I want to model radial crack growth using XFEM under contact load. I defined Crack under contact area and Selected XFEM from interaction module. I specify contact property and I selected Fracture criterion and “VCCT “with mixed mode “BK”. I calculated critical energy release rate for my material. Crack is located at glass material and I put 6 for mode I and Mode II and I put for Mode III (1) and exponent (n=1); however, when I run contact modeling my crack doesn’t propagate and abaqus gives an error when job was running. Would you please help me to model crack propagate under contact load.