“…For a quantitative approach to this problem, it is evident that it is crucial to find an efficient method for calculating the applied mode II SIFs at subsurface crack tips, because this is the basic parameter for assessing both propagation threshold and growth rate. In the literature, an analytical solution is given only for a crack in an infinite body, known as the Griffith crack; other few solutions 19,20,26 have been obtained by numerical methods, such as body forces, boundary element (BEM) and FEM. The main drawback of these methods is that the model construction is usually complex and time consuming; therefore, in a crack growth process, only a few load cycles can be simulated in a reasonable time, as each step requires a re‐construction of the model.…”