“…and an initial crack set S u(⋅,t 0 ) . To render the model Equation (1) well-defined for each instant of time, the displacement field u needs to be chosen in a suitable function space of discontinuous functions, see Chambolle et al 51,52 To treat the Francfort-Marigo model numerically, it is customary to discretize the quasi-static evolution 53 in time by a backward Euler scheme, that is, to subdivide the time interval [t 0 , t 1 ] into K+1 increasing time instants t 0 = t 0 , t 1 , t 2 , … , t K−1 , t K = t 1 and to minimize, for each k = 1, 2, … , K, the functional (1) subjected to the boundary condition u k = u 0 (⋅, t k ) on Ω and the constraint S u k ⊇ S u k−1 , where u k and u k−1 denote minimizers of the functional (1) corresponding to the time steps t k and t k−1 , respectively.…”