CC BY
99
NAUČNI ČLANCI SCIENTIFIC PAPERS
COMPARISON OF DIFFERENT LATERAL ACCELERATION AUTOPILOTS FOR A SURFACE-TO-SURFACE MISSILE
Ćuk V. Danilo, University of Belgrade, Faculty of Mechanical Engineering, Department for Weapon Systems, Belgrade, Mandić D. Slobodan, Ministry of Defence of the Republic of Serbia, Military Technical Institute, Division for Rocket Dynamics, Belgrade
UDC: 623.462.2 ; 629.7.051.5
FIELD: Mechanical engineering (Rocket technique)
Summary:
This paper presents a comparison of three lateral acceleration autopilots for a surface-to-surface missile: three-loop conventional acceleration autopilot, and gamma-dot and three-loop acceleration autopilot based upon the inverse-dynamic control. The surface-to-surface missile motion is described by nonlinear differential equations whose parameters change rapidly over a very wide range due to variable velocity and altitude. The requirement for the accurate controlling of the missile in such an environment represents a challenge for the autopilot designer. The brief review of the calculation of the autopilot gains is given using the concept of the “point” stability for the linear time-varying system with “frozen” dynamic coefficients. The method of the inverse-dynamic control is presented in the next section for two types of the autopilots: gamma-dot and acceleration autopilot. Both of them require the design of the estimators for the variables used as inputs to the control law. Finally, six-degree-of-freedom simulation results of the missile response to the demanded command on the typical ballistic trajectory are presented. The comparison of three autopilots considers the steady state errors and the sensitivity of the response to the highly variable environment. It was shown that the inverse-dynamic control can be very effective in the controlling of the surface-to-surface missile.
7
Ćuk D. i dr., Comparison of different lateral accelarion autopilots for a surface-to-surface missile, pp. 7-28
VOJNOTEHNIČKI GLASNIK (MILITARY TECHNICAL COURIER), 2011, Vol. LIX, No. 3
Key words: Autopilot, lateral Acceleration Autopilot, Gamma-dot Autopilot, Inverse-Dynamic Control, Guided Surface-to-Surface Missile, ”Six-Degree-Of-Freedom” Model.
Introduction
The aim of this paper is to present the main results of the compari-s on of different autopilot designs for a surface-to-surface missile (SSM). The free flight rocket designed to achieve the range about 50 km was modified to a guided weapon with reduced dispersion. The missile carries the flight computer, inertial measurement unit (IMU) with three rate gyros and three accelerometers, autopilot, and control section. The free flight rolling rocket was redesigned to canard a configuration whose control and guidance sections were stabilized against the rolling motion due to the application of the IMU and the strapdown navigation algorithm.
The purpose of the autopilot is to control the modified SSM during the whole flight, i.e., at low-, medium- and high altitude flight conditions. This is a difficult problem because the aerodynamic controls are very sensitive to the high variation of the air density, Mach number, angle of attack, and other missile dynamic properties needed for the autopilot design.
Typically, the design of the autopilot for the missile is based on the concept of the point stability technique. Most of the autopilots have a fixed structure whose gains are scheduled upon flight conditions such as dynamic pressure and the angle of attack. Therefore, tedious process of generating numerous aerodynamic transfer functions is required. After this step, the "frozen-point” stability is done to develop autopilot gains for the purpose of gain scheduling. This method of the autopilot design was described in many references such as Refs. [1], [2] and [3]. The authors of this paper have developed the computer codes for the complete numerical linearization of six-degree-of-freedom (6-DOF) model [4] and the autopilot design [5]. Different reference trajectories (ballistic, constant manoeuvre, straight line) were employed to generate aerodynamic transfer functions. One of the approaches in preparing the linear missile model based upon 6-DOF concept is shown in Ref. [6].
Many modern control methods with the application to the autopilot designs were published in papers, Refs. [7] - [10]. The method of the robust nonlinear inverse-dynamic control, or feedback linearization, [11] solves the problem of the synthesis of the autopilot for the missile having a highly varying nature. In addition to large variations in aerodynamics, mass, and inertia properties of a SSM, which occur during the boost phase, several constraints were faced by the autopilot designers. The problem of a small manoeuvre capability at high altitudes can be solved partially by allowance of small static instability (the centre of mass is behind
the centre of pressure). The Mach number is changed during the boost phase from the zero value to that greater than 3, when the static instability appears. The dynamic instability and dispersion for a free-flight rocket with extended range was analyzed in details in Ref. [12]. The robust control using the inverse-dynamic method can be applied to the SSM independently of the time instant on the trajectory even in the case of the dynamic instability [13].
In this paper, the results of different autopilot designs are discussed. The three-loop conventional acceleration autopilot is studied first. Using the concept of point stability, all necessary expressions for the autopilot’s gains are rewritten from the published papers. The next section includes the method of the inverse-dynamic control and its application to two types of autopilots: gamma-dot and three-loop lateral acceleration autopilots. The numerical simulation results using a 6-DOF model will be presented to compare the SSM responses to the given commands on the trajectory for all three autopilot designs. Finally, some conclusions are given at the end of this paper.
The autopilot designs, presented in this paper, can be applied to both the missiles with Lambert guidance [14] and those without a thrust terminating mechanism [15].
Three-loop Conventional Acceleration Autopilot
Classical control techniques have dominated missile autopilot designs over the past several decades. Most missile autopilots employ acceleration and rate feedback with the proportional and integral control to stabilize the statically unstable missile and to track the guidance commands (Fig. 1).
The aerodynamic transfer functions of the pitch rate and normal acceleration to the controls deflections in Fig. 1 are, respectively
Aq , ) _ Kq (TgS + l)
Tn2 s 2 + 2gnTns +1
(1)
(2)
where
T 2 _ ______________I
П "» / _ _
n aft -(mwUk - zwMq)
1
(3)
C£>
Ćuk D. i dr., Comparison of different lateral accelarion autopilots for a surface-to-surface missile, pp. 7-28
VOJNOTEHNIČKI GLASNIK (MILITARY TECHNICAL COURIER), 2011, Vol. LIX, No. 3
Kq
znmw zwmn
mc,zw - mwUk
(4)
2grTr
zn
Uk(mw -zwmn)
= mgZn
Uk(mw -z^tnq)
(5)
Uk « V - missile velocity, zw, - derivatives of linear acceleration with respect to the parameter in the subscript, mw, mq, mn- derivatives of angular acceleration with respect to the parameter in the subscript, Ref. [16]. If the reference trajectory is ballistic, Aq = q , Aaz = az, An = n .
Figure 1 - Three-loop conventional acceleration autopilot
It is assumed that the acceleration is at the mass centre (CM), Aa'z = Aaz. If the accelerometer is away from the CM, the transfer function Aaz / Aq(s) has the same form, but the parameters T.f and 2grTr should be modified as it was shown in [1] and [3]. From Fig. 1, the
overall transfer function of the lateral acceleration to demanded acceleration can be developed as:
Aaz
Aazd
(s ) = K
+
T2 s2 + 2gyTys +1
ЦКТ + t2 j s2 +(K_ + 2g0T0
s +1
(7)
where
Kо =-KrK,Kq (1 + KacUk )
(8)
Ćuk D. i dr., Comparison of different lateral accelarion autopilots for a surface-to-surface missile, pp. 7-28
VOJNOTEHNIČKI GLASNIK (MILITARY TECHNICAL COURIER), 2011, Vol. LIX, No. 3
Autopilot Design Using Inverse-Dynamic Control
Theoretical Background
This section follows the description of the inverse-dynamic control given in Ref. [13]. The missile’s nonlinear equations of motion can be separated into a nonlinear homogeneous term plus the term which is the linear function of the control vector
x = f (x) + G (x)u (16)
where x is a n x 1 state vector and u is the control vector of the dimension m x 1; G(x) is n x m matrix of control sensitivity. We assume that the m x 1 output vector y is a linear function of the state vector
y = Hx (17)
The matrix H is selected by the designer. The derivative of the chosen output vector is:
y = Hx = H[f(x) + G(x)u] @ f* (x) + G* (x)u (18)
where f* (x) and G*(x) have dimensions m x 1 and m x m , respectively. Denoting the derivative of y as v
У = v (19)
the inverting control law is obtained
u = G*-1( x)[ v - f* (x)] (20)
The control law includes a model of missile’s homogeneous dynamics in the feedback loop and the inverse of its control effects in the forward loop (Fig. 2). The final step is to design a controller for the system described by Eq. (19).
Figure 2 - Inverse dynamic control
The matrix G* (x) = HG (x) is often not invertible because some elements of y are not linear functions ofu. The repeated differentiation eventually reveals a linear relationship. Suppose that this can be applied to the first element of y, i.e., y1
V = hif (x) (21) where h is the first row of H. The second derivative of y1 is
VI = hjf (x) = h ^[f (x) + G( x)u] (22)
dx
If the necessary linear relationship between u and y is established,
the first rows of f* (x) and G* (x) are redefined as h1 f(x) f(x) and
dx
hi df(x) G(x), respectively.
dx
The second element of y is checked then. If Eq. (22) does not produce a linear relationship, the next derivative of У1 is generated. We denote the vector of suitable derivatives of y by y(d), and f* (x) and G* (x) are defined row by row. Hence,
y(d) = f* (x) + G* (x)u (23)
where G* (x) is invertible for linearly independent controls. If v represents the desired values of y(d), the inverting control is described by Eq. (20). The application of the inverse control can be complex because the evaluation of f* (x) and G* (x) requires that a full, d - differentiable model of the missile be included in the control system. It is worth to say that the inverse-dynamic control can be applied to both nonlinear and linear time-varying systems.
Gamma - Dot Autopilot
The gamma-dot autopilot is used to control the velocity vector turn rate, у. It is useful in the case of SSM with proportional navigation because it provides explicit control of the missile velocity vector turn rate required by the guidance law. If the missile has the pitch and yaw rate gyros for measuring body rates, an estimator is implemented to reconstruct the missile gamma-dots. The error in estimating the gamma-dots depends on the errors in estimating missile’s aerodynamic characteristics. In the ca-
СТз>
Ćuk D. i dr., Comparison of different lateral accelarion autopilots for a surface-to-surface missile, pp. 7-28
VOJNOTEHNIČKI GLASNIK (MILITARY TECHNICAL COURIER), 2011, Vol. LIX, No. 3
se of full strapdown navigation systems, the gamma-dots and all the other parameters required by the autopilot are generated using a navigation algorithm. Since the pitch and yaw channels are identical, only the former will be described.
Using notations in Fig. 3, we can develop a simplified model for missile’s motion. The velocity vector turn rate is given by
(Z + Zc)cosa F -X .
Y =-------c-----+ —----sin a
mV mV
(24)
where Z , Zc Fx, X, m are normal force due to angle-of-attack and controls deflection, thrust, aerodynamic axial force, and missile mass, respectively.
Substituting the following approximations
Z = Zaa, Zc = П M = Maa, Mc = M^n,
Fx - X 1 1 ma
ax = —----, sina^a, cosa^ 1, n =----—a
m mn
into Eq. (24), we get
f = faa
where
fa
dy
da
:(
m
m„
) + -
(25)
(26) (27)
di>
z = ~a = vz z = —
m m
m = ■
Ma V МП
■ = Vmw, mn =—
L,
w n
We select the following state variables:
X1 = у x2 = q
Using the kinematical relationship
a = q-y
after differentiating Eqs. (29) and (30), we get
x = fa ( q - f aa ) + f a a
X2 = Щаa + mvn
(28)
(29)
(30)
(31)
(32)
(33)
Now, we can implement the method of the inverse-dynamic control. The estimated quantity у will be a controlled variable. After an additional differentiating step in Eq. (32) and introducing the new state variables
Ух = X = Y (34)
У2 = X (35)
we get the transformed system of differential equations
II V; (36)
У 2 = V (37)
where
V = fa maa + mn- fa (q - fa- faa + 2fa (q - faa) + faU (38)
The new system described by Eqs. (36) and (37) is a linear, time-invariant system with a new control input v. The deflection of control surfaces may be easily determined from Eq. (38) in terms of the new control input v and the parameters of missile motion.
The control law for the new time-invariant system can be defined to satisfy some design requirements
v = C
Tifc - Ух) + j ifc - Ух )dt ] - C2 У 2 - C1Ух
(39)
where &c is the demanded value of the velocity vector turn rate and yx = у. The coefficients C0,Cx,C2 and the time constant т may be
dr>
Ćuk D. i dr., Comparison of different lateral accelarion autopilots for a surface-to-surface missile, pp. 7-28
VOJNOTEHNIČKI GLASNIK (MILITARY TECHNICAL COURIER), 2011, Vol. LIX, No. 3
defined through the synthesis of the autopilot. The application of Laplace transformation in Eqs. (36), (37) and (39) and solving for y1 = y in terms
of demanded input yc gives the transfer function of the overall time-invariant system
X (s) = _------£(££+1)------------= 3 CoTs + D--------- (40)
Yc s3 + C2 s2 + (C1 + rC0 )s + C0 s3 + C2 s £ + C£ + C0
Equation (38) can be solved for the deflection of controls as
П = K7[T(j'c ~f) + j(YC - f)dt -KA-Kua + Псог(t) (41)
where
C
K =■
famn
K = C2 ~ fa
Г mn
Ka = — (C + C - Krfamn - KyTfamn )
(42)
(43)
(44)
Icor (t)
1
famn
fa-fa(-C2 +3fa)
a
2 fa
famn
q
(45)
The block-diagram of the gamma-dot autopilot is shown in Fig. 4. The deflection of controls depends on the measured value of the pitch
rate q, estimated values of y,a, and the additional correction цоог(t) which is the consequence of the variable velocity vector turn rate derivative fa(t). If fa = fa= 0, the correction to the deflection of controls is equal to zero.
The coefficients C0,C{,C2 are defined by the placement of the closed-loop poles, while the time constant т is selected to compensate for the aperiodic time lag in the autopilot loop.
Estimator Figure 4 - Gamma-dot autopilot
Three-Loop Acceleration Autopilot
The inverse-dynamic control will be applied to the three-loop acceleration autopilot as well. Instead of using the velocity turn rate у as a state
variable, the normal acceleration is introduced
an = -VY = -Vfaa = gaa = X (46)
Since the second state variable remains the pitch rate (x2 = q), the simplified state space model takes the following form:
X = gaa + Saa (47)
x2 = maa + mn (48)
an
a=q-y=q+у
(49)
The controlled variable is normal acceleration, y = an. Taking two
differentiations on the normal acceleration and defining new variables in the transformed system, y = an,y2 = an, we get again:
У = У2 (50)
У 2 = v where
(51)
Ćuk D. i dr., Comparison of different lateral accelarion autopilots for a surface-to-surface missile, pp. 7-28
VOJNOTEHNIČKI GLASNIK (MILITARY TECHNICAL COURIER), 2011, Vol. LIX, No. 3
v=2£al q+у 1+ga+ga
maa + mn-gJ q + у 1-gact
(52)
Hence, we get the time-invariant system and the design of the optimal autopilot can be done easily. Again, after two differentiations the linear function with respect to the deflection of controls in the expression for a new control input v, Eq. (52), is obtained.
The integral control law was applied
v = C0 {(anc - У d - C2У2 - C1У1 (53)
The transfer function of the normal acceleration y1 = an to the demanded value of the normal acceleration anc follows from Eqs. (50), (51) and (53):
a
a
(s) =
Co
s + C2 s + CiS + Cq
(54)
The D.C. gain is equal to unity, and the zero steady-state error occurs. From Eqs. (50) and (51), with the new controlled variable given by Eq. (53), the solution of Eq. (52) for the controls deflection n gives:
n = Kac { (anc - an )dt - Krq - Kaa + Vcor (t) (55)
where
K
ac
cq
Vfamn
(56)
Kr = —(C2 - fa)
m
Ka =---(C1 + ma - Krfamn)
m„
(57)
(58)
Icor (t)
1 г
g m
о a n
ga + ga (Krmn -fa) gafa
a-^&^ q g m
о a n
(59)
The control law of the three-loop acceleration autopilot generated by the inverse-dynamic control is similar to that of the gamma-dot autopilot. In order to obtain a linear time-invariant system, instead of the pitch angle feedback as in the case of the conventional acceleration autopilot, we need the angle-of-attack feedback [see Eq. (55)]. The autopilot is desig-
ned by using the estimator for the angle-of-attack. The compensation of the rapid changes in ga is introduced via the corrected value of the control deflection ncor (t), Eq. (59).
Numerical Examples and Simulation
To determine the response of the missile having different autopilot designs, the six-degree-of-freedom model and numerical simulation were applied. The missile configuration used in the simulation study is an artillery rocket of the range about 50 km. The aspect ratio of the body is l / d = 17,8. The canard section is built into the nose section inside the diameter of the cylindrical body at the distance of l / d = 1,55 from the missile tip. The fin-stabilized configuration has four pop-out fins without cant for generating the rolling motion. The roll autopilot is used to stabilize the canard section with the IMU against the rolling motion. The rocket motor having the total impulse of 33500 dNs imparts the velocity over Mach 3 to the missile to achieve the range over 50 km. The motor has two levels of the thrust 27800/8000 dN with burn times of 0.17/3.80 sec. The missile weight, the mass centre location from the missile tip, the roll and the transverse inertia radius before and after burn is 390/225 kg, 2.65/2.47 m,
0.094/0.098 and 1.2/1.4 m, respectively. The elevation angle of 54.19 deg was chosen to achieve the range of 50 km. The apogee of the trajectory is 19600 m. The aerodynamic data for the 6-DOF model were generated by numerical simulation and wind tunnel tests.
The missile configuration is statically unstable in the time interval t e (3,9;5,1) s when a high value of the Mach number is achieved. Two control points were selected for the autopilot design: t = 10 s (after the burn out time when the missile becomes statically stable again), and t = 60 s (near to the apogee of the trajectory).
Using the "frozen’’ aerodynamic transfer functions parameters for t =10 s and the desired autopilot’s dynamics given by
a>c = Ъ<ап = 3 x5,35 «15 rad/s, ge = 0,6, Te = 0,1 s (60)
the gains of the three-loop conventional autopilot were calculated by Eqs.
(11) - (15):
Kac = 0,012 s/m, Kt = 3,58 1/s, Kr = 0,647 s, kac = 0,91 (61)
The responses of the autopilot to the demanded unit step input for t = 3,5,10, and 15 s are shown in Fig. 5. The time varying transfer functions parameters have important effects on the autopilot’s responses
Ćuk D. i dr., Comparison of different lateral accelarion autopilots for a surface-to-surface missile, pp. 7-28
VOJNOTEHNIČKI GLASNIK (MILITARY TECHNICAL COURIER), 2011, Vol. LIX, No. 3
which are characterized by an increased dynamic errors if the time is different from the point selected for the design.
Figure 5 - Responses of a conventional autopilot designed for t=10 s
The response of the missile to the given command of ±1 g by using the 6-DOF model for the gains in Eq. (61) is shown in Figs. 6 and 7 for the normal acceleration and the deflection of controls and angle-of-at-tack, respectively. A more realistic numerical simulation by the 6-DOF model shows that the satisfactory autopilot response was not obtain even in the vicinity of the chosen control point, t = 10 s. The additional adjustment of the autopilot is required. The change of the gain Kac to the value of Kac = 0.02 gives the response with the reduced steady state error, but
with the increased overshoot, as shown in Fig. 8. The concept of the „frozen" point stability in the design of the three-loop conventional acceleration autopilot cannot give good results for the wide interval of the missile flight.
The identical procedure was repeated for the control point near to the apogee of the trajectory, t = 60 s. Since the natural frequency of the missile is reduced to a>n = 2 rad/s, the desired values of the autopilot dynamics were chosen to be a>c = 6 rad/s, те = 0,5 s, ge = 0,6 . The new values of the autopilot gains were obtained asKac = 0,057 s/m, Kt = 2,49 1/s, Kr = 2,93 s, kac = 0,96. The results of the numerical -simulation by the 6-DOF model in Fig. 9 show again that the additional
adjustments of the gains are useful. Since the manoeuvre capability of the SSM is highly reduced at the apogee, the demanded acceleration was ±0.1 g.
The gamma-dot autopilot was designed using the inverse dynamic control for the same desired dynamic characteristics as those for a conventional autopilot at t=10 s, Eq. (60). Based on these values, the coefficients in the transfer function Eq. (40) were calculated as
т = те = 0.1 s, C0 = 260, C1 = 87, C2 = 16 . Eqs. (42) - (45) were used for
the calculation of the autopilot gains and the correction of the controls deflection due to time varying dynamics needed for the control law, defined by Eq. (41). Fig. 10 illustrates the autopilot response to a step command
±1 deg/s.
In comparison to the conventional autopilot, the gamma-dot autopilot shows better performances because the desired dynamics was achieved without steady state errors. For the same desired dynamics as in the case of the conventional autopilot at t=60 s, the response to the demanded value of the velocity vector turn rate of ±0,1 deg/s is shown in Fig. 11. The comparison of the diagrams in Figs. 9 and 11 proves some advantages of the gamma-dot over the conventional autopilot: faster response and lower values of dynamic errors in tracking the demanded velocity vector turn rate.
Figure 6 - Three-loop conventional acceleration autopilot (t = 10 s, Kac = 0.012)
Ćuk D. i dr., Comparison of different lateral accelarion autopilots for a surface-to-surface missile, pp. 7-28
VOJNOTEHNIČKI GLASNIK (MILITARY TECHNICAL COURIER), 2011, Vol. LIX, No. 3
10
0 5 10 15 20 25
Time, sec
Figure 7 - Three-loop conventional acceleration autopilot (t = 10 s, Kac = 0.012 s/m, 5m = ц)
Figure 8 - Three-loop conventional acceleration autopilot (t = 10 s, Kac = 0.02 s/m)
5
-5
-10
Time, sec
Figure 9 - Three-loop conventional acceleration autopilot (t = 60 s, Kac = 0.057 s/m)
1,5
о
тр
1,0 0,5 0,0 -0,5 -1,0 -1,5
0 5 10 15 20 25
Л f 1 X- 50M
rL 1 n
Vs г
1
1
—°— input
output \P
Time, sec
Figure 10 - Gamma-dot autopilot (design point t = 10 s)
0,15
0,10
0,05
{/)
~В)
0,00
-0,05
-0,10
-0,15-|-----.------------.-------------------------------------------
45 50 55 60 65 70
- input
- output
X-50M
Time, sec
Figure 11 - Gamma-dot autopilot (design point t = 60 s)
The design of the three-loop acceleration autopilot was done using the inverse-dynamic control and having the identical requirements at the control points t = 10 s and t = 60 s as in the examples of the conventional and gamma-dot autopilot. The corresponding responses in Figs. 12 and 13 show that good performances can be achieved with this autopilot. The slower responses were obtained in the comparison with the gamma-dot autopilot because the control law Eq. (55) has only integral of the difference between the demanded signal and the output signal. In spite of time varying missile dynamics, the quality of the transient response does not change with the flight time.
Ćuk D. i dr., Comparison of different lateral accelarion autopilots for a surface-to-surface missile, pp. 7-28
VOJNOTEHNIČKI GLASNIK (MILITARY TECHNICAL COURIER), 2011, Vol. LIX, No. 3
1,5
1,0
3)
^ 0,5
о
2
*> 0,(
<D
О
0 ГО
1 -0,5
ф
га
-1 -1,0 -1,5
0 5 10 15 20
Time, sec
Figure 12 - Three-loop acceleration autopilot designed by the inverse-dynamic control (t = 10 s )
25
0,15
0,10
0,05
0,00
-0,05
-0,10
-0,15
45
50
55 60
Time, sec
65
70
Figure 13 - Three-loop acceleration autopilot designed by the inverse-dynamic control (t = 60 s)
Conclusion
Three types of lateral acceleration autopilots for a surface-to-surface missile were analyzed: three-loop conventional acceleration autopilot, and gamma-dot and three-loop acceleration autopilot based upon the application of the inverse-dynamic control. The design of the first one employs the concept of the point stability for the system with "frozen” aerodynamic transfer functions. After the calculation of the autopilot gains for the full envelope of the missile flight, the additional adjustment of the gains should be done to compensate for the time varying dynamics of
C2D
the SSM. Contrary to this method, the inverse-dynamic control is well suited for the controlling of a missile with the rapidly time-varying dynamics as in the case of a surface-to-surface missile. However, its application requires the design of the estimators for the velocity vector turn rate and the angle-of-attack. The synthesis of such an autopilot begins with the choice of the pre-filter and the compensator for the linear time-invariant system. After this step, the calculation of the controls deflection in terms of the defined new control input and the missile measured or estimated parameters of motion should be done. The gamma-dot and three-loop acceleration autopilot, designed by the inverse-dynamic control, demonstrate excellent performances - fast response and zero steady state errors in a rapidly time varying environment of the SSM. This was proved by the numerical simulation results using a 6-DOF model for the missile flight to the range of 50 km.
There are still some questions to be investigated and solved. It is recommended that research be continued by the design of the required estimators, including the study of the influence of the missile dynamics uncertainties on the autopilot response.
References
[1] Nesline, F.W., Nesline, M.L., How Autopilot Requirements Constrain the Aerodynamic Design of Homing Missiles, AIAA Guidance and Control Conference, TA7, 1983, pp. 716-730.
[2] Gazzina, A., How to Control Missile Airframes: Methodology and Limitation, AGARD-CP-451, Stability and Control of Tactical Missile Systems, 1989, pp. 13.1-13.8.
[3] Ćuk, D., Mandić, S., Pitch Autopilot Design for a Missile with Highly Non-Stationary Dynamical Parameters, Scientific Technical Review, Military Technical Institute, Belgrade, 3/2000.
[4] Ćuk, D., Ćurčin, M., Mandić, S., GMTC_3D: Guided Missile Trajectory Calculation - Three Degree of Freedom Motion Model, Theoretical Manual, Military Technical Institute, Belgrade, 2002.
[5] Ćuk, D., Ćurčin, M., Mandić, S., Autopilot Design, Theoretical Manual, Military Technical Institute, Belgrade, 2004.
[6] Bar-on, J. R., Adams, R. J., Linearization of a Six-Degree-of-Freedom Missile for Autopilot Analysis, J. Guidance, 1/1998, pp. 184-187.
[7] Arrow, A., Williams, D. E., Comparison of Classical and Modern Autopilot Design and Analysis Techniques, J. Guidance, 2/1989, pp. 220-227.
[8] Puri, N. N., Design of an Optimal-Adaptive Digital Autopilot, J. Spacecraft, 7/1970, pp. 1172-1176.
[9] Buschek, H., Full Envelope Missile Autopilot Design Using Gain Scheduled Robust Control, Journal of Guidance, Control, and Dynamics, 1/1999, pp. 115-128.
СЖ>
Ćuk D. i dr., Comparison of different lateral accelarion autopilots for a surface-to-surface missile, pp. 7-28
VOJNOTEHNIČKI GLASNIK (MILITARY TECHNICAL COURIER), 2011, Vol. LIX, No. 3
[10] Zhu, J. J., Mickle, M. C., Missile Autopilot Design using a New Linear Time-Varying Control Technique, Journal of Guidance, Control, and Dynamics, 1/1997, pp. 150-157.
[11] Gratt, H. J., McCowan, W. L., Feedback Linearization Autopilot Design for the Advanced Kinetic Energy Missile Boost Phase, Journal of Guidance, Control, and Dynamics, 5/1995, pp. 945-950.
[12] Ćuk, D., Influence of Range Extension on Dynamic Stability for Artillery Rockets with Wrap Around Fins, Military Technical Courier (Vojnotehnički glasnik), Ministry of Defence of the Republic of Serbia, vol. 55, No. 3, pp. 296-307, ISSN 0042-8469, UDC 623+355/359, Beograd 2007.
[13] Stengel, R. F., Flight Dynamics, Princeton University Press, 2004.
[14] Ćuk, D., Choice of the Rotational Factor of the Thrust Vector for the Ballistic Missile with Lambert Guidance, Military Technical Courier (Vojnotehnički glasnik), Ministry of Defence of the Republic of Serbia, vol. 56, No. 2, pp. 133-146, ISSN 0042-8469, UDC 623+355/359, Beograd 2008.
[15] Gregoriou, G., CEP Calculation for a Rocket with Different Control System, J. Guidance, 3/1988, pp. 193-198.
[16] Garnel, P., Guided Weapon Systems, Pergamon Press, New York, 1980.
POREĐENJE RAZLIČITIH AUTOPILOTA ZA UPRAVLJANJE NORMALNIM UBRZANJEM RAKETE ZEMLJA-ZEMLJA
OBLAST: Mašinstvo (Raketna tehnika)
Sažetak
Ovaj rad prikazuje poređenje tri autopilota za upravljanje normalnim ubrzanjem rakete zemlja-zemlja: konvencionalnog autopilota sa tri povrat-ne sprege, i autopilota za upravljanje ugaonom brzinom vektora brzine i normalnim ubrzanjem koji su zasnovani na primeni inverznog dinamičkog upravljanja. Kretanje rakete zemlja-zemlja opisuje se nelinearnim diferen-cijalnim jednačinama čiji se parametri menjaju u širokom opsegu zbog promene brzine i visine leta. Zahtev za tačnim upravljanjem rakete u ta-kvom okruženju predstavlja izazov za konstruktora autopilota. Dat je sažet pregled proračuna faktora pojačanja konvencionalnog autopilota prime-nom metode lokalne stabilnosti za linearni nestacionarni sistem sa „zamr-znutim“ dinamičkim koeficijentima. Metoda inverznog dinamičkog upravljanja prikazana je u narednom odeljku za dva tipa autopilota: za upravljanje ugaonom brzinom vektora brzine i normalnim ubrzanjem. Oba zahtevaju sintezu estimatora za veličine koje se koriste kao ulazne za zakon upra-vljanja raketom. Na kraju, pomoću modela „6-stepeni slobode kretanja“ daju se rezultati simulacije odgovora rakete na zahtevanu komandu pri le-tu po tipičnoj balističkoj putanji. Poređenjem tri autopilota razmatraju se greške u stacionarnom stanju i osetljivost odgovora na izrazito promenljivo okruženje. Pokazano je da inverzno dinamičko upravljanje može da bude veoma efikasno u upravljanju raketom zemlja-zemlja.
Uvod
Nevođena raketa, koja je projektovana za domet od oko 50 km, mo-difikovana je u vođenu sa smanjenim rasturanjem. Za tako modifikovanu raketu neophodna je primena autopilota za upravljanje normalnim ubrza-njem kako bi se dobio zadovoljavajući odgovor rakete na svim visinama leta. Posle pregleda publikovanih radova iz oblasti raketne tehnike, daje se predmet ovog rada koji obuhvata analizu rezultata odgovora tri tipa autopilota rakete zemlja-zemlja: konvencionalnog autopilota ubrzanja sa tri povratne sprege i dva autopilota koji se zasnivaju na primeni inverznog dinamičkog upravljanja.
Konvencionalni autopilot ubrzanja sa tri povratne sprege
Data je sažeta sinteza autopilota ubrzanja sa tri povratne sprege za statički nestabilnu raketu. Primenom kvazistacionarne metode, sračunava-ju se faktori pojačanja autopilota u zavisnosti od dinamičkih parametara le-ta rakete.
Sinteza autopilota primenom inverznog dinamičkog upravljanja
Posle opisa metode inverznog dinamičkog upravljanja, razmatrana su dva tipa autopilota: autopilot za upravljanje ugaonom brzinom vektora brzine i autopilot za upravljanje normalnim ubrzanjem pakete. Određena je struktura zakona upravljanja i uvedeni određeni estimatori veličina sta-nja. Izvedeni su svi potrebni izrazi za proračun faktora pojačanja i korekci-ju komande zbog nestacionarnosti parametara leta rakete.
Numerički primeri i simulacija
Verifikacija odabranih autopilota izvršena je numeričkom simulacijom odgovora rakete na datu komandu pomoću modela „6 stepeni slobode kretanja“. Pokazano je da primena konvencionalnog autopilota ubrzanja ne daje zadovoljavajuće rezultate zbog povećane statičke greške i preba-čaja koji su funkcije parametara leta rakete. Suprotno tome, autopiloti koji su zasnovani na inverznom dinamičkom upravljanju daju željene odgovo-re nezavisno od promenljivih parametara leta rakete.
Zaključak
Sumirane su dinamičke osobine tri tipa autopilota: konvencionalnog autopilota ubrzanja sa tri povratne sprege, autopilota za upravljanje ugao-nom brzinom tangente na putanju i autopilota za upravljanje normalnim ubrzanjem rakete. Poslednja dva zasnivaju se na primeni savremene metode upravljanja poznate kao inverzno dinamičko upravljanje. Faktori poja-čanja konvencionalnog autopilota dodatno se podešavaju i verifikuju po-moću modela „6 stepeni slobode kretanja“. Suprotno tome, metoda in-verznog dinamičkog upravljanja veoma je pogodna za dinamičke objekte sa jako promenljivim parametrima kakva je i raketa zemlja-zemlja. Zakon
(tT)
Ćuk D. i dr., Comparison of different lateral accelarion autopilots for a surface-to-surface missile, pp. 7-28
VOJNOTEHNIČKI GLASNIK (MILITARY TECHNICAL COURIER), 2011, Vol. LIX, No. 3
upravljanja određuje se u prvom koraku za linearan stacionaran sistem, da bi se u drugom prolazu odredio otklon upravljačkih krila u zavisnosti od novog upravljačkog vektora i promenljivih dinamičkih parametara leta ra-kete. Primena inverznog dinamičkog upravljanja zahteva sintezu određe-nih estimatora veličina stanja rakete, koji bi bili predmet budućih istraživa-nja autopilota za rakete zemlja-zemlja.
Ključne reči: autopilot, autopilot normalnog ubrzanja, autopilot ugaone brzine vektorom brzine, inverzno dinamičko upravljanje, vođena raketa zemlja-zemlja, model „6 stepeni slobode kretanja“.
Datum prijema članka: 26. 01.2011.
Datum dostavljanja ispravki rukopisa: 10. 02. 2011.
Datum konačnog prihvatanja članka za objavljivanje: 12. 02. 2011.
