分离上机
- 格式：doc
- 大小：62.50 KB
- 文档页数：4
分离工程上机实验报告
一、实验目的：
通过上机编程熟悉等温闪蒸过程计算，通过上机编程熟悉多组分
精馏塔的简捷算法的过程。
二、实验内容：
1.试问苯（A）和正庚烷（B）混合物（系统组成zA=0.5）在333.15K和42500Pa
条件下，是否存在闪蒸现象，如果存在闪蒸，计算汽液相组成和汽化比。
已知 正庚烷的安托因常数为
A=20.7665 B=2911.32K C=-56.51K
苯的安托因常数为
A=20.7936 B=2788.15K C=-52.36K。
分别应用范拉方程和马克勒斯方程计算活度系数
安托因方程：㏑Pis=A-B/（T+C）
范拉方程：㏑Y1=A/(1+AX1/BX2)
2
㏑Y2=B/(1+BX2/AX1)2
(式中A=0.31559，B=0.402295)
马克勒斯方程：㏑Y1=X22[A+(B-A)X1]
㏑Y2=X12[B+(A-B)X2]
(式中A=0.312251，B=0.398646)
2：多组分精馏塔的简捷算法
苯（B）－甲苯（T）－二甲苯（X）－异丙苯（C）的混合物送入精馏塔分离，
进料组成为：ZB=0.2, ZT=0.3, ZX=0.1, ZC=0.4（摩尔分数）。相对挥发度数据：
B=2.25,T=1.00,X=0.33,C =0.21。分离要求：馏出液中异丙苯不大
于0.15％;釜液中甲苯不大于0.3％（摩尔分数）。计算最小理论板和全回流下的
物料分配。
三、实验报告：
第一题：
解：
编码：Private Sub Command1_Click()
Dim A1 As Single, B1 As Single, C1 As Single, A2 As Single, B2 As Single, C2 As Single,
A As Single, B As Single
Dim p As Single, p1 As Single, p2 As Single
Dim r1 As Single, r2 As Single
Dim z1 As Single, z2 As Single
Dim g1 As Single, g2 As Single
A = 0.315559: B = 0.402295
z1 = 0.5: z2 = 0.5
p = 42500: A1 = 20.7665: B1 = 2911.32: C1 = -56.51: A2 = 20.7936: B2 = 2788.15: C2
= -52.36
p1 = Exp(A1 - (B1 / 333.15 + C1))
p2 = Exp(A2 - (B2 / 333.15 + C2))
r1 = Exp(A / (1 + A * z1 / B * z2) ^ 2)
r2 = Exp(B / (1 + B * z2 / A * z1) ^ 2)
k1 = r1 * p1 / p
k2 = r2 * p2 / p
g1 = k1 * z1 + k2 * z2: g2 = z1 / k1 + z2 / k2
If g1 > 1 And g2 > 1 Then
lable11.Caption = "闪蒸问题成立"
Else
Label1.Caption = "闪蒸问题不成立"
End If
If g1 > 1 And g2 > 1 Then
Do
Ψ = Ψ1
γ1 = Exp(A / (1 + A * X1 / B * X2) ^ 2)
γ2 = Exp(B / (1 + B * X2 / A * X1) ^ 2)
K1 = γ1 * p1 / p
K2 = γ2 * p2 / p
f = (K1 - 1) * z1 / (1 + Ψ * (K1 - 1)) + (K2 - 1) * z2 / (1 + Ψ * (K2 -
1))
g = -((K1 - 1) ^ 2 * z1 / (1 + Ψ * (K1 - 1)) ^ 2 + (K2 - 1) ^ 2 * z2 / (1
+ Ψ * (K2 - 1)) ^ 2)
X1 = z1 / (1 + Ψ * (K1 - 1)): X2 = z2 / (1 + Ψ * (K2 - 1))
y1 = K1 * z1 / (1 + Ψ * (K1 - 1)): y2 = K2 * z2 / (1 + Ψ * (K2 - 1))
Ψ1 = Ψ - f / g
Loop While Abs(Ψ1 - Ψ) >= 0.0001
Text5.Text = Format$(Ψ1, "0.000000"): Text1.Text = Format$(X1, "0.000000"):
Text2.Text = Format$(X2, "0.000000"):
Text3.Text = Format$(y1, "0.000000"): Text4.Text = Format$(y2, "0.000000"):
End If
End Sub
Private Sub Command2_Click()
End
End Sub
End Sub
运行界面：
第二题
解：
编码：Private Sub Command1_Click()
Dim d1 As Single, d2 As Single, d3 As Single, d4 As Single, d31 As Single
Dim w1 As Single, w2 As Single, w3 As Single, w4 As Single
Dim α1 As Single, α2 As Integer, α3 As Single, α4 As Single
Dim D As Single, W As Single, G As Single, N As Single
α1 = 2.25: α2 = 1: α3 = 0.33: α4 = 0.21
w1 = 0: d1 = 20: w3 = 9: d3 = 1
Do
d3 = d31
D = 50 + d3 - 0.003 * (100 - D) + 0.0015 * D
W = 100 - D
d2 = 30 - 0.003 * W
d4 = 0.0015 * D
w2 = 0.003 * W
w4 = 40 - 0.0015 * D
G = (d2 / w2) / (d4 / w4)
N = Log(G) / Log(α2 / α4)
d31 = 10 * (d4 / w4) * (α3 / α4) ^ N / (1 + (d4 / w4) * (α3 / α4) ^ N)
w3 = 10 - d31
Loop While Abs(d31 - d3) >= 0.0001
Text1.Text = Format$(N, "0.00"): Text2.Text = Format$(d1, "0.0000"): Text3.Text = Format$(d2,
"0.0000")
Text4.Text = Format$(d31, "0.0000"): Text5.Text = Format$(d4, "0.0000"): Text6.Text =
Format$(w1, "0.0000")
Text7.Text = Format$(w2, "0.0000"): Text8.Text = Format$(w3, "0.0000"): Text9.Text =
Format$(w4, "0.0000")
End Sub
Private Sub Command2_Click()
End
End Sub
Private Sub Form_Load()
Print "以100摩尔进料为计算基准。根据题意定甲苯为轻关键组分，异丙苯为重关键组分."
Print
Print "从相对挥发度的大小可以看出，二甲苯为中间组分。在做物料衡算时，要根据它的"
Print
Print "相对挥发度与轻，重关键组分相对挥发度的比例初定在馏出液和釜液中的分配比。"
Print
Print "f1=20,f2=30,f3=10,f4=40,F=100,"
End Sub
运行界面如下图：
