化工原理精馏塔设计及流程原理图 第11页

化工原理精馏塔设计及流程原理图 第11页
    fn=s-1; 
    fnd=(101.3/760)*(x0*B[0]*exp(A[0]-B[0]/(T+C[0]))/(pa*(T+C[0])*(T+C[0]))+(1-x0)*B[1]*exp(A[1]-B[1]/(T+C[1]))/(pa*(T+C[1])*(T+C[1])));
   }
 while(fabs(s-1)>0.001);
 printf("T=%f ",T);
 printf("pa=%f ",pa);
for(i=0;i<2;i++)
 {printf("p=%f ",p[i]);
printf("K=%f\n",K[i]);}
 a=K[0]/K[1];
 printf("a=%f ",a);
 
 
}

main()
{double f,c,c20,uf,u,A,At,Ad,b,Dc,D,hb,hw,lw,how,bs,bc,bd,Aa,x,r,fi,d0,Ao,uo,uo1,t,ev,co,g,ho,
test(316,1721.3,0.98);
a1=a;
do{
a2=a;
a=(a1+a2)/2;
xd=0.98;xw=0.02;qmd=0.5364;qmw=0.2939;qmf=0.8303;dl=470;dv=28,Ht=0.45;
xe=0.65;ye=a*xe/(1+(a-1)*xe);
Rmin=(xd-ye)/(ye-xe);
R=Rmin*1.6;
Y[0]=xd;
for(j=0;;j++)
{
X[j]=Y[j]/(a+(1-a)*Y[j]);
Y[j+1]=R*X[j]/(R+1)+xd/(R+1);
if(X[j]<0.65) {Nf=(int)(j/0.6)+1;break;}
}
for(;;j++)
{
X[j]=Y[j]/(a+(1-a)*Y[j]);
Y[j+1]=(R+1.5238)*X[j]/(R+1)-0.5238*xw/(R+1); 
if(X[j]<0.02) break;        
}
pb=1721.3+0.461*(int)((j-1)/0.6);
aa=a;
test(324,pb,0.02);
}
while(fabs(a-a2)>0.0001);
a=aa;
if((10*(j-1))%6!=0) N=(int)((j-1)/0.6)+1;
else N=10*(j-1)/6;
qml[0]=R*qmd;
qmv[0]=(R+1)*qmd;
qml[1]=qml[0]+qmf;
qmv[1]=qmv[0];
printf("a=%f ",a);
printf("R=%f ",R);
printf("Nf=%d ",Nf);
printf("N=%d\n",N);
for(i=0;i<2;i++)
{printf("qml=%f ",qml[i]);
 printf("qmv=%f ",qmv[i]);
 printf("\n");
}
Flv=(qml[0]/qmv[0])*sqrt(dv/dl);
printf("Flv=%f\n",Flv);

f=4.5;c20=0.061;b=0.15;
c=c20*pow((f/20),0.2);
uf=c*sqrt((dl-dv)/dv);
u=uf*0.7;
A=qmv[0]/(dv*u);
At=A/(1-b);
Dc=sqrt(4*At/3.14159);

printf("c=%f ",c);
printf("uf=%f ",uf);
printf("Dc=%f\n",Dc);

D=1.8;
At=0.25*3.14159*D*D;
A=(1-b)*At;
Ad=At-A;
u=qmv[0]/(A*dv);
m=u/uf;
b=Ad/At;

printf("At=%f ",At);
printf("A=%f ",A);
printf("Ad=%f ",Ad);
printf("u=%f ",u);
printf("m=%f ",m);
printf("b=%f\n",b);

bs=0.080;bc=0.050;d0=0.007;
bd=0.21*D;
r=D/2-bc;
x=D/2-(bd+bs);
辣文论文网
printf("Aa=%f ",Aa);
printf("t=%f ",t);
printf("fi=%f ",fi);
printf("Ao=%f ",Ao);
printf("uo=%f ",uo);
printf("n=%d\n",n);

d=0.003;
hb=0.045; 
hw=0.040;
lw=0.82*D;
how=0.00284*pow(3600*qml[0]/(dl*lw),0.666667);
printf("d=%f ",d);
printf("hb=%f ",hb);
printf("hw=%f ",hw);
printf("lw=%f ",lw);
printf("how=%f\n",how);

printf("开始校核\n");
Hf=2.5*(hw+how);
ev=0.0057/f*pow(u/(Ht-Hf),3.2);
printf("ev=%f\n",ev);

 << 上一页  [11] [12] [13] 下一页

化工原理精馏塔设计及流程原理图 第11页下载如图片无法显示或论文不完整，请联系qq752018766