地震時ダムの流体-構造連成
地震時ダムの流体-構造連成の理論基礎
現象の概要
地震時にダムで流体-構造連成が重要になるのはなぜですか？
地震動でダム堤体が振動すると、背後の貯水池の水が動圧としてダム面に作用する。この動水圧がダムの応答を増大させ、場合によっては堤体の安全性を脅かす。Westergaard（1933）の古典的な研究以来、ダムの耐震設計では動水圧の考慮が必須だ。
支配方程式
貯水池の水はどう扱うんですか？
水を非粘性・非圧縮性流体と仮定するとラプラス方程式になるが、地震応答では圧縮性を考慮して波動方程式を使うことが多い。
$c$ は水中の音速（約1440 m/s）だ。ダム面での境界条件は、
$\ddot{u}_n$ はダム面の法線方向加速度だ。
ダム堤体は構造力学の運動方程式で記述する。
$\{F_{eq}\}$ は地震慣性力、$\{F_{hydro}\}$ は動水圧による荷重ベクトルだ。
Westergaardの付加質量法はまだ使われていますか？
簡易手法として現在も使われている。剛体ダムに対する動水圧分布を付加質量に換算する。
$H$ は水深、$y$ は水面からの深さだ。ただし、ダムの柔軟性や貯水池の有限長さの効果は考慮できないので、詳細評価にはFEM-BEM連成が必要だ。
ダムと水の「固有振動数対決」——地震波が引き起こす増幅の罠
コンクリートダムは一見「超剛体」に見えますが、貯水池の水と一体になると固有振動数が大きく変わります。たとえば空のダムと満水のダムでは、流体慣性の追加効果（付加質量）で固有振動数が20〜40%も低下することが理論的に示されています。1971年のサン・フェルナンド地震では、あるロックフィルダムが設計地震力を超える揺れを受けたにもかかわらず崩壊を免れましたが、後の調査で「満水状態が振動数を下げたため主要地震波の卓越周期とずれた」という皮肉な幸運が明らかになりました。理論を理解しないまま設計すると、満水か空かだけで耐震性評価がまるで変わってしまうんです。
地震時ダムの流体-構造連成の数値計算手法
ダム堤体は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. 堤体の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解析に使えるソフトウェアは何がありますか？
主要ツールを整理しよう。
| ツール | 流体モデル | 特徴 |
|---|---|---|
| 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が圧倒的に得意。「ダムのどの現象を見たいか」によってツールを使い分けるのが賢い選択で、欧米の大規模ダム評価プロジェクトでは両者を組み合わせた事例が増えています。
地震時ダムの流体-構造連成の先端研究
亀裂内への水圧浸透は、亀裂開口幅に応じた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では設定方法が全然違うので要注意。「圧力が振動していておかしい」と思ったら、まず境界条件を疑いましょう。
関連トピック
なった
詳しく
報告