In this work, the problem of resonance caused by some gravitational potentials due to Mercury and a third body, namely the Sun, together with some non-gravitational perturbations, specifically coronal mass ejections and solar wind in addition to radiation pressure, are investigated. Some simplifying assumptions without loss of accuracy are employed. The considered force model is constructed. Then the Delaunay canonical set is introduced. The Hamiltonian of the problem is obtained then it is expressed in terms of the Deluanay canonical set. The Hamiltonian is re-ordered to adopt it to the perturbation technique used to solve the problem. The Lie transform method is surveyed. The Hamiltonian is doubly averaged. The resonance capture is investigated. Finally, some numerical simulations are illustrated and are analyzed. Many resonant inclinations are revealed.

Frozen orbit is an attractive option for orbital design owing to its characteristics (its argument of pericenter and eccentricity are kept constant on an average). Solar sails are attractive solutions for massive and expensive missions. However, the solar radiation pressure effect represents an additional force on the solar sail that may greatly affect its orbital behavior in the long run. Thus, this force must be included as a perturbation force in the dynamical model for more accuracy. This study shows the calculations of initial conditions for a lunar solar sail frozen orbit. The disturbing function of the problem was developed to include the lunar gravitational field that is characterized by uneven mass distribution, third body perturbation, and the effect of solar radiation. An averaging technique was used to reduce the dynamical problem to a long period system. Lagrange planetary equations were utilized to formulate the rate of change of the argument of pericenter and eccentricity. Using the reduced system, frozen orbits for the Moon sail orbiter were constructed. The resulting frozen orbits are shown by two 3Dsurface (semimajor, eccentricity, inclination) figures. To simplify the analysis, we showed inclination–eccentricity contours for different values of semi-major axis, argument of pericenter, and values of sail lightness number.