跳到主要内容

早稲田大学 基幹理工学研究科 電子物理システム学専攻 2023年8月実施 力学 (その1) [1]

Author​

Miyake

Description​

問題の要約 — 大学公表問題

問1​

質量を無視できる半径 aa の円環に質量 mm の質点を固定し、円環を直線 y=−ay=-a に沿って滑らずに転がす。初めの中心は原点、質点は (0,−a)(0,-a) にある。下図では右向きの転がり角を θ\theta とする。重力による位置エネルギーを mgymgy とし、摩擦による散逸は考えない。

  1. 質点の位置を θ\theta で表す。
  2. θ\theta を一般化座標とするラグランジアンを求める。
  3. 共役運動量とハミルトニアンを求める。
  4. 正準方程式を書く。

問2​

長さ aa の質量のない棒の上端から 2a/32a/3 の位置に質量 mm の質点を固定する。棒の両端はそれぞれ壁と床に接したまま滑り、鉛直な壁との角を φ\varphi とする。重力は下向き、摩擦は無視する。

  1. 質点の位置を φ(t)\varphi(t) で表す。
  2. ラグランジアンを求める。
  3. φ\varphi の運動方程式を導く。

円環と壁に沿って滑る棒の座標を示す独自の模式図

题目描述​

【问题 1】

在半径为 aa、质量可忽略的圆环上固定质量为 mm 的质点,使圆环沿直线 y=−ay=-a 无滑动滚动。初始圆心为原点、质点位于 (0,−a)(0,-a),向右滚动角为 θ\theta,重力势能为 mgymgy,忽略摩擦耗散。

  1. 用 θ\theta 写出质点坐标,几何关系为 X=aθX=a\theta、x=aθ−asin⁡θx=a\theta-a\sin\theta、y=−acos⁡θy=-a\cos\theta。
  2. 求 Lagrangian L(θ,θ˙)L(\theta,\dot\theta)。
  3. 求共轭动量 pθp_\theta 和 Hamiltonian H(θ,pθ)H(\theta,p_\theta)。
  4. 写出 Hamilton 正则方程。

【问题 2】

长度为 aa 的无质量杆上,距上端 2a/32a/3 处固定质量为 mm 的质点。杆的两端分别贴着竖直墙面与水平地面无摩擦滑动,杆与竖直墙面的夹角为 φ\varphi,重力向下。接触点坐标满足 X=asin⁡φX=a\sin\varphi、Y=acos⁡φY=a\cos\varphi,质点坐标为 x=2X/3x=2X/3、y=Y/3y=Y/3。

  1. 用 φ(t)\varphi(t) 表示质点坐标。
  2. 求 L(φ,φ˙)L(\varphi,\dot\varphi)。
  3. 用 Euler–Lagrange 方程导出 φ\varphi 的运动方程。

Kai​

問 1​

(1)​

円環の中心の x 座標は

X=aθ \begin{aligned} X = a \theta \end{aligned}

であり、質点の座標は

x=X−asin⁡θ=aθ−asin⁡θ,y=−acos⁡θ \begin{aligned} x &= X - a \sin \theta = a \theta - a \sin \theta ,\\ y &= - a \cos \theta \end{aligned}

である。

(2)​

x˙=aθ˙−aθ˙cos⁡θ,y˙=aθ˙sin⁡θ \begin{aligned} \dot{x} &= a \dot{\theta} - a \dot{\theta} \cos \theta , \\ \dot{y} &= a \dot{\theta} \sin \theta \end{aligned}

なので、求めるラグランジアンは

L(θ,θ˙)=12m(x˙2 +y˙2)−mgy=ma2θ˙2(1−cos⁡θ)+mgacos⁡θ \begin{aligned} L \left( \theta, \dot{\theta} \right) &= \frac{1}{2} m \left( \dot{x}^2\ + \dot{y}^2 \right) - mgy \\ &= ma^2 \dot{\theta}^2 ( 1 - \cos \theta ) + mga \cos \theta \end{aligned}

である。

(3)​

以下のルジャンドル変換と正準方程式は 1−cos⁡θ>01-\cos\theta>0 の範囲で成り立つ。θ=2kπ\theta=2k\pi では速度の係数が零になり、この座標による変換は特異になる。

pθ=∂L∂θ˙=2ma2θ˙(1−cos⁡θ),H(θ,pθ)=θ˙pθ−L=2ma2θ˙2(1−cos⁡θ)−ma2θ˙2(1−cos⁡θ)−mgacos⁡θ=ma2θ˙2(1−cos⁡θ)−mgacos⁡θ=pθ24ma2(1−cos⁡θ)−mgacos⁡θ \begin{aligned} p_\theta &= \frac{\partial L}{\partial \dot{\theta}} \\ &= 2ma^2 \dot{\theta} ( 1 - \cos \theta ) ,\\ H \left( \theta, p_\theta \right) &= \dot{\theta} p_\theta - L \\ &= 2ma^2 \dot{\theta}^2 ( 1 - \cos \theta ) - ma^2 \dot{\theta}^2 ( 1 - \cos \theta ) - mga \cos \theta \\ &= ma^2 \dot{\theta}^2 ( 1 - \cos \theta ) - mga \cos \theta \\ &= \frac{p_\theta^2}{4ma^2 ( 1 - \cos \theta )} - mga \cos \theta \end{aligned}

(4)​

θ˙=∂H∂pθ=pθ2ma2(1−cos⁡θ),p˙θ=−∂H∂θ=pθ2sin⁡θ4ma2(1−cos⁡θ)2−mgasin⁡θ \begin{aligned} \dot{\theta} &= \frac{\partial H}{\partial p_\theta} = \frac{p_\theta}{2ma^2(1 - \cos \theta)} ,\\ \dot{p}_\theta &= - \frac{\partial H}{\partial \theta} = \frac{p_\theta^2 \sin \theta}{4ma^2 (1 - \cos \theta)^2} - mga \sin \theta \end{aligned}

問 2​

(1)​

時刻 tt において、棒と床が接触している点の x 座標を X(t)X(t) とし、 棒と壁が接触している点の y 座標を Y(t)Y(t) とすると、

X(t)=asin⁡φ(t),  Y(t)=acos⁡φ(t) \begin{aligned} X(t) = a \sin \varphi (t) , \ \ Y(t) = a \cos \varphi (t) \end{aligned}

なので、

x(t)=23X(t)=23asin⁡φ(t),y(t)=13Y(t)=13acos⁡φ(t) \begin{aligned} x(t) &= \frac{2}{3} X(t) = \frac{2}{3} a \sin \varphi (t) ,\\ y(t) &= \frac{1}{3} Y(t) = \frac{1}{3} a \cos \varphi (t) \end{aligned}

である。

(2)​

x˙(t)=23aφ˙(t)cos⁡φ(t),y˙(t)=−13aφ˙(t)sin⁡φ(t) \begin{aligned} \dot{x} (t) &= \frac{2}{3} a \dot{\varphi}(t) \cos \varphi (t) ,\\ \dot{y} (t) &= - \frac{1}{3} a \dot{\varphi}(t) \sin \varphi (t) \end{aligned}

なので、求めるラグランジアンは

L(φ,φ˙)=12m(x˙2+y˙2)−mgy=118ma2φ˙2(4cos⁡2φ+sin⁡2φ)−13mgacos⁡φ=118ma2φ˙2(3cos⁡2φ+1)−13mgacos⁡φ \begin{aligned} L \left( \varphi, \dot{\varphi} \right) &= \frac{1}{2} m \left( \dot{x}^2 + \dot{y}^2 \right) - mgy \\ &= \frac{1}{18} m a^2 \dot{\varphi}^2 \left( 4 \cos^2 \varphi + \sin^2 \varphi \right) - \frac{1}{3} mga \cos \varphi \\ &= \frac{1}{18} m a^2 \dot{\varphi}^2 \left( 3 \cos^2 \varphi + 1 \right) - \frac{1}{3} mga \cos \varphi \end{aligned}

である。

(3)​

ddt∂L∂φ˙=ddt19ma2φ˙(3cos⁡2φ+1)=19ma2(φ¨(3cos⁡2φ+1)−6φ˙2cos⁡φsin⁡φ),∂L∂φ=−13ma2φ˙2cos⁡φsin⁡φ+13mgasin⁡φ \begin{aligned} \frac{d}{dt} \frac{\partial L}{\partial \dot{\varphi}} &= \frac{d}{dt} \frac{1}{9} m a^2 \dot{\varphi} \left( 3 \cos^2 \varphi + 1 \right) \\ &= \frac{1}{9} m a^2 \left( \ddot{\varphi} \left( 3 \cos^2 \varphi + 1 \right) - 6 \dot{\varphi}^2 \cos \varphi \sin \varphi \right) ,\\ \frac{\partial L}{\partial \varphi} &= - \frac{1}{3} ma^2 \dot{\varphi}^2 \cos \varphi \sin \varphi + \frac{1}{3} mga \sin \varphi \end{aligned}

なので、求める運動方程式は

φ¨(3cos⁡2φ+1)−6φ˙2cos⁡φsin⁡φ=−3φ˙2cos⁡φsin⁡φ+3gasin⁡φ \begin{aligned} \ddot{\varphi} \left( 3 \cos^2 \varphi + 1 \right) - 6 \dot{\varphi}^2 \cos \varphi \sin \varphi &= - 3 \dot{\varphi}^2 \cos \varphi \sin \varphi + 3 \frac{g}{a} \sin \varphi \end{aligned}
∴  φ¨(3cos⁡2φ+1)−3φ˙2cos⁡φsin⁡φ−3gasin⁡φ=0 \begin{aligned} \therefore \ \ \ddot{\varphi} \left( 3 \cos^2 \varphi + 1 \right) - 3 \dot{\varphi}^2 \cos \varphi \sin \varphi - 3 \frac{g}{a} \sin \varphi &= 0 \end{aligned}

である。