Ellipsirata

Suljetun, ellipsimäisen kiertoradan muodon, koon ja suunnan määrittämiseen tarvitaan kuusi riippumatonta muuttujaa, esimerkiksi tieto kappaleen paikka- ja nopeusvektorista tietyllä hetkellä (kaksi kolmiulotteisen avaruuden vektoria muodostaa kuusi muuttujaa). Ellipsiradan tapauksessa käytetään yleensä Johannes Keplerin elementtejä: ellipsin isoakselin puolikkaan pituus, eksentrisyys, inklinaatio eli kulma johonkin vertailutasoon verrattuna, nousevan solmun pituus, periapsin argumentti, periapsisaika.

Periapsis on kiertoradan matalin kohta ja vastaavasti apoapsis on korkein. Gravitaatiovoimien takia kappale kulkee sitä suuremmalla nopeudella mitä lähempänä se on kiertämäänsä kappaletta.

Ellipsiradan lähin piste on peri- ja kaukaisin ap-/apo-. Lisäksi kuvataan se kappale, josta peri tai ap/apo on otettu:

• Maa: perigeum ja apogeum (esimerkiksi satelliitin ellipsiradan lähin ja kaukaisin piste)
• Aurinko: periheli ja apheli (esimerkiksi planeetan Aurinkoa lähin ja kaukaisin piste)
• Tähti: periastroni ja apoastroni
• Mielivaltainen taivaankappale: periapsis ja apoapsis.

Tavallinen on tapaus, jossa hyvin pienimassainen kappale kiertää valtavan massiivista kappaletta. Tällöin pienen kappaleen vaikutus isompaan voidaan useasti unohtaa ja kirjoittaa yksinkertaisempia yhtälöitä. Kiertoaika lasketaan seuraavalla kaavalla:

${\displaystyle T={2\pi \over {\sqrt {G\cdot M}}}a^{3 \over {2}}}$,

jossa

• ${\displaystyle T}$ = kiertoaika
• ${\displaystyle a}$ = kappaleiden keskietäisyys
• ${\displaystyle G}$ = gravitaatiovakio
• ${\displaystyle M}$ = suuren kappaleen massa

Tämä yhtälö selittää Keplerin aikoinaan huomaaman Keplerin kolmannen lain, jonka mukaan planeettojen kiertoaikojen neliöt suhtautuvat toisiinsa kuin niiden isoakselin puolikkaiden kuutiot.

${\displaystyle v={\sqrt {2GM\left({1 \over {r}}-{1 \over {2a}}\right)}}}$

• ${\displaystyle v}$ = kappaleen ratanopeus
• ${\displaystyle r}$ = kappaleen etäisyys
• ${\displaystyle a}$ = kappaleen keskietäisyys
• ${\displaystyle G}$ = gravitaatiovakio
• ${\displaystyle M}$ = keskuskappaleen massa

Kahdella samamassaisella kappaleella ei voi olla samaa rataa ja eri ratanopeutta, sillä suuremman ratanopeuden omaavalla kappaleella on enemmän liike-energiaa, jolloin kappale pyrkii keskipakovoiman ansiosta poistumaan vetovoimakentästä. Tällöin rata suurenee.

Ellipsiradan soikeutta kuvataan eksentrisyydellä e, joka on ympyrälle 0. Ellipsiradalle e on aina alle 1. Planeetoille on monesti ympyrämäinen rata, jossa karkeasti e on suunnilleen 0 (ehkä e <0,1).

Planeetan etäisyys keskuskappaleesta ellipsiradalla, kun planeetan kulmaetäisyys perihelistä tunnetaan, on[1][2][3]

${\displaystyle r={\frac {a(1-e^{2})}{1+e\cos(\theta -\omega )}}}$

missä

• ${\displaystyle r}$ = planeetan etäisyys keskuskappaleesta
• ${\displaystyle a}$ = planeetan radan isoakselin puolikas
• ${\displaystyle e}$ = planeetan radan eksentrisyys
• ${\displaystyle \theta -\omega }$ = perihelin ja planeetan sijainnin välinen kulma

Planeetan kulmaetäisyys ν perihelistä ja säde ${\displaystyle r}$ saadaan laskettua, kun tunnetaan periheliaika ja planeetan nykyinen olinaika. Historiallisista syistä käytetään nelivaiheista menetelmää, jossa yhdessä tarvitaan iterointia.

1. Lasketaan planeetan keskianomalia ${\displaystyle M}$ kaavasta

${\displaystyle M={\frac {2\pi t}{P}}}$

2. Lasketaan eksentrinen anomalia ${\displaystyle E}$ numeerisesti iteroiden Kaplerin yhtälöä:

${\displaystyle \ M=E-\epsilon \cdot \sin E}$

3. Lasketaan todellinen anomalia ν yhtälöstä:

${\displaystyle \tan {\frac {\nu }{2}}={\sqrt {\frac {1+\epsilon }{1-\epsilon }}}\cdot \tan {\frac {E}{2}}}$

4. Lasketaan etäisyys auringosta eli heliosentrinen etäisyys ${\displaystyle r}$ :

${\displaystyle r={\frac {p}{1+\epsilon \cdot \cos \nu }}}$

Ja myös[4]

${\displaystyle r={\frac {k^{2}/(G(m_{1}+m_{2}))}{1+e\cos(\theta -\omega )}},}$

missä ${\displaystyle k}$ on vakio. Vakio saadaan kaavasta[5]

${\displaystyle k={\sqrt {aG(m_{1}+m_{2})|1-e^{2}|}}}$,

missä ${\displaystyle m_{1}}$ ja ${\displaystyle m_{2}}$ ovat kappaleiden massat.

Planeetan paikkavektori ${\displaystyle {\bar {r}}}$ lasketaan kaavasta[6]

${\displaystyle {\bar {r}}=a(\cos(E)-e){\bar {i}}+a{\sqrt {1-e^{2}}}\sin(E){\bar {j}}}$

eli

${\displaystyle r_{i}=a(\cos(E)-e)}$

ja

${\displaystyle r_{j}=a{\sqrt {1-e^{2}}}\sin(E)}$

ja paikkavektorin pituus eli etäisyys Auringosta on[7]

${\displaystyle \vert {\bar {r}}\vert =a(1-e\cos(E))}$

Planeetan todellinen eli luonnollinen anomalia ${\displaystyle \nu }$ lasketaan jakamalla paikkavektorin komponentit ${\displaystyle r_{j}}$ ja ${\displaystyle r_{i}}$ keskenään ja ottamalla tästä osamäärästä käänteinen tangentti[8].

${\displaystyle \nu =\arctan \left({\frac {r_{j}}{r_{i}}}\right)}$

eli silloin

${\displaystyle \nu =\arctan \left({\frac {\cos(E)-e}{\sin(E){\sqrt {1-e^{2}}}}}\right)}$

jossa

• ${\displaystyle \nu }$ = todellinen anomalia
• ${\displaystyle a}$ = radan isoakselin puolikas
• ${\displaystyle E}$ = eksentrinen anomalia
• ${\displaystyle e}$ = radan eksentrisyys
• ${\displaystyle i}$,${\displaystyle j}$ = paikkavektorin komponentit, ${\displaystyle j}$ isoakselin suuntainen

Toinen tapa lausua planeetan paikkavektori on käyttää ellipsin isoakselia ${\displaystyle {\bar {a}}}$ ja pikkuakselia ${\displaystyle {\bar {b}}}$[9]

${\displaystyle {\bar {r}}=(\cos(E)-e){\bar {a}}+\sin(E){\bar {b}}}$

Ellipsiradalla on kaltevuus eli inklinaatio ${\displaystyle i}$, joka tavallisesti ilmoitetaan asteina. Esim kaksoistähden inklinaatio ilmoitetaan taivaanpallon pintaa vastaan.

Kun kappale kiertää ellipsiradalla toista kappaletta, sen energiatasapaino vaihtelee potentiaalienergian ja liike-energian välillä. Vuoroin jompikumpi on suurempi, mutta ellipsiradan ominaisrataenergia on vakio. Jos kappale on lähellä keskuskappalettaan, sen potentiaalienergia on pieni mutta liike-energia suuri, kauempana päinvastoin. Ympyrämäisellä radalla potentiaalienergia ja liike-energia ovat tasapainossa.

${\displaystyle {v^{2} \over {2}}-{\mu \over {r}}=-{\mu \over {2a}}=\epsilon <0}$

• ${\displaystyle v\,}$ kiertävän kappaleen ratanopeus
• ${\displaystyle r\,}$ kiertävän kappaleen etäisyys keskuskappaleesta.
• ${\displaystyle a\,}$ on radan isoakselin puolikas,
• ${\displaystyle \mu \,}$ on vakiogravitaatioparametri, μ = G*M
• ${\displaystyle G}$ = gravitaatiovakio
• ${\displaystyle M}$ = massa

• Keplerin lait
• Eksentrinen anomalia
• Paraabelirata
• Hyperbelirata