- 追加された行はこの色です。
- 削除された行はこの色です。
#rst2hooktail_source
============================================================
空気抵抗がある時のモンキーハンティング
============================================================
空気抵抗がある時のモンキーハンティングについて書きます。
銃弾の初期位置を原点 $(x,y)=(0,0)$ に、サルの初期位置を $(X,Y)=(X_0,Y_0)$ とします。
そこには一様に重力 $-g$ がかかっていて、空気抵抗は速度に比例してかかるものとします。
そして、サルの質量を $M$ ,銃弾の質量を $m$ とします。
すると、サルの運動方程式は、
<tex>
M \ddot{X} &= 0 \tag{##} \\
M \ddot{Y} &= -k \dot{Y} -mg \tag{##}
</tex>
速度 $(\dot{X},\dot{Y})$ は、
<tex>
M \dot{X} &= 0 \tag{##}
</tex>
はいいとして、Y方向は、
は良いとして、Y方向は、
<tex>
M \dfrac{d}{dt} ( \dot{Y} + \dfrac{Mg}{k} ) &= -k( \dot{Y} + \dfrac{Mg}{k} ) \tag{##}
</tex>
<tex>
\dfrac{(d/dt)( \dot{Y} + \dfrac{Mg}{k} )}{\dot{Y} + \dfrac{Mg}{k}} &= -\dfrac{k}{M} \tag{##}
</tex>
時間0からtまで定積分すると、Y方向への初期速度は0より、
<tex>
\log \left| \dfrac{\dot{Y} + \dfrac{Mg}{k}}{ 0 + \dfrac{mg}{k}} \right| &= -\dfrac{kt}{M} \tag{##}
\log \left| \dfrac{\dot{Y} + \dfrac{Mg}{k}}{ 0 + \dfrac{Mg}{k}} \right| &= -\dfrac{kt}{M} \tag{##}
</tex>
<tex>
\dot{Y} + \dfrac{Mg}{k} &= \dfrac{Mg}{k} e^{-\dfrac{kt}{M}} \tag{##}
</tex>
<tex>
\dot{Y} &= \dfrac{Mg}{k} e^{-\dfrac{kt}{M}} - \dfrac{Mg}{k} \tag{##}
</tex>
となります。
銃弾の運動方程式は、同様に解くと、初期速度 $(\dot{x},\dot{y})=(\dot{x}_0,\dot{y}_0)$ として、
<tex>
\dot{x} &= \dot{x}_0 e^{-\dfrac{k}{m}t} \tag{##} \\
\dot{y} &= (\dot{y}_0 + \dfrac{mg}{k}) e^{-\dfrac{k}{m}t} - \dfrac{mg}{k} \tag{##}
\dot{x} &= \dot{x}_0 e^{-\dfrac{kt}{m}} \tag{##} \\
\dot{y} &= (\dot{y}_0 + \dfrac{mg}{k}) e^{-\dfrac{kt}{m}} - \dfrac{mg}{k} \tag{##}
</tex>
ここで、落下するサルから銃弾の速度変化を見ます。つまり、
<tex>
\dfrac{\dot{y}-\dot{Y}}{\dot{x}-\dot{X}} = \dfrac{(\dot{y}_0 + \dfrac{mg}{k}) e^{-\dfrac{k}{m}t} - \dfrac{mg}{k} - (\dfrac{Mg}{k} e^{-\dfrac{kt}{M}} - \dfrac{Mg}{k}) }{\dot{x}_0 e^{-\dfrac{k}{m}t}} \tag{##}
\dfrac{\dot{y}-\dot{Y}}{\dot{x}-\dot{X}} = \dfrac{(\dot{y}_0 + \dfrac{mg}{k}) e^{-\dfrac{kt}{m}} - \dfrac{mg}{k} - (\dfrac{Mg}{k} e^{-\dfrac{kt}{M}} - \dfrac{Mg}{k}) }{\dot{x}_0 e^{-\dfrac{kt}{m}}} \tag{##}
</tex>
となり、よく分からない量になってしまいますが、サル、銃弾の両者の質量 $M=m$ が等しいときには、
<tex>
\dfrac{\dot{y}-\dot{Y}}{\dot{x}-\dot{X}} &= \dfrac{(\dot{y}_0 + \dfrac{mg}{k}) e^{-\dfrac{k}{m}t} - \dfrac{mg}{k} - (\dfrac{mg}{k} e^{-\dfrac{kt}{m}} - \dfrac{mg}{k}) }{\dot{x}_0 e^{-\dfrac{k}{m}t}} \\
&= \dfrac{\dot{y}_0 e^{-\dfrac{k}{m}t}}{\dot{x}_0 e^{-\dfrac{k}{m}t}} \\
\dfrac{\dot{y}-\dot{Y}}{\dot{x}-\dot{X}} &= \dfrac{(\dot{y}_0 + \dfrac{mg}{k}) e^{-\dfrac{kt}{m}} - \dfrac{mg}{k} - (\dfrac{mg}{k} e^{-\dfrac{kt}{m}} - \dfrac{mg}{k}) }{\dot{x}_0 e^{-\dfrac{kt}{m}}} \\
&= \dfrac{\dot{y}_0 e^{-kt/m}}{\dot{x}_0 e^{-kt/m}} \\
&= \dfrac{\dot{y}_0}{\dot{x}_0} \\
&= \dfrac{Y_0}{X_0} \tag{##}
</tex>
最後の等式は、銃弾の初期速度がサルの方向を向いていたということを示しています。
つまり、この時銃弾は、サルから見ると真っ直ぐに近づいてくることになり、サルに当たります。
ただし、これは銃弾は初期速度によって、一定値を超えることができないので [*]_ 、
銃弾の $x$ 座標が $X_0$ を超えることができる時です。
それでは、今日はこの辺で。お疲れ様でした。
.. [*] 式 $(9)$ を積分すれば、<tex> x = \dfrac{m\dot{x}_0}{k}(1- e^{-kt/m}) </tex>となるので、
.. [*] 式 $(9)$ を積分すれば、 $ x = \dfrac{m\dot{x}_0}{k}(1- e^{-kt/m}) $ となるので、
つまりこの時、銃弾は $x=\dfrac{m \dot{x}_0}{k}$ を超えられません。
@@author:クロメル@@
@@accept:2013-06-07@@
@@category:力学@@
@@id:monkeyInAir@@