matlab作图法计算精馏理论板数.docVIP

MATLAB图解法计算精馏塔理论板数 中文摘要：双组份精馏是化工生产中重要的单元操作，运用matlab【1】的强大功能，通过绘制精馏段和提留段操作线方程以及进料q线方程，运用图解法计算理论板数和进料板位置，使得求解精确、简洁。 中文关键词：matlab 精馏 理论板 图解法 Title：Graphical Method MATLAB number of theoretical distillation Abstract：Two-component distillation is an important chemical production unit operations, using the power of matlab, by drawing the rectifying section operating line segment and the retention equation and the feed q line equation, the use of graphical method to calculate the theoretical plate number and location of the feed plate Makes solving the accurate, concise 朗读 显示对应的拉丁字符的拼音 ? 字典 Keywords：Matlab Distillation Theoretical plate Graphic 引 言 二元精馏塔的计算可用图解法计算理论板数。通过函数求解曲线绘制，采用图解法计算理论板数。在工程计算中matlab语言拥有大量的命令集和可用函数集可以完成各种计算和数据处理，集数值计算和图形处理功能于一身，形式简单，易于掌握。 1 问题叙述：常压操作的的连续精馏塔，分离含二硫化碳0.53（摩尔分数）的二硫化碳—四氯化碳混合液，要求塔顶产品中含二硫化碳不低于0.986，塔底产品中含二硫化碳不高于0.0220.操作回流比为3.8，计算进料液相分率为1.550时的理论板层数和加料板位置。 常压下的二硫化碳和四氯化碳的气液平衡数据【2】 液相中二硫化碳摩尔分率x 气相中二硫化碳摩尔分率y 液相中二硫化碳摩尔分率x 气相中二硫化碳摩尔分率y 0 0.0296 0.0615 0.1106 0.1435 0.2580 0 0.0823 0.1555 0.2660 0.3325 0.4950 0.3908 0.5318 0.6630 0.7574 0.8604 1.0 0.6340 0.7470 0.8290 0.8790 0.9320 1.0 2 数学模型： 精馏段操作线方程： 提馏段操作线方程： 相平衡方程： q线方程： 3解决问题方法： 1、列出已知条件 2、用line绘出对角线，拟合数据用plot绘制平衡线 3、利用fplot绘制精馏段操作线 4、利用fplot绘制进料线 5、求3,4两线交点f，采用数值解法，调用fzero 6、连接f与塔底组成点c绘制提馏段操作线 7、绘制阶梯，计算理论板数 框图 具体程序【1，2】如下： %图解法确定精馏塔理论板层数 function fdistillation %图解法确定理论板层数 clear,clc %1.已知数据 xD = 0.986;%塔顶组成 xW = 0.0220%塔釜组成 xF = 0.53;%进料组成 R = 3.8;%回流比 q = 1.550;%进料热状况 vledata = [%汽液平衡数据 0 0.0296 0.0615 0.1106 0.1435 0.2580 0.3908 0.5318 0.6630 0.7574 0.8604 1.0 0 0.0823 0.1555 0.2660 0.3325 0.4950 0.6340 0.7470 0.8290 0.8790 0.9320 1.0]; %2 绘制平衡线 plot(vledata(1,:),vledata(2,:),*) hold on p = polyfit(vledata(1,:),vledata(2,:),3); x1 = [0.0 :0.01:1.0]; y1 = polyval(p,x1); plot(x1,y1) hold on %3 绘制对角线 line([0.0,1.0],[0.0,1.0]) % 4 绘制精馏段操作线 fplot(@recope,[0 xD],[],[],r,R,xD) hold on % 5 绘制进料线方程 fplot(@qline,[