CC BY
134
УДК 62.5; 519.6
СИНТЕЗ СИСТЕМЫ УПРАВЛЕНИЯ СПУСКОМ КОСМИЧЕСКОГО АППАРАТА В АТМОСФЕРЕ МАРСА
Дивеев Асхат Ибрагимович, д.т.н.,
Северцев Николай Алексеевич, д.т.н.,
Шмалько Елизавета Юрьевна, к.т.н.
Вычислительный центр им. А. А. Дородницына РАН ул. Вавилова, 40, Москва, Россия, 119333
Рассмотрена задача спуска космического аппарата в атмосфере Марса. Требования времени таковы, что следующее поколение космических кораблей, исследующих Марс, должны с высокой точностью доставлять максимально нагруженные научные лаборатории на Марс. Для синтеза оптимального управления предложен метод сетевого оператора. Управление обеспечивает минимум максимальной перегрузки в процессе спуска и попадание в заданную терминальную точку при вариации начальных условий.
Ключевые слова: спуск космического аппарата, синтез системы управления, сетевой оператор
Исследование и осмысление основных тенденций спуска в атмосфере Марса и последующей посадки космического аппарата (КА) на поверхность планеты является основополагающим шагом к последующему повышению качества полетов к "красной" планете и, быть может, высадки людей на поверхность Марса.
Как показала практика, 66% предпринятых миссий на Марс оказались провальными. И, несмотря на то, что вход КА в атмосферу, торможение и посадка его на поверхность Марса длятся в общей сложности около 6 минут (а весь полет составляет порядка двух лет), данный этап является самым сложным с точки зрения безопасности и точности полета.
Очевидно, что на сегодняшний день не создано достаточно надежной системы спуска и посадки КА на Марс, позволяющей безопасно в смысле перегрузки и целостности аппарата и с требуемой точностью посадить аппарат в запланированное место посадки. Актуальным так же остается вопрос об увеличении массы спускаемого аппарата, что представляет собой сложную задачу в условии крайне высокой разреженности атмосферы Марса. [1]. Кроме того, параметры атмосферы Марса значительно меняются в зависимости от времени года, времени суток, наличию пылевых бурь и др.
Значительная удаленность Марса в свою очередь затрудняет управление с Земли процессом спуска космического аппарата. Все вышеуказанные трудности, а также высокая стоимость приближенного моделирования марсианской среды накладывают определенные ограничения на исследования в области спуска на Марс и выработку оптимальных управлений.
Одной из основных проблем, являющихся причиной неудачных посадок заключается в отсутствии точных знаний о начальном состоянии спускаемого космического аппарата. Траектория спуска и его основные параметры, точность посадки и максимальная перегрузка, существенно зависят от начального угла входа в атмосферу и начальной высоты.
Преодолеть указанные трудности в данной работе предлагается путем ухода от построения оптимального управления как функции времени, так как оно не учитывает поведения системы уже в процессе функционирования и влияния этого поведения на дальнейшее состояние всей системы. В свою очередь предлагается использовать синтез системы управления как функции от состояния объекта. Именно такое поведенческое описание систем в динамике является приоритетным на сегодняшний день, поскольку максимально учитывает работу объекта в реальности.
Математическая модель спускаемого космического аппарата представляет собой систему дифференциальных уравнений четвертого порядка.
dx1
dt
— x
2,
(1)
dx2
dt
XX2 V
YX4
V
cos(u)
g 0 RM X1
(x2 + x3 f
dx3
dt
— x
4 ,
(2)
(3)
dx4 — — Xх4 + Y —cos(u) dt V V
g 0 rm x3 (x2 + x2 )3/2
(4)
где x1, x3 - координаты центра масс космического аппарата в ортогональной системе координат, x2, x4 - компоненты скорости космического аппарата, и
- управление, и — < и £ и +, Rm - экваториальный радиус Марса, g0 -ускорение свободного падения на поверхности Марса,
V=4
xf + x| , X =
CxPV2S Y = CyrV2S
2m 2m
m - масса космического аппарата, S - площадь поверхности сопротивления комического аппарата, Сx, Сy - аэродинамические коэффициенты,
Р = Р 0 е,
r0 - плотность атмосферы на поверхности Марса, l - коэффициент разреженности атмосферы,
Н — x 2 + x3 — Rm .
Заданы начальные xi(0) — 0 условия: xi(0) — 0, x2(0) = V0 cos( J(0)), x3(0) — H(0) + Rm , x4(0) — V0 sin( J(0)),
где J(0) - угол наклона траектории при входе в атмосферу,
J(0) = Jj, Jj = Jо + jDJ, j є {0,1,..., Kj}, (5)
H (0) - расстояние до поверхности планеты, с которого начинается вход в атмосферу,
H (0) = H, Hi = H 0 + (2i - 1)DH , і є {0,1,..., Kh }. (6)
Цель управления описывается следующими функционалами:
J = max G + bE(tf)®min, (7)
0<t £tf J
где b - весовой коэффициент, G - перегрузка
G =
V + V} + V}
g 0
V
x
V
y
-XV - У^щ)-7^0ЩЇ
V V (x2 + x2 )32
Xх4 + Y^ccs^) g0RM}
V V
(x12 + x3 )
32
E (t f) - терминальная ошибка
Vz = У sin(u),
E(tf) = L(tf) - Lf + H(tf) - Hf L = RM arctg
f \ X1
w
V x3
h=+X7—Rm
Управление ищем в виде
и = j(L, H ,V, G, q), (8)
где j(x, q) - искомое математическое выражение, q - вектор параметров, не
меняющихся в процессе управления.
Особенностью рассматриваемой задачи является наличие множества начальных значений (5), (6). При различных начальных значениях (5), (6) значения функционала (7), различны.
Система управления (9) должна для каждого начального значения из множеств (5), (6) обеспечивать такое управление, которое бы давало оптимальные значения функционалов, т.е. такие значения, которые могли бы быть получены при решении задачи оптимального управления.
Для решения такой задачи синтеза используем два этапа. На первом этапе решаем задачи оптимального управления для каждого начального значения. Получаем массив оптимальных значений функционалов J (i, j) , где i, j определяет начальные значения (5). (6).
На втором этапе решаем задачу синтеза (8) с функционалом
J
K J к
Z Z
j=0 i=0
H J bjzJihA J(i, j)
® min
(9)
При решении задач оптимального управления используем метод аппроксимации кривыми Безье [2]. При решении задачи синтеза используем вычислительный метод сетевого оператора [3,4].
В таблице 1 приведены результаты решения задачи оптимального
управления при различных начальных значениях.
При моделировании использовались следующие числовые значения параметров: g0 = 3,71 • 10-3 км / с2, RM = 3396 км, р0 = 0,02 кг / м3, 1 = 0,093, S = 9,61625 м2, m = 930 кг, Сх = 1.68, Су = 0.3024, V0 = 4.5 км / с, Hf = 10 км, Lf = 600 км, t + = 300 с, u- = -p, и + = 0, b = 0.25.
Таблица 1.
J0 H 0 max G E (tf)
- 0,28 155 18,7794 20,82
- 0,28 160 18,5722 43,27
- 0,28 165 18,3602 61,41
- 0,3 155 21,0018 32,60
- 0,3 160 20,8197 12,27
- 0,3 165 20,6339 7,26
- 0,32 155 23,1743 77,62
- 0,32 160 23,0138 59,32
- 0,32 165 22,8499 40,88
Литература
1. Braun R.D., Manning R.M. Mars Exploration Entry, Descent and Landing Challenges//IEEE Aerospace Conference, Big Sky Montana, March 2006. Paper №0076.
2. Дивеев А.И., Шмалько Е.Ю. Метод аппроксимации кривыми Безье для решения задачи оптимального управления посадкой космического аппарата// Труды ИСА РАН. Динамика неоднородных систем/ Под ред. чл.- корр. РАН Ю.С. Попкова. М.: ИСА РАН, КомКнига, 2007. Том 31(1). С. 8-13.
3. Дивеев А.И., Шмалько Е.Ю. Многокритериальный структурнопараметрический синтез системы управления спуском космического аппарата на основе метода сетевого оператора// Вестник Российского университета дружбы народов. Серия инженерные исследования (информационные технологии и управление). 2008, № 4. С. 86 - 93.
4. Дивеев А.И. Метод сетевого оператора. М.: Изд-во ВЦ РАН. 2010. 178 с.
SYNTHESIS OF CONTROL SYSTEM FOR SPACECRAFT DECENT IN
MARTIAN ATMOSPHERE
A.I. Diveev, N.A. Severtsev, E.Y. Shmalko
Spacecraft decent in Martian atmosphere is considered. The next generation of Mars exploration landers must precisely deliver massive scientific payloads to sites of interest. Network operator method is used for optimal control synthesis. It delivers minimum to maximal G-force and landing in the defined area allowing some variations in initial conditions.
Key words: spacecraft landing, synthesis of control system, network operator
