地震時ダムの流体-構造連成
理論と物理
現象の概要
地震時にダムで流体-構造連成が重要になるのはなぜですか？
地震動でダム堤体が振動すると、背後の貯水池の水が動圧としてダム面に作用する。この動水圧がダムの応答を増大させ、場合によっては堤体の安全性を脅かす。Westergaard（1933）の古典的な研究以来、ダムの耐震設計では動水圧の考慮が必須だ。
支配方程式
貯水池の水はどう扱うんですか？
水を非粘性・非圧縮性流体と仮定するとラプラス方程式になるが、地震応答では圧縮性を考慮して波動方程式を使うことが多い。
$c$ は水中の音速（約1440 m/s）だ。ダム面での境界条件は、
$\ddot{u}_n$ はダム面の法線方向加速度だ。
ダム堤体は構造力学の運動方程式で記述する。
$\{F_{eq}\}$ は地震慣性力、$\{F_{hydro}\}$ は動水圧による荷重ベクトルだ。
Westergaardの付加質量法はまだ使われていますか？
簡易手法として現在も使われている。剛体ダムに対する動水圧分布を付加質量に換算する。
$H$ は水深、$y$ は水面からの深さだ。ただし、ダムの柔軟性や貯水池の有限長さの効果は考慮できないので、詳細評価にはFEM-BEM連成が必要だ。
ダムと水の「固有振動数対決」——地震波が引き起こす増幅の罠
コンクリートダムは一見「超剛体」に見えますが、貯水池の水と一体になると固有振動数が大きく変わります。たとえば空のダムと満水のダムでは、流体慣性の追加効果（付加質量）で固有振動数が20〜40%も低下することが理論的に示されています。1971年のサン・フェルナンド地震では、あるロックフィルダムが設計地震力を超える揺れを受けたにもかかわらず崩壊を免れましたが、後の調査で「満水状態が振動数を下げたため主要地震波の卓越周期とずれた」という皮肉な幸運が明らかになりました。理論を理解しないまま設計すると、満水か空かだけで耐震性評価がまるで変わってしまうんです。
各項の物理的意味
- 構造-熱連成項：温度変化による熱膨張が構造変形を誘発し、変形が温度場に影響する。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常の例】夏に線路のレールが伸びて隙間が狭くなる——温度上昇→熱膨張→応力発生の典型例。電子基板がはんだ付け後に反るのも、異なる材料の熱膨張率差による。エンジンのシリンダーブロックは高温部と低温部の温度差で熱応力が発生し、最悪の場合亀裂に至る。
- 流体-構造連成（FSI）項：流体圧力・せん断力が構造を変形させ、構造変形が流体領域を変化させる双方向の相互作用。【日常の例】強風で吊り橋のケーブルが振動する（渦励振）——風の力が構造を揺らし、揺れた構造が風の流れを変え、さらに振動が増幅する。心臓の血流と血管壁の弾性変形、航空機の翼のフラッタ（空力弾性不安定性）も典型的なFSI問題。片方向のみの連成で済む場合もあるが、変形が大きい場合は双方向連成が必須。
- 電磁-熱連成項：ジュール発熱 $Q = J^2/\sigma$ が温度上昇を引き起こし、温度変化が電気抵抗を変化させるフィードバックループ。【日常の例】電気ストーブのニクロム線は電流が流れると発熱（ジュール熱）して赤くなる——温度が上がると抵抗が変わり、電流分布も変化する。IHクッキングヒーターの渦電流発熱、送電線の温度上昇による弛み増加もこの連成の例。
- データ転写項：異なる物理場間のメッシュ不一致を補間で解決。【日常の例】天気予報で「気温のデータ」と「風のデータ」を合わせて体感温度を計算するとき、それぞれの観測地点が異なれば補間が必要——CAEの連成解析でも、構造メッシュとCFDメッシュは一般に一致しないため、界面でのデータ転写（補間）精度が結果の信頼性に直結する。
仮定条件と適用限界
- 弱連成仮定（片方向連成）：一方の物理場が他方に影響するが逆は無視可能な場合に有効
- 強連成が必要なケース：FSIでの大変形、電磁-熱連成での温度依存性が強い場合
- 時間スケールの分離：各物理場の特性時間が大きく異なる場合、サブサイクリングで効率化可能
- 界面条件の整合性：連成界面でのエネルギー・運動量保存が数値的に満たされることを確認
- 適用外ケース：3つ以上の物理場が同時に強く連成する場合、モノリシック手法が必要になることがある
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 熱膨張係数 $\alpha$ | 1/K | 鋼: 約12×10⁻⁶、アルミ: 約23×10⁻⁶ |
| 連成界面力 | N/m²（圧力）またはN（集中力） | 流体側と構造側で力の釣り合いを確認 |
| データ転写誤差 | 無次元（%） | 補間精度はメッシュ密度比に依存。5%以下が目安 |
数値解法と実装
ダム堤体はFEM（ソリッド要素）、貯水池は音響流体要素（acoustic element）で離散化するのが標準的だ。AbaqusのAC3D系要素やAnsysのFLUID30要素が使われる。
連成系の運動方程式は以下の形になる。
$R$ は連成マトリクスで、流体-構造界面での面積積分から構成される。
時間積分はどうするんですか？
Newmark-β法（$\beta = 0.25, \gamma = 0.5$）やHHT-α法が標準的だ。地震波の入力は加速度時刻歴を基盤面に設定する。サンプリング間隔は通常0.01秒だが、高周波成分を考慮する場合は0.005秒以下が必要だ。
貯水池の無限延長の処理
貯水池は上流側に無限に広がりますよね。どう処理するんですか？
ソンマーフェルト放射条件（無反射境界）を上流端に設定する。Ansysではimpedance boundary、AbaqusではNon-Reflecting Boundary Conditionとして実装されている。あるいは、ダム面から5〜10倍の水深相当の距離にモデル端を設定し、吸収境界を配置する方法もある。
ウェスターガード式の「過信」——単純公式が招く設計ミス
ダムの地震動水圧を計算するとき、今でも教科書に載っている「ウェスターガード式（1933年）」が使われることがあります。この式は直立平面ダムに鉛直方向地震加速度が作用すると仮定した解析解で、計算が非常に簡単です。しかし問題は、アーチダムや台形断面ダムには適用誤差が出ること、そして「ダム-水-岩盤の連成振動」を無視していること。実際、ある重力ダムの解析で、ウェスターガード式による動水圧はFEM連成解析の結果より30%以上過大評価していたケースがあります。「簡単な式だから安全側」と思いがちですが、過大評価は不必要な補強コストにつながる——数値解法の必要性はコスト削減にもあります。
モノリシック法
全物理場を1つの連立方程式系として同時に解く。強い連成に対して安定だが、実装が複雑でメモリ消費が大きい。
パーティション法（分離反復法）
各物理場を独立に解き、界面でデータ交換。実装が容易で既存ソルバーを活用可能。弱い連成に適する。
界面データ転写
最近傍法（最も簡単だが精度低い）、射影法（保存的）、RBF補間（メッシュ非一致に強い）。保存性と精度のバランスが重要。
サブイタレーション
各連成ステップ内で十分な反復を行い、界面条件の整合性を確保。残差基準は各物理場の典型値に基づいてスケーリング。
Aitken緩和
連成反復の緩和係数を自動調整。過緩和による発散を防止し、収束を加速する適応的手法。
安定性条件
added mass効果（流体-構造連成で構造密度≈流体密度の場合）に注意。不安定な場合はロビン型界面条件やIQN-ILS法を適用。
Aitken緩和のたとえ
Aitken緩和は「シーソーのバランス取り」に似ている。一方が強く押しすぎると反対側が跳ね上がり、その反動でまた強く押しすぎる——この振動を抑えるために、押す力を自動的に調整するのがAitken緩和。連成反復が振動して収束しないとき、前回の修正量を見て次の修正量を自動調整する適応的手法。
実践ガイド
典型的なアーチダムの場合、
1. 堤体のFEモデル作成（ソリッド要素、20節点六面体が望ましい）
2. 基礎岩盤のモデル化（ダム高の2〜3倍の範囲）
3. 貯水池の音響流体モデル（ダム面から上流方向に3〜5倍の水深距離）
4. 流体-構造界面の定義（tie constraint）
5. 地震入力波の設定（基盤面にuniform excitation or deconvolved motion）
基礎岩盤のモデル化が重要な理由は何ですか？
ダムの応答はダム-基礎-貯水池の3者連成で決まる。基礎岩盤の質量効果（inertia effect）と放射減衰（radiation damping）を無視すると応答を過大評価する。mass-less foundationの仮定は保守的だが、非現実的に大きな応力を予測する場合がある。
材料モデル
コンクリートダムの材料モデルは何を使うんですか？
線形弾性でスクリーニング解析を行い、必要に応じてコンクリート損傷塑性モデル（Concrete Damaged Plasticity: CDP）に切り替える。AbaqusのCDPモデルは圧縮・引張の損傷変数 $d_c, d_t$ で剛性低下を表現する。
| パラメータ | 典型値（マスコンクリート） |
|---|---|
| Young率 | 25〜35 GPa |
| Poisson比 | 0.18〜0.20 |
| 密度 | 2400 kg/m³ |
| 引張強度 | 2〜4 MPa |
| 圧縮強度 | 20〜40 MPa |
| 減衰比 | 5% |
現場技術者が語る「水位1m差の恐怖」——実務でのダムFSI
ダム管理の現場で「地震のときは水位を下げておけ」とよく言われますが、これには定量的な根拠があります。水位が設計満水位より10m下がるだけで、動水圧荷重は面積比で大幅に減少し、ダム底部の引張応力が15〜25%改善されるケースが報告されています。ただし「水を下げると下流への洪水リスクが消えるのでは？」と思いきや、急激な水位操作は吸い出し現象やパイピング（水みち形成）を誘発するリスクがある。実務では「地震前に段階的に水位を下げる」判断基準を事前に決めておくことが重要で、日本の主要ダムでは地震対応マニュアルに具体的な水位管理手順が盛り込まれています。
解析フローのたとえ
風船を膨らませたことがありますか？ あの瞬間、実は高度な流体-構造連成が起きています。内部の空気圧（流体）がゴム壁（構造）を押し広げ→広がった壁が内部の圧力分布を変え→変わった圧力がさらに壁を変形させる…このキャッチボールを計算ステップごとに繰り返すのがFSI解析です。
初心者が陥りやすい落とし穴
「片方向連成で十分でしょ？」——この判断ミスが連成解析で最も危険です。構造の変形が微小なら確かに片方向で足りますが、心臓弁の開閉のように変形が流路を大きく変える場合、片方向では全く話になりません。目安は「変形量が代表長さの1%を超えるか」。超えるなら双方向連成は必須です。片方向で済ませてしまった場合、結果が「もっともらしいけど実は大間違い」になる——これが最も怖いパターンです。
境界条件の考え方
連成界面のデータ交換は「国境の出入国管理」と同じです。各国（物理場）には独自の法律（支配方程式）がありますが、国境（界面）で人や物（力・温度・変位）のやり取りを正確に管理しないと、両国の経済（エネルギーバランス）が崩壊します。メッシュが一致していない場合の補間は「通訳」のようなもの——誤訳（補間誤差）が小さいほど良い結果が得られます。
ソフトウェア比較
ツール比較
ダムの地震FSI解析に使えるソフトウェアは何がありますか？
主要ツールを整理しよう。
| ツール | 流体モデル | 特徴 |
|---|---|---|
| Abaqus | 音響要素(AC3D) | CDPモデル対応。ダム解析の実績豊富 |
| Ansys Mechanical | FLUID30/220要素 | 音響-構造連成。大規模並列に強い |
| DIANA FEA | 音響要素 | コンクリート特化。Westergaard自動設定 |
| LS-DYNA | ALE流体 or SPH | 越流・破壊を含む極限解析向き |
| EACD-3D | BEM（音響） | 米陸工兵隊開発。ダム専用コード |
DIANA FEAはダム解析に特化しているんですか？
DIANA（Displacement Analyzer）はオランダTNO由来のソフトで、コンクリート構造の非線形解析に強みがある。アーチダムの3D解析でICOLDベンチマーク問題の参照解として広く使われている。
規格や基準への準拠はどうなっていますか？
ICOLD（国際大ダム会議）のBulletinやUSACE EM 1110-2-6051が参照基準だ。日本では国土交通省の河川構造物設計指針に従う。大規模ダムの耐震設計では動的FEM-FSI解析が標準になりつつある。
「音響要素 vs. ALE流体」——ダムFSIでソフト選びが難しい理由
ダムの地震FSI解析でよく議論になるのが「音響流体要素で十分か、それともNavier-Stokes全解析が必要か」という問いです。貯水池の水は通常、圧縮性が小さく流速も低いため、音響近似（音響要素）で地震応答は十分に捉えられます。ただし越流や漏水、スロッシング（水面の大波）を同時に評価したい場合はALE法やSPH法の出番になる。AbaqusとOpenFOAMを比較すると、Abaqusは音響要素が堅牢で減衰設定も豊富ですが、スロッシング計算はOpenFOAMが圧倒的に得意。「ダムのどの現象を見たいか」によってツールを使い分けるのが賢い選択で、欧米の大規模ダム評価プロジェクトでは両者を組み合わせた事例が増えています。
選定で最も重要な3つの問い
- 「何を解くか」：地震時ダムの流体-構造連成に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
亀裂内への水圧浸透は、亀裂開口幅に応じたcubic flow法則で記述される。
$w$ は亀裂開口幅だ。AbaqusのXFEMと流体要素の連成でこの種の解析が可能だ。
確率論的地震安全性評価
地震の不確かさをどう考慮するんですか？
Fragility Curve（フラジリティ曲線）を作成する。地震強度指標（PGA等）に対する損傷超過確率を表す曲線で、数十〜数百ケースの動的解析が必要だ。ラテン超方格サンプリングや代理モデル（Kriging等）で計算コストを削減する。
機械学習による応答予測
AIを使ったダム応答の予測研究もあるんですか？
LSTMやCNNで地震波入力からダム応答（頂部変位、最大応力）を直接予測する研究が発表されている。学習データにはFSI解析の結果を使い、数千ケースの地震波に対する応答データベースを構築する。リアルタイムのダム健全性モニタリングへの応用が期待されている。
ダム亀裂と水圧の「デスマッチ」——XFEM研究の最前線
大地震でダム堤体にひびが入ったとき、最も怖いのはそこに貯水池の水圧が染み込むことです。亀裂先端に水圧が作用すると応力拡大係数が跳ね上がり、亀裂が一気に進展する——いわゆる「水力破砕（ハイドロリックフラクチャリング）」と同じメカニズムが起きます。近年はXFEM（拡張有限要素法）で亀裂経路を追跡しながら、亀裂内部の水流と水圧を連成させる手法が研究されています。カナダのSteep Rock Lakeダムの事後解析では、この水圧連成を考慮したモデルが実際の亀裂パターンを95%の精度で再現できたと報告されており、先端技術が実ダムの評価に直結し始めています。
トラブルシューティング
音響要素の異常応答
貯水池の圧力に非物理的な振動が出るんですが。
よくある問題だ。原因と対策をまとめよう。
| 原因 | 症状 | 対策 |
|---|---|---|
| 無反射境界の不備 | 上流端からの反射波 | ダム面から5H以上離す。impedance boundary追加 |
| 自由表面条件の不整合 | 表面での圧力振動 | p=0境界条件の確認。重力予圧なしの設定確認 |
| メッシュ密度不足 | 高周波数成分の解像不良 | 対象最高周波数の1/6波長以下に |
| 時間刻みが粗い | 高周波ノイズ | $\Delta t < T_{min}/20$ |
基礎岩盤モデルの影響
基礎岩盤をmass-lessにするかどうかで結果が大きく変わるんですが、どちらが正しいんですか？
massless foundationは保守的（応答が大きくなる方向）だが、非現実的に大きな応答を予測することがある。USACE EM 1110-2-6051では放射減衰を含む基礎モデルの使用を推奨している。Lysmer-Kuhlemeyer粘性境界やPMLを基礎岩盤の外周に設定して放射減衰を表現する。
結果の検証指標はありますか？
ICOLDのベンチマーク問題（Theme A: Cotter Dam、Theme C: Morrow Point Dam等）の結果と比較するのが有効だ。堤頂部の加速度応答スペクトルや基本固有振動数が参照値と合っているかを確認する。
「反射波が戻ってきた！」——貯水池の無反射境界設定で詰まったら
ダム地震FSI解析でよく起きるトラブルが「貯水池上流端からの圧力波の反射」です。本来、水面から遠ざかる方向に伝わる波は無限遠に消えるはずですが、計算領域を有限にするとそこで反射してしまい、非物理的な振動が発生します。これを「反射波アーティファクト」と呼びます。ある実務プロジェクトで、上流端をただ「固定壁」にしてしまったため、ダムの応力波形が実測の3倍近くになるという大失敗があったそうです。対策はインピーダンス境界（吸収境界条件）の設定ですが、AbaqusとAnsysでは設定方法が全然違うので要注意。「圧力が振動していておかしい」と思ったら、まず境界条件を疑いましょう。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——地震時ダムの流体-構造連成の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告