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.