Abaqus非线性迭代平均力和第二步的K怎么求的? 300

浏览:1041 回答:6

在研究非线性迭代如何做到和Abaqus完全一致,就做了一个简单的非线性弹簧,有两个地方Abaqus不一样,一直没有研究出来为什么?请教一下。

非线性弹簧本构关系如下:

SpringMat.png

也就是说U=1时,拉力为595,U=2时拉力为1128.

在弹簧右端加一个力为5000,初始增量为0.1

第一个增量步:

Inc1.png

查看msg,第一个迭代步得到的AVERAGE FORCE =188,第一个力增量=500,

第一个问题:那么这个188怎么从500得到呢?


第二个增量步:

第一个增量步一个迭代步就收敛了,第二个增量步开始

Inc2.png

Abaqus非线性迭代平均力和第二步的K怎么求的?的图4

在第二个迭代步,力增量为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结果。

Abaqus非线性迭代平均力和第二步的K怎么求的?的图5SpringNonLinearCAE.rar

两个问题能解决任何一个都非常非常感谢。

邀请回答 我来回答

全部回答

(3)
默认 最新
fan+CAE
大佬,这个问题有答案了吗,小弟目前也在查找相关的资料,可以分享下吗?
2020年9月8日
评论 点赞
xunai97

第二个问题是否是0.16*595+(0.92-0.16)*533=500

2020年6月18日
评论 1 点赞
刘阳宇东

x=0-1,y=595*x

x=1-2,y=533*x+62

第二个问题:力等于1000的时候对应的位移是(1000-62)/533=1.7598,而图中两步位移加起来为0.84+0.92=1.76,二者相等。这就是算出0.92的原因。

“那么Abaqus计算的K应该就是500/0.92=543.4783”这里的K是割线刚度了,拿割线刚度和切线刚度比当然是不对的。

2019年12月17日
已采纳 评论 2 点赞

没解决?试试专家一对一服务

换一批