Numerical integration of non‐linear elastic multi‐body systems

Abstract

SUMMARYThis paper is concerned with the modelling of nonlinear elastic multi-body systems discretized using the finite element method. The formulation uses Cartesian co-ordinates 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, non-linear, 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 proposed in which the equations of motion are discretized so that they imply conservation of the total energy 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 non-linear elastic multi-body systems. Examples of the procedure are presented.

Performance of time‐stepping schemes for discrete models in fracture dynamic analysis

“…This scheme is close to being one of the best time integration schemes for this kind of discrete model. However, in the presence of geometric non-linearities (global behaviour or elastic beam behaviour with large displacement), the scheme can sometimes lead to numerical problems like generation of energy rather than its dissipation [29]. The last proposed scheme has been developed in order to avoid this kind of non-physical behaviour.…”

Simulation of Flexible Multibody Systems Using Linear Graph Theory

“…Hence, the preprocessor generates look-up tables that can be exploited when needed during the dynamic simulation of FMS. To enforce the constraints at the position and velocity levels, an energy algorithm proposed by (Bauchau et al, 1995) has also been implemented.…”

A new family of generalized‐α time integration algorithms without overshoot for structural dynamics

“…Hulbert and Hughes [15] proved that HHT-method suffers from the overshooting of velocity. Bauchau et al [16] highlighted the overshooting consequences in the responses of HHTmethod when it is applied to a stiff dynamical system. Later, it was also proved that CHmethod and even the whole family of G-algorithm exhibit the overshoot property [14].…”