双剪弹塑性编的子程序为何单轴拉正常单轴压不正常啊？

lxw917

材料子程序编好了，验算单轴受拉正常（双线性应力应变关系），但受压不正常----看上去应力应变还是弹性关系。用的双剪强度理论，定义受拉屈服强度/受压屈服强度=0.1。 我自己编了个主程序用应变增量{-ε νε νε 0 0 0} 来调用调试，应力刚好超出屈服极限时流动矢量算出来是{-0.1 0.5 0.5 0 0 0}，这个矢量是不是不对？第一项是不是应该是-1.0，第二三项则要远远小于第一项的绝对值？算出来应该从总弹性应力增量中因塑性变形扣除的应力是一个正数，按理应该是个负数啊？哪里出了问题呢？ 好急啊，怎么也检查不出错误所在啊。