騒音予測シミュレーション — FEM-BEM連成からSEAまで

ざっくり言うと2段階だ。まず構造のFEM振動解析で表面速度分布を求める。次にBEM（境界要素法）で、その振動する表面から放射される音圧分布を計算する。

一番わかりやすいのは自動車のロードノイズだね。タイヤ→サスペンション→車体パネルと振動が伝わり、車室内に音が放射される。この「構造→音響」の連鎖を追うのがFEM-BEM連成の定番用途だ。ほかにも、コンプレッサーの騒音、家電のモーター音、建設機械の環境騒音なんかもこの手法で予測する。

概念はシンプルだけど、周波数帯域で手法を使い分ける必要がある。低周波（〜500 Hz）はFEM-BEM、中高周波（500 Hz〜）はSEA（統計エネルギー解析）、空力騒音にはFW-H法と、目的に応じたツールの棲み分けがある。この全体像をまず頭に入れておくと、後の話が理解しやすくなるよ。

出発点は波動方程式だ。音圧 $p(\mathbf{x}, t)$ に対して：

ここで $c$ は音速（空気中で約343 m/s）。これを周波数領域に変換すると、角周波数 $\omega$ で振動する定常音場に対するHelmholtz方程式が得られる：

その通り。波長 $\lambda = c/f$ に対して1波長あたり最低6要素（2次要素なら3要素）が経験則だ。例えば1000 Hzなら $\lambda \approx 0.34$ m だから、要素サイズは約57 mm以下にする必要がある。周波数が上がるとメッシュが爆発的に増えるのが騒音予測の宿命だね。

無限剛壁に埋め込まれた振動面（バッフル面）からの音響放射にはRayleigh積分が使える。面の法線速度 $v_n(\mathbf{y})$ から任意の点 $\mathbf{x}$ の音圧を直接求める：

ここで $r = |\mathbf{x} - \mathbf{y}|$ は音源点から受音点までの距離、$\rho_0$ は空気密度。$e^{-jkr}/r$ はグリーン関数で、点音源からの球面波を表す。スピーカーの放射パターンの計算なんかはこれで十分だ。

そこで登場するのがKirchhoff-Helmholtz積分方程式とBEMだ。任意の閉曲面からの放射を扱える、もっと一般的な定式化だよ。

Helmholtz方程式のグリーンの定理を適用すると、Kirchhoff-Helmholtz積分方程式が得られる。外部領域（放射問題）では：

$G(\mathbf{x},\mathbf{y}) = \dfrac{e^{-jk|\mathbf{x}-\mathbf{y}|}}{4\pi|\mathbf{x}-\mathbf{y}|}$ が3Dの自由空間グリーン関数。$c(\mathbf{x})$ は境界上で $1/2$、外部領域内で $1$ となる係数だ。$\partial p/\partial n = j\omega\rho_0 v_n$ だから、構造FEMで求めた表面速度 $v_n$ を代入すると音圧 $p$ が求まる。

そう。外部放射問題ではFEMだと音が遠方まで伝わるから、巨大な空気領域をメッシュ化したうえに無反射境界（PML等）を設定する必要がある。BEMなら表面メッシュだけで済むし、放射条件（Sommerfeld条件）が自動的に満たされる。自動車のパスバイ騒音のような外部放射問題ではBEMが圧倒的に有利だ。

主要な指標は2つ。まず音圧レベル（SPL）は受音点での音圧を基準値と比較したもの：

次に音響パワーレベル（$L_W$）は音源自体の放射エネルギーの指標で、観測位置に依存しない：

音響パワー $W$ は振動面の音響インテンシティ $I$ を面積積分して求める：

そうだ。規格でいうと、環境騒音の規制値はたいていSPLで規定されるけど、騒音源の性能比較には $L_W$ を使う。ISO 3744やISO 3745で音響パワーの測定法が定められているから、シミュレーション結果との比較はこの枠組みで行うことが多い。

基本は同じだけど、騒音予測特有のポイントがある。強制振動の周波数応答解析で、表面の法線速度 $v_n(\omega)$ を出力するのがゴールだ。FEMの運動方程式は：

変位 $\mathbf{u}$ が求まれば表面速度は $v_n = j\omega\, \mathbf{u} \cdot \mathbf{n}$ で得られる。実務で重要なのは減衰モデルの選び方だ。Rayleigh減衰（$\mathbf{C} = \alpha\mathbf{M} + \beta\mathbf{K}$）は手軽だけど、広い周波数帯域では減衰比が非物理的に変動する。自動車のNVH解析では実測モーダル減衰を周波数ごとに設定するのが一般的だよ。

その通り。減衰比が実際の2倍になるだけで、共振ピークのSPLが6 dBも変わることがある。だから構造FEMの段階で実験モーダル解析との相関（MAC値）をしっかり確認しておくことが、騒音予測の信頼性を決めるんだ。

表面 $S$ を $N$ 個の境界要素に分割して、積分方程式を離散化する。行列形式で書くと：

$\mathbf{p}$ は節点音圧ベクトル、$\mathbf{q} = \partial p/\partial n = j\omega\rho_0 v_n$ は法線方向の音圧勾配。構造FEMから $v_n$ がわかっているから $\mathbf{q}$ は既知で、この連立方程式を解けば表面音圧 $\mathbf{p}$ が求まる。表面音圧がわかれば、任意の外部点（フィールドポイント）の音圧も積分で計算できる。

いい質問だ。BEMの行列は密行列（フル行列）になる。全ての節点が全ての節点と相互作用するからね。$N$ 個の節点なら $N \times N$ のフル行列で、メモリは $O(N^2)$、求解は直接法で $O(N^3)$。これが大規模問題でのBEMのボトルネックで、後で話すFMM-BEMで解決するんだ。

片方向連成（弱連成）は最もシンプルで、実務の大半はこれだ。構造FEMで表面速度を求め、それをBEMの境界条件としてそのまま渡す。音圧が構造に及ぼす反力（音響負荷）を無視する。

双方向連成（強連成）は、音響負荷が構造の振動に影響する場合に必要だ。薄い鋼板パネルが水中で振動するケースなどが典型例。連成方程式は：

$\mathbf{L}$ が音圧→構造力のカップリング行列。実務では、空気中の騒音予測ならほとんど片方向で十分。空気の密度は構造に比べて3桁小さいから、音響負荷の影響は無視できることが多い。水中音響や車室内の閉空間音響では双方向が必要になるケースがあるよ。

FMM（Fast Multipole Method）が革命的だった。遠方の要素群からの寄与をまとめて計算することで、行列-ベクトル積を $O(N^2)$ から $O(N\log N)$ に落とせる。メモリも $O(N)$ になるから、従来BEMでは不可能だった100万節点クラスの問題が解けるようになった。

その通り。Siemens Simcenter 3DやFFT ACTRAN、ESI VA Oneといった商用ツールはFMM-BEMを標準搭載している。OpenBEMやFastBEMのようなオープンソースも選択肢にはあるけど、大規模産業問題では商用ツールの安定性と前後処理の完成度がまだ差がある。

SEA（Statistical Energy Analysis）の出番だ。構造や音場を「サブシステム」に分割し、各サブシステムのエネルギーの流れをパワーバランス方程式で記述する：

サブシステム間のパワーフローは連成損失係数 $\eta_{ij}$ を使って $P_{ij} = \omega\eta_{ij}(E_i/n_i - E_j/n_j)$ と表す。$n_i$ はモード密度、$E_i$ は全エネルギー。個々のモードを追わず、統計的な平均エネルギーで記述するから、数千Hzの高周波でも数十のサブシステムで車両全体をモデル化できる。

SEAは空間的に平均化されたエネルギーしか出せないから、「車室内のどの座席が一番うるさいか」みたいな空間分布の議論はできない。あくまで「このパネルからこの空間へどれだけエネルギーが流れるか」の傾向把握が目的だ。精度は連成損失係数の信頼性に強く依存する。実測で同定するか、解析的な理論値を使うかでかなり変わるよ。

空力騒音にはFfowcs Williams-Hawkings（FW-H）方程式を使う。非定常CFD（通常はLESかDES）で計算した流体場の圧力変動から、遠方での音圧を求める音響アナロジー法だ。

第1項がモノポール（厚み騒音）、第2項がダイポール（荷重騒音）。「ret」は遅延時間での評価。自動車のサイドミラー風切り音、航空機のランディングギア騒音、風車のブレード騒音がこの手法の典型的な適用対象だ。Ansys FluentやOpenFOAMにFW-Hソルバーが組み込まれている。

そうだ。LESで十分な時間幅の非定常データを溜めないと統計的に安定したスペクトルが得られない。自動車のサイドミラー周りの空力騒音解析で、数百コア×数日の計算はざらだ。それでもFW-H法の利点は、CFDメッシュと音響メッシュを分離できること。音の伝播を再計算せずに、受音点位置だけ変えて遠方場の音圧を即座に評価できるんだ。

FEM-BEM連成による外部放射騒音予測の標準フローはこうだ：

車体の各パネル（ルーフ、フロア、ドアなど）が受音点のSPLにどれだけ寄与しているかを分解する手法だ。「どのパネルを制振処理すれば一番効くか」を定量的に判断できるから、対策のコストパフォーマンスを最適化するのに不可欠。BEMの積分をパネルごとに分けて計算するんだよ。

最重要ルールは「6要素/波長」ルールだ。最大解析周波数 $f_{\max}$ に対して：

その通り。実務では「FEM-BEMは2000 Hz程度まで、それ以上はSEA」という使い分けが多い。中間帯域をカバーするハイブリッドFEM-SEAも最近のトレンドだ。

めちゃくちゃ重要だ。車室内の騒音予測で吸音材のインピーダンスが間違っていると、中高周波のSPLが5〜10 dBずれることもある。多孔質材料のモデルにはBiot理論やJohnson-Champoux-Allardモデルが使われるけど、必要なパラメータ（流れ抵抗率、空隙率、トーチュオシティなど）の取得が結構大変なんだ。メーカーが材料データベースを提供してくれる場合もあるけど、自社で測定するのがベストだよ。

外部BEMで使うKirchhoff-Helmholtz積分方程式は、じつは物体内部の音場の固有周波数で行列が特異になるんだ。物理的には外部音場には関係ないんだけど、数値的にはその周波数で解が汚れる。対策としてBurton-Miller定式化（通常の積分方程式とその法線微分を組み合わせる）が標準で、これを使えばほぼ完全に除去できる。最近の商用BEMソルバーはデフォルトでBurton-Millerが有効になっていることが多いよ。

騒音予測（NVH: Noise, Vibration, Harshness）には専用の音響ソルバーが必要だ。主なツールはこうなっている：

いいところに気づいたね。Simcenter 3Dはサイバネット、ActranはMSCソフトウェア（ヘキサゴン）、VA OneはESI Japan、COMSOLは直販体制で、それぞれ日本語のテクニカルサポートがある。ライセンス費用だけでなく「困ったときに電話で相談できるか」は実務で死活的に重要だよ。

まさにその「ミッドフリークエンシー問題」を解決するのがハイブリッドFEM-SEAだ。長波長の振動を決定論的にFEMで解き、短波長のモードを確率的にSEAで処理する。同じモデルの中で「低周波はFEM精度」「高周波はSEA効率」を両立させる。

基本的なアイデアは、FEMのサブシステム（確定的な振動場を持つ部分）とSEAのサブシステム（ランダムな短波長場を持つ部分）を「拡散場の相反性」で結合する。Langley & Cotoni (2004) の定式化が基盤で、Simcenter 3DとVA Oneが実装している。自動車のダッシュボード周りのように、剛性の高い骨格部材（FEM）と柔らかいパネル（SEA）が混在する領域で特に威力を発揮するよ。

サロゲートモデルは自動車OEMで実用レベルに近づいている。設計初期段階でFEM-BEMを何百ケースも回すのはコスト的に不可能だから、少数のFEM-BEM結果で学習したNNを最適化に使うのは理にかなっている。PINNはまだ学術論文が中心だけど、ポテンシャルは非常に大きい。「5年後にはCAEエンジニアの必須スキルになっている」と予測する人もいる。

最近の騒音予測のゴールは「dBを下げる」から「音質を設計する」に変わりつつある。EV（電気自動車）は内燃機関のマスキング音がないため、以前は気にならなかったモーターの高周波純音やインバーターのスイッチングノイズが際立つようになった。

そこで使われるのが心理音響指標だ：

CAEでは、BEMやSEAで計算した受音点の音圧時刻歴・スペクトルをポスト処理で心理音響指標に変換する。Simcenter Testlab、HEAD acoustics ArtemiS、ANSYSのSound Analysis Toolboxなどがこの変換機能を持っている。EVの開発では「走行中のシャープネスを1.5 acum以下」みたいな目標値が設定されることもあるんだ。

10 dBのずれは珍しくないんだけど、原因は大きく3つに分類できる：

その通り。「構造がダメなら音響は絶対に合わない」が鉄則だ。逆に、構造FEMと実測の相関がしっかり取れていれば、BEMのSPLは3〜5 dBの精度で一致するのが普通だよ。

典型的なメッシュ不足の症状だ。6要素/波長を下回ると数値分散誤差が急激に増大し、位相誤差が蓄積してSPLが大きく狂う。対策は：

外部BEMの非一意性問題（fictitious eigenfrequency problem）だね。物体内部の音響固有周波数で行列の条件数が悪化し、解が信頼できなくなる。対策は2つ：

最近のSimcenter 3DやActranはBurton-Millerがデフォルトだから、通常は問題にならない。ただし古いバージョンや自作BEMコードでは要注意。「なんか変なスパイクが出る」と思ったら、まずBurton-Millerの有効化を確認してみるといい。

SEAの精度は連成損失係数（CLF: Coupling Loss Factor）の質で決まると言っても過言じゃない。よくある問題は：

対策としては、パワー注入法（PIM: Power Injection Method）で実測CLFを同定するのが最も確実だ。各サブシステムに順番に振動を入力し、エネルギー応答行列から逆算する。手間はかかるけど、これをやるとSEA結果の信頼性が格段に上がるよ。

騒音予測は構造FEMの精度が土台だから、まずはモーダル相関をしっかりやることから始めるといい。構造が合えば音響は自然とついてくる。あと、周波数帯域で手法を切り替える判断力は実務で鍛えるしかないから、最初はFEM-BEMの低周波問題で経験を積むことを勧めるよ。分からないことがあったらいつでも聞いてくれ。