Artificial satellites orbits in 2:1 resonance: GPS constellation

In this work, the resonance problem in the artificial satellites motion is studied. The development of the geopotential includes the zonal harmonics J₂₀ and J₄₀ and the tesseral harmonics J₂₂ and J₄₂. Through an averaging procedure and successive Mathieu transformations, the order of dynamical system is reduced and the final system is solved by numerical integration. In the simplified dynamical model, three critical angles are studied. The half-width of the separatrix is calculated through a linearized model which describes the behavior of the dynamical system in a neighborhood of each critical angle. Through the resonance overlap criterion the possible regular and irregular motions are investigated by the time behavior of the semi-major axis, argument of perigee and eccentricity. The largest Lyapunov exponent is used as tool to verify the chaotic motion.