ラム波のソースを表示

'''ラム波'''（[[英語|英]]：''Lamb waves''）とは均質・等方な[[弾性体]]の薄板中を伝搬する[[弾性波]]伝搬モードの一種である。
周波数により位相速度が変化する速度分散性および無数の伝搬モードが存在する多モード性を有する。
1917年に[[Horace Lamb]]により[[分散関係]]が導かれた<ref>Lamb, H. "On Waves in an Elastic Plate." Proc. Roy. Soc. London, Ser. A 93, 114&ndash;128, 1917.</ref>。

ラム波は均質・等方な板中を伝搬するガイド波の中でも、その振動面が板表面に対して垂直であるものを指す。なお、振動面が板表面に対して平行なものは、'''SH波による板波'''と呼ばれている。

== 基礎方程式 ==
=== 分散関係式 ===
伝搬方向に対して垂直な方向の長さが無限であることを仮定する。このとき、真空中に置かれた厚さ2h、縦波速度c<sub>L</sub>、横波速度c<sub>T</sub>の[[等方性媒質|等方弾性体]]の板を伝搬するラム波の[[分散関係]](角周波数ωと波数kの関係)は以下のように表される。ここで波数kとは、ラム波の伝搬方向に対する波長をλとするとき、k=λ/(2π)で表される量である。

{{Indent|<math>\Omega_{\rm S}(\omega,k)\equiv\left(q^2-k^2\right)^2\left(\frac{\sin(qh)}{q}\right)\cos(ph)+4k^2p\sin(ph)\cos(qh)=0,</math>}}
{{Indent|<math>\Omega_{\rm A}(\omega,k)\equiv\left(q^2-k^2\right)^2\left(\frac{\sin(ph)}{p}\right)\cos(qh)+4k^2q\sin(qh)\cos(ph)=0.</math>}}

この式は特にレイリー=ラム周波数方程式と呼ばれている。ただし式中のp,qはそれぞれ縦波および横波の面外方向の波数成分を表し、

{{Indent|<math>p=\sqrt{\left(\frac{\omega}{c_{\rm L}}\right)^2-k^2},</math>}}
{{Indent|<math>q=\sqrt{\left(\frac{\omega}{c_{\rm T}}\right)^2-k^2},</math>}}

であり、実数あるいは純虚数をとる。レイリー=ラム周波数方程式は実数の角周波数ωに対して、一般に実数あるいは複素数の波数kを解に持つ。kが実数の場合、ラム波は長距離を伝搬しうるため、伝搬モードと呼ばれる。kが複素数あるいは純虚数の場合、伝搬に従いその振幅が指数関数的に減衰するため、非伝搬モードとなる。レイリー=ラム周波数方程式を満たすという意味ではいずれもラム波であるが、文脈によってはラム波伝搬モードをラム波と呼んでいることもある。レイリー=ラム周波数方程式はその板厚あるいは周波数の無限大の極限として[[レイリー波]]の分散関係式を含む。特に板厚が弾性波の波長に対して同程度であるときに上述したラム波としての性質（速度分散性・多モード性）が顕著に現れる。

=== 変位分布と対称・反対称モード ===
板の中立面上の一点を原点とし、伝搬方向をx方向、板厚方向をzとするとき、定数A, B, C, Dを用いて、ラム波の変位場は以下の式で表される。

{{Indent|<math>u_x=\operatorname{Re} \left \lbrack \left \lbrace kA\cos(pz)+qB\cos(qz)+kC\sin(pz)+qD\sin(qz)\right \rbrace \exp\left \lbrace i(kx-\omega t)\right \rbrace \right \rbrack,</math>}}
{{Indent|<math>u_z=\operatorname{Re} \left \lbrack i\left \lbrace pA\sin(pz)-kB\sin(qz)-pC\cos(pz)+kD\cos(qz)\right \rbrace \exp\left \lbrace i(kx-\omega t)\right \rbrace \right \rbrack ,</math>}}

A,B,C,Dは振幅を表す任意定数であり、

{{Indent|<math>
\begin{pmatrix}
2kp\sin(ph) & (q^2-k^2)\sin(qh) \\
(q^2-k^2)\cos(ph) & -2kq\cos(qh) 
\end{pmatrix}
\begin{pmatrix}
A \\ B 
\end{pmatrix}
= 0,</math>}}
{{Indent|<math>
\begin{pmatrix}
2kp\cos(ph) & (q^2-k^2)\cos(qh) \\
(q^2-k^2)\sin(ph) & -2kq\sin(qh) 
\end{pmatrix}
\begin{pmatrix}
C \\ D 
\end{pmatrix}
= 0,</math>}}

を満たす。上の分散関係式におけるΩ<sub>S</sub>=0を満たす場合、A,Bがともに0になる自明解以外の解をとりうる（A,Bがいずれかあるいは双方が非ゼロをとりうる）。このとき、変位場はz=0の面に対して対称となるため、このような変位場を有するラム波伝搬モードは特に'''対称モード'''と呼ばれる。

同様にΩ<sub>A</sub>=0を満たす場合、C,Dがともに0になる自明解以外の解をとりうる（C,Dがいずれかあるいは双方が非ゼロをとりうる）。このとき、変位場はz=0の面に対して反対称となるため、このような変位場を有するラム波伝搬モードは特に'''反対称モード'''と呼ばれる。

=== 位相速度・群速度と分散曲線 ===
ラム波の'''[[位相速度]]'''c<sub>p</sub>はc<sub>p</sub>=ω/kと表される。例えば単一周波数のラム波が薄板中を伝搬しているときに、表面における位相が伝搬方向に対して進む速さはこの位相速度c<sub>p</sub>に一致する。また、ラム波の角周波数ωを波数kの関数ω(k)とみなすことにより、c<sub>g</sub>=∂ω/∂kを定義できる。このc<sub>g</sub>はラム波の'''[[群速度]]'''と呼ばれ、ラム波の波束が伝搬する速度に対応する。

[[Image:Sym_asym_sigma0.27_und_0.34_edited2.svg|thumb|位相速度に関する分散曲線。青あるいは赤の線はそれぞれ対称および反対称モードを示している。また実線と点線は[[ポアソン比]]σの異なる2種類の材料の分散曲線を示している。]]

横軸に周波数を取り、縦軸に位相速度あるいは群速度をとった図は'''分散曲線'''と呼ばれ、ラム波の性質を知る上で重要な役割を果たす。図は位相速度に関する分散曲線を示している。横軸、縦軸はそれぞれ角周波数ωと板厚d(=2h)、位相速度c<sub>p</sub>（図中ではv）を横波速度c<sub>T</sub>（図中ではv<sub>s</sub>）で除することにより正規化している。青あるいは赤の線はそれぞれ対称および反対称モードを示している。また実線と点線は[[ポアソン比]]の異なる2種類の材料の分散曲線を示している。

== 性質 ==
有限の大きさの振動子から弾性波を励起した場合、伝搬距離方向以外の方向に拡散する。このため弾性波は一般に距離に従って減衰する。ラム波は力学的なエネルギーが板中に閉じ込められるため、無限媒質中を伝搬する縦波や横波に比べてこの効果による減衰が小さい。なお、このことは境界にそって伝搬する性質を有するガイド波の一般的な特徴である。

ラム波は特別な場合を除き、一般に強い [[速度分散性]]を有する。このため、広帯域な波を励起した場合伝搬に従い波形が変わる。また、単一の周波数でも複数のモードが存在しうる。

=== 極限 ===
0次の対称・反対称モードは任意の周波数でレイリー=ラム周波数方程式を満たし、伝搬モードとなる。0次モードの位相速度および群速度は周波数fあるいは板厚d(=2h)の無限大の極限において、レイリー波の伝搬速度に一致する。また、1次以上の対称・反対称モードの位相速度および群速度は、周波数fあるいは板厚d(=2h)の無限大の極限において横波の伝搬速度に一致する。

1次以上の対称・反対称モードはある周波数以上でのみレイリー=ラム周波数方程式を満たす。そこでレイリー=ラム周波数方程式を満たす解のうち、波数kのk→0の極限における周波数を'''カットオフ周波数'''と呼ぶ<ref>Graff, K. F. "Wave Motion in Elastic Solids," "Chapter 8 Wave propagation in plates and rods," Dover, New York, 1975</ref>。周波数が（おおよそ）このカットオフ周波数を超えるたびにラム波の伝搬モードが増えるため、実用上・解析上重要な周波数である。k→0の極限では対称モード、反対称モードの角周波数ω<sub>cutoff</sub>はそれぞれ、

{{Indent|<math>\cos\left(\frac{\omega_{\rm cutoff} h}{c_{\rm L}}\right)\sin\left(\frac{\omega_{\rm cutoff} h}{c_{\rm T}}\right)=0,</math>}}
{{Indent|<math>\sin\left(\frac{\omega_{\rm cutoff} h}{c_{\rm L}}\right)\cos\left(\frac{\omega_{\rm cutoff} h}{c_{\rm T}}\right)=0,</math>}}

を満たす。
なお、あるモードにおけるカットオフ周波数はあくまでk→0の極限であるため、厳密にはカットオフ周波数ω<sub>cutoff</sub>がこのラム波伝搬モードにおける最小の周波数になるとは限らない。つまり、同じラム波伝搬モードのなかで最も小さい周波数をω<sub>cr</sub>とするとき、ω<sub>cr</sub> < ω<sub>cutoff</sub>となる場合が存在する。

<!--　=== 定点 ===

== 歴史 ==

== 応用 ==
--->

== 脚注 ==
{{脚注ヘルプ}}
{{Reflist}}

== 参考文献 ==
* Viktorov, I. A. "Rayleigh and Lamb Waves: Physical Theory and Applications," Plenum Press, New York, 1967.
* Rose, J. L., "Ultrasonic Waves in Solid Media," Cambridge University Press, 1999.

== 関連項目 ==
* [[2022年のフンガ・トンガ噴火]]

{{デフォルトソート:らむは}}
[[Category:振動と波動]]
[[Category:応用力学]]
[[Category:固体力学]]
[[Category:弾性]]