甲醇—水汽液平衡.doc

。ropt1.01670 1.01772 1.02178 1.02687 1.03195 1.03703 1.04212 R/Rmin1.000 1.0011.0051.011.0151.021.025年总费用1344560.331329387.811318482.11314378131418813144461314831.5%2.37%1.22%0.39%0.07%0.06%0.08%0.11%ropt1.04720 1.05228 1.05533 1.05737 1.06754 1.07770 1.08787 R/Rmin1.031.0351.0381.041.051.061.07年总费用1316669131808413189261319446132355513271491331896%0.25%0.36%0.42%0.46%0.77%1.05%1.41%ropt1.09804 1.10820 1.11837 1.22004 1.32171 1.42338 1.52505 R/Rmin1.081.091.11.21.31.41.5年总费用1336073134051613484191399738145355015096671566271%1.73%2.06%2.67%6.57%10.67%14.94%19.25%ropt1.626721.728391.830061.931732.033400R/Rmin1.61.71.81.9200年总费用162348516806581741322179836518565011557936.51607201.2%36.56%41.37%46.48%51.27%56.16%31.05%35.19%附录二 甲醇水汽液平衡数据(摩尔组成)txytxy100.000.000.00075.300.400.72996.400.020.13473.100.500.77993.500.040.23471.200.600.82591.200.060.30469.300.700.87089.300.080.36567.600.800.91587.700.100.41866.000.900.95884.400.150.51765.000.950.97981.700.200.57964.501.001.00078.000.300.665附录一甲醇水系统的主要物理性质附录三 优化设计程序源代码优化程序定义全局变量Dim J1#, J2#, J3#, J4#, JJ#Dim N#, R#, Ropt#Dim lilunbanshu#, jinliaoweizhi%, tajing#, chukouwendu#, chuanremianji#, zongtagao#, tiliuduanbanshu#, jinliuduanbanshu#Dim XF#, F#, q#, XD#, D#, td#, rD#, po#, u#, Rmin#, t1#, Cw#, Cp#, SI#, HETP#Dim Co#, HA#, f1#, f2#, a#, b#, FL#, #, #, bo#, Fc#优化所需参数Public Sub Form_Load()XF = 0.3151: F = 402.34: q = 1XD = 0.982: D = 128.97: td = 64.93: rD = 35373.48: po = 101.3u = 5.4464: Rmin = 1.0167t1 = 20: Cw = 0.0002: Cp = 4.1875: Co = 0.03: cpa = 15674.4HETP = 0.462: HA = 6f1 = 1: f2 = 6.5: a = 487: b = 0.72: SI = 3.73FL = 6.22: = 7200: = 7860: bo = 0.005: Fc = 0.125Text1.Text = 402.34Text2.Text = 0.3151Text3.Text = 128.97Text4.Text = 0.982Text5.Text = 35373.48Text6.Text = 64.93Text7.Text = 1Text8.Text = 1.0167Text9.Text = 7200Text10.Text = 3.73Text11.Text = 0.125Text12.Text = 6.22Text13.Text = 0.005Text14.Text = 7860Text15.Text = 5.4464Text16.Text = 0.462Text17.Text = 6Text18.Text = 15674.4Text19.Text = 0.0002Text20.Text = 4.1875Text21.Text = 20Text22.Text = 2000Text23.Text = 1Text24.Text = 6.5Text25.Text = 487Text26.Text = 0.72Text27.Text = 0.03Text28.Text = 1.01Text29.Text = 2Text30.Text = 0.0001Text31.Text = Text32.Text = Text33.Text = Text34.Text = Text35.Text = Text36.Text = Text37.Text = Text38.Text = Text39.Text = Text40.Text = Text41.Text = Text42.Text = Text43.Text = Text44.Text = Text45.Text = Text46.Text = Text47.Text = End Sub主程序Private Sub Command1_Click() 菲波拿契法求RoptDim Aa#, Bb#, W#(1 To 50), i%, K%, N#, M%, R1#, R2#, #Dim JJ1#, JJ2#Aa = 1.01 * Rmin: Bb = 2 * Rmin 搜索区间Aa,BbW(1) = 1: W(2) = 2: W(3) = 3: i = 1: = 0.0001Do While W(i + 2) = (Bb - Aa) / ) i = i + 1 W(i + 2) = W(i) + W(i + 1)Loop R1 = Aa + (Bb - Aa) * W(i) / W(i + 2): JJ1 = j(R1)N = i + 2: K = 1: M = 0Do While K N - 1 If M = 0 Then R2 = Aa + (Bb - Aa) * W(N - K) / W(N - K + 1) JJ2 = j(R2) Else R1 = Aa + (Bb - Aa) * W(N - K - 1) / W(N - K + 1) JJ1 = j(R1) End If If JJ1 JJ2 Then Bb = R2: R2 = R1: JJ2 = JJ1: M = 1 Else Aa = R1: R1 = R2: JJ1 = JJ2: M = 0 End If K = K + 1LoopR = (Aa + Bb) / 2Ropt = RJJ = j(R)Text31.Text = RoptText32.Text = RminText45.Text = Ropt / RminText33.Text = lilunbanshuText34.Text = zongtagaoText40.Text = J1Text41.Text = J2Text42.Text = J3Text43.Text = J4Text44.Text = JJText37.Text = tajingText38.Text = chukouwenduText39.Text = chuanremianjiText46.Text = Ropt * DText47.Text = (Ropt + 1) * DText35.Text = tiliuduanbanshu * HETPText36.Text = jinliuduanbanshu * HETPEnd SubJ函数Public Function j(R#) As DoubleCall jjj1(R#, J1#)Call jjj2(R#, J2#)Call jjj3(R#, J3#)Call jjj4(R#, J4#)j = J1 + J2 + J3 + J4End Function求J1Public Sub jjj1(R#, J1#)Dim DT#, H#, Ws#, CH#Call tabanshu(R#, N#)DT = Sqr(R + 1) * D * 22.4 / (3600 * 0.785 * u) * (273 + td) / 273 * 101.3 / po)H = N * HETP + HAWs = 3.14 * DT * (H + 0.8116 * DT) * bo * 为碳钢的密度CH = FL * Exp(6.95 + 0.1808 * Log(Ws) + 0.02468 * (Log(Ws) 2 + 0.0158 * H / DT)J1 = SI * (Fc + 0.06) * CHtajing = DTzongtagao = HEnd Sub求J2Public Sub jjj2(R#, J2#)Dim xx1#, xx0#, CD#, ff#, df#, t2#, AD#, KD#KD = 2000: xx1 = 70Do 牛顿迭代法求冷却水最佳出口温度t2 xx0 = xx1 CD = 1.3 * SI * a * b * f1 * f2 * Fc * (R + 1) * D * rD / (td - t1) (b - 1) / KD b ff = -Cw * / Cp + CD * (xx0 - 1) / xx0 / Log(xx0) (1 - b) * (xx0 - 1 - Log(xx0) df = CD * (xx0 - 1) / xx0 / Log(xx0) (2 - b) * (b - 1) * (xx0 - 1 - Log(xx0) 2 / (xx0 - 1) 2 + Log(xx0) xx1 = xx0 - ff / dfLoop Until Abs(xx1 - xx0) XF Then Y(i + 1) = R * X(i) / (R + 1) + XD / (R + 1) n1 = i + 1 + (X(i) - XF) / (X(i) - X(i + 1) Else Y(i + 1) = W * (X(i) - XW) / V If X(i) XW Then Exit Do End If i = i + 1 xx(i) = (Y(i) / (3.3874 * (1 - Y(i) (1 / 0.7977) X(i) = xx(i) / (1 + xx(i) Loop N = i - 1 + (X(i - 1) - XW) / (X(i - 1) - X(i) lilunbanshu = N tiliuduanbanshu = n1 jinliuduanbanshu = N - n1End Sub调整ROPT程序：定义全局变量Dim J1#, J2#, J3#, J4#, JJ#Dim N#, R#, Ropt#Dim lilunbanshu#, jinliaoweizhi%, tajing#, chukouwendu#, chuanremianji#, zongtagao#, tiliuduanbanshu#, jinliuduanbanshu#Dim XF#, F#, q#, XD#, D#, td#, rD#, po#, u#, Rmin#, t1#, Cw#, Cp#, SI#, HETP#Dim Co#, HA#, f1#, f2#, a#, b#, FL#, #, #, bo#, Fc#优化所需参数Public Sub Form_Load()XF = 0.3151: F = 402.34: q = 1XD = 0.982: D = 128.97: td = 64.93: rD = 35373.48: po = 101.3u = 5.4464: Rmin = 1.0167t1 = 20: Cw = 0.0002: Cp = 4.1875: Co = 0.03: cpa = 15674.4HETP = 0.462: HA = 6f1 = 1: f2 = 6.5: a = 487: b = 0.72: SI = 3.73FL = 6.22: = 7200: = 7860: bo = 0.005: Fc = 0.125Text1.Text = 402.34Text2.Text = 0.3151Text3.Text = 128.97Text4.Text = 0.982Text5.Text = 35373.48Text6.Text = 64.93Text7.Text = 1Text8.Text = 1.0167Text9.Text = 7200Text10.Text = 3.73Text11.Text = 0.125Text12.Text = 6.22Text13.Text = 0.005Text14.Text = 7860Text15.Text = 5.4464Text16.Text = 0.462Text17.Text = 6Text18.Text = 15674.4Text19.Text = 0.0002Text20.Text = 4.1875Text21.Text = 20Text22.Text = 2000Text23.Text = 1Text24.Text = 6.5Text25.Text = 487Text26.Text = 0.72Text27.Text = 0.03Text28.Text = 1.01Text29.Text = 2Text30.Text = 0.0001Text31.Text = Text32.Text = Text33.Text = Text34.Text = Text35.Text = Text36.Text = Text37.Text = Text38.Text = Text39.Text = Text40.Text = Text41.Text = Text42.Text = Text43.Text = Text44.Text = Text45.Text = Text46.Text = Text47.Text = End Sub主程序Private Sub Command1_Click() 菲波拿契法求RoptR = Text31.TextRopt = RJJ = j(R)Text32.Text = RminText45.Text = Ropt / RminText33.Text = lilunbanshuText34.Text = zongtagaoText40.Text = J1Text41.Text = J2Text42.Text = J3Text43.Text = J4Text44.Text = JJText37.Text = tajingText38.Text = chukouwenduText39.Text = chuanremianjiText46.Text = Ropt * DText47.Text = (Ropt + 1) * DText35.Text = tiliuduanbanshu * HETPText36.Text = jinliuduanbanshu * HETPEnd SubJ函数Public Function j(R#) As DoubleCall jjj1(R#, J1#)Call jjj2(R#, J2#)Call jjj3(R#, J3#)Call jjj4(R#, J4#)j = J1 + J2 + J3 + J4End Function求J1Public Sub jjj1(R#, J1#)Dim DT#, H#, Ws#, CH#Call tabanshu(R#, N#)DT = Sqr(R + 1) * D * 22.4 / (3600 * 0.785 * u) * (273 + td) / 273 * 101.3 / po)If DT 1 ThenDT = Int(DT * 10 + 1) / 10ElseDT = Int(DT * 5 + 1) * 0.2End IfH = N * HETP + HAWs = 3.14 * DT * (H + 0.8116 * DT) * bo * 为碳钢的密度CH = FL * Exp(6.95 + 0.1808 * Log(Ws) + 0.02468 * (Log(Ws) 2 + 0.0158 * H / DT)J1 = SI * (Fc + 0.06) * CHtajing = DTzongtagao = HEnd Sub求J2Public Sub jjj2(R#, J2#)Dim xx1#, xx0#, CD#, ff#, df#, t2#, AD#, KD#KD = 2000: xx1 = 70Do 牛顿迭代法求冷却水最佳出口温度t2 xx0 = xx1 CD = 1.3 * SI * a * b * f1 * f2 * Fc * (R + 1) * D * rD / (td - t1) (b - 1) / KD b ff = -Cw * / Cp + CD * (xx0 - 1) / xx0 / Log(xx0) (1 - b) * (xx0 - 1 - Log(xx0) df = CD * (xx0 - 1) / xx0 / Log(xx0) (2 - b) * (b - 1) * (xx0 - 1 - Log(xx0) 2 / (xx0 - 1) 2 + Log(xx0) xx1 = xx0 - ff / dfLoop Until Abs(xx1 - xx0) XF Then Y(i + 1) = R * X(i) / (R + 1) + XD / (R + 1) n1 = i + 1 + (X(i) - XF) / (X(i) - X(i + 1) Else Y(i + 1) = W * (X(i) - XW) / V If X(i) XW Then Exit Do End If i = i + 1 xx(i) = (Y(i) / (3.3874 * (1 - Y(i) (1 / 0.7977) X(i) = xx(i) / (1 + xx(i) Loop N = i - 1 + (X(i - 1) - XW) / (X(i - 1) - X(i) lilunbanshu = N tiliuduanbanshu = n1 jinliuduanbanshu = N - n1End Sub目 录1 前言-12 方案论证2.1 精馏塔类型-12.2 精馏压力-12.3 进料方式-12.4 填料类型-22.5 加热方式-22.6 塔材料类型-23 数学模型的建立3.1 精馏塔塔体年投资折旧费及维修费用-33.2 冷凝器年运转费用-43.3 直接蒸汽加热费用-53.4 填料年折旧费- -54 数学模型的求解4.1 数学模型决策变量分析-54.2 主要工艺参数的求解-54.2.1 塔径的计算-54.2.2 塔板数的计算-6 4.2.2.1 相平衡关系的表示-6 4.2.2.2 N的计算-64.2.3 冷凝器年运转费用的计算- - -74.2.3.1 冷却水用量及冷凝器传热面积的计算- -74.2.3.2 冷凝器冷却水最佳出口温度的确定-74.2.4 直接加热蒸气费用的计算-84.3 数学模型的求解- -84.3.1 单变量最优化方法- -84.3.2 优化设计程序框图- -8 4.3.2.1 函数调用关系-9 5 优化设计计算5.1 数据预处理-10 5.1.1 进塔物料的计算-10 5.1.2 塔顶蒸气温度的计算-10 5.1.3 等板高度的计算-10 5.1.4 产品汽化潜热的计算-11 5.1.5 最