this can be ordered along with an additional penalty term in the
phase-field evolution equation7, or alternatively via
a history variable \(H\) which ensures the irreversibility condition to
prevent crack healing. This variable,\(\ H\), replaces the quantity\(\psi_{e}\) in Eq.(27).37