Non‐linear dynamics of three‐dimensional rods: Exact energy and momentum conserving algorithms

Abstract

The long-term dynamic response of non-linear geometrically exact rods under-going finite extension, shear and bending, accompanied by large overall motions, is addressed in detail. The central objective is the design of unconditionally stable time-stepping algorithms which exactly preserve fundamental constants of the motion such as the total linear momentum, the total angular momentum and, for the Hamiltonian case, the total energy. This objective is accomplished in two steps. First, a class of algorithms is introduced which conserves linear and angular momentum. This result holds independently of the definition of the algorithmic stress resultants. Second, an algorithmic counterpart of the elastic constitutive equations is developed such that the law of conservation of total energy is exactly preserved. Conventional schemes exhibiting no numerical dissipation, symplectic algorithms in particular, are shown to lead to unstable solutions when the high frequencies are not resolved. Compared to conventional schemes there is little, if any, additional computational cost involved in the proposed class of energy-momentum methods. The excellent performance of the new algorithm in comparison to other standard schemes is demonstrated in several numerical simulations. J. C. SIMO. N. TARNOW AND M. DOBLARE class of mechanical systems incorporates finite rotations and large deformations without restrictions placed on the allowable flexibility, and furnishes the canonical model problem for nonlinear structural dynamics.Rod models of the type considered here are fairly classical. In fact, the local balance of momentum equations are essentially contained in the works of Euler, Clebch, Maxwell, Kirchhoff and others (see e.g. Reference 2). Extensions to incorporate effects such as transverse shear deformation and warping distortion have been addressed by a number of authors from different perspectives; see References 3-7 and references therein. The key point from a numerical analysis perspective concerns the specific choice of parametrization employed in the mathematical description of the kinematics of the rod. The parametrization adopted here, introduced in Reference 5, renders a form of the momentum equations which strongly resembles the classical Euler equations of rigid body dynamics and is well suited both for mathematical and numerical analysis (see References 8-10 and the work of Mielke.") The same parametrization is adopted in Reference 12; alternative approaches are discussed in References 13 and 14.The preceding non-linear rod model can be shown to define an infinite dimensional Hamiltonian system.' Energy-momentum schemes for finite dimensional Hamiltonian systems with a linear configuration space have been considered by a number of authors, see e.g. References 16-18 for the case of classical particle mechanics. However, considerable difficulties arise when the configuration manifold is non-linear, as in the case of rigid body mechanic^,'^*^^ or for the infinite dimensional case which is the situation of interest here...

A constrained generalised- $$\alpha $$ α method for coupling rigid parallel chain kinematics and elastic bodies

“…The GEMM from Kuhl and Crisfield [23] is based on the energymomentum method developed by Simo and colleagues (see Refs. [31][32][33] for examples), which preserves energy and momentum exactly within each time step. However, due to convergence problems resulting from higher frequencies in stiff structural dynamics, controllable numerical dissipation is often desirable to prevent non-physical responses in these analyses.…”

Consistent structural linearisation in flexible-body dynamics with large rigid-body motion

“…[33]. The resulting ordinary differential equations can be solved by any standard time-integration method [34][35][36][37].…”

Dual extremum principles for geometrically exact finite strain beams

“…We consider the geometrically exact finite strain beam theory introduced by Reissner (1973Reissner ( , 1981 [34,35] and Antman (1974) [1] and further extended by Simo and co-workers (1985,1986,1991,1995) [43,45,46,44]. Since these pioneering works, considerable progress has been made on the geometrically exact analysis of three-dimensional framed structures, from both theoretical and numerical points of view, see e.g.…”