24小时热门版块排行榜

>论坛更新日志 (2106)
>虫友互识 (69)
>文献求助 (30)
>考博 (25)
>导师招生 (24)
>基金申请 (16)
>休闲灌水 (16)
>海外博后 (14)
>硕博家园 (13)
>找工作 (12)
>论文投稿 (10)
>人文社科 (5)
>博后之家 (5)
>有机交流 (3)
>生物科学 (3)
>考研 (3)

返回列表

yuyi2009

铜虫 (小有名气)

应助: 2 (幼儿园)
金币: 60.6
散金: 29
红花: 4
帖子: 289
在线: 174.3小时
虫号: 1103053
注册: 2010-09-19
性别: GG
专业: 计算机网络

[求助] VB入门解惑已有1人参与

刚开始学习VB，还不太熟悉，帮我看个题，谢啦~就是在主函数里面调用UNIFAC模型，不知道在哪里调试！
Option Explicit
Private Sub Command1_Click()
Dim NG As Integer, NC As Integer, T As Single, XC() As Variant, RK() As Variant, QK() As Variant, A() As Single, V() As Single
NG = 4
NC = 2
T = 298
XC = Array(0.1, 0.9)
RK = Array(0.9011, 0.6744, 0.1, 2.8266)
QK = Array(0.848, 0.54, 1.2, 2.472)
A(1, 1) = 0
A(1, 2) = 0
A(1, 3) = 986.5
A(1, 4) = 526.5
A(2, 1) = 0
A(2, 2) = 0
A(2, 3) = 986.5
A(2, 4) = 526.5
A(3, 1) = 156.4
A(3, 2) = 156.4
A(3, 3) = 0
A(3, 4) = -202.1
A(4, 1) = 50.49
A(4, 2) = 50.49
A(4, 3) = -25.87
A(4, 4) = 0
V(1, 1) = 1
V(1, 2) = 1
V(1, 3) = 1
V(1, 4) = 0
V(2, 1) = 0
V(2, 2) = 0
V(2, 3) = 0
V(2, 4) = 1
Call UNIFAC
End Sub
Sub UNIFAC(NG As Integer, NC As Integer, T As Single, XC() As Single, RK() As Single, QK() As Single, A() As Single, V() As Single, GAMMA() As Single)
'This program calculates liquid phase activity coefficients using UNIFAC method
'NC--number of components in the mixture
'NG-number of functional group in the mixture
'QK(k)-the area parameter of functional group
'RK(k)-the volume parameter of functional group
'A(k,m)-the interactional parameter of functional group
'V(i,k)-the number of "'k" functional group in the "i" components
'T-system temperature,K
'XC(i)-mol fraction of "i" components
'GAMMA(i)--activity coefficients of "i" components

Dim i As Integer, K As Integer, M As Integer, n As Integer
Dim SUM1 As Single, SUM2 As Single, SUM3 As Single, SUM4 As Single
Dim Q(10) As Single, R(10) As Single, THETA1(10) As Single, THETA2(10, 10) As Single, PHI(10, 10) As Single

SUM1 = 0
SUM2 = 0
SUM3 = 0
For i = 1 To NC
Q(i) = 0
R(i) = 0
For K = 1 To NG
   Q(i) = Q(i) + V(i, K) * QK(K)
   R(i) = R(i) + V(i, K) * RK(K)
Next K
SUM1 = SUM1 + XC(i) * (5 * (R(i) - Q(i)) - R(i) + 1) '计算SUM(x(i)*l(i))
SUM2 = SUM2 + XC(i) * Q(i)
SUM3 = SUM3 + XC(i) * R(i)
Next i

For i = 1 To NC
GAMMA(i) = Log(R(i) / SUM3) + 5 * Q(i) * Log(Q(i) * SUM3 / R(i) / SUM2)
GAMMA(i) = Exp(GAMMA(i) + 5 * (R(i) - Q(i)) - R(i) + 1 - R(i) * SUM1 / SUM3)
Next i

SUM1 = 0
For i = 1 To NC
SUM2 = 0
SUM3 = 0
For K = 1 To NG
      SUM2 = SUM2 + V(i, K) * XC(i)
      SUM3 = SUM3 + V(i, K)
Next K
SUM1 = SUM1 + SUM2
SUM4 = 0
For K = 1 To NG
      THETA2(i, K) = V(i, K) / SUM3
      SUM4 = SUM4 + QK(K) * THETA2(i, K)
Next K
For K = 1 To NG
      THETA2(i, K) = QK(K) * THETA2(i, K) / SUM4
Next K
Next i

SUM2 = 0
For K = 1 To NG
SUM3 = 0
For i = 1 To NC
   SUM3 = SUM3 + V(i, K) * XC(i)
Next i
THETA1(K) = SUM3 / SUM1
SUM2 = SUM2 + QK(K) * THETA1(K)
Next K

For K = 1 To NG
THETA1(K) = QK(K) * THETA1(K) / SUM2
For M = 1 To NG
      PHI(K, M) = Exp(-A(K, M) / T)
Next M
Next K

For K = 1 To NG
SUM1 = 0
SUM2 = 0
For M = 1 To NG
      SUM1 = SUM1 + THETA1(M) * PHI(M, K)
      SUM3 = 0
      For n = 1 To NG
         SUM3 = SUM3 + THETA1(n) * PHI(n, M)
      Next n
      SUM2 = SUM2 + THETA1(M) * PHI(K, M) / SUM3
Next M
Q(K) = QK(K) * (1 - Log(SUM1) - SUM2)
Next K

For i = 1 To NC
SUM4 = 0
For K = 1 To NG
      SUM1 = 0
      SUM2 = 0
      For M = 1 To NG
         SUM1 = SUM1 + THETA2(i, M) * PHI(M, K)
         SUM3 = 0
         For n = 1 To NG
            SUM3 = SUM3 + THETA2(i, n) * PHI(n, M)
         Next n
         SUM2 = SUM2 + THETA2(i, M) * PHI(K, M) / SUM3
      Next M
      SUM4 = SUM4 + V(i, K) * (Q(K) - QK(K) * (1 - Log(SUM1) - SUM2))
Next K
GAMMA(i) = GAMMA(i) * Exp(SUM4)

Next i

End Sub

End Sub

回复此楼