我正在尝试在 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.
我想问您是否可以帮助我解决这个问题（将初始条件定义为每个变量的先前值）。