Loading [MathJax]/jax/output/CommonHTML/jax.js

User Tools

Site Tools


module:mathematik:nkoerperproblem

This is an old revision of the document!


N-Körper-Problem

Quellen

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=0a=0dvdt=0v=const

(2.) F=ma

(3.) Fab=Fba

Gravitationsgesetz

Fm1m2=Gm1m2r3m1m2rm1m2

Fm1m2=Gm1m2r2m1m2erm1m2

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) , dLdt=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 : (xy)=(+cos(α)sin(α)+sin(α)+cos(α))(xy)

Rücktransformation um α von CCSCCS : (xy)=(+cos(α)+sin(α)sin(α)+cos(α))(xy)

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:

du(t)dt=ddt(cos(φ(t))sin(φ(t)))=(dφ(t)dtsin(φ(t))+dφ(t)dtcos(φ(t)))=dφ(t)dt(sin(φ(t))+cos(φ(t)))=dφ(t)dtw(t)

Kurzschreibweise: ˙u(t)=˙φ(t)w(t)

<h5>Einheitsvektor Geschwindigkeit</h5>

Einheitsvektor Geschwindigkeit: w(t)=(wxwy)

Drehung u(t) um +π2:u(t)\rarrw(t)

(wxwy)=(+cos(π2)sin(π2)+sin(π2)+cos(π2))(uxuy)

(wxwy)=(01+10)(uxuy)

(wxwy)=(01+10)(cos(φ(t))sin(φ(t)))

w(t)=(sin(φ(t))cos(φ(t)))

Ableitungen des Geschwindigkeits-Einheitsvektors w(t) nach der Zeit:

dw(t)dt=ddt(sin(φ(t))+cos(φ(t)))=(dφ(t)dtcos(φ(t))dφ(t)dtsin(φ(t)))=dφ(t)dt(cos(φ(t))sin(φ(t)))=dφ(t)dtu(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):=dr(t)dt=:˙r(t)

<h5>Beschleunigung</h5>

als Zeit-Ableitung der Geschwindigkeit nach der Zeit: a(t):=dv(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)=dr(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!!!

Ergebnisse

module/mathematik/nkoerperproblem.1638902670.txt.gz · Last modified: 2021/12/07 20:44 (external edit)