CC BY
23
UDC 629.764.7
Siberian Journal of Science and Technology. 2017, Vol. 18, No. 2, P. 314-322
DEVELOPMENT OF LAUNCH VEHICLE CONTROL ALGORITHM FOR THE INITIAL PART OF THE TRAJECTORY USING THE ACOR METHOD
A. Sh. Altshuler1, V. T. Bobronnikov2, M. V. Trifonov2*
1Khrunichev State Research and Production Space Center 18, Novozavodskaya Str., Moscow, 121087, Russian Federation 2Moscow Aviation Institute (National Research University) 4, Volokolamskoe Shosse, A-80, GSP-3, Moscow, 125993, Russian Federation
E-mail: [email protected]
The control system of a launch vehicle (LV) at the initial phase offlight at altitudes between 0 and 300 meters is the object of investigation in the article. At this phase the flame of the LV jet engine is the cause of a negative impact on facilities of the launch complex. This effect can be reduced by displacing the jet flame in a certain radial direction with increased resistance of the facilities using a specially developed vehicle motion control program.
The purpose of the article is to develop an algorithm for the controller of the LV motion control system on the considered "displacement phase" of the trajectory that provides an implementation of such a program.
To solve the problem, we developed the modified version of the Letov's method of analytical construction of regulators (ACOR). The peculiarity of the modified statement of the problem solved in the work is that the controlled output vector of the system depends explicitly not only on the LV state vector, but also on the control variable.
The quality of control is evaluated using a quadratic terminal-integral optimality criterion. This kind of criterion allows to trace with the specified accuracy the preliminary calculated program for supporting the required position of the trace of the LV jet flame on the launching plane, and also to ensure the vertical position and the zero angular velocity of the vehicle at the end of the displacement phase.
To solve the problem of constructing the algorithm, a special linearized model of the LV motion has been developed. The results of simulating the controlled motion of a launch vehicle with the use of the algorithm confirm the operability and demonstrate efficiency of the developed optimal regulator of the LV control system at the displacement phase under consideration. Calculation results show that angular position of the LV at the end of the displacement phase is close to vertical, the angle of the engine nozzle deflection is within the permissible limits and the deviation of the current position of the jet flame track from the program value does not exceed 0.5 meters.
Keywords: launch vehicle, facilities of the launch complex, jet flame, displacement phase, optimal regulator, quadratic criterion.
Сибирский журнал науки и технологий. 2017. Т. 18, № 2. С. 314-322
РАЗРАБОТКА АЛГОРИТМА УПРАВЛЕНИЯ ДВИЖЕНИЕМ РАКЕТЫ-НОСИТЕЛЯ НА НАЧАЛЬНОМ УЧАСТКЕ ПОЛЕТА С ИСПОЛЬЗОВАНИЕМ МЕТОДА АКОР
А. Ш. Альтшулер1, В. Т. Бобронников2, М. В. Трифонов2*
Государственный космический научно-производственный центр имени М. В. Хруничева Российская Федерация, 121087, г. Москва, ул. Новозаводская, 18 2Московский авиационный институт (национальный исследовательский университет) Российская Федерация, 125993, г. Москва, Волоколамское шоссе, 4, А-80, ГСП-3
E-mail: [email protected]
В качестве объекта исследования рассматривается система управления движением РН на начальном участке полета 0-300 метров. На этом участке газодинамическая струя двигателя РН оказывает негативное воздействие на сооружения стартового комплекса. Это воздействие может быть снижено за счет увода струи в радиальном направлении от точки старта по заранее заданной программе управления.
Целью исследования является разработка оптимального регулятора СУ РН на этапе увода, обеспечивающего реализацию программы увода.
Предлагается модифицированный вариант решения задачи аналитического конструирования регуляторов А. М. Летова. Особенностью постановки задачи является то, что вектор выхода системы зависит и от вектора состояния РН, и от управляющей переменной, и качество управления в задаче оценивается с использованием квадратичного терминально-интегрального критерия. Такой вид критерия позволяет с заданной точностью отслеживать требуемую программу изменения следа газодинамической струи РН на стартовой плоскости и обеспечить заданное положение РН в конце этапа увода.
Задача решается с использованием линеаризованных уравнений движения РН. Результаты моделирования движения РН подтверждают работоспособность и эффективность разработанного оптимального регулятора системы управления РН. Расчеты показывают, что движение РН на этапе увода близко к вертикальному, угол отклонения сопла двигателя находится в допустимом диапазоне и величина рассогласования текущего положения следа струи с заданным программным не превышает 0,5 метра.
Ключевые слова: ракета-носитель, сооружения стартового комплекса, газодинамическая струя, этап увода, оптимальный регулятор, квадратичный критерий.
Introduction. At the initial stage of the launch vehicle (LV) flight at altitudes between 200-300 meters one of the flight control system's tasks is to displace the launch vehicle's jet engine flame from the launch complex (LC). This process is of great importance as the jet engine flame has a negative impact on the LC [1].
There are several methods to reduce negative thermal effect on LC. Some of them are:
1. Methods of "passive" protection: to use heatproof materials in the constructions of the LC or to use additional measures such as foam or water to cool the launching complex elements [2; 3].
2. To use some modified algorithms of a LV flight control to displace LV jet flame in specially organized direction and so on.
For example, at the LV of ground start "ZENIT" to cool jet engine flame and to protect LC facilities water is used [2]. While at the LV of sea launch "ZENIT-3SL" the technology of jet flame deviation to certain direction by means of modified algorithms of the LV flight control is used [4-6].
Two approaches can be proposed to perform controlled displacement of a LV trajectory using the vehicle flight control system (CS):
1. Using pre-calculated motion program of a LV flight at the displacement phase (DP) taking into account initial and terminal (boundary) states of the vehicle and the subsequent executing of this program with use of the CS of the vehicle.
2. Solving the synthesis problem of flight control taking into account the current and required terminal LV states.
In this paper we consider the second approach - synthesis of the LV fligth control algorithm. Solving the problem of jet engine flame displacement, the position of flame trace on the launch pad is one of output parameters of the LV flight control system.
To implement the controlled DP of the LV flight, the program of flame flow location on the launch pad in the direction of displacement lp(h) as a function the LV altitude h has to be preliminarily calculated. This program is fulfilled by the LV CS with flight control algorithm (regulator), specially developed for this phase of the LV flight.
In each time of DP the jet engine flame trace must be located inside a circular area of a certain size. The center of the area must be located on the line specifying the direction of the displacement. Size of the area has to be calculated in advance as a function of time for every particular LV and LC. Any high-rise structures should not be located in the direction of displacement.
Generally the model of LV motion at the initial phase of LV flight is described by nonlinear equations. But
while solving the problem of LV jet flame displacement from launching complex constructions this model can be replaced by linearized equations with time variant parameters.
The problem of finding the optimal control of the system described by linear equations can be solved using the method of analytical construction of regulators (ACOR). The first scientific papers on this theme belong to A. M. Letov, published in 1960-s [7; 8]. In foreign publications this method is known as the tracking task method.
The aim of this work is to develop the method of construction of optimal regulator of a LV control system for the DP that provides the executing the program of jet flame trace displacement.
The problem is solved for linear equations describing the launch vehicle motion in the direction of displacement. The new modification of the ACOR method is proposed. The modified version of the method allows finding the solution of the problem when the system output is dependent explicitly on input variable.
Statement of the problem. While starting the LV nearest facility of the LC is the cable-filling tower. In order to prevent the LV collision with the tower the program of jet flame trace displacement on the launch pad as a function of trajectory altitude is used (fig. 1).
0 50 100 150 200 250 ЗОВ
fe m
Fig. 1. Preset program changes the LV jet flame trace location as a function of flight altitude
Рис. 1. Заданная программа изменения положения следа оси ГД струи РН в зависимости от высоты полета
Double-stage character of this program is stipulated with the fact that after the lift-off the LV should move vertically to avoid hitting the cable-filling tower, then softly move in the direction of displacement on preset distance from starting point not setting engine plume close to the tower.
The control variable in the jet engine flame displacement problem is the engine nozzle deviation angle 5.
The task is to develop an optimal LV motion control law taking into account LV parameters and preset program of jet flame trace location at the launch pad, to choose the structure and gains of the LV CS regulator for executing the preset program of the jet flame trace displacement.
Efficiency of control is evaluated by quadratic terminal integral criteria. The terminal part of criteria is used to provide preset final angular state of the LV at the end of DP - vertical orientation and zero rate of LV. This requirement is caused by the fact that after completing the DP, the LV must continue its motion according to the regular pitch program. The integral part is some kind of a "penalty" for deviation of current output parameter lj, from its preset program value l3 [10].
The model of LV motion at the phase of controlled displacement. Parameters describing motion of the LV in the displacement plane are depicted in fig. 2.
Fig. 2. Parameters of LV motion at the displacement phase in the starting frame
Рис. 2. Параметры движения РН на этапе увода в стартовой системе координат
Fig. 2 shows the main parameters of LV motion [11]: OL is the LV jet flame displacement direction, RA is the main vector of external forces, acting on LV, С is LV center of mass, F is LV pressure center.
Assumption:
1) LV motion is performed in the plane of displacement maneuver;
2) LV is considered as a solid body;
3) the Earth gravitational field is homogeneous;
4) the Earth is plain and nonrotating;
5) atmospheric density during the DP is constant, wind disturbances are not considered;
6) LV mass and inertia moments during the DP are constant;
7) LV aerodynamic characteristics are: Cx = const, Cy = C"a , Cd = const , where a - angle of attack; Cx -aerodynamic drag force coefficient; Cy - aerodynamic normal force coefficient; C - derivative of aerodynamic normal force coefficient; Cd - coefficient of pressure center;
8) LV engine thrust during the DP is constant;
9) the dynamics of engine thrust vector control actuator is neglected.
LV motion equation. Several reference frames are used in the article for description of LV motion (fig. 2): the starting frame OXCYCZC, the body frame CXYZ, the flow frame CXaYaZa [12]. The beginning O of starting frame is located in the crossing point of LV longitudinal axis, when LV is located on the launch pad (LP), with horizontal plane on zero level of launch pad. Axis OZC form the right-handed frame OXCYCZC.
With enumerated assumptions the LV motion in the displacement maneuver plane is described by the following differential equations:
Vh =(P sin (S +8 ) -mg-Xa sinS + Ya cosS)/m, V =(P cos (S + 8) - Ya sin S-Xa cos S)/m, ra = (-P sin 8 (xT -xB) + Ya (Xd - xT)) /I, l = V, h = Vh, S = ra,
where m - LV mass; I - moment of inertia relative to lateral axis; Vh h Vl - projections of velocity
on flow axes h and l; P - engine thrust; S - pitch angle; 9 = arctg (Vh/Vl) - trajectory inclination angle; Xa = CXSpV2 /2 - aerodynamic longitudinal force; Ya = C^a SpV2/2 - aerodynamic normal force, a = S-9; h - LV center of mass (COG) flight altitude; l - distance from starting point to LV COG projection on the launch pad plane; ra - angular rate of LV relative to lateral axis; xT - distance from nozzle edge plane of sus-tainer engine to COG; xB - distance from nozzle edge plane to center of deviation of sustainer engine chamber; xd - distance from nozzle edge plane of sustainer engine to LV pressure center: xd = L(1 - Cd); L - characteristic length; S - characteristic area; p - atmospheric density; 5 - angle of deviation of sustainer engine chamber; g - gravitational acceleration.
Simplifying and linearization of LV motion model. Considering values 5, a, AS = S-n/2, A9 = 9-tc/2 small, we will get:
sin S «1, cos S « -AS, sin 8 « 8, cos 8 «1, sin (S+8 )«1, cos (S + 8) « -AS-8, Vh = V sin 9 « V, Vt = V cos 9«-V A9.
With these assumptions the system (1) transformed to the following two systems:
V> = ( p - mg - CxS pV2/2) /m, h = V.
V = (-PДâ - P8- Cj SpV2 а /2 + CXSpV2ДЗ /2)/m, / = V,
CD = (-P8(xT — хВ ) + CjaaS(xd — xT )pV2 /2)//, ДЗ = ю,
(3)
where a = AS - A8, A9 = - V^V.
The system (3) can be rewritten in the following form:
V =
i
(CX - Ca)SpV2 - 2P ДЗ CJSpVV P
2m
V--8,
2m m
I = V,
. = CY(xd - xt)SPV2 ДЗ + 21
Ca (xd - xt )SP£V - P(Xt - xB )
(4)
8, ДЗ = ю .
2I I
System (2) is independent from the system (4). It can be integrated with initial conditions: V(0) = 0 and h(0) = h0 = xT to get explicit expressions for V and h as functions of time.
Several transformations of the first equation in the system (2) give the following results:
m -V = 1__Cx P^ V 2;
P - mg 2(P - mg)
CxPS
a =-, b = -
P - mg 2(P - mg)
adV
aV = 1 -bV
adV
J1 - bV2 =f
= dt,
a 1 + 4b V t C ln-;=-= t + C.
24b i-4b V As V(0) = 0, so C = 0 and
V (t ) t;
4b a
h(t ) = h0 + f Vdt = h0 + a ln ch—t. J b a
(5)
so
lp = l - (xT - xB )cos З-
- (h - (xT - xB ) sin З) ctg (З + 8).
At small angles 5 and AS = S -—
lp = l + A (AS + 8) - (xT - xB) 8 .
So the final formula for calculation of output lp is:
lp = l + AAS + hB8,
(6)
(7)
where h is the known function of time (5); hB = h - xT + xB -altitude of sustainer engine chamber center of deviation above the starting plane.
Preset coordinate of LV jet flow trace on the starting plane is denoted as l3(hc), where hc - altitude of nozzle edge plane of the sustainer engine above the starting plane.
hc = h - xT. (8)
Formulas (5), (8) and l3(hc) can be used to recalculate the preset coordinate of LV jet flame trace along x axis as a function of flight time.
Motion equation of LV in vector-matrix form. System (4) can be written in normal Cauchy form:
dx(t) = A(t) x(t) + B(t )u(t),
dt
(9)
where x = \Vl ю l ДЗ]Т isstate vector;
A =
a11 0 0 a14
a21 0 0 a24
1 0 0 0
0 1 0 0
- dynamic matrix with elements:
One of system's output parameter is the coordinate of jet flow trace lp along axis OL.
Coordinates sustainer chamber center of deviations on the on the starting frame axes are described as follows:
11 ik = l ( Xt Xb ) cos S,
h]^K = h - (xT - xB) sin S.
Equation of the straight-line H = H(L), that goes through the sustainer chamber center of deviations and is parallel to thrust vector:
H - hjK = (L -Ijjk)tg (S + 8).
Coordinate of jet flame trace on the starting plane can be evaluated with H = 0:
lp = V - hm ctg (S+8^
a11 =-Ca SpV /2m, a14 =[(CX -Ca) S pV2 -2P a21 = Ca (xd - xT )SpV/21, a24 = C^a(xd - xT )SpV2/2/;
B = [bi b2 0 0] - vector with elements:
b =-P / m, b2 =-P(xt - xB)/1; u = [5] - control variable.
The system output vector is described by the algebraic
equation:
>>(t) = C (t) x(t) + D(t)u(t), (10)
where y = [V ю lp Д^]Т ters;
C =
1 0 0 0
0 1 0 0
0 0 1 h
0 0 0 1
vector of output parame-
D = [0 0 hB 0]T.
vector
T
of output parameters is
Preset
z = [0 0 l3 0]T .
Method of analytical construction of regulators (ACOR). Let's consider linear time variable observable and controllable dynamic system (9); (10) with the initial conditions x(t0) = x0 [9; 13]. In this system x(t) -«-dimensional state vector, u(t) - r-dimensional vector of
control variables, y(t) - m-dimensional vector of outputs, z(t) - m-dimensional vector of demanded outputs.
Vector equation of tracking error at the time t is calculated as follows:
e(t) = y(t) - z(t). It is necessary to find the control minimizing the criteria:
J (e, u) = 2 eT (tK )Fe(tK ) +
1 1к
— J (eT (t )Q(t )e(t ) + uT (t)R(t)u(t) ) dt.
(11)
In (11) F, R, Q - preset "weight" matrices, control variables u(t) are not constrained, tk - preset final time.
It is shown [9; 13], that the solution of the problem if D = 0 gives to the optimal control law in the form:
u (t) = R-1(t) BT (t) [ g (t) - K (t) x(t)],
where K(t) is «-dimensional quadratic symmetric positive definite matrix, that obey to the matrix differential Riccatti equation:
dK (t)
dt
■ = K (t)B(t)R-1 (t)BT (t)K(t) -
dt
- [ K (t ) B(t ) R-1 (t ) BT - AT (t ) ] g (t ) - CT (t )Q(t ) z (t )
with terminal condition:
g (tk ) = CT (tk )Fz (tk ).
But such statement of the problem is not general. In the technical task considered in this article, output vector y depends both on state vector x and input variable u = l3.
To find the solution in the such case the ACOR problem has to be reformulated in the following more general form: to find the optimal control law u of linear dynamic system (9), (10) with D(t)u(t) term, that minimizes criteria (11). The solution in this case we will call as solution of the modified ACOR problem.
Solution of modified ACOR problem. According to the Pontryagin's principle of maximum the optimal control of linear system (9), (10) must minimize Hamiltonian [14]:
1 T
H(p, x, u, t) = z - Cx - Du) x
x Q (z - Cx - Du)+ uTRu + pT Ax + pT Bu, where vector-function p(t) satisfies the vector equation
dpT = aH
dt dx
or
dp dt
= CTQ (z-Cx - Du )-ATp.
(12)
Control u minimizing H can be fined from the necessary optimality condition:
— = xTCTQD - zTQD + uT (R + DTQD) + pT B = 0. du v '
Solving this equation we get
u =
(r + DTQD) 1 (DTQz - DTQCx - BTp). (13)
Inserting (13) into (9) h (12) results in the following equations
dx
— = Lx - Mp + Nz, dt
— = -Ux-LTp + Wz, dt
(14)
(15)
-K (t) A(t) - AT (t )K (t) - CT (t )Q(t )C (t)
with terminal condition:
K(tk) = CT (tk)FC(tk);
g(t) is «-dimensional column vector that is solution of linear vector differential equation
dg (t)
where L = A - B (R + DTQD)-1 DTQC;
M = B (R + DTQD BT; N = B (R + DTQD )-1 DTQ; U = CTQC - CQD (R + DTQD DTQC; W = CTQ - CTQD (R + DTQD DT Q;
[ R(t) + DT (t )Q(t) D(t) ] - positively defined symmetric matrix.
The (14) and (15) jointly form nonhomogeneous linear system of differential equations with variables x, p. The solution of this system should satisfy boundary conditions
x(t0) = x0,
p(tK) = CT (tK)FC(tK)x(tK) - CT (tK)Fz(tK). (16)
Let's represent the vector-functionp in the form
p = K(t) x - g(t), (17)
where K(t) is square matrix of « x « size; g(t) - «-dimensional vector.
Make equation for determining K(t) and g(t). For this we differentiate both its parts of (17) by time
dp dK Tjr dx dg
— =-x + K---—.
dt dt dt dt
Taking into account (14), (15) and (17)
dK
-Vx - LT (Kx - g) + Wz = — x +
+KLx - KM (Kx - g) + KNz -
dg dt '
(18)
By making coefficients of terms with x equal to free terms in both parts of (18), we get the equations
dK t
— = -LL k - kl + kmk - и, dt
(19)
dg_ dt
■■(KM - L )g + (KN - W )z. (20)
0
From the boundary condition (16) and formulae (17) follows that matrix K(t) and vector g(t) satisfy terminal conditions:
K (tk ) = CT (tk ) FC (tk ), g (tk ) = CT (tk )Fz (tk ).
(21)
Resulting optimal control u* as function of state vector x is described by the following expression
u =(R + DTQD) (BTg + DTQz - (BTK + DTQC) x). (22)
When D(t) = 0, the equations (19), (20) and (22) coincide with the known equations [9; 13]:
— = -ATK - KA + KBR~lBTK - CTQC , dt
K (tk) = CT (tk)FC (tk), dg = (KBR-lBT - AT )g - CTQz , g (tk ) = CT (tk)Fz (tk), u = R- (BTg - BT Kx).
So, the solving of modified problem ACOR for the linear non-stationary dynamic system (9), (10) let to get the optimal control u* (22) with criteria (11).
Numerical solution of the problem. Consider hypothetic LV moving at altitudes from 0 to 300 meters. Given data are the following parameters of the LV: m, P, I, Cx, Cy \ Cd, L, S, Xd.
Vertical velocity and altitude of the LV are calculated by (5). The "weight" matrixes F, R, Q (11) have the following forms:
F =
fii 0
0 f22
0 0
0 0
0 0
0 0
fas 0
0 f44
Q =
0 0 0
0 0 0
0 q33 0
0 0 0
R = [r].
J =
2 (( (t* )+ f22®2 (tk )+ fas (lp (tk )- Ç (tk ))f
f44AS2 (tk ) + Ij (qaa ( (t) - /ç(t) )2 + r52 (t)) dt.
(23)
If to put (19), (20) with boundary conditions (21) b (22) we get the optimal control in the example task:
u* = KPV1 + Kfflra + K1l + KsAS + u3. (24)
Here the terms of (22) and (23) can be compared in the following way:
[kp1 kw k1 KA»]T =
= -(R +Dt (t)QD(t))-1 (BTK(t) +DT (t)QC(t)), (25)
uQ = (R + DT (t )QD(t) )-1 (BTg (t) + DT (t)Qz (t)).
In example task the input variable u is nozzle deviation angle of the LV sustain engine 5, which is the function of flight time. The preset program 53(t) is calculated together with the preset pitch angle program calculation: 53 = -KasASq. So optimal program (24) can be rewritten in the following form:
8(t) = Kp,(t)p + Ka(t)® + K(t)1 + Kas(t)(AS-ASq). (26)
The first stage of the example task numerical solution is calculation of regulator coefficients by integrating of (19) and (20) in reverse time with boundary conditions (21). The input known data are the elements of matrix C(tk) and vector z(tk).
The second stage is simulation of LV CS with regulator coefficients Kv, Km, K,, KAS calculated in advance.
Several results of numerical solution of the problem are presented below by plots of the LV motion parameters on displacement phase as functions of vehicle flight time (fig. 3-6).
If to state that u(tk) = 0 as state vector x at the end of the displacement phase doesn't affect on the solution, so the criteria (11) has the following scalar form:
In general case "weight" matrixes F, R, Q are unknown and have to be found heuristically. At this technical example elements of "weight" matrixes in criteria (23) were chosen as follows:
/„ = 0.01 (c2/m2), f22 = 5 (c2), ¡33 = 0.05 (1/m2), /44 = 10, q33 = 0.02 (1/m2), r = 6.
It should be pointed out that "weight" coefficient r is the parameter characterizing the "expenditure" of rudder. In the considered task this parameter is the measure of sus-tainer engine chamber deviation. Consequently at low values of parameter r the deviation angle rises. "Weight" coefficient q33 can be interpreted as a "penalty" for deviation of the LV current jet flame trace from preset program value.
Fig. 3. LV altitude of flight h(t) Рис. 3. Высоты полета ЦМ РН
Fig. 7 demonstrates the LV CS regulator coefficients. Fig. 8 shows, that sustainer engine nozzle deviation angle does not exceed 2.5 degrees. After start LV control system commands to deviate the combustion chamber and to hold the current position of LV jet flame on the preset trajectory on starting plane. At 20-30 meters altitude there is some deviation of executed location of LV jet flame from preset location because LV CS provides conditions for further pitch program executing.
Fig. 9 presents the executed and preset positions of LV jet flame trace on the starting plane.
Fig. 4. LV velocity of flight V(t) Fig. 7. LV CS regulator coefficients
Рис. 4. Скорость ЦМ РН Рис. 7. Коэффициенты регулятора СУ РН
Fig. 5. LV horizontal speed Vi (t) Рис. 5. Горизонтальная скорость ЦМ РН
Fig. 8. LV sustainer engine nozzle deviation angle 5(f) Рис. 8. Угол отклонения сопла МД РН
1.5 1 0 5 Я Щ -1 -1 Ш ш -7 m
10 12
Fig. 6. Increments of LV pitch angle 9(f) Рис. 6. Приращение угла тангажа РН
Fig. 9. Preset and executed positions of LV jet flame trace on the starting plane
Рис. 9. Заданное и текущее положения следов струи РН на стартовой плоскости
Fig. 10 show that tracking error does not exceed 0.5 meter. It proves the efficiency of developed algorithms for LV controlled motion displace jet flame from starting complex facilities with tolerated accuracy.
■ю-1-1-1-i-1-
0 50 100 150 200 250 300
h, m
Fig. 10. The difference between preset and executed positions of LV jet flame trace on the starting plane
Рис. 10. Величина рассогласования между заданным
и текущим положениями следов струи РН на стартовой плоскости
Conclusion. Having solved problem ACOR for the system (9), (10) with quadratic terminal integral criteria (11), the optimal control of linear nonstationary dynamic system was developed. Dependency of system output vector on input variable was taken into account in optimal control algorithm (22). The solution is valid for the task of displacement of LV jet flame from the starting complex facilities during the preset displacement program execution. It is demonstrated that the optimal regulator (26) and the calculated coefficients (25) satisfy the specified LV CS requirements. It means that the deviations of executed positions of LV jet flame from preset ones does not exceed 0.5 meter. The CS at the end of the displacement maneuver is holding the LV trajectory close to vertical. The value of sustainer engine chamber deviations do not exceed permissible limits and stays in ±2.5 degree interval.
The problem solved is worth to continue by analysis of three-dimensional path of LV and by the LV jet flame trace preset program tacking accuracy analysis taking into account horizontal wind in atmosphere as a random function of altitude using methods of statistical dynamics [15; 16].
References
1. Al'tshuler A. Sh., Volodin V. D. [Space rocket movement control in the initial flight path, taking into consideration the requirements to reduce the gas-dynamic impact of engine jets on the launching facilities]. Avia-
kosmicheskaya tekhnika i tekhnologiya. 2007, No. 2, P. 3-8 (In Russ.).
2. But A. B., Gonchar A. G., Sokolova T. T. et al. [Development of effective systems of water giving in starting constructions for decrease in gasdynamic loadings].
Kosmonavtika i raketostroenie. 2009, No. 4, P. 11-18 (In Russ.).
3. Stromskiy I. V. Kosmicheskie porty mira [Space ports of the world]. Moscow, Mashinostroenie Publ., 1996, 113 p.
4. Degtyarev A. V., Degtyareva E. A., Novikov A. V. et al. [The results of the analysis of the dynamics of start space rocket "Zenit-3SL" with the afloat sea launch platform]. Aviatsionno-kosmicheskaya tekhnika i tekhnologiya. 2013, No. 9, P. 25-31 (In Russ.).
5. Dyad'kin A. A. [Aerogas dynamics of the spacerocket complex "SEA LAUNCH"]. Kosmicheskaya tekhnika i tekhnologii. 2014, No. 2, P. 14-31 (In Russ.).
6. Legostaev V. P. [Start from the surface of the ocean]. Polet. 1999, No. 2, P. 3-14 (In Russ.).
7. Letov A. M. [Analytical design of controllers I]. Avtomatika i telemekhanika. 1960, No. 4, P. 436-441 (In Russ.).
8. Dmitrievskiy A. A., Lysenko L. N. Prikladnye zad-achi teorii optimal'nogo upravleniya dvizheniem bespi-lotnykh letatel'nykh apparatov [Application-oriented tasks of the theory of optimal control of unmanned aerial vehicles motion]. Moscow, Mashinostroenie Publ., 1978, 328 p.
9. Athans M., Falb P. Optimal Control: An Introduction to the Theory and Its Applications. New York : Dover Publications, 2006. 894 p.
10. Letov A. M. Dinamika poleta i upravlenie [Flight dynamics and control]. Moscow, Nauka Publ., 1969, 360 p.
11. Du W. Dynamic modeling and ascent flight control of Ares-I Crew Launch Vehicle: Dissertation of the Doctor of Philosophy. Ames : Iowa State University. 2010, 167 p. Available at: http://lib.dr.iastate.edu/cgi/ viewcontent.cgi?article=2578&context=etd (accessed 16.04.2017).
12. Lebedev A. A., Chernobrovkin L. S. Dinamika poleta bespilotnykh LA [Flight dynamics unmanned aerial vehicles]. Moscow, Mashinostroenie Publ., 1973, 616 p.
13. Afanas'ev V. N. Teoriya optimal'nogo upravleniya nepreryvnymi dinamicheskimi sistemami [Theory of optimal control of continuous dynamic systems]. Moscow, Moscow St. Univ. Publ., 2011, 168 p.
14. Ivanov V. A., Faldin N. V. Teoriya optimal'nykh sistem avtomaticheskogo upravleniya [Theory of optimal systems of automatic control]. Moscow, Nauka Publ., 1981,336 p.
15. Bobronnikov V. T., Kozorez D. A., Krasil'shchikov M. N. et al. Statisticheskaya dinamika i optimizatsiya upravleniya letatel'nykh apparatov [Statistical dynamics and optimization control aerial vehicle]. Moscow, Al'yans Publ., 2013, 468 p.
16. Bobronnikov V. T., Trifonov M. V. [Technique of the statistical analysis motion of 1st stage LV taking into account random wind loads]. VestnikMAI. 2014, No. 1, P. 33-42 (In Russ.).
Библиографические ссылки
1. Альтшулер А. Ш., Володин В. Д. Управление движением ракеты космического назначения на начальном участке полета с учетом требований по снижению газодинамического воздействия струй двигателей на сооружения стартового комплекса // Авиакосмическая техника и технология. 2007. № 2. C. 3-8.
2. Создание эффективных систем водоподачи в стартовых сооружениях для снижения газодинамических нагрузок / А. Б. Бут [и др.] // Космонавтика и ракетостроение. 2009. № 4. С. 11-18.
3. Стромский И. В. Космические порты мира. М. : Машиностроение, 1996. 113 с.
4. Результаты анализа динамики старта РКН «Зе-нит-38Ь» с находящейся на плаву морской стартовой платформы / А. В. Дегтярев [и др.] // Авиационно-космическая техника и технология. 2013. № 9. С. 25-31.
5. Дядькин А. А. Аэрогазодинамика ракетно-космического комплекса «Морской старт» // Космическая техника и технологии. 2014. № 2. C. 14-31.
6. Легостаев В. П. Старт с поверхности океана // Полет. 1999. № 2. С. 3-14.
7. Летов А. М. Аналитическое конструирование регуляторов I // Автоматика и телемеханика. 1960. № 4. С. 436-441.
8. Дмитриевский А. А., Лысенко Л. Н. Прикладные задачи теории оптимального управления движением беспилотных летательных аппаратов. М. : Машиностроение, 1978. 328 с.
9. Äthans M., Falb P. Optimal Control: An Introduction to the Theory and Its Applications. New York : Dover Publications, 2006. 894 p.
10. Летов А. М. Динамика полета и управление. М. : Наука, 1969. 360 с.
11. Du W. Dynamic modeling and ascent flight control of Ares-I Crew Launch Vehicle : Dissertation of the Doctor of Philosophy. Ames : Iowa State University, 2010. 167 p.
12. Лебедев А. А., Чернобровкин Л. С. Динамика полета беспилотных ЛА. М. : Машиностроение, 1973. 616 с.
13. Афанасьев В. Н. Теория оптимального управления непрерывными динамическими системами. М. : Физический факультет МГУ, 2011. 168 с.
14. Иванов В. А., Фалдин Н. В. Теория оптимальных систем автоматического управления. М. : Наука, 1981. 336 с.
15. Статистическая динамика и оптимизация управления летательных аппаратов / В. Т. Бобронни-ков [и др.] ; под ред. М. Н. Красильщикова, В. В. Малышева. М. : Альянс, 2013. 468 c.
16. Бобронников В. Т., Трифонов М. В. Методика статистического анализа движения первой ступени ракеты-носителя с учетом случайных ветровых нагрузок // Вестник МАИ. 2014. № 1. C. 33-42.
© Altshuler A. Sh., Bobronnikov V. T., Trifonov M. V.., 2017
