On the formulation of high-frequency dissipative time-stepping algorithms for nonlinear dynamics. Part II: second-order methods

We present in this paper the formulation of a new high-frequency dissipativ e time-stepping algorithm for nonlinear elastodynamics that is second-order accurate in time, The new scheme exhibits unconditional energy dissipation and momentum conservation (and thus the given name of energy-dissipative, momentum-conserving second-order scheme (EDMC-2)), leading also to the cons ervation of the relative equilibria of the underlying physical system. The unconditional character of these properties applies not only with respect t o the time step size but, equally important, with respect to the considered elastic potential. Moreover, the dissipation properties are fully controll ed through an algorithmic parameter, reducing to existing fully conserving schemes, if desired, The design of the new algorithm is described in detail , including a complete analysis of its dissipation/conservation properties in the fully nonlinear range of finite elasticity. To motivate the differen t constructions that lead to the dissipative properties of the final scheme , the same arguments are used first in the construction of new linear time- stepping algorithms for the system of linear elastodynamics, including firs t- and second-order schemes. The new schemes exhibit a rigorous decay of th e physical energy, with the second-order schemes formulated in a general tw o-stage formula accommodating the aforementioned control parameter in the d issipation of the energy. A complete spectral analysis of the new schemes i s presented in this linear range to evaluate their different numerical prop erties. In particular, the dissipative character of the proposed schemes in the high-frequency range is fully demonstrated. In fact, it is shown that the new second-order scheme is L-stable. Most remarkably, the extension of these ideas to the nonlinear range is accomplished avoiding the use of mult i-stage formulae, given the freedom gained in using general nonlinear relat ions, while preserving the conservation laws of the momenta and the corresp onding relative equilibria. Several representative numerical simulations ar e presented in the context of nonlinear elastodynamics to evaluate the perf ormance of the newly proposed schemes. (C) 2001 Elsevier Science B.V. All r ights reserved.