资源描述：
《连续多组分三元）精馏塔的模拟计算.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、1.连续多组分（三元）精馏塔的模拟计算functionConDistillclearallclcglobalFz1z2z3Ralpha1alpha2alpha3M1MNMNtNfV1VDLL1WF=40;%进料流量,kmol/hrR=5;%回流比z1=0.6;%苯的进料组成（摩尔分率）z2=0.25;%甲苯的进料组成（摩尔分率）z3=1-z1-z2;%苯乙烯的进料组成（摩尔分率）M1=75;%塔顶冷凝器中的滞液量(kmol)M=10;%塔板上的滞液量(kmol)MN=150;%塔釜中的滞液量(kmol)q=1;%饱和进料tf=35;dt
2、=1;%相对挥发度alpha1=2.75;alpha2=1;alpha3=0.4;Nt=10;%塔板总数Nf=5;%进料位置V1=150;%从塔釜蒸发上来的蒸汽流量(kmol/hr)%精馏段V=V1-(1-q)*F;D=V/(R+1);L=V-D;%提馏段L1=L+F;W=L1-V1;%初始化x1和x2---开车时塔内所有板上的x1和x2分别与进料的z1和z2相同x1=z1*ones(1,Nt);x2=z2*ones(1,Nt);[t,y]=ode45(@DistMassBalances,[0:dt:tf],[x1x2])%输出结果x1
3、=y(:,1:Nt);%苯的液相组成（摩尔分率）x2=y(:,Nt+1:2*Nt);%甲苯的液相组成（摩尔分率）x3=1-x1-x2;%苯乙烯的液相组成（摩尔分率）plot(t,x1(:,1),'r-',t,x2(:,1),'k--',t,x3(:,1),'b:',t,x1(:,end),'r.-',t,x2(:,end),'k-.',t,x3(:,end),'b.--')xlabel('Time(h)')ylabel('x_1_1,x_1_2,x_1_3,x_N_1,x_N_2,x_N_3')title('塔顶和塔釜产品从进料开始直至
4、稳态的动态浓度曲线')legend('x_1_1','x_1_2','x_1_3','x_N_1','x_N_2','x_N_3')%稳态图figureplate=1:Nt;plot(plate,x1(end,:),'r.-',plate,x2(end,:),'k.--',plate,x3(end,:),'b.:')xlabel('塔板')ylabel('稳态时苯,甲苯,苯乙烯的组成')title('稳态时精馏塔的浓度曲线')legend('苯','甲苯','苯乙烯')%---------------------------------
5、---------------------------------functiondydt=DistMassBalances(t,y)%物料平衡方程组globalFz1z2z3Ralpha1alpha2alpha3M1MNMNtNfV1VDLL1Wx1=y(1:Nt);%组分1（苯）x2=y(Nt+1:2*Nt);%组分2（甲苯）x3=1-x1-x2;%组分3（苯乙烯）%气相平衡denom=alpha1*x1+alpha2*x2+alpha3*x3;y1=alpha1*x1./denom;y2=alpha2*x2./denom;%对塔顶
6、冷凝器（i=1）i=1;dx1dt(i)=(V*y1(i+1)-(L+D)*x1(i))/M1;dx2dt(i)=(V*y2(i+1)-(L+D)*x2(i))/M1;%精馏段（i=2~Nf-1）fori=2:Nf-1dx1dt(i)=(L*(x1(i-1)-x1(i))+V*(y1(i+1)-y1(i)))/M;dx2dt(i)=(L*(x2(i-1)-x2(i))+V*(y2(i+1)-y2(i)))/M;end%进料板（i=Nf）i=Nf;dx1dt(i)=(F*z1+L*x1(i-1)-L1*x1(i)+V1*y1(i+1)-V
7、*y1(i))/M;dx2dt(i)=(F*z2+L*x2(i-1)-L1*x2(i)+V1*y2(i+1)-V*y2(i))/M;%提馏段（Nf+1~Nt-1）fori=Nf+1:Nt-1dx1dt(i)=(L1*(x1(i-1)-x1(i))+V1*(y1(i+1)-y1(i)))/M;dx2dt(i)=(L1*(x2(i-1)-x2(i))+V1*(y2(i+1)-y2(i)))/M;end%塔釜（i=Nt）i=Nt;dx1dt(i)=(L1*x1(i-1)-V1*y1(i)-W*x1(i))/MN;dx2dt(i)=(L1*x2
8、(i-1)-V1*y2(i)-W*x2(i))/MN;dydt=[dx1dtdx2dt]';2.恒定滞液量的双组分间歇精馏的回流比控制functionBatchDistill%假设：%1.沸腾上升的蒸汽流量
