УДК 621.454.2 DOI: 10.18287/2541-7533-2019-18-4-129-145
РЕШЕНИЕ ЗАДАЧИ МИНИМАКСНОГО АДАПТИВНОГО УПРАВЛЕНИЯ РАСХОДОМ ТОПЛИВА ЖИДКОСТНОЙ РАКЕТЫ-НОСИТЕЛЯ
© 2019
А. Ф. Шориков профессор, доктор физико-математических наук;
Уральский федеральный университет, г. Екатеринбург; a.f.shorikov@urfu.ru
В. И. Калёв ведущий инженер-конструктор отдела управления движением;
НПО автоматики им. академика Н.А. Семихатова, г. Екатеринбург; persona@npoa.ru
Предлагается математическая формализация и метод решения задачи минимаксного (гарантированного) адаптивного терминального управления расходом топлива двигательной установки (ДУ) жидкостной ракеты-носителя (ЖРН). Исходная дискретно-непрерывная нелинейная модель объекта управления линеаризуется вдоль заданной опорной фазовой траектории и аппроксимируется линейной дискретной многошаговой динамической системой, включающей в себя вектор состояния (фазовый вектор) системы, вектор управления и вектор внешнего возмущения, описывающий погрешность формирования аппроксимирующей модели. При заданных геометрических ограничениях на фазовый вектор, векторы управления и возмущения для аппроксимирующей системы формулируется основная задача минимаксного адаптивного терминального управления расходом топлива ДУ ЖРН, состоящая из решения ряда вспомогательных задач минимаксного программного терминального управления. При решении каждой из них используется аппарат построения и анализа обобщённых областей достижимости аппроксимирующей линейной дискретной динамической системы, который реализуется с помощью модификации общего рекуррентного алгебраического метода. Для рассматриваемой задачи минимаксного адаптивного терминального управления расходом топлива ДУ ЖРН предлагается метод её решения и соответствующий численный алгоритм, сводящийся к реализации конечной последовательности только одношаговых алгебраических и оптимизационных операций. Эффективность предлагаемого решения исследуемой задачи демонстрируется и подтверждается на примере компьютерного моделирования управления процессом расхода топлива ДУ третьей ступени ЖРН «Союз-2-1б».
Ракета-носитель; двигательная установка; управление расходом топлива; адаптивное управление; минимаксный результат; гарантированное управление; робастное управление.
Цитирование: Шориков А.Ф., Калёв В.И. Решение задачи минимаксного адаптивного управления расходом топлива жидкостной ракеты-носителя // Вестник Самарского университета. Аэрокосмическая техника, технологии и машиностроение. 2019. Т. 18, № 4. С. 129-145. DOI: 10.18287/2541-7533-2019-18-4-129-145
Введение
Одной из основных задач, решаемых системой управления отечественных ЖРН, является задача оптимизации терминального управления расходом топлива. В качестве объекта управления в ней рассматривается совокупность ДУ и топливных баков ЖРН. Перед системой управления ставится цель: полностью и синхронно израсходовать рабочие запасы компонентов топлива (окислителя и горючего) к заданному моменту времени. Достижение этой цели может быть формализовано как минимизация некоторого критерия качества, оценивающего отклонение фазовых координат объекта управления от их желаемых значений в финальный момент времени.
Известно, например [1-3], что при моделировании объектов управления ракетнокосмической техники информация об априори неопределённых параметрах системы (начальном множестве фазовых состояний, реализации возмущающего воздействия), как правило, может задаваться в виде соответствующих ограниченных множеств возможных значений этих параметров.
В работе рассматривается исходная нелинейная дискретно-непрерывная математическая модель объекта управления, которая аппроксимируется с помощью линейной дискретной управляемой динамической системы, в которой учитывается возмущающее воздействие (возмущение), описывающее погрешность формирования аппроксимирующей модели или внешнее воздействие (помеху). Множества, ограничивающие возможные значения фазового вектора и вектора возмущения, рассматриваются как выпуклые, замкнутые и ограниченные многогранники (с конечным числом вершин) в соответствующих конечномерных векторных пространствах. Далее для краткости будем писать «многогранники», подразумевая все указанные выше свойства. Множество допустимых управляющих воздействий является конечным множеством соответствующего векторного пространства. Задача оптимизации управления расходом топлива ЖРН формулируется как задача минимаксного адаптивного (по принципу обратной связи) терминального управления, минимизирующего возможные наихудшие (максимальные) значения выбранного критерия качества, соответствующие наихудшим реализациям допустимых значений возмущения.
Основываясь на результатах работ [4;5], для исходной дискретно-непрерывной динамической системы решается задача минимаксного адаптивного терминального управления расходом топлива ЖРН в выбранном классе допустимых стратегий адаптивного управления посредством сведения её к последовательному решению вспомогательных задач минимаксного программного терминального управления сформированной линейной дискретной динамической системой. Решение каждой вспомогательной задачи предлагается осуществлять с помощью рекуррентного алгоритма, сводящего исходную многошаговую задачу к реализации конечной последовательности решения только одношаговых задач линейного и выпуклого математического программирования, конечных систем линейных алгебраических равенств и неравенств, а также к операциям над выпуклыми множествами. Для решения вспомогательных задач используется общий рекуррентный алгебраический метод построения и анализа областей достижимости (ОД) линейных дискретных управляемых динамических систем, разработанный в [4;5]. Для компьютерной реализации этого метода разработан и создан программный комплекс, описание и возможности которого представлены в [4;5]. Модификации алгоритмов численной реализации этого метода представлены в [6].
Математическая модель расхода топлива ДУ ЖРН
На промежутке времени \j30,df ^ рассматривается нелинейная дискретнонепрерывная математическая модель [7], описывающая динамику работы ДУ ЖРН. Полагается, что в0 - время начала управляемого режима работы ДУ, 0f - время выключения ДУ. Управляющее воздействие (управление) в данной модели из-за особенностей конструкции тракта управления расходом топлива ДУ ЖРН реализуется в заданные дискретные моменты времени {90,в1, ...,вт ч] <z.^90,9f ], где 6T = 9f, в соответствие которым может быть поставлен целочисленный набор (целочисленный промежуток времени) 0,T-1 = {0,1,...,T-1] {T е N). Здесь T - количество заданных моментов времени, в которые планируется осуществлять реализацию управления; N - множество всех натуральных чисел [8]. За счёт управления u (t) на целочисленном промежутке времени 0, T -1 регулируется угол поворота дросселя ath (t) в соответствии с рекуррентным соотношением:
ah (t + 1) = а* (t) + cou (t), tе О, Г -1, ath (0) = 0, (1)
где c0 - известный коэффициент привода дросселя. Важно отметить, что значение угла поворота дросселя ath (в) при в е\вг ,9t+1), t е 0, T -1 имеет фиксированное значение (не изменяется).
Поворот дросселя на угол ath (в) изменяет значение коэффициента соотношения расходов окислителя и горючего [2] следующим образом:
Кт (в) = K + &K + c,а.,(в, ве[в„,в/], (2)
где К - номинальное значение коэффициента соотношения расходов компонентов топлива; ДК - известный параметр выставки дросселя в номинальное положение; с, - коэффициент эффективности дросселя.
Выражение для вычисления пустотной тяги ДУ ЖРН имеет вид
P, (в) = P + с2 (Km(в)-К)2 + Сз (Km (в)-К), ве 00], (3)
где P - номинальное значение пустотной тяги ДУ; с2, с3 - коэффициенты рабочего режима ДУ ЖРН.
Удельный импульс тяги ДУ РН в пустоте вычисляется согласно формуле
I, (в) = I +с4 ( (в) - к) + cs (Кт (в) - К), ве[в„,в, ], (4)
где I - номинальное значение удельного импульса в пустоте; с4, с5 - коэффициенты рабочего режима ДУ ЖРН.
Значения массовых расходов окислителя и горючего вычисляются посредством следующих соотношений:
то (в)
mf (в)
P (в) Кт (в)
L (в ( + Кт (в) ,
P (в
Is (в ( + Кт (в) ,
ве[в0,вf ] , ве\в0,в/ ].
(5)
Текущие массы компонентов топлива в баках, рассчитанные на выполнение полётной программы ЖРН, вычисляются по формулам:
М„(в) = Mо.о +ДМо -вто(0)d0,
•°0
Mf (в) = М + ДМf -f mf (в) dO,
•°0
где M0o,M0f - соответственно номинальные рабочие массы окислителя и горючего,
заправляемые в топливные баки для реализации программного режима тяги; AMo, AMf - погрешности заправки компонентов топлива в баки.
Расчётная опорная траектория движения ЖРН реализуется при известном режиме изменения тяги ДУ, которому соответствует заданная программа изменения рабочих масс компонентов топлива в баках:
M'f(e) = M„о -вт;', m'f = ~~~K,
I + IK
M'f (0) = M„,-em'f, mf =
(7)
В связи с порогово-дискретным принципом измерения уровней компонентов топлива в баках [1-3; 8] и, как следствие, управлением, формируемым в эти дискретные моменты времени |в0,в1,^,вг-х] в бортовой цифровой вычислительной машине, линеаризованная вдоль опорной траектории (7) непрерывная составляющая исходной модели дискретизируется. Процедура линеаризации и дискретизации исходной модели динамики расхода топлива ДУ ЖРН описана в работе [7].
Рассмотрим на целочисленном промежутке времени 0,T = {0,1,...,T}(T е N) линейную дискретную динамическую систему, аппроксимирующую исходную математическую модель. Следующие рекуррентные соотношения позволяют рассчитывать значения массовых расходов компонентов топлива в ДУ ЖРН:
то (t +1) = то (t) + au (t) + ylvl (t), mo (0) = mrf + aAK,
mf (t +1) = mf (t) + /u (t) + y2v2 (t), mf (0) = mrf + /AK, ^
где t е 0, T -1; a, /3 - коэффициенты, сформированные в результате разложения в ряд Тейлора уравнений исходной модели; u (t) - скалярное управление; v (t), v2 (t) - неконтролируемые возмущения, вызванные погрешностью формирования модели (аналогичным образом могут рассматриваться любые другие возникающие в данной системе возмущающие воздействия, например, изменение температуры компонентов топлива в полёте, нештатные ситуации в работе ДУ); ух, у2 - известные коэффициенты, оцениваемые на этапе моделирования.
Рекуррентные уравнения для вычисления значений масс окислителя и горючего в топливных баках ДУ ЖРН имеют вид:
Mo (t +1) = Mo (t)-AT(t)mo (t), M„ (0) =M,o +AM.,
Mt(t +1) = Mt(t)-AT(t)mf (t), Mf (0) = M0f + AMf, ( )
где AT (t) = в(+j -в(, t е 0, T -1 - шаг дискретизации исходной системы (время между двумя соседними управлениями).
Тогда динамика объекта управления (8), (9) описывается дискретной многошаговой управляемой системой в векторно-матричном виде
x (t +1) = Ax (t) + Bu (t) + Dv(t), x (0) = x0. (10)
Здесь х (t)е М4 - фазовый вектор системы (здесь и далее, для n е N, R” - n-мерное
векторное пространство векторов-столбцов), на который наложено фазовое ограничение по безотказной работе ДУ вида:
х(t)е X1 (t)е М4, t е 0, T, (11)
где Xj (t) - многогранник в векторном пространстве М4; u (t) е М1 - управление, стеснённое заданным ограничением:
u (t)е U1 (t)е М1, t е 0, T -1, (12)
где U (t) - конечный набор векторов в М1, определяющий все возможные значения управления в момент времени t ;
v(t) е М2 - вектор неконтролируемых возмущений (помех), удовлетворяющий известному ограничению
v(t)е V1 (t)е М2, t е 0, T-1, (13)
где Vj (t) - многогранник в векторном пространстве М2;
A е М4x4, B е М4х1, D е М4х2 - матрицы состояния системы, управления и возмущения соответственно; х (0) = х0 - заданное начальное значение фазового вектора.
Для фиксированного целочисленного промежутка времени г, T ^ 0, T (г< T) и ограничения (12) определим множество всех допустимых реализаций программных управлений u () = {u (t )} в виде следующего конечного множества:
U(г,Т) = ju ( ) | u ( ) = ju(t)}te— е М1х(г-г), Vt е г, T -1, u (t)е U1 (t)J. (14)
Аналогично для целочисленного промежутка времени г, T и ограничения (13) определим множество всех допустимых программных реализаций вектора возмущений v () = {v ( )}„ 1 соотношением
V(r,Т) = {v()1 v() = {v(t)}te— е М2(T-г)Vt е г,Т-1,v(t) е V1 {t)}. (15)
Набор w(г) = {г, х(г)} е 0, Tх Х1 (г) (w(0) = w0 = {0,х0}) назовём т-позицией дискретной динамической системы (10) - (13), а множество
W (г) = {г} х Xj (г) (W (0) = {w (0) = w0 | w0 = {0, х0}}) - множеством всех допустимых т-позиций.
Качество рассматриваемого процесса управления на промежутке времени г, T ^ 0, T будем оценивать выпуклым терминальным функционалом
у— : W (r)x U (г, T)x V (г, T) ^ К1, значения которого на реализациях наборов (w(),u(),v( )) G W(r)x U(г, T )x V (г, T) определяются следующим образом:
у-(((т),и()v(-)) = \\х(T)-xd\l =Ф(x(T)) . (16)
Здесь и( ) = {u(t)}tG—,v( ) = {v(t)}tG—; x(T) = x(t;г,T,x0,и( ),v(•)) - финальное фазовое состояние движения (фазовой траектории) динамической системы (10); xd g К4 -заданный вектор, определяющий желаемое финальное фазовое состояние рассматриваемой системы; || • ||4 — евклидова норма в пространстве К4; Ф:Ж4 ^Ж1 — выпуклый функционал, определённый на финальных фазовых состояниях объекта управления.
Постановка задачи минимаксного адаптивного управления
Сформулируем вспомогательную задачу минимаксного программного терминального управления [4] аппроксимирующей системой (10) - (13), (16).
Задача 1. Для заданного целочисленного промежутка времени г, T ^ 0, T (г < T) и
реализации т-позиции w(г) = {г,x(г)} g W(г) в управляемой динамической системе
(10) - (13), (16) требуется найти множество U(e) (г, T, w(r))e U(г, T) программных
управлений и(e) ( ) = {u(e) (t)} _g U (г , T), удовлетворяющих условию минимакса [4]:
U(e) (г,T,w (r)) = {u(e) () 1 u(e} () = {u(e) (t)}tG____ G U (r, T),
)(г,T,w(г) ) = г (w(,u(e)(),v(e)()) = max г (w(,uWQ,v() ) = (17)
= min max v— (w(г),u(•),v(•))},
u (•)GU(r,r ('(•)GV(r,T )
которое будем называть множеством минимаксных программных управлений, а число c(ye) (г, T, w (г)) будем называть минимаксным или оптимальным гарантированным результатом для этой задачи.
Сформулируем содержательно задачу минимаксного адаптивного терминального управления рассматриваемой динамической системой (10) - (13), (16).
На промежутке времени 0, T необходимо сформировать управление u ()={u (t )g U1(t )L системой (10) с ограничениями (11) - (13) как стратегию адаптивного (по принципу обратной связи) терминального управления, использующую полную информацию о её t-позиции w(t) = {t, x(t)} g W (t) в каждый момент времени
t g 0, T -1 так, чтобы к концу процесса управления функционал , определяемый согласно (16), принимал наименьшее возможное значение при условии, что значения возмущения v( ) = {v(t) __ g V (0, T) могут быть наихудшими, то есть максимизирующими данный функционал.
Для формализации приведённой цели адаптивного управления введём ряд определений. Допустимой стратегией адаптивного терминального управления расходом
топлива ДУ ЖРН Ua для дискретной динамической системы (10) - (13), (16) на промежутке времени 0,T будем называть отображение Ua : W(г)^ U (г), которое каждому моменту времени ге 0,T-1 и реализации т-позиции w(г) = {г,х(г^е W (г) назначает
множество Ua ( w (г ))=U1 ( г) управлений и(г)е Ux (г) . Множество всех допустимых стратегий адаптивного терминального управления расходом топлива ДУ ЖРН обозначим как ua.
Стратегией минимаксного адаптивного терминального управления системой (10) - (13), (16) на промежутке времени 0, T будем называть такую реализацию специфической стратегии U^ = U^^ (w (г)) е ua, w(г) = |г,х(г)| е W(г), г е 0,T-1, из класса всех допустимых стратегий адаптивного терминального управления расходом топлива ДУ ЖРН Ua, которая определяется следующими соотношениями:
т-позиций w(г) = |г, х(е) (г) j е W (г)
1)для всех ге 0,T-1 и
( (0) = {0, Х0 j = w0) пусть
(w(e) (г)) = u(e) (w(e) (г)) , (18)
где множество U(e)( w(e) (г)) определяется согласно (17) из решения Задачи 1, а именно:
(т)) = {u(e) (г) | и(e) (г) е Uj (г),и(e) (•) = {u(t)j ____е U(e) (г,T, w(e) (r))j;
для допустимых на промежутке 0,г реализаций управления
ur (-) = k(tе U (0,г) и возмущения VrO = {Vr(tе V (0,г) :
х(е) (г) = х (г;0,г,х0,иг (•), v (•));
2)для всех ге 0,T-1 и т-позиций w(г) £ W(г) пусть
Ui''( w (г)) = U, (г), (19)
и пусть u(e) (•) = {uae) (r)j _е U(0,T) и V.e| (') = {v^1 Dj е V(,T) есть реализа-
ции соответственно управления и возмущения на промежутке времени 0,T, которые сформированы в результате использования стратегии U(ae) е Ua на этом промежутке времени, причём значения u^ (T -1) и v^ (T -1) удовлетворяют соотношению (17) при г = T -1 и можно вычислить следующее число:
Уё )=r*J ( «a)(■), (w(e)(T - ^ Uie)(T - ^ Va}(T - !)) =
-1) Ут-гг (((т - (, ((т -1, v (T -!)) = ф (xW (т)) = (20)
((;0,т, (оА^ (-),va (•))) = А (0,т, wo )
= min
u(T-l}eUr (T-1) v(T
= Ф
где w(e) (T -1) = {T -1, xW (T - l)}e W (T -1) , x(e) (T -1) = x (t -1; 0, T, (,, A (•), Va (•)) .
Сформулируем следующую многошаговую задачу минимаксного адаптивного терминального управления расходом топлива ДУ ЖРН для аппроксимирующей дискретной управляемой динамической системы (10) - (13), (16).
Задача 2. Для заданного промежутка времени 0, T и начальной позиции w0 = {0, x0}e W (0) в дискретной динамической системе (10) - (13), (16) требуется найти стратегию минимаксного адаптивного терминального управления U(ae) = U(ae) (w (t)) e U*a , W (t) = {t, X (t)} e W (t) , t e 0, T-1, которая удовлетворяет соотношениям (18), (19), и число АУ( 0, T, w0) - минимаксный (оптимальный гарантированный) результат, соответствующий реализации этой стратегии, путём реализации конечного числа только одношаговых операций.
Нетрудно убедиться, что оптимальный гарантированный результат, полученный в результате решения Задачи 2, будет, по крайней мере, не хуже оптимального гарантированного результата, обеспечиваемого посредством решения Задачи 1. Отметим, что хотя Задача 2 является основной, но её решение и формализация базируются на решении Задачи 1.
Алгоритмы решения Задачи 1 и Задачи 2
При решении сформулированной Задачи 1 используются результаты работ [4-7], а также аппарат построения и анализа ОД линейных управляемых динамических систем, широко применяемый в теоретических и прикладных задачах [4-6; 9-11].
Введём необходимое определение обобщённой ОД для системы (10) - (13) [4].
Определение 1. Обобщённой ОД линейной дискретной управляемой динамической системы (10) - (13) при фиксированном допустимом программном управлении
u* () = {u '(t)}te_ e U( t, T) на момент времени T, соответствующей набору
(X(t) , u* (•)), где X (0) = {x0} , называется множество вида:
G (X (t), u* О;т Иx (T_)KT )e Г’x ('+1) = ^(') x (') + *(')u'(t) + D (t )v (t ),(21)
t e t,T-1, x(r) e X(r), v(t) e V1 (t)}.
В силу принятых предположений о том, что множества в ограничении (13) относятся к классу многогранников, в [4; 5] было показано, что множество, описывающее обобщённую ОД такой динамической системы, также будет принадлежать к этому классу.
Приведём описание общего рекуррентного алгебраического метода [4; 5] для построения обобщённых ОД рассматриваемой аппроксимирующей системы.
Алгоритм 1. (Общий рекуррентный алгебраический метод построения ОД).
1
. Инициализация: X (г) = {х (г)] , и * () = {u* (t ) __ е U (г, T ) .
2. Для всех t ег, T -1 последовательно выполнить
V (t) = Vertex ( V1 (t));
X (t +1) = ( A (t) X (t)® B (t) u* (t ))0 D (t )V (t);
X (t +1) = RemoveRedundancy (X (t +1)).
3. Закончить.
Операции Vertex и RemoveRedundancy в алгоритме обозначают соответственно выделение вершин приведённого внутри скобок многогранника и решение задачи нахождения крайних точек для приведённого внутри скобок множества точек (формирование неизбыточного вершинного описания многогранника), которая решается как задача линейного математического программирования способом, предложенным в [6]. Операция суммирования множеств в данном алгоритме понимается как геометрическая сумма этих множеств и обозначена знаком ® .
Учитывая сказанное выше и результаты работ [4-6], алгоритм решения Задачи 1 минимаксного программного терминального управления расходом топлива ДУ ЖРН можно свести к следующей алгоритмической последовательности.
Алгоритм 2. (Минимаксное программное терминальное управление).
1. Упорядочение по возрастанию натурального индекса j конечного множества
U(_ T), состоящего из N допустимых программных управлений
и(1) () = {u(j)(t)j е U(_T) (формирование множества U(_T) = {u(j)(t)] ^ ).
2. Построение обобщённых ОД G_,X(г),и(j)();T) при фиксированных допустимых управлениях и( 1 () для всех j е 1, N (Алгоритм 1).
3. Формирование подмножества U _, T ) е U _, T) разрешённых программных
управлений U()() = {u()(t)j е U(г,T), i е 1,M посредством удаления из множества
U(г,T) допустимых программных управлений u(j) () = {u(j) (t)j е U(_T), обобщённые ОД которых нарушают ограничения (11) для любого t ег +1,T, и упорядочение по возрастанию натурального индекса i множества U (г, T) .
4. Формирование фасетного описания (в виде систем линейных алгебраических равенств и неравенств) многогранников G(г,{х(г)],u()();T), i е1,M , например, с помощью модификации метода двойного описания [12], предложенной в [4].
5. Оптимизация функционала (16) на множестве G(г,{х(r)j,U()();T) для всех
i е 1, M методами выпуклого математического программирования, то есть нахождение следующего значения функционала:
г = (w (( ЧУ v(')()
JX (T )-
= max_
v(')Gv(r,r'
г (w (г), u ()() v ()):
max
XA\ = X'
х(Т )ес(
(T )-
d 114 -
где х (Т ) = х (т ;т, Т, x (г), u^Q, v (•)); x(e)(T ) = x (t ;t, T, x (r), uwQ, v(e) (•)) (решается, например, с помощью метода Зойтендейка [13]).
6. Нахождение множества U( % Т, w (г)) минимаксных программных терми-
нальных управлений и гарантированного результата управления е/е)(т, Т, w (г)) по средством решения следующей конечной задачи дискретной оптимизации:
U(e) (т, Т, w (г)) = {u(e) (•) | u(e) (•) е U (т, Т),
„mUfr) «Шт/т(w ()u v ())=й? т=т}.
(23)
Используя Алгоритм 2, опишем решение основной Задачи 2 также в виде алгоритма. Пусть в дискретной динамической системе (10) - (13), (16) на промежутке времени 0, т реализовались управление ur(-) = {ur(t)}^_1 е U (0, т) и возмущение
vr (•) = {vr (t )}te_- е V (о, г) . Тогда для те 0, Т -1 минимаксная адаптивная терминаль-
ная стратегия U(ae) е U формируется путём реализации следующей конечной последовательности операций.
Алгоритм 3. (Минимаксное адаптивное терминальное управление).
1. Для допустимых на промежутке 0, г реализаций управления
urO = {ur(t)},ет е U (0,т) и возмущения vrO = {vr(tе V (0,т) формирование w(e) (г) = {г, x(e) (г)} е W (г) ( w(e) (0) = {0, Х0} = W0),
'(r) = x (t;0,t,x0,utQ,vtO).
т-позиции w (t) = 'Jt,x (г)е W(г) (w (0) = {0,x0 } = w0), где
x( v
2. Решение Задачи 1, то есть вычисление множества минимаксных программных терминальных управлений U(e)(r, Т, w(e) (г)) вида (22) и гарантированного результата
е{ге)(т, Т, w(e) (г)) = У~Т на промежутке времени г, Т (Алгоритм 2).
3. Вычисление множества U(ae) (w(e) (г ))=U й w|*)(r)e u, (г) по формуле (18).
4. Выбор любого управления u^ (г) е U(ae) (w(e) (г)) .
5. Вычисление по уравнениям нелинейной системы (1), (2) (т+1)-позиции системы
ww(r + 1) = {г +1, xw(r +1)}, xw(r +1) = x (г + 1;т,т +1, x0, u[(;)(z), va (г)),
где va (г) е V (г) есть допустимая реализация возмущения в момент времени т (либо
формирование её любым другим доступным способом). Например, если процедура реализуется на реальной ЖРН, то информация о (т+1)-позиции может быть получена путём обработки информации с датчиков уровней и расходомерных датчиков окислителя и горючего.
6. Осуществление перехода к пункту 1, если т< Т — 2, в противном случае - к пункту 7.
7. Вычисление согласно (20) оптимального гарантированного (минимаксного) результата решения Задачи 2 ее (0, Т, w0), соответствующего реализациям управления
u\
0 = {u!T)(t)[е0Т—j е U (0 Т) и возмущения Va () = {« (t)},е0Т— е V (0, Т), которые яв-
ляются результатом использования стратегии минимаксного адаптивного терминально-
Отметим, что предложенный алгоритм решения Задачи 2 - минимаксного адаптивного терминального управления расходом топлива ДУ ЖРН - сводится к реализации только конечного числа операций - решения задач: конечной дискретной оптимизации, конечного числа систем линейных алгебраических равенств и неравенств, операций над выпуклыми многогранниками.
Предлагаемый метод решения Задачи 1 и Задачи 2 апробирован на модельном примере, иллюстрирующем решение задачи минимаксного адаптивного терминального управления расходом топлива ДУ третьей ступени ЖРН «Союз-2-1б».
Работа ДУ ЖРН рассматривается на отрезке времени \О0,в1 ], где в0 - время начала установившегося режима функционирования ДУ, вх - время выключения ДУ. В топливные баки окислителя и горючего заправлено нештатное количество топлива -Mo.0 + AMo и M0f +AMf соответственно. Турбонасосный агрегат и дроссельная заслонка изначально выставлены на отклонённый от номинального коэффициент соотношения массовых расходов компонентов топлива K + AK. Вал привода дросселя в каждый момент времени может занимать одно из девяти фиксированных положений, составляющих конечное множество в ограничении управляющего воздействия (12). Требуется к моменту времени вх (терминальное управление) одновременно израсходовать окислитель и горючее в топливных баках. При этом для обеспечения безопасной работы ДУ запрещается отклонять коэффициент соотношения массовых расходов компонентов топлива больше чем на 7% от номинального значения K (ограничение (11)).
Для решения этой задачи используется исходная нелинейная дискретнонепрерывная математическая модель (1) - (6), достаточно точно описывающая динамику установившегося режима работы ДУ. Значения параметров данной модели приведены в табл. 1.
Таблица 1. Значения параметров нелинейной системы
I, с P, тс K ЬЙ о о M 0 f , кг AMo, кг AM f , кг AK
360 30.24 2.5 18000 7200 100 -50 -0.075
C0 c1 C2 C3 C4 C5 ^, с вх, с
10 0.01 -6700 900 -15 -12 0 300
В соответствие исходной нелинейной дискретно-непрерывной модели ставится новая модель, которая формируется путём линеаризации относительно опорной траектории (7) и дискретизации с шагом, равным заданному числу управляющих воздействий (T = 4). При этом предполагается, что моменты времени переключения управляющего воздействия в системе (1) - (6) совпадают с целочисленными значениями
промежутка времени 0, T -1 = 0,3. Тогда сформированная аппроксимирующая модель (10) имеет следующий вид:
го управления
на промежутке времени 0, T, где
Модельный пример
х(t +1) = Ax(t) + Bu(t) + Dv(t), t e 0,3.
Здесь x(t) = {x (t), x2 (t), x3 (t), x4 (t)}e E4; x1 (t) - массовый расход окислителя; x2 (t) - масса окислителя в баке; x3 (t) - массовый расход горючего; x4 (t) - масса горючего в баке; фазовый вектор системы с учётом условия безотказной работы ДУ (11) для всех t = 0,4 ограничен многогранником вида
X1 (t) = {57.5 < x1 (t) < 62.5, - 30 < x2 (t) < 20000, 23 < x3 (t) < 25, - 30 < x4 () < 8000}; u (t) - скалярное управляющее воздействие, обозначающее фиксированные положения вала привода дроссельной заслонки, принимающее свои значения из конечного множества U (t) = {1; 0.75; 0.5; 0.25; 0; - 0.25; - 0.5; -0.75;-1}, Vt 0,3; v(t) - погрешность формирования аппроксимирующей модели, принимающая свои значения из множества: V (t) = {v(t) | - 0.06 < v1 (t) < 0.06, - 0.06 < v2 (t) < 0.06} e E2, Vt e 0,3; начальное фазовое
состояние x(0) = x0 = (59.134; 18100; 24.385; 7150)T (расчётная изначальная несинхронность ^кон = 12.9 с); матрицы A, B и D принимают значения:
Г 1 0 0 0 > Г 1.064 > Г1 0 Л
-75 1 0 0 0 0 0
A = 0 0 1 0 , в = -0.534 , D = 0 -1
v 0 0 -75 17 v 0 7 v 0 0 7
Качество процесса управления оценивается значениями выпуклого терминального функционала в финальный момент времени T = 4 :
ф((T)) = V(x, (T)-60)2 +((T)) +( (T)-24)2 +(x4 (T)) +( (T)-Kx, (T))2,
в котором первые четыре слагаемые под квадратным корнем обозначают отклонение от номинального (желаемого) финального состояния xd = (60; 0; 24; 0)T, а пятое слагаемое
отвечает за «одновременность» окончания компонентов топлива.
Сначала рассмотрим решение задачи в постановке минимаксного программного терминального управления для упрощённой системы (Задача 1), которое реализовано с помощью Алгоритма 2. В качестве решения сформировано множество минимаксных
программных терминальных управлений U( e)(0, T, w (0)) , состоящее из одного допустимого программного управления u(e)( ) = {u(e) (t)} _={1,1, - 0.75, -0.25}, с оптимальным гарантированным результатом решения Задачи 1 не хуже, чем (0,T, w0) = 115.79. При использовании найденного управления u(e) (•) в исходной модели финальное фазовое состояние принимает значение
x(6) = (60.26; 22.6; 23.86; 5)T, в котором значение показателя качества равно ф(x(61)) = 25.295. Таким образом, если использовать минимаксное программное управление u(e) (), в баках окислителя и горючего останется по 22.6 кг и 5 кг соответственно, а несинхронность окончания компонентов топлива составит ^tKm = 0.17 с.
Далее рассмотрим решение задачи в постановке адаптивного минимаксного терминального управления исходной системой (Задача 2), реализованное с помощью Алгоритма 3. Результатом применения минимаксной адаптивной стратегии U^ стало
множество управлений U^ (t)), состоящее из допустимого на промежутке време-
ни 0,4 управления u(e) (•) = {u^ (w (t)) _= {1, 0.5, 0.25, -1}, гарантирующего результат решения Задачи 2 не хуже, чем (0, T, w0 ) = 22.658 . При подстановке найденного управления (•) в исходную модель финальное фазовое состояние принимает значе-
ние x(в1) = (60; 20.3; 24; 4.5)T, в котором показатель качества равен ф(x (в )) = 22.657.
Таким образом, при использовании стратегии минимаксного адаптивного управления в баках окислителя и горючего остаётся по 20.3 и 4.5 кг соответственно, а несинхронность составит Кон = 0.15 с.
На рис. 1 представлены проекции фазовых траекторий исходной модели, порождённые найденным минимаксным программным терминальным управлением u(e) ( ) и
стратегией минимаксного адаптивного терминального управления U^). В табл. 2 сведены результаты применения найденных минимаксных программного и адаптивного управлений к исходной модели (1) - (6), а также их результаты для наихудшего случая в упрощённой модели. На рис. 2 показаны графики изменения параметра несинхронности Atкон при решении Задачи 1 и Задачи 2.
Из результатов моделирования можно сделать вывод, что стратегия минимаксного адаптивного терминального управления при её использовании для управления процессом расхода топлива ДУ третьей ступени ЖРН «Союз-2-1б» обеспечивает оптимальный гарантированный (минимаксный) результат не хуже, чем результат управления в соответствующей линейной дискретной системе при возможных наихудших реализациях возмущений (погрешностей аппроксимации модели).
Рис. 1. Фазовые траектории (слева - бак окислителя, справа - горючего) при решении Задачи 1 и Задачи 2: штриховая линия - опорная траектория; штрихпунктирная линия - траектория исходной системы при минимаксном программном терминальном управлении; сплошная линия - траектория исходной системы при минимаксном адаптивном терминальном управлении
Таблица 2. Результаты моделирования
Управление * (Т ), кг/c *2 (Т ) , кг *3 (Т), кг/c *4 (Т) , кг ф (* (т )) Кон
Опорная траектория О II > 60 0 24 0 0 0
Упрощённая 60.08 47.7 24.09 -22.2 115.79 1.71
модель Mie)0 59.82 20.3 23.91 4.5 22.658 0.15
Исходная »ио 60.26 22.6 23.86 5 25.295 0.17
модель '£’(■) 60 20.3 24 4.5 22.657 0.15
О 50 100 150 200 250 300
С
Рис. 2. Изменение параметра несинхронности А(кон: штриховая линия — минимаксное программное управление; сплошная линия — минимаксное адаптивное управление
Заключение
Исследована нелинейная дискретно-непрерывная управляемая система, описывающая динамику расхода топлива ДУ ЖРН, которая путём линеаризации относительно опорной траектории и дискретизации преобразована в соответствующую ей линейную дискретную управляемую динамическую систему. Задача оптимизации управления расходом топлива ДУ ЖРН сформулирована как задача минимаксного адаптивного терминального управления линейной дискретной динамической системой с выпуклым функционалом качества.
Сформулированная многошаговая задача минимаксного адаптивного терминального управления расходом топлива ДУ ЖРН решается путём реализации конечной последовательности решений только одношаговых задач линейного и выпуклого математического программирования, дискретной оптимизации, конечных систем линейных алгебраических равенств и неравенств, а также выполнения операций над выпуклыми множествами. На основе описанного метода решения рассматриваемой задачи разработан соответствующий численный алгоритм и получены результаты моделирования для задачи терминального управления расходом топлива ДУ третьей ступени ЖРН «Союз-2-1б».
Результаты моделирования демонстрируют как эффективность общего рекуррентного алгебраического метода построения ОД линейных дискретных управляемых динамических систем, используемого при решении вспомогательной задачи минимаксного программного терминального управления, так и предлагаемый численный алгоритм для решения задачи минимаксного адаптивного терминального управления. Результаты позволяют сделать вывод о возможности применения предлагаемой методики при разработке систем управления расходом топлива ДУ ЖРН.
Библиографический список
1. Петров Б.Н. Избранные труды. Т. 2. Управление авиационными и космическими аппаратами. М.: Наука, 1983. 328 с.
2. Челомей В.Н. Пневмогидравлические системы двигательных установок с жидкостными ракетными двигателями. М.: Машиностроение, 1978. 289 с.
3. Завадский В.К., Иванов В.П., Каблова Е.Б., Кленовая Л.Г. Системы управления расходованием топлива (назначение, принципы построения, алгоритмы управления) // Сборник трудов Всероссийской конференции «Актуальные проблемы ракетнокосмической техники (V Козловские чтения)» (11-15 сентября 2017 г., Самара). Т. 1. Самара: СамНЦ РАН, 2017. С. 243-254.
4. Шориков А.Ф. Минимаксное оценивание и управление в дискретных динамических системах. Екатеринбург: Уральский университет, 1997. 242 с.
5. Тюлюкин В.А., Шориков А.Ф. Об одном алгоритме построения области
достижимости линейной управляемой системы // В сб.: «Негладкие задачи
оптимизации и управление». Свердловск: УрО АН СССР, 1988. С. 55-61.
6. Шориков А.Ф., Булаев В.В., Горанов А.Ю., Калёв В.И. Аппроксимация областей достижимости нелинейных дискретных управляемых динамических систем // Вестник Бурятского государственного университета. Математика, информатика. 2018. № 1. С. 52-65. DOI: 10.18101/2304-5728-2018-1-52-65
7. Шориков А.Ф., Калёв В.И. Формирование линейной дискретной динамической модели для решения задачи оптимального терминального управления расходом топлива ракеты-носителя // Труды Пятой Международной научной конференции «Информационные технологии и системы» (24-28 февраля 2016 г., Банное). Челябинск: Челябинский государственный университет, 2016. С. 61-66.
8. Бобровников Г.Н., Катков А.Г. Методы измерения уровня. М.: Машиностроение, 1977. 168 с.
9. Брайсон А., Хо Ю-Ши. Прикладная теория оптимального управления. М.: Мир, 1972. 544 с.
10. Красовский Н.Н. Теория управления движением. М.: Наука, 1968. 476 с.
11. Красовский Н.Н., Субботин А.И. Позиционные дифференциальные игры. М.: Наука, 1974. 456 с.
12. Черников С.Н. Линейные неравенства. М.: Наука, 1968. 488 с.
13. Зойтендейк Г. Методы возможных направлений. М.: Изд-во иностранной литературы, 1963. 176 с.
SOLVING THE PROBLEM OF MINIMAX CLOSED-LOOP CONTROL
OF LIQUID-PROPELLANT LAUNCH VEHICLE FUEL CONSUMPTION CONTROL
© 2019
A. F. Shorikov Professor, Doctor of Science (Phys. & Math.);
Ural Federal University, Yekaterinburg, Russian Federation; a.f.shorikov@urfu.ru
V. I. Kalev Leading Engineer of the Department of Motion Control;
Scientific and Production Association of Automatics named after Academician
N.A. Semikhatov, Yekaterinburg, Russian Federation;
persona@npoa.ru
The paper provides mathematical formalization and a method of solving the problem of minimax (guaranteed) closed-loop terminal control of fuel consumption of a liquid-propellant launch vehicle power plant. The initial discrete-continuous nonlinear model of the controlled object is linearized along the given reference phase path and is approximated by a linear discrete-time multistep dynamical system. The approximated system includes the state vector, the control vector and the disturbance vector that defines the error of formation of the approximated model. Taking into account the geometrical constrains of control and disturbance vectors in the approximated system, we formulate the main problem of minimax closed-loop terminal control of propellant consumption of the launch vehicle’s propulsion system. This problem consists in solving a number of auxiliary tasks of minimax open-loop terminal control. To solve each of these tasks we use an instrument of development and analysis of generalized attainability domains of the approximated linear discrete dynamical system. These techniques are implemented by modifying the general recurrent algebraic method. To solve the problems under consideration we propose an approach and an appropriate numerical algorithm that is reduced to the implementation of a finite sequence of only one-step algebraic and optimization operations. The efficiency of the proposed approach to solving the problem under consideration is demonstrated and verified by a computer simulation example. This simulation example consists in controlling the process of propellant consumption for “Soyuz-2.1b” launch vehicle’s third stage propulsion system.
Launch vehicle; propulsion system; closed-loop control; propellant consumption control; minimax; guaranteed control; robust control.
Citation: Shorikov A.F., Kalev V.I. Solving the problem of minimax closed-loop control of liquid-propellant launch vehicle fuel consumption control. Vestnik of Samara University. Aerospace and Mechanical Engineering. 2019. V. 18, no. 4. P. 129-145. DOI: 10.18287/2541-7533-2019-18-4-129-145
References
1. Petrov B.N. Izbrannye trudy. T. 2. Upravlenie aviatsionnymi i kosmicheskimi apparatami [Selectas. Vol. 2. Air- and spacecraft control]. Moscow: Nauka Publ., 1983. 328 p.
2. Chelomey V.N. Pnevmogidravlicheskie sistemy dvigatel'nykh ustanovok s zhidkostnymi raketnymi dvigatelyami [Pneumatic / hydraulic systems of liquid-propellant engine power plants]. Moscow: Mashinostroenie Publ., 1978. 289 p.
3. Zavadskiy V.K., Ivanov V.P., Kablova E.B., Klenovaya L.G. Propellant-consumption control systems (mission, design concept, algorithms for systems). Sbornik trudov Vserossiyskoy konferentsii «Aktual'nye problemy raketno-kosmicheskoy tekhniki (V Kozlovskie chteniya)» (September, 11-15, 2017, Samara). Samara: Samarskiy Nauchnyy Tsentr RAN Publ., 2017. P. 243-254. (In Russ.)
4. Shorikov A.F. Minimaksnoe otsenivanie i upravlenie v diskretnykh dinamicheskikh sistemakh [Minimax estimation and control in discrete dynamical systems]. Ekaterinburg: Ural University Publ., 1997. 242 p.
5. Tyulyukin V.A., Shorikov A.F. Ob odnom algoritme postroeniya oblasti dostizhimosti lineynoy upravlyaemoy sistemy. V sb.: «Negladkie Zadachi Optimizatsii i Upravlenie». Sverdlovsk: UrO AN SSSR Publ., 1988. P. 55-61. (In Russ.)
6. Shorikov A.F., Bulaev V.V., Goranov A.Yu., Kalev V.I. Approximation of attainability domains of nonlinear discrete-time controlled dynamical systems. BSU Bulletin. Mathematics, Informatics. 2018. No. 1. P. 52-65. DOI: 10.18101/2304-5728-2018-1-52-65. (In Russ.)
7. Shorikov A.F., Kalev V.I. Linear discrete-time dynamical model forming for solving optimal terminal fuel consumption problem of launch vehicle. Proceedings of the Fifth International Scientific Conference «Information Technologies and Systems» (February, 2428, 2016, Bannoe). Chelyabinsk: Chelyabinsk State University Publ., 2016. P. 61-66. (In Russ.)
8. Bobrovnikov G.N., Katkov A.G. Metody izmereniya urovnya [Methods of level measurement]. Moscow: Mashinostroenie Publ., 1977. 168 p.
9. Bryson A.E., Ho Yu-Chi. Applied optimal control. New York: Routledge Publ., 1975. 496 p.
10. Krasovskiy N.N. Teoriya upravleniya dvizheniem [Motion control theory]. Moscow: Nauka Publ., 1968. 476 p.
11. Krasovskiy N.N., Subbotin A.I. Pozitsionnye differentsial'nye igry [Positional differential games]. Moscow: Nauka Publ., 1974. 456 p.
12. Chernikov S.N. Lineynye neravenstva [Linear inequalities]. Moscow: Nauka Publ., 1968. 488 p.
13. Zoutendijk G. Methods of feasible directions: a study in linear and non-linear programming. Amsterdam: Elsevier Publ., 1960. 126 p.