chrono::fea::ChLoadContactSurfaceMesh Class Reference

Description

Class for applyings loads to a contact mesh as a cluster of forces operating on the nodes of the underlying finite elements.

It is useful for doing cosimulation: one can pass this object's vertex & faces to an external software (ex. CFD) that in turn will perform collision detection with its entities, compute forces, send forces back to Chrono via this object. Note, this is based on a cluster of std::vector< std::shared_ptr<ChLoadXYZnode> >, but the class itself could bypass all methods of ChLoadXYZnode and directly implement a more efficient LoadIntLoadResidual_F, however this is left in this way for didactical reasons.

#include <ChLoadContactSurfaceMesh.h>

Inheritance diagram for chrono::fea::ChLoadContactSurfaceMesh:
Collaboration diagram for chrono::fea::ChLoadContactSurfaceMesh:

Public Member Functions

 ChLoadContactSurfaceMesh (std::shared_ptr< ChContactSurfaceMesh > cmesh)
 
void OutputSimpleMesh (std::vector< ChVector<>> &vert_pos, std::vector< ChVector<>> &vert_vel, std::vector< ChVector< int >> &triangles)
 A –> B Get the collision mesh in a pointer-less way, where vertex are given in a vector of xyz points, and triangles are given as indexes to the three vertexes in that vector. More...
 
void InputSimpleForces (const std::vector< ChVector<>> vert_forces, const std::vector< int > vert_ind)
 A <– B Set the forces to the nodes in a pointer-less way, where forces are given as a vector of xyz vectors and indexes to the referenced vertex, as obtained by OutputSimpleMesh. More...
 
void SetContactMesh (std::shared_ptr< ChContactSurfaceMesh > mmesh)
 Set the contact mesh (also resets the applied nodes)
 
std::shared_ptr
< ChContactSurfaceMesh
GetContactMesh ()
 Get the contact mesh.
 
std::vector< std::shared_ptr
< ChLoadXYZnode > > & 
GetForceList ()
 Access the list of applied forces, so you can add new ones by using push_back(), remove them, count them, etc. More...
 
virtual int LoadGet_ndof_x ()
 Gets the number of DOFs affected by this load (position part)
 
virtual int LoadGet_ndof_w ()
 Gets the number of DOFs affected by this load (speed part)
 
virtual void LoadGetStateBlock_x (ChState &mD)
 Gets all the current DOFs packed in a single vector (position part)
 
virtual void LoadGetStateBlock_w (ChStateDelta &mD)
 Gets all the current DOFs packed in a single vector (speed part)
 
virtual void LoadStateIncrement (const ChState &x, const ChStateDelta &dw, ChState &x_new) override
 Increment a packed state (ex. More...
 
virtual int LoadGet_field_ncoords ()
 Number of coordinates in the interpolated field, ex=3 for a tetrahedron finite element or a cable, = 1 for a thermal problem, etc. More...
 
virtual void ComputeQ (ChState *state_x, ChStateDelta *state_w)
 Compute Q, the generalized load. More...
 
virtual void ComputeJacobian (ChState *state_x, ChStateDelta *state_w, ChMatrix<> &mK, ChMatrix<> &mR, ChMatrix<> &mM)
 Compute jacobians. More...
 
virtual bool IsStiff ()
 Report if this is load is stiff. More...
 
virtual void CreateJacobianMatrices ()
 Create the jacobian loads if needed, and also set the ChVariables referenced by the sparse KRM block. More...
 
virtual void LoadIntLoadResidual_F (ChVectorDynamic<> &R, const double c)
 Adds the internal loads Q (pasted at global nodes offsets) into a global vector R, multiplied by a scaling factor c, as R += forces * c.
 
virtual void InjectKRMmatrices (ChSystemDescriptor &mdescriptor)
 Tell to a system descriptor that there are item(s) of type ChKblock in this object (for further passing it to a solver) Basically does nothing, but inherited classes must specialize this. More...
 
virtual void KRMmatricesLoad (double Kfactor, double Rfactor, double Mfactor)
 Adds the current stiffness K and damping R and mass M matrices in encapsulated ChKblock item(s), if any. More...
 
- Public Member Functions inherited from chrono::ChLoadBase
ChLoadJacobiansGetJacobians ()
 Access the jacobians (if any, i.e. if this is a stiff load)
 
virtual void Update ()
 Update: this is called at least at each time step. More...
 

Additional Inherited Members

- Protected Attributes inherited from chrono::ChLoadBase
ChLoadJacobiansjacobians
 

Member Function Documentation

virtual void chrono::fea::ChLoadContactSurfaceMesh::ComputeJacobian ( ChState state_x,
ChStateDelta state_w,
ChMatrix<> &  mK,
ChMatrix<> &  mR,
ChMatrix<> &  mM 
)
virtual

Compute jacobians.

Not needed when forces are constant, btw.

Parameters
state_xstate position to evaluate jacobians
state_wstate speed to evaluate jacobians
mKresult dQ/dx
mRresult dQ/dv
mMresult dQ/da

Implements chrono::ChLoadBase.

virtual void chrono::fea::ChLoadContactSurfaceMesh::ComputeQ ( ChState state_x,
ChStateDelta state_w 
)
virtual

Compute Q, the generalized load.

Parameters
state_xstate position to evaluate Q
state_wstate speed to evaluate Q

Implements chrono::ChLoadBase.

virtual void chrono::fea::ChLoadContactSurfaceMesh::CreateJacobianMatrices ( )
virtual

Create the jacobian loads if needed, and also set the ChVariables referenced by the sparse KRM block.

Implements chrono::ChLoadBase.

std::vector<std::shared_ptr<ChLoadXYZnode> >& chrono::fea::ChLoadContactSurfaceMesh::GetForceList ( )

Access the list of applied forces, so you can add new ones by using push_back(), remove them, count them, etc.

Note that if you add nodes, these should belong to the referenced mesh.

virtual void chrono::fea::ChLoadContactSurfaceMesh::InjectKRMmatrices ( ChSystemDescriptor mdescriptor)
virtual

Tell to a system descriptor that there are item(s) of type ChKblock in this object (for further passing it to a solver) Basically does nothing, but inherited classes must specialize this.

Reimplemented from chrono::ChLoadBase.

void chrono::fea::ChLoadContactSurfaceMesh::InputSimpleForces ( const std::vector< ChVector<>>  vert_forces,
const std::vector< int >  vert_ind 
)

A <– B Set the forces to the nodes in a pointer-less way, where forces are given as a vector of xyz vectors and indexes to the referenced vertex, as obtained by OutputSimpleMesh.

NOTE! do not insert/remove nodes from the collision mesh between the OutputSimpleMesh-InputSimpleForces pair!

Parameters
vert_forcesarray of forces (absolute xyz forces in [N])
vert_indarray of indexes to vertexes to whom you apply forces
virtual bool chrono::fea::ChLoadContactSurfaceMesh::IsStiff ( )
virtual

Report if this is load is stiff.

If so, InjectKRMmatrices will provide the jacobians of the load.

Implements chrono::ChLoadBase.

virtual void chrono::fea::ChLoadContactSurfaceMesh::KRMmatricesLoad ( double  Kfactor,
double  Rfactor,
double  Mfactor 
)
virtual

Adds the current stiffness K and damping R and mass M matrices in encapsulated ChKblock item(s), if any.

The K, R, M matrices are added with scaling values Kfactor, Rfactor, Mfactor.

Reimplemented from chrono::ChLoadBase.

virtual int chrono::fea::ChLoadContactSurfaceMesh::LoadGet_field_ncoords ( )
virtual

Number of coordinates in the interpolated field, ex=3 for a tetrahedron finite element or a cable, = 1 for a thermal problem, etc.

Implements chrono::ChLoadBase.

virtual void chrono::fea::ChLoadContactSurfaceMesh::LoadStateIncrement ( const ChState x,
const ChStateDelta dw,
ChState x_new 
)
overridevirtual

Increment a packed state (ex.

as obtained by LoadGetStateBlock_x()) by a given packed state-delta. Compute: x_new = x + dw. Ex. this is called by the BDF numerical differentiation routine that computes jacobian in the default ComputeJacobian() fallback, if not overriding ComputeJacobian() with an analytical form.

Implements chrono::ChLoadBase.

void chrono::fea::ChLoadContactSurfaceMesh::OutputSimpleMesh ( std::vector< ChVector<>> &  vert_pos,
std::vector< ChVector<>> &  vert_vel,
std::vector< ChVector< int >> &  triangles 
)

A –> B Get the collision mesh in a pointer-less way, where vertex are given in a vector of xyz points, and triangles are given as indexes to the three vertexes in that vector.

Similarly to Wavefront .OBJ meshes. Note, indexes are 0-based. These vectors can be later sent to another computing node that computes, say, CFD forces on the mesh.

Parameters
vert_posarray of vertexes (absolute xyz positions)
vert_velarray of vertexes (absolute xyz velocities, might be useful)
trianglesarray of triangles (indexes to vertexes, ccw)