Научная статья на тему 'Алгоритм оптимального адаптивного терминального управления расходом топлива жидкостной двигательной установки ракеты-носителя. Часть II'

Алгоритм оптимального адаптивного терминального управления расходом топлива жидкостной двигательной установки ракеты-носителя. Часть II Текст научной статьи по специальности «Математика»

CC BY
116
11
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ОПТИМАЛЬНОЕ ПРОГРАММНОЕ И АДАПТИВНОЕ УПРАВЛЕНИЕ / ТЕРМИНАЛЬНОЕ УПРАВЛЕНИЕ / РАСХОД ТОПЛИВА / РАКЕТА-НОСИТЕЛЬ / ОБЛАСТИ ДОСТИЖИМОСТИ / OPTIMAL CLOSED-LOOP CONTROL / TERMINAL CONTROL / PROPELLANT CONSUMPTION / LAUNCH VEHICLE / REACHABLE SETS

Аннотация научной статьи по математике, автор научной работы — Шориков Андрей Федорович, Калёв Виталий Игоревич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Шориков Андрей Федорович, Калёв Виталий Игоревич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

THE OPTIMAL CLOSED-LOOP TERMINAL PROPELLANT CONSUMPTION CONTROL ALGORITHM FOR LIQUID PROPULSION SYSTEM OF LAUNCH VEHICLE. PART II

In this paper, we propose the algorithm of launch vehicle’s liquid propulsion system propellant consumption control problem solving, in which this problem is formulated as the optimal closedloop terminal control problem for corresponding linear discrete-time dynamical system. To achieve this, initial nonlinear continuous model, which describes the plant dynamics, is linearized along given reference trajectory and then is discretized in accordance to control process requirements. As the result of this, we have formed approximating model represented as the vector-matrix system of recurrence equations. The constraints of state vector and control vector are also taken into account in this approximating system, and we assume that the constraints are convex, closed and limited polyhedra with finite number of vertices in corresponding vector spaces. The optimal open-loop and closed-loop terminal control problems are formulated for generated linear discrete-time dynamical system. For solving the problems the optimal closed-loop terminal control recurrence algorithm is developed, which consists in the solving of number of optimal open-loop terminal control problems for linear discrete-time dynamical system. Based on solution of finite sequence of optimal open-loop terminal control problem for linear approximate model we provide recurrence optimal closed-loop terminal control algorithm for initial nonlinear dynamical system. In the implementation of proposed optimal open-loop control algorithm we use the tool of reachable sets computation and analysis based on general algebraic recurrence approach for linear discrete-time dynamical systems and some modifications of this approach intended for reducing the computational complexity and consequently for increasing the operation speed. The performance of developed optimal closed-loop terminal control algorithm is presented on numerical model example of launch vehicle’s first stage liquid propulsion system propellant consumption closed-loop control optimization. This work is a second part of a paper.

Текст научной работы на тему «Алгоритм оптимального адаптивного терминального управления расходом топлива жидкостной двигательной установки ракеты-носителя. Часть II»

УДК 681.511 DOI: 10.14529/^сг190203

АЛГОРИТМ ОПТИМАЛЬНОГО АДАПТИВНОГО ТЕРМИНАЛЬНОГО УПРАВЛЕНИЯ РАСХОДОМ ТОПЛИВА ЖИДКОСТНОЙ ДВИГАТЕЛЬНОЙ УСТАНОВКИ РАКЕТЫ-НОСИТЕЛЯ. ЧАСТЬ II

А.Ф. Шориков1, В.И. Калёв1'2

1 Уральский федеральный университет им. первого Президента России Б.Н. Ельцина, г. Екатеринбург, Россия,

2 АО «НПО автоматики им. академика Н.А. Семихатова», г. Екатеринбург, Россия

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

Данная работа является второй частью статьи.

Ключевые слова: оптимальное программное и адаптивное управление, терминальное управление, расход топлива, ракета-носитель, области достижимости.

Введение

В данной статье на основании работ [1-10] для решения задачи оптимального адаптивного управления дискретной динамической системой предлагается рекуррентный алгоритм, который сводит исходную многошаговую задачу к решению конечной последовательности задач оптимального программного терминального управления в форме Задачи 1. В свою очередь, каждая Задача 1 сводится к реализации конечной последовательности одношаговых операций решения задач линейного и выпуклого математического программирования, алгебраических операций над выпуклыми компактными многогранниками, систем линейных алгебраических операций и вспомогательных одношаговых краевых задач для линейных управляемых систем.

Действие предлагаемых в работе алгоритмов реализовано в программной среде МА^АВ и проиллюстрировано на модельном примере. Работа продолжает первую часть этой статьи, термины и обозначения из которой используются здесь без дополнительных пояснений.

1. Алгоритм формирования множества оптимальных программных терминальных управлений

Пусть в дискретной динамической системе

х^ +1) = А^)х^) + В^)u(t), х(0) = х0 (1)

с ограничениями фазового состояния и управления

х^) е X1(t) с М4, t е 0Т; (2)

ы(0 е U1(t) с М1, t е 0,Т -1, (3)

качество управления в которой на временном интервале т,Т с 0,Т оценивается функционалом у^Г ( ^(т), ы(-) ) = | |х(Т) - |4 =Ф( х(Т)), (4)

сформирована т-позиция w(^e) (т) = {т, х(е) (т)} е W(т) (w(e) (0) = w0 ). Тогда для вычисления множества оптимальных программных управлений и(е)(т,Т,{х(т)}, (Т)) с и(т,Т) разработан метод

Ут,Т

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

1) пошагового построения прямых областей достижимости G + (т, х(е)(т);^ , Vtет + 1,Т на основе модификации общего рекуррентного алгебраического метода [10] путем решения рекуррентной последовательности задач линейного математического программирования, выполнения одношаговых операций над компактными многогранниками и описания их с помощью вершин и конечных систем алгебраических равенств и неравенств;

2) оптимизации функционала у—т ^(т),ы(-)), определяемого (4), на множестве G + (т, х(е)(т);Т), т. е. формирования множества Х(е)_ (Т) с М4 оптимальных финальных состояний системы (1)-(4)

Ут,Т

путем решения следующей задачи выпуклого математического программирования [6, 11]:

Х(еЛ (Т) = | х(е) (Т) | х(е) (Т) е G + (т, х(е) (т);Т),

'т,Т

c<e)(x,T, w(x)) = = min x(e) (T) - x,

x,T x (T )eG + (x,x(e)(x);T )

Де),

3) пошагового построения обратных областей достижимости О _ (Т, Х(~е^(Т); t),

Ут,Т

Vt е {Т _ 1,Т _ 2, ..., т +1} на основе модификации общего рекуррентного алгебраического метода [10] путем решения рекуррентной последовательности задач линейного математического программирования, выполнения одношаговых операций над компактными многогранниками и описания их с помощью вершин и конечных систем алгебраических равенств и неравенств;

4) последовательного построения пересечений сформированных прямых и обратных областей достижимости путем реализации конечной последовательности одношаговых операций перехода от вершинного описания областей достижимости к описанию соответствующими системами линейных алгебраических уравнений и неравенств и наоборот:

О(т, х(т); t) = О + (т, х(т);I) П О_(Т,Х(1 (Т);t), Vt е т,Т _ 1;

т,Т

О(т, х(т); т) = {х(т)}; О(т, х(т);Т) = Х(1 (Т);

5) решения вспомогательных краевых задач на промежутке времени т,Т для системы (1)-(4) с краевыми условиями х(т) и х(Т) = х^е)(Т) е Х^ (Т), т. е. нахождения следующего множества

Ут,Т

оптимальных программных управлений путем реализации конечной последовательности алгебраических операций:

U(e)(x,T,{x(x)},X<1 (T)) = {u(e)(-) | u(e)(t) e Ui(t), Vt ex,T -1,

U,t

x (t +1) = ^(t)x (t) + B(t)u(e) (t) e G(x, x(x); t +1), x(e)(x) e X1 (x)}.

4

2. Алгоритм формирования оптимального адаптивного терминального управления

Пусть в дискретной динамической системе (1)-(4) сформирована т-позиция w(e) (т) = {т,х(е) (т)} е W(т) (w(e) (0) = w0 ).

Тогда для те 0,Т-1 оптимальная адаптивная стратегия управления ^^ = U1) е формируется на основе применения метода прямых и обратных конструкций [3, 6] путем реализации следующей последовательности операций.

1. Решения Задачи 1, т. е. вычисления множества

^е) (ТТ, w(т)) = ^е) (ТТ, {х(т)}, X™ (Т)) е ЩТТ)

' т,Т

и числа сУе)(т,Т, w(т)) = фТу на отрезке времени т,Т .

2. На основании формул (13)-(15) первой части статьи вычисления множества управлений ^^Чт)) = U1%(е)(т)) е U1(т).

3. Выбора любого управления и(е)(т) е Uae)(w(e)(т)) .

4. Вычисления по уравнениям нелинейного объекта (т+1)-позиции: w(e) (Т +1) = {т +1, х(е) (т +1)}, х(е) (т +1) = х(т +1; т, т +1, х(е) (т),и(е) (т)).

5. Если (т +1) < Т -1, то осуществить переход на п. 1, иначе - на п. 6.

6. Расчет числа с1еУ(0,Т,w0) = у^-уи(е)( )) - оптимального результата решения Задачи 2, соответствующего реализации управления и(е)( ) = {и(е)^)}?е0тт-1 е U(0,T), полученного посред-

тт(е) тт*

ством применения оптимальной адаптивной стратегии ^ = а е Ua .

В работах [3, 6] доказано, что для линейных дискретных систем вида (1)-(3) предложенный алгоритм адаптивного управления, основанный на решении вспомогательных задач оптимального программного терминального управления методом прямых и обратных конструкций, является оптимальным алгоритмом адаптивного управления для терминального функционала вида (4).

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

3. Численный модельный пример решения задачи оптимального адаптивного терминального управления ДУ РН

Покажем эффективность предложенного подхода и соответствующих алгоритмов на численном примере, моделирующем решение задачи оптимального адаптивного терминального управления расходом топлива ДУ РН [8, 9]. Рассматривается нелинейная непрерывная управляемая система, описывающая динамику ДУ РН на промежутке времени [0, уравнениями:

т (Л = (Р + ОЦЦ (О2 + 02Ц (} ))(К + и и))

^О (') О ,

а+Съи ^ )2+сии Ш1+к+и ^))

т, ^) =_Р + с!и (t )2 + с2и (t)_,

1 (I + с3и ^ )2 + с4и (t))(1 + К + и ^)) М0 ^) = М°° +АИо -Г'1 т0 (t)dt,

(6)

Мг ^) = М° + АМГ - Г1 тг ^)dt, с параметрами, заданными в табл. 1.

Таблица 1

Параметры исходной нелинейной системы

I, c P, кг K, б/р M0, кг M0, кг AMo, кг AM, кг

320 528 000 8/3 120 000 45 000 -500 300

AK, б/р C1, б/р C2, б/р c3, б/р c4, б/р c5, б/р t1, с

-1/25 8/75 -6700 9000 -15 -12 100

Исходной нелинейной непрерывной системе поставлена в соответствие аппроксимирующая линейная дискретная динамическая система, полученная посредством линеаризации вдоль опорной траектории вида траектории:

m

ref (t) =

m

ref

(t) =

P ■ K I +1 ■ K' P ; I +1 ■ K'

Mf (t) = M0 -Mf (t) = M0 -

P ■ K I +1 ■ K

_P_

I +1 ■ K

-t,

-t,

(7)

и дискретизации. Таким образом, аппроксимирующая рекуррентная векторно-матричная система, соответствующая (1), имеет вид

х^ +1) = Л(0х^) + ), t е 0,15, где х(0 = {x1(t),х2(0,х3^),x4(t)}е М4; х^) - расход окислителя ДУ РН; х2(0 - масса окислителя ДУ; xз(t) - расход горючего ДУ РН; х4(0 - масса горючего ДУ РН; и({) е М1 - вектор управления; х(0) = х0 = (1192.5,119 500, 454, 45 300)г - начальный фазовый вектор объекта управления; матрицы Л(^) и В(^) стационарные и имеют следующие значения:

A(t) =

Г 1 0 0 0 > Г 17 >

-6,25 1 0 0 0

, B(t) = , Vt е 0,15

0 0 1 0 -12

ч 0 0 -6,25 1V ч 0 V

Фазовые ограничения ДУ РН для всех t е 0,16 имеют вид

х1^) е [1176; 1224], |х2(0| < 120 000, х3^) е [432; 468], |х4(0| < 46 000.

Качество процесса управления жидкостной ДУ РН оценивается значением выпуклого терминального функционала в момент времени Т = 16, соответствующий моменту времени t\ = 100 для начальной нелинейной непрерывной системы, значения которого вычисляются по формуле

Ф(х(Т)) = ^(1200 - х1(Т))2 + х2(Т)2 + К2 (450- х3(Т))2 + К2х4(Т)2 ,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

где первый и третий члены под квадратным корнем обозначают отклонение расходов компонентов топлива от требуемых значений, а второй и четвертый члены - количество топлива в баках жидкостной ДУ РН.

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

и(е)(-) = {и1)(-)} ^, У/ е 1,к: и1)(-) = = {и^ 1 ЧО}^^, которое порождает соответствующие

оптимальные траектории x(e, 1)(-) = х(• ;0,16,х0,u(e, 1)(-)), ] е 1,к . Финальное состояние таких траекторий - х^е'1 )(16) = х(16;0,16, х0, и -^(О), для всех у е 1, к есть вычисленный вектор

Де, j )

(16) = (1200.1, 38.97, 450.03, 10.5), который доставляет функционалу качества Ф минималь-

ное значение Ф(е) =Ф(x(e)(16)) = 47,988. Далее выбираются четыре программных управления

Де)

и (е,/), / = 1,4 и подставляются в исходную систему (5), (6) для сравнения результатов реализации процедур программного и адаптивного управления.

Решение задачи оптимального адаптивного терминального управления нелинейной системой реализовано с помощью алгоритма, приведенного во второй главе данной работы, и описывается множеством управлений и(ае)(-) = {и!?^Ч-)}/^ , У/ е1,d : и^е^)(-) = {и1е,/Чх^))}^^. Для сравнения полученных результатов реализации программного и адаптивного управления в рассматриваемом примере были выбраны и подставлены пошагово в систему (5), (6) четыре стратегии адаптивного управления и^-7) (х^)), / е 1,4 . Проекции соответствующих траекторий представлены на рис. 1.

Рис. 1. Проекции фазовых траекторий при управлениях:

а) и(еД) и ие'1); б) и

(еД);

/е'2) и и(е'2);

; в) и(е '3) и ие'3); г) и

(е,3) ;

/е'4) и и(е'4)

В табл. 2 приведены значения функционала Ф для траекторий, соответствующих реализациям оптимального программного терминального управления линейной и нелинейной системами, а также адаптивного терминального управления исходной нелинейной системой.

Период времени моделирования для одного варианта решения задачи оптимального программного терминального управления жидкостной ДУ РН составляет 0,4 с. Весь период времени моделирования решения задачи оптимального адаптивного терминального управления составляет 10 с. Из табл. 2 видно, что для исходной нелинейной системы результат реализации оптимального адаптивного терминального управления лучше, чем результат реализации оптимального программного терминального управления.

Таблица 2

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

Управление xi(T) xi(T) x3(T) x4(T Ф^(Т)

Опорная траектория u(t) = 0 1200 0 450 0 0

Нелинейная система u fe1)(t) 1200,1 39,1 450 10,6 48,261

x(t)) 1200,2 32,7 449,9 8,5 39,837

u (e'2)(t) 1200,1 39,6 450 12,1 51,021

uae,2)( x(t)) 1200,1 36,3 450 9,5 44,228

u (e'3)(t) 1200 38,7 450 9,6 46,339

uae,3)( x(t)) 1200,1 31,1 450 8,2 37,935

u (e'4)(t) 1200,1 39,9 450 10,6 48,906

u{:,4)( x(t)) 1200 36,6 450 10,1 45,372

Линейная система uj }(t) 1200,1 39 450 10,5 47,988

Численное моделирование примера решения задач оптимального программного и адаптивного терминального управления жидкостной ДУ РН реализовано в среде MATLAB R2014a на ПК с Intel© Core™ i7-3770 CPU @ 3.4 GHz, 8 Gb RAM.

Заключение

В статье динамическая система управления расходом топлива жидкостной ДУ РН описана как нелинейная непрерывная дифференциальная система. Задача терминального управления расходом топлива жидкостной ДУ РН сформулирована как задача оптимального адаптивного терминального управления соответствующей аппроксимирующей линейной дискретной динамической системой. Для вычисления оптимального адаптивного терминального управления ДУ РН с нелинейной динамикой предложена процедура, основанная на реализации решения конечной рекуррентной последовательности вспомогательных задач оптимального программного терминального управления соответствующей аппроксимирующей линейной дискретной системой. Подробно описаны алгоритмы решения основной и вспомогательной задач оптимального управления. Показано, что задачу оптимального адаптивного терминального управления можно свести к решению конечного числа задач линейного и выпуклого математического программирования, выполнения одношаговых операций над компактными многогранниками и описания их с помощью вершин и конечных систем алгебраических равенств и неравенств, решения вспомогательных краевых задач. Результаты компьютерного моделирования показали применимость общего рекуррентного алгебраического метода построения областей достижимости и метода прямых и обратных конструкций [3] к решению задачи оптимального адаптивного терминального управления расходом топлива жидкостной ДУ РН.

Работа выполнена при финансовой поддержке РФФИ (проект № 18-01-00544).

Литература

1. Бортовые терминальные системы управления (принципы построения и элементы теории) / Б.Н. Петров, Ю.П. Портнов-Соколов, А.Я. Андриенко, В.П. Иванов. - М.: Машиностроение, 1983. - 200 с.

2. Красовский, Н.Н. Теория управления движением / Н.Н. Красовский. - М.: Наука, 1968. -476 с.

3. Шориков, А.Ф. Минимаксное оценивание и управление в дискретных динамических системах /А.Ф. Шориков. - Екатеринбург: Изд-во Урал. гос. ун-та, 1997. - 242 с.

4. Шориков, А.Ф. Формирование линейной дискретной динамической модели для решения задачи оптимального терминального управления расходом топлива ракеты-носителя / А.Ф. Шориков, Калёв В.И. // Информационные технологии и системы: тр. 5-й Междунар. науч. конф. -2016. - С. 61-66.

5. Калёв, В.И. Моделирование задачи терминального управления расходом топлива жидкостных ракет / В.И. Калёв, А.Ф. Шориков // Известия высших учебных заведений. Физика. - 2016. -Т. 59, № 8-2. - С. 45-48.

6. Тюлюкин, В.А. Об одном алгоритме построения области достижимости линейной управляемой системы /В.А. Тюлюкин, А.Ф. Шориков //Негладкие задачи оптимизации и управление. -Свердловск: УрО АН СССР. - 1988. - С. 55-63.

7. Тюлюкин, В.А. Алгоритм решения задачи терминального управления для линейной дискретной системы /В.А. Тюлюкин, А.Ф. Шориков //Автоматика и телемеханика. - 1993. - № 4. -С. 115-127.

8. Шориков, А.Ф. Алгоритм решения задачи оптимального терминального управления в линейных дискретных динамических системах / А.Ф. Шориков // Информационные технологии в экономике: теория, модели и методы: сб. науч. тр. - Екатеринбург: Изд-во Урал. гос. экон. ун-та, 2005. - С. 119-138.

9. Шориков, А.Ф. Описание библиотеки компьютерных программ для моделирования решения задачи апостериорного минимаксного оценивания /А.Ф. Шориков, В.А. Тюлюкин //Известия Урал. гос. экон. ун-та. - 1999. - № 2. - С. 36-49.

10. Аппроксимация областей достижимости нелинейных дискретных управляемых динамических систем /А.Ф. Шориков, В.В. Булаев, А.Ю. Горанов, В.И. Калев //Вестник БГУ. Математика, информатика. - 2018. - № 1. - С. 52-65.

11. Зойтендейк Г. Методы возможных направлений / Г. Зойтендейк. - М. : Изд-во иностр. лит., 1963. - 176 с.

Шориков Андрей Федорович, д-р физ.-мат. наук, профессор кафедры прикладной математики Уральского энергетического института, Уральский федеральный университет им. первого Президента России Б.Н. Ельцина, г. Екатеринбург; afshorikov@mail.ru.

Калёв Виталий Игоревич, аспирант кафедры прикладной математики Уральского энергетического института, Уральский федеральный университет им. первого Президента России Б.Н. Ельцина; инженер-конструктор, АО «НПО автоматики им. академика Н.А. Семихатова», г. Екатеринбург; v.i.kalev@urfu.ru.

Поступила в редакцию 25 января 2019 г.

DOI: 10.14529/ctcr190203

THE OPTIMAL CLOSED-LOOP TERMINAL PROPELLANT CONSUMPTION CONTROL ALGORITHM FOR LIQUID PROPULSION SYSTEM OF LAUNCH VEHICLE. PART II

A.F. Shorikov1, afshorikov@mail.ru, V.I. Kalev1'2, v.i.kalev@urfu.ru

1 Ural Federal University named after the first President of Russia B.N. Yeltsin, Ekaterinburg, Russian Federation,

2 JSC "Scientific and Production Association of automatics named after academician N.A. Semikhatov", Ekaterinburg, Russian Federation

In this paper, we propose the algorithm of launch vehicle's liquid propulsion system propellant consumption control problem solving, in which this problem is formulated as the optimal closed-loop terminal control problem for corresponding linear discrete-time dynamical system. To achieve this, initial nonlinear continuous model, which describes the plant dynamics, is linearized along given reference trajectory and then is discretized in accordance to control process requirements. As the result of this, we have formed approximating model represented as the vector-matrix system of recurrence equations. The constraints of state vector and control vector are also taken into account in this approximating system, and we assume that the constraints are convex, closed and limited polyhedra with finite number of vertices in corresponding vector spaces. The optimal open-loop and closed-loop terminal control problems are formulated for generated linear discrete-time dynamical system. For solving the problems the optimal closed-loop terminal control recurrence algorithm is developed, which consists in the solving of number of optimal open-loop terminal control problems for linear discrete-time dynamical system. Based on solution of finite sequence of optimal open-loop terminal control problem for linear approximate model we provide recurrence optimal closed-loop terminal control algorithm for initial nonlinear dynamical system. In the implementation of proposed optimal open-loop control algorithm we use the tool of reachable sets computation and analysis based on general algebraic recurrence approach for linear discrete-time dynamical systems and some modifications of this approach intended for reducing the computational complexity and consequently for increasing the operation speed. The performance of developed optimal closed-loop terminal control algorithm is presented on numerical model example of launch vehicle's first stage liquid propulsion system propellant consumption closed-loop control optimization.

This work is a second part of a paper.

Keywords: optimal closed-loop control, terminal control, propellant consumption, launch vehicle, reachable sets.

References

1. Petrov B.N., Portnov-Sokolov A.Yu., Andrienko A.Ya., Ivanov V.P. Bortovye terminalnye sistemy upravleniya [On-board Terminal Control Systems]. Moscow, Mashinostroenie Publ., 1983. 200 p.

2. Krasovskiy N.N. Teoriya upravleniya dvizheniem (Lineynye sistemy) [Theory of Motion Control (Linear Systems)]. Moscow, Nauka Publ., 1968. 476 p.

3. Shorikov A.F. Minimaksnoe otsenivanie i upravlenie v diskretnykh dinamicheskikh sistemakh [Minimax Estimation and Control in Discrete-Time Dynamical Systems]. Ekaterinburg, Ural State University Publ., 1997. 242 p.

4. Shorikov A.F., Kalev V.I. [Linear Discrete-Time Dynamical Model Forming for Solving Optimal Terminal Fuel Consumption Problem of Launch Vehicle]. Informatsionnye tekhnologii i sistemy. Proceedings of 5th International Conference, 2016, pp. 61-66. (in Russ.)

5. Kalev V.I., Shorikov A.F. [Fuel Consumption Terminal Control Problem Statement for Liquid-Propellant Rockets]. Russian Physics Journal, 2016, vol. 59, no. 8-2, pp. 45-48. (in Russ.)

6. Tyulyukin V.A., Shorikov A.F. [About One Algorithm of Creation of Area of Approachibility of the Linear Operated System]. Optimization and Control Nonsmooth Problems, 1988, pp. 55-63. (in Russ.)

7. Tyulyukin V.A., Shorikov A.F. [The Solution Algorithm of Terminal Control Problem for Linear Discrete-Time System]. Automatics and Telemechanics, 1993, no. 4, pp. 115-127. (in Russ.)

8. Shorikov A.F. [Algorithm of the Solution of a Problem of Optimum Terminal Control in Linear Discrete Dynamic Systems]. Information Technologies in Economy: Theory, Models and Methods: Collection of Scientific Works, Ekaterinburg, 2005, pp. 119-138. (in Russ.)

9. Shorikov A.F., Tyulyukin V.A. [Description of Library of Computer Programs for Modeling of the Solution of a Problem of a Posteriori Minimax Estimation]. News the Ural State Econ. Univ., 1999, no 2, pp. 36-49. (in Russ.)

10. Shorikov A.F., Bulaev V.V., Goranov A.Yu., Kalev V.I., [Approximation of Attainability Domains of Nonlinear Discrete-Time Controlled Dynamical Systems]. Bulletin of the Buryat State University. Mathematics, Informatics, 2018, no. 1, pp. 52-65. (in Russ.)

11. Zoutendijk G. Metody vozmozhnykh napravleniy [Methods of Feasible Direction]. Moscow, Inostrannaya literatura Publ, 1963. 176 p.

Received 25 January 2019

ОБРАЗЕЦ ЦИТИРОВАНИЯ

FOR CITATION

Шориков, А.Ф. Алгоритм оптимального адаптивного терминального управления расходом топлива жидкостной двигательной установки ракеты-носителя. Часть II / А.Ф. Шориков, В.И. Калёв // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2019. - Т. 19, № 2. -

С. 29-37. DOI: 10.14529/ctcr190203

Shorikov A.F., Kalev V.I. The Optimal Closed-Loop Terminal Propellant Consumption Control Algorithm for Liquid Propulsion System of Launch Vehicle. Part II. Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics,, 2019, vol. 19, no. 2, pp. 29-37. (in Russ.) DOI: 10.14529/ctcr190203

i Надоели баннеры? Вы всегда можете отключить рекламу.