跳到主要内容

東北大学 理学研究科 物理学専攻 2018年8月実施 問題2(力学)

Author​

Miyake

Description​

下図に示すばねと振子からなる装置を考えよう.O を原点として水平方向に xx 軸,鉛直上向きに yy 軸をとる.ばね(ばね定数 kk)の一端は壁に固定され,反対側の端に質量 2m2m の質点 P が取り付けられている.ばねと質点 P は,xx 軸に沿ってたわむことなく,摩擦なしに運動する.さらに,長さ ll の剛体棒と質点 Q(質量 mm)からなる振子が,質点 P に取り付けられ,質点 P を支点として xyxy 平面内で自由に振れるようになっている.ばねが自然長をとるとき,質点 P は原点 O にあるとし,剛体棒が yy 軸となす角を図のように θ\theta とする.また,重力加速度を gg とし,ばねと剛体棒の質量は無視できるものとして以下の問いに答えよ.

ばねと振子からなる装置

  1. 質点 P が座標 (x,0)(x,0) にあり,剛体棒の角度が θ\theta であるとき,質点 Q の座標 (X,Y)(X,Y),及びその速度 (X˙,Y˙)(\dot X,\dot Y) を,l,x,x˙,θ,θ˙l,x,\dot x,\theta,\dot\theta のうち必要なものを用いて表せ.

  2. この装置のラグランジアン LL を求めよ.

  3. x,x˙,θ,θ˙x,\dot x,\theta,\dot\theta が微小であるとして,問2)で求めた LL を,x,x˙,θ,θ˙x,\dot x,\theta,\dot\theta の2次まで展開せよ.

以下の問いでは,ばね定数が k=2mglk=\dfrac{2mg}{l} の場合を考える.また,s=lθ, Ω=gls=l\theta,\ \Omega=\sqrt{\dfrac gl} とする.

  1. 問3)の LL を用いて,xx と ss が満たす2個の運動方程式を,それぞれ s,s¨,x,x¨,Ωs,\ddot s,x,\ddot x,\Omega のうち必要なものを用いて表せ.

  2. 問4)の運動方程式の解を,x=Asin⁡ωt, s=Bsin⁡ωtx=A\sin\omega t,\ s=B\sin\omega t の形に仮定する.ここで,A,B,ωA,B,\omega は定数とする.自明でない解((A,B)=(0,0)(A,B)=(0,0) 以外の解)が存在するように ω\omega の値 ω1,ω2\omega_1,\omega_2 を求め,それぞれ Ω\Omega を用いて表せ.ただし,ω1>ω2>0\omega_1>\omega_2>0 とする.

  3. 問5)で求めた ω1,ω2\omega_1,\omega_2 それぞれに対して,AA と BB の比を求め,θ\theta が最大となる時刻での質点 P と Q の位置を解答用紙の図中に書き入れよ.

题目描述​

考虑上图所示由弹簧和摆组成的装置。以 O 为原点,水平方向为 xx 轴,竖直向上为 yy 轴。弹簧(劲度系数 kk)的一端固定在墙上,另一端连接质量为 2m2m 的质点 P。弹簧和质点 P 沿 xx 轴无弯曲、无摩擦地运动。此外,由长度为 ll 的刚杆和质点 Q(质量 mm)组成的摆连接于质点 P,以 P 为支点在 xyxy 平面内自由摆动。弹簧处于自然长度时 P 位于原点 O;刚杆与 yy 轴的夹角按图定义为 θ\theta。重力加速度为 gg,忽略弹簧和刚杆的质量,回答下列问题。

1) 当 P 的坐标为 (x,0)(x,0)、刚杆角度为 θ\theta 时,用 l,x,x˙,θ,θ˙l,x,\dot x,\theta,\dot\theta 中必要的量表示 Q 的坐标 (X,Y)(X,Y) 及速度 (X˙,Y˙)(\dot X,\dot Y)。

2) 求该装置的拉格朗日量 LL。

3) 设 x,x˙,θ,θ˙x,\dot x,\theta,\dot\theta 都是微小量,将问2)所得 LL 展开至这些量的二阶。

以下考虑 k=2mglk=\dfrac{2mg}{l},并令 s=lθ, Ω=gls=l\theta,\ \Omega=\sqrt{\dfrac gl}。

4) 利用问3)中的 LL,分别用 s,s¨,x,x¨,Ωs,\ddot s,x,\ddot x,\Omega 中必要的量表示 x,sx,s 满足的两个运动方程。

5) 假设问4)的运动方程的解具有 x=Asin⁡ωt, s=Bsin⁡ωtx=A\sin\omega t,\ s=B\sin\omega t 的形式,其中 A,B,ωA,B,\omega 为常数。求存在非平凡解(即 (A,B)≠(0,0)(A,B)\ne(0,0))的两个 ω\omega 值 ω1,ω2\omega_1,\omega_2,分别用 Ω\Omega 表示,且 ω1>ω2>0\omega_1>\omega_2>0。

6) 对问5)求出的 ω1,ω2\omega_1,\omega_2,分别求 AA 与 BB 的比值,并在答题纸的图中标出 θ\theta 最大时 P、Q 的位置。

Kai​

1)​

X=x+lsin⁡θY=−lcos⁡θX˙=x˙+lθ˙cos⁡θY˙=lθ˙sin⁡θ\begin{aligned} X &= x + l \sin \theta \\ Y &= - l \cos \theta \\ \dot{X} &= \dot{x} + l \dot{\theta} \cos \theta \\ \dot{Y} &= l \dot{\theta} \sin \theta \end{aligned}

2)​

ラグランジアンは、(運動エネルギー) - (位置エネルギー) であるから、

L=122mx˙2+12m(X˙2+Y˙2)−12kx2−mgY=12m(3x˙2+2lx˙θ˙cos⁡θ+l2θ˙2)−12kx2+mglcos⁡θ\begin{aligned} L &= \frac{1}{2} 2m \dot{x}^2 + \frac{1}{2} m \left( \dot{X}^2 + \dot{Y}^2 \right) - \frac{1}{2} k x^2 - mgY \\ &= \frac{1}{2} m \left( 3 \dot{x}^2 + 2 l \dot{x} \dot{\theta} \cos \theta + l^2 \dot{\theta}^2 \right) - \frac{1}{2} k x^2 + mgl \cos \theta \end{aligned}

である。

3)​

問2) で求めた LL は x,x˙,θ,θ˙x, \dot{x}, \theta, \dot{\theta} の2次までで、

L=12m(3x˙2+2lx˙θ˙+l2θ˙2)−12kx2+mgl(1−12θ2)\begin{aligned} L &= \frac{1}{2} m \left( 3 \dot{x}^2 + 2 l \dot{x} \dot{\theta} + l^2 \dot{\theta}^2 \right) - \frac{1}{2} k x^2 + mgl \left( 1 - \frac{1}{2} \theta^2 \right) \end{aligned}

である。

4)​

問3) で求めた LL について、次のように計算する:

ddt∂L∂x˙=ddt(3mx˙+mlθ˙)=3mx¨+mlθ¨=3mx¨+ms¨,∂L∂x=−kx=−2mΩ2x,ddt∂L∂θ˙=ddt(mlx˙+ml2θ˙)=mlx¨+ml2θ¨=mlx¨+mls¨,∂L∂θ=−mglθ=−mgs.\begin{aligned} \frac{d}{dt} \frac{\partial L}{\partial \dot{x}} &= \frac{d}{dt} \left( 3m \dot{x} + ml \dot{\theta} \right) = 3m \ddot{x} + ml \ddot{\theta} = 3m \ddot{x} + m \ddot{s} ,\\ \frac{\partial L}{\partial x} &= -kx = -2m \Omega^2 x ,\\ \frac{d}{dt} \frac{\partial L}{\partial \dot{\theta}} &= \frac{d}{dt} \left( ml \dot{x} + ml^2 \dot{\theta} \right) = ml \ddot{x} + ml^2 \ddot{\theta} = ml \ddot{x} + ml \ddot{s} ,\\ \frac{\partial L}{\partial \theta} &= -mgl \theta = -mgs . \end{aligned}

よって、運動方程式は次のようになる:

{3mx¨+ms¨=−2mΩ2xmlx¨+mls¨=−mgs∴    {3x¨+s¨=−2Ω2xx¨+s¨=−Ω2s\begin{aligned} & \begin{cases} 3m \ddot{x} + m \ddot{s} = -2m \Omega^2 x \\ ml \ddot{x} + ml \ddot{s} = -mgs \end{cases} \\ \therefore \ \ \ \ & \begin{cases} 3 \ddot{x} + \ddot{s} = - 2 \Omega^2 x \\ \ddot{x} + \ddot{s} = - \Omega^2 s \end{cases} \end{aligned}

5)​

問4) で求めた運動方程式に x=Asin⁡ωt,s=Bsin⁡ωtx=A \sin \omega t, s=B \sin \omega t を代入すると、

{−3ω2Asin⁡ωt−ω2Bsin⁡ωt=−2Ω2Asin⁡ωt−ω2Asin⁡ωt−ω2Bsin⁡ωt=−Ω2Bsin⁡ωt∴    {(3ω2−2Ω2)A+ω2B=0ω2A+(ω2−Ω2)B=0\begin{align} & \begin{cases} -3 \omega^2 A \sin \omega t - \omega^2 B \sin \omega t = -2 \Omega^2 A \sin \omega t \\ - \omega^2 A \sin \omega t - \omega^2 B \sin \omega t = - \Omega^2 B \sin \omega t \end{cases} \nonumber \\ \therefore \ \ \ \ & \begin{cases} \left( 3 \omega^2 - 2 \Omega^2 \right) A + \omega^2 B = 0 \\ \omega^2 A + \left( \omega^2 - \Omega^2 \right) B = 0 \end{cases} \tag{1} \end{align}

となるから、自明でない解をもつための条件は、

0=∣3ω2−2Ω2ω2ω2ω2−Ω2∣=(ω+2Ω)(ω−2Ω)(2ω+Ω)(2ω−Ω)\begin{gathered} 0 = \begin{vmatrix} 3 \omega^2 - 2 \Omega^2 & \omega^2 \\ \omega^2 & \omega^2 - \Omega^2 \end{vmatrix} \\ = (\omega + \sqrt{2} \Omega) (\omega - \sqrt{2} \Omega) (\sqrt{2} \omega + \Omega) (\sqrt{2} \omega - \Omega) \end{gathered}

である。

ω1>ω2>0\omega_1 \gt \omega_2 \gt 0 であるから、

ω1=2Ω,  ω2=Ω2\begin{aligned} \omega_1 = \sqrt{2} \Omega , \ \ \omega_2 = \frac{\Omega}{\sqrt{2}} \end{aligned}

を得る。

6)​

(i) ω=ω1\omega = \omega_1 を式(11)に代入して、次を得る:

AB=−12.\begin{aligned} \frac{A}{B} = - \frac{1}{2} . \end{aligned}

(ii) ω=ω2\omega = \omega_2 を式(11)に代入して、次を得る:

AB=1.\begin{aligned} \frac{A}{B} = 1 . \end{aligned}

図の位置は微小振動の精度で次のように表される。s=b>0s=b>0 を最大変位とすると、θ=b/l\theta=b/l であり、

  • ω1=2Ω\omega_1=\sqrt2\Omega:P=(−b/2,0)P=(-b/2,0)、Q=(b/2,−l)Q=(b/2,-l)。支点と振り子の水平変位は逆向き。
  • ω2=Ω/2\omega_2=\Omega/\sqrt2:P=(b,0)P=(b,0)、Q=(2b,−l)Q=(2b,-l)。両者の水平変位は同じ向き。

ここで Y=−l+O(b2/l)Y=-l+O(b^2/l) と近似している。

Spring and pendulum: positions at maximum positive angle

Apparatus

High frequency: opposite motion

Low frequency: same direction

Mode sketches: s=lθ=b>0s=l\theta=b>0; coordinates use the small-angle approximation. Horizontal displacements are enlarged.

Reference​