在求解非线性动力系统周期解时,遇到刚性的系统怎么办呢?

胡晓宇 · 发表于 2009-3-23 15:22

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

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

x

大家有没有好的建议,或者方法,指点一些资料或者书也行.

讨论求知 · 发表于 2009-3-23 19:39

请问什么叫刚性系统，和非线性系统有什么区别

无水1324 · 发表于 2009-3-24 08:25

看一下 matlab的主页，找一下stiff ode，一般刚性方程可以用 ode15s，就是s结尾的算法进行计算

胡晓宇 · 发表于 2009-3-24 10:17

这里的刚性系统就是方程Jacobian矩阵的条件数很大，为病态，在用newton-raphson求解时，解得变化很大。

胡晓宇 · 发表于 2009-3-24 10:23

谢谢无水的指点。刚性系统的积分是可以的，比如隐式法或者你说的matlab中s的算法。但是，求周期解是一个初值问题，积分一个周期后，在用newton-raphson方法求解时，由非线性方程组得到的解往往在各个维上的数值差别比较大，那么下一步迭代就不容易处理了，往往直接就发散了。是由于我的方程无量刚化有问题么？

无水1324 · 发表于 2009-3-24 12:17

对方程量纲处理，一个主要的目的就是降低方程的刚性，你再仔细看一下量纲一化是否出了问题

胡晓宇 · 发表于 2009-3-24 13:48

我做的是滚动轴承-转子系统的动力学分析.无量刚化选取两个变量:长度是以轴承游隙;时间就是wt

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

[稳定性与分岔] 在求解非线性动力系统周期解时,遇到刚性的系统怎么办呢?

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

回复楼主胡晓宇的帖子

回复沙发讨论求知的帖子

回复板凳无水1324 的帖子

回复 5楼胡晓宇的帖子

回复 6楼无水1324 的帖子

[稳定性与分岔] 在求解非线性动力系统周期解时,遇到刚性的系统怎么办呢?

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

回复 楼主 胡晓宇 的帖子

回复 沙发 讨论求知 的帖子

回复 板凳 无水1324 的帖子

回复 5楼 胡晓宇 的帖子

回复 6楼 无水1324 的帖子

回复楼主胡晓宇的帖子

回复沙发讨论求知的帖子

回复板凳无水1324 的帖子

回复 5楼胡晓宇的帖子

回复 6楼无水1324 的帖子