固有振動
固有振動(こゆうしんどう、characteristic vibration, normal mode)とは対象とする振動系が自由振動を行う際、その振動系に働く特有の振動のことである。このときの振動数を固有振動数という。
振動の速さは単位時間に起こる往復運動の回数で表され、この回数を振動数または周波数という。単位はHz(ヘルツ)である。
振動の1回の往復運動は円運動1周に対応していて、振動の速さは単位時間におこなわれる円運動の回転角で表される。これを角振動数という。角振動数は振動数に1周の角度2π(rad)をかけて定義される。単位はrad/sである。
質量mの物体を一端を固定したばね定数kのばねの他端に取り付けて、摩擦の無い水平面上に置く。
右向きを正にx軸をとり、ばねが自然長の時の物体の位置を0とする。
物体を正の向きに移動させるとばねが伸び、負の向きに移動させるとばねは縮む。
いずれもばねはフックの法則に従うため、物体の変位をx、物体がばねから受ける力をFとすると
が成り立つ。また物体の加速度をxの時間tによる2階微分で表すと、
である。
(1-1)と(1-2)から
を得る。この2階微分方程式を解くと一般解は
となる。ただしA,\omega,\phiは定数で\omega =\sqrt{k\over m} である。
このときのωがばね-質量系の固有角振動数である。
単振り子は微小振動をしているとき水平面内で単振動をしているとみなすことができる。おもり(質点とみなす)の質量をm、糸の長さをℓとする。糸が鉛直線となす角度θが十分小さいとき、水平方向にx軸をとると変位は
物体の加速度をxの時間tによる2階微分で表すと、ニュートンの運動方程式は
である。(2-1)、(2-2)、(2-3)から
を得る。この2階微分方程式を解くと一般解は
となる。ただしA,\omega,\phiは定数で\omega =\sqrt{g\over l} である。このときのωが単振り子の固有角振動数である。
線密度ρ(kg/m)で張力T(N)で引っ張られている弦に関して、v=\sqrt {T/\rho}とおくと
の波動方程式を得る。この波動方程式を解くと、
このような各y_n(x,t)を基準モードという。また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は
(3-1)においてn=1,2,3の基準モードは右図のような振動を示す。
\omega_n={n\pi v\over l}={n\pi \over l}\sqrt {T\over \rho}
である。
空気の密度をρ(g/㎥)、体積弾性率をK(N/㎡)、v=\sqrt {K/\rho}とする。ここでは開口で実際に生じる開口端補正を無視して考える。
の波動方程式を得る。この波動方程式を解くと、
また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は
\omega_n={(2n-1)\pi v\over 2l}={(2n-1)\pi \over 2l}\sqrt {K\over \rho}
である。
の波動方程式を得る。この波動方程式を解くと、
また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は
\omega_n={n\pi v\over l}={n\pi \over l}\sqrt {K\over \rho}
である。
(1-2)と(1-5)から
(1-6)式でm\omega^2 =kを満足していれば解であることがいえる。
(2-4)と(2-6)から
(2-7)式で\omega^2 ={g\over 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方向成分は次のように表すことができる。
したがってy方向の力F_yは
ここでT\sin\theta_{(x+\delta x)}にテイラー級数展開を適用すると
δxは微小であるため2次以上の項を無視できる。よって
(3-2)を(3-1)に代入すると、
θ十分に小さいとき\sin\theta\approx \tan\thetaと近似できる。また\tan\theta=\frac{\partial y}{\partial x}と置き換えられるから
線分\delta sの質量は\rho \delta sであるから[ニュートンの運動方程式は
δyが小さいから\delta s\approx \delta x ,さらにv=\sqrt {T\over \rho}とおくと
の波動方程式を得る。
波動方程式を解くために、変数分離法を用いる。
関数y(x,t)がxの関数X(x)とtの関数T(t)の積の形で表されると仮定して
とおく。(3-5)を(3-4)に代入して整理し、両辺をX(x)T(t)でわると
このとき左辺はxのみの関数、右辺はtのみの関数であり、xとtは独立変数である。両辺が等しいということは両辺の値が定数であるということになる。この定数をKとおくと(3-6)から
と書きかえられる。
xについての方程式\frac{\partial}^2X(x)}{\partial x^2}-KX(x)=0… (3-7)を解く。
ⅰ)K=0のとき
となる。この微分方程式の一般解はX(x)=ax+bである。
ⅱ)K>0のとき
実数の定数k用いてK=k^2とすると
と表される。ここでX(x)=e^{\alpha x}とおくと、\frac{\partial}^2X(x)}{\partial x^2}={\alpha}^2e^{\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_1e^{kx}+C_2e^{-kx}も解である。k=\sqrt Kから
ⅲ)K<0のとき
実数の定数k用いてK=-k^2とすると
と表される。ここでX(x)=e^{\alpha x}とおくと、\frac{\partial}^2X(x)}{\partial x^2}={\alpha}^2e^{\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_1e^{ikx}+C_2e^{-ikx}も解である。k=\sqrt {-K}から
C_1(\cos{\sqrt{-K}x}+i\sin{\sqrt{-K}x})+C_2(\cos{\sqrt{-K}x}-i\sin{\sqrt{-K}x})=C_3\cos{\sqrt{-K}x}+C_4\sin{\sqrt{-K}x}(C_3=C_1+C_2,C_4=iC_1-iC_2はそれぞれ定数)
ⅰ)~ⅲ)から
K=0のとき…X(x)=ax+b… (3-11)
K>0のとき…X(x)=C_1e^{\sqrt Kx}+C_2e^{-\sqrt Kx} … (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)を与えると
(3-12)に条件(3-14)を与えると
(3-13)に条件(3-14)を与えると
X(x)=0は弦が振動していない様子を表すので、振動する弦の解は
である。
tについての方程式\frac{\partial}^2T(t)}{\partial t^2}-Kv^2T(t)=0… (3-8)を解く。xについての微分方程式を解いたとき、導いた解はK<0のときであった。よってここでもK<0のときのみを考える。実数の定数kを用いてK=-k^2とすると(3-8)は
と表される。この2階微分方程式を解くと一般解は
となる。ただし、C_5,\omega_n,\phi_nは定数で、\omega_n=kv={n\pi v\over l}である。
(3-15)、(3-17)から
また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は
断面積Sの円筒の中の空気の振動を考える。空気の密度をρ[g/㎥]、空気のx軸方向の変位をy(x,t)とする。大気圧をP_0とすると、位置xにおける圧力はP_0+\delta P(x,t)と表される。
この円筒の中のxとx+δxの微小部分について考える。空気が振動していないとき微小部分の体積はV=Sδxである。空気が振動したときの体積の変化は
と表される。空気の体積と圧力の間には
の関係が成り立つ。ここでKは体積弾性率である。(4-1)を(4-2)に代入すると
δx→0で
空気の断面にはそれぞれ圧力がはたらいている。xにおける断面にはたらく力は
x+δxにおける断面にはたらく力は
また微小部分の質量はm=\rho S\delta xであり、ニュートンの運動方程式を整理すると
x→0で
(4-3),(4-5)より
v=\sqrt {K\over \rho}とおくと
の波動方程式を得る。
「弦に関する波動方程式の解法」と同様にして変数分離法で波動方程式を解いていくと、xについての方程式は次の解を得る。
K=0のとき…X(x)=ax+b… (4-7)
K>0のとき…X(x)=C_1e^{\sqrt Kx}+C_2e^{-\sqrt Kx} … (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つまり{\partial y(l,t) \over \partial x}=0。したがって管の満たすべき条件は
y(0,t)=0 and {\partial y(l,t) \over \partial x}=0… (4-10)
である。(4-7)に条件(4-10)を与えると
(4-8)に条件(4-10)を与えると
(4-9)に条件(4-10)を与えると
X(x)=0は気柱が振動していない様子を表すので、振動する気柱の解は
である。また、「弦に関する波動方程式の解法」と同様にしてtについての方程式を解くと、
となる。ただし、C_5,\omega_n,\phi_nは定数で、\omega_n=kv={(2n-1)\over 2l}\pi vである。したがって
また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は
ここでは開口で実際に生じる開口端補正を無視して解きすすめる。両端が開口で長さlの管について考えると、両端開口による条件は
{\partial y(0,t) \over \partial x}=0 and {\partial y(l,t) \over \partial x}=0… (4-15)
である。(4-7)に条件(4-15)を与えると
(4-8)に条件(4-15)を与えると
(4-9)に条件(4-15)を与えると
X(x)=0は気柱が振動していない様子を表すので、振動する気柱の解は
である。また、「弦に関する波動方程式の解法」と同様にしてtについての方程式を解くと、
となる。ただし、C_5,\omega_n,\phi_nは定数で、\omega_n=kv={n \over l}\pi vである。したがって
また各y(x,t)は線形微分方程式の解であるから、それらの和もまた解である。したがって一般解は