よみもの:解析力学

解析力学（かいせきりきがく、analytical mechanics）は、主に直交座標を用いて記述されたニュートン力学を数学的に洗練し、様々な状況を簡単に数式で表すためにオイラー、ラグランジュ、ハミルトン等によって作られた力学の体系である。古典力学の集大成であり、解析力学的な考え方は統計力学、量子力学、相対性理論などの現代物理学に引き継がれている。

ここはそんな解析力学を一から導出してしまおう、という趣旨で書いてみた記事である。最終的には、ここにわかりやすい教科書のようなものが出来上がればいいな。

序論

一般化座標と一般化速度の定義

ある座標系に存在する全ての質点の位置座標$q_i$をずらーっと並べて書いた集合を$\{q_i\}$と表し、これを一般化座標と呼ぶ(たとえば、3次元直交座標系にN個の質点がある場合$\{q_i\}$の成分の数は3N個となる)。

一般化座標の各成分を時間で微分したものの集合、すなわち$\{{\dot{q}}_i\}$を一般化速度と呼ぶ。

このとき、各$q_i$および${\dot{q}}_i$はそれぞれ時刻$t$の関数であるから、それぞれ

$q_i=q_i(t)$

${\dot{q}}_i={\dot{q}}_i(t)$

と表せる。

とりあえず、ある時刻の$(\{q_i\})$を状態と呼び、$\{q_i(t)\}$の関数形を過程と呼ぶことにする。

保存力とポテンシャル

直交座標系$(\{x_j\})$において、位置$(\{x_j\})$と時刻$t$の関数$U=U(\{x_j\},t)$を用いて、力$f_j$が $f_j=-\frac{\partial U}{\partial x_j}$と表せるとき、この力を保存力と呼び、

そしてこのときの$U$をポテンシャルと呼ぶ。

また、保存力でない力$f_j=s_j$を非保存力と呼ぶ

ラグランジュ形式

解析力学の基礎には、物理学に様々な形で登場する変分原理の一つであるハミルトンの原理が登場する。変分原理は他にも、光に適用できるフェルマーの原理や、電磁気学におけるディリクレの原理などがある。このような考え方は量子力学の基礎にもなっている。ここでは、ニュートン力学とハミルトンの原理から、ニュートンの運動方程式よりも適用範囲の広いラグランジュ方程式を導出する。

ハミルトンの原理

非保存力が存在しない系が状態A(時刻:$t_A$)から状態B(時刻:$t_B$)に変化するとき、

$I=\underset{t_A}{\overset{t_B}{\int }}L(\{q_i\},\{{\dot{q}}_i\},t)dt$

とおき、いまのところ謎の関数$L=L(\{q_i\},\{{\dot{q}}_i\},t)$に含まれる$\{q_i\}$の各成分を、微小な$\delta q_i=\delta q_i(t)$だけ変化させるとき、$\{{\dot{q}}_i\}$の各成分も微小に$\delta {\dot{q}}_i=\frac{d}{dt}\delta q_i(t)$だけ変化することになる。

ただし、状態A、状態Bそれぞれのときの$\{q_i\}$はすでに決めたため、$\delta q_i(t_A)=\delta q_i(t_B)=0$である。

このとき、$L$の変化量は

$\delta L=L(\{q_i+\delta q_i\},\{{\dot{q}}_i+\delta {\dot{q}}_i\},t)-L(\{q_i\},\{{\dot{q}}_i\},t)$

である。さらに、このときの$I$の変化量は

$\delta I=\underset{t_A}{\overset{t_B}{\int }}L(\{q_i+\delta q_i\},\{{\dot{q}}_i+\delta {\dot{q}}_i\},t)dt-\underset{t_A}{\overset{t_B}{\int }}L(\{q_i\},\{{\dot{q}}_i\},t)dt$

である。以上の2つの式より

$\delta I=\underset{t_A}{\overset{t_B}{\int }}\delta Ldt$

である。 このとき、"$\delta I=0$となるような過程のみが実現する"と仮定する。

そして、この仮定をハミルトンの原理と呼び、$I$を作用と呼び、

ハミルトンの原理をみたす$L$をラグランジアンと呼ぶ。

ラグランジュ方程式

さて、$\{q_i\}$のある1成分$q_i$のみを$\delta q_i$だけ微小変化させるとき、${\dot{q}}_i$も微小に$\delta {\dot{q}}_i$だけ微小変化することになる。よって、このとき$L$は、

$\delta L=\frac{\partial L}{\partial q_i}\delta q_i+\frac{\partial L}{\partial {\dot{q}}_i}\delta {\dot{q}}_i$ だけ変化することになる。よって、

$\delta I=\underset{t_A}{\overset{t_B}{\int }}(\frac{\partial L}{\partial q_i}\delta q_i+\frac{\partial L}{\partial {\dot{q}}_i}\delta {\dot{q}}_i)dt$

となる。よって、

$\delta I=\underset{t_A}{\overset{t_B}{\int }}\frac{\partial L}{\partial q_i}\delta q_{i}dt+\underset{t_A}{\overset{t_B}{\int }}\frac{\partial L}{\partial {\dot{q}}_i}\delta {\dot{q}}_{i}dt$

である。ここで第2項にのみ部分積分を用いると、

$\delta I=\underset{t_A}{\overset{t_B}{\int }}\frac{\partial L}{\partial q_i}\delta q_{i}dt+[\frac{\partial L}{\partial {\dot{q}}_i}\delta q_i{]}_{t_A}^{t_B}-\underset{t_A}{\overset{t_B}{\int }}\frac{d}{dt}(\frac{\partial L}{\partial {\dot{q}}_i})\delta q_{i}dt$

となる。第1項と第3項をまとめ、第2項を展開すると、

$\delta I=\underset{t_A}{\overset{t_B}{\int }}(\frac{\partial L}{\partial q_i}-\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_i})\delta q_{i}dt+(\frac{\partial L}{\partial {\dot{q}}_i}(t_B)\delta q_i(t_B)-\frac{\partial L}{\partial {\dot{q}}_i}(t_A)\delta q_i(t_A))$

となる。ここで$\delta q_i(t_A)=\delta q_i(t_B)=0$だったから、後半部分が0となり、

$\delta I=\underset{t_A}{\overset{t_B}{\int }}(\frac{\partial L}{\partial q_i}-\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_i})\delta q_{i}dt$

となる。ここで、$\delta q_i=\delta q_i(t)$がどんな関数形だったとしても$\delta I=0$となるためには

$\frac{\partial L}{\partial q_i}-\frac{d}{dt}\frac{\partial L}{\partial {\dot{q}}_i}=0$

となっていればよい。この式をラグランジュ方程式と呼ぶ。

ラグランジュ方程式は、ハミルトンの原理と同じことを言っている式である。

ラグランジアンの具体的な形

直交座標を用いてある系の運動エネルギーの合計$T$を表すと、

$T=\sum _j\frac{1}{2}{m}_j{\dot{x_j}}^2$

となる。ここで$m_j$は、各${\dot{x}}_j$に対応する質量である。この式から、

$\frac{\partial T}{\partial {\dot{x}}_j}=m_j{\dot{x}}_j$

である。また、$T=T(\{{\dot{x}}_j(\{q_i\},\{{\dot{q}}_i\})\})$と表せることから、

$\frac{\partial T}{\partial {\dot{q}}_i}=\sum _j\frac{\partial {\dot{x}}_j}{\partial {\dot{q}}_i}\frac{\partial T}{\partial {\dot{x}}_j}$

である。また、$x_j=x_j(\{q_i(t)\})$と表せることから、

${\dot{x}}_j=\sum _i\frac{\partial x_j}{\partial q_i}{\dot{q}}_i$

である。よって

$\frac{\partial {\dot{x}}_j}{\partial {\dot{q}}_i}=\frac{\partial x_j}{\partial q_i}$

となる。以上のことから、

$\frac{\partial T}{\partial {\dot{q}}_i}=\sum _j\frac{\partial x_j}{\partial q_i}{m}_j{\dot{x}}_j$

である。ここで、$p_i=\frac{\partial T}{\partial {\dot{q}}_i}$とおき、これを一般化運動量と呼ぶ。すると、

$p_i=\sum _j\frac{\partial x_j}{\partial q_i}{m}_j{\dot{x}}_j$

となる。両辺を時間で微分すると、

${\dot{p}}_i=\frac{d}{dt}(\sum _j\frac{\partial x_j}{\partial q_i}{m}_j{\dot{x}}_j)$

となる。ここで、各$m_i$は時間変化しないから、

${\dot{p}}_i=\sum _j{m}_j\frac{d}{dt}({\dot{x}}_j\frac{\partial x_j}{\partial q_i})$

となる。ここで、時間微分の部分を展開してシグマを2つに分けると、

${\dot{p}}_i=\sum _j\frac{\partial x_j}{\partial q_i}{m}_j{\ddot{x}}_j+\sum _i\frac{\partial {\dot{x}}_j}{\partial q_i}{m}_j{\dot{x}}_j$

また、ニュートン力学の運動方程式より、

$f_j=m_j{\ddot{x}}_j$

である。以上のことから、

${\dot{p}}_i=\sum _j\frac{\partial x_j}{\partial q_i}{f}_j+\sum _j\frac{\partial {\dot{x}}_j}{\partial q_i}\frac{\partial T}{\partial {\dot{x}}_j}$

である。よって

${\dot{p}}_i=\sum _j\frac{\partial x_j}{\partial q_i}{f}_j+\frac{\partial T}{\partial q_i}$

である。ここで、$F_i=\sum _j\frac{\partial x_j}{\partial q_i}{f}_j$とおき、これを一般化力と呼ぶ。すると、

${\dot{p}}_i=F_i+\frac{\partial T}{\partial q_i}$

となる。ここで、力$f_i$が保存力$-\frac{\partial U}{\partial x_j}$と非保存力$s_j$で構成されているとき、

$F_i=\sum _j\frac{\partial x_j}{\partial q_i}(-\frac{\partial U}{\partial x_j}+s_j)$

であるから

$F_i=-\sum _j\frac{\partial x_j}{\partial q_i}\frac{\partial U}{\partial x_j}+\sum _j\frac{\partial x_j}{\partial q_i}{s}_j$

となり

$F_i=-\frac{\partial U}{\partial q_i}+\sum _j\frac{\partial x_j}{\partial q_i}{s}_j$

となる。ここで、非保存力$s_j$による一般化力$\sum _j\frac{\partial x_j}{\partial q_i}{s}_j$を$S_i$と表すと、

$F_i=-\frac{\partial U}{\partial q_i}+S_i$

となる。以上のことから、

${\dot{p}}_i=-\frac{\partial U}{\partial q_i}+S_i+\frac{\partial T}{\partial q_i}$

である。第1項と第3項をまとめると、

${\dot{p}}_i=\frac{\partial (T-U)}{\partial q_i}+S_i$

となる。ここで、一般化運動量の定義から

$\frac{d}{dt}\frac{\partial T}{\partial {\dot{q}}_i}=\frac{\partial (T-U)}{\partial q_i}+S_i$

である。また、$U=U(\{x_i(\{q_i\})\},t)$と表せるから、

$\frac{\partial U}{\partial {\dot{q}}_i}=0$

である。以上のことから、

$\frac{d}{dt}(\frac{\partial T}{\partial {\dot{q}}_i}-\frac{\partial U}{\partial {\dot{q}}_i})=\frac{\partial (T-U)}{\partial q_i}+S_i$

となり、よって

$\frac{d}{dt}\frac{\partial (T-U)}{\partial {\dot{q}}_i}-\frac{\partial (T-U)}{\partial q_i}=S_i$

である。ここで、非保存力$s_j$がすべて0のとき、非保存力の一般化力$S_i$もすべて0であるから、

$\frac{d}{dt}\frac{\partial (T-U)}{\partial {\dot{q}}_i}-\frac{\partial (T-U)}{\partial q_i}=0$

となる。これはラグランジュ方程式そのものであるから、$T-U$はラグランジアンであるといえる。

そこで、$L=T-U$とおくと、

$\frac{\partial L}{\partial {\dot{q}}_i}=\frac{\partial T}{\partial {\dot{q}}_i}-\frac{\partial U}{\partial {\dot{q}}_i}$

であり、$\frac{\partial U}{\partial {\dot{q}}_i}=0$であったから、

$\frac{\partial L}{\partial {\dot{q}}_i}=\frac{\partial T}{\partial {\dot{q}}_i}$

であり、

$\frac{\partial L}{\partial {\dot{q}}_i}=p_i$

であることがわかる。

ハミルトン形式

ニュートン力学の適用範囲を広げたラグランジュ形式であったが、ここでは、ラグランジュ方程式からさらに適用範囲を広げたハミルトンの正準方程式を導出する。その際、ルジャンドル変換という数学的テクニックが暗に用いられているが、そんなものは知らなくても理解できる。

ハミルトンの正準方程式

以下では、非保存力による一般化力$S_i$がない場合について考える。このとき、

$H=\sum _i{p}_i{\dot{q}}_i-L$

とおき、この$H$をハミルトニアンと呼ぶ。すると、$H$の微小変化$dH$は、

$dH=\sum _i(d{p}_i{\dot{q}}_i+p_{i}d{\dot{q}}_i)-dL$

であり、よって

$dH=\sum _i(d{p}_i{\dot{q}}_i+p_{i}d{\dot{q}}_i)-\sum _i\frac{\partial L}{\partial {\dot{q}}_i}d{\dot{q}}_i$

となる。ここで、前節より

$\frac{\partial L}{\partial {\dot{q}}_i}=p_i$

だったから、

$dH=\sum _i(d{p}_i{\dot{q}}_i+\frac{\partial L}{\partial {\dot{q}}_i}{\dot{q}}_i)-\sum _i\frac{\partial L}{\partial {\dot{q}}_i}d{\dot{q}}_i$

であり、

$dH=\sum _i{\dot{q}}_{i}d{p}_i$

である。よって

${\dot{q}}_i=\frac{\partial H}{\partial q_i}$

である。

また、ハミルトニアンの定義より、

$\frac{\partial H}{\partial q_i}=-\frac{\partial L}{\partial q_i}$

であり、また

$\frac{\partial L}{\partial {\dot{q}}_i}=p_i$

であったから

${\dot{p}}_i=-\frac{\partial H}{\partial q_i}$

となる。以上で

${\dot{q}}_i=\frac{\partial H}{\partial p_i}$

${\dot{p}}_i=-\frac{\partial H}{\partial q_i}$

の2式が導けた。この2式をハミルトンの正準方程式と呼ぶ。

参考文献

• よくわかる解析力学(前野昌弘 著)
• EMANの解析力学(広江克彦 著)