热力学计算.docx

4100型非增压柴油机热力学计算源程序#includc graphics.h#include conio.h#include dos.h#include stdio.h#include math.hvoid draw(nl,Pa,s 1 ,n2,Pb,l,p,Pz)float n 1 ,Pa,s 1 ,n2,Pb,l,p,Pz;(float k,j;int gdriver,gniode,i,x,y,r,L=500,B=50,u,z;gdriver=DETECT;initgraph(&gdriver,&gmode,);setbkcolor( 15);setcolor(8);cleardevice();line(20,55,20,410);line(20,410,600,410);line(2O.55,l 7,65);line(20,55,23,65);line(600,410,590,407);linc(600,410,590,413);gotoxy(2.4);printf(p);gotoxy(2,5);printf(/);gotoxy(2,6);printf(Ma);k=(float)(L/s l);u=L;z=(int)(B*Pa);gotoxy( 1,2);printf(rn 叩 1-13125);x=L;r=(int)(B*Pa*pow(float)L/x,n 1);u=x;z=r;for(x=L;x=(int)(k);x-)r=(int)(B*Pa*pow(float)L/x,n 1);line(u,400-z,x,400-r);u=x,z=r;line(x+1,400-r,(int)(p*L/s 1 -1 ),(int)(400-B *Pz);j=(float)(L/l);u=L;z=(int)(B*Pb);moveto(505,415);outtext(V* IOe-4/(m*m*m);y=L;r=(int)(B *Pb*pow(float)L/y,n2);u=y;z=r;for(y=L;y=(int)j;y-)r=(int)(B*Pb*pow(float)L/y,n2);line(u,400-z,y.400-r);u=y;z=r;gotoxy(6,4)；prinlf(z,);moveto(505,375);outtext(nb);movcto(25,250);outtcxt(c);mo veto(505,400) ;ou ttext( a);arc(L,400-(int)(B*(Pa+(Pb-Pa)/2),270,450,(int)(B*(Pb-Pa)/2);getch();closegraph();main() float Hu=42705,gc=0.86,gh=0.13,go=0.01;float Pz=6.0,a= 1.6,sz=0.70,sb=0.86,Rr=0.03,Tr=750,Qin=0.84,s= 16,P0=0.09,T0=293;float Tl = l 0,sc= 1.11 ,a2= 19.26,b=0.0025,o=0.98;floatTa,Pa,Qv,n 1 ,n2,m,n,Tc,Pc,L(),L,gi,B,B0,Xz,Bz,Cvpmx,Cppmx,Cvpmc,Cvpmz,Cvpmb,Cv l,Cv 11,Cv2,Va,Tz,Tb,Tk,k,l,Pb,Pk,Pi 1 ,Pi,Qi,Pmi,ge,Qe,Pe,Cm,Hu 1 ,Cppmz,p,Bb;printf(”*n”)Ta=(T0+Tl +sc*Rr*Tr)/( 1 +Rr);Pa=0.87*P0;Qv=s*T0*Pa/(s-1 )*Ta*P0*( 1 +Rr);printf(Ta=%fKn,Ta);printf(Qv=%fn.Qv);n 1 = 1.32;do m=nl;nl=l+8.315/(a2+b*Ta*(l +pow(s,(m-1); while(fabs(m-nl )0.0001);piintf(nl=%fnnl);Tc=Ta*pow(s,n 1 -1);Pc=Pa*pow(s,nl);L0=(gc/12+gh/4-go/32)/0.21;L=a*L0;B0=l+0.0639/a;B=(B0+Rr)/(l+Rr);Xz=sz/sb;Bz=l+(B0-l)*Xz/(l+Rr);printf(,Tc=%fKnPc=%fMPan,Tc,Pc);Cv2=19.26+0.0025*Tc;Cvll=20.47+0.0036*Tc;Cvpmc=(l .064+Rr)*Cv 11 +(a-1)*(1 +Rr)*Cv2)/(a*( 1 +Rr)+0.064); Tb=293;Cvl = l 9.26+0.0()25*Tb;Cv 11 =20.47+0.0036*Tb;Cvpmb=(1.064+Rr)*Cv 11 +(a-1)*(1 +Rr)*Cv l)/(a*(l +Rr)+0.064);Hui =Hu+L*( 1 +Rr)*(B *C vpmb-C v 1)*293;printfCHu l=%fkJ/kgn,Hu 1);k=Pz/Pc;Tb=897;Cvl = 19.26+0.0025 *Tb:Cv 11 =20.47+0.0036叮 b;Cvpmb=(L064+Rr)*Cv II +(a-1)*(1 +Rr)*Cv l)/(a*(l +Rr)+0.064);n=1800;do(Tz=n;Cvl = 19.26+0.0025 叮 z;Cv 11 =20.47+0.0036*Tz;Cvpmz=( 1.064*Xz+Rr)*Cv 11 +(a*(l +Rr)-(Xz+Rr)*Cv l)/(a*(l +Rr)+0.064*Xz);Cppmz=Cvpmz+8.315;n=(sz*Hu l/(a*L0)+(Cv2+8.315*k)*Tc+Rr*(Cvpmc+8.315*k)*Tc)/(Bz*( 1 +Rr)*Cppmz); while(fabs(Tz-n)0.0001);printf(,Tz=%fKn,Tz);p=B*Tz/(k*Tc);l=s/p;Tb=1000;do m=Tb;Cvl = l 9.26+0.0025 *Tb;Cv 11 =20.47+0.0036*Tb;n2=8.315*(Tz-Tb)/(sb-sz)*Hu 1/( l+Rr)*L*B)+Tz*Cvpmz-Cvpmb*Tb)+1;Tb=Tz*B/(Bz*pow(I,n2-1);while(fabs(m-Tb)0.01);printf(Tb=%fKnn2=%An,Tb,n2);Pb=Pz/pow(l,n2);printf(Pb=%fMPan.Pb);Pil =Pc*(k*(p- l)+k*p*(l- l/pow(l,n2- l)/(n2- l)-(l-l/pow(s,nl- l)/(n 1-l)/(s-1);Pi=Pil*o;gi=l 2.54e6*Qv*P0/(28.96*L*T0*Pi);Qi=3600.0*l 000.0/(gi*Hu);Pc=Pi*Qm;gc=gi/Qm;Qc=Qi*Qm;printfCPil =%fMPanPi=%fMPangi=%fg/kw.hnQi=%fnH,Pil,Pi,gi,Qi);pnntf(Pe=%fMPan,.Pe);printf(ge=%fg/kw.hnQe=%fn,ge,Qe); printf(”*#*n”) gctch();draw(n 1 ,Pa,s.n2,Pb,l,p,Pz,Pc);