Marcel Moldenhauer

Fleurianne Bertrand, Marcel Moldenhauer, Gerhard Starke

A stress equilibration procedure for hyperelastic material models is proposed andanalyzed in this paper. Based on the displacement-pressure approximation computed with a stable finite element pair, it constructs, in a vertex-patch-wise manner, an $H(div)$-conforming approximation to the first Piola-Kirchhoff stress. This is done in such a way that its associated Cauchy stress is weakly symmetric in the sense that its anti-symmetric part is zero tested against continuous piecewise linear functions. Our main result is the identification of the subspace of test functions perpendicular to the range of the local equilibration system on each patch which turn out to be rigid body modes associated with the current configuration. Momentum balance properties are investigated analytically and numerically and the resulting stress reconstruction is shown to provide improved results for surface traction forces by computational experiments.

Fleurianne Bertrand, Marcel Moldenhauer, Gerhard Starke

The nonconforming triangular piecewise quadratic finite element space by Fortin and Soulie can be used for the displacement approximation and its combination with discontinuous piecewise linear pressure elements is known to constitute a stable combination for incompressible linear elasticity computations. In this contribution, we extend the stress reconstruction procedure and resulting guaranteed a posteriori error estimator developed by Ainsworth, Allendes, Barrenechea and Rankin \cite{AinAllBarRan:12} and by Kim \cite{Kim:12a} to linear elasticity. In order to get a guaranteed reliability bound with respect to the energy norm involving only known constants, two modifications are carried out: (i) the stress reconstruction in next-to-lowest order Raviart-Thomas spaces is modified in such way that its anti-symmetric part vanishes in average on each element; (ii) the auxiliary conforming approximation is constructed under the constraint that its divergence coincides with the one for the nonconforming approximation. An important aspect of our construction is that all results hold uniformly in the incompressible limit. Local efficiency is also shown and the effectiveness is illustrated by adaptive computations involving different Lamé parameters including the incompressible limit case.