请教多自由一阶非线性微分方程组数值解法?

咕噜噜 · 发表于 2007-11-28 14:38

用电脑抓屏，在画图板里面存为jpeg文件上传

nymy · 发表于 2007-11-28 15:04

Equations

[ 本帖最后由咕噜噜于 2007-11-28 15:07 编辑 ]

nymy · 发表于 2007-11-28 15:07

咕噜噜

看到我的方程了么?

其实除了14个要求的变量以外,要改变的参数就是A1到A7了.

不过现在即使用最基本的参数,最后得到的都是复数的解,不能用ode得到实数的解.

咕噜噜 · 发表于 2007-11-28 15:08

你的程序呢，给我看看

nymy · 发表于 2007-11-28 15:10

程序就是调用ode函数了,运行没有问题的.

nymy · 发表于 2007-11-28 15:11

现在的问题其实就是这个方程能不能用ode求解.

nymy · 发表于 2007-11-28 15:12

另外就是如何判断这个方程就是刚性方程了? 如果是刚性方程了,应该怎么办? 有没有类似的例子?

咕噜噜 · 发表于 2007-11-28 15:28

通常参数问题会造成刚度方程，但不可能一直都是刚性方程
没有你的程序和参数我不敢说

nymy · 发表于 2007-11-28 15:58

咕噜噜 · 发表于 2007-11-28 16:40

如果是刚性方程用newmark方法自己编写程序试试看

nymy · 发表于 2007-11-28 19:48

谢了！

无水1324 · 发表于 2007-11-28 22:20

还是改用其他的方法吧

nymy · 发表于 2007-11-29 14:39

不知道该用什么方法，郁闷死了，而且一点也不懂，各位指点下吧？

无水1324 · 发表于 2007-11-29 15:29

newmark，论坛有类似的帖子你可以看一下

nymy · 发表于 2007-11-29 16:06

---------------------------------------------------

newmark方法是针对多自由度的线性系统的简化方法。一般的结构振动问题都是线性的或弱非线性的，用此法可以做到“多快好省”。但是如果系统是高度非线性的，您就老老实实用Runge Kutta法吧。如果自由度很多，考虑一下动态子结构之类的缩减方法。
------------------------------------------------------

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

[计算数学] 请教多自由一阶非线性微分方程组数值解法?

Equations

回复 #26 nymy 的帖子

回复 #28 nymy 的帖子