Next Article in Journal
Control of Cowl Shock/Boundary Layer Interaction in Supersonic Inlet Based on Dynamic Vortex Generator
Previous Article in Journal
Testing of the N2O/HDPE Vortex Flow Pancake Hybrid Rocket Engine with Augmented Spark Igniter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Dynamic Cross-Collaborative Interception Algorithm Based on GTSMC and Virtual Geometry

1
School of Astronautics, Beijing Institute of Technology, Beijing 100081, China
2
Shanghai Institute of Mechanical and Electrical Engineering, Shanghai 201109, China
3
Aerospace College, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(8), 728; https://doi.org/10.3390/aerospace10080728
Submission received: 10 March 2023 / Revised: 7 June 2023 / Accepted: 14 June 2023 / Published: 20 August 2023

Abstract

:
In the model (m:n), to improve the autonomous collaborative interception capability for air vehicle, a new autonomous cross-collaborative interception algorithm based on GTSMC (Global Terminal Sliding Mode Control) and real-time virtual geometry is proposed in this paper. Firstly, the conception of an autonomous cross-collaboration is defined and the multi-air vehicle for the multi- object interception problem is formulated. Then, this paper presents the dynamic situation assessment function, which considers the real-time flight status and cooperative status of the air vehicle during the interception of the object. At the same time, this paper states the condition of whether the air vehicle is in a cooperative state and proves it. After completing the dynamic situation assessment, and considering the dynamic of the air vehicles, a new controller is designed by using GTSMC and the idea of backstepping method. Simultaneously, this paper gives a stability analysis of the closed-loop system by using Lyapunov theory. Finally, to demonstrate the effectiveness of the proposed algorithm, several simulation cases which consider different interception scenarios are given. The simulation results show that the new collaborative interception algorithm can provide better autonomous cross-collaborative interception capability and higher accuracy.

1. Introduction

Recently, in the wake of developments in aerospace science and technology, the research on cooperative interception has become a hotspot of multi-air vehicle to multi-object interception research all over the world. Facing the problem of multi-air vehicles and multi-objects scenarios by using cooperative interception strategy can greatly enhance the interception effectiveness [1] and improve the overall effectiveness [2]. In a dynamic environment, real-time assessment and object autonomous allocation play a crucial role. To solve the abovementioned problem, many previous works have been considered. For instance, ref. [3] designed a new WTA algorithm for a multiagent TMD scenario by considering the relative distance, angle of direction, and velocity aspects, as well as cooperative linear guidance laws for the object and the allocated defenders. Considering the terminal time and angle constraints, the look angle shaping technique is proposed in [4]. Based on Multi-source information fusion, especially the devices’ performance and the energy, the work in [5,6] investigates the object allocation problem. Considering the impact time factors, a time-to-go prediction formula and a desired error dynamic are proposed in [7]. In the work of [8], a two-stage cooperative strategy to address the optimal decentralized three-dimensional cooperative guidance problem is proposed. This method considers the problem of the optimal consensus control under directed and periodical switching topologies during the first phase. To solve the cooperative salvo attack problem, a new fixed finite-time consensus algorithm based on sliding mode control is proposed in [9]. In this approach, the cooperative salvo problem was transformed to a consensus problem over a cycle digraph, which aids in designing a guidance strategy through the selection of heterogeneous gains. Furthermore, in [10], a new cooperative guidance approach by using the receding horizon control technique is proposed to coordinate the impact time of a group of interceptors against the stationary object. To realize the simultaneous attack, the cooperative guidance problem for multiple interceptors with fixed and switching directed communication topologies was investigated in [11]. Considering the object with unknown maneuverability, a cooperative salvo guidance law by using asymptotic consensus over undirected graphs was presented in [12]. In the case of salvo attacks which require interceptors to hit the object simultaneously, a cooperative guidance algorithm based on the impact-time-control guidance (ITCG) law was proposed in [13]. For this cooperative guidance scheme, the coordination algorithms and local guidance laws were merged. Considering the desired impact angles, an integrated distributed cooperative guidance and control law for multi-air vehicle to attack a single object is proposed in [14]. Considering the condition of the actuation failures, a new fault-tolerant cooperative guidance under partial actuator capability was presented in [15]. However, this method mainly focuses on the stationary object, instead of the maneuvering objects. To counter the aircraft object, a high-order SMC control law was designed in [16]. However, the higher order of the system is, the more complex the computational effort of the high-order sliding mode surface will become. This is detrimental to the air vehicle’s ability to track the object in real time. In the work of [17], a novel cooperative guidance law for the task of salvo attack based on efficient information was investigated. In this approach, the effective time of the efficient information was regarded as an important parameter to indicate the ability of an air vehicle to attack a given stationary object. At the same time, this parameter was estimated by modifying the current time-to-go estimation method in the presence of a lateral acceleration constraint, while these guidance law often assume that the object is stationary. Considering the condition of communication delay, a new salvo attack guidance law designing for multiple interceptors against a stationary object is presented in [18]. In [19], the finite time simultaneous attack problem was investigated, when intercepting the maneuvering object with unknown acceleration. However, this guidance strategy uses radial acceleration to achieve consensus in time-to-go estimates, which may be difficult to implement in practice. By introducing the bias term into both pitch and yaw channels, a three-dimensional PNG based impact time control guidance law was proposed in [20]. However, this PNG guidance law used linear engagement kinematics or estimated time-to-go, which may generate large errors. In [21,22,23], it was found that many impact time control guidance laws based on SMC methods have complicated structures, which is very stressful to deal with the look angle constraint. Moreover, to satisfy impact time constraint, guidance gains or parameters are often tuned by trial and error, or by using an optimization routine, which can make on-line calculations less efficient. In addition, a fixed-time nonlinear circular guidance law that satisfies the impact time constraint is proposed in [24]. Using the geometric principle that the length of a circular arc connecting the air vehicle and the object can be analytically calculated, the exact expression of time-to-go can be obtained. Thus, the impact time error can be shaped to zero and the air vehicle can intercept the object at the desired time, which is crucial in a salvo attack. To improve the multiple-air vehicle cooperative attack capability and penetration capability, two three-dimensional impact-angle-constrained cooperative guidance strategies against the maneuvering object was proposed in [25]. In the work of [26], which assumes that the speed of all interceptors is the same, a fixed-time cooperative guidance algorithm by utilizing optimal control theory to achieve the cooperative attack on a stationary object was investigated. Based on the basis of proportional guidance, a two-stage cooperative guidance algorithm was proposed in [27]. In this approach, the first stage was to adjust the lead angle and the relative distance, and when the initial condition of the second stage was satisfied, the proportional guidance algorithm was selected to complete the cooperative attack on the stationary object. The work of [28] investigated the problem of multi-air vehicle cooperative attack with impact angle constraints by using proportional guidance. However, this approach often uses the linearized cooperative guidance model, and the guidance model is often nonlinear in real intercepting scenarios. Thus, this research is very suitable for attacking a stationary or slow-moving object. In [29], a cooperative guidance algorithm in three-dimensional space is presented. This approach both can satisfy the stationary object and constant-value maneuvering object. However, this method has a fatal disadvantage that the stability analysis was not given. To hit the stationary ground objects in a specified direction, a nonlinear impact angle control guidance law based on Lyapunov stability theory is proposed in [30]. While this method just considers how to intercept the stationary ground objects. When intercepting the maneuvering objects, this approach has not been investigated. By introducing a novel alternative strategy to ensure the pointwise solvability of SDRE, a new three-dimensional guidance law was designed in [31]. However, this research did not take any attack constraint into account. Based on the fusion of finite-time SDRE approach and integral SMC method, a finite and fixed time convergent impact angle guidance was derived in [32]. The main idea of this approach is that it converts the impact angle problem to a tracking problem. However, the existence of SDRE solutions was not guaranteed and numerical methods were involved, which restricted the implementation. Therefore, designing an effective and precise guidance law with impact angle constraint based on the SDRE technique is worthy of further investigation. Considering the interception of the maneuverable object without knowing the object’s information, a new adaptive backstepping integrated guidance law was presented in [33].
In summary, many studies have been conducted to solve the problem of multi-air vehicle to multi-object interception. However, several problems are not fully considered, for instance, most previous work often considers the particle model, instead of considering the dynamics of the air vehicle. Secondly, most current collaborative algorithms often focus on multi-to-one problems in a static environment. On the one hand, multi-to-multi situation is rarely researched, and, on the other hand, most collaborative methods do not consider the problem of switching the objects autonomously when considering multi-to-multi system. Thirdly, considering the cooperation strategy, the traditional situation assessment model only considers the relative position and velocity between the object and air vehicle. However, the cooperation state between the object and air vehicle is rarely considered. Finally, there is no method to determine whether the air vehicle is in a cooperative relationship or not. Thus, a new cross-collaborative interception algorithm based on GTSMC and attack geometry is proposed. This paper mainly focusses on the problem of multi-air vehicle interception of multi-object, especially in a dynamic interception scenario. The main contributions of this study are as follows:
  • The paper gives the definition of the cross-collaborative interception autonomously. When considering the problem of multi-air vehicle to multi-object interception in a dynamic environment, the paper takes the dynamics of the air vehicle into account, instead of using the particle model.
  • The paper establishes a new dynamic situation assessment model. The model consists of the flight statues and the cooperation state between each air vehicle and the object. At the same time, according to the attack geometry, the paper gives the cooperative conditions of the air vehicle and proves it.
  • Based on GTSMC and the idea of backstepping method, the paper designs a new controller to intercept the assigned object. Then, Using the Lyapunov theory, the closed-loop system was proved to be stable.
  • Finally, several simulation cases which consider different dynamic interception scenarios are given to demonstrate the effectiveness of the proposed algorithm.
Finally, the remainder of this paper is organized as follows: the introduction is stated in Section 1. Then, formulation of the problem is given in Section 2. The dynamic situation assessment, which considers the flight statues and the cooperation status, is designed in Section 3. In Section 4, the controller is designed by combining GTSMC and backstepping method. At the same time, in this section the stability of the closed-loop system is proved using Lyapunov theory. The simulation results and analysis are presented in Section 5. Finally, the conclusions are summarized in Section 6.

2. Problem Formulation

Equations of Kinematics and Air vehicle Dynamics

Before establishing the dynamic of the air vehicle–object, the following assumptions and definitions will be considered when analyzing and designing the interception law.
Definition 1.
Cross-collaborative interception: In a multi-air vehicle and multi-object interception scenario, the air vehicle can automatically switch to intercept the object, according to the flight statues situation of the air vehicle and the object and the cooperation status of the air vehicle when intercepting the same object.
Assumption 1.
The speed of the air vehicle and the object is constant.
Assumption 2.
For all objects, the paper defines that they are all unpowered at the terminal phase of their trajectories.
According to the above assumptions, the paper will consider the cross-cooperative guidance problem between n air vehicles and m objects in two-dimensional space, which is shown in Figure 1. In this interception scenario, the coordinate for each interceptor and each object is defined as M x m i , y m i , i = 1 , 2 , 3 , , n and T x t j , y t j , j = 1 , 2 , 3 , , m respectively. r i j denotes the relative distance between the i t h air vehicle and the j t h object. v t j , θ t j and a t j denotes the velocity, the flight path angle and the normal acceleration of the j t h object, respectively. Similarly, v m i , θ m i and a m i represent the velocity, the flight path angle, and the normal acceleration of the i t h air vehicle, respectively. q i j is the LOS angle between the i t h air vehicle and the j t h object. O m i x b i is the body axis of the i t h air vehicle. At last, ϑ i is the pitch angle of the i t h air vehicle.
According to all above conditions, the air vehicle–object engagement dynamics are defined as follows.
x ˙ t j = v t j cos θ t j y ˙ t j = v t j sin θ t j θ ˙ t j = a t j / v t j x ˙ m i = v m i cos θ m i y ˙ m i = v m i sin θ m i θ ˙ m i = a m i / v m i
r ˙ i j = v t j cos θ t j q i j v m i cos θ m i q i j r i j q ˙ i j = v t j sin θ t j q i j v m i sin θ m i q i j
Then, differentiating Equation (2), with respect to time, the LOS dynamics can be obtained as:
q ¨ i j = 2 r ˙ q ˙ i j r i j + a t j cos ϕ t j r i j a m i cos ϕ m i r i j
where ϕ t j and ϕ m are defined as ϕ t j = θ t j q i j , ϕ m i = θ m i q i j .
At the same time, the nonlinear dynamic model of the air vehicle considered in this paper can be written as follows.
a m i = n y g = 57.3 q s c y α α m g cos θ m i ϑ ˙ m i = ω z J z ω ˙ z = 57.3 q s l m z α α + q s l 2 m z ω z v m i ω z + 57.3 q s l m z δ z δ z α m i = ϑ m i θ m i n ˙ y = 57.3 q s c y α m g ω z + g sin θ m i v m i 57.3 q s c y α m v m i n y
According to Equation (4), m , ω z and J z denote the mass, the pitch angular rate, and the moment of inertia of the air vehicle. g is the gravitational acceleration. At the same time, q = ρ v m 2 / 2 . s denotes the reference area and l is the reference length. δ z denotes the rudder deflection angle. In addition, c y α denotes the lift force derivative with respect to α . m z α , m z ω z , and m z δ z denotes the pitch moment.

3. Dynamic Situation Assessment

In this section, the paper will design the dynamic situation assessment function to improve the collaborative performance, according to the dynamic situation. To reflect the real interception scenario, we will consider two main factors. These are the flight status of the air vehicle and the potential cooperation state between the air vehicle and the object. In terms of the first factor, we will mainly consider the relative position, radial relative velocity, and the acceleration between a pair of an air vehicle and its object. In terms of the second factor, this paper mainly considers the intercepting triangle between the cooperation air vehicle and the object.

3.1. Flight Status Situation Assessment Function

In this section, we will construct a situation assessment model between the air vehicle and the object. Firstly, a pair of the i t h air vehicle and the j t h object is defined as χ m i , T j in this paper.
As is widely known, during the interception process between multi-air vehicle and multi-object, the relative position between the air vehicle and the object is the primary consideration. In this paper, the position advantage function between the air vehicle and the object is defined as follows.
χ r i j m i , T j = e r i j σ r
According to Equation (5), it is clearly that the smaller relative distance between the i t h air vehicle and the j t h object is, the better for the i t h air vehicle intercepts.
In addition to the relative distance factor, the radial relative velocity between the i t h air vehicle and the j t h object is also important for the interception. Thus, the radial velocity advantage function can be defined as follows.
χ r v m i , T j = s i g m o i d v m i cos θ m i q i j v t j cos θ t j q i j σ r v 1 2
where σ r v denotes a positive constant. The radial velocity advantage function reflects the radial approach characteristics between the air vehicle and the object. The faster they approach, the better it is better for the air vehicle. Additionally, we also conclude that when χ r v m i , T j is positive, the relative distance decreases and the air vehicle is approaching the object. When χ r v m i , T j is negative, the air vehicle is way off object.
Besides the above factors, we also consider the variation of q i j between the i t h air vehicle and the j t h object in this paper. As is widely known, the smaller q ˙ i j is, the smoother the interception trajectory is. This is very advantageous to the interceptor. Thus, the q i j advantage function is defined as follows.
χ q i j m i , T j = cos p i 2 + s i g m o i d q i j 1 2

3.2. Cooperation Status Function

In this section, the paper will consider the cooperation state between each air vehicle when intercepting the same object during the process.
Theorem 1.
The air vehicles are in a state of cooperation if there is at least one intersection point between the air vehicles and the object in the direction of velocity.
Proof. 
As previously stated, the paper will take the two air vehicles M i , M i + 1 and one object T j as an example and use these two air vehicles and one object to form a triangle. In this triangle, the velocities of these are v M 1 , v M 2 and v T . The types of intercepting triangle at a certain moment are shown in Figure 2. □
Figure 2 shows some types of intercepting triangle. Now, we will take one of them as an example to prove Theorem 1. As shown in Figure 3, after connecting the air vehicle and the object, the paper rotates the connecting line q i j degrees to the right as follows.
According to Figure 3, we can establish the following line equation between the i t h air vehicle and the j t h object.
y M i = tan q i j θ m i x m i y T j = tan p i q i j + θ T j x T j x , 0 < q i j < p i
Assuming that there is an intersection point in the direction of velocity, we will get the following equation.
tan p i q i j + θ T j x T j x = tan q i j θ m i x m i
Further,
tan p i q i j + θ T j x T j + tan q i j θ m i x m i tan p i q i j + θ T j = x
where x T j , x m i , θ m i and q i j are known. Thus, we can find a value, x , which will satisfy Equation (10). Thus, the air vehicles are in a state of cooperation.
Similarly, we can prove the other relationship between the i + 1 t h air vehicle and the j t h object. Thus, the assessment model of potential cooperation is defined as:
g T j m i , m i + 1 = tanh S T j _ m i _ m i + 1
where S T j _ m i _ m i + 1 is the area of the triangle formed by m i , m i + 1 is the air vehicle and the T j object.
According to Equations (5)–(11), the total dynamic assessment function for the air vehicle and the object is defined as follows.
J = i n α 1 χ r i j m i , T j + α 2 χ r v m i , T j + α 3 χ q i j m i , T j + α c j = 1 m g T j m i , m i + 1
where α 1 , α 2 , α 3 , α c are parameters.

4. Controller Design

Instead of using the particle model, the paper will consider the dynamic model. So, after completing the dynamic situation assessment, the air vehicle can get the suitable interception object for itself. Thus, this paper will design the appropriate control law to intercept the assigned object in this section. Before designing the control law, this paper will establish the IGC model in state space first.

4.1. Dynamic Model in State Space

To establish the dynamic mode in state space, the relevant parameters can be defined as follows:
p 1 = 57.3 q s c y α m v m , p 2 = 57.3 q s c y δ z m v m , p 3 = 57.3 q s l m z α J z p 4 = q s l 2 m z ω z J z v m , p 5 = 57.3 q s l m z δ z J z , p 6 = g sin θ m v m
When intercepting the object, this paper will pay more attention to how to intercept the object successfully. Therefore, this paper will focus on changing the air vehicle–object distance. Thus, we will define the following state vector.
X = x 1 , x 2 , x 3 , x 4 , x 5 = Δ e i j , Δ e ˙ i j , n y m i , ω z m i
where Δ e i j is defined as Δ e i j = x m i x T j + y m i y T j . The control input u = δ m i .
According to the above equations, the dynamic model with strict-feedback state equation can be written as follows.
x ˙ 1 = x 2 x ˙ 2 = f 2 x 2 + b 2 x 3 x ˙ 3 = f 3 x 3 + b 3 x 4 x ˙ 4 = f 4 x 4 + b 4 u
where
f 2 x 2 = 0 f 3 x 3 = p 6 p 1 x 3 f 4 x 4 = p 4 x 4 + p 3 α , b 2 = cos θ m i sin θ m i b 3 = p 1 v m / g b 4 = p 5

4.2. GTSMC Controller Design Based on Backstepping

As is known, the conventional SMC method often chooses a linear sliding mode surface which has a fatal disadvantage: the tracking errors cannot converge to zero in a finite time; meanwhile, in this paper, we hope that the air vehicle–object distance can converge to zero within a designated time, when intercepting the object. Thus, a new modified global terminal sliding mode surface will be designed.
Step 1: According to Equation (15), x 1 , x 2 are related to the Δ e , Δ e ˙ . To ensure the air vehicle–object distance Δ e = 0 and Δ e ˙ = 0 , the new sliding mode surface is designed as follows:
s 1 = c 1 x 1 + c 2 x 1 ς / μ + c 3 x 2
where x 1 , x 2 are the state variable. c 1 , c 2 , c 3 are positive constants. ς and μ are positive odd integers which satisfy 1 < ς / μ < 2 .
The reaching law is defined as follows:
s ˙ 1 = k l s 1 k s s 1 λ / η
where λ and η are all positive odd constants. k l > 0 , k s > 0 determine the reaching speed.
Then, when differentiating Equation (17) is combined with Equations (15) and (18), we can get the following equation.
x 3 d = c 3 b 2 1 c 1 x 2 + c 2 ς / μ x 1 ς μ / μ x 2 + k l s 1 + k s s 1 λ / η
Step 2: To track the desired x 3 d , and avoid the “differential explosion” phenomenon, the paper defines the following equation.
x ˙ 3 = x ˙ 3 d + k 3 x 3 d x 3
Then, combining this with Equation (15), the desired pitch angular rate can be written as follows.
x 4 d = b 3 1 x ˙ 3 d + k 3 x 3 d x 3 f 3 x 3
Step 3: Similarly, we can get the desired deflection angle as follows.
u = b 4 1 x ˙ 4 d + k 4 x 4 d x 4 f 4 x 4

4.3. Stability Analysis

In this section, the paper will give the stability of the closed-loop system by using Lyapunov theory. Firstly, the paper defines the following tracking error vector.
E = e s , e x 3 , e x 4 = s 1 , x 3 x 3 d , x 4 x 4 d
Then, by defining Lyapunov function as follows:
V E = V e s + V e x 3 + V e x 4 = 1 2 s 1 2 + 1 2 e x 3 2 + 1 2 e x 4 2
Differentiating Equation (24), yields
V e s = s 1 s ˙ 1 = s 1 c 1 x 2 + c 2 ς / μ x 1 ς μ / μ x 2 + c 3 x ˙ 2 = s 1 c 1 x 2 + c 2 ς / μ x 1 ς μ / μ x 2 + c 3 b 2 x 3 = s 1 c 1 x 2 + c 2 ς / μ x 1 ς μ / μ x 2 + c 3 b 2 e x 3 + x 3 d = s 1 c 1 x 2 + c 2 ς / μ x 1 ς μ / μ x 2 + c 3 b 2 e x 3 + c 3 b 2 x 3 d = s 1 c 3 b 2 e x 3 k l s 1 k s s 1 λ / η k l s 1 2 + c 3 b 2 s 1 2 + e x 3 2 / 4
Similarly,
V ˙ e x 3 = e x 3 e ˙ x 3 = e x 3 x ˙ 3 x ˙ 3 d = e x 3 f 3 x 3 + b 3 x 4 x ˙ 3 d = e x 3 f 3 x 3 + b 3 e x 4 + x 4 d x ˙ 3 d = e x 3 b 3 e x 4 + k 3 x 3 d x 3 = e x 3 b 3 e x 4 k 3 e x 3 k 3 e x 3 2 + b 3 e x 3 2 4 + e x 4 2
V ˙ e x 4 = e x 4 e ˙ x 4 = e x 4 x ˙ 4 x ˙ 4 d = e x 4 f 4 x 4 + b 4 x 5 x ˙ 4 d = e x 4 f 4 x 4 + b 4 x 5 d x ˙ 4 d = e x 4 k 4 x 4 d x 4 = k 4 e x 4 2
According to Equations (24)–(27), V ˙ e can be obtained as follows.
V ˙ e = V ˙ e s + V ˙ e x 3 + V ˙ e x 4 = s 1 s ˙ 1 + e x 3 e ˙ x 3 + e x 4 e ˙ x 4 k l s 1 2 + c 3 b 2 s 1 2 + e x 3 2 / 4 + k 3 e x 3 2 + b 3 e x 3 2 4 + e x 4 2 + k 4 e x 4 2 = k l c 3 b 2 s 1 2 k 3 c 3 b 2 4 b 3 4 e x 3 2 k 4 b 4 b 3 e x 4 2 = e s , e x 3 , e x 4 k l c 3 b 2 k 3 c 3 b 2 4 b 3 4 k 4 b 4 b 3 e s , e x 3 , e x 4 T
where Q is a positive definite matrix. By adjusting the control parameters, we can get the following equation.
V ˙ e = V ˙ e s + V ˙ e x 3 + V ˙ e x 4 = s 1 s ˙ 1 + e x 3 e ˙ x 3 + e x 4 e ˙ x 4 e 2 Q
Hence, the closed-loop system is asymptotically stable.

5. Simulation

In this section, this paper will give several examples to illustrate the effectiveness of the algorithm. The aerodynamic parameters for the air vehicle are defined as follows:
p 1 = 3.1166 , p 2 = 0.2337 , p 3 = 82.6918 p 4 = 0.9749 , p 5 = 128.6316
Considering the real properties of the air vehicle, the control constraint is set as δ z c max 20 . All these simulation step sizes are 0.01 s.

5.1. Case I

In this section, we will demonstrate the autonomous and collaborative capability of the air vehicle. So, the paper will use four air vehicles to intercept two objects. The initial parameters if the air vehicle and the object are shown in Table 1 and Table 2.
The simulation results are shown in Figure 4, Figure 5, Figure 6 and Figure 7.
Figure 4 shows the pursuit trajectory with the two objects and four air vehicles. It is obvious that all the objects can be intercepted. In addition, the 2 t h air vehicle and 3 t h air vehicle do not approach the object. This is due to the 1 t h air vehicle and 4 t h air vehicle having intercepted the object. Further analysis reveals that the 1 t h air vehicle and 2 t h air vehicle have an initial air vehicle–object distance advantage over the 2 t h object, while there is at least one intersection point between the 1 t h air vehicle, the 2 t h air vehicle, and the 2 t h object in the direction of velocity. Thus, the 1 t h air vehicle and the 2 t h air vehicle have no cooperative interception relationship with the 2 t h object. Instead, there is at least one intersection point between the 3 t h air vehicle, the 4 t h air vehicle, and the 2 t h object in the direction of velocity. Additionally, they did not have the distance advantage over the 1 t h object. Thus, the 3 t h air vehicle and the 4 t h air vehicle have a cooperative interception relationship with the 2 t h object.
As can be seen in Figure 5, which shows the rudder deflection angle for each air vehicle, it is clear that the changing of the rudder deflection angle is relatively small during the whole interception process, which also reflects the controller designed in this paper as having a better control performance.
As shown in Figure 6, which gives the pitch angle variation curve for each air vehicle during the whole interception process, it is clear that the variation of each curve is relatively smooth, which indirectly reflects the air vehicle attitude changing stably.
As can been seen in Figure 7, which presents the variation of the air vehicle–object distance for each air vehicle, it is obvious that the distance of the 2 t h air vehicle and the 3 t h air vehicle are not converged at 0. The reason for this phenomenon is that the 2 t h air vehicle and the 1 t h air vehicle are in a cooperative state, and the 1 t h object is intercepted by the 1 t h air vehicle. Similarly, the 3 t h air vehicle and the 4 t h air vehicle are in a cooperative state, and the 2 t h object is intercepted by the 4 t h air vehicle.

5.2. Case II

To further verify the autonomous cross-interception capability of the method, this paper will take the interception process of vertical plane as an example. In this case, the objects do not appear at the same time as the 2 t h object appears, which is three seconds later than the 1 t h object. The initial parameters of the air vehicles and the objects are shown in Table 3 and Table 4. The cooperative interception relationship between the air vehicle and the objects is shown in Table 5. Additionally, all these simulations are shown in Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13.
As can be seen in Figure 8 and Figure 9, which show the pursuit trajectory with two objects and four air vehicles at different times, it is clear that all the objects can be intercepted. In addition, the 2 t h air vehicle and 3 t h air vehicle do not approach the object. This is due to the 1 t h air vehicle and 4 t h air vehicle having intercepted the object. Furthermore, based on the dynamic situation assessment model, all these air vehicles are in a collaborative relationship to intercept the 1 t h object before 3 s, while with the 2 t h object appears and the motion of these air vehicles, m 3 , m 4 , does not satisfy the conditions of cooperative interception of the 1 t h object. Instead, these two air vehicles are in a collaborative relationship in the interception of the 2 t h object.
As shown in Figure 10, which presents the changing of the air vehicle–object distance for each air vehicle, it is obvious that the distance of the 2 t h air vehicle and the 3 t h air vehicle are not converged to 0. The reason for this phenomenon is that the 2 t h air vehicle and the 1 t h air vehicle are in a cooperative state and the 1 t h object is intercepted by the 1 t h air vehicle. Similarly, the 3 t h air vehicle and the 4 t h air vehicle are in a cooperative state and the 2 t h object is intercepted by the 4 t h air vehicle. In addition, it can also be concluded that the air vehicle–object distance of the 3 t h air vehicle and the 4 t h air vehicle changes suddenly at 3 s. This is due to the air vehicle changing its interception object.
As shown in Figure 11 and Figure 12, which show the variation of LOS angle and LOS rate for each air vehicle, it is clear that the LOS angle and LOS rate for the 3 t h air vehicle and the 4 t h air vehicle changes suddenly at 3 s. This is also because these two air vehicles change the interception object from the 1 t h object to the 2 t h object.
As shown in Figure 13 and Figure 14, which give the variation of the radial relative velocity and the normal relative velocity for each air vehicle, it is clear that the radial relative velocity for the 3 t h air vehicle and the 4 t h air vehicle changes suddenly at 3 s. This also demonstrates that these two air vehicles change their interception object from the 1 t h object to the 2 t h object.
As shown in Table 5, which shows the cooperation state between the objects and air vehicles at different times, it is clear that before the 2 t h object appears, all air vehicles are in a collaborative state to intercept the 1 t h object. After 3 s, m 1 , m 2 are in a collaborative relationship between 1 t h object, and m 3 , m 4 are in a collaborative relationship between 2 t h object.
Thus, according to all these simulations, we can conclude that the proposed method in this paper can provide a better autonomous cross-cooperation interception ability during multiple-air vehicle and multiple-object interception. This also verifies the statement in Definition 1.

6. Conclusions

To solve the problem of multi-object interception with multiple air vehicles, a new dynamic cross-collaborative interception algorithm based on GTSMC and virtual geometry is discussed in this paper. Firstly, the paper formulates the problem of multi-air vehicle to multi-object interception and gives the definition of cross-collaborative interception. Instead of considering the particle model, the paper takes the dynamics of the air vehicle into account. To evaluate the situation in real time, the paper designs the dynamic situation assessment model, which includes the flight statues and cooperation statues of the air vehicle. Further, the paper defines the condition of the air vehicle, whether the air vehicle is in a cooperative state, and proves it. After introducing the dynamic situation assessment model, the paper designs a new controller by combining GTSMC with the idea of backstepping method. Then, the paper proves the stability of system by using Lyapunov theory. Finally, several simulation cases which consider different interception scenarios are given to illustrate the effectiveness of the method. The simulation results show that the new collaborative interception algorithm can provide a better autonomous cross-collaborative interception capability. Even so, there are two important issues that need to be addressed. The first is how to avoid a collision between the air vehicles, and the second is establishing a dynamic situation assessment model with an overload constraint. These two schemes are actively important in future research.

Author Contributions

Conceptualization, K.N., X.B., X.C. and J.Y.; methodology, K.N., X.B., X.C., J.Y. and H.L.; formal analysis, K.N. and X.C.; data curation, H.L.; writing—original draft preparation, K.N., X.B., X.C. and J.Y.; writing—review and editing, K.N., X.B., X.C., J.Y. and H.L. All authors have read and agreed to the published version of the manuscript.

Funding

The authors are grateful for the support of the National Natural Science Foundation of China (No. 61673212).

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the flight control lab of Beijing University of Technology for their support for this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shi, L.; Zhu, Z.H.; Tang, S.; Yan, X.D. Hybrid cooperative guidance law for active aircraft defense against a guided missile. J. Guid. Control Dyn. 2018, 41, 531–537. [Google Scholar]
  2. Zhang, Y.A.; Wang, X.L.; Wu, H.L. Impact time control guidance law with field of view constraint. Aerosp. Sci. Technol. 2014, 39, 361–369. [Google Scholar] [CrossRef]
  3. Shalumov, V.; Shima, T. Weapon-target-allocation strategies in multiagent target-missile-defender engagement. J. Guid. Control Dyn. 2017, 40, 2452–2464. [Google Scholar] [CrossRef]
  4. Tekin, R.; Erer, K.S. Impact time and angle control against moving targets with look angle shaping. J. Guid. Control Dyn. 2020, 43, 1020–1025. [Google Scholar] [CrossRef]
  5. Ou, A.; Ziqian, Z. A method of threat assessment based on MADM and results of situation assessment in air-to-air combat. Fire Control Radar Technol. 2006, 35, 64–67. [Google Scholar]
  6. Kong, S.; Zhang, H. Target Track Estimation and Threat Assessment Based on Multi-Source Information Fusion. Ph.D. Thesis, China Aerospace Science and Technology Corporation, Beijing, China, 2017. [Google Scholar]
  7. Tahk, M.; Shim, S.; Hong, S.; Choi, H.; Lee, C. Impact time control based on time-to-go prediction for sea-skimming antiship missiles. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 2043–2052. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Tang, S.; Guo, J. Two-stage cooperative guidance strategy using a prescribed-time optimal consensus method. Aerosp. Sci. Technol. 2020, 100, 105641. [Google Scholar] [CrossRef]
  9. Kumar, S.R.; Mukherjee, D. Cooperative salvo guidance using finite-time consensus over directed cycles. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 1504–1514. [Google Scholar] [CrossRef]
  10. Zhao, J.; Zhou, R. Distributed three-dimensional cooperative guidance via receding horizon control. Chin. J. Aeronaut. 2016, 29, 972–983. [Google Scholar] [CrossRef]
  11. Qilun, Z.; Xiwang, D.; Liang, Z.; Chen, B.; Jian, C.; Zhang, R. Distributed cooperative guidance for multiple missiles with fixed and switching communication topologies. Chin. J. Aeronaut. 2017, 30, 1570–1581. [Google Scholar]
  12. Wei, X.; Yang, J. Cooperative guidance laws for simultaneous attack against a target with unknown maneuverability. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2019, 233, 2518–2535. [Google Scholar] [CrossRef]
  13. Shiyu, Z.; Rui, Z. Cooperative guidance for multi-missile salvo attack. Chin. J. Aeronaut. 2008, 21, 533–539. [Google Scholar] [CrossRef]
  14. Wang, X.H.; Tan, C.P. Distributed cooperative controller design considering guidance loop and impact angle. J. Frankl. Inst. 2018, 355, 6927–6946. [Google Scholar] [CrossRef]
  15. Li, G.F.; Wu, Y.J.; Xu, P.Y. Adaptive fault-tolerant cooperative guidance law for simultaneous arrival. Aerosp. Sci. Technol. 2018, 82, 243–251. [Google Scholar] [CrossRef]
  16. Yamasaki, T.; Balakrishnan, S.N.; Takano, H. Integrated guidance and autopilot design for a chasing UAV via high-order sliding modes. J. Frankl. Inst. 2012, 349, 531–558. [Google Scholar] [CrossRef]
  17. Guo, J.G.; Wang, X.M.; Zhou, J. Efficient information-based cooperative guidance law of multi-missiles. Trans. Inst. Meas. Control 2019, 41, 2651–2665. [Google Scholar]
  18. He, S.M.; Kim, M.; Song, T.; Lin, D.F. Three-dimensional salvo attack guidance considering communication delay. Aerosp. Sci. Technol. 2018, 73, 1–9. [Google Scholar] [CrossRef]
  19. Wei, X.; Yang, J. Finite time simultaneous attack for a maneuvering target with unknown acceleration. Trans. Inst. Meas. Control 2019, 41, 1849–1860. [Google Scholar] [CrossRef]
  20. He, S.; Lin, D. Three-dimensional optimal impact time guidance for antiship missiles. J. Guid. Control Dyn. 2019, 42, 941–948. [Google Scholar] [CrossRef]
  21. Hu, Q.; Han, T.; Xin, M. Sliding-mode impact time guidance law design for various target motions. J. Guid. Control Dyn. 2019, 42, 136–148. [Google Scholar] [CrossRef]
  22. Hou, Z.; Yang, Y.; Liu, L.; Wang, Y. Terminal sliding mode control-based impact time and angle constrained guidance. Aerosp. Sci. Technol. 2019, 93, 105142. [Google Scholar] [CrossRef]
  23. Mukherjee, D.; Kumar, S.R. Field-of-View Constrained Impact Time Guidance Against Stationary Targets. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 3296–3306. [Google Scholar] [CrossRef]
  24. Li, X.; Chen, W.; Chen, Z.; Wang, T.; Shi, H. Fixed-Time Circular Impact-Time Guidance with Look Angle Constraint. Aerospace 2022, 9, 356. [Google Scholar] [CrossRef]
  25. You, H.; Chang, X.L.; Zhao, J.F.; Wang, S.H.; Zhang, Y.H. Three-dimensional impact-angle-constrained cooperative guidance strategy against maneuvering target. ISA Trans. 2023, in press. [Google Scholar] [CrossRef] [PubMed]
  26. Chun, W.; Dong, W.; Jia, W. Impact-angle-constrained cooperative guidance for salvo attack. J. Guid. Control Dyn. 2022, 45, 684–703. [Google Scholar]
  27. He, S.; Wang, W.; Lin, D.; Lei, H. Consensus-based two-stage salvo attack guidance. IEEE Trans. Aerosp. Electron. Syst. 2017, 54, 1555–1566. [Google Scholar] [CrossRef]
  28. Li, B.; Lin, D.F.; Wang, H. Finite time convergence cooperative guidance law based on graph theory. Optik 2016, 127, 10180–10188. [Google Scholar] [CrossRef]
  29. Zhao, J.; Zhou, R.; Dong, Z. Three-dimensional cooperative guidance algorithms against stationary and maneuvering targets. Chin. J. Aeronaut. 2015, 28, 1104–1120. [Google Scholar] [CrossRef]
  30. Cheng, Z.; Liu, L.; Wang, Y. Lyapunov-based switched-gain impact angle control guidance. Chin. J. Aeronaut. 2018, 31, 765–775. [Google Scholar] [CrossRef]
  31. Lin, L.-G.; Xin, M. Missile guidance law based on new analysis and design of SDRE scheme. J. Guid. Control Dyn. 2019, 42, 853–868. [Google Scholar] [CrossRef]
  32. Kumar, S.R.; Maity, A. Finite-horizon robust suboptimal control-based impact angle guidance. IEEE Trans. Aerosp. Electron. Syst. 2019, 56, 1955–1965. [Google Scholar] [CrossRef]
  33. He, S.; Song, T.; Lin, D. Impact angle constrained integrated guidance and control for maneuvering target interception. J. Guid. Control Dyn. 2017, 40, 2652–2660. [Google Scholar] [CrossRef]
Figure 1. 2D guidance geometry of multiple-air vehicle and multiple-object.
Figure 1. 2D guidance geometry of multiple-air vehicle and multiple-object.
Aerospace 10 00728 g001
Figure 2. Different types of the intercepting triangle under two air vehicles and one object.
Figure 2. Different types of the intercepting triangle under two air vehicles and one object.
Aerospace 10 00728 g002
Figure 3. The rotation of the interception triangle.
Figure 3. The rotation of the interception triangle.
Aerospace 10 00728 g003
Figure 4. The pursuit trajectory with two objects and four air vehicles.
Figure 4. The pursuit trajectory with two objects and four air vehicles.
Aerospace 10 00728 g004
Figure 5. The rudder deflection angle of each air vehicle.
Figure 5. The rudder deflection angle of each air vehicle.
Aerospace 10 00728 g005
Figure 6. The pitch angle of each air vehicle.
Figure 6. The pitch angle of each air vehicle.
Aerospace 10 00728 g006
Figure 7. The distance curve of each air vehicle.
Figure 7. The distance curve of each air vehicle.
Aerospace 10 00728 g007
Figure 8. The pursuit trajectory with two objects and four air vehicles before 3 s.
Figure 8. The pursuit trajectory with two objects and four air vehicles before 3 s.
Aerospace 10 00728 g008
Figure 9. The pursuit trajectory with two objects and four air vehicles.
Figure 9. The pursuit trajectory with two objects and four air vehicles.
Aerospace 10 00728 g009
Figure 10. The relative distance curve for each air vehicle.
Figure 10. The relative distance curve for each air vehicle.
Aerospace 10 00728 g010
Figure 11. The LOS angle for each air vehicle.
Figure 11. The LOS angle for each air vehicle.
Aerospace 10 00728 g011
Figure 12. The radial relative speed angle for each air vehicle.
Figure 12. The radial relative speed angle for each air vehicle.
Aerospace 10 00728 g012
Figure 13. The radial relative velocity for each air vehicle.
Figure 13. The radial relative velocity for each air vehicle.
Aerospace 10 00728 g013
Figure 14. The normal relative velocity for each air vehicle.
Figure 14. The normal relative velocity for each air vehicle.
Aerospace 10 00728 g014
Table 1. Parameters of the air vehicles.
Table 1. Parameters of the air vehicles.
Air VehicleM1M2M3M4
(xtj, ytj)/km(0, 2.5)(0, 2)(0, 0)(2, 0)
θ /deg70703070
Table 2. Parameters of the objects.
Table 2. Parameters of the objects.
ObjectT1T2
(xtj, ytj)/km(6.9, 4)(6, 3.5)
θ /deg175−135
Table 3. Parameters of the air vehicles.
Table 3. Parameters of the air vehicles.
Air VehicleM1M2M3M4
(xtj, ytj)/km(0.5, 0)(1, 0)(1.7, 0)(2.5, 0)
θ /deg90909090
Table 4. Parameters of the objects.
Table 4. Parameters of the objects.
ObjectT1T2
(xtj, ytj)/km(2.1, 5)(0.6, 5)
θ /deg−135−315
Table 5. The cooperative interception relationship.
Table 5. The cooperative interception relationship.
T1T2
before 3 sM1, M2, M3, M4-
after 3 sM1, M2M3, M4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Niu, K.; Bai, X.; Chen, X.; Yu, J.; Liu, H. A Dynamic Cross-Collaborative Interception Algorithm Based on GTSMC and Virtual Geometry. Aerospace 2023, 10, 728. https://doi.org/10.3390/aerospace10080728

AMA Style

Niu K, Bai X, Chen X, Yu J, Liu H. A Dynamic Cross-Collaborative Interception Algorithm Based on GTSMC and Virtual Geometry. Aerospace. 2023; 10(8):728. https://doi.org/10.3390/aerospace10080728

Chicago/Turabian Style

Niu, Kang, Xu Bai, Xi Chen, Jianqiao Yu, and Haiying Liu. 2023. "A Dynamic Cross-Collaborative Interception Algorithm Based on GTSMC and Virtual Geometry" Aerospace 10, no. 8: 728. https://doi.org/10.3390/aerospace10080728

APA Style

Niu, K., Bai, X., Chen, X., Yu, J., & Liu, H. (2023). A Dynamic Cross-Collaborative Interception Algorithm Based on GTSMC and Virtual Geometry. Aerospace, 10(8), 728. https://doi.org/10.3390/aerospace10080728

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