Alvin Chen,N. Sukumar
Abstract:In this paper, we propose a robust low-order stabilization-free virtual element method on quadrilateral meshes for linear elasticity that is based on the stress-hybrid principle. We refer to this approach as the Stress-Hybrid Virtual Element Method (SH-VEM). In this method, the Hellinger$-$Reissner variational principle is adopted, wherein both the equilibrium equations and the strain-displacement relations are variationally enforced. We consider small-strain deformations of linear elastic solids in the compressible and near-incompressible regimes over quadrilateral (convex and nonconvex) meshes. Within an element, the displacement field is approximated as a linear combination of canonical shape functions that are $\textit{virtual}$. The stress field, similar to the stress-hybrid finite element method of Pian and Sumihara, is represented using a linear combination of symmetric tensor polynomials. A 5-parameter expansion of the stress field is used in each element, with stress transformation equations applied on distorted quadrilaterals. In the variational statement of the strain-displacement relations, the divergence theorem is invoked to express the stress coefficients in terms of the nodal displacements. This results in a formulation with solely the nodal displacements as unknowns. Numerical results are presented for several benchmark problems from linear elasticity. We show that SH-VEM is free of volumetric and shear locking, and it converges optimally in the $L^2$ norm and energy seminorm of the displacement field, and in the $L^2$ norm of the hydrostatic stress.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a robust low - order stabilization - free virtual element method (SF - VEM) on quadrilateral meshes for the stress - hybrid virtual element method (SH - VEM) of linear elastic materials under compressive and nearly incompressible conditions. Specifically, the paper aims to overcome the volumetric locking and shear locking problems that occur in traditional finite element methods when dealing with nearly incompressible materials.
### Background and Motivation
1. **Volumetric Locking and Shear Locking**:
- Traditional low - order fully integrated displacement - based finite elements are prone to volumetric locking when the Poisson ratio \( \nu \) is close to 0.5, resulting in inaccurate calculation results.
- For bending - dominated problems, the finite element method will produce spurious shear strains, leading to the shear locking phenomenon and making the elements too stiff.
2. **Development of the Virtual Element Method**:
- The introduction of the virtual element method (VEM) provides a new way to solve the locking problems. By using techniques such as the mixed variational principle, mixed formulation, B - bar method, and selective reduced integration strategy, VEM performs well in dealing with constrained problems.
- However, existing stabilization - free virtual element methods (such as SF - VEM) still have volumetric locking problems under nearly incompressible conditions.
### Main Contributions of the Paper
1. **Stress - Hybrid Virtual Element Method (SH - VEM)**:
- This paper proposes a stress - hybrid virtual element method based on the Hellinger - Reissner variational principle. In this method, the displacement field and the stress field are independent, the displacement field is continuous, and the stress field (discontinuous) is represented by a linear combination of symmetric tensor polynomials.
- By using a five - parameter extended stress field in each element and applying the stress transformation equation, SH - VEM can effectively solve linear elastic problems on quadrilateral meshes and avoid volumetric locking and shear locking.
2. **Numerical Implementation**:
- The paper details the numerical implementation process of SH - VEM, including the definition of the virtual element space, the construction of projection operators, and the calculation of element stiffness matrices and force vectors.
- Through a series of standard test problems, the effectiveness and robustness of SH - VEM are verified, especially when the Poisson ratio is close to 0.5.
### Mathematical Model and Method
1. **Hellinger - Reissner Variational Principle**:
- Starting from the Hellinger - Reissner functional, the weak form of the linear elastic problem is constructed. The equilibrium equation and the strain - displacement relationship in the variational principle are enforced through the weak form.
- By introducing stress assumptions and displacement assumptions, the relationship between stress coefficients and nodal displacements is derived, and finally a displacement formula containing only nodal displacements is obtained.
2. **Virtual Element Space and Projection Operators**:
- The virtual element space is defined, and the energy projection operator and the stress - hybrid projection operator are constructed.
- The integrity and rotation invariance of the stress field are ensured by using a rotating coordinate system and a stress transformation equation.
### Numerical Results
1. **Benchmark Test Problems**:
- The paper shows the performance of SH - VEM in multiple benchmark test problems, including thin cantilever beams, Cook diaphragms, infinite plate hole tension, hollow pressure cylinders, and flat punches.
- The results show that SH - VEM has the optimal convergence rate in the \( L^2 \) norm and the energy semi - norm, and also performs well in the \( L^2 \) norm of the hydrostatic stress.
2. **Stability Analysis**:
- Through eigenvalue analysis, the stability of SH - VEM on quadrilateral elements of different shapes, including convex and non - convex elements, is verified.
- The results show that SH - VEM does not have spurious zero - energy modes and can remain stable even when the element shape is severely distorted.
### Summary
The SH - VEM proposed in this paper successfully solves the linear elastic problem by combining the stress - hybrid technique and the virtual element method.