tuxing.m

内燃机转子仿真

load 'bneiranji';
load 'bneiranji1';
load 'bneiranji2';
load 'bneiranji3';
load 'bneiranji4';
load 'bneiranji5';
figure(1);
subplot(6,1,1);
plot(t*(0.6/h),xyz(2,:));
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
xlim([7200 7920]);
title('额定转速下自由端扭角位移时域图');
subplot(6,1,2);
plot(t*(0.6/h),xyz1(2,:));
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
xlim([7200 7920]);
title('额定转速e=0.1自由端扭角位移时域图');
subplot(6,1,3);
plot(t*(0.6/h),xyz2(2,:));
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
xlim([7200 7920]);
title('额定转速e=0.2自由端扭角位移时域图');
subplot(6,1,4);
plot(t*(0.6/h),xyz3(2,:));
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
xlim([7200 7920]);
title('额定转速e=0.3自由端扭角位移时域图');
subplot(6,1,5);
plot(t*(0.6/h),xyz4(2,:));
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
xlim([7200 7920]);
title('额定转速e=0.4自由端扭角位移时域图');
subplot(6,1,6);
plot(t*(0.6/h),xyz5(2,:));
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
xlim([7200 7920]);
title('额定转速e=0.45自由端扭角位移时域图');

figure(2);
subplot(6,1,1);
plot(t*(0.6/h),xyz(9,:));
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
xlim([7200 7920]);
title('额定转速下第九轴段扭角位移时域图');
subplot(6,1,2);
plot(t*(0.6/h),xyz1(9,:));
xlim([7200 7920]);
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
title('额定转速e=0.1第九轴段扭角位移时域图');
subplot(6,1,3);
plot(t*(0.6/h),xyz2(9,:));
xlim([7200 7920]);
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
title('额定转速e=0.2第九轴段扭角位移时域图');
subplot(6,1,4);
plot(t*(0.6/h),xyz3(9,:));
xlim([7200 7920]);
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
title('额定转速e=0.3第九轴段扭角位移时域图');
subplot(6,1,5);
plot(t*(0.6/h),xyz4(9,:));
xlim([7200 7920]);
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
title('额定转速e=0.4第九轴段扭角位移时域图');
subplot(6,1,6);
plot(t*(0.6/h),xyz5(9,:));
xlim([7200 7920]);
ylabel('扭角位移(rad)');
xlabel('曲轴转角(度)');
title('额定转速e=0.45第九轴段扭角位移时域图');