HAUPTTEIL: Der Lagrange
Lassen Sie die Bewegungsgleichungen und die Euler-Lagrange-Gleichungen mit null rechten Seiten ausdrücken
Q¨k+ω2kQk− 2Q˙Q∑JGkj _Q˙J−Q¨Q−Q˙2Q2∑JGkj _QJ−Q˙2Q2∑jℓ _Gjk _Gjℓ _Qℓ= 0(01a)
MQ¨+∂v( q)∂Q−1Q∑k , j( -1 _)k + jωkωJQkQJ= 0(01b)
DDT(∂L∂Q˙k) −∂L∂Qk= 0(02a)
DDT(∂L∂Q˙) −∂L∂Q= 0(02b)
WoL ( q,Q˙,Qk,Q˙k)
der Lagrange.
Wir gehen zu den folgenden Definitionen über, um die große Menge an Variablen und Indizes mittels komprimierter vereinfachter Ausdrücke zu handhaben:
Q≡def⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢Q1Q2⋮Qk⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥Q˙≡def⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢Q˙1Q˙2⋮Q˙k⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥Q¨≡def⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢Q¨1Q¨2⋮Q¨k⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥(03)
G≡def⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢0G21⋮Gk 1⋮G120⋮Gk2 _⋮G13G23⋮Gk 3⋮⋯⋯⋮⋯⋮G1 kG2 k⋮0⋮⋯⋯⋮⋯⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥= −GT(04)
Ω ( q)≡def⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢ω10⋮0⋮0ω2⋮0⋮⋯⋯⋮⋯⋮00⋮ωk⋮⋯⋯⋮⋯⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥=πQ⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢10⋮0⋮02⋮0⋮⋯⋯⋮⋯⋮00⋮k⋮⋯⋯⋮⋯⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥=ΩT( q)(05)
ϕ ( q,Q˙)≡defQ˙Q(06)
Wir definieren unten auch den reellen Skalar, so etwas wie das innere Produkt von reellen Vektoren
< Q , P >≡def∑kQkPk(07)
Unter diesen Definitionen und unter Verwendung der Gleichungen (A-01), siehe HILFSABSCHNITT, haben wir die folgenden Ausdrücke (08) anstelle der Bewegungsgleichungen (01) und (09) anstelle von (02):
Q¨+Ω2( q) Q − 2 ϕ ( q,Q˙) GQ˙−ϕ˙( q,Q˙) G Q +ϕ2( q,Q˙)G2Q = 0(08a)
MQ¨+∂v( q)∂Q−1Q<Ω2( q) Q , Q > +1Q<Ω2( q) Q , G Q > = 0(08b)
DDT(∂L∂Q˙) −∂L∂Q= 0(09a)
DDT(∂L∂Q˙) −∂L∂Q= 0(09b)
während die Lagrange-Funktion des Systems, siehe Gleichung (3.1) in der Frage, unter Verwendung der Gleichungen (A-02) ausgedrückt wird als
L ( q,Q˙, Q ,Q˙) =12<Q˙,Q˙> −12<Ω2( q) Q , Q > +12MQ˙2−V _( q) − ϕ < G Q ,Q˙> −12ϕ2<G2Q , Q >
-------------------------------------------------- ---------(10)
und in noch kompakterer Form
L ( q,Q˙, Q ,Q˙) =12(∥∥ϕ GQ − _Q˙∥∥2−∥ Ω Q ∥2) +12MQ˙2−V _(10')
Wir werden versuchen, die Lagrange-Funktion Schritt für Schritt durch ein Trial-and-Error-Verfahren zu erstellen.
Wir erwarten also, dass der 1. Term von Gleichung (08a) aus einem Lagrange-Teil stammtL1(Q˙)
so dass durch (09a)
DDT(∂L1∂Q˙) =Q¨⟹∂L1∂Q˙=Q˙(11)
Aus der Regel (A-3.d), siehe HILFSABSCHNITT,
L1
Ist
L1(Q˙) =12<Q˙,Q˙>(12)
seit
∂( <Q˙,Q˙> )∂Q˙= 2Q˙(13)
Für den 2. Term von Gleichung (08a) erwarten wir einen Lagrange-AnteilL2( q, Q )
so dass durch (09a)
−∂L2∂Q=Ω2( q) F(14)
So
L2( q, Q ) = −12<Ω2( q) Q , Q >(15)
da aus der Regel (A-3.c) und der symmetrischen (genauer: diagonalen) Matrix
Ω2=(Ω2)T
∂( <Ω2Q , Q > )∂Q= [Ω2+(Ω2)T] Q =2Ω2Q(16)
Aber als Lagrange-Teil
L2( q, Q )
ist eine Funktion von
Q
außerdem erzeugt es Elemente in den Bewegungsgleichungen, wenn es in den 2. Term von (09b) eingefügt wird:
−∂L2∂Q= +12∂( <Ω2Q , Q > )∂Q= + < Ω∂Ω∂QQ , Q >=−1Q<Ω2( q) Q , Q >(17)
das ist genau der 3. Term in Gleichung (08b).
Andererseits sind die ersten beiden Terme von (08b) die eines Teilchens, das sich in einem Potential bewegt, also stammen sie aus einem Lagrange-TeilL3( q,Q˙)
:
L3( q,Q˙) =12MQ˙2−V _( q)(18)
Dieser Teil
L3( q,Q˙)
wenn in (9a) eingesetzt, ergibt sich nichts (kein Term in Bewegungsgleichungen). Nun, in (08a) geben die Hälfte des 3. Terms und des 4. Terms an
− ϕ ( q,Q˙) GQ˙−ϕ˙( q,Q˙) G Q =DDT( − ϕ G Q )(19)
wir erwarten also einen Lagrange-Anteil
L4( q,Q˙, Q ,Q˙)
so dass durch (09a)
∂L4∂Q˙= − ϕ ( q,Q˙) GQ _(20)
das ist
L4( q,Q˙, Q ,Q˙) =−ϕ ( q,Q˙) < G Q ,Q˙>(21)
Aber wegen der Antisymmetrie von
G
, dieser Teil kann auch ausgedrückt werden als
L4( q,Q˙, Q ,Q˙) =+ϕ ( q,Q˙) < GQ˙, Q >(22)
also einfügen in den 2. Term von (09a)
−∂L4∂Q= − ϕ ( q,Q˙) GQ˙(23)
das ist die andere Hälfte des 3. Terms in (08a). Das bedeutet, dass
L4
, wenn in (09a) eingefügt, erzeugt den 3. und 4. Term von (08a)
DDT(∂L4∂Q˙) −∂L4∂Q= − 2 ϕ ( q,Q˙) GQ˙−ϕ˙( q,Q˙) GQ _(24)
Die Ausgabe der Einfügung von
L4
in (09b) würde später zusammen mit untersucht
L5
. Der 5. Term von (08a) kann aus einem Lagrange-Teil stammen
L5( q,Q˙, Q ,Q˙)
so dass durch (09a)
−∂L5∂Q=ϕ2( q,Q˙)G2Q(25)
So
L5( q,Q˙, Q ,Q˙) =−12ϕ2<G2Q , Q >(26)
da von (A-03.c) und die Symmetrie von
G2
∂( <G2Q , Q > )∂Q= (G2+(G2)T) Q =2G2Q(27)
Es kann bewiesen werden, siehe A PROOF ABSCHNITT, dass die Summe
L45=L4+L5
L45( q,Q˙, Q ,Q˙) =L4+L5= −12ϕ2<G2Q , Q >−ϕ ( q,Q˙) < G Q ,Q˙>(28)
wenn in (09b) eingefügt, entsteht der 4. Term von (08b)
DDT(∂L45∂Q˙) −∂L45∂Q= +1Q<Ω2( q) Q , G Q >(29)
In Gleichung (30) unten summieren wir die gefundenen Lagrange-Teile und die endgültige Lagrange-Funktion ist
L ( q,Q˙, Q ,Q˙) =12<Q˙,Q˙>L1−12<Ω2( q) Q , Q >L2+12MQ˙2−V _( q)L3− ϕ < G Q ,Q˙>L4−12ϕ2<G2Q , Q >L5
-------------------------------------------------- ---------(30)
identisch mit der in der Veröffentlichung angegebenen Gleichung (10).
Gleichungen (31) sind die Bewegungsgleichungen (08) mit geschweiften Klammern unter Punkten, die von den Lagrange-Termen angegeben sindLM
Diese Artikel stammen von:
Q¨L1+Ω2( q) FL2− 2 ϕ ( q,Q˙) GQ˙−ϕ˙( q,Q˙) GQ _L4+ϕ2( q,Q˙)G2QL5= 0(31a)
MQ¨+∂v( q)∂QL3−1Q<Ω2( q) Q , Q >L2+1Q<Ω2( q) Q , G Q >L4+L5= 0(31b)
Beachten Sie, dass die kanonischen Impulse
P ,p
konjugieren zu
Q ,q
jeweils sind
PP=∂L∂Q˙=Q˙−Q˙QGQ _=∂L∂Q˙= mQ˙−1Q< G Q , P >(32a)(32b)
wo für den Beweis von (32b)
P=∂L∂Q˙= mQ˙−1Q< GQ , _Q˙> −Q˙Q2<G2Q , Q >= mQ˙−1Q< GQ , _Q˙> +Q˙Q2< G Q , G Q >= mQ˙−1Q< GQ , _Q˙−Q˙QGQ _P> = mQ˙−1Q< G Q , P >(32b')
Die Gleichungen (32a) und (32b) sind identisch mit (3.3) bzw. (3.4) der unten angegebenen Veröffentlichung
PkP=Q˙k−Q˙Q∑JGkj _QJ= mQ˙−1Q∑jk _Gkj _PkQJ(3.3)(3.4)
HILFSTEIL : Komprimierte vereinfachte Ausdrücke und partielle Differenzierungsregeln
Gleichungen (A-01) sind nützlich für die Umwandlung der Bewegungsgleichungen von der Form (01) in die Form (08) :
ω2kQk=[Ω2( q) Q ]k∑JGkj _QJ=[ GQ ] _k∑JGkj _Q˙J=[ GQ˙]k∑jℓ _Gjk _Gjℓ _Qℓ= −∑ℓ(∑JGkj _Gjℓ _)Qℓ= −∑ℓ(G2)kℓ _Qℓ= −(G2F )kQ¨Q−Q˙2Q2=DDT(Q˙Q) =Dϕ ( q,Q˙)DT=ϕ˙( q,Q˙)∑k , j( -1 _)k + jωkωJQkQJ= <Ω2( q) Q , Q > − <Ω2( q) Q , G Q >(A-01.a)(A-01.b)(A-01.c)(A-01.d)(A-01.e)(A-01.f)
Der Beweis von (A-01.f) läuft wie folgt ab
∑k , j( -1 _)k + jωkωJQkQJ=∑kω2kQ2k<Ω2( q) Q , Q >+∑k , j ≠ k( -1 _)k + jωkωJQkQJ− <Ω2( q) Q , G Q >(A-01.f')
seit
∑k , j ≠ k( -1 _)k + jωkωJQkQJ=(πQ)2∑k , j ≠ k( -1 _)k + jkj _QkQJ=(πQ)2∑k , j ≠ k( -1 _)k + j2 k jJ2−k2Gkj _J2−k22QkQJ=12∑k , jGkj _(ω2J−ω2k)QkQJ=−12∑J(ω2JQJ)[Ω2( q) Q ]J∑kGjk _Qk[ GQ ] _J−12∑k(ω2kQk)[Ω2( q) Q ]k∑JGkj _QJ[ GQ ] _k= − <Ω2( q) Q , G Q >
-------------------------------------------------- ---------(A-01.f„ “)
Die Gleichungen (A-02) und (A-03) sind nützlich für die Umwandlung der Lagrange-Funktion von der Form (3.1), siehe fragliche Gleichung, in die Form (10) und für die schrittweise Konstruktion dieser Lagrange-Funktion aus der Bewegungsgleichungen (08) :
∑kQ˙2k= <Q˙,Q˙> =∥∥Q˙∥∥2∑kω2k( q)Q2k= <Ω2Q , Q >=<Ω Q ,ΩTQ >=<Ω Q ,Ω Q >=∥ Ω Q ∥2∑j , kGkj _Q˙kQJ= < GQ , _Q˙> = − < GQ˙, Q >∑j , k , lGkj _Gkℓ _QℓQJ= − <G2Q , Q >=< G Q , G Q >=∥ GQ ∥ _2(A-02.a)(A-02.b)(A-02.c)(A-02.d)
Die Gleichungen (A-02.c) und (A-02.d) werden jeweils wie folgt bewiesen
∑j , kGkj _Q˙kQJ=∑k(∑JGkj _QJ)Q˙k=∑k[ GQ ] _k[Q˙]k=< GQ , _Q˙> = < Q ,GTQ˙> = < Q , − GQ˙> = − < GQ˙, Q >(A-02.c')
∑j , k , lGkj _Gkℓ _QℓQJ=∑k(∑JGkj _QJ) (∑ℓGkℓ _Qℓ) =∑k[ GQ ] _k[ GQ ] _k= < G Q , G Q > = <GTG Q , Q >=−<G2Q , Q >(A-02.d')
Die nachstehenden Gleichungen (A-03) sind gewissermaßen partielle Differentiationsregeln einer Skalarfunktion einer VektorvariablenS
in Bezug auf diese Variable. Die Skalarfunktionen sind normalerweise Skalarprodukte und der variable Vektor ist esS = QoderQ˙
. Im FolgendenA , R
sind Vektoren undF
lineare Transformation alle unabhängig vom variablen VektorS
. NormalerweiseF =Ω,Ω2, G ,G2
:
∂( < A , S > )∂S=∂( < S , A > )∂S= A∂( < R , F S > )∂S=∂( <FTR , S > )∂S=FTR∂( < F S , S > )∂S= ( F +FT) S∂( < S , S > )∂S= 2 S(A-03.a)(A-03.b)(A-03.c)(A-03.d)
(A-03.b) ist ein Sonderfall von (A-03.a) mit
A =FTR
und (A-03.d) ist ein Spezialfall von (A-03.c) mit
F = Ich
.
Eine im folgenden Abschnitt nützliche Identität ist
GT= − G⟹< G S , S > = 0 ,für jeden reellen Vektor S(A-04)
seit
< G S , S > = < S ,GTS >=< S , ( - G ) S >=-< G S , S >(A-04')
A BEWEISABSCHNITT: Beweis von Gleichung (29) bei gegebener Gleichung (28).
Wir werden Gleichung (29) aus (28) beweisen, wobei die beiden Gleichungen hier der Einfachheit halber wiederholt werden
DDT(∂L45∂Q˙) −∂L45∂Q= +1Q<Ω2( q) Q , G Q >(29)
Wo
L45( q,Q˙, Q ,Q˙)≡def−12ϕ2( q,Q˙) <G2Q , Q >−ϕ ( q,Q˙) < G Q ,Q˙>(28)
−∂L45∂Q= ϕ∂ϕ∂Q<G2Q , Q >+∂ϕ∂Q< GQ , _Q˙>= (Q˙Q)∂(Q˙Q)∂Q<G2Q , Q >+∂(Q˙Q)∂Q< GQ , _Q˙>
So
−∂L45∂Q= ( -Q˙2Q3) <G2Q , Q >+ ( −Q˙Q2) < G Q ,Q˙>(B-01)
Jetzt
∂L45∂Q˙∂L45∂Q˙= − ϕ∂ϕ∂Q˙<G2Q , Q >−∂ϕ∂Q˙< GQ , _Q˙>⟹= ( -Q˙Q2) <G2Q , Q >+ ( −1Q) < G Q ,Q˙>(B-02)
Differenzieren (B-02) bzgl
T
DDT(∂L45∂Q˙) = ( -Q¨Q− 2Q˙2Q3) <G2Q , Q >+ ( −Q˙Q2) <G2Q˙, Q >+ ( -Q˙Q2) <G2F ,Q˙> + (Q˙Q2) < G Q ,Q˙> + ( −1Q)< GQ˙,Q˙>= 0 , siehe (A-04)+ ( -1Q) < G Q ,Q¨>(B-03)
Hinzufügen von (B-01) und (B-03)
DDT(∂L45∂Q˙) −∂L45∂Q=( -Q¨Q−Q˙2Q3) <G2Q , Q >+ ( −2Q˙Q2) <G2F ,Q˙> + ( −1Q) < G Q ,Q¨> =+1Q< (Q¨−Q˙2Q2) G Q +2Q˙QGQ˙−Q¨, G Q > = +1Q<ϕ˙G Q +2ϕ GQ˙−Q¨=Ω2( q) F +ϕ2G2Q , siehe (08a), G Q >+1Q<Ω2( q) F +ϕ2G2Q , G Q >=+1Q<Ω2( q) Q , G Q > +ϕ2Q<G2Q , GQ > _= 0 , siehe (A-04)
So
DDT(∂L45∂Q˙) −∂L45∂Q= +1Q<Ω2( q) Q , G Q >(B-04)
QED.
udrv