行星轮系振动问题请教！

fineshang

我正在分析一个振动问题，列出的是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进行求解，最后输出
   （1）X与Y的轨迹图
   （2）Xsp(t)=Rs×φ2-Rp×φ3-Rc×φ7 与 X'sp(t)=Rs×φ'2-Rp×φ'3-Rc
×φ'7
   （3）Xcl(t)=Rc×(φ7-φ8) 与 X'cl(t)=Rc×(φ'7-φ'8)
   （4）Xsp(t) 与 X'sp(t)的相平面图
程序编好了，一运行，结果前几千次迭代全为零，后面就显示NaN，曲线也输不
出来，本人百思不得其原因何在，请教各位高手帮帮忙，给分析一下。是不是因
为没有量纲一化的缘故呢？？？？？？？？？？？
非常感谢！！！！
程序很简单，就是输入质量矩阵，刚度矩阵，阻尼矩阵等，调用一次ode45!

caichengtao

自己设计过行星轮（体积小，减速比大），但是对其特性还没有研究过！期待你的成果！

wenzihui

建议搂主把程序贴出来看看，这样无法作出判断：）

xautmcf

建议你看一下Lin Jian和G。Parker的论文

咕噜噜

程序编好了，一运行，结果前几千次迭代全为零，后面就显示NaN，曲线也输不出来？
把你的程序贴出来看看，没程序怎么帮你分析是怎么回事

无水1324

这个问题很棘手的，你最好检查一下参数方程，是否有点错误

yejet

初步估计是你的刚度矩阵，阻尼矩阵有问题

[ 本帖最后由 yejet 于 2008-10-14 20:25 编辑 ]

huaijuliu

LZ消失咯？ 请贴出原程序

ttwwooblueyes

这刚度矩阵似乎写起来挺费劲的，得好好检查下；
初始条件也得检查下；
特征值解过了吗？建议先解一下，对系统有个初步的了解再分析结果。如果系统数值刚度大的话也不一定要用ode45.