Anisotropic discretization- and model-error estimation in solid mechanics by local Neumann problems

Abstract

First, a survey of existing residuum-based error-estimators and error-indicators is given. Generally, residual error estimators (which have at least upper bound in contrast to indicators) can be locally computed from residua of equilibrium and stress-jumps at element interfaces using Dirichlet or Neumann conditions for element patches or individual elements (REM). Another equivalent method for error estimation can be derived from a posteriori computed improved boundary tractions which provide exact equilibrium of elements. They are computed from local Neumann problems and yield improved local solution by testing with higher test functions. This method is called Posterior Equilibrium Method (PEM) and yields better efficiency indices than REM for the investigated problem, even in case of locking. Also, the important feature of model-adaptivity via hierarchical error estimation can only be realized by this strategy. Global Neumann-type estimators and examples are presented for elastic and elastoplastic deformations, controlling equilibrium in tangential space.