解析力学から様々な理論分野を組み立てる 第2部 第5章 §2.5.6 ★★ 弦の素片に Newton の運動方程式を立てる導出と場の Euler–Lagrange 方程式からは同じ式が得られ、端の条件とエネルギーの流れも一致する
弦の素片の運動方程式
解析力学を使わない力学・振動波動の授業では、弦の波動方程式を次のように出す。張力 T s T_{\mathrm s} T s で張った弦の、x x x と x + Δ x x+\Delta x x + Δ x の間の小さな部分(素片)を考える。素片の両端で弦が水平から角 θ \theta θ だけ傾いていると、張力の縦の成分は T s sin θ ≈ T s tan θ = T s ∂ x ϕ T_{\mathrm s}\sin\theta\approx T_{\mathrm s}\tan\theta=T_{\mathrm s}\partial_x\phi T s sin θ ≈ T s tan θ = T s ∂ x ϕ である。素片には右端から + T s ∂ x ϕ ( x + Δ x ) +T_{\mathrm s}\partial_x\phi(x+\Delta x) + T s ∂ x ϕ ( x + Δ x ) 、左端から − T s ∂ x ϕ ( x ) -T_{\mathrm s}\partial_x\phi(x) − T s ∂ x ϕ ( x ) の縦の力がかかるので、Newton の運動方程式は
ρ Δ x ∂ t 2 ϕ = T s [ ∂ x ϕ ( x + Δ x ) − ∂ x ϕ ( x ) ] = T s ∂ x 2 ϕ Δ x + O ( Δ x 2 ) (2.5.29) \rho\,\Delta x\;\partial_t^2\phi=T_{\mathrm s}\bigl[\partial_x\phi(x+\Delta x)-\partial_x\phi(x)\bigr]=T_{\mathrm s}\,\partial_x^2\phi\,\Delta x+\mathcal{O}(\Delta x^2)\tag{2.5.29} ρ Δ x ∂ t 2 ϕ = T s [ ∂ x ϕ ( x + Δ x ) − ∂ x ϕ ( x ) ] = T s ∂ x 2 ϕ Δ x + O ( Δ x 2 ) ( 2.5.29 )
で、Δ x → 0 \Delta x\to0 Δ x → 0 で ρ ∂ t 2 ϕ = T s ∂ x 2 ϕ \rho\,\partial_t^2\phi=T_{\mathrm s}\partial_x^2\phi ρ ∂ t 2 ϕ = T s ∂ x 2 ϕ になる。これは式 (2.5.15) の Y Y Y を T s T_{\mathrm s} T s にしたものである(T s T_{\mathrm s} T s の添字 s は運動エネルギー・温度の T T T と分けるため)。
この導出を場の言葉で読み直す。張った弦の位置エネルギーは「伸びた長さ×張力」で、素片の長さは 1 + ( ∂ x ϕ ) 2 Δ x \sqrt{1+(\partial_x\phi)^2}\,\Delta x 1 + ( ∂ x ϕ ) 2 Δ x だから、位置エネルギーの密度は
T s ( 1 + ( ∂ x ϕ ) 2 − 1 ) = T s 2 ( ∂ x ϕ ) 2 + O ( ( ∂ x ϕ ) 4 ) T_{\mathrm s}\Bigl(\sqrt{1+(\partial_x\phi)^2}-1\Bigr)=\frac{T_{\mathrm s}}{2}(\partial_x\phi)^2+\mathcal{O}\bigl((\partial_x\phi)^4\bigr) T s ( 1 + ( ∂ x ϕ ) 2 − 1 ) = 2 T s ( ∂ x ϕ ) 2 + O ( ( ∂ x ϕ ) 4 )
である。式 (2.5.1) の − Y 2 ( ∂ x ϕ ) 2 -\frac Y2(\partial_x\phi)^2 − 2 Y ( ∂ x ϕ ) 2 は、この展開の最初の項を残したものになる。授業の導出の sin θ ≈ tan θ \sin\theta\approx\tan\theta sin θ ≈ tan θ は、L \mathcal L L のこの展開を2次で打ち切ることに当たる。
近似の裏付け
型 B(極限) | 小さいもの 傾き ∣ ∂ x ϕ ∣ → 0 \lvert \partial_x\phi\rvert \to0 ∣ ∂ x ϕ ∣ → 0 (変位の勾配が小さい) | 証明 条件付き(方程式の係数の誤差は O ( ∣ ∂ x ϕ ∣ 2 ) \mathcal{O}(\lvert \partial_x\phi\rvert ^2) O (∣ ∂ x ϕ ∣ 2 ) 。長い時間にわたる解の誤差の評価は要確認) | 行き先 近似の数理 第2章
並べて見る
授業での書き方 この章の書き方 出発点 素片にかかる張力の縦成分 L = ρ 2 ( ∂ t ϕ ) 2 − Y 2 ( ∂ x ϕ ) 2 \mathcal L=\frac\rho2(\partial_t\phi)^2-\frac Y2(\partial_x\phi)^2 L = 2 ρ ( ∂ t ϕ ) 2 − 2 Y ( ∂ x ϕ ) 2 運動方程式 素片の Newton の式 (2.5.29) 場の Euler–Lagrange 方程式 (2.5.12) 自由端 質量の無い輪に縦の力がかからない:T s ∂ x ϕ = 0 T_{\mathrm s}\partial_x\phi=0 T s ∂ x ϕ = 0 端の項 ∂ L ∂ ( ∂ x ϕ ) = 0 \frac{\partial \mathcal L}{\partial (\partial_x\phi)}=0 ∂ ( ∂ x ϕ ) ∂ L = 0 (式 (2.5.20) ) ばねつきの端 輪にかかる力の釣り合い 端の Lagrangian − K e 2 ϕ ( ℓ ) 2 -\frac{K_{\mathrm e}}2\phi(\ell)^2 − 2 K e ϕ ( ℓ ) 2 を足す エネルギー 運動と伸びのエネルギーを足す H = π ∂ t ϕ − L \mathcal H=\pi\,\partial_t\phi-\mathcal L H = π ∂ t ϕ − L エネルギーの流れ 左の部分が右の部分にする仕事率 S = − Y ∂ x ϕ ∂ t ϕ \mathcal S=-Y\partial_x\phi\,\partial_t\phi S = − Y ∂ x ϕ ∂ t ϕ 3次元への拡張 膜・音波ごとに力の釣り合いを立て直す L \mathcal L L の ∂ x \partial_x ∂ x を ∇ \nabla ∇ にする
自由端の条件は、授業では「端の輪は質量が無いので、縦の力の合力が 0 0 0 でなければ無限の加速度を受ける」という議論で出す。この章では、端の変位を自由に動かす変分から機械的に出た。エネルギーの流れも、授業では「点 x x x で左の部分が右の部分に及ぼす縦の力 − T s ∂ x ϕ -T_{\mathrm s}\partial_x\phi − T s ∂ x ϕ と、その点の速度 ∂ t ϕ \partial_t\phi ∂ t ϕ の積」として出す。これは式 (2.5.27) の S \mathcal S S と同じである。2つの書き方から得られる物理は同じである。場の書き方の利点は、膜・音波・電磁場・量子力学の波動関数へ移るとき、毎回力の釣り合いを考え直さずに L \mathcal L L を取り替えるだけで済むことにある。
例題 2.5.6 進む波が運ぶ仕事率 ★
オリジナル 例題 2.5.1 の弦(ρ = 0.10 k g m − 1 \rho=0.10\,\mathrm{kg\,m^{-1}} ρ = 0.10 kg m − 1 、Y = 10 N Y=10\,\mathrm N Y = 10 N )の上を、進む波 ϕ = A cos ( k x − ω t ) \phi=A\cos(kx-\omega t) ϕ = A cos ( k x − ω t ) (A = 1.0 × 10 − 3 m A=1.0\times10^{-3}\,\mathrm m A = 1.0 × 1 0 − 3 m 、ω = 126 s − 1 \omega=126\,\mathrm s^{-1} ω = 126 s − 1 )が右へ進む。
(1) 点 x x x で左の部分が右の部分にする仕事率を、力×速度として求め、S \mathcal S S と一致することを確かめよ。
(2) その時間平均を求め、平均のエネルギー密度×波の速さに等しいことを示せ。数値はいくらか。
解答 (1) 点 x x x で左の部分が右の部分に及ぼす縦の力は − Y ∂ x ϕ -Y\partial_x\phi − Y ∂ x ϕ (右の部分が上へ傾いて ∂ x ϕ > 0 \partial_x\phi>0 ∂ x ϕ > 0 なら、左の部分は右の部分を下へ引く)、その点の速度は ∂ t ϕ \partial_t\phi ∂ t ϕ なので、仕事率は − Y ∂ x ϕ ∂ t ϕ = S -Y\partial_x\phi\,\partial_t\phi=\mathcal S − Y ∂ x ϕ ∂ t ϕ = S である。この波では
S = − Y ( − A k sin θ ) ( A ω sin θ ) = Y A 2 k ω sin 2 θ , θ = k x − ω t . \mathcal S=-Y(-Ak\sin\theta)(A\omega\sin\theta)=YA^2k\omega\sin^2\theta,\qquad \theta=kx-\omega t . S = − Y ( − A k sin θ ) ( A ω sin θ ) = Y A 2 k ω sin 2 θ , θ = k x − ω t .
(2) sin 2 \sin^2 sin 2 の平均 1 / 2 1/2 1/2 と Y k = ρ v 2 k = ρ v ω Yk=\rho v^2k=\rho v\omega Y k = ρ v 2 k = ρ v ω を使って S ‾ = 1 2 Y A 2 k ω = 1 2 ρ v ω 2 A 2 \overline{\mathcal S}=\frac12YA^2k\omega=\frac12\rho v\omega^2A^2 S = 2 1 Y A 2 k ω = 2 1 ρ v ω 2 A 2 。例題 2.5.1 の平均のエネルギー密度 1 2 ρ ω 2 A 2 \frac12\rho\omega^2A^2 2 1 ρ ω 2 A 2 に v v v を掛けたものに等しく、エネルギーが波の速さで運ばれていることを示す。数値は S ‾ = 7.90 × 10 − 4 J m − 1 × 10 m s − 1 = 7.90 × 10 − 3 W \overline{\mathcal S}=7.90\times10^{-4}\,\mathrm{J\,m^{-1}}\times10\,\mathrm{m\,s^{-1}}=7.90\times10^{-3}\,\mathrm W S = 7.90 × 1 0 − 4 J m − 1 × 10 m s − 1 = 7.90 × 1 0 − 3 W 。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com