Virtual elements on agglomerated finite elements to increase the critical time step in elastodynamic simulations

N. Sukumar,Michael R. Tupek

DOI: https://doi.org/10.48550/arXiv.2110.00514

2022-05-31

Abstract:In this paper, we use the first-order virtual element method (VEM) to investigate the effect of shape quality of polyhedra in the estimation of the critical time step for explicit three-dimensional elastodynamic finite element (FE) simulations. Low-quality finite elements are common when meshing realistic complex components, and while tetrahedral meshing technology is generally robust, meshing algorithms cannot guarantee high-quality meshes for arbitrary geometries or for non-water-tight computer-aided design models. For reliable simulations on such meshes, we consider FE meshes with tetrahedral and prismatic elements that have badly-shaped elements$-$tetrahedra with dihedral angles close to $0^\circ$ and $180^\circ$, and slender prisms with triangular faces that have short edges$-$and agglomerate such `bad' elements with neighboring elements to form a larger polyhedral virtual element. On each element, the element-eigenvalue inequality is used to estimate the critical time step. For a suite of illustrative finite element meshes with $\epsilon$ being a mesh-coordinate parameter that leads to poor mesh quality, we show that adopting VEM on the agglomerated polyhedra yield critical time steps that are insensitive as $\epsilon \rightarrow 0$. The significant reduction in solution time on meshes with agglomerated virtual elements vis-$à$-vis tetrahedral meshes is demonstrated through explicit dynamics simulations on a tapered beam.

Numerical Analysis

What problem does this paper attempt to address?

This paper aims to solve the problem of overly small time steps due to poor mesh quality in explicit three - dimensional elastodynamic finite - element simulations. Specifically, when using low - quality tetrahedral meshes (for example, slender tetrahedra with dihedral angles close to 0° or 180°) for simulations, these poorly - shaped elements will significantly reduce the critical time step, making the simulation impractical. To solve this problem, the author adopts the first - order Virtual Element Method (VEM) to study the influence of polyhedral shape quality on estimating the critical time step by aggregating these low - quality elements to form larger polyhedral virtual elements.

The key points of the paper include:
1. **Problem background**: When dealing with actual components of complex geometries, high - quality mesh generation is a challenge. Although tetrahedral mesh techniques are generally robust, existing mesh algorithms cannot guarantee the generation of high - quality meshes for arbitrary geometries or non - closed Computer - Aided Design (CAD) models.
2. **Solution**: By using VEM, low - quality tetrahedral elements can be aggregated into larger polyhedral virtual elements. This method can not only increase the critical time step but also significantly reduce the solution time while maintaining computational accuracy.
3. **Theoretical basis**: VEM is based on a variational framework and constructs a discretized bilinear form by defining virtual basis functions and projection operators. These basis functions do not need to be explicitly calculated when implementing the method, thus providing the ability to handle arbitrary - shaped polygonal and polyhedral meshes.
4. **Experimental verification**: Through a series of numerical experiments on finite - element meshes with different mesh parameters, it is proved that after using VEM to aggregate polyhedral elements, the sensitivity of the critical time step to mesh quality is greatly reduced. In particular, in an explicit dynamic simulation of a three - dimensional tapered beam, it is shown that using polyhedral virtual elements can significantly accelerate the computational speed compared to traditional tetrahedral elements.

In conclusion, this paper proposes an effective method. By aggregating low - quality elements to form polyhedral virtual elements, it solves the problem of an overly small critical time step caused by mesh quality problems, thereby improving the efficiency and feasibility of simulations of complex geometries.