水力計(jì)算程序

附錄一水力計(jì)算程序#include<stdio.h>#include"math.h”voidmain(){FILE*fp;intn,i,j,d,y,nt[4];charPGT[4][30];floatp1,p2,p3,pQ,pk,pkcp,YB,YB1,q,z,zkcp,DB,DH,x,o,K,L,Lkcp,Tcp,Tkcp,Tcp1,Tkcp1,T0,TQ;floatQ, Kcp,Cp,Cpkcp,Cp1,Cp2,aL,aL1,aL2,aLkcp,Di,Dik,Di1,Di2,Vs,mm,C,a,b,c,p1min,p2max,pjmax,pjmin;floatQ1,Q2,Q3,Q4,Q5,Q6,Q7,L1,L2,q1,q2,q3,q4,q6,q7,N,l,l1,l2,Tcp2,Tcp3,Tcp4,Tcp5,L3,L4,L5,L6,L7,L8,L9;floatR, k,A[4],G[4],zB,H,N0[4],Niso[4],Nsite[4],g[4],Ncomp[4],HRiso[4],HRsite[4],qF[4],Hl,QF[4],Qsumf,nn;printf("輸入設(shè)計(jì)壓力P1(Mp)\n");/*輸入一系列參數(shù)值*/scanf("%f”,&p1);printf("輸入壓比YB\n");scanf("%f",&YB);printf("輸入管外徑DH(mm)\n");scanf("%f",&DH);a=p1*DH/2/552/0.72;/*計(jì)算壁厚*/printf("輸出壁厚％%.1f(mm)\n",a);printf("輸入化整后的壁厚\n");scanf("%f",&b);K=0.005;DB=DH-2*b;/*計(jì)算內(nèi)徑*/p2=p1/YB;p3=0.666667*(p1+p2*p2/(p1+p2));/*計(jì)算平均壓力*/c=p3*DH/2/552/0.72;/*計(jì)算等強(qiáng)度壁厚*/printf("等強(qiáng)度設(shè)計(jì)壁厚％%.1f(mm)\n",c);pk=2.1;/*出站最小壓力*/pkcp=0.666667*(p1+pk*pk/(p1+pk));/*末段平均壓力*/z=100.0/(100.0+0.113*pow(10.0*p3,1.15));zkcp=100.0/(100.0+0.113*pow(10.0*pkcp,1.15));/*壓縮因子*/x=1.05*0.077*pow(K,0.2)/(0.9025*pow(DB,0.2));/*水力摩阻系數(shù)*/printf("水力摩阻系數(shù)x=%.3f\n",x);o=0.595;/*相對(duì)密度*/Tcp=293;Tkcp=290;Q=280;/*末段輸氣量*/q=Q/350*100;L=(105.113*pow(DB*0.001,2.5)/q)*(105.113*pow(DB*0.001,2.5)/q)*((p1*p1-p2*p2)/(o*x*z*Tcp));Lkcp=(105.113*pow(DB*0.001,2.5)/q)*(105.113*pow(DB*0.001,2.5)/q)*((p1*p1-pk*pk)/(o*x*zkcp*Tkcp))/2;/*0.計(jì)算末端平均溫度和站間距*/Kcp=1.75;T0=285;TQ=299;Cp=1.695+0.001838*Tcp+1960000*(p3-0.1)/pow(Tcp,3);Cpkcp=1.695+0.001838*Tkcp+1960000*(pkcp-0.1)/pow(Tkcp,3);aL=0.225*Kcp*DH*0.001*L/(q*o*Cp);aLkcp=0.225*Kcp*DH*0.001*Lkcp/(q*o*Cpkcp);Tcp1=T0+(TQ-T0)*(1-exp(-aL))/aL-Di*(p1*p1-p2*p2)*(1-(1-exp(-aL))/aL)/2/aL/p3;Tkcp1=T0+(TQ-T0)*(1-exp(-aLkcp))/aLkcp-Dik*(p1*p1-pk*pk)*(1-(1-exp(-aLkcp))/aLkcp)/2/aLkcp/pkcp;Di=1/Cp*(980000/Tcp1/Tcp1-1.5);Dik=1/Cpkcp*(980000/Tkcp1/Tkcp1-1.5);while(fabs(Tcp1-Tcp)>0.001&&fabs(Tkcp1-Tkcp)>0.001){Tcp=Tcp1;Tkcp=Tkcp1;L=pow((105.113*pow(DB*0.001,2.5)/q),2)*(p1*p1-p2*p2)/(o*x*z*Tcp);Lkcp=pow((105.113*pow(DB*0.001,2.5)/q),2)*(p1*p1-pk*pk)/(o*x*zkcp*Tkcp)/2;Cp=1.695+0.001838*Tcp+1960000*(p3-0.1)/pow(Tcp,3);Cpkcp=1.695+0.001838*Tkcp+1960000*(pkcp-0.1)/pow(Tkcp,3);aL=0.225*Kcp*DH*0.001*L/(q*o*Cp);aLkcp=0.225*Kcp*DH*0.001*Lkcp/(q*o*Cpkcp);Tcp1=T0+(TQ-T0)*(1-exp(-aL))/aL-Di*(p1*p1-p2*p2)*(1-(1-exp(-aL))/aL)/2/aL/p3;Tkcp1=T0+(TQ-T0)*(1-exp(-aLkcp))/aLkcp-Dik*(p1*p1-pk*pk)*(1-(1-exp(-aLkcp))/aLkcp)/2/aLkcp/pkcp;Di=1/Cp*(980000/Tcp1/Tcp1-1.5);Dik=1/Cpkcp*(980000/Tkcp/Tkcp-1.5);}Tcp=Tcp1;Tkcp=Tkcp1;printf("z=%.3f\n”,z);printf("zkcp=%.3f\n”,zkcp);printf("末段平均溫度％%.1f(K)\n",Tkcp);printf("計(jì)算段平均溫度％%.1f(K)\n”,Tcp);printf("末段平均壓力為Pkcp=%.2fMP\n",pkcp);printf("計(jì)算段平均壓力為Pcp=%.2fMP\n",p3);printf("末段長(zhǎng)度Lkcp=%.0fkm\n”,Lkcp);printf("末端站間距L=%.0fkm\n”,L);C=3.1415926/4*pow(DB/1000,2)*293*Lkcp/0.101325/Tkcp/zkcp*pow(10,3);p2max=sqrt(p1*p1-o*Lkcp*x*zkcp*Tkcp*q*q/105.113/105.113/pow(DB/1000,5));p1min=sqrt(pk*pk+o*Lkcp*x*zkcp*Tkcp*q*q/105.113/105.113/pow(DB/1000,5));printf("末站最小壓比=%.2f\n",p1min/p2);printf("p1min=%.1fMP\n”,p1min);printf("p2min=%.1fMP\n”,p2max);pjmax=0.666667*(p1+pow(p2max,2)/(p1+p2max));printf("pjmax=%.1fNm3\n”,pjmax);pjmin=2.000/3.000*(p1min+pow(pk,2)/(pk+p1min));printf("pjmin=%.1fNm3\n”,pjmin);Vs=C*(pjmax-pjmin);mm=Vs/q/10000;/*計(jì)算儲(chǔ)氣量，調(diào)峰量*/printf("儲(chǔ)氣量Vs=%.1fNm3\n",Vs);printf("儲(chǔ)氣率mm=%.1f%%\n”,mm);/*1.末端相關(guān)計(jì)算*/i=0;l=1350;while(l-Lkcp-i*L>L){i=i+1;}L1=l-Lkcp-i*L;/*計(jì)算L1*/printf("末段壓氣站數(shù)i=%d\n",i+1);printf("距后分氣點(diǎn)距離為L(zhǎng)1=%.0fkm\n",L1);Q3=35;Q1=Q;Q2=Q1+Q3;q1=Q1/350*100;q2=Q2/350*100;q3=Q3/350*100;N=o*x*z*Tcp/pow(105.113,2)/pow(DB/1000,5);L2=(p1*p1-p2*p2-L1*N*q1*q1)/N/q2/q2;/*計(jì)算L2*/printf("后分氣點(diǎn)前一段管長(zhǎng)L2=%.0f(km)\n”,L2);/*2.中間段相關(guān)計(jì)算*/Q4=Q1+Q3;Tcp4=293;q4=Q4/350*100;L3=(105.113*pow(DB*0.001,2.5)/q4)*(105.113*pow(DB*0.001,2.5)/q4)*((p1*p1-p2*p2)/(o*x*z*Tcp4));T0=285;TQ=299;Cp1=1.695+0.001838*Tcp4+1960000*(p3-0.1)/pow(Tcp4,3);aL1=0.225*Kcp*DH*0.001*L3/(q4*o*Cp1);Di1=1/Cp1*(980000/Tcp4/Tcp4-1.5);Tcp5=T0+(TQ-T0)*(1-exp(-aL1))/aL1-Di1*(p1*p1-p2*p2)*(1-(1-exp(-aL1))/aL1)/2/aL1/p3;while(fabs(Tcp4-Tcp5)>0.001)/*計(jì)算L4*/{Tcp4=Tcp5;L4=pow((105.113*pow(DB*0.001,2.5)/q4),2)*(p1*p1-p2*p2)/(o*x*z*Tcp4);Cp1=1.695+0.001838*Tcp4+1960000*(p3-0.1)/pow(Tcp4,3);aL1=0.225*Kcp*DH*0.001*L4/(q4*o*Cp1);Tcp5=T0+(TQ-T0)*(1-exp(-aL1))/aL1-Di1*(p1*p1-p2*p2)*(1-(1-exp(-aL1))/aL1)/2/aL1/p3;Di1=1/Cp1*(980000/Tcp5/Tcp5-1.5);}printf("輸出中間段站間距L4=%.0f(km)\n",L4);/*計(jì)算中間段站間距*/j=0;l1=700;while(l1-L2-j*L4>L4){j=j+1;}L5=l1-L2-j*L4;/*計(jì)算L5*/printf("輸出中間段溫度Tcp4=%.1f(k)\n",Tcp4);printf("中間段壓氣站數(shù)j=%d\n",j+1);printf("輸出中間段站間距L4=%.0f(km)\n",L4);printf("距離前分氣點(diǎn)距離L5=%.0fkm\n",L5);Q5=35;Q6=Q4+Q5;q6=Q6/350*100;q4=Q4/350*100;N=o*x*z*Tcp4/pow(105.113,2)/pow(DB/1000,5);L6=(p1*p1-p2*p2-L5*N*q4*q4)/N/q6/q6;/*計(jì)算L6*/printf("輸出前分氣點(diǎn)前一段管長(zhǎng)L6=%.0f(km)\n”,L6);/*3.首段相關(guān)計(jì)算*/Q7=Q6;Tcp2=292;q7=Q7/350*100;L7=(105.113*pow(DB*0.001,2.5)/q7)*(105.113*pow(DB*0.001,2.5)/q7)*((p1*p1-p2*p2)/(o*x*z*Tcp2));T0=285;TQ=299;Cp2=1.695+0.001838*Tcp2+1960000*(p3-0.1)/pow(Tcp2,3);aL2=0.225*Kcp*DH*0.001*L7/(q7*o*Cp2);Di2=1/Cp2*(980000/Tcp2/Tcp2-1.5);Tcp3=T0+(TQ-T0)*(1-exp(-aL2))/aL2-Di2*(p1*p1-p2*p2)*(1-(1-exp(-aL2))/aL2)/2/aL2/p3;while(fabs(Tcp2-Tcp3)>0.001)/*計(jì)算L8*/Tcp2=Tcp3;L8=pow((105.113*pow(DB*0.001,2.5)/q7),2)*(p1*p1-p2*p2)/(o*x*z*Tcp2);Cp2=1.695+0.001838*Tcp2+1960000*(p3-0.1)/pow(Tcp2,3);aL2=0.225*Kcp*DH*0.001*L8/(q7*o*Cp2);Tcp3=T0+(TQ-T0)*(1-exp(-aL2))/aL2-Di2*(p1*p1-p2*p2)*(1-(1-exp(-aL2))/aL2)/2/aL2/p3;Di2=1/Cp2*(980000/Tcp3/Tcp3-1.5);}/*計(jì)算首段平均溫度和站間距*/d=0;l2=3500;while(l2-L6-d*L8>L8){d=d+1;}L9=l2-L6-d*L8;/*計(jì)算L9*/printf("輸出首段溫度Tcp2=%.1f(k)\n",Tcp2);printf("首段壓氣站數(shù)d=%d\n",d+1);printf("輸出首段站間距L8=%.0f(km)\n",L8);y=i+j+d+3;/*壓氣站總站數(shù)計(jì)算*/printf("總壓氣站數(shù)y=%d\n",y);printf("首戰(zhàn)距離起點(diǎn)距離L9=%.0fkm\n",L9);pQ=sqrt(p2*p2+o*L9*x*z*Tcp2*q7*q7/pow(105.113,2)/pow(DB/1000,5));/*起點(diǎn)壓力計(jì)算*/YB1=pQ/4.1;/*起點(diǎn)壓比*/printf("起點(diǎn)升壓YB1=%.2f\n",YB1);/*4.燃?xì)夂牧坑?jì)算*/R=848/17.241;/*基本參數(shù)*/Hl=37484.3;k=1.5;A[1]=q;/*流量賦值*/A[2]=q4;A[3]=q7;zB=100.0/(100.0+0.113*pow(10.0*p2,1.15));/*壓縮機(jī)入口壓縮因子*/H=k/(k-1)*zB*R*293*(pow(YB,(k-1)/k)-1);/*計(jì)算多變能頭*/for(n=1;n<=3;n++){G[n]=A[n]*pow(10,10)*1.033/24/3600/R/293;/*計(jì)算質(zhì)量流量*/N0[n]=G[n]*H/102/0.8;/*計(jì)算單站功率*/printf（"輸出第％d段的壓氣站計(jì)算功率％%.0f（kw）\n”,n,N0[n]）;printf（"輸入燃?xì)廨啓C(jī)型號(hào)\n"）;scanf（"%s",PGT[n]）;/*壓縮機(jī)型號(hào)*/printf（"輸入第％d段的壓氣站在ISO下的額定功率（kw）\n",n）;scanf（"%f",&Niso[n]）;Nsite[n]=Niso[n]*1.04*0.952*0.985*0.995;/*計(jì)算可用功率*/nt[n]=N0[n]/Nsite[n]+1;/*計(jì)算一站燃-壓組數(shù)*/g[n]=G[n]/nt[n];/*每臺(tái)壓縮機(jī)排量*/Ncomp[n]=g[n]*H/102/0.8+100;/*壓縮機(jī)實(shí)際功率*/printf（"輸入燃?xì)廨啓C(jī)ISO下的額定熱耗率\n"）;scanf（"%f",&HRiso[n]）;HRsite[n]=HRiso[n]*0.995*1.0098*1.018;/*熱耗率*/qF[n]=HRsite[n]/Hl;/*單位時(shí)間燃?xì)夂牧?/QF[n]=qF[n]*24*365*0.875*nt[n]*Ncomp[n];/*單站燃?xì)饽旰牧?/}Qsumf=QF[1]*（i+1）+QF⑵*（j+1）+QF[3]*（d+1）;/*一種方案燃?xì)饪偤牧?/nn=Qsumf/350/pow（10,6）;/*全線燃?xì)夂牧空伎傒斄康谋壤?/printf（"全線燃?xì)夂牧空伎傒斄康谋壤齨n=%.1f%%\n”,nn）;fp=fopen（"E:\\初算計(jì)算結(jié)果.xls","a"）;/*創(chuàng)建文件到E盤(pán)*/if（fp==NULL）{printf（"創(chuàng)建文件失??！\n"）;return;}else{printf（"在E盤(pán)創(chuàng)建XLs文件記錄結(jié)果\n"）;fprintf（fp,”設(shè)計(jì)壓力（Mpa）\t壓比\t管道外徑（mm）\t管道壁厚（mm）\t內(nèi)涂層（mm）\t水力摩阻系數(shù)\t末端壓縮系數(shù)\t末段平均溫度（K）\t末段長(zhǎng)度（km）\t末段平均壓力（MPa）\t計(jì)算段平均壓力（MPa）\t調(diào)峰量\t計(jì)算段壓縮系數(shù)\t平