Energy decaying scheme for nonlinear elastic multi-body systems

This paper is concerned with the modeling of nonlinear elastic multi-body systems discretized using the finite element method. The formulation uses Cartesian coordinates to represent the position of each elastic body with respect to a single inertial frame. The kinematic constraints among the various bodies of the system are enforced via the Lagrange multiplier technique. The resulting equations of motion are stiff, nonlinear, differential-algebraic equations. The integration of these equations presents a real challenge as most available techniques are either numerically unstable or present undesirable high frequency oscillations of a purely numerical origin. An approach is proposd in which the equations of motion are discretized so that they imply an energy decay inequality for the elastic components of the system, whereas the forces of constraint are discretized so that the work they perform vanishes exactly. The combination of these two features of the discretization guarantees the stability of the numerical integration process for nonlinear elastic multi-body systems and provides high frequency numerical dissipation. Examples of the procedure are presented and compared with other available methodologies.