A Correction Method for Crack Area Overestimation in Phase-Field Fracture

TL;DR

This paper presents a novel, mesh-independent correction method for overestimated crack areas in phase-field fracture models, improving accuracy in 2D and 3D simulations by leveraging energy equipartition principles.

Contribution

A new correction framework based on energy equipartition that accurately estimates crack areas, overcoming limitations of existing mesh-dependent and skeletonization methods.

Findings

The method accurately corrects crack area overestimation in benchmark tests.

It is applicable to complex 3D geometries with curvilinear cracks.

Validation shows improved accuracy over existing correction strategies.

Abstract

Phase-field fracture models are known to overestimate the crack area, a discrepancy that compromises the accuracy of fracture predictions. This issue stems from the diffuse crack representation and numerical artifacts, such as strain localization, where the phase-field variable artificially saturates across finite elements. Existing correction strategies, including mesh-dependent factors and skeletonization algorithms, have significant limitations. Mesh-based corrections are often unreliable for unstructured meshes, while skeletonization can be complex and inaccurate for intricate crack topologies, especially in three dimensions. This paper introduces a novel and robust framework to correct this overestimation. Our approach is founded on the principle of energy equipartition, where the energy contributions from the phase-field and its gradient are equal as the length-scale parameter approaches zero. Since numerical artifacts primarily affect the phase-field term while leaving the gradient term largely unperturbed, we propose that the crack area can be accurately approximated as twice the gradient-dependent energy. This method is inherently mesh-independent and readily applicable to the entire domain, including 3D simulations. The proposed methodology is validated against benchmarks with analytical solutions and compared with established methods like skeletonization to demonstrate its accuracy. It is then applied to complex geometries with curvilinear crack paths and evaluated in a three-dimensional simulation.