Stable generalized finite element method (SGFEM) for three-dimensional crack problems

Abstract

This paper proposes a stable generalized finite element method (SGFEM) for the linear 3D elasticity problem with cracked domains. Conventional material-independent branch functions serve as singular enrichments. We prove that the proposed SGFEM with the geometric enrichment scheme yields the optimal order of convergence in the energy norm, O(h), for fully 3D elasticity planar crack problems; h is the mesh parameter. To improve the conditioning of SGFEM, two stability techniques have been employed, namely, (a) a cubic polynomial has been used as the PU (partition of unity), instead of the standard FE hat-functions, to address the possible almost linear dependence between the PU functions and the enrichments, and (b) a local principal component analysis (LPCA) has been implemented to address the local bad conditioning produced by multi-fold enrichments at a node. The scaled condition number for the proposed SGFEM is shown to be O(h^{-2}) (same as that of a standard Finite Element Method), for various relative positions of crack surface and mesh. The robustness of the scaled condition number for the proposed SGFEM, with respect to the relative positions of the crack-surface and the element boundaries, has been observed numerically. The numerical experiments for both the planar and fully 3D planar crack problems are presented to show the efficiency of the proposed SGFEM.