Graphical generation of periodic orbits of tschauner-hempel equations

T HE relative motion of a satellite (follower) with respect to the reference satellite (leader) in a circular orbit is described by autonomous nonlinear differential equations. The linearized equations at the origin are known as Hill–Clohessy–Wiltshire (HCW) equations [1]. The in-plane motion and the out-of plane motion are independent. The latter is a simple sinusoidal motion. The former has also periodic solutions under the so-called Clohessy–Wiltshire condition. Periodic solutions are given by sinusoidal functions, and the relative orbits are ellipses. These orbits are parametrized by two parameters representing size and phase. Because of this simple nature, relative periodic orbits are used for formation flying [2]. Periodic orbits of the nonlinear relative dynamics are also characterized by the parameters of the follower’s elliptic orbit and its inclination angle in an inertial reference frame, and initial conditions for periodic orbits are explicitly given [3]. If the reference orbit is elliptic, the equations of relative motion involve the true anomaly and the radius of the orbit, which are periodic functions. The linearized equations of motion at the origin are known as Tschauner–Hempel (TH) equations. The in-plane motion and the out-of plane motion remain independent. However, the derivation of the state transition matrix in this case is not immediate. In earlier studies [4–6], the true anomaly is used as independent variable, because this simplifies the resulting equations and the transition matrix can be obtained for them. Rendezvous problems have been extensively studied using this representation [7,8]. The out-of-plane motion is a sinusoidal motion, and the inplane motion has periodic solutions [2,9]. The condition, which generalizes the CW condition, for periodic solutions is given in [9]. Periodic orbits of the TH equations are also used for formation flying [2,9,10]. However, they are more complicated compared to ellipses of the HCW equations, and their shape changes with eccentricity. Some of them are very much distorted. Hence, it is not intuitively clear how to design relative orbits. The direct study of the relative dynamics using time as independent variable is given in [11], and a series expansion of the state transition matrix in the eccentricity is obtained in [12]. A closed form of the state transition matrix is recently derived in [13],where a brief history of elliptic rendezvous is also found. The characterization of periodic orbits of the nonlinear relative dynamics as functions of time is extended to the elliptic case, and initial conditions for periodic orbits are explicitly given [14]. The general solution of the relative dynamics based on the orbital elements is obtained in [15], and the relative motion invariant manifold, where the solution lies, is determined. This technical note is concernedwith the generation of all periodic solutions of the in-plane motion of the TH equations. For this purpose, the state transition matrix of [16], which is given as a function of the true anomaly, is used. The solution is parametrized by four constants, and the condition for periodic solution follows from this. As in the case of the HCWequations, size and phase parameters are introduced. The main contributions of this note are summarized as follows. First, the region where all periodic solutions lie is determined by two ellipses, whose size depends on the eccentricity of the leader orbit. Given size and phase parameters, a graphical method to determine the initial point of a periodic orbit at a given true anomaly is proposed. Varying the value of the true anomaly, the whole family of periodic orbits is generated. The time domain versions of the region and periodic orbits are determined by multiplying the radius of the leader orbit. This clarifies how the shape of a periodic orbit changes with eccentricity and why the shapes of some orbits are distorted. As the eccentricity goes to zero, the region of periodic motion shrinks to a single ellipse, which is a periodic orbit of the HCW equations. The size parameter gives the distance between the leader and the follower, and the phase determines the position of the follower. As in [15], the maximum and minimum distances to the leader can be calculated. Hence, our method is useful to design relative orbits with size and position specifications. This note is organized as follows. Section II reviews the equations of relativemotion along an elliptic orbit. Section III gives a geometric method to draw periodic orbits of the TH equations. Finally, Section IV gives conclusions.