tecktom.com

授業で習った形との対応

2026年10月15日6分で読めますAI生成教科書『虎の巻』第2部 第1章 L=T-V の例題

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第1章 §2.1.6 ★★ 換算質量の引き算・Kepler の保存則・Euler 方程式・連成振動の行列を、この章の LL からの導き方と並べる

解析力学を使わない力学の授業では、この章の4つの系を Newton の運動方程式から次のように扱う。

振動子:mx¨=−kx+F(t)m\ddot x=-kx+F(t) を立て、強制振動は特解と一般解の和で解く。周期的な力 F0cos⁡ωtF_0\cos\omega t への応答の振幅が ω→ω0\omega\to\omega_0 で大きくなる共鳴を主に扱う。この章の式 (2.1.8) は、一般の時刻依存の力について、終わったときのエネルギーを Fourier 成分1つで書いたもので、共鳴の式と同じ内容を時間領域で表している。

2体:m1r¨1=Fm_1\ddot{\bm{r}}_1=\bm{F}, m2r¨2=−Fm_2\ddot{\bm{r}}_2=-\bm{F}(作用反作用)を足して重心の等速運動、それぞれを質量で割って引いてr¨=(1/m1+1/m2)F\ddot{\bm{r}}=(1/m_1+1/m_2)\bm{F}、つまり mrr¨=Fm_{\mathrm r}\ddot{\bm{r}}=\bm{F} を出す。この章では LL の運動エネルギーを座標変換 (2.1.1) で書き直しただけで同じ分離が出て、作用反作用は VV が r1−r2\bm{r}_1-\bm{r}_2 だけによることの帰結になる(第1部第10章の対応の節)。

Kepler 問題:面積速度一定(角運動量)とエネルギーの2つの保存則から動径方向の式を作り、u=1/ru=1/r に変えて Binet の式を積分して円錐曲線を出す。この章の Runge–Lenz ベクトル (2.1.10) は、その積分の積分定数(離心率と近点の向き)を、初めから保存量として持っていた量である。

剛体:角運動量の時間変化がトルクに等しい dLdt=N\frac{\mathrm{d} \bm{L}}{\mathrm{d} t}=\bm{N} を、剛体に固定した主軸の成分で書く。剛体に固定した座標系は回っているので、L=(I1ω1,I2ω2,I3ω3)\bm{L}=(I_1\omega_1,I_2\omega_2,I_3\omega_3) の成分の時間微分に ω×L\bm{\omega}\times\bm{L} が足され、

I1ω˙1=(I2−I3)ω2ω3+N1,I2ω˙2=(I3−I1)ω3ω1+N2,I3ω˙3=(I1−I2)ω1ω2+N3(2.1.25)I_1\dot\omega_1=(I_2-I_3)\omega_2\omega_3+N_1,\quad I_2\dot\omega_2=(I_3-I_1)\omega_3\omega_1+N_2,\quad I_3\dot\omega_3=(I_1-I_2)\omega_1\omega_2+N_3\tag{2.1.25}

という Euler 方程式になる。この章では Euler 角を一般化座標にとって LL (2.1.14) を書いたので、Euler 方程式を通らずに保存量 pφp_\varphi, pψp_\psi, EE を直接得た。Euler 方程式は角速度の成分(一般化座標の時間微分ではない量)についての式で、LL から出すには剛体の回転の群の上で変分をとる別の道具が要る。

連成振動:ばねでつないだ質点の運動方程式を直交座標で mx¨=−Kxm\ddot{\bm{x}}=-K\bm{x} と行列で書き、x=ae−iωt\bm{x}=\bm{a}\mathrm{e}^{-\mathrm{i}\omega t} を入れて det⁡(K−mω2)=0\det(K-m\omega^2)=0 を解く。質量がそろっていれば普通の固有値問題で、この章の一般化固有値問題 (2.1.21) はそれを、質量がそろっていない場合や、振り子の角度のような直交座標でない一般化座標の場合へ広げたものである。剛性行列が VV の2階微分で、質量行列が TT の係数であるという形で書くと、どちらの場合も同じ手順で扱うことができる。

例題 2.1.6 本を投げ上げると中間の軸のまわりで裏返る ★★

オリジナル外からのトルクのない剛体(N=0\bm{N}=0)が、主慣性モーメント I1<I2<I3I_1<I_2<I_3 のうち I2I_2 の軸のまわりに角速度 ω0\omega_0 で回っている。

  1. (1) Euler 方程式 (2.1.25) を (ω1,ω2,ω3)=(η1,ω0,η3)(\omega_1,\omega_2,\omega_3)=(\eta_1,\omega_0,\eta_3)(η1\eta_1, η3\eta_3 は小さい)のまわりで1次まで展開し、ずれが eλt\mathrm{e}^{\lambda t} で増えることを示して λ\lambda を求めよ。
  2. (2) 縦 21 cm21\,\mathrm{cm}・横 15 cm15\,\mathrm{cm}・厚さ 3 cm3\,\mathrm{cm}、質量 0.6 kg0.6\,\mathrm{kg} の一様な直方体(本)を、横の辺に平行な軸のまわりに毎秒 33 回転で投げ上げる。λ\lambda と、ずれが 10001000 倍になるまでの時間を求めよ。
解答

(1) 1次までで ω˙2=0\dot\omega_2=0、

I1η˙1=(I2−I3)ω0η3,I3η˙3=(I1−I2)ω0η1.I_1\dot\eta_1=(I_2-I_3)\omega_0\eta_3,\qquad I_3\dot\eta_3=(I_1-I_2)\omega_0\eta_1 .

2つ目を微分して1つ目を入れると η¨3=λ2η3\ddot\eta_3=\lambda^2\eta_3、λ2=ω02(I3−I2)(I2−I1)I1I3\lambda^2=\omega_0^2\dfrac{(I_3-I_2)(I_2-I_1)}{I_1I_3}。I1<I2<I3I_1<I_2<I_3 なので λ2>0\lambda^2>0 で、ずれは eλt\mathrm{e}^{\lambda t} で増える。最大または最小の軸のまわりの回転では λ2<0\lambda^2<0 になり、ずれは振動するだけで増えない。

(2) 直方体の重心のまわりの主慣性モーメントは、辺 aa に平行な軸について M(b2+c2)/12M(b^2+c^2)/12 などで、縦の辺の軸 I1=0.6(0.152+0.032)/12=1.17×10−3I_1=0.6(0.15^2+0.03^2)/12=1.17\times10^{-3}、横の辺の軸 I2=0.6(0.212+0.032)/12=2.25×10−3I_2=0.6(0.21^2+0.03^2)/12=2.25\times10^{-3}、厚さの向きの軸 I3=0.6(0.212+0.152)/12=3.33×10−3I_3=0.6(0.21^2+0.15^2)/12=3.33\times10^{-3}(単位 kg m2\mathrm{kg\,m^2})。横の辺の軸が中間の軸である。ω0=2π×3=18.85 s−1\omega_0=2\pi\times3=18.85\,\mathrm{s^{-1}} で

λ=18.85(3.33−2.25)(2.25−1.17)1.17×3.33=10.3 s−1,\lambda=18.85\sqrt{\frac{(3.33-2.25)(2.25-1.17)}{1.17\times3.33}}=10.3\,\mathrm{s^{-1}},

10001000 倍になる時間は ln⁡1000/λ=0.67 s\ln1000/\lambda=0.67\,\mathrm s、およそ2回転である。投げ上げた本が空中で1回裏返るように見えるのはこの不安定のためで、Euler 方程式を数値積分すると、3秒のうちに ω2\omega_2 の符号が反転する。保存量の言葉では、エネルギー 12∑Iaωa2\frac12\sum I_a\omega_a^2 と角運動量の大きさ ∑Ia2ωa2\sum I_a^2\omega_a^2 がどちらも保存し、中間の軸の回転はその2つの等高面が交わる鞍点にあたる。

木下の勉強ステータス

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

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

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