gneiranjincrackok.m

内燃机转子仿真

    Ta1=2*Ta1+Ta1*Ta1;
  end
  T1=I+Ta1;
  for ix=1:2
      f(ix)=0;
  end
      
      sumf1=aa0;
      sumf2=aa0;
      sumf3=aa0;
      sumf4=aa0;
      sumf5=aa0;
      sumf6=aa0;
     
      for iw=0.5:0.5:48
sumf1=sumf1+aa11(2*iw)*cos(iw*w*(t(:,i)-tt0(1)+0.5*h+0.5*h*xi(ii)))+aa22(2*iw)*sin(iw*w*(t(:,i)-tt0(1)+0.5*h+0.5*h*xi(ii)));
sumf2=sumf2+aa11(2*iw)*cos(iw*w*(t(:,i)-tt0(2)+0.5*h+0.5*h*xi(ii)))+aa22(2*iw)*sin(iw*w*(t(:,i)-tt0(2)+0.5*h+0.5*h*xi(ii)));
sumf3=sumf3+aa11(2*iw)*cos(iw*w*(t(:,i)-tt0(3)+0.5*h+0.5*h*xi(ii)))+aa22(2*iw)*sin(iw*w*(t(:,i)-tt0(3)+0.5*h+0.5*h*xi(ii)));
sumf4=sumf4+aa11(2*iw)*cos(iw*w*(t(:,i)-tt0(4)+0.5*h+0.5*h*xi(ii)))+aa22(2*iw)*sin(iw*w*(t(:,i)-tt0(4)+0.5*h+0.5*h*xi(ii)));
sumf5=sumf5+aa11(2*iw)*cos(iw*w*(t(:,i)-tt0(5)+0.5*h+0.5*h*xi(ii)))+aa22(2*iw)*sin(iw*w*(t(:,i)-tt0(5)+0.5*h+0.5*h*xi(ii)));
sumf6=sumf6+aa11(2*iw)*cos(iw*w*(t(:,i)-tt0(6)+0.5*h+0.5*h*xi(ii)))+aa22(2*iw)*sin(iw*w*(t(:,i)-tt0(6)+0.5*h+0.5*h*xi(ii)));

end
f(3)=sumf1;
f(4)=sumf2;
f(5)=sumf3;
f(6)=sumf4;
f(7)=sumf5;
f(8)=sumf6;    

for ix=9:10;
    f(ix)=0;
end

%%%%%%%%%%%%%%%%%me

%%%%%%%%%%%%%%%%%%%%me
F=[o2;f];

vf=vf+we(ii)*T1*F;
end
vk(:,1)=T*vk+0.5*h*vf;    

x=vk(1:n,1);
%vk(1:n,1)=vk(1:n,1)-vk(1:n,1)/n;
 zhuanzhix=x';

%xyz(:,i)=vk(1:n);
%cxyz1d2e7(:,i)=vk(1:n);

if kic==0;
v(:,i)=vk(:);
vv(1:n,i)=inv(M)*(vk(n+1:2*n)-0.5*C*vk(1:n));
aaii(1:n,i)=inv(M)*(f-C*vv(1:n,i)-K*v(1:n,i));
elseif kic==1;
    v1(:,i)=vk(:);
vv1(1:n,i)=inv(M)*(vk(n+1:2*n)-0.5*C*vk(1:n));
aaii1(1:n,i)=inv(M)*(f-C*vv1(1:n,i)-K*v1(1:n,i));
elseif kic==2;
    v2(:,i)=vk(:);
vv2(1:n,i)=inv(M)*(vk(n+1:2*n)-0.5*C*vk(1:n));
aaii2(1:n,i)=inv(M)*(f-C*vv2(1:n,i)-K*v2(1:n,i));
elseif kic==3;
    v3(:,i)=vk(:);
vv3(1:n,i)=inv(M)*(vk(n+1:2*n)-0.5*C*vk(1:n));
aaii3(1:n,i)=inv(M)*(f-C*vv3(1:n,i)-K*v3(1:n,i));
elseif kic==4;
    v4(:,i)=vk(:);
vv4(1:n,i)=inv(M)*(vk(n+1:2*n)-0.5*C*vk(1:n));
aaii4(1:n,i)=inv(M)*(f-C*vv4(1:n,i)-K*v4(1:n,i));
elseif kic==5;
    v5(:,i)=vk(:);
vv5(1:n,i)=inv(M)*(vk(n+1:2*n)-0.5*C*vk(1:n));
aaii5(1:n,i)=inv(M)*(f-C*vv5(1:n,i)-K*v5(1:n,i));
end
for ll=1:(n-1)
         xiangduix(ll)=zhuanzhix(ll)-zhuanzhix(ll+1);
         xdnj(ll)=xiangduix(ll)*gangdu(ll+1);
end
if kic==0;
Txyz(2:n,i)=xdnj';
xyz(2:n,i)=xiangduix;
elseif kic==1;
 xdnj1=xdnj';  
 xiangduix1=xiangduix;
Txyz1(2:n,i)=xdnj1;
xyz1(2:n,i)=xiangduix1;
elseif kic==2;
xdnj2=xdnj';  
xiangduix2=xiangduix;
Txyz2(2:n,i)=xdnj2;
xyz2(2:n,i)=xiangduix2;
elseif kic==3;
xdnj3=xdnj';
xiangduix3=xiangduix;
Txyz3(2:n,i)=xdnj3;
xyz3(2:n,i)=xiangduix3;
elseif kic==4;
    xdnj4=xdnj';
xiangduix4=xiangduix;
Txyz4(2:n,i)=xdnj4;
xyz4(2:n,i)=xiangduix4;
elseif kic==5;
    xdnj5=xdnj';
xiangduix5=xiangduix;
Txyz5(2:n,i)=xdnj5;
xyz5(2:n,i)=xiangduix5;
end

Ffff1(:,i)=F;
if oo>=2*(k+1)*pi
        k=k+1;
    end



end;



if kic==0
save('bneiranji','v','xyz','vv','Txyz','aaii','t','nn','Ts','h');
elseif kic==1
save('bneiranji1','v1','xyz1','vv1','Txyz1','aaii1','t','nn','Ts','h');

elseif kic==2
save('bneiranji2','v2','xyz2','vv2','Txyz2','aaii2','t','nn','Ts','h');
elseif kic==3
save('bneiranji3','v3','xyz3','vv3','Txyz3','aaii3','t','nn','Ts','h');

elseif kic==4
save('bneiranji4','v4','xyz4','vv4','Txyz4','aaii4','t','nn','Ts','h');
elseif kic==5
save('bneiranji5','v5','xyz5','vv5','Txyz5','aaii5','t','nn','Ts','h');
else
    disp('nosave');
end

if exn==1
mxdjwy12=max(xyz1(2,0.0*nn+1:nn-0));
mxdjwy23=max(xyz1(3,0.0*nn+1:nn-0));
mxdjwy34=max(xyz1(4,0.0*nn+1:nn-0));
mxdjwy45=max(xyz1(5,0.0*nn+1:nn-0));
mxdjwy56=max(xyz1(6,0.0*nn+1:nn-0));
mxdjwy67=max(xyz1(7,0.0*nn+1:nn-0));
mxdjwy78=max(xyz1(8,0.0*nn+1:nn-0));
mxdjwy89=max(xyz1(9,0.0*nn+1:nn-0));
mxdjwy910=max(xyz1(10,0.0*nn+1:nn-0));

mixdjwy12=min(xyz1(2,0.0*nn+1:nn-0));
mixdjwy23=min(xyz1(3,0.0*nn+1:nn-0));
mixdjwy34=min(xyz1(4,0.0*nn+1:nn-0));
mixdjwy45=min(xyz1(5,0.0*nn+1:nn-0));
mixdjwy56=min(xyz1(6,0.0*nn+1:nn-0));
mixdjwy67=min(xyz1(7,0.0*nn+1:nn-0));
mixdjwy78=min(xyz1(8,0.0*nn+1:nn-0));
mixdjwy89=min(xyz1(9,0.0*nn+1:nn-0));
mixdjwy910=min(xyz1(10,0.0*nn+1:nn-0));




 
 
wn12=pi*(0.0001)^3/16;
wn23=pi*(0.21)^3/16;
wn34=pi*(0.21)^3/16;
wn45=pi*(0.21)^3/16;
wn56=pi*(0.21)^3/16;
wn67=pi*(0.21)^3/16;
wn78=pi*(0.21)^3/16;
wn89=pi*(0.3)^3/16;
wn910=pi*(0.22)^3/16;


xyz12=Txyz(2,0.0*nn+1:nn-0);
maxTxyz12=max(Txyz(7,0.0*nn+1:nn-0));
Tmax12=max(xyz12);
Tmin12=min(xyz12);
yingli12=xyz12/wn12/1e6;
maxyl12=max(yingli12);
minyl12=min(yingli12);
feng12=maxyl12-minyl12;
Tfeng12=Tmax12-Tmin12;

xyz23=Txyz(3,0.0*nn+1:nn-0);
maxTxyz23=max(Txyz(7,0.0*nn+1:nn-0));
Tmax23=max(xyz23);
Tmin23=min(xyz23);
yingli23=xyz23/wn23/1e6;
maxyl23=max(yingli23);
minyl23=min(yingli23);
feng23=maxyl23-minyl23;
Tfeng23=Tmax23-Tmin23;

xyz34=Txyz(4,0.0*nn+1:nn-0);
maxTxyz34=max(Txyz(7,0.0*nn+1:nn-0));
Tmax34=max(xyz34);
Tmin34=min(xyz34);
yingli34=xyz34/wn34/1e6;
maxyl34=max(yingli34);
minyl34=min(yingli34);
feng34=maxyl34-minyl34;
Tfeng34=Tmax34-Tmin34;

xyz45=Txyz(5,0.0*nn+1:nn-0);
maxTxyz45=max(Txyz(7,0.0*nn+1:nn-0));
Tmax45=max(xyz45);
Tmin45=min(xyz45);
yingli45=xyz45/wn45/1e6;
maxyl45=max(yingli45);
minyl45=min(yingli45);
feng45=maxyl45-minyl45;
Tfeng45=Tmax45-Tmin45;

xyz56=Txyz(6,0.0*nn+1:nn-0);
maxTxyz56=max(Txyz(7,0.0*nn+1:nn-0));
Tmax56=max(xyz56);
Tmin56=min(xyz56);
yingli56=xyz56/wn56/1e6;
maxyl56=max(yingli56);
minyl56=min(yingli56);
feng56=maxyl56-minyl56;
Tfeng56=Tmax56-Tmin56;

xyz67=Txyz(7,0.0*nn+1:nn-0);
maxTxyz67=max(Txyz(7,0.0*nn+1:nn-0));
Tmax67=max(xyz67);
Tmin67=min(xyz67);
yingli67=xyz67/wn67/1e6;
maxyl67=max(yingli67);
minyl67=min(yingli67);
feng67=maxyl67-minyl67;
Tfeng67=Tmax67-Tmin67;

%%%%%%%%%%
xyz78=Txyz(8,0.0*nn+1:nn-0);
maxTxyz78=max(Txyz(8,0.0*nn+1:nn-0));
Tmax78=max(xyz78);
Tmin78=min(xyz78);
yingli78=xyz78/wn78/1e6;
maxyl78=max(yingli78);
minyl78=min(yingli78);
feng78=maxyl78-minyl78;
Tfeng78=Tmax78-Tmin78;

xyz89=Txyz(9,0.0*nn+1:nn-0);
maxTxyz89=max(Txyz(9,0.0*nn+1:nn-0));
Tmax89=max(xyz89);
Tmin89=min(xyz89);
yingli89=xyz89/wn89/1e6;
maxyl89=max(yingli89);
minyl89=min(yingli89);
feng89=maxyl89-minyl89;
Tfeng89=Tmax89-Tmin89;

xyz910=Txyz(10,0.0*nn+1:nn-0);
maxTxyz910=max(Txyz(7,0.0*nn+1:nn-0));
Tmax910=max(xyz910);
Tmin910=min(xyz910);
yingli910=xyz910/wn910/1e6;
maxyl910=max(yingli910);
minyl910=min(yingli910);
feng910=maxyl910-minyl910;
Tfeng910=Tmax910-Tmin910;


%TTmm(:,1)=[Tmax12;Tmax23;Tmax34;Tmax45;Tmax56;Tmax67;Tmax78;Tmax89;Tmax910];
YLmax(:,1)=[0;maxyl23;maxyl34;maxyl45;maxyl56;maxyl67;maxyl78;maxyl89;maxyl910];
YLfeng(:,1)=[0;feng23;feng34;feng45;feng56;feng67;feng78;feng89;feng910];
Tfeng(:,1)=[0;Tfeng23;Tfeng34;Tfeng45;Tfeng56;Tfeng67;Tfeng78;Tfeng89;Tfeng910];
mTi=[0;maxTxyz23;maxTxyz34;maxTxyz45;maxTxyz56;maxTxyz67;maxTxyz78;maxTxyz89;maxTxyz910];
mxdi=[0;mxdjwy23;mxdjwy34;mxdjwy45;mxdjwy56;mxdjwy67;mxdjwy78;mxdjwy89;mxdjwy910];
mixdi=[0;mixdjwy23;mixdjwy34;mixdjwy45;mixdjwy56;mixdjwy67;mixdjwy78;mixdjwy89;mixdjwy910];


%%%%%%%%%%%%%%%%%%%%%%曲线%%%%%%%%%%%%%

if tu1==0
figure(1)
plot(t,yingli12,'g');title('第12轴段变化曲线');


figure(2)
plot(t,yingli23,'g');title('第23轴段变化曲线');


figure(3)
plot(t,yingli34,'g');title('第34轴段变化曲线');
figure(4)
plot(t,yingli45,'r');title('第45轴段变化曲线');
figure(5)
plot(t,yingli56,'g');title('第56轴段变化曲线');
figure(6)
plot(t,yingli67,'r');title('第67轴段变化曲线');
figure(7)
plot(t,yingli78,'g');title('第78轴段变化曲线');
figure(8)
plot(t,yingli89,'r');title('第89轴段变化曲线');
figure(9)
plot(t,yingli910,'r');title('第910轴段变化曲线');

elseif tu1==1
figure(1)
plot(t,xyz(2,:),'g');title('第12轴段变化曲线');


figure(2)
plot(t,xyz(3,:),'g');title('第23轴段变化曲线');


figure(3)
plot(t,xyz(4,:),'g');title('第34轴段变化曲线');
figure(4)
plot(t,xyz(5,:),'r');title('第45轴段变化曲线');
figure(5)
plot(t,xyz(6,:),'g');title('第56轴段变化曲线');
figure(6)
plot(t,xyz(7,:),'r');title('第67轴段变化曲线');
figure(7)
plot(t,xyz(8,:),'g');title('第78轴段变化曲线');
figure(8)
plot(t,xyz(9,:),'r');title('第89轴段变化曲线');
figure(9)
plot(t,xyz(10,:),'r');title('第910轴段变化曲线');
figure(10)
plot(t,xyz(11,:),'r');title('第1011轴段变化曲线');



elseif tu1==2
figure(1)
plot(t,Txyz(2,:),'g');title('第12轴段变化曲线');


figure(2)
plot(t,Txyz(3,:),'g');title('第23轴段变化曲线');


figure(3)
plot(t,Txyz(4,:),'g');title('第34轴段变化曲线');
figure(4)
plot(t,Txyz(5,:),'r');title('第45轴段变化曲线');
figure(5)
plot(t,Txyz(6,:),'g');title('第56轴段变化曲线');
figure(6)
plot(t,Txyz(7,:),'r');title('第67轴段变化曲线');
figure(7)
plot(t,Txyz(8,:),'g');title('第78轴段变化曲线');
figure(8)
plot(t,Txyz(9,:),'r');title('第89轴段变化曲线');
figure(9)
plot(t,Txyz(10,:),'r');title('第910轴段变化曲线');


elseif tu1==3
    figure(1)
plot(t,xyz1(2,:),'g');title('第12轴段变化曲线');


figure(2)
plot(t,xyz1(3,:),'g');title('第23轴段变化曲线');


figure(3)
plot(t,xyz1(4,:),'g');title('第34轴段变化曲线');
figure(4)
plot(t,xyz1(5,:),'r');title('第45轴段变化曲线');
figure(5)
plot(t,xyz1(6,:),'g');title('第56轴段变化曲线');
figure(6)
plot(t,xyz1(7,:),'r');title('第67轴段变化曲线');
figure(7)
plot(t,xyz1(8,:),'g');title('第78轴段变化曲线');
figure(8)
plot(t,xyz1(9,:),'r');title('第89轴段变化曲线');
figure(9)
plot(t,xyz1(10,:),'r');title('第910轴段变化曲线');


end
end