Mapped finite element methods: High‐order approximations of problems on domains with cracks and corners

Abstract

Summary Linear elasticity problems posed on cracked domains, or domains with re‐entrant corners, yield singular solutions that deteriorate the optimality of convergence of finite element methods. In this work, we propose an optimally convergent finite element method for this class of problems. The method is based on approximating a much smoother function obtained by locally reparameterizing the solution around the singularities. This reparameterized solution can be approximated using standard finite element procedures yielding optimal convergence rates for any order of interpolating polynomials, without additional degrees of freedom or special shape functions. Hence, the method provides optimally convergent solutions for the same computational complexity of standard finite element methods. Furthermore, the sparsity and the conditioning of the resulting system is preserved. The method handles body forces and crack‐face tractions, as well as multiple crack tips and re‐entrant corners. The advantages of the method are showcased for four different problems: a straight crack with loaded faces, a circular arc crack, an L‐shaped domain undergoing anti‐plane deformation, and lastly a crack along a bimaterial interface. Optimality in convergence is observed for all the examples. A proof of optimal convergence is accomplished mainly by proving the regularity of the reparameterized solution. Copyright © 2016 John Wiley & Sons, Ltd.

The efficiency of thermal fatigue in regolith generation on small airless bodies

“…The numerical model presented here is designed to efficiently model the fatigue growth of cracks in regolith over a several thousands of years, while using timesteps as small as 15 minutes. More sophisticated fracture codes, such as the mapped finite element method (Chiaramonte et al, 2017), have been developed recently.…”

A method to compute mixed‐mode stress intensity factors for nonplanar cracks in three dimensions

“…Because displacement gradients converge with only order 1/2 for the standard FEM, the final convergence rates of the stress intensity factors are very slow. Meanwhile, high-order methods, including MFEM 63 (which only has been used for 2D problems) and certain forms of XFEM, will produce faster converging gradients. This, in turn, will produce faster converging stress intensity factors, as can be seen from the continuity estimate (47).…”

An algorithm for the simulation of curvilinear plane‐strain and axisymmetric hydraulic fractures with lag using the universal meshes

“…To resolve this, we need to use a higher‐order method. Chiaramonte et al have introduced a technique, which maps the linear elasticity problem to a new solution domain, “removing” the stress singularity around the crack tip. In doing so, the method recovers the optimal convergence rates for high‐order elements.…”