Variational integrators and the Newmark algorithm for conservative and dissipative mechanical systems

The purpose of this work is twofold. First, we demonstrate analytically that the classical Newmark family as well as related integration algorithms are variational in the sense of the Veselov formulation of discrete mechanics. Such variational algorithms are well known to be symplectic and momentum preserving and to often have excellent global energy behaviour. This analytical result is veri ed through numerical examples and is believed to be one of the primary reasons that this class of algorithms performs so well. Second, we develop algorithms for mechanical systems with forcing, and in particular, for dissipative systems. In this case, we develop integrators that are based on a discretization of the Lagrange dʹAlembert principle as well as on a variational formulation of dissipation. It is demonstrated that these types of structured integrators have good numerical behaviour in terms of obtaining the correct amounts by which the energy changes over the integration run. C