Next Article in Journal
Experimental and Numerical Study on Residual Strength of Honeycomb Sandwich Composite Structure after Lightning Strike
Next Article in Special Issue
Constrained Motion Planning of 7-DOF Space Manipulator via Deep Reinforcement Learning Combined with Artificial Potential Field
Previous Article in Journal
Separation Reliability Analysis for the Low-Shock Separation Nut with Mechanism Motion Failure Mode
Previous Article in Special Issue
Bifurcations of Periodic Orbits in the Gravitational Field of Irregular Bodies: Applications to Bennu and Steins
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Finite-Time Distributed Leaderless Cooperative Guidance Law for Maneuvering Targets under Directed Topology without Numerical Singularities

1
School of Automation Science and Electrical Engineering, Beihang University, Beijing 100191, China
2
Science and Technology on Complex System Control and Intelligent Agent Cooperation Laboratory, Beijing 100074, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(3), 157; https://doi.org/10.3390/aerospace9030157
Submission received: 16 January 2022 / Revised: 25 February 2022 / Accepted: 2 March 2022 / Published: 13 March 2022
(This article belongs to the Special Issue Recent Advances in Spacecraft Dynamics and Control)

Abstract

:
A new distributed leaderless cooperative guidance algorithm is suggested for multi-directional saturation strikes against maneuvering targets. First, the finite-time disturbance observer (FDO) is used to estimate the target’s unknown maneuvers. After that, the guidance laws along the line-of-sight (LOS) direction and the LOS normal direction are designed separately based on the finite-time consensus theory. The proposed guidance law could satisfy both impact time and LOS angle constraints. Numerical singularities due to feedback linearization are avoided by nonlinear stability analysis. Finally, the effectiveness of the proposed cooperative guidance law is proved by numerical simulation arithmetic.

1. Introduction

In the last decade, cooperative guidance has received significant attention for its effectiveness in increasing the rate of defense penetration and aircraft-target interception. Well-coordinated multiple vehicles tend to be more efficient. They can raise the effectiveness of defense penetration by attacking the target from numerous different directions simultaneously. Cooperative guidance against maneuvering targets allows for a multi-vehicle-to-target interception posture. As a result, the target’s cost of maneuvering to escape increases, and the interception success rate improves.
The cooperative guidance law can be functionally divided into impact angle cooperation and impact time cooperation. The information exchanged between the vehicles was not considered in the initial study. Instead, the cooperative attack was achieved through predefined impact angles and impact times [1,2,3,4]. The limitation of this method is that it is difficult to find reasonable predetermined values. If the time and angle of impact are not set properly, there is a risk of increasing energy consumption and even resulting in guidance failure.
The current research focuses on cooperative guidance. The remaining flight time and impact angle are selected as coordination variables. Then the distributed cooperative guidance law (DCGL) was established by local communication based on consensus theory. Among the existing literature on cooperative guidance, the most extensive research has been conducted on cooperative time-coordinated guidance against a stationary target [5,6,7,8,9,10,11]. Wang, et al. [5] proposed a two-step guidance algorithm. At the first step, a DCGL was developed based on the consensus theory of the second-order system. During the second step, a simultaneous attack was achieved using a proportional guidance method. Zhao, et al. [6] proposed a new cooperative guidance method that reduces the computational effort by triggering only at a specific time. Jiang, et al. [7] simultaneously attacking with multiple constraints was realized based on the backward horizon control (RHC) algorithm. Zhang, et al. [8] investigated the issue of optimal DCGL for stationary targets under directed topologies. A two-stage guidance method was developed to optimize energy consumption while ensuring simultaneous collision against the target. Chen, et al. [9] proposed a cooperative guidance law for a vehicle with thrust control. This guidance law enables a coordinated attack under a hit angle constraint considering the velocity constraint. Li, et al. [10] investigated the issue of simultaneous arrival of multiple interceptors with effective partial actuators. A fault-tolerant cooperative guidance method was proposed, where an adaptive method was devised to handle uncertainties. Simultaneous arrivals within a fixed time interval under actuator failure conditions were achieved. Teng, et al. [11] suggested a new cooperative guidance method that achieves simultaneous hits in multiple directions without radial velocity measurement.
Currently, there are limited results of cooperative guidance for maneuvering targets [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26]. Nikusokhan, et al. [12], it was hypothesized that the linearization condition of the engagement dynamics could be satisfied. The measurement information of the target acceleration, which is difficult to obtain, is directly used in [12,13,14,15,16,17,18]. Wang, et al. [19] presented a three-dimensional DCGL for several vehicles to strike a maneuvering target at predetermined impact angles. Yu, et al. [20] studied the design and analysis of DCGL against maneuvering targets. An extended state observer is first utilized to evaluate the target’s maneuver. On this basis, a cooperative guidance law that enables a head-on saturation attack is proposed based on the leader-follower model. Wang, et al. [21] proposed a DCGL for hypersonic vehicles that solves the simultaneous arrival problem in the presence of uncontrollable velocity. Chen, et al. [22] proposed a three-dimensional nonlinear DCGL that enables multiple vehicles to simultaneously attack a maneuvering target at a predetermined LOS angle. Liu, et al. [23] investigated robust differential games and their application in cooperative guidance. The suggested guidance method is able to avoid input saturation while synchronizing arrival times.
At present, the following problems still exist with cooperative guidance for maneuvering targets. The first issue is how to estimate the maneuver of the target. In [12,13,14,15,16,17] need to obtain the acceleration of the target directly, which is usually difficult to measure directly by sensors. In the literature [20,21], it also needs to be assumed that the target acceleration is constant or slowly varying. The second problem is the large tangential acceleration command for the vehicle. In the literature [24,25], the acceleration signal is made singular at the end of the guidance due to the feedback linearization method. In practice, the tangential acceleration of the aircraft cannot respond to excessive commands. In the literature [19,22,26], both the missile-target distance and the rate of the distance are needed to achieve consensus. This is unnecessary and could lead to larger energy consumption. Finally, in terms of communication topology, the approach in the literature [18] is centralized, while the approach in the literature [19,22,24,25] can only be applied to undirected topologies.
This paper investigates the cooperative guidance issue for maneuvering targets and proposes a new guidance law with the following advantages:
  • The suggested guidance strategy could be used for the maneuvering target. At the same time, the guidance method does not need direct access to the acceleration measurement information of the target and does not need the acceleration of the target to be constant or slowly varying.
  • The design and analysis of the guidance law are conducted directly on the nonlinear model, avoiding the disadvantages of numerical singularities and excessive guidance commands due to feedback linearization. It is also shown that the guidance law is finite-time converged based on the homogeneous system stability theory.
  • The proposed guidance method is distributed and only requires neighborhood information. At the same time, he suggested that the method can be used to a directed communication topology through a rigorous stability analysis.
The rest of the paper is structured as follows. In Section 2, some necessary preliminaries are presented. In Section 3, the studied cooperative guidance problem is described. In Section 4, the major results, including the design of the distributed guidance algorithm and the stability analysis, are elaborated. In Section 5, the simulation verification on the proposed distributed guidance algorithm is performed. Lastly in Section 6, conclusions are drawn.

2. Preliminaries

Lemma 1.
Under Assumption 1, there exists a positive column vector θ   =   θ 1 , , θ N T such that θ i a i j = θ j a j i for all i , j = 1 , N . Denote Θ = diag θ 1 , , θ N . Then, the matrix Θ L is symmetric.
Lemma 2.
[27] Consider the system
x ˙ ( t )   =   f ( t , x ( t ) ) , x ( 0 )   =   x 0 .
Assume that there is a continuous positive definite function V ( x ) : D for any real number c > 0 and α ( 0 , 1 ) , the following inequality holds
V ˙ ( x )   +   c V α ( x ) 0 , x N \ { 0 } .
Then, the origin is a finite-time stable equilibrium of (1), and the settling time is
T x 0     1 c ( 1 α ) V ( 1 α ) x 0 .
Lemma 3.
[28] Let ξ 1 , ξ 2 , , ξ N 0 and 0 < p 1 . Then
i = 1 N ξ i p     i = 1 N ξ i p .
Definition 1.
Select the following system
x ˙ = f ( x ) , f ( 0 ) = 0 , x ( 0 ) = x 0 , x m
where f ( x )   =   f 1 ( x ) , f 2 ( x ) , , f m ( x ) T : U 0 m is a continuous vector field in an open neighborhood U 0 around the origin. Let r 1 , r 2 , , r m m with r i > 0 ,   i = 1 , 2 , , m . f ( x ) is called homogeneous of degree κ with respect to r 1 , r 2 , , r m if for any given ε > 0 , f i ε r 1 x 1 , ε r 2 x 2 , , ε r m x m   =   ε κ + r i f i ( x ) ,   i = 1 , 2 , , m ,   x m . System (5) is called homogeneous if f ( x ) is homogeneous.
Lemma 4.
Consider system (5) with f ( x ) as a continuous homogeneous vector field of degree κ < 0 with respect to r 1 , r 2 , , r m . If the system (5) is globally asymptotically stable, for any initial, x ( 0 ) m it converges to origin in finite time.

3. Problem Description

In this section, the problem of cooperative guidance is described. The problem we consider is that several unmanned aerial vehicles (UAVs) attack a target in an arbitrary maneuver. The following is a typical assumption when considering cooperative guidance problems in engineering practice.
The following assumptions are given:
Assumption 1.
The topology of directed communication between agents is strongly connected and detail balanced.
Assumption 2.
The disturbance satisfies the following conditions:
w ˙ r i     w 1 w ˙ q i     w 2
where w r i and w q i denote the disturbance caused by the target maneuver in the direction of the line of sight and the component normal to the line of sight, respectively. Meanwhile w 1 and w 2 are known constants respectively.
Assumption 3.
The seeker and autopilot dynamics of the missiles are fast enough in comparison with the guidance loop.
The geometric relationships in the guidance process are illustrated in Figure 1 where m i denotes the i th UAV, and T denotes the target.
The terms q i ,   η i ,   θ i stand for the LOS angle, leading angle, flight path angle, respectively. From the geometric relations in Figure 1, we can obtain that
η i = q i θ i .
We note that a case η i   0 , π , ( i . e . ,   θ i ( t ) q i ( t ) ) can be considered as an asymmetry of a case η i   0 , π , ( i . e . ,   θ i ( t ) q i ( t ) ) . The case η i t     0 , π is treated in this paper. The case pursued is made by
r ˙ i = V T cos η T i V i cos η i r i q ˙ i = V i sin η i V T sin η T i θ ˙ i = n i y / V i V ˙ i = n i x
where r i is the distance from m i to the target; V i represents the velocity of m i , while n i x as well as n i y are the overloads of m i in its velocity frame, tuning the magnitude and direction of V i , and η T i , V T stands for the leading angle and the velocity of the target.
Taking the derivative of r ˙ i we can obtain that
r ¨ i = V ˙ T cos η T i V T sin η T i η ˙ T i V ˙ i cos η i + V i sin η i η ˙ i = V ˙ T cos η T i + V T sin η T i θ ˙ T V ˙ i cos η i V i sin η i θ ˙ i + r i q ˙ i 2 .
Denote that
w r i = V ˙ T cos η T i + V T sin η T i θ ˙ T u r i = V ˙ i cos η i + V i sin η i θ ˙ i
then we can obtain further
r ¨ i = w r i u r i + r i q ˙ i 2 .
The remaining flight time t g o i of the vehicle can be roughly estimated as
t g o i = r i r ˙ i .
Taking the derivative of (10) and substituting (9) we can obtain
t ˙ g o i = 1 + r i 2 q ˙ i 2 r ˙ i 2 r i r ˙ i 2 u r i + r i r ˙ i 2 w r i .
Taking the derivative of q ˙ i , we can obtain
q ¨ i = 1 r i 2 ( r i ( t ) ( V ˙ i sin η i + V i cos η i η ˙ i V ˙ T sin η T i   V T cos η T i η ˙ T i ) r i r ˙ i q ˙ i ) .
Denote that
w q i = V ˙ T sin η T i V T cos η T i θ ˙ T
and
u q i = V ˙ i sin η i V i cos η i θ ˙ i .
Note that η ˙ i = q ˙ i θ ˙ i , and η ˙ T i = q ˙ i θ ˙ T then we can get
q ¨ i = 1 r i ( u q i w q i ) 2 q ˙ i r ˙ i r i .
Illustrative Example. The purpose of the cooperative attack is to enable multiple UAVs to achieve a simultaneous attack on the target while ensuring their cooperation in terms of hitting angles. Simultaneous strikes can be achieved if and only if the remaining flight time achieves consensus. From (10), we can find that the remaining flight time of the UAV is determined by r i and r ˙ i . At the same time, we can see that if lim t T e n d q ˙ i = 0 then the UAV m i is able to hit the target. Then we are able to obtain the following definition.
Definition 2.
If the following equations hold simultaneously, the multi-UAVs system achieves a coordinated attack in finite time.
lim t T e n d ( ( q i c i ) ( q j c j ) ) = 0 lim t T e n d q ˙ i = lim t T e n d q ˙ j = 0 lim t T e n d ( t g o i t g o j ) = 0 i , j N
where c i and c j are default constants and T e n d is a limited adjustment time.

4. Main Result

This section proposes a sliding-mode FDO to evaluate the interference caused by target maneuvers in finite time. After that, based on the finite-time consensus theory, u r i and u q i are designed to enable simultaneous multi-directional attacks against the target, respectively.

4.1. FDO Design

Inspired by [29], the FDO is designed as
w ^ ˙ q i = α w ˜ q i 2 p q + β w ˜ q i p q + w 1 sgn w ˜ q i w ^ ˙ r i = α w ˜ r i 2 p q + β w ˜ r i p q + w 2 sgn w ˜ r i
then the observation error are:
w ˜ q i = w q i w ^ q i = r i ( q ¨ i 2 q ˙ i t g o i + 1 r i u q i ) w ^ q i w ˜ r i = w r i w ^ r i = r ˙ i x 3 i ( t ˙ g o i 1 + q ˙ i 2 t g o i 2 + t g o i r ˙ i u r i ) w ^ r i .
Theorem 1.
If Assumption 2 holds, and the parameters satisfy that α > 0 ,   β > 0 , q > p > 0 , then the disturbance observation error w ˜ q i and w ˜ r i are able to converge in finite time.
Proof of Theorem 1.
Define the Lyapunov function as:
V 1 i = 1 2 w ˜ q i 2 + 1 2 w ˜ r i 2 .
Taking the time derivative of the above function, we can get
V ˙ 1 i = w ˜ q i w ˜ ˙ q i + w ˜ r i w ˜ ˙ r i = w ˜ q i · ( w ˙ q i α w ˜ q i 2 p q β w ˜ q i p q w 1 sgn w ˜ q i ) +   w ˜ r i · ( w ˙ r i α w ˜ r i 2 p q β w ˜ r i p q w 2 sgn w ˜ r i ) w ˜ q i · ( α w ˜ q i 2 p q β w ˜ q i p q + w 1 ( 1 sgn w ˜ q i ) ) +   w ˜ r i · ( α w ˜ r i 2 p q β w ˜ r i p q + w 2 ( 1 sgn w ˜ r i ) ) α ( w ˜ q i 2 ) 3 q p 2 q β ( w ˜ q i 2 ) p + q 2 q α ( w ˜ r i 2 ) 3 q p 2 q β ( w ˜ r i 2 ) p + q 2 q .
According to Lemma 3 we can obtain
V ˙ 1 i 2 α ( V 1 i 2 ) 3 q p 2 q 2 β ( V 1 i 2 ) p + q 2 q = 2 ( α V 1 i q p q + β ) V 1 i p + q 2 q .
We can obtain α V 1 i ( p q ) / q > 0 from α > 0 ,   q > p > 0 then (18) can be rewritten as:
V ˙ 1 i 2 β V ( p + q ) / 2 p .
According to the finite-time stability theory in Lemma 2, if α > 0 ,   β > 0 ,   q > p > 0 , we can obtain that V 1 i 0 in a finite time T 1 . By the structure of V 1 i we can get that w ˜ q i 0 in a finite time T 1 , where T 1 q β ( q p ) V q p 2 q ( 0 ) . □

4.2. Impact Angle Cooperation Part

According to Definition 2, we get the purpose of impact angle cooperation is
lim t T e n d ( ( q i ( t ) c i ) ( q j ( t ) c j ) ) = 0 lim t T e n d q ˙ i ( t ) = lim t T e n d q ˙ j ( t ) = 0 .
This means that the purpose of the LOS angle cooperation is to make the LOS angles of different UAVs into a predefined sequence. The cooperative guidance law is also supposed to make the LOS angular rate of each UAV asymptotically equal to zero.
Let
θ 1 i = q i c i θ 2 i = q ˙ i .
It follows from Equation (12) that its derivative satisfies
θ ˙ 1 i = θ 2 i θ ˙ 2 i = 2 r ˙ i r i θ 2 i + 1 r i ( u q i w q i ) .
and the guidance law is designed as follows
u q i = 2 i r ˙ i θ 2 i + w ^ q i +   r i ( k 1 j = 1 N a i j sig ( θ 1 i θ 1 j ) α 1 + k 2 sig ( θ 2 i ) α 2 ) .
Lemma 5.
Under Assumption 1, there exists an appositive column vector τ = ( ς 1 , ς 2 , ς N ) T , such that ς i a i j = ς j a j i for all i , j = 1 , N . Denote ς = diag ( ς 1 , ς 2 , ς N ) then the matrix ς L is symmetric.
Theorem 2.
If the parameters in the guidance law (23) are satisfied k 1 > 0 ,   k 2 > 0 ,   0 < α 1 < 1 ,   α 2 = 2 α 1 1 + α 1 , then the impact angle coordination condition (20) is achieved in finite-time T 2 .
Proof of Theorem 2.
Substituting (23) into (22) yields
θ ˙ 1 i = θ 2 i θ ˙ 2 i = w ˜ q i r i k 1 j = 1 N a i j sig ( θ 1 i θ 1 j ) α 1 k 2 sig ( θ 2 i ) α 2 .
Let
ξ 1 = M θ 1 ξ 2 = θ 2
with θ 1 = ( θ 11 , θ 12 , , θ 1 N ) T , θ 2 = ( θ 21 , θ 22 , , θ 2 N ) T and M = I 1 χ T of which χ T = ( χ 1 , χ 2 , , χ N ) satisfies χ T L = 0 and i = 1 N χ i = 1 . According to the algebraic properties of M , we have ξ 1 = 0 and ξ 2 = 0 if and only if θ 11 = θ 12 = = θ 1 N , θ 21 = θ 22 = = θ 2 N = 0 respectively. From (24) and (25) we obtain
ξ ˙ 1 = M ξ 2 ξ ˙ 2 = k 1 sig ( L ξ 1 ) α 1 k 2 sig ( ξ 2 ) α 2 + w ˜ q
where w ˜ q = d i a g ( w ˜ q 1 r 1 , w ˜ q 2 r 2 , , w ˜ q N r N ) . Consider a Lyapunov function candidate
V 2 = k 1 0 L ξ 1 ς sig ( s ) α 1 d s + 1 2 ξ 2 T ( ς L ) ξ 2
where ς = diag ( ς 1 , ς 2 , ς N ) > 0 . Noting that L ξ 1 and sig ( L ξ 1 ) α 1 have the same component sign, we are able to obtain 0 L ξ 1 ς sig ( s ) α 1 d s > 0 for any ξ 1 0 . Furthermore, from ς L = ( ς L ) T we know that ξ 2 T ( ς L ) ξ 2 > 0 for any ξ 2 0 . Therefore, we obtain V 2 as a positive definition. Then the time derivative of the above function yields
V ˙ 2 = k 1 ( L ξ ˙ 1 ) T ς sig ( L ξ 1 ) α 1 + ξ 2 T ς L ξ ˙ 2 = k 1 ξ 2 T L T ς sig ( L ξ 1 ) α 1 k 1 ξ 2 T ς L sig ( L ξ 1 ) α 1 k 2 ξ 2 T ς L sig ( ξ 2 ) α 2 +   ξ 2 T ς L w ˜ q = k 2 ξ 2 T ς L sig ( ξ 2 ) α 2 + ξ 2 T ς L w ˜ q .
From the definition of q ˙ i in (8) and note that both V i and V T are bounded, we can get that ξ 2 is bounded in [ 0 , T 1 ) . Note that w ˜ q is also bounded in [ 0 , T 1 ) we can obtain that ξ 2 T ς L w ˜ q is bounded in [ 0 , T 1 ) and its upper bound is assumed to be w ¯ q . It follows from (28) that:
V ˙ 2 k 2 ξ 2 T ς L sig ( ξ 2 ) α 2 + w ¯ q w ¯ q .
Thus, V 2 is bounded in [ 0 , T 1 ) . When t T 1 , we get that w ˜ q = 0 . Then we can receive that
V ˙ 2 k 2 ξ 2 T ς L sig ( ξ 2 ) α 2 0 .
It can be seen that V ˙ 2 = 0 implies ξ 1 = ξ 2 = 0 . It can be seen from LaSalle’s invariance theorem that the system (26) can reach globally asymptotically stable for its zero equilibrium.
Next, we will prove that the system dynamics have a negative degree of homogeneity. Let ζ = ( ξ 11 , ξ 12 , , ξ 1 N , ξ 21 , ξ 22 , , ξ 2 N ) T R 2 N while the derivative of ζ is ζ ˙ = f ( ζ ) . Consider the dilation ( s 1 , s 1 , s 2 , s 2 ) , s 1 > 0 , s 2 > 0 and homogeneity κ then we can gain
f i ( ε s 1 ξ 11 , , ε s 1 ξ 1 N , ε s 2 ξ 21 , , ε s 2 ξ 2 N ) = ε s 2 ξ 2 i j = 1 N χ j ε s 2 ξ 2 j = ε s 2 ( ξ 2 i j = 1 N χ j ξ 2 j )
f N + i ( ε s 1 ξ 11 , , ε s 1 ξ 1 N , ε s 2 ξ 21 , , ε s 2 ξ 2 N ) = k 1 1 r i sig ( i = 1 N a i j ( ε s 1 ξ 1 i ε s 1 ξ 1 j ) ) α 1 k 2 1 r i sig ( ε s 2 ξ 2 i ) α 2 .
We can find for every ε > 0 we have f i ε s 1 ξ 1 , ε s 2 ξ 2 , , ε s 2 N ξ 2 N   = ε κ + s i f i ( x ) , i = 1 , 2 , , 2 N . By Lemma 4, we have κ < 0 holds by setting 0 < α 1 < 1 ,   α 2 = 2 α 1 1 + α 1 . According to Lemma 4, the system (26) can achieve global finite-time stability. Then the impact angle coordination condition (20) could be achieved in finite-time. We have completed the proof. □

4.3. Impact Time Cooperation Part

By Definition 2, the goal of time cooperation is
lim t T e n d ( t g o i t g o j ) = 0 .
According to this goal, we design time cooperative guidance law as
u r i = w ^ r i k 3 sig ( j = 1 N a i j ( t g o i t g o j ) ) α 3 .
Theorem 3.
If the parameters in the guidance law (32) are satisfied k 3 > 0 ,   0 < α 3 < 1 , then the impact angle coordination condition (31) could be achieved in finite time.
Proof of Theorem 3.
Substituting (32) into the system (11) we get
t ˙ g o i = 1 + q ˙ i 2 t g o i 2 k 3 r i r ˙ i 2 j = 1 N a i j sig ( t g o i t g o j ) α 3 t g o i r ˙ i w ˜ r i .
Let x = ( t g o 1 , t g o 2 , , t g o N ) T then (33) can be rewritten as
x ˙ = 1 + Φ k 3 r sig ( L x ) α 3 w ˜ r
where Φ = diag ( q ˙ 1 2 t g o 1 2 , q ˙ 2 2 t g o 2 2 , , q ˙ N 2 t g o N 2 ) , w ˜ r = ( t g o 1 r ˙ 1 w ˜ r 1 , t g o 2 r ˙ 2 w ˜ r 2 , , t g o N r ˙ N w ˜ r N ) and r = diag ( r 1 r ˙ 1 2 , r 2 r ˙ 2 2 , r N r ˙ N 2 ) . Let ζ = M x . According to the algebraic properties of M , we have ζ = 0 if and only if t g o 1 = t g o 2 = = t g o N , respectively. Take the time derivative of the above function yields:
ζ ˙ = k 3 r sig ( L ζ ) α 3 + M Φ M w ˜ r .
Consider a Lyapunov function candidate V 3 = 1 2 ζ T ς L ζ . It is convenient to obtain that V 3 is positive definite. Its derivative along (35) satisfies
V ˙ 3 = k 3 ζ T ς L M r sig ( L ζ ) α 3 ζ T ς L w ˜ r + ζ T ς L Φ .
Based on the definition of ζ and Φ , we are able to obtain that V 3 is bounded in t [ 0 , T 2 ) . From Theorems 1 and 2, we are able to get w ˜ r = 0 and Φ = 0 in t T 2 . Then (36) can be reformulated as
V ˙ 3 = k 3 ζ T L T ς M r sig ( L ζ ) α 3 .
Notice that r > 0 while L ζ and sig ( L ζ ) α 3 has the same sign component-wise, we can obtain V ˙ 3 0 and it can be seen that V ˙ 3 = 0 indicates that ζ = 0 . It can be seen from LaSalle’s invariance theorem that the system (35) can reach globally asymptotically stable for its zero equilibrium.
Using the same analysis as in Theorem 2, we are able to obtain that the system (35) has a negative degree of homogeneity when 0 < α 3 < 1 . Then the impact time coordination condition (31) can be achieved in a finite time. □
Remark 1.
For comparison, it is helpful to review the literature’s engagement models and guidance laws in References [24,25]. In Ref. [24] the derivative of t g o is given as:
t ˙ g o i = 1 + r i 2 q ˙ i 2 r ˙ i 2 r i r ˙ i 2 u r i .
A three-dimensional version is given in the literature [25] as
t ˙ g o = 1 + r i 2 r ˙ i 2 q ˙ ε 2 + r i 2 r ˙ i 2 q ˙ β 2 cos 2 q ε r i r ˙ i 2 u r i + d r .
By feedback linearization, the guidance law was designed in the literature [24,25] as
u r i = r i q ˙ i 2 r ˙ i 2 r i u ˜ r i
and
u r i = r i q ˙ ε i 2 + r i q ˙ β i 2 cos 2 q ε i r ˙ i 2 r i u ˜ r i
t ˙ g o i can be rewritten as t ˙ g o i = 1 + u ˜ ri in this way. References [24,25] designed consensus-based DCGL to achieve simultaneous arrival based on this model. However, the overload command has critical singularities when r i 0 . Due to the singularities, a r i would diverge to infinity when r i 0 . The guidance law (32) does not use feedback linearization and therefore avoids singularities above.

5. Numerical Example

We verify the performance of the DCGL by a numerical example of attacking a maneuvering target by four UAVs. In which the speed of the target is 300 m/s. The target’s normal acceleration is set to 20 cos ( 0.2 t ) m / s 2 . The communication topology of the vehicles is shown in Figure 2.
The initial conditions of multiple UAVs are shown in Table 1.
For the FDO, the parameters are set to be α = 1 ,   β = 1 ,   p = 1 ,   q = 2 ,   w 1 = w 2 = 10 . The parameters of the guidance law in the normal direction of LOS are set as k 1 = 0.5 ,   k 2 = 0.08 ,   α 1 = 0.5 ,   α 2 = 0.667 . Further, the guidance law parameters for the LOS direction are chosen as k 3 = 2 ,   α 3 = 0.5 .

5.1. Simulation of Pursuit and Head-On Attacks

We set up two sets of simulations, where group A is a chase attack and group B is a head-on attack. The initial position of the target in group A is (12,000, 0) and the initial orientation is the initial heading angle of the target is 0°. Group B target’s initial position is (25,000, 0) the initial heading angle of the target is 180°.
The simulation curves of the suggested cooperative guidance law, including the trajectories of the vehicles and the target, time-to-go of four UAVs to attack a maneuvering target, impact angle, LOS angular rate, radial relative velocities, actual and estimated values of external disturbance, and tangential and normal acceleration command, are shown in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10.
In Figure 3 and Figure 4, it can be seen that the four UAVs achieve a simultaneous attack on a maneuvering target. In addition, it can be observed in Figure 4 that the four UAVs have different remaining flight times at moment zero. However, using the designed finite-time consensus protocol (34), the remaining flight time reaches the same value after approximately 10 seconds. Figure 5 shows that the relative velocity between the UAVs and the target converges to fixed values after some time. At the same time, the relative velocities of the different UAVs and the target do not converge to the same value, which avoids unnecessary energy consumption. Figure 6 and Figure 7 indicate that under the angular cooperative guidance law (25), the LOS angle converges to the desired sequence after some time, while the LOS angle rate converges to zero.
As shown in Figure 8 and Figure 9, the proposed FDO (15) has a favorable performance in estimating the disturbance caused by the unknown maneuvers of the target.
Figure 10 and Figure 11 show that both the tangential acceleration commands and normal acceleration commands are smooth.

5.2. Comparison with Feedback Linearization Methods

We compared the tangential acceleration commands of the two methods for the initial condition A. From the Figure 12 we can see that the required overload command of the proposed method is significantly smaller. From the right-hand image, we can see that the required acceleration at the initial moment is large at around 4 g . At the same time, the acceleration required for the proposed method is less than 1 g . This is due to the fact that r i 2 q ˙ i 2 r ˙ i 2 needs to be compensated directly in the feedback linearization method. Usually, q ˙ i does not converge to zero at the initial moment and r i is large at the same time, therefore resulting in a large overload command. At the same time, we can see that the feedback linearization method has a significantly higher acceleration command at the moment of hit. This is because the acceleration command includes r ˙ i 2 r i u ˜ r i . This leads that the overload command has critical singularities when r i 0 . Through comparison, it can be seen that the proposed method in this paper has a smaller overload command and is more promising for application.
Figure 12. Comparison of overload commands.
Figure 12. Comparison of overload commands.
Aerospace 09 00157 g012

6. Conclusions

This paper proposed a new cooperative guidance law that can be applied to targets with large maneuvers. We propose a new sliding mode FDO to estimate the rapidly changing acceleration of the target. Meanwhile, the guidance law was designed in two parts: angular cooperation and time cooperation to satisfy impact time and impact angle constraints. We adopt the nonlinear guidance law design to avoid numerical singularity caused by feedback linearization. A rigorous derivation demonstrates that the proposed guidance law can achieve a finite-time cooperative attack under directed communication topology. Numerical simulations show that the proposed guidance law could achieve a cooperative attack on a highly maneuverable target while avoiding the high tangential overload command of the UAVs. Future work will attempt to consider autopilot dynamics in the design of cooperative guidance law for actual engagement missions.

7. Annexes

A directed graph G = ( V , E ) is developed for N agents to represent the interactions between agents, the vertex set and the edge set. Moreover, edges are an ordered pair of vertices V , implying that agent j can receive information from agent i. If a directed edge from i to j exists, then i would be defined as the parent node, and j would be defined as the child node, the neighbors of node i are represented by N i = j V | v j , v i E and N i are the neighbor numbers of agent i.
The adjacency matrix A associated with G is defined such that a i j = 1 if i j and node i is neighboring to node j, while a i j = 0 otherwise. For a directed graph G , a i j a j i and therefore the matrix A is asymmetric. The Laplacian matrix of the graph associated with adjacency matrix A is given as L = l i j , where l i i = j = 1 N a i j and l i j = a i j ,   i j .

Author Contributions

Writing—original draft X.D.; Writing—review & editing Z.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Technology Innovation 2030-Key Project of “New Generation Artificial Intelligence” under Grant 2020AAA0108200, the National Natural Science Foundation of China under Grants 61873011, 61922008, 61973013 and 61803014, the Defense Industrial Technology Development Program under Grant JCKY2019601C106, the Innovation Zone Project under Grant 18-163-00-TS-001-001-34, the Foundation Strengthening Program Technology Field Fund under Grant 2019-JCJQ-JJ-243, and the Fund from Key Laboratory of Dependable Service Computing in Cyber Physical Society under Grant CPSDSC202001.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jeon, I.S.; Lee, J.I.; Tahk, M.J. Impact-time-control guidance law for anti-ship missiles. IEEE Trans. Control Syst. Technol. 2006, 14, 260–266. [Google Scholar] [CrossRef]
  2. Kim, M.; Jung, B.; Han, B.; Lee, S.; Kim, Y. Lyapunov-based impact time control guidance laws against stationary targets. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 1111–1122. [Google Scholar] [CrossRef]
  3. Saleem, A.; Ratnoo, A. Lyapunov-based guidance law for impact time control and simultaneous arrival. J. Guid. Control. Dyn. 2016, 39, 164–173. [Google Scholar] [CrossRef]
  4. Tekin, R.; Erer, K.S.; Holzapfel, F. Polynomial shaping of the look angle for impact-time control. J. Guid. Control. Dyn. 2017, 40, 2668–2673. [Google Scholar] [CrossRef]
  5. Wang, Y.; Dong, S.; Ou, L.; Liu, L. Cooperative control of multi-missile systems. IET Control Theory Appl. 2014, 9, 441–446. [Google Scholar] [CrossRef]
  6. Zhao, J.; Zhou, R. Obstacle avoidance for multi-missile network via distributed coordination algorithm. Chin. J. Aeronaut. 2016, 29, 441–447. [Google Scholar] [CrossRef] [Green Version]
  7. Jiang, H.; An, Z.; Chen, S.; Xiong, F. Cooperative guidance with multiple constraints using convex optimization. Aerosp. Sci. Technol. 2018, 79, 426–440. [Google Scholar] [CrossRef]
  8. Yao, Z.A.; St, B.; Jie, G.B. Two-stage cooperative guidance strategy using a prescribed-time optimal consensus method. Aerosp. Sci. Technol. 2020, 100, 105641. [Google Scholar]
  9. Chen, Y.; Wang, J.; Shan, J.; Xin, M. Cooperative guidance for multiple powered missiles with constrained impact and bounded speed. J. Guid. Control. Dyn. 2021, 44, 1–17. [Google Scholar] [CrossRef]
  10. Li, G.; Wu, Y.; Xu, P. Adaptive fault-tolerant cooperative guidance law for simultaneous arrival. Aerosp. Sci. Technol. 2018, 82, 243–251. [Google Scholar] [CrossRef]
  11. Lyu, T.; Li, C.; Guo, Y.; Ma, G. Three-dimensional finite-time cooperative guidance for multiple missiles without radial velocity measurements. Chin. J. Aeronaut. 2019, 32, 241–251. [Google Scholar] [CrossRef]
  12. Nikusokhan, M.; Nobahari, H. Closed-form optimal cooperative guidance law against random step maneuver. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 319–336. [Google Scholar] [CrossRef]
  13. Zhao, J.; Zhou, R.; Dong, Z. Three-dimensional cooperative guidance laws against stationary and maneuvering targets. Chin. J. Aeronaut. 2015, 28, 1104–1120. [Google Scholar] [CrossRef] [Green Version]
  14. Zhao, J.; Zhou, R. Unified approach to cooperative guidance laws against stationary and maneuvering targets. Nonlinear Dyn. 2015, 81, 1635–1647. [Google Scholar] [CrossRef]
  15. Zhao, Q.L.; Dong, X.W.; Song, X.; Ren, Z. Cooperative time-varying formation guidance for leader-following missiles to intercept a maneuvering target with switching topologies. Nonlinear Dyn. 2019, 95, 129–141. [Google Scholar] [CrossRef]
  16. Zhou, J.; Lü, Y.; Li, Z.; Yang, J. Cooperative guidance law design for simultaneous attack with multiple missiles against a maneuvering target. J. Syst. Sci. Complex. 2018, 31, 287–301. [Google Scholar] [CrossRef]
  17. Shaferman, V.; Shima, T. Cooperative optimal guidance laws for imposing a relative intercept angle. J. Guid. Control Dyn. 2015, 38, 1395–1408. [Google Scholar] [CrossRef]
  18. Shaferman, V.; Shima, T. Cooperative differential games guidance laws for imposing a relative intercept angle. J. Guid. Control Dyn. 2017, 40, 2465–2480. [Google Scholar] [CrossRef]
  19. Wang, X.H.; Tan, C.P. 3-D impact angle constrained distributed cooperative guidance for maneuvering targets without angular-rate measurements. Control Eng. Pract. 2018, 78, 142–159. [Google Scholar] [CrossRef]
  20. Yu, J.; Dong, X.; Li, Q.; Ren, Z. Distributed cooperative encirclement hunting guidance for multiple flight vehicles system. Aerosp. Sci. Technol. 2019, 95, 105475. [Google Scholar] [CrossRef]
  21. Wang, X.; Zhang, Y.; Liu, D.; He, M. Three-dimensional cooperative guidance and control law for multiple reentry missiles with time-varying velocities. Aerosp. Sci. Technol. 2018, 80, 127–143. [Google Scholar] [CrossRef]
  22. Chen, Z.Y.; Chen, W.C.; Liu, X.M.; Cheng, J. Three-dimensional fixed-time robust cooperative guidance law for simultaneous attack with impact angle constraint. Aerosp. Sci. Technol. 2021, 110, 106523. [Google Scholar] [CrossRef]
  23. Liu, F.; Dong, X.; Li, Q.; Ren, Z. Robust multi-agent differential games with application to cooperative guidance. Aerosp. Sci. Technol. 2021, 111, 106568. [Google Scholar] [CrossRef]
  24. Lin, M.; Ding, X.; Wang, C.; Cui, N. Fixed-time cooperative guidance law with angle constraint for multiple missiles against maneuvering target. IEEE Access 2020, 1, 99. [Google Scholar]
  25. Song, J.; Song, S.; Xu, S. Three-dimensional cooperative guidance law for multiple missiles with finite-time convergence. Aerosp. Sci. Technol. 2017, 67, 193–205. [Google Scholar] [CrossRef]
  26. Zhai, S.; Wei, X.; Yang, J. Cooperative guidance law based on time-varying terminal sliding mode for maneuvering target with unknown uncertainty in simultaneous attack. J. Frankl. Inst. 2020, 357, 11914–11938. [Google Scholar] [CrossRef]
  27. Bhat, S.; Bernstein, D. Finite time stability of continuous autonomous systems. SIAM J. Control Optim. 2000, 38, 751–766. [Google Scholar] [CrossRef]
  28. Zuo, Z.; Tie, L. Distributed robust finite-time nonlinear consensus protocols for multi-agent systems. Int. J. Syst. Sci. 2016, 47, 1366–1375. [Google Scholar] [CrossRef]
  29. Ran, M.; Wang, Q.; Dong, C. Stabilization of a class of nonlinear systems with actuator saturation via active disturbance rejection control. Automatica 2016, 63, 302–310. [Google Scholar] [CrossRef]
Figure 1. The geometric relationship between multiple UAVs and targets.
Figure 1. The geometric relationship between multiple UAVs and targets.
Aerospace 09 00157 g001
Figure 2. Communication topology among UAVs.
Figure 2. Communication topology among UAVs.
Aerospace 09 00157 g002
Figure 3. Trajectories of four UAVs to attack a maneuvering target simultaneously.
Figure 3. Trajectories of four UAVs to attack a maneuvering target simultaneously.
Aerospace 09 00157 g003
Figure 4. Time-to-go for four UAVs to attack a maneuvering target.
Figure 4. Time-to-go for four UAVs to attack a maneuvering target.
Aerospace 09 00157 g004
Figure 5. Radial relative velocities between four UAVs and a maneuvering target.
Figure 5. Radial relative velocities between four UAVs and a maneuvering target.
Aerospace 09 00157 g005
Figure 6. The impact angle of the UAVs.
Figure 6. The impact angle of the UAVs.
Aerospace 09 00157 g006
Figure 7. LOS angular rate of the UAVs.
Figure 7. LOS angular rate of the UAVs.
Aerospace 09 00157 g007
Figure 8. Actual and estimated values of external disturbance in the direction of LOS.
Figure 8. Actual and estimated values of external disturbance in the direction of LOS.
Aerospace 09 00157 g008
Figure 9. Actual and estimated values of external disturbance in the normal direction of LOS.
Figure 9. Actual and estimated values of external disturbance in the normal direction of LOS.
Aerospace 09 00157 g009
Figure 10. Normal overload commands.
Figure 10. Normal overload commands.
Aerospace 09 00157 g010
Figure 11. Tangential overload commands.
Figure 11. Tangential overload commands.
Aerospace 09 00157 g011
Table 1. Initial conditions of multiple UAVs.
Table 1. Initial conditions of multiple UAVs.
UAVInitial Position (m,m)Initial Velocity/m/sHeading Angle/°Desired Angle/°
1(0, −2000)5901515
2(0, −1000)6001010
3(0, 1000)610−105
4(0, 2000)590−150
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dong, X.; Ren, Z. Finite-Time Distributed Leaderless Cooperative Guidance Law for Maneuvering Targets under Directed Topology without Numerical Singularities. Aerospace 2022, 9, 157. https://doi.org/10.3390/aerospace9030157

AMA Style

Dong X, Ren Z. Finite-Time Distributed Leaderless Cooperative Guidance Law for Maneuvering Targets under Directed Topology without Numerical Singularities. Aerospace. 2022; 9(3):157. https://doi.org/10.3390/aerospace9030157

Chicago/Turabian Style

Dong, Xiaofei, and Zhang Ren. 2022. "Finite-Time Distributed Leaderless Cooperative Guidance Law for Maneuvering Targets under Directed Topology without Numerical Singularities" Aerospace 9, no. 3: 157. https://doi.org/10.3390/aerospace9030157

APA Style

Dong, X., & Ren, Z. (2022). Finite-Time Distributed Leaderless Cooperative Guidance Law for Maneuvering Targets under Directed Topology without Numerical Singularities. Aerospace, 9(3), 157. https://doi.org/10.3390/aerospace9030157

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