A posteriori error estimation for extended finite elements by an extended global recovery

This contribution presents an extended global derivative recovery for enriched finite element methods (FEMs), such as the extended FEM along with an associated error indicator. Owing to its simplicity, the proposed scheme is ideally suited to industrial applications. The procedure is based on global minimization of the L2 norm of the difference between the raw strain field (C-1) and the recovered (C0) strain field. The methodology engineered in this paper extends the ideas of Oden and Brauchli (Int. J. Numer. Meth. Engng 1971; 3) and Hinton and Campbell (Int. J. Numer. Meth. Engng 1974; 8) by enriching the approximation used for the construction of the recovered derivatives (strains) with the gradients of the functions employed to enrich the approximation employed for the primal unknown (displacements). We show linear elastic fracture mechanics examples, both in simple two-dimensional settings, and for a three-dimensional structure. Numerically, we show that the effectivity index of the proposed indicator converges to unity upon mesh refinement. Consequently, the approximate error converges to the exact error, indicating that the error indicator is valid. Additionally, the numerical examples suggest a novel adaptive strategy for enriched approximations in which the dimensions of the enrichment zone are first increased, before standard h- and p-adaptivities are applied; we suggest to coin this methodology e-adaptivity. Copyright © 2008 John Wiley & Sons, Ltd.