fineshang 发表于 2007-8-31 16:28

行星轮系振动问题请教!

我正在分析一个振动问题,列出的是10自由度二阶微分方程组,形如:
{X''}+{X'}+{X}={P}
其中=diag{J1,J2,J3,J4,J5,J6,J7,J8,m,m}
    {X}=(φ1,φ2,φ3,φ4,φ5,φ6,φ7,φ8,X,Y)
    为10×10阶的常数矩阵
    为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 发表于 2007-9-3 00:08

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

wenzihui 发表于 2007-9-3 20:57

建议搂主把程序贴出来看看,这样无法作出判断:)

xautmcf 发表于 2008-10-10 19:16

楼主!!

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

咕噜噜 发表于 2008-10-11 10:35

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

无水1324 发表于 2008-10-11 21:01

回复 5楼 咕噜噜 的帖子

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

yejet 发表于 2008-10-14 20:24

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

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

huaijuliu 发表于 2009-12-25 18:47

LZ消失咯? 请贴出原程序

ttwwooblueyes 发表于 2009-12-25 22:45

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

402144999 发表于 2014-7-26 22:31

页: [1]
查看完整版本: 行星轮系振动问题请教!