《内燃机课程设计程序》由会员分享，可在线阅读，更多相关《内燃机课程设计程序（7页珍藏版）》请在金锄头文库上搜索。
1、内燃机课程设计程序内燃机课程设计程序%(1)活塞加速度 J 随曲轴转角 的变化曲线(J-);cleara1=0:1:720;S=0.078r=S/2L=0.121n=100w=2*pi*nJ=r*w2*(cos(a1/180*pi)+r/L*cos(2*a1/180*pi)figure,plot(a1,J);grid on;title(活塞加速度 J 随曲轴转角 的变化曲线(J-);xlabel(a();ylabel(加速度 J(m/s2);%(2)往复质量惯性力 Pj 随曲轴转角 的变化曲线(Pj-);cleara2=0:1:720;m01=0.21;m1=0.09;m2=0.268;m02
2、=m1+m2;S=0.078;r=S/2;L=0.121;n=100;w=2*pi*nL1=m2*L/m02;L2=m1*L/m01;mj=m01+L2/L*m02;Fj=-mj*r*w2*(cos(a2/180*pi)+r/L*cos(2*a2/180*pi);figure;plot(a2,Fj);grid on;title(往复质量惯性力 Pj 随曲轴转角 的变化曲线(Pj-);xlabel(a();ylabel(往复质量惯性力 Pj(N);clearPa=0.98/10;n1=1.38;Vs=0.3/1000;Vc=Vs/8.3;n=100;D=0.070;S=0.078;r=S/2;L
3、=0.121;w=2*pi*n;m01=0.21;m1=0.09;m2=0.268;m02=m1+m2;L1=m2*L/m02;L2=m1*L/m01;mj=m01+L2/L*m02;%压缩a1=180:1:360;Bta=asin(r/L*sin(a1/180*pi);A1=1-cos(a1/180*pi);A2=1-cos(a1/180*pi.*2)A3=r/(4*L)A4=A2*A3A5=A1+A4A6=pi*D2/4*rVx1=A6*A5;V1=Vc+Vx1e1=(Vc+Vx1).(Vc+Vs);Pc=Pa*e1.n1;Pc1=Pa*9.31.38F11=Pc*pi*D2/4*106F
4、12=-mj*r*w2*(cos(a1/180*pi)+r/L*cos(2*a1/180*pi)F1=F11+F12figure,plot(a1,F1);FL1=F1./(cos(Bta)figure,plot(a1,FL1);Fn1=F1.*tan(Bta);figure,plot(a1,Fn1);Ft1=FL1.*sin(a1/180*pi+Bta);figure,plot(a1,Ft1);Fk1=FL1.*cos(a1/180*pi+Bta);figure,plot(a1,Fk1);%连线g=360:1:372;Btg=asin(r/L*sin(g/180*pi);G=785.9,100
5、0,1500,2000,3000,5000,10000,13000,15000,17800,17900,17900,17920figure,plot(g,G);FL5=G./(cos(Btg)figure,plot(g,FL5);Fn5=G.*tan(Btg);figure,plot(g,Fn5);Ft5=FL5.*sin(g/180*pi+Btg);figure,plot(g,Ft5);Fk5=FL5.*cos(g/180*pi+Btg);figure,plot(g,Fk5);%膨胀b1=372:1:540;Btb=asin(r/L*sin(b1/180*pi);B1=1-cos(b1/18
6、0*pi);B2=1-cos(b1/180*pi.*2);B3=r/(4*L);B4=B2*B3;B5=B1+B4;B6=pi*D2/4*r;Vx2=B6*B5;V2=Vc+Vx2;e2=(Vc).(Vc+Vx2);n2=1.24;zgp=3.52;%最高压力升高比;Pz=zgp*Pc1;Pb=(e2.n2).Pz;Pb1=Pz/9.31.24;F21=Pb*pi*D2/4*106F22=-mj*r*w2*(cos(b1/180*pi)+r/L*cos(2*b1/180*pi)F2=F21+F22;figure,plot(b1,F2);FL2=F2./(cos(Btb);figure,plot
7、(b1,FL2);Fn2=F2.*tan(Btb);figure,plot(b1,Fn2);Ft2=FL2.*sin(b1/180*pi+Btb);figure,plot(b1,Ft2);Fk2=FL2.*cos(b1/180*pi+Btb);figure,plot(b1,Fk2);%进气c1=0:1:180;Btc=asin(r/L*sin(c1/180*pi);C1=1-cos(c1/180*pi);C2=1-cos(c1/180*pi.*2)C3=r/(4*L)C4=C2*C3C5=C1+C4C6=pi*D2/4*rVx3=C6*C5V3=Vc+Vx3P0=0.09996;aa1=(0.
8、095-P0)/Vs;bb1=0.095-aa1*Vc;Py1=aa1.*V3+bb1;F31=Py1*pi*D2/4*106F32=-mj*r*w2*(cos(c1/180*pi)+r/L*cos(2*c1/180*pi)F3=F31+F32;figure,plot(c1,F3);FL3=F3./(cos(Btc)figure,plot(c1,FL3);Fn3=F3.*tan(Btc);figure,plot(c1,Fn3);Ft3=FL3.*sin(c1/180*pi+Btc);figure,plot(c1,Ft3);FL3=F3./(cos(Btc);Fk3=FL3.*cos(c1/18
9、0*pi+Btc);%排气d1=540:1:720;Btd=asin(r/L*sin(d1/180*pi);D1=1-cos(d1/180*pi);D2=1-cos(d1/180*pi.*2)D3=r/(4*L)D4=D2*D3D5=D1+D4D6=pi*D2/4*rVx4=D6*D5V4=Vc+Vx4;P0=0.09996;aa4=(Pb1-P0)/Vs;bb1=0.095-aa4*Vc;Py4=aa4.*V4+bb1;F41=Py4*pi*D2/4*106F42=-mj*r*w2*(cos(d1/180*pi)+r/L*cos(2*d1/180*pi)F4=F41+F42;figure,p
10、lot(d1,F4);FL4=F4./(cos(Btd)figure,plot(d1,FL4);Fn4=F4.*tan(Btd);figure,plot(d1,Fn4);Ft4=FL4.*sin(d1/180*pi+Btd);figure,plot(d1,Ft4);Fk4=FL4.*cos(d1/180*pi+Btd);figure,plot(d1,Fk4);figure,plot(a1,F1,b1,F2,g,G,c1,F3,d1,F4);grid on;title(作用在活塞上的合力 P 随曲轴转角 的变化曲线(P-);xlabel();ylabel(合力(N);figure;plot(a1
11、,FL1,b1,FL2,g,FL5,c1,FL3,d1,FL4);grid on;title(连杆力 Pl 随曲轴转角 的变化曲线（Pl-）);xlabel();ylabel(连杆力(N);figure,plot(a1,Fn1,b1,Fn2,g,Fn5,c1,Fn3,d1,Fn4);grid on;title(气缸侧压力 PN 随曲轴转角 的变化曲线（PN-）);xlabel();ylabel(气缸侧压力 PN(N);figure,plot(a1,Ft1,b1,Ft2,g,Ft5,c1,Ft3,d1,Ft4);grid on;title(作用在曲柄销上的切向力 T 随曲轴转角 的变化曲线（T-）);xlabel();ylabel(切向力 T(N);figure,plot(a1,Fk1,b1,Fk2,g,Fk5,c1,Fk3,d1,Fk4);grid on;title(作用在曲柄销上的径向力 K 随曲轴转角 的变化曲线（K-）);xlabel();ylabel(径向力 K(N);
