计算机辅助简捷法求理论版层数.doc

某大学化工原理课程论文 题 目： 计算机辅助简捷法求理论版层数 院 （系）： 资源与环境工程学院 班 级： 化工09-1 姓 名： 张* 学 号： 09号 指导教师： 熊* 教师职称： 教授 目录计算机辅助计算简捷法求理论版层数3摘要：3前言：4设计方法与步骤：5计算5计算最小回流比：7计算理论塔板数N：10界面辅助设计11结语：11参考文献12计算机辅助计算简捷法求理论版层数基于 Gilliland图与Liddle回归方程的简捷法求理论版层数的计算机辅助计算程序的设计摘要： 本文用VB程序实现了简捷法求理论版层数问题，基于Liddle的Gilliland图的回归方程，体现了计算机在实际化工设计领域的应用。 关键词： VB，简捷法求理论版层数，吉利兰/Gilliland，李德/Liddle，计算机辅助计算，化工原理。Abstract:This paper use VB program to achieve a simple method for the number of theory plate, which is based on Liddles regression equation (regression of Gilliland graph). This program reflects the computer applying in the field of the chemical engineering design.Key words:VB, simple method for the version of the theory of plies, Gilliland, Liddle, computer aided design, principles of chemical engineering.前言： 化工设计中精馏塔的理论塔板数的计算是一项繁琐而又重要的工作，即使是运用简捷法计算理论塔板数，也有很大的计算工作量。随着计算机技术突飞猛进的发展，化工行业中用计算机来辅助设计必将能大大减少设计人员的计算量，提高工作效率。在非计算机专业的本科教育中，大部分工科学生都以vb为基础计算机编程语言，所以本文以vb为编程语言，结合李德回归方程设计出简捷法求理论版层数的程序。设计方法与步骤：李德根据吉利兰图得到的回归方程：Y=0.545827-0.591422X+0.002743/X其中 ; 其中R为回流比，N为理论塔板数，为最小回流比，为最小理论塔板数。所以与为重要计算步骤。所以将vb程序分为3部分，也就是三个计算步骤。设计vb操作界面如下：程序运行界面计算计算芬斯克方程来计算： （不包括再沸器）vb可识别的线性计算公式为： = (Log(Xd / (1 - Xd) * (1 - Xw) / Xw) / (Log(Am) 1如此，将，输入到相应文本框中并赋值给相应变量，点击按钮即可得出 的值。计算最小理论踏板数的主要代码如下：Private Sub Command1_Click() If Text1.Text = Or Text2.Text = Or Text1.Text = Then If Text1.Text = Then MsgBox 请输入塔顶产品易挥发组分含量 Xd 的值 ElseIf Text3.Text = 请输入am值 Then MsgBox 请输入平均挥发度am值 ElseIf Text2.Text = Then MsgBox 请输入塔底产品易挥发组分含量 Xw 的值 End If Else Xd = CDbl(Text1.Text) Xw = CDbl(Text2.Text) Am = CDbl(Text3.Text) Ntem = (Log(Xd / (1 - Xd) * (1 - Xw) / Xw) / (Log(Am) - 1 Nmin = Int(Ntem) + 1 不包括再沸器 Text4.Text = Ntem 即为Nmin End IfEnd Sub计算计算最小回流比：同理，将最小回流比的解析法计算公式：转换为线性形式：Rmin = (1 / (Am - 1) * (Xd / Xq - Am * (1 - Xd) / (1 - Xq)其中，值的算法分两种，用option控件来选择计算方式：1，已知，直接输入2，由q线方程和气液平衡方程计算Xq,Yq的值其中，第二种方法计算时用此公式： = (-BB + (BB 2 - 4 * AA * CC) 0.5) / (2 * AA)进行计算。AA、BB、CC分别为一元二次方程的系数，他们的值由q线方程和气液平衡方程联立消去y后得到的一元二次方程得到，分别赋值计算，从而简化输入：AA = q * (a - 1) / (a - 1)BB = (q / (q - 1) - a - Xf)CC = -Xf / (q - 1)计算方法赋给按钮的点击事件，代码如下：Private Sub Command5_Click() If Text11.Text = Or Text5.Text = Then MsgBox 请输入Xf和进料热状况q 的值 Else a = Text3.Text q = Text5.Text Xf = Text11.Text AA = q * (a - 1) / (a - 1) BB = (q / (q - 1) - a - Xf) CC = -Xf / (q - 1) Xq = (-BB + (BB 2 - 4 * AA * CC) 0.5) / (2 * AA) Yq = a * Xq / (1 + (a - 1) * Xf) Text6.Text = Xq Text7.Text = Yq Command3.Enabled = True End IfEnd Sub 得到了Xq,Yq的值，紧接着，计算最小回流比：方法类似，代码如下：Private Sub Command3_Click() If Text6.Text = Or Text7.Text = Then If Text6.Text = Then MsgBox 请输入或计算 Xq 的值 ElseIf Text7.Text = 请输入am值 Then MsgBox 请输入Yq值 End If Else Xq = Text6.Text Yq = Text7.Text Rmin = (1 / (Am - 1) * (Xd / Xq - Am * (1 - Xd) / (1 - Xq) 数字的自然对数的运算vb自身就有。 Text8.Text = Rmin Command4.Enabled = True End IfEnd Sub 计算计算Xq，Yq与最小回流比Rmin计算理论塔板数N：计算理论塔板数需要知道回流比R，可输入R与最小回流比Rmin的比之，然后计算得出最终结果：代码和界面如下：Private Sub Command4_Click()If Text9.Text = ThenMsgBox 请输入回流比为最小回流比的倍数Elser = Rmin * Text9.TextX = (r - Rmin) / (r + 1)Y = 0.545827 - 0.591422 * X + 0.002743 / XN = Int(2 * Y + Ntem) / (1 - Y) + 1Text10.Text = NEnd ifEnd Sub计算理论板层数N 界面辅助设计 接下来，进行软件的菜单的简单编辑以及帮助提示等信息的加入。 菜单的文件选项下简单加入清除数据和退出两个菜单选项，文件中添加使用说明和联系方式等信息。如图：菜单-文件菜单-帮助结语： 实际上，本文原本想做成一键式计算程序，但是考虑到教学等方面的方便，将整体计算分为了多步，将关键量的计算分别进行计算，体现出了计算过程，便于理解简洁法求最少理论塔板数的方法。本论文是在我的导师熊楚安老师的亲切关怀和悉心指导下完成的。他严肃的科学态度，严谨的治学精神，精益求精的工作作风，深深地感染和激励着我。熊老师不仅在学业上给我以精心指导，同时还在思想、生活上给