Рассматривается проблема получения дифференциальных уравнений, определяющих вековые возмущения орбитальных элементов в многопланетной системе в случае, когда центральная звезда теряет свою массу изотропно, а массы планет могут изменяться анизотропно, что приводит к появлению реактивных сил. В качестве модели многопланетной системы используется классическая задача (n + 1) тел переменной массы, когда n тел движутся вокруг центральной звезды по квазиэллиптическим непересекающимся орбитам и взаимодействуют друг с другом в соответствии с законом всемирного тяготения. Предполагается, что массы тел изменяются с различными скоростями, причем законы изменения масс считаются произвольными заданными функциями времени. Получены дифференциальные уравнения движения тел в оскулирующих элементах апериодического движения по квазиконическим сечениям, соответствующие планетарным уравнениям Лагранжа. Обсуждается алгоритм вычисления возмущающих функций в виде степенных рядов по малым параметрам и получение дифференциальных уравнений, определяющих вековые возмущения орбитальных элементов. Все необходимые символьные вычисления выполняются с использованием системы компьютерной алгебры Wolfram Mathematica.
Индексирование
Scopus
Crossref
Higher Attestation Commission
At the Ministry of Education and Science of the Russian Federation