|
OOFEM 3.0
|
#include <libeam3dnl2.h>
Public Member Functions | |
| LIBeam3dNL2 (int n, Domain *d) | |
| virtual | ~LIBeam3dNL2 () |
| void | computeLumpedMassMatrix (FloatMatrix &answer, TimeStep *tStep) override |
| void | computeConsistentMassMatrix (FloatMatrix &answer, TimeStep *tStep, double &mass, const double *ipDensity=NULL) override |
| void | computeStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) override |
| int | computeNumberOfDofs () override |
| void | giveDofManDofIDMask (int inode, IntArray &answer) const override |
| double | computeVolumeAround (GaussPoint *gp) override |
| int | computeGlobalCoordinates (FloatArray &answer, const FloatArray &lcoords) override |
| const char * | giveInputRecordName () const override |
| const char * | giveClassName () const override |
| void | initializeFrom (InputRecord &ir, int prio) override |
| void | postInitialize () override |
| Performs post initialization steps. | |
| Element_Geometry_Type | giveGeometryType () const override |
| void | drawRawGeometry (oofegGraphicContext &gc, TimeStep *tStep) override |
| void | drawDeformedGeometry (oofegGraphicContext &gc, TimeStep *tStep, UnknownType) override |
| void | computeStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) override |
| void | giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override |
| integrationDomain | giveIntegrationDomain () const override |
| MaterialMode | giveMaterialMode () override |
| void | saveContext (DataStream &stream, ContextMode mode) override |
| void | restoreContext (DataStream &stream, ContextMode mode) override |
| Public Member Functions inherited from oofem::NLStructuralElement | |
| NLStructuralElement (int n, Domain *d) | |
| virtual | ~NLStructuralElement () |
| Destructor. | |
| int | giveGeometryMode () |
| void | computeFirstPKStressVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) |
| void | computeCauchyStressVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) |
| void | computeStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) override |
| void | computeInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep) override |
| void | computeStiffnessMatrix_withIRulesAsSubcells (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) |
| void | giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override |
| void | giveInternalForcesVector_withIRulesAsSubcells (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override |
| virtual void | computeDeformationGradientVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) |
| double | computeCurrentVolume (TimeStep *tStep) |
| void | initializeFrom (InputRecord &ir, int priority) override |
| void | giveInputRecord (DynamicInputRecord &input) override |
| const char * | giveClassName () const 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 | computeMassMatrix (FloatMatrix &answer, TimeStep *tStep) |
| virtual void | giveMassMtrxIntegrationgMask (IntArray &answer) |
| void | computeStiffnessMatrix_withIRulesAsSubcells (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) |
| virtual void | computeLumpedInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep) |
| void | computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer) override |
| 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) |
| 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) |
| virtual int | computeNumberOfGlobalDofs () |
| int | computeNumberOfPrimaryMasterDofs () |
| virtual bool | computeGtoLRotationMatrix (FloatMatrix &answer) |
| 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 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 FEInterpolation * | giveInterpolation () const |
| virtual FEInterpolation * | giveInterpolation (DofIDItem id) 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 int | giveIntegrationRuleLocalCodeNumbers (IntArray &answer, IntegrationRule &ie) |
| int | giveRegionNumber () |
| virtual void | initializeYourself (TimeStep *timeStepWhenICApply) |
| virtual bool | isActivated (TimeStep *tStep) |
| virtual bool | isCast (TimeStep *tStep) |
| 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 () |
| virtual int | testElementExtension (ElementExtension ext) |
| int | giveGlobalIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) |
| virtual double | giveLengthInDir (const FloatArray &normalToCrackPlane) |
| virtual double | giveCharacteristicLength (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 bool | computeLocalCoordinates (FloatArray &answer, const FloatArray &gcoords) |
| virtual void | giveLocalCoordinateSystemVector (InternalStateType isttype, FloatArray &answer) |
| virtual void | computeMidPlaneNormal (FloatArray &answer, const GaussPoint *gp) |
| 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 | 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 | 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. | |
| virtual Interface * | giveInterface (InterfaceType t) |
| std::string | errorInfo (const char *func) const |
| Returns string for prepending output (used by error reporting macros). | |
Private Attributes | |
| double | l0 |
| Initial length. | |
| FloatArray | q |
| Quaternion at the center (last equilibrated). | |
| FloatArray | tempQ |
| Temporary quaternion at the center. | |
| StateCounterType | tempQCounter |
| Time stamp of temporary centre quaternion. | |
| int | referenceNode |
| Reference node. | |
Static Private Attributes | |
| static ParamKey | IPK_LIBeam3dNL2_refnode |
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 |
| Protected Attributes inherited from oofem::NLStructuralElement | |
| int | nlGeometry =0 |
| Flag indicating if geometrical nonlinearities apply. | |
| 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 inherited from oofem::NLStructuralElement | |
| static ParamKey | IPK_NLStructuralElement_nlgeoflag |
This class implements a 3-dimensional Linear Isoparametric Mindlin theory beam element, with reduced integration. Geometric nonlinearities are taken into account. Based on Element due to Simo and Vu-Quoc, description taken from Crisfield monograph. Similar to Libeam3dNL, but rotational update is done using quaternions.
Definition at line 56 of file libeam3dnl2.h.
| oofem::LIBeam3dNL2::LIBeam3dNL2 | ( | int | n, |
| Domain * | d ) |
Definition at line 60 of file libeam3dnl2.C.
References l0, oofem::NLStructuralElement::NLStructuralElement(), oofem::Element::numberOfDofMans, q, referenceNode, tempQ, and tempQCounter.
|
inlinevirtual |
Definition at line 76 of file libeam3dnl2.h.
|
inlineoverrideprotectedvirtual |
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 120 of file libeam3dnl2.h.
References OOFEM_ERROR.
|
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 699 of file libeam3dnl2.C.
References oofem::StructuralElement::computeBodyLoadVectorAt(), oofem::CS_Area, oofem::Element::giveCrossSection(), and oofem::FloatArray::times().
|
inlineoverridevirtual |
Computes consistent mass matrix of receiver using numerical integration over element volume. Mass matrix is computed as \( M = \int_V N^{\mathrm{T}} \rho N dV \), where \( N \) is displacement approximation matrix. The number of necessary integration points is determined using this->giveNumberOfIPForMassMtrxIntegration service. Only selected degrees of freedom participate in integration of mass matrix. This is described using dof mass integration mask. This mask is obtained from this->giveMassMtrxIntegrationgMask service. The nonzero mask value at i-th position indicates that i-th element DOF participates in mass matrix computation. The result is in element local coordinate system.
| answer | Mass matrix. |
| tStep | Time step. |
| mass | Total mass of receiver. |
Reimplemented from oofem::StructuralElement.
Definition at line 79 of file libeam3dnl2.h.
References computeLumpedMassMatrix().
|
overrideprotectedvirtual |
Computes large strain 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::NLStructuralElement.
Definition at line 443 of file libeam3dnl2.C.
References OOFEM_ERROR.
|
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 437 of file libeam3dnl2.C.
References oofem::StructuralCrossSection::give3dBeamStiffMtrx(), 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 606 of file libeam3dnl2.C.
References computeLength(), oofem::GaussPoint::giveWeight(), and OOFEM_ERROR.
|
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 425 of file libeam3dnl2.C.
References oofem::Element::giveCrossSection(), oofem::Element::integrationRulesArray, and oofem::CrossSection::setupIntegrationPoints().
|
overridevirtual |
Computes the global coordinates from given element's local coordinates.
| answer | Requested global coordinates. |
| lcoords | Local coordinates. |
Reimplemented from oofem::Element.
Definition at line 570 of file libeam3dnl2.C.
References oofem::FloatArray::at(), oofem::DofManager::giveCoordinate(), oofem::Element::giveNode(), and oofem::FloatArray::resize().
|
overrideprotectedvirtual |
Computes the length (zero for all but 1D geometries)
Reimplemented from oofem::Element.
Definition at line 484 of file libeam3dnl2.C.
References oofem::DofManager::giveCoordinate(), oofem::Element::giveNode(), and l0.
Referenced by computeEdgeVolumeAround(), computeLumpedMassMatrix(), computeVolumeAround(), and giveLocalCoordinateSystem().
|
overrideprotectedvirtual |
Returns transformation matrix from global coordinate system to local element coordinate system for element load vector components. If no transformation is necessary, answer is empty matrix (default);
| answer | Transformation matrix. |
Reimplemented from oofem::StructuralElement.
Definition at line 655 of file libeam3dnl2.C.
References oofem::FloatMatrix::at(), giveLocalCoordinateSystem(), oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().
|
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 684 of file libeam3dnl2.C.
References oofem::FloatMatrix::clear().
|
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 504 of file libeam3dnl2.C.
References oofem::FloatMatrix::at(), computeLength(), oofem::CS_Area, oofem::CrossSection::give(), oofem::Element::giveCrossSection(), oofem::StructuralElement::giveStructuralCrossSection(), oofem::Element::integrationRulesArray, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().
Referenced by computeConsistentMassMatrix().
|
overrideprotectedvirtual |
Computes interpolation matrix for element unknowns. The order and meaning of unknowns is element dependent.
| iLocCoord | Local coordinates. |
| answer | Interpolation matrix evaluated at gp. |
Reimplemented from oofem::StructuralElement.
Definition at line 519 of file libeam3dnl2.C.
References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().
|
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 83 of file libeam3dnl2.h.
|
protected |
Computes the normalized quaternion from the given rotation matrix.
| answer | Computed quaternion. |
| R | Source rotation matrix. |
Definition at line 182 of file libeam3dnl2.C.
References oofem::FloatArray::at(), oofem::FloatMatrix::at(), and oofem::FloatArray::resize().
Referenced by postInitialize().
|
protected |
Evaluates the rotation matrix for large rotations according to Rodrigues formula for given pseudovector psi.
| answer | Result. |
| psi | Pseudovector. |
Definition at line 91 of file libeam3dnl2.C.
References oofem::FloatMatrix::add(), oofem::FloatMatrix::at(), oofem::FloatMatrix::beProductOf(), oofem::FloatArray::computeNorm(), computeSMtrx(), oofem::FloatArray::giveSize(), OOFEM_ERROR, oofem::FloatMatrix::resize(), S, oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().
Referenced by computeTempCurv().
|
protected |
Computes rotation matrix from given quaternion.
| answer | Returned rotation matrix. |
| q | Input quaternion |
Definition at line 161 of file libeam3dnl2.C.
References oofem::FloatMatrix::at(), q, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::times().
Referenced by computeStiffnessMatrix(), computeStrainVector(), computeTempCurv(), drawDeformedGeometry(), and giveInternalForcesVector().
|
protected |
Evaluates the S matrix from given vector vec.
| answer | Assembled result. |
| vec | Source vector. |
Definition at line 72 of file libeam3dnl2.C.
References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatArray::giveSize(), OOFEM_ERROR, and oofem::FloatMatrix::resize().
Referenced by computeRotMtrx(), computeStiffnessMatrix(), and computeXMtrx().
|
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 338 of file libeam3dnl2.C.
References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatMatrix::beProductOf(), oofem::FloatMatrix::beProductTOf(), oofem::FloatMatrix::clear(), computeConstitutiveMatrixAt(), computeRotMtrxFromQuaternion(), computeSMtrx(), computeStrainVector(), computeStressVector(), computeXdVector(), computeXMtrx(), oofem::IntegrationRule::getIntegrationPoint(), oofem::Element::giveDefaultIntegrationRulePtr(), l0, tempQ, oofem::FloatMatrix::times(), and updateTempQuaternion().
|
overridevirtual |
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 230 of file libeam3dnl2.C.
References oofem::FloatArray::at(), computeRotMtrxFromQuaternion(), computeTempCurv(), computeXdVector(), l0, oofem::FloatArray::resize(), tempQ, and updateTempQuaternion().
Referenced by computeStiffnessMatrix(), and giveInternalForcesVector().
|
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 450 of file libeam3dnl2.C.
References oofem::StructuralCrossSection::giveGeneralizedStress_Beam3d(), and oofem::StructuralElement::giveStructuralCrossSection().
Referenced by computeStiffnessMatrix(), and giveInternalForcesVector().
|
protected |
Compute the temporary curvature at the centre to the state identified by given solution step.
| answer | Computed curvature. |
| tStep | Solution step identifying reached state. |
Definition at line 734 of file libeam3dnl2.C.
References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatArray::beProductOf(), oofem::FloatMatrix::beProductOf(), oofem::FloatArray::beTProductOf(), oofem::FloatArray::beVectorProductOf(), computeRotMtrx(), computeRotMtrxFromQuaternion(), oofem::FloatArray::computeSquaredNorm(), oofem::Element::computeVectorOf(), oofem::IntegrationRule::getIntegrationPoint(), oofem::Element::giveDefaultIntegrationRulePtr(), oofem::GaussPoint::giveMaterialStatus(), oofem::FloatArray::giveSize(), oofem::StructuralMaterialStatus::giveStrainVector(), l0, q, oofem::FloatArray::resize(), and oofem::FloatArray::times().
Referenced by computeStrainVector().
|
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 554 of file libeam3dnl2.C.
References computeLength(), and oofem::GaussPoint::giveWeight().
|
protected |
Computes x_21' vector for given solution state.
| answer | Returned x_21'. |
| tStep | Determines solution state. |
Definition at line 320 of file libeam3dnl2.C.
References oofem::FloatArray::at(), oofem::Element::computeVectorOf(), oofem::DofManager::giveCoordinate(), oofem::Element::giveNode(), and oofem::FloatArray::resize().
Referenced by computeStiffnessMatrix(), computeStrainVector(), and computeXMtrx().
|
protected |
Computes X matrix at given solution state.
| answer | Returned X matrix. |
| tStep | Determines solution state. |
Definition at line 260 of file libeam3dnl2.C.
References oofem::FloatMatrix::at(), computeSMtrx(), computeXdVector(), oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().
Referenced by computeStiffnessMatrix(), and giveInternalForcesVector().
|
overridevirtual |
Reimplemented from oofem::Element.
Definition at line 891 of file libeam3dnl2.C.
References oofem::FloatMatrix::at(), computeRotMtrxFromQuaternion(), gc, oofem::Element::giveNode(), oofem::Node::giveUpdatedCoordinate(), l0, OOFEG_DEFORMED_GEOMETRY_LAYER, OOFEG_DEFORMED_GEOMETRY_WIDTH, and q.
|
overridevirtual |
Reimplemented from oofem::Element.
Definition at line 865 of file libeam3dnl2.C.
References gc, oofem::DofManager::giveCoordinate(), oofem::Element::giveNode(), OOFEG_RAW_GEOMETRY_LAYER, and OOFEG_RAW_GEOMETRY_WIDTH.
|
inlineoverridevirtual |
Reimplemented from oofem::Element.
Definition at line 91 of file libeam3dnl2.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 564 of file libeam3dnl2.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 588 of file libeam3dnl2.C.
References oofem::IntArray::at(), OOFEM_ERROR, and oofem::IntArray::resize().
|
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 94 of file libeam3dnl2.h.
|
inlineoverridevirtual |
Implements oofem::FEMComponent.
Definition at line 90 of file libeam3dnl2.h.
References _IFT_LIBeam3dNL2_Name.
|
inlineoverridevirtual |
Returns integration domain for receiver, used to initialize integration point over receiver volume. Default behavior is taken from the default interpolation.
Reimplemented from oofem::Element.
Definition at line 104 of file libeam3dnl2.h.
References oofem::_Line.
|
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 285 of file libeam3dnl2.C.
References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatArray::beProductOf(), computeRotMtrxFromQuaternion(), computeStrainVector(), computeStressVector(), computeXMtrx(), oofem::IntegrationRule::getIntegrationPoint(), oofem::Element::giveDefaultIntegrationRulePtr(), oofem::GaussPoint::giveMaterialStatus(), oofem::StructuralMaterialStatus::giveStressVector(), tempQ, and updateTempQuaternion().
|
overrideprotectedvirtual |
Returns local coordinate system of receiver Required by material models with ortho- and anisotrophy. Returns a unit vectors of local coordinate system at element stored row-wise. If local system is equal to global one, set answer to empty matrix and return zero value.
Reimplemented from oofem::Element.
Definition at line 618 of file libeam3dnl2.C.
References oofem::FloatArray::at(), oofem::FloatMatrix::at(), computeLength(), oofem::DofManager::giveCoordinate(), oofem::FEMComponent::giveDomain(), oofem::Domain::giveNode(), oofem::Element::giveNode(), length(), referenceNode, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().
Referenced by computeLoadGToLRotationMtrx(), and postInitialize().
|
inlineoverridevirtual |
Returns material mode for receiver integration points. Should be specialized.
Reimplemented from oofem::Element.
Definition at line 105 of file libeam3dnl2.h.
|
overrideprotectedvirtual |
Initializes receivers state to new time step. It can be used also if current time step must be re-started. Default implementation invokes initForNewStep member function for all defined integrationRules. Thus all state variables in all defined integration points are re initialized.
Reimplemented from oofem::Element.
Definition at line 724 of file libeam3dnl2.C.
References oofem::Element::initForNewStep(), q, and tempQ.
|
overridevirtual |
Reimplemented from oofem::Element.
Definition at line 457 of file libeam3dnl2.C.
References oofem::FEMComponent::domain, oofem::NLStructuralElement::initializeFrom(), IPK_LIBeam3dNL2_refnode, oofem::FEMComponent::number, PM_UPDATE_PARAMETER, and referenceNode.
|
overridevirtual |
Performs post initialization steps.
Reimplemented from oofem::Element.
Definition at line 466 of file libeam3dnl2.C.
References oofem::FloatMatrix::beTranspositionOf(), computeQuaternionFromRotMtrx(), oofem::FEMComponent::domain, giveLocalCoordinateSystem(), IPK_LIBeam3dNL2_refnode, oofem::NLStructuralElement::nlGeometry, oofem::FEMComponent::number, OOFEM_ERROR, PM_ELEMENT_ERROR_IFNOTSET, oofem::Element::postInitialize(), q, and referenceNode.
|
overridevirtual |
Restores the receiver state previously written in stream.
| stream | Input stream. |
| mode | Determines amount of info available in stream (state, definition, ...). |
| throws | an ContextIOERR exception if error encountered. |
Reimplemented from oofem::Element.
Definition at line 853 of file libeam3dnl2.C.
References oofem::CIO_OK, q, oofem::Element::restoreContext(), and THROW_CIOERR.
|
overridevirtual |
Stores receiver state to output stream.
| stream | Output stream. |
| mode | Determines amount of info required in stream (state, definition, ...). |
| throws | an ContextIOERR exception if error encountered. |
Reimplemented from oofem::Element.
Definition at line 841 of file libeam3dnl2.C.
References oofem::CIO_OK, q, oofem::Element::saveContext(), and THROW_CIOERR.
|
protected |
Updates the temporary triad at the centre to the state identified by given solution step. The attribute tempQ is changed to reflect new state and tempQCounter is set to solution step counter to avoid multiple updates.
| tStep | Solution step identifying reached state. |
Definition at line 121 of file libeam3dnl2.C.
References oofem::FloatArray::at(), oofem::Element::computeVectorOf(), oofem::TimeStep::giveSolutionStateCounter(), q, tempQ, and tempQCounter.
Referenced by computeStiffnessMatrix(), computeStrainVector(), giveInternalForcesVector(), and updateYourself().
|
overrideprotectedvirtual |
Updates element state after equilibrium in time step has been reached. Default implementation updates all integration rules defined by integrationRulesArray member variable. Doing this, all integration points and their material statuses are updated also. All temporary history variables, which now describe equilibrium state are copied into equilibrium ones. The existing internal state is used for update.
| tStep | Time step for newly reached state. |
Reimplemented from oofem::Element.
Definition at line 708 of file libeam3dnl2.C.
References q, tempQ, updateTempQuaternion(), and oofem::Element::updateYourself().
|
staticprivate |
Definition at line 72 of file libeam3dnl2.h.
Referenced by initializeFrom(), and postInitialize().
|
private |
Initial length.
Definition at line 60 of file libeam3dnl2.h.
Referenced by computeLength(), computeStiffnessMatrix(), computeStrainVector(), computeTempCurv(), drawDeformedGeometry(), and LIBeam3dNL2().
|
private |
Quaternion at the center (last equilibrated).
Definition at line 62 of file libeam3dnl2.h.
Referenced by computeRotMtrxFromQuaternion(), computeTempCurv(), drawDeformedGeometry(), initForNewStep(), LIBeam3dNL2(), postInitialize(), restoreContext(), saveContext(), updateTempQuaternion(), and updateYourself().
|
private |
Reference node.
Definition at line 70 of file libeam3dnl2.h.
Referenced by giveLocalCoordinateSystem(), initializeFrom(), LIBeam3dNL2(), and postInitialize().
|
private |
Temporary quaternion at the center.
Definition at line 64 of file libeam3dnl2.h.
Referenced by computeStiffnessMatrix(), computeStrainVector(), giveInternalForcesVector(), initForNewStep(), LIBeam3dNL2(), updateTempQuaternion(), and updateYourself().
|
private |
Time stamp of temporary centre quaternion.
Definition at line 68 of file libeam3dnl2.h.
Referenced by LIBeam3dNL2(), and updateTempQuaternion().