Abstract:In this work, we present a first-order stabilization-free virtual element method (SFVEM) for three-dimensional hyperelastic problems. Different from the conventional virtual element method, which necessitates additional stabilization terms in the bilinear formulation, the method developed in this work operates without the need for any stabilization. Consequently, it proves highly suitable for the computation of nonlinear problems. The stabilization-free virtual element method has been applied in two-dimensional hyperelasticity and three-dimensional elasticity problems. In this work, the format will be applied to three-dimensional hyperelasticity problems for the first time. Similar to the techniques used in the two-dimensional stabilization-free virtual element method, the new virtual element space is modified to allow the computation of the higher-order projection of the gradient. This paper reviews the calculation process of the traditional projection operator; and describes in detail how to calculate the high-order projection operator for three-dimensional problems. Based on this high-order projection operator, this paper extends the method to more complex three-dimensional nonlinear problems. Some benchmark problems illustrate the capability of the stabilization-free VEM for three-dimensional hyperelastic problems.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to apply the Stabilization - Free Virtual Element Method (SFVEM) in three - dimensional hyperelastic problems. Specifically, the traditional virtual element method needs to introduce additional stabilization terms when dealing with nonlinear problems to ensure that the element stiffness matrix has the correct rank. This increases the complexity and computational cost of the method. In this paper, a virtual element method without stabilization is proposed. By increasing the order of the strain polynomial to accurately describe the energy within the element, the use of stabilization terms is avoided. This method not only improves the computational efficiency but also enhances the adaptability of the method to complex - shaped elements and high - order elements.
The main contributions of the paper are as follows:
1. **Stabilization - free term**: For the first time, a first - order stabilization - free virtual element method applicable to three - dimensional hyperelastic problems is proposed.
2. **High - order projection operator**: The high - order \(L^2\) projection operator is used to construct the stabilization - free virtual element space, enabling the method to handle higher - order nonlinear problems.
3. **Numerical verification**: The effectiveness and accuracy of this method are verified through several benchmark problems, showing its advantages in solving complex nonlinear problems.
### Overview of the paper structure
1. **Introduction**: The basic concepts of the virtual element method (VEM) and its application background in different fields are introduced.
2. **Governing equations of the finite elasticity problem**: The basic governing equations of the finite elasticity problem are reviewed, including the displacement field, deformation gradient, right Cauchy - Green deformation tensor, Green - Lagrange strain tensor, etc.
3. **Three - dimensional stabilization - free virtual element space**: The basic idea and mathematical framework of the three - dimensional stabilization - free virtual element method are described in detail, including the definitions and calculation methods of the local lifting space, \(H^1\) projection operator and high - order \(L^2\) projection operator.
4. **Projection operator**: The calculation methods of the \(H^1\) projection operator and the \(L^2\) projection operator are derived in detail, especially how to perform integration on polyhedral elements.
5. **Application of the stabilization - free virtual element method in three - dimensional hyperelastic problems**: Based on the above theory, the specific implementation of the stabilization - free virtual element method is constructed and its application in nonlinear problems is discussed.
6. **Numerical examples**: The effectiveness and accuracy of the method are verified through several benchmark problems (such as Cook's membrane problem, punch problem, cylinder torsion problem).
7. **Conclusion**: The main findings of the research and future work directions are summarized.
### Mathematical formulas
- **Deformation gradient**:
\[
F=\text{Grad}\phi = I+\nabla u
\]
- **Right Cauchy - Green deformation tensor**:
\[
C = F^T\cdot F
\]
- **Green - Lagrange strain tensor**:
\[
E=\frac{1}{2}(C - I)
\]
- **Equilibrium equation**:
\[
-\text{Div}P = f
\]
- **First Piola - Kirchhoff stress tensor**:
\[
S = F^{-1}\cdot P
\]
- **Strain energy function** (Neo - Hookean model):
\[
\Psi=\frac{\mu}{2}(I_C - 3)-\mu\ln J+\frac{\lambda}{2}(\ln J)^2
\]
- **Second Piola - Kirchhoff stress tensor**:
\[
S=\frac{\partial\Psi}{\partial C}=\mu(I - C^{-1})+\lambda(\ln J)C^{-1}
\]
### Conclusion
This paper successfully proposes a stabilization - free virtual element method applicable to three - dimensional hyperelastic problems. Through the introduction of the high - order \(L^2\) projection operator, the problem of requiring stabilization terms in the traditional virtual element method is effectively solved. The numerical experimental results show that this method is effective in dealing with complex nonlinear problems.