r-k calculation.txt

来自「基于VB语言编写的热力学性质计算的应用软件（数据库操作）」· 文本 代码 · 共 18 行

    
    R=8.314
    Tc1=33.2
    Pc1=1.297
    a = 0.42748 * R ^ 2 * Tc1 ^ 2.5 / (Pc1 * 1000000)
    b = 0.08664 * R * Tc1 / (Pc1 * 1000000)
       p=
       V0=0.00020402
       T0 = P * V0/ R
       Do While 1   
         Cv=26.88+4.347*T0/1000-0.3265*T0^2/1000000
         V=(2687.402-Cv*(T0-298.15))/P
         T1 = P * (V - b) / R + a * (V - b) / (Sqr(T0) * V * (V + b) * R) '迭代方程式
         If Abs(T0 - T1) < 0.00001 Then Exit Do
         T0 = T1
       Loop
       t = T1