甲醇水汽液平衡

ropt 1.01670 1.01772 1.02178 1.02687 1.03195 1.03703 1.04212 R/Rmin 1.000 1.001 1.005 1.01 1.015 1.02 1.025

年总费

用1344560.33 1329387.8

1

1318482

.1

1314831.

5

% 2.37% 1.22% 0.39% 0.07% 0.06% 0.08% 0.11% ropt 1.04720 1.05228 1.05533 1.05737 1.06754 1.07770 1.08787 R/Rmin 1.03 1.035 1.038 1.04 1.05 1.06 1.07

年总费

用

% 0.25% 0.36% 0.42% 0.46% 0.77% 1.05% 1.41% ropt 1.09804 1.10820 1.11837 1.22004 1.32171 1.42338 1.52505 R/Rmin 1.08 1.09 1.1 1.2 1.3 1.4 1.5

年总费

用

% 1.73% 2.06% 2.67% 6.57% 10.67% 14.94% 19.25% ropt 1.62672 1.72839 1.83006 1.93173 2.0334 0 0

R/Rmin 1.6 1.7 1.8 1.9 2 0 0

年总费用1557936.

5

1607201.

2

% 36.56% 41.37% 46.48% 51.27% 56.16% 31.05% 35.19%

附录二甲醇—水汽液平衡数据(摩尔组成)

t x y t x y

100.00 0.00 0.000 75.30 0.40 0.729 96.40 0.02 0.134 73.10 0.50 0.779 93.50 0.04 0.234 71.20 0.60 0.825 91.20 0.06 0.304 69.30 0.70 0.870 89.30 0.08 0.365 67.60 0.80 0.915 87.70 0.10 0.418 66.00 0.90 0.958 84.40 0.15 0.517 65.00 0.95 0.979 81.70 0.20 0.579 64.50 1.00 1.000 78.00 0.30 0.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 = 1

XD = 0.982: D = 128.97: td = 64.93: rD = 35373.48: po = 101.3

u = 5.4464: Rmin = 1.0167

t1 = 20: Cw = 0.0002: Cp = 4.1875: Co = 0.03: cpa = 15674.4

HETP = 0.462: HA = 6

f1 = 1: f2 = 6.5: a = 487: b = 0.72: SI = 3.73

FL = 6.22: θ= 7200: ρ= 7860: bo = 0.005: Fc = 0.125

Text1.Text = 402.34

Text2.Text = 0.3151

Text3.Text = 128.97

Text4.Text = 0.982

Text5.Text = 35373.48

Text6.Text = 64.93

Text7.Text = 1

Text8.Text = 1.0167

Text9.Text = 7200

Text10.Text = 3.73

Text11.Text = 0.125

Text12.Text = 6.22

Text13.Text = 0.005

Text14.Text = 7860

Text15.Text = 5.4464

Text16.Text = 0.462

Text17.Text = 6

Text18.Text = 15674.4

Text19.Text = 0.0002

Text20.Text = 4.1875

Text21.Text = 20

Text22.Text = 2000

Text23.Text = 1

Text24.Text = 6.5

Text25.Text = 487

Text26.Text = 0.72

Text27.Text = 0.03

Text28.Text = 1.01

Text29.Text = 2

Text30.Text = 0.0001

Text31.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() '菲波拿契法求Ropt

Dim 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.0001

Do While W(i + 2) <= ((Bb - Aa) / ε)

i = i + 1