質量 $m$ の物体が空気抵抗を受けながら落下する運動について考えます。鉛直上向きを $y$ 軸とします。物体には下向きの重力 $mg$ と、速度に比例する上向きの空気抵抗 $bv$ が働きます。以下の問いに答えます。 (i) 物体が満たす運動方程式を立ててください。 (ii) $t=0$ で $v=v_0$ を満たす運動方程式の解が $v(t) = \left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t} + \frac{mg}{b}$ となることを確かめてください。 (iii) 十分時間が経過したとき、速度が一定の速度に漸近することを示し、終端速度を求めてください。

応用数学微分方程式運動方程式力学終端速度指数関数
2025/5/22

1. 問題の内容

質量 mm の物体が空気抵抗を受けながら落下する運動について考えます。鉛直上向きを yy 軸とします。物体には下向きの重力 mgmg と、速度に比例する上向きの空気抵抗 bvbv が働きます。以下の問いに答えます。
(i) 物体が満たす運動方程式を立ててください。
(ii) t=0t=0v=v0v=v_0 を満たす運動方程式の解が v(t)=(v0mgb)ebmt+mgbv(t) = \left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t} + \frac{mg}{b} となることを確かめてください。
(iii) 十分時間が経過したとき、速度が一定の速度に漸近することを示し、終端速度を求めてください。

2. 解き方の手順

(i) 運動方程式の導出
物体の運動方程式は、ニュートンの運動の第二法則 F=maF = ma から導かれます。
鉛直下向きを正とすると、物体に働く力は重力 mgmg と空気抵抗 bv-bv です。したがって、運動方程式は次のようになります。
ma=mgbvma = mg - bv
加速度 aa は速度 vv の時間微分 dvdt\frac{dv}{dt} なので、運動方程式は以下のようになります。
mdvdt=mgbvm\frac{dv}{dt} = mg - bv
(ii) 解の確認
与えられた解 v(t)=(v0mgb)ebmt+mgbv(t) = \left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t} + \frac{mg}{b} が運動方程式を満たすことを確認します。まず、v(t)v(t) を時間 tt で微分します。
dvdt=(v0mgb)(bm)ebmt=bm(v0mgb)ebmt\frac{dv}{dt} = \left(v_0 - \frac{mg}{b}\right) \left(-\frac{b}{m}\right) e^{-\frac{b}{m}t} = -\frac{b}{m}\left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t}
これを運動方程式 mdvdt=mgbvm\frac{dv}{dt} = mg - bv に代入します。
m(bm(v0mgb)ebmt)=mgb((v0mgb)ebmt+mgb)m\left(-\frac{b}{m}\left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t}\right) = mg - b\left(\left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t} + \frac{mg}{b}\right)
b(v0mgb)ebmt=mgb(v0mgb)ebmtmg-b\left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t} = mg - b\left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t} - mg
b(v0mgb)ebmt=b(v0mgb)ebmt-b\left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t} = -b\left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t}
左辺と右辺が一致するので、与えられた v(t)v(t) は運動方程式を満たします。
次に初期条件 t=0t=0v=v0v=v_0 を満たすことを確認します。
v(0)=(v0mgb)ebm(0)+mgb=(v0mgb)(1)+mgb=v0mgb+mgb=v0v(0) = \left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}(0)} + \frac{mg}{b} = \left(v_0 - \frac{mg}{b}\right)(1) + \frac{mg}{b} = v_0 - \frac{mg}{b} + \frac{mg}{b} = v_0
初期条件も満たすので、与えられた v(t)v(t) は運動方程式の解です。
(iii) 終端速度の導出
十分時間が経過したとき、tt \to \inftyv(t)v(t) がどうなるかを調べます。
limtv(t)=limt[(v0mgb)ebmt+mgb]\lim_{t \to \infty} v(t) = \lim_{t \to \infty} \left[\left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t} + \frac{mg}{b}\right]
tt \to \infty のとき、ebmt0e^{-\frac{b}{m}t} \to 0 なので、
limtv(t)=(v0mgb)(0)+mgb=mgb\lim_{t \to \infty} v(t) = \left(v_0 - \frac{mg}{b}\right)(0) + \frac{mg}{b} = \frac{mg}{b}
したがって、十分時間が経過すると速度は mgb\frac{mg}{b} に漸近します。これが終端速度です。

3. 最終的な答え

(i) 運動方程式: mdvdt=mgbvm\frac{dv}{dt} = mg - bv
(ii) 解の確認: 与えられた v(t)=(v0mgb)ebmt+mgbv(t) = \left(v_0 - \frac{mg}{b}\right)e^{-\frac{b}{m}t} + \frac{mg}{b} は運動方程式と初期条件を満たします。
(iii) 終端速度: mgb\frac{mg}{b}

「応用数学」の関連問題

長さ $L$、質量 $m$ の一様な棒ABの一端Aが壁に蝶番で固定され、他端Bが糸で壁に繋がれています。糸と水平面のなす角は30度です。Aにはたらく力の水平方向と鉛直方向の成分をそれぞれ $F_x$ ...

力学モーメント力のつりあい剛体物理
2025/5/22

半径2の円周上を運動する質点A, Bについて、時刻tにおける位置ベクトル$r^A(t)$と$r^B(t)$が与えられています。 $r^A(t) = 2 (cos(\frac{\pi t}{3} - \...

ベクトル円運動角速度加速度軌跡微分
2025/5/22

半径2の円周上を運動する質点A, Bについて、時刻tにおける位置がそれぞれ以下のように与えられています。 $r^A(t) = 2(\cos(\frac{\pi t}{3} - \frac{\pi}{6...

ベクトル円運動角速度加速度運動解析微分
2025/5/22

2階線形同次微分方程式 $4y'' - 12y' + 9y = 0$ の一般解を求め、初期条件 $x = 0$ のとき $y = 1$, $y' = 2$ を満たす解を、選択肢の中から選びます。

微分方程式線形微分方程式初期条件一般解特性方程式
2025/5/22

半径2の円周上を運動する質点AとBについて、時刻 $t$ における位置ベクトルがそれぞれ $r^A(t) = 2(\cos(\frac{\pi t}{3} - \frac{\pi}{6})i + \s...

ベクトル円運動角速度加速度接線成分法線成分軌跡
2025/5/22

半径2の円周上を運動する質点AとBについて、時刻 $t$ におけるそれぞれの位置が与えられています。 Aの位置ベクトル: $r^A(t) = 2(\cos(\frac{\pi t}{3} - \fra...

ベクトル円運動角速度加速度軌跡微分
2025/5/22

10 kgとP [kg]の物体がロープで吊るされ、図のような位置で静止している。ロープAB, BC, CDの張力の大きさとPの値を求めよ。

力学力の釣り合いベクトル三角関数
2025/5/22

デルタ結線された三相交流回路に、線間電圧 $E_{ab} = 200\angle 0^\circ$, $E_{bc} = 200\angle -120^\circ$, $E_{ca} = 200\an...

電気回路三相交流インピーダンス複素数オームの法則
2025/5/22

質量5.0kgの小球が軽い糸で天井から吊るされている。小球は水平方向に力$F$で押され、糸が天井と30°の角をなして静止している。このときの糸の張力$T$を求める。重力加速度は9.8m/s²とする。

力学ベクトル力のつりあい三角関数
2025/5/22

水平な地面から小球を斜方投射した。重力加速度の大きさを $9.8 \ m/s^2$ とする。初速度の水平成分の大きさは $4.9 \ m/s$、鉛直成分の大きさは $9.8 \ m/s$ である。 (...

力学物理斜方投射運動重力
2025/5/22