Motion of a satellite with a variable mass distribution in a central field of Newtonian attraction

Within the framework of the so-called satellite approximation, configurations of the relative equilibrium are built and their stability is analyzed. In this case the elliptic Keplerian motion of the satellite/the spacecraft tight group mass center is predefined. The attitude motion of the system does not influence its orbital motion. The principal central axes of inertia are assumed to move as a rigid body. Simultaneously masses of the body can redistribute in a way such that the values of moments of inertia can change. Thus, all configurations can perform pulsing motions changing it own dimensions.
One obtains a system of equations of motion for such a compound satellite. It turns out that the resulting system of equations is similar to the well-known equation of V.V.Beletsky for the satellite in elliptic orbit planar oscillations. We use true anomaly as an independent variable as it is in the Beletsky equation. It turned out that there are planar pendulum-like librations of the whole system which may be regarded as perturbations of the mathematical pendulum.
One can introduce action-angle variables in this case and can construct the dynamics of mappings over the non-autonomous perturbation period. As a result, one is able to apply the well-known Moser theorem on an invariant curve for twisting maps of annulus. After that one can get a general picture of motion in the case of the system planar oscillations. So, the whole description in the paper splits into two topics: (a) general dynamical analysis of the satellite planar attitude motion using KAM theory; (b) construction of periodic solutions families depending on the perturbation parameter and rising from equilibrium as the perturbation value grows. The latter families depend on the parameter of the perturbation and are absent in the non-perturbed problem.