DOCSLIB.ORG

The Perturbations of a Hyperbolic Orbit by an Oblate Planet

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
The Perturbations of a Hyperbolic Orbit by an Oblate Planet Technical Report No. 32-731 The Perturbations of a Hyperbolic Orbit by an Oblate Planet Carl G. Sauer, Jr. OTS PRICE I XEROX MICROFILM $ JET PROPULSION LABORATORY C-AtlFORNiA 1NSTlTUTE OF TECHNOLOGY PASADENA.CALIFORNIA January 15, 1963 NATIONAL AERONAUTICSAND SPACE ADMINISTRATION CONTRACTNo. NAS 7-100 Tecbnical Report No. 32-731 The Perturbations of a Hyperbolic Orbit by an Oblate Planet Carl G. Sauer, Jr. C. Gates, Chief Systems Analysis Section JET PROPULSION LABORATORY CALIFORNIA INSTITUTE OF TECHNOLOGY PASADENA,CALIFORNIA January 15, 1963 Copyright @ 1963 Jet Propulsion laboratory California Institute of Technology The material in this Report appeared in the ARS Journal, Vol. 32, No. 5, pp. 714-717. JPL TECHNICAL REPORT NO. 32-131 CONTENTS 1. Introduction ..................... 1 II. The Perturbative Potential ................2 111. Variation of the Orbital Elements ............. 2 IV. Example of an Earth-Escape Mission ............4 References ....................... 8 1. Example of an Earth-escape mission ..............5 2 . Asymptotic value of the orbital elements ............5 3 . Asymptotic perturbations of the hyperbolic excess velocity vector ....6 4 . Calculated and observed perturbations at 240 min .........6 5 . Position and velocity at T = 240 min ..............7 FIGURES 1. Perturbation of the semi-transverse axis ............5 2 . Perturbation of the eccentricity ................5 3 . Perturbation of the time of peri-focal passage .......... 6 4 . Perturbation of the inclination ................6 5 . Perturbation of the argument of perigee .............6 6 . Perturbation of the longitudeof the ascending node ........7 7. True anomaly vs time (e = 1.25) ...............7 JPL TECHNICAL REPORT NO. 32-131 The perturbations of the hyperbolic orbital elements of a vehicle in the gravitational field of an oblate planet are derived as functions of the initial osculating elements. Assumptions are made that atmos- pheric drag is absent and that the gravitational potential of the planet may be represented adequately by the principal term and the second harmonic. An example of an Earth-escape mission is presented in which a comparison is made between calculated orbital perturbations and results from a numerical integration of the equations of motion. 1. INTRODUCTION Although the variations of the orbital elements of Earth perturbations of an interplanetary probe by an oblate satellites have been derived numerous times, nevertheless planet at encounter. a treatment of the variation of the orbital elements of a vehicle for an escape mission is lacking except for the 'Since 1960, when this report was originally prepared, there have case of perturbed hyperbolic motion in the equatorial been two additional papers published concerning the perturbing plane (see Ref. 1 ) There are several instances in which effects of an oblate Planet upon a vehicle on a hYFrbolic tra- jectory (Ref. 2 and 3). Of particular interest is Ref. 3 by G. Hori the Of this paper are Of interest, the first being in which the Von-&ipel was used to derive the pefiw&d the perturbations Of the escape from Earth Of motion. The present paper uses the method of Poisson-LaGrange, a lunar or interplanetary probe, and the second being the as does Ref. 2. 1 JPL TECHNICAL REPORT NO. 32-131 II. PERTURBATIVE POTENTIAL If an axially symmetric density distribution of the planet where is assumed, the gravitational potential function U may be 1 P, (sin+) = -(3sin, 4 - 1 expanded in terms of surface zonal harmonics of the form 2 (3) (Ref. 4 and 5). The distance from the gravicenter T and gravicenter lati- tude + expressed as functions of the osculating elements of the vehicle are where T = p/(l + ecosv) (4) J, = the coefficients of the zonal hormonics P, ( sin 4 ) and k2 = GM = gravitational constant of attracting body r = radial distance from gravicenter in equatorial radii sin4 = sinisin (W + V) (5) 0 = gravicenter latitude where In an analysis of the perturbations of the bound orbit p = h2/k2= ( -u) (e' - 1) = semilatus rectum of an Earth satellite, it is necessary to consideral several of a = semitransverse axis the harmonic terms to describe accurately the motion of h = angular momentum the vehicle over an extended period of time (Ref. 6 and e = eccentricity i = inclination to the equatorial plane 7). However, only the second harmonic J2 contributes significantly to the perturbations of the elements of a non- w = argument of perigee bound or hyperbolic orbit, since the time spent in the v = trueanomaly vicinity of the planet is short and the effects of the higher The disturbing function R, when expressed in terms of order harmonics are negligible. the osculating elements of the orbit, becomes The perturbing potential R due to the second harmonic R = ],(k2/2pS)(1 + BCOS v)3 [l - -sin23. i J2 is given by 2 + 3 ;z sin2 i cos (2, + 2~)](6) 111. VARIATION OF ORBITAL ELEMENTS The solution of the differential equations representing E = - k2/2a = undisturbed energy the variations of the orbital elements with respect to time is derived by a method analogous to that employed by h, = h cos i = component of angular momentum along Kozai (Ref. 8) in investigations of elliptic motion. The polar axis orbital elements E, h, h,, T, O, and sz are employed, since they both form a canonical set and are also well behaved T = time of perifocal passage for eccentricities near unity. The elements not previously defined are sz = longitude of the ascending node 2 ~ ~~ JPL TECHNICAL REPORT NO. 32-131 The rate of change of the orbital elements with time n( t - T), where the mean motion n is a function of the may be expressed in terms of the differential coefficients element E through the harmonic law n = ( 2E):’’2k-2.The of R as variable A.1 corresponds to the mean anomaly for hyper- bolic motion and is a function of both the true anomaly and the eccentricity, and hence the disturbing function involves the element h both explicitly and also implicitly through the true anomaly. dddt = -?R/ah (9) Since an axially symmetric density distribution is as- sumed, the disturbing function R is independent of 0,and Eq. ( 12) becomes dhz/dt = 0 (14) and h, = h,(O) = const (15) Since the orbital elements of the vehicle change only the component of angular momentum along the polar slightly due to the effects of the oblateness, the pertur- axis being unaffected by the oblateness. bations of the first order may be derived by considering the elements appearing in R on the right side of Eq. (7) The perturbations of the remaining five orbital ele- through (12) constant. If the elements are assumed con- ments are given in Eq. ( 16) through (20) where M and stant, the true anomaly may be regarded as a known appearing in Eq. (20) are the mean anomaly and mean function of time, and the independent variable may be motion, respectively. These expressions for the perturba- transformed from time to the true anomaly by tions are to be evaluated between the limits v0 and V. SE = R1’ JVO h sin2 i Sh = J, - 3e cos + ”) + 3 cos (2w + 2”) + e cos (20 + 3”) (17) 4p? I:, cos i 60 = -I., -i6(” + e sin “) - [3e sin (2, + v) + 3 sin (2” + 2”) + e sin (2~+ 3v)] - 4pX 1 >:, e 6” + cos is0 = J, 3e3 sin (2w - (12 - 21e’) sin (2w - 36e sin (2, 2v) - 48e v) + + v) + (28 + lle?) sin (2, + 3“) - 18e sin (2, + 4”) - 3e2 sin (20~+ 5~) The perturbations of the semitransverse axis a, eccen- d/dt = h/r2 d/dv (13) tricity e, and orbital inclination i are found from The disturbing function involves the element E and ‘T implicity through the true anomaly from the relation M = 3 JPL TECHNICAL REPORT NO. 32-131 ah and 4h being the right ascension and declination of the hyperbolic excess velocity vector and Vh being the angle between the hyperbolic excess velocity vector and peri- 8i = cot i (6h/h) (23) apsis, The expressions for the variations of the elements E, The angle Vh is related to the eccentricity of the hyper- h, h,, and 0 are well behaved for eccentricities near unity; bola by however, the expression for the variation of the time of perifocal passage breaks down when e = 1. Fortunately, cos Vh = - l/e (28) in most cases of interest, the change in the time of peri- focal passage is of not great importance. The perturbation of the magnitude of the hyperbolic excess velocity vector is given by Of additional importance in the analysis of interplan- etary escape trajectories is the velocity vector of the 8vh = 8E/Vh (29) vehicle at a great distance from Earth. The magnitude of this hyperbolic excess velocity vector is given by the whereas the perturbations of the right ascension and dec- following function of the energy: lination are given by The unit vector S in the direction of the hyperbolic excess velocity vector and also in the direction of the outgoing asymptote of the hyperbola has direction cosines of the form cos i sin (w + VI,) sin i cos (0 + Vh) 84h = (80 f 8Vh) cos $‘h si + cos 4h S, = cos C#3h sin ah = sin 0 cos (W + vIl) + The perturbation of the angle between S and periapsis cos 0 sin (W + vII)cos i (26) is given by S, = sin 41,= sin (W + Vh) sin i (27) IV. EXAMPLE OF AN EARTH-ESCAPE MISSION As an example of the application of the deriviation of 6 as functions of the true anomaly. The time along the the perturbations, an Earth-escape mission with the initial trajectory may be determined from Fig. 7 as a function of parameters shown in Table 1 was used. The perturbations the true anomaly. The maximum true anomaly is approxi- of the elements a, e, i, T, W, and R were calculated from mately 2.5 rad; in the following discussion the time of Eq.
Load more

Recommended publications
  • AAS 13-250 Hohmann Spiral Transfer with Inclination Change Performed
    AAS 13-250 Hohmann Spiral Transfer With Inclination Change Performed By Low-Thrust System Steven Owens1 and Malcolm Macdonald2 This paper investigates the Hohmann Spiral Transfer (HST), an orbit transfer method previously developed by the authors incorporating both high and low- thrust propulsion systems, using the low-thrust system to perform an inclination change as well as orbit transfer. The HST is similar to the bi-elliptic transfer as the high-thrust system is first used to propel the spacecraft beyond the target where it is used again to circularize at an intermediate orbit. The low-thrust system is then activated and, while maintaining this orbit altitude, used to change the orbit inclination to suit the mission specification. The low-thrust system is then used again to reduce the spacecraft altitude by spiraling in-toward the target orbit. An analytical analysis of the HST utilizing the low-thrust system for the inclination change is performed which allows a critical specific impulse ratio to be derived determining the point at which the HST consumes the same amount of fuel as the Hohmann transfer. A critical ratio is found for both a circular and elliptical initial orbit. These equations are validated by a numerical approach before being compared to the HST utilizing the high-thrust system to perform the inclination change. An additional critical ratio comparing the HST utilizing the low-thrust system for the inclination change with its high-thrust counterpart is derived and by using these three critical ratios together, it can be determined when each transfer offers the lowest fuel mass consumption.
    [Show full text]
  • Appendix a Orbits
    Appendix A Orbits As discussed in the Introduction, a good ¯rst approximation for satellite motion is obtained by assuming the spacecraft is a point mass or spherical body moving in the gravitational ¯eld of a spherical planet. This leads to the classical two-body problem. Since we use the term body to refer to a spacecraft of ¯nite size (as in rigid body), it may be more appropriate to call this the two-particle problem, but I will use the term two-body problem in its classical sense. The basic elements of orbital dynamics are captured in Kepler's three laws which he published in the 17th century. His laws were for the orbital motion of the planets about the Sun, but are also applicable to the motion of satellites about planets. The three laws are: 1. The orbit of each planet is an ellipse with the Sun at one focus. 2. The line joining the planet to the Sun sweeps out equal areas in equal times. 3. The square of the period of a planet is proportional to the cube of its mean distance to the sun. The ¯rst law applies to most spacecraft, but it is also possible for spacecraft to travel in parabolic and hyperbolic orbits, in which case the period is in¯nite and the 3rd law does not apply. However, the 2nd law applies to all two-body motion. Newton's 2nd law and his law of universal gravitation provide the tools for generalizing Kepler's laws to non-elliptical orbits, as well as for proving Kepler's laws.
    [Show full text]
  • Astrodynamics
    Politecnico di Torino SEEDS SpacE Exploration and Development Systems Astrodynamics II Edition 2006 - 07 - Ver. 2.0.1 Author: Guido Colasurdo Dipartimento di Energetica Teacher: Giulio Avanzini Dipartimento di Ingegneria Aeronautica e Spaziale e-mail: [email protected] Contents 1 Two–Body Orbital Mechanics 1 1.1 BirthofAstrodynamics: Kepler’sLaws. ......... 1 1.2 Newton’sLawsofMotion ............................ ... 2 1.3 Newton’s Law of Universal Gravitation . ......... 3 1.4 The n–BodyProblem ................................. 4 1.5 Equation of Motion in the Two-Body Problem . ....... 5 1.6 PotentialEnergy ................................. ... 6 1.7 ConstantsoftheMotion . .. .. .. .. .. .. .. .. .... 7 1.8 TrajectoryEquation .............................. .... 8 1.9 ConicSections ................................... 8 1.10 Relating Energy and Semi-major Axis . ........ 9 2 Two-Dimensional Analysis of Motion 11 2.1 ReferenceFrames................................. 11 2.2 Velocity and acceleration components . ......... 12 2.3 First-Order Scalar Equations of Motion . ......... 12 2.4 PerifocalReferenceFrame . ...... 13 2.5 FlightPathAngle ................................. 14 2.6 EllipticalOrbits................................ ..... 15 2.6.1 Geometry of an Elliptical Orbit . ..... 15 2.6.2 Period of an Elliptical Orbit . ..... 16 2.7 Time–of–Flight on the Elliptical Orbit . .......... 16 2.8 Extensiontohyperbolaandparabola. ........ 18 2.9 Circular and Escape Velocity, Hyperbolic Excess Speed . .............. 18 2.10 CosmicVelocities
    [Show full text]
  • Optimisation of Propellant Consumption for Power Limited Rockets
    Delft University of Technology Faculty Electrical Engineering, Mathematics and Computer Science Delft Institute of Applied Mathematics Optimisation of Propellant Consumption for Power Limited Rockets. What Role do Power Limited Rockets have in Future Spaceflight Missions? (Dutch title: Optimaliseren van brandstofverbruik voor vermogen gelimiteerde raketten. De rol van deze raketten in toekomstige ruimtevlucht missies. ) A thesis submitted to the Delft Institute of Applied Mathematics as part to obtain the degree of BACHELOR OF SCIENCE in APPLIED MATHEMATICS by NATHALIE OUDHOF Delft, the Netherlands December 2017 Copyright c 2017 by Nathalie Oudhof. All rights reserved. BSc thesis APPLIED MATHEMATICS \ Optimisation of Propellant Consumption for Power Limited Rockets What Role do Power Limite Rockets have in Future Spaceflight Missions?" (Dutch title: \Optimaliseren van brandstofverbruik voor vermogen gelimiteerde raketten De rol van deze raketten in toekomstige ruimtevlucht missies.)" NATHALIE OUDHOF Delft University of Technology Supervisor Dr. P.M. Visser Other members of the committee Dr.ir. W.G.M. Groenevelt Drs. E.M. van Elderen 21 December, 2017 Delft Abstract In this thesis we look at the most cost-effective trajectory for power limited rockets, i.e. the trajectory which costs the least amount of propellant. First some background information as well as the differences between thrust limited and power limited rockets will be discussed. Then the optimal trajectory for thrust limited rockets, the Hohmann Transfer Orbit, will be explained. Using Optimal Control Theory, the optimal trajectory for power limited rockets can be found. Three trajectories will be discussed: Low Earth Orbit to Geostationary Earth Orbit, Earth to Mars and Earth to Saturn. After this we compare the propellant use of the thrust limited rockets for these trajectories with the power limited rockets.
    [Show full text]
  • Flight and Orbital Mechanics
    Flight and Orbital Mechanics Lecture slides Challenge the future 1 Flight and Orbital Mechanics AE2-104, lecture hours 21-24: Interplanetary flight Ron Noomen October 25, 2012 AE2104 Flight and Orbital Mechanics 1 | Example: Galileo VEEGA trajectory Questions: • what is the purpose of this mission? • what propulsion technique(s) are used? • why this Venus- Earth-Earth sequence? • …. [NASA, 2010] AE2104 Flight and Orbital Mechanics 2 | Overview • Solar System • Hohmann transfer orbits • Synodic period • Launch, arrival dates • Fast transfer orbits • Round trip travel times • Gravity Assists AE2104 Flight and Orbital Mechanics 3 | Learning goals The student should be able to: • describe and explain the concept of an interplanetary transfer, including that of patched conics; • compute the main parameters of a Hohmann transfer between arbitrary planets (including the required ΔV); • compute the main parameters of a fast transfer between arbitrary planets (including the required ΔV); • derive the equation for the synodic period of an arbitrary pair of planets, and compute its numerical value; • derive the equations for launch and arrival epochs, for a Hohmann transfer between arbitrary planets; • derive the equations for the length of the main mission phases of a round trip mission, using Hohmann transfers; and • describe the mechanics of a Gravity Assist, and compute the changes in velocity and energy. Lecture material: • these slides (incl. footnotes) AE2104 Flight and Orbital Mechanics 4 | Introduction The Solar System (not to scale): [Aerospace
    [Show full text]
  • A Strategy for the Eigenvector Perturbations of the Reynolds Stress Tensor in the Context of Uncertainty Quantification
    Center for Turbulence Research 425 Proceedings of the Summer Program 2016 A strategy for the eigenvector perturbations of the Reynolds stress tensor in the context of uncertainty quantification By R. L. Thompsony, L. E. B. Sampaio, W. Edeling, A. A. Mishra AND G. Iaccarino In spite of increasing progress in high fidelity turbulent simulation avenues like Di- rect Numerical Simulation (DNS) and Large Eddy Simulation (LES), Reynolds-Averaged Navier-Stokes (RANS) models remain the predominant numerical recourse to study com- plex engineering turbulent flows. In this scenario, it is imperative to provide reliable estimates for the uncertainty in RANS predictions. In the recent past, an uncertainty estimation framework relying on perturbations to the modeled Reynolds stress tensor has been widely applied with satisfactory results. Many of these investigations focus on perturbing only the Reynolds stress eigenvalues in the Barycentric map, thus ensuring realizability. However, these restrictions apply to the eigenvalues of that tensor only, leav- ing the eigenvectors free of any limiting condition. In the present work, we propose the use of the Reynolds stress transport equation as a constraint for the eigenvector pertur- bations of the Reynolds stress anisotropy tensor once the eigenvalues of this tensor are perturbed. We apply this methodology to a convex channel and show that the ensuing eigenvector perturbations are a more accurate measure of uncertainty when compared with a pure eigenvalue perturbation. 1. Introduction Turbulent flows are present in a wide number of engineering design problems. Due to the disparate character of such flows, predictive methods must be robust, so as to be easily applicable for most of these cases, yet possessing a high degree of accuracy in each.
    [Show full text]
  • Spacecraft Guidance Techniques for Maximizing Mission Success
    Utah State University DigitalCommons@USU All Graduate Theses and Dissertations Graduate Studies 5-2014 Spacecraft Guidance Techniques for Maximizing Mission Success Shane B. Robinson Utah State University Follow this and additional works at: https://digitalcommons.usu.edu/etd Part of the Mechanical Engineering Commons Recommended Citation Robinson, Shane B., "Spacecraft Guidance Techniques for Maximizing Mission Success" (2014). All Graduate Theses and Dissertations. 2175. https://digitalcommons.usu.edu/etd/2175 This Dissertation is brought to you for free and open access by the Graduate Studies at DigitalCommons@USU. It has been accepted for inclusion in All Graduate Theses and Dissertations by an authorized administrator of DigitalCommons@USU. For more information, please contact [email protected]. SPACECRAFT GUIDANCE TECHNIQUES FOR MAXIMIZING MISSION SUCCESS by Shane B. Robinson A dissertation submitted in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY in Mechanical Engineering Approved: Dr. David K. Geller Dr. Jacob H. Gunther Major Professor Committee Member Dr. Warren F. Phillips Dr. Charles M. Swenson Committee Member Committee Member Dr. Stephen A. Whitmore Dr. Mark R. McLellan Committee Member Vice President for Research and Dean of the School of Graduate Studies UTAH STATE UNIVERSITY Logan, Utah 2013 [This page intentionally left blank] iii Copyright c Shane B. Robinson 2013 All Rights Reserved [This page intentionally left blank] v Abstract Spacecraft Guidance Techniques for Maximizing Mission Success by Shane B. Robinson, Doctor of Philosophy Utah State University, 2013 Major Professor: Dr. David K. Geller Department: Mechanical and Aerospace Engineering Traditional spacecraft guidance techniques have the objective of deterministically min- imizing fuel consumption.
    [Show full text]
  • Elliptical Orbits
    1 Ellipse-geometry 1.1 Parameterization • Functional characterization:(a: semi major axis, b ≤ a: semi minor axis) x2 y 2 b p + = 1 ⇐⇒ y(x) = · ± a2 − x2 (1) a b a • Parameterization in cartesian coordinates, which follows directly from Eq. (1): x a · cos t = with 0 ≤ t < 2π (2) y b · sin t – The origin (0, 0) is the center of the ellipse and the auxilliary circle with radius a. √ – The focal points are located at (±a · e, 0) with the eccentricity e = a2 − b2/a. • Parameterization in polar coordinates:(p: parameter, 0 ≤ < 1: eccentricity) p r(ϕ) = (3) 1 + e cos ϕ – The origin (0, 0) is the right focal point of the ellipse. – The major axis is given by 2a = r(0) − r(π), thus a = p/(1 − e2), the center is therefore at − pe/(1 − e2), 0. – ϕ = 0 corresponds to the periapsis (the point closest to the focal point; which is also called perigee/perihelion/periastron in case of an orbit around the Earth/sun/star). The relation between t and ϕ of the parameterizations in Eqs. (2) and (3) is the following: t r1 − e ϕ tan = · tan (4) 2 1 + e 2 1.2 Area of an elliptic sector As an ellipse is a circle with radius a scaled by a factor b/a in y-direction (Eq. 1), the area of an elliptic sector PFS (Fig. ??) is just this fraction of the area PFQ in the auxiliary circle. b t 2 1 APFS = · · πa − · ae · a sin t a 2π 2 (5) 1 = (t − e sin t) · a b 2 The area of the full ellipse (t = 2π) is then, of course, Aellipse = π a b (6) Figure 1: Ellipse and auxilliary circle.
    [Show full text]
  • Richard Dalbello Vice President, Legal and Government Affairs Intelsat General
    Commercial Management of the Space Environment Richard DalBello Vice President, Legal and Government Affairs Intelsat General The commercial satellite industry has billions of dollars of assets in space and relies on this unique environment for the development and growth of its business. As a result, safety and the sustainment of the space environment are two of the satellite industry‘s highest priorities. In this paper I would like to provide a quick survey of past and current industry space traffic control practices and to discuss a few key initiatives that the industry is developing in this area. Background The commercial satellite industry has been providing essential space services for almost as long as humans have been exploring space. Over the decades, this industry has played an active role in developing technology, worked collaboratively to set standards, and partnered with government to develop successful international regulatory regimes. Success in both commercial and government space programs has meant that new demands are being placed on the space environment. This has resulted in orbital crowding, an increase in space debris, and greater demand for limited frequency resources. The successful management of these issues will require a strong partnership between government and industry and the careful, experience-based expansion of international law and diplomacy. Throughout the years, the satellite industry has never taken for granted the remarkable environment in which it works. Industry has invested heavily in technology and sought out the best and brightest minds to allow the full, but sustainable exploitation of the space environment. Where problems have arisen, such as space debris or electronic interference, industry has taken 1 the initiative to deploy new technologies and adopt new practices to minimize negative consequences.
    [Show full text]
  • Perturbation Theory in Celestial Mechanics
    Perturbation Theory in Celestial Mechanics Alessandra Celletti Dipartimento di Matematica Universit`adi Roma Tor Vergata Via della Ricerca Scientifica 1, I-00133 Roma (Italy) ([email protected]) December 8, 2007 Contents 1 Glossary 2 2 Definition 2 3 Introduction 2 4 Classical perturbation theory 4 4.1 The classical theory . 4 4.2 The precession of the perihelion of Mercury . 6 4.2.1 Delaunay action–angle variables . 6 4.2.2 The restricted, planar, circular, three–body problem . 7 4.2.3 Expansion of the perturbing function . 7 4.2.4 Computation of the precession of the perihelion . 8 5 Resonant perturbation theory 9 5.1 The resonant theory . 9 5.2 Three–body resonance . 10 5.3 Degenerate perturbation theory . 11 5.4 The precession of the equinoxes . 12 6 Invariant tori 14 6.1 Invariant KAM surfaces . 14 6.2 Rotational tori for the spin–orbit problem . 15 6.3 Librational tori for the spin–orbit problem . 16 6.4 Rotational tori for the restricted three–body problem . 17 6.5 Planetary problem . 18 7 Periodic orbits 18 7.1 Construction of periodic orbits . 18 7.2 The libration in longitude of the Moon . 20 1 8 Future directions 20 9 Bibliography 21 9.1 Books and Reviews . 21 9.2 Primary Literature . 22 1 Glossary KAM theory: it provides the persistence of quasi–periodic motions under a small perturbation of an integrable system. KAM theory can be applied under quite general assumptions, i.e. a non– degeneracy of the integrable system and a diophantine condition of the frequency of motion.
    [Show full text]
  • New Closed-Form Solutions for Optimal Impulsive Control of Spacecraft Relative Motion
    New Closed-Form Solutions for Optimal Impulsive Control of Spacecraft Relative Motion Michelle Chernick∗ and Simone D'Amicoy Aeronautics and Astronautics, Stanford University, Stanford, California, 94305, USA This paper addresses the fuel-optimal guidance and control of the relative motion for formation-flying and rendezvous using impulsive maneuvers. To meet the requirements of future multi-satellite missions, closed-form solutions of the inverse relative dynamics are sought in arbitrary orbits. Time constraints dictated by mission operations and relevant perturbations acting on the formation are taken into account by splitting the optimal recon- figuration in a guidance (long-term) and control (short-term) layer. Both problems are cast in relative orbit element space which allows the simple inclusion of secular and long-periodic perturbations through a state transition matrix and the translation of the fuel-optimal optimization into a minimum-length path-planning problem. Due to the proper choice of state variables, both guidance and control problems can be solved (semi-)analytically leading to optimal, predictable maneuvering schemes for simple on-board implementation. Besides generalizing previous work, this paper finds four new in-plane and out-of-plane (semi-)analytical solutions to the optimal control problem in the cases of unperturbed ec- centric and perturbed near-circular orbits. A general delta-v lower bound is formulated which provides insight into the optimality of the control solutions, and a strong analogy between elliptic Hohmann transfers and formation-flying control is established. Finally, the functionality, performance, and benefits of the new impulsive maneuvering schemes are rigorously assessed through numerical integration of the equations of motion and a systematic comparison with primer vector optimal control.
    [Show full text]
  • Coordinate Systems in Geodesy
    COORDINATE SYSTEMS IN GEODESY E. J. KRAKIWSKY D. E. WELLS May 1971 TECHNICALLECTURE NOTES REPORT NO.NO. 21716 COORDINATE SYSTElVIS IN GEODESY E.J. Krakiwsky D.E. \Vells Department of Geodesy and Geomatics Engineering University of New Brunswick P.O. Box 4400 Fredericton, N .B. Canada E3B 5A3 May 1971 Latest Reprinting January 1998 PREFACE In order to make our extensive series of lecture notes more readily available, we have scanned the old master copies and produced electronic versions in Portable Document Format. The quality of the images varies depending on the quality of the originals. The images have not been converted to searchable text. TABLE OF CONTENTS page LIST OF ILLUSTRATIONS iv LIST OF TABLES . vi l. INTRODUCTION l 1.1 Poles~ Planes and -~es 4 1.2 Universal and Sidereal Time 6 1.3 Coordinate Systems in Geodesy . 7 2. TERRESTRIAL COORDINATE SYSTEMS 9 2.1 Terrestrial Geocentric Systems • . 9 2.1.1 Polar Motion and Irregular Rotation of the Earth • . • • . • • • • . 10 2.1.2 Average and Instantaneous Terrestrial Systems • 12 2.1. 3 Geodetic Systems • • • • • • • • • • . 1 17 2.2 Relationship between Cartesian and Curvilinear Coordinates • • • • • • • . • • 19 2.2.1 Cartesian and Curvilinear Coordinates of a Point on the Reference Ellipsoid • • • • • 19 2.2.2 The Position Vector in Terms of the Geodetic Latitude • • • • • • • • • • • • • • • • • • • 22 2.2.3 Th~ Position Vector in Terms of the Geocentric and Reduced Latitudes . • • • • • • • • • • • 27 2.2.4 Relationships between Geodetic, Geocentric and Reduced Latitudes • . • • • • • • • • • • 28 2.2.5 The Position Vector of a Point Above the Reference Ellipsoid . • • . • • • • • • . .• 28 2.2.6 Transformation from Average Terrestrial Cartesian to Geodetic Coordinates • 31 2.3 Geodetic Datums 33 2.3.1 Datum Position Parameters .
    [Show full text]