文件名称：自控工程减压塔课程设计
文件大小：331KB
文件格式：DOC
更新时间：2014-04-07 14:15:57
减压塔
自控工程课程设计 程序： #include“stdio.h” main() { function y1 clear all clc global F z1 z2 R alpha1 alpha2 M1 MN M Nt Nf V1 V D L L1 W F=25;%进料流量,kmol/h R=5;%回流比 z1=0.4;%子塔段一塔顶油分的进料组成（摩尔分率） z2=0.6;%减二线油分的进料组成（摩尔分率） M1=75;%冷凝器的塔板滞液量(kmol) M=10;%塔板的滞液量(kmol) MN=150;%塔底的滞液量(kmol) q=1;%饱和进料 tf=10; dt=1; %相对挥发度 alpha1=2.49; alpha2=1; Nt=8;%塔板总数 Nf=2;%进料位置 V1=1.5;%塔底进入蒸汽流量(kmol/h) %精馏段 V=V1-(1-q)*F; D=V/(R+1); L=V-D; %提馏段 L1=L+F; W=L1-V1; x1=z1*ones(1,Nt); x2=z2*ones(1,Nt); [t,y]=ode45(@DistMassBalances,[0:dt:tf],[x1 x2]) %输出结果 x1=y(:,1:Nt);%子塔段四塔顶油分的液相组成（摩尔分率） x2=y(:,Nt+1:2*Nt);%减二线油分的液相组成（摩尔分率） plot(t,x1(:,1),'r-') xlabel('Time(h)') title('减二线油组分从进料开始直至稳态的动态浓度曲线') %稳态图 figure plate=1:Nt; plot(plate,x1(end,:),'r.-') xlabel('塔板') ylabel('稳态时减二线油的组成') title('稳态时子塔段四的浓度曲线') legend('减二线油') function dydt=DistMassBalances(t,y)%物料平衡方程组 global F z1 z2 R alpha1 alpha2 M1 MN M Nt Nf V1 V D L L1 W x1=y(1:Nt);%组分一（子塔段四塔顶油） x2=y(Nt+1:2*Nt);%组分二（减二线油） %气相平衡 denom=alpha1*x1+alpha2*x2; y1=alpha1*x1./denom; y2=alpha2*x2./denom; %对塔顶冷凝器（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;哈尔滨工业大学工学硕士学位论文 -46- %精馏段（i=2~Nf-1） for i=2:Nf-1 dx1dt(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*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） for i=Nf+1:Nt-1 dx1dt(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(i-1)-V1*y2(i)-W*x2(i))/MN; dydt=[dx1dt dx2dt]';}