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],
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
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.
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
As usual for sliding mode controllers, it is shown in the next subsection that
converges to zero in finite time, and maintains equal to zero thereafter. On the sliding surface, that is, when
, the coupled error behaves according to the linear system
which has all poles in the left plane and, thereafter, is exponentially asymptotically stable for project parameters
.
The proposed control law for
i-th UAV is
where
and
are, respectively, a smooth signal and a filtered signal of the control law, computed by
Equation (
35) defines a low pass filter with cutoff frequency
that converts a chattering signal
to a smooth signal
. 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
. The neighborhood information is contained in
, defined in Equation (
31), which is a function of
from Equation (
30), which is a function of the own local error
and the neighborhood errors
,
.
Remark 1. The variables and define the natural frequency and damping factor of the 2nd order local sliding surface 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 and the same gain 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,
. It is supposed that the controller has no access to
but there are known upper bounds
,
and
on the magnitude of the components of
and upper bounds
,
, and
on the derivatives of the components of
, that is,
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
,
, and
on the components of the disturbance signal in the wind frame
,
and from the upper bounds
,
and
for the
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
for the
i-th UAV angular velocity
and define the bounds vectors
and
. From Equation (
11), it can be seen that
The upper bounds of each component of
are
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]
where
contains two components. The first,
, is the derivative of the disturbance
, as seen by the wind frame. The second,
, 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,
is used to represent the sum of these components in the inertial frame.
For the bounds
,
, and
, it is obtained
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
, virtual control inputs
, and disturbances
from all UAVs of the fleet results in respectively
,
, and
, all
vectors. In this way, the dynamics of the fleet of UAVs is given by concatenating Equation (
29) as
Similarly, the error and coupled error in
x axis are
vectors given by
and
which are related by
where ⊗ denotes the Kronecker product and matrix
is given by Equation (
24). The concatenation of the
n UAVs sliding surfaces
is obtained as
The proposed sliding mode control law is written as
where
and
are computed by
with
,
, and
.
To analyze the fleet stability, the following Lyapunov functional candidate is proposed
Note that, since and are a positive definite matrix, and are also a positive definite matrix, so is always positive for .
By using Equations (
45), (
48) and (
49), the sliding surface given by Equation (
47) can be rewritten as
Since
is constant, the derivative of Equation (
52) is
By deriving Equation (
53) and after using Equation (
50),
is rewritten to
The upper bounds of the disturbance and its derivative are given, respectively, by
and
, which are computed by, respectively, Equations (
41) and (
44). It is shown in [
24] that
. By using these upper bounds in Equation (
55), it can be seen that
By choosing
satisfying
for some arbitrarily chosen constant
, it is obtained
where
is the 1-norm of
. Using the fact that the 1-norm is greater than the Euclidean norm of the same vector, then
which means that
and, therefore,
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 and the output of the filter . This has two main implications: - 1.
Since it is shown here that , it follows that . In this way, 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 , if 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
and
are chosen to give the same weight for the UAV own error and for each of its relative errors. The choice
and
provide a critically damped sliding surface with natural frequency
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
.
A fleet with a non-rectilinear 3D trajectory is described, which is defined by the virtual leader path given by
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
is defined as
from Equation (
16) and the clearance vectors
related to the virtual leader are
The initial position of each UAV is defined as
The initial velocity of each UAV is defined as
The disturbance is simulated as
From Equation (
65), the magnitude of the upper bound vector
of
is computed as
. The magnitude of the upper bound vector of
is computed also from Equation (
65) as
.
The upper bound of each component of
is computed by Equation (
41) resulting in
. By simulation experiments it is verified that
rad/s is an upper bound for the angular velocity amplitude; the upper bound in
is computed by Equation (
44), resulting in
. By choosing
, it is obtained from Equation (
57) that
.
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,
rad/s. By using the same upper bound
and by choosing the same
, it is computed
. 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
,
, and
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
,
, and
for each
i-th UAV, which is generated by adding the smooth
control signal and
, obtained by filtering the chattering signal
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.