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 verified 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.
|Original language||English (US)|
|Number of pages||31|
|Journal||International Journal for Numerical Methods in Engineering|
|State||Published - Dec 2000|
ASJC Scopus subject areas
- Numerical Analysis
- Applied Mathematics