logo

物理学におけるオイラー‐ラグランジュ方程式 📂古典力学

物理学におけるオイラー‐ラグランジュ方程式

概要

この記事は、読者が古典力学のカテゴリー、具体的にはラグランジュ力学とハミルトンの最小作用の原理を読んだという前提で書かれています。可能な限り重複する記法や内容も再説明しますが、説明されていない記法はリンクを参照してください。

運動の経路にわたるラグランジアンの積分を作用と言い、JJで表されます。位置をyyとすると、

J=t1t2Ldt=t1t2L(y(t), y(t), t)dt J =\int_{t_{1}} ^{t_2} Ldt=\int_{t_{1}}^{t_2}L\big( y^{\prime}(t),\ y(t),\ t \big)dt

この場合、実際の物体の運動経路に対してこの作用が最小値を持つことがハミルトンの原理です。

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

作用が極値(極大または極小)を持つ実際の運動経路のラグランジアンLLは、次の方程式を満たします。

Lyddt(Ly)=0 \dfrac{ \partial L}{\partial y}-\dfrac{d}{dt} \left( \dfrac{\partial L}{\partial y^{\prime}} \right)=0

この方程式はオイラー・ラグランジュ方程式Euler-Lagrangeと言います。

説明

重要な点は、その逆が成り立たないことです。極大や極小で微分すると00となりますが、微分して00であったとしても、それが極大や極小を意味するわけではありません。微分して00であれば、変曲点の可能性もあります。同様に、実際の運動経路はオイラー・ラグランジュ方程式を満たしますが、オイラー・ラグランジュ方程式を満たすすべての経路が実際の運動経路であるわけではありません。

証明

5D21C02B1.jpg

可能な任意の運動経路をy(α,t)=y(0,t)+αη(t)y(\alpha, t)=y(0,t)+\alpha\eta (t)としましょう。ここでy(0,t)y(0,t)は実際の運動経路です。また、すべての経路の始点と終点は同じなので、η(t1)=η(t2)=0\eta (t_{1})=\eta (t_2)=0が成り立ちます。そして、任意の運動経路に対するラグランジアンとそれに対する作用を次のように考えます。

L=L(y(α,t), y(α,t), t) L = L \big( y(\alpha, t)^{\prime},\ y(\alpha, t),\ t \big)

J(α)=t1t2L(y(α,t), y(α,t), t)dt J(\alpha) = \int_{t_{1}}^{t_2} L \big( y(\alpha, t)^{\prime},\ y(\alpha, t),\ t\big) dt

y(α,t)y(\alpha, t)α=0\alpha=0のとき実際の運動経路となり、実際の運動経路のとき作用は極値を持つので、α=0\alpha=0で微分すると00となります。

Jαα=0=0 \dfrac{\partial J}{\partial \alpha} \bigg|_{\alpha=0}=0

この式は、教科書や他の資料で簡単にδJ=0\delta J=0として表されます。(δ=α(\delta=\frac{\partial}{\partial \alpha}) 実際に上記の微分を計算すると、以下の式が得られます。

Jα=αt1t2L(y(α,t), y(α,t), t)dt=t1t2αL(y(α,t), y(α,t), t)dt=t1t2(Lyyα+Lyyα+Lttα)dt \begin{align*} \dfrac{\partial J}{\partial \alpha} &= \dfrac{\partial}{\partial \alpha} \int_{t_{1}}^{t_2} L \big( y(\alpha, t)^{\prime},\ y(\alpha, t),\ t \big) dt \\ &= \int_{t_{1}}^{t_2} \dfrac{\partial}{\partial \alpha} L \big( y(\alpha, t)^{\prime},\ y(\alpha, t),\ t \big) dt \\ &= \int_{t_{1}}^{t_2} \left( \dfrac{\partial L}{\partial y^{\prime}}\dfrac{\partial y^{\prime}}{\partial \alpha}+\dfrac{\partial L}{\partial y}\dfrac{\partial y}{\partial \alpha} +\dfrac{\partial L}{\partial t}\dfrac{\partial t}{\partial \alpha} \right) dt \end{align*}

このとき、yy^{\prime}{}^{\prime}は、ttに対する微分を意味することに注意してください。y(α,t)=y(0,t)+αη(t)y(\alpha, t)=y(0,t)+\alpha\eta (t)であるため、y(α,t)=y(0,t)+αη(t)y^{\prime}(\alpha, t)=y^{\prime}(0,t)+\alpha\eta^{\prime}(t)が成り立ちます。したがって、以下のようになります。

yα=η(t),yα=η(t),tα=0 \dfrac{\partial y^{\prime}}{\partial \alpha}=\eta^{\prime}(t), \quad \dfrac{\partial y}{\partial \alpha}=\eta (t),\quad \dfrac{\partial t}{\partial \alpha}=0

ttα\alphaは独立変数なので、微分すると00となります。したがって、以下のようになります。

Jα=t1t2(Lyη(t)+Lyη(t))dt \dfrac{\partial J}{\partial \alpha} = \int_{t_{1}}^{t_2} \left( \dfrac{\partial L}{\partial y^{\prime}}\eta^{\prime}(t)+\dfrac{\partial L}{\partial y}\eta (t) \right) dt

最初の項だけを部分積分で計算すると、その結果は以下のようになります。

t1t2Lyη(t)dt=Lyη(t)]t1t2t1t2ddt(Ly)η(t)dt=t1t2ddt(Ly)η(t)dt \begin{align*} \int_{t_{1}}^{t_2} \dfrac{\partial L}{\partial y^{\prime}}\eta^{\prime}(t) dt &= \left. \dfrac{\partial L}{\partial y^{\prime}}\eta (t) \right]_{t_{1}}^{t_2}-\int_{t_{1}}^{t_2} \dfrac{d}{dt}\left( \dfrac{\partial L}{\partial y^{\prime}}\right) \eta (t) dt \\ &= -\int_{t_{1}}^{t_2} \dfrac{d}{dt}\left( \dfrac{\partial L}{\partial y^{\prime}}\right) \eta (t) dt \end{align*}

定積分項が消える理由は、η(t1)=η(t2)=0\eta (t_{1})=\eta (t_2)=0だからです。これを元の式に代入し、η(t)\eta (t)に関して整理すると、次の式が得られます。

Jα=t1t2(ddtLy+Ly)η(t)dt \dfrac{\partial J}{\partial \alpha} = \int_{t_{1}}^{t_2} \left(- \dfrac{d}{dt}\dfrac{\partial L}{\partial y^{\prime}}+\dfrac{\partial L}{\partial y} \right) \eta (t)dt

しかし、Jα\dfrac{\partial J}{\partial \alpha}の値は、経路が実際の運動経路であるとき、(α=0(\alpha=0、つまり))00となります。したがって、α=0\alpha=0とすると、以下の式が得られます。

Jα=t1t2(ddtLy+Ly)η(t)dt=0 \dfrac{\partial J}{\partial \alpha} = \int_{t_{1}}^{t_2} \left(- \dfrac{d}{dt}\dfrac{\partial L}{\partial y^{\prime}}+\dfrac{\partial L}{\partial y} \right) \eta (t)dt=0

この式は、実際の経路に対するラグランジアンの作用を表し、あらゆるη(t)\eta (t)に対して常に成立しなければなりません。つまり、どのようなη(t)\eta (t)に対しても積分値が00と出なければなりません。これが可能であるためには、かけられた括弧内の値が00でなければならないということです。したがって、実際の経路に対するラグランジアンは、以下の方程式を満たします。

Lyddt(Ly)=0 \dfrac{ \partial L}{\partial y}-\dfrac{d}{dt} \left( \dfrac{\partial L}{\partial y^{\prime}} \right)=0

参照