The paper presents a numerical implementation of a novel approximation scheme for Griffith’s theory of brittle fracture recently proposed in Schmidt et al. (2009). The total potential energy of a brittle body (including bulk and surface terms) is variationally approximated by a family of functionals, depending on a small penalty parameter ". The (two-field) approximating functionals have as arguments the displacement field u and an eigendeformation field. The latter describes the regions of the body under high strains, where fracture will occur. Gamma-convergence of such a family of functionals to Griffith’s energy has been proved in Schmidt et al. (2009). Here we investigate numerical examples for the quasi-static crack propagation in mixe...