Orbit Characteristics

Similar documents
Lecture 22: Gravitational Orbits

Lecture 15 - Orbit Problems

Central force motion/kepler problem. 1 Reducing 2-body motion to effective 1-body, that too with 2 d.o.f and 1st order differential equations

F = ma. G mm r 2. S center

Astromechanics. 10. The Kepler Problem

MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. Principles of Space Systems Design

Celestial Mechanics II. Orbital energy and angular momentum Elliptic, parabolic and hyperbolic orbits Position in the orbit versus time

L03: Kepler problem & Hamiltonian dynamics

MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. Principles of Space Systems Design

Spacecraft Dynamics and Control

!! r r θ! 2 r!φ 2 sinθ = GM / r 2 r!! θ + 2!r θ! r!φ 2 sinθ cosθ = 0 r!! φ sinθ + 2r θ!!φ cosθ + 2!r!φ sinθ = 0

Lecture D30 - Orbit Transfers

Universal Gravitation

Chapter 8. Orbits. 8.1 Conics

Motion under the Influence of a Central Force

Astromechanics. 6. Changing Orbits

Use conserved quantities to reduce number of variables and the equation of motion (EOM)

Celestial Mechanics I. Introduction Kepler s Laws

We start by noting Kepler's observations (before Newton) regarding planetary motion.

Chapter 13: universal gravitation

Physics 115/242 The Kepler Problem

A. Correct! These are the corresponding rectangular coordinates.

ASE 366K Spacecraft Dynamics

Orbital Mechanics! Space System Design, MAE 342, Princeton University! Robert Stengel

Flight and Orbital Mechanics

Satellite Communications

Orbital Mechanics MARYLAND

A Unified Method of Generating Conic Sections,

Radial Acceleration. recall, the direction of the instantaneous velocity vector is tangential to the trajectory

ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS

Conic Sections in Polar Coordinates

OPTIMAL MANEUVERS IN THREE-DIMENSIONAL SWING-BY TRAJECTORIES

Previous Lecture. Orbital maneuvers: general framework. Single-impulse maneuver: compatibility conditions

Extending the Patched-Conic Approximation to the Restricted Four-Body Problem

Chapter 13. Gravitation

Chapter 13. Universal Gravitation

1 Summary of Chapter 2

REVIEW OF KEY CONCEPTS

Spacecraft Dynamics and Control

Lecture 1: Oscillatory motions in the restricted three body problem

GRAVITATION. F = GmM R 2

History. Geocentric model (Ptolemy) Heliocentric model (Aristarchus of Samos)

Celestial Mechanics Lecture 10

ISIMA lectures on celestial mechanics. 1

Teaching Mathematics in Twenty First Century: An Integrated Approach

Gravitation. chapter 9

Copyright 2009, August E. Evrard.

Explanation: The escape velocity and the orbital velocity for a satellite are given by

Orbital Mechanics MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design

AST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements

Two-Body Problem. Central Potential. 1D Motion

Lambert's Problem: given r 0, r 1, t; solve for the conic parameters.

The Heliocentric Model of Copernicus

Celestial Orbits. Adrienne Carter Ottopaskal Rice May 18, 2001

VISUAL PHYSICS ONLINE

Honors Precalculus Chapter 8 Summary Conic Sections- Parabola

(b) The period T and the angular frequency ω of uniform rotation are related to the cyclic frequency f as. , ω = 2πf =

Seminar 3! Precursors to Space Flight! Orbital Motion!

Notes on Planetary Motion

Orbital Mechanics MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design

Lecture 8. Kepler's IInd: Angular Momentum

Satellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN

Chapter 13. Gravitation

Precalculus Conic Sections Unit 6. Parabolas. Label the parts: Focus Vertex Axis of symmetry Focal Diameter Directrix

Patch Conics. Basic Approach

Distance and Midpoint Formula 7.1

Gravitation and the Waltz of the Planets

Gravitation and the Waltz of the Planets. Chapter Four

TWO KINDS Evenly spaced parallel lines measure distance from a line.

Lecture Outline. Chapter 13 Gravity Pearson Education, Inc. Slide 13-1

Intermediate Math Circles Wednesday, April 5, 2017 Problem Set 8

10 Orbit and Constellation Design Selecting the Right Orbit

Gravitation and the Motion of the Planets

PHY321 Homework Set 10

9-4 Ellipses. Write an equation of each ellipse. 1. ANSWER: ANSWER:

Newton s Gravitational Law

A space probe to Jupiter

Astrodynamics (AERO0024)

Gravitation & Kepler s Laws

Copyright 2008 Pearson Education, Inc., publishing as Pearson Addison-Wesley.

APPLICATION OF THE HÉNON S ORBIT TRANSFER PROBLEM TO MANEUVER A SATELLITE IN A CONSTELLATION

A = 6561 times greater. B. 81 times greater. C. equally strong. D. 1/81 as great. E. (1/81) 2 = 1/6561 as great Pearson Education, Inc.

Chapter 3. Algorithm for Lambert's Problem

Astrodynamics (AERO0024)

Gravitation and Central-force motion

A Study of the Close Approach Between a Planet and a Cloud of Particles

Chapter 4. Integrated Algorithm for Impact Lunar Transfer Trajectories. using Pseudo state Technique

Central Force Problem

MINIMUM IMPULSE TRANSFERS TO ROTATE THE LINE OF APSIDES

AST111 PROBLEM SET 2 SOLUTIONS. RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch J2000).

Orbital Mechanics MARYLAND. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design

Theorem of the Keplerian kinematics

Introduction to Computer Graphics (Lecture No 07) Ellipse and Other Curves

Satellite meteorology

18. Kepler as a young man became the assistant to A) Nicolaus Copernicus. B) Ptolemy. C) Tycho Brahe. D) Sir Isaac Newton.

Stellar Dynamics and Structure of Galaxies

Course Overview/Orbital Mechanics

Physics Lecture 03: FRI 29 AUG

Lecture 1a: Satellite Orbits

Orbital Mechanics MARYLAND. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design

Transcription:

Orbit Characteristics We have shown that the in the two body problem, the orbit of the satellite about the primary (or vice-versa) is a conic section, with the primary located at the focus of the conic section. Hence the orbit is either an ellipse, parabola, or a hyperbola, depending on the eccentricity. For conic sections we have the following classifications: e < 1 ellipse e = 1 parabola (An exception this relation is that all rectilinear orbits have e = 1, and angular momentum = 0) e > 1 hyperbola Of main interest for Earth centered satellites (Geocentric satellites) and Sun centered satellites (Heliocentric satellites) are elliptic orbits. However when we go from one regime to another such as leaving the Earth and entering into an interplanetary orbit then we mus deal with hyperbolic orbits. Or if we approach a planet form a heliocentric orbit, hyperbolic orbits are of interest. Parabolic orbits, on the other hand are more theoretical than practical and simply define the boundary between those orbits which are periodic and hang around, (elliptic orbits), and those orbits which allow one to escape from the system (hyperbolic orbits). So one might say that the parabolic orbit is the minimum energy orbit that allows escape. In the following, we will determine the properties of each of these types of orbits and write some equations that are applicable to all orbits and other equations that are applicable only to the type of orbit of interest. It was indicated earlier that the constant e (eccentricity) was related to some other constants (angular momentum and energy that we developed earlier). We can now determine this relation by using the expression we found for r and substitute it into the energy equation. The energy equation is given by V & µ r En 0r % r 0 & µ r En (1) However, since the energy is a constant over the whole orbit, we can evaluate the energy equation a any convenient point we choose. If you pick an arbitrary point in the orbit, the result will be the same. To reduce the algebra required, we can select the periapse location and evaluate Eq. 1 at periapse. At periapse, the radial component of velocity is zero, hence we eliminate a term. Also a periapse, we can note that r becomes, h r r min /µ 1 % e cos< /0 < 0 h /µ 1% e () r From the angular momentum relation 0 h, we can replace 0 with h/r. Equation 1

becomes, h r 4 & µ r En (3) We can replace r with Eq. (r at periapse) to get h (1 % e ) (h /µ ) & µ(1& e) h /µ En (4) If we simplify we end up with the desired result, e 1 % h En µ (5) Equation 5 relates the eccentricity of the orbit to its angular momentum and energy. It is valid for any type of orbit. It is clear from Eq. 5 that the following classification can occur. En < 0 e < 1 orbit is ellipse En = 0 e = 1 orbit is parabola En > 0 e > 1 orbit is a hyperbola These are important relations that should be committed to memory. Eq. 5 is also an important equation to know. Summary of General Results for All Orbits (elliptic, parabolic, and hyperbolic) The Equation of the orbit: r(<) h /µ 1 % e cos < p 1 % ecos < (6) where p = h /µ = Orbit Parameter r = radial distance from center of attracting body (focus) e = eccentricity < = true anomaly (measured from point of closest approach (apoapse) Energy Equation:

& µ r V r % V V & µ r En (7) where: V r 0r, V r 0< Angular Momentum: h r 0< rv rvcosn const (8) where N (or () is the flight path angle ( angle velocity vector makes with local horizontal). We can develop other relations as follows. Radial Velocity in Terms of True Anomaly The radial velocity is just 0r, which can be obtained by differentiating Eq. 6 with respect to time. pesin< 0< r 0< esin< 0r (9) (1 % e cos< ) 1 % e cos< By multiplying and dividing by p = h /µ we can rewrite Eq. 9 as, 0r 1 p r 0< e sin< µ p e sin< (10) From the basic definition of the flight path angle and from Eq. 9 we have, Tan N V r V 0r r 0< e sin< 1 % ecos< (11) We can also define the semi-latus rectum as the value of r when < = B/. Hence the semi-latus rectum = p = h /. Finally we can note that since cos < = cos -<, that the orbit is symmetric about the point of closest approach ( or about the major axis of the conic section).

PROPERTIES OF THE ORBITS Parabolic Orbit (e = 1, En = 0) The parabolic orbit serves as a boundary between the elliptic (periodic) orbits and the hyperbolic (escape) orbits. It is the orbit of least energy that allows escape. The orbit equation becomes, r(<) h /µ 1 % cos< p 1 % cos< (1) Further the periapse distance = r(0) = p/ = h / µ. The energy equation becomes. V & µ r 0 Y V escape µ r (13) Equation 13 defines the escape velocity, which is the minimum velocity to escape the two body system at the given radius r. Note that the speed required for escape is independent of it direction! The flight path angle in a parabolic orbit is given by, tan N or sin< 1 % cos< sin</ cos</ 1 % cos </ & sin </ tan</ (14) N < 0< cosn µr r Also it is easy to show ( use V cos N = r ) that p. r Elliptic Orbits (e < 1, En, 0) Recall that the radius position is measured from the focus of the ellipse (not the center). By evaluating the orbit equation at values of the true anomaly of 0 and B, we can determine the closest approach (periapse distance, ) and the furthest approach (apoapse distance, r a). The sum of these two comprise the major axis of the ellipse. Of interest to us is half the distance or the semi-major axis, usually designated by the symbol a. Hence a is the distance from the center of the ellipse to the periapse and apoapse. We can set < = 0 and B to obtain and r a respectively, h

r(0) r min h /µ 1%e, r(b) r max r a h /µ 1 & e (15) Then the major axis, a = r a +, and we have, a h /µ 1 1 % e % 1 1 & e h /µ 1 & e (16) or h /µ a (1& e ). Hence the orbit equation for the ellipse becomes; r(<) a (1& e ) 1 % e cos< (17) Also it follows: r = a ( 1 - e ) p r = a ( 1 + e ) a a - r = a e p Energy in an Elliptic Orbit e r a & r a % ba 1 & e (semi-minor axis) We can eliminate the eccentricity from Eq. 16b using Eq. 5, h /µ a (1& e ) a 1 & 1 % h En & ah En µ µ (18) or En & µ a (19) Hence the energy equation takes the form (for an elliptic orbit) V & µ r &µ a (0)

Additionally we have the flight path angle relation tan N e sin< 1 % e cos< (1) and in terms of r, cos N r 0< V a (1& e ) r & r a () An interesting set of relations happen at the point on the orbit at the end of the semiminor axis. At that point, the following relations can be found: r = a, cos < = -e, tan N = e, sin N = e. 1 & e Period of an Elliptic Orbit An elliptic orbit is the only type of orbit that has a period. We can determine the period by recalling the angular momentum equation and noting that half the angular momentum is the aerial rate, da dt 1 r 0< h (3) Then integrating both sides over one period or once around the orbit we have A T p da m m 0 0 h dt h T p B ab (area of ellipse) (4) But h µ a (1& e ), and b a (1& e ), and the (1 - e ) terms cancel leading to the period given as T p B a 3 µ (5)

We can see that both the energy and the period of the orbit depend only on the size of the orbit and not on the shape (e). It is useful to define the mean angular rate, n = B / T. With this definition we can write p a form of Kepler s third law, the square of the period of an orbit is proportional to the cube of the orbit size. n a 3 µ (6) Properties of Elliptic Orbits in Term s of Periapse and Apoapse Distances From geometry we noted previously that a = r a +. Then the energy equation can be written as V & µ r & µ r a % En (7) If we evaluate Eq. 7 at periapse and apoapse, we can determine the velocities at these points. The results, after very little algebra, are, V p µ r a 1 % r a V a µ r a 1 % r a (8) Since h = r V = r V we can calculate h in terms of r and r, that gives p p a a p a h µ r a 1 % r a µ ra 1 % r a µ r a % r a (9) Circular Orbit - As a Special Case of Elliptic Orbit (e = 0, En & µ ) r c With e = 0, the orbit equation gives, r p h /µ a const r c (30)

From the energy equation, V c & µ r c & µ r c Y V c µ r c (31) V c is called the circular speed and is defined at every radius r as orbit. V c µ r regardless of the Hyperbolic Orbits (e > 1, En > 0) From the orbit equation, it is clear that the radial distance will go to infinity when the denominator term goes to zero. Then r 4 when ( 1 + e cos < ) = 0 < min max cos &1 & 1 e (3) Equation 3 puts limits on the true anomaly of an hyperbolic orbit. Just as in an elliptic orbit, we can calculate the closest approach at < = 0. We can also calculate the apoapse distance by letting < = B. However the result is negative and represents the closest approach of the other branch of the hyperbola, one that has no meaning in our orbit. However formally we can then note that h /µ 1 & e, r B h /µ 1 % e &( % a) (33) Substituting for r in the above equation yields the hyperbolic orbit result, p h /µ a ( e & 1) (34) and the corresponding hyperbolic orbit equation r(<) a (e & 1) 1 % e cos < (35) e 1 % h En If we substitute for in Eq. 34 we find the hyperbolic orbit energy, µ

En µ a (36) The corresponding energy equation is V & µ r µ a (37) One of the consequences of Eq. 37 is that at infinity, the velocity is no longer zero and is given by V 4 µ a (38) and is defined as the hyperbolic excess velocity. We can also determine the flight path angle in terms of < or in terms of r in a similar manner as we did for elliptic orbits. The results are tan N esin < 1 % e cos < (39) and cos N h rv a ( e & 1) r % r a (40) ` Finally, the hyperbolic orbit has a property that no other orbit has. A vehicle traveling the length of the orbit will arrive coming from some point at infinity, and then fly by through the closest approach point (apoapse), and then leave, going to some point at infinity. The approach direction comes in from a direction of -< 4 and leave in the direction of +< 4. The angle through which the vehicle turns is called the turning angle, *. This turning angle can be determined from the properties of the hyperbola. * < 4 & B Y * < 4 & B (41)

Then sin * sin < 4 cos B & cos < 4 sin B &cos < 4 1. Consequently the turning angle is e given by sin * 1 e (4) The properties characteristic of all orbits are presented in this section. Equations which apply to all orbits are given and those which are applicable to specific orbits are presented.