求助：多自由度微分方程的求解

luoluo · 发表于 2007-6-28 18:43

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

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

x

我建立了一个四自由度的振动方程:
[M]{x''}+[C(t)]{x'}+[K(t)]{x}={F}
其中的刚度和阻尼都是时变的，对于这类方程，该如何求解呢？
用matlab自带的ode45可以求解吗？
哪位朋友有类似方程的matlab源码啊,可否共享一下，先行谢过了！

风花雪月 · 发表于 2007-7-2 10:59

可以的，搜索论坛有很多例子

octopussheng · 发表于 2007-7-2 13:42

将方程写成一阶微分方程组的形式，然后用ode求解，论坛有很多这方面的例子，别太偷懒啊，呵呵呵呵

fineshang · 发表于 2007-8-14 20:57

我现在做一个行星减速器的振动分析，需要解二阶时变系数的振动微分方程组，类似楼上的，然后通过一定工具绘出时域、频域特性曲线等。哪位高手知道请明示，我没用过MATLAB。或加QQ12700731

linqus · 发表于 2007-8-16 07:13

楼主的问题用newmark－β等方法求解该是没问题的，
只要在每个time step中更新C，K即可。
建议尝试。

探讨。

fineshang · 发表于 2007-8-20 20:46

请问楼上的朋友，什么是newmark－β法，可以解决以上提出的问题么，推荐一本着方面的书给我吧，我急需解决这个问题，谢谢啦！！！也请做过这个问题的朋友与我联系，QQ12700731，多多沟通交流啊！！

风花雪月 · 发表于 2007-9-3 10:11

原帖由 fineshang 于 2007-8-20 20:46 发表

 登录/注册后可看大图

请问楼上的朋友，什么是newmark－β法，可以解决以上提出的问题么，推荐一本着方面的书给我吧，我急需解决这个问题，谢谢啦！！！也请做过这个问题的朋友与我联系，QQ12700731，多多沟通交流啊！！

大部分结构力学方面的书籍都会有这一算法的说明
比如王焕定的《结构力学》
刘晶波的《结构动力学》等等

luoluo · 发表于 2007-9-13 09:25

可以提供一个时变参数多自由度微分方程求解的范例吗？谢谢了

风花雪月 · 发表于 2007-10-12 09:47

原帖由 luoluo 于 2007-9-13 09:25 发表

 登录/注册后可看大图

可以提供一个时变参数多自由度微分方程求解的范例吗？谢谢了

时变参数和定参数求解的程序没有太大区别，增加个传递变量即可

无水1324 · 发表于 2007-10-12 10:51

将时变量直接写出就可以了。与其他常数一样

songjianjiang · 发表于 2007-10-12 23:02

将矩阵解耦,应该prony方法可以解决这个问题,matlab中办法就不知道了

风花雪月 · 发表于 2007-12-6 09:25

原帖由 songjianjiang 于 2007-10-12 23:02 发表

 登录/注册后可看大图

将矩阵解耦,应该prony方法可以解决这个问题,matlab中办法就不知道了

有很多方程解耦是很困难的

su200703 · 发表于 2007-12-7 23:21

原帖由 fineshang 于 2007-8-20 20:46 发表

 登录/注册后可看大图

请问楼上的朋友，什么是newmark－β法，可以解决以上提出的问题么，推荐一本着方面的书给我吧，我急需解决这个问题，谢谢啦！！！也请做过这个问题的朋友与我联系，QQ12700731，多多沟通交流啊！！

C++源代码：http://www.chinavib.com/forum/thread-49261-1-1.html

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

[其他相关] 求助：多自由度微分方程的求解

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

一样的问题请示

请教！

回复 #8 luoluo 的帖子