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

动力学方程求出的值特大

请问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没给齐资讯, 楼上诸君都只能猜测:@)
页: [1]
查看完整版本: 动力学方程求出的值特大