Терминальное управление субгиперзвуковой первой ступенью аэрокосмической системы

УДК 629.78 В01: 10.18287/2412-7329-2016-15-2-9-21

ТЕРМИНАЛЬНОЕ УПРАВЛЕНИЕ СУБГИПЕРЗВУКОВОЙ ПЕРВОЙ СТУПЕНЬЮ АЭРОКОСМИЧЕСКОЙ СИСТЕМЫ

© 2016 В. Л. Балакин, А. В. Ковалёв

Самарский национальный исследовательский университет имени академика С. П. Королёва

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

Аэрокосмическая система, субгиперзвуковая первая ступень, летательный аппарат, терминальное управление, алгоритм, коэффициент подъёмной силы, угол атаки, возмущения плотности атмосферы.

Введение

Одним из вариантов первой ступени аэрокосмической системы (АэрКС) для запуска космических аппаратов (КА) на околоземную орбиту рассматривается субгиперзвуковой летательный аппарат (ЛА), развивающий на активном участке движения с работающей двигательной установкой (ДУ) скорость, несколько меньшую гиперзвуковой (число Маха равно 4) [1; 2]. Отличительной особенностью данного ЛА-MPV (MIPCC-powered vehicle - аппарат с охлаждаемыми двигателями; Mass Injection Pre-Compressor Cooling - система дополнительной подачи кислорода и охлаждения двигателя), входящего в аэрокосмическую систему RASCAL (Responsive Access, Small Cargo, Affordable Launch - доступная система для запуска малых грузов по требованию, США), является пассивный (с выключенной ДУ) трансатмосферный режим полёта

до высоты 62,8 км [1]. На указанной высоте происходит отделение и дальнейший самостоятельный полёт второй (ракетной) ступени системы RASCAL.

Отметим, что в проектах аэрокосмических / авиационно-космических систем со сверхзвуковыми или гиперзвуковыми первыми ступенями, например [3 - 5], конечная высота не превышает 30 км. Поэтому возникает вопрос об управляемости ЛА на трансатмосферных высотах для обеспечения заданных конечных условий движения первой ступени.

Будем рассматривать пассивное движение ЛА от момента выключения ДУ до момента отделения от него ракетной ступени. Начальные и конечные условия рассматриваемого пассивного движения заданы и известна номинальная программа управления углом атаки (коэффициентом аэродинамической подъёмной силы), полученная для условий невозмущённого

Цитирование: Балакин В.Л., Ковалёв А.В. Терминальное управление субгиперзвуковой первой ступенью аэрокосмической системы // Вестник Самарского государственного аэрокосмического университета имени академика С.П. Королёва (национального исследовательского университета). 2016. Т. 15, № 2. С. 9-21. БО!: 10.18287/2412-7329-2016-15-2-9-21

движения. При наличии возмущений атмосферы (отклонений плотности от стандартных значений) конечные условия движения ЛА будут отличаться от заданных. Поэтому необходимо проводить коррекции угла атаки (коэффициента аэродинамической подъёмной силы) с целью уменьшения отклонений конечных условий возмущённого движения от заданных значений. Для этого требуется разработать алгоритм формирования командного управления углом атаки (коэффициентом аэродинамической подъёмной силы).

Целью данной работы является разработка, исследование работоспособности и оценка эффективности алгоритма терминального управления.

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

Рассмотрим модель пассивного движения ЛА в вертикальной плоскости, которая представляет собой традиционную систему дифференциальных уравнений, записанных в траекторной системе координат:

Программу управления коэффициентом подъёмной силы, которую можно использовать как номинальную для терминального управления при возмущённом движении, определим как кусочно-непрерывную функцию.

Примем следующую номинальную программу управления коэффициентом подъёмной силы [6], которая была предложена на основании полученных результатов решения задачи оптимизации [7]:

\С J, если t < tп ;

Сном _ I Уи1- 11

"" \Суа2, есЛи t ^ tП •

Здесь индекс «ном» означает «номинальная»; Суа1, Суа 2 - соответственно начальное и конечное значение коэффициента; /П - время переключения коэффициента со значения Суа1 на значение Суа2.

На коэффициент подъёмной силы Суа и на определяющий его угол атаки а согласно [1] наложены ограничения:

dV dt _ C -р ^ ха

de dt II V 1 f C • уа V

dh dt _ У sine.

р-У2 • S

2m р-У2 • S 2m

- g • sin в,

- g • cose

М/

Omn — а — amax'

С — C — C

Уа min Уа Уа max'

ГДе amin _-Ю°

«max _ 20'

Здесь ^ - время, секунды; V - скорость, М ; в - угол наклона траектории; И -

высота, м; Сха - коэффициент лобового сопротивления; Суа - коэффициент аэродинамической подъёмной силы (в дальнейшем подъёмной силы), М - число Маха; р = f (И) - плотность воздуха на заданной высоте, кг/ 3; £ - площадь крыла, / м

м2; g - ускорение свободного падения,

Суа min И

C определяются соответственно ми-

уа max А

нимальным amin и максимальным amax

min max

углами атаки и зависят от числа М.

Коэффициент Суа1 принимается равным значению, соответствующему окончанию активного участка движения ЛА. Это позволит при дальнейшем движении (увеличении высоты с уменьшением скорости) иметь значение нормальной скоростной перегрузки пуа, меньшее максимально допустимой nya • Для максимизации конечной скорости ЛА в конце участка пассивного движения коэффициент Суа2 примем равным нулю согласно решениям [7].

Подбором времени переключения tП можно обеспечить заданное значение конечного угла наклона траектории вк (принято равным 20° согласно [1]) и определить соответствующие ему конечные значения высоты hк и скорости Vк. В дальнейшем будем считать их заданными конечными условиями движения ЛА: вк *, h * V *

"к ' " к

Начальные условия пассивного движения ЛА имеют следующие значения: V = 1174 м/с; во = 38,4°; к = 26800 м [1].

В соответствии с результатами, полученными в [6], определим значения параметров программы управления, обеспечивающие выполнение конечных условий движения, близких к аппарату MPV:

= 0,374; t п = 14с.

Для расчётов используется универсальная модель отклонений плотности атмосферы [8]:

Р(к) = Рст (к) + ^ар(к),

стандартной атмосферой, величина ; -нормально распределённая случайная величина с нулевым математическим ожиданием и дисперсией, равной единице. Среднее квадратическое отклонение ор(к) определяется соотношением:

ар(к) = аро  ехр(-0,15Ю-3 • к),

р0

где плотность атмосферы р(И) представляется как нормальная случайная функция высоты к, математическое ожидание рст(к) вычисляется в соответствии со

где ор0 = 0,05 кг/м - средняя величина, соответствующая условиям на нулевой высоте [8].

В дальнейшем при моделировании возмущённого движения ЛА использовались два предельных значения нормально распределённой случайной величины ; : -5,5 и 5,5, соответствующих предельно «разреженной» и предельно «плотной» атмосфере.

Результаты моделирования возмущённого движения ЛА с номинальной программой коэффициента подъёмной силы приведены в табл. 1. Возмущения атмосферы представлены коэффициентом ; кк - конечная высота; Vк - конечная скорость; - конечное время; Дкк, АУк -соответственно отклонения конечной высоты и конечной скорости от заданных значений.

Таблица 1. Результаты моделирования при номинальной программе управления

; кк, м Vк, м/с 4, с Дкк, м ДVк, м/с

0 63275,10 677,99 66,56 0 0

-5,5 61471,63 722,12 63,33 -1803,47 44,12

5,5 64886,06 634,81 69,53 1610,96 -43,17

Из табл. 1 можно сделать следующие выводы. При «разреженной» атмосфере (; = -5,5) ЛА при выполнении заданного конечного условия по углу наклона траектории не достигает заданной конечной высоты и поэтому имеет большую конечную скорость. Можно сказать, что ЛА «не хватает» подъёмной силы из-за низкой плотности атмосферы. При «плотной» атмосфере (; = 5,5) ЛА при выполнении заданного конечного условия по углу наклона траектории, напротив,

имеет конечную высоту больше заданной и поэтому имеет меньшую конечную скорость. Можно сказать, что ЛА имеет «лишнюю» подъёмную силу из-за высокой плотности атмосферы.

Поэтому при терминальном управлении следует ожидать увеличения коэффициента подъёмной силы при «разреженной» атмосфере и уменьшения коэффициента подъёмной силы при «плотной» атмосфере.

Алгоритм командного управления

Для компенсации влияния возмущений на конечные условия движения ЛА предлагается использовать алгоритм терминального управления по определению командного значения коэффициента подъёмной силы Суа2.

Для принятой программы управления (2) в общем случае имеется три параметра (Суа1, р Суа2), изменяя которые

можно компенсировать влияние действия возмущений на конечные значения трёх контролируемых фазовых координат: высоты Ик, скорости Vк и угла наклона траектории вк. Поскольку при использовании терминального управления они прогнозируются в результате интегрирования дифференциальных уравнений движения (1), то одна из фазовых координат всегда имеет заданное конечное значение, т. к. это является условием окончания интегрирования. Поэтому в регулировании нуждаются две фазовые координаты.

Примем, что определение командного управления Суа2 начинается после времени переключения /П , поскольку это позволит не увеличивать значение нормальной скоростной перегрузки пуа свыше максимально допустимой п™х из-за

возможного увеличения коэффициента подъёмной силы при терминальном управлении.

Таким образом, два параметра номинальной программы управления, а именно Суа1 и /П, остаются неизменными

при терминальном управлении. Поэтому в дальнейшем при терминальном управлении будем использовать один оставшийся параметр, а именно Суа2, при двух имеющихся краевых условиях.

Примем, что при формировании командного управления будут заданы конечное значение угла наклона траектории вк * как условие окончания интегрирования уравнений движения (1) и конечное

значение высоты Ик * как регулируемая фазовая координата.

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

На каждом шаге управления (коррекции) г длительностью Ат алгоритм формирования командного значения коэффициента подъёмной силы Суа2 предполагает следующую последовательность действий.

1. Прогнозирование движения до выполнения заданного конечного условия движения, а именно заданного значения конечного угла наклона траектории вк *.

2. Формирование командного управления: при + Ат > решается одноточечная краевая задача определения Суа2, обеспечивающего выполнение заданного конечного условия движения по высоте.

3. Прогнозирование движения с полученным управлением.

4. Принятие решения о необходимости следующей итерации для определения командного управления. Итерация проводится, если отклонение значения конечной высоты от заданного Ик * превышает установленную величину.

Одноточечная краевая задача решается методом Ньютона.

Результаты моделирования терминального управления

Для проверки работоспособности и эффективности предложенного алгоритма управления проведено моделирование возмущённого движения при идеальной навигации, ориентации и стабилизации ЛА.

Моделирование проводилось путём решения системы уравнений (1) при возмущениях плотности атмосферы согласно (3). Окончание интегрирования определялось достижением заданного конечного

значения угла наклона траектории вк *. Шаг по времени при формировании командного управления принят постоянным и равным 10 с.

На этапе проверки работоспособности предложенного алгоритма проводилась только одна итерация для определения командного значения с использова-

нием метода Ньютона, т.е. выполнялись вышеуказанные действия 1 - 3.

Результаты проведённого моделирования приведены на рис. 1 - 6.

На рис. 1 - 4 показаны отклонения регулируемой фазовой координаты - конечного значения высоты.

Рис. 1. Значение \ЛИК\ до коррекции при £ = -5,5

Рис. 2. Значение \ЛИк\ при £= -5,5 после коррекции

1800

1600

1400

1200

1000

800

600

400

200

20

350

200

150

100

50

30 40 50 60

Рис. 3. Значение \ЛИк\ до коррекции при £ = 5,5

70

I с

20 30 40 50 60

Рис. 4. Значение \ЛИк\ при £ = 5,5 после коррекции

70 и с

Суа2

1,0 -

0,9 -

0.8 -

0.7 -

0.6 -

0.5 -

0,4 -

0.3 -0.2 0.1 0.0

20

30

град. 25 -

40

50

60

20

15

10

70 С

20

30

40

50

60

70 г. С

б

Рис. 5. Командные значения Суа2 и а2 после одной итерации при £ = -5,5: а - коэффициент Суа2; б - угол а2;

ограничение по углу атаки

а

а2

Суа2 град.

20 30 40 50 60 70 С 20 30 40 50 60 70 £, с

0,0

-0.1

-0,3

-0,4

-4 -6 -8

-10 -

-12

б

а

Рис. 6. Командные значения Cya2 и а2 после одной итерации при % = 5,5: а - коэффициент Суа2; б - угол а2; ч-т-гт-г _ ограничение по углу атаки

Из рис. 1 - 4 следует, что алгоритм работоспособен, поскольку на каждом шаге управления происходит уменьшение прогнозируемого отклонения конечной высоты от заданного значения за счёт коррекции коэффициента подъёмной силы Суа2.

На рис. 5, 6 показаны командные значения коэффициента подъёмной силы Суа2 и соответствующие ему углы атаки а2.

Из рис. 5, 6 следует, как и отмечалось выше, что командные значения коэффициента подъёмной силы Суа2 увеличиваются при «разреженной» атмосфере (рис. 5, % = -5,5) и уменьшаются при «плотной» атмосфере (рис. 6, % = 5,5). Однако на некоторых шагах углы атаки а2, соответствующие полученным согласно пункту 2 алгоритма терминального управления командным значениям коэффициента подъёмной силы Суа2, выходят на ограничение, как при % = -5,5 (рис. 5, б), так и при % = 5,5 (рис. 6, б). Это объясняется тем, что прогнозируемые отклонения высоты слишком «велики» для имеющих-

ся ограниченных величинами ат1П или атах возможностей по изменению коэффициента подъёмной силы Суа2.

Поэтому на втором этапе оценки эффективности алгоритма на шаге управления проводилось несколько итераций для определения командного значения коэффициента подъёмной силы Суа2. Число итераций определялось достижением заданного значения отклонения конечной высоты, которое было принято равным 1 м.

Результаты проведённого моделирования приведены на рис. 7 - 12.

На рис. 7 - 10 показаны отклонения регулируемой фазовой координаты - конечного значения высоты. На рис. 11, 12 показаны командные значения коэффициента подъёмной силы Суа2 и соответствующие ему углы атаки а2. На рис. 13 приведено число итераций при определении командного значения коэффициента подъёмной силы Суа2 на шагах управления.

2000

20 30 4'."' 50 60 70 с

Рис. 7. Значение \ЛИк\ до первой итерации при £ = -5,5

0.7

20 30 40 50 60 70 с

Рис. 8. Значение \ЛИк\ после последней итерации при £ = -5,5

1800

20 30 40 50 60 70

Рис. 9. Значение \ЛИк\ до первой итерации при £ = 5,5

250

200

150

100

50

Суа2

1.0 -

0,5 -

0,8 -

0,7 ~

0,6 -

0,5 -

0,4 -

0,3 -

0,2 -

ОД -

О .Я -

20 30 40 50 60 70 с

Рис. 10. Значение \ЛИк\ после последней итерации при % = 5,5

а2 град. 25

20

15

10

20 30 40 50 60 70 С 20 30 40 50 60 70С

а б

Рис. 11. Командные значения Суа2 и а2 после нескольких итераций при % = -5,5: а - коэффициент Суа2; б - угол а2; - ограничение по углу атаки

Суа2

20 30 40 50 60 70 * с

0.0

-0.1

-0.2

-0,3

-0.4

Л2 град.

20 30 40 50 60 70 ¡г с

О

-2 -4 -6 -8

-ю 4 -12

б

Рис. 12. Командные значения Суа2 и а2 после нескольких итераций при % = 5,5: а - коэффициент Суа2; б - угол а2; тпп-* - ограничение по углу атаки

а

N

3 -

2 -

О -I------

20 30 40 50 60 65 Т, С

а

Рис. 13. Количество итераций при

При Е, = -5,5 происходит значительное уменьшение прогнозируемого отклонения конечной высоты (рис. 8) за счёт проведения нескольких итераций на первых четырёх шагах управления (рис. 13, а). Значение угла атаки а2 и соответствующее ему командное значение коэффициента подъёмной силы Суа2 выходят на ограничение только на последнем шаге управления (рис. 11).

При Е, = 5,5 на первых двух шагах управления значение угла атаки а2 и соответствующее ему командное значение коэффициента подъёмной силы Суа2 при первой же итерации выходят на ограничение (рис. 12). Поэтому на этих шагах управления прогнозируемое отклонение конечной высоты остаётся значительным (рис.10). На следующих шагах управления происходит значительное уменьшение прогнозируемого отклонения конечной высоты (рис. 10) за счёт проведения нескольких итераций (рис. 13, б).

Заключение

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

В случае предельно «разреженной» атмосферы при номинальной программе управления имеется недостаток высоты при избытке скорости. При терминальном управлении отклонения конечной высоты

.V

3-----

2-------

1-------

О -I------

20 30 40 50 60 65 С

б

коррекции: а - £= -5,5; б - £ = 5,5

уменьшаются за счёт увеличения командного коэффициента подъёмной силы. По сути, происходит перевод кинетической энергии ЛА в потенциальную. В конце траектории полёта из-за снижения управляемости ЛА в разреженных слоях атмосферы командное значение коэффициента подъёмной силы выходит на ограничение, определяемое максимальным углом атаки атах. Прогнозируемые отклонения конечной скорости уменьшаются до величины нескольких метров.

В случае предельно «плотной» атмосферы при номинальной программе управления имеется избыток высоты при недостатке скорости. При терминальном управлении отклонения конечной высоты уменьшаются за счёт отрицательного командного коэффициента подъёмной силы. В начале траектории командное значение коэффициента подъёмной силы выходит на ограничение из-за достаточно «жёсткого» по сравнению с максимальным (атах = 20 °) ограничения на минимальное

значение угла атаки (атт =-10°). Прогнозируемое отклонение конечной скорости снижается до -25 м/с. Дальнейшее увеличение конечной скорости ЛА невозможно, поскольку заданное конечное значение скорости Ук * определено для минимального лобового сопротивления, а при командном управлении лобовое сопротивление увеличивается.

Возможным совершенствованием ного конечного значения угла наклона

алгоритма терминального управления яв- траектории вк * или для уменьшения от-

ляется придание ему адаптиюж свойств, клонения конечного угла наклона траек-

а именно: коррекция коэффициента подъ- тории при достижении заданного конеч-

ёмной силы для уменьшения отклонения ного значения высоты И конечной высоты при достижении задан-

*

Библиографический список

1. Young D.A., Olds J.R. Responsive Access Small Cargo Affordable Launch (RASCAL) Independent Performance Evaluation // 13 th International Space Planes and Hypersonics Systems and Technologies Conference (May 2005, Capua, Italy). 23 p. http://hdl.handle.net/1853/8372. DOI:10.2514/6.2005-3241

2. Urschel P.H., Cox T.H. Launch Condition Deviations of Reusable Launch Vehicle Simulations in Exo-Atmospheric Zoom Climbs // AIAA Atmospheric Flight Mechanics Conference and Exhibit (August 2003, Austin, United States). DOI: 10.2514/6.2003-5544

3. Нечаев Ю.Н. Силовые установки гиперзвуковых и воздушно-космических летательных аппаратов. М.: Академия космонавтики им. К.Э. Циолковского, 1996. 214 с.

4. Hague C.N., Siegenthaler E., Rothman J. Enabling responsive space: F-15 microsatellite launch vehicle // Aerospace Conference Proceedings. 2003. V. 6. P. 2703-2708. DOI: 10.1109/AERO.2003.1235195

5. Балашов В.В., Бузулуков В.М., Давидсон Б.Х., Смирнов А.В. Возможности использования сверхзвуковых самолётов-носителей для запуска малых, мини- и микроспутников // Тезисы докладов XXXVIII чтений, посвящённых разработке научного наследия и развитию идей К.Э. Циолковского. Калуга: Государственный музей истории космонавтики им. К.Э. Циолковского, 2003. С. 61-63.

6. Балакин В.Л., Потапов В.И. Номинальная программа управления движением сверхзвукового самолёта-носителя // Вестник Самарского государственного аэрокосмического университета имени академика С.П. Королёва (национального исследовательского университета). 2011. № 6 (30). С. 15-21.

7. Потапов В.И. Программы управления и траектории движения сверхзвуковой первой ступени авиационно-космической системы // Вестник Самарского государственного аэрокосмического университета имени академика С.П. Королёва (национального исследовательского университета). 2010. № 1 (21). С. 75-83.

8. Школьный Е.П., Майборода А. Атмосфера и управление движением летательных аппаратов. Л.: Гидрометеоиздат, 1973. 308 с.

Информация об авторах

Балакин Виктор Леонидович, доктор технических наук, профессор, профессор кафедры автоматических систем энергетических установок, Самарский национальный исследовательский университет имени академика С.П. Королёва. E-mail: bala-kin@ssau.ru. Область научных интересов: динамика и управление движением летательных аппаратов.

Ковалев Андрей Викторович, магистрант, Самарский национальный исследовательский университет имени академика С.П. Королёва. Область научных интересов: динамика и управление движением летательных аппаратов.

TERMINAL CONTROL OF AEROSPACE SYSTEM SUBHYPERSONIC FIRST STAGE

© 2016 V. L. Balakin, A. V. Kovalyov Samara National Research University, Samara, Russian Federation

The problem of forming command control of the subhypersonic first stage of an aerospace system in climb is considered. Passive motion of spacecraft in conditions of maximum atmospheric density disturbance is analyzed. Achieving the prescribed value of the angle of climb is the terminal motion condition. Terminal height is a controlled value. An algorithm of terminal control for the formation of command value of aerodynamic lift coefficient is proposed. The Newton method with one or more iterations at the correction step is used in determining command control. The serviceability and efficiency of the algorithm compensating the influence of variations of atmospheric density on the preset terminal altitude condition of spacecraft motion are analyzed. The results of simulating spacecraft motion with terminal control for maximally «rarefied» and maximally «dense» atmosphere are discussed. Possible improvement of the terminal control algorithm is suggested.

Aerospace system, subhypersonic first stage, space vehicle, terminal control, algorithm, lift coefficient, angle of attack, atmospheric density disturbance.

References

1. Young D.A., Olds J.R. Responsive Access Small Cargo Affordable Launch (RASCAL) Independent Performance Evaluation. 13th International Space Planes and Hypersonics Systems and Technologies Conference (May 2005, Capua, Italy). 23 p. Available at: http://hdl.handle.net/1853/8372. DOI:10.2514/6.2005-3241.

2. Urschel P.H., Cox T.H. Launch Condition Deviations of Reusable Launch Vehicle Simulations in Exo-Atmospheric Zoom Climbs. AIAA Atmospheric Flight Mechanics Conference and Exhibit (August 2003, Austin, United States). DOI: 10.2514/6.2003-5544

3. Nechaev Yu.N. Silovye ustanovki giperzvukovykh i vozdushno-kosmicheskikh le-tatel'nykh apparatov [Power plants of hypersonic and space-air vehicles] Moscow: Akademi-ya Kosmonavtiki im. K.E. Tsiolkovskogo Publ., 1996. 214 p.

4. Hague C.N., Siegenthaler E., Rothman J. Enabling responsive space: F-15 microsatellite launch vehicle. Aerospace Conference Proceedings. 2003. V. 6. P. 2703-2708. DOI: 10.1109/AERO.2003.1235195

5. Balashov V.V., Buzulukov V.M., Davidson B.Kh., Smirnov A.V. Vozmozhnosti ispol'zovaniya sverkhzvukovykh samoletov-nositeley dlya zapuska malykh, mini- i mikrosputnikov. Tezisy dokladov XXXVIII chteniy, posvyashchennykh razrabotke nauchnogo naslediya i razvitiyu idey K.E. Tsiolkovskogo. Kaluga: Gosudarstvennyy muzey istorii kosmonavtiki im. K.E. Tsiolkovskogo Publ., 2003. P. 61-63.

6. Balakin V.L., Potapov V.I. Nominal motion control program in a supersonic carrier aircraft. Vestnik of the Samara State Aerospace University. 2011. No. 6

7. Potapov V.I. Control programs and motion trajectories of hypersonic first stage of an aerospace system. Vestnik of the Samara State Aerospace University. 2010. No. 1

8. Shkol'nyy E.P., Mayboroda A. Atmosfera i upravlenie dvizheniem letatel'nykh apparatov [Atmosphere and aircraft motion control]. Leningrad: Gidrometeoizdat Publ., 1973. 308 p.

Citation: Balakin V.L., Kovalyov A.V. Terminal control of aerospace system subhypersonic first stage. Vestnik of the Samara State Aerospace University. 2016. V. 15, no. 2. P. 9-21. DOI: 10.18287/2412-7329-2016-15-2-9-21

About the authors

Balakin Victor Leonidovich, Doctor of Science (Engineering), Professor, Professor of the Department of Automatic Systems of Power Plants, Samara National Research University, Samara, Russian Federation. E-mail: balakin@ssau.ru. Area of Research: flying vehicle dynamics and motion control.

Kovalyov Andrei Viktorovich, postgraduate student, Samara National Research University, Samara, Russian Federation. Area of Research: flying vehicle dynamics and motion control.