Gamma-Convergence of A Variational Cohesive Phase-Field Model
Please login to view abstract download link
Reproducing the key macroscopic features of fracture behavior under multiaxial stress states is essential for accurate modeling. Experimental evidence indicates that three intrinsic macroscopic material properties govern fracture nucleation in elastic materials: elasticity, strength, and fracture toughness (or critical energy release rate). Among these, strength remains the most often misunderstood, as it is not a single scalar quantity but rather a full surface in stress space. The flexibility in defining this strength envelope in phase-field models poses significant challenges, especially under complex loading conditions. Existing models in the literature often fail to capture both the qualitative shape and the quantitative fit of experimentally observed strength surfaces. To address this limitation, a recent work by Vicentini et al. introduces a new energy functional within a cohesive phase-field framework, specifically designed to control the shape of elastic domains. This model introduces an internal variable to describe nonlinear elasticity. Notably, the strength is decoupled from the regularization parameter, that is not interpreted as a material length scale, as often done in literature, but rather as a purely variational tool. The proposed functional allows for a rigorous variational framework, enabling the use of tools from the calculus of variations. We investigate the Gamma-convergence of the model to a sharp cohesive fracture energy in the one and two dimensional setting, using a finite finite element discrete formulation and exploiting the strong localization of the damage variable. Notably, unlike classical models where the elastic and fracture energies converge independently, this model exhibits a coupling of all energy terms. The limiting cohesive energy arises from the combined asymptotic behavior of the elastic energy (concentrated in a single element), the fracture energy, and the plastic potential, while the remaining elastic energy converges separately. We also present numerical simulations exploring the sensitivity of the model to mesh anisotropy, offering insight into both its theoretical robustness and its practical implementation.