《柴油机压缩膨胀过程matlab》由会员分享，可在线阅读，更多相关《柴油机压缩膨胀过程matlab（2页珍藏版）》请在金锄头文库上搜索。
1、Matlab 调用x,y=ode45(fun_compress,220.5,310.5,320);function dY=fun_compress(x,y)以下编辑成m函数T1=340+273 ;T2=290+273 ;T3=240+273;n=1860;S=0.195;R=287.11;w=n*6;x0=220.5;J=0.28;y0=305;%连杆曲柄比D=0.170;p0=2.5*105;ep=13.1;%压比A1=(D/2)2*pi;V0=pi*D2/4*(S/(ep-1)+S/2*(1+1/J)-cos(pi/180*x0)-1/J*sqrt(1-(J2)*(sin(pi/180*x
2、0)2);m=p0*V0/(y0*R);A2=(D/2)2*pi;Sc=S/(ep-1);%高度余隙V=pi*D2/4*(S/(ep-1)+S/2*(1+1/J)-cos(pi/180*x)-1/J*sqrt(1-(J2)*(sin(pi/180*x)2);B=pi2*D2*S/(8*180)*(sin(pi*x/180)+J/2*sin(pi*2*x/180)/sqrt(1-(J2)*(sin(pi/180*x)2);%气缸工作容积dedaoshuCm=n*S/30;%活塞平均速度p=m*R*y/V;ag=7.8*10(-3)*Cm(1/3)*sqrt(p*y);U=S/2*(1+1/J)-
3、cos(pi*x/180)-1/J*sqrt(1-J2*(sin(pi*x/180)2);A3=pi*D*(Sc+U);Qw=ag*(A1*(-y+T1)+A2*(-y+T2)+A3*(-y+T3)/w; %周壁传热 A1活塞，A2气缸盖cv=144.55*(-3)*0.0975*(10)(-6)*(y-273.15)(2)+2*7.768*(10)(-4)*(y-273.15)+489.6*(10)(-2);dY=(Qw-p*B)/(m*cv);function dY=fire1(x,y)T1=340+273 ;T2=290+273 ;T3=240+273;n=1860;S=0.195;R=
4、287.11;w=n*6;x0=220.5;J=0.28;y0=305;%连杆曲柄比D=0.170;p0=2.5*105;ep=13.1;%压比Cm=n*S/30;%活塞平均速度A1=(D/2)2*pi;V0=pi*D2/4*(S/(ep-1)+S/2*(1+1/J)-cos(pi/180*x0)-1/J*sqrt(1-(J2)*(sin(pi/180*x0)2);m=p0*V0/(y0*R);A2=(D/2)2*pi;U=S/2*(1+1/J)-cos(pi*x/180)-1/J*sqrt(1-J2*(sin(pi*x/180)2);Sc=S/(ep-1);%高度余隙A3=pi*D*(Sc+
5、U);V=pi*(0.17)2/4*(0.195)/12.5+(0.195)/2*(1+1/0.28)-cos(pi/180*x)-1/(0.28)*(1-(0.28)2*(sin(pi/180*x)2)(1/2);%瞬时气缸容积p=m*R*y/V;ag=7.8*10(-3)*Cm(1/3)*sqrt(p*y);%Eicheiberg 瞬时平均换热系数Qw=ag*(A1*(-y+T1)+A2*(-y+T2)+A3*(-y+T3)/w;%气缸周壁传热mx=12.5*10(-3)+0.51434*10(-3)*(1-(exp(-6.908*(x-342)/70)4)/100;%理论空气数14.3
6、总质量=初始空气质量+气缸内以燃烧的质量a=12.5*10(-3)/(14.3*0.51434*10(-3)*(1-(exp(-6.908*(x-352)/70)4)/100);%燃烧起始角352;Cv=144.55*(-3)*(0.0975+0.0485/a0.75)*(10)(-6)*(y-273.15)(2)+2*(7.768+3.36/a0.8)*(10)(-4)*(y-273.15)+(489.6+46.4/a0.03)*(10)(-2);DV=pi2*D2*S/(8*180)*(sin(pi*x/180)+J/2*sin(pi*2*x/180)/sqrt(1-(J2)*(sin(p
7、i/180*x)2);%气缸容积随曲周转角的变化率Dx=6.908*4/70*(x-342)/70)(3)*exp(-6.908*(x-342)/70)4);% 燃烧品质指术3 持续70%缸内燃料的瞬时放热率 循环喷油量0.51434*10(-3) Hu=42692*103u=0.14455*(-(0.0975+0.0485/(a(0.75)*(y-273.15)3*10(-6)+(7.768+3.36/(a(0.8)*(y-273.15)2*10(-4)+(489.6+46.4/(a(0.03)*(y-273.15)*10(-2)+1356.8);dY=1/(mx*Cv)*(Qw-p*DV+
8、0.51434*10(-3)*(43000*103-u)*Dx);function dY=PZ(x,y)T1=340+273 ;T2=290+273 ;T3=240+273;n=1860;S=0.195;R=287.11;w=n*6;x0=220.5;J=0.28;y0=305;%连杆曲柄比D=0.170;p0=2.5*105;ep=13.1;%压比A1=(D/2)2*pi;V0=pi*D2/4*(S/(ep-1)+S/2*(1+1/J)-cos(pi/180*x0)-1/J*sqrt(1-(J2)*(sin(pi/180*x0)2);m=p0*V0/(y0*R);mx=m+0.51434*1
9、0(-3)*(1-(exp(-6.908*(422-352)/70)4);%理论空气数14.3 总质量=初始空气质量+气缸内以燃烧的质量a=12.5*10(-3)/(14.3*mx);Cv=144.55*(-3)*(0.0975+0.0485/a0.75)*(10)(-6)*(y-273.15)(2)+2*(7.768+3.36/a0.8)*(10)(-4)*(y-273.15)+(489.6+46.4/a0.03)*(10)(-2);A2=(D/2)2*pi;Sc=S/(ep-1);%高度余隙V=pi*D2/4*(S/(ep-1)+S/2*(1+1/J)-cos(pi/180*x)-1/J*
10、sqrt(1-(J2)*(sin(pi/180*x)2);B=pi2*D2*S/(8*180)*(sin(pi*x/180)+J/2*sin(pi*2*x/180)/sqrt(1-(J2)*(sin(pi/180*x)2);%气缸工作容积dedaoshuCm=n*S/30;%活塞平均速度p=m*R*y/V;ag=7.8*10(-3)*Cm(1/3)*sqrt(p*y);U=S/2*(1+1/J)-cos(pi*x/180)-1/J*sqrt(1-J2*(sin(pi*x/180)2);A3=pi*D*(Sc+U);Qw=ag*(A1*(-y+T1)+A2*(-y+T2)+A3*(-y+T3)/w; %周壁传热 A1活塞，A2气缸盖dY=(Qw-p*B)/(mx*Cv);
