автореферат диссертации по авиационной и ракетно-космической технике, 05.07.09, диссертация на тему:Синтез алгоритмов управления космическими аппаратами с учетом требований безопасности проведения динамических операций

кандидата технических наук
Старков, Александр Владимирович
город
Москва
год
2012
специальность ВАК РФ
05.07.09
цена
450 рублей
Диссертация по авиационной и ракетно-космической технике на тему «Синтез алгоритмов управления космическими аппаратами с учетом требований безопасности проведения динамических операций»

Автореферат диссертации по теме "Синтез алгоритмов управления космическими аппаратами с учетом требований безопасности проведения динамических операций"

На правах рукописи УДК 629.783

Старков Александр Владимирович

СИНТЕЗ АЛГОРИТМОВ УПРАВЛЕНИЯ КОСМИЧЕСКИМИ АППАРАТАМИ С УЧЕТОМ ТРЕБОВАНИЙ БЕЗОПАСНОСТИ ПРОВЕДЕНИЯ ДИНАМИЧЕСКИХ ОПЕРАЦИЙ

Специальность 05.07.09 «Динамика, баллистика, управление движением летательных аппаратов»

АВТОРЕФЕРАТ

диссертации на соискание учёной степени кандидата технических наук

1 и ЯНВ 2013

МОСКВА

005048074

Работа выполнена на кафедре «Системный анализ и управление» Федерального государственного бюджетного образовательного учреждения высшего профессионального образования «Московский авиационный институт (национальный исследовательский университет)» (МАИ).

Научный руководитель: Малышев Вениамин Васильевич, доктор технических

наук, профессор, заведующий кафедрой «Системный анализ и управление» Московского авиационного института (национального исследовательского университета) (МАИ)

Официальные оппоненты: Падалко Сергей Николаевич, доктор технических наук,

профессор, заместитель заведующего кафедрой «Прикладная информатика» Московского авиационного института (национального исследовательского университета) (МАИ)

Суханов Константин Георгиевич, кандидат технических наук, доцент, ведущий специалист центра 128 Федерального государственного унитарного предприятия научно-производственное объединение имени С.А. Лавочкина, (ФГУП «НПО им. С.А. Лавочкина», 141400, г. Химки, Московская область, ул. Ленинградская, д. 24).

Ведущая организация: Федеральное государственное унитарное предприятие

«Центральный научно-исследовательский институт машиностроения» (ФГУП ЦНИИмаш, 141070, Московская область, г. Королёв, ул. Пионерская, дом 4)

Защита состоится «13» декабря 2012 г. в 13.00 часов на заседании диссертационного совета Д 212.125.12 в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Московский авиационный институт (национальный исследовательский университет)» (МАИ) по адресу: 125993, г. Москва, А-80, ГСП-3, Волоколамское шоссе, д.4.

С диссертацией можно ознакомиться в библиотеке Федерального государственного бюджетного образовательного учреждения высшего профессионального образования «Московский авиационный институт (национальный исследовательский университет)» (МАИ) по адресу: 125993, г. Москва, А-80, ГСП-3, Волоколамское шоссе, д.4.

Автореферат разослан «12» ноября 2012г.

Отзывы, заверенные печатью, просим высылать по адресу: 125993, г. Москва, ГСП-3, А-80, Волоколамское шоссе, д.4, Ученый совет МАИ.

/ ^

Ученый секретарь диссертационного совета

Д 212.125.12, к.т.н., доцент В.В. Дарнопых

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. В результате человеческой деятельности количество управляемых и неуправляемых космических объектов в околоземном пространстве с каждым годом будет расти, а вероятность столкновения космических объектов друг с другом - увеличиваться. Столкновения и разрушения космических аппаратов (КА) приводят к образованию новых крупных фрагментов и осколков. По приблизительным оценкам, в околоземном пространстве находится более девяти тысяч объектов с размерами более 10 сантиметров, десятки тысяч с размерами 1-10 сантиметров и сотни тысяч с размерами менее 1 сантиметра. В последнее время «космический мусор» стал представлять все большую опасность для вновь запускаемых КА, но приемлемые способы борьбы с ним пока не разработаны. Поэтому потребность в обеспечении безопасного проведения динамических операций в околоземном пространстве очевидна. Уменьшение угрозы столкновения КА с другими космическими объектами представляет собой сложную организационно-техническую задачу, решению одной из подзадач которой посвящена данная работа.

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

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

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

Актуальность и новизна работы следует из анализа развития специализированного программного обеспечения, как на мировом рынке, так и во внутрикорпоративном сегменте (продукты Satellite Tool Kit, Bernese Software, «Orbit Determination Toolkit», «Free Flyer» и др.), который выявил тенденцию их

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

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

2) единого, стандартизованного языка высокоуровневого описания объектов моделирования, входной и выходной информации, исполнительной макропрограммы с поддержкой механизмов межпрограммного взаимодействия;

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

Анализ публикаций на данную тему и результаты работ, выполненных авторами: В.А. Бартеневым, М.Н. Красилыциковым, A.A. Лебедевым, JI.H. Лысенко, В.В. Малышевым, В.Н. Почукаевым, М.Ф. Решетневым, Г.М. Чернявским, Р. Шенноном, позволил определить основные особенности построения алгоритмов управления космическим аппаратом при проведении динамических операций, а также пути исследований, направленных на определение общих принципов и методов решения поставленной задачи. В работах этих авторов излагались теоретические основы и методы проектирования КА, исследовались проблемы навигационно-баллистического обеспечения полета, модели движения и принципы управления КА, работа целевой аппаратуры и служебных систем. Вместе с тем остается еще ряд задач, связанных с созданием унифицированных средств автоматизации отработки динамических операций КА.

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

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

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

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

построения программно-моделирующего комплекса (ПМК) для отработки динамических операций.

Для достижения поставленной цели решаются научно-технические задачи:

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

2. Разрабатываются математические модели бортовых систем КА и внешней среды.

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

4. Определяются принципы построения и технический облик программно-моделирующего комплекса.

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

Методы исследования. В качестве методологической основы используется системный подход. На нем основывается принцип создания единой базовой подсистемы - типового унифицированного программного средства, методический и алгоритмический аппарат которого позволяет определять движение широкого множества КА ближнего и среднего космоса, гибко реагируя при этом на изменение параметров орбит КА, допустимых значений показателей качества и различных моделей измерений. Основными методами исследования, используемыми в работе, являются методы динамики полета, теории управления, статистические методы обработки данных, методы оптимального управления и обработки информации. При программной реализации математического обеспечения используются методы объектно-ориентированного программирования и мультизадачность операционных систем Windows, Linux, системы управления реляционными базами данных, а также средства обеспечения доступа в сетях Интернет и Интранет.

Научная новизна результатов. В работе получены следующие результаты, обладающие новизной и научной значимостью:

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

2. Определены математические модели проведения орбитальных коррекций КА при выполнении динамических операций.

3. Сформирована адаптированная к требованиям программно-моделирующего комплекса (ПМК) математическая модель бортовых систем, внешних воздействий и контура управления.

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

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

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

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

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

Результаты диссертационной работы используются в учебном процессе кафедры «Системный анализ и управление» МАИ и рабочей деятельности ФГУП ЦНИИмаш, ФГУП «НПО им. С.А. Лавочкина». Имеются соответствующие акты внедрения от МАИ, ФГУП ЦНИИмаш, ФГУП «НПО им. С.А. Лавочкина».

Достоверность результатов подтверждается использованием

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

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

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

3. Методика разработки базы данных программно-моделирующего комплекса.

4. Результаты использования программно-моделирующего комплекса для отработки безопасного проведения динамических операций на высокоэллиптических и круговых геостационарных орбитах (ГСО). Апробация работы и публикации. Результаты работы докладывались и

получили одобрение на научно-технических конференциях: 16-ой Международной научной конференции «Системный анализ, управление и навигация» (Евпатория, 2011 г.), 17-ой Международной научной конференции «Системный анализ, управление и навигация» (Евпатория, 2012 г.)

Основные результаты диссертационной работы опубликованы в трех статьях [1-3] в журналах, входящих в рекомендованный ВАКом Минобрнауки России перечень изданий, и в четырех работах [4-7] в сборниках тезисов докладов на научно-технических конференциях.

Структура и объём работы. Диссертация состоит из введения, четырёх глав, заключения, списка литературы из 51 наименования. Текст диссертации изложен на 218 машинописных страницах, включает 80 рисунков и 12 таблиц.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

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

• ввод спутника в заданное место (слот) орбитальной структуры после выведения носителем на промежуточную орбиту;

• коррекция относительного положения спутников, когда искажение орбитальной структуры приводит к снижению целевой эффективности системы;

• поддержание относительного положения спутников в заданных пределах при проведении экспериментов;

• перевод спутника из одного слота структуры в другой, например - при замене вышедшего из строя спутника резервным;

• удержание спутника в окрестности номинального положения слота в течение длительных промежутков времени (например, удержание спутника на ГСО, удержание долготы восходящего узла высокоапогейной синхронной орбиты);

• сближение с другим КА для инспекции;

• выполнение маневров уклонения на фоне решения КА его основной задачи;

• регулирование относительного положения спутников при проведении экспериментов и т.п.

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

1. Задаче приведения вектора начального состояния КА группировки в заданную область пространства.

2. Задаче выполнения маневров в заданной области пространства в течение срока активного существования.

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

Вторая группа требований это требования к программному обеспечению, которое должно:

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

• обеспечить возможность назначения пользователем критерия оптимальности или взвешенного обобщенного критерия с коэффициентами предпочтения;

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

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

• обеспечить «прозрачность» (понятность) генерируемых решений для конечного потребителя.

Требования, предъявляемые к средствам отработки динамических операций КА для исключения возникновения опасных ситуаций, позволили сформулировать две взаимосвязанные между собой задачи синтеза алгоритмов управления КА с учетом безопасности проведения динамических операций и формирования методики построения программно-моделирующего комплекса (ПМК) для их отработки.

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

Далее рассматривается постановка задачи управления при проведении следующих динамических операций:

• обеспечение начальной орбиты с заданными параметрами (задача выведения);

• перевод на промежуточную и возвращение на опорную орбиту при выполнении маневра уклонения;

• удержание КА на расчетной орбите.

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

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

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

Ук+г = Ук + [0* + (1 + Ик)ик№к + Ме* + (1)

1 = вк + (1 + и*)«* + 2Ь*Д1* + ; ь =

где к -индекс, соответствующий моменту начала проведения маневра; Ы— количество коррекций; ук - отклонение географической долготы восходящего узла (ГДВУ) от требуемого значения (в град.); -Эк - скорость изменения ук или скорость дрейфа ГДВУ (в град./зв.сут.); ик - корректирующее воздействие или приращение скорости дрейфа, обусловленное работой корректирующей двигательной установки (в град./зв.сут.); цк - мультипликативная ошибка реализации корректирующего воздействия ик,%Ук, - случайные аддитивные возмущения; Ьг - постоянные на интервале времени Д^ коэффициенты, вычисляемые по формулам Ьг = .»(А*^) > ъ = + ^ д

»0 уо

гдег, я - направления радиус-вектора КА и трансверсали соответственно, Ддгг, АдГ5 - проекции гравитационного возмущающего ускорения; &дтзг, ~

проекции осредненного на интервале Д Ск возмущающего ускорения от гравитационных полей Луны и Солнца.

Статистические характеристики случайных факторов Ик^ук>^к считаются известными:

МЫ = MU = МЫ = °< Ч&] = М1Ш =

J = = al

где М[.] - математическое ожидание

Система уравнений плоского движения в матричном виде:

хк+1 = Акхк + (1 + ßk)Bkuk + Dk + ffc, (2)

где к - число маневров; хк - n-мерный вектор состояния системы; ик - ш-мерный вектор управления; Ак - квадратная пхп матрица;Вк - прямоугольная пхт матрица; Dk - и-мерный неслучайный вектор; цк - мультипликативная ошибка управления; - случайный вектор ошибок прогноза вектора состояния хк+1,

= (уAV. (, - (fy.f.,)' А* = Q, "{*)• в» =

В качестве характеристики конечной точности принимается величина У = min{Ufc) M[x%+1ÁxN+1], (3)

где X - симметричная матрица с известными элементами АцЛ^ Лгг-

Таким образом, задача оптимизации состоит в определении такого алгоритма коррекций [ик(хк), к = 1 ,N}, который позволит перевести КА из произвольного начального положения х0 в некоторое конечное xN+1 и обеспечит минимум критерия (3). Отметим, что предложенная модель движения содержит два типа аддитивных возмущений - систематические неслучайные Dk и случайные Этим она отличается от традиционных моделей, содержащих только случайные аддитивные возмущения. Решение поставленной задачи возможно в стохастической и минимаксной постановках.

В стохастической постановке оптимальное управление - линейная функция: ик = ~Lkxk - dk, (4)

где коэффициенты обратной связи Lk и смещение dk вычисляются с помощью рекуррентных соотношений:

Lk ~ Yk Afc+i^k. dk = Yk(4+A + Gk+i), У* = (1 +

Afc = Ak^k+l^k ~ ЬкУк^к,

Gk = Ak^k+lDk + ÄkGk+1 ~ Ykdktfo

Ck = cft+i + sp[F¿Ák+1FkEk] + D¡Ak+1Dk - ykd\ + 2G¡+1Dk с граничными условиями на правом конце: Лдг+i = Я; Слг+1 = 0; cN+1 = 0 .

В минимаксной постановке, когда задача состоит в поиске оптимального управления {ufc(;efc)}, которое обеспечивает минимакс критерия

ик = arg minUfc max^ Щ(х%+1Лхы+1)], (5)

оптимальное управление также является линейной функцией:

ик = = ~{ßTkXk+iAkxk + BlÄk+1Dk + Gfc+1ßfc)/(ßfcAfc+1ßfc),(6)

которое с точностью до слагаемого сг| совпадает с решением задачи в стохастической постановке.

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

вычисляются по формуле (4) при N=1. При этом интервалы между коррекциями ДСвыбираются из условия

< ДЛтах,

где ДЛтах - заданная величина, характеризующая интервал удержания по средней долготе.

Для этого случая получены скалярные соотношения для вычисления коэффициентов обратной связи и смещения:

¿у = (Д£ + Х2)(1 + о2У\ М2 + 2ХгДС + ХьГ1,

= (1 + ст2Г\ (7)

й = [Д Ь2(ЪГ + МО + ХгД*(ЗМС+Ьг) + 2Ь5Д«Х1](1 + (Т2)-\М2 + 2Х2Д£ + хО-1, где

V - V — ¿21

Л2 з > Л1 — 1 >

а индекс к, равный 1, у всех величин опущен (для сокращения записи).

Считая М постоянным, выражение для управления имеет вид

и = -1уу -Цд-й (8)

С учетом (8) и предполагая постоянство интервалов между коррекциями Аь, эволюция долготы и скорости дрейфа может быть описана следующим матричным уравнением:

М[*к+1] = Ахк + С, (9)

где

'\-Lybt

:ЬrДt+ЬsДt2 -2ЬЛД{ -й

Из соотношения (9) находим координаты КА х* в установившемся режиме: х* = (Е - Л)~1С. (10)

Равенство (10) показывает, что вектор состояния КА в установившемся режиме зависит лишь от интервалов времени между коррекциями ДС и требуемой амплитуды колебаний долготы восходящего узла и скорости дрейфа, определяемой параметрами матрицы Лд,+1. Такой подход позволяет на фоне режима удержания проводить динамические операции по незначительному изменению в нужный момент времени долготы восходящего узла орбиты КА и дальнейшему возвращению ее к расчетному значению без проведения дополнительных сложных расчетов.

А = С =

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

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

Для формирования терминального алгоритма предлагается заменить дифференциальные уравнения движения соотношениями, связывающими текущее г* и конечное Т состояния, при условии, что весь импульс скорости прикладывается в текущий момент времени:

S, (Г) = S'" (St,Vk, Atk ) + (1 + ni)AVk{Pl)k Atk, V, (T) = V (lk ,Vk,Ah) + (l + Л )AVk (P,)k, где l = x,y,z, Atk=T-tk;Sh V\- проекции векторов положения и скорости на оси выбранной системы координат,^ -проекции единичного вектора направления действия управляющей силы на оси выбранной системы координат (ориентация двигательной установки); AVk(P,)k -вектор проекций полного приращение скорости КА, вызванного работой двигательной установки на интервале Ык, на оси выбранной системы координат; Ser, V" - векторные функции прогноза состояния КА на интервал времени Aik, вычисляемые по кеплеровской теории; (1+) -мультипликативная ошибка, вызванная работой двигательной установки, причем

мкЫ

Определение пространственной ориентации корректирующей двигательной установки (ДУ) в этом случае сводится к задаче минимизации линейной функции F = F0+2^->min, (11)

где к- номер коррекции,

RT = Ry, Rz\

я,

F0 = [S(T) -S'Y [S(T)-s'] + bf +b\, è,=(l + at2) AVkIV0i Ь2=(\ + ст2к)А1кАУк/г0, S*, V - векторы требуемых положения и скорости КА, А К*-приращение характеристической скорости за время Д tk. Решение задачи определяет пространственную ориентацию вектора тяги:

и* г

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

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

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

Основным функциональным назначением ПМК является:

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

• организация дискретно-событийного взаимодействия моделей в модельном и в реальном времени с возможностью розыгрыша отказов и помех;

• управление процессом моделирования, как в диалоговом, так и в пакетном режиме работы, написание сценариев моделирования;

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

• регистрация и обработка результатов моделирования, взаимодействие с бортовыми регистраторами.

Предлагаемая методика включает в себя:

• определение состава и структуры ПМК в целом

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

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

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

В работе на основе технической декомпозиции КА предлагается состав и иерархическая структура программно-моделирующего комплекса (см. рис. 1).

На верхнем уровне находится ядро ПМК, которое взаимодействует с базой данных (БД), хранящей исходные данные для моделирования. Заложенный в нем принцип объектно-ориентированного программирования позволяет легко модифицировать

. „ исходное ядро в

Рис. 1. Структурная схема программного комплекса г ,

зависимости от фазы

полета или постановки решаемой задачи, не меняя при этом сущности

моделирования. В состав ядра ПМК входят управляющий модуль, графический

модуль и интегратор.

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

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

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

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

ЯДРО

Управляющий модуль

Интегратор

0

Модели бортовых систем КА

¡-ая фаза л о лета

Графический модуль

С

Модель внешней среды

/-ал Фаза полета

Модуль доступа к базе данных

Я

База данных 1-ая фаза полета

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

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

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

выполнены в соответствии со стандартами,

поддерживаемыми различными ОС. Более подробно архитектура управляющего модуля представлена на рис. 2. Модель операций

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

представляет собой

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

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

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

Рис. 2. Архитектура управляющего модуля

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

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

• интеллектуальное хранилище данных моделирования;

• обеспечение межпрограммного взаимодействия.

При таком подходе работа с объектом моделирования с точки зрения разработчика программного обеспечения состоит из нескольких этапов.

На первом этапе с использованием созданной автором системой команд (языком описания): определяется объектно-ориентированная иерархическую структуру объекта моделирования; формируются исходные и текущие данные для проведения моделирования; описываются методы обработки данных; организуется сценарий межпрограммного взаимодействия.

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

На третьем этапе проводится моделирование, состоящее из представления объекта моделирования с необходимыми программными модулями в памяти компьютера (компиляции проекта) и непосредственно самого моделирования с обеспечением возможностей «быстрого» обмена данными и событиями.

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

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

• Проведена отработка средств выведения КА с промежуточной круговой на рабочую орбиту.

• Проведена отработка средств выведения КА на круговую орбиту.

• Проведена отработка средств проведения динамических операций на высокоэллиптических орбитах.

• Проведена отработка средств проведения динамических операций на круговой геостационарной орбите.

Отработка средств выведения КА с промежуточной круговой на рабочую орбиту проиллюстрирована на примере выведения КА системы ГЛОНАСС. Исследуется участок траектории разгонного блока «Фрегат» с низкой круговой орбиты высотой Яг 180км до рабочей круговой орбиты с высотой 19135км. На этом участке проводятся два включения двигательной установки разгонного блока. Первый активный участок формирует эллиптическую орбиту с высотой апогея равной высоте рабочей орбиты КА ГЛОНАСС, а второй переводит КА на почти круговую рабочую орбиту. На первом участке циклограмма выведения задана жестко, т.е. обратной связи по измерениям параметров траектории нет. На втором участке может использоваться как жесткая циклограмма выведения, так и обобщенный терминальный алгоритм управления. На рис. 3-5 представлен характер изменения по времени орбитальных параметров как при отсутствии возмущение, так и при наличии предельных возмущений.

\п«„.

Р

Рис. 3. Изменение радиуса апоцентра по времени.

Рис. 4. Изменение радиуса перицентра по времени

Рис. 5. Изменение эксцентриситета по времени

Ниже представлены разбросы конечных параметров орбит при использовании на втором участке как жесткой циклограммы выведения (см. рис. 6), так и обобщенного терминального алгоритма (см. рис. 7)

Жесткая циклограмма выведения

гМВМ1гя,КХ1 * ^

•*- V '» * -

Алгоритм терминального управления

*т «Г«*

И00» 2М0*

»«04 Э2006

* г * .V'

Рис. 6. Разбросы радиусов перицентра и апоцентра при жесткой циклограмме выведения

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

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

На рис. 8-11 показаны результаты изменения периода обращения и перигея орбиты с учетом и без учета проведения маневров:

.....4.....;..... ......;......1.....[-■- г* '

.....}■■■■].....}■•- -

.....!......>■■■ X : ;

.....!......Г' "!..........'""' ..........1......

.....!......г ■ "!......Г.....1..... .....1......!Ч"'......

.....|.....Г' .....|.....1 ' ; [ ]ч

й

од

Рис. 8. Эволюция периода обращения

Рис. 9. Изменение периода обращения после проведения коррекции

""ГТ~

Рис. 10. Эволюция высоты перигея Рис. 11. Изменение высоты перигея после

проведения коррекции

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

__независимо, причем алгоритм

-0.1»|

' > V *р*а./|«.с)гг.

-9.1

т»я су.

У'Р'Я-

8.1";

II

коррекции наклонения - заранее выбранная программа

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

представляет собой

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

средней долготы. Начало очередного сеанса коррекций средней долготы, их продолжительность рассчитывается в соответствии с

алгоритмом, описанным в четвертой главе. Проведено моделирование процесса

коррекций долготы восходящего узла спутника на стационарной орбите: г=42164км, е^Ю'4, £2=30°, ¡=4°, Хтр~45°. На рис. 12 представлены фазовые портреты в координатах: долгота восходящего узла - скорость дрейфа для разных начальных условий и рассчитанной по формуле (10) точкой сходимости процесса управления (у*, в*).

Рис. 12. Зависимость скорости дрейфа о от долготы восходящего узла у.

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

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

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

1. Сформированы требования к алгоритмическому и программному обеспечению средств отработки динамических операций.

2. Предложены математические модели бортовых систем КА и внешней среды.

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

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

5. Создан программно-моделирующий комплекс для отработки алгоритмов и программ проведения динамических операций КА.

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

1. Вывод КА на расчетную высокоэллиптическую и круговую орбиты, как штатными алгоритмами выведения, так и с помощью терминального метода.

2. Выполнение маневров на высокоэллиптической орбите.

3. Выполнение маневров на геостационарной орбите.

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

ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

Публикации в изданиях из рекомендованного ВАК Минобрнауки России перечня:

1. Малышев В.В., Старков A.B., Федоров A.B. Методика разработки программно-моделирующего комплекса для отработки средств проведения динамических операций космических аппаратов. Электронный журнал «Труды МАИ», 2012, №57.

2. Малышев В.В., Старков A.B., Федоров A.B. Программно-моделирующий комплекс для отработки средств проведения динамических операций космических аппаратов. Вестник компьютерных и информационных технологий, 2012, №9, Москва: Издательский дом «СПЕКТР», с. 7-12.

3. Малышев В.В., Старков A.B., Федоров A.B. Методика актуализации банков данных объектов космического назначения. Вестник Московского Авиационного Института, 2012, №4. Москва: Изд-во МАИ, с. 56-64.

Другие публикации:

4. Малышев В.В., Старков A.B., Федоров A.B. Моделирующий комплекс для отработки средств безопасного проведения динамических операций КА. Тезисы докладов XVII международной конференции «Системный анализ, управление и навигация», 2012, Москва: Изд-во МАИ-ПРИНТ, с. 79.

5. Старков A.B., Шмигирилов С.Ю. Применение трёхмерной визуализации в задачах моделировании движения КА. Тезисы докладов XVII международной конференции «Системный анализ, управление и навигация», 2012, Москва: Изд-во МАИ-ПРИНТ, с. 61.

6. Малышев В.В., Старков A.B. Автоматизация процессов обработки данных в современных геоинформационных системах. Тезисы докладов XVI международной конференции «Системный анализ, управление и навигация», 2011, Москва: Изд-во МАИ-ПРИНТ, с. 96.

7. Малышев В.В., Соколов A.A., Старков A.B. Использование современных образовательных технологий для подготовки специалистов в области применения результатов космической деятельности. Тезисы докладов XVI международной конференции «Системный анализ, управление и навигация», 2011, Москва: Изд-во МАИ-ПРИНТ, с. 110.

Множительный центр МАИ (НИУ) Заказ от 09, Н. 2012 г. Тираж/¿/¿у* экз.

Оглавление автор диссертации — кандидата технических наук Старков, Александр Владимирович

Введение.

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

1.1. Причины возникновения опасных ситуаций с КА.

1.2. Основные направления, пути и средства повышения надёжности и безопасности управления КА.

1.3. Требования, предъявляемые к средствам отработки проведения динамических операций КА.

1.3.1. Требования к алгоритмическому обеспечению.

1.3.2. Требования к программно-моделирующему комплексу.

1.4. Постановка задач исследования.

1.5. Выводы по главе 1.

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

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

2.1.1. Системы координат и уравнения движения КА.

2.1.2. Невозмущенное движение.

2.1.3. Влияние нецентральности гравитационного поля Земли.

2.1.4. Влияние гравитационных полей Луны и Солнца.

2.1.5. Влияние давления солнечного света.

2.1.6. Расчет координат Солнца.

2.1.7. Расчет координат Луны.

2.1.8. Модель вектора тяги.

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

2.2.1. Методика изменения высоты перигея.

2.2.2. Методика изменения периода обращения.

2.3. Математическая постановка задачи оптимизации управления.

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

2.4.1. Стохастический подход.

2.4.1. Гарантирующий подход.

2.4.3. Совмещение задач удержания и выполнения динамической операции в окрестности опорной орбиты.

2.5. Терминальный алгоритм управления движением.

2.6. Выводы по главе 2.

Глава 3. Методика построения и технический облик программно-моделирующего комплекса для отработки средств проведения динамических операций КА.

3.1. Состав системы ориентации и управления движением.

3.2. Состав и структура программно-моделирующего комплекса.

3.3. Диспетчер событий и глобальных переменных.

3.3.1. Организация данных и дескрипторы.

3.3.2. Учет ограничений и моделирование случайных возмущений.

3.3.3. Программный интерфейс.

3.3.6. Пример.

3.4. Основные алгоритмы проведения расчётов в управляющем модуле и организации взаимосвязи между компонентами программного комплекса.

3.4.1. Учет частоты срабатывания бортовых систем КА.

3.4.2. Методика интегрирования уравнений движения. Определение текущего шага интегрирования.

3.4.3. Организация работы управляющего модуля с использованием диспетчера событий и глобальных переменных.

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

3.5.1. Архитектура базы данных ПМК.

3.5.2. Соглашения об элементах реляционной модели данных.

3.5.3. Реляционная модель данных.

3.5.4. Средства контроля доступа.

3.5.5. Организация связи с другими модулями программно-моделирующего комплекса.

3.5.6. Интерфейс пользователя базы данных.

3.6. Выводы по главе 3.

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

4.1. Использование программно-моделирующего комплекса для отработки средств выведения КА на высокоэллиптическую орбиту.

4.2. Использование программно-моделирующего комплекса для отработки средств выведения КА на круговую орбиту.

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

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

4.5. Выводы по главе 4.

Введение 2012 год, диссертация по авиационной и ракетно-космической технике, Старков, Александр Владимирович

Актуальность темы. В результате человеческой деятельности количество управляемых и неуправляемых космических объектов в околоземном пространстве с каждым годом будет расти, а вероятность столкновения космических объектов друг с другом - увеличиваться. Столкновения и разрушения космических аппаратов (КА) приводят к образованию новых крупных фрагментов и осколков. По приблизительным оценкам, в околоземном пространстве находится более девяти тысяч объектов с размерами более 10 сантиметров, десятки тысяч с размерами 1-10 сантиметров и сотни тысяч с размерами менее 1 сантиметра. В последнее время «космический мусор» стал представлять все большую опасность для вновь запускаемых КА, но приемлемые способы борьбы с ним пока не разработаны. Поэтому потребность в обеспечении безопасного проведения динамических операций в околоземном пространстве очевидна. Уменьшение угрозы столкновения КА с другими космическими объектами представляет собой сложную организационно-техническую задачу, решению одной из подзадач которой посвящена данная работа.

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

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

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

Актуальность и новизна работы следует из анализа развития специализированного программного обеспечения, как на мировом рынке, так и во внутрикорпоративном сегменте (продукты Satellite Tool Kit, Bernese Software, «Orbit Determination Toolkit», «Free Flyer» и др.), который выявил тенденцию их развития по пути универсализации. Производитель стремится включить в свой продукт все больше функций и инструментов решения задач. Это неизбежно приводит к значительному усложнению и снижению эффективности их применения конечным пользователем. Ключевая цель обеспечения функциональной полноты программных средств остается нереализованной. Разумным шагом на пути к преодолению этого недостатка является специализация решений на основе концепции открытой архитектуры, которая предлагает строить ПМК из следующих компонент:

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

2) единого, стандартизованного языка высокоуровневого описания объектов моделирования, входной и выходной информации, исполнительной макропрограммы с поддержкой механизмов межпрограммного взаимодействия;

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

Анализ публикаций на данную тему и результаты работ, выполненных авторами: В А. Бартеневым, М.Н. Красилыциковым, А А. Лебедевым, J1.H. Лысенко, В.В. Малышевым, В.Н. Почукаевым, М.Ф. Решетневым, Г.М. Чернявским, Р. Шенноном, позволил определить основные особенности построения алгоритмов управления космическим аппаратом при проведении динамических операций, а также пути исследований, направленных на определение общих принципов и методов решения поставленной задачи. В работах этих авторов излагались теоретические основы и методы проектирования КА, исследовались проблемы навигационно-баллистического обеспечения полета, модели движения и принципы управления КА, работа целевой аппаратуры и служебных систем. Вместе с тем остается еще ряд задач, связанных с созданием унифицированных средств автоматизации отработки динамических операций КА.

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

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

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

Целью диссертационной работы является совершенствование алгоритмического и программного обеспечения за счет синтеза алгоритмов управления КА с учетом требований безопасности и разработке методики построения программно-моделирующего комплекса (ПМК) для отработки динамических операций.

Для достижения поставленной цели решаются научно-технические задачи:

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

2. Разрабатываются математические модели бортовых систем КА и внешней среды.

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

4. Определяются принципы построения и технический облик программно-моделирующего комплекса.

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

Методы исследования. В качестве методологической основы используется системный подход. На нем основывается принцип создания единой базовой подсистемы - типового унифицированного программного средства, методический и алгоритмический аппарат которого позволяет определять движение широкого множества КА ближнего и среднего космоса, гибко реагируя при этом на изменение параметров орбит КА, допустимых значений показателей качества и различных моделей измерений. Основными методами исследования, используемыми в работе, являются методы динамики полета, теории управления, статистические методы обработки данных, методы оптимального управления и обработки информации. При программной реализации математического обеспечения используются методы объектно-ориентированного программирования и мультизадачность операционных систем Windows, Linux, системы управления реляционными базами данных, а также средства обеспечения доступа в сетях Интернет и Интранет.

Научная новизна результатов. В работе получены следующие результаты, обладающие новизной и научной значимостью:

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

2. Определены математические модели проведения орбитальных коррекций К А при выполнении динамических операций.

3. Сформирована адаптированная к требованиям программно-моделирующего комплекса (ПМК) математическая модель бортовых систем, внешних воздействий и контура управления.

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

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

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

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

Результаты диссертационной работы используются в учебном процессе кафедры «Системный анализ и управление» МАИ и рабочей деятельности ФГУП ЦНИИмаш, ФГУП «НПО им. С.А. Лавочкина». Имеются соответствующие акты внедрения от МАИ, ФГУП ЦНИИмаш, ФГУП «НПО им. С.А. Лавочкина».

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

Основные положения, выносимые на защиту.

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

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

3. Методика разработки базы данных программно-моделирующего комплекса.

4. Результаты использования программно-моделирующего комплекса для отработки безопасного проведения динамических операций на высокоэллиптических и круговых геостационарных орбитах (ГСО). Апробация работы и публикации. Результаты работы докладывались и получили одобрение на научно-технических конференциях: 16-ой Международной научной конференции «Системный анализ, управление и навигация» (Евпатория, 2011 г.), 17-ой Международной научной конференции «Системный анализ, управление и навигация» (Евпатория, 2012 г.)

Содержание работы.

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

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

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

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

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

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

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

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

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

Заключение диссертация на тему "Синтез алгоритмов управления космическими аппаратами с учетом требований безопасности проведения динамических операций"

4.5. Выводы по главе 4.

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

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

• вывод КА на расчетную высокоэллиптическую и круговую орбиты как штатными алгоритмами выведения, так и с помощью терминального метода;

• выполнение маневров на высокоэллиптической орбите;

• выполнение маневров на круговой орбиты.

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

3. Для повышения точности выведения на заключительном этапе необходимо использование терминальных алгоритмов наведения с обратной связью. Они позволяют компенсировать разбросы начальных условий и других случайных параметров за счет использования информации о реальном векторе состояния, получаемой либо с наземных, либо бортовых навигационных средств, в том числе с использованием спутниковой навигационной системы ГЛОНАСС.

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

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

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

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

Заключение

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

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

1. Сформированы требования к алгоритмическому и программному обеспечению средств отработки динамических операций.

2. Предложены математические модели бортовых систем КА и внешней среды.

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

Практическая значимость работы заключается в том, что автор

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

2. Создал программно-моделирующий комплекс для отработки средств проведения динамических операций КА.

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

1. Вывод КА на расчетную высокоэллиптическую и круговую орбиты как штатными алгоритмами выведения, так и с помощью терминального метода.

2. Выполнение маневров на высокоэллиптической орбите.

3. Выполнение маневров на геостационарной орбите.

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

Библиография Старков, Александр Владимирович, диссертация по теме Динамика, баллистика, дистанционное управление движением летательных аппаратов

1. Liou J.-C., Matney М.J., Anz-Meador P.D., Kessler D., Janscn M., Theall J.R. The New NASA Orbital Debris Engineering Model ORDEM2000. NASA/TP-2002-210780, May, 2002

2. Moknight D.S. Lorenson G. Collision matrix for low earth orbit satellites/ Journal of Spacecraalf and Rockets -26-№2-1989

3. Walker R., Martin C., Stokes H., Wilkinson J., Sdunnus H., Hauptmann S., Beltrami P., Klinkrad HL Update of the ESA Space Debris Mitigation Handbook. Executive Summary July 2002

4. Аксенов Е.П. Специальные функции в небесной механике. М.: Наука. Гл. ред. физ.-мат. лит., 1986.

5. Алексеев К.Б., Бебенин Г.Г., Ярошевский В.А. Маневрирование космических аппаратов. -М.: Машиностроение, 1970.

6. Аппазов Р.Ф., Сытин О.Г. Методы проектирования траекторий носителей и спутников Земли. -М.: Наука, 1987.

7. Артюхин Ю.П., Каргу Л.И., Симаев В.Л. Системы управления космических аппаратов, стабилизированных вращением. М.: Наука, 1979.

8. Бебенин Г.Г., Скребушевский Б.С, Соколов Г.А. Системы управления полетом космических аппаратов. -М.: Машиностроение, 1978.

9. Гродзовский Г.Л., Иванов Ю.Н., Токарев В.В. Механика космического полета с малой тягой. М.: Наука, 1966.

10. Гущин В.Н. Основы устройства космических аппаратов: Учебник для вузов.-М. Машиностроение, 2003.

11. Дейт К.Дж. Введение в системы баз данных = Introduction to Database Systems. — 8-е изд. —М.: «Вильяме», 2006.

12. Иванова Г.С. Технология программирования: Учебник для вузов. / Г.С. Иванова. М.: Изд-во МГТУ им. Н.Э. Баумана, 2002.

13. Касаев К.С. Новые наукоемкие технологии в технике: Энциклопедия. Том 10. Системный подход к сложным техническим объектам. М., 1997.

14. Когаловский М.Р. Перспективные технологии информационных систем. — М.: ДМК Пресс; Компания АйТи, 2003.

15. Когаловский М.Р. Энциклопедия технологий баз данных. — М.: Финансы и статистика, 2002.

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

17. Лебедев A.A. Введение в анализ и синтез систем: учебное пособие. М.: Изд-воМАИ, 2001.

18. Лебедев A.A. Ведение в анализ и синтез систем. Учебное пособие. Изд-во МАИ, 2001.

19. Лебедев A.A., Аджимамудов Г.Г. Основы синтеза систем летательных аппаратов: Учебник. -М., МАИ, 1996.

20. Лебедев A.A., Аджимамудов A.A., Баранов В.Н., Бобронников В.Т., Кр асильщиков М.Н., Малышев В.В., Нестеренко О.П. Основы синтеза систем летательных аппаратов. -М.: Машиностроение, 1996.

21. Лебедев A.A., Баранов В.Н., Бобронников В.Т., Красильщиков М.Н., Малышев В.В., Нестеренко О.П. Основы синтеза систем ЛА. М.: Машиностроение, 1987.

22. Лебедев A.A., Бобронников В.Т., Красильщиков М.Н., Малышев В.В. Статистическая динамика управляемого полета. — М.: Машиностроение, 1978.

23. Лебедев A.A., Бобронников В.Т., Красильщиков М.Н., Малышев В.В.

24. Статистическая динамика и оптимальное управление летательных аппаратов. -М.: Машиностроение, 1985.

25. Лебедев A.A., Красильщиков М.Н., Малышев В.В. Оптимальное управление движением космических летательных аппаратов. — М.: Машиностроение, 1974.

26. Лебедев A.A., Нестеренко О.П. Космические системы наблюдения. Синтез и моделирование -М.: Машиностроение, 1991.

27. Лебедев A.A., Соколов В.Б. Встреча на орбите. — М.: Машиностроение, 1969.

28. Лысенко Л.Н. Наведение и навигация баллистических ракет: Учеб.пособие. -М.: Изд-во МГТУ им. Н.Э. Баумана, 2007.

29. Максимов Г.Ю. Теоретические основы разработки космических аппаратов. М.: Наука. Главная редакция физико-математической литературы, 1980.

30. Малышев В.В. Методы оптимизации в задачах системного анализа и управления: Учебное пособие. -М.: Изд-во МАИ-ПРИНТ, 2010.

31. Малышев В.В., Бобронников В.Т., Карп К.А., Нестеренко О.П., Федоров A.B. Пакеты прикладных программ для оптимизации и управления авиационно-космичесими комплексами. В сб. «Системный анализ в технике-3». М.: Изд-во МАИ, 1994.

32. Малышев В.В., Карп К.А., Федоров A.B. Автоматизация решения вероятностных задач. М.: Изд-во МАИ, 1995.

33. Малышев В.В., Кибзун А.И. Анализ и синтез высокоточного управления летательным аппаратом. -М.: Машиностроение, 1987.

34. Малышев В.В., Кузнецов А.К., Федоров A.B. Достаточные условия оптимальности в задачах управления движением. М.: Изд-во МАИ, 1990.

35. Малышев В.В., Старков A.B., Федоров A.B. Методика актуализации банков данных объектов космического назначения. // «Вестник Московского Авиационного Института». Том 19. Выпуск №4. М.: Изд-во МАИ, 2012г.

36. Малышев В.В., Старков A.B., Федоров A.B. Методика разработки программно-моделирующего комплекса для отработки средств проведения динамических операций космических аппаратов. Электронный журнал «Труды МАИ». Выпуск № 57, 2012 г.

37. Малышев В.В., Старков A.B., Федоров A.B. Программно-моделирующий комплекс для отработки средств проведения динамическихопераций космических аппаратов. Вестник компьютерных и информационных технологий. №9, 2012 г.

38. Малышев В.В., Федоров A.B. Минимаксный подход в задачах определения гарантирующей стратегии управления. Доклад на 4-й конференции.по динамике и управлению. Труды МАИ, 1976.

39. Малышев В.В., Федоров A.B. Управление движением спутников космической системы // сборник трудов II Международной конференции «Малые спутники. Новые технологии. Минитюаризация. Области эффективного применения в XXI веке». М.: ЦНИИМАШ, 2000.

40. Назаренко А.И., Скребушевский Б.В. Эволюция и устойчивость спутниковых систем. -М.: Машиностроение, 1981.

41. Назваренко А.И. Построение высотно-широтного распределения объектов в околоземном космическом пространстве. Проблема загрязнения космоса (космический мусор). Москва, Космосинформ, 1993.

42. Назваренко А.И., Чернявский Г.М. Моделирование загрязнения околоземного пространства. Столкновения в околоземном пространстве (космический мусор). Москва, Космосинформ, 1995.

43. Нариманов Г.С., Тихонравов М.К. Основы теории полета космических аппаратов. — М.: Машиностроение, 1972.

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

45. Салмин В.В. Оптимизация космических перелетов с малой тягой. М.: Машиностроение, 1987.

46. Скребушевский Б.С. Управление полетом беспилотных космических аппаратов. — М.: «Владмо», 2003.

47. Соловьев В.А., Лысенко Л.Н., Любинский В.Е. Управление космическими полетами. Учебное пособие в 2-х частях. М.: Изд-во МГТУ им.Н.Э.Баумана, 2009.

48. Сухой Ю.Г. Коррекции орбит геостационарных спутников: В 3-х частях. Часть 1. Особенности управления орбитальным движением и возмущения орбит геостационарных спутников: Пособие для специалистов. М: Издательство «Спутник+», 2011.

49. Федоров A.B. Применение минимаксного подхода в задачах оптимизации управления дискретными системами. В сб. «Новые методы высокоточного управления и оценивания и управления JIA». М.: Изд-во МАИ, 1986.

50. Чернявский Г.М., Бартенев В.А. Орбиты спутников связи. М., «Связь», 1978.

51. Чернявский Г.М., Бартенев В.А., Малышев В.А. Управление орбитой стационарного спутника. — М.: Машиностроение, 1984.

52. Эльясберг П.Е. Введение в теорию полета искусственных спутников Земли. 2-е изд. - М.: Эдиториал УРСС, 2011.