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

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

Захваткин Михаил Витальевич

ОПРЕДЕЛЕНИЕ И ПРОГНОЗИРОВАНИЕ

ПАРАМЕТРОВ ДВИЖЕНИЯ КОСМИЧЕСКОГО АППАРАТА С УЧЕТОМ ВОЗМУЩЕНИЙ, ВЫЗВАННЫХ РАБОТОЙ БОРТОВЫХ СИСТЕМ

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

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

Москва — 2014

2 О ФЕБ 2014

005545338

005545338

Работа выполнена в Федеральном государственном бюджетном учреждении науки Институте прикладной математики им. М.В. Келдыша Российской академии наук.

Защита состоится 25 марта 2014 года в

П00

часов на заседании диссертационного совета Д 002.024.01, созданного на базе Института прикладной математики им. М.В. Келдыша РАН, расположенного по адресу: 125047, Москва, Миусская пл., 4.

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

Автореферат разослан 11 февраля 2014 года.

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

Научный руководитель:

д.ф.-м.н., профессор Сазонов Виктор Васильевич

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

д.т.н., Хуторовский Захарий Наумович, начальник сектора ОАО «МАК ВЫМПЕЛ»

д.т.н., Назиров Равиль Равильевич, зам. директора ФБГУН ИКИ РАН

Федеральное государственное унитарное предприятие «Центральный научно-исследовательский институт машиностроения» (ФГУП ЦНИИмаш)

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

д.ф.-м.н.

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

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

Актуальность работы

Фундаментальные космические исследования являются приоритетным направлением Федеральной космической программы России до 2030 г. В настоящее время на рабочей орбите находится космический аппарат (КА) «Спектр-Р». В текущем десятилетии планируется развертывание внеатмосферных обсерваторий «Спектр-РГ», «Спектр-УФ», «Гамма-400», «Спектр-М», а также запуск орбитальных зондов и посадочных аппаратов для исследования Луны по проекту «Луна-Глоб». Задачи управления движением и навигационной привязки проводимых научных экспериментов предъявляют высокие требования к точности определения и прогнозирования параметров движения этих КА.

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

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

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

го в июле 2011 г. в соответствии с проектом создания высокоорбитального наземно-космического радиоинтерферометра «Радиоастрон». Большое удаление КА от Земли обеспечивает более высокое разрешения по сравнению с наземными интерферометрами и космическими интерферометрами на низких орбитах. Реализовать это преимущество можно только при условии высокоточного определения параметров базы интерферометра, а следовательно, параметров движения КА. Основные требования к точности определения параметров движения КА составляют по положению ± 600 м, по скорости ± 2 см/с, по ускорению ±10-8 м/с2. Выполнить эти требования невозможно без учета возмущений, обусловленных работой бортовых систем аппарата. Возникла необходимость разработки параметризованной модели давления солнечной радиации, учитывающей ориентацию аппарата и свойства его поверхности, а также методики определения параметров движения аппарата путем совместной обработки внешнетраекторных и телеметрических данных [1].

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

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

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

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

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

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

Научную новизну работы составляют:

• Параметризованная модель силы и момента светового давления, использованная для повышения точности определения параметров орбиты КА «Спектр-Р».

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

• Метод долгосрочного прогнозирования движения КА с учетом прогноза изменения кинетического момента и разгрузок ЭМИО.

• Математическая модель видимого блеска КА «Спектр-Р», согласующаяся с моделью влияния светового давления на движение аппарата.

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

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

Практическая значимость

Модели и методики, разработанные в диссертации, использованы в рамках баллистико-навигационного обеспечения КА «Спектр-Р» — основного элемента международного проекта «Радиоастрон». Результаты, полученные в работе, позволили решить ряд важных прикладных задач, среди которых:

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

2. Долгосрочное прогнозирование движения КА «Спектр-Р», позволяющее избежать больших продольных ошибок в расчетном движении КА и соответствующих им временных ошибок при планировании работы наземно-космического интерферометра.

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

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

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

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

2. Метод использования скоростей вращения маховиков ЭМИО в качестве измерений внешнего возмущающего момента при уточнении параметров движения КА. Расчетная величина возмущающего момента зависит от уточняемых параметров модели светового давления.

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

4. Математическая модель видимого блеска КА «Спектр-Р», которая согласуется с моделью влияния светового давления на движение аппарата.

Апробация работы и публикации

Результаты работы изложены в 5 печатных работах в изданиях, рекомендованных ВАК [1—5].

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

• V Международная научная конференция «Наблюдение околоземных космических объектов». Исследования первых оптических наблюдений КА «Спектр-Р», 9-12 ноября 2011 г., Москва, Россия.

• Математические модели и моделирование в лазерно-плазменных процессах и передовых научных технологиях ЬРрМЗ. Высокоточное определение орбиты КРТ миссии «Радиоастрон», 28 мая - 2 июня 2012 г., Петро-вац, Черногория.

• Radioastron International Science Council 2012. Possibility of determination of Earth's and Moon's gravitational fields from Radioastron orbital data, 18-20 июня 2012 г., Пущино, Россия.

• Radioastron International Science Council 2013. Orbit determination of the Radioastron mission, 20-21 июня 2013 г., Москва, Россия.

• Семинар сектора № 2 отдела № 5 Института прикладной математики им. М.В. Келдыша РАН. Руководитель: д.ф.-м.н. А.Г. Тучин, 19 июня, 2 ноября 2013 г., Москва, Россия.

• Расширенный семинар отдела .X5 5 «Механика космического полета и управление движением» Института прикладной математики им. М.В. Келдыша РАН. Руководитель: проф. Ю.Ф. Голубев, 5 декабря 2013 г., Москва, Россия.

• НТС подсекции № 4 секции № 4 НТС ФГУП ЦНИИмаш. Руководитель: И.И. Олейников, 10 декабря 2013 г., Королев, Россия.

Структура и объем диссертации

Диссертация состоит из введения, четырех глав, заключения и библиографии. Общий объем диссертации составляет 122 страницы, из них 117 страниц текста, включая 46 рисунков и 13 таблиц. Библиография включает 44 наименования на 5 страницах.

Содержание работы по главам

Введение

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

Глава 1. Модель движения

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

В разделе 1.2 описывается используемая в работе модель движения центра масс КА. Штатный полет происходит по пассивной траектории, прерывающейся сеансами проведения разгрузок ЭМИО, сопровождающихся включением двигателей стабилизации (ДС). Разгрузки происходят несколько раз за сутки и представляют собой попеременное включение ДС с целью погасить суммарный кинетический момент аппарата вместе с маховиками. Процесс длится 1-3 минуты, за которые происходит несколько десятков включений ДС. Работа двигателей телеметрируется, для каждого включения во время г-й разгрузки известно время длительность работы т^ и расход топлива Amf\ Номинальная тяга ДС определяется длительностью включения, поэтому приращение скорости КА в результате одного включения ДС является известной функцией длительности включения и расхода топлива Av<?} = Avlj){Am{p ,rlj)). Суммарное приращение скорости КА в результате г-й разгрузки представляется в виде одного вектора Avi, обозначаемого импульсом разгрузки. Средневзвешенный момент времени приложения импульса и его измеренное значение определяются выражениями

у-ЛГ . (j),(j) N

Е^л^ 'it

и)

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

На пассивных участках траектории учитываются следующие силы, действующие на движение КА: гравитация Земли, содержащая центральную и нецентральную части геопотенциала, гравитация Луны, Солнца и планет; изменение гравитационного поля вследствие деформации Земли под действием притяжения Луны и Солнца (т.н. «твердые приливы»); давление солнечного излучения; атмосферное торможение (на участках ниже 1.5 тыс. км над поверхностью Земли); влияние альбедо Земли, а также дополнительное возмущающее ускорение, обеспечивающее учет эффектов общей теории относительности. Гравитационное поле Земли представлено разложением геопотенциала по сферическим функциям в соответствии с моделью EGM-96 до гармоники 75x75. Для получения координат Луны, Солнца и планет используются таблицы, основанные на теории движения DE421. Для описания влияния приливных сил используется модель, учитывающая деформацию Земли в направлении возмущающего тела (Луны или Солнца) в форме гармоники второго порядка. Расчет плотности атмосферы осуществляется на основе модели, рекомендованной ГОСТ Р 25645.166-2004. Учет переизлученного от Земли света производится согласно модели разбиения поверхности планеты на 160 участков с постоянными коэффициентами альбедо. Поправки для учета эффектов общей теории относительности вычислялись с использованием модели движения в центральном невращающемся гравитационном поле Земли.

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

Таким образом была построена модель силы и момента светового давления, действующих на КА «Спектр-Р». Для этого поверхность аппарата разбита на три части: космический радиотелескоп (КРТ), центральный блок и панели солнечных батарей [2]. Описание модели поверхности К А дается в разделе 1.3. Для описания свойств поверхности КРТ и центрального блока используются коэффициенты а\ и fii, поскольку эти элементы К А покрыты одинаковой многослойной изоляцией. Отражающая способность панелей солнечных батарей задается коэффициентом второй коэффициент не вводится. Из-за того, что панели всегда ориентированы практически нормально

Рис. 1. Модель поверхности КА «Спектр-Р», пример возникающей тени и конечноэлементное представление поверхности КРТ

по отношению к солнечному потоку, направление действующей силы определено с достаточной точностью, а параметры светового давления скоррелиро-ваны. Для определенности полагается цъ = 0. Для учета тени, возникающей на КРТ, поверхность антенны разбивается на множество элементов (рис. 1), центр каждого из которых проверялся на затенение от центрального блока и панелей солнечных батарей. Разработанная модель описывает силу светового давления, действующую на КА «Спектр-Р», при помощи трех параметров qi, Mi и

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

По известному движению аппарата вокруг центра масс и скоростям вращения ДМ из уравнений движения получается оценка возмущающего момента, зависящего, в том числе, от неизвестных коэффициентов светового давления. Рассматривается случай поддержания аппаратом постоянной ориентации в инерциальном пространстве, преобладающий в движении КА «Спектр-Р» и актуальный для будущих космических обсерваторий. Если на интервале (ii, ¿2) движение происходит вдали от гравитирующих тел, накап-

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

-----= М5р(8,аь .. .,ат,т,.. .,Цт), (2)

12 — ъ 1

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

момент инерции, скорость вращения г'-го ДМ, Мяр — момент сил светового давления, я — единичный вектор направления солнечных лучей в связанной с аппаратом системе координат.

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

N Лдс

~ X] Г3 х > 0. 3 = • • •> мдс, (3)

!=1 ¿=1

где г^ — радиус вектор центра сопла ДС в связанной с аппаратом системе, е,- — направление истечения рабочего тела ДС. В общем случае выражение (3) без введения дополнительных условий не обеспечивает однозначности перевода накопленного кинетического момента в импульс разгрузки. Эти дополнительные условия рассмотрены в главе 4.

Глава 2. Моделирование траекторных измерений

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

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

Разделы с 2.3 по 2.7 посвящены построению расчетных значений траекторных измерений штатной радиотехнической системы, беззапросных до-

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

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

К = аЦЕ - е • ет) + а2уе ■ ет, (4)

где Е — единичная матрица, е — направление тяги двигателя, сг„ — ошибка определения величины вектора, о^ — ошибка в направлениях, ортогональных е. В случае если направления всех Д^'^ совпадают, ковариационная матрица импульса разгрузки Д^ представляется в более простом виде. Такая ситуация имеет место в случае КА «Спектр-Р», который компенсирует накопленный кинетический момент вдоль одной из осей связанной системы координат моментным образом, а в плоскости двух других осей — силовым образом при помощи ДС с одинаковым направлением тяги.

Раздел 2.9 содержит описание измерений внешнего возмущающего момента и соответствующих им ковариационных матриц. Скорость изменения кинетического момента системы ДМ равна внешнему возмущающему моменту, если К А не вращается (рис. 2). Вдали от Земли и других гравитирующих тел основным источником возмущающего момента является световое давление. Измеренное значение момента в этом случае соответствует левой части выражения (2), расчетное значение в правой части зависит от коэффициентов светового давления.

Глава 3. Определение орбиты

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

Рис. 2. Телеметрия КУДМ КА «Спектр-Р». Кинетический момент ДМ линейно возрастает на интервалах поддержания аппаратом постоянной ориентации.

КА «Спектр-Р», построенный по реальным траекторным и телеметрическим данным.

Постановка задачи определения параметров движения дается в разделе 3.1. В условии задачи световое давление, действующее на КА, зависит от неизвестных параметров {q;,^}^. На интервале уточнения движение аппарата возмущено импульсами разгрузок ДМ {¿.¡, Av,}™=1, помимо этого в течение интервала аппарат N раз находился в неизменной ориентации в инерциальной СК вдали от притягивающих тел. Набор уточняемых параметров включает в себя вектор состояния КА в начальный момент времени внутри интервала уточнения, 2т параметров светового давления и векторы импульсов разгрузок

Q = {Xo(io),ai,. ■ - ,ат, fii,..., цт, Avb ..., Avn}.

Уточнение параметров производится при помощи траекторных измерений Ф и данных, полученных в результате обработки телеметрической информации, включающих измерения импульсов разгрузок {Av?}"=1 и рассогласования измеренных и расчетных величин возмущающего момента полученные из (2). Определяется вектор параметров Q. минимизирующий функционал [3; 4]

N

Ф =(Ф0 - ФС)ТР(Ф0 - Фс) + £

i=l

П

+ £(Ду° - Avj)TPj(Av° - Av,:), (5)

i=i

где индексы о и с вектора траекторных измерений обозначают измеренные и расчетные значения соответственно, Р — весовая матрица траекторных измерений, Рг и Р®р — весовые матрицы измерений импульсов разгрузок и возмущающих моментов, полученные в разделах 2.8 и 2.9.

В разделе 3.2 описывается алгоритм поиска неизвестного вектора С^. Определение параметров осуществляется при помощи последовательных приближений, полученных методом обобщенных касательных Ньютона. Поправки к приближению вектора параметров С2 на каждом шаге определяются решением системы нормальных уравнений

матрица А и правые части В которой выражаются через входящие в (5) рассогласования, весовые матрицы и частные производные <9ФС/<9(3, д^/дЦ и д^/дС^.

Частные производные траекторных измерений Фс рассчитываются при помощи производных измерений по текущему вектору состояния, описанных во второй главе, и производных текущего вектора состояния по уточняемым параметрам <9ФС/<9(3 = ЭФС/<9Х ■ дХ/дЦ. Последние получаются интегрированием дополнительных 36 4- 12т уравнений в вариациях для определения производных по Х0 и параметрам светового давления. Производные по импульсам разгрузки рассчитываются следующим образом [3]:

где Ё - последние 3 столбца единичной матрицы 6x6.

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

В разделе 3.3 разработанная методика используется для уточнения параметров движения КА «Спектр-Р» на двух интервалах: с 20.02.2013 по 10.04.2013 и с 10.04.2013 по 30.05.2013. Для оценки эффективности разработанного метода орбита аппарата на выбранных интервалах уточняется четырьмя способами, различающимися моделями движения КА и набором уточняемых параметров.

1. В первом варианте возмущение от светового давления рассчитывается по упрощенной модели и зависит только от одного уточняемого параметра ■к. Разгрузки ДМ не учитываются.

2. Во втором варианте возмущение от светового давления также рассчитывается по упрощенной модели с одним уточняемым параметром. Разгрузки ДМ учитываются, но не уточняются, а используются измеренные значения.

А Д<2 = В,

(6)

№ Световое давление Разгрузки 1Т2 Дг, км Ау, мм/с

1 Упрощенное, х Не учитываются 12.43677 9.18588 71.71 288.1

2 Упрощенное, х Учитываются 4.72914 6.78832 36.76 113.3

3 Сложное, ai.ui, <*2 Учитываются 1.20896 0.63767 7.57 8.9

4 Сложное, оц, , а2 Уточняются 0.36210 0.31607 0.21 2.3

Таблица 1. Безразмерные СКО траекторных измерений, полученные в результате

уточнения орбиты КА «Спектр-Р» на интервалах 20.02.2013-10.04.2013 (01) и 10.04.2013-30.05.2013 (02)- Рассогласование полученных решений в момент времени 10.04.13 00:00:00 11Т по положению (Дг) и по скорости (Дг>).

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

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

При уточнении ошибка измерения дальности полагалась равной 100 м, ошибка запросных измерений радиальной скорости была установлена на уровне 10 мм/с, беззапросных измерений — 5 мм/с. Предполагалось, что ошибки оптических изхмерений прямого восхождения и склонения равны одной угловой секунде и не имеют между собой корреляции. Результаты согласования траекторных измерений приведены в таблице 1. Как видно из таблицы, постепенное улучшение модели движения приводит к улучшению согласования измерений; учет сложного светового давления и уточнение импульсов разгрузок улучшает согласование измерений в 29-34 раза на выбранных интервалах уточнения. Рассогласования измеренных и расчетных величин траекторных измерений, полученные с использованием сложной модели светового давления и уточнением импульсов разгрузок, приведены на рисунке 3.

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

Для каждого из способов уточнения полученные на двух интервалах параметры светового давления согласуются между собой. В случае использования сложной модели светового давления и уточнении импульсов разгрузок коэффициенты светового давления уточнились в следующих интервалах: а\ = 0.86 — 0.87, /XI = 0.08 — 0.13, аг = 0.07 — 0.09, что хорошо соотносится со

Наклонная дальность. О-С

——* ** ш т чФН^ ** j I

Склонение, О-С

+

1

+ 3,,. 1... - * ■

+

Склонение, О-С

А

+

Н Китаб ЕЗ Монды Ужгород 1 1 Звенигород

ЕИЭ Евпатория fcs! Краснодар Е=1 Благовещенск

I Китаб Е I Евпатория I

I Н06, Мейхилп I ) Краснодар I

I Ужгород 1

] 062, Австралия

Рис. 3. Рассогласования траекторных измерений КА «Спектр-Р»

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

Глава 4. Прогнозирование параметров движения

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

Раздел 4.1 посвящен расчету будущих возмущений движения аппарата, наибольшая неопределенность которых обусловлена разгрузками. Для учета будущих разгрузок прогнозируются времена их проведения и величины соответствующих импульсов. Из накопленного ^ДМ кинетического момента при помощи (3) и дополнительного условия Yli=i Pi -S> min рассчитывается величина импульса оптимальной разгрузки. Рассчитанные таким образом

импульсы хорошо согласуются с измеренными значениями, полученными из телеметрии ДС (см. рис. 4).

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

Проверка описанной методики прогноза с учетом возмущений от разгрузок ДМ проводится на КА «Спектр-Р» и описывается в разделе 4.2. Прогноз движения осуществляется на временной интервал с 10.04.2013 по 30.05.2013. При этом используются параметры, полученные в результате уточнения на интервале с 20.02.2013 по 10.04.2013 (см. главу 3). По известным параметрам светового давления и ориентации КА на интервале прогнозирования было рассчитано накопление кинетического момента ДМ и будущие разгрузки. В Рис. 4. Соотношение величин одном случае прогнозировались вы-

импульсов оптимальных и фактических нужденные разгрузки, возникающие разгрузок 08.2011-09.2013 При выходе кинетического момента из

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

Построено три варианта прогноза. Первый рассчитывался без учета разгрузок и сложного светового давления. Второй и третий — с использованием вынужденных и регулярных разгрузок соответственно. На рисунке 5 изображено отклонение прогнозного положения КА от фактического в продольном направлении, составляющее основную часть ошибки. Максимальное отклонение прогноза, рассчитанного по упрощенной модели, не использующей результаты этой работы, составило 387 км, прогноза с использованием вынужденных разгрузок — 93 км, с использованием регулярных разгрузок — 52 км. Таким образом, предложенная методика позволила улучшить точность прогноза положения КА более чем в 7 раз на исследуемом интервале.

Условие проведения разгрузки при достижении критического кии. момента также имело место, но на практике не достигалось.

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

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

Видимый блеск определяется плотностью светового потока, отраженного от поверхности КА в сторону наблюдателя. Подход, используемый в динамической модели, недостаточно точно описывает распределение зеркально отраженного света. В разработанной модели блеска интенсивность зеркально отраженного света пропорциональна соэ^ в', где в' — угол между направлением на наблюдателя и направлением зеркального отражения [5]. Отражение также разделятся на зеркальное и диффузное, их долевое соотношение описывается коэффициентами а и /г поверхности.

Модель видимого блеска построена для КА «Спектр-Р». Разработанная модель обеспечивает в 1.5 раз лучшее согласование попутных фотометрических измерений, чем упрощенная модель сферы, не учитывающая ориентацию аппарата. Такое согласование достигается с использованием коэффициентов а и ¡х, полученных в ходе уточнения орбиты. При использовании уточненных значений ^ = 0.57219, = 0.17931 согласование измерений улучшается в 1.9 раз по сравнению с простой моделью.

Заключение

Отмечается:

На основе разработанной модели силы и момента светового давления, учитывающей форму и характеристики поверхности КА, была решена задача определения орбиты КА «Спектр-Р» и навигационной привязки интерферо-метрических измерений с ошибками в 20-30 раз меньшими, чем при исполь-

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

Учет ориентации КА на интервале прогнозирования и предстоящих разгрузок ДМ позволяет существенно повысить точность прогноза движения аппарата. В приведенном в качестве примера случае ошибка прогнозирования положения КА «Спектр-Р» на 50 суток вперед уменьшилась в 7 раз (с 387 км до 52 км).

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

Список литературы

1. Кардахиев Я., Хартов В., ..., Захваткин М. [и др.] Космическая миссия «Радио-астрон». Первые результаты // Вестник НПО им. С.А.Лавочкина. — 2012. — Т. 3, № 14. - С. 4-21. - ISSN 2075-6941.

2. Кардашев Н., Хартов В., ..., Захваткин М. [и др.] "Радиоастрон" - телескоп размером 300 000 км: основные параметры и первые результаты наблюдений // Астрономический Журнал. — 2013. — Т. 90, N« 3. — С. 179—222. — DOI: 10.7868/ S000462991303002X.

3. Боровин Г., Захваткин М., Степанъянц В., Тучин А. [и др.] Определение параметров орбиты и маневра космического аппарата при заданном времени приложения импульса // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. — 2012. — Спец. выпуск № 4 «Математическое моделирование». — С. 7G—86.

4. Боровин Г., Захваткин М., Степанъянц В., Тучин А. [и др.] Идентификация маневров, выполненных двигателями малой тяги космического аппарата // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. — 2012. — Спец. выпуск № 3 «Математическое моделирование». — С. 27—36.

5. Захваткин М. Моделирование видимого блеска космического аппарата «Спектр-Р» для планирования астрометрических наблюденнй // Наука и образование. — 2013. — Май. — № 3. — С. 10. — DOI: 10.7463/0513.0571011.

Подписано в печать 16.01.2014. Формат 60x84/16. Усл. печ. л. 1,3. Тираж 100 экз. Заказ П-60. ИПМ им. М.В. Келдыша РАН. 125047, Москва, Миусская пл., 4

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ УЧРЕЖДЕНИЕ НАУКИ ИНСТИТУТ ПРИКЛАДНОЙ МАТЕМАТИКИ ИМ. М.В.КЕЛДЫША РОССИЙСКОЙ АКАДЕМИИ НАУК

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

04201456421

ЗАХВАТКИН МИХАИЛ ВИТАЛЬЕВИЧ

ОПРЕДЕЛЕНИЕ И ПРОГНОЗИРОВАНИЕ ПАРАМЕТРОВ ДВИЖЕНИЯ КОСМИЧЕСКОГО

АППАРАТА С УЧЕТОМ ВОЗМУЩЕНИЙ, ВЫЗВАННЫХ РАБОТОЙ БОРТОВЫХ СИСТЕМ

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

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

Научный руководитель д.ф.-м.н., профессор Сазонов В. В.

Содержание

Введение 4

1 Модель движения 12

1.1 Системы координат и времени...................12

1.2 Движение центра масс........................16

1.2.1 Модель гравитационных возмущений ...........17

1.2.2 Модель негравитационных возмущений..........22

1.2.3 Сбалансированность модели.................29

1.3 Модель поверхности К А.......................30

1.4 Движение вокруг центра масс ...................34

1.4.1 Возмущающие моменты...................36

1.4.2 Режим движения с неизменной ориентацией.......38

1.4.3 Разгрузка маховиков.....................40

2 Моделирование траекторных измерений 43

2.1 Решение светового уравнения....................43

2.2 Задержки распространения в среде ................45

2.3 Радиотехнические измерения дальности..............46

2.4 Радиотехнические измерения радиальной скорости........49

2.5 Беззапросные измерения радиальной скорости..........50

2.6 Лазерные измерения дальности...................53

2.7 Оптические астрометрические измерения.............55

2.8 Измерения импульсов разгрузок..................57

2.9 Измерения возмущающих моментов................59

3 Определение орбиты 64

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

3.2 Алгоритм решения..........................66

3.3 Уточнение орбиты К А «Спектр-Р».................70

3.3.1 Модели движения ......................71

3.3.2 Траекторные измерения...................72

3.3.3 Телеметрическая информация ...............75

3.3.4 Результаты уточнения орбиты ...............76

3.3.5 Проверка лазерными измерениями.............91

4 Прогнозирование параметров движения 93

4.1 Моделирование будущих возмущений...............94

4.1.1 Расчет времен разгрузок...................95

4.1.2 Расчет величин импульсов разгрузок ...........96

4.2 Результаты прогнозирования движения К А «Спектр-Р».....99

4.3 Прогноз видимого блеска......................107

4.3.1 Модель блеска ........................108

4.3.2 Моделирование видимого блеска КА «Спектр-Р».....111

Заключение 114

Список рисунков 116

Список таблиц 117

Литература 118

Введение

Одним из приоритетных пунктов Стратегии развития космической деятельности России до 2030 года являются фундаментальные космические исследования. К способам реализации подобных исследований в текущем десятилетии в первую очередь относится развертывание внеатмосферных астрофизических обсерваторий, а также запуск орбитальных зондов и посадочных аппаратов для исследование Луны. Большинство космических аппаратов (КА), призванных решать описанные исследовательские задачи, разрабатывается в НПО им. С. А. Лавочкина, к ним относится запущенная астрофизическая обсерватория «Спектр-Р» (проект Радиоастрон[30; 36]), а также предусмотренные федеральной космической программой миссии «Спектр-РГ», «Луна-Глоб», «Спектр-УФ», «Гамма-400» и «Спектр-М». Перечисленные КА имеют различную научную нагрузку и требования к баллистико-навигационному обеспечению. Обсерватории на базе аппаратов «Спектр-Р» и «Спектр-М» в числе прочего могут работать в режиме наземно-космического интерферометра. Для успешной корреляции научных данных, полученных в таком режиме, требуется высокая точность знания движения этих аппаратов. Доставка КА «Луна-Глоб» на рабочую орбиту спутника Луны требует уточнения орбиты перелета и расчета необходимых параметров маневра в сжатые сроки. КА «Спектр-РГ» должен быть доставлен в окрестность точки либрации 1/2 системы Солнце-Земля и удерживаться на квазиустойчивой орбите при помощи реактивных двигателей. От точности баллистико-навигационного обеспечения (БНО) аппарата зависит величина корректирующих импульсов и, как следствие, время его активного существования.

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

ентации, основанные на работе маховнчных электромеханических исполнительных органов (ЭМИО), управление которыми позволяет компенсировать внешние возмущающие моменты и менять ориентацию аппарата в пространстве. С течением времени маховики ЭМИО достигают ограничений по скорости вращения, вследствие чего возникает необходимость разгрузки системы — значительного сокращения угловой скорости маховиков и компенсации суммарного кинетического момента КА при помощи реактивных двигателей стабилизации (ДС). Поскольку двигатели являются частью базового модуля, то в общем случае они не позволяют реализовать моментную схему разгрузки, поэтому каждая разгрузка маховиков создает небольшое возмущение в движении центра масс. Для запущенного в 2011 году КА «Спектр-Р» средняя величина приращения скорости в результате разгрузки составляет 5 мм/с. Игнорирование подобных возмущений при определении движения КА даже на коротких интервалах приводит к ошибкам, выходящим за рамки технического задания.

Влияние светового давления на движение КА зависит от его ориентации и свойств поверхности. Эта зависимость должна быть отражена в динамической модели аппарата. Точность модели светового давления особенно важна для космических обсерваторий «Спектр-Р» и «Спектр-М», оборудованных десятиметровыми параболическими антеннами для наблюдений в радиодиапазоне. Отношение эффективной площади сечения к массе этих аппаратов достигает 0.03 — 0.07 м2/кг, а соответствующее возмущение от светового давления всего на 1-2 порядка меньше гравитационного возмущения. Помимо этого, в зависимости от ориентации относительно Солнца световое давление создает момент относительно центра масс, действие которого компенсируется работой ЭМИО. Таким образом, световое давление оказывает влияние на возмущения от разгрузок маховиков.

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

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

Тема данной работы, посвященной проблеме навигации космических аппаратов, определению и прогнозированию параметров их движения с учетом как внешних возмущений, так и возмущений, обусловленных работой бортовых систем аппарата, является весьма актуальной. Актуальность темы подтверждена двухлетним опытом работ по баллистико-навигационному обеспечению полета и навигационной привязки научных экспериментов космического аппарата «Спектр-Р», запущенного в июле 2011 г. в соответствии с проектом создания высокоорбитального наземно-космического радиоинтерферометра «Радиоастрон». Большое удаление К А от Земли обеспечивает более высокое разрешения по сравнению с наземными интерферометрами и космическими интерферометрами на низких орбитах. Реализовать это преимущество можно только при условии высокоточного определения параметров базы интерферометра, а следовательно, параметров движения КА. Основные требования к точности определения параметров движения КА составляют по положению ± 600 м, по скорости ± 2 см/с, по ускорению ±10~8 м/с2. Выполнить эти требования невозможно без учета возмущений, обусловленных работой бортовых систем аппарата.

В рамках поставленной в работе цели была решена задача разработки параметризованной модели светового давления К А, учитывающей его форму и характеристики поверхности. Данный подход имеет ряд преимуществ по сравнению с эмпирическим и аналитическим, основанным на точных априорных оценках формы и поверхности КА, поскольку с одной стороны позволяет избавиться от ошибок априорного определения параметров и эффектов, связанных со старением материалов поверхности, с другой — будучи основанным на физическом взаимодействии потока излучения и поверхности КА, позволяет описать динамику не менее точно, чем эмпирические модели, зачастую лишенные этого свойства. Эффективность подхода была продемонстрирована в работах [2; 11; 21]. Однако, описанные в работах модели используют простую форму поверхности КА, что делает их неприменимыми для аппаратов сложной формы или со сложной структурой возникающей тени.

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

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

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

Научную новизну работы составляют: параметризованная модель силы и момента светового давления, метод использования измеренных скоростей вращения маховиков и разгрузок ЭМИО при уточнении параметров движения КА, метод долгосрочного прогнозирования движения КА с учетом прогноза изменения кинетического момента и разгрузок ЭМИО, а также математическая модель видимого блеска КА «Спектр-Р», согласующаяся с моделью влияния светового давления на движение аппарата.

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

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

• Метод использования скоростей вращения маховиков ЭМИО в качестве измерений внешнего возмущающего момента при уточнении параметров движения К А. Расчетная величина возмущающего момента зависит от уточняемых параметров модели светового давления.

• Метод долгосрочного прогнозирования движения КА с учетом разгрузок ЭМИО. Прогноз изменения кинетического момента, накапливаемого ДМ, рассчитывается при заданных параметрах светового давления и программе ориентации аппарата. Расчет возмущений от будущих разгрузок строится на основании прогноза накопления ДМ кинетического момента.

• Математическая модель видимого блеска К А «Спектр-Р», которая согласуется с моделью влияния светового давления на движение аппарата.

Модели и методики, разработанные в диссертации, использованы в рамках баллистико-навигационного обеспечения К А «Спектр-Р» — основного элемента международного проекта «Радиоастрон». Результаты, полученные в работе, позволили решить ряд важных прикладных задач, основной из которых является апостериорное определение параметров (реконструкция) движения К А «Спектр-Р», обеспечивающее высокоточную навигационную привязку научных измерений. Точность полученных реконструированных орбит позволила осуществить корреляцию интерферометрических измерений наземно-космического интерферометра. Долгосрочное прогнозирование движения КА «Спектр-Р» с использованием разработанного метода позволило избежать больших продольных ошибок в расчетном движении КА и соответствующих им временных ошибок при планировании работы наземно-космического интерферометра. Расчетное значение видимого блеска аппарата используется для планирования работы оптических наблюдательных средств. Использование этих данных позволило расширить набор телескопов, способных наблюдать КА, а также предсказывать моменты времени, неблагоприятные для проведения наблюдений.

Результаты, представленные в диссертации, докладывались автором на симпозиуме Radioastron International Science Council (RISC) (Пущино, июнь 2012 г.), на 5-й Международной конференции «Наблюдение околоземных космических объектов» (Москва, ноябрь 2011 г.) и на съезде RISC (Москва, август 2013 г.).

Основные результаты по теме диссертации изложены в 5 статьях, опубликованных печатных изданиях из перечня ВАК [32; 33; 35; 37; 38], 4 из них в соавторстве. Из работ, выполненных в соавторстве, в диссертацию включены только результаты, полученные автором.

Диссертация состоит из введения, четырех глав и заключения, содержащего основные результаты, полученные в ходе исследования. Полный объем диссертации составляет 122 страницы с 46 рисунками и 13 таблицами. Список литературы содержит 44 наименования.

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

Построение модели светового давления производится на примере КА «Спектр-Р». Для аппарата вводится упрощенная форма поверхности, согласно которой рассчитывается, какие участки поверхности освещены Солнцем в зависимости от ориентации. Для элементов поверхности К А вводятся параметры, при помощи которых производится расчет силы и момента светового давления.

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

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

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

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

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