物理/運動の法則の応用
提供: Internet Web School
(→回転運動の方程式 ) |
(→固定軸の周りの剛体の回転運動の方程式) |
||
820 行: | 820 行: | ||
==== 固定軸の周りの剛体の回転運動の方程式==== | ==== 固定軸の周りの剛体の回転運動の方程式==== | ||
+ | 回転運動の運動方程式から、任意の軸の周りの回転運動の方程式が簡単に導出できる。<br/> | ||
+ | z軸周りの場合を例にとり、説明する。<br/> | ||
+ | z軸周りの回転力はT→ez=→N⋅→ezなので、<br/> | ||
+ | 回転運動の方程式から<br/> | ||
+ | $T_{\vec e_z}=\vec N \cdot \vec e_z | ||
+ | =\frac{d\vec L}{dt} \cdot \vec e_z$<br/> | ||
+ | この式の右辺に,L=∑i→ri×→pi を代入すると<br/> | ||
+ | 右辺<br/> | ||
+ | =d∑i→ri×→pidt⋅→ez 微分の加法性から <br/> | ||
+ | =(∑i d→ri×→pidt)⋅→ez 内積の加法性から <br/> | ||
+ | =∑i(d→ri×→pidt⋅→ez) ベクトル積の性質8から <br/> | ||
+ | =∑i(dridt×→pi+→ri×d→pidt | ||
==== てこの原理と力のモーメント==== | ==== てこの原理と力のモーメント==== |
2014年12月25日 (木) 10:30時点における版
物理 > 運動の法則の応用
目次[非表示] |
解説
運動の3法則、万有引力の法則と力の法則を用いると、分子から銀河まであらゆる物体の運動を求めることが出来きる。
その正しさは地上の物体や人工衛星、惑星の運動などで確かめられている。
しかし、もっとはるかかなたの宇宙でもこれ等の法則は正しいのだろうか。
天体観測は、世界各地で行われ、年々新しい発見がされているが、現在のところ、この理論が間違っていることを示す観測結果は、得られていない。
そこで、これらの法則は宇宙の全体を支配しているものと、現在は信じられている。
運動の3法則からはエネルギー保存則や運動量保存則などの重要な保存則を導く事が出来る。
これらの保存則は、色々な運動を調べるとき、大変役立つ。これらについては次節で学ぶ。
質点の色々な運動
最初に最も簡単な運動から考える。
それは質点とみなせる物体の運動である。
質点の落体運動
地球上の物体は高いところから落とすと、時間とともに速度を増しながら落下する。
質点とみなせる物体の落下運動を、運動法則と力の法則を用いて、解析しよう。
質点の質量をmとすると、そこに作用する重力による力は、
真下(厳密には地球の重心;後で学ぶ)の方向・向きに大きさMgである。
落下の向きを負にした落下方向の一次元座標を考えると、重力加速度は−gで、質点mに作用する力は−mgである。
落下の加速度をαと置くと、運動の第2法則よりmα=−mg.
ゆえに質点の落下加速度αは負の重力加速度−gに等しい。
tで微分して−gとなる関数は−gt+cなので、質点の速度は−gt+cである。
ここでcは定数で、初期時刻0における質点の速度であり、初期速度と呼ばれる。
微分して−gt+cとなる関数を求めれば質点の位置x(t)=−12gt2+ct+dが得られる。
ここで、dは定数で初期時刻0での質点の位置(高さ)である。
これはガリレオが明らかにした落体法則である。
参考文献;
投射体の運動
質点を地面に対して角度θ(ラジアン)、速さuで投げたときの、質点はどのような運動を行うだろうか。
ガリレオは、慣性法則と落体の法則を組み合わせて利用して、放物線を描いて飛ぶことを発見した。
ニュートン力学を用いれば、運動の第2法則と質点に働く力(重力)から、以下のように、この運動を導ける。
適切な座標系をいれる
質点が投げ出された場所を原点とし、飛んでいく方向に地面と水平に引いた半直線をx軸の正の側に、地面と直角で上方に向かう半直線をy軸の正の側とする座標を定める。図参照。
質点に作用する力を求める
空気抵抗を無視すれば、質点に作用する力は、地球からの重力だけである。この力は、質点の質量をM,重力加速度をgとすると、質点の位置に関係なく常に、→F=(o,−Mg)である。
運動の第2法則から質点の運動方程式をつくる
質点の位置ベクトルを→r=(x,y)で表すと
運動方程式は、M(d2/dt2)→r(t)=→Fである。
座標成分表示すると
M(d2/dt2)x(t)=0, M(d2/dt2)y(t)=−Mg
運動の初期状態の指定
投げ上げた瞬間を時刻t=0とおくと、質点の初期位置は→r(0)=(0,0), 初期速度は→v(0)=(ucosθ,usinθ)
運動方程式を初期状態を使って解く
(1)x成分の式を解く
M(d2/dt2)x(t)=0は、M(d/dt)vx(t)=0なので(d/dt)vx(t)=0。 tで微分して零となるtの関数は定数なのでaと書くと、vx(t)=a
速度の定義より、(d/dt)x(t)=vxなので、(d/dt)x(t)=a. tで微分してaとなるのはat+b(bは未知定数)なので、x(t)=at+b
初期条件から、a=vx(0)=ucosθ, またx(0)=a0+b=0なのでb=0。
故に、x(t)=(ucosθ)t
(2)y成分の式を解く
M(d2/dt2)y(t)=−Mgは、(d/dt)vy(t)=−g tで微分して−gとなる関数は−gt+c(cは未知定数)なので、
vy(t)=−gt+c 故に(d/dt)y(t)=−gt+c
tで微分して−gt+cとなる関数は、−12gt2+ct+dなので、y(t)=−12g2t+ct+d
初期速度の条件から、c=−g0+c=vy(0)=usinθ d=−12g0+c0+d=y(0)=0
故に、y(t)=−12gt2+(usinθ)t
(3)運動の軌跡(xとyとの関係式)を求める
x(t)の式からt=x(t)/(ucosθ)
これをy(t)=−12gt2+(usinθ)tに代入すると
y(t)=(−g/2u2cos2θ)x2(t)+(tanθ)x(t)
これは上に凸な放物線である。
参考文献は
- ウィキブックス(高等学校理科 物理I 運動とエネルギー)の2.4.1 ニュートン方程式
惑星運動
前述のようにケプラーは、火星と太陽の観測データをユークリッド幾何学を巧みに利用して分析し次の惑星運動の3法則を発見した。
惑星運動の3法則を運動の第2法則と万有引力の法則から導く
この3法則は、運動の第2法則と万有引力の法則から導くことが出来るが少し難しい数学が必要である。大学で学ぶ。
惑星の軌道を太陽を中心とする円運動に限定すると、高校の数学の知識で3法則を導ける。
この場合ケプラーの第一法則は、仮定から、明白なので、第二法則から始める。
ケプラーの第2法則の導出
第二法則は、太陽と惑星を結ぶ動径の単位時間に掃く面積が一定であることを主張する。円運動のばあい、これは等速円運動であることと同じである。
そこで等速円運動であることを導こう。
太陽と惑星は質点として扱い、質量をそれぞれM,mとする。
惑星の軌道面をxy平面にし、太陽をその原点にとる。円運動の半径をr, 太陽と時刻tにおける惑星を結ぶ線分が、x軸となす角度をθ=θ(t)とおく。
惑星Pの位置;→r(t)=r(cosθ(t),sinθ(t))
惑星の速度;→v(t)=d→r(t)/dt=r(dcosθ(t)/dt,dsinθ(t)/dt)
=r(−sinθ(t)dθ(t)dt,cosθ(t)dθ(t)dt)
=rdθ(t)dt(−sinθ(t),cosθ(t))
惑星の加速度;→α(t)=d→v(t)/dt=r(d2θ(t)/dt2)(−sinθ(t),cosθ(t))
+r(dθ(t)/dt)(−cosθ(t)dθ(t)dt,−sinθ(t)dθ(t)dt)
=r(d2θ(t)/dt2)(−sinθ(t),cosθ(t))−r(dθ(t)dt)2(cosθ(t),sinθ(t))
惑星に働く力;万有引力の法則より、太陽の方向に向いた、大きさGMm/r2の力なので
→F(t)=−(GMm/r2)(cosθ(t),sinθ(t))
と表せる。
この力が、惑星の運動を変化させ、上述の加速度を生じさせたのだから、運動の第2法則m→α(t)=→F(t)より、
mr(d2θ(t)/dt2)(−sinθ(t),cosθ(t))−mr(dθ(t)dt)2(cosθ(t),sinθ(t)
=−(GMm/r2)(cosθ(t),sinθ(t))
変形すると、
mr(d2θ(t)/dt2)(−sinθ(t),cosθ(t))
=(mr(dθ(t)dt)2−GMm/r2)(cosθ(t),sinθ(t)) −−−−−−(1)
(−sinθ(t),cosθ(t)) と(cosθ(t),sinθ(t))は直交するベクトルなので、(1)式が成立する必要十分条件は、
d2θ(t)/dt2=0−−−−−−(2),
mr(dθ(t)dt)2−GMm/r2=0−−−−−−(3)
である。
(2)式から、角速度ω(t)=dθ(t)dt=ω0(定数)が
(3)式から、mr(dθ(t)dt)2=GMm/r2が
得られる。
これらより、惑星は等角速度
ω0=±√GM/r3 −−−−−− (4)
で太陽の周りを回転することが分かり、ケプラーの第2法則が得られた。
ケプラーの第3法則の導出
惑星が太陽の周りを一周する時間T(周期という)は、T=2π/ω0なので、(4)式より、
T=2π/√GM/r3=2π√r3/GM,
故にT2=4π2r3/GM,
T2/r3=4π2/GM
これは軌道が円の場合のケプラーの第3法則である。
万有引力の法則を,ケプラーの法則と運動の第2法則から導く
惑星が太陽の周りを円運動しているとき、太陽が惑星に及ぼしている力を計算する。
ケプラーの第2法則より、円運動する惑星は角速度一定である。これをω0とする。
太陽の位置を原点とし円の半径をrとすると、この惑星の加速度は→α(t)=−r(dθ(t)/dt)2(cosθ(t),sinθ(t))=−rω20(cosθ(t),sinθ(t)) 。これは、太陽にむかう大きさrω20のベクトル。
運動の第2法則より、惑星に働く力→Fは、太陽の方向に、大きさmrω20
ここで、m は惑星の慣性質量である。
ω20をrの関数で表すためケプラーの第3法則と用いる。
惑星の公転周期Tと円の半径rの間にはT2/r3=C,C;定数
T=2π/ω0なので
(2π/ω0)2/r3=C∴ω20=4π2/(Cr3)
それゆえ、力の大きさは
mrω20=4π2Cmr2
さらに、太陽の質量Mがk倍になると、質量Mの太陽がk個あり、それぞれが惑星に上記の力を与えると考えられる。
すると惑星に働く力はk倍になるので力の比例部分4π2Cは太陽の質量Mに比例することが分かる。
比例定数をGとおくと、4π2C=GM
従って惑星に働く力の大きさは、太陽の方向に、
GMmr2=GmMr2
これは万有引力の法則である。
(注)この式は万有引力の法則の式と同じだが、質量mは、慣性質量であり、対称性から太陽の質量Mも慣性質量と考えられる。
しかしニュートンは重力を生む質量は、慣性質量と完全には一致しない可能性もあると考え、重力質量という概念を生みだしと思われる。
既述のように、多くの実験の結果、両質量は同一であると考えられている。
重量質量を使わず、慣性質量だけを用いても、ニュートン力学を構成することが出来る。これを提唱する物理学者もいる。
それには万有引力の法則のかわりに、次の法則を採用すればよい。
外力が働かないときは、どんな2質点も、お互いに相手に向かって,
加速度運動して近ずく。両者の加速度は、両者の距離の2乗r2に反比例し、それぞれの慣性質量の比に反比例する。
式で書くと、
質点1の慣性質量と加速度の大きさをm1,α1
質点2の慣性質量と加速度の大きさをm2,α2
とすると、m1α1=m2α2、m1∝1/r2,m2∝1/r2
この法則と運動法則により2質点間に働く力(万有引力)を求めると、
ニュートンの万有引力の法則と同じ式だが、質量は慣性質量になり、
重量質量を用いずニュートン力学が構成できる。
振り子と単振動
- ウィキペディア(単振動)の「振り子」の項を見てください。
質点のつり合い
質点に力F1,,Fnが作用し、質点が静止したまま(あるいは等速直線運動)であるとき、それらの力は釣り合っているという。
釣り合いの条件は、F1+ +Fn=0です(運動の第2法則と力の合成則から導出できる)。
仕事とエネルギー
仕事
物体に力を加えて動かす時、力はこの物体に仕事をするという。仕事(の量)は力の大きさと動かした距離の積に比例する。
正確には、加えられる力→F が一定で、力の向きに対して角度θ[rad] だけ傾いている直線上を→s移動したとき、
仕事W は、
W=‖→F‖‖→s‖cosθ
で定義する。
ここで任意のヴェクトル\vec aに対して、\|\vec{a}\|は、そのベクトルの大きさ\sqrt{\sum_{i}a_i^2}である。
特に、この式において\theta=0(すなわち \cos\theta = 1)とすると
「加えられる力が一定であり力の方向が運動の方向と一致している場合」になり、
W=\|\vec F\|\|\vec s\| である。
また、\theta=\pi/2(\cos\theta = 0)のとき、W = 0となる。
すなわち、力が運動の方向と直角方向にはたらいている場合、その力は仕事をしない。
W=\|\vec F\|(\|\vec s\| \cos\theta)と表現すると、
仕事は、力の方向に\|\vec s\| \cos\thetaだけ動かしたときの仕事に等しいことが分かる。
W=(\|\vec F\| \cos\theta)\|\vec s\|と表現すると、
仕事は、
大きさ\|\vec F\| \cos\theta の\vec s方向の力を加えて、\vec sだけ動かしたときの仕事に等しい
ことが分かる。
- ウィキペディア(仕事)を参照のこと。
仕事の内積を用いた表現
内積は、仕事の記述や計算に便利な数学の概念である。
内積の定義と仕事の内積表現
ベクトル\vec a,\vec bの内積 \vec a \cdot \vec b は、\|\vec{a}\|\|\vec{b}\|
\cos\thetaで定義する。
ここで、\thetaは、ベクトル\vec a,\vec bのなす角(0\le \theta \le \pi )である。
- ウィキブックス(高等学校数学B ベクトル) の1.1.6~ 1.1.8を参照のこと。
ウィキブックスでは2次元のベクトルを中心にして説明しているが、
3次元ベクトルの場合にも、成り立つように修正することは容易である。
例えば、ベクトル\vec a = (a _1,a _2,a_3)の長さは、\|\vec a\|= \sqrt {a _1^2 +a _2^2+a _3^2},
ベクトルの内積は、この長さを使えば、全く同じ式で良い。
内積を使った 仕事の表現
内積 \cdot を用いると、
物体に力\vec{F}を加えて、\vec{PQ}(P点からQ点まで)動かした時の力のなす仕事は、
W=\vec{F}\cdot\vec{PQ} と表せる。
内積の性質
仕事は、前述のように内積で表現できるので、内積の性質を調べておくと、仕事について考察する時に役に立つ。
\vec a,\vec b,\vec cが、すべて同じ次元(2か3)のベクトルとし、 \alphaは実数とする。
(1)\vec a \cdot \vec b =\vec b \cdot \vec a
(2)\vec a \cdot \vec b =\sum_{i}a_ib_i 、
ここでa_1,b_1はそれぞれ\vec a,\vec bのx座標成分、同様に、添え字2はy座標成分、3はz座標成分
直交座標系はどんなものでも良い。しかしすべてのベクトルは同じ座標系で座標成分表示しなければならない。
(3)(\vec a +\vec b) \cdot \vec c =\vec a \cdot \vec c+\vec b \cdot \vec c
(4)(\alpha \vec a)\cdot \vec b =\vec a \cdot (\alpha \vec b)=\alpha (\vec a \cdot \vec b)
が成り立つ。
(証明)
(1)は、内積の定義から明らか。
(2);次の三角形の余弦定理を利用する。
三角形の第2余弦定理;
図のような\triangle {ABC}を考える。
頂点A,B,Cの対辺の長さをそれぞれa,b,cとし、\angle{ACB}=\thetaとする。
すると、c^2=a^2+b^2-2ab\cos\theta
余弦定理の証明;頂点Aから対辺BCにおろした垂線の足をHとする。
ピタゴラスの定理により、
c^2=\overline{BH}^2+\overline{AH}^2。\qquad 右辺の第2項に、再び、ピタゴラスの定理を適用して、
=\overline{BH}^2+(b^2-\overline{CH}^2) \qquad \overline{BH}=a-\overline{CH}を代入すると、
=(a-\overline{CH})^2+(b^2-\overline{CH}^2)=a^2+b^2-2a\overline{CH},\quad \overline{CH}=b\cos\thetaなので、代入すると
=a^2+b^2-2ab\cos\theta
証明終わり。
(2)の証明
ベクトル\vec a と\vec b を、
始点が点Cである有向線分で表現し、その終点をB,Cで表す。
すると\vec a=\vec{CB}, \vec b=\vec{CA}である。
ベクトル\vec c=\vec a-\vec bを導入すると、
\vec c=\vec a-\vec b=\vec{CB}-\vec{CA}=\vec{CB}+\vec{AC}=\vec{AB}
3角形\triangle {ABC}を考え、第2余弦定理を適用しよう。
\angle{ACB}=\thetaとおく。すると、
\|\vec c\|^2=\|\vec a\|^2+\|\vec b\|^2-2\|\vec a\|\|\vec b\|\cos{\theta}
=\|\vec a\|^2+\|\vec b\|^2-2\vec a \cdot \vec bが得られる。
この式を変形して\vec a \cdot \vec bだけを左辺に置くと、
\vec a \cdot \vec b=(\|\vec a\|^2+\|\vec b\|^2-\|\vec c\|^2)/2 。
\vec c=\vec{AB}=\vec{AC}+\vec{CB}=-\vec b+\vec aなので、
\vec a \cdot \vec b=(\|\vec a\|^2+\|\vec b\|^2-\|\vec a-\vec b\|^2)/2
この右辺を、ベクトルの直交座標成分で表すと、次式が得られる。
\vec a \cdot \vec b=(\sum_{i}a_i^2+\sum_{i}b_i^2-\sum_{i}(a_i-b_i)^2 )/2
=\sum_{i}a_i b_i
(2)の証明終わり。
(性質3)の証明;ある一つの直交座標系をさだめ、両辺を、性質(2)を利用して、座標成分であらわす。両辺が等しいことが分かる。
(性質4)の証明;同様に、3つの式を、座標成分表示すれば、みな等しいことが、簡単に分かる。
物体が曲線運動するときの仕事量の求め方
力を受けた時の物体の運動は直線とは限らないが、運動の軌跡を細かく区切って眺めると、線分に近いので、物体の変位は、ごく短い線分をつなぎ合わせたものと考える。すると各線分毎に仕事を計算しそれをたせば、全体の仕事量を求めることができる。
エネルギー
物質の持っている仕事をする能力をエネルギーという。
- エネルギー(ウィキペディア)の自然科学の項を参照のこと。
仕事の単位
仕事の定義W=\|\vec F\|\|\vec s\| \cos\thetaから、仕事の単位は、力の大きさ\|\vec F\|の単位と長さ\|\vec s\|の単位を掛けたものになる( \cos\theta は無単位なので )。
MKSA単位系では、力の大きさの単位はN(ニュートン)、長さの単位はm(メートル)なので、仕事の単位はNm となる。
これをJ(ジュール)と呼ぶ。J=Nmである。
質点系の運動
2個以上の質点が集まって出来ている系を質点系という。
質点系というときは、各質点は密集していても、離れ離れでも良い。互いに固着しようが、自由に動けようが構わない。
すべての物質は、分子の集合と考えたり、細分化して極小部分に分け、それらの集合と考えれば、十分な精度で、質点系とみなすことができる。
そのため質点系の運動の法則を、ニュートンの運動法則から導出すれば、その応用範囲は非常に広い。
質点系の運動と重心
系の任意の2つの質点間には作用・反作用の法則を満たす力が働いていてもよい。
この力を質点系の”内力”という。
質点系の各質点に外部から力(外力という)が加わる時、この質点系はどんな運動をするだろうか。
質点系の各質点の位置を\vec{r_i}、質量をm_i とし、
質点m_i に作用する外力を\vec{f_i}、
m_i に、他の質点m_j から作用する内力を\vec{f_{ij}}とする(i,j=1 \ldots N)。
すると、各質点に対して、運動の第2法則により、
d (m_i \vec{v_i})/dt=\vec{f_i}+\sum_{j\neq i}\vec{f_{ij}} \qquad ここで\vec{v_i}=d\vec{r_i}/dt、
各ベクトルを自由ベクトルとみなしてi=1 \ldots Nについて加え合わせると、\vec{f_{ij}}+\vec{f_{ji}}=0なので、
\frac{d^2}{dt^2} \sum_i{ m_i \vec{r_i}} =\frac{d}{dt} \sum_i{ m_i \vec{v_i}} =\sum_i{\vec{f_i}}
が得られる。
質点系の全質量M= \sum_i{m_i} と質点系に働く全外力\vec{F}= \sum_i{\vec{f_i}} を用いて書きなおすと、
M\frac{d^2}{dt^2}(\sum_i{ m_i \vec{r_i}}/M)= \vec{F}
質点系の重心\vec{R}を \quad \vec{R}=\sum_i{ m_i \vec{r_i}}/M で定義すると、
M\frac{d^2}{dt^2}\vec R= \vec{F}
この式は、力\vec{F}をうける質量Mの質点の運動方程式と同じである。
以下の解説も参考にしてください。
複雑にみえる運動も重心の運動をみれば簡単である
体操選手の運動は、跳躍や着地などで空中をまいながら、回転や体の屈伸、ひねりなどを行う。大変複雑である。
しかし、導出した質点系の重心の運動法則から、体の重心の運動は、投射体の運動であり、放物線をえがいて移動することが分かる。
空中に飛び出た瞬間の速度(速さと方向・向き)で、その軌跡は完全に決まってしまうのである。
剛体の運動とつり合い
剛体
剛体(Rigid body)とは、
質点系であって、それらの、どの2質点の間の距離も変わらない,特殊な系のことを言う。
どの2質点の間の距離も変わらなければ変形は起こらない。
固くて変形しにくい物体を理想化した概念である。
剛体の運動
剛体は変形しない質点系なので、その運動は、重心の運動と、重心の周りの回転運動を合成したものになる。
重心の運動は前の節で説明したように、質点の運動と同じように簡単に扱える。
重心の周りの回転運動について解析するには、少し難しい数学が必要になる。
- ウィキペディア(剛体の力学)を参照のこと。
このテキストでは、固定軸の周りの回転運動を中心に、 剛体運動の初歩と釣合の条件について学ぶ。
固定軸のまわりの回転運動
剛体が、剛体の中を通る固定軸の周りを回転する運動(車輪の回転など)を考える。
応用も考え、回転軸は重心を通らなくてもよいように一般化しておく。
(注)なお、軸が動かないようにするためには軸受が必要である。
工夫しても回転時に軸は軸受から多少の摩擦力を受け、回転にブレーキがかかる。
しかし、これは無視出来るほど小さいと仮定する。
すると軸が受ける力は、軸の変動を防ぎ、固定軸の周りの運動に限定させる作用を持ち、
回転を遅める作用は持たないことになる。
回転運動の表示法
固定軸まわりの剛体の運動はどのように表示したらよいだろうか。
・剛体の位置を表す変数;回転角
剛体が幾ら回転したか分かるように、剛体の、回転軸上にない一点P_sに印を付ける。
次に、角度を測る基準線をきめるため、座標系を決めよう。
P_sから固定軸へ垂線をひき、その足を原点Oとし,固定軸をz座標とする(静止した)3次元直交座標O-xyzを考える。
剛体が固定軸の周りを回転すると、印P_sはxy平面上を、原点Oを中心に円を描いて動くことになる。
その位置ベクトル\vec{OP_s}がx軸の正方向となす角度\phiを、回転角と呼ぶ。図参照。
但し、x軸から反時計回りの角を正にする。
また一回転した後ならば、一回転の角2\piを加え、逆周りに一回転した後なら2\piを引き、
角度だけでなく回転数も分かるようにする。
回転角が指定されると、点P_sの位置が決まる。
それだけでなく剛体は変形しないので、剛体のすべての点の位置がきまる。
そこで回転角\phiの時間変化\phi=\phi (t)を明らかにすれば、剛体の回転運動は定まる。
固定軸のまわりの回転運動において回転角の果たし役割は、質点の運動において質点の位置が果たし役割に対応していることが分かる。
・回転の角速度と角加速度
\phi=\phi (t)を時間で微分したd\phi (t)/dtを回転の角速度と呼ぶ。
直観的には、時刻tの瞬間の、回転の速さ(回転角の時間に対する変化率)を表す。
さらにもう一回時間微分したd^2\phi (t)/dt^2を回転の角加速度と呼ぶ。
回転力(トルク)
質点の運動に倣って、剛体に作用する力によって、その位置(=回転角)がどう変化するかの法則を導出したい。
しかし、剛体の回転の場合、ある方向の力は、剛体の回転に全く関係しない。
例えば、回転軸から放射状にでる半直線方向の力は全く回転の変化に寄与しない。
そこで剛体の回転を変化させる力とはなにかという問題から考察する必要が起こる。
質点運動における力の定義(力と運動量の変化の関係)や力と仕事の関係など力の係っている式のなかから、
剛体の回転運動に容易に拡張出来るものを選び、その式から、回転に関する力を求めることを試みる。
力の定義からは、回転運動への拡張を、推測することは難しい。
力と仕事の関係の考察をしてみよう。
力と仕事の関係からの考察
適当な直交座標系をさだめ、ベクトルは、座標成分で表示する。
質点に、一定の力\vec F=(F_x,F_y,F_z)を作用させて、x軸方向に変位させる。
質点はこの軸の上でしか動けないように拘束され、摩擦はないと仮定する。
質点の変位ベクトルは一次元の変数xを使って\vec s=(x,0,0)と表せる。
すると力のなす仕事は、W=\vec F \cdot (x,0,0)=F_{x}xである。
逆に物体に一定の力を加え、x軸上でxだけ変位させた時の仕事Wが分かれば、質点を動かした力は
F_x=W/x
で求められる。
F_y,F_zは、質点をx軸上で動かすことには全く寄与せず、
x軸に拘束された質点を動かす力は、F_xなのである。
固定軸まわりの回転もその変位は一次元の変数である回転角度で表わせるので、
これに倣って、
W/回転した角度
を、回転にかんする力であると考える。これを回転力と呼ぶ。トルクともいう。
この方針を実行して回転力を具体的に求めよう。
剛体に力を加え微小角動かす時の、力のなす仕事の算出
図4.1のように剛体の任意の一点P(x,y,z)を考える。
z座標の上方からxy平面を見下ろしているので、z座標は点になりOと書いてある。
まず一点P(x,y,z)に力\vec F=(F_{x},F_{y},F_{z})が作用して、微小角\Delta\thetaだけ回転したときの
仕事\Delta Wを計算し回転力を求めよう。
P点から回転軸(z軸)に垂線を下ろし、その足をO'=(0,0,z)とする。
\vec{O'P}の長さをr、x軸となす角を\theta(ラジアン)と置く。
この角度は、
剛体につけた印の位置ベクトル\vec{OP_s}がx軸となす回転角\phiと
このベクトルと\vec{O'P}(をxy平面に平行移動したベクトル)の間の角の和である。
後者は、剛体なので、運動しても変わらない定数である。そこで、\theta=\phi+定数,と書ける。
剛体がz軸の周りを微小角\Delta\theta回転して、点Pが図の点Qに移動したとする。
すると角\angle OPQはほぼ直角(=\pi /2)で\vec{PQ}の長さPQは、PQ=r(\Delta\theta)。
\vec{PQ}のx成分とy成分は、図4-1中に示したように、それぞれ、-QR=-PQ*y/r、PR=PQ*x/r。
PQ=r(\Delta\theta)を代入すると、
\vec{PQ}_x=-y(\Delta\theta)、\vec{PQ}_y=x(\Delta\theta)、\vec{PQ}_z=0
点P(x,y,z)に作用する力\vec{F}=(F_{x},F_{y},F_{z})が、物体を\vec{PQ}だけ動かしたので、
その仕事は、\Delta W=\vec{F} \cdot \vec{PQ}(内積)。
この右辺を内積の性質を用いて座標成分で表すと、
F_{x}*(-y)\Delta\theta+F_{y} x\Delta\theta+F_{z}* 0
=(xF_{y}-yF_{x})*\Delta\theta
z軸まわりの回転力の導出
ゆえに、力\vec{F}のz軸まわりの回転力(トルク)T_\vec{e_z}は\Delta W/\Delta\theta=xF_{y}-yF_{x}
z軸まわりの回転力(トルク)の性質
(1)力\vec{F}のz軸まわりの回転力は,\vec{F}_zには関係しない。
言いかえるとz軸を固定軸とする剛体にz軸の方向の力を加えても、z軸の周りの回転は起こらない。
(2)剛体の1点P(x,y,z)に作用する力\vec Fを考える。
点P(x,y,z)からz軸に下ろした垂線の足をO'(0,0,z)と書く。
力\vec Fを、,
\vec{O'P}方向の成分\vec F_rと、
z軸まわりの回転によりPの描く、O'を中心とする回転円の(左回りの)接線方向の成分\vec F_t
および、これら2成分に直交する成分(z軸と平行)
に分解する(図参照)。この時、
・力\vec F_rのz軸まわりの回転力は、零である。
すなわち、動径方向の力は回転に寄与しない。
・力\vec Fのz軸まわりの回転力は、\vec F_tのz軸まわりの回転力に等しい。
数式で表すと、xF_{y}-yF_{x}=x(F_t)_{y}-y(F_t)_{x}
(3)剛体に作用する力の作用点を、力の作用線上で動かす限り、回転力は変化しない。
ここで、力の作用線とは、力の作用点を通り、力の方向と重なる直線のこと。
これらはいずれも直観と合致する。
証明は、試みてほしい。
他の軸の周りの回転力
力\vec{F}のx軸、y軸まわりの回転力も同様に計算できる。結果は、
x軸まわりの回転力;yF_{z}-zF_{y}=y(F_t)_{z}-z(F_t)_{y}
y軸まわりの回転力;zF_{x}-xF_{z}=z(F_t)_{x}-x(F_t)_{z}
原点まわりの力のモーメント
位置ベクトル\vec r=(x,y,z)の剛体の点Pに作用する力\vec Fの原点まわりの力のモーメントを、
\vec N=(x軸まわりのトルク、y軸まわりのトルク、z軸まわりのトルク)で定義する。
数式で書くと、
\vec N=(yF_{z}-zF_{y},zF_{x}-xF_{z},xF_{y}-yF_{x}),
ベクトル積と力のモーメントのベクトル積表示
以上の結果は、ベクトル積(クロス積ともいう)を用いると簡潔、正確に表現でき、
回転運動の性質を調べるのが容易になる。
3次元ベクトル\vec a,\vec b のベクトル積\vec a \times \vec bとは、3次元ベクトルであり,
大きさは\vec a,\vec b を2辺とする平行四辺形の面積に等しく、
方向はこの四辺形に垂直で、向きは、(\vec a,\vec b,\vec a \times \vec b)が右手系をなすように定めたものである。
次の項で説明するベクトル積の性質6を用いると、
位置ベクトル\vec rの点に作用する\vec F の
原点まわりの力のモーメントは、\vec N = \vec r \times \vec F
x軸まわりの回転力(トルク)は、\vec N \cdot \vec e_x と表せることが分かる。
y軸とz軸周りの回転力も、それぞれ
\vec N \cdot \vec e_y ,\quad \vec N \cdot \vec e_z で
表せる。
ベクトル積の性質
力のモーメントやトルクの性質を調べるには、ベクトル積の性質についての知識が必要になる。
\vec{a}, \vec{b}, \vec{c}を3次元ベクトル
\alphaを実数とする。
すると次の性質が成り立つ。
性質1. \quad \vec{a} を, \vec{c} と垂直な成分 \vec{a_\perp} と,平行な成分\vec{a_\parallel} の和に分解するとき、
\quad \vec{a} \times \vec{c}= \vec{a_\perp} \times \vec{c}
\quad \vec{a_\parallel} \times \vec{c}= 0
性質2. \quad \vec{a} \times \vec{b}= -\vec{b} \times \vec{a}
性質3. \quad (\vec{a}+ \vec{b})\times \vec{c}= \vec{a} \times \vec{c} + \vec{b} \times \vec{c}
性質3の系. \quad \vec{a} \times (\vec{b}+ \vec{c})= \vec{a} \times \vec{b} + \vec{a} \times \vec{c}
\quad (\vec{a}+ \vec{b}+\vec{c})\times \vec{d}=\vec{a}\times \vec{d}+\vec{b}\times \vec{d}+\vec{c}\times \vec{d}
性質4. \quad (\alpha\vec{a})\times \vec{b}= \alpha(\vec{a} \times \vec{b})= \vec{a}\times (\alpha\vec{b})
性質5.\quad (\vec{e_1},\vec{e_2}, \vec{e_3}) を
それぞれ大きさ(長さ)1で互いに直交し、右手系をなす、ベクトル(右手系をなす正規直交基底)とする。
この時、
\quad \vec{e_1} \times \vec{e_2} = \vec{e_3}, \quad
\vec{e_2} \times \vec{e_3} = \vec{e_1}, \quad
\vec{e_3} \times \vec{e_1} = \vec{e_2}
性質6.ベクトル\vec a, \vec bを,性質5で用いた基底 (\vec{e_1},\vec{e_2}, \vec{e_3}) で決まる座標の座標成分で表示しておく。
すると\vec a \times \vec b=(a_yb_z-a_zb_y,a_zb_x-a_xb_z,a_xb_y-a_yb_x)
性質7. \quad (\vec{a} \times \vec{b})\cdot \vec{c}= (\vec{c} \times \vec{a})\cdot\vec{b} =(\vec{b} \times \vec{c})\cdot\vec{a}
性質8. \quad \vec{a(t)} と \vec{b(t)} を,tにかんして微分可能な、ベクトルに値をとる関数とする。すると、
\quad \vec{a(t)} \times \vec{b(t)} は、tにかんして微分可能で、
\quad \frac{d}{dt}(\vec{a(t)} \times \vec{b(t)}) =(\frac{d}{dt}\vec{a(t)} )\times \vec{b(t)}+\vec{a(t)}\times (\frac{d}{dt}\vec{b(t)})
証明
性質1の証明;ベクトル積の定義から、容易に示せる。
2つのベクトルの作る平行四辺形の面積と方向・向きを考えれば良い。
性質2の証明;2つのベクトルを入れ替えても、それらが作る平行四辺形の面積は変わらず、この四辺形に直交する直線の方向も変わらない。しかし、右手系をなす方向は、逆向きになる。ベクトル積の定義から、 \qquad \vec{a} \times \vec{b}= -\vec{b} \times \vec{a} が示せた。
性質3の証明;
この証明には少し工夫が必要である。
ベクトル積の性質の中でも、もっとも大切なものなので、詳しく説明しよう。
① \vec{a}, \quad \vec{b} と \vec{c} が直交する場合。図参照のこと
・議論をやさしくするため、ベクトルを、空間の原点O を始点とする有向線分で代表させる。
・ \vec{c} と直交しO を通る平面をHとする。
・仮定より \vec{a},\quad \vec{b}は、ともに平面H上のベクトルである。
・\vec{a} \times \vec{c} ,\quad \vec{b} \times \vec{c}も、
ベクトル積の定義により、共に \vec{c} と直交するので、H上のベクトルである。
これら四つのベクトルはすべて平面H上にあるので、今後の議論はこの平面上で進める。
ⅰ)\vec{a} \times \vec{c}, \vec{b} \times \vec{c} の張る平行四辺形は,
\vec{a}, \vec{b}の張る平行四辺形を、\| \vec{c}\|倍し,原点周りに90度回転したものになることを、示そう。
・\vec{a} \times \vec{c} は、ベクトル積の定義から、 \vec{a} と直交する。
そのため、\vec{a} を平面H上で、原点まわりに、90度右回りか、左回りすれば、方向と向きが一致する。
・\vec{b} \times \vec{c} も、同様に考え、\vec{b} を平面H上で、原点まわりに、90度右回りか、左回りすれば、方向と向きが一致することが分かる。
・どちら周りの回転になるかは、ベクトル積の定義によって決まるが、
後者の回転の向きが、前者の回転の向きと一致することが分かる。
・\vec{a}\times \vec{c} の大きさは、
\|\vec{a}\times \vec{c}\|=\|\vec{a}\|\|\vec{c}\|\cos\pi/2=\|\vec{a}\|\|\vec{c}\| なので、\vec{a} の大きさの\|\vec{c}\|倍になる。
同様に、\vec{b}\times \vec{c} の大きさは、\vec{a} の大きさの\|\vec{c}\|倍になる。
・以上の結果より、所望の結果は示された。
ⅱ) \qquad (\vec{a}+ \vec{b})\times \vec{c}= \vec{a} \times \vec{c} + \vec{b} \times \vec{c}を示そう。
・ ⅰ)と同じ議論により、
(\vec{a}+ \vec{b}) \times \vec{c}は\vec{a}, \vec{b}の張る平行四辺形の対角線を、原点周りに90度、同じ向きに回転させ、\|\vec{c}\|倍させたものであることが分かる。
・すると、ⅰ)で示したことから、(\vec{a}+ \vec{b}) \times \vec{c}は
\vec{a} \times \vec{c}, \vec{b} \times \vec{c} の張る平行四辺形の対角線\vec{a} \times \vec{c}+\vec{b} \times \vec{c} に等しいことが分かる。
・以上で①が示せた。
② 一般の場合。
性質1より、\perp を\vec{c}と垂直な成分を表すとすると、 (\vec{a}+ \vec{b})\times \vec{c}= (\vec{a}+ \vec{b})_\perp \times \vec{c} \qquad \qquad \qquad (1)
(\vec{a}+ \vec{b})_\perp =\vec{a}_\perp +\vec{b}_\perpなので、(1)式は、
= (\vec{a}_\perp +\vec{b}_\perp) \times \vec{c}
①より、
= \vec{a}_\perp \times \vec{c}+\vec{b}_\perp\times \vec{c}=\vec{a} \times \vec{c}+\vec{b} \vec{c} \qquad 性質3の証明終わり。
性質3の系の証明;
性質2より、
\vec{a} \times (\vec{b}+ \vec{c})= -(\vec{b}+ \vec{c})\times \vec{a}
性質3より、
= -(\vec{b} \times \vec{a}+ \vec{c} \times \vec{a})
再び性質2より、
=\vec{a} \times \vec{b} + \vec{a} \times \vec{c} \quad 前半の証明終わり
性質2より、
\quad (\vec{a}+ \vec{b}+\vec{c})\times \vec{d}=(\vec{a}+ \vec{b})\times \vec{d}+\vec{c})\times \vec{d}
再び性質2より、
=\vec{a}\times \vec{d}+\vec{b}\times \vec{d}+\vec{c}\times \vec{d}
\quad証明終わり。
性質4の証明;実数\alpha が正、零、負の場合に分けて考える。いずれの場合にも
ベクトル積の定義とベクトルと実数の積の性質から、容易に証明できる。
性質5の照明;ベクトル積と(e_1,e_2,e_3) の定義から明らかである。
性質6の証明;\vec a=a_x\vec{e_x}+a_y\vec{e_y}+a_z\vec{e_z},
\vec b=b_x\vec{e_x}+b_y\vec{e_y}+b_z\vec{e_z}と表せるので、
\vec a \times \vec b=(a_x\vec{e_x}+a_y\vec{e_y}+a_z\vec{e_z})\times \vec b
性質3の系から
=a_x\vec{e_x}\times \vec b
+a_y\vec{e_y}\times \vec b
+a_z\vec{e_z}\times \vec b \qquad (1)
式(1)の第1項
a_x\vec{e_x}\times \vec b
に
\vec b=b_x\vec{e_x}+b_y\vec{e_y}+b_z\vec{e_z}
を代入して、性質3の系を使って変形すると、
a_x\vec{e_x}\times \vec b
=a_x\vec{e_x}\times b_x\vec{e_x}
+a_x\vec{e_x}\times b_y\vec{e_y}
+a_x\vec{e_x}\times b_z\vec{e_z} \qquad (2)
性質4と性質5を使うと、
a_x\vec{e_x}\times b_x\vec{e_x}
=a_x b_x\vec{e_x}\times \vec{e_x}
=\vec 0 。
同様の計算を行うと、
a_x\vec{e_x}\times b_y\vec{e_y}
=a_x b_y\vec{e_x}\times \vec{e_y}
=a_x b_y\vec{e_z}
a_x\vec{e_x}\times b_z\vec{e_z}
=a_x b_z\vec{e_x}\times \vec{e_z}
=-a_x b_z\vec{e_y}
式(2)にこれらを代入して、
a_x\vec{e_x}\times \vec b
=a_x b_y\vec{e_z} - a_x b_z\vec{e_y} \qquad (3)
式(1)の第2項、第3項も同様に計算すると、
a_y\vec{e_y}\times \vec b
=a_y b_z\vec{e_x} - a_y b_x\vec{e_z} \qquad (4)
a_z\vec{e_z}\times \vec b
=a_z b_x\vec{e_y} - a_z b_y\vec{e_x} \qquad (5)
式(3),(4),(5) を、式 (1)に代入すると、
\vec a \times \vec b
=a_x b_y\vec{e_z} - a_x b_z\vec{e_y}
+a_y b_z\vec{e_x} - a_y b_x\vec{e_z}
+a_z b_x\vec{e_y} - a_z b_y\vec{e_x}
=(a_y b_z - a_z b_y)\vec{e_x}
+(a_z b_x - a_x b_z)\vec{e_y}
+(a_x b_y - a_y b_x)\vec{e_z}
性質6の証明終わり。
性質7の証明;
\quad (\vec{a} \times \vec{b})\cdot \vec{c}= (\vec{c} \times \vec{a})\cdot\vec{b}を証明しよう。
残りも、同様に証明出来るので各自試みてください。
右手系をなす一つの直交座標を決める。
3つのベクトルを、この座標の成分で表示して、性質6と内積の性質を使えば、左右が等しいことが証明できる。
概略をスケッチしよう。
\quad (\vec{a} \times \vec{b})\cdot \vec{c}
=(a_yb_z-a_zb_y,a_zb_x-a_xb_z,a_xb_y-a_yb_x)
\cdot (c_x,c_y,c_z)
=(a_yb_z-a_zb_y)c_x+(a_zb_x-a_xb_z)c_y+(a_xb_y-a_yb_x)c_z
\quad (\vec{c} \times \vec{a})\cdot\vec{b}も、これと同じように計算する。
これら両式を整頓すると、同じものであることが分かる。
性質7の証明終わり。
性質8の証明;
性質8. \quad \vec{a(t)} と \vec{b(t)} を,tにかんして微分可能な、ベクトルに値をとる関数とする。すると、
\quad \vec{a(t)} \times \vec{b(t)} は、tにかんして微分可能で、
\quad \frac{d}{dt}(\vec{a(t)} \times \vec{b(t)})
\quad =(\frac{d}{dt}\vec{a(t)} )\times \vec{b(t)}+\vec{a(t)}\times (\frac{d}{dt}\vec{b(t)})
すでにこのテキストで紹介した、ベクトル値関数の微分の定義を用いて証明する。
\quad \frac{d}{dt}(\vec{a(t)} \times \vec{b(t)})
=\lim_{\delta t \to 0}
(\vec a(t+\delta t)\times \vec{b(t+\delta t)}- \vec{a(t)} \times \vec{b(t)})/\delta t \qquad (1)
この極限が存在し、
\frac{d}{dt}\vec{a(t)} \times \vec{b(t)}+\vec{a(t)}\times \frac{d}{dt}\vec{b(t)}
になることを示せば性質8は証明できたことになる。
極限の計算が進むよう、右辺の式の分母は変形しよう。
関数の積の微分公式の証明と同じ技巧を用いる。
\vec a(t+\delta t)\times \vec{b(t+\delta t)}
- \vec{a(t)} \times \vec{b(t)}
= \vec a(t+\delta t)\times \vec{b(t+\delta t)}
-\vec a(t)\times \vec{b(t+\delta t)}
+\vec a(t)\times \vec{b(t+\delta t)}
- \vec{a(t)} \times \vec{b(t)}
ベクトル積の性質3を利用すると、
= \left(\vec a\left(t+\delta t\right) -\vec a\left(t\right)\right)
\times
\vec b\left(t+\delta t\right)
+\vec a\left(t\right)\times \left(\vec b\left(t+\delta t\right)- \vec b\left(t\right)\right)
この式を式(1)の右辺の分子の項に代入し整頓すると
\quad \frac{d}{dt}(\vec{a(t)} \times \vec{b(t)})
=\lim_{\delta t \to 0}
\frac{\vec a(t+\delta t)\times \vec b(t+\delta t)- \vec a(t) \times \vec b(t)}
{\delta t}
=\lim_{\delta t \to 0}
\frac{\left(\vec a\left(t+\delta t\right) -\vec a\left(t\right)\right)
\times
\vec b\left(t+\delta t\right)
+\vec a\left(t\right)\times \left(\vec b\left(t+\delta t\right)- \vec b\left(t\right)\right) }
{\delta t}
ベクトル積の性質4を使い、
=\lim_{\delta t \to 0}\left(
\frac{\vec a(t+\delta t) -\vec a(t)}{\delta t}
\times
\vec b\left(t+\delta t\right)
+
\vec a(t)\times \frac{\vec b(t+\delta t)- \vec b(t)}
{\delta t}
\right)
極限の性質を使って、
=\lim_{\delta t \to 0}
\frac{\vec a(t+\delta t) -\vec a(t)}{\delta t}
\times
\lim_{\delta t \to 0}\vec b(t+\delta t)
+
\vec a(t)\times
\lim_{\delta t \to 0}\frac{\vec b(t+\delta t)- \vec b(t)}{\delta t}
式中の極限は、\vec a,\vec bが、微分可能なので存在し、
\lim_{\delta t \to 0} \frac{\vec a(t+\delta t) -\vec a(t)}{\delta t}
=\frac{d\vec a(t)}{dt}
\lim_{\delta t \to 0} \frac{\vec b(t+\delta t) -\vec b(t)}{\delta t}
=\frac{d\vec b(t)}{dt}
また、\lim_{\delta t \to 0}\vec b(t+\delta t)=\vec b(t) なので、
所望の結果が得られた。性質8の証明終わり。
力のモーメントの性質
ベクトル積の性質が分かったところで、再び、力のモーメントの考察に戻る。
剛体の一点 P に加えられた力 \vec F の、原点周りの力のモーメントは、
\vec N= \vec{OP} \times \vec F で定義した。
すると、
x軸まわりの回転力(トルク)は、\vec N \cdot \vec e_x 、
y軸とz軸周りの回転力も、それぞれ
\vec N \cdot \vec e_y ,\vec N \cdot \vec e_z
で表せることは、すでに説明した。
ところが、もっと一般に、どんな軸の周りの回転力も、\vec N から得られる。
定理;
\vec eを、原点を始点とする大きさ1の任意のベクトルとする。
すると、
\quad \vec N \cdot \vec eは、力\vec{F}の\vec e軸の周りの回転力になる。
証明;
9つに分けて示す。
ⅰ)準備
図のように、剛体の点 P から、\vec e 軸に垂線を下ろし、その足を Q とする。
力 \vec F のもとで、剛体が \vec e を固定軸にして、
微小時間に、微小角\delta \phi だけ回転したとする。
このとき、P が移った先を、P' とする。
ⅱ)回転角 \delta \phi が微小なので、
この回転中の P の軌跡(円弧の微小部分)は、有向線分\vec{PP'} で精度高く、近似できる。
ⅲ)この間に力 \vec F がなした仕事 \delta W は、\delta W=\vec{PP'} \cdot \vec F
この仕事を、回転角\delta \phiで割ると、力の \vec e 軸周りの回転力が得られる。そこで、\vec{PP'} を、この定理で与えられている諸量を使って表現し、これを用いて、仕事を計算しよう。
ⅳ)有向線分\vec{PP'}の方向を求める。
\vec{PP'} は、\vec e 軸と垂直でQ を通る平面H上にあり、
Qを中心とする円の弧の微小部分をなすので、線分QP と直交する。\vec{PP'}\perp QP
また、\vec e 軸と垂直でQ を通る平面H上にあるので、
\vec{PP'}は\vec e 軸とも直交し、従って線分OQと直交する。\vec{PP'}\perp OQ
ゆえに、\vec{PP'} は、3点O,Q,Pを通る平面 OQP と直交する。
すると、\vec{PP'} は、平面 OQP 上のすべての線分と直交する。
ゆえに、\vec{PP'}\perp \vec e,\quad \vec{PP'}\perp \vec{OP}
これで、\vec{PP'} の方向は、求まった。
ⅴ)有向線分\vec{PP'} の向き
点 P は、\vec e 軸の周りを右周りに回転するので、その向きは、
\vec e \times \vec{OP} と同じ向きである。
ⅵ)\vec{PP'} の大きさ。
\vec{PP'}は、 Q を中心とする、半径 \| \vec{QP} \| の円弧の一部なので、
その中心角\delta \phi を用いて、\| \vec{PP'}\|=\|\vec{QP}\|\delta \phi
ⅶ)ⅳ)、ⅴ)、ⅵ)から
\vec{PP'}=\frac {\vec e \times \vec r}{\|\vec e \times \vec r \|}\|\vec{QP}\|\delta \phi
ⅷ)\vec{PP'}=\vec e \times \vec r \delta \phiが成り立つ。
なぜなら、
\|\vec e \times \vec r \|= \|\vec e \|\|\vec r \|\sin \theta =\|\vec r \|\sin \theta =\| \vec{QP} \| ,ここで \theta は\vec e と\vec r の間の角。
この式をⅶ)で得られた式に代入すれば、所望の結果が得られる。
ⅸ)\delta W=\vec{PP'} \cdot \vec F
=(\vec e \times \vec r \delta \phi) \cdot \vec F
=(\vec e \times \vec r) \cdot \vec F \delta \phi
=(\vec r \times \vec F)\cdot \vec e \delta \phi
ⅹ) T_\vec e = \frac{\delta W}{\delta \phi} =(\vec r \times \vec F)\cdot \vec e =\vec N \cdot \vec e
定理の証明終わり。
(注)剛体が固定軸の周りでなく、自由に回転するときでも、
ある瞬間には、ある軸の周りの回転になっている。
力のモーメントは、どんな軸周りの回転力の情報も含んでいることが証明されたので、
回転運動一般に有効な概念であることが分かる。
剛体の複数個所に作用する力の回転力
次に剛体の多くの点に力を加えたときの回転力を求めよう。
力の作用点をP_i(x_i,y_i,z_i)、力を\vec{F_i}\quad (i=1,2,,,n)とする。
これらの力のもとで剛体がz軸まわりを\Delta\thetaだけ微小回転するときの、各力のなす仕事の合計は、
(\sum_{i=1}^{n}(x_{i}(F_{i})_{y}-y_{i}(F_{i})_{x}))*\Delta\theta
従って、作用点P_i(x_i,y_i,z_i)の力、\vec{F_i}\quad (i=1,2,,,n)のz軸まわりの回転力は、
\sum_{i=1}^{n}(x_{i}(F_{i})_{y}-y_{i}(F_{i})_{x})
この結果は、力が一か所に作用する場合と全く同じ考えで導ける。
同様に、x軸まわりとy軸まわりの回転力も導ける。
従って、力\vec{F_i}\quad (i=1,2,,,n)の原点まわりの力のモーメントは、
\vec N=\sum_{i=1}^{n}\vec N_i
=\sum_{i=1}^{n}(y_{i}(F_{i})_{z}-z_{i}(F_{i})_{y},z_{i}(F_{i})_{x}-x_{i}(F_{i})_{z},x_{i}(F_{i})_{y}-y_{i}(F_{i})_{x})
回転運動の方程式
\vec N が、あらゆる回転軸にかんする回転力を表現していることがわかった。
力Fと運動量の変化の関係をあたえるニュートンの運動方程式(第2法則)を変形して、回転力\vec Nにかんする方程式を導こう。
直交右手座標系O-xyz を定める。原点 O は、考察対象に都合のよい点を選ぶ。
剛体をN個の(質点と考えてよい)微小部分P^i(i=1 \cdots N)に分け、
その質量をm_i、位置ベクトルを\vec{r}^i(x_i,y_i,z_i)とする。
P_iが外部から受ける力を\vec {F}^i、
P_i が剛体の他の部分P_j(j\neq i) から受ける力(内力)を\vec {F}^{ij}とおく。
後者は、剛体が変形しないよう、剛体の原子間に働かせる力に起因する。
この原子間の力は、原子の電荷による電気力と、
原子同士が接近しすぎたときに作用する量子力学的力により生じる。
作用・反作用の法則(運動の第3法則)から、\vec F^{ij}=-\vec F^{ji} 。
さらに、剛体の2点間に働く内力の方向は、
その2点を結ぶ直線の方向と同じだと、仮定する。
各質点のニュートンの運動方程式
各質点ごとに、ニュートンの運動方程式を立てると、
m_i\frac{d^2\vec r^i}{dt^2}=\vec F^i+\sum_{j\neq i}\vec F^{i,j} \quad (i=1 \cdots N)
これを変形して
\vec F^i=m_i\frac{d^2\vec r^i}{dt^2}-\sum_{j\neq i}\vec F^{i,j} \quad (i=1 \cdots N) \qquad (1)
この式から、
力\vec F^iの回転力\vec N^i=\vec r^i \times \vec F^iにかんする式を導こう。
\vec N^i=\vec r^i \times \vec F^iにかんする式の誘導
式(1)の両辺に左側から、\vec r^i のベクトル積を施すと、
\vec N^i=\vec r^i \times \vec F^i
=\vec r^i \times (m_i\frac{d^2\vec r^i}{dt^2}-\sum_{j\neq i}\vec F^{i,j})
(i=1 \cdots N)
ベクトル積の性質3と性質4により、
=m_i\vec r^i \times \frac{d^2\vec r^i}{dt^2}-\sum_{j\neq i}\vec r^i \times\vec F^{i,j} \qquad (2)
ここで、ベクトル積の性質8より
\frac{d}{dt}(\vec r^i \times \frac{d \vec r^i}{dt})
= \frac{d \vec r^i}{dt} \times \frac{d \vec r^i}{dt}
+\vec r^i \times \frac{d^2\vec r^i}{dt^2}
=\vec r^i \times \frac{d^2\vec r^i}{dt^2}
なので、
\vec N^i=m_i\frac{d}{dt}(\vec r^i \times \frac{d \vec r^i}{dt})
-\sum_i\vec r^i\times \vec F^{i,j}
= \frac{d}{dt}(\vec r^i \times m_i\frac{d \vec r^i}{dt})
-\sum_{j\neq i}\vec r^i\times \vec F^{i,j} \qquad (3)
質点P_iの運動量を\vec P^iと書くと、
P^i=m_i\vec v^i=m_i\frac{d\vec r^i}{dt}なので、
\vec N^i=\frac{d}{dt}(\vec r^i \times \vec P^i)
-\sum_{j\neq i} \vec r^i \times \vec F^{i,j}
定義;角運動量(運動量のモーメントともいう)
質点の位置ベクトルを\vec r、運動量を\vec pと書くとき、
\vec l=\vec r \times \vec pを,この質点の角運動量と呼ぶ。
これを用いると、
\vec N^i=\frac{d\vec l^i}{dt}-\sum_{j\neq i} \vec r^i \times \vec F^{i,j}
回転の運動方程式の導出
故に、
\vec N=\sum_i\vec N^i=\frac{d\sum_i \vec l^i}{dt}-\sum_i \sum_{j\neq i} \vec r^i \times \vec F^{i,j} \qquad (4)
ここで、
\sum_i \sum_{j\neq i} \vec r^i \times \vec F^{i,j}=\sum \sum_{i<j}\vec r^i \times \vec F^{i,j}+\sum \sum_{i>j}\vec r^i \times \vec F^{i,j} \qquad (5)
式(4)の右辺の第2項の上付き添え字i,jを、それぞれ、j'と i'でおきかえられるので、
\sum \sum_{i>j}\vec r^i \times \vec F^{i,j}
=\sum \sum_{j'>i'}\vec r^{j'} \times \vec F^{j',i'}
内力は作用反作用の法則が適用できると仮定しているので、
\vec F^{j',i'}=-\vec F^{i',j'} 。この式を上の式の右辺に代入すると、
\sum \sum_{i>j}\vec r^i \times \vec F^{i,j}
=-\sum \sum_{j'>i'}\vec r^{j'} \times \vec F^{i',j'}
この式の右辺の和をとる変数i',j' を i,j におきかえると、
\sum \sum_{i>j}\vec r^i \times \vec F^{i,j}=-\sum \sum_{i<j}\vec r^{j} \times \vec F^{i,j}
この式を、式(5)の右辺の第2項に代入して整頓すると、
\sum_i \sum_{j\neq i} \vec r^i \times \vec F^{i,j}
=\sum \sum_{i<j}(\vec r^i - \vec r^{j}) \times \vec F^{i,j}
さらに、内力に関する第2の仮定により、\vec r^i - \vec r^{j} と\vec F^{i,j}は同じ方向なので、ベクトル積の定義より、この項は、零となることが分かる。
故に、式(4)の右辺の第2項は零となり、
\vec N=\frac{d\sum_i \vec l^i}{dt} \qquad (6)
が得られる。全角運動量を\vec L =\sum_i \vec l^i とおけば、
式(6)は、次のように書ける。
命題;回転運動の関するオイラーの運動方程式
剛体の内力に上述の2つの仮定を付ける。このとき、
剛体に作用する全ての外部力の原点周りの力のモーメント\vec N=\sum_i\vec N^i=\sum_i\vec r^i \times \vec F^iと、
全角運動量\vec L =\sum_i \vec l^i =\sum_i \vec r^i \times \vec p^iの間には、
\vec N=\frac{d\vec L}{dt} \qquad (7)
この命題の導出までは詳しく述べたが、本テキストではこれ以上は深入りしない。
この先にも興味がある方は、次の記事をご覧ください。
固定軸の周りの剛体の回転運動の方程式
回転運動の運動方程式から、任意の軸の周りの回転運動の方程式が簡単に導出できる。
z軸周りの場合を例にとり、説明する。
z軸周りの回転力はT_{\vec e_z}=\vec N \cdot \vec e_zなので、
回転運動の方程式から
T_{\vec e_z}=\vec N \cdot \vec e_z
=\frac{d\vec L}{dt} \cdot \vec e_z
この式の右辺に,L=\sum_i \vec r^i \times \vec p^i を代入すると
右辺
=\frac{d\sum_i \vec r^i \times \vec p^i}{dt}\cdot \vec e_z \qquad 微分の加法性から
=(\sum_i \frac{d \vec r^i \times \vec p^i}{dt})\cdot \vec e_z \qquad 内積の加法性から
=\sum_i(\frac{d \vec r^i \times \vec p^i}{dt}\cdot \vec e_z) \qquad ベクトル積の性質8から
=\sum_i(\frac{dr^i}{dt}\times \vec p^i+\vec r^i \times \frac{d\vec p^i}{dt}
てこの原理と力のモーメント
てこの原理については、
作用線の定理
剛体の場合、作用線に沿って力の作用点を移動しても、力の作用は変わらない。何故かは、考えてみましょう。
剛体のつり合い
いくつかの力が作用し、剛体が静止したままであるか、
重心が等速直線運動を続け、重心の周りの回転が変化しない場合に、剛体(に作用している力)は釣り合っているという。