Hi Maria,
You probably haven't change any parameters in the model. (rho_ref, psi_ref, young_modulus) Originally, they were chosen only to assure convergence (which is poor anyway), they have no actual meaning. You can tweak them yourself for particular application. The results you are producing are fine, they just need better postprocessing, because unloaded discontinuity results in negative density. We should further look into it.
Load the attached file in paraview to see how it looks with IsoVolume filter.
You can load file in preview: In top menu, you have File -> Load State.