Go Back   Orbiter-Forum > Articles > A3. Orbiter Tutorials
Register Blogs Orbinauts List Social Groups FAQ Projects Mark Forums Read

Article Tools
george7378's Avatar
george7378 george7378 is offline

Default So, you want to calculate orbits?
by george7378 03-04-2012, 02:09 PM

Well, do you? If so, I hope this can show you how to calculate and therefore simulate an orbit yourself! By this I mean you will be able to calculate a satellite's altitude and velocity at a chosen time, among other things, like OrbitMFD shows. Let's assume we're orbiting Earth, and that there are no other forces acting on the satellite other than the Earth's gravity. So, what do we do with this assumption?

1. Decide on an apogee and perigee.

Let's have a nice sensible 400 x 1400 km orbit.

2. Calculate some constants for the orbit.

These constants are a)Apogee and perigee radii b) Semi-major axis, c) Eccentricity, d) Mean motion and e) Orbital period.

a) Apogee and perigee radii are simply the apogee and perigee + the radius of Earth:

R_a = apogee+6378 = 7778km
R_p = perigee+6378 = 6778km

b) Semi-major axis is the distance from the centre of your orbit to one of the extreme ends:

SMa = \frac{R_a + R_p} {2} = 7278km

c) Eccentricity is how close (or far away, depending on how you look at it) your orbit is to being a circle:

e = \frac{R_a - R_p}{R_a + R_p} = 0.069

d) Mean motion is a measure of how quickly a satellite moves in its orbit (The constant GM is the gravitational constant multiplied by Earth's mass, and equals 398600.4):

M_m = \sqrt{\frac{GM}{SMa^3}} = 1.02 x 10^-3

e) Orbital period is how long it takes to complete one orbit:

T = 2\pi\sqrt{\frac{SMa^3}{GM}} = 6179 seconds

3. Decide on a time for which you want to calculate the altitude and velocity.

By 'a time', I mean 'how far into the orbit do you want to be?' - eg: if t=0, you are at perigee. If t=T/2, you are halfway round the orbit and at apogee. If you choose anything in between, you will calculate for that position. So this basically means 'how many seconds past the perigee do you want to be?' - you can't have a negative number, or a number larger than the orbital period.

We will calculate for the apogee. Therefore:

t = T/2 = 3089.5 seconds.

4. Calculate some anomalies (angles).

These are simply angles between different points in the orbit. They are a) Mean anomaly, b) Eccentric anomaly and c) True anomaly. To understand them, it's best to see a diagram like this one:

These angles are in radians!

a) Mean anomaly:

M = M_m \times t = 3.15

b) To find the Eccentric anomaly, you have to iterate. The formula you need to iterate is:

E = E_0 - \frac{E_0-esin(E_0)-M}{1-ecos(E_0)}

To start, make E_0 = M

If you're not faimilar with iteration, it means that you have to take an initial value (in this case, M, your Mean anomaly) and substitute it into the given equation. Then, you replace your initial value with your result of your first calculation, and do the same thing again and again until the result you get doesn't change. In this case, because you are at apogee, E = M, so you get:

E = 3.15

c) True anomaly:

v = cos^{-1}(\frac{cos(E)-e}{1-ecos(E)}) = 3.15

As you can see, all three anomalies are equal at this point, meaning that our method is working (since we are calculating for the apogee where they all ~ 3.15 radians. The actual value is 3.141 radians, but the reason we are slightly out is because of rounding in our previous answers. No biggie at all).

5. Calculate the final result!

To arrive at the final result of altitude and velocity, you need one more interim value. The flightpath:

F_p = tan^{-1}(\frac{esin(v)}{1+ecos(v)}) = -6.2 x 10^-4

A little note: This number should be zero at apogee, but again, it's because of rounding. This is not a big deal, and will still give us a perfectly accurate answer.

Finally, we can now calculate our results!

alt = (SMa \times 1-ecos(E))-6378 = 1402km.

vel = \sqrt{ GM( \frac{2}{alt+6378}-\frac{1}{SMa})} = 6.91km/s

As you can see, there are minor errors of ~ 2km in altitude, and hence a tiny error in velocity, but this is all because of rounding in the constants and prior results. If you use a program to solve these equations, you will get far more accurate values.

It is of course inconvenient to solve these equations manually, but if used as part of a program which continues to solve as time increases, a very accurate orbit simulation can be made!

I have created a program which runs on these equations alone, and which simulates an Earth-orbiting satellite in the way we just have. It's called Orbit-Vis, and is available here:


Thanks for reading!
Reply With Quote
Views 31572 Comments 0
Total Comments 0



  Orbiter-Forum > Articles > A3. Orbiter Tutorials

Article Tools

Posting Rules
BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts
Forum Jump

All times are GMT. The time now is 01:36 AM.

Quick Links Need Help?

About Us | Rules & Guidelines | TOS Policy | Privacy Policy

Orbiter-Forum is hosted at Orbithangar.com
Powered by vBulletin® Version 3.8.6
Copyright ©2000 - 2017, Jelsoft Enterprises Ltd.
Copyright 2007 - 2012, Orbiter-Forum.com. All rights reserved.