带初始不规则间隙分布的接触问题

想分析受力后轧辊的变形情况，但中间辊带有特殊的三次方曲线（如下图，轧辊为360度回转体，下图仅为轧辊中间截面），引起接触面不规则的间隙分布。

如果不调整从面节点坐标，运算时总是不收敛；但要是运算时调整从面节点坐标，辊形就没意义了，不知道怎么样在不调整中间辊辊形的基础上，使接触稳定的建立起来，并能运算收敛。
注：上图中间辊辊形夸张化了，实际中间辊辊形的辊径差最大值仅为0.591325mm。
谢谢大家的帮助！~
inp文件：