基于Excel精馏塔理论塔板数图解法.docVIP

基于Excel的精馏塔理论塔板数的图解法 1引言 精馏是一种重要的传质单元操作，精馏塔理论塔板数的计算是其重要的内容, 无论是对设计过程还是操作过程都有重要的意义。理论塔板数的计算方法有逐板计算法和图解法。其中逐板计算法是计算理论塔板数的基本方法, 结果准确、清晰, 不仅可以计算出塔板数和进料板位置, 而且可以得到每一块塔板上的气液组成, 但此法计算过程繁琐, 工作量大。用图解法比较直观、简便，但手工图解法准确性较差。用Excel的图表功能, 可以比较方便地解决这一问题。本文给出用图解法求精馏塔理论塔板数的通用程序, 程序采用VBA编制，该程序具有较好的通用性，对于图解法求精馏塔理论塔板数的不同问题, 只需在表格中输入体系的气液平衡数据以及分离要求，即可计算出精馏塔理论塔板数，加料板位置及每一块理论塔板上的汽液组成，并且计算结果是以图形的方式表示，计算结果较为直观。 2原理和步骤 2.1有关方程 (1)气液相平衡方程[1] y=f(x) （1） 在计算过程中主要是由气相的组成来计算相应板上的液相组成，上述方程应变换为以气相组成来表示液相组成的方程。 x=g(y) （2） 通常情况下已知x-y气液平衡数据，这时，对于给定的y，可用一元三点插值公式来求取x的数值。 （3） 其中xk, xk+1, xk+2和yk, yk+1, yk+2分别为最接近y的三个插值结点和相应的函数值，通过程序中的函数即可采用插值法进行计算。 (2)精馏段操作线方程 （4） (3)提馏段操作线方程 （5） (4)q线方程 （6） 2.2有关计算和图表绘制 根据分离要求(塔顶、塔底和进料组成xD，xW，xF)、进料状态q值和回流比R值，进行以下计算： (1)根据精馏段操作线方程和q线方程，求出两条线的交点d，坐标为(xd，yd) （7） （8） (2)根据d点(xd，yd)和D点(xD，xD)求出精馏段操作线方程 y=(xD-yd)/(xD-xd)*(x-xd)+yd （9） (3) 根据d点(xd，yd)和W点(xW，xW)求出提馏段操作线方程 y=(xW-yd)/(xW -xd)*(x-xd)+yd （10） (4)求出精馏段操作线在y轴上的截距b b=xD/(R+1) （11） (5)连接b点(0，b)和D点(xD，xD)画出精馏段操作线 (6)连接F点(xF，xF)和d点(xd，yd)画出q线 (7)连接d点(xd，xd)和W点(xW，xW)画出提馏段操作线 (8)根据平衡线，精馏段操作线和提馏段操作线求解每层塔板上的气液平衡数据，进一步求解出精馏段塔板数和提馏段塔板数。 3 VBA程序 精馏塔理论板数的求解过程采用表格中的宏进行计算，将宏与表格中的图形进行关联，用鼠标单击图形即可调用宏进行计算，计算过程中需要的原始数据通过表格输入，计算结果输出到表格中，然后进一步画出图形。 表中CDEFNOP七列为原始数据，GHIPKLAB八列为计算出的数据，表格中C，D两列数据为对角线，E，F两列数据为平衡线，N，O, P 三列数据为操作条件，分离要求等原始数据和计算出的d点坐标，G，H两列数据为精馏段操作线，I，J两列数据为q线以及提馏段操作线，K，L两列数据为阶梯线，A，B两列数据为每层塔板上的气液平衡组成。 表格中的宏采用Visual Basic 编辑器进行编制。此法制成的应用程序操作简便，通用性强，并且把复杂的手工操作交给程序来完成[2,3]，计算程序如下所示： Sub hong() Dim xcz(100) As Variant Dim ycz(100) As Variant Dim xph(100) As Variant Dim yph(100) As Variant Dim u(100) As Variant Dim v(100) As Variant For i=1 To 100 u(i)=ActiveSheet.Cells(i, 6) v(i)=ActiveSheet.Cells(i, 5) If u(1)=1 And u(i)=0 Then nph=i If u(1)=0 And u(i)=1 T