Аннотация:Предложена модель Солнечной системы, состоящая из N однородных вязкоупругих шаров взаимодействующих по закону Всемирного тяготения. С каждым шаром связывается система координат, относительно которой количество движения шара и интеграл от ротора поля перемещений равны нулю. Деформации шаров описываются линейной теорией вязкоупругости. Под действием собственного вращения и приливных гравитационных сил шаровая планета меняет свою форму: происходит «сплющивание» планеты в направлении вектора ее угловой скорости и формирование приливных горбов по линиям, соединяющим центр планеты с центрами остальных планет. Предполагается, что приливные гравитационные, центробежные и упругие силы приводят к малому изменению шаровой формы планеты. Методом разделения движений получена система 6N обыкновенных дифференциальных уравнений, описывающая поступательно-вращательные движения всех шаров. В предложенной модели имеет место диссипация энергии, источником которой являются внутренне вязкие силы каждой планеты. В результате деформаций планет в законе всемирного тяготения появляются малые консервативные поправки. Система уравнений допускает первый интеграл – закон сохранения момента количеств движения относительно барицентра системы. Доказано, что стационарные движения системы соответствуют ее вращению с постоянной угловой скоростью относительно барицентра как твердого тела, когда центры всех деформированных шаров расположены в плоскости, ортогональной вектору момента количеств движения, а планеты находятся в стационарном деформированном состоянии. Все стационарные вращения неустойчивы, если число шаров N>2. В случае двух шаров возможно существование двух стационарных вращений, одно из которых с меньшим расстоянием между центрам неустойчиво, а второе с большим расстоянием – устойчиво.