DOCSLIB.ORG

Satellite Orbit and Ephemeris Determination Using Inter Satellite Links

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Satellite Orbit and Ephemeris Determination Using Inter Satellite Links Satellite Orbit and Ephemeris Determination using Inter Satellite Links by Robert Wolf Vollständiger Abdruck der an der Fakultät für Bauingenieur- und Vermessungswesen der Universität der Bundeswehr München zur Erlangung des akademischen Grades eines Doktors der Ingenieurswissenschaften (Dr.-Ing.) eingereichten Dissertation. (© 2000) Inter Satellite Links Page ii R. Wolf Inter Satellite Links Abstract Global navigation satellite systems like GPS, GLONASS or the future systems like Galileo require precise orbit and clock estimates in order to provide high positioning performance. Within the frame of this Ph. D. thesis, the theory of orbit determination and orbit computation is reviewed and a new approach for precise orbit and ephemeris determination using inter- satellite links is developed. To investigate the achievable accuracy, models of the various perturbing forces acting on a satellite have been elaborated and coded in a complex software package, allowing system level performance analysis as well as detailed evaluation of orbit prediction and orbit estimation algorithms. Several satellite constellations have been simulated, involving nearly all classes of orbit altitude and the results are compared. The purpose of orbit determination in a satellite navigation system is the derivation of ephemeris parameters which can be broadcast to the user community (or the other satellites) and allow easy computation of the satellites position at the desired epoch. The broadcast ephemeris model of both today existing satellite navigation systems, GPS and GLONASS are investigated, as well as two new models developed within this thesis, which are derivates of the GLONASS model. Furthermore, the topic of autonomous onboard processing is addressed. A conceptual design for an onboard orbit estimator is proposed and investigated with respect to the computational load. The algorithms have been implemented. The main benefits of ISL onboard processing, especially with respect to the great potential to ephemeris and clock state monitoring are investigated using complex simulations of failure scenarios. By simulating several types of non-integrity cases, it is showed that one single fault detection mechanism is likely to be insufficient. Within the algorithm design of the onboard processor, a reasonable combination of fault detection mechanisms is presented, covering different fault cases. Zusammenfassung Globale Navigationssysteme wie GPS, GLONASS oder zukünftige Systeme wie Galileo erfordern die hochpräzise Bestimmung der Orbital- und Uhrenparameter, um hohe Navigationsgenauigkeit bieten zu können. Im Rahmen dieser Dissertation wurde die Theorie der Orbitprädiktion und der Orbitbestimmung erörtert und ein neuer Ansatz für die präzisen Orbitbestimmung mit Hilfe von Intersatelliten-Messungen entwickelt. Um die erreichbare Genauigkeit und Präzision der Orbitbestimmung zu untersuchen, wurden mathematische Modelle der zahlreiche Orbitstörungen erarbeitet und in einem komplexen Software-Paket implemetiert. Dieses bietet die Möglichkeit für Systemstudien von Satellitennavigations- Systemen beliebiger Orbitklassen, sowie zur detaillierten Untersuchung spezieller Fragestellungen der Orbitprädiktion und -bestimmung. Eine Reihe von Simulationen mit existierenden sowie fiktiven Satelliten-Navigations-Systemen wurden durchgeführt, deren Ergebnisse in dieser Arbeit präsentiert werden. Die präzise Orbitbestimmung in einem SatNav-System ist kein Selbstzweck, sondern dient lediglich der Bestimmung der Ephemeridenparameter, die - vom Satellite gesendet - es dem Nutzer-Empfänger erlauben, mit Hilfe einfacher Berechnungen die Position des Satelliten zu ermitteln. Die Ephemeridenformate beider existierender SatNav-Systeme - GPS und R. Wolf Page iii Inter Satellite Links GLONASS - wurden untersucht und mit zwei weiteren Formaten verglichen, die im Rahmen dieser Arbeit entwickelt wurden. Desweiteren wurde das Thema der bordautonomen Verarbeitung von Messungen behandelt. Ein konzeptuelles Design für einen Onboard-Prozessor wurde vorgeschlagen und die Algorithmen implementiert. Dabei erfolgte eine Abschätzung der benötigten Prozessorleistung. Einer der Hauptvorteile der bordautonomen Verarbeitung von Intersatellitenmessungen, die Möglichkeit zur Überwachung der Integrität der Ephemeriden und Uhrenparameter, wurde in komplexen Simulationen untersucht. Durch die Simulation verschiedener Fehlerfälle wurde gezeigt, das kein Detektionsmechanismus allein, wohl aber eine sinnvolle Kombination solcher Mechanismen, zur bordautonomen Integritätsüberwachung geeignet sind. Die Ergebissen werden hier präsentiert. Page iv R. Wolf Inter Satellite Links Table of Contents 1 INTRODUCTION............................................................................................................................ 1 2 ISL OBSERVATION MODEL....................................................................................................... 4 2.1 DERIVATION OF THE RANGE EQUATION .................................................................................. 4 2.2 DERIVATION OF THE RANGE RATE EQUATION ........................................................................ 5 3 STATE ESTIMATION.................................................................................................................... 7 3.1 LINEARIZATION OF DYNAMIC AND OBSERVATION MODEL .................................................... 7 3.2 STATE VECTOR........................................................................................................................... 8 3.3 STATE TRANSITION AND TRANSITION MATRIX ..................................................................... 10 3.4 LEAST SQUARES BATCH ESTIMATION .................................................................................... 13 3.4.1 WEIGHTED LEAST SQUARES ............................................................................................... 14 3.4.2 INTRODUCING APRIORI STATISTIC INFORMATION............................................................... 15 3.5 KALMAN FILTERING ................................................................................................................ 15 3.5.1 REAL TIME ESTIMATION ..................................................................................................... 15 3.5.2 FILTERING TO EPOCH........................................................................................................... 16 3.5.3 FILTER STRUCTURES ........................................................................................................... 17 4 ORBIT COMPUTATION.............................................................................................................19 4.1 ANALYTICAL SOLUTION .......................................................................................................... 19 4.1.1 KEPLER ORBITS ................................................................................................................... 20 4.1.2 ACCOUNTING FOR SECULAR PERTURBATIONS.................................................................... 22 4.2 NUMERICAL INTEGRATION OF THE EQUATIONS OF MOTION............................................... 23 4.2.1 EARTH’S GRAVITY .............................................................................................................. 26 4.2.1.1 Computation of Legendre Polynomials and Functions................................................. 27 4.2.1.2 Normalisation................................................................................................................ 28 4.2.1.3 Computation of Gravity ................................................................................................ 29 4.2.2 THIRD BODY ATTRACTION.................................................................................................. 32 4.2.3 SOLAR PRESSURE ................................................................................................................ 32 4.2.4 AIR DRAG ............................................................................................................................ 33 4.2.5 SOLID EARTH TIDES ............................................................................................................ 35 4.2.6 OCEAN TIDES....................................................................................................................... 36 4.2.7 EARTH ALBEDO ................................................................................................................... 37 4.2.8 VEHICLE THRUST ................................................................................................................ 37 4.3 FORCE MODEL ERRORS........................................................................................................... 39 4.3.1 EARTH'S GRAVITY ............................................................................................................... 39 4.3.2 THIRD BODY ATTRACTION (DIRECT TIDAL EFFECTS) ........................................................ 47 4.3.3 SOLAR RADIATION PRESSURE ............................................................................................. 50 4.3.4 AIR DRAG ............................................................................................................................ 52 4.3.5 OTHER PERTURBATIONS...................................................................................................... 55 4.3.6 NUMERICAL ERRORS ..........................................................................................................
Load more

Recommended publications
  • Analysis of Perturbations and Station-Keeping Requirements in Highly-Inclined Geosynchronous Orbits
    ANALYSIS OF PERTURBATIONS AND STATION-KEEPING REQUIREMENTS IN HIGHLY-INCLINED GEOSYNCHRONOUS ORBITS Elena Fantino(1), Roberto Flores(2), Alessio Di Salvo(3), and Marilena Di Carlo(4) (1)Space Studies Institute of Catalonia (IEEC), Polytechnic University of Catalonia (UPC), E.T.S.E.I.A.T., Colom 11, 08222 Terrassa (Spain), [email protected] (2)International Center for Numerical Methods in Engineering (CIMNE), Polytechnic University of Catalonia (UPC), Building C1, Campus Norte, UPC, Gran Capitan,´ s/n, 08034 Barcelona (Spain) (3)NEXT Ingegneria dei Sistemi S.p.A., Space Innovation System Unit, Via A. Noale 345/b, 00155 Roma (Italy), [email protected] (4)Department of Mechanical and Aerospace Engineering, University of Strathclyde, 75 Montrose Street, Glasgow G1 1XJ (United Kingdom), [email protected] Abstract: There is a demand for communications services at high latitudes that is not well served by conventional geostationary satellites. Alternatives using low-altitude orbits require too large constellations. Other options are the Molniya and Tundra families (critically-inclined, eccentric orbits with the apogee at high latitudes). In this work we have considered derivatives of the Tundra type with different inclinations and eccentricities. By means of a high-precision model of the terrestrial gravity field and the most relevant environmental perturbations, we have studied the evolution of these orbits during a period of two years. The effects of the different perturbations on the constellation ground track (which is more important for coverage than the orbital elements themselves) have been identified. We show that, in order to maintain the ground track unchanged, the most important parameters are the orbital period and the argument of the perigee.
    [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]
  • A Hot Subdwarf-White Dwarf Super-Chandrasekhar Candidate
    A hot subdwarf–white dwarf super-Chandrasekhar candidate supernova Ia progenitor Ingrid Pelisoli1,2*, P. Neunteufel3, S. Geier1, T. Kupfer4,5, U. Heber6, A. Irrgang6, D. Schneider6, A. Bastian1, J. van Roestel7, V. Schaffenroth1, and B. N. Barlow8 1Institut fur¨ Physik und Astronomie, Universitat¨ Potsdam, Haus 28, Karl-Liebknecht-Str. 24/25, D-14476 Potsdam-Golm, Germany 2Department of Physics, University of Warwick, Coventry, CV4 7AL, UK 3Max Planck Institut fur¨ Astrophysik, Karl-Schwarzschild-Straße 1, 85748 Garching bei Munchen¨ 4Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA 5Texas Tech University, Department of Physics & Astronomy, Box 41051, 79409, Lubbock, TX, USA 6Dr. Karl Remeis-Observatory & ECAP, Astronomical Institute, Friedrich-Alexander University Erlangen-Nuremberg (FAU), Sternwartstr. 7, 96049 Bamberg, Germany 7Division of Physics, Mathematics and Astronomy, California Institute of Technology, Pasadena, CA 91125, USA 8Department of Physics and Astronomy, High Point University, High Point, NC 27268, USA *[email protected] ABSTRACT Supernova Ia are bright explosive events that can be used to estimate cosmological distances, allowing us to study the expansion of the Universe. They are understood to result from a thermonuclear detonation in a white dwarf that formed from the exhausted core of a star more massive than the Sun. However, the possible progenitor channels leading to an explosion are a long-standing debate, limiting the precision and accuracy of supernova Ia as distance indicators. Here we present HD 265435, a binary system with an orbital period of less than a hundred minutes, consisting of a white dwarf and a hot subdwarf — a stripped core-helium burning star.
    [Show full text]
  • Electric Propulsion System Scaling for Asteroid Capture-And-Return Missions
    Electric propulsion system scaling for asteroid capture-and-return missions Justin M. Little⇤ and Edgar Y. Choueiri† Electric Propulsion and Plasma Dynamics Laboratory, Princeton University, Princeton, NJ, 08544 The requirements for an electric propulsion system needed to maximize the return mass of asteroid capture-and-return (ACR) missions are investigated in detail. An analytical model is presented for the mission time and mass balance of an ACR mission based on the propellant requirements of each mission phase. Edelbaum’s approximation is used for the Earth-escape phase. The asteroid rendezvous and return phases of the mission are modeled as a low-thrust optimal control problem with a lunar assist. The numerical solution to this problem is used to derive scaling laws for the propellant requirements based on the maneuver time, asteroid orbit, and propulsion system parameters. Constraining the rendezvous and return phases by the synodic period of the target asteroid, a semi- empirical equation is obtained for the optimum specific impulse and power supply. It was found analytically that the optimum power supply is one such that the mass of the propulsion system and power supply are approximately equal to the total mass of propellant used during the entire mission. Finally, it is shown that ACR missions, in general, are optimized using propulsion systems capable of processing 100 kW – 1 MW of power with specific impulses in the range 5,000 – 10,000 s, and have the potential to return asteroids on the order of 103 104 tons. − Nomenclature
    [Show full text]
  • AFSPC-CO TERMINOLOGY Revised: 12 Jan 2019
    AFSPC-CO TERMINOLOGY Revised: 12 Jan 2019 Term Description AEHF Advanced Extremely High Frequency AFB / AFS Air Force Base / Air Force Station AOC Air Operations Center AOI Area of Interest The point in the orbit of a heavenly body, specifically the moon, or of a man-made satellite Apogee at which it is farthest from the earth. Even CAP rockets experience apogee. Either of two points in an eccentric orbit, one (higher apsis) farthest from the center of Apsis attraction, the other (lower apsis) nearest to the center of attraction Argument of Perigee the angle in a satellites' orbit plane that is measured from the Ascending Node to the (ω) perigee along the satellite direction of travel CGO Company Grade Officer CLV Calculated Load Value, Crew Launch Vehicle COP Common Operating Picture DCO Defensive Cyber Operations DHS Department of Homeland Security DoD Department of Defense DOP Dilution of Precision Defense Satellite Communications Systems - wideband communications spacecraft for DSCS the USAF DSP Defense Satellite Program or Defense Support Program - "Eyes in the Sky" EHF Extremely High Frequency (30-300 GHz; 1mm-1cm) ELF Extremely Low Frequency (3-30 Hz; 100,000km-10,000km) EMS Electromagnetic Spectrum Equitorial Plane the plane passing through the equator EWR Early Warning Radar and Electromagnetic Wave Resistivity GBR Ground-Based Radar and Global Broadband Roaming GBS Global Broadcast Service GEO Geosynchronous Earth Orbit or Geostationary Orbit ( ~22,300 miles above Earth) GEODSS Ground-Based Electro-Optical Deep Space Surveillance
    [Show full text]
  • Handbook of Satellite Orbits from Kepler to GPS Michel Capderou
    Handbook of Satellite Orbits From Kepler to GPS Michel Capderou Handbook of Satellite Orbits From Kepler to GPS Translated by Stephen Lyle Foreword by Charles Elachi, Director, NASA Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California, USA 123 Michel Capderou Universite´ Pierre et Marie Curie Paris, France ISBN 978-3-319-03415-7 ISBN 978-3-319-03416-4 (eBook) DOI 10.1007/978-3-319-03416-4 Springer Cham Heidelberg New York Dordrecht London Library of Congress Control Number: 2014930341 © Springer International Publishing Switzerland 2014 This work is subject to copyright. All rights are reserved by the Publisher, whether the whole or part of the material is concerned, specifically the rights of translation, reprinting, reuse of illustrations, recitation, broadcasting, reproduction on microfilms or in any other physical way, and transmission or information storage and retrieval, electronic adaptation, computer software, or by similar or dissimilar methodology now known or hereafter developed. Exempted from this legal reservation are brief excerpts in connection with reviews or scholarly analysis or material supplied specifically for the purpose of being entered and executed on a computer system, for exclusive use by the purchaser of the work. Duplication of this pub- lication or parts thereof is permitted only under the provisions of the Copyright Law of the Publisher’s location, in its current version, and permission for use must always be obtained from Springer. Permis- sions for use may be obtained through RightsLink at the Copyright Clearance Center. Violations are liable to prosecution under the respective Copyright Law. The use of general descriptive names, registered names, trademarks, service marks, etc.
    [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]
  • Orbit and Spin
    Orbit and Spin Overview: A whole-body activity that explores the relative sizes, distances, orbit, and spin of the Sun, Earth, and Moon. Target Grade Level: 3-5 Estimated Duration: 2 40-minute sessions Learning Goals: Students will be able to… • compare the relative sizes of the Earth, Moon, and Sun. • contrast the distance between the Earth and Moon to the distance between the Earth and Sun. • differentiate between the motions of orbit and spin. • demonstrate the spins of the Earth and the Moon, as well as the orbits of the Earth around the Sun, and the Moon around the Earth. Standards Addressed: Benchmarks (AAAS, 1993) The Physical Setting, 4A: The Universe, 4B: The Earth National Science Education Standards (NRC, 1996) Physical Science, Standard B: Position and motion of objects Earth and Space Science, Standard D: Objects in the sky, Changes in Earth and sky Table of Contents: Background Page 1 Materials and Procedure 5 What I Learned… Science Journal Page 14 Earth Picture 15 Sun Picture 16 Moon Picture 17 Earth Spin Demonstration 18 Moon Orbit Demonstration 19 Extensions and Adaptations 21 Standards Addressed, detailed 22 Background: Sun The Sun is the center of our Solar System, both literally—as all of the planets orbit around it, and figuratively—as its rays warm our planet and sustain life as we know it. The Sun is very hot compared to temperatures we usually encounter. Its mean surface temperature is about 9980° Fahrenheit (5800 Kelvin) and its interior temperature is as high as about 28 million° F (15,500,000 Kelvin).
    [Show full text]
  • Multisatellite Determination of the Relativistic Electron Phase Space Density at Geosynchronous Orbit: Methodology and Results During Geomagnetically Quiet Times Y
    JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110, A10210, doi:10.1029/2004JA010895, 2005 Multisatellite determination of the relativistic electron phase space density at geosynchronous orbit: Methodology and results during geomagnetically quiet times Y. Chen, R. H. W. Friedel, and G. D. Reeves Los Alamos National Laboratory, Los Alamos, New Mexico, USA T. G. Onsager NOAA, Boulder, Colorado, USA M. F. Thomsen Los Alamos National Laboratory, Los Alamos, New Mexico, USA Received 10 November 2004; revised 20 May 2005; accepted 8 July 2005; published 20 October 2005. [1] We develop and test a methodology to determine the relativistic electron phase space density distribution in the vicinity of geostationary orbit by making use of the pitch-angle resolved energetic electron data from three Los Alamos National Laboratory geosynchronous Synchronous Orbit Particle Analyzer instruments and magnetic field measurements from two GOES satellites. Owing to the Earth’s dipole tilt and drift shell splitting for different pitch angles, each satellite samples a different range of Roederer L* throughout its orbit. We use existing empirical magnetic field models and the measured pitch-angle resolved electron spectra to determine the phase space density as a function of the three adiabatic invariants at each spacecraft. Comparing all satellite measurements provides a determination of the global phase space density gradient over the range L* 6–7. We investigate the sensitivity of this method to the choice of the magnetic field model and the fidelity of the instrument intercalibration in order to both understand and mitigate possible error sources. Results for magnetically quiet periods show that the radial slopes of the density distribution at low energy are positive, while at high energy the slopes are negative, which confirms the results from some earlier studies of this type.
    [Show full text]
  • Instructor's Guide for Virtual Astronomy Laboratories
    Instructor’s Guide for Virtual Astronomy Laboratories Mike Guidry, University of Tennessee Kevin Lee, University of Nebraska The Brooks/Cole product Virtual Astronomy Laboratories consists of 20 virtual online astronomy laboratories (VLabs) representing a sampling of interactive exercises that illustrate some of the most important topics in introductory astronomy. The exercises are meant to be representative, not exhaustive, since introductory astronomy is too broad to be covered in only 20 laboratories. Material is approximately evenly divided between that commonly found in the Solar System part of an introductory course and that commonly associated with the stars, galaxies, and cosmology part of such a course. Intended Use This material was designed to serve two general functions: on the one hand it represents a set of virtual laboratories that can be used as part or all of an introductory astronomy laboratory sequence, either within a normal laboratory setting or in a distance learning environment. On the other hand, it is meant to serve as a tutorial supplement for standard textbooks. While this is an efficient use of the material, it presents some problems in organization since (as a rule of thumb) supplemental tutorial material is more concept oriented while astronomy laboratory material typically requires more hands-on problem-solving involving at least some basic mathematical manipulations. As a result, one will find material of varying levels of difficulty in these laboratories. Some sections are highly conceptual in nature, emphasizing more qualitative answers to questions that students may deduce without working through involved tasks. Other sections, even within the same virtual laboratory, may require students to carry out guided but non-trivial analysis in order to answer questions.
    [Show full text]
  • Positioning: Drift Orbit and Station Acquisition
    Orbits Supplement GEOSTATIONARY ORBIT PERTURBATIONS INFLUENCE OF ASPHERICITY OF THE EARTH: The gravitational potential of the Earth is no longer µ/r, but varies with longitude. A tangential acceleration is created, depending on the longitudinal location of the satellite, with four points of stable equilibrium: two stable equilibrium points (L 75° E, 105° W) two unstable equilibrium points ( 15° W, 162° E) This tangential acceleration causes a drift of the satellite longitude. Longitudinal drift d'/dt in terms of the longitude about a point of stable equilibrium expresses as: (d/dt)2 - k cos 2 = constant Orbits Supplement GEO PERTURBATIONS (CONT'D) INFLUENCE OF EARTH ASPHERICITY VARIATION IN THE LONGITUDINAL ACCELERATION OF A GEOSTATIONARY SATELLITE: Orbits Supplement GEO PERTURBATIONS (CONT'D) INFLUENCE OF SUN & MOON ATTRACTION Gravitational attraction by the sun and moon causes the satellite orbital inclination to change with time. The evolution of the inclination vector is mainly a combination of variations: period 13.66 days with 0.0035° amplitude period 182.65 days with 0.023° amplitude long term drift The long term drift is given by: -4 dix/dt = H = (-3.6 sin M) 10 ° /day -4 diy/dt = K = (23.4 +.2.7 cos M) 10 °/day where M is the moon ascending node longitude: M = 12.111 -0.052954 T (T: days from 1/1/1950) 2 2 2 2 cos d = H / (H + K ); i/t = (H + K ) Depending on time within the 18 year period of M d varies from 81.1° to 98.9° i/t varies from 0.75°/year to 0.95°/year Orbits Supplement GEO PERTURBATIONS (CONT'D) INFLUENCE OF SUN RADIATION PRESSURE Due to sun radiation pressure, eccentricity arises: EFFECT OF NON-ZERO ECCENTRICITY L = difference between longitude of geostationary satellite and geosynchronous satellite (24 hour period orbit with e0) With non-zero eccentricity the satellite track undergoes a periodic motion about the subsatellite point at perigee.
    [Show full text]
  • Spacecraft Trajectories in a Sun, Earth, and Moon Ephemeris Model
    SPACECRAFT TRAJECTORIES IN A SUN, EARTH, AND MOON EPHEMERIS MODEL A Project Presented to The Faculty of the Department of Aerospace Engineering San José State University In Partial Fulfillment of the Requirements for the Degree Master of Science in Aerospace Engineering by Romalyn Mirador i ABSTRACT SPACECRAFT TRAJECTORIES IN A SUN, EARTH, AND MOON EPHEMERIS MODEL by Romalyn Mirador This project details the process of building, testing, and comparing a tool to simulate spacecraft trajectories using an ephemeris N-Body model. Different trajectory models and methods of solving are reviewed. Using the Ephemeris positions of the Earth, Moon and Sun, a code for higher-fidelity numerical modeling is built and tested using MATLAB. Resulting trajectories are compared to NASA’s GMAT for accuracy. Results reveal that the N-Body model can be used to find complex trajectories but would need to include other perturbations like gravity harmonics to model more accurate trajectories. i ACKNOWLEDGEMENTS I would like to thank my family and friends for their continuous encouragement and support throughout all these years. A special thank you to my advisor, Dr. Capdevila, and my friend, Dhathri, for mentoring me as I work on this project. The knowledge and guidance from the both of you has helped me tremendously and I appreciate everything you both have done to help me get here. ii Table of Contents List of Symbols ............................................................................................................................... v 1.0 INTRODUCTION
    [Show full text]