我正在尝试在ABAQUS中实现一个子程序。
这是一个非常简单的非线性弹性模型，其中杨氏模量取决于平均压力，详见E=3*(1-2*poisson)*p/kap（其中，poisson=0.3是泊松系数和kap=0.005膨胀指数）。对于sigma11,22和33，初始应力为1e5 Pa。
当我运行子程序时，它给出了E=3*(1-2*0.3)*(3*1e5/3)/0.05的线性行为（这是用初始应力计算的杨氏模量）。如果所有组件的初始压力为0，则所有计算都为0，因为E=3*(1-2*0.3)*(3*0/3)/0.05=0。
我想问你是否可以帮助我解决这个问题（将初始条件定义为每个变量的先前值）。