CC BY
111
УДК 623.54:623.451.08
КОРРЕКЦИЯ ТРАЕКТОРИИ ПЕРЕХОДА КОСМИЧЕСКОГО АППАРАТА НА КРУГОВУЮ ОРБИТУ СПУТНИКА ПРИ ИСПОЛЬЗОВАНИИ ТОРМОЖЕНИЯ В АТМОСФЕРЕ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ ЕЕ ПАРАМЕТРОВ
В.П. Казаковцев, В.В. Корянов, Зо Мин Тун
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация
e-mail: kafsm3@bmstu.ru; vkoryanov@mail.ru; eaglez.lion@gmail.com
Современный этап развития мировой космонавтики характеризуется усложнением решаемых в условиях космического пространства задач, требующих для этих целей применения достаточно громоздких вспомогательных систем и средств. Значительная масса топлива расходуется при переходе космического аппарата с пролетной траектории на орбиту спутника планеты. Использование промежуточной эллиптической орбиты, у которой перицентр находится в верхних слоях атмосферы планеты, позволяет уменьшить расход топлива. Рассмотрен вопрос использования атмосферы планеты для торможения космического аппарата при его переходе на орбиту спутника в условиях неопределенности знания параметров атмосферы.
Ключевые слова: орбита, пространственное движение, космический аппарат, торможение в атмосфере, тепловой поток, скоростной напор.
CORRECTION OF SPACECRAFT TRAJECTORY
DURING ITS TRANSITION INTO CIRCULAR SATELLITE ORBIT
USING ATMOSPHERIC BRAKING WITH UNCERTAIN ATMOSPHERIC
PARAMETERS
V.P. Kazakovtsev, V.V. Koryanov, Zaw Min Tun
Moscow State Technical University, Moscow, Russian Federation. e-mail: kafsm3@bmstu.ru; vkoryanov@mail.ru; eaglez.lion@gmail.com
The current stage of space exploration is characterized by complication of the problems being solved in space that require quite cumbersome auxiliary systems and tools. A considerable mass of fuel is consumed during the spacecraft transition from the planet gravitation assist trajectory into the planet satellite orbit. The use of an intermediate elliptical orbit with the pericenter located in the planet's upper atmosphere can reduce the fuel consumption. The article is devoted to the use of the planet's atmosphere for slowing down the spacecraft during its transition into the satellite orbit with uncertain atmospheric parameters.
Keywords: orbit, spatial motion, spacecraft, braking in atmosphere, heat flow, dynamic pressure.
Обсуждение проблемы. В работах [1-3] показано, что использование торможения космического аппарата (КА) в верхних слоях атмосферы планеты при его переводе на орбиту спутника приводит к экономии топлива на 20-30%. Однако реализация движения КА в верхних слоях атмосферы сталкивается с большими трудностями. В работах [4, 5] рассмотрены некоторые проблемы стратегии и реализации торможения КА в верхних слоях атмосферы. Приведены особенности
38 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
такого торможения в космических миссиях Odyssey и Mars Global Surveyor. Торможение скорости КА осуществлялось в основном за счет энергии солнечных батарей. При этом максимальная допустимая температура нагрева поверхности солнечных батарей не должна была превышать 175 °С. Температура поверхности солнечных батарей во многом определяется суммарным количеством теплоты, подведенным к ней за время движения КА в верхних слоях атмосферы. Суммарное количество теплоты является интегральной характеристикой от теплового потока. Кроме того, на панели солнечных батарей действует силовая нагрузка, пропорциональная скоростному напору. Поэтому в процессе движения КА в верхних слоях атмосферы необходимо контролировать как тепловой поток, так и скоростной напор.
Неточность знания параметров атмосферы значительно затрудняет выбор высоты перицентра переходной орбиты для торможения КА [6]. Например, плотность в верхних слоях атмосферы Марса для одной и той же высоты может сильно отличатся от плотности атмосферы основной модели. В табл. 1 показаны значения плотности в атмосфере Марса, взятые из работы [7].
Таблица 1
Значения плотности для различных моделей атмосферы Марса
Высота H, км Плотность атмосферы, кг-м 3
основная максимальная минимальная
140 1,27 • 10-9 2,03-10-8 1,1710-10
120 V О О 1 со 2,24-10-7 1,29-10-9
100 1,54-10-7 2,47-10-6 1,42-10-8
80 О О 1 Об 2,72-10-5 1,57-10-7
Рассмотрим основные моменты движения КА, перицентр орбиты которого находится в атмосфере, а апоцентр — на достаточно большом расстоянии от планеты.
Математическая модель движения КА. Выберем математическую модель пространственного движения КА в соответствии с работами [8-10], сформированную при следующем допущении: планета постоянного радиуса (R) вращается с постоянной угловой скоростью (ш):
—^=тш2(cos sin в— sin ^ cos в sin ф) cos ^—
— gR sin в + дт cos в sin ф — Ббq;
de
dt
1
V
тш2 (cos cos в + sin sin в sin ф) cos ^ — gR cos в—
— gT sin в sin ф + 2Vw cos ^ cos ф
-1 V cos в
+---------
т
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 5 39
dф dt
1
V cos в
ru2 cos p sin в sin ф+
+ 2Vu(cos p sin ф sin в — sin p cos в)+ gT cos ф
-1 V tg p cos в cos ф r
dh т/ • д
— = V sin в;
dt
dp V cos в sin ф dt r
dA V cos в cos ф dt
(1)
r cos p
где F — скорость КА; r — радиус; A — высота; в — угол наклона вектора скорости; ф — угол поворота траектории; p, A — широта и долгота положения КА; q = 0,5pV2 — скоростной напор; р — плотность атмосфе-
ры; S — площадь миделевого сечения; S6 =
SCx
m
— баллистический
коэффициент; Cx — аэродинамический коэффициент силы лобового
ГА ,1 + 0,002948(1 — 3 sin2 p)R2
сопротивления; m - масса КА; gR = р
0,005126 cos2 pr\ 2
)r2 — радиальная составляющая ускорения тяго-
R2
тения; gT = — р
(0,002948 sin2 pR2 0,002563 sin2 pr2
R2
r2 — транс-
версальная составляющая ускорения тяготения; р — гравитационный параметр поля тяготения планеты.
Определение условий учета неопределенности параметров атмосферы. Температура поверхности тела при наличии отвода теплоты излучением определяется из условия баланса теплоты [11]:
qw — eoTW = qx,
(2)
2
где qw — тепловой поток, подводимый к поверхности от газа; qx — тепловой поток, отводимый внутрь материала (во многих задачах при оценке теплового баланса считают qx = 0); Tw — температура поверхности материала; е — излучательная способность материала; о — постоянная Стефана-Больцмана.
При движении КА со скоростью порядка 3,5... 4,5 км/с в атмосфере плотностью 10-7... 10-10 кг-м-3 можно считать подводимый тепловой поток конвективным. Для конвективного теплового потока в критической точке поверхности, согласно работе [12], можно записать следующее выражение:
АконЯ
утн
(3)
40 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
где AkohS = (4,5... 5,3) • 104, rH = 1, p0 ,V±k — начальные значения плотности и скорости; n = 0,5, m = 3,0 ... 3,2.
В настоящей работе будем использовать формулу (3) для оценки подводимого теплового потока в следующем виде:
qT = k± p0,5 V3. (4)
Рассмотрим пример расчета движения КА по орбите вокруг Марса при следующих параметрах: высота перицентра ha = 100 км, высота апоцентра hA = 500 ... 10000 км, модель атмосферы — основная. Результаты расчетов приведены на рис. 1 и 2.
Отметим, что максимальные значения скоростного напора qmax и теплового потока qt max возникают одновременно в момент прохождения КА перицентра орбиты. Изменение qmax соответствует изменению
0,5 1,5 -3 -1 -2
qt,KT м -с ;д,кг-м с
0,8 0,6 0,4 0,2
0 100 200 300 U с
а
0,5 1,5 -3 -1 -2
qt, кг • м • с ; q, кг -м -с
0,4 0,3 0,2 ОД
0 200 400 600 800 U с
6
Рис. 1. Графики изменения относительного значения теплового потока qt и апоцентра hA = 10000 (а) и 500 км (б):
1 — максимальные значения теплового потока qt max и 2 — максимальные значения скоростного напора qmax
-1 -2
4тах>КГ-М -С
2.0 I------------------------------------------
1,5 - ^____________________________
1.0 - ^ '
0,5 -
_______I_______I_____I_______I_______I______
0 2000 4000 6000 8000 10000 hA, км
Рис. 2. График изменения максимального скоростного напора qmax в зависимости от высоты апоцентра орбиты hA при ^п = 100 км
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. №5 41
qtmax- Таким образом, контролируя максимальный скоростной напор одновременно контролируем максимальный тепловой поток.
При фиксированной высоте перицентра наибольший максимальный скоростной напор будет в начале движения КА по переходной орбите. Поэтому выбирать начальное значение высоты перицентра надо по параметрам орбиты с максимальным значением высоты апоцентра.
При наличии акселерометра, расположенного по продольной оси КА, определяется осевая перегрузка. По максимальному значению перегрузки вычисляется максимальное значение скоростного напора при каждом витке прохождении КА атмосферного участка орбиты. Это значение qmax сравнивается с заранее рассчитанным допустимым значением максимального скоростного напора qmaxfl0n и таким образом определяется отклонение максимального скоростного напора от его допустимого значения
^qmax
qmax qmaxдоп•
(5)
Если Aqmax > 0, то определяется значение корректирующего импульса A Va в апоцентре последующего витка орбиты. Заранее по параметрам номинальной траектории и для основной атмосферы рассчитываются
/ 9К \
элементы матрицы частных производных —---- для каждого витка
Vdqmax' i
орбиты КА. Далее определяется изменение высоты движения КА в перицентре орбиты
A К
дЪп_
^qmax
Aqmax-
(6)
Также заранее рассчитываются элементы матрицы частных произ-для каждого витка переходной орбиты. Корректирую-
i
щий импульс скорости в апоцентре рассчитывается по формуле
водных
dVA
A
dh
AVa = djA Ahn. dhn
(7)
Методика определения корректирующих импульсов скорости в апоцентре переходной орбиты в условиях неопределенности параметров атмосферы планеты. 1. Выбираем параметры номинальной орбиты (высоты перицентра и апоцентра), которые обеспечивают заданное время перехода КА на требуемую орбиту спутника планеты.
2. На основании расчетов параметров теплового нагрева панелей солнечных батарей для скоростей в перицентре номинальной орбиты определяем максимальные значения теплового потока и скоростного напора. Выбираем допустимое значение максимального скоростного
напора qmaxдоп-
42 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
3. По параметрам номинальной орбиты и основной атмосферы
dhn
планеты для каждого витка рассчитываем элементы матриц —--- и
dVA
дд„
A
dh-п
4. Для атмосферы планеты с максимальной плотностью по значению дтахдоп уточняем начальное значение высоты перицентра hH1 для первого витка движения КА по орбите.
5. На каждом i-м витке орбиты определяем максимальный скоростной напор qmax *. По формуле (5) рассчитываем Agmax а
6. Вычисляем по формуле (6) необходимое изменение высоты перицентра A hnA
7. По формуле (7) определяем корректирующий импульс скорости в апоцентре орбиты — AVAi и проводим коррекцию орбиты.
8. Повторяем действия по п.п. 5-7 до момента достижения заданной высоты апоцентра, равной высоте требуемой орбиты спутника планеты.
Таким образом, в процессе движения КА по переходной орбите на каждом витке корректируется высота перицентра по допустимому максимальному скоростному напору и тем самым выполняется температурный режим на поверхности панелей солнечных батарей.
Пример использования методики учета неопределенности параметров атмосферы. Рассмотрим пример расчета коррекций высоты перицентра орбиты КА при его переводе на круговую орбиту спутника Марса с использованием торможения в атмосфере. Расчет тепловой нагрузки на панели солнечных батарей, которые используются для торможения в атмосфере Марса, показал, что значение скоростного напора в перицентре должно удовлетворять условию дпв < qmax = = 1,49 ккал-м-2с-1. Это верхний предел зоны допустимого максимального скоростного напора. Выберем нижний предел зоны, удовлетворяющей нашей задаче: дпн = 1,20 ккал-м-2с-1.
Номинальная орбита для основной атмосферы Марса: hn = 100 км, hA = 10000 км.
Предположим, что реальная плотность атмосферы несколько больше и на высоте 100 км составляет р100 = 2 • 10-7. Поэтому выбираем начальное значение высоты перицентра hn0 = 116 км. Далее по предложенной методике рассчитываем корректирующие импульсы для осторожного понижения на 2 км высоты перицентра.
Результаты расчетов по предложенной методике приведены в табл. 2 и на рис. 3
На первом витке дп1 < дпн, поэтому вычислен и проведен корректирующий импульс для понижения на 2 км высоты перицентра. Далее до пятого витка проводились корректирующие импульсы, пока значение
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 5 43
Результаты расчета коррекций орбиты
Номер витка АУаь ,мс-1 Ki, км — 1 — 9 qm, кгм 1 •с 9
1 0 116 0,516
2 0,26 114 0,656
3 0,26 112 0,833
4 0,26 110 1,06
5 0,26 108 1,35
6 0 108 1,35
7 0 108 1,35
Рис. 3. Диаграмма изменения высоты перицентра орбиты в функции от номера витка
скоростного напора не вошло в заданную зону. Таким образом, для данной случайным образом выбранной атмосферы проведены корректирующие импульсы. Требуемая высота перицентра оказалась равной 108 км.
Выводы. 1. Разработана методика коррекции параметров переходной орбиты КА с использованием торможения в атмосфере планеты для случая неопределенности знания параметров атмосферы.
2. Методика позволяет учитывать влияние торможения в атмосфере на тепловые потоки, температуру и скоростной напор, определять необходимые корректирующие импульсы скорости в апоцентре орбиты для сохранения температурного режима на поверхности панелей солнечных батарей.
3. Анализ тепловых потоков и скоростных напоров на этапе движения КА в верхних слоях атмосферы показал адекватность изменения их максимальных значений. Это позволило определять значения корректирующих импульсов только по скоростному напору в перицентре орбиты.
4. Приведен пример расчета корректирующих импульсов скорости в апоцентре орбиты при неопределенности знания параметров атмосферы. Показана работоспособность предложенной методики.
ЛИТЕРАТУРА
1. Казаковцев В.П., Сухенко А.В. Об использовании предварительного тормозного
импульса для выведения космического аппарата на орбиту спутника Марса //
Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2005. № 4. С. 3-11.
44 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
2. Иванов В.М., Соколов Н.Л.Оптимальное управление космическим аппаратом за счет изменения вектора тяги двигательной установки при проведении межорбитальных маневров // Космонавтика и ракетостроение. 2014. № 2. С. 80-88.
3. Голомазов М.М., Литвинов Л.А., Литвинов И.А., Иванков А.А., Финченко В.С. Численное моделирование обтекания спускаемых аппаратов при входе в атмосферу планеты // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2011. № 4. С. 42-53.
4. David A. Spencer and Robert Tolson. Airobraking Cost and Risk Decisions // Journal of Spacecraft and Rockets. 2007. Vol. 44. No. 6. R 1285-1293.
5. Smith John C., Bell Julia L. 2001 Mars odyssey airobraking // Journal of Spacecraft and Rockets. 2005. Vol. 42. No. 3. Р. 406-415.
6. Казаковцев В.П., Корянов В.В. Динамика движения спускаемого аппарата с надувным тормозным устройством в атмосфере планеты // Естественные и технические науки. 2013. № 6. С. 266-270.
7. Иванов Н.М., Мартынов А.И. Движение космических летательных аппаратов в атмосферах планет. М.: Наука, 1985. 384 с.
8. Новицкий А.В., Зеленцов В.В. Одна задача маневрирования в атмосфере Марса // Управляющие системы и машины. 1999. № 1. С. 3-9.
9. Новицкий А.В., Зеленцов В.В. Оптимизация атмосферного маневра при выведении искусственного спутника Марса // Проблемы управления и информатики. 1999. №4. С. 143-157.
10. Ярошевский В.А. Вход в атмосферу космических летательных аппаратов. М.: Наука, 1988. 336 с.
11. Цаплин С.В., Болычев С.А., Романов А.Е. Теплообмен в космосе. Самара: Самарский университет, 2013. 53 с.
12. Андреевский В.В. Динамика спуска космических аппаратов на Землю. М.: Машиностроение, 1970. 235 с.
REFERENCES
[1] Kazakovsev V. R, Sukhenko A.V. Use of the Preliminary Brake Impetus for Launching Spacecraft into Orbit of a Satellite of Mars. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2005, no. 4, pp. 3-11 (in Russ.).
[2] Ivanov V.M., Sokolov N.L. Optimal control of a spacecraft due to the change of the trust vector of the propulsion system during interorbital maneuvers. Kosmonavtika i raketostroenie [Cosmonautic and Rocket Engineering], 2014, no. 2, pp. 80-88 (in Russ.).
[3] Golomazov M.M., Litvinov L.A., Litvinov I.A., Ivankov A.A., Finchenko V.S. Numerical simulation of flow past descent vehicles during planetary entry. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2011, no. 4, pp. 42-53 (in Russ.).
[4] Spencer David A., Tolson Robert. Airobraking Cost and Risk Decisions. J. of Spacecraft and Rockets, 2007, vol. 44, no. 6, pp. 1285-1293.
[5] Smith John C., Bell Julia L. 2001 Mars Odyssey Airobraking. J. of Spacecraft and Rockets, 2005, vol. 42, no. 3, pp. 406-415.
[6] Kazakovsev V.P., Karanov V.V. Dynamics of motion of descent spacecraft with inflatable break device in the planet’s atmosphere. Estestvennye i tekhnicheskie nauki [Natural and Technical Sciences], 2013, no. 6, pp. 266-270 (in Russ.).
[7] Ivanov N.M., Martynov A.I. Dvizhenie kosmicheskikh letatel’nykh apparatov v atmosferakh planet [Motion of spacecraft in planetary atmospheres]. Мoscow, Nauka Publ., 1985. 384 p.
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 5 45
[8] Novitsky V, Zelentsov V. One Task of Maneuvering in the Atmosphere of Mars. Upravlyayushchie sistemy i mashiny [Control Systems and Machines], 1999, no. 1, pp. 3-9 (in Russ.).
[9] Novitsky V., Zelentsov V.V. Optimization of atmospheric maneuver during insertion of an artificial satellite of Mars. Problemy upravleniya i informatiki [Problems of Control and Informatics], 1999, no. 4, pp. 143-157 (in Russ.).
[10] Yaroshevsky V.A. Vkhod v atmosferu kosmicheskikh letatel’nykh apparatov [The Entry of Spacecraft]. Moscow, Nauka Publ., 1988. 336 p.
[11] Saplin S.V., Bolychev S.A., Romanov A.E. Teploobmen v kosmose [Heat Transfer in Space]. Samara, Samara University Publ., 2013. 53 p.
[12] Andreevskiy V.V. Dinamika spuska kosmicheskikh apparatov na Zemlyu [Dynamics of Spacecraft Descent to the Ground]. Moscow, Mashinostroenie Publ., 1970. 235 p.
Статья поступила в редакцию 30.03.2015
Казаковцев Виктор Поликарпович — д-р техн. наук, профессор кафедры “Динамика и управление полетом ракет и космических аппаратов” МГТУ им. Н.Э. Баумана. МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Kazakovtsev V.P. — D.Sc. (Eng.), Professor, Department of Dynamics, Ballistics and Flight Vehicles Movement Control, Bauman Moscow State Technical University.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Корянов Всеволод Владимирович — канд. техн. наук, доцент кафедры “Динамика и управление полетом ракет и космических аппаратов” МГТУ им. Н.Э. Баумана. МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Koryanov V.V. — Ph.D. (Eng.), Associate Professor, Department of Dynamics, Ballistics and Flight Vehicles Movement Control, Bauman Moscow State Technical University. Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Зо Мин Тун — аспирант кафедры “Динамика и управление полетом ракет и космических аппаратов” МГТУ им. Н.Э. Баумана.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Zaw Min Tun — Ph.D. student, Department of Dynamics, Ballistics and Flight Vehicles Movement Control, Bauman Moscow State Technical University.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Просьба ссылаться на эту статью следующим образом:
Казаковцев В.П., Корянов В.В., Зо Мин Тун. Коррекция траектории перехода космического аппарата на круговую орбиту спутника при использовании торможения в атмосфере в условиях неопределенности ее параметров // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2015. № 5. C. 38-46.
Please cite this article in English as:
Kazakovtsev V.P., Koryanov V.V., Zaw Min Tun. Correction of spacecraft trajectory during its transition into circular satellite orbit using atmospheric braking with uncertain atmospheric parameters. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2015, no. 5, pp. 38-46.
46 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
