導入 摩擦のないワイヤー上を移動するビーズ。ワイヤーはビーズに反作用力C を及ぼし、ビーズをワイヤー上に留めている。この場合の非拘束力N は重力である。ワイヤー上のビーズの初期位置によって、異なる動きが生じることに注目してほしい。 単振り子。棒は剛体であるため、おもりの位置は方程式f ( x , y ) = 0 に従って拘束され、拘束力C は棒の張力です。この場合も、非拘束力Nは重力です。 ニュートンの法則と力の概念は、機械システムについて教える際の一般的な出発点です。[ 5 ] この方法は多くの問題には有効ですが、他の問題ではアプローチが悪夢のように複雑になります。[ 6 ] 例えば、真珠が内部で滑っている水平面上を転がるトーラスの運動を計算する場合、トーラスの 角速度 やトーラスに対する真珠の運動などの時間的に変化する拘束力により、ニュートンの方程式でトーラスの運動を決定することが困難になります。[ 7 ] ラグランジュ力学は、力ではなくエネルギーを基本要素として採用し、[ 5 ] より複雑な問題に取り組むことができる、より抽象的な方程式につながります。[ 6 ]
特に、ラグランジュのアプローチは、各物体の位置と速度に対して独立した一般化座標 を設定することであり、これによりラグランジアン(系の全運動エネルギーから位置エネルギーを引いたもの)の一般形を書き下すことが可能になり、粒子のすべての可能な運動経路にわたってこれを合計すると「作用」の式が得られ、それを最小化して一般化された方程式のセットが得られます。この合計量は、粒子が実際にたどる経路に沿って最小化されます。この選択により、拘束力が結果として得られる一般化された方程式系 に入る必要がなくなります。特定の瞬間に粒子に対する拘束の影響を直接計算しないため、方程式の数が少なくなります。[ 7 ]
さまざまな物理システムにおいて、質量のある物体の大きさや形状が無視できる場合、それを点粒子として扱うことは有用な簡略化です。 質量 m 1 、m 2 、 ...、m N のN 個の点粒子のシステムでは、各粒子は位置ベクトル r 1 、r 2 、 ...、r N を持ちます。デカルト座標で 十分な場合が多く、r 1 = ( x 1 、 y 1 、 z 1 ) 、r 2 = ( x 2 、 y 2 、 z 2 ) などとなります。3次元空間 では、各位置ベクトルは点の位置を一意に定義するために3 つの座標を必要とするため、システムの構成を一意に定義するために 3 N 個の座標があります。これらはすべて粒子の位置を特定するための空間内の特定の点です。空間内の一般的な点は r = ( x 、 y 、 z ) と書きます。各粒子の速度は 、粒子が運動経路に沿って移動する速さであり、位置の 時間微分であるため、 v 1 = d r 1 d t 、 v 2 = d r 2 d t 、 … 、 v N = d r N d t 。 {\displaystyle \mathbf {v} _{1}={\frac {d\mathbf {r} _{1}}{dt}},\mathbf {v} _{2}={\frac {d\mathbf {r} _{2}}{dt}},\ldots ,\mathbf {v} _{N}={\frac {d\mathbf {r} _{N}}{dt}}.} ニュートン力学では、運動方程式は ニュートンの法則 によって与えられます。第二法則「正味の力は 質量×加速度 に等しい」、 ∑ F = m d 2 r d t 2 、 {\displaystyle \sum \mathbf {F} =m{\frac {d^{2}\mathbf {r} }{dt^{2}}},} これは各粒子に適用されます。3次元のN 粒子系の場合、粒子の位置に関する3N個の2次常微分方程式を解く必要があり ます。
ラグランジアン ラグランジュ力学では、力の代わりにシステムのエネルギーを使用します。ラグランジュ力学の中心となる量は ラグランジアン であり、これはシステム全体のダイナミクスを要約する関数です。一般的に、ラグランジアンはエネルギーの単位を持ちますが、すべての物理システムに対して単一の表現はありません。物理法則に合致する正しい運動方程式を生成する関数は、ラグランジアンとして扱うことができます。ただし、幅広い応用分野に対して一般的な表現を構築することは可能です。電磁場のない粒子システムの非相対論的ラグランジアンは、 [ 8 ]で与えられます。 L = T − V 、 {\displaystyle L=TV,} どこ T = 1 2 ∑ k = 1 N m k v k 2 {\displaystyle T={\frac {1}{2}}\sum _{k=1}^{N}m_{k}v_{k}^{2}} はシステムの 全運動エネルギー であり、の運動エネルギーの合計Σ に等しい。 N {\displaystyle N} 粒子。各粒子にはラベルが付けられています。k {\displaystyle k} 質量を持つm k 、 {\displaystyle m_{k},} また、v k 2 = v k · v k は速度の大きさの二乗であり、速度とそれ自身との内積に等しい。 [ 9 ]
運動エネルギーT はシステムの運動のエネルギーであり、速度v k のみの関数であり、位置r k や時間t の関数ではないため、T = T ( v 1 , v 2 , ...) となります。
系のポテンシャルエネルギー V は 、粒子間の相互作用エネルギー、つまり、外部からの影響と合わせて、ある粒子が他のすべての粒子からどれだけのエネルギーを得ているかを反映しています。保存力( ニュートン重力 など) の場合、V は粒子の位置ベクトルのみの関数なので、V = V ( r 1 , r 2 , ...) となります。適切なポテンシャル (電磁ポテンシャル など) から導出できる非保存力の場合、速度も現れ、V = V ( r 1 , r 2 , ..., v 1 , v 2 , ...) となります。外部場や外部駆動力が時間とともに変化する場合、ポテンシャルも時間とともに変化するので、最も一般的にはV = V ( r 1 , r 2 , ..., v 1 , v 2 , ..., t ) となります。
数式による定式化(有限粒子系の場合)ラグランジアンの同等かつより数学的に厳密な定義は以下のとおりである。[ 10 ] 3次元空間におけるN 個の粒子 の系の場合、その系の構成空間は滑らかな多様体である。 Q ⊆ R 3 N \displaystyle Q\subseteq \mathbb {R} ^{3N}} 各構成q ∈ Q {\displaystyle q\in Q} これは、特定の時点における各粒子の空間的位置を指定するものであり、多様体は、システムに対する制約によって許容されるすべての構成から構成される。
ラグランジアンは滑らかな関数です。 L : T Q × R → R 、 {\displaystyle L:TQ\times \mathbb {R} \to \mathbb {R} ,} どこT Q ⊆ R 6 N {\displaystyle TQ\subseteq \mathbb {R} ^{6N}} は構成空間の接束 です。つまり、の各要素はT Q {\displaystyle TQ} は粒子の位置と速度の両方を表し、タプルとして記述できます。( q 1 、 … 、 q N 、 q ˙ 1 、 … q ˙ N ) {\displaystyle (q_{1},\ldots ,q_{N},{\dot {q}}_{1},\ldots {\dot {q}}_{N})} とq 私 q_i そしてq ˙ 私 \displaystyle {\dot {q}}_{i}} i 番目の粒子の位置と速度をそれぞれ指定する。時間依存性t ∈ R \displaystyle t\in \mathbb {R} } これにより、ラグランジアンは時間依存の力やポテンシャルを記述することができる。
システムの軌跡は 滑らかな関数である q : [ t 0 、 t 1 ] → Q 、 {\displaystyle q:[t_{0},t_{1}]\to Q,} 時間の経過に伴う構成の変化を記述する。その速度は時間微分である。q ˙ ( t ) ∈ T q ( t ) Q {\displaystyle {\dot {q}}(t)\in T_{q(t)}Q} のq ( t ) {\displaystyle q(t)} 、そしてそのペア( q ( t ) 、 q ˙ ( t ) ) ∈ T Q {\displaystyle (q(t),{\dot {q}}(t))\in TQ} したがって、は任意のバンドルの要素であるt {\displaystyle t} したがって、軌道の 作用汎関数は、経路に沿ったラグランジアンの積分として定義できる。 S [ q ] = ∫ t 0 t 1 L ( q ( t ) 、 q ˙ ( t ) 、 t ) d t 。 {\displaystyle S[q]=\int _{t_{0}}^{t_{1}}L(q(t),{\dot {q}}(t),t)\,dt.}
ラグランジュ力学における運動法則は、与えられた2つの配置間のすべての軌道のうち、実際に系がたどる軌道は作用汎関数の臨界点(多くの場合、局所最小値であるが、必ずしもそうとは限らない)でなければならないという仮定から導き出される。これにより、オイラー・ラグランジュ方程式 が得られる(下記も参照)。
運動方程式 質量を持つ粒子のシステムの場合m 1 、 … 、 m N {\displaystyle m_{1},\ldots ,m_{N}} 運動エネルギー は次のようになります。 T ( q 、 q ˙ ) = 1 2 ∑ 私 = 1 N m 私 ‖ q ˙ 私 ‖ 2 、 {\displaystyle T(q,{\dot {q}})={\frac {1}{2}}\sum _{i=1}^{N}m_{i}\|{\dot {q}}_{i}\|^{2},} どこq ˙ 私 ∈ R 3 \displaystyle {\dot {q}}_{i}\in \mathbb {R} ^{3}} は粒子i の速度です。
位置エネルギー V ( q 、 t ) {\displaystyle V(q,t)} これは構成(そして場合によっては時間)のみに依存し、通常は保存力から生じる。
標準ラグランジアン は、以下の差分で与えられる。 L ( q 、 q ˙ 、 t ) = T ( q 、 q ˙ ) − V ( q 、 t ) 。 {\displaystyle L(q,{\dot {q}},t)=T(q,{\dot {q}})-V(q,t).}
この定式化は、保存系と時間依存系の両方を対象としており、連続系(場)、制約系、および曲面構成空間を持つ系への一般化の基礎となる。
時間変動制約や外部影響により、T またはV 、あるいはその両方が時間に対して明示的に依存する場合、ラグランジアンL ( r 1 , r 2 , ... v 1 , v 2 , ... t )は時間依存と なります。ポテンシャルエネルギーも運動エネルギーも時間に依存しない場合、ラグランジアンL ( r 1 , r 2 , ... v 1 , v 2 , ...)は時間に対して明示的に独立 となります。いずれの場合も、ラグランジアンは一般化座標を通して暗黙的に時間依存性を持ちます。
これらの定義を用いると、ラグランジュ方程式 は次のようになる[ 11 ]
ラグランジュ方程式 ∂ L ∂ r k − d d t ∂ L ∂ r ˙ k + ∑ 私 = 1 C λ 私 ∂ f 私 ∂ r k = 0 、 {\displaystyle {\frac {\partial L}{\partial \mathbf {r} _{k}}}-{\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L}{\partial {\dot {\mathbf {r} }}_{k}}}+\sum _{i=1}^{C}\lambda _{i}{\frac {\partial f_{i}}{\partial \mathbf {r} _{k}}}=0,}
ここで、k = 1, 2, ..., N は 粒子をラベル付けし、各制約方程式f i にはラグランジュ乗数 λ i があり、 ∂ ∂ r k ≡ ( ∂ ∂ x k 、 ∂ ∂ y k 、 ∂ ∂ z k ) 、 ∂ ∂ r ˙ k ≡ ( ∂ ∂ x ˙ k 、 ∂ ∂ y ˙ k 、 ∂ ∂ z ˙ k ) {\displaystyle {\frac {\partial }{\partial \mathbf {r} _{k}}}\equiv \left({\frac {\partial }{\partial x_{k}}},{\frac {\partial }{\partial y_{k}}},{\frac {\partial }{\partial z_{k}}}\right),\quad {\frac {\partial }{\partial {\dot {\mathbf {r} }}_{k}}}\equiv \left({\frac {\partial }{\partial {\dot {x}}_{k}}},{\frac {\partial }{\partial {\dot {y}}_{k}}},{\frac {\partial }{\partial {\dot {z}}_{k}}}\right)} はそれぞれ、指定された変数に関する偏微分ベクトル ∂/∂ の略記です(ベクトル全体に関する微分ではありません)。 [ nb 1 ] 各上付きドットは時間微分 の略記です。この手順では、位置座標と乗数に 3 N 個の連立 2 階微分方程式とC 個の 拘束方程式があるため、ニュートンの法則と比較して解く方程式の数が 3 N から3 N + C に増加します。ただし、粒子の位置座標とともに解くと、乗数から拘束力に関する情報を得ることができます。拘束方程式を解くことで座標を消去する必要はありません。
ラグランジアンでは、位置座標と速度成分はすべて独立変数 であり、ラグランジアンの導関数は通常の微分規則 に従ってこれらに関して個別に取られます(たとえば、粒子2のz 速度成分v z ,2 = dz 2 / dt で定義されるL の偏微分は、単に∂ L /∂ v z ,2 です。速度成分を対応する座標z 2に関連付けるために 、 面倒な連鎖規則 や全微分を使用する必要はありません)。
各制約方程式において、1つの座標は他の座標から決定されるため冗長です。したがって、独立座標の数は n = 3 N − C となります。各位置ベクトル、ひいては位置座標を一般化座標と時間の関数として表すことで、各位置ベクトルを共通の n 個の 一般化座標 のセットに変換できます。これは、n タプルq = ( q 1 , q 2 , ... q n ) として表すと便利です。 r k = r k ( q 、 t ) = ( x k ( q 、 t ) 、 y k ( q 、 t ) 、 z k ( q 、 t ) 、 t ) 。 {\displaystyle \mathbf {r} _{k}=\mathbf {r} _{k}(\mathbf {q} ,t)={\big (}x_{k}(\mathbf {q} ,t),y_{k}(\mathbf {q} ,t),z_{k}(\mathbf {q} ,t),t{\big )}.}
ベクトルq はシステムの構成空間 内の点です。一般化座標の時間微分は一般化速度と呼ばれ、各粒子について、その速度ベクトルの変換、つまり時間に関する位置の 全微分は次のようになります。 q ˙ j = d q j d t 、 v k = ∑ j = 1 n ∂ r k ∂ q j q ˙ j + ∂ r k ∂ t 。 {\displaystyle {\dot {q}}_{j}={\frac {\mathrm {d} q_{j}}{\mathrm {d} t}},\quad \mathbf {v} _{k}=\sum _{j=1}^{n}{\frac {\partial \mathbf {r} _{k}}{\partial q_{j}}}{\dot {q}}_{j}+{\frac {\partial \mathbf {r} _{k}}{\partial t}}.}
このv k を考慮すると、一般化座標における 運動エネルギーは、一般化速度、一般化座標、および時間に依存します。位置ベクトルが時間変動制約により明示的に時間に依存する場合、T = T ( q 、 q ˙ 、 t ) 。 {\displaystyle T=T(\mathbf {q} ,{\dot {\mathbf {q} }},t).}
これらの定義により、オイラー・ラグランジュ方程式 、[ 12 ] [ 13 ]
ラグランジュ方程式 d d t ( ∂ L ∂ q ˙ j ) = ∂ L ∂ q j {\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial L}{\partial {\dot {q}}_{j}}}\right)={\frac {\partial L}{\partial q_{j}}}}
これらは変分法の 数学的結果であり、力学でも使用できます。ラグランジアンL ( q , d q /d t , t ) に代入すると、システムの運動方程式 が得られます。方程式の数は、ニュートン力学と比較して、3 N からn = 3 N − C 個の一般化座標における連立2階微分方程式に減少しています。これらの方程式には拘束力は全く含まれておらず、非拘束力のみを考慮すれば十分です。
運動方程式には偏微分が 含まれるが、偏微分の結果は依然として粒子の位置座標に関する常微分方程式である。 全時間微分 d/d tは、多くの場合、 陰関数微分 を伴う。どちらの方程式もラグランジアンに関しては線形であるが、一般的には座標に関しては非線形連立方程式である。
拡張機能 既に述べたように、この形式のL は多くの重要なクラスのシステムに適用できますが、すべての場合に適用できるわけではありません。相対論的ラグランジュ力学 の場合、特殊相対性理論 (ローレンツ変換の下でのスカラー) または一般相対性理論 (4 スカラー) と整合する関数で全体を置き換える必要があります。[ 14 ] 磁場が存在する場合、ポテンシャルエネルギーの式を再記述する必要があります。また、散逸力 (摩擦など) の場合、エネルギー損失を考慮するために、ラグランジュ関数に加えて、しばしば「 レイリー散逸関数 」と呼ばれる別の関数を導入する必要があります。[ 15 ]
粒子の 1 つ以上は、それぞれ 1 つ以上のホロノミック拘束条件を受ける可能性があります。このような拘束条件は、 f ( r , t ) = 0 の形式の方程式で表されます。システム内の拘束条件の数がC である場合、各拘束条件は、 f 1 ( r , t ) = 0、f 2 ( r , t ) = 0、 ...、f C ( r , t ) = 0 という方程式を持ち、 それぞれが任意の粒子に適用できます。粒子kが拘束条件 i を受ける場合、f i ( r k , t ) = 0 と なります。任意の瞬間において、拘束された粒子の座標は互いにリンクしており、独立ではありません。拘束方程式は、粒子が移動できる許容経路を決定しますが、各瞬間における粒子の位置や速度は決定しません。非ホロノミック拘束条件は、 粒子の速度、加速度、または位置の高階微分に依存します。ラグランジュ力学は、制約条件が存在する場合、それがすべてホロノミックであるシステムにのみ適用できます 。非ホロノミック制約の3つの例は次のとおりです。[ 16 ] 制約方程式が非積分可能な場合、制約条件に不等式が含まれる場合、または制約条件に摩擦などの複雑な非保存力が含まれる場合。非ホロノミック制約には特別な処理が必要であり、ニュートン力学 に戻るか、他の方法を使用する必要がある場合があります。[ 17 ]
ニュートン力学からラグランジュ力学へ
ニュートンの法則アイザック・ニュートン (1642年 - 1727年)簡略化のため、ニュートンの法則は一般性をほとんど損なうことなく、1つの粒子について説明できます(N 個の粒子からなる系の場合、これらの式はすべて系内の各粒子に適用されます)。質量mが一定の粒子の 運動方程式は、1687年の ニュートンの第2法則 であり、現代のベクトル表記では 次のようになります。F = m 1 、 {\displaystyle \mathbf {F} =m\mathbf {a} ,} ここで、a はその加速度、F はその物体 に 作用する合力です。質量が変化する場合には、運動量の時間微分を取るように方程式を一般化する必要があります。3 次元空間では、このベクトル方程式には 3 つの成分があるため、これは 3 つの連立 2 階常微分方程式のシステムであり、解く必要があります。解は、 t = 0のときのr とv の初期条件を満たす 、時刻t における粒子の位置ベクトルrです 。
ニュートンの法則はデカルト座標では使いやすいが、デカルト座標は必ずしも便利ではなく、他の座標系では運動方程式が複雑になることがある。曲線座標 ξ = ( ξ 1 , ξ 2 , ξ 3 ) の法則は、 テンソルインデックス表記では 「ラグランジュ形式」 [ 18 ] [ 19 ] となる。F 1 = m ( d 2 ξ 1 d t 2 + Γ 1 b c d ξ b d t d ξ c d t ) = g 1 k ( d d t ∂ T ∂ ξ ˙ k − ∂ T ∂ ξ k ) 、 ξ ˙ 1 ≡ d ξ 1 d t 、 {\displaystyle F^{a}=m\left({\frac {\mathrm {d} ^{2}\xi ^{a}}{\mathrm {d} t^{2}}}+\Gamma ^{a}{}_{bc}{\frac {\mathrm {d} \xi ^{b}}{\mathrm {d} t}}{\frac {\mathrm {d} \xi ^{c}}{\mathrm {d} t}}\right)=g^{ak}\left({\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial T}{\partial {\dot {\xi }}^{k}}}-{\frac {\partial T}{\partial \xi ^{k}}}\right),\quad {\dot {\xi }}^{a}\equiv {\frac {\mathrm {d} \xi ^{a}}{\mathrm {d} t}},} ここで、F a は 粒子に作用する合力のa 番目の反変成分 であり、Γ a bc は第 2 種の クリストッフェル記号 である。T = 1 2 m g b c d ξ b d t d ξ c d t {\displaystyle T={\frac {1}{2}}mg_{bc}{\frac {\mathrm {d} \xi ^{b}}{\mathrm {d} t}}{\frac {\mathrm {d} \xi ^{c}}{\mathrm {d} t}}} は粒子の運動エネルギーであり、g bc は曲線座標系の 計量テンソル の共変成分 です。添え字a 、b 、c はそれぞれ 1、2、3 の値をとります。曲線座標は一般化座標とは異なります。
ニュートンの法則をこのような形で表現するのは複雑すぎるように思えるかもしれませんが、利点があります。クリストッフェル記号で表される加速度成分は、運動エネルギーの導関数を評価することで回避できます。粒子に作用する合力がゼロの場合、F = 0 であれば、 粒子は加速せず、一定の速度で直線上を移動します。数学的には、微分方程式の解は測地線 、つまり空間内の2点間の極値の長さの曲線です(これらは最短経路、つまり最小になる場合もありますが、必ずしもそうとは限りません)。平坦な3次元実空間では、測地線は単なる直線です。したがって、自由粒子の場合、ニュートンの第2法則は測地線方程式と一致し、自由粒子は移動できる極値軌道である測地線をたどると述べています。粒子に力F ≠ 0 が作用する場合、 粒子は作用する力によって加速し、自由であればたどるはずの測地線から逸れます。ここで平坦な 3D 空間で与えられた量を適切に 4D曲がった時空 に拡張すると、上記のニュートンの法則の形式はアインシュタイン の一般相対性理論 にも引き継がれ、その場合、自由粒子は曲がった時空の中で、もはや通常の意味での「直線」ではない測地線に沿って移動します。[ 20 ]
しかし、粒子に作用する合計合力F を知る必要があり、そのためには合力非拘束力N と合力拘束力C が必要となります。 F = C + N 。 {\displaystyle \mathbf {F} =\mathbf {C} +\mathbf {N} .}
拘束力は一般的に時間に依存するため、複雑になる可能性がある。また、拘束が存在する場合、曲線座標は独立ではなく、1つ以上の拘束方程式によって関連付けられる。
拘束力は、運動方程式から消去して非拘束力のみを残すか、拘束方程式を運動方程式に含めることによって含めることができる。
ダランベールの原理ジャン・ダランベール (1717年 - 1783年)曲線上に閉じ込められた質量m の粒子に対する拘束力C と仮想変位δ r 。結果として生じる非拘束力はN である。
解析力学 における基本的な結果の一つは、静的平衡を 理解するために1708年にジャック・ベルヌーイ によって導入され、 1743年にダランベールによって動的問題を解決するために発展させられた ダランベールの原理 である。[ 21 ] この原理は、N 個の粒子に対して仮想仕事、すなわち仮想変位に沿った仕事が、δ r k {\displaystyle \delta \mathbf {r} _{k}} 、はゼロです: [ 9 ] ∑ k = 1 N ( N k + C k − m k 1 k ) ⋅ δ r k = 0. {\displaystyle \sum _{k=1}^{N}(\mathbf {N} _{k}+\mathbf {C} _{k}-m_{k}\mathbf {a} _{k})\cdot \delta \mathbf {r} _{k}=0.}
仮想変位 、δ r k {\displaystyle \delta \mathbf {r} _{k}} は、定義上、ある瞬間に システムに作用する拘束力と一致するシステム構成の微小変化であり、[ 22 ] つまり、拘束力が拘束された運動を維持するような変化です。これらは、粒子を加速および移動させるために粒子に作用する結果として生じる拘束力と非拘束力によって引き起こされる、システム内の実際の変位とは異なります。[ 注2 ] 仮想仕事は 、任意の力(拘束力または非拘束力)に対して仮想変位に沿って行われる仕事です。
拘束力はシステム内の各粒子の運動に垂直に作用して拘束を維持するため、システムに作用する拘束力による仮想仕事の総量はゼロです。[ 23 ] [ nb3 ] ∑ k = 1 N C k ⋅ δ r k = 0 、 {\displaystyle \sum _{k=1}^{N}\mathbf {C} _{k}\cdot \delta \mathbf {r} _{k}=0,} となることによって ∑ k = 1 N ( N k − m k 1 k ) ⋅ δ r k = 0. {\displaystyle \sum _{k=1}^{N}(\mathbf {N} _{k}-m_{k}\mathbf {a} _{k})\cdot \delta \mathbf {r} _{k}=0.}
このように、ダランベールの原理により、運動方程式から拘束力を排除し、適用された非拘束力のみに集中することができます。[ 23 ] [ 24 ] 示された形式は、座標の選択にも依存しません。ただし、変位が任意の座標系で運動方程式を設定するために容易に使用することはできません。δ r k {\displaystyle \delta \mathbf {r} _{k}} 制約方程式によってN 個の項が結び付けられている場合があり、その場合、N個の個々の項を0に設定することはできません。そこで、個々の項が0の場合に限り総和が0となるような、互いに独立した座標系を探します。各項を0に設定することで、最終的に運動方程式を分離することができます。
ダランベールの原理に基づく運動方程式粒子k に制約がある場合、位置の座標r k = ( x k , y k , z k ) は制約方程式によって結び付けられているため、仮想変位 δ r k = ( δx k , δy k , δz k )も同様に結び付けられます。一般化座標は独立しているため、一般化座標の仮想変位に変換することでδ r k の複雑さを回避できます。これらは全微分 と同じ形式で関連付けられています。[ 9 ] δ r k = ∑ j = 1 n ∂ r k ∂ q j δ q j 。 {\displaystyle \delta \mathbf {r} _{k}=\sum _{j=1}^{n}{\frac {\partial \mathbf {r} _{k}}{\partial q_{j}}}\delta q_{j}.}
これは仮想的な変位、つまりある 瞬間の拘束条件に沿った変位であるため、時間に関する偏微分に時間増分を乗じたものは存在しない。
上記のダランベールの原理における最初の項は、仮想変位δ r k に沿って非拘束力N k が行う仮想仕事であり、一般性を失うことなく、一般化された力 の定義によって一般化された類似物に変換することができる。Q j = ∑ k = 1 N N k ⋅ ∂ r k ∂ q j 、 {\displaystyle Q_{j}=\sum _{k=1}^{N}\mathbf {N} _{k}\cdot {\frac {\partial \mathbf {r} _{k}}{\partial q_{j}}},} となることによって ∑ k = 1 N N k ⋅ δ r k = ∑ k = 1 N N k ⋅ ∑ j = 1 n ∂ r k ∂ q j δ q j = ∑ j = 1 n Q j δ q j 。 {\displaystyle \sum _{k=1}^{N}\mathbf {N} _{k}\cdot \delta \mathbf {r} _{k}=\sum _{k=1}^{N}\mathbf {N} _{k}\cdot \sum _{j=1}^{n}{\frac {\partial \mathbf {r} _{k}}{\partial q_{j}}}\delta q_{j}=\sum _{j=1}^{n}Q_{j}\delta q_{j}.}
これは一般化座標への変換の半分です。残りは加速度項を一般化座標に変換することですが、これはすぐには明らかではありません。ニュートンの第2法則のラグランジュ形式を思い出すと、運動エネルギーの一般化座標と速度に関する偏微分を求めることで、目的の結果が得られます。[ 9 ] ∑ k = 1 N m k 1 k ⋅ ∂ r k ∂ q j = d d t ∂ T ∂ q ˙ j − ∂ T ∂ q j 。 {\displaystyle \sum _{k=1}^{N}m_{k}\mathbf {a} _{k}\cdot {\frac {\partial \mathbf {r} _{k}}{\partial q_{j}}}={\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial T}{\partial {\dot {q}}_{j}}}-{\frac {\partial T}{\partial q_{j}}}.}
さて、ダランベールの原理は、要求どおり一般化座標で表され、 ∑ j = 1 n [ Q j − ( d d t ∂ T ∂ q ˙ j − ∂ T ∂ q j ) ] δ q j = 0 、 {\displaystyle \sum _{j=1}^{n}\left[Q_{j}-\left({\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial T}{\partial {\dot {q}}_{j}}}-{\frac {\partial T}{\partial q_{j}}}\right)\right]\delta q_{j}=0,} これらの仮想変位δq j は 独立でゼロではないため、係数をゼロに等しくすることができ、その結果、ラグランジュ方程式 [ 25 ] [ 26 ] または一般化された運動方程式 [ 27 ]が得られます 。 Q j = d d t ∂ T ∂ q ˙ j − ∂ T ∂ q j {\displaystyle Q_{j}={\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial T}{\partial {\dot {q}}_{j}}}-{\frac {\partial T}{\partial q_{j}}}}
これらの式は、非拘束力に対する ニュートンの法則と同等です。この式における一般化された力は、非拘束力のみから導出されます。拘束力はダランベールの原理から除外されており、求める必要はありません。一般化された力は、ダランベールの原理を満たす限り、非保存力となる可能性があります。[ 28 ]
オイラー・ラグランジュ方程式とハミルトンの原理システムが進化するにつれて、q は 構成空間 を通る経路をたどります(一部のみ表示)。システムがたどる経路 (赤) は、システムの構成 (δ q ) の小さな変化の下で定常作用 (δ S = 0) を持ちます。[ 29 ] 速度に依存する非保存力の場合、位置と速度に依存するポテンシャルエネルギー関数Vを見つけることができる かもしれません。一般化された力 Q i が ポテンシャルV から導出できる場合、[ 30 ] [ 31 ] Q j = d d t ∂ V ∂ q ˙ j − ∂ V ∂ q j 、 {\displaystyle Q_{j}={\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial V}{\partial {\dot {q}}_{j}}}-{\frac {\partial V}{\partial q_{j}}},} ラグランジュ方程式と等式を結び、ラグランジアンをL = T − V と定義すると、第 2 種のラグランジュ方程式 、すなわちオイラー・ラグランジュ 運動 方程式が得られる。∂ L ∂ q j − d d t ∂ L ∂ q ˙ j = 0. {\displaystyle {\frac {\partial L}{\partial q_{j}}}-{\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L}{\partial {\dot {q}}_{j}}}=0.}
しかし、オイラー・ラグランジュ方程式は、図に示すようにポテンシャルが見出せる場合にのみ 、非保存力を考慮することができます。非保存力の場合、常にポテンシャルが見出せるとは限らず、ラグランジュ方程式はポテンシャルを含まず、一般化された力のみを扱うため、オイラー・ラグランジュ方程式よりも一般的です。
オイラー・ラグランジュ方程式も変分法 から導かれる。ラグランジアンの 変分は δ L = ∑ j = 1 n ( ∂ L ∂ q j δ q j + ∂ L ∂ q ˙ j δ q ˙ j ) 、 δ q ˙ j ≡ δ d q j d t ≡ d ( δ q j ) d t 、 {\displaystyle \delta L=\sum _{j=1}^{n}\left({\frac {\partial L}{\partial q_{j}}}\delta q_{j}+{\frac {\partial L}{\partial {\dot {q}}_{j}}}\delta {\dot {q}}_{j}\right),\quad \delta {\dot {q}}_{j}\equiv \delta {\frac {\mathrm {d} q_{j}}{\mathrm {d} t}}\equiv {\frac {\mathrm {d} (\delta q_{j})}{\mathrm {d} t}},} これはL の全微分 に似た形式ですが、微分が仮想変位とその時間微分に置き換えられ、仮想変位の定義に従って時間増分はありません。時間に関する部分積分により、 δq j の時間微分を∂ L /∂(d q j /d t ) に変換でき、その過程で d( δq j )/d tを δq j と交換することで、独立した仮想変位をラグランジアンの微分から因数分解できます。 ∫ t 1 t 2 δ L d t = ∫ t 1 t 2 ∑ j = 1 n ( ∂ L ∂ q j δ q j + d d t ( ∂ L ∂ q ˙ j δ q j ) − d d t ∂ L ∂ q ˙ j δ q j ) d t = ∑ j = 1 n [ ∂ L ∂ q ˙ j δ q j ] t 1 t 2 + ∫ t 1 t 2 ∑ j = 1 n ( ∂ L ∂ q j − d d t ∂ L ∂ q ˙ j ) δ q j d t 。 {\displaystyle {\begin{aligned}\int _{t_{1}}^{t_{2}}\delta L\,\mathrm {d} t&=\int _{t_{1}}^{t_{2}}\sum _{j=1}^{n}\left({\frac {\partial L}{\partial q_{j}}}\delta q_{j}+{\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial L}{\partial {\dot {q}}_{j}}}\delta q_{j}\right)-{\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L}{\partial {\dot {q}}_{j}}}\delta q_{j}\right)\,\mathrm {d} t\\&=\sum _{j=1}^{n}\left[{\frac {\partial L}{\partial {\dot {q}}_{j}}}\delta q_{j}\right]_{t_{1}}^{t_{2}}+\int _{t_{1}}^{t_{2}}\sum _{j=1}^{n}\left({\frac {\partial L}{\partial q_{j}}}-{\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L}{\partial {\dot {q}}_{j}}}\right)\delta q_{j}\,\mathrm {d} t.\end{aligned}}}
さて、条件δq j ( t 1 ) = δq j ( t 2 ) = 0がすべてのj について成り立つ場合、積分されていない項はゼロになります。さらに、δL の時間積分全体がゼロである場合、δq j は独立であり、定積分がゼロになる唯一の方法は被積分関数がゼロに等しい場合であるため、 δq j の各係数もゼロでなければなりません。すると、運動方程式が得られます。これはハミルトンの原理 で要約できます。 ∫ t 1 t 2 δ L d t = 0. {\displaystyle \int _{t_{1}}^{t_{2}}\delta L\,\mathrm {d} t=0.}
ラグランジアンの時間積分は、アクション と呼ばれる別の量であり、[ 32 ]のように定義されます。 S = ∫ t 1 t 2 L d t 、 {\displaystyle S=\int _{t_{1}}^{t_{2}}L\,\mathrm {d} t,} これは関数 です。t 1 から t 2 までのすべての時間に対するラグランジュ関数を入力として受け取り、スカラー値を返します。その次元は[ 角運動量 ] 、[ エネルギー ] ・[ 時間 ] 、または[ 長さ ] ・[ 運動量 ] と同じです。この定義では、ハミルトンの原理は次のようになります。 δ S = 0. {\displaystyle \delta S=0.}
粒子が加えられた力に応じて加速すると考える代わりに、初期時刻と最終時刻に構成空間内の経路の終点が固定された状態で、粒子が静止した作用で経路を選択すると考えることもできる。ハミルトンの原理は、いくつかの作用原理 の1つである。[ 33 ]
歴史的に、粒子が力の作用を受けてたどる最短経路を見つけるという考えは、1696年にジャン・ベルヌーイ が解決した最速降下曲線問題 や、同時期にライプニッツ 、ダニエル・ベルヌーイ 、ロピタル 、そして翌年にニュートン が解決した力学の問題への変分法の最初の応用を促しました。 [ 34 ] ニュートン自身も変分法の考え方に沿って考えていましたが、発表しませんでした。[ 34 ] これらの考えは、フェルマー 、モーペルテュイ 、オイラー 、ハミルトン などの力学の変分原理につながりました。
ハミルトンの原理は、拘束方程式を座標に関する1階微分方程式の線形結合 という特定の形式にできる場合、非ホロノミック拘束 にも適用できます。結果として得られる拘束方程式は、1階微分方程式に並べ替えることができます。[ 35 ] ここではその説明は省略します。
ラグランジュ乗数と制約条件 ラグランジアンL は、 N 個の 粒子に対して、デカルト座標r k で変化させることができます。 ∫ t 1 t 2 ∑ k = 1 N ( ∂ L ∂ r k − d d t ∂ L ∂ r ˙ k ) ⋅ δ r k d t = 0. {\displaystyle \int _{t_{1}}^{t_{2}}\sum _{k=1}^{N}\left({\frac {\partial L}{\partial \mathbf {r} _{k}}}-{\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L}{\partial {\dot {\mathbf {r} }}_{k}}}\right)\cdot \delta \mathbf {r} _{k}\,\mathrm {d} t=0.}
ハミルトンの原理は、 L の 座標が独立でない場合でも依然として有効です (ここではr k ) 。ただし、拘束条件は依然としてホロノミックであると仮定されています。[ 36 ] いつものように、すべてのk に対して端点は固定されていますδ r k ( t 1 ) = δ r k ( t 2 ) = 0。δ r k は独立ではないため、δ r k の 係数 を単純にゼロに等しくすることはできません。代わりに、ラグランジュ乗数 法を使用して拘束条件を含めることができます。各拘束方程式f i ( r k , t ) = 0に i = 1, 2, ..., C のラグランジュ乗数λ i を乗じ、その結果を元のラグランジアンに加えると、新しいラグランジアンが得られます。 L ′ = L ( r 1 、 r 2 、 … 、 r ˙ 1 、 r ˙ 2 、 … 、 t ) + ∑ 私 = 1 C λ 私 ( t ) f 私 ( r k 、 t ) 。 {\displaystyle L'=L(\mathbf {r} _{1},\mathbf {r} _{2},\ldots ,{\dot {\mathbf {r} }}_{1},{\dot {\mathbf {r} }}_{2},\ldots ,t)+\sum _{i=1}^{C}\lambda _{i}(t)f_{i}(\mathbf {r} _{k},t).}
ラグランジュ乗数は時間t の任意の関数ですが、座標r k の関数ではないため、乗数は位置座標と同等の地位にあります。この新しいラグランジアンを変化させて時間に関して積分すると、次のようになります。 ∫ t 1 t 2 δ L ′ d t = ∫ t 1 t 2 ∑ k = 1 N ( ∂ L ∂ r k − d d t ∂ L ∂ r ˙ k + ∑ 私 = 1 C λ 私 ∂ f 私 ∂ r k ) ⋅ δ r k d t = 0. {\displaystyle \int _{t_{1}}^{t_{2}}\delta L'\mathrm {d} t=\int _{t_{1}}^{t_{2}}\sum _{k=1}^{N}\left({\frac {\partial L}{\partial \mathbf {r} _{k}}}-{\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L}{\partial {\dot {\mathbf {r} }}_{k}}}+\sum _{i=1}^{C}\lambda _{i}{\frac {\partial f_{i}}{\partial \mathbf {r} _{k}}}\right)\cdot \delta \mathbf {r} _{k}\,\mathrm {d} t=0.}
導入された乗数は、 r k が独立ではないにもかかわらず、δ r k の係数がゼロになるように見つけることができます。運動方程式は次のようになります。前述の解析から、この積分の解を得ることは、次の記述と同等です。 ∂ L ′ ∂ r k − d d t ∂ L ′ ∂ r ˙ k = 0 ⇒ ∂ L ∂ r k − d d t ∂ L ∂ r ˙ k + ∑ 私 = 1 C λ 私 ∂ f 私 ∂ r k = 0 、 {\displaystyle {\frac {\partial L'}{\partial \mathbf {r} _{k}}}-{\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L'}{\partial {\dot {\mathbf {r} }}_{k}}}=0\quad \Rightarrow \quad {\frac {\partial L}{\partial \mathbf {r} _{k}}}-{\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L}{\partial {\dot {\mathbf {r} }}_{k}}}+\sum _{i=1}^{C}\lambda _{i}{\frac {\partial f_{i}}{\partial \mathbf {r} _{k}}}=0,} これらは第一種ラグランジュ方程式 である。また、新しいラグランジアンに対する λ iオイラー・ラグランジュ方程式は拘束方程式を返す。 ∂ L ′ ∂ λ 私 − d d t ∂ L ′ ∂ λ ˙ 私 = 0 ⇒ f 私 ( r k 、 t ) = 0. {\displaystyle {\frac {\partial L'}{\partial \lambda _{i}}}-{\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L'}{\partial {\dot {\lambda }}_{i}}}=0\quad \Rightarrow \quad f_{i}(\mathbf {r} _{k},t)=0.}
r k 座標のみの関数であるポテンシャルエネルギーV の勾配によって与えられる保存力の場合、ラグランジアンL = T − V を代入すると次のようになる 。∂ T ∂ r k − d d t ∂ T ∂ r ˙ k ⏟ − F k + − ∂ V ∂ r k ⏟ N k + ∑ 私 = 1 C λ 私 ∂ f 私 ∂ r k = 0 、 {\displaystyle \underbrace {{\frac {\partial T}{\partial \mathbf {r} _{k}}}-{\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial T}{\partial {\dot {\mathbf {r} }}_{k}}}} _{-\mathbf {F} _{k}}+\underbrace {-{\frac {\partial V}{\partial \mathbf {r} _{k}}}} _{\mathbf {N} _{k}}+\sum _{i=1}^{C}\lambda _{i}{\frac {\partial f_{i}}{\partial \mathbf {r} _{k}}}=0,} 運動エネルギーの導関数を(負の)合力とみなし、ポテンシャルの導関数を非拘束力とみなすと、拘束力は次のようになる。 C k = ∑ 私 = 1 C λ 私 ∂ f 私 ∂ r k 、 {\displaystyle \mathbf {C} _{k}=\sum _{i=1}^{C}\lambda _{i}{\frac {\partial f_{i}}{\partial \mathbf {r} _{k}}},} これにより、拘束力は拘束方程式とラグランジュ乗数を用いて明示的に表される。
ラグランジアンの性質
非一意性 与えられたシステムのラグランジアンは一意ではありません。ラグランジアンLにゼロでない定数 a を掛け、任意の定数b だけシフトすると、新しいラグランジアンL ′ = aL + b はL と同じ運動を記述します。上記のように、与えられた時間間隔[ t st , t fin ] にわたる軌道q と固定された終点P st = q ( t st ) およびP fin = q ( t fin ) に制限すると、同じシステムを記述する 2 つのラグランジアンは、関数f ( q , t ) の「全時間微分」だけ異なる可能性があります。[ 37 ] L ′ ( q 、 q ˙ 、 t ) = L ( q 、 q ˙ 、 t ) + d f ( q 、 t ) d t 、 {\displaystyle L'(\mathbf {q} ,{\dot {\mathbf {q} }},t)=L(\mathbf {q} ,{\dot {\mathbf {q} }},t)+{\frac {\mathrm {d} f(\mathbf {q} ,t)}{\mathrm {d} t}},} どこd f ( q 、 t ) d t {\textstyle {\frac {\mathrm {d} f(\mathbf {q} ,t)}{\mathrm {d} t}}} 手段∂ f ( q 、 t ) ∂ t + ∑ 私 ∂ f ( q 、 t ) ∂ q 私 q ˙ 私 。 {\textstyle {\frac {\partial f(\mathbf {q} ,t)}{\partial t}}+\sum _{i}{\frac {\partial f(\mathbf {q} ,t)}{\partial q_{i}}}{\dot {q}}_{i}.}
ラグランジアンL とL ′ はどちらも同じ運動方程式を生成します[ 38 ] [ 39 ]。 対応するアクションS とS ′ は次のように関連付けられています。 S ′ [ q ] = ∫ t st t フィン L ′ ( q ( t ) 、 q ˙ ( t ) 、 t ) d t = ∫ t st t フィン L ( q ( t ) 、 q ˙ ( t ) 、 t ) d t + ∫ t st t フィン d f ( q ( t ) 、 t ) d t d t = S [ q ] + f ( P フィン 、 t フィン ) − f ( P st 、 t st ) 、 {\displaystyle {\begin{aligned}S'[\mathbf {q} ]&=\int _{t_{\text{st}}}^{t_{\text{fin}}}L'(\mathbf {q} (t),{\dot {\mathbf {q} }}(t),t)\,dt\\&=\int _{t_{\text{st}}}^{t_{\text{fin}}}L(\mathbf {q} (t),{\dot {\mathbf {q} }}(t),t)\,dt+\int _{t_{\text{st}}}^{t_{\text{fin}}}{\frac {\mathrm {d} f(\mathbf {q} (t),t)}{\mathrm {d} t}}\,dt\\&=S[\mathbf {q} ]+f(P_{\text{fin}},t_{\text{fin}})-f(P_{\text{st}},t_{\text{st}}),\end{aligned}}} 最後の 2 つの成分f ( P fin 、t fin ) とf ( P st 、t st )はq に依存しない。
一般化座標q の集合が与えられたとき、点変換 Q = Q ( q , t ) に従ってこれらの変数を新しい一般化座標Qの集合に変更すると、 q = q ( Q , t ) として可逆となり、新しいラグランジアンL ′ は新しい座標の関数となり、制約についても同様です。 L ′ ( Q 、 Q ˙ 、 t ) = L ( q ( Q 、 t ) 、 q ˙ ( Q 、 Q ˙ 、 t ) 、 t ) 、 ϕ j ′ ( Q 、 t ) = ϕ j ( q ( Q 、 t ) 、 t ) {\displaystyle {\begin{aligned}L'(\mathbf {Q} ,{\dot {\mathbf {Q} }},t)&=L(\mathbf {q} (\mathbf {Q} ,t),{\dot {\mathbf {q} }}(\mathbf {Q} ,{\dot {\mathbf {Q} }},t),t),\\\phi _{j}'(\mathbf {Q} ,t)&=\phi _{j}(\mathbf {q} (\mathbf {Q} ,t),t)\end{aligned}}} また、偏微分に関する連鎖律 により、ラグランジュ方程式はこの変換の下で不変である。[ 40 ]
d d t ∂ L ′ ∂ Q ˙ 私 = ∂ L ′ ∂ Q 私 + ∑ j λ j ∂ ϕ j ′ ∂ Q 私 。 {\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L'}{\partial {\dot {Q}}_{i}}}={\frac {\partial L'}{\partial Q_{i}}}+\sum _{j}\lambda _{j}{\frac {\partial \phi '_{j}}{\partial Q_{i}}}.}
証拠 座標変換の場合Q 私 = Q 私 ( q 、 t ) {\displaystyle Q_{i}=Q_{i}(\mathbf {q} ,t)} 、 我々は持っています d Q 私 = ∑ k ∂ Q 私 ∂ q k d q k + ∂ Q k ∂ t d t 、 {\displaystyle d{Q_{i}}=\sum _{k}{\frac {\partial Q_{i}}{\partial q_{k}}}d{q_{k}}+{\frac {\partial Q_{k}}{\partial t}}dt,} これはつまりQ 私 ˙ = ∑ k ∂ Q 私 ∂ q k ( q 、 t ) q ˙ k + ∂ Q k ∂ t ( q 、 t ) {\displaystyle {\dot {Q_{i}}}=\sum _{k}{\frac {\partial Q_{i}}{\partial q_{k}}}(\mathbf {q} ,t)\,{\dot {q}}_{k}+{\frac {\partial Q_{k}}{\partial t}}(\mathbf {q} ,t)} これはつまり∂ Q 私 ˙ ∂ q ˙ k = ∂ Q 私 ∂ q k {\displaystyle {\frac {\partial {\dot {Q_{i}}}}{\partial {\dot {q}}_{k}}}={\frac {\partial Q_{i}}{\partial q_{k}}}} 。
したがって、次のことも言える。 ∂ Q 私 ˙ ∂ q j = ∑ k ∂ 2 Q 私 ∂ q j ∂ q k ( q 、 t ) q ˙ k + ∂ 2 Q k ∂ q j ∂ t ( q 、 t ) {\displaystyle {\frac {\partial {\dot {Q_{i}}}}{\partial q_{j}}}=\sum _{k}{\frac {\partial ^{2}Q_{i}}{\partial q_{j}\partial q_{k}}}(\mathbf {q} ,t)\,{\dot {q}}_{k}+{\frac {\partial ^{2}Q_{k}}{\partial q_{j}\partial t}}(\mathbf {q} ,t)} 同様に: d d t ( ∂ Q 私 ∂ q j ) = ∑ k ∂ 2 Q 私 ∂ q k ∂ q j ( q 、 t ) q ˙ k + ∂ 2 Q k ∂ t ∂ q j ( q 、 t ) {\displaystyle {\frac {d}{dt}}\left({\frac {\partial {Q_{i}}}{\partial q_{j}}}\right)=\sum _{k}{\frac {\partial ^{2}Q_{i}}{\partial q_{k}\partial q_{j}}}(\mathbf {q} ,t)\,{\dot {q}}_{k}+{\frac {\partial ^{2}Q_{k}}{\partial t\partial q_{j}}}(\mathbf {q} ,t)} これは、d d t ( ∂ Q 私 ∂ q k ) = ∂ Q ˙ 私 ∂ q k {\displaystyle {\frac {d}{dt}}\left({\frac {\partial Q_{i}}{\partial q_{k}}}\right)={\frac {\partial {\dot {Q}}_{i}}{\partial q_{k}}}} 導出された2つの関係式は証明に利用できる。
一般化座標の初期セットにおけるオイラー・ラグランジュ方程式から出発して、次の式が得られます。 d d t ∂ L ∂ q ˙ 私 − ∂ L ∂ q 私 − ∑ j λ j ∂ ϕ j ∂ q 私 = 0 ∑ k ( d d t ( ∂ L ∂ Q ˙ k ∂ Q ˙ k ∂ q ˙ 私 ) − ∂ L ∂ Q k ∂ Q k ∂ q 私 − ∂ L ∂ Q ˙ k ∂ Q ˙ k ∂ q 私 − ∑ j λ j ∂ ϕ j ∂ Q k ∂ Q k ∂ q 私 ) = 0 ∑ k ( d d t ( ∂ L ∂ Q ˙ k ) ∂ Q ˙ k ∂ q ˙ 私 + ∂ L ∂ Q ˙ k d d t ( ∂ Q ˙ k ∂ q ˙ 私 ) − ∂ L ∂ Q k ∂ Q ˙ k ∂ q ˙ 私 − ∂ L ∂ Q ˙ k d d t ( ∂ Q k ∂ q 私 ) − ∑ j λ j ∂ ϕ j ∂ Q k ∂ Q k ∂ q 私 ) = 0 ∑ k ( d d t ( ∂ L ∂ Q ˙ k ) − ∂ L ∂ Q k − ∑ j λ j ∂ ϕ j ∂ Q k ) ∂ Q k ∂ q 私 = 0 {\displaystyle {\begin{aligned}{\frac {d}{dt}}{\frac {\partial L}{\partial {\dot {q}}_{i}}}-{\frac {{\partial }L}{\partial q_{i}}}-\sum _{j}\lambda _{j}{\frac {\partial \phi _{j}}{\partial q_{i}}}=0\\\sum _{k}\left({\frac {d}{dt}}\left({\frac {\partial L}{\partial {\dot {Q}}_{k}}}{\frac {\partial {\dot {Q}}_{k}}{\partial {\dot {q}}_{i}}}\right)-{\frac {{\partial }L}{\partial Q_{k}}}{\frac {\partial Q_{k}}{\partial {q}_{i}}}-{\frac {{\partial }L}{\partial {\dot {Q}}_{k}}}{\frac {\partial {\dot {Q}}_{k}}{\partial {q}_{i}}}-\sum _{j}\lambda _{j}{\frac {\partial \phi _{j}}{\partial Q_{k}}}{\frac {\partial Q_{k}}{\partial q_{i}}}\right)=0\\\sum _{k}\left({\frac {d}{dt}}\left({\frac {\partial L}{\partial {\dot {Q}}_{k}}}\right){\frac {\partial {\dot {Q}}_{k}}{\partial {\dot {q}}_{i}}}+{\frac {\partial L}{\partial {\dot {Q}}_{k}}}{\frac {d}{dt}}\left({\frac {\partial {\dot {Q}}_{k}}{\partial {\dot {q}}_{i}}}\right)-{\frac {{\partial }L}{\partial Q_{k}}}{\frac {\partial {\dot {Q}}_{k}}{\partial {\dot {q}}_{i}}}-{\frac {{\partial }L}{\partial {\dot {Q}}_{k}}}{\frac {d}{dt}}\left({\frac {\partial Q_{k}}{\partial q_{i}}}\right)-\sum _{j}\lambda _{j}{\frac {\partial \phi _{j}}{\partial Q_{k}}}{\frac {\partial Q_{k}}{\partial q_{i}}}\right)=0\\\sum _{k}\left({\frac {d}{dt}}\left({\frac {\partial L}{\partial {\dot {Q}}_{k}}}\right)-{\frac {\partial L}{\partial Q_{k}}}-\sum _{j}\lambda _{j}{\frac {\partial \phi _{j}}{\partial Q_{k}}}\right){\frac {\partial Q_{k}}{\partial q_{i}}}=0\\\end{aligned}}}
からの変化以来q → Q {\displaystyle q\rightarrow Q} が可逆であるため、オイラー・ラグランジュ方程式の形式は不変である、すなわち、 d d t ∂ L ∂ Q ˙ 私 − ∂ L ∂ Q 私 − ∑ j λ j ∂ ϕ j ∂ Q 私 = 0. {\displaystyle {\frac {d}{dt}}{\frac {\partial L}{\partial {\dot {Q}}_{i}}}-{\frac {{\partial }L}{\partial Q_{i}}}-\sum _{j}\lambda _{j}{\frac {\partial \phi _{j}}{\partial Q_{i}}}=0.}
循環座標と保存運動量 ラグランジアンの重要な特性は、保存量が そこから容易に読み取れることである。座標q i に「正準的に共役」な一般化運動量 は次のように定義される。 p 私 = ∂ L ∂ q ˙ 私 。 {\displaystyle p_{i}={\frac {\partial L}{\partial {\dot {q}}_{i}}}.}
ラグランジアンL が 何らかの座標q i に依存しない 場合、オイラー・ラグランジュ方程式から直ちに次のことが導かれる。 p ˙ 私 = d d t ∂ L ∂ q ˙ 私 = ∂ L ∂ q 私 = 0 {\displaystyle {\dot {p}}_{i}={\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L}{\partial {\dot {q}}_{i}}}={\frac {\partial L}{\partial q_{i}}}=0} 積分すると、対応する一般化運動量は定数、つまり保存量となることがわかる。これはネーターの定理 の特殊な場合である。このような座標は「巡回座標」または「無視可能な座標」と呼ばれる。
例えば、システムにはラグランジアンがあるかもしれない。 L ( r 、 θ 、 s ˙ 、 z ˙ 、 r ˙ 、 θ ˙ 、 ϕ ˙ 、 t ) 、 {\displaystyle L(r,\theta ,{\dot {s}},{\dot {z}},{\dot {r}},{\dot {\theta }},{\dot {\phi }},t),} ここで、r とz は直線に沿った長さ、s はある曲線に沿った弧長、θ とφ は角度です。z、s、φ は速度は含まれているにもかかわらず、ラグランジアンには含まれていないことに注意してください 。 すると運動 量は p z = ∂ L ∂ z ˙ 、 p s = ∂ L ∂ s ˙ 、 p ϕ = ∂ L ∂ ϕ ˙ 、 {\displaystyle p_{z}={\frac {\partial L}{\partial {\dot {z}}}},\quad p_{s}={\frac {\partial L}{\partial {\dot {s}}}},\quad p_{\phi }={\frac {\partial L}{\partial {\dot {\phi }}}},} これらはすべて保存量です。各一般化運動量の単位と性質は、対応する座標によって異なります。この場合、p z はz 方向の並進運動量、p s は曲線s に沿った並進運動量、p φ は 角度φ が測定される平面内の角運動量です。システムの運動がどれほど複雑であっても、すべての座標と速度は、これらの運動量が保存されるように変化します。
エネルギー ラグランジアンが与えられた場合L 、 {\displaystyle L,} 対応する力学系のハミルトニアンは、定義 により、 H = ( ∑ 私 = 1 n q ˙ 私 ∂ L ∂ q ˙ 私 ) − L 。 {\displaystyle H={\biggl (}\sum _{i=1}^{n}{\dot {q}}_{i}{\frac {\partial L}{\partial {\dot {q}}_{i}}}{\biggr )}-L.} 一般化座標が自然座標である場合、つまり位置ベクトルを表す際に明示的な時間依存性がない場合、この量はエネルギーと等価になります。r = r ( q 1 、 ⋯ 、 q n ) {\displaystyle \mathbf {r} =\mathbf {r} (q_{1},\cdots ,q_{n})} 。 から: T = m 2 v 2 = m 2 ∑ 私 、 j ( ∂ r → ∂ q 私 q ˙ 私 ) ⋅ ( ∂ r → ∂ q j q ˙ j ) = m 2 ∑ 私 、 j 1 私 j q ˙ 私 q ˙ j {\displaystyle T={\frac {m}{2}}v^{2}={\frac {m}{2}}\sum _{i,j}\left({\frac {\partial {\vec {r}}}{\partial q_{i}}}{\dot {q}}_{i}\right)\cdot \left({\frac {\partial {\vec {r}}}{\partial q_{j}}}{\dot {q}}_{j}\right)={\frac {m}{2}}\sum _{i,j}a_{ij}{\dot {q}}_{i}{\dot {q}}_{j}} ∑ k = 1 n q ˙ k ∂ L ∂ q ˙ k = ∑ k = 1 n q ˙ k ∂ T ∂ q ˙ k = m 2 ( 2 ∑ 私 、 j 1 私 j q ˙ 私 q ˙ j ) = 2 T {\displaystyle \sum _{k=1}^{n}{\dot {q}}_{k}{\frac {\partial L}{\partial {\dot {q}}_{k}}}=\sum _{k=1}^{n}{\dot {q}}_{k}{\frac {\partial T}{\partial {\dot {q}}_{k}}}={\frac {m}{2}}\left(2\sum _{i,j}a_{ij}{\dot {q}}_{i}{\dot {q}}_{j}\right)=2T} H = ( ∑ 私 = 1 n q ˙ 私 ∂ L ∂ q ˙ 私 ) − L = 2 T − ( T − V ) = T + V = E {\displaystyle H=\left(\sum _{i=1}^{n}{\dot {q}}_{i}{\frac {\partial L}{\partial {\dot {q}}_{i}}}\right)-L=2T-(T-V)=T+V=E} どこ1 私 j = ∂ r ∂ q 私 ⋅ ∂ r ∂ q j {\displaystyle a_{ij}={\frac {\partial \mathbf {r} }{\partial q_{i}}}\cdot {\frac {\partial \mathbf {r} }{\partial q_{j}}}} は、導出のために定義された対称行列です。
任意の時刻t において、エネルギーは構成空間 座標の変化q → Q に対して不変である。すなわち(自然座標を使用) E ( q 、 q ˙ 、 t ) = E ( Q 、 Q ˙ 、 t ) 。 {\displaystyle E(\mathbf {q} ,{\dot {\mathbf {q} }},t)=E(\mathbf {Q} ,{\dot {\mathbf {Q} }},t).} この結果に加えて、以下の証明は、このような座標変換の下で導関数が∂ L / ∂ q ˙ 私 {\displaystyle \partial L/\partial {\dot {q}}_{i}} 線形形式の係数として変化する。
証拠 座標変換Q = F ( q ) の場合、次のようになります。 d Q = F * ( q ) d q 、 {\displaystyle d\mathbf {Q} =F_{*}(\mathbf {q} )d\mathbf {q} ,} どこF * ( q ) {\displaystyle F_{*}(\mathbf {q} )} これはベクトル空間の 接写像 です{ ∑ 私 = 1 n q ˙ 私 ⋅ ( ∂ ∂ q 私 | q ) | q ˙ 私 ∈ R } {\displaystyle \left\{\sum _{i=1}^{n}{\dot {q}}_{i}\cdot \left(\left.{\frac {\partial }{\partial q_{i}}}\right|_{\mathbf {q} }\right)\ {\biggl |}\ {\dot {q}}_{i}\in \mathbb {R} \right\}} ベクトル空間へ { ∑ 私 = 1 n Q ˙ 私 ⋅ ( ∂ ∂ Q 私 | F ( q ) ) | Q ˙ 私 ∈ R } 、 {\displaystyle \left\{\sum _{i=1}^{n}{\dot {Q}}_{i}\cdot \left(\left.{\frac {\partial }{\partial Q_{i}}}\right|_{F(\mathbf {q} )}\right)\ {\biggl |}\ {\dot {Q}}_{i}\in \mathbb {R} \right\},} そして F * ( q ) = ( ∂ F 私 ∂ q j | q ) 私 、 j = 1 n {\displaystyle \textstyle F_{*}(\mathbf {q} )=\left(\left.{\frac {\partial F_{i}}{\partial q_{j}}}\right|_{\mathbf {q} }\right)_{i,j=1}^{n}} はヤコビアンです。座標ではq ˙ 私 {\displaystyle {\dot {q}}_{i}} そしてQ ˙ 私 、 {\displaystyle {\dot {Q}}_{i},} 以前の式d Q {\displaystyle d\mathbf {Q} } 形式はQ ˙ = F * ( q ) q ˙ 。 {\displaystyle {\dot {\mathbf {Q} }}=F_{*}(\mathbf {q} ){\dot {\mathbf {q} }}.} 積の法則を用いた微分後、 d Q ˙ = G ( q 、 q ˙ ) d q + F * ( q ) d q ˙ 、 {\displaystyle d{\dot {\mathbf {Q} }}=G(\mathbf {q} ,{\dot {\mathbf {q} }})d\mathbf {q} +F_{*}(\mathbf {q} )d{\dot {\mathbf {q} }},} どこ G ( q 、 q ˙ ) d q = 定義 d ( F * ( q ) ) q ˙ = ( ∑ k = 1 n ∂ 2 F 私 ∂ q j ∂ q k | q d q k ) 私 、 j = 1 n q ˙ = ( ∑ j = 1 n q ˙ j ∑ k = 1 n ∂ 2 F 私 ∂ q j ∂ q k | q d q k ) 私 = 1 、 … 、 n T = ( ∑ k = 1 n d q k ∑ j = 1 n ∂ 2 F 私 ∂ q j ∂ q k | q q ˙ j ) 私 = 1 、 … 、 n T = ( ∑ j = 1 n ∂ 2 F 私 ∂ q j ∂ q k | q q ˙ j ) 私 、 k = 1 n d q 。 {\displaystyle {\begin{aligned}G(\mathbf {q} ,{\dot {\mathbf {q} }})d\mathbf {q} &\,{\stackrel {\text{def}}{=}}\,d(F_{*}(\mathbf {q} )){\dot {\mathbf {q} }}=\left(\sum _{k=1}^{n}{\frac {\partial ^{2}F_{i}}{\partial q_{j}\partial q_{k}}}{\biggl |}_{\mathbf {q} }dq_{k}\right)_{i,j=1}^{n}{\dot {\mathbf {q} }}=\left(\sum _{j=1}^{n}{\dot {q}}_{j}\sum _{k=1}^{n}{\frac {\partial ^{2}F_{i}}{\partial q_{j}\partial q_{k}}}{\biggl |}_{\mathbf {q} }dq_{k}\right)_{i=1,\ldots ,n}^{T}\\&=\left(\sum _{k=1}^{n}dq_{k}\sum _{j=1}^{n}{\frac {\partial ^{2}F_{i}}{\partial q_{j}\partial q_{k}}}{\biggl |}_{\mathbf {q} }{\dot {q}}_{j}\right)_{i=1,\ldots ,n}^{T}=\left(\sum _{j=1}^{n}{\frac {\partial ^{2}F_{i}}{\partial q_{j}\partial q_{k}}}{\biggl |}_{\mathbf {q} }{\dot {q}}_{j}\right)_{i,k=1}^{n}d\mathbf {q} .\end{aligned}}}
ベクトル表記では、 d L ( Q 、 Q ˙ 、 t ) = ∂ L ∂ Q d Q + ∂ L ∂ Q ˙ d Q ˙ + ∂ L ∂ t d t = ( ∂ L ∂ Q F * ( q ) + ∂ L ∂ Q ˙ G ( q 、 q ˙ ) ) d q + ∂ L ∂ Q ˙ F * ( q ) d q ˙ + ∂ L ∂ t 。 {\displaystyle {\begin{aligned}dL(\mathbf {Q} ,{\dot {\mathbf {Q} }},t)&={\frac {\partial L}{\partial \mathbf {Q} }}d\mathbf {Q} +{\frac {\partial L}{\partial {\dot {\mathbf {Q} }}}}d{\dot {\mathbf {Q} }}+{\frac {\partial L}{\partial t}}dt\\&=\left({\frac {\partial L}{\partial \mathbf {Q} }}F_{*}(\mathbf {q} )+{\frac {\partial L}{\partial {\dot {\mathbf {Q} }}}}G(\mathbf {q} ,{\dot {\mathbf {q} }})\right)d\mathbf {q} +{\frac {\partial L}{\partial {\dot {\mathbf {Q} }}}}F_{*}(\mathbf {q} )d{\dot {\mathbf {q} }}+{\frac {\partial L}{\partial t}}.\end{aligned}}}
一方で、 d L ( q 、 q ˙ 、 t ) = ∂ L ∂ q d q + ∂ L ∂ q ˙ d q ˙ + ∂ L ∂ t d t 。 {\displaystyle dL(\mathbf {q} ,{\dot {\mathbf {q} }},t)={\frac {\partial L}{\partial \mathbf {q} }}d\mathbf {q} +{\frac {\partial L}{\partial {\dot {\mathbf {q} }}}}d{\dot {\mathbf {q} }}+{\frac {\partial L}{\partial t}}dt.}
先に述べたように、ラグランジアンは構成空間座標の選択に依存しない。L ( Q 、 Q ˙ 、 t ) = L ( q 、 q ˙ 、 t ) 。 {\displaystyle L(\mathbf {Q} ,{\dot {\mathbf {Q} }},t)=L(\mathbf {q} ,{\dot {\mathbf {q} }},t).} このことから、d L ( Q 、 Q ˙ 、 t ) = d L ( q 、 q ˙ 、 t ) 、 {\displaystyle dL(\mathbf {Q} ,{\dot {\mathbf {Q} }},t)=dL(\mathbf {q} ,{\dot {\mathbf {q} }},t),} そして ∂ L ∂ Q ˙ F * ( q ) = ∂ L ∂ q ˙ 。 {\displaystyle {\frac {\partial L}{\partial {\dot {\mathbf {Q} }}}}F_{*}(\mathbf {q} )={\frac {\partial L}{\partial {\dot {\mathbf {q} }}}}.} これは、各q 、 {\displaystyle \mathbf {q} ,} q ˙ 、 {\displaystyle {\dot {\mathbf {q} }},} そしてt 、 {\displaystyle t,} ∑ 私 = 1 n ∂ L ∂ q ˙ 私 d q ˙ 私 {\displaystyle \textstyle \sum \limits _{i=1}^{n}{\frac {\partial L}{\partial {\dot {q}}_{i}}}d{\dot {q}}_{i}} は、係数が明確に定義された線形形式です。∂ L ∂ q ˙ 私 {\displaystyle \textstyle {\frac {\partial L}{\partial {\dot {q}}_{i}}}} は反変1テンソルです。等式の両辺をq ˙ {\displaystyle {\dot {\mathbf {q} }}} そして上記の式を用いてQ ˙ {\displaystyle {\dot {\mathbf {Q} }}} 収量 Q ˙ ∂ L ∂ Q ˙ = q ˙ ∂ L ∂ q ˙ 。 {\displaystyle {\dot {\mathbf {Q} }}{\frac {\partial L}{\partial {\dot {\mathbf {Q} }}}}={\dot {\mathbf {q} }}{\frac {\partial L}{\partial {\dot {\mathbf {q} }}}}.} エネルギーの不変性E {\displaystyle E} 続く。
保全 ラグランジュ力学では、システムが閉じて いるのは、そのラグランジュ関数がL {\displaystyle L} エネルギーは時間に直接依存しない。エネルギー保存の法則 によれば、エネルギーはE {\displaystyle E} 閉鎖系の運動積分は、運動の積分 である。
より正確には、q = q ( t )を極値 とする。(言い換えれば、q は オイラー・ラグランジュ方程式を満たす。)この極値に沿ってL の全時間微分を取り、EL 方程式を用いると、次のようになる。 d L d t = q ˙ ∂ L ∂ q + q ¨ ∂ L ∂ q ˙ + ∂ L ∂ t − ∂ L ∂ t = d d t ( ∂ L ∂ q ˙ ) q ˙ + q ¨ ∂ L ∂ q ˙ − L ˙ − ∂ L ∂ t = d d t ( ∂ L ∂ q ˙ q ˙ − L ) = d H d t {\displaystyle {\begin{aligned}{\frac {dL}{dt}}&={\dot {\mathbf {q} }}{\frac {\partial L}{\partial \mathbf {q} }}+{\ddot {\mathbf {q} }}{\frac {\partial L}{\partial \mathbf {\dot {q}} }}+{\frac {\partial L}{\partial t}}\\-{\frac {\partial L}{\partial t}}&={\frac {d}{dt}}\left({\frac {\partial L}{\partial \mathbf {\dot {q}} }}\right){\dot {\mathbf {q} }}+{\ddot {\mathbf {q} }}{\frac {\partial L}{\partial \mathbf {\dot {q}} }}-{\dot {L}}\\-{\frac {\partial L}{\partial t}}&={\frac {d}{dt}}\left({\frac {\partial L}{\partial \mathbf {\dot {q}} }}\mathbf {\dot {q}} -L\right)={\frac {dH}{dt}}\end{aligned}}}
ラグランジアンL が明示的に時間に依存しない場合、∂ L /∂ t = 0 となり、H は 粒子の時間発展によって変化せず、実際には運動の積分となる。 H ( q ( t ) 、 q ˙ ( t ) 、 t ) = 時間の定数 。 {\displaystyle H(\mathbf {q} (t),{\dot {\mathbf {q} }}(t),t)={\text{constant of time}}.} したがって、選択された座標が自然座標であれば、エネルギーは保存される。
運動エネルギーと位置エネルギー これらの状況すべてにおいて、[ 41 ] 定数 E = T + V {\displaystyle E=T+V} は系の全エネルギーです。運動エネルギーと位置エネルギーは系が進化するにつれて変化しますが、系の運動はそれらの合計、つまり全エネルギーが一定になるようにします。これは非常に有用な簡略化です。なぜなら、エネルギーE は積分定数であり、この問題では任意の定数として扱われるため、このエネルギー関係から速度を積分して座標を求めることができる可能性があるからです。
機械的類似性 ポテンシャルエネルギーが座標の同次関数であり、時間に依存しない場合、 [ 42 ] すべての位置ベクトルが同じ非ゼロ定数α でスケーリングされ、r k ′ = α r k となるので、 V ( α r 1 、 α r 2 、 … 、 α r N ) = α N V ( r 1 、 r 2 、 … 、 r N ) {\displaystyle V(\alpha \mathbf {r} _{1},\alpha \mathbf {r} _{2},\ldots ,\alpha \mathbf {r} _{N})=\alpha ^{N}V(\mathbf {r} _{1},\mathbf {r} _{2},\ldots ,\mathbf {r} _{N})} そして、時間が係数β でスケーリングされ、t ′ = βt となると、速度v k は係数α / β でスケーリングされ、運動エネルギーT は( α / β ) 2 でスケーリングされる。ラグランジアン全体が同じ係数でスケーリングされている場合、 α 2 β 2 = α N ⇒ β = α 1 − N 2 。 {\displaystyle {\frac {\alpha ^{2}}{\beta ^{2}}}=\alpha ^{N}\quad \Rightarrow \quad \beta =\alpha ^{1-{\frac {N}{2}}}.}
長さと時間がスケーリングされているため、システム内の粒子の軌跡は、サイズが異なる幾何学的に類似した経路をたどります。元の軌跡で時間t で移動した長さl は、新しい軌跡で 時間t ′ で移動した新しい長さl ′に対応し、その比率は次のようになります。 t ′ t = ( l ′ l ) 1 − N 2 。 {\displaystyle {\frac {t'}{t}}=\left({\frac {l'}{l}}\right)^{1-{\frac {N}{2}}}.}
相互作用する粒子 与えられたシステムにおいて、2 つのサブシステムA とB が相互作用しない場合、システム全体のラグランジアンL は、サブシステムのラグランジアンL A とL B の合計になります。 [ 37 ] L = L A + L B 。 {\displaystyle L=L_{A}+L_{B}.}
相互作用する場合、これは不可能です。状況によっては、システムのラグランジアンL を、相互作用しないラグランジアンの和と、相互作用に関する情報を含む 別のラグランジアンL ABに分離できる場合があります。 L = L A + L B + L A B 。 {\displaystyle L=L_{A}+L_{B}+L_{AB}.}
これは、相互作用のないラグランジアンを運動エネルギーのみとみなし、相互作用のあるラグランジアンを系の全ポテンシャルエネルギーとみなすことで、物理的に説明できるかもしれない。また、相互作用が無視できる極限の場合、L AB は ゼロに近づき、上記の相互作用のないケースに帰着する。
相互作用しない2つ以上のサブシステムへの拡張は簡単です。全体のラグランジアンは、各サブシステムの個別のラグランジアンの合計になります。相互作用がある場合は、相互作用ラグランジアンを追加できます。
例 以下の例では、ラグランジュの第2種方程式を力学的な問題に適用する。
保守勢力 質量m の粒子は、スカラーポテンシャル の勾配 ∇から導かれる保存力 の影響下で運動する。 F = − ∇ V ( r ) 。 {\displaystyle \mathbf {F} =-{\boldsymbol {\nabla }}V(\mathbf {r} ).}
粒子の数が増えると、上記の結果に従って、全運動エネルギーはすべての粒子の運動エネルギーの合計となり、ポテンシャルはすべての座標の関数となる。
デカルト座標 粒子のラグランジアンは次のように書ける。 L ( x 、 y 、 z 、 x ˙ 、 y ˙ 、 z ˙ ) = 1 2 m ( x ˙ 2 + y ˙ 2 + z ˙ 2 ) − V ( x 、 y 、 z ) 。 {\displaystyle L(x,y,z,{\dot {x}},{\dot {y}},{\dot {z}})={\frac {1}{2}}m({\dot {x}}^{2}+{\dot {y}}^{2}+{\dot {z}}^{2})-V(x,y,z).}
粒子の運動方程式は、 x 座標 に対してオイラー・ラグランジュ方程式を適用することによって求められる。 d d t ( ∂ L ∂ x ˙ ) = ∂ L ∂ x 、 {\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial L}{\partial {\dot {x}}}}\right)={\frac {\partial L}{\partial x}},} デリバティブを含む ∂ L ∂ x = − ∂ V ∂ x 、 ∂ L ∂ x ˙ = m x ˙ 、 d d t ( ∂ L ∂ x ˙ ) = m x ¨ 、 {\displaystyle {\frac {\partial L}{\partial x}}=-{\frac {\partial V}{\partial x}},\quad {\frac {\partial L}{\partial {\dot {x}}}}=m{\dot {x}},\quad {\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial L}{\partial {\dot {x}}}}\right)=m{\ddot {x}},} したがって m x ¨ = − ∂ V ∂ x 、 {\displaystyle m{\ddot {x}}=-{\frac {\partial V}{\partial x}},} y 座標とz 座標 についても同様です。方程式をベクトル形式でまとめると、次のようになります。 m r ¨ = − ∇ V {\displaystyle m{\ddot {\mathbf {r} }}=-{\boldsymbol {\nabla }}V} これは、保存力が作用する粒子に対するニュートンの運動第二法則である。
二体中心力問題 質量m 1 とm 2 の 2 つの物体が、位置ベクトルr 1 とr 2 を持ち、引力中心ポテンシャル V によって互いの周りを公転している。位置座標をそのまま用いてラグランジアンを書き下すこともできるが、次のように 2 体問題を 1 体問題に変換する確立された手順がある。ヤコビ座標 、物体間の距離r = r 2 − r 1 、重心の 位置R = ( m 1 r 1 + m 2 r 2 )/( m 1 + m 2 ) を導入する。するとラグランジアンは次のようになる[ 44 ] [ 45 ] [ nb 4 ] L = 1 2 M R ˙ 2 ⏟ L cm + 1 2 μ r ˙ 2 − V ( | r | ) ⏟ L 関係 {\displaystyle L=\underbrace {{\frac {1}{2}}M{\dot {\mathbf {R} }}^{2}} _{L_{\text{cm}}}+\underbrace {{\frac {1}{2}}\mu {\dot {\mathbf {r} }}^{2}-V(|\mathbf {r} |)} _{L_{\text{rel}}}} ここで、M = m 1 + m 2 は全質量、μ = m 1 m 2 /( m 1 + m 2 ) は換算質量 、V は 半径方向の力のポテンシャルであり、分離の大きさ | r | = | r 2 − r 1 | のみに依存します。ラグランジアンは、重心 項L cm と相対運動 項L rel に分割されます。
R に関するオイラー・ラグランジュ方程式は単純に M R ¨ = 0 、 {\displaystyle M{\ddot {\mathbf {R} }}=0,} これは、質量中心が一定の速度で直線上を移動することを意味する。
相対運動は分離の大きさにのみ依存するため、極座標( r , θ ) を使用し、r = | r | とするのが理想的です。 L 関係 = 1 2 μ ( r ˙ 2 + r 2 θ ˙ 2 ) − V ( r ) 、 {\displaystyle L_{\text{rel}}={\frac {1}{2}}\mu \left({\dot {r}}^{2}+r^{2}{\dot {\theta }}^{2}\right)-V(r),} したがって、θ は対応する保存(角)運動量を持つ周期座標である。 p θ = ∂ L 関係 ∂ θ ˙ = μ r 2 θ ˙ = ℓ 。 {\displaystyle p_{\theta }={\frac {\partial L_{\text{rel}}}{\partial {\dot {\theta }}}}=\mu r^{2}{\dot {\theta }}=\ell .}
動径座標r と角速度d θ /d t は時間とともに変化するが、 ℓ が一定である場合に限る。r に関するラグランジュ方程式は次のようになる 。μ r θ ˙ 2 − d V d r = μ r ¨ 。 {\displaystyle \mu r{\dot {\theta }}^{2}-{\frac {dV}{dr}}=\mu {\ddot {r}}.}
この方程式は、共回転 参照系、つまり、質量が減少するにつれて回転して静止しているように見える系でニュートンの法則を使用して得られた動径方程式と同一です。この動径方程式から角速度d θ /d t を消去すると、[ 46 ] μ r ¨ = − d V d r + ℓ 2 μ r 3 。 {\displaystyle \mu {\ddot {r}}=-{\frac {\mathrm {d} V}{\mathrm {d} r}}+{\frac {\ell ^{2}}{\mu r^{3}}}.} これは、質量μ の粒子が内向きの中心力−d V /d r と、この文脈では(ラグランジュ)遠心力と呼ばれる2番目の外向きの力を受ける1次元問題の運動方程式です( 遠心力#用語のその他の使用法を 参照)。 F c f = μ r θ ˙ 2 = ℓ 2 μ r 3 。 {\displaystyle F_{\mathrm {cf} }=\mu r{\dot {\theta }}^{2}={\frac {\ell ^{2}}{\mu r^{3}}}.}
もちろん、完全に一次元的な定式化の範囲内にとどまる場合、ℓは 外部からの外向きの力の何らかの強制パラメータとしてのみ現れ、角運動量としてのその解釈は、一次元の問題の起源となったより一般的な二次元の問題に依存する。
共回転座標系でニュートン力学を用いてこの方程式を導出する場合、その解釈は、座標系自体の回転によるその座標系における遠心力として明らかです。一般化座標( r 、θ )を直接使用し、座標系について全く考えずにラグランジュ形式に従うだけでこの方程式を導出する場合、遠心力は 極座標の使用 から生じるものとして解釈されます。ヒルデブランドが言うように: [ 47 ]
「このような量は真の物理的な力ではないため、慣性力 と呼ばれることが多い。それらの存在の有無は、対象となる特定の問題ではなく、選択された座標系によって 決まる。」特に、デカルト座標系を選択した場合、遠心力は消滅し、定式化には中心力のみが含まれる。中心力は、曲線運動に対する向心力を提供する。
仮想的な力が座標の選択に由来するというこの見解は、ラグランジュ法の利用者によってしばしば表明される。この見解は、参照系が(おそらく無意識のうちに)座標の選択によって選択されるため、ラグランジュ法では自然に生じる。例えば、慣性系と非慣性系におけるラグランジアンの比較については[ 48 ]を参照。また、「全」および「更新」ラグランジュ定式化に関する議論については [ 49 ] を参照。残念ながら、この「慣性力」の使用は、ニュートン力学における慣性力の概念と矛盾する。ニュートン力学では、慣性力は座標系の選択ではなく、観測系の加速(それが慣性系ではないという事実)に由来する。物事を明確にするために、ラグランジュ慣性力を 一般化 慣性力と呼び、ニュートンベクトル慣性力と区別するのが最も安全である。つまり、ヒルデブランドが( 155ページで)「我々は常に 一般化された 力、速度、加速度、運動量を扱う。簡潔にするため、『一般化された』という形容詞は頻繁に省略する」と述べている箇所には従うべきではない。
システムのラグランジアンは一意ではないことが知られている。ラグランジアンの形式論では、ニュートンの仮想力は、仮想力が消える別のラグランジアンの存在によって特定することができ、システムの対称性を利用することによって見つかる場合もある。[ 50 ]
非保守勢力を含めるための拡張
散逸力 散逸 (すなわち非保存系)も、自由度をある方法で倍増させることで定式化された有効ラグランジアンで扱うことができる。[ 51 ] [ 52 ] [ 53 ] [ 54 ]
より一般的な定式化では、力は保存力と粘性力の両方である可能性がある。Fi から 適切 な変換 が見つかる場合、レイリーは 次の形式の散逸関数 D を使用することを提案している。 [ 55 ] D = 1 2 ∑ j = 1 m ∑ k = 1 m C j k q ˙ j q ˙ k 、 {\displaystyle D={\frac {1}{2}}\sum _{j=1}^{m}\sum _{k=1}^{m}C_{jk}{\dot {q}}_{j}{\dot {q}}_{k},} ここで、C jk は物理システムの減衰係数に関連する定数ですが、必ずしもそれらと等しいとは限りません。D がこの ように定義される場合、[ 55 ] Q j = − ∂ V ∂ q j − ∂ D ∂ q ˙ j {\displaystyle Q_{j}=-{\frac {\partial V}{\partial q_{j}}}-{\frac {\partial D}{\partial {\dot {q}}_{j}}}} そして d d t ( ∂ L ∂ q ˙ j ) − ∂ L ∂ q j + ∂ D ∂ q ˙ j = 0. {\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial L}{\partial {\dot {q}}_{j}}}\right)-{\frac {\partial L}{\partial q_{j}}}+{\frac {\partial D}{\partial {\dot {q}}_{j}}}=0.}
電磁気 テスト粒子とは、質量 と電荷が 非常に小さく、外部系への影響が無視できると仮定される粒子のことです。多くの場合、質量と電荷以外の特性を持たない、仮想的な単純化された点粒子です。電子 やアップクォーク などの実際の粒子はより複雑で、ラグランジアンにさらに項が含まれます。場は非保存ポテンシャルを形成するだけでなく、これらのポテンシャルは速度に依存することもあります。
電荷 q を持つ荷電粒子が 電磁場 と相互作用する場合のラグランジアンは、速度依存ポテンシャルの典型的な例である。電気スカラーポテンシャル ϕ = ϕ ( r , t ) と磁気ベクトルポテンシャル A = A ( r , t )は、 電場 E = E ( r , t ) と磁場 B = B ( r , t ) から次のように定義される。 E = − ∇ ϕ − ∂ A ∂ t 、 B = ∇ × A 。 {\displaystyle \mathbf {E} =-{\boldsymbol {\nabla }}\phi -{\frac {\partial \mathbf {A} }{\partial t}},\quad \mathbf {B} ={\boldsymbol {\nabla }}\times \mathbf {A} .}
電磁場中の質量のある荷電試験粒子のラグランジアン L = 1 2 m r ˙ 2 + q r ˙ ⋅ A − q ϕ 、 {\displaystyle L={\tfrac {1}{2}}m{\dot {\mathbf {r} }}^{2}+q\,{\dot {\mathbf {r} }}\cdot \mathbf {A} -q\phi ,} これは最小結合 と呼ばれます。これは、ラグランジアンが運動エネルギーから位置エネルギーを引いたものであるという一般的な経験則 が間違っている良い例です。オイラー・ラグランジュ方程式と組み合わせると、 ローレンツ力 法則 が得られます。m r ¨ = q E + q r ˙ × B {\displaystyle m{\ddot {\mathbf {r} }}=q\mathbf {E} +q{\dot {\mathbf {r} }}\times \mathbf {B} }
ゲージ変換 の下で: A → A + ∇ f 、 ϕ → ϕ − f ˙ 、 {\displaystyle \mathbf {A} \rightarrow \mathbf {A} +{\boldsymbol {\nabla }}f,\quad \phi \rightarrow \phi -{\dot {f}},} ここで、f ( r , t ) は空間と時間の任意のスカラー関数であり、前述のラグランジュ変換は次のように表される。 L → L + q ( r ˙ ⋅ ∇ + ∂ ∂ t ) f = L + q d f d t 、 {\displaystyle L\rightarrow L+q\left({\dot {\mathbf {r} }}\cdot {\boldsymbol {\nabla }}+{\frac {\partial }{\partial t}}\right)f=L+q{\frac {df}{dt}},} これは依然として同じローレンツ力法則を生み出す。
正準運動量 (位置r に共役)は、運動量に A 場からの寄与(位置運動量として知られる) を加えたものであることに注意してください。p = ∂ L ∂ r ˙ = m r ˙ + q A 。 {\displaystyle \mathbf {p} ={\frac {\partial L}{\partial {\dot {\mathbf {r} }}}}=m{\dot {\mathbf {r} }}+q\mathbf {A} .}
この関係は、量子力学 および量子場理論における 最小結合 処方でも使用されます。この式から、正準運動量 p は ゲージ不変ではないため、測定可能な物理量ではないことがわかります。ただし、rが周期的である場合(つまり、ラグランジアンが位置 r に依存しない場合)、これはϕ 場とA 場が一様である場合に発生します。この場合、ここで与えられた正準運動量p は保存運動量ですが、測定可能な物理的運動量m v は保存されません。
その他の文脈と表現 ラグランジュ力学の考え方は、物理学の他の分野にも数多く応用されており、変分法から得られる一般化された結果を取り入れることができる。
オイラー・ラグランジュ方程式は、一般化座標ではなく一般化運動量を用いて定式化することもできます。一般化座標ラグランジアンL ( q , d q /d t , t ) にルジャンドル変換を施すと、元のラグランジアンを用いた一般化運動量ラグランジアンL ′ ( p , d p /d t , t ) と、一般化運動量を用いたオイラー・ラグランジュ方程式が得られます。どちらのラグランジアンも同じ情報を含んでおり、どちらを用いても系の運動を解くことができます。実際には、一般化座標の方が一般化運動量よりも扱いやすく解釈しやすいです。
一般化座標の高階微分 一般化座標の導関数を1次のみに制限する数学的な理由は存在しません。高階導関数を含むラグランジアンに対して修正されたEL方程式を導出することは可能です(詳細はオイラー・ラグランジュ方程式 を参照)。しかし、物理的な観点からは、1次を超える時間導関数を含めることには障害があり、これはオストログラツキーによる非退化高階導関数ラグランジアンの正準形式の構築によって示唆されています(オストログラツキー不安定性を参照)。
光学 ラグランジュ力学は幾何光学 にも応用でき、媒質中の光線に変分原理を適用することで、EL方程式を解くことにより光線がたどる経路の方程式が得られる。
ラグランジュ力学は、特殊相対性理論 と一般相対性理論の 両方で定式化できます。ラグランジュ力学のいくつかの特徴は相対論的理論にも引き継がれますが、他の点ではすぐに困難が生じます。特に、EL方程式は同じ形式をとり、周期座標と保存運動量との関係は依然として適用されますが、ラグランジアンは修正する必要があり、単に粒子の運動エネルギーから位置エネルギーを引いたものではありません。また、多粒子系を明らかに共変的な 方法で扱うことは容易ではなく、特定の参照系を選択すれば可能になる場合もあります。
古典場理論 ラグランジュ力学では、一般化座標はシステムの構成を定義する離散的な変数の集合を形成します。古典場理論 では、物理システムは離散的な粒子の集合ではなく、3次元空間の領域上で定義された連続場ϕ ( r , t )です。この場にはラグランジュ密度が 関連付けられています。L ( ϕ 、 ∇ ϕ 、 ϕ ˙ 、 r 、 t ) {\displaystyle {\mathcal {L}}(\phi ,\nabla \phi ,{\dot {\phi }},\mathbf {r} ,t)} は、位置r および時刻t における場とその空間微分および時間微分によって定義されます。粒子の場合と同様に、非相対論的応用においては、ラグランジアン密度は場の運動エネルギー密度から位置エネルギー密度を差し引いたものとなります(これは一般には当てはまらず、ラグランジアン密度は「逆算」する必要があります)。ラグランジアンは、3 次元空間におけるラグランジアン密度の 体積積分となります。 L ( t ) = ∫ L d 3 r {\displaystyle L(t)=\int {\mathcal {L}}\,\mathrm {d} ^{3}\mathbf {r} } ここで、d 3 r は 3D微分 体積要素 です。ラグランジアンは時間の関数です。ラグランジアン密度は場を介して暗黙的に空間依存性を持ち、明示的な空間依存性を持つ場合もありますが、これらは積分で除去され、ラグランジアンの変数として時間だけが残ります。