The computation of bounds for linear-functional outputs of weak solutions to the two-dimensional elasticity equations

We present a method for the computation of upper and lower bounds for linear-functional outputs of the exact solutions to the two dimensional elasticity equations. The method can be regarded as a generalization of the well known complementary energy principle. The desired output is cast as the supremum of a quadratic-linear convex functional over an infinite dimensional domain. Using duality the computation of an upper bound for the output of interest is reduced to a feasibility problem for the complementary, or dual, problem. In order to make the problem tractable from a computational perspective two additional relaxations that preserve the bounding properties are introduced. First, the domain is triangulated and a domain decomposition strategy is used to generate a sequence of independent problems to be solved over each triangle. The Lagrange multipliers enforcing continuity are approximated using piecewise linear functions over the edges of the triangulation. Second, the solution of the adjoint problem is approximated over the triangulation using a standard Galerkin finite element approach. A lower bound for the output of interest is computed by repeating the process for the negative of the output. Reversing the sign of the computed upper bound for the negative of the output yields a lower bound for the actual output. The method can be easily generalized to three dimensions. However, a constructive proof for the existence of feasible solutions for the outputs of interest is only given in two dimensions. The computed bound gaps are found to converge optimally, that is, at the same rate as the finite element approximation. An attractive feature of the proposed approach is that it allows for a data set to be generated that can be used to certify and document the computed bounds. Using this data set and a simple algorithm, the correctness of the computed bounds can be established without recourse to the original code used to compute them. In the present paper, only computational domains whose boundary is made up of straight sided segments and polynomially varying loads are considered. Two examples are given to illustrate the proposed methodology.

[1]  Pierre Ladevèze,et al.  Error Estimate Procedure in the Finite Element Method and Applications , 1983 .

[2]  D. Arnold,et al.  Mixed Finite Elements for Elasticity in the Stress-Displacement Formulation , 2002 .

[3]  Anthony T. Patera,et al.  A hierarchical duality approach to bounds for the outputs of partial differential equations , 1998 .

[4]  D. Kelly,et al.  The self‐equilibration of residuals and complementary a posteriori error estimates in the finite element method , 1984 .

[5]  Philippe Destuynder,et al.  Explicit error bounds in a conforming finite element method , 1999, Math. Comput..

[6]  B. D. Veubeke.  Displacement and equilibrium models in the finite element method , 1965 .

[7]  M. Geradin,et al.  B.M. Fraeijs de Veubeke Memorial Volume of Selected Papers , 1980 .

[8]  J. Peraire,et al.  A posteriori finite element bounds for linear-functional outputs of elliptic partial differential equations , 1997 .

[9]  Douglas N. Arnold,et al.  Mixed finite elements for elasticity , 2002, Numerische Mathematik.

[10]  J. Tinsley Oden,et al.  Practical methods for a posteriori error estimation in engineering applications , 2003 .

[11]  Anthony T. Patera,et al.  Bounds for Linear–Functional Outputs of Coercive Partial Differential Equations : Local Indicators and Adaptive Refinement , 1998 .

[12]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[13]  R. Rannacher,et al.  A feed-back approach to error control in finite element methods: application to linear elasticity , 1997 .

[14]  Serge Prudhomme,et al.  On goal-oriented error estimation for elliptic problems: application to the control of pointwise errors , 1999 .

[15]  R. Bank,et al.  Some a posteriori error estimators for elliptic partial differential equations , 1985 .

[16]  Antonio Huerta,et al.  Computing Bounds for Linear Functionals of Exact Weak Solutions to Poisson's Equation , 2004, SIAM J. Numer. Anal..

[17]  Pierre Ladevèze,et al.  ERROR ESTIMATION AND MESH OPTIMIZATION FOR CLASSICAL FINITE ELEMENTS , 1991 .

[18]  Pedro Díez,et al.  Recovering lower bounds of the error by postprocessing implicit residual a posteriori error estimates , 2003 .

[19]  Jaime Peraire,et al.  Computing Bounds for Linear Functionals of Exact Weak Solutions to the Advection-Diffusion-Reaction Equation , 2005, SIAM J. Sci. Comput..

[20]  J. Oden,et al.  A Posteriori Error Estimation in Finite Element Analysis , 2000 .

[21]  Ivo Babuška,et al.  Guaranteed computable bounds for the exact error in the finite element solution—Part II: bounds for the energy norm of the error in two dimensions† , 2000 .