Содержание к диссертации
Введение 3
Глава 1. Условия бесстолкновительного движения двойного астероида.
1. Уравнения движения и условия существования системы двойного асте
роида 25
Глава 2. Исследование устойчивости равновесий систем с односторонним ограничением по первым приближениям.
1 .Определение системы ОДУ с односторонним ограничением, её решений
и условий устойчивости неподвижных точек 51
Глава 3. Исследование устойчивости равновесий механических систем с односторонними голономными и дифференциальными связями по первому и второму приближению.
1.Общие соображения.-Модельная задача о космическом лифте как при
мер системы с (с,/) = 0 82
Иллюстративные примеры систем с односторонними дифференциальными связями 93
Задача о движении твёрдого тела с неподвижной точкой, на которое наложена односторонняя дифференциальная связь 99
Заключение 114
Список литературы 116
Введение к работе
Диссертация посвящена исследованию систем с ограничениями на значения фазовых переменных или с неудерживающими связями (в механических терминах). Это необходимо для изучения систем, состоящих из нескольких твёрдых тел, чьё взаимное расположение ограничено тем или иным образом. В работе рассматривается задача двух взаимно-гравитирующих твёрдых тел, и строится теория, позволяющая исследовать устойчивость состояний, возникающих из-за наличия односторонних связей в механике или ограничений общего вида для соответствующих образом определённых систем дифференциальных уравнений с неравенствами. Полученные теоретические результаты поясняются в иллюстративных примерах.
Для системы двойного астероида найдены достаточные условия на величину механической энергии и кинетического момента, обеспечивающие вечное бесстолкновительное движение тел на конечном расстоянии друг от друга.
Для систем обыкновенных дифференциальных уравнений общего вида с неравенством, ограничивающим допустимую для динамики область фазового пространства, получены условия (не)устойчивости неподвижных точек, возникающих из-за наличия неравенства. Эти и родственные результаты применяются для исследования механических систем с неудерживающими связями.
Небесномеханические системы из нескольких твёрдых тел активно исследуются в современной механике. Наиболее точный анализ возможен в случае систем вида "тело-точка", которые позволяют точно (численно) находить как возможные стационарные движения (D.J. Scheeres [14]), так и статистические характеристики типа плотности выпадения частиц на тело (пересечение орбиты точки и поверхности тела), как это было сделано в работе О.О. Василько-вой[37]. Такие задачи представляют интерес для исследования процесса формирования астероидов, но не всегда пригодны для изучения динамики конкретных систем.
Большинство известных автору современных работ, посвященных системам двойных астероидов, используют численное моделирование, опирающееся либо на экспериментальные данные (D.J. Scheeres и др.[9]), либо на авторские представления о типичности рассматриваемых тел. Тела могут приближаться многогранниками (D.J. Scheeres, E.G. Fahnestock [11]), либо иметь симметрии (W.S. Кооп и др. [13]). Эффективными оказываются модели, задаваемые телами относительно простой формы - эллипсоидами и сферами. Они появляются при желании аналитически исследовать динамику систем многих притягивающихся твёрдых тел (В.В.Белецкий, А.В. Родников [33], [4], В.В.Белецкий [23]), либо динамику в окрестности таких систем (К Gabern и др. [12], В.В. Белецкий [25]). Отдельно стоит упомянуть работы (В.В. Видякин [38]), где налагаемые на тела требования симметрии являются относительно слабыми, но находятся только решения, соответствующие классическим для случая точечных тел. Представления потенциала в таких ситуациях, полученные ещё Г.Н. Дубошиным [42], успешно применяются для численного моделирования (О.О. Василькова [37]). К сожалению, численное интегрирование уравнений не позволяет проводить качественный анализ систем общего вида.
В силу аналитической сложности получать общие результаты в полной задаче двух твёрдых тел произвольного вида представляется затруднительным (сложности имеют место уже для одного тела в поле ньютоновского центра [22] - предельный случай задачи двух твёрдых тел состоящей из очень массивного однородного шара и тела произвольной формы). Современные теоретические исследования, открытые книгой Г.Н. Дубошина [42] оперируют модельными системами тел простой формы. В общем случае до сих пор проводился лишь грубый качественный анализ. В полной постановке для изолированной системы двух твёрдых тел условия бесстолкновительности рассматривались в работах D.J. Scheeres[18], [19], в которых, однако, энергия и кинетический момент системы выступали как независимые параметры, которые приходилось подбирать. Ниже пороговые значения параметров будут явно вычисляться с использовани-
ем величины моментов инерции тел системы, что позволит сформулировать пригодные к практическому применению условия обеспечения сохранности системы, то есть устойчивости по Лагранжу некоторого класса её движений.
Допустим, что по тем или иным причинам нельзя ограничить взаимное расположение тел без постороннего вмешательства. Традиционным способом ограничения относительного движения является создание лёгких, но жёстких пространственных конструкций, удерживающих вместе части системы. Модельные задачи, открывающие эти исследования рассматривают базовую модель - две материальные точки, соединённые невесомым нерастяжимым стержнем и называемые гантелью, находящейся в гравитационном поле. Эти модели являются также дальнейшим упрощением, вводимым при описании твёрдых тел сложной формы (В.В. Белецкий [27], [32]). Даже классическая задача трёх тел (В.И. Арнольд и др.[20], А.П. Маркеев [51]) может быть обобщена для случаев подобного рода (WJ. Robinson [15], В.В. Белецкий [27]). Они допускают, например, исследование возмущений, вносимых атмосферой (В.В. Белецкий, М.Л. Пивоваров [31]), а также интересные инженерные применения (А.В. Родников [16], [17], М.А. Муницина [52]) для работы на околопланетной орбите.
Второй тип систем - с нежёстким закреплением частей. Характерной чертой в этом случае является малая масса удерживающих механизмов, имеющих при этом большую жёсткость на растяжение. Перспективные применения подобных конструкций рассматривались в (Белецкий В.В. [26], Белецкий В.В. Левин Е.М. [29]). Для их моделирования удобно применять теорию систем с односторонними связями и ударами (А.П. Иванов [43], В.В. Козлов, Д,В. Трещёв [46], В. Brogliato [6]). В частности, для них строятся методики исследования устойчивости. В книге [46] рассмотрен вопрос о реализации связей подобного рода, и конструкции, применяемые для их построения, эффективно используются для изучения модельных систем, таких, например, как биллиард Биркгофа. В книге [43] излагаются вопросы об импульсивном движении, возникающем при
наличии односторонних связей, приводятся авторские обобщения метода функций Ляпунова для исследования устойчивости положений равновесия и периодических движений систем с односторонними связями. В целом представляется, что существуют два основных способа определения решений - с помощью вариационных принципов ([46], В.В. Румянцев [60]) либо конструктивно [43]. Последний будет применяться ниже, так как он наиболее удобен для исследования конкретных систем. Существуют также два основных подхода для исследования устойчивости - прямой метод Ляпунова (использовался для исследования систем с импульсами у D. Bainov и др. в [1], [2], [3], разработан А.П. Ивановым для механических систем [43], [44]) и метод точечных отображений (СП. Горбиков [41]).
При формальном описании можно по-разному определять поведение системы, когда "включается" ограничение. Простейший и уже традиционный путь - зафиксировать моменты скачков, как это сделано в [54], [53] А.Д. Мышкисом, и в [58] Н.А. Перестюком и A.M. Самойленко. Он не подходит, однако, для механических систем с неудерживающими связями (см, напр. С.Н. Березинская и др. [34]). Введение дополнительной функции, которая будет использоваться для изменения векторного поля, определяющего систему при "включённой реакции" отличает предлагаемый ниже подход от используемого в [40],[39]. Выстраиваемая в работе теория допускает скачки фазовых переменных и работает для систем произвольной конечной размерности, что отличает работу от близкой по духу [45], где системы двумерны и при выходе на связь система к ней "прилипает". В отличие от [45] нерегулярные точки границы связи и неподвижные точки системы, попавшие на связь, не исследовались в работе, так как в случае общего положения они изолированы и такие ситуации пропадают при малой деформации границы связи. Единообразное описание позволяет отталкиваться от общей базы при описании систем, механическими аналогами которых являются системы с неудерживающими голономными и неголономными (в терминах [34]) связями. Полученные теоремы близки по смыслу к теоремам Ляпу-
нова об исследовании устойчивости положения равновесия по первому приближению (см. Н.Г. Четаев [65]). Для демонстрации взяты примеры, близкие к рассматриваемым в работах С.Н. Березинской и др. [34], И.И.Косенко и С.Я. Степанова [48], книги Ю.И. Неймарка и Н.А. Фуфаева [56].
Рассмотрим кратко основные моменты.
Астероид 243 Ида со спутником Дактиль
Последние исследования космического пространства показали наличие большого количества систем естественных космических тел, которые состоят из нескольких близко расположенных элементов - системы астероидов, состоящие из двух и более тел. Из последних результатов интересно открытие тройной системы Сильвия-Ром-Ремул [10].
Астероид Сильвия со спутниками Рем и Ромул (анимация NASA).
Наблюдения показывают, что элементы таких систем находятся на конечном расстоянии друг от друга и соударения в них, как правило, отсутствуют.
Полное рассмотрение подобных задач очень трудоёмко и выходит за рамки классической механики. Однако представляется возможным использовать подобные идеи в рукотворных конструкциях и исследовать соответствующие модели в рамках динамики космического полёта. Элементы системы при этом можно моделировать с помощью абсолютно твёрдых тел. Разделение системы на независимые модули упрощает её монтаж и ремонт, а также может способствовать повышению безопасности эксплуатации подобных систем, когда, например, источник энергии космической станции может представлять угрозу при аварии. В последнем случае его удаление сильно упрощается из-за отсутствия непосредственного контакта с остальными модулями станции.
В главе 1 для развития общих методов исследования подобных систем представлена задача по нахождению условий, при которых система, состоящая из двух произвольных твёрдых тел, не требует фиксирования частей друг относительно друга, но при этом её динамика не приводит к соударению тел или их разлёта на слишком большое расстояние. Это можно назвать устойчивостью по
Лагранжу движений, удовлетворяющих найденным условиям. Задача ставится в рамках небесной механики, но, разумеется, может применяться и в динамике космического полёта [57]. Ещё один аспект - это задача о невыходе на границу неудерживающеи голономнои связи, задаваемой как условие касания поверхностей тел.
В качестве модели будут рассматриваться взаимно-гравитирующие твёрдые тела Т1г Т2 с центрами масс Сх, С2, массами т^, т2 и центральными тензорами инерции J и А соответственно. Пусть (ег, е2, е3), (а, /?, у) - правые тройки, задающие вмороженные оси С&щ и C2afiy первого и второго тел. Внешние силы отсутствуют.
Можно показать, что динамика системы описывается уравнениями:
Г + СО X г
Ц(І7 + СО X V) =
ди дг
г , т dU , п- dU , dU , ди
/со + сох/со = ах — + рх — + ух — + гх—-
да г 9Р ду дг
П = СЛ-гСт(К0 -/со -цг xv) а = (ft - со) х а р = (ft - со) х р
у = (ft - со) х у С = [а р у]
(0.1)
в которых введены обозначения:
г=СгС2, v = І аг <»>, ft - абсолютные угловые скорости первого и второго тел,
m1m2
ц = приведенная масса системы.
mi+mz
U — /Т хТ . _ о oiCTi(Pi)t72(P2)^Pi^P2 " силовая функция гравитационного
dr взаимодействия тел, К0 = /со + Ail + \ir X — = k - кинетический момент сис-
темы относительно ее центра масс.
Все производные в (0.1) записаны в осях, вмороженных в первое тело. Тогда Л = СТАС - центральный тензор инерции второго тела в своих главных осях. Уравнения (0.1) замкнуты, но для завершения описания нужно ещё определить ориентацию первого тела в абсолютном пространстве. Это делается путём интегрирования уравнений Пуассона
Гйег
'—=(ОХЄі
dt г
5=wxe2 (0.2)
сіє*
после решения (0.1).
Уравнения (0.1) проще стандартных, получаемых из общих теорем динамики, если ставится задача поиска относительных равновесий. Теорема об изменении кинетического момента второго тела заменена более простым уравнением для вектора Ко. Это приводит к существенному усложнению соотношений, описывающих изменение ориентации второго тела относительно первого, и не позволяет говорить об упрощении при наличии динамики в связанных осях. Но выведенные уравнения имеют одно заметное преимущество: общие теоремы динамики, выписанные для каждого из тел, требуют одновременного отслеживания ориентации двух реперов, а не одного, как сделано выше. Усложнение в формулах Пуассона спасает от рассмотрения локальных производных в двух подвижных системах отсчёта, либо от нескольких соотношений с использованием непостоянных матриц при записи в абсолютных осях.
Г.Н. Дубошин в книге [42] получил уравнения движения системы многих твёрдых тел. Их скалярный вид удобен для применения методов усреднения или при наличии симметрии, тогда как рассматриваемые в данной главе больше подходят для качественных оценок из-за компактности их векторной записи. В уравнениях из [42] для случая двух твёрдых тел
также присутствуют выражения теорем об изменении кинетического момента каждого из тел.
В качестве демонстрации относительно несложно можно получить условия относительного равновесия системы - оба тела вращаются как твёрдое целое вокруг их совместной главной центральной оси инерции. Непосредственное нахождение таких решений затруднено общим видом силовой функции U.
Для грубого качественного описания динамики рассматриваемой системы предлагается найти условия, при которых тела вечно остаются на конечном расстоянии друг от друга и не сталкиваютсяГВ работе [18] была предложена методика исследования и получена оценка, существенно и независимо использующая модуль вектора К0. Ниже мы откажемся от его повсеместного использования и найдём правило вычисления порогового значения для величины кинетического момента системы, необходимой для существования системы двойного астероида.
Методика нахождения необходимых и достаточных условий, предложенная в [18] и развиваемая ниже, восходит к классической работе С. Смейла [61]. В силу сложного устройства полного потенциала взаимодействия тел предлагается получить приближённые оценки границ области возможности движения (ОВД).
Как известно, в случае наличия у механической системы интеграла (энергии) вида Т + V = h = const, где Т - положительно - определённая квадратичная функция скоростей, ОВД определяется условием V < h. Используя дополнительные первые интегралы, можно установить нижнюю границу значений функции Т, считая координаты заданными параметрами. Это позволяет уточнить, где именно могут быть движения системы, при её нахождении на совместном уровне интегралов.
Потенциал взаимодействия задан интегралом, что затрудняет аналити-чес кое исследование. Вводятся две оценки потенциальной энергии- сни-
зу И, < V и сверху V < V*. Выполнение неравенства V* < h влечёт V < h и гарантирует наличие движения в рассматриваемой точке конфигурационного многообразия. Неравенство V* > h напротив, обеспечивает V > h и обеспечивает невозможность движения с заданной энергией. Комбинируя эти условия можно оценивать положение границ ОВД.
Используя интеграл кинетической энергии для оценки снизу, предыдущие рассуждения в рассматриваемой задаче приводят к системе неравенств.
1 к^ 1
~ 6+1 - & 2(i*+62)2 ~ 5^1 - Х (0,3)
Здесь 5- безразмерное расстояние между центрами масс тел, к2 - безразмерный квадрат кинетического момента системы, і* - безразмерная сумма максимальных главных центральных моментов инерции тел, % - безразмерная механическая энергия системы. В них использованы оценки потенциальной энергии сверху и снизу соответственно. Получены следующие результаты:
Утверждение 1. Пусть значения параметров к2и і* таковы, что при некотором х > 1 выполнены условия
і) K2=*JpH
' х2-1
on 2^ а*+*2)3
2) KZ > —; —.
7 2хО;-1)2
Тогда существует промежуток значений безразмерной энергии х Є [Хі/Хг] < 0, при которых все движения системы, начавшиеся на расстояниях между центрами масс, достаточно близкими к х (или больших х), будут ограниченными и бесстолкновительными.
Утверждение 2. Пусть значения параметров к2 и і* таковы, что при некотором х > 1 выполнены условия:
1)к2=Нії!>!
' х-1
2) к > 0^ .2.
,2 ^ (і*+*2)3
гх(х-\)2
Тогда существует промежуток значений безразмерной энергии [#*, 0) < 0, при которых движения системы, начавшиеся на расстояниях не меньших х, будут ограниченными и бесстолкновительными.
Утверждение 1 может иметь место только при значении і* безразмерной суммы максимальных из главных центральных моментов инерции тел из
интервала 0, Л0 J. Утверждение 2 может работать при любых значе
ниях *, но оно требует больших величин кинетического момента систе
мы для обеспечения бесстолкновительности, нежели утверждение 1.
