Andrea Lamperti,Massimiliano Cremonesi,Umberto Perego,Alessandro Russo,Carlo Lovadina
Abstract:A recent mixed formulation of the Virtual Element Method in 2D elastostatics, based on the Hu-Washizu variational principle, is here extended to 2D elastodynamics. The independent modeling of the strain field, allowed by the mixed formulation, is exploited to derive first order quadrilateral Virtual Elements (VEs) not requiring a stabilization (namely, self-stabilized VEs), in contrast to the standard VEs, where an artificial stabilization is always required for first order quads. Lumped mass matrices are derived using a novel approach, based on an integration scheme that makes use of nodal values only, preserving the correct mass in the case of rigid-body modes. In the case of implicit time integration, it is shown how the combination of a self-stabilized stiffness matrix with a self-stabilized lumped mass matrix can produce excellent performances both in the compressible and quasi-incompressible regimes with almost negligible sensitivity to element distortion. Finally, in the case of explicit dynamics, the performances of the different types of derived VEs are analyzed in terms of their critical time-step size.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to construct self - stable quadrilateral virtual elements by the hybrid virtual element method (VEM) based on the Hu - Washizu variational principle in two - dimensional linear elastodynamics. Specifically, the paper focuses on how to achieve self - stabilization of first - order quadrilateral virtual elements without introducing artificial stabilization terms, and particularly emphasizes the self - stabilization treatment of the mass matrix. The paper proposes two types of self - stable elements: one contains internal moment degrees of freedom, and the other only uses nodal displacement degrees of freedom. In addition, the paper also explores the performance of these self - stable elements in explicit and implicit time - integration schemes, especially in terms of quasi - incompressible material behavior and element distortion sensitivity.
### Analysis of the Core Problems in the Paper
1. **Achievement of Self - Stabilization**:
- **Self - Stabilization of the Stiffness Matrix**: By selecting an appropriate strain model, the paper enables the quadrilateral virtual element to naturally avoid non - physical zero - energy modes (such as hourglass modes), thus achieving self - stabilization.
- **Self - Stabilization of the Mass Matrix**: The paper proposes a new integration scheme that uses only nodal values for integration, thereby directly generating a diagonal lumped mass matrix and ensuring that it has the correct mass distribution in rigid - body modes.
2. **Application of Low - Order Quadrilateral Virtual Elements**:
- **Advantages of Low - Order Elements**: In practical engineering applications, low - order quadrilateral elements are widely used because of their simplicity and efficiency. Through self - stabilization techniques, the paper further improves the performance of these elements, especially when dealing with highly nonlinear and dynamic problems.
- **Explicit and Implicit Time Integration**: The paper analyzes the performance of self - stabilized virtual elements in explicit and implicit time - integration schemes respectively, especially in terms of critical time step and eigenfrequency analysis.
3. **Numerical Verification**:
- **Convergence and Accuracy**: The paper verifies the convergence and accuracy of different types of self - stabilized virtual elements through multiple numerical experiments, especially in quasi - incompressible problems.
- **Eigenfrequency Analysis**: The paper also conducts eigenfrequency analysis on different types of self - stabilized elements and compares their performance in explicit time - integration schemes.
### Key Technologies and Methods
1. **Hu - Washizu Variational Principle**:
- Based on the Hu - Washizu variational principle, the paper independently models displacement, strain, and stress fields, thereby deriving the control equations in weak form.
- By introducing an appropriate strain model, self - stabilization of the stiffness matrix is achieved.
2. **Self - Stabilization of the Mass Matrix**:
- A new integration scheme is proposed, which uses only nodal values for integration, thereby directly generating a diagonal lumped mass matrix.
- The correctness of the mass matrix in rigid - body modes is ensured, and it can be directly used in explicit time - integration schemes.
3. **Numerical Experiments**:
- The performance of self - stabilized virtual elements is verified through multiple numerical experiments, especially their performance when dealing with highly nonlinear and dynamic problems.
### Conclusion
The paper successfully realizes self - stabilized quadrilateral virtual elements based on the Hu - Washizu variational principle and verifies their effectiveness and superiority in two - dimensional linear elastodynamics problems through numerical experiments. In particular, the self - stabilization mass matrix technique proposed in the paper not only improves the performance of the elements but also simplifies the implementation of explicit time - integration schemes. These achievements provide important theoretical and technical support for the further promotion of the virtual element method in engineering applications.