Abstract:In this paper we analyze a virtual element method for the two dimensional elasticity spectral problem allowing small edges. Under this approach, and with the aid of the theory of compact operators, we prove convergence of the proposed VEM and error estimates, where the influence of the Lamé constants is presented. We present a series of numerical tests to assess the performance of the method where we analyze the effects of the Poisson ratio on the computation of the order of convergence, together with the effects of the stabilization term on the arising of spurious eigenvalues.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: in the case where small edges are allowed to exist, use the Virtual Element Method (VEM) to solve the two - dimensional elastic spectral problem. Specifically, the author aims to:
1. **Analyze and prove the convergence of the method**: By introducing the Virtual Element Method that allows small edges, the author hopes to prove the convergence and error estimates of this method when dealing with elastic spectral problems. This includes analyzing the influence of Lamé constants.
2. **Handle the discretization of complex geometries**: One of the significant advantages of the Virtual Element Method is that it can handle complex geometric shapes, such as domains with cracks or non - convex regions. The author hopes to further expand the application range of VEM by relaxing the restrictions on the side lengths of polygonal meshes.
3. **Study the influence of Poisson ratio and stabilization terms**: The author also focuses on the influence of Poisson ratio on the order of convergence and the influence of stabilization terms on the appearance of spurious eigenvalues. Through a series of numerical experiments, evaluate how these factors affect the accuracy of calculation results.
### Main problem description
The model problem in the paper is to solve the following system of equations:
\[
\begin{cases}
\text{div}(\sigma(w)) = -\rho \kappa w & \text{in } \Omega, \\
w = 0 & \text{on } \partial \Omega,
\end{cases}
\]
where $\Omega \subset \mathbb{R}^2$ is a bounded, open and convex region with boundary $\partial \Omega$. $\sigma(w)$ represents the stress tensor, $\rho$ is the density, $\kappa$ is the eigenvalue, and $w$ is the displacement field.
### Variational form
The variational form of this problem is: find $(\kappa, w) \in \mathbb{R} \times H_0^1(\Omega)$ such that
\[
a(w, v) = \kappa b(w, v) \quad \forall v \in H_0^1(\Omega),
\]
where
\[
a(u, v) = \int_{\Omega} \sigma(u) : \varepsilon(v) \, dx, \quad b(u, v) = \int_{\Omega} \rho u \cdot v \, dx.
\]
### Core idea of the Virtual Element Method
To solve the above problem, the author introduced the Virtual Element Method (VEM) and analyzed it under the following assumptions:
- Assumption A1: Each polygon $E \in \mathcal{T}_h$ is a star - shaped region with respect to a certain ball $B_E$, with the center of the ball being $x_E$ and the radius being $\rho_E \geq \gamma h_E$.
Through this assumption, the author defined local and global virtual spaces and introduced discrete bilinear forms $a_h(\cdot, \cdot)$ and $b_h(\cdot, \cdot)$, thereby constructing a discrete eigenvalue problem.
### Numerical experiments
The author verified the effectiveness of the method through a series of numerical experiments, including:
- Calculate the lowest four eigenfrequencies on different types of meshes.
- Study the influence of different Poisson ratios ($\nu$) and stabilization terms on the calculation results.
- Compare the effects of different stabilization terms.
- Test the performance of the method on non - convex domains (such as L - shaped domains).
These experimental results show that even when small edges are allowed, the Virtual Element Method can still effectively solve the elastic spectral problem and has good convergence and stability.