Calculation of stress intensity factors for an interfacial crack between dissimilar anisotropic media, using a hybrid element method and the mutual integral

Using Beom and Atluri's complete eigen-function solutions for stresses and displacements near the tip of an interfacial crack between dissimilar anisotropic media, a hybrid crack tip finite-element is developed. This element, as well as a mutual integral method are used to determine the stress intensity factors for an interfacial crack between dissimilar anisotropic media. The hybrid element has, for its Galerkin basis functions, the eigen-function solutions for stresses and displacements embedded within it. The "mutual integral" approach is based on the application of the path-independent J integral to a linear combination of two solutions: one, the problem to be solved, and the second, an "auxiliary" solution with a known singular solution. A comparison with exact solutions is made to determine the accuracy and efficiency of both the methods in various mixed mode interfacial crack problems. The size of the hybrid element was found to have very little effect on the accuracy of the solution: an acceptable numerical solution can be obtained with a very coarse mesh by using a larger hybrid element. An equivalent domain integral method is used in the application of the "mutual" integral instead of the line integral method. It is shown that the calculated mutual integral is domain independent. Therefore, the mutual integral can be evaluated far away from the crack-tip where the finite element solution is more accurate. In addition, numerical examples are given to determine the stress intensity factors for a delamination crack in composite lap joints and at plate-stiffener interfaces.