neiranjinmk1.m

内燃机转子仿真

XDNJ78=XDNJ(7,:);
yingli78=XDNJ78/wn78/1e6;
maxyl78=max(yingli78);
minyl78=min(yingli78);
feng78=maxyl78-minyl78;

XDNJ89=XDNJ(8,:);
yingli89=XDNJ89/wn89/1e6;
maxyl89=max(yingli89);
minyl89=min(yingli89);
feng89=maxyl89-minyl89;

XDNJ910=XDNJ(9,:);
yingli910=XDNJ910/wn910/1e6;
maxyl910=max(yingli910);
minyl910=min(yingli910);
feng910=maxyl910-minyl910;

XDNJ1011=XDNJ(10,:);
yingli1011=XDNJ1011/wn1011/1e6;
maxyl1011=max(yingli1011);

XDNJ1112=XDNJ(11,:);
yingli1112=XDNJ1112/wn1112/1e6;
maxyl1112=max(yingli1112);


XDNJ1213=XDNJ(12,:);
yingli1213=XDNJ1213/wn1213/1e6;
maxyl1213=max(yingli1213);
minyl1213=min(yingli1213);
feng1213=maxyl1213-minyl1213;

XDNJ1314=XDNJ(13,:);
yingli1314=XDNJ1314/wn1314/1e6;
maxyl1314=max(yingli1314);
minyl1314=min(yingli1314);
feng1314=maxyl1314-minyl1314;

XDNJ1415=XDNJ(14,:);
yingli1415=XDNJ1415/wn1415/1e6;
maxyl1415=max(yingli1415)
minyl1415=min(yingli1415);
feng1415=maxyl1415-minyl1415;;

XDNJ1516=XDNJ(15,:);
yingli1516=XDNJ1516/wn1516/1e6;
maxyl1516=max(yingli1516);
minyl1516=min(yingli1516);
feng1516=maxyl1516-minyl1516;

%%%%%%%%%%%%%%%%%%%%%%曲线%%%%%%%%%%%%%
wn12=pi*(0.21)^3/16;
wn23=pi*(0.21)^3/16;
wn34=pi*(0.23)^3/16;
wn45=pi*(0.23)^3/16;
wn56=pi*(0.23)^3/16;
wn67=pi*(0.23)^3/16;
wn78=pi*(0.23)^3/16;
wn89=pi*(0.23)^3/16;
wn910=pi*(0.23)^3/16;
wn1011=pi*(0.23)^3/16;%联轴器内外之间的连接轴
wn1112=pi*(0.23)^3/16;%离合器主从之间的连接轴
wn1213=pi*(0.314^4-0.18^4)/0.314/16;%一号传动轴
wn1314=pi*(0.314^4-0.18^4)/0.314/16;%2号传动轴
wn1415=pi*(0.25^4-0.15^4)/0.25/16;%电轴1
wn1516=pi*(0.38^4-0.064^4)/0.38/16;%两个电枢之间

xyz12=xyz(2,:);
yingli12=xyz12/wn12/1e6;
maxyl12=max(yingli12);
minyl12=min(yingli12);
feng12=maxyl12-minyl12;

xyz23=xyz(3,:);
yingli23=xyz23/wn23/1e6;
maxyl23=max(yingli23);
minyl23=min(yingli23);
feng23=maxyl23-minyl23;

xyz34=xyz(4,:);
yingli34=xyz34/wn34/1e6;
maxyl34=max(yingli34);
minyl34=min(yingli34);
feng34=maxyl34-minyl34;

xyz45=xyz(5,:);
yingli45=xyz45/wn45/1e6;
maxyl45=max(yingli45);
minyl45=min(yingli45);
feng45=maxyl45-minyl45;

xyz56=xyz(6,:);
yingli56=xyz56/wn56/1e6;
maxyl56=max(yingli56);
minyl56=min(yingli56);
feng56=maxyl56-minyl56;

xyz67=xyz(7,:);
yingli67=xyz67/wn67/1e6;
maxyl67=max(yingli67);
minyl67=min(yingli67);
feng67=maxyl67-minyl67;

xyz78=xyz(8,:);
yingli78=xyz78/wn78/1e6;
maxyl78=max(yingli78);
minyl78=min(yingli78);
feng78=maxyl78-minyl78;

xyz89=xyz(9,:);
yingli89=xyz89/wn89/1e6;
maxyl89=max(yingli89);
minyl89=min(yingli89);
feng89=maxyl89-minyl89;

xyz910=xyz(10,:);
yingli910=xyz910/wn910/1e6;
maxyl910=max(yingli910);
minyl910=min(yingli910);
feng910=maxyl910-minyl910;

xyz1011=xyz(11,:);
yingli1011=xyz1011/wn1011/1e6;
maxyl1011=max(yingli1011);

xyz1112=xyz(12,:);
yingli1112=xyz1112/wn1112/1e6;
maxyl1112=max(yingli1112);


xyz1213=xyz(13,:);
yingli1213=xyz1213/wn1213/1e6;
maxyl1213=max(yingli1213);
minyl1213=min(yingli1213);
feng1213=maxyl1213-minyl1213;

xyz1314=xyz(14,:);
yingli1314=xyz1314/wn1314/1e6;
maxyl1314=max(yingli1314);
minyl1314=min(yingli1314);
feng1314=maxyl1314-minyl1314;

xyz1415=xyz(15,:);
yingli1415=xyz1415/wn1415/1e6;
maxyl1415=max(yingli1415);
minyl1415=min(yingli1415);
feng1415=maxyl1415-minyl1415;;

xyz1516=xyz(16,:);
yingli1516=xyz1516/wn1516/1e6;
maxyl1516=max(yingli1516);
minyl1516=min(yingli1516);
feng1516=maxyl1516-minyl1516;

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

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


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


figure(3)
plot(d*i,yingli34,'g');title('第34轴段变化曲线');
figure(4)
plot(d*i,yingli45,'r');title('第45轴段变化曲线');
figure(5)
plot(d*i,yingli56,'g');title('第56轴段变化曲线');
figure(6)
plot(d*i,yingli67,'r');title('第67轴段变化曲线');
figure(7)
plot(d*i,yingli78,'g');title('第78轴段变化曲线');
figure(8)
plot(d*i,yingli89,'r');title('第89轴段变化曲线');
figure(9)
plot(d*i,yingli910,'r');title('第910轴段变化曲线');
figure(10)
plot(d*i,yingli1011,'r');title('第1011轴段变化曲线');
figure(11)
plot(d*i,yingli1112,'r');title('第1112轴段变化曲线');

figure(12)
plot(d*i,yingli1213,'g');title('第1213轴段变化曲线');
figure(13)
plot(d*i,yingli1314,'r');title('第1314轴段变化曲线');
figure(14)
plot(d*i,yingli1415,'r');title('第1415轴段变化曲线');
figure(15)
plot(d*i,yingli1516,'g');title('第1516轴段变化曲线');

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轴段变化曲线');
figure(11)
plot(t,xyz(12,:),'r');title('第1112轴段变化曲线');    
figure(12)
plot(t,xyz(13,:),'g');title('第1213轴段变化曲线');
figure(13)
plot(t,xyz(14,:),'r');title('第1314轴段变化曲线');
figure(14)
plot(t,xyz(15,:),'r');title('第1415轴段变化曲线');
figure(15)
plot(t,xyz(16,:),'g');title('第1516轴段变化曲线');
elseif tu1==2
figure(1)
plot(t,xy(1,:),'g');title('第12轴段变化曲线');


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


figure(3)
plot(3*t/d,xy(3,:),'g');title('第34轴段变化曲线');
figure(4)
plot(3*t/d,xy(4,:),'r');title('第45轴段变化曲线');
figure(5)
plot(3*t/d,xy(5,:),'g');title('第56轴段变化曲线');
figure(6)
plot(3*t/d,xy(6,:),'r');title('第67轴段变化曲线');
figure(7)
plot(3*t/d,xy(7,:),'g');title('第78轴段变化曲线');
figure(8)
plot(3*t/d,xy(8,:),'r');title('第89轴段变化曲线');
figure(9)
plot(3*t/d,xy(9,:),'r');title('第910轴段变化曲线');
figure(10)
plot(3*t/d,xy(10,:),'r');title('第1011轴段变化曲线');
figure(11)
plot(3*t/d,xy(11,:),'r');title('第1112轴段变化曲线'); 
figure(12)
plot(3*t/d,xy(12,:),'r');title('第1112轴段变化曲线'); 
figure(13)
plot(3*t/d,xy(13,:),'r');title('第1112轴段变化曲线'); 
figure(14)
plot(3*t/d,xy(14,:),'r');title('第1112轴段变化曲线'); 
figure(15)
plot(3*t/d,xy(15,:),'r');title('第1112轴段变化曲线'); 
figure(16)
plot(3*t/d,xy(16,:),'r');title('第1112轴段变化曲线'); 
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轴段变化曲线');
figure(10)
plot(t,xyz1(11,:),'r');title('第1011轴段变化曲线');
figure(11)
plot(t,xyz1(12,:),'r');title('第1112轴段变化曲线');  
figure(12)
plot(t,xyz1(13,:),'r');title('第1112轴段变化曲线'); 
figure(13)
plot(t,xyz1(14,:),'r');title('第1112轴段变化曲线'); 
figure(14)
plot(t,xyz1(15,:),'r');title('第1112轴段变化曲线'); 

end