========================================
抵抗力のある落下運動 2
========================================

抵抗力のある落下運動_ では抵抗力が速度に比例する運動を考えました。
そこでは終端速度が $v_{\infty} = \frac{mg}{\kappa}$ となることを学びました。

ここでは抵抗力が速度の二乗に比例する場合(慣性抵抗と呼ばれています)にどのような運動になるかを見ていきます。

.. _抵抗力のある落下運動: http://www12.plala.or.jp/ksp/mechanics/resistdown/

------------------------------
運動方程式
------------------------------

落下運動に限らず、重力下で慣性抵抗を受けながら運動する物体の運動方程式は、次のようになります。

<tex>
m\frac{d\bm{v}}{dt} = m\bm{g}-\kappa |\bm{v}|^2\frac{\bm{v}}{|\bm{v}|} \tag{#def(init)}
</tex>

この記事では話を簡単にするために、鉛直方向の運動のみを扱うことにします。
つまり落下運動または鉛直投げ上げということになります。
このとき (#ref(init)) は、

<tex>
m\frac{dv}{dt} = mg - \kappa v^2 \frac{v}{|v|} \tag{#def(init2)}
</tex>

となります。ここで $m$ は物体の質量、 $g$ は重力加速度、 $\kappa > 0$ は空気抵抗の比例係数になります。

------------------------------
落下運動
------------------------------

落下時の様子を絵に描くと次図のようになります。落下運動なので $v > 0$ で考えます(軸を下向き正に撮っていることに注意！)

.. figure:: co_vsquare_fig01.png

            抵抗のある場合の落下

運動方程式 (#ref(init2)) は $v > 0$ より

<tex>
m\frac{dv}{dt} = mg - \kappa v^2 \ \ (v > 0) \tag{#def(eq01)}
</tex>

となります。抵抗力の符号は $-$ 、つまり抵抗力は上向きに働くことになりますね。

速度の時間変化を求めてみることにしましょう。(#ref(eq01))の両辺を $m$ で割って、式を整理します。

<tex>
\frac{dv}{dt} & = g - \frac{\kappa}{m} v^2 \\
\frac{dv}{dt} & = g(1-k^2 v^2) \ \ \left(\frac{\kappa}{mg} \equiv k^2 \right) 
</tex>
<tex>
\frac{1}{1-k^2v^2} dv = g dt \tag{#def(eq02)}
</tex>

(#ref(eq02))を積分すれば速度変化を求めることができます。
どうすれば積分を実行できるでしょうか。ここでは部分分数分解を利用することにします。

<tex>
\frac{1}{1-k^2v^2} dv = g dt
</tex>
<tex>
\frac{1}{2}\left(\frac{1}{1-kv} + \frac{1}{1+kv}\right) dv = g dt \tag{#def(eq02-2)}
</tex>

両辺を積分します。

<tex>
\int \left( \frac{1}{1-kv} + \frac{1}{1+kv} \right) dv = 2g \int dt
</tex>
<tex>
- \ln |1-kv| + \ln |1+kv| = 2g k t + 2C. \tag{#def(eq03)}
</tex>

ここで $2C$ は積分定数です。 $2C$ と置いたのは後々のためです。

<tex>
\ln \left| \frac{1+kv}{1-kv} \right| & = 2gkt + 2C\\
\left| \frac{1+kv}{1-kv} \right| & = e^{2(gkt+C)} \tag{#def(eq04)}
</tex>

式 (#ref(eq04)) は分母の $1-kv$ の正負によって場合分けが必要です。

計算練習だと思って手を動かしてみましょう。

<tex>
v(t) & = \frac{1}{k} \ \frac{\pm e^{2(gkt+C)} - 1}{1 \pm e^{2(gkt+C)}} \\
     & = \frac{1}{k} \ \frac{\pm e^{gkt+C} - e^{-(gkt+C)}}{\pm e^{(gkt+C)} + e^{-(gkt+C)}} \tag{#def(eq04-2)}
</tex>

ここで $\pm$ は $0 < v < \frac{1}{k}$ のとき $+$ 、 $\frac{1}{k} < v$ のとき $-$ をとります。
定数 $k = \sqrt{\frac{\kappa}{mg}}$ を元に戻してやると、

<tex>
v(t) & = \sqrt{\frac{mg}{\kappa}} \frac{\pm e^{\sqrt{\frac{g\kappa}{m}}t+C} - e^{-(\sqrt{\frac{g\kappa}{m}}t+C)}}{\pm e^{(\sqrt{\frac{g\kappa}{m}}t+C)} + e^{-(\sqrt{\frac{g\kappa}{m}}t+C)}} \tag{#def(eq06)}
</tex>

となります。
式を見やすくするために $\sqrt{\frac{mg}{\kappa}} \equiv v_{\infty}$ 、 $\sqrt{\frac{m}{g\kappa}} \equiv \tau$ と置くことにします。
(#ref(eq06))式を書き直すと、

<tex>
v(t) & = v_{\infty} \frac{\pm e^{\frac{t}{\tau}+C}-e^{-(\frac{t}{\tau}+C)}}{\pm e^{\frac{t}{\tau}+C}+e^{-(\frac{t}{\tau}+C)}}\\
     & = \begin{cases} v_{\infty} \tanh \left( \frac{t}{\tau} + C \right) \ \  &(0 < v < v_{\infty}) \\
                v_{\infty} \left( \tanh \left( \frac{t}{\tau} + C \right) \right)^{-1} \ \ & (v_{\infty} < v)
\end{cases} \tag{#def(eq-rakka-final)}
</tex>

こうして $v(t)$ の時間変化を得ることができました。



初期条件として $v(t=0) = 0$ をとってやることにしましょう。
(#ref(eq-rakka-final)) で $t=0$ 、 $v=0$ としてやると、 $C=0$ が得られます。

したがって、 $v(t=0) = 0$ を初期条件にとったとき、

<tex>
v(t) = v_{\infty} \tanh \left( \frac{t}{\tau} \right) \tag{#def(eq08)}
</tex>

となります。

このときの速度の変化をグラフに書くと次のようになります。

.. figure:: co_vsquare_fig02.png

            速度の変化(落下運動)

速度は時間が経過すると $v_{\infty}$ へと漸近していく様子がわかります。

~~~~~~~~~~~~~~~~~~~~~~~~~~~~
問い 1.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~

抵抗力の大きさ( $\kappa$ )によって、終端速度 $v_{\infty}$ はどう変わるでしょうか。
また終端速度に達するまでの時間スケールを表す緩和時間 $\tau$ はどうなるでしょうか。


~~~~~~~~~~~~~~~~~~~~~~~~~~~~
問い 2.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~

式 (#ref(eq-rakka-final)) で $v(t=0) > v_{\infty}$ とすると、どのような v-t グラフになるでしょうか。



----------------------------
鉛直投げ上げ
----------------------------

おまけとして鉛直投げ上げをした場合の運動について考えてみます。やはり軸を下向き正にとっていることに注意して下さい。投げ上げなので、 $ v < 0 $ の場合を考えることになります。

.. figure:: co_vsquare_fig03.png

            抵抗のある場合の投げ上げ

運動方程式 (#ref(init2)) は $v < 0$ より次のようになります。

<tex>
m\frac{dv}{dt} = mg + \kappa v^2 \ \ (v < 0, \ \kappa > 0) \tag{#def(eq09)}
</tex>

(#ref(eq01))との違いは、抵抗力につく符号だけです。今度は $+$ なので抵抗力は下向きにかかることになります。
(#ref(eq01))と同様にして解いていくことにしましょう。

<tex>
\frac{dv}{dt} = g( 1 + k^2 v^2) \\
\frac{1}{1+k^2v^2} dv = g dt
</tex>

積分しましょう。

<tex>
\int \frac{1}{1+k^2v^2} dv = g \int dt \tag{#def(eq09-2)}
</tex>

左辺の積分について考えましょう。
$kv = \tan \theta \ ( -\frac{\pi}{2} < \theta < 0 )$ と置換すると $dv = \frac{1}{k}\frac{1}{\cos^2 \theta} d\theta$ となりますので、

<tex>
\frac{1}{k} \ \int \frac{1}{1+\tan^2 \theta} \frac{1}{\cos^2 \theta} d\theta = g \int dt
</tex>
<tex>
\int \cos^2 \theta \frac{1}{\cos^2 \theta} d\theta = gk \int dt
</tex>
<tex>
\int d\theta = gk \int dt.
</tex>

積分を実行すると、

<tex>
\theta = gkt + C
</tex>
<tex>
\tan^{-1} \left( kv \right) = gkt + C
</tex>
<tex>
v(t) = \frac{1}{k} \tan \left( gkt + C \right).
</tex>

$C$ は積分定数です。 $-\frac{\pi}{2} < \theta < 0$ でしたから、 $-\frac{\pi}{2} < gkt + C < 0$ です。
先ほど定義した $v_{\infty}$ と $\tau$ を用いて書くと、

<tex>
v(t) & = v_{\infty} \tan \left( \frac{t}{\tau} + C \right) \ \ \ \rm{where} \ -\frac{\pi}{2}  < \frac{t}{\tau}+C < 0 \tag{#def(eq10)}
</tex>

となります。

初期条件として、 $v(t=0) / v_{\infty} = -1$ をとってみましょう。
$C = -\frac{\pi}{4}$ となりますので、(#ref(eq10))は

<tex>
v(t) & = v_{\infty} \tan \left( \frac{t}{\tau} - \frac{\pi}{4} \right) \tag{#def(eq11)}
</tex>

となります。
$\frac{t}{\tau} = \frac{\pi}{4}$ で速度が $0$ となり、あとは上で考えた落下運動へと移行します。

この様子をグラフにすると、次のようになります。赤線が速度変化を表しています。

.. figure:: co_vsquare_fig04.png

            速度の変化(速度が 0 になると、最初に考えた落下運動へと移行する)


---------------------
別解
---------------------

「落下運動」のセクションでは部分分数分解を用いて積分を、「鉛直投げ上げ」では置換積分を行いました。
積分の形は下のように $\pm$ が違うだけです。

<tex>
\int \frac{1}{1\pm k^2v^2} dv = g \int dt
</tex>

部分分数分解による方法、または置換積分による方法、どちらかだけで解けないものでしょうか。
そのほうが解き方を覚えるのも楽ですよね。

~~~~~~~~~~~~~~~~~~~
落下運動
~~~~~~~~~~~~~~~~~~~

<tex>
\int \frac{1}{1 - k^2 v^2} dv = g \int dt \tag{#def(eq-rakkaomake-1)}
</tex>

まず、落下運動を置換積分で解けないか考えてみます。
結果は(#ref(eq08))のようになることがすでに分かっていて、 $\tanh x$ が出てくるのでした。

そういえば $\sinh x$ 、 $\cosh x$ には $\cosh^2 x - \sinh^2 x = 1$ という関係があり、三角関数とよく似ています。
注目すべきは、両辺を $\cosh^2 x$ で割れば、 $1-\tanh^2 x = \frac{1}{\cosh^2 x}$ という関係が得られることです。
$kv = \tanh \theta$ と置換してやると、うまく行きそうな気になってきませんか？やってみましょう。

と、ここで注意が必要です。
$-1 < \tanh x < 1$ なので、全ての $0 < v$ にたいして $kv = \tanh \theta$ と置換するわけにはいきません。
$0 < kv < 1$ と $1 < kv$ で場合分けが必要です。

我々は落下運動を既に解いて、結果が (#ref(eq-rakka-final)) となることを知っています。なので $0 < v < \frac{1}{k} \ (=v_{\infty})$ では $kv = \tanh \theta$ 、 $v_{\infty} < v$ では $kv = \frac{1}{\tanh \theta}$ と置いてみることにします。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
$0 < v < v_{\infty}$ の場合
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

$dv = \frac{1}{k} \ \frac{1}{\cosh^2 \theta} \ d\theta$ となります。
(#ref(eq-rakkaomake-1)) は、

<tex>
\int \frac{1}{1-\tanh^2 \theta} \ \frac{1}{\cosh^2 \theta} d\theta = kg \int dt
</tex>
<tex>
\int d\theta = gk \int dt
</tex>

となります。積分を実行すると

<tex>
\theta = kgt + C
</tex>

となります。 $\theta$ を元に戻すと

<tex>
\tanh^{-1} (kv) = kgt + C
</tex>
<tex>
v(t) & = \frac{1}{k} \tanh \left( kgt + C \right)\\
     & = v_{\infty} \tanh \left( \frac{t}{\tau} + C \right) \tag{#def(eq-rakka-case1)}
</tex>

となりました。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
$v_{\infty} < v$ の場合
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

$dv = - \frac{1}{k} \ \frac{1}{\tanh^2 \theta} \ \frac{1}{\cosh^2 \theta} \ d\theta$ となります。
(#ref(eq-rakkaomake-1)) は、

<tex>
\int \frac{1}{1-\frac{1}{\tanh^2 \theta}} \ \frac{-1}{\tanh^2 \theta} \ \frac{1}{\cosh^2 \theta} d\theta = kg \int dt
</tex>
<tex>
\int \frac{\tanh^2 \theta}{1-\tanh^2 \theta} \ \frac{1}{\tanh^2 \theta} \ \frac{1}{\cosh^2 \theta} d\theta = kg \int dt
</tex>
<tex>
\int d\theta = gk \int dt
</tex>

となります。積分を実行すると

<tex>
\theta = kgt + C
</tex>

となります。 $\theta$ を元に戻すと

<tex>
\tanh^{-1} \left(\left(kv\right)^{-1}\right) = kgt + C
</tex>
<tex>
v(t) & = \frac{1}{k} \left( \tanh \left( kgt + C \right) \right)^{-1}\\
     & = v_{\infty} \left( \tanh \left( \frac{t}{\tau} + C \right) \right)^{-1} \tag{#def(eq-rakka-case2)}
</tex>

となりました。

式 (#ref(eq-rakka-case1))、(#ref(eq-rakka-case2)) の結果を合わせると、

<tex>
v(t) & = \begin{cases} v_{\infty} \tanh \left( \frac{t}{\tau} + C \right) \ \  &(0 < v < v_{\infty}) \\
                v_{\infty} \left( \tanh \left( \frac{t}{\tau} + C \right) \right)^{-1} \ \ & (v_{\infty} < v)
\end{cases}
</tex>

となり、(#ref(eq-rakka-final)) と一致しました！

~~~~~~~~~~~~~~~~~~~
鉛直投げ上げ
~~~~~~~~~~~~~~~~~~~

では鉛直投げ上げの場合を部分分数分解を用いて積分できるでしょうか。
やってみましょう。
複素数を用いて、無理矢理にでも部分分数分解してやると

<tex>
\int \frac{1}{1 + k^2 v^2} dv = g \int dt \tag{#def(eq-nageomake-1)}
</tex>
<tex>
\frac{1}{2} \int \left( \frac{1}{1+ikv} + \frac{1}{1-ikv} \right) dv = g \int dt \tag{#def(eq-nageomake-2)}
</tex>

となります。積分すると

<tex>
\ln \left( \frac{1+ikv}{1-ikv} \right) = 2i ( gkt + C)
</tex>

となります。ここで  $C$ は積分定数です。 $v$ について整理してやると

<tex>
v(t) & = \frac{1}{ik} \ \frac{e^{igkt+C} - e^{-(igkt+C)}}{e^{(igkt+C)} + e^{-(igkt+C)}} .
</tex>

$\sin x = \frac{e^{ix} - e^{-ix}}{2i}$ 、 $\cos x = \frac{e^{ix} + e^{-ix}}{2}$ の関係を用いてやれば

<tex>
v(t) & = \frac{1}{k} \tan \left( gkt + C \right)
</tex>

が得られます。 $v_{\infty}$ 、 $\tau$ を用いて書き換えると、

<tex>
v(t) & = v_{\infty} \tan \left( \frac{t}{\tau} + C \right)
</tex>

となり (#ref(eq10)) と一致しました！

@@author: CO@@
@@accept: 2006-06-20@@
@@category: 力学@@
@@id: vsquare@@