A novel mesh regularization approach based on finite element distortion potentials: Application to material expansion processes with extreme volume change

TL;DR

This work proposes a novel mesh regularization approach allowing to restore a non-distorted high-quality mesh in an adaptive manner without the need for expensive re-meshing procedures.

Abstract

The accuracy of finite element solutions is closely tied to the mesh quality. In particular, geometrically nonlinear problems involving large and strongly localized deformations often result in prohibitively large element distortions. In this work, we propose a novel mesh regularization approach allowing to restore a non-distorted high-quality mesh in an adaptive manner without the need for expensive re-meshing procedures. The core idea of this approach lies in the definition of a finite element distortion potential considering contributions from different distortion modes such as skewness and aspect ratio of the elements. The regularized mesh is found by minimization of this potential. Moreover, based on the concept of spatial localization functions, the method allows to specify tailored requirements on mesh resolution and quality for regions with strongly localized mechanical deformation and mesh distortion. In addition, while existing mesh regularization schemes often keep the boundary nodes of the discretization fixed, we propose a mesh-sliding algorithm based on variationally consistent mortar methods allowing for an unrestricted tangential motion of nodes along the problem boundary. Especially for problems involving significant surface deformation (e.g., frictional contact), this approach allows for an improved mesh relaxation as compared to schemes with fixed boundary nodes. To transfer data such as tensor-valued history variables of the material model from the old (distorted) to the new (regularized) mesh, a structure-preserving invariant interpolation scheme for second-order tensors is employed, which has been proposed in our previous work and is designed to preserve important mechanical properties of tensor-valued data such as objectivity and positive definiteness... {continued see pdf}

Figures (26)

• Figure 1: Notation and kinematics to depict the interaction between two deformable bodies.
• Figure 2: Geometry with sharp edge
• Figure 3: Illustration of edge vectors defined to formulate included angle and edge constraints for a hexahedral element.
• Figure 4: Illustration of constraints: (a) initial geometry (b) resulting geometry when elemental constraints with $l^{1}_{\text{r}}=l^{2}_{\text{r}}=l^{3}_{\text{r}}$ and $\theta^{mn}_r=\pi/2$ are applied on (a), and (b) when $l^{1}_{\text{r}}\neq l^{2}_{\text{r}}\neq l^{3}_{\text{r}}$ with $\theta^{mn}_r=\pi/2$ are applied on (a).
• Figure 5: Illustration of mesh sliding interface surfaces with auxiliary boundary $\Gamma^{(2)}_m$.

Theorems & Definitions (8)

• Remark 2.3.1
• Remark 2.4.1
• Remark 3.1.1
• Remark 3.1.2
• Remark 3.1.3
• Remark 3.1.4
• Remark 3.1.5
• Remark 3.1.6