物理/質点の運動と質点系
提供: Internet Web School
(→振り子と単振動) |
(→振り子) |
||
185 行: | 185 行: | ||
==== 振り子 ==== | ==== 振り子 ==== | ||
*[[wikipedia_ja:自由振動|ウィキペディア(単振動)]] | *[[wikipedia_ja:自由振動|ウィキペディア(単振動)]] | ||
- | + | 上記の記事の運動方程式の導出は不正確なので、簡単に説明する。 | |
- | + | ||
[[File:GENPHY00010203-01.jpg|right|frame|図 振り子]] | [[File:GENPHY00010203-01.jpg|right|frame|図 振り子]] | ||
- | 図のように振り子の支点を通り水平な直線をx軸、 | + | 図のように振り子の支点を通り水平な直線をx軸、<br/> |
- | 原点を通り鉛直な直線をy軸にとる。y軸と振り子の腕のなす角度を | + | 原点を通り鉛直な直線をy軸にとる。y軸と振り子の腕のなす角度をθ、<br/> |
- | θ | + | 腕の長さをl とする。<br/> |
+ | 振り子の腕の重さは無視できると仮定し、振り子のおもりは、質量mの質点と考え運動方程式を導こう。<br/> | ||
+ | 振り子のおもり(質量mの質点と考える)の運動の軌道は、 | ||
+ | x-y平面内の、支点Oを中心とする半径l の円周上に拘束されている。 | ||
+ | |||
==== 単振動 ==== | ==== 単振動 ==== | ||
*[[wikipedia_ja:自由振動|ウィキペディア(自由振動]]中の単振動を | *[[wikipedia_ja:自由振動|ウィキペディア(自由振動]]中の単振動を |
2015年6月29日 (月) 03:39時点における版
目次[非表示] |
質点の運動と質点系
運動の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)で表すと
運動方程式は、Md2dt2→r(t)=→Fである。
座標成分表示すると
Md2dt2x(t)=0, Md2dt2y(t)=−Mg
これらの式の両辺を、Mで割ると
d2dt2x(t)=0, d2dt2y(t)=−g(1)
運動の初期状態の指定
投げ上げた瞬間を時刻t=0とおくと、
質点の初期位置は→r(0)=(x(0),y(0))=(0,0),
初期速度は→v(0)=(vx(0),vy(0))=(ucosθ,usinθ)
運動方程式を初期状態を使って解く
(1)x成分の式を解く
d2dt2x(t)=0は、速度の定義vx:=ddtx(t)から
ddtvx(t)=0と書ける。
上の式を解こう。
tで微分して零となる関数は定数なのでaと書くと、vx(t)=a (1)
ここで、速度の初期条件からvx(0)=ucosθなので
a=ucosθ(2)
次に(1)式を解こう。
tで微分してaとなるのはat+b(bは未知定数)なので、
x(t)=at+b
ここで、ⅹの初期条件から、x(0)=0 なのでb=0
故に、x(t)=at
(2)式から、x(t)=ucosθtが得られる。
(2)y成分の式を解く
d2dt2y(t)=−gは、速度を用いて表現するとddtvy(t)=−g
tで微分して−gとなる関数は−gt+c(cは未知定数)なので、
vy(t)=−gt+c(3)
速度の定義から(3)式は、ddty(t)=−gt+c
tで微分して−gt+cとなる関数は、−12gt2+ct+dなので、
y(t)=−12gt2+ct+d(4)
次に位置と速度の初期条件から係数c、dを決めよう。
(3)式を用いると、初期速度の条件vy(0)=usinθから、
usinθ=vy(0)=−g0+c=c 故にc=usinθ(5)
(4)式のt=0を代入すると yの初期値y(0)=0から、
0=−12g02+c0+d=d故にd=0(6)
(4)式に(5),(6)式を代入して、
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質点間に働く力(万有引力)を求めると、
ニュートンの万有引力の法則と同じ式だが、質量は慣性質量になり、
重量質量を用いずニュートン力学が構成できる。
振り子と単振動
振り子
上記の記事の運動方程式の導出は不正確なので、簡単に説明する。
図のように振り子の支点を通り水平な直線をx軸、
原点を通り鉛直な直線をy軸にとる。y軸と振り子の腕のなす角度をθ、
腕の長さをl とする。
振り子の腕の重さは無視できると仮定し、振り子のおもりは、質量mの質点と考え運動方程式を導こう。
振り子のおもり(質量mの質点と考える)の運動の軌道は、
x-y平面内の、支点Oを中心とする半径l の円周上に拘束されている。
単振動
- ウィキペディア(自由振動中の単振動を
見てください。
質点のつり合い
質点に力F1,,Fnが作用し、質点が静止したまま(あるいは等速直線運動)であるとき、それらの力は釣り合っているという。
釣り合いの条件は、F1+ +Fn=0です(運動の第2法則と力の合成則から導出できる)。
質点系の運動
2個以上の質点が集まって出来ている系を質点系という。
質点系というときは、各質点は密集していても、離れ離れでも良い。互いに固着しようが、自由に動けようが構わない。
すべての物質は、分子の集合と考えたり、細分化して極小部分に分け、それらの集合と考えれば、十分な精度で、質点系とみなすことができる。
そのため質点系の運動の法則を、ニュートンの運動法則から導出すれば、その応用範囲は非常に広い。
質点系の運動と重心
系の任意の2つの質点間には相互に及ぼしあう力が働いていてもよい。
この相互作用力は作用・反作用の法則から、
大きさと方向は同じで向きは逆である。
この力を質点系の”内力”という。
質点系の各質点に外部から力(外力という)が加わる時、この質点系はどんな運動をするだろうか。
質点系の各質点の位置を→ri、質量をmiとし、
質点mi に作用する外力を→fi、
mi に、他の質点mjから作用する内力を→fijとする(i,j=1…N)。
すると、各質点に対して、運動の第2法則により、
ddt(mi→vi)=→fi+∑j≠i→fij ここで→vi=d→ri/dt、
各ベクトルを自由ベクトルとみなしてi=1…Nについて加え合わせると、
d2dt2∑imi→ri=∑iddtmi→vi=∑i→fi+∑i∑j≠i→fij
=∑i→fi+∑i<j→fij+∑j<i→fij
上の式の第3項の変数名iとjを入れ替えると
=∑i→fi+∑i<j→fij+∑i<j→fji
剛体の内力についての仮定から、→fji=−→fijなので、
=∑i→fi+∑i<j(→fij−→fij)=∑i→fi
故に、
d2dt2∑imi→ri=∑i→fi
が得られる。
質点系の全質量M=∑imiと質点系に働く全外力→F=∑i→fiを用いて書きなおすと、
Md2dt2(∑imi→riM)=→F
質点系の重心→Rを →R=∑imi→riM で定義すると、
Md2dt2→R=→F
この式は、力→Fをうける質量Mの質点(位置ベクトル→R)の運動方程式と同じである。
以下の解説も参考にしてください。
複雑にみえる運動も重心の運動をみれば簡単である
体操選手の運動は、跳躍などで空中をまいながら、回転や体の屈伸、ひねりなどを行う。大変複雑で美しい。
しかし、導出した質点系の重心の運動法則から、体の重心の運動は、投射体の運動であり、放物線をえがいて移動することが分かる。
空中に飛び出た瞬間の重心の位置と速度(速さと方向・向き)で、その軌跡は完全に決まってしまうのである。
複雑さ・美しさは、重心周りの体の姿勢・ひねりとその変化、それに伴う重心周りの回転変化によって、もたらされる。
回転運動については、
[|「2.5 剛体と回転力」]と
[|「2.6 剛体の回転運動と釣合い」]で調べる。
ガリレイの相対性原理
どのような慣性座標系で観測しても力学の法則は同じであるという原理。
一つの慣性系にたいして等速度並進運動(注)する観測系を考えると、
力の働いてない物体はやはり、等速度運動するので慣性系であり、
運動の第2、第3法則、万有引力の法則、力の合成則が成立することを主張している。
(注) 座標系S′(O′−x′1x′2x′3)の任意の点(x′1,x′2,x′3)が、
座標系S(O−x1x2x3)からみると、
すべて同じ速度(d→OO′dt(t)で移動すること。
言い換えると、S'系の各座標軸上の点(1,0,0),(0,1,0),(0,0,1)が、
S系からみると皆、同じ速度(d→OO′dt(t)で移動すること。
重要な原理なので、「2.8 ガリレイ変換とガリレイの相対性原理」で詳しく考察しよう。
仕事
物体に力を加えて動かす時、力はこの物体に仕事をするという。
仕事(の量)は力の大きさと動かした距離の積に比例する。
正確には、加えられる力→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 Fの\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
\vec a \cdot(\vec b+\vec c) =\vec a \cdot \vec b+\vec a \cdot \vec c
(4)(\alpha \vec a)\cdot \vec b =\vec a \cdot (\alpha \vec b)=\alpha (\vec a \cdot \vec b)
が成り立つ。
(5)\|\vec a \cdot \vec b\| \leq \|\vec a\|\|\vec b\|
(6)ノルムの性質;\|\vec a + \vec b\| \leq \|\vec a\| + \|\vec b\|
証明は、本テキストの「8章 物理数学の8.1平面と空間のベクトル」にある。
力が変動したり、物体の移動が曲線であるときの仕事
物体に作用する万有引力は、
物体の位置により決まる。
力の場
質点がどこにあろうが、
その位置\vec xに応じて力\vec{F}(\vec x)が作用する空間を力の場という。
万有引力は位置\vec xについての連続なベクトル値関数である。
力の場の力のおこなう仕事の近似値
この力を受けて運動する物体は曲線を描いて運動することが多い(彗星、惑星の運動など)。
物体がさまざまな理由で運動を拘束され曲線を描いて動くときもある。
こうした場合、力のなす仕事をどのように決めたらよいだろうか。
運動の軌跡(向きつき曲線C)を、細かく区切って
n個の向きつき小曲線部分{C^n}_iに分ける(i=1,2,,,n)。
各小部分{C^n}_iの両端の位置ベクトルを、{\vec{x}^n}_{i-1},{\vec{x}^n}_iと書く。図参照のこと。
この分割に名前を付け\Delta^n =\{{C^n}_i \mid i=1,2,,,n\}と書く。
nを増やして分割を細かくすると、
連続的に変化する力は、各小部分{C^n}_i上でほぼ一定となる。
{C^n}_iの形状もほぼ線分とみなせるようになる。
そこで、i番目の小部分上の一点{\xi}^{n}_iを選び、
この小部分上では、力を\vec{F}({\xi}^{n}_i)で近似する。
また、i番目の小部分を向きつき線分
{\vec{x}^n}_i-{\vec{x}^n}_{i-1}
で近似する。
すると、この小部分の移動で力の行う仕事は
W_{i}({\Delta}^n,{\xi}^{n}_i)
=\vec{F}({\xi}^{n}_i)\cdot ({\vec{x}^n}_i-{\vec{x}^n}_{i-1})
で近似できる。
これをすべて加えた
W({\Delta}^n,\{{\xi}^{n}_i\}_{i=1}^{n})
=\sum_{i=1}^{n}W_{i}({\Delta}^n,{\xi}^{n}_i)
=\sum_{i=1}^{n}\vec{F}({\xi}^{n}_i)\cdot ({\vec{x}^n}_i-{\vec{x}^n}_{i-1})\qquad (1)
は、力のなした仕事の近似値を与えると考えられる。
線積分可能性の定義
ある数Wが存在して、
分割数nをどんどん増やし、小部分の長さを零に近づくて行く。
この時、分割の仕方\Delta^nや\{{\xi}^{n}_i\}_{i=1}^{n}の選び方によらず
W=\lim_{n\to \infty}W({\Delta}^n, \{{\xi}^{n}_i\}_{i=1}^{n})
となるとき、
ベクトル値関数\vec{F}(\vec x)は、曲線Cに沿って線積分可能という。
このとき、Wのことを、\vec{F}(\vec x)の曲線Cに沿った線積分といい、
W=\int_{C}\vec{F}(\vec x)\cdot{d \vec x}
と記す。
曲線に沿った仕事の定義
力\vec{F}(\vec x)が、向きつき曲線Cに沿って線積分可能であるとき、
W=\lim_{n\to \infty}W({\Delta}^n, \{{\xi}^{n}_i\}_{i=1}^{n})=\int_{C}\vec{F}(\vec x)\cdot{d \vec x}
を、力\vec{F}(\vec x)によって質点をCに沿って動かすときの、
力の行う行う仕事と定義する。
どのような条件があると線積分が可能になるか、とか、
この計算をこれ以上具体的に進めるには、
空間中の向きのついた曲線を数式で表現しなければならない。
連続な力の場のなす仕事
命題;連続な力の場のなす仕事
質点が連続な力の場から力\vec{F}(\vec x)を受け、
点Pから点Q まで曲線Cに沿って移動したとする。
この軌跡(曲線Cに、点Pから点Qへむけた向きをいれたもの)が、
実数のある閉区間I=[a,b]上で定義された、
滑らかなベクトル値関数\vec x=\vec{x}(p),(a\leq p \leq b)で表せるならば、
力のなす仕事W=\int_{C}\vec{F}(\vec x)\cdot{d \vec x}は、
W=\int_{a}^{b}\vec{F}(\vec x(p))\cdot \frac{d \vec x(p)}{dp}dp
(注)滑らかな関数とは、I=[a,b]上で連続で、
開区間I^{\circ}=(a,b)で微分可能で導関数が連続なもの。
我々が経験するすべての曲線は、
曲線を表示するための、パラメータ区間上のベクトル値関数であらわせる。
この関数が滑らかというのが、この命題の課す制約である。
証明;
曲線Cにそった、力の線積分が可能であり、Wとなることを示す。
曲線の分割を得るため、パラメータの分割を行う。
パラメータaからbまでを、n個の小区間
[p^{n}_{i-1},p^{n}_i](i=1,2,,,n)に分割する。ここで(p^{n}_0=a<p^{n}_{1}<p^{n}_{2},,,,<p^{n}_{n}=b)
この分割を\Delta^{n}:=\{[p^{n}_{i-1},p^{n}_i]\mid i=1,2,3,,,n\}と書く。
n個の小区間の長さの最大値を|\Delta^{n}|と記す。|\Delta^{n}|:=\max_{i}(p^{n}_i-p^{n}_{i-1})
\vec x^{n}_i:=\vec{x}(p^{n}_i)(i=1,2,,,,,n)をCの分割の分点とすることで、
Cのn分割C^{n}_i,(i=1,2,,,n)が得られる。これを\Delta'^nと書く。
任意の\eta^{n}_i \in [p^{n}_{i-1},p^{n}_i]を選び
\xi^{n}_i:=\vec{x}(\eta^{n}_i)
を、C^{n}_i上の力の代表点とする。
W({\Delta'}^n,\{{\xi}^{n}_i\}_{i=1}^{n})
=\sum_{i=1}^{n}W_{i}({\Delta'}^n,{\xi}^{n}_i)
=\sum_{i=1}^{n}\vec{F}({\xi}^{n}_i)\cdot ({\vec{x}^n}_i-{\vec{x}^n}_{i-1})
=\sum_{i=1}^{n}\vec{F}(\vec{x}(\eta^{n}_i))\cdot (\vec{x}(p^{n}_i)-\vec{x}(p^{n}_{i-1}))
が、パラメータの分割を細かくしていくと、
分割\Delta や \eta^{n}_iの選び方に無関係な定数に収束することを示せばよい。
p^{n}_i- p^{n}_{i-1}が小さいので、微分の定義から
\frac{\vec{x}(p^{n}_i)- \vec{x}(p^{n}_{i-1})}{p^{n}_i)- p^{n}_{i-1}}
\approx \frac{d\vec{x}}{dp} (p^{n}_{i-1})
そこで、
\vec{x}(p^{n}_i)- \vec{x}(p^{n}_{i-1})
\approx
\frac{d\vec{x}}{dp} (p^{n}_{i-1})(p^{n}_i)- p^{n}_{i-1})
故に、
W^{n}_i \approx \vec{F}(p^{n}_{i-1})\cdot \frac{d\vec{x}}{dp} (p^{n}_{i-1})(p^{n}_i)- p^{n}_{i-1})
=(\vec{F}\cdot\frac{d\vec{x}}{dp})(p^{n}_{i-1}) (p^{n}_i)- p^{n}_{i-1})
力がパラメータ区間[a,b]の間に行う仕事W^nは、
W^n=\sum_{i=1}^{n}W^{n}_i
\approx
\sum_{i=1}^{n}(\vec{F}\cdot \frac{d\vec{x}}{dp})((p^{n}_{i-1}) (p^{n}_i)- p^{n}_{i-1})
この式は、nが大きくなるほど、等式に近くなる。
最後の項\sum_{i=1}^{n}(\vec{F}\cdot\frac{d\vec{x}}{dp})((p^{n}_{i-1}) (p^{n}_i)- p^{n}_{i-1})は、
[a,b]上で定義された連続な実数値関数
f(p):=(\vec{F}\cdot\frac{d\vec{x}}{dp})(p)=\vec{F}(p)\cdot\frac{d\vec{x}}{dp}(p)
の、パラメータ区間の分割\Delta^{n}:=\{[p^{n}_{i-1},p^{n}_i]\mid i=1,2,3,,,n\}に対応する
リーマン和である。
そこで、nを無限に大きくしていくと、最後の式は
W:=\lim_{n\to \infty}W^n=\int_{a}^{b}f(p)dp
=\int_{a}^{b}(\vec{F}\cdot\frac{d\vec{x}}{dp})(p)
=\int_{a}^{b}\vec{F}(p)\cdot \frac{d\vec{x}}{dp}(p)dp
に収束する(8章の8.3 積分に証明あり)。
(注)この証明は、
「この式は、nが大きくなるほど、等式に近くなる」
ことの証明がないので、不正確である。
きちんとした証明は、8章物理数学の8.3積分で与える。
仕事の単位
仕事の定義W=\|\vec F\|\|\vec s\| \cos\thetaから、仕事の単位は、力の大きさ\|\vec F\|の単位と長さ\|\vec s\|の単位を掛けたものになる( \cos\theta は無単位なので )。
MKSA単位系では、力の大きさの単位はN(ニュートン)、長さの単位はm(メートル)なので、仕事の単位はNm となる。
これをJ(ジュール)と呼ぶ。J=Nmである。