Next Article in Journal
Numerical Modeling of Thermal Behavior during Lunar Soil Drilling
Next Article in Special Issue
Flight Tracking Control for Helicopter Attitude and Altitude Systems Using Output Feedback Method under Full State Constraints
Previous Article in Journal
Nonintrusive Aerodynamic Shape Optimisation with a POD-DEIM Based Trust Region Method
Previous Article in Special Issue
In-Plane Libration Suppression of a Two-Segment Tethered Towing System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Full State Constrained Flight Tracking Control for Helicopter Systems with Disturbances

1
School of Automation and Information Engineering, Xi’an University of Technology, Xi’an 710048, China
2
School of Sciences, Xi’an Technological University, Xi’an 710021, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Aerospace 2023, 10(5), 471; https://doi.org/10.3390/aerospace10050471
Submission received: 7 March 2023 / Revised: 21 April 2023 / Accepted: 15 May 2023 / Published: 17 May 2023
(This article belongs to the Special Issue Advanced Motion Planning and Control in Aerospace Applications)

Abstract

:
In this paper, a full state-constrained anti-disturbance dynamic surface control method is proposed for six-degree-of-freedom unmanned helicopter systems under full state constraints and disturbances. Firstly, due to the underactuated characteristics of six-degree-of-freedom unmanned helicopter systems, an input–output feedback linearization method is used to transform the complex nonlinear systems into facilitated-control nonlinear ones. Based on the transformed systems, the nonlinear disturbance-observer-based control, backstepping control and Barrier Lyapunov function methods are used to construct the flight controller via uniting the state constraint control and dynamic surface control technologies. Then, Lyapunov stability theory is adopted for analysing the closed-loop tracking error systems, which confirms that the tracking errors are bounded under the proposed flight control scheme. Finally, a simulation in the MATLAB/Simulink environment verifies that the unmanned helicopter system can constrain all states under the action of the designed controller, with good dynamic performance.

1. Introduction

Unmanned aerial vehicles (UAVs) are common flying machines which are widely used in many fields of application, such as the military, agriculture, and others, and have demonstrated strong practical abilities. Indeed, UAVs can be divided into fixed-wing UAVs and unmanned helicopters according to their different flight principles. The unmanned helicopter has numerous advantages, including vertical take-off and landing, low-speed flying, fixed point hovering, and so on, compared to the fixed-wing UAV. Over the past few decades, with the development of composite materials, control theory, navigation and communication systems, unmanned helicopters have entered a new stage of development and received more and more attention [1,2]. However, unmanned helicopter systems have the characteristics of nonlinearity, strong coupling and being under-actuated, which make them vulnerable to the influence of disturbance during flight, which means the controller design becomes a difficult problem to overcome [3]. In recent years, many control methods have been used for the controller design of unmanned helicopter systems. In [4], an input-shaping and model-following control method was proposed to study the adaptive tracking control problem of helicopter hanged flight. In [5], the authors used the linear quadratic regulator control method to study the approach and landing of the helicopter on a limited platform, such as a moving vessel. In [6], the data-driven control problem was discussed for unmanned helicopters without model knowledge, in which the proposed control method was based on model-free adaptive control and the integration and variance approach, and the controller was designed using the input and output data. In [7], a control scheme was proposed by combining a state feedback linear controller and an H controller for reducing the adverse effects of helicopter icing conditions. The disturbance issues were mentioned for unmanned helicopter systems in parts of the above articles, which resulted in dynamic disturbances that can seriously affect the flight safety of helicopters. Therefore, it is necessary to study the advanced anti-disturbance control theory.
The anti-disturbance control is an important research direction in control theory domain. Based on the development of anti-disturbance control theory, many anti-disturbance control methods had been proposed, for example, the disturbance-observer-based control [8], active-disturbance rejection control [9] and H robust control [10], and so on. Among these excellent anti-disturbance control technologies, the disturbance-observer-based control scheme is an effective method to reject the influence of disturbances, in which the disturbance observer is used to estimate the disturbances, whose estimation is applied in the feedforward controller, and ensures that the controller achieves good control performance. Hence, the disturbance-observer-based control scheme is widely used in various control systems. In [11], stochastic systems were investigated under multiple heterogeneous disturbances, which included parts of the known information disturbances, such as being H -norm bounded and induced by white noise, and a random disturbance observer was proposed via combining it with the H control method to achieve disturbance suppression and attenuation. In [12], a shipboard platform system with unknown disturbance was studied, and the anti-disturbance controller was designed by uniting the nonlinear disturbance observer and vectorial backstepping technique. In [13], for the anti-disturbance attitude control problem of spacecraft, the extended state observer and backstepping control methods were employed to suppress disturbances. In [14], the author studied the anti-disturbance control problem of UAVs with unknown disturbance by using a sliding-mode control method. In [15], the disturbance-suppression problem of a space robot arm in a space environment was studied with the sliding-mode control method. Unmanned helicopters fly in the full-disturbances environment. Many disturbances affect the flight quality of helicopters, such as gusts, turbulence and wind shear. For these reasons, the anti-disturbance control is adopted based on the disturbance-observer-based control scheme to solve the robust control issue for helicopter systems. In our past works [16,17], the multiple random disturbances were discussed for helicopter nonlinear systems using composite anti-disturbance control schemes. Therefore, the anti-disturbance control issues are solved for helicopters, which would promote the control performance of the flight controller.
Furthermore, in order to further enhance the flying safety and stability for helicopters, it is necessary to constrain the intermediate control variables when accomplishing satisfactory performance and disturbance suppression ability. The constraint control technology has a wide range of applications in control theory, and there are many areas in practical engineering whose control variables need to be maintained within specified limits. In [18], a full state-constrained prescribed performance control problem was investigated through the adaptive fuzzy control method. In [19], the authors used the backstepping control and Barrier Lyapunov function methods to solve the state constraint problem of nonlinear switching systems. In [20], for the adaptive tracking control problem with asymmetric full state constraints, a control objective with bounded error was achieved by using the unified barrier function and the command filtered backstepping control methods. In [21], a state-constraint control problem of a class of electromagnetic active suspension systems was studied with a cotangent nonlinear state function with the neural network adaptive control method. In [22], a full state-constrained control problem for nonholonomic systems with parameter uncertainty was discussed via using the adaptive control method. In [23], the adaptive tracking control issue was studied for the interconnected nonlinear stochastic system with full state constraints. The above research results provide many new ideas for anti-disturbance control and state-constrained control, which are the new approach to improve the flying dynamic performance of helicopters.
Based on the above studies, in this paper, a flight-tracking controller is designed for six-degree-of-freedom unmanned helicopter systems with full state constraints and disturbances. To enhance the practicability of the proposed controller, the intermediate control variables are constrained for the helicopter while the flight-tracking control task is achieved. Firstly, the six-degrees-of-freedom unmanned helicopter system is transformed by the input–output feedback linearization theory, and the tracking error system of the unmanned helicopter is obtained; secondly, a nonlinear disturbance observer is designed to estimate the system disturbance, and the disturbance estimate is used in the subsequent controller design; Then, in the process of controller design, the nonlinear disturbance observer, backstepping control, dynamic surface and Barrier Lyapunov function methods are combined to construct the full state-constrained flight-tracking control scheme. Under the conditions of full state constraints, the states of the tracking error systems are within the expected bounded range. Finally, the simulation is carried out in the MATLAB/Simulink environment to verify the effectiveness of the proposed flight controller.
This paper is organized as follows: in Section 2, we show the modeling process of helicopter systems; the full state-constrained anti-disturbance flight-tracking controller is constructed for helicopter systems in Section 3; the stability analysis is shown in Section 4 for helicopter tracking error systems; in Section 5, we give the numerical simulation to verify the validity of proposed flight control scheme; and Section 6 is the conclusions.

2. Unmanned Helicopter Models

According to the flight characteristics of an unmanned helicopter, the unmanned helicopter is considered as a six-degree-of-freedom rigid body. By considering the aerodynamic forces and moments during flying, the six-degree-of-freedom system models are established according to the Newton–Euler equations, which are shown as follows [17]:
P ˙ ( t ) = V ( t ) V ˙ ( t ) = g e 3 + 1 m R ( t ) f ( t ) + d 1 ( t ) Ω ˙ ( t ) = H ( t ) w ( t ) J w ˙ ( t ) = w ( t ) × J w ( t ) + τ ( t ) + d 2 ( t )
where P ( t ) = [ x t y ( t ) z ( t ) ] T denote the positions of the helicopter corresponding to the three coordinate axes, V ( t ) = [ v x ( t ) v y ( t ) v z ( t ) ] T denote the velocities of the helicopter, Ω ( t ) = [ ϕ ( t ) θ ( t ) ψ ( t ) ] T are the roll, pitch and yaw angle of the helicopter, respectively, and w ( t ) = [ p ( t ) q ( t ) r ( t ) ] T are the three angular velocity components of the bodyframe coordinate system relative to the terrestrial coordinate system. m and g are the mass of helicopter and gravitational acceleration, J = d i a g { j x x j y y j z z } is the inertia matrix of the helicopter, d 1 ( t ) and d 2 ( t ) denote the interference aerodynamic forces and moments of the helicopter during flight. R ( t ) denotes the rotation matrix of the bodyframe coordinate system relative to the terrestrial coordinate system:
R t = c θ c ψ s θ s ϕ c ψ c ϕ s ψ c ϕ s θ c ψ + s ϕ s ψ c θ s ψ s ϕ s θ s ψ + c ϕ c ψ c ϕ s θ s ψ s ϕ c ψ s θ s ϕ c θ c ϕ c θ
where s ( · ) , c ( · ) denote s i n ( · ) and c o s ( · ) , H ( t ) is attitude matrix:
H t = 1 s ϕ s θ c θ C ϕ s θ c θ 0 c ϕ s ϕ 0 s ϕ c θ c ϕ c θ
w ( t ) × is a cross-product matrix:
w t × = 0 r t q t r t 0 p t q t p t 0
f ( t ) and τ ( t ) are the sums of external forces and moments on the centroid of the helicopter, respectively:
f t = s a c b T M t c a s b T M t T T t c a c b T M t
τ t = z m b t T M t z t T T t + c m b b t z m a t T M t + c m a a t x t T T t Q m t
where T M ( t ) and T T ( t ) are main rotor thrust and tail rotor thrust, respectively, a ( t ) and b ( t ) are the longitudinal and lateral flapping angles, c m a and c m b are the main rotor pitch and roll moment intensity factors, h m = [ x m y m z m ] T and h t = [ x t y t z t ] T are the relative distances from the main rotor and tail rotor to the centroid of the helicopter on the bodyframe coordinate system, respectively. Q m t = c m q T M 1.5 t + d m q is the torque generated by the main rotor, where c m q and d m q are the main rotor torque factors.
The unmanned helicopter is flying with low velocity in the cruising flight phase. The flapping angles a ( t ) and b ( t ) are small, which satisfy s a a t , s b b t , c a = c b 1 , and in this case, we assume: a t T M t 0 , b t T M t T T t . Therefore, the sum of external forces f ( t ) is re-expressed as f t = 0 0 T M t T , and in this paper T M t T T t a t b t T is considered as the system input.
Before designing the controller, in order to simplify the controller design process, the following Assumption and Lemmas are necessary for building the flight controller, which are given as:
Assumption 1
([24]). The disturbance aerodynamic forces and aerodynamic moments satisfy | | d ˙ 1 t | | < D 1 , | | d ˙ 2 t | | < D 2 , where D 1 and D 2 denotes the unknown boundary.
Lemma 1
([25]). For any constant ϵ > 0 and the appropriate dimensional vectors or matrices X and Y , the following inequality holds:
X T Y + Y T X ε X T X + ε 1 Y T Y
Lemma 2
([26]). For any constant k b and real variable z ( t ) , the following inequality holds while z ( t ) < k b :
ln k b 2 k b 2 z 2 ( t ) z 2 ( t ) k b 2 z 2 ( t )

3. Flight-Tracking Controller Design

3.1. System Transformation

In this paper, we design a state-constrained anti-disturbance flight controller for unmanned helicopter systems (1) and implement the asymptotic tracking of the expected tracking signals under the proposed controller. The control outputs are P d t and ψ d t , which denote the helicopter’s expected position trajectories and yaw angle trajectory, respectively. In fact, the control inputs of systems (1) are T M ( t ) , T T ( t ) , a ( t ) and b ( t ) , which are four-dimensional variables. Apparently, the helicopter systems (1) are the six-degree-of-freedom rigid body system models, which are typical underactuated nonlinear systems. In order to achieve the control objective, the input–output feedback linearization method is used to obtain new equivalent systems via the systems (1). The flight-tracking controller is designed under the new systems.
Similar to the analysis process of literature [17], the systems (1) are extended to the following systems by adding two new variables:
P ˙ t = V t V ˙ t = g e 3 1 m R t T M t + d 1 t T ˙ M t = T M 1 t T ˙ M 1 t = T M 2 t Ω ˙ t = H t w t J w ˙ t = w t × Jw t + τ t + d 2 t
The inputs and outputs of new systems (4) are [ T M 2 t τ T t ] T and [ P d T t ψ d t ] T . Then, based on (4), a group of new variables are constructed from the systems (4) via using the input–output feedback linearization methods [17], which are given by:
x 1 t = P t P d t x 2 t = V t P ˙ d t x 3 t = g e 3 e 3 m R t T M t P ¨ d t x 4 t = e 3 m R t T M 1 t e 3 m R t w × t T M t P d t x 5 t = ψ t ψ d t x 6 t = α 1 t w t ψ ˙ d t
where α 1 t = 0 sin ϕ t sec θ t cos ϕ t sec θ t . Moreover, new control inputs are constructed by using the input–output feedback linearization method according to [ T M 2 t τ T t ] T [17], which are given by
u 1 t = T M t m R t β 1 J 1 w t × Jw t 2 T M 1 t m R t w t × e 3 T M t m R t w t × w t × e 3 T M 2 t m R t e 3 T M t m R t β 1 J 1 τ t P d 4 t u 2 ( t ) = α 1 ( t ) J 1 w t × Jw t + c ϕ c θ ϕ ˙ t + s ϕ s θ c θ 2 θ ˙ t q t s ϕ c θ ϕ ˙ t c ϕ s θ c θ 2 θ ˙ t r t + α 1 t J 1 τ t ψ ¨ d t
Indeed, according to [17], x 1 t x 2 t x 3 t x 4 t x 5 t x 6 t T have practical physical significance, which are the helicopter position tracking errors, velocity tracking errors, accelerations generated by controllable join force, force variation rates generated by controllable join force, yaw angle tracking error and yaw angle rate error, respectively.
From (4) and (5), the new transformed systems are obtained, which are given as:
x ˙ 1 t = x 2 t x ˙ 2 t = x 3 t + d 1 t x ˙ 3 t = x 4 t x ˙ 4 t = u 1 t T M t m R t β 1 J 1 d 2 t x ˙ 5 t = x 6 t x ˙ 6 t = u 2 t + α 1 t J 1 d 2 t
where
β 1 = 0 1 0 1 0 0 0 0 0
According to the above process, in this paper, the flight-tracking controller is designed based on the systems (6), u 1 ( t ) and u 2 ( t ) are the control inputs of systems (6) and x 1 t and x 5 t are the outputs of systems (6). In what follows, the tracking controllers u 1 ( t ) and u 2 ( t ) are constructed such that the control outputs x 1 t and x 5 t converge to arbitrary bounded scopes, and the other states are constrained within reasonable domains.

3.2. Nonlinear Disturbance Observer Design

In order to suppress the influence of disturbance aerodynamic forces d 1 ( t ) and aerodynamic moments d 2 ( t ) on the helicopter systems, the nonlinear disturbance observers are used to estimate the disturbances, and substitute the disturbance estimates into the feedforward controller. For the disturbance aerodynamic forces d 1 ( t ) , the disturbance estimator errors d ˜ 1 t = d 1 t d ^ 1 t are defined and the disturbance observer is designed as follows:
d ^ 1 t = δ 1 t + L 1 x 2 t δ ˙ 1 t = L 1 x 3 t + d ^ 1 t
where d ^ 1 t are the disturbance estimates, L 1 = d i a g l 1 , l 1 , l 1 is the disturbance observer gain matrix and δ 1 t denotes the constructive function. The dynamics of the errors in the disturbance estimates are:
d ˜ ˙ 1 ( t ) = d ˙ 1 ( t ) d ˙ ^ 1 ( t ) = d ˙ 1 ( t ) [ δ ˙ 1 ( t ) + L 1 x ˙ 2 ( t ) ] = d ˙ 1 ( t ) [ L 1 x 3 ( t ) L 1 d ^ 1 ( t ) + L 1 x 3 ( t ) + L 1 d 1 ( t ) ] = L 1 d ˜ 1 ( t ) + d ˙ 1 ( t )
Note that (4) and (6) are equivalent systems. The disturbance observer of d 2 ( t ) is designed on the basis of systems (4) as follows:
d ^ 2 t = δ 2 t + L 2 JW t δ ˙ 2 t = L 2 w t × Jw t + τ t + d ^ 2 t
where d ^ 2 t are the disturbance estimates, L 2 = d i a g l 2 , l 2 , l 2 is the disturbance observer gain matrix, and δ 2 t denotes the constructive function. The dynamics of the errors in the disturbance estimates are:
d ˜ ˙ 2 ( t ) = d ˙ 2 ( t ) d ˙ ^ 2 ( t ) = d ˙ 2 ( t ) [ δ ˙ 2 ( t ) + L 2 J w ˙ ( t ) ] = d ˙ 2 ( t ) [ L 2 w ( t ) × J w ( t ) L 2 τ ( t ) L 2 d ^ 2 ( t ) L 2 w ( t ) × J w ( t ) + L 2 τ ( t ) + L 2 d 2 ( t ) = L 2 d ˜ 2 ( t ) + d ˙ 2 ( t )
Based on the designed disturbance observers, in what follows the disturbance estimates are introduced into the process of the state-constrained tracking controller.

3.3. State-Constrained Backstepping Controller Design

In this section, a full state-constrained anti-disturbance controller is designed based on backstepping control, and such that the expected trajectories of position P d ( t ) and yaw angle ψ d ( t ) are tracked under the error state’s constrained condition. The control flow graph of the unmanned helicopter is shown in Figure 1:
The control design mainly includes the design of a nonlinear disturbance observer and a composite anti disturbance tracking controller based on the combination of disturbance compensation and backstepping. The observer obtains disturbance estimates for feedforward compensation, handles the impact of disturbance on the closed-loop control system and derives a reference feedback controller based on the backstepping method to achieve the goal of asymptotically converging the system state. The two work together to achieve the effect of anti-disturbance tracking control. Even if the disturbance degrades to zero, the composite controller can maintain the performance of the reference feedback control and still maintain the asymptotic stability of the system.
Combining the backstepping control scheme with the dynamic surface control method, the first-order filters are introduced to estimate the virtual control inputs, and the problem of “differential explosion” is solved for the nonlinear controller. The specific design process is as follows:
Step 1: Define the tracking error as e 1 t = x 1 t . To ensure accurate tracking, the expected position P d ( t ) with adequate tracking control performance, the maximum helicopter position tracking error values are set as k b 1 i , and e 1 i 0 k b 1 i , for i { 1 , 2 , 3 } .
Define the Barrier Lyapunov function candidate V 1 ( t ) as:
V 1 t = 1 2 i = 1 3 ln k b 1 i 2 k b 1 i 2 e 1 i 2 t
From (6), the derivative of (11) is given as
V ˙ 1 t = e 1 T t C 1 t e ˙ 1 t = e 1 T t C 1 t x 2 t
where C 1 t = d i a g 1 k b 11 2 e 11 2 t 1 k b 12 2 e 12 2 t 1 k b 13 2 e 13 2 t . To ensure the nonnegative definiteness of (12), the virtual control inputs α 1 ( t ) are chosen as:
α 1 t = k 1 A 1 t e 1 t ε 1 2 C 1 t e 1 t
where k 1 1 k b 1 i 2 e 1 i 2 for i = 1 , 2 , 3 is controller gain, ϵ 1 > 0 is the positive constant, and A 1 t = C 1 1 t .
To solve the “differential explosion” problem that arises during the backstepping controller designed process, a first-order filter with time constant λ 1 is constructed, and its outputs are substituted for the virtual control inputs α 1 ( t ) . The initial values of this filter are the same as the initial values of the virtual control inputs. The filter is designed as follows:
λ 1 x ˙ 2 c t + x 2 c t = α 1 t , x 2 c 0 = α 1 0
Define the errors between the filter outputs and the virtual control inputs as γ 2 t = x 2 c t α 1 t . Then, the dynamic of γ 2 t is given by:
γ ˙ 2 t = x ˙ 2 c t α ˙ 1 t = γ 2 t λ 1 α ˙ 1 t
Step 2: Define the errors of the first dynamic surface as e 2 t = x 2 t x 2 c t . We choose the Barrier Lyapunov function candidate V 2 ( t ) as:
V 2 t = 1 2 i = 1 3 ln k b 2 i 2 k b 2 i 2 e 2 i 2 t
where k b 2 i > 0 , for i { 1 , 2 , 3 } , are the maximum tracking errors of e 2 i t . The derivative of (16) along with (6) is shown as:
V ˙ 2 t = e 2 T t C 2 t e ˙ 2 t = e 2 T t C 2 t x 3 t + d 1 t + γ 2 t λ 1
where C 2 t = d i a g 1 k b 21 2 e 21 2 t 1 k b 22 2 e 22 2 t 1 k b 23 2 e 23 2 t .
Similar with the Step 1, in order to guarantee the nonnegative definiteness of V ˙ 2 t , the virtual control inputs α 2 ( t ) are chosen as:
α 2 t = k 2 A 2 t e 2 t A 2 t C 1 t e 1 t ε 2 2 C 2 t e 2 t ε 3 2 C 2 t e 2 t d ^ 1 t γ 2 t λ 1
where k 2 1 k b 2 i 2 e 2 i 2 for i = 1 , 2 , 3 is the controller gain, and ϵ 2 > 0 and ϵ 3 > 0 are positive constants, A 2 t = C 2 1 t .
We design the second first-order filter as:
λ 2 x ˙ 3 c t + x 3 c t = α 2 t , x 3 c 0 = α 2 0
and define the errors between the filter outputs and the virtual control inputs as γ 3 t = x 3 c t α 2 t
The dynamic of γ 3 t is given by:
γ ˙ 3 t = x ˙ 3 c t α ˙ 2 t = γ 3 t λ 2 α ˙ 2 t
Step 3: The error of the second dynamic surface is defined as e 3 t = x 3 t x 3 c t .
The Barrier Lyapunov function candidate V 3 ( t ) is chosen as:
V 3 t = 1 2 i = 1 3 ln k b 3 i 2 k b 3 i 2 e 3 i 2 t
where k b 3 i > 0 , i = { 1 , 2 , 3 } , are the maximum acceleration errors of the helicopter-controllable join force. The derivative of V 3 t along with (6) is given as:
V ˙ 3 t = e 3 T t C 3 t e ˙ 3 t = e 3 T t C 3 t x 4 t + γ 3 t λ 2
where C 3 t = d i a g 1 k b 31 2 e 31 2 t 1 k b 32 2 e 32 2 t 1 k b 33 2 e 33 2 t .
To ensure the nonnegative definiteness of V ˙ 3 t , the virtual control inputs α 3 ( t ) are chosen as:
α 3 t = k 3 A 3 t e 3 t A 3 t C 2 t e 2 t ε 4 2 C 3 t e 3 t γ 3 t λ 2
where k 3 1 k b 3 i 2 e 3 i 2 for i = 1 , 2 , 3 is controller gain, ϵ 4 > 0 is a positive constant and A 3 t = C 3 1 t .
We design the third first-order filter as:
λ 3 x ˙ 4 c t + x 4 c t = α 3 t , x 4 c 0 = α 3 0
and define the errors between the filter outputs and the virtual control inputs as γ 4 t = x 4 c t α 3 t . Then, the dynamic of this error γ 4 t is shown as
γ ˙ 4 t = x ˙ 4 c t α ˙ 3 t = γ 4 t λ 3 α ˙ 3 t
Step 4: The errors of the third dynamic surface are defined as e 4 t = x 4 t x 4 c t . The Barrier Lyapunov function candidate V 4 ( t ) is chosen as:
V 4 t = 1 2 i = 1 3 ln k b 4 i 2 k b 4 i 2 e 4 i 2 t
where k b 4 i , i { 1 , 2 , 3 } are the maximum force variation rate errors from the helicopter-controllable join force. The derivative of V 4 t along with (6) is given by:
V ˙ 4 t = e 4 T t C 4 t e ˙ 4 t = e 4 T t C 4 t u 1 t T M t m R t J 1 d 2 t + γ 4 t λ 3
where C 4 t = d i a g 1 k b 41 2 e 41 2 t 1 k b 42 2 e 42 2 t 1 k b 43 2 e 43 2 t .
To ensure the nonnegative definiteness of V ˙ 4 t , the virtual control inputs u 1 ( t ) are chosen as:
u 1 ( t ) = k 4 A 4 t e 4 t A 4 t C 3 t e 3 t ε 5 2 T M t m R t J 1 T M t m J 1 R T t C 4 t e 4 t + T M t m R t J 1 d ^ 2 t γ 4 t λ 3
where k 4 1 k b 4 i 2 e 4 i 2 , for i = 1 , 2 , 3 , is controller gain and ϵ 5 > 0 is a positive constant, A 4 t = C 4 1 t .
Step 5: Define the error term e 5 t = x 5 t and define the Barrier Lyapunov function candidate V 5 ( t ) as:
V 5 t = 1 2 ln k b 5 2 k b 5 2 e 5 2 t
where k b 5 > 0 is the maximum helicopter yaw angle tracking error. From (6), the derivative of is shown as:
V ˙ 5 t = e 5 t e ˙ 5 t k b 5 2 e 5 2 t = e 5 t x 6 t k b 5 2 e 5 2 t
To ensure the nonnegative definiteness V ˙ 5 t , the virtual control input α 5 ( t ) is chosen as:
α 5 t = k 5 k b 5 2 e 5 2 t e 5 t ε 6 e 5 t k b 5 2 e 5 2 t
where k 5 1 k b 5 2 e 5 2 is controller gain, and ϵ 6 > 0 is a positive constant. The fourth first-order filter is designed as:
λ 4 x ˙ 6 c t + x 6 c t = α 5 t , x 6 c 0 = α 5 0
Then, the error between the filter output and the virtual control input is defined as γ 6 t = x 6 c t α 5 t and the dynamic of γ 6 t is given by:
γ ˙ 6 t = x ˙ 6 c t α ˙ 5 t = γ 6 t λ 4 α ˙ 5 t
Step 6: The error of the fourth dynamic surface is defined as e 6 t = x 6 t x 6 c t and the Barrier Lyapunov function candidate V 6 ( t ) is chosen as:
V 6 t = 1 2 ln k b 6 2 k b 6 2 e 6 2 t
where k b 6 > 0 is the maximum helicopter yaw angle rate error. The derivative of V 6 t is:
V ˙ 6 t = e 6 t e ˙ 6 t k b 6 2 e 6 2 t = e 6 t u 2 t + α 1 t J 1 d 2 t + γ 6 t λ 4 k b 6 2 e 6 2 t
To ensure the nonnegative definiteness of V ˙ 6 t , the control inputs u 2 ( t ) are chosen as:
u 2 t = k 6 k b 6 2 e 6 2 t e 6 t k b 6 2 e 6 2 t k b 5 2 e 5 2 t e 5 t ε 7 2 k b 6 2 e 6 2 t α 1 t J 1 J 1 α 1 T t e 6 t α 1 t J 1 d ^ 2 t γ 6 t λ 4
where k 6 1 k b 6 2 e 6 2 is controller gain, and ϵ 7 > 0 is the positive constant.
Remark 1.
As assumed in many other studies [27], Ω D is a positive invariant set of the virtual control laws α ˙ 1 ( t ) , α ˙ 2 ( t ) , α ˙ 3 ( t ) , α ˙ 5 ( t ) , V : D R is a continuous differentiable function satisfying V ˙ ( α ˙ i ) 0 for i = 1 , 2 , 3 , 5 within Ω, E is the set of all points within Ω satisfying V ˙ ( α ˙ i ) = 0 , and M is the maximal invariant set within E.
Remark 2.
In this paper, an anti-disturbance flight controller is constructed for the helicopter systems via the disturbance observer, backstepping controller and BLF methods. During the helicopter flying, all of the tracking errors of helicopter systems are constrained in limitative scopes, which have many advantages for improving the dynamic performance of helicopter systems. Firstly, the helicopter would be safer under our proposed control scheme. For example, the variable x 3 ( t ) denotes the driving force acceleration of the helicopter, which is constrained in a limitative scope would avoid helicopter overturning. Secondly, the dynamic performance of helicopter would be better using our designed controller. For instance, the variable x 1 ( t ) is the position tracking error, which is constrained in a limitative scope would avoid the big overshoot, and enhance the flying performance.

4. Stability Analysis

Theorem 1.
Consider the tracking error systems of unmanned helicopter (6) with disturbances, design the nonlinear disturbance observers (7) and (9), the virtual control laws (13), (18), (23) and (31), and tracking controller (28) and (36). For given some parameters k b j i > 0 , j = 1 , 2 , 3 and i = 1 , 2 , 3 , 4 , k b h > 0 , h = 5 , 6 , the tracking errors of unmanned helicopter systems (6) converge to the arbitrary bounded scopes with full state constraints. If there exist some parameters ϵ j > 0 , j = 1 , 2 . . . , 13 , k i > 0 , i = 1 , 2 , 6 , such that the following conditions hold: k 1 1 k b 1 i 2 e 1 i 2 , k 2 1 k b 2 i 2 e 2 i 2 , k 3 1 k b 3 i 2 e 3 i 2 , k 4 1 k b 4 i 2 e 4 i 2 , k 5 1 k b 5 2 e 5 2 , k 6 1 k b 6 2 e 6 2 , k ¯ 7 = 2 λ 1 ε 1 1 ε 10 , k ¯ 8 = 2 λ 2 ε 2 1 ε 11 , k ¯ 9 = 2 λ 3 ε 4 1 ε 12 , k ¯ 10 = 2 λ 4 2 ε 6 1 2 ε 13 , k ¯ 11 = 2 l 1 ε 3 1 ε 8 , k ¯ 12 = 2 l 2 ε 5 1 ε 7 1 ε 9 .
Proof. 
According to the above discussion, choose the Lyapunov function candidate V ( t ) as:
V ( t ) = V 1 t + V 2 t + V 3 t + V 4 t + V 5 t + V 6 t + 1 2 γ 2 T t γ 2 t + 1 2 γ 3 T t γ 3 t + 1 2 γ 4 T t γ 4 t + 1 2 γ 6 T t γ 6 t + 1 2 d ˜ 1 T ( t ) d ˜ 1 ( t ) + 1 2 d ˜ 2 T ( t ) d ˜ 2 ( t )
Combining (12), (17), (22), (27) and (30), we have:
V ˙ ( t ) = V ˙ 1 t + V ˙ 2 t + V ˙ 3 t + V ˙ 4 t + V ˙ 5 t + V ˙ 6 t + γ 2 T t γ ˙ 2 t + γ 3 T t γ ˙ 3 t + γ 4 T t γ ˙ 4 t + γ 6 t γ ˙ 6 t + d ˜ 1 T t ˜ ˙ d 1 t + d ˜ 2 T t ˜ ˙ d 2 t = k 1 e 1 T t e 1 t + e 1 T t C 1 t γ 2 t ε 1 2 e 1 T t C 1 t C 1 t e 1 t k 2 e 2 T t e 2 t + e 2 T t C 2 t γ 3 t + e 2 T t C 2 t d ˜ 1 t ε 2 2 e 2 T t C 2 t C 2 t e 2 t ε 3 2 e 2 T t C 2 t C 2 t e 2 t k 3 e 3 T t e 3 t + e 3 T t C 3 t γ 4 t ε 4 2 e 3 T t C 3 t C 3 t e 3 t k 4 e 4 T t e 4 t T M t m e 4 T t C 4 t R t J 1 d ˜ 2 t ε 5 2 T M 2 t m 2 e 4 T t C 4 t R t J 1 J 1 R T t C 4 t e 4 t k 5 e 5 T t e 5 t + e 5 t γ 6 t k b 5 2 e 5 2 t ε 6 e 5 T t e 5 t k b 5 2 e 5 2 t 2 k 6 e 6 T t e 6 t + e 6 t α 1 t J 1 d ˜ 2 t k b 6 2 e 6 2 t ε 7 2 e 6 T t e 6 t k b 6 2 e 6 2 t 2 α 1 t J 1 J 1 α 1 T t + γ 2 T t γ ˙ 2 t + γ 3 T t γ ˙ 3 t + γ 4 T t γ ˙ 4 t + γ 6 T t γ ˙ 6 t l 1 d ˜ 1 T t d ˜ 1 t + d ˜ 1 T t d ˙ 1 t l 2 d ˜ 2 T t d ˜ 2 t + d ˜ 2 T t d ˙ 2 t
According to Lemma 1 and Young’s inequality, there exist ε 1 > 0 , ε 2 > 0 , ε 3 > 0 , ε 4 > 0 , ε 5 > 0 , ε 6 > 0 , ε 7 > 0 , ε 8 > 0 , ε 9 > 0 such that:
e 1 T t C 1 t γ 2 t ε 1 2 e 1 T t C 1 t C 1 T t e 1 t + ε 1 1 2 γ 2 T t γ 2 t e 2 T t C 2 t γ 3 t ε 2 2 e 2 T t C 2 t C 2 T t e 2 t + ε 2 1 2 γ 3 T t γ 3 t e 2 T t C 2 t d ˜ 1 t ε 3 2 e 2 T t C 2 t C 2 T t e 2 t + ε 3 1 2 d ˜ 1 T t d ˜ 1 t e 3 T t C 3 t γ 4 t ε 4 2 e 3 T t C 3 t C 3 T t e 3 t + ε 4 1 2 γ 4 T t γ 4 t T M t m e 4 T t C 4 t R t J 1 d ˜ 2 t ε 5 2 T M 2 t m 2 e 4 T t C 4 t R t J 1 J 1 R T t C 4 t e 4 t + ε 5 1 2 d ˜ 2 T t d ˜ 2 t e 5 t γ 6 t k b 5 2 e 5 2 t ε 6 e 5 T t e 5 t k b 5 2 e 5 2 t 2 + ε 6 1 γ 6 T t γ 6 t e 6 t α 1 t J 1 d ˜ 2 t k b 6 2 e 6 2 t ε 7 2 e 6 T t e 6 t k b 6 2 e 6 2 t 2 α 1 t J 1 J 1 α 1 T t + ε 7 1 2 d ˜ 2 T t d ˜ 2 t d ˜ 1 T t d ˙ 1 t ε 8 2 d ˜ 1 T t d ˜ 1 t + ε 8 1 2 d ˙ 1 T t d ˙ 1 t d ˜ 2 T t d ˙ 2 t ε 9 2 d ˜ 2 T t d ˜ 2 t + ε 9 1 2 d ˙ 2 T t d ˙ 2 t
According to (15), (20), (25) and (33), using the Lemma 1 and Young’s inequality, there exist ε 10 > 0 , ε 11 > 0 , ε 12 > 0 , ε 13 > 0 such that:
γ 2 T t γ ˙ 2 t γ 2 T t γ 2 t λ 1 + ε 10 2 γ 2 T t γ 2 t + ε 10 1 2 α ˙ 1 T t α ˙ 1 t γ 3 T t γ ˙ 3 t γ 3 T t γ 3 t λ 2 + ε 11 2 γ 3 T t γ 3 t + ε 11 1 2 α ˙ 2 T t α ˙ 2 t γ 4 T t γ ˙ 4 t γ 4 T t γ 4 t λ 3 + ε 12 2 γ 4 T t γ 4 t + ε 12 1 2 α ˙ 3 T t α ˙ 3 t γ 6 t γ ˙ 6 t γ 6 T t γ 6 t λ 4 + ε 13 γ 6 T t γ 6 t + ε 13 1 α ˙ 5 T t α ˙ 5 t
Moreover, from Lemma 2, the following inequalities hold:
1 2 i = 1 3 ln k b 1 i 2 k b 1 i 2 e 1 i 2 t e 1 T t C 1 t e 1 t 1 2 i = 1 3 ln k b 2 i 2 k b 2 i 2 e 2 i 2 t e 2 T t C 2 t e 2 t 1 2 i = 1 3 ln k b 3 i 2 k b 3 i 2 e 3 i 2 t e 3 T t C 3 t e 3 t 1 2 i = 1 3 ln k b 4 i 2 k b 4 i 2 e 4 i 2 t e 4 T t C 4 t e 4 t 1 2 ln k b 5 2 k b 5 2 e 5 2 t e 5 T t e 5 t k b 5 2 e 5 2 t 1 2 ln k b 6 2 k b 6 2 e 6 2 t e 6 T t e 6 t k b 6 2 e 6 2 t
Hence, we have:
V ˙ t k 1 e 1 T t e 1 t k 2 e 2 T t e 2 t k 3 e 3 T t e 3 t k 4 e 4 T t e 4 t k 5 e 5 T t e 5 t k 6 e 6 T t e 6 t 1 2 k ¯ 7 γ 2 T t γ 2 t 1 2 k ¯ 8 γ 3 T t γ 3 t 1 2 k ¯ 9 γ 4 T t γ 4 t 1 2 k ¯ 10 γ 6 T t γ 6 t 1 2 k ¯ 11 d ˜ 1 T t d ˜ 1 t 1 2 k ¯ 12 d ˜ 2 T t d ˜ 2 t + D
where D = ε 10 1 2 α ˙ 1 T t α ˙ 1 t + ε 11 1 2 α ˙ 2 T t α ˙ 2 t + ε 12 1 2 α ˙ 3 T t α ˙ 3 t + ε 13 1 α ˙ 5 T t α ˙ 5 t + ε 8 1 2 d ˙ 1 T t d ˙ 1 t + ε 9 2 d ˜ 2 T t d ˜ 2 t + ε 9 1 2 d ˙ 2 T t d ˙ 2 t , for i = 1 , 2 , 3
By defining λ = min k 1 , k 2 , k 3 , k 4 , k 5 , k 6 , k ¯ 7 , k ¯ 8 , k ¯ 9 , k ¯ 10 , k ¯ 11 , k ¯ 12 , (38) can be written as:
V ˙ 2 λ V + D
Thus the unmanned helicopter tracking error systems (6) can converge to the expected bounded range. □

5. Numerical Simulation

In order to verify the state-constrained flight-tracking controller designed in this paper. The parameters presented in Table 1 are selected as the physical parameters of unmanned helicopter and simulated by MATLAB/Simulink [28]:
The initial states of the unmanned helicopter are:
P ( 0 ) = [ 2 2 1.5 ] T m , ψ ( 0 ) = 0.1 r a d
The expected position and yaw of the unmanned helicopter are:
P d = [ 10 c o s ( 0.2 t ) 10 s i n ( 0.2 t ) 5 + 0.1 t ] T m , ψ d = 0.1 s i n ( 0.1 t ) r a d
In order to verify the anti-disturbance ability of proposed control method, the disturbances of the unmanned helicopter are set as:
d 1 = [ 1.5 sin ( 0.5 t ) 0.8 sin ( 0.4 t ) 1.3 sin ( 0.6 t ) ] T , d 2 = [ 0.5 sin ( 0.45 t ) 0.6 sin ( 0.5 t ) 0.6 sin ( 0.4 t ) ] T .
The simulation results of the unmanned helicopter are shown in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10. The position tracking trajectories of the unmanned helicopter are given in Figure 2, which is a 3D figure, in which the helicopter could achieve the position-tracking task with good control performances. Since the unmanned helicopter is under the action of the full state-constrained flight-tracking controller, the velocities of the helicopter, accelerations generated by controllable join force and force variation rates generated by controllable join force are all constrained. Therefore, the positions of all three axes track the desired trajectories asymptotically with a relatively gentle trend, and maintain the tracking states in the subsequent process. The Figure 3 presents the position error curves on the three axes of the unmanned helicopter, which converge to the desired bounded range under the action of state constraint controller. Moreover, the maximum errors of position variables do not exceed 1 × 10 2 m. Obviously, the helicopter has high control accuracy with the constructed anti-disturbance flight controller. The Figure 4 shows the unmanned helicopter attitude angular curves. All these helicopter intermediate control variables are also bounded and kept within a reasonable range which means, though the attitude angles are not the control target variables, these important variables are also in the bounded scopes. Hence, the results from this figure declare the helicopter could safely fly under our full state-constrained control scheme. The control inputs of the unmanned helicopter are shown in Figure 5. To suppress the influence of disturbances, the control inputs display the corresponding waves with the feasible scopes. Figure 6 and Figure 7 represent the disturbance forces and disturbance moments d 1 ( t ) , d 2 ( t ) and their estimated curves.
In order to demonstrate the advantages of our full state-constrained anti-disturbance control scheme, a typical nonlinear control method is applied for designing the flight-tracking controller without using the state-constrained control. The typical nonlinear control method is constructed by combining the nonlinear disturbance observer technology and backstepping control method, which does not adopt the full state-constrained control scheme. The simulation results are shown in Figure 8, Figure 9 and Figure 10, which represent the comparison of the unmanned helicopter position, velocity and attitude angle under the action of the full state-constrained flight-tracking controller and the backstepping controller, respectively, in which the helicopter systems have the same initial conditions. From the figures, it can be seen that the positions, velocities and attitude angles of the unmanned helicopter can be ensured to be constrained within the desired range under the action of the full state-constrained flight-tracking controller to avoid safety problems. Apparently, the curves of the positions, velocity and attitude angles have the smaller overshoot using our flight controller, which implies the helicopter has better control dynamic performance under the proposed flight control scheme.

6. Conclusions

In this paper, the anti-disturbance flight-tracking control problem is studied for unmanned helicopter under full state constraints and disturbances. Firstly, the input–output feedback linearization method is adopted for the unmanned helicopter system to reduce the complexity of controller design. Secondly, two nonlinear disturbance observers are used to estimate these system disturbances. Thirdly, a Barrier Lyapunov function method is combined with the backstepping control technology to construct the full state-constrained anti-disturbance flight-tracking controller. Then, the virtual control laws and their derivatives are estimated by the dynamic surface method to reduce the degree of “differential explosion” caused by nonlinear controller for the high-order system. Finally, the system stability is guaranteed by using the semi-global stability theory, and the effectiveness of designed controller is verified by numerical simulation.

Author Contributions

Conceptualization, Y.L. methodology, Y.L. and Y.H.; software, Y.H.; validation, Y.L. and D.L.; formal analysis, D.L.; writing—original draft preparation, Y.L. and Y.H.; writing—review and editing, D.L.; supervision, H.L.; funding acquisition, Y.L. and D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China grant numbers 62103327 and 12101481, and the Project funded by China Postdoctoral Science Foundation grant numbers 2021MD703879 and 2022T150524.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Raptis, I.A.; Valavanis, K.P. Linear and Nonlinear Control of Small-Scale Unmanned Helicopters; Springer: London, UK, 2011. [Google Scholar]
  2. Ren, B.; Ge, S.S.; Chen, C.; Fua, C.H.; Lee, T. Modeling, Control and Coordination of Helicopter Systems; Springer: New York, NY, USA, 2013. [Google Scholar]
  3. Huang, M.; Xu, D.; He, L.; Zhu, Q.; Wang, L. Development overview and key technologies of high speed hybrid helicopter with single main rotor. J. Aerosp. Power 2021, 36, 1156–1168. [Google Scholar]
  4. Zhang, S.; Ma, Z.; Zhang, H. Research on slung-load swing control method in helicopter suspension flight. Adv. Aeronaut. Sci. Eng. 2022, 13, 57–63. [Google Scholar]
  5. Sebastian, T.; Janusz, N.; Przemyslaw, B. Helicopter Control During Landing on a Moving Confined Platform. IEEE Access 2020, 8, 107315–107325. [Google Scholar]
  6. Xian, B.; Gu, X.; Pan, X. Data driven adaptive robust attitude control for a small size unmanned helicopter. Mech. Syst. Signal Process. 2022, 177, 109205. [Google Scholar] [CrossRef]
  7. Hendra, G.H.; Yoonsoo, K. Flight envelope estimation for helicopters under icing conditions via the zonotopic reachability analysis. Aerosp. Sci. Technol. 2020, 102, 105859. [Google Scholar]
  8. Uwe, N. Design of a disturbance observer based on an already existing observer without disturbance model. AT-Automatisierungstechnik 2022, 70, 134–141. [Google Scholar]
  9. Wang, Y.; Liu, J.; Chen, Z.; Sun, M.; Sun, Q. On the stability and convergence rate analysis for the nonlinear uncertain systems based upon active disturbance rejection control. Int. J. Robust Nonlinear Control 2020, 30, 5728–5750. [Google Scholar] [CrossRef]
  10. Wang, K.; Liu, X.; Jing, Y. Robust finite-time H-infinity congestion control for a class of AQM network systems. Neural Comput. Appl. 2021, 33, 3105–3112. [Google Scholar] [CrossRef]
  11. Dong, L.; Wei, X.; Hu, X.; Zhang, H.; Han, J. Disturbance observer-based elegant anti-disturbance saturation control for a class of stochastic systems. Int. J. Control 2020, 93, 2859–2871. [Google Scholar] [CrossRef]
  12. Hu, X.; Gong, Q. Composite disturbance observer-based backstepping stabilisation for 3-DOF ship-borne platform with unknown disturbances. Ships Offshore Struct. 2022; early Access. [Google Scholar] [CrossRef]
  13. Hou, L.; Sun, H. Anti-disturbance attitude control of flexible spacecraft with quantized states. Aerosp. Sci. Technol. 2020, 99, 105760. [Google Scholar] [CrossRef]
  14. Zhang, X.; Zhu, Z.; Yi, Y. Anti-Disturbance Fault-Tolerant Sliding Mode Control for Systems with Unknown Faults and Disturbances. Electronics 2021, 10, 1487. [Google Scholar] [CrossRef]
  15. Yao, W.; Liu, L.; Guo, Y.; Guo, J. Anti-disturbance control of free floating flexible space manipulator. J. Nanjing Univ. Sci. Technol. 2021, 45, 63–70. [Google Scholar]
  16. Li, Y.; Chen, M.; Ge, S.S.; Li, D. Anti-disturbance control for attitude and altitude systems of the helicopter under random disturbances. Aerosp. Sci. Technol. 2020, 96, 105561. [Google Scholar] [CrossRef]
  17. Li, Y.; Chen, M.; Shi, P.; Li, T. Stochastic anti-disturbance flight control for helicopter systems with switching disturbances under Markovian parameters. IEEE Trans. Aerosp. Electron. Syst. 2023, 1–13. [Google Scholar] [CrossRef]
  18. Shen, L.; Wang, H.; Yue, H. Prescribed Performance Adaptive Fuzzy Control for Affine Nonlinear Systems With State Constraints. IEEE Trans. Fuzzy Syst. 2022, 30, 5351–5360. [Google Scholar] [CrossRef]
  19. Liu, L.; Liu, Y.; Chen, A.; Tong, S.; Chen, C. Integral Barrier Lyapunov function-based adaptive control for switched nonlinear systems. Sci. China Inf. Sci. 2020, 63, 132203. [Google Scholar] [CrossRef]
  20. Gao, Z.; Liu, D.; Qian, M. Decentralised adaptive tracking control for the interconnected nonlinear systems with asymmetric full state dynamic constraints. Int. J. Control 2022, 95, 2840–2853. [Google Scholar] [CrossRef]
  21. Shao, X.; Ye, D. Neuroadaptive deferred full-state constraints control without feasibility conditions for uncertain nonlinear EASSs. J. Frankl. Inst. Eng. Appl. Math. 2022, 359, 2810–2832. [Google Scholar] [CrossRef]
  22. Zhang, Z.; Zhang, S. Adaptive stabilization of state-constrained uncertain nonholonomic system via dynamic surface control. Int. J. Adapt. Control Signal Process. 2021, 35, 1649–1661. [Google Scholar] [CrossRef]
  23. Hua, C.; Meng, R.; Li, K. Hybrid threshold strategy-based adaptive tracking control for nonlinear stochastic systems with full-state constraints. Int. J. Robust Nonlinear Control 2020, 30, 6189–6206. [Google Scholar] [CrossRef]
  24. Xiong, H.; Zhang, H.; Qin, K.; Wang, S.; Li, B. Disturbance observer based fixed-tim integrated attitude-position control for the UAV. Manuf. Autom. 2023, 45, 149–155. [Google Scholar]
  25. Chen, M.; Chen, W. Disturbance-observer-based robust control for time delay uncertain systems. Int. J. Control. Autom. Syst. 2010, 8, 445–453. [Google Scholar] [CrossRef]
  26. Ren, B.; Ge, S.; Tee, K.; Lee, T. Adaptive neural control for output feedback nonlinear systems using a barrier Lyapunov function. IEEE Trans. Neural Netw. 2010, 21, 1339–1345. [Google Scholar]
  27. Chen, M.; Tao, G.; Jiang, B. Dynamic surface control using neural networks for a class of uncertain nonlinear systems with iuput saturation. IEEE Trans. Neural Netw. Learn. Syst. 2015, 26, 2086–2097. [Google Scholar] [CrossRef]
  28. Sun, X.; Fang, Y.; Sun, N. Backstepping-based adaptive attitude and height control of a small-scale unmanned helicopter. Control. Theory Appl. 2012, 29, 381–388. [Google Scholar]
Figure 1. The position-tracking trajectories of the unmanned helicopter.
Figure 1. The position-tracking trajectories of the unmanned helicopter.
Aerospace 10 00471 g001
Figure 2. Three−dimensional image of position tracking trajectories of the unmanned helicopter.
Figure 2. Three−dimensional image of position tracking trajectories of the unmanned helicopter.
Aerospace 10 00471 g002
Figure 3. Position tracking error curves of the unmanned helicopter.
Figure 3. Position tracking error curves of the unmanned helicopter.
Aerospace 10 00471 g003
Figure 4. Attitude angle curves of the unmanned helicopter.
Figure 4. Attitude angle curves of the unmanned helicopter.
Aerospace 10 00471 g004
Figure 5. Control input curves of the unmanned helicopter.
Figure 5. Control input curves of the unmanned helicopter.
Aerospace 10 00471 g005
Figure 6. The curves of disturbance forces and their estimates.
Figure 6. The curves of disturbance forces and their estimates.
Aerospace 10 00471 g006
Figure 7. The curves of disturbance moments and their estimates.
Figure 7. The curves of disturbance moments and their estimates.
Aerospace 10 00471 g007
Figure 8. Position curves of the unmanned helicopter.
Figure 8. Position curves of the unmanned helicopter.
Aerospace 10 00471 g008
Figure 9. Velocity curves of the unmanned helicopter.
Figure 9. Velocity curves of the unmanned helicopter.
Aerospace 10 00471 g009
Figure 10. Attitude angle curves of the unmanned helicopter.
Figure 10. Attitude angle curves of the unmanned helicopter.
Aerospace 10 00471 g010
Table 1. Parameters of unmanned helicopter.
Table 1. Parameters of unmanned helicopter.
Parameter (Unit)Parameter DescriptionParameter (Unit)Parameter Description
m = 8 kgquality of the helicopter J x x = 0.26 kg · m 2 moment of rotation
J y y = 0.35 kg · m 2 moment of rotation J z z = 0.29 kg · m 2 moment of rotation
C m a = 107 ( N · m / rad ) pitch moment intensity factor C m b = 199 ( N · m / rad ) rolling moment intensity factor
C M Q = 0.0044 ( M · N 1 2 ) main rotor torque factor D M Q = 0.6304 ( M · N 1 2 ) main rotor torque factor
x m = 0 mdistance between the center of the main rotor and the x-axis of the helicopter’s center of gravity z m = 0.284 mdistance between the center of the main rotor and the z-axis of the helicopter’s center of gravity
x t = 0.915 mdistance between the center of the tail and the x-axis of the helicopter’s center of gravity z t = 0.104 mdistance between the center of the tail and the z-axis of the helicopter’s center of gravity
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Y.; Huang, Y.; Liu, H.; Li, D. Full State Constrained Flight Tracking Control for Helicopter Systems with Disturbances. Aerospace 2023, 10, 471. https://doi.org/10.3390/aerospace10050471

AMA Style

Li Y, Huang Y, Liu H, Li D. Full State Constrained Flight Tracking Control for Helicopter Systems with Disturbances. Aerospace. 2023; 10(5):471. https://doi.org/10.3390/aerospace10050471

Chicago/Turabian Style

Li, Yankai, Yulong Huang, Han Liu, and Dongping Li. 2023. "Full State Constrained Flight Tracking Control for Helicopter Systems with Disturbances" Aerospace 10, no. 5: 471. https://doi.org/10.3390/aerospace10050471

APA Style

Li, Y., Huang, Y., Liu, H., & Li, D. (2023). Full State Constrained Flight Tracking Control for Helicopter Systems with Disturbances. Aerospace, 10(5), 471. https://doi.org/10.3390/aerospace10050471

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop