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

российская академия наук

Ордена Ленина ИНСТИТУТ ПРИКЛАДНОЙ МАТЕМАТИКИ им. М.В. Келдыша

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

БАБКИН ЕВГЕНИЙ ВЯЧЕСЛАВОВИЧ

ОПРЕДЕЛЕНИЕ ВРАЩАТЕЛЬНОГО ДВИЖЕНИЯ ОРБИТАЛЬНЫХ СТАНЦИЙ И АНАЛИЗ МИКРОГРАВИТАЦИОННОЙ ОБСТАНОВКИ ПРИ ПРОВЕДЕНИИ КОСМИЧЕСКИХ ЭКСПЕРИМЕНТОВ

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

АВТОРЕФЕРАТ

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

чех 11

Москва-2004

Работа выполнена в Ракетно-космической корпорации "Энергия" им. СП. Королева.

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

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

доктор физико-математических наук, профессор В.В. Сазонов доктор технических наук, профессор М.Ю. Беляев

доктор физико-математических наук, профессор М.Ю. Овчинников

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

кандидат технических наук СВ. Симаков

Центральный научно-исследовательский институт машиностроения

Защита диссертации состоится " " 2004 г. на засе-

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

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

Автореферат разослан" " 2004 г.

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

доктор физико-математических наук Т.А. Полилова

Диссертационного совета

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

Актуальность проблемы. На начальном этапе развертывания Международной космической станции (МКС) и в последние полтора года полета орбитального комплекса (ОК) "Мир" длительное хранение этих объектов на орбите без экипажа осуществлялось в режиме неуправляемого движения относительно центра масс. Тип движения выбирался из условия достаточно большого среднесуточного энергосъема с солнечных батарей объекта при достаточно малой абсолютной угловой скорости последнего. Мониторинг (систематическое наблюдение) вращательного движения сводился к регулярной реконструкции этого движения на сравнительно коротких интервалах времени по результатам обработки данных измерений бортовых датчиков ориентации. Реконструкция управляемого и неуправляемого вращательного движения космического аппарата необходима также для расчета квазистатической компоненты остаточных микроускорений на его борту. Мониторинг остаточных микроускорений проводился на ОК "Мир" и систематически проводится на МКС для планирования выполнения космических экспериментов с гравитационно-чувствительными системами.

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

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

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

Апробация работы. Результаты диссертации докладывались на VIII Всероссийском съезде по теоретической и прикладной механике (г. Пермь, 23-29 августа 2001 г.), XXXVI Чтениях К.Э. Циолковского (г. Калуга, 18-20 сентября 2001 г.), семинаре «Механика невесомости и гравитационно -чувствительные системы» Института проблем механики РАН (рук. В.И. Полежаев и В.В. Сазонов), семинаре по механике космического полета механике- национальная

I БИБЛИОТЕКА

I С.П<тер&рг и^

• ,_ О»

ко-математического факультета МГУ (рук. В.В. Белецкий, К.Г. Григорьев и В.В. Сазонов) и научно-технических конференциях РКК Энергия.

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

Диссертация состоит из введения, трех глав и заключения.

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

Первая глава диссертации содержит описание методов реконструкции неуправляемого вращательного движения ОК "Мир" и МКС, а также анализ наиболее интересных результатов, полученных в процессе мониторинга движения этих объектов в 1999 - 2001 гт. Движение МКС определялось по измерениям трехосного датчика угловой скорости, движение ОК "Мир" - по измерениям трехосного магнитометра. Для описания неуправляемого вращательного движения обеих орбитальных станций использовалась единая математическая модель, построенная при следующих допущениях. Станция считалась твердым телом, геоцентрическое движение центра масс которого было кеплеровым эллиптическим. Элементы этого движения определялись по данным радиоконтроля орбиты. В уравнениях движения станции относительно центра масс учитывались гравитационный и восстанавливающий аэродинамический моменты. Последний рассчитывался в предположении, что станция имеет форму сферы, центр которой смещен относительно ее центра масс. Координаты центра сферы и ее радиус (точнее, три комбинации этих четырех величин и минимального главного центрального момента инерции станции) служили параметрами согласования - определялись в процессе обработки данных измерений наряду с начальными условиями движения. Плотность набегающего на станцию аэродинамического потока рассчитывалась согласно модели верхней атмосферы, определенной ГОСТ 22721-77. Учитывался перекос между главными центральными осями инерции станции, в которых записывались уравнения ее вращательного движения, и осями строительной системы координат, в которой интерпретировались измерения датчиков ориентации.

Данные измерений, выполненные на временном интервале длиной несколько часов, обрабатывались совместно методом наименьших квадратов с помощью интегрирования уравнений движения станции относительно центра масс. Число данных, относящихся к обрабатываемому интервалу, обычно оказывалось весьма большим. В такой ситуации не было необходимости использовать все данные непосредственно при построении сглаживающего их решения уравнений движения. Имеющаяся информация предварительно сжималась. При этом выполнялась оценка ее качества, и заведомо неточные измерения отбрасывались. Сжатие выполнялось с помощью построения функций, аппроксимирующих данные измерений. Эти функции строились в виде дискретных рядов Фурье. На решениях уравнений движения станции задавалсяЛ функционал, который представлял собой сумму квадратов разно-

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

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

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

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

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

Недостаток цилиндрической прецессии заключается в ее неустойчивости при достаточно малых абсолютных значениях угловой скорости вращения тела вокруг оси симметрии. Граничные по устойчивости значения угловой скорости тем выше по модулю, чем более вытянуто тело вдоль оси симметрии. Указанное обстоятельство делает движения, близкие к цилиндрической прецессии, практически непригодными для вытянутых спутников из-за необходимости закрутки с большой угловой скоростью. Для ОК «Мир» указанные граничные значения угловой скорости оказались достаточно малыми - около 0.14 градУс. При реализации этих движений начальные значения угловой скорости ОК были примерно равны 0.2 градУс.

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

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

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

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

Перейдем к МКС. На этой станции в 1999 г. применялся режим гравитационной ориентации вращающегося спутника. В этом режиме станция вращалась с угловой скоростью 0.2 - 0.4 град./с вокруг продольной оси, направленной приблизительно вдоль местной вертикали. Для применения такого режима необходимо выполнение трех условий. Во-первых, спутник должен иметь специфический центральный эллипсоид инерции: большая и средняя полуоси этого эллипсоида должны мало отличаться друг от друга и быть существенно больше (в три и более раз) малой полуоси. Во-вторых, приложенный к спутнику гравитационный момент должен существенно превышать другие действующие на спутник механические моменты. В-третьих, орбита спутника должна быть близка к круговой. Если перечисленные условия выполнены, то существуют движения спутника, близкие к так называемой конической прецессии осесимметричного твердого тела на круговой орбите под действием гравитационного момента. При этом в случае малой угловой скорости закрутки спутника вокруг продольной оси отклонения последней от местной вертикали будут малыми. Такой режим применялся на орбитальных комплексах «Салют-6» - «Космос-1267», «Салют-7» - «Кос-мос-1443» и «Салют-7» - «Космос-1686». МКС в конфигурации 1999 г. также удовлетворяла всем перечисленным условиям.

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

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

Во второй главе диссертации описывается методика расчета реальной квазистатической составляющей микроускорения по информации о враща-

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

Микроускорением точки, закрепленной на корпусе космического аппарата, называется разность между напряженностью гравитационного поля в этой точке и абсолютным ускорением последней. Возникновение остаточных микроускорений на борту искусственного спутника Земли в неуправляемом полете обусловлено несколькими причинами. Главные из них: 1) движение спутника относительно центра масс, 2) градиент гравитационного поля, 3) аэродинамическое торможение, 4) упругие колебания конструкции. Микроускорения измеряются акселерометрами, их показания на крупногабаритных спутниках типа ОК «Мир» и МКС естественным образом разбиваются на две составляющие - высокочастотную и квазистатическую (низкочастотную). Квазистатическая составляющая имеет спектр в диапазоне от нуля до нескольких тысячных долей герца и обусловлена первыми тремя указанными выше причинами, которые могут реализоваться для спутника - твердого тела. Высокочастотная составляющая имеет спектр в диапазоне выше нескольких сотых долей герца и вызвана упругими колебаниями конструкции спутника и функционированием его бортовых устройств.

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

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

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

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

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

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

В третьей главе описывается методика обработки данных измерений высокочастотных акселерометров, и приводятся результаты определения вибрационной компоненты микроускорения на ОК «Мир» и МКС.

Во время выполнения экспериментов с установкой ДАКОН вибрационная компонента микроускорения на ОК "Мир" измерялась с помощью де-вятиканальной системы акселерометров ВМ-09 (Виброметр Многоканальный, 9 - максимальное число каналов), разработанной в НПО "Композит". Амплитудный диапазон ускорений, измерявшихся этой системой, - 0.0001 -0.01м/с2, частотный диапазон - от нескольких десятых долей герца до 100 Гц. Скорость выборки прибора - 200 отсчетов в секунду, длина непрерывного сеанса измерений — не более 10 мин. Снимаемые данные передавались на Землю в режиме непосредственной передачи во время пролета ОК над приемным пунктом телеметрической информации. При обработке этих данных на Земле решались следующие задачи: 1) выделение циклических трендов из данных измерений, 2) оценка спектральной плотности составляющей этих данных с непрерывным спектром, 3) низкочастотная фильтрация данных измерений. Полученные данные использовались в задачах математического моделирования экспериментов, выполненных на борту ОК.

На российском сегменте МКС вибрационная составляющая микроускорения измеряется с помощью 10 трехкомпонентных акселерометров ИМУ-128. Пять из них расположены на Функциональном грузовом блоке (ФГБ) и пять - на СМ. Для увеличения точности измерений в области малых амплитуд измерения по каждой компоненте передаются по двум телеметрическим каналам. По одному каналу передается полный диапазон амплитуд до 0.01 м/с2, по другому - в 10 раз меньший — до 0.001 м/с2. Полученные данные измерений представляют собой отрезки временных рядов длительностью в несколько минут. Обработка этих данных включала решение двух задач: 1) выделение циклических трендов, 2) оценивание спектральной плотности ряда остатков, получающегося из исходного ряда удалением значимых цикли-

ческих трендов. Как оказалось, вибрационная составляющая микроускорения на ФГБ не превышала 0.01 м/с2 в диапазоне частот от 0 до 20 Гц. Во время спокойной обстановки на СМ вибрационная составляющая микроускорения также не превосходила указанного предела. При включенном компрессоре системы кондиционирования воздуха уровень микроускорений на СМ повышался более чем на порядок.

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

ПОЛОЖЕНИЯ, ВЫНОСИМЫЕ НА ЗАЩИТУ.

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

2. Разработана и реализована в виде программ на персональном компьютере методика расчета реальных квазистатических микроускорений на ОК «Мир» и МКС по телеметрической информации о вращательном движении. Методика предназначена для контроля микроускорений во время проведения космических экспериментов с гравитационно-чувствительными системами. Она использовалась при анализе результатов экспериментов с датчиком конвекции ДАКОН и установкой ALICE, выполненных на ОК «Мир» в 1998 — 1999 гг. Точность методики подтверждена совпадением результатов расчетов с данными измерений акселерометра MAMS на МКС л

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

СПИСОК РАБОТ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ

1. Е.В. Бабкин, М.Ю. Беляев, Н.И. Ефимов, В.В. Сазонов, В.М. Стажков.

Неуправляемое вращательное движение станции «Мир». Космические

исследования, 2001, т. 39, № 1, с. 27-42.

2. В.В. Сазонов, М.Ю. Беляев, Н.И. Ефимов, В.М. Стажков, Е.В. Бабкин.

Определение квазистатической компоненты микроускорения на станции «Мир». Космические исследования, 2001, т. 39, № 2, с. 136-147. '

0-4884

3. В.И. Ветлов, СМ. Новичкова, В.В. Сазонов, Н В. Матвеев, Е.В. Бабкин. Режим гравитационной ориентации Международной космической станции. Космические исследования, 2001, т. 39, №4, с. 436-448.

4. Е.В. Бабкин, М.Ю. Беляев, НИ. Ефимов, В.В. Сазонов, В.М. Стажков. Неуправляемое вращательное движение орбитальной станции «Мир» в последние месяцы ее полета. Космические исследования, 2003, т. 41, № 3, с. 285-294.

5. Е.В. Бабкин, М.Ю. Беляев, Н.И. Ефимов, В.В. Сазонов, В.М. Стажков Неуправляемое вращательное движение орбитальной станции «Мир» в 1999 г. Препринт ИПМ им. М.В. Келдыша РАН, 2000, № 15.

6. В.И. Ветлов, СМ. Новичкова, ВВ. Сазонов, Н.В. Матвеев, Е.В. Бабкин. Режим гравитационной ориентации Международной космической станции. Препринт ИПМ им. М.В. Келдыша РАН, 2000, № 24.

7. Е.В. Бабкин, М.Ю. Беляев, Н И. Ефимов, С.С. Обыденников, В.В. Сазонов, В.М. Стажков. Первые результаты определения микроускорений на российском сегменте международной космической станции: Препринт ИПМ им. М.В. Келдыша РАН, 2001, № 83.

8. Е.В. Бабкин, М.Ю. Беляев, Н И. Ефимов, В.В. Сазонов, В.М. Стажков. Определение квазистатической компоненты микроускорения, возникающего на борту Международной космической станции. Препринт ИПМ им. М.В. Келдыша РАН, 2003, №11.

9. Е.В. Бабкин, М.Ю. Беляев, Н.И. Ефимов, В.М. Стажков, В.В. Сазонов. Неуправляемое вращательное движение орбитальной станции "Мир". Материалы VIII Всероссийского съезда по теоретической и прикладной механике, Пермь, 23-29 августа 2001 г.

10.Е.В. Бабкин, М.Ю. Беляев, Н.И. Ефимов, С.С. Обыденников, ВВ. Сазонов, В.М. Стажков, В.В. Измерение и моделирование микроускоре-, ний на МКС. Материалы XXXVI Чтений, посвященных разработке научного наследия и развития идей К.Э. Циолковского, Калуга, 18-20 сентября 2001 г.

И П М. 3 а I а з № 14. Тираж! 00 экс.

Введение.

1. Мониторинг неуправляемого вращательного движения орбитальных станций.

1.1. Задачи, решаемые мониторингом неуправляемого движения.

1.2. Уравнения вращательного движения станции.

1.3. Определение движения по измерениям бортовых магнитометров.

1.4. Определение движения по телеметрической информации.

1.5. Учет априорной информации о параметрах математической модели.

1.6. Закрутка станции «Мир» вокруг нормали к плоскости орбиты.

1.7. Вращение станции «Мир» вокруг оси максимального момента инерции, совершающей колебания относительно нормали к плоскости орбиты.

1.8. Движение станции «Мир», близкое регулярной прецессии Эйлера.

1.9. Двухосная закрутка станции «Мир» в плоскости орбиты.

1.10. Режим гравитационной ориентации МКС.

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

2.1. Использование информации о вращательном движении станции для определения микроускорений на ее борту.

2.2. Аппроксимация вращательного движения станции с помощью сплайнов.

2.3. Низкочастотная составляющая в телеметрических значениях угловой скорости.

2.4. Аппроксимация вращательного движения станции с решением вариационной задачи.

2.5. Микроускорения в неуправляемом движении станции.

2.6. Определение квазистатической составляющей микроускорения на борту МКС по телеметрическим измерениям.

2.7. Сглаживание данных измерений низкочастотного акселерометра MAMS.

2.8. Сравнение данных измерения с результатами вычисления.

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

3.1. Измерение микроускорений на орбитальной станции «Мир» с помощью системы акселерометров ВМ-09.

3.2. Интерпретация данных измерений системы ВМ-09.

3.3. Спектральный анализ данных измерений бортовых акселерометров.

3.4. Спектральный анализ данных измерений системы ВМ-09.

3.5. Фильтрация данных измерений системы ВМ-09.

3.6. Измерение вибрационной составляющей микроускорения на МКС акселерометрами ИМУ-128.

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

Реконструкция управляемого и неуправляемого вращательного движения орбитального комплекса необходима также для расчета квазистатической компоненты остаточных микроускорений на его борту. Знание микроускорений необходимо при анализе космических экспериментов с гравитационно-чувствительными системами. Мониторинг остаточных микроускорений проводился на ОК "Мир" и систематически проводится на МКС для планирования космических экспериментов с такими системами.

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

Диссертация содержит три главы. В первой главе описаны методы реконструкции неуправляемого вращательного движения ОК "Мир" и МКС, а также представлены наиболее интересные результаты, полученные в процессе мониторинга движения этих объектов в 1999 - 2001 гг. Движение МКС определялось по измерениям трехосного датчика угловой скорости, движение ОК "Мир" - по измерениям трехосного магнитометра. Для описания неуправляемого вращательного движения обеих орбитальных станций использовалась единая математическая модель, построенная при следующих допущениях. Станция считалась твердым телом, геоцентрическое движение центра масс которого было кеплеровым эллиптическим. Элементы этого движения определялись по данным радиоконтроля орбиты. В уравнениях движения станции относительно центра масс учитывались гравитационный и восстанавливающий аэродинамический моменты. Последний рассчитывался в предположении, что станция имеет форму сферы, центр которой смещен относительно ее центра масс. Координаты центра сферы и ее радиус (точнее, три комбинации этих четырех величин и минимального центрального момента инерции станции) служили параметрами согласования - определялись в процессе обработки данных измерений наряду с начальными условиями движения. Плотность набегающего на станцию аэродинамического потока рассчитывалась в рамках допущений, определяемых ГОСТ 22721-77. Учитывался перекос между главными центральными осями инерции станции, в которых записывались уравнения ее вращательного движения, и осями строительной системы координат, в которой интерпретировались измерения датчиков ориентации.

Данные измерений, выполненные на временном интервале в несколько часов, обрабатывались совместно методом наименьших квадратов с помощью интегрирования уравнений движения станции относительно центра масс. Число данных, относящихся к обрабатываемому интервалу, обычно оказывалось весьма большим. В такой ситуации не было необходимости использовать все данные непосредственно при построении сглаживающего их решения уравнений движения. Имеющаяся информация предварительно сжималась. При этом выполнялась оценка ее качества, и заведомо неточные измерения отбрасывались. Сжатие выполнялось с помощью построения функций, аппроксимирующих данные измерений. Эти функции строились в виде дискретных рядов Фурье. На решениях уравнений движения станции задавался функционал, который представлял собой сумму квадратов разностей между значениями функций, аппроксимирующих данные измерений, и расчетными аналогами этих данных, вычисляемыми вдоль решения уравнений движения. Указанная сумма вычислялась по моментам времени, образующим равномерную сетку с шагом 1-2 мин внутри обрабатываемого интервала. Минималь такого функционала считалась реконструкцией фактического движения станции на обрабатываемом интервале. Минимизация функционала выполнялась методом Гаусса-Ньютона по начальным условиям движения станции и параметрам математической модели. Использовались два варианта минимизации. В одном из них центральный тензор инерции станции, отнесенный к осям строительной системы координат, был фиксирован, в другом пять из шести характеристик этого тензора уточнялись. А именно, уточнялись три угла, задающих положение главных центральных осей инерции станции в строительной системе координат, и две безразмерных комбинации трех главных центральных моментов инерции станции.

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

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

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

Тензор инерции ОК "Мир" в конфигурации 1999 г. был близок к осесим-метричному. Два главных центральных момента инерции станции почти совпадали между собой и заметно превосходили третий момент. Отношение минимального и максимального из этих моментов составляло примерно 0.62, отношение среднего и максимального - 0.91. Указанное обстоятельство позволило использовать для беспилотного полета ОК неуправляемое движение относительно центра масс, близкое так называемой цилиндрической прецессии осе-симметричного твердого тела на круговой орбите. Это движение представляло собой вращение ОК вокруг своей продольной оси (оси минимального момента инерции), совершавшей малые колебания относительно нормали к плоскости орбиты.

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

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

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

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

Перейдем к МКС. На этой станции в 1999 г. применялся режим гравитационной ориентации вращающегося спутника. В этом режиме станция вращалась вокруг продольной оси, направленной приблизительно вдоль местной вертикали. Для применения такого режима необходимо выполнение трех условий. Во-первых, спутник должен иметь специфический центральный эллипсоид инерции: большая и средняя полуоси этого эллипсоида должны мало отличаться друг от друга и быть существенно больше (в три и более раз) малой полуоси. Во-вторых, приложенный к спутнику гравитационный момент должен существенно превышать другие действующие на спутник механические моменты. В-третьих, орбита спутника должна быть близка к круговой. Если перечисленные условия выполнены, то существуют движения спутника, близкие к так называемой конической прецессии осесимметричного твердого тела на круговой орбите под действием гравитационного момента. При этом в случае малой угловой скорости закрутки спутника вокруг продольной оси отклонения последней от местной вертикали будут малыми. Такой режим применялся на орбитальных комплексах «Салют-6» - «Космос-1267», «Салют-7» - «Космос-1443» и «Салют-7» - «Космос-1686». МКС в конфигурации 1999 г. также удовлетворяла всем перечисленным условиям.

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

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

Во второй главе описывается методика расчета реальной квазистатической составляющей микроускорения по информации о вращательном движении орбитальной станции, и приводятся результаты таких расчетов для ОК «Мир» и МКС. Напомним, что микроускорением точки, закрепленной на корпусе космического аппарата, называется разность между напряженностью гравитационного поля в этой точке и абсолютным ускорением последней. Возникновение остаточных микроускорений на борту искусственного спутника Земли в неуправляемом полете обусловлено несколькими причинами. Главные из них: 1) движение спутника относительно центра масс, 2) градиент гравитационного поля, 3) аэродинамическое торможение, 4) упругие колебания конструкции. Микроускорения измеряются акселерометрами, их показания на крупногабаритных спутниках типа ОК «Мир» и МКС естественным образом разбиваются на две составляющие - высокочастотную и квазистатическую (низкочастотную). Квазистатическая составляющая имеет спектр в диапазоне от нуля до нескольких тысячных долей герца и обусловлена первыми тремя указанными выше причинами, которые могут реализоваться для спутника - твердого тела. Высокочастотная составляющая имеет спектр в диапазоне выше нескольких сотых долей герца и вызвана упругими колебаниями конструкции спутника и функционированием его бортовых устройств.

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

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

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

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

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

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

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

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

В третьей главе рассматриваются данные измерений двух акселерометров. Один из них, ВМ-09 (Виброметр Многоканальный, 9 - максимальное число каналов), был разработан в НПО "Композит" и установлен на ОК «Мир». Фактически это была система из шести однокомпонентных акселерометров. Она использовалась для измерения вибрационной составляющей микроускорения во время выполнения экспериментов с установкой ДАКОН. Амплитудный диапазон ускорений, измерявшихся этой системой, - 0.0001 - 0.01м/с2, частотный диапазон - от нескольких десятых долей герца до 100 Гц. Скорость выборки прибора - 200 отсчетов в секунду, длина непрерывного сеанса измерений -не более 10 мин. Снимаемые данные передавались на Землю в режиме непосредственной передачи во время пролета ОК над приемным пунктом телеметрической информации. При обработке этих данных на Земле решались следующие задачи: 1) выделение циклических трендов из данных измерений, 2) оценка спектральной плотности составляющей этих данных с непрерывным спектром, 3) низкочастотная фильтрация данных измерений. Полученные данные использовались в задачах математического моделирования экспериментов, выполненных на борту ОК.

Второй прибор, данные измерений которого анализируются в третьей главе, это трехкомпонентный акселерометр ИМУ-128. На российском сегменте МКС установлено 10 таких приборов. Пять из них расположены на Функциональном грузовом блоке (ФГБ) и пять - на СМ. Для увеличения точности измерений в области малых амплитуд измерения по каждой компоненте передаются по двум телеметрическим каналам. По одному каналу передается полный диапазон амплитуд до 0.01 м/с2, по другому - в 10 раз меньший - до 0.001 м/с2. Полученные данные измерений представляют собой отрезки временных рядов длительностью в несколько минут. Обработка этих данных включала решение двух задач: 1) выделение циклических трендов, 2) оценивание спектральной плотности ряда остатков, получающегося из исходного ряда удалением значимых циклических трендов. Анализ данных измерений акселерометров ИМУ-128 позволил получить первые оценки уровня вибрационных микроускорений на Российском сегменте МКС. Вибрационная составляющая микроускорения на ФГБ в диапазоне частот от 0 до 20 Гц не превышала 0.01 м/с2. Во время спокойной обстановки на СМ вибрационная составляющая микроускорения также не превосходила указанного предела. При включенном компрессоре системы кондиционирования воздуха уровень микроускорений на СМ повышался более чем на порядок.

Основные результаты диссертации, выносимые на защиту, сформулированы в Заключении.

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

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

2. Разработана и реализована в виде программ на персональном компьютере методика расчета реальных квазистатических микроускорений на ОК «Мир» и МКС по телеметрической информации о вращательном движении. Методика предназначена для контроля микроускорений во время проведения космических экспериментов с гравитационно-чувствительными системами. Она использовалась при анализе результатов экспериментов с датчиком конвекции ДАКОН и установкой ALICE, выполненных на ОК «Мир» в 1998 - 1999 гг. Точность методики подтверждена совпадением результатов расчетов с данными измерений акселерометра MAMS на МКС.

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

Заключение.

1. Сарычев В.А., Сазонов В.В. Гравитационная ориентация вращающегося спутника. Космические исследования, 1981, т. 19, вып. 4, с. 499-512.

2. Сазонов В.В., Петров АЛ. Эволюция режима гравитационной ориентации вращающегося спутника под действием непотенциального аэродинамического момента. Космические исследования, 1987, т. 25, вып. 4, с. 508-522.

3. Ветлов В.И., Сазонов В.В., Сарычев В.А. Влияние демпфирования на режим гравитационной ориентации вращающегося спутника. Известия АН СССР. Механика твердого тела, 1990, вып. 1, с. 3-11.

4. Костенко И.К., Ветлов В.И., Нырков А.Г., Сарычев В.А., Сазонов В.В. Режим обобщенной гравитационной ориентации на орбитальных комплексах "Салют-6" "Космос-1267" и "Салют-7" - "Космос-1443". Космические исследования, 1986, т. 24, вып. 1, с. 46-51.

5. Ветлов В.И., Новичкова С.М., Сазонов В.В. Исследование режима гравитационной ориентации вращающегося спутника. Препринт Института прикладной математики РАН, 1995, N 24.

6. Сарычев В.А., Сазонов В.В., Беляев М.Ю., Ефимов Н.И., Лапшина И. Л. Определение пассивного вращательного движения орбитальной станции "Мир" по измерениям напряженности геомагнитного поля. Космические исследования, 1995, т. 33, N1, с. 12-19.

7. Модель верхней атмосферы для баллистических расчетов. ГОСТ 2272177. М., Изд-во стандартов, 1978.

8. Сарычев В.А., Сазонов В.В., Беляев М.Ю., Ефимов Н.И. Повышение точности определения вращательного движения орбитальных станций "Салют-6" и "Салют-7" по данным измерений. Космические исследования, 1991, т. 29, N 3, с. 375-389.

9. Черноусько Ф.Л. Об устойчивости регулярной прецессии спутника. Прикладная математика и механика, 1963, т. 28, вып. 1, с. 155-157.

10. Сарычев В.А., Беляев М.Ю., Сазонов В.В., Тян Т.Н. Определение микроускорений на орбитальных комплексах "Салют-6" и "Салют-7". Космические исследований, 1986, т. 24, N 3, с. 337-344.

11. Сазонов В.В., Комаров М.М., Беляев М.Ю., Зыков С.Г., Стажков В.М. Оценка квазистатической компоненты микроускорения на борту искусственного спутника Земли. Препринт ИПМ им. М.В.Келдыша РАН, 1995, N45.

12. Сазонов В.В., Зыков С.Г. Расчет и измерение низкочастотной составляющей микроускорения на борту ИСЗ. Препринт ИПМ им. М.В.Келдыша РАН, 1996, N 31.

13. Сазонов В.В., Чебуков С.Ю., Абрашкин В.И., Казакова А.Е., Зайцев А.С. Анализ низкочастотных микроускорений на борту ИСЗ "Фотон-11". Препринт ИПМ им. М.В.Келдыша РАН, 1999, N 33.

14. Белецкий В.В. Движение искусственного спутника относительно центра масс. М., Наука, 1965.

15. Беляев М.Ю., Ефимов Н.И., Сазонов В.В. Определение ориентации орбитального комплекса "Мир" по показаниям оптического звездного датчика. Космические исследования, 1995, т. 33, N 4, с. 395-402.

16. Бабкин Е.В., Беляев М.Ю., Ефимов Н.И., Сазонов В.В., Стажков В.М. Неуправляемое вращательное движение орбитальной станции "Мир" в 1999 г. Препринт ИПМ им. М.В.Келдыша РАН, 2000, N 15.

17. Беляев М.Ю., Зыков С.Г., Рябуха С.Б., Сазонов В.В., Сарычев В.А., Стажков В.М. Математическое моделирование и измерение микроускорений на орбитальной станции "Мир"// Известия Академии наук. Механика жидкости и газа. 1994. N5. С. 5-14.

18. Сазонов В.В., Комаров М.М., Беляев М.Ю., Ефимов Н.И., Зыков С.Г. Оценка квазистатической компоненты микроускорения на борту орбитальной станции "Мир" по показаниям оптического звездного датчика. Космические исследования, 1996, т. 34, N 5, с. 491-499.

19. Полежаев В.И., Никитин С.А., Сазонов В.В. Измерение квазистатической составляющей микроускорения на борту ИСЗ с помощью датчика конвекции// Космические исследования 2001. Т39.№ 2.

20. Сазонов В.В., Беляев М.Ю., Ефимов Н.И., Стажков В.М., Бабкин Е.В. Определение квазистатической компоненты микроускорения на станции «Мир». Космические исследования, 2001, т. 39, № 2, с. 136-147.

21. Бабкин Е.В., Беляев М.Ю., Ефимов Н.И., Сазонов В.В., Стажков В.М. Неуправляемое вращательное движение станции «Мир». Космические исследования, 2001, т. 39, № 1, с. 27-42.

22. Reinsch С.Н. Smoothing by spline functions. Numerische mathematik, 1975, B. 24, No. 5, S. 383-393.

23. Хеннан Э. Анализ временных рядов. М., Наука, 1964.

24. Теребиж В.Ю. Анализ временных рядов в астрофизике. М., Наука, 1992.

25. Отнес Р., Эноксон Л. Прикладной анализ временных рядов. М., Мир, 1982.

26. Ланцош К. Практические методы прикладного анализа. М., Физматгиз, 1961.

27. Хемминг Р.В. Цифровые фильтры. М., Недра, 1987.

28. Бабкин Е.В., Беляев М.Ю., Ефимов Н.И., Обудёнников С.С., Сазонов В.В., Стажков В.М. Первые результаты определения микроускорений на российском сегменте международной космической станции: Препринт ИПМ им. М.В. Келдыша РАН, 2001, №83

29. B.JI. Левтов, В.В. Романов, А.И. Иванов, С.Б. Рябуха, С.С. Сазонов. Результаты летно-космических испытаний виброзащитной платформы ВЗП-1К на станции "Мир". Препринт ИПМ им. М.В.Келдыша РАН, 2000, N 46.

30. В.И. Ветлов, С.М. Новичкова, В.В. Сазонов, Н.В. Матвеев, Е.В. Бабкин. Режим гравитационной ориентации Международной космической станции. Космические исследования, 2001, т. 39, №4, с.436-448.

31. В.В. Голубков, И.Г. Хацкевич. Определение ориентации ИСЗ по заданной системе измерений. Космические исследования, 1969, т. 7, №4, с.510-521.

32. В.В. Белецкий, В.В. Голубков, Е.А. Степанова, И.Г. Хацкевич. Результаты определения ориентации спутника «Протон-2» и описание его движения относительно центра масс. Космические исследования, 1969, т. 7, №4, с.522-533.

33. Е.В. Бабкин, М.Ю. Беляев, Н.И. Ефимов, В.В. Сазонов, В.М. Стажков. Неуправляемое вращательное движение орбитальной станции «Мир» в последние месяцы ее полета. Космические исследования, 2003, т. 41, № 3, с. 285-294.

34. В.И. Ветлов, С.М. Новичкова, В.В. Сазонов, Н.В. Матвеев, Е.В. Бабкин. Режим гравитационной ориентации Международной космической станции. Препринт ИПМ им. М.В. Келдыша РАН, 2000, № 24.

35. Е.В. Бабкин, М.Ю. Беляев, Н.И. Ефимов, В.В. Сазонов, В.М. Стажков. Определение квазистатической компоненты микроускорения, возникающего на борту Международной космической станции. Препринт ИПМ им. М.В. Келдыша РАН, 2003, №11.

36. Е.В. Бабкин, М.Ю. Беляев, Н.И. Ефимов, В.М. Стажков, В.В. Сазонов. Неуправляемое вращательное движение орбитальной станции "Мир". Материалы VIII Всероссийского съезда по теоретической и прикладной механике, Пермь, 23-29 августа 2001 г.

37. Микроускорения п2, п3, |п| (1(Гб м/с2)t-t0 (мин) Рис. 25. Микроускорения на интервале 7

38. Микроускорения пх,п2,щ\г\ (ХО^м/с1)

39. Компоненты кватерниона qQ,qi,q2,q353.9t-tQ (мин)26.9539 80.8t-t0 (мин)107.8134.7а) (б)

40. Рис. 26. Момент t = 0 соответствует 19:06:23 ДМВ 17.07.1998.

41. Угловые скорости тх,со2,ео3 (10"3с 1) Угловые ускорения сь1,сй2,а3(\0~6с'2)t-t0( мин) /-/0(мин)в) (г)

42. Рис. 26. Момент t = 0 соответствует 19:06:23 ДМВ 17.07.1998.

43. Микроускорения щ,п2>пгМ (XQT6 м!с*)

44. Компоненты кватерниона <7о > Ч\ > (h > Чг47 48.4321 15.7-0.60.0839 38.6 -8.8 -96.1-103.40.024.022.122.122.144.266.344.266.344.266.300 22.1 44.2 66.3t-t0 (мин)88.41.4110.5110.5110.5110.522.1442 66.3t-t0 (мин)88.4110.5110.5110.5110.5а) (б)

45. Рис. 27 Момент t = 0 соответствует 08:22:32 ДМВ 14.10.1998.2.1210 1.6322 1.14Э4 0.6546 0.1(38

46. Угловые скорости , а>2, со3 (10"3 с-1)Л

47. Угловые ускорения fy, <02, й>3 (10-6 с~2 )-0.3230