The symmetric Galerkin boundary element method is an approach characterized by the use of the standard continuous formulation centered on integral equations, based on the combined use of single-layer and double-layer sources, so that the integral operator will be symmetric. The difficulty in this method is the higher order singularity of the kernel functions.

In this paper, the authors examine singular integrals in the two-dimensional elastodynamic symmetric Galerkin boundary element method. The weak singular double integrals in the space domain are evaluated by splitting them into two parts: the singular part is evaluated with an analytic approach, and the regular part is computed by direct numerical evaluation. The authors also provide some numerical examples in order to demonstrate the accuracy of the proposed symmetric Galerkin boundary element method.