Abstract
Numerical methods have proven to be reliable tools in generating long-term bounded relative trajectories in nontrivial dynamical environments. Yet, none of the existing procedures have a systematic approach to come up with all of the bounded relative trajectories associated with a variety of satellites in Earth orbit. In this paper, such a systematic procedure is developed based on the assumption that any trajectory, except at critical inclinations, is either periodic or quasi periodic in the four-dimensional Routh reduced system describing the motion of a mass particle about an axisymmetric body. This allows key design parameters to be identified, such as the nodal period and the node drift per nodal period, which are later used to constrain a differential corrector scheme looking for families of quasi-periodic invariant tori that yield bounded relative motion in the Earth-centered inertial frame. Numerical simulations show that the computed solutions are good candidates for cluster flight missions at a large range of altitude and inclination values.