| 查看: 3932 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | ||
zkunF金虫 (小有名气)
|
[求助]
动力有限元隐式求解,牛顿迭代法不收敛
|
|
|
大家好,现有一问题求助各位虫友,希望大家帮忙解决,先行谢过。 问题如下: 在编写动力有限元程序时,运用牛顿迭代法求解动力方程,具体过程参照了庄茁老师译著《连续体和结构的非线性有限元》。 !------------------------------------------- 基本方程: K△d = r ; d(v+1) = d(v)+△d (v为迭代次数) 收敛准则: r的2范数残差准则,记Res_0=||r(0)||2;Res_v=||r(v)||2; 收敛判断——相对误差Res_v/Res_0<ERROR 在求解过程中,前两步算完,相对误差大致与时间步长一个量级,再进行迭代,相对误差基本保持不变,也就是说参数r(v)基本不变,每一迭代步的△ d变化很小,但并不逼近0. 例如,时间步长取1.0e-4 Iteration step=1 相对误差 = 2.0e-3 Iteration step=2 相对误差 = 3.0e-4 . . . Iteration step=v 相对误差 = 3.0e-4 另,当把时间步长调小后(小2个数量级以上),也是前两步算完,相对误差就基本保持不变了。 !------------------------------------------- 这个问题困扰蛮久,也看过数值分析关于牛顿法在某些情况下的收敛问题,尝试过修正牛顿法等,效果并不好。切刚K是按照参考书上推导而来的,检查过应该没有问题。 希望万能的虫友们不吝相助,谢谢! |
» 猜你喜欢
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
为什么资助数各大高校都创新高,自己申请怎么就这么难
已经有15人回复
科研人应该花精力去思考如何解决问题,而不是去凝练问题
已经有14人回复
要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够
已经有17人回复
基金系统什么内容也没有
已经有10人回复
学科评审组评审是指会评吗?
已经有5人回复
面上合作单位盖章
已经有9人回复
面上函评意见出来了,像什么等级?
已经有17人回复
小白求助 投论文要求的highlights应该如何写
已经有4人回复

3楼2018-10-07 13:41:45
|
本帖内容被屏蔽 |
2楼2016-10-13 18:17:56
4楼2023-03-08 19:11:17










回复此楼
想问问你解决了吗?
30