This is an old revision of the document!
−Table of Contents
N-Körper-Problem
Quellen
Open Hard- & Software [ DokuWiki WebSites MediaWiki NextCloud ]
Prof.Dr. Edmund Weitz, HAW Hamburg: Kepler, Newton und das Dreikörperproblem
Vorgaben
Newtonsche Axiome
- Kräftefreie Masse bewegt sich mit konstanter Geschwindigkeit.</li>
- Kraft gleich Masse mal Beschleunigung.</li>
- Kraft erzeugt Gegenkraft - Kraft auf Masse A erzeugt Gegenkraft mit gleichem Betrag und umgekehrter Richtung.
(1.) →F=→0⇒→a=→0⇒d→vdt=→0⇒→v=const
(2.) →F=m⋅→a
(3.) →Fab=−→Fba
Gravitationsgesetz
→Fm1m2=Gm1m2r3m1m2⋅→rm1m2
→Fm1m2=Gm1m2r2m1m2⋅→erm1m2
Annahmen
- - Zwei Massen bewegen sich immer in einer Ebene
- - Es gibt nur Kräfte in Richtung der Verbindungsvektoren der Massen
- - Drehimpuls bleibt erhalten und damit konstant:
→L=C(→r×→v) , d→Ldt=→0
Definitionen
Definition - Punkt in Polarkoordinaten: PPCS=P(φ(t),r(t))
Definition - Punkt in Kartesischen Koordinaten: PCCS=P(x(t),y(t)).
Definition - Umrechnung von PCS nach CCS: PCCS=(xy)=(rcos(φ)rsin(φ))\larrPPCS
PPCS=(φr)=(arctan(yx)√x2+y2)\larrPCCS
Definition - Drehung des CCS um Ursprung mit Winkel: Drehung um α von CCS\rarrCCS′ : (x′y′)=(+cos(α)−sin(α)+sin(α)+cos(α))(xy)
Rücktransformation um α von CCS′⇒CCS : (xy)=(+cos(α)+sin(α)−sin(α)+cos(α))(x′y′)
Definition - Lokales Planeten-(Hilfs-)Koordinatensystem
Definition zweier senkrechter Einheitsvektoren →u(t) und →w(t) :
Bild!
Die Einheitsvektoren →u(t) und →w(t) unterteilen: - den Ortsvektor →r(t) in Betrag ‖→r(t)‖=r(t) und Richtung →u(t) - den Geschwindigkeitsvektor →v(t) in Betrag ‖→v(t)‖=v(t) und Richtung →w(t)
<h5>Einheitsvektor Ort</h5>
Einheitsvektor Ort: →u(t)=(uxuy) →u(t)=→r(t)r(t)=(cos(φ(t))sin(φ(t)))
Ableitungen des Orts-Einheitsvektors →u(t) nach der Zeit:
d→u(t)dt=ddt(cos(φ(t))sin(φ(t)))=(−dφ(t)dtsin(φ(t))+dφ(t)dtcos(φ(t)))=dφ(t)dt(−sin(φ(t))+cos(φ(t)))=dφ(t)dt→w(t)
Kurzschreibweise: ˙→u(t)=˙φ(t)→w(t)
<h5>Einheitsvektor Geschwindigkeit</h5>
Einheitsvektor Geschwindigkeit: →w(t)=(wxwy)
Drehung →u(t) um +π2:→u(t)\rarr→w(t)
(wxwy)=(+cos(π2)−sin(π2)+sin(π2)+cos(π2))(uxuy)
(wxwy)=(0−1+10)(uxuy)
(wxwy)=(0−1+10)(cos(φ(t))sin(φ(t)))
→w(t)=(−sin(φ(t))cos(φ(t)))
Ableitungen des Geschwindigkeits-Einheitsvektors →w(t) nach der Zeit:
d→w(t)dt=ddt(−sin(φ(t))+cos(φ(t)))=(−dφ(t)dtcos(φ(t))−dφ(t)dtsin(φ(t)))=dφ(t)dt(−cos(φ(t))−sin(φ(t)))=−dφ(t)dt→u(t)
Kurzschreibweise: ˙→w(t)=−˙φ(t)→u(t)
<h4>Definition - Variablen Geschwindigkeit und Beschleunigung</h4>
<h5>Geschwindigkeit</h5>
als Zeit-Ableitung des Ortes nach der Zeit: →v(t):=d→r(t)dt=:˙→r(t)
<h5>Beschleunigung</h5>
als Zeit-Ableitung der Geschwindigkeit nach der Zeit: →a(t):=d→v(t)dt=:˙→v(t)=¨→r(t)
—
Herleitungen
Herleitung: Bewegungsgleichung der Massen
Aufspaltung des Ortsvektor nach Produkt von Betrag und Richtung: →r(t)=‖→r(t)‖⋅→u(t)=r(t)⋅→u(t)
Grund: bei Berechnung von Geschwindigkeit und Beschleunigung als Ableitung(en) des Ortsvektors nach der Zeit werden Betrag und Richtung (durch Produktregel) separiert: →v(t)=d→r(t)dt=ddt(r(t)⋅→u(t))=˙r(t)⋅→u(t)+r(t)⋅˙→u(t)
mit ˙→u(t)=˙φ(t)→w(t) folgt: →v(t)=˙r(t)→u(t)+r(t)˙φ(t)→w(t)
ebenso die Beschleunigung: →a(t)=˙→v(t)=ddt[˙r(t)→u(t)+r(t)˙φ(t)→w(t)]
→a(t)=ddt[˙r(t)→u(t)]+ddt[r(t)˙φ(t)→w(t)]
ddt[˙r(t)→u(t)]=¨r(t)→u(t)+˙r(t)˙→u(t) und ddt[r(t)˙φ(t)→w(t)]=˙r(t)˙φ(t)→w(t)+r(t)¨φ(t)→w(t)+r(t)˙φ(t)˙→w(t)
→a(t)=[¨r(t)→u(t)+˙r(t)˙→u(t)]+[˙r(t)˙φ(t)→w(t)+r(t)¨φ(t)→w(t)+r(t)˙φ(t)˙→w(t)]
Mit ˙→u(t)=˙φ(t)→w(t) und ˙→w(t)=−˙φ(t)→u(t) folgt:
→a(t)=[¨r(t)→u(t)+˙r(t)˙φ(t)→w(t)]+[˙r(t)˙φ(t)→w(t)+r(t)¨φ(t)→w(t)−r(t)˙φ(t)˙φ(t)→u(t)]
→a(t)=¨r(t)→u(t)+2˙r(t)˙φ(t)→w(t)+r(t)¨φ(t)→w(t)−r(t)˙φ2(t)→u(t)
→a(t)=→u(t)[¨r(t)−r(t)˙φ2(t)]+→w(t)[2˙r(t)˙φ(t)+r(t)¨φ(t)]
Trick: ddt[a2(t)˙b(t)]=2a(t)˙a(t)˙b(t)+a2(t)¨b(t)=a(t)[2˙a(t)˙b(t)+a(t)¨b(t)]
Damit: 1a(t)ddt[a2(t)˙b(t)]=a(t)a(t)[2˙a(t)˙b(t)+a(t)¨b(t)]=2˙a(t)˙b(t)+a(t)¨b(t)
Damit: [2˙r(t)˙φ(t)+r(t)¨φ(t)]=1r(t)ddt[r2(t)˙φ(t)]
Damit: →a(t)=→u(t)[¨r(t)−r(t)˙φ2(t)]+→w(t)[1r(t)ddt[r2(t)˙φ(t)]
Da Kraft und damit Beschleunigung nur radial wirkt, muss die →w(t)-Komponente Null sein:
0=1r(t)ddt[r2(t)˙φ(t)] und damit muss r2(t)˙φ(t) zeitlich konstant sein!!!