你在这里
Riks方法中的弧长因子
星期五,2012-07-13 14:42 -Staszkinson
你好,
如果有人能向我解释如何在Abaqus中校准Static, Riks方法中的弧长增量因子,我将不胜感激。
对于静态风险步骤中的默认输入,我得到的最大LPF因子等于147(这个值可能是正确的,预期为147.8)。不幸的是,我对增量的数量(12个增量的结果)不满意。问题是:如何校准弧长增量因子(初始、最小、最大弧长增量和估计的总弧长)以获得更准确的结果。
我发现了一些东西,原话是:
在严重非线性的情况下,Riks方法可能不收敛或收敛到先前得到的平衡构型。通常这种问题是由于弧长太大造成的。如果您怀疑这就是问题所在,请在消息文件中查找所使用的典型弧长,并以较小的最大弧长返回分析结果。
典型的弧长是多少?
附件包含lpf弧长图。
附件 | 大小 |
---|---|
graph.JPG | 28.69 KB |
论坛:
免费的标签:
弧长法相当
由于力和位移都是未知的,弧长法非常独特。
为了理解这个方案,想象一下你被蒙住眼睛,沿着一堵弯曲的墙行进(不允许偏离曲线)。所以,在任何时候,你拿起你的棍子,试着轻敲你的前面。在这种情况下,手杖的水平投影为搜索半径,一步所走的距离为弧长。
弧长也是如此。在abaqus中,您可以给出初始增量,总弧长(默认为1,大多数情况下都有效)和最大弧长增量(如果我没记错的话,默认为1E+36)。因此,指定最大弧长(我总是从1开始,并密切关注状态文件)。在状态文件中,将打印LPF和LPF的增量。如果你发现它太小,增加最大弧长增量。你可以把它解释为
LPF*应用负载=当前负载
LPF inc * Applied Load = Load增量。
希望它能帮助你达成目标。
库马尔
--
世界从0开始,以0前进,但不想要0。
谢谢你的
谢谢你非常好的解释。
我必须同意,我发现总弧长在很多情况下等于1是正确的值。
现在我可以自由地改变因子,我或多或少知道我在做什么。
Ł。
弧长法
请告诉我E Ramm的弧长法的算法。我正在使用以下类型的算法,但我没有得到正确的结果----
你{0},λ= 0
Fext
N个步数
Kg(刚度u{0})
u1=(Kg)^-1 * text '
[L D] =低密度脂蛋白(公斤)
m =侦破(D)
arc_length =规范(u1)
λ= 0.1 * arc_length /√(u1’* u1 + 1)
如果m < 0
λ=λ
结束
u =λ* u1
迭代:
公斤(u {0} + u)
弗林特(u {0} + u)
耐火的= Fint-lamda * Fext
你=发票(公斤)*耐火的”
* Fext u1 =发票(公斤)
del_lamda = - (u ' * ur) / (u ' * u1)
del_ur =你+ del_lamda * u1
u = = u + del_ur
λ=λ+ del_lamda
如果规范(即< 10 ^ 5)
结束
u {0} = {0} + u
如果λ> = 2
结束
所以请帮助我……急切地等待解决方案。
谢谢
Md Rushdie Ibne Islam
引用请求
亲爱的“irushdie”,
你们有拉姆的电子版的研究报告吗?如果有,你可以发邮件给我,以便我帮助你吗?我的邮箱是gpapazafeiropoulos@yahoo.gr。
提前非常感谢你。
乔治
弧长法
在以下链接中,可以找到实现弧长方法的MATLAB源代码(Crisfield, 1981和Fafard & Massicotte, 1993):
http://en.pudn.com/downloads599/sourcecode/math/detail2447809_en.html
http://www.mathworks.com/matlabcentral/fileexchange/44352-arc-length-method
_______________________________________________
乔治Papazafeiropoulos
中尉,基础设施工程师,希腊空军
土木工程师,硕士,博士研究生,台大
电子邮件:gpapazafeiropoulos@yahoo.gr
网站:http://users.ntua.gr/gpapazaf/