Description

Base class for most structural elements of 'beam' type.

#include <ChElementBeam.h>

Inheritance diagram for chrono::fea::ChElementBeam:
Collaboration diagram for chrono::fea::ChElementBeam:

Public Member Functions

virtual void EvaluateSectionDisplacement (const double eta, ChVector3d &u_displ, ChVector3d &u_rotaz)=0
 Gets the xyz displacement of a point on the beam line, and the rotation RxRyRz of section plane, at abscissa 'eta'. More...
 
virtual void EvaluateSectionFrame (const double eta, ChVector3d &point, ChQuaternion<> &rot)=0
 Gets the absolute xyz position of a point on the beam line, and the absolute rotation of section plane, at abscissa 'eta'. More...
 
virtual void EvaluateSectionForceTorque (const double eta, ChVector3d &Fforce, ChVector3d &Mtorque)=0
 Gets the force (traction x, shear y, shear z) and the torque (torsion on x, bending on y, on bending on z) at a section along the beam line, at abscissa 'eta'. More...
 
virtual void EvaluateSectionStrain (const double eta, ChVector3d &StrainV)=0
 Gets the axial and bending strain of the ANCF "cable" element.
 
double GetMass ()
 The full mass of the beam, (with const. section, density, etc.)
 
double GetRestLength ()
 The rest length of the bar.
 
void SetRestLength (double ml)
 Set the rest length of the bar (usually this should be automatically done when SetupInitial is called on beams element, given the current state, but one might need to override this, ex for precompressed beams etc).
 
- Public Member Functions inherited from chrono::fea::ChElementGeneric
ChKRMBlockKstiffness ()
 Access the proxy to stiffness, for sparse solver.
 
virtual void EleIntLoadResidual_F (ChVectorDynamic<> &R, const double c) override
 Add the internal forces (pasted at global nodes offsets) into a global vector R, multiplied by a scaling factor c, as R += forces * c This default implementation is SLIGHTLY INEFFICIENT.
 
virtual void EleIntLoadResidual_Mv (ChVectorDynamic<> &R, const ChVectorDynamic<> &w, const double c) override
 Add the product of element mass M by a vector w (pasted at global nodes offsets) into a global vector R, multiplied by a scaling factor c, as R += M * w * c This default implementation is VERY INEFFICIENT.
 
virtual void EleIntLoadLumpedMass_Md (ChVectorDynamic<> &Md, double &error, const double c) override
 Adds the lumped mass to a Md vector, representing a mass diagonal matrix. More...
 
virtual void EleIntLoadResidual_F_gravity (ChVectorDynamic<> &R, const ChVector3d &G_acc, const double c) override
 Add the contribution of gravity loads, multiplied by a scaling factor c, as: R += M * g * c This default implementation is VERY INEFFICIENT. More...
 
virtual void ComputeGravityForces (ChVectorDynamic<> &Fg, const ChVector3d &G_acc) override
 Compute the gravitational forces. More...
 
virtual void ComputeMmatrixGlobal (ChMatrixRef M) override
 Calculate the mass matrix, expressed in global reference. More...
 
virtual void InjectKRMMatrices (ChSystemDescriptor &descriptor) override
 Register with the given system descriptor any ChKRMBlock objects associated with this item.
 
virtual void LoadKRMMatrices (double Kfactor, double Rfactor, double Mfactor) override
 Compute and load current stiffnes (K), damping (R), and mass (M) matrices in encapsulated ChKRMBlock objects. More...
 
virtual void VariablesFbLoadInternalForces (double factor=1.) override
 Add the internal forces, expressed as nodal forces, into the encapsulated ChVariables.
 
virtual void VariablesFbIncrementMq () override
 Add M*q (internal masses multiplied current 'qb').
 
- Public Member Functions inherited from chrono::fea::ChElementBase
virtual unsigned int GetNumNodes ()=0
 Get the number of nodes used by this element.
 
virtual unsigned int GetNumCoordsPosLevel ()=0
 Get the number of coordinates in the field used by the referenced nodes. More...
 
virtual unsigned int GetNumCoordsPosLevelActive ()
 Get the actual number of active degrees of freedom. More...
 
virtual unsigned int GetNodeNumCoordsPosLevel (unsigned int n)=0
 Get the number of coordinates from the specified node that are used by this element. More...
 
virtual unsigned int GetNodeNumCoordsPosLevelActive (unsigned int n)
 Get the actual number of active coordinates from the specified node that are used by this element. More...
 
virtual std::shared_ptr< ChNodeFEAbaseGetNode (unsigned int n)=0
 Access the nth node.
 
virtual void GetStateBlock (ChVectorDynamic<> &mD)=0
 Fill the D vector with the current field values at the nodes of the element, with proper ordering. More...
 
virtual void ComputeNodalMass ()
 Compute element's nodal masses.
 
virtual void ComputeKRMmatricesGlobal (ChMatrixRef H, double Kfactor, double Rfactor=0, double Mfactor=0)=0
 Set H as the stiffness matrix K, scaled by Kfactor. More...
 
virtual void ComputeInternalForces (ChVectorDynamic<> &Fi)=0
 Compute the internal forces. More...
 
virtual void Update ()
 Update, called at least at each time step. More...
 
virtual void EleDoIntegration ()
 This is optionally implemented if there is some internal state that requires integration.
 

Protected Attributes

double mass
 
double length
 
- Protected Attributes inherited from chrono::fea::ChElementGeneric
ChKRMBlock Kmatr
 

Member Function Documentation

◆ EvaluateSectionDisplacement()

virtual void chrono::fea::ChElementBeam::EvaluateSectionDisplacement ( const double  eta,
ChVector3d u_displ,
ChVector3d u_rotaz 
)
pure virtual

Gets the xyz displacement of a point on the beam line, and the rotation RxRyRz of section plane, at abscissa 'eta'.

Note, eta=-1 at node1, eta=+1 at node2.

Implemented in chrono::fea::ChElementBeamANCF_3243, chrono::fea::ChElementBeamANCF_3333, chrono::fea::ChElementBeamTaperedTimoshenko, chrono::fea::ChElementBeamEuler, chrono::fea::ChElementBeamIGA, chrono::fea::ChElementCableANCF, and chrono::fea::ChElementBeamTaperedTimoshenkoFPM.

◆ EvaluateSectionForceTorque()

virtual void chrono::fea::ChElementBeam::EvaluateSectionForceTorque ( const double  eta,
ChVector3d Fforce,
ChVector3d Mtorque 
)
pure virtual

Gets the force (traction x, shear y, shear z) and the torque (torsion on x, bending on y, on bending on z) at a section along the beam line, at abscissa 'eta'.

Note, eta=-1 at node1, eta=+1 at node2. Results are not corotated, and are expressed in the reference system of beam.

Implemented in chrono::fea::ChElementBeamANCF_3243, chrono::fea::ChElementBeamANCF_3333, chrono::fea::ChElementBeamTaperedTimoshenko, chrono::fea::ChElementBeamIGA, chrono::fea::ChElementBeamEuler, chrono::fea::ChElementCableANCF, and chrono::fea::ChElementBeamTaperedTimoshenkoFPM.

◆ EvaluateSectionFrame()

virtual void chrono::fea::ChElementBeam::EvaluateSectionFrame ( const double  eta,
ChVector3d point,
ChQuaternion<> &  rot 
)
pure virtual

Gets the absolute xyz position of a point on the beam line, and the absolute rotation of section plane, at abscissa 'eta'.

Note, eta=-1 at node1, eta=+1 at node2.

Implemented in chrono::fea::ChElementBeamANCF_3243, chrono::fea::ChElementBeamANCF_3333, chrono::fea::ChElementBeamTaperedTimoshenko, chrono::fea::ChElementBeamIGA, chrono::fea::ChElementBeamEuler, and chrono::fea::ChElementCableANCF.


The documentation for this class was generated from the following file:
  • /builds/uwsbel/chrono/src/chrono/fea/ChElementBeam.h