Direct computation of 3-d stress intensity factors of straight and curved planar cracks with the p-version finite element method and contour integral method

Abstract

This paper presents direct computations of 3-D fracture parameters including stress intensity factors (SIFs) and T-stress for straight and curved planar cracks with the p-version finite element method (P-FEM) and contour integral method (CIM). No excessive singular element or enrichment function is required for the computation. To demonstrate the accuracy and efficiency of the proposed approaches, several benchmark numerical models of 3-D planar straight and curved cracks with single and mixed-mode fractures are considered and analyzed: a through thickness edge straight crack in a homogeneous material, a through thickness inclined straight crack, a penny-shaped crack embedded in a cube and a central ellipse shaped crack in a homogeneous cube. Numerical results are analyzed and compared with analytical solutions and those reported by the extended finite element method (XFEM) and the scaled boundary finite element method (SBFEM) in the available literature. Numerical experiments show the accuracy, robustness and effectiveness of the present method.