Abstract:Assumed stress hybrid methods are known to improve the performance of standard displacement-based finite elements and are widely used in computational mechanics. The methods are based on the Hellinger-Reissner variational principle for the displacement and stress variables. This work analyzes two existing 4-node hybrid stress quadrilateral elements due to Pian and Sumihara [Int. J. Numer. Meth. Engng, 1984] and due to Xie and Zhou [Int. J. Numer. Meth. Engng, 2004], which behave robustly in numerical benchmark tests. For the finite elements, the isoparametric bilinear interpolation is used for the displacement approximation, while different piecewise-independent 5-parameter modes are employed for the stress approximation. We show that the two schemes are free from Poisson-locking, in the sense that the error bound in the a priori estimate is independent of the relevant Lame constant $\lambda$. We also establish the equivalence of the methods to two assumed enhanced strain schemes. Finally, we derive reliable and efficient residual-based a posteriori error estimators for the stress in $L^{2}$-norm and the displacement in $H^{1}$-norm, and verify the theoretical results by some numerical experiments.
What problem does this paper attempt to address?
The problems that this paper attempts to solve are as follows: In the plane linear elasticity problem, the traditional four - node displacement - based finite element method (such as the bilinear isoparametric element) has poor results under a coarse grid when dealing with bending problems, and the "Poisson locking" phenomenon will occur when approaching the incompressible limit (Poisson ratio ν → 0.5). Specifically, the upper bound of the error estimate depends on the Lamé constant λ, resulting in a decrease in the accuracy of the numerical solution.
To solve these problems, the paper focuses on two hybrid finite element methods based on assumed stress: the PS finite element proposed by Pian and Sumihara (Pian and Sumihara, 1984) and the ECQ4 finite element proposed by Xie and Zhou (Xie and Zhou, 2004). These methods use the Hellinger - Reissner variational principle, which allows the simultaneous solution of displacement and stress variables. The main objective of the paper is to establish the uniform convergence and posterior error estimates of these two four - node assumed - stress hybrid quadrilateral elements and verify the theoretical results.
### Main contributions
1. **No Poisson locking**: It is proved that these two elements do not depend on the Lamé constant λ in the a priori error estimate, thus avoiding the Poisson locking problem.
2. **Equivalence**: The equivalence relationship between these two hybrid finite element methods and two assumed enhanced strain methods (EAS) is established.
3. **Posterior error estimate**: Reliable and effective residual - type posterior error estimates for stress in the L2 - norm and displacement in the H1 - norm are derived, and the theoretical results are verified by numerical experiments.
### Mathematical model
Consider a bounded open set \(\Omega\subset\mathbb{R}^2\) with its boundary \(\Gamma = \Gamma_D\cup\Gamma_N\), where \(\text{meas}(\Gamma_D)>0\). The plane linear elasticity model can be expressed as:
\[
\begin{cases}
-\text{div}\sigma = f & \text{in }\Omega, \\
\sigma = C\varepsilon(u) & \text{in }\Omega, \\
u|_{\Gamma_D} = 0,\quad\sigma n|_{\Gamma_N} = g,
\end{cases}
\]
where:
- \(\sigma\in\mathbb{R}^{2\times2}_{\text{sym}}\) represents the symmetric stress tensor field,
- \(u\in\mathbb{R}^2\) represents the displacement field,
- \(\varepsilon(u)=(\nabla u + (\nabla u)^T)/2\) represents the strain,
- \(f\in\mathbb{R}^2\) represents the body force density,
- \(g\in\mathbb{R}^2\) represents the surface force,
- \(n\) represents the unit outer normal vector,
- \(C\) represents the elastic modulus tensor, defined as:
\[
C\varepsilon(u)=2\mu\varepsilon(u)+\lambda\text{tr}(\varepsilon(u))I,
\]
where \(I\) is a 2×2 identity matrix, and \(\mu\) and \(\lambda\) are Lamé parameters respectively. For the plane strain problem:
\[
\mu=\frac{E}{2(1 + \nu)},\quad\lambda=\frac{E\nu}{(1 + \nu)(1 - 2\nu)},
\]
For the plane stress problem:
\[
\mu=\frac{E}{2(1 + \nu)},\quad\lambda=\frac{E\nu}{(1 + \nu)(1-\nu)},
\]
where \(0 < \nu < 0.5\) is the Poisson ratio and \(E\) is the Young's modulus.
Through these methods, the paper...