一个简单的非线性弹性的umat不收敛 请大佬指教

清欢

大佬们，我最近在写一个图里的子程序，因为是弹性的嘛， z方向应力我就按照公式那样更新 只不过用STRAN(2)+DSTRAN(2)代替公式里z方向的正应变 DDSDDE中的DDSDDE(3，3)我是把公式关于z方向正应变求导，再使用STRAN(2)+DSTRAN(2)代替求导结果里z方向的正应变，这样做法应该是对的吧？

但是现在遇到一个问题。写完了子程序我想看看对不对嘛，就在abaqus中以单压的形式进行仿真，获得应力应变曲线来对比，当我建出来上方压头并用压头位移为边界条件时，仿真可以收敛，在计算出大概的力以后，我尝试着去掉压头，在构件上表面施加压力，大小由上述力换算得到，这时候发现不收敛了。我一开始以为是ddsdde在应变太小的时候有问题，因为压头位移为边界条件时，应力应变和公式曲线一致的，说明应力更新不会出错。为了收敛，我将umat中应变0.01之前改为线性弹性段，但是abaqus算到0.01之后就报错 还是不收敛。想知道这个问题怎么解决file:///C:/Users/Administrator/Desktop/fiberbed/%E6%9C%AC%E6%9E%84%E5%85%B3%E7%B3%BB.jpg

清欢

大佬们麻烦来看一下