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

運動方程式

/*ニュートン力学を勉強する意味 残念ながら、ニュートン力学ではミクロな世界や極端に重い物体の付近や光速に近い物理現象を説明することはできない。しかし、身の回りで起きている雨滴の落下から、惑星の運動までかなり広い範囲の運動を説明することができる。更に量子力学を勉強するために、エネルギーや運動量について知っておく必要があるし、相対性理論は、ニュートン力学に補正を入れた形をしているため、まずは、ニュートン力学を勉強することが必要なのだ。*/
本来物体は大きさがあり、物体の運動には重心や回転などの要素があるが、ニュートン力学で扱う物体は全て質量をもつ点として扱う。ニュートン力学第1法則の準備として位置と速度をベクトルで定義しておく。
図のような座標系\(K\)で原点から見た点\(P\)の位置を \[\boldsymbol{x}:=\begin{pmatrix}x\\y\\z\end{pmatrix}\] と定める。空間上で場所を表すには、縦、横、高さの3つを指定する必要がある。地図のような平面では、縦(緯度)、横(経度)の2つを指定する必要があるので、\(\boldsymbol{x}=(x,y)\)と表せる。さらに点\(P\)が動いている場合、\(P\)の位置は時間の関数となるので、 \[\boldsymbol{x}(t)=\begin{pmatrix}x(t)\\y(t)\\z(t)\end{pmatrix}\] と表す。点\(P\)の速度を \[\boldsymbol{v}(t):=\frac{d}{dt}\boldsymbol{x}(t)\] と定義する。\(\boldsymbol{v}(t)=\dot{\boldsymbol{x}}(t)\)とも表す。位置をベクトルで表したため、各成分の時間微分でも速度を表すことができる。ニュートン力学のゴールは、\(\boldsymbol{x}=\boldsymbol{x}(t)\)位置を時間で表すことである。これによって何時何分に物体がどこにあるということを知ることができる。しかしそれが難しい場合は妥協案として、何かのパラメータ\(s\)を用いて、\(\boldsymbol{x}=\boldsymbol{x}(s)\)のように軌道(軌跡)を求めることもある。例として、一定の速度、\(\boldsymbol{v}_0=(v_x,v_y,v_z)\)で進む物体の運動の解は、 \[\boldsymbol{v}_0=\frac{d}{dt}\boldsymbol{x}(t)\] \(x\)成分についてみると、 \[v_x=\frac{d}{dt}x(t)\] 両辺を\(t\)で積分すると、 \[\int v_xdt=\int\frac{dx(t)}{dt}dt\] \[x(t)=v_xt+c_{x}\] \(c_x\)は積分定数、\(y,z\)成分にいても同様なので、\(\boldsymbol{x}(t)=\boldsymbol{v}t+\boldsymbol{c}\)の等速直線運動になる。(静止も\(\boldsymbol{v}=0\)の等速直線運動とみなせる。)普段の経験として我々が見ている光景では、空気抵抗や摩擦があるせいで、ぱっと見力が加わっていない状態でも速度が徐々に減衰して見えるが、理想的には、力を受けいていない物体は静止もしくは等速直線運動をするのである。しかし、力を受けいていない物体が等速直線運動するためには、その物体が置かれている環境(系)も静止、等速直線運動している必要がある。いきなり車(系)が急発進した場合、車の中のものが、静止していたのに発進した方とは逆向きに動き出すことがある。この現象の解決方法は車の外から見ることで、外から見ると、ものが急に動き出すわけではなく、車が急発進して、ものがその場に置いて行かれているように見えるはずである。そもそも力がはたらかない物体が等速直線運動するためには、系にも制限がかかっている。このニュアンスでニュートン力学第一法則と言っている教本が多いように思える。
**ニュートン力学第1法則**
力を受けていない物体は静止または等速直線運動する。
 他の物体や、力を受けていない物体が静止または等速直線運動するという法則であるが、真の意味としては力を受けいていない物体が等速直線運動するためには、加速度系であってはならないということである。等速直線運動する系を慣性系\(K\)という。座標系\(K\)から見て\(\boldsymbol{r_0}\)離れていて一定の速度\(\boldsymbol{V}\)で移動する座標系を\(K'\)とする。\(K\)での物体の位置を\(\boldsymbol{x}(t)\)、\(K'\)での物体の位置を\(\boldsymbol{x'}(t)\)とすると、 \[\boldsymbol{x}(t)=\boldsymbol{x}'(t)+\boldsymbol{r_0}+\boldsymbol{V}t\tag{1}\] と表すことができる。両辺の時間微分は \[\dot{\boldsymbol{x}}(t)=\dot{\boldsymbol{x}}'(t)+\boldsymbol{V}\tag{2}\]
上の図は\(K\)から見て物体の速度が\(\boldsymbol{V}\)であったときである。(2)式に\(\dot{\boldsymbol{x}}(t)=\boldsymbol{V}\)を代入すると、 \[\boldsymbol{V}=\dot{\boldsymbol{x}}'(t)+\boldsymbol{V},\ \dot{\boldsymbol{x}}'(t)=\boldsymbol{0}\] つまり\(K'\)から見ると物体が静止して見える。次に\(K\)から見て物体の任意の一定の速度\(\boldsymbol{v}\)で動いているとき、(2)式に\(\dot{\boldsymbol{x}}(t)=\boldsymbol{v}\)を代入すると、 \[\boldsymbol{v}=\dot{\boldsymbol{x}}'(t)+\boldsymbol{V},\ \dot{\boldsymbol{x}}'(t)=\boldsymbol{v}-\boldsymbol{V}\] \(K'\)から見て、物体の速度が\(\boldsymbol{v}-\boldsymbol{V}\)の等速直線運動をしているように見える。\(K\)が慣性系であるとき、一定の速度\(\boldsymbol{V}\)で移動する座標系\(K'\)も第1法則を満たす。一般に慣性系\(K\)に対して(1)式のような座標変換で得た座標系も慣性系となる。
 次に第2法則の準備として運動量と加速度を定義する。 物体の質量\(m\)、速度を\(\boldsymbol{v}(t)\)とする。 物体の運動量を次のように定義する。 \[\boldsymbol{p}(t):=m\boldsymbol{v}(t)\] 加速度を次のように定義する。 \[\boldsymbol{a}(t):=\frac{d}{dt}\boldsymbol{v}(t)\] \(\boldsymbol{a}(t)=\dot{\boldsymbol{v}}(t)=\ddot{\boldsymbol{x}}(t)\)のような表記もされる。
**ニュートン力学第2法則**
慣性系において物体に加わる力\(\boldsymbol{F}\)を次のように定義する。 \[\boldsymbol{F}(t):=\frac{d}{dt}\boldsymbol{p}(t)\]
上記の方程式を運動方程式という。運動量\(p\)を質量\(m\)と、速度\(v\)の積で定義したのは、より高速であるまたはより質量の大きい物体の運動量を時間変化させるためにはより大きな力が必要になるようにするためである。質量は時間に依らない場合、 \[\boldsymbol{F}(t)=\frac{d}{dt}(m\boldsymbol{v}(t))=m\frac{d}{dt}\boldsymbol{v}(t)=m\boldsymbol{a}(t)\] と書き換えることもできる。力とは質量と加速度の積でも表すことができる。第1法則の時に紹介した(2)式の座標変換ガリレイ変換という。運動方程式にガリレイ変換をしてみよう。\(\boldsymbol{v}(t)\to \boldsymbol{v}(t)+\boldsymbol{V}\)と変換した運動方程式は \[\boldsymbol{F}(t)=m\frac{d}{dt}(\boldsymbol{v}(t)-\boldsymbol{V})\] となる。少し丁寧に計算すると、 \[\boldsymbol{F}(t)=m\lim_{\varDelta t\to 0}\frac{(\boldsymbol{v}(t+\varDelta t)-\boldsymbol{V})-(\boldsymbol{v}(t)-\boldsymbol{V})}{\varDelta t}\] \[\boldsymbol{F}(t)=m\lim_{\varDelta t\to 0}\frac{\boldsymbol{v}(t+\varDelta t)-\boldsymbol{v}(t)}{\varDelta t}\] \[\boldsymbol{F}(t)=m\frac{d}{dt}\boldsymbol{v}(t)\] 運動方程式はガリレイ変換に対して不変である。座標の取り方によって式の形は変わるかも知れないが、運動方程式は任意の慣性系で成り立つことを意味する。 \[\boldsymbol{a}(t)=\frac{d}{dt}\boldsymbol{v}(t)=\frac{\boldsymbol{F}(t)}{m}\] と変形すれば同じ力を受けている物体の加速度は質量に反比例することがわかる。
**ニュートン力学第3法則**
物体1物体2の間に相互作用がはたらくとき、物体1から物体2にはたらく力を\(\boldsymbol{F}_{1\to 2}\)と物体2から物体1にはたらく力を\(\boldsymbol{F}_{2\to 1}\)は同一作用線上で逆向きにはたらく。 \[\boldsymbol{F}_{1\to 2}=-\boldsymbol{F}_{2\to 1}\]
例をだすと、地上に接している物体の重力と垂直抗力や、質量のある物体どうしにはたらく万有引力などがある。より一般にある系で、\(n\)個の物体が相互作用しあう場合、 \[\sum_{i=1,i\neq1}^{n-1}\boldsymbol{F}_{1\to i}+\sum_{i=1,i\neq 2}^{n-1}\boldsymbol{F}_{2\to i}+\sum_{i=1,i\neq 3}^{n-1}\boldsymbol{F}_{3\to i}+\cdots+\sum_{i=1}^{n-1}\boldsymbol{F}_{n\to i}=\boldsymbol{0}\] となる。