On the formulation of high-frequency dissipative time-stepping algorithms for nonlinear dynamics. Part I: low-order methods for two model problems and nonlinear elastodynamics

Citation
F. Armero et I. Romero, On the formulation of high-frequency dissipative time-stepping algorithms for nonlinear dynamics. Part I: low-order methods for two model problems and nonlinear elastodynamics, COMPUT METH, 190(20-21), 2001, pp. 2603-2649
Citations number
42
Categorie Soggetti
Mechanical Engineering
Journal title
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING
ISSN journal
00457825 → ACNP
Volume
190
Issue
20-21
Year of publication
2001
Pages
2603 - 2649
Database
ISI
SICI code
0045-7825(2001)190:20-21<2603:OTFOHD>2.0.ZU;2-4
Abstract
We present in this paper the development of a class of time-stepping algori thms for nonlinear elastodynamics that exhibits the controllable numerical dissipation in the high-frequency range required for the robust solution of the resulting numerically stiff systems. To motivate and illustrate better the developments in this general case, we present first the formulation an d analysis of these methods for two simple model problems. Namely, we consi der a nonlinear elastic spring/mass system and a simplified model of thin e lastic beams. As it is discussed in detail in this paper, the conservation by the numerical algorithm of the momenta and corresponding relative equili bria of these characteristic Hamiltonian systems with symmetry is of the ma in importance. These conservation properties lead for a fixed and finite ti me step to a correct qualitative picture of the phase space where the discr ete dynamics takes place, even in the presence of the desired and controlle d numerical dissipation of the energy. This situation is contrasted with tr aditional "dissipative" numerical schemes, which are shown through rigorous analyses to not only lose their dissipative character in the general nonli near range, but also the aforementioned conservation properties, thus leadi ng to a qualitatively distorted approximation of the phase dynamics. The ke y for a successful algorithm in this context is the incorporation of the nu merical dissipation in the internal modes of the motion while not affecting the group motions of the systems. The algorithms presented in this work ac complish these goals. The focus in this first part is given to first-order methods. Representative numerical simulations, ranging from applications in nonlinear structural dynamics to nonlinear continuum three-dimensional ela stodynamics, are presented in the context of the finite element method to i llustrate these ideas and results. (C) 2001 Elsevier Science B.V. All right s reserved.