Spacecraft Formation Control Using Analytical Integration of Gauss’ Variational Equations

Spacecraft Formation Control Using Analytical Integration of Gauss’ Variational Equations

SPACECRAFT FORMATION CONTROL USING ANALYTICAL INTEGRATION OF GAUSS’ VARIATIONAL EQUATIONS Mohamed Khalil Ben Larbi, Enrico Stoll Institute of Space Systems Technische Universitaet Braunschweig Hermann-Blenk-Str. 23, 38108 Braunschweig, Germany ABSTRACT Apollo space program, which had the Lunar Excursion Mod- ule and the Command and Service Module being assembled This paper derives a control concept for far range Formation in orbit. The purpose is not to correct the Earth relative orbit Flight (FF) applications assuming circular reference orbits. itself during this maneuver, but rather to adjust and control The paper focuses on a general impulsive control concept for the relative orbit between two vehicles. The relative distance FF which is then extended to the more realistic case of non- is decreased to zero in a very slow and controlled manner dur- impulsive thrust maneuvers. The control concept uses a de- ing the docking maneuver [1]. scription of the FF in relative orbital elements (ROE) instead of the classical Cartesian description since the ROE provide a The modern-day focus of s/c formation flying has ex- direct access to key aspects of the relative motion and are par- tended to maintain a formation of various s/c. Several forma- ticularly suitable for relative orbit control purposes and col- tion flying missions are currently operating or in the deisgn lision avoidance analysis. Although Gauss’ variational equa- stage: synthetic aperture interferometers for Earth obser- tions have been first derived to offer a mathematical tool for vation (e.g. TanDEM-X/TerraSAR-X), dual s/c telescopes processing orbit perturbations, they are suitable for several (e.g. XEUS) and laser interferometer for the detection of different applications. If the perturbation acceleration is due gravitational waves (e.g. LISA). It became obvious that to a control thrust, Gauss’ variational equations show the ef- the formation flying concept overcomes significant techni- fect of such a control thrust on the keplerian orbital elements. cal challenges and even avoids financial limitations. Indeed, Integrating the Gauss’ variational equations offers a direct re- the distribution of sensors and payloads among several s/c lation between velocity increments in the local vertical lo- allows higher redundancy, flexibility, and new applications cal horizontal (LVLH) frame and the subsequent change of that would not be achievable with a single s/c [2]. Recently, keplerian orbital elements. For proximity operations, these formation flying with non-cooperative objects has emerged equations can be generalized from describing the motion of as a new focus, motivated through the increasing number of single spacecraft to the description of the relative motion of such objects especially in low Earth Orbit. The exploitation two spacecraft. This will be shown for impulsive and finite- of satellite-based resources has led to an ever increasing num- duration maneuvers. Based on that, an analytical tool to es- ber of non-cooperative objects such as defunct satellites and timate the error induced through impulsive planning is pre- rocket upper stages termed as space debris. A collisional sented. The resulting control schemes are simple and effective cascading effect has been postulated by Kessler in the late and thus, also suitable for on-board implementation. Simula- seventies [3] and seems more real than ever since 2009 when tions show that the proposed concept improves the timing of two intact artificial satellites, Iridium 33 (operational) and the thrust maneuver executions and thus reduces the residual Kosmos 2251 (out of service) collided distributing debris error of the formation control. across thousands of cubic kilometers. Using the the NASA long-term orbital debris projection model (LEGEND), Liou Index Terms— Gauss’ variational equations, relative or- showed that, besides already implemented mitigation mea- bital elements, impulsive thrust, continuous thrust sures, the annual active debris removal (ADR) of 5-10 prior- itized objects from orbit, is required in order to stabilize the LEO environment [4]. Those results have been backed by 1. INTRODUCTION several other studies [5] and meanwhile are widely accepted in the space debris community. The theory of spacecraft (s/c) formation flying has become the focus of considerably extensive research and development This evolution has led to a substantial research effort to effort during the last decades. Earlier design techniques ad- develop a theory that could simply and explicitly address dressed rendezvous and docking missions such as those of the the relative motion and collision avoidance issue in control design. Thus the development of the upcoming theories, 2. DYNAMICS OF RELATIVE MOTION such as the relative orbital elements (ROE) and Eccentric- ity/Inclination (E/I) vector separation, originally developed First of all we define some notations adopted in this paper. for collocation of geostationary satellites [6], were conducted. The motion of a single s/c orbiting the Earth is described in the Earth Centered Inertial (ECI) frame. The relative mo- The control of satellite formation is performed by the ac- tion of two s/c orbiting the Earth is described in the radial- tivation of on-board thrusters. Typically, impulsive control tangential-normal (RTN) frame (Fig. 1). (very short duration thrust) is preferred to continuous control The ∆(:) operator indicates arithmetic differences be- (finite-duration thrust). This is due to historical limitation on tween absolute Cartesian or Orbital parameters. The δ(:) propulsion technologies, typical payload requirements espe- operator indicates the relative Cartesian position and veloc- cially for scientific missions, and the simplicity of impulsive ity in the RTN Frame. It refers generally to a non-linear planning often allowing pure analytical maneuver design. Re- combination of the absolute Cartesian/orbital parameters. cent advances in propulsion and computer technologies sug- The s/c about which all other s/c motions are referenced gest a deeper study of continuous planning. This approach is called the Client and is denoted with subscript (• ). The is not only justified by the precision of continuous maneuvers 1 second s/c, referred to as Servicer, is to fly in formation with planning (the finite duration is explicitly addressed and no im- the Client and is denoted with the subscript (• ). Absolute pulsive assumptions are made) but also because of the typical 2 Cartesian and orbital parameters without subscript are to be advantages of low thrust propulsion systems, such as reduced understood as Client parameters. mass, limited required power, and variable exhaust velocity. A vast amount of literature exists on formation recon- 2.1. Linearized Equations of Relative Motion figuration maneuvers with impulsive thrust mainly modeled with Clohessy-Wiltshire and Lawden’s equations of relative The Client position vector in the Earth-centered inertial frame motion. The Gauss’ variational equations (GVE) of motion (ECI) is noted r1. The relative orbit will be described in the however offer an ideal mathematical framework for design- rotating local orbital frame RTN in terms of the Cartesian co- T ing impulsive control laws [7]. These equations have been ordinate vector δr = x y z . We introduce dimension- extensively used in the last decades for absolute orbit keeping of single s/c, but have only recently been exploited for forma- tion flying control in LEO as introduced by Schaub et. al [8] and subsequently by Vaddi et al. [9], Breger et al. [10] and D’amico [11]. The reason for such slow development is that GVE describe the effect of control acceleration on the time derivative of the Keplerian orbital elements which were nor- mally used to parametrize the motion of a single s/c but not the relative motion of a formation [10]. In this paper, we build upon the previously mentioned ref- erences and give the following original contributions. Firstly, a comprehensive literature survey of the relative motion parametrization and Gauss’ variational equations for relative Fig. 1. Illustration of a s/c formation in RTN frame motion is presented demonstrating the convenience of ROE and GVE-based maneuver planning as opposed to Cartesian less spatial coordinates and a dimensionless time τ via the parametrization. Secondly, Gauss’ variational equation for equations: relative motion using finite duration thrust are derived for δr T a specific set of ROE. Thirdly, an impulsive maneuver plan δρ = = x¯ y¯ z¯ and dτ = ndt ; (1) based on[2] is extended to the general case of non-zero rel- r1 ative semi major axis and finally translated to the case of finite-duration thrust. The paper is organized as follows. δr T In section 2, an overview of the theory of relative motion, δρ = = x¯ y¯ z¯ (2a) r including ROE, is presented. In section 3, an overview of 1 the GVE and their application in relative orbit control is dτ = ndt ; (2b) presented and the integrated GVE are derived. Section 4 is with n the mean motion. The differentiation with respect to dedicated to the study of the effects of impulsive and finite the independent variable τ is written here as duration thrust on ROE. Subsequently the impulsive and fi- nite duration maneuver schemes are derived in section 5 and d () ()0 ≡ : (3) verified via numerical simulation. dτ The dimensionless state vector is then noted The amplitudes and phases of the in-plane and out-of-plane relative motion oscillations are T δρ δρ 0 δx = : (4) q q 2 2 2 2 c34 = c3 + c4 c56 = c5 + c6 (8) The Clohessy-Wiltshire (CW) equations of motion take a very c c elegant numerically advantageous form if written in a non- ' = arctan 4 θ = arctan 6 : (9) dimensional form [12]: c3 c5 We can easily see from (7) that the Servicer moves in an x¯00 − 2¯y0 − 3¯x = 0 (5a) y¯00 + 2¯x0 = 0 (5b) z¯00 +z ¯ = 0 : (5c) Note that these equations of motion are valid only if : • the Client orbit is circular, • δρ r1, • the Client and Servicer s/c have a pure Keplerian mo- tion.

View Full Text

Details

  • File Type
    pdf
  • Upload Time
    -
  • Content Languages
    English
  • Upload User
    Anonymous/Not logged-in
  • File Pages
    11 Page
  • File Size
    -

Download

Channel Download Status
Express Download Enable

Copyright

We respect the copyrights and intellectual property rights of all users. All uploaded documents are either original works of the uploader or authorized works of the rightful owners.

  • Not to be reproduced or distributed without explicit permission.
  • Not used for commercial purposes outside of approved use cases.
  • Not used to infringe on the rights of the original creators.
  • If you believe any content infringes your copyright, please contact us immediately.

Support

For help with questions, suggestions, or problems, please contact us