bryant 发表于 2008-5-13 10:17

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

本人需要做jeffcott裂纹转子的振动响应分析，首先要分析无裂纹时的转子振动响应，但得出的图形是一个衰减的图形，跟理论不相符，请大家指点指点。
微分方程：
X''+(2ξ/Ω)X'+(1/Ω^2)X=(1/Ω^2)+ecosωt
Y''+(2ξ/Ω)Y'+(1/Ω^2)Y=esinωt
主函数：
global e xi Omega omega;
e=0.1;xi=0.01;Omega=0.5;omega=100;
=ode45('odex',,);
=ode45('odey',,);
plot(x(:,1),y(:,1));
子函数：
function dx=odex(t,x)
global e xi Omega omega
dx=zeros(2,1);
dx(1)=x(2);
dx(2)=e*cos(omega*t)+1/Omega^2-2*xi/Omega*x(2)-1/Omega^2*x(1);

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

whecust 发表于 2008-5-27 16:48

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

gh688 发表于 2008-5-27 16:56

回复 2楼 的帖子

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

whecust 发表于 2008-5-28 09:48

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

whecust 发表于 2008-5-29 09:20

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

gh688 发表于 2008-5-29 22:30

回复 5楼 的帖子

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

whecust 发表于 2008-5-31 15:50

回复 6楼 的帖子

固有频率是0.5？可能lz这个值太小了吧
页: [1]
查看完整版本: 求助：关于jeffcott转子振动响应