水撃(ウォーターハンマー)解析 戻る
配管・流体工学

水撃(ウォーターハンマー)解析シミュレーター

管路長・内径・流速・弁閉鎖時間を調整して Joukowsky 式に基づき最大水撃圧・圧力波速度・臨界閉鎖時間をリアルタイム計算。圧力波形の変化と弁閉鎖時間の安全判定を可視化し、プラント配管の耐圧設計に活用できます。

管路・流体条件

— 判定中 —
最大水撃圧 ΔP (MPa)
圧力波速度 a (m/s)
臨界閉鎖時間 Tr (s)
静水圧 P₀ (MPa)
最大管内圧 Pmax (MPa)
tc / Tr 比

弁位置での圧力-時間波形(簡易矩形波近似)。赤点線 = 許容圧力(PN16=1.6MPa)目安。

弁閉鎖時間 tc に対する最大水撃圧 ΔP の変化。縦点線 = 現在の tc、横点線 = 臨界閉鎖時間 Tr。

管路長 L に対する臨界閉鎖時間 Tr と最大水撃圧(同一 V₀)の変化。

水撃(ウォーターハンマー)の仕組み — 会話で理解する

🧑‍🎓
「水撃」って水道の蛇口を急に閉めると「ドン!」っていう音がするやつですよね。でも工場配管で問題になるのはなぜですか?
🎓
そう、あれが水撃だ。家庭の蛇口なら音がする程度だけど、工場や発電所の大口径配管では流速も大きくて、弁を急閉鎖すると何MPaもの圧力衝撃波が発生することがある。これで配管の破断・継手の損傷・ポンプの逆転が起きる。過去に水力発電所や化学プラントで大事故が起きているんだ。
🧑‍🎓
Joukowsky の式 ΔP = ρaΔV、この「a」は何ですか?音速と同じですか?
🎓
水中の音速に近いけど、正確には「管内の圧力波速度」で、液体の体積弾性係数 K と管壁の弾性 E の両方が影響する。純水の a は約 1480 m/s、でも鋼管内では管壁の弾性で少し遅くなって 1200〜1400 m/s 程度になる。ゴム管や HDPE 管だとさらに遅くて 200〜500 m/s 程度になるため水撃圧が下がる。
🧑‍🎓
「臨界閉鎖時間」って何ですか?これより遅く閉めれば大丈夫なんですか?
🎓
Tr = 2L/a は圧力波が管路を往復する時間だ。弁閉鎖時間 tc > Tr なら「緩閉鎖」で、反射波が戻ってくる前に部分的に弁が開いたままなので水撃圧が下がる。tc < Tr だと「急閉鎖」で Joukowsky の最大値 ΔP = ρaΔV が発生する。例えば L = 100m、a = 1400m/s なら Tr = 0.14s。これより遅く閉めれば完全水撃は避けられる。
🧑‍🎓
水撃を防ぐには弁をゆっくり閉めるだけでいいんですか?他に対策はありますか?
🎓
緩閉鎖が最も基本的な対策だ。他にはサージタンク(管路に設けた開放タンクで圧力を吸収)、エアチャンバー(密閉空気タンクが圧力衝撃を緩和)、安全弁・リリーフバルブ(過圧時に逃がす)、フライホイール付きポンプ(急停止防止)なども使われる。複合的に対策することが多い。
🧑‍🎓
水撃の「負圧」って何ですか?水が「引っ張られる」んですか?
🎓
正圧波(ΔP > 0)の後に負圧波が来て、管内圧力が大気圧以下に下がることがある。これが「カラム分離(column separation)」を引き起こすと、水が気化・空洞化して、その後再び衝突するときに「二次水撃」として非常に大きな衝撃が発生する。長い管路や起伏のある配管で特に問題になる危険な現象だ。
🧑‍🎓
MOC(方法of Characteristics)解析ってよく聞きますが、このシミュレーターとどう違うんですか?
🎓
このツールは Joukowsky 式を使った「単純な1次元解析」だ。実際の設計ではMOC(特性曲線法)という数値解法を使って、管路の各区間での圧力・流速の時間変化を逐次計算する。分岐管・ポンプ・弁の特性曲線を組み込んで複雑な管網全体を精密にシミュレーションできる。AFT Impulse や HAMMER(Bentley)などの専用ソフトウェアがある。

理論メモ — Joukowsky の水撃式と圧力波速度

急閉鎖時(tc ≤ 2L/a)の最大水撃圧上昇(Joukowsky 式):

\[\Delta P = \rho \cdot a \cdot \Delta V\]

弾性管内の圧力波速度(Halliwell 式):

\[a = \sqrt{\frac{K/\rho}{1 + (K \cdot D)/(E \cdot e)}}\]

ここで K = 液体の体積弾性係数、E = 管材のヤング率、D = 内径、e = 管壁肉厚。

臨界(往復)時間と緩閉鎖時の水撃圧(Allievi 近似):

\[T_r = \frac{2L}{a}, \quad \Delta P_{slow} \approx \Delta P_{Joukowsky} \cdot \frac{T_r}{t_c} \quad (t_c > T_r)\]

よくある質問

臨界閉鎖時間 Tr = 2L/a です。tc ≤ Tr で急閉鎖(最大 ΔP = ρaΔV 発生)、tc > Tr で緩閉鎖(ΔP が 1/Tr×tc 倍に低減)。L=200m、a=1400m/s なら Tr=0.29s。これより遅く弁を閉めれば最大水撃圧を抑制できます。実設計では Tr の 2〜3 倍以上の閉鎖時間を目安とすることが多いです。
鋼管(E=210GPa): 1200〜1400 m/s、鋳鉄管: 1000〜1200 m/s、コンクリート管: 1000〜1200 m/s、HDPE 管(E=1GPa): 200〜400 m/s、ゴム管: 100〜300 m/s。弾性の低い材料ほど a が小さく、Joukowsky 式から水撃圧も下がります。HDPE や可撓管を使う水撃低減設計もあります。
負圧波が大気圧以下(絶対圧で液体の蒸気圧以下)になったとき、液体が気化して空洞(vapor cavity)が形成される現象です。その後空洞が崩壊するときに「二次水撃」として非常に大きな衝撃圧が発生します。長い管路、起伏のある管路、高速ポンプ停止時に特に危険です。MOC 解析での空洞モデルによるシミュレーションが必要です。
ポンプが急停止すると逆流防止弁が閉じることで弁閉鎖水撃と同様の現象が起きます。ΔV はポンプ定格流速と同程度になります。対策としてフライホイール(慣性力で停止を緩やかにする)、緩閉逆止弁、バイパス配管などが使われます。ポンプ保護のためにチェックバルブの動特性も MOC 解析に組み込む必要があります。
管路に設置した密閉容器に空気を封入したものです。圧力波が来ると空気が圧縮されてエネルギーを吸収し、圧力ピークを低減します。また負圧時には空気が膨張して液体を補給し、カラム分離を防止します。設計では空気容積・初期圧力・位置の最適化が必要で、圧縮性ガスの状態方程式(多変過程)を考慮します。
JISB 8265(圧力容器の構造)や JEAG 4601(原子力)では過渡荷重として水撃を考慮します。上水道ではWSP(日本水道協会規格)で弁操作時間に制限が設けられています。石油プラントでは API 651 等の国際規格、発電所では ASME B31.1 が参照されます。設計では MOC 解析結果の最大圧力に安全率をかけてフランジ定格圧力(ANSI PN)を選定します。