New Nonconforming Elements for Linear Strain Gradient Elastic Model

Hongliang Li,Pingbing Ming,Zhong-ci Shi

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

2018-09-04

Abstract:Based on a new H$^2-$Korn's inequality, we propose new nonconforming elements for the linear strain gradient elastic model. The first group of elements are H$^1-$conforming but H$^2-$nonconforming. The tensor product NTW element [Tai:2001] and the tensor product Specht triangle are two typical representatives. The second element is based on Morley's triangle with a modified elastic strain energy. We proved new interpolation error estimates for all these elements, which are key to prove uniform rates of convergence for the proposed elements. Numerical results are reported and they are consistent with the theoretical prediction.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop new non - conforming finite element methods for the linear strain - gradient elasticity model. Specifically, based on a new \(H^2 - \)Korn inequality, the author proposes new non - conforming finite elements suitable for the linear strain - gradient elasticity model. These new elements aim to simplify the complex constructions in existing methods and improve computational efficiency and accuracy.

### Problem Background

1. **Strain - Gradient Elasticity Theory**:
- Strain - gradient elasticity theory is an extension of the classical elasticity theory. It describes the strong size effect of heterogeneous materials by introducing high - order strain tensors and microscopic parameters.
- Compared with the classical Mindlin model, the linear strain - gradient elasticity model proposed by Aifantis et al. contains only one material parameter, thus simplifying the model and successfully eliminating the strain singularity in the brittle crack - tip field.

2. **Challenges in Numerical Simulation**:
- One of the main challenges in computer - simulating the strain - gradient model is the existence of a large number of material parameters in the model, which leads to complex numerical solutions.
- Although traditional \(C^1\) finite elements (such as Argyris triangles) can discretize this model, the number of degrees of freedom is huge. Especially in three - dimensional problems, high - order polynomial basis functions need to be used, and the computational cost is extremely high.

### Solution

To overcome the above difficulties, the author proposes a new non - conforming finite element method:

1. **New \(H^2 - \)Korn Inequality**:
- The author proves a new \(H^2 - \)Korn inequality and its discrete version, which provides a theoretical basis for designing robust strain - gradient finite elements.
- The innovation of the discrete \(H^2 - \)Korn inequality lies in removing the gradient - field jump term, thus avoiding the gradient - field degrees of freedom along each edge and greatly simplifying the element construction.

2. **Two Types of New Non - conforming Elements**:
- **The First Type of Elements**: \(H^1 - \)conforming but \(H^2 - \)non - conforming. These elements include NTW elements and Specht triangle elements. They have 9 degrees of freedom on each element and locally adopt quartic polynomials.
- **The Second Type of Elements**: Based on the modified Morley triangle. Although the original Morley triangle diverges in the scalar strain - gradient model, it can converge uniformly after appropriately modifying the elastic strain energy.

3. **Interpolation Error Estimation**:
- The author proves new interpolation error estimates for all the proposed elements, which is the key to proving the uniform convergence of the proposed elements.

4. **Numerical Results**:
- Numerical experimental results are provided to verify the correctness of the theoretical predictions.

### Conclusion

By introducing the new \(H^2 - \)Korn inequality and designing the corresponding non - conforming finite elements, the author successfully simplifies the numerical solution process of the strain - gradient elasticity model, improves the computational efficiency, and ensures the accuracy of the numerical solution.