The discrete energy-momentum method. Conserving algorithms for nonlinear elastodynamics

In the absence of external loads or in the presence of symmetries (i.e., translational and rotational invariance) the nonlinear dynamics of continuum systems preserves the total linear and the total angular momentum. Furthermore, under assumption met by all classical models, the internal dissipation in the system is non-negative. The goal of this work is the systematic design of conserving algorithms that preserve exactly the conservation laws of momentum and inherit the property of positive dissipation forany step-size. In particular, within the specific context of elastodynamics, a second order accurate algorithm is presented that exhibits exact conservation of both total (linear and angular) momentum and total energy. This scheme is shown to be amenable to a completely straightforward (Galerkin) finite element implementation and ideally suited for long-term/large-scale simulations. The excellent performance of the method relative to conventional time-integrators is conclusively demonstrated in numerical simulations exhibiting large strains coupled with a large overall rigid motion.

[1]  R. D. Vogelaere,et al.  Methods of Integration which Preserve the Contact Transformation Property of the Hamilton Equations , 1956 .

[2]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[3]  J. Marsden,et al.  Elementary classical analysis , 1974 .

[4]  V. Arnold.  Mathematical Methods of Classical Mechanics , 1974 .

[5]  Donald Greenspan,et al.  Energy and momentum conserving methods of arbitrary order for the numerical integration of equations of motion , 1975 .

[6]  Thomas J. R. Hughes,et al.  FINITE-ELEMENT METHODS FOR NONLINEAR ELASTODYNAMICS WHICH CONSERVE ENERGY. , 1978 .

[7]  J. K. Knowles.  On the dissipation associated with equilibrium shocks in finite elasticity , 1979 .

[8]  R. Ruth.  A Can0nical Integrati0n Technique , 1983, IEEE Transactions on Nuclear Science.

[9]  K. Feng.  Difference schemes for Hamiltonian formalism and symplectic geometry , 1986 .

[10]  The Hamiltonian Structure of Nonlinear Elasticity: The Convective Representation of Solids, Rods, and Plates, , 1987 .

[11]  J. M. Sanz-Serna,et al.  Runge-kutta schemes for Hamiltonian systems , 1988 .

[12]  Jerrold E. Marsden,et al.  The Hamiltonian structure of nonlinear elasticity: The material and convective representations of solids, rods, and plates , 1988 .

[13]  F. Lasagni.  Canonical Runge-Kutta methods , 1988 .

[14]  J. Marsden,et al.  Lie-Poisson Hamilton-Jacobi theory and Lie-Poisson integrators , 1988 .

[15]  C. Scovel.  Symplectic Numerical Integration of Hamiltonian Systems , 1991 .

[16]  J. C. Simo,et al.  Unconditionally stable algorithms for rigid body dynamics that exactly preserve energy and momentum , 1991 .

[17]  J. C. Simo,et al.  Nonlinear stability of the time-discrete variational problem of evolution in nonlinear heat conduction, plasticity and viscoplasticity , 1991 .

[18]  J. C. Simo,et al.  Exact energy-momentum conserving algorithms and symplectic schemes for nonlinear dynamics , 1992 .

[19]  J. C. Simo,et al.  On a stress resultant geometrically exact shell model. Part VI: Conserving algorithms for non‐linear dynamics , 1992 .