油气集输程序文本

《油气集输程序文本》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、#include#includeintmain(){inti,temp=1,flag=0;doublesum=0,sumC=0,Mpj=0;doublez[8],w[8]={0};doubleC[8]={0,25.9,6.6,12.3,12.4,12.4,16.5,913.9};intM[8]={0,16,30,44,58,72,86,472};for(i=1;i<=7;i++)sum+=C[i];sumC=sum;for(i=1;i<=7;i++){w[i]=C[i]/sum;Mpj+=w[i]/M[i];}Mpj=1

2、/Mpj;for(i=1;i<=7;i++){z[i]=(Mpj*w[i])/M[i];}doublesumx=0,Mid,sumz=0;doubley[8],sumy=0;doubleWL[8]={0},sumWL=0,Wg[8]={0},sumWg=0;doublex[8]={0};doubleK1[8]={0,100,17,5.2,2.4,0.8,0.23,0.018};doubleK2[8]={0,255,43,13,4.95,1.65,0.6,0.0375};for(flag=0;flag<2;flag++){doubleL1=0,V2=1;s

3、umx=0;while(fabs(sumx-1)>1e-6){Mid=(L1+V2)/2;sumx=0;for(i=1;i<=7;i++){x[i]=z[i]/(Mid+(1-Mid)*K1[i]);sumx+=x[i];}if(sumx-1>0)L1=Mid;elseV2=Mid;}printf("在第%d级分离后，组成为:L=%f,V=%f",flag+1,Mid,1-Mid);sumy=0;for(i=1;i<=7;i++){y[i]=K1[i]*x[i];sumy+=y[i];}doublePg;sum=0;for(i=1;i<7;i++

4、)sum+=y[i]*M[i];sum=sum+y[7]*114;Pg=sum/22.4;printf("在第%d级分离后，气相密度为:%f",flag+1,Pg);temp--;sumWL=0;sumWg=0;for(i=1;i<=7;i++){WL[i]=(x[i]*Mid*C[i])/z[i];sumWL+=WL[i];Wg[i]=C[i]-WL[i];sumWg+=Wg[i];}printf("CiXiYiZiWLiWgi");for(i=1;i<=6;i++){printf("");printf("C%d%f%f%f%f%f%f

5、",i,C[i],x[i],y[i],z[i],WL[i],Wg[i]);}printf("");printf("C7+%f%f%f%f%f%f",C[7],x[7],y[7],z[7],WL[7],Wg[7]);printf("");printf("合计%f%f%f%f%f%f",sumC,sumx,sumy,sumz,sumWL,sumWg);if(flag==0){for(i=1;i<=7;i++){K1[i]=K2[i];z[i]=x[i];sumz+=z[i];}}}sumy=0;}