Premium
Unconditionally stable algorithms for rigid body dynamics that exactly preserve energy and momentum
Author(s) -
Simo J. C.,
Wong K. K.
Publication year - 1991
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.1620310103
Subject(s) - angular momentum , rigid body dynamics , rigid body , norm (philosophy) , mathematics , classical mechanics , dynamics (music) , total angular momentum quantum number , momentum (technical analysis) , point (geometry) , mathematical analysis , physics , geometry , finance , political science , acoustics , law , economics
We show that, for rigid body dynamics, the mid‐point rule formulated in body co‐ordinates exactly conserves energy and the norm of the angular momentum for incremental force‐free motions, but fails to conserve the direction of the angular momentum vector. Further, we show that the mid‐point rule formulated in the spatial representation is, in general, physically and geometrically meaningless. An alternative algorithm is developed which exactly preserves energy, and the total spatial angular momentum in incremental force‐free motions. The implicit version of this algorithm is unconditionally stable and second order accurate. The explicit version conserves exactly angular momentum in incremental force‐free motions. Numerical simulations are presented which illustrate the excellent performance of the proposed procedure, even for incremental rotations over 65 degrees. The procedure is directly applicable to transient dynamic calculations of geometrically exact rods and shells.