跳到主要内容

京都大学 理学研究科 物理学・宇宙物理学専攻 2019年度 問題III-1

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

自然長 aa、ばね定数 kk のばねで質点を一直線につなぐ。まず図1の両端固定、全長 4a4a の鎖の3質点を考え、左から位置を a+x1,2a+x2,3a+x3a+x_1,2a+x_2,3a+x_3、質量を m,M,mm,M,m とする。ω0=2k/m\omega_0=\sqrt{2k/m}、β=m/M\beta=m/M とおく。

原題 図1:両端固定の3質点ばね系

  1. 両外側の質点を x1=x3=0x_1=x_3=0 に固定したとき、初期条件 x2(0)=X,x˙2(0)=Vx_2(0)=X,\dot x_2(0)=V に対する運動を求めよ。
  2. 全質点を自由に動かした運動方程式を求め、xj=Ajcos⁡ωtx_j=A_j\cos\omega t として、小さい順の基準角振動数 ωs,ωc,ωl\omega_s,\omega_c,\omega_l の二乗を求めよ。
  3. β=1\beta=1 で各モードの A2,A3A_2,A_3 を A1A_1 で表せ。
  4. 0≤β≤20\leq\beta\leq2 で ωs2/ω02,ωc2/ω02,ωl2/ω02\omega_s^2/\omega_0^2,\omega_c^2/\omega_0^2,\omega_l^2/\omega_0^2 を図示し、β=0,1\beta=0,1 の値を示せ。

次に図2のように十分多数の質点を同じばねで一直線につなぎ、中央の n=0n=0 のみ質量 MM、残りは mm とする。平衡位置は nana、変位を xnx_n とする。

原題 図2:中央の質量が異なるばね鎖

  1. 中央を角振動数 ω\omega で強制振動させる。xn=Re⁡[Bneiωt]x_n=\operatorname{Re}[B_ne^{i\omega t}]、Bn=B0bnB_n=B_0b^n を質量 mm の運動方程式へ代入し、bb と ∣b∣≠1|b|\neq1 となる条件を求めよ。∣b∣=1|b|=1 は進行波、∣b∣<1|b|<1 は局在波に対応する。
  2. 強制力を止める。M<mM<m では局在した基準振動が存在する。左右で ∣n∣|n| とともに減衰する解を選び、x−1,x+1x_{-1},x_{+1} を x0,bx_0,b で表せ。中央の運動方程式からその角振動数を求め、m→∞m\to\infty の極限も求めよ。

题目描述​

如上方图1,三质点以相同弹簧连接两固定端,弹簧自然长 aa、劲度 kk,质量为 m,M,mm,M,m,平衡位置为 a,2a,3aa,2a,3a。定义 ω0=2k/m,β=m/M\omega_0=\sqrt{2k/m},\beta=m/M。

  1. 固定两侧质点,求中央质点给定初位移和初速度下的运动。
  2. 放开三质点,求方程及三个本征频率。
  3. 等质量时求三个本征振型。
  4. 绘制频率平方随质量比的曲线。
  5. 对上方图2中仅中央质量不同的无限长弹簧链,以 Bn=B0bnB_n=B_0b^n 求传播因子及非传播频率条件。
  6. 放开中央质点,求 M<mM<m 时左右衰减的局域模频率及 m→∞m\to\infty 极限。

Kai​

(1)–(3)​

(1) Mx¨2=−2kx2M\ddot x_2=-2kx_2。Ω=2k/M\Omega=\sqrt{2k/M} とおけば

x2=Xcos⁡Ωt+VΩsin⁡Ωt.\boxed{x_2=X\cos\Omega t+\frac V\Omega\sin\Omega t.}

(2) 運動方程式は

mx¨1=k(x2−2x1),Mx¨2=k(x1−2x2+x3),mx¨3=k(x2−2x3).m\ddot x_1=k(x_2-2x_1),\quad M\ddot x_2=k(x_1-2x_2+x_3),\quad m\ddot x_3=k(x_2-2x_3).

z=ω2/ω02z=\omega^2/\omega_0^2 とすると固有方程式は (z−1)[z2−(1+β)z+β/2]=0(z-1)[z^2-(1+\beta)z+\beta/2]=0。したがって

ωs2ω02=1+β−1+β22,ωc2ω02=1,ωl2ω02=1+β+1+β22.\boxed{\frac{\omega_s^2}{\omega_0^2}=\frac{1+\beta-\sqrt{1+\beta^2}}2,\quad \frac{\omega_c^2}{\omega_0^2}=1,\quad \frac{\omega_l^2}{\omega_0^2}=\frac{1+\beta+\sqrt{1+\beta^2}}2.}

(3) β=1\beta=1 では順に

(A1,A2,A3)∝(1,2,1), (1,0,−1), (1,−2,1).\boxed{(A_1,A_2,A_3)\propto(1,\sqrt2,1),\ (1,0,-1),\ (1,-\sqrt2,1).}

(4)​

質量比と三つの基準振動数

β=0\beta=0 で値は (0,1,1)(0,1,1)、β=1\beta=1 で (1−1/2,1,1+1/2)(1-1/\sqrt2,1,1+1/\sqrt2)。上下の枝は単調増加し、中間の枝は一定である。β=0\beta=0 は M→∞M\to\infty の極限を意味する。

(5)​

n≠0n\neq0 では mx¨n=k(xn−1−2xn+xn+1)m\ddot x_n=k(x_{n-1}-2x_n+x_{n+1})。z=ω2/ω02z=\omega^2/\omega_0^2 とすると

b+b−1=2(1−z),b=1−z±(1−z)2−1.b+b^{-1}=2(1-z),\qquad \boxed{b=1-z\pm\sqrt{(1-z)^2-1}.}

実周波数では 0≤z≤20\leq z\leq2 で ∣b∣=1|b|=1、z>2z>2 で互いに逆数の負の実根を持つ。従って

∣b∣≠1 ⟺ ω2>2ω02.\boxed{|b|\neq1\ \Longleftrightarrow\ \omega^2>2\omega_0^2.}

(6)​

z>2z>2 では b=1−z+(1−z)2−1b=1-z+\sqrt{(1-z)^2-1} が −1<b<0-1<b<0 を満たす。左右に減衰させるには Bn=B0b∣n∣B_n=B_0b^{|n|} とすればよく、x−1=x+1=bx0x_{-1}=x_{+1}=bx_0。中央の式から z=β(1−b)z=\beta(1-b)。これと (5) を連立すると

b=−12β−1,ω2=2β22β−1ω02.\boxed{b=-\frac1{2\beta-1},\qquad \omega^2=\frac{2\beta^2}{2\beta-1}\omega_0^2.}

局在条件 ∣b∣<1|b|<1 は β>1\beta>1、すなわち M<mM<m と一致する。M,kM,k を固定して m→∞m\to\infty とすると、周囲は動かなくなり

ω⟶2k/M.\boxed{\omega\longrightarrow\sqrt{2k/M}.}

Reference​