Contact Convergence Issue

[grnlan03]

I have been having trouble with a thermal-stress model that will not converge. I believe my problem is somewhere in the contact interaction I have applied to the model called “internal contact.” The model will converge with a low density mesh with the contact interaction on but it will not converge at a higher density mesh with the interaction on. It will also converge when the interaction is turned off for both the low density and high density mesh.

I have attached python script file that will create the higher density mesh model with the interaction on i.e. “Internal Contact.” I’m hoping someone can look at it and give some insights.

Thanks