最全空气物性计算-源代码-VBA-加载模块.docx
91 浏览量
2024-03-10
10:56:40
上传
评论
收藏 19KB DOCX 举报
'空气
'常数
Const R = 287.11 'J/kg.K，气体常数
Const M = 0.02896 'kg/mol，空气的分子量
'Const Pcr = 5 'MPa，临界压力
'Const roucr = 536 'kg/m3，临界密度
'Const tb = -185.86 '℃，沸点
'Const tm = -189.37 '℃，熔点
'Const cp0 = 5198 'J/kg.K
'Const cv0 = 3121 'J/kg.K
Const h0 = -5059 'J/kg 温度 315.56℃，压力 0.1013MPa 时空气的焓
Const s0 = 2607.17 'J/kg.K 温度 315.56℃，压力 0.1013MPa 时空气的熵
Const T0 = 273.15 'K
Const p0 = 100000 'Pa
'Const rou0 = 0.1762 'kg/m3
'空气状态方程
Function f_vpT_Air(V, P, T) 'p-MPa,T-K,v-m3/kg
Dim A0, a, B0, b, c
A0 = 157.204
a = 0.000666782
B0 = 0.0015922
b = -0.00038018
c = 1498.62
'p = p * 1000000
Dim f
f = ((R * T * (1 - c / (V * T ^ 3))) / (V ^ 2)) * (V + B0 * (1 - b / V))
f = f - (1 - a / V) * A0 / V ^ 2
f_vpT_Air = f - P
'f = f
On Error Resume Next
End Function
'A.1. 比容，m3/kg
Function V_PT_Air(P, T)
Dim vi(10000)
vi(1) = 0.01
vi(2) = 2
WuCha = 0.00001
If f_vpT_Air(vi(1), P, T) <> 0 And f_vpT_Air(vi(2), P, T) Then
剩余12页未读，继续阅读
资源评论
