2つのボールが衝突する場合の、各ボールの運動を求めたい。
ボール同士の衝突では、ダランベールの原理と弾性衝突条件を再考する必要がある
右図のように、2つのボールの半径をそれぞれ
、中心座標を
とおく。また、両方のボールの物理量をまとめて大文字で書くことにする。例えば、位置と速度はそれぞれ
と書く。知りたいのは、互いに衝突する2つのボールの運動
を計算する方法である。
これまでの議論で扱えないのは、衝突する瞬間のみである。よって新たに知る必要があるのは、前章と同様、「入射速度
から反射速度
を導く衝突公式」である。今回も、最もよく跳ねかえる極限である弾性衝突を考える。衝突の瞬間以外は、個々のボールに対して第2章で述べたニュートンの運動方程式を使って計算できる。
第4章で考えた「壁とボールの衝突」の場合、弾性衝突公式を導くには、2つの条件:速度変化の方向を決める「ダランベールの原理」(=拘束力は拘束面に垂直)と、弾性衝突を表す「弾性衝突条件」(=壁との相対速度の大きさが変化しない)を連立すればよかった。この方法はしかし、今回の問題にはそのままでは使えない。実際、ダランベールの原理は「各ボールについての速度変化の方向」を与えるので、残る未知数は「各ボールの速度変化の大きさ」の2つ。一方、弾性衝突条件は1つだけであり、条件が1つ足りず、
の6成分全てを決めることができない。よって、ダランベールの原理を拡張する必要がありそうである。また、弾性衝突条件についても、どのような形で表されるか今の段階では分からない。
以上により、
を求めるには、衝突での速度の変化を与える弾性衝突公式が必要であり、そのためには、ボール同士の衝突における「ダランベールの原理」と「弾性衝突条件」を導く必要がある。これらを踏まえてこの章では、以下の3つの節に分けて考えることにする:
ボール同士の衝突におけるダランベールの原理ボール同士の衝突における弾性衝突条件ボールの運動の計算
なお、ボールは球対称なものを考える。(一般の形状の物体の衝突では、物体の回転も変化する。そういった場合については、第10章以降で扱う剛体の理論が必要である。)
6.1ボール同士の衝突におけるダランベールの原理
この節では、ボール同士の弾性衝突におけるダランベールの原理()を導く。この式を変形することで、衝突前後の速度に関する条件式()が得られる(未知数
を決めるには、次節の弾性衝突条件が必要である)。
拘束条件()から拘束力 の向きを求めたい
第4章で見たように、ダランベールの原理は、拘束条件
から拘束力
の方向を求めるものである。ボール同士の衝突においても、拘束条件
を書き下しておこう。これは、
間の距離が
以上:
とすればよい。
また、2つのボールがあるので、運動方程式は2本になるが、これらもまとめて以下の形になる:
は衝突時に働く拘束力、
はそれ以外の重力などの力(=外力)である。(行列
のサイズは
である。)
知りたいのは、ダランベールの原理、即ち
拘束条件から拘束力の向きを求める法則
である。
各々のボールに対するダランベールの原理()→だが条件が足りない
まず、ボールが1つの場合のダランベールの原理は第4章で既に述べたので、ボール1に適用してみよう。ダランベールの原理によると、ボール1が受ける拘束力
は、衝突面に垂直(右図)、即ち、拘束条件
を
で偏微分したもの
に比例する:
は未定乗数、
は
を正規化したもの。第2式の最後の変形(絶対値の微分)は、第4章の公式【4.3-注1】を使った。
同様に、ボール2についてダランベールの原理を適用すると、未定乗数を
として
となる。これら
をまとめて書くと、以下のようになる:
この式の未知数は
の2つである。しかし、冒頭で述べた通り、未知数が2つあると、弾性衝突条件と連立した時に未知数が残ってしまい、衝突後の速度が一意的に決まらない。従って、これら
に対する関係式が、他にもまだあるはずである。細かいことを言うと、そもそも衝突後の速度が、質量や速度といった既出の変数だけから決定できない可能性もある。その場合はもちろん未知数が残るのが正しい。今は、未知数が残らないことが実験的に知られている、という想定で話を進めている。
ボール同士の衝突におけるダランベールの原理:式()
2つの未定乗数
の間に条件を課すためには、各ボールに働く拘束力
を関係づける新しい法則が必要である。そのような法則として、よく知られているものがある。作用・反作用の法則である。これは、「接触する2つの物体が互いに作用しあう力は、大きさが等しく向きが逆になる」という法則である(以下の【6.1-注1】)。これにより、拘束力は
を満たすことになる(右図)。従って、式()より
とならなければならない。
よって結局、式()の拘束力
は、1つの未定乗数
を用いて書けることになる:
これが、ボール同士の衝突におけるダランベールの原理である。拘束力が、拘束条件の微分
に比例する形になっているが、これは、ボールが1つだけの場合のダランベールの原理と全く同じ形である。実はこれが一般の拘束条件で成り立つことを、第8章で示す。
ダランベールの原理()が得られたので、衝突前後の速度、それぞれ
に対する関係式の形にしたい。考え方は第4章と同じよである。まず、衝突は一瞬なので、運動方程式()において外力
の影響は無視できる。そのため、
が変化する方向は、拘束力
の方向に一致する。即ち、衝突の前後での速度の関係式は、以下のように書ける:
よってあと1つ条件があれば、未知数
が決まり、衝突後の速度
が求まる。その条件とは、弾性衝突条件となるはずである。なお、作用・反作用の法則は、弾性衝突かどうかとは無関係に常に成り立つので、式()も弾性衝突でなくても成り立つ(ただし横方向の摩擦力が無視できるとする)。
【6.1-注1】作用・反作用の法則:式()
2つの物体
が互いに接触しているとする(右図)。この時、「
が
から受ける接触力
」と「
が
から受ける接触力
」は、大きさが同じで向きが逆になる:
これを作用・反作用の法則という。
補足
作用・反作用の法則により、複合物体の運動方程式から拘束力を消すことができる。例えば右図のように、物体
が、くっついたままそれぞれ外力
を受けながら運動するとする。この時、
それぞれの運動方程式:(
は共通の加速度)
について辺々和をとると、以下のようになる:
作用・反作用の法則()により
部分が消えるので、この式は、質量
を持つ1つの物体に、力
がかかっている場合の運動方程式とみなせる。よって、拘束力は、うまいこと消えて考える必要がなくなるのである。もし、作用・反作用の法則が成り立っていないと、物体の運動を考える際、物体の内部でどのような力がかかっているのかを考えなければならず、面倒なことになる。
6.2ボール同士の衝突における弾性衝突条件
この節では、弾性衝突条件()を導出する。これと、前節で求めた速度の変化の式()を連立すれば、未知数
が決まり、求めたかった弾性衝突公式()が得られる。
重心系 から見ると、速度の大きさの変化は等しい:式()
弾性衝突は、「最もよく跳ね返る」という状況を表すわけであるが、この時、衝突の前後で「運動の勢い」は変化しないはずである。=衝突を繰り返すたびに、両方のボールがどんどん加速していったり減速していったりすることはないボールが1つだけの場合であれば、「運動の勢い」とは速度の大きさであるが、2つのボールがある場合、「運動の勢い」とはどのように定量化されるのだろうか。
そのための取っ掛かりとして、2つのボールが同じ質量である場合を考えてみよう。まず、衝突前の速度の大きさが同じで正面衝突する場合:
直感的に、「弾性衝突後の速度の大きさは変化しない」というのが弾性衝突条件になると予想できる。次に、式()を満たさない場合であるが、この場合であっても、適当な慣性系
の観測者から見ることにより、式()が成り立つようにできる。式()は各成分3つの条件を課しているが、
の取り方も速度3成分の自由度があるので数があっているよって、弾性衝突条件が求められそうである。ということは、質量が異なる場合であっても、上手く慣性系
を選ぶことによって、弾性衝突条件が得られるのではないだろうか。
そこで着目するのが、「作用・反作用の法則が成り立つ系」に存在する保存則である。即ち、重心
(以下の【6.2-注1】の式())の移動速度
:
が、一定となる(同注の式()):
重心速度
は衝突の前後で変わらないので、重心とともに等速度で移動する慣性系
(=重心系)をとることができる。
から見ると、重心位置
は静止してる(右図)。重力などの外力があると式()は成り立たなくなるが、衝突の瞬間の微小時間を考えているので、外力は無視できる。
重心系
から見た重心速度
はゼロを保つ、即ち、衝突の前後での式()の値はそれぞれ:(添え字
は入射速度、
は反射速度を表す)
となる。ここから質量
を消して、速度だけの関係式にすると以下のようになる:
具体的な導出手順:「両式の分母を払う」→「ボール2に対する量を右辺に移項する」→「両辺の絶対値をとって辺々割る」。これにより、重心系
から見れば、衝突による速度の大きさの変化の割合は、両方のボールで等しいことが分かる(ともに
倍になる)。
【6.2-注1】重心
2つのボール
の重心
は、右図のように、
間を
に内分する点として定義される:
もっと粒子数が多い場合にも、容易に拡張できる:
補足
2つの粒子に働く力をそれぞれ
とする。これらが
を満たすならば、重心速度
は一定となる:
【導出】定義式()を2階微分して運動方程式
を代入するだけである:
加速度
がゼロなので、速度
は一定となる。
弾性衝突条件:式()
さて、上式()により、重心系
から見た時、衝突による速度の大きさの変化の比率
は、両方のボールで等しい(右図)。よって、この系において、衝突時にあり得る可能性は、「どちらのボールも衝突前より遅くなる/速くなる/変わらない」のいずれかである。弾性衝突に対応するのは、もちろん「運動の勢い」が変化しないもの、即ち、「変わらない」である(
)。
以上により、重心系
での弾性衝突条件としては、「ボール1の速度の大きさが変わらない」ことを要求すればよい:
(式()が成り立つので、どちらか一方のボールについてだけでよい。)これを本来の慣性系に戻せば、弾性衝突条件が得られる:
導出の手順:「ガリレイ変換
(第5章の【5.1-注2】)の時間微分
を式()に代入する」→「重心の定義式()を用いて
を消去する」。この式()が、求めたかった弾性衝突条件である。
弾性衝突条件()は、相対速度の大きさが衝突の前後で変化しないことを意味している。これは、動く壁との衝突の場合の弾性衝突条件(第5章の【5.2-注2】)と同じ形になっている。壁は「質量が非常に大きい物体」だとみなせるので、今考えている問題の特殊な場合(=一方のボールの質量が非常に大きい場合)に対応しているわけである。
ボール同士の衝突における弾性衝突公式:式()
以上により、「ボール同士の弾性衝突の公式」を得るための条件式がそろった。具体的には、先ほど求めた弾性衝突条件()、および、前節で求めたダランベールの原理に対応する式()を連立して、未知数
を求めればよい。実際に計算すると、以下の【6.2-注2】のようになる。
【6.2-注2】ボール同士の弾性衝突の公式
2つのボールが弾性衝突した時、衝突後の速度
は、衝突前の速度
を用いて、以下のように計算できる:(
は
を正規化したもの)
導出の概要
式()の上側2式は、式()をボールごとに分けて書いたものである。後は、未知数
を求めればよい。式()を式()(を辺々2乗したもの)と連立して、
を消去すれば、単なる
の2次方程式となる。2次方程式なので解の候補は2つあるが、第4章の【4.1-注3】と同様に、一方の解は
となる。よって、もう一方の解が正しい解である。この
を愚直に計算すると、式()の第3式が得られる。
6.3ボールの運動 の計算
ボール同士の弾性衝突の公式()が得られたので、弾性衝突する2つのボールの運動
を計算することができる。この節では、具体的な計算手法についてまとめた後、いくつかの例について数値計算を行う。
の計算方法
弾性衝突する2つのボールの運動
を数値計算によって求める手順は、弾性衝突公式が変わる以外は前章までと同じである。まず、衝突の瞬間(=拘束条件
が成り立つ時刻
)まで、運動方程式
(
は重力などの外力)を用いて計算を進める。その後、衝突点
において、弾性衝突公式()を使って衝突後の速度
を計算する。その後は、その
を初期値として、再び
を用いて、次の衝突まで計算をすすめる。以上を繰り返せばよい。
【例題1】一方のボールが静止している場合
静止したボールに「様々な質量のボール」を衝突させる場合、数値計算を行うと右図のようになる。。
【例題2】重力下での衝突
球の内部で重力を受けて運動する場合、数値計算を行うと右図のようになる。