Derivation of the integral equation for general 3D crack problems was examined based on the theory of body force method. In the present analysis, stress intensity factors (SIFs) along a front of arbitrary shaped 3D planar crack are obtained directly only by solving simultaneous equations expressing a boundary condition. The crack surface is discretized using number of triangular elements and the variation of the force doublet embedded in each triangle is assumed at constant. The derived boundary integral equation was transformed into a set of simultaneous equations and was solved computationally. In order to improve the accuracy of the numerically examined boundary integral, a polar transformation scheme combined with Tayler expansion of the fundamental solutions is introduced. Not only a single crack problem but also an interference among coplanar cracks can be calculated using the unique program developed in this research. It was verified that as the number of triangular elements increases, the evaluated SIF converges to the reference solution.
展开▼
