有哪位大神做过梁单元的混凝土材料子程序UMAT啊

我是大四一毕业生，毕业论文就是这个题目。梁单元简化了很多过程。UMAT的编制也本质上就是给定一个应变，求出曲线上相应的应力，以及该点的切线。
我用一个简单的模型来测试子程序。就是一个直杆，一端固定，另一端位移加载。然后做某一点的应力应变曲线，发现曲线上某些地方相邻两点的应力相同，但是应变不同，这点很难解释为什么会出现这种情况。
不知道是子程序的位移，还是建模过程中的参数设置有问题。
希望坐过站这方面的大神能指点一二。
下面是应力应变曲线（我是把数据提取出来在origin里面画的）局部放大图，横轴是应变轴，纵轴是应力轴。