質量 $m$ の小球を初速度 $v_0$ で鉛直上向きに投げ上げたときの運動を考える問題です。小球は速度に比例する空気抵抗 $f = -bv$ を受けます。 (a) 鉛直上向きに $y$ 軸を取り、小球の運動方程式を立てる。 (b) 運動方程式を解き、$t$ 秒後のボールの速度 $v(t)$ を求め、そのグラフを描く。また、投げ上げ直後の速度の近似式を求め、そのグラフを描く。 (c) 加速度の時間変化を求め、そのグラフを描く。また、投げ上げ直後の近似式を求め、そのグラフを描く。 (d) 小球が最高点に達するまでの時間を求める。 (e) 特別な初速度 $v_0 = \frac{mg}{b}$ で小球を投げ上げた場合、空気抵抗がある場合の最高点 $h$ と、空気抵抗がない場合の最高点 $h'$ を求め、その比 $h/h'$ を求め、この $v_0$ がどう特別なのかを考察する。

応用数学運動方程式微分方程式空気抵抗力学積分
2025/4/23

1. 問題の内容

質量 mm の小球を初速度 v0v_0 で鉛直上向きに投げ上げたときの運動を考える問題です。小球は速度に比例する空気抵抗 f=bvf = -bv を受けます。
(a) 鉛直上向きに yy 軸を取り、小球の運動方程式を立てる。
(b) 運動方程式を解き、tt 秒後のボールの速度 v(t)v(t) を求め、そのグラフを描く。また、投げ上げ直後の速度の近似式を求め、そのグラフを描く。
(c) 加速度の時間変化を求め、そのグラフを描く。また、投げ上げ直後の近似式を求め、そのグラフを描く。
(d) 小球が最高点に達するまでの時間を求める。
(e) 特別な初速度 v0=mgbv_0 = \frac{mg}{b} で小球を投げ上げた場合、空気抵抗がある場合の最高点 hh と、空気抵抗がない場合の最高点 hh' を求め、その比 h/hh/h' を求め、この v0v_0 がどう特別なのかを考察する。

2. 解き方の手順

(a) 運動方程式:
ニュートンの運動方程式より、
mdvdt=mgbvm \frac{dv}{dt} = -mg - bv
(b) 運動方程式の解:
運動方程式を解くために、変数分離を行う。
dvdt=gbmv\frac{dv}{dt} = -g - \frac{b}{m}v
dvg+bmv=dt\frac{dv}{g + \frac{b}{m}v} = -dt
両辺を積分する。
dvg+bmv=dt\int \frac{dv}{g + \frac{b}{m}v} = \int -dt
mblng+bmv=t+C\frac{m}{b} \ln |g + \frac{b}{m}v| = -t + C
初期条件 t=0t = 0v=v0v = v_0 を代入して積分定数 CC を求める。
mblng+bmv0=C\frac{m}{b} \ln |g + \frac{b}{m}v_0| = C
したがって、
mblng+bmv=t+mblng+bmv0\frac{m}{b} \ln |g + \frac{b}{m}v| = -t + \frac{m}{b} \ln |g + \frac{b}{m}v_0|
lng+bmv=bmt+lng+bmv0\ln |g + \frac{b}{m}v| = -\frac{b}{m}t + \ln |g + \frac{b}{m}v_0|
g+bmv=(g+bmv0)ebmtg + \frac{b}{m}v = (g + \frac{b}{m}v_0)e^{-\frac{b}{m}t}
bmv=(g+bmv0)ebmtg\frac{b}{m}v = (g + \frac{b}{m}v_0)e^{-\frac{b}{m}t} - g
v(t)=(v0+mgb)ebmtmgbv(t) = (v_0 + \frac{mg}{b})e^{-\frac{b}{m}t} - \frac{mg}{b}
投げ上げ直後の近似:
ebmt1bmte^{-\frac{b}{m}t} \approx 1 - \frac{b}{m}t
v(t)(v0+mgb)(1bmt)mgb=v0gtv(t) \approx (v_0 + \frac{mg}{b})(1 - \frac{b}{m}t) - \frac{mg}{b} = v_0 - gt
グラフは省略。
(c) 加速度の時間変化:
a(t)=dv(t)dt=bm(v0+mgb)ebmta(t) = \frac{dv(t)}{dt} = -\frac{b}{m}(v_0 + \frac{mg}{b})e^{-\frac{b}{m}t}
a(t)=(g+bmv0)ebmta(t) = - (g + \frac{b}{m}v_0)e^{-\frac{b}{m}t}
投げ上げ直後の近似:
a(t)(g+bmv0)(1bmt)gbmv0a(t) \approx - (g + \frac{b}{m}v_0)(1 - \frac{b}{m}t) \approx -g - \frac{b}{m}v_0
グラフは省略。
(d) 最高点に達するまでの時間:
最高点では v(t)=0v(t) = 0 なので、
0=(v0+mgb)ebmtmgb0 = (v_0 + \frac{mg}{b})e^{-\frac{b}{m}t} - \frac{mg}{b}
ebmt=mg/bv0+mg/be^{-\frac{b}{m}t} = \frac{mg/b}{v_0 + mg/b}
bmt=lnmg/bv0+mg/b-\frac{b}{m}t = \ln \frac{mg/b}{v_0 + mg/b}
t=mblnmg/bv0+mg/b=mblnv0+mg/bmg/bt = -\frac{m}{b} \ln \frac{mg/b}{v_0 + mg/b} = \frac{m}{b} \ln \frac{v_0 + mg/b}{mg/b}
(e) v0=mgbv_0 = \frac{mg}{b} の場合:
空気抵抗がある場合:
v(t)=(mgb+mgb)ebmtmgb=2mgbebmtmgbv(t) = (\frac{mg}{b} + \frac{mg}{b})e^{-\frac{b}{m}t} - \frac{mg}{b} = \frac{2mg}{b}e^{-\frac{b}{m}t} - \frac{mg}{b}
v=dydtv = \frac{dy}{dt} より、
y(t)=v(t)dt=(2mgbebmtmgb)dt=2m2gb2ebmtmgbt+Cy(t) = \int v(t) dt = \int (\frac{2mg}{b}e^{-\frac{b}{m}t} - \frac{mg}{b}) dt = -\frac{2m^2g}{b^2}e^{-\frac{b}{m}t} - \frac{mg}{b}t + C
t=0t=0y=0y=0 より、
0=2m2gb2+C0 = -\frac{2m^2g}{b^2} + C
C=2m2gb2C = \frac{2m^2g}{b^2}
y(t)=2m2gb2ebmtmgbt+2m2gb2y(t) = -\frac{2m^2g}{b^2}e^{-\frac{b}{m}t} - \frac{mg}{b}t + \frac{2m^2g}{b^2}
最高点では v(t)=0v(t) = 0 より、(d)の結果より、
t=mblnmg/b+mg/bmg/b=mbln2t = \frac{m}{b} \ln \frac{mg/b + mg/b}{mg/b} = \frac{m}{b} \ln 2
h=y(mbln2)=2m2gb2ebmmbln2mgbmbln2+2m2gb2=2m2gb212m2gb2ln2+2m2gb2=m2gb2(1ln2)h = y(\frac{m}{b} \ln 2) = -\frac{2m^2g}{b^2}e^{-\frac{b}{m}\frac{m}{b} \ln 2} - \frac{mg}{b}\frac{m}{b} \ln 2 + \frac{2m^2g}{b^2} = -\frac{2m^2g}{b^2} \frac{1}{2} - \frac{m^2g}{b^2} \ln 2 + \frac{2m^2g}{b^2} = \frac{m^2g}{b^2}(1 - \ln 2)
空気抵抗がない場合:
v(t)=v0gt=mgbgtv(t) = v_0 - gt = \frac{mg}{b} - gt
v=0v = 0 のとき、 t=mgbgt = \frac{mg}{bg}
t=mbt = \frac{m}{b}
y(t)=v(t)dt=mgbt12gt2+Cy(t) = \int v(t) dt = \frac{mg}{b}t - \frac{1}{2}gt^2 + C
y(0)=0y(0) = 0 より C=0C = 0
h=y(mb)=mgbmb12g(mb)2=m2gb212m2gb2=12m2gb2h' = y(\frac{m}{b}) = \frac{mg}{b} \frac{m}{b} - \frac{1}{2}g (\frac{m}{b})^2 = \frac{m^2g}{b^2} - \frac{1}{2} \frac{m^2g}{b^2} = \frac{1}{2} \frac{m^2g}{b^2}
hh=m2gb2(1ln2)12m2gb2=2(1ln2)\frac{h}{h'} = \frac{\frac{m^2g}{b^2}(1 - \ln 2)}{\frac{1}{2} \frac{m^2g}{b^2}} = 2(1 - \ln 2)
v0=mgbv_0 = \frac{mg}{b} は、終端速度と等しい。

3. 最終的な答え

(a) 運動方程式: mdvdt=mgbvm \frac{dv}{dt} = -mg - bv
(b) v(t)=(v0+mgb)ebmtmgbv(t) = (v_0 + \frac{mg}{b})e^{-\frac{b}{m}t} - \frac{mg}{b}, v(t)v0gtv(t) \approx v_0 - gt
(c) a(t)=(g+bmv0)ebmta(t) = - (g + \frac{b}{m}v_0)e^{-\frac{b}{m}t}, a(t)gbmv0a(t) \approx -g - \frac{b}{m}v_0
(d) t=mblnv0+mg/bmg/bt = \frac{m}{b} \ln \frac{v_0 + mg/b}{mg/b}
(e) h=m2gb2(1ln2)h = \frac{m^2g}{b^2}(1 - \ln 2), h=12m2gb2h' = \frac{1}{2} \frac{m^2g}{b^2}, hh=2(1ln2)\frac{h}{h'} = 2(1 - \ln 2)

「応用数学」の関連問題

3点$(1, 1)$, $(2, 5)$, $(3, 3)$を近似する回帰直線の式を、最小二乗法を用いて求める。

最小二乗法回帰直線線形代数統計
2025/6/2

アルゴンガスを初期体積 $2000 \ cm^3$ から最終体積 $500 \ cm^3$ まで圧縮し、同時に初期温度 $300 \ K$ から最終温度 $400 \ K$ まで加熱した時の、モルエン...

熱力学エントロピー積分対数
2025/6/2

グラフが与えられており、横軸は$V_0$、縦軸は$\theta_{max}$、パラメータは$R$です。$R = 0$の場合のグラフの概形と、ある$V_0$に対する$\theta_{max}$の値がどの...

グラフ解釈物理モデル
2025/6/2

DNAの10塩基対の長さが3.4 nmであるとき、長さが102 μmの2本鎖DNA分子に含まれるヌクレオチドの数を求める。

単位変換比例計算生物学
2025/6/2

ある工場で製品Aと製品Bを製造する。製品Aを1トン作るには原料Pが2トン、原料Qが4トン必要。製品Bを1トン作るには原料Pが6トン、原料Qが2トン必要。1ヶ月あたり、原料Pは140トン、原料Qは120...

線形計画法最適化不等式グラフ
2025/6/1

H国とF国の2国を考え、H国のみに市場が存在する。H国の需要曲線は、$D_H = -2p + 20$ である。最初に、H国の価格が9、F国の価格が8の場合を考える。次に、H国の価格が2、F国の価格が1...

経済学需要曲線価格弾力性収入利益
2025/6/1

ある工場で製品Aと製品Bを生産します。製品A, Bを1トン生産するのに必要な原料P, Qの量と製品A, Bの価格が表で与えられています。この工場へ1日に供給できる原料Pが最大9トン、原料Qが最大8トン...

線形計画法最適化制約条件グラフ
2025/6/1

質量 $m$ [kg] の物体が軽い糸1で天井から吊るされ、水平方向に糸2で引かれている。糸1が天井となす角は $\theta$ [°] である。重力加速度の大きさを $g$ [m/s²] とする。 ...

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

質量3.0kgの物体が糸1と糸2によって引っ張られている。糸1は鉛直方向から30°の角度をなし、糸2は水平方向に引っ張られている。重力加速度の大きさは9.8m/s²、$\sqrt{3} = 1.7$と...

力学ベクトル力の釣り合い三角関数物理
2025/6/1

ボールを高さ20の位置から45度または30度の角度で発射した時に、ボールが地面に落下するまでの水平距離を計算し、45度と30度どちらの角度で発射した方が水平距離が長いかを判断する。

力学放物運動水平距離角度物理
2025/6/1