UMAT定义的正交异性弹塑性本构模型，为何独立单元测试通过但模拟多单元时会出现不同程度的收敛困难？

一、理论方面

正交各向异性材料的基于UMAT的弹塑性本构模型二次开发：1）弹性阶段，正交异性刚度矩阵；2）塑性阶段，Hill强度准则，各向同性指数硬化模型。

二、程序方面

1)弹性试算；2）塑性修正，通过N-R迭代法。

三、遇到的问题

在模拟边长为 1mm 的 立方体块 时，首先基于独立单元进行调试，各个正交方向单轴加载均没有问题；但是随着单元数量的增加（4个单元、125个单元、1000个单元），程序的收敛性越来越差，

如果不采用UMAT，而是将相同的CAE模型（几何、边界等），将强化曲线座位塑性数据直接输入到ABAQUS的塑性定义中，对于不同数量的单元数量则具有相同的收敛性：

四、提问

请问各位老师，为什么使用UMAT的时候，在独立单元可以调试通过，可以模拟到应变达 0.30，而在单元较多的时候，却无法模拟到最终，会提前终止，请问这是什么原因造成的