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

日本の高校で行われている物理の授業では、物体の落下については空気抵抗を無視する場合だけ扱っています。
発展的な話題として空気抵抗が取り上げられたとしても、落下し始めてからたくさん時間が経った後の速度を求めるくらいしかやりません。 [*]_

だけど空気抵抗があるときの落下運動というのは、Newtonも著書「プリンキピア」に書いたくらい歴史ある問題です。

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

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

.. [*] 空気抵抗を考えているときに落下し始めてから“たくさん時間が経った後の速度”のことを「終端速度」と呼んだりします。これを求めることしかしないってのは、文科省の学習指導要領解説に“空気抵抗については定性的な話題に留めること”と記載されている影響が強いのでしょう。

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



状況設定
------------------------------
落下している物体に働く空気抵抗は、いろいろなパターンを考えることができます。
だけど、解析しやすいものといえば、速さの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>

となります。
力 $F$ の部分は重力 $mg$ と空気抵抗 $|f|=k|v|$ を足し算したものです。
力の向きに注意しましょう。 重力は下向きだから、正。
空気抵抗は速度と逆向きにかかります。だから $kv$ の前にマイナス符号がついています。なお $a$ という記号は物体の加速度を表しています。



たくさん時間が経った後の速度（終端速度）
------------------------
落下し始めたときは、空気抵抗の大きさと重力の大きさは違うはずで、物体は加速か減速をすることになるでしょう。
だけど、たくさん時間が経った後には、空気抵抗と重力がちょうどつり合う状態になるはずです。それを式で表すと、たくさん時間が経った後では、
<tex>
0=mg-kv_\infty
</tex>
ということ。 $v_\infty$ というのは、“たくさん(= $\infty$ )時間が経った”という意味を込めて、 $v$ に小さく $\infty$ を付けて書きました。

この式で $kv_\infty$ を左辺に移項して、両辺を $k$ で割り算すれば、たくさん時間が経った後の速度がわかります。

<tex>
v_\infty=\frac{mg}{k}
</tex>
ですね。
物体はこの速度に達したら、それ以上は加速も減速もせず、等速で落下を続けることになるでしょう。
(これが、雨粒が驚異的な速さで落下してこないことの説明になります。)

さて、それではこれだけで満足せずに、速度が始めの状態から $v=\frac{mg}{k}$ という状態へどう変化していくのか、その過程を求めてみましょう。



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

<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}
</tex>

ですね。ここに現れる極限操作の部分、 $\displaystyle \lim_{n\rightarrow\infty} (1-\frac{kt}{m}\frac{1}{n})^n $ って、どこか見覚えがありませんか?
ここで、自然定数 $e$ という書き方を思い出しましょう。
次のように定義される極限値を、特別に $e$ (自然定数)という文字で書き表すという約束事をしているのでした。

<tex>
 e\equiv \lim_{N\rightarrow\infty}\bigl( 1+\frac{1}{N}\bigr)^N
</tex>

この極限値にそっくりの形が出てきていますね。 $\displaystyle -\frac{kt}{m}\frac{1}{n}$ をまとめて、 $\displaystyle +\frac{1}{N}$ と思ってしまいましょう。 $\displaystyle n=-\frac{kt}{m}N$ ということです。すると、

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

となります。このことを考えると、

<tex>
 v(t)
 &=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@@

.. _微分方程式を使って解きます: ../../resistdown/index.html
.. _運動方程式: ../../maf/index.html

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