Energy and momentum conserving algorithms for rigid body contact

Energy-momentum conserving methods are developed for rigid body dynamics with contact. Because these methods are unconditionally stable, they are not time step dependent and, hence, are well suited for incorporation into structural mechanics finite element codes. Both penalty and Lagrange multiplier methods are developed herein and are the extension of the energy-momentum conserving integration schemes for rigid bodies given by Simo and Wong [1]