A priori and a posteriori error analyses of a pseudostress-based mixed formulation for linear elasticity

Abstract

In this paper we present the a priori and a posteriori error analyses of a non-standard mixed finite element method for the linear elasticity problem with non-homogeneous Dirichlet boundary conditions. More precisely, the approach introduced here is based on a simplified interpretation of the pseudostress-displacement formulation originally proposed in Arnold and Falk (1988), which does not require symmetric tensor spaces in the finite element discretization. In addition, physical quantities such as the stress, the strain tensor of small deformations, and the rotation, are computed through a simple postprocessing in terms of the pseudostress variable. Furthermore, we also introduce a second element-by-element postprocessing formula for the stress, which yields an optimally convergent approximation of this unknown with respect to the broken ℍ(div)-norm. We apply the classical Babuška-Brezzi theory to prove that the corresponding continuous and discrete schemes are well-posed. In particular, Raviart-Thomas spaces of order k≥0 for the pseudostress and piecewise polynomials of degree ≤k for the displacement can be utilized. Moreover, we remark that in the 3D case the number of unknowns behaves approximately as 9 times the number of elements (tetrahedra) of the triangulation when k=0. This factor increases to 12.5 when one uses the classical PEERS. Next, we derive a reliable and efficient residual-based a posteriori error estimator for the mixed finite element scheme. Finally, several numerical results illustrating the performance of the method, confirming the theoretical properties of the estimator, and showing the expected behaviour of the associated adaptive algorithm, are provided.

En este artículo presentamos los análisis de error a priori y a posteriori de un método de elementos finitos mixtos no estándar para el problema de la elasticidad lineal con condiciones de contorno de Dirichlet no homogéneas. Más precisamente, el enfoque introducido aquí se basa en una interpretación simplificada de la formulación de desplazamiento de pseudoesfuerzo propuesta originalmente en Arnold y Falk (1988), que no requiere espacios tensoriales simétricos en la discretización de elementos finitos. Además, las cantidades físicas como el esfuerzo, el tensor de deformación de pequeñas deformaciones y la rotación, se calculan a través de un simple procesamiento posterior en términos de la variable de pseudoesfuerzo. Además, también introducimos una segunda fórmula de postprocesamiento elemento por elemento para el estrés, que produce una aproximación óptimamente convergente de esta incógnita con respecto a la forma rota. Aplicamos la teoría clásica de Babuška-Brezzi para demostrar que los esquemas continuos y discretos correspondientes están bien planteados. En particular, se pueden utilizar espacios de orden de Raviart-Thomas para el pseudoesfuerzo y polinomios de grado por partes para el desplazamiento. Además, observamos que en el caso 3D el número de incógnitas se comporta aproximadamente como 9 veces el número de elementos (tetraedros) de la triangulación cuando. Este factor aumenta a 12.5 cuando uno usa los PEERS clásicos. A continuación, derivamos un estimador de error a posteriori confiable y eficiente basado en residuos para el esquema de elementos finitos mixtos. Finalmente, se proporcionan varios resultados numéricos que ilustran el rendimiento del método, confirman las propiedades teóricas del estimador y muestran el comportamiento esperado del algoritmo adaptativo asociado.