Three-dimensional fracture simulation with a single-domain,direct boundary element formulation

An efficient hypersingular boundary integral equation method for three-dimensional fracture mechanics was presented in a previous paper. The details of the numerical implementation of this method are further discussed herein. In particular, an algorithm for achieving the required differentiability of the crack surface displacement function is discussed. To illustrate the utility of the method, computational results for several strongly interacting multiple-crack geometries are presented. The calculated stress intensity factors are in excellent agreement with those obtained by an approximate analytical method due to Kachanov and Laures.