Analytical investigations of quasi-circular frozen orbits in the Martian   gravity field
  Frozen orbits are always important foci of orbit design because of their valuable characteristics that their eccentricity and argument of pericentre remain constant on average. This study investigates quasi-circular frozen orbits and examines their basic nature analytically using two different methods. First, an analytical method based on Lagrangian formulations is applied to obtain constraint conditions for Martian frozen orbits. Second, Lie transforms are employed to locate these orbits accurately, and draw the contours of the Hamiltonian to show evolutions of the equilibria. Both methods are verified by numerical integrations in an 80\times80 Mars gravity field. The simulations demonstrate that these two analytical methods can provide accurate enough results. By comparison, the two methods are found well consistent with each other, and both discover four families of Martian frozen orbits: three families with small eccentricities and one family near the critical inclination. The results also show some valuable conclusions: for the majority of Martian frozen orbits, argument of pericentre are kept at 270 degrees because J3 has the same sign with J2; while for a minority of ones with low altitude and low inclination, argument of pericentre are possible to be kept at 90 degrees because of the effect of the higher degree odd zonals; for the critical inclinations cases, argument of pericentre can also be kept at 90 degrees. It is worthwhile to note that there exist some special frozen orbits with extremely small eccentricity, which could provide much convenience for reconnaissance. Finally, the stability of Martian frozen orbits is estimated based on the trace of the monodromy matrix. The analytical investigations can provide good initial conditions for numerical correction methods in the more complex models. 
Xiaodong Liu, Hexi Baoyin, Xingrui Ma
astro-ph.EP math-ph math.MP physics.space-ph