<<

Appendix to Section 3: A brief overview of Geodetics

MAE 5540 - Propulsion Systems

• Navigation Geeks do Calculations in Geocentric (spherical) Coordinates

• Map Makers Give Data in Terms ! of Geodetic (elliptical) Coordinates • Need to have some idea how to relate one to another -- science of geodesy MAE 5540 - Propulsion Systems How Does the Earth Vary with ?

k z ? 2 ? 2 REarth = x + y2 + z = r + z r R " : y j

! r 2 + z 2 = 1 j 2 Req Req 1 - eEarth y x i "a" "b" r Prime ! x

MAE 5540 - Propulsion Systems How Does the vary with Latitude? r 2 + z 2 = 1 2 Req Req 1 - eEarth

! 2 2 2 2 2 1 - eEarth r + z = Req 1 - eEarth ! z 2 2 2 2 z 2 r Req = r + 2 = r 1 + 2 1 - eEarth 1 - eEarth

Rearth z " r 2 2 2 z 2 2 r = Rearth cos " r = tan "

MAE 5540 - Propulsion Systems How Does the Earth Radius vary with Latitude? 2 2 Req 2 tan ! 2 = cos ! 1 + 2 = Rearth 1 - eEarth

2 2 2 1 - eEarth cos ! + sin ! 2 = 1 - eEarth Inverting .... 2 2 2 2 2 2 cos ! + sin ! - eEarth cos ! 1 - eEarth cos ! 2 = 2 1 - eEarth 1 - eEarth R earth ! = Req 2 1 - eEarth 2 2 1 - eEarth cos !

MAE 5540 - Propulsion Systems Earth Radius vs Geocentric Latitude

R earth ! = Polar Radius: 6356.75170 km R 2 eq 1 - eEarth Equatorial Radius: 6378.13649 km 2 2 1 - eEarth cos !

2 2 2 e = 1 - b = a - b = Earth a a2

6378.13649 2 6378.136492 - = 0.08181939 6378.13649

MAE 5540 - Propulsion Systems Earth Radius vs Geocentric Latitude (concluded)

6380.

6375.

6370. Radius, Km 6365.

6360.

6355. -100 -50 0 50 100 Geocentric Latitude, deg. MAE 5540 - Propulsion Systems Earth Radius … alternate formula

• Earth radius as Function of Latitude

MAE 5540 - Propulsion Systems What is the mean radius of the earth?

RE = Req 2 ! 2 dA = " [ RE cos ( !) ] 1 - eEarth ! 2 2 1 - eEarth cos !

RE cos(!) • IAU Convention: ! Based on Earth's RE d! cos (!) ! RE ! Volume: !

3 4 " R = V 3 Emean E

2 RE cos(!) x RE d! cos (!) dV = " [ ] !

MAE 5540 - Propulsion Systems What is the Earth's Mean Radius? (continued) • Earth's () Volume

! 2 V = R cos 3d = E ! E" " " -! 2

! 2 1 - e2 3/2 R3 Earth cos3 d = eq ! 2 2 " " 1 - eEarth cos " -! 2

3 4 ! 2 R 1- e eq 3 earth

MAE 5540 - Propulsion Systems What is the Earth's Mean Radius? (continued) • Based on Volume 3 Ellipsoid Volume: 4 ! 1- e2 R 3 eq

3 Sphere Volume: 4 ! R 3 sphere " R R = 1- 0.081819392 1/6 6378.13649 = 6371.0002 km sphere # mean • Mean Radius we have been using is for a Sphere with same volume as the Earth ME = $E VE "gravitational radius"

MAE 5540 - Propulsion Systems Earth Radius vs Geocentric Latitude (concluded)

“Gaussian Surface” 6380.

6375.

6370. Radius, Km 6365.

6360.

6355. -100 -50 0 50 100 Geocentric Latitude, deg. MAE 5540 - Propulsion Systems Geocentric vs Geodetic Coordinates

Mean North (Celestial) Pole

• Map makers Geocentric vs Geodetic C2 ooordinates define a new Req 1 - e(earth) latitude which is Mean Greenwich the angle that Meridian RE normal to the ! Earth's surface ! makes with the Geocenter R respect to the !' eq equatorial plane " Req • Geodetic latitude

MAE 5540 - Propulsion Systems Geocentric vs Geodetic Coordinates (contined)

Mean •Since the Earth is North (Celestial) Elliptical only along the Pole z-axis ... geodetic and geocentric 2 h are Req 1 - e(earth) identical Mean Greenwich Meridian RE • Altitude is an ! extension of the Geocenter ! line of latitude Req geodetic) !' " Req Equator MAE 5540 - Propulsion Systems Geocentric vs Geodetic Coordinates (contined)

•Complex nonlinear Mean equations describe North (Celestial) relationship between Pole geocentric and geodetic latitude 2 h Req 1 - e(earth) Mean • Derivation requires Greenwich Meridian RE Extensive ! Knowledge of Spherical ! Trigonometry Geocenter R !' eq " Req Equator MAE 5540 - Propulsion Systems Geocentric vs Geodetic Coordinates (contined)

Geocentric Cartesian Coordinates * rtarget xtarget = R' +h cos !' cos " !' h

y = R' +h cos ' sin z target !' ! "

2 !' ztarget = R' 1 - e +h sin !' ! !' earth R' ' Req ! R' = !' 2 1 - e sin2 !' !' earth rtargetr

* We would be here all week if I try to derive this "Radius of " Req MWAhEit m55o4re0, S- tPeprohepnu Als.i, oand S Hyasetreinmgs, Edward A., Jr., FORTRAN Program for the Analysis of Ground Based Range Tracking Data--Usage and Derivations, NASA TM 104201, December, 1992 Geocentric vs Geodetic Coordinates (contined)

Geocentric Polar Coordinates rtarget R = x 2 + y 2 + z 2 target target target target h

z -1 ytarget target = tan " xtarget

!' ! z -1 target R' !target = tan 2 2 !' xtarget + ytarget !' rtargetr

"" Req MAE 5540 - Propulsion Systems Geocentric vs Geodetic Coordinates (concluded)

Inverse Relationships, non-linear no direct solution

2 2 xtarget + ytarget h = - R' cos !' !' y = tan-1 target "target xtarget

ztarget 1 ' = -1 # ! target tan 2 2 R' xtarget + ytarget 1 - e2 !' earth R' + h !' target MWAhEit m55o4re0, S- tPeprohepnu Als.i, oand S Hyasetreinmgs, Edward A., Jr., FORTRAN Program for the Analysis of Ground Based Range Tracking Data--Usage and Derivations, NASA TM 104201, December, 1992 Pulling it all together • Given geodetic coordinates -- compute geocentric i) Compute geocedtric cartesian coordinates Range, runway threshholds, radar antennae, beacon ....

x = R' +h cos !' cos " target !'

y = R' +h cos !' sin " target !'

2 ztarget = R' 1 - e +h sin !' !' earth Req R' = !' 2 1 - e sin2 !' MAE 5540 - Propulsion Systems earth Pulling it all together (continued)

• Given geodetic coordinates -- compute geocentric ii) Compute Geocentric polar coordinates next

2 2 2 Rtarget = xtarget + ytarget + ztarget

-1 ytarget target = tan ! xtarget

-1 ztarget "target = tan 2 2 xtarget + ytarget

MAE 5540 - Propulsion Systems Pulling it all together (concluded) • Given geocentric (usually x,y,z) coordinates -- compute geodetic Req R' = ' 2 ! 1 - e sin2 !' GPS, INS, TLE's .... earth

No explicit solution: 2 2 xtarget + ytarget requires h = - R' 1) series expansion solution, cos !' !' 2) numerical iteration, -1 ytarget 3) or a special solution "target = tan called "Ferrari's method"** xtarget

z -1 target 1 !'target = tan # x 2 + y 2 R' target target 1 - e2 !' earth R' + h !' target **NASA Technical Paper 3430, Whitmore and Haering, FORTRAN Program forAnalyzing Ground-Based Tracking Data: Usage and Derivations, Version 6.2, 1995 MAE 5540 - Propulsion Systems Numerical Example

• Edwards Air Force Base, Radar Site #34 Numerical example !' = 34.96081°

" = –117.91150° h = 2563.200 ft

• Find corresponding geocentric cartesian and polar coordinates

MAE 5540 - Propulsion Systems Numerical Example (cont'd)

• Compute Local Radius of Curvature

Req R' = = !' 2 1 - e sin2 !' earth

6378.13649 km = 1 - 0.08181939 sin 34.96081 " # 2 180

6392.187109 km

MAE 5540 - Propulsion Systems Numerical Example (cont'd)

• Compute X and Y (geocentric) r = R' +h cos ' = target !' !

6392.1871 + 2536.2 " 3. 048 " 10-4 cos 34.96081 " # = 180

5239.3131 km

xtarget = rtarget cos $ = ytarget = rtarget sin $ =

5239.3131 km " -117.91150 " # = 5216.0074 km " sin -117.91150 " # = cos 180 180

-2452.5602 km -4629.83218 km

MAE 5540 - Propulsion Systems Numerical Example (cont'd)

• Compute z (geocentric)

2 ztarget = R' 1 - e +h sin !' = !' earth

6392.1871 km 0.081819392 2536.2 " 3. 048 " 10-4 sin 34.96081 " # = 1 - + 180

3638.7480 km

MAE 5540 - Propulsion Systems Numerical Example (cont'd) • Compute Geocentric Polar Coordinates

2 2 2 Rtarget = 2452.5602 + 4629.83218 + 3638.7480 =

6378.94104km y = tan-1 target = !target xtarget

180 # -1 -4629.83218 = -117.9115° " tan -2452.5602 z -1 target = $target = tan 2 2 xtarget + ytarget

180 # -1 3638.7480 = 34.7803° " tan 5239.3131 MAE 5540 - Propulsion Systems Numerical Example (cont'd) • Compute Local Earth Radius and Geocentric Above

RE = Req = ! 2 1 - eEarth 2 2 1 - eEarth cos !

6378.13649 km = 2 1 - 0.08181939 0.081819392 2 34.7083 " # 1 - cos 180

6364.23 km D = R - R = geoid target E!

6378.94104 - 6364.23 km = 14.7 km = 48228.25 ft

MAE 5540 - Propulsion Systems Numerical Example (concluded)

• Comparison

Geodetic Geocentric !' = 34.96081° ! = 34.7803°

" = –117.91150° " = –117.91150°

h = 2563.200 ft Dgeoid = 48228.25 ft

• Earth Oblateness is NOT trivial, and in the REAL World -- it must be accounted for

MAE 5540 - Propulsion Systems Appendix II: Rigorous Derivation of Realizable Launch Inclination

• Launch Initial Conditions Greenwich Meridian • Position: ", Latitude z #, Longitude (Inertial) Launch meridian h, Altitude

Equator

y Re

"

x

Vernal Equinox !g

95 MAE 5540 - Propulsion Systems • Launch Initial Conditions

North Pole Greenwich Meridian z Launch meridian • Launch Initial Conditions (Inertial Coordinates) Equator

y Re

"

x

Vernal Equinox !g Sidereal hour angle 96 MAE 5540 - Propulsion Systems Sidereal Hour Angle

Ignore for now!

97 MAE 5540 - Propulsion Systems Computing the Hour Angle

98 MAE 5540 - Propulsion Systems