第10講の導入
前々回、常微分方程式の解法についてお話ししました。ここでは、2階定数係数型の同次線型常微分方程式の解法の具体例として単振動と減衰振動について学んでいきましょう。必要に応じて第08講の内容を復習しながら読み進めていってください。
フックの法則とバネによる単振動
ここでは一番簡単な単振動の例としてバネによる単振動を議論します。まずはフックの法則の復習から始めましょう。外力を受けないバネにはたらく復元力の性質は、経験則として以下のようにまとめられます。これらをフックの法則と言います。
- フックの法則
バネの復元力Fは自然長からの変位xに比例してF=−kxで与えられる。ここで、k(>0)はバネの形状などから決まる係数でバネ定数と呼ばれる。
さて、ニュートンの運動方程式にF=−kxを代入すると、バネにつながれた質量mの物体の運動方程式は以下で与えられます。
md2xdt2=−kx
これは第08講で説明した2階の常微分方程式の形なので、この常微分方程式は解くことが出来ます。ω=√k/mとすれば、任意定数A、Bを用いて
x=Acosωt+Bsinωt
と書けます。あるいは、振幅aと初期位相δを用いて
x=acos(ωt+δ)
と書くことも出来ます。このとき、これらの文字はA=acosδ、B=−asinδで結ばれています。この運動は2次元平面内を一定の速さで回転する円運動の1つの軸への射影であるとみなすことも出来ます。
先ほどのωは振動数と呼ばれ、運動状態は周期T=2π/ωであらわされます。バネの場合は
T=2π√mk
です。これは高校物理で暗記させられていた人もいるかもしれません。実はこのような式変形の下で導出される式なのです。
振り子による単振動
別の例として、振り子による単振動を考えてみましょう。長さlのひもに質量mのおもりをくっつけた単振り子の運動を考えます。ひもは伸び縮みせず、重さが無視できるほど軽いとします。運動方程式はふれ角θを用いて次のように書くことが出来ます。
maθ=m(r¨θ+2˙r˙θ)=ml¨θ=−mgsinθ
ここで、ふれ角が十分小さい場合の運動を考えます。|θ|≪1のときの近似式sinθ≃θを利用すると、運動方程式は次のように書くことが出来ます。
¨θ=−glθ
これは単振動の運動方程式に他なりません。ω2=g/lとおくと、振り子の周期は
T=2π√lg
となります。
減衰振動
現実のバネや振り子は永遠に単振動を繰り返すことはなく、振幅はゆっくりと減少して最後には止まります。この減衰振動を記述するために、バネ定数kのバネにつながれた質量mの物体に速度に比例した抵抗力−2h˙xがはたらく状況を考えましょう。このとき、運動方程式は次のように書くことが出来ます。
m¨x=−kx−2h˙x
この式はγ=h/m、ω20=k/mとおくと次のように書き換えることが出来ます。
¨x+2γ˙x+ω20x=0
従って、これを解くことが次の目標です。解をx=e−γtf(t)と仮定すると方程式は
{¨f−2γ˙f+γ2f+2γ(˙f−γf)+ω20f}e−γt
これを整理すると¨f+(ω20−γ2)f=0が得られます。
減衰振動
粘性抵抗が小さくγ<ω0であるとき、ω1=√ω20−γ2とおくとfの式は¨f=−ω21fとなって単振動の方程式の形になるのでf=acos(ω1t+δ)となります。従って、x=ae−γtcos(ω1t+δ)となります。この解の振幅は緩和時間1/γで指数関数的に減少します。これを減衰振動と言います。
過減衰
粘性抵抗が大きくω0<γの場合はσ=√γ2−ω20とおくとfの式は¨f=σ2fとなります。このときの一般解はf=a1eσt+a2e−γtとなります。従ってx==a1e−(γ−σ)t+a2e−(γ+σ)t。このような運動は過減衰と言います。
臨界減衰
減衰振動と過減衰の間の状態、ω0=γのとき、¨f=0となるのでf=a+btとなります。従って、x=−γt(a+bt)となります。このような運動を臨界減衰と言います。
第10講のまとめ
最後に一言コメントをして第10講を終わりにしましょう。ここでは古典力学を具体例として単振動や減衰振動を解説しました。しかし、単振動や減衰振動は古典力学以外の分野でも様々なところで顔を出します。例えば、コイルとコンデンサーで構成されるLC 回路を流れる電流は単振動の形であらわすことができます。また、インスリンによる血糖値の制御は減衰振動の形であらわせることが知られています。このように、古典力学で学んだ事柄は他分野でも盛んに応用されているのです。