KEYENCE KVスクリプトで、外気冷房判断
・外気冷房判断のために、温度 - 相対湿度、露点温度からエンタルピを求めます。
※絶対湿度 [kg/m3] 換算の間違いを修正 （2015/07/11)
※絶対湿度 [kg/m3] の計算式を変更 （2015/07/12)
※露点温度 [℃DP] から相対湿度 [%RH] の計算を追加 （2015/07/12)
※定数の間違い等あれば、指摘して下さい。
'' ***********************************************
'' 温度 - 相対湿度から、露点温度、エンタルピを求める
'' 温度 - 露点温度から相対湿度を求める
'' KEYENCE KVスクリプト
'' ***********************************************
'' 温度
dm0.df = 25.0
'' 絶対温度
dm4.df = dm0.df + 273.15
'' 水の飽和水蒸気圧 es (Pa)
'' [Sonntagの式]
'' 未使用
dm10.df = exp(-6096.9385 / dm4.df + 21.2409642 - 0.02711193 * dm4.df + _
dm4.df * dm4.df * 0.00001673952 + _
2.433502 * log(dm4.df))
'' [Tetensの近似式]
dm10.df = 6.11 * 10^(7.5 * dm0.df / (dm0.df + 237.3)) * 100.0
'' メモ：実数対応べき乗計算の例： 4^5 = exp(5 * log(4))
'dm10.df = 6.11 * exp(7.5 * dm0.df / (dm0.df + 237.3) * log(10)) * 100.0
'' [Goff-Gratch + Okada]
'' こちらを使用(-30℃～35℃では、精度がある。JISの飽和水蒸気圧表に近い）
dm14.df = exp(1.809378 + 7.266115E-2 * dm0.df - 3.003879E-4 * dm0.df * dm0.df + _
1.181765E-6 * dm0.df * dm0.df * dm0.df + _
- 3.863083E-9 * dm0.df * dm0.df * dm0.df * dm0.df) * 100.0
'' [Wexler-Hylandの式]
'' 未使用
'dm14.df = exp(-5800.2206 / dm4.df + 1.3914993 - 4.8640239E-2 * dm4.df _
' + 4.1764768E-5 * dm4.df * dm4.df _
' - 1.4452093E-8 * dm4.df * dm4.df * dm4.df _
' + 6.5459673 * log(dm4.df))
'' 相対湿度(%RH)
dm20.df = 50.0
dm24.df = dm20.df / 100.0
'' 水蒸気圧 e (Pa)
dm28.df = dm14.df * dm24.df
'' 露点温度(℃DP)
'' 未使用
'' 計算用テンポラリ変数
'dm32.df = log(dm28.df / 611.213)
'if dm32.df >= 0 then
' dm40.df = 13.715 * dm32.df + dm32.df * dm32.df * 0.84262 + _
' dm32.df * dm32.df * dm32.df * 0.019048 + _
' 0.0078158 * dm32.df * dm32.df * dm32.df * dm32.df
'else
' dm40.df = 13.7204 *dm32.df + dm32.df * dm32.df * 0.736631 + _
' dm32.df * dm32.df * dm32.df * 0.0332136 + _
' dm32.df * dm32.df * dm32.df * dm32.df * 0.000778591
'end if
'' 露点温度(℃DP)
'' こちらを使用
dm60.df = 237.3 / (7.5 / log10((dm28.df / 100) / 6.1078) -1)
'' 上の式から、水蒸気圧を逆算
dm64.df = 10^(7.5 * dm60.df / (237.3 + dm60.df)) * 6.1078 * 100
'' 逆算した相対湿度(%RH)
dm68.df = dm64.df / dm14.df * 100
'' 水の分子量 = 18.015 [g/mol]
'' 空気の平均分子量 = 28.967 [g/mol]
'' 空気の全圧 = 海抜0mの大気圧（101.325[kPa]）
'' 絶対湿度 [kg/kgDA]
dm44.df = 18.015 * dm28.df / (28.967 * (101325 - dm28.df))
'' 乾き空気の比重量は 1.293 [kg/m3] 湿り空気は一般的に 1. 2 [kg/m3]
'' 絶対湿度 [kg/m3]
'dm48.df = dm44.df * 1.2
'' 絶対湿度 [kg/m3]
'' こちらを使用
dm48.df = 2.16679 * dm28.df / dm4.df / 1000
'' 空気の比熱（定圧比熱容量）= 1.006 [kJ/kg･Ｋ]
'' 水の蒸発熱（0℃基準）= 2501 [kJ/kg]
'' 蒸気の比熱 = 1.805 [kJ/kg･Ｋ]
'' 比エンタルピ [kJ/kgD.A.]
dm52.df = 1.006 * dm0.df + (2501.0 + 1.805 * dm0.df) * dm44.df
'' [kcal/kgDA]は、4.18605を乗じて[kJ/kgDA]となる
'' 比エンタルピ [kcal/kgD.A.]
dm56.df = dm52.df / 4.18605