OOFEM 3.0
Loading...
Searching...
No Matches
oofem::Beam2d Class Reference

#include <beam2d.h>

Inheritance diagram for oofem::Beam2d:
Collaboration diagram for oofem::Beam2d:

Public Member Functions

 Beam2d (int n, Domain *aDomain)
virtual ~Beam2d ()
void computeConsistentMassMatrix (FloatMatrix &answer, TimeStep *tStep, double &mass, const double *ipDensity=NULL) override
void computeInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep) override
void computeLumpedInitialStressMatrix (FloatMatrix &answer, TimeStep *tStep) override
void computeStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep) override
int giveLocalCoordinateSystem (FloatMatrix &answer) override
void giveInternalForcesVector (FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0) override
void giveEndForcesVector (FloatArray &answer, TimeStep *tStep)
int testElementExtension (ElementExtension ext) override
InterfacegiveInterface (InterfaceType) override
FEInterpolationgiveInterpolation () const override
FEInterpolationgiveInterpolation (DofIDItem id) const override
int computeNumberOfDofs () override
int computeNumberOfGlobalDofs () override
void giveDofManDofIDMask (int inode, IntArray &) const override
int giveNumberOfInternalDofManagers () const override
DofManagergiveInternalDofManager (int i) const override
void giveInternalDofManDofIDMask (int i, IntArray &answer) const override
void giveBoundaryLocationArray (IntArray &locationArray, const IntArray &bNodes, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) override
void giveBoundaryLocationArray (IntArray &locationArray, const IntArray &bNodes, const IntArray &dofIDMask, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) override
double computeVolumeAround (GaussPoint *gp) override
void printOutputAt (FILE *file, TimeStep *tStep) override
Element_Geometry_Type giveGeometryType () const override
const char * giveClassName () const override
const char * giveInputRecordName () const override
void initializeFrom (InputRecord &ir, int priority) override
void postInitialize () override
 Performs post initialization steps. Called after all components are created and initialized.
void drawRawGeometry (oofegGraphicContext &gc, TimeStep *tStep) override
void drawDeformedGeometry (oofegGraphicContext &gc, TimeStep *tStep, UnknownType) override
void computeStrainVectorInLayer (FloatArray &answer, const FloatArray &masterGpStrain, GaussPoint *masterGp, GaussPoint *slaveGp, TimeStep *tStep) override
int giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
Public Member Functions inherited from oofem::BeamBaseElement
 BeamBaseElement (int n, Domain *d)
virtual ~BeamBaseElement ()
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 computeLumpedMassMatrix (FloatMatrix &answer, TimeStep *tStep)
virtual void giveMassMtrxIntegrationgMask (IntArray &answer)
void computeStiffnessMatrix_withIRulesAsSubcells (FloatMatrix &answer, MatResponseMode rMode, 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)
virtual void computeStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep)
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)
StructuralCrossSectiongiveStructuralCrossSection ()
 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
Elementoperator= (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 int giveNumberOfDofs ()
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 IntArraygiveBodyLoadList () const
const IntArraygiveBoundaryLoadList () 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 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< IntegrationRulegiveBoundaryEdgeIntegrationRule (int order, int boundary)
virtual std::unique_ptr< IntegrationRulegiveBoundarySurfaceIntegrationRule (int order, int boundary)
int giveDofManagerNumber (int i) const
const IntArraygiveDofManArray () const
void addDofManager (DofManager *dMan)
DofManagergiveDofManager (int i) const
NodegiveNode (int i) const
virtual ElementSidegiveSide (int i) const
virtual const FEInterpolationgetGeometryInterpolation () const
virtual MaterialgiveMaterial ()
int giveMaterialNumber () const
CrossSectiongiveCrossSection ()
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 IntegrationRulegiveDefaultIntegrationRulePtr ()
int giveNumberOfIntegrationRules ()
virtual IntegrationRulegiveIntegrationRule (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)
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 int computeGlobalCoordinates (FloatArray &answer, const FloatArray &lcoords)
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 IntArraygivePartitionList () const
void setPartitionList (IntArray &pl)
virtual double predictRelativeComputationalCost ()
virtual double giveRelativeSelfComputationalCost ()
virtual double predictRelativeRedistributionCost ()
IntArraygiveBodyLoadArray ()
 Returns array containing load numbers of loads acting on element.
IntArraygiveBoundaryLoadArray ()
 Returns array containing load numbers of boundary loads acting on element.
void initializeFrom (InputRecord &ir, int priority) override
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 IntArraygiveSharedEdgeIDs () const
const IntArraygiveSharedSurfaceIDs () const
Public Member Functions inherited from oofem::FEMComponent
 FEMComponent (int n, Domain *d)
virtual ~FEMComponent ()=default
 Virtual destructor.
DomaingiveDomain () 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::LayeredCrossSectionInterface
 LayeredCrossSectionInterface ()
Public Member Functions inherited from oofem::Interface
 Interface ()
 Constructor.
virtual ~Interface ()

Protected Member Functions

void computeBoundaryEdgeLoadVector (FloatArray &answer, BoundaryLoad *load, int edge, CharType type, ValueModeType mode, TimeStep *tStep, bool global=true) override
void computeBmatrixAt (GaussPoint *, FloatMatrix &, int=1, int=ALL_STRAINS) override
void computeNmatrixAt (const FloatArray &iLocCoord, FloatMatrix &) override
bool computeGtoLRotationMatrix (FloatMatrix &answer) override
void computeConstitutiveMatrixAt (FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep) override
void computeStressVector (FloatArray &answer, const FloatArray &strain, GaussPoint *gp, TimeStep *tStep) override
void computeBodyLoadVectorAt (FloatArray &answer, Load *load, TimeStep *tStep, ValueModeType mode) override
double giveKappaCoeff (TimeStep *tStep)
double computeLength () override
double givePitch ()
void computeGaussPoints () override
MaterialMode giveMaterialMode () override
int giveNumberOfIPForMassMtrxIntegration () override
bool hasDofs2Condense ()
Protected Member Functions inherited from oofem::BeamBaseElement
virtual void computeLocalForceLoadVector (FloatArray &answer, TimeStep *tStep, ValueModeType mode)
virtual void computePointLoadVectorAt (FloatArray &answer, Load *load, TimeStep *tStep, ValueModeType mode, bool global=true)
virtual void giveEdgeDofMapping (IntArray &answer, int iEdge) const
virtual void giveSurfaceDofMapping (IntArray &answer, int iSurf) const
virtual double computeEdgeVolumeAround (GaussPoint *gp, int iEdge)
virtual double computeSurfaceVolumeAround (GaussPoint *gp, int iSurf)
virtual int computeLoadGToLRotationMtrx (FloatMatrix &answer)
virtual int computeLoadLEToLRotationMatrix (FloatMatrix &answer, int iEdge, GaussPoint *gp)
virtual int computeLoadLSToLRotationMatrix (FloatMatrix &answer, int iSurf, GaussPoint *gp)
void condense (FloatMatrix *stiff, FloatMatrix *mass, FloatArray *load, IntArray *what)
virtual void setupIRForMassMtrxIntegration (IntegrationRule &iRule)

Protected Attributes

double kappa
double pitch
double length
DofManagerghostNodes [2]
int numberOfCondensedDofs
 number of condensed DOFs
Protected Attributes inherited from oofem::StructuralElement
std::unique_ptr< FloatArrayinitialDisplacements
 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.
Domaindomain
 Link to domain object, useful for communicating with other FEM components.

Static Protected Attributes

static FEI2dLineLin interp_geom
static FEI2dLineHermite interp_beam
static ParamKey IPK_Beam2d_dofsToCondense
 [optional] DOFs to condense

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

Detailed Description

This class implements a 2-dimensional beam element with cubic lateral displacement, quadratic rotations, and linear longitudinal displacements and geometry. This is an exact displacement approximation for beam with no nonnodal loading.

This class is not derived from linear beam or truss element, because it does not support any material nonlinearities (if should, stiffness must be integrated)

Definition at line 62 of file beam2d.h.

Constructor & Destructor Documentation

◆ Beam2d()

◆ ~Beam2d()

oofem::Beam2d::~Beam2d ( )
virtual

Definition at line 83 of file beam2d.C.

References ghostNodes.

Member Function Documentation

◆ computeBmatrixAt()

void oofem::Beam2d::computeBmatrixAt ( GaussPoint * gp,
FloatMatrix & answer,
int lowerIndx = 1,
int upperIndx = ALL_STRAINS )
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.

Parameters
gpIntegration point for which answer is computed.
answerGeometric matrix of receiver.
lowerIndxIf 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).
upperIndxIf 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 105 of file beam2d.C.

References oofem::FloatMatrix::at(), computeLength(), oofem::FEMComponent::domain, giveKappaCoeff(), oofem::GaussPoint::giveNaturalCoordinate(), kappa, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by computeStiffnessMatrix().

◆ computeBodyLoadVectorAt()

void oofem::Beam2d::computeBodyLoadVectorAt ( FloatArray & answer,
Load * load,
TimeStep * tStep,
ValueModeType mode )
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.

Parameters
answerComputed load vector due to body load
loadBody load which contribution is computed.
tStepTime step.
modedetermines the response mode

Reimplemented from oofem::StructuralElement.

Definition at line 495 of file beam2d.C.

References oofem::CS_Area, oofem::Element::giveCrossSection(), and oofem::FloatArray::times().

◆ computeBoundaryEdgeLoadVector()

void oofem::Beam2d::computeBoundaryEdgeLoadVector ( FloatArray & answer,
BoundaryLoad * load,
int edge,
CharType type,
ValueModeType mode,
TimeStep * tStep,
bool global = true )
overrideprotectedvirtual

Computes the contribution of the given load at the given boundary edge. In general, the answer should include only relevant DOFs at the edge. The related is giveBoundaryLocationArray method, which should return corresponding code numbers..

Note
Elements which do not have an contribution should resize the vector to be empty.
Parameters
answerRequested contribution of load (in Global c.s.).
loadLoad to compute contribution from.
edgeEdge number.
typeType of the contribution.
modeDetermines mode of answer.
tStepTime step when answer is computed.
globalif true (default) then contribution is in global c.s., when false then contribution is in element local c.s.
Todo
Change the load type to "BoundaryEdgeLoad" maybe?

Reimplemented from oofem::Element.

Definition at line 449 of file beam2d.C.

References oofem::FloatArray::clear(), oofem::Element::computeGlobalCoordinates(), computeGtoLRotationMatrix(), computeLength(), oofem::StructuralElement::computeLoadGToLRotationMtrx(), computeNmatrixAt(), oofem::Load::computeValues(), oofem::BoundaryLoad::giveCoordSystMode(), oofem::Element::giveDefaultIntegrationRulePtr(), oofem::Load::giveFormulationType(), N, OOFEM_ERROR, oofem::FloatArray::plusProduct(), and oofem::FloatArray::rotatedWith().

◆ computeConsistentMassMatrix()

void oofem::Beam2d::computeConsistentMassMatrix ( FloatMatrix & answer,
TimeStep * tStep,
double & mass,
const double * ipDensity = NULL )
overridevirtual

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.

Parameters
answerMass matrix.
tStepTime step.
massTotal mass of receiver.

Reimplemented from oofem::StructuralElement.

Definition at line 577 of file beam2d.C.

References oofem::FloatMatrix::at(), computeLength(), oofem::CS_Area, oofem::Element::giveCrossSection(), giveKappaCoeff(), oofem::StructuralElement::giveStructuralCrossSection(), oofem::Element::integrationRulesArray, kappa, oofem::FloatMatrix::resize(), oofem::FloatMatrix::symmetrized(), and oofem::FloatMatrix::zero().

◆ computeConstitutiveMatrixAt()

void oofem::Beam2d::computeConstitutiveMatrixAt ( FloatMatrix & answer,
MatResponseMode rMode,
GaussPoint * gp,
TimeStep * tStep )
overrideprotectedvirtual

Computes constitutive matrix of receiver. Default implementation uses element cross section giveCharMaterialStiffnessMatrix service.

Parameters
answerConstitutive matrix.
rModeMaterial response mode of answer.
gpIntegration point for which constitutive matrix is computed.
tStepTime step.

Implements oofem::StructuralElement.

Definition at line 196 of file beam2d.C.

References oofem::StructuralElement::giveStructuralCrossSection().

Referenced by computeStiffnessMatrix(), and giveKappaCoeff().

◆ computeGaussPoints()

void oofem::Beam2d::computeGaussPoints ( )
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.

See also
IntegrationRule

Reimplemented from oofem::Element.

Definition at line 132 of file beam2d.C.

References oofem::Element::giveCrossSection(), oofem::Element::integrationRulesArray, and oofem::Element::numberOfGaussPoints.

◆ computeGtoLRotationMatrix()

bool oofem::Beam2d::computeGtoLRotationMatrix ( FloatMatrix & answer)
overrideprotectedvirtual

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.

Parameters
answerComputed rotation matrix.
Returns
Nonzero if transformation is necessary, zero otherwise.

Reimplemented from oofem::Element.

Definition at line 210 of file beam2d.C.

References oofem::FloatMatrix::at(), oofem::FloatMatrix::beProductOf(), computeNumberOfGlobalDofs(), ghostNodes, givePitch(), hasDofs2Condense(), pitch, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by computeBoundaryEdgeLoadVector().

◆ computeInitialStressMatrix()

void oofem::Beam2d::computeInitialStressMatrix ( FloatMatrix & answer,
TimeStep * tStep )
overridevirtual

Computes initial stress matrix for linear stability problem. Default implementation is not provided. Please note, that initial stress matrix depends on normal forces of element, corresponding engineering model must take this into account.

Parameters
answerComputed initial stress matrix.
tStepTime step.

Reimplemented from oofem::StructuralElement.

Definition at line 627 of file beam2d.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), computeLength(), giveEndForcesVector(), giveKappaCoeff(), kappa, oofem::min(), N, oofem::FloatMatrix::resize(), oofem::FloatMatrix::symmetrized(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().

◆ computeLength()

double oofem::Beam2d::computeLength ( )
overrideprotectedvirtual

◆ computeLumpedInitialStressMatrix()

void oofem::Beam2d::computeLumpedInitialStressMatrix ( FloatMatrix & answer,
TimeStep * tStep )
overridevirtual

Computes lumped initial stress matrix of receiver.

Parameters
answerLumped initial stress matrix.
tStepTime step.

Reimplemented from oofem::StructuralElement.

Definition at line 669 of file beam2d.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), computeLength(), giveEndForcesVector(), N, oofem::FloatMatrix::resize(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().

◆ computeNmatrixAt()

void oofem::Beam2d::computeNmatrixAt ( const FloatArray & iLocCoord,
FloatMatrix & answer )
overrideprotectedvirtual

Computes interpolation matrix for element unknowns. The order and meaning of unknowns is element dependent.

Parameters
iLocCoordLocal coordinates.
answerInterpolation matrix evaluated at gp.

Reimplemented from oofem::StructuralElement.

Definition at line 145 of file beam2d.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), computeLength(), oofem::FEMComponent::domain, giveKappaCoeff(), kappa, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by computeBoundaryEdgeLoadVector().

◆ computeNumberOfDofs()

int oofem::Beam2d::computeNumberOfDofs ( )
inlineoverridevirtual

Computes or simply returns total number of element's local DOFs. Must be defined by particular element.

Returns
Number of local DOFs of element.

Reimplemented from oofem::Element.

Definition at line 101 of file beam2d.h.

◆ computeNumberOfGlobalDofs()

int oofem::Beam2d::computeNumberOfGlobalDofs ( )
inlineoverridevirtual

Computes the total number of element's global dofs. The transitions from global c.s. to nodal c.s. should NOT be included.

Returns
Total number of global DOFs of element.

Reimplemented from oofem::Element.

Definition at line 102 of file beam2d.h.

Referenced by computeGtoLRotationMatrix().

◆ computeStiffnessMatrix()

void oofem::Beam2d::computeStiffnessMatrix ( FloatMatrix & answer,
MatResponseMode rMode,
TimeStep * tStep )
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).

Parameters
answerComputed stiffness matrix (symmetric).
rModeResponse mode.
tStepTime step.

Reimplemented from oofem::StructuralElement.

Definition at line 179 of file beam2d.C.

References oofem::FloatMatrix::beProductOf(), oofem::FloatMatrix::clear(), computeBmatrixAt(), computeConstitutiveMatrixAt(), computeLength(), oofem::Element::giveDefaultIntegrationRulePtr(), oofem::FloatMatrix::plusProductSymmUpper(), and oofem::FloatMatrix::symmetrized().

◆ computeStrainVectorInLayer()

void oofem::Beam2d::computeStrainVectorInLayer ( FloatArray & answer,
const FloatArray & masterGpStrain,
GaussPoint * masterGp,
GaussPoint * slaveGp,
TimeStep * tStep )
overridevirtual

Computes full 3D strain vector in element layer. This function is necessary if layered cross section is specified..

Parameters
answerFull layer strain vector.
masterGpStrainGeneralized strain at master gauss point.
masterGpElement integration point.
slaveGpSlave integration point representing particular layer.
tStepTime step.

Implements oofem::LayeredCrossSectionInterface.

Definition at line 274 of file beam2d.C.

References oofem::FloatArray::at(), oofem::CS_BottomZCoord, oofem::CS_TopZCoord, oofem::Element::giveCrossSection(), oofem::GaussPoint::giveNaturalCoordinate(), and oofem::FloatArray::resize().

◆ computeStressVector()

void oofem::Beam2d::computeStressVector ( FloatArray & answer,
const FloatArray & strain,
GaussPoint * gp,
TimeStep * tStep )
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.

Parameters
answerStress vector.
strainStrain vector.
gpIntegration point.
tStepTime step.

Implements oofem::StructuralElement.

Definition at line 203 of file beam2d.C.

References oofem::StructuralElement::giveStructuralCrossSection().

◆ computeVolumeAround()

double oofem::Beam2d::computeVolumeAround ( GaussPoint * gp)
overridevirtual

Returns volume related to given integration point. Used typically in subroutines, that perform integration over element volume. Should be implemented by particular elements.

See also
GaussPoint
Parameters
gpIntegration point for which volume is computed.
Returns
Volume for integration point.

Reimplemented from oofem::Element.

Definition at line 267 of file beam2d.C.

References computeLength(), and oofem::GaussPoint::giveWeight().

◆ drawDeformedGeometry()

void oofem::Beam2d::drawDeformedGeometry ( oofegGraphicContext & gc,
TimeStep * tStep,
UnknownType type )
overridevirtual

Reimplemented from oofem::Element.

Definition at line 722 of file beam2d.C.

References gc, oofem::Element::giveNode(), OOFEG_DEFORMED_GEOMETRY_LAYER, and OOFEG_DEFORMED_GEOMETRY_WIDTH.

◆ drawRawGeometry()

void oofem::Beam2d::drawRawGeometry ( oofegGraphicContext & gc,
TimeStep * tStep )
overridevirtual

Reimplemented from oofem::Element.

Definition at line 696 of file beam2d.C.

References gc, oofem::Element::giveNode(), OOFEG_RAW_GEOMETRY_LAYER, and OOFEG_RAW_GEOMETRY_WIDTH.

◆ giveBoundaryLocationArray() [1/2]

void oofem::Beam2d::giveBoundaryLocationArray ( IntArray & locationArray,
const IntArray & bNodes,
const IntArray & dofIDMask,
const UnknownNumberingScheme & s,
IntArray * dofIds = NULL )
inlineoverridevirtual

Reimplemented from oofem::Element.

Definition at line 132 of file beam2d.h.

References oofem::Element::giveLocationArray().

◆ giveBoundaryLocationArray() [2/2]

void oofem::Beam2d::giveBoundaryLocationArray ( IntArray & locationArray,
const IntArray & bNodes,
const UnknownNumberingScheme & s,
IntArray * dofIds = NULL )
inlineoverridevirtual

Returns the location array for the boundary of the element. Only takes into account nodes in the bNodes vector.

Reimplemented from oofem::Element.

Definition at line 128 of file beam2d.h.

References oofem::Element::giveLocationArray().

◆ giveClassName()

const char * oofem::Beam2d::giveClassName ( ) const
inlineoverridevirtual
Returns
Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 141 of file beam2d.h.

◆ giveDofManDofIDMask()

void oofem::Beam2d::giveDofManDofIDMask ( int inode,
IntArray & answer ) const
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.

Parameters
inodeMask is computed for local dofmanager with inode number.
answerMask for node.

Reimplemented from oofem::Element.

Definition at line 295 of file beam2d.C.

◆ giveEndForcesVector()

◆ giveGeometryType()

Element_Geometry_Type oofem::Beam2d::giveGeometryType ( ) const
inlineoverridevirtual

Returns the element geometry type. This information is assumed to be of general interest, but it is required only for some specialized tasks.

Returns
Geometry type of element.

Implements oofem::Element.

Definition at line 139 of file beam2d.h.

◆ giveInputRecordName()

const char * oofem::Beam2d::giveInputRecordName ( ) const
inlineoverridevirtual
Returns
Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 142 of file beam2d.h.

References _IFT_Beam2d_Name.

◆ giveInterface()

Interface * oofem::Beam2d::giveInterface ( InterfaceType t)
overridevirtual

Interface requesting service.

See also
InterfaceType
Returns
Requested interface if implemented, otherwise NULL.

Reimplemented from oofem::FEMComponent.

Definition at line 94 of file beam2d.C.

References oofem::LayeredCrossSectionInterface::LayeredCrossSectionInterface(), and oofem::LayeredCrossSectionInterfaceType.

◆ giveInternalDofManager()

DofManager * oofem::Beam2d::giveInternalDofManager ( int i) const
inlineoverridevirtual

Returns i-th internal element dof manager of the receiver

Parameters
iInternal number of DOF.
Returns
DOF number i.

Reimplemented from oofem::Element.

Definition at line 105 of file beam2d.h.

References ghostNodes, oofem::FEMComponent::number, and OOFEM_ERROR.

◆ giveInternalDofManDofIDMask()

void oofem::Beam2d::giveInternalDofManDofIDMask ( int inode,
IntArray & answer ) const
inlineoverridevirtual

Returns internal 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.

Parameters
inodeMask is computed for local dofmanager with inode number.
answermask for node.

Reimplemented from oofem::Element.

Definition at line 114 of file beam2d.h.

References ghostNodes, oofem::FEMComponent::number, and OOFEM_ERROR.

◆ giveInternalForcesVector()

void oofem::Beam2d::giveInternalForcesVector ( FloatArray & answer,
TimeStep * tStep,
int useUpdatedGpRecord = 0 )
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.

Parameters
answerInternal nodal forces vector.
tStepTime step.
useUpdatedGpRecordIf 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 426 of file beam2d.C.

Referenced by giveEndForcesVector().

◆ giveInterpolation() [1/2]

FEInterpolation * oofem::Beam2d::giveInterpolation ( ) const
overridevirtual
Returns
Interpolation of the element geometry, or NULL if none exist.

Reimplemented from oofem::Element.

Definition at line 90 of file beam2d.C.

References interp_geom.

◆ giveInterpolation() [2/2]

FEInterpolation * oofem::Beam2d::giveInterpolation ( DofIDItem id) const
inlineoverridevirtual

Returns the interpolation for the specific dof id. Special elements which uses a mixed interpolation should reimplement this method.

Parameters
idID of the dof for the for the requested interpolation.
Returns
Appropriate interpolation, or NULL if none exists.

Reimplemented from oofem::Element.

Definition at line 99 of file beam2d.h.

◆ giveIPValue()

int oofem::Beam2d::giveIPValue ( FloatArray & answer,
GaussPoint * gp,
InternalStateType type,
TimeStep * tStep )
overridevirtual

Returns the integration point corresponding value in full form.

Parameters
answerContain corresponding integration point value, zero sized if not available.
gpIntegration point to check.
typeDetermines the type of internal variable.
tStepTime step.
Returns
Nonzero if o.k, zero otherwise.
Todo
Which "error type" should be used? Why are there several? I don't see the point of this enum when there could be different function calls just as well (and different IST values)

Reimplemented from oofem::Element.

Definition at line 504 of file beam2d.C.

References oofem::FloatArray::at(), oofem::GaussPoint::giveMaterialStatus(), oofem::StructuralMaterialStatus::giveStrainVector(), oofem::StructuralMaterialStatus::giveStressVector(), and oofem::FloatArray::resize().

◆ giveKappaCoeff()

◆ giveLocalCoordinateSystem()

int oofem::Beam2d::giveLocalCoordinateSystem ( FloatMatrix & answer)
overridevirtual

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.

Returns
nonzero if answer computed, zero value if answer is empty, i.e. no transformation is necessary.

Reimplemented from oofem::Element.

Definition at line 359 of file beam2d.C.

References oofem::FloatMatrix::at(), givePitch(), pitch, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

◆ giveMaterialMode()

MaterialMode oofem::Beam2d::giveMaterialMode ( )
inlineoverrideprotectedvirtual

Returns material mode for receiver integration points. Should be specialized.

Returns
Material mode of element.

Reimplemented from oofem::Element.

Definition at line 171 of file beam2d.h.

◆ giveNumberOfInternalDofManagers()

int oofem::Beam2d::giveNumberOfInternalDofManagers ( ) const
inlineoverridevirtual
Returns
Number of internal DOF managers of element.

Reimplemented from oofem::Element.

Definition at line 104 of file beam2d.h.

References ghostNodes.

◆ giveNumberOfIPForMassMtrxIntegration()

int oofem::Beam2d::giveNumberOfIPForMassMtrxIntegration ( )
inlineoverrideprotectedvirtual

Return desired number of integration points for consistent mass matrix computation, if required.

Returns
Number of integration points for mass matrix.

TODO this is without the jacobian and density

Reimplemented from oofem::StructuralElement.

Definition at line 172 of file beam2d.h.

◆ givePitch()

double oofem::Beam2d::givePitch ( )
protected

◆ hasDofs2Condense()

bool oofem::Beam2d::hasDofs2Condense ( )
inlineprotected

Definition at line 174 of file beam2d.h.

References ghostNodes.

Referenced by computeGtoLRotationMatrix().

◆ initializeFrom()

void oofem::Beam2d::initializeFrom ( InputRecord & ir,
int priority )
overridevirtual

◆ postInitialize()

void oofem::Beam2d::postInitialize ( )
overridevirtual

◆ printOutputAt()

void oofem::Beam2d::printOutputAt ( FILE * file,
TimeStep * tStep )
overridevirtual

Prints output of receiver to stream, for given time step. This is used for output into the standard output file.

Parameters
fileFile pointer to print to.
tStepTime step to write for.

Reimplemented from oofem::FEMComponent.

Definition at line 546 of file beam2d.C.

References oofem::Element::computeVectorOf(), giveEndForcesVector(), oofem::Element::giveLabel(), oofem::FEMComponent::giveNumber(), and oofem::Element::integrationRulesArray.

◆ testElementExtension()

int oofem::Beam2d::testElementExtension ( ElementExtension ext)
inlineoverridevirtual

Tests if the element implements required extension. ElementExtension type defines the list of all available element extensions.

Parameters
extExtension to be tested.
Returns
Nonzero if extension supported.
See also
ElementExtension

Reimplemented from oofem::Element.

Definition at line 94 of file beam2d.h.

References oofem::Element_EdgeLoadSupport.

Member Data Documentation

◆ ghostNodes

DofManager* oofem::Beam2d::ghostNodes[2]
protected

Ghost nodes are used to introduce additional DOFs at element. These are needed as we actually do not want to condense selected DOFs, but rather allocate an extra equation to these. This allows to get cooresponding DOFs directly from the global system, avoiding the need to postprocess local displacements at element.

Definition at line 72 of file beam2d.h.

Referenced by Beam2d(), computeGtoLRotationMatrix(), giveInternalDofManager(), giveInternalDofManDofIDMask(), giveNumberOfInternalDofManagers(), hasDofs2Condense(), postInitialize(), and ~Beam2d().

◆ interp_beam

FEI2dLineHermite oofem::Beam2d::interp_beam
staticprotected

Definition at line 77 of file beam2d.h.

◆ interp_geom

FEI2dLineLin oofem::Beam2d::interp_geom
staticprotected

Definition at line 76 of file beam2d.h.

Referenced by giveInterpolation().

◆ IPK_Beam2d_dofsToCondense

ParamKey oofem::Beam2d::IPK_Beam2d_dofsToCondense
staticprotected

[optional] DOFs to condense

Definition at line 79 of file beam2d.h.

Referenced by initializeFrom(), and postInitialize().

◆ kappa

double oofem::Beam2d::kappa
protected

◆ length

double oofem::Beam2d::length
protected

Definition at line 65 of file beam2d.h.

Referenced by Beam2d(), and computeLength().

◆ numberOfCondensedDofs

int oofem::Beam2d::numberOfCondensedDofs
protected

number of condensed DOFs

Definition at line 74 of file beam2d.h.

Referenced by Beam2d(), and postInitialize().

◆ pitch

double oofem::Beam2d::pitch
protected

Definition at line 65 of file beam2d.h.

Referenced by Beam2d(), computeGtoLRotationMatrix(), giveLocalCoordinateSystem(), and givePitch().


The documentation for this class was generated from the following files:

This page is part of the OOFEM-3.0 documentation. Copyright Copyright (C) 1994-2025 Borek Patzak Bořek Patzák
Project e-mail: oofem@fsv.cvut.cz
Generated at for OOFEM by doxygen 1.15.0 written by Dimitri van Heesch, © 1997-2011