我解这个微分方程组解了七个多月还没解出来！

siyanger · 发表于 2006-7-7 06:57

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

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

x

有做过星型或行星齿轮动力学方面的数值计算的吗？
我解一个行星微分方程组解了七个多月还没解出来！总得不到收敛的结果，查了资料说可能是重频率导致不收敛？不知道哪位研友做过？
要不然我都毕业了，真得很伤心哪！

[ 本帖最后由 siyanger 于 2006-7-7 07:48 编辑 ]

toes · 发表于 2006-7-7 09:01

你用什么方法求解的啊？
要不把方程组具体形式传上来大家看看。

7个月。。。:L 楼主不会换个方向啊？

siyanger · 发表于 2006-7-7 09:34

换方向了，近似解析法，数值法，matlab,fortrain都用过了，结果都是有的，可是都是发散的，不收敛．我都快疯了！

siyanger · 发表于 2006-7-7 09:49

方程的形式就是一般的微分方程：ＭＸ＇＇＋ＣＸ＇＋ＫＸ＝Ｆ（t）
这里的ＭＣＫ为１６Ｘ１６的矩阵。Ｍ是对角阵，Ｃ，Ｋ为对称阵。
如果ＭＣＫ都为定常的话，用wilson直接积分法可以算出结果来．
可是因为要做的是Ｋ为时变的（时间的显函数），弹性力那一项为ＫＦ（Ｘ），其中Ｆ（Ｘ）为非线性函数（分段函数），所以想用数值积分法算一下，但数值法现在线性都出不来，非线性更别说了。

有人说是因为系数病态，有人说是因为系统有重频，所以不收敛。

用ＭＸ＇＇＋ＫＸ＝０（ＭＫ为定常矩阵）算系统的固有频率，１６个自由度有３组重频．有时改变一下Ｋ中元素的系数，还会有固有频率是虚数的情况，这样就更不可能了啊，无阻尼系统的固有频率总是大于等于０的吧？

zhuofeng · 发表于 2006-7-7 11:38

同情中，感觉还是方法不对头

jbs0645 · 发表于 2006-7-9 12:54

原帖由 siyanger 于 2006-7-7 09:49 发表
方程的形式就是一般的微分方程：ＭＸ＇＇＋ＣＸ＇＋ＫＸ＝Ｆ（t）
这里的ＭＣＫ为１６Ｘ１６的矩阵。Ｍ是对角阵，Ｃ，Ｋ为对称阵。
如果ＭＣＫ都为定常的话，用wilson直接积分法可以算出结果来．
可是因为要做 ...

我想问一下，刚度是时变的，而x也是时间的函数，所以我觉得尽管k可能一直是大于0的函数在微分方程的求解过程中也是有可能求不出频率的吗？楼主觉得呢？

[ 本帖最后由无水1324 于 2008-3-5 20:21 编辑 ]

siyanger · 发表于 2006-7-9 16:47

楼上的，我没懂你什么意思。你说得频率是指固有频率吗？固有频率也可求的。
我现在用数值法求K是定常的也没求出来昵（直接积分法可以求线性的，但非线性的不好做，所以想用RK法等数值积分法做一下。），别说K是时变且F（X）是非线性的了。换了好几种方法，花了N多时间，也没能成功。不知道问题出在哪里。

fineshang · 发表于 2008-3-2 20:41

请问，现在问题解决了么？我刚刚开始，也遇到了类似问题，也是百思不得其解。欢迎讨论

lin20080111 · 发表于 2008-3-5 16:10

本帖最后由 VibInfo 于 2016-5-5 14:30 编辑

原帖由 fineshang 于 2008-3-2 20:41 发表
请问，现在问题解决了么？我刚刚开始，也遇到了类似问题，也是百思不得其解。欢迎讨论

LZ一年前提的问题现在应该不会再回来看了
K或M时变的问题不能求固有频率，一般只能通过求解系统响应来分析问题
以前做过一个K或M时变气弹颤振问题，就是不能通过求固有频率求解临界马赫数，最后取不同的马赫数根据系统响应情况趋近了临界马赫数

无水1324 · 发表于 2008-3-5 20:22

请问这里可以用一些什么方法

WSYcxl · 发表于 2008-3-8 20:51

我也很想知道，碰到这种时变的问题，该如何求其频率。

meteor · 发表于 2008-3-12 23:13

可以试试龙格-库塔方法

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

[其他相关] 我解这个微分方程组解了七个多月还没解出来！

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

回复楼主的帖子

评分

回复 5楼的帖子

[其他相关] 我解这个微分方程组解了七个多月还没解出来！

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

回复 楼主 的帖子

评分

回复 5楼 的帖子

回复楼主的帖子

回复 5楼的帖子