,
1.Research Institute of Pilotless Aircraft,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,P.R.China;
2.Key Laboratory of Unmanned Aerial Vehicle Technology,Ministry of Industry and Information Technology,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,P.R.China;
3.College of Computer Science and Technology,Nanjing University of Aeronautics and Astronautics,Nanjing 211106,P.R.China
Abstract:To deal with the high nonlinearities and strong couplings in the transition stage of tail-sitter aircraft,an adaptive gain-scheduling controller is proposed by combining the guardian maps theory and H∞control theory.This method is applied to track the flight-path angle of the transition stage of tail-sitter aircraft,and compared with the linear quadratic regulator(LQR)method based on traditional gain scheduling.Simulation results show that the controller based on the guardian maps theory can autonomously schedule the appropriate control parameters and accomplish the stable transition.Besides,the proposed method shows better tracking performance than the LQR method based on traditional gain scheduling.
Key words:tail-sitter aircraft;transition stage;guardian maps;adaptive gain-scheduling controller
Tail-sitter aircraft is a kind of aircraft that can“sit”on the ground through tail wing[1],which makes it possible to take-off and landing vertically.During the flight,the pitch attitude of aircraft will be adjusted to 90°and the aircraft will transform into the high-speed forward flight mode.Tail-sitter aircraft combines the advantages of both rotor aircraft and fixed-wing aircraft:It can hover,maneuver around and switch to level flight in order to achieve high-speed cruise and rapid transfer in fixed-wing flight mode.Owing to the ability to switch between hovering and level flight mode,tail-sitter aircraft overcomes not only the shortcomings of fixed-wing aircraft limited by the environment but also the weaknesses of rotor aircraft with short battery life and slow speed.Compared with other vertical takeoff and landing aircraft,the main advantage of tailsitter aircraft is its simple mechanical structure.There is no need to increase the mechanism to change its power propulsion direction.In addition,it is also competitive in low failure rate,easy maintenance and so on[2-3].Because of its advantages in task performance and structure,tail-sitter aircraft has attracted considerable attention in scientific research,civil and military fields.
The transition process of tail-sitter aircraft is a unique flight mode.In transition mode,tail-sitter aircraft will switch between horizontal cruise mode and vertical flight mode,which is also a key process for such aircraft to acquire vertical take-off and landing capability and efficient horizontal cruise capability.The difficulty lies in the strong nonlinearity and dramatic dynamic characteristic change in the whole process.It is a common way to apply gain-scheduling control strategy to deal with this kind of issue.
Osborne et al.analyzed the aerodynamic characteristics of tail-sitter aircraft in detail,used gain scheduling technology in the controller design,carried out linearization analysis for the operating points of the aircraft transition flight and designed linear quadratic regulator(LQR)linear feedback control law for each operating point.Then the flight envelope was segmented and the appropriate gain scheduling strategy was designed according to the characteristics of the segmented flight envelope[4].Batailla et al.[5]established the linear parameter varying(LPV)model to describe the transition flight process of tail-sitter aircraft and proposed a ginscheduling controller by interpolating linear controllers for each linear time invariant(LTI)for longitudinal transition control.Besides,the preset transition trajectory is generated by using the neural network generator,and the genetic algorithm with variable environment is used to optimize the neural network instruction generator for transition trajectory optimization,so as to find the most suitable transition trajectory[6].Knoebel et al.[7]used gain scheduling strategy of LQR controller based on the LPV model in simulation,and there are always some errors in tracking instructions in simulation experiments.Kokume et al.[8]mainly studied the controller architecture which can transform horizontal flight to hover state based on fixed-wing unmanned aerial vehicles(UAVs).The dynamic inversion controller andH∞controller were designed for UAV's nonlinear motion equation and the numerical simulation was carried out.Then the author studied the mathematical models of single-propeller fixed-wing aircraft in vertical take-off and vertical transition stages,designed the controllers and validated with LQR andH∞respectively.However,the traditional fixedwing simulation model did not consider the effect of propeller wake on aircraft in transition and vertical states[9].In Ref.[10],a robust formation control method was proposed to achieve the aggressive time-varying formation subject to nonlinear dynamics and uncertainties.In Ref.[11],altitude and position controllers were designed separately based on a nominal controller approach,and robust compensators were adopted to restrain the disturbances.
Since the gain scheduling method is based on each operating point of the nonlinear model,the selection of the operating point has strong subjectivity.A new operating point should be added if the dynamic characteristics of two operating points distinguish from each other.In addition,even if the controller designed at each working point performs well,control system performance during the switching process cannot be guaranteed.In another word,the global stability of the gain-scheduling controller has no theoretical proof.In order to ensure global stability and improve the conservativeness of traditional gain scheduling controller,the guardian maps theory proposed by Saydy was used[12-15].In this paper,a linear parameter varying model of the tail-sitter aircraft is established based on the Jacobian method.The nominalH∞theory is adopted to design the controller at the initial point,and the guardian maps theory is used to analyze the stability boundary of the controller.Then a new controller is designed at the stability boundary,with the stability boundary calculated.The above process is repeated until the whole envelope is covered by the stability boundary.Finally,the proposed controller is used to track the flightpath angle during the tail-sitter aircraft transition process.Simulation results indicate that the controller based on the guardian maps theory can schedule the appropriate control parameters automatically and achieve satisfactory tracking performance,as well as the robustness performance,which is superior over the traditional gain scheduling method based on LQR tracking performance.
The tail-sitter aircraft studied in this paper is shown in Fig.1(a).Similar to Ref.[16],the conventional fixed-wing aircraft model is modified and the tail seat bracket is installed to achieve vertical take-off and landing.The rotors of the tail-sitter are located at the nose of the fuselage,which provides thrust for the vertical takeoff and balances the drag during the cruise.Elevons,which are located at the trailing edge of the fixed-wing,are used to control the pitch angle. The schematic diagram of the tail-sitter flight mode is shown in Fig.1(b).The goal of the flight control system during the vertical-to-horizontal transition is to complete the mode transformation by changing flight speed and pitch angle while keeping the yaw angle and roll angle stable[17].According to the aircraft characteristics during the transition process,the mathematical model is decoupled,and the longitudinal kinematics model is used to design the transition controller.

Fig.1 Flight mode of tail-sitter aircraft
As shown in Fig.2,by decomposing the ground speed vectorVgof the aircraft,we can get the course angle between the horizontal component and the earth-surface inertial reference frame,and the flight path angle between the horizontal component and the ground speed vectorVg.

Fig.2 Flight path angle and course angle
Due to the large angle maneuver of tail-sitter aircraft during the transition process,the small-angle mathematical model will lead to singular phenomena for the large angle maneuver with pitch angle over 45°.Therefore,the horizontal Euler angle(θ)is used to describe the small-angle mathematical model,and the vertical Euler angle(θv)is used to parameterize the attitude matrix when the pitch angle exceeds±45°.According to the motion mode of tail-sitter aircraft during transition flight,the mathematical model is decoupled based on level and non-sideslip flight state,and the equation of flight path angle is added to its longitudinal motion equation system to extract its longitudinal kinematics model.
During the transition between horizontal and vertical flight modes,the vertical Euler angle and the horizontal Euler angle need to be transformed.The transformation formula from horizontal Euler angle to vertical Euler angle is given as[18]

Since the horizontal Euler angle is singular when the pitch angle is 90°in vertical flight mode,the original aircraft mathematical model is transformed into a vertical Euler angle coordinate system by using the transformation formula from the horizontal Euler angle to vertical Euler angle.According to Eq.(1),the longitudinal model of the transition stage can be described by two sets of Euler angles.The longitudinal dynamic equations of the aircraft in the horizontal/vertical Euler angle inertial coordinate system are expressed as

whereTis the engine thrust,Dthe drag,Lthe lift,Jyythe moment of inertia,andMthe pitching moment.In the longitudinal model of aircraft,the flight path angleγis the main variable describing whether the aircraft is flying horizontally during the transition process. Therefore,during the controller design,the flight path angle is selected as the mission evaluation criteria.
Different from nonlinear systems,LPV modeling uses linear systems with varying parameters.There are three approaches to build the LPV model:State transformation,function substitution,and Jacobian linearization.In this paper,the Jacobian linearization method is applied to establish the LPV model of tail-sitter aircraft.By linearizing the nonlinear model at the equilibrium point,a set of linearized models of the tail-sitter aircraft during the transition process can be obtained.Matrix elements are fitted into the polynomial form by means of the linear fitting method.
The range of flight speed of the tailstock aircraft studied in this paper isVa∈(0,10)m/s.AssumingVa=7 m/s before the vertical transition,the flight state and input are decided according to the flight path angle of the aircraft.The states and control inputs at the equilibrium point are obtained by trimming the longitudinal model,then the nonlinear model is linearized to obtain the corresponding state-space model.Because of the high similarity between the linearized system and original non-linear system near the equilibrium point,the resulting LPV model is suitable for the controller design during the transition process.The data listed in Table 1 are the state variables and their inputs at the selected operating points,andAnandBnare the corresponding linearized state-space matrix.
Five representative operating points are selected as the input for linear parameter varying modeling. And the flight-path angleγis taken as the scheduling parameter, and elements in the state-space matrix are fitted as the polynomial functions of the scheduling parameter.In this paper,the maximum polynomial degree is set to be 3,and the detailed polynomials are listed in Appendix A.

where

In order to verify the rationality of the established LPV model,the trajectory angle response curves of the LPV model and the non-linear model under the same step input are compared at a given trajectory angle(γ=89°),as shown in Fig.3.

Table 1 Trim states at different operation points
According to the response results shown in Fig.3,the proposed LPV model of the tail-sitter aircraft matches well with the nonlinear one.

Fig.3 Step response comparison diagram
Basically,guardian maps[17]are scalar-valued maps defined on the set ofn×nreal matrices(ornth-order polynomials)that take nonzero values on the set of“stable”matrices(or polynomials)and vanish on its boundary.The definition of guardian maps is mapping the matrix Rn×nto a complex plane with mappingv,ifdenotes closure of the set,andν(A)=0,thenνguards.The formula is expressed as

The typical stability region consists of the following constraints:Region stability margin constraint,natural frequency constraint,and damping ratio constraint,as shown in Fig.4.
If the guardian maps corresponding to respective regions are given byv1(A),v2(A),…,vm(A),the guardian map for the intersection region can be written as

Synthesize the properties of the guardian maps and make all the eigenvalues in the open left half-plane of the complex plane.Then in the unit circle,we obtain the guardian maps which can remain the characteristic root of the system in the range of the stability margin bigger thanσand the natural frequency smaller thanωn.

According to the synthesis property of guardian maps,the guardian map of sector region shown in Fig.5 can be obtained as


Fig.4 Basic regions

Fig.5 Classical stable region
For the unknown parameterrranging fromrmintormaxin matrixA,ifA(r0)is stable withΩregions andνΩis the guardian map of it,the maximum allowable range of parameterrin the guardian maps region can be obtained by the following procedure:Solve the value ofand the maximum real rootwhich has a minimum value0 in the range of[rmin,r0),-r=rifhas no real root in[rmin,r0),also solve the value of the minimum real rootwhich has a interval upper boundνΩ[A(r) ]=0 in the range of(r0,rmax],r=rmaxifhas no real root in(r0,rmax],which means the maximum stable(or unstable)interval includingr0is.
The improvement ofH∞controller[18]design based on guardian maps theory mainly includes two aspects:One is to make the closed-loop system stable with respect to guardian maps area while satisfying the performance ofH∞;the other is to divide LPV model into several stable and controllable subsystems which can meet the performance requirements in the range of the scheduling parameters,at the same time to complete the stable operation in the whole varying process.
2.2.1 H∞controller design
Considering the system as shown in Fig.6[18]

whereX∈Rnis the state vector of the system,W∈Rqthe external input,U∈Rpthe control inputs,Z∈Rrthe control outputs,Y∈Rmthe measurable output,andA,B*,C*,D*are all real matrices with corresponding dimensions.

Fig.6 Block diagram of the control system
For a full-state feedback controllerU=KX,the closed-loop system is expressed as

In order to stabilize the matrixAclin Eq.(9)and ensure that there exist symmetric matrix when given the necessary and sufficient condition of performance indexthe linear matrix inequality must satisfy

The most used stable-mapping regions in guardian maps are represented by linear matrix inequality(LMI)regions as follows:
(1)When the guardian mapping region is a left half-plane mapping region with stability margin|σ|,whereσis the length of the imaginary axis of the stable region,andL=2σ,M=1,the necessary and sufficient condition for matrixXclto make the LMI region stable is

(2)When the guardian mapping region is a circle region of taking the origin as the center of the circle with the radius less thanωn,the necessary and sufficient condition for matrixXclto make the LMI region stable is

(3)When the guardian mapping region is a sector region with an angle less than 2θ,the necessary and sufficient condition for matrixXclto make the LMI region stable is

If there exists a symmetric positive definite matrixXclsatisfying Eqs.(10)—(13),the closed-loop system hasH∞performance while maintaining the stability of the mapping regionΩ.
2.2.2 Switching control design based on guardian maps
For the improvement of gain scheduling control strategy,we can refer to the guardian maps theory.Firstly,the lower boundaryr0of the LPV model is set to the initial value of the variable parameter model,and the control lawK0satisfying Eqs.(10—13)is calculated[19].
From the equation0,we can calculate the scheduling parameter interval[]of the controllerK0that makes the system to maintain the stability in the guardian map regionΩ.We regard the system model of the parameters in the interval as a subsystem.At this point,the system model of the controllerK0in the parameter range can achieve the control effect which satisfies the performance requirements.Next,is regarded as the initial state of the adjacent second subsystem to solve the controllerK1,then this process is repeated until the scheduling parametersThrough this process,the LPV model can be divided into several subsystems which are stable and controllable and performance requirements satisfied in the range of scheduling parameters;each controller is regarded as the controller within the corresponding scheduling parameters,and the controller is switched on the boundary of each parameter so as to complete the stable operation of the whole changing process.The switching process ofH∞control law based on guardian maps is shown in Fig.7.

Fig.7 Flow chart of robust adaptive H∞controller
In the transitional control of tail-sitter aircraft,considering the problem of speed and flight path angle,the original system state is extended,with the speed and flight path angle errors of the aircraft added into the system state.The expected output value isYrand the error between the actual output and the expected one is expressed as[20]

The integral quantity of the error is

Adding the error integration into the state equation of the LPV model yields

where

The full-state feedback control law is designed for the augmented system,which makesu=Kxaug,and the closed-loop matrix of the system is expressed asAaug+BaugK. The poles of the closed-loop system of tail-sitter aircraft are placed in a stable region. Then,the control parameters,which satisfy the stability and robustness performance index,could be calculated by solving the above LMI conditions for a positive definite symmetric matrixXcl.
The proposed method is applied to design a tracking controller on flight path angle of the tail-sitter aircraft.Before the vertical to horizontal transition,the initial velocity of the aircraft is 7 m/s,and the flight path angle is 90°.After the vertical to horizontal transition,the aircraft will cruise with a constant velocity and zero flight path angle.Then,a horizontal to vertical transition will be presented.In the whole above process,the velocity is assumed to be constant.
With the predefined stability region(|σ|≥2,ωn≤6,ξ≥0.707),the initial controller parameters are calculated by solving Eqs.(10)—(13)with the closed-loop system given by Eq.(17).The resulting initial controllerK0is given by

Then,the stability boundary of the initial controller is calculated by settingresulting inr∈(0,1.13). The stable scheduling parameters of the system are used to calculate the controller again until the boundary of the parameter covers the whole control process.Finally,two sets of controllers are obtained,and the parameter scheduling range that makes the system about the guardian mapping regionΩstable is(-1.336 1,0.53),(0.280 4,2.214).

The controller parameters switch at the midpoint of the intersectional stability regions.We apply the control parameters into five equilibrium points to observe the pole distribution.The pole distribution of the aircraft in open and closed-loop cases is shown in Fig.8.According to Fig.8,the adaptiveH∞controller based on guardian maps can obtain predictable flight quality.
To compare the proposed adaptiveH∞controller based on guardian maps with the traditional gain scheduling LQR controller,the LQR control law is designed at five equilibrium points for the established LPV model,and the gain scheduling control technology is applied to the nonlinear simulation.

Fig.8 Pole distribution comparison diagram
We put forward an adaptiveH∞controller based on guardian maps and the simulation results are compared with the traditional gain scheduling LQR controller. The linear speedVais set to,the initial flight path angle is 90°,the pitch is 90°,the input command is flight path tracking,the slope is-1(°/s),the ramp signal is limited by 90°,the input speed command is 7,and the simulation time is 100 s.The simulation results are shown in Fig.9.

Fig.9 Command tracking comparison diagram
According to the simulation results of Fig.9,when LQR is used to track the flight path angle,there is a certain lag and overshoot in flight path angle tracking and a certain fluctuation of aircraft attitude at the beginning and end of the transition,and airspeed tracking also has some errors.Meanwhile,there is no lag and overshoot in flight path angle tracking when we use adaptiveH∞controller based on guardian maps,and airspeed tracking has no errors either.Moreover,the control by LQR cannot guarantee the flight speed to be maintained in the axial direction,which will produce large components in the other two axes.But this phenomenon will not occur in the guardian maps control.In addition,the computer running time is 30 min when LQR is used for tracking simulation,while theH∞gain scheduling controller designed by guardian maps theory runs for 2 min when tracking simulation is carried out.
Furthermore,the robustness of the proposed controller is verified.Assuming the uncertainties associated with aerodynamic coefficients are within 20%on account of external disturbance,the flight path angle and velocity tracking results of Monte Carlo simulation are shown in Figs.10,11.

Fig.10 Command tracking with 20%uncertainty
According to Fig.11,the throttle angle and elevator deflection at the input end ofH∞gain scheduling controller based on the guardian maps theory is within the acceptable range.And theH∞gain scheduling controller based on the guardian maps theory can effectively control the transition of the tailstock air-sitter and guarantee the satisfactory flight quality.

Fig.11 Control input under 20%uncertainty
A guardian map based adaptive robustH∞controller is proposed for the tail-sitter control during the transition phase,which simultaneously satisfies the stability and robust performance requirements.Then,the proposed method is applied to the longitudinal nonlinear model of tail-sitter aircraft to track the flight path angle at fixed airspeed.Compared with the traditional LQR based gain scheduling method,simulation results indicate that the gain scheduling controller designed by the proposed method achieves satisfactory tracking performance in the transition region.In addition,different from LQR,the proposed method guarantees the global stability of the resulting scheduling controller.
Appendix A


Transactions of Nanjing University of Aeronautics and Astronautics2020年2期