オイラー＝ラグランジュ方程式

 

オイラー＝ラグランジュ方程式

解析力学とは力学の計算の手法として生み出されたものの一つであり、ニュートンの運動方程式のように本質的に新たなことが出てくるわけではない。

 

ラグランジアン

まずはラグランジアンの導入です。ラグランジアン$L$は以下のように表します。

ラグランジアンとは　 $L=T-U$　（運動エネルギー $-$ 位置エネルギー）という量です。

$L=L(q(t) , \dot q(t) , t)$　のように依存する変数を明示して書くこともある。

 

$q(t)$  ……具体的に書くと$q_{1}(t) , q_{2}(t) , \cdots  , q_{N}(t)$で$N$自由度の力学変数を表している。

例えば、$N=3$のとき、

直交座標では$q_{1}(t)=x(t) , q_{2}(t)=y(t) ,q_{3}(t)=z(t)$

極座標では $q_{1}(t)=r(t) , q_{2}(t)=\theta(t) , q_{3}(t)=\phi(t)$

$\dot q(t)$ 　……　これは上で述べた力学変数の時間微分です。 

$t$　　……　時間です。

最小作用の原理

解析力学の根本となる原理で、$t_{1}$～$t_{2}$ における運動は、$L$をラグランジアンとして

$S[q]=\displaystyle\int_{t_{1}}^{t_{2}}L dt$　を最小（停留）するような運動になる。

 

オイラー＝ラグランジュ方程式

上の最小作用の原理から以下の式が導かれる。

$\displaystyle\frac{\partial L}{\partial q_{i}}-\displaystyle\frac{d}{dt }\displaystyle\frac{\partial L}{\partial \dot q_{i}}=0$

これを、オイラー＝ラグランジュ方程式といい、運動を決定する方程式となる。

証明

$\delta S[q]=S[q+\delta q]-S[q]=$$\displaystyle\int_{t_{1}}^{t_{2}}L(q+\delta q , \dot q+\delta\dot q , t)dt-\displaystyle\int_{t_{1}}^{t_{2}}L(q , \dot q , t)dt$

 

$L(q+\delta q , \dot q+\delta\dot q , t)=L(q , \dot q , t)+\displaystyle\frac{\partial L(q , \dot q , t)}{\partial q_{i}}\delta q_{i}+\displaystyle\frac{\partial L(q , \dot q , t)}{\partial \dot q_{i}}\delta\dot q_{i}$ 

という展開公式（多変数のテイラー展開公式）を使用すると（第四項以降は無視）

 

$=\displaystyle\int_{t_{1}}^{t_{2}}\displaystyle\frac{\partial L(q , \dot q , t)}{\partial q_{i}}\delta q_{i}dt+\displaystyle\int_{t_{1}}^{t_{2}}\displaystyle\frac{\partial L(q , \dot q , t)}{\partial \dot q_{i}}\displaystyle\frac{d}{dt}\delta q_{i}dt$ 　（$\delta\dot q_{i}=\displaystyle\frac{d}{dt}\delta q_{i}$）

$=\displaystyle\int_{t_{1}}^{t_{2}}\displaystyle\frac{\partial L}{\partial q_{i}}\delta q_{i}dt+\biggl[\displaystyle\frac{\partial L}{\partial \dot q_{i}}\delta q_{i}\biggr]_{t_{1}}^{t_{2}}-\displaystyle\int_{t_{1}}^{t_{2}}\displaystyle\frac{d}{dt}\displaystyle\frac{\partial  L}{\partial \dot q}\delta q_{i}dt$　（第二項部分積分）

$=\displaystyle\int_{t_{1}}^{t_{2}}\biggl(\displaystyle\frac{\partial L}{\partial q_{i}}-\displaystyle\frac{d}{dt}\displaystyle\frac{\partial  L}{\partial \dot q}\biggr)\delta q_{i}dt+\biggl[\displaystyle\frac{\partial L}{\partial \dot q_{i}}\delta q_{i}\biggr]_{t_{1}}^{t_{2}}$　　（式を整理した）

 

ここで$\delta q_{1}=\delta q_{2}=0$ なので（微小変化で力学変数は変化しない。）最後の項は0になる。よって最小作用の原理により、この$\delta S[q]$が$0$になることから

 

$\displaystyle\frac{\partial L}{\partial q_{i}}-\displaystyle\frac{d}{dt }\displaystyle\frac{\partial L}{\partial \dot q_{i}}=0$