Convolution quadrature method‐based symmetric Galerkin boundary element method for 3‐d elastodynamics

The boundary integral equations in 3‐d elastodynamics contain convolution integrals with respect to the time. They can be performed analytically or with the convolution quadrature method. The latter time‐stepping procedure's benefit is the usage of the Laplace‐transformed fundamental solution. Therefore, it is possible to apply this method also to problems where analytical time‐dependent fundamental solutions might not be known. To obtain a symmetric formulation, the second boundary integral equation has to be used which, unfortunately, requires special care in the numerical implementation since it involves hypersingular kernel functions. Therefore, a regularization for closed surfaces of the Laplace‐transformed elastodynamic kernel functions is presented which transforms the bilinear form of the hypersingular integral operator to a weakly singular one. Supplementarily, a weakly singular formulation of the Laplace‐transformed elastodynamic double layer potential is presented. This results in a time domain boundary element formulation involving at least only weakly singular integral kernels. Finally, numerical studies validate this approach with respect to different spatial and time discretizations. Further, a comparison with the wider used collocation method is presented. It is shown numerically that the presented formulation exhibits a good convergence rate and a more stable behavior compared with collocation methods. Copyright © 2008 John Wiley & Sons, Ltd.