On penalization in variational phase-field models of brittle fracture

Abstract

Irreversible evolution is one of the central concepts as well as implementation challenges of both the variational approach to fracture by Francfort and Marigo (1998) and its regularized counterpart by Bourdin, Francfort and Marigo (2000, 2007 and 2008), which is commonly referred to as a phase-field model of brittle fracture. Irreversibility of the crack phase-field imposed to prevent fracture healing leads to a constrained minimization problem, whose optimality condition is given by a variational inequality. In our study, the irreversibility is handled via penalization. Provided the penalty constant is well-tuned, the penalized formulation is a good approximation to the original one, with the advantage that the induced equality-based weak problem enables a much simpler algorithmic treatment. We propose an analytical procedure for deriving the optimal penalty constant, more precisely, its lower bound, which guarantees a sufficiently accurate enforcement of the crack phase-field irreversibility. Our main tool is the notion of the optimal phase-field profile, as well as the Γ-convergence result. It is shown that the explicit lower bound is a function of two formulation parameters (the fracture toughness and the regularization length scale) but is independent of the problem setup (geometry, boundary conditions etc.), the formulation ingredients (degradation function, tension–compression split etc.), as well as the discretization (mesh size). The optimally-penalized formulation is tested for three benchmark problems, including one with available analytical solution. We also compare our results with those obtained by using the alternative irreversibility technique based on the notion of the history field by Miehe et al. (2010).