2-D time-harmonic BEM for solids of general anisotropy with application to eigenvalue problems

We present the direct formulation of the two-dimensional boundary element method (BEM) for time-harmonic dynamic problems in solids of general anisotropy. We split the fundamental solution, obtained by Radon transform, into static singular and dynamics regular parts. We evaluate the boundary integrals for the static singular part analytically and those for the dynamic regular part numerically over the unit circle. We apply the developed BEM to eigenvalue analysis. We determine eigenvalues of full non-symmetric complex-valued matrices, depending non-linearly on the frequency, by first reducing them to the generalized linear eigenvalue problem and then applying the QZ algorithm. We test the performance of the QZ algorithm thoroughly in comparison with the FEM solution. The proposed BEM is not only a strong candidate to replace the FEM for industrial eigenvalue problems, but it is also applicable to a wider class of two-dimensional time-harmonic problems.