Next Article in Journal
A Conflict Resolution Strategy at a Taxiway Intersection by Combining a Monte Carlo Tree Search with Prior Knowledge
Next Article in Special Issue
High-Speed Virtual Flight Testing Platform for Performance Evaluation of Pitch Maneuvers
Previous Article in Journal
Study on the Transition to True North in Air Navigation
Previous Article in Special Issue
Modeling and Application of Out-of-Cabin and Extra-Vehicular Dynamics of Airdrop System Based on Kane Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coordinated Control of Quadrotor Suspension Systems Based on Consistency Theory

1
College of Marine Electrical Engineering, Dalian Maritime University, Dalian 116026, China
2
Key Laboratory of Technology and System for Intelligent Ships of Liaoning Province, Dalian 116026, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Aerospace 2023, 10(11), 913; https://doi.org/10.3390/aerospace10110913
Submission received: 21 August 2023 / Revised: 16 October 2023 / Accepted: 23 October 2023 / Published: 26 October 2023
(This article belongs to the Special Issue Flight Dynamics, Control & Simulation)

Abstract

:
This paper designs a cooperative control method for the multi-quadrotor suspension system based on consistency theory and realizes the cooperative formation trajectory tracking control of the multi-quadrotor suspension system by designing a consistent formation cooperative algorithm of virtual piloting and a nonlinear controller. First, a new quadrotor suspension system model is established based on the traditional quadrotor model using the Newton–Euler method. This model can accurately reflect the influence of the load on the quadrotor while obtaining the swing of the load. Then, the vertical and horizontal positions are designed separately based on the quadrotor motion characteristics, and the formation algorithm based on the virtual pilot consistency theory ensures that the final convergence of each position is consistent. An integral backstepping controller and an integral backstepping sliding mode controller are designed for quadrotor position, attitude, and load swing control to achieve accurate and fast quadrotor trajectory tracking control while reducing load swing. The stability of all the controllers is demonstrated using Lyapunov functions. Finally, a multi-quadrotor suspension system formation cooperative simulation experiment is designed to verify the designed control method.

1. Introduction

The quadrotor suspension system is a versatile technology with diverse applications [1,2], including freight transport [3,4,5,6], agriculture, forestry, crop protection, and environmental sensing [7]. This system is highly regarded for its exceptional flight flexibility, cost effectiveness, and ease of loading and unloading [8,9,10,11]. Nevertheless, single quadrotor suspension systems have their inherent limitations, characterized by a constrained load capacity and diminished efficiency when dealing with expansive areas or multiple loads [12,13]. To surmount these limitations, the adoption of multi-quadrotor suspension systems is gaining momentum. Cooperative formation control has emerged as an enticing solution across various sectors, such as transportation and detection, capturing substantial interest among researchers [14,15,16,17,18]. This approach harnesses the combined capabilities of multiple quadrotors, elevating their performance and efficiency across an extensive spectrum of applications.
While there has been substantial research into the control of single quadrotor suspension systems [11,19,20,21,22,23,24,25,26], many of these methods exhibit limitations. Much of the existing research predominantly concentrates either on the load [19,20,21] or solely on the quadrotor [11,22,23]. Although this specialized approach ensures effective control of the individual component being studied, it often sacrifices the overall system’s dynamics. This can lead to undesirable outcomes such as horizontal position oscillations of the quadrotor or significant swings in the suspended load [24,25,26]. To address these issues comprehensively, there is a pressing need for integrated control strategies that simultaneously consider both the quadrotor and the suspended load. This approach aims to optimize the performance and stability of the entire system.
Multi-quadrotor suspension systems operating in coordinated formations offer several advantages, including enhanced mission efficiency, extended operational reach, and the capacity to establish networks for task execution [27]. Various formation control methods govern these systems, each with its unique strengths and approaches. Prominent among these methods are the Leader–Follower Method [28], which designates a leader quadrotor to guide the formation, albeit with potential limitations on flexibility. The Virtual Structure Method [29] introduces a virtual framework within the formation, allowing for the maintenance of desired relative positions while affording greater formation geometry flexibility. The Behavior-Based Method [30] empowers each quadrotor to exhibit specific behaviors within the formation, enhancing adaptability to dynamic environments. The Graph Theory-Based Method [31] leverages graph theory to model communication relationships among quadrotors, transforming formation control challenges into consistency and system stability problems. In recent years, there has been a notable upsurge in research centered on graph theory-based methods. These techniques, by representing communication relationships as edges in a graph, provide a potent framework for modeling multi-quadrotor systems. Through matrix transformations, they facilitate the conversion of formation control challenges into consistency and system stability problems, thereby enhancing the precision and adaptability of multi-quadrotor formation control strategies.
Consistency theory, firmly rooted in graph theory-based methodologies, stands as a fundamental concept within the domain of multi-quadrotor formation control. Its primary objective revolves around the assurance that quadrotors, operating within a formation, maintain a steadfast and unchanging relative spatial configuration as they navigate in perfect harmony and synchrony [32]. This theory assumes a central role in safeguarding the intended formation shape and preserving the relative positions of quadrotors. Consequently, it fortifies the overarching unity and stability of the formation. The bedrock of consistency theory lies in the cultivation of synchronized motion and seamless coordination among quadrotors. As a result, it profoundly enhances the efficiency and effectiveness of multi-quadrotor systems across a diverse spectrum of applications [33].
Reference [34] introduces a novel consistent formation control algorithm based on the super-twisting algorithm, notable for its robustness and effectiveness in addressing parameter uncertainties. Reference [35] combines sliding mode control and consistency theory to suppress time-varying disturbances in the position loop, ensuring that quadrotors adhere closely to desired trajectories while maintaining a fixed formation. Reference [36] tackles consistent formation control in discrete time, employing system identification techniques to obtain dynamic models for altitude, horizontal, and yaw motions of quadrotors. A consistent formation controller is designed using the Riccati equation and parameter uncertainty theorem to ensure the stable flight of multi-quadrotor formations. In contrast to the extensive research on consistent formation control of quadrotors, the field of formation flight involving quadrotor suspension systems remains relatively unexplored. Recognizing the quadrotor as the central element of the system, we propose the application of established consistent formation control techniques to quadrotor suspension systems. This extension opens new possibilities for achieving precise and coordinated control within multi-quadrotor suspension setups.
Compared to traditional quadrotor formation flight control, achieving formation with a quadrotor suspension system introduces heightened complexity and challenges [37,38,39]. This complexity arises from various factors, including the mutual influence between the horizontal position of the quadrotor and the oscillations of the suspended load. The introduction of two additional degrees of freedom associated with the load’s swing angles necessitates innovative underdrive system control schemes. Furthermore, effective control must manage load oscillations to prevent collisions between adjacent loads. In cooperative formation control, the unique characteristics of cable-suspended loads mandate separate handling of the horizontal direction and longitudinal height. In summary, achieving cooperative formation control in quadrotor suspension systems poses a formidable challenge due to its highly coupled, under-driven, and complex nonlinear nature. Innovative control strategies and approaches are imperative to address these challenges effectively.
Recognizing the intricate nature of quadrotor suspension systems and building upon the foundational insights from previous research that have underpinned the control of individual quadrotor suspension systems [40,41,42], and drawing inspiration from the field of quadrotor formation control, this paper introduces an innovative cooperative control methodology firmly rooted in the principles of consistency theory. The primary contributions and unique aspects of this study are outlined as follows:
  • Formation Control with a Virtual Leader: In response to the intricacies inherent to the operation of quadrotor suspension systems, this research introduces a novel formation controller, strategically incorporating a virtual leader. Deliberately designed to address both the horizontal and vertical aspects of the system, this controller places a primary emphasis on preserving the coherency of quadrotor relative positions. Notably, the design enhances the anticipated relative distances among quadrotors, effectively reducing the potential for collisions and, consequently, promoting a heightened degree of safety and synchronization within the formation.
  • Integrated Control for Position, Attitude, and Swing Angles: A sophisticated control strategy is meticulously devised to govern the position, attitude, and swing angles within the quadrotor suspension system. This intricate control paradigm harnesses the seamless integration of three pivotal components: the position controller, the swing angle controller, and the attitude controller. Notably, the control of swing angles in the suspension system is achieved by translating output signals from the position and swing angle controllers into inputs for the attitude controller. This holistic methodology ensures the precise tracking of quadrotor trajectories while mitigating the oscillations induced by load dynamics.
  • Integral Type Backstepping Sliding Mode Controller: To address the challenge posed by oscillations in the position of quadrotors induced by suspended loads, an integral-type backstepping sliding mode controller is introduced. This controller offers a resilient solution to mitigate oscillations through the skillful integration of integral control techniques within the backstepping sliding mode framework. An innovative adaptation involves substituting the conventional sign function with the hyperbolic tangent function. This alteration not only minimizes jitter in controller outputs, but also fortifies the overall robustness of the system, resulting in enhanced performance and reliability.
In summary, this paper introduces a comprehensive cooperative control framework for quadrotor suspension systems, grounded in consistency theory. The novel contributions include formation control with a virtual leader to maintain relative position coherence, an integrated control strategy covering position, attitude, and swing angles, and the implementation of an integral-type backstepping sliding mode controller with a hyperbolic tangent function for improved stability and robustness. Collectively, these innovations expand the frontiers of cooperative control for quadrotor suspension systems, tackling their inherent complexities, and laying the groundwork for more effective and reliable operations.
The specific chapters are organized as follows. Section 2 presents the dynamic modeling of the quadrotor suspension system. Section 3 shows the cooperative formation controller design and stability proof. Section 4 presents the experimental simulation results to verify the designed control scheme. Section 5 concludes the paper.

2. Quadrotor Suspension System Dynamic

The quadrotor suspension system includes a quadrotor body with six degrees of freedom and a suspended load with two degrees of freedom, which cannot simply be treated as the drive of the load or the load as a disturbance to the quadrotor, so a new quadrotor suspension system model must be established. The model is shown in Figure 1.
In Figure 1, a common cross-type quadrotor is shown, defining the inertial coordinate system O I X I , Y I , Z I . The center of mass of the quadrotor is considered as the origin of the quadrotor body coordinate system [5], x points in the direction of the tail of the aircraft, y points to the right of the aircraft, the z vertical x o y plane points above the aircraft body, α is the angle between the cable and the x o z plane, β is the angle between the cable and the y o z plane, and the cable length is L. The conversion of inertial and airframe coordinate systems is described by Euler angles, the rotation around the x-axis is the roll angle ϕ , the rotation around the y-axis is the pitch angle θ , the rotation around the z-axis is the yaw angle ψ . The rotation matrix is
R B I = c θ c ψ s ϕ s θ c ψ c ϕ s ψ s ϕ s ψ c ϕ s θ c ψ c θ s ψ s ϕ s θ s ψ + c ϕ c ψ s ϕ c ψ c ϕ s θ s ψ s θ s ϕ c θ c ϕ c θ ,
where s means sin and c means cos.
We define the position of the quadrotor as ξ = x , y , z T , the attitude of the quadrotor as Θ = ϕ , θ , ψ T , the swing angle of the load as δ = α , β T . The cable attachment point is considered to be the center of mass of the quadrotor, and the cable is considered to be tight at all times. Considering that the load is always under the quadrotor, that is, 90 < α , β < 90 .
The system is driven by a quadrotor motor, the rotor generates lift T i and reaction-torque Q i , i = 1 , 2 , 3 , 4 , and the control signal U B is
U B = U 1 U 2 U 3 U 4 = T 1 + T 2 + T 3 + T 4 T 3 T 4 T 2 T 1 Q 1 + Q 2 Q 3 Q 4 .
Since the load cable is assumed to be tied to the quadrotor’s center of mass, the load is assumed to have no effect on the quadrotor’s attitude angle, so we model the attitude dynamics using the Newton–Euler equation.
I Q ς ˙ + ς × ( I Q ς ) = M a ,
where I Q = d i a g ( I x , I y , I z ) is the rotational inertia and ς = p , q , r T is the attitude angular velocity. Using the principle of small-angle approximation, we let ς = Θ , M a = l U 2 , l U 3 ; U 4 T is external moments, l is the quadrotor arm length.
The positions of the quadrotor and the swing angle are coupled, so we use the Euler–Lagrange function to model the dynamics.
The position of the load ξ L = x L , y L , z L T can be obtained from the quadrotor position and the swing angle, with cable length
ξ L = ξ + L ρ ,
where ρ = ρ 1 ρ 2 ρ 3 = sin β sin α cos β cos α cos β .
The Euler–Lagrange function is as follows, and q = x , y , z , α , β T ,
L ( q , q ˙ ) = T ( q , q ˙ ) V ( q ) ,
where T ( q , q ˙ ) = T Q + T L is system kinetics composed of quadrotor kinetic energy T Q and load kinetic energy T L , V ( q ) is the potential energy of the system.
The kinetic energy of the system is
T ( q , q ˙ ) = 1 2 m Q ξ ˙ T ξ ˙ + 1 2 m L ξ ˙ L T ξ ˙ L ,
where m Q , m L are quadrotor and load mass.
System potential energy is
V ( q ) = m Q g z + m L g z L .
We obtain the Euler–Lagrange equations for system
d d t L ( q , q ˙ ) q ˙ L ( q , q ˙ ) q = F g ,
where F g = F p T , M s T T is the internal and external forces and moments on the system, and F p = R B I 0 , 0 , U 1 0 , 0 , ( m Q + m L ) g T , M s = 0 , 0 T . Substituting q = x , y , z , α , β T , we can obtain the quadrotor suspension system position and the swing angle model.
In summary, we can obtain the dynamics model of the quadrotor suspension system.
ϕ ¨ = l U 2 I x + I y I z I x θ ˙ ψ ˙ θ ¨ = l U 3 I y + I z I x I y ϕ ˙ ψ ˙ ψ ¨ = U 4 I z + I x I y I z ϕ ˙ θ ˙ x ¨ = a 1 U 1 m L L b 1 m Q + m L y ¨ = a 2 U 1 m L L b 2 m Q + m L z ¨ = a 3 U 1 m L L b 3 m Q + m L g α ¨ = c 1 ( cos α a 2 + sin α a 3 ) U 1 m Q L cos β β ¨ = c 2 ( cos β a 1 + sin β sin α a 2 + sin β cos α a 3 ) U 1 m Q L ,
where a 1 = sin ϕ sin ψ cos ϕ sin θ cos ψ a 2 = sin ϕ cos ψ cos ϕ sin θ sin ψ a 3 = cos ϕ cos θ , b 1 = 2 ρ 1 t 2 b 2 = 2 ρ 2 t 2 b 3 = 2 ρ 3 t 2 , c 1 = 2 m Q L sin β β ˙ α ˙ c 2 = m Q L sin β cos β 2 .

3. Cooperative Formation Controller Design

The cooperative formation controller in this paper is divided into a trajectory generation part based on consistency theory and a nonlinear controller part based on integral type backstepping sliding mode, and the control block diagram is shown in Figure 2.
In Figure 2, the desired trajectory for each quadrotor is computed using the first-order consistent formation system, employing virtual piloting. The control parameters u x , u y , u α , and u β are derived from both the position controller and the swing angle controller. Notably, for the attitude solution, the influence of roll and pitch on position is disregarded, that is, ϕ d = ( u y + u β ) cos ψ ( u x + u α ) sin ψ θ d = ( u y + u β ) sin ψ ( u x + u α ) cos ψ . The desired attitude is obtained, and the force and moment of the system are obtained by the attitude angle controller and position controller to realize the trajectory tracking control of the quadrotor suspension system.
The control structure designed in this paper has the advantage of converting the control quantity of the swing angle controller to the attitude control input, thus avoiding the problem of start failure caused by the conversion to the position control input, and at the same time realizing the accurate control of the quadrotor position and the swing reduction control of the suspended load.

3.1. Collaborative Trajectory Generation

In the context of ensuring the consistent formation control of multi-quadrotor suspension systems, a fundamental understanding of algebraic graph theory is indispensable. Let us define an oriented graph G = V , E , where the vertex set is denoted as V = v 1 , v 2 , , v n , representing a collection of n nodes, and the edge set is E V × V , signifying the adjacency relationships between these nodes. The index set for the nodes is denoted as I 1 , 2 , , n . To describe the communication relationships among the nodes, we introduce the adjacency matrix A = a i j R n × n . Here, the element a i j takes on the value > 0 to indicate that node j can receive information from node i, while a i j = 0 signifies the absence of such information flow. Furthermore, the Laplacian matrix of G is defined as L A = l i j R n × n , with diagonal elements l i i = j = 1 n a i j and off-diagonal elements l i j = a i j .
Theorem 1 
([43]). For a communication structure directed graph G consisting of n intelligent agents, a sufficient necessary condition for this system to converge to average consistency, which means that the states or behaviors of these agents collectively tend to reach a specific consistent state or value, is that graph G is a strongly connected equilibrium graph or a connected undirected graph.
The first-order consistent formation control algorithm with navigator designed in this paper, i.e., the consistency of the quadrotor positions, thus considers the agents as single-integral systems with the equation of state
ξ ˙ i = u i , i = 1 , 2 , , n ,
where ξ i R is the position state of agent i, and u i R is the control input.
Theorem 2 
([43]). In any initial state, the control algorithm u i enables the multi-agent system (10) to satisfy
lim t ξ i ξ j = 0 , i = 1 , 2 , , n .
Then, the multi-intelligent body system (10) is said to be able to achieve asymptotic consistency.
The design consistency formation algorithm is
u i = h ξ i ξ i * ξ 0 j = 1 n a i j ξ i ξ i * ξ j ξ j * ,
where h denotes the communication relationship between the virtual navigator and the follower, h = 1 is the ability of the follower to receive information from the navigator, ξ * R is the deviation matrix related to system formation, and ξ 0 R is the position state of the navigator. From Theorem 2, it can be seen that the designed algorithm can make the intelligent body agree while maintaining the relative distance from the navigator.
We propose a virtual pilot and six quadrotor followers. Each follower can obtain the information of the pilot, and the directed communication topology is shown in Figure 3.
The communication graph in Figure 3 is a strongly connected graph containing a directed spanning tree. The adjacency matrix is A = 0 1 1 1 1 1 1 0 0 1 0 0 0 0 0 0 0 1 0 1 0 0 1 0 0 0 1 0 0 0 1 0 0 0 0 0 0 0 0 1 0 1 0 0 0 1 0 0 0 , the Laplacian matrix is L = 6 1 1 1 1 1 1 0 1 1 0 0 0 0 0 0 2 1 0 1 0 0 1 0 2 0 1 0 0 0 1 0 1 0 0 0 0 0 0 1 2 1 0 0 0 1 0 0 1 . From Theorem 1, the system is consistent.

3.2. Controller Design

The controller design of the quadrotor suspension system is divided into three modules: quadrotor attitude, quadrotor position, and suspension swing angle, where attitude is the inner loop and position and suspension swing angle are the outer loop. In the following, k , λ , μ , h and ε are all constants greater than zero. First, the integral type backstepping controller is designed. We take the roll angle ϕ as an example.
Step 1: Define the error
e ϕ = ϕ d ϕ ,
where ϕ d is the desired roll angle, deriving Equation (13) and obtaining the error of the roll angle velocity
e ˙ ϕ = ϕ ˙ d ϕ ˙ = ϕ ˙ d w ϕ
where w ϕ is the virtual input volume.
Step 2: Design the Lyapunov function of e ϕ :
V 1 = 1 2 e ϕ 2 .
The derivative is obtained as
V ˙ 1 = e ϕ e ˙ ϕ = e ϕ ( ϕ ˙ d w ϕ ) .
Take w ϕ as the virtual control input, i.e., the actual roll speed, and define its expected value as w ϕ d to make V ˙ 1 0 ; let w ϕ d = k ϕ e ϕ + ϕ ˙ d .
Step 3: Add an integration term to the desired virtual input to increase robustness and reduce error:
w ϕ d = k ϕ e ϕ + ϕ ˙ d + λ ϕ χ ϕ ,
where χ ϕ = 0 t e ϕ ( τ ) d τ .
Step 4: Define the virtual input error
e ϕ ˙ = w ϕ d w ϕ .
Substituting in w ϕ d , we obtain
w ϕ = k ϕ e ϕ + ϕ ˙ d + λ ϕ χ ϕ e ϕ ˙ .
The derivative of (18) can be obtained as
e ˙ ϕ ˙ = w ˙ ϕ d w ˙ ϕ = k ϕ e ˙ ϕ + ϕ ¨ d + λ ϕ e ϕ ϕ ¨ ,
where ϕ ¨ is a roll angle acceleration model. It can be obtained from (19),
V ˙ 1 = e ϕ ( ϕ ˙ d w ϕ ) = k ϕ e ϕ 2 e ϕ λ ϕ χ ϕ + e ϕ e ϕ ˙ .
To make V ˙ 1 0 , it is also needed to make e ϕ ˙ , χ ϕ tend to zero.
Step 5: Design Lyapunov functions for e ϕ , e ϕ ˙ , χ ϕ :
V 2 = 1 2 e ϕ 2 + 1 2 e ϕ ˙ 2 + 1 2 λ ϕ χ ϕ 2 .
Differentiate (22)
V ˙ 2 = e ϕ e ˙ ϕ + e ϕ ˙ e ˙ ϕ ˙ + λ ϕ χ ϕ χ ˙ ϕ .
Substitute (20) and (21) into (23):
V ˙ 2 = k ϕ e ϕ 2 + e ϕ e ϕ ˙ + e ϕ ˙ k ϕ e ˙ ϕ + ϕ ¨ d + λ ϕ e ϕ l U 2 + θ ˙ ψ ˙ ( I y I z ) I x .
To make V ˙ 2 0 , design the control variable U 2 as
U 2 = I x l 1 k ϕ 2 + λ ϕ e ϕ + k ϕ + k ϕ ˙ e ϕ ˙ k ϕ λ ϕ χ ϕ + ϕ ¨ d θ ˙ ψ ˙ ( I y I z ) I x .
The proof of the stability of the integral type backstepping controller is as follows:
Proof. 
Taking U 2 into (24), we can obtain
V ˙ 2 = k ϕ e ϕ 2 k ϕ ˙ e ϕ ˙ 2 0 .
It can be concluded that the designed control law is asymptotically stable. □
Similarly, the control laws for pitch angle and yaw angle can be obtained, along with stability proof.
U 3 = I y l 1 k θ 2 + λ θ e θ + k θ + k θ ˙ e θ ˙ k θ λ θ χ θ + θ ¨ d ϕ ˙ ψ ˙ ( I z I x ) I y ,
U 4 = I z 1 k ψ 2 + λ ψ e ψ + k ψ + k ψ ˙ e ψ ˙ k ψ λ ψ χ ψ + ψ ¨ d ϕ ˙ θ ˙ ( I x I y ) I z ,
where θ d , ψ d represent the desired pitch angle and yaw angle, respectively, and
e θ = θ d θ , χ θ = 0 t e θ ( τ ) d τ e θ ˙ = k θ e θ + θ ˙ d + λ θ χ θ θ ˙ e ψ = ψ d ψ , χ ψ = 0 t e ψ ( τ ) d τ e ψ ˙ = k ψ e ψ + ψ ˙ d + λ ψ χ ψ ψ ˙ .
Similarly, the swing angle controller is designed in the same manner.
u α = m Q L cos β 1 k α 2 + λ α e α + α ¨ d + k α + k α ˙ e α ˙ k α λ α χ α + c 1 U 1 u β = m Q L 1 k β 2 + λ β e β + β ¨ d + k β + k β ˙ e β ˙ k β λ β χ β + c 2 U 1 ,
where α ¨ d and β ¨ d are the desired values, and the definitions of each variable that appear are as follows: e α = α d α , χ α = 0 t e α ( τ ) d τ e α ˙ = k α e α + α ˙ d + λ α χ α α ˙ s α = μ α e α + e α ˙ , e β = β d β , χ β = 0 t e β ( τ ) d τ e β ˙ = k β e β + β ˙ d + λ β χ β β ˙ s β = μ β e β + e β ˙ .
The position of the quadrotor is susceptible to disturbances and load. Therefore, an integral backstepping sliding mode controller is designed by combining sliding mode control with backstepping control to address this issue. Taking the altitude z of the quadrotor as an example, the first four steps are the same as mentioned above.
Step 1: Define the error term e z = z d z and the virtual input w z .
Step 2: Design a Lyapunov function V 1 = 1 2 e z 2 for the error term e z such that the desired virtual input is w z d = k z e z + z ˙ d .
Step 3: Add an integral term to the function w z d = k z e z + z ˙ d + λ z χ z , χ z = 0 t e z ( τ ) d τ .
Step 4: Define the error term for the virtual input e z ˙ = w z d w z = k z e z + z ˙ d + λ z χ z z ˙ .
Step 5: Based on the integral backstepping method, add sliding mode control and define the sliding surface.
s z = μ z e z + e z ˙ = μ z + k z e z + e ˙ z + λ z χ z .
Taking the derivative of s z ,
s ˙ z = μ z e ˙ z + e ˙ z ˙ = μ z + k z e z ˙ k z e z λ z χ z + λ z e z + z ¨ d a 3 U 1 m L L b 3 m Q + m L + g .
Step 6: Design Lyapunov functions for e z , χ z , s z :
V 2 = 1 2 e z 2 + 1 2 λ ϕ χ ϕ 2 + 1 2 s z 2 .
Taking the derivative of (32),
V ˙ 2 = e z e ˙ z + λ z χ z χ ˙ z + s z s ˙ z = k z e z 2 + e z e z ˙ + s z μ z + k z e z ˙ k z e z λ z χ z + λ z e z + z ¨ d a 3 U 1 m L L b 3 m Q + m L + g .
To make the controller stable and reduce chattering in sliding mode control, use the hyperbolic tangent function instead of the sign function and design the control law.
U 1 = m Q + m L μ z + k z e z ˙ k z e z λ z χ z + λ z e z + z ¨ d + g + h z s z + ε z tanh ( s z ) + m L L b 3 a 3 ;
ε z is the sliding mode switching gain.
The proof of the stability of the integral backstepping sliding mode controller is as follows:
Proof. 
Substituting (34) into (33),
V ˙ 2 = k z e z 2 + e z e z ˙ + s z h z s z ε z tanh ( s z ) k z e z 2 + e z e z ˙ h z s z 2 ε z s z .
We let Q be a positive definite matrix.
Q = k z + h z μ z 2 h z μ z 1 2 h z μ z 1 2 h z .
The determinant of Q is positive, which means the control parameters need to satisfy h z k z + μ z 1 4 > 0 . We can obtain
e z e z ˙ T k z + h z μ z 2 h z μ z 1 2 h z μ z 1 2 h z e z e z ˙ = k z e z 2 e z e z ˙ + h z s z 2 .
We substitute into (35)
V ˙ 2 e z e z ˙ T Q e z e z ˙ ε z s z 0 .
Therefore, it can be proved that the control system is stable. □
Similarly, horizontal position controllers and swing angle controllers can be designed. It should be noted that the horizontal position controller does not directly obtain the desired roll and pitch angles, but obtains the desired rotation matrix, that is, u x = a 1 u y = a 2 , to avoid the calculation of inverse sine and greatly reduce the computational complexity. The design process and stability analysis are not repeated here, and the results are obtained.
u x = m Q + m L μ x + k x e x ˙ k x e x λ x χ x + λ x e x + x ¨ d + h x s x + ε x tanh ( s x ) + m L L b 1 U 1 u y = m Q + m L μ y + k y e y ˙ k y e y λ y χ y + λ y e y + y ¨ d + h y s y + ε y tanh ( s y ) + m L L b 2 U 1 ,
where x ¨ d and y ¨ d are the desired values, and the definitions of each variable that appear are as follows: e x = x d x , χ x = 0 t e x ( τ ) d τ e x ˙ = k x e x + x ˙ d + λ x χ x x ˙ s x = μ x e x + e x ˙ , e y = y d y , χ y = 0 t e y ( τ ) d τ e y ˙ = k y e y + y ˙ d + λ y χ y y ˙ s y = μ y e y + e y ˙ .

4. Simulation Results

This paper employs Simulink simulation experiments to validate the proposed control scheme. The model parameters for the quadrotor suspension system are shown in Table 1.
The initial positions of the six quadrotors are ξ 0 = 0 1 1 2 0 2 0 1 1 2 2 2 0 0 0 0 0 0 T , with position deviation matrix ξ * = 0 1.5 1.5 3 0 3 2 3 3 / 2 3 / 2 3 3 3 0 0 0 0 0 0 T . The trajectory of the virtual leader is 0 , t > 0 0 , t > 0 3 , t > 0 , 0.3 ( t 10 ) , t > 10 2 sin 0.3 ( t 10 ) , t > 10 3 , t > 10 , with the first 10 s being the rising stage. The control system parameters are shown in Table 2.
The following are the results of the simulation experiments.
The expected trajectory of the quadrotor generated by the consensus formation algorithm shown in Figure 4. It can be observed from the figure that the designed formation algorithm can effectively complete the formation task under the guidance of the virtual leader, generating a formation trajectory that brings the quadrotors’ horizontal and vertical positions from the initial position to the expected formation trajectory. Additionally, the algorithm forms the corresponding formation according to the designed deviation matrix and maintains consistency among the quadrotors.
Figure 5, Figure 6 and Figure 7 illustrate the flight trajectories of six quadrotor suspension systems controlled by the designed method, backstepping control, and PID control, respectively. It is evident that all three methods enable the quadrotors to reach the desired altitude and successfully track their respective planned formation paths. Consequently, they achieve coordinated control of the multi-quadrotor suspension system. Nevertheless, variations in position control accuracy exist among these methods.
Figure 8 and Figure 9 illustrate the position errors of the quadrotors using the three control methods. It is clear that, when compared to the control without the swing angle controller, the method proposed in this paper achieves stability within 5 s, while the backstepping method takes 10 s, and PID control requires 15 s. Furthermore, the approach outlined in this paper exhibits a smaller overshoot. Once stability is reached, the positions of the quadrotors are influenced by load swinging due to the absence of a swing angle controller. As a result, periodic errors persist even after stability is achieved, and this issue remains unresolved with the backstepping method, while PID control induces more significant oscillations. However, the design presented in this paper significantly mitigates the impact of load swinging, preventing the repetitive oscillation of quadrotor positions.
Figure 10 illustrates the swinging of the suspended load under the three control methods. It is evident that, with the swing angle controller in place, the suspended load experiences smaller swings when the quadrotor initiates horizontal movement and quickly achieves stability. Minimal swings occur as the quadrotor’s position undergoes periodic motion, with the absence of high-frequency oscillations. This design effectively mitigates the risks associated with excessive and frequent load swings.
Figure 11 shows the control of the quadrotor’s attitude angles. It can be observed that the control method designed in this paper can accurately and quickly track the desired attitude angles.

5. Conclusions

To solve the problem of coordinated control of quadrotor suspension systems, this paper proposes a nonlinear control method based on a consistency algorithm and backstepping sliding mode control. An accurate eight-degree-of-freedom dynamic model of the quadrotor suspension system is established using the Lagrangian equation. A consistency formation algorithm with a virtual leader is designed to plan the formation flight trajectory of multiple quadrotor suspension systems. To address the issue of the quadrotor position being sensitive to load swinging, an integral backstepping sliding mode controller is designed for the quadrotor position channel, while integral backstepping controllers are designed for the quadrotor attitude and load swinging. By converting the output values of the outer-loop position and swinging controllers into the inner-loop attitude reference values, the control method not only achieves quadrotor trajectory tracking control, but also suppresses load swinging. Finally, simulation experiments using Simulink are conducted to verify the effectiveness of the proposed control method, which outperforms the backstepping control method in terms of control accuracy and load swinging suppression.

Author Contributions

Conceptualization, G.W. and X.C.; methodology, X.C.; software, X.C.; validation, D.M.; formal analysis, Y.F.; investigation, X.C.; resources, G.W.; data curation, D.M.; writing—original draft preparation, X.C.; writing—review and editing, X.C.; visualization, D.M.; supervision, D.M.; project administration, G.W.; funding acquisition, Y.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “the Nature Science Foundation of China” (grand number 61976033, 51609033), “the Nature Science Foundation of Liaoning Province of China” (grand number 20180520005), “the Fundamental Research Funds for the Central Universities” (grand numbers 3132021106, 3132020110), “the Liaoning Province Doctoral Research Start-up” (grand number 2022-BS-094), “the China Postdoctoral Science Foundation” (grand number 2022M710569) and “the postgraduate education and teaching reform project of Dalian Maritime University” (grand number YJG2022713).

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Acknowledgments

This work was supported by “the Nature Science Foundation of China” (grand number61976033, 51609033), “the Nature Science Foundation of Liaoning Province of China” (grand number 20180520005), “the Fundamental Research Funds for the Central Universities” (grand numbers 3132021106, 3132020110), “the Liaoning Province Doctoral Research Start-up” (grand number 2022-BS-094), “the China Postdoctoral Science Foundation” (grand number 2022M710569) and “the postgraduate education and teaching reform project of Dalian Maritime University” (grand number YJG2022713).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jiang, B.; Li, B.; Zhou, W.; Lo, L.-Y.; Chen, C.-K.; Wen, C.-Y. Neural Network Based Model Predictive Control for a Quadrotor UAV. Aerospace 2022, 9, 460. [Google Scholar] [CrossRef]
  2. Mariani, M.; Fiori, S. Design and Simulation of a Neuroevolutionary Controller for a Quadcopter Drone. Aerospace 2023, 10, 418. [Google Scholar] [CrossRef]
  3. Potter, J.J.; Adams, C.J.; Singhose, W. A planar experimental remote-controlled helicopter with a suspended load. IEEE/ASME Trans. Mechatronics 2015, 20, 2496–2503. [Google Scholar] [CrossRef]
  4. Duan, D.; Wang, Z.; Li, J.; Zhang, C.; Wang, Q. Stabilization control for unmanned helicopter-slung load system based on active disturbance rejection control and improved sliding mode control. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2021, 235, 1803–1816. [Google Scholar] [CrossRef]
  5. Yang, S.; Xian, B.; Cai, J.; Wang, G. Finite-time convergence control for a quadrotor unmanned aerial vehicle with a slung load. IEEE Trans. Ind. Inform. 2023; early access. [Google Scholar]
  6. Outeiro, P.; Cardeira, C.; Oliveira, P. Control Architecture for a Quadrotor Transporting a Cable-Suspended Load of Uncertain Mass. Drones 2023, 7, 201. [Google Scholar] [CrossRef]
  7. Si, P.; Fu, Z.; Shu, L.; Yang, Y.; Huang, K.; Liu, Y. Target-barrier coverage improvement in an insecticidal lamps internet of UAVs. IEEE Trans. Veh. Technol. 2022, 71, 4373–4382. [Google Scholar] [CrossRef]
  8. Sreenath, K.; Lee, T.; Kumar, V. Geometric control and differential flatness of a quadrotor UAV with a cable-suspended load. In Proceedings of the 52nd IEEE Conference on Decision and Control, Firenze, Italy, 10–13 December 2013; pp. 2269–2274. [Google Scholar]
  9. Vahdanipour, M.; Khodabandeh, M. Adaptive fractional order sliding mode control for a quadrotor with a varying load. Aerosp. Sci. Technol. 2019, 86, 737–747. [Google Scholar] [CrossRef]
  10. Wang, Z.; Zhao, T. Based on robust sliding mode and linear active disturbance rejection control for attitude of quadrotor load UAV. Nonlinear Dyn. 2022, 108, 3485–3503. [Google Scholar] [CrossRef]
  11. Roy, K.R.; Waghmare, L.M.; Patre, B.M. Dynamic modeling and displacement control for differential flatness of quadrotor UAV slung-load system. Int. J. Dyn. Control. 2023, 11, 637–655. [Google Scholar] [CrossRef]
  12. Yang, P.; Zhang, A.; Bi, W.; Li, M. Cooperative group formation control for multiple quadrotors system with finite-and fixed-time convergence. ISA Trans. 2023, 138, 186–196. [Google Scholar] [CrossRef]
  13. Doakhan, M.; Kabganian, M.; Azimi, A. Robust adaptive control for formation-based cooperative transportation of a payload by multi quadrotors. Eur. J. Control. 2023, 69, 100763. [Google Scholar] [CrossRef]
  14. Zhang, K.; Yang, Z.; Başar, T. Multi-agent reinforcement learning: A selective overview of theories and algorithms. In Handbook of Reinforcement Learning and Control; Springer: Berlin/Heidelberg, Germany, 2021; pp. 321–384. [Google Scholar]
  15. Oroojlooy, A.; Hajinezhad, D. A review of cooperative multi-agent deep reinforcement learning. Appl. Intell. 2023, 53, 13677–13722. [Google Scholar] [CrossRef]
  16. Du, Z.; Negenborn, R.R.; Reppa, V. Cooperative multi-agent control for autonomous ship towing under environmental disturbances. IEEE/CAA J. Autom. Sin. 2021, 8, 1365–1379. [Google Scholar] [CrossRef]
  17. Mukras, S.M.S.; Omar, H.M. Development of a 6-DOF testing platform for multirotor flying vehicles with suspended loads. Aerospace 2021, 8, 355. [Google Scholar] [CrossRef]
  18. Ding, F.; Sun, C.; Ai, Y.; Huang, J. Sliding Mode Control for Quadrotor-Slung Load Transportation System with State Constraints. In Proceedings of the 2022 IEEE International Conference on Cyborg and Bionic Systems (CBS), Wuhan, China, 24–26 March 2023; pp. 368–373. [Google Scholar]
  19. Yu, G.; Cabecinhas, D.; Cunha, R.; Silvestre, C. Aggressive maneuvers for a quadrotor-slung-load system through fast trajectory generation and tracking. Auton. Robot. 2022, 46, 499–513. [Google Scholar] [CrossRef]
  20. Lv, Z.; Zhao, Q.; Li, S.; Wu, Y. Finite-time control design for a quadrotor transporting a slung load. Control. Eng. Pract. 2022, 122, 105082. [Google Scholar] [CrossRef]
  21. Lv, Z.; Wu, Y.; Sun, X.-M.; Wang, Q.-G. Fixed-time control for a quadrotor with a cable-suspended load. IEEE Trans. Intell. Transp. Syst. 2022, 23, 21932–21943. [Google Scholar] [CrossRef]
  22. Salih, Z.; Saleh, M.H. Attitude and Altitude Control of Quadrotor Carrying a Suspended Payload using Genetic Algorithm. J. Eng. 2022, 28, 25–40. [Google Scholar] [CrossRef]
  23. Sun, H.; Gu, X.; Luo, S.; Liang, Y.; Bai, J. Robust stabilization technique for a quadrotor slung-load system using sliding mode control. J. Phys. Conf. Ser. 2022, 2232, 012013. [Google Scholar] [CrossRef]
  24. Chandra, A.; Lal, P.P.S. Higher Order Sliding Mode Controller for a Quadrotor UAV with a Suspended Load. IFAC-PapersOnLine 2022, 55, 610–615. [Google Scholar] [CrossRef]
  25. Wang, J.; Yuan, X.; Zhu, B. Geometric control for trajectory-tracking of a quadrotor UAV with suspended load. IET Control. Theory Appl. 2022, 16, 1271–1281. [Google Scholar] [CrossRef]
  26. Omar, H.M.; Mukras, S. Integrating anti-swing controller with px4 autopilot for quadrotor with suspended load. J. Mech. Sci. Technol. 2022, 36, 1511–1519. [Google Scholar] [CrossRef]
  27. Yan, D.; Zhang, W.; Chen, H. Design of a multi-constraint formation controller based on improved MPC and consensus for quadrotors. Aerospace 2022, 9, 94. [Google Scholar] [CrossRef]
  28. Wang, P.K.C. Navigation strategies for multiple autonomous mobile robots moving in formation. J. Robot. Syst. 1991, 8, 177–195. [Google Scholar] [CrossRef]
  29. Lewis, M.A.; Tan, K.H. High precision formation control of mobile robots using virtual structures. Auton. Robot. 1997, 4, 387–403. [Google Scholar] [CrossRef]
  30. Balch, T.; Arkin, R.C. Behavior-based formation control for multirobot teams. IEEE Trans. Robot. Autom. 1998, 14, 926–939. [Google Scholar] [CrossRef]
  31. Desai, J.P.; Ostrowski, J.P.; Kumar, V. Modeling and control of formations of nonholonomic mobile robots. IEEE Trans. Robot. Autom. 2001, 17, 905–908. [Google Scholar] [CrossRef]
  32. Liu, R.; Qu, B.; Wei, T.; Zhang, L.; Yan, L.; Chai, X. Research on UAV Formation Obstacle Avoidance Based on Consistency Control. In Proceedings of the 2023 IEEE 12th Data Driven Control and Learning Systems Conference (DDCLS), Xiangtan, China, 12–14 May 2023; pp. 155–160. [Google Scholar]
  33. Yu, H.; Ning, L. Coordinated Obstacle Avoidance of Multi-AUV Based on Improved Artificial Potential Field Method and Consistency Protocol. J. Mar. Sci. Eng. 2023, 11, 1157. [Google Scholar] [CrossRef]
  34. Rojo-Rodriguez, E.G.; Garcia, O.; Ollervides, E.J.; Zambrano-Robledo, P.; Espinoza-Quesada, E.S. Robust consensus-based formation flight for multiple quadrotors. J. Intell. Robot. Syst. 2019, 93, 213–226. [Google Scholar] [CrossRef]
  35. Wang, J.; Ma, X.; Li, H.; Tian, B. Self-triggered sliding mode control for distributed formation of multiple quadrotors. J. Frankl. Inst. 2020, 357, 12223–12240. [Google Scholar] [CrossRef]
  36. Steinleitner, A.; Ballam, R.; McFadyen, A. Practical consensus-based formation control for quadrotor systems. In Proceedings of the 2022 International Conference on Unmanned Aircraft Systems (ICUAS), Dubrovnik, Croatia, 21–24 June 2022; pp. 1510–1519. [Google Scholar]
  37. Liu, Y.; Li, Y. Application of Inverse Optimal Formation Control for Euler-Lagrange Systems. IEEE Trans. Intell. Transp. Syst. 2023, 24, 5655–5662. [Google Scholar] [CrossRef]
  38. Hwang, C.L.; Abebe, H.B. Generalized and heterogeneous nonlinear dynamic multiagent systems using online RNN-based finite-time formation tracking control and application to transportation systems. IEEE Trans. Intell. Transp. Syst. 2021, 23, 13708–13720. [Google Scholar] [CrossRef]
  39. Omar, H.M.; Akram, R.; Mukras, S.M.S.; Mahvouz, A.A. Recent advances and challenges in controlling quadrotors with suspended loads. Alex. Eng. J. 2023, 63, 253–270. [Google Scholar] [CrossRef]
  40. Chen, X.; Zhao, Y.; Fan, Y. Adaptive Integral Backstepping Control for a Quadrotor with Suspended Flight. In Proceedings of the 2020 5th International Conference on Automation, Control and Robotics Engineering (CACRE), Dalian, China, 19–20 September 2020; 2020; pp. 226–234. [Google Scholar]
  41. Fan, Y.; Guo, H.; Han, X.; Chen, X. Research and Verification of Trajectory Tracking Control of a Quadrotor Carrying a Load. Appl. Sci. 2022, 12, 1036. [Google Scholar] [CrossRef]
  42. Fan, Y.; Chen, X.; Zhao, Y.; Song, B. Nonlinear control of quadrotor suspension system based on extended state observer. Acta Autom. Sin. 2023, 49, 1758–1770. [Google Scholar]
  43. Olfati-Saber, R.; Fax, J.A.; Murray, R.M. Consensus and cooperation in networked multi-agent systems. Proc. IEEE 2007, 95, 215–233. [Google Scholar] [CrossRef]
Figure 1. Quadrotor suspension system model.
Figure 1. Quadrotor suspension system model.
Aerospace 10 00913 g001
Figure 2. Cooperative formation and backstepping sliding mode control structure.
Figure 2. Cooperative formation and backstepping sliding mode control structure.
Aerospace 10 00913 g002
Figure 3. Quadrotor suspension system communication topology.
Figure 3. Quadrotor suspension system communication topology.
Aerospace 10 00913 g003
Figure 4. Consistent formation-generated quadrotor desired trajectories.
Figure 4. Consistent formation-generated quadrotor desired trajectories.
Aerospace 10 00913 g004
Figure 5. 3D trajectories under proposed method.
Figure 5. 3D trajectories under proposed method.
Aerospace 10 00913 g005
Figure 6. 3D trajectories through Backstepping controller without swing angle controller.
Figure 6. 3D trajectories through Backstepping controller without swing angle controller.
Aerospace 10 00913 g006
Figure 7. 3D trajectories through PID controller without swing angle controller.
Figure 7. 3D trajectories through PID controller without swing angle controller.
Aerospace 10 00913 g007
Figure 8. Position control error of 6 quadrotors throughout the process.
Figure 8. Position control error of 6 quadrotors throughout the process.
Aerospace 10 00913 g008
Figure 9. Control error of 6 quadrotors after 20 s.
Figure 9. Control error of 6 quadrotors after 20 s.
Aerospace 10 00913 g009
Figure 10. Comparison of swing angles under different control methods.
Figure 10. Comparison of swing angles under different control methods.
Aerospace 10 00913 g010
Figure 11. Desired and real attitude angle control.
Figure 11. Desired and real attitude angle control.
Aerospace 10 00913 g011
Table 1. Model Parameters for the Quadrotor Suspension System.
Table 1. Model Parameters for the Quadrotor Suspension System.
ParameterValue
m Q 1.8 kg
m L 0.2 kg
l 0.2 m
I x 0.03 kg·m 2
I y 0.03 kg·m 2
I z 0.04 kg·m 2
g 9.81 m/(s 2 )
L 0.3 m
Table 2. Control System Model Parameters.
Table 2. Control System Model Parameters.
ParameterValue
k ϕ , θ , k ψ , k x , y , k z , k α , β 53.3, 5.8, 1.45, 7.8, 1
k ϕ ˙ , θ ˙ , k ψ ˙ , k α ˙ , β ˙ 5.5, 3.5, 2
λ ϕ , θ , λ ψ , λ x , y , λ z , λ α , β 8.5, 0.01, 0.015, 20.55, 0.1
μ x , y , μ z 0.5, 20
h x , y , h z 0.8, 1
ε x , y , ε z 2, 5
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, X.; Fan, Y.; Wang, G.; Mu, D. Coordinated Control of Quadrotor Suspension Systems Based on Consistency Theory. Aerospace 2023, 10, 913. https://doi.org/10.3390/aerospace10110913

AMA Style

Chen X, Fan Y, Wang G, Mu D. Coordinated Control of Quadrotor Suspension Systems Based on Consistency Theory. Aerospace. 2023; 10(11):913. https://doi.org/10.3390/aerospace10110913

Chicago/Turabian Style

Chen, Xinyu, Yunsheng Fan, Guofeng Wang, and Dongdong Mu. 2023. "Coordinated Control of Quadrotor Suspension Systems Based on Consistency Theory" Aerospace 10, no. 11: 913. https://doi.org/10.3390/aerospace10110913

APA Style

Chen, X., Fan, Y., Wang, G., & Mu, D. (2023). Coordinated Control of Quadrotor Suspension Systems Based on Consistency Theory. Aerospace, 10(11), 913. https://doi.org/10.3390/aerospace10110913

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