Effects of the smoothness of partitions of unity on the quality of representation of singular enrichments for GFEM/XFEM stress approximations around brittle cracks

Abstract

The convergence rates of the conventional generalized/extended finite element method (GFEM/XFEM) in crack modeling are similar to the convergence rates of the finite element method (FEM) (Laborde et al., 2005 and Béchet et al., 2005) unless the crack tip enrichment functions are applied in a subdomain with fixed dimension, independently of the mesh parameter h and also demanding some special care for blending along the transition zone (Chahine et al., 2006 and Taracón et al., 2009). Thus, to improve convergence rates, more degrees of freedom (DOF) are generated due to the larger quantity of enriched nodes. This work seeks to identify and understand the advantages of better capturing the information provided by singular enrichments over mesh-based smooth partitions of unity (PoU). Such PoU with higher regularity can be built through the so-called Ck-GFEM framework, following Duarte et al. (2006), based on a moving least square of degree zero and considering mesh-based smooth weighting functions associated with arbitrary polygonal clouds. The purpose herein is to investigate some possible advantages of mesh-based smooth PoU for modeling discontinuities and singularities, in two-dimensional problems of linear elastic fracture mechanics, in such a fashion that the discretization error associated to stress discontinuities inherent in standard C0-continuous GFEM/XFEM approximations is eliminated. The procedure shares features similar to the standard FEM regarding domain partition and numerical integration but, as neither the PoU nor the enrichment functions are defined in natural domains, the integrations are performed using only global coordinates. The approximation capabilities of Ck-GFEM discretizations with different patterns of singular enrichment distribution are investigated by analyzing the convergence rates of the h and p versions, considering global measures in terms of strain energy and L2-norm of displacements. The effects on stability are also verified by analyzing the evolution of the condition number. The effect of smoothness on conditioning is investigated and the eigenvalues distributions are used to identify the several aspects involved: the smoothness of the PoU, the different types of enrichment functions and the pattern of enrichment. The performance of the smooth approximations is compared to the C0 counterparts built using conventional C0 FEM-based PoU. It is shown that smoothness in the presence of extrinsically applied singular enrichments is important, as the information provided by such enrichment is better captured. In addition, there are no stress jumps around singularities, reducing error propagation beyond the neighborhood of the singularity.