跳到主要内容

早稲田大学 基幹理工学研究科 電子物理システム学専攻 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\thetax=aθasinθx=a\theta-a\sin\thetay=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\varphiY=acosφY=a\cos\varphi,质点坐标为 x=2X/3x=2X/3y=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=Xasinθ=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(1cosθ)+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)

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

pθ=Lθ˙=2ma2θ˙(1cosθ),H(θ,pθ)=θ˙pθL=2ma2θ˙2(1cosθ)ma2θ˙2(1cosθ)mgacosθ=ma2θ˙2(1cosθ)mgacosθ=pθ24ma2(1cosθ)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)

θ˙=Hpθ=pθ2ma2(1cosθ),p˙θ=Hθ=pθ2sinθ4ma2(1cosθ)2mgasinθ \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(4cos2φ+sin2φ)13mgacosφ=118ma2φ˙2(3cos2φ+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)

ddtLφ˙=ddt19ma2φ˙(3cos2φ+1)=19ma2(φ¨(3cos2φ+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}

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

φ¨(3cos2φ+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}
  φ¨(3cos2φ+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}

である。