跳到主要内容

大阪大学 理学研究科 物理学専攻 2019年8月実施 問題1

Author

Miyake

Description

題意の要約(暫定)

質量 mm、長さ ll の単振り子を考える。重力加速度を gg とし、振れ角 θ\theta は十分小さいとする。初期条件は

θ(0)=α,θ˙(0)=0.\theta(0)=\alpha,\qquad\dot\theta(0)=0.

である。

【I:微小振動】

  1. 運動方程式から θ(t)\theta(t) と角振動数 ω\omega を求めよ。
  2. 振動の力学的エネルギー EE を求めよ。
  3. θ2\theta^2 の項まで残す近似の下で、糸の張力 T(t)T(t) を求めよ。
  4. 一周期にわたる平均張力 T\langle T\rangle を求めよ。

【II:振り子の長さの緩やかな変化】

振り子の長さを緩やかに δl\delta l だけ変化させ、微小変化の一次の項のみを考える。

  1. 角振動数の変化 δω\delta\omega を求めよ。
  2. 平均張力による仕事から、力学的エネルギーの変化 δE\delta E を求めよ。
  3. この緩やかな変化に対し、
    E+δEω+δω=Eω+O(δl2),\frac{E+\delta E}{\omega+\delta\omega} =\frac E\omega+O(\delta l^2),
    すなわち E/ωE/\omega が断熱不変量であることを示せ。

【III:位相空間】

  1. 微小振動する単振り子の Lagrangian を書き、θ\theta に共役な運動量 pp を求めよ。
  2. エネルギー EE を固定したとき、位相平面 (θ,p)(\theta,p) における楕円軌道の方程式を書け。
  3. 作用積分 pdθ,\oint p\,d\theta, を求め、振り子の長さを緩やかに変えるときに不変であることを説明せよ。

题目描述

考虑质量为 mm 的单摆,摆长为 ll,重力加速度为 gg,摆角 θ\theta 很小;初始条件为 θ(0)=α,θ˙(0)=0.\theta(0)=\alpha,\qquad\dot\theta(0)=0.

【I:小振动】

  1. 由运动方程求 θ(t)\theta(t) 与角频率 ω\omega
  2. 求振动的机械能 EE
  3. 在保留至 θ2\theta^2 的近似下,求摆线张力 T(t)T(t)
  4. 求一个周期内的平均张力 T\langle T\rangle

【II:缓慢改变摆长】

将摆长缓慢改变 δl\delta l,只保留一阶小量。

  1. 求角频率变化 δω\delta\omega
  2. 由平均张力所做功求机械能变化 δE\delta E
  3. 证明在这一缓慢变化下
    E+δEω+δω=Eω+O(δl2),\frac{E+\delta E}{\omega+\delta\omega} =\frac E\omega+O(\delta l^2),
    E/ωE/\omega 为绝热不变量。

【III:相空间】

  1. 写出小角度单摆的 Lagrangian,并求与 θ\theta 共轭的动量 pp
  2. 在固定能量 EE 下,写出相平面 (θ,p)(\theta,p) 中轨道的椭圆方程。
  3. 计算作用量 pdθ,\oint p\,d\theta, 并说明它在摆长缓慢变化时保持不变。

Kai

I.

(1)

θ¨=dθ˙/dt\ddot{\theta} = d \dot{\theta} / dt と書くと、運動方程式は次のようになる:

mlθ¨=mgsinθ  θ¨=glsinθglθ\begin{aligned} ml \ddot{\theta} &= - mg \sin \theta \\ \therefore \ \ \ddot{\theta} &= - \frac{g}{l} \sin \theta \\ &\simeq - \frac{g}{l} \theta \end{aligned}

よって、 ω=g/l\omega = \sqrt{g/l} とおいて、

θ(t)=αcosωt=αcos(glt)\begin{aligned} \theta (t) &= \alpha \cos \omega t \\ &= \alpha \cos \left( \sqrt{\frac{g}{l}} t \right) \end{aligned}

が与えられた初期条件を満たす解である。

(2)

θ˙(t)=αωsinωt\dot{\theta}(t) = - \alpha \omega \sin \omega t の最大値は αω\alpha \omega であるから、

E=12ml2α2ω2=12mglα2\begin{aligned} E &= \frac{1}{2} m l^2 \alpha^2 \omega^2 \\ &= \frac{1}{2} mgl \alpha^2 \end{aligned}

である。

(3)

糸の方向の運動方程式より、

mlθ˙2=Tmgcosθ  T=mlθ˙2+mgcosθmlθ˙2+mg12mgθ2\begin{aligned} ml \dot{\theta}^2 &= T - mg \cos \theta \\ \therefore \ \ T &= ml \dot{\theta}^2 + mg \cos \theta \\ &\simeq ml \dot{\theta}^2 + mg - \frac{1}{2} mg \theta^2 \end{aligned}

がわかる。

(4)

(1), (3) より、

T=mlα2ω2sin2ωt+mg12mgα2cos2ωt=mgα2sin2ωt+mg12mgα2cos2ωt=mg(α2+443α24cos2ωt)  T=α2+44mg\begin{aligned} T &= ml \alpha^2 \omega^2 \sin^2 \omega t + mg - \frac{1}{2} mg \alpha^2 \cos^2 \omega t \\ &= mg \alpha^2 \sin^2 \omega t + mg - \frac{1}{2} mg \alpha^2 \cos^2 \omega t \\ &=mg \left( \frac{\alpha^2 + 4}{4} - \frac{3\alpha^2}{4} \cos 2 \omega t \right) \\ \therefore \ \ \langle T \rangle &= \frac{\alpha^2 + 4}{4} mg \end{aligned}

がわかる。

II.

(5)

δω=gl+δlω=ω(1+δll)1/2ωω(1δl2l)ω=δl2lω\begin{aligned} \delta \omega &= \sqrt{\frac{g}{l + \delta l}} - \omega \\ &= \omega \left( 1 + \frac{\delta l}{l} \right)^{-1/2} - \omega \\ &\simeq \omega \left( 1 - \frac{\delta l}{2l} \right) - \omega \\ &= - \frac{\delta l}{2l} \omega \end{aligned}

(6)

δW=Tδl=α2+44mgδl\begin{aligned} \delta W &= - \langle T \rangle \delta l \\ &= - \frac{\alpha^2 + 4}{4} mg \delta l \end{aligned}

また、 δW=δEmgδl\delta W = \delta E - mg \delta l より、

δE=δW+mgδl=α24mgδl\begin{aligned} \delta E &= \delta W + mg \delta l \\ &= - \frac{\alpha^2}{4} mg \delta l \end{aligned}

である。

(7)

(2), (5), (6) の結果を使って、 δl\delta l の1次までで次のように計算できる:

E+δEω+δω=Eω(1+δEE)(1+δωω)1Eω(1+δEEδωω)=Eω(1δl2l+δl2l)=Eω\begin{aligned} \frac{E + \delta E}{\omega + \delta \omega} &= \frac{E}{\omega} \left(1 + \frac{\delta E}{E} \right) \left(1 + \frac{\delta \omega}{\omega} \right)^{-1} \\ &\simeq \frac{E}{\omega} \left(1 + \frac{\delta E}{E} - \frac{\delta \omega}{\omega} \right) \\ &= \frac{E}{\omega} \left(1 - \frac{\delta l}{2l} + \frac{\delta l}{2l} \right) \\ &= \frac{E}{\omega} \end{aligned}

III.

(8)

L=12ml2θ˙2mgl(1cosθ)  p=Lθ˙=ml2θ˙\begin{aligned} L &= \frac{1}{2} ml^2 \dot{\theta}^2 - mgl (1 - \cos \theta) \\ \therefore \ \ p &= \frac{\partial L}{\partial \dot{\theta}} \\ &= ml^2 \dot{\theta} \end{aligned}

(9)

E=12ml2θ˙2+mgl(1cosθ)12ml2θ˙2+12mglθ2=p22ml2+12mglθ2  p22ml2E+θ22Emgl=1\begin{aligned} E &= \frac{1}{2} ml^2 \dot{\theta}^2 + mgl (1 - \cos \theta) \\ &\simeq \frac{1}{2} ml^2 \dot{\theta}^2 + \frac{1}{2} mgl \theta^2 \\ &= \frac{p^2}{2ml^2} + \frac{1}{2} mgl \theta^2 \\ \therefore \ \ &\frac{p^2}{2ml^2E} + \frac{\theta^2}{\frac{2E}{mgl}} = 1 \end{aligned}

これは楕円の方程式である。

(10)

(9) より、

pdθ=π2ml2E2Emgl=2πElg=2πEω\begin{aligned} \oint p d \theta &= \pi \sqrt{2ml^2E} \sqrt{\frac{2E}{mgl}} \\ &= 2 \pi E \sqrt{\frac{l}{g}} \\ &= 2 \pi \frac{E}{\omega} \end{aligned}

となるが、 (7) より、これは糸の長さが δl\delta l 変化する前後で変わらない。