1

高校数学で空気抵抗ありの落下運動
山本明＠物理のかぎプロジェクト

2005-08-05

日本の高校で行われている物理の授業では，物体の落下については空気抵抗を無視する場合だけ扱っています．発展的な話題として空気抵抗が取り上げられたとしても，落下し始めてからたくさん時間が経った後の速度を求めるくらいしかやりません．*1
だけど空気抵抗があるときの落下運動というのは，Newton も著書「プリンキピア」に書いたくらい歴史ある問題です．
そして高校で習う数学の知識だけでも，速さに比例する空気抵抗の問題を解くことは可能なのです．せっかくならば、たくさん時間が経過した後の速度を求めるだけでなく，時間とともに速度がどんな具合に増えていくのかってことも求めてみましょう．
ちなみに，同じ問題を大学では 微分方程式を使って解きます ．本質的にやっていることは一緒なのですが，微分方程式を使う方が機械的にできて，楽かもしれません．ここでは，高校で習う数学の知識だけで求める…というお遊びです．*2

状況設定
落下している物体に働く空気抵抗は，いろいろなパターンを考えることができます．だけど，解析しやすいものといえば、速さの 1 乗に比例する空気抵抗です．そして速さに比例する空気抵抗は，物体の速さが遅いときには，実験結果とよく一致することが知られています．
そこで、ここでも速さの 1 乗に比例する空気抵抗を考えたいと思います．比例定数を k (k > 0)と書くことにしましょう．空気抵抗の大きさ |f | は，|f | = k|v|ということになります． v と書いたのは物体の速度です．
位置も速度も力の向きも，すべて下方向を正の方向だと決めましょう．反対に上方向を正だと決めてもいいのですが，これは単に趣味の問題．このページ内では，下向きを正とします．
*1

空気抵抗を考えているときに落下し始めてから たくさん時間が経った後の速度 のことを「終端速度」と呼んだりします．これを求めることしかしないってのは，文科省の学習指導要領解説に 空気抵抗については定性的な話題に留めること と記載されている影響が強いのでしょう．
*2 なお高校で習う数学の知識といっているのは，「等比数列」の知識と「極限」の操作です．それらに関する解説は，他をご参照ください．

2

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

運動方程式は…
考えている物体の質量を m と書くことにしますと，Newton の 運動方程式 ma = F はma = mg − kvとなります．力 F の部分は重力 mg と空気抵抗 |f | = k|v| を足し算したものです．力の向きに注意しましょう． 重力は下向きだから，正．空気抵抗は速度と逆向きにかかります．だから kv の前にマイナス符号がついています．なお a という記号は物体の加速度を表しています．

たくさん時間が経った後の速度（終端速度）落下し始めたときは，空気抵抗の大きさと重力の大きさは違うはずで、物体は加速か減速をすることになるでしょう．だけど，たくさん時間が経った後には，空気抵抗と重力がちょうどつり合う状態になるはずです．それを式で表すと，たくさん時間が経った後では，

0 = mg − kv∞ということ． v∞ というのは， たくさん (= ∞ ) 時間が経ったという意味を込めて， v に小さく ∞ を付けて書きました．

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

v∞ =mgkですね．物体はこの速度に達したら，それ以上は加速も減速もせず，等速で落下を続けることになるでしょう．(これが，雨粒が驚異的な速さで落下してこないことの説明になります．)
さて，それではこれだけで満足せずに，速度が始めの状態から v = mgk という状態へどう変化し— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —3

高校数学で空気抵抗ありの落下運動ていくのか，その過程を求めてみましょう．

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

ma = mg − kv
運動方程式は，物体がどんな運動を行うのか教えてくれます．この方程式を変形して， v は時間が経つにつれどう変化するのか求めましょう．

加速度って…
運動方程式の左辺に加速度 a が現れています．加速度とは，ある時間が経って速度がどれだけの割合で変化したかを表す数字ですね．つまり，時刻 t = 0 のときの速度を v0 ，時刻 t = ∆t のときの速度を v1 とすると，その間の加速度は，a=v1 − v0∆tとなります．これを時刻 t = 0 での加速度として， a0 と書くことにしましょう．

a0 =v1 − v0∆tということです．さらに同じように，時刻 t = ∆t のときの加速度を a1 ，時刻 t = 2∆t のときの加速度を a2 ，時刻 t = 3∆t のときの加速度を a3 …というように考えていきましょう．
時刻は物体を落下させ始めてから測ることにします． n という文字は自然数を表すこととして，
時刻 t = n∆t での状態まで、考えることにします．

v2 − v1∆tv3 − v2a2 =∆t..
.
vn+1 − vnan =∆ta1 =∆t はかなり小さな値だと考えると，それぞれの時刻で運動方程式 ma = mg − kv は同じ形で成— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —4

高校数学で空気抵抗ありの落下運動立し続けます．ということは，

ma0 = mg − kv0
ma1 = mg − kv1
ma2 = mg − kv2..
.

man = mg − kvn⇓
v1 − v0= mg − kv0∆tv2 − v1m= mg − kv1∆tv3 − v2m= mg − kv2∆t..
.
vn+1 − vnm= mg − kvn∆tmとなります．

運動方程式を解きやすく変形
右辺に mg の部分があると計算を進めにくいので、 Vn = vn −しょう．上の式に vn = Vn +mgを代入してみます．すると，kmgという書き方をしてみま
kV1 − V0= −kV0∆t
V2 − V1m= −kV1∆t
V3 − V2m= −kV2∆t..
.
Vn+1 − Vn= −kVnm∆tm

いくらか，すっきりした形になりました．ここでさらに，両辺に ∆t を掛け算して，分母から ∆tをなくしてしまいましょう．

mV1 − mV0 = −kV0 ∆tmV2 − mV1 = −kV1 ∆tmV3 − mV2 = −kV2 ∆t..
.

mVn+1 − mVn = −kVn ∆t— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —5

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

さらに少し整理すると，

mV1 = (m − k∆t)V0mV2 = (m − k∆t)V1mV3 = (m − k∆t)V2..
.

mVn+1 = (m − k∆t)Vnとなりますね．数列を勉強した方は，この式を見て，等比数列になっていることに気づけるのではないでしょうか．

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

k∆t)V0mkV2 = (1 − ∆t)V1mkV3 = (1 − ∆t)V2m..
.
V1 = (1 −Vn+1 = (1 −k∆t)Vnm

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

V2 = (1 −k∆t)2 V0mということがわかりますね．これをさらに 3 つ目の式に代入してみると，

V3 = (1 −k∆t)3 V0mとなります．さあ，これを最後の式まで繰り返していったらどうなるか…というと，Vn+1 = (1 −k∆t)n+1 V0mとなります．
…って，そもそも求めたかったのは， t = n∆t の時の速度 Vn でしたから， Vn+1 は行き過ぎですね．ひとつ戻って Vn は，Vn = (1 −と求まります．ちなみに vn = Vn +k∆t)n V0mmgとしていたのだから，
kvn = (1 −mgk∆t)n V0 +mk— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —6

高校数学で空気抵抗ありの落下運動となりますね．この段階でも，もう一般の時刻での速度がわかった!と考えることもできます．*3
だけどせっかくなので、もう少し現代風の見慣れた形にしていきましょう．

最後に極限操作を
いま速度を求めたい時刻を t = n∆t としていますね．つまり， ∆t =tということ．ある時刻
nt までを n 等分して，解析しているということです．
その結果として，時刻 t での速度 v(t) = vn は，vn = (1 −k t nmg) V0 +mnkとなるのでした．Newton の運動方程式は， ∆t が小さいときに成立すると考えられているわけです．ですから，ここで n → ∞ の極限をとりましょう．

kt 1 nmg) V0 +mnkkt 1 n mg
= V0 lim (1 −) +n→∞
mnkv(t) = lim (1 −n→∞kt 1 n) って，どこか見覚えがありませんか? こn→∞mn
こで、自然定数 e という書き方を思い出しましょう．次のように定義される極限値を，特別に eですね．ここに現れる極限操作の部分， lim (1 −(自然定数) という文字で書き表すという約束事をしているのでした．
¡1 ¢Ne ≡ lim 1 +N →∞N
この極限値にそっくりの形が出てきていますね． −しょう． n = −kt
N ということです．すると，mlim (1 −n→∞

1kt 1をまとめて， +と思ってしまいまmnNkt 1 n1 kt) = lim (1 + )− m N
N →∞
mn
Nknt
1 N o− m= lim (1 + )N →∞Nk−mt=eとなります．このことを考えると，

v(t) = V0 lim (1 −n→∞k= V0 e− m t +kt 1 n mg) +mnkmgkmgとしていたことも思い出して，k¡mg ¢ − k t mgv(t) = v0 −e m +kkという結果になります．最後に V0 = v0 −という式を得ます．これで、すっかり現代風の書き方になりました．

— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —7

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

mg
図は v0 < mg
k のとき．だんだん落下速度は速くなり， k に近づいていく．

速度は指数関数 e−t で、mgという値に近づいていくのです．
k

なにをやっていたのか
微分方程式を解いて求めた速度の式と同じ式を，微分方程式を解かずに求めることができました．これはすごいことだ!と思うでしょうか?
ミソは n → ∞ という操作にあります．

ここで行った解法では， n → ∞ という操作をしても，等比数列で求めた vn = · · · の式は成立

dv= F というものも，dtvn+1 − vn
= F という関係が ∆t → 0 の極限でも有限の (微小な) 時間間隔で成り立っていた m∆tし続ける…という信念を抱いています．そして，Newton の運動方程式: m成り立ち続けるという信念を用いて，書き表したものなのです．

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

Newton がその著書「プリンキピア」で示した空気抵抗がある場合の解は，だいたいこのような形だったそうです．
— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —