Error estimation for the polygonal finite element method for smooth and singular linear elasticity

• Published in: Computers & mathematics with applications (1987) 2021-06, Vol.92, p.109-119
• Main Authors: González-Estrada, Octavio A., Natarajan, Sundararajan, Ródenas, Juan José, Bordas, Stéphane P.A.
• Format: Article
• Language: English
• Subjects: Approximation; Engineering Sciences; Equilibrium conditions; Error analysis; Error estimation; Exact solutions; Finite element analysis; Finite element method; Laplace interpolants; Mathematical analysis; Mechanics; Polygonal finite element method; Polygons; Recovery; Singularity; Statical admissibility; Stress distribution; Structural mechanics
• ISSN: 0898-1221; 1873-7668
• DOI: 10.1016/j.camwa.2021.03.017

A recovery-based error indicator developed to evaluate the quality of polygonal finite element approximations is presented in this paper. Generalisations of the finite element method to arbitrary polygonal meshes have been increasingly investigated in the last years, as they provide flexibility in meshing and improve solution accuracy. As any numerical approximation, they have an induced error which has to be accounted for in order to validate the approximate solution. Here, we propose a recovery type error measure based on a moving least squares fitting of the finite element stress field. The quality of the recovered field is improved by imposing equilibrium conditions and, for singular problems, splitting the stress field into smooth and singular parts. We assess the performance of the error indicator using three problems with exact solution, and we also compared the results with those obtained with standard finite element meshes based on simplexes. The results indicate good values for the local and global effectivities, similar to the values obtained for standard approximations, and are always within the recommended range. •Recovery-based error indicator for polygonal FEM.•Enhanced moving least squares fitting, adapted to polygonal meshes.•Error estimates in energy norm with very good values for local and global effectivities.•Recovery procedure prepared for smooth and singular linear elasticity problems.