Next Article in Journal
Study on the Influence of Measuring AE Sensor Type on the Effectiveness of OLTC Defect Classification
Next Article in Special Issue
Optimization and Communication in UAV Networks
Previous Article in Journal
Shifting the Specificity of E. coli Biosensor from Inorganic Arsenic to Phenylarsine Oxide through Genetic Engineering
Previous Article in Special Issue
Priority-Based Data Collection for UAV-Aided Mobile Sensor Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Decentralized Low-Chattering Sliding Mode Formation Flight Controller for a Swarm of UAVs

by
Thiago F. K. Cordeiro
1,
João Y. Ishihara
2 and
Henrique C. Ferreira
2,*
1
Gama Campus, University of Brasília, Brasília 72444-240, Brazil
2
Department of Electrical Engineering, Darcy Ribeiro Campus, University of Brasília, Brasília 70910-900, Brazil
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(11), 3094; https://doi.org/10.3390/s20113094
Submission received: 30 March 2020 / Revised: 20 April 2020 / Accepted: 7 May 2020 / Published: 30 May 2020
(This article belongs to the Special Issue Optimization and Communication in UAV Networks)

Abstract

:
In this paper, a nonlinear robust formation flight controller for a swarm of unmanned aerial vehicles (UAVs) is presented. It is based on the virtual leader approach and is capable of achieving and maintaining a formation with time-varying shape. By using a decentralized architecture, the local controller in each UAV uses information only from the UAV itself, its neighbors, and from the virtual leader. Also, a synchronization control objective provides a mechanism to weight between the fleet achieving the desired formation shape, that is, achieving the desired relative position between the UAVs, and each UAV achieving its desired absolute position. The use of a combination of a sliding mode controller and a low pass filter reduces the usual chattering effect, providing a smooth control signal while maintaining robustness. Simulation results show the effectiveness of the proposed decentralized controller.

1. Introduction

The use of an unmanned aerial vehicle (UAV) swarm brings several advantages in search and rescue, disaster monitoring, aerial mapping, traffic monitoring, reconnaissance missions, and surveillance [1,2,3]. A swarm of UAVs provides system redundancy, reconfiguration ability, and structure flexibility, being more effective, flexible, robust, and reliable than single vehicles [4,5]. The formation control is a critical task of attempting cooperation among UAVs. In general, a formation control problem is to find a coordination scheme to enable UAVs to reach and maintain some desired, possibly time-varying formation or group configuration [6].
In the view of communication networks, the existing formation control approaches can be classified into the centralized method, where a single controller is used to control the whole team based on the information from the whole team [7] and the decentralized method, where each team member generates its own control based on local information from its neighbors [1,2,4,8,9,10,11]. Centralized formation control can be a good strategy for a small team of UAVs. When considering a team with a large number of UAVs, the need for greater computational capacity and a large communication bandwidth would mandate a decentralized formation control [4].
The main structures considered for formation control of UAVs swarm are leader-follower, behavioral, and virtual structure/virtual leader [12,13]. In the leader-follower approach [3,6,10,14], a common leader is chosen and the rest of the agents are assigned as followers. The group leader broadcasts its position information to the followers who then begin to follow the leader at an offset. In the behavioral approach [15,16], several desired behaviors are prescribed for agents in this approach. Such desired behaviors may include cohesion, collision avoidance, obstacle avoidance. In the virtual structure/virtual leader approach [4,9,11], the entire formation is treated as a single rigid body. The virtual structure can evolve as a whole in a given direction with some given orientation and can maintain a rigid geometric relationship among multiple vehicles. In the virtual leader approach, the leader is a known virtual entity and its information can be made available for each agent software.
There are several control techniques used in UAV formation control, based on distinct premises and aiming to achieve distinct objectives. A common approach is to use a nonlinear dynamic inversion (NLDI) which, via nonlinear functions, encapsulates the nonlinear system in a box with virtual inputs/outputs that behaves as a linear system. This linearized system act as a set of double integrators, that then is controlled by any linear or nonlinear technique, such as pole placement [14,17], H control [9], differential game approach [10] or sliding mode control (SMC) [8,18]. There are two main approaches when using the NLDI fixed-wing UAV formation flight control, related to in which frame the whole formation is described. One is to choose a global frame, such as north-east-down [9,10], and the virtual inputs are accelerations in north/south, east/west and up/down directions. Other is to use a leader related frame [8,14,17,18], and the virtual inputs accelerates toward forward/backward, left/right and above/below the leader.
In References [14,17], classic controllers were designed after applying of NLDI procedure in the nonlinear dynamics of UAVs formation flight. In Reference [10], a differential game approach is used to achieve an optimal controller that weights between minimizing the terminal position and velocity error of each UAV and minimizing the control effort. Another option, as in Reference [2], is the model predictive control (MPC), which can be used to compute an optimal control output to achieve formation control while avoiding obstacle and dealing with actuator saturation. It is however computationally expensive, since it reevaluates, at each time instant, the optimal control output over a finite time horizon. In Reference [2], the computational cost is partially reduced by maintaining the previously computed control output and reevaluating only when certain trigger events indicates that the control output must be changed, which works well in steady maneuvers, such as straight level flight or in constant-rate turns. In all of these approaches [2,10,14,17], the project does not account for the effect of disturbances or model uncertainty. Robust [8,9,11,18,19] and adaptive [1] approaches are appropriate for tackling this problem, where robust approaches usually has fast response but has a high control effort and/or chattering, whereas adaptive approach has a slower convergence, but uses a smoother control signal. In this way, the robust approach is recommended if precision of the formation is more important than control effort. In Reference [9], the proposed H linear controller is robust to noises, disturbances, and delays in communication between the UAVs. The SMC is an interesting technique since it, ideally, can completely compensates the effects of model uncertainties and bounded disturbances. As disadvantage, it provides a discontinuous, chattering control signal whose source is a signum (sign) function [11,19,20]. A possible solution is to change this signum function to a saturation (sat) function, as in our previous work for UAV formation flight [18], but this generates a trade-off between precision and chattering. Another solution is the use of a second-order SMC (SOSMC) [21,22], which uses the integral of the chattering signal as control input of a plant. A generalization of the second-order SMC is the low-pass filter (LPF) [23,24,25,26]. The integral, or the more general low-pass filter, is included as part of the plant model, which improves precision. For example, in Reference [23] an architecture with sliding mode control and low pass filter is proposed for synchronous position control for multiple robotic manipulator systems. However, its control law involves the computation of derivatives whose order is higher than the plant order, which can difficult implementation for example in an embedded system. In Reference [24] an attitude controller for the reentry of a space vehicle based on low pass filter SMC architecture is proposed. Differently from, for example, Reference [23], the LPF in Reference [24] is used to filter only the signal component that contains a signum function, while bypassing the smooth component of the control signal directly to the plant. This approach, compared to the approach in Reference [23], avoids the computation of higher-order derivatives.
In this paper, using an LPF-based SMC approach, a decentralized controller for a time-varying synchronized formation of multiple UAVs with a virtual leader is proposed. It is considered that each UAV is subject to unknown bounded disturbance. The computation of higher-order derivatives is not required. This is achieved by decomposing the control signal in smooth and in chattering components and filtering only the chattering component. Compared with Reference [24], which considers a single space vehicle, the proposed controller considers a synchronized and decentralized formation of multiple UAVs. In a synchronized formation, multiple UAVs simultaneously converge to desired positions. In comparison with Reference [23], which uses LPF SMC for synchronized position control for multiple robotic manipulator systems in a ring-link communication topology, the proposed controller not require computation of higher-order derivatives, a more general information exchange topology is adopted, and the problem of UAV formation flying is considered. Different from our previous work [18], the proposed controller uses an LPF for chattering attenuation. The finite-time convergence to a linear sliding surface is proven by introduction of an appropriated Lyapunov function candidate and simulation results show the effectiveness of the proposed control architecture.
To use the LPF-based SMC, the upper bound of the disturbance must be known. Most of this disturbance is better described in the wind frame of the aircraft. For example, a model uncertainty can affect the lift force computation. The difference between the true and computed lift forces is equivalent to a disturbance force applied in the lift direction. Similar discussion can be made of the thrust, drag and side forces. If the NLDI linearizes the system in the leader’s wind frame [8,14,17,18], this upper bound can be used directly, under the assumption that the leader’s wind frame is similar enough to the followers’ wind frames, as the fleet in formation flies approximated to the same direction. If, however, the NLDI linearizes the system in a global frame [9,10], the wind-frame-described upper bound must be translated to the global frame. The equations to translate the disturbance upper bound to the global frame are developed in this paper.
The main contribution is now summarized. A formation flight controller is developed that includes, in a single controller, the following characteristics:
  • uses the robust sliding mode control technique [8,11,18,19];
  • has low-chattering with low degradation in performance by the use of a low-pass filter modelled as a plant component [23,24,25,26];
  • uses a variant of the LPF SMC that is mathematically and computationally simpler than the usual approach, and removes computation of higher-order derivatives [24];
  • is a multi-agent decentralized/synchronous approach [18,23].
It is worth noting that each individual characteristic of the controller listed above has already been developed in other papers but, to the best knowledge of the authors, there is no controller that includes all characteristics in a single controller. It is also worth noting that, to include all characteristics in a single controller, an appropriate Lyapunov that unifies theses characteristics is developed.
As a second contribution, a set of equations that translate the disturbance upper bound and is derivative upper bound from the wind frame to the global frame is proposed. These equations are used in the proposed LPF-SMC, but can be used, with minor modifications, in most fixed-wing formation SMC or SOSMC that are described in a global frame.
The remainder of this paper is organized as follows. Section 2 defines the problem, presents the mathematical models for the UAVs, formation flight, and communication graph. Section 3 presents the proposed controller, equations to compute the disturbance’s upper bound, and proves the stability of the controller. Section 4 evaluates the proposed controller by simulation against an unfiltered SMC, where it is shown that the controller significantly reduces its chattering without significantly reducing its performance, and Section 5 concludes this paper.

2. Preliminaries

In this section, models for an individual UAV and for a fleet formation are presented.

2.1. UAV Model

The dynamics relating the input and output of the i-th vehicle in a fleet of n UAVs can be described by using the so-called point-mass aircraft model. It assumes a non-rotating flat Earth with a constant gravitational acceleration g. This model provides adequate precision to aircraft guidance and control problems and for short-range trajectory planning. It also assumes that the intensity of the wind is mild such that the airflow can be considered aligned with the vehicle fuselage, that is, that the angle-of-attack and sideslip angle are null, which are reasonable suppositions to cruise flight and coordinate maneuvers. Under these assumptions, as depicted in Figure 1, the drag force D i ( t ) , generated by the airflow, is aligned to the fuselage, pointing backward, whereas the lift force L i ( t ) is perpendicular to the fuselage. It is assumed that the propulsive system provides a thrust force vector T i ( t ) aligned with the fuselage/airflow and in the opposite direction to D i ( t ) and that the net angular momentum generated by the propulsive system is null. Finally, it is assumed that the vehicle mass m i is (approximately) constant, i.e., the propulsive system is electric or if fuel-based, that the consumption is small compared to the total vehicle mass. These simplified suppositions are commonly used in the literature [9,10,27,28]. To achieve extra precision in more aggressive maneuvers, the effect of the angle-of-attack and sideslip angle should be included in the model [14]. In general, these quantities require a dedicated sensory system—which is usually not present in small UAVs—to these angles be measured. In this work, the angle-of-attack and sideslip angle are allowed to be unmeasured but it is supposed that their effect can be incorporated in the model as a bounded disturbance b i ( t ) .
The state vector of the point-mass model of the i-th UAV is composed of its position vector p i ( t ) = [ p x i ( t ) p y i ( t ) p z i ( t ) ] T described in the inertial Cartesian reference frame NED (North-East-Down) and by its velocity, described in a spherical coordinate system composed by the ground speed V i ( t ) , flight path angle γ i ( t ) and course angle χ i ( t ) . By rotating the reference frame first by χ i ( t ) around the z axis and after by γ i ( t ) around the rotated y axis, the i-th aircraft wind frame is obtained. Since it is assumed that the aircraft velocity vector is aligned with its fuselage, χ i ( t ) and γ i ( t ) are equivalent respectively to the yaw ψ i ( t ) and pitch θ i ( t ) attitude angles. The point-mass model includes also the roll attitude angle ϕ i ( t ) , which is a rotation of the fuselage around the direction of the velocity vector. The definition of the roll, pitch and yaw angles can be seen in [29]. Figure 1 shows the i-th UAV and its vectors and attitude angles.
The state change given by the derivative of p i ( t ) is computed as
p ˙ i ( t ) = p ˙ x i ( t ) p ˙ y i ( t ) p ˙ z i ( t ) = R i ( t ) V i ( t ) 0 0 = V i ( t ) cos γ i ( t ) cos χ i ( t ) cos γ i ( t ) sin χ i ( t ) sin γ i ( t ) ,
where
R i ( t ) = cos γ i ( t ) cos χ i ( t ) sin χ i ( t ) sin γ i ( t ) cos χ i ( t ) cos γ i ( t ) sin χ i ( t ) cos χ i ( t ) sin γ i ( t ) sin χ i ( t ) sin γ i ( t ) 0 cos γ i ( t )
is a rotation matrix that rotates from the wind frame to the reference frame with angular velocity
ω i ( t ) = [ χ ˙ i ( t ) sin γ i ( t ) γ ˙ i ( t ) χ ˙ i ( t ) cos γ i ( t ) ] T ,
and the variables V i ( t ) , χ i ( t ) , and γ i ( t ) can be computed by
tan χ i ( t ) = p ˙ y i ( t ) p ˙ x i ( t ) , sin γ i ( t ) = p ˙ z i ( t ) V i ( t ) , V i 2 ( t ) = p ˙ x i 2 ( t ) + p ˙ y i 2 ( t ) + p ˙ z i 2 ( t ) .
The UAV dynamics is described by [9]
V ˙ i ( t ) = T i ( t ) D i ( t ) m i g sin γ i ( t ) + b t i ( t ) , χ ˙ i ( t ) = L i ( t ) sin ϕ i ( t ) m i V i ( t ) cos γ i ( t ) + b ψ i ( t ) V i ( t ) cos γ i ( t ) , γ ˙ i ( t ) = L i ( t ) cos ϕ i ( t ) m i V i ( t ) g cos γ i ( t ) V i ( t ) + b θ i ( t ) V i ( t ) ,
where the disturbance signal b i ( t ) = [ b t i ( t ) b θ i ( t ) b ψ i ( t ) ] encompasses model approximations, parameter uncertainty, and disturbances in acceleration, generated by several sources, such as wind. It is supposed that b i ( t ) and b ˙ i ( t ) are unknown but with known bounds. The subscripts t, θ , and ψ from the elements of b i ( t ) means respectively thrust, pitch, and yaw. The thrust force magnitude T i is a function of the engine throttle; D i ( t ) is the magnitude of the drag force; the lift force magnitude L i is a function of several parameters, such as air density and aircraft speed, and is adjusted mainly by changing the elevator position and ϕ i is adjusted by a combination of aileron and rudder positions. The variables T i ( t ) , L i ( t ) , and ϕ i ( t ) are the control inputs at the i-th UAV. It can be seen that Equation (5) presents a singularity when V i ( t ) = 0 . Fixed-wing vehicles must maintain non-null airspeed to maintain its lift. Assuming null or mild wind speed, the ground speed V i ( t ) is also non-null and this singularity does not occur. It can be seen that cos γ i ( t ) = 0 also presents a singularity in Equation (5). This occurs only when the UAV is flying exactly in an up or down direction. However, this is not an achieved state, except in highly acrobatic vehicles.
By defining the load factor n i ( t ) as [28]
n i ( t ) L i ( t ) m i g ,
and defining the following virtual control input
Γ i ( t ) = a t i ( t ) a ψ i ( t ) a θ i ( t ) = T i ( t ) D i ( t ) m i g n i ( t ) sin ϕ i ( t ) g n i ( t ) cos ϕ i ( t ) ,
the dynamics Equation (5) can be rewritten as
V ˙ i ( t ) = a t i ( t ) g sin γ i ( t ) b t i ( t ) , χ ˙ i ( t ) = a ψ i ( t ) b ψ i ( t ) V i ( t ) cos γ i ( t ) , γ ˙ i ( t ) = a θ i ( t ) g cos γ i ( t ) b θ i ( t ) V i ( t ) .
By deriving p ˙ i ( t ) from Equation (1), and applying some manipulations, is obtained that
p ¨ i ( t ) = R i ( t ) Γ i ( t ) + b i ( t ) + g ,
where g = [ 0 0 g ] T is the gravitational acceleration vector. By defining
τ i ( t ) = [ τ x i ( t ) τ y i ( t ) τ z i ( t ) ] T R i ( t ) Γ i ( t ) + g ,
d i ( t ) = [ d x i ( t ) d y i ( t ) d z i ( t ) ] T , R i ( t ) b i ( t ) ,
the dynamics are finally rewritten as
p ¨ i ( t ) = τ i ( t ) + d i ( t ) ,
where τ i ( t ) R 3 is a virtual controller input and d i ( t ) R 3 is the virtual disturbance described in the reference frame.
For the controller design, the model given by Equation (12) will be used. Once the virtual control signals are known, the original variables can be obtained. Since for any rotation matrix, R i 1 ( t ) = R i T ( t ) , the virtual input Γ i ( t ) can always be obtained from τ i ( t ) by
Γ i ( t ) = R i T ( t ) τ i ( t ) g ,
and then T i ( t ) , n i ( t ) , ϕ i ( t ) can be obtained from Equation (7), which can finally be used as the input of an inner loop controller that actuates over the engine and control surfaces [8].

2.2. UAV Formation

It is considered a formation of UAVs with a virtual leader scheme. The virtual leader is designated here as the 0-th UAV, and consists of a virtual point with a position p 0 ( t ) = [ p x 0 ( t ) p y 0 ( t ) p z 0 ( t ) ] T in space, known by all UAVs, which describes a smooth trajectory as a function of time. The results of this work can also be used for a non-virtual leader configuration by assuming that the leader UAV can broadcast its position to all followers UAVs.
The fleet formation is planned by the generation of the desired position p i d ( t ) = [ p x i d ( t ) p y i d ( t ) p z i d ( t ) ] T for the i-th UAV which is described as
p i d ( t ) = p 0 ( t ) + p ˜ i ( t ) ,
where p ˜ i ( t ) = [ p ˜ x i ( t ) p ˜ y i ( t ) p ˜ z i ( t ) ] T is the desired (time-varying) clearance, which is described in the reference frame.
To achieve a formation shape that rotates with the leader is interesting to describe the desired clearance p ˜ i ( t ) in the leader’s wind frame or any other frame related to the leader as
p ˜ i ( t ) = R r ( t ) p ˜ i r ( t ) ,
where p ˜ i r ( t ) is the clearance vector described in a leader’s frame, such as wind, and the formation rotation matrix R r ( t ) rotates from the leader’s frame to the reference frame. For example, by defining R r ( t ) = R 0 ( t ) (see Equation (2) with i = 0 ), it is achieved formation description aligned with the (virtual) leader’s trajectory as in, for example, Reference [8,9]. If, instead, R r ( t ) is defined as
R r ( t ) = R χ ( t ) cos χ 0 ( t ) sin χ 0 ( t ) 0 sin χ 0 ( t ) cos χ 0 ( t ) 0 0 0 1 ,
it is achieved a formation description aligned with the horizontal projection of the (virtual) leader’s trajectory, used in, for example, References [14,17].
Another option is to describe the formation using a leader’s frame defined by the attitude Euler angles yaw ψ 0 ( t ) , pitch θ 0 ( t ) , and roll ϕ 0 ( t ) . This can be useful, for example, for maneuvers involving close interaction between the leader and the followers, such as to a boom-receptacle automatic aerial refueling. In this case, R r ( t ) = R b ( t ) , where [29]
R b ( t ) cos ψ 0 cos θ 0 cos ψ 0 sin θ 0 sin ϕ 0 sin ψ 0 cos ϕ 0 cos ψ 0 sin θ 0 cos ϕ 0 + sin ψ 0 sin ϕ 0 sin ψ 0 cos θ 0 sin ψ 0 sin θ 0 sin ϕ 0 + cos ψ 0 cos ϕ 0 sin ψ 0 sin θ 0 cos ϕ 0 cos ψ 0 sin ϕ 0 sin θ 0 cos θ 0 sin ϕ 0 cos θ 0 cos ϕ 0 .
The derivatives of p i d ( t ) in Equation (14) can be computed as
p ˙ i d ( t ) = p ˙ 0 ( t ) + p ˜ ˙ i ( t ) ,
p ¨ i d ( t ) = p ¨ 0 ( t ) + p ˜ ¨ i ( t ) .
Using the Theorem of Coriolis [29], the derivatives of p ˜ i ( t ) in Equation (15) can be computed as
p ˜ ˙ i ( t ) = R r ( t ) p ˜ ˙ i r ( t ) + ω r ( t ) × p ˜ i r ( t ) ,
p ˜ ¨ i ( t ) = R r ( t ) p ˜ ¨ i r ( t ) + 2 ω r ( t ) × p ˜ ˙ i r ( t ) + ω ˙ r ( t ) × p ˜ i r ( t ) + ω r ( t ) × [ ω r ( t ) × p ˜ i r ( t ) ] ,
where ω r ( t ) is the angular velocity between the rotating leader’s frame and the reference frame and is given by
ω r ( t ) = leader s gyro measurements , if R r ( t ) = R b ( t ) , [ χ ˙ 0 sin γ 0 γ ˙ 0 χ ˙ 0 cos γ 0 ] T , if R r ( t ) = R 0 ( t ) , [ 0 0 χ ˙ 0 ] T , if R r ( t ) = R χ ( t ) .
It is worth noting that when using the non-virtual leader’s body frame, the angular velocity ω r ( t ) is the body angular velocity, which can be directly measured by a gyro sensor at the leader. By using a non-virtual leader and any of the wind frame variants, the ground velocity obtained from a GPS sensor or from a navigation algorithm must be used. When using a virtual leader approach, its trajectory is smooth, pre-known, and artificially generated, in a way that ω r ( t ) can be pre-computed analytically or numerically with arbitrary precision depending on how the trajectory is created.

2.3. Communication Graph

Each follower UAV can exchange data with their neighbors. The communication network is represented by an undirected graph, which means that, if an i-th UAV receives data from a j-th UAV, this means that the j-th UAV receives data from the i-th UAV. The set of the UAVs that are neighbors of the i-th UAV is defined as N i .
The Laplacian matrix L represents the connectivity between the UAVs
L i j = a i j , if j i and j N i , k N i a i k , if j = i , 0 , otherwise ,
where a i j = 0 means that there is no communication between the i-th and j-th UAVs and a i j > 0 means that there is a communication link between the i-th and j-th UAVs, and the value of a i j is used as a weight to the control algorithm that is developed in this paper. If all UAVs are reachable, that is, if someone starts from any UAV and can achieve any-other UAV via the communication links, L is semidefinite positive.
In decentralized controllers, the weight that is given to the information present in the own i-th UAV is also described, which is given by the diagonal matrix Λ . The matrix H includes both the own weight and the neighborhood weight. These matrices are given by
H = Λ + L ,
Λ = diag ( [ λ 1 λ n ] ) .
Note that since λ 1 , , λ n > 0 and L is semidefinite positive, the matrix H is invertible.

2.4. Formation Tracking and Synchronization Errors

The tracking error of each aircraft e i ( t ) = [ e x i ( t ) e y i ( t ) e z i ( t ) ] T R 3 , relative to a desired position in the reference frame, is defined as
e i ( t ) p i ( t ) p i d ( t ) .
The synchronization error Δ e i j ( t ) = [ Δ e x i j ( t ) Δ e y i j ( t ) Δ e z i j ( t ) ] T R 3 , which can be seen as a relative position error between the UAVs, is defined as
Δ e i j ( t ) e i ( t ) e j ( t ) = p i ( t ) p ˜ i ( t ) p j ( t ) p ˜ j ( t ) .
It can be seen that Δ e i j ( t ) can be computed without knowing the leader’s position. However, since the computation of p ˜ i ( t ) and p ˜ j ( t ) in Equation (15) can be chosen to be dependent on the leader’s flight direction or attitude angles, it is assumed here that the leader’s data is available to all UAVs.
It is assumed that each i-th UAV can communicate only with a correspondent set of neighbor UAVs, N i { 1 , 2 , , n } . The communication graph is assumed to be undirected, connected, not change with time, and previously known. Each UAV receives the tracking error information of other UAVs in the fleet only through its neighbors (as, for example, in the simulation in Section 4). The virtual leader can be seen as an extra node in the graph, that connects to every other UAV in a directed way, from leader to each follower.
The coupled error at i-th UAV is defined as the weighted sum of its tracking error and the synchronization error with respect to its neighbors, that is,
e i c ( t ) = [ e x i c ( t ) e y i c ( t ) e z i c ( t ) ] T λ i e i ( t ) + j N i a i j Δ e i j ( t ) = λ i e i ( t ) + j = 1 n a i j Δ e i j ( t ) ,
in which λ i > 0 weights its own tracking error and a i j > 0 weights the error difference between the neighbor UAV j of the UAV i. In the last equality in Equation (28), if j N i then a i j = 0 . The synchronization control objective is to make the coupled errors approach to zero.

2.5. A Componentwise Formation Description

It is supposed that each component of d i ( t ) is independent of each other which implies that each component of p ¨ i ( t ) is independent of each other. In this way, the controller design is simplified since the description of only one axis is sufficient. A controller policy can be developed to a single axis and then it can be directly applied to the other two.
The one-dimensional dynamics from axis l = x , y , z , of the reference frame, is obtained from Equation (12) as
p ¨ l i ( t ) = τ l i ( t ) + d l i ( t ) .
Accordingly, the coupled tracking-synchronization error is obtained from Equation (28) as
e l i c ( t ) = λ i e l i ( t ) + j N i a i j e l i ( t ) e l j ( t ) .

3. Proposed Controller

Here, a synchronous sliding mode controller is proposed. Figure 2 shows the proposed control structure. It achieves robustness against model uncertainty and disturbance. The chattering is attenuated by the use of a low pass filter (LPF).
To achieve synchronization, each UAV uses tracking errors of its neighbors to compute a sliding surface in the coupled error space. The sliding surface at the i-th UAV for the l axis is defined as
s l i ( t ) = e ¨ l i c ( t ) + k d e ˙ l i c ( t ) + k p e l i c ( t ) .
As usual for sliding mode controllers, it is shown in the next subsection that s l i ( t ) converges to zero in finite time, and maintains equal to zero thereafter. On the sliding surface, that is, when s l i ( t ) = 0 , the coupled error behaves according to the linear system
e ¨ l i c ( t ) + k d e ˙ l i c ( t ) + k p e l i c ( t ) = 0 ,
which has all poles in the left plane and, thereafter, is exponentially asymptotically stable for project parameters k d , k p > 0 .
The proposed control law for i-th UAV is
τ l i ( t ) = τ l i s ( t ) + τ l i f ( t ) ,
where τ l i s ( t ) and τ l i f ( t ) are, respectively, a smooth signal and a filtered signal of the control law, computed by
τ l i s ( t ) = p ¨ l i d ( t ) k d e ˙ l i ( t ) k p e l i ( t ) ,
τ ˙ l i f ( t ) + ξ i τ l i f ( t ) = u l i ( t ) ,
u l i ( t ) = sign ( s l i ( t ) ) η .
Equation (35) defines a low pass filter with cutoff frequency ξ i > 0 that converts a chattering signal u l i ( t ) to a smooth signal τ l i f ( t ) . The parameter η must be chosen by the designer to guarantee the stability of the overall system.
The proposed control law given by Equations (33)–(36) contains only information from the virtual leader (or from a broadcasting non-virtual leader), from the own i-th UAV, and from its neighborhood N i . The neighborhood information is contained in s l i ( t ) , defined in Equation (31), which is a function of e l i c ( t ) from Equation (30), which is a function of the own local error e l i ( t ) and the neighborhood errors e l j ( t ) , j N i .
Remark 1.
The variables k p and k d define the natural frequency and damping factor of the 2nd order local sliding surface s l i of the i-th UAV from Equation (31). As can be seen in [20], these gains also define a control bandwidth, which must be sufficiently small to account for, for example, to actuator dynamics. Since it is chosen the same gain k p and the same gain k d to all UAVs, it means that they have sliding surfaces that share the same control bandwidth. This is reasonable if all UAVs have similar physical, actuator, and aerodynamic characteristics. However, if there are distinct UAVs, the constants must be chosen to respect the control bandwidth of the UAV with the slowest dynamics.

3.1. Disturbance Model

Measurement or computation errors and the effect of non-modeled dynamics are incorporated in the dynamics model, given by Equation (12), as a disturbance signal described in the reference frame, d i = [ d x i d y i d z i ] T . It is supposed that the controller has no access to d i but there are known upper bounds Δ x i , Δ y i and Δ z i on the magnitude of the components of d i and upper bounds Δ ˜ x i , Δ ˜ y i , and Δ ˜ z i on the derivatives of the components of d i , that is,
| d l i ( t ) | Δ l i , | d ˙ l i ( t ) | Δ ˜ l i , l = { x , y , z } .
These upper bounds are used to define the value of η in Equation (36), as explained in Section 3.2. As a contribution of this paper is shown that the upper bounds on the components in the reference frame coordinates can be computed from the upper bounds δ t i , δ θ i , and δ ψ i on the components of the disturbance signal in the wind frame b i ( t ) ,
| b t i ( t ) | δ t i , | b θ i ( t ) | δ θ i , | b ψ i ( t ) | δ ψ i ,
and from the upper bounds δ ˜ t i , δ ˜ θ i and δ ˜ ψ i for the
| b ˙ t i ( t ) | δ ˜ t i , | b ˙ θ i ( t ) | δ ˜ θ i , | b ˙ ψ i ( t ) | δ ˜ ψ i .
The wind frame components of the disturbances are more naturally obtained, for example, in description of imprecision in the computation of drag or thrust forces. Assume that there is an upper bound Ω i for the i-th UAV angular velocity ω i and define the bounds vectors δ i [ δ t i δ θ i δ ψ i ] T and δ ˜ i [ δ ˜ t i δ ˜ θ i δ ˜ ψ i ] T . From Equation (11), it can be seen that
| d l i ( t ) | d i ( t ) = R i ( t ) b i ( t ) = b i ( t ) δ i .
The upper bounds of each component of d i are
Δ x i = Δ y i = Δ z i = δ i .
Since Equation (11) involves two frames in which one rotates related to the other, its derivative is obtained by using the Theorem of Coriolis [29]
d i ˙ ( t ) = R i ( t ) b ˙ i ( t ) + ω i ( t ) × b i ( t ) ,
where d i ˙ ( t ) contains two components. The first, b ˙ i ( t ) , is the derivative of the disturbance b i ( t ) , as seen by the wind frame. The second, ω i ( t ) × b i ( t ) , is generated by the rotation of the wind frame related to the inertial frame. See that a constant disturbance in the wind frame is a varying disturbance in the inertial frame, because of its rotation. Finally, R i ( t ) is used to represent the sum of these components in the inertial frame.
For the bounds δ i , δ ˜ i , and Ω i , it is obtained
| d ˙ l i ( t ) | d i ˙ ( t ) b ˙ i ( t ) + ω i ( t ) × b i ( t ) δ ˜ i + Ω i δ i .
In this way,
Δ ˜ x i = Δ ˜ y i = Δ ˜ z i = δ ˜ i + Ω i δ i .
Equations (40) and (44) provide the upper bounds to the proposed controller.

3.2. Stability Proof

To analyze the overall fleet behavior, all local variables must be concatenated in vectors. Concatenating the positions p i , virtual control inputs τ i ( t ) , and disturbances d i ( t ) from all UAVs of the fleet results in respectively P ( t ) = [ p 1 T ( t ) p n T ( t ) ] T , τ ( t ) = [ τ 1 T ( t ) τ n T ( t ) ] T , and D ( t ) = [ d 1 T ( t ) d n T ( t ) ] T , all R 3 n vectors. In this way, the dynamics of the fleet of UAVs is given by concatenating Equation (29) as
P ¨ ( t ) = τ + D ( t ) .
Similarly, the error and coupled error in x axis are R n vectors given by E ( t ) = [ e 1 T ( t ) e n T ( t ) ] T and E c ( t ) = [ e 1 c T ( t ) e n c T ( t ) ] T which are related by
E c ( t ) = ( H I 3 ) E ( t ) ,
where ⊗ denotes the Kronecker product and matrix H is given by Equation (24). The concatenation of the n UAVs sliding surfaces S ( t ) = [ s 1 T ( t ) s n T ( t ) ] T is obtained as
S ( t ) = E ¨ c ( t ) + k d E ˙ c ( t ) + k p E c ( t ) = ( H I 3 ) E ¨ ( t ) + k d E ˙ ( t ) + k p E ( t ) .
The proposed sliding mode control law is written as
τ ( t ) = τ s ( t ) + τ f ( t ) ,
where τ s ( t ) and τ f ( t ) are computed by
τ s ( t ) = P ¨ d ( t ) k d E ˙ ( t ) k p E ( t ) ,
τ ˙ f ( t ) + Ξ τ f ( t ) = U ( t ) ,
U ( t ) = diag η | s l i ( t ) | S ( t ) ,
with U ( t ) = [ u 1 T ( t ) u n T ( t ) ] T , u i ( t ) = [ u x i ( t ) u y i ( t ) u z i ( t ) ] T , and Ξ diag ( [ ξ 1 ξ n ] ) I 3 R 3 n × 3 n .
To analyze the fleet stability, the following Lyapunov functional candidate is proposed
V ( t ) = 1 2 S T ( t ) ( H I 3 ) 1 S ( t ) .
Note that, since H and H I 3 are a positive definite matrix, H 1 and ( H I 3 ) 1 are also a positive definite matrix, so V ( t ) is always positive for S ( t ) 0 .
By using Equations (45), (48) and (49), the sliding surface given by Equation (47) can be rewritten as
S ( t ) = ( H I 3 ) P ¨ ( t ) P ¨ d ( t ) + k d E ˙ ( t ) + k p E ( t ) = ( H I 3 ) τ s ( t ) + τ f ( t ) + D ( t ) P ¨ d ( t ) + k d E ˙ ( t ) + k p E ( t ) = ( H I 3 ) D ( t ) + τ f ( t ) .
Since ( H I 3 ) 1 is constant, the derivative of Equation (52) is
V ˙ ( t ) = S T ( t ) ( H I 3 ) 1 S ˙ ( t ) .
By deriving Equation (53) and after using Equation (50), V ˙ ( t ) is rewritten to
V ˙ ( t ) = S T ( t ) D ˙ ( t ) + τ ˙ f ( t ) = S T ( t ) D ˙ ( t ) S T ( t ) Ξ τ f ( t ) + S T ( t ) U ( t ) = i = 1 n s i T ( t ) d ˙ i ( t ) ξ i s i T ( t ) τ i f ( t ) + s i T ( t ) u i = i = 1 n l = { x , y , z } ( s l i d ˙ l i ξ i s l i τ l i f | s l i | η ) i = 1 n l = { x , y , z } | s l i | ( | d ˙ l i | + ξ i | τ l i f | η ) .
The upper bounds of the disturbance and its derivative are given, respectively, by Δ l i | d l i ( t ) | and Δ ˜ l i | d ˙ l i ( t ) | , which are computed by, respectively, Equations (41) and (44). It is shown in [24] that | τ l i f ( t ) | | d l i ( t ) | Δ l i . By using these upper bounds in Equation (55), it can be seen that
V ˙ ( t ) i = 1 n l = { x , y , z } | s l i | ( Δ ˜ l i + ξ i Δ l i η ) .
By choosing η satisfying
η Δ ˜ l i + ξ i Δ l i + ϵ , i { 1 , , n } , , l { x , y , z } ,
for some arbitrarily chosen constant ϵ > 0 , it is obtained
V ˙ ( t ) i = 1 n l = { x , y , z } | s l i ( t ) | ϵ = ϵ i = 1 n l = { x , y , z } | s l i ( t ) | = ϵ S ( t ) 1 ,
where S ( t ) 1 is the 1-norm of S ( t ) . Using the fact that the 1-norm is greater than the Euclidean norm of the same vector, then
V ˙ ( t ) ϵ S ( t ) ,
which means that V ( t ) and, therefore, S ( t ) go to zero in finite time [20]. On the sliding surface, the system behaves as a stable linear system given by Equation (32) and the error converges asymptotically to zero.
Remark 2.
Note that the sliding surface given by Equation (47), when rewritten in Equation (53), is a function only of the disturbance D ( t ) and the output of the filter τ f ( t ) . This has two main implications:
1. 
Since it is shown here that S ( t ) 0 , it follows that τ f ( t ) D ( t ) . In this way, τ f ( t ) estimates and compensates disturbances. Since the effect of airflow is not aligned to the fuselage is a disturbance, the presence of a disturbance compensation shows that the wind effect can be neglected in the initial model if this effect has known bounds.
2. 
If the disturbance is null at t = 0 , S ( 0 ) = 0 if τ f ( 0 ) = 0 and the system already starts in sliding condition. Similarly, if the known disturbance upper bound is relatively small, the system starts near the sliding surface and converges fast to the sliding surface.

4. Simulation

In this section, a simulation is made to show the effectiveness of the proposed controller. A scenario of 5 UAVs with communication links described by Figure 3 is used.
The matrices
L = 2 1 1 0 0 1 3 1 1 0 1 1 3 0 1 0 1 0 1 0 0 0 1 0 1
and Λ = I 5 are chosen to give the same weight for the UAV own error and for each of its relative errors. The choice k p = 0.5 and k d = 0.0625 provide a critically damped sliding surface with natural frequency ω n = 0.25 rad/s. These gains are chosen relatively small, as a way to limit the maximum commanded acceleration, even if the UAVs are initially far from their desired position. The low pass filters are settled such that Ξ = I 5 I 3 .
A fleet with a non-rectilinear 3D trajectory is described, which is defined by the virtual leader path given by
p x 0 ( t ) = 80 + 45 t [ m ] , p y 0 ( t ) = 20 cos ( 0.1 t ) [ m ] , γ 0 ( t ) = π 36 rad , ( z 0 ( 0 ) = 100 m ) .
For easy visualization, a time-varying formation is considered, whose horizontal projection in the reference frame has a V-shape and the altitude has time-varying oscillation. Accordingly, the formation rotation matrix R r is defined as R χ from Equation (16) and the clearance vectors p ˜ i r ( t ) related to the virtual leader are
p ˜ 1 r ( t ) = 0 0 10 sin ( 0.1 t ) , p ˜ 2 r ( t ) = 40 40 10 sin ( 0.1 t + 2 π / 5 ) , p ˜ 3 r ( t ) = 40 40 10 sin ( 0.1 t + 4 π / 5 ) , p ˜ 4 r ( t ) = 80 80 10 sin ( 0.1 t + 6 π / 5 ) , p ˜ 5 r ( t ) = 80 80 10 sin ( 0.1 t + 8 π / 5 ) .
The initial position of each UAV is defined as
p 1 ( 0 ) = 60 0 100 , p 2 ( 0 ) = 20 30 100 , p 3 ( 0 ) = 50 20 100 , p 4 ( 0 ) = 10 50 100 , p 5 ( 0 ) = 20 80 100 .
The initial velocity of each UAV is defined as
p ˙ 1 ( 0 ) = 50 5 0 , p ˙ 2 ( 0 ) = 40 10 0 , p ˙ 3 ( 0 ) = 45 10 0 , p ˙ 4 ( 0 ) = 40 5 0 , p ˙ 5 ( 0 ) = 45 0 0 .
The disturbance is simulated as
b i ( t ) = 0.2 cos ( 0.5 t ) cos ( 0.5 t ) cos ( 0.5 t ) T , i { 1 , 2 , 3 , 4 , 5 } .
From Equation (65), the magnitude of the upper bound vector δ i of b i ( t ) is computed as δ i = 0.35 . The magnitude of the upper bound vector of δ ˜ i is computed also from Equation (65) as δ ˜ i = 0.17 .
The upper bound of each component of d i ( t ) is computed by Equation (41) resulting in Δ x i = Δ y i = Δ z i = 0.35 . By simulation experiments it is verified that Ω i = 0.17 rad/s is an upper bound for the angular velocity amplitude; the upper bound in d ˙ i ( t ) is computed by Equation (44), resulting in Δ ˜ x i = Δ ˜ y i = Δ ˜ z i = 0.23 . By choosing ϵ = 0.42 , it is obtained from Equation (57) that η = 1 .
The system is implemented using an ode4 Runge-Kutta solver, with a fixed-step size of 1 ms. Since it is impossible to perfectly simulate the effect of a chattering input signal in a continuous differential equation, the controller output is evaluated at 10 ms time steps and maintained constant between time intervals.
For comparison purposes, the unfiltered synchronous formation flight controller presented in Reference [18] is also simulated. It is configured to be as similar as possible to the proposed controller. The first order sliding surface is defined with the same natural frequency as the proposed controller, that is, ω n = 0.25 rad/s. By using the same upper bound Δ x i = Δ y i = Δ z i = 0.35 and by choosing the same ϵ = 0.42 , it is computed η = 0.77 . Other parameters are exactly the same as the proposed controller.
Figure 4 shows the desired trajectory for each UAV in black, and the trajectory achieved by each UAV in distinct colors. Square and ‘*’ markers show respectively the desired and achieved positions in specific and equally spaced time instants. When a ‘*’ is inside the square, the UAV is in its desired position.
Figure 5 shows the formation flight error components e x i , e y i , and e z i for each i-th UAV for both controllers. Figure 6, shows the coupled error of each i-th UAV, which is given by Equation (46) for both controllers. It can be seen that, for both controllers, the system rapidly enters in sliding mode, the coupled errors slide in the prescribed linear sliding surface and achieve the performance described by the linear system that defines the sliding surface. It can also be seen that the error converges to zero, which shows that both controllers completely compensate for the added input disturbance.
Figure 7 shows the controller output τ x i , τ y i , and τ z i for each i-th UAV, which is generated by adding the smooth τ i s control signal and τ i f , obtained by filtering the chattering signal U i in the proposed controller, or is the unfiltered control signal in the controller from Reference [18]. As can be seen, the proposed control output is smooth, whereas the control output from the unfiltered SMC chatters.

5. Conclusions

A decentralized architecture for synchronous formation flight of UAVs based on sliding mode control with a low pass filter was proposed. The use of the SMC technique provides robustness to disturbances, in a way that the system slides in the prescribed sliding surface even in the presence of disturbances. The LPF virtually removes the chattering while maintaining the convergence to a null error in steady-state. In the proposed architecture only the chattering component of the control signal is filtered. As a result, the controller has a simpler expression when compared to recent results of the literature, such as in [23]. Also, it is presented an equation that is used to compute the upper bounds in the disturbance and in its derivative to a formation described in a global frame. This equation assumes that the upper bounds are known in the wind frame of each follower UAV. It is proved that the proposed controller is stable, achieving a prescribed sliding surface in finite time.
For future work, more realistic models for UAV and wind gusts can be implemented. Also, it is desired to implement other SOSMC, such as presented in References [21,22], in the context of the synchronous formation flight.

Author Contributions

Conceptualization, J.Y.I. and H.C.F.; Formal analysis, T.F.K.C. and J.Y.I.; Investigation, T.F.K.C.; Methodology, T.F.K.C., J.Y.I. and H.C.F.; Software, T.F.K.C.; Supervision, H.C.F.; Visualization, T.F.K.C.; Writing—original draft, T.F.K.C., J.Y.I. and H.C.F. All authors have read and agree to the published version of the manuscript.

Funding

This work was supported in part by National Council for Scientific and Technological Development—CNPq, under grants 312528/2017-5/PQ and 460311/2014-0 and by District Federal Research Foundation—FAPDF, under grant 0193.001153/2015.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zheng, Z.; Qian, M.; Li, P.; Yi, H. Distributed Adaptive Control for UAV Formation with Input Saturation and Actuator Fault. IEEE Access 2019, 7, 144638–144647. [Google Scholar] [CrossRef]
  2. Cai, Z.; Wang, L.; Zhao, J.; Wu, K.; Wang, Y. Virtual target guidance-based distributed model predictive control for formation control of multiple UAVs. Chin. J. Aeronaut. 2019. [Google Scholar] [CrossRef]
  3. Zhao, H.; Wu, S.; Wen, Y.; Liu, W.; Wu, X. Modeling and Flight Experiments for Swarms of High Dynamic UAVs: A Stochastic Configuration Control System with Multiplicative Noises. Sensors 2019, 19, 3278. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Liao, F.; Teo, R.; Wang, J.L.; Dong, X.; Lin, F.; Peng, K. Distributed Formation and Reconfiguration Control of VTOL UAVs. IEEE Trans. Control Syst. Technol. 2017, 25, 270–277. [Google Scholar] [CrossRef]
  5. Park, C.; Cho, N.; Lee, K.; Kim, Y. Formation Flight of Multiple UAVs via Onboard Sensor Information Sharing. Sensors 2015, 15, 17397–17419. [Google Scholar] [CrossRef]
  6. He, L.; Bai, P.; Liang, X.; Zhang, J.; Wang, W. Feedback formation control of UAV swarm with multiple implicit leaders. Aerosp. Sci. Technol. 2018, 72, 327–334. [Google Scholar] [CrossRef]
  7. Brandão, A.S.; Sarcinelli-Filho, M. On the Guidance of Multiple UAV using a Centralized Formation Control Scheme and Delaunay Triangulation. J. Intell. Robot. Syst. 2016, 84, 397–413. [Google Scholar] [CrossRef]
  8. Singh, S.N.; Zhang, R.; Chandler, P.; Banda, S. Decentralized nonlinear robust control of UAVs in close formation. Int. J. Robust Nonlinear Control 2003, 13, 1057–1078. [Google Scholar] [CrossRef]
  9. Rezaee, H.; Abdollahi, F.; Talebi, H.A. H based motion synchronization in formation flight with delayed communications. IEEE Trans. Ind. Electron. 2014, 61, 6175–6182. [Google Scholar] [CrossRef]
  10. Lin, W. Distributed UAV formation control using differential game approach. Aerosp. Sci. Technol. 2014, 35, 54–62. [Google Scholar] [CrossRef]
  11. Li, Z.; Xing, X.; Yu, J. Decentralized output-feedback formation control of multiple 3-DOF laboratory helicopters. J. Frankl. Inst. 2015, 352, 3827–3842. [Google Scholar] [CrossRef]
  12. Ren, W. Consensus strategies for cooperative control of vehicle formations. IET Control Theory A 2007, 1, 505–512. [Google Scholar] [CrossRef]
  13. Oh, K.K.; Park, M.C.; Ahn, H.S. A survey of multi-agent formation control. Automatica 2015, 53, 424–440. [Google Scholar] [CrossRef]
  14. Cordeiro, T.F.K.; Ferreira, H.C.; Ishihara, J.Y. Non linear controller and path planner algorithm for an autonomous variable shape formation flight. In Proceedings of the 2017 International Conference on Unmanned Aircraft Systems (ICUAS), Miami, FL, USA, 13–16 June 2017. [Google Scholar] [CrossRef]
  15. Liu, W.; Gao, Z. A distributed flocking control strategy for UAV groups. Comput. Commun. 2020, 153, 95–101. [Google Scholar] [CrossRef]
  16. Qiu, H.; Duan, H. A multi-objective pigeon-inspired optimization approach to UAV distributed flocking among obstacles. Inf. Sci. 2020, 509, 515–529. [Google Scholar] [CrossRef]
  17. Campa, G.; Gu, Y.; Seanor, B.; Napolitano, M.R.; Pollini, L.; Fravolini, M.L. Design and flight-testing of non-linear formation control laws. Control Eng. Pract. 2007, 15, 1077–1092. [Google Scholar] [CrossRef]
  18. Cordeiro, T.F.K.; Ferreira, H.C.; Ishihara, J.Y. Robust and Synchronous Nonlinear Controller for Autonomous Formation Flight of Fixed Wing UASs. In Proceedings of the 2019 International Conference on Unmanned Aircraft Systems (ICUAS), Atlanta, GA, USA, 11–14 June 2019. [Google Scholar] [CrossRef]
  19. Yu, J.; Dong, X.; Li, Q.; Ren, Z. Time-varying formation tracking for high-order multi-agent systems with switching topologies and a leader of bounded unknown input. J. Frankl. Inst. 2018, 355, 2808–2825. [Google Scholar] [CrossRef]
  20. Slotine, J.J.; Li, W. Applied Nonlinear Control; Prentice Hall: Englewood Cliffs, NJ, USA, 1990. [Google Scholar]
  21. Mei, K.; Ding, S. Second-order sliding mode controller design subject to an upper-triangular structure. IEEE Trans. Syst. Man Cybern. Syst. 2018. [Google Scholar] [CrossRef]
  22. Ding, S.; Park, J.H.; Chen, C.C. Second-order sliding mode controller design with output constraint. Automatica 2020, 112, 108704. [Google Scholar] [CrossRef]
  23. Zhao, D.; Li, C.; Zhu, Q. Low-pass-filter-based position synchronization sliding mode control for multiple robotic manipulator systems. Proc. Inst. Mech. Eng. I J. Syst. Control Eng. 2011, 225, 1136–1148. [Google Scholar] [CrossRef]
  24. Chong, S.; Sheng, Y.; Xiangyuan, Z.; Liu, X. An improved chattering-free sliding mode control with finite time convergence for reentry vehicle. In Proceedings of the 2016 IEEE Chinese Guidance, Navigation and Control Conference (CGNCC), Nanjing, China, 12–14 August 2016. [Google Scholar] [CrossRef]
  25. Balamurugan, S.; Venkatesh, P.; Varatharajan, M. Fuzzy sliding-mode control with low pass filter to reduce chattering effect: An experimental validation on Quanser SRIP. Sādhanā 2017, 42, 1693–1703. [Google Scholar] [CrossRef] [Green Version]
  26. Xu, B.; Li, J.; Yang, Y.; Wu, H.; Postolache, O. A Novel Sliding Mode Control with Low-Pass Filter for Nonlinear Handling Chain System in Container Ports. Complexity 2020, 2020, 7254503. [Google Scholar] [CrossRef]
  27. Han, T.; Guan, Z.H.; Wu, Y.; Zheng, D.F.; Zhang, X.H.; Xiao, J.W. Three-dimensional containment control for multiple unmanned aerial vehicles. J. Frankl. Inst. 2016, 353, 2929–2942. [Google Scholar] [CrossRef]
  28. Anderson, J.D., Jr. Introduction to Flight, 3rd ed.; Mcgraw-Hill: New York, NY, USA, 1988. [Google Scholar]
  29. Stevens, B.L.; Lewis, F.L. Aircraft Control and Simulation; Wiley-Interscience: New York, NY, USA, 2003. [Google Scholar]
Figure 1. The i-th Unmanned Aerial Vehicle (UAV), its force and velocity vectors, and attitude angles.
Figure 1. The i-th Unmanned Aerial Vehicle (UAV), its force and velocity vectors, and attitude angles.
Sensors 20 03094 g001
Figure 2. Block diagram of the control structure.
Figure 2. Block diagram of the control structure.
Sensors 20 03094 g002
Figure 3. Five UAVs and their undirected communication links. The virtual leader is not shown here. All UAVs have access to the virtual leader’s trajectory information.
Figure 3. Five UAVs and their undirected communication links. The virtual leader is not shown here. All UAVs have access to the virtual leader’s trajectory information.
Sensors 20 03094 g003
Figure 4. Desired trajectory and UAV position.
Figure 4. Desired trajectory and UAV position.
Sensors 20 03094 g004
Figure 5. Position errors for all UAVs.
Figure 5. Position errors for all UAVs.
Sensors 20 03094 g005
Figure 6. Coupled errors for all UAVs.
Figure 6. Coupled errors for all UAVs.
Sensors 20 03094 g006
Figure 7. The controller output for all UAVs.
Figure 7. The controller output for all UAVs.
Sensors 20 03094 g007

Share and Cite

MDPI and ACS Style

Cordeiro, T.F.K.; Ishihara, J.Y.; Ferreira, H.C. A Decentralized Low-Chattering Sliding Mode Formation Flight Controller for a Swarm of UAVs. Sensors 2020, 20, 3094. https://doi.org/10.3390/s20113094

AMA Style

Cordeiro TFK, Ishihara JY, Ferreira HC. A Decentralized Low-Chattering Sliding Mode Formation Flight Controller for a Swarm of UAVs. Sensors. 2020; 20(11):3094. https://doi.org/10.3390/s20113094

Chicago/Turabian Style

Cordeiro, Thiago F. K., João Y. Ishihara, and Henrique C. Ferreira. 2020. "A Decentralized Low-Chattering Sliding Mode Formation Flight Controller for a Swarm of UAVs" Sensors 20, no. 11: 3094. https://doi.org/10.3390/s20113094

APA Style

Cordeiro, T. F. K., Ishihara, J. Y., & Ferreira, H. C. (2020). A Decentralized Low-Chattering Sliding Mode Formation Flight Controller for a Swarm of UAVs. Sensors, 20(11), 3094. https://doi.org/10.3390/s20113094

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