±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 720  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 0854Çóµ÷¼Á +3 assdll 2026-04-03 3/150 2026-04-03 14:32 by fxue1114
[¿¼ÑÐ] 320Çóµ÷¼Á +3 ũҵ¹¤³ÌÓëÐÅÏ¢¼ 2026-04-03 3/150 2026-04-03 11:40 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] ר˶ 351 086100 Ò²ÊÇ¿¼µÄ²Ä¿Æ»ù ±¾¿ÆÒ²ÊDzÄÁÏ +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 295Çóµ÷¼Á +7 Ô¸ÂÃ;ÓÀԶ̹Ȼ 2026-04-02 7/350 2026-04-03 08:22 by fangshan711
[¿¼ÑÐ] 346Çóµ÷¼Á +5 Ö£³ÏÀÖ 2026-04-02 5/250 2026-04-02 16:38 by SZW_UJN
[¿¼ÑÐ] Ò»Ö¾Ô¸Ö£´ó²ÄÁϹ¤³Ì290Çóµ÷¼Á +20 Youth_ 2026-03-30 20/1000 2026-04-02 14:48 by 5896
[¿¼ÑÐ] 08¹¤¿ÆÇóµ÷¼Á290·Ö +5 1314Åõ»¨ 2026-04-02 8/400 2026-04-02 13:16 by ÇÇßÕßÕßÕ
[¿¼ÑÐ] ²ÄÁÏ»¯¹¤340Çóµ÷¼Á +5 jhx777 2026-03-30 5/250 2026-04-02 12:45 by smileboy2006
[¿¼ÑÐ] 292·Ö£¬²ÄÁÏÓ뻯¹¤£¬ÉêÇëµ÷¼Á +19 ³ÌÇçÖ® 2026-04-01 23/1150 2026-04-02 11:59 by ³ÌÇçÖ®
[¿¼ÑÐ] Ò»Ö¾Ô¸346ÉϺ£´óѧÉúÎïѧ +3 ÉϺ£´óѧ346µ÷¼Á 2026-04-01 3/150 2026-04-02 08:36 by w³æ³æ123
[¿¼ÑÐ] Ò»Ö¾Ô¸±±¾©¿Æ¼¼£¬085601×Ü·Ö305Çóµ÷¼Á +9 °ëÉú¹Ï£¡ 2026-04-01 11/550 2026-04-02 08:28 by Wang200018
[¿¼ÑÐ] 302Çóµ÷¼ÁÒ»Ö¾Ô¸±±º½070300£¬±¾¿ÆÖ£´ó»¯Ñ§ +8 Ê¥ÈÕ¶úÂüÌõ 2026-04-01 11/550 2026-04-02 07:40 by chemdavid
[¿¼ÑÐ] Çóµ÷¼Á +4 ͼ¼ø212 2026-03-30 5/250 2026-04-01 15:32 by ͼ¼ø212
[¿¼ÑÐ] 085410È˹¤ÖÇÄÜ ³õÊÔ316·Ö Çóµ÷¼Á +3 ²ÐÐÇ·÷Êï 2026-03-31 3/150 2026-04-01 11:09 by СÐÜraider
[¿¼ÑÐ] 08¹¤¿Æ£¬295£¬½ÓÊÜ¿çרҵµ÷¼Á +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by ÄæË®³Ë·ç
[¿¼ÑÐ] 267Çóµ÷¼Á +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by ̽123
[¿¼ÑÐ] 0855»úе³õÊÔ280Çóµ÷¼Á +3 kazenotori 2026-03-31 3/150 2026-04-01 10:08 by JourneyLucky
[¿¼ÑÐ] 080200ѧ˶£¬»úе¹¤³Ìרҵ277·Ö£¬Çó´ø×ߣ¡ +4 Æ¿×ÓPZ 2026-03-31 4/200 2026-03-31 20:16 by vgtyfty
[¿¼ÑÐ] 315Çóµ÷¼Á +6 akie... 2026-03-28 7/350 2026-03-31 16:48 by asdfzly
[¿¼ÑÐ] 332Çó92µ÷¼Á +8 ½¶½¶123 2026-03-28 8/400 2026-03-29 10:46 by ÖÜè÷µ¤
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û