Effect of 3rd-degree gravity harmonics and Earth perturbations on lunar artificial satellite orbits

In a previous work we studied the effects of (I) the J (2) and C (22) terms of the lunar potential and (II) the rotation of the primary on the critical inclination orbits of artificial satellites. Here, we show that, when 3rd-degree gravity harmonics are taken into account, the long-term orbital behavior and stability are strongly affected, especially for a non-rotating central body, where chaotic or collision orbits dominate the phase space. In the rotating case these phenomena are strongly weakened and the motion is mostly regular. When the averaged effect of the Earth's perturbation is added, chaotic regions appear again for some inclination ranges. These are more important for higher values of semi-major axes. We compute the main families of periodic orbits, which are shown to emanate from the 'frozen eccentricity' and 'critical inclination' solutions of the axisymmetric problem ('J (2) + J (3)'). Although the geometrical properties of the orbits are not preserved, we find that the variations in e, I and g can be quite small, so that they can be of practical importance to mission planning.