Содержание к диссертации
Введение
Глава 1. Математические модели углового движения космического аппарата
1.1. Динамические и кинематические уравнения Эйлера 15
1.2. Кватернионные уравнения движения осесимметричного космического аппарата в осцилляторной и нормальной формах 16
Глава 2. Построение аналитических решений кватернионных уравнений движения космического аппарата в частных случаях его движения
2.1. Аналитическое решение уравнений движения космического аппарата в случае его движения по инерции 23
2.2. Аналитическое решение уравнений движения космического аппарата в случае коллинеарности момента внешних сил кинетическому моменту .29
Глава 3. Построение прогнозного движения космического аппарата
3.1. Постановка задачи построения обобщенного прогнозного движения осесимметричного космического аппарата 36
3.2. Решение задачи 37
3.3. Аналитические частные решения задачи 40
3.4. Численное решение задачи построения прогнозного движения осесимметричного космического аппарата 46
3.4.1. Описание алгоритма численного решения 46
3.4.2. Результаты численного решения 48
3.5. Постановка задачи построения прогнозного движения космического аппарата произвольной динамической конфигурации 51
3.6. Решение задачи 52
3.7. Численное решение задачи построения прогнозного движения космического аппарата произвольной конфигурации 53
3.7.1 Описание алгоритма численного решения 53
3.7.2. Результаты численного решения 56
Глава 4. Построение законов оптимального изменения вектора кинетического момента осесимметричного космического аппарата
4.1. Построение законов оптимального изменения вектора кинетического момента осесимметричного космического аппарата с использованием уравнений, записанных во вращающейся системе координат 61
4.1.1. Постановка задачи 61
4.1.2. Необходимые условия оптимальности 65
4.1.3. Условия трансверсальности 69
4.1.4. Анализ задачи 72
4.1.5. Аналитическое решение задачи 76
4.2. Построение оптимальных законов изменения вектора кинетического момента осесимметричного космического аппарата с использованием уравнений, записанных в инерциальной системе координат 86
4.2.1. Постановка задачи 86
4.2.2. Необходимые условия оптимальности 88
4.2.3. Условия трансверсальности 92
4.2.4. Анализ задачи 93
Глава 5. Оптимальное управление угловым движением осесимметричного космического аппарата 97
5.1. Решение задачи оптимального управления угловым движением осесимметричного космического аппарата с использованием кватернионных уравнений движения в нормальной форме 97
5.1.1. Постановка задачи 97
5.1.2. Необходимые условия оптимальности 100
5.1.3. Условия трансверсальности 103
5.1.4. Анализ задачи 106
5.2. Описание алгоритма численного решения задачи оптимального управления угловым движением осесимметричного космического аппарата 109
5.3. Результаты численного решения 112
Заключение 119
Библиографический список 121
- Кватернионные уравнения движения осесимметричного космического аппарата в осцилляторной и нормальной формах
- Аналитическое решение уравнений движения космического аппарата в случае его движения по инерции
- Описание алгоритма численного решения
- Необходимые условия оптимальности
Введение к работе
Задачам динамики и управления угловым движением космического аппарата посвящено большое число публикаций как в России, так и за рубежом. Однако сложность стоящих здесь проблем, отсутствие общих аналитических решений и трудности численного решения дифференциальных краевых задач, к которым сводятся задачи оптимального управления пространственным движением космического аппарата, продолжают оставлять эту проблематику актуальной.
Успех в решении задач динамики и управления движением космического аппарата во многом зависит от выбранной модели движения космического аппарата. Угловое движение космического аппарата, рассматриваемого как твердое тело, описывается двумя группами уравнений: динамическими уравнениями Эйлера и кинематическими уравнениями, записанными в тех или иных кинематических параметрах (углах Эйлера-Крылова, направляющих косинусах, параметрах Эйлера (Родрига-Гамильтона), Кейли-Клейна). В работах В.Н. Кошлякова [40-45], В.Н. Бранеца [15-19], И.П. Шмыглевского [15-18], Ю.Н. Челнокова [79-91], Д.В. Лебедева [51, 118], В. Wie [100], Н.Л. Стрелковой [76, 77] и других авторов показано, что использование в качестве кинематических параметров таких невырождающихся параметров как параметры Эйлера повышает эффективность аналитического исследования и численного решения многих задач динамики и управления угловым движением космического аппарата. При этом удобным математическим аппаратом оказывается аппарат кватернионов Гамильтона.
В большинстве работ по динамике и управлению угловым движением космического аппарата, использующих в качестве параметров ориентации параметры Эйлера, применяется модель углового движения твердого тела, состоящая из динамических уравнений Эйлера и кватернионного кинематического уравнения. Кроме этой модели для решения задач динамики и управления успешно используется другая модель твердого тела, имеющая вид системы дифференциальных уравнений второго порядка относительно параметров Эйлера (В.Н. Кошляков [40-45], Ю.Н. Челноков [79-83]). Эти уравнения движения, записанные в связанной с твердым телом системе координат, удобны для решения ряда задач динамики твердого тела и космического аппарата (например, для решения задач устойчивости и стабилизации углового движения), однако они оказываются неудобными для построения оптимальных (программных) траекторий и управлений угловым движением космического аппарата с использованием принципа максимума.
В диссертационной работе для решения задач динамики и управления угловым движением осесимметричного космического аппарата - твердого тела используются две формы (осцилляторная и нормальная) новых кватернионных уравнений движения осесимметричного твердого тела, записанных в специальной системе координат, вращающейся с абсолютной угловой скоростью, коллинеарной вектору кинетического момента твердого тела. Эти уравнения, предложенные Ю.Н. Челноковым [81,82,87], отличаются от уравнений, записанных в связанной с твердым телом системе координат, большей простотой и компактностью. Их использование, как показано в диссертационной работе, облегчает аналитическое и численное решение изучаемых задач.
Одной из важных задач теории управления угловым движением космического аппарата является задача построения прогнозного движения -одного из основных программных угловых движений космического аппарата, имеющего смысл движения по инерции, удовлетворяющего заданным краевым условиям по угловому положению космического аппарата. Известное решение этой задачи основывается на аналитических решениях уравнений движения твердого тела в случае Эйлера. Такой подход особенно эффективен в случае осесимметричного космического аппарата и изучался в ряде работ (М.В. Невского [52] и др.). В диссертационной работе решается задача построения обобщенного прогнозного движения осесимметричного космического аппарата, когда программный управляющий момент не равен нулю, а коллинеарен вектору кинетического момента космического аппарата. Частным случаем этой задачи является задача построения обычного прогнозного движения. Решение задачи основывается на построенном в работе аналитическом решении кватернионных дифференциальных уравнениях движения осесимметричного космического аппарата - твердого тела в случае коллинеарности главного вектора моментов внешних сил и вектора кинетического момента твердого тела. Этот случай интегрируемости уравнений движения осесимметричного твердого тела впервые был установлен (в классических переменных) Б.А. Смольниковым [74]. Построенная в работе новая, кватернионная форма аналитического решения уравнений движения осесимметричного твердого тела для указанного случая удобна для построения обобщенного прогнозного движения космического аппарата.
В работе также рассматривается новое решение задачи построения обычного прогнозного движения космического аппарата с произвольным распределением масс, основанное на решении дифференциальной краевой задачи, описываемой кватернионным дифференциальным уравнением, связывающим кватернион ориентации космического аппарата с его вектором кинетического момента.
Другой актуальной задачей теории оптимального управления угловым движением осесимметричного космического аппарата, связанной с созданием высокоточных систем ориентации, использующих в качестве исполнительных устройств вращающиеся маховики, является задача построения оптимальных (программных) законов изменения вектора кинетического момента осесимметричного космического аппарата. В диссертационной работе эта задача, по-видимому, впервые решена для двух основных функционалов качества с использованием принципа максимума и кватернионного дифференциального уравнения, связывающего кватернион ориентации космического аппарата с его вектором кинетического момента.
Отметим, что если в задаче построения обычного прогнозного движения космического аппарата с произвольным распределением масс используется кватернионное дифференциальное уравнение, v связывающее кватернион ориентации космического аппарата с его вектором кинетического момента, записанное в системе координат, связанной с космическим аппаратом, то в задаче построения оптимальных законов изменения вектора кинетического момента осесимметричного космического аппарата используется другое кватернионное дифференциальное уравнение, также связывающее кватернион ориентации космического аппарата с его вектором кинетического момента, но записанное в вышеуказанной специальной системе координат.
Большое количество работ посвящено изучению задачи оптимального по быстродействию управления угловым движением космического аппарата -твердого тела, однако эта задача в общей постановке до настоящего времени далека от завершения. Среди работ, посвященных решению этой задачи в различных постановках с использованием аппарата кватернионов, отметим работы В.Н. Бранеца [15-19], И.П. Шмыглевского [15-18], М.Б. Чертока [19], Ю.В. Казначеева [19], Д.В. Лебедева [51], В.В. Маланина [56,57], Н.А. Стрелковой [76,77], А.Н. Сиротина [70-72], А.В. Молоденкова [34,60,61], Я.Г. Сапункова [61].
Обширный обзор литературы [5-9,11,14,17,21,25,19,47,48,49,51,58,63, 70,72,75,92-96,99] по задаче оптимального в смысле быстродействия переориентацией космического аппарата содержится в монографии [58]. В этой монографии отмечается, что при исследовании динамических задач оптимального по быстродействию управления движением объектов, когда в качестве управляющего воздействия выступает вектор моментов внешних сил, основная аналитическая трудность заключается в рассмотрении оптимального по быстродействию управления ориентацией твердого тела. Обычно решение данной задачи разбивается на две подзадачи: торможение (гашение) вращений и переориентация твердого тела (начальное и конечное положения тела заданы и являются состояниями покоя). В работе [58] отмечается, что если торможение вращений достаточно полно исследовано [1, 2, 3, 4, 10, 14, 36, 53, 54, 73, 74, 92, 97], то задача оптимального по быстродействию управления переориентацией твердого тела в полной постановке является задачей нерешенной, тем более далека от завершения общая задача синтеза оптимальных управлений ориентацией твердого тела (с одновременным гашением угловой скорости).
В диссертационной работе задача оптимального управления угловым движением осесимметричного космического аппарата исследуется при произвольных краевых условиях. В качестве функционала качества рассматривается более общий комбинированный функционал, характеризующий расход времени и импульса модуля управляющего момента за промежуток времени управляемого движения на перевод космического аппарата из начального в конечное состояние. Для решения задачи используется новая кватернионная модель углового движения космического аппарата - твердого тела, повышающая эффективность аналитического исследования и численного решения задачи оптимального управления.
Каждая из исследуемых в диссертации задач представляет самостоятельный интерес. Методологическая общность их решения заключается в использовании нового класса кватернионных моделей движения осесимметричного космического аппарата - твердого тела.
Целями и основными задачами работы являются: - получение решений задач динамики и управления угловым движением осесимметричного космического аппарата с использованием новых кватернионных моделей углового движения твердого тела;
- разработка алгоритмов и программ численного решения краевых алгебраических и дифференциальных задач, к которым сводятся задачи построения прогнозного движения и задачи оптимального управления угловым движением космического аппарата.
В работе решаются следующие задачи:
- построение аналитических решений уравнений неуправляемого и управляемого движений осесимметричного космического аппарата в частных случаях его движения;
- построение обобщенного прогнозного движения осесимметричного космического аппарата и обычного прогнозного движения космического аппарата произвольной динамической конфигурации;
- получение оптимальных законов изменения вектора кинетического момента осесимметричного космического аппарата;
- построение оптимальных управлений и траекторий углового движения осесимметричного космического аппарата.
В первой главе диссертации рассматриваются известные классические и кватернионные уравнения углового движения осесимметричного космического аппарата - твердого тела вокруг центра масс, находящегося под действием произвольного внешнего момента, в том числе две формы (осцилляторная и нормальная) новых кватернионных уравнений движения осесимметричного космического аппарата, записанных в специальной системе координат, вращающейся с абсолютной угловой скоростью, коллинеарной вектору кинетического момента космического аппарата.
Во второй главе строятся аналитические решения кватернионных дифференциальных уравнений движения осесимметричного космического аппарата в двух частных случаях его движения: в случае движения космического аппарата по инерции и в случае коллинеарности главного вектора моментов внешних сил и вектора кинетического момента космического аппарата, удобные для решения задач построения прогнозного движения космического аппарата.
В третьей главе диссертационной работы решается задача построения обобщенного прогнозного движения осесимметричного космического аппарата, когда программный управляющий момент не равен нулю, а коллинеарен вектору кинетического момента космического аппарата. Показано, что решение задачи построения обобщенного прогнозного движения осесимметричного космического аппарата сводится к решению системы двух скалярных алгебраических трансцендентных уравнений. Построены аналитические решения задачи в случае сферической симметрии космического аппарата и в двух частных случаях разворотов космического аппарата (на любой угол) вокруг определенным образом ориентированных в пространстве осей. В этой главе также рассмотрено новое решение задачи построения обычного прогнозного движения космического аппарата с произвольным распределением масс, основанное на решении кватернионной дифференциальной краевой задачи. Разработаны алгоритмы и программа численного решения задачи построения прогнозного движения.
Четвертая глава диссертационной работы посвящена задаче построения оптимальных законов изменения вектора кинетического момента осесимметричного космического аппарата с использованием кватернионных уравнений, записанных во вращающейся специальной или инерциальной системе координат. С помощью принципа максимума [67] и кватернионного дифференциального уравнения, связывающего вектор кинетического момента осесимметричного космического аппарата с кватернионом ориентации космического аппарата, сформулированы (для двух функционалов качества) дифференциальные краевые задачи для построения оптимальных законов изменения вектора кинетического момента космического аппарата. Построены общие аналитические решения дифференциальных уравнений краевых задач, образующих системы девяти нелинейных дифференциальных уравнений. Показано, что решение дифференциальных краевых задач сводится к решению двух скалярных алгебраических трансцендентных уравнений. Получены, как явные функции времени, зависимости для кватерниона ориентации, вектора абсолютной угловой скорости и вектора кинетического момента космического аппарата, описывающие оптимальное управляемое движение космического аппарата. Разработаны эффективные алгоритмы и программа численного решения задач построения прогнозного движения космического аппарата для любых заданных начальной и конечной ориентации космического аппарата.
В пятой главе предложено новое решение задачи оптимального управления угловым движением осесимметричного космического аппарата с использованием новых кватернионных уравнений в нормальной форме. Задача оптимального управления угловым движением осесимметричного космического аппарата исследуется при произвольных краевых условиях. В качестве функционала качества рассматривается комбинированный функционал, характеризующий расход времени и импульса модуля управляющего момента за промежуток времени управляемого движения на перевод космического аппарата из начального в конечное состояние. Построены новые кватернионные дифференциальные уравнения краевой задачи принципа максимума для построения оптимальных управлений (проекций вектора управляющего момента) и траекторий углового движения осесимметричного космического аппарата. Найдено несколько первых интегралов задачи, установлен частный случай интегрируемости дифференциальных уравнений краевой задачи. Разработаны кватернионные алгоритмы и программа численного решения задачи оптимального управления угловым движением осесимметричного космического аппарата для произвольных краевых условий. Построены примеры численного решения задачи.
Научная новизна: 1. Построено кватернионное решение дифференциальных уравнений движения осесимметричного космического аппарата в случае коллинеарности главного вектора моментов внешних сил и вектора кинетического момента космического аппарата в форме, удобной для построения прогнозного движения космического аппарата.
2. Получены новые кватернионные алгебраические и дифференциальные уравнения краевых задач для построения обобщенного прогнозного движения осесимметричного космического аппарата и обычного прогнозного движения космического аппарата с произвольным распределением масс. Построены частные аналитические решения краевых задач. Разработаны алгоритмы численного решения краевых задач.
3. Найдено общее аналитическое решение задачи построения оптимальных законов изменения вектора кинетического момента осесимметричного космического аппарата для двух основных функционалов качества.
4. Получены новые кватернионные дифференциальные уравнения краевой задачи принципа максимума для построения оптимальных управлений и траекторий углового движения осесимметричного космического аппарата (минимизируется время управляемого движения и импульс модуля управляющего момента за это время), построено частное аналитическое решение задачи, найдено несколько первых интегралов дифференциальных уравнений краевой задачи, построены алгоритмы численного решения краевой задачи.
Достоверность результатов обеспечивается корректностью математической постановки задач, строгостью применяемых методов решения, сопоставлением полученных результатов с результатами других авторов. На защиту выносятся:
1. Кватернионное аналитическое решение дифференциальных уравнений движения осесимметричного космического аппарата в случае коллинеарности главного вектора моментов внешних сил и вектора кинетического момента осесимметричного космического аппарата.
2. Решения трех задач теории управления угловым движением космического аппарата:
- задачи построения прогнозного движения,
- задачи построения оптимальных законов изменения вектора кинетического момента осесимметричного космического аппарата,
- задачи оптимального управления угловым движением осесимметричного космического аппарата, построенные с использованием новых кватернионных моделей динамики твердого тела.
3. Геометрические интерпретации управляемого углового движения осесимметричного космического аппарата.
4. Алгоритмы численного решения задач построения прогнозного движения космического аппарата (обобщенного прогнозного движения осесимметричного космического аппарата, обычного прогнозного движения космического аппарата с произвольным распределением масс) и задачи оптимального управления угловым движением осесимметричного космического аппарата.
Практическая значимость. Полученные законы оптимального управления и траектории углового движения космического аппарата могут быть использованы в качестве законов программных управлений и программных траекторий при построении систем ориентации космических аппаратов, использующих в качестве исполнительных устройств вращающиеся маховики и (или) реактивные двигатели. Разработанные алгоритмы и программы могут быть использованы для построения программных траекторий и управлений угловым движением космического аппарата - твердого тела и математического моделирования управляемого движения космического аппарата.
Использование результатов. Результаты, полученные в диссертационной работе, были использованы в лаборатории механики, навигации и управления ИПТМУ РАН (г. Саратов, 2003-2004 гг.) при выполнении работ по заданию президиума РАН (тема № 0120.0403260 «Разработка кватернионных и бикватернионных моделей и методов механики твердого тела, методов пространства состояний в задачах динамики и управления движением») и проекта РФФИ № 02-01-00988 «Кватернионные модели и методы в пространственных нелинейных задачах оптимального управления движением космических аппаратов».
Апробация работы. Основные результаты работы докладывались: на научных конференциях механико-математического факультета Саратовского государственного университета «Актуальные проблемы математики и механики» (2001-2004); на XXXIV и XXXV постоянно действующем научно-техническом семинаре «Проблемы теории, конструкции, проектирования и эксплуатации ракет» (г. Саратов, 2002, 2003); на научных семинарах Института
Кватернионные уравнения движения осесимметричного космического аппарата в осцилляторной и нормальной формах
Четвертая глава диссертационной работы посвящена задаче построения оптимальных законов изменения вектора кинетического момента осесимметричного космического аппарата с использованием кватернионных уравнений, записанных во вращающейся специальной или инерциальной системе координат. С помощью принципа максимума [67] и кватернионного дифференциального уравнения, связывающего вектор кинетического момента осесимметричного космического аппарата с кватернионом ориентации космического аппарата, сформулированы (для двух функционалов качества) дифференциальные краевые задачи для построения оптимальных законов изменения вектора кинетического момента космического аппарата. Построены общие аналитические решения дифференциальных уравнений краевых задач, образующих системы девяти нелинейных дифференциальных уравнений. Показано, что решение дифференциальных краевых задач сводится к решению двух скалярных алгебраических трансцендентных уравнений. Получены, как явные функции времени, зависимости для кватерниона ориентации, вектора абсолютной угловой скорости и вектора кинетического момента космического аппарата, описывающие оптимальное управляемое движение космического аппарата. Разработаны эффективные алгоритмы и программа численного решения задач построения прогнозного движения космического аппарата для любых заданных начальной и конечной ориентации космического аппарата.
В пятой главе предложено новое решение задачи оптимального управления угловым движением осесимметричного космического аппарата с использованием новых кватернионных уравнений в нормальной форме. Задача оптимального управления угловым движением осесимметричного космического аппарата исследуется при произвольных краевых условиях. В качестве функционала качества рассматривается комбинированный функционал, характеризующий расход времени и импульса модуля управляющего момента за промежуток времени управляемого движения на перевод космического аппарата из начального в конечное состояние. Построены новые кватернионные дифференциальные уравнения краевой задачи принципа максимума для построения оптимальных управлений (проекций вектора управляющего момента) и траекторий углового движения осесимметричного космического аппарата. Найдено несколько первых интегралов задачи, установлен частный случай интегрируемости дифференциальных уравнений краевой задачи. Разработаны кватернионные алгоритмы и программа численного решения задачи оптимального управления угловым движением осесимметричного космического аппарата для произвольных краевых условий. Построены примеры численного решения задачи.
Научная новизна: 1. Построено кватернионное решение дифференциальных уравнений движения осесимметричного космического аппарата в случае коллинеарности главного вектора моментов внешних сил и вектора кинетического момента космического аппарата в форме, удобной для построения прогнозного движения космического аппарата. 2. Получены новые кватернионные алгебраические и дифференциальные уравнения краевых задач для построения обобщенного прогнозного движения осесимметричного космического аппарата и обычного прогнозного движения космического аппарата с произвольным распределением масс. Построены частные аналитические решения краевых задач. Разработаны алгоритмы численного решения краевых задач. 3. Найдено общее аналитическое решение задачи построения оптимальных законов изменения вектора кинетического момента осесимметричного космического аппарата для двух основных функционалов качества. 4. Получены новые кватернионные дифференциальные уравнения краевой задачи принципа максимума для построения оптимальных управлений и траекторий углового движения осесимметричного космического аппарата (минимизируется время управляемого движения и импульс модуля управляющего момента за это время), построено частное аналитическое решение задачи, найдено несколько первых интегралов дифференциальных уравнений краевой задачи, построены алгоритмы численного решения краевой задачи. Достоверность результатов обеспечивается корректностью математической постановки задач, строгостью применяемых методов решения, сопоставлением полученных результатов с результатами других авторов. На защиту выносятся: 1. Кватернионное аналитическое решение дифференциальных уравнений движения осесимметричного космического аппарата в случае коллинеарности главного вектора моментов внешних сил и вектора кинетического момента осесимметричного космического аппарата. 2. Решения трех задач теории управления угловым движением космического аппарата: - задачи построения прогнозного движения.
Аналитическое решение уравнений движения космического аппарата в случае его движения по инерции
Показано, что решение задачи построения обобщенного прогнозного движения осесимметричного космического аппарата сводится к решению системы четырех алгебраических (трансцендентных) уравнений. Построены аналитические решения задачи в случае сферической симметрии космического аппарата и в двух частных случаях разворотов осесимметричного космического аппарата (на любой угол) вокруг определенным образом ориентированных в пространстве осей.
Сформулирована дифференциальная краевая задача для построения обычного прогнозного движения космического аппарата с произвольным распределением масс, основанная на интегрировании кватернионного уравнения, связывающего кватернион ориентации космического аппарата с вектором кинетического момента космического аппарата.
Разработаны алгоритмы и программа численного решения задач построения прогнозного движения космического аппарата для любых заданных начальной и конечной ориентации космического аппарата. Установлено, что в качестве эффективного начального приближения для численного решения задач служит аналитическое решение задачи в случае сферически симметричного космического аппарата или в случае разворота осесимметричного космического аппарата вокруг его оси симметрии.
Построение законов оптимального изменения вектора кинетического момента осесимметричного космического аппарата Построение законов оптимального изменения вектора кинетического момента осесимметричного космического аппарата с использованием уравнений, записанных во вращающейся системе координат
Исследуем задачу управления угловым движением космического аппарата, рассматриваемого как осесимметричное твердое тело, когда управляющей вектор - функцией является вектор кинетического момента космического аппарата. Эта задача представляет не только теоретический, но и практический интерес в связи с построением систем управления пространственной ориентацией космического аппарата с использованием управляющих маховиков. Для построения такого рода систем необходимо знание оптимальных в том или ином смысле законов изменения вектора кинетического момента космического аппарата, позволяющих решать задачи переориентации космического аппарата. Отметим, что эти законы могут быть реализованы не только с помощью управляющих маховиков, но и с помощью моментов внешних сил или интегралов от них, создаваемых реактивными импульсными газовыми двигателями.
Будем рассматривать угловое движение космического аппарата, используя уравнения движения, записанные в системе координат Y , вращающейся с абсолютной угловой скоростью, коллинеарной вектору кинетического момента космического аппарата. Движение базиса Y относительно системы координат X будем задавать кватернионом z, который связан с кватернионом X, определяющим ориентацию космического аппарата в системе координат X , операцией вращения (1.2.14):
Описание алгоритма численного решения
Построены общие аналитические решения кватернионных дифференциальных уравнений движения осесимметричного космического аппарата в случае движения по инерции и в случае коллинеарности главного вектора моментов внешних сил вектору кинетического момента космического аппарата в форме, удобной для построения прогнозного движения осесимметричного космического аппарата. Для этих случаев движения космического аппарата получены в виде явных функций времени компактные зависимости, описывающие законы изменения кватерниона ориентации и вектора абсолютной угловой скорости космического аппарата. 2. Показано, что решение задачи построения обобщенного прогнозного движения осесимметричного космического аппарата сводится к решению системы двух скалярных алгебраических трансцендентных уравнений. Построены аналитические решения задачи в случае сферической симметрии космического аппарата и в двух частных случаях разворотов космического аппарата (на любой угол) вокруг определенным образом ориентированных в пространстве осей. Сформулирована дифференциальная краевая задача для построения обычного прогнозного движения космического аппарата с произвольным распределением масс, основанная на использовании кватернионного дифференциального уравнения, связывающего кватернион ориентации космического аппарата с вектором кинетического момента космического аппарата. Разработаны эффективные алгоритмы и программа численного решения задач построения прогнозного движения космического аппарата для любых заданных начальной и конечной ориентации космического аппарата. Установлено, что начальные приближения, найденные с помощью построенных частных аналитических решений, обеспечивают высокую сходимость итерационных вычислительных процессов. 3. С помощью принципа максимума и кватернионного дифференциального уравнения, связывающего вектор кинетического момента осесимметричного космического аппарата с кватернионом ориентации космического аппарата, сформулированы (для двух функционалов качества) дифференциальные краевые задачи для построения оптимальных законов изменения вектора кинетического момента космического аппарата. Построены общие аналитические решения дифференциальных уравнений краевых задач, образующих системы девяти нелинейных дифференциальных уравнений. Показано, что решение дифференциальных краевых задач сводится к решению двух скалярных алгебраических трансцендентных уравнений. Получены, как явные функции времени, зависимости для кватерниона ориентации, вектора абсолютной угловой скорости и вектора кинетического момента космического аппарата, описывающие оптимальное управляемое движение космического аппарата. 4. Построены новые кватернионные дифференциальные уравнения краевой задачи принципа максимума для построения оптимальных управлений (проекций вектора управляющего момента) и траекторий углового движения осесимметричного космического аппарата (минимизируется комбинированный в отношении времени и энергетических затрат функционал качества). Найдено несколько первых интегралов задачи, установлен частный случай интегрируемости дифференциальных уравнений краевой задачи, имеющий место при выполнении условий коллинеарности векторов MY , LY , Ф и п у . Разработаны кватернионные алгоритмы и программа численного решения задачи оптимального управления угловым движением осесимметричного космического аппарата для произвольных краевых условий. Построены примеры численного решения задачи. 5. Даны наглядные геометрические интерпретации управляемого углового движения осесимметричного космического аппарата.
Необходимые условия оптимальности
Большое количество работ посвящено изучению задачи оптимального по быстродействию управления угловым движением космического аппарата -твердого тела, однако эта задача в общей постановке до настоящего времени далека от завершения. Среди работ, посвященных решению этой задачи в различных постановках с использованием аппарата кватернионов, отметим работы В.Н. Бранеца [15-19], И.П. Шмыглевского [15-18], М.Б. Чертока [19], Ю.В. Казначеева [19], Д.В. Лебедева [51], В.В. Маланина [56,57], Н.А. Стрелковой [76,77], А.Н. Сиротина [70-72], А.В. Молоденкова [34,60,61], Я.Г. Сапункова [61].
Обширный обзор литературы [5-9,11,14,17,21,25,19,47,48,49,51,58,63, 70,72,75,92-96,99] по задаче оптимального в смысле быстродействия переориентацией космического аппарата содержится в монографии [58]. В этой монографии отмечается, что при исследовании динамических задач оптимального по быстродействию управления движением объектов, когда в качестве управляющего воздействия выступает вектор моментов внешних сил, основная аналитическая трудность заключается в рассмотрении оптимального по быстродействию управления ориентацией твердого тела. Обычно решение данной задачи разбивается на две подзадачи: торможение (гашение) вращений и переориентация твердого тела (начальное и конечное положения тела заданы и являются состояниями покоя). В работе [58] отмечается, что если торможение вращений достаточно полно исследовано [1, 2, 3, 4, 10, 14, 36, 53, 54, 73, 74, 92, 97], то задача оптимального по быстродействию управления переориентацией твердого тела в полной постановке является задачей нерешенной, тем более далека от завершения общая задача синтеза оптимальных управлений ориентацией твердого тела (с одновременным гашением угловой скорости).
В диссертационной работе задача оптимального управления угловым движением осесимметричного космического аппарата исследуется при произвольных краевых условиях. В качестве функционала качества рассматривается более общий комбинированный функционал, характеризующий расход времени и импульса модуля управляющего момента за промежуток времени управляемого движения на перевод космического аппарата из начального в конечное состояние. Для решения задачи используется новая кватернионная модель углового движения космического аппарата - твердого тела, повышающая эффективность аналитического исследования и численного решения задачи оптимального управления. Каждая из исследуемых в диссертации задач представляет самостоятельный интерес. Методологическая общность их решения заключается в использовании нового класса кватернионных моделей движения осесимметричного космического аппарата - твердого тела. Целями и основными задачами работы являются: - получение решений задач динамики и управления угловым движением осесимметричного космического аппарата с использованием новых кватернионных моделей углового движения твердого тела; - разработка алгоритмов и программ численного решения краевых алгебраических и дифференциальных задач, к которым сводятся задачи построения прогнозного движения и задачи оптимального управления угловым движением космического аппарата. В работе решаются следующие задачи: - построение аналитических решений уравнений неуправляемого и управляемого движений осесимметричного космического аппарата в частных случаях его движения; - построение обобщенного прогнозного движения осесимметричного космического аппарата и обычного прогнозного движения космического аппарата произвольной динамической конфигурации; - получение оптимальных законов изменения вектора кинетического момента осесимметричного космического аппарата; - построение оптимальных управлений и траекторий углового движения осесимметричного космического аппарата. В первой главе диссертации рассматриваются известные классические и кватернионные уравнения углового движения осесимметричного космического аппарата - твердого тела вокруг центра масс, находящегося под действием произвольного внешнего момента, в том числе две формы (осцилляторная и нормальная) новых кватернионных уравнений движения осесимметричного космического аппарата, записанных в специальной системе координат, вращающейся с абсолютной угловой скоростью, коллинеарной вектору кинетического момента космического аппарата. Во второй главе строятся аналитические решения кватернионных дифференциальных уравнений движения осесимметричного космического аппарата в двух частных случаях его движения: в случае движения космического аппарата по инерции и в случае коллинеарности главного вектора моментов внешних сил и вектора кинетического момента космического аппарата, удобные для решения задач построения прогнозного движения космического аппарата. В третьей главе диссертационной работы решается задача построения обобщенного прогнозного движения осесимметричного космического аппарата, когда программный управляющий момент не равен нулю, а коллинеарен вектору кинетического момента космического аппарата. Показано, что решение задачи построения обобщенного прогнозного движения осесимметричного космического аппарата сводится к решению системы двух скалярных алгебраических трансцендентных уравнений.
Похожие диссертации на Кватернионное решение задач динамики и управления угловым движением осесимметричного космического аппарата
