wie löst man das Differentialgleichungssystem für dieses Teilchen?

Ich versuche, dieses Problem zu lösen

Ein Teilchen der Masse m bewegt sich unter der Wirkung der Schwerkraft auf der Innenfläche eines Rotationsparaboloids X 2 + j 2 = A z was reibungsfrei angenommen wird. Erhalten Sie die Bewegungsgleichungen.

Der Lagrange in Polarkoordinaten, unter der Annahme, dass die Schwerkraft in Richtung der negativen z-Achse liegt, ist

L = 1 2 M ( ρ ˙ 2 + ρ 2 φ ˙ 2 + z ˙ 2 ) M G z ; Q ¨ D Q ˙ D T D 2 Q D T 2
Einschränkungen auferlegte Bedingung impliziert
F = ρ 2 A z = 0
Lagrange-Gleichungen für dieses System sind
D D T ( L Q ˙ ich ) L Q ich = λ F Q ich
λ ist ein Lagrange-Multiplikator. Dann zu Q 1 = ρ , Q 2 = z , Q 3 = φ
M ρ ¨ M ρ φ ˙ 2 = 2 ρ λ
M z ¨ + M G = A λ
D D T ( ρ 2 φ ˙ ) = 0
ρ 2 A z = 0
ρ [ 0 , + ) , φ [ 0 , 2 π ) , z ( , + )
Ich weiß nicht, wie ich dieses Gleichungssystem lösen soll, und das Wichtigste ist, es zu bestimmen λ .

Aus dem Gleichungssystem könnte ich das ableiten

ρ 2 φ ˙ = C 0
M ρ ¨ M C 0 2 ρ 3 = 2 ρ λ
A z φ ˙ = C 0
z ˙ φ ˙ + z φ ¨ = 0


in kartesischen Koordinaten

M j ¨ = 2 j λ
M X ¨ = 2 X λ
M z ¨ + M G = A λ
X 2 + j 2 A z = 0

Die Einheiten des zweiten Terms deiner ersten eingerahmten Gleichung stimmen nicht mit dem Rest der Gleichung überein. Also würde ich die beteiligten Derivate noch einmal überprüfen.
ok korrigiert, danke...
Finden φ ˙ aus der dritten Gl und ersetzen Sie es in der 1..

Antworten (2)

Sie scheinen erfolgreich eine Reihe von Bewegungsgleichungen für das Teilchen gefunden zu haben. Ja, Sie könnten einige Variablen eliminieren, aber wenn das Problem dies nicht verlangt, besteht wirklich keine Notwendigkeit.

Wenn Sie das lösen wollen, werden die Dinge etwas komplexer. Wir wissen, dass dies ein energie- und drehimpulserhaltendes System ist, also wissen wir, dass die Bewegungsgleichungen vereinfacht werden können

ρ 2 φ ˙ = C 0 M 2 ( ρ ˙ 2 + ρ 2 φ ˙ 2 + z ˙ 2 ) + M G z = E ρ 2 A z = 0.
Um dies aus den Lagrange-Bewegungsgleichungen abzuleiten, multiplizieren Sie die ρ Gleichung durch ρ ˙ , Die z Gleichung durch z ˙ , und addiere sie zusammen. Es vereinfacht sich alles zu der oben angegebenen Energieeinsparung.

Einsetzen von für φ ˙ Und z gibt

M 2 [ ( 1 + 4 ρ 2 A 2 ) ρ ˙ 2 + C 0 2 ρ 2 ] + M G A ρ 2 = E .
Dies ist eine trennbare Gleichung erster Ordnung, also im Prinzip lösbar. Wir können viele der Konstanten loswerden, indem wir verwenden ρ ( T ) = ( A / 2 ) R ( T 2 G / A ) :
( 1 + R 2 ) R ˙ 2 + k R 2 + R 2 = ϵ ; k 8 C 0 2 A 3 G , ϵ 4 E M G A
Obwohl es merklich sauberer ist als zuvor, ist es immer noch eine Art heißes Durcheinander. So ist das Leben in Dynamik. Trennen und Verwenden von a u -Substitution gibt
1 2 R 0 2 R 2 1 + u ϵ u u 2 k D u = T T 0
Dies ist eine Lösung, wenn auch eine implizite, mit der möglicherweise schwer zu arbeiten ist.

Findest du normalerweise nicht λ , aber verwenden Sie Ihre Einschränkung, um den unphysikalischen Freiheitsgrad zu eliminieren (3 Variablen in Ihrer Parametrisierung, 2 Freiheitsgrade im System).

Genauer gesagt, differenzieren Sie die Einschränkung zweimal, stecken Sie sie in die z Bewegungsgleichung, entsprechend multiplizieren, um die gleichzusetzen ρ Gleichung. Du erhältst

M ρ ¨ M ρ ϕ ˙ 2 = 2 ρ M A ( G + 2 A ( ρ ˙ 2 + ρ ρ ¨ ) )

was dazu führt

ρ ¨ ρ ϕ ˙ 2 + 2 G A ρ + 4 ρ A 2 ( ρ ˙ 2 + ρ ρ ¨ ) = 0.

Dies ist die Bewegungsgleichung der Lagrange-Funktion (Beseitigung z von Anfang an, M = 1 )

L = 1 2 ( ρ ˙ 2 + ρ 2 ϕ ˙ 2 + 4 A 2 ( ρ ρ ˙ ) 2 ) G A ρ 2

so wie es sein sollte.

Was die Lösung der Gleichungen betrifft, weiß ich es nicht.