Recovery strategies, a posteriori error estimation, and local error indication for second-order G/XFEM and FEM

Author(s)

Bento, Murilo H. C. ; Proenca, Sergio P. B. ; Duarte, C. Armando

• Document type: Journal article
• Source: INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING; v. 124, n. 13, p. 38-pg., 2023-03-23.

Abstract

This article presents a computationally efficient and straightforward to implement a posteriori error estimator for second-order G/XFEM and FEM approximations. The formulation is based on the recently proposed block-diagonal Zienkiewicz-Zhu (ZZ-BD) a posteriori error estimator. The focus is on linear elastic fracture mechanics (LEFM) problems but the proposed error estimator formulation is general and can be adapted to other types of problems such as those involving material interfaces. The proposed ZZ-BD error estimator is based on a new strategy to recover stress fields from second-order G/XFEM and FEM approximations. This recovery procedure involves locally weighted L-2 projections of raw stresses on an approximation space for discontinuous and singular stress fields. The basis functions for these stress approximations are defined using a low-order partition of unity together with polynomial, discontinuous, and singular recovery enrichment functions. These singular enrichments are provided by the gradient of G/XFEM enrichments for LEFM problems and are, thus, available in most G/XFEM implementations. They also lead to a more computationally efficient stress recovery procedure than the one in the original ZZ-BD error estimator. The adopted G/XFEM are optimally convergent with the error in the energy norm of 0(h(2)). Numerical experiments show that the error of the recovered stress field converges at the same rate as the discretization error. They also show that the effectivity index of the error estimator is close to the unity and that it provides good local error indicators for mesh adaptivity algorithms. Second-order FEM for elasticity problems is a special case of the G/XFEM adopted here. This is exploited to define a block-diagonal ZZ error estimator for the FEM as a special case of the proposed estimator for the G/XFEM. This brings the computational efficiency, simplicity, and accuracy of the ZZ-BD error estimator, originally developed for the G/XFEM, to second-order standard FEM. (AU)

• FAPESP's process: 19/00435-3 - Hp-adaptive formulation for the global-local generalized finite element method and applications
• Grantee: Murilo Henrique Campana Bento
• Support Opportunities: Scholarships in Brazil - Doctorate
• FAPESP's process: 20/14463-6 - New developments to the hp-adaptive global-local generalized FEM and applications to 3-D fracture mechanics problems
• Grantee: Murilo Henrique Campana Bento
• Support Opportunities: Scholarships abroad - Research Internship - Doctorate