Jeffcott转子振动微分方程求解
m*X''+c1*X'+k1*X=m*e*(Z'*Z'*cosZ+Z''*sinZ);m*Y''+c1*Y'+k1*X=m*e*(Z'*Z'*sinZ-Z''*cosZ)+m*g;
(J+m*e*e)*Z''+c2*Z'+K2*Z=m*e;
式中X,Y,Z为未知量,用ODE-45的话不知道怎么表示该方程啊?大家有什么办法解这个微分方程吗?急着用啊,帮帮兄弟!
[ 本帖最后由 venture 于 2007-1-1 10:14 编辑 ] 设
x1=x
x2=x'
x3=y
x4=y'
代入上述方程就能整理出
x'=f(x)形式的标准方程 x'=f(x)形式的标准方程
设
x1=x
x2=x'
x3=y
x4=y'
x5=z
x6=z'
x1'=x2
x2'=....x6' 怎么含有x6'呢 该怎么整理呢
回复 #3 zhubaowei7 的帖子
是不是需要对原方程进行化简 变系数非线性微分方程组,只能用数值解法,用ode45或者ode23吧回复 #5 咕噜噜 的帖子
他这里主要是右边的二阶导数怎么处理的问题 我看先直接写进去求解一下试试,看看能不能求解!直觉应该是可以的!呵呵回复 #1 huzhaolin 的帖子
以前看到这类转子好像没有出现这个问题回复 #1 huzhaolin 的帖子
把参数贴一下吧!
页:
[1]