離心式壓縮機(jī)的熱力計(jì)算(C++)

1、離心式壓縮機(jī)的熱力計(jì)算 #include #include void main()double p0=200000,T0=350,ps_=380000,n=60000,G=0.25,H=0.635,c2r,c1=c2r=90,o1=0.98,a=0.05,n2=1.55,n3=1.9,n5=1.8;/壓頭系數(shù)取H=0.635,葉輪進(jìn)口空氣速度c1=90m/s,總壓恢復(fù)系數(shù)o1=0.98,/ D1h/D2=0.3 double sc=ps_/p0; /1.壓縮機(jī)增壓比 cout1.壓縮機(jī)增壓比sc=scendl;double lcs=1004.5*T0*(pow(sc,0.4/1.4)-1);

2、/2.等熵過程的壓縮功 cout2.等熵過程的壓縮功lcs=lcsendl; double u2=sqrt(lcs/H); /3,H=0.635葉輪外徑圓周速度coutu2=u2endl;double T1=T0-c1*c1/2/1005; /4.c1=90,葉輪進(jìn)口處溫度 coutT1=T1endl;double p1=o1*p0*(pow(T1/T0,3.5); /5進(jìn)氣管總壓恢復(fù)系數(shù)為0.98，葉輪進(jìn)口空氣壓力coutp1=p1endl;double j1=p1/T1/287; /6.葉輪進(jìn)口密度coutj1=j1endl; double D2=60*u2*1000/3.14159/n;

3、 /7.葉輪外徑coutD2=D2endl;double D1h=0.3*D2; /8,D1h/D2=0.3葉輪進(jìn)口葉根直徑coutD1h=D1hendl;double D1t=sqrt(D1h*D1h+4*G*1000000/3.1415/j1/c1); /9.葉輪進(jìn)口葉尖直徑coutD1t=D1tendl; coutD1t/D2=D1t/D2 一般范圍在0.500.65之間endl; double D1=sqrt(D1h*D1h/2+D1t*D1t/2) ; /10.葉輪進(jìn)口平均直徑coutD1=D1endl;double u1t=D1t/D2*u2; /11.葉輪進(jìn)口葉尖圓周速度cout

4、u1t=u1tendl; double w1t=sqrt(c1*c1+u1t*u1t); /12. 葉輪進(jìn)口葉尖處氣流相對(duì)速度coutw1t=w1tendl; double Mw1t=w1t/20.1/sqrt(T1); /13.葉輪進(jìn)口葉尖馬赫數(shù)coutMw1t=Mw1t 一般希望Mw1t取值范圍0.800.92endl;double u1=D1/D2*u2; /14.葉輪進(jìn)口平均直徑處的圓周速度coutu1=u1endl;double w1=sqrt(c1*c1+u1*u1); /15. .葉輪進(jìn)口平均直徑處的相對(duì)速度coutw1=w1endl;double zi=10+0.03*D2;

5、/葉輪葉片數(shù)coutzi=ziendl;double h=1/(1+2*3.1415/3/zi*(1/(1-D1*D1/D2/D2); /16.葉輪滑動(dòng)系數(shù)couth=hendl; double T2=T1+(h+a-h*h/2)*u2*u2/1005; /17.葉輪出口溫度coutT2=T2endl;double p2=p1*pow(T2/T1,n2/(n2-1); /18.葉輪出口壓力 coutp2=p2endl;double j2=p2/287/T2; /19.葉輪出口空氣密度coutj2=j2endl;double c2u=h*u2; /20.葉輪出口周向分速度coutc2u=c2ue

6、ndl;double c2=sqrt(c2u*c2u+c2r*c2r); /21.葉輪出口氣流絕對(duì)速度coutc2=c2endl;double a2=atan(c2r/c2u); /22.葉輪出口氣流角couta2=a2endl;double w2=sqrt(c2*c2+u2*u2-2*c2*u2*cos(a2); /23.葉輪出口相對(duì)速度coutw2=w2endl;coutw2/w1=w2/w1判斷葉輪中的相對(duì)運(yùn)動(dòng)是否是減速的endl; double l2=G*1000000/3.1415/D2/j2/c2r; /24.葉輪出口軸向?qū)挾萩outl2=l2endl;double D3=1.08

7、*D2; /25.無葉擴(kuò)壓器出口直徑coutD3=D3endl;double c3=D2/D3*c2; /26. 無葉擴(kuò)壓器出口出口氣流速度coutc3=c3endl;double T3=T2+(c2*c2-c3*c3)/2/1005; /27. 無葉擴(kuò)壓器出口溫度coutT3=T3endl; double p3=p2*pow(T3/T2,n3/(n3-1); /28.取n3=1.9, 無葉擴(kuò)壓器出口壓力coutp3=p3endl;double Mc3=c3/sqrt(T3)/20.1; /29.驗(yàn)算coutMc3=Mc3 是否=0.800.95endl;double a3g=3.1415/1

8、80*25; double a3=a3g;double a4g=a3g+3.1415/180*20;double a4=a4g;double D4=D3*cos(a3)/cos(a4); /30.直葉片擴(kuò)壓器出口直徑coutD4=D4endl;coutD4/D3=D4/D3 是否=1.21.4endl;double r4=D4/2;double r3=D3/2;double L=(r4-r3)/sin(a3g+a4g)/2); /31.擴(kuò)壓器流道直線長度coutL=Lendl;double aav=3.14159/18; /平均擴(kuò)張角取10coutaav=aavendl;double zd=3

9、.1415*(r4*sin(a4g)-r3*sin(a3g)/tan(aav/2)/L; /32.擴(kuò)壓器葉片數(shù)coutzd=zdendl;double n4=1.6,c4=111.3; /葉片擴(kuò)壓器多變指數(shù)n4=1.6,葉片擴(kuò)壓器出口速度c4=111.3/33.直葉片擴(kuò)壓器出口溫度double T4_=T3*pow(c3*D3*sin(a3)/c4/D4/sin(a4),(n4-1); coutT4_=T4_endl;double T4=T3+(c3*c3-c4*c4)/2/1005;coutT4=T4endl;coutT4-T4_的絕對(duì)值為fabs(T4-T4_)是否小于0.5endl; d

10、ouble p4=p3*pow(T4/T3,n4/(n4-1); /34. 葉片擴(kuò)壓器出口壓力coutp4=p4endl;double j4=p4/287/T4; /35. 葉片擴(kuò)壓器出口空氣密度coutj4=j4endl;double c4u=c4*cos(a4); /36. 葉片擴(kuò)壓器出口周向分速coutc4u=c4uendl;double K=c4u*r4/1000; /37.系數(shù)KcoutK=Kendl;double C=720*j4*K*3.1415/G; /38.系數(shù)CcoutC=Cendl;double jcmax=360*1000/C+sqrt(720*1000*r4/C);

11、/39.蝸殼出口圓截面直徑coutjcmax=jcmaxendl;double c5=G*1000000/3.1415/jcmax/jcmax/j4; /40.蝸殼出口氣流速度coutc5=c5endl;double T5=T4+(c4*c4-c5*c5)/2/1005; /41. 蝸殼出口空氣溫度coutT5=T5endl;double p5=p4*pow(T5/T4,n5/(n5-1); /42. 蝸殼出口壓力coutp5=p5endl;double Ts=T5+c5*c5/2010; /43. 蝸殼出口氣流總溫coutTs=Tsendl;double ps=p5*pow(Ts/T5,3.

12、5);/44 /44. 蝸殼出口氣流總壓coutps=psendl;cout與所要求的出口總壓的相對(duì)誤差fabs(ps-ps_)/ps_小于1%，滿足要求endl;double lc=1005*(Ts-T0); /45.壓縮機(jī)的實(shí)際耗功coutlc=lcendl;double Nc=G*lc/1000; /46.壓縮機(jī)消耗功率coutNc=Ncendl;double ncs=lcs/lc; /47.壓縮機(jī)絕熱效率cout壓縮機(jī)絕熱效率為ncsendl;結(jié)果：1.壓縮機(jī)增壓比sc=1.92.等熵過程的壓縮功lcs=70764.73.葉輪外徑圓周速度u2=333.8274.葉輪進(jìn)口處溫度T1=34

13、5.975.進(jìn)氣管總壓恢復(fù)系數(shù)為0.98，葉輪進(jìn)口空氣壓力p1=1882156.葉輪進(jìn)口密度j1=1.895547.葉輪外徑D2=106.2618.葉輪進(jìn)口葉根直徑D1h=31.87829.葉輪進(jìn)口葉尖直徑D1t=53.6853D1t/D2=0.505224 一般范圍在0.500.65之間10.葉輪進(jìn)口平均直徑D1=44.149411.葉輪進(jìn)口葉尖圓周速度u1t=168.65712. 葉輪進(jìn)口葉尖處氣流相對(duì)速度w1t=191.16813.葉輪進(jìn)口葉尖馬赫數(shù)Mw1t=0.511329 一般希望Mw1t取值范圍0.800.9214.葉輪進(jìn)口平均直徑處的圓周速度u1=138.69915.葉輪進(jìn)口平均

14、直徑處的相對(duì)速度w1=165.34葉輪葉片數(shù)zi=13.187816.葉輪滑動(dòng)系數(shù)h=0.83896717.葉輪出口溫度T2=405.5218.葉輪出口壓力p2=29446419.葉輪出口空氣密度j2=2.5301120.葉輪出口周向分速度c2u=280.0721.葉輪出口氣流絕對(duì)速度c2=294.17522.葉輪出口氣流角a2=0.31092623.葉輪出口相對(duì)速度w2=104.832w2/w1=0.63404判斷葉輪中的相對(duì)運(yùn)動(dòng)是否是減速的 /小于1為減速24.葉輪出口軸向?qū)挾萳2=3.2888925.無葉擴(kuò)壓器出口直徑D3=116.88726. 無葉擴(kuò)壓器出口出口氣流速度c3=267.4

15、3227. 無葉擴(kuò)壓器出口溫度T3=412.99228.無葉擴(kuò)壓器出口壓力p3=30603629.驗(yàn)算Mc3=0.654706 是否=0.800.9530.直葉片擴(kuò)壓器出口直徑D4=149.812D4/D3=1.28169 是否=1.21.431.擴(kuò)壓器流道直線長度L=28.703平均擴(kuò)張角取10aav=0.17453332.擴(kuò)壓器葉片數(shù)zd=35.362133.直葉片擴(kuò)壓器出口溫度T4_=442.16T4=442.411T4-T4_的絕對(duì)值為0.251063是否小于0.534. 葉片擴(kuò)壓器出口壓力p4=36767535. 葉片擴(kuò)壓器出口空氣密度j4=2.8957236. 葉片擴(kuò)壓器出口周向分

16、速c4u=78.702837.系數(shù)K=5.8953338.系數(shù)C=15445239.蝸殼出口圓截面直徑j(luò)cmax=21.017340.蝸殼出口氣流速度c5=62.214541. 蝸殼出口空氣溫度T5=446.64842. 蝸殼出口壓力p5=37564643.蝸殼出口氣流總溫Ts=448.57444. 蝸殼出口氣流總壓ps=381345與所要求的出口總壓的相對(duì)誤差0.00353975小于1%，滿足要求45.壓縮機(jī)的實(shí)際耗功lc=99066.846.壓縮機(jī)消耗功率Nc=24.7667壓縮機(jī)絕熱效率為0.714312Press any key to continue向心渦輪級(jí)的熱力計(jì)算#includ

17、e #include void main()double p0=360000, p3=110000 ,T0=360, G=0.25,D1=0.12,x=0.97, y=0.88, a1g=0.314159,B2=0.610865,u=0.4,zn=27;/預(yù)取D1=120mm,噴嘴環(huán)速度系數(shù)x=0.97,葉輪速度系數(shù)y=0.88/噴嘴環(huán)出口葉片幾何角a1g=18，葉輪出口相對(duì)氣流角B2=35,u=0.4 double a1=a1g; /double zi=0.03*(33-a1/3.14159*180)*(33-a1/3.14159*180)+12; /double modf(double z

18、i_,double *zi); double zi=18; cout 1.葉輪葉片數(shù)zi=ziendl; double M=1-2/zi;cout2.速度比M=Mendl;double i1=cos(a1)*u/M*cos(a1)*u/M; double i1_=cos(a1)*u/M*u/M; double i2=1/y/cos(B2)*1/y/cos(B2)-1; double i3=2*cos(a1)*cos(a1)/M-1; double pt=x*x*(i1*i2+i3)/(1+x*x*(i1_*i2+i3); cout3.反力度pt=ptendl;double x0=x*sqrt(

19、1-pt)*cos(a1)/M;cout4.渦輪速比X0=x0endl;double s=p0/p3;cout5.渦輪膨脹比s=sendl;double vhs=1005*T0*(1-1/(pow(s,0.286);cout6.渦輪總等焓降vhs=vhsendl;double Cs=sqrt(2*vhs);cout7.等熵速度Cs=Csendl;double C1=x*sqrt(1-pt)*Cs; cout8.噴嘴環(huán)出口氣流速度c1=C1endl;double k1_=pow(1/s,0.286);double k2_=pt+(1-pt)*k1_;double k3_=pow(k2_,3.5)

20、;double p1=p0*k3_; cout9.噴嘴環(huán)出口氣流壓力p1=p1endl;double T1=T0-x*x/1005*(1-pt)*vhs;cout10.噴嘴環(huán)出口氣流溫度T1=T1endl;double j1=p1/287/T1;cout11.噴嘴環(huán)出口氣流密度j1=j1endl;double n=1/(1-log10(T1/T0)/log10(p1/p0);cout12.氣流在噴嘴環(huán)中的多變指數(shù)n=nendl;double pcrn1=pow(2/(n+1),n/(n-1);cout13.臨界壓力比pcrn=pcrn1endl;double p11=p1/p0;cout14.

21、噴嘴環(huán)出口壓力與渦輪進(jìn)口壓力之比p11=p11endl;coutpcrn1是否小于p11endl;double u1=x0*Cs;cout15.葉輪進(jìn)口圓周速度u1=u1endl;double N1=60*u1/3.1415/D1;cout16.渦輪轉(zhuǎn)速N1=N1endl;double C1u=M*u1;cout17.葉輪進(jìn)口絕對(duì)速度的周向分量C1u=C1uendl;double C1r=C1*sin(a1);cout18.葉輪進(jìn)口絕對(duì)速度的徑向分量C1r=C1rendl;double w1=sqrt(C1*C1+u1*u1-2*C1*u1*cos(a1);cout19.葉輪進(jìn)口相對(duì)速度w1=

22、w1endl;double w1u=u1-C1u;cout20.葉輪進(jìn)口相對(duì)速度的周向分量w1=w1uendl;double B1=1.5708+atan(w1u/C1r);cout21.葉輪進(jìn)口氣流角B1=B1endl;double u2=u*u1;cout22.葉輪出口周向速度u2=u2endl;double w2s=sqrt(pt*Cs*Cs+w1*w1-u1*u1+u2*u2);cout23.葉輪出口理論相對(duì)速度w2s=w2sendl;double w2=y*w2s;cout24.葉輪出口實(shí)際相對(duì)速度w2=w2endl;double o=w2/w1; cout25.氣流在葉輪中的加速因

23、子o=oendl;double C2=sqrt(w2*w2+u2*u2-2*w2*u2*cos(B2);cout26.葉輪出口絕對(duì)速度c2=C2endl;double a2=w2*sin(B2/C2);cout27.葉輪出口氣流角的sin值a2endl;double T2=T1-pt*vhs/1005+(1-y*y)*w2s*w2s/2/1005;cout28.葉輪出口溫度T2=T2endl;double j2=p3/T2/287;cout29.葉輪出口空氣密度j2=j2endl;double en=(1-pt)*(1-x*x);cout30.噴嘴環(huán)相對(duì)能量損失en=enendl;double

24、 ei=(1-y*y)*(w2s/Cs)*(w2s/Cs);cout31.葉輪相對(duì)能量損失ei=eiendl;double ec2=(C2/Cs)*(C2/Cs);cout32.余速相對(duì)能量損失ec2=ec2endl;double nu=1-(en+ei+ec2);cout33.輪周效率nu=nuendl;double Ndf=736*1.18*j1*D1*D1*(pow(u1/100,3);cout34.輪盤摩擦損失功率Ndf=Ndfendl;double edf=Ndf/G/vhs;cout35.輪盤摩擦相對(duì)能量損失edf=edfendl;double zb=0.023; double n

25、_d=1-1.3*zb;cout36.間隙漏氣修正系數(shù)n_d=n_dendl;double nt=(nu-edf)*n_d;cout37.渦輪絕熱效率nt=ntendl;double vT0=vhs/1005;cout38.渦輪理論溫降vT0=vT0endl;double vT=vT0*nt;cout39.渦輪實(shí)際溫降vT=vTendl;double T3=T0-vT;cout40.渦輪出口速度T3=T3endl;double Nt=G*vhs*nt;cout41.渦輪功率Nt=Ntendl;double Nte=Nt*0.97;cout42.渦輪有效功率Nte=Nteendl;double

26、j3=p3/T3/287;cout43.渦輪出口空氣密度j3=j3endl;double Q=G/j3;cout44.渦輪出口空氣體積流量Q=Qendl;double w=3.1415*N1/30;cout45.葉輪角速度w=wendl;double Ns=w*sqrt(Q)/(pow(vhs,0.75);cout46.渦輪比轉(zhuǎn)速Ns=Nsendl;double Fmin=sqrt(287*T0)*G*1000000/p0/sqrt(7*(pow(p1/p0,2/n)-(pow(p1/p0,(n+1)/n);cout47.噴嘴環(huán)最小截面積Fmin=Fminendl;double D1n=D1*

27、1000+2;cout48.噴嘴環(huán)內(nèi)徑D1n=D1nendl;double F1n=G*1000000/C1/sin(a1)/0.81/j1;/第49步cout49.噴嘴環(huán)出口總面積F1n=F1nendl;double ln=F1n/3.1415/D1n;cout50.噴嘴環(huán)出口葉片高度ln=lnendl;double bn=Fmin/zn/ln;cout51.噴嘴喉部寬度bn=bnendl;double h1=1-1.5*zi/3.14/D1/1000;cout52.葉輪進(jìn)口阻塞系數(shù)h1=h1endl;double l1=G*1000/3.1415/D1/C1r/j1/h1;coutl1=l

28、1endl;coutl1?ln l1=l1 ln=lnendl; l1=1.18*ln;cout53.葉輪進(jìn)口葉片高度l1=l1endl;double D2=u*D1*1000;cout54.葉輪出口平均直徑D2=D2endl;double D2h=0.25*D1*1000; /預(yù)取D2h/D1=0.25cout55.葉輪出口葉根直徑D2h=D2hendl;double F2=3.1415/2*(D2*D2-D2h*D2h);cout56.葉輪出口總面積F2=F2endl;double D2t=sqrt(D2*D2+2*F2/3.1415);cout57.葉輪出口葉尖直徑D2t=D2tendl

29、;double l2=(D2t-D2h)/2;cout58.葉輪出口葉片高度l2=l2endl;double F2e=G*1000000/C2/j2;/sina2=1cout59.葉輪出口有效面積F2e=F2eendl;double h2=F2e/F2;cout60.葉輪出口的阻塞系數(shù)h2=h2endl;double v2=(F2-F2e)*sin(B2-0.0349066)/l2/zi*2;cout61.葉輪出口葉片尾緣厚度v2=v2endl;double lav=(l1+l2)/2;cout62.葉輪葉片高度lav=lav 0.23/lav與第36步預(yù)取值比較endl;結(jié)果1.葉輪葉片數(shù)z

30、i=182.速度比M=0.8888893.反力度pt=0.5292054.渦輪速比X0=0.7121095.渦輪膨脹比s=3.272736.渦輪總等焓降vhs=1040477.等熵速度Cs=456.1728.噴嘴環(huán)出口氣流速度c1=303.619.噴嘴環(huán)出口氣流壓力p1=21635710.噴嘴環(huán)出口氣流溫度T1=314.1411.噴嘴環(huán)出口氣流密度j1=2.3997512.氣流在噴嘴環(huán)中的多變指數(shù)n=1.3654213.臨界壓力比pcrn=0.53417714.噴嘴環(huán)出口壓力與渦輪進(jìn)口壓力之比p11=0.6009920.534177是否小于0.600992 /小于15.葉輪進(jìn)口圓周速度u1=324.84416.渦輪轉(zhuǎn)速N1=51702.1 /小于9000017.葉輪進(jìn)口絕對(duì)速度的周向分量C1u=288.75118.葉輪進(jìn)口絕對(duì)速度的徑向分量C1r=93.820719.葉輪進(jìn)口相對(duì)速度w1=100