地震時ダムの流体-構造連成 — 数値解法と実装
離散化手法
ダムと貯水池をどう離散化するんですか？
ダム堤体は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倍の水深相当の距離にモデル端を設定し、吸収境界を配置する方法もある。
リバティ船の溶接割れ——連成問題の教訓
第二次世界大戦中、アメリカは「リバティ船」を溶接で大量生産し、戦争の物流を支えました。しかし約1,500隻のうち約400隻に船体の亀裂が発生。原因は溶接残留応力と低温脆性の連成——溶接時の急激な温度変化が残留応力を生み、北大西洋の冷たい海水で鋼材が脆くなり、亀裂が伝播したのです。現代の溶接シミュレーションは、この「温度→残留応力→破壊」の連鎖を予測できます。
離散化手法の詳細解説
空間離散化における手法選択が数値精度・安定性・計算コストに与える影響を詳述する。
モノリシック法
全物理場を1つの連立方程式系として同時に解く。強い連成に対して安定だが、実装が複雑でメモリ消費が大きい。
パーティション法（分離反復法）
各物理場を独立に解き、界面でデータ交換。実装が容易で既存ソルバーを活用可能。弱い連成に適する。
界面データ転写
最近傍法（最も簡単だが精度低い）、射影法（保存的）、RBF補間（メッシュ非一致に強い）。保存性と精度のバランスが重要。
マトリクスソルバーの選定指針
問題規模と特性に応じた最適なソルバー選択のガイドライン。
| ソルバー種別 | 詳細・推奨条件 |
|---|---|
| ブロック対角前処理 | 各物理場の前処理を対角ブロックとして使用。実装が容易でスケーラビリティが良い。 |
| フィールド分割法 | 速度-圧力分割（流体）、変位-温度分割（熱-構造）等。物理的な意味に基づくブロック分割が安定性に寄与。 |
時間積分法と収束判定
ソルバー内部の制御パラメータと収束判定基準について記述する。
サブイタレーション
各連成ステップ内で十分な反復を行い、界面条件の整合性を確保。残差基準は各物理場の典型値に基づいてスケーリング。
Aitken緩和
連成反復の緩和係数を自動調整。過緩和による発散を防止し、収束を加速する適応的手法。
安定性条件
added mass効果（流体-構造連成で構造密度≈流体密度の場合）に注意。不安定な場合はロビン型界面条件やIQN-ILS法を適用。
数値解法の直感的理解
連成ソルバーのイメージ
モノリシック手法は「バンドの全員が同じ楽譜を見て同時に演奏する」——完璧な同期が取れるが、楽譜（連立方程式系）が巨大で複雑。パーティション手法は「各パートが別々に練習して合わせる」——個々の練習（ソルバー）は簡単だが、合わせ（データ交換）のタイミングと精度が課題。
Aitken緩和のたとえ
Aitken緩和は「シーソーのバランス取り」に似ている。一方が強く押しすぎると反対側が跳ね上がり、その反動でまた強く押しすぎる——この振動を抑えるために、押す力を自動的に調整するのがAitken緩和。連成反復が振動して収束しないとき、前回の修正量を見て次の修正量を自動調整する適応的手法。
連成解析の安定性やデータ転写の精度は、マルチフィジックスの永続的な課題です。 — Project NovaSolverはこの課題に正面から取り組んでいます。
地震時ダムの流体-構造連成の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →