水力计算程序

附录一水力计算程序#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("输入设计压力P1(Mp)\n");/*输入一系列参数值*/scanf("%f”,&p1);printf("输入压比YB\n");scanf("%f",&YB);printf("输入管外径DH(mm)\n");scanf("%f",&DH);a=p1*DH/2/552/0.72;/*计算壁厚*/printf("输出壁厚％%.1f(mm)\n",a);printf("输入化整后的壁厚\n");scanf("%f",&b);K=0.005;DB=DH-2*b;/*计算内径*/p2=p1/YB;p3=0.666667*(p1+p2*p2/(p1+p2));/*计算平均压力*/c=p3*DH/2/552/0.72;/*计算等强度壁厚*/printf("等强度设计壁厚％%.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));/*水力摩阻系数*/printf("水力摩阻系数x=%.3f\n",x);o=0.595;/*相对密度*/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.计算末端平均温度和站间距*/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("计算段平均温度％%.1f(K)\n”,Tcp);printf("末段平均压力为Pkcp=%.2fMP\n",pkcp);printf("计算段平均压力为Pcp=%.2fMP\n",p3);printf("末段长度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;/*计算储气量，调峰量*/printf("储气量Vs=%.1fNm3\n",Vs);printf("储气率mm=%.1f%%\n”,mm);/*1.末端相关计算*/i=0;l=1350;while(l-Lkcp-i*L>L){i=i+1;}L1=l-Lkcp-i*L;/*计算L1*/printf("末段压气站数i=%d\n",i+1);printf("距后分气点距离为L1=%.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;/*计算L2*/printf("后分气点前一段管长L2=%.0f(km)\n”,L2);/*2.中间段相关计算*/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)/*计算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=0;l1=700;while(l1-L2-j*L4>L4){j=j+1;}L5=l1-L2-j*L4;/*计算L5*/printf("输出中间段温度Tcp4=%.1f(k)\n",Tcp4);printf("中间段压气站数j=%d\n",j+1);printf("输出中间段站间距L4=%.0f(km)\n",L4);printf("距离前分气点距离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;/*计算L6*/printf("输出前分气点前一段管长L6=%.0f(km)\n”,L6);/*3.首段相关计算*/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)/*计算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);}/*计算首段平均温度和站间距*/d=0;l2=3500;while(l2-L6-d*L8>L8){d=d+1;}L9=l2-L6-d*L8;/*计算L9*/printf("输出首段温度Tcp2=%.1f(k)\n",Tcp2);printf("首段压气站数d=%d\n",d+1);printf("输出首段站间距L8=%.0f(km)\n",L8);y=i+j+d+3;/*压气站总站数计算*/printf("总压气站数y=%d\n",y);printf("首战距离起点距离L9=%.0fkm\n",L9);pQ=sqrt(p2*p2+o*L9*x*z*Tcp2*q7*q7/pow(105.113,2)/pow(DB/1000,5));/*起点压力计算*/YB1=pQ/4.1;/*起点压比*/printf("起点升压YB1=%.2f\n",YB1);/*4.燃气耗量计算*/R=848/17.241;/*基本参数*/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));/*压缩机入口压缩因子*/H=k/(k-1)*zB*R*293*(pow(YB,(k-1)/k)-1);/*计算多变能头*/for(n=1;n<=3;n++){G[n]=A[n]*pow(10,10)*1.033/24/3600/R/293;/*计算质量流量*/N0[n]=G[n]*H/102/0.8;/*计算单站功率*/printf（"输出第％d段的压气站计算功率％%.0f（kw）\n”,n,N0[n]）;printf（"输入燃气轮机型号\n"）;scanf（"%s",PGT[n]）;/*压缩机型号*/printf（"输入第％d段的压气站在ISO下的额定功率（kw）\n",n）;scanf（"%f",&Niso[n]）;Nsite[n]=Niso[n]*1.04*0.952*0.985*0.995;/*计算可用功率*/nt[n]=N0[n]/Nsite[n]+1;/*计算一站燃-压组数*/g[n]=G[n]/nt[n];/*每台压缩机排量*/Ncomp[n]=g[n]*H/102/0.8+100;/*压缩机实际功率*/printf（"输入燃气轮机ISO下的额定热耗率\n"）;scanf（"%f",&HRiso[n]）;HRsite[n]=HRiso[n]*0.995*1.0098*1.018;/*热耗率*/qF[n]=HRsite[n]/Hl;/*单位时间燃气耗量*/QF[n]=qF[n]*24*365*0.875*nt[n]*Ncomp[n];/*单站燃气年耗量*/}Qsumf=QF[1]*（i+1）+QF⑵*（j+1）+QF[3]*（d+1）;/*一种方案燃气总耗量*/nn=Qsumf/350/pow（10,6）;/*全线燃气耗量占总输量的比例*/printf（"全线燃气耗量占总输量的比例nn=%.1f%%\n”,nn）;fp=fopen（"E:\\初算计算结果.xls","a"）;/*创建文件到E盘*/if（fp==NULL）{printf（"创建文件失败！\n"）;return;}else{printf（"在E盘创建XLs文件记录结果\n"）;fprintf（fp,”设计压力（Mpa）\t压比\t管道外径（mm）\t管道壁厚（mm）\t内涂层（mm）\t水力摩阻系数\t末端压缩系数\t末段平均温度（K）\t末段长度（km）\t末段平均压力（MPa）\t计算段平均压力（MPa）\t调峰量\t计算段压缩系数\t平