Motion of deformable planets and stability of their stationary rotations

A mechanical system consisting from N deformable spheres interacting according to the law of gravity is considered as a model of planetary system. Deformations of the viscoelastic spheres are described according to the model of the theory of elasticity of small deformations, the Kelvin-Voigt model of viscous forces, and occur under the action of gravitational fields and fields of centrifugal forces. Approximate equations describing motions of the centers of mass of the spheres and their rotations relative to the centers of mass are constructed by the method of separation of motions on the basis of solving quasistatic problems of the theory of viscoelasticity with allowance made for smallness of sphere deformations. Using the first integral of conservation of the angular momentum of the system relative to its center of mass, the expression for the changed potential energy is obtained with the use of the Routh method. An investigation of stationary rotations is carried out, and it is shown that all of them are unstable, if the number of planets is more than two.