3D elastodynamic crack analysis by a non-hypersingular BIEM

A boundary integral equation method (BIEM) is presented for 3D elastodynamic crack analysis. The method is based on a non-hypersingular BIE formulation, where the unknown quantities are the crack opening dispacements and their derivatives. The numerical scheme applied here uses a constant shape function for elements away from the crack front, and a square-root crack-tip shape function for elements near the crack front to describe the proper behavior of the unknown quantities at the crack front. A collocation method is applied to convert the non-hypersingular BIEs to a system of linear algebraic equations which are solved numerically. For several geometrical configurations, numerical results are presented for both the elastodynamic stress intensity factors and the scattering cross section. They are in good agreement with those obtained by other authors.