求助：关于jeffcott转子振动响应

bryant

本人需要做jeffcott裂纹转子的振动响应分析，首先要分析无裂纹时的转子振动响应，但得出的图形是一个衰减的图形，跟理论不相符，请大家指点指点。
微分方程：
X''+(2ξ/Ω)X'+(1/Ω^2)X=(1/Ω^2)+ecosωt
Y''+(2ξ/Ω)Y'+(1/Ω^2)Y=esinωt
主函数：

1. global e xi Omega omega;
   
2. e=0.1;xi=0.01;Omega=0.5;omega=100;
   
3. [t,x]=ode45('odex',[0:0.05:50],[0 0]);
   
4. [t,y]=ode45('odey',[0:0.05:50],[0 0]);
   
5. plot(x(:,1),y(:,1));

复制代码

子函数：

1. function dx=odex(t,x)
   
2. global e xi Omega omega
   
3. dx=zeros(2,1);
   
4. dx(1)=x(2);
   
5. dx(2)=e*cos(omega*t)+1/Omega^2-2*xi/Omega*x(2)-1/Omega^2*x(1);

复制代码

1. function dy=odey(t,y)
   
2. global e xi Omega omega
   
3. dy=zeros(2,1);
   
4. dy(1)=y(2);
   
5. dy(2)=e*sin(omega*t)-2*xi/Omega*y(2)-1/Omega^2*y(1);

复制代码

whecust

lz请问你的Ω指的是什么意思啊？

gh688

他这里的Ω应该指的是转子的转速

whecust

如果Ω指的是转子的转速，那么ω代表什么那？

whecust

而且Omega=0.5，不太可能是转子的转速吧

gh688

按照楼主的公式，ξ—就是阻尼比，Ω—是sqrt(k/m)，w才是转速

whecust

固有频率是0.5？可能lz这个值太小了吧