力の等価変換

質点の力学ではすべての運動はニュートンの運動の3法則で計算できると説明されていますが, 剛体の力学ではある点の回りのモーメントの計算が必要になります. そのためには,まずモーメントの計算では力が束縛ベクトルである(大きさと向きが同じでも 作用点によって働きが異なる)ことを再認識しなければなりません. ここでは,作用点を明示した束縛ベクトルの表現と計算例を示します.

等価な力とは?

簡単な例として,質量 m_1, m_2 の質点を軽い棒でつないだ下図の剛体を考えます. \bm F_i, \bm F_{ik} はそれぞれ質量 m_i の質点に働く外力,質量 m_k の質点から 質量 m_i の質点に働く内力です.

pulsar-BoundVector-Fig1.gif

2質点をつないだ剛体

質量 m_1, m_2 の質点の位置ベクトルをそれぞれ \bm r_1, \bm r_2 とすると,内力は

\bm F_{12} + \bm F_{21} = \bm 0
(\bm r_1 - \bm r_2) \times \bm F_{12} = \bm 0

という性質をもっているので,運動方程式

m_i \frac{\mathrm d^2 \bm r_i}{\mathrm d t^2}= \bm F_i + \sum_{k \neq i}\bm F_{ik}

から

m_1 \frac{\mathrm d^2 \bm r_1}{\mathrm d t^2} + m_2 \frac{\mathrm d^2 \bm r_2}{\mathrm d t^2} = \bm F_1 + \bm F_2
\bm r_1 \times m_1 \frac{\mathrm d^2 \bm r_1}{\mathrm d t^2} + \bm r_2 \times m_2 \frac{\mathrm d^2 \bm r_2}{\mathrm d t^2} + = \bm r_1 \times \bm F_1 + \bm r_2 \times \bm F_2

という内力を含まない式が得られ,この微分方程式を解くことによって \bm r_1, \bm r_2 を求めることができます [*] .剛体の運動を考えるときモーメントを計算するのは, 前記の (\bm r_1 - \bm r_2) \times \bm F_{12} = \bm 0 という制約によるものです.

[*]重心 \bm r_{\mathrm G}(m_1 + m_2) \bm r_{\mathrm G} = m_1 \bm r_1 + m_2 \bm r_2 で定義すると,質点の運動方程式と類似の (m_1 + m_2) \frac{\mathrm d^2 \bm r_{\mathrm G}}{\mathrm d t^2} = \bm F_1 + \bm F_2 が得られます.

この剛体に働く力がつり合っているということは \bm F_1 + \bm F_2 = \bm 0, \bm r_1 \times \bm F_1 + \bm r_2 \times \bm F_2 = \bm 0 であることを意味します. このことから, \bm F_1, \cdots , m, \cdots , F_m, -\bm F_{m+1}, \cdots , -\bm F_n が つり合っているとき,すなわち

\sum_{i=1}^m \bm F_i = \sum_{i=m+1}^n \bm F_i
\sum_{i=1}^m \bm r_i \times \bm F_i = \sum_{i=m+1}^n \bm r_i \times \bm F_i

であるとき, \bm F_1, \cdots , \bm F_m\bm F_{m+1}, \cdots , \bm F_n は 等価であると定義します.ここで \bm r_i\bm F_i の作用点の位置ベクトルを 表しています.このとき,任意の \bm r_0 について

\sum_{i=1}^m (\bm r_i - \bm r_0) \times \bm F_i = \sum_{i=m+1}^n (\bm r_i - \bm r_0) \times \bm F_i

となることを確認してください.

束縛ベクトルの表現

\bm r_i を作用点とする力 \bm F_i\bm F_i(\bm r_i) と表しても, 通常は \bm F_i(\bm r_i) = \bm F_k(\bm r_k)\bm F_i(\bm r_i)\bm F_k(\bm r_k) の大きさと向きが等しいと解釈されます.このため, 新しい記号 \bm F_i[\bm r_i] を用いて

\sum_{i=1}^m \bm F_i[\bm r_i] = \sum_{i=m+1}^n \bm F_i[\bm r_i]

の意味を

\sum_{i=1}^m \bm F_i = \sum_{i=m+1}^n \bm F_i
\sum_{i=1}^m \bm r_i \times \bm F_i = \sum_{i=m+1}^n \bm r_i \times \bm F_i

で定義します.この定義からただちに (\bm r_1 - \bm r_2) \times \bm F_0 ならば \bm F_0[\bm r_1] = \bm F_0[\bm r_2] (力はその作用線上を移動しても働きは変わらない) という性質が導かれます [†]

[†](\bm r - \bm r_1) \times \bm F_0 となる点 \bm r の集合を力 \bm F_0[\bm r_1] の作用線といいます.

二つの力の合力

与えられた \bm F_1[\bm r_1], \bm F_2[\bm r_2] に対して

\bm F_1[\bm r_1] + \bm F_2[\bm r_2] = \bm F_3[\bm r_3]

となる \bm F_3[\bm r_3] が存在するとき, \bm F_3[\bm r_3]\bm F_1[\bm r_1], \bm F_2[\bm r_2] の合力といいます. \bm F_1[\bm r_1], \bm F_2[\bm r_2] が 同一平面上にあれば, \bm F_1 + \bm F_2 = \bm 0, \bm r_1 \neq \bm r_2 の場合を除いて,合力が存在します.

平行でない2力

\bm F_1[\bm r_1], \bm F_2[\bm r_2] が平行でないときは,これらのベクトルは 1次独立ですから \bm r_1 - \bm r_2 = a \bm F_1 - b \bm F_2 となる a, b が 存在し,

\bm F_3 = \bm F_1 + \bm F_2
\bm r_3 = \bm r_1 - a \bm F_1 = \bm r_2 - b \bm F_2

と表せます( \bm r_3 は2力の作用線の交点).

例えば \bm F_1[\bm r_1] = (F, 0)[(2r, 0)], \bm F_2[\bm r_2] = (0, F)[(0, 2r)] のとき,

\left( \begin{array}{cc} 2r & -2r \\ \end{array} \right)= \left( \begin{array}{cc} a & -b \\ \end{array} \right)\left(\begin{array}{cc} F & 0 \\ 0 & F \\ \end{array} \right)

ですから, \bm F_3[\bm r_3] = (F, F)[(0, 0)] となります.ただし, (0, 0) は 剛体上にないので, \bm F_3[\bm r_3] の作用線上にある剛体上の点 \bm r_4=(r, r) を求め, \bm F_3[\bm r_3] と等価な \bm F_3[\bm r_4] = (F, F)[(r, r)] を 合力とする方が自然でしょう [‡]

pulsar-BoundVector-Fig2.gif

(F, 0)[(2r, 0)] + (0, F)[(0, 2r)]

[‡]任意の r_a, r_b, r_c について (F, F)[(0, 0)] = (F, 0)[(r_a, r_b)] + (0, F)[(r_c, r_a)] が成立します.

平行な2力

\bm F_2 = c \bm F_1 \hspace{1zw} (c \neq -1) であれば,

\bm r_1 \times \bm F_1 + \bm r_2 \times \bm F_2= (\bm r_1 + c \bm r_2) \times \bm F_1

ですから,

\bm F_3 = (1 + c)\bm F_1
\bm r_3 = \frac{\bm r_1 + c \bm r_2}{1 + c}

です.

例えば \bm F_1[\bm r_1] = (0, -F)[(0, 0)], \bm F_2[\bm r_2] = (0, -2F)[(3r, 0)] のとき,

\bm r_3 = \frac{(0, 0)+2(3r, 0)}{1+2} = (2r, 0)

で, \bm F_3[\bm r_3] = (0, -3F)[(2r, 0)] となります.

pulsar-BoundVector-Fig3.gif

(0, -F)[(0, 0)] + (0, -2F)[(3r, 0)]

偶力

\bm F_1 + \bm F_2 = \bm 0, \bm r_1 \neq \bm r_2 のとき, \bm F_1[\bm r_1], \bm F_2[\bm r_2] の合力は存在しません.このような力の組 \bm F_1[\bm r_1], \bm F_2[\bm r_2] を偶力といい,以下では \bm F_0[\bm r_1] - \bm F_0[\bm r_2]\bm F_0[\bm r_1, \bm r_2] と略記します.

偶力のモーメントが自由ベクトルであること,すなわち任意の \bm r_0 について

\bm r_1 \times \bm F_0 + \bm r_2 \times (-\bm F_0)= (\bm r_1 - \bm r_0) \times \bm F_0 + (\bm r_1 - \bm r_0) \times (-\bm F_0)

が成立することや,偶力 \bm F_0[\bm r_1, \bm r_2] と 偶力 c^{-1} \bm F_0[c(\bm r_1 - \bm r_2), \bm 0] のモーメントが等しいこと等が 容易に確かめられます.

3次元の力の等価変換

\bm F_1[\bm r_1], \bm F_2[\bm r_2] が同一平面上になければ これらの合力は存在しませんが,任意の \bm r_0 について

\bm F_1[\bm r_1] + \bm F_2[\bm r_2]= \bm F_3[\bm r_0] + \bm F_4[\bm r_4, \bm r_0]

が成立するように \bm F_3, \bm F_4, \bm r_4 を選ぶことができます.

pulsar-BoundVector-Fig4.gif

\bm F_i[\bm r_i] の等価変換

例えば, \bm F_i[\bm r_i]

\bm F_i[\bm r_i] = \bm F_i[\bm r_0] + \bm F_i[\bm r_i] - \bm F_i[\bm r_0]= \bm F_i[\bm r_0] + \bm F_i[\bm r_i, \bm r_0]

と変形して加算した式

\sum_{i=1}^2 \bm F_i[\bm r_i]= \left( \sum_{i=1}^2 \bm F_i \right)[\bm r_0] + \sum_{i=1}^2 \bm F_i[\bm r_i, \bm r_0]

の右辺第2項は偶力の和ですから,自由ベクトルであるこれらのモーメントの和

\bm N = \sum_{i=1}^2 (\bm r_i - \bm r_0) \times \bm F_i

は容易に計算でき, (\bm r_4 - \bm r_0) \times \bm F_4 = \bm N を満足する \bm F_4, \bm r_4 が求められます [§]

[§]\bm N\bm F_3 と平行でないとき, \bm F_3 と平行な偶力の成分を \bm F_3 に加算し(合力が存在します),これと \bm F_3 と直交する偶力 (モーメントは \bm F_3 と平行)との和に変形できます.

上式で \bm r_0 が任意の点でよかったことを思い出すと,多数の3次元ベクトル \bm F_1[\bm r_1], \cdots , \bm F_n[\bm r_n] についても

\sum_{i=1}^n \bm F_i[\bm r_i]= \bm F_{n+1}[\bm r_0] + \bm F_{n+2}[\bm r_{n+2}, \bm r_0]

となる \bm F_{n+1}, \bm F_{n+2}, \bm r_{n+2} が存在することが分かります.

一様な密度の剛体に働く力

重力

\bm F_0[\bm r_0 - \bm r_1] + \bm F_0[\bm r_0 + \bm r_2] = 2 \bm F_0[\bm r_0] ですから,密度が一様で,点 \bm r_0 に関して対称な剛体に働く重力は,重心を 作用点とし,大きさと向きが力の総和に等しい力と等価であることが分かります.

例えば,両端が \bm r_1, \bm r_2 にある質量 M の棒に働く重力は,棒を 2 n 個の微小領域に等分して点対称な位置にある微小領域の対を考えると,各対の合力の 作用点はすべて棒の中心になり,これらの合力の総和が棒に働く重力と等しくなります.

固定点から受ける力

剛体が固定点の回りを運動しているとき,通常は固定点回りの角運動量の時間微分が 固定点回りの外力のモーメントの総和に等しいとして運動方程式を解きます.この場合, 固定点に働いている力を計算する必要はありませんが,上記の説明と関連付けるために 固定点 \bm r_0 に働く力が \bm F_0[\bm r_0], その他の外力(重力は重心に集中 しているとして扱います)が \bm F_i[\bm r_i] (i = 1, \cdots , n) である場合 について考えましょう.

この剛体に働く力は

\bm F[\bm r] = \left( \sum_{i=0}^n \bm F_i \right) [\bm r]+ \sum_{i=0}^n \bm F_i[\bm r_i, \bm r]

と等価ですから, \bm r = \bm r_0 とおくと \bm r_0 回りの角運動量を求める 方程式から \bm F_0 が消えてしまいます.つまり, \bm F_0 を求める必要がなければ 運動量を求める方程式は考えなくてよいのです.

数値例について

「二つの力の合力」では,数値例に \bm F_1[\bm r_1] = (0, -F)[(0, 0)], \bm F_2[\bm r_2] = (0, -2F)[(3r, 0)] のような表現を用いました.この点について 補足します.

通常,物理量を表す文字は単位を含んでいます.例えば,速さ v で時間 t だけ 移動したときの移動距離は vt であるといったとき, vt の単位は自由に選べます. v = 72 \mathrm{km/h}, t = 3 \mathrm{s} ならば vt = 60 \mathrm{m} です. なお,速さ v1 \mathrm{s} 間移動したときの移動距離は v ではなく (1 \mathrm{s})v です [¶] .上記の (0, -2F)[(3r, 0)] のような数値例は (0, -2 \mathrm{N})[(3 \mathrm{m}, 0)] のような表現より分かりやすいと思います.

[¶]0 \mathrm{km/h} = 0 \mathrm{m/s} だから v = 0 の右辺に単位は不要ですが, 温度 T = 0^{\circ} \mathrm C の単位は省略できません. (蛇足ですが,熱容量や比熱に使われる温度差の 1^{\circ} \mathrm C274.15 \mathrm K でなく 1 \mathrm K です. [\mathrm J] でない [\mathrm N \cdot \mathrm m], 無次元でない [\mathrm m / \mathrm m] 等, 単位については苦し紛れの区別がいくつかあります.)

あとがき

ここで述べた内容はすべて私が学生のときに学んだ教科書[1](古書,ISBNなし)に書かれて います.既知の事柄を珍奇な記号を使って書き換えただけだと感じた方が多いと思いますが, モデリング,記法,プログラミング等に関心のある高校生・大学生諸君のために所属学会の 研究会で発表した内容を補足し,数値例も付け加えて紹介しました.教科書・参考書を超えて 自分でいろいろ工夫してほしいというのが筆者の願いです.

参考文献

  1. 国井修二郎,千田香苗著,『 力学?』( 丸善, 1958, 80-86)