Next Article in Journal
Sand Net Device to Control the Meanders of a Coastal River: The Case of the Authie Estuary (France)
Previous Article in Journal
Maritime Cyber Security: A Global Challenge Tackled through Distinct Regional Approaches
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distributed Swarm Control Algorithm of Multiple Unmanned Surface Vehicles Based on Grouping Method

College of Intelligent Systems Science and Engineering, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2021, 9(12), 1324; https://doi.org/10.3390/jmse9121324
Submission received: 30 October 2021 / Revised: 15 November 2021 / Accepted: 21 November 2021 / Published: 23 November 2021
(This article belongs to the Section Ocean Engineering)

Abstract

:
This paper addresses the distributed swarm control problem of multiple unmanned surface vehicles (USVs) in Euclidean space with virtual leader. Firstly, to investigate the topology of the neighborhood relations between vehicles, a new time-variant topology structure is proposed. Secondly, to research the dynamic properties of the group for the case where the number of virtual leader is different, a grouping method based on cosine similarity is proposed. Thirdly, to ensure the high effeciency of information transmission and the reduction of costs, a distributed swarm control algorithm is proposed, which is mainly composed of three parts: gradient descent term, velocity consensus term and navigational feedback term. To analyze the stability of system, the concept of translation framework is introduced. Based on the properties of the Hamiltonian and LaSalle invariance principle, the stability of multiple USVs swarm motion is proved. Finally, simulation results illustrate the effectiveness of the proposed methods.

1. Introduction

Swarming is ubiquition in nature, and it is a form of collective behaviour of multiple interacting agents with a common group objective. Swarm is defined as multiple autonomous agents moving cooperatively to fulfil global objective of a scientific or technological mission. This term is often observed in nature. Each autonomous agent is modelled as a particle and characterised by its position and a function describing its dynamics [1,2,3]. Understanding the mechanisms and operational principles in them can provide useful ideas for developing formation control, distributed cooperative control and coordination of multiple mobile autonomous agents/robots [4].
Over the past several years, the theory and application of multi-agent systems have attracted great attention. Specifically, the multi-agent distributed control not only has received much attention, but also has made significant research progress on account of its extensive applications [5,6,7,8]. Multiple USVs, which is one of the most actively studied topics within the realm of multi-agents systems, generally aims to drive USVs to achieve prescribed constraints on their states [9]. The studies focus on designing appropriate protocols and algorithms so that the multiple USVs can achieve various forms of agreement.
Keeping a swarming group connected is a connectivity problem. In representative literatures [10,11] always establishes group stabilization by assuming that the vehicles are connected through a time-invariant topology. In [12], the vessels are assumed interconnected through a directed topology rather than bidirectional. The connectivity of the group topology relies on the actual distances between vehicles and may be destroyed during group evolution. Efforts should be made to preserve connectivity of the topology while achieving the desired flocking objective. The disadvantage of time-invariant topology is that requiring the number of multiple USVs is fixed. For example, a group vehicles disperse when facing multi-task, and a flock splits when encountering predators, the topology needs to be redesigned.
The second concern is the communication constraints of these strategies. In recent years, a leader-follower system was investigated in many literatures [11,13,14]. A theoretical framework for the design and analysis of flocking algorithm is proposed in [15]. However, these swarming algorithms proposed in thses literatures require substantial communication energy consumptions during the flocking process, it is impossible to apply this algorithm to engineering applications such as distributed sensing in mobile sensor networks. Under the circumstances, the information flow cause substantial increase in the load, which bring enormous pressure on the server. Multiple virtual leaders become another choice, which can reduce the communication loads and realize the information diffluence. Considering that only partial agents can be informed of the virtual leaders, according to the extent of social distancing, a group composed of N mobile agents is divided into multiple different subgroups in [16]. The complex flocking of multiple inertial agents is analyzed via the decomposition approach in [17]. In [18], an algorithm solving the flocking problem for multiple groups of mobile agents in heterogeneous networks is studied. To the best of our knowledge, few results have discussed the method for a group of multiple USVs to flocking in any specified number of groups by cosine similarity method.
The third part focus on the swarm control algorithm. With the deepening research on multi-USVs system, many algorithms have been proposed in literatures [19,20,21,22]. In [19,20], graph theory has been used to investigate the linear consensus problem. Local attractive/repulsive potentials have been used to define the interactive forces between neighboring agents to deal with the separation and cohesion problems in [21,22]. In flocking control, the potential function approach is used to keep a swarming group cohesive and separate. Furthermore, consensus algorithms have been studies extensively in the context of aooperative control of multi-agent systems [23,24,25]. Both the sonsensus algorithm and the potentional function approach obtain information via neighbor-to-neighbor communication. In [26], swarming robots make their motion decisions based on partial information of the entire flocking network. If all members have global information or know the desired trajectory, the cohesive force generated by this global information can keep the flocking group connected. There are also many researches on flocking algorithm. Such as, semi-flocking algorithm [27], local flocking algorithm [28], a individual-based alignment/repulsion algorithm [29], etc. Meanwhile, ref. [11] theoretically proved that provably convergent flocking, formation control, and path following can be integrated in a single control architecture. In the above literature, the formation controller or flocking algorithm use consensus method and distributed approaches for the navigation of multi-USVs. Consensus algorithms and their variants are applied to formation control and swarm control [23]. This means that there is no essential difference between formation control and swarm control algorithm. These algorithms are relevant to the one proposed in this note in the sense that the emergent flocking phenomenon results form the alignment mechanism.
To sum up what is discussed above, firstly, motivated by the importance of network connectivity in multiple USVs, a time-variant topology is designed. It is dependent on proximity nets with interaction range of USVs cluster and α - g r i d s with a class of USV satisfying some algebraic constraints. Secondly, to reduce the communication load and realize the information diffluence, a grouping method for a group of multiple USVs to flocking in any specified number of groups based on cosine similarity is proposed. Thirdly, swarming-based algorithms have several advantages that make them suitable for use in sensor management. First, they are completely distributed algorithms. They are highly compatible with the distributed nature of sensor management in sensor networks. Second, in swarming-based algorithms, each vehicle needs to communicate only with its neighbors. Using swarming-based algorithms for sensor management requires only local communication between sensors. Third, in flocking-based algorithms, vehicles apply simple flocking rules, using this type of algorithms for sensor management has low computation overhead for the sensors. Furthermore, they are highly flexible and scalable. To ensure the high effeciency of information transmission and the reduction of costs, a distributed swarm control algorithm is proposed. Furthermore, the concept of translation framework is introduced. The stability of multiple USVs swarm motion is proved by the Hamiltonian and LaSalle invariance principle.
In summary, the main contributions of this paper are summarized as follows.
  • To achieve multiple USVs connectivity, a time-variant topology is designed, which is dependent on proximity nets and α - g r i d s . The connectivity of vehicles are bidirectional.
  • To reduce the communication loads and realize the information diffluence, a grouping method for a group of multiple USVs to flocking in any specified number of groups based on cosine similarity is proposed.
  • To ensure the high effeciency of information transmission and the reduction of costs, a distributed swarm control algorithm based on gradient descent, velocity consensus term and navigational feedback is proposed.
This paper is organized as follows. A table of notations and some used variables in the paper is presented (Table 1). Section 2 describes problem description and some necessary preliminaries. Section 3 designs the model of USVs, virtual leader, and grouping method. Section 4 proposes the distributed swarm control algorithm and proves the stability of the USVs swarming. Section 5 gives the simulation results and comparison results to verify effectiveness of the proposed algortihm. Section 6 concludes this paper.

2. Problem Description and Preliminaries

2.1. Problem Description

Multiple USVs system is a complex system. In this paper, we investigate the dynamic properties of the group for the case where the number of virtual leader is different and the topology of the neighborhood relations between vehicles is dynamic. Under the premise of learning preparatory knowledge, a grouping method based on cosine similarity is proposed. It is grouped according to the number of virtual leaders. Then, a disributed swarm control algorithm is propsoed. It is mainly composed of gradient descent term, velocity consensus term and navigational feedback term. The structure of the proposed grouping method and distributed swarm control algorithm for multiple USVs system is shown in Figure 1.

2.2. Algebraic Graph Theory

A directed graph G = ( V , E ) consists of a vertex set V = { 1 , 2 , , n + 1 } and the set of edges E { V × V \ ( i , i ) } . The graph is in general directed and has no self-loops, it means that E { ( i , j ) : i , j V , j i } .
The adjacency matrix A = [ a i j ] n × n of a graph is called to be undirected if ( i , j ) E ( j , i ) E (i.e., a i j = a j i = 1 ); otherwise, A is directed. The adjacency matrix A is a matrix with nonzero elements, in which a i j 0 ( i , j ) E . For an undirected graph G , The adjacency matrix A is symmetric, A T = A . The Laplacian matrix L associates with the graph G is define as L = D A where D = d i a g { d 1 , d 2 , , d n } with d i = j = 1 n a i j . Laplacian matrix L always has a right eigenvector of 1 n = ( 1 , 1 , , 1 ) T associated with eigenvalue λ 1 = 0 . It satisfies the following assumption.
Assumption 1.
In the graph theory, there is at least one path from the virtual leader navigator to each vehicle.

2.3. Proximity Nets

In low-speed applications, the damping can be assumed to be linear. Its quadratic velocity term can be ignored. Ship models are usually reduced-order models for control of the horizontal palne motions (surge, sway and yaw).
The set of neighbors of ith USV is described as
N i = { j V : ( i , j ) E } = { j V : a i j 0 } .
Define the interaction radius r > 0 , it is the interaction range between any two USVs. As shown in Figure 2, an open-loop circle neighborhood is defined as [30]
N i = { j V : X j X i < r } ,
where · represents the Euclidean norm in R 2 , X i : = ( x i , y i , ψ i ) T R 3 denotes the position and yaw angle of ith USV for all i V , X i is the center, r is the radius.
The set of edges E X is denoted as
E X = { ( i , j ) V × V : X j X i < r , i j } .
The vector X = c o l ( X 1 , X 2 , , X n ) R 3 × n is called the configuration of all USVs of the graph G . A framework is a pair ( G , X ) consists of a graph G and the configuration X of its USVs. The proximity net G ( X ) = ( V , E X ) can be defined by V and E X . The framework ( G ( X ) , X ) is called a proximity structure.
Assumption 2.
The interaction radius r of all USVs is the same. The topology of the USV wireless sensor network is a proximity net G with a interaction range r. All proximity nets are bidirectional graphs.

2.4. Geometric Structure: α - g r i d s

In order to describe the movement process of multiple USVs swarm more clearly, we use g r i d to represent the geometry of all USVs swarm.
First, we find a set of USVs with the following properties form all USVs, where each USV is equally distanced form all of its neighbors on a proximity net G ( X ) . This kind of USV can be described as a set of solutions satisfying the following algebraic constraints:
X j X i = d , j N i ( X ) .
To describe the set of solutions X of another class different form X, which are very close to the set of solutions satisfying Formula (4), the set of inequalities is shown as
δ X j X i d δ , ( i , j ) E X ,
and refer to these solutions satisfying the above inequalities constraints as quasi α - g r i d . The α - g r i d and quasi α - g r i d are shown in Figure 3 and Figure 4, respectively.
Next, to quantitatively describe the different degrees of configuration X and α - g r i d , we introduce a variable: deviation energy function (DEF). It is defined as
E ( X ) = 1 ( | E X | + 1 ) i = 1 n j N i ω ( X j X i d ) ,
where ω ( z ) = z 2 is called potential energy function (PEF).
For the swarm with n USVs, the DEF can be regarded as a nonsmooth PEF. The α - g r i d s are the solutions with PEF equal to zero. Moreover, the α - g r i d s are also the global minimum value of the E ( x ) .
For the quasi α - g r i d X with an side length uncertainty of δ , the deviation energy function is defined as
E ( x ) | E x | ( | E x | + 1 ) δ 2 δ 2 = ε 2 d 2 , ε 1 .

2.5. Nonnegative Mapping Function: λ -Norm

In order to construct the smooth collective potential energy function (CPEF) of the USVs swarm and the planar adjacency matrix of proximity structure, we define a concept of nonnegative mapping function: λ -norm.
The λ -norm of a vector is defined as
z λ = 1 σ [ 1 + σ z 2 1 ] ,
where σ > 0 is a constant, λ -norm is a map R 3 R 0 .
The gradient λ σ ( z ) = z λ is given as
λ σ ( z ) = z 1 + σ z 2 = z 1 + σ z λ .
The constant σ of the λ -norm remains fixed in this paper. The mapping z λ is differentiable everywhere, but it is not differentiable at z = 0 . Since the model of multi-USVs swarm is discrete, we can construct a smooth continuous CPEF of USVs swarm by using the property of λ -norm.
The scalar function ρ h ( z ) is a convex function, which has the following properties: it is continuous and smooth in the interval ( 0 , 1 ) . In this paper, we use convex function to construct smooth PEF with finite cut-off and smooth adjacency matrix. The convex function is as follows [15]:
ρ h ( z ) = 1 , z [ 0 , h ) 1 2 [ 1 + c o s ( z h 1 h π ) ] , z [ h , 1 ] 0 , o t h e r w i s e ,
where h ( 0 , 1 ) .
Through the convex function ρ h ( z ) , we define a concept of planar adjacency matrix A ( X ) . The elements in A ( X ) are defined as
a i j ( X ) = ρ h ( X j X i λ r α ) [ 0 , 1 ] , j i ,
where r α = r λ and a i i ( X ) = 0 for all i and X.

2.6. Collective Potential Energy Function

The collective potential energy function of multiple USVs swarm is a nonnegative function V, which is a mapping: R 3 n R 0 . In this paper, the CPEF is a smooth form of the DEF with the PEF that has a finite cut-off. The most common way to construct a CPEF is to multiply a convex function by using PEF.
First, define a smooth attractive/repulsive potential function ω ( z ) : R 0 R 0 . It has the following properties: when z = d , it has a global minimum, and it has a finite cut-off r.
Next, define the function φ ( X ) is defined as:
φ ( X ) = 1 2 i n j i n ω ( X j X i ) ,
where the φ ( X ) is not differentiable at the singular configuration where two different nodes coincide. In order to solve this problem, we introduce λ -norm into Formula (4), and its algebraic constraints can be written as
X j X i λ = d α , j N i ( X ) ,
where d α = d λ .
The end, according to Formulas (12) and (13), the smooth CPEF of multiple USVs swarm is obtained as
V ( X ) = 1 2 i n j i n ω α ( X j X i λ ) ,
where
ω α ( z ) = d α z ϕ α ( s ) d s , ϕ α ( z ) = ρ h ( z / r α ) ϕ ( z d α ) , ϕ ( z ) = 1 2 [ ( a 1 + a 2 ) λ 1 ( z + a 3 ) + ( a 1 a 2 ) ] ,
and the λ 1 ( z ) = z / 1 + z 2 , and the parameters in ϕ ( z ) are: 0 < a 1 a 2 , a 3 = | a 1 a 2 | / 4 a 1 a 2 , so as to sure ϕ ( 0 ) = 0 .
Through the above analysis, it can be concluded that the CPEF V ( X ) has the following properties: it has a finite cut-off at r α = r λ ; when z = d α , it has a global minimum.

3. Modeling of USV, Virtual Leader, and Grouping Method

3.1. Modeling of USV

In low-speed applications, the damping can be assumed to be linear. Its quadratic velocity term can be ignored. Ship models are usually reduced-order models for control of the horizontal palne motions (surge, sway and yaw). The 3 degrees of freedom (DOFs) kinematic equations of the USV can be expressed in vector form as [31]
η ˙ = J η υ ,
M υ ˙ + D ( υ ) υ = τ + τ ,
where J η : = R ( ψ ) is the rotation matrix, it is given as
R ( ψ ) = cos ( ψ ) sin ( ψ ) 0 sin ( ψ ) cos ( ψ ) 0 0 0 1 ,
with the properties: R ( ψ ) = 1 and R T ( ψ ) R ( ψ ) = I 3 × 3 . η : = [ x , y , ψ ] T is the position and yaw angle in the earth-fixed frame X E O E Y E (see Figure 5). υ : = [ u , v , r ] T is the velocity vector in the body-fixed frame X B O B Y B , where u, v and r represent the surge, sway and yaw angular velocities, respectively. The system inertia matrix M is positive definite and constant, where M = M T R 3 × 3 and M ˙ = 0 . The damping matrix D is also symmetric and positive define, that is D R 3 × 3 . τ : = [ τ 1 , τ 2 , τ 3 ] T is the control input, which is produced by the propellers. τ : = τ w i n d + τ w a v e + τ o c e a n is total environment disturbance by the wind, waves, and ocean currents, respectively.
Through the derivation of the above formula, we can rewritten the Formulas (16) and (17) as
X ˙ = θ , θ ˙ = U ,
where X = η , θ = η ˙ , U = M η 1 ( η ) D η ( υ , η ) η ˙ + M η 1 ( η ) [ τ η + τ η ] .
Assuming that the USVs swarm is composed of n USVs, the motion state equation of the ith USV is described as
X i ˙ = θ i , θ i ˙ = U i ,
where X i = η i R 3 n ( i = 1 , 2 , , n ) is the position and yaw angle vectors, θ i = η ˙ i R 3 n is the velocity vector, U i = M i η i 1 ( η i ) D i η i ( υ i , η i ) η i ˙ + M i η i 1 ( η i ) [ τ i η i + τ i η i ] R 3 n is the control input acting on the ith USV.

3.2. Modeling of Virtual Leader: γ -USV

To get better effect of navigational feedback term in USVs swarm algorithm (see Formula (28) for details), the motion state of virtual leader also conforms to Formula (20). Therefore, the motion state of virtual leader ( γ -USV, γ = 1 , , m ) is as follows [32,33]:
X γ ˙ = θ γ , θ γ ˙ = U γ ,
where X γ , θ γ , and U γ R 3 are the position and yaw angle vector, velocity vector, and control input of the virtual leader, respectively. ( X i ( 0 ) , θ i ( 0 ) ) = ( X d , θ d ) . A static USV has a fixed static that is equal to ( X d , θ d ) for all time. The design of U γ for virtual leader is part of tracking control design for a group of USVs. For example, the choice of U γ 0 leads to that the virtual leader moves along a straight line with a desired velocity θ γ . Based on expression of u i γ (see Formula (29)), a secondary objective of an α - g r i d is to track a dynamic/static virtual leader ( γ -USV).

3.3. Grouping Method

In order to expand the search scope of multiple USVs cluster, it is not enough to increase the number of USVs. By increasing the number of virtual leaders, multi-USV can be grouped and searched from different directions, so as to expand the search scope. Under the circumstances, a multi-USV swarm grouping method based on cosine similarity is proposed.
Suppose that all USVs are in a rectangular region, the coordinates of the lower left corner of the rectangular is Q ( x q , y q ) , and the ith USV initial positions are T i ( x i ( 0 ) , y i ( 0 ) ) , as shown in Figure 6. Connect point Q and T i to form a line segment Q T i . The angle θ i is formed by the line segment Q T i and the positive direction of the x axis. Counter clockwise is positive. Tangent angle θ i is as follows:
θ i = a t a n ( ( y i ( 0 ) y q ) / ( x i ( 0 ) x q ) ) ,
where i = 1 , 2 , , n .
The design process of this method is as follows.
Step 1. According to Formula (22), the tangent angle of the ith USV is θ 1 , θ 2 ,..., θ i ,..., and θ n , respectively. i = 1 , 2 , , n .
Step 2. Sort them by the order of tangent angle. The ordered array satisfies the following condition:
θ k 1 θ k 2 θ k n ,
The index k i of θ k i (i = 1,..., n) represents the ith number in an ordered array from the step 2.
Step 3. USVs were grouped according to the number of virtual leaders. If we divide n by m, the quotient is W and the remainder is c. The value W meets the following condition
W = n / m ,
where · denotes rounding down number.
There are following cases to explain.
1 If c = 0 , all USVs are evenly distributed, and the number of USVs in each group is equal. The subscripts of γ th group USV satisfy the following conditions:
k ( γ 1 ) W + 1 , k ( γ 1 ) W + 2 , , k ( γ 1 ) W + W ,
where γ = 1 , 2 , , m .
2 If c 0 , the subscripts of γ th group USV are as follows:
(1) For the ( m c ) groups, the subscripts of USV are as follows:
k ( γ 1 ) W + 1 , k ( γ 1 ) W + 2 , , k ( γ 1 ) W + W ,
where γ = 1 , 2 , , m c .
(2) For the rest of c groups, the subscripts of USV are as follows:
k ( γ 1 ) ( W + 1 ) + 1 , k ( γ 1 ) ( W + 1 ) + 2 , , k ( γ 1 ) ( W + 1 ) + W + 1 ,
where γ = m c + 1 , , m .

4. USVs Distributed Swarm Control Algorithm

In this section, the distributed swarm control problem for multiple unmanned surface vehicles is addressed by grouping method, distributed swarm control algorithm and other control technology. To accurately estimate environmental disturbances and guarantee estimation error converges in finite time, a finite time disturbance observer is designed. In order to show the control algorithm more intuitively, the schematic diagram is presented in Figure 7.

4.1. USVs Distributed Swarm Control Algorithm

Based on the flocking algorithm proposed in [15], heading angle and yaw angular velocity are considered in the flocking algorithm. In this case, a distributed swarm control algorithm (DSCA) is proposed.
In USVs swarm, a control input u i is used to denote the input signal of the ith USV, which consists of three terms. The input signal u i is as follows:
u i = f i g + f i v + f i γ ,
where f i g , f i v and f i γ are gradient-based term, velocity consensus term and navigational control law, respectively. The specific expressions are as follows:
f i g = p i j N i ϕ α ( X j X i λ ) n i j , f i v = p i j N i a i j ( X ) ( θ j θ i ) , f i γ = u i γ = f i γ ( X i , θ i , X γ , θ γ ) , = p i ( b 1 ( X i X γ ) b 2 ( θ i θ γ ) ) , ( b 1 , b 2 > 0 ) ,
where p i = 1 if ith USV is informed; otherwise, p i = 0 . n i j = λ σ ( X j X i ) = ( ( X j X i ) / 1 + σ X j X i 2 ) is a vector along the line connecting X j to X i . σ ( 0 , 1 ) is a fixed parameter of the λ -norm. the pair ( X γ , θ γ ) is the state of the virtual leader. The parameters b 1 > 0 and b 2 > 0 represents the navigational feedback weight coefficients of X and θ , respectively.
Substituting (29) into (28), the expression (28) is written as
u i = p i ( j N i ϕ α ( X j X i λ ) n i j + j N i a i j ( X ) ( θ j θ i ) b 1 ( X i X γ ) b 2 ( θ i θ γ ) ) , ( b 1 , b 2 > 0 ) .

4.2. USVs Swarm Dynamics

From the Formulas (20) and (30), the collective dynamics equation of the USVs swarm is as follows:
X i ˙ = θ i , θ i ˙ = V ( X i ) L ^ ( X i ) θ i + f i γ ( X i , θ i , X γ , θ γ ) ,
where V ( X ) is a smooth CPEF (see Formula (14)), L ^ ( X ) is the 3-D Laplacian of the proximity net G ( X ) . L ^ = L I m . The L ^ has the following property:
z T L ^ z = 1 2 ( i , j ) E a i j z j z i 2 , z R 3 n ,
where z = c o l ( z 1 , z 2 , , z n ) and z i R 3 n for all i.
When f i γ 0 , the USVs swarm dynamics whose motion form is Formula (31) is a energy consumption system with Hamiltonian (see the Formula (41) for details):
H ξ ( X , θ ) = E ξ ( X ) + 1 2 i = 1 n θ i 2 .
Consider that a translational frame is centered at X c , i.e., the center of mass of all USVs, let A v e ( z ) = ( 1 / n ) i = 1 n z i denote the average of the z i , z = c o l ( z 1 , z 2 , , z n ) . Let X c = A v e ( X ) and θ c = A v e ( θ ) denote the position and yaw angle vector and velocity vector of the origin of the translational frame, respectively. Therefore, X ˙ c ( t ) = θ c ( t ) and θ ˙ c ( t ) = A v e ( u ( t ) ) .
Lemma 1.
([15]) Suppose that the navigational feedback f γ ( X , θ ) : = f i γ ( X , θ , X γ , θ γ ) is linear, i.e., there exists a decomposition of f γ ( X , θ ) in the following form:
f i γ ( X , θ , X γ , θ γ ) = g ( X , θ ) + 1 n h ( X c , θ c , X γ , θ γ ) ,
where g ( X , θ ) = p i b 1 X p i b 2 θ , 1 n = ( 1 , 1 , , 1 ) T , h ( X c , θ c , X γ , θ γ ) = p i b 1 ( X c X γ ) p i b 2 ( θ c θ γ ) , and the pair ( X γ , θ γ ) is the state of the virtual leader.
Substituting (34) to (31), the expression (31) can be consisted of n second-order systems
X ˙ = θ θ ˙ = V ( X ) L ^ ( X ) θ + g ( X , θ )
and one second-order systems
X c ˙ = θ c θ c ˙ = h ( X c , θ c , X γ , θ γ ) .

4.3. Stability Analysis

According to the Lemma 1, we can decompose the stable USVs swarm motion into the combination of the following two forms of stability properties: (1) the stability of certain equilibria of the structural dynamics, and (2) the stability of the expected equilibrium points in translational dynamics. Then, the stability of USVs swarm is equivalent to the stability of these two forms of systems.
According to the Formula (31), the structural dynamics of the multiple USVs cluster has the following forms:
X ˙ = θ , θ ˙ = E ξ ( X ) D d ( X ) θ ,
where
E ξ ( X ) = V ( X ) + ξ J ( X ) ,
J ( X ) = 1 2 i = 1 n X i 2 ,
D d ( X ) = b 2 I 3 + L ^ ( X ) ,
where E λ ( X ) denotes the aggregate potential energy function (APEF), J ( X ) is the moment of inertia of all USVs, ξ = b 1 > 0 is the navigational feedback parameter, D d ( X ) is a positive definite matrix, b 2 > 0 .
Before proving the stability of the system, we introduce a concept: the structural Hamiltonian of the system [15]. Hamiltonian is the sum of the potential energy of all particles plus the kinetic energy of the particles related to the system. In this paper, its expression is as follows:
H ξ ( X , θ ) = E ξ ( X ) + K ( θ ) ,
where K ( θ ) = ( 1 / 2 ) i θ i 2 is the kinetic energy function (KEF). According to Formulas (14), (38) and (39), the expression (41) is written as
H ξ ( X , θ ) = 1 2 i n j i n ω α ( X j X i λ ) + 1 2 i n ξ X i 2 + 1 2 i n θ i 2 = 1 2 i n ( j i n ω α ( X j X i λ ) + ξ X i 2 + θ i 2 ) .
We also need to introduce two other concepts: cohesion of the USV group and cluster.
Definition 1.
Cohesion of the USV group. If there exists a circle with X c ( t ) = A v e ( X ( t ) ) as the center and R > 0 as the radius, which contains all the USV for all time t [ t 0 , t f ] , (i.e., R > 0 : X ( t ) R , t [ t 0 , t f ] ), we called the USV group is cohension.
Definition 2.
Swarms. If the proximity net G ( X ( t ) ) is connected over [ t 0 , t f ) ( t f > t 0 ) , the group of α- g r i d s is called a swarm over the interval [ t 0 , t f ) .
Lemma 2.
The local minima of CPEF V ( X ) is a α- g r i d , and vice versa.
Proof. 
Define a ε -neighborhood of X as
N ε ( X ) = { X R 3 n : X i X i < ε , i V }
where X = c o l ( X 1 , X 2 , , X n ) . If there exists a ε -neighborhood N ε ( X ) of X such that V ( X ) V ( X ) for all X N ε ( X ) , we called the configuration X is a local minima of V ( X ) .
The V ( X ) can be decomposed into two terms: a graph-induced PEF by the proximity net G ( X ) and an integer factor of h 0 = ω α ( r α ) , the expression of V ( X ) is as follows:
V ( X ) = V G ( X ) + k ( X ) h 0 ,
where
V G ( X ) = 1 2 ( i , j ) E ( X ) ω α ( X j X i λ ) ,
k ( X ) = 1 2 ( n ( n 1 ) | E ( X ) | ) I
with I = { 0 , 1 , 2 , , n ( n 1 ) } . For a fixed k I , we have
V ( X ) = V G ( X ) + k h 0 k h 0 .
Any configuration X that achieves the equality in Formula (46) is a local minima of V ( X ) and satisfies V G ( X ) = 0 . ω α ( z ) = 0 if and only if z = d α . Therefore, the configuration X satisfy the following constraints:
X j X i λ = d α X j X i = d , j N i ( X ) .
It means that X is a α - g r i d . The proof of the converse is rather similar and is omitted.
This completes the proof. □
Remark 1.
The environmental disturbances consist of low-frequency part and high-frequency part. In this paper, only the low-frequency part is considered during the control process.Regardless of the rotation and translation of X i , the CPEF V ( X ) remains invariant. Therefore, the V ( X ) has no isolated local minimum.
Theorem 1.
Consider that a group of USVs adopt USVs swarm algorithm of Formula (30) with b 1 , b 2 > 0 and satisfy structural dynamics of Formula (37). Let Ω c = { ( X T , θ T ) R 6 n : = H ξ ( X , θ ) c } be a level-set of the Hamiltonian H ξ ( X , θ ) of Formula (41) such that for any solution starting in Ω c . Suppose that their initial inertia J ( X ( 0 ) ) and the kinetic energy K ( θ ( 0 ) ) are finite. Then, the following statements hold. (i) The group of USVs remain cohesive for all t 0 . (ii) Almost every solution of Formula (37) asymptotically converges to an equilibrium point ( X λ , 0 ) , where X λ is a local minimum of function E ξ ( X ) . (iii) All USVs move asymptotically at the same speed. (iv) It is supposed that the initial structural energy of USVs swarm system is less than ( k + 1 ) c with c = M ω α ( 0 ) and k Z + . Then, at most k distinct pairs of different α-grids may collide. When k = 0 , there is no interaction collision in the USVs swarm system.
Proof. 
According to the symmetry of ω ( z ) = z 2 and adjacency matrix A ( t ) , we can get the following equation:
ω α ( X j i λ ) X j i = ω α ( X j i λ ) X j = ω α ( X j i λ ) X i
with ω α ( X j i λ ) = ω α ( X j X i λ ) = ω α ( X i X j λ ) .
Due to the USVs swarm is an energy consuming system, the time derivative of Formula (41) is obtained as
H ˙ ξ ( X , θ ) = θ T ( b 2 I 3 + L ^ ( X ) ) θ = b 2 ( θ T θ ) θ T L ^ ( X ) θ 0 .
According to the Formula (50), H ξ ( X , θ ) is a monotonically decreasing function for all ( X , θ ) , such that
H ξ ( X ( t ) , θ ( t ) ) H 0 : = H ξ ( X ( 0 ) , θ ( 0 ) ) < + ,
where
H 0 : = H ξ ( X ( 0 ) , θ ( 0 ) ) = V ( X ( 0 ) ) + ξ J ( X ( 0 ) ) + K ( θ ( 0 ) ) .
According to the Formulas (51) and (52) and the assumption that initial inertia J ( X ( 0 ) ) and the kinetic energy K ( θ ( 0 ) ) are finite, we have
E ξ ( X ( t ) ) H 0 , t 0 E ξ ( X ( t ) ) = V ( X ( t ) ) + ξ 2 X T ( t ) X ( t ) . ξ > 0
From the Formula (53), the following equation is obtained:
V ( X ( t ) ) + ξ 2 X T ( t ) X ( t ) H 0 X T ( t ) X ( t ) 2 ( H 0 V ( X ( t ) ) ) ξ X T ( t ) X ( t ) 2 H 0 ξ . t 0
According to the same proof process, yields
K ( θ ( t ) ) = 1 2 i n θ i ( t ) 2 H 0 , t 0 θ T ( t ) θ ( t ) 2 H 0 . t 0
According to the Formulas (50)–(55), we can draw the following conclusion: the USVs swarm can achieve aggregation for all t 0 , and keep in a circle neighborhood with X c as the center and R = 2 H 0 / ξ as the radius.
Therefore, conclusion (i) is valid.
From to the Formula (50), H ξ ( X , θ ) is a monotonically decreasing function for all ( X , θ ) . Hence, H ξ ( X , θ ) c for all t 0 that implies the Ω c = { ( X T , θ T ) R 6 n : = H ξ ( X , θ ) c } is an invariant set.
From the LaSalle’s invariance principle [34], all the solutions of Formula (41) starting in Ω c , will converge to the largest invariant set in S. The expression of S is as follows:
S = { ( X T , θ T ) R 6 n : = H ˙ ξ ( X , θ ) = 0 } .
Thus, all the solutions are in the set Ω c . Any solution ( X ( t ) , θ ( t ) ) of the USVs swarm algorithm is uniquely mapped to a solution ( X ( t ) , θ ( t ) ) of the multiple USVs system.
According to the Formula (55), the velocity will be bounded as
K ( θ ( t ) ) = 1 2 i n θ i ( t ) 2 H 0 c . t 0 θ T ( t ) θ ( t ) 2 H 0 2 c . t 0
This ensures that the upper bound of uniform convergence of velocity is θ ( t ) 2 c .
The Formulas (53)–(57) guarantee that the solution of the USVs swarm is bounded. Let z = c o l ( X , θ ) , we have
z ( t ) 2 = X T ( t ) X ( t ) + θ T ( t ) θ ( t ) 2 ( 1 ξ + 1 ) H 0 R 2 + 2 c : = C < + ,
where C > 0 is a constant.
From the Formula (50), we draw the following conclusions: if θ = 0 , H ˙ ξ ( X , θ ) = 0 ; otherwise, H ˙ ξ ( X , θ ) < 0 . Therefore, the USVs swarm are always in motion for all t 0 . Thus, the velocities of the group of USVs match in the moving frame, or θ 1 = θ 2 = = θ n . When i θ i = 0 , we can get θ i = 0 ( i = 1 , 2 , , n ) . In the reference frame, the velocities of USVs are asymptotically consistent. Meanwhile, X 1 = X 2 = = X n . Thus, E ξ ( X ξ ) = 0 . Therefore, almost every solution of the multiple USVs system asymptotically converges to the equilibrium point z ξ = ( X ξ , 0 ) , where X ξ is a local minima of E ξ ( X ) . According to the Lemma 2, the local minima of CPEF V ( X ) is a α - g r i d . These prove that (ii) is holds.
We now prove (iii) of Theorem 1. Through the above proof of the (i) and (ii), we can draw the following conclusion: the velocities of all USVs are asymptotically consistent. Thus, the velocity of each USV asymptotically consistent in the reference frame.
In order to prove (iv), we use the counter evidence method to prove.
According to the information in the (iv), we can conclude that the following inequality holds:
H 0 < ( k + 1 ) c = ( k + 1 ) ω α ( 0 ) , k Z + .
Thus, there are at most k distinct pairs of USVs that collide at a given time t 1 0 .
Next, we make the following assumption: there must be at least k + 1 distinct pairs of USVs collide with each other at a given time t 1 . Therefore, we get the following conclusion: the Hamiltonian H 0 of the USVs swarm is at least ( k + 1 ) ω α ( 0 ) at a given time t = t 1 . Then, the following inequality holds:
H 0 ( k + 1 ) ω α ( 0 ) ( k + 1 ) c = ( k + 1 ) ω α ( 0 ) H 0 ( k + 1 ) c .
Obviously, the Formulas (59) and (60) are contradictory. Therefore, our hypothesis is not tenable. Therefore, no more than k distinct pairs of the USVs can possibly collide at any time t > 0 .
Next, we prove that: there is no collision between USV when k = 0 . When k = 0 , the Hamiltonian H 0 < c = ω α ( 0 ) . Assuming that there exists two USV colliding at t = t 1 , the two USVs are labeled as n 1 th USV and n 2 th USV, respectively. i.e., X n 1 ( t 1 ) = X n 2 ( t 1 ) . From the Formula (44), we can get the following results:
V ( X ( t ) ) = 1 2 i n i j n ω α ( X j X i λ ) = 1 2 i V \ { n 1 , n 2 } n j V \ { i , n 1 , n 2 } n ω α ( X j X i λ ) + ω α ( X n 1 ( t ) X n 2 ( t ) λ ) ω α ( X n 1 ( t ) X n 2 ( t ) λ ) .
Substituting t = t 1 into Formula (61), the following inequality holds:
V ( X ( t 1 ) ) ω α ( 0 ) = c .
However, from the Formula (33), the inequality is shown as follows:
V ( X ( t ) ) H ( X ( t ) , θ ( t ) ) K ( θ ( t ) ) V ( X ( t ) ) H ( X ( t ) , θ ( t ) ) c < c .
Obviously, the Formulas (62) and (63) are contradictory. Therefore, our hypothesis is untenable. Thus, no two USVs collide with k = 0 at any time t 0 .
This completes the proof. □

5. Simulation Results

In simulations, the model of surface ship Cybership II is used in [35]. The time-varying environmental disturbances are modeled as the white Gaussian noise. Numerical simulation and comparison results are given to show the effectiveness of the proposed DSCA. We consider twenty follower vehicles ( n = 20 ) under different number virtual leader ( γ USV), i.e., case 1, case 2 and comparison group, respectively.
The parameters of DSCA are choose as follows. The sensing radius is choose at r = 5 m , the desired distance d = 3 m , λ = 0.1 for the λ -norm, h = 0.2 for the bump function ρ h ( z ) . After a lot of simulation testing and verification, the navigation feedback gain are set as b 1 = 0.01 and b 2 = 0.15 , respectively.

5.1. Case 1: 20 USVs and a Virtual Leader

In this subsection, the simulation results are given to verify the performance of proposed DSCA of the USVs with a virtual leader.
The initial positions of 20 USVs and a virtual leader are randomly distribute in the range of abscissa and ordinate [−60 m, 60 m] × [−80 m, 60 m], and heading angle of the USVs is also randomly distributed in [ π r a d , π r a d ]. The initial parameters of the virtual leader are X γ ( 0 ) = [ x ( 0 ) m , y ( 0 ) m , ψ ( 0 ) r a d ] T , where ( x ( 0 ) , y ( 0 ) ) denotes the position of the initial moment, and ψ ( 0 ) = π / 4 , and θ γ ( 0 ) = [ 2 m / s , 2 m / s , 0 m / s 2 ] T .
The simulation results are shown in Figure 8 and Figure 9. Figure 8 shows the topology structure of 20 USVs and a virtual leader at different times. As shown in Figure 8a, it represnets the position of 20 USVs and a virtual leader at the initial time. Figure 8b shows that 19 USVs are clustered at 30 s . Figure 8c represents that 20 USVs are clustered at 40 s . Figure 8d represents that 20 USVs can always maintain cluster navigation. During sail, the topology structure composed of multiple USVs remained stable. Figure 9a represents the trajectories of 20 USVs and a virtual leader. Figure 9b shows the positions and headings (X). Figure 9c represents the tracking errors. Figure 9d shows the USVs speeds and velocities mismatch. The simulation results show that the group of USVs remain cohesive and all USVs move asymptotically at the same speed.
Since the initial conditions (x, y positions, heading and velocity) of the USVs are different, there are small deviations from vessel to vessel. Apparently, the topology show that no collisions occur during the whole swarming process. Through the comparison of Figure 8b–d, it can be seen that the topology structure is time-varying. Quasi α - g r i d is unstable. After a long enough time, quasi α - g r i d will eventually form α - g r i d . It shows the effective of the proposed time-variant topology structure. Through the simulation results, it can be found that under the proposed DSCA, the positions converge uniformly without collision, the heading errors tend to zero, and velocities consensus are also achieved.

5.2. Case 2: 20 USVs and Three Virtual Leaders

In this subsection, the simulation results are given to verify the performance of proposed DCSA with multi-USV under three virtual leaders.
The initial positions of the 20 USVs are distribute in the range of abscissa and ordinate [−100 m, 50 m] × [−100 m, 50 m], and heading angle of the USVs is also randomly distributed in [ π r a d , π r a d ].
Different parameters are set for the initial state of the three virtual leaders. The parameters of the first virtual leader are X γ 1 ( 0 ) = [ 20 m , 80 m , π / 12 r a d ] T and θ γ 1 ( 0 ) = [ 1.5 m / s , 1.5 m / s , 0 m / s 2 ] T . The parameters of the second virtual leader are X γ 2 ( 0 ) = [ 20 m , 20 m , π / 4 r a d ] T and θ γ 2 ( 0 ) = [ 2 m / s , 2 m / s , 0 m / s 2 ] T . The parameters of the third virtual leader are X γ 3 ( 0 ) = [ 80 m , 10 m , 5 π / 12 r a d ] T and θ γ 3 ( 0 ) = [ 2.5 m / s , 2.5 m / s , 0 m / s 2 ] T .
Figure 10a represents the trajectories of 20 USVs and three virtual leaders. Figure 10b shows the positions and headings (X). Figure 10c represents the tracking errors. Figure 10d shows the USVs speeds and velocities mismatch. The simulation results show that the grouping method based on cosine similarity is effectiveness, the group of USVs remain cohesive and all USVs move asymptotically at the same speed.

5.3. Comparison Group

The comparison group is carried out by using a distributed convex optimization proposed in [36].
The controller is designed by using only local interaction information rather than global bidirectional topology, and the case of multi-USV grouping is not considered in [36]. In simulations, the control inputs and environment distribution of two methods are both same. The simulation results are shown in Figure 11 and Figure 12.
Figure 11 shows the topology structure of 20 USVs and a virtual leader at different times. As shown in Figure 11a, it represnets the position of 20 USVs and a virtual leader at the initial time. Figure 11b shows that 12 USVs are clustered at 30 s. Figure 11c represents that 17 USVs are clustered at 50 s. Figure 11d represents that 20 USVs are clustered at 100 s. Figure 12a represents the trajectories of 20 USVs and a virtual leader. Figure 12b shows the positions and headings (X). Figure 12c represents the tracking errors. Figure 12d shows the USVs speeds and velocities mismatch.
After 50 experiments, the number of USVs achieving clusters by using defferent algorithms are shown in Table 2. Through the comparative analysis of the data in Table 2, it is shown that the number of USV obtained by method distributed swarm control algorithm is more than that obtained by distributed convex optimization. It can be proved that the effectiveness of the proposed algorithm.
In order to study the influence of environmental disturbances on the algorithm, we use a finite time disturbance observer to estimate the environmental disturbances. The comparative simulation is carried out by using a modular design approach proposed in [37].
The simulation result is shown in Figure 13. Taking the first USV as an example, Figure 13 depicts that the time-varying environmental disturbances value and estimates obtained by using modular design approach and finite time disturbance observer.
Through the analysis of Figure 13, it is shown that the more accurate extimation of the environmantal disturbances is obtained by using finite time disturbance observer than modular design approach. This means that the proposed control scheme is robust against disturbances.

6. Conclusions

In this paper, we investigate the dynamic properties of the group for the case where the number of virtual leader is different and the topology of the neighborhood relations between vehicles is dynamic. Firstly, to achieve multiple USVs connectivity, a new time-variant topology structure is proposed, which is based on proximity nets and α - g r i d s . Secondly, to reduce the communication load and realize the information diffluence, a grouping method based on cosine similarity is proposed. Thirdly, to ensure the high effeciency of information transmission and the reduction of costs, a distributed swarm control algorithm is proposed. To analyze the stability of system, the concept of translation framework is introduced. The stability of multiple USVs system is proved by using the Hamiltonian and LaSalle invariance principle. Simulation results show the effectiveness of the proposed methods. However, this paper does not consider the problem of static/dynamic obstacles. These problems will be studied in the next step.

Author Contributions

Conceptualization, G.X. and X.S.; methodology, X.S.; software, X.S. and X.X.; formal analysis, X.S. and X.X.; investigation, X.S.; supervision, G.X. and X.X.; funding acquisition, G.X. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the 7th Generation Ultra Deep Water Drilling Unit Innovation Project.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
USVunmanned surface vehicle
DSCAdistributed swarm control algorithm
DEFdeviation energy function
PEFpotential energy function
CPEFcollective potential energy function
APEFaggregate potential energy function
KEFkinetic energy function

References

  1. Singh, Y. Cooperative Swarm Optimisation of Unmanned Surface Vehicles. Ph.D. Dissertation, University of Plymouth, Plymouth, UK, 2019. [Google Scholar]
  2. Mina, T.; Singh, Y.; Min, B.C. A Novel Double Layered Weighted Potential Field Framework for Multi-USV Navigation towards Dynamic Obstacle Avoidance in a Constrained Maritime Environment. In Proceedings of the OCEANS 2019 MTS/IEEE SEATTLE, Seattle, WA, USA, 27–31 October 2019; pp. 1–9. [Google Scholar]
  3. Mina, T.; Singh, Y.; Min, B.-C. Maneuvering Ability-Based Weighted Potential Field Framework for Multi-USV Navigation, Guidance, and Control. Mar. Technol. Soc. J. 2020, 19, 40–58. [Google Scholar] [CrossRef]
  4. Huang, Y.; Ding, H.; Zhang, Y.; Wang, H.; Cao, D.; Xu, N.; Hu, C. A Motion Planning and Tracking Framework for Autonomous Vehicles Based on Artificial Potential Field Elaborated Resistance Network Approach. IEEE Trans. Ind. Electron. 2020, 67, 1376–1386. [Google Scholar] [CrossRef]
  5. Shi, P.; Yan, B. A Survey on Intelligent Control for Multiagent Systems. Systems 2005, 51, 161–175. [Google Scholar] [CrossRef]
  6. Liang, X.; Qu, X.; Wang, N.; Li, Y.; Zhang, R. Swarm control with collision avoidance for multiple underactuated surface vehicles. Ocean Eng. 2019, 41, 106516. [Google Scholar] [CrossRef]
  7. Ning, B.; Han, Q.L.; Zuo, Z.; Jin, J.; Zheng, J. Collective Behaviors of Mobile Robots Beyond the Nearest Neighbor Rules With Switching Topology. IEEE Trans. Cybern 2018, 418, 1577–1590. [Google Scholar] [CrossRef]
  8. Jia, Y.; Wang, L. Leader–Follower Flocking of Multiple Robotic Fish. IEEE/ASME Trans. Mechatronics 2015, 20, 1372–1383. [Google Scholar] [CrossRef]
  9. Messina, F.; Vasilakos, A.; Meo, P.D. Introduction to the special section on Recent trends in flocking control and communication for Unmanned vehicles. Comput. Electr. Eng. 2019, 80, 1–10. [Google Scholar] [CrossRef]
  10. Jing, G.; Wang, L. Multiagent Flocking With Angle-Based Formation Shape Control. IEEE Trans. Autom. Control 2020, 65, 817–823. [Google Scholar] [CrossRef]
  11. Reyes, L.A.V.; Tanner, H.G. Flocking, Formation Control, and Path Following for a Group of Mobile Robots. IEEE Trans. Control. Syst. Technol. 2015, 23, 1268–1282. [Google Scholar] [CrossRef]
  12. Xia, G.; Sun, C.; Zhao, B.; Xue, J. Cooperative Control of Multiple Dynamic Positioning Vessels with Input Saturation Based on Finite-time Disturbance Observer. Int. J. Control 2019, 17, 370–379. [Google Scholar] [CrossRef]
  13. Jia, Y.; Li, Q.; Qiu, S. Distributed Leader-Follower Flight Control for Large-Scale Clusters of Small Unmanned Aerial Vehicles. IEEE Access 2018, 6, 32790–32799. [Google Scholar] [CrossRef]
  14. Ibuki, T.; Wilson, S.; Yamauchi, J.; Fujita, M.; Egerstedt, M. Optimization-Based Distributed Flocking Control for Multiple Rigid Bodies. IEEE Robot. Autom. Lett. 2020, 5, 1891–1898. [Google Scholar] [CrossRef]
  15. Olfati-Saber, R. Flocking for multi-agent dynamic systems: Algorithms and theory. IEEE Trans. Autom. Control 2006, 51, 401–420. [Google Scholar] [CrossRef] [Green Version]
  16. Wei, H.; Chen, X.-B. Flocking for Multiple Subgroups of Multi-Agents With Different Social Distancing. IEEE Access 2020, 8, 164705–164716. [Google Scholar] [CrossRef]
  17. Chen, S.; Pei, H.; Lai, Q.; Yan, H. Multitarget Tracking Control for Coupled Heterogeneous Inertial Agents Systems Based on Flocking Behavior. IEEE Trans. Syst. Man Cybern. Syst. 2019, 49, 2605–2611. [Google Scholar] [CrossRef]
  18. Jing, G.; Zheng, Y.; Wang, L. Flocking of multi-agent systems with multiple groups. Int. J. Control 2014, 87, 2573–2582. [Google Scholar] [CrossRef]
  19. Olfati-Saber, R.; Murray, R.M. Consensus problems in networks of agents with switching topology and time-delays. IEEE Trans. Autom. Control 2004, 49, 1520–1533. [Google Scholar] [CrossRef] [Green Version]
  20. Ren, W.; Beard, R.W. Consensus seeking in multi-agent systems under dynamically changing interaction topologies. IEEE Trans. Autom. Control 2005, 50, 655–661. [Google Scholar] [CrossRef]
  21. Gazi, V.; Passino, K.M. Stability analysis of swarms. IEEE Trans. Autom. Control 2003, 48, 692–697. [Google Scholar] [CrossRef] [Green Version]
  22. Moreau, L. Stability of multiagent systems with time-dependent communication links. IEEE Trans. Autom. Control 2007, 50, 169–182. [Google Scholar] [CrossRef]
  23. Dai, J.; Zhu, S.; Chen, C.; Guan, X. Connectivity-Preserving consensus algorithms for Multi-agent Systems. IFAC Proc. Vol. 2011, 44, 5675–5680. [Google Scholar] [CrossRef]
  24. Fax, J.A.; Murray, R.M. Information flow and cooperative control of vehicle formations. IEEE Trans. Robot. Autom. 2004, 49, 1465–1476. [Google Scholar] [CrossRef] [Green Version]
  25. Tanner, H.G.; Jadbabaie, A.; Pappas, G.J. Flocking in fixed and switching networks. IEEE Trans. Autom. Control 2007, 52, 863–868. [Google Scholar] [CrossRef]
  26. Chen, Y.; Wang, Z.; Liang, J. Automatioc dynamic flocking in mobile actuator sensor networks by central Voronoi tessellations. In Proceedings of the IEEE International Conference Mechatronics and Automation, Niagara Falls, ON, Canada, 29 July–1 August 2005; pp. 1630–1635. [Google Scholar]
  27. Semnani, S.H.; Basir, O.A. Semi-flocking algorithm for motion control of mobile sensors in large-scale surveillance systems. IEEE Trans. Cybern 2015, 45, 129–137. [Google Scholar] [CrossRef] [PubMed]
  28. Pei, H.; Chen, S.; Lai, Q. A local flocking algorithm of multi-agent dynamic systems. Int. J. Control 2015, 88, 2242–2249. [Google Scholar] [CrossRef]
  29. Zhang, H.-T.; Zhai, C.; Chen, Z. A General Alignment Repulsion Algorithm for Flocking of Multi-Agent Systems. Trans. Autom. Control 2011, 56, 430–435. [Google Scholar] [CrossRef]
  30. Zhu, J.; Lu, J.; Yu, X. Flocking of Multi-Agent Non-Holonomic Systems With Proximity Graphs. IEEE Trans. Circuits Syst. I Regul. Pap. 2013, 60, 199–210. [Google Scholar] [CrossRef]
  31. Fossen, T.I.; Strand, J.P. Passive nonlinear observer design for ships using Lyapunov methods: Full-Scale experiments with a supply vessel. Automatica 1999, 35, 3–16. [Google Scholar] [CrossRef]
  32. Su, H.; Wang, X.; Lin, Z. Flocking of Multi-Agents With a Virtual Leader. IEEE Trans. Autom. Control 2009, 54, 293–307. [Google Scholar] [CrossRef]
  33. Shi, H.; Wang, L.; Chu, T. Flocking of multi-agent systems with a dynamic virtual leader. Int. J. Control 2008, 82, 43–58. [Google Scholar] [CrossRef]
  34. Sun, Y. Stability analysis of flocking for multi-agent dynamic systems. Nonlinear Anal. Real World Appl. 2013, 14, 1075–1081. [Google Scholar] [CrossRef]
  35. Skjetne, R.; Fossen, T.I.; Kokotović, P.V. Adaptive maneuvering, with experiments, for a model ship in a marine control laboratory. Automatica 2005, 41, 289–298. [Google Scholar] [CrossRef]
  36. Zhang, Q.; Gong, Z.; Yang, Z.; Chen, Z. Distributed Convex Optimization for Flocking of Nonlinear Multi-agent Systems. Int. J. Control. Autom. Syst. 2019, 17, 1177–1183. [Google Scholar] [CrossRef]
  37. Peng, Z.; Wang, D.; Wang, J. Cooperative Dynamic Positioning of Multiple Marine Offshore Vessels: A Modular Design. IEEE/ASME Trans. Mechatronics 2016, 21, 1210–1221. [Google Scholar] [CrossRef]
Figure 1. The structure of the proposed grouping method and distributed swarm control algorithm for multiple USVs system.
Figure 1. The structure of the proposed grouping method and distributed swarm control algorithm for multiple USVs system.
Jmse 09 01324 g001
Figure 2. The ith USV and its circle neighborhood (Ni represents the set of neighbors of ith USV, i = 1, 2, …, 6).
Figure 2. The ith USV and its circle neighborhood (Ni represents the set of neighbors of ith USV, i = 1, 2, …, 6).
Jmse 09 01324 g002
Figure 3. Structure of α - g r i d s : regular triangle and square.
Figure 3. Structure of α - g r i d s : regular triangle and square.
Jmse 09 01324 g003
Figure 4. Quasi α - g r i d : multiple regular polygons.
Figure 4. Quasi α - g r i d : multiple regular polygons.
Jmse 09 01324 g004
Figure 5. Earth-fixed frame and body-fixed frame.
Figure 5. Earth-fixed frame and body-fixed frame.
Jmse 09 01324 g005
Figure 6. USVs are distributed in the rectangular region.
Figure 6. USVs are distributed in the rectangular region.
Jmse 09 01324 g006
Figure 7. The diagram of the proposed grouping method, distributed swarm control algorithm and other control technology for multiple USVs.
Figure 7. The diagram of the proposed grouping method, distributed swarm control algorithm and other control technology for multiple USVs.
Jmse 09 01324 g007
Figure 8. Motion instants of swarming/automated rendezvous using DSCA of Case 1. The blue ship shape and “+” represent the vessel and virtual leader, respectively. (a) t = 0.1 s. (b) t = 30 s. (c) t = 40 s. (d) t = 100 s.
Figure 8. Motion instants of swarming/automated rendezvous using DSCA of Case 1. The blue ship shape and “+” represent the vessel and virtual leader, respectively. (a) t = 0.1 s. (b) t = 30 s. (c) t = 40 s. (d) t = 100 s.
Jmse 09 01324 g008
Figure 9. The simulation results using DSCA of Case 1. (a) Trajectories. (b) Positions and headings. (c) Tracking errors. (d) USV speeds and velocity mismatch.
Figure 9. The simulation results using DSCA of Case 1. (a) Trajectories. (b) Positions and headings. (c) Tracking errors. (d) USV speeds and velocity mismatch.
Jmse 09 01324 g009
Figure 10. The simulation results using DSCA of Case 2. (a) Trajectories. (b) Positions and headings. (c) Tracking errors. (d) USV speeds and velocity mismatch.
Figure 10. The simulation results using DSCA of Case 2. (a) Trajectories. (b) Positions and headings. (c) Tracking errors. (d) USV speeds and velocity mismatch.
Jmse 09 01324 g010
Figure 11. Motion instants of swarming/automated rendezvous by using proposed method in [36]. The blue ship shape and “+” represent the vessel and virtual leader, respectively. (a) t = 0.1 s. (b) t = 30 s. (c) t = 50 s. (d) t = 100 s.
Figure 11. Motion instants of swarming/automated rendezvous by using proposed method in [36]. The blue ship shape and “+” represent the vessel and virtual leader, respectively. (a) t = 0.1 s. (b) t = 30 s. (c) t = 50 s. (d) t = 100 s.
Jmse 09 01324 g011
Figure 12. The simulation results by using proposed method in [36]. (a) Trajectories. (b) Positions and headings. (c) Tracking errors. (d) USV speeds and velocity mismatch.
Figure 12. The simulation results by using proposed method in [36]. (a) Trajectories. (b) Positions and headings. (c) Tracking errors. (d) USV speeds and velocity mismatch.
Jmse 09 01324 g012
Figure 13. The time-varying environmental disturbances value and estimates obtained by using modular design approach and finite time disturbance observer.
Figure 13. The time-varying environmental disturbances value and estimates obtained by using modular design approach and finite time disturbance observer.
Jmse 09 01324 g013
Table 1. Notations and variables used in this paper.
Table 1. Notations and variables used in this paper.
VariableDefinition
A \ B Set of elements belonging to A but not belonging to B
· Absolute value of a scalar
· Euclidean norm
R n × n n × n dimensional Euclidean Space
d i a g { a i } A block-diagonal matrix with a i being the ith diagonal element
( · ) T Transpose of a matrix
( · ) 1 Inverse of a matrix
Kronecker product of matrix
A R n × n Adjacency matrix defined as A = [ a i j ] n × n with a i j = a j i
d i Defined as d i = j = 1 n a i j
D R n × n Degree matrix defined as D = d i a g { d 1 , d 2 , , d n }
L R n × n Laplacian matrix defined as L = D A
Table 2. The number of USVs achieving clusters via different algorithms ( n = 20 ).
Table 2. The number of USVs achieving clusters via different algorithms ( n = 20 ).
t = 0 st = 30 st = 50 st = 100 s
distributed swarm control algorithm0192020
distributed convex optimization0121720
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xia, G.; Sun, X.; Xia, X. Distributed Swarm Control Algorithm of Multiple Unmanned Surface Vehicles Based on Grouping Method. J. Mar. Sci. Eng. 2021, 9, 1324. https://doi.org/10.3390/jmse9121324

AMA Style

Xia G, Sun X, Xia X. Distributed Swarm Control Algorithm of Multiple Unmanned Surface Vehicles Based on Grouping Method. Journal of Marine Science and Engineering. 2021; 9(12):1324. https://doi.org/10.3390/jmse9121324

Chicago/Turabian Style

Xia, Guoqing, Xianxin Sun, and Xiaoming Xia. 2021. "Distributed Swarm Control Algorithm of Multiple Unmanned Surface Vehicles Based on Grouping Method" Journal of Marine Science and Engineering 9, no. 12: 1324. https://doi.org/10.3390/jmse9121324

APA Style

Xia, G., Sun, X., & Xia, X. (2021). Distributed Swarm Control Algorithm of Multiple Unmanned Surface Vehicles Based on Grouping Method. Journal of Marine Science and Engineering, 9(12), 1324. https://doi.org/10.3390/jmse9121324

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