On Numerical Integration of Discontinuous Approximations in Partition of Unity Finite Elements

Abstract

This contribution presents two advances in the formulation of discontinuous approximations in finite elements. The first method relies on Schwarz-Christoffel mapping for integration on arbitrary polygonal domains [1]. When an element is split into two subdomains by a piecewise continuous discontinuity, each of these polygonal domains is mapped onto a unit disk on which cubature rules are utilized. This suppresses the need for the usual two-level isoparametric mapping. The second method relies on strain smoothing applied to discontinuous finite element approximations. By writing the strain field as a non-local weighted average of the compatible strain field, integration on the surface of the finite elements is transformed into boundary integration, so that the usual subdivision into integration cells is not required, an isoparametric mapping is not needed and the derivatives of the shape (enrichment) functions do not need to be computed. Results in fracture mechanics and composite materials are presented and both methods are compared in terms of accuracy and simplicity. The interested reader is referred to [1,6,13] for more details and should contact the authors to receive a version of the MATLAB codes used to obtain the results herein.