Abaqus非线性迭代平均力和第二步的K怎么求的? 300
浏览:1041 回答:6
在研究非线性迭代如何做到和Abaqus完全一致,就做了一个简单的非线性弹簧,有两个地方Abaqus不一样,一直没有研究出来为什么?请教一下。
非线性弹簧本构关系如下:
也就是说U=1时,拉力为595,U=2时拉力为1128.
在弹簧右端加一个力为5000,初始增量为0.1
第一个增量步:
查看msg,第一个迭代步得到的AVERAGE FORCE =188,第一个力增量=500,
第一个问题:那么这个188怎么从500得到呢?
第二个增量步:
第一个增量步一个迭代步就收敛了,第二个增量步开始
在第二个迭代步,力增量为500,位移增量为0.92,那么Abaqus计算的K应该就是500/0.92=543.4783。
但这个与本构关系不一致。
因为位移肯定在2以内,2以内的本构关系算出的K为
0-1s内,K=595,
1-2s内,K=1128-595=533。
第二个问题:Abaqus的K与595和533都不一致,它怎么算出来的?
附上6.12-1版本的cae文件和msg,注意,无法分析成功,但通过msg可以查看前两个分析步的Abaqus结果。
两个问题能解决任何一个都非常非常感谢。