油氣儲(chǔ)運(yùn)網(wǎng)

 找回密碼
 立即注冊(cè)

QQ登錄

QQ快速登錄

一步迅速開始

pdms培訓(xùn)&項(xiàng)目數(shù)字化360網(wǎng)
查看: 6742|回復(fù): 4
打印 上一主題 下一主題

[軟件分享] C++油氣集輸程序

[復(fù)制鏈接]
跳轉(zhuǎn)到指定樓層
樓主
發(fā)表于 2013-5-3 09:42:48 | 只看該作者 回帖獎(jiǎng)勵(lì) |倒序?yàn)g覽 |閱讀模式
#include <stdio.h>
#include <math.h>
int main()
{
    int i,temp=1,flag=0;
        double sum=0,sumC=0,Mpj=0;
        double z[8],w[8]={0};
        double C[8]= {0,25.9,6.6,12.3,12.4,12.4,16.5,913.9};
        int M[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/Mpj;
        for(i=1;i<=7;i++)
        {
                z[i]=(Mpj*w[i])/M[i];
        }
    double sumx=0,Mid,sumz=0;
        double y[8],sumy=0;
    double WL[8]={0},sumWL=0,Wg[8]={0},sumWg=0;
        double x[8]={0};
        double K1[8]={0,100,17,5.2,2.4,0.8,0.23,0.018};
        double K2[8]={0,255,43,13,4.95,1.65,0.6,0.0375};
        for(flag=0;flag<2;flag++)
        {
                double L1=0,V2=1;
                sumx=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;
                        else
                                V2=Mid;

                }
                printf("在第%d級(jí)分離后,組成為:  L=%f,  V=%f\n\n",flag+1,Mid,1-Mid);

    sumy=0;
        for(i=1;i<=7;i++)
        {
                y[i]=K1[i]*x[i];
                sumy+=y[i];
        }
    double Pg;
        sum=0;
        for(i=1;i<7;i++)
                sum+=y[i]*M[i];
        sum=sum+y[7]*114;
        Pg=sum/22.4;
        printf("在第%d級(jí)分離后,氣相密度為: %f \n\n",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("         Ci       Xi        Yi        Zi       WLi      Wgi\n");
        for(i=1;i<=6;i++)
        {
                printf("\n");
                printf(" C%d %f %f %f %f %f %f\n",i,C[i],x[i],y[i],z[i],WL[i],Wg[i]);
        }
        printf("\n");
        printf(" C7+%f %f %f %f %f %f\n",C[7],x[7],y[7],z[7],WL[7],Wg[7]);
        printf("\n");
        printf("合計(jì)%f %f %f %f %f %f\n\n",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;
}


油氣集輸程序.rar

913 Bytes, 下載次數(shù): 23

售價(jià): 3 金幣  [記錄]


油氣儲(chǔ)運(yùn)網(wǎng) - 論壇版權(quán)1、本主題所有言論和圖片純屬會(huì)員個(gè)人意見,與本論壇立場(chǎng)無關(guān)
2、本站所有主題由該帖子作者發(fā)表,該帖子作者與油氣儲(chǔ)運(yùn)網(wǎng)享有帖子相關(guān)版權(quán)
3、其他單位或個(gè)人使用、轉(zhuǎn)載或引用本文時(shí)必須同時(shí)征得該帖子作者和油氣儲(chǔ)運(yùn)網(wǎng)的同意
4、帖子作者須承擔(dān)一切因本文發(fā)表而直接或間接導(dǎo)致的民事或刑事法律責(zé)任
5、本帖部分內(nèi)容轉(zhuǎn)載自其它媒體,但并不代表本站贊同其觀點(diǎn)和對(duì)其真實(shí)性負(fù)責(zé)
6、如本帖侵犯到任何版權(quán)問題,請(qǐng)立即告知本站,本站將及時(shí)予與刪除并致以最深的歉意
7、油氣儲(chǔ)運(yùn)網(wǎng)管理員和版主有權(quán)不事先通知發(fā)貼者而刪除本文

沙發(fā)
發(fā)表于 2015-1-8 09:05:54 | 只看該作者
正需要,非常感謝樓主。
板凳
發(fā)表于 2017-4-22 09:21:26 | 只看該作者
謝謝樓主分享!
5#
發(fā)表于 2017-6-1 17:21:05 | 只看該作者
謝謝分享,下載學(xué)習(xí)
您需要登錄后才可以回帖 登錄 | 立即注冊(cè)  

本版積分規(guī)則

掃碼訪問手機(jī)版

QQ|關(guān)于我們|網(wǎng)站地圖|油氣儲(chǔ)運(yùn)網(wǎng) ( 魯ICP備11007657號(hào)-3 )

GMT+8, 2024-12-23 22:09 , Processed in 0.025749 second(s), 14 queries , Gzip On, MemCached On.

Powered by Discuz! X3.5

Copyright © 2001-2023, Tencent Cloud.

快速回復(fù) 返回列表