空気抵抗の微分方程式の解き方は？と検索してこのページにたどり着いた方は、おそらく力学の授業や受験勉強、あるいは大学のレポート課題で行き詰まっているのではないでしょうか。
空気抵抗が絡む運動方程式は、重力だけの単純な自由落下と違い、速度に依存する項が加わることで一気に難易度が上がります。
教科書を読んでも、途中の積分計算でつまずいてしまい、結局どこで何をしているのかわからなくなった経験がある方も多いはずです。
特に速度に比例する空気抵抗と速度の2乗に比例する空気抵抗では、解法の道筋がまったく異なります。
片方は指数関数、もう片方は双曲線関数が登場するため、混同してしまうとどちらの計算をしているのか迷子になりがちです。
さらに終端速度という概念も登場し、パラシュートや雨粒の落下運動といった身近な現象を理解するうえで欠かせない知識となります。
この記事では、空気抵抗の微分方程式の解き方について、運動方程式の立て方から丁寧に解説していきます。
変数分離法の使い方や双曲線関数を用いた積分計算、終端速度の求め方まで、順を追って確認していただければ、空気抵抗を含む運動方程式に対する苦手意識はきっと解消されるはずです。
数式だけでなく、パラシュートやスカイダイビングといった具体例も交えながら解説していきますので、最後まで読んでみてください。
空気抵抗の微分方程式の解き方の結論
それでは、空気抵抗の微分方程式の解き方について、結論からお伝えしていきます。
結論として、空気抵抗を含む運動方程式は、抵抗力が速度に比例するか、速度の2乗に比例するかによって解法が大きく変わります。
速度比例モデルでは変数分離法を使い、指数関数を含む解が得られます。
一方でv²比例モデルでは非線形の項が現れるため、双曲線正接関数、いわゆるtanhを使った解が導かれます。
どちらのモデルでも、時間が十分に経過すると速度が一定値に収束していきます。
この一定値こそが終端速度と呼ばれるものです。
結論を先取りしてしまえば、空気抵抗の微分方程式を解く作業とは、抵抗力の形を見極めたうえで正しい積分公式を選び出す作業にほかなりません。
空気抵抗の微分方程式は、抵抗力の形によって解法が変わるという点がもっとも重要なポイントです。
速度比例なら指数関数、速度の2乗比例なら双曲線関数というキーワードをまず押さえておきましょう。
この2つの型さえ見分けられれば、あとは機械的な積分計算で答えにたどり着けます。
速度に比例する空気抵抗の場合の結論
質量mの物体が速度vで運動し、空気抵抗が速度に比例してkvと表せる場合を考えてみましょう。
このときの運動方程式は、m dv/dt イコール mg マイナス kv という形になります。
この微分方程式を変数分離法で解くと、速度vは時間tの関数として指数関数的に終端速度へ近づいていく式になります。
具体的には v イコール (mg/k)(1 マイナス e の マイナスkt/m 乗) という解が得られます。
この式からもわかるように、時間が経つほど指数の項がゼロに近づき、速度はmg/kという一定値に落ち着いていきます。
低速で運動する小さな物体、たとえば粘性の高い液体中を沈む球などは、このモデルでよく近似されます。
速度の2乗に比例する空気抵抗の場合の結論
次に、空気抵抗が速度の2乗に比例してkv²と表せる場合です。
この場合の運動方程式はm dv/dt イコール mg マイナス kv²となり、非線形微分方程式になります。
変数分離法自体は使えるものの、積分の過程で双曲線関数が登場する点が特徴的です。
最終的な解は v イコール √(mg/k) かける tanh(√(gk/m) かける t) という形にまとまります。
tanhは値が1に収束する関数なので、時間が経つと速度は√(mg/k)という終端速度に近づいていきます。
空気中を落下する多くの物体、たとえば人間や雨粒、パラシュートなどは、こちらのモデルで近似されることが一般的です。
終端速度の求め方の結論
終端速度を求めたいだけであれば、実は微分方程式をすべて解く必要はありません。
加速度がゼロになる瞬間、つまりdv/dtがゼロになる条件を使えば、代数的に終端速度を導けます。
速度比例モデルであればmgイコールkvより、終端速度はmg/kです。
v²比例モデルであればmgイコールkv²より、終端速度は√(mg/k)となります。
この方法は計算量が少なく、テストや演習で時間がないときに非常に役立つテクニックです。
ただし途中の速度の変化を知りたい場合は、やはり微分方程式そのものを解く必要があります。
空気抵抗を含む運動方程式の立て方
続いては、空気抵抗を含む運動方程式の立て方について確認していきます。
微分方程式を正しく解くためには、まず物理的な状況を数式に落とし込む力が欠かせません。
ここでつまずいてしまうと、そのあとの積分計算がどれだけ正確でも答えが合わなくなってしまいます。
ニュートンの運動方程式の基本形
すべての出発点はニュートンの第二法則です。
質量mの物体にはたらく力の合計Fは、質量と加速度の積に等しく、F イコール ma と表されます。
加速度aは速度vの時間微分dv/dtで表せるため、F イコール m dv/dt と書き換えられます。
この式が、空気抵抗を含むあらゆる運動方程式の土台になります。
逆にいえば、この基本形さえ理解していれば、抵抗力の種類が変わっても対応できるということです。
重力と空気抵抗の力のつり合い
落下運動を考える場合、物体にはたらく力は主に重力と空気抵抗の2つです。
重力は下向きにmg、空気抵抗は運動方向と逆向き、つまり上向きにはたらきます。
下向きを正とすると、合力はmg マイナス 抵抗力という形で表現できます。
この抵抗力の部分に、速度比例のkvを入れるか、速度の2乗比例のkv²を入れるかで、その後の展開がまったく違うものになります。
どちらの抵抗力を採用するかは、物体の大きさや速さ、流体の性質によって決まる点も覚えておきましょう。
微分方程式として表す方法
力のつり合いの式をF イコール m dv/dt に代入すると、微分方程式が完成します。
速度比例モデルの場合
m dv/dt イコール mg マイナス kv
v²比例モデルの場合
m dv/dt イコール mg マイナス kv²
ここでkは空気抵抗の大きさを決める比例定数で、物体の形状や空気の粘性などによって変わる値です。
この2つの式こそが、これから解いていく空気抵抗の微分方程式の正体です。
式の形さえ正しく立てられれば、あとは数学の積分計算の話になりますので、物理と数学の役割分担を意識しておくと理解が進みやすくなります。
速度に比例する空気抵抗の微分方程式の解き方
続いては、速度に比例する空気抵抗の微分方程式を、実際に手を動かして解いていきます。
変数分離法による解法
m dv/dt イコール mg マイナス kv という式を、変数分離法で解いていきます。
変数分離法とは、vを含む項とtを含む項を左右に分けてから、両辺を積分する手法です。
この式を整理すると、dv/(mg マイナス kv) イコール dt/m という形に変形できます。
左辺はvについて、右辺はtについて積分すればよいだけなので、比較的シンプルに解き進められます。
この段階で符号を間違えるミスが多いので、マイナスの扱いには注意しながら計算を進めましょう。
一般解の導出過程
両辺を積分すると、マイナス(1/k)ln|mg マイナス kv| イコール t/m プラス C という式が得られます。
初期条件として、時刻t イコール0のとき速度v イコール0とすると、積分定数Cを決定できます。
この条件を代入して整理していくと、最終的に v イコール (mg/k)(1 マイナス e の マイナスkt/m 乗) という一般解にたどり着きます。
この式の形を見れば、tがゼロのときにvがゼロになり、tが無限大に近づくとvがmg/kに収束することが一目でわかるはずです。
指数関数の肩にマイナスの符号がつくという点が、この解の大きな特徴です。
具体的な計算例
質量1キログラムの物体が、比例定数k イコール0.2で自由落下するケースを考えてみましょう。
重力加速度をgイコール9.8とすると、終端速度はmg/kイコール49メートル毎秒と計算できます。
この値を先ほどの一般解に代入すれば、任意の時刻における速度を求められます。
たとえばt イコール5のとき、指数部分はe の マイナス1乗となり、速度はおよそ31メートル毎秒程度まで加速していることがわかります。
このように、パラメータを具体的な数値に置き換えることで、微分方程式の解が持つ意味をより実感しやすくなります。
計算問題として出題される際も、まずはこの一般解の形を覚えておくと見通しよく解答できます。
速度の2乗に比例する空気抵抗の微分方程式の解き方
続いては、速度の2乗に比例する空気抵抗の微分方程式について確認していきます。
非線形微分方程式としての特徴
m dv/dt イコール mg マイナス kv² という式は、vの2乗が含まれているため非線形微分方程式に分類されます。
非線形とはいっても、変数分離自体は可能なので、落ち着いて計算すれば十分に解けるタイプの式です。
ただし積分の過程で登場する関数の形が、速度比例モデルとはまったく異なる点に注意が必要です。
ここで焦って速度比例モデルと同じように対数関数で処理しようとすると、計算が破綻してしまいます。
双曲線関数を使った解法
まず終端速度をv∞イコール√(mg/k)と置いておきます。
すると元の式は dv/dt イコール g(1 マイナス v²/v∞²) という形に書き換えられます。
この式を変数分離すると、dv/(1 マイナス v²/v∞²) イコール g dt となります。
左辺の積分は、双曲線正接関数の逆関数であるアークタンハイパボリックの形になることが知られています。
部分分数分解を使って対数関数だけで処理することも可能ですが、双曲線関数を使ったほうが式が簡潔にまとまります。
積分の計算手順
積分を実行すると v∞ かける artanh(v/v∞) イコール gt プラス C という式が得られます。
初期条件v イコール0、t イコール0を代入するとCイコール0となります。
これをvについて解くと、v イコール v∞ かける tanh(gt/v∞) という解にたどり着きます。
tanh関数はtが大きくなるにつれて1に近づく性質を持つため、この解も時間経過とともにv∞、つまり終端速度に収束していきます。
v²比例モデルの解は指数関数ではなく双曲線関数で表されるという点が、v比例モデルとの決定的な違いです。
グラフに描いてみると、速度比例モデルよりも初期の立ち上がりが急で、終端速度に近づいたあとの動きが緩やかになる傾向があります。
終端速度とは何か、その求め方
続いては、ここまで何度も登場してきた終端速度について、もう少し深く確認していきます。
終端速度の定義と物理的意味
終端速度とは、重力と空気抵抗がつり合い、加速度がゼロになったときの速度のことです。
加速度がゼロということは、それ以上速度が増加も減少もしない状態を意味します。
落下し続ける物体は、最終的にこの終端速度に近づいていき、理論上は無限の時間をかけて到達する形になります。
実際にはごく短時間でほぼ終端速度に達するため、日常の感覚としては一定の速さで落ちているように見えます。
このため終端速度は、実用上その物体が最終的にたどり着く速さと考えて差し支えありません。
速度比例モデルにおける終端速度
速度比例モデルでは、dv/dtイコール0を代入すると0イコールg マイナス(k/m)vとなります。
これを整理すると、終端速度はv∞イコールmg/kと求められます。
質量が大きいほど終端速度は大きくなり、比例定数kが大きいほど終端速度は小さくなる関係が見て取れます。
v²比例モデルにおける終端速度
v²比例モデルでも同様に、dv/dtイコール0を代入すると0イコールg マイナス(k/m)v²となります。
これを解くと、終端速度はv∞イコール√(mg/k)です。
速度比例モデルとは異なり、質量やkの平方根が関わってくる点が特徴です。
| モデル | 運動方程式 | 終端速度 | 解の形 |
|---|---|---|---|
| 速度比例モデル | m dv/dt イコール mg マイナス kv | mg/k | 指数関数 |
| 速度の2乗比例モデル | m dv/dt イコール mg マイナス kv² | √(mg/k) | 双曲線関数 |
表にまとめると、2つのモデルの違いがより明確に整理できるのではないでしょうか。
問題を解くときは、まずどちらのモデルが与えられているかを確認し、この表を思い出しながら計算方針を立てるとスムーズです。
落下運動における空気抵抗の応用例
続いては、空気抵抗の微分方程式が実際の落下運動でどのように使われているのか、具体例とあわせて確認していきます。
パラシュートの落下運動
パラシュートは表面積が大きいため、空気抵抗の影響を強く受ける代表的な例です。
パラシュートが開いた瞬間は速度が大きいため、空気抵抗も非常に強くはたらき、急激に減速していきます。
その後は終端速度に近づき、地面に着くころにはほぼ一定の速度で降下する状態になります。
この一連の流れは、v²比例モデルの微分方程式を使うことでかなり正確に説明できます。
パラシュートの設計者は、この終端速度を安全な着地速度に収まるよう、比例定数kにあたる面積や形状を調整しています。
雨粒の落下運動
雨粒も空気抵抗を強く受けながら落下する身近な存在です。
小さな雨粒は速度比例モデルに近い挙動を示すことが多く、比較的低い終端速度で安定します。
一方で大きな雨粒になると空気抵抗の影響がより複雑になり、v²比例モデルで近似されることが一般的です。
雨粒の大きさによって終端速度が異なるため、大粒の雨ほど地面に強く当たるという現象にもつながっています。
スカイダイビングの物理
スカイダイバーが自由落下する際の速度も、空気抵抗の微分方程式で説明できる現象のひとつです。
体を大の字に広げると空気抵抗が大きくなり、終端速度は時速およそ200キロメートル程度に抑えられます。
逆に頭から突っ込むような姿勢を取ると空気抵抗が小さくなり、終端速度はさらに大きくなります。
パラシュートを開くタイミングも、この終端速度に達したあとの安定した状態を利用して計算されています。
スカイダイビングの姿勢による速度の違いは、まさに比例定数kの大小が終端速度を左右する好例です。
空気抵抗の微分方程式を理解していると、こうしたスポーツの物理現象も数式で説明できるようになります。
身近な現象と数式を結びつけて考える習慣が、微分方程式への理解をより深めてくれるはずです。
まとめ
ここまで、空気抵抗の微分方程式の解き方について、運動方程式の立て方から具体的な計算過程まで解説してきました。
速度に比例する空気抵抗の場合は変数分離法によって指数関数を含む解が得られ、速度の2乗に比例する場合は変数分離のあとに双曲線関数を含む解が得られるという違いを押さえておきましょう。
どちらのモデルでも、時間が十分に経過すると速度は終端速度と呼ばれる一定値に収束していきます。
終端速度だけを求めたいのであれば、微分方程式をすべて解かずとも、加速度がゼロになる条件から代数的に導ける点も便利なポイントです。
パラシュートや雨粒、スカイダイビングといった身近な落下運動も、こうした微分方程式によって説明できることがおわかりいただけたのではないでしょうか。
空気抵抗の微分方程式の解き方は？という疑問を持って調べ始めた方も、この記事を通じて運動方程式の立式から解法、終端速度の意味まで、一通りの流れをつかんでいただけたなら幸いです。
実際に手を動かして計算してみることで、公式の暗記ではなく理解として身についていくはずですので、ぜひ紙とペンを使って一度ご自身で導出してみてください。
