オイラー=ラグランジュ方程式とは?公式・導出・具体例と解き方を徹底解説

目次
オイラー=ラグランジュ方程式とは?公式・導出・具体例と解き方を徹底解説
オイラー=ラグランジュ方程式とは?公式・導出・具体例と解き方を徹底解説
@ creator • Click to Play Video Inline
🎵 オイラー=ラグランジュ方程式とは?公式・導出・具体例と解き方を徹底解説

オイラー=ラグランジュ方程式は、物体がどのように運動するのかを、エネルギーと変分法を使って導き出すための基本公式です。ニュートンの運動方程式が力と加速度の関係から運動を記述するのに対し、この方程式は「作用積分」と呼ばれる量に注目し、物理的な運動を支配する微分方程式を導きます。

公式そのものはシンプルですが、導出には偏微分と全微分の違い、作用積分の変分、境界条件など、解析力学の重要な概念が登場します。2026年時点でも大学の力学教育で扱われる基礎的な枠組みであり、自由粒子や振り子の運動から、より複雑な力学系や場の理論まで幅広く応用されています。

📌 【この記事の重要ポイントまとめ】
  • 基本公式: d/dt(∂L/∂q̇)−∂L/∂q=0 が、オイラー=ラグランジュ方程式の標準形。
  • 導出の原理: ハミルトンの原理に基づき、作用積分の一次変分をゼロにすることで運動方程式を導く。
  • 具体的な応用: 自由粒子では等速直線運動、単振り子では角度に関する非線形微分方程式が得られる。

【基本から理解】オイラー=ラグランジュ方程式とは何か?

オイラー=ラグランジュ方程式とは、ある積分型の汎関数が停留値を取る条件から、その汎関数を定める関数が満たすべき微分方程式を導く公式です。解析力学では、これを利用して物体の運動方程式を求めます。

一般化座標を q、時間微分を q̇=dq/dt、ラグランジアンを L とすると、1つの一般化座標についての基本公式は次のとおりです。

d/dt(∂L/∂q̇)−∂L/∂q=0

複数の自由度を持つ系では、一般化座標 q₁、q₂、…、qₙのそれぞれについて、次の方程式を立てます。

d/dt(∂L/∂q̇ᵢ)−∂L/∂qᵢ=0 (i=1, 2, …, n)

ここで重要なのは、ラグランジアン L が座標 q、速度 q̇、時間 t の関数であることです。通常の古典力学では、運動エネルギーを T、ポテンシャルエネルギーを V とすると、

L=T−V

と定義します。ただし、これは代表的な力学系で使われる形であり、より一般的な理論では、ラグランジアンの構造が異なる場合もあります。

ニュートンの第二法則 F=ma では、物体に働く力を把握して運動を求めます。一方、ラグランジュ力学では、適切な座標系を選び、エネルギーから L を構成して方程式を導きます。拘束条件のある系や、角度などの一般化座標を使う系では、後者の方法によって計算を整理しやすくなります。

当時のメディア報道・掲載写真
【検証資料 1】当時のメディア報道・掲載写真(出典:i.ytimg.com)

【公式の意味】ラグランジアン・作用積分・変分法の基本原理

オイラー=ラグランジュ方程式を理解するには、ラグランジアン、作用積分、変分法という3つの概念を押さえる必要があります。

1.ラグランジアンの定義

ラグランジアン L は、力学系の状態を座標と速度によって記述する関数です。標準的な保存力系では L=T−V と表されます。たとえば、質量 m の粒子が直線上を動く場合、運動エネルギーは T=½m q̇² です。位置 q に依存するポテンシャルエネルギーを V(q) とすれば、L=½m q̇²−V(q) となります。

2.作用積分の定義

作用 S は、ラグランジアンを時間について積分した量です。

S[q]=∫ₜ₁ᵗ² L(q(t), q̇(t), t) dt

積分区間は時刻 t₁から t₂までです。作用は、時刻ごとの座標や速度の値だけでなく、その間の運動経路全体によって決まります。そのため、作用は通常の関数ではなく、関数を入力として数値を返す「汎関数」として扱われます。

SI単位では、ラグランジアンがジュール、時間が秒なので、作用の単位はジュール秒です。

3.変分法の基本原理とハミルトンの原理

変分法は、関数そのものを少し変化させたときに、積分などの汎関数がどのように変化するかを調べる数学的手法です。通常の微分法では変数を少し動かして関数の変化を調べますが、変分法では、ある関数を別の近接した関数へと変化させます。

古典力学におけるハミルトンの原理では、適切な条件のもとで、実際の運動経路は作用 S が停留値を取る経路として特徴づけられます。これを数式で表すと、

δS=0

となります。ここで δS は作用の一次変分です。「作用が必ず最小になる」という意味ではありません。停留値には極小値、極大値、鞍点などが含まれます。「最小作用の原理」という呼び名だけで、常に最小値を取ると解釈しないことが大切です。

カナダのゲルフ大学の公開物理教材や、オランダのデルフト工科大学(TU Delft)の公開教材でも、作用積分と境界条件を用いてオイラー=ラグランジュ方程式を導く方法が説明されています。MIT OpenCourseWareの力学講義資料でも、一般化座標や単振り子への適用が扱われています。これらに共通するのは、力を個別に列挙する代わりに、系のエネルギーと運動経路の性質から運動方程式を組み立てるという考え方です。

【導出方法を解説】作用積分の変分から公式が出るまで

オイラー=ラグランジュ方程式の導出では、実際の運動経路と、その近くにある仮想的な経路を比較します。ポイントは、作用の変化を計算したあと、部分積分によって時間微分を整理することです。

ステップ1:近接した経路を設定する

実際の経路を q(t) とし、比較する経路を次のように定めます。

qε(t)=q(t)+εη(t)

ε は十分に小さい実数、η(t) は経路の変化を表す滑らかな関数です。両端の時刻で位置を固定するため、境界条件を

η(t₁)=η(t₂)=0

とします。

ステップ2:作用積分を変化させる

変化させた経路に対する作用を S[qε] とします。ε で微分し、ε=0 で評価すると、作用の一次変分は次の形になります。

δS=∫ₜ₁ᵗ²{(∂L/∂q)η+(∂L/∂q̇)η̇}dt

これは、ラグランジアンが座標 q と速度 q̇ の両方に依存することから得られます。η̇ は η の時間微分です。

ステップ3:部分積分を使う

第2項の積分を部分積分すると、次のように変形できます。

∫ₜ₁ᵗ²(∂L/∂q̇)η̇ dt
=[(∂L/∂q̇)η]ₜ₁ᵗ²−∫ₜ₁ᵗ²{d/dt(∂L/∂q̇)}η dt

境界条件 η(t₁)=η(t₂)=0 により、最初の境界項はゼロになります。その結果、作用の一次変分は次の形に整理されます。

δS=∫ₜ₁ᵗ²{∂L/∂q−d/dt(∂L/∂q̇)}η dt

ステップ4:任意の変分について成立する条件を使う

ハミルトンの原理から δS=0 です。この条件が、端点でゼロになる任意の適切な変分 η に対して成立するためには、積分の被積分関数が区間の内部でゼロでなければなりません。したがって、

∂L/∂q−d/dt(∂L/∂q̇)=0

となり、整理すると、

d/dt(∂L/∂q̇)−∂L/∂q=0

が得られます。これがオイラー=ラグランジュ方程式です。

この導出で見落としやすいのが、部分積分と境界条件の役割です。作用を変分するだけでは、公式の形はまだ完成しません。境界項を消去し、変分関数が区間の内部で任意に選べることを利用することで、局所的な微分方程式が得られます。

活動歴および当時の関連ビジュアル記録
【検証資料 2】活動歴および当時の関連ビジュアル記録(出典:st-note.com)

【具体例で理解】自由粒子と単振り子の解き方

公式の意味をつかむには、ラグランジアンを定義し、偏微分を計算し、時間について全微分を取るという手順を具体例で確認するのが効果的です。

例1:自由粒子の運動

質量 m の粒子が一次元を運動し、外力もポテンシャルエネルギーもないとします。座標を x とすると、ラグランジアンは運動エネルギーだけで構成されます。

L=½m ẋ²

まず、速度について偏微分します。

∂L/∂ẋ=m ẋ

次に、位置 x について偏微分します。L は x に直接依存しないため、

∂L/∂x=0

です。これらをオイラー=ラグランジュ方程式に代入すると、

d/dt(m ẋ)=0

質量 m が一定なので、

m ẍ=0、すなわち ẍ=0

が得られます。加速度がゼロであるため、粒子は静止するか、一定の速度で直線運動を続けます。ニュートンの第一法則および第二法則と一致する結果です。

例2:単振り子の運動

長さ ℓ、質量 m の振り子を考えます。糸は伸びず、空気抵抗はなく、重力加速度を g とします。角度 θ は鉛直下向きの平衡位置から測るものとします。

運動エネルギーとポテンシャルエネルギーは、

T=½mℓ²θ̇²

V=mgℓ(1−cosθ)

と表せます。したがって、ラグランジアンは、

L=½mℓ²θ̇²−mgℓ(1−cosθ)

です。角速度 θ̇ について偏微分すると、

∂L/∂θ̇=mℓ²θ̇

これを時間について全微分すると、

d/dt(∂L/∂θ̇)=mℓ²θ̈

一方、角度 θ について偏微分すると、

∂L/∂θ=−mgℓ sinθ

となります。これらを基本公式に代入して整理すると、

mℓ²θ̈+mgℓ sinθ=0

両辺を mℓ² で割れば、

θ̈+(g/ℓ)sinθ=0

という振り子の運動方程式が得られます。振り子の糸に働く張力を未知の力として個別に求めなくても、一般化座標 θ とエネルギーから運動方程式を導ける点が、ラグランジュ力学の利点です。

さらに、振れ角が十分に小さい場合には sinθ≈θ と近似できます。そのとき、

θ̈+(g/ℓ)θ=0

となり、単振動の微分方程式に変わります。角振動数は ω=√(g/ℓ)、周期は T₀=2π√(ℓ/g)です。この周期の式は、小振幅近似の範囲で成立します。

2つの例を比べると、公式を適用する流れが共通していることが分かります。

項目自由粒子単振り子
一般化座標xθ
ラグランジアン½m ẋ²½mℓ²θ̇²−mgℓ(1−cosθ)
運動方程式ẍ=0θ̈+(g/ℓ)sinθ=0
特徴加速度がゼロ重力による復元運動

実際に解くときの手順

問題演習では、次の順番で計算すると、微分の取り違えや符号ミスを減らせます。

  1. 運動を表す一般化座標 q を選ぶ。
  2. 運動エネルギー T とポテンシャルエネルギー V を求める。
  3. ラグランジアン L=T−V を作る。
  4. ∂L/∂q̇ と ∂L/∂q を別々に計算する。
  5. ∂L/∂q̇ を時間について全微分する。
  6. オイラー=ラグランジュ方程式に代入し、運動方程式を整理する。
  7. 初期条件や境界条件を使い、必要に応じて微分方程式を解く。

【成立条件と境界条件】公式を適用できる範囲と注意点

オイラー=ラグランジュ方程式は広く使える公式ですが、どのような問題にも同じ形を無条件で適用できるわけではありません。導出の前提と、得られた方程式を解くための条件を区別する必要があります。

1.ラグランジアンの滑らかさ

標準的な導出では、L が座標や速度について必要な偏微分を持ち、時間についても適切に微分できることを仮定します。変分関数も十分に滑らかで、部分積分が正当化できる必要があります。関数に不連続性や特異点がある場合には、区間を分割するなど、問題に応じた追加の取り扱いが必要になることがあります。

2.境界条件と固定端点

先ほどの導出では、始点と終点の位置を固定し、η(t₁)=η(t₂)=0 としました。この条件によって部分積分の境界項が消えます。

ただし、境界条件が常に固定端点とは限りません。端点の位置が自由な問題では、境界項をゼロとするための別の条件が必要になることがあります。場の理論や最適化問題でも、境界条件の設定によって得られる解が変わります。

3.停留作用と解の一意性は別の問題

δS=0 は、実際の運動が作用の停留条件を満たすことを表します。しかし、この条件だけで解の一意性や安定性、作用の最小性が自動的に保証されるわけではありません。解の存在や一意性を調べるには、ラグランジアンの性質や境界条件などをさらに検討する必要があります。

4.拘束条件と非保存力

糸の長さが一定の振り子のような系では、適切な一般化座標を選ぶことで拘束条件を座標の定義に組み込めます。一方、一般化座標だけでは表現しにくい拘束条件がある場合には、ラグランジュの未定乗数などを使う方法があります。

また、空気抵抗や摩擦のような非保存力が存在する場合、単純な L=T−V だけでは必要な力学を完全に表現できないことがあります。その場合は、一般化力を右辺に加える方法や、散逸関数などを導入する方法を検討します。

さらに、ラグランジアンの速度に関する二階の偏微分を並べた行列が特異な場合には、通常の形の運動方程式から加速度を一意に求められないことがあります。このような系では、拘束や制約を含めた解析が必要です。

公の場での発言・インタビュー報道記録
【検証資料 3】公の場での発言・インタビュー報道記録(出典:i.redd.it)

【よくある誤解を整理】偏微分と全微分、解析力学の使い方

オイラー=ラグランジュ方程式の計算で特に間違えやすいのが、偏微分と全微分の使い分けです。公式の左辺には2種類の微分が登場しますが、それぞれの意味は異なります。

偏微分 ∂L/∂q と ∂L/∂q̇

偏微分は、複数の変数を持つ関数について、ほかの変数を固定して特定の変数だけを変化させたときの変化率です。L=½m q̇²−V(q) の場合、速度について偏微分すると ∂L/∂q̇=m q̇、座標について偏微分すると ∂L/∂q=−dV/dq となります。

全微分 d/dt(∂L/∂q̇)

一方、d/dt は、時間の経過に伴う変化を表します。∂L/∂q̇ は q や q̇、場合によっては t にも依存するため、その時間微分には連鎖律が関係します。

たとえば、L=½m q̇²−V(q) なら、

d/dt(∂L/∂q̇)=d/dt(m q̇)=m q̈

です。速度についての偏微分を計算したあと、さらに時間について微分することを忘れないようにしましょう。

よくある誤解

  • 作用は必ず最小になる:正しくは、基本原理が要求するのは停留条件です。
  • L はいつでも T−V と書ける:標準的な力学では有用ですが、一般の系では追加の項や別の構造が必要になることがあります。
  • 公式を使えば運動の解まで完成する:公式から得られるのは運動方程式であり、具体的な軌道を求めるには初期条件や境界条件、微分方程式の解法が必要です。
  • 偏微分と全微分は同じ:前者は変数ごとの依存性を調べ、後者は時間に沿った変化を計算します。

解析力学の基礎を学ぶ段階では、まず自由粒子、ばね振動、単振り子のような例で計算手順を身につけるのが有効です。その後、一般化座標、拘束条件、対称性、保存則へと学習を進めると、ラグランジュ力学の利点が見えやすくなります。

特に、座標の選択によって計算量が大きく変わる問題では、この方法の価値が際立ちます。直交座標では複雑な拘束条件も、角度や距離などの一般化座標を使えば簡潔に表現できる場合があります。さらに、ラグランジアンの対称性から運動量やエネルギーの保存則を導く考え方は、解析力学を深く理解するうえで重要です。

【オイラー=ラグランジュ方程式】に関するよくある質問(FAQ)

Q1:オイラー=ラグランジュ方程式の公式を簡単に覚える方法は?
A1:d/dt(∂L/∂q̇)−∂L/∂q=0 と覚えましょう。まず速度 q̇ について偏微分し、その結果を時間について全微分してから、座標 q についての偏微分を引きます。複数の一般化座標がある場合は、座標ごとに同じ計算を繰り返します。

Q2:オイラー=ラグランジュ方程式はニュートンの運動方程式と何が違う?
A2:ニュートンの運動方程式は力と加速度の関係から運動を求めます。オイラー=ラグランジュ方程式は、ラグランジアンと作用の停留条件から運動方程式を導きます。適切な条件のもとで両者は同じ物理的な運動を記述し、拘束条件や一般化座標がある問題ではラグランジュ力学のほうが計算を整理しやすいことがあります。

Q3:オイラー=ラグランジュ方程式を導出する際、なぜ境界条件が必要なの?
A3:標準的な導出では、始点と終点の位置を固定し、変分関数が両端でゼロになると仮定します。これによって部分積分で現れる境界項が消え、積分の内部にある式から運動方程式を導けます。端点が自由な場合には、別の境界条件を考慮する必要があります。

【プロの結論】まとめ:公式の暗記から導出と応用の理解へ

オイラー=ラグランジュ方程式は、作用積分の変分とハミルトンの原理を通じて、力学系の運動方程式を導くための基本的な道具です。公式は d/dt(∂L/∂q̇)−∂L/∂q=0 という簡潔な形ですが、その背後にはラグランジアンの定義、変分法、部分積分、境界条件という一連の論理があります。

理解を確実にするには、公式だけを暗記するのではなく、自由粒子や単振り子を使って、ラグランジアンの設定から運動方程式が得られるまでを自分で計算することが重要です。その際、偏微分と時間についての全微分を区別し、作用の停留条件と境界条件の関係を意識してください。

さらに、適用条件を把握しておけば、保存力系だけでなく、拘束条件のある力学系や、より高度な解析力学の問題にも対応しやすくなります。公式を正しく使うことに加えて、なぜその式が成立するのか、どの条件のもとで有効なのかを説明できることが、解析力学を身につけるうえでの確かな判断基準です。 (出典: オイラー ラグランジュ 方程式(Yahoo!ニュース))

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