A second-order analytic analysis for frozen orbits

An analytical method for studying characteristics of frozen orbits using mean orbital elements is presented. By removing some of the approximations of Brouwer and Kosai (1959) and using the method of integration by parts, an analytic solution is obtained accurate to the second order in the sense of expansion of the perturbing potential of an aspherical approximation. The solution applies to orbits at and near the critical inclination as well as to frozen orbits and, of course, to unfrozen orbits. The analysis reveals that a frozen orbit can be interpreted as an orbit in which movements of the long-periodic part of the rate of change of periapsis cancels the secular part exactly at any time.