Traditional analytical theories of celestial mechanics are not well-adapted when dealing with highly elliptical orbits. On the one hand, analytical solutions are quite generally expanded into power series of the eccentricity and so limited to quasi-circular orbits. On the other hand, the time-dependency due to the motion of the third body (e.g. Moon and Sun) is almost always neglected. We propose several tools to overcome these limitations. Firstly, we have expanded the third-body disturbing function into a finite polynomial using Fourier series in multiple of the satellite's eccentric anomaly (instead of the mean anomaly) and involving Hansen-like coefficients. Next, by combining the classical Brouwer-von Zeipel procedure and the time-dependent Lie-Deprit transforms, we have performed a normalization of the expanded Hamiltonian in order to eliminate all the periodic terms. One of the benefits is that the original Brouwer solution for J2 is not modified. The main difficulty lies in the fact that the generating functions of the transformation must be computed by solving a partial differential equation, involving derivatives with respect to the mean anomaly, which appears implicitly in the perturbation. We present a method to solve this equation by means of an iterative process. Finally we have obtained an analytical tool useful for the mission analysis, allowing to propagate the osculating motion of objects on highly elliptical orbits (e>0.6) over long periods efficiently with very high accuracy, or to determine initial elements or mean elements. Comparisons between the complete solution and the numerical simulations will be presented.

When dealing with satellites orbiting a central body on a highly elliptical orbit, it is necessary to consider the effect of gravitational perturbations due to external bodies. Indeed, these perturbations can become very important as soon as the altitude of the satellite becomes high, which is the case around the apocentre of this type of orbit. For several reasons, the traditional tools of celestial mechanics are not well adapted to the particular dynamic of highly elliptical orbits. On the one hand, analytical solutions are quite generally expanded into power series of the eccentricity and therefore limited to quasi-circular orbits [17, 25]. On the other hand, the time-dependency due to the motion of the third-body is often neglected. We propose several tools to overcome these limitations. Firstly, we have expanded the disturbing function into a finite polynomial using Fourier expansions of elliptic motion functions in multiple of the satellite's eccentric anomaly (instead of the mean anomaly) and involving Hansen-like coefficients. Next, we show how to perform a normalization of the expanded Hamiltonian by means of a time-dependent Lie transformation which aims to eliminate periodic terms. The difficulty lies in the fact that the generator of the transformation must be computed by solving a partial differential equation involving variables which are linear with time and the eccentric anomaly which is not time linear. We propose to solve this equation by means of an iterative process.