|
OOFEM 3.0
|
#include <quad1mindlinshell3d.h>
Public Member Functions | |
| Quad1MindlinShell3D (int n, Domain *d) | |
| FEInterpolation * | giveInterpolation () const override |
| FEInterpolation * | giveInterpolation (DofIDItem id) const override |
| MaterialMode | giveMaterialMode () override |
| int | testElementExtension (ElementExtension ext) override |
| void | computeStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) override |
| void | giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override |
| const char * | giveInputRecordName () const override |
| const char * | giveClassName () const override |
| void | initializeFrom (InputRecord &ir, int priority) override |
| int | computeNumberOfDofs () override |
| int | computeNumberOfGlobalDofs () override |
| void | giveDofManDofIDMask (int inode, IntArray &) const override |
| Element_Geometry_Type | giveGeometryType () const override |
| void | computeMidPlaneNormal (FloatArray &answer, const GaussPoint *gp) override |
| double | giveCharacteristicLength (const FloatArray &normalToCrackPlane) override |
| double | computeVolumeAround (GaussPoint *gp) override |
| void | computeLumpedMassMatrix (FloatMatrix &answer, TimeStep *tStep) override |
| void | computeMassMatrix (FloatMatrix &answer, TimeStep *tStep) override |
| int | giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override |
| virtual void | computeLCS () |
| bool | computeGtoLRotationMatrix (FloatMatrix &answer) override |
| Interface * | giveInterface (InterfaceType it) override |
| void | SPRNodalRecoveryMI_giveSPRAssemblyPoints (IntArray &pap) override |
| void | SPRNodalRecoveryMI_giveDofMansDeterminedByPatch (IntArray &answer, int pap) override |
| int | SPRNodalRecoveryMI_giveNumberOfIP () override |
| SPRPatchType | SPRNodalRecoveryMI_givePatchType () override |
| Public Member Functions inherited from oofem::StructuralElement | |
| StructuralElement (int n, Domain *d) | |
| virtual | ~StructuralElement () |
| Destructor. | |
| void | giveCharacteristicMatrix (FloatMatrix &answer, CharType, TimeStep *tStep) override |
| void | giveCharacteristicVector (FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep) override |
| virtual void | computeConsistentMassMatrix (FloatMatrix &answer, TimeStep *tStep, double &mass, const double *ipDensity=NULL) |
| virtual void | giveMassMtrxIntegrationgMask (IntArray &answer) |
| void | computeStiffnessMatrix_withIRulesAsSubcells (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) |
| virtual void | computeInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep) |
| virtual void | computeLumpedInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep) |
| void | computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer) override |
| virtual void | giveInternalForcesVector_withIRulesAsSubcells (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) |
| int | giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override |
| virtual void | computeResultingIPTemperatureAt (FloatArray &answer, TimeStep *tStep, GaussPoint *gp, ValueModeType mode) |
| virtual void | computeResultingIPEigenstrainAt (FloatArray &answer, TimeStep *tStep, GaussPoint *gp, ValueModeType mode) |
| virtual void | computeNmatrixAt (const FloatArray &iLocCoord, FloatMatrix &answer) |
| void | updateBeforeNonlocalAverage (TimeStep *tStep) override |
| virtual void | giveNonlocalLocationArray (IntArray &locationArray, const UnknownNumberingScheme &us) |
| virtual void | addNonlocalStiffnessContributions (SparseMtrx &dest, const UnknownNumberingScheme &s, TimeStep *tStep) |
| int | adaptiveUpdate (TimeStep *tStep) override |
| void | updateInternalState (TimeStep *tStep) override |
| void | updateYourself (TimeStep *tStep) override |
| int | checkConsistency () override |
| void | giveInputRecord (DynamicInputRecord &input) override |
| const char * | giveClassName () const override |
| int | giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *tStep) override |
| void | showSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep) override |
| Shows sparse structure. | |
| void | showExtendedSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep) override |
| Shows extended sparse structure (for example, due to nonlocal interactions for tangent stiffness). | |
| void | computeLoadVector (FloatArray &answer, BodyLoad *load, CharType type, ValueModeType mode, TimeStep *tStep) override |
| void | computeBoundarySurfaceLoadVector (FloatArray &answer, BoundaryLoad *load, int boundary, CharType type, ValueModeType mode, TimeStep *tStep, bool global=true) override |
| void | computeBoundaryEdgeLoadVector (FloatArray &answer, BoundaryLoad *load, int boundary, CharType type, ValueModeType mode, TimeStep *tStep, bool global=true) override |
| virtual void | computeEdgeNMatrix (FloatMatrix &answer, int boundaryID, const FloatArray &lcoords) |
| computes edge interpolation matrix | |
| virtual void | computeSurfaceNMatrix (FloatMatrix &answer, int boundaryID, const FloatArray &lcoords) |
| StructuralCrossSection * | giveStructuralCrossSection () |
| Helper function which returns the structural cross-section for the element. | |
| virtual void | createMaterialStatus () |
| Public Member Functions inherited from oofem::Element | |
| Element (int n, Domain *aDomain) | |
| Element (const Element &src)=delete | |
| Element & | operator= (const Element &src)=delete |
| virtual | ~Element () |
| Virtual destructor. | |
| void | giveLocationArray (IntArray &locationArray, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) const |
| void | giveLocationArray (IntArray &locationArray, const IntArray &dofIDMask, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) const |
| virtual void | giveBoundaryLocationArray (IntArray &locationArray, const IntArray &bNodes, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) |
| virtual void | giveBoundaryLocationArray (IntArray &locationArray, const IntArray &bNodes, const IntArray &dofIDMask, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) |
| virtual int | giveNumberOfDofs () |
| virtual int | giveNumberOfInternalDofManagers () const |
| virtual DofManager * | giveInternalDofManager (int i) const |
| virtual void | setInternalDofManager (int num, std::unique_ptr< DofManager > dm) |
| virtual double | giveCharacteristicValue (CharType type, TimeStep *tStep) |
| virtual void | computeTangentFromSurfaceLoad (FloatMatrix &answer, BoundaryLoad *load, int boundary, MatResponseMode rmode, TimeStep *tStep) |
| virtual void | computeTangentFromEdgeLoad (FloatMatrix &answer, BoundaryLoad *load, int boundary, MatResponseMode rmode, TimeStep *tStep) |
| const IntArray & | giveBodyLoadList () const |
| const IntArray & | giveBoundaryLoadList () const |
| void | computeVectorOf (ValueModeType u, TimeStep *tStep, FloatArray &answer) |
| void | computeVectorOf (const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false) |
| void | computeBoundaryVectorOf (const IntArray &bNodes, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false) |
| void | computeVectorOf (PrimaryField &field, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false) |
| void | computeVectorOfPrescribed (ValueModeType u, TimeStep *tStep, FloatArray &answer) |
| void | computeVectorOfPrescribed (const IntArray &dofIDMask, ValueModeType type, TimeStep *tStep, FloatArray &answer) |
| int | computeNumberOfPrimaryMasterDofs () |
| virtual bool | giveRotationMatrix (FloatMatrix &answer) |
| virtual bool | computeDofTransformationMatrix (FloatMatrix &answer, const IntArray &nodes, bool includeInternal) |
| virtual void | giveInternalDofManDofIDMask (int inode, IntArray &answer) const |
| virtual void | giveElementDofIDMask (IntArray &answer) const |
| virtual double | computeVolumeAreaOrLength () |
| Computes the volume, area or length of the element depending on its spatial dimension. | |
| double | computeMeanSize () |
| virtual double | computeVolume () |
| virtual double | computeArea () |
| virtual double | computeLength () |
| virtual IntArray | giveBoundaryEdgeNodes (int boundary, bool includeHierarchical=false) const |
| virtual IntArray | giveBoundarySurfaceNodes (int boundary, bool includeHierarchical=false) const |
| virtual IntArray | giveBoundaryNodes (int boundary) const |
| virtual std::unique_ptr< IntegrationRule > | giveBoundaryEdgeIntegrationRule (int order, int boundary) |
| virtual std::unique_ptr< IntegrationRule > | giveBoundarySurfaceIntegrationRule (int order, int boundary) |
| int | giveDofManagerNumber (int i) const |
| const IntArray & | giveDofManArray () const |
| void | addDofManager (DofManager *dMan) |
| DofManager * | giveDofManager (int i) const |
| Node * | giveNode (int i) const |
| virtual ElementSide * | giveSide (int i) const |
| virtual const FEInterpolation * | getGeometryInterpolation () const |
| virtual Material * | giveMaterial () |
| int | giveMaterialNumber () const |
| CrossSection * | giveCrossSection () |
| int | getActivityTimeFunctionNumber () |
| void | setActivityTimeFunctionNumber (int funcIndx) |
| void | setMaterial (int matIndx) |
| virtual void | setCrossSection (int csIndx) |
| virtual int | giveNumberOfDofManagers () const |
| void | setNumberOfDofManagers (int i) |
| Sets number of element dof managers. | |
| virtual int | giveNumberOfNodes () const |
| void | setDofManagers (const IntArray &dmans) |
| void | setDofManager (int id, int dm) |
| void | setBodyLoads (const IntArray &bodyLoads) |
| void | setIntegrationRules (std ::vector< std ::unique_ptr< IntegrationRule > > irlist) |
| virtual integrationDomain | giveIntegrationDomain () const |
| virtual int | giveIntegrationRuleLocalCodeNumbers (IntArray &answer, IntegrationRule &ie) |
| int | giveRegionNumber () |
| virtual void | initializeYourself (TimeStep *timeStepWhenICApply) |
| virtual bool | isActivated (TimeStep *tStep) |
| virtual bool | isCast (TimeStep *tStep) |
| virtual void | initForNewStep () |
| virtual Element_Geometry_Type | giveEdgeGeometryType (int id) const |
| Returns the receiver edge geometry type. | |
| virtual Element_Geometry_Type | giveSurfaceGeometryType (int id) const |
| Returns the receiver surface geometry type. | |
| virtual int | giveSpatialDimension () |
| virtual int | giveNumberOfBoundarySides () |
| Returns number of boundaries (entities of element_dimension-1: points, edges, surfaces). | |
| virtual int | giveNumberOfEdges () const |
| virtual int | giveNumberOfSurfaces () const |
| virtual int | giveDefaultIntegrationRule () const |
| virtual IntegrationRule * | giveDefaultIntegrationRulePtr () |
| int | giveNumberOfIntegrationRules () |
| virtual IntegrationRule * | giveIntegrationRule (int i) |
| std::vector< std ::unique_ptr< IntegrationRule > > & | giveIntegrationRulesArray () |
| int | giveGlobalIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) |
| virtual double | giveLengthInDir (const FloatArray &normalToCrackPlane) |
| double | giveCharacteristicLengthForPlaneElements (const FloatArray &normalToCrackPlane) |
| double | giveCharacteristicLengthForAxisymmElements (const FloatArray &normalToCrackPlane) |
| virtual double | giveCharacteristicSize (GaussPoint *gp, FloatArray &normalToCrackPlane, ElementCharSizeMethod method) |
| virtual double | giveParentElSize () const |
| virtual int | computeGlobalCoordinates (FloatArray &answer, const FloatArray &lcoords) |
| virtual bool | computeLocalCoordinates (FloatArray &answer, const FloatArray &gcoords) |
| virtual int | giveLocalCoordinateSystem (FloatMatrix &answer) |
| virtual void | giveLocalCoordinateSystemVector (InternalStateType isttype, FloatArray &answer) |
| virtual int | adaptiveMap (Domain *oldd, TimeStep *tStep) |
| virtual int | mapStateVariables (Domain &iOldDom, const TimeStep &iTStep) |
| virtual int | adaptiveFinish (TimeStep *tStep) |
| void | updateLocalNumbering (EntityRenumberingFunctor &f) override |
| template<class T> | |
| void | ipEvaluator (T *src, void(T ::*f)(GaussPoint *gp)) |
| Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f. | |
| template<class T, class S> | |
| void | ipEvaluator (T *src, void(T ::*f)(GaussPoint *, S &), S &_val) |
| Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f as well as additional array. | |
| virtual void | drawYourself (oofegGraphicContext &gc, TimeStep *tStep) |
| virtual void | drawAnnotation (oofegGraphicContext &gc, TimeStep *tStep) |
| virtual void | drawRawGeometry (oofegGraphicContext &gc, TimeStep *tStep) |
| virtual void | drawDeformedGeometry (oofegGraphicContext &gc, TimeStep *tStep, UnknownType) |
| virtual void | drawScalar (oofegGraphicContext &gc, TimeStep *tStep) |
| virtual void | drawSpecial (oofegGraphicContext &gc, TimeStep *tStep) |
| virtual void | giveLocalIntVarMaxMin (oofegGraphicContext &gc, TimeStep *tStep, double &emin, double &emax) |
| virtual int | giveInternalStateAtSide (FloatArray &answer, InternalStateType type, InternalStateMode mode, int side, TimeStep *tStep) |
| int | giveLabel () const |
| int | giveGlobalNumber () const |
| void | setGlobalNumber (int num) |
| elementParallelMode | giveParallelMode () const |
| void | setParallelMode (elementParallelMode _mode) |
| Sets parallel mode of element. | |
| virtual elementParallelMode | giveKnotSpanParallelMode (int) const |
| int | packUnknowns (DataStream &buff, TimeStep *tStep) |
| int | unpackAndUpdateUnknowns (DataStream &buff, TimeStep *tStep) |
| int | estimatePackSize (DataStream &buff) |
| const IntArray * | givePartitionList () const |
| void | setPartitionList (IntArray &pl) |
| virtual double | predictRelativeComputationalCost () |
| virtual double | giveRelativeSelfComputationalCost () |
| virtual double | predictRelativeRedistributionCost () |
| IntArray * | giveBodyLoadArray () |
| Returns array containing load numbers of loads acting on element. | |
| IntArray * | giveBoundaryLoadArray () |
| Returns array containing load numbers of boundary loads acting on element. | |
| void | initializeFinish () override |
| void | postInitialize () override |
| Performs post initialization steps. | |
| void | saveContext (DataStream &stream, ContextMode mode) override |
| void | restoreContext (DataStream &stream, ContextMode mode) override |
| void | printOutputAt (FILE *file, TimeStep *tStep) override |
| virtual const IntArray | giveLocation () |
| virtual void | recalculateCoordinates (int nodeNumber, FloatArray &coords) |
| void | setSharedEdgeID (int iedge, int globalID) |
| void | setSharedSurfaceID (int isurf, int globalID) |
| const IntArray * | giveSharedEdgeIDs () const |
| const IntArray * | giveSharedSurfaceIDs () const |
| Public Member Functions inherited from oofem::FEMComponent | |
| FEMComponent (int n, Domain *d) | |
| virtual | ~FEMComponent ()=default |
| Virtual destructor. | |
| Domain * | giveDomain () const |
| virtual void | setDomain (Domain *d) |
| int | giveNumber () const |
| void | setNumber (int num) |
| virtual void | initializeFrom (InputRecord &ir) |
| virtual void | printYourself () |
| Prints receiver state on stdout. Useful for debugging. | |
| std::string | errorInfo (const char *func) const |
| Returns string for prepending output (used by error reporting macros). | |
| Public Member Functions inherited from oofem::ZZNodalRecoveryModelInterface | |
| ZZNodalRecoveryModelInterface (Element *element) | |
| Constructor. | |
| virtual bool | ZZNodalRecoveryMI_computeNValProduct (FloatMatrix &answer, InternalStateType type, TimeStep *tStep) |
| virtual void | ZZNodalRecoveryMI_computeNNMatrix (FloatArray &answer, InternalStateType type) |
| Public Member Functions inherited from oofem::Interface | |
| Interface () | |
| Constructor. | |
| virtual | ~Interface () |
| Public Member Functions inherited from oofem::SPRNodalRecoveryModelInterface | |
| SPRNodalRecoveryModelInterface () | |
| Constructor. | |
Protected Member Functions | |
| void | computeGaussPoints () override |
| void | computeBodyLoadVectorAt (FloatArray &answer, Load *load, TimeStep *tStep, ValueModeType mode) override |
| void | computeBmatrixAt (GaussPoint *gp, FloatMatrix &answer, int=1, int=ALL_STRAINS) override |
| void | computeStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) override |
| void | computeStressVector (FloatArray &answer, const FloatArray &strain, GaussPoint *gp, TimeStep *tStep) override |
| void | computeConstitutiveMatrixAt (FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep) override |
| void | giveEdgeDofMapping (IntArray &answer, int iEdge) const override |
| double | computeEdgeVolumeAround (GaussPoint *gp, int iEdge) override |
| int | computeLoadLEToLRotationMatrix (FloatMatrix &answer, int iEdge, GaussPoint *gp) override |
| void | computeVectorOfUnknowns (ValueModeType mode, TimeStep *tStep, FloatArray &shellUnknowns, FloatArray &drillUnknowns) |
| virtual void | computePointLoadVectorAt (FloatArray &answer, Load *load, TimeStep *tStep, ValueModeType mode, bool global=true) |
| virtual void | giveSurfaceDofMapping (IntArray &answer, int iSurf) const |
| virtual double | computeSurfaceVolumeAround (GaussPoint *gp, int iSurf) |
| virtual int | computeLoadGToLRotationMtrx (FloatMatrix &answer) |
| virtual int | computeLoadLSToLRotationMatrix (FloatMatrix &answer, int iSurf, GaussPoint *gp) |
| virtual int | giveNumberOfIPForMassMtrxIntegration () |
| void | condense (FloatMatrix *stiff, FloatMatrix *mass, FloatArray *load, IntArray *what) |
| virtual void | setupIRForMassMtrxIntegration (IntegrationRule &iRule) |
Protected Attributes | |
| std::vector< FloatArray > | lnodes |
| Cached nodal coordinates in local c.s.,. | |
| FloatMatrix | lcsMatrix |
| Cached coordinates in local c.s.,. | |
| bool | reducedIntegrationFlag |
| Flag controlling reduced (one - point) integration for shear. | |
| Protected Attributes inherited from oofem::StructuralElement | |
| std::unique_ptr< FloatArray > | initialDisplacements |
| Initial displacement vector, describes the initial nodal displacements when element has been casted. | |
| Protected Attributes inherited from oofem::Element | |
| int | numberOfDofMans |
| Number of dofmanagers. | |
| IntArray | dofManArray |
| Array containing dofmanager numbers. | |
| int | material |
| Number of associated material. | |
| int | crossSection |
| Number of associated cross section. | |
| IntArray | bodyLoadArray |
| IntArray | boundaryLoadArray |
| std::vector< std ::unique_ptr< IntegrationRule > > | integrationRulesArray |
| FloatMatrix | elemLocalCS |
| Transformation material matrix, used in orthotropic and anisotropic materials, global->local transformation. | |
| int | activityTimeFunction |
| Element activity time function. If defined, nonzero value indicates active receiver, zero value inactive element. | |
| int | globalNumber |
| int | numberOfGaussPoints |
| elementParallelMode | parallel_mode |
| Determines the parallel mode of the element. | |
| IntArray | partitions |
| IntArray | globalEdgeIDs |
| IntArray | globalSurfaceIDs |
| Protected Attributes inherited from oofem::FEMComponent | |
| int | number |
| Component number. | |
| Domain * | domain |
| Link to domain object, useful for communicating with other FEM components. | |
Static Protected Attributes | |
| static FEI2dQuadLin | interp |
| static IntArray | shellOrdering = { 1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, 20, 21, 22, 23 } |
| Ordering for the normal shell stiffness (everything but the out-of-plane rotations). | |
| static IntArray | drillOrdering = { 6, 12, 18, 24 } |
| Ordering for the drilling dofs (the out-of-plane rotations). | |
| static ParamKey | IPK_Quad1MindlinShell3D_reducedIntegration |
Additional Inherited Members | |
| Static Public Attributes inherited from oofem::Element | |
| static ParamKey | IPK_Element_mat |
| static ParamKey | IPK_Element_crosssect |
| static ParamKey | IPK_Element_nodes |
| static ParamKey | IPK_Element_bodyload |
| static ParamKey | IPK_Element_boundaryload |
| static ParamKey | IPK_Element_lcs |
| static ParamKey | IPK_Element_partitions |
| static ParamKey | IPK_Element_remote |
| static ParamKey | IPK_Element_activityTimeFunction |
| static ParamKey | IPK_Element_nip |
This class implements an quadrilateral four-node shell element, using Mindlin plate theory. Each node has 3 degrees of freedom (out-of-plane displacement, in-plane rotations). This type of element exhibit strong shear locking (thin plates exhibit almost no bending). No reduced integration is used, as it causes numerical problems.
Loading types supported;
Reference: Robert Cook, David Malkus, Michael Plesha Concepts and Applications of Finite Element Analysis - Third edition ISBN: 0-471-84788-7
Definition at line 69 of file quad1mindlinshell3d.h.
| oofem::Quad1MindlinShell3D::Quad1MindlinShell3D | ( | int | n, |
| Domain * | d ) |
Definition at line 65 of file quad1mindlinshell3d.C.
References lnodes, oofem::Element::numberOfDofMans, oofem::Element::numberOfGaussPoints, reducedIntegrationFlag, oofem::SPRNodalRecoveryModelInterface::SPRNodalRecoveryModelInterface(), oofem::StructuralElement::StructuralElement(), and oofem::ZZNodalRecoveryModelInterface::ZZNodalRecoveryModelInterface().
|
overrideprotectedvirtual |
Computes the geometrical matrix of receiver in given integration point. The product of this matrix (assembled at given integration point) and element displacement vector is element strain vector. If lowerIndx and upperIndx parameters are specified, answer is formed only for strains within this interval. This will affects the size of answer.
| gp | Integration point for which answer is computed. |
| answer | Geometric matrix of receiver. |
| lowerIndx | If specified, answer is formed only for strain with index equal and greater than lowerIndx. This parameter has default value 1 (answer is formed from first strain). |
| upperIndx | If specified, answer is formed only for strain with index less and equal than upperIndx. This parameter has default value ALL_STRAINS (answer is formed for all strains). |
Implements oofem::StructuralElement.
Definition at line 188 of file quad1mindlinshell3d.C.
References oofem::FloatMatrix::giveDeterminant(), giveGeometryType(), oofem::GaussPoint::giveNaturalCoordinates(), interp, lnodes, oofem::FloatArray::normalize(), OOFEM_WARNING, reducedIntegrationFlag, oofem::FloatMatrix::resize(), oofem::sqr(), and oofem::FloatMatrix::zero().
Referenced by computeStiffnessMatrix(), computeStrainVector(), and giveInternalForcesVector().
|
overrideprotectedvirtual |
Computes the load vector due to body load acting on receiver, at given time step. Default implementation computes body load vector numerically as \( l=\int_V N^{\mathrm{T}} f \rho\;\mathrm{d}V \) using default integration rule. Result is transformed to global c.s.
| answer | Computed load vector due to body load |
| load | Body load which contribution is computed. |
| tStep | Time step. |
| mode | determines the response mode |
Reimplemented from oofem::StructuralElement.
Definition at line 105 of file quad1mindlinshell3d.C.
References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::FloatArray::beProductOf(), oofem::BodyLoadBGT, oofem::FloatArray::clear(), oofem::Load::computeComponentArrayAt(), computeVolumeAround(), oofem::CS_Thickness, oofem::ForceLoadBVT, oofem::CrossSection::give(), oofem::GeneralBoundaryCondition::giveBCGeoType(), oofem::GeneralBoundaryCondition::giveBCValType(), oofem::Element::giveCrossSection(), oofem::FloatArray::giveSize(), oofem::StructuralElement::giveStructuralCrossSection(), oofem::Element::integrationRulesArray, interp, lcsMatrix, OOFEM_ERROR, oofem::FloatArray::resize(), and oofem::FloatArray::zero().
|
overrideprotectedvirtual |
Computes constitutive matrix of receiver. Default implementation uses element cross section giveCharMaterialStiffnessMatrix service.
| answer | Constitutive matrix. |
| rMode | Material response mode of answer. |
| gp | Integration point for which constitutive matrix is computed. |
| tStep | Time step. |
Implements oofem::StructuralElement.
Definition at line 333 of file quad1mindlinshell3d.C.
References oofem::StructuralCrossSection::give3dShellStiffMtrx(), and oofem::StructuralElement::giveStructuralCrossSection().
Referenced by computeStiffnessMatrix().
|
overrideprotectedvirtual |
Computes volume related to integration point on local edge.
| gp | edge integration point |
| iEdge | edge number |
Reimplemented from oofem::StructuralElement.
Definition at line 575 of file quad1mindlinshell3d.C.
References giveGeometryType(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::GaussPoint::giveWeight(), interp, and lnodes.
|
overrideprotectedvirtual |
Initializes the array of integration rules member variable. Element can have multiple integration rules for different tasks. For example structural element family class uses this feature to implement transparent support for reduced and selective integration of some strain components. Must be defined by terminator classes.
Reimplemented from oofem::Element.
Definition at line 91 of file quad1mindlinshell3d.C.
References computeLCS(), oofem::Element::giveCrossSection(), oofem::Element::integrationRulesArray, oofem::Element::numberOfGaussPoints, and oofem::CrossSection::setupIntegrationPoints().
|
overridevirtual |
Returns transformation matrix from global c.s. to local element c.s., i.e. \( r_l =T r_g \). If no transformation is necessary then answer is empty matrix and zero value is returned.
| answer | Computed rotation matrix. |
Reimplemented from oofem::Element.
Definition at line 643 of file quad1mindlinshell3d.C.
References lcsMatrix, oofem::FloatMatrix::resize(), oofem::FloatMatrix::setSubMatrix(), and oofem::FloatMatrix::zero().
|
virtual |
Definition at line 619 of file quad1mindlinshell3d.C.
References oofem::cross(), oofem::Element::giveNode(), lcsMatrix, lnodes, and oofem::normalize().
Referenced by computeGaussPoints().
|
overrideprotectedvirtual |
Returns transformation matrix from local edge c.s to element local coordinate system of load vector components. Necessary, because integration must be done in local coordinate system of entity (edge or surface). If no transformation is necessary, answer is empty matrix (default);
| answer | Computed rotation matrix. |
| iEdge | Edge number. |
| gp | Integration point (point, where transformation is computed, useful for curved edges). |
Reimplemented from oofem::StructuralElement.
Definition at line 594 of file quad1mindlinshell3d.C.
References oofem::FloatMatrix::at(), oofem::Element::giveNode(), interp, length(), oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().
|
overridevirtual |
Computes lumped mass matrix of receiver. Default implementation returns lumped consistent mass matrix. Then returns lumped mass transformed into nodal coordinate system. The lumping procedure zeroes all off-diagonal members and zeroes also all diagonal members corresponding to non-displacement DOFs. Such diagonal matrix is then rescaled, to preserve the element mass. Requires the computeNmatrixAt and giveMassMtrxIntegrationgMask services to be implemented.
| answer | Lumped mass matrix. |
| tStep | Time step. |
Reimplemented from oofem::StructuralElement.
Definition at line 501 of file quad1mindlinshell3d.C.
References computeVolumeAround(), oofem::CrossSection::give(), oofem::StructuralElement::giveStructuralCrossSection(), oofem::Element::integrationRulesArray, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().
Referenced by computeMassMatrix().
|
inlineoverridevirtual |
Computes mass matrix of receiver. Default implementation returns consistent mass matrix and uses numerical integration. Returns result of this->computeConsistentMassMatrix service, transformed into nodal coordinate system. Requires the computeNmatrixAt and giveMassMtrxIntegrationgMask services to be implemented.
| answer | Mass matrix. |
| tStep | Time step. |
Reimplemented from oofem::StructuralElement.
Definition at line 119 of file quad1mindlinshell3d.h.
References computeLumpedMassMatrix().
|
overridevirtual |
Computes mid-plane normal of receiver at integration point. Only for plane elements in space (3d) (shells, plates, ....).
| answer | The mid plane normal. |
| gp | The integration point at which to calculate the normal. |
Reimplemented from oofem::Element.
Definition at line 474 of file quad1mindlinshell3d.C.
References oofem::cross(), oofem::DofManager::giveCoordinates(), oofem::Element::giveNode(), and oofem::norm().
|
inlineoverridevirtual |
Computes or simply returns total number of element's local DOFs. Must be defined by particular element.
Reimplemented from oofem::Element.
Definition at line 107 of file quad1mindlinshell3d.h.
|
inlineoverridevirtual |
Computes the total number of element's global dofs. The transitions from global c.s. to nodal c.s. should NOT be included.
Reimplemented from oofem::Element.
Definition at line 108 of file quad1mindlinshell3d.h.
|
overridevirtual |
Computes numerically stiffness matrix of receiver. Default implementation computes element stiffness using \( K=\int_v B^{\mathrm{T}} D B \mathrm{d}V \) formulae, where \( B \) is element geometric matrix and \( D \) is material stiffness matrix. No geometrical nonlinearity is taken into account. NUmerical integration procedure uses integrationRulesArray for numerical integration. Support for reduced or selected integration is implemented. The individual integration rules are assumed to correspond to different terms from which the overall matrix is assembled.
If there is one integration rule, the full integration of all coefficients is performed. Otherwise, integration is performed using following rules. Each integration rule can specify start and end strain index of strain vector components for which is valid. It is necessary to ensure that these start and end indexes, dividing geometrical matrix into blocks, are not overlapping and that each strain component is included.
Then stiffness matrix is obtained as summation of integrals \( I_{ij}=\int_v B^{\mathrm{T}}_i D_{ij} B_j \mathrm{d}V \) where \( B_i \) is i-th block of geometrical matrix and \( D_{ij} \) is corresponding constitutive sub-matrix. The geometrical matrix is obtained using computeBmatrixAt service and the constitutive matrix is obtained using computeConstitutiveMatrixAt service. The \( I_{ij} \) integral is evaluated using such integration rule, which is valid for i-th or j-th block and has smaller number of integration points.
For higher numerical performance, only one half of stiffness matrix is computed and answer is then symmetrized. Therefore, if element matrix will be generally nonsymmetric, one must specialize this method. Finally, the result is transformed into global coordinate system (or nodal coordinate system, if it is defined).
| answer | Computed stiffness matrix (symmetric). |
| rMode | Response mode. |
| tStep | Time step. |
Reimplemented from oofem::StructuralElement.
Definition at line 414 of file quad1mindlinshell3d.C.
References oofem::FloatMatrix::assemble(), oofem::FloatMatrix::beProductOf(), computeBmatrixAt(), computeConstitutiveMatrixAt(), computeVolumeAround(), oofem::CS_DrillingStiffness, drillOrdering, oofem::CrossSection::give(), oofem::StructuralElement::giveStructuralCrossSection(), oofem::Element::integrationRulesArray, interp, oofem::FloatMatrix::plusDyadSymmUpper(), oofem::FloatMatrix::plusProductSymmUpper(), oofem::FloatMatrix::resize(), shellOrdering, oofem::FloatMatrix::symmetrized(), and oofem::FloatMatrix::zero().
|
overrideprotectedvirtual |
Compute strain vector of receiver evaluated at given integration point at given time step from element displacement vector. The nature of strain vector depends on the element type.
| answer | Requested strain vector. |
| gp | Integration point where to calculate the strain. |
| tStep | Time step. |
Reimplemented from oofem::StructuralElement.
Definition at line 350 of file quad1mindlinshell3d.C.
References oofem::FloatArray::beProductOf(), computeBmatrixAt(), and computeVectorOfUnknowns().
|
overrideprotectedvirtual |
Computes the stress vector of receiver at given integration point, at time step tStep. The nature of these stresses depends on the element's type.
| answer | Stress vector. |
| strain | Strain vector. |
| gp | Integration point. |
| tStep | Time step. |
Implements oofem::StructuralElement.
Definition at line 326 of file quad1mindlinshell3d.C.
References oofem::StructuralCrossSection::giveGeneralizedStress_Shell(), and oofem::StructuralElement::giveStructuralCrossSection().
|
protected |
Definition at line 340 of file quad1mindlinshell3d.C.
References oofem::FloatArray::beSubArrayOf(), oofem::Element::computeVectorOf(), drillOrdering, and shellOrdering.
Referenced by computeStrainVector(), and giveInternalForcesVector().
|
overridevirtual |
Returns volume related to given integration point. Used typically in subroutines, that perform integration over element volume. Should be implemented by particular elements.
| gp | Integration point for which volume is computed. |
Reimplemented from oofem::Element.
Definition at line 492 of file quad1mindlinshell3d.C.
References giveGeometryType(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::GaussPoint::giveWeight(), interp, and lnodes.
Referenced by computeBodyLoadVectorAt(), computeLumpedMassMatrix(), computeStiffnessMatrix(), and giveInternalForcesVector().
|
overridevirtual |
Returns the size of element in the given direction, in some cases adjusted (e.g. if the direction is perpendicular to a planar element). Required by material models relying on the crack-band approach to achieve objectivity with respect to the mesh size.
| normalToCrackPlane | Normal to the expected crack band. |
Reimplemented from oofem::Element.
Definition at line 485 of file quad1mindlinshell3d.C.
References oofem::Element::giveLengthInDir().
|
inlineoverridevirtual |
Reimplemented from oofem::Element.
Definition at line 104 of file quad1mindlinshell3d.h.
|
overridevirtual |
Returns dofmanager dof mask for node. This mask defines the dofs which are used by element in node. Mask influences the code number ordering for particular node. Code numbers are ordered according to node order and dofs belonging to particular node are ordered according to this mask. If element requests dofs using node mask which are not in node then error is generated. This masking allows node to be shared by different elements with different dofs in same node. Elements local code numbers are extracted from node using this mask. Must be defined by particular element.
| inode | Mask is computed for local dofmanager with inode number. |
| answer | Mask for node. |
Reimplemented from oofem::Element.
Definition at line 467 of file quad1mindlinshell3d.C.
|
overrideprotectedvirtual |
Assembles edge dof mapping mask, which provides mapping between edge local DOFs and "global" element DOFs. Mask can be imagined as local edge code numbers used to localize local edge DOFs to element DOFs.
| answer | Edge DOF mapping. |
| iEdge | Edge number. |
Reimplemented from oofem::StructuralElement.
Definition at line 558 of file quad1mindlinshell3d.C.
References OOFEM_ERROR.
|
inlineoverridevirtual |
Returns the element geometry type. This information is assumed to be of general interest, but it is required only for some specialized tasks.
Implements oofem::Element.
Definition at line 110 of file quad1mindlinshell3d.h.
Referenced by computeBmatrixAt(), computeEdgeVolumeAround(), and computeVolumeAround().
|
inlineoverridevirtual |
Implements oofem::FEMComponent.
Definition at line 103 of file quad1mindlinshell3d.h.
References _IFT_Quad1MindlinShell3D_Name.
|
overridevirtual |
Interface requesting service.
Reimplemented from oofem::FEMComponent.
Definition at line 658 of file quad1mindlinshell3d.C.
References oofem::SPRNodalRecoveryModelInterface::SPRNodalRecoveryModelInterface(), oofem::SPRNodalRecoveryModelInterfaceType, oofem::ZZNodalRecoveryModelInterface::ZZNodalRecoveryModelInterface(), and oofem::ZZNodalRecoveryModelInterfaceType.
|
overridevirtual |
Returns equivalent nodal forces vectors. Useful for nonlinear analysis. Default implementation computes result as \( F=\int_v B^{\mathrm{T}} \sigma \mathrm{d}V \), where \( \sigma \) is the real element stress vector obtained using computeStressVector service (if useUpdatedGpRecord=0) or (if useUpdatedGpRecord=1) from integration point status. The geometric matrix is obtained using computeBmatrixAt service. Integration is performed using default integration rule, which should produce always valid results, assuming that strains used for computation of stresses are valid.
| answer | Internal nodal forces vector. |
| tStep | Time step. |
| useUpdatedGpRecord | If equal to zero, the stresses in integration points are computed (slow but safe), else if nonzero the stresses are taken directly from integration point status (should be derived from StructuralMaterialStatus) (fast, but engineering model must ensure valid status data in each integration point). |
Reimplemented from oofem::StructuralElement.
Definition at line 363 of file quad1mindlinshell3d.C.
References oofem::FloatArray::add(), oofem::FloatArray::assemble(), oofem::FloatArray::beProductOf(), computeBmatrixAt(), computeVectorOfUnknowns(), computeVolumeAround(), oofem::CS_DrillingStiffness, oofem::FloatArray::dotProduct(), drillOrdering, oofem::CrossSection::give(), oofem::StructuralCrossSection::giveGeneralizedStress_Shell(), oofem::StructuralMaterialStatus::giveStressVector(), oofem::StructuralElement::giveStructuralCrossSection(), oofem::Element::integrationRulesArray, interp, oofem::FloatArray::plusProduct(), oofem::FloatArray::resize(), shellOrdering, and oofem::FloatArray::zero().
|
overridevirtual |
Reimplemented from oofem::Element.
Definition at line 77 of file quad1mindlinshell3d.C.
References interp.
|
overridevirtual |
Returns the interpolation for the specific dof id. Special elements which uses a mixed interpolation should reimplement this method.
| id | ID of the dof for the for the requested interpolation. |
Reimplemented from oofem::Element.
Definition at line 84 of file quad1mindlinshell3d.C.
References interp.
|
overridevirtual |
Returns the integration point corresponding value in full form.
| answer | Contain corresponding integration point value, zero sized if not available. |
| gp | Integration point to check. |
| type | Determines the type of internal variable. |
| tStep | Time step. |
Reimplemented from oofem::Element.
Definition at line 521 of file quad1mindlinshell3d.C.
References oofem::FloatArray::at(), oofem::StructuralElement::giveIPValue(), oofem::GaussPoint::giveMaterialStatus(), oofem::StructuralMaterialStatus::giveStrainVector(), oofem::StructuralMaterialStatus::giveStressVector(), and oofem::FloatArray::resize().
|
inlineoverridevirtual |
Returns material mode for receiver integration points. Should be specialized.
Reimplemented from oofem::Element.
Definition at line 96 of file quad1mindlinshell3d.h.
|
overridevirtual |
Reimplemented from oofem::Element.
Definition at line 458 of file quad1mindlinshell3d.C.
References oofem::Domain::elementPPM, oofem::FEMComponent::giveDomain(), oofem::Element::initializeFrom(), oofem::Element::IPK_Element_activityTimeFunction, oofem::FEMComponent::number, PM_UPDATE_PARAMETER, and reducedIntegrationFlag.
|
overridevirtual |
Implements oofem::SPRNodalRecoveryModelInterface.
Definition at line 681 of file quad1mindlinshell3d.C.
References oofem::IntArray::at(), oofem::Element::giveNode(), oofem::FEMComponent::giveNumber(), OOFEM_ERROR, and oofem::IntArray::resize().
|
inlineoverridevirtual |
Implements oofem::SPRNodalRecoveryModelInterface.
Definition at line 131 of file quad1mindlinshell3d.h.
References oofem::Element::numberOfGaussPoints.
|
inlineoverridevirtual |
Implements oofem::SPRNodalRecoveryModelInterface.
Definition at line 132 of file quad1mindlinshell3d.h.
References oofem::SPRPatchType_2dxy.
|
overridevirtual |
Implements oofem::SPRNodalRecoveryModelInterface.
Definition at line 672 of file quad1mindlinshell3d.C.
References oofem::IntArray::at(), oofem::Element::giveNode(), oofem::FEMComponent::giveNumber(), and oofem::IntArray::resize().
|
inlineoverridevirtual |
Tests if the element implements required extension. ElementExtension type defines the list of all available element extensions.
| ext | Extension to be tested. |
Reimplemented from oofem::Element.
Definition at line 97 of file quad1mindlinshell3d.h.
References oofem::Element_EdgeLoadSupport.
|
staticprotected |
Ordering for the drilling dofs (the out-of-plane rotations).
Definition at line 86 of file quad1mindlinshell3d.h.
Referenced by computeStiffnessMatrix(), computeVectorOfUnknowns(), and giveInternalForcesVector().
|
staticprotected |
Definition at line 81 of file quad1mindlinshell3d.h.
Referenced by computeBmatrixAt(), computeBodyLoadVectorAt(), computeEdgeVolumeAround(), computeLoadLEToLRotationMatrix(), computeStiffnessMatrix(), computeVolumeAround(), giveInternalForcesVector(), giveInterpolation(), and giveInterpolation().
|
staticprotected |
Definition at line 88 of file quad1mindlinshell3d.h.
|
protected |
Cached coordinates in local c.s.,.
Definition at line 77 of file quad1mindlinshell3d.h.
Referenced by computeBodyLoadVectorAt(), computeGtoLRotationMatrix(), and computeLCS().
|
protected |
Cached nodal coordinates in local c.s.,.
Definition at line 75 of file quad1mindlinshell3d.h.
Referenced by computeBmatrixAt(), computeEdgeVolumeAround(), computeLCS(), computeVolumeAround(), and Quad1MindlinShell3D().
|
protected |
Flag controlling reduced (one - point) integration for shear.
Definition at line 79 of file quad1mindlinshell3d.h.
Referenced by computeBmatrixAt(), initializeFrom(), and Quad1MindlinShell3D().
|
staticprotected |
Ordering for the normal shell stiffness (everything but the out-of-plane rotations).
Definition at line 84 of file quad1mindlinshell3d.h.
Referenced by computeStiffnessMatrix(), computeVectorOfUnknowns(), giveInternalForcesVector(), and oofem::REGISTER_Element().