Abstract A frozen orbit is beneficial for observation owing to its stationary apsidal line. The traditional gravitational field model of orbits only considers the main zonal harmonic terms J 2 and limited high-order terms, which cannot meet stringent demands all missions. In this study, expanded 15 Hamiltonian canonical form described by Delaunay variables used. coefficients Earth are chosen as...