该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
2元物系泡点和露点计算的通用程序!
老师给我一道例题:要求不必编计算程序,要编通用程序
计算任意2种物系组成的泡点和露点!
(1)含有80%(mol)醋酸乙酯(A)和20%乙醇(E)的二元物系,液相活度系数用Van Laar方程计算, =0.144, =0.170。试计算在101.3kPa压力下的露点温度和泡点温度。
我把露点温度计算程序发一下,泡点一样,只是多加几步就出来了,
#include
#include
main()
{ int i,k;
float x1,x2,a[6],tb,P,p[2],x,A1,A2,m,r,y1,y2,c,d;
printf("p,x1,x2\n");
scanf("%f,%f,%f",&P,&x1,&x2);
for(i=0;i<6;i++)
{ printf("a[%d]=",i);
scanf("%f",&a[i]);
printf("%f",a[i]);
printf(" %d\n",i); }
printf("A1,A2\n");
scanf("%f,%f",&A1,&A2);
m=((A1*x1)/(A2*x2))*((A1*x1)/(A2*x2));
r=exp(A1/(1+m));
printf("r=%f\n",r);
do
{printf("tb=");
scanf("%f",&tb);
p[0]=exp(a[0]-a[1]/(tb+273.15-a[2]));
printf("p[0]=%f\n",p[0]);
p[1]=exp(a[3]-a[4]/(tb+273.15-a[5]));
printf("p[1]=%f\n",p[1]);
x=(P-p[0])/(p[1]-p[0]);
printf("x=%f\n",x);
c=x/x2;
d=fabs(c-1);
printf("d=%f\n",d);
}while(d>0.01);
y2=r*x*(p[1]/P);
y1=(1-x)*(p[0]/P);
printf("tb=%f\nx1=%f,y1=%f x2=%f,y2=%f\n",tb,1-x,y1,x,y2);
getch();
}