Processing math: 72%

物理/運動の法則の応用2

提供: Internet Web School

(版間での差分)
( 回転の運動方程式の導出  )
(慣性モーメントの計算(1)リーマン和の極限を求める方法)
 
(間の22版分が非表示)
27 行: 27 行:
*[[wikipedia_ja:質点|ウィキペディア(質点系の力学)]]
*[[wikipedia_ja:質点|ウィキペディア(質点系の力学)]]
====複雑にみえる運動も重心の運動をみれば簡単である  ====
====複雑にみえる運動も重心の運動をみれば簡単である  ====
-
体操選手の運動は、跳躍や着地などで空中をまいながら、回転や体の屈伸、ひねりなどを行う。大変複雑である。<br/>
+
体操選手の運動は、跳躍などで空中をまいながら、回転や体の屈伸、ひねりなどを行う。大変複雑で美しい。<br/>
しかし、導出した質点系の重心の運動法則から、体の重心の運動は、投射体の運動であり、放物線をえがいて移動することが分かる。<br/>
しかし、導出した質点系の重心の運動法則から、体の重心の運動は、投射体の運動であり、放物線をえがいて移動することが分かる。<br/>
-
空中に飛び出た瞬間の速度(速さと方向・向き)で、その軌跡は完全に決まってしまうのである。
+
空中に飛び出た瞬間の重心の位置と速度(速さと方向・向き)で、その軌跡は完全に決まってしまうのである。
==  剛体の運動とつり合い==
==  剛体の運動とつり合い==
35 行: 35 行:
剛体(Rigid body)とは、<br/>
剛体(Rigid body)とは、<br/>
-
質点系であって、それらの、どの2質点の間の距離も変わらない,特殊な系のことを言う。<br/>
+
質点系であって、それらの、どの2質点の間の距離も変わらない,特殊な質点系のことを言う。<br/>
どの2質点の間の距離も変わらなければ変形は起こらない。<br/>
どの2質点の間の距離も変わらなければ変形は起こらない。<br/>
固くて変形しにくい物体を理想化した概念である。<br/>
固くて変形しにくい物体を理想化した概念である。<br/>
485 行: 485 行:
質点の速度 ˙x=dx(t)dt <===>剛体のe軸周りの角速度˙ϕ(t);<br/>
質点の速度 ˙x=dx(t)dt <===>剛体のe軸周りの角速度˙ϕ(t);<br/>
質点の運動量 m˙x <===>  剛体の角運動量I˙ϕ;<br/>
質点の運動量 m˙x <===>  剛体の角運動量I˙ϕ;<br/>
 +
運動方程式F=m¨x <===> Te=I¨ϕ
という、対応関係があることが分かる。<br/>
という、対応関係があることが分かる。<br/>
この節で得た固定軸まわりの回転運動の方程式から、<br/>
この節で得た固定軸まわりの回転運動の方程式から、<br/>
-
もしN=0 ならば、あらゆる軸まわりの回転力が零なので、<br/>
+
もしN=0 ならば、任意の軸まわりの回転力が零なので、<br/>
-
剛体の軸まわりの角速度が一定となることが分かる。<br/>
+
剛体の任意の軸まわりの角加速度が零、角速度が一定となることが分かる。<br/>
=====剛体の回転の運動エネルギー  =====
=====剛体の回転の運動エネルギー  =====
516 行: 517 行:
ここでIは、軸まわりの、振り子の慣性質量。
ここでIは、軸まわりの、振り子の慣性質量。
-
=====慣性モーメントの計算1(一次元の剛体) =====
+
=====剛体の慣性モーメントの計算(一次元の剛体) =====
-
細長い棒Vの、軸まわりの慣性モーメントを具体的に計算しよう。<br/>
+
剛体Vは、ごく細く、まっすぐな棒で,<br/>
-
Vは、ごく細い一様な質量密度のまっすぐな棒とする。<br/>
+
長さl、質量密度(単位長さあたりの質量)は一定でρとする。<br/>
-
長さをl、質量をMとすると、単位長さ当たりの質量は$\rho=M/l$<br/>
+
棒の左端からl1の場所Oを通り、棒に直交する軸まわりの慣性モーメントを具体的に計算しよう。<br/>
-
 
+
-
軸は棒と直角で、左端からl1の場所Oを通るとする。<br/>
+
Oを原点とし、棒と同じ方向の数直線を考え、これを座標系として採用。<br/>
Oを原点とし、棒と同じ方向の数直線を考え、これを座標系として採用。<br/>
V=[a=l1,b=ll1]と表現する。<br/>
V=[a=l1,b=ll1]と表現する。<br/>
529 行: 528 行:
I=imi(ˆri)2で定義した。但しx0=a,xn=b<br/>
I=imi(ˆri)2で定義した。但しx0=a,xn=b<br/>
-
======'''慣性モーメントの近似式とリーマン和''' ======
+
======'''剛体の分割と慣性モーメントの近似式・リーマン和''' ======
Viの質量miは、Viの長さxixi1に質量密度ρを掛ければ得られるので<br/>
Viの質量miは、Viの長さxixi1に質量密度ρを掛ければ得られるので<br/>
-
mi=ρ(xixi1)となり、問題なく定まる。<br/>
+
mi=ρ(xixi1)であり、<br/>
 +
I=iρ(ˆri)2(xixi1)<br/>
 +
と書ける。<br/>
しかし、剛体V=[a,b]をいくら細かく分割しても、<br/>
しかし、剛体V=[a,b]をいくら細かく分割しても、<br/>
各小区間Vi=[xi1,xi]は大きさ(長さ)をもつので、<br/>
各小区間Vi=[xi1,xi]は大きさ(長さ)をもつので、<br/>
540 行: 541 行:
imi(ˆri)2=iρ(ξi)2(xixi1)<br/>
imi(ˆri)2=iρ(ξi)2(xixi1)<br/>
で近似される。<br/>
で近似される。<br/>
-
これは分割の仕方と分割小区間の代表点ξiの選び方によって変化する。<br/>
+
 
そこで、この分割を<br/>
そこで、この分割を<br/>
Δ={Vi=[xi1,xi]i=1,2,,,,n}
Δ={Vi=[xi1,xi]i=1,2,,,,n}
547 行: 548 行:
で,慣性モーメントの近似式を表すことにする。<br/>
で,慣性モーメントの近似式を表すことにする。<br/>
すると、<br/>
すると、<br/>
-
分割Δとその小区間Viから代表点ξiを選んだときの<br/>
 
'''慣性モーメントの近似式'''は、<br/>
'''慣性モーメントの近似式'''は、<br/>
-
IΔ(ξ1,,,ξn)=ni=1ρ(ξi)2(xixi1)<br/>
+
$I^{\Delta}(\xi_1,,,\xi_n)=\sum_{i=1}^{n}\rho(\xi_i)^2(x_i-x_{i-1})\qquad  (1)$<br/>
と書ける。<br/>
と書ける。<br/>
 +
この値は分割の仕方と分割小区間の代表点ξi(Vi)の選び方によって変化する。<br/>
 +
Viの中で、原点に最も近い点ξmi(Vi)(i=1,2,n)にとると<br/>
 +
最小値
 +
Im(Δ):=ni=1ρ(ξmi)2(xixi1)<br/>
 +
をとり、<br/>
 +
Viの中で、原点に最も遠い点ξMi(Vi)(i=1,2,n)にとると<br/>
 +
最大値 IM(Δ):=ni=1ρ(ξMi)2(xixi1)<br/>
 +
を取る。<br/>
関数y=f(x)=ρx2を使って表現すれば、<br/>
関数y=f(x)=ρx2を使って表現すれば、<br/>
IΔ(ξ1,,,ξn)=ni=1f(ξi)(xixi1)<br/>
IΔ(ξ1,,,ξn)=ni=1f(ξi)(xixi1)<br/>
-
数学の分野では、この和の式は、<br/>
+
であり、Im(Δ)IΔ(ξ1,,,ξn)IM(Δ)を満たす。<br/>
-
分割Δ={Vi=[xi1,xi]i=1,2,,,,n}Viの代表点ξi,(i=1,2,,,n)に関する関数y=f(x)の'''リーマン和'''と呼ばれ、<br/>
+
質量密度が場所で変わるときは、、関数はy=f(x)=ρ(x)x2になり、<br/>
-
$I^{f,\Delta}(\xi_1,,,\xi_n)$<br/>
+
剛体の重心を求めるときは、後述するように、別の関数が現れる。<br/>
-
と表示する。<br/>
+
そこで、数学の分野では、一般の関数y=f(x)にたいして<br/>
-
関数y=f(x)は任意にきめられるので、どの関数のリーマン和か分かるように、明示してある。<br/>
+
If,Δ(ξ1,,,ξn)=ni=1f(ξi)(xixi1)(2)<br/>
-
 
+
を求め、分割Δ={Vi=[xi1,xi]i=1,2,,,,n}Viの代表点ξi,(i=1,2,,,n)に関する関数y=f(x)の'''リーマン和'''と呼ぶ。<br/>
 +
その最小値$I_{m}(f,\Delta)I_{M}(\Delta)$も,同様に定義される。<br/><br/>
 +
Im(f,Δ)If,Δ(ξ1,,,ξn)IM(Δ)(3)
 +
慣性モーメントの近似式(1)は、関数$y=f(x)=\rho(x) x^2$にたいするリーマン和である。
======慣性モーメントの近似式の意味 ======
======慣性モーメントの近似式の意味 ======
568 行: 579 行:
y=f(x)=ρx2を代入すれば、考察対象の剛体の慣性モーメントの話になる。<br/>
y=f(x)=ρx2を代入すれば、考察対象の剛体の慣性モーメントの話になる。<br/>
リーマン和<br/>
リーマン和<br/>
-
IΔ(ξ1,,,ξn)=ni=1f(ξi)(xixi1)<br/>
+
$I^{f,\Delta}(\xi_1,,,\xi_n)=\sum_{i=1}^{n}f(\xi_i)(x_i-x_{i-1})$<br/>
は、y=f(x)のグラフを、棒グラフで近似したときの棒グラフの作る面積(各角柱の面積和)であることが分かる。図参照。<br/>
は、y=f(x)のグラフを、棒グラフで近似したときの棒グラフの作る面積(各角柱の面積和)であることが分かる。図参照。<br/>
-
$y=f(x)$のグラフとx軸およびy軸と平行な直線x=ax=bで囲まれる部分の面積を近似している。<br/>
+
また、$I_{m}(f,\Delta)$は一点鎖線でしめす、小さいほうの長方形の和であり、<br/>
-
分割を限りなく細かくしていくとき、この値が分割や代表点の選び方に関係ない数に収束するならば、
+
$I_{Mm}(f,\Delta)$は点線でしめす、大きいほうの長方形の和である。
-
その極限値は、$y=f(x)$のグラフとx軸およびy軸と平行な直線x=ax=bで囲まれる部分の面積と考えられる。<br/>
+
-
======可積分条件======
 
-
この節は、区間上で定義された関数の
 
-
積分可能な条件を紹介する。<br/>
 
-
大学の教養コース程度の数学を使うが、テキスト中で理解できるように説明する。<br/>
 
-
興味のない方は、とばしてください。<br/>
 
-
<br/>
 
-
'''(1)準備;集合論の初歩'''<br/>
+
If,Δ(ξ1,,,ξn)は、y=f(x)のグラフとx軸およびy軸と平行な直線x=ax=bで囲まれる部分の面積$S$を近似している。<br/>
-
 
+
また、$I_{m}(f,\Delta)\leq S \leq I_{Mm}(f,\Delta)  \qquad (4)  $<br/>
-
集合論の初歩の知識は前提にして記述するので、<br/>
+
-
なじみのない方は、下記を参考に、<br/>
+
-
集合の素朴な定義、集合の和集合や共通集合、集合の包含関係などについて学習してほしい。
+
-
*[[wikipedia_ja:集合 |ウィキペディア(集合)]]
+
-
 
+
-
'''(2)区間上の関数のリーマン和'''<br/>
+
-
'''定義;リーマン和 ''' <br/>
+
-
区間V=[a,b]で定義され、実数に値をとる関数y=f(x)を考える。<br/>
+
-
この区間の分割<br/>
+
-
Δ={Vi=[xi1,xi]i=1,2,,,n},x0=a,xn=b<br/>
+
-
と、その代表点ξiVi(i=1,2,,,n)に関する、y=f(x)のリーマン和とは、
+
-
$I^{f,\Delta}(\xi_1,,,\xi_n)=$<br/>
+
-
if(ξi)v(Vi)=if(ξi)(xixi1) <br/>
+
-
で定義する。<br/>
+
-
リーマン和は、y=f(x)のグラフを、棒グラフで近似したときの棒グラフの作る面積(各角柱の面積和)であることが分かる。図参照。<br/>
+
-
y=f(x)のグラフとx軸およびy軸と平行な直線x=ax=bで囲まれる部分の面積を近似している。<br/>
+
-
 
+
-
 
+
-
'''(3)リーマン和の不足リーマン和と過剰リーマン和による評価'''<br/>
+
-
$V$を分割して得られた小区間Vi=[xi1,xi]を考える。<br/>
+
-
関数y=f(x)をこの小区間上に限定した時、<br/>
+
-
関数は、この区間上の点で最大値と最小値をとると仮定する(注参照)。<br/>
+
-
関数の最大値max{f(x)xVi}と最小値min{f(x)xVi}を、それぞれ、m(f;Vi),M(f;Vi)と書く。<br/>
+
-
(注) 区間上で最大値、最小値を取らない関数では、<br/>
+
-
[[wikipedia_ja:有界函数 |有界な関数]]でありさえすれば、常に存在し、しかも最大値、最小値と殆ど同じ性質をもつ、<br/>
+
-
上限、下限に置き換えれば以後の、議論は成り立つ。上限、下限については「(5)不足リーマン和の上限と過剰リーマン和の下限」で説明する。<br/><br/>
+
-
 
+
-
すると、Viの任意の点ξ に対して、<br/>
+
-
m(f;Vi)f(ξ)M(f;Vi)  <br/>
+
-
故に、<br/>
+
-
'''補題1'''<br/>
+
-
ⅰ)どのような代表点{ξi}i,(ξiVi,i=1,2,,,n)に対しても<br/>
+
-
$I_{m}(f,\Delta):=I^{f,\Delta}(\xi_{1}^m,,,\xi_{n}^m)
+
-
=\sum_i m(f;V_i)v(V_i)$<br/>
+
-
$\leq  
+
-
I^{f,\Delta}(\xi_1,,,\xi_n)=\sum_i f(\xi_i)v(V_i)$<br/>
+
-
$\leq
+
-
\sum_i M(f;V_i)v(V_i)
+
-
=I_{M}(f,\Delta)=I^{f,\Delta}(\xi_{1}^M,,,\xi_{n}^M)       \qquad   (1)$ <br/>
+
-
そこで、Im(f,Δ)を'''Δ)に関するfの'''不足リーマン和'''、IM(Δ)を'''過剰リーマン和'''と呼ぶ。<br/>
+
-
ⅱ)Im(f,Δ)=minξiVi,i=1,2,,,nIf,Δ(ξM1,,,ξMn) <br/>
+
-
IM(f,Δ)=maxξiVi,i=1,2,,,nIf,Δ(ξM1,,,ξMn) <br/>
+
-
証明は明らかなので省略。<br/>
+
-
 
+
-
'''(4)分割の細分とリーマン和の評価式'''<br/>
+
-
'''定義;分割の細分'''<br/>
+
-
Vの分割Δが分割Δの細分というのは、<br/>
+
-
Δの分点の集合{x0,x1,,,,xn}が、<br/>
+
-
Δの分点の集合{x0,x1,,,,xn}に真に含まれることと定義する。<br/>
+
-
記号でかけば、$\{x_0,x_1,,,,x_n\}\subset \{x'_0,x'_1,,,,x'_{n'}\},
+
-
\{x_0,x_1,,,,x_n\}\neq \{x'_0,x'_1,,,,x'_{n'}\}$。
+
-
記号では、ΔΔと記す。<br/><br/>
+
-
 
+
-
 
+
-
'''補題2'''<br/>
+
-
Δ Δという分割に対し、<br/>
+
-
$I_{m}(f,\Delta)
+
-
\leq I_{m}(f,\Delta')
+
-
\leq I_{M}(f,\Delta')
+
-
\leq I_{M}(f,\Delta)      \qquad  (2)$
+
-
<br/>
+
-
が成り立つ。    <br/>
+
-
(証明)<br/>
+
-
Δの小区間Vi=[xi1,xi]が分割Δでは、<br/>
+
-
{Vj=[xi1,xj],Vj+1=[xj,xi]}の2つに分割されたとする。<br/><br/>
+
-
すると、区間上の関数の最大値と最小値の定義から、<br/>
+
-
m(f;Vi)m(f;Vj)  m(f;Vi)m(f;Vj+1)<br/>
+
-
M(f;Vi)M(f;Vj)  M(f;Vi)M(f;Vj+1)<br/>
+
-
これらから、命題は成立することが分かる。<br/><br/>
+
-
 
+
-
'''(5)不足リーマン和の上限と過剰リーマン和の下限'''
+
-
補題2から、分割の細分を繰り返していくと、その分割に対応する、<br/>
+
-
不足リーマン和は、広義増加(増加するか、同じ値にとどまる)し、<br/>
+
-
過剰リーマン和は、広義減少する。<br/>
+
-
分割を細かくしていったとき、これらの極限が一致すれば、補題1から、<br/>
+
-
リーマン和の極限値は、代表点に無関係に、定まることになる。<br/>
+
-
 
+
-
そこで色々な分割に対応する不足リーマン和のなかの最大値と過剰リーマン和の最小値を求めることが、重要になる。しかし一般にはこれらは存在しないことが示せる。<br/>
+
-
そこで最大値に近い性質を持つ上限と最小値に近い下限という概念を利用する。<br/>
+
-
'''定義;上界と下界'''<br/>
+
-
Rを、全ての実数を要素とする集合とし、<br/>
+
-
Aをその部分集合とする。<br/>
+
-
実数uAの上界(upper bound)とは、<br/>
+
-
任意のaAに対して、auがなりたつこと。<br/>
+
-
実数lAの下界(lower bound)とは、<br/>
+
-
任意のaAに対して、laがなりたつこと。<br/>
+
-
UAAの上界をすべて集めた集合、<br/>
+
-
LAAの上界をすべて集めた集合とする。<br/>
+
-
UAが空集合でない(すなわち、Aの上界が少なくとも一つ存在する)とき、<br/>
+
-
Aは'''上に有界'''であるといい、<br/>
+
-
LAの時、Aは'''下に有界'''であるという。<br/>
+
-
上に有界で、下にも有界な集合(R)は、'''有界'''という。<br/>
+
-
 
+
-
'''実数の連続の公理'''<br/>
+
-
以下の性質は、色々な極限の存在の根拠を与えるもので、実数の持つ最も重要な性質の一つである。<br/>
+
-
ARとする。<br/>
+
-
もし、UAならば、UAは、最小元を持つ。<br/>
+
-
これをAの'''上限(supremum)'''あるいは'''最小上界(least upper bound)'''という。<br/>
+
-
もし、LAならば、LAは、最大元を持つ。<br/>
+
-
これをAの'''下限(infimum)'''あるいは'''最大下界(greatest lower bound)'''というという。<br/><br/>
+
-
 
+
-
補題3<br/>
+
-
uA(R) の上限となるための必要十分条件は、<br/>
+
-
ⅰ)uAの上界。すなわち任意のaAにたいしてau   <br/>
+
-
ⅱ)x<uである任意のxAの上界ではない。すなわち、x<aとなるaAが存在。<br/>
+
-
ⅲ)Aが最大値を持つ場合には、上限は最大値と一致する。<br/>
+
-
同様に、lA の下限となるための必要十分条件は、<br/>
+
-
ⅰ)lAの下界。すなわち任意のaAにたいしてla   <br/>
+
-
ⅱ)l<xである任意のxAの下界ではない。すなわち、a<xとなるaAが存在。<br/>
+
-
ⅲ)Aが最小値を持つ場合には、下限は最小値と一致する。<br/><br/>
+
-
A の上限をsupA、下限をinfAと書く。<br/><br/>
+
-
証明は、上限、下限の定義から、明らかなので省略する。<br/>
+
-
例;A=(0,1)のとき、supA=1,infA=0。<br/>
+
-
これらは、ともにAの要素でないので、<br/>
+
-
上限1はAの最大元(最大値)ではなく、下限0はAの最小元(最小値)ではない。<br/>
+
-
A=[0,1]のとき、supA=1,infA=0。<br/>
+
-
これらは、ともにAの要素なので、<br/>
+
-
上限は最大限であり、下限は最小限となる。<br/>
+
-
 
+
-
補題4.
+
-
ABRで、Bは有界集合とする。<br/>
+
-
このとき、infBinfAsupAsupB<br/>
+
-
証明は容易である。<br/><br/>
+
-
 
+
-
関数y=f(x)が連続でない時は、区間上で最大値や最小値を取らないことがある。
+
-
この場合も考慮して、最大値を上限に、最小値を下限に置き換えて、m(f;Vi)=inf{f(x)xVi},M(f;Vi)=sup{f(x)xVi}で定義すれば、
+
-
有界関数に対して、これらは常に定義され、今までの議論はすべて成り立つ。
+
-
 
+
-
 
+
-
''' 2つの分割の共通の細分 '''<br/>
+
-
分割Δの分点の集合{xjj=1,2,,,m}と、<br/>
+
-
分割Δ の分点の集合{xjj=1,2,,,n}の<br/>
+
-
和集合{xjj=1,2,,,m}{xjj=1,2,,,n}を分点とする分割をΔΔと書く。<br/>
+
-
すると新しい分割は<br/>
+
-
ΔΔΔ と
+
-
ΔΔΔ<br/>
+
-
を満たす。<br/>
+
-
これを用いると、不足リーマン和の上限s(f)と過剰リーマン和の下限S(f)が存在することが証明できる。<br/>
+
-
 
+
-
補題5<br/>
+
-
fを区間V=[a,b]で定義され実数値をとる有界関数(すなわち、{f(x)xV}Rの有界集合)とする。<br/>
+
-
V=[a,b]の分割を全て集めて作った集合をD(V)と書く。<br/>
+
-
すると、<br/>
+
-
ⅰ)任意のΔ,ΔD(V)に対して、
+
-
Im(f,Δ)IM(f,Δ))<br/>
+
-
ⅱ)集合{Im(f,Δ)ΔD(V)}は上に有界、<br/>
+
-
集合{IM(f,Δ)ΔD(V)}は下に有界<br/>
+
-
ⅲ)s(f):=sup{Im(f,Δ)ΔD(V)}
+
-
<br/>
+
-
S(f):=inf{IM(f,Δ)ΔD(V)}は存在し、<br/>
+
-
s(f) S(f) <br/>
+
-
証明;<br/>
+
-
ⅰ)ΔΔΔ なので、補題2から、<br/>
+
-
$I_m(f,\Delta)
+
-
\leq I_m(f,\Delta\vee {\Delta}') 
+
-
\leq I_M(f,\Delta\vee {\Delta}') 
+
-
\leq I_M(f,{\Delta}') $<br/>
+
-
ⅱ)1)で証明した不等式で、分割Δ は固定する。<br/>
+
-
すると全ての分割 Δに対して、Im(f,Δ)IM(f,Δ))なので<br/>
+
-
集合{Im(f,Δ)ΔD(V)}は、上界IM(f,Δ))を持ち、上に有界である。<br/>
+
-
後者も同様にして下に有界であることが示せる。<br/>
+
-
ⅲ)従って、実数の連続性の公理から、<br/>
+
-
集合{Im(f,Δ)ΔD(V)}は上限s(f)をもち、<br/>
+
-
集合{IM(f,Δ)ΔD(V)}は下限S(f)をもつ。<br/>
+
-
上限は、上界の中の最小値なので、<br/>
+
-
s(f)IM(f,Δ)<br/>
+
-
この式は任意のΔについて成立するので、<br/>
+
-
s(f)は、集合{IM(f,Δ)ΔD(V)}の下界である。<br/>
+
-
下限S(f)は、下界のなかの最大値なのでs(f)S(f)を得る。<br/>
+
-
 
+
-
 
+
-
 
+
-
'''(6)分割を細かくしていくと不足リーマン和はその上限sに、過剰リーマン和はその下限Sに収束する'''<br/>
+
-
この命題を正確に述べるには、
+
-
まず、分割Δ={Vi=[xi1,xi]i=1,2,,,n}の大きさを、きちんと定める必要がある。<br/>
+
-
定義;Δ={Vi=[xi1,xi]i=1,2,,,n}の大きさd(Δ)とは、<br/>
+
-
この分割で得られた小区間の長さの、最大値で定義する。記号で書くと<br/>
+
-
d(Δ)=max{xixi1i=1,2,,,n}
+
-
 
+
-
定理(ダルブー;Darboux)<br/>
+
-
V=[a,b]<br/>
+
-
fを、Vで定義され、実数に値を取る有界関数とする。<br/>
+
-
このとき、<br/>
+
-
ⅰ)limd(Δ)0Im(f,Δ)=s(f)<br/>
+
-
ⅱ)limd(Δ)0IM(f,Δ)=S(f)<br/>
+
-
証明;<br/>
+
-
ⅰ)を示す。ⅱ)は同じようにして証明できるので略す。<br/>
+
-
これを示すには、どんなに小さい正の実数ϵに対しても、それに応じた、小さい正の実数δϵを適切に選べば、分割の大きさがδϵより小さい、どんな分割Δも、s(f)Im(f,Δ)<ϵであることを示せばよい。以下に、数段階に分けて、これを証明する。<br/>
+
-
 
+
-
$\quad 1)<br/>
+
-
$上限の性質(補題3)から、<br/>
+
-
ある分割<br/>
+
-
D={VDi=[xDi1,xDi]i=1,2,,,n}inD(V)<br/>
+
-
が存在して、<br/>
+
-
s(f)Im(f,D)<ϵ2(1)<br/>
+
-
今後このDを使って、証明を進める。<br/>
+
-
 
+
-
2)<br/>
+
-
分割Dの小区間VDiの長さ(xDixDi1)(i=1,2,,,n)
+
-
最小値をeとおくと<br/>
+
-
e=minni=1(xDixDi1) <br/>
+
-
eに比べて非常に小さい大きさを持つ分割、<br/>
+
-
Δ={VΔi=[xΔi1,xΔi]i=1,2,,,N}
+
-
<br/>
+
-
d(Δ)=maxi=1,2,,,N(xΔixΔi1)e<br/>
+
-
<br/>
+
-
を考える。<br/>
+
-
もし、DΔならば補題2より、Im(f,D)Im(f,Δ)、<br/>
+
-
すると$\mathscr{s}(f)-I_m(f,\Delta)\leq \mathscr{s}(f)-I_m(f,D) \leq
+
-
\frac{\epsilon}{2}\leq \epsilon$ <br/>
+
-
通常、分割Δは、Dの細分になっていない。<br/>
+
-
この場合は、いくつか(高々n-1個)のΔの小区間が、Dの小区間には含まれず、<br/>
+
-
Dの分点xDi(i=1,2,,,n1)をまたぐことになる。図参照のこと。<br/>
+
-
議論を簡単にするため、Dの分点xDi(i=1,2,,,n1)が全て、Δの小区間によって跨がれていると仮定し、議論を進める。<br/>
+
-
他のケースでも、証明はおなじようにできるので、このように仮定しても何の問題も起こらない。<br/>
+
-
Dの分点xDiを跨ぐΔの小区間をVΔmiとする(i=1,2,,,n-1)。<br/>
+
-
3)  <br/>
+
-
2つの分割DΔからΔ:=DΔを作る。<br/>
+
-
すると<br/>
+
-
Δ={VΔ1,VΔ2,,,,,,,,,VΔm11,<br/>
+
-
[xΔm11,xD1],[xD1,xΔm1],<br/>
+
-
VΔm1+1,VΔm1+2,,,,,,,,,VΔm21,<br/>
+
-
[xΔm21,xD2],[xD2,xΔm2],<br/>
+
-
VΔm2+1,VΔm2+2,,,VΔm31,<br/><br/>
+
-
,,,,,,,,,<br/><br/>
+
-
VΔmn1+1,VΔmn1+2,,,,,,,,,VΔN} (2)<br/>
+
-
と書ける。<br/>
+
-
 
+
-
ΔΔで、 DΔ なので、<br/>
+
-
Im(f,Δ)Im(f,Δ), Im(f,D)Im(f,Δ)<br/>
+
-
後者の式から、<br/>
+
-
0s(f)Im(f,Δ)s(f)Im(f,D)<br/>
+
-
この式と(1)式から、<br/>
+
-
0s(f)Im(f,Δ)<ϵ2<br/>
+
-
そこで、<br/>
+
-
d(Δ)0ならば、Im(f,Δ)Im(f,Δ)<ϵ2
+
-
<br/>
+
-
が示せれば、<br/>
+
-
0s(f)Im(f,Δ)<br/>
+
-
=(s(f)Im(f,Δ)+(Im(f,ΔIm(f,Δ)ϵ<br/>
+
-
が示され、証明が終わる。<br/>
+
-
$\quad 4)$ <br/>
+
-
Im(f,Δ)=Ni=1m(f;VΔi)v(VΔi)
+
であり、<br/>
であり、<br/>
-
(2)式から、<br/>
+
$I_{m}(f,\Delta)$は面積を下から評価し、<br/>
-
$I_m(f,{\Delta}')$<br/>
+
$I_{M}(f,\Delta)$は面積を上から評価していることがわかる。<br/>
-
$=\sum_{i\notin \{m_1,m_2,,,,m_{n-1}\}} m(f;V^{\Delta}_i)v(V^{\Delta}_i)$<br/>
+
分割を限りなく細かくしていくとき、<br/>
-
$+\sum_{k=1}^{n-1} m(f;[x^{\Delta}_{m_{k}-1},x^{D}_k])v([x^{\Delta}_{m_{k}-1},x^{D}_k])$<br/>
+
リーマン和が分割や代表点の選び方に関係ない数に収束するならば、<br/>
-
$+\sum_{k=1}^{n-1} m(f;[x^{D}_k,x^{\Delta}_{m_{k}}])v([x^{D}_k,x^{\Delta}_{m_{k}}])$<br/>
+
その極限値は、<br/>
-
なので、<br/>
+
$y=f(x)$のグラフとx軸およびy軸と平行な直線$x=a$$x=b$で囲まれる部分の面積<br/>
-
Im(f,Δ)Im(f,Δ)<br/>
+
と考えられる。<br/>
-
=n1k=1m(f;[xΔmk1,xDk])v([xΔmk1,xDk])<br/>
+
もし、分割Δを細かくしていくとき<br/>
-
$+\sum_{k=1}^{n-1} m(f;[x^{D}_k,x^{\Delta}_{m_{k}}])v([x^{D}_k,x^{\Delta}_{m_{k}}])$<br/>
+
$I_{m}(f,\Delta)$$I_{M}(f,\Delta)$が同じ値に収束することが示せれば、<br/>
-
$-\sum_{i\in \{m_1,m_2,,,,m_{n-1}\}} m(f;V^{\Delta}_i)v(V^{\Delta}_i)$<br/>
+
(3)式と(4)式から、リーマン和は、関数のグラフの作る面積$S$に収束することが分かった。
-
関数はV上で有界なので、適切に正の実数Mを選ぶと、xVの要素ならば<br/>
+
-
|f(x)|Mが成立する。<br/>
+
-
すると$|m(f;[x^{\Delta}_{m_{k}-1},x^{D}_k])|, |m(f;[x^{D}_k,x^{\Delta}_{m_{k}}])|
+
-
\leq M$<br/>
+
-
が成り立つ。また<br/>
+
-
$v(V^{\Delta}_{m_k})
+
-
=v([x^{\Delta}_{m_{k}-1},x^{D}_k])+v([x^{D}_k,x^{\Delta}_{m_{k}}])$で、<br/>
+
-
$v(V^{\Delta}_i)\leq d(\Delta) $<br/>
+
-
なので<br/>
+
-
$|I_m(f,{\Delta}')-I_m(f,\Delta)|\leq 2M\sum_{i\in \{m_1,m_2,,,,m_{n-1}\}} v(V^{\Delta}_i)\leq 2M(n-1)d(\Delta)$<br/>
+
-
そこで、<br/>
+
-
δϵ=ϵ4Mn
+
-
と選べば、<br/>
+
-
$d(\Delta)\leq \delta_{\epsilon}\Delta$も、<br/>
+
-
$0\leq I_m(f,{\Delta}')-I_m(f,\Delta)|\leq \frac{\epsilon}{2}$<br/>
+
-
を満たすことが証明できた。証明終わり。<br/>
+
 +
======可積分の定義と積分 ======
 +
「分割を細かくしていくとき、リーマン和が収束する」ということは、<br/>
 +
面積を決める上で決定的に重要がことなので、<br/>
 +
可積分という名を付けて、数学的に厳密に定義する。
 +
このためにはまず、分割の大きさを定める必要がある。<br/>
 +
'''定義:分割の大きさ'''<br/>
 +
分割 Δ={Vi=[xi1,xi]i=1,2,,,,n}の大きさとは、<br/>
 +
d(Δ):=maxi=1,2,n(xixi1)<br/>
-
  '''(7)可積分条件'''<br/>
+
'''定義:可積分と積分'''<br/>
-
定理;可積分条件 <br/>
+
fを、有界閉区間V上で定義され、実数の値をとる関数とする。<br/>
-
V=[a,b]<br/>
+
-
fを、Vで定義され、実数に値を取る有界関数とする。<br/>
+
-
次の条件のうち1つが成立すれば、残り2つは成立する(互いに同値という)。
+
-
ⅰ)fV上で(リーマン)可積分<br/>
+
-
ⅱ)limd(Δ)0(IM(f,Δ)Im(f,Δ))=0<br/>
+
-
ⅲ)S(f)=s(f)<br/>
+
-
<br/>
+
-
証明<br/>
+
-
ⅰ)を仮定する。ⅱ)が成立することを示そう。<br/>
+
-
fの積分値をαとおくと、可積分の定義から、<br/>
+
-
任意のϵ>0に対して、δ>0が存在して、<br/>
+
-
d(Δ)<δである任意の分割と、その分割の任意の代表点ξi,(i=1,2,,,)に対し,<br/>
+
-
|If,Δ(ξ1,,,ξn)α|<12ϵ<br/>
+
-
が成立する。<br/>
+
-
変形すると<br/>
+
-
$\alpha-\frac{1}{2}\epsilon
+
-
<I^{f,\Delta}(\xi_1,,,\xi_n)
+
-
<\alpha+\frac{1}{2}\epsilon \qquad (1)  $<br/>
+
-
ここで、補題1のⅱ)から、<br/>
+
-
inf{ξi}If,Δ(ξ1,,,ξn)=Im(f,Δ)<br/>
+
-
sup{ξi}If,Δ(ξ1,,,ξn)=IM(f,Δ)<br/>
+
-
なので、<br/>
+
-
(1)式から、<br/>
+
-
$\alpha-\frac{1}{2}\epsilon
+
-
\leq
+
-
I_{m}(f,\Delta)
+
-
\leq
+
-
I_{M}(f,\Delta)
+
-
\leq
+
-
\alpha+\frac{1}{2}\epsilon$<br/>
+
-
これより、任意のϵ>0に対して、δ>0が存在して、<br/>
+
-
d(Δ)<δ(0IM(f,Δ)Im(f,Δ)ϵ)<br/>
+
-
ⅱ)が示せた。<br/>
+
-
ⅱ)を仮定する。 ⅲ)が成り立つことを示す。<br/>
+
-
$I_{m}(f,\Delta)
+
もし、ある実数$I$が存在して、<br/>
-
\leq
+
どんな分割$\Delta=\{V_i=[x_{i-1},x_i] \mid i=1,2,,,,n\} $と代表点$\xi_i\in V_i(i=1,2,\cdots ,n)$であっても、<br/>
-
\mathscr{s}(f):=\sup_{\Delta}I_{m}(f,\Delta)
+
$\lim_{d(\Delta) \to 0}I^{f,\Delta}(\xi_1,,,\xi_n)=I$<br/>
-
\leq
+
が成り立つ時、<br/>
-
\mathscr{S}(f):=\inf_{\Delta}I_{M}(f,\Delta)
+
f$V$上で(リーマン)可積分であるという。<br/>
-
\leq
+
このとき、$I$ f$V$上での(リーマン)積分といい、<br/>
-
I_{M}(f,\Delta)$<br/>
+
$I=\int_{V}f=\int_{V} f(x)dx$<br/>
-
なので、<br/>
+
などと書く。
-
$0
+
-
\leq
+
-
\mathscr{S}(f)-\mathscr{s}(f)
+
-
\leq
+
-
I_{M}(f,\Delta)-I_{m}(f,\Delta)$<br/>
+
-
故に、分割を細かくしていき、極限をとると、<br/>
+
-
$0
+
-
\leq
+
-
\mathscr{S}(f)-\mathscr{s}(f)
+
-
\leq
+
-
\lim_{d(\Delta)\to 0}(I_{M}(f,\Delta)-I_{m}(f,\Delta))$<br/>
+
-
ⅱ)が成立するので、<br/>
+
-
$=0$<br/>
+
-
ⅲ)が示せた。<br/>
+
-
ⅲ)を仮定する。 $\alpha=\mathscr{S}(f)=\mathscr{s}(f)$とおく。<br/>
+
-
ⅰ)が成り立つことを示そう。<br/>
+
-
補題1のⅰ)から、どのような分割Δと、その代表点$\{\xi_i\}_{i}, (\xi_i \in V_i,i=1,2,,,n)$に対しても<br/>
+
-
$I_{m}(f,\Delta)
+
-
\leq
+
-
I^{f,\Delta}(\xi_1,,,\xi_n)
+
-
\leq
+
-
I_{M}(f,\Delta)$<br/>
+
-
ここで、ダルブーの定理から、<br/>
+
-
$\lim_{d(\Delta) \to 0}I_{m}(f,\Delta)=\mathscr{s}(f)=\alpha$,<br/>
+
-
$\lim_{d(\Delta) \to 0}I_{M}(f,\Delta)=\mathscr{S}(f)=\alpha$<br/>
+
-
が成り立つので、<br/>
+
-
$\lim_{d(\Delta) \to 0}I^{f,\Delta}(\xi_1,,,\xi_n)=\alpha$ <br/>
+
-
が成り立つ。<br/>
+
-
ⅰ)が示せた。<br/>
+
-
 
+
-
  '''(8)有限個の点を除いて連続な閉区間上の関数は積分可能'''<br/>
+
-
色々な関数のグラフを書くとつながっているところを、跳んでいるところが出来る。<br/>
+
-
y=Xのグラフはずっとつながっている。<br/>
+
-
x<0のとき  y=0, 0xのとき  y=1 で関数を定義すると、x=0のところで
+
-
そのグラフは跳んでいる。<br/>
+
-
連続や不連続は関数の非常に重要な性質であり、それを調べることはとても豊かな知識をもたらす。<br/>
+
-
しかし正確に議論するには、連続とは何かをきちんと定義する必要がある。<br/>
+
-
関数の連続性の定義;<br/>
+
-
 
+
-
実数値関数 $f(x)$ がある点''' $x_0$で連続'''であるとは、<br/>
+
-
xx0 に限りなく近づくならば、$f(x)f(x_0)\lim_{x\to x_0} f(x) = f(x_0)$と記す。<br/>
+
-
 
+
-
これはイプシロン-デルタ論法(ε-δ論法)を用いれば次のように定式化できる。<br/>
+
-
(小さな)正の数 ε が任意に与えられたとき、<br/>
+
-
(小さな)正の数 δ をうまくとってやれば、<br/>
+
-
$x_0$ と δ 以内の距離にあるどんな x に対しても、<br/>
+
-
$f(x)$ $f(x)$ の差が ε より小さいようにすることができる。<br/>
+
-
 
+
-
関数 $f(x)I$ で連続であるとは、<br/>
+
-
$I$ に属するそれぞれの点において連続であることを言う。<br/>
+
-
 
+
-
定理 <br/>
+
-
有界閉区間上$V=[a,b]f$は、V上で可積分である。<br/>
+
-
略証;<br/>
+
-
有界閉区間上の連続関数は[[wikipedia_ja:一様連続 |一様連続]]なので、<br/>
+
-
任意の$\epsilon>0\delta>0$が存在して、<br/>
+
-
|xx|δを満たすVの任意の2点に対して、<br/>
+
-
$|f(x)-f(x')|<  \frac{\epsilon}{b-a}$<br/>
+
-
が成立する。<br/>
+
-
$V=[a,b]\Delta$を細かくして、<br/>
+
-
d(Δ)<δ<br/>
+
-
を満たすようにする。<br/>
+
-
すると、その分割によって得られた小区間Vi(i=1,2,,,n)の長さは、<br/>
+
-
全てδより小さくなるので、<br/>
+
-
$\sup \{f(x)\mid x\in V_i\}-\inf\{f(x)\mid x\in V_i\}<\frac{\epsilon}{b-a}$<br/>
+
-
M(f;Vi),m(f;Vi)の定義から<br/>
+
-
$M(f;V_i)-m(f;V_i)<\frac{\epsilon}{b-a}, (i=1,2,,,n)$
+
-
これを用いると、<br/>
+
-
$I_M(f,\Delta)-I_m(f,\Delta)=\sum_{i=1}^{n} M(f;V_i)v(V_i)-\sum_i m(f;V_i)v(V_i)$<br/>
+
-
$=\sum_i(M(f;V_i)- m(f;V_i))v(V_i)
+
-
\leq
+
-
\sum_i \frac{\epsilon}{b-a}v(V_i)
+
-
=\frac{\epsilon}{b-a}\sum_{i=1}^{n}v(V_i)
+
-
=\frac{\epsilon}{b-a}(b-a)
+
-
=\epsilon$<br/>
+
-
故に、<br/>
+
-
任意のϵ>0に対して、δ>0が存在して、<br/>
+
-
d(Δ)<δを満たす任意の分割Δにたいして、<br/>
+
-
IM(f,Δ)Im(f,Δ)ϵが示せた。<br/>
+
-
$\mathscr{S}(f)-\mathscr{s}(f)\leq
+
-
I_M(f,\Delta)-I_m(f,\Delta)$<br/>
+
-
なので<br/>
+
-
S(f)s(f) ϵ<br/>
+
-
が任意のϵ>0にたいして成立する。故に<br/>
+
-
S(f)=s(f)<br/>
+
-
可積分条件のⅲ)が示せた。証明終わり。<br/><br/>
+
-
定理の系;有限個の不連続点をもつ、有界閉区間上の関数は積分可能である。<br/>
+
====== 積分の性質 ======
====== 積分の性質 ======
1,015 行: 647 行:
従って(1)式の左辺の極限I(cf+dg) も存在して、右辺の極限と一致する。
従って(1)式の左辺の極限I(cf+dg) も存在して、右辺の極限と一致する。
証明終わり。<br/><br/>
証明終わり。<br/><br/>
-
 
======慣性モーメントの計算(1)リーマン和の極限を求める方法======
======慣性モーメントの計算(1)リーマン和の極限を求める方法======
Vは、先述の、ごく細い一様な質量密度ρ=M/lのまっすぐな棒で、<br/>
Vは、先述の、ごく細い一様な質量密度ρ=M/lのまっすぐな棒で、<br/>
1,021 行: 652 行:
原点を通りこの棒と直交する軸のまわりの(この棒の)慣性モーメントを、<br/>
原点を通りこの棒と直交する軸のまわりの(この棒の)慣性モーメントを、<br/>
リーマン和の極限を取って求めよう。<br/>
リーマン和の極限を取って求めよう。<br/>
-
区間V=[l1,ll1]をn等分して得られる点列,<br/>
+
区間V=[l1,ll1]をn(2)等分して得られる点列,<br/>
xn0=l1,xn1=l1+l/n,xni=l1+i(l/n),,,xnn=ll1<br/>
xn0=l1,xn1=l1+l/n,xni=l1+i(l/n),,,xnn=ll1<br/>
を分点とする分割をΔnと記す。すると、<br/>
を分点とする分割をΔnと記す。すると、<br/>
xnixni1=l/n,(i=1,2,,,n), d(Δn)=l/nであり、<br/>
xnixni1=l/n,(i=1,2,,,n), d(Δn)=l/nであり、<br/>
Δn={Vnj=[xnj1,xnj]j=1,2,,,n}  <br/>
Δn={Vnj=[xnj1,xnj]j=1,2,,,n}  <br/>
-
$\{{\Delta}^n \mid n=2,,,,n\}\lim_{n\to\infty} d({\Delta}^n)=\lim_{n\to\infty}\frac{l}{n}=0$を満たす。<br/>
+
$\{{\Delta}^n \mid n=2,3,,,\}$という分割の列は、<br/>
 +
limnd(Δn)=limnln=0を満たす。<br/>
y=f(x)=ρx2がリーマン可積分であることを認めれば、<br/>
y=f(x)=ρx2がリーマン可積分であることを認めれば、<br/>
可積分の定義から、どんな代表点ξnjVnjを選んでも、<br/>
可積分の定義から、どんな代表点ξnjVnjを選んでも、<br/>
-
limnIΔn(ξn1,ξn2,,,ξnn)=Iとなる。<br/>
+
$\lim_{n\to \infty}I^{f,{\Delta}^n}({{\xi}^n}_1,{{\xi}^n}_2,,,{{\xi}^n}_n)=I$となる。<br/>
そこで、代表点をξnj=xnj=l1+j(l/n)(n=2,,,),(j=1,2,,,,n)と選ぶ。<br/>
そこで、代表点をξnj=xnj=l1+j(l/n)(n=2,,,),(j=1,2,,,,n)と選ぶ。<br/>
関数y=f(x)=ρx2を用いると、
関数y=f(x)=ρx2を用いると、
分割Δnを用いた慣性モーメントの近似値は次のようになる。<br/>
分割Δnを用いた慣性モーメントの近似値は次のようになる。<br/>
-
$I^{{\Delta}^n}({x^n}_1,{x^n}_2,,,,{x^n}_n)
+
$I^{f,{\Delta}^n}({x^n}_1,{x^n}_2,,,,{x^n}_n)
=\sum_j f({x^n}_j)1/n
=\sum_j f({x^n}_j)1/n
=\sum_j f(-{l_1}+j(1/n))\frac{l}{n}
=\sum_j f(-{l_1}+j(1/n))\frac{l}{n}
1,044 行: 676 行:
=M(l12l1ln+1n+l26n+1n2n+1n)<br/>
=M(l12l1ln+1n+l26n+1n2n+1n)<br/>
故に、<br/>
故に、<br/>
-
$I=\lim_{n\to \infty}I^{{\Delta}^n}({x^n}_1,{x^n}_2,,,,{x^n}_n)
+
$I=\lim_{n\to \infty}I^{f,{\Delta}^n}({x^n}_1,{x^n}_2,,,,{x^n}_n)
  =\frac{M}{3}(l^2-3{l_1}l+3{l_1}^2)$<br/><br/>
  =\frac{M}{3}(l^2-3{l_1}l+3{l_1}^2)$<br/><br/>
   
   
1,072 行: 704 行:
[a,b]f=F(b)F(a)<br/>
[a,b]f=F(b)F(a)<br/>
上記の条件を満たす関数Fを、fの'''原始関数'''という。<br/>
上記の条件を満たす関数Fを、fの'''原始関数'''という。<br/>
-
 
(注)関数Fは、V上でしか定義されていないので、<br/>
(注)関数Fは、V上でしか定義されていないので、<br/>
端点a,bでは、通常の微分は定義できない。そこで、<br/>
端点a,bでは、通常の微分は定義できない。そこで、<br/>
1,108 行: 739 行:
I=[a,b]ρx2=13ρ(b3a3)<br/>
I=[a,b]ρx2=13ρ(b3a3)<br/>
ρ=M/l,a=l1,b=ll1を代入して、整頓すると、<br/>
ρ=M/l,a=l1,b=ll1を代入して、整頓すると、<br/>
-
=M3(l23l1l+3l12) <br/><br/>
+
=M3(l23l1l+3l12)  
-
 
+
======  重心の計算への応用 ======
-
 
+
質量密度が場所により変わる、長さlのごく細い棒Vの重心を求めてみよう。<br/>
-
======可積分条件======
+
考えやすくするため、
-
この節は、区間上で定義された関数の
+
棒の一端を原点にし、他端がx軸の正の位置にくるように座標系Oxをいれる。<br/>
-
積分可能な条件を紹介する。<br/>
+
この座標系で剛体はV=[0,l]と書ける。
-
大学の教養コース程度の数学を使うが、テキスト中で理解できるように説明する。<br/>
+
Vを小区間$V_i=[{x^n}_{i-1},{x^n}_{i}],i=1,2.\cdots,n),{x^n}_n=0,{x^n}_n=l\Delta$と記す)し、これらの小区間を質点とみなせば、その重心は、<br/>
-
興味のない方は、とばしてください。<br/>
+
$\vec{R}=\sum_i{ m_i \vec{r_i}}/M $ <br/>
 +
で定義された(1.1.1節参照)。ここでmi は第i質点の質量、$M=\sum_{i} m_i$、
 +
riは第i質点の位置ベクトル。 <br/>
 +
ベクトルを座標成分表示すると、この問題では一次元なので、
 +
$R=\sum_{i} m_i r_i/M M=\sum_{i} m_i$ <br/>
 +
しかし、実際には Viは、質点ではないので、
 +
位置ベクトルは、定まらない。<br/>
 +
またその質量も密度が一定ならば、$m_i=\rho ({x^n}_{i}- {x^n}_{i-1})$できまるが、<br/>
 +
密度が変化するならば、定まらない。<br/>
 +
そこで、各小区間 Viの代表点ξi(Vi)を選び<br/>
 +
$m_i=\rho(\xi_i)({x^n}_{i}- {x^n}_{i-1}),r_i=\xi_i $<br/>
 +
で近似する。<br/>
 +
すると分割Δと代表点{ξi}ni=1に対応する、
 +
質量Mと重心Gの近似値は、それぞれ<br/>
 +
Mρ,Δ(ξ1,,,ξn)<br/>
 +
=iρ(ξi)v(Vi)=iρ(ξi)(xixi1) <br/>
<br/>
<br/>
-
'''(1)準備;集合論の初歩'''<br/>
+
$G^{f,\Delta}(\xi_1,,,\xi_n)$<br/>
-
 
+
$=\sum_i f(\xi_i)v(V_i)=\sum_i f(\xi_i)(x_i-x_{i-1})$ <br/>
-
集合論の初歩の知識は前提にして記述するので、<br/>
+
ここで、$f(x)=\frac{1}{M}\rho(x)x$<br/>
-
なじみのない方は、下記を参考に、<br/>
+
もし、関数$\rho(x)$が積分可能ならば、分割Δを細かくしていけば
-
集合の素朴な定義、集合の和集合や共通集合、集合の包含関係などについて学習してほしい。
+
$M=\lim_{d(\Delta) \to 0}M^{\rho,\Delta}(\xi_1,,,\xi_n)$<br/>
-
*[[wikipedia_ja:集合 |ウィキペディア(集合)]]
+
$=\int_{[0,l]}\rho$<br/>
-
 
+
もし関数$\rho$の原始関数$P$が存在する($\rho(x)=\frac{dP(x)}{dx}$)ならば<br/>
-
'''(2)区間上の関数のリーマン和'''<br/>
+
$=P(l)-P(0)$<br/>
-
'''定義;リーマン和 ''' <br/>
+
もし関数f(x)も積分可能ならば、分割Δを細かくしていけば
-
区間V=[a,b]で定義され、実数に値をとる関数y=f(x)を考える。<br/>
+
$G=\lim_{d(\Delta) \to 0}G^{f,\Delta}(\xi_1,,,\xi_n)$<br/>
-
この区間の分割<br/>
+
$=\int_{[0,l]}f$<br/>
-
Δ={Vi=[xi1,xi]i=1,2,,,n},x0=a,xn=b<br/>
+
もし関数fの原始関数$F$が存在する($f(x)=\frac{dF(x)}{dx}$)ならば<br/>
-
と、その代表点ξiVi(i=1,2,,,n)に関する、y=f(x)のリーマン和とは、
+
$=F(l)-F(0)$<br/>
-
$I^{f,\Delta}(\xi_1,,,\xi_n)=$<br/>
+
例;$\rho(x)=\rho_0$ならば、$P(x)=\rho_0 x$なので  $M=\rho_0 l$<br/>
-
if(ξi)v(Vi)=if(ξi)(xixi1) <br/>
+
また、$f(x)=\frac{1}{M}\rho_{0}x$となるので$F(x)=\frac{1}{2M}\rho_{0}x^2$
-
で定義する。<br/>
+
となり、$G=F(l)-F(0)=\frac{l}{2}$<br/>
-
リーマン和は、$y=f(x)$のグラフを、棒グラフで近似したときの棒グラフの作る面積(各角柱の面積和)であることが分かる。図参照。<br/>
+
例;$\rho(x)=x$ならば、$P(x)=\frac{x^2}{2}$なので $M=\frac{l^2}{2}$<br/>
-
$y=f(x)x=ax=b$で囲まれる部分の面積を近似している。<br/>
+
このとき$f(x)=\frac{1}{M}\rho(x)x=\frac{1}{M}x^2$なので<br/>
-
 
+
$F(x)=\frac{1}{3M}l^3$である。$G=F(l)-F(0)=\frac{2l}{3}$<br/>
-
 
+
-
'''(3)リーマン和の不足リーマン和と過剰リーマン和による評価'''<br/>
+
-
Vを分割して得られた小区間$V_i=[x_{i-1},x_i]$を考える。<br/>
+
-
関数y=f(x)をこの小区間上に限定した時、<br/>
+
-
関数は、この区間上の点で最大値と最小値をとると仮定する(注参照)。<br/>
+
-
関数の最大値$max\{f(x)\mid x\in V_i\}min\{f(x)\mid x\in V_i\}m(f;V_i),M(f;V_i)$と書く。<br/>
+
-
(注) 区間上で最大値、最小値を取らない関数では、<br/>
+
-
[[wikipedia_ja:有界函数 |有界な関数]]でありさえすれば、常に存在し、しかも最大値、最小値と殆ど同じ性質をもつ、<br/>
+
-
上限、下限に置き換えれば以後の、議論は成り立つ。上限、下限については「(5)不足リーマン和の上限と過剰リーマン和の下限」で説明する。<br/><br/>
+
-
 
+
-
すると、Viの任意の点$\xi$ に対して、<br/>
+
-
$m(f;V_i)\leq f(\xi) \leq M(f;V_i)$  <br/>
+
-
故に、<br/>
+
-
'''補題1'''<br/>
+
-
ⅰ)どのような代表点$\{\xi_i\}_{i}, (\xi_i \in V_i,i=1,2,,,n)$に対しても<br/>
+
-
$I_{m}(f,\Delta):=I^{f,\Delta}(\xi_{1}^m,,,\xi_{n}^m)
+
-
=\sum_i m(f;V_i)v(V_i)$<br/>
+
-
$\leq
+
-
I^{f,\Delta}(\xi_1,,,\xi_n)=\sum_i f(\xi_i)v(V_i)$<br/>
+
-
$\leq
+
-
\sum_i M(f;V_i)v(V_i)
+
-
=I_{M}(f,\Delta)=I^{f,\Delta}(\xi_{1}^M,,,\xi_{n}^M)       \qquad  (1)$ <br/>
+
-
そこで、$I_{m}(f,\Delta)\Delta)fI_{M}(\Delta)$を'''過剰リーマン和'''と呼ぶ。<br/>
+
-
ⅱ)$I_{m}(f,\Delta)=\min_{\xi_i \in V_i,i=1,2,,,n}I^{f,\Delta}(\xi_{1}^M,,,\xi_{n}^M) <br/>
+
-
$I_{M}(f,\Delta)=\max_{\xi_i \in V_i,i=1,2,,,n}I^{f,\Delta}(\xi_{1}^M,,,\xi_{n}^M) <br/>
+
-
証明は明らかなので省略。<br/>
+
-
 
+
-
'''(4)分割の細分とリーマン和の評価式'''<br/>
+
-
'''定義;分割の細分'''<br/>
+
-
Vの分割Δが分割$\Delta$の細分というのは、<br/>
+
-
$\Delta\{x_0,x_1,,,,x_n\}$が、<br/>
+
-
Δの分点の集合{x0,x1,,,,xn}に真に含まれることと定義する。<br/>
+
-
記号でかけば、$\{x_0,x_1,,,,x_n\}\subset \{x'_0,x'_1,,,,x'_{n'}\},
+
-
\{x_0,x_1,,,,x_n\}\neq \{x'_0,x'_1,,,,x'_{n'}\}$。
+
-
記号では、ΔΔと記す。<br/><br/>
+
-
 
+
-
 
+
-
'''補題2'''<br/>
+
-
Δ Δという分割に対し、<br/>
+
-
$I_{m}(f,\Delta)
+
-
\leq I_{m}(f,\Delta')
+
-
\leq I_{M}(f,\Delta')
+
-
\leq I_{M}(f,\Delta)      \qquad  (2)$
+
-
<br/>
+
-
が成り立つ。    <br/>
+
-
(証明)<br/>
+
-
Δの小区間$V_i=[x_{i-1},x_i]{\Delta}'$では、<br/>
+
-
$\{V'_j=[x_{i-1},x'_j],V'_{j+1}=[x'_j,x_i]\}$の2つに分割されたとする。<br/><br/>
+
-
すると、区間上の関数の最大値と最小値の定義から、<br/>
+
-
$m(f;V_i)  \leq  m(f;V'_j)$ $\quad m(f;V_i)  \leq  m(f;V'_{j+1})$<br/>
+
-
$M(f;V_i)  \geq  M(f;V'_j)$  $\quad M(f;V_i)  \geq  M(f;V'_{j+1})$<br/>
+
-
これらから、命題は成立することが分かる。<br/><br/>
+
-
 
+
-
'''(5)不足リーマン和の上限と過剰リーマン和の下限'''
+
-
補題2から、分割の細分を繰り返していくと、その分割に対応する、<br/>
+
-
不足リーマン和は、広義増加(増加するか、同じ値にとどまる)し、<br/>
+
-
過剰リーマン和は、広義減少する。<br/>
+
-
分割を細かくしていったとき、これらの極限が一致すれば、補題1から、<br/>
+
-
リーマン和の極限値は、代表点に無関係に、定まることになる。<br/>
+
-
 
+
-
そこで色々な分割に対応する不足リーマン和のなかの最大値と過剰リーマン和の最小値を求めることが、重要になる。しかし一般にはこれらは存在しないことが示せる。<br/>
+
-
そこで最大値に近い性質を持つ上限と最小値に近い下限という概念を利用する。<br/>
+
-
'''定義;上界と下界'''<br/>
+
-
Rを、全ての実数を要素とする集合とし、<br/>
+
-
Aをその部分集合とする。<br/>
+
-
実数uAの上界(upper bound)とは、<br/>
+
-
任意のaAに対して、auがなりたつこと。<br/>
+
-
実数lAの下界(lower bound)とは、<br/>
+
-
任意のaAに対して、laがなりたつこと。<br/>
+
-
UAAの上界をすべて集めた集合、<br/>
+
-
LAAの上界をすべて集めた集合とする。<br/>
+
-
UAが空集合でない(すなわち、Aの上界が少なくとも一つ存在する)とき、<br/>
+
-
Aは'''上に有界'''であるといい、<br/>
+
-
LAの時、Aは'''下に有界'''であるという。<br/>
+
-
上に有界で、下にも有界な集合(R)は、'''有界'''という。<br/>
+
-
 
+
-
'''実数の連続の公理'''<br/>
+
-
以下の性質は、色々な極限の存在の根拠を与えるもので、実数の持つ最も重要な性質の一つである。<br/>
+
-
ARとする。<br/>
+
-
もし、UAならば、UAは、最小元を持つ。<br/>
+
-
これをAの'''上限(supremum)'''あるいは'''最小上界(least upper bound)'''という。<br/>
+
-
もし、LAならば、LAは、最大元を持つ。<br/>
+
-
これをAの'''下限(infimum)'''あるいは'''最大下界(greatest lower bound)'''というという。<br/><br/>
+
-
 
+
-
補題3<br/>
+
-
uA(R) の上限となるための必要十分条件は、<br/>
+
-
ⅰ)uAの上界。すなわち任意のaAにたいしてau   <br/>
+
-
ⅱ)$x<uxAx<aa\in A$が存在。<br/>
+
-
ⅲ)Aが最大値を持つ場合には、上限は最大値と一致する。<br/>
+
-
同様に、lA の下限となるための必要十分条件は、<br/>
+
-
ⅰ)lAの下界。すなわち任意のaAにたいしてla   <br/>
+
-
ⅱ)l<xである任意のxAの下界ではない。すなわち、a<xとなるaAが存在。<br/>
+
-
ⅲ)Aが最小値を持つ場合には、下限は最小値と一致する。<br/><br/>
+
-
A の上限をsupA、下限をinfAと書く。<br/><br/>
+
-
証明は、上限、下限の定義から、明らかなので省略する。<br/>
+
-
例;$A=(0,1)\sup A=1,\inf A=0$。<br/>
+
-
これらは、ともにAの要素でないので、<br/>
+
-
上限1はAの最大元(最大値)ではなく、下限0はAの最小元(最小値)ではない。<br/>
+
-
A=[0,1]のとき、supA=1,infA=0。<br/>
+
-
これらは、ともにAの要素なので、<br/>
+
-
上限は最大限であり、下限は最小限となる。<br/>
+
-
 
+
-
補題4.
+
-
$A \subset B \subset {\bf R}B$は有界集合とする。<br/>
+
-
このとき、infBinfAsupAsupB<br/>
+
-
証明は容易である。<br/><br/>
+
-
 
+
-
関数$y=f(x)$が連続でない時は、区間上で最大値や最小値を取らないことがある。
+
-
この場合も考慮して、最大値を上限に、最小値を下限に置き換えて、$m(f;V_i)=\inf\{f(x)\mid x\in V_i\},M(f;V_i)=\sup \{f(x)\mid x\in V_i\}$で定義すれば、
+
-
有界関数に対して、これらは常に定義され、今までの議論はすべて成り立つ。
+
-
 
+
-
 
+
-
''' 2つの分割の共通の細分 '''<br/>
+
-
分割$\Delta\{x_j \mid j=1,2,,,m\}$と、<br/>
+
-
分割Δ の分点の集合{xjj=1,2,,,n}の<br/>
+
-
和集合{xjj=1,2,,,m}{xjj=1,2,,,n}を分点とする分割をΔΔと書く。<br/>
+
-
すると新しい分割は<br/>
+
-
ΔΔΔ と
+
-
ΔΔΔ<br/>
+
-
を満たす。<br/>
+
-
これを用いると、不足リーマン和の上限$\mathscr{s}(f)\mathscr{S}(f)$が存在することが証明できる。<br/>
+
-
 
+
-
補題5<br/>
+
-
$fV=[a,b](\{f(x)\mid x\in V\}$Rの有界集合)とする。<br/>
+
-
V=[a,b]の分割を全て集めて作った集合をD(V)と書く。<br/>
+
-
すると、<br/>
+
-
ⅰ)任意の$\Delta,{\Delta}'\in \mathscr{D}(V)$に対して、
+
-
$I_m(f,\Delta) \leq I_M(f,{\Delta)}')$<br/>
+
-
ⅱ)集合{Im(f,Δ)ΔD(V)}は上に有界、<br/>
+
-
集合{IM(f,Δ)ΔD(V)}は下に有界<br/>
+
-
ⅲ)$\mathscr{s}(f):=\sup\{I_m(f,\Delta) \mid \Delta \in \mathscr{D}(V)\}$と
+
-
<br/>
+
-
S(f):=inf{IM(f,Δ)ΔD(V)}は存在し、<br/>
+
-
s(f) S(f) <br/>
+
-
証明;<br/>
+
-
ⅰ)ΔΔΔ なので、補題2から、<br/>
+
-
$I_m(f,\Delta)
+
-
\leq I_m(f,\Delta\vee {\Delta}') 
+
-
\leq I_M(f,\Delta\vee {\Delta}') 
+
-
\leq I_M(f,{\Delta}') $<br/>
+
-
ⅱ)1)で証明した不等式で、分割Δ は固定する。<br/>
+
-
すると全ての分割 Δに対して、Im(f,Δ)IM(f,Δ))なので<br/>
+
-
集合{Im(f,Δ)ΔD(V)}は、上界IM(f,Δ))を持ち、上に有界である。<br/>
+
-
後者も同様にして下に有界であることが示せる。<br/>
+
-
ⅲ)従って、実数の連続性の公理から、<br/>
+
-
集合{Im(f,Δ)ΔD(V)}は上限s(f)をもち、<br/>
+
-
集合{IM(f,Δ)ΔD(V)}は下限S(f)をもつ。<br/>
+
-
上限は、上界の中の最小値なので、<br/>
+
-
s(f)IM(f,Δ)<br/>
+
-
この式は任意のΔについて成立するので、<br/>
+
-
s(f)は、集合{IM(f,Δ)ΔD(V)}の下界である。<br/>
+
-
下限S(f)は、下界のなかの最大値なのでs(f)S(f)を得る。<br/>
+
-
 
+
-
 
+
-
 
+
-
'''(6)分割を細かくしていくと不足リーマン和はその上限sに、過剰リーマン和はその下限Sに収束する'''<br/>
+
-
この命題を正確に述べるには、
+
-
まず、分割Δ={Vi=[xi1,xi]i=1,2,,,n}の大きさを、きちんと定める必要がある。<br/>
+
-
定義;Δ={Vi=[xi1,xi]i=1,2,,,n}の大きさd(Δ)とは、<br/>
+
-
この分割で得られた小区間の長さの、最大値で定義する。記号で書くと<br/>
+
-
d(Δ)=max{xixi1i=1,2,,,n}
+
-
 
+
-
定理(ダルブー;Darboux)<br/>
+
-
V=[a,b]<br/>
+
-
fを、Vで定義され、実数に値を取る有界関数とする。<br/>
+
-
このとき、<br/>
+
-
ⅰ)$\lim_{d(\Delta) \to 0}I_m(f,\Delta)=\mathscr{s}(f)$<br/>
+
-
ⅱ)$\lim_{d(\Delta) \to 0}I_M(f,\Delta)=\mathscr{S}(f)$<br/>
+
-
証明;<br/>
+
-
ⅰ)を示す。ⅱ)は同じようにして証明できるので略す。<br/>
+
-
これを示すには、どんなに小さい正の実数$\epsilon\delta_{\epsilon}\delta_{\epsilon}\Delta\mathscr{s}(f)-I_m(f,\Delta)<\epsilon$であることを示せばよい。以下に、数段階に分けて、これを証明する。<br/>
+
-
 
+
-
$\quad 1)<br/>
+
-
$上限の性質(補題3)から、<br/>
+
-
ある分割<br/>
+
-
$D=\{{V^D}_i=[{x^D}_{i-1},{x^D}_i] \mid i=1,2,,,n\}in \mathscr{D}(V)$<br/>
+
-
が存在して、<br/>
+
-
$\mathscr{s}(f)-I_m(f,D)<\frac{\epsilon}{2} \qquad (1)$<br/>
+
-
今後このDを使って、証明を進める。<br/>
+
-
 
+
-
2)<br/>
+
-
分割Dの小区間VDiの長さ(xDixDi1)(i=1,2,,,n)
+
-
最小値をeとおくと<br/>
+
-
e=minni=1(xDixDi1) <br/>
+
-
eに比べて非常に小さい大きさを持つ分割、<br/>
+
-
$\Delta=\{V^{\Delta}_i=[{x^{\Delta}}_{i-1},{x^{\Delta}}_i] \mid i=1,2,,,N\}$、
+
-
<br/>
+
-
d(Δ)=maxi=1,2,,,N(xΔixΔi1)e<br/>
+
-
<br/>
+
-
を考える。<br/>
+
-
もし、DΔならば補題2より、$I_m(f,D) \leq I_m(f,\Delta)$<br/>
+
-
すると$\mathscr{s}(f)-I_m(f,\Delta)\leq \mathscr{s}(f)-I_m(f,D) \leq
+
-
\frac{\epsilon}{2}\leq \epsilon$ <br/>
+
-
通常、分割$\Delta$は、Dの細分になっていない。<br/>
+
-
この場合は、いくつか(高々n-1個)の$\DeltaD$の小区間には含まれず、<br/>
+
-
Dの分点${x^D}_i(i=1,2,,,n-1)$をまたぐことになる。図参照のこと。<br/>
+
-
議論を簡単にするため、Dの分点${x^D}_i(i=1,2,,,n-1)\Delta$の小区間によって跨がれていると仮定し、議論を進める。<br/>
+
-
他のケースでも、証明はおなじようにできるので、このように仮定しても何の問題も起こらない。<br/>
+
-
Dの分点xDiを跨ぐΔの小区間をVΔmiとする(i=1,2,,,n-1)。<br/>
+
-
3)  <br/>
+
-
2つの分割DΔからΔ:=DΔを作る。<br/>
+
-
すると<br/>
+
-
Δ={VΔ1,VΔ2,,,,,,,,,VΔm11,<br/>
+
-
[xΔm11,xD1],[xD1,xΔm1],<br/>
+
-
VΔm1+1,VΔm1+2,,,,,,,,,VΔm21,<br/>
+
-
[xΔm21,xD2],[xD2,xΔm2],<br/>
+
-
VΔm2+1,VΔm2+2,,,VΔm31,<br/><br/>
+
-
,,,,,,,,,<br/><br/>
+
-
VΔmn1+1,VΔmn1+2,,,,,,,,,VΔN} (2)<br/>
+
-
と書ける。<br/>
+
-
 
+
-
ΔΔで、 DΔ なので、<br/>
+
-
Im(f,Δ)Im(f,Δ), Im(f,D)Im(f,Δ)<br/>
+
-
後者の式から、<br/>
+
-
0s(f)Im(f,Δ)s(f)Im(f,D)<br/>
+
-
この式と(1)式から、<br/>
+
-
$0 \leq \mathscr{s}(f)-I_m(f,{\Delta}')<\frac{\epsilon}{2}$<br/>
+
-
そこで、<br/>
+
-
「$d(\Delta) \to 0 I_m(f,{\Delta}')-I_m(f,\Delta)<\frac{\epsilon}{2}$
+
-
<br/>
+
-
が示せれば、<br/>
+
-
$0 \leq \mathscr{s}(f)-I_m(f,\Delta)$<br/>
+
-
$=(\mathscr{s}(f)-I_m(f,{\Delta}')+(I_m(f,{\Delta}'-I_m(f,\Delta)\leq \epsilon$<br/>
+
-
が示され、証明が終わる。<br/>
+
-
$\quad 4)$  <br/>
+
-
$I_{m}(f,\Delta)=\sum_{i=1}^{N} m(f;V^{\Delta}_i)v(V^{\Delta}_i)$
+
-
であり、<br/>
+
-
(2)式から、<br/>
+
-
$I_m(f,{\Delta}')$<br/>
+
-
$=\sum_{i\notin \{m_1,m_2,,,,m_{n-1}\}} m(f;V^{\Delta}_i)v(V^{\Delta}_i)$<br/>
+
-
$+\sum_{k=1}^{n-1} m(f;[x^{\Delta}_{m_{k}-1},x^{D}_k])v([x^{\Delta}_{m_{k}-1},x^{D}_k])$<br/>
+
-
$+\sum_{k=1}^{n-1} m(f;[x^{D}_k,x^{\Delta}_{m_{k}}])v([x^{D}_k,x^{\Delta}_{m_{k}}])$<br/>
+
-
なので、<br/>
+
-
Im(f,Δ)Im(f,Δ)<br/>
+
-
$=\sum_{k=1}^{n-1} m(f;[x^{\Delta}_{m_{k}-1},x^{D}_k])v([x^{\Delta}_{m_{k}-1},x^{D}_k])$<br/>
+
-
$+\sum_{k=1}^{n-1} m(f;[x^{D}_k,x^{\Delta}_{m_{k}}])v([x^{D}_k,x^{\Delta}_{m_{k}}])$<br/>
+
-
i{m1,m2,,,,mn1}m(f;VΔi)v(VΔi)<br/>
+
-
関数はV上で有界なので、適切に正の実数Mを選ぶと、xVの要素ならば<br/>
+
-
|f(x)|Mが成立する。<br/>
+
-
すると$|m(f;[x^{\Delta}_{m_{k}-1},x^{D}_k])|, |m(f;[x^{D}_k,x^{\Delta}_{m_{k}}])|
+
-
\leq M$<br/>
+
-
が成り立つ。また<br/>
+
-
$v(V^{\Delta}_{m_k})
+
-
=v([x^{\Delta}_{m_{k}-1},x^{D}_k])+v([x^{D}_k,x^{\Delta}_{m_{k}}])$で、<br/>
+
-
v(VΔi)d(Δ)<br/>
+
-
なので<br/>
+
-
$|I_m(f,{\Delta}')-I_m(f,\Delta)|\leq 2M\sum_{i\in \{m_1,m_2,,,,m_{n-1}\}} v(V^{\Delta}_i)\leq 2M(n-1)d(\Delta)$<br/>
+
-
そこで、<br/>
+
-
δϵ=ϵ4Mn
+
-
と選べば、<br/>
+
-
d(Δ)δϵをみたすどのような分割Δも、<br/>
+
-
$0\leq I_m(f,{\Delta}')-I_m(f,\Delta)|\leq \frac{\epsilon}{2}$<br/>
+
-
を満たすことが証明できた。証明終わり。<br/>
+
-
 
+
-
 
+
-
  '''(7)可積分条件'''<br/>
+
-
定理;可積分条件 <br/>
+
-
$V=[a,b]$<br/>
+
-
fを、Vで定義され、実数に値を取る有界関数とする。<br/>
+
-
次の条件のうち1つが成立すれば、残り2つは成立する(互いに同値という)。
+
-
ⅰ)fV上で(リーマン)可積分<br/>
+
-
ⅱ)$\lim_{d(\Delta) \to 0}(I_M(f,\Delta)-I_m(f,\Delta))=0$<br/>
+
-
ⅲ)$\mathscr{S}(f)=\mathscr{s}(f)$<br/>
+
-
<br/>
+
-
証明<br/>
+
-
ⅰ)を仮定する。ⅱ)が成立することを示そう。<br/>
+
-
fの積分値をαとおくと、可積分の定義から、<br/>
+
-
任意のϵ>0に対して、δ>0が存在して、<br/>
+
-
d(Δ)<δである任意の分割と、その分割の任意の代表点ξi,(i=1,2,,,)に対し,<br/>
+
-
$|I^{f,\Delta}(\xi_1,,,\xi_n)-\alpha |<\frac{1}{2}\epsilon$<br/>
+
-
が成立する。<br/>
+
-
変形すると<br/>
+
-
$\alpha-\frac{1}{2}\epsilon
+
-
<I^{f,\Delta}(\xi_1,,,\xi_n)
+
-
<\alpha+\frac{1}{2}\epsilon \qquad (1)  $<br/>
+
-
ここで、補題1のⅱ)から、<br/>
+
-
$\inf_{\{\xi_i\}}I^{f,\Delta}(\xi_1,,,\xi_n) =I_{m}(f,\Delta)$<br/>
+
-
$\sup_{\{\xi_i\}}I^{f,\Delta}(\xi_1,,,\xi_n) =I_{M}(f,\Delta)$<br/>
+
-
なので、<br/>
+
-
(1)式から、<br/>
+
-
$\alpha-\frac{1}{2}\epsilon
+
-
\leq
+
-
I_{m}(f,\Delta)
+
-
\leq
+
-
I_{M}(f,\Delta)
+
-
\leq
+
-
\alpha+\frac{1}{2}\epsilon$<br/>
+
-
これより、任意のϵ>0に対して、δ>0が存在して、<br/>
+
-
$d(\Delta)<\delta \implies (0\leq I_{M}(f,\Delta)-I_{m}(f,\Delta)\leq \epsilon)$<br/>
+
-
ⅱ)が示せた。<br/>
+
-
ⅱ)を仮定する。 ⅲ)が成り立つことを示す。<br/>
+
-
 
+
-
$I_{m}(f,\Delta)
+
-
\leq
+
-
\mathscr{s}(f):=\sup_{\Delta}I_{m}(f,\Delta)
+
-
\leq
+
-
\mathscr{S}(f):=\inf_{\Delta}I_{M}(f,\Delta)
+
-
\leq
+
-
I_{M}(f,\Delta)$<br/>
+
-
なので、<br/>
+
-
$0
+
-
\leq
+
-
\mathscr{S}(f)-\mathscr{s}(f)
+
-
\leq
+
-
I_{M}(f,\Delta)-I_{m}(f,\Delta)$<br/>
+
-
故に、分割を細かくしていき、極限をとると、<br/>
+
-
$0
+
-
\leq
+
-
\mathscr{S}(f)-\mathscr{s}(f)
+
-
\leq
+
-
\lim_{d(\Delta)\to 0}(I_{M}(f,\Delta)-I_{m}(f,\Delta))$<br/>
+
-
ⅱ)が成立するので、<br/>
+
-
$=0$<br/>
+
-
ⅲ)が示せた。<br/>
+
-
ⅲ)を仮定する。 $\alpha=\mathscr{S}(f)=\mathscr{s}(f)$とおく。<br/>
+
-
ⅰ)が成り立つことを示そう。<br/>
+
-
補題1のⅰ)から、<br/>
+
-
どのような分割Δと、<br/>
+
-
その代表点$\{\xi_i\}_{i} (\xi_i \in V_i,i=1,2,,,n)$に対しても<br/>
+
-
$I_{m}(f,\Delta)
+
-
\leq
+
-
I^{f,\Delta}(\xi_1,,,\xi_n)
+
-
\leq
+
-
I_{M}(f,\Delta)$<br/>
+
-
ここで、ダルブーの定理から、<br/>
+
-
$\lim_{d(\Delta) \to 0}I_{m}(f,\Delta)=\mathscr{s}(f)=\alpha$,<br/>
+
-
$\lim_{d(\Delta) \to 0}I_{M}(f,\Delta)=\mathscr{S}(f)=\alpha$<br/>
+
-
が成り立つので、<br/>
+
-
limd(Δ)0If,Δ(ξ1,,,ξn)=α <br/>
+
-
が成り立つ。<br/>
+
-
ⅰ)が示せた。<br/>
+
-
 
+
-
  '''(8)有限個の点を除いて連続な閉区間上の関数は積分可能'''<br/>
+
-
色々な関数のグラフを書くとつながっているところを、跳んでいるところが出来る。<br/>
+
-
y=Xのグラフはずっとつながっている。<br/>
+
-
x<0のとき  y=0, 0xのとき  y=1 で関数を定義すると、x=0のところで
+
-
そのグラフは跳んでいる。<br/>
+
-
連続や不連続は関数の非常に重要な性質であり、それを調べることはとても豊かな知識をもたらす。<br/>
+
-
しかし正確に議論するには、連続とは何かをきちんと定義する必要がある。<br/>
+
-
関数の連続性の定義;<br/>
+
-
 
+
-
実数値関数 $f(x)x_0$で連続'''であるとは、<br/>
+
-
$xx_0f(x)f(x_0)\lim_{x\to x_0} f(x) = f(x_0)$と記す。<br/>
+
-
 
+
-
これはイプシロン-デルタ論法(ε-δ論法)を用いれば次のように定式化できる。<br/>
+
-
(小さな)正の数 ε が任意に与えられたとき、<br/>
+
-
(小さな)正の数 δ をうまくとってやれば、<br/>
+
-
x0 と δ 以内の距離にあるどんな x に対しても、<br/>
+
-
f(x)f(x) の差が ε より小さいようにすることができる。<br/>
+
-
 
+
-
関数 f(x) がある区間I で連続であるとは、<br/>
+
-
I に属するそれぞれの点において連続であることを言う。<br/>
+
-
 
+
-
定理 <br/>
+
-
有界閉区間上V=[a,b]で定義され、実数に値を取る連続関数fは、V上で可積分である。<br/>
+
-
略証;<br/>
+
-
有界閉区間上の連続関数は[[wikipedia_ja:一様連続 |一様連続]]なので、<br/>
+
-
任意のϵ>0に対して、δ>0が存在して、<br/>
+
-
|xx|δを満たすVの任意の2点に対して、<br/>
+
-
$|f(x)-f(x')|<  \frac{\epsilon}{b-a}$<br/>
+
-
が成立する。<br/>
+
-
V=[a,b]の分割Δを細かくして、<br/>
+
-
d(Δ)<δ<br/>
+
-
を満たすようにする。<br/>
+
-
すると、その分割によって得られた小区間$V_i(i=1,2,,,n)$の長さは、<br/>
+
-
全てδより小さくなるので、<br/>
+
-
$\sup \{f(x)\mid x\in V_i\}-\inf\{f(x)\mid x\in V_i\}<\frac{\epsilon}{b-a}$<br/>
+
-
$M(f;V_i),m(f;V_i)$の定義から<br/>
+
-
$M(f;V_i)-m(f;V_i)<\frac{\epsilon}{b-a}, (i=1,2,,,n)$
+
-
これを用いると、<br/>
+
-
$I_M(f,\Delta)-I_m(f,\Delta)=\sum_{i=1}^{n} M(f;V_i)v(V_i)-\sum_i m(f;V_i)v(V_i)$<br/>
+
-
$=\sum_i(M(f;V_i)- m(f;V_i))v(V_i)
+
-
\leq
+
-
\sum_i \frac{\epsilon}{b-a}v(V_i)
+
-
=\frac{\epsilon}{b-a}\sum_{i=1}^{n}v(V_i)
+
-
=\frac{\epsilon}{b-a}(b-a)
+
-
=\epsilon$<br/>
+
-
故に、<br/>
+
-
任意のϵ>0に対して、δ>0が存在して、<br/>
+
-
d(Δ)<δを満たす任意の分割Δにたいして、<br/>
+
-
IM(f,Δ)Im(f,Δ)ϵが示せた。<br/>
+
-
$\mathscr{S}(f)-\mathscr{s}(f)\leq
+
-
I_M(f,\Delta)-I_m(f,\Delta)$<br/>
+
-
なので<br/>
+
-
S(f)s(f) ϵ<br/>
+
-
が任意のϵ>0にたいして成立する。故に<br/>
+
-
S(f)=s(f)<br/>
+
-
可積分条件のⅲ)が示せた。証明終わり。<br/><br/>
+
-
定理の系;有限個の不連続点をもつ、有界閉区間上の関数は積分可能である。<br/>
+
 +
===== 2次元以上の物体の慣性モーメントについて =====
====  てこの原理と力のモーメント====
====  てこの原理と力のモーメント====
図のように剛体の棒の中間に支点Oがあり、<br/>
図のように剛体の棒の中間に支点Oがあり、<br/>
1,556 行: 814 行:
も参照のこと。
も参照のこと。
-
====  作用線の定理====
+
====  剛体に働く力の作用線====
-
剛体の場合、作用線に沿って力の作用点を移動しても、力の作用は変わらない。何故かは、考えてみましょう。<br/>
+
力が作用する点を着力点といい、<br/>
-
*[[wikipedia_ja:作用線の定理|ウィキペディア(作用線の定理)]]
+
着力点を通り力のベクトルと方向が等しい直線を、力の作用線という。<br/>
 +
剛体に働く力は、その着力点をかえると、一般には、剛体の運動への効果が異なってしまう。<br/>
 +
しかし、力のベクトル和と、力のモーメント和が不変となるように力の着力点を移動したり力の合成をすることは、<br/>
 +
剛体の運動には全く影響がでないので、許される。<br/>
 +
例えば、力の着力点をその作用線にそってうごかしたり、<br/>
 +
同じ着力点をもつ複数の力を、それらのベクトル和に置き換えることは許される。。
 +
<br/>
-
=== 剛体のつり合い===
+
===剛体のつり合い====  
いくつかの力が作用し、剛体が静止したままであるか、<br/>
いくつかの力が作用し、剛体が静止したままであるか、<br/>
-
重心が等速直線運動を続け、重心の周りの回転が変化しない場合に、剛体(に作用している力)は釣り合っているという。
+
重心Gが等速直線運動(静止も含む)を続け、<br/>
-
<br/>
+
重心の周りの回転が変化しない(回転しないままか、同じ回転を続ける)場合に、<br/>
-
==  気体や液体の圧力と浮力==
+
剛体(に作用している力)は釣り合っているという。<br/>
 +
重心が等速直線運動を行うのは、<br/>
 +
剛体に作用する外力のベクトル和が0になることであり、その場合に限る。<br/>
 +
これについては、「1.1.1  質点系の運動と重心」で説明した。<br/>
 +
重心周りの回転が変化しないのは、重心まわりの外力のモーメントの総和が0になることであり、この場合に限る。これについては、「1.2.3.5 固定軸の周りの剛体の回転運動の方程式」で説明した。<br/>
 +
定理;剛体のつり合い<br/>
 +
剛体に、外力F1,F2,,,,Fnがはたらいている。<br/>
 +
このとき、次の条件は同等である。<br/>
 +
ⅰ)剛体は釣り合っている。<br/>
 +
ⅱ)外力のベクトル和が零で、重心Gまわりの外力のモーメントの和が零。<br/>
 +
ⅲ)外力のベクトル和が零で、任意の固定点Pまわりの外力のモーメントの和が零。<br/>
 +
証明;条件ⅰ)とⅱ)が同等であることは、すでに、説明した。<br/>
 +
条件ⅱ)とⅲ)の同等性を示そう。<br/>
 +
外力の和が零であるという条件の下で、<br/>
 +
「任意の固定点Pまわりの外力のモーメントの和NPは常に等しい」<br/>
 +
ことを示せば良い。<br/>
 +
外力Fiの作用点をPi(i=1,2,n)とする。<br/>
 +
すると、<br/>
 +
Pまわりの外力のモーメントの和NPは<br/>
 +
NP=ni=1PPi×Fi(1)
 +
 
 +
任意の点Qまわりの外力のモーメントの和NQは<br/>
 +
NQ=ni=1QPi×Fi(2)<br/>
 +
 
 +
PPi=PQ+QPiを(1)式に代入すると<br/>
 +
$\vec N_{P}=\sum_{i=1}^{n}\vec{PP_i}\times \vec{F^i}
 +
=\sum_{i=1}^{n}(\vec{PQ}+\vec{QP_i})\times \vec{F^i}$<br/>
 +
ベクトル積の性質から、<br/>
 +
$=\sum_{i=1}^{n}(\vec{PQ}\times \vec{F^i}+\vec{QP_i}\times \vec{F^i})
 +
=\vec{PQ}\times \sum_{i=1}^{n}\vec{F^i}+\sum_{i=1}^{n}\vec{QP_i}\times \vec{F^i}$<br/>
 +
仮定と(2)式から、<br/>
 +
=PQ×0+NQ=NQ<br/>
 +
故に、NP=NQ<br/>
 +
証明終わり。  
 +
 
 +
==  気体と液体の圧力 ==
 +
この節では気体や液体を、<br/>
 +
分子や原子という粒子から構成されるという微視的立場でなく、<br/>
 +
巨視的に捉え空間的に滑らかな連続体であるとみなす。<br/>
 +
連続体の内部の微小部分に働く力を考え、其の釣合いについて考え、<br/>
 +
圧力の性質を導く。
 +
 
 +
=== 気体や液体とは何か。===
 +
*[[wikipedia_ja:気体 |ウィキペディア(気体)]]
 +
*[[wikipedia_ja:液体 |ウィキペディア(液体)]]
 +
=== 気体と液体の特徴===
 +
気体と液体は体積の変化には抵抗するが、<br/>
 +
形の変化には、抵抗しない。(ただし非常に速い変化には抵抗する)。<br/>
 +
但し、気体の体積変化への抵抗は小さく、液体は非常に大きい。<br/>
 +
=== 静止気体と液体の圧力===
 +
気体や液体は、その表面または内部に任意の面を考えると、その面で2分される部分は、<br/>
 +
互いに他を押している。それらは大きさ・方向は等しく、逆向きである(作用反作用の法則)<br/>
 +
単位面積当たりのこの力を'''応力'''とよぶ。<br/>
 +
その発生は、重力の存在と前述の気体や液体の特徴(形の変化に抵抗しない)に起因する。<br/>
 +
この力の性質を、気体・液体の特徴から導こう。
 +
====応力は面に垂直に働く ====
 +
説明は便宜上、液体の語で述べる。<br/>
 +
命題1:<br/>
 +
静止した液体(気体)の表面あるいは内部に任意のなめらかな面(注参照)を考える。<br/>
 +
この面上の応力は、常にこの面に直角に働く。<br/>
 +
面と常に直角に働く応力を、'''圧力'''と呼ぶ。<br/>
 +
(注)面のどの一点においても、その点にごく近い面の部分だけをみれば、平面とみなせる曲面のこと。
 +
理由;<br/>
 +
もし、ある面上のある一点Pの周辺の微小面部分(Sと書く)で、押し合う力がこの面と平行な成分を持つとする。<br/>
 +
Sは仮定より、平面(の一部)と考えてよい。<br/>
 +
 
 +
図のように、面部分Sとそれと平行な平面の一部S’から作られる、<br/>
 +
非常に薄い液体の板状部分Vを考える。<br/>
 +
 
 +
するとVがSを通して液体から受ける力の総和FSは、面Sと平行な成分をもつ。<br/>
 +
面SとS’は、非常に近いので、<br/>
 +
Sを挟んで押し合う力と、S’を挟んで押し合う力は、単位面積当たり、ほぼ等しいと考えてよい。<br/>
 +
すると、VがS’を通して液体から受ける力FSは、<br/>
 +
Sを通して受ける力と大きさと方向はほぼ同じで、逆向きになる。<br/>
 +
FSの面Sと平行な成分も、FSのSと平行な成分と大きさはおなじで、逆向きになる。<br/>
 +
液体は自由に形を変えられるので、VのS面とS’面は逆方向に動いてしまい、<br/>
 +
静水という条件に反してしまう。<br/>
 +
従って、<br/>
 +
「ある面上のある一点Pの周辺の微小面部分Sで、押し合う力がこの面と平行な成分を持つ」<br/>
 +
という仮定はあり得ないことが示された。<br/><br/>
 +
 
 +
命題2<br/>
 +
どの面にも直角に働く応力(圧力)は、どの点でも面の方向によらず一定の強さ(大きさ)をもつ。<br/>
 +
証明;<br/>
 +
液体中の任意の点をOとする。<br/>
 +
Oを原点とする、直交右手系Oxyzを定める。<br/>
 +
Oを通る任意の面Hをとる。<br/>
 +
O点における、<br/>
 +
この面における圧力pとxy平面における圧力pz、yz平面、zx平面における圧力px,py<br/>
 +
が等しいことを示そう。<br/>
 +
平面Hと平行でO点の近くを通る平面Hが<br/>
 +
x軸、y軸、z軸と交わる点をそれぞれ、<br/>
 +
A(αa,0,0),B(0,αb,0),C(0,0,αc)とおく。図参照。<br/>
 +
四面体OABCの外部の液体が、<br/>
 +
OBCを押す力をFx,OCAを押す力をFy,OABを押す力をFz,ACBを押す力をF<br/>
 +
とおく。<br/>
 +
四面体内の液体が静止しているので、<br/>
 +
Fx+Fy+Fz+F=0(1)<br/>
 +
が成り立つ。<br/>
 +
この式を圧力で表示しよう。<br/>
 +
limα0Fx|OBC|=pxなので、<br/>
 +
αが十分小さければ<br/>
 +
Fx=|OBC|px=12|αbαc|px<br/>
 +
故に、2Fx=OB×OCpx=αbαcpxex<br/>
 +
同様に2Fy=OC×OApy=αcαapyey、<br/>
 +
2vecFz=OA×OBpz=αaαbpzez<br/>
 +
2F=AC×ABp=p(αa,0,αc)×(αa,αb,0)<br/>
 +
これらを(1)式に代入して<br/>
 +
pxex+pyey+pzez+AC×ABp=0(2)<br/>
 +
これを計算すると、<br/>
 +
(α2bc(pxp),α2ca(pyp),α2ab(pzp))=0<br/>
 +
これより、p=px=py=pz   証明終わり。<br/>
 +
 
 +
命題3<br/>
 +
ⅰ)一様な重力のもとで静止している気体・液体内では、同一水平面上での圧力の大きさは一定である。 <br/>
 +
ⅱ)もし液体の密度ρが圧力によって変化しないならば、 <br/>
 +
深さl1の水平面H1上の圧力p1と  <br/>
 +
深さl2(l2>l1)の水平面H2上の圧力p2には  <br/>
 +
次の関係が成り立つ。 <br/>
 +
p2=p1+ρg(l2l1)  <br/><br/>
 +
 
 +
図示した液体部分Vが静止しているので、Vに作用する力の総和が零になっている。<br/>
 +
このことから、この命題は容易に証明できる。<br/>
 +
 
 +
命題4 アルキメデスの原理<br/>
 +
*[[wikipedia_ja:アルキメデスの原理|ウィキペディア(アルキメデスの原理)]] 
 +
 
 +
 
 +
====気体の圧力と大気圧====
 +
気体は圧力が増すと縮むので、命題3のⅱ)の結論は成立しない。<br/>
 +
大気は静止していると仮定し、地表の大気圧から高度zでの大気圧を求めてみよう。
 +
地表の一点を原点とし、鉛直上方をz軸の正方向になる座標Oxyzをいれる。<br/>
 +
図のように、下底面が高さz、上底面が高さz+hの、単位断面積の角柱Vを考える。<br/>
 +
その部分の気体が受ける力の和は零となるので、<br/>
 +
次式が成り立つ。<br/>
 +
p(z+h)+mg=p(z)(1)<br/>
 +
ここで<br/>
 +
p(z)は高さzの地点の大気圧(命題3のⅰ)から、高度が同じ水平面上で圧力は一定)、<br/>
 +
mVの質量。Vの体積hと平均質量密度ρの積。<br/>
 +
圧力が大きいと空気は縮み質量密度は高くなるので、両者の関係を求めねばならない。<br/>
 +
空気体積の変動にともなう温度変化がないとすると、<br/>
 +
ボイルの法則(3章1節 熱とエネルギー参照)から、<br/>
 +
pVm=ccは温度だけに依存する数)<br/>
 +
質量密度ρ=mVを代入すると、<br/>
 +
pρ=c,ゆえに、ρ=pc<br/>
 +
1cを、cとおくと、<br/>
 +
ρ=cp  (2)<br/>
 +
この質量密度と圧力の関係を用いると、<br/>
 +
m=hρhcp(z)(hが小さいほど差は少なくなる)<br/>
 +
この式を(1)式に代入して、<br/>
 +
p(z+h)+cgp(z)hp(z)、変形すると <br/>
 +
p(z+h)p(z)hcgp(z)。これより<br/>
 +
dp(z)dz=limh0p(z+h)p(z)h=cgp(z)<br/>
 +
を得る。これを積分して<br/>
 +
p(z)=p0ecgz<br/>
 +
を得る。<br/>
 +
ここでp0は、地表での圧力、eは[[wikipedia_ja:ネイピア数 |ネイピア数]]である。<br/>
 +
地表での質量密度がρ0ならば,(2)式から、<br/>
 +
c=ρ0p0
 +
 
 +
=== 圧力の単位      ===
 +
圧力は、単位面積当たりの力なので、その単位は面積の単位m2と力の単位Nから得られる。<br/>
 +
Pa=N/m2=kgm1s2<br/>
 +
が圧力の単位で、パスカルと呼ばれる。

2015年2月5日 (木) 05:58 時点における最新版

目次

[非表示]

解説

この節では複数の質点が集まって作る質点系と、硬くて形を変えない質点系である剛体の運動を、運動法則を用いて解析しよう。

質点系の運動

2個以上の質点が集まって出来ている系を質点系という。
質点系というときは、各質点は密集していても、離れ離れでも良い。互いに固着しようが、自由に動けようが構わない。
すべての物質は、分子の集合と考えたり、細分化して極小部分に分け、それらの集合と考えれば、十分な精度で、質点系とみなすことができる。
そのため質点系の運動の法則を、ニュートンの運動法則から導出すれば、その応用範囲は非常に広い。

質点系の運動と重心

系の任意の2つの質点間には作用・反作用の法則を満たす力が働いていてもよい。
この力を質点系の”内力”という。  
質点系の各質点に外部から力(外力という)が加わる時、この質点系はどんな運動をするだろうか。
質点系の各質点の位置をri、質量をmiとし、
質点mi に作用する外力をfi
mi に、他の質点mjから作用する内力をfijとする(i,j=1N)。
すると、各質点に対して、運動の第2法則により、
d(mivi)/dt=fi+jifij  ここでvi=dri/dt
各ベクトルを自由ベクトルとみなしてi=1Nについて加え合わせると、fij+fji=0なので、
d2dt2imiri=ddtimivi=ifi
が得られる。
質点系の全質量M=imiと質点系に働く全外力F=ifiを用いて書きなおすと、
Md2dt2(imiri/M)=F
質点系の重心R R=imiri/M で定義すると、
Md2dt2R=F
この式は、力Fをうける質量Mの質点の運動方程式と同じである。
以下の解説も参考にしてください。

複雑にみえる運動も重心の運動をみれば簡単である  

体操選手の運動は、跳躍などで空中をまいながら、回転や体の屈伸、ひねりなどを行う。大変複雑で美しい。
しかし、導出した質点系の重心の運動法則から、体の重心の運動は、投射体の運動であり、放物線をえがいて移動することが分かる。
空中に飛び出た瞬間の重心の位置と速度(速さと方向・向き)で、その軌跡は完全に決まってしまうのである。

剛体の運動とつり合い

剛体

剛体(Rigid body)とは、
質点系であって、それらの、どの2質点の間の距離も変わらない,特殊な質点系のことを言う。
どの2質点の間の距離も変わらなければ変形は起こらない。
固くて変形しにくい物体を理想化した概念である。

剛体の運動 

剛体は変形しない質点系なので、その運動は、重心の運動と、重心の周りの回転運動を合成したものになる。
重心の運動は前の節で説明したように、質点の運動と同じように簡単に扱える。
重心の周りの回転運動について解析するには、少し難しい数学が必要になる。

このテキストでは、固定軸の周りの回転運動を中心に、 剛体運動の初歩と釣合の条件について学ぶ。

固定軸のまわりの回転運動 

剛体が、剛体の中を通る固定軸の周りを回転する運動(車輪の回転など)を考える。
応用も考え、回転軸は重心を通らなくてもよいように一般化しておく。
(注)なお、軸が動かないようにするためには軸受が必要である。
工夫しても回転時に軸は軸受から多少の摩擦力を受け、回転にブレーキがかかる。
しかし、これは無視出来るほど小さいと仮定する。
すると軸が受ける力は、軸の変動を防ぎ、固定軸の周りの運動に限定させる作用を持ち、
回転を遅める作用は持たないことになる。

回転運動の表示法 

固定軸まわりの剛体の運動はどのように表示したらよいだろうか。
・剛体の位置を表す変数;回転角
剛体が幾ら回転したか分かるように、剛体の、回転軸上にない一点Psに印を付ける。
次に、角度を測る基準線をきめるため、座標系を決めよう。
Psから固定軸へ垂線をひき、その足を原点Oとし,固定軸をz座標とする(静止した)3次元直交座標Oxyzを考える。
剛体が固定軸の周りを回転すると、印Psはxy平面上を、原点Oを中心に円を描いて動くことになる。
その位置ベクトルOPsがx軸の正方向となす角度ϕを、回転角と呼ぶ。図参照。
但し、x軸から反時計回りの角を正にする。
また一回転した後ならば、一回転の角2πを加え、逆周りに一回転した後なら2πを引き、
角度だけでなく回転数も分かるようにする。
回転角が指定されると、点Psの位置が決まる。
それだけでなく剛体は変形しないので、剛体のすべての点の位置がきまる。
そこで回転角ϕの時間変化ϕ=ϕ(t)を明らかにすれば、剛体の回転運動は定まる。
固定軸のまわりの回転運動において回転角の果たし役割は、質点の運動において質点の位置が果たし役割に対応していることが分かる。
・回転の角速度と角加速度
ϕ=ϕ(t)を時間で微分したdϕ(t)/dtを回転の角速度と呼ぶ。
直観的には、時刻tの瞬間の、回転の速さ(回転角の時間に対する変化率)を表す。
さらにもう一回時間微分したd2ϕ(t)/dt2を回転の角加速度と呼ぶ。

回転力(トルク) 

質点の運動に倣って、剛体に作用する力によって、その位置(=回転角)がどう変化するかの法則を導出したい。
しかし、剛体の回転の場合、ある方向の力は、剛体の回転に全く関係しない。 例えば、回転軸から放射状にでる半直線方向の力は全く回転の変化に寄与しない。
そこで剛体の回転を変化させる力とはなにかという問題から考察する必要が起こる。
質点運動における力の定義(力と運動量の変化の関係)や力と仕事の関係など力の係っている式のなかから、
剛体の回転運動に容易に拡張出来るものを選び、その式から、回転に関する力を求めることを試みる。
力の定義からは、回転運動への拡張を、推測することは難しい。
力と仕事の関係の考察をしてみよう。

力と仕事の関係からの考察 

適当な直交座標系をさだめ、ベクトルは、座標成分で表示する。
質点に、一定の力F=(Fx,Fy,Fz)を作用させて、x軸方向に変位させる。
質点はこの軸の上でしか動けないように拘束され、摩擦はないと仮定する。
質点の変位ベクトルは一次元の変数xを使ってs=(x,0,0)と表せる。
すると力のなす仕事は、W=F(x,0,0)=Fxxである。 
逆に物体に一定の力を加え、x軸上でxだけ変位させた時の仕事Wが分かれば、質点を動かした力は
Fx=W/x
で求められる。
Fy,Fzは、質点をx軸上で動かすことには全く寄与せず、
x軸に拘束された質点を動かす力は、Fxなのである。
固定軸まわりの回転もその変位は一次元の変数である回転角度で表わせるので、
これに倣って、
W/回転した角度 
を、回転にかんする力であると考える。これを回転力と呼ぶ。トルクともいう。

この方針を実行して回転力を具体的に求めよう。

剛体に力を加え微小角動かす時の、力のなす仕事の算出 

図4.1のように剛体の任意の一点P(x,y,z)を考える。
z座標の上方からxy平面を見下ろしているので、z座標は点になりOと書いてある。

図4.1 ☆☆キャプションはココに書いて下さい☆☆

まず一点P(x,y,z)に力F=(Fx,Fy,Fz)が作用して、微小角Δθだけ回転したときの
仕事ΔWを計算し回転力を求めよう。
P点から回転軸(z軸)に垂線を下ろし、その足をO=(0,0,z)とする。
OPの長さをr、x軸となす角をθ(ラジアン)と置く。
この角度は、
剛体につけた印の位置ベクトルOPsがx軸となす回転角ϕ
このベクトルとOP(をxy平面に平行移動したベクトル)の間の角の和である。
後者は、剛体なので、運動しても変わらない定数である。そこで、θ=ϕ+定数,と書ける。
剛体がz軸の周りを微小角Δθ回転して、点Pが図の点Qに移動したとする。
すると角OPQはほぼ直角(=π/2)でPQの長さPQは、PQ=r(Δθ)

PQのx成分とy成分は、図4-1中に示したように、それぞれ、QR=PQy/rPR=PQx/r
PQ=r(Δθ)を代入すると、
PQx=y(Δθ)PQy=x(Δθ)PQz=0
P(x,y,z)に作用する力F=(Fx,Fy,Fz)が、物体をPQだけ動かしたので、
その仕事は、ΔW=FPQ(内積)。
この右辺を内積の性質を用いて座標成分で表すと、
Fx(y)Δθ+FyxΔθ+Fz0
=(xFyyFx)Δθ

z軸まわりの回転力の導出 

ゆえに、力Fのz軸まわりの回転力(トルク)TezΔW/Δθ=xFyyFx に等しい。
これより、ΔW=TezΔθが得られる。
この式と、直線上に拘束された質点の運動における、力と仕事の関係式(  節  項)と対比させると、
Tez は、拘束された直線の上を動かすときに、働いた力の成分が対応し、
Δθ は、変位量   に対応していることが分かる。

z軸まわりの回転力(トルク)の性質

(1)力Fのz軸まわりの回転力は,Fzには関係しない。
言いかえるとz軸を固定軸とする剛体にz軸の方向の力を加えても、z軸の周りの回転は起こらない。
(2)剛体の1点P(x,y,z)に作用する力Fを考える。
P(x,y,z)からz軸に下ろした垂線の足をO(0,0,z)と書く。 力Fを、, OP方向の成分Frと、
z軸まわりの回転によりPの描く、Oを中心とする回転円の(左回りの)接線方向の成分Ft
および、これら2成分に直交する成分(z軸と平行)
に分解する(図参照)。この時、
・力Frのz軸まわりの回転力は、零である。
すなわち、動径方向の力は回転に寄与しない。  
・力Fのz軸まわりの回転力は、Ftのz軸まわりの回転力に等しい。
数式で表すと、xFyyFx=x(Ft)yy(Ft)x
(3)剛体に作用する力の作用点を、力の作用線上で動かす限り、回転力は変化しない。
ここで、力の作用線とは、力の作用点を通り、力の方向と重なる直線のこと。

これらはいずれも直観と合致する。
証明は、試みてほしい。

他の軸の周りの回転力

Fのx軸、y軸まわりの回転力も同様に計算できる。結果は、
x軸まわりの回転力;yFzzFy=y(Ft)zz(Ft)y
y軸まわりの回転力;zFxxFz=z(Ft)xx(Ft)z


原点まわりの力のモーメント

位置ベクトルr=(x,y,z)の剛体の点Pに作用する力Fの原点まわりの力のモーメントを、
N=(x軸まわりのトルク、y軸まわりのトルク、z軸まわりのトルク)で定義する。
数式で書くと、
N=(yFzzFy,zFxxFz,xFyyFx),

ベクトル積と力のモーメントのベクトル積表示

以上の結果は、ベクトル積(クロス積ともいう)を用いると簡潔、正確に表現でき、
回転運動の性質を調べるのが容易になる。
3次元ベクトルa,b のベクトル積a×bとは、3次元ベクトルであり,
大きさはa,b を2辺とする平行四辺形の面積に等しく、
方向はこの四辺形に垂直で、向きは、(a,b,a×b)が右手系をなすように定めたものである。
すると、ベクトル積に関して以下の8つの命題が成り立つ。

ベクトル積にかんする命題    

以下に述べる全ての命題で、
a,b,cは3次元ベクトル
αは実数とする。
命題1.
a を, cと垂直な成分a と,平行な成分a の和に分解するとき、
a×c=a×c
a×c=0

命題2.
a×b=b×a

命題3
(αa)×b=α(a×b)=a×(αb) 

命題4.
(a+b)×c=a×c+b×c 

命題4の系  
   a×(b+c)=a×b+a×c
(a+b+c)×d=a×d+b×d+c×d

命題5. 
(e1,e2,e3)
それぞれ大きさ(長さ)1で互いに直交し、右手系をなす、ベクトル(右手系をなす正規直交基底)とする。
この時、
e1×e2=e3,e2×e3=e1,e3×e1=e2
 

命題6.
ベクトルa,bを,命題5で用いた基底(e1,e2,e3) で決まる座標の座標成分で表示しておく。
するとa×b=(aybzazby,azbxaxbz,axbyaybx) 

命題7.
(a×b)c=(c×a)b=(b×c)a

命題8.
a(t)b(t)を,tにかんして微分可能な、ベクトルに値をとる関数とする。すると、
a(t)×b(t) は、tにかんして微分可能で、
ddt(a(t)×b(t))=(ddta(t))×b(t)+a(t)×(ddtb(t))

これらの証明は、後節で扱う。 ベクトル積に関しては以下を参照のこと。

 力のモーメントのベクトル積表示   

ベクトル積の命題6を用いると、
位置ベクトルrの点に作用するF
原点まわりの力のモーメントは、N=r×F
x軸まわりの回転力(トルク)は、Nex と表せることが分かる。
y軸とz軸周りの回転力も、それぞれ 
Ney ,Nezで 表せる。

力のモーメントの性質

もっと一般に、どんな軸の周りの回転力も、N から得られる。
定理;
eを、原点を始点とする大きさ1の任意のベクトルとする。
すると、
Neは、力Fe軸の周りの回転力になる。式で書くと、Te=Ne 
この式を、回転力の定義に基づいて言い換えると、
Fのもとで、剛体をe軸の右まわりに角度ϕだけ回転させたとき、 Fのなす仕事Wは、W=Teϕ=(Ne)ϕ
証明;
9つに分けて示す。
ⅰ)準備 
図のように、剛体の点 P から、e 軸に垂線を下ろし、その足を Q とする。
力 F のもとで、剛体が e を固定軸にして、
微小時間に、微小角δϕ だけ回転したとする。
このとき、P が移った先を、P とする。
ⅱ)回転角 δϕ が微小なので、
この回転中の P の軌跡(円弧の微小部分)は、有向線分PP で精度高く、近似できる。
ⅲ)この間に力 F がなした仕事 δW は、δW=PPF
この仕事を、回転角δϕで割ると、力の e 軸周りの回転力が得られる。そこで、PP を、この定理で与えられている諸量を使って表現し、これを用いて、仕事を計算しよう。
ⅳ)有向線分PPの方向を求める。
PP は、e 軸と垂直でQ を通る平面H上にあり、
Qを中心とする円の弧の微小部分をなすので、線分QP と直交する。PPQP
また、e 軸と垂直でQ を通る平面H上にあるので、 PPe 軸とも直交し、従って線分OQと直交する。PPOQ
ゆえに、PP は、3点O,Q,Pを通る平面 OQP と直交する。
すると、PP は、平面 OQP 上のすべての線分と直交する。
ゆえに、PPe,PPOP 
これで、PP の方向は、求まった。
ⅴ)有向線分PP の向き 
点 P は、e 軸の周りを右周りに回転するので、その向きは、 e×OP と同じ向きである。
ⅵ)PP の大きさ。
PPは、 Q を中心とする、半径 QP の円弧の一部なので、 その中心角δϕ を用いて、PP=QPδϕ
ⅶ)ⅳ)、ⅴ)、ⅵ)から  PP=e×re×rQPδϕ
ⅷ)PP=e×rδϕが成り立つ。
なぜなら、
e×r=ersinθ=rsinθ=QP ,ここで θ はer の間の角。
この式をⅶ)で得られた式に代入すれば、所望の結果が得られる。
ⅸ)δW=PPF=(e×rδϕ)F=(e×r)Fδϕ=(r×F)eδϕ
ⅹ)Te=δWδϕ=(r×F)e=Ne
定理の証明終わり。

(注)剛体が固定軸の周りでなく、自由に回転するときでも、
ある瞬間には、ある軸の周りの回転になっている。
力のモーメントは、どんな軸周りの回転力の情報も含んでいることが証明されたので、
  回転運動一般に有効な概念であることが分かる。

剛体の複数個所に作用する力の回転力 

次に剛体の多くの点に力を加えたときの回転力を求めよう。
力の作用点をPi(xi,yi,zi)、力をFi(i=1,2,,,n)とする。
これらの力のもとで剛体がz軸まわりをΔθだけ微小回転するときの、各力のなす仕事の合計は、
(ni=1(xi(Fi)yyi(Fi)x)Δθ
従って、作用点Pi(xi,yi,zi)の力Fi(i=1,2,,,n)の全体がもつz軸まわりの回転力は、
Tez=ni=1Tiez=ni=1(xi(Fi)yyi(Fi)x) ここでTiezは力Fiのz軸まわりの回転力。

同様に、x軸まわりとy軸まわりの回転力も、それぞれ
Tex=ni=1Tiex=ni=1(yi(Fi)zzi(Fi)y)
Tey=ni=1Tiey=ni=1(zi(Fi)xxi(Fi)z)
Fiの原点周りに力のモーメントNiNi=(Tiex,Tiey,Tiez)で定義した。
全ての力の原点周りの力のモーメントも、同様に
N=(Tex,Tey,Tez)で定義する。すると、
N=(ni=1Tiex,ni=1Tiey,ni=1Tiez)=ni=1Ni
全ての力の原点周りの力のモーメントNも、上述の定理と同様の定理(定理の系と呼ぶ)が成り立つ。
定理の系
Nを剛体に作用する全ての力のモーメントとし、
eを、原点を始点とする大きさ1の任意のベクトルとする。
すると、
Neは、力Fe軸の周りの回転力になる。
式で書くと、Te=Ne 
この式を、回転力の定義に基づいて言い換えると、
Fi(i=1,2,,,n)のもとで、剛体をe軸の右まわりに角度ϕだけ回転させたとき、
これらの力のなす仕事Wは、W=Teϕ=(Ne)ϕ

この系は、内積の性質を使えば、定理から、容易に導かれる。

質点系に作用する重力のモーメント

n個の質点系を考える。
第i質点の質量をmi、位置ベクトルをriとする。
鉛直上方をz軸の正方向とする直交座標系0xyzをいれる。
この質点系に作用する重力の原点周りのモーメントNを求めよう。
第i質点に働く重力は、
fi=(0,0,mig)
なので、
N=ni=1ri×fi=ni=1ri×(0,0,mig)
=ni=1(miri×(0,0,g))=(ni=1miri)×(0,0,g))
すでに学んだことから、この質点系の重心は、
R=ni=1miri)M
であった。ここで、 M=ni=1mi 。
これを用いて、モーメントを書きなおすと、
N=MR×(0,0,g)=R×(0,0,Mg)
となる。
これは、質点系の重心の位置に質点系の全質量が集中している時の、 原点周りの重力のモーメントに等しい。

回転運動の方程式 

N が、あらゆる回転軸にかんする回転力を表現していることがわかった。
Fと運動量の変化の関係をあたえるニュートンの運動方程式(第2法則)を変形して、
回転力Nにかんする方程式を導こう。
直交右手座標系Oxyz を定める。原点 O は、考察対象に都合のよい点を選ぶ。

剛体をN個の(質点と考えてよい)微小部分Pi(i=1N)に分け、
その質量をmi、位置ベクトルをri(xi,yi,zi)とする。
Piが外部から受ける力をFi
Pi が剛体の他の部分Pj(ji) から受ける力(内力)をFijとおく。
後者は、剛体が変形しないよう、剛体の原子間に働かせる力に起因する。
この原子間の力は、原子の電荷による電気力と、
原子同士が接近しすぎたときに作用する量子力学的力により生じる。
作用・反作用の法則(運動の第3法則)から、Fij=Fji 。
さらに、剛体の2点間に働く内力の方向は、
その2点を結ぶ直線の方向と同じだと、仮定する。

各質点のニュートンの運動方程式  

各質点ごとに、ニュートンの運動方程式を立てると、
mid2ridt2=Fi+jiFi,j(i=1N)
これを変形して
Fi=mid2ridt2jiFi,j(i=1N) (1) 
この式から、
Fiの回転力Ni=ri×Fiにかんする式を導こう。

Ni=ri×Fiにかんする式の誘導  

式(1)の両辺に左側から、ri のベクトル積を施すと、
Ni=ri×Fi=ri×(mid2ridt2jiFi,j) (i=1N)
ベクトル積の性質3と性質4により、
=miri×d2ridt2jiri×Fi,j (2)
ここで、ベクトル積の性質8より
ddt(ri×dridt)=dridt×dridt+ri×d2ridt2=ri×d2ridt2
なので、 Ni=middt(ri×dridt)iri×Fi,j=ddt(ri×midridt)jiri×Fi,j(3)
質点Piの運動量をPiと書くと、
Pi=mivi=midridtなので、
Ni=ddt(ri×Pi)jiri×Fi,j
定義;角運動量(運動量のモーメントともいう)
質点の位置ベクトルをr、運動量をpと書くとき、
l=r×pを,この質点の角運動量と呼ぶ。
これを用いると、
Ni=dlidtjiri×Fi,j

 回転の運動方程式の導出  

故に、
N=iNi=dilidtijiri×Fi,j(4)
ここで、
ijiri×Fi,j=i<jri×Fi,j+i>jri×Fi,j(5)
式(4)の右辺の第2項の上付き添え字i,jを、それぞれ、j'と i'でおきかえられるので、
i>jri×Fi,j=j>irj×Fj,i
内力は作用反作用の法則が適用できると仮定しているので、
Fj,i=Fi,j 。この式を上の式の右辺に代入すると、
i>jri×Fi,j=j>irj×Fi,j
この式の右辺の和をとる変数i',j' を i,j におきかえると、
i>jri×Fi,j=i<jrj×Fi,j
この式を、式(5)の右辺の第2項に代入して整頓すると、
ijiri×Fi,j=i<j(rirj)×Fi,j
さらに、内力に関する第2の仮定により、rirj とFi,jは同じ方向なので、ベクトル積の定義より、この項は、零となることが分かる。
故に、式(4)の右辺の第2項は零となり、
N=dilidt(6)
が得られる。全角運動量をL=iliとおけば、
式(6)は、次のように書ける。

命題;回転運動の関するオイラーの運動方程式
剛体の内力に上述の2つの仮定を付ける。このとき、
剛体に作用する全ての外部力の原点周りの力のモーメントN=iNi=iri×Fiと、
全角運動量L=ili=iri×piの間には、
N=dLdt (7)

この命題の導出までは詳しく述べたが、本テキストではこれ以上は深入りしない。
この先にも興味がある方は、次の記事をご覧ください。

固定軸の周りの剛体の回転運動の方程式

回転運動の運動方程式から、任意の軸の周りの回転運動の方程式が簡単に導出できる。
z軸周りの場合を例にとり、説明する。
z軸周りの回転力はTez=Nezなので、
回転運動の方程式から
Tez=Nez=dLdtez
この式の右辺に,L=iri×pi を代入すると
右辺
=diri×pidtez 微分の加法性から
=(idri×pidt)ez 内積の加法性から
=i(dri×pidtez) ベクトル積の性質8から
=i(dridt×pi+ri×dpidt)ez pi=midridtを代入し、ベクトル積の性質を用いると、
=i(ri×d2ridt2)ez
故に、
Tez=i(miri×d2ridt2)ez(1)
剛体はz軸の周りを回転するので、
その各点Pi(位置ベクトルri=OPi)は、
z軸と直交する平面上を、z軸を中心とする円を描いて運動する。
この拘束条件を考慮して、
時刻tの位置ベクトルri(t)の座標成分を書きなおすと、
ri(t)=(xi,yi,zi)=(ˆricosθi(t),ˆrisinθi(t),zi)(2)
ここでˆriは、点Piとz軸との距離、
θ(t)は、ri(t)をxy平面に正射影した像がx軸となす角度である。図参照。
剛体につけておいた印Psの位置ベクトルOPs
xy平面に正射影した像がx軸となす角(回転角)ϕを用いると、
θi(t)=ϕ(t)+ϕi (3)
ϕiは、Piごとに決まる、定数)と書ける。

式(1)の右辺を、式(2)を利用して、変形すると、
=imi((ˆricosθi(t),ˆrisinθi(t),zi)׈ri(cosθi˙θi2sinθi¨θi,sinθi˙θi2+cosθi¨θi,0))ez
=imiˆri((ˆricosθi(t),ˆrisinθi(t),zi)×(cosθi˙θi2sinθi¨θi,sinθi˙θi2+cosθi¨θi,0))3
ベクトル積の性質6より、
=imiˆri
(ˆricosθi(t)(sinθi(t)˙θi(t)2+cosθi(t)¨θi(t))ˆrisinθi(t)(cosθi(t)˙θi(t)2sinθi(t)¨θi(t))
=\sum_i m_i(\hat{r}_i)^2\ddot{\theta_i}(t)
ここで、\theta_i(t)=\phi(t)+\phi_iを代入すると
=(\sum_i m_i(\hat{r}_i)^2)\ddot{\phi}(t)
以上により、
T_{\vec e_z}=\vec N \cdot \vec e_z 
=(\sum_i m_i(\hat{r}_i)^2)\ddot{\phi}(t)
が得られた。 I=\sum_i m_i(\hat{r}_i)^2とおくと、この式は
T_{\vec e_z}=I \ddot{\phi}(t) \qquad (4)
と書ける。ここでIを、剛体の軸まわりの慣性モーメントと呼ぶ。
これがz軸を固定軸とする剛体の回転運動の運動方程式である。
原点を始点とする任意の回転軸\vec{e},\|\vec{e}\|=1まわりの回転の方程式も同様に得られる。

この方程式の変数\phi は、一次元のスカラーなので、
質点がなめらかに拘束され、直線上を運動するときの運動方程式
F=m\ddot{x}
と、対比させる。すると、
質点に作用する力 F  <===> 剛体に作用する回転力T_{\vec e}=\vec N \cdot \vec e
質点の質量 m     <===> 剛体の\vec{e}軸まわりの慣性モーメント
\qquad \qquad \qquad \qquad \qquad \qquad I=\sum_i m_i(\hat{r}_i)^2\hat{r}_iは質量m_i\vec{e}軸を延長した直線との距離
質点の位置変数 x(t)  <===> 剛体の\vec{e}軸周りの回転角変数\phi(t)
質点の速度 \dot{x}=\frac{dx(t)}{dt} <===>剛体の\vec{e}軸周りの角速度\dot{\phi}(t);
質点の運動量 m\dot{x} <===> 剛体の角運動量I\dot{\phi};
運動方程式F=m\ddot{x} <===> T_{\vec e}=I\ddot{\phi}

という、対応関係があることが分かる。
この節で得た固定軸まわりの回転運動の方程式から、
もし\vec N=0 ならば、任意の軸まわりの回転力が零なので、
剛体の任意の軸まわりの角加速度が零、角速度が一定となることが分かる。

剛体の回転の運動エネルギー  

剛体の各微小部分(質量m_i)の速度を v_iと書くと、
その運動エネルギーは \frac{1}{2}m_i {v_i}^2,(i=1 \cdots n)なので、
剛体全体の運動エネルギーは、K=\sum_{i}\frac{1}{2}m_i {v_i}^2
回転運動している各微小部分の速度は、v_i=\hat{r}_i\dot{\phi}と書けるので、
K=\sum_{i}\frac{1}{2}m_i {\hat{r}_i}^2 {\dot{\phi} }^2=\frac{1}{2}I{\dot{\phi} }^2,\qquad (5)

物理振り子

剛体は、重心を通らない水平軸の周りで、重力の作用を受け振動する。
これを物理振り子、あるいは実体振り子という。

水平回転軸をx軸とし、鉛直上方をz軸の正方向とし、yz平面が剛体の重心を通る座標系を考え、
回転軸とこの平面の交点を原点O、重心をGと記す。図参照。
回転はなめらかで摩擦力は無視できるとする。
すると、回転軸から、この剛体が受ける力は、剛体をこの軸に支える作用を持つだけで、剛体の振動に何の影響も与えない。
そこで、剛体にかかる力は、重力だけと考えて良い。
重力の原点周りの力のモーメント\vec Nは、
剛体の重心\vec Rに、剛体の全質量Mがあるとしたときの
重力の原点周りのモーメントに等しいことが分かっている。 故に、
\vec N=\vec R \times (0,0,-Mg)=(-R_{2}Mg,R_{1}Mg,0)
x軸まわりの力のモーメントは、
\vec N \cdot \vec e_{x}=-R_{2}Mg=-Mg\|\vec{OG}\|\sin\phi
従って、回転の運動方程式は
I\frac{d^{2}\phi}{dt^2}=-Mg\|\vec{OG}\|\sin\phi
ここでIは、軸まわりの、振り子の慣性質量。

剛体の慣性モーメントの計算(一次元の剛体) 

剛体Vは、ごく細く、まっすぐな棒で,
長さl、質量密度(単位長さあたりの質量)は一定で\rhoとする。
棒の左端からl_1の場所Oを通り、棒に直交する軸まわりの慣性モーメントを具体的に計算しよう。
Oを原点とし、棒と同じ方向の数直線を考え、これを座標系として採用。
V=[a=-{l_1},b=l-l_1]と表現する。
剛体Vの慣性モーメントは、
剛体を質点とみなせるほど細かい部分V_i=[x_{i-1},x_i],(i=1,2,,,n)に分割して、
V_iの質量m_iと、V_iOとの距離\hat{r}_iを用いて、
I=\sum_i m_i(\hat{r}_i)^2で定義した。但しx_0=a,x_n=b

剛体の分割と慣性モーメントの近似式・リーマン和 

V_iの質量m_iは、V_iの長さx_i-x_{i-1}に質量密度\rhoを掛ければ得られるので
m_i=\rho (x_i-x_{i-1})であり、
I=\sum_i \rho (\hat{r}_i)^2(x_i-x_{i-1})
と書ける。
しかし、剛体V=[a,b]をいくら細かく分割しても、
各小区間V_i=[x_{i-1},x_i]は大きさ(長さ)をもつので、
原点との距離\hat{r}_iは、一つに定まらない。
そこで、各小区間V_iから、代表点\xi_iを選びだし、その点の原点からの距離|\xi_i|、(\xi_i絶対値)を、\hat{r}_iとみなす。
すると、慣性モーメントIの式は
\sum_i m_i(\hat{r}_i)^2=\sum_i \rho(\xi_i)^2(x_i-x_{i-1})
で近似される。

そこで、この分割を
\Delta=\{V_i=[x_{i-1},x_i] \mid i=1,2,,,,n\} と表し、 I^{\Delta}(\xi_1,,,\xi_n)
で,慣性モーメントの近似式を表すことにする。
すると、

慣性モーメントの近似式は、
I^{\Delta}(\xi_1,,,\xi_n)=\sum_{i=1}^{n}\rho(\xi_i)^2(x_i-x_{i-1})\qquad (1)
と書ける。
この値は分割の仕方と分割小区間の代表点\xi_i(\in V_i)の選び方によって変化する。
V_iの中で、原点に最も近い点{{\xi}^m}_i(\in V_i)(i=1,2\cdots,n)にとると
最小値 I_{m}(\Delta):=\sum_{i=1}^{n}\rho({{\xi}^m}_i)^2(x_i-x_{i-1})
をとり、
V_iの中で、原点に最も遠い点{{\xi}^M}_i(\in V_i)(i=1,2\cdots,n)にとると
最大値 I_{M}(\Delta):=\sum_{i=1}^{n}\rho({{\xi}^M}_i)^2(x_i-x_{i-1})
を取る。
関数y=f(x)=\rho x^2を使って表現すれば、
I^{\Delta}(\xi_1,,,\xi_n)=\sum_{i=1}^{n}f(\xi_i)(x_i-x_{i-1})
であり、I_{m}(\Delta)\leq I^{\Delta}(\xi_1,,,\xi_n)\leq I_{M}(\Delta)を満たす。
質量密度が場所で変わるときは、、関数はy=f(x)=\rho(x) x^2になり、
剛体の重心を求めるときは、後述するように、別の関数が現れる。
そこで、数学の分野では、一般の関数y=f(x)にたいして
I^{f,\Delta}(\xi_1,,,\xi_n)=\sum_{i=1}^{n}f(\xi_i)(x_i-x_{i-1})\qquad (2)
を求め、分割\Delta=\{V_i=[x_{i-1},x_i] \mid i=1,2,,,,n\} V_iの代表点\xi_i,(i=1,2,,,n)に関する関数y=f(x)リーマン和と呼ぶ。
その最小値I_{m}(f,\Delta)と最大値I_{M}(\Delta)も,同様に定義される。

I_{m}(f,\Delta)\leq I^{f,\Delta}(\xi_1,,,\xi_n)\leq I_{M}(\Delta)\qquad (3) 慣性モーメントの近似式(1)は、関数y=f(x)=\rho(x) x^2にたいするリーマン和である。

慣性モーメントの近似式の意味 

今後、関数y=f(x)は、V=[a,b]で定義された有界関数として、 議論を進める。
有界関数とは、十分大きな正数Mを選べば、
V=[a,b]の全ての点xに対して、|f(x)| \leq Mとなること。
y=f(x)=\rho x^2を代入すれば、考察対象の剛体の慣性モーメントの話になる。
リーマン和
I^{f,\Delta}(\xi_1,,,\xi_n)=\sum_{i=1}^{n}f(\xi_i)(x_i-x_{i-1})
は、y=f(x)のグラフを、棒グラフで近似したときの棒グラフの作る面積(各角柱の面積和)であることが分かる。図参照。
また、I_{m}(f,\Delta)は一点鎖線でしめす、小さいほうの長方形の和であり、
I_{Mm}(f,\Delta)は点線でしめす、大きいほうの長方形の和である。


I^{f,\Delta}(\xi_1,,,\xi_n)は、y=f(x)のグラフとx軸およびy軸と平行な直線x=ax=bで囲まれる部分の面積Sを近似している。
また、I_{m}(f,\Delta)\leq S \leq I_{Mm}(f,\Delta)  \qquad (4)  
であり、
I_{m}(f,\Delta)は面積を下から評価し、
I_{M}(f,\Delta)は面積を上から評価していることがわかる。
分割を限りなく細かくしていくとき、
リーマン和が分割や代表点の選び方に関係ない数に収束するならば、
その極限値は、
y=f(x)のグラフとx軸およびy軸と平行な直線x=ax=bで囲まれる部分の面積
と考えられる。
もし、分割\Deltaを細かくしていくとき
I_{m}(f,\Delta)I_{M}(f,\Delta)が同じ値に収束することが示せれば、
(3)式と(4)式から、リーマン和は、関数のグラフの作る面積Sに収束することが分かった。

可積分の定義と積分 

「分割を細かくしていくとき、リーマン和が収束する」ということは、
面積を決める上で決定的に重要がことなので、
可積分という名を付けて、数学的に厳密に定義する。 このためにはまず、分割の大きさを定める必要がある。
定義:分割の大きさ
分割 \Delta=\{V_i=[x_{i-1},x_i] \mid i=1,2,,,,n\} の大きさとは、
d(\Delta):=max_{i=1,2,\cdots n}(x_i-x_{i-1})

定義:可積分と積分
fを、有界閉区間V上で定義され、実数の値をとる関数とする。

もし、ある実数Iが存在して、
どんな分割\Delta=\{V_i=[x_{i-1},x_i] \mid i=1,2,,,,n\} と代表点\xi_i\in V_i(i=1,2,\cdots ,n)であっても、
\lim_{d(\Delta) \to 0}I^{f,\Delta}(\xi_1,,,\xi_n)=I
が成り立つ時、
fV上で(リーマン)可積分であるという。
このとき、IfV上での(リーマン)積分といい、
I=\int_{V}f=\int_{V} f(x)dx
などと書く。

 積分の性質 

定理(積分の線形性)
f, g \quadを、区間I上で定義された、任意の実数値関数であり、
c, d \quadを任意の実数とする。
このとき、
(1)f,\quad g \quadI上で可積分ならば、cf+dg \quadI上で可積分
(2)このとき、 \int_{I}(cf+dg)=c\int_{I}f+d\int_{I}g

証明;リーマン和の定義から、区間Iの任意の分割\Delta=\{I_1,,,,I_n\} と 分割区間の任意の代表点\xi\in V_i(i=1,2,,,,n) (\xiV_iに含まれる意)に対して、
S(cf+dg,\Delta,\{\xi_i\}_{i=1}^{n}) =cS(f,\Delta,\{\xi_i\}_{i=1}^{n}) +dS(g,\Delta,\{\xi_i\}_{i=1}^{n}) \qquad (1)
f,\quad g \quadは可積分なので、その定義から、
\lim_{d(\Delta) \to 0}S(f,\Delta,\{\xi_i\})=\int_{I}f
\lim_{d(\Delta) \to 0}S(g,\Delta,\{\xi_i\})=\int_{I}g
(1)式の両辺の極限\lim_{d(\Delta) \to 0} をとろう。
右辺の極限
=\lim_{d(\Delta) \to 0} \left(cS(f,\Delta,\{\xi_i\}_{i=1}^{n}) +dS(g,\Delta,\{\xi_i\}_{i=1}^{n}\right)
極限の性質から、
=c\lim_{d(\Delta) \to 0}S(f,\Delta,\{\xi_i\}_{i=1}^{n}) +d\lim_{d(\Delta) \to 0}S(g,\Delta,\{\xi_i\}_{i=1}^{n}
=c\int_{I}f+d\int_{I}g
従って(1)式の左辺の極限 \int_{I}(cf+dg) も存在して、右辺の極限と一致する。 証明終わり。

慣性モーメントの計算(1)リーマン和の極限を求める方法

Vは、先述の、ごく細い一様な質量密度\rho=M/lのまっすぐな棒で、
座標系を入れて、V=[a=-{l_1},b=l-l_1]と表現しておく。
原点を通りこの棒と直交する軸のまわりの(この棒の)慣性モーメントを、
リーマン和の極限を取って求めよう。
区間V=[-{l_1},l-l_1]をn(\geq 2)等分して得られる点列,
{x^n }_0=-l_1, {x^n }_1={-l_1}+l/n, {x^n }_i={-l_1}+i(l/n),,,{x^n }_n=l-{l_1}
を分点とする分割を{\Delta}^nと記す。すると、
{x^n }_i-{x^n }_{i-1}=l/n,\quad(i=1,2,,,n), d({\Delta}^n)=l/nであり、
{\Delta}^n=\{{V^n}_j=[{x^n}_{j-1},{x^n}_j] \mid j=1,2,,,n\}  
\{{\Delta}^n \mid n=2,3,,,\}という分割の列は、
\lim_{n\to\infty} d({\Delta}^n)=\lim_{n\to\infty}\frac{l}{n}=0を満たす。
y=f(x)=\rho x^2がリーマン可積分であることを認めれば、
可積分の定義から、どんな代表点{{\xi}^n}_j\in {V^n}_jを選んでも、
\lim_{n\to \infty}I^{f,{\Delta}^n}({{\xi}^n}_1,{{\xi}^n}_2,,,{{\xi}^n}_n)=Iとなる。

そこで、代表点を{{\xi}^n}_j={x^n}_j=-l_1+j(l/n) \quad (n=2,,,),(j=1,2,,,,n)と選ぶ。
関数y=f(x)=\rho x^2を用いると、 分割\Delta^nを用いた慣性モーメントの近似値は次のようになる。
I^{f,{\Delta}^n}({x^n}_1,{x^n}_2,,,,{x^n}_n) =\sum_j f({x^n}_j)1/n =\sum_j f(-{l_1}+j(1/n))\frac{l}{n} =\rho\sum_{j=1}^{n} (-{l_1}+j(1/n))^2\frac{l}{n}
ここで、 \sum_{j=1}^{n} j=\frac{1}{2}n(n+1),\quad \sum_{j=1}^{n} j^2=\frac{1}{6}n(n+1)(2n+1)(注参照)を利用して、この式を計算すると、
=\rho l ({l_1}^2-{l_1}l\frac{n+1}{n}+\frac{l^2}{6} \frac{n+1}{n} \frac{2n+1}{n})
\rho=M/lなので、
=M({l_1}^2-{l_1}l\frac{n+1}{n}+\frac{l^2}{6} \frac{n+1}{n} \frac{2n+1}{n})
故に、
I=\lim_{n\to \infty}I^{f,{\Delta}^n}({x^n}_1,{x^n}_2,,,,{x^n}_n) =\frac{M}{3}(l^2-3{l_1}l+3{l_1}^2)

(注)S_{1}:=\sum_{j=1}^{n} j=\frac{1}{2}n(n+1)の証明
(j+1)^{2}-j^{2}=2i+1 なので、両辺のj=1,2,,,n に関する和を取る。
左辺の和は\sum_{j=1}^{n}((j+1)^{2}-j^{2})=(n+1)^{2}-1
右辺の和は\sum_{j=1}^{n}(2j+1)=2\sum_{j=1}^{n}j+n=2S_{1}+n
故に、(n+1)^{2}-1=2S_{1}+n (n+1)^{2}-1-n=2S_{1} S_{1}=\frac{1}{2}\left((n+1)^{2}-1-n\right)=\frac{1}{2}n(n+1)
S_{2}:=\sum_{j=1}^{n} j^2=\frac{1}{6}n(n+1)(2n+1)の略証
(j+1)^{3}-j^{3}=3j^2+3j+1なので、この両辺のj=1,2,,,nに関する和を取る。
左辺の和は(n+1)^{3}-1、右辺の和は3S_{2}+3S_{1}+n,故に3S_{2}+3S_{1}+n=(n+1)^{3}-1


慣性モーメントの計算(2)原始関数を利用する方法

積分可能な関数の積分をリーマン和の極限から求める計算は煩雑であり、複雑な形状の剛体の慣性モーメントを求めるにはふさわしくない。
次の定理が強力な計算法を提供する。

定理
V=[a,b]を数直線上の区間、
fV上可積分な実数値関数
とする。
もしFが、
V上で微分可能で
全てのVの点xで、\frac{d}{dx}F(x)=f(x)
を満たす関数ならば(注参照)、
\int_{[a,b]}f=F(b)-F(a)
上記の条件を満たす関数Fを、f原始関数という。
(注)関数Fは、V上でしか定義されていないので、
端点a,bでは、通常の微分は定義できない。そこで、
\frac{d}{dx}F(a):=\lim_{h \to 0,h\geq 0}\frac{F(a+h)-F(a)}{h}
\frac{d}{dx}F(b):=\lim_{h \to 0,h\leq 0}\frac{F(b+h)-F(b)}{h}
と定義する。
証明;
区間[a,b]の任意の分割
\Delta=\{V_i=[x_{i-1},x_i]\mid 1 \leq i \leq n,x_0=a,x_n=b\}
に対して、
代表点を\xi_i\in V_i(\xi_iV_iの点の意)とすると、 fのリーマン和は
I^{f,\Delta}(\xi_1,,,\xi_n)
=\sum_i f(\xi_i)(x_i-x_{i-1})
小区間V_i=[x_{i-1},x_i]での関数Fの平均勾配
\frac{F(x_i)-F(x_{i-1})}{x_i-x_{i-1}}
は、平均値の定理により、
V_i=[x_{i-1},x_i]の中のある一点\eta_iにおけるy=F(x)の接線の勾配
\frac{d}{dt}F(\eta_i)に等しので、
\frac{F(x_i)-F(x_{i-1})}{x_i-x_{i-1}}=\frac{d}{dt}F(\eta_i)=f(\eta_i)
故に、f(\eta_i)(x_i-x_{i-1})=F(x_i)-F(x_{i-1})
そこで、各小区間V_iの代表点を\eta_i,(i=1,2,,,n)と選べば、
I^{f,\Delta}(\eta_1,,,\eta_n)
=\sum_i f(\eta_i)(x_i-x_{i-1})
=\sum_{i=1}^{n}\left(F(x_i)-F(x_{i-1})\right)
=F(x_n)-F(x_0)=F(b)-F(a)
fは可積分なので、 \int_{[a,b]}f=\lim_{d(\Delta)\to 0}I^{f,\Delta}(\eta_1,,,\eta_n)
=\lim_{d(\Delta)\to 0}(F(b)-F(a))=F(b)-F(a)
証明終わり。

さて、慣性モーメントを求めたい剛体では、
f(x)=\rho x^2なので、その原始関数は、
F(x)=\frac{1}{3}\rho x^3
従って、慣性モーメントは、定理を適用して、
I=\int_{[a,b]}\rho x^2=\frac{1}{3}\rho (b^3-a^3)
\rho=M/l,a=-l_1,b=l-l_1を代入して、整頓すると、
=\frac{M}{3}(l^2-3l_{1}l+3{l_1}^2)

重心の計算への応用 

質量密度が場所により変わる、長さlのごく細い棒Vの重心を求めてみよう。
考えやすくするため、 棒の一端を原点にし、他端がx軸の正の位置にくるように座標系O-xをいれる。
この座標系で剛体はV=[0,l]と書ける。 Vを小区間V_i=[{x^n}_{i-1},{x^n}_{i}],i=1,2.\cdots,n),{x^n}_n=0,{x^n}_n=lに分割(分割\Deltaと記す)し、これらの小区間を質点とみなせば、その重心は、
\vec{R}=\sum_i{ m_i \vec{r_i}}/M
で定義された(1.1.1節参照)。ここでm_i は第i質点の質量、M=\sum_{i} m_i\vec{r_i}は第i質点の位置ベクトル。
ベクトルを座標成分表示すると、この問題では一次元なので、 R=\sum_{i} m_i r_i/M M=\sum_{i} m_i
しかし、実際には V_iは、質点ではないので、 位置ベクトルは、定まらない。
またその質量も密度が一定ならば、m_i=\rho ({x^n}_{i}- {x^n}_{i-1})できまるが、
密度が変化するならば、定まらない。
そこで、各小区間 V_iの代表点\xi_i (\in V_i)を選び
m_i=\rho(\xi_i)({x^n}_{i}- {x^n}_{i-1}),r_i=\xi_i
で近似する。
すると分割\Deltaと代表点\{\xi_i\}_{i=1}^{n}に対応する、 質量Mと重心Gの近似値は、それぞれ
M^{\rho,\Delta}(\xi_1,,,\xi_n)
=\sum_i \rho(\xi_i)v(V_i)=\sum_i \rho(\xi_i)(x_i-x_{i-1})

G^{f,\Delta}(\xi_1,,,\xi_n)
=\sum_i f(\xi_i)v(V_i)=\sum_i f(\xi_i)(x_i-x_{i-1})
ここで、f(x)=\frac{1}{M}\rho(x)x
もし、関数\rho(x)が積分可能ならば、分割\Deltaを細かくしていけば M=\lim_{d(\Delta) \to 0}M^{\rho,\Delta}(\xi_1,,,\xi_n)
=\int_{[0,l]}\rho
もし関数\rhoの原始関数Pが存在する(\rho(x)=\frac{dP(x)}{dx})ならば
=P(l)-P(0)
もし関数f(x)も積分可能ならば、分割\Deltaを細かくしていけば G=\lim_{d(\Delta) \to 0}G^{f,\Delta}(\xi_1,,,\xi_n)
=\int_{[0,l]}f
もし関数fの原始関数Fが存在する(f(x)=\frac{dF(x)}{dx})ならば
=F(l)-F(0)
例;\rho(x)=\rho_0ならば、P(x)=\rho_0 xなので M=\rho_0 l
また、f(x)=\frac{1}{M}\rho_{0}xとなるのでF(x)=\frac{1}{2M}\rho_{0}x^2 となり、G=F(l)-F(0)=\frac{l}{2}
例;\rho(x)=xならば、P(x)=\frac{x^2}{2}なので M=\frac{l^2}{2}
このときf(x)=\frac{1}{M}\rho(x)x=\frac{1}{M}x^2なので
F(x)=\frac{1}{3M}l^3である。G=F(l)-F(0)=\frac{2l}{3}

 2次元以上の物体の慣性モーメントについて 

てこの原理と力のモーメント

図のように剛体の棒の中間に支点Oがあり、
この点をとおり、図面に垂直な軸の周りを自由に回転する装置を梃子(てこ)と呼ぶ。

てこの原理

梃子の端A_1に力\vec{f^1}が作用し、他端A_2に力\vec{f^2}が作用して、
つりあう(静止し続ける)とき、2つの力の間にはどのような関係があるだろうか。
棒は軽くて無視できるとして考察する。
軸周りに静止し続けるということは、
固定軸まわりの運動方程式(1.4.3.5節)から、
梃子に働く外力\vec{f^1}, \vec{f^2}の、回転軸まわり回転力が零であることを意味する。
Oを原点、回転軸をz軸,梃子の棒をx軸とする、直交座標系O-xyzを導入すると、
\vec{OA_1}=(-l_1,0,0),\quad \vec{OA_2}=(l_2,0,0)
\vec{f^1}=({f^1}_x,{f^1}_y,{f^1}_z),\quad \vec{f^2}=({f^2}_x,{f^2}_y,{f^2}_z)
と表現できる。
そこで、1.4.3.2.3節(z軸まわりの回転力の導出)から
z軸まわりのトルク(回転力)は T_{\vec{e_z}}=-l_{1}{f^1}_y+l_{2}{f^2}_y
となる。
従って
つりあい条件は、
l_{1}{f^1}_y=l_{2}{f^2}_y
これをてこの原理という。
l_{2}l_{1}に比べて、非常に大きくとれば、
少しの力{f^2}_yで非常に大きな力{f^1}_yと釣り合わせることが出来ることが分かる。
てこの原理については、

も参照のこと。

剛体に働く力の作用線

力が作用する点を着力点といい、
着力点を通り力のベクトルと方向が等しい直線を、力の作用線という。
剛体に働く力は、その着力点をかえると、一般には、剛体の運動への効果が異なってしまう。
しかし、力のベクトル和と、力のモーメント和が不変となるように力の着力点を移動したり力の合成をすることは、
剛体の運動には全く影響がでないので、許される。
例えば、力の着力点をその作用線にそってうごかしたり、
同じ着力点をもつ複数の力を、それらのベクトル和に置き換えることは許される。。

剛体のつり合い

いくつかの力が作用し、剛体が静止したままであるか、
重心Gが等速直線運動(静止も含む)を続け、
重心の周りの回転が変化しない(回転しないままか、同じ回転を続ける)場合に、
剛体(に作用している力)は釣り合っているという。
重心が等速直線運動を行うのは、
剛体に作用する外力のベクトル和が0になることであり、その場合に限る。
これについては、「1.1.1  質点系の運動と重心」で説明した。
重心周りの回転が変化しないのは、重心まわりの外力のモーメントの総和が0になることであり、この場合に限る。これについては、「1.2.3.5 固定軸の周りの剛体の回転運動の方程式」で説明した。
定理;剛体のつり合い
剛体に、外力\vec{F^1},\vec{F_2},,,,\vec{F_n}がはたらいている。
このとき、次の条件は同等である。
ⅰ)剛体は釣り合っている。
ⅱ)外力のベクトル和が零で、重心Gまわりの外力のモーメントの和が零。
ⅲ)外力のベクトル和が零で、任意の固定点Pまわりの外力のモーメントの和が零。
証明;条件ⅰ)とⅱ)が同等であることは、すでに、説明した。
条件ⅱ)とⅲ)の同等性を示そう。
外力の和が零であるという条件の下で、
「任意の固定点Pまわりの外力のモーメントの和\vec N_{P}は常に等しい」
ことを示せば良い。
外力\vec{F^i}の作用点を\vec {P_i}(i=1,2,\cdots n)とする。
すると、
Pまわりの外力のモーメントの和\vec N_{P}
\vec N_{P}=\sum_{i=1}^{n}\vec{PP_i}\times \vec{F^i} \qquad \qquad (1)

任意の点Qまわりの外力のモーメントの和\vec N_{Q}
\vec N_{Q}=\sum_{i=1}^{n}\vec{QP_i}\times \vec{F^i} \qquad \qquad (2)

\vec{PP_i}=\vec{PQ}+\vec{QP_i}を(1)式に代入すると
\vec N_{P}=\sum_{i=1}^{n}\vec{PP_i}\times \vec{F^i} =\sum_{i=1}^{n}(\vec{PQ}+\vec{QP_i})\times \vec{F^i}
ベクトル積の性質から、
=\sum_{i=1}^{n}(\vec{PQ}\times \vec{F^i}+\vec{QP_i}\times \vec{F^i}) =\vec{PQ}\times \sum_{i=1}^{n}\vec{F^i}+\sum_{i=1}^{n}\vec{QP_i}\times \vec{F^i}
仮定と(2)式から、
=\vec{PQ}\times 0 + \vec N_{Q}=\vec N_{Q}
故に、\vec N_{P}=\vec N_{Q}
証明終わり。  

気体と液体の圧力

この節では気体や液体を、
分子や原子という粒子から構成されるという微視的立場でなく、
巨視的に捉え空間的に滑らかな連続体であるとみなす。
連続体の内部の微小部分に働く力を考え、其の釣合いについて考え、
圧力の性質を導く。

 気体や液体とは何か。

 気体と液体の特徴

気体と液体は体積の変化には抵抗するが、
形の変化には、抵抗しない。(ただし非常に速い変化には抵抗する)。
但し、気体の体積変化への抵抗は小さく、液体は非常に大きい。

 静止気体と液体の圧力

気体や液体は、その表面または内部に任意の面を考えると、その面で2分される部分は、
互いに他を押している。それらは大きさ・方向は等しく、逆向きである(作用反作用の法則)
単位面積当たりのこの力を応力とよぶ。
その発生は、重力の存在と前述の気体や液体の特徴(形の変化に抵抗しない)に起因する。
この力の性質を、気体・液体の特徴から導こう。

応力は面に垂直に働く 

説明は便宜上、液体の語で述べる。
命題1:
静止した液体(気体)の表面あるいは内部に任意のなめらかな面(注参照)を考える。
この面上の応力は、常にこの面に直角に働く。
面と常に直角に働く応力を、圧力と呼ぶ。
(注)面のどの一点においても、その点にごく近い面の部分だけをみれば、平面とみなせる曲面のこと。 理由;
もし、ある面上のある一点Pの周辺の微小面部分(Sと書く)で、押し合う力がこの面と平行な成分を持つとする。
Sは仮定より、平面(の一部)と考えてよい。

図のように、面部分Sとそれと平行な平面の一部S’から作られる、
非常に薄い液体の板状部分Vを考える。

するとVがSを通して液体から受ける力の総和\vec F_Sは、面Sと平行な成分をもつ。
面SとS’は、非常に近いので、
Sを挟んで押し合う力と、S’を挟んで押し合う力は、単位面積当たり、ほぼ等しいと考えてよい。
すると、VがS’を通して液体から受ける力\vec F_{S'}は、
Sを通して受ける力と大きさと方向はほぼ同じで、逆向きになる。
\vec F_{S'}の面Sと平行な成分も、\vec F_SのSと平行な成分と大きさはおなじで、逆向きになる。
液体は自由に形を変えられるので、VのS面とS’面は逆方向に動いてしまい、
静水という条件に反してしまう。
従って、
「ある面上のある一点Pの周辺の微小面部分Sで、押し合う力がこの面と平行な成分を持つ」
という仮定はあり得ないことが示された。

命題2
どの面にも直角に働く応力(圧力)は、どの点でも面の方向によらず一定の強さ(大きさ)をもつ。
証明;
液体中の任意の点をOとする。
Oを原点とする、直交右手系O-xyzを定める。
Oを通る任意の面Hをとる。
O点における、
この面における圧力pとxy平面における圧力p_z、yz平面、zx平面における圧力p_x,p_y
が等しいことを示そう。
平面Hと平行でO点の近くを通る平面H'
x軸、y軸、z軸と交わる点をそれぞれ、
A(\alpha a,0,0),B(0,\alpha b,0),C(0,0,\alpha c)とおく。図参照。
四面体OABCの外部の液体が、
\triangle{OBC}を押す力を\vec F^x,\triangle{OCA}を押す力を\vec F^y,\triangle{OAB}を押す力を\vec F^z,\triangle{ACB}を押す力を\vec F
とおく。
四面体内の液体が静止しているので、
\vec F^x+\vec F^y +\vec F^z+\vec F=0 \qquad (1)
が成り立つ。
この式を圧力で表示しよう。
\lim_{\alpha \to 0}\frac{\|\vec F^x\|}{|\triangle{OBC}|}=p_xなので、
\alpha が十分小さければ
\|\vec F^x\|=|\triangle{OBC}| p_x=\frac{1}{2}|\alpha b \alpha c|p_x
故に、2\vec F^x=\vec{OB}\times \vec{OC}p_x=\alpha b \alpha c p_x\vec{e_x}
同様に2\vec F^y=\vec{OC}\times \vec{OA}p_y=\alpha c \alpha a p_y\vec{e_y}
2vec F^z=\vec{OA}\times \vec{OB}p_z=\alpha a \alpha b p_z\vec{e_z}
2\vec F=\vec{AC}\times \vec{AB}p=p(-\alpha a,0,\alpha c)\times (-\alpha a,\alpha b,0)
これらを(1)式に代入して
p_x\vec{e_x}+p_y\vec{e_y}+p_z\vec{e_z}+\vec{AC}\times \vec{AB}p=0 \quad (2)
これを計算すると、
\left({\alpha}^{2}bc(p_x-p),{\alpha}^{2}ca(p_y-p),{\alpha}^{2}ab(p_z-p)\right)=0
これより、p=p_x=p_y=p_z   証明終わり。

命題3
ⅰ)一様な重力のもとで静止している気体・液体内では、同一水平面上での圧力の大きさは一定である。 
ⅱ)もし液体の密度\rhoが圧力によって変化しないならば、 
深さl_1の水平面H_1上の圧力p_1と  
深さl_2 \quad(l_2>l_1)の水平面H_2上の圧力p_2には  
次の関係が成り立つ。 
p_2=p_1+\rho g(l_2-l_1)  

図示した液体部分Vが静止しているので、Vに作用する力の総和が零になっている。
このことから、この命題は容易に証明できる。

命題4 アルキメデスの原理


気体の圧力と大気圧

気体は圧力が増すと縮むので、命題3のⅱ)の結論は成立しない。
大気は静止していると仮定し、地表の大気圧から高度zでの大気圧を求めてみよう。 地表の一点を原点とし、鉛直上方をz軸の正方向になる座標O-xyzをいれる。
図のように、下底面が高さz、上底面が高さz+hの、単位断面積の角柱Vを考える。
その部分の気体が受ける力の和は零となるので、
次式が成り立つ。
p(z+h)+mg=p(z) \qquad \qquad (1)
ここで
p(z)は高さzの地点の大気圧(命題3のⅰ)から、高度が同じ水平面上で圧力は一定)、
mVの質量。Vの体積hと平均質量密度\rhoの積。
圧力が大きいと空気は縮み質量密度は高くなるので、両者の関係を求めねばならない。
空気体積の変動にともなう温度変化がないとすると、
ボイルの法則(3章1節 熱とエネルギー参照)から、
p\frac{V}{m}=ccは温度だけに依存する数)
質量密度\rho=\frac{m}{V}を代入すると、
\frac{p}{\rho}=c,ゆえに、\rho=\frac{p}{c}
\frac{1}{c}を、cとおくと、
\rho=cp  \qquad \qquad (2)
この質量密度と圧力の関係を用いると、
m=h\rho \approx hcp(z)(hが小さいほど差は少なくなる)
この式を(1)式に代入して、
p(z+h)+cgp(z)h\approx p(z) 、変形すると 
\frac{p(z+h)-p(z)}{h} \approx -cgp(z) 。これより
\frac{dp(z)}{dz}=\lim_{h\to 0}\frac{p(z+h)-p(z)}{h}= -cgp(z)
を得る。これを積分して
p(z)=p_{0}e^{-cgz}
を得る。
ここでp_{0}は、地表での圧力、eネイピア数である。
地表での質量密度が\rho_{0}ならば,(2)式から、
c=\frac{\rho_{0}}{p_{0}}

圧力の単位

圧力は、単位面積当たりの力なので、その単位は面積の単位m^2と力の単位Nから得られる。
Pa=N/m^2=kg\cdot m^{-1}\cdot s^{-2}
が圧力の単位で、パスカルと呼ばれる。

個人用ツール