Premium
Higher‐order extended FEM for weak discontinuities – level set representation, quadrature and application to magneto‐mechanical problems
Author(s) -
Kästner M.,
Müller S.,
Goldmann J.,
Spieler C.,
Brummund J.,
Ulbricht V.
Publication year - 2012
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.4435
Subject(s) - classification of discontinuities , finite element method , extended finite element method , numerical integration , boundary value problem , smoothing , bilinear interpolation , computation , representation (politics) , mathematics , quadrature (astronomy) , mathematical analysis , algorithm , structural engineering , physics , engineering , statistics , optics , politics , political science , law
SUMMARY In this article, we present the application of bilinear and biquadratic extended FEM (XFEM) formulations to model weak discontinuities in magnetic and coupled magneto‐mechanical boundary value problems. For properly resolving the location of curved interfaces and the discontinuous physical behaviour, the major part of the contribution is devoted to review and develop methods for level set representation of curved interfaces and numerical integration of the weak form in higher‐order XFEM formulations. In order to reduce the complexity of the representation of curved interfaces, an element local approach that allows for an automated computation of the level set values and also improves the compatibility between the level set representation and the integration subdomains is proposed. Integration rules for polygons and strain smoothing are applied in conjunction with biquadratic elements and compared with curved integration subdomains. Eventually, a coupled magneto‐mechanical demonstration problem is modelled and solved by XFEM. For demonstration purposes, a magneto‐mechanical coupling due to magnetic stresses is considered. Errors and convergence rates are analysed for the different level set representations and numerical integration procedures as well as their dependence on the ratio of material parameters at an interface. Copyright © 2012 John Wiley & Sons, Ltd.