固有振動

太鼓の表面における固有振動

固有振動（こゆうしんどう、テンプレート:Lang-en）とは、ある系が自由振動を行う際に現れる、いくつかの特定の振動形式のことである。固有振動の振動数を固有振動数という。

代表的な振動系の固有振動

ばね‐質量系の固有振動

ばね‐質量系の振動

質量mの物体を一端を固定したばね定数kのばねの他端に取り付けて、摩擦の無い水平面上に置く。 右向きを正にx軸をとり、ばねが自然長の時の物体の位置を0とする。 物体を正の向きに移動させるとばねが伸び、負の向きに移動させるとばねは縮む。 いずれもばねはフックの法則に従うため、物体の変位をx、物体がばねから受ける力をFとすると テンプレート:Indent が成り立つ。また物体の加速度をxの時間tによる2階微分で表すと、 ニュートンの運動方程式は テンプレート:Indent である。 (1-1)と(1-2)から テンプレート:Indent を得る。この2階微分方程式を解くと一般解は テンプレート:Indent となる。ただし$A,\omega ,\varphi$は定数で$\omega =\sqrt{k/m}$である。 このときのωがばね-質量系の固有角振動数である。

単振り子の固有振動

単振り子の様子

単振り子は微小振動をしているとき水平面内で単振動をしているとみなすことができる。おもり(質点とみなす)の質量をm、糸の長さをℓとする。糸が鉛直線となす角度θが十分小さいとき、水平方向にx軸をとると変位は テンプレート:Indent 水平方向の力は テンプレート:Indent 物体の加速度をxの時間tによる2階微分で表すと、ニュートンの運動方程式は テンプレート:Indent である。(2-1)、(2-2)、(2-3)から テンプレート:Indent テンプレート:Indent を得る。この2階微分方程式を解くと一般解は テンプレート:Indent となる。ただし$A,\omega ,\varphi$は定数で$\omega =\sqrt{g/l}$である。このときのωが単振り子の固有角振動数である。

弦の固有振動

線密度ρ(kg/m)で張力T(N)で引っ張られている弦に関して、$v=\sqrt{T/\rho }$とおくと テンプレート:Indent の波動方程式を得る。この波動方程式を解くと、 テンプレート:Indent このような各$y_n(x,t)$を基準モードという。また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は テンプレート:Indent (3-1)においてn=1,2,3の基準モードは右図のような振動を示す。

n=1のとき第1調和振動

n=2のとき第2調和振動

n=3のとき第3調和振動

またこの系における固有角振動数は

${\omega }_n=\frac{n\pi v}{l}=\frac{n\pi }{l}\sqrt{\frac{T}{\rho }}$

である。

気柱の固有振動

空気の密度をρ(g/㎥)、体積弾性率をK(N/㎡)、$v=\sqrt{K/\rho }$とする。ここでは開口で実際に生じる開口端補正を無視して考える。

一端が閉口で他端が開口の管

テンプレート:Indent の波動方程式を得る。この波動方程式を解くと、 テンプレート:Indent また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は テンプレート:Indent この系における固有角振動数は

${\omega }_n=\frac{(2n-1)\pi v}{2l}=\frac{(2n-1)\pi }{2l}\sqrt{\frac{K}{\rho }}$

である。

両端が開口の管

テンプレート:Indent の波動方程式を得る。この波動方程式を解くと、 テンプレート:Indent また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は テンプレート:Indent この系における固有角振動数は

${\omega }_n=\frac{n\pi v}{l}=\frac{n\pi }{l}\sqrt{\frac{K}{\rho }}$

である。

付録

(1-4)式が(1-3)式の解であることの証明

テンプレート:Indent テンプレート:Indent (1-2)と(1-5)から テンプレート:Indent (1-6)式で$m{\omega }^2=k$を満足していれば解であることがいえる。

(2-5)式が(2-4)式の解であることの証明

テンプレート:Indent テンプレート:Indent (2-4)と(2-6)から テンプレート:Indent (2-7)式で${\omega }^2=\frac{g}{l}$を満足していれば解であることがいえる。

弦に関する波動方程式

振動する弦の微小部分

波動方程式の導出

　線密度ρ(kg/m)で張力T(N)で引っ張られている弦がXY平面上にあるとする。その弦のxとx+δxの微小部分について考える。位置xにおける弦の接線とx軸のなす角を${\theta }_x$、位置x+δxにおける弦の接線とx軸のなす角を${\theta }_{x+\delta x}$とすると張力$T_A$と$T_B$のx方向成分、y方向成分は次のように表すことができる。 テンプレート:Indent テンプレート:Indent テンプレート:Indent テンプレート:Indent 　したがってy方向の力$F_y$は テンプレート:Indent ここで$T\sin {\theta }_{(x+\delta x)}$にテイラー級数展開を適用すると テンプレート:Indent

δxは微小であるため2次以上の項を無視できる。よって テンプレート:Indent (3-2)を(3-1)に代入すると、 テンプレート:Indent θ十分に小さいとき$\sin \theta \approx \tan \theta$と近似できる。また$\tan \theta =\frac{\partial y}{\partial x}$と置き換えられるから テンプレート:Indent 線分$\delta s$の質量は$\rho \delta s$であるから[ニュートンの運動方程式は テンプレート:Indent δyが小さいから$\delta s\approx \delta x$ ,さらに$v$=$\sqrt{\frac{T}{\rho }}$とおくと テンプレート:Indent の波動方程式を得る。

波動方程式の解法

波動方程式を解くために、変数分離法を用いる。 関数y(x,t)がxの関数X(x)とtの関数T(t)の積の形で表されると仮定して テンプレート:Indent とおく。(3-5)を(3-4)に代入して整理し、両辺をX(x)T(t)でわると テンプレート:Indent このとき左辺はxのみの関数、右辺はtのみの関数であり、xとtは独立変数である。両辺が等しいということは両辺の値が定数であるということになる。この定数をKとおくと(3-6)から テンプレート:Indent テンプレート:Indent と書きかえられる。

• xについての方程式$\frac{d^{2}X(x)}{d{x}^2}-KX(x)=0$… (3-7)を解く。

ⅰ)K＝0のとき

テンプレート:Indent となる。この微分方程式の一般解は$X(x)=ax+b$である。

ⅱ)K>0のとき

実数の定数k用いて$K=k^2$とすると テンプレート:Indent と表される。ここで$X(x)=e^{\alpha x}$とおくと、$\frac{d^{2}X(x)}{d{x}^2}={\alpha }^2{e}^{\alpha x}$なので(3-9)は$({\alpha }^2-k^2)X(x)=0$と書きかえられる。X(x)は任意の関数であるから${\alpha }^2-k^2=0$を考える。つまり$\alpha =\pm k$である。したがって解は$X(x)=e^{kx}$と$X(x)=e^{-kx}$であり、またその線形結合の$X(x)=C_1{e}^{kx}+C_2{e}^{-kx}$も解である。$k=\sqrt{K}$から テンプレート:Indent

ⅲ)K<0のとき

実数の定数k用いて$K=-k^2$とすると テンプレート:Indent と表される。ここで$X(x)=e^{\alpha x}$とおくと、$\frac{d^{2}X(x)}{d{x}^2}={\alpha }^2{e}^{\alpha x}$なので(3-10)は$({\alpha }^2+k^2)X(x)=0$と書きかえられる。X(x)は任意の関数であるから${\alpha }^2+k^2=0$を考える。つまり$\alpha =\pm ik$である。したがって解は$X(x)=e^{ikx}$と$X(x)=e^{-ikx}$であり、またその線形結合の$X(x)=C_1{e}^{ikx}+C_2{e}^{-ikx}$も解である。$k=\sqrt{-K}$から テンプレート:Indent オイラーの公式を適用すると テンプレート:Indent($C_3=C_1+C_2,C_4=i{C}_1-i{C}_2$はそれぞれ定数)

ⅰ)～ⅲ)から

K=0のとき…$X(x)=ax+b$… (3-11)

K>0のとき…$X(x)=C_1{e}^{\sqrt{K}x}+C_2{e}^{-\sqrt{K}x}$ … (3-12)

K<0のとき…$X(x)=C_3\cos \sqrt{-K}x+C_4\sin \sqrt{-K}x$… (3-13)

両端固定の長さ$l$の弦について考えると、両端固定による条件は

$y(0,t)=0$ and $y(l,t)=0$… (3-14)

(3-11)に条件(3-14)を与えると テンプレート:Indent (3-12)に条件(3-14)を与えると テンプレート:Indent (3-13)に条件(3-14)を与えると テンプレート:Indent $X(x)=0$は弦が振動していない様子を表すので、振動する弦の解は テンプレート:Indent である。

• tについての方程式$\frac{d^{2}T(t)}{d{t}^2}-K{v}^{2}T(t)=0$… (3-8)を解く。xについての微分方程式を解いたとき、導いた解はK<0のときであった。よってここでもK<0のときのみを考える。実数の定数kを用いて$K=-k^2$とすると(3-8)は

テンプレート:Indent と表される。この2階微分方程式を解くと一般解は テンプレート:Indent となる。ただし、$C_5$,${\omega }_n$,${\varphi }_n$は定数で、${\omega }_n=kv=\frac{n\pi v}{l}$である。

(3-15)、(3-17)から テンプレート:Indent また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は テンプレート:Indent

気柱に関する波動方程式

波動方程式の導出

断面積Sの円筒の中の空気の振動を考える。空気の密度をρ[g/㎥]、空気のx軸方向の変位をy(x,t)とする。大気圧を$P_0$とすると、位置xにおける圧力は$P_0+\delta P(x,t)$と表される。

気柱の変位

この円筒の中のxとx+δxの微小部分について考える。空気が振動していないとき微小部分の体積はV=Sδxである。空気が振動したときの体積の変化は テンプレート:Indent と表される。空気の体積と圧力の間には テンプレート:Indent の関係が成り立つ。ここでKは体積弾性率である。(4-1)を(4-2)に代入すると テンプレート:Indent δx→0で テンプレート:Indent

気柱にはたらく圧力

空気の断面にはそれぞれ圧力がはたらいている。xにおける断面にはたらく力は テンプレート:Indent x+δxにおける断面にはたらく力は テンプレート:Indent したがって微小部分にはたらく力は テンプレート:Indent

また微小部分の質量は$m=\rho S\delta x$であり、ニュートンの運動方程式を整理すると テンプレート:Indent x→0で テンプレート:Indent (4-3),(4-5)より テンプレート:Indent $v$=$\sqrt{\frac{K}{\rho }}$とおくと テンプレート:Indent の波動方程式を得る。

波動方程式の解法

「弦に関する波動方程式の解法」と同様にして変数分離法で波動方程式を解いていくと、xについての方程式は次の解を得る。

K=0のとき…$X(x)=ax+b$… (4-7)

K>0のとき…$X(x)=C_1{e}^{\sqrt{K}x}+C_2{e}^{-\sqrt{K}x}$ … (4-8)

K<0のとき…$X(x)=C_3\cos \sqrt{-K}x+C_4\sin \sqrt{-K}x$… (4-9)

一端が閉口で他端が開口の管の場合

ここでは開口で実際に生じる開口端補正を無視して解きすすめる。左端が閉口で右端が開口な長さ$l$の管について考えると、左端が閉口による条件は$y(0,t)=0$、右端が開口による条件は$P(l,t)=0$つまり$\frac{\partial y(l,t)}{\partial x}=0$。したがって管の満たすべき条件は

$y(0,t)=0$ and $\frac{\partial y(l,t)}{\partial x}=0$… (4-10)

である。(4-7)に条件(4-10)を与えると テンプレート:Indent (4-8)に条件(4-10)を与えると テンプレート:Indent (4-9)に条件(4-10)を与えると テンプレート:Indent $X(x)=0$は気柱が振動していない様子を表すので、振動する気柱の解は テンプレート:Indent である。また、「弦に関する波動方程式の解法」と同様にしてtについての方程式を解くと、 テンプレート:Indent となる。ただし、$C_5$,${\omega }_n$,${\varphi }_n$は定数で、${\omega }_n=kv=\frac{(2n-1)}{2l}\pi v$である。したがって テンプレート:Indent また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は テンプレート:Indent

両端が開口の管の場合

ここでは開口で実際に生じる開口端補正を無視して解きすすめる。両端が開口で長さ$l$の管について考えると、両端開口による条件は

$\frac{\partial y(0,t)}{\partial x}=0$ and $\frac{\partial y(l,t)}{\partial x}=0$… (4-15)

である。(4-7)に条件(4-15)を与えると テンプレート:Indent (4-8)に条件(4-15)を与えると テンプレート:Indent (4-9)に条件(4-15)を与えると テンプレート:Indent $X(x)=0$は気柱が振動していない様子を表すので、振動する気柱の解は テンプレート:Indent である。また、「弦に関する波動方程式の解法」と同様にしてtについての方程式を解くと、 テンプレート:Indent となる。ただし、$C_5$,${\omega }_n$,${\varphi }_n$は定数で、${\omega }_n=kv=\frac{n}{l}\pi v$である。したがって テンプレート:Indent また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は テンプレート:Indent

参考文献

• テンプレート:Cite book - 原タイトル：The Physics of Musical Instruments
• テンプレート:Cite book
• 気柱の振動

関連項目

• 共振
• 共鳴
• 固有関数