波動・音響の基礎 — 波動方程式・音響圧・NVH・構造音響連成まで
波動とは何か — CAEにおける波動現象
波動とは「エネルギーと情報が媒質を通じて伝播する現象」です。CAEの世界では、音響解析（NVH）・超音波非破壊検査・衝撃波・電磁波伝播と、非常に幅広い物理現象を支配する共通の数学フレームワークです。
NVH解析って「Noise, Vibration, Harshness」ですよね。振動解析とは別物ですか？ どう違うんでしょう。
ざっくり言うと、振動解析は構造（固体）の変形を見て、NVHはそこから生まれる音と乗り心地の問題まで包括的に扱う。自動車のNVHで言えば、エンジンの燃焼力がボディ構造を振動させ（Vibration）、その振動が車室内の空気を揺らして音（Noise）になり、路面凸凹が急に来る衝撃感（Harshness）まで含む。CAEで言えば、振動は構造FEM、Noiseは音響FEM/BEMや統計的エネルギー解析（SEA）が担当する。つまり波動・音響の理論は特に「振動が音になる」メカニズムの理解に不可欠なんだ。
1.1 波動の分類
| 種類 | 媒質の変位方向 | 例 |
|---|---|---|
| 縦波（Longitudinal） | 伝播方向と同じ（疎密波） | 空気中の音、固体中の圧縮波（P波） |
| 横波（Transverse） | 伝播方向と垂直 | 電磁波、固体中のせん断波（S波） |
| 表面波（Surface wave） | 表面に沿って減衰 | 地震のRayleigh波、超音波探傷 |
| 曲げ波（Flexural） | 板・梁の面外変位 | 車体パネルの振動、騒音放射の主役 |
1.2 波数・波長・周波数・位相速度の関係
$c$ は位相速度、$f$ は周波数（Hz）、$\lambda$ は波長（m）、$\omega = 2\pi f$ は角周波数（rad/s）、$k = 2\pi/\lambda$ は波数（rad/m）です。音速 $c_0 \approx 340\,\text{m/s}$（20°C空気）で 1000 Hz の音波の波長は $\lambda = 0.34\,\text{m}$。これは自動車のボディパネルと同じオーダーで、構造と音響が強く連成する理由の一つです。
1.3 分散性と非分散性
空気中の音波は非分散（位相速度が周波数によらず一定）。一方、板の曲げ波は分散性があり、高周波ほど速く伝播します（$c_b \propto \sqrt{f}$）。この分散性がコインシデンス効果（後述）の原因です。
波動方程式
2.1 1次元波動方程式の導出
弦の微小要素の運動方程式から（張力 $T$、線密度 $\rho_L$）：
音響圧 $p'$ も同形の方程式に従います（連続式と運動量保存から導出）：
2.2 D'Alembert解
1次元波動方程式の一般解は「右向き進行波」と「左向き進行波」の重ね合わせです：
調和波（正弦波）の場合：$u(x,t) = A e^{i(kx - \omega t)} + B e^{i(-kx - \omega t)}$。複素指数表示が数値計算では便利で、最終的に実部を取ります。
2.3 双曲型方程式と熱方程式の本質的違い
波動方程式と熱伝導方程式って見た目が似てますよね。どう違うんですか？
見た目は似てるけど物理的に全然違う。波動方程式は $\partial^2 u/\partial t^2 = c^2 \nabla^2 u$ で2階の時間微分がある（双曲型）。熱方程式は $\partial T/\partial t = \alpha \nabla^2 T$ で1階の時間微分（放物型）。双曲型は情報が有限の速度 $c$ で伝播する――音を叩いてから相手に届くまでに時間差がある。放物型は情報が瞬時に無限速で伝わる――熱伝導の厳密解では系のどこかを加熱した瞬間に無限遠でも温度が変化する（ただし微小な変化）。数値解法も違って、双曲型はCFL条件が出てくるし、放物型は陰解法で安定的に解ける。
2.4 3次元波動方程式（Helmholtz方程式）
時間高調波（角周波数 $\omega$ の定常振動）を仮定すると $p'(x,t) = \hat{p}(x) e^{-i\omega t}$ と分離でき、Helmholtz方程式が得られます：
これが音響有限要素法（Acoustic FEM）や境界要素法（BEM）の支配方程式です。
音響の基礎
3.1 音速と媒質
$B$ は体積弾性率です。空気中の音速は温度で変わり、$c_0 \approx 331.3 + 0.606\,T_C\,\text{[m/s]}$（$T_C$ は摂氏）。20°Cで約 340 m/s、1000°Cの燃焼ガス中では約 600 m/s になります。
| 媒質 | 音速 (m/s) | 特記事項 |
|---|---|---|
| 空気（20°C） | 343 | NVH解析の基本 |
| 水（20°C） | 1481 | 超音波検査・水中音響 |
| 鋼（縦波） | 5900 | 超音波探傷・衝撃波 |
| コンクリート | 3000〜4000 | 建築音響・構造健全性評価 |
| CFRP（面内） | 8000超 | 積層方向依存性大 |
3.2 音響インピーダンスと反射
特性音響インピーダンス（固有音響抵抗）：
空気：$Z_0 \approx 410\,\text{Rayl}$、水：$Z_0 \approx 1.5 \times 10^6\,\text{Rayl}$。インピーダンス差が大きいほど反射率が高い。これが空気–水界面で音の99.9%が反射する理由であり、超音波探傷でカプラント（グリス等）を使う理由でもあります。
3.3 音圧レベルと音響パワーレベル
基準音圧 $p_\text{ref} = 20\,\mu\text{Pa}$（人間の可聴限界に相当）。音響パワーレベル：
自動車の設計仕様書に「エンジン音 65dB以下」みたいに書いてあるんですが、CAEでこの dB を計算するにはどうするんですか？
まず Ansys Mechanical や Abaqus の音響解析で車室内の音圧分布（Pa）を周波数ごとに計算して、運転席位置の $p_\text{rms}$ を取り出す。それを上の式で dB に変換する。ただし「人間がどう感じるか」はさらに A特性フィルタ（dBA）をかけて評価するのが普通だ。低周波のドスンという音と高周波の甲高い音は同じ dB でも聴こえ方が違うから、聴感補正した dBA で比較する。CFD で計算した流れ場のノイズ源から Ffowcs Williams-Hawkings（FWH）方程式で遠距離の音圧を計算するアプローチもある。排気音やエアコンファン騒音の解析ではよく使う。
定在波と共鳴
4.1 1次元定在波
固定端（$x = 0$）と自由端（$x = L$）を持つ管の音響固有モード：
両端固定（または両端開口）の管：
例えば自動車の排気管（長さ $L = 1.5\,\text{m}$、$c_0 = 340\,\text{m/s}$）の1次固有周波数は $f_1 = 340/(4 \times 1.5) \approx 57\,\text{Hz}$。エンジンの燃焼周波数がこれに一致するとマフラーが共鳴して排気音が増大します。
4.2 3次元矩形室の固有周波数
$m, n, p$ はそれぞれ $x, y, z$ 方向のモード次数。室内音響（コンサートホール・録音スタジオ・自動車車室）の設計では低周波域のモード分布がこれで決まります。典型的な乗用車の車室（$L_x \approx 1.8$, $L_y \approx 1.2$, $L_z \approx 1.0\,\text{m}$）の最低次固有周波数は 90 Hz 程度。
4.3 コインシデンス効果と音響透過損失
板の曲げ波速度 $c_b = (2\pi f)^{1/2} (B/\rho_s)^{1/4}$（$B$ は曲げ剛性）が音速 $c_0$ と一致する周波数をコインシデンス周波数と言います：
コインシデンスって実務で何が問題になるんですか？
この周波数で遮音性能（透過損失TL）が谷を形成して急激に悪化する。例えば 2 mm 厚の鋼板パネルのコインシデンス周波数は約 6 kHz。これが自動車の高速走行中の風切り音が車内に入りやすい周波数帯と重なると問題になる。対策は板厚の変更（$f_c$ をシフト）、非対称積層材の使用、制振シートによる損失ファクター増加など。Ansys の Acoustic 解析でパネルの TL を周波数ごとに計算して、コインシデンスによる谷がどこに出るかを設計段階で確認できる。
散乱・回折・反射
5.1 Snellの法則（音響版）
媒質1（音速 $c_1$）から媒質2（音速 $c_2$）への音波入射：
$c_2 > c_1$ のとき全反射が起きる（全反射臨界角 $\theta_c = \arcsin(c_1/c_2)$）。水中から鋼板への音波入射では全反射が起きやすく、超音波探傷の設計に影響します。
5.2 遮音壁の挿入損失
音源と受音点の間に障壁を置いたときの挿入損失（Maekawa の近似式）：
$\delta$ はフレネル数の分子（音源・障壁頂点・受音点の経路差）、$N$ はフレネル数。道路沿いの防音壁設計や工場内騒音対策でよく使われる近似式です。ただし高精度が必要な場合は音響BEMや回折積分を使います。
5.3 室内音響と残響時間
Sabineの残響時間：
$V$ は室体積（m³）、$\alpha_i$ は各面の吸音率、$S_i$ は面積（m²）。コンサートホールでは $T_{60} = 1.8 \sim 2.2\,\text{s}$、会話が主の会議室では $0.4 \sim 0.6\,\text{s}$ が目安です。
構造振動と音響放射（構造音響連成）
6.1 振動する構造から音が放射するメカニズム
表面が速度 $v_n$ で振動する構造は周囲の流体を押して音響圧 $p$ を生成します。これが構造音響連成（Vibroacoustics）です。連成の強さは音響インピーダンス比で決まります：
この比が小さい（重い構造・軽い流体）ほど連成が弱く（空気と鋼板）、大きい（軽い構造・重い流体）ほど強い（水中の薄板）。自動車のNVHでは連成が弱いため、構造解析→音響解析の一方向連成（弱連成）で十分なことが多いです。
6.2 放射効率
放射効率 $\sigma$ は「実際の放射音響パワー」と「全面が同位相で振動する剛体ピストンが放射するパワー」の比です：
低次の振動モードは放射効率が低く（体積速度が相殺する）、モード次数が上がるとコインシデンス周波数付近で $\sigma \to 1$ に近づきます。
6.3 SEA（統計的エネルギー解析）の概念
高周波域（波長がシステムの特性寸法より十分小さい）では、個々のモードではなく「サブシステム間のエネルギーフロー」を統計的に扱うSEAが有効です：
$P_{12}$ はサブシステム1から2へのエネルギーフロー、$\eta_{12}$ は連成損失因子、$E_1, E_2$ は各サブシステムのエネルギーです。車体の1000 Hz 以上の高周波騒音解析、船舶の機械室騒音、飛行機客室の騒音解析でSEAが標準的に使われます。
じゃあ FEM と SEA、どっちを使えばいいかはどうやって決めるんですか？
目安として「波長がシステムの代表寸法の1/6より大きい」低周波域はFEM、逆に小さい高周波域はSEA。境目は解析対象によるけど、自動車の場合だと大体 500〜1000 Hz あたりが境界線になることが多い。その境界域（mid-frequency）が一番難しくて、FEMは計算コストが爆発し、SEAは統計仮定が怪しくなる。そこを狙い打ちしたハイブリッド FEM/SEA 手法が VA One（ESI Group）や Actran などのソフトに実装されている。実務ではドライブサイクルの主要周波数帯（どの周波数が問題か）をまず特定して、それに合わせて手法を選ぶ。
数値音響解析
7.1 有限要素法（FEM）による音響解析
Helmholtz 方程式の弱形式から音響FEMの離散方程式が得られます：
$\mathbf{M}_a$：音響質量行列（体積積分）、$\mathbf{C}_a$：吸音境界からの減衰行列、$\mathbf{K}_a$：音響剛性行列。構造との連成は：
$\mathbf{L}$ が構造-音響連成マトリクス（界面での法線方向連成）です。Ansys Mechanical では FLUID30/FLUID220/FLUID221 要素、Abaqus では AC2D4/AC3D4 要素が音響要素として使えます。
7.2 境界要素法（BEM）の利点
BEMは境界（表面）のみをメッシュ化するため、無限音響領域（外部音場）の自動処理が得意です。Sommerfeld放射条件を自動的に満足するため、マフラー外の放射音・屋外の騒音分布計算に向いています。ただし非対称・稠密な行列になり、計算コストは $O(N^2) \sim O(N^3)$（高速多重極法で改善可能）。
7.3 非反射境界条件（PML・Absorbing BC）
FEMで無限領域を扱うには、計算領域の外縁で人工的に波を吸収する完全整合層（PML: Perfectly Matched Layer）が使われます。PML 内で：
$\sigma_i$ は境界での吸収係数。Ansys Acoustics・Comsol・OpenFOAM の音響モジュールで利用可能です。
実務でのNVH解析のポイント
8.1 周波数帯域と解析手法の選択
| 周波数帯域 | 推奨手法 | 主な対象 |
|---|---|---|
| 〜200 Hz | 構造FEM + 音響FEM | アイドル振動・ブーミング音・低周波こもり音 |
| 200〜1000 Hz | FEM + SEA ハイブリッド | ロードノイズ・ギアノイズ・排気音 |
| 1000 Hz〜 | SEA / 統計手法 | 風切り音・吸音材評価・高周波透過損失 |
8.2 減衰材・吸音材のモデリング
粘弾性制振シート（CLD: Constrained Layer Damping）の損失係数：
実務ではObata式などの解析解か、FEMで複素剛性（$E^* = E(1 + i\eta)$）を使って損失を評価します。
8.3 伝達経路解析（TPA）とCAE連携
TPA（Transfer Path Analysis）は「音源から受音点へのどの経路がどれだけ寄与しているか」を定量化する手法です。CAEでは周波数応答関数（FRF）マトリクスを計算し、実測の力入力との積で合成音圧を予測します：
$H_{pk}$ は経路 $k$ から評価点 $p$ への伝達関数、$F_k$ は経路 $k$ への入力力。CAEで FRF を計算して実測の力入力から車内音圧を合成すると、「どの締結部を強化すれば最も効果的か」が設計段階で分かります。
TPAって実際の自動車開発でも使うんですか？
かなり標準的なツールになっている。量産車開発では試作車が出来たらまずTPAで「どのマウントからどれくらいの音が来ているか」を実測で洗い出す。最近はソースTPA（OTPA）という実走行データだけから解析する手法も普及してきた。CAEとの連携では、試作前に車体FEMで各経路の伝達関数を計算しておき、サプライヤーから貰ったパワートレインの振動データを入力として車内音圧を予測する。目標値を超える経路を見つけたら構造補強か防振マウントの変更を設計にフィードバックする、という流れが一般的だ。
関連インタラクティブツール
理論を手を動かして確認しよう
- 1自由度系周波数応答ツール — 固有振動数・減衰比から応答倍率カーブをリアルタイム生成
- オイラー座屈ツール — 構造の不安定現象を可視化
- 梁のたわみ計算ツール — 振動解析の基礎となる梁理論を確認
なった
詳しく
報告