Large deformation spherical (hertz) contact simulation

Hello Everyone,

I am simulating contact of a sphere with a soft elastomer (Hertz contact). Unlike the hertz contact, I will like to simulate large deformations. My model is a 2D axisymmetric model and the load is displacement controlled. The current model fails at total time/freq ~ 0.84. I have tried refining my mesh and have reduced the increment size and increased the number of increments. Currently, step my parameters are:

 Max. number of increments: 200000 and Increment size:0.005/1E-5/0.01

 Nlgeom is On because of the large deformation.

Please, does anyone know why the model is not running to completion? I will greatly appreciate any suggestion.

Thank you.

Uyi 

>> I have attached the inp file. I am only simulating the loading step in this case.

files/JKR_Contact.txt