±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 721  |  »Ø¸´: 12

zyj8119

ľ³æ (ÖøÃûдÊÖ)


[½»Á÷] ¡¾ÇóÖú¡¿Çë¸ßÊÖÖ¸³öһϴ˳ÌÐòÓжàÉÙ´íÎó£¿

CODE:
#include
#include

#define M 16.043
#define omiga 0.012
#define Tc  190.6
#define Pc 45.99
#define R 83.14
#define epsilon 1-pow(2,0.5)
#define sigma 1+pow(2,0.5)
#define        alpha (1+(0.37464+1.54226*omiga-0.26992*pow(omiga,2)*pow(1-pow((T/Tc),(0.5))),2)
#define   a ((pow(0.45724*R,2)*pow(Tc,2)/Pc)*alpha
#define   b 0.07779*R*Tc/Pc
#define   beta b*P/(R*T)
#define   q a/(b*R*T)
int fugacity(float P, int T)
{   int s,I,N;
    double phi,f,P1;
    s=sizeof(P);
        for(int k=0;k<=s;k++){
    double z0(k)=M(s,1)=0;
        double z1(k)=N(s,1)=1;
    if((z0(k)-z1(k)).le.1e-6){
     z0(k)=z1(k);
        z1(k)=1+beta-q*beta*(z0(k)-beta)/((z0(k)+epsilon*beta)*(z0(k)+sigma*beta));
        }
        }
    I=(1/(sigma-epsilon))*log((z1(k)+sigma*beta)/(z1(k)+epsilon*beta));
    phi=exp(z1(k)-1-log(z1(k)-beta)-q*I);
    f=phi*P1;
        return f;
}

void main()
{
        fugacity(1000,298);
}

»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

» ÇÀ½ð±ÒÀ²£¡»ØÌû¾Í¿ÉÒԵõ½:

²é¿´È«²¿É¢½ðÌù

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

imyourkobe

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)


¡ï
ajian04(½ð±Ò+1):лл²ÎÓë½»Á÷~ 2010-11-24 20:48:45
ºÇºÇ£¬´íÎóºÜ¶àµÄ°¡¡£Äãµ½µ×ÒªÎÊiʲô£¬ÊÇÐèÒª°ïÖúÐ޸ĴíÎó£¿
2Â¥2010-11-24 18:35:29
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

imyourkobe

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)


zyj8119(½ð±Ò+5):ÄÇô³ÌÐòÓ¦¸ÃÔõô¸ÄÄØ£¿¸ßÊÖ£¿ 2010-11-24 19:09:15
ajian04:лл²ÎÓë½»Á÷~ 2010-11-24 20:48:53
ÆäÖУºz0(k),z1(k)Òª¶¨ÒåΪÊý×é
double z0[100],z1[100];
double z0(k)=M(s,1)=0;³õʼ»¯ÒªÖØÐ¸Äд¡£
ÀýÈ磺
for(int i=0;i z0=0
3Â¥2010-11-24 18:39:23
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zyj8119

ľ³æ (ÖøÃûдÊÖ)


ajian04:лл²ÎÓë½»Á÷~ 2010-11-24 20:48:59
ÒýÓûØÌû:
Originally posted by imyourkobe at 2010-11-24 18:35:29:
ºÇºÇ£¬´íÎóºÜ¶àµÄ°¡¡£Äãµ½µ×ÒªÎÊiʲô£¬ÊÇÐèÒª°ïÖúÐ޸ĴíÎó£¿

Ï£ÍûÄú°ïæÐÞ¸ÄÏ£¬ÎÒÊÇÐÂÊÖ¡£¡£¡£¡£¡£
4Â¥2010-11-24 19:10:39
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

379631818

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)


¡ï
zyj8119(½ð±Ò+1):лл²ÎÓë
ajian04:лл²ÎÓë½»Á÷~ 2010-11-24 20:49:04
zyj8119(½ð±Ò+2):ÊÇÒ»¸öËãÒݶȵijÌÐò¡£¡£¡£ 2010-11-24 21:05:36
µ½µ×ÒªËãʲôµÄ£¬Â¥Ö÷»áµÄÌ«¶àÁË
5Â¥2010-11-24 20:43:54
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

imyourkobe

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)


#defineÖÐÉÙÁË ¶ÔTºÍPµÄ¶¨Òå¡£
ÆäËûºÃ¸ÄµÄ¡£
6Â¥2010-11-24 21:11:12
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zyj8119

ľ³æ (ÖøÃûдÊÖ)


ÒýÓûØÌû:
Originally posted by imyourkobe at 2010-11-24 21:11:12:
#defineÖÐÉÙÁË ¶ÔTºÍPµÄ¶¨Òå¡£
ÆäËûºÃ¸ÄµÄ¡£

µ«ÊÇTºÍPÊÇÖ÷³ÌÐò¶¨ÒåµÄ°¡¡£
7Â¥2010-11-24 21:19:15
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

imyourkobe

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)


ÒýÓûØÌû:
Originally posted by zyj8119 at 2010-11-24 21:19:15:

µ«ÊÇTºÍPÊÇÖ÷³ÌÐò¶¨ÒåµÄ°¡¡£

Äã ÉÏÃæ
#define   beta b*P/(R*T)
#define   q a/(b*R*T)

Óõ½ÁËT ºÍP£¬Õâ˵Ã÷TºÍPÊdz£Á¿¡£
ÓëÖ÷³ÌÐòÖÐ񻂿·Ö¿ª¡£
8Â¥2010-11-24 21:23:31
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zyj8119

ľ³æ (ÖøÃûдÊÖ)


ÒýÓûØÌû:
Originally posted by imyourkobe at 2010-11-24 21:23:31:


Äã ÉÏÃæ
#define   beta b*P/(R*T)
#define   q a/(b*R*T)

Óõ½ÁËT ºÍP£¬Õâ˵Ã÷TºÍPÊdz£Á¿¡£
ÓëÖ÷³ÌÐòÖÐ񻂿·Ö¿ª¡£

ÄǸÃÔõô¸ÄÄØ£¿
9Â¥2010-11-24 21:29:57
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

imyourkobe

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)


zyj8119(½ð±Ò+50):¶àлÐֵܣ¡£¡£¡£¡£¡£¡£¡£¡ 2010-11-24 21:50:00
#include
#include

#define M 16.043
#define omiga 0.012
#define Tc  190.6
#define T  190.6
#define Pc 45.99
#define P 45.99
#define R 83.14
#define epsilon 1-pow(2,0.5)
#define sigma 1+pow(2,0.5)
#define        alpha 1+0.37464+1.54226*omiga-0.26992*pow(omiga,2)*pow(1-pow(T/Tc,0.5),2)
#define   a ((pow(0.45724*R,2)*pow(Tc,2)/Pc)*alpha
#define   b 0.07779*R*Tc/Pc
#define   beta b*P/(R*T)
#define   q a/(b*R*T)

int fugacity(float PP, int TT)
{
        int s,I,N;
double phi,f,P1;
//s=sizeof(P);
double z0[1000]={0};
double z1[1000]={1};
for(int k=0;k<=PP;k++)
{

    if((z0[k]-z1[k])<0.000001)
        {
                z0[k]=z1[k];
        z1[k]=1+beta-q*beta*(z0[k]-beta)/(z0[k]+epsilon*beta)*(z0[k]+sigma*beta));
        }
}
I=1/(sigma-epsilon)*log((z1[k]+sigma*beta)/(z1[k]+epsilon*beta));
phi=exp(z1[k]-1-log(z1[k]-beta)-q*I));
f=phi*P1;


return f;
}

void main()
{
        fugacity(1000,298);
}
10Â¥2010-11-24 21:48:09
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zyj8119

ľ³æ (ÖøÃûдÊÖ)


ÒýÓûØÌû:
Originally posted by imyourkobe at 2010-11-24 21:48:09:
#include
#include

#define M 16.043
#define omiga 0.012
#define Tc  190.6
#define T  190.6
#define Pc 45.99
#define P 45.99
#define R 83.14
#define epsilon ...

Ã²ËÆ»¹ÊÇÓм¸¸öWARNING:
--------------------Configuration: bbb - Win32 Debug--------------------
Compiling...
bbb.cpp
C:\Program Files\Microsoft Visual Studio\Common\MSDev98\My Projects\bbb.cpp(35) : warning C4244: '=' : conversion from 'double' to 'int', possible loss of data
C:\Program Files\Microsoft Visual Studio\Common\MSDev98\My Projects\bbb.cpp(40) : warning C4244: 'return' : conversion from 'double' to 'int', possible loss of data
C:\Program Files\Microsoft Visual Studio\Common\MSDev98\My Projects\bbb.cpp(21) : warning C4101: 's' : unreferenced local variable
C:\Program Files\Microsoft Visual Studio\Common\MSDev98\My Projects\bbb.cpp(21) : warning C4101: 'N' : unreferenced local variable
C:\Program Files\Microsoft Visual Studio\Common\MSDev98\My Projects\bbb.cpp(37) : warning C4700: local variable 'P1' used without having been initialized

bbb.obj - 0 error(s), 5 warning(s)
11Â¥2010-11-24 21:51:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zyj8119

ľ³æ (ÖøÃûдÊÖ)


ÒýÓûØÌû:
Originally posted by imyourkobe at 2010-11-24 21:48:09:
#include
#include

#define M 16.043
#define omiga 0.012
#define Tc  190.6
#define T  190.6
#define Pc 45.99
#define P 45.99
#define R 83.14
#define epsilon ...

Ã²ËÆÔËÐнá¹ûʲôҲûÓа¡¡£¡£¡£¡£
12Â¥2010-11-24 22:23:06
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

imyourkobe

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)


zyj8119(½ð±Ò+12):·Ç³£¸Ðл£¡£¡£¡£¡£¡ 2010-11-25 00:35:11
ÒýÓûØÌû:
Originally posted by zyj8119 at 2010-11-24 22:23:06:

Ã²ËÆÔËÐнá¹ûʲôҲûÓа¡¡£¡£¡£¡£

Ðֵܰ¡¡£ÄãÕâ¸ö³ÌÐòÖ÷º¯Êý¼ÓÉÏÊä³öº¯Êýprintf("%d",º¯Êý·µ»ØÖµ);
ÁíÍ⣬ÉÏÃæ¶¨ÒåµÄPºÍTµÄ³£Á¿µÄÖµÄãÒªÐÞ¸Äһϰ¡¡£
13Â¥2010-11-25 00:16:21
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ zyj8119 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 321Çóµ÷¼Á +10 ÈÏÕæÇóÉÏѧ 2026-04-02 10/500 2026-04-03 16:17 by lijunpoly
[¿¼ÑÐ] 357Çóµ÷¼Á +3 1050389037 2026-04-03 3/150 2026-04-03 15:51 by ÒÅÍüÏûʧµÄž™
[¿¼ÑÐ] Çó²ÄÁϵ÷¼Á£¬Ò»Ö¾Ô¸Ö£ÖÝ´óѧ289·Ö +14 ˶ÐǸ° 2026-04-03 14/700 2026-04-03 15:48 by Çç¿Õ210210
[¿¼ÑÐ] 0856£¬269·ÖÇóµ÷¼Á +14 ÓÐѧÉϾÍÐÐÇóÇóÁ 2026-03-30 17/850 2026-04-03 15:12 by ×ÁçíØ¼
[¿¼ÑÐ] 326·ÖÇóµ÷¼Á +3 ÓÚÊǺõÄØ 2026-04-01 5/250 2026-04-03 14:23 by ÓÚÊǺõÄØ
[¿¼ÑÐ] Çó²ÄÁϵ÷¼Á Ò»Ö¾Ô¸Äϲý´óѧ 328·Ö +5 yyy..... 2026-04-03 5/250 2026-04-03 13:46 by °ÙÁéͯ888
[¿¼ÑÐ] 321Çóµ÷¼Á +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[¿¼ÑÐ] 282Çóµ÷¼Á +5 ºôÎü¶¼ÊǼõ·Ê 2026-03-31 5/250 2026-04-03 12:03 by 1753564080
[¿¼ÑÐ] 313Çóµ÷¼Á +3 ¡«Î¢Î¢Á¹¡« 2026-04-03 3/150 2026-04-03 11:25 by à£à£à£0119
[¿¼ÑÐ] Çóµ÷¼Á22408 288·Ö +5 new382 2026-04-02 5/250 2026-04-03 09:13 by ×íÔÚ·çÀï
[¿¼ÑÐ] Çóµ÷¼Á£¡ÉúÎïÓëҽҩר˶ +4 Äæ×ªÂ½ÏÈÉú 2026-04-01 5/250 2026-04-03 08:33 by Jaylen.
[¿¼ÑÐ] Çóµ÷¼ÁÇóµ÷¼Á +7 121. 2026-04-02 7/350 2026-04-02 19:15 by dick_runner
[¿¼ÑÐ] 337Çóµ÷¼Á +11 ¡¶Ê÷¡· 2026-03-29 11/550 2026-04-02 10:20 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 11408 321·ÖÇóµ÷¼Á +3 huchun12138 2026-03-30 4/200 2026-04-01 22:48 by guanxin1001
[¿¼ÑÐ] °²È«¹¤³Ì 285 Çóµ÷¼Á +3 Xinyu56 2026-04-01 4/200 2026-04-01 21:50 by ¾²¾²¾²¾²¾²¾²¾²¾
[¿¼ÑÐ] Çóµ÷¼Á£¬Ò»Ö¾Ô¸±±ÁÖʳƷÓëÓªÑø095500£¬301·Ö£¬ÒѹýÁù¼¶£¬ÓпÆÑо­Àú +4 ¿ìÀÖ´¢Ðî¹Þ 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[¿¼ÑÐ] 333Çóµ÷¼Á +4 °¢¿ÆÒÝ 2026-03-31 4/200 2026-04-01 09:11 by jp9609
[¿¼ÑÐ] 080200ѧ˶£¬»úе¹¤³Ìרҵ277·Ö£¬Çó´ø×ߣ¡ +4 Æ¿×ÓPZ 2026-03-31 4/200 2026-03-31 20:16 by vgtyfty
[¿¼ÑÐ] Ò»Ö¾Ô¸Äϲý´óѧ324Çóµ÷¼Á +6 hanamiko 2026-03-29 6/300 2026-03-31 16:35 by hypershenger
[¿¼ÑÐ] 085601Ò»Ö¾Ô¸Î÷±±¹¤Òµ´óѧ³õÊÔ346 +4 085601³õÊÔ346 2026-03-30 4/200 2026-03-31 07:47 by jp9609
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û