泡点压力及气相组成.txt

Private Sub Form_Load()
Rem calculating the bubble points for ideal systems           '程序说明
Dim N As Integer, i As Integer, T As Single, x(10) As Single
Dim A(10) As Single, B(10) As Single, C(10) As Single, ps(10) As Double, pc, yy As Double, y(10) As Double, Tl As String
Const R# = 8.314                                          '定义R为常数8.314
Open "D:\VB\3in.dat" For Input As #6          '打开带路径的输入数据文件
Input #6, Tl                                              '读入标题
Input #6, N, T                                     '读入组分数N，温度T/K，压力p/kPa，阻尼系数
Open "D:\VB\OutTb.dat" For Output As #8        '打开带路径的输出数据文件
Print #8, Tab(6); Tl                                      '在输出文件中打印标题
Print #8, Tab(12); "输入数据"
Print #8, Tab(2); "N="; Format(N, "#0"); "  T/K="; Format(T, "####0.00")  '在输出文件中打印N,T
For i = 1 To N
Input #6, x(i), A(i), B(i), C(i)                          '利用循环语句读入xi，安托尼参数Ai,Bi,Ci
Print #8, Tab(2); "x"; Format(i, "#0"); "="; Format(x(i), "#0.00"); "  A"; Format(i, "#0"); "="; Format(A(i), "#00.0000"); _
          "  B"; Format(i, "#0"); "="; Format(B(i), "####0.00"); "  C"; Format(i, "#0"); "="; Format(C(i), "###0.00")
Next i
Close #6                                                  '关闭输入文件
Print #8, Tab(12); "计算结果"
pc = 0#                                                   'pc为计算的压力
For i = 1 To N
ps(i) = Exp(A(i) - B(i) / (T + C(i)))                     '计算饱和蒸汽压
pc = pc + ps(i) * x(i)                                    '计算压力
Next i
  Print #8, Tab(10); "  p/Pa="; Format(pc, "####0.000")            '输出迭代次数和相应的泡点温度
yy = 0#
For i = 1 To N                                            '对气相组成归一化处理
y(i) = ps(i) * x(i) / pc
yy = yy + y(i)
Next i