tecktom.com

この章で選ぶ LL:4つの系の T−VT-V

2026年10月15日11分で読めますAI生成教科書『虎の巻』第2部 第1章 L=T-V の例題
目次
  1. 4つの Lagrangian
  2. 2体:重心と相対座標
  3. 剛体:重心の並進と回転(Koenig の分解)
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第1章 §2.1.1 ★★★ 4つの系それぞれについて、一般化座標と L=T−VL=T-V を書き下す。2体は重心と相対座標に、剛体は重心の並進と回転に分ける

第1部では、Lagrangian L=T−VL=T-V を選ぶ理由(第1部第5章)、Euler–Lagrange 方程式(第1部第4章)、対称性から保存量を取り出す Noether の定理(第1部第10章)、Legendre 変換で作る Hamiltonian と相空間(第1部第7章・第1部第8章)を、主に1粒子の簡単な例で作ってきた。第2部はこの道具を古典の場と相対論へ広げる部である。その前にこの章で、道具を力学の代表的な4つの系に続けて当て、「LL を書く → 循環座標と保存量を探す → 1自由度に落とす → 平衡点のまわりで2次まで展開する」という手順を一通り使えるようにする。

選ぶ系は次の4つである。第一に1次元の振動子で、平衡点のまわりのどんな系も2次まで展開すればこの形になるので、手順の行き着く先として最初に置く。外から時刻に依存する力を加えると LL が時刻を含み、Hamiltonian が保存しなくなる例にもなる。第二に中心力を及ぼし合う2体で、重心と相対運動に分けると1体の問題になり、回転対称性から角運動量が保存して平面運動になる。Kepler 問題では、回転対称性とは別にもう1つのベクトルの保存量(Runge–Lenz ベクトル)があり、これが軌道の閉じる理由になっている。第三に剛体で、運動エネルギーを重心の並進と重心のまわりの回転に分け、回転の部分を慣性テンソルで書く。重力の中の対称なこまでは Euler 角のうち2つが循環座標になり、残る1自由度の有効ポテンシャルから歳差運動の角速度が決まる。第四に多自由度の小振動で、質量行列と剛性行列を同時に対角化する基準振動を作り、直線形の分子 CO2_2 の振動数を計算して観測値と比べる。

最後に、Newton の運動方程式から同じ結果を出す授業の書き方(換算質量の引き算・Euler 方程式・連成振動の行列)と、この章の書き方を並べる。

4つの Lagrangian

この章で選ぶ Lagrangian を先に並べる。どれも L=T−VL=T-V の形で、違うのは一般化座標の選び方と、TT・VV の具体形だけである。

系一般化座標L=T−VL=T-V
振動子変位 xxm2x˙2−m2ω02x2+F(t) x\dfrac m2\dot x^2-\dfrac m2\omega_0^2x^2+F(t)\,x
中心力の2体重心 R\bm{R}、相対位置 r\bm{r}M2R˙ 2+mr2r˙ 2−V(∣r∣)\dfrac M2\dot{\bm{R}}^{\,2}+\dfrac{m_{\mathrm r}}2\dot{\bm{r}}^{\,2}-V(\lvert \bm{r}\rvert )
剛体重心 R\bm{R}、向きを決める3つの角M2R˙ 2+12∑a=13Iaωa2−V\dfrac M2\dot{\bm{R}}^{\,2}+\dfrac12\sum_{a=1}^3I_a\omega_a^2-V
小振動平衡点からのずれ ηi\eta_i12∑i,jMijη˙iη˙j−12∑i,jKijηiηj\dfrac12\sum_{i,j}M_{ij}\dot\eta_i\dot\eta_j-\dfrac12\sum_{i,j}K_{ij}\eta_i\eta_j

振動子の F(t)F(t) は外から加える時刻の関数の力で、−F(t)x-F(t)x を時刻に依存するポテンシャルとして VV に含めた。MM は全質量、mrm_{\mathrm r} は換算質量、IaI_a は主慣性モーメント、ωa\omega_a は剛体に固定した軸の方向の角速度の成分である。小振動の MijM_{ij} と KijK_{ij} は質量行列と剛性行列で、2.1.2節で作る。表の2行目と3行目の形は自明でなく、どちらも「重心を原点に選ぶと交差項が消える」という同じ事実から出るので、次の2つの小節で続けて導く。

2体:重心と相対座標

質量 m1m_1, m2m_2 の2粒子が、距離だけによるポテンシャル V(∣r1−r2∣)V(\lvert \bm{r}_1-\bm{r}_2\rvert ) で相互作用するとする。新しい一般化座標として

R=m1r1+m2r2M,r=r1−r2,M=m1+m2(2.1.1)\bm{R}=\frac{m_1\bm{r}_1+m_2\bm{r}_2}{M},\qquad \bm{r}=\bm{r}_1-\bm{r}_2,\qquad M=m_1+m_2\tag{2.1.1}

を選ぶと、r1=R+m2Mr\bm{r}_1=\bm{R}+\frac{m_2}{M}\bm{r}, r2=R−m1Mr\bm{r}_2=\bm{R}-\frac{m_1}{M}\bm{r} である。運動エネルギーに代入すると R˙⋅r˙\dot{\bm{R}}\cdot\dot{\bm{r}} の交差項の係数はm1m2M−m2m1M=0\frac{m_1m_2}{M}-\frac{m_2m_1}{M}=0 で消え、

L=M2R˙ 2+mr2r˙ 2−V(∣r∣),mr=m1m2m1+m2(2.1.2)L=\frac M2\dot{\bm{R}}^{\,2}+\frac{m_{\mathrm r}}2\dot{\bm{r}}^{\,2}-V(\lvert \bm{r}\rvert ),\qquad m_{\mathrm r}=\frac{m_1m_2}{m_1+m_2}\tag{2.1.2}

となる。1次元のばねで同じ分離をしたのが第1部第5章の総合問題で、3次元でも成分ごとに同じ計算になる。LL は R\bm{R} を含まないので R\bm{R} の3成分はすべて循環座標で、全運動量 MR˙M\dot{\bm{R}} が保存する(空間の並進対称性の Noether の保存量)。相対運動は、質量 mrm_{\mathrm r} の1粒子が固定した中心のまわりで V(r)V(r) を受けて動く問題になる。

剛体:重心の並進と回転(Koenig の分解)

剛体を NN 個の質点 mam_a の集まりとし、各質点の位置を重心 R\bm{R} と重心からの位置 ba\bm{b}_a に分けて ra=R+ba\bm{r}_a=\bm{R}+\bm{b}_a と書く。重心の定義から ∑amaba=0\sum_am_a\bm{b}_a=0 である。運動エネルギーは

T=12∑ama(R˙+b˙a)2=M2R˙ 2+R˙⋅ddt(∑amaba)+12∑amab˙a 2=M2R˙ 2+12∑amab˙a 2(2.1.3)T=\frac12\sum_am_a\bigl(\dot{\bm{R}}+\dot{\bm{b}}_a\bigr)^2=\frac M2\dot{\bm{R}}^{\,2}+\dot{\bm{R}}\cdot\frac{\mathrm{d} }{\mathrm{d} t}\Bigl(\sum_am_a\bm{b}_a\Bigr)+\frac12\sum_am_a\dot{\bm{b}}_a^{\,2} =\frac M2\dot{\bm{R}}^{\,2}+\frac12\sum_am_a\dot{\bm{b}}_a^{\,2}\tag{2.1.3}

と、重心の並進と重心に対する運動の和に分かれる(Koenig の分解)。中央の交差項は ∑amaba=0\sum_am_a\bm{b}_a=0 の時間微分なので消える。この分解は剛体でなくても成り立ち、2体の式 (2.1.2) はその N=2N=2 の場合にあたる(b1=m2Mr\bm{b}_1=\frac{m_2}{M}\bm{r}, b2=−m1Mr\bm{b}_2=-\frac{m_1}{M}\bm{r} を入れると 12mrr˙ 2\frac12m_{\mathrm r}\dot{\bm{r}}^{\,2} になる)。

剛体では、重心に対する各点の運動は全体の回転だけである。ある瞬間の角速度ベクトルを ω\bm{\omega} とすると b˙a=ω×ba\dot{\bm{b}}_a=\bm{\omega}\times\bm{b}_a と書け、∣ω×b∣2=∣ω∣2∣b∣2−(ω⋅b)2\lvert \bm{\omega}\times\bm{b}\rvert ^2=\lvert \bm{\omega}\rvert ^2\lvert \bm{b}\rvert ^2-(\bm{\omega}\cdot\bm{b})^2 を使うと

12∑amab˙a 2=12∑i,j=13Iij ωiωj,Iij=∑ama(∣ba∣2δij−ba,iba,j)(2.1.4)\frac12\sum_am_a\dot{\bm{b}}_a^{\,2}=\frac12\sum_{i,j=1}^3I_{ij}\,\omega_i\omega_j,\qquad I_{ij}=\sum_am_a\bigl(\lvert \bm{b}_a\rvert ^2\delta_{ij}-b_{a,i}b_{a,j}\bigr)\tag{2.1.4}

となる。IijI_{ij} を重心のまわりの慣性テンソルと呼ぶ。IijI_{ij} は実対称行列なので直交行列で対角化でき、対角化する座標軸を慣性主軸、対角成分 I1I_1, I2I_2, I3I_3 を主慣性モーメントと呼ぶ。剛体に固定した主軸の方向の角速度の成分を ω1\omega_1, ω2\omega_2, ω3\omega_3 とすれば、回転の運動エネルギーは12(I1ω12+I2ω22+I3ω32)\frac12(I_1\omega_1^2+I_2\omega_2^2+I_3\omega_3^2) である。ωa\omega_a を一般化座標の時間微分で書く式は2.1.4節で作る。

剛体の1点が固定されているときは、重心の代わりにその固定点を原点に選んでもよい。固定点が動かないので、ra=ba\bm{r}_a=\bm{b}_a と書いて同じ計算をすれば並進の項がない T=12∑IijωiωjT=\frac12\sum I_{ij}\omega_i\omega_j になり、慣性テンソルは固定点のまわりで計算したものになる。

例題 2.1.1 細い棒の振り子の LL を2通りに書く ★★

オリジナル長さ ℓ\ell、質量 MM の一様な細い棒の一端を天井の点 O に留め、鉛直面内で振らせる。鉛直下向きからの角度を θ\theta とする。

  1. (1) O から距離 ss の位置の小片 Mℓds\frac M\ell\mathrm{d} s の速さを書き、棒全体で積分して TT を求めよ。
  2. (2) 式 (2.1.3) の分解(重心の並進+重心のまわりの回転)で TT を求め、(1) と一致することを確かめよ。
  3. (3) L=T−VL=T-V を書き、小振動の角振動数を求めよ。ℓ=1 m\ell=1\,\mathrm m, g=9.81 m s−2g=9.81\,\mathrm{m\,s^{-2}} で、同じ長さの単振り子(質点)と比べよ。
解答

(1) O から距離 ss の点は半径 ss の円を角速度 θ˙\dot\theta で回るので、速さは sθ˙s\dot\theta である。

T=∫0ℓ12⋅Mℓ (sθ˙)2 ds=Mθ˙22ℓ⋅ℓ33=12⋅Mℓ23 θ˙2.T=\int_0^\ell\frac12\cdot\frac M\ell\,(s\dot\theta)^2\,\mathrm{d} s=\frac{M\dot\theta^2}{2\ell}\cdot\frac{\ell^3}{3}=\frac12\cdot\frac{M\ell^2}{3}\,\dot\theta^2 .

Mℓ2/3M\ell^2/3 は O のまわりの慣性モーメントで、式 (2.1.4) の固定点を原点に選んだ形にあたる。

(2) 重心は O から ℓ/2\ell/2 の点で、速さは ℓ2θ˙\frac\ell2\dot\theta。重心のまわりの慣性モーメントは ∫−ℓ/2ℓ/2Mℓu2 du=Mℓ212\int_{-\ell/2}^{\ell/2}\frac M\ell u^2\,\mathrm{d} u=\frac{M\ell^2}{12}、角速度は同じ θ˙\dot\theta である。

T=M2(ℓ2θ˙)2+12⋅Mℓ212θ˙2=12(Mℓ24+Mℓ212)θ˙2=12⋅Mℓ23θ˙2.T=\frac M2\Bigl(\frac\ell2\dot\theta\Bigr)^2+\frac12\cdot\frac{M\ell^2}{12}\dot\theta^2=\frac12\Bigl(\frac{M\ell^2}{4}+\frac{M\ell^2}{12}\Bigr)\dot\theta^2=\frac12\cdot\frac{M\ell^2}{3}\dot\theta^2 .

(1) と一致する。Mℓ23=M(ℓ/2)2+Mℓ212\frac{M\ell^2}{3}=M(\ell/2)^2+\frac{M\ell^2}{12} は、固定点のまわりの慣性モーメントが「重心のまわりの慣性モーメント+全質量×距離の2乗」になるという関係(平行軸の定理)の1例である。

(3) 重心の高さは −ℓ2cos⁡θ-\frac\ell2\cos\theta なので V=−Mgℓ2cos⁡θV=-\frac{Mg\ell}{2}\cos\theta、

L=Mℓ26θ˙2+Mgℓ2cos⁡θ.L=\frac{M\ell^2}{6}\dot\theta^2+\frac{Mg\ell}{2}\cos\theta .

Euler–Lagrange 方程式は Mℓ23θ¨=−Mgℓ2sin⁡θ\frac{M\ell^2}{3}\ddot\theta=-\frac{Mg\ell}{2}\sin\theta。sin⁡θ≈θ\sin\theta\approx\theta とすると角振動数はω=Mgℓ/2Mℓ2/3=3g2ℓ\omega=\sqrt{\dfrac{Mg\ell/2}{M\ell^2/3}}=\sqrt{\dfrac{3g}{2\ell}} で、ℓ=1 m\ell=1\,\mathrm m では 3.84 s−13.84\,\mathrm{s^{-1}} である。単振り子の g/ℓ=3.13 s−1\sqrt{g/\ell}=3.13\,\mathrm{s^{-1}} より速く、長さ 2ℓ/32\ell/3 の単振り子と同じになる。棒は O の近くの部分が小さい半径で回っていて慣性モーメントが Mℓ2M\ell^2 より小さく、その分だけ速く振れる。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ