you are correct stabilization was decreased and the unexpected deformation was minimized. Also, interestingly enough probing # contacting (in the actual model) showed that the material displaces after number contacting changes to 1 which is exactly as expected. Implementing this in the actual complicated model involved increasing the number of steps but still worked. I have high hopes of success.

I have (in the past) implemented an "artificial" gravity field, I did and do not understand what this will actually do. What does it add to the solver or can you briefly explain this idea a little for me?