音響-構造相互作用 — 理論と支配方程式
概要
先生、音響-構造相互作用って、つまり音と構造が互いに影響し合うってことですか？
その通りだ。構造物が振動すると周囲の空気や流体を押して音波が発生する。逆に音圧が構造に荷重として作用して振動を引き起こす。この双方向のカップリングが音響-構造相互作用（ASI）だよ。
具体的にはどんな場面で問題になるんですか？
自動車の車室内騒音が代表例だね。エンジンの振動がボディパネルを励振し、パネルが車室内の空気を圧縮・膨張させて乗員が聞く音になる。他には航空機の胴体透過騒音、建築物の遮音設計、潜水艦のソナードーム設計などがある。
支配方程式
構造側と音響側、それぞれどんな方程式で記述するんですか？
構造側は通常の弾性体の運動方程式だ。
右辺の $[A]\{p\}$ が音圧から構造への荷重項で、$[A]$ は連成マトリクスだ。音響側はヘルムホルツ方程式の時間領域版、つまり波動方程式になる。
この2つはどこで繋がるんですか？
連成界面での境界条件で繋がる。音響側から見ると、界面の法線方向速度が構造の変位速度と一致する必要がある。
これを行列形式で書くと、連成系全体は非対称な方程式系になるんだ。
非対称なのが厄介そうですね。
そうなんだ。NastranではFLUID要素（CAABSF, CHEXAなど）を使うと内部でこの非対称連成系を組んでくれる。AbaqusではAcoustic-Structural連成のための*TIE制約が用意されているよ。
周波数領域の定式化
振動騒音って周波数で議論することが多いですよね？
調和応答を仮定すると $\{u\} = \{\hat{u}\}e^{i\omega t}$、$\{p\} = \{\hat{p}\}e^{i\omega t}$ として周波数領域に変換できる。
周波数ごとに解くので、周波数応答関数（FRF）を効率よく求められるんだ。
この定式化だと、どのモードがどの周波数で問題になるか見やすいですね。
まさにそう。特に構造の固有振動数と音響キャビティの固有振動数が近い場合、連成共振が起きて音圧が増幅される。車室内のブーミング問題がこの典型だよ。
実務上のポイント
解析する上で気をつけるべきことは何ですか？
まず音響メッシュの要素サイズだ。最高解析周波数 $f_{max}$ に対して、要素あたり最低6節点が波長内に入るようにする。
空気中で $c = 343$ m/s、$f_{max} = 1000$ Hz なら $h \leq 57$ mm になる。二次要素ならもう少し粗くてもよいが、一次要素では守った方がいい。
なるほど、低次要素だとかなり細かいメッシュが必要になりますね。全体像がだいぶ見えてきました。
心臓シミュレーション——究極のFSI問題
人間の心臓は1日に約10万回拍動し、血液を全身に送り出します。この過程は流体（血液）-構造（心筋・弁）-電気（刺激伝導系）の3場連成問題。心臓のデジタルツインの構築は連成解析の「聖杯」と呼ばれ、世界中の研究者が挑戦しています。実現すれば、手術のシミュレーションや薬の効果予測が患者ごとにカスタマイズできるようになります。
各項の物理的意味
- 構造-熱連成項：温度変化による熱膨張が構造変形を誘発し、変形が温度場に影響する。$\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%以下が目安 |
数値例：熱-構造連成（アルミ板, α=23×10⁻⁶/K, 温度差ΔT=100K, 拘束あり）
自由熱膨張 ε=αΔT = 23×10⁻⁶×100 = 0.23% 完全拘束時の熱応力 σ=EαΔT = 70000×23×10⁻⁶×100 = 161 MPa
連成手法別の計算コスト比較（同一問題）：
片方向で済むなら計算コストは1/8！ ただし温度変形が流体領域を大きく変える場合は強連成が必須。「過剰品質」と「不足」の見極めが実務のカギです。
連成解析の安定性やデータ転写の精度は、マルチフィジックスの永続的な課題です。 — Project NovaSolverはこの課題に正面から取り組んでいます。
Project NovaSolver — CAE実務の課題に向き合う研究開発
「音響-構造相互作用をもっと効率的に解析できないか？」——私たちは実務者の声に耳を傾け、既存ワークフローの改善を目指す次世代CAEプロジェクトに取り組んでいます。具体的な機能はまだ公開前ですが、開発の進捗をお届けします。
進捗通知を受け取る →