楽しい科学(理論)チャンネル

単振動の一般解

今回はx軸上でフックの法則従う物体の運動についてその運動を見ていきたい。僕が一番好きな運動。2通りの解法をお見せしよう。いろいろな教本に乗っている解をx(t)=eλtと仮定する方法と、エネルギー保存則を使い方法だ。後者の解法は、楕円の方程式を極座標で表すことで自然に解けるためおすすめの解法だ。
 単振動を簡単な例を出すと、片方が固定されたばねをちょっとだけ引っ張って離すとびよーん、びよーんとなる。これがまさに単振動なのだ。実際は空気抵抗や摩擦により運動は減衰していくが今回はずっと運動を続けるものとする。この物体にはたらく力Fは原点から遠ざかるに比例して大きくなる。その比例係数をkとしている。なんではたらく力がマイナスなんだ、と思うかもしれない。符号はx軸正の方向を+とすると、x>0のときFx軸負(引っ張る)の方向に力がはたらく。x<0のときFx軸正の方向(押し返す)に力がはたらく。
 話がそれたが運動方程式を解いていこう。この物体の運動方程式は md2dt2x(t)=kx(t) 普通に積分して解ける形ではない。xを2回微分してxがでてきている。指数関数f(x)=eλx+μが使えそう。f(x)を微分してみよう。 ddxf(x)=aeλx+μ もういっちょ d2dx2f(x)=λ2eλx+μ=a2f(x) 運動方程式とかなり似た形になった。今回は解をx(t)=eλt+μと仮定しよう。これを運動方程式に代入する。 d2dt2x(t)=kmx(t) d2dt2eλx+μ=kmeλx+μ a2eλx+μ=kmeλx+μ (λ2+km)eλx+μ=0 実数範囲では解けないですな。eλx+μは常には0でないので λ=±ikm x(t)=eikmt+μ, eikmt+μ 解が複数ある時、その線形結合も運動方程式を満たすので、 x(t)=c1eikmt+μ+c2eikmt+μ x(t)=c1eμeikmt+c2eμeikmt うわー思ったより複雑、さらにオイラーの公式eiα=cosα+isinαを使う。 x(t)=c1eμcoskmt+ic1eμsinkmt +c2eμcoskmtic2eμsinkmt x(t)=(c1eμ+c2eμ)coskmt+i(c1eμc2eμ)sinkmt c1eμ+c2eμ, i(c1eμc2eμ), kmtに依らない定数であるため、これらをa,b,ωと再定義すると x(t)=acosωt+bsinωtこれが単振動の一般解とゆうやつである。A=a2+b2, δ=arctanba となるような定数A,δをとると、一般解は x(t)=a2+b2(aa2+b2cosωt+ba2+b2sinωt) =A(cosδcosωt+sinδsinωt)=Acos(ωtδ) このようにもかける。 x(t)=Acos(ωtδ) 今回の記事は短いので、おまけで力学的エネルギー保存則から一般解も示そう。 U(x)=12kx2 のような、ポテンシャルエネルギーを考える。 ddxU(x)=kx=F(x) を満たすので保存力である。仕事は W=x0xF(x)dx=x0xdU(x)dxdx=x0xdU(x) (1)W=U(x0)U(x)=12kx0212kx2 F=mv˙として仕事を計算する mdvdtdx=mdvdtdxdtdt=mdvdtvdt=12mdv2dtdt mdvdtdx=12mdv2 であるから、 W=x0xF(x)dx=x0xdvdtdx=12mv(x0)v(x)dv2 (2)W=12mv212mv02 (1)、(2)式から、 W=12mv212mv02=12kx0212kx2 12mv2+12kx2=E 力学的エネルギー保存則を示した。 v22Em2+x22Ek2=1 このように式変形すると、単振動の力学的エネルギー保存則は(x,v)平面上で楕円の周の方程式となる。すごい! (3)x=2Ekcosϕ (4)v=2Emsinϕ とかける媒介変数ϕが存在するはずだ。(3)式をtで微分すると、 ddtx=2Eksinϕddtϕ (4)式を代入すると、 ddtx=mkvddtϕ dx/dt=vなので mkϕ˙=1 両辺をtで積分する。 ϕ+ϕ0=km t k/m=ωとすると、単振動の一般解は(3)式に代入して、 x(t)=2mEωcos(ωt+ϕ0) となる。2mEω=A, ϕ0=δとすれば、微分方程式をごりごり解いた時の一般解と同じになる。