Баллистико-навигационное проектирование полетов к Луне, планетам и малым телам Солнечной системы тема автореферата и диссертации по механике, 01.02.01 ВАК РФ

004618882

На правах рукописи

ТУЧИН Андрей Георгиевич

Баллистико-навигационное проектирование полётов к Луне, планетам и малым телам Солнечной системы

Специальность 01.02.01 - Теоретическая механика

Автореферат диссертации на соискание учёной степени доктора физико-математических наук

1 3 ЯН В 2071

Москва-2010

004618882

Работа выполнена в Учреждении Российской академии наук Институте прикладной математики им. М.В. Келдыша РАН.

Официальные оппоненты:

Академик РАН, доктор физико-математических наук, профессор Маров Михаил Яковлевич, Учреждение Российской академии наук Институт геохимии и аналитической химии им. В.И. Вернадского РАН.

Доктор технических наук, профессор Почукаев Владимир Николаевич, Федеральное государственное унитарное предприятие «Центральный научно-исследовательский институт машиностроения».

Доктор технических наук, профессор Константинов Михаил Сергеевич, Государственное образовательное учреждение высшего профессионального образования «Московский авиационный институт (государственный технический университет)».

Ведущая организация:

Учреждение Российской академии наук Институт космических исследований

Защита диссертации состоится 15 февраля 2011 г. в 11 час. на заседании диссертационного совета Д 002.024.01 при Институте прикладной математики им. М.В. Келдыша РАН по адресу 125047, Москва, Миусская пл., 4, конференц-зал.

С диссертацией можно ознакомиться в библиотеке Института прикладной математики им. М.В. Келдыша РАН.

Автореферат разослан « »2010г.

РАН.

Учёный секретарь диссертационного совета доктор физико-математических наук

Полилова Т.А.

1. Общая характеристика работы

Диссертация посвящена теоретико-механическим вопросам проектирования полётов космических аппаратов (КА) к планетам Солнечной системы и их естественным спутникам. Разработан метод проектирования квазисинхронных орбит КА вокруг Фобоса, предназначенных для обеспечения посадки КА на его поверхность. Разработаны методы определения параметров движения КА по траекторным измерениям на фоне работы двигательной установки; разработаны методы расчёта точности приведения КА к планете-цели при использовании химических и электроракетных двигателей с учётом ошибок прогноза движения КА, ошибок исполнения манёвров и коррекций, и различных возмущений, вносимых в движение центра масс КА. Разработанные методы опираются на опыт баллистико-навигационного обеспечения полётов к Венере: КА «Венера-15,16», «Вега-1,2», к Марсу: КА «Фобос-2» и нашли применение в проекте «Фобос-Грунт», при проектировании полётов к Луне в проектах «Луна-Ресурс» и «Луна-Глоб», к Венере (проект «Венера-Д»), в систему Юпитера с целью посадки КА на его естественный спутник Европу. Это определяет актуальность и практическую значимость диссертации.

В диссертации использованы методы небесной механики, асимптотического интегрирования дифференциальных уравнений, математической статистики и теории управления движением.

Цели диссертации

Цель работы состоит в разработке теоретико-механических и математических методов, обеспечивающих баллистику и навигацию в проекте «Фобос-Грунт», и распространению этих методов на решение задач навигации и управления в перспективных отечественных проектах полётов к Луне (проекты «Луна-Ресурс» и «Луна-Глоб»), к Венере (проект «Венера-Д») и в систему Юпитера (проект «Лаплас»),

Основные результаты

1. Разработан метод проектирования квазисинхронных орбит КА вокруг Фобоса, позволяющий вычислить параметры орбиты КА, которая в заданное время проходит над заданной долготой поверхности Фобоса и обладает свойством минимального дрейфа, т.е. проходит заданную долготу примерно на одинаковом расстоянии от поверхности Фобоса. Метод позволил предварительно вычислить набор таблиц, каждая из которых соответствует значению максимального удаления от поверхности Фобоса. На основе этих таблиц простой вычислительный алгоритм определяет параметры орбиты КА, которая обладает свойством минимального дрейфа и в заданное время проходит над заданной долготой поверхности Фобоса.

2. Разработаны методы бортовой навигации, обеспечивающие сближение с малым небесным телом и посадку на его поверхность. Методы предназначены для их реализации на бортовом компьютере, предполагают использование цифровой модели поверхности небесного тела, используют измерения дальности и скорости в направлении лучей навигационных приборов и позволяют работать на фоне работающих двигателей. Алгоритмы и бортовые программы навигационного обеспечения участка посадки проекта «Фобос-Грунт» разработаны на основе этих методов.

3. Разработан метод определения параметров движения КА в условиях воздействия шума. Метод основан на представлении модели движения КА в виде композиции опорного движения и движения относительно опорного. Движение относительно опорного представляет собой случайный процесс с нулевым средним и описывается линейным стохастическим дифференциальным уравнением. Оценка вычисляется из условия минимизации функционала, содержащего квадраты взвешенных невязок между измеренными значениями и их расчётными моделями, а также взвешенные квадраты модулей векторов суммарных возмущений между измерениями. Метод применён в алгоритмах определения параметров

движения КА на фоне работы двигательной установки. На основе этих методов разработаны алгоритмы определения параметров движения КА с электроракетными двигательными установками (ЭРДУ).

4. Разработан метод оценки точности приведения КА в заданную область в пространстве целевых параметров с учётом ошибок исполнения и навигационных ошибок. В качестве целевых параметров могут быть выбраны параметры картинной плоскости при решении задачи приведения КА к планете или Луне; высота над поверхностью, вертикальная и горизонтальная составляющая вектора скорости при решении задачи обеспечения условий начала сеанса посадки на поверхность Луны; параметры, определяющие геометрию гало-орбиты, при решении задачи перелёта на гало-орбиту. Метод использован при оценке точности приведения КА в проектах «Фобос-Грунт», «Луна-Ресурс», «Луна-Глоб», «Венера-Д».

Научная новизна

1. Новизна разработанного метода проектирования квазисинхронных орбит КА вокруг Фобоса состоит в том, что усреднённая система учитывает дрейф эллипса, по которому совершается облёт КА вокруг Фобоса, в направлении, которое ортогонально орбитальному движению Фобоса. Наличие этого уравнения позволило найти начальные приближения для орбит с минимальным дрейфом, а затем построить искомые семейства орбит.

2. Новизна методов бортовой навигации на участке посадки состоит в том, что найдены такие модели движения КА и построены такие модели измерений, которые, с одной стороны, обладают простотой, позволяющей реализовать их на бортовом компьютере, а с другой стороны, точностью, позволяющей решить поставленную задачу.

3. Научная новизна метода оценки параметров движения КА в условиях воздействия шума состоит в предложенной модели движения, представляющей собой композицию детерминированного опорного движения и случайного процесса, который описывается линейным стохастическим

дифференциальным уравнением, а также в методе, который обеспечивает одновременное получение оценки вектора состояния КА и векторов суммарных возмущений между измерениями.

Апробация работы

Результаты, представленные в диссертации, докладывались автором на 16-м симпозиуме IFAC по автоматическому управлению (Санкт-Петербург, июнь 2004 г); на 17-м (Москва, июнь 2003г) и 18-м (Мюнхен, октябрь 2004г.) Международных симпозиумах по динамике космического полёта; на Общероссийском семинаре «Современные методы навигации и управления движением» (Институт проблем Управления РАН, 31 марта 2009 г.); на семинаре «Солнечная система и смежные проблемы физики и механики» (ИПМ им. М.В. Келдыша РАН, 18 марта 2008 г.).

Публикации

Список основных публикаций по теме диссертации приведён в третьем разделе автореферата. Список состоит из 18 работ; 7 работ из этого списка опубликованы в изданиях, рекомендованных ВАК; 5 из них - в соавторстве. Из работ, выполненных с соавторами, в диссертацию включены только результаты, полученные автором.

Структура и объём работы

Диссертация состоит из введения, пяти глав, заключения, списка использованных источников. Работа изложена на 238 страницах, содержит 28 рисунков, 33 таблицы. Список использованных источников включает 141 наименование.

2. Краткое описание содержания работы

В главе 1 рассмотрены задачи проектирования квазисинхронных орбит КА вокруг Фобоса для решения задачи посадки на его поверхность. Результаты, представленные в этой главе, описаны в работах: [1; 2; 8; 10; 13].

Моделирование показало [3; 5], что требуется четыре дня полёта КА на квазисинхронной орбите (КСО) для определения параметров движения с необходимой точностью. Поэтому интервал пребывания на КСО оценивается пятью сутками. К номинальной КСО предъявляется ряд требований. За время пребывания на КСО с борта КА должно быть получено телевизионное изображение предполагаемого района посадки и проверена работоспособность лазерного высотомера автономной системы посадки. В момент начала сеанса посадки, при получении изображения района посадки и проверки работоспособности лазерного высотомера, удаление от поверхности Фобоса не должно превосходить 60 км и должны выполняться условия по освещённости и радиосвязи со станциями слежения в Евпатории и Уссурийске. Если в штатном сеансе работы автономной системы посадка не будет осуществлена и КА останется на прежней орбите, должна быть обеспечена возможность повторения попытки посадки на Фобос.

Выбор КСО предполагается проводить по следующей схеме. Из условий освещённости, обеспечения радиосвязи со станциями слежения, возможно, других условий выбирается точка над поверхностью Фобоса и время, в которое КСО должна пройти через эту точку. Из семейства квазисинхронных орбит, проходящих через заданную точку в заданное время, выбирается такая КСО, которая при других прохождениях долготы выбранной точки удаляется от неё на минимальное расстояние. Это требование должно обеспечить выполнение условия по максимальному удалению КА от поверхности Фобоса при каждом прохождении над районом посадки. Тем самым возможность телевизионной съемки предполагаемого участка посадки и условия проверки работоспособности лазерного высотомера будут обусловлены только условиями освещённости и радиосвязи со станциями слежения.

Специалисты Института геохимии и аналитической химии им. В.И. Вернадского РАН указали несколько интересных для исследования точек поверхности Фобоса. Координаты точек посадки задаются а системе

координат, центр которой находится в центре фигуры Фобоса. Положение точек местности определяется сферическими координатами. Широты считаются положительными к северу от экватора. Долгота измеряется в экваториальной плоскости к западу от нулевого меридиана. Сеанс посадки целесообразно начинать с упреждением достижения КА долготы точки посадки. Величина этого упреждения может находиться в диапазоне от 1° до 60°. В момент начала сеанса посадки должны быть обеспечены условия освещённости и радиосвязи с наземными станциями слежения. По моменту времени прохождения над заданной долготой всегда можно определить момент времени прохождения долготы 270°. Поэтому рассмотрим задачу построения такой КСО, которая проходит в заданное время над точкой с долготой 270° и удалена от поверхности на расстояние в диапазоне от 50 до 60 км. Потребуем от этой КСО, чтобы при других прохождениях КА долготы 270° его удаление от начальной точки было минимальным.

Движение по КСО представлено как облёт Фобоса по дрейфующему эллипсу. Такое описание уже использовалось в работах М. Хенона1, Д. Бенеста2, А.Ю. Когана3, М.Л. Лядова и М.А. Вашковьяка4 для получения качественных результатов и развито в диссертации в целях получения инструмента для проектирования посадки на Фобос. Большая полуось этого эллипса ориентирована вдоль орбитального движения Фобоса. КА обегает эллипс в направлении, обратном к орбитальному движению Фобоса, с периодом, близким к периоду обращения Фобоса вокруг Марса. Дрейф

1 Henon М. Numerical Exploration of the Restricted Problem // VI. Hill's Case:

Non-Periodic Orbits. Astronomy and Astrophysics. 1970. V. 24. № 9. P. 24-36.

3 Benest D. Libration effects for retrograde satellites in the restricted three-body problem // Cel. Mech. 1976. V. 13. № 2. P. 203-215.

3 Коган А.Ю. Далекие спутниковые орбиты в ограниченной круговой задаче трёх тел // Космич. исслед. 1988. Т. 26. № 6. С. 813-818.

4 Лидов М.Л., Вашковьяк М.А. О квазиспутниковых орбитах для эксперимента по уточнению гравитационной постоянной II Письма в АЖ. 1994. Т. 20. № 3. С. 229-240.

эллипса также проходит вдоль оси, ориентированной по орбитальному движению Фобоса.

Для обеспечения успешной посадки нужна такая КСО, дрейф которой был бы минимален. На рис. 1. показаны КСО с обычным и минимальным дрейфом. Эти орбиты показаны во вращающейся СК, в центре которой находится Фобос. Ось X направлена по линии визирования Марс - Фобос. Ось У ортогональна оси X, лежит в плоскости орбиты Фобоса и направлена в сторону его орбитального движения. Так как нулевой меридиан Фобоса постоянно обращен в сторону Марса, точкам, расположенным на отрицательной части оси У, в этой системе координат соответствуют подспутниковые точки с долготой, близкой к 270°.

У, КМ

У, КМ

X, КМ

X, КМ

Рис. 1 — Квазисинхронные орбиты с обычным (а) и минимальным (б) дрейфом

Найденные КСО с минимальным дрейфом обладают тем важным свойством, что прохождение над точкой посадки всегда происходит на одной

и той же высоте. Это позволит при необходимости, если процесс посадки не был начат, начать его в одном из следующих прохождений над районом посадки, повторить посадку, произвести телевизионную съёмку района посадки, выполнить измерения высоты до поверхности Фобоса при его облёте.

Быстродействие алгоритма расчёта КСО было получено за счёт использования начального приближения, формируемого из эволюционных уравнений.

При поиске КСО в качестве упрощённой модели движения использовалась плоская эллиптическая задача Хилла в безразмерных переменных. Плоская эллиптическая задача Хилла обладает тем свойством, что траектории, полученные в ней, сохраняют свою структуру при переходе к полной модели сил. Круговая задача Хилла таким свойством не обладает.

Эволюционное движение системы плоского кругового движения задачи Хилла в условиях удаления тела нулевой массы от тела меньшей массы на расстояния, значительно превышающие радиус сферы Хилла, изучено в работе Бенеста2. В этой работе получены эволюционные уравнения движения центра эллипса вдоль оси Y. Пространственный случай эволюционного движения в аналогичных условиях исследован в работе Когана3. В этих работах с целью изучения эволюции вдоль оси У при упрощении системы уравнений была исключена составляющая, описывающая эволюционное движение вдоль X. Тем самым, в них отброшено множество движений, представляющих практический интерес. Это множество движений рассмотрено в диссертации.

Решая задачу построения системы эволюционных уравнений для плоской эллиптической задачи Хилла, обладающих заданным свойством, нужно специальным образом выбрать фазовые переменные. В диссертации в качестве таких переменных выбраны: большая полуось эллипса, по которому происходит движение КА (А); разность фаз между движением КА вокруг

У

Рис. 2 - Траектории движения точки нулевой массы (КА) и центра перемещающегося эллипса. Мгновенные положения эллипсов, по которым движется К А, показаны для двух моментов времени, когда и + <р = О (а) и и + ф = Зтт/2 (б). Центры этих эллипсов находятся в точках внутреннего кольца, из которых начинаются стрелки

Фобоса и Фобоса вокруг V. (<р); полуось эллипса, по которому движется центр эллипса движения КА (а), разность фаз между движением центра эллипса и дьижением КА по эллипсу (с,). Геометрическая интерпретация

выбранных переменных показана на рисунке 2. Центр внешнего эллипса находится на внутреннем эллипсе. Ищутся такие начальные условия, при которых движение КА происходит в узкой кольцеобразной области. Предварительный анализ показал, что при искомых движениях не только КА, но и центр эллипса движется в узкой кольцеобразной области. На рис. 2 - это внешняя и внутренняя кольцеобразные области. КА имеет обратное движение по отношению к движению Фобоса вокруг Марса^

КА движется по перемещающемуся эллипсу, центр которого движется в противофазе по отношению к КА. Поэтому КА всегда находится в области пересечения перемещающегося эллипса и внешнего кольца. Большая полуось эллипса практически не меняется. Перемещение эллипса происходит за счёт движения его центра в пределах внутреннего кольца. При этом значения а и А близки к своим средним значениям, а значение близко к л.

Усреднение системы дифференциальных уравнений в указанных выше переменных позволяет получить уравнение, связывающее средние значения полуосей внутреннего и внешнего эллипсов. Одним из основных элементов при построении системы эволюционных уравнений является представление

интеграла — — через фазовые параметры. В диссертации найдено новое 2п I г

представление этого интеграла через выбранные фазовые переменные а, Ъ, А, Ф и свободный параметр к:

где К(к) - полный эллиптический интеграл первого рода, к - его модуль,

1 ги»

(1)

¿ = -№¡=5111^, \Ч2=СО$£,.

Усреднённые уравнения движения имеют вид:

. ___1__л__| 3 | 1 ( 4 д! | 2 31

акпд^' + 4к + 2к{АтдАт+ аткдат

т ~т \ т т т т

А 2 91 ■ 2 31

Для поиска начального приближения потребуем, чтобы =

При к = 2 условие на преобразуется в уравнение, связывающее А и а:

ЮАа2 -5Ага +—К

'А

2У

л

Корни уравнения имеют действительные значения при:

= 0. (3)

(4)

Уравнение (3) определяет зависимость а(А), которая используется при построении начального приближения. Рассмотрим алгоритм формирования начального приближения.

Пусть заданы: хо, Уо, и. По переменным \ и ср выполняется сканирование. На каждом шаге сканирования по значениям и, х0, у0, ^ и <р определяются а и А. Только в случае, если А\ < Да, формируется начальное приближение

и вычисляется ширина кольца.

Сформированное приближение используется численным алгоритмом, определяющим начальные условия. Численный алгоритм построен на основе комбинированного применения методов градиентного спуска и золотого сечения.

В результате вычислены таблицы начальных условий для максимальных удалений от поверхности Фобоса в диапазоне от 50 до 200 км.

Для начальных удалений 50, 55 и 60 км кольцеобразные области, в которых находятся КСО, имеют ширину: 5.3, 6.7 и 8.3 км соответственно. Облёт КА вокруг Фобоса на таких орбитах происходит быстрее по отношению к периоду обращения Фобоса вокруг Марса. Средняя

относительная разность угловых скоростей находится в диапазонах: 0.215— 0.234, 0.182-0.199 и 0.154-0.170 для КСО с удалениями 50, 55 и 60 км соответственно.

Глава 2 посвящена алгоритмам навигации и управления в схеме посадки на Фобос. Результаты, представленные в этой главе, описаны в работах [4] и [11]. В этой главе содержится модель движения КА относительно Фобоса, которая достаточно проста для бортовых вычислений, обеспечивает требуемую точность и позволяет рассчитывать коррекции и решать навигационные задачи. Далее строится алгоритм оценки вектора состояния КА по измерениям лазерного высотомера и доплеровского измерителя скорости и дальности с использованием цифровой модели поверхности Фобоса, разработанной специалистами Института геохимии и аналитической химии им. В.И. Вернадского РАН. Ниже приводится анализируемая в диссертации схема посадки и условия ее выполнения.

Для обеспечения условий посадки КА должен быть приведён в точку, которая находится на высоте не более 60 км над предполагаемым районом посадки. Должны быть выполнены следующие условия при подготовке и на интервале посадки:

- получение телевизионного изображения района посадки за несколько суток до начала сеанса посадки;

- прогноз движения КА относительно Фобоса на момент схода с КСО с ошибками, не превосходящими 3 км по положению и 1 м/с по скорости;

- возможность повторения сеанса посадки, если он не был начат в назначенное время;

- реализация схода с КСО, сближение и мягкая посадка в рамках тяговооруженности и запасов топлива КА;

- проверка работоспособности основных бортовых систем, обеспечивающих посадку, до её начала;

- нахождение угла Солнце - Фобос - КА в диапазоне от 20° до 70° в течение сеанса посадки;

- радиовидимость КА со станций слежения в Уссурийске и Медвежьих Озёрах;

- радиовидимость Земли в допустимом диапазоне углов привода остронаправленной антенны (ОНА).

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

- бесплатформенный инерциальный блок (БИБ);

- звёздный прибор (БОКЗ-МФ);

- лазерный высотомер (ЛВ);

- доплеровский измеритель скорости и дальности (ДИСД);

- телевизионная система.

Бесплатформенный инерциальный блок имеет в своем составе три акселерометра и три датчика угловой скорости. Измерения приращений кажущейся скорости в направлении осей чувствительности акселерометров и приращения углов поворота вокруг осей чувствительности датчиков угловой скорости используются для определения ориентации, угловых скоростей и вектора состояния КА методом счисления пути.

Звёздный прибор обеспечивает определение ориентации КА по звёздам с высокой точностью. Параметры ориентации, получаемые звёздным прибором, используются для коррекции параметров ориентации, получаемых методом счисления пути.

Лазерный высотомер обеспечивает измерения расстояния до поверхности по четырём лучам. Он выполняет навигационные измерения на участке от момента схода с КСО до высоты 500 м. Доплеровский измеритель скорости и дальности выполняет измерения дальности до подстилающей поверхности в направлении четырёх лучей и проекции вектора скорости КА относительно Фобоса на направления этих лучей. Доплеровская система

работает с высоты 3 км. Измерения лазерного высотомера и доплеровской системы используются навигационной задачей, входящей в бортовой комплекс управления посадкой, для определения вектора состояния КА, которым корректируется вектор состояния, получаемый методом счисления пути. Измерения дальностей по лучам, получаемые ЛВ и ДИСД, используются также для определения нормали к подстилающей поверхности.

Основной задачей телевизионной системы является определение района в окрестности номинальной точки посадки, в которой рельеф поверхности максимально подходит для посадки. Телевизионная камера позволяет получать навигационные измерения, в том числе скорость движения КА относительно поверхности. Эти измерения являются резервными на случай отказа доплеровской системы.

Схема посадки на Фобос (рис. 3) потребовала анализа траектории с использованием четырёх участков:

- схода с КСО;

- перелёта с КСО в точку, расположенную над районом посадки;

- вертикального спуска;

- прецизионного торможения.

Участок перелёта с КСО в точку, расположенную над районом посадки, начинается с манёвра схода с КСО и завершается в момент попадания в заданную точку. В ходе перелёта предусмотрены коррекции траектории.

При движении КА на участке вертикального спуска предложено использовать простой метод компенсации горизонтальной составляющей скорости. Если она превосходит пороговое значение, включается двигатель, который её компенсирует. При этом в горизонтальной плоскости накапливается ошибка по положению. Поэтому скорость спуска должна быть достаточно велика, чтобы за время спуска не накопилась большая величина ошибки.

Рис. 3 - Схема управляемой посадки в СК, связанной с фигурой Фобоса

На участке прецизионного торможения постепенно гасится вертикальная составляющая скорости до величины, с которой допускается соприкосновение с поверхностью: 1.5-2 м/с. При этом боковая составляющая вектора скорости не должна превосходить ] м/с.

Управление ориентацией КА должно быть организовано следующим образом. До начала сеанса посадки выполняется разворот, обеспечивающий заданную ориентацию. Она определяется из условий:

- «захвата» поверхности Фобоса лучами ЛВ;

- освещённости солнечных батарей в ходе и после посадки;

- работы привода ОНА, обеспечивающей связь с Землей.

В ходе посадки необходимо обеспечить совмещение средней нормали к поверхности Фобоса с продольной осью (ОХ) связанной системы координат

У часто к торможения

Импульс перехода на участок: .......вертикального спуска

2-я кор

1-я коррекция рекция

Сход с КСО

(СК) КА. Это позволяет привести КА на поверхность Фобоса так, чтобы его ось ОХ была направлена по нормали к поверхности, и при этом избежать больших разворотов, которые могут создать сложности в работе ЛВ, ДИСД и привода ОНА.

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

- по длительности всей операции посадки;

- по максимальным значениям остаточных величин вертикальной и горизонтальной составляющих скорости в момент посадки (порядка единиц м/с);

- по высоте, начиная с которой нельзя включать двигатели, обеспечивающие тягу по вертикали, чтобы не испортить оптические условия наблюдения поверхности;

- по энергетическим затратам.

С точки зрения надёжности работы бортовой аппаратуры, чем короче сеанс посадки, тем лучше, так как в сеансе посадки используется большое число аппаратных средств и надёжность всей системы падает с увеличением времени. Основная часть времени сеанса посадки приходится на перелёт с КСО в точку, расположенную над выбранным районом.

С точки зрения затрат характеристической скорости, чем больше это время, тем меньше суммарные затраты на импульс схода с КСО и импульс гашения продольной составляющей скорости в момент перехода в заданную точку.

Разработанная система управления посадкой содержит два варианта схода с КСО: при прохождении долготы траверза точки посадки и с упреждением прохождения долготы точки посадки. В любом из этих вариантов импульс схода с КСО в номинальном случае обеспечивает приведение КА в прицельную точку, находящуюся на заданной высоте (не выше 10 км) над точкой предполагаемой посадки. С этой прицельной точки

начинается участок вертикального спуска. Моделирование процесса посадки показало, что длительность перелёта от момента схода с КСО до начала участка вертикального спуска составляет -30 мин.

При получении от телевизионной системы достоверных данных о неровности рельефа подстилающей поверхности допускается горизонтальный манёвр КА в направлении более благоприятной точки посадки.

Выбор точки посадки, варианта схода с КСО, высоты приведения над точкой посадки, расчёт импульса схода с КСО выполняется на Земле. Перед выполнением сеанса посадки на борт передаются следующие данные, необходимые для обеспечения посадки:

- момент схода с КСО и длительность перелёта от момента схода до прицельной точки;

- импульс схода и положение прицельной точки;

- заданная ориентация КА, векторы состояния КА и Фобоса на момент схода с КСО;

- моменты возможных коррекций на интервале перелёта.

Алгоритм решения навигационной задачи и алгоритм расчёта коррекций используют уравнения движения КА относительно Фобоса. Уравнения записаны в его орбитальной СК, центр которой совпадает с центром масс Фобоса. Ось 0x1 направлена по линии визирования Марс - Фобос. Ось 0x2 ортогональна оси Ох\, лежит в плоскости орбиты Фобоса и направлена в сторону его орбитального движения. Ось Охз дополняет систему до правой. Орбитальная СК Фобоса удобна ещё и тем, что переход из неё в СК, связанную с фигурой Фобоса, выполняется линейным преобразованием, коэффициенты которого постоянны на интервале посадки. Коэффициенты этого линейного преобразования вычисляются на Земле при подготовке сеанса посадки. Уравнения движения КА относительно Фобоса в нормированных переменных имеют вид:

( ^ V Iх У

= + р

сЬ Л ¿2х2 _

<Л>2 Л |х|

+ А:яр3 («г, соэи + я2 вши),

= -2—- - + кар3 (-а, эт -о + аг соб и), (5)

= -рх3 (есови +1) - рц-^ + £ар3а3, «о |х|

где е, р, с - эксцентриситет, фокальный параметр и интеграл площадей орбиты Фобоса; а = (а„а2,а3)Т - вектор кажущегося ускорения в орбитальной

1 и р3

СК Фобоса, р =-; ц = ———; ка = цм и (хр|1 - гравитационные

1 + есозг> Им+Н-рь с

параметры Марса и Фобоса.

Векторы положения хР[; и скорости х№ КА в орбитальной СК связаны с

векторами положения и скорости в нормированных переменных следующими

простыми соотношениями:

с

Р

1 ейс

---ьеБти-х

р сЬ

(6)

Для расчёта импульса схода с КСО и решения навигационной задачи необходим расчёт частных производных текущего вектора состояния КА по начальному вектору состояния. Для этого на бортовой машине нужно интегрировать уравнения в вариациях, что позволяет делать достаточно простой вид системы уравнений (5).

Расчёт импульса коррекции производится итерационно. На каждом шаге итерационного процесса уточняется вектор приращения скорости, обеспечивающий приведение КА в точку, заданную в орбитальной СК Фобоса. Для этого используется матрица, связывающая в линейном приближении невязки по компонентам положения на момент г с поправками к вектору скорости на момент времени (0, которыми компенсируются невязки по положению. Матрица вычисляется с использованием результатов интегрирования уравнений в вариациях. Для расчёта коррекции необходимо

не более трёх итераций. Первая итерация позволяет определить вектор приращения скорости с точностью до единиц сантиметров в сек, а третья - с точностью до долей миллиметров в сек.

Навигационная задача (определение вектора состояния КА по измерениям лазерного высотомера и доплеровской системы) решается классическим методом наименьших квадратов. Измеряемыми функциями являются расстояние до поверхности Фобоса вдоль заданного направления и проекция вектора скорости на заданное направление. Алгоритм решения навигационной задачи основан на минимизации функции, представляющей собой сумму квадратов взвешенных невязок расчётных у""'0 и измеренных \|/(теа5 значений:

Z(x) = (х - Х0)Т с;1 (х- х0) + £ vv, (vr - vr10)2, (?)

1=1

где х0 - априорное значение уточняемого вектора, Со - ковариационная матрица, соответствующая х0, Wj - веса измеряемых параметров.

Поиск минимума функции приводит к решению нелинейной системы уравнений методом Ньютона, который сводится к серии последовательных приближений (не более 5). При программной реализации этого алгоритма в БВК расчёты распределены по тактам управления (не более 80 по 0.1 с).

Основная сложность при поиске минимума функции (7) состоит в нахождении точки пересечения луча с цифровой моделью поверхности Фобоса, что необходимо при получении расчётных значений, соответствующих измерениям JIB и ДИСД. Основная идея алгоритма поиска пересечения луча, направленного от КА, с поверхностью Фобоса, показана на рис. 4. Пересечение луча с цифровой поверхностью Фобоса (ЦПФ) ищется итерационно. Первая точка ть, на луче определяется как пересечение луча с референц-эллипсоидом (РЭ). Вторая и последующие точки rbi+1,/^l вычисляются при выполнении следующих действий:

- на ЦПФ определяется точка г51., лежащая на нормали к РЭ, проходящей через точку гь,;

- через найденную точку г5, проводится плоскость, ортогональная нормали к РЭ;

- находится точка гь (+, пересечения луча с этой плоскостью.

Алгоритмы и методы, описанные в главе 2, были реализованы в результате совместной работы специалистов ИПМ им. М.В. Келдыша и НПО им. С.А. Лавочкина при непосредственном участии автора диссертации в системе управления КА «Фобос-Грунт», разрабатываемой в НПО им. С.А. Лавочкина.

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

и фоновой. Основная задача работает в рамках такта управления КА. Фоновая задача выполняет расчёты, распределённые по тактам управления. Тестирование программного обеспечения управления посадкой выполнялось на созданном для этой цели стенде. Стенд содержит макет бортовой машины и персональный компьютер для моделирования внешней среды и систем КА. В систему управления на макете бортовой машины поступают измерения, получаемые моделями измерительных приборов: БИБ, ЛВ, ДИСД, БОКЗ-МФ. Система управления определяет номера включаемых двигателей. В персональном компьютере моделируется работа двигательной установки с учетом переходных процессов и определяются векторы тяги и момента, создаваемые двигателями. Далее интегрируются уравнения движения центра масс КА, Фобоса и уравнения движения КА вокруг центра масс. Полученные параметры используются для моделирования выхода измерительных приборов.

Разработанные модели и алгоритмы навигации, коррекции и посадки достаточно универсальны, и их можно использовать при создании систем посадки КА на малые спутники планет и астероиды. Полученный опыт может быть использован при разработке системы посадки на Луну и другие небесные тела.

В главе 3 рассмотрены вопросы определения параметров движения КА по результатам траекторных измерений при наличии немоделируемых действующих ускорений.

Задача определения параметров движения космического аппарата (КА) является одной из основных задач, решаемых в ходе управления его полётом. При решении этой задачи часто возникает ситуация, в которой определение параметров движения КА надо выполнять на фоне работы двигателей. В качестве примера можно привести следующие задачи: контроль участка выведения КА на орбиту искусственного спутника; оперативная оценка исполнения импульсов по измерениям наземных средств на фоне работы

двигательной установки; определение параметров движения КА с электроракетной двигательной установкой (ЭРДУ).

К вопросам, которые рассматриваются в диссертации, в первую очередь относятся проблемы решения задач навигации и управления полётом с включённой ЭРДУ. Фактическое ускорение, создаваемое ЭРДУ, отличается от модели этого ускорения, заложенного в расчёты. Имеются ошибки величины и ориентации вектора тяги ЭРДУ в пространстве.

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

Результаты, представленные в данной главе, описаны в [7;12;18].

При решении указанных выше задач применяются различные модели шума и, соответственно, используются различные методы и алгоритмы оценки вектора состояния. Могут применяться комбинированные методы осреднения.

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

- метод наименьших квадратов;

- расширенный фильтр Калмана.

Многолетняя практика показала, что метод наименьших квадратов является очень надёжным методом определения параметров. Оценка получается в результате поиска минимума функционала, представляющего собой сумму взвешенных невязок между измеренными значениями и их расчётными аналогами. При этом расчётные аналоги функционально зависят от уточняемых параметров. Особенностью метода является то, что модель

движения КА должна быть достаточно точной. Не допускается наличие больших возмущений, которые не задаются в виде зависимостей от уточняемых параметров. Поэтому применение метода наименьших квадратов в случае неизвестных интервалов работы двигателя малой тяги в данном случае вызывает определённые проблемы. Метод наименьших квадратов может быть успешно использован на этапе, когда уже имеется оценка моментов включения и выключения двигателя и создаваемого им ускорения.

Расширенный фильтр Калмана является эмпирическим расширением фильтра Калмана для линейных систем на нелинейный случай. Расширенный фильтр Калмана предполагает наличие неизвестного фазового шума, воздействующего на систему. Характеристики шума задаются его ковариационной матрицей. Особенностью этого метода является то, что текущие значения оцениваемых параметров должны находиться в достаточно близкой окрестности относительно их истинных значений. Кроме того, метод требует, чтобы не было длительных интервалов времени, в которых нет измерений.

В диссертации для оценки вектора состояния в условиях воздействия шума разработан метод, который обеспечивает оценку вектора состояния и суммарных возмущений между измерениями. Предполагается, что фазовый шум, воздействующий на систему, является белым или может быть получен из белого шума формирующими фильтрами. Пример использования таких фильтров для формирования шума, вызванного ошибками исполнения программы ЭРДУ, рассмотрен в пятой главе.

Суть предложенного метода состоит в том, что оценка определяется из условия минимизации функционала, зависящего как от невязок измеренных значений и их расчётных аналогов, так и от величин определяемых возмущений. Минимизация функционала выполняется итерационно. На каждом шаге итерационного процесса определяется поправка к искомым параметрам. Определение поправки производится из условия минимума

функционала для линейной системы. Поиск минимума функционала для линейной системы приводит к двум последовательностям рекуррентных формул. Первая последовательность рекуррентных формул идёт от первого измерения к последнему измерению и позволяет определить поправку на момент последнего измерения. Вторая последовательность рекуррентных формул идёт от последнего значения к первому и позволяет восстановить возмущения. Первая последовательность рекуррентных формул эквивалентна рекуррентным формулам фильтра Калмана для линейной системы. Вторая последовательность рекуррентных формул в литературе называется сглаживанием.

Представим модель движения К А на интервале в виде суммы

опорного движения и движения относительно опорного.

хЛ0 = хв(0 + х,(0, (8)

где

^ яг

+ В(0<Лу, (9)

х=хо(0

х - вектор состояния;

Р(?,х) ~ вектор-функция;

В(?) - матрица, описывающая воздействие шума на систему;

- случайный процесс с независимыми приращениями, производная которого представляет белый шум £,(?) с нулевым математическим ожиданием и заданной матрицей интенсивности 0(/). Начальные условия для (8) задаются априорным вектором х0 и его

ковариационной матрицей Р0.

Математическая модель в форме (8) проще, чем модель в форме

нелинейного стохастического векторного дифференциального уравнения. При

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

реализациями этой модели, имеют такие же математические ожидания, как и

для модели, представляемой стохастическим дифференциальным уравнением.

В моменты времени производятся измерения, модели

которых описываются функциями Ч^/Ч',,...,1?,,. Измеряемыми функциями

могут быть наземные или бортовые траекторные измерения. В случае наемных траекторных измерений - это измерения наклонной дальности, радиальной скорости или интерферометрические измерения. В случае бортовых траекторных измерений используют показания акселерометров, углы, определяющие направление на небесное тело (планету, спутник планеты или астероид), получаемое по телевизионному изображению, а также высоту над небесным телом. Измеренное значение функции Ч*, обозначим через (¥,) • Для каждого момента времени справедливо

(Ю)

где Г),- - случайный вектор, имеющий нулевое математическое ожидание и ковариационную матрицу Я,.

Запись в качестве параметра х(-) функции означает, что функция Ч*, зависит не от мгновенного значения вектора состояния, а от функции х(?), которая может быть представлена в виде (8).

Зависимость между векторами состояния хД/,-) в дискретные моменты времени может быть выражена разностным уравнением:

'м

= + Ф ('м.ОЫО-МФ 1 Ф(';+1,т)ВШт)Л, (11)

Ч

где Ф(/,т) - переходная матрица линейной модели уравнений (9).

'м

Обозначим случайный вектор | Ф(?;+1,т)В(т)4(т) г/т как V,. Этот

случайный вектор имеет нулевое математическое ожидание и ковариационную матрицу 'м

<1,= / Ф(?,+1,т)В(т)д(т)Вт(т)Фт(/,+1,т) с1т. (12)

ч

Положим = Тогда вектор начальных условий

х1,(гд,) = хД^) = х(глг) однозначно определяет значения вектор-функции

Х£ О \

Хг

х^ (/) в дискретных точках: . Однако при вычислении значений

функции Ч'Д/^хД-)) нужно знать зависимость х.,(7) в окрестности каждого момента времени . Представим эту зависимость в виде: = если / =

,(<) + Ф(Г(/,) - х0 (/,)) + уы, если < / < 1 < / < ЛГ; (13) > (0 + Ф(Л*0 )(хл ('о) - хс(го))> если < Г < Г,.

Таким образом, построена параметрическая зависимость где q

- вектор уточняемых параметров, состоящий из компонент векторов

Критерием качества оценки является функционал, содержащий квадрат взвешенного отклонения априорно заданного вектора состояния от его расчётного значения, а также квадраты взвешенных невязок измерений и взвешенных модулей векторов суммарных возмущений между измерениями. Этот функционал можно представить в виде:

^ = (хл ('о,'ч) ■- ■хо )Т Ро'' (х А{(0, ч)- х0)] +

, ; т 04)

1 Л'-]

42>№У, .

1=0

Если минимум функционала (14) искать методом Ньютона, то поправки каждого шага итерации минимизируют квадратичную форму, полученную из (14) заменой нелинейных зависимостей линейными членами ряда Тейлора. Это означает, что на шаге итерации решается задача оптимальной оценки состояния линейной системы.

Матрицы А(,)(г) и Н;5' (г), ; = этой системы вычисляются по

следующим формулам:

1 ' Эх

-*,('.ч("4)' ' 5)1

(15)

4-Л

Обозначим

-^.(^хД^^)), / = 1,Ф, =Ф(/,+1,/,), / = 0,1,...,Л^—1.

Значения поправок 5х( / = 0,для которых достигается минимум квадратичной формы

2

дм

-)ТК, -НЙ5х,+и) + ад

/=0 ^ ^

(16)

при ограничениях:

5х,.+|1Л = Ф,(1)5х^ + ^Л,, Ф, = Ф(гм,*,),/ = 0Д...,#-1 (17)

удовлетворяют системе уравнений:

бх0Л = Р0(Ф«)\ + х0-х^>,

= (18) = (ФЙ)Т + Я,;1,{ъм -НЙ5х(+и), для / = 0,...,ЛГ -1,

Эта система уравнений распадается на две системы рекуррентных соотношений. Одна система связывает априорную оценку с конечной поправкой на момент последнего измерения. Происходит продвижение по измерениям вперёд. Эта рекуррентная система совпадает с алгоритмом фильтра Калмана. В результате применения первой системы рекуррентных соотношений получается поправка на момент последнего измерения. Однако чтобы вычислить остаточные невязки и перейти к следующему шагу итерационного процесса, нужно получить оценки векторов шума. Это обеспечивает другая рекуррентная система, которая продвигается от последнего измерения к первому и осуществляет сглаживание. Она позволяет определить оценки векторов состояния (поправок) на момент времени по информации на интервале от ^ до 4-, а также восстановить значения векторов шума.

Во втором разделе главы показано, что в случае отсутствия

немоделируемых ускорений алгоритм получения оценки совпадает с классическим алгоритмом5.

Если шум близок к постоянным систематическим воздействиям, целесообразно применять алгоритм, позволяющий оценивать средние значения приращений этих воздействий на мерном интервале. Алгоритм оценки вектора состояния и средних значений возмущений рассмотрен в третьем разделе главы.

В тех случаях, когда точность и состав измеряемых функций не позволяют оценить параметры шума, целесообразно применять метод мешающих параметров. Если в качестве оцениваемого вектора состояния выбрать вектор состояния на конец мерной базы, то неучтённый шум будет приводить к увеличивающимся ошибкам модели по мере перемещения от конца мерной базы к её началу. Суть метода мешающих параметров состоит в учёте этой нарастающей ошибки модели в весовой матрице измерений. Алгоритм оценки вектора состояния с использованием метода мешающих параметров рассмотрен в четвёртом разделе главы. Рассмотрены варианты алгоритма для двух типов возмущений: белого шума и случайных величин, постоянных на всем интервале. Алгоритм оценки вектора состояния с использованием метода мешающих параметров целесообразно применять при решении задач оценки точности определения параметров движения КА. Этот алгоритм позволяет оценить воздействие шума при приближённых представлениях о его статистических характеристиках.

В пятом разделе главы приведено обоснование алгоритма, использованного в первом разделе. Рассмотрен алгоритм оценки вектора состояния и возмущений дискретной динамической системы. Задача в такой постановке решается на каждом шаге итерационного процесса в алгоритме оценки вектора состояния и суммарных воздействий возмущений между

5 Аким Э.Л., Энеев Т.М. Определение параметров движения космического летательного аппарата по данным траекторных измерений // Космические исследования. 1963. Вып. 1. Т.1.

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

В главе 4 рассмотрены вопросы баллистики и навигации в проектах полётов к Луне, планетам и малым телам Солнечной системы, в разработке которых принимал участие автор диссертации. Цель главы состоит в изучении особенностей проектов с точки зрения баллистики и навигации и сведению частных задач проектов к некоторой обобщённой задаче, методы решения которой ([3], [5]) разработаны в интересах проекта «Фобос-Грунт».

Глава содержит анализ проблем баллистики и навигации в проектах полётов к Луне, и планетам Солнечной системы. Рассмотрены отечественные проекты «Фобос-Грунт», полёты к Луне, к Венере и Юпитеру, а также проект "8о1аг ОгЬкег".

Схема полёта КА «Фобос-Грунт» включает четыре последовательных

этапа:

- старт космического комплекса и выведение КА на отлётную от Земли траекторию;

- перелёт КА от Земли к Марсу, завершаемый переходом на орбиту искусственного спутника Марса (ИСМ);

- орбитальное движение КА в сфере действия Марса и маневрирование, обеспечивающее тесное сближение КА с Фобосом, посадку КА на его поверхность, забор образцов грунта, старт с Фобоса возвращаемого аппарата (ВА) и его полёт по околомарсианской орбите ожидания;

- старт с орбиты ожидания возвращаемого КА и перелёт к Земле,

завершаемый входом в её атмосферу и посадкой в заданный район территории нашей страны.

Первоначально проект «Фобос-Грунт» был ориентирован на применение ракеты-носителя (РН) среднего класса типа «Союз», исходя из чего были разработаны основные проектно-компоновочные решения по КА и соответствующая им схема полёта. Но в процессе дальнейшей работы возникло предложение - и оно было принято - доставить на начальную трёхсуточную орбиту ИСМ в качестве попутной полезной нагрузки китайский спутник «УН-1». Однако массовый резерв КА «Фобос-Грунт» при этом оказался недостаточным, и в связи с этим начата проработка запуска с помощью более мощной РН «Зенит-28ЬБ». Выбор схемы проведения коррекций, оценка затрат характеристической скорости и точность приведения КА к Марсу рассмотрены в пятой главе.

Решение задачи сближения с Фобосом начинается после перехода КА на орбиту наблюдения. Плоскость орбиты этой околокруговой орбиты совпадает с плоскостью орбиты Фобоса, а радиальное расстояние больше на 500 км. На этой орбите происходит уточнение орбитальных параметров движения Фобоса и его гравитационного параметра. При полёте по этой орбите окончательно выбираются дата и время посадки и тем самым параметры КСО, с которой будет выполняться сход и посадка ПМ на Фобос. После этого при необходимости выполняется манёвр «фазирования», который должен обеспечить оптимальные условия перехода с орбиты наблюдения на КСО. Выполненные расчёты и моделирование показали, что затраты характеристической скорости на переход с орбиты наблюдения на КСО не должны различаться более, чем на 5 м/с по сравнению с переходом Гомана между круговыми орбитами, радиусы которых равны радиусам орбиты наблюдения и орбиты Фобоса соответственно. Перед манёвром перехода на КСО начинается этап интенсивных траекторных измерений и определения орбиты КА и Фобоса. Выполненные в [3] и [5] оценки показали, что для этого

требуется мерная база протяжённостью 23 дня при условии, что возмущения, вносимые служебными системами ориентации и стабилизации в движение центра масс КА, на этом интервале будут сведены к минимуму. В этом случае полученная ошибка прогноза относительного движения на момент перехода с орбиты наблюдения на КСО не превосходит 1.5 км по относительному положению и 0.4 м/с по относительной скорости, ошибка определения гравитационного параметра Фобоса не превышает 10% от его величины. Двухимпульсным манёвром КА переходит с орбиты наблюдения на КСО. Ошибки исполнения каждого импульса этого двухимпульсного манёвра приводят, с одной стороны, к искажению параметров КСО, а с другой стороны, к увеличению ошибки прогноза относительного положения КА и Фобоса. При этом искажение параметров КСО со временем возрастает, а относительная ошибка уменьшается за счет выполнения новых бортовых и наземных траекторных измерений и определения параметров движения по ним. Искажение параметров КСО может быть частично компенсировано за счёт смещения во времени начала сеанса посадки относительно проектного (номинального) момента времени. Для определения и прогнозирования параметров относительного движения КА и Фобоса на КСО используется следующий состав навигационных измерений: доплеровские измерения с измерительных пунктов в Уссурийске и Медвежьих Озёрах (по одному сеансу каждого измерительного пункта на каждом из 10 витков); три сеанса бортовых телевизионных наблюдений Фобоса; три сеанса измерения расстояния с борга КА до Фобоса с помощью лазерного высотомера. Определение параметров относительного движения КА и Фобоса по указанному составу измерений на 4-суточном интервале позволяют прогнозировать относительное движение КА и Фобоса на сутки вперёд с ошибками 3 км по положению и 0.5 м/с по скорости. Эти орбитальные данные движения КА по КСО используются для расчёта параметров манёвра перехода КА с КСО на траекторию спуска и посадки на поверхность Фобоса. Результаты расчёта закладываются с Земли

на борт КА и реализуются в автономном режиме в сеансе спуска и посадки КА на поверхность Фобоса.

Старт с Фобоса представляет собой сложную динамическую операцию. Эскизным проектом предусмотрено, что он должен выполняться автоматически в заданное время. Это время задается перед началом сеанса посадки и может быть изменено в сеансе связи с КА, находящемся на поверхности Фобоса. При старте с Фобоса возвращаемый аппарат (ВА) не должен повредить станцию, остающуюся на поверхности Фобоса, должна быть исключена возможность столкновения с Фобосом. При выполнении таких динамических операций возрастает ошибка прогноза движения ВА. Поэтому должна быть предусмотрена возможность поиска ВА наземными станциями слежения после его старта с Фобоса.

После старта с Фобоса ВА выходит на орбиту, именуемую базовой. С учётом места посадки КА в качестве базовой орбиты выбрана круговая орбита, высота которой на 300-350 км ниже высоты орбиты Фобоса. К сожалению, ВА не может долго находиться на базовой орбите из-за затенений солнечных батарей. Поэтому В А будет переведён на промежуточную орбиту, на которой он должен безопасно существовать до начала манёвров для отлёта к Земле.

Основной задачей, решаемой на участке полёта возвращаемого аппарата к Земле, является приведение спускаемого аппарата на заданный полигон посадки на территории страны и надёжное обнаружение этого аппарата. Схема полёта предполагает проведение пяти коррекций траектории.

В главе рассмотрены также два проекта полёта к Луне, предусмотренные Федеральной космической программой России: проект «Луна-Ресурс» и проект «Луна-Глоб», которые предусматривают посадку на Луну. Для выполнения таких проектов необходимо выполнение следующих операций:

- приведение КА к Луне;

- торможение у Луны и переход КА на околокруговую полярную орбиту ИСЛ с высотой около 100 км над поверхностью;

- проведение на околокруговой орбите с высотой 100 км траекторных измерений, определение по ним параметров движения КА, расчёт и выполнение двух связанных коррекций, обеспечивающих выполнение посадки в заданную область;

- проведение траекторных измерений на околокруговой орбите ИСЛ после выполнения коррекций, определение параметров движения по этим измерениям, расчёт и выполнение орбитального манёвра, обеспечивающего переход на предпосадочную орбиту с высотами апоселения и периселения 100 и 18 км соответственно, с угловой дальностью 8° - 10е от периселения до точки посадки;

- проведение траекторных измерений на предпосадочной орбите ИСЛ, определение параметров движения и расчёт уставок на посадку;

- торможение и приведение КА в точку, расположенную примерно на высоте 2 км над точкой посадки, движение КА на участке прецизионного торможения, достижение поверхности Луны со скоростью, не превосходящей 3 м/с по вертикальной составляющей и 1 м/с — по горизонтальной.

Вопросы оценки точности приведения КА к Луне рассмотрены в пятой

главе.

Научные задачи другого проекта - «Венера-Д» предъявляют высокие требования к навигационному обеспечению. На этапе подлёта к Венере КА должен выполнить ряд динамических операций: выполнить коррекцию траектории, чтобы обеспечить точное наведение зондов, далее отделить их и выполнить манёвр перехода на пролётную траекторию с целью последующего торможения и перехода на орбиту искусственного спутника Венеры. При этом КА должен опередить зонды так, чтобы в момент их входа в атмосферу уже находиться на рабочей орбите и обеспечить приём информации от зондов. В

качестве рабочей орбиты КА предполагается полярная орбита с периодом около одних земных суток и высотой перицентра около 300 км.

Целью проекта полёта к Юпитеру «Лаплас» является посадка космического аппарата на поверхность естественного спутника Юпитера Европы. Перелёт к Юпитеру планируется осуществить с использованием ЭРДУ. Различные варианты параметров схемы перелёта были рассмотрены в работе6. Вопросы баллистики и навигации этого проекта рассмотрены в работе [18]. В пятой главе диссертации на примере этого проекта рассмотрены вопросы навигации и управления при постоянно работающей электроракетной двигательной установке.

На примере проекта «Solar Orbiter» в диссертации рассмотрены вопросы управления и навигации при проведении гравитационных манёвров.

Для получения требований по точности к управлению и навигации на участке гравитационного манёвра было выполнено моделирование полёта КА на участке от второго пролёта Земли (2018/08/08) до третьего пролёта Венеры (2018/10/09). В результате определены требуемые точности приведения КА к Венере, необходимые для успешного выполнения гравитационного манёвра при её втором пролёте 2018/10/09. Цель этого гравитационного манёвра состоит в увеличении наклонения КА по отношению к эклиптике и уменьшении периода до величины, которая составляет 2/3 от периода обращения Венеры вокруг Солнца. Такие элементы орбиты КА должны обеспечить гравитационный манёвр при следующем пролёте Венеры через три оборота КА вокруг Солнца. В результате гравитационного манёвра скорость КА должна быть уменьшена до расчётной величины. Для этого КА должен пройти через точку пересечения орбит Венеры и КА перед Венерой. Высота перицентра пролётной гиперболы должна составлять 1122 км. В результате гравитационного манёвра ICA сообщается импульс 4.148 км/с.

6 Федотов Г.Г. Об использовании гравитационного потенциала Земли для межпланетных полётов // Космические исследования, 2007, том 45, № 2, С. 169-179.

Проекция вектора импульса на направление скорости в гелиоцентрической СК составляет -3.903 км/с.

Из-за ошибок приведения КА к Венере изменение его скорости в результате гравитационного манёвра произойдет не на расчётную величину. Из-за этого не будут выполняться условия, необходимые для исполнения гравитационного манёвра при следующем пролёте. Поэтому ошибки исполнения гравитационного манёвра должны быть компенсированы коррекцией, выполняемой двигательной установкой КА. Величина характеристической скорости, необходимой для выполнения этой коррекции, зависит от ошибок исполнения гравитационного манёвра, которые, в свою очередь, зависят от ошибок приведения. Ошибки приведения рассматривались в картинной плоскости Венеры. Были выполнены расчёты разбросов дополнительных затрат характеристической скорости, необходимой для компенсации ошибок исполнения гравитационного манёвра, в зависимости от ошибок приведения. В случае если предельная ошибка приведения не превосходит 50 км, требуется не более 40 м/с дополнительной характеристической скорости на её исправление. Это составляет 1% от сообщаемого КА импульса и соответствует средним ошибкам исполнения импульсов, выполняемых двигателями КА. В случае если при выполнении гравитационного маневра предельные ошибки приведения в картинной плоскости не превосходят 10 км, для их компенсации потребуется не более 15 м/с характеристической скорости. В свою очередь это потребует точности навигации, при которой ошибки в картинной плоскости не превзойдут величину 5 км.

Конкретные задачи проектов, рассмотренных в этой главе, соответствуют следующей обобщённой задаче.

После выполнения динамических операций имеется область ошибок приведения КА к цели, например, приведение КА к планете или выход над заданной точкой спутника планеты в заданное время для начала сеанса

посадки. В связи с этим следует определить требования к средствам наземных траекторных измерений: точности, которые должны обеспечивать средства, и их территориальное размещение. Далее следует выбрать интервалы проведения навигационных сеансов, число коррекций, моменты проведения коррекций с учётом ошибок их исполнения и навигационных ошибок.

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

Например, увеличение интервала выполнения траекторных измерений позволяет уменьшить навигационную ошибку, но приводит к увеличению затрат характеристической скорости. Заблаговременное проведение коррекции при подлёте к планете-цели, позволяет сократить затраты характеристической скорости, но может привести к большим ошибкам прилёта, т. к. с течением времени возрастает отклонение от идеальной траектории, получаемое от ошибок исполнения коррекции. Следует также отметить эффект уменьшения навигационных ошибок по мере приближения к планете-цели, который обусловлен повышением информативности наземных траекторных измерений при достижении КА сферы действия планеты. Предлагаемый метод гарантированной оценки точности решения обобщённой баллистической задачи рассмотрен в главе 5.

Глава 5 посвящена методу анализа точности приведения КА к цели и оценки энергетических затрат. Рассмотрены вопросы навигации и управления при постоянно работающей ЭРДУ.

Требования по точности приведения определяются задачами проектов. Характерные примеры рассмотрены в главе 4. Фактическая точность приведения КА к планете-цели определяется ошибками прогнозирования параметров движения, ошибками исполнения манёвров или коррекций, а также неучтёнными возмущениями в движении центра масс КА.

Результаты, представленные в данной главе, получены на основе общего метода и описаны в разделах анализа точности определения и прогнозирования параметров движения КА в проекте «Фобос-Грунт» [3] и [5]. Этот метод применим для анализа точности определения и прогнозирования параметров движения КА в полётах к Луне, планетам и малым телам Солнечной системы.

Метод основан на комбинированном применении метода ковариационного анализа и метода статистических испытаний. Для оценки навигационных ошибок применяется ковариационный анализ. Влияние ошибок исполнения определяется методом статистических испытаний. Пусть заданы моменты проведения коррекций и исходное множество ошибок, например, после выхода на отлётную траекторию от Земли или после выхода на орбиту искусственного спутника Луны или планеты. Выполним расчёт навигационных ошибок на момент проведения первой коррекции с использованием опорной траектории. Свяжем каждую реализацию ошибки по целевым параметрам с отклонением вектора скорости от вектора скорости опорной траектории на момент задания исходного множества ошибок. Разобьём множество ошибок на ячейки по значениям ошибок целевых параметров. Каждую ячейку будем представлять одной или несколькими реализациями. Выполним прогноз движения КА на момент проведения первой коррекции по выбранным представителям каждой ячейки, далее выполним расчёт коррекции и исказим его навигационными ошибками и ошибками исполнения. В результате получим множество ошибок на момент проведения первой коррекции и максимальные затраты характеристической скорости на её проведение. Далее переходим ко второй коррекции и последующим коррекциям. После выполнения расчётов по последней коррекции получаем суммарные максимальные затраты характеристической скорости и оценку ошибок приведения. Выполняя расчёты по описанному выше методу, получаем характеристики рассматриваемой схемы проведения коррекций:

точность приведения и затраты характеристической скорости. Далее, варьируя числом коррекций и моментами их исполнения, получаем приемлемую схему.

Для применения метода необходимо иметь функциональную связь ошибок по целевым параметрам с отклонением вектора скорости от его значения по опорной траектории. Например, при решении задачи приведения КА к Луне или планете - эта функциональная связь определяется зависимостью значений параметров картинной плоскости от компонент вектора скорости. При решении задач приведения КА на гало-орбиту вокруг точки ¿2 эта зависимость определяется зависимостью от вектора скорости параметров, определяющих геометрию гало-орбиты. Далее работа метода показана на задачах приведения КА к Луне или планете.

При построении моделей и выполнении расчётов предполагалось, что для определения параметров движения в случае КА с ЖРД используется классический метод, описанный в работе5, а в случае КА с ЭРДУ - методы, описанные в третьей главе. Измеряемыми параметрами являются наклонная дальность и радиальная скорость.

При условии идеального исполнения управляющих импульсов имеются следующие источники ошибок прогноза параметров движения КА:

- случайная составляющая ошибок измерений радиальной скорости и наклонной дальности;

- систематические смещения измеренных величин радиальной скорости и наклонной дальности за время сеанса измерений;

- неопределённость знания положения наземного измерительного пункта;

- неопределённость знания параметров ионосферы и тропосферы в принятых моделях измерений для случаев работы с учётом и без учёта калибровки по сигналам навигационных спутников;

- неопределённость коэффициента светового давления.

Модели движения КА с ЖРД и ЭРДУ отличаются моделями шума, воздействующего на движение КА. В случае ЖРД - этот шум обусловлен

наличием небольших по величине (1(Г9 - 10~7 мм/с2) негравитационных ускорений, вызванных неконтролируемыми реактивными силами. В случае ЭРДУ - это ошибки по направлению (Г) и величине вектора тяги (5% от величины).

Ковариационная матрица ошибок определения вектора состояния вычисляется по формуле:

К(0 = К„ (0 + ДКЬ (0 + ЛК„(г) + ДК1(г) + ДК,(/), (19)

где

к0(г) - ковариационная матрица ошибок, обусловленная неучтёнными ускорениями и случайными составляющими ошибок измерений;

ДКЬ(/) - поправка на ошибки, обусловленные систематическими составляющими ошибок измерений;

ДКа(г) - поправка на ошибки, обусловленные ошибками привязки фазового центра антенны;

ДК,^) - поправка на ошибки, обусловленные тропосферной составляющей ошибок измерений;

ДК,(7) - поправка на ошибки, обусловленные ионосферной составляющей ошибок измерений.

Для расчёта ковариационной матрицы К0(?) используется метод мешающих параметров, который рассмотрен в третьей главе диссертации. В качестве оцениваемого вектора состояния выбран вектор состояния на конец мерной базы. Неучтённый шум приводит к увеличению ошибок модели по мере перемещения от конца мерной базы к ее началу. Нарастающие ошибки модели учитываются в весовой матрице измерений.

Движение КА в условиях неопределённости возмущающих ускорений описывается рассмотренной в третьей главе моделью (8), (9).

Матрица поправок ДКЬ (7Л,) на ошибки, обусловленные систематическими составляющими ошибок измерений, вычисляется как разность между ковариационной матрицей ошибок оценки вектора состояния, обусловленных систематической и случайной составляющими ошибок

измерений, и ковариационной матрицей ошибок оценки вектора состояния, обусловленных только случайными составляющими ошибок измерений.

Ковариационная матрица оценки вектора состояния КА, обусловленная систематическими и случайными ошибками измерений, вычисляется с использованием расширенного вектора состояния. Этот вектор наряду с вектором состояния КА содержит систематические составляющие ошибок измерений. При расчёте ковариационной матрицы ошибок определения расширенного вектора состояния используется априорная информация о СКО систематических ошибок измерений.

Аналогично вычисляются матрицы ЛКа(/), ЛКД?) и ДК.; (г). Специфика состоит в функциональной зависимости измеренных значений от компонент расширенного вектора состояния, а также в представлении априорной информации.

При выборе схемы проведения коррекций рассматривались ошибки (невязки между фактическими и заданными значениями) по положению в картинной плоскости (Д^, Лг[) и ошибка по времени приведения (Дт). Направления осей и г| картинной плоскости определяются единичными векторами е^и е^:

е х V V

К XV.! |У.|

еп - единичный вектор, направленный от центра планеты (или Луны) к её северному полюсу.

Сопоставление схем маневрирования между собой проводится по соответствующим им максимально возможным затратам характеристической скорости и максимальным ошибкам приведения КА к планете.

Оценка максимальных затрат характеристической скорости и максимальных ошибок приведения выполняется следующим методом. Трёхмерное пространство ошибок: Ат] и Дт разделяется на ячейки -

прямоугольные параллелепипеды в соответствии с заданным шагом.

Множество возможных ошибок после каждой коррекции покрывалось набором этих ячеек. Если в некоторую ячейку попала хотя бы одна точка, то при оценке затрат на коррекцию нужно рассматривать все точки этой ячейки. Если на этапе коррекции к нужно устранить ошибки, оцениваемые

вектором Г = Дт), Дт)Т, то поставим в соответствие этому вектору множество Уг в пространстве ошибок. Это множество получается после расчёта корректирующего импульса и формирования множества ошибок, обусловленных ошибками исполнения и ошибками прогноза.

Обозначим как Ек результат построения покрытия прямоугольными параллелепипедами множества ошибок до проведения к-й коррекции.

Ек=[}иг (21)

где

£/, - прямоугольный параллелепипед с индексом_/';

- множество индексов, соответствующих параллелепипедам, в которых содержатся ошибки.

Для каждого и), у е Л построим множество Жк .:

(22)

Множество РГк . состоит из прямоугольных параллелепипедов, в которые попали ошибки, полученные при исправлении ошибок {т,ш = 1,...М. Эти ошибки обусловлены ошибками исполнения и ошибками прогноза. В качестве векторов исправляемых ошибок (т,т = 1,...М можно выбрать вершины параллелепипеда £/ и его среднюю точку.

Рассмотрим объединение множеств Шкдля индексов у е Ук. Построим покрытие этого множества системой параллелепипедов и обозначим

множество индексов параллелепипедов, которые вошли в покрытие, как Jм. В результате получим множество Еы :

¿V, = ■ (23)

1

Оценка максимальных затрат характеристической скорости получается автоматически, если при построении У^ определить максимальную характеристическую скорость, а затем выбрать максимальное значение по всем параллелепипедам ир

В последнем разделе главы рассмотрены три примера расчёта точности приведения КА, оснащённых жидкостными реактивными двигателями: к Марсу (на примере проекта «Фобос-Грунт»), к Луне (на примере проекта «Луна-Глоб») и к Венере (на примере проекта «Венера-Д») для КА, а также пример расчёта точности приведения к Юпитеру КА, оснащенного ЭРДУ.

3. Заключение

Основные результаты работы были перечислены в начале автореферата. В заключение приведённого изложения содержания диссертации следует упомянуть, что оно является попыткой обобщения многолетнего опыта практической работы в трёх направлениях прикладной небесной баллистики:

- в выполнении проектных баллистических исследований, позволяющих определять облик космического объекта, сроки его запуска и параметры баллистического обеспечения будущего полёта;

- в определении параметров траектории полёта КА после его успешного запуска и в прогнозировании будущего движения с гарантированной оценкой точности получаемого прогноза;

- в построении манёвров управления траекторией полёта на всех её этапах.

4. Основные публикации по теме диссертации

(Публикации 1-7 из основного перечня ВАК)

1. Тучин А.Г. Квазисинхронные орбиты и их использование для сближения космического аппарата с Фобосом. // Космические исследования. Т. 45. № 2, 2007. - С.144-149.

2. Тучин А.Г. Проектирование квазисинхронных орбит космического аппарата вокруг Фобоса для решения задачи посадки на его поверхность. Н Космические исследования. Т. 46, № 6, 2008. -С.536-546.

3. Аким Э.Л., Заславский Г.С., Морской И.М., Рузский Е.Г., СтепаньянцВ.А., Тучин А.Г. Доставка на землю реликтового вещества с Фобоса - проект «Фобос-Грунт». Баллистика, навигация и управление полётом. // Астрономический вестник, 2010, том 44, № 1, С.29-40.

4. Аким ЭЛ., Архангельский Р.Н., ЗайкоЮ.К., Лавренов С.М., Порошин А.Л., Е. Г. Рузский, В. А. Степаньянц, Тучин А.Г., Тучин Д.А., Федотов В.П., Ярошевский B.C. Баллистика, навигация и управление движением КА на этапе его посадки на поверхность Фобоса. // Космические исследования, 2009, Т. 47, № 4, С.332-342.

5. Аким Э.Л., Заславский Г.С., Морской И.М., Степаньянц В.А., Тучин А.Г. Баллистика, навигация и управление полётом космического аппарата в проекте «Фобос-Грунт» // Известия РАН. Теория и системы управления. 2002. №5, С. 153-161.

6. Akim E.L., Stcpaniants V.A., Vlasova Z.P., Tuchin A.G. Orbit determination for a radar mapping satellite of Venus. Earth, Moon, and Planets (ISSN 0167-9295), vol. 50-51, July-Dec. 1990, p. 559-578.

7. Аким Э.Л., Астахов А.П., Бакитько P.B., Полыциков В.П., Степаньянц В.А., Тучин А.Г., Тучин Д.А., Федотов В.П., Ярошевский B.C. Автономная навигационная система околоземного космического аппарата. // Известия РАН. Теория и системы управления. 2009, №2. - С.139-156.

8. Аким Э.Л., Тучин А.Г. Проект «Фобос-Грунт» и прикладная математика // Будущее прикладной математики: Лекции для молодых исследователей. Поиски и открытия / Под ред. Г.Г. Малинецкого. - М.: Книжный дом «Либроком», 2009, С.363 -390. ISBN 978-5-397-00638-5.

9. Akim E.L., Popov G.A., Tuchin A.G. Mechanics and Motion Control of a Space Vehicle in the Project of Relict Substance Delivery on Earth (The Project "Phobos-Grunt"). // Proc. 16th IF AC Symposium on automatic control in Aerospace. 2004. Preprints. Saint Petersburg. Russia. V. 1.

10. Akim E.L., BotkinA.V, Stepaniants V.A., ShishovV.A., RuzskiyE.G., Tuchin A.G. Orbit Selection, Navigation and Maneuvers before the Landing on the Phobos Surface for Phobos Sample Return Project. // Proc. I7lh ISSFD. 16-20 June 2003. Moscow. Russia. V. 1.

11. AkimE.L., Stepaniants V. A., ShishovV.A., RuzskiyE.G., Tuchin A.G. Ballistics, navigation and motion control of the SC on stages of the Phobos surface

approaching and landing. // Proc. 18th ISSFD. 11-15 October 2004. Munich. Germany.

12. Тучин А.Г. Определение параметров движения КА по результатам измерений при наличии шума в динамической системе. // Препринт ИПМ им. М.В. Келдыша РАН, 2004, № 2.

13. Тучин А.Г. Проектирование квазисинхронных орбит КА вокруг Фобоса для решения задачи посадки на его поверхность. // Препринт ИПМ им. М.В. Келдыша РАН, 2008, №15.

14. Аким Э.Л., Горохова A.A., Киселева И.П., Степаньянц В.А., Тучин А.Г. Локальная обработка измерений радиосистемы межпланетных космических аппаратов. // Препринт ИПМ им. М.В. Келдыша РАН, 2002, №11.

15. АкимЭ.Л., Горохова A.A., Киселева И.П., Степаньянц В.А., Тучин А.Г. Небесно-механическая интерпретация измерений радиосистемы межпланетных космических аппаратов «Квант-Д». // Препринт ИПМ им. М.В. Келдыша РАН, 2002, № 2.

16. Квашнин А.Г., Тучин А.Г. Синтез вычислительной схемы алгоритма оценивания параметров динамической системы. // Препринт ИПМ им. М.В. Келдыша РАН, 1991, № 26.

17. Квашнин А.Г., Тучин А.Г. Построение вычислительной схемы алгоритма оценивания параметров динамической системы. // Автоматизация проектирования, информатика. - 1993. - Вып. 4. - С.88-94.

18. АкимЭ.Л., Сазонов В.В., Степаньянц В. А., Тучин А.Г. Управление и навигация в перелёте КА с малой тягой от Земли к Юпитеру. // Препринт ИПМ им. М.В. Келдыша РАН, 2010, № 36.

Подписано в почать 17.11.2010г. Заказ № 94. Формат бумаги 60X90 1/16. Тираж 70 экз.

Отпечатано в Институте прикладной математики РАН Москва, Миусская пл. 4

Обозначения и сокращения.

Введение

Глава 1 Проектирование квазисинхронных орбит КА вокруг

Фобоса для решения задачи посадки на его поверхность

1.1 Численный алгоритм построения множества КСО с минимальным дрейфом

1.1.1 Формализация постановки задачи

1.1.2 Упрощённая модель движения КА относительно Фобоса.

1.1.3 Область перебора

1.1.4 Переменные перебора.

1.1.5 Численный анализ

1.1.6 Результаты численного анализа.

1.1.7 Сравнение характеристик движения в полной и упрощённой моделях.

1.2 Начальное приближение для расчёта множества КСО

1.2.1 Упрощённые уравнения движения.

1.2.2 Предварительный численный анализ. Постановка задачи аналитических исследований

1.2.3 Подготовка системы (1.19-1), (1.19-2) к усреднению?.

1.2.4 Частный случай плоского кругового движения задачи Хилла (р = 1).

1.2.5 Вычисление значения интеграла Г— о ' д1 д

1.2.6 Вычисление производных —, —, —. сЫ да д^

1.2.7 Усреднённые уравнения движения.

1.2.8 Условие для поиска начального приближения начальных условий.

Выбор параметра К

1.2.9 Алгоритм формирования начального приближения.

1.2.10 Вычислительная процедура поиска минимума.

1.2.11 Результаты расчётов

1.2.12 Алгоритм проектирования КСО.

Глава 2 Баллистика, навигация и управление движением КА на этапе его посадки на поверхность Фобоса.

2.1 Сближение с Фобосом и посадка на его поверхность: общая схема.

2.1.1 Сближение с Фобосом

2.1.2 Условия посадки.

2.1.3 Навигационные приборы обеспечения посадки.

2.1.4 Схема управляемой посадки.

2.1.5 Управление включением двигателей.

2.2 Алгоритмы управления движением КА.

2.2.1 Бортовые алгоритмы навигации и управления.

2.2.2 Уравнения движения КА относительно Фобоса.

2.2.3 Расчёт коррекций на участке перелёта от момента схода с КСО до точки начала вертикального спуска.

2.2.4 Определение вектора состояния КА по измерениям лазерного высотомера и доплеровской системы.

2.3 Реализация алгоритмов посадки в среде операционной системы реального времени.

Глава 3 Определение параметров движения КА по результатам измерений при наличии немоделируемых ускорений.

3.1 Алгоритм оценки вектора состояния и суммарных воздействий возмущений между измерениями.

3.1.1 Постановка задачи.

3.1.2 Линейный случай

3.1.3 Нелинейный случай.

3.1.4 Проверка качества измерений с использованием приведённого среднеквадратического отклонения

3.2 Алгоритм оценки вектора состояния в случае отсутствия шума.

3.3 Алгоритм оценки вектора состояния и средних значений приращений возмущений

3.4 Алгоритм оценки вектора состояния с использованием метода мешающих параметров

3.4.Г Мешающие параметрь1 в форме белого шума

3.4.2 Мешающие параметры в форме случайных величин, постоянных на всем интервале.л.

3.5 Оценка вектора состояния и возмущений дискретной динамической системы и свойства этих оценок.

Глава 4 Вопросы баллистики и навигации в проектах полётов к Дуне, планетам и малым телам Солнечной системы

4.1 Проблемы баллистикй и навигации в проектах полётов к Луне, точке Ь2 системы Солнце - Земля, планетам и малым телам

Солнечной системы.

4.1.1 Проект «Фобос-Грунт»

4.1.2 Российские проекты полёта к Луне

4.1.3 Полёт к точке ¿2 системы Солнце — Земля.

4.1.4 Проект полёта к Венере с целью доставки на её поверхность долгоживущей станции, баллонов в атмосферу Венеры н выхода орбитальной станции на орбиту искусственного спутника Венеры.

4.1.5 Проект полёта к Юпитеру с целью посадки КА на.поверхность Европы, естественного спутника Юпитера (проект «Лаплас»).

4.1.6 Особенности баллйстико-навигационного обеспечения полётов к дальним планетам.

4.1.7 Подготовка и проведение гравитационного маневра КА «Розетта».

4.1.8 Проект полёта к Солнцу.

4.2 Методы обеспечения навигации и управления КА дальнего космоса.

4.2.1 Управление и навигация на участках перелёта к планетам с использованием ЖРД

4.2.2 Управление и навигация на участках перелёта к планетам с использованием ЭРДУ.

4.2.3 Управление и навигация при проведении гравитационных маневров .148'

4.2.4 Автономная навигация по изображениям планет на участках перелёта от Земли к планете-цели

4.2.5 Управление и навигация на орбитах искусственных спутников планет.

4.3 Обобщённая баллистическая задача.

Глава 5 Анализ точности приведения КА-к Луне и планетам.

5.1 Методика и вычислительный алгоритм оценки точности определения и прогнозирования параметров движения КА по наземным траекторным измерениям.

5.1.1 Расчёт ковариационной матрицы ошибок определения вектора состояния, обусловленных неучтёнными ускорениями и случайными составляющими ошибок измерений.

5.1.2 Поправки на ошибки, обусловленные систематическими составляющими ошибок измерений.

5.1.3 Поправки на ошибки,.обусловленные ошибками привязки фазового t центра антенны

5.1.4 Поправки на ошибки, обусловленные тропосферными составляющими ошибок измерений.

5.1.5 Поправки на ошибки^-обусловленные ионосферными составляющими ошибок измерений.

5.1.6 Прогнозирование ковариационной матрицы.

5.2 Методика выбора схемы-проведения коррекций.при перелётах к Луне и планетам.

5.3 Примеры расчётов определения точности приведения КА к Луне и планетам с использованием жидкостных реактивных двигателей.

5.3.1 Определение точности приведения КА к Марсу (на примере проекта «Фобос-Грунт»)

5.3.2 Определение точности приведения КА к Луне (на примере проекта «Луна-Глёб»).

5.3.3 Определение точности приведения КА к Венере (на примере проекта «Венера-Д»).

5.4 Определение точности приведения КА с ЭРДУ к Юпитеру в проекте «Лаплас».

5.4.1 Моделирование ойределения орбиты.

5.4.2 Модель ошибок ЭРДУ

5.4.3 Расширенная система.

5.4.4 Оценка ошибок определения вектора состояния и параметров линейных функций

5.4.5 Технология управления КА на участке перелёта.

5.4.6 Расчёт ошибок приведения КА к Юпитеру и оценки-дополнительного расхода топлива.

5.5 Требования по точности наземных и бортовых траекторных измерений.

5.5.1 Требования к точности наземных траекторных измерений

5.5.2 Требования к точности бортовых акселерометров

5.5.3 Требования к точности бортовых измерений направления на планеты, их спутники и астероиды.

Диссертация посвящена теоретико-механическим вопросам проектирования полетов космических аппаратов (КА) к планетам Солнечной системы и их естественным спутникам. Разработан метод -проектирования квазисинхронных орбит КА вокруг Фобоса, предназначенных для обеспечения посадки КА на его поверхность. Разработаны .методы, определения параметров движения КА по траекторным измерениям на фоне работы двигательной установки; разработаны методы расчета точности приведения КА к планете цели при использовании химических и электроракетных двигателей с учетом ошибок прогноза движения КА, ошибок исполнения маневров и коррекций и различных возмущений, вносимых в движение центра масс КА. разработанные методы опираются на опыт баллистико-навигационного обеспечения полетов к Венере КА «Венера-15,16», «Вега-1,2», к Марсу КА «Фобос-2» и >;нашли применение в проекте «Фобос-Грунт», при проектировании полетов к Луне в проектах «Луна-Ресурс» и «Луна-Епоб», к Венере (проект «Венера-Д»), в систему Юпитера с целью посадки КА на его естественный спутник Европу. Это определяет актуальность и практическую значимость диссертации.

Цель работы состоит в разработке теоретико-механических и математических методов, обеспечивающих баллистику и навигацию, в проекте «Фобос-Грунт» и распространению этих методов на решение задачь навигации и управления в перспективных отечественных проектах полётов к Луне (проекты «Луна-Ресурс» и «Луна-Глоб»), к Венере (проект «Венера-Д»), и в систему Юпитера (проект «Лаплас»).

Результаты, представленные в диссертации, докладывались автором на 16-м симпозиуме ШАС по автоматическому управлению (Санкт-Петербург, июнь 2004 г); на 17-м (Москва, июнь 2003г) и 18-м (Мюнхен, октябрь 2004г.) Международных симпозиумах по динамике космического полёта, на Общероссийском семинаре «Современные методы навигации и управления движением» (Институт проблем Управления РАН, 31 марта 2009 г.), на семинаре «Солнечная система и смежные проблемы физики и механики» (ИПМ им. М.В. Келдыша РАН, 18 марта 2008 г.).

10

Список основных публикаций пр теме--диссертации состоит из 18 работ; 7 работ из этого списка опубликованы в изданиях, рекомендованных.ВАК; 5 из них в соавторстве. Из работ, выполненных с соавторами, в диссертацию включены только результаты, полученные автором.

Основные результатыгполученные в диссертации, приведены.в,заключении.

В первой главе рассмотрены методы проектирования квазисинхронных орбит КА вокруг Фобоса для решения задачи посадки на его поверхность. Результаты, представленные в этой главе, описаны в работах: [112-ТисЫп, 2003; 57-Тучин, 2007; 60 - Тучин, 2008; 61 - Тучин, 2009; 59 - Тучин, 2008]

Одним из основных этапов проекта «Фобос-Грунт» является этап непосредственного сближения КА с Фобосом. На этом, этапе должны быть получены навигационные данные, необходимые для успешной, работы автономной системы посадки. Для решения этой задачи необходима такая орбита КА, двигаясь по которой КА не удалялся-бь^более, чем на 70 км. от Фобоса. К классу таких орбит в терминологии ограниченйой задачи трёх тел относятся- квазиспутниковые орбиты [37- Лидов, 1994]. Это орбиты с обратным движением, охватывающие тело меньшей массы и расположенные вне сферы его действия. Ниже будем называть такие орбиты квазисинхронными. Двигаясь по; такой орбите, КА облетает Фобос за интервал времени, который примерно равен: периоду обращения Фобоса вокруг Марса.

Моделирование показало, что требуются четыре дня полёта КА .на КСО для определения параметров движения с необходимой точностью [67 - Тучин, 2002; 66-Тучин, 2010]. Поэтому интервал пребывания на КСО оценивается пятью сутками. К номинальной КСО предъявляется ряд .требований. За время пребывания на КСО с борта КА должно быть получено телевизионное изображение предполагаемого района посадки и проверена работоспособность лазерного высотомера автономной системы посадки. В момент начала сеанса посадки, при получении изображения района посадки и проверки работоспособности лазерного высотомера, удаление от поверхности Фобоса не-^должно превосходить 60 км и должны выполняться условия по освещённости и радиосвязи со станциями слежения в Евпатории и Уссурийске. Если в штатном сеансе работы автономной системы посадка не будет осуществлена и КА останется на прежней орбите, должна быть обеспечена возможность повторения попытки посадки на Фобос.

Выбор КСО предполагается проводить по следующей схеме.''Из условий* освещенности, обеспечения-радиосвязи со станциями слежения; возможно, других условий выбирается точка над поверхностью Фобоса и момент времени, в который КСО должна пройти через эту точку. В первой главе диссертации разработан метод, позволяющий найти такую квазисинхронную орбиту, проходящюю через заданную точку в заданное время, которая при других прохождениях долготы выбранной точки удаляется от неё на минимальное расстояние. Это свойство обеспечивает выполнение условия по максимальному удалению КА от поверхности Фобоса при каждом прохождении над районом посадки. Тем самым возможность телевизионной съёмки предполагаемого участка посадки и условия проверки работоспособности лазерного высотомера будут обусловлены только условиями ^.освещённости и радиосвязи со станциями слежения.

Вторая глава посвящена алгоритмам навигации и управления в схеме посадки на Фобос. Результаты, представленные* в этой главе, описаны в работах [114-Tuchin, 2004; 63-Тучин, 2009]. В этой главе.содержится модель движения КА относительно Фобоса, которая достаточно проста для бортовых вычислений, обеспечивает требуемую точность и позволяет.-рассчитывать, коррекции и решать навигационные задачи. Далее строится алгоритм оценки вектора состояния КА по измерениям лазерного высотомера и доплеровского измерителя скорости и дальности с использованием цифровой модели поверхности Фобоса, разработанной специалистами Института геохимии и аналитической химии им. В.И. Вернадского РАН. В этой главе анализируется схема посадки и условия ее выполнения.

В третьей главе рассмотрены вопросы определения параметров движения КА по результатам траекторных измерений при наличии немоделируемых действующих ускорений.

Задача определения параметров движения космического аппарата (КА) является одной из основных задач, решаемых в ходе управления его полётом. При решении этой задачи часто возникает ситуация,- в которой определение параметров движения КА надо выполнять на фоне работы двигателей. В качестве примера можно привести следующие задачи: контроль-)£частка, выведения КА на орбиту искусственного спутника; оперативная, оценка исполнения импульсов* по измерениям, наземных средств на фоне ^работы двигательной' установки; определение параметров движения КА с электроракетной двигательной установкой11

ОРДУ).

К вопросам; которые рассматриваются в диссертации, в, первую очередь, относятся проблемы решения задач навигации исправления полётом с включённым ЭРДУ. Фактическое ускорение, создаваемое ЭРДУ, отличается от модели этого5 ускорения, заложенного' в расчёты. Имеются ошибки, величины и ориентации» вектора тяги ЭРДУ в пространстве.

Решение указанных, выше задач основано на< применении« моделей динамических систем, в- которых помехи имеются^не только в измерениях, но и влияют на поведение самого объекта. Такие модели обычно исследуются в рамках линейных моделей в общей теории систем. Применение этих методов в задачах определения движения КА требует развития соответствующих нелинейных моделей и учёта особенностей уравнений динамики, и измеряемых данных.

Результаты, представленные в данной1 главе, описаны »в работе [62-Тучин, 2009] в части алгоритмов определения'параметров-движения КА и в работе [58-Тучин, 2004].

При решении указанных выше задач применяются различные модели шума и, соответственно, используются/ различные- методы и алгоритмы- оценки вектора* состояниям Могут применяться комбинированные методы осреднения:

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

- метод наименьших квадратов;

- расширенный фильтр Калмана.

Многолетняя практика показала, что метод-наименьших'квадратов является очень надёжным методом определения параметров. Оценка получается в результате поиска минимума функционала, представляющего собрй сумму взвешенных невязок между измеренными-значениями и их расчётными аналогами. При этом расчётные аналоги функционально зависят от уточняемых-параметров. Особенностью метода является« то, что модель, движения КА должна ;.быть достаточно^ точной. Не допускается^ наличие "больших возмущений, которые не задаются- в. виде зависимостей от уточняемых параметров. Поэтому применение метода наименьших квадратов в случае неизвестных интервалов! работы двигателя« малой тяги в данном случае вызывает определённые проблемы. Метод,,наименьших квадратов может быть успешно использован на этапе, когда уже имеется оценка моментов включения и выключения двигателями создаваемого им ускорения.

Расширенный фильтр Калмана является эмпирическим расширением фильтра Калмана для линейных систем на -нелинейный ^случай. Расширенный фильтр Калмана предполагает наличие неизвестного фазового шума, воздействующего на систему. Характеристики шума задаются, его ковариационной- матрицей. Особенностью этош метода является то, что текущие значения оцениваемых параметров- должны 1 находиться в^достаточно близкой окрестности относительно их истинных значений. Кроме того, метод требует, чтобы не было длительных интервалов времени, в которых нет измерений.

В1 диссертации для оценки вектора состояния в условия воздействия* шума разработан метод, который обеспечивает оценку вектора состояния и суммарные возмущения между измерениями. Предполагается, что фазовый шум, воздействующий на систему, является белым-или-может быть получен из белого шума формирующими фильтрами. Пример использования формирующих фильтров для формирования шума; вызванного ошибками- исполнения программы ЭРДУ, рассмотрен в пятой главе.

Суть предложенного метода состоит в-том, что оценка - определяется из условия минимизации функционала, зависящего как - от невязок измеренных значений и их расчётных.аналогов, так и от величин определяемых возмущений. Минимизация функционала выполняется итерационно. На каждом шаге итерационного процесса определяется поправка к искомым параметрам.

Определение поправки производится из условия* минимума функционала для линейной системы. Поиск минимума функционала для линейной системы приводит к двум* последовательностям рекуррентных" формул. Первая последовательность И рекуррентных формул идёт от первого измерения, к последнему измерению и позволяет определить? поправку на момент последнего измерения. Вторая-последовательность рекуррентных формул идёт от последнего значения-к первому и позволяет восстановить, возмущения: Первая^ последовательность рекуррентных формул эквивалентна- рёкуррентным -формулам фильтра: Калмана для^ линейной' системы.

В четвёртой главе рассмотрены^ вопросы баллистики и навигации в проектах полётов.к Луне, планетам и-малым телам Солнечной системы, в которых принимал участие автор1диссертации. Цель главы состоит в изучении особенностей проектов с точки зрения баллистики и навигации и сведению частных задач проектов к некоторой обобщённой задаче.

Четвертая глава'содержит анализ пр о б л ем^балл и сти ки и навигации в проектах полётов к Луне, точке Ь2 системы Солнце - Земля и планетам Солнечной системы. Рассмотрены, отечественные проекты «Фобос-Грунт», полёты к Луне, в окрестность точки Ь2 системы .Солнце - Земля, к Венере и Юпитеру, а .также проекты «Кассини», «New Horizons» и «Пионер-10, -11» (США).

Конкретные задачи проектов, рассмотренных в четвёртой' главе, соответствуют следующей обобщённой задаче.

После- выполнения динамических операций имеется область ошибок приведения КАк цели, например, приведение КА к планете или выход над заданной точкой* спутника планеты в, заданное* время для начала сеанса* посадки. В1'связи* с этим следует определить требования к .средствам наземных траекторных измерений: их точности и их размещение на оси времени полёта. Далее следует выбрать интервалы проведения навигационных сеансов, число коррекций, моменты проведения коррекций с учётом ошибок их исполнения и навигационных ошибок.

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

Например, увеличение интервала выполнения траекторных измерений позволяет уменьшить навигационную ошибку,^ но приводит к увеличению затрат характеристической скорости. Заблаговременное, проведение коррекции при подлёте к планете-цели, позволяет сократить затраты характеристической скорости, но может привести к большим ошибкам прилёта, т. к. с течением времени возрастает отклонение от идеальной траектории, получаемое от ошибок исполнения коррекции. Следует также отметить эффект уменьшения навигационных ошибок по мере приближения к планете-цели, который обусловлен повышением информативности наземных траекторных „измерений при достижении КА сферы действия планеты. ?

Предлагаемый метод гарантированной оценки точности решения обобщённой баллистической задачи рассмотрен в пятой главе.

Результаты, представленные в пятой главе, описаны в разделах анализа точности определения и прогнозирования параметров движения КА в проекте «Фобос-Грунт» [61 -.Тучин, 2009; -66 - Тучин, 2010; 67 - Тучин, 2002, 113 — ТисЫп, 2004]. Эти методы применимы для ^анализа точности определения и прогнозирования параметров' движения КА в полётах к Луне, планетам и малым телам Солнечной системы и были использованы ".при подготовке научно-технических отчётов [122 - НТО'5-14-04; 124 ^НТО £-04:06; 132 -НТО 5-31-07] по проекту «Фобос-Грунт»; [128-НТО 5-06-08; 135 - НТО 5-02-093 по проектам полёта к Луне; [141 - НТО 5-013-09] по проекту «Венера-Д».

Фактическая точность- приведения КА к-планере-цели-определяется ошибками прогнозирования параметров движения; ошибками исполнения маневров или коррекций, а также неучтёнными возмущениями "-выдвижении центра масс КА.

Предложенный в пятой главе методгарантированной оценки точности решения обобщённой баллистической задачи основан, на комбинированном применении метода ковариационного анализа^метода статистических испытаний. Для оценки навигационных, ошибок применяется ковариационный анализ. Влияние ошибок исполнения определяется методом статистических испытаний. Пусть заданы моменты проведения коррекций и исходное множество ошибок, например, после выхода на отлётную траекторию от Земли или после выхода на орбиту искусственного спутника Луны или планеты. „Выполним расчет навигационных ошибок, на момент проведения- первой коррекции с .использованием опорной траектории. Свяжем каждую реализацию ошибки по целевым параметрам с отклонением вектора скорости от вектора скорости опорной траектории на момент задания исходного множества ошибок., Разобьём множество ошибок на ячейки по значениям ошибок целевых параметров*. Каждую ячейку будем представлять одной или несколькими реализациями. Выполним прогноз движения. КА на момент проведения первой коррекциишо выбранным представителям каждой ячейки, далее выполним расчет коррекции и исказим его навигационными ошибками и ошибками исполнения. В результате получено множество ошибок на момент проведения первой коррекции и максимальные затраты характеристической скорости на её проведение. Далее переходим ко- второй коррекции, и последующим коррекциям. После выполнения расчётов по последней коррекции получаем суммарные максимальные затраты характеристической скорости и оценку ошибок приведения. Выполняя расчёты по описанному выше методу, получаем характеристики рассматриваемой схемы проведения коррекций: точность приведения и затраты характеристической скорости. Далее, варьируя-числом коррекций-и моментами их исполнения, получаем приемлемую схему.

Для применения метода необходимо иметь функциональную связь ошибок по целевым параметрам с отклонением вектора скорости от его значения по опорной-траектории. Например, при решении задачи приведения КА к Луне или планете -эта функциональная' связь определяется зависимостью значений параметров картинной плоскости от компонент вектора скорости. При решении задач приведения КА на гало-орбиту вокруг точки Ь2 эта зависимость определяется зависимостью от вектора скорости параметров, определяющих геометрию гало-орбиты. Работа метода показана на задачах приведения КА к Луне или планете.

Основные результаты диссертации состоят в следующем:

1. Разработан метод проектирования- квазисинхронных орбит КА вокруг Фобоса; позволяющий, вычислить параметры .орбиты КА, которая в заданное время, проходит над заданной- долготой^ поверхности Фобоса, и обладает свойством минимального дрейфа, т.е. проходит заданную долготу примерно на одинаковом- расстоянии- от поверхности Фобоса. Метод позволил предварительно вычислить набор таблиц, каждая- из которых соответствует значению максимального удаления от поверхности Фобоса. На основе этих таблиц простой вычислительный алгоритм определяет параметры^ орбиты КА, которая^обладает свойством минимального дрейфа и в заданное время проходит над заданной, долготой-поверхности Фобоса.

2'. Разработаны, методы бортовой навигации, обеспечивающие-сближение с малым-небесным телом и посадку на его поверхность. Методы предназначены для их реализации на бортовом компьютере, предполагают использование цифровой модели поверхности небесного тела, используют измерения дальности и скорости в направлении лучей навигационных приборов и позволяют работать на фоне работающих двигателей. Алгоритмы и бортовые программы навигационного обеспечения участка посадки проекта «Фобос-Грунт» разработаны на основе этих методов.

3. Разработан метод определения параметров движения КА в условиях воздействия шума. Метод основан на представлении модели движения- КА в виде композиции опорного движения-и движения относительно опорного. Движение относительно опорного представляет собой случайный процесс с нулевым средним и описывается линейным стохастическим дифференциальным уравнением. Оценка вычисляется из условия минимизации функционала, содержащего квадраты взвешенных невязок между измеренными значениями и их расчётными моделями, а также взвешенные нормы векторов суммарных возмущений между измерениями. Метод применён в алгоритмах определения- параметров движения КА на фоне работы двигательной установки. На основе этих методов разработаны, алгоритмы» определения» параметров движения КА с электроракетными двигательными установками (ЭРДУ).

4. Разработан' метод оценки; точности приведения КА в заданную область, в пространстве: целевых параметров с учётом; ошибок исполнения и навигационных ошибок. В качестве целевых параметров могут быть выбраны параметры картинной плоскости при решении; задачи приведения КА к планете или Луне; высота над поверхностью, вертикальная и горизонтальная составляющая^ вектора скорости при решении задачи обеспечения- условий начала, сеанса посадки на поверхность Луны; параметры, определяющие: геометрию гало-орбиты, при решении задачи перелёта на гало-орбиту. Метод использован при оценке точности приведения КА в проектах «Фобос-Грунт», «Луна-Ресурс», «Луна-Глоб»^ «Венера-Д».

Научная новизна разработанных методов.заюпочаетсяш следующем: 1. Новизна разработанного метода проектирования квазисинхронных орбит КА вокруг Фобоса состоит в том, что усреднённая система учитывает дрейф эллипса; по которому совершается; облёт К А вокруг Фобоса, в направлении,- которое ортогонально, орбитальному движению Фобоса. Наличие этого; уравнения* позволило найти начальные приближения для-орбит с минимальным дрейфом, а затем построить искомые семейства орбит., 2. Новизна методов бортовой навигации на участке посадки состоит в том, что найдены, такие модели движения КА и построены такие модели измерений; которые, с одной; стороны, обладают простотой, позволяющей реализовать их на бортовом компьютере, а с другой стороны, точностью, позволяющей решить поставленную задачу.

3. Научная новизна метода оценки параметров движения КА в условиях воздействия шума состоит в предложенной модели движения; представляющей собой композицию детерминированного опорного движения и случайного процесса, который описывается линейным стохастическим дифференциальным уравнением, а также в методе, который обеспечивает одновременное получение оценки вектора состояния КА и векторов суммарных возмущений между измерениями.

Заключение

1. Г.: Абрамович; С.К;,, Агеева Т.Д., Аким Э:Л. и др. Баллистика и навигация автоматических межпланетных станций" «Венера-9»- и «Венера-10». // Космические исследования;-1976. Т. 14. Ш 51— С. 667-676;

2. Абрамович С.К., Агеева Т.Д., Аким Э Л. и др. Баллистика и навигация автоматических межпланетных станций «Венера-11» и «Венера-12». // Космические исследования. — 1979. Т. 17. № 5. - С. 670-680.

3. Агаджанов П.А., Дулевич В.Е., Коростелев A.A. Космические траекторные измерения: Радиотехнические методы; измерений и математическая обработка данных. М.: Советское радио, 1969. - 504 с.

4. Агеева Т.Д., Аким, ЭЛ.,. Иванов Н.М. и; др. Баллистика и навигация автоматических; межпланетных станций «Венера-13» и «Вёнера-14» // Космические исследования. — 1983. Т. 21. № 2. - С. 160-169.

5. Аким Э.Л. Определение поля тяготения Луны по движению искусственного спутника Луны «Луна-10» // ДАН СССР; 1966. -Т. 170, № 4. -- С.799-802.

6. Аким Э.Л;, Бажинов И.К., Павлов В.П., Почукаев В.Н. Поле тяготения Луны и движение её искусственных спутников. М.: Машиностроение,. 1984.-288 с.

7. Аким Э.Л., Власова З.Г1. Исследования гравитационного поля Луны по ' данным измерений траекторий советских' ИСЛ. // Космические исследования. 1983. - Т. 21, Вып. 4.

8. Аким Э.Л;, Власова 3.П. Модель гравитационного поля Луны по наблюдениям за движением: ее искусственных спутников Луны «Луна-10, -12, -14, -19 и -22» //ДАН СССР. 1977. - Т. 235, № 1.-С.38-41.

9. Аким Э.Л;, Горохова A.A., Киселева И.П., Степаньянц В.А., Тучин А.Г. Небесно-механическая интерпретация и первичная обработка измерений КА «Гранат» и «Интербол». // Препринт ИПМ им. М:В. Келдыша РАН. -1996.-№83.-21 е.

10. Аким Э.Л., Горохова A.A., Степаньянц В.А., Тучин А.Г. Контроль траектории выведения по данным траекторных измерений. // Препринт ИПМ им. М:В. Келдыша РАН. 1999. - № 10. - 26 с.

11. АкимЭ.Л., Степаньянц В. А. Численная теория движения Венеры по результатам радиолокационных наблюдений // ДАН СССР. 1979. -Т. 233. -С. 314-318.

12. Аким Э.Л., Энеев Т.М. Определение параметров движения космического летательного аппарата по данным траекторных измерений // Космические исследования. — 1963—Т. 1, Вып. 1. — С.5-50

13. АокиМ. Оптимизация стохастических систем. М.: Наука, 1967. - 424 с.

14. Бакитько Р.В., Васильев М.Б., Винницкий A.C. и др. Радиосистемы межпланетных космических аппаратов. М.: Радио и связь, 1993. - 328 с.

15. Берке У. Пространство-время, геометрия, космология. — М.: Мир, 1985. — 414 с.

16. Брумберг В.А., Финкелыптейн А.М. Проблема сопоставления результатов наблюдений и теоретических данных при проверке релятивистских эффектов в Солнечной системе. // Журнал экспериментальной и теоретической физики. — 1979. Т. 76, Вып. 5.

17. Брумберг В.А. Релятивистская небесная механика. — М.: Наука, 1972. — 382 с.

18. Брумберг В.А. Релятивистские эффекты при радиолокационных, оптических и радиоинтерферометрических измерениях. // Астрономический журнал. 1981. - Т. 58; Вып. 1. - С. 181-183

19. Градштейн И.С., Рыжик И.М. Таблицы интегралов сумм, рядов и произведений. М,: Наука, 1971. - 1108 с.

20. Дубошин Г.Н. Небесная механика. Аналитические и качественные методы. -М: Наука, 1968. 456 с.

21. Дубошин Г.Н. Небесная механика. Основные задачи и методы. М:

22. Наука, 1975.-800 с. 221 Егоров В :А.,ГусевЛ.И. Динамика* перелётов* между Землей и Луной.

23. М.,«Наука», 1980. 543 с. 231 ЗасоваШ®;, М0роз!В:И;, Линкин BtM: и др. Строение атмосферы Венеры от поверхности до 100 т. // Космические исследования. - 2006.- TV 44: № 4: С 381-400i

24. Засова Л:В;, Мороз'В:№, Формизано В. и др. Исследование Венеры с помощью ИК. Фурье-спектрометра на Венере-15 и Планетного Фурье-спектрометра на борту Венеры-Экспресс. // Космические исследования:.- 2006: Т. 441 № 4: - С. 365-380.

25. Зельманов А.Л., Агаков В.Г. Элементы общей теории относительности.- Mi: Наука;, 19892-240:с.

26. Иванов H.M., Поляков B.C. Наведение автоматических межпланетных станций:- М;: «Машиностроение», 1987. 312 с.

27. Квашнин А.Г., Тучин А.Г. Баллистическое обеспечение радиоинтер-ферометрических измерений. // Препринт ИПМ им. М.В. Келдыша РАН.- 1992.-№98.-26 с.

28. Кислик М.Д. Релятивистские возмущения при астрометрических наблюдениях планет. // Письма в АЖ. 1980:.- Т. 6, № 12.

29. Коган А.Ю. Далёкие спутниковые орбиты в ограниченной круговой задаче трёх тел. // Космические исследования. — 1988. Т. 26, № 6. - С. 813-818.

30. Колосов М.А., Арманд H.A., Яковлев О.И. Распространение радиоволн при космической связи. — М.: Связь, 1969; — 155 с.

31. Константинов М.С., Попов Г.А., Федотов Г.Г. Оценка использования солнечной электрореактивной установки для выведения спутника Юпитера. // Космические исследования. 2002. - Т. 40, № 2. - С. 201208.

32. Котельников В.А., Барсуков B.JL, Аким Э.Л. и др. Атлас поверхности Венеры. — М.: Главное управление геодезии и картографии при Совете Министров СССР, 1989. 328 с.

33. Красовский H.H. Теория управления движением. М.: Наука, 1968. -476 с.

34. Ли Р. Оптимальные оценки; определение характеристик и управление. — М.: Наука, 1966. 176 с.

35. ЛидовМ.Л., Вашковьяк М.А. О квазиспутниковых орбитах в ограниченной эллиптической задаче трёх тел. // Письма в АЖ. 1994'. — Т. 20, №10.-С. 781-795.

36. ЛидовМ.Л:, Вашковьяк М.А. О квазиспутниковых орбитах для эксперимента по > уточнению гравитационной постоянной. // Письма в АЖ. 1994. - Т. 20, № 3. - С. 229-240.

37. ЛидовМ.Л., Вашковьяк М.А. Теория возмущений и анализ эволюции квазиспутниковых орбит в ограниченной задаче трёх тел. // Космические исследования. 1993. - Т. 31, № 2. - С. 75-99.

38. Лидов М.Л., Рабинович В.Ю. Исследование семейств пространственных периодических орбит задачи трёх тел. // Космические исследования. -1979.-Т. 17. №3.

39. ЛидовМ.Л. Метод построения семейств пространственных периодических орбит в задаче Хилла. // Космические исследования. — 1982. Т. 20, № 6. - С. 787-807.

40. ЛидовМ.Л. Об' одном семействе пространственных периодических орбит около Луны и планет. // ДАН СССР. 1977. - Т. 233, № 6. - С. 1068-1071.

41. Льюнг Л. Идентификация систем. Теория для пользователя. М.: Наука, 1991.-432 с.

42. Маров М.Я. Планеты Солнечной .системы. М: Наука, 1986. - 320 с.

43. Молотов Е.П. Наземные радиотехнические системы управления космическими аппаратами. М.: ФИЗМАТ ЛИТ, 2004. - 256 с.

44. МюррейК., ДермоттС. Динамика солнечной системы М.: ФИЗМАТЛИТ, 2009. - 588 с.

45. Острём К. Введение в стохастическую теорию управления. М.: Мир, 1973.-322 с.

46. Охоцимский Д.Е., Сихарулидзе Ю.Г. Основы механики космического' полёта. М., Наука, 1990. - 448 с.

47. Платонов А.К. Исследование свойств корректирующих маневров в межпланетных полётах // Космические исследования. — 1966. Т. 4, Вып. 5.-С. 670-693.

48. Платонов А.К. Оптимальные свойства корректирующих маневров при использовании двигателя с ограниченной тягой // Космические исследования. 1967. — Т. 5, Вып. 2.

49. Растригин Л.А. Системы экстремального управления. -М.: Наука, 1974. -632 с.

50. Решетнев М.Ф., Лебедев A.A., Бартенев В.А. и др. Управление и навигация искусственных спутников Земли на околокруговых орбитах. — М.: Машиностроение, 1988, — 336 стр.

51. Себехей В. Теория орбит. Ограниченная задача трёх тел. М: Наука, 1982.-656 с.

52. Сейделмен П. Спутники Марса. Mi: Мир, 1981". - 100 с;

53. Скворцов A.B. Триангуляция Делоне и её применения. Томск.: Изд-во Томского ун-та, 2002. — 128 с.

54. Титов Д.В., Сведхем X., МкКой Д. и др. Венера-Экспресс: научные задачи, аппаратура и сценарий миссии. // Космические исследования. -2006. Т. 44, № 4. - С. 349-364.

55. ИПМ им. М.В; Келдыша РАН. 2004. - № 2'. - 32 с. '1.59; Тучин А.Г. Проектирование квазисинхронных орбит КА вокруг Фобоса1.для решения задачи посадки на его поверхность. // Препринт ИПМ им.j М.В. Келдыша РАН. 2008 - №15. - 30 с.

56. Будущее прикладной математики: Лекции для молодых1.исследователей: Поиски и:открытия / Под ред. Г.Г. Малинецкого. М.:

57. Книжный дом «Либроком», 2009. С.363-390.1., 62. Тучин АлГ., Аким Э.Л., Астахов А.П., Бакитысо Р.В:, Польщиков В.П.,

58. Степаньянц В.А., Тучин Д.А., Ярошевский В;С. Автономная навигации• онная система околоземного космического аппарата. //Известия: РАН;.

59. Теория и системы управления. 2009. - №2 - С. 139-156.

60. Тучин А.Г., Аким Э.Л., Архангельский P.H., Зайко Ю.К., Лавренов С.

61. М., Порошин А.Л., Рузский Е.Г., Степаньянц В:А., Тучин Д.А., ,

62. Федотов-В:П., Ярошевский В;G. Баллистика, навигация и управление движением КА на этапе его посадки на поверхность Фобоса: // Космические исследования. 2009. - Т. 47, № 4. - С.ЗЗ2-342.

63. Тз^чин А.Г., Аким Э.Л., Горохова A.A., Киселева И.П;, Степаньянц В;А. Локальная обработка измерений радиосистемы межпланетных космических аппаратов. // Препринт ИПМ им. М.В. Келдыша РАН, 2002.-№ 11.-20 с.

64. Тучин А.Г., Аким Э.Л., Заславский F.C., Морской И.М., Степаньянц B.A. Баллистика, навигация и управление полётом космического аппарата в проекте «Фобос-Грунт» // Известия РАН! Теория и системы управления. 2002. - № 5, - С. 153-161.

65. Тучин А.Г., Аким Э.Л., Сазонов BIB:, Степаньянц В.А. Управление и навигация в перелёте КА с малой тягой от Земли к Юпитеру. // Препринт ИПМ им. М.В. Келдыша РАН: 2010. - № 36.

66. Тучин A.F., Квашнин А.Г. Построение- вычислительной схемы алгоритма: оценивания параметров; динамической системы. // Автоматизация проектирования, информатика. 1993; — Вып. 4. — С.88-94.

67. Тучин А.Г., Квашнин А.Г. Синтез вычислительной схемы алгоритма оценивания; параметров динамической системы. // Препринт ИПМ им. М.В. Келдыша РАН. 1991. - № 126 - 19 с.

68. Тучин А.Г., Шатаев В.В. Реализация реляционной СУБД с языком данных типа SQL. // Препринт ИПМ им. М.В. Келдыша РАН. 1990. -№ 54.-28 с.

69. Чёрный Ф:Б. Распространение радиоволн. М.: Советское радио, 1972; 464 с. ' ■ ■ '

70. Шишов В.А. Модель движения/ Фобоса, ш методика уточнения; параметров; в: проекте; «Фобос-Грунт». // Препринт И11М им. М.В. Келдыша РАИ. 2008. - № 10.-27 с.

71. Шишов В.А. Определение параметров движения КА и Фобоса в проекте «Фобос-Грунт». // Астрономический вестник. 2008. - Т. 42, № 4. -С.341-350.

72. AkimiE.L., Stepaniants V.A., Tuchin;A.G: Tracking of the launch-vehicle* during the insertion to Earth orbit // RBCM J. of the Braz. Soc. Mecanical Sciences.- Vol; XXL - Special lssue-1999;-PL387-3991

73. Akim E.L., Zaslavsky G.S., Zharov V.G., Chernov A.V. Interplanetary Flight Control With Electric Engine In View Of Thrust Errors // The Proceedings of the 18th International Symposium on Space Flight Dynamics, 11-15 October 2004, Munich, Germany.

74. Berman A.L. A New Approach to the Evaluation and Prediction of Wet Tropospheric Zenith Wet Refraction. JPL DEEP SPACE Progress Report 4225.

75. Berman A.L. The Prediction of Zenith Range Refraction From Surface Measurements of Meteorological Parameters. JPL Technical1 Report 32-1602/ 1976.-40 p.

76. AmoozegarF. Applications and Operation Concepts of Large Transmit Phased Array of Parabolic Reflectors. I IIEEEAC paper #1570, 2006.

77. Benest D. Libration effects for retrograde satellites in the restricted three-body problem. // CelestiabMechanics. 1976. -V. 13, № 2. - P. 203-215.

78. Bhaskaran S. The application* of noncoherent Doppler data types for deep space navigation. // TDAProgress Report 42-121. JPL, 1995.

79. Carpenter J.D., Folta D.C., Moreau М.С., Quinn D.A. Libration Point Navigation Concepts Supporting the Vision for Space Exploration. // AIAA 2004-4747, AIAA/AAS < Astrodynamics Specialist Conference and Exhibit, Providence, Rhode Island, August 2004.

80. CrouseB., ZanettiR., D'SouzaC., Spanos P.D. Autonomous optical lunar navigation // AAS 09-122, Savannah, GA, February 2009.

81. D'Souza C., Clark F. Linear covariance analysis techniques applied to Orion cislunar operations // AAS-08-100, Galveston, TX, January 2008.

82. Estefan J.A., Folkner W.M, Sensitivity of Planetary Cruise Navigation to Earth Orientation Calibration Errors. TDA Progress Report 42-123.

83. Ferri P., Accomazzo A., Montagnon E., Morales J. Rosetta in the year of the swing-bys // IAC-07-A3.5.02. Доклад на 58-м международном конгрессе по астронавтике. Индия. 24-28 сентября 2007.

84. Getchius J., Crain Т., D'Souza С. Optical navigation for the Orion vehicle. // AAS-08-105, Galveston, TX, January 2008.

85. HauptG., KasdinN., KeiserG., Parkinson В. An optimal recursive iterative algorithm for discrete nonlinear least-square estimation. // AIAA-95-3218. AIAA Guidance, navigation, and control conference. August 7-10, 1995. Baltimore, P. 404-417

86. Henon M. Numerical Exploration of the Restricted Problem I IVI. Hill's Case: Non-Periodic Orbits. Astronomy and Astrophysics. 1970. - V. 24, № 9. - P. 24-36.

87. HofmannP., von Heise-Rotenburg R., Schulte W., Widani C. Sample preparation and handling system & payload* aspects for the european Exomars mission // IAC-07-A3.3.04. Доклад на 58-м международном конгрессе по астронавтике. Индия. 24-28 сентября 2007.

88. Hoofs R.M.T., TitovD., SvedhemH., KoschnyD., Witass О., Tanco I. Venus-Express science observations experience at Venus // IAC-07-A3.2.03. Доклад на 58-м международном конгрессе по астронавтике. Индия. 24-28 сентября 2007.

89. Jones В.A. Surface feature navigation in low lunar orbit // AAS 08-104, Galveston, TX, January 2008.

90. McNutt R.L., Solomon S.C., Grant D.G., Finnegan E.J., Bedini P.D. The MESSENGER Mission on to Mercury: Status after the Venus Flybys // IAC-07-A3.2.02. Доклад на 58-м международном конгрессе по астронавтике. Индия. 24-28 сентября 2007.

91. Mitchell R.T. The Cassini Mission Exploring Saturn // IAC-07-A3.2.01. Доклад на 58-м международном конгрессе по астронавтике. Индия. 2428 сентября 2007.

92. Moyer. T.D. Formulation for Observed and Computed Values of Deep Space Network Data Types for Navigation. Jet Propulsion Laboratory California Institute of Technology. - 2000.

93. Opperman B.D.L. Reconstructing ionospheric TEC over South Africa using signals from a Regional GPS network. A thesis submitted in fulfillment of the requirements for the degree Doctor of philosophy of Rhodes University. 28 Nov. 2007.

94. Osenar M., Clark F., D'SouzaC. Performance of an automated feature tracking lunar navigation system // AAS 08-101, Galveston, TX, January 2008.

95. Seidelmann Р.К., Abalakin V.K. et all Report of the IAU/IAG Working

96. Group on cartographic coordinates and rotational elements of the planets and\satellites: 2000 // Celestial Mechanics and Dynamical Astronomy, 2002, 82(1).

97. Shor V.A. The motion of the Martian satellites // Celestial Mechanics, Vol. 12, 1975.

98. Takayuki Sugano and Kosuke Heki. High resolution lunar gravity anomaly map from the lunar prospector line-of-sight acceleration data. // Earth Planets Space, 56, 2004.

99. Chao C.C. The Tropospheric Calibration Model for Mariner Mars 1971. // JPL Technical Report 32-1587. -P.61-76.

100. Thurman S.W. Deep space navigation with differenced data types. Part 1: Differenced range information content // The telecommunications and data, acquisition progress report 42-103, JPL, 15, 1990. P. 47-60.

101. Turyshev S.G., Nieto M.M., Anderson J.D. Lessons Learned from' the Pioneers 10/11 for a Mission to Test the Pioneer Anomaly. // Adv. Space Res. 39(2), 2007.-P. 291-296.

102. Turyshev S:G., Nieto M;M., Anderson J;Di Study of the Pioneer anomaly: A problem set. // American J. Phvs. 73(11), 2005. P. 1033-1044.118; Mendes V.B., LangLey R.B. Revised Refractive Index Formulae and Their

103. Анализ и дальнейшее развитие принципов обеспечения навигации и управления: КА дальнего космоса. : отчёт о НИР / ИГ1М им. М.В. Келдыша РАН. -М., 2009. -78 с.- Инв. № 5-016-09.

104. Баллистико-навигационное обеспечение выведения КА «Фобос-Грунт» на траекторию-отлёта от Зёмли; : отчёт о;НИР;У ИИМ им. М.В. Келдыша РАН М;, 2006: - 116 е. - Инв; №5^-04-06:

105. Баллистико-навигационное обеспечение заключительного этапа полета орбитальной станции «МИР» (БЦ ИПМ). : отчёт о НИР / ИПМ им; М.В ; Келдыша РАН: М., 2001.- Инв;.№ 5-01-02.

106. Баллистико-навигационное обеспечение проекта «Спектр-РГ». : отчёт о НИР / ИПМ им; М;В; Келдыша PAHL -M;, 2009. 137 с. - Инв; Ж 5-0509. ,

107. Баллистико-навигационное обеспечение старта с Фобоса и выхода на базовую орбиту. : отчёт о НИР / ИПМ им. М.В. Келдыша РАН. М., 2005, - 116 с. -Инв. № 5-10-05.

108. Навигационное обеспечение проекта «Луна-Глоб1». : отчёт о НИР / ИПМ им. М.В. Келдыша РАН. M., 2008! -120 с, - Инв; № 5-06-08-.

109. Оптимизация баллистического г построения экспедиции с луноходом и окололунным орбитальным аппаратом (Проект «Луна — Ресурс»). : отчёт о НИР / ИПМ им; М.В. Келдыша РАН. М., 2009. - 105 с. - Инв. № 5018-09.

110. Проработка предложению по схеме: полёта к Юпитеру и проведению работ по баллистической поддержке миссии: к Юпитеру и Европе на участке:, перелёта Земля Юпитер. : отчёт о НИР / ИПМ: им. М.В. Келдыша РАН. -М., 2009. - 58 е.- Ипв. № 5-012-09. ,

111. Разработка алгоритмов математической! модели сближения и посадки: КА на Фобос: с; квазисинхронной орбиты: : отчёт о НИР / ИПМ им: М.В: Келдыша РАИ. - М:, 2005. - 99 с. - Инв. № 5-07-05.

112. Разработка программных модулей основных алгоритмических блоков для модели сближения и посадки КА на поверхность Фобоса с использованием цифровых моделей Фобоса. : отчёт о НИР / ИПМ им. М.В. Келдыша РАН. М., 2005. - 80 с. - Инв. № 5-11-05.

113. Разработка технических предложений по схеме перелёта Земля-Венера, маневрирования в сфере действия гравитационных сил Венеры и баллистико-навигационному обеспечению полета. : отчёт о НИР / ИПМ им. М.В. Келдыша РАН. М., 2009. - 68 с. - Инв. № 5-013-09.