In this paper, a numerical procedure, incorporated with the finite element method, is developed for calculation of the mixed-mode stress intensity factors along a 3D curved crack with non-planar surfaces. The approach is formulated by modifying the concep