动态非线性的拓扑优化

大家好。我现在用optistruct11.0.201，做动态非线性的拓扑优化。用ESLM方法将动态非线性转化成静态线性，目前动态载荷、材料非线性的调试都已经算通，可是一旦加入接触，优化就不能计算了。用的求解器是nlgeom。出现的错误提示是：内部程序出错，“in file "tabscr.F",at location#1022”