Three alternative approaches, namely the extended/generalized finite element method (XFEM/GFEM), the scaled boundary finite element method (SBFEM) and phase field methods, are surveyed and compared in the context of linear elastic fracture mechanics (LEFM). The purpose of the study is to provide a critical literature review, emphasizing on the mathematical, conceptual and implementation particularities that lead to the specific advantages and disadvantages of each method, as well as to offer numerical examples that help illustrate these features.

Summary: A variation of the extended finite element method for three‐dimensional fracture mechanics is proposed. It utilizes a novel form of enrichment and point‐wise and integral matching of displacements of the standard and enriched elements in order to achieve higher accuracy, optimal convergence rates, and improved conditioning for two‐dimensional and three‐dimensional crack problems. A bespoke benchmark problem is introduced to determine the method's accuracy in the general three‐dimensional case where it is demonstrated that the proposed approach improves the accuracy and reduces the number of iterations required for the iterative solution of the resulting system of equations by 40% for moderately refined meshes and topological enrichment. Moreover, when a fixed enrichment volume is used, the number of iterations required grows at a rate which is reduced by a factor of 2 compared with standard extended finite element method, diminishing the number of iterations by almost one order of magnitude. Copyright © 2015 John Wiley & Sons, Ltd.

A numerical scheme is proposed for the detection of multiple cracks in three dimensional (3D) structures. The scheme is based on a variant of the extended finite element method (XFEM) and a hybrid optimizer solution. The proposed XFEM variant is particularly well-suited for the simulation of 3D fracture problems, and as such serves as an efficient solution to the so-called forward problem. A set of heuristic optimization algorithms are recombined into a multiscale optimization scheme. The introduced approach proves effective in tackling the complex inverse problem involved, where identification of multiple flaws is sought on the basis of sparse measurements collected near the structural boundary. The potential of the scheme is demonstrated through a set of numerical case studies of varying complexity.