CC BY
27
При малых углах отклонения двигателей от оси Z эффективность приложения тяги по долготе низка, а уменьшение затрат топлива ДУОС на компенсацию моментов невелико, поэтому суммарная масса увеличивается по сравнению даже с исходным вариантом (рис. 4). При увеличении угла отклонения до 30° эффективность тяги по наклонению (по оси Z) остается еще достаточной (87 % от идеальной), а эффективность тяги по долготе становится уже приемлемой (50 %), затраты топлива на компенсацию моментов заметно уменьшаются. При дальнейшем увеличении угла отклонения до 45° эффективность приложения тяги по наклонению падает до 70 %, и этот фактор перевешивает повышение эффективности тяги по долготе и уменьшение запаса топлива на компенсацию возмущающих моментов. Таким образом, при компоновке с «универсальным» назначением двигателей применительно к изделию типа «Экспресс-АМ», оптимальным является размещение двигателей в плоскости YOZ под углом примерно 30° к оси ^ При таком размещении наблюдается наилучший баланс между эффективностью тяги по наклонению и долготе, относительно невелики затраты топлива на компенсацию возмущающих моментов. За счет сокращения количества двигателей и массы СПУ достигается заметное преимущество (22-23 кг) в сравнении с базовым вариантом с 8-ю специализированными двигателями. С учетом уменьшения мас-
сы кабелей (~18 кг) и конструктивных элементов для размещения блоков коррекции на корпусе КА (~7 кг) преимущество может достигать ощутимой величины 47-48 кг.
Предложены и рассмотрены критерии эффективности размещения двигателей коррекции орбиты на геостационарных КА (оценка «по топливу», оценка «по суммарной массе»). Указанные критерии весьма наглядны и позволяют с приемлемыми трудозатратами сравнивать между собой различные варианты компоновки как для случая специализации двигателей по выполняемой задаче, так для случая универсального назначения двигателей. Получаемое с помощью оптимизации по предложенным критериям преимущество в массе, по сравнению с базовым вариантом, может составлять величину порядка 47-48 кг.
Библиографические ссылки
1. Effects of interaction of electric thruster jets and spacecraft solar battareies / M. P. Burgasov, A. B. Nadiradzer, A. A. Chirov [at al.] // European Space Power Conference ; Graz, Austria. 23-27Aug. 1993. (ESAWPP-054,Aug. 1993).
2. Разработка и поставка пакета прикладных программ для расчета воздействия струй СПД и ТКД на элементы КА : отчет о НИР. Тема№ 208-97040. Моск. гос. авиационный ин-т. М., 2001.
Yu. M. Yermoshkin, V. M. Urusov
CRITERIA OF ORBITAL CONTROL THRUSTER OPTIMALALLOCATION ON BOARD OF GEOSTATIONARY SPACECRAFT
During designing of new, especially long-life geostationary spacecrafts we face a problem of rational orbital control thruster allocation. In this paper we consider the possible criteria for solution of such problem and the results of their application.
Keywords: spacecraft, thruster, propulsion, fuel, working body.
© Ермошкин Ю. М., Урусов В. М., 2011
УДК 629.76.78.015
И. Н. Карцан, В. С. Пономарев, Е. С. Жукова, С. В. Литошик
АСТРОНОМИЧЕСКИЕ ИСТОЧНИКИ ПЕРВИЧНОЙ НАВИГАЦИОННОЙ ИНФОРМАЦИИ ДЛЯ АВТОНОМНОГО ОПРЕДЕЛЕНИЯ ДВИЖЕНИЯ ОБЪЕКТОВ*
Анализируются астрономические источники первичной навигационной информации для автономного определения движения объектов.
Ключевые слова: зенитное расстояние звезды, сжатие земного сфероида, вектор напряженности.
Астрономические источники ПНП, а также методы и ные на использовании этой информации, в значительной средства решения задачи навигации в космосе, основан- степени определяют успех выполнения всех задач, кото-
* Работа выполнена в рамках реализации ФЦП «Научные и научно-педагогические кадры России» на 2009-2013 гг. ГК № 14.740.11.0329 от 17.02.2010 г.
рые ставятся перед космическими системами самого различного назначения. В настоящее время практически нет таких космических аппаратов и кораблей, в составе бортовых систем которых не были бы представлены устройства навигации и ориентации, основывающиеся на использовании информации, получаемой от астрономических источников.
Для целей навигации и ориентации в космическом пространстве могут использоваться все поля, возбуждаемые астрономическими объектами: электромагнитное (во всем возможном диапазоне длин волн от ультрафиолета до радиодиапазона), магнитное, гравитационное. При этом первичной навигационной информацией могут являться результаты измерений любых физических величин, обусловленных полями, характеристики которых зависят от координат и скорости КА. К числу таких величин можно отнести величину и градиенты полей, интенсивность излучения, углы положения астрономических ориентиров относительно некоторых направлений, удаленность от поверхности близких небесных тел, измеряемую дальномерами [1].
Измерение положения звезды относительно направления текущего радиус-вектора.
Пусть г - радиус-вектор, определяющий положение КА в некоторой системе координат OXYZ (рис. 1), а я0 -единичный вектор, направленный из начала координат в точку расположения неподвижной звезды. Обозначим буквой в угол между г и я0. Величина этого угла определяется простым выражением
0 , 1.0 . 0 xa + yb + zc
(1)
s0 в системе
OXYZ.
Измеряя зенитное расстояние второй звезды, построим вторую коническую поверхность. Две линии пересечения (в общем случае) полученных поверхностей позволяют определить два направления возможного положения радиус-вектора r. С одним из этих направлений совпадает истинное положение r.
Уравнения прямых, соответствующих двум направлениям r, определяются в результате совместного решения двух уравнений, отличающихся только коэффициентами a0, b°, c°. Для нахождения положения КА необходимо измерить третье зенитное расстояние и построить третью коническую поверхность, однако выбор из двух направлений того, которое соответствует положению КА, обычно не представляет затруднений.
При r << да для реализации этого состава измерений необходимо располагать бортовыми устройствами, воспроизводящими направление радиус-вектора r иосуще-ствляющими пеленгацию звезд. При полетах вблизи небесных тел такими устройствами могут являться приборная вертикаль и автоматический секстант. При этом все измерительное устройство целесообразно осуществлять в виде единого блока, позволяющего избежать промежуточных вычислений по приведению составляющих r и s0 к единой системе отсчета.
Чувствительность измеряемого угла в к изменению координат x, y, z определяется производными, полученными дифференцированием уравнения (1):
cos в-a01 Ir sin в,
dx"
ß = arccos —
r i 2 . 2 . 2 V “
(x + y + z )
где x, y, z - составляющие радиус-вектора r в системе OXYZ , a°, b°, c° - составляющие единичного вектора
% = (y C0Sß-b°УГSinß’ % = ( ~ C0S ß-c° )/r Sin ß’
(2)
Связь между приращением дв угла в и малым отклонением дг радиус-вектора г определяется формулой [1]
пЪг
5ß=-
(3)
где п - единичный вектор, лежащий в плоскости измерений и перпендикулярный вектору г. Плоскость измерений образуется векторами я0 и г, проходящими через начало координат.
Измерение видимого углового диаметра Земли. Если геометрический центр сферической Земли расположен в точке, определяемой радиус-вектором гп (рис. 2), то ее видимый угловой диаметр d находится с помощью выражения
d = 2у = 2arcsin — =
Р
= 2arcsin-
Рис. 1. Схема измерения зенитного расстояния
Зная положение звезды и измерив угол в, который называют зенитным расстоянием, можно лишь указать, что КА находится на конической поверхности с вершиной в точке О, образованной вращением луча вокруг вектора я0 лучом под углом в. На одной из образующих этой поверхности лежит вектор г.
[(( - x)2 + (Уп - У)2 +(л - z)2 ]
(4)
где ^п - радиус Земли; хп, уп, гп - его составляющие по осям X, Y, Z ; р = гп - г - радиус-вектор, определяющий взаимное положение центров КА и Земли; р - модуль радиус-вектора.
Если начало системы координат OXYZ совмещено с геометрическим центром О Земли, то гп = 0 и выражение (4) принимает вид
r
d = 2агс8Іп
(5)
§у = -
0,95
5г =
0,95
(Г2 - Я ) Я
погрешности от сжатия земного сфероида
§у = ■
0,243Я2е2
5г =
0,243Я2е2
дd/ ду = -2
у-
~Уп d
^ г
д^д1 = -2 ^.
Р 2
, тОг
оу =----------tgy ,
(7)
I 2 . 2 . 2 У'
(х + у + г )
При измерениях видимого углового диаметра реальных планет вследствие сжатия и отклонения рельефа от правильной сферической поверхности появляются методические погрешности, детерминированный учет которых представляет практически неразрешимую задачу.
где т - единичный вектор, направленный из центра Земли в точку наблюдения (вдоль вектора г).
Измерение положения звезды относительно линии видимого горизонта. Пусть начало системы координат OXYZ совмещено с центром сферической Земли О, имеющим радиус Яп (рис. 3). Обозначим: г - радиус-вектор положения КА; я0 - единичный вектор направления на бесконечно удаленную звезду; в - угол между г и я0 в точке наблюдения.
Тогда угол Ф между направлением на звезду и линией видимого горизонта определится выражением
Ф = л-в~У, (8)
где в - определяется равенством (1), а у выглядит следующим образом:
. *п
у = агсзіп— = агсзіп
г
Я
п/ Г
у
(9)
Рис. 2. Схема измерения видимого углового диаметра Земли
Однако при наличии некоторых усредненных характеристик распределения рельефа такой учет может быть произведен в статистическом плане. Так, например, в результате обработки данных об истинных фигурах Земли и Луны могут быть получены статистические оценки погрешностей измерения видимых угловых радиусов 5у этих небесных тел, определяемых расстоянием до их центров 5г [1]. С вероятностью 90 % эти погрешности имеют следующие значения: погрешности от отклонений рельефа поверхности Земли относительно уровня сфероида (без учета рефракции)
1 / 2 -п2\!
г (г - ЯЭ)
Следует иметь в виду, что ошибки измерения угла у вблизи Земли в реальных условиях будут значительно выше приведенных за счет размытости самого горизонта светящейся атмосферой.
Производные, характеризующие чувствительность угла d к изменению координат, имеют вид
д<і/дх = -2 —,
Р2 2
(6)
Если центр Земли совмещен с началом системы координат, то в формулах (6) необходимо учесть, что р = г, хп = уп = гп = 0. В общем случае для сферической Земли радиусам Яп, связь между приращениями 0г вектора г иугла 5у = 5(d/2) определяется формулой
Рис. 3. Измерение положения звезды относительно видимого горизонта
В большинстве случаев звезды Земли не являются сферическими, и для точных измерений необходимо было бы учитывать истинное значение Яп в точке горизонта под звездой. Однако такой учет ведет к значительному усложнению всех вычислений. Поэтому расчет по формуле (9) допустим только для измерений, производимых на значительных удалениях от планет с относительно небольшим сжатием. На очень малых удалениях от поверхности Земли некоторое повышение точности может быть достигнуто в результате учета радиуса Земли в точке пересечения радиус-вектора г с ее поверхностью. Так, например, если высота полета КА над поверхностью Земли составляет 150-180 км, то при совмещении оси Z с полюсом мира радиус Земли в горизонте может быть принят равным ее радиусу Яп в подспутниковой точке. Если при этом считать, что фигура Земли соответствует эллипсоиду вращения с экваториальным радиусом Я3 и сжатием е, то радиус подспутниковой точки имеет вид
г
откуда
Y = arcsin
(x2 + y2 + z2)
(10)
1
дФ =— (n - mtgy )r,
(11)
-Г
Связь между приращениями радиус-вектора г и измеряемого угла Ф может быть найдена дифференцированием уравнения (8).
дФ _ дв ду
В результате получим ~х __"дХ~~дХ ’
дФ __5р_дУ дФ __5р_5у где производные
ду ду ду ’ дг дг дг ’
5р 5р 5р вычисляются по формулам (2).
дх ’ ду ’ дг
ду
Производные — вычисляются по формулам (6) с
дд
заменой р на г. Малые приращения угла Ф связаны с приращениями вектора г уравнениями
где п - единичный вектор, лежащий в плоскости, образуемой векторами г и я0, и перпендикулярный вектору г; т - единичный вектор в точке наблюдения, направленный по вектору г. Таким образом, данный состав измерений по своим информативным свойствам эквивалентен одновременному измерению зенитного расстояния в звезды и углового диаметра d Земли. Методические ошибки, обусловленные сжатием и рельефом Земли, при измерениях угла Ф могут быть учтены по тем же формулам, что и для измерений видимого углового диаметра d Земли.
Измерение угла между звездой и известным ориентиром на поверхности Земли. Совместим начало системы координат OXYZ с геометрическим центром Земли О (рис. 4). Пусть г -радиус-вектор положения КА; гоп -радиус-вектор, определяющий положение ориентира на поверхности Земли; я0 - единичный вектор, направленный на бесконечно удаленную звезду.
Тогдаугол ß между векторами я0 и р = r - гоп определится выражением
ß = arccos У-Х-К + (y-y.И + (z-z-pК . (12)
[(X - Хор К +(y - Уор К + (z - Z* )2 ]
Производные, определяющие связь между малыми приращениями угла ß и координат x, y, z, имеют следующий вид:
dß/dx = [(x - Хор) cos ß - «оР]/Р2 sin ß, dß/5y = [(У - Уор)cos ß - boP]/Р2 sin ß 3ß/3x = [(z - Zор )cos ß-Cop]/p2sln ß, (13)
где p -радиус-вектор, направленный из точки наблюдения к точке на поверхности Земли. Связь между малыми приращениями 5ß и 5r определяется выражением
dß = щЪг/р , (14)
где n1 - единичный вектор, перпендикулярный вектору р и лежащий в плоскости, образуемой векторами я0 и р.
Возможности использования для навигации некоторых физических полей. Помимо электромагнитного и гравитационного полей астрономические объекты имеют другие поля, которые принципиально могут явиться источниками ПНП. К их числу относятся магнитное и электростатическое поля, а также поле атмосферы. Однако к настоящему времени только магнитное поле Земли изучено в такой степени, что оно может быть практически использовано для навигации в космосе.
В первом приближении магнитное поле Земли может быть представлено в виде поля однородного намагниченного шара или поля центрального диполя. Начало вектора магнитного момента Мэтого поля совмещено с геометрическим центром эллипсоидальной Земли, и он наклонен под углом 11°30' к ее оси вращения. Вектор напряженности T этого поля (как поля диполя) описывается уравнением
T =
rM M - 3 — r
r2
(15)
Рис. 4. Измерение положения звезды относительно ориентира на поверхности
где г - радиус-вектор точки наблюдения в геоцентрической системе координат.
Фактическое магнитное поле имеет более сложную структуру, причем меняющуюся во времени. Причинами отклонения реального поля Земли от дипольного и его изменения во времени являются мировые и региональные аномалии, изменение интенсивности и распределения токовых систем внутри тела Земли и токовых систем в атмосфере.
При среднем значении напряженности магнитного поля Земли у ее поверхности, равном 3 -104 у (У -гамма-единица, соответствующая 105Э ), вековые вариации магнитного момента от изменения внутренних токовых систем достигают 20 у в год. При этом систематический дрейф момента к западу составляет 0°15' в год.
Атмосферные токовые системы в спокойное время создают магнитное поле интенсивности 20 у ...70 у на расстоянии трех радиусов Земли (г « 20 000 км) и 10 у ...20 у на ее поверхности. Это поле подвержено наиболее ин-
тенсивным и быстрым изменениям, обусловленным процессами в атмосфере и активностью Солнца. При этом могут наблюдаться глобальные возмущения, охватывающие всю Землю, с изменениями интенсивности в несколько тысяч гамм при продолжительности в несколько суток; месячные возмущения с изменением интенсивности в сотни гамм при продолжительности в несколько часов; наконец, высокочастотные флюктуации интенсивности в несколько гамм и с периодом до нескольких минут. Статистические характеристики всех возмущений магнитного поля Земли в настоящее время достаточно хорошо изучены. Это позволяет создать математическую модель этого поля, относительно точно описывающую реальное поле.
Математическое ожидание Т0 вектора напряженности как среднее по множеству многолетних измерений, описывается разложением его составляющих Тх, Ту, Тг в геоцентрической системе по шаровым функциям [1]:
т‘=££ т, =И
(k coskX + i" sinkx) j—j + (" coskX + e" sinkx)j -r-
( sin(X — П coskxj|^—J + (k sinkX — e" coskx))
d \^P" (cosQ)]
dQ
d \^P" (cos Q)]
Tx = I (k coskX + i" sinkx)" + 1))—"j + (" coskX + e" sinkx)kf — 1 |—[ " (-----), (16)
Г1!d [p(cos0)]
R J J dQ ’
где Rз - радиус сферической поверхности Земли; X -географическая долгота; 0 = 90° — ф (ф - географическая широта); Pk (cos 0) - присоединенные полиномы Лежандра 1-го рода; i", i", E", e" - нормированные коэффициенты, зависящие от внутренних и атмосферных источников магнитного поля Земли. В настоящее время эти коэффициенты определены с относительно высокой точностью.
Первая гармоника этого разложения (" = 1) соответствует полю диполя. Точность знания входящих в уравнения (16) коэффициентов позволяет вычислять поля с учетом до 10-й гармоники (" = 10). При этом расчетные значения составляющих Г0 отличаются от фактически измеренных не более чем на 0,1...1,5 %при r < R3. С увеличением удаленности более 5 R3 ошибки быстро возрастают.
Помимо математического ожидания при расчете магнитного поля Земли необходимо учитывать случайные составляющие 5T, в большинстве своем носящие нестационарный характер как во времени, так и в пространстве. В настоящее время с помощью магнитометрической аппаратуры уже накоплена значительная информация о статистических характеристиках магнитного поля Земли на различных высотах от ее поверхности и о вариациях этих характеристик во времени. Это создает определенные предпосылки к тому, что в будущем при достаточном уточнении таких характеристик могут быть созданы условия для практической реализации навигационных систем искусственных спутников Земли, основанных на магнитометрических измерениях.
В основу построения таких навигационных систем могут быть доложены две основные схемы: с измерением полного вектора напряженности T и измерением лишь его модуля.
В первом случае магнитометрический измерительный блок, состоящий из трех магнитометров с ортогонально расположенными осями чувствительности, дол-
жен устанавливаться на площадке, стабилизированной в инерциальном пространстве или системе координат, связанной с Землей. Оси магнитометров ориентируются (с точностью до инструментальных и методических ошибок) по осям инерциальной или связанной с Землей системами координат. При этом каждый магнитометр измеряет составляющую вектора Т по соответствующей оси. Если бы магнитное поле было детерминированным и описывалось уравнением (13), знание составляющих вектора позволило бы однозначно определять место его измерения, т. е. находить координаты КА. Однако уже из-за регулярных вариаций, учитываемых в уравнении (14), по одному измерению место измерения определяется неоднозначно, поскольку вектор Т по величине и направлению может иметь одно и то же значение во многих точках. Случайные вариации приводят к тому, что количество точек становится весьма большим и даже может стремиться к бесконечности. Техническая реализация навигационной системы, базирующейся на измерениях полного вектора Т, является весьма сложной проблемой из-за необходимости материализации стабилизированной платформы. Особенно трудно создать на подвижном объекте платформу, связанную с фиксированной относительно Земли системой координат.
Навигационные системы, базирующиеся на измерениях модуля вектора напряженности Т, технически реализуются гораздо проще, так как они не требуют наличия стабилизированной платформы. Измерительный блок, состоящий, как и в первом случае, из трех магнитометров с ортогонально расположенными осями, может быть установлен на корпусе объекта. В любой момент времени и при любой ориентации измерительного блока модуль вектора Т рассчитывается по формуле
Т _( + Ту + Т/ )1/2, (17)
где Тх, Ту, Тг - измеренные составляющие вектора Тпо осям измерительного блока.
Однако определение места по измеренным значениям модуля вектора становится весьма сложной задачей, так как любому его значению соответствует поверхность весьма сложной конфигурации. Неоднозначность определения места как по измерениям полного вектора напряженности Т , так и его модуля делает неизбежным вывод, что навигационные магнитометрические системы могут работать только в режиме уточнения места наблюдения. При этом предварительная оценка места должна быть достаточно точной, точность знания магнитного поля - весьма высокой, а навигационная задача должна решаться статистическими методами, по результатам обработки большого количества измерений.
Исходя из выше сказанного, использование астрономических источников первичной информации для автономного определения движения объектов является весьма актуальной для космической отрасли.
Библиографическая ссылка
1. Портфирьев Л. Ф., Смирнов В. В., Кузнецов В. И. Аналитические методы оценки точности автономных методов определения орбит. М.: Машиностроение, 1987.
I. N. Kartsan, V. S. Ponomarev, E. S. Zhukova, S. V. Litoshik
ASTRONOMICAL SOURCES OF PRIMARY NAVIGATION INFORMATION FOR AUTONOMOUS DEFINITION OF THE OBJECT MOTION
Akalyzed astrokomical sources of primary "avigatio" ikformatio" for autokomous determikatio" of the motio" of objects.
Keywords: a"tiaircraft dista"ce of a star, compressio" of a terrestrial spheroid, a" i"te"sity vector.
© Карцан И. Н., Пономарев В. С., Жукова Е. С., Литошик С. В., 2011
УДК 629.7.017.1
А. В. Кацура, В. А. Лавренов, А. А. Рябин
ПРИМЕНЕНИЕ МЕТОДОВ НЕРАЗРУШАЮЩЕГО КОНТРОЛЯ ДЛЯ ВЫЯВЛЕНИЯ КОРРОЗИОННЫХ ПОРАЖЕНИЙ ЭЛЕМЕНТОВ КОНСТРУКЦИЙ ЛЕТАТЕЛЬНЫХ АППАРАТОВ
Использование этих методов на различных этапах эксплуатации обеспечивает достоверность оценки технического состояния летательных аппаратов (ЛА).
Ключевые слова: техническая диагностика, неразрушающий контроль, коррозионное поражение.
Применение методов неразрушающего контроля для исходного материала, заготовок, полуфабрикатов и готовых деталей на этапах производства авиационной техники способствует выпуску продукции высокого качества. Использование методов неразрушающего контроля на различных этапах эксплуатации при выполнении периодических и целевых технических осмотров обеспечивает предотвращение лётных происшествий, полную выработку назначенных и межремонтных ресурсов и их обоснованное увеличение при сохранении заложенного уровня надёжности, повышение достоверности оценки технического и коррозионного состояния высоконагруженных элементов конструкции.
Контроль технического состояния выполняется для получения информации о фактическом состоянии летательного аппарата.
Внезапный отказ, обусловленный повреждением или усталостным разрушением элемента, обычно легко выявляется при визуальном контроле.
Коррозия считается актуальной проблемой эксплуатации. В связи с этим по требованиям ИКАО - программы по сохранению целостности конструкции ЛА -обязательно включают программы по предупреждению и контролю уровня коррозии. В частности, программа коррозионной защиты содержит описание методов обнаружения, удаления и предупреждения коррозии. Основным средством обнаружения повреждения вследствие коррозии является визуальный осмотр, для недоступных зон применяются методы неразрушающего контроля [1].
Коррозионные поражения элементов, возникающие в процессе эксплуатации, требуют своевременного обнаружения и устранения. Неудаленные очаги коррозии
при дальнейшей эксплуатации развиваются, как правило, без торможения во времени.
Своевременное обнаружение коррозионных поражений при визуальных и инструментальных осмотрах зависит от возможности доступа к осматриваемой зоне и величины (относительного размера) поражения.
Схемы и рисунки элементов конструкции, которые необходимо контролировать по коррозионному состоянию при выполнении технического обслуживания, методы удаления очагов коррозии и восстановления защитных покрытий приведены разработчиком ЛА в руководстве по защите от коррозии летательного аппарата в эксплуатации.
Коррозионные поражения элементов, находящихся в легкодоступных местах, также выявляются при визуальных осмотрах.
Для осмотра труднодоступных зон необходимо использовать переносные лампы, поворотные зеркала, перископические дефектоскопы, эндоскопы и другие инструменты.
При осмотре силовых элементов конструкции с целью контроля коррозионного состояния необходимо через определенные периоды времени применять инструментальные методы неразрушающего контроля.
Для обеспечения полного использования возможностей методов неразрушающего контроля должны выполняться следующие условия:
1) в процессе разработки летательного аппарата и эксплуатационно-технической документации должно быть предусмотрено применение методов неразрушающего контроля для проверки технического состояния частей и узлов ЛА;
2) при проектировании летательного аппарата необходимо предусмотреть контролепригодность конструк-
