Lunar Frozen Orbits

Based upon the derivation of the equations of frozen orbits about the Earth, an analytical formulation is presented for lunar frozen orbits (LFO) in a wide altitude range. Global examination is made of combinations of inclination, semimajor axis, and eccentricity to zero out the mean eccentricity and argument of periapsis drift. Included in this development are the effects of both the lunar gravity potential and third body perturbations. Unlike the Earth frozen orbits, with the major perturbation due to the uniform 52 through J5 gravity terms, lunar frozen orbits are restricted by the distribution of the lunar mascons. With accurate modeling, they exhibit a predictive pattern that drives the overall mission design. In addition to a derivation for analytical expressions, a fill1 numerical simulation is used for verification and examination of secular drift from all third body effects. This presentation addresses a selection of usable lunar orbits that reduce or eliminate the need for stationkeeping while maintaining usable science or communication/navigation orbit parameters.