理论塔板数的计算.doc

首先要得到相平衡方程和精馏段、提馏段方程，再根据逐板计算求得精馏塔的理论塔板数。源程序：#include #include #include static double R,L,F,W,T0,P,E;static double Xf,Xd,Xw,q,Xe;static double P1,P2;/td和tb的变量static double K1;static double K2;static double X1;static double X2;static double Y1;static double Y2;static double A1;static double B1;static double C1;static double A2;static double B2;static double C2;static double F1;static double F2;double Pressure(double A,double B,double C,double T)double temp1=101325.0/760.0;double temp2=A*1.0-B*1.0/(C+T)*1.0;double temp3=pow(10.0,temp2);double result=temp1*temp3;return result;double Ftd(double y,double P1,double P2)double temp1=pow(K1,2);double temp2=pow(C1+T0,2);double temp3=pow(K2,2);double temp4=pow(C2+T0,2);double temp5=temp1*P*temp2;double temp6=temp3*P*temp4;double temp7=-1.0*y*P1*B1*log(10);double temp8=-1.0*(1-y)*P2*B2*log(10);double result=temp7/temp5+temp8/temp6;return result;double Ftb(double x,double P1,double P2)double temp2=pow(C1+T0,2);double temp4=pow(C2+T0,2);double temp5=P*temp2;double temp6=P*temp4;double temp7=1.0*x*P1*B1*log(10);double temp8=1.0*(1-x)*P2*B2*log(10);double result=temp7/temp5+temp8/temp6;return result;double Tdm(double y)double Td=T0;double P1;double P2;KTdm:P1=Pressure(A1,B1,C1,Td); P2=Pressure(A2,B2,C2,Td); K1=P1*1.0/P; K2=P2*1.0/P; X1=y*1.0/K1; X2=(1-y)*1.0/K2; F1=X1+X2-1; F2=Ftd(y,P1,P2); if(fabs(F1)=E) Td=Td-F1*1.0/F2; goto KTdm; return Td;double Tbm(double x)double Tb=T0;double P1;double P2;KTbm:P1=Pressure(A1,B1,C1,Tb); P2=Pressure(A2,B2,C2,Tb); K1=P1/P; K2=P2/P; Y1=x*K1; Y2=(1-x)*K2; F1=Y1+Y2-1; F2=Ftb(x,P1,P2); if(fabs(F1)=E) Tb=Tb-F1*1.0/F2; goto KTbm; return Tb;void main()double Td;double Tb;double x;double y;int n=0;int m=0;int N,M;int Nt;int flag;char end;printf(Please input Rn);scanf(%lf,&R);printf(Please input Ln);scanf(%lf,&L);printf(Please input Fn);scanf(%lf,&F);printf(Please input Wn);scanf(%lf,&W);printf(Please input Pn);scanf(%lf,&P);printf(Please input T0n);scanf(%lf,&T0);printf(Please input En);scanf(%lf,&E);printf(Please input Xfn);scanf(%lf,&Xf);printf(Please input Xdn);scanf(%lf,&Xd);printf(Please input Xwn);scanf(%lf,&Xw);printf(Please input qn);scanf(%lf,&q);printf(Please input Xen);scanf(%lf,&Xe);printf(Please input A1n);scanf(%lf,&A1);printf(Please input B1n);scanf(%lf,&B1);printf(Please input C1n);scanf(%lf,&C1);printf(Please input A2n);scanf(%lf,&A2);printf(Please input B2n);scanf(%lf,&B2);printf(Please input C2n);scanf(%lf,&C2);/-/-y=Xd;flag=1;while(flag)flag=0;/-框框-x=0.0;Td=Tdm(y);P1=Pressure(A1,B1,C1,Td);P2=Pressure(A2,B2,C2,Td);K1=P1/P;K2=P2/P;Kup:Y1=x*K1;Y2=(1.0-x)*K2;if(fabs(Y1+Y2-1.0)=E)x=1.0*(x+E);goto Kup;/-n=n+1;if(x=Xe)flag=1;y=R*x/(R+1.0)+Xd/(R+1.0);N=n-1;/-/-x=Xw;flag=1;while(flag)flag=0;/-框框-y=1;Tb=Tbm(x);P1=Pressure(A1,B1,C1,Tb);P2=Pressure(A2,B2,C2,Tb);K1=P1/P;K2=P2/P;Kdown:X1=y/K1; X2=(1-y)/K2; if(fabs(X1+X2-1)=E) y=y-E; goto Kdown; /- m=m+1; if(xXf) fla