Rectangular Mixed Elements for Elasticity with Weakly Imposed symmetry Condition

Gerard Awanou

DOI: https://doi.org/10.48550/arXiv.1012.1906

2011-03-04

Abstract:We present new rectangular mixed finite elements for linear elasticity. The approach is based on a modification of the Hellinger-Reissner functional in which the symmetry of the stress field is enforced weakly through the introduction of a Lagrange multiplier. The elements are analogues of the lowest order elements described in Arnold, Falk and Winther [ Mixed finite element methods for linear elasticity with weakly imposed symmetry. Mathematics of Computation 76 (2007), pp. 1699--1723]. Piecewise constants are used to approximate the displacement and the rotation. The first order BDM elements are used to approximate each row of the stress field.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to construct a mixed finite element method with weakly imposed symmetry conditions on a rectangular grid for linear elasticity problems. Specifically, the author aims to develop a new rectangular mixed finite element method in which the symmetry of the stress field is weakly processed by introducing Lagrange multipliers. The keys of this method are as follows:

1. **Weakening symmetry conditions**: Traditional finite element methods require the stress tensor to be strictly symmetric, which is very difficult to achieve at the discrete level and needs to satisfy the Brezzi stability condition. To simplify this requirement, the author adopts a modified Hellinger - Reissner functional, in which the symmetry of the stress tensor is weakly imposed through Lagrange multipliers.

2. **Construction of low - order elements**: Based on the work of Arnold, Falk and Winther, the author proposes the lowest - order mixed finite elements suitable for rectangular grids. These elements use piecewise constants to approximate displacement and rotation, and use first - order BDM (Brezzi - Douglas - Marini) elements to approximate each row of the stress tensor.

3. **Simplifying computational complexity**: Compared with previous methods, such as the Morley element, the newly proposed elements are simpler in three - dimensional cases, and the rotational degrees of freedom can be eliminated by static condensation, thereby reducing the computational burden.

4. **Ensuring stability and convergence**: To ensure the stability and convergence of the proposed finite element method, the author verifies the Brezzi condition and proves the error estimate formula:
\[
||\sigma - \sigma_h||_{H(\text{div})}+||u - u_h||_{L^2}+||\gamma - \gamma_h||_{L^2}\leq Ch||u||_3
\]
where $\sigma$ and $u$ represent the stress tensor and displacement vector in the continuous solution respectively, and $\sigma_h$, $u_h$, $\gamma_h$ are the corresponding discrete solutions, and $h$ is the mesh size.

In summary, the main objective of this paper is to develop a mixed finite element method with weakly imposed symmetry conditions on a rectangular grid to simplify the complexity in traditional methods and ensure the stability and convergence of numerical solutions.