境界要素法（BEM）による音響解析
境界要素法（BEM）による音響の理論基礎
音響BEMとは
先生、音響BEMって何ですか？
Helmholtz方程式を境界積分方程式に変換し、表面メッシュだけで音場を計算する手法。体積メッシュが不要なので、外部音場（無限領域）を自然に扱えるのが最大の利点。
支配方程式
音響場のHelmholtz方程式:
$k = \omega/c$: 波数。これをグリーンの定理で境界積分形に変換:
$G$: 自由空間グリーン関数:
$c(\mathbf{x})$: 境界上で$1/2$、領域内で$1$、領域外で$0$。
FEMと比べて何がいいんですか？
まとめ
BEMの起源は1960年代の弾性理論
境界要素法の数学的基盤を築いたのはスタンフォード大学のJaswon（1963年）とSympson。彼らは弾性問題の積分方程式を離散化する手法を発表したが、当時は音響への応用など想定していなかった。音響BEMへの転用はその10年後、Schenckが1968年にIBMの技術報告として発表したCHEFS法が先駆けとなり、後のNASTRANや市販ソルバーに影響を与えた。
境界要素法（BEM）による音響の数値計算手法
音響BEMの離散化
境界積分方程式をどうやって数値的に解くんですか？
表面を三角形 or 四角形要素で離散化。音圧$p$と法線速度$v_n$を節点値で近似。
$[H]$, $[G]$: 影響行列。各要素間のグリーン関数の積分値。
特異積分の処理
BEMの実装で最も技術的に難しいのが特異積分（$r \to 0$でグリーン関数が発散）の処理:
- 弱い特異性（$1/r$）: 極座標変換で正則化
- 強い特異性（$1/r^2$）: Cauchy主値で定義
- 超特異積分: Hadamard有限部分積分
非唯一性問題
外部BEMでは、閉じた境界の内部固有振動数で解が非唯一になる。対策:
- CHIEF法: 内部点に追加方程式を配置（過剰決定系）
- Burton-Miller法: 通常のBIEと法線微分BIEを線形結合。理論的に完全
Burton-Miller法が標準。超特異積分の処理が必要だが、確実に非唯一性を除去できる。
FMM（高速多重極法）
密行列で大規模問題が解けないのが弱点と言っていましたが…
FMM（Fast Multipole Method）で解決。遠方の要素群をまとめて近似計算:
- メモリ: $O(N^2) \to O(N)$
- 計算量: $O(N^2) \to O(N\log N)$
100万要素規模のBEMが実用的になった。Actran, COMSOL等に実装済み。
まとめ
FMMがBEMの計算量をO(N²)→O(N log N)に
従来のBEMは節点数Nに対してO(N²)のメモリと計算量が必要で、大規模モデルは事実上不可能だった。1987年にGreengardとRokhlinが発表した高速多重極法（Fast Multipole Method）が革命をもたらし、自動車ボディ規模（数十万節点）の音響解析も現実的になった。NuancesのVirtualLab AcousticsやLMS Sysnoise 5.xはFMM-BEMをいち早く実装した製品として知られる。
境界要素法（BEM）による音響の実務適用
音響BEMの実務
エンジン放射音、タイヤ騒音、変圧器騒音、排気管の音響放射が典型的な適用先。
解析フロー
1. 構造の振動解析 — FEMで表面振動速度$v_n$を取得
2. BEM表面メッシュ作成 — FEMメッシュから抽出 or 独立作成
3. 境界条件設定 — $v_n$（Neumann条件）を設定
4. BEM求解 — 表面音圧$p$を計算
5. 音場評価 — 任意の観測点での音圧、放射パワーを計算
実務チェックリスト
BEMメッシュのガイドライン
| 最高周波数 [Hz] | 音速340m/sでの波長 [m] | 要素サイズ上限 [mm] |
|---|---|---|
| 500 | 0.68 | 113 |
| 1000 | 0.34 | 57 |
| 2000 | 0.17 | 28 |
| 5000 | 0.068 | 11 |
| 10000 | 0.034 | 5.7 |
5kHz以上だとメッシュがかなり細かくなりますね。
そう。高周波ではBEMのメッシュコストが急増するので、SEAへの切り替えを検討する。
車室内音響はBEMとFEMの使い分けが鍵
自動車の車室内（キャビン）騒音解析では、1kHz以下の低周波はFEM、1kHz以上の高周波はBEMが有利とされる。Audi AG が2000年代初頭に発表したベンチマーク研究では、BEM-FEM連成で実験値と3dB以内の精度を1500Hz帯域まで達成した事例が報告されている。ただし実務では境界条件の設定（特に吸音材の複素インピーダンス値）が精度を左右する最大の課題だ。
境界要素法（BEM）による音響のソフトウェア比較
ツール
Actranが圧倒的な印象ですが…
自動車NVHではActranがデファクト。航空宇宙でも広く使われている。COMSOLは小規模な研究・検証用途に便利。Siemens系はNastranユーザーとの親和性が高い。
Sysnoise買収がActranの音響BEM機能を強化
音響BEMの老舗ソルバーLMS Sysnoise（ベルギーLMS International製）は2000年代に業界標準の地位を確立したが、2012年にSiemensがLMSを買収しSimcenter Testlabブランドへ統合された。一方、FFT社（現Siemens傘下）のActranは音響有限要素と境界要素のハイブリッドで差別化を図り、航空機エンジンナセルの透過損失解析に強みを持つ。OptiStructやNASTRANもBEMモジュールを後年追加している。
境界要素法（BEM）による音響の先端研究
先端技術
IGA-BEMってすごいですね。メッシュ不要になるんですか？
CADのNURBS表面をそのまま使うので、「メッシュ生成」という工程がなくなる。特に曲面が多い形状（自動車ボディ、航空機胴体）で威力を発揮する。ただしまだ研究段階で、商用実装は限定的だ。
水中音響BEMは潜水艦設計で機密扱いだった
水中放射音の境界要素解析は1970〜80年代の米海軍DARPA プログラムの中で急速に発展した。潜水艦のスクリュー放射音低減は軍事機密だったため、民間への技術移転は冷戦終結後の1990年代まで遅れた。現在はABAQUS/Acoustics や Actran が水中音響BEMを標準機能として提供しており、洋上風力発電基礎の水中騒音評価にも使われている。
境界要素法（BEM）による音響のトラブル対応
音響BEMのトラブル
| 症状 | 原因 | 対策 |
|---|---|---|
| 特定周波数で音圧が発散 | 非唯一性（内部共鳴） | Burton-Miller法を有効化 |
| 音圧結果が全面で反転 | 法線方向が内向き | 全要素の法線を外向きに修正 |
| 結果が振動解析と不整合 | FEM→BEMの速度マッピングミス | メッシュ対応関係を確認。補間精度を検証 |
| 計算が遅い/メモリ不足 | 密行列のまま解いている | FMMを有効化。メッシュ数を確認 |
| 近距離場の音圧が不正確 | 観測点が要素面上 or 極近傍 | 観測点を要素サイズの2〜3倍離す |
| 高周波で結果が不安定 | メッシュが粗い（6要素/波長未満） | メッシュ細分化。またはSEAに切り替え |
法線方向を間違えるとそんなに影響があるんですか？
法線が反転すると、境界積分方程式の符号が変わる。音圧の正負が逆転したり、Burton-Miller法が正しく機能しなくなる。BEMモデルの最初の確認事項として法線チェックを必ず行うこと。
CHIEF法がなければBEMは固有振動数で破綻する
音響BEMには「不正則振動数（irregular frequency）」問題がある。構造物の特定の固有振動数でBEM方程式が特異行列になり、解が発散する現象だ。Schenck（1968年）が提案したCHIEF（Combined Helmholtz Integral Equation Formulation）法は内部点に追加方程式を配置することで問題を回避し、現在も多くの商用ソルバーの標準手法となっている。CHIEF点の配置を誤ると再び不安定になる点に注意が必要だ。
なった
詳しく
報告