|
楼主 |
发表于 2011-11-11 15:17
|
显示全部楼层
回复 1 # bingkafei7155 的帖子
clc
clear
y=[];
f=0.5:0.5:10; % 转速参数
k=0;
for i=1:length(f)
disp(f(i))
period=2*pi/9;
b=f(i); %转速的参数传递
k=k+1;
step=period/100;
tspan=0:step:500*period;
options=odeset('RelTol',1.0e-9);
y0=[0.1;0.2;0.1;0.1];
[t,y]=ode15s('jeff_fun118',tspan,y0,options,b);
plot(f(i),y(30000:100:end,1),'r.','markersize',1) %画分叉图
pause(1);
hold on
end
save y.dat y -ascii
函数 function dy=jeff_fun118(t,y,flag,W)
NB=9;a=0;R1=28.262;r1=18.788;S1=0;S2=0;cc=0.630;C=0;
for j=1:1:NB;
Sita(j)=2*pi*(j-1)/NB+t*r1/(r1+R1);
Deltasita(j)=(y(1)*cos(Sita(j))+y(3)*(sin(Sita(j))))*cos(a)-1;
if Deltasita(j)<=0
Deltasita(j)=0;
end
Sa=(Deltasita(j)).^(3/2)*cos(Sita(j));
Sb=(Deltasita(j)).^(3/2)*sin(Sita(j));
S1=Sa+S1;
S2=Sb+S2;
end
Kp=283.4;
Qx=Kp*S1;
Qy=Kp*S2;
dy(1)=y(2);
dy(3)=y(4);
dy(2)=-cc*y(2)/W+C*cos(t)-Qx./W^2+1/W^2;
dy(4)=-cc*y(4)/W+C*sin(t)-Qy./W^2;
t(end,1)
dy=[dy(1);dy(2);dy(3);dy(4)];
上面的有关系数都是参考一个文献上的数值并经过无量纲化处理过的,这是分叉图程序,但是和原文献的图差距很大,哪位大虾可否指导一下啊?万分感谢啊,参考文献是《转子-滚动轴承非线性动力学研究》韩宝财 |
|