Invariant relative orbits for satellite constellations: A second order theory

Working with the mean orbit elements, the secular drift of the longitude of the ascending node and the sum of the argument of perigee and mean anomaly are set equal between two neighboring orbits. By having both orbits drift at equal angular rates on the average, they will not separate over time due to the influence of the perturbative effects of the asphericity of the Earth, as is considered in this work. The problem is stated. The expressions for the time rate of change of the longitude of the ascending node and the sum of the argument of perigee and mean anomaly in terms of the Delaunay canonical elements are obtained. The expressions for the second order conditions that guarantee that the drift rates of two neighboring orbits are equal on the average are derived.