Smoothed finite element method for stress intensity factor estimation: benefits and limitations

I.V.Rokach

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

2019-03-27

Abstract:The aim of this study was to check how efficient can be smoothed finite element method (FEM) for solution of the linear fracture mechanics problems. Accuracy of stress intensity factor (SIF) computation were investigated using three types of smoothed FEM (namely ES-FEM, NS-FEM and $\alpha$-FEM). Two types of simple benchmark problems (uniaxial tension of the specimens with central and edge cracks) have been considered. SIF values were calculated by virtual crack extension and modified crack closure integral methods on almost uniform meshes. Influence of the mesh density and type of the method used for SIF determination were investigated too. Conclusions: 1) $\alpha$-FEM has the highest accuracy among all tested methods. 2) In the case of the simplest linear triangles, ES-FEM accuracy is same as in the case of traditional FEM when second order (quadratic) elements are used. 3) NS-FEM always overestimates SIF values. Unfortunately, solutions obtained using this method or $\alpha$-FEM contain unphysical oscillations of displacements and stress near the crack tip.

Computational Engineering, Finance, and Science,Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to evaluate the efficiency of the smoothed finite element method (S - FEM) when using the finite element method (FEM) to solve linear fracture mechanics problems. Specifically, the paper focuses on the accuracy of stress intensity factor (SIF) calculations. Three types of S - FEM (i.e., ES - FEM, NS - FEM, and α - FEM) were used in the study, and SIF values were calculated on an almost uniform mesh by the virtual crack extension method and the modified crack closure method. In addition, the influence of mesh density and the type of method used to determine SIF was also studied.

### Research Background and Objectives
- **Background**: Calculating the stress intensity factor (SIF) is an important topic in applying the finite element method (FEM) to solve engineering problems. Newly emerging finite element methods are usually immediately tested for their applicability in solving fracture mechanics problems.
- **Objectives**: To evaluate the efficiency of the smoothed finite element method (S - FEM) in solving linear fracture mechanics problems, especially the accuracy of SIF calculations.

### Main Research Contents
1. **Method Introduction**:
- **ES - FEM**: Edge - smoothed finite element method, which reduces the jumps in stress and strain by averaging on the common edges of adjacent elements.
- **NS - FEM**: Node - smoothed finite element method, which reduces the total stiffness by averaging on the common nodes of the initial mesh.
- **α - FEM**: A method that combines traditional FEM and NS - FEM, and optimizes the results by adjusting the parameter α.
2. **Test Problems**:
- Use two simple benchmark problems: uniaxial tensile specimens with a central crack and an edge crack.
- Calculate SIF values and compare the accuracy of different methods.
3. **Experimental Design**:
- **Mesh Density**: The influence of different mesh densities on SIF calculations was studied.
- **SIF Calculation Method**: Use the virtual crack extension method (VCE) and the modified crack closure method (MCCI) to calculate SIF values.

### Main Conclusions
1. **α - FEM**: Among all the tested methods, α - FEM has the highest accuracy.
2. **ES - FEM**: For the simplest linear triangular element, the accuracy of ES - FEM is the same as that of the traditional FEM when using quadratic elements.
3. **NS - FEM**: It always overestimates SIF values and produces unphysical displacement and stress oscillations near the crack tip.

### Significance and Applications
- **Practical Applications**: α - FEM can obtain high - precision SIF values using the simplest linear triangular element on a very rough mesh, which is particularly useful for topology optimization problems that require fast calculations.
- **Limitations**: NS - FEM has instability near the crack tip, which limits its application in practical engineering.

Through these studies, the paper provides important references and guidance for the application of the smoothed finite element method in fracture mechanics problems.