Смульский И.И. Осесимметричная задача гравитационного взаимодействия N-тел// Математическое моделирование. - 2003, т. 15, № 5, с. 27-36.
УДК 521.11

И.И.Смульский

ОСЕСИММЕТРИЧНАЯ ЗАДАЧА ГРАВИТАЦИОННОГО

ВЗАИМОДЕЙСТВИЯ N-ТЕЛ
1. Введение. До недавнего времени существовало только одно точное решение задачи ньютоновского гравитационного взаимодействия тел. Для случая двух тел она была решена 300 лет назад И.Ньютоном. Десять лет назад Б.Эльмабсут [1] нашел решение задачи для одного случая осесимметрично расположенных на плоскости n-тел, которые равномерно обращаются по окружности вокруг оси симметрии, т.е. при наличии у тел только одной определенной трансверсальной скорости. Совсем недавно Е.А.Гребенников [2] получил решения этой задачи при наличии у тел произвольной радиальной и трансверсальной скоростей. Это же решение получено в нашей работе [3], но другим методом. Если в первой работе используется анализ дифференциальных уравнений движения n-тел, то во второй рассчитывается сила воздействия всех тел на одно, а затем решается дифференциальное уравнение движения одного тела под воздействием суммарной силы. Ниже в более совершенном виде приводится это решение, а также временные характеристики движения n-тел.

В общем случае гравитационного взаимодействия точечных n-тел, тело с массой m_j воздействует на тело с массой m_i (см. фиг.1) силой

, (1)

где G – гравитационная постоянная;

- вектор-расстояние от j-того тела до i–того.

Сила воздействия всех тел на i-тое тело будет равна

. (2)

Тогда ускорение i- того тела согласно 2-ому закону механики запишется

. (3)

где – радиус-вектор i-того тела в инерциальной системе координат.

Уравнения (3) представляют систему 3n обыкновенных нелинейных дифференциальных уравнений второго порядка, которые при задании 3n начальных значений координат n-тел и 3n начальных значений компонент их скоростей представляют задачу n-тел. В общем случае аналитически решена задача только для двух тел, и результаты ее приведены во многих источниках, в том числе [3] и [4]. Ниже представлено решение этой задачи для n-тел, осесимметрично расположенных на плоскости.

2. Сила воздействия на каждое из n-тел. Рассмотрим взаимодействие точечных n-тел, равномерно расположенных по окружности радиуса r (фиг. 2) с угловым интервалом ₁ = 2/n и массы которых одинаковы: m_j = m₁. Тело с массой m_j, отстоящее от тела с массой m₁ на угол _j = 2(j-1)/n, в соответствии с (1) будет воздействовать на него силой, проекция которой на ось x равна

. (4)

Треугольник ОВА – равнобедренный, поэтому 2 + _j = , откуда получаем:

. (5)

Из этого треугольника определяется расстояние АВ = 2r cos Тогда сила (4) запишется так:

, (6)

где j в зависимости от номера тела может изменяться от 2 до n.

С учетом того, что m_j = m₁, просуммируем по индексу j силы воздействия (6) тел от m₂ до m_n на первое тело:

, (7)

где

. (8)

Аналогично (4), сумма проекций на ось y сил воздействия всех тел на первое тело будет:

.

Поэтому силой F_1x полностью определяется воздействие всех тел на тело m₁. Так как эта сила направлена к общему центру масс и обратно-пропорциональна квадрату расстояния r тела до него, то в векторном виде она запишется так:

. (9)

Воздействие всех (n-1) тел на любое тело m_i, будет аналогично воздействию на тело m₁. Поэтому выражением (9) определяется сила воздействия всех тел на тело m_i, где – радиус-вектор тела m_i относительно центра масс O, а функция f(n) – для всех тел одинакова.

Так как тела расположены осесимметрично, то центр их масс находится в начале координат (см. фиг. 2) и, как видно из (8), сила воздействия всех тел на любое из них определяется только расстоянием r до общего центра масс О. Сомножитель f(n) определяется осевой симметрией n тел. В процессе движения число тел не изменяется и отсутствуют причины для нарушения осевой симметрии, т.е. f(n) = const и не зависит от времени.

При сопоставлении силы (8) с силой всемирного тяготения (1) видно, что (8) отличается только этим постоянным коэффициентом f(n) = const. Следовательно, движение каждого из n тел будет происходить под действием центральной силы (8) аналогично движению одного тела, находящегося под воздействием другого. Траектории каждого из тел будут коническими сечениями, в фокусе которых находится центр масс О.

Если дополнительно поместить (n + 1)-ое тело с массой m₀ в центр масс О (фиг. 2), то оно на каждое из n тел будет воздействовать силой

. (10)

Тогда суммарная сила воздействия на каждое из n тел с массой m_j = m₁ будет равна сумме сил (8) и (10):

. (11)

Сила воздействия всех n тел на центральное тело m₀ в силу осевой симметрии равна нулю и оно остается неподвижным.

Из сопоставления (11) с (1) видно, что воздействие совокупности тел на каждое тело с массой m_i эквивалентно воздействию на это тело некоторого тела с массой , находящегося в центре масс O всей совокупности тел.

Как уже отмечалось, сила (11) зависит от расстояния r тела m_i от центра масс O (см. рис. 2). Так как центр масс системы (n+1) тел движется неускоренно, то ускорение тела в инерциальной системе, связанной с ним, будет

, (12)

где

. (13)

Дифференциальное уравнение (12) движения тела m_i относительно центра масс идентично дифференциальному уравнению движения тела m₁ относительно тела с массой , т.е. движение каждого из m_i тел полностью определяется классической проблемой двух тел [4]. В наших работах [3, 5, 6] рассмотрен общий случай задачи двух тел, когда силы зависят не только от относительного расстояния между телами, но и от скорости. Представление временных параметров движения является не тривиальным, поэтому целесообразно их описать более подробно.

3. Уравнения движения. Выразим ускорение в левой части (12) через ускорения в полярной системе координат: на направление радиус-вектора и на направление , перпендикулярное к нему. Составляющую ускорения на перпендикулярное направление называют трансверсальной. Здесь – единичный вектор, а  = – угловая скорость поворота радиус-вектора . Радиальная компонента ускорения будет равна

. (14)

После подстановки этих величин в (12) получаем

. (15)

Приравнивая выражения в правой и левой частях этого равенства при соответствующих векторах и , получаем два уравнения, из которых первое дает:

.

В процессе движения момент количества движения единицы массы, т.е. кинематический момент количества движения h

(16)

остается неизменным. Здесь = r- трансверсальная скорость.

Приравнивая коэффициенты при векторе в (15) с учетом (16) получаем второе уравнение движения:

. (17)

Так как радиальная скорость , то уравнение (17) запишется в виде

. (18)

Интегрируя (18) c учетом начальных условий

при t = t₀ , r = r₀ , _

получаем выражение для радиальной скорости

. (20)

Так как

, (21)

отсюда можно получить интегралы для траектории и для времени движения по траектории

; (22)

где радиальная скорость определяется формулой (20).

В точке траектории, называемой перицентром, взаимодействующие объекты сближаются на минимальное расстояние r_p, при котором их радиальная скорость равна нулю, а трансверсальная принимает наибольшее значение. Преобразуем уравнения для радиальной скорости (20) и траектории (22) к безразмерному виду, приведенному к параметрам в перицентре:

, (23)

, (24)

где

; (25)

– относительный радиус; – относительная радиальная скорость; ₁ – параметр траектории.

Выражение (24) при радиальной скорости (23) заменой легко интегрируется, и при граничном условии при  = 0, получаем уравнение траектории

. (26)

Выражение (26) в полярной системе координат при ₁ = -1 представляет окружность, при –0.5 > – эллипсы, при ₁ = -0.5 – параболу, – гиперболы, при ₁ = 0 – прямую. Итак, траектории (26) описывают сечения конуса плоскостью и обычно [4] приводятся в зависимости от двух параметров:

, (27)

где – фокальный параметр; – эксцентриситет.

Традиционные параметры в соответствии с (26) – (27) выражаются через параметр траектории ₁ так:

, . (28)

Следует отметить, что при взаимодействии заряженных частиц по закону Кулона траектории также описываются выражением (27). Однако при больших скоростях движения необходимо учитывать магнитную составляющую силы, поэтому траектории будут уже другими. Они зависят не только от параметра ₁, но и от скорости частицы в перицентре [3, 5, 6].

4. Время движения по траектории. Интеграл для времени (22) зависит от вида траектории. В соответствии с [3] приведем (22) не к точке перицентра r_p, а к произвольной точке траектории r₀:

,

где

, , , . (29)

В соответствии с определением, связь между параметрами траектории следующая

. (30)

Результат интегрирования зависит от параметра

. (31)

При A > 0 траектория будет гиперболой и время движения по ней запишется:

, (32)

где

, (33)

– радиальная скорость при = 1.

При A < 0 траектория будет эллипсом и время движения по ней в результате интегрирования получаем в виде:

. (34)

При A = 0 траектория является параболой и время движения по ней будет:

. (35)

В случае чисто радиального движения, при h = 0, введем другие приведенные параметры

, (36)

с учетом которых время движения в результате интегрирования получаем в виде:

, (37)

где в соответствии с (20) радиальная скорость

. (38)

Итак, время движения по траектории при осесимметричном взаимодействии многих тел, как и при взаимодействии двух тел, подразделяется на четыре случая движения, три из которых имеют место при наличии трансверсальной составляющей скорости: 1) гиперболический (32); 2) эллиптический (34); 3) параболический (35). Траектории при этом определяются (26) при , и , соответственно. В четвертом случае радиального движения время и закон движения задаются соотношением (37), а траекторией является проходящая через центр прямая линия.

5. Конкретный пример осесимметричного взаимодействия. Если задано число n периферийных тел с массой m₁ и масса центрального тела m₀, в соответствии с (13) определяется ключевой параметр взаимодействия μ₁. Затем согласно (25) рассчитывается параметр траектории α₁ или в соответствии с (29) .

Радиальная скорость вычисляется согласно (23), а трансверсальная . Все величины здесь нормированы к параметрам в перицентре υ_p и r_p. Время движения тел по траекториям определяется в соответствии с (32), (34), (35) и (38), где параметры нормированы к величинам в произвольной точке r_0,, в которой тело имеет скорость υ_t0 и υ_r0 относительно центра О (фиг.2). Движение каждого из n тел происходит при сохранении кинематического момента h = const. Если считать, что на фиг. 2 тела изображены в перицентре, то линия апсид тела m₁ проходит по оси x. Апоцентр его будет расположен на отрицательных значениях оси x. Линия апсид тела m_i наклонена к оси x под углом _0i = 2(i-1)/n. Уравнение траектории тела с номером i определяется выражением (26) при замене φ на φ_i . Центральное тело с массой m₀ будет неподвижно находиться в т. О.

В качестве примера рассмотрим движение четырех тел (n = 4), которые на радиусе r₀ имеют радиальную скорость υ_r0 = 0, при последовательном увеличении υ_t0 от 0 до некоторого наибольшего значения. Вариации траекторий, полученные в [3] в результате численного интегрирования уравнений (3), представлены на фиг. 3. Покажем, что они описываются точными решениями.

В соответствии с (8) величина f(n) = 0.25 + . Согласно (13) параметр взаимодействия μ₁ = - Gm₁(0,25 + ), а параметр траектории α₁⁰ = Gm₁(0.25 + )/(r₀υ_t0). На фиг. 3 представлены траектории в декартовых координатах (ось x – горизонтальна), которые смещены относительно оси симметрии и нормированы к первоначальному расстоянию между ближайшими телами. Связь координат x,y с полярными координатами (r, φ) следующая:

x/a = 0.5 + (r/a)cosφ; y/a = 0.5 + (r/a)sinφ. (39)

Заменяя на в соответствии с (30), выразим уравнение траектории (26) через r/a

, (40)

где i = 1, 2, 3, 4 – номера тел, а φ_0i – угол расположения линии апсид. Так как на фиг.3, в отличие от фиг.2, горизонтальная ось x расположена между телами, то

φ_0i = π/4 + (i-1) π/2. (41)

Согласно (40) при (φ - φ_0i) = 0 частицы находятся в перицентре. При (φ - φ_0i) = π частицы будут в апоцентре

. (42)

При выборе в качестве начального радиуса радиус апоцентра r₀ = r_a из (42) получаем выражение для радиуса перицентра

. (43)

Так как в точке r₀ в соответствии с (29) трансверсальная скорость υ_t0 = определяется величиной α₁⁰ при постоянных μ₁ и r₀, то разные значения υ_t0 будем задавать вариацией α₁⁰. В случае чисто радиального движения υ_t0 = 0 и h = 0. Тогда радиальная скорость изменяется по закону (38), а трансверсальная υ_t = h/r=0, т.е. все четыре тела движутся по радиальным линиям (40) начиная с до соприкосновения в начале координат (см.фиг.3,1). Закон движения, т.е. зависимость r(t) определяется соотношением (37).

При небольшой трансверсальной скорости, определяемой α₁⁰ = -17.17 радиус перицентра, согласно (43) r_p / r₀ = 0.03 , а параметр траектории α₁ = -0.515. То есть орбита будет вытянутым эллипсом. Согласно (40) при (φ - φ_0i) = π частицы будут находиться в апоцентре , с увеличением угла (φ - φ_0i) до 2π будут двигаться к центру и достигают перицентров, как показано на фиг. 3,2.

При увеличении трансверсальной скорости (α₁⁰ = -3.37) тела движутся по более наполненным эллипсам и в соответствии с (43) проходят на большем расстоянии от центра r_p = 0.174r₀ (см.фиг.3,3). Последовательное увеличение трансверсальной скорости : α₁⁰ = -1.88; -1.058 приводит к дальнейшему увеличению r_p = 0.362r₀ и r_p = 0.896r_0,, т.е. уменьшается эксцентриситет орбит (см. фиг. 3,4 и 3,5), соответственно. И при α₁⁰ = -1 радиус перицентра, согласно (44) r_p = r₀, т.е. все тела движутся по окружности (см.рис.2,6) вокруг общего центра масс. Параметр траектории α₁ = -1, т.е. в этом случае совпадает с α₁⁰.

При дальнейшем увеличении скорости υ_t0, которая соответствует α₁⁰ = -0.559, согласно (43) r_p = 8.47 r₀. Так как r_p > r₀, то эта точка уже является апоцентром. Тела уже не сближаются (см. фиг.3,7), а удаляются друг от друга и при φ – φ₀ = π достигают апоцентра r_a = 8.47 r₀. Как и в предыдущем случае, параметр траектории α₁⁰ отнесен к начальному радиусу r₀, который в данном случае является радиусом перицентра, поэтому α₁⁰ совпадает с α₁. С увеличением υ_t0, которая соответствует α₁⁰ = α₁ = -0.5, радиус апоцентра согласно (43) возрастает и стремится к бесконечности, т.е. тела движутся по параболе (см. фиг. 3,8). При еще большей трансверсальной скорости (см. α₁⁰ = -0.47, на фиг. 3,9) тела движутся по гиперболам.

Представленные на фиг.3 траектории отражают всевозможные типы движения при осесимметричном взаимодействии не только четырех тел, но и любого числа n тел. При приближенном аналитическом или численном решении системы уравнений (3) эллиптические орбиты получаются незамкнутые, т.е. они поворачиваются и перицентр прецессирует. С увеличением точности счета величина прецессии уменьшается, но полностью избавиться от нее невозможно. Поэтому при приближенном решении задачи гравитационного взаимодействия тел всегда будет получаться прецессия орбит, даже в случаях, когда ее нет.

Следует отметить, что при взаимодействии наэлектризованных тел (некулоновское взаимодействие) орбиты не являются коническими сечениями и их перицентры прецессируют [3, 5] несмотря на то, что движения проходят при сохранении кинематического момента h.

Итак, осесимметричная задача n-тел сводится к задаче взаимодействия двух тел. Движение каждого из n-тел происходит по неизменной в пространстве траектории, которая принадлежит к семейству кривых конического сечения и определяется параметром ₁. Практическое значение данной задачи заключается в возможности определения реальной погрешности вычислительных алгоритмов проблемы многих тел.

Литература

1. Elmabsout B. Sur l'existence de certaines configurations d'equillibre relatif dans le probleme des n corps. //Celestial Mechanics and Dynamical Astronomy. 1988, v.4, N1, p.131-151.

2. Гребеников Е.А. Существование точных симметричных решений в плоской ньютоновской проблеме многих тел//Математическое моделирование.-1998.-Т.10, N8. - С.74-80.

3. Смульский И.И. Теория взаимодействия.- Новосибирск: Из-во Новосибирского ун-та, ННЦ ОИГГМ СО РАН. 1999. - 294 с.

4. Справочное руководство по небесной механике и астродинамике. – Под рук. Г.Н. Дубошина. Изд. 2-е, доп. и перераб. – М.: Наука.-1976.-862 с.

5. Смульский И.И. Траектории при взаимодействии двух тел, зависящем от относительных расстояния и скорости // Мат. моделирование.– 1995.– Т.7.– N.7.– С. 111-125.

6. Смульский И.И. Электромагнитное и гравитационное воздействие (нерелятивистские трактаты).– Новосибирск: ВО “Наука”. Сиб. издат. фирма, 1994.– 224 с.

г.н.с. Института криосферы Земли СО РАН, д. ф.-м. н.

10.06.02 г. Смульский Иосиф Иосифович

625000, Тюмень, а/я 1230, ИКЗ СО РАН

Тел. (8-345-2)-27-35-18 Fax:(8-345-2)-25-11-53

Ф

иг. 1. Сила гравитационного воздействия тела m_j на тело m_i.
Фиг. 2. Взаимодействие осесимметрично расположенных n-тел.


Фиг. 3. Траектории движения осесимметрично расположенных четырех тел при вариации начальной трансверсальной скорости , определяемой для случаев 2÷9 параметром ₁⁰ (в скобках указан параметр траектории ₁):

1 – = 0; 2 – -17.17 (-0.515); 3 – -3.37 (-0.587); 4 – -1.88 (-0.681); 5 – -1.058 (-0.948); 6 – -1 (-1); 7 – -0.559 (-0.559); 8 – -0.5 (-0.5); 9 – -0.47 (-0.47).

Список принятых обозначений
F – сила воздействия на тело;

G – гравитационная постоянная;

h = υ_t r – кинематический момент количества движения;

m – масса частицы или тела;

n – число взаимодействующих тел;

– радиус-вектор тела относительно центра масс;

r_p – радиус перицентра;

– приведенное расстояние между телами;

– радиальная и трансверсальная скорости тел;

– радиальная и трансверсальная скорости тела в начальной точке r₀;

, – параметр траектории;

– параметр траектории при отнесении к точке с радиусом r₀;

– эксцентриситет траектории;

– константа гравитационного взаимодействия;

УДК 521.11

И.И.Смульский

ОСЕСИММЕТРИЧНАЯ ЗАДАЧА ГРАВИТАЦИОННОГО

ВЗАИМОДЕЙСТВИЯ N-ТЕЛ

Автореферат

Рассмотрено гравитационное взаимодействие n-тел, осесимметрично расположенных на плоскости. Рассчитана сила воздействия на каждое тело. Показано, что она направлена к центру масс и обратно пропорциональна расстоянию до него. Приведены зависимости для скорости, траектории и времени движения по траектории для всех возможных случаев движения. В частном случае движения четырех тел полученные точные решения подтверждены численными решениями этой задачи. Точное решение дает ответы на ряд вопросов, которые не могут быть получены в результате приближенных решений. Кроме того, оно позволяет тестировать приближенные методы решения проблемы многих тел.

Ил. 3. Библиогр. 6 назв.

Joseph J. Smulsky

The Axisymmetric Problem of Gravitational

Interaction of N-bodies
It is considered the gravitational interaction of n-bodies, which axially symmetric located on a plane. The force of action on each body is calculated. It is shown, that force is directed to a center of masses and is inversely to a distance to it. The expressions for a velocity, trajectory and time of moving on trajectory are received for all possible cases of interaction. The obtained exact solutions are confirmed by numerical solutions of this task in the specific case of four bodies. The exact solution gives the answers to a many of problems, which cannot be obtained by approximate solutions. Besides it allows to test approximate methods of a solution of a many bodies problem.