1-4-3. はねかえり係数

Coefficient of Restitution

はねかえり係数の定義と基本性質

定義

はねかえり係数(反発係数)

ある直線上を運動する物体Aの速度が vA\boldsymbol{v}_{\mathrm{A}} となった瞬間に,同じ直線上を運動する速度 vB\boldsymbol{v}_{\mathrm{B}} の物体Bに衝突したとする。その直後の物体A,Bの速度がそれぞれ vA\boldsymbol{v}'_{\mathrm{A}}vB\boldsymbol{v}'_{\mathrm{B}} となったとき,物体Aと物体Bの間の衝突の性質を表すはねかえり係数 ee は次のように定義される。 ただし,空気抵抗などの影響は考えなくてよいものとする。

e=vBvAvBvA=vBvAvBvAe = \dfrac{|\boldsymbol{v}'_{\mathrm{B}} - \boldsymbol{v}'_{\mathrm{A}}|}{|\boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}}|} = -\dfrac{\boldsymbol{v}'_{\mathrm{B}} - \boldsymbol{v}'_{\mathrm{A}}}{\boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}}}
説明

はねかえり係数は,衝突の前後で2つの物体が互いに「近づく速さ」と「遠ざかる速さ」の比を表す無次元量である。

分母の vBvA\boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}} は衝突前にAから見たBの相対速度であり,物体が衝突するためには互いに近づいている必要がある。一方,分子の vBvA\boldsymbol{v}'_{\mathrm{B}} - \boldsymbol{v}'_{\mathrm{A}} は衝突後の相対速度であり,衝突後は互いに遠ざかるため,必ず衝突前とは逆向きのベクトルとなる。この向きの反転を補正し,比を正の値として定義するために数式にはマイナスの符号がつけられている。

公式

壁に衝突する物体

ある直線上を運動する小球の速度が v\boldsymbol{v} となった瞬間,静止している壁に垂直に衝突し,その直後の小球の速度が v\boldsymbol{v}' となった。このとき,小球と壁の間のはねかえり係数 ee は次のように求める。 ただし,空気抵抗は考えなくてよいものとする。

e=vve = \dfrac{|\boldsymbol{v}'|}{|\boldsymbol{v}|}
導出

静止している壁を一方の物体とみなしてはねかえり係数の定義式を適用する。 壁は衝突の前後を通じて静止し続けているため,壁の速度は常に 0\mathbf{0} である。これを定義式に代入すると,

e=v0v0=vve = \dfrac{|\boldsymbol{v}' - \mathbf{0}|}{|\boldsymbol{v} - \mathbf{0}|} = \dfrac{|\boldsymbol{v}'|}{|\boldsymbol{v}|}

となり,壁などの固定された巨大な物体との衝突においては,単なる衝突前後の速度の大きさ(速さ)の比としてはねかえり係数を計算できることが導かれる。

公式

床と斜めに衝突する物体

水平右向きに xx 軸,鉛直上向きに yy 軸をとる。小球の速度が v1=(v1x,v1y)\boldsymbol{v}_1 = (v_{1x}, v_{1y}) となった瞬間になめらかな水平な床に衝突したとき,その直後の小球の速度 v2\boldsymbol{v}_2 は次のように求める。 ただし,小球と床とのはねかえり係数を ee とし,空気抵抗は考えなくてよいものとする。

v2=(v1x,ev1y)\boldsymbol{v}_2 = (v_{1x}, -e v_{1y})
導出

小球が床に衝突する瞬間において,床から小球に対してはねかえす力(垂直抗力およびその力積)がはたらく。 小球に対して運動方程式を立て,成分ごとに分割して考察する。床からの力を FF,小球の質量を mm とすると,

{mdvxdt=0x 軸方向)mdvydt=Fmgy 軸方向)\begin{cases} m\dfrac{dv_x}{dt} = 0 & \text{($x$ 軸方向)} \\ m\dfrac{dv_y}{dt} = F - mg & \text{($y$ 軸方向)} \end{cases}

水平方向( xx 成分)には床からの力がはたらかない(床はなめらかであると仮定されているため摩擦力は生じない)ため,加速度はゼロである。したがって,xx 方向の速度は衝突前後で変化せず等速運動を保つ。ゆえに衝突後の xx 成分を v2xv_{2x} とすると,v2x=v1xv_{2x} = v_{1x} が成り立つ。

鉛直方向( yy 成分)については,床からの非常に大きな力 FF によって速度が反転する。床の yy 軸方向の速度は常にゼロであるから,はねかえり係数の定義式を yy 成分にのみ適用することができる。

e=0v2y0v1y=v2yv1ye = \dfrac{|0 - v_{2y}|}{|0 - v_{1y}|} = -\dfrac{v_{2y}}{v_{1y}}

これを v2yv_{2y} について解くと,v2y=ev1yv_{2y} = -e v_{1y} となる。 水平成分は維持され,鉛直成分のみが e-e 倍されるため,衝突後の速度ベクトルは v2=(v1x,ev1y)\boldsymbol{v}_2 = (v_{1x}, -e v_{1y}) となることが導出される。

衝突後の速度と換算質量

定義

換算質量

nn 個の質量 m1,m2,,mnm_1, m_2, \dots, m_n からなる系において,その相互作用を単一の物体の運動として取り扱う際に現れる有効的な質量を換算質量 μ\mu といい,次の式を満たすものとして定義される。

1μ=1m1+1m2++1mn\dfrac{1}{\mu} = \dfrac{1}{m_1} + \dfrac{1}{m_2} + \dots + \dfrac{1}{m_n}
導出

物理学において,2つの物体が互いに力を及ぼし合いながら運動する問題(いわゆる二体問題)を解析する際,それぞれの物体の絶対的な運動を追跡するよりも,一方から見た他方の「相対運動」として1つの方程式にまとめる方が見通しが良くなることが多い。 この相対運動の方程式を立てたとき,本来の質量の代わりに現れる係数が換算質量である。2物体の場合は上の定義式を変形し,和分の積の形として次のように表される。

μ=m1m2m1+m2\mu = \dfrac{m_1 m_2}{m_1 + m_2}

換算質量を用いることで,複雑な二体問題を「質量 μ\mu の1つの物体が,固定された中心のまわりを運動する」という単純な一体問題に置き換えることが可能になる。

公式

衝突後の物体の速度

ある直線上を運動する質量 mAm_{\mathrm{A}} の物体Aの速度が vA\boldsymbol{v}_{\mathrm{A}} となった瞬間に,同じ直線上を運動する速度 vB\boldsymbol{v}_{\mathrm{B}},質量 mBm_{\mathrm{B}} の物体Bに衝突した。 物体Aと物体Bとのはねかえり係数を ee とし,空気抵抗は考えないものとする。


衝突直後の物体Aの速度

衝突直後の物体Aの速度 vA\boldsymbol{v}'_{\mathrm{A}} は次のように求める。

vA=mAvA+mBvBmBe(vAvB)mA+mB\boldsymbol{v}'_{\mathrm{A}} = \dfrac{m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}} - m_{\mathrm{B}}e(\boldsymbol{v}_{\mathrm{A}} - \boldsymbol{v}_{\mathrm{B}})}{m_{\mathrm{A}} + m_{\mathrm{B}}}
衝突直後の物体Bの速度

衝突直後の物体Bの速度 vB\boldsymbol{v}'_{\mathrm{B}} は次のように求める。

vB=mAvA+mBvBmAe(vBvA)mA+mB\boldsymbol{v}'_{\mathrm{B}} = \dfrac{m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}} - m_{\mathrm{A}}e(\boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}})}{m_{\mathrm{A}} + m_{\mathrm{B}}}
導出

2つの物体の衝突において,外部からの力がはたらかない閉じた系であるとみなせるため,運動量保存則が成り立つ。

mAvA+mBvB=mAvA+mBvBm_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}} = m_{\mathrm{A}}\boldsymbol{v}'_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}'_{\mathrm{B}}

またはねかえり係数の定義式より,衝突後の相対速度を衝突前の相対速度を用いて表すことができる。

e=vBvAvBvA    vBvA=e(vBvA)=e(vAvB)e = -\dfrac{\boldsymbol{v}'_{\mathrm{B}} - \boldsymbol{v}'_{\mathrm{A}}}{\boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}}} \implies \boldsymbol{v}'_{\mathrm{B}} - \boldsymbol{v}'_{\mathrm{A}} = -e(\boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}}) = e(\boldsymbol{v}_{\mathrm{A}} - \boldsymbol{v}_{\mathrm{B}})

この関係式を vB\boldsymbol{v}'_{\mathrm{B}} について解くと,vB=vA+e(vAvB)\boldsymbol{v}'_{\mathrm{B}} = \boldsymbol{v}'_{\mathrm{A}} + e(\boldsymbol{v}_{\mathrm{A}} - \boldsymbol{v}_{\mathrm{B}}) となる。これを運動量保存則の式に代入して,衝突後の物体Bの速度を消去する。

mAvA+mBvB=mAvA+mB{vA+e(vAvB)}m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}} = m_{\mathrm{A}}\boldsymbol{v}'_{\mathrm{A}} + m_{\mathrm{B}} \left\{ \boldsymbol{v}'_{\mathrm{A}} + e(\boldsymbol{v}_{\mathrm{A}} - \boldsymbol{v}_{\mathrm{B}}) \right\}

右辺を展開して vA\boldsymbol{v}'_{\mathrm{A}} でくくる。

mAvA+mBvB=(mA+mB)vA+mBe(vAvB)m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}} = (m_{\mathrm{A}} + m_{\mathrm{B}})\boldsymbol{v}'_{\mathrm{A}} + m_{\mathrm{B}}e(\boldsymbol{v}_{\mathrm{A}} - \boldsymbol{v}_{\mathrm{B}})

これを vA\boldsymbol{v}'_{\mathrm{A}} について解くことで,物体Aの衝突後の速度が導出される。

vA=mAvA+mBvBmBe(vAvB)mA+mB\boldsymbol{v}'_{\mathrm{A}} = \dfrac{m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}} - m_{\mathrm{B}}e(\boldsymbol{v}_{\mathrm{A}} - \boldsymbol{v}_{\mathrm{B}})}{m_{\mathrm{A}} + m_{\mathrm{B}}}

同様に,はねかえり係数の式を vA=vBe(vAvB)=vB+e(vBvA)\boldsymbol{v}'_{\mathrm{A}} = \boldsymbol{v}'_{\mathrm{B}} - e(\boldsymbol{v}_{\mathrm{A}} - \boldsymbol{v}_{\mathrm{B}}) = \boldsymbol{v}'_{\mathrm{B}} + e(\boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}}) と変形して運動量保存則の式に代入することで,物体Bの衝突後の速度も対称的な形として導出される。

vB=mAvA+mBvBmAe(vBvA)mA+mB\boldsymbol{v}'_{\mathrm{B}} = \dfrac{m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}} - m_{\mathrm{A}}e(\boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}})}{m_{\mathrm{A}} + m_{\mathrm{B}}}

衝突とエネルギーの損失

公式

衝突における運動エネルギーの減少

ある直線上で,質量 mAm_{\mathrm{A}} の物体Aが,質量 mBm_{\mathrm{B}} の物体Bに衝突した。衝突直前に物体Aから見た物体Bの相対速度を vABv_{\mathrm{AB}},はねかえり係数を ee とする。 このとき,衝突によって失われた系全体の運動エネルギー ΔE\Delta E は,換算質量 μ=mAmBmA+mB\mu = \dfrac{m_{\mathrm{A}}m_{\mathrm{B}}}{m_{\mathrm{A}} + m_{\mathrm{B}}} を用いて次のように求めることができる。

ΔE=12μ(1e2)vAB2\Delta E = \dfrac{1}{2}\mu (1 - e^2){v_{\mathrm{AB}}}^2
はねかえり係数の範囲

物理的に運動エネルギーが衝突によって増加することはないため,はねかえり係数 ee は常に次の不等式を満たす。

0e10 \le e \le 1
導出

物体A,Bの衝突前の速度をそれぞれ vA\boldsymbol{v}_{\mathrm{A}}, vB\boldsymbol{v}_{\mathrm{B}},系全体の重心の速度を vG\boldsymbol{v}_G とする。重心速度の定義より,

vG=mAvA+mBvBmA+mB\boldsymbol{v}_G = \dfrac{m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}}}{m_{\mathrm{A}} + m_{\mathrm{B}}}

衝突の直前と直後では外力がはたらかないため運動量保存則が成り立ち,重心速度 vG\boldsymbol{v}_G は衝突前後で一定(不変)である。 ここで,静止している観測者からではなく,重心とともに動く観測者から見た座標系(重心系)を考える。重心系における物体A,Bの速度をそれぞれ uA\boldsymbol{u}_{\mathrm{A}}, uB\boldsymbol{u}_{\mathrm{B}} とすると,これらは絶対速度から重心速度を引くことで得られる。

{uA=vAvG=vAmAvA+mBvBmA+mB=mB(vAvB)mA+mBuB=vBvG=vBmAvA+mBvBmA+mB=mA(vBvA)mA+mB\begin{cases} \boldsymbol{u}_{\mathrm{A}} = \boldsymbol{v}_{\mathrm{A}} - \boldsymbol{v}_G = \boldsymbol{v}_{\mathrm{A}} - \dfrac{m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}}}{m_{\mathrm{A}} + m_{\mathrm{B}}} = \dfrac{m_{\mathrm{B}}(\boldsymbol{v}_{\mathrm{A}} - \boldsymbol{v}_{\mathrm{B}})}{m_{\mathrm{A}} + m_{\mathrm{B}}} \\ \boldsymbol{u}_{\mathrm{B}} = \boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_G = \boldsymbol{v}_{\mathrm{B}} - \dfrac{m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}}}{m_{\mathrm{A}} + m_{\mathrm{B}}} = \dfrac{m_{\mathrm{A}}(\boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}})}{m_{\mathrm{A}} + m_{\mathrm{B}}} \end{cases}

相対速度ベクトルを vAB=vBvA\boldsymbol{v}_{\mathrm{AB}} = \boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}} とおくと,上の結果は uA=mBvABmA+mB\boldsymbol{u}_{\mathrm{A}} = -\dfrac{m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{AB}}}{m_{\mathrm{A}} + m_{\mathrm{B}}} および uB=mAvABmA+mB\boldsymbol{u}_{\mathrm{B}} = \dfrac{m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{AB}}}{m_{\mathrm{A}} + m_{\mathrm{B}}} と表される。

次に,衝突前の系全体の運動エネルギー KK を計算する。総質量を M=mA+mBM = m_{\mathrm{A}} + m_{\mathrm{B}} とおき,絶対速度 vA=vG+uA\boldsymbol{v}_{\mathrm{A}} = \boldsymbol{v}_G + \boldsymbol{u}_{\mathrm{A}}vB=vG+uB\boldsymbol{v}_{\mathrm{B}} = \boldsymbol{v}_G + \boldsymbol{u}_{\mathrm{B}} を運動エネルギーの式に代入して展開する。

K=12mAvA2+12mBvB2=12mAvG+uA2+12mBvG+uB2=12(mA+mB)vG2+(mAuA+mBuB)vG+12mAuA2+12mBuB2\begin{aligned} K &= \dfrac{1}{2} m_{\mathrm{A}} |\boldsymbol{v}_{\mathrm{A}}|^2 + \dfrac{1}{2} m_{\mathrm{B}} |\boldsymbol{v}_{\mathrm{B}}|^2 \\ &= \dfrac{1}{2} m_{\mathrm{A}} |\boldsymbol{v}_G + \boldsymbol{u}_{\mathrm{A}}|^2 + \dfrac{1}{2} m_{\mathrm{B}} |\boldsymbol{v}_G + \boldsymbol{u}_{\mathrm{B}}|^2 \\ &= \dfrac{1}{2} (m_{\mathrm{A}} + m_{\mathrm{B}}) |\boldsymbol{v}_G|^2 + (m_{\mathrm{A}}\boldsymbol{u}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{u}_{\mathrm{B}}) \cdot \boldsymbol{v}_G + \dfrac{1}{2} m_{\mathrm{A}} |\boldsymbol{u}_{\mathrm{A}}|^2 + \dfrac{1}{2} m_{\mathrm{B}} |\boldsymbol{u}_{\mathrm{B}}|^2 \end{aligned}

ここで,第2項の括弧内 mAuA+mBuBm_{\mathrm{A}}\boldsymbol{u}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{u}_{\mathrm{B}} は「重心系における全運動量」を表す。重心系から見れば質量の中心は常に原点に静止しているため,この全運動量はゼロになる。 残った第3項と第4項に先ほど求めた uA\boldsymbol{u}_{\mathrm{A}}uB\boldsymbol{u}_{\mathrm{B}} を代入して整理する。

K=12MvG2+12mAmBvABmA+mB2+12mBmAvABmA+mB2=12MvG2+12(mAmB2+mA2mB(mA+mB)2)vAB2=12MvG2+12(mAmB(mA+mB)(mA+mB)2)vAB2=12MvG2+12μvAB2\begin{aligned} K &= \dfrac{1}{2} M |\boldsymbol{v}_G|^2 + \dfrac{1}{2} m_{\mathrm{A}} \left| -\dfrac{m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{AB}}}{m_{\mathrm{A}} + m_{\mathrm{B}}} \right|^2 + \dfrac{1}{2} m_{\mathrm{B}} \left| \dfrac{m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{AB}}}{m_{\mathrm{A}} + m_{\mathrm{B}}} \right|^2 \\ &= \dfrac{1}{2} M |\boldsymbol{v}_G|^2 + \dfrac{1}{2} \left( \dfrac{m_{\mathrm{A}}{m_{\mathrm{B}}}^2 + {m_{\mathrm{A}}}^2 m_{\mathrm{B}}}{(m_{\mathrm{A}} + m_{\mathrm{B}})^2} \right) |\boldsymbol{v}_{\mathrm{AB}}|^2 \\ &= \dfrac{1}{2} M |\boldsymbol{v}_G|^2 + \dfrac{1}{2} \left( \dfrac{m_{\mathrm{A}} m_{\mathrm{B}}(m_{\mathrm{A}} + m_{\mathrm{B}})}{(m_{\mathrm{A}} + m_{\mathrm{B}})^2} \right) |\boldsymbol{v}_{\mathrm{AB}}|^2 \\ &= \dfrac{1}{2} M |\boldsymbol{v}_G|^2 + \dfrac{1}{2} \mu |\boldsymbol{v}_{\mathrm{AB}}|^2 \end{aligned}

この結果(ケーニヒの定理)は,系全体の運動エネルギーが「重心自体の並進運動エネルギー 12MvG2\frac{1}{2} M |\boldsymbol{v}_G|^2 」と「重心から見た相対運動エネルギー 12μvAB2\frac{1}{2} \mu |\boldsymbol{v}_{\mathrm{AB}}|^2 」の和に分離できることを示している。

衝突後の系の運動エネルギー KK' についても全く同様の計算が成り立つ。外力がはたらかないため重心速度 vG\boldsymbol{v}_G は不変であり,変化するのは相対速度のみである。はねかえり係数の定義より,衝突後の相対速度の大きさは衝突前の ee 倍になる( vAB=evAB|\boldsymbol{v}'_{\mathrm{AB}}| = e |\boldsymbol{v}_{\mathrm{AB}}| )。 したがって,衝突によって失われたエネルギー ΔE\Delta E は,相対運動エネルギーの差として計算される。

ΔE=KK=(12MvG2+12μvAB2)(12MvG2+12μvAB2)=12μvAB212μ(evAB)2=12μ(1e2)vAB2\begin{aligned} \Delta E &= K - K' \\ &= \left( \dfrac{1}{2} M |\boldsymbol{v}_G|^2 + \dfrac{1}{2} \mu |\boldsymbol{v}_{\mathrm{AB}}|^2 \right) - \left( \dfrac{1}{2} M |\boldsymbol{v}_G|^2 + \dfrac{1}{2} \mu |\boldsymbol{v}'_{\mathrm{AB}}|^2 \right) \\ &= \dfrac{1}{2} \mu |\boldsymbol{v}_{\mathrm{AB}}|^2 - \dfrac{1}{2} \mu (e |\boldsymbol{v}_{\mathrm{AB}}|)^2 \\ &= \dfrac{1}{2} \mu (1 - e^2) |\boldsymbol{v}_{\mathrm{AB}}|^2 \end{aligned}
はねかえり係数の範囲の導出

外部からエネルギーの供給がない通常の衝突において,系全体の運動エネルギーが増加することはない。つまり,失われたエネルギー ΔE\Delta E は必ずゼロ以上( ΔE0\Delta E \ge 0 )でなければならない。 換算質量 μ\mu および相対速度の2乗 vAB2|\boldsymbol{v}_{\mathrm{AB}}|^2 はともに正の値をとるため,ΔE0\Delta E \ge 0 を満たすためには (1e2)0(1 - e^2) \ge 0 でなければならない。

1e20    e211 - e^2 \ge 0 \implies e^2 \le 1

また,はねかえり係数 ee は速さの比として定義されているため,常にゼロ以上の実数である。これらの条件を組み合わせることで,

0e10 \le e \le 1

という,はねかえり係数が取り得る値の絶対的な範囲が数学的に証明される。 とくに e=1e=1 のとき ΔE=0\Delta E = 0 となり,運動エネルギーが完全に保存される衝突(弾性衝突)となる。

落体の衝突と連続バウンド

公式

落体の衝突と撃力近似

ある直線上を鉛直下向きに運動する落体Aの速度が vA\boldsymbol{v}_{\mathrm{A}} となった瞬間に,同じ直線上を運動する速度 vB\boldsymbol{v}_{\mathrm{B}} の落体Bに衝突し,その直後の速度がそれぞれ vA\boldsymbol{v}'_{\mathrm{A}}vB\boldsymbol{v}'_{\mathrm{B}} となった。 このとき,重力の影響があるにも関わらず,はねかえり係数 ee は通常の定義式通りに次のように求めることができる。

e=vBvAvBvAe = \dfrac{|\boldsymbol{v}'_{\mathrm{B}} - \boldsymbol{v}'_{\mathrm{A}}|}{|\boldsymbol{v}_{\mathrm{B}} - \boldsymbol{v}_{\mathrm{A}}|}
導出

空気抵抗は考えず,重力などの外力がはたらく状況での衝突を考える。物体A,Bそれぞれに対して,衝突している最中の運動方程式を立てる。物体Aが物体Bから受けた力を FF,鉛直上向きを正として重力加速度を gg とすると,

{mAdvAdt=FmAg(物体A)mBdvBdt=FmBg(物体B)\begin{cases} m_{\mathrm{A}}\dfrac{d\boldsymbol{v}_{\mathrm{A}}}{dt} = F - m_{\mathrm{A}}g & \text{(物体A)} \\ m_{\mathrm{B}}\dfrac{d\boldsymbol{v}_{\mathrm{B}}}{dt} = -F - m_{\mathrm{B}}g & \text{(物体B)} \end{cases}

となる。衝突にかかる極めて短い時間を Δt\Delta t とし,両辺をこの時間で積分して力積を求める。

{mAvAmAvA=FdtmAgdtmBvBmBvB=FdtmBgdt\begin{cases} m_{\mathrm{A}}\boldsymbol{v}'_{\mathrm{A}} - m_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} = \displaystyle \int F dt - \displaystyle \int m_{\mathrm{A}}g dt \\ m_{\mathrm{B}}\boldsymbol{v}'_{\mathrm{B}} - m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}} = -\displaystyle \int F dt - \displaystyle \int m_{\mathrm{B}}g dt \end{cases}

ここで,衝突の瞬間に物体間で及ぼし合う力(撃力) FF は非常に大きく,それに比べて重力 mgmg は極めて小さい。衝突時間 Δt\Delta t が十分に短ければ,重力が及ぼす力積の項 mgdt\displaystyle \int mg dt は,撃力による力積 Fdt\displaystyle \int F dt に対して無視できるほど小さくなる(撃力近似)。 重力の力積をゼロとみなして2式を足し合わせると,

mAvA+mBvB=mAvA+mBvBm_{\mathrm{A}}\boldsymbol{v}_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}_{\mathrm{B}} = m_{\mathrm{A}}\boldsymbol{v}'_{\mathrm{A}} + m_{\mathrm{B}}\boldsymbol{v}'_{\mathrm{B}}

となり,外力である重力が存在しているにも関わらず,衝突の瞬間においては近似的に運動量保存則が成り立つことがわかる。これと同時に,速度変化が純粋に内部の相互作用のみによって決定されるとみなせるため,はねかえり係数の定義式もそのまま適用することが正当化される。

公式

はねかえり係数と高さ

物体を高さ hh から自由落下させて床に衝突させ,はね上がった最高点の高さを hh' とする。このとき,物体と床とのはねかえり係数 ee は次のように求めることができる。

e=hhe = \sqrt{\dfrac{h'}{h}}
導出

物体の質量を mm,重力加速度を gg とする。空気抵抗を無視し,力学的エネルギー保存則を用いて衝突直前および衝突直後の速度を計算する。 床を高さの基準にとる。

落下を開始してから床に衝突する直前までの運動について,エネルギー保存則を立てる。衝突直前の速度の大きさを vv とすると,

mgh+12m02=mg0+12mv2    v2=2ghmg h + \dfrac{1}{2}m \cdot 0^2 = mg \cdot 0 + \dfrac{1}{2}mv^2 \implies v^2 = 2gh

次に,床に衝突した直後から最高点に達するまでの運動について,エネルギー保存則を立てる。衝突直後の速度の大きさを vv' とすると,最高点では速度がゼロになるため,

12m(v)2+mg0=12m02+mgh    (v)2=2gh\dfrac{1}{2}m(v')^2 + mg \cdot 0 = \dfrac{1}{2}m \cdot 0^2 + mg h' \implies (v')^2 = 2gh'

床(静止した壁)との衝突においてはねかえり係数は e=vve = \dfrac{v'}{v} で定義される。ここで求めた速度の2乗の比をとると,

e2=(v)2v2=2gh2gh=hhe^2 = \dfrac{(v')^2}{v^2} = \dfrac{2gh'}{2gh} = \dfrac{h'}{h}

両辺の平方根をとることで,e=hhe = \sqrt{\dfrac{h'}{h}} が導出される。この式は,はねかえり係数を速度の測定なしに高さの比だけで求めることができる実用的な関係式である。

公式

床と複数回衝突する物体(連続バウンド)

時刻 t=0t=0 に,物体を高さ h0h_0 から自由落下させ,床との衝突を繰り返した。反発係数を ee(ただし 0<e<10 < e < 1),重力加速度を gg とし,空気抵抗は考えないものとする。

nn 回衝突後に物体のはね上がる高さ hnh_n

物体が nn 回目の衝突をした後にはね上がる最高点の高さ hnh_n は次のように求める。

hn=h0e2nh_n = h_0 e^{2n}衝突と衝突の間の滞空時間 tnt_n

高さ hn1h_{n-1} から落下し始めてから,nn 回目の衝突を経て高さ hnh_n に達するまでの時間(1サイクルにかかる時間)を tnt_n とすると,次のように求める。

tn=2h0g(en1+en)t_n = \sqrt{\dfrac{2h_0}{g}} (e^{n-1} + e^n)物体の運動が止まる時刻 TT

バウンドを繰り返し,最終的に物体の運動が完全に停止する時刻 TT は次のように求める。

T=2h0g1+e1eT = \sqrt{\dfrac{2h_0}{g}} \cdot \dfrac{1+e}{1-e}
導出
はね上がる高さ hnh_n の導出

前項の「はねかえり係数と高さの関係」より,任意の kk 回目の衝突において,衝突前後の最高点の高さの比は e=hkhk1e = \sqrt{\dfrac{h_k}{h_{k-1}}} を満たす。これを変形すると,

hk=e2hk1h_k = e^2 h_{k-1}

となる。これは,数列 {hn}\{ h_n \} が初項 h0h_0,公比 e2e^2 の等比数列であることを示している。したがって,nn 回目の衝突後の高さは初項から公比を nn 回掛けたものとなり,

hn=h0(e2)n=h0e2nh_n = h_0 (e^2)^n = h_0 e^{2n}

が導かれる。


滞空時間 tnt_n の導出

ある高さ hh から自由落下して床に達するまでの時間は,等加速度直線運動の公式 h=12gt2h = \dfrac{1}{2}gt^2 より t=2hgt = \sqrt{\dfrac{2h}{g}} と計算される。また,床から上方に打ち出され,最高点 hh に達するまでの時間も対称性から全く同じである。

ここで,高さ h0h_0 から落下し,1回目の衝突を経て高さ h1h_1 まで上昇するのに要する時間 t1t_1 を考える。 落下に要する時間を t1at_{1a},上昇に要する時間を t1bt_{1b} とすると,

t1a=2h0g,t1b=2h1g=2(e2h0)g=e2h0gt_{1a} = \sqrt{\dfrac{2h_0}{g}}, \quad t_{1b} = \sqrt{\dfrac{2h_1}{g}} = \sqrt{\dfrac{2(e^2 h_0)}{g}} = e\sqrt{\dfrac{2h_0}{g}}

したがって,t1t_1 はこれらを足し合わせたものとなる。

t1=t1a+t1b=2h0g+e2h0g=(1+e)2h0gt_1 = t_{1a} + t_{1b} = \sqrt{\dfrac{2h_0}{g}} + e\sqrt{\dfrac{2h_0}{g}} = (1+e)\sqrt{\dfrac{2h_0}{g}}

一般の nn 回目のサイクル(高さ hn1h_{n-1} から落下し,hnh_n まで上昇する時間) tnt_n についても,落下時間は hn1h_{n-1} によって,上昇時間は hnh_n によって決まる。

tn=2hn1g+2hng=2(h0e2(n1))g+2(h0e2n)gt_n = \sqrt{\dfrac{2h_{n-1}}{g}} + \sqrt{\dfrac{2h_n}{g}} = \sqrt{\dfrac{2(h_0 e^{2(n-1)})}{g}} + \sqrt{\dfrac{2(h_0 e^{2n})}{g}}

これを整理すると,

tn=en12h0g+en2h0g=2h0g(en1+en)t_n = e^{n-1} \sqrt{\dfrac{2h_0}{g}} + e^n \sqrt{\dfrac{2h_0}{g}} = \sqrt{\dfrac{2h_0}{g}} (e^{n-1} + e^n)

となり,各サイクルにかかる時間が導出される。これは初項 t1t_1,公比 ee の等比数列をなしている。


停止時刻 TT の導出

物体がバウンドを停止するというのは,無限回の衝突を終えた極限( nn \to \infty )の状態を意味する。このとき 0<e<10 < e < 1 であるため,はね上がる高さの極限は limnhn=0\lim_{n \to \infty} h_n = 0 となり,確かに運動は停止する。 運動が完全に止まるまでの総時間 TT は,すべてのサイクルにかかる時間 tnt_n を無限回足し合わせた無限級数の和として計算される。

T=n=1tn=n=12h0g(en1+en)T = \sum_{n=1}^\infty t_n = \sum_{n=1}^\infty \sqrt{\dfrac{2h_0}{g}} (e^{n-1} + e^n)

この級数は,初項 t1=(1+e)2h0gt_1 = (1+e)\sqrt{\dfrac{2h_0}{g}},公比 ee の無限等比級数である。公比の絶対値が 11 より小さいため,この無限等比級数は収束し,その和の公式 初項1公比\dfrac{\text{初項}}{1 - \text{公比}} を用いることができる。

T=(1+e)2h0g1e=2h0g1+e1eT = \dfrac{(1+e)\sqrt{\dfrac{2h_0}{g}}}{1 - e} = \sqrt{\dfrac{2h_0}{g}} \cdot \dfrac{1+e}{1-e}

これにより,物体は無限回のバウンドを繰り返すにも関わらず,有限の時刻 TT において完全に静止するという結論が数学的に導かれる。