甲醇—水汽液平衡.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.11314831.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年总费用%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年总费用%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年总费用1557936.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,Bb]W(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 Sub'J函数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。