Adaptive backstepping sliding mode control for heavy-weight airdrop operations

Abstract. This paper proposes an adaptive backstepping sliding mode flight control method that is compatible with heavyweight cargo airdrop. The goal is to maintain the plane states during cargo extraction process, in the presence of uncertainties of both constant and time-varying types, as well as matched and unmatched types. A backstepping sliding mode flight control law with parameter adaptation is presented based on the plane-cargo dynamics in strict-feedback form. The control approach consists in having an adaptation law that approximates the disturbance and uncertain aerodynamic function, which is separated from the complex nonlinearities. Also, the adaptation algorithm with projection can bound the estimated function. This ensures the robustness of the controller against time-varying disturbance and uncertainty. The convergence performance and robustness property of the control law are proved by the Lyapunov theory. The control effect is evaluated on a transport plane performing a maximum load airdrop task in a number of simulation scenarios.


Introduction
Heavyweight cargo airdrop is a main function of a transport plane, and it has been widely used in modern military tasks and humanitarian aid [1][2][3].The successful release and landing of heavy payloads greatly depends on the operations of the flight control system of the plane, which is imperative to reject sudden and large disturbances caused by the heavy cargo and also to accommodate various system uncertainties [4,5].
Over recent years, some achievements have been reported in developing advanced flight control laws that are compatible with the heavyweight cargo airdrop.Based on the linear system at a given trimming position, [6] investigated the L 1 adaptive control approach subject to system uncertainties, [7] and [8] made a research on the use of active disturbance rejection control technique, and [9] proposed a robust control approach to achieve specified handing qualities.A key problem is that thus linear design controllers may cause unsatisfactory performance in the event that the cargo becomes increasingly heavy.In such an event, the plane dynamics would deviate far from the operating point in cargo extraction, thus leading to a highly nonlinear system.It is interesting to linearize the model at a train of operating points throughout the airdrop process [9], but this is really a tedious work and cannot fundamentally solve the aforementioned problem.Many nonlinear control methods have been proposed to overcome the shortcomings of the linear design controllers.The theoretically established feedback linearization method is the one that is most widely applied [10].
A nonlinear system can be transformed into a constant linear system rather than linear approximations via nonlinear feedback and exact state transformations.However, precise knowledge of the plant model must be known for achieving perfect feedback linearization.It is generally not the case for the airdrop flight controller design, as the complex nonlinear aerodynamic characteristics are very difficult to ascertain and model precisely [5,6,11].Moreover, aerodynamic data obtained from wind tunnel tests always contain a certain degree of uncertainty.Sliding mode control (SMC) is an efficient method to deal with model deficiencies.Based on the feedback linearization model of the airdrop process, a preliminary research on the use of linear SMC techniques to stabilize the plane speed and pitch attitude is introduced in [12], and an iterative SMC strategy that integrates the merits of the global sliding mode and integral sliding mode is proposed in [11].In spite of their interesting ideas, these methods require the upper bounds of the uncertainties to specify the control gains to satisfy the requirement of stability and robustness.However, the complex uncertainties are always unavailable.In this case, the control gains are usually required to be selected large enough to operate correctly under a variety of conditions, which is generally a very conservative strategy [13], and might result in severe chattering phenomena that even damage systems and actuators [14].
Additionally, the aforementioned approaches cannot handle unmatched uncertainties.This also limits their implementations in the design of airdrop flight controllers from a purely practical perspective.Backstepping control is a recursive method for stabilizing systems with both matched and unmatched uncertainties.Combining with the sliding mode method, authors of [15] designed a backstepping SMC law, which solves the unmatched uncertain control problem for cargo airdrop, but their approach still has the conservative problem as mentioned above.As the motivation is the non-requirement of the bounds on the complex nonlinear uncertainties, the idea of employing adaptation algorithms to approximate the uncertainties have attracted much attention.Papers [16][17][18] developed a class of online learning algorithm which is compatible with uncertainties estimation of complex systems.Authors of [19,20] proposed a novel estimation method based on swarm intelligence optimization algorithms, and this method also provide a feasible solution to deal with uncertainties of complex plants.Papers [21,22] introduced neural networks to parameterize the complex uncertainties in designing adaptive backstepping controller, so that the update laws can adapt the network weights.Paper [23] separated the uncertain parameters from the complex nonlinearities and directly estimated the uncertain parameters using adaptation algorithms.The design procedure, as well as the performance analysis of such an approach, is relatively easy when compared with that using the neural networks method.
The main motivation for this work is to present a flight controller design that can accommodate large changes in plane dynamics and reject uncertainties of both constant and time-varying types, as well as matched and unmatched types.The contributions include 1) a flight controller design that inherits the merits of the SMC and backstepping methods, thus solving the unmatched control problem of cargo airdrop; 2) the employment of adaptation techniques to approximate the disturbance and uncertain aerodynamic function, which is separated from the complex nonlinearities; 3) the formation of adaptation laws using the projection operator to bound the estimated function [24,25], which theoretically ensures the robustness of the controller against time-varying uncertainties; and 4) a proof of the convergence performance and robustness property of the control approach based on the Lyapunov theory.

Plane-cargo dynamics
The airdrop process starts as the plane adjusts to a steady flight at the target area.The cargo moves along the rail and is pulled out of the cabin with the help of the extraction system, as shown in Fig. 1.We assumed that the cargo is a particle and the rail coincides with the longitudinal axis of the plane.Thus, the cargo only affects the lateral dynamics of the plane.The disturbances that the cargo imposes on the plane are the contact force , friction force , and additional moment .The aerodynamic forces imposing on the plane include and , which are the lift and drag forces.The thrust and pitch aerodynamic moment are denoted by and , respectively.The pitch flight is characterized by the airspeed , climb angle , attack angle (AOA) , and pitch angle , respectively.
The plane dynamics during airdrop are as follows [11]: where stands for the pitch rate; stands for the mass of the plane; and stands for the pitch moment of inertia.The lift and drag forces are found by: where stands for the dynamic pressure; stands for the wing area; stands for the elevator deflection; * and * represent the drag and lift coefficients, respectively.is given by: where * represent the pitch aerodynamic moment coefficients and stands for the mean aerodynamic chord.
The thrust force is given by: where is the throttle opening ranging from 0 to 100 % and is the maximal thrust.The expressions of , , and derived from the cargo dynamics.As shown in Fig. 2, the cargo dynamics can be obtained as: where stands for the mass of the cargo; and are the -acceleration and -acceleration of the cargo in the wind-axes frame, respectively; stands for the coefficient of cargo to friction; and is the extraction force found by = with denoting the extraction ratio.Using the theorem of composition of particles' acceleration, it follows that the acceleration of the cargo is the sum of its implicated acceleration, Coriolis acceleration, and relative acceleration [11]: where is the distance of point to point .

Fig. 2. Forces analysis of the cargo
Substituting Eqs. ( 12) and ( 13) into Eqs.( 9) and ( 10) yields: = cos + sin + sin − cos + sin + + cos From Fig. 1, the disturbance moment can be obtained as: Remark 1. Eqs. ( 1)-( 3) are the acceleration equation, centripetal acceleration equation, and angular acceleration equation of the plane.Eq. ( 4) is the pitch angle kinematic equation.It is observed from Eqs. ( 1)-( 3) and ( 14)-( 17) that the plane-cargo dynamics form a strongly nonlinear system subject to the coupling of the plane and cargo states.Also, the system contains various uncertainties, such as aerodynamic data perturbation.Readers can refer to [11] for detailed discussions about the model.

Control law and convergence performance
As shown in Fig. 3, the overall control system contains three feedback loops.In the outer loop, an altitude controller is designed by the use of the proportion-differentiation (PD) method.This loop generates a command for the pitch angle controller in the second layer.The controlled variables of the inner loop are [ , ], and the pitch rate command is generated by pitch angle controller.The altitude and airspeed commands are characterized by and , which are defined as those of the trimming position.The design of the pitch angle, pitch rate and speed controllers is the main work.3. Autopilot control architecture with three layers of feedback First, consider the first subsystem in Eq. ( 18), which is rewritten here for convenience: We define: where ( ) = and ( ) is the tracking error of ( ).We can design the control law of Eq. (39) as: it then follows from Eq. ( 49) that: which further leads to: Combining Eqs. ( 60), (61), and (64), the ( ) dynamics is such that: Since (0) ≤ , it follows that ( ) ≤ for every ≥ 0. Note that if + ‖ ‖ < ( ), then + ‖ ‖ < .This completes the proof.Remark 2. Such a projection-based adaptation law can bound the estimated function.This theoretically ensures the robustness of the controller against time-varying disturbance and uncertainty.
Remark 3. It follows from Eq. ( 48) that one can obtain the arbitrary desired tracking performance by increasing the adaptation gain Γ.However, a big gain yields large control power, which can also lead to control signal oscillations.Thus, the adaptation gain needs to be tuned using trial and error.

Simulation results
We simulate a 24,955 kg transport plane airdrop with an 8,000 kg cargo for example.The cargo is initially fixed at the center of gravity of the plane.The parameters of the proposed adaptive backstepping SMC are set as = 1, = 0.5, = 1, = 0.001, and Γ = 0.5 after experimental tuning.The compact sets are conservatively set to = 0.3, Θ = { = ( ) × ∈ × : ∈ [−2,2]}.The outer-loop altitude-hold PD control parameters are selected as = 0.05 and = 0.02.To test the requirements of task performance and flight safety, the evaluation criteria for the heavyweight airdrop can be given as [11]: being the stalling AOA.The performance and robustness of the controller is first tested in the presence of constant aerodynamic uncertainty.As shown in Fig. 4, for all the three scenarios, the altitude and speed are maintained at the trimming position and fully stabilized at about 10 s and 6 s, respectively.The final pitch angle becomes smaller compared to that of the trimming position due to the loss of heavy weight.Note that the observed change in value of the pitch angle is less than 5 deg.The criteria for a successful drop have all been met for the three scenarios.Fig. 5 illustrates that the throttle opening and elevator deflection remain within the limits and they do not occur severe chattering phenomenon.Overall, the pitch up movement of the plane caused by the releasing of the payload is effectively restrained through regulating the throttle and elevator appropriately.(t) = 0.01sin(2t) rad/s, P = 0; (t) = 0, P = 15%sin(2t)[C L0 , C L , C D0 , C D , C m0 , C m , C mq ] T performance evaluation criteria.The elevator deflection and throttle opening change periodically to reject the time-varying uncertainty, and there is no severe chattering, as presented in Fig. 7. Responses comparison of the dropping process with the proposed method and the SMC method [11] (in the presence of 0.01 rad/s pitch rate disturbance and 15 % aerodynamic coefficients uncertainty) We further make a comparison of control performance between the proposed control method and the SMC method in [11], in the presence of 0.01 rad/s pitch rate disturbance and 15 % aerodynamic coefficients uncertainty.The sliding mode controller is recalled from [11] as following: where , , , > 0 are exponentially approaching constants; , are robust terms which need to be chosen larger than the bounds of the system uncertainties; , and are positive constants; , , and are sliding manifolds with = [ , ] ; is the boundary layer thickness; and (⋅) is defined as: (t) = 0.01sin(2t) rad/s, P = 0; (t) = 0, P = 15%sin(2t with sgn(⋅) being the sign function.We set , = 0.1 and = 2 after referring to [11].The outer-loop altitude-hold controller parameters are set as that of the proposed method.Fig. 8 illustrates that the altitude and speed are well maintained within 10 seconds for both of these two controllers, in the presence of uncertainties, and the results meet the task performance evaluation criteria.However, the curves of the throttle opening and speed imply that the SMC method might lead to severe chattering.This is unfavorable for real applications.While using the proposed controller, the undesired chattering is reduced significantly.This is because the projection-based adaptation technique is employed to estimate the uncertainties, and the uncertainties can be compensated by the estimated value in the control input.To summarize, the proposed method not only guarantees tracking performance but also yields a moderate control behavior.

Conclusions
This paper proposes a flight control approach that integrates backstepping sliding mode control with parameter adaptation for the heavyweight cargo airdrop.The controller is capable of stabilizing the plane dynamics in cargo extraction, while being robust to uncertainties of both constant and time-varying types, as well as matched and unmatched types.The method uses projection-based adaptation strategies to achieve robustness against uncertainties, and this overcomes the conservation of the SMC method that relies on the knowledge of the bounds on the complex uncertainties.The convergence and stability properties of the controller are proved by the Lyapunov theory.The control performance is evaluated on a 24,955 kg transport plane airdrop with an 8,000 kg payload.The responses of the plane satisfy the airdrop task performance criteria in the presence of pitch rate disturbance and aerodynamic uncertainty.The application of this research can be used to achieve higher levels of performance and safety in practical airdrop tasks.
The control gains are tuned in a trial and error way, which is a tedious work.The investigation of developing a criterion to specify the control parameters is interesting.This will be one of our future studies.Also, the well-known problems such as input constraints and dead zones have not been discussed in the current work.While a flight control system is actually implemented using this approach, these problems will of course have to be considered.

Fig. 4 .
Fig. 4. The responses of the plane in cargo extraction in the presence of constant uncertainty

Fig. 5 .Fig. 6 .
Fig. 5. Curves of the elevator and throttle in the presence of constant uncertainty

Fig. 7 .
Fig. 7. Curves of the elevator and throttle in the presence of time-varying uncertainty )[C L0 , C L , C D0 , C D , C m0 , C m , C mq ] T