The Earthquake Engineering Online Archive

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

Armero, Francisco; Romero, Ignacio

UCB/SEMM-1999/05, Dept. of Civil and Environmental Engineering, University of California, Berkeley, 1999-02, 66 pages (500/C23/1999/05)

This report presents the development of a class of time-stepping algorithms for nonlinear elastodynamics that exhibits a controllable energy dissipation in the high-frequency range, thus allowing the elimination of the modeling/discretization errors that are known to accumulate in this range of frequencies. To motivate and illustrate better the developments in this general case, the formulation and analysis of these methods are first presented for two simple model problems. Namely, a nonlinear elastic spring/mass system and a simplified model of thin elastic beams are considered. As discussed in detail, the conservation by the numerical algorithm of the momenta and corresponding relative equilibria of these characteristic Hamiltonian systems with symmetry is of main importance. These conservation properties lead for a fixed and finite time step to a correct qualitative picture of the phase space where the discrete dynamics take place, even in the presence of the desired and controlled numerical dissipation of the energy. This situation is contrasted with traditional "dissipative" numerical schemes, which are shown through rigorous analyses to not only lose their dissipative character in the general nonlinear range, but also the aforementioned conservation properties, thus leading to a qualitatively distorted approximation of the phase dynamics. The key for a successful algorithm in this context is the incorporation of the numerical dissipation in the internal modes of the motion while not affecting the group motions of the system. The algorithms presented in this work accomplish these goals. The focus in this first part is on low-order methods. Representative numerical simulations, ranging from applications in nonlinear structural dynamics to nonlinear continuum three-dimensional elastodynamics, are presented in the context of the finite element method to illustrate these ideas and results.

Available online: http://nisee.berkeley.edu/documents/SEMM/SEMM-99-05.pdf (6 MB)