==============================
高校数学で空気抵抗ありの落下運動
==============================

日本の高校で行われている物理の授業では、物体の落下については空気抵抗を無視する場合だけ扱っています。
発展的な話題としても、空気抵抗がある場合の終端速度を求めるというだけです。

だけど高校で習う数学の知識だけでも、速度に比例する空気抵抗の問題を解くことは可能です。
せっかくならば、終端速度(落下後たくさん時間が経った後の速度)を求めるだけでなく、時間とともに速度がどんな具合に増えていくのか、求めてみましょう。

ちなみに、同じ問題を大学では微分方程式を使って解きます。
本質的にやっていることは一緒なのですが、微分方程式を使う方が機械的にできて、楽かもしれません。
ここでは、高校数学の範囲だけで求める…というお遊びです。 [*]_ 

.. [*] なお高校で習う数学の知識といっているのは、「等比数列」の知識と「極限」の操作です。それらに関する解説は、他をご参照ください。



状況設定
------------------------------

落下している物体に働く空気抵抗は、いろいろなパターンを考えることができます。
だけど、解析しやすいものといえば、速度の1乗に比例する空気抵抗です。そして速度に比例する空気抵抗は、物体の速度が遅いときには、実験結果とよく一致することが知られています。

そこで、ここでも速度の1乗に比例する空気抵抗を考えたいと思います。
比例定数を $k\,(k>0)$ と書くことにしましょう。力の大きさ $|F|$ は、

<tex>
|F|=k|v|
</tex>

ということになります。 $v$ と書いたのは物体の速度です。

位置も速度も力の向きも、すべて下方向を正の方向だと決めましょう。反対に上方向を正だと決めてもいいのですが、これは単に趣味の問題。
このページ内では、下向きを正とします。

.. image:: yamamoto-rakka-fig2.png

運動方程式は…
^^^^^^^^^^^^^^^^
考えている物体の質量を $m$ と書くことにしますと、Newtonの運動方程式 $ma=F$ は

<tex>
ma=mg-kv
</tex>

となります。力の向きに注意しましょう。 $mg$ というのは重力です。
空気抵抗は速度と逆向きにかかります。だから $kv$ の前にマイナス符号がついています。なお $a$ という記号は物体の加速度を表しています。


速度の時間変化を求める
------------------------
状況設定に従って、運動方程式を求めました。

<tex>
ma=mg-kv
</tex>

運動方程式は、物体がどんな運動を行うのか教えてくれます。この方程式を変形して、 $v$ は時間が経つにつれどう変化するのか求めましょう。

加速度って…
^^^^^^^^^^^^^^
運動方程式の左辺に加速度 $a$ が現れています。
加速度とは、ある時間が経って速度がどれだけの割合で変化したかを表す数字ですね。
つまり、時刻 $t=0$ のときの速度を $v_0$ 、時刻 $t=\Delta t$ のときの速度を $v_1$ とすると、その間の加速度は、

<tex>
a=\frac{v_1-v_0}{\Delta t}
</tex>

となります。これを時刻 $t=0$ での加速度として、 $a_0$ と書くことにしましょう。

<tex>
a_0  = \frac{v_1-v_0}{\Delta t}
</tex>

ということです。
さらに同じように、時刻 $t=\Delta t$ のときの加速度を $a_1$ 、時刻 $t=2\Delta t$ のときの加速度を $a_2$ 、時刻 $t=3\Delta t$ のときの加速度を $a_3$ …というように考えていきましょう。

時刻は物体を落下させ始めてから測ることにします。 $n$ という文字は自然数を表すこととして、時刻 $t=n\Delta t$ での状態まで、考えることにします。

<tex>
a_1 &= \frac{v_2-v_1}{\Delta t}
\\
a_2 &= \frac{v_3-v_2}{\Delta t}
\\
    & \vdots
\\
a_n &= \frac{v_{n+1}-v_n}{\Delta t}
</tex>

$\Delta t$ はかなり小さな値だと考えると、それぞれの時刻で運動方程式 $ma=mg-kv$ は同じ形で成立し続けます。ということは、

<tex>
ma_0 &= mg-kv_0
\\
ma_1 &= mg-kv_1
\\
ma_2 &= mg-kv_2
\\
    & \vdots
\\
ma_n &= mg-kv_n
\\
\\
&\Downarrow
\\
\\
m \frac{v_1-v_0}{\Delta t}
 &= mg-kv_0
\\
m \frac{v_2-v_1}{\Delta t}
 &= mg-kv_1
\\
m \frac{v_3-v_2}{\Delta t}
 &= mg-kv_2
\\
    & \vdots
\\
m \frac{v_{n+1}-v_n}{\Delta t}
 &= mg-kv_n
</tex>

となります。

運動方程式を解きやすく変形
^^^^^^^^^^^^^^^^^^^^^^^^^^^
右辺に $mg$ の部分があると計算を進めにくいので、 $\displaystyle V_n=v_n-\frac{mg}{k}$ という書き方をしてみましょう。
上の式に $\displaystyle v_n=V_n+\frac{mg}{k}$ を代入してみます。すると、

<tex>
m \frac{V_1-V_0}{\Delta t}
 &= -kV_0
\\
m \frac{V_2-V_1}{\Delta t}
 &= -kV_1
\\
m \frac{V_3-V_2}{\Delta t}
 &= -kV_2
\\
    & \vdots
\\
m \frac{V_{n+1}-V_n}{\Delta t}
 &= -kV_n
</tex>

いくらか、すっきりした形になりました。
ここでさらに、両辺に $\Delta t$ を掛け算して、
分母から $\Delta t$ をなくしてしまいましょう。

<tex>
m {V_1}-m{V_0} &= -k{V_0}\Delta t
\\
m {V_2}-m{V_1} &= -kV_1{\Delta t}
\\
m {V_3}-m{V_2} &= -kV_2{\Delta t}
\\
    & \vdots
\\
m {V_{n+1}}-m{V_n} &= -kV_n{\Delta t}
</tex>

さらに少し整理すると、

<tex>
m {V_1} &= (m-k\Delta t) V_0
\\
m {V_2} &= (m-k{\Delta t})V_1
\\
m {V_3} &= (m-k{\Delta t})V_2
\\ 
    & \vdots
\\
m {V_{n+1}} &= (m-k{\Delta t})V_n
</tex>

となりますね。数列を勉強した方は、この式を見て、等比数列になっていることに気づけるのではないでしょうか。

等比数列の一般項を求める
^^^^^^^^^^^^^^^^^^^^^^^^^^
いま求めたばかりの式、両辺を $m$ で割ります。

<tex>
 {V_1} &= (1-\frac{k}{m}\Delta t) V_0
\\
 {V_2} &= (1-\frac{k}{m}{\Delta t})V_1
\\
 {V_3} &= (1-\frac{k}{m}{\Delta t})V_2
\\
    & \vdots
\\
 {V_{n+1}} &= (1-\frac{k}{m}{\Delta t})V_n
</tex>

これらの式で、1つ目の式を2つ目の式に代入してみましょう。すると、

<tex>
 {V_2} = (1-\frac{k}{m}{\Delta t})^2 V_0
</tex>

ということがわかりますね。
これをさらに3つ目の式に代入してみると、

<tex>
 {V_3} = (1-\frac{k}{m}{\Delta t})^3 V_0
</tex>

となります。
さあ、これを最後の式まで繰り返していったらどうなるか…というと、

<tex>
 V_{n+1} = (1-\frac{k}{m}{\Delta t})^{n+1} V_0
</tex>

となります。

…って、そもそも求めたかったのは、 $t=n\Delta t$ の時の速度 $V_n$ でしたから、 $V_{n+1}$ は行き過ぎですね。ひとつ戻って $V_n$ は、

<tex>
 V_n = (1-\frac{k}{m}{\Delta t})^n V_0
</tex>

と求まります。
ちなみに $\displaystyle v_n=V_n+\frac{mg}{k}$ としていたのだから、

<tex>
 v_n = (1-\frac{k}{m}{\Delta t})^n V_0 +\frac{mg}{k}
</tex>

となりますね。
この段階でも、もう一般の時刻での速度がわかった!と考えることもできます。 [*]_ 

だけどせっかくなので、もう少し現代風の見慣れた形にしていきましょう。

.. [*] Newtonがその著書「プリンキピア」で示した空気抵抗がある場合の解は、だいたいこのような形だったそうです。

最後に極限操作を
^^^^^^^^^^^^^^^^^
いま速度を求めたい時刻を $t=n\Delta t$ としていますね。
つまり、 $\displaystyle \Delta t=\frac{t}{n}$ ということ。
ある時刻 $t$ までを $n$ 等分して、解析しているということです。

その結果として、時刻 $t$ での速度 $v(t)=v_n$ は、

<tex>
 v_n = (1-\frac{k}{m}\frac{t}{n})^n V_0 +\frac{mg}{k}
</tex>

となるのでした。
Newtonの運動方程式は、 $\Delta t$ が小さいときに成立すると考えられているわけです。ですから、ここで $n\rightarrow\infty$ の極限をとりましょう。

<tex>
 v(t)
 &=\lim_{n\rightarrow\infty}
 (1-\frac{kt}{m}\frac{1}{n})^n V_0 +\frac{mg}{k}
\\
 &=V_0\lim_{n\rightarrow\infty}
 (1-\frac{kt}{m}\frac{1}{n})^n  +\frac{mg}{k}
\\
 &=V_0 e^{-\frac{k}{m}t} +\frac{mg}{k}
</tex>

となります。最後に $\displaystyle V_0=v_0-\frac{mg}{k}$ としていたことも思い出して、

<tex>
 v(t)
 =\bigl(v_0-\frac{mg}{k}\bigr) e^{-\frac{k}{m}t} +\frac{mg}{k}
</tex>

という式を得ます。これで、すっかり現代風の書き方になりました。

速度は指数関数 $e^{-t}$ で、 $\displaystyle \frac{mg}{k}$ という値に近づいていきます。

.. figure:: yamamoto-rakka-fig1.png

 図は $v_0<\frac{mg}{k}$ のとき。
 だんだん落下速度は速くなり、 $\frac{mg}{k}$ に近づいていく。


なにをやっていたのか
---------------------
微分方程式を解いて求めた速度の式と同じ式を、微分方程式を解かずに求めることができました。これはすごいことだ!と思うでしょうか?

ミソは $n\rightarrow\infty$ という操作にあります。

ここで行った解法では、 $n\rightarrow\infty$ という操作をしても、等比数列で求めた $v_n=\cdots$ の式は成立し続ける…という信念を抱いています。
そして、Newtonの運動方程式: $\displaystyle m\frac{dv}{dt}=F$ というものも、有限の(微小な)時間間隔で成り立っていた $\displaystyle m\frac{v_{n+1}-v_{n}}{\Delta t}=F$ という関係が $\Delta t\rightarrow 0$ の極限でも成り立ち続けるという信念を用いて、書き表したものなのです。

ここで行った解法も微分方程式を用いる解法も、同じ信念に基づいているので、どちらも同じ形になるのです。



@@reference: 山本義隆, 古典力学の形成 --ニュートンからラグランジュへ, 日本評論社, 1997, 57-64, 4535782431@@


@@author: 山本明@@
@@accept: 2005-08-05@@
@@category: 力学@@