音響モーダル解析
音響モーダルの理論基礎
音響モーダル解析とは
先生、音響モーダル解析って構造の振動ではなく「空気の振動」を解析するんですか？
そう。閉空間（車室、部屋、ダクト等）内の空気の固有振動数とモード形状を求める。構造の固有振動が「骨の振動」なら、音響モーダルは「肉（空気）の振動」だ。
支配方程式
音響場の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の固有モードが数千個現れ、サビン式の統計的正確性を個別モードレベルで確認できる。
音響モーダルの数値計算手法
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以下とする。
音響モーダルの実務適用
音響モーダルの実務
音響モーダルは実務でどう使いますか？
自動車の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%削減した。
音響モーダルのソフトウェア比較
音響解析のツール
音響解析にはどんなツールが使えますか？
ActranやVirtual.Labが音響専用ツールなんですね。
選定ガイド
MSC Nastranの音響解析モジュール
MSC Nastranは構造と音響を統合した連成解析（SOL 103/108）でNVH解析の業界標準だ。BMW・Mercedes-Benzなど欧州自動車メーカーが内装パネルとキャビン音響の連成モード解析に採用。1モデル200万DOFの固有値解析を数時間で完了できる。2019年にSol 108の並列化が強化され従来比3倍の高速化を達成した。
音響モーダルの先端研究
音響解析の先端研究
音響解析の最前線を教えてください。
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以下にすること。また硬壁での反射の際にインピーダンス境界条件の設定が漏れていないか確認する。
関連トピック
なった
詳しく
報告