Evaluating fracture energy predictions using phase-field and gradient-enhanced damage models for elastomers

TL;DR

This paper compares phase-field and a new gradient-enhanced damage model for predicting fracture energy in elastomers, addressing numerical challenges and proposing a novel stretch-based approach that outputs energy release rates directly.

Contribution

It introduces a novel stretch-based GED model for elastomer fracture, offering an alternative to phase-field methods with direct energy release rate computation.

Findings

Artificial viscosity stabilizes phase-field simulations but affects energy release rate measurements.

The stretch-based GED model can directly output energy release rates during crack propagation.

Numerical issues remain in the new GED approach, similar to strain-based formulations.

Abstract

Recently, the phase field method has been increasingly used for brittle fractures in soft materials like polymers, elastomers, and biological tissues. When considering finite deformations to account for the highly deformable nature of soft materials, the convergence of the phase-field method becomes challenging, especially in scenarios of unstable crack growth. To overcome these numerical difficulties, several approaches have been introduced, with artificial viscosity being among the most widely utilized. This study investigates the energy release rate due to crack propagation in hyperelastic nearly-incompressible materials and compares the phase-field method and a novel gradient-enhanced damage (GED) approach. First, we simulate unstable loading scenarios using the phase-field method, which leads to convergence problems. To address these issues, we introduce artificial viscosity to stabilize the problem and analyze its impact on the energy release rate utilizing a domain J-integral approach giving quantitative measurements during crack propagation. It is observed that the measured energy released rate during crack propagation does not comply with the imposed critical energy release rate, and shows non-monotonic behavior. In the second part of the paper, we introduce a novel stretch-based GED model as an alternative to the phase-field method for modeling crack evolution in elastomers. It is demonstrated that in this method, the energy release rate can be obtained as an output of the simulation rather than as an input which could be useful in the exploration of rate-dependent responses, as one could directly impose chain-level criteria for damage initiation. We show that while this novel approach provides reasonable results for fracture simulations, it still suffers from some numerical issues that strain-based GED formulations are known to be susceptible to.