Quasi-static crack propagation with a Griffith criterion using a variational discrete element method

Frédéric Marazzato,Alexandre Ern,Laurent Monasse

DOI: https://doi.org/10.1007/s00466-021-02102-5

2021-10-21

Abstract:A variational discrete element method is applied to simulate quasi-static crack propagation. Cracks are considered to propagate between the mesh cells through the mesh facets. The elastic behaviour is parametrized by the continuous mechanical parameters (Young modulus and Poisson ratio). A discrete energetic cracking criterion coupled to a discrete kinking criterion guide the cracking process. Two-dimensional numerical examples are presented to illustrate the robustness and versatility of the method.

Numerical Analysis,Computational Engineering, Finance, and Science

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop a method based on the variational discrete element method (VDEM) combined with the Griffith criterion to simulate quasi - static crack propagation. Specifically, the paper aims to:

1. **Propose a new numerical method**: Use the variational discrete element method to simulate the crack propagation process in materials. This method can simulate crack propagation through mesh faces between mesh elements.

2. **Introduce the energy release rate and kink criteria**: In order to guide the crack propagation process, the paper introduces the discrete energy release rate approximation and the discrete kink criterion. These criteria are used to determine the propagation path and direction of the crack.

3. **Verify the effectiveness and robustness of the method**: Through two - dimensional numerical experiments, verify the effectiveness and robustness of the proposed method in simulating quasi - static crack propagation.

### Background and Motivation

In the fields of materials science and engineering, crack propagation is an important research topic. The traditional discrete element method (DEM) performs well in simulating granular materials, soil and rock mechanics, but has limitations when dealing with crack propagation in continuous media. Therefore, a new method needs to be developed to more accurately simulate the crack propagation process.

### Method Overview

The method proposed in the paper is based on the following points:

- **Variational discrete element method (VDEM)**: This method can handle polygonal meshes and can derive the discrete form of elastic problems through variational principles.
- **Griffith criterion**: This criterion is based on the concept of energy release rate and is used to judge whether a crack will propagate at a certain location.
- **Discrete approximation of energy release rate**: Calculate the approximate value of the energy release rate at each vertex to guide crack propagation.
- **Kink criterion**: Used to determine the direction of crack propagation, select the internal mesh face that maximizes the local elastic energy density as the next fracture face.

### Numerical Experiments

The paper verifies the effectiveness and robustness of the proposed method through a series of two - dimensional numerical experiments. The experimental results show that this method can accurately simulate the crack propagation process under different conditions.

### Conclusion

The paper successfully develops a quasi - static crack propagation simulation method based on the variational discrete element method and the Griffith criterion. This method can not only handle polygonal meshes, but also effectively guide the crack propagation process through the energy release rate and kink criteria. The numerical experimental results verify the effectiveness and robustness of this method, providing strong support for future research.

### Formula Display

The formulas involved in the paper are as follows:

- **Linearized strain tensor**:
\[
\varepsilon(u):=\frac{1}{2}(\nabla u + (\nabla u)^T)\in\mathbb{R}^{2\times2}
\]

- **Stress tensor**:
\[
\sigma(u):=C:\varepsilon(u)
\]
where \(C\) is a fourth - order stiffness tensor.

- **Energy release rate**:
\[
G(y):=\frac{1 - \nu^2}{E}(K_1(y)^2 + K_2(y)^2)
\]
where \(K_1(y)\) and \(K_2(y)\) are stress intensity factors (SIF), defined as:
\[
K_1(y):=\lim_{y'\to y}\frac{\sigma_{nn}(y')}{\sqrt{2\pi d(y,y')}}
\]
\[
K_2(y):=\lim_{y'\to y}\frac{\sigma_{n\tau}(y')}{\sqrt{2\pi d(y,y')}}
\]

- **Discrete energy release rate**:
\[
G_h(v):=\max_{F\in F^{\Gamma}_v(t_k,m)}\max_{F'\in F^i_v(t_k,m)}\pi n_F