开发工具:
文件大小: 331kb
下载次数: 0
上传时间: 2011-03-18
详细说明： 自控工程课程设计 程序： #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]';} ...展开收缩
(系统自动生成,下载前可以参看下载内容)
下载文件列表
相关说明
- 本站资源为会员上传分享交流与学习,如有侵犯您的权益,请联系我们删除.
- 本站是交换下载平台，提供交流渠道，下载内容来自于网络，除下载问题外，其它问题请自行百度。
- 本站已设置防盗链，请勿用迅雷、QQ旋风等多线程下载软件下载资源，下载后用WinRAR最新版进行解压.
- 如果您发现内容无法下载，请稍后再次尝试；或者到消费记录里找到下载记录反馈给我们.
- 下载后发现下载的内容跟说明不相乎，请到消费记录里找到下载记录反馈给我们，经确认后退回积分.
- 如下载前有疑问，可以通过点击"提供者"的名字，查看对方的联系方式，联系对方咨询.
