Evaluating Fracture Energy Predictions Using Phase-Field and Gradient-Enhanced Damage Models for Elastomers

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 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.