Parametric enrichment adaptivity by the extended finite element method

An adaptive method within the extended finite element method (XFEM) framework which adapts the enrichment function locally to the physics of a problem, as opposed to polynomial or mesh refinement, is presented. The method minimizes a local residual and determines the parameters of the enrichment function. We consider an energy form and a 'strong' form of the residual as error measures to drive the algorithm. Numerical examples for boundary layers and solid mechanics problems illustrate that the procedure converges. Moreover, when only the character of the solution is known, a good approximation is obtained in the area of interest. It is also shown that the method can be used to determine the order of singularities in solutions. Copyright (c) 2007 John Wiley & Sons, Ltd.