Study of the strong prolongation equation for the construction of statically admissible stress fields: Implementation and optimization

This paper focuses on the construction of statically admissible stress fields (SA-fields) for a posteriori error estimation. In the context of verification, the recovery of such fields enables to provide strict upper bounds of the energy norm of the discretization error between the known finite element solution and the unavailable exact solution. The reconstruction is a difficult and decisive step insofar as the effectiveness of the estimator strongly depends on the quality of the SA-fields. This paper examines the strong prolongation hypothesis, which is the starting point of the Element Equilibration Technique (EET). We manage to characterize the full space of SA-fields satisfying the prolongation equation so that optimizing inside this space is possible. The computation exploits topological properties of the mesh so that implementation is easy and costs remain controlled. In this paper, we describe the new technique in details and compare it to the classical EET and to the flux-free technique for different 2D mechanical problems. The method is explained on first degree triangular elements, but we show how extensions to different elements and to 3D are straightforward.