matlab作圖法計(jì)算精餾理論板數(shù)

1、MATLAB圖解法計(jì)算精餾塔理論板數(shù) 中文摘要：雙組份精餾是化工生產(chǎn)中重要的單元操作，運(yùn)用matlab【1】的強(qiáng)大功能，通過繪制精餾段和提留段操作線方程以及進(jìn)料q線方程，運(yùn)用圖解法計(jì)算理論板數(shù)和進(jìn)料板位置，使得求解精確、簡(jiǎn)潔。中文關(guān)鍵詞：matlab 精餾 理論板 圖解法Title：Graphical Method MATLAB number of theoretical distillationAbstract：Two-component distillation is an important chemical production unit operations, using the p

2、ower 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朗讀顯示對(duì)應(yīng)的拉丁字符的拼音 字典Keywords：Mat

3、lab Distillation Theoretical plate Graphic引 言 二元精餾塔的計(jì)算可用圖解法計(jì)算理論板數(shù)。通過函數(shù)求解曲線繪制，采用圖解法計(jì)算理論板數(shù)。在工程計(jì)算中matlab語(yǔ)言擁有大量的命令集和可用函數(shù)集可以完成各種計(jì)算和數(shù)據(jù)處理，集數(shù)值計(jì)算和圖形處理功能于一身，形式簡(jiǎn)單，易于掌握。1 問題敘述：常壓操作的的連續(xù)精餾塔，分離含二硫化碳0.53（摩爾分?jǐn)?shù)）的二硫化碳四氯化碳混合液，要求塔頂產(chǎn)品中含二硫化碳不低于0.986，塔底產(chǎn)品中含二硫化碳不高于0.0220.操作回流比為3.8，計(jì)算進(jìn)料液相分率為1.550時(shí)的理論板層數(shù)和加料板位置。常壓下的二硫化碳和四氯化碳的氣

4、液平衡數(shù)據(jù)【2】液相中二硫化碳摩爾分率x氣相中二硫化碳摩爾分率y液相中二硫化碳摩爾分率x氣相中二硫化碳摩爾分率y00.02960.06150.11060.14350.258000.08230.15550.26600.33250.49500.39080.53180.66300.75740.86041.00.63400.74700.82900.87900.93201.02 數(shù)學(xué)模型：精餾段操作線方程：提餾段操作線方程：相平衡方程：q線方程：3解決問題方法：1、列出已知條件2、用line繪出對(duì)角線，擬合數(shù)據(jù)用plot繪制平衡線3、利用fplot繪制精餾段操作線4、利用fplot繪制進(jìn)料線5、求3,4

5、兩線交點(diǎn)f，采用數(shù)值解法，調(diào)用fzero6、連接f與塔底組成點(diǎn)c繪制提餾段操作線7、繪制階梯，計(jì)算理論板數(shù)給定，R，F(xiàn)框圖繪制對(duì)角線，擬合氣液平衡數(shù)據(jù)繪制平衡線繪制精餾段操作線方程繪制進(jìn)料線方程聯(lián)立兩直線求交點(diǎn)繪制提餾段操作線繪制階梯，求理論板數(shù)具體程序【1，2】如下：%圖解法確定精餾塔理論板層數(shù)function fdistillation%圖解法確定理論板層數(shù)clear,clc%1.已知數(shù)據(jù)xD = 0.986;%塔頂組成xW = 0.0220%塔釜組成xF = 0.53;%進(jìn)料組成R = 3.8;%回流比q = 1.550;%進(jìn)料熱狀況vledata = %汽液平衡數(shù)據(jù)0 0.0296 0

6、.0615 0.1106 0.1435 0.2580 0.3908 0.5318 0.6630 0.7574 0.8604 1.00 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

7、 繪制對(duì)角線line(0.0,1.0,0.0,1.0)% 4 繪制精餾段操作線fplot(recope,0 xD,'r',R,xD)hold on % 5 繪制進(jìn)料線方程fplot(qline,xF 1,'g',q,xF)hold on % 6 與進(jìn)料線的交點(diǎn)xc = fzero(cross,0.5,optimset('fzero'),R,xD,q,xF);yc = recope(xc,R,xD);% 7 繪制提餾段操作線line(xc,xW,yc,xW,'Color','c')% 8 繪制梯級(jí)x0 = xD;y0

8、 = xD;n=0;while 1n = n+1;x0,y0 = tri(x0,y0,R ,xD,p ,xc ,yc,xW,n);if x0<xWbreak;endend% 9修飾s = sprintf('理論板數(shù)為%d',n);text(0.1,0.9,s)xlim(0,1),ylim(0,1),title('圖解法確定精餾塔理論板數(shù)'),xlabel('x'),ylabel('y')%-精餾段操作線-function y = recope(x,R,xD)y = R*x/(R+1)+xD/(R+1);%-進(jìn)料線方程-fun

9、ction y = qline(x,q,xF)y = q*x/(q-1)-xF/(q-1);%-交點(diǎn)-function y = cross(x,R,xD,q,xF)y = recope(x,R,xD)-qline(x,q,xF);%-平衡關(guān)系逆運(yùn)算-function y = vlei(x,p,yn)y = polyval(p,x)-yn;%-繪制梯級(jí)-function x,y = tri(x0,y0,R,xD,p,xc,yc,xW,n)x = fzero(vlei,x0,optimset('fzero'),p,y0);if x>xcy = recope(x,R,xD);e

10、lsey = yc + (yc-xW)*(x-xc)/(xc-xW);endline (x0,x,x,y0,y0,y,'Color','m'),text(x,y0,int2str(n),'VerticalAlignment','bottom')grid on結(jié)果如下：由圖可以看出需要12塊理論板結(jié)論：應(yīng)用 matlab求解二元精餾塔理論板數(shù)，方法簡(jiǎn)單實(shí)用，可以方便快捷的得到計(jì)算結(jié)果。參考文獻(xiàn)1 田文德，王曉紅.化工過程計(jì)算機(jī)基礎(chǔ).北京：化工過程計(jì)算機(jī)應(yīng)用基礎(chǔ).北京：化學(xué)工業(yè)出版社20072 陳敏恒，叢德滋，方圖南等.化工原理.第3

11、版.北京：化學(xué)工業(yè)出版社200600.01000000000000000.02000000000000000.04000000000000000.06000000000000000.08000000000000000.1000000000000000.1400000000000000.1800000000000000.2000000000000000.2500000000000000.3000000000000000.3500000000000000.4000000000000000.4500000000000000.5000000000000000.5500000000000000.6000

12、000000000000.6500000000000000.7000000000000000.7500000000000000.8000000000000000.8500000000000000.8940000000000000.9000000000000000.950000000000000100.1100000000000000.1750000000000000.2730000000000000.3400000000000000.3920000000000000.4300000000000000.4820000000000000.5130000000000000.5250000000000000.5510000000000000.5750000000000000.5950000000000000.6140000000000000.635000000000