前期 大問〔3〕

解法の指針

〔3〕は前半(問1)が等速円運動の位置・速度・加速度をベクトル成分で追う問題、後半(問2)がばねでつながれた多数の物体の単振動(1 次元格子振動)で、隣どうしの位相差から運動方程式を立て $\omega$ と $\lambda$ の関係(分散関係)を導く問題です。円運動を直線に投影すると単振動になる、という一本の糸で前後がつながっています。

Point

この設問での円運動は時計回りで、$t=0$ に $\vec{r}(0)=(0,\,r)$(=$y$ 軸の正の位置)から出発する点に注意。角 $\omega t$ は $y$ 軸から時計回り($+x$ 側)に測る。したがって $x$ 成分が $\sin$、$y$ 成分が $\cos$ になります(反時計回り・$x$ 軸出発の教科書公式とは $\sin$/$\cos$ が入れ替わる)。

使う道具

この大問で問われている 15 個の空欄(ア〜ソ) 答えはすべて $r,\ \omega,\ t,\ d,\ \lambda,\ m,\ k$ のみで表す(間隔は $d$、$x_n=nd$)。

問1 等速円運動の位置・速度・加速度(ア〜コ)

直感的理解
半径 $r$ の円周上を時計回りに回る点。出発点は真上 $(0,r)$。「位置 → 速度 → 加速度」は、いずれも大きさは一定($r,\ r\omega,\ r\omega^2$)で、向きだけが同じ角速度 $\omega$ で時計回りに回り続ける 3 本の矢印です。速度は位置より $90^\circ$、加速度は速度よりさらに $90^\circ$ 進んだ(=位置と逆向き・中心向きの)ベクトルになります。

位置ベクトル(ア・イ)

立式:時計回りで、時刻 0 に $\vec{r}(0)=(0,\,r)$。角 $\omega t$ を $y$ 軸から時計回り($+x$ 側)に測ると、$x$ 成分は $\sin$、$y$ 成分は $\cos$ で書ける:

$$x(t)=r\sin(\omega t),\qquad y(t)=r\cos(\omega t)$$

確認:$t=0$ で $(x,y)=(0,\,r)$。少し時間がたつと $x\gt 0$ に動き出す=時計回り。$x^2+y^2=r^2$ で半径 $r$ の円上にある。

初速度(ウ・エ)

立式:$\Delta t$ が十分小さいとき $\vec{r}(\Delta t)-\vec{r}(0)=\vec{v}(0)\,\Delta t$。$\sin\theta\fallingdotseq\theta,\ \cos\theta\fallingdotseq 1$($\theta=\omega\Delta t$ が微小)を使うと:

$$\vec{r}(\Delta t)-\vec{r}(0)=\big(r\sin(\omega\Delta t)-0,\ r\cos(\omega\Delta t)-r\big)\fallingdotseq\big(r\omega\Delta t,\ 0\big)$$

両辺を $\Delta t$ で割ると初速度が得られる:

$$v_x(0)=r\omega,\qquad v_y(0)=0$$

初速度は $+x$ 方向・大きさ $r\omega$。真上の点で時計回りなら、確かに右向きに動き出す。

速度(オ・カ)

立式:$\vec{v}(0)$ と $\vec{v}(t)$ のなす角は $\omega t$(速度ベクトルも時計回りに同じだけ回る)。$\vec{v}(0)=(r\omega,\,0)$ を時計回りに $\omega t$ 回すと:

$$v_x(t)=r\omega\cos(\omega t),\qquad v_y(t)=-r\omega\sin(\omega t)$$

検算(微小変化から):$x(t)=r\sin\omega t$ の $\Delta t$ 変化は $r\omega\cos(\omega t)\,\Delta t$、$y(t)=r\cos\omega t$ の変化は $-r\omega\sin(\omega t)\,\Delta t$。上式と一致する。大きさは $\sqrt{v_x^2+v_y^2}=r\omega$(一定)。

初期加速度と加速度(キ・ク・ケ・コ)

立式:同様に $\vec{v}(\Delta t)-\vec{v}(0)=\vec{a}(0)\,\Delta t$。$\vec{v}(t)$ を微小時間で展開すると:

$$\vec{v}(\Delta t)-\vec{v}(0)=\big(r\omega\cos(\omega\Delta t)-r\omega,\ -r\omega\sin(\omega\Delta t)\big)\fallingdotseq\big(0,\ -r\omega^2\Delta t\big)$$

$\Delta t$ で割って初期加速度:

$$a_x(0)=0,\qquad a_y(0)=-r\omega^2$$

加速度(ケ・コ):$\vec{a}(0)=(0,\,-r\omega^2)$($-y$ 方向=中心向き)を時計回りに $\omega t$ 回すと:

$$a_x(t)=-r\omega^2\sin(\omega t),\qquad a_y(t)=-r\omega^2\cos(\omega t)$$

これは $\vec{a}(t)=-\omega^2\,\vec{r}(t)$ に等しく、常に中心を向く向心加速度(大きさ $r\omega^2$)。

答え:
ア・イ:$x(t)=r\sin(\omega t),\ \ y(t)=r\cos(\omega t)$
ウ・エ:$v_x(0)=r\omega,\ \ v_y(0)=0$
オ・カ:$v_x(t)=r\omega\cos(\omega t),\ \ v_y(t)=-r\omega\sin(\omega t)$
キ・ク:$a_x(0)=0,\ \ a_y(0)=-r\omega^2$
ケ・コ:$a_x(t)=-r\omega^2\sin(\omega t),\ \ a_y(t)=-r\omega^2\cos(\omega t)$
別解:微分で一気に求める

$x=r\sin\omega t,\ y=r\cos\omega t$ を $t$ で微分すると速度、もう 1 回で加速度:

$$\vec{v}=\frac{d\vec{r}}{dt}=\big(r\omega\cos\omega t,\ -r\omega\sin\omega t\big)$$ $$\vec{a}=\frac{d\vec{v}}{dt}=\big(-r\omega^2\sin\omega t,\ -r\omega^2\cos\omega t\big)=-\omega^2\vec{r}$$

本問は微分を使わず「微小変化 = 速度$\times\Delta t$」で誘導しているが、結果は同じ。$\vec{a}=-\omega^2\vec{r}$ が向心加速度の本質。

Point

回転の向き(時計回り)と出発点 $(0,r)$ を取り違えると $\sin$/$\cos$ や符号を全部間違える。「$t=0$ でどこにいて、次の瞬間どちらへ動くか」を必ず 1 点だけ確かめて向きを固定するのがコツ。

問2 格子上の単振動:隣どうしの回転角(サ・シ)

直感的理解
間隔 $d$ で並ぶ物体たちの変位が、波長 $\lambda$ の正弦波の上に乗って $+x$ 向きに進む。波が右へ進むので、左の物体($P-1$)は先に同じ振動をし(位相が進む)、右の物体($P+1$)は遅れる。位置が $d$ ずれるごとに位相は $\dfrac{2\pi}{\lambda}d$ だけずれる。基準の $P$ の回転角が $\omega t$ なら、隣は $\omega t\pm\dfrac{2\pi d}{\lambda}$。

立式:物体は間隔 $d$、位置 $x_n=nd$。変位が波長 $\lambda$ の正弦波に乗り $+x$ 向きに進むので、$P$ 番目を基準(回転角 $\omega t$、$u_P=r\sin(\omega t)$)にとると、$n$ 番目の変位は

$$u_n(t)=r\sin\!\left(\omega t-\frac{2\pi}{\lambda}(x_n-x_P)\right)=r\sin\!\left(\omega t-\frac{2\pi d}{\lambda}(n-P)\right)$$

波が $+x$ に進むとき、位置が大きいほど「同じ位相に達するのが遅れる」ので、位相(回転角)に $-\dfrac{2\pi}{\lambda}(x_n-x_P)$ が付く。

代入($P-1$ 番目):$n=P-1$ なら $x_{P-1}-x_P=-d$。

$$u_{P-1}(t)=r\sin\!\left(\omega t-\frac{2\pi d}{\lambda}\cdot(-1)\right)=r\sin\!\left(\omega t+\frac{2\pi d}{\lambda}\right)$$

左の物体は位相が $\dfrac{2\pi d}{\lambda}$ だけ進む。よって サ $=\omega t+\dfrac{2\pi d}{\lambda}$。

代入($P+1$ 番目):$n=P+1$ なら $x_{P+1}-x_P=+d$。

$$u_{P+1}(t)=r\sin\!\left(\omega t-\frac{2\pi d}{\lambda}\right)$$

右の物体は位相が $\dfrac{2\pi d}{\lambda}$ だけ遅れる。よって シ $=\omega t-\dfrac{2\pi d}{\lambda}$。

答え:
サ($P-1$ 番目の回転角):$\displaystyle \omega t+\frac{2\pi d}{\lambda}$
シ($P+1$ 番目の回転角):$\displaystyle \omega t-\frac{2\pi d}{\lambda}$
補足:位相差の見積もり

1 波長 $\lambda$ で位相が $2\pi$ 進む。間隔 $d$ に対応する位相は比例配分で $\dfrac{2\pi}{\lambda}\times d=\dfrac{2\pi d}{\lambda}$。

例:$\lambda=4d$ なら位相差 $\dfrac{\pi}{2}$($90^\circ$)。$\lambda=2d$ なら位相差 $\pi$($180^\circ$)で隣どうしが逆位相。$\lambda\gg d$ なら位相差はほぼ 0 で、隣どうしほとんど同じ動き(=なめらかな長波長の波)。

Point

「$+x$ 向きに進む波」では、後方($x$ が小さい側)が位相が進む。$u_n=r\sin(\omega t-\frac{2\pi}{\lambda}x)$ の $x$ の前の符号がマイナスであることと結びつけて覚えると、サ・シ の符号を間違えない。

問2 加速度・左右ばねの合力・分散関係(ス・セ・ソ)

直感的理解
$P$ 番目の物体を単独で見れば、変位 $u_P=r\sin\omega t$ の単振動なので加速度は $-\omega^2 u_P$。一方、$P$ を動かす力は左右のばねが引っぱる(押す)合力で、これは隣との変位の差で決まる。両者を運動方程式で結ぶと、$\omega$ と $\lambda$ の関係(分散関係)が出てくる。和積・半角公式で $\sin^2\!\dfrac{\pi d}{\lambda}$ の形にまとめるのがゴール。

ス:$P$ 番目の加速度

立式・計算:変位 $u_P=r\sin(\omega t)$ を時間で 2 回微分(単振動なので加速度は $-\omega^2\times$変位):

$$a_P=\frac{d^2u_P}{dt^2}=-r\omega^2\sin(\omega t)=-\omega^2 u_P$$

よって ス $=-r\omega^2\sin(\omega t)$。

セ:左右ばねからの合力の係数

立式:ばね定数 $k$。ばねの力は「隣との変位の差」で決まる。左のばねは $k(u_{P-1}-u_P)$、右のばねは $k(u_{P+1}-u_P)$。合力は

$$F=k(u_{P-1}-u_P)+k(u_{P+1}-u_P)=k\big(u_{P-1}+u_{P+1}-2u_P\big)$$

代入(サ・シ を利用):

$$u_{P-1}+u_{P+1}=r\sin\!\left(\omega t+\frac{2\pi d}{\lambda}\right)+r\sin\!\left(\omega t-\frac{2\pi d}{\lambda}\right)$$

和積公式 $\sin(A{+}B)+\sin(A{-}B)=2\sin A\cos B$($A=\omega t,\ B=\frac{2\pi d}{\lambda}$)より:

$$u_{P-1}+u_{P+1}=2r\sin(\omega t)\cos\!\left(\frac{2\pi d}{\lambda}\right)=2u_P\cos\!\left(\frac{2\pi d}{\lambda}\right)$$

これを合力に戻すと:

$$F=k\left[\,2u_P\cos\!\left(\frac{2\pi d}{\lambda}\right)-2u_P\right]=-2k\left[1-\cos\!\left(\frac{2\pi d}{\lambda}\right)\right]u_P$$

半角公式 $1-\cos\theta=2\sin^2\dfrac{\theta}{2}$($\theta=\frac{2\pi d}{\lambda}$)を使うと:

$$F=-4k\sin^2\!\left(\frac{\pi d}{\lambda}\right)u_P=\left[-4k\sin^2\!\left(\frac{\pi d}{\lambda}\right)\right]\times r\sin(\omega t)$$

問題文は「合力 = セ $\times\,r\sin(\omega t)$」だから、セ $=-4k\sin^2\!\left(\dfrac{\pi d}{\lambda}\right)$。

ソ:$\omega^2$ と $\lambda$ の関係(分散関係)

立式:運動方程式 $m a_P=F$ に ス と セ を代入する:

$$m\big(-r\omega^2\sin\omega t\big)=\left[-4k\sin^2\!\left(\frac{\pi d}{\lambda}\right)\right]r\sin\omega t$$

両辺の $-r\sin(\omega t)$ を約分して整理:

$$m\omega^2=4k\sin^2\!\left(\frac{\pi d}{\lambda}\right)\quad\Longrightarrow\quad \omega^2=\frac{4k}{m}\sin^2\!\left(\frac{\pi d}{\lambda}\right)$$
答え:
ス($P$ の加速度):$\displaystyle -r\omega^2\sin(\omega t)$
セ(合力の係数):$\displaystyle -4k\sin^2\!\left(\frac{\pi d}{\lambda}\right)$
ソ(分散関係):$\displaystyle \omega^2=\frac{4k}{m}\sin^2\!\left(\frac{\pi d}{\lambda}\right)$
別解:合力を「差分の差分」でとらえる

$F=k(u_{P+1}-u_P)-k(u_P-u_{P-1})$ と書くと、これは「右のばねの伸び」から「左のばねの伸び」を引いた形(隣接差の差)。$u_n=r\sin(\omega t-\frac{2\pi d}{\lambda}(n-P))$ を代入し、$\sin$ の差を積和公式で整理しても同じ結果に至る:

$$F=-2k\Big[1-\cos\tfrac{2\pi d}{\lambda}\Big]u_P=-4k\sin^2\!\Big(\tfrac{\pi d}{\lambda}\Big)u_P$$

物理的には「隣との差が大きいほど(=波長が短いほど)復元力が強く、$\omega$ が大きい」ことを表している。

補足:長波長極限と最大振動数

長波長 $\lambda\gg d$:$\dfrac{\pi d}{\lambda}$ が微小なので $\sin\dfrac{\pi d}{\lambda}\fallingdotseq\dfrac{\pi d}{\lambda}$。

$$\omega^2\fallingdotseq\frac{4k}{m}\left(\frac{\pi d}{\lambda}\right)^2\ \Longrightarrow\ \omega\fallingdotseq\frac{2\pi d}{\lambda}\sqrt{\frac{k}{m}}$$

$\omega\propto\dfrac{1}{\lambda}$ となり、波の速さ $v=\dfrac{\lambda}{2\pi}\omega=d\sqrt{\dfrac{k}{m}}$(一定)。これはふつうの音波(非分散)にあたる。

最短波長 $\lambda=2d$:$\sin^2\dfrac{\pi d}{\lambda}=\sin^2\dfrac{\pi}{2}=1$ で $\omega^2=\dfrac{4k}{m}$、$\omega_{\max}=2\sqrt{\dfrac{k}{m}}$。隣どうしが逆位相で振動する状態で、これ以上短い波は格子上に存在できない。

Point

ばねの力は「その物体の変位」ではなく隣との変位の差で決まる。ここを $-ku_P$ とだけ書くと誤り。和積・半角公式で $\sin^2\!\dfrac{\pi d}{\lambda}$ にまとめる流れは、格子振動(フォノン)の分散関係そのもの。

まとめ:円運動 → 単振動 → 格子振動

直感的理解
等速円運動を軸に投影すると単振動になる。この単振動を間隔 $d$ で多数並べ、ばねでつなぐと、位相のずれた振動が正弦波として伝わる(格子振動)。問1 の $\vec{a}=-\omega^2\vec{r}$ と、問2 の $a_P=-\omega^2u_P$ は同じ「振動の方程式」で、円運動と格子波が一本の理屈でつながっている。

要点の整理

数値でチェック(分散関係)

例として $k=1.0\times10^{2}\ \text{N/m}$、$m=1.0\times10^{-2}\ \text{kg}$ の格子を考える。最短波長 $\lambda=2d$ のとき、$\sin^2\!\dfrac{\pi d}{\lambda}=\sin^2\dfrac{\pi}{2}=1$ なので最大角振動数は:

$$\omega_{\max}=2\sqrt{\frac{k}{m}}=2\sqrt{\frac{1.0\times10^{2}}{1.0\times10^{-2}}}=2\times100=2.0\times10^{2}\ \text{rad/s}$$

すなわち最大角振動数は約 200 rad/s(周期にすると約 0.031 s)。一方 $\lambda=4d$ なら $\sin^2\!\dfrac{\pi d}{\lambda}=\sin^2\dfrac{\pi}{4}=\dfrac{1}{2}$ なので:

$$\omega=\sqrt{\frac{4k}{m}\cdot\frac{1}{2}}=\sqrt{\frac{2k}{m}}=\sqrt{\frac{2\times1.0\times10^{2}}{1.0\times10^{-2}}}\fallingdotseq 1.4\times10^{2}\ \text{rad/s}$$

波長を短くする(隣との差を大きくする)ほど $\omega$ が大きくなることが数値でも確認できる。

Point

「円運動の投影=単振動」という見方を軸に、$\vec{a}=-\omega^2\vec{r}$ から $\omega^2=\dfrac{4k}{m}\sin^2\!\dfrac{\pi d}{\lambda}$ まで一直線。符号(時計回り・$+x$ 伝搬)「ばねの力=変位の差」の 2 点さえ外さなければ完答できる大問。