|
正在做行星轮系的动特性分析!求助!
我正在分析一个振动问题,列出的是10自由度二阶微分方程组,形如:
[m]{X''}+[C]{X'}+[K]{X}={P}
其中[m]=diag{J1,J2,J3,J4,J5,J6,J7,J8,m,m}
{X}=(φ1,φ2,φ3,φ4,φ5,φ6,φ7,φ8,X,Y)
[C]为10×10阶的常数矩阵
[K]为10×10阶的时变矩阵,其中每个元素均是时间t的函数,Kij=Kij(t)
{P}=diag{Tin,0,0,0,0,0,0,0,0,-Tout}
想用数值解法(四阶龙格库塔法),借助matlab进行求解,程序编好了,一运行,结果前几千次迭代全为零,后面就显示NaN,曲线也输不出来,本人百思不得其原因何在,请教各位高手帮帮忙,给分析一下。是不是因为没有量纲一化的缘故呢???????????
非常感谢!!!!
程序很简单,就是输入质量矩阵,刚度矩阵,阻尼矩阵等,调用一次ode45!
刚度数量级为8次,阻尼数量级为3次,质量矩阵数量级为 —3 次,这个有问题吗?
高手能不能跟我联系,我好讨教下,谢谢!QQ12700731,请注明“‘振动’‘! |
|