vumat杆件模型遇到了问题，前辈们能不能帮帮忙？ 100

我的模型载荷载是从x方向单向加载（荷载一直增加），一根简单的杆件模型 stateNew(i，1)=stateOld（i，1）+strainInc（i，1）这样计算出来的不是x方向的应变吗？ 模型一共四段直线，每一段都是按照线弹性来定义的，只不过弹性模量不一样 我用应力判断程序进入到哪一个区间，每个区间都是一个do-enddo循环 我的问题是： 在每个do循环我都会输出stressNEW（i，1）和stateNew(i，1)（有程序图），他们的值非常奇怪： 1、在开始计算的时候，二者都出现了负值，但是我是单项拉伸不应该出现负值！ 是vumat的计算机制的问题么？ 2、在计算的中间，我想输出1方向的应变：stateNew(i，1)，可是数值一直不对，要么是0，要么是负值 让我很崩溃（有log结果图） 另外：在计算的时候，应力直接从0跳到10^7级别，我的增量步设置的很小，应该怎么控制应力变化？ 虽然我需要的数据数值不对，但是 ABAQUS+vumat计算的结果是对的（有图） 我应该怎么调整我的程序？虽然程序不难，但是作为编程小白，很吃力，在此请教前辈们了！！！