Advanced
Moving Mass Actuated Reentry Vehicle Control Based on Trajectory Linearization
Moving Mass Actuated Reentry Vehicle Control Based on Trajectory Linearization
International Journal of Aeronautical and Space Sciences. 2013. Sep, 14(3): 247-255
Copyright ©2013, The Korean Society for Aeronautical Space Science
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/bync/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : July 09, 2013
  • Accepted : September 11, 2013
  • Published : September 30, 2013
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Xiao-Long Su
Jian-Qiao Yu
jianqiao@bit.edu.cn
Ya-Fei Wang
Lin-lin Wang

Abstract
The flight control of re-entry vehicles poses a challenge to conventional gain-scheduled flight controllers due to the widely spread aerodynamic coefficients. In addition, a wide range of uncertainties in disturbances must be accommodated by the control system. This paper presents the design of a roll channel controller for a non-axisymmetric reentry vehicle model using the trajectory linearization control (TLC) method. The dynamic equations of a moving mass system and roll control model are established using the Lagrange method. Nonlinear tracking and decoupling control by trajectory linearization can be viewed as the ideal gain-scheduling controller designed at every point along the flight trajectory. It provides robust stability and performance at all stages of the flight without adjusting controller gains. It is this “plug-and-play” feature that is highly preferred for developing, testing and routine operating of the re-entry vehicles. Although the controller is designed only for nominal aerodynamic coefficients, excellent performance is verified by simulation for wind disturbances and variations from -30% to +30% of the aerodynamic coefficients.
Keywords
Nomenclature
  • A(t),B(t),Az(t)= state-space system matrices
  • x(t)= state vector
  • u(t)= input vector
  • y(t)= output vector
PPT Slide
Lager Image
  • = nominal state
  • ȳ(t) = nominal output trajectories
  • ū(t) = nominal control
PPT Slide
Lager Image
  • = state error
  • ulc=u-ū= tracking error control input
  • M= mass of maneuvering re-entry vehicle (MaRV) exclusive of moving mass
  • m = mass of moving-mass element
  • V= velocity of MaRV
  • pb= relative position of mass with respect to body coordinate system
  • F= net aerodynamic force on twobody system
  • G= gravitational force on two-body system
  • ρ= air density
  • S= characteristic area
  • α= angle of attack
  • β= sideslip angle
  • Cx0,Cxα2,Cxβ2= resistance coefficients
  • Cy0,Cyα= lift coefficients
  • Czo,Czβ= lateral force coefficients
  • ω = angular velocity
  • ϑ= pitch angle
  • ψ= yaw angle
  • γ= roll angle
PPT Slide
Lager Image
  • = the direction-cosine matrix from body coordinate system to ground coordinate systemg
  • J= moment of inertia of MaRV with respect to body coordinate system
  • Jm= moment of inertia of moving mass with respect to body coordinate system
  • Ma= net aerodtnamic moment about MaRV's center of mass
  • L= characteristic length
PPT Slide
Lager Image
  • = roll moment coefficients
PPT Slide
Lager Image
  • = yaw moment coefficients
PPT Slide
Lager Image
  • = pitch moment coefficients
  • μ= reduced mass parameter
  • ρi,k(t)(k=1, 2) = PD spectrum
  • ζi,ζ= constant damping
  • ωni(t),ω(t) = time-varting bandwidth
  • λ1,λ2= time-varying paramenters
  • Superscripts:
  • g = ground coordinate system
  • b = body coordinate system
1. introduction
Increasing emphasis has been placed on the need for maneuvering re-entry vehicle (MaRV) designs, since future missions for atmospheric re-entry vehicles are facing the problem of a complex environment, short action time, severe weight and volume constraints on actuation and instrumentation. The simplicity of a moving-mass roll control system (MMRCS), combined with its unique ability to provide roll control from within the MaRV's protective shell, make it an attractive alternative to more traditional aerodynamic or thruster-based roll control systems [1] . The purpose of this paper is to present the roll controller using trajectory linearization control (TLC) method that can handle the uncertainties in disturbances and modeling of many modern control problems as exemplified by the controller for a MaRV.
The governing equations of motion of a coupled MaRVmoving mass two-body system are derived using the Lagrange method [2 , 4 , 5] . The mathematical model has a clear physical meaning and is free from force analysis. Classical control theories, such as PID, can barely meet the needs of MMRCS due to the nonlinearity, coupling and time-varying characteristics of the mathematical model. So modern control methods, such as optimum control [1] , quadratic programming [5] , and H control [6] are widely used in designing MMRCS. Although the results of the modern control methods are very attractive, the methods are unsuitable for engineering application. For example, an accurate mathematical model is essential for most modern control methods. However, the uncertainties in disturbances and modeling of an actual system may lead to degradation or failure in the controller.
TLC is an effective nonlinear control method and it has been successfully applied in the control systems of missiles [7 , 8 , 10] , robots [12] and X33 vehicle [13] . The design procedure of TLC consists of the design of two controller subsections. The first one is designed to put the vehicle on the desired trajectory by inverting the nonlinear plant. The second one is a PD-spectrum assignment controller that exponentially stabilizes the linearized tracking error dynamics. This method provides closed-loop global exponential stability without disturbances and gains the maximum robustness against disturbances. The original mathematical model of the nonlinear system is not suitable for analysis, because TLC is based on affine nonlinear systems. So the simplified roll channel dynamic equation derived from the original mathematical model is used to design the controller. Although the controller is designed only for nominal aerodynamic coefficients, excellent performance is verified by simulation for wind disturbances and variations from -30% to +30% of the aerodynamic coefficients.
2. Trajectory Linearization
Suppose the nonlinear system is described by
PPT Slide
Lager Image
where x ( t ) ∈ ℝ n , u ( t ) ∈ ℝ p , y ( t ) ∈ ℝ m , f(x(t)) , g(x(t)) are sufficiently smooth known vector fields of time that are bounded, and have bounded, continuous derivatives up to (n - 1) times. h(x(t) ) is a smooth known function. Let
PPT Slide
Lager Image
ȳ ( t ), ū ( t ) be the nominal state, nominal output trajectories and nominal control satisfying
PPT Slide
Lager Image
Define the state errors and the tracking error control input by
PPT Slide
Lager Image
ulc = u-ū. Then the tracking error dynamics are governed by
PPT Slide
Lager Image
Asymptotic tracking can then be achieved by a 2 Degree-of-Freedom (DOF) controller consisting of: (i) a dynamic inverse I/O mapping of the plant to compute the nominal control function ū for any given nominal output trajectory ȳ ( t ), of which there is a detailed discussed in Ref [11] about the pseudo-inverse and the non-minimum phase case and (ii) a tracking error stabilizing control law ulc to account for modeling simplifications and uncertainties, disturbances and excitation of internal dynamics [9 , 14] . For the unperturbed system of Equation (3), exponential stability is the strongest robustness with respect to all kinds of perturbations, and it guarantees finite gain boundedinput- bounded-output stability. The structure of TLC control is illustrated in Fig. 1 .
Since nominal state
PPT Slide
Lager Image
and nominal input ū ( t ) can be regarded as additional time-varying paramenters of Equation (3), we can rewrite Equation (3) as
PPT Slide
Lager Image
where
PPT Slide
Lager Image
Assumption 1 Let x =0 be an exponentially stable equilibrium point of the nominal system (4), where F : [0, ∞)× D →ℝ n is continuously cifferentiable, D = { e ∈ℝ n |║ e ║˂ r 0 } and the Jaccobian matrix [ ∂F/∂e ] is bounded and Lipshitz on D , uniformly in t . There exists a nominal control law ū and a time-varying feedback control law ulc such that ė = F ( t , e ) is locally esponentially stable.
With the assumption that the tracking errors e are small by performance requirement, the tracking error dynamics can be linearized along the nominal trajectory as
PPT Slide
Lager Image
where
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
Nonlinear tracking system configuration
Assumption 2 The system (5), ( A ( t ), B ( t )) is uniformly completely controllable.
Suppose the linearized error dynamics (5) satisfy Assumption 1 and Assumption 2. Then, there exists a LTV state feedback
PPT Slide
Lager Image
that can exponentially stabilize the system (5) at the origin by assigning to the close-loop system the desired PD spectrum [11] , where
PPT Slide
Lager Image
According to Theorem 3.11 in Ref [15] , nonlinear error dynamic along the nominal trajectory is also exponentially stable at the origin. Thus, the system can be exponentially stabilized along the nominal trajectory.
The detailed design procedure and theory for PD-spectrum assignment is presented in Ref [11] , along with guidelines on the selection of the closed-loop PD-spectrum. According to Theorem 3.1-5.2 in Ref [11] .
PPT Slide
Lager Image
where
PPT Slide
Lager Image
If the subsystem
PPT Slide
Lager Image
is a second-order system, the PD spectrum ρi,k ( t )( k =1, 2) of the i-th second-order system is designed as
PPT Slide
Lager Image
where ζi is the constant damping and ωni ( t ) is the time-varying bandwidth. Then the time-varying coefficients of the second-order system are
PPT Slide
Lager Image
3. Governing Equations of Motion
The realization of the MaRV-moving mass two-body system is shown in Fig. 2 , and it consists of a cone-shaped body. The moving mass is allowed to translate with respect to the MaRV, but is not allowed to rotate with respect to the MaRV.
The system translational dynamics are given by Equation (10).
PPT Slide
Lager Image
The aerodynamic force model is given by Equation (11).
PPT Slide
Lager Image
The kinematic equations of attitude when the system is rolling against the centroid of the shell are given by Equation (12).
PPT Slide
Lager Image
Coordinate-frame definitions
PPT Slide
Lager Image
The system rotational dynamics are given by Equation (13).
PPT Slide
Lager Image
where
PPT Slide
Lager Image
and the reduced-mass parameter is given by
  • μ=mM/(M+m)
The aerodynamic moment model is given by Equation (14).
PPT Slide
Lager Image
where
PPT Slide
Lager Image
Equation (10) and Equation (13) describe the mathematical model of the MaRV-moving mass two-body system. See Ref [4] and [5] for the detailed derivation for the governing equations of motion.
The roll channel dynamic equation is derived according to Ref [4] .
Let l = q =0 and Js = J + μJm for brevity, then
PPT Slide
Lager Image
Where
PPT Slide
Lager Image
From Equation (12) we can get
PPT Slide
Lager Image
Derivate both sides of Equation (16) and substitute it into Equation (15). Then Equation (15) can be rewritten as
PPT Slide
Lager Image
With further consolidation, Equation (17) is rewritten as
PPT Slide
Lager Image
where
PPT Slide
Lager Image
Analysis shows that the system rotational dynamic equation is non-linear, coupled and time-varying. There are also numbers of disturbling moments during the re-entry. However, the sidely used classical PD comtrol theory cannot meet the needs of MMRCS. This paper presents the attitude controller for for the roll channel using TLC.
PPT Slide
Lager Image
Response of roll angle
4. MaRV Controller Design
The controller is based on the roll channel dynamic equation (18) and the desired equation ignoring the disturbance term is rewritten as
PPT Slide
Lager Image
According to the design philosophy of TLC, it is essential to get the nominal control instruction of the system. The nominal control instruction of the system is the control instruction of the vehicle’s roll angle and roll angular velocity, namely
PPT Slide
Lager Image
The nominal position of the moving-mass is
PPT Slide
Lager Image
At the same time, to ensure causality causality,
PPT Slide
Lager Image
are obtained by the following pseudo differentiator
PPT Slide
Lager Image
where ωdiff is the bandwidth of the low pass filter. Various factors should be comprehensively considered in choosing ωdiff , so the low pass filter can eliminate high frequency noise and allow the given control instruction.
Define the state error of system as
PPT Slide
Lager Image
According to the design philosophy of TLC, the linearized matrix of tracking error dynamics e along
PPT Slide
Lager Image
is
PPT Slide
Lager Image
where
PPT Slide
Lager Image
  • a2(t)=(ω2sinγ+ω3cosγ)tanϑ
If the desired closed-loop dynamic behavior is
PPT Slide
Lager Image
Then according to
PPT Slide
Lager Image
The expression of K(t) is
PPT Slide
Lager Image
where
  • k11(t)=0
  • k12(t)=0
  • k21(t)=(-λ1-a1(t))/gf
  • k22(t)=(-λ2-a2(t))/gf
The control input of the system is
PPT Slide
Lager Image
Tracking error of roll angle
PPT Slide
Lager Image
According to the pre-established PD-spectrum theory, the time-varying parameters are
  • λ1=ω2(t)
PPT Slide
Lager Image
5. Simulation
A numerical simulation of the full, nonlinear 6-DOF equations of motion is used to examine the time response of the TLC for the given roll command.
The initial conditions for the simulation are: initial speed V 0 =7000 m / s , initial height h 0 =50 km , initial flight path angle r 0 =10°. The damping and time-varying bandwidth of the linear time-varying regulator are ζ =0.8 and ω ( t )=50 for all the parameters used in the controller. The bandwidth of the low pass filter ωdiff is 10. The lateral position limit of the moving mass is ±0.5m.
The time histories of the roll angle and tracking error are shown Fig. 3 and Fig. 4 . As can be seen from the plot. the roll response is very quick with little overshoot. The maximum peak overshoot is about 0.7 degree, or 1.75% of the 40-degree commanded roll angle. Also, the tracking error is exponentially stabilized as time goes on.
The envelope values of wind speed with a 99% probability are shown in Table 1 according to Ref [16] . Simulations are performed at the same given roll command. Fig. 5 and Fig. 6 show the responses and tracking errors of roll angle with wind disturbances.
The controller is stable when there are wind disturbances.
Envelope values of wind speed with a 99% probability
PPT Slide
Lager Image
Envelope values of wind speed with a 99% probability
PPT Slide
Lager Image
Responses of roll angle with wind disturbances
PPT Slide
Lager Image
Tracking errors of roll angle with wind disturbances
The maximum peak overshoot of all curves is about 0.7 degree, or 1.75% of the 40-degree commanded roll angle. Also, the tracking errors all follow the same trend when exponentially stabilized.
Considering ±30% variations in aerodynamic coefficients and ±10% variations in atmospheric density, the simulations are performed at the same given roll command. Fig. 7 and Fig. 8 show the responses and tracking errors of roll angle in various aerodynamic coefficients.
Obviously, the controller is still stable when there are variations in aerodynamic coefficients. The maximum peak overshoot of all curves is about 0.9 degree, or 2.25% of the 40-degree commanded roll angle. Also, the tracking errors
PPT Slide
Lager Image
Responses of roll angle in ±30% variations
PPT Slide
Lager Image
Tracking errors of roll angle in ±30% variations
all follow the same trend when exponentially stabilized.
6. Conclusion
This paper presented a nonlinear, time-varying controller design for an MaRV using the trajectory linearization method. The nonlinearity, coupling and time-varying characteristics of the MaRV pose great challenges to the controller and TLC provides a satisfactory solution for the MMRCS. The controller structure exhibits considerable inherent robustness and decoupling capability without high actuator activity, providing a useful framework to deal with MaRV problems. Simulation shows that the controller is capable of dealing with different instructions. Although the controller is designed only for nominal aerodynamic coefficients, excellent performance is verified for wind disturbances and ±30% variations of the aerodynamic coefficients. It is this “plug-and-play” feature that is highly preferential for developing, testing and routine operating of the re-entry vehicles.
Future research plans include improving controller performance by: (i) using a nonlinear observer to take advantage of the ignored disturbance term ds for a better output-feedback and (ii) ωi ( t ) is a time-varying coefficient and time variation bandwidth (TVB) method should be taken into account. In particular, (i) should prove effective in overall tracking performance.
References
Petsopoulos Thomas , Regan Frank J. , Barlow Jewel “Moving-mass roll control system for fixed-trim reentry vehicle” Journal of Spacecraft and Rockets 33 (1) 54 - 60    DOI : 10.2514/3.55707
Menon P. K. , Sweriduk G. D. , Ohlmeyer E. J. , Malyevac D. S. “Integrated guidance and control of moving-mass actuated kinetic warheads” Journal of Guidance, control, and Dynamics 27 (1) 118 - 126    DOI : 10.2514/1.9336
Page J. A. , Rogers R. O. “Guidance and control of maneuvering reentry vehicles” Decision and Control including the 16th Symposium on Adaptive Processes and A Special Symposium on Fuzzy Set Theory and Applications, 1977 IEEE Conference on, Vol. 16 659 - 664
Zixing LI , Gaofeng LI “Moving centroid reentry vehicle modeling and active disturbance rejection roll control” Acta Aeronautica et Astronautica Sinica (in Chinese) 33 (11) 2121 - 2129
Vaddi Sesha Sai. 2006 “Moving mass actuated missile control using convex optimization techniques” AIAA Guidance, Navigation, and Control Conference and Exhibit Keyston, Colorado    DOI : 10.2514/6.2006-6575
Menon P. K. , Vaddi S. S. , Ohlmeyer Ernest J. 2006 “Finite-horizon robust integrated guidance-control of a moving-mass actuated kinetic warhead” AIAA Guidance, Navigation, and Control Conference and Exhibit Keyston, Colorado    DOI : 10.2514/6.2006-6787
Mickle M. Chris , Zhu J. Jim 1998 “A nonlinear roll-yaw missile autopilot based on plant inversion and PD-spectral assignment” Decision and Control, Proceedings of the 37th IEEE Conference on, Vol. 4 4679 - 4684
Mickle M. Chris , Zhu J. Jim 1997 “Nonlinear missile planar autopilot design based on pd-spectrum assignment” Decision and Control, Proceedings of the 36th IEEE Conference on, Vol. 4 3914 - 3919
Huang Rui , Mickle M. Christopher , Zhu J. Jim 2003 “Nonlinear time-varying observer design using trajectory linearization” American Control Conference, Proceedings of the 2003 Vol. 6 4772 - 4778
Zhu J. Jim , Mickle M. Christopher “Missile autopilot design using a new linear time-varying control technique” Journal of guidance, control, and dynamics 20 (1) 150 - 157    DOI : 10.2514/2.4009
Zhu J. Jim 1997 “PD-spectral theory for multivariable linear time-varying systems” Decision and Control, Proceedings of the 36th IEEE Conference on, Vol. 4 3908 - 3913
Liu Y. , Wu X. , Jim Zhu, J. , Lew J. 2003 “Omni-directional mobile robot controller design by trajectory linearization.” American Control Conference, Proceedings of the 2003 Vol. 6 3423 - 3428
Zhu J. Jim , Banker D. , Hall Chahes E. 2000 “X-33 ascent flight control design by trajectory linearization-a singular perturbation approach.” AIAA Guidance, Navigation, and Control Conference and Exhibit Denver, Colorado    DOI : 10.2514/6.2000-4159
Zhu Liang “New trajectory linearization control for nonlinear systems undergoing harmonic disturbance” Systems Engineering and Electronics, Journal of 20 (3) 571 - 576
Kahlil Hassan K. 1996 Nonlinear Systems 2nd ed. Prmtice-Hall Upper Saddle River, NI
Johnson Dale Leroy 1993 “Terrestrial environment (climatic) criteria guidelines for use in aerospace vehicle development, 1993 revision.” NASA-TM-4511