tecktom.com

この章のまとめ

2026年10月15日9分で読めますAI生成教科書『虎の巻』第1部 第4章 Euler–Lagrange 方程式と境界条件(数学だけ)

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第4章 §1.4.7 ★★★ Euler–Lagrange 方程式・第一積分・自然境界条件の3つで、端の条件まで含めた変分問題を解ける

汎関数 S[q]=∫L(q,q˙,t) dtS[q]=\int L(q,\dot q,t)\,\mathrm{d} t の第1変分を、連鎖律と部分積分で ϵ[∂L/∂q˙ η]+ϵ∫(∂L/∂q−ddt∂L/∂q˙)η dt\epsilon[\partial L/\partial\dot q\,\eta]+\epsilon\int(\partial L/\partial q-\frac{\mathrm{d} }{\mathrm{d} t}\partial L/\partial\dot q)\eta\,\mathrm{d} t の形に書いた。積分の中の括弧が汎関数微分で、端を固定した問題ではそれを 00 とおいた Euler–Lagrange 方程式が停留の条件になる。ddt\frac{\mathrm{d} }{\mathrm{d} t} は経路に沿った微分で、∂2L/∂q˙2≠0\partial^2L/\partial\dot q^2\neq0 なら方程式は2階になり、条件が2つ要る。自由度が複数あれば成分ごとに1本ずつ方程式が立つ。LL が座標 qjq_j を含まなければ ∂L/∂q˙j\partial L/\partial\dot q_j が、時刻を含まなければエネルギー関数 h=∑q˙i∂L/∂q˙i−Lh=\sum\dot q_i\partial L/\partial\dot q_i-L が解に沿って一定で、方程式の階数が1つ下がる。LL に時刻の関数の全微分 dΛ/dt\mathrm{d}\Lambda/\mathrm{d} t を足しても方程式は変わらない。端の値を自由にすると、部分積分の端の項から自然境界条件 ∂L/∂q˙=0\partial L/\partial\dot q=0 が出て、手放した端の値の代わりの条件になる。端に費用 Ψ\Psi があれば条件は ∂L/∂q˙+Ψ′=0\partial L/\partial\dot q+\Psi'=0 に変わる。回転面の面積と吊るした鎖からは懸垂線 acosh⁡((x−b)/a)a\cosh((x-b)/a) が、降下時間からはサイクロイドが、どちらもエネルギー関数の第一積分から出た。懸垂面の例では、停留する曲線が2本ある場合と1本もない場合があった。次の章(第1部第5章)では、物理の作用として L=T−VL=T-V を選ぶ理由と、「最小作用」が停留であることの意味を扱う。

この章の道具

  • 汎関数微分:δSδq(t)=∂L∂q−ddt∂L∂q˙\dfrac{\delta S}{\delta q(t)}=\frac{\partial L}{\partial q}-\frac{\mathrm{d} }{\mathrm{d} t}\frac{\partial L}{\partial \dot q}。Euler–Lagrange 方程式は ddt∂L∂q˙i−∂L∂qi=0\frac{\mathrm{d} }{\mathrm{d} t}\frac{\partial L}{\partial \dot q_i}-\frac{\partial L}{\partial q_i}=0。
  • 第一積分:∂L/∂qj=0\partial L/\partial q_j=0 なら ∂L/∂q˙j\partial L/\partial\dot q_j が一定。∂L/∂t=0\partial L/\partial t=0 なら h=∑iq˙i∂L/∂q˙i−Lh=\sum_i\dot q_i\partial L/\partial\dot q_i-L が一定。
  • L→L+dΛ(q,t)/dtL\to L+\mathrm{d}\Lambda(q,t)/\mathrm{d} t で方程式は変わらない。
  • 自由な端では ∂L/∂q˙=0\partial L/\partial\dot q=0(自然境界条件)。端の費用 Ψ\Psi があれば ∂L/∂q˙+Ψ′=0\partial L/\partial\dot q+\Psi'=0。

総合問題 終点の高さを決めない最速降下線 ★★★

オリジナル原点から静かに放した玉が、摩擦なしに曲線 y(x)y(x)(yy は下向き)に沿って滑り、鉛直な直線 x=Xx=X(X>0X>0)のどこかに着くまでの時間を最も短くする曲線を探す。終点の高さ y(X)y(X) は決めない。重力加速度を gg とし、玉の速さは v=2gyv=\sqrt{2gy} とする。

  1. (1) (使う道具:第1部第3章 の汎関数/1.4.1節)降下時間 T[y]\mathcal T[y] を汎関数として書き、汎関数微分 δT/δy(x)\delta\mathcal T/\delta y(x) を y,y′,y′′y,y',y'' で表せ。
  2. (2) (使う道具:1.4.2節)エネルギー関数の第一積分を使い、停留曲線が y(1+y′2)=2ry(1+y'^2)=2r を満たすことを示せ。
  3. (3) (使う道具:1.4.3節)x=Xx=X での自然境界条件を書け。
  4. (4) (使う道具:1.4.5節)停留曲線のサイクロイドの rr と終点の高さ、降下時間を XX と gg で表せ。X=1 mX=1\,\mathrm m, g=9.81 m s−2g=9.81\,\mathrm{m\,s^{-2}} での値を求め、例題 1.4.5 の終点 (1 m,1 m)(1\,\mathrm m,1\,\mathrm m) のサイクロイドの時間と比べよ。
  5. (5) (使う道具:第1部第3章 の離散化の考え方/近似の数理 第4章)比べる曲線を原点から (X,Y)(X,Y) への直線だけに限り、YY を動かして時間を最小にせよ。(4) と比べて何が言えるか。
解答

(1) T[y]=∫0X1+y′22gy dx\mathcal T[y]=\displaystyle\int_0^X\frac{\sqrt{1+y'^2}}{\sqrt{2gy}}\,\mathrm{d} x、y(0)=0y(0)=0。L=1+y′2 (2gy)−1/2L=\sqrt{1+y'^2}\,(2gy)^{-1/2} について

∂L∂y=−1+y′22y2gy,∂L∂y′=y′2gy1+y′2.\frac{\partial L}{\partial y}=-\frac{\sqrt{1+y'^2}}{2y\sqrt{2gy}},\qquad \frac{\partial L}{\partial y'}=\frac{y'}{\sqrt{2gy}\sqrt{1+y'^2}} .

後者を xx で微分すると、ddxy′1+y′2=y′′(1+y′2)3/2\frac{\mathrm{d} }{\mathrm{d} x}\dfrac{y'}{\sqrt{1+y'^2}}=\dfrac{y''}{(1+y'^2)^{3/2}} と ddx(2gy)−1/2=−y′2y(2gy)−1/2\frac{\mathrm{d} }{\mathrm{d} x}(2gy)^{-1/2}=-\dfrac{y'}{2y}(2gy)^{-1/2} から

ddx∂L∂y′=12gy(y′′(1+y′2)3/2−y′22y1+y′2).\frac{\mathrm{d} }{\mathrm{d} x}\frac{\partial L}{\partial y'}=\frac{1}{\sqrt{2gy}}\Bigl(\frac{y''}{(1+y'^2)^{3/2}}-\frac{y'^2}{2y\sqrt{1+y'^2}}\Bigr).

よって式 (1.4.4) より

δTδy(x)=−12gy(1+y′22y−y′22y1+y′2+y′′(1+y′2)3/2)=−12gy (1+y′2)3/2(1+y′22y+y′′).\frac{\delta\mathcal T}{\delta y(x)}=-\frac{1}{\sqrt{2gy}}\Bigl(\frac{\sqrt{1+y'^2}}{2y}-\frac{y'^2}{2y\sqrt{1+y'^2}}+\frac{y''}{(1+y'^2)^{3/2}}\Bigr) =-\frac{1}{\sqrt{2gy}\,(1+y'^2)^{3/2}}\Bigl(\frac{1+y'^2}{2y}+y''\Bigr).

(括弧の第1・2項は (1+y′2)−y′22y1+y′2=12y1+y′2\frac{(1+y'^2)-y'^2}{2y\sqrt{1+y'^2}}=\frac{1}{2y\sqrt{1+y'^2}} にまとまる。)停留の条件は 2yy′′+1+y′2=02yy''+1+y'^2=0 である。

(2) LL は xx を陽に含まないので、1.4.5節と同じ計算で y′∂L/∂y′−L=−1/(2gy1+y′2)y'\partial L/\partial y'-L=-1/(\sqrt{2gy}\sqrt{1+y'^2}) が一定となり、これを −1/4gr-1/\sqrt{4gr} とおくと y(1+y′2)=2ry(1+y'^2)=2r を得る。これが (1) の2階の方程式の第一積分になっていることは、y(1+y′2)y(1+y'^2) を xx で微分すると y′(1+y′2)+2yy′y′′=y′(1+y′2+2yy′′)y'(1+y'^2)+2yy'y''=y'(1+y'^2+2yy'') となり、(1) の停留条件の下で 00 になることからも確かめられる。

(3) 右端の値が自由なので、自然境界条件 (1.4.11) は

∂L∂y′∣x=X=y′2gy1+y′2∣x=X=0,すなわちy′(X)=0.\frac{\partial L}{\partial y'}\Big|_{x=X}=\frac{y'}{\sqrt{2gy}\sqrt{1+y'^2}}\Big|_{x=X}=0,\qquad\text{すなわち}\qquad y'(X)=0 .

停留曲線は終点で水平になる。

(4) サイクロイド (1.4.18) の傾きは y′=cot⁡(θ/2)y'=\cot(\theta/2) なので、y′=0y'=0 は θ=π\theta=\pi、すなわちサイクロイドの最も低い点である。X=r(π−sin⁡π)=πrX=r(\pi-\sin\pi)=\pi r から r=X/πr=X/\pi。終点の高さは y(X)=2r=2X/πy(X)=2r=2X/\pi、時間は式 (1.4.19) から

T∗=πrg=πXg.\mathcal T_*=\pi\sqrt{\frac rg}=\sqrt{\frac{\pi X}{g}} .

X=1 mX=1\,\mathrm m では r≈0.318 mr\approx0.318\,\mathrm m、終点の深さ 2/π≈0.637 m2/\pi\approx0.637\,\mathrm m、T∗=π/9.81 s≈0.566 s\mathcal T_*=\sqrt{\pi/9.81}\,\mathrm s\approx0.566\,\mathrm s である。例題 1.4.5 の終点を (1 m,1 m)(1\,\mathrm m,1\,\mathrm m) に固定したサイクロイドは 0.583 s0.583\,\mathrm s だった。終点を固定した曲線は終点を自由にした曲線の一部なので、自由な問題の最小値はそれより大きくならず、実際に約 3%3\% 短い。

(5) 原点から (X,Y)(X,Y) への直線では、長さ ℓ=X2+Y2\ell=\sqrt{X^2+Y^2}、斜面に沿った加速度 gY/ℓgY/\ell で、ℓ=12(gY/ℓ)T2\ell=\frac12(gY/\ell)\mathcal T^2 から

T(Y)=2(X2+Y2)gY.\mathcal T(Y)=\sqrt{\frac{2(X^2+Y^2)}{gY}} .

(X2+Y2)/Y=X2/Y+Y(X^2+Y^2)/Y=X^2/Y+Y は Y=XY=X で最小値 2X2X をとるので、最良の直線は 45∘45^\circ の直線で、T=4X/g\mathcal T=\sqrt{4X/g} である。X=1 mX=1\,\mathrm m では 0.639 s0.639\,\mathrm s で、例題 1.4.5 の直線と同じ値になる。比は T∗/4X/g=π/4≈0.886\mathcal T_*/\sqrt{4X/g}=\sqrt{\pi/4}\approx0.886 で、直線に限った最小値はサイクロイドより約 13%13\% 長い。比べる曲線を直線という狭い集まりに限ると、その中の最小値は本当の最小値の上界になる。これは、探す関数の空間を有限次元の部分空間に限る A 型(射影)の近似で、上界になることを近似の数理 第4章が Rayleigh–Ritz の方法として扱っている。直線の族は1つの数 YY で表されるので、関数の変分問題が1変数の普通の極値問題に置き換わった。

木下の勉強ステータス

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

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

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