Потенциальная функция в интегральный критерии устойчивости (экстремальное свойство) синхронных движений для небесномеханических орбитальных систем
Как нетрудно видеть, свободные орбитальные системы, описанные в п. 1 § 6 гл. 12, соответствуют рассматриваемым в небесной механике системам, состоящим из некоторого числа к свободных твердых тел В і, ..., Bh (в простейшем случае — материальных точек), движущихся по замкнутым траекториям относительно некоторого центрального тела В0 (в простейшем случае — также материальной точки),, обладающего обычно значительно большей массой (рис. 44). Вращательные координаты тел <р« (f + а») при отсутствии взаимодействия между телами (в более общем случае — между соответствующими степенями свободы) определены вследствие автономности с точностью до произвольных начальных фаз а, (см. равенства (6.1)1 гл. 12). Основные уравнения (4.6) или (5.8) Рис. 44. гл. 12, из которых определяются зна
чения этих фаз в возможных синхронных движениях, имеют вид
f,(«!.............. = + = 0 (*=l,...,fc) (2.1)
(как и в § 6 гл. 12, далее рассматриваем случай, когда ф*' ~ ф* — орбитальные координаты, число которых равно числу тел к). Здесь Ло — усредненная за период Т = 2л/и синхронного движения функция Лагранжа системы связей, вычисленная на порождающих траекториях, а А, — усредненные неконсервативные, обобщенные силы, являющиеся в данном случае диссипативными силами, е, — так называемая крутизна частотной характеристики.
Предположим сначала, что существует потенциал усредненных неконсервативных сил В, т. е. функция B — B(aі, ..., ак), удовлетворяющая соотношениям
Тогда, как указывается в § 5 гл. 12, за потенциальную функцию D может быть принято выражение
где о = sgn е„ — величина, определяемая характером анизохро - низма объектов, который предполагается одинаковым для всех s. Согласно интегральному критерию, устойчивые синхронные движения могут соответствовать точкам грубых минимумов функции D но разностям фаз а, — ак.
Для рассматриваемой орбитальной системы кинетическая энергия несомых связей Т111) и потенциальная энергия несущих связей Па) рарны нулю, а потенциальная энергия несомых связей П(Н) = — С7(П), где t/(1I) — потенциал сил тяготения. Поэтому, согласно сказанному в § 5 гл. 12, в случае малости колебаний несущего тела выражение (2.3) для потенциальной функции представится в форме
D = [<(гт)> - <(t/‘n>)> - В]с, (2.4)
где Тт — кинетическая энергия центрального (несущего) тела В0.
В случае, когда все тела представляют собой свободные тяготеющие точечные массы, имеем (см., например, [121J)
л л
(2.5) |
І7®> і/У. У
Asj = 1/ (*, — Xjf + (ys — у jf + (zs — Zjf;
T 1 _ ,, „ 2/3 „ 1 ms 2/3—4/3
"•з------- 2 * (/^o®«) і (Д jig) 3 (M>) s » £2 g)
os = sgn es = — 1.
Здесь / — постоянная тяготения, Д. і — расстояние между массами те и тj, ft, — полная энергия массы т, в ее эллиптическом кемеровском движении вокруг неподвижного притягивающего центрального тела, принимаемого также за точечную массу то, со,— частота обращения массы т„ в указанном движении (эту частоту в небесной механике называют средним движением и обозначают через п,), штрих при знаке суммы указывает на пропуск слагаемого, соответствующего j — s. Поскольку в данном случае о, = — 1, т. е. объекты являются мягко анизохронными, то из формулы (2.4) получаем
В = <(1/<и>)> - <(21(1,)> + В. (2.7)
Именно такая функция должна минимизироваться (по фазам а. — а,) на устойчивых синхронных движениях тел т,.
Если пренебречь величинами <(21(1>)> и В, что соответствует предположению о значительно большей величине массы центрального тела то по сравнению с массами Тої, ..., m*, а также
о малости диссипативных сил, то получится
(2.8) |
D^<{UlII')>,
Рис. 45. |
т. е. мы приходим к заключению о «приближенной минимальности» среднего потенциала взаимного тяготения масс в устойчивых синхронных движениях. Таким образом, эвристический принцип, предлагаемый в работе [324] (см. также [1061), действительно вытекает при соответствующих предположениях из ранее установленного в работах [40—42] интегрального критерия устойчивости синхронных движений.
Напомним, что для рассмотренной свободной орбитальной системы условия устойчивости, выраженные интегральным критерием, являются лишь необходимыми. Эти условия будут и достаточными, если вместо свободной рассмотреть соответствующую несвободную систему, т. е. систему, в которой массы т, движутся, например, по жестким эллиптическим каркасам (рис. 45). Конечно, такая несвободная система является достаточно грубой моделью планетных систем хотя бы потому, что уменьшение числа степеней свободы каждой массы т, до одной приводит к смене характера анизохронизма объекта.
Действительно, в случае точечной массы те, движущейся в поле тяготения неподвижной центральной массы т0 по заданной кривои, полная энергия
hs |
fm0ms
so
где Fs0 и rs0 — соответственно модуль скорости и расстояние меяеду массами тп, и тп0 в некоторый начальный момент времени. Поскольку теперь среднее движение со» не зависит от г6о, а определяется только начальной скоростью F„0, то
dh |
es = |
de>r |
(to. |
■m. V. |
(о. |
Нетрудно показать, что das/dVBQ> 0, т. e. что со* увеличивается при увеличении F«o. Поэтому теперь, в отличие от (2.6), es > 0 и о, = sgn ее = 1, т. е. объекты являются яїєстко анизохронными. Также жестко анизохронными являются твердые тела, вращающиеся вокруг некоторых фиксированных осей |ha = - j/eCOs,
ея =
В результате для несвободной (каркасной) системы на устойчивых синхронных движениях будет приближенно минимизироваться уже не средний потенциал сил тяготения <(UlIl))>, а средняя потенциальная энергия
<(П(П>)> = —<Ш(П)»,
т. е. средний потенциал <(Ї7П)> на таких движениях будет иметь (также приближенно) не минимум, а максимум.
Заметим далее, что устойчивости описанной несвободной орбитальной системы соответствует, по терминологии А. М. Ляпунова, так называемая условная устойчивость соответствующей свободной орбитальной системы [178].