声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 7315|回复: 25

[编程技巧] newmark积分算法步长选择

[复制链接]
发表于 2016-4-15 11:24 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
求解,最近做非线性这一块,用newmark积分算法求解方程,步长选择0.1和0.01 出来的相图差异很大,很困惑,请教newmark算法积分步长的选择的原则是什么有没有有经验的小伙伴,在网上搜索了挺久也没能搜到。 求助

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2016-4-15 13:13 | 显示全部楼层
转载一篇文章看看:

前段时间做了一个瞬态动力学的项目分析,对于时间步长的选择也进行过一段时间的研究,现把我学到的一些资料和自己的体会融合贴出,请大家指教。
一、瞬态动力学分析包含数值为时间函数的载荷,在载荷/时间曲线上的每个“拐角”都应该作为一个荷载步。瞬态分析求解的精度决定积分时间步的大小,要想计算出最优时间步长,应当遵循下列5个准则:
1、解算响应频率时,时间步长应当足够小以能求解出结构的响应。
结构的动力学响应可以看做是各阶模态响应的组合,时间步长应当小到能解出对整体响应有贡献的最高阶模态。对于NEWMARK时间积分,已经发现当时间步长取值20倍最高频率时会产生比较合理精度的解,即积分时间步长ITS=1/20F。
2、解算所加荷载/时间关系曲线时,时间步长应当小到足以跟随载荷函数
响应总是倾向滞后于施加的载荷,特别是对于阶跃载荷。要跟随阶跃载荷,ITS要小到和1/80F相近。
3、解算接触频率时。
在设计接触(碰撞)问题时,时间步长应当小到足以捕捉到两个接触表面之间的动量传递,如果违反,将发生明显的能量损失,从而碰撞将不是完全弹性的,积分时间步长可由接触频率FC确定
ITS=1/NFC
FC=(K/M)**0.5/(2*PI)
K 是间隙刚度,M 是作用在间隙上的有效质量,N 是每周的点数,要使能量损失最小,每周至少要取30个点,即N=30;如果要得到加速度结果,可能N取值更大些;对于缩减法和模态叠加法,N至少需要7才能确定求解的稳定性。
4、解算波的传播时,则时间步长应当小到当波在单元之间传播时足以捕捉到波,所以要根据单元大小来确定步长。
5、解算非线性时,大部分问题要求满足前面四个准则的时间步长就能捕捉到非线性行为。
还有少数例外的情况,例如当结构在载荷作用下趋于刚化,则必须求解被激活的高阶模态。
二、自动时间步长。
在用合适的准则计算出时间步长后,在具体分析中应该用最小的步长值。当然可以采用自动时间步长来让ANSYS决定在求解中什么时间增大或减小时间步长,这样的好处是可以减少子步的总数,还可以大大减少可能需要重新分析的次数,如果有非线性,还可以适当增加载荷并在达不到收敛时回溯到先前收敛的解(二分法),可以用AUTOTS激活自动时间步长。要避免使用极小的时间步长,尤其是在建立初始条件时,这样可能会导致数值障碍。
在下列情况不适宜用自动时间步长
1、只是在结构的局部有动力学行为的问题。这时系统部件的低频能量部分远远高于高频部分。
2、受恒定激励的问题,如地震荷载。在这种情况下当不同频率被激活时时间步长趋于连续变化。
3、运动学问题,在这种情况下刚体运动对响应频率的贡献占主导地位。
 楼主| 发表于 2016-4-15 14:59 | 显示全部楼层
glise 发表于 2016-4-15 13:13
转载一篇文章看看:

谢谢!  
 楼主| 发表于 2016-4-15 15:01 | 显示全部楼层
glise 发表于 2016-4-15 13:13
转载一篇文章看看:

老师您有没有经验用newmark 算过类似方程 步长是如何确定的

点评

一般newmark法步长主要考虑收敛性比较多 做响应计算式需要满足Nyquist采样定理  详情 回复 发表于 2016-4-15 15:40
发表于 2016-4-15 15:40 | 显示全部楼层
a9303 发表于 2016-4-15 15:01
老师您有没有经验用newmark 算过类似方程 步长是如何确定的

一般newmark法步长主要考虑收敛性比较多

做响应计算式需要满足Nyquist采样定理
 楼主| 发表于 2016-4-15 16:06 | 显示全部楼层
frogfish 发表于 2016-4-15 15:40
一般newmark法步长主要考虑收敛性比较多

做响应计算式需要满足Nyquist采样定理

老师我的两个积分步长求解出来都是收敛 是不是积分步长越小就越精确点,因为我的程序是循环步长选0.01的话要运行5个多小时

点评

正常情况下应该是这样的,所以步长的选择还需要结合计算效率来考虑  详情 回复 发表于 2016-4-15 16:09
发表于 2016-4-15 16:09 | 显示全部楼层
a9303 发表于 2016-4-15 16:06
老师我的两个积分步长求解出来都是收敛 是不是积分步长越小就越精确点,因为我的程序是循环步长选0.01的 ...

正常情况下应该是这样的,所以步长的选择还需要结合计算效率来考虑
 楼主| 发表于 2016-4-15 16:37 | 显示全部楼层
感谢!!
发表于 2016-7-13 22:57 | 显示全部楼层

楼主肯定是做非线性。

我也遇到了,积分步长不够细的话,甚至会发散!

我也没有办法,这个玩意儿真难。
发表于 2016-7-14 09:14 | 显示全部楼层
楼主可以参考一下论坛里这个专题:
http://forum.vibunion.com/forum. ... d=53&fromop=all
 楼主| 发表于 2016-7-14 14:53 | 显示全部楼层
mxlzhenzhu 发表于 2016-7-13 22:57
楼主肯定是做非线性。

我也遇到了,积分步长不够细的话,甚至会发散!

对的 的确在做非线性这块,理论上步长越小精度越高,但是还要考虑运算时间 。这两点要综合下,
 楼主| 发表于 2016-7-14 14:54 | 显示全部楼层
minoz 发表于 2016-7-14 09:14
楼主可以参考一下论坛里这个专题:
http://forum.vibunion.com/forum.php?mod=collection&action=view&cti ...

谢谢,关注过这个专题。

点评

基本是常见问题 里面都有  详情 回复 发表于 2016-7-14 15:41
发表于 2016-7-14 15:41 | 显示全部楼层
a9303 发表于 2016-7-14 14:54
谢谢,关注过这个专题。

基本是常见问题  里面都有
发表于 2016-7-24 10:07 | 显示全部楼层
LZ可以关注下 显式求解 和 隐式求解的区别,
貌似Newmark属于隐式求解法, 只要能收敛, 步长不是问题...
发表于 2016-7-25 09:03 | 显示全部楼层
我也有同样的问题,最近做非线性的时候,发现步长的不同队相图的影响很大,那么到底该如何选择呢

点评

sd
如果用n-r迭代呢  详情 回复 发表于 2016-11-15 08:40
我估计吧Newmark用来做土木工程比较好,目前MATLAB里面ode45都是联合四阶和五阶RK法,可能这就是原因吧;很多非线性振动问题终究也还是简单的ODE。  发表于 2016-10-6 16:16
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-3-29 22:15 , Processed in 0.066185 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表