冷凍サイクル計算プログラムソフト（自然冷媒）
操作説明書
日本冷凍空調学会

冷媒技術分科会

目次
１．はじめにP2

２．プログラムの実行

P2

３．各シートの使い方

P33.1）「飽和表」シート

P33.2）「基本サイクル」シート

P43.3）「EPR」シート

P63.4）「二段」シート

P83.5）「ｐ−ｈ線図」シート

P103.6）「Ｔ−ｓ線図」シート

P113.7）「検算飽和」シート

P123.8）「検算過熱」シート

P13

４．計算上の注意

P14

５．折込関数（計算精度ほか）

P14

６．使用上の注意

P15

原著文献

P15

付録

折込関数の応用ほか

A1．圧力基準飽和蒸気表

P16A2．二段圧縮一段膨張冷凍サイクル

P18

A3．折込関数（Function）一覧

P27

− 1 −

１．はじめに
このプログラムソフトは、一段圧縮、二段圧縮、蒸発圧力調整弁の冷凍サイクル計算ができます（ただし二酸化炭素は一段圧縮のみ掲載）。冷媒の状態式は Excel の関数として組み込ま
れています。各関数を使えば、ユーザが独自の冷凍サイクル･シミュレーションなどを Excel
上に作ることができます。付録は自作ソフトを作る際の参考用です。
本プログラムソフト扱う冷媒は表１に示す自然冷媒です。状態式は本学会が推奨するものに基づきます。HC 冷媒は Miyamoto and Watanabe の状態式 １ , ２ , ３ によります。この式はASHRAE Handbook６にも採用され、ｐ−ｈ線図、飽和表などが掲載されています。アンモニ
アは本学会冷媒蒸気表委員会,SI 化小委員会７が採択した Baehr and Tillner-Roth の状態式４によります。二酸化炭素は Span and Wanger の状態式 ５ によります。この式は ASHRAE
Handbook８にも採用されています。
本ソフトは Microsoft 社の Excel 上で作動します。従って、計算環境として Excel がインストールされたパソコンが必要です。Excel のバージョンは Excel97 以上を推奨します。
表１
No
１
２
３
４
５

冷凍サイクル計算プログラムソフトの冷媒（自然冷媒）

ファイル名
R 290(プロパン).xls
R 600(ブタン).xls
R 600a(イソブタン).xls
R 717(アンモニア).xls
R 744(二酸化炭素).xls

種類
HC
HC
HC
無機物
無機物

分子式
CH3CH2CH3
CH3CH2CH2CH3(CH3)2CHCH3
NH3
CO2

使用例ほか
冷凍冷蔵
冷凍冷蔵
冷凍冷蔵
冷凍冷蔵

原著
文献１
文献２
文献３
文献４
文献５

２．プログラムの実行
計算したい冷媒のファイルをダブルクリックすると Excel が起動し、下記ダイアログが表示されます。ここでは「マクロを有効にする」をクリックします。配布されている CD 上のファ
イルはウィルスに汚染されていないことを確認済です。

図１．マクロを有効にするダイアログ「マクロを無効にする」をクリックすると、以前の計算結果が表示されるのみで、新しい計
算はできません。
なお、上記のダイアログが出ず直ちに Excel が立上る場合があります。これは、使っている
Excel のセキュリティレベルが高(H)に設定されている場合です。本ソフトはレベル高(H)に対応していないので、その時は [ツール]→[マクロ]→[セキュリティ]の順にクリックして レベル
中(M)を選択してください。レベル低(L)は勧められません。

− 2 −３．各シートの使い方
Excel が立上ると、表２に示すシートが表示されます。
表２

プログラムに含まれるシート一覧

シート名
シートの機能「飽和表」
温度基準・飽和蒸気表「基本サイクル」 一段圧縮冷凍サイクルの計算 （含：ポリトロープ圧縮）「EPR」
蒸発圧力調整弁付きの理論冷凍サイクル計算「二段」
二段圧縮一段膨張の理論冷凍サイクル計算「p-h 線図」「二段」「EPR」「基本サイクル」のデータを書込みできます「T-ｓ線図」「二段」「EPR」「基本サイクル」のデータを書込みできます「検算飽和」
飽和蒸気域の計算誤差を示すシート「検算過熱」
過熱蒸気域の計算誤差を示すシート

詳細説明
3.1)
3.2)
3.3)
3.4)
3.5)
3.6)
3.7)
3.8)

注）計算では各部圧力損失は無視します。なお、理論冷凍サイクルは等エントロピー圧縮です。
3.1) 「飽和表」シート

図２「飽和表」シート（R 290 の例）

この表は温度基準で、１℃きざみで計算したものです。各セルには計算に使用した折込関数がそのまま残してあり、適宜応用することができます。例えばもっと細かい温度きざみや、圧力基準の飽和表などに応用可能です。ただし、変更内容を間違えると元の設定に戻らなくなるので、注意して行ってください。なお、元の設定に戻らなくなった場合には、配布されているCD から再度インストールし直せば復帰可能です。
折込関数の行にも関数名を記入してあります。関数の詳細は５章を参照してください。
− 3 −3.2) 「基本サイクル」シート（一段圧縮冷凍サイクル）

図３ 「基本サイクル」シート（R 290 の例）

このシートを使って、一段圧縮の冷凍サイクル

圧力 （MPa）を計算できます。右図に各点の番号を示します。

２

３

１は圧縮機吸込み、２は圧縮機吐出し、３は膨張
弁入口、４は蒸発器入口です。

ｖ１

シートの黄色ハッチング部分に数値を入力し１

４てください。列ごとに、圧縮機吐出圧力〜成績係

h４

数までの項目を計算します。

h１h２

比エンタルピ （kJ/kg）

E 列データは「ｐ−ｈ線図」、「T−ｓ線図」の図４ 一段圧縮 冷凍サイクル

シートに赤い線で書込まれます。二酸化炭素については F 列のデータも「線図」と連動します。それ以外の列は「線図」と連動しません。
２３log(p)

３

２

２

T

２

３
３

５
４

１

1

５

1

４

図５ p-h 線図h

図６ T-s 線図

− 4 −

１’

s

「基本サイクル」シート８行目の圧縮機吐出温度の入力では、実際の圧縮機と理論冷凍サ
イクルとを区別できます。実際の圧縮機で計算したい場合は、実機の圧縮機吐出し温度を［E8］セルに入力してください。すると、ポリトロープ圧縮（ pv n = 一定 ）による計算を行います。理論冷凍サイクルで計算する場合は、［E8］セルをブランクにしてください。こちらは、等エントロピー圧縮の計算を行います。
なお、理論圧縮・実際圧縮共に、サイクル各部の圧力損失は無いものとしています。
冷凍サイクルを描いていないｐ−ｈ線図が必要な場合は、E 列の凝縮温度、蒸発温度（E4、E5 セル）を消去します。ただし、「基本サイクル」シートの他に、「EPR」「二段」シートでも線図に書込んでいるので、そちらのデータも同様に消去する必要があります。
データを消去したとき、いくつかのセルに＃VALUE が表示される場合がありますが、プ
ログラムの動作に問題はありません。再びデータをインプットすれば正常に計算します。
「基本サイクル」シートの形は冷媒によって変えてあります。二酸化炭素については圧力を入力します。それ以外の冷媒では温度基準で入力します。
「基本サイクル」シートは計算式が埋めこんであり、誤操作によって式を壊してしまわな
いように保護をかけてあります。従って、黄色ハッチング部分以外には書込みができません。
Excel 上級者ならば、保護を外し独自の冷凍サイクルなどに応用展開が可能です。ただし、元に戻らなくなる恐れがあるので、初級者は保護を外さないことをお薦めします。
補足）シート保護の外し方
追記や修正などをする場合は、下記方法でシート保護を外してから行ってください。
１）「ツール」→「保護」→「シート保護の解除」の順にクリックします。すると、シート保護が解除されます。解除パスワードは省略しており不要です。
２）修正が完了したら、再度シート保護をかけておく事を推奨します。
「ツール」→「保護」→「シートの保護」の順にクリックすると保護がかかります。
そのとき図７のダイアログが出ます。必要ならばパスワードを設定してください。設定
したパスワードは保護の解除時に必要になるので、憶えておく必要があります。
ダイアログの保護対象は全てチェックを付けたままを推奨します。
なお、同様の方法でブックの保護もかけられますが、動作が不安定になる場合がありますのでお勧めしません。
（ブックとは各シートを含む Excel ファイ
ル全体を意味します）

図７

− 5 −

シート保護のダイアログ

3.3) 「EPR」シート（ 蒸発圧力調整弁（EPR）を用いたサイクル）
注）二酸化炭素にはこのシートはありません。

図 8 「EPR」シート（R 290 の例）

図９に示す２つの蒸発器を持ち、そのどちらか又

2は両方に蒸発圧力調整弁が設置されている冷凍サイ

COMP

クルを計算できます。計算できるのは理論冷凍サイクルです。（等エントロピー圧縮で圧力損失は無視）
黄色ハッチング部にデータを入力してください。
蒸発温度や蒸発器出口の過熱度は蒸発器ごとに入力

COND

6

3

EPR

5 EVAPⅠ 7
3

1
8

EVAPⅡ

4

図 9 EPR 冷凍サイクルが必要です。EPR の設定圧力に相当する温度を蒸発
温度のセルに入力してください。二酸化炭素では、蒸発温度のかわりに蒸発圧力で入力します。
なお、EVAPⅠを EPR 側とすることを推奨します。
両方の蒸発器共に EPR 付とする場合は、同じ様に設定する蒸発温度を入力してください。
次に、２つの蒸発器に流れる冷媒流量比を入力します。EVAPⅠ側流量／全体流量を流量比としており、流量比＝１とすると EVAPⅡの無い蒸発器Ⅰだけの EPR サイクルになります。
計算した EPR サイクルはｐ−ｈ線図、Ｔ−ｓ線図に緑色の線で書込まれます。書込まれる− 6 −のはＥ列データだけです。書込んだサイクルは下図の形をしています。図９と図 10〜11 は記号を合わせてあるので各点の対応状況を確認してください。

log(p)

２

３

T

2

3

５
４

4

７

7

８１６

図 10 p-h 線図

h

4

5

5

１

7

図 11 T-s 線図s各状態点を表３に示します。ダッシュ付きの番号は直近の飽和線上の点を意味します。
表３
点
１
２
３
４

EPR（蒸発圧力調整弁）冷凍サイクルの状態点

状態点
圧縮機吸込み（点８､点６の混合）圧縮機吐出し膨張弁入口
EVAPⅠ入口

点
５
６
７
８

状態点
EVAPⅠ出口、EPR 入口
EPR 出口EVAPⅡ入口EVAPⅡ出口

冷凍サイクルを描いていない線図が必要な場合は、E 列の凝縮温度、蒸発温度（E4、E5セル）を消去してください。二酸化炭素については F 列も消去してください。同様に、「基本サイクル」「二段」シートのデータも消去する必要があります。その際に、いくつかのセルに＃VALUE が表示される場合がありますが、プログラムの動作に問題はありません。再び
データを書込めば正常に計算します。
「EPR」シートにも計算式が埋めこんであり、誤操作によって式を壊してしまわないように保護をかけてあります。従って、黄色ハッチング部分以外に書込みはできません。Excel 上
級者ならば、保護を外し独自の冷凍サイクルなどに応用展開が可能です。ただし、元に戻らなくなる恐れがあるので、初級者は保護を外さないことを推奨します。
入力する際の注意事項は 3.2)項と同じですので、そちらも参照してください。

− 7 −3.4）「二段」シート（二段圧縮一段膨張 理論冷凍サイクル）
注）二酸化炭素にはこのシートはありません。

図 12 「二段」シート（R 290 の例）

図５に示す、二段圧縮一段膨張の冷凍サイクルを計算できます。中間冷却器および膨張弁を２つ

４

COND

５

です。なお計算できるのは理論冷凍サイクルです３６（等エントロピー圧縮で圧力損失は無視）。
他の入力シートと同様に、黄色ハッチング部に

２

中間冷却器

持ち、同一冷媒がサイクル内を循環するサイクル

入力します。
中間冷却器 CF`（コンタクトファクター）等については、付録Ａ２章を参照してください。計算

１

EVAP

８

７

図 13 二段圧縮一段膨張冷凍サイクル

方法を説明してあります。
計算した二段圧縮冷凍サイクルはｐ−ｈ線図、Ｔ−ｓ線図に青い線で書込まれます。ただし、書き込まれるのはＥ列データだけです。書込んだサイクルは図 14、図 15 の形をしています。
両図は、図 13 と記号を合わせてあるので、各点の対応を確認してください。

− 8 −log(p)

７ ５
5
6

4

T

４

4

5

６

3

３

7
6

２

4

5

2
3

6

3

1

8
８

1

1

8

8

1

s

h
図 14

p-h 線図

図 15

T-s 線図各状態点を表４に示します。ダッシュ付きの番号は直近の飽和線上の点を意味します。
表４
点
１
２
３
４

二段圧縮一段膨張冷凍サイクルの状態点

状 態 点
低段圧縮機吸込み
低段圧縮機吐出し
高段圧縮機吸込み（点２､点６の混合）高段圧縮機吐出し

点
５
６
７
８

状 態 点
高段用膨張弁入口
中間冷却器液側入口
低段用膨張弁入口
蒸発器入口

冷凍サイクルを描いていない線図が必要な場合は、［C5］、［C6］、［C7］セルのデータを消去してください。同様に、「基本サイクル」「EPR」シートのデータも消去する必要があります。その際に、いくつかのセルに＃VALUE が表示される場合がありますが、プログラムの動作に問題はありません。再びデータを書込めば正常に計算します。
「二段」シートにも計算式が埋めこんであり、誤操作によって式を壊してしまわないように保護をかけてあります。従って、黄色ハッチング部分以外に書込みはできません。Excel 上級者ならば、保護を外し独自の冷凍サイクルなどに応用展開が可能です。ただし、元に戻らなくなる恐れがあるので、初級者は保護を外さないことを推奨します。
入力する際の注意事項は 3.2)項 3.3)項と同じですので、そちらも参照してください。

− 9 −3.5）「ｐ−ｈ線図」シート
温度目盛

℃

臨界点

液相線

気相線
図 16 「ｐ−ｈ線図」シート

等エントロピー線

等温線

等比体積線

（R 290 の例）

ｐ−ｈ線図は縦軸を圧力（対数目盛）、横軸を比エンタルピーとする線図です。気相線、液相線上に温度を２℃刻みで目盛ってあります。なお、本プログラムソフトでは圧縮液域の等温度線や二相域の等乾き度線は省略しています。
上記ｐ−ｈ線図には、各シートの入力データに基づく冷凍サイクルが表示されています。
赤い線は「基本サイクル」シートデータによる一段圧縮冷凍サイクル、緑の線は「EPR」シ
ートデータによる蒸発圧力調整弁の理論冷凍サイクル、青い線は「二段」シートデータによる二段圧縮一段膨張の理論冷凍サイクルを示します。
冷凍サイクルの線の消し方は、3.2)〜3.4)の説明を参照してください。
Ctrl

と

g

のボタンを同時に押すと、保存・出力用の線図をイメージデータとして作ることができます。イメージデータにはプログラムが付きませんので、ファイルを小さくできます。また、詳細情報は添付されませんので、情報を保護することができます。
なお、「ｐ−ｈ線図」シートも誤操作予防のために保護をかけてあります。Excel 上級者ならば、保護を外して独自のｐ−ｈ線図を作ることも可能です。

− 10 −

3.6）「Ｔ−ｓ線図」シート
例として、R290 のＴ−ｓ線図を示します。

臨界点

等比体積線

液相線

等乾き度線

気相線

図 17 「T−s 線図」シート

等圧線
等エンタルピー線（R 290 の例）

T−ｓ線図は縦軸を温度（絶対温度）、横軸を比エントロピーとする線図です。等圧線、等エンタルピー線、等比体積線などが書込まれています。この線図ではカルノーサイクルが長方形になります。種々冷凍サイクルを書込んだときに、長方形に近いものほど高効率になります。
上記Ｔ−ｓ線図には、各シートの入力データに基づく冷凍サイクルが示されています。
赤い線は「基本サイクル」シートデータによる一段圧縮冷凍サイクル、緑の線は「EPR」シ
ートデータによる蒸発圧力調整弁の理論冷凍サイクル、青い線は「二段」シートデータによる二段圧縮一段膨張の理論冷凍サイクルを示します。
冷凍サイクルの線の消し方は、3.2)〜3.4)の説明を参照してください。
Ctrl

と

g

のボタンを同時に押すと、保存・出力用の線図をイメージデータとして作ることができます。イメージデータにはプログラムが付きませんので、ファイルを小さくできます。また、詳細情報は添付されませんので、情報を保護することができます。
なお、「Ｔ−ｓ線図」シートも誤操作予防のために保護をかけてあります。Excel 上級者ならば、保護を外して独自のＴ−ｓ線図を作ることも可能です。

− 11 −

3.7）「検算飽和」シート（飽和蒸気域の計算誤差を示すシート）

例として、R 290 の「検算飽和」シートを示します。
このシートは飽和蒸気域での本ソフトの計算誤差を示しており、原著または原データに記載された数値と本プログラムソフトの計算値とを対比して誤差を計算しています。表中、原デー
タの列は原著または原データに記載された数値で、本ソフトの列は本プログラムソフトの計算値です。差の列は誤差を示しており、［本ソフト計算値］／［原著または原データ］−１

を％で表示しています。
本シートに示された範囲で、誤差はほとんど 0.00％以下と小さいことが確認できています。
ただし、原データの有効桁数が小さい場合などでは、四捨五入のために誤差の％値が大きくなっていますので、その領域のデータは丸め誤差があることに注意してください。
また下表にて、飽和蒸気の比体積 v" にて−150℃の場合などでは誤差が少し大きくなっていますが、その理由は分かっていません。このような場合には、本ソフトの計算値はシートに示された誤差を含むとして取扱ってください。

図 18 「検算飽和」シート（R 290 の例）− 12 −

3.8）「検算過熱」シート（過熱蒸気域の計算誤差を示すシート）

注）HC 冷媒（R 290、R 600、R 600a）にはこのシートはありません。
例として、R 717 の「検算過熱」シートを示します。
このシートは加熱蒸気域での本ソフトの計算誤差を示しており、原著または原データに記載された数値と本プログラムソフトの計算値とを対比して誤差を計算しています。表中、原デー
タの列は原著または原データに記載された数値で、本ソフトの列は本プログラムソフトの計算値です。差の列は誤差を示しており、［本ソフト計算値］／［原著または原データ］−１

を％で表示しています。
本シートに示された範囲で、誤差はほとんど 0.00％程度と小さいことが確認できています。
ただし、原データの有効桁数が小さい場合などでは、四捨五入のために誤差の％値が大きくなっていますので、その領域のデータは丸め誤差があることに注意してください。
また下表のように、誤差が±0.03％程度になっている場合がありますが、その理由は分かっていません。ほとんどは丸め誤差であると思われますが、このような場合には本ソフトの計算
値はシートに示された誤差を含むとして取扱ってください。

図 19 「検算過熱」シート（R 717 の例）− 13 −

４．計算上の注意
このプログラムソフトでは、ニュートン法などの逐次近似法を多く用いています。
逐次近似法は求める解を挟むある範囲内で計算を繰返すので、場合によって式の適用範囲を超えて計算してしまう場合があります。従って、臨界点の近くや低温の領域などでは計算できない場合が発生します。線図に示された範囲や飽和蒸気表に示された範囲は、計算できた結果を表示したものです。なるべくこの範囲で計算してください。ただし、臨界点付近の領域では表示されている範囲内であっても計算できない等の例外があります。
臨界点のどの位近くまで計算できるかは冷媒によって異なります。これはプログラムの動作を軽くする目的で、使用しているアルゴリズム毎に計算できない範囲が異なるためです。計算できない場合は、プログラムは強制的に中断されます。したがって、臨界点近くの領域の熱力学性質を検討する場合は、原著文献を参照してください。
計算結果が無限大になってしまう場合は、Excel が計算を中断します。この場合、該当セルに #VALUE が表示されます（Excel での中断なので、図 20 の警告は出ません）。入力したデータに誤りがあったのが原因です。適切なデータを再入力すれば正常な状態に戻ります。
計算が収束しない場合は、Excel が延々と計算を続けるのを避けるために、プログラムが自分で計算を中断します。この場合は図 20 の警告表示が出ます。 ＯＫをクリックすると、以前の状態に戻ります。計算できない領域のデータを入力したことなどが主な原因です。適切なデータを再入力すれば正常な状態に戻ります。

図 20 警告ダイアログ

５．折込関数（計算精度ほか）
折込関数は、「飽和表」シートに組み込まれた関数以外にもありますので、以下に紹介します。これらは Excel の関数と同じように使えます。ただし、その冷媒のプログラム（Excel ではブックと言う）に各関数は付属しており、当該ブック内のシートで使えますが、他の Excel
ファイル（ブック）では動作しません。
冷媒が違えば同じ名前の関数でも、状態方程式も定数値も違うので、計算結果は同じにはなりません。従って、関数はブック内のみで使えるようになっています。
プログラムの計算結果と原著文献に掲載された値とを対比した表を、各冷媒の「検算飽和」、「検算過熱」シートに表示してあります。計算精度や計算可能範囲の確認に使ってください。
付録の折込関数（Function）一覧（P27）も参照ください。

− 14 −

表 5）

自然冷媒

（R 290、R 600、R 600a、R 717、R 744、）

求める状態量

関数名

（単位）

変 数

飽和圧力

PsT(T)（MPa）

T=温度（K）

飽和温度

TsP(P)（K）

P=圧力（MPa）

圧力

PVT(V,T)（MPa）

V=比体積（m3/kg）、T=温度（K）

比体積

VPT(P,T)（m3/kg）

P=圧力（MPa）、T=温度（K）

比ｴﾝﾀルﾋﾟｰ

HVT(V,T)（kJ/kg）

比ｴﾝﾄﾛﾋﾟｰ

SVT(V,T)（kJ/(kg･K))

V=比体積（m3/kg）、T=温度（K）

温度

TPS(P,S)

（K）

飽和
状態

蒸気
(気相)

P=圧力（MPa）、S=比ｴﾝﾄﾛﾋﾟｰ(ｋJ/(kg･K))

飽和
液体
(液相)

比ｴﾝﾀﾙﾋﾟｰ

HfT(T)（ｋJ/kg)比ｴﾝﾄﾛﾋﾟｰ

SfT(T)（kJ/(kg･K))

比体積

VfT(T)（m3/kg）

圧縮

比ｴﾝﾀﾙﾋﾟｰ

HfPT(P,T)（ｋJ/kg）

液体

比ｴﾝﾄﾛﾋﾟｰ

SfPT(P,T)（kJ/(kg･K))

(液相)

比体積

VfPT(P,T)（m3/kg）

T=飽和液の温度（K）

P=圧力（MPa）、T=温度（K）

６．使用上の注意
このプログラムソフトを使用することによって、ユーザが受ける可能性がある全ての損害等について、配布元、開発者は賠償請求の責を免れるものとします。
このプログラムは日本冷凍空調学会 冷媒技術分科会が、作成、監修、検査、承認したものです（開発担当 冷媒技術分科会 宇田川義紘）。
平成１７年４月日

配布元 社団法人日本冷凍空調学会

文献
１ H.Miyamoto, K.Watanabe: Thermodynamic Property Model or Fluid-Phase Propane,
International Journal of Thermophysics, Vol.21, No.5,pp.1045-1072, 2000
２ H.Miyamoto, K.Watanabe: Thermodynamic Property Model or Fluid-Phase n-Butane,
International Journal of Thermophysics, Vol.22, No.2, pp.459-475, 2001
３ H.Miyamoto, K.Watanabe: Thermodynamic Property Model or Fluid-Phase Isobutane,
International Journal of Thermophysics, Vol.23, No.2, pp.477-499, 2002
４ H.D.Baehr, R.Tillner-Roth: Thermodynamic Properties of Environmentally Acceptable
Refrigerants,ISBN 3-540-58693-8,Springer-Verg,1995
５ R.Span, W.Wanger: A New Equation of State for Carbon Dioxide, J.Phys.Che.Ref.Data,
1996
６ ASHRAE Handbook, Fundamentals, pp.20.44-20.49, 2005
７ 小口幸成,高石吉登,東之弘,矢田直之,並木孝行：冷凍,Vol.71,No.826, pp.826- ,1996-8
８ ASHRAE Handbook, Fundamentals, pp.20.38- , 2005

− 15 −

（付録）

折込関数の応用

本文５章で説明した折込関数を使って、専用のデータ処理などを行うことができます。その例として、圧力基準飽和表と二段圧縮理論冷凍サイクルのプログラム作成例を示します。二段
圧縮は本体に織り込み済みですが、ここで説明する方法で作ったものです。
Excel の操作方法についてはある程度知っているものとして、細かい説明は省略します。不明な点は Excel のヘルプを見てください。なお、データ処理を追加していくと、ファイルが重くなり計算時間が長くなります。パソコンの性能との兼ね合いでデータ処理する量を決めてください。
A１．圧力基準飽和表
新しいシートに圧力基準の飽和表を作る例を説明します。
A１.１

新しいシートを挿入

まず新しいシートを挿入してください。［挿入］→［ワークシート］の順にクリックすれば新しいシート Sheet1 ができます。このシートに表のタイトルや項目などを書き込んでください。シートのタブ名も必要により変更してください。
次に、A 列に計算したい圧力の値を適宜書き込んでください。サンプルでは 0.1 から 1.5MPa
まで 0.1MPa きざみで書き込みましたが、使い方に従って自由に決めてください。書込んだサ
ンプルを下図に示します。赤い字や線は説明用ですので書込み不要です。

この範囲に関数を書き込みます− 16 −

A１.２

各セルへ関数書き込み

まず B7 のセルに=TsP(A7)-273.15と書込み Enter してください。瞬時に計算して–42.4 と表示されます。以下同様に各セルにそれぞれの関数を書込みます。
計算精度は原著文献に従ってください。R 290 はおおむね有効数字５桁で表示されています
が、サンプルは４桁にしました。計算精度は原著文献に従うことを推奨しますが、文献が手元にない場合等には有効数字５桁以下にしてください。
どの関数も絶対温度で作ってあるので、273.15 のセルシウス度(℃)への温度換算は必ず行ってください。
B7
C7
D7
E7
F7
G7
H7

=TsP(A7)-273.15=VfT(B7+273.15)=VPT(A7,B7+273.15)=HfT(B7+273.15)=HVT(D7,B7+273.15)=SfT(B7+273.15)=SVT(D7,B7+273.15)

飽和温度
飽和液体の比体積
飽和蒸気の比体積
飽和液体の比エンタルピー
飽和蒸気の比エンタルピー
飽和液体の比エントロピー
飽和蒸気の比エントロピー

これで B7 から H7 まで一通り計算がおわりました。この B7〜H7 のセルを他のセルにコピーすれば飽和表が完成します。出来たサンプルを下図に示します。罫線など適宜追加して読みやすいようにしてください。

− 17 −

A２．二段圧縮の理論冷凍サイクル
本プログラムソフトで取扱うことができる二段圧縮の理論冷凍サイクルを図１、図２に示します。中間冷却器を用いる一段膨張サイクルで、各膨張弁が各圧縮機吸込状態での過熱度を制御しています。
計算シートの黄色ハッチングされたセルにデータを入力すれば、自動的にこの二段圧縮理論
冷凍サイクルを計算し、ｐ−ｈ線図上に冷凍サイクルを描きます。
冷凍サイクルとしての線図が不要な場合は、圧力（Pk、Pm、Po）の入力データを消去してください。その際、いくつかのセルに＃VALUE と表示されますが、動作上に問題はありません(３.２項 参照)。
圧力(MPa)

４

凝縮器

Pk
中間冷却器

３

５

６

２

１

蒸発器

８

７ ５

Pm

９

Po

８

４

ｖ3
６

３

ｖ１

１
比エンタルピー、h (kJ/kg)

CF’

７

２

(CF’=100%)

図１

冷凍サイクル図２

ｐ−ｈ線図

プログラムを使う場合、コンタクトファクター、CF の意味と使い方が必要ですので以下に説明します。

h 7 を求める式は、二段圧縮一段膨張に特有な以下の式です。
h 7 = h 2 − R ×( h 3 − h 6 )ここで、R は冷媒流量比（高段流量／低段流量）です。
Ref. 冷凍空調技術冷凍編、p.51、日本冷凍空調学会、東京、(H12-7-31)
二段圧縮の計算には、この式の他に中間冷却器の性能が必要です。ここでは中間冷却器の性能に次の CF を使っています。
CF ＝1−( t 7 − t m )／( t 5 − t m )t m は中間圧力 P m に対応する飽和温度です。

CF ＝0％は熱交換器無し、CF ＝熱交換が 100％(完全に)実現する場合になります。但し CF
は t m を用いて求めているので、厳密なコンタクトファクターではありません。厳密な CF 値との関係は実機で確認が必要です。
流量比 R には適用範囲があります。最大は R max ＝ ( h 2 − h 9 )／ ( h 3 − h 5 ) であり、t 7 ＝ t mの時、すなわち CF ＝100％で起ります。最小は R min ＝ ( h 2 − h 5 )／ ( h 3 − h 5 ) であり、高段の圧縮機吸込み状態での過熱度が液の直接冷却のみで制御される時、すなわち t 7 ＝ t 5 の時、
CF ＝0％で起ります。CF はこの R max と R min との間でどの辺りかを示すものです。
図２にて、CF は点７が点５から左にどれだけ進むかを示しています。CF は最大でも 100%
（点９の上）です。CF 入力値は 50〜90％程が適当です。実機と照合してください。
以下に、二段圧縮の計算プログラムの作り方を説明します。
− 18 −

A２.１

新しいシートを挿入

新しいシートに下図のように文字および罫線を書いてください。赤枠部分には後で計算式を書きます。データ入力セルには黄色ハッチングをします。入力可能はこのセルだけです。

ここに関数や数式を書き込みます

数式を書き込みます
数式を書き込みます

A２.２

関数と数式の書込み

各セルに書き込む数式は下表の通りです。青の式は該当するセルのコピーでも作れます。
A

B

C

D

E

F

G

H

②

計算結果

圧力

飽和温度

温度

比体積

比ｴﾝﾀﾙﾋﾟｰ比ｴﾝﾄﾛﾋﾟｰ

ｋJ/kg

KJ/(kg･K)=HVT(F10,=SVT(F10,

8
9

MPa

10 １

℃

℃

ｍ3/kg

低圧段圧縮機=C4=TsP(C10)=D10+G6=VPT(C10,

低圧段圧縮機=C5=TsP(C11)=TPS(C11,=VPT(C11,

高圧段圧縮機=C5=TsP(C12)

吸込
11

吐出

−273.15

12 ３
吸込

高圧段圧縮機=C6

吐出=TsP(C13)

H11)−273.15 E11+273.15)

E11+273.15)=HVT(F12,=TPS(C13,−273.15
凝縮器出口=C6=TsP(C14)

主膨張弁入口=C6=TsP(C15)=VPT(C13,=H13=HVT(F13,=HfPT(E14− 19 −=SVT(F12,E12+273.15)=D14−G5−273.15=H10

E12+273.15) E12+273.15)
E13+273.15)=TfH(H15)−273.15=VPT(C12,

H13)−273.15 E13+273.15)−273.1515 ７=HVT(F11,=D12+G5−273.1513 ４
14 ５

E10+273.15) E10+273.15) E10+273.15)−273.15

２=VfT(E14+273.15)+273.15=VfPT(C15,=G11−F18*

E15+273.15)(G12−G14)

E11 セルおよび E13 セルでは、比エントロピーから圧縮機吐出状態での温度を求めています。
E15 セルの TfH(H)はあまり使わないので、本文５章折込関数の項では紹介していませんが、液体の比エンタルピーから飽和液体温度を求める関数です（付録 A3 章を参照してください）。
G15 セルは二段圧縮一段膨張の場合に固有の数式です。
③計算結果まとめ
17
18

冷媒流量比

の各セルには、次の式を書き込んでください。
D
E
F
R min
R max
R 計算結果=(G11−G14)/(G12−G14) =(G11-HfPT(C12,D12=D18+G3*(E18−D18)+273.15)) /(G12−G14)

19

kJ/kg

D
高圧段動力

20

比ｴﾝﾀﾙﾋﾟｰ差= G13−G12

A２.３

E
低圧段動力
=G11−G10

F
合計動力

G
冷凍効果

H
COP

=F18*D20+F20 = G10−G17= G20/F20

冷凍サイクル計算の完成形

見やすくするために、罫線を書き加えてあります。有効数字は３〜４桁にしましたが、これで実用上十分と思います。必要ならばもう１桁加えてください。

（参考）このシートは R 290 のものですが、他の冷媒でも同様です。
なお、数式や文字が記入されたセルを保護するには、入力セル(黄色)のロックを外してから、シートの保護をかけます。詳細は Excel のヘルプを参照ください。
− 20 −

A２.４ ｐ−ｈ線図への書き込み準備
二段圧縮の理論サイクルを、右図の外殻の線（A 線）と高圧段側の膨張蒸発（B 線）とに分割して記入します。一筆書きもできます。興味があれば試してみてください。

圧力 (MPa)

Pk

７

Po

４

６
８

３
A線

る点数を追加します。プログラムに折込んだ二段圧縮は点数を追加してあり、線図へ

ｖ3

B線

Pm

このサンプルでは、圧縮行程は直線です。
等エントロピー線の曲線にするには計算す

５

２

ｖ１

１
比エンタルピ-、 h (kJ/kg)

二段圧縮理論冷凍サイクルの書込みデータも別シートで作っています。
圧力軸（縦軸）の対数は自然対数（Excel で LN、VBA で LOG）です。

まず、冷凍サイクル各線のデータ準備から始めます。Sheet2 に下図のようにデータ表を追加してください。各セルに文字および数式を書いてデータ表を作ります。
参照先をクリックして数式を書けば簡単です。赤い矢印はその時の参照先を示しています。
示したのは A 線（外郭の線）の参照です。他は倣ってください。
セルに数式が入力されると直ちに計算が行われ、そのセルに計算結果が表示されます。

− 21 −

A２.５ ｐ−ｈ線図への書込み（A 線）
準備が終ったので、これからｐ−ｈ線図上に理論冷凍サイクルを書き込みます。Excel の図
形描画のテクニックを多く使いますので、慣れていない方は注意して進めてください。
まず、シートを「線図」に切替えてください。「線図」に一段圧縮のサイクルが書かれていたら消してください。後で消しても OK です。消し方は本文 3.2）項を見てください。
２番目にシートの保護を外します。［ツール(T)］→［保護(P)］→［シート保護の解除(P)］の順にクリックすれば、保護は解除されます。
３番目に元のデータを呼び出します。線図のどこかをクリックして図を選択します。図の４辺にハンドル ■ が付けば選択されています。そのまま［グラフ(C)］→［元のデータ(S)］の順にクリックすれば、元のデータのダイアログが現れます（次項で説明）。
2.5.1 「元のデータ」ダイアログ（A 線のデータ入力と X の値の入力準備）

このダイアログでは、ｐ−ｈ線図の比エンタルピー軸（横軸）が「X の値」、圧力軸（縦軸）が「Y の値」と表現されています。以下の順序で操作してください。
１）タブ「系列」をクリックします。
２）「追加(A)」をクリックします。すると、系列７が追加されます。但し、系列の番号が７まであるかどうかは状況によります。
３）系列７のままでは不便なので、「名前(N)」のボックスをクリックして、A 線と記入します。系列７の名前もA 線 に変ります。
４）「X の値(X)」の右端のマークをクリックしてください。ダイアログが下図のように小さくなり、シートの切替えやデータ入力がより簡単にできるようになります。

− 22 −

A2.5.2

X の値（比エンタルピー）の入力

Sheet2 に切替え、B24 から H24 までドラッグすると X の値が入力できます。
下図は X が入力された後の状態を表示しています。小さくなったダイアログが、入力する
位置に重なってしまったときは、ダイアログを適宜動かしてください。

A2.5.3

Y の値の入力準備
前項（A2.5.2）のダイアログにて、右端のマークをクリックすると、基本ダイアロ
グ（下図）に戻ります。
「Y の値(Y)」右端のマークをクリックしてください。ダイアログは一番下の図のように小さくなり、シートの切替えが可能になります。

前の操作の結果を確認できます。
A線
X の値が入っているセル− 23 −

A2.5.4

Y の値（LOG（圧力））の入力

Sheet2 に切替えて、B25 から H25 までドラッグすると Y の値が入力されます。
ここでは背景（Sheet2 のこと）込みで説明しています。この右端のマークをクリックすると、「元のデータ」のダイアログに戻ります。

最後に OK をクリックす
ると、「線図」シートに戻ります。シートには理論冷凍サイク
ルが書かれています。但し、その図は下図の様に歪んでいる
場合があります。

ハンドル

この線は理論冷凍サイクルです。ただし、歪んだ曲線になっており修正が必要です。修正す
るために線上をクリックします。（詳細は次項で）− 24 −

A2.5.5

A 線（外郭線）の書式設定

歪んでいる A 線の上をクリックして線を選択します。きちんとクリックしないと他の線が選択されるので注意してください。選択された線上に表示されるハンドル ■ を確認して、目的の線か否かを確認してください。なお、線が非表示の場合もあります。その時はハンドルだ
けが頼りです。
そこで、［書式(O)］→［選択したデータの系列(E)］の順にクリックすると、次のダイアロ
グが出ます。以下、順々に操作してください。

１）［パターン］を選ぶ。
２）色(C)は自動から、青に変える（又は好みの色）
３）太さ(W)は更に太めに変える（又は好みに）
４）スムージング(M)のチェ
ックは外す。
歪んだ曲線になった原因はこのチェックが有りのためです。
５）OK をクリックすればダイアログが消えて完成です。完成した A 線（外郭線）を下図に示します。
６）サイクル上の各点に記号をつけたければ、マーカーの各項目で行います（説明は略します）。

− 25 −

A２.６ ｐ−ｈ線図への書込み（高圧段側の膨張・蒸発：B 線）

B 線の書き方は、A 線（外郭線）と同様です。ここでは操作のエッセンスと、結果のｐ−ｈ線図を説明します。詳細は A2.５の説明を見てください。
１）「元のデータ」を呼び出す。
(５項 参照)２）「追加(A)」をクリックする。
３）B 線

と書き込む

（左図はこの状態）
４）Sheet2 のセル B27 から D27を入力する。
５）Sheet2 のセル B28 から D28を入力する。
６）OK すれば線ができる。
７）書式設定で線の形を整える。

Ａ２.７

二段圧縮冷凍サイクルのまとめ

以上で二段圧縮の計算とそのｐ−ｈ線図への書込みが完成します。上記方法を応用すれば種々の用途にも使えると思います。独自の冷凍サイクル開発などに応用してください。

− 26 −

A3．折込関数（Function）一覧

Function 一覧表 （３）
R 744

R 717

R 290

R 600

R 600a

CO2

NH3

プロパン

ブタン

イソブタン

132.36405.51
11.36114
224.500.0044543
1120.53.5571
195.5
0.00609
0.064196.675369.825
4.24709218.5
0.0045767
558.8
2.061985.48

151.975425.1253.796
227.840.0043890
694.9
2.3652134.87134.667
407.8173.640
224.360.0044571
636.4
2.2327113.56

PsT(T)
TsP(P)TsV(V)
TsH1(H)
TsH2(H)PsT(T)
TsP(P)TsV(V)
TsH1(H)
TsH2(H)

30.9782tc℃
304.1282
Tc
K
7.3773pc MPa467.60ρｃ kg/m3
vｃ m3/kg 0.0021386-174.53
hｃ kJ/kg
sｃ kJ/kg･K -1.3054
216.592
K
Tt
0.51795
Pt MPaρt" kg/m3 13.7614

S1
S2
S3
S4 -1
-2

PsT(T)
TsP(P)TsV(V)
TsH1(H)
TsH2(H)PsT(T)
TsP(P)TsV(V)
TsH1(H)
TsH2(H)PsT(T)
TsP(P)TsV(V)
TsH1(H)
TsH2(H)-24.5℃
49.2℃
73℃
136℃
119.1℃h s max -69.72 kJ/kg 1491.1 kJ/kg 631.2 kJ/kg 748.6 kJ/kg 687.2 kJ/kg

TsS(S)TsS(S)

原著文献に示された定数値

三重点

S4 境界

S5 -1
-2
-3
G1
G2
G3
G4
G5
G6
G7
G8
G9
L1
L2
L3
L4
L5
L1 -1
L2 -1
L3 -1
L5 -1

臨界点

TsS(S)TsS(S)TsS(S)

この関数は、飽和ガスの計算用です。（気相線上）
気相線上のエンタルピ、エント
ロピから飽和温度を求めます。
ＴｓＨには、低温（TsH1）と高温(TsH2)の２種類があります。

ｐ

TsH2

液相線

PVT(V,T)
TPV(P,V)
VPT(P,T)
HPT(V,T)
TPH(P,H)
PHT(H,T)
SPT(V,T)TPS(P,S)

PVT(V,T)
TPV(P,V)
VPT(P,T)
HPT(V,T)
TPH(P,H)
PHT(H,T)
SPT(V,T)TPS(P,S)

PVT(V,T)
TPV(P,V)
VPT(P,T)
HPT(V,T)
TPH(P,H)
PHT(H,T)
SPT(V,T)TPS(P,S)

PVT(V,T)
TPV(P,V)
VPT(P,T)
HPT(V,T)
TPH(P,H)
PHT(H,T)
SPT(V,T)TPS(P,S)

PVT(V,T)
TPV(P,V)
VPT(P,T)
HPT(V,T)
TPH(P,H)
PHT(H,T)
SPT(V,T)TPS(P,S)

VfT(T)
HfT(T)
SfT(T)

VfT(T)
HfT(T)
SfT(T)

VfT(T)
HfT(T)
SfT(T)

VfT(T)
HfT(T)
SfT(T)

VfT(T)
HfT(T)
SfT(T)

飽和液の計算用です。

TfH(H)
VfPT(P,T)
HfPT(P,T)
SfPT(P,T)
TfPH(P,H)

TfH(H)
VfPT(P,T)
HfPT(P,T)
SfPT(P,T)
TfPH(P,H)

TfH(H)
VfPT(P,T)
HfPT(P,T)
SfPT(P,T)
TfPH(P,H)

TfH(H)
VfPT(P,T)
HfPT(P,T)
SfPT(P,T)
TfPH(P,H)

TfH(H)
VfPT(P,T)
HfPT(P,T)
SfPT(P,T)
TfPH(P,H)

圧縮液（過冷却液）の計算用。
超臨界のガス側にも使えます。
（但し境界は自ら確認ください。）− 27 −

TsH1

気相線

ｈ
ガス側の計算用です。

ｐ
圧縮液

ガス側

ｈ