A study on time schemes for DRBEM analysis of elastic impact wave

Abstract The precise integration and differential quadrature methods are two new unconditionally stable numerical schemes to approximate time derivative with more than the second order accuracy. Recent studies showed that compared with the Houbolt and Newmark methods, they produced more accurate solutions with large time step for the problems where response is primarily dominated by low and intermediate frequency modes. This paper aims to investigate these time schemes in the context of the dual reciprocity BEM (DRBEM) formulation of various shock-excited scalar elastic wave problems, where high modes have important affect on traction response. The Houbolt method was widely recommended in many literatures for such DRBEM dynamic formulations. However, this study found that the damped Newmark algorithm was the most efficient and accurate for impact traction analysis in conjunction with the DRBEM. The precise integration and differential quadrature methods are shown inapplicable for such shock-excited problems due to the absence of numerical damping. On the other hand, we also found that to achieve the same order of accuracy, the differential quadrature method required much less computing effort than the precise integration method due to the use of the Bartels–Stewart algorithm solving the resulting Lyapunov matrix analogization equation.

[1]  Wen Chen,et al.  New Insights in Boundary-only and Domain-type RBF Methods , 2002, ArXiv.

[2]  Demosthenes Polyzos,et al.  Some Studies on Dual Reciprocity BEM for Elastodynamic Analysis , 1996 .

[3]  M. Dokainish,et al.  A survey of direct time-integration methods in computational structural dynamics—I. Explicit methods , 1989 .

[4]  Jin Yeon Cho,et al.  From the trapezoidal rule to higher-order accurate and unconditionally stable time-integration method for structural dynamics , 1997 .

[5]  Frederic Ward Williams,et al.  A high precision direct integration scheme for structures subjected to transient dynamic loading , 1995 .

[6]  G. Dahlquist.  A special stability problem for linear multistep methods , 1963 .

[7]  Demosthenes Polyzos,et al.  Three-dimensional structural vibration analysis by the Dual Reciprocity BEM , 1998 .

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

[9]  W. Mansur,et al.  Time discontinuous linear traction approximation in time-domain BEM scalar wave propagation analysis , 1998 .

[10]  C. Bert,et al.  Differential Quadrature Method in Computational Mechanics: A Review , 1996 .

[11]  C. S. Chen,et al.  Some comments on the use of radial basis functions in the dual reciprocity method , 1998 .

[12]  M. Tanaka,et al.  Coupling dual reciprocity BEM and differential quadrature method for time-dependent diffusion problems , 2001 .

[13]  Wen Chen,et al.  DUAL RECIPROCITY BEM APPLIED TO TRANSIENT ELASTODYNAMIC PROBLEMS WITH DIFFERENTIAL QUADRATURE METHOD IN TIME , 2001 .

[14]  P. W. Partridge,et al.  The dual reciprocity boundary element method , 1991 .

[15]  C. Brebbia,et al.  A new approach to free vibration analysis using boundary elements , 1983 .

[16]  Edward L. Wilson,et al.  Numerical methods in finite element analysis , 1976 .

[17]  Richard H. Bartels,et al.  Algorithm 432 [C2]: Solution of the matrix equation AX + XB = C [F4] , 1972, Commun. ACM.

[18]  T. Hughes,et al.  Collocation, dissipation and [overshoot] for time integration schemes in structural dynamics , 1978 .

[19]  H. R. Busby,et al.  Practical inverse analysis in engineering , 1997 .

[20]  K. Park.  An Improved Stiffly Stable Method for Direct Integration of Nonlinear Structural Dynamic Equations , 1975 .

[21]  Dimitri E. Beskos,et al.  Transient dynamic elastoplastic analysis by the dual reciprocity BEM , 1993 .