求救，非线性弹簧模拟螺栓！！

十万火急，请大侠救命啊！！！！！:t(（Shell model)我尝试了两种方法：
1. 首先我用参考点RP 确定了螺栓的位置，然后用弹簧连接各螺栓点。每个螺栓用两个弹簧表示，自由度分别是（1,1）和（2,2）。 RP和面耦合，各组件设置了Hard contact。 这种方法可以工作，但是老师说不准确
2. 这是老师的方法。 与上不同的是，我在INP文件里把RP改成了面上的节点。然后用非线性弹簧连接它们，如上。不考虑接触。但是这种不能收敛。
小弟附上了个简单文件的inp文件，希望大神们赐教！！！！！！！！！！！！！！