Saturday, September 22, 2012

Reference: 13 Steps to Compute Classic Orbital Elements

13 Steps to Compute Classic Orbital Elements
Seq.NameSymbolEquation
1.Distance
r
2.Speed
v
3.Radial velocity
vr
4.Specific angular momentum
h
5.Magnitude of specific angular momentum
6.Inclinationi
7.Node line
N
8.Magnitude of N
9.Right Ascension (RA) of ascending node
Ω
10.Eccentricity vector
11.Eccentricity
12.Argument of perihelion
13.True anomaly.
ν
To compute the six classic orbital elements, some experts suggest 13 steps as described in textbook published by Embry Riddle University, Space Mechanics by Howard D. Curtis, pages 6.10-12, Sep 2003.

Auxilliary Circle is the circumcircle of an orbit. The circle's centercoincides with the orbit's center, and the circle's radius equals the orbit's semimajor axis.

EccentricPosition in OrbitDist fm SunTrueMean
AnomalyX-valueY-valueR-valueAnomalyAnomaly
(Deg)(AU)(AU)(AU)(Deg)(Deg)
01.470.000.6500.00.0
151.420.320.67827.86.7
301.270.610.76053.414.0
451.040.860.89075.722.3
600.741.061.06094.632.2
750.381.181.258110.544.0
900.001.221.470123.957.9
105-0.381.181.682135.574.0
120-0.741.061.880145.892.2
135-1.040.862.050155.1112.3
150-1.270.612.180163.8134.0
165-1.420.322.262172.0156.7
180-1.470.002.290180.0180.0
195-1.42-0.322.262188.0203.3
210-1.27-0.612.180196.2226.0
225-1.04-0.862.050204.9247.7
240-0.74-1.061.880214.2267.8
255-0.38-1.181.682224.5286.0
2700.00-1.221.470236.1302.1
2850.38-1.181.258249.5316.0
3000.74-1.061.060265.4327.8
3151.04-0.860.890284.3337.7
3301.27-0.610.760306.6346.0
3451.42-0.320.678332.2353.3
3601.470.000.650360.0360.0
Givena×cos(E)b×sin(E)SQRT((x - c)2 + y2)cos-1((x-c)/r)E - ν
TRUE ANOMALY: Θ = 360 - cos-1((x-c)/r) after E>180°
  • E=Eccentric Anomaly: Angle originates from elliptical center and starts with line to perihelion.
  • a=semi-major axis. For Apollo, a = 1.47 AU
  • b=semi-minor axis. For Apollo, b = 1.22 AU
  • c=focal length. For Apollo, c = 0.82 AU, distance of Sun from orbit center.
  • ν=True Anomaly: Angle originates from Sun and starts with line to perihelion.
Source material: Fundamentals of Astrodynamics, by Roger R. Bate, Donald D. Mueller, Jerry E. White.  Observing enclosed figures described above, determine areas of certain figures.

ECCENTRIC ANOMALY. If we assume this figure to be a circular sector defined by points, C, P, and ACIR. Points are shown in figure in above table and defined below.
  1. C is center of Auxiliary Circle which coincides with orbital center.
  2. P is perihelion of orbit.
  3. ACIR is Apollo's position superimposed upon Auxiliary Circle as described in preceding text.
  4. This area is a circular sector and easily computed.

TRUE ANOMALY. The outline of this figure vaguely resembles previous figure but with some significant differences. The three points which define the border of this area are shown in above figure and described below.
  1. S, Sun, is a focus of the orbit and the apex of the angle ASP.
  2. P is perihelion of orbit.
  3. AORB, Apollo's position in orbit.
  4. This area is not so easy to compute.
    • The enclosing arc is a portion of an ellipse, thus not circular.
    • The angle apex is at a focus, not at the center.
Thus, computing the enclosed area seems far from simple. However, Kepler discovered a straightforward way to do it.

MEAN ANOMALY. Observing the above figure which shows both of above, it seems the area of eccentric anomaly must always contain area bounded by true anomaly. Intuitively, it seems there must be a circular sector with area less then that of eccentric anomaly which equals area of the true anomaly figure, a non-circular sector. Once we determine this circular sector, then computing its area would be a simple task.
In fact, Kepler discovered this circular sector and named it, Mean Anomaly. It's computed as follows.
  • M = E - e sin E
  • M = mean anomaly (radians)
  • E = eccentric anomaly (radians)
  • e = eccentricity
Purpose of determining enclosed areas: Recall Kepler's 2nd Law of orbital motion.

MEAN MOTION. As the positional vector, r, sweeps out areas during Apollo's orbit, areas per time can be related to motion of Apollo throughout it's elliptical orbit.


0 Comments:

Post a Comment

<< Home