質量 $m$ の物体が、初期位置 $(x,y) = (0, h)$ から水平方向に初速度 $v_0$ で射出される。重力加速度を $g$ (鉛直上向きが正) とし、以下の問いに答える。 (1) 空気抵抗がない場合、任意の時刻 $t$ における $x$ 方向と $y$ 方向の速度を求める。運動方程式を微分方程式の形で書き、解くこと。 (2) (1) の場合、任意の時刻 $t$ における物体の位置を求める。 (3) 粘性抵抗 (比例係数 $\gamma$) が働く場合、任意の時刻 $t$ における $x$ 方向と $y$ 方向の速度を求める。 (4) (3) の場合、任意の時刻 $t$ における物体の位置を求める。 (5) 慣性抵抗 (比例係数 $\gamma$) が働く場合、任意の時刻 $t$ における $x$ 方向と $y$ 方向の速度を求める。

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

1. 問題の内容

質量 mm の物体が、初期位置 (x,y)=(0,h)(x,y) = (0, h) から水平方向に初速度 v0v_0 で射出される。重力加速度を gg (鉛直上向きが正) とし、以下の問いに答える。
(1) 空気抵抗がない場合、任意の時刻 tt における xx 方向と yy 方向の速度を求める。運動方程式を微分方程式の形で書き、解くこと。
(2) (1) の場合、任意の時刻 tt における物体の位置を求める。
(3) 粘性抵抗 (比例係数 γ\gamma) が働く場合、任意の時刻 tt における xx 方向と yy 方向の速度を求める。
(4) (3) の場合、任意の時刻 tt における物体の位置を求める。
(5) 慣性抵抗 (比例係数 γ\gamma) が働く場合、任意の時刻 tt における xx 方向と yy 方向の速度を求める。

2. 解き方の手順

(1) 空気抵抗がない場合:
* xx 方向の運動方程式: mdvxdt=0m \frac{dv_x}{dt} = 0。初期条件は vx(0)=v0v_x(0) = v_0
解は vx(t)=v0v_x(t) = v_0
* yy 方向の運動方程式: mdvydt=mgm \frac{dv_y}{dt} = -mg。初期条件は vy(0)=0v_y(0) = 0
解は vy(t)=gtv_y(t) = -gt
(2) 空気抵抗がない場合の位置:
* xx 方向の速度 vx(t)=dxdt=v0v_x(t) = \frac{dx}{dt} = v_0 を時間 tt で積分する。初期条件は x(0)=0x(0) = 0
x(t)=v0tx(t) = v_0 t
* yy 方向の速度 vy(t)=dydt=gtv_y(t) = \frac{dy}{dt} = -gt を時間 tt で積分する。初期条件は y(0)=hy(0) = h
y(t)=h12gt2y(t) = h - \frac{1}{2}gt^2
(3) 粘性抵抗がある場合:
* xx 方向の運動方程式: mdvxdt=γvxm \frac{dv_x}{dt} = -\gamma v_x。初期条件は vx(0)=v0v_x(0) = v_0
dvxvx=γmdt\frac{dv_x}{v_x} = -\frac{\gamma}{m} dt より積分して lnvx=γmt+C\ln v_x = -\frac{\gamma}{m} t + C。初期条件より C=lnv0C = \ln v_0
vx(t)=v0eγmtv_x(t) = v_0 e^{-\frac{\gamma}{m}t}
* yy 方向の運動方程式: mdvydt=mgγvym \frac{dv_y}{dt} = -mg - \gamma v_y。初期条件は vy(0)=0v_y(0) = 0
dvydt=gγmvy\frac{dv_y}{dt} = -g - \frac{\gamma}{m} v_y。これを解くと、
vy(t)=mgγ(1eγmt)v_y(t) = -\frac{mg}{\gamma} (1 - e^{-\frac{\gamma}{m}t})
(4) 粘性抵抗がある場合の位置:
* xx 方向の位置 x(t)x(t)vx(t)=dxdt=v0eγmtv_x(t) = \frac{dx}{dt} = v_0 e^{-\frac{\gamma}{m}t} を時間 tt で積分する。初期条件は x(0)=0x(0) = 0
x(t)=0tv0eγmtdt=v0[mγeγmt]0t=mv0γ(1eγmt)x(t) = \int_0^t v_0 e^{-\frac{\gamma}{m}t'} dt' = v_0 \left[ -\frac{m}{\gamma} e^{-\frac{\gamma}{m}t'} \right]_0^t = \frac{m v_0}{\gamma} (1 - e^{-\frac{\gamma}{m}t})
* yy 方向の位置 y(t)y(t)vy(t)=dydt=mgγ(1eγmt)v_y(t) = \frac{dy}{dt} = -\frac{mg}{\gamma} (1 - e^{-\frac{\gamma}{m}t}) を時間 tt で積分する。初期条件は y(0)=hy(0) = h
y(t)=0tmgγ(1eγmt)dt+h=mgγ[t+mγeγmt]0t+h=mgγ(t+mγeγmtmγ)+h=hmgtγm2gγ2(eγmt1)y(t) = \int_0^t -\frac{mg}{\gamma} (1 - e^{-\frac{\gamma}{m}t'}) dt' + h = -\frac{mg}{\gamma} \left[ t' + \frac{m}{\gamma} e^{-\frac{\gamma}{m}t'} \right]_0^t + h = -\frac{mg}{\gamma} (t + \frac{m}{\gamma} e^{-\frac{\gamma}{m}t} - \frac{m}{\gamma}) + h = h - \frac{mgt}{\gamma} - \frac{m^2g}{\gamma^2} (e^{-\frac{\gamma}{m}t} - 1)
(5) 慣性抵抗がある場合:
* xx 方向の運動方程式: mdvxdt=γvx2m \frac{dv_x}{dt} = -\gamma v_x^2。初期条件は vx(0)=v0v_x(0) = v_0
dvxvx2=γmdt\frac{dv_x}{v_x^2} = -\frac{\gamma}{m} dt より積分して 1vx=γmt+C-\frac{1}{v_x} = -\frac{\gamma}{m}t + C。初期条件より C=1v0C = -\frac{1}{v_0}
1vx=γmt+1v0=γv0t+mmv0\frac{1}{v_x} = \frac{\gamma}{m}t + \frac{1}{v_0} = \frac{\gamma v_0 t + m}{mv_0}
vx(t)=mv0γv0t+m=v01+γv0mtv_x(t) = \frac{mv_0}{\gamma v_0 t + m} = \frac{v_0}{1 + \frac{\gamma v_0}{m} t}
* yy 方向の運動方程式: mdvydt=mgγvy2m \frac{dv_y}{dt} = -mg - \gamma v_y^2。初期条件は vy(0)=0v_y(0) = 0
dvydt=gγmvy2\frac{dv_y}{dt} = -g - \frac{\gamma}{m} v_y^2dvyg+γmvy2=dt\int \frac{dv_y}{g + \frac{\gamma}{m} v_y^2} = -\int dt を解く。
vy(t)=mgγtanh(γgmt)v_y(t) = -\sqrt{\frac{mg}{\gamma}} \tanh \left( \sqrt{\frac{\gamma g}{m}}t \right)

3. 最終的な答え

(1) 空気抵抗がない場合:
vx(t)=v0v_x(t) = v_0
vy(t)=gtv_y(t) = -gt
(2) 空気抵抗がない場合の位置:
x(t)=v0tx(t) = v_0 t
y(t)=h12gt2y(t) = h - \frac{1}{2}gt^2
(3) 粘性抵抗がある場合:
vx(t)=v0eγmtv_x(t) = v_0 e^{-\frac{\gamma}{m}t}
vy(t)=mgγ(1eγmt)v_y(t) = -\frac{mg}{\gamma} (1 - e^{-\frac{\gamma}{m}t})
(4) 粘性抵抗がある場合の位置:
x(t)=mv0γ(1eγmt)x(t) = \frac{m v_0}{\gamma} (1 - e^{-\frac{\gamma}{m}t})
y(t)=hmgtγm2gγ2(eγmt1)y(t) = h - \frac{mgt}{\gamma} - \frac{m^2g}{\gamma^2} (e^{-\frac{\gamma}{m}t} - 1)
(5) 慣性抵抗がある場合:
vx(t)=v01+γv0mtv_x(t) = \frac{v_0}{1 + \frac{\gamma v_0}{m} t}
vy(t)=mgγtanh(γgmt)v_y(t) = -\sqrt{\frac{mg}{\gamma}} \tanh \left( \sqrt{\frac{\gamma g}{m}}t \right)

「応用数学」の関連問題

一定の速度で走る列車が、長さ250mのトンネルを通過するのに20秒、長さ550mの鉄橋を通過するのに35秒かかる。列車の長さを求める。

速さ距離方程式線形方程式
2025/6/20

電車が時速72kmでトンネルを通過する。トンネルの長さは800m、電車の長さは200mである。電車がトンネルを通過するのに何秒かかるかを求める。

速度距離時間物理
2025/6/20

記録タイマーで等加速度直線運動を測定した結果から、以下の3つの問いに答える。 (1) OA間の平均の速さを求める。 (2) 打点Oの時刻を0とし、横軸に時間$t$[s]、縦軸に速さ$v$[m/s]をと...

等加速度直線運動平均速度v-tグラフ加速度
2025/6/20

片持ち梁の中央 (B点) と先端 (C点) にそれぞれ $2P$ と $P$ の集中荷重がかかっている。 1. 壁からの反力 $R_A$ と固定モーメント $M_A$ を求める。

構造力学力学モーメントせん断力曲げモーメント片持ち梁
2025/6/20

長さ $l$ (m) の片持ち梁に、単位長さあたり $q$ (N) の等分布荷重が作用している。壁に生じる支持反力と固定モーメントを求める。

構造力学片持ち梁等分布荷重支持反力固定モーメント力学
2025/6/20

単純支持梁に、A点から距離$l$の位置に荷重$P$、距離$2l$の位置に荷重$2P$が作用している。A点とB点における反力$R_A$、$R_B$を求めよ。

力学構造力学モーメント力の釣り合い
2025/6/20

Aさんが家から2520m離れた駅まで自転車で向かいました。出発して9分後に自転車が故障し、そこから歩いて駅まで行ったところ、予定より7分遅れて到着しました。自転車の速さが歩く速さの2.4倍であるとき、...

速さ距離方程式文章問題
2025/6/20

ある財の需要曲線が $X = -3P + 2400$ で、供給曲線が $X = 5P$ である。ここで、$X$ は数量、$P$ は価格である。政府により上限価格が 150 に規制されたとき、取引量、消...

経済学需要と供給価格規制消費者余剰生産者余剰死荷重
2025/6/20

単純支持はりに集中荷重 $P = 200N$ が加えられている。せん断力図 (SFD) と曲げモーメント図 (BMD) を描く問題です。支点Aから荷重までの距離は0.3m、荷重から支点Bまでの距離は0...

力学構造力学せん断力図曲げモーメント図静力学
2025/6/20

ある財の市場において、需要曲線が $X = -P + 200$ で表される。当初、参入規制があり、供給曲線は $X = (1/2)P - 10$ であった。その後、参入規制が撤廃され、供給曲線は $X...

経済学需要曲線供給曲線均衡価格総余剰消費者余剰生産者余剰
2025/6/20