In a two dimensional setting, we present a result of convergence of an alternate minimization scheme applied to a phase field model of fracture with non-interpenetration. Our analysis is based on the study of suitable gradient flows of the phase field energy, which connect all the states of the algorithm. The limit evolutions are described in terms of parametrized BV-solutions. This is a joint work with M. Negri.