Next Article in Journal
Improved Crop Biomass Algorithm with Piecewise Function (iCBA-PF) for Maize Using Multi-Source UAV Data
Next Article in Special Issue
Task Assignment of UAV Swarms Based on Deep Reinforcement Learning
Previous Article in Journal
Blade Twist Effects on Aerodynamic Performance and Noise Reduction in a Multirotor Propeller
Previous Article in Special Issue
Autonomous Maneuver Decision-Making of UCAV with Incomplete Information in Human-Computer Gaming
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Path-Following Control of Small Fixed-Wing UAVs under Wind Disturbance

1
School of Aerospace Engineering, North University of China, Taiyuan 030051, China
2
College of Mechatronics Engineering, North University of China, Taiyuan 030051, China
*
Author to whom correspondence should be addressed.
Drones 2023, 7(4), 253; https://doi.org/10.3390/drones7040253
Submission received: 28 February 2023 / Revised: 23 March 2023 / Accepted: 6 April 2023 / Published: 7 April 2023

Abstract

:
Aiming at the problems of low following accuracy and weak anti-disturbance ability in the three-dimensional path-following control of small fixed-wing Unmanned Aerial Vehicles (UAV), a Globally Stable Integral Sliding Mode Radial Basis Function S-Plane (GSISM+RBF S-Plane) controller is designed. The controller adopts the inner and outer loop mode, the outer loop adopts the Globally Stable Integral Sliding Mode (GSISM) control, and the inner loop adopts the S-Plane control. At the same time, the unknown disturbance in the model is estimated via an RBF neural network. Firstly, the outer loop controller is designed based on the GSISM, and its stability is proved using the Lyapunov theory. Then, the S-Plane controller is designed for the instruction signal of the inner loop. Considering the complexity of the derivation in the S-Plane controller, a second-order differentiator is introduced. Finally, considering the problem of external wind disturbance, the controller is modeled, studied, and processed in order to better reflect the impact of real external wind on UAV path following. Finally, the Globally Stable Sliding Mode (GSSM) control and Globally Stable Integral Sliding Mode S-Plane (GSISM S-Plane) control are used for a comparative experiment. The simulation results show that the designed GSISM+RBF S-Plane controller can accurately track the ideal path compared with the GSSM and GSISM S-Plane controller, and it has good control performance and anti-disturbance performance.

1. Introduction

With the development of electronic information technology, communication, and intelligent control, UAVs are becoming increasingly intelligent. UAVs play an important role in the military and civil fields, especially in the service fields [1]. Viewed from the angle of technology, UAVs can be mainly classified into fixed-wing UAVs, rotor UAVs, and unmanned helicopters. Fixed-wing UAVs have a large load, long range, and fast flight speed. They have been widely used in communication relays, enemy attacks, swarm operations, reconnaissance and surveillance, forest fire control, disaster relief and reduction, sightseeing tourism, agricultural plant protection, etc. [2]. Among the various subsystems of UAVs, the flight control system is an important component in determining whether the UAV can successfully complete a scheduled task. Trajectory tracking and path following are important components of the UAV flight control system. Trajectory tracking requires the UAV to reach a predetermined place within a predetermined time according to the ideal trajectory, while path following has no specific requirements for the arrival time of a UAV following a path. The research in reference [3] also shows that path following is more valuable than trajectory tracking in small fixed-wing UAVs.
UAV path following involves geometric methods and control methods. Geometric methods mainly include the pure pursuit method, line of sight method [4], vector field method [5], nonlinear guidance rate [6], L1 guidance rate, and a combination of the line-of-sight method and pure pursuit method. The control method mainly includes PID control [7], model predictive control [8,9], sliding mode control [10,11,12], LQR control [13] backstepping control [14,15], linear quadratic regulator control, etc. The control principle of the pure pursuit method, line of sight method, and nonlinear guidance rate [16,17] aim to provide a virtual point on the ideal path. The UAV aligns its course with the virtual target point by adjusting the yaw angle so as to follow the ideal path. The main idea of the vector field method is to construct a vector field around the desired path, direct the UAV to move into the vector field, and ensure that the tracking error uniformly converges to 0. The L1 guidance rate is an algorithm that chooses a reference point on the current path and calculates the horizontal expected acceleration according to the horizontal velocity of the UAV and its distance from the reference point. The geometric method has the advantages of a simple structure and easy implementation. It is mostly used for UAVs following a two-dimensional curved path. It is mainly applicable for UAV cruise flights at a certain altitude, as it is hard to apply for the UAV’s whole flight phase. The control method is mainly applied in the controller design based on the error or error change rate of the two-dimensional or three-dimensional path following. It has a wide range of applications and can be used to design controllers according to specific requirements. Compared to PID control and model predictive control, sliding mode control can overcome the uncertainty of the system and has the advantages of a fast response, insensitivity to parameter changes and disturbances, and a simple structure, and thus has great advantages in path-tracking control algorithms for unmanned aerial vehicles. The main disadvantage of sliding mode control is that chattering can occur when state trajectory reaches the sliding mode surface, which is also the main factor limiting the practical application of sliding mode control. Integral sliding mode control (ISMC) can make its initial state on the sliding mode surface by design so as to eliminate the static error, effectively suppress the vibration, and improve the robustness of the system.
Liu et al. designed the S-Plane controller [18]. It features a simple structure, easy parameter adjustment, and a strong disturbance rejection capability. Zhao et al. introduced S-Plane control into the longitudinal motion control of a UAV [19]. It showed a good ability to resist external disturbance and modeling uncertainty. When Dong et al. introduced S-Plane control into a micro-USV control system, the simulation results showed that the S-Plane control had a good motion control performance [20]. Li et al. applied S-Plane to the control of an ocean observation platform to aid it in resisting the influence of ocean currents. Simulink and tests obtained good results [21].
In the movement of small fixed-wing UAVs, external wind disturbance is a critical and non-negligible disturbance factor. Zhao et al. took 20% to 50% of the UAV flight speed as the wind velocity [22]. Brezoescu et al. took the gust and crosswind disturbance as a fixed value [23]. Zhang et al. regarded 5 m/s fixed wind from west to east in the UAV path trajectory as an external wind disturbance [24]. In the above literature, the external wind disturbance was set to a certain value. This method is easy to operate in simulation, but it does not accurately reflect the impact of real external wind disturbance on UAVs. However, Wei et al. expressed the gust, random wind, and gradual wind of nature with corresponding equations, which truly reflects the composition of wind disturbance [25]. This paper will expand on this research on the composition of wind disturbance so as to more truly and accurately reflect the impact of external wind disturbance on UAV path following.
In this paper, based on the global stability theory, we propose a Globally Stable Integral Sliding Mode Radial Basis Function S-Plane (GSISM+RBF S-Plane) control method by combining the integral sliding mode control and the S-Plane control. It has the characteristics of strong disturbance rejection, resistance to modeling uncertainty, and robustness. In addition, the integral sliding mode control can improve the control accuracy. Compared with the Globally Stable Sliding Mode (GSSM) control and GSISM S-Plane control, GSISM+RBF S-Plane control has the characteristics of high tracking accuracy, a strong anti-disturbance ability, and smooth cutting into the ideal path. The specific contributions of this paper are as follows.
Aiming at the problems of low following accuracy and weak anti-disturbance ability of fixed-wing UAVs in three-dimensional path following, an S-Plane control algorithm with a strong anti-disturbance performance is introduced, and a radial basis function neural network is used to estimate the unknown disturbance in the model, and a GSISM+RBF S-Plane control algorithm is proposed. This allows the fixed-wing UAV to cut into the ideal path smoothly with a certain radian, and it has high following accuracy and anti-disturbance performance under external wind.
Considering the influence of wind on small fixed-wing UAVs, the composition of natural wind is modeled and studied, and these wind models are used to simulate the actual flight environment of fixed-wing UAVs in a simulation so as to improve the control performance of fixed-wing UAVs.
In order to verify the performance of the designed GSISM S-Plane control algorithm, the simulation is carried out considering a spatial straight line, spiral line, and special space curve, in which the special space curve includes taxiing, takeoff, cruise, and landing phases of the fixed-wing UAVs.
The remainder of this paper is organized as follows: UAV modeling and S-Plane control are described in Section 2. The design of the fixed-wing UAV controller is presented in Section 3. A comparative simulation and simulation verification are described in Section 4. Concluding remarks are summarized in Section 5.

2. UAV Modeling and S-Plane Control

2.1. UAV Modeling

The positional relationship between a UAV and an ideal path in a three-dimensional space is shown in Figure 1.
In Figure 1, ( x , y , z ) represents the position coordinates of the fixed-wing UAV at a certain point in space, and ( x d , y d , z d ) represents the position coordinates of the fixed-wing UAV on the ideal path in three-dimensional space. ( x e , y e , z e ) represents the error between the current position ( x , y , z ) of the fixed-wing UAV and the corresponding position ( x d , y d , z d ) on the ideal path. The motion model of fixed-wing UAV in three-dimensional space is shown in Equation (1) [26]:
x ˙ = v cos ψ cos γ + d x y ˙ = v sin ψ cos γ + d y z ˙ = v sin γ + d z
where ψ is the yaw angle, γ is the flight path angle, and v is the flight speed of the fixed-wing UAV. d x , d y , and d z are the unknown disturbances of each coordinate axis, respectively.

2.2. S-Plane Control

The S-Plane control model is shown in Equation (2):
u = 2 1 + exp [ ( k 11 e + k 12 e ˙ ) ] 1.0
where e is the error of the input signal, e ˙ is the error change rate of the input signal, and k 11 and k 12 represent the adjustment coefficients of the input signal error and error change rate, respectively. A diagram of its three-dimensional surface is shown in Figure 2.
S-Plane control is similar in structure to proportional differential control. Therefore, in the adjustment of k 11 and k 12 ,we can refer to the adjustment mode of its parameters in proportional differential control. PD control is linear control, while S-Plane control is a kind of nonlinear control based on the fuzzy control table. Thus, it also has the advantages of strong robustness and anti-disturbance of fuzzy control.
In the path following a fixed-wing UAV, the mathematical model of the ideal path is ζ = f ( t ) . It can be seen from Section 2 that q   ( x d ,   y d ,   z d ) is a point on the ideal path. The current position of the fixed-wing UAV is p   ( x ,   y ,   z ) . The purpose of path following is to make lim t p q = 0 by designing a control law so as to achieve x following x d , y following y d , z following z d . That is, the error of each coordinate axis in the three-dimensional space tends to be 0. The error equation is shown in Equation (3).
x e = x x d y e = y y d z e = z z d
The derivative of error can be obtained through the derivation of Equation (3) according to Equation (1). The derivative of error is shown in Equation (4).
x ˙ e = v cos ψ cos γ x ˙ d y ˙ e = v sin ψ cos γ y ˙ d z ˙ e = v sin γ z ˙ d
In order to facilitate the design of late control law, we make a simple replacement, such as Equation (5).
v cos ψ cos γ = U 1 v sin ψ cos γ = U 2 v sin γ = U 3

3. Controller Design

3.1. GSISM S-Plane Controller and Stability Analysis

In order to design the three-dimensional path following the controller of the fixed-wing UAV, a theorem from the literature [27] is adopted.
Lemma 1: The trivial solution of the system
η ˙ = α tanh ( k η )
is globally asymptotically stable for any selected constants α ,   k > 0 , and for each η ( 0 ) , the state η ( t ) converges exponentially to zero.
Proof: It is sufficient to prove that when t , η 0 . That is, when Equation (6) is asymptotically stable, Lyapunov is defined as V = η 2 / 2 .
Then,
V ˙ = η η ˙ = α η tanh ( k η )
since x tanh ( x ) = x e x e x e x + e x 0 , k η tanh ( k η ) 0 ; thus, V ˙ 0
if, and only if, η = 0 , V ˙ = 0
which completes the proof [28].
The above theorem is also used in reference [29] for the controller design of a wheeled mobile robot. The wheeled mobile robot is eventually able to track the ideal path with a certain accuracy. However, compared with fixed-wing UAVs, wheeled mobile robots have the problem of having fewer degrees of freedom. In this paper, based on the literature [27], a GSISM+RBF S-Plane control method is proposed to improve the path-following accuracy of a fixed-wing UAV. The specific control law is designed as follows:
Firstly, the x -axis direction is taken to design the control law, and the GSISM plane is taken as
s 1 = x e + c 1 0 t x e d t + a 1 tanh ( p 1 x e ) , c 1 , a 1 , p 1 0
The time derivative of s ˙ 1 is given by
s ˙ 1 = x ˙ e + c 1 x e + a 1 p 1 ( 1 / cosh ( p 1 x e ) ) 2 = x ˙ x ˙ d + c 1 x e + a 1 p 1 ( 1 / cosh ( p 1 x e ) ) 2 = u 1 + d x x ˙ d + c 1 x e + a 1 p 1 ( 1 / cosh ( p 1 x e ) ) 2
The control law in the x -axis direction can be designed as
u 1 = x ˙ d c 1 x e a 1 p 1 ( 1 / cosh ( p 1 x e ) ) 2 d x k 1 s 1
Because d x is unknown, d x is estimated using an RBF neural network. The formula of the RBF neural network is
h ( x ) = exp ( x c j 2 2 b j 2 )
Design d x = d 1 * T h ( x ) + ε x
where x is the input of the neural network, j is the number of nodes in the hidden layer of the network, h ( x ) is the output part of the Gaussian function, d 1 * is the ideal weight, and ε x is the approximation error.
In the x direction, the actual output of the RBF neural network is
d ^ x = d ^ 1 T h ( x )
Then, Equation (10) can be rewritten as
u 1 = x ˙ d c 1 x e a 1 p 1 ( 1 / cosh ( p 1 x e ) ) 2 d ^ x k 1 s 1
Substitute Equation (13) into Equation (9) to obtain
s ˙ 1 = u 1 + d x x ˙ d + c 1 x e + a 1 p 1 ( 1 / cosh ( p 1 x e ) ) 2 = ( x ˙ d c 1 x e a 1 p 1 ( 1 / cosh ( p 1 x e ) ) 2 d ^ x k 1 s 1 )   + d x * x ˙ d + c 1 x e + a 1 p 1 ( 1 / cosh ( p 1 x e ) ) 2 = d ^ x + d x k 1 s 1 = d ˜ x k 1 s 1
d ˜ x = d x d ^ x = d 1 * T h ( x ) + ε x d ^ 1 T h ( x ) = d ˜ x T h ( x ) + ε x
Define the Lyapunov function as
L 1 = 1 2 s 1 2 + 1 2 γ 1 d ˜ x T d ˜ x
The time derivative of L 1 is given by
L 1 = s 1 s ˙ 1 + γ 1 d ˜ x T d ˜ ˙ x = s 1 ( d ˜ x k 1 s 1 ) γ 1 d ˜ x T d ^ ˙ x = s 1 ( d ˜ x T h ( x ) ε x k 1 s 1 ) γ 1 d ˜ x T d ^ ˙ x =   d   ˜ x T ( s 1 h ( x ) + γ 1 d ^ ˙ x ) k 1 s 1 2 s 1 ε x
The adaptive rate of the designed x axis is
d ^ ˙ x = 1 γ 1 s 1 h ( x )
Then, L ˙ 1 = s 1 ( k 1 s 1 ε x ) = k 1 s 1 2 ε x s 1 ; when ε x is small enough, that is ε x 0 , L ˙ 1 0 .
Similarly, the GSISM function in the y direction can be designed as
s 2 = y e + c 2 0 t y e d t + a 2 tanh ( p 2 y e ) , c 2 , a 2 , p 2 0
The control law in the y -axis direction can be designed as
u 2 = y ˙ d c 2 y e a 2 p 2 ( 1 / cosh ( p 2 y e ) ) 2 d y k 2 s 2
Finally, the GSISM function in the z direction can be designed as
s 3 = z e + c 3 0 t z e d t + a 3 tanh ( p 3 z e ) , c 3 , a 3 , p 3 0
The control law in the z -axis direction can be designed as
u 3 = z ˙ d c 3 z e a 3 p 3 ( 1 / cosh ( p 3 z e ) ) 2 d z k 3 s 3
When L ˙ i 0 ( i = 1 , 2 , 3 ) , s i 0 ( i = 1 , 2 , 3 ) , according to the invariance principle of L a S a l l e , the closed-loop system is asymptotically stable.

3.2. Controller Structure

The structure of the fixed-wing UAV path-following controller proposed in this paper is shown in Figure 3.
GSISM+RBF S-Plane controller adopts the double-loop structure of inner and outer loop control. The running process of the GSISM+RBF S-Plane controller is as follows. Firstly, an ideal flight path is obtained through the track generator. Then, the ideal path is expressed in the form of a parametric equation and transmitted to the GSISM controller. Flight speed v , flight path angle γ , and yaw angle ψ are obtained under the action of the GSISM controller, where the flight speed v is transmitted to the motion model through the outer loop control. The flight path angle γ and yaw angle ψ are transmitted to the motion model after passing through the S-Plane control in the inner loop. The actual UAV flight path coordinate is obtained after the action of the motion model. A part of the signal is output, and another part is fed back to the input signal to participate in the calculation of the next step. By repeating this process, the fixed-wing UAV can follow the ideal path. When the flight path angle γ and yaw angle ψ are transmitted to the inner loop control, the error of the flight path angle γ and yaw angle ψ is obtained under the action of the feedback signal. Then, the error change rate is obtained under the action of the differentiator, and the obtained error and error change rate are transmitted to the S-Plane controller. Part of the signal is output to the kinematic model after the action of the S-Plane controller; another part of the signal is fed back to the input signal of the inner loop control to participate in the calculation of the next step in the inner loop control. In addition, the RBF neural network takes the error of the input expected position and output actual position as the input signal and its output signal is input to the GSISM controller as the estimated value of unknown disturbance. In this way, the following of the intermediate command signal can be realized. The intermediate command signal can be obtained by Equation (5). The yaw angle ψ of path following is
ψ = arctan ( U 2 U 1 )
The flight path angle γ of the path to be followed can be obtained according to Equations (5) and (23).
γ = arctan ( U 3 U 2 sin ψ )
The velocity v of the path to be followed is obtained according to Equations (5), (23), and (24).
v = U 3 sin γ
The input signal in S-Plane control is the error and error change rate. Integral explosion in the solution of the error change rate occurs easily. Therefore, the derivation process of the intermediate instruction signal is realized using the second-order differential method [28]. The second-order differentiator is shown in Equation (26):
x ˙ 1 = x 2 x ˙ 2 = 2 R 2 ( x 1 n ( t ) ) R x 2 y = x 2
where n ( t ) is the input signal, x 1 is the following of the signal, and x 2 is the estimation for the first-order derivative of the signal. The differentiator is an integral chain operation, which can suppress the noise in signal derivation.

3.3. GSISM+RBF S-Plane Control Algorithm

The path-following control algorithm based on GSISM+RBF S-Plane model for the fixed-wing UAV is shown as Algorithm 1.
Algorithm 1: GSISM+RBF S-Plane Controller
Outer loop input: expected path point q   ( x d ,   y d ,   z d ) . Current location point p   ( x ,   y ,   z ) . GSISM+RBF S-Plane control gain. d x ,   d y ,   d z .
Inner loop input: intermediate command signal ψ , γ , S-Plane gain k 11 , k 12 , k 21 , k 22 .
1: calculating the error and error change rate through Equations (3) and (4).
2: ( s 1 , s 2 , s 3 ) ( x e , y e , z e )
3: calculate u 1 , u 2 , u 3 through Equations (10), (20), and (22)
4: calculate ψ , γ , v through Equations (23)–(25)
5: calculating ψ ˙ and γ ˙ through Equation (26).
6: ( e ψ , e γ ) ( ψ ψ , γ γ )
7: ψ 2 / ( 1 + exp [ ( k 11 e ψ + k 12 e ˙ ψ ) ] ) 1.0
8: γ 2 / ( 1 + exp [ ( k 21 e γ + k 22 e ˙ γ ) ] ) 1.0
9: inner loop return: ψ , γ
10: ( x , y , z ) ( v cos ψ cos γ , v sin ψ cos γ , v sin γ )
11: outer loop return: x , y , z

4. Simulation and Results Analysis

4.1. Wind Disturbance Modeling

In the design of a small fixed-wing UAV controller, the factor of wind has a great influence on the performance of the designed controller. In the field of wind power generation, in-depth research has been carried out on the composition and performance of wind. Therefore, we can refer to these research results to study the impact of external wind disturbance on small fixed-wing UAVs. Wei et al. hypothesize that the wind occurring in nature is composed of basic wind, gust wind, gradual wind, and random wind and that the random combinations of these four types of wind represent most wind types occurring in nature [25]. However, the wind speed used in wind power generation is generally low, which is not fully applicable to wind disturbance modeling in the control process of a fixed-wing UAV. So, based on relevant research results, this paper further studies a wind disturbance model affecting the control of a fixed-wing UAV; adjustment factors, including gust, gradual wind, and random wind, are added to adjust the wind speed according to different needs.
(1)
Basic Wind
Basic wind describes an average change in the external wind. It is also the best wind speed for UAV flights. In most cases, the wind occurring in nature can be described as basic wind. We treat it as a constant when we simulate it.
V c = k ,   k   is   a   constant
(2)
Gust Wind
Gust is a wind whose wind direction remains unchanged and whose wind speed increases suddenly for a short time. It reflects a sudden change in the wind. The model of gust wind can be expressed as follows:
V g = 0 ,   t < t 1 k g V g max 2 1 cos ( 2 π ( t t 1 t g ) ) 0 ,   t > t 1 + t g ,   t 1 t t 1 + t g
where V g is the wind speed of the gust, V g max is the highest wind speed of the gust, t 1 is the start time of the gust, t g is the duration of the gust, and k g is the gust wind speed adjustment factor.
(3)
Gradual Wind
Gradual wind reflects a gradual nature in the wind, and its model can be expressed as follows:
V w = 0 ,   others k w V w max t t 2 t 3 t 2 ,   t 2 t t 3 k w V w max ,   t 3 < t t 3 + T
where V w is the wind speed of the gradual wind, V w max is the maximum wind speed of the gradual wind, t 2 is the start time of the gradual wind, t 3 is the end time of the gradual wind, T is the holding time of the gradual wind, and k w is the wind speed adjustment factor of the gradual wind.
(4)
Random Wind
Random wind reflects arbitrariness, uncertainty, and randomness of wind, and its model is shown as follows:
V n = k n V n max R a m ( 1 , 1 ) cos ( ω + φ )
where k n is the adjustment factor of random wind speed, V n is the random wind speed, V n max is the maximum random wind speed, R a m ( 1 , 1 ) is a random number between −1 and 1, and ω is the average distance of the wind speed fluctuation. The general value is 0.5 π 2 π ; φ is a random quantity uniformly distributed on 0 2 π . After modeling the above four kinds of wind, we can represent the influence of real external wind on a small fixed-wing UAV flight in any combination. This is of great significance for improving the control performance of small fixed-wing UAVs. The following are the simulation results obtained by randomly selecting a combination of several winds.
Bring V c = 2   m / s , k g = 0.5   m / s , V g max = 6   m / s , t 1 = 1 , t g = 7 , k w = 0.5 , V w max = 6 into the gust and gradual wind formula above. The results for a combination of basic wind, gradual wind, and gust can be obtained via the MATLAB simulation, as shown in Figure 4.
The simulation of the random wind function is complex, so we approximately estimate the arbitrary wind model of nature through the random number module, transfer Fcn module, and gain module in Simulink. Bring k g = 0.125   m / s , V g max = 8   m / s , t 1 = 15   m / s , t g = 3   m / s into the random wind and gust formula above. The simulation of a specific random wind and gust combination model is shown in Figure 5.

4.2. Semi-Physical Simulation System

The hardware in the loop simulation system used in this paper is shown in Figure 6. It consists of four parts, namely, a ground control station, a switch, a UAV, and a ROS simulation platform. During operation, communication between the ROS simulation platform, UAV, and ground control station is achieved through the switch; the ROS simulation platform sends the simulation results to the ground control station through the switch-and-subscribe method. Then, the ground control station sends the control command to the UAV’s airborne computer through the switch to control the flight of the UAV. The simulation principle of the ROS simulation platform is as follows. In Figure 6, the inner part of the red dotted box is the publisher creation process, and the rest is the subscriber creation process. During the creation of the publisher, when the blank message module outputs a ROS message (bus signal) to the bus assignment module, the GSISM+RBF S-Plane controller module designed in this paper generates a UAV path following specific points, which are relayed via a ROS message under the action of the bus assignment module. Then, the bus assignment module passes the ROS message containing the UAV path points onto the publish module. After receiving the bus signal, the publish module sends it to the topic of “/location.” During the creation process of the subscriber, the subscribe module subscribes to the messages sent to the “/my_topic” by the publisher and extracts the coordinates of the UAV path points from the message, and then transfers these to the ground control station. The UAV model includes the shape and size design of the UAV, mainly to maintain its similarity to the actual UAV. It also includes a dynamic model based on Newtonian mechanics, a sensor model considering sensors such as cameras and radars, and a controller model based on control theory and UAV control characteristics.

4.3. Simulation Test and Results Analysis

In order to verify the following accuracy and anti-disturbance performance of the GSISM+RBF S-Plane controller designed in this paper, the GSSM controller in reference [29] and the GSISM S-Plane controller are used for a comparative experiment. In addition, the control effect under different flight paths following a spatial straight line, a spiral line, and a special space curve is verified. The coefficient of two S-Plane controllers is shown in Table 1.
Where the unknown disturbances in Equation (1) are:
d x = 5 sin ( 0.5 t ) d y = 5 cos ( 0.5 t ) d z = 5 sin ( 0.5 t )
The coefficient of the RBF function is shown in Table 2.
The parameters of the GSISM function in adaptive estimation are a 1 i = c 1 i = p 1 i = 0.015   ( i = 1 , 2 , 3 ) . The GSISM+RBF S-Plane control coefficients for the UAV when tracking the spatial straight line, spiral line, and special spatial curve are shown in Table 3.

4.3.1. Spatial Straight Line Path-Following Simulation

The parameter expression when the ideal path followed by the fixed-wing UAV is a spatial straight line is shown by Equation (32):
ζ : x d = 50 + 4 t y d = 50 + 7 t z d = 50 + 5 t
The initial position of a fixed-wing UAV during takeoff when following a spatial straight line is calculated in [ 20 , 20 , 35 ] . Under the action of control laws, the fixed-wing UAV completes the flight process along the spatial straight line. A diagram of the spatial straight line followed by the fixed-wing UAV in the simulation is shown in Figure 7. The error diagram of each coordinate is shown in Figure 8. As can be seen from Figure 8, when the fixed-wing UAV flies in a spatial straight line, the GSISM+RBF S-Plane controller can cut into the ideal path smoothly and quickly compared with the GSSM controller and GSISM S-Plane control. Figure 9 shows the error diagram of each coordinate axis of the spatial straight line. It can be seen from Figure 9 that in the whole flight process, especially when a disturbance is added in at around 10–20 s, the error between GSISM+RBF S-Plane control and the desired path is basically 0, which is smaller than that achieved with GSSM control and GSISM S-Plane control. Therefore, the GSISM+RBF S-Plane controller has the advantages of high following accuracy and strong anti-disturbance ability.

4.3.2. Spiral Line Path following Simulation

When the ideal flight path of the fixed-wing UAV is a circle, in order to better reflect the following performance of the designed controller in practice, the simulation is carried out using a spiral instead of a circle. The parametric equation of the spiral line is shown in Equation (33).
ζ : x d = 100 cos ( t / 2 ) y d = 100 sin ( t / 2 ) z d = 20 + 5 t
When flying along the spiral line, the initial position of UAV takeoff is [ 0 , 0 , 0 ] in the GSISM+RBF S-Plane controller. The simulation results for the fixed-wing UAV following the spiral line flight path are shown in Figure 10.
It can be seen from Figure 10 that GSISM+RBF S-Plane control, GSSM control, and GSISM S-Plane control had good control performances and could track the ideal path with high accuracy. However, compared with GSSM control and GSISM S-Plane control, GSISM+RBF S-Plane control could cut into the ideal path with a certain radian. GSSM control and GSISM S-Plane control cannot track with high accuracy in the initial stage and only cut into the ideal path after a few seconds. Figure 11 is the path-following diagram of each coordinate axis when the UAV follows the spiral line. Figure 12 shows the error value between the ideal path and the actual flight path of each coordinate axis when the UAV follows the spiral line. In the whole flight segment of Figure 12, the error of GSISM+RBF S-Plane control is around 0. It can be seen that the path-following accuracy is particularly high. In addition, we added the wind disturbance designed in Section 4 at around 15–25 s into the spiral line path-following simulation. It can be seen from Figure 11 and Figure 12 that when there is a gust disturbance at 17–23 s, the fluctuation of the UAV flying along the ideal path is greater. Compared with GSSM control and GSISM S-Plane control, GSISM+RBF S-Plane control leads to less fluctuation in gust disturbance. It can be seen that its performance against external wind disturbance is superior.
Figure 13 and Figure 14 present input and output diagrams of the flight path angle and yaw angle in S-Plane control. S-Plane control has good anti-disturbance performance. It can be seen from Figure 13 and Figure 14 that the output value of the command signal is significantly different from the input value after being controlled by the S-Plane, regardless of the flight path angle or yaw angle. The corresponding value is large before entering the S-Plane controller, but after the action of the S-Plane controller, it is basically normalized to [ 1 , 1 ] . Although the intermediate command signal changes to a certain extent when passing through the S-Plane control, the path-following accuracy of the fixed-wing UAV does not change; this also shows that the GSISM+RBF S-Plane controller designed in this paper has good control performance.

4.3.3. Following Simulation of Special Space Curve

When a UAV flies in real space, its path can usually be synthesized by Section 4.3.1 and Section 4.3.2 spatial curves. In order to more realistically simulate the accuracy of the actual path-following performance of the UAV, the special space curve in reference [30] is used to simulate the whole flight process of a UAV in this paper. Its parametric equation is shown in Table 4.
Figure 15 presents a diagram of the special space curve path. In Figure 15, ① is the taxiing phase on the runway, ② is the takeoff phase of the fixed-wing UAV, ③ and ⑤ are the cruise flight phase, ④ is the flight phase under special space conditions, ⑥ is the landing phase, and ⑦ is the taxiing phase after the aircraft lands. When following along a special space curve, the initial position of the UAV when taxiing is [ 0 , 0 , 1 ] . When following a special space curve path, it can be seen from Figure 15 that the GSISM+RBF S-Plane control, the GSSM control, and the GSISM S-Plane control can follow the ideal path with high accuracy for the whole flight phase, showing good control performance. We then add the wind disturbance designed in Section 4 (20–70 s) into the special space curve path simulation in Figure 15. It can be seen from Figure 15, Figure 16 and Figure 17 that when wind disturbance is added, the disturbance of GSSM control and GSISM S-Plane control is more severe than that of GSISM+RBF S-Plane control. When there is no disturbance in Figure 17, the error of GSSM control and GSISM S-Plane control is about 1, while the error of GSISM+RBF S-Plane control basically tends to 0. It can be seen from the diagram of the special space curve path of each coordinate axis in Figure 16 that GSISM+RBF S-Plane control can smoothly track the ideal path, and the external wind disturbance has little impact on it. GSSM control and GSISM S-Plane control are more sensitive to external wind disturbance, showing large fluctuations. Therefore, it can be seen that the GSISM+RBF S-Plane controller designed in this paper has high control performance and the ability to resist external wind disturbance.
Figure 18 and Figure 19 present diagrams of the input and output for the flight path angle and yaw angle of the S-Plane controller in the special spatial curve path-following phase. It can be seen from Figure 18 and Figure 19 that except for the period when wind disturbance is added, the output signals at other times are normalized to the range of [ 1 , 1 ] under the action of the S-Plane controller. Although the output value changes compared with the input value, it does not affect the following accuracy of the UAV. On the contrary, it improves its anti-disturbance ability. It also shows that the S-Plane controller has good anti-disturbance performance.

4.3.4. Simulation Results Analysis

In order to better prove the anti-interference performance of the designed controller, the maximum perturbation values of the x, y, and z directions in simulations are shown in Table 3, Table 4 and Table 5.
As Shown in Table 5, Table 6 and Table 7, when tracking spatial straight, helical, or spatially special curves, the GSISM+RBF S-Plane controller has the smallest perturbation values in the presence of external wind disturbance. Compared to the GSISM+RBF S-Plane controller, the perturbation values of the GSISM S-Plane controller and GSISM S-Plane controller are higher. In other words, the GSISM+RBF S-Plane controller has a superior anti-wind disturbance performance.

5. Conclusions

Aiming at the path-tracking control problem of small fixed-wing UAVs, a GSISM+ RBF S-plane controller is designed in this paper, and its control performance is verified by tracking spatial straight lines, spiral lines, and special spatial curves. The following conclusions are reached:
(1)
The proposed controller can track the ideal path with high accuracy and a smooth cut into curved paths; thus, it has good control accuracy and anti-disturbance performance under external wind disturbance;
(2)
The proposed controller is based on the inner and outer loop control idea, which has the characteristics of a simple structure, easy realization, and practical application for fixed-wing UAV control;
(3)
The S-plane control can adjust the input signal well, and the second-order differentiators have certain advantages in suppressing the integral explosion problem during signal derivation; thus, the proposed controller demonstrates excellent performance in anti-interference.

Author Contributions

Conceptualization, G.Z. and P.C.; methodology, P.C. and G.Z.; validation, Z.C. and J.L.; formal analysis, P.C. and J.L.; writing—original draft preparation, P.C.; writing—review and editing, Q.Y. and Z.C.; project administration, J.L.; funding acquisition, P.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the key research and development program of Shanxi Province (202202020101001); the National Natural Science Foundation of China under Grant (51909245, 62003314); the Fundamental Research Program of Shanxi Province (202103021224187), and the Postgraduate Science and Technology Project of NUC (20221876).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fari, S.; Wang, X.M.; Roy, S.; Baldi, S. Addressing unmodeled path-following dynamics via adaptive vector field. A UAV test case. IEEE Trant. Aero Elect. System. 2020, 56, 1613–1622. [Google Scholar] [CrossRef]
  2. Chen, P.Y.; Zhang, G.B.; Guan, T.; Yuan, M.N.; Shen, J. The motion controller based on neural network s-plane model for fixed-wing UAVs. IEEE Access 2021, 9, 93927–93936. [Google Scholar] [CrossRef]
  3. Sujit, P.B.; Saripalli, S.; Sousa, J.B. Unmanned aerial vehicle path following a survey and analysis of algorithms for fixed-wing unmanned aerial vehicles. IEEE Control Syst. Mag. 2014, 34, 42–59. [Google Scholar]
  4. Cui, Z.Y.; Wang, Y. Nonlinear adaptive line-of-sight path following control of unmanned aerial vehicles considering sideslip amendment and system constraints. Math. Probl. Eng. 2020, 2020, 4535698. [Google Scholar] [CrossRef] [Green Version]
  5. Miao, J.M.; Wang, S.P.; Tomovic, M.M.; Zhao, Z.P. Compound line-of-sight nonlinear path following control of underactuated marine vehicles exposed to wind, waves, and ocean currents. Nonlinear Dyn. 2017, 89, 2441–2459. [Google Scholar] [CrossRef]
  6. Jung, W.; Lim, S.; Lee, D.; Bang, H. Unmanned aircraft vector field path following with arrival angle control. J. Intell. Robot Syst. 2016, 84, 311–325. [Google Scholar] [CrossRef]
  7. Capello, E.; Guglieri, G.; Ristorto, G. Guidance and control algorithms for mini-UAV autopilots. Aircr. Eng. Aerosp. Technol. 2017, 89, 133–144. [Google Scholar] [CrossRef]
  8. Christoph, H.; Dirk, A. Model predictive trajectory following for a ground vehicle in a heterogeneous rendezvous with a fixed-wing aircraft. IFAC-PapersOnLine 2020, 53, 15693–15698. [Google Scholar]
  9. Yang, J.; Liu, C.J.; Coombes, M.; Yan, Y.D.; Chen, W.H. Optimal Path Following for Small Fixed-Wing UAVs Under Wind Disturbances. IEEE Trans. Control. Syst. Technol. 2021, 29, 996–1008. [Google Scholar] [CrossRef]
  10. Shah, M.Z.; Samar, R.; Bhatti, A.I. Lateral track control of UAVs using the sliding mode approach: From design to flight testing. Trans. Inst. Meas. Control. 2015, 37, 457–474. [Google Scholar] [CrossRef]
  11. Castaneda, H.; Salas-Pena, O.S.; Leon-Morales, J. Extended observer based on adaptive second order sliding mode control for a fixed wing UAV. ISA Trans. 2017, 66, 226–232. [Google Scholar] [CrossRef]
  12. Dehghani, M.A.; Menhaj, M.B. Integral sliding mode formation control of fixed-wing unmanned aircraft using seeker as a relative measurement system. Aerosp. Sci. Technol. 2016, 58, 318–327. [Google Scholar] [CrossRef]
  13. Benkhoud, K.; Bouallègue, S. Dynamics modeling and advanced metaheuristics based LQG controller design for a Quad Tilt Wing UAV. Int. J. Dyn. Control. 2018, 6, 630–651. [Google Scholar] [CrossRef]
  14. Julian, C.; Bernardo, M.; Fatiha, N. Modeling and adaptive backstepping control for TX-1570 UAV path following. Aerosp. Sci. Technol. 2014, 39, 342–351. [Google Scholar]
  15. Muslimov, T.Z.; Munasypov, R.A. Consensus-based cooperative control of parallel fixed-wing UAV formations via adaptive backstepping. Aerosp. Sci. Technol. 2021, 109, 106416. [Google Scholar] [CrossRef]
  16. Whang, I.H.; Cho, S. LQR gain-schedule controller for vertical line following. Electron. Lett. 2010, 46, 991–992. [Google Scholar] [CrossRef]
  17. Cho, N.; Kim, Y.; Park, S. Three-Dimensional Nonlinear Differential Geometric Path-Following Guidance Law. J. Guid. Control Dyn. 2015, 38, 2366–2385. [Google Scholar] [CrossRef]
  18. Liu, X.M.; Xu, Y.R. S control of automatic underwater vehicles. Ocean Eng. 2001, 3, 81–84. [Google Scholar]
  19. Zhao, X.C.; Yuan, M.N.; Cheng, P.Y.; Xin, L.; Yao, L.B. Robust H-infinity/S-plane controller of longitudinal control for UAVs. IEEE Access 2019, 7, 91367–91374. [Google Scholar] [CrossRef]
  20. Dong, Z.P.; Wan, L.; Song, Y.F.; Li, Y.M. Design of control system for micro-usv based on adaptive expert s plane algorithm. Shipbuild. China 2017, 58, 178–188. [Google Scholar]
  21. Li, Y.; An, L.; Jiang, Y.Q.; He, J.Y.; Cao, J.; Guo, H.D. Dynamic positioning test for removable of ocean observation platform. Ocean Eng. 2018, 153, 112–121. [Google Scholar] [CrossRef]
  22. Zhao, S.L.; Wang, X.K.; Zhang, D.B.; Shen, L.C. Curved path following control for fixed-wing unmanned aerial vehicles with control constraint. J. Intell. Robot Syst. 2018, 89, 107–119. [Google Scholar] [CrossRef]
  23. Brezoescu, A.; Espinoza, T.; Castillo, P.; Lozano, R. Adaptive trajectory following for a fixed-wing UAV in presence of crosswind. J. Intell. Robot Syst. 2013, 69, 257–271. [Google Scholar] [CrossRef] [Green Version]
  24. Zhang, J.M.; Li, Q.; Cheng, N.; Liang, B. Nonlinear path-following method of fixed-wing unmanned aerial vehicles. J. Zhejiang Univ.-Sci. C Comput. Electron. 2013, 14, 125–132. [Google Scholar] [CrossRef]
  25. Wei, Y.L.; Han, S.X.; Shi, S.Z. The modelling and simulation of the combined wind speed in the wind power system. Renew. Energy Resour. 2010, 28, 18–20. [Google Scholar]
  26. Lugo-Cárdenas, I.; Salazar, S.; Lozano, R. Lyapunov based 3D path following kinematic controller for a fixed wing UAV. IFAC-PapersOnLine 2017, 50, 15946–15951. [Google Scholar] [CrossRef]
  27. Ailon, A.; Zohar, I. Controllers for trajectory tracking and string-like formation in Wheeled Mobile Robots with bounded inputs. In Proceedings of the Melecon 2010—2010 15th IEEE Mediterranean Electrotechnical Conference, Wuhan, China, 26–28 April 2010; pp. 1563–1568. [Google Scholar]
  28. Wang, X.H.; Liu, J.K. Differentiator Design and Application—Signal Filtering and Differentiation; Publishing House of Electronics Industry: Beijing, China, 2010. [Google Scholar]
  29. Liu, J.K. Sliding Mode Control Design and Matlab Simulation: The Design Method of Advanced Control System; Tsinghua University Press: Beijing, China, 2015. [Google Scholar]
  30. Cardoso, D.N.; Esteban, S.; Raffo, G.V. A new robust adaptive mixing control for trajectory tracking with improved forward flight of a tilt-rotor UAV. ISA Trans. 2021, 110, 86–104. [Google Scholar] [CrossRef]
Figure 1. Position relationship between fixed-wing UAV and ideal path.
Figure 1. Position relationship between fixed-wing UAV and ideal path.
Drones 07 00253 g001
Figure 2. Spatial structure of S-Plane control.
Figure 2. Spatial structure of S-Plane control.
Drones 07 00253 g002
Figure 3. GSISM S-Plane controller structure.
Figure 3. GSISM S-Plane controller structure.
Drones 07 00253 g003
Figure 4. Simulation results for a basic wind, gradual wind, and gust combination.
Figure 4. Simulation results for a basic wind, gradual wind, and gust combination.
Drones 07 00253 g004
Figure 5. Simulation results of a random wind and gust combination.
Figure 5. Simulation results of a random wind and gust combination.
Drones 07 00253 g005
Figure 6. Schematic diagram of hardware in the loop simulation system.
Figure 6. Schematic diagram of hardware in the loop simulation system.
Drones 07 00253 g006
Figure 7. Spatial straight line path-following diagram.
Figure 7. Spatial straight line path-following diagram.
Drones 07 00253 g007
Figure 8. Spatial straight line path-following diagram for each coordinate axis.
Figure 8. Spatial straight line path-following diagram for each coordinate axis.
Drones 07 00253 g008
Figure 9. Error diagram of each coordinate axis of the spatial straight line.
Figure 9. Error diagram of each coordinate axis of the spatial straight line.
Drones 07 00253 g009
Figure 10. Spiral line path-following diagram.
Figure 10. Spiral line path-following diagram.
Drones 07 00253 g010
Figure 11. Spiral line path-following diagram for each coordinate axis.
Figure 11. Spiral line path-following diagram for each coordinate axis.
Drones 07 00253 g011
Figure 12. Error diagram of each coordinate axis of the spiral line.
Figure 12. Error diagram of each coordinate axis of the spiral line.
Drones 07 00253 g012
Figure 13. Input and output diagram of pitch angle in S-Plane control.
Figure 13. Input and output diagram of pitch angle in S-Plane control.
Drones 07 00253 g013
Figure 14. Input and output diagram of yaw angle in S-Plane control.
Figure 14. Input and output diagram of yaw angle in S-Plane control.
Drones 07 00253 g014
Figure 15. Diagram of special space curve path-following simulation.
Figure 15. Diagram of special space curve path-following simulation.
Drones 07 00253 g015
Figure 16. Diagram of special space curve path-following simulation for each coordinate axis.
Figure 16. Diagram of special space curve path-following simulation for each coordinate axis.
Drones 07 00253 g016
Figure 17. Error diagram of each coordinate axis of special space curve.
Figure 17. Error diagram of each coordinate axis of special space curve.
Drones 07 00253 g017
Figure 18. Diagram of input and output for flight path angle of S-Plane control in special space curve simulation.
Figure 18. Diagram of input and output for flight path angle of S-Plane control in special space curve simulation.
Drones 07 00253 g018
Figure 19. Diagram of input and output of yaw angle of S-Plane control in special space curve simulation.
Figure 19. Diagram of input and output of yaw angle of S-Plane control in special space curve simulation.
Drones 07 00253 g019
Table 1. The coefficient of two S-Plane controllers.
Table 1. The coefficient of two S-Plane controllers.
k 11 k 12 k 21 k 22
two S-Plane controllers40.01100.1
Table 2. The coefficient of the RBF function.
Table 2. The coefficient of the RBF function.
c i j b j γ 1 γ 2 γ 3
RBF function 0.1 0.5 0.5 0.25 0.25 0 0 0.25 0.25 0.5 0.5 15505050
Table 3. The GSISM+RBF S-Plane control coefficients for the UAV when tracking the spatial straight line, spiral line, and special spatial curve.
Table 3. The GSISM+RBF S-Plane control coefficients for the UAV when tracking the spatial straight line, spiral line, and special spatial curve.
k1k2k3c1c2c3a1p1a2p2a3p3
spatial straight line0.50.01650.0932020400.251.120.2
spiral line730300.200.30.320.213.41020.2
special spatial curve4598101520.20.51.10.20.2
Table 4. Parametric equation of spatial space curve.
Table 4. Parametric equation of spatial space curve.
t x d y d z d
t < 10 tt1
10 t < 30 2 t 10 + 10 1.5 t 10 + 10 0.5 t 10 + 1
30 t < 39.2 235 23 t 30 + 50 50 23 t 30 20 11
39.2 t < 84.44 144 cos 15 144 t 39.2 π 2 + 144 144 sin 15 144 t 39.2 π 2 + 144 11 + 6 sin 2 π 22.62 t 39.2
84.44 t < 93.64 50 23 t 84.44 144 235 23 ( t 84.44 ) 11
93.64 t < 113.64 1.5 t 93.64 20 50 2 t 93.64 11 0.5 t 93.64
113.64 t < 123.64 123.64 t 123.64 t 1
Table 5. The maximum perturbation values of the x direction.
Table 5. The maximum perturbation values of the x direction.
GSISM+RBF S-PlaneGSSMGSISM S-Plane
Spatial straight line0.5−1−3
Spiral line−0.8−15
Special space curve−0.4−1.2−2
Table 6. The maximum perturbation values of the y direction.
Table 6. The maximum perturbation values of the y direction.
GSISM+RBF S-PlaneGSSMGSISM S-Plane
Spatial straight line−0.4−4−2.2
Spiral line0.5−4.8−3
Special space curve35−4
Table 7. The maximum perturbation values of the z direction.
Table 7. The maximum perturbation values of the z direction.
GSISM+RBF S-PlaneGSSMGSISM S-Plane
Spatial straight line−1.52−5
Spiral line510−10
Special space curve−1.73−9
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

Chen, P.; Zhang, G.; Li, J.; Chang, Z.; Yan, Q. Path-Following Control of Small Fixed-Wing UAVs under Wind Disturbance. Drones 2023, 7, 253. https://doi.org/10.3390/drones7040253

AMA Style

Chen P, Zhang G, Li J, Chang Z, Yan Q. Path-Following Control of Small Fixed-Wing UAVs under Wind Disturbance. Drones. 2023; 7(4):253. https://doi.org/10.3390/drones7040253

Chicago/Turabian Style

Chen, Pengyun, Guobing Zhang, Jiacheng Li, Ze Chang, and Qichen Yan. 2023. "Path-Following Control of Small Fixed-Wing UAVs under Wind Disturbance" Drones 7, no. 4: 253. https://doi.org/10.3390/drones7040253

APA Style

Chen, P., Zhang, G., Li, J., Chang, Z., & Yan, Q. (2023). Path-Following Control of Small Fixed-Wing UAVs under Wind Disturbance. Drones, 7(4), 253. https://doi.org/10.3390/drones7040253

Article Metrics

Back to TopTop