音響モーダル解析
理論と物理
音響モーダル解析とは
先生、音響モーダル解析って構造の振動ではなく「空気の振動」を解析するんですか？
そう。閉空間（車室、部屋、ダクト等）内の空気の固有振動数とモード形状を求める。構造の固有振動が「骨の振動」なら、音響モーダルは「肉（空気）の振動」だ。
支配方程式
音響場のHelmholtz方程式：
$p$ は音圧、$c$ は音速。FEMで離散化すると：
構造の固有値問題 $([K] - \omega^2 [M])\{u\} = \{0\}$ と全く同じ形！
そう。未知数が変位 $u$ ではなく音圧 $p$ に変わるだけ。同じLanczos法で解ける。
車室の音響モード
自動車の車室内の代表的な音響モード：
| モード | 振動数 | 特徴 |
|---|---|---|
| 1次縦モード | 80〜120 Hz | 前後方向の定在波 |
| 1次横モード | 200〜300 Hz | 左右方向 |
| 1次高さモード | 300〜400 Hz | 上下方向 |
1次縦モードが80 Hzって、エンジンの回転振動と近いですよね。
4気筒エンジンの2次（主要次数）が3000 rpmで100 Hz。車室の1次音響モードと共振する可能性がある。これがブーミングノイズの原因。NVH設計でこの共振を回避するのが重要。
まとめ
音響モーダル解析を整理します。
要点：
- 閉空間の空気の固有振動 — Helmholtz方程式の固有値問題
- 構造の固有値問題と同じ形 — Lanczos法で解ける
- 車室のブーミングノイズ — 音響モードと構造振動の共振
- NVH設計の基礎 — 音響モードの把握が騒音対策の第一歩
コンサートホールの残響設計
1900年に音響学の父サビン（ハーバード大）は残響時間T=0.161V/Aの式を確立した。ボストンシンフォニーホールの設計にこの式が使われ、1.8秒の残響時間が実現した。現在FEMで室内音響モードを計算すると20Hz〜5000Hzの固有モードが数千個現れ、サビン式の統計的正確性を個別モードレベルで確認できる。
各項の物理的意味
- 慣性項（質量項）：$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか？ あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
- 剛性項（弾性復元力）：$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね？ あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか？ 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い：「剛性が高い＝強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
- 外力項（荷重項）：体積力 $f_b$（重力など）と表面力 $f_s$（圧力、接触力など）。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」（体積力）、タイヤが路面を押す力は「表面だけにかかる力」（表面力）。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗：荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
- 減衰項：レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか？ いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら？ 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 変位 $u$ | m（メートル） | mm入力時は荷重・弾性率もMPa/N系に統一すること |
| 応力 $\sigma$ | Pa（パスカル）= N/m² | MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意 |
| 歪み $\varepsilon$ | 無次元（m/m） | 工学歪みと対数歪みの区別に注意（大変形時） |
| 弾性率 $E$ | Pa | 鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意 |
| 密度 $\rho$ | kg/m³ | mm系ではtonne/mm³（= 10⁻⁹ tonne/mm³ for 鋼） |
| 力 $F$ | N（ニュートン） | mm系ではN、m系ではNで統一 |
数値解法と実装
FEMでの音響解析
音響場のFEM要素はどんなものですか？
音響要素は1自由度（音圧 $p$）の要素。形状は構造要素と同じ（四面体、六面体等）だが、節点変数が変位ではなく音圧。
| ソルバー | 音響要素 | 備考 |
|---|---|---|
| Nastran | CAERO（パネル法）or FLUID | 流体要素 |
| Abaqus | AC3D4, AC3D8 | 音響四面体/六面体 |
| Ansys | FLUID30, FLUID220 | 音響要素。構造連成対応 |
音響要素で車室内をメッシュするんですね。
車室の空間を音響要素で充填する。壁面（ボディパネル）は構造要素。構造と音響の界面で流体-構造連成（FSI）を定義する。
構造-音響連成
連成の固有値問題：
$[A]$ が構造-音響の連成行列。
構造の変位と音響の音圧が連成するんですね。
パネルが振動すると音場に音圧が発生し、音圧が構造に力を及ぼす。この双方向連成を解くことで、構造振動→車室内騒音の伝達が予測できる。
まとめ
音響モーダルの数値手法、整理します。
要点：
- 音響要素（音圧自由度）で閉空間をメッシュ — AC3D4/8（Abaqus）, FLUID30（Ansys）
- 構造-音響連成 — 界面で変位と音圧が連成
- 連成固有値問題 — 構造と音響の同時固有値
- NVH解析のコアツール — 車室内騒音の予測
FEM音響モード解析の境界条件
音響モード解析では空気をポテンシャル流体要素（圧力が未知数）でモデル化し、壁面を剛体として固定、開口部を自由端（P=0）とする。有限要素の最小メッシュサイズは最高評価周波数の波長λの1/6以下が必要で、1000Hzの空気中音波（λ=340mm）では要素サイズ57mm以下とする。
線形要素（1次要素）
節点間を線形補間。計算コストは低いが、応力の精度が低い。せん断ロッキングに注意（低減積分やB-bar法で緩和）。
2次要素（中間節点付き）
曲線的な変形を表現可能。応力精度が大幅に向上するが、自由度は約2〜3倍に増加。推奨：応力評価が重要な場合。
完全積分 vs 低減積分
完全積分：過剰拘束（ロッキング）のリスク。低減積分：アワーグラスモード（零エネルギーモード）のリスク。適材適所で選択。
アダプティブメッシュ
誤差指標（ZZ推定量等）に基づく自動細分化。応力集中部の精度を効率的に向上。h法（要素分割）とp法（次数増加）がある。
ニュートン・ラフソン法
非線形解析の標準的手法。接線剛性マトリクスを毎反復更新。収束半径内で2次収束するが、計算コストが高い。
修正ニュートン・ラフソン法
接線剛性マトリクスを初期値または数反復毎に更新。各反復のコストは低いが、収束速度は線形的。
収束判定基準
力の残差ノルム: $||R|| / ||F_{ext}|| < \epsilon$（一般に $\epsilon = 10^{-3}$〜$10^{-6}$）。変位増分ノルム: $||\Delta u|| / ||u|| < \epsilon$。エネルギーノルム: $\Delta u \cdot R < \epsilon$
荷重増分法
全荷重を一度に負荷せず、小刻みに増加させる。弧長法（Riks法）は荷重-変位関係の極値点を越えて追跡可能。
直接法 vs 反復法のたとえ
直接法は「連立方程式を筆算で正確に解く」方法——確実だが大規模問題では時間がかかりすぎる。反復法は「当て推量を繰り返して正解に近づく」方法——最初は大雑把な答えだが、反復するたびに精度が上がる。辞書で言葉を探すとき、最初のページから順番に探す（直接法）より、見当をつけて開き、前後に調整する（反復法）方が効率的なのと同じ原理。
メッシュの次数と精度の関係
1次要素は「定規で曲線を近似する」——直線の折れ線で表現するため精度に限界がある。2次要素は「フレキシブルカーブ」——曲線的な変化を表現でき、同じメッシュ密度でも格段に精度が向上する。ただし、1要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
音響モーダルの実務
音響モーダルは実務でどう使いますか？
自動車のNVH開発が最大の適用。
メッシュの要件
音響要素のサイズは波長の1/6以下（二次要素）。$f_{max} = 500$ Hzの場合：
要素サイズ: $0.68 / 6 \approx 0.11$ m = 110 mm。
音響の要素は構造より粗いんですね。
音波の波長は構造の弾性波の波長より長いから、メッシュは粗くて済む。ただし高周波（1000 Hz以上）では細かいメッシュが必要。
実務チェックリスト
吸音材の効果もFEMで評価できるんですか。
吸音材はインピーダンス境界条件として音響面に設定する。吸音係数 $\alpha$ またはインピーダンス $Z$ をパラメータとして入力。
自動車室内音響モードのNVH解析
乗用車の室内音響空間（体積約3m³）は20〜200Hzに約100個の音響固有モードを持つ。フロアやルーフの構造共振と一致すると「ブーミング」が発生する。BMW 5シリーズでは1990年代末に室内音響FEMを設計プロセスに組み込み、試作段階でのNVH手直しを40%削減した。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
音響解析のツール
音響解析にはどんなツールが使えますか？
ActranやVirtual.Labが音響専用ツールなんですね。
選定ガイド
MSC Nastranの音響解析モジュール
MSC Nastranは構造と音響を統合した連成解析（SOL 103/108）でNVH解析の業界標準だ。BMW・Mercedes-Benzなど欧州自動車メーカーが内装パネルとキャビン音響の連成モード解析に採用。1モデル200万DOFの固有値解析を数時間で完了できる。2019年にSol 108の並列化が強化され従来比3倍の高速化を達成した。
選定で最も重要な3つの問い
- 「何を解くか」：音響モーダル解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
音響解析の先端研究
音響解析の最前線を教えてください。
SEA（統計的エネルギー解析）
高周波（1000 Hz以上）ではFEMのメッシュが膨大になる。SEAはモード密度が高い領域のエネルギーの流れを統計的に予測する。FEMとSEAのハイブリッド法が中周波帯の標準。
メタマテリアルによる音響制御
音響メタマテリアル（人工的に設計された周期構造）でバンドギャップ（特定周波数の音波を遮断）を実現する研究。FEMでユニットセルの分散関係を計算し、バンドギャップを設計。
AIによるNVH最適化
構造-音響連成のFEM結果をニューラルネットワークで学習し、構造のパラメータから車室内騒音をリアルタイム予測。設計空間の探索を高速化。
まとめ
音響解析の先端研究、まとめます。
流体-構造連成音響固有値解析
音響-構造連成問題では音響域と構造域の固有値を同時に解く「coupled FEM」が必要だ。圧力pと変位uの非対称連成行列の固有値計算には特殊なソルバー（unsymmetric Lanczos）が必要で、純音響・純構造の固有値とは異なるhybridモードが新たに現れる。潜水艦の筐体振動に伴う放射音予測で実用化されている。
トラブルシューティング
音響モーダルのトラブル
音響モーダル解析でよくあるトラブルは？
音響モードと構造モードが連成しない
界面の連成定義が間違っている。確認：
- 構造面と音響面がTIE / *TIE / FSI Interface で正しく接続されているか
- 法線方向が正しいか（音響面の法線が空間内側を向いていること）
音響振動数が理論値と合わない
矩形部屋の音響モード理論解：
FEMと理論値が合わない場合、音速 $c$ と密度 $\rho$ の設定を確認。
メッシュが粗すぎる
高周波の音響モードが出ない場合、メッシュが粗い。要素サイズ < $\lambda_{min}/6$ を確認。
まとめ
音響モーダルのトラブル対処、整理します。
音響モードのメッシュ細分化の目安
音響FEM解析で実測とモード周波数が5%以上ずれる場合、多くはメッシュが粗すぎる。評価周波数の波長に対してメッシュサイズλ/6より粗いと数値分散（数値的音速低下）が生じる。1000Hzまで評価するなら空気要素のサイズを50mm以下にすること。また硬壁での反射の際にインピーダンス境界条件の設定が漏れていないか確認する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——音響モーダル解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告