ttlwsj.m

内燃机转子仿真

clear all
 load 'pyjsdata'
 load 'datasj1'
 
 

maxylall=max(YLfengall) ;
 x=[400:50:1100];
 www=[400:25:1100];
 xxxx(1:length(www))=49;
 figure(1);
 subplot(3,1,1);
 plot(www,Tfengall(2,:)/2,'k',x,Mfuzhi(2,:),'-.k');
 title('第2轴段的最大力矩曲线');
 xlabel('转速(r.p.m)');
 ylabel('幅值(Nm)');
 h = legend('精细时程积分','文献',2);
  %subplot(7,1,2);
 %plot(www,Tfengall(3,:)/2,x,Mfuzhi(3,:),'-.');
  %subplot(7,1,3);
% plot(www,Tfengall(4,:)/2,x,Mfuzhi(4,:),'-.');
  %subplot(7,1,4);
% plot(www,Tfengall(5,:)/2,x,Mfuzhi(5,:),'-.');
  %title('第5轴段的最大力矩曲线');
 %xlabel('横轴转速(r.p.m)');
 %ylabel('幅值(Nm)');
% h = legend('精细时程积分','文献',2);
 % subplot(7,1,5);
% plot(www,Tfengall(6,:)/2,x,Mfuzhi(6,:),'-.');
 % title('第6轴段的最大力矩曲线');
 %xlabel('横轴转速(r.p.m)');
 %ylabel('幅值(Nm)');
 %h = legend('精细时程积分','文献',2);
 subplot(3,1,2);
 plot(www,Tfengall(7,:)/2,'k',x,Mfuzhi(7,:),'-.k');
   title('第7轴段的最大力矩曲线');
 xlabel('转速(r.p.m)');
 ylabel('幅值(Nm)');
 h = legend('精细时程积分','文献',2);
 % subplot(7,1,7);
 %plot(www,Tfengall(8,:)/2,x,Mfuzhi(8,:),'-.');
   subplot(3,1,3);
 plot(www,Tfengall(9,:)/2,'k',x,Mfuzhi(9,:),'-.k');
   title('第9轴段的最大力矩曲线');
 xlabel('转速(r.p.m)');
 ylabel('幅值(Nm)');
 h = legend('精细时程积分','文献',2);
 
 figure(2);
plot(www,maxylall/2,'k',x,Tfuzhi,'.-k',www,xxxx,'--k');
title('内燃机轴系最大应力');
 xlabel('转速(r.p.m)');
 ylabel('幅值(Mpa)');
  h = legend('精细时程积分','文献','许用应力',2);