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

応用数学微分方程式力学運動空気抵抗終端速度
2025/5/22

1. 問題の内容

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

2. 解き方の手順

(i) 運動方程式を立てる。
物体に働く力は、重力 mgmg (下向き) と空気抵抗 bv-bv (上向き) です。したがって、運動方程式は
mdvdt=mgbvm\frac{dv}{dt} = mg - bv
となります。
(ii) 与えられた解が運動方程式を満たすことを確認する。
v(t)=(v0mgb)ebmt+mgbv(t) = (v_0 - \frac{mg}{b})e^{-\frac{b}{m}t} + \frac{mg}{b} を時間 tt で微分すると、
dvdt=(v0mgb)(bm)ebmt=bm(v0mgb)ebmt\frac{dv}{dt} = (v_0 - \frac{mg}{b}) (-\frac{b}{m}) e^{-\frac{b}{m}t} = -\frac{b}{m} (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t}
これを運動方程式に代入すると、
mdvdt=m(bm(v0mgb)ebmt)=b(v0mgb)ebmtm \frac{dv}{dt} = m (-\frac{b}{m} (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t}) = -b (v_0 - \frac{mg}{b}) e^{-\frac{b}{m}t}
一方、
mgbv=mgb((v0mgb)ebmt+mgb)=mgb(v0mgb)ebmtmg=b(v0mgb)ebmtmg - bv = mg - b ((v_0 - \frac{mg}{b})e^{-\frac{b}{m}t} + \frac{mg}{b}) = mg - b (v_0 - \frac{mg}{b})e^{-\frac{b}{m}t} - mg = -b (v_0 - \frac{mg}{b})e^{-\frac{b}{m}t}
したがって、mdvdt=mgbvm\frac{dv}{dt} = mg - bv が成り立ちます。
また、t=0t = 0 のとき、
v(0)=(v0mgb)e0+mgb=v0mgb+mgb=v0v(0) = (v_0 - \frac{mg}{b})e^0 + \frac{mg}{b} = v_0 - \frac{mg}{b} + \frac{mg}{b} = v_0
となり、初期条件を満たします。
(iii) 終端速度を求める。
十分時間が経過したとき、tt \to \inftyebmt0e^{-\frac{b}{m}t} \to 0 となるので、
limtv(t)=limt[(v0mgb)ebmt+mgb]=0+mgb=mgb\lim_{t \to \infty} v(t) = \lim_{t \to \infty} [(v_0 - \frac{mg}{b})e^{-\frac{b}{m}t} + \frac{mg}{b}] = 0 + \frac{mg}{b} = \frac{mg}{b}
したがって、終端速度は mgb\frac{mg}{b} です。

3. 最終的な答え

(i) 運動方程式: mdvdt=mgbvm\frac{dv}{dt} = mg - bv
(ii) 与えられた解は運動方程式と初期条件を満たす。
(iii) 終端速度: mgb\frac{mg}{b}

「応用数学」の関連問題

一様な板ABが支柱C, Dによって水平に支えられており、その上に重さ24Nのおもりが置かれている。板が支柱C, Dから受ける力の大きさをそれぞれ$N_C$ [N], $N_D$ [N]とする。おもりは...

力学モーメント力のつり合いグラフ
2025/5/22

長さ $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