音響放射パワー

Q: 先生、音響放射パワーって何ですか？

振動する構造表面が周囲の空気を振動させて音を出す。その音のパワーを定量化する指標だ。 $$ W = \int_S \frac{1}{2}\text{Re}(p \cdot v_n^*)\,dS $$ $p$: 表面の音圧、$v_n$: 表面の法線方向振動速度、$*$: 複素共役。表面全体で音響インテンシティを積分したもの。

Q: 振動していれば必ず音が出るんですか？

振動しても音にならないことがある。それを表すのが放射効率 $\sigma_{rad}$。 $$ \sigma_{rad} = \frac{W}{\rho c S \langle v_n^2 \rangle} $$ $\rho c$: 空気の特性インピーダンス、$S$: 放射面積、$\langle v_n^2 \rangle$: 面平均振動速度二乗値。 $\sigma_{rad} = 1$: 完全放射体（ピストン振動） $\sigma_{rad} $\sigma_{rad} > 1$: まれだが、共振時に発生することがある

Q: 音響短絡って何ですか？

隣り合う領域が逆位相で振動すると、片方の正圧がもう片方の負圧で打ち消される。板のモード振動で顕著で、低周波ほど放射効率が低い理由だ。

カテゴリ: 構造解析 | 統合版 2026-04-06

CAE visualization for acoustic radiation theory - technical simulation diagram

音響放射パワー

音響放射パワーの理論基礎

音響放射とは

🧑‍🎓

先生、音響放射パワーって何ですか？

🎓

振動する構造表面が周囲の空気を振動させて音を出す。その音のパワーを定量化する指標だ。

$$ W = \int_S \frac{1}{2}\text{Re}(p \cdot v_n^*)\,dS $$

$p$: 表面の音圧、$v_n$: 表面の法線方向振動速度、$*$: 複素共役。表面全体で音響インテンシティを積分したもの。

放射効率

🧑‍🎓

振動していれば必ず音が出るんですか？

🎓

振動しても音にならないことがある。それを表すのが放射効率 $\sigma_{rad}$。

$$ \sigma_{rad} = \frac{W}{\rho c S \langle v_n^2 \rangle} $$

$\rho c$: 空気の特性インピーダンス、$S$: 放射面積、$\langle v_n^2 \rangle$: 面平均振動速度二乗値。

$\sigma_{rad} = 1$: 完全放射体（ピストン振動）
$\sigma_{rad} < 1$: 放射が弱い（音響短絡が発生）
$\sigma_{rad} > 1$: まれだが、共振時に発生することがある

🧑‍🎓

音響短絡って何ですか？

🎓

隣り合う領域が逆位相で振動すると、片方の正圧がもう片方の負圧で打ち消される。板のモード振動で顕著で、低周波ほど放射効率が低い理由だ。

臨界周波数と放射

🎓

放射効率は周波数によって大きく変わる:

$f < f_c$（コインシデンス周波数以下）: $\sigma_{rad} \ll 1$。音響短絡が支配的
$f = f_c$: $\sigma_{rad}$が急上昇
$f > f_c$: $\sigma_{rad} \approx 1$。全面が効率的に放射

まとめ

🎓

音響放射パワー $W$ — 振動面が放出する音のエネルギー

放射効率 $\sigma_{rad}$ — 振動がどれだけ効率的に音になるか

音響短絡 — 低周波で逆位相振動が音を打ち消す

コインシデンス周波数以上で$\sigma_{rad} \approx 1$

Coffee Break よもやま話

放射効率の概念はRayleigh卿の1877年論文が起源

音響放射効率（radiation efficiency）の理論的基盤はLord Rayleighが1877年に著した「The Theory of Sound」第2巻に遡る。彼は無限平板の「ピストンモデル」を用いて放射インピーダンスを導出したが、有限構造への拡張は20世紀まで待たねばならなかった。Wallace（1972年）が有限矩形平板の放射効率を解析的に求め、これが現代の構造音響理論の礎になっている。

音響放射パワーの数値計算手法

音響放射の計算手法

🧑‍🎓

音響放射パワーをFEMで計算するにはどうすればいいですか？

🎓

主に3つのアプローチがある。

1. FEM連成解析

🎓

構造FEM + 音響FEMの連成。音響領域を有限要素でモデル化し、外部にPML（完全整合層）を配置。

利点: 近距離場も遠距離場も計算可能
欠点: 音響領域のメッシュが大きくなる

2. BEM（境界要素法）

🎓

構造表面のメッシュだけで外部音場を計算。無限領域をそのまま扱える。

$$ c(\mathbf{x})p(\mathbf{x}) = \int_S \left(G\frac{\partial p}{\partial n} - p\frac{\partial G}{\partial n}\right)dS $$

$G$: 自由空間グリーン関数。表面の振動速度を境界条件として、音圧場を計算。

3. Rayleigh積分（平面近似）

🎓

無限バッフルに取り付けられた振動面を仮定:

$$ p(\mathbf{x}) = \frac{i\omega\rho}{2\pi}\int_S \frac{v_n(\mathbf{y})e^{ikr}}{r}\,dS(\mathbf{y}) $$

平面パネルの放射パワーの迅速な推定に便利。計算コストはBEMより小さい。

まとめ

🎓

FEM連成: 閉空間内部の放射に有効

BEM: 外部放射の標準手法。メッシュは表面のみ

Rayleigh積分: 平面バッフル仮定。高速な初期評価

Coffee Break よもやま話

近接場音響ホログラフィは逆問題解析の傑作

Near-field Acoustic Holography（NAH）は1985年にMaynardらがJASAに発表した革新的手法で、音圧マイクロフォンアレイの測定から放射面の振動分布を再構成できる。NAHはフォードやBMWのNVH部門が1990年代に実験ベンチ装備として導入し、エンジンカバーや燃料タンクの寄与音源を非接触で特定した。現在では音響カメラ（例：gfai techのSoundCam）として小型化・商品化されている。

音響放射パワーの実務適用

放射パワー解析の実務

🎓

自動車のエンジン放射騒音、家電のモーター騒音、変圧器のハム音が典型的な適用例。

解析フロー

🎓

1. 構造の振動解析 — モーダル or 周波数応答で表面振動速度を取得

2. 放射パワー計算 — BEM or Rayleigh積分で音場を計算

3. 放射効率の評価 — $\sigma_{rad}$を周波数ごとにプロット

4. 音圧レベル — 観測点の音圧をdB(A)で評価

実務チェックリスト

🎓

[ ] 構造モデルの固有振動数が実験と合っているか

[ ] BEM表面メッシュが$\lambda_{min}/6$以下か

[ ] 観測点の距離が正しいか（通常1m）

[ ] 放射パワーの結果がdB(A)で正しく重み付けされているか

[ ] エッジ部（音響短絡が起きやすい箇所）のメッシュが十分か

放射パワーレベル（SWL）

🎓

$$ L_W = 10\log_{10}\frac{W}{W_0} \quad [\text{dB}] $$

$W_0 = 10^{-12}$ W（基準音響パワー）。$L_W = 80$ dBなら$W = 10^{-4}$ W（0.1 mW）。

Coffee Break よもやま話

電動モーターの放射音は電磁力が支配的

EVモーターの放射騒音は燃焼エンジンと異なり、電磁トルクリップルによる電磁力が主因となる。トヨタ初代プリウス（1997年）開発時、エンジン停止時のモーター音が車室内に明瞭に聞こえることが設計後期まで未認識だった。現在はMaxwell応力テンソルを用いた電磁-構造-音響の3連成解析がデファクトとなり、JMAG＋ABAQUS＋Actranの組み合わせが国内大手サプライヤーで採用されている。

音響放射パワーのソフトウェア比較

ツール

🎓

ツール	手法	特徴
Actran（FFT/MSC）	FEM+BEM	音響放射の業界標準。FEM-BEM連成
Virtual.Lab（Siemens）	FEM+BEM	Nastranとのシームレス連携
LMS Sysnoise（Siemens）	BEM	音響BEMの古典。現在はSimcenter Acousticsに統合
COMSOL Acoustics	FEM+BEM	小〜中規模。マルチフィジックス連成
OpenBEM	BEM	オープンソース。教育・研究用

🧑‍🎓

どれを使えばいいですか？

🎓

外部放射（エンジン、モーター）: BEMベース → Actran or Virtual.Lab

内部放射（車室内）: FEMベース → Nastran + Actran

研究・教育: COMSOL or OpenBEM

Coffee Break よもやま話

Actranはエアバスのエンジン音響認証に使われる

MSC Actran（旧FFT Actran）はエアバスA380のエンジンナセル透過損失の型式認証解析に採用された実績を持つ。EASA（欧州航空安全機関）の騒音規制ICAO Annex 16 Chapter 14への適合確認にActranが活用され、ナセルライナーの吸音設計が最適化されている。国内では川崎重工・三菱重工がActranをMRJ（現SpaceJet）の機体放射音解析に使用したことを技術報告で公表している。

音響放射パワーの先端研究

先端技術

🎓

音響ホログラフィ（NAH） — 近距離場の音圧測定から振動面の速度分布を逆推定。音源同定に使用

IBEM（間接BEM） — 直接BEMでは困難な薄板構造の両面放射を効率的に計算

FMM-BEM（高速多重極BEM） — 従来BEMの$O(N^2)$を$O(N\log N)$に高速化。大規模モデル対応

AIによる放射パワー予測 — 構造振動パターンから放射パワーを即座に推定するサロゲートモデル

🧑‍🎓

FMM-BEMってそんなに速いんですか？

🎓

10万要素のBEMモデルで、従来法なら数日かかる計算が数時間で終わる。自動車の車体全体の音響放射解析が現実的になった。

Coffee Break よもやま話

音響メタマテリアルで放射音を操作する研究が加速

2000年代以降、局所共振型音響メタマテリアル（Liu et al., Science 2000）の登場で、構造物の放射音を能動的に制御する研究が加速した。MITのHunt研究室は2018年、厚さ1mm以下のメタ表面を鋼板に貼付するだけで特定周波数帯の放射音を15dB低減した実証実験を発表。商用化はまだ途上だが、航空機の胴体パネルへの適用研究がAirbusとともに進行中だ。

音響放射パワーのトラブル対応

音響放射のトラブル

🎓

症状	原因	対策
放射パワーがゼロ	振動速度の法線成分がゼロ or FSI界面未定義	速度出力の方向を確認。FSI界面を正しく設定
放射パワーが実測の10倍以上	BEMの内部共鳴（虚固有振動数）	CHIEF法 or Burton-Miller法で非唯一性を除去
特定周波数でスパイク	構造共振 + 高放射効率モード	制振処理で減衰追加。モード形状を確認
遠距離場の音圧が不正確	観測点がBEMメッシュに近すぎる	観測点をメッシュから十分離す（要素サイズの3倍以上）

🧑‍🎓

BEMの内部共鳴問題って何ですか？

🎓

外部BEMでは、閉じた表面の内部音響固有振動数で解が非唯一になる。物理的には意味がないが数値的にスパイクが出る。Burton-Miller法（法線微分方程式を重みつき結合）で解決するのが標準。

Coffee Break よもやま話

音圧測定とシミュレーションのズレは境界反射が原因

無響室が確保できない実環境で測定した放射音圧とシミュレーション結果が10dB以上乖離するケースが多い。主因は壁面・床面の反射音が測定値に混入することで、IEC 61672に準拠した測定でも低周波（200Hz以下）では反射の影響が無視できない。CAEとの比較には「自由音場補正」が必須であり、半無響室であれば床面の鏡像音源を追加したBEMモデルで補正する手法が確立されている。

音響放射パワー

音響放射パワーの理論基礎

音響放射とは

放射効率

臨界周波数と放射

まとめ

放射効率の概念はRayleigh卿の1877年論文が起源

音響放射パワーの数値計算手法

音響放射の計算手法

1. FEM連成解析

2. BEM（境界要素法）

3. Rayleigh積分（平面近似）

まとめ

近接場音響ホログラフィは逆問題解析の傑作

音響放射パワーの実務適用

放射パワー解析の実務

解析フロー

実務チェックリスト

放射パワーレベル（SWL）

電動モーターの放射音は電磁力が支配的

音響放射パワーのソフトウェア比較

ツール

Actranはエアバスのエンジン音響認証に使われる

音響放射パワーの先端研究

先端技術

音響メタマテリアルで放射音を操作する研究が加速

音響放射パワーのトラブル対応

音響放射のトラブル

音圧測定とシミュレーションのズレは境界反射が原因

関連トピック

関連する分野