名古屋大学 理学研究科 物理学教室 2019年8月実施 物理学 [I]
Author
Miyake
Description
題意の要約。
原点 O O O を軸とする半径 R R R の円柱に、質量を無視できる長さ 2 π R 2\pi R 2 π R の伸びない糸を巻く。一端を円柱上の S S S 、他端を質量 m m m の質点 P P P に結ぶ。初め S S S と P P P は y y y 軸上の S 0 S_0 S 0 にあり、t = 0 t=0 t = 0 に質点へ + y +y + y 方向の初速 v 0 > 0 v_0>0 v 0 > 0 を与える。糸が張ったままほどける範囲を考え、ポテンシャルエネルギーはない。接点を C C C 、時計回りを正として θ = ∠ S O C \theta=\angle SOC θ = ∠ SOC とする。
円柱を固定した場合:
O C → = ( R sin θ , R cos θ ) \overrightarrow{OC}=(R\sin\theta,R\cos\theta) OC = ( R sin θ , R cos θ ) から C P → \overrightarrow{CP} CP と O P → \overrightarrow{OP} OP を求める。
t > 0 t>0 t > 0 の速度成分を R , θ , θ ˙ R,\theta,\dot\theta R , θ , θ ˙ で表す。
Lagrangian を求める。
θ \theta θ の Euler–Lagrange 方程式を求める。
方程式を解いて θ ( t ) = 2 v 0 t / R \theta(t)=\sqrt{2v_0t/R} θ ( t ) = 2 v 0 t / R を導く。
O O O まわりの角運動量の大きさと、その保存性を理由とともに示す。
質点の運動エネルギーの保存性を理由とともに示す。
次に、円柱が軸まわりに摩擦なく回転できる場合を考える。時計回りの回転角を ψ = ∠ S 0 O S \psi=\angle S_0OS ψ = ∠ S 0 OS 、慣性モーメントを I > 0 I>0 I > 0 とする。円柱は初め静止し、質点の初期位置と初速度は前と同じとする。
円柱の回転と質点の並進を含む Lagrangian を求める。
θ , ψ \theta,\psi θ , ψ の Euler–Lagrange 方程式を求める。
ψ \psi ψ の方程式から保存量を求める。
θ ˙ > 0 \dot\theta>0 θ ˙ > 0 に注意し、円柱の回転方向を理由とともに示す。
質点の運動エネルギーの保存性を理由とともに示す。
题目描述
半径为 R R R 的圆柱上缠有质量可忽略、长度 2 π R 2\pi R 2 π R 的不可伸长细线,线端质点 P P P 的质量为 m m m 。一端系在圆柱上的 S S S ,另一端接质点。初始时 S , P S,P S , P 位于 y y y 轴上的 S 0 S_0 S 0 ,给质点沿 + y +y + y 的初速度 v 0 > 0 v_0>0 v 0 > 0 。令 θ = ∠ S O C \theta=\angle SOC θ = ∠ SOC ,以顺时针为正,接触点为 C C C ,圆心为 O;考虑细线保持张紧并展开的范围,忽略势能与耗散。
固定圆柱(问题 1—7) :
求 O P → \overrightarrow{OP} OP 关于 θ \theta θ 的坐标,说明轨迹为圆的渐开线。
求质点速度 ( v x , v y ) (v_x,v_y) ( v x , v y ) 。
写出 Lagrangian。
推导 θ \theta θ 的 Euler–Lagrange 运动方程。
在初始条件 t = 0 t=0 t = 0 时 θ = 0 , v x = 0 , v y = v 0 \theta=0,\ v_x=0,\ v_y=v_0 θ = 0 , v x = 0 , v y = v 0 下,求解第4问的微分方程,导出
θ ( t ) = 2 v 0 t R \theta(t)=\sqrt{\frac{2v_0t}{R}} θ ( t ) = R 2 v 0 t
。
求质点关于 O 的角动量大小,并判断是否守恒。其大小可写为 m 2 R v 0 3 t m\sqrt{2Rv_0^3t} m 2 R v 0 3 t 。
判断质点动能是否守恒。
可转动圆柱(问题 8—12) :现在圆柱可绕 O 转动,顺时针转角为 ψ = ∠ S 0 O S \psi=\angle S_0OS ψ = ∠ S 0 OS 、转动惯量为 I > 0 I>0 I > 0 ;质点相对圆柱展开角仍为 θ \theta θ 。
求质点位置、速度和含圆柱转动能的 Lagrangian。
写出关于 θ , ψ \theta,\psi θ , ψ 的两条 Euler–Lagrange 方程。
利用 ψ \psi ψ 是循环坐标求相应守恒广义动量。
在 t = 0 t=0 t = 0 时 θ = 0 , ψ ˙ = 0 \theta=0,\dot\psi=0 θ = 0 , ψ ˙ = 0 且随后 θ , θ ˙ > 0 \theta,\dot\theta>0 θ , θ ˙ > 0 的条件下,判断圆柱转动方向并说明理由;ψ ˙ < 0 \dot\psi<0 ψ ˙ < 0 对应逆时针转动。
判断此时质点自身的动能是否守恒,并说明能量转移。
Kai
問 1.
C P → = R θ ( − cos θ , sin θ ) ∴ O P → = O C → + C P → = R ( sin θ − θ cos θ , cos θ + θ sin θ ) \begin{aligned}
\overrightarrow{\mathrm{CP}}
&= R \theta (- \cos \theta, \sin \theta)
\\
\therefore \ \
\overrightarrow{\mathrm{OP}}
&= \overrightarrow{\mathrm{OC}} + \overrightarrow{\mathrm{CP}}
\\
&= R (\sin \theta - \theta \cos \theta, \cos \theta + \theta \sin \theta)
\end{aligned} CP ∴ OP = Rθ ( − cos θ , sin θ ) = OC + CP = R ( sin θ − θ cos θ , cos θ + θ sin θ )
問 2.
( v x , v y ) = R θ θ ˙ ( sin θ , cos θ ) \begin{aligned}
(v_x, v_y)
= R \theta \dot{\theta} (\sin \theta, \cos \theta)
\end{aligned} ( v x , v y ) = Rθ θ ˙ ( sin θ , cos θ )
問 3.
L = m 2 ( v x 2 + v y 2 ) = m 2 R 2 θ 2 θ ˙ 2 \begin{aligned}
\mathcal{L}
&= \frac{m}{2} \left( v_x^2 + v_y^2 \right)
= \frac{m}{2} R^2 \theta^2 \dot{\theta}^2
\end{aligned} L = 2 m ( v x 2 + v y 2 ) = 2 m R 2 θ 2 θ ˙ 2
問 4.
d d t ∂ L ∂ θ ˙ = d d t m R 2 θ 2 θ ˙ = m R 2 ( 2 θ θ ˙ 2 + θ 2 θ ¨ ) ∂ L ∂ θ = m R 2 θ θ ˙ 2 \begin{aligned}
\frac{\mathrm{d}}{\mathrm{d}t}
\frac{\partial \mathcal{L}}{\partial \dot{\theta}}
&= \frac{\mathrm{d}}{\mathrm{d}t} m R^2 \theta^2 \dot{\theta}
= m R^2 \left( 2 \theta \dot{\theta}^2 + \theta^2 \ddot{\theta} \right)
\\
\frac{\partial \mathcal{L}}{\partial \theta}
&= m R^2 \theta \dot{\theta}^2
\end{aligned} d t d ∂ θ ˙ ∂ L ∂ θ ∂ L = d t d m R 2 θ 2 θ ˙ = m R 2 ( 2 θ θ ˙ 2 + θ 2 θ ¨ ) = m R 2 θ θ ˙ 2
なので、オイラー-ラグランジュ方程式より、
θ θ ˙ 2 + θ 2 θ ¨ = 0 \begin{aligned}
\theta \dot{\theta}^2 + \theta^2 \ddot{\theta} = 0
\end{aligned} θ θ ˙ 2 + θ 2 θ ¨ = 0
問 5.
t > 0 t>0 t > 0 で θ > 0 \theta>0 θ > 0 だから、問4を θ \theta θ で割ると
d d t ( θ θ ˙ ) = 0. \frac{d}{dt}(\theta\dot\theta)=0. d t d ( θ θ ˙ ) = 0.
初速度の条件より R θ θ ˙ = v 0 R\theta\dot\theta=v_0 Rθ θ ˙ = v 0 。これを積分し、θ ( 0 ) = 0 \theta(0)=0 θ ( 0 ) = 0 を使うと
θ 2 = 2 v 0 t R , θ ( t ) = 2 v 0 t R . \theta^2=\frac{2v_0t}{R},\qquad
\theta(t)=\sqrt{\frac{2v_0t}{R}}. θ 2 = R 2 v 0 t , θ ( t ) = R 2 v 0 t .
与えられた θ ( t ) = 2 v 0 t / R \theta(t) = \sqrt{2 v_0 t / R} θ ( t ) = 2 v 0 t / R が
運動方程式と初期条件を満たすことは次のようにして確かめられる。
まず、
θ ˙ ( t ) = v 0 2 R t , θ ¨ ( t ) = − v 0 8 R t 3 \begin{aligned}
\dot{\theta} (t) = \sqrt{\frac{v_0}{2Rt}}
, \ \
\ddot{\theta} (t) = - \sqrt{\frac{v_0}{8Rt^3}}
\end{aligned} θ ˙ ( t ) = 2 Rt v 0 , θ ¨ ( t ) = − 8 R t 3 v 0
より、運動方程式 θ θ ˙ 2 + θ 2 θ ¨ = 0 \theta \dot{\theta}^2 + \theta^2 \ddot{\theta} = 0 θ θ ˙ 2 + θ 2 θ ¨ = 0
を満たすことがわかる。
次に、
( v x , v y ) = v 0 ( sin 2 v 0 t R , cos 2 v 0 t R ) \begin{aligned}
(v_x, v_y)
= v_0 \left(\sin \sqrt{\frac{2v_0t}{R}}, \cos \sqrt{\frac{2v_0t}{R}} \right)
\end{aligned} ( v x , v y ) = v 0 ( sin R 2 v 0 t , cos R 2 v 0 t )
より、初期条件( t = 0 t=0 t = 0 で θ = 0 , v x = 0 , v y = v 0 \theta = 0, \ v_x = 0, \ v_y = v_0 θ = 0 , v x = 0 , v y = v 0
)も t → 0 + t\to0^+ t → 0 + の極限で満たすことがわかる。θ ˙ \dot\theta θ ˙ 自体はこの極限で発散するが、質点の速度は有限である。
問 6.
点 P の位置を ( x , y ) (x,y) ( x , y ) とすると、求める角運動量の大きさは、
L = m ∣ x v y − y v x ∣ = m R 2 θ 2 ∣ θ ˙ ∣ \begin{aligned}
L
&= m \left| x v_y - y v_x \right|
\\
&= m R^2 \theta^2 \left| \dot{\theta} \right|
\end{aligned} L = m ∣ x v y − y v x ∣ = m R 2 θ 2 θ ˙
であるが、問5で与えられた解を使うと、
L = m 2 R v 0 3 t \begin{aligned}
L
&= m \sqrt{2R v_0^3 t}
\end{aligned} L = m 2 R v 0 3 t
を得る。
これは t t t に依存するので保存量ではない。
(糸の一端がy軸上に固定されていることから、この系は中心軸に関する回転対称性を持たない。)
問 7.
問 5 より、 v x 2 + v y 2 = v 0 2 v_x^2 + v_y^2 = v_0^2 v x 2 + v y 2 = v 0 2 なので、運動エネルギーは保存する。
(摩擦や非弾性衝突などがなく、ポテンシャルエネルギーも考えていないので、運動エネルギーが保存しない理由がない。)
問 8.
問 1, 2, 3 と同じように考えて、
O C → = R ( sin ( θ + ψ ) , cos ( θ + ψ ) ) C P → = R θ ( − cos ( θ + ψ ) , sin ( θ + ψ ) ) O P → = O C → + C P → = R ( sin ( θ + ψ ) − θ cos ( θ + ψ ) , cos ( θ + ψ ) + θ sin ( θ + ψ ) ) v ⃗ = R ( ψ ˙ cos ( θ + ψ ) + θ ( θ ˙ + ψ ˙ ) sin ( θ + ψ ) , − ψ ˙ sin ( θ + ψ ) + θ ( θ ˙ + ψ ˙ ) cos ( θ + ψ ) ) ∣ v ⃗ ∣ 2 = R 2 ( ψ ˙ 2 + θ 2 ( θ ˙ + ψ ˙ ) 2 ) L = 1 2 I ψ ˙ 2 + 1 2 m ∣ v ⃗ ∣ 2 = 1 2 ( I + m R 2 ) ψ ˙ 2 + 1 2 m R 2 θ 2 ( θ ˙ + ψ ˙ ) 2 \begin{aligned}
\overrightarrow{\mathrm{OC}}
&= R \left( \sin(\theta+\psi), \cos(\theta+\psi) \right)
\\
\overrightarrow{\mathrm{CP}}
&= R \theta \left(-\cos(\theta+\psi), \sin(\theta+\psi) \right)
\\
\overrightarrow{\mathrm{OP}}
&= \overrightarrow{\mathrm{OC}} + \overrightarrow{\mathrm{CP}}
\\
&= R \left( \sin(\theta+\psi) - \theta \cos(\theta+\psi),
\cos (\theta + \psi) + \theta \sin (\theta + \psi) \right)
\\
\vec{v}
&= R \left(
\dot{\psi} \cos(\theta+\psi) + \theta (\dot{\theta}+\dot{\psi}) \sin(\theta+\psi),
-\dot{\psi} \sin(\theta+\psi) + \theta (\dot{\theta}+\dot{\psi}) \cos(\theta+\psi)
\right)
\\
\left| \vec{v} \right|^2
&= R^2 \left( \dot{\psi}^2 + \theta^2 (\dot{\theta}+\dot{\psi})^2 \right)
\\
\mathcal{L}
&= \frac{1}{2} I \dot{\psi}^2 + \frac{1}{2} m \left| \vec{v} \right|^2
\\
&= \frac{1}{2} \left( I + mR^2 \right) \dot{\psi}^2
+ \frac{1}{2} mR^2 \theta^2 \left( \dot{\theta} + \dot{\psi} \right)^2
\end{aligned} OC CP OP v ∣ v ∣ 2 L = R ( sin ( θ + ψ ) , cos ( θ + ψ ) ) = Rθ ( − cos ( θ + ψ ) , sin ( θ + ψ ) ) = OC + CP = R ( sin ( θ + ψ ) − θ cos ( θ + ψ ) , cos ( θ + ψ ) + θ sin ( θ + ψ ) ) = R ( ψ ˙ cos ( θ + ψ ) + θ ( θ ˙ + ψ ˙ ) sin ( θ + ψ ) , − ψ ˙ sin ( θ + ψ ) + θ ( θ ˙ + ψ ˙ ) cos ( θ + ψ ) ) = R 2 ( ψ ˙ 2 + θ 2 ( θ ˙ + ψ ˙ ) 2 ) = 2 1 I ψ ˙ 2 + 2 1 m ∣ v ∣ 2 = 2 1 ( I + m R 2 ) ψ ˙ 2 + 2 1 m R 2 θ 2 ( θ ˙ + ψ ˙ ) 2
問 9.
d d t ∂ L ∂ θ ˙ = d d t m R 2 θ 2 ( θ ˙ + ψ ˙ ) = m R 2 ( 2 θ θ ˙ ( θ ˙ + ψ ˙ ) + θ 2 ( θ ¨ + ψ ¨ ) ) ∂ L ∂ θ = m R 2 θ ( θ ˙ + ψ ˙ ) 2 d d t ∂ L ∂ ψ ˙ = d d t ( ( I + m R 2 ) ψ ˙ + m R 2 θ 2 ( θ ˙ + ψ ˙ ) ) = ( I + m R 2 ) ψ ¨ + 2 m R 2 θ θ ˙ ( θ ˙ + ψ ˙ ) + m R 2 θ 2 ( θ ¨ + ψ ¨ ) ∂ L ∂ ψ = 0 \begin{aligned}
\frac{\mathrm{d}}{\mathrm{d} t}
\frac{\partial \mathcal{L}}{\partial \dot{\theta}}
&= \frac{\mathrm{d}}{\mathrm{d} t}
mR^2 \theta^2 \left( \dot{\theta} + \dot{\psi} \right)
\\
&= mR^2 \left( 2 \theta \dot{\theta} (\dot{\theta}+\dot{\psi})
+ \theta^2 (\ddot{\theta}+\ddot{\psi}) \right)
\\
\frac{\partial \mathcal{L}}{\partial \theta}
&= mR^2 \theta \left( \dot{\theta} + \dot{\psi} \right)^2
\\
\frac{\mathrm{d}}{\mathrm{d} t}
\frac{\partial \mathcal{L}}{\partial \dot{\psi}}
&= \frac{\mathrm{d}}{\mathrm{d} t}
\left( \left( I+mR^2 \right) \dot{\psi}
+ mR^2 \theta^2 \left( \dot{\theta} + \dot{\psi} \right) \right)
\\
&= \left( I+mR^2 \right) \ddot{\psi}
+ 2mR^2 \theta \dot{\theta} \left( \dot{\theta} + \dot{\psi} \right)
+ mR^2 \theta^2 \left( \ddot{\theta} + \ddot{\psi} \right)
\\
\frac{\partial \mathcal{L}}{\partial \psi}
&= 0
\end{aligned} d t d ∂ θ ˙ ∂ L ∂ θ ∂ L d t d ∂ ψ ˙ ∂ L ∂ ψ ∂ L = d t d m R 2 θ 2 ( θ ˙ + ψ ˙ ) = m R 2 ( 2 θ θ ˙ ( θ ˙ + ψ ˙ ) + θ 2 ( θ ¨ + ψ ¨ ) ) = m R 2 θ ( θ ˙ + ψ ˙ ) 2 = d t d ( ( I + m R 2 ) ψ ˙ + m R 2 θ 2 ( θ ˙ + ψ ˙ ) ) = ( I + m R 2 ) ψ ¨ + 2 m R 2 θ θ ˙ ( θ ˙ + ψ ˙ ) + m R 2 θ 2 ( θ ¨ + ψ ¨ ) = 0
なので、 θ , ψ \theta, \psi θ , ψ に関するオイラー-ラグランジュ方程式はそれぞれ次のようになる:
2 θ θ ˙ ( θ ˙ + ψ ˙ ) + θ 2 ( θ ¨ + ψ ¨ ) = θ ( θ ˙ + ψ ˙ ) 2 ( I + m R 2 ) ψ ¨ + 2 m R 2 θ θ ˙ ( θ ˙ + ψ ˙ ) + m R 2 θ 2 ( θ ¨ + ψ ¨ ) = 0 \begin{aligned}
2 \theta \dot{\theta} \left( \dot{\theta}+\dot{\psi} \right)
+ \theta^2 \left( \ddot{\theta}+\ddot{\psi} \right)
&= \theta \left( \dot{\theta} + \dot{\psi} \right)^2
\\
\left( I+mR^2 \right) \ddot{\psi}
+ 2mR^2 \theta \dot{\theta} \left( \dot{\theta} + \dot{\psi} \right)
+ mR^2 \theta^2 \left( \ddot{\theta} + \ddot{\psi} \right)
&= 0
\end{aligned} 2 θ θ ˙ ( θ ˙ + ψ ˙ ) + θ 2 ( θ ¨ + ψ ¨ ) ( I + m R 2 ) ψ ¨ + 2 m R 2 θ θ ˙ ( θ ˙ + ψ ˙ ) + m R 2 θ 2 ( θ ¨ + ψ ¨ ) = θ ( θ ˙ + ψ ˙ ) 2 = 0
問 10.
問 9 からわかるように、
d d t ( ( I + m R 2 ) ψ ˙ + m R 2 θ 2 ( θ ˙ + ψ ˙ ) ) = 0 \begin{aligned}
\frac{\mathrm{d}}{\mathrm{d} t}
\left( \left( I+mR^2 \right) \dot{\psi}
+ mR^2 \theta^2 \left( \dot{\theta} + \dot{\psi} \right) \right)
= 0
\end{aligned} d t d ( ( I + m R 2 ) ψ ˙ + m R 2 θ 2 ( θ ˙ + ψ ˙ ) ) = 0
つまり、
( I + m R 2 ) ψ ˙ + m R 2 θ 2 ( θ ˙ + ψ ˙ ) \begin{aligned}
\left( I+mR^2 \right) \dot{\psi}
+ mR^2 \theta^2 \left( \dot{\theta} + \dot{\psi} \right)
\end{aligned} ( I + m R 2 ) ψ ˙ + m R 2 θ 2 ( θ ˙ + ψ ˙ )
は保存量である。
問 11.
t = 0 t=0 t = 0 において θ = 0 , ψ ˙ = 0 \theta = 0, \ \dot{\psi}=0 θ = 0 , ψ ˙ = 0 であるから、
問 10 の保存量は t ≥ 0 t \geq 0 t ≥ 0 において 0 0 0 である:
( I + m R 2 ) ψ ˙ + m R 2 θ 2 ( θ ˙ + ψ ˙ ) = 0 ∴ ψ ˙ = − m R 2 θ 2 θ ˙ I + m R 2 + m R 2 θ 2 \begin{aligned}
\left( I+mR^2 \right) \dot{\psi}
+ mR^2 \theta^2 \left( \dot{\theta} + \dot{\psi} \right)
= 0
\\
\therefore \ \
\dot{\psi}
= - \frac{mR^2 \theta^2 \dot{\theta}}
{I + mR^2 + mR^2 \theta^2}
\end{aligned} ( I + m R 2 ) ψ ˙ + m R 2 θ 2 ( θ ˙ + ψ ˙ ) = 0 ∴ ψ ˙ = − I + m R 2 + m R 2 θ 2 m R 2 θ 2 θ ˙
t > 0 t \gt 0 t > 0 において θ > 0 , θ ˙ > 0 \theta \gt 0, \ \dot{\theta} \gt 0 θ > 0 , θ ˙ > 0 であるから
ψ ˙ < 0 \dot{\psi} \lt 0 ψ ˙ < 0 であり、時計回りを正とした ψ \psi ψ が減少するため、円柱は反時計回りに回転することがわかる。
問 12.
t = 0 t=0 t = 0 での質点の運動エネルギー ( 1 / 2 ) m v 0 2 (1/2)mv_0^2 ( 1/2 ) m v 0 2
の一部は円柱の運動エネルギーとなるので、質点の運動エネルギーは保存しない。
Reference