Simulate with complex geometries and complex physics
MassSpringDashpot
movement of a mass-spring-damper system coupled (by joints) with other %MOVE_rigid% or %MOVE_MassSpringDashpot% geometries
A mass-spring-damper requires two BND_point definitions. Those two points initially mark the beginning and end location of the spring, and form the connection bonds to other %MOVE_rigid% or %MOVE_MassSpringDashpot% geometries. In addition, optional triangles and lines can be used to compute the hydrodynamical forces on the mass-spring-damper system.
If the optional non-zero Spring Force and Damping Force are provided, then the total forces are calculated as follows:
\begin{align}{\bf F_{spring}} = {\bf F_{spring}} + {\bf K} \cdot {\bf dx}, \quad {\bf F_{damping}} = {\bf F_{damping}} + {\bf C} \cdot {\bf dv}\end{align}
\( {\bf F_{spring}}\) is the total spring force, \( {\bf F_{damping}}\) is the total damping force, \( {\bf K}\) is the springConstant, and \( {\bf C}\) is the dampingConstant.
\( {\bf dx}\) is the change in spring length which can be obtained via the function xCOG(4) or xCOG(5), see xCOG() for details.
\( {\bf dv}\) is the change in length rate which can be obtained via the function vCOG(4) or vCOG(5), see vCOG() for details.
begin_alias{}
"MSD" = " ... MOVE$MOVE_MSD$ ... "
end_alias
begin_boundary_elements{}
BND_point &MSD& x1 y1 z1 # most preferably, these points should link to other BND_point items of other MOVE-statements (%MOVE_rigid%, ... )
BND_point &MSD& x2 y2 z2 # otherwise, the spring-damper system would hang in the air
end_boundary_elements
MOVE($MOVE_MSD$) = ( %MOVE_MassSpringDashpot%, mass, precompressionForce, springConstant, dampingConstant , OPTIONAL: Spring Force, OPTIONAL: Damping Force )