动力学方程求出的值特大

sytibm · 发表于 2012-3-22 16:44

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

x

请问m*d2x+k*x=q这个方程k是关于x的函数，可是用龙格库塔求出x的值太大了，都10的几百次方了？谁知道错误在哪啊？谢谢。

VibrationMaster · 发表于 2012-3-22 19:53

检查k是否会随x变大而变成负的. 那样的话,方程有问题

玉林 · 发表于 2012-3-23 20:01

方程不收敛？

sytibm · 发表于 2012-3-27 15:26

回复 2 # VibrationMaster 的帖子

转换后d2x=m'*q-(m'*k)*x；在matlab中m'*k运行几步会出现NaN

玉林 · 发表于 2012-3-31 21:59

方程不收敛，先用特征值计算，看看是不是有的根的实部大于0.

dw04116 · 发表于 2012-4-5 16:32

这类数值算法，得如履薄冰，我在用newmark法解个非线性问题，各种问题纠结中……

无水1324 · 发表于 2012-4-6 18:21

方程不收敛，应该是参数有问题，注意阻尼的符号

ChaChing · 发表于 2012-4-13 00:14

回复 1 # sytibm 的帖子

个人水平有限, 但感觉LZ没给齐资讯, 楼上诸君都只能猜测:@)

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

[非线性振动] 动力学方程求出的值特大