Lagrange Formalismus - Berlin

Inhaltsverzeichnis
Lagrange Formalismus
Frank Essenberger
FU Berlin
10.Oktober 2006
1 Oszillatoren 1
1.1 Fadenpendel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Stabpendel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3 U-Rohr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2 Perlen 5
2.1 Perle auf rotierenden Stab . . . . . . . . . . . . . . . . . . . . . . 5
2.2 Perle auf Schraubenlinie . . . . . . . . . . . . . . . . . . . . . . . 7
2.3 Perle fällt durch Kugel . . . . . . . . . . . . . . . . . . . . . . . . 8
1 Oszillatoren
1.1 Fadenpendel
Wir nehmen an eine Punktmasse der Größe m hänge an einem masselosen Faden
der Länge S und zwar im Kraftfeld der Erde. Wir haben eine Zwangsbedingung,
nämlich dass das Seil gespannt seinen soll und sich nicht dehnen kann:
−→ r (t) · ←− r (t) = S 2 .
Wenn man den Ursprung des Koordinantensystems in den Punkt legt wo
das Pendel befestigt ist ist ergibt sich:
x(t) = S · sin(ϕ(t)) ⇒ ˙x(t) = S · cos(ϕ(t)) ˙ϕ(t)
y(t) = S · cos(ϕ(t)) ⇒ ˙y(t) = −S · sin(ϕ(t)) ˙ϕ(t).
1

Abbildung 1: Fadenpendel
Durch die eine generalisierte Kordinate reduziert sich die Anzahl der Variabeln
von zwei ( −→ ( )
x(t)
r (t) = ) auf eine ϕ(t). Nun muss die Lagrange Gleichung
y(t)
bestimmt werden. Dies ist hier ganz einfach:
T = 1
2 mv2
V = −mgh
L = T − V = 1
2 m ˙−→ r (t) · ˙←− r (t) − mgy(t)
= 1
2 m[ ˙x(t) 2 + ˙y(t) 2] −mgy(t)
= 1
2 mS2 ˙ϕ(t) 2[ sin(ϕ(t)) 2 + cos(ϕ(t)) 2] +mgS cos(ϕ(t))
L = 1
2 mS2 ˙ϕ(t) 2 − mgS cos(ϕ(t)).
Nun bilden wir die Euler Lagrange Gleichung für die generalisierte Kordinate
ϕ(t):
0 = d [ ∂
dt ∂ ˙ϕ L] − ∂
∂ϕ L
0 = d
[
∂ [1
dt ∂ ˙ϕ 2 mS2 ˙ϕ(t) 2 − mgS cos(ϕ(t)) ]]
− ∂ [1
∂ϕ 2 mS2 ˙ϕ(t) 2 + mgS cos(ϕ(t)) ]
0 = d [ ] 2
mS ˙ϕ(t) +mgS sin(ϕ(t))
dt
0 = S ¨ϕ(t) + g sin(ϕ(t)).
Für kleine Auslenkungen ist sin(ϕ(t)) ≈ ϕ(t) und so ergibt sich:
0 = ¨ϕ(t) + g
S ϕ(t)).
Die Lösung sieht man sofort:
√ √
g
g
ϕ(t) = A cos( t) + B sin(
S S t).
Die Anpassung an Randbedingungen, die die Konstanten festlegen bleibt
dem Leser überlassen.
2

1.2 Stabpendel
Dieses Pendel besteht aus einem steifen Stab der Länge S und der Masse m
(homogen verteilt). Auch er ist fest befestigt und in homogenen Schwerefeld
der Erde oder einem anderen homogenen Feld. Unsere Zwangsbedingung ist
diesmal, dass der Stab steif ist. So lässt sich −→ r (t) wieder nur mit einer Kordinate
ausdrücken:
−→ r (t) =
( S ′ sin(ϕ(t))
S ′ cos(ϕ(t))
)
.
Abbildung 2: Stabpendel
Nun muss wieder T und V bestimmt werden. Dazu erstmal ˙r ′ (t) 2 bestimmen.
r wurde deshalb gestrichen, da es vom Abstand S ′ zum Ursprung abhängt:
r ˙−→ ′
)
(t) =
.
( S ′ cos(ϕ(t)) ˙ϕ(t)
−S ′ sin(ϕ(t)) ˙ϕ(t)
˙r ′ (t) 2 = S ′2 [cos(ϕ(t)) 2 ˙ϕ(t) 2 + sin(ϕ(t)) 2 ˙ϕ(t) 2 ] = S ′2 ˙ϕ(t) 2
Damit wird T zu:
Und V ist:
T =
=
∫ S
dS
0
∫ S
dS
0
′ 1 m
2 S ˙r′ (t) 2
= 1
6 mS2 ˙ϕ(t) 2 .
′ 1 m
2 S S′2 ˙ϕ(t) 2
∫ S
′ m
V = − g dS
0 S S′ cos(ϕ(t))
= − 1
Sgm cos(ϕ(t)).
2
Wir bilden wieder L = T − V und bilden so wie im ersten Beispiel die
3

Eulerlagrange Gleichung:
0 = d [ ∂
dt ∂ ˙ϕ L] − ∂
∂ϕ L
0 = d [ ∂ 1
dt ∂ ˙ϕ 6 mS2 ˙ϕ(t) 2 + 1
2 Sgm cos(ϕ(t))] − ∂ [1
∂ϕ 6 mS2 ˙ϕ(t) 2 + 1
Sgm cos(ϕ(t))]
2
0 = d [ ∂ 1
dt ∂ ˙ϕ 6 mS2 ˙ϕ(t) 2] − ∂ [1
Sgm cos(ϕ(t))]
∂ϕ 2
0 = 1
3 mS2 ¨ϕ(t) + 1
Sgm sin(ϕ(t))
2
0 = ¨ϕ(t) + 3 g
2 S sin(ϕ(t)).
Für kleine Auslenkungen ist sin(ϕ(t)) ≈ ϕ(t) und so ergibt sich:
0 = ¨ϕ(t) + 3 g
2 S ϕ(t).
Die Lösung sieht man sofort:
√ √
3g
3g
ϕ(t) = A cos( t) + B sin(
2S 2S t).
1.3 U-Rohr
Hier befindet sich in einem U-Rohr der Dicke A eine Flüssigkeit der Dichte ρ.
Diese soll sich reibungsfrei und inkompressibel sich im Rohr bewegen können.
Das heißt die Geschwindigkeit der Oberfläche ist gleich der Geschwindigkeit
überall im Rohr.Unsere generalisiertekordinate ist die Höhe h(t). Wir müssen
nun wieder V und T bestimmen.
Abbildung 3: U-Rohr
Um ein Volumenelement um ∆y anzuheben muss ∆V an sich anheben und
auch noch auf der anderen Seite ∆V um ∆y aus der Ruhelage heben. Damit
folgt:
∫ y
V = − dh2 · hρAg
o
= −h 2 ρAg
4

Und die kinetische Energie ist:
T = 1
2 m˙ h 2 .
Nun wieder Die Euler Lagrangegleichung bestimmenL = T − V :
0 = d [ ∂
dt ∂ ˙ h L] − ∂
∂h L
0 = d [ ∂
dt ∂ ˙ 1
h 2 m˙ h 2 − 2 · mgh ] − ∂ [1
∂h 2 m˙ h 2 − gh 2 ρA ]
0 = d [ ]
mh ˙ +2 · ghρA
dt
0 = ¨ h + 2 · ghρA
0 = ¨ h + 2 · ghρA
m
0 = ¨ h + 2 · ghρA
ρV
0 = ¨ h + 2g
· h
l
Die Lösung sieht man wieder leicht:
√ √
2g
2g
h(t) = A cos( t) + B sin(
l l t).
2 Perlen
Ganz häufige Aufgaben sind irgendwelche gleitenden Perlen auf Stäben, Ringen
oder ähnlichen. Dabei wird die Reibung meistens nicht berücksichtigt.
2.1 Perle auf rotierenden Stab
Angenommen eine Perle sei auf einen Stab gesteckt, der sich mit konstanter
Winkelgeschwindigkeit in der X-Y Ebene bewegt.
Abbildung 4: Gleitende Perle auf rotirendem Stab
5

Man nimmt dabei an, dass die X-Y Ebene eine Äquipotentialfläche ist. Wir
führen also wieder passende generalisierte kordinaten ein:
x(t) = l(t) cos(ϕ(t))
y(t) = l(t) sin(ϕ(t)).
Hier ist alles nun sehr leicht. V ist einfach:
Und die kinetische Energie ist:
V = const.
T = 1
2 m[ ˙x(t)2 + ˙y(t) 2 ]
= 1
2 m[l(t)2 sin(ϕ(t)) 2 ˙ϕ(t) 2 + l(t) 2 cos(ϕ(t)) 2 ˙ϕ(t) 2 + ˙ l(t) 2 sin(ϕ(t)) 2 + ˙ l(t) 2 cos(ϕ(t)) 2 ]
= 1
2 [ml(t)2 ˙ϕ(t) 2 + ˙ l(t) 2 ]
Damit ergibt sich als Lagrangegleichung für l(t):
0 = d [ ∂
dt ∂ ˙ h L] − ∂
∂h L
0 = d [ ∂
dt ∂ ˙ 1
l 2 [ml(t)2 ˙ϕ(t) 2 + ˙ l(t) 2 − const] ] + ∂ [1
∂l 2 [ml(t)2 ˙ϕ(t) 2 + ˙ l(t) 2 − const ] ]
0 = d
l(t) + ml(t) ˙ϕ(t)2
dt
0 = ¨l(t) + ml(t) ˙ϕ(t) 2 .
Da der Stab sich mit konstanter Winkelgeschwindigkeit dreht ist ϕ(t) = ωt
also wird die Gleichung leicht:
0 = ¨ l(t) + mω 2 l(t).
Die Lösund sieht man wieder sofort:
l(t) = Ae −√mω2 √
t mω2 + Be
t.
Für eine ruhende Perle im Abstand l0 zum Zeitpunkt t=0 ergibt sich (l(0) =
l0 ˙ l(0) = 0):
l0 = A + B
∧ 0 = √ mω 2 B − √ mω 2 A.
Aus dem LGS sieht man sofort, dass A = B und A = l0
2 ist. Damit
vereinfacht sich diese Lösung der DGL zu:
l(t) = l0 cosh( √ mω 2 t).
6

Die Bahnkurve sieht dann so aus:
Abbildung 5: Bahnkurve der Perle
2.2 Perle auf Schraubenlinie
Diese Perle der Masse m soll nun im Schwerefeld der Erde reibungsrei entlang
einer Schraubenlinie rutschen.
Abbildung 6: Perle auf Schraubenlinie
Bei einer Windung soll die Perle dabei die Strecke S zurücklegen. Die Wahl
der Kordinaten liegt wieder nahe. In den Zylinderkordinaten lassen sich die
Zwangsbedingungen:
x(t) 2 + y(t) 2 = R
z(t) = ϕ(t)
leicht implementieren. So wird aus −→ r (t) = ( x(t) y(t) z(t) ) ein Vektor der
nurnoch von einer variablen abhängt: −→ r (t) = ( R cos(ϕ(t)) R sin(ϕ(t)) Sϕ(t) ) .
7

Die kinetische und Potentielle Energie zu bestimmen ist nun Routine:
V = − mgz
= − mgSϕ(t).
Die kinetische energie ist nur hinzuschreiben:
T = 1
= 1
2 m
[
2 m[ ˙x(t)2 + ˙y(t) 2 + ˙z(t) 2 ]
[R cos(ϕ(t)) ˙ϕ(t)] 2 + [R sin(ϕ(t)) ˙ϕ(t)] 2 + [S ˙ϕ(t)] 2
= 1
2 m[R2 ˙ϕ(t) 2 + S 2 ˙ϕ(t) 2 ]
= 1
2 m ˙ϕ(t)2 [R 2 + S 2 ].
Nun wieder Die Euler Lagrangegleichung bestimmenL = T − V :
0 =
0 =
0 =
d [ ∂
dt ∂ ˙ h L] − ∂
∂h L
[
d ∂ 1
[
dt ∂ ˙ϕ 2 m ˙ϕ(t)2 [R 2 + S 2 ]
] − mgSϕ(t) ] − ∂
[
1
∂ϕ 2 m ˙ϕ(t)2 [R 2 + S 2 ]
] − mgSϕ(t)
d
dt [m ˙ϕ(t)2 [R 2 + S 2 ]] + ∂
∂ϕ mgSϕ(t)
0 = m ¨ϕ(t) 2 [R 2 + S 2 ]] + mgS
gS
[R2 + S2 ] .
⇒ ¨ϕ(t) = −
Nun muss man nurnoch zweimal überintegrieren um ϕ(t) zu erhalten:
ϕ(t) =
−gS
2[R 2 + S 2 ] · t2 − ˙ϕo · t − ϕ0.
Für eine perle die einfach nur fällt wäre R = 0 und S = 1 so ergäbe sich:
so wie erwartet.
2.3 Perle fällt durch Kugel
ϕ(t) = z(t) = − 1
2 g · t2 − ˙zo · t − z0.
Dies ist ein Beispiel in dem Nur das Potential zu finden ein Problem darstellt.
Wenn man es hat kann man sich ganz einfach nach Schema f zum Ergebniss
durchhangeln. Angenommen man hatt durch eine Kugel einen schacht gebohrt,
durch den man eine Kugel fallen lässt. Wichtig für die Aufgabe ist, dass keine
äußeren Kräfte wirken, und nur die Kugel das Potential erzeugt. Gesucht ist
nun also das Potential im inneren einer Kugel der Masse M.
8
]

Abbildung 7: Perle fliegt durch Kugel
Wenn man annimmt, das die Dichtefunktion ρ( −→ x ) aus konzentrischen Kugelschalen
gleicher Dichte aufgebaut ist, dann muss das Feld im inneren der
Kugel symetrisch sein. Wir wollen uns auf den einfacheren Fall ρ( −→ x ) = ρ0 be-
schränken. Im weiteren ist | −→ x | = r und | −→ x ′ | = r ′ . Dann ist das Potential
definiert als:
φ( −→ ∫
x ) = γ dx
V
′3 ρ(−→ x ′ )
| −→ x − −→ x ′ ∫
= γ
| V
dx ′3 ρ0Θ[RKugel − r ′ ]
√∑
(xi − x ′ .
i )2
Wenn man nun den Schacht entlang der z-Achse gebohrt hat, ist nur −→ ( ) x =
0 0 x3 interessant, damit vereinfacht sich die Gleichung zu:
φ( ( 0 0
)
x3 ) = γρ0
∫
= γρ0
= γρ0
∫
V
V
∫ ∞
0
dx ′3 Θ[RKugel − r ′ ]
√∑
(xiδ3i − x ′ ∫
= γρ0
i )2 V
dx ′3
Θ[RKugel − r ′ ]
dx ′3
√ r 2 + [r ′ cos(θ)] 2 − 2r ′ cos(θ)x3
dr ′ r ′2
∫ 2π
0
dφ ′
∫ π
0
Θ[RKugel − r ′ ]
√ r 2 + x ′2
3 − 2x ′ 3 x3
Θ[RKugel − r ′ ]
dθ ′ sin(θ ′ ) √
r2 + [r ′ cos(θ ′ )] 2 − 2r ′ cos(θ ′ )x3
Nun nuzt man die Θ[RKugel − r ′ ] aus und Substituiert noch cos(θ ′ ) mit y’:
φ( −→ x ) = γρ0
∫ R Kugel
0
dr ′ r ′2
∫ 2π
0
dφ ′
∫ 1
−1
dy ′
1
√ r 2 + [r ′ y ′ ] 2 − 2r ′ y ′ x3
Das Integral ist sehr schwer zu lösen, deshalb bedienen wir uns der Piosson-
9
.
.

Gleichung und nutzen die Symetrie des Problems aus:
φ( −→ x ) =
∫
γ dx
V
′3 ρ(−→ x ′ )
∇ 2 φ( −→ x ) = ∇ 2 γ
∫
V
| −→ x − −→ x ′ |
dx ′3 ρ(−→ x ′ )
| −→ x − −→ x ′ = γ
|
∫
dx
V
′3 ρ( −→ x ′ )∇ 2
1
| −→ x − −→ x ′ .
|
Nun die δ −F unktion einsetzen und das Konservative Gradientenfeld ←− ∇φ =
− ←− E ( −→ x ) ausnutzen:
∇ 2 φ( −→ x ) =
∫
γ dx
V
′3 ρ( −→ x ′ )∇ 2 1
| −→ x − −→ x ′ −→ ←−
∇ · ∇φ =
|
∫
γ dx ′3 ρ( −→ x ′ )δ( −→ x − −→ x ′ )
V
− −→ ∇ ←− E ( −→ x ) = γρ( −→ x ). (1)
Für unser Problem ergibt sich also:
− −→ ∇ ←− E ( −→ x ) = γρ( −→ x ).
∫
dV ′ − −→ ∇ ←− E ( −→ ∫
x ) = dV ′ γρ0
v
v
∫
− d
S
−→ S ′ · ←− E ( −→ x )| |
−→
x |=r0 = 4
3 πr3 0γρ0.
Da, dass Problem Kugelsymetrisch ist, muss das Feld auf einer Kugeloberfläche
mit Radius r0 = const. sein und in richtung der Flächennormalen −→ n zeigen.
So ergibt sich:
∫
−E(r0) dS
S
−→ n ←− n = 4
3 πr3 0γρ0
−E(r0)4πr 2 0 = 4
3 πr3 0γρ0
Damit ergibt sich für das Potential:
⇒ E(r) = − γρ0
3 r.
φ( −→ x ) = γρ0
6 r2 = Γr 2 .
Wir gehen nun davon aus, dass der Schacht für die Perle ganz dünn ist und
das Potential unverändert bleibt. Dann ergibt sich als L = V − T :
L = 1
2 m ˙r(t)2 − Γr 2 .
10

Wenn nun der Schacht auf der z-Achse liegt ist die Bewegung der Perle auf
−→ x = ( 0 0 x3(t) ) festgelegt und es ergibt sich:
L = 1
2 m ˙z(t)2 − Γz 2 .
Nun wieder Die Euler Lagrangegleichung bestimmen:
0 = d [ ∂
dt ∂ ˙z L] − ∂
∂z L
0 = d [ ∂ 1
dt ∂ ˙z 2 m ˙z(t)2 − Γz 2] − ∂ [1
∂z 2 m ˙z(t)2 − Γz 2]
0 = d [ ]
m ˙z(t) +2Γz
dt
0 = m¨z(t) + 2Γz.
Die Lösung sieht man wieder schnell:
Mit der Kreisfrequenz ω =
z(t) = A cos(ωt) + B sin(ωt).
√
2Γ
m = √ γρ0
3m
11