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

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 \(\Gamma \) 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.