波の基本性質 — 波長・振動数・位相速度から音響CAEまで
関連トピック
1. 波とは — エネルギーの伝搬
超音波探傷検査（UT）って非破壊試験で使いますよね？ 波のどんな性質を利用してるんですか？
反射と音速の差を利用している。超音波を材料内部に送信して、き裂や空洞など「音響インピーダンスが異なる境界」で反射して戻ってくる時間を計測する。伝播速度（鋼で約5900 m/s）が分かれば反射源までの距離が計算できる。MRIや魚群探知機も同じ原理だよ。
波は「物質自体が移動するのではなく、エネルギーと位相（振動パターン）が伝わる現象」です。水面の波紋を見ると水の粒子は上下に動くだけで横には移動しませんが、波のパターン（エネルギー）は遠くに広がります。
2. 波の分類と基本用語
横波と縦波
- 横波（Transverse Wave）：振動方向が波の進行方向に垂直。光波・電磁波・地震のS波・水面波。固体中のせん断波もこれ。
- 縦波（Longitudinal Wave）：振動方向が波の進行方向に平行（疎密波）。音波・地震のP波・超音波の圧縮波。
波の基本用語
- 振幅 $A$：最大変位の大きさ
- 波長 $\lambda$ (m)：同位相の隣り合う点の間の距離
- 振動数（周波数）$f$ (Hz)：1秒間に通過する波の数
- 角振動数 $\omega = 2\pi f$ (rad/s)
- 波数 $k = 2\pi/\lambda$ (rad/m)：空間的な角振動数
- 位相速度 $v = f\lambda = \omega/k$ (m/s)
| 媒体・波の種類 | 波速の目安 |
|---|---|
| 音波（空気中、20℃） | 343 m/s |
| 音波（水中） | 1480 m/s |
| 縦波（鋼中） | 5920 m/s |
| 地震P波（地殻） | 5000〜7000 m/s |
| 地震S波（地殻） | 3000〜4000 m/s |
| 光（真空中） | 3×10⁸ m/s |
3. 波の基本式（正弦波）
$+x$ 方向に伝播する正弦波：
位相の意味：$kx - \omega t = \text{const}$ を満たす点が同位相面（波面）となり、時間とともに $+x$ 方向に速度 $v = \omega/k$ で移動します。これが位相速度の定義です。
1次元波動方程式（Helmholtz方程式の時間域版）：
ここで $T$ は張力、$\mu$ は単位長さあたりの質量です。この偏微分方程式を有限要素で離散化したものがFEM音響解析の基礎です。
4. 重ね合わせの原理と干渉
線形波動方程式では、複数の波の解を足し合わせた関数もまた解になります（重ね合わせの原理）：
建設的干渉（同位相）
2つの波が同位相（位相差 $\delta = 0$）のとき：
振幅が2倍になります。
相殺的干渉（逆位相）と消音の原理
逆位相（位相差 $\delta = \pi$）のとき：
完全消音。これがアクティブノイズキャンセリング（ANC）の原理です。マイクで取り込んだ雑音の逆位相波をスピーカーで出力することで騒音を打ち消します。自動車の排気音キャンセル（アクティブマフラー）にも応用されています。
唸り（ビート）
わずかに異なる周波数 $f_1$, $f_2$ の2つの波の合成：
唸りの周波数：$f_{beat} = |f_1 - f_2|$。楽器のチューニングで2つの音が一致するとビートが消えることを利用します。機械振動でも、わずかに異なる回転数の2つのモーターが干渉して唸りが生じる現象があります。
5. 定在波（Standing Wave）
同じ振幅・周波数を持つ2つの進行波が逆方向に進むとき定在波が生じます：
振幅ゼロの点（節、Node）と振幅最大の点（腹、Antinode）が固定されます。
両端固定弦の共鳴周波数
$n=1$：基本音（基本波）、$n=2,3,\ldots$：倍音（高調波）。弦楽器の音階はこの共鳴周波数の組み合わせです。
室内音響と定在波
矩形室内の定在波（Room Mode）は3次元に広がります：
室内音響設計ではこのRoom Modeが低周波での不均一な音圧分布（ブーミーな場所・聞こえにくい場所）の原因となります。FEM音響解析でMode形状を計算して吸音材の配置を最適化します。
6. 反射・屈折・回折
スネルの法則（屈折）
光の屈折と同じ関係が音波・弾性波にも成立します。
全反射
$v_2 > v_1$ のとき、入射角が臨界角 $\theta_c = \arcsin(v_1/v_2)$ を超えると全反射が起きます。超音波探傷での斜角探傷法（斜めに超音波を入射して横波モードに変換）はこの原理を応用しています。
回折と波長
波の回折（障害物の後ろに回り込む現象）は、障害物のサイズが波長と同程度以下のときに顕著です。音波（波長数cm〜数m）は壁の角で回り込み、光（波長数百nm）は同じ障害物を通りません。建物の防音壁（遮音壁）設計では回折効果が重要で、FEM/BEM音響解析で壁の高さや形状を最適化します。
7. ドップラー効果
音源または観測者が移動するとき、観測される振動数が変化します：
符号：上符号は観測者が音源に近づく場合（+）または音源が観測者に近づく場合（−）。
工学的応用
- 超音波流速計（ドップラー流速計）：流体中の粒子や赤血球からの反射波の周波数シフトから流速を計測。血流計測・パイプ流量計に使用。
- 気象レーダー：雨粒からの電波の周波数シフトで風速を計測。
- 自動車の速度測定レーダー（オービス）：マイクロ波のドップラーシフトで車速を計測。
8. 波の分散（Dispersion）
位相速度が周波数に依存する現象を分散（Dispersion）と言います。
梁の曲げ波の分散関係
弾性梁中の曲げ波の位相速度は周波数の平方根に比例します：
高周波数ほど速く伝播します。これは棒波（縦波）が分散なし（$v = \sqrt{E/\rho}$）であるのと対照的です。
構造音響解析では、薄板や梁の曲げ波の分散特性が放射音に大きく影響します。板の臨界周波数（coincidence frequency）以上では放射効率が急激に上がるため、FEM-BEM連成解析での確認が重要です。
9. FEM/BEMによる波動・音響解析
音場の支配方程式はHelmholtz方程式（周波数領域の波動方程式）：
ここで $p$ は音圧振幅、$k$ は音波数（Acoustic Wave Number）です。
FEM音響解析（有限要素法）
閉じた内部音場（室内・車室・マフラー内部）に適しています。要素サイズは最短波長の6分の1程度が目安：$\Delta x \leq \lambda_{min}/6 = c/(6f_{max})$。
車の室内騒音をFEMで解析するとき、どんなソフトを使うんですか？ 周波数によって違うんですか？
使い分けがあるんだ。1000Hz以下はFEM（Actran, NASTRAN, Abaqus Acoustics）が主流で、1000〜3000Hzはハイブリッド手法、3000Hz以上はSEA（統計的エネルギー解析、AutoSEAなど）か光線追跡法を使うのが実務の標準だ。低周波では波の干渉パターンが重要で、高周波では統計的に扱ったほうが効率的という理由だよ。BEM（境界要素法）は外部音場（放射音）に使う。
BEM音響解析（境界要素法）
外部音場（エンジン音の放射・建物外部の騒音など）に適しています。FEMと違い、内部を要素分割する必要がなく（境界のみ）、無限遠の放射条件を自動的に満たせる特長があります。ただし周波数ごとに大規模な密行列計算が必要です。
FEM-BEM連成解析（構造-音響連成）
振動する構造物（エンジン・電動モーター・パネル）が発する放射音を予測する場合、構造FEMと音響BEMを連成させます：
- 構造FEMで表面の振動速度を計算
- 音響BEMで表面速度から遠方音圧を計算
- あるいは音響FEMで車室内の音圧分布を計算
吸音材のモデリング
内装材などの吸音材は複素インピーダンス境界条件として近似します：
ここで $Z_s$ は表面音響インピーダンスです。インピーダンス管試験（ISO 10534）で実測した値をFEMに入力します。