質量 $m$ のおもりが長さ $L$ の棒で吊るされた振り子を考える。支点は鉛直下方に等速 $V_0$ で下降しており、時刻 $t$ での支点の座標 $y_0$ は $y_0 = V_0 t$ で表される。おもりは速度に比例する抵抗力(大きさ $c\vec{v}$)を受ける。 (1) $c=0$ のとき、振り子の振れ角 $\theta$ は小さく、単振動であるとして、振り子の周期 $T_0$ を求めよ。次に、$y_0$ が等加速度 $a$ で運動し、$y_0 = \frac{1}{2}at^2 + V_0 t$ であるとして、振り子の周期 $T_a$ を求めよ。 (2) 振り子の振れ角 $\theta$ は小さく、$\sin\theta \approx \theta$, $\cos\theta \approx 1$ と近似する。運動方程式から次の式を導き、$P$, $Q$, $R$ にあてはまる文字式を求めよ。 $\frac{d^2\theta}{dt^2} = P\theta + Q\frac{d\theta}{dt} + R$ (3) $0 < c < c_0$ のとき、運動は減衰振動となる。$c_0$ を求めよ。 (4) $c = m\sqrt{\frac{g}{L}}$, $V_0 = \frac{1}{2}\sqrt{gL}$ とする。時刻 $t=0$ のとき、おもりは $(0, L)$ の位置にあり、速度は $(U_0, -V_0)$ であったとする。$\Theta(t)$ を求めよ。

応用数学力学振り子微分方程式減衰振動
2025/7/18
はい、承知いたしました。以下に解答を記載します。

1. 問題の内容

質量 mm のおもりが長さ LL の棒で吊るされた振り子を考える。支点は鉛直下方に等速 V0V_0 で下降しており、時刻 tt での支点の座標 y0y_0y0=V0ty_0 = V_0 t で表される。おもりは速度に比例する抵抗力(大きさ cvc\vec{v})を受ける。
(1) c=0c=0 のとき、振り子の振れ角 θ\theta は小さく、単振動であるとして、振り子の周期 T0T_0 を求めよ。次に、y0y_0 が等加速度 aa で運動し、y0=12at2+V0ty_0 = \frac{1}{2}at^2 + V_0 t であるとして、振り子の周期 TaT_a を求めよ。
(2) 振り子の振れ角 θ\theta は小さく、sinθθ\sin\theta \approx \theta, cosθ1\cos\theta \approx 1 と近似する。運動方程式から次の式を導き、PP, QQ, RR にあてはまる文字式を求めよ。
d2θdt2=Pθ+Qdθdt+R\frac{d^2\theta}{dt^2} = P\theta + Q\frac{d\theta}{dt} + R
(3) 0<c<c00 < c < c_0 のとき、運動は減衰振動となる。c0c_0 を求めよ。
(4) c=mgLc = m\sqrt{\frac{g}{L}}, V0=12gLV_0 = \frac{1}{2}\sqrt{gL} とする。時刻 t=0t=0 のとき、おもりは (0,L)(0, L) の位置にあり、速度は (U0,V0)(U_0, -V_0) であったとする。Θ(t)\Theta(t) を求めよ。

2. 解き方の手順

(1) c=0c=0y0=V0ty_0 = V_0 t の場合:
振り子の運動方程式は、
mLd2θdt2=mgsinθmL\frac{d^2\theta}{dt^2} = -mg\sin\theta
θ\theta が小さいとき sinθθ\sin\theta \approx \theta なので、
d2θdt2=gLθ\frac{d^2\theta}{dt^2} = -\frac{g}{L}\theta
これは単振動の式であり、角振動数 ω=gL\omega = \sqrt{\frac{g}{L}} である。周期 T0T_0
T0=2πω=2πLgT_0 = \frac{2\pi}{\omega} = 2\pi\sqrt{\frac{L}{g}}
次に、y0=12at2+V0ty_0 = \frac{1}{2}at^2 + V_0 t の場合:
慣性系で考えると、重力加速度は g=g+ag' = g + a となる。したがって、
Ta=2πLg+aT_a = 2\pi\sqrt{\frac{L}{g+a}}
(2) 減衰振動の場合:
運動方程式は、
mLd2θdt2=mgsinθcLdθdtmL\frac{d^2\theta}{dt^2} = -mg\sin\theta - cL\frac{d\theta}{dt}
y0=V0ty_0 = V_0t で、座標系が下降しているので、慣性力は働かない。
θ\theta が小さいとき sinθθ\sin\theta \approx \theta なので、
mLd2θdt2+cLdθdt+mgθ=0mL\frac{d^2\theta}{dt^2} + cL\frac{d\theta}{dt} + mg\theta = 0
d2θdt2=cmdθdtgLθ\frac{d^2\theta}{dt^2} = -\frac{c}{m}\frac{d\theta}{dt} - \frac{g}{L}\theta
したがって、P=gLP = -\frac{g}{L}, Q=cmQ = -\frac{c}{m}, R=0R = 0
(3) 減衰振動の条件:
d2θdt2+cmdθdt+gLθ=0\frac{d^2\theta}{dt^2} + \frac{c}{m}\frac{d\theta}{dt} + \frac{g}{L}\theta = 0
この微分方程式の解は、
θ(t)=Aeλt\theta(t) = Ae^{\lambda t}
とおくと、λ2+cmλ+gL=0\lambda^2 + \frac{c}{m}\lambda + \frac{g}{L} = 0
λ=cm±(cm)24gL2\lambda = \frac{-\frac{c}{m} \pm \sqrt{(\frac{c}{m})^2 - 4\frac{g}{L}}}{2}
減衰振動の条件は、(cm)2<4gL(\frac{c}{m})^2 < 4\frac{g}{L} である。
c2<4m2gLc^2 < 4m^2\frac{g}{L}
c<2mgLc < 2m\sqrt{\frac{g}{L}}
したがって、c0=2mgLc_0 = 2m\sqrt{\frac{g}{L}}
(4) 特定の条件でのΘ(t)\Theta(t)
c=mgLc = m\sqrt{\frac{g}{L}}, V0=12gLV_0 = \frac{1}{2}\sqrt{gL} のとき、
λ=cm±(cm)24gL2=gL±gL4gL2=gL±i3gL2\lambda = \frac{-\frac{c}{m} \pm \sqrt{(\frac{c}{m})^2 - 4\frac{g}{L}}}{2} = \frac{-\sqrt{\frac{g}{L}} \pm \sqrt{\frac{g}{L} - 4\frac{g}{L}}}{2} = \frac{-\sqrt{\frac{g}{L}} \pm i\sqrt{3\frac{g}{L}}}{2}
λ=12gL±i32gL\lambda = -\frac{1}{2}\sqrt{\frac{g}{L}} \pm i\frac{\sqrt{3}}{2}\sqrt{\frac{g}{L}}
したがって、ω=32gL\omega = \frac{\sqrt{3}}{2}\sqrt{\frac{g}{L}}
Θ(t)=e12gLt[Acos(ωt)+Bsin(ωt)]\Theta(t) = e^{-\frac{1}{2}\sqrt{\frac{g}{L}}t} [A\cos(\omega t) + B\sin(\omega t)]
初期条件 Θ(0)=0\Theta(0) = 0 より、A=0A=0
Θ(t)=Be12gLtsin(ωt)\Theta(t) = Be^{-\frac{1}{2}\sqrt{\frac{g}{L}}t} \sin(\omega t)
dΘdt(0)=U0L=B32gLB=2U03Lg\frac{d\Theta}{dt}(0) = \frac{U_0}{L} = \frac{B\sqrt{3}}{2}\sqrt{\frac{g}{L}} \Rightarrow B = \frac{2 U_0}{\sqrt{3}}\sqrt{\frac{L}{g}}
Θ(t)=2U03Lge12gLtsin(32gLt)\Theta(t) = \frac{2U_0}{\sqrt{3}}\sqrt{\frac{L}{g}} e^{-\frac{1}{2}\sqrt{\frac{g}{L}}t} \sin(\frac{\sqrt{3}}{2}\sqrt{\frac{g}{L}}t)

3. 最終的な答え

(1) T0=2πLgT_0 = 2\pi\sqrt{\frac{L}{g}}, Ta=2πLg+aT_a = 2\pi\sqrt{\frac{L}{g+a}}
(2) P=gLP = -\frac{g}{L}, Q=cmQ = -\frac{c}{m}, R=0R = 0
(3) c0=2mgLc_0 = 2m\sqrt{\frac{g}{L}}
(4) Θ(t)=2U03Lge12gLtsin(32gLt)\Theta(t) = \frac{2U_0}{\sqrt{3}}\sqrt{\frac{L}{g}} e^{-\frac{1}{2}\sqrt{\frac{g}{L}}t} \sin(\frac{\sqrt{3}}{2}\sqrt{\frac{g}{L}}t)

「応用数学」の関連問題

角度 $\theta = 30^\circ$ のなめらかな斜面において、高さ $h = 0.50$ m の位置に質量 $m = 3.0$ kg の物体を置く。手を離すと、物体は斜面を滑り降りる。物体が...

力学エネルギー保存則摩擦
2025/7/18

99mTcに関する以下の3つの問題を解きます。 (1) 99mTcの壊変定数を求め、壊変によって放出される放射線の名前を答えます。 (2) 2.6 ngの99mTcを含む注射剤の放射能を求めます。 (...

放射線半減期指数関数壊変定数放射能
2025/7/18

問題は、発散するベクトル場 $\vec{A}$ に対して、$\mathrm{div} \vec{A} = \rho$ が成り立つ理由を説明せよ、というものです。ここで、$\mathrm{div} \v...

ベクトル解析発散ガウスの発散定理連続の式電磁気学
2025/7/18

片持ち梁の先端に集中荷重 $P$ が作用する場合のせん断力図 (SFD) と曲げモーメント図 (BMD) を描き、最大曲げ応力 $\sigma_{max}$ と最大せん断応力 $\tau_{max}$...

構造力学せん断力図曲げモーメント図応力
2025/7/18

一定の時間間隔でパルスを発生させる音源があり、観測者の運動状態に応じて観測されるパルスの振動数が変化するドップラー効果に関する問題です。空欄 【23】 から 【27】 に当てはまる選択肢を選びます。

物理ドップラー効果速度振動数
2025/7/18

与えられたポテンシャル $U(x, y, z)$ に対して、保存力 $\vec{F}$ の各成分 $F_x, F_y, F_z$ を求める問題です。 保存力 $\vec{F}$ はポテンシャル $U$...

勾配偏微分ポテンシャル保存力
2025/7/18

問題は、ベクトルに関する恒等式と全微分の式を証明することです。 a) $A \cdot (B \times C) = (A \times B) \cdot C$ を証明する。 b) $dz = \na...

ベクトルベクトル解析全微分勾配
2025/7/18

質量 $m$ のおもりが長さ $L$ の棒で吊るされた振り子を考える。支点は鉛直下方に等速 $V_0$ で下降しており、時刻 $t$ での支点の座標 $y_0$ は $y_0 = V_0 t$ で表さ...

力学振り子微分方程式減衰振動
2025/7/18

質量 $m$ のおもりが長さ $L$ の棒でぶら下げられた振り子について、以下の問題を解きます。 (1) 抵抗がない場合 ($c=0$) の振り子の周期 $T_0$ を求めます。次に、支点の $y$ ...

物理振り子運動方程式微分方程式減衰振動
2025/7/18

質量 $m$ のおもりが長さ $L$ の棒で吊り下げられた振り子を考えます。支点は鉛直下方に等速 $V_0$ で下降し、時刻 $t$ での支点の座標 $y_0$ は $y_0 = V_0t$ で表され...

力学振り子微分方程式減衰振動
2025/7/18