Vumat梁单元计算后材料的力学响应为什么与内置本构不一样？ 10

算例是金属棒被套在一段由梁单元组成的环内部，对金属棒施加Y方向位移，对环产生力直至拔断环。

金属棒是离散刚体，除开U2方向位移外的其他一切自由度均为固定，环是B31梁单元组成的，环是弹塑性本构，想要模拟塑料丝被金属拉断的效果，左边是我自己编写的弹塑性vumat，右边是abaqus内置弹塑性本构，现在的问题就是子程序与自带本构力学响应完全不一样，这两个算例只改了本构，其他边界条件，接触等等都没变，这是什么原因导致的呢？