Эволюция движения механических систем с бесконечным числом степеней свободы
На правах рукописи
Шатина Альбина Викторовна ЭВОЛЮЦИЯ ДВИЖЕНИЯ МЕХАНИЧЕСКИХ СИСТЕМ С БЕСКОНЕЧНЫМ ЧИСЛОМ СТЕПЕНЕЙ СВОБОДЫ Специальность: 01.02.01 – теоретическая механика
АВТОРЕФЕРАТ
диссертации на соискание ученой степени доктора физико-математических наук
Москва – 2007
Работа выполнена на кафедре Высшей математики Московского государственного института радиотехники, электроники и автоматики (технического университета)
Научный консультант:
доктор физико-математических наук, профессор Вильке Владимир Георгиевич.
Официальные оппоненты:
член-корреспондент РАН, доктор физико-математических наук, профессор Белецкий Владимир Васильевич;
доктор физико-математических наук, профессор Баркин Юрий Владимирович;
доктор физико-математических наук, профессор Марков Юрий Георгиевич.
Ведущая организация:
Вычислительный центр им. А.А. Дородницына Российской академии наук.
Защита диссертации состоится 18 мая 2007 года в 16 часов на заседании диссертационного совета Д 501.001.22 при Московском государственном университете им. М.В. Ломоносова по адресу: 119992, Москва, Ленинские горы, МГУ, механико-математический факультет, аудитория 16-10.
С диссертацией можно ознакомиться в библиотеке механико математического факультета МГУ (Главное здание, 14 этаж).
Автореферат разослан2007 года
Ученый секретарь диссертационного совета кандидат физико-математических наук, доцент В.А. Прошкин
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы. Диссертация посвящена исследованию эволюции поступательно-вращательного движения механических систем с бесконечным числом степеней свободы, моделируемых либо вязкоупругими телами, либо твердыми телами с жестко прикрепленными к ним вязкоупругими элементами.
Вопрос о влиянии внутреннего вязкого трения на поступательно вращательное движение деформируемого тела возник много лет назад, прежде всего в связи с изучением приливной эволюции движения планет Солнечной системы. В небесной механике для описания движений естественных и искусственных тел, как правило, используются простейшие модели классической механики – материальная точка и абсолютно твердое тело, а приливная теория базируется на ряде гипотез относительно величины приливных горбов и их расположения относительно вращающейся планеты.
Первые фундаментальные работы по изучению приливной эволюции в системе “планета-спутник” были выполнены в конце XIX века небесным механиком и космогонистом Дж.Г. Дарвиным. В 60-е годы XX века новая научная информация о планетах и спутниках, полученная с помощью космических аппаратов и радиолокационной астрономии возродила интерес к приливной теории. Макдональд Г., Гойлдрах П., Пил С., Каула У. и ряд других авторов провели более детальное исследование эффектов приливного трения.
Приливные силы оказывают влияние не только на орбиты спутников, но также и на их вращательные движения, и являются важнейшим диссипативным фактором, приводящим произвольное первоначальное движение тела к захвату в резонансный режим движения.
Ряд важных результатов по приливной эволюции вращательного движения небесных тел был получен Белецким В.В. (1975, 1978).
Исследования по влиянию деформируемости и внутренней вязкости материала на движение тела стали также актуальны в связи с попыткой объяснения расхождений между теоретическим результатами и данными наблюдений в динамике Земли, с появлением искусственных спутников и обнаружением новых эффектов в их угловых движениях, обусловленных упругими свойствами.
При изучении динамики систем, содержащих вязкоупругие элементы большой жесткости, широко применяется асимптотический метод, предложенный Черноусько Ф.Л. (1978, 1980). Получаемые с помощью этого метода уравнения движения имеют вид уравнений динамики твердого тела с дополнительными слагаемыми, обусловленными внутренней упругостью и диссипацией, и описывают движение системы, устанавливающееся после затухания собственных упругих колебаний и вызванное внешними силами и силами инерции.
Для исследования эволюции движения механических систем, содержащих вязкоупругие элементы большой жесткости, используется асимптотический метод разделения движений и усреднения (Вильке В.Г., 1983), который позволяет перейти от уравнений в бесконечномерных банаховых пространствах к конечной системе обыкновенных дифференциальных уравнений, неизвестными в которой являются переменные “действие” невозмущенной задачи.
Одно из важных направлений, посвященных изучению динамики распределенных систем, представлено в монографиях Вильке В.Г. (1986, 1997), где хорошо развитые методы классической аналитической механики обобщаются на случай систем с бесконечным числом степеней свободы. Для ряда моделей рассматриваемых систем доказаны теоремы существования и единственности обобщенных решений уравнений движения. Рассмотрены задачи о движении вязкоупругого тела в центральном ньютоновском поле сил, системы вязкоупругих тел, взаимодействующих по закону всемирного тяготения. Эти и другие задачи являются модельными задачами в теории приливов на новом ее этапе развития, когда на смену моделей абсолютно твердого тела и материальной точки приходит модель пространственного вязкоупругого тела.
Цель диссертации состоит в развитии и углублении методов исследования эволюции движения систем с бесконечным числом степеней свободы и применении этих и ранее известных методов к исследованию диссипативной эволюции движения естественных и искусственных тел в задачах небесной механики и механики космического полета.
Основные результаты диссертации и их научная новизна.
• Для исследования эволюции поступательно-вращательного движения вязкоупругого тела, а также твердого тела с жестко прикрепленными к нему вязкоупругими элементами, предложен асимптотический метод, сочетающий в себе метод разделения движений и метод Крылова-Боголюбова для систем с быстрыми и медленными переменными. Этот метод в предположении больших коэффициентов жесткости и демпфирования вязкоупругой среды, а также наличия малых внешних периодических возмущений, позволяет получить приближенную систему обыкновенных дифференциальных уравнений, описывающую медленную диссипативную эволюцию поступательно-вращательного движения изучаемой механической системы.
• Получена эволюционная система уравнений и проведен ее анализ в задаче о поступательно-вращательном движении вязкоупругого шара в центральном ньютоновском поле сил. Указанная система уравнений описывает взаимное изменение параметров орбиты и вращательного движения вязкоупругого шара и может быть использована при изучении приливной эволюции движения планет Солнечной системы.
• Исследована эволюция поступательно-вращательного движения вязкоупругого шара в ограниченной круговой задаче трех тел, являющейся модельной задачей небесной механики. Рассмотрен случай, когда массивные тела, моделируемые материальными точками с массами 1 и µ ( µ 1), и центр масс вязкоупругого шара движутся в одной плоскости, а вращение вязкоупругого шара относительно центра масс происходит вокруг нормали к этой плоскости. Показано, что эволюция движения вязкоупругого шара разбивается на три этапа, характеризующихся различными временами. Получены уравнения, описывающие движение вязкоупругого шара на каждом из этих этапов в классе квазикруговых орбит.
• Исследована эволюция движения двойной планеты, моделируемой материальной точкой массы и однородным деформируемым вязкоупругим шаром массы m, в гравитационном поле неподвижного центра – материальной точки единичной массы µ ( µ m 1) в классе квазикруговых орбит. Рассмотрен случай, когда движение центра масс вязкоупругого шара и материальной точки массы µ происходит в одной плоскости, проходящей через притягивающий центр, а вращение шара относительно центра масс происходит вокруг нормали к этой плоскости. Построены фазовые портреты, позволяющие выявить взаимосвязь между изменением радиуса орбиты центра масс двойной планеты и расстоянием между составляющими ее элементами.
• Получена эволюционная система уравнений, описывающая поступательно-вращательное движение сферически симметричного спутника с гибкими вязкоупругими стержнями в центральном ньютоновском поле сил. Найдены многообразия стационарных решений и исследована их устойчивость на основе уравнений в вариациях.
• Исследована эволюция вращательного движения динамически симметричного спутника с гибкими вязкоупругими стержнями на круговой орбите. Получены приближенные уравнения, описывающие поведение рассматриваемой механической системы на этапах “быстрой” и “медленной” диссипативной эволюции. Для этапа “медленной” диссипативной эволюции построены фазовые портреты.
• Рассмотрен пространственный вариант задачи о двойной планете.
Методом разделения движения построена возмущенная система уравнений, описывающая поступательно-вращательное движение планеты и спутника с учетом возмущений, вызываемых упругостью и диссипацией. Проведен анализ квазистатических деформаций вязкоупругой планеты, моделируемой однородным изотропным вязкоупругим шаром. Получено квазистатическое решение задачи теории упругости в случае, когда планета моделируется механической системой, состоящей из абсолютно твердой невесомой сферы, к которой с внешней стороны жестко прикреплена вязкоупругая сферическая оболочка, а внутри имеется подвижное внутреннее ядро.
Все результаты диссертации являются новыми.
Методы исследования. В работе используются методы аналитической механики, метод разделения движений, применяемый к механическим системам, содержащим деформируемые элементы большой жесткости (Черноусько Ф.Л. (1980), Вильке В.Г. (1983)), метод усреднения, асимптотический метод Крылова-Боголюбова для систем с быстрыми и медленными переменными (Арнольд В.И., Козлов В.В., Нейштадт А.И. (1985)).
Достоверность результатов. Все результаты диссертационной работы получены на основе сформулированных в ней гипотез и строго обоснованы.
Практическая ценность работы. Результаты диссертационной работы, полученные в задачах о поступательно-врашательном движении вязкоупругого шара в центральном ньютоновском поле сил, в задаче о двойной планете и в ограниченной круговой задаче трех тел, могут найти применение в приливной теории движения планет и их спутников и геодинамике. Результаты исследования поступательно-вращательного движения спутника с гибкими вязкоупругими стержнями в центральном ньютоновском поле сил и вращательного движения спутника на круговой орбите могут быть использованы на этапе предварительного проектирования больших космических конструкций. Методические приемы, использованные при выводе точных уравнений движения в рамках рассматриваемой модели и получении приближенной системы уравнений, описывающей поступательно-вращательное движение, могут быть полезны при исследовании сложных механических систем.
Апробация работы. Результаты, изложенные в диссертации, докладывались на III-ем и IV-ом Международных симпозиумах по классической и небесной механики (Великие Луки, 1998, 2001), на Международном научном симпозиуме по проблемам механики деформируемых тел, посвященном 90-летию со дня рождения А.А.
Ильюшина (Москва, 2001), на Восьмом Всероссийском съезде по теоретической и прикладной механике (Пермь, 2001), на VIII-ой Международной конференции Устойчивость, управление и динамика твердого тела (Донецк, 2002). Результаты диссертации докладывались на научных семинарах механико-математического факультета МГУ им.
Ломоносова под руководством акад. Румянцева В.В. и проф. Карапетяна А.В. (1999, 2002, 2007), под руководством чл.-корр. РАН Белецкого В.В. и проф. Голубева Ю.Ф. (2001, 2002), под руководством проф. Вильке В.Г. и проф. Самсонова В.А. (2001);
на научном семинаре в ГАИШ им.
Штернберга П.К. под руководством проф. Баркина Ю.В. (2001), на научном семинаре кафедры теоретической механики Московского энергетического института под руководством проф. Мартыненко Ю.Г.
(2001), на 52-ой, 53-ей, 54-ой научно-технических конференциях МИРЭА (Москва, 2003, 2004, 2005), на семинаре под руководством проф. Шевалье Л. и Паскаль М. в университете Пьера и Мари Кюри (Париж, Франция, 2001).
Публикации. Основные результаты диссертации изложены в статьях (из них 10 – в рецензируемых отечественных журналах), список которых приведен в конце автореферата. В совместных работах с Вильке В.Г. авторы внесли равный вклад и несут равную ответственность за полученные результаты.
Структура и объем диссертации. Диссертация состоит из введения, семи глав, заключения и списка литературы. Работа содержит 34 рисунка. Общий объем диссертации составляет 250 страниц. Список цитируемой литературы включает 141 наименование.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность темы диссертации, сформулирована цель работы, приведен обзор литературы соответствующей тематики и кратко изложено содержание диссертации.
Глава 1. Асимптотические методы в механике систем с бесконечным числом степеней свободы.
В первой главе диссертации формулируется постановка задачи о поступательно-вращательном движении механической системы с бесконечным числом степеней свободы, выводятся уравнения движения и излагаются асимптотические методы исследования. Под механической системой с бесконечным числом степеней свободы понимается либо вязкоупругое тело, либо твердое тело с жестко прикрепленными к нему вязкоупругими элементами.
Постановка задачи и вывод уравнений движения в форме уравнений Рауса из вариационного принципа Д’Аламбера-Лагранжа составляют содержание §1.1.
Предполагается, что рассматриваемая механическая система занимает область = 1 2 в евклидовом пространстве 3, причем 1 это область, занимаемая твердым телом, а 2 вязкоупругим в естественном недеформированном состоянии (в случае вязкоупругого тела 1 = ). Помимо инерциальной системы координат OXYZ вводятся система осей Кенига C1 2 3 и подвижная система координат Cx1 x 2 x 3 с началом в центре масс механической системы в деформированном состоянии. Положение точки M относительно инерциальной системы координат OXYZ определяется радиус-вектором R (r, t ) = R C (t ) + (t )(r + u(r, t )), (1) где r, u(r, t ) – вектор упругого смещения, обращающийся в нуль для радиус-вектор точки C, – оператор точек твердого тела, R C – перехода от системы координат Cx1 x 2 x 3 к системе осей Кёнига C1 2 3.
Радиус-вектор R C и “связанная” с деформируемой системой система координат Cx1 x 2 x 3 определяются однозначно по заданному векторному полю R (r, t ) следующими условиями:
u(r, t ) (r)dx = 0, rot u (r)dx = 0, R C (t ) = R (r, t ) (r )dx, (2) m 2 m где масса всей системы, распределенная плотность, dx = dx1 dx 2 dx 3.
Предполагается, что деформируемая часть рассматриваемой механической системы однородна, изотропна и имеет постоянную плотность, а относительные перемещения точек при деформациях u(r, t ) малы.
Потенциальная энергия упругих деформаций задается в соответствии с линейной теорией упругости квадратичным функционалом E [u ], который по предположению содержит множителем “большой” параметр N, характеризующий жесткость упругой среды.
Вводится малый параметр = N 1. При = 0 вектор-функция u(r, t ) полагается равной нулю, и механическая система представляет собой абсолютно твердое тело. Функционал внутренних диссипативных сил D[u] = E [u] D[u] определяется соотношением где &, – коэффициент внутреннего вязкого трения, т.е. рассматривается модель Кельвина-Фойгта.
Кинетическая энергия системы определяется равенством:
[ ] 1 &2 1 1 R C + (r + u ) + u dx, () = 1(), & & T= R dx = & (3) 2 а функционал потенциальной энергии внешних сил полагается заданным в виде:
= [R,, µ ], [R, + 2, µ ] = [R,, µ ], (4) где (t ) = 0 t + (0), 0 = const, µ малый параметр.
Вариационный принцип Д’Аламбера-Лагранжа имеет вид:
RRdx + ( ) + ( uE [u] + uD[u], u ) && R, R + & & rot udx = 0, + 1 udx + 2 u B, (5) 2 где 1 = 1 (t ), 2 = 2 (t ) – неопределенные множители Лагранжа, порожденные условиями (2).
Конфигурационным пространством системы является прямое произведение M B, где M – шестимерное дифференцируемое многообразие с локальными координатами q1,..., q 6, а B – банахово пространство векторных функций u(r, t ). Обобщенные координаты q1,..., q 6 определяют положение точки C в системе координат OXYZ и ориентацию подвижной системы координат Cx1 x 2 x 3 относительно осей Кенига C1 2 3 :
R C = R C (q1, q 2, q 3 ), = (q 4, q5, q 6 ), откуда следует, что T = T [q, q, u, u ], = [q, u,, µ ].
& & Уравнения движения выводятся из вариационного принципа Д’Аламбера-Лагранжа (5) и выписываются в форме уравнений Рауса, причем канонические переменные (p, q ), где p = q T – обобщенные & импульсы, используются для описания поступательно-вращательного движения деформированной системы, а лагранжевы переменные u i (r, t ) (i = 1,2,3) для описания деформаций.
Функционал Рауса определяется равенством R = pi qi T + + E [u] & (q, q ) (p, q ), & i = а уравнения движения имеют вид:
R R pk =, qk =, k = 1,...,6, & & (6) q k p k d uR + uR + uD[u ] + 1, u + 2 rot udx = 0 u B & (7) & & dt 2 Система уравнений (6)–(7) представляет собой сложную систему интегродифференциальных уравнений в банаховом пространстве, содержащую малые параметры и µ. Малый параметр отвечает за возмущения, связанные с нежесткостью рассматриваемой механической системы, а малый параметр µ – с малыми возмущениями, индуцируемыми полем внешних сил. Предполагается, что при = 0, µ = 0 задача является интегрируемой, а потому в целях использования асимптотических методов от переменных (p, q ) осуществляется переход к переменным (I, ), которые в невозмущенной задаче (т.е. при = 0, µ = 0 ) являются переменными “действие-угол”. В классе рассматриваемых в данной диссертации задач в роли таких переменных выступают переменные Андуайе-Делоне.
В §1.2 дано описание переменных Андуайе-Делоне, выписаны соотношения, позволяющие выразить функционал Рауса в этих переменных, т.е. представить его в виде R = R[I,, u, u,, µ ]. Показано, & что уравнения (6) в переменных (I, ) сохраняют свою структуру, т.е.
имеют вид:
R R & Ik =, k =, k = 1,..., 6.
& (8) k I k и совместно с уравнениями (7) и условиями (2) образуют полную систему уравнений, определяющую движение данной механической системы.
На следующем шаге к полученной системе уравнений применяется метод разделения движений, описание которого содержится в §1.3.
При = 0, когда жесткость упругой среды бесконечно велика, перемещения u(r, t ) обращаются в нуль, и рассматриваемая механическая система движется как твердое тело. Функционал Рауса в этом случае имеет вид:
R[I,, 0,0, µ ] = R0 (I ) + µK [I,,, µ ], а уравнения (0.8) выглядят следующим образом:
K K R & Ik = µ, k = k ( I ) + µ, k = 1,...,6, k (I ) = &. (9) k I k I k Уравнение (7) после умножения обеих его частей на N становится сингулярно возмущенным, т.е. содержащим малый параметр при старшей производной по времени. Из уравнения (7) определяется вектор-функция u, описывающая вынужденные колебания точек вязкоупругой среды, обусловленные действием внешних сил и сил инерции, которые называются квазистатическими. Указанное решение представляется в виде:
u = u1 (r, I,,,, µ ) + O( ( ) 2 ), (10) где = max k (I (0)).
k & После линеаризации правых частей уравнений (8) по u и u и & подстановки в них вместо u и u найденного частного решения (10) получается замкнутая система обыкновенных дифференциальных уравнений относительно переменных, описывающих поступательно вращательное движение механической системы, и учитывающая возмущения, связанные с внутренней упругостью и диссипацией.
Математическое обоснование и оценка погрешности метода разделения движений в динамике систем с упругими и диссипативными элементами в случае, когда эта система имеет конечное число степеней свободы, изложено в работах Черноусько Ф.Л., Шамаева А.С. (1983), Синицына Е.В. (1991) на основе метода пограничных функций для систем сингулярно возмущенных обыкновенных дифференциальных уравнений (Васильева А.Б., Бутузов В.Ф. (1973)), а в работе Сидоренко В.В. (1995) на основе метода теории интегральных многообразий (Митропольский Ю.А., Лыкова О.Б. (1973), Стрыгин В.В., Соболев В.А. (1988)). В работе Шатиной А.В. (1990) получено асимптотическое представление решения линейного дифференциального уравнения второго порядка в банаховом пространстве с малым параметром при старшей производной и неоднородном члене.
В §1.4 излагаются асимптотические методы исследования полученной в §1.3 системы обыкновенных дифференциальных уравнений, которая имеет вид:
K + Fk (I,,,, µ ), k = 1,..., n, & I k = µ k K l = l (I ) + µ + Gl (I,,,, µ ), l = 1,..., m.
& (11) I l В рассматриваемой постановке задачи k, l = 1,..., 6. Однако, вследствие применения метода разделения движений правые части получаются независящими от некоторых угловых переменных и, следовательно, уравнения для этих переменных отделяются от системы (11). Порядок системы уравнений (11) в ряде случаев также понижается из-за наличия первых интегралов.
При = 0 система уравнений (11) становится стандартной в смысле применения к ней метода усреднения (Арнольд В.И. (1978)). Переменные “действие” в этом случае не эволюционируют, т.к. при отсутствии диссипации энергии, вызываемой внутренним вязким трением упругой среды, данная система является консервативной.
Если K 0, то система уравнений (11) исследуется методом усреднения. Усреднение по быстрым угловым переменным в предположении отсутствия резонансов приводит к замкнутой системе обыкновенных дифференциальных уравнений относительно переменных действие и медленных угловых переменных.
Если K (I,,, µ ) 0, то система уравнений (11) рассматривается как система с малым параметром µ при фиксированном значении параметра. Рассматривается ситуация, когда при µ = 0 уравнения (11) выглядят следующим образом:
& I = F (I ), k = 1,..., n, (12) k k l = l (I ) + Gl 0 (I ), l = 1,..., m.
& (13) Система уравнений образует замкнутую систему (12) обыкновенных дифференциальных уравнений относительно переменных действие. Предполагается, что она обладает притягивающим множеством аттрактором, к которому стремятся почти все решения этой системы.
Это притягивающее множество содержится во множестве стационарных решений, определяемых системой уравнений:
Fk 0 (I ) = 0, k = 1,..., n.
Решение системы уравнений (11) по аналогии с методом методом Крылова-Боголюбова для систем с быстрыми и медленными переменными (Арнольд В.И., Козлов В.В., Нейштадт А.И. (1985)) ищется в виде:
I k = J k + µN k1 (J,,, ) + µ 2 N k 2 (J,,, ) + L, k = 1,..., n, l = l + µM l1 (J,,, ) + µ 2 M l 2 (J,,, ) +L, l =1,..., m, (14) где J = ( J 1,..., J n ), = (, 1,..., m ), функции N ks, M lp 2 -периодичны по угловым переменным 1,..., m, и имеют по этим переменным нулевое среднее.
1,..., m Как функции времени переменные и J 1,..., J n определяются следующей системой дифференциальных уравнений:
& J = A (J,, ) + µA (J,, ) + µ 2 A (J,, ) +L, k =1,..., n, k k0 k1 k l = Bl 0 (J, ) + µBl1 (J,, ) + µ 2 Bl 2 (J,, ) +L, l =1,..., m, & (15) где Bl 0 (J, ) = l (J )+ Gl 0 (J ).
Ak 0 (J,, ) = Fk 0 (J ), (16) µ Полагая в асимптотически устойчивого -окрестности ~ Ak 0 (J,, ) = µAk 0 (J ) и стационарного решения системы (12) подставляя разложения (14), (15) в уравнения (11), приравниванием коэффициентов при одинаковых степенях малого параметра µ можно последовательно определить неизвестные функции Akj, N kj, Blj, M lj (k = 1,..., n, l = 1,..., m, j = 1,2,...) при условии отсутствия резонансных соотношений.
Если Ak1 (J,, ) 0, k = 1,..., n, то приближенные уравнения, описывающие эволюцию переменных действие на этапе медленной диссипативной эволюции, имеют вид:
& J = F (J ) + µF (J ), k =1,..., n. (17) k k0 k Если же Ak1 (J,, ) 0, k = 1,..., n, то следует найти функции Ak 2 (J,, ), определяющие эволюцию переменных действие во втором приближении по малому параметру µ. Функции Ak 2 (J,, ) второго приближения по малому параметру µ необходимо искать также в случае, когда система уравнений (17) обладает асимптотически устойчивым стационарным решением. В µ 2 - окрестности этого стационарного решения необходимо учитывать следующее приближение по µ. Тогда эволюционная система уравнений примет вид:
& J = ( F (J ) + µF (J ) + µ 2 F (J )), k =1,..., n. (18) k k0 k1 k Система уравнений (12) описывает “быструю” диссипативную эволюцию переменных действие, а приближенная система дифференциальных уравнений (17) (или (18)) описывает “медленную” диссипативную эволюцию этих переменных.
Глава 2. Эволюция движения вязкоупругого шара в центральном ньютоновском поле сил.
Вторая глава диссертации посвящена исследованию эволюции движения вязкоупругого шара в центральном ньютоновском поле сил.
В §2.1 формулируется постановка задачи, выводятся уравнения движения в форме уравнений Рауса с использованием переменных Андуайе-Делоне, составляющих гамильтонову часть переменных и описывающих поступательно-вращательное движение шара.
Обобщенные координаты, описывающие деформацию шара, составляют лагранжеву часть переменных. В данной задаче поле внешних сил центральное, и, следовательно, сохраняется вектор момента количеств движения G 0 относительно притягивающего центра, который равен сумме векторов моментов количеств движения в орбитальном движении и в движении относительно центра масс: G 0 = G + K. Инерциальная система координат OXYZ выбирается так, что ее начало совпадает с притягивающим центром, а ось OZ направлена по вектору G 0.
Деформированное состояние шара описывается классической теорией упругости малых деформаций, в частности функционал потенциальной энергии упругих деформаций имеет вид:
E [u ] = 1 ( I 2 2 II )dx, 1 0, 0 2 3, (1 ) 2(1 2 ) 1 =, 2 =, 2(1 + )(1 2 ) 1 u i u j 3 I = eii, II = (eii e jj eij ), eij = + 2 x j x i i =1 i j где – модуль упругости Юнга, – коэффициент Пуассона.
Потенциальная энергия внешнего и внутреннего гравитационных полей представляется функционалами fm(r, u ) dx = 1 = dx,, (R С + (r + u ) ) r где = fm1, f – гравитационная постоянная, m1 – масса притягивающего центра, m - масса вязкоупругого шара, r0 – его радиус в естественном недеформированном состоянии.
В §2.2 методом разделения движений и усреднения осуществляется построение приближенной системы обыкновенных дифференциальных уравнений относительно переменных “действие”, описывающих диссипативную эволюцию поступательно-вращательного движения вязкоупругого шара в центральном поле сил.
В соответствии с рассматриваемой в первой главе моделью жесткость деформируемого шара предполагается большой, и вводится малый параметр, обратно пропорциональный модулю упругости Юнга.
При = 0 соответствующая система уравнений описывает движение абсолютно твердого шара в центральном ньютоновском поле сил. Эта система является интегрируемой, а движение, ею описываемое, таково:
центр масс шара движется по кеплеровской орбите как материальная точка с массой, равной массе шара;
при этом шар равномерно вращается вокруг оси, неизменно ориентированной в инерциальной системе координат. Рассматриваются финитные движения, когда кеплеровская орбита в невозмущенном движении представляет собой эллипс. Это невозмущенное движение используется в качестве порождающего для определения вынужденных колебаний вязкоупругого шара.
Решение задачи о деформации шара при отсутствии массовых сил было получено Томсоном с использованием сферических функций.
Специальный вид массовых сил в исследуемой задаче позволяет получить ее решение в виде суммы однородных сферических функций третьего и первого порядков (Ляв А. (1935), Лейбензон Л.С. (1942)). Функции, описывающие вынужденные колебания вязкоупругого шара, используются для формирования правых частей “возмущенной” системы уравнений, описывающей его поступательно-вращательное движение.
Далее применяется метод усреднения (рассматривается нерезонансный случай).
Полученные в результате указанной процедуры эволюционные уравнения представляют собой замкнутую систему обыкновенных дифференциальных уравнений шестого порядка относительно переменных действие:
( ) ( ) 1 + z 2 3 9 & I 1 = 4 1 e 2 cos 2 F1 (e) 2 zF2 (e) 4 1 e, 2 ( ) ( ) 9 2 1 + 3 z & I 2 = 4 1 e 2 F1 (e) 2 zF2 (e) 4 1 e, 2 ( ) ( ) 9 2 1 3 & I 3 = 4 1 e 2 ( z cos i + cos 1 ) F1 (e) 2 cos iF2 (e) 4 1 e 4, 2 ( ) ( ) 6 3 L = 4 1 e 2 zF2 (e) 2 F3 (e) 4 1 e &, ( ) ( ) 9 2 3 & G = 4 1 e 2 zF1 (e) 2 F2 (e) 4 1 e 4, & & H = I 3, (19) где 4r07 (1 + )(9 + 13) 2m3 I I = 18 D 2, D 2 =, 4 =, 2 = 2, cos 1 = 3, 105(5 + 7) L3 A I I cos 2 = 1, cos i = H / G, z = cos( 1 + i ), e 2 = 1 G 2 L2, I 3 15 45 F1 (e) = 1 + 3e 2 + e 4, F2 (e) = 1 + e 2 + e 4 + e 6, 8 2 8 31 255 4 185 6 25 F3 (e) = 1 + e 2 + e+ e + e.
2 8 16 Здесь I 2 = K, K – вектор момента количеств движения шара относительно центра масс, I 1, I 3 - проекции вектора K на оси Cx 3 и C соответственно, L, G, H - переменные Делоне, описывающие орбиту центра масс. Переменная G – модуль вектора момента количеств движения G в орбитальном движении, H – проекция вектора G на ось OZ, e – эксцентриситет орбиты центра масс шара, i – наклонение орбиты, cos i = H / G.
Уравнения (19) имеют три первых интеграла. Два из них являются следствием закона сохранения момента количеств движения G 0 ;
третий является следствием сферической симметрии и соответствует тому факту, что в усредненных уравнениях угол между осью Cx 3 связанной с шаром подвижной системы координат, и вектором кинетического момента K относительно точки С остается неизменным. Эти уравнения описывают приливную эволюцию поступательно-вращательного движения вязкоупругой планеты в поле притягивающего центра.
В §2.3 проводится анализ эволюционной системы уравнений.
Находятся стационарные решения системы и исследуется их устойчивость на основе уравнений в вариациях. Показано, что в стационарном движении центр масс движется по круговой орбите, ортогональной вектору G 0, ось вращения шара направлена по нормали к плоскости орбиты, при этом деформированный шар обращен одной стороной к притягивающему центру. Если модуль вектора G удовлетворяет условию G 0 4 3 A 2 m 3, (20) где A – момент инерции недеформированного шара относительно его диаметра, – гравитационная постоянная, m – масса шара, то имеются две стационарные орбиты. Если G 0 = 3 A 2 m 3, то такая орбита одна.
В случае, когда G 0 3 A 2 m 3, стационарных решений эволюционная система не имеет. В случае двух стационарных орбит устойчивым является стационарное движение, соответствующее орбите большего радиуса, а стационарное движение с меньшим радиусом орбиты неустойчиво. Этот результат соответствует полученному ранее в работе Вильке В.Г. (1980) при анализе векторных эволюционных уравнений относительно вектора кинетического момента вращательного движения K и радиус-вектора центра масс шара.
Следует отметить, что в случае существования одной стационарной орбиты ее радиус равен 1,2r0, где r0 радиус шара, а в случае двух стационарных орбит радиус меньшей орбиты меньше этого значения.
Однако эволюционные уравнения выводятся в предположении, что радиус шара много меньше расстояния от его центра масс до притягивающего центра. Это предположение исключает из рассмотрения орбиты с радиусами порядка r0. В изолированной системе планета – притягивающий центр (Солнце) условие (20) выполнено для всех планет Солнечной системы. Таким образом, можно утверждать, что диссипативная эволюция поступательно-вращательного движения планеты в центральном ньютоновском поле сил протекает так, что ее движение стремится к стационарному.
Тот факт, что приливные силы приводят первоначальное вращение тела к захвату в резонансный режим движения, когда период обращения по орбите и период вращения вокруг оси совпадают, был установлен Белецким В.В. (1965) при исследовании вращательного движения твердого тела под действием приливных сил, момент которых определялся феноменологически. В Солнечной системе почти стационарное движение осуществляет Луна относительно Земли, спутники Юпитера Ио, Европа, Ганимед, Каллисто, ряд других спутников планет.
В силу полученной эволюционной системы уравнений найдены производные по времени от наклонения орбиты i (угла между векторами G 0 и G ), угла 1 между векторами G 0 и K, и эксцентриситета орбиты e. Показано, что существует класс движений, когда вращение вязкоупругого шара относительно центра масс происходит вокруг нормали к плоскости орбиты, совпадающей с вектором G 0. Если i 0, то di 0, и, следовательно, наклонение орбиты в процессе эволюции dt уменьшается. Эволюция угла 1 имеет более сложный характер. В частности, в случае обратных вращений (когда 2 1 + i ) угол уменьшается, а в случае прямых вращений ( 0 1 + i 2 ) может происходить увеличение угла 1, если угловая скорость вращения планеты вокруг своей оси превышает ее среднее движение по орбите.
Эффект возможного увеличения угла 1 с последующим его уменьшением до нуля описан в работах Белецкого В.В. (1975) при исследовании приливной эволюции вращательного движения планет.
Показано, что уравнения движения допускают класс орбит с нулевым эксцентриситетом. Если e 0, то в случае обратных вращений e уменьшается, а в случае прямых вращений может происходить увеличение эксцентриситета при достаточно большой угловой скорости собственного вращения по сравнению с орбитальной.
В §2.4 рассмотрен частный случай задачи о движении вязкоупругого шара в центральном ньютоновском поле сил, когда центр масс шара движется в плоскости OXY, а вращение относительно центра масс происходит вокруг нормали к плоскости орбиты. Правые части “возмущенной” системы уравнений в этом случае зависят лишь от одной угловой переменной – истинной аномалии, по которой и производится усреднение. Проблема резонансов в данном случае отсутствует.
Построены фазовые портреты в плоскости переменных ( L, I ), где I =| K | в случае прямого вращения и I = | K | в случае обратного вращения, а L =| G | (1 e 2 ) 0.5, а также в плоскости переменных ( I, e).
Третья глава диссертации посвящена исследованию эволюции движения вязкоупругого шара в ограниченной круговой задаче трех тел.
Рассматривается следующая постановка задачи. Два массивных тела, моделируемых материальными точками M 1 и M 2 с массами равными единице и µ, где µ 1, движутся по круговым орбитам под действием сил ньютоновского притяжения вокруг общего центра масс O в плоскости OXY. При этом OM 2 = b, OM 1 = µb, а угол между осью и радиус-вектором точки изменяется по OX OM 2 M закону (t ) = 3 t (1 + µ ) + (0), 3 = fb 3, где f – гравитационная постоянная. Третье тело – это вязкоупругий шар массы m, m µ.
Предполагается, что третье малое по массе тело не влияет на движение первых двух и движется в гравитационном поле, порожденном первыми двумя телами. Его центр масс C движется в плоскости OXY, а вращение относительно центра масс происходит вокруг нормали к этой плоскости, R – радиус-вектор точки C.
Потенциальная энергия внешних гравитационных полей задается функционалом:
fdx µfdx =, (R 1 + (r + u )) (R 2 + (r + u )) 2 где R 1 = M 1C = R + µbe, R 2 = M 2 C = R be, e = (cos, sin,0).
В §3.1 формулируется постановка задачи и выводятся уравнения движения в форме уравнений Рауса с использованием канонических переменных Андуайе-Делоне, описывающих поступательно вращательное движение вязкоупругого шара, и лагранжевых переменных u i (r, t ), i = 1, 2,3, описывающих его деформации. Конфигурационное пространство данной механической системы есть прямое произведение трехмерного дифференцируемого многообразия M и банахова пространства векторных функций u(r, t ). Полученная система уравнений движения содержит малый параметр, обратно пропорциональный модулю упругости шара, и малый параметр µ. При = соответствующие уравнения описывают движение абсолютно твердого шара в ограниченной круговой задаче трех тел. К этим уравнениям можно применить метод усреднения. Переменные “действие”, определяющие орбиту центра масс шара и момент количеств движения относительно центра масс, в усредненных уравнениях в этом случае не изменяются.
В §3.2 методом разделения движений и асимптотическим методом, аналогичным методу Крылова-Боголюбова для систем с быстрыми и медленными переменными, строится приближенная система уравнений, описывающая диссипативную эволюцию переменных “действие”.
Равенству µ = 0 соответствует задача о движении вязкоупругого шара в центральном ньютоновском поле сил, рассмотренная в §2.4. Эта задача имеет в качестве притягивающего множества стационарное решение – гравитационно-стабилизированное вращение вязкоупругого шара на круговой орбите. Показано, что уравнения движения допускают класс квазикруговых орбит, т.е. орбит с нулевым эксцентриситетом. В этом классе движений строится приближенная система уравнений относительно переменных действие, описывающих изменение радиуса орбиты центра масс вязкоупругого шара и его вращательного движения.
Показано, что эволюция движения вязкоупругого шара разбивается на три этапа, характеризующихся различными временами. На этапе “быстрой” диссипативной эволюции, протекающей со скоростью порядка (нулевое приближение по малому параметру µ ), уравнения, описывающие эволюцию переменных действие, имеют вид:
& & I = 18D2 2 2 ( 2 1 ), = 18D2 2 2 ( 2 1 ), 4 (21) где I =| K | в случае прямого вращения шара и I = | K | в случае обратного вращения шара, K – вектор момента количеств движения шара относительно центра масс, f 2m3 f I = fm R, 1 =, 2 =, 3 =, b A A - момент инерции недеформированного шара относительно диаметра.
Система уравнений (21) имеет асимптотически устойчивое стационарное решение, соответствующее режиму гравитационной стабилизации шара на круговой орбите. В µ -окрестности указанного стационарного решения учитывается первое приближение по малому параметру µ.
Приближенные уравнения, описывающие эволюцию переменных I, с учетом членов первого приближения по малому параметру µ, имеют вид (для ясности изложения здесь сохранены обозначения для переменных действие):
& I = 18D2 2 2 {( 2 1 ) + 2µ2 a(q)}, & = 18D2 2 2 {( 2 1 ) + 2µ 2 a(q)}, (22) 1 q cos 1 b (1 2q cos + q 2 ) 3 где a (q ) = d, q=.
2 R При q = 1 интеграл, определяющий функцию a (q ), расходится. При q 1, т.е. для внешних орбит вязкоупругого шара, функция a (q ) принимает положительные значения, а при q 1, т.е. для внутренних орбит вязкоупругого шара, функция a (q ) принимает отрицательные значения.
Из системы (22) следует, что и в первом приближении по малому параметру µ сохраняется момент количеств движения G0 относительно & точки O : G0 = 0, G0 = I +. Система уравнений (22) также как и система (21), обладает асимптотически устойчивым стационарным решением, соответствующим движению шара, при котором его центр масс движется по орбите радиуса R 01, отличного от R0, с угловой скоростью 1 = 2 (1+ 2 µa (q )).
Таким образом, на первом этапе “медленной” диссипативной эволюции, протекающей со скоростью порядка µ, орбита вязкоупругого шара незначительно смещается к орбите тела массы µ, т.е.
для внутренних орбит вязкоупругого шара радиус квазикруговой орбиты увеличивается, а для внешних – уменьшается, стремясь к новому своему стационарному значению R01.
В µ 2 -окрестности этого асимптотически устойчивого стацио нарного решения учитывается следующее приближение по малому параметру µ. Получено уравнение, описывающее эволюцию момента & количеств движения G0 в виде: G0 =18µ 2 D2 2 2 (q)b(q). Показано, что эволюционная система уравнений с учетом членов второго порядка по малому параметру µ стационарных решений не имеет. На втором этапе “медленной” диссипативной эволюции, протекающей со скоростью порядка µ 2, орбита вязкоупругого шара удаляется от орбиты тела массы µ : для внутренних планет радиус квазикруговой орбиты уменьшается, а для внешних – увеличивается.
Полученные на основе приближенных уравнений результаты справедливы для областей изменения радиусов квазикруговых орбит, исключающих окрестность точки b, в которой могут находиться точки либрации, так как в этом случае оказываются некорректными процедура усреднения по быстрой угловой переменной и разложение в ряд одного из гравитационных потенциалов.
В четвертой главе диссертации изучается эволюция поступательно вращательного движения вязкоупругой планеты при взаимодействии приливов, порождаемых на ней центральным телом и спутником.
Двойная планета, моделируемая материальной точкой массы µ и однородным деформируемым вязкоупругим шаром массы m, движется в гравитационном поле неподвижного центра с массой, равной ( µ m 1 ). Предполагается, что движение центра масс деформируемой планеты и материальной точки происходит в одной плоскости, проходящей через притягивающий центр, а вращение вязкоупругой планеты относительно ее центра масс происходит вокруг нормали к этой плоскости. При этом радиус орбиты центра масс двойной планеты много больше взаимного расстояния между составляющими ее элементами.
В §4.1 выводятся уравнения движения в форме уравнений Рауса.
Канонические переменные Андуайе-Делоне используются для описания орбиты центра масс двойной планеты и орбит вязкоупругого шара и материальной точки массы µ относительно их общего центра масс, а также для описания вращательного движения вязкоупругого шара.
Функции, описывающие деформации шара, составляют лагранжеву часть переменных. Показано, что уравнения движения допускают класс квазикруговых орбит, т.е. орбит с нулевыми эксцентриситетами.
Дальнейшее исследование проводится в этом классе движений. Указанная упрощенная модель позволяет выявить наиболее важные закономерности эволюционных процессов, происходящих в рассматриваемой механической системе.
В соответствии с общей постановкой задачи, сформулированной в первой главе, жесткость вязкоупругого шара предполагается большой и вводится малый параметр, обратно пропорциональный модулю упругости. При = 0 шар представляет собой абсолютно твердое тело, и исходная задача распадается на три, не связанные друг с другом задачи:
задачу о движении центра масс двойной планеты относительно притягивающего центра, задачу о движении связки двух планет относительно их общего центра масс и задачу о вращении шара относительно собственного центра масс. Рассеяние энергии при деформациях вязкоупругого шара вызывает “перекрестные” связи между этими независимыми задачами и является причиной эволюции рассматриваемой системы.
В §4.2 методом разделения движений и усреднения осуществляется построение приближенной системы эволюционных уравнений, которая представляет собой систему обыкновенных дифференциальных уравнений 3-го порядка относительно переменных действие, описывающих изменение расстояния от притягивающего центра до центра масс двойной планеты, расстояния от центра масс вязкоупругого шара до материальной точки массы µ и вращения вязкоупругого шара.
Указанная система уравнений имеет вид:
& 1 = 18 2 D21 (1 3 ), µ 2 & 2 = 18 2 D 2 ( 2 3 ), (m + µ ) & & & I = 1 2. (23) Переменные 1, 2 выражаются через радиус орбиты центра масс двойной планеты и расстояние между центром масс R1 R деформируемого шара и спутником с помощью соотношений:
2 2 µm R1 =, R2 =, mr =, f 0 = f (m + µ ).
1 m+µ f (m + µ ) 2 f 0 mr Переменная I равна модулю момента количеств шара относительно центра масс, взятого со знаком плюс в случае прямого вращения, и со знаком минус в случае обратного вращения. Величины 1, 2, 3, входящие в правые части уравнений (23), имеют вид:
f 2 (m + µ ) 3 f 02 m r I 1 =, 2 =, 3 =, 3 3 A 1 где A - момент инерции недеформированного шара относительно диаметра, f - универсальная гравитационная постоянная. Уравнения (23) имеют первый интеграл: I + 1 + 2 = G 0, который является следствием закона сохранения момента количеств движения рассматриваемой механической системы относительно притягивающего центра.
В §4.3 на основе полученной в §4.2 системы уравнений проводится качественное исследование эволюции движения системы. Находятся стационарные решения, и исследуется их устойчивость на основе уравнений в вариациях. Показано, что в зависимости от величины модуля момента количеств движения G0 существует либо два, либо одно, либо ни одного стационарного решения. В стационарном движении (в случае его существования) конфигурация двойной планеты неизменна в орбитальной системе координат, связанной с движением ее центра масс.
Все стационарные решения являются неустойчивыми.
~ Построены фазовые портреты в плоскости переменных ( 1, 2 ), (m + µ ) 4/ ~ где 2 = 2. Соответствующая система уравнений имеет вид:
µm k ~ & 1 = 12 [ 3 G 0 + 1 + n 2 ], 1 k & ~ ~ 2 = ~12 [ ~ 3 G0 + 1 + n 2 ], 2 µ где = 18 2 D2 f 8 (m + µ )12 A 1, =, m(m + µ ) 2 / µm k = Af 2 (m + µ ) 3.
n=, (m + µ ) 4/ Движение по фазовым траекториям происходит таким образом, что, в конце концов, либо двойная планета приближается к притягивающему центру, либо спутник (тело массы µ ) падает на планету. Численным методом получена фазовая траектория для двойной планеты Земля-Луна, движущейся в поле притяжения Солнца, в рамках рассматриваемой в этой главе постановки задачи. Эволюция движения системы Земля-Луна происходит так, что в настоящее время радиусы орбит R1 и R увеличиваются, а вращение Земли вокруг собственной оси замедляется.
После того как угловая скорость вращения Земли сравняется с орбитальной в ее движении относительно Луны (период обращения Луны вокруг Земли при этом согласно полученным результатам увеличится в 1,476 раз и составит 40,32 современных суток) расстояние между Луной и Землей начнет уменьшаться. Следует отметить, что часть интегральной кривой, соответствующая уменьшению R 2, расположена вблизи кривой, на которой угловая скорость вращения Земли совпадает с орбитальной в ее движении относительно Луны. Подобную картину эволюции Земли и Луны описывал Дарвин и позже Макдональд. По последним оценкам увеличение суток за счет приливного трения составляет 0,003 с за 100 лет (Сидоренков Н.С., 2004).
В пятой главе диссертации изучается эволюция поступательно вращательного движения симметричного спутника с гибкими вязкоупругими стержнями в центральном ньютоновском поле сил.
Спутник представляет собой динамически симметричное твердое тело, по оси симметрии которого расположена пара гибких вязкоупругих стержней. При отсутствии деформаций в стержнях главные центральные моменты инерции спутника равны между собой, т.е. его центральный эллипсоид инерции – сфера. Диссипация энергии происходит за счет внутреннего трения в материале стержней. Жесткость стержней предполагается большой, и вводится малый параметр, обратно пропорциональный жесткости упругих элементов.
В §5.1 формулируется постановка задачи и выводятся уравнения движения в форме уравнений Рауса. Используется линейная теория изгиба тонких нерастяжимых стержней. Канонические переменные Андуайе-Делоне описывают поступательно-вращательное движение спутника, а лагранжевы переменные u1 ( s, t ), u 2 ( s, t ) – деформации стержней (здесь s – координата точки стержня по оси симметрии спутника, t – время). Так как поле внешних сил центральное, то имеет место закон сохранения момента количеств движения относительно притягивающего центра.
В §5.2 методом разделения движений и усреднения строится приближенная система обыкновенных дифференциальных уравнений относительно переменных “действие” и медленной угловой переменной g – долготы перигелия от восходящего узла, описывающая эволюцию движения рассматриваемой механической системы.
В §5.3 определяются многообразия стационарных решений усредненной системы уравнений, и исследуется их устойчивость на основе уравнений в вариациях. Стационарные точки эволюционной системы уравнений соответствуют стационарным движениям системы, на которых диссипация энергии равна нулю. Для одного класса стационарных движений спутника его центр масс описывает эллиптическую орбиту, стержни перпендикулярны плоскости орбиты, а скорость вращения вокруг оси симметрии, совпадающей со стержнями, произвольна. Описанное стационарное движение устойчиво практически только для прямых вращений спутника с угловой скоростью большей среднего движения по орбите.
Другой класс стационарных решений соответствует движению центра масс по круговой орбите, когда угловая скорость вращения спутника совпадает с орбитальной, а, вообще говоря, деформированные стержни имеют произвольную ориентацию в орбитальной системе координат. Так же, как и в задаче о поступательно-вращательном движении вязкоупругого шара в центральном ньютоновском поле сил, рассмотренной в главе 2, в зависимости от величины модуля момента количеств движения G 0 число стационарных орбит может быть равно нулю, одному или двум. В случае двух стационарных орбит устойчивым является стационарное решение, соответствующее орбите большего радиуса.
В окрестности этого стационарного решения условие отсутствия резонансов нарушается, поэтому его исследование проводится на основе возмущенной системы уравнений, полученной в результате применения метода разделения движений. Точные уравнения движения допускают такой класс движений, при котором эксцентриситет и наклонение орбиты равны нулю, а ось вращения ортогональна плоскости орбиты. В указанном классе движений проводится исследование устойчивости стационарных решений, соответствующих режиму гравитационной стабилизации спутника на круговой орбите. Стационарные движения, соответствующие такому положению спутника, при котором стержни лежат в плоскости орбиты и направлены по радиус-вектору центра масс, неустойчивы. А стационарные решения, соответствующие расположению стержней по касательной к орбите, устойчивы.
В шестой главе диссертации изучается эволюция движения симметричного спутника с гибкими вязкоупругими стержнями на круговой орбите. Спутник моделируется осесимметричным твердым телом, по оси симметрии которого расположена пара гибких вязкоупругих стержней.
В §6.1 формулируется постановка задачи и выводятся уравнения движения. Для определения радиус-вектора точки стержня в связанной со спутником системе координат используется линейная теория изгиба тонких нерастяжимых стержней. Вводится малый параметр, пропорциональный обратной величине жесткости упругих элементов, и малый параметр µ, пропорциональный квадрату орбитальной угловой скорости. Для описания вращательного движения спутника относительно центра масс используются канонические переменные Андуайе.
В §6.2 методом разделения движений осуществляется построение “возмущенной” системы уравнений относительно переменных Андуайе.
При µ = 0 указанная система уравнений описывает быструю диссипативную эволюцию вращательного движения симметричного твердого тела с двумя вязкоупругими стержнями относительно центра масс, причем правые части этих уравнений не зависят от угловых переменных. Переменные I 2, I 3 ( I 2 – модуль вектора кинетического момента K вращательного движения, I 3 – проекция вектора K на нормаль к плоскости орбиты) в этих уравнениях не эволюционируют, а скорость эволюции переменной I 1 ( I 1 – проекция вектора K на ось симметрии недеформированного спутника) имеет порядок. В случае динамически вытянутого спутника ( A C ) значение переменной I уменьшается и стремится к нулю, а в случае динамически сжатого спутника ( A C ) значение переменной I 1 увеличивается и стремится совпасть с I 2. Эволюционная система уравнений, описывающая быструю диссипативную эволюцию, имеет вид:
( A C) ( I 2 I 1 ) I 1, I 2 = 0, I 3 = 0 (d 1 0 ).
& & & I 1 = d 1 2 AC В §6.3 изучается этап медленной диссипативной эволюции в случае A C. Асимптотическим методом, аналогичным методу Крылова Боголюбова для систем с быстрыми и медленными переменными, в окрестности аттрактора I 1 = 0, I 2 = I 2 (0), I 3 = I 3 (0) получены уравнения, описывающие эволюцию переменных I 2, I 3 с учетом периодических возмущений, вызываемых переменной I 1. Скорость эволюции этих переменных имеет порядок µ 2. С сохранением обозначений для переменных действие эволюционная система уравнений имеет вид:
&2 = 4 n1 2 ( I 2 )1 + I 3 2 I 3, I I I 4 &3 = n1 8 2 ( I 2 ) I 3 + 3 I 3 6 I 3 5, I I I2 I2 µ 2 d 1 2 0 C 2 A 2, 2 ( I 2 ) = I 2 A, = где n1 =, = fm1, f 16 R гравитационная постоянная, m1 масса притягивающего центра.
Построен фазовый портрет в плоскости переменных x = I 3 I 2 1 и y = I 2 A 1, который аналогичен фазовому портрету, полученному Белецким В.В.
(1975) при анализе вращательного движения сферически симметричного твердого тела под действием приливного момента. Система дифференциальных уравнений относительно переменных x и y имеет единственное стационарное решение x = 1, y = ( – орбитальная угловая скорость), которое соответствует гравитационной стабилизации спутника в орбитальной системе координат. Соответствующая стационарная точка в плоскости ( x, y ) является устойчивым узлом. На основе неусредненной возмущенной системы уравнений получен следующий результат: положение равновесия, при котором ось симметрии спутника направлена по радиус-вектору его центра масс, устойчиво, а положение равновесия, при котором ось симметрии спутника направлена по касательной к орбите, неустойчиво.
В §6.4 изучается этап медленной диссипативной эволюции в случае A C. Точные уравнения движения допускают такой класс решений, когда вектор кинетического момента K направлен по оси симметрии ( I 1 = I 2 ). В этом классе движений асимптотическим методом Крылова Боголюбова строится приближенная система уравнений, описывающая эволюцию переменных I 2, I 3. Скорость эволюции переменных I 2, I 3 имеет порядок µ 2. Соответствующие уравнения имеют вид:
2 & = n 1 I 3 1 + 3 I 3 I 2 4 I 3, 2 I2 I 2 I2 C I 2 & = 4 n 1 I 3 I 3 1 + I 3, 2 I3 I 2 C I где n 2 = µ 2 d1 2 0.
Построен фазовый портрет системы в плоскости переменных ( x, y ).
Анализ эволюционной системы уравнений позволяет сделать вывод, что в случае A C движение спутника стремится к прямому вращению вокруг оси симметрии, направленной по нормали к плоскости орбиты. Значение предельной угловой скорости определяется начальными условиями.
В 7-ой главе диссертации рассмотрен пространственный вариант задачи о двойной планете: исследуется поступательно-вращательное движение вязкоупругой планеты в гравитационном поле притягивающего центра и спутника, которые моделируются материальными точками.
Спутник и планета движутся относительно общего центра масс, которые в свою очередь совершают движение относительно притягивающего центра.
В §7.1 из вариационного принципа Даламбера-Лагранжа выводится точная система уравнений движения рассматриваемой механической системы в рамках линейной модели теории вязкоупругости.
В §7.2 методом разделения движений строится возмущенная система уравнений, описывающая поступательно-вращательное движение планеты и спутника с учетом возмущений, вызываемых упругостью и диссипацией.
В §7.3 проводится анализ квазистатических деформаций вязкоупругой планеты, моделируемой однородным изотропным вязкоупругим шаром. Деформированное состояние вязкоупругой планеты, вызываемое полем гравитационных сил и сил инерции, описывается вектором упругого смещения u = u 1 = (u 10 + u 11 + u 120 u 120 + u 130 u 130 ), & & (24) где 2 fm u 10 = 2 3 [d 1 r 2 + d 2 r02 ]r, 3 r { } u 11 = a1 2 r 2 (, r ) 2 r + (a 2 r 2 + a 3 r02 ) 2 r (, r ), 1 1 6 3 {a 1 r 1 (,r) r + (a r + a r )1 r (,r) }, 3 fµ u 120 = 2 2 2 1 2 2 30 2 3 6 2 R {a 1 r 1 (,r) r + (a r + a r )1 r (,r) }, 3 f = 2 2 2 u 130 1 1 2 30 1 3 6 2 Q - вектор угловой скорости системы координат, интегральным образом связанной с шаром, Q1 = Q 1, R2 = R 2, Q 1 = OD, R 2 = DF, D - центр масс деформируемого шара, F - материальная точка, моделирующая спутник, 1, 2 - единичные векторы, направленные по Q1 и R соответственно, коэффициенты a i, d j, i = 1,2,3, j = 1,2 являются дробно рациональными функциями коэффициента Пуассона, = 1, модуль Юнга, r0 - радиус шара в естественном недеформированном состоянии, - его плотность. Зависимость от времени вектора упругого смещения u = u (r, t ) осуществляется через, 1, 2 из Q1, R 2, невозмущенной системы уравнений, описывающей поступательно вращательное движение абсолютно твердого шара в гравитационном поле притягивающего центра и спутника.
Получено уравнение поверхности вращающейся планеты, определяемое слагаемыми u 10, u 11 в правой части (24). Найдены главные центральные моменты инерции вращающейся планеты в виде:
A = A0 [1 + 2 r02 f1 ( ) fmr01 f 2 ( )], C = A0 [1 + 2 r02 f 3 ( ) fmr01 f 2 ( )], 4(1 2 )(4 3 ) 95 3 13 2 167 + где f 1 ( ) =, f 2 ( ) =, 35(1 )(5 + 7) 35(1 ) 2(25 3 39 2 61 + 59) f 3 ( ) =, A0 = mr02, m - масса планеты.
35(1 )(5 + 7) Получены значения модуля вектора упругого смещения, позволяющие оценить величину приливных горбов, создаваемых притягивающим центром и спутником и определяемых слагаемыми u 120, u 130 в формуле (24).
В §7.4 найдены стационарные движения и исследована их устойчивость на основе уравнений в вариациях. Показано, что в стационарном движении, когда диссипация энергии равна нулю, центр масс планеты и спутник движутся по круговым орбитам относительно притягивающего центра, располагаясь с ним на одной прямой.
Вязкоупругая планета в стационарном движении неподвижна в орбитальной системе координат. Показано, что указанное стационарное движение является неустойчивым.
В §7.5 в качестве модели вязкоупругой планеты рассматривается система, состоящая из абсолютно твердой невесомой сферы, к которой с внешней стороны жестко прикреплена вязкоупругая сферическая оболочка, а внутри имеется подвижное внутреннее ядро. Получено квазистатическое решение задачи теории упругости, которое описывает деформации упругого слоя планеты под действием гравитационных сил и сил инерции. Показано, что действие внешних гравитационных полей в рассматриваемом приближении вызывает центрально симметричные деформации упругого сферического слоя планеты, а колебания внутреннего ядра нарушают эту симметрию.
В заключении сформулированы основные результаты диссертационной работы.
ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ 1) Шатина А.В. Об асимптотических свойствах решений одного класса механических систем с бесконечным числом степеней свободы // Вестник МГУ, сер. 1, матем.-мех., №4, 1990, с. 85-89.
2) Вильке В.Г., Шатина А.В. Эволюция движения симметричного спутника с гибкими вязкоупругими стержнями в центральном ньютоновском поле сил // Космические исследования, т.37, вып.3, 1999, с.
289-295.
3) Вильке В.Г., Шатина А.В. Эволюция движения вязкоупругого шара в ограниченной круговой задаче трех тел // Прикладная математика и механика, т.64, вып.5, 2000, с. 772-782.
4) Шатина А.В. Эволюция движения вязкоупругого шара в центральном ньютоновском поле сил // Космические исследования, 2001, т.39, №3, с.303-315.
5) Вильке В,Г., Шатина А.В. Эволюция движения двойной планеты // Космические исследования, 2001, т.39, №3, с. 316-323.
6) Вильке В.Г., Шатина А.В. Эволюция движения вязкоупругого шара в ограниченной круговой задаче трех тел // Нелинейная механика. М.:
Наука, Физматлит, 2001, с. 385-401.
7) Шатина А.В. Эволюция вращательного движения симметричного спутника с гибкими вязкоупругими стержнями // Космические исследования, 2002, т.40 №2, с. 178-192.
8) Вильке В.Г., Шатина А.В. Медленная диссипативная эволюция движения вязкоупругого шара в ограниченной круговой задаче трех тел //Прикладная математика и механика, т.66, вып.5, 2002, с. 782-792.
9) Шатина А.В. Эволюция поступательно-вращательного движения вязкоупругого шара // Механика твердого тела, ISSN 0321-1975, 2002, Вып.32, с. 194-202.
10) Шатина А.В. Эволюция движения вязкоупругой планеты // Новая геометрия природы. Труды международной научной конференции (август Т. Астрономия. Образование.
25-сентябрь 5, 2003) III.
Естественнонаучная философия. Казань: Казанский гос. университет им.
В.И. Ульянова-Ленина, 2003, с. 166-170.
11) Шатина А.В. Быстрая и медленная диссипативная эволюция в механических системах, содержащих вязкоупругие элементы // Известия РАН, Механика твердого тела, №2, 2004, с. 14-23.
12) Вильке В.Г., Шатина А.В. О поступательно-вращательном движении вязкоупругого шара в гравитационном поле притягивающего центра и спутника // Космические исследования, 2004, т. 42, №1, с 95-106.
13) Шатина А.В. О деформациях планеты, содержащей подвижное внутреннее ядро, в гравитационном поле центрального тела и спутника // Известия РАН, Механика твердого тела, №1, 2005, с. 3-12.