24小时热门版块排行榜

>论坛更新日志 (3447)
>虫友互识 (396)
>休闲灌水 (238)
>文献求助 (230)
>导师招生 (174)
>考博 (87)
>考研 (75)
>硕博家园 (54)
>博后之家 (53)
>论文投稿 (46)
>基金申请 (39)
>公派出国 (35)
>论文道贺祈福 (32)
>教师之家 (32)
>招聘信息布告栏 (28)
>找工作 (27)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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

回复此楼