Linear and Nonlinear Formulation of Phase Field Model with Generalized Polynomial Degradation Functions for Brittle Fractures

Tabiei, Ala

DOI: https://doi.org/10.1007/s10338-024-00501-8

2024-06-12

Acta Mechanica Solida Sinica

Abstract:The classical phase field model has wide applications for brittle materials, but nonlinearity and inelasticity are found in its stress–strain curve. The degradation function in the classical phase field model makes it a linear formulation of phase field and computationally attractive, but stiffness reduction happens even at low strain. In this paper, generalized polynomial degradation functions are investigated to solve this problem. The first derivative of degradation function at zero phase is added as an extra constraint, which renders higher-order polynomial degradation function and nonlinear formulation of phase field. Compared with other degradation functions (like algebraic fraction function, exponential function, and trigonometric function), this polynomial degradation function enables phase in [0, 1] (should still avoid the first derivative of degradation function at zero phase to be 0), so there is no convergence problem. The good and meaningful finding is that, under the same fracture strength, the proposed phase field model has a larger length scale, which means larger element size and better computational efficiency. This proposed phase field model is implemented in LS-DYNA user-defined element and user-defined material and solved by the Newton–Raphson method. A tensile test shows that the first derivative of degradation function at zero phase does impact stress–strain curve. Mode I, mode II, and mixed-mode examples show the feasibility of the proposed phase field model in simulating brittle fracture.

materials science, multidisciplinary,mechanics

What problem does this paper attempt to address?

The paper aims to address the issues of nonlinearity and irreversibility in the classical phase-field model when dealing with brittle materials. Specifically, the paper focuses on the following points:

1. **Problems of the Classical Phase-Field Model**:
- The classical phase-field model exhibits nonlinearity and irreversibility in the stress-strain curve.
- The degradation function leads to stiffness reduction at low strains.

2. **Research Objectives**:
- Propose an improved degradation function to address the above issues, enabling the model to better describe the linear elastic behavior of brittle materials.
- Study polynomial degradation functions of different orders and explore their impact on the phase-field model.
- Propose a new nonlinear phase-field model with better computational efficiency and higher length scales.

3. **Specific Methods**:
- Introduce the first derivative as an additional constraint to construct higher-order polynomial degradation functions.
- Compare the effects of different degradation functions (such as algebraic fractional functions, exponential functions, and trigonometric functions).
- Validate the feasibility of the newly proposed phase-field model in simulating brittle fracture through numerical simulations.

In summary, the main objective of this paper is to improve the phase-field model by introducing new degradation functions, enabling it to better reproduce the linear elastic behavior when dealing with brittle materials and enhancing computational efficiency.