Explicit time-domain approaches based on numerical Green's functions computed by finite differences - The ExGA family

The present paper describes a new family of time stepping methods to integrate dynamic equations of motion. The scalar wave equation is considered here; however, the method can be applied to time-domain analyses of other hyperbolic (e.g., elastodynamics) or parabolic (e.g., transient diffusion) problems. The algorithms presented require the knowledge of the Green's function of mechanical systems in nodal coordinates. The finite difference method is used here to compute numerically the problem Green's function; however, any other numerical method can be employed, e.g., finite elements, finite volumes, etc. The Green's matrix and its time derivative are computed explicitly through the range [0,@Dt] with either the fourth-order Runge-Kutta algorithm or the central difference scheme. In order to improve the stability of the algorithm based on central differences, an additional matrix called step response is also calculated. The new methods become more stable and accurate when a sub-stepping procedure is adopted to obtain the Green's and step response matrices and their time derivatives at the end of the time step. Three numerical examples are presented to illustrate the high precision of the present approach.

[1]  Francesco Ubertini,et al.  A methodology for the generation of low‐cost higher‐order methods for linear dynamics , 2003 .

[2]  John C. Houbolt,et al.  A Recurrence Matrix Solution for the Dynamic Response of Elastic Aircraft , 1950 .

[3]  T. C. Fung,et al.  Higher‐order accurate time‐step‐integration algorithms by post‐integration techniques , 2002 .

[4]  Kumar K. Tamma,et al.  Virtual-pulse time integral methodology, a new approach for computational dynamics: Part 1 - Theory for linear structural dynamics , 1992 .

[5]  Jintai Chung,et al.  A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method , 1993 .

[6]  William Weaver,et al.  Structural dynamics by finite elements , 1987 .

[7]  J. Z. Zhu,et al.  The finite element method , 1977 .

[8]  W. Daniel.  Analysis and implementation of a new constant acceleration subcycling algorithm , 1997 .

[9]  Kumar K. Tamma,et al.  AN OVERVIEW OF RECENT ADVANCES AND EVALUATION OF THE VIRTUAL‐PULSE (VIP) TIME‐INTEGRAL METHODOLOGY FOR GENERAL STRUCTURAL DYNAMICS PROBLEMS: COMPUTATIONAL ISSUES AND IMPLEMENTATION ASPECTS , 1996 .

[10]  Mario Paz,et al.  Structural Dynamics: Theory and Computation , 1981 .

[11]  P. Morse,et al.  Methods of theoretical physics , 1955 .

[12]  Delfim Soares,et al.  A time domain FEM approach based on implicit Green's functions for non‐linear dynamic analysis , 2005 .

[13]  W. Zhong,et al.  A Precise Time Step Integration Method , 1994 .

[14]  Delfim Soares,et al.  An efficient time/frequency domain algorithm for modal analysis of non-linear models discretized by the FEM , 2003 .

[15]  W. Thomson.  Theory of vibration with applications , 1965 .

[16]  Jintai Chung,et al.  Explicit time integration algorithms for structural dynamics with optimal numerical dissipation , 1996 .

[17]  Kumar K. Tamma,et al.  The time dimension: A theory towards the evolution, classification, characterization and design of computational algorithms for transient/ dynamic applications , 2000 .

[18]  J. Oden,et al.  High-Order Taylor-Galerkin Methods for Linear Hyperbolic Systems , 1995 .

[19]  Gregory M. Hulbert,et al.  Automatic time step control algorithm for structural dynamics , 1995 .

[20]  J. T. Oden,et al.  High-order Taylor-Galerkin and adaptive h-p methods for second-order hyperbolic systems: Application to elastodynamics , 1993 .

[21]  Kumar K. Tamma,et al.  A new unified theory underlying time dependent linear first‐order systems: a prelude to algorithms by design , 2004 .

[22]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

[23]  K. Tamma,et al.  A robust self‐starting explicit computational methodology for structural dynamic applications: Architecture and representations , 1990 .

[24]  M. Géradin,et al.  Mechanical Vibrations: Theory and Application to Structural Dynamics , 1994 .

[25]  Thomas J. R. Hughes,et al.  Improved numerical dissipation for time integration algorithms in structural dynamics , 1977 .

[26]  Webe João Mansur,et al.  Time-Domain Bem Analysis for the 2d Scalar Wave Equation: Initial Conditions Contributions to Space and Time Derivatives , 1996 .

[27]  P. Smolinski.  Subcycling integration with non-integer time steps for structural dynamics problems , 1996 .

[28]  K. Bathe.  Finite Element Procedures , 1995 .

[29]  J.A.M. Carrer,et al.  A fourth order finite difference method applied to elastodynamics: Finite element and boundary element formulations , 2004 .

[30]  T. Fung.  Unconditionally stable higher-order Newmark methods by sub-stepping procedure , 1997 .

[31]  Wanxie Zhong,et al.  On a New Time Integration Method for Solving Time Dependent Partial Differential Equations , 1996 .

[32]  Kumar K. Tamma,et al.  Time discretized operators. Part 1: towards the theoretical design of a new generation of a generalized family of unconditionally stable implicit and explicit representations of arbitrary order for computational dynamics , 2003 .

[33]  Jintai Chung,et al.  A new family of explicit time integration methods for linear and non‐linear structural dynamics , 1994 .

[34]  R. Clough,et al.  Dynamics Of Structures , 1975 .

[35]  T. C. Fung,et al.  A PRECISE TIME-STEP INTEGRATION METHOD BY STEP-RESPONSE AND IMPULSIVE-RESPONSE MATRICES FOR DYNAMIC PROBLEMS , 1997 .

[36]  Robert L. Taylor,et al.  Higher derivative explicit one step methods for non‐linear dynamic problems. Part I: Design and theory , 1990 .

[37]  G. Hahn.  A modified Euler method for dynamic analyses , 1991 .