基于热力学参数计算模型设计与实现.pdfVIP

#includestdio.h

#includestdlib.h

#includemath.h

floats(floatp,floatt)

{

floatsx;

sx=1.693*log(t)-0.4795*log(p)-2.9347+0.0005442*t+(0.02747-4.495*p)*pow((t/100),-

4.1)+(0.005731-2.51*pow(10,4)*pow(p,3))*pow((t/100),-14.5);

return(sx);

}

floath(floatp,floatt)

{

floathx;

hx=2018.24+1.693*t+0.0002721*pow(t,2)+(3.634-594.5*p)*pow((t/100),-3.1)+(0.6156-2.69

6*pow(10,6)*pow(p,3))*pow((t/100),-13.5);

return(hx);

}

floatt(floatp,floatsm)

{

floatt=273.149,s;

do

{

t=t+0.001;

s=1.693*log(t)-0.4795*log(p)-2.9347+0.0005442*t+(0.02747-4.495*p)*pow((t/100),-4.1)+(0

.005731-2.51*pow(10,4)*pow(p,3))*pow((t/100),-14.5);

}while(abs(10000*(s-sm))/10000.00.0001);

return(t);

}

voidcal()

{

floatn,x,p1,p2,t1,h1,h2,h3,h5,s1,s2,s3,s5,nex;

floatpa,pb,ta,tb,t2,ha,hb,sa,s6,s7,h6,h7;

printf(nexp1(MPa)t1(℃)p2(MPa)pa(MPa)ta(℃)\n);

scanf(%f%f%f%f%f%f,nex,p1,t1,p2,pa,ta);

t1=t1+273.15;

ta=ta+273.15;

p1=p1*10;

pa=pa*10;

h1=h(p1,t1);

ha=h(pa,ta);

#includestdio.h#includestdlib.h#includemath.h浮点s(浮点p,浮点t){浮点sx;

sx=1.693*log(t)‑0.4795*log(p)‑2.9347+0.0005442*t+(0.02747‑4.495*p)*pow((t/100),‑4.1)+(

0.005731‑2.51*pow(10,4)*pow(p,3))*pow((t/100),‑14.5);返回（SX）；}浮点h(浮点

p,浮点t){浮点hx;hx=2018.24+1.693*t+0.0002721*pow(t,2)+

(3.634‑594.5*p)*pow((t/100),‑3.1)+(0.6156‑2.69

6*pow(10,6)*pow(p,3))*pow((t/100),‑13.5);返回（hx）；}浮点t(浮点p,浮点sm){浮

点t=273.149,s;做{t=t+0.001;s=1.693*log(t)‑0.4795*log(p)‑2.9347+0.0005442*t+

(0.02747‑4.495*p)*pow((t/100),‑4.1)+(0.

005731‑2.51*pow(10,4)*pow(p,3))*pow((t/100),‑14.5);}

while(abs(10000*(s‑sm))/10000.00.0001);返回（t）；}voidcal(){float

n,x,p1,p2,t1,h1,h2,h3,h5,s1,s2,s3,s5,nex;浮动pa,pb,ta,tb,t2,ha,hb,sa,s6,s7,h6,h7;

printf(nexp1(MPa)t1(℃)p2(MPa)pa(MPa)ta(℃)\n);scanf(%f%f%f%f%f%f,

nex,p1,t1,p2,pa,ta);t1=t1+273.15；塔=塔+273.15;