音響解析 — CAE用語解説
音響解析
音響解析って「音の解析」ですよね。車の車内の静粛性とかスピーカーの指向性とか？ FEMで音が計算できるって驚きです。
理論と物理
音響波動の基本方程式
音響解析の最も基本的な支配方程式は何ですか？流体の運動方程式からどう導かれるんですか？
線形音響の基礎は、連続の式、運動量保存則、状態方程式から導かれる波動方程式です。非粘性・等エントロピー流れを仮定し、圧力や密度を平均値と微小擾動に分けて線形化すると、次のヘルムホルツ方程式が得られます。
「線形化」というのは具体的にどういう操作をしているんですか？なぜ非線形項を無視できるのですか？
良い着眼点だ。例えば、圧力を
音速
音速は媒質の体積弾性率と密度で決まる。気体では
数値解法と実装
FEMによる離散化と境界条件
ヘルムホルツ方程式を有限要素法（FEM）で解く場合、どのように離散化するんですか？特に波数kが含まれる項はどう扱う？
ガラーキン法を用いて弱形式を導く。支配方程式に重み関数
境界積分の
それが境界条件の本質だ。硬い壁（完全反射）では、法線方向粒子速度が0、つまり
メッシュサイズはどのように決めればいいですか？周波数によって変える必要がありますか？
音響FEMでは、波長を十分に分解できるメッシュが必要だ。経験則として、解析対象の最高周波数
実践ガイド
室内音響解析のワークフロー
会議室のような室内の音響特性（残響時間など）をCAEで評価する場合、どのような手順で進めればいいですか？
典型的なワークフローはこうだ。1) 3D CADで室内形状を作成。2) 壁・天井・床の材料ごとに音響インピーダンス（吸音率）を設定。JIS A 1409やISO 354規格に基づく実測値を使う。3) メッシュ生成。先ほどのルールに従い、例えば500 Hzまで解析なら要素サイズ~70mm以下。4) 音源を定義。点音源や面振動源が使える。5) 周波数応答解析を実行。125, 250, 500, 1000, 2000, 4000 Hzの1/1オクターブバンド中心周波数が一般的。6) 後処理で音圧分布、特定点の周波数応答、そして残響時間（Sabineの式やエネルギー減衰曲線から算出）を評価する。
吸音率の設定で注意すべき点は？「吸音率0.5」をそのまま入力していいんですか？
絶対にダメだ。吸音率は実数だが、音響インピーダンスは一般に複素数だ。吸音率αからインピーダンスZを逆算する近似式はあるが、特に低周波数では誤差が大きい。実務では、材料メーカーが提供する「複素特性インピーダンス」または「複素反射係数」の周波数データを直接入力すべきだ。例えば、音響材料メーカーのデータシートには、125-4000Hzの範囲で実部と虚部が掲載されている。COMSOLの「インピーダンス境界条件」では、この実数部と虚数部をテーブル入力できる。
解析結果の残響時間が実測と合わない場合、まずどこを疑えばいいですか？
以下のチェックリストを順に確認する。1) **材料データ**: インピーダンス値が正しい周波数で入力されているか。特に低周波数（125Hz）のデータは不確実性が高い。2) **メッシュ**: 最高周波数に対して十分に細かいか。粗いメッシュは減衰を過大評価する。3) **漏れ**: ドアの隙間など、音エネルギーが抜ける経路をモデル化しているか。4) **家具の影響**: 会議室の机や椅子は音を散乱・吸収する。これらを等価吸音面積として追加する必要があるかもしれない。5) **ソルバー設定**: 減衰（もしあれば）の設定や、求解周波数点数が十分か。
ソフトウェア比較
主要ソフトウェアの特徴と適用範囲
Ansys、COMSOL、Abaqusで音響解析を行う場合、それぞれどういう特徴や得意分野があるんですか？
大きく分かれる。**Ansys**は「Mechanical + ACT Acoustics」または「Harmonic Acoustics」が主流。構造振動との連成（FSI）が強く、自動車の車内騒音や家電の振動騒音に強い。大規模モデルの並列計算性能に優れる。**COMSOL Multiphysics**は「Pressure Acoustics」モジュールが核で、その名の通り多物理場連成が本領。音-構造-熱-流れの連成を一つの環境で柔軟に設定できる。研究開発向け。**Abaqus**は「Abaqus/Standard」のAcoustic要素を用いる。構造解析との統合がシームレスで、タイヤの路面騒音や水中音響など、固体と流体の大変形を伴う複雑な接触問題にも対応できる。
無料や低価格のオープンソースソフトウェア（OSS）ではどうですか？実用レベルで使えますか？
限定的だが可能だ。**CalculiX**や**Code_Aster**といったFEMソルバーには音響要素がある。前処理・後処理には**Salome-Meca**や**Gmsh**、**ParaView**を組み合わせる。しかし、商用ソフトのような高度な吸収境界条件（PML）や、複雑な周波数依存インピーダンスの設定、効率的な周波数掃引ソルバーは自分で実装するか、既存のマクロを探す必要がある。教育や小規模な等方性モデルの検証には有用。実務で製品開発に使うには、サポートや検証済み機能の面で商用ソフトに分がある。
「BEM（境界要素法）を使った音響解析」というのも聞きます。FEMと何が違うんですか？
根本的な違いは、FEMが領域全体をメッシュ分割するのに対し、BEMは境界面のみをメッシュ化する点だ。これにより、無限遠の音場（放射音）を自然に表現できるのが最大の利点。自動車の外側放射騒音や、スピーカーの遠方音場解析に適する。代表的な商用BEMソルバーは、**Altair Acoustics**（旧LMS Virtual.Lab）、**FFT Actran**（FEM/BEM両方可能）、**COMET**など。欠点は、方程式の係数行列が密行列になるため、大規模モデルでは計算コストが高くなること。近年は高速BEM（Fast Multipole BEMなど）で克服されつつある。
トラブルシューティング
よくある収束エラーと対策
周波数応答解析で、特定の周波数（例えば275Hz）で解が発散したり、異常に大きな音圧値が出たりします。原因は何ですか？
それはほぼ間違いなく「固有モードの共振」だ。音響キャビティには、ドラムの膜のような固有振動数（音響モード）が存在する。減衰が十分にモデル化されていないと、その周波数で応答が理論上無限大に発散する。まず、その周波数付近の固有値解析を実行し、固有周波数と一致するか確認せよ。対策は3つ。1) 現実的な減衰（材料内部損失や境界での吸収）をモデルに追加する。2) 構造と連成解析を行う（現実の壁は完全剛体ではない）。3) 結果を評価する際、その共振周波数を避けるか、平均化して評価する。Ansysでは「Modal Damping」や「Rayleigh Damping」の設定が有効だ。
メッシュを細かくしたら、むしろ結果がおかしくなり、ソルバーが「Matrix is singular」というエラーを出しました。なぜですか？
それは「数値的ロック」または「スピューリアスモード」の典型的な症状だ。音響FEMでは、圧力のみを変数とする定式化（圧力公式）を用いることが多い。これは非圧縮性流れのFEMで起こる「体積ロック」に似て、メッシュを極端に細かくすると、要素レベルで
PML（完全整合層）を設定したのに、境界で明らかに反射が起きているように見えます。考えられる原因は？
PMLの設定ミスだ。以下の点をチェックせよ。1) **層の厚さ**: 一般的に、対象周波数の最低波長の1/2〜1波長分の厚さが必要。500Hzで空気中なら、約34cm以上。2) **減衰係数の分布**: 内部領域からPML境界に向かって滑らかに増加する関数（二次関数など）が使われる。急峻すぎると反射する。3) **PMLの配置**: 波の進行方向に対して垂直に配置する必要がある。複雑な形状や散乱体では、PMLを球状や円筒状に配置する。4) **メッシュ**: PML領域内でも、少なくとも3〜5層の要素が必要。COMSOLの「スケーリング座標系」やAnsysの「PML要素タイプ」の設定をデフォルトから変えていないか確認すること。
構造-音響連成解析で、流体（空気）のメッシュと構造メッシュのインターフェースはどう扱うべきですか？一致させる必要がありますか？
必ずしも一致させる必要はないが、整合性のある方法で結合しなければならない。主流は2つの手法。1) **一致メッシュ**: インターフェースで節点を共有する。最も単純で精度が高いが、メッシュ作成の柔軟性に欠ける。2) **非一致メッシュ**: それぞれ独立したメッシュを用い、拘束方程式（MPC）やラグランジュ乗数法で「圧力と法線方向変位の連続条件」を課す。Ansysの「FSIインターフェース」やAbaqusの「TIE」拘束がこれに該当。非一致メッシュを使う場合、音響側のメッシュが構造側より粗すぎると、高周波数のエネルギー伝達を正しく捉えられないので注意。インターフェース近くの音響メッシュは、構造メッシュのサイズと同程度か、それより細かくすることを推奨する。
関連トピック
なった
詳しく
報告