音響アナロジー — CAE用語解説

いい質問だ。音響アナロジーは、音響（音波）の伝播を記述する支配方程式が、他の物理現象の方程式と数学的に同じ形をしていることを利用する手法だ。つまり、音響解析用のソフトウェアや知見を、全く別の分野の問題に「流用」できる。例えば、ヘルムホルツ方程式

本質は「波動現象」にある。薄板の横振動を支配するのは板の曲げ波動方程式だ。例えば、四辺単純支持の矩形板の固有振動数は

対応表を覚えると便利だ。音響系では、音圧p [Pa]、体積速度U [m³/s]、音響インピーダンスZ_a [Pa·s/m³]が基本変数だ。これが機械系では、力F [N]、速度v [m/s]、機械インピーダンスZ_m [N·s/m]に対応する。電気系では、電圧V [V]、電流I [A]、電気インピーダンスZ_e [Ω]だ。このアナロジーを使えば、スピーカーの電気-機械-音響結合系を一つの等価回路でモデル化できる。

その通り、支配方程式を離散化すると同じ形式の行列方程式になる。音響のヘルムホルツ方程式をガラーキン法で弱形式化すると、次の線形システムが得られる。

これもアナロジーが成り立つ。音響の「剛体壁」（速度=0）は、構造の「固定端」（変位=0）のディリクレ条件に対応する。音響の「開放端」（音圧≈0）は、構造の「自由端」（力=0）のニューマン条件に近い。さらに、吸音材を表す「インピーダンス境界条件」

それが大きな違いだ。構造解析では形状の曲率に応じてメッシュを細かくするが、音響では解析最高周波数における波長が基準になる。経験則として、1波長あたり6～10個以上の要素が必要だ。例えば、空気中で最高周波数1000Hzを解析する場合、波長は約0.34mだから、要素サイズは最大で0.034～0.056mにしなければならない。一方、構造の固有振動解析では、モード形状を捉えられれば良いので、これほど細かいメッシュは必要ない場合が多い。

まず、支配的な伝播経路を特定せよ。1) 構造振動が直接放射する「構造放射音」か、2) 内部流体の振動が壁を伝って放射される「伝達音」か。例えば、電気洗濯機のモータ騒音は1が主で、エアコンのダクト音は2が主だ。これで、連成界面（構造と音響領域が接する面）をどこに設定するかが決まる。最初のチェックリストは、「音源の特定」、「連成界面の定義」、「解析対象周波数範囲の設定（例: 20-2000Hz）」の3点だ。

音響領域の材料には、密度と音速の2つが必須だ。常温の空気なら、密度ρ0=1.225 kg/m³、音速c0=340 m/sだ。水ならρ0=1000 kg/m³、c0=1500 m/sと大きく変わる。構造側は、ヤング率、ポアソン比、密度、損失係数（減衰）が必要だ。実務では、カタログ値やJIS規格（例えばJIS K 6920-2のプラスチック物性）を参照する。忘れがちなのが「減衰」で、鋼は損失係数η=0.001程度だが、ゴムはη=0.1以上と大きく、結果に与える影響が大きい。

良い着眼点だ。これが実務上のハードルの一つだ。構造メッシュが粗く、音響メッシュが細かい場合、界面で変位や力を補間（マッピング）する必要がある。多くのCAEソフト（Abaqus, Ansys）は自動でこの処理を行うが、マッピング誤差が生じる。対策は、可能な限り界面でメッシュサイズと節点位置を一致させることだ。それが無理なら、構造側の界面メッシュを、音響側の要求する要素サイズ（先述の波長基準）に合わせて局所細分化する。

コアとなる有限要素法の理論は同じだが、ソルバーの実装とユーザインタフェースに違いがある。Ansysは、Acoustics専用の「ACT」アドオンがあり、Harmonic AcousticsやModal Acousticsなど、解析タイプが細分化されている。一方、Abaqusでは、「Acoustic, Steady-State」や「Acoustic, Transient」といったステップ定義で音響解析を設定する。音響要素は「AC3D8」などと明示的に選択する。最大の違いはCOMSOL Multiphysicsで、「Pressure Acoustics」物理場を選択し、他の物理場（固体力学、流体流れ）とGUI上で直感的に連成設定ができる点だ。

これが各社の技術の見せ所だ。3者とも「無限要素」または「完全整合層（PML）」を実装している。Ansysでは「Infinite Acoustic」要素をメッシュの外縁に追加する。Abaqusでは「Acoustic Infinite」要素を使う。COMSOLは「Perfectly Matched Layer」ドメインを設定するのが一般的だ。PMLは吸収効率が高く、特に斜め入射波に対して優れるが、メッシュが複雑化する。商用ソフトでは、Ansys Mechanicalの「Acoustic PML」機能が2019R3以降で強化され、設定が簡便になった。

歴史的に、専用ソフトは中高周波数帯（統計的エネルギー解析: SEA）や、非線形音響、ファン騒音などの高度なモデリングに強かった。しかし、汎用CAEも追いついている。例えば、Ansys 2023R1では「Harmonic BEM」ソルバーが導入され、大規模な外部放射問題の計算効率が向上した。一方、Actranは自動車業界でタイヤ騒音やマフラー音の解析で強く、豊富な実験検証済みのテンプレートを持つ。選択は「解析対象の周波数帯」と「必要な物理モデルの専門性」による。1kHz以下の低中音域で構造連成が主なら、汎用CAEで十分なことが多い。

それはほぼ間違いなく「非物理的共振」、いわゆる「フィクティシャスモード」または「スプリアスモード」だ。音響領域のメッシュが粗すぎて、波長を十分に表現できていない場合に発生する。先ほど言った「1波長に6要素以上」はこれを防ぐための経験則だ。特に、低音圧領域（節）でメッシュが粗いと、音圧が計算上ゼロに固定され、剛性が過大評価されて固有周波数が高く出てしまう。対策は、メッシュを細かくするか、要素次数を上げる（1次要素から2次要素へ）ことだ。

音圧を変数とするFEMでは、全体の音圧レベルが一意に定まらない「剛体モード」が存在する。つまり、領域全体の音圧が一様に上下するモード（ゼロ周波数モード）に対応する。これが質量行列を特異にする。構造解析で言えば、物体が浮いている状態だ。これを防ぐには、少なくとも1つの節点で音圧を固定（ディリクレ条件）する必要がある。実務では、開放端境界（p=0）を1箇所設けるか、あるいは「圧力基準点」を設定する機能をソフトウェアで使う。COMSOLでは「Point Constraint」、Ansysでは「Acoustic Pressure」の固定支持がそれに当たる。

いくつかの戦略がある。1) 部分構造合成法（CMS）を使う。構造側を超要素（減次モデル）化する。Abaqusの「Substructure」、Ansysの「CMS」機能を使えば良い。2) 対称性・周期性を利用する。ファンやフィンのように周期的な構造なら、1ピッチ分だけモデル化し、周期性境界条件を適用する。3) ソルバー設定を見直す。連成問題は非対称行列になることが多く、直接法ソルバー（スパースソルバー）が向く。Ansysでは「Sparse Direct」、Abaqusでは「Direct」を選択する。反復法ソルバーは収束しないか、時間がかかる場合がある。