(3) 粘性抵抗が働く場合、物体の任意の時刻 $t$ における $x$ 方向と $y$ 方向の速度を求めよ。ただし、粘性抵抗の比例係数を $\gamma_1$ とする。 (4) (3)の場合に、物体の任意の時刻 $t$ における位置を求めよ。 (5) 慣性抵抗が働く場合、物体の任意の時刻 $t$ における $x$ 方向と $y$ 方向の速度を求めよ。ただし、慣性抵抗の比例係数を $\gamma_2$ とする。

応用数学微分方程式力学粘性抵抗運動
2025/6/23

1. 問題の内容

(3) 粘性抵抗が働く場合、物体の任意の時刻 tt における xx 方向と yy 方向の速度を求めよ。ただし、粘性抵抗の比例係数を γ1\gamma_1 とする。
(4) (3)の場合に、物体の任意の時刻 tt における位置を求めよ。
(5) 慣性抵抗が働く場合、物体の任意の時刻 tt における xx 方向と yy 方向の速度を求めよ。ただし、慣性抵抗の比例係数を γ2\gamma_2 とする。

2. 解き方の手順

(3) 粘性抵抗が働く場合:
xx方向の運動方程式は、mdvxdt=γ1vxm \frac{dv_x}{dt} = -\gamma_1 v_x である。初期条件を vx(0)=v0v_x(0) = v_0 とすると、
dvxvx=γ1mdt\frac{dv_x}{v_x} = -\frac{\gamma_1}{m} dt
両辺を積分して、
v0vxdvxvx=0tγ1mdt\int_{v_0}^{v_x} \frac{dv_x'}{v_x'} = \int_0^t -\frac{\gamma_1}{m} dt'
lnvxv0=γ1mt\ln \frac{v_x}{v_0} = -\frac{\gamma_1}{m} t
vx(t)=v0eγ1mtv_x(t) = v_0 e^{-\frac{\gamma_1}{m} t}
yy方向の運動方程式は、mdvydt=mgγ1vym \frac{dv_y}{dt} = -mg - \gamma_1 v_y である。
dvydt=gγ1mvy\frac{dv_y}{dt} = -g - \frac{\gamma_1}{m} v_y
dvyg+γ1mvy=dt\frac{dv_y}{g + \frac{\gamma_1}{m} v_y} = -dt
0vydvyg+γ1mvy=0tdt\int_{0}^{v_y} \frac{dv_y'}{g + \frac{\gamma_1}{m} v_y'} = \int_0^t -dt'
mγ1lng+γ1mvyg=t\frac{m}{\gamma_1} \ln \frac{g + \frac{\gamma_1}{m} v_y}{g} = -t
lng+γ1mvyg=γ1mt\ln \frac{g + \frac{\gamma_1}{m} v_y}{g} = -\frac{\gamma_1}{m} t
g+γ1mvyg=eγ1mt\frac{g + \frac{\gamma_1}{m} v_y}{g} = e^{-\frac{\gamma_1}{m} t}
g+γ1mvy=geγ1mtg + \frac{\gamma_1}{m} v_y = g e^{-\frac{\gamma_1}{m} t}
γ1mvy=g(eγ1mt1)\frac{\gamma_1}{m} v_y = g(e^{-\frac{\gamma_1}{m} t} - 1)
vy(t)=mgγ1(eγ1mt1)v_y(t) = \frac{mg}{\gamma_1} (e^{-\frac{\gamma_1}{m} t} - 1)
(4) (3)の場合:
xx方向の位置 x(t)x(t) を求めるには、vx(t)v_x(t) を積分する。初期位置を x(0)=0x(0) = 0 とする。
x(t)=0tvx(t)dt=0tv0eγ1mtdt=v0[mγ1eγ1mt]0t=mv0γ1(1eγ1mt)x(t) = \int_0^t v_x(t') dt' = \int_0^t v_0 e^{-\frac{\gamma_1}{m} t'} dt' = v_0 \left[ -\frac{m}{\gamma_1} e^{-\frac{\gamma_1}{m} t'} \right]_0^t = \frac{m v_0}{\gamma_1} (1 - e^{-\frac{\gamma_1}{m} t})
yy方向の位置 y(t)y(t) を求めるには、vy(t)v_y(t) を積分する。初期位置を y(0)=hy(0) = h とする。
y(t)=0tvy(t)dt+h=0tmgγ1(eγ1mt1)dt+h=mgγ1[mγ1eγ1mtt]0t+h=mgγ1(mγ1eγ1mtt+mγ1)+h=hmgtγ1+m2gγ12(1eγ1mt)y(t) = \int_0^t v_y(t') dt' + h = \int_0^t \frac{mg}{\gamma_1} (e^{-\frac{\gamma_1}{m} t'} - 1) dt' + h = \frac{mg}{\gamma_1} \left[ -\frac{m}{\gamma_1} e^{-\frac{\gamma_1}{m} t'} - t' \right]_0^t + h = \frac{mg}{\gamma_1} \left( -\frac{m}{\gamma_1} e^{-\frac{\gamma_1}{m} t} - t + \frac{m}{\gamma_1} \right) + h = h - \frac{mgt}{\gamma_1} + \frac{m^2 g}{\gamma_1^2} (1 - e^{-\frac{\gamma_1}{m} t})
(5) 慣性抵抗が働く場合:
xx方向の運動方程式は、mdvxdt=γ2(vx)2m \frac{dv_x}{dt} = -\gamma_2 (v_x)^2 である。
dvxvx2=γ2mdt\frac{dv_x}{v_x^2} = -\frac{\gamma_2}{m} dt
v0vxdvx(vx)2=0tγ2mdt\int_{v_0}^{v_x} \frac{dv_x'}{(v_x')^2} = \int_0^t -\frac{\gamma_2}{m} dt'
[1vx]v0vx=γ2mt\left[ -\frac{1}{v_x'} \right]_{v_0}^{v_x} = -\frac{\gamma_2}{m} t
1vx+1v0=γ2mt-\frac{1}{v_x} + \frac{1}{v_0} = -\frac{\gamma_2}{m} t
1vx=1v0+γ2mt=m+γ2v0tmv0\frac{1}{v_x} = \frac{1}{v_0} + \frac{\gamma_2}{m} t = \frac{m + \gamma_2 v_0 t}{m v_0}
vx(t)=mv0m+γ2v0tv_x(t) = \frac{m v_0}{m + \gamma_2 v_0 t}
yy方向の運動方程式は、mdvydt=mgγ2(vy)2m \frac{dv_y}{dt} = -mg - \gamma_2 (v_y)^2 である。
dvydt=gγ2mvy2\frac{dv_y}{dt} = -g - \frac{\gamma_2}{m} v_y^2
dvyg+γ2mvy2=dt\frac{dv_y}{g + \frac{\gamma_2}{m} v_y^2} = -dt
初期条件 vy(0)=0v_y(0) = 0 とする。
0vydvyg+γ2m(vy)2=0tdt\int_{0}^{v_y} \frac{dv_y'}{g + \frac{\gamma_2}{m} (v_y')^2} = - \int_0^t dt'
0vydvyg+γ2m(vy)2=mgγ2arctan(vyγ2mg)=t\int_{0}^{v_y} \frac{dv_y'}{g + \frac{\gamma_2}{m} (v_y')^2} = \sqrt{\frac{m}{g\gamma_2}} \arctan(v_y\sqrt{\frac{\gamma_2}{mg}}) = -t
arctan(vyγ2mg)=tgγ2m\arctan(v_y\sqrt{\frac{\gamma_2}{mg}}) = -t\sqrt{\frac{g\gamma_2}{m}}
vyγ2mg=tan(tgγ2m)v_y\sqrt{\frac{\gamma_2}{mg}} = \tan(-t\sqrt{\frac{g\gamma_2}{m}})
vy=mgγ2tan(tgγ2m)v_y = -\sqrt{\frac{mg}{\gamma_2}} \tan(t\sqrt{\frac{g\gamma_2}{m}})

3. 最終的な答え

(3) 粘性抵抗が働く場合:
vx(t)=v0eγ1mtv_x(t) = v_0 e^{-\frac{\gamma_1}{m} t}
vy(t)=mgγ1(eγ1mt1)v_y(t) = \frac{mg}{\gamma_1} (e^{-\frac{\gamma_1}{m} t} - 1)
(4) (3)の場合:
x(t)=mv0γ1(1eγ1mt)x(t) = \frac{m v_0}{\gamma_1} (1 - e^{-\frac{\gamma_1}{m} t})
y(t)=hmgtγ1+m2gγ12(1eγ1mt)y(t) = h - \frac{mgt}{\gamma_1} + \frac{m^2 g}{\gamma_1^2} (1 - e^{-\frac{\gamma_1}{m} t})
(5) 慣性抵抗が働く場合:
vx(t)=mv0m+γ2v0tv_x(t) = \frac{m v_0}{m + \gamma_2 v_0 t}
vy(t)=mgγ2tan(tgγ2m)v_y(t) = -\sqrt{\frac{mg}{\gamma_2}} \tan(t\sqrt{\frac{g\gamma_2}{m}})

「応用数学」の関連問題

小球が水平面に対して30°の角度で速度 $v$ で衝突し、跳ね返った直後の速度のx成分 $v'_x$ とy成分 $v'_y$ を、反発係数 $\frac{1}{3}$ を用いて、$v$ で表す問題です...

力学衝突ベクトル物理
2025/6/23

与えられた7つの2階常微分方程式の一般解を求める問題です。特に、(5)から(7)については、括弧内に与えられた斉次解を利用することが推奨されています。

常微分方程式微分方程式解法一般解
2025/6/23

与えられた問題はDSB(両側波帯)変調に関する3つの設問です。 1) DSB変調波 $f_{dsb}(t)$ を、搬送波 $V_c \sin(\omega_c t)$、情報信号 $s(t)$ および変...

信号処理フーリエ変換DSB変調周波数スペクトル
2025/6/23

質量0.50kgの球を軽い糸でつるし、糸の上端を持って鉛直上向きに球を引き上げる問題です。 (1) 球にはたらく重力の大きさを求めます。 (2) 球が一定の速さで上昇する場合と、鉛直下向きの加速度で上...

力学運動方程式重力張力加速度
2025/6/23

質量 $0.50 \mathrm{kg}$ の球に働く重力の大きさを求める問題です。重力加速度を $g = 9.8 \mathrm{m/s^2}$ とします。

物理力学重力質量重力加速度
2025/6/23

質量5.0kgの物体に、大きさ20Nの力が角度$\theta$で加えられています。ここで、$\tan \theta = \frac{3}{4}$です。物体の加速度の向きと大きさを求める問題です。

力学運動方程式ベクトル三角関数加速度
2025/6/23

質量 5.0 kg の物体が、傾斜角 $\theta$ ($\sin \theta = \frac{3}{5}$, $\cos \theta = \frac{4}{5}$) の斜面上に置かれている。物...

力学運動方程式加速度重力張力
2025/6/23

以下の式を証明する。 $(a \times b) \cdot (a \times b) + (a \cdot b)^2 = |a|^2 |b|^2$

ベクトル内積外積ベクトル三重積正射影
2025/6/23

問題2-7は、初速度 $v_0$、角度 $\theta$ で打ち出された物体が放物線を描くこと、すなわち、物体の軌跡が $y = Ax^2 + Bx$ の形で表されることを示す問題です。

力学放物運動軌跡物理
2025/6/23

あるプロ野球球団の1年間の応援グッズの売上が、売場Aから売場Eまでまとめられています。売場Aから売場Dまでの売上高(万円)と、メガホン、タオル、Tシャツの売上個数が分かっています。売場Eの売上高を推測...

線形代数連立方程式回帰分析モデル化
2025/6/23