Evolution of almost circular orbits of satellites under the action of noncentral gravitational field of the Earth and lunisolar perturbations

Based on the results of paper [1] by G.V. Mozhaev, joint perturbations produced by nonsphericity of the Earth and by attraction of the Moon and the Sun are investigated using the method of averaging. Arbitrary number of spherical harmonics was taken into account in the force function of the Earth’s gravitational filed, and only the principal term was retained in the perturbing function of the Sun. In the perturbing function of the Moon two parallactic terms were considered in addition to the dominant term. The flight altitude was chosen in such a way that perturbations produced by the Sun and Moon would have the second order of smallness relative to the polar oblateness of the Earth. As a result, the formulas for calculation of satellite coordinates are derived that give a high precision on long time intervals.