请问用Newmark求解微分方程时初值是怎么确定的

wangrong631 · 发表于 2011-4-30 17:30

ansys里求解就可以选择newm法

qiandashuai · 发表于 2011-5-1 19:34

本帖最后由 qiandashuai 于 2011-5-1 19:51 编辑

对于转子动力学方程来说，即使考虑非线性油膜力，Newmark方法中初值的影响也不像其他一些系统那样那么明显。对于普通转子系统来说，相空间中平衡点、吸引子的类别和个数相对来说较少，尤其是当转子转速给定时。Newmark算法本身的稳定性也是比较好的，初值对算法本身的影响应该也比较小。
另外，不知道楼主为什么要把Newmark的算法参数取为0.6和0.45？出于何种考虑？

qiandashuai · 发表于 2011-5-1 19:55

另，看了一下楼主程序，plot(x(1,:),x(1,:))画出来的不是直线又能是啥

zhoujoe · 发表于 2011-5-2 09:12

回复 17 # qiandashuai 的帖子

老师你好，去这两个参数是为了让其收敛的
取0.5和0.25的时候结果发散

zhoujoe · 发表于 2011-5-2 09:50

回复 18 # qiandashuai 的帖子

这个地方是小失误，应该是plot(x(1,:),x(2,:))这个画出来的看着也是直线呀

hiv5 · 发表于 2011-5-4 21:42

你说的圆或椭圆，是指相平面图？

mayuanzhuo · 发表于 2011-12-19 21:02

回复 11 # meiyongyuandeze 的帖子

你好，请问用newmark求解非线性系统的话，应该怎么处理，怎么迭代，谢谢！

402144999 · 发表于 2014-7-26 18:34

提示: 作者被禁止或删除内容自动屏蔽

dollfish000 · 发表于 2014-11-1 00:10

虽然过去很多年，但是这样编程好累，我现在计算汽轮机转子，模化后接近300个单元1500个自由度，如果按照下载的文档里的方法，那得累坏。消耗5个体能下载……我手贱了

haobaobaos · 发表于 2014-12-26 11:18

学习了谢谢

账号		自动登录	找回密码
密码			我要加入

402144999 402144999 当前离线积分 22 头像被屏蔽	发表于 2014-7-26 18:34 \| 显示全部楼层提示: 作者被禁止或删除内容自动屏蔽
402144999 402144999 当前离线积分 22 头像被屏蔽
	回复支持反对使用道具举报显身卡

[动力学和稳定性] 请问用Newmark求解微分方程时初值是怎么确定的