chunhuo420 发表于 2007-6-6 08:18:59

如何移植弧长法到我的问题上

有限元中,有个弧长法能够追踪 荷载和位移曲线,并能得到极限荷载,
而我做的,如图,也是类似 荷载-位移曲线,但不是荷载-位移关系,
不过P-a同样存在收敛问题,请问能不能通过 P-a 曲线,得到 P的极限值。。
目前的做法是,逐渐增加P,当a趋向于无穷大,或者 曲线曲率最大时,对应的P即为极限值,
但是,通过不断增加P,存在两个缺点,1、P过大时,可能不收敛;2、P过小,计算步数太多。
后来看了弧长法,好像能解决这个问题,请问大家如何将弧长法移植到我的问题上。

hillyuan 发表于 2007-6-6 12:32:56

Why don't you 逐渐增加a,当P不变时对应的P即为极限值. It's much simple

chunhuo420 发表于 2007-6-6 14:12:44

原帖由 hillyuan 于 2007-6-6 12:32 发表 http://www.simwe.com/forum/images/common/back.gif
Why don't you 逐渐增加a,当P不变时对应的P即为极限值. It's much simple


谢谢你的回答,a是由改变P计算得到的,
并不能先改变a

chunhuo420 发表于 2007-6-7 15:46:43

我顶

baibing 发表于 2007-6-8 21:24:52

建议采用位移控制法 DCM

happysimwe 发表于 2007-6-24 11:07:13

在计算非线性问题的时候,单一的载荷控制或者位移控制都不是有效的。弧长控制能够极好地稳定程序运行。
1) 计算当前的切线刚度阵KT,给出初始dP
2) 位移增量da = KT dP
3) 计算弧长 l , 是||da||, ||dP||的无量纲化后的函数
4) 根据设定的弧长来线性调整da, dP
5) 重新迭代,必要时更新KT

适当地预测可以提高效率。弧长控制法可以跟踪屈曲过程。
页: [1]
查看完整版本: 如何移植弧长法到我的问题上