Abstract:It is known that the 3-node hybrid triangular element MiSP3 and 4-node hybrid quadrilateral element MiSP4 presented by Ayad, Dhatt and Batoz (Int. J. Numer. Meth. Engng 1998, 42: 1149-1179) for Reissner-Mindlin plates behave robustly in numerical benchmark tests. These two elements are based on Hellinger-Reissner variational principle, where continuous piecewise linear/isoparametric bilinear interpolations, as well as the mixed shear interpolation/projection technique of MITC family, are used for the approximations of displacements, and piecewise-independent equilibrium modes are used for the approximation of bending moments/shear stresses. We show that the MiSP3 and MiSP4 elements are uniformly stable with respect to the plate thickness and thus free from shear-locking.
What problem does this paper attempt to address?
This paper attempts to solve the shear locking problem encountered in the numerical calculation of the Reissner - Mindlin plate model. Specifically, when the plate thickness \(t\) approaches zero, the traditional low - order finite element discretization method will lead to a serious decline in the accuracy of the numerical solution and cannot accurately simulate the behavior of thin and thick plates. To solve this problem, the author analyzes two hybrid finite element methods: MiSP3 (triangular elements) and MiSP4 (quadrilateral elements), and proves the uniform stability of these elements at different thicknesses, thus avoiding the shear locking phenomenon.
### Main research content:
1. **Problem background**:
- The Reissner - Mindlin plate model is the main two - dimensional model for calculating the bending behavior of three - dimensional thick/thin plates.
- When the plate thickness \(t\) approaches zero, the standard low - order finite element method will have a shear locking problem, resulting in inaccurate numerical solutions.
2. **Existing methods**:
- Many researchers overcome the shear locking problem by modifying the variational formula or introducing shear stress interpolation techniques.
- Among them, the MITC (Mixed Interpolated Tensorial Components) series of elements is a commonly used method, but it has some limitations, especially in the application of triangular elements.
3. **Improved methods**:
- Ayad, Dhatt and Batoz proposed an improved MiSP3 and MiSP4 elements, based on the Hellinger - Reissner variational principle and combined with the shear stress interpolation technique of the MITC series.
- These elements use continuous piecewise linear/bilinear interpolation to approximate the displacement and use independent equilibrium modes to approximate the bending moment and shear stress.
4. **Main contributions**:
- The author proves that MiSP3 and MiSP4 elements have uniform stability at different plate thicknesses and can effectively avoid the shear locking phenomenon.
- The stability and convergence of these elements are verified through theoretical analysis and numerical experiments.
### Key conclusions:
- **Uniform stability**: MiSP3 and MiSP4 elements remain uniformly stable at different thicknesses and will not have shear locking.
- **Error estimate**: The author derives the error estimates of these elements and proves their effectiveness in numerical calculations.
### Mathematical expressions:
The mathematical formulas involved in the article include:
- **Helmholtz decomposition**:
\[
\tau=\text{grad} s+\text{curl} q, \quad(s, q) \in H^1_0(\Omega) \times \hat{H}^1(\Omega)
\]
where,
\[
\hat{H}^1(\Omega):=\left\{q \in H^1(\Omega): \int_\Omega q \, dx = 0\right\}
\]
- **Error estimate**:
\[
| \| M - M_h \| |_{h,1}+| \| (w - w_h, \beta - \beta_h) \| |_{h,2} \lesssim h (\| M \|_1+\| \beta \|_2+\| w \|_2+\| r \|_2+\| p \|_1+t \| p \|_2) \lesssim h \| g \|_0
\]
Through these analyses, the author successfully solves the shear locking problem encountered in the numerical calculation of the Reissner - Mindlin plate model and provides a reliable numerical tool for practical engineering applications.