There are three considerations when solving a 3D-elasticity problem for cracked bodies when the goal is the accurate computation of Stress Intensity Factors (SIFs):

1. Approximability of the exact solution: The best rate of convergence is obtained by p-extension with a mesh graded in geometric progression towards the crack front with a grading factor of 0.15.
2. Element mapping: Distorted elements decrease the quality of the computed stress field which is used in the extraction of SIFs by means of the Contour Integral Method (CIM).
3. Path-independence: When computing SIFs by the CIM for corner and thru-thickness cracks in 3D, the SIFs are not independent of the radius of the integration path. It is necessary to use the stress field in the close vicinity of the crack front. This point is addressed in the following text and practical guidelines are provided.