Stress intensity factors for three-dimensional surface cracks using enriched finite elements

The analysis of three-dimensional crack problems using enriched crack tip elements is examined in this paper. It is demonstrated that the enriched nite element approach is a very e ective technique for obtaining stress intensity factors for general three-dimensional crack problems. The in uence of compatibility, integration, element shape function order, and mesh re nement on solution convergence is investigated to ascertain the accuracy of the numerical results. It is shown that integration order has the greatest impact on solution accuracy. Sample results are presented for semi-circular surface cracks and compared with previously obtained solutions available in the literature. Good agreement is obtained between the di erent numerical solutions, except in the small zone near the free surface where previously published results have often neglected the change in the stress singularity at the free surface. The enriched crack tip element appears to be particularly e ective in this region, since boundary conditions can be easily imposed on the stress intensity factors to accurately represent the correct free surface condition.