Next Article in Journal
A Spatio-Temporal Spotting Network with Sliding Windows for Micro-Expression Detection
Previous Article in Journal
Performance Evaluation of Single-Carrier and Orthogonal Frequency Divison Multiplexing-Based Autoencoders in Comparison with Low-Density Parity-Check Encoder
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Fuzzy Inference System-Based Three-Dimensional Resolution Algorithm for Collision Avoidance of Fixed-Wing UAVs

Department of Aerospace Engineering & Engineering Mechanics, University of Cincinnati, Cincinnati, OH 45221, USA
*
Author to whom correspondence should be addressed.
Electronics 2023, 12(18), 3946; https://doi.org/10.3390/electronics12183946
Submission received: 23 August 2023 / Revised: 8 September 2023 / Accepted: 13 September 2023 / Published: 19 September 2023
(This article belongs to the Section Systems & Control Engineering)

Abstract

:
Fixed-wing Unmanned Aerial Vehicles (UAVs) cannot fly at speeds lower than critical stall speeds. As a result, hovering during a potential collision scenario, like with rotary-wing UAVs, is impossible. Moreover, hovering is not an optimal solution for Collision Avoidance (CA), as it increases mission time and is innately fuel-inefficient. This work proposes a decentralized Fuzzy Inference System (FIS)-based resolution algorithm that modulates the point-to-point mission path while ensuring the continuous motion of UAVs during CA. A simplified kinematic guidance model with coordinated turn conditions is considered to control the UAVs. The model employs a proportional-derivative control of commanded airspeed, bank angle, and flight path angle. The commands are derived from the desired path, characterized by airspeed, heading, and altitude. The desired path is, in turn, obtained using look-ahead points generated for the target point. The FIS aims to mimic human behavior during collision scenarios, generating modulation parameters for the desired path to achieve CA. Notably, it is also scalable, which makes it easy to adjust the algorithm parameters, as per the required missions, and factors specific to a given UAV. A genetic algorithm was used to optimize FIS parameters so that the distance traveled during the mission was minimized despite path modulation. The proposed algorithm was optimized using a pairwise conflict scenario. The effectiveness of the algorithm was evaluated through a Monte Carlo simulation of random conflict scenarios involving multiple UAVs operating in a confined space.

1. Introduction

Global air traffic is expected to grow annually at a rate of 4.1% [1], and global advanced air mobility of automated aircraft at lower altitudes is expected to grow at a compounded annual growth rate of 22.45% [2]. Unmanned Aerial Vehicles (UAVs) would be an integral part of this future air traffic, used for surveillance, military, transportation, cargo, and many more applications. Thus, given the inevitability of potential collisions among aircraft, Collision Avoidance (CA) algorithms are vital for all UAVs, including fixed-wing UAVs. Fixed-wing UAVs carry many advantages over rotary-wing UAVs like greater speed and endurance, higher payload capacity for the same endurance, lower maintenance, lower noise levels, etc., which could be essential for many emergency, reconnaissance, and cargo missions. However, the formulation of CA algorithms is challenging for fixed-wing UAVs. Unlike rotary-wing UAVs, fixed-wing UAVs cannot hover, or even slow down below critical stall speeds during potential collision scenarios. Hovering is also a non-optimal solution for CA as it increases mission time and is innately fuel-inefficient. Thus, this work proposes a decentralized Fuzzy Inference System (FIS)-based resolution algorithm that modulates the point-to-point mission path while ensuring the continuous motion of UAVs during CA. This algorithm, ensuring continuous motion, could likewise provide an efficient solution for CA in rotary-wing UAVs, eliminating the need for hover-based solutions in numerous potential collision scenarios.
Aircraft CA System (ACAS) II is the most common method of CA in crewed aircraft. It is implemented as Traffic CAS (TCAS) II, which provides traffic advisory for warning and resolution advisory for appropriate CA maneuvers to the pilots [3]. However, TCAS II has several limitations, encompassing the provision of resolutions exclusively within the vertical plane like climb, descent, or level flight, and these resolutions are designed for specific encounter geometries; furthermore, TCAS II operates deterministically [4]. The review and analysis of TCAS in [5] also state that expanding the horizontal resolution algorithm will offer superior collision detection and resolution performance. Similarly, the newly proposed method of ACAS X and its variation ACAS Xu for unmanned aircraft encourage the design of resolutions in both vertical and horizontal directions using probabilistic models [6,7]. Therefore, this study proposes to generate resolutions in both vertical and horizontal directions.
A vast number of research studies and algorithms have been proposed for the CA of aircraft, widely classified by Kuchar et al. [8] as prescribed, optimized, force field, or manual. With the advancement of navigation and control technologies for UAVs, CA algorithms can be further classified into various domains, including path planning, conflict resolution, potential function, geometric guidance, and motion planning [9,10]. However, it is essential to note that many of these approaches primarily cater to rotary-wing UAVs capable of stopping and hovering, rendering them less suitable for fixed-wing UAVs. For fixed-wing aircraft, a range of geometric approaches, such as forming circular arcs or Three-Dimensional (3D) trajectories [11,12,13], velocity modulation in the horizontal plane [14], velocity obstacle concepts [15,16,17,18], differential geometry concepts [19], and the use of collision cones [20,21], have been proposed. While these methods provide some applicability to fixed-wing aircraft, they may suffer from limitations such as high computational demands and a lack of intuitive understanding. Furthermore, some concepts involve potential/force fields modified or reformulated with UAV’s physical constraints for better real-life performance [22,23,24], or enhanced and optimized using evolutionary methods to avoid local minima problems [22,25,26]. Unfortunately, many force field methods, like the optimized path planning-based approach of the 3D vector field histogram [27] are predominantly suited for rotary-wing UAVs and may not translate well to fixed-wing UAV scenarios. While path planning-based methods offer high accuracy, some of them, like combining differential game problems with tree-based path planning [28], utilizing reinforcement learning for UAV guidance [29], optimizing flight trajectory with a bank-turn mechanism [30], employing the Radau-pseudospectral approach [31], applying probabilistic methods in collision detection [32], or generating an automated distributed policy for multi-robot motion planning [33], may face significant computational challenges when implemented in real time for small UAVs. Analytical formulations such as the speed approach [34] and the use of buffered Voronoi Cells for path planning [35], though computationally efficient, may lack intuitive understanding and require specialized expertise to implement effectively.
Humans, on the other hand, are innately capable of avoiding collision intuitively [36,37]. So, in recent times, fuzzy logic that emulates the human way of decision making with the linguistic characterization of numerical variables has been gaining popularity in devising CA algorithms in UAVs. This research also uses an FIS to avoid computational complexity while also providing an intuitive understanding of the approach. Just like pilots who perform CA maneuvers manually based on their visual perception and flight instrument information about incoming static or dynamic obstacles, the proposed algorithm also takes into account the sensor readings for incoming intruders to modulate the current path and avoid collisions. Previously, fuzzy logic for CA was predominantly researched in surface vehicles [38,39,40,41]. Now, research on the use of fuzzy logic for CA in UAVs is also gaining traction. For instance, a fuzzy-based aircraft CA system capable of generating an alert of a potential midair collision while taking control if no preventive action is taken within a specified time was proposed by Younas et al. [42]. Choi et al. [43] used an FIS to improve the performance of their enhanced potential field-based CA. Likewise, several basic avoidance methods have been devised to avoid collision using fuzzy logics [44,45,46]. However, most of them are suitable for rotary-wing UAVs. Cook et al. [47] used a fuzzy logic-based approach to help mitigate the risk of collisions among aircraft, including fixed-wing and quadcopter, using separation assurance and CA techniques.
Although similar to the CA logic of [47], which directly generates low-level control input like turn rate using intuitively defined FISs, this work proposes a 3D Collision Avoidance Fuzzy Inference System (3D-CAFIS) to modulate the path variables—heading, velocity, and altitude of the UAVs. Thus, the UAV performs both horizontal and vertical maneuvers. The modulation is based on the relative states of the UAV in consideration, the own-ship, with respect to the nearest intruder or obstacle point. Moreover, a Genetic Algorithm (GA) was used to optimize the FISs such that the distance traveled during the mission is minimal despite path modulation while ensuring separation during CA is above the appropriate threshold. FIS optimization was conducted using a training scenario of three basic pairwise conflicts. The algorithm’s validity can be confirmed using formal methods, which will be similar to that outlined in [47], as the basic approach of defining the FISs for this study was similar. However, due to the complexity of the algorithm with many variables, the effectiveness of the algorithm is evaluated through multiple simulations of randomly generated UAV missions within a closed functional space.
The paper begins by introducing the system architecture and the UAV model utilized in this study in Section 2. Next, Section 3 provides a step-by-step explanation of the proposed algorithm. Section 4 presents the findings from our simulation studies. Lastly, Section 5 offers a summary of this work and addresses future works.

2. System Modeling

UAVs operate in the N-frame. Hereafter, ‘own-ship’ refers to the UAV currently executing the CA algorithm, while ‘intruders’ denote other UAVs in the airspace. The own-ship’s body frame is designated as the O-frame and the intruder’s as the I-frame. A UAV’s position in the N-frame is defined by N ρ , and its attitude is represented by the 3-2-1 set of Euler angles ( Λ = ψ , θ , ϕ ), as shown in Figure 1. For simplicity, one omits the N-frame expression.
Furthermore, the state of a fixed-wing UAV can be comprehensively defined with parameters including airspeed ( ν air ), ground speed ( ν ground ), wind velocity ( ν wind ), α , β , χ , and γ , as shown in Figure 2.
The UAV’s attitude can also be represented in the G-frame, obtained using the 3-2-1 set of Euler angles ( G Λ = χ , γ , ϕ ). Consequently, the UAV state space vector considered is p = ρ , ν air , G Λ , ϕ ˙ . Note that ν air is the norm of ν air .
A basic system architecture for the waypoint control and CA of a UAV is depicted in Figure 3.
In summary, the significance of each block is explained below:
  • Waypoint Follower (Section 3.1) determines the desired path along with the desired change variables based on set waypoints.
  • Collision Detection (CD) (Section 3.2) checks whether the own-ship is on a collision course with an intruder using range-bearing sensors and modified TCAS logic.
  • 3D-CAFIS (Section 3.3) generates modulation parameters based on relative states obtained from CD using the proposed FIS tree. These parameters modulate the desired path obtained from the waypoint follower, thus avoiding potential collisions.
  • Controller (Section 3.4) calculates control commands based on the path, either from the waypoint follower or, if CA is initiated, from the 3D-CAFIS, to control the UAV.
  • UAV Plant represents the kinematic mathematical model of a fixed-wing UAV that simulates UAV motion.
  • Sensors consist of a global positioning system and an inertial measurement unit for determining the position, velocity, and angular velocity of the own-ship, along with a range-bearing sensor that measures the relative position of intruders in airspace.
  • State Estimator utilizes filtering techniques to estimate the position and attitude of the UAV. This is not within the scope of this research.
Here, UAVs are controlled via a proportional-derivative control of ν air , γ , and ϕ . The subsequent kinematic guidance model for the UAV is given by [48]
ρ ˙ = C N / G G ν ground ,
ν ˙ air = k ν air ν air c ν air ,
γ ˙ = k γ γ c γ ,
χ ˙ = g cos ( χ ψ ) ν ground tan ϕ ,
ϕ ¨ = k ϕ ϕ c ϕ + k ϕ ˙ ϕ ˙ ,
where G ν ground = ν ground , 0 , 0 is the ground speed vector in the G-frame, while k γ , k ν air , k ϕ , and k ϕ ˙ are gain values tuned to achieve smoother flight maneuvers. Note that under wind-less conditions, ν air = ν ground , ψ = χ , and θ = γ . For simplicity, UAVs are assumed to be in a coordinated-turn condition with zero side-slip and zero angles of attack. The relationship between the course angle and bank angle is expressed in Equation (4). Also, ϕ c , γ c , and ν air c serve as the command variables for controlling the UAV.

3. Proposed Method

3.1. Waypoint Follower

The waypoint follower is a path generator that calculates the desired variables, fully defining a trajectory for the UAV during a mission to a waypoint. This approach draws inspiration from pure pursuit control [28], incorporating the look-ahead point position vector ( ρ l ) and l into the calculation of the desired variables. Given a waypoint position vector ( ρ w ) in the N-frame, its relative position vector with respect to the own-ship is expressed as follows:
ρ w / O = ρ w ρ O .
Then, the look-ahead point position vector is obtained via
ρ l = ρ + l ρ ^ w / O ,
and the desired path requisites are calculated as follows:
χ d = tan 1 ρ 2 l ρ 2 ρ 1 l ρ 1 ,
Δ χ d = χ d χ ,
h d = ρ 3 l ,
Δ h d = h d h .
The desired airspeed ( ν air d ) is determined based on the mission requirements. In this research, it is assumed to be equal to an optimal airspeed ( ν air opt ), which gradually decreases quadratically as the target point is approached. Thus, the desired path is characterized by ν air d , χ d , Δ χ d , h d , and Δ h d .

3.2. Collision Detection (CD)

The CD begins with intruder detection, where the intruder’s presence is initially identified using the intruder’s 3D spherical coordinates ( r , η , ϵ ) acquired from the onboard range-bearing sensor. These coordinates are then transformed into Cartesian form to establish the separation between the own-ship and the intruder ( O s = O ρ I / O ). The horizontal and vertical views of two potentially colliding UAVs—one own-ship in the O-frame and another intruder in the I-frame—are depicted in Figure 4. Figure 4a shows the top-down view of the two UAVs, while Figure 4b offers a side view of the UAVs on a vertical plane. Here, ν O and ν I are the velocities of the own-ship and intruder, respectively, and their relative velocity is expressed as ν O / I . Also, ψ M is the angle formed by s with respect to n ^ 1 , while χ O / I is the relative course angle of the own-ship with respect to the intruder.
The relative velocity and separation between the own-ship and the intruder play a crucial role in CD. These values are essential for applying the CD logic provided by TCAS II, which assesses whether the own-ship is on a potential collision course with an intruder. To begin, the relative velocity of the intruder with respect to the own-ship, as observed in the N-frame, can be determined using the transport theorem, and this is expressed in the O-frame as follows:
N O s ˙ = O O s ˙ + O ω O / N × O s .
Here, O O s ˙ is determined with successive range-bearing measurements as follows:
O O s ˙ = O s t O s t Δ t Δ t ,
where Δ t is the time interval between two successive measurements. Also, the angular velocity in the O-frame is derived from the inertial measurement unit.
Now, one introduces an M-frame achieved by rotating the N-frame about n ^ 3 by ψ M . The horizontal and vertical separations (x and z) and their respective closure rates ( ν x and ν z ) can be determined from the relative separation and relative velocity when expressed in the M-frame, respectively, as shown below:
M s = C M / N C N / O O s = x , y , z ,
N M ν O / I = C M / N C N / O N O s ˙ = ν x , ν y , ν z .
Note that z = h , and thus, z d = h d and Δ z d = Δ h d .
The CD is confirmed if the time to the Closest Point of Approach (CPA) is less than the CPA thresholds ( τ th ) in both the horizontal range ( τ x ) and vertical range ( τ z ) directions. These calculations are carried out as described below [49]:
τ x = ν x / x ,
τ z = ν z / z .
The thresholds are selected using TCAS II logic, which considers the altitude of operation and the sensitivity level [3]. However, since TCAS is designed for large aircraft, and this study focuses on small UAVs, the thresholds are determined based on the specific UAV and sensor limitations. Similarly, χ O / I and ψ M are crucial relative variables for CA. They assist in determining inputs, namely Δ ψ M and sign Δ χ O / I , for the 3D-CAFIS. In this context, Δ ψ M indicates the anticipated position of the own-ship with respect to the intruder at the CPA, while sign Δ χ O / I denotes the direction of relative motion in comparison to the current course angle, thereby distinguishing between collisions in the same direction and opposite directions. These values are calculated as follows:
Δ ψ M = χ O / I ψ M ,
sign ( Δ χ O / I ) = sign ( χ O / I χ O ) .
However, in cases of slower closure rates, this method could trigger CD alerts even when a significant separation exists. Therefore, when dealing with slower closure rates, a comparison is made between the horizontal and vertical separations against the distance modification threshold (DMOD) and altitude threshold (ZTHR) for the small UAVs, respectively. A collision is assumed if either the horizontal separation is below the threshold ( x th ) or the vertical separation is below the threshold ( z th ).

3.3. Three-Dimensional Collision Avoidance Fuzzy Inference System (3D-CAFIS)

3.3.1. Fuzzy Tree

A fuzzy tree is designed to obtain path modulation parameters from the relative variables determined in CD, as shown in Figure 5. The blocks FIS_H1 and FIS_H2 represent FISs that yield modulation parameters Δ χ ca and Δ ν x ca , respectively, for horizontal maneuvers. Similarly, the blocks FIS_V1 and FIS_V2 are FISs that yield modulation parameters Δ h ca and Δ ν z ca , respectively, for vertical maneuvers. Ultimately, the outputs Δ χ ca and Δ h ca serve as inputs to the FIS_W block, yielding Υ x and Υ z which are weights for horizontal and vertical maneuvers, respectively. Note that there is an input to the fuzzy tree, Δ z d * , which is a modified version of Δ z d designed to induce complementary altitude adjustments between potentially colliding UAVs coming from opposite directions. Its calculation is as follows:
Δ z d * = sign ( Δ χ O / I ) Δ z d
The inputs of each FIS block are normalized using a factor determined according to the potential range of the actual input values, which takes into account factors such as sensor limits, UAV performance limits, and so on. The outputs of the fuzzy tree are derived by multiplying the normalized outputs of the FISs by their respective normalization factors. These factors are determined based on the expected range that the respective outputs can span over the sampling time. The relative variables utilized as inputs for the FIS tree, as well as the modulation parameters for outputs, are presented in Table 1, along with their corresponding normalization factor and significance.
For this study, each of the normalized inputs and outputs of the FISs is fuzzified into three Membership Functions (MFs)—two trapezoidal ones at each end and one triangular MF in the middle, as shown in Figure 6. The decision to use three MFs serves a dual purpose. Firstly, it provides simplicity in emulating three fundamental directions or maneuvers associated with human behavior during CA, such as left, center, and right or up, in front, and down. Secondly, the MFs for all inputs and outputs can be symmetrically arranged around the midpoints of their respective ranges, reducing complexity when optimizing using a GA. The specifics of each normalized input and output of the FISs are provided in Table 2.
The rule base is developed based on basic human intuition regarding potential collision scenarios. For instance, in FIS_H1:
  • IF x ^ is Close, meaning that the own-ship and the intruder are nearly colliding, AND Δ ψ ^ M is Left, meaning that the own-ship is positioned to the left of the intruder at CPA, AND Δ χ ^ d is Left, meaning that the desired look-ahead point is to the left of the own-ship, THEN Δ χ ^ ca is Go-Left, signifying that the own-ship’s path is to be modulated to the left.
  • IF x ^ is Far, meaning that the own-ship and the intruder are still considerably apart, AND Δ ψ ^ M is Center, meaning that the own-ship is nearly colliding head-on with the intruder at CPA, AND Δ χ ^ d is Center, meaning that the desired look-ahead point is directly ahead of the own-ship, THEN Δ χ ^ ca is Continue, indicating that the own-ship’s path should not be significantly modulated.
Similarly, in FIS_V2:
  • IF ν ^ z indicates an Upward motion, implying that the own-ship is ascending relative to the intruder, AND z ^ denotes Above, meaning that the intruder is positioned higher than the own-ship, THEN Δ ν ^ z ca signifies Descend, meaning that the own-ship’s vertical velocity component should be adjusted to descend.
  • IF ν ^ z is a Downward motion, meaning that the own-ship is descending relative to the intruder, AND z ^ signifies In-line, meaning that the intruder is directly ahead of the own-ship, THEN Δ ν ^ ) z ca signifies Climb, implying that the own-ship’s vertical velocity component should be adjusted to climb.
Thus, a total of 81 rules were defined—27 rules for FIS_H1 and FIS_V1 each and 9 rules for FIS_H2, FIS_V2, and FIS_W each.
It is also worth noting that, to simulate a real-life environment where perfectly symmetric conflict scenarios are unlikely, slight randomness is introduced to the relative variables Δ ψ M and z. Similarly, for special ‘Top-Down’ scenarios in which the intruder is directly above or below the own-ship and both of them are moving in the same direction towards a potential collision, the modulation parameters for vertical maneuvers, Δ χ ca and Δ h ca , are adjusted as follows:
Δ ν z ca = T D factor × Δ ν z ca ,
Δ h ca = T D factor × Δ h ca ,
where T D factor is an arbitrary factor that can be tuned for this special ‘Top-Down’ case. An instance of this scenario is illustrated by the pairwise conflict involving UAVs 5 and 6 in Figure 9a.

3.3.2. Desired Path Modulation

The desired path obtained from the waypoint follower is modulated using modulation parameters obtained from 3D-CAFIS. First, χ c and h c are calculated as
χ c = χ d + Υ x Δ χ ca ,
h c = h d + Υ z Δ h ca .
Now, the modulation of the desired airspeed begins with the determination of the desired velocity in the M-frame as follows:
M ν x d , ν y d , ν z d = C M / N C N / O O ν air d , 0 , 0 .
Then, ν x d and ν z d are modulated to obtain the following commanded horizontal and vertical closure rates:
ν x c = ν x d + Υ x Δ ν x ca ,
ν z c = ν z d + Υ z Δ ν z ca .
To avoid extreme deviation from the desired path, one sets ν y c = ν y d . The commanded airspeed ( ν air c ) is obtained from the norm of modulated velocity given below:
ν air c = C M N ν x c , ν y c , ν z c .

3.4. Controller

The controller block determines the primary control commands: γ c , ϕ c , and ν air c . It is important to note that when no potential collisions are detected, the values are such that ν air c = ν air d , χ c = χ d , and h c = h d . However, in the presence of a potential collision, γ c and ϕ c are determined from the path characterized by ν air c , χ c , and h c obtained from the desired path modulation as follows:
cos ( χ ψ ) ν ground tan ϕ c = k χ χ c χ ,
ν air sin ( γ c ) = min max k z h c h , ν air , ν air .
Here, k h and k χ are gain values tuned to achieve smoother bank and flight path angles. Lastly, γ c , ϕ c , and ν air c are the control commands used to control the UAV toward the desired waypoint while avoiding collision. Equation (30) represents a saturation function designed to limit the excessive vertical speed of the UAV. However, before being transmitted to the flight controller, these control commands are subjected to limits set by the physical and aerodynamic limitations of respective UAVs.
The algorithm, presented and discussed in Section 3.1, Section 3.2, Section 3.3 and Section 3.4, is summarized in Figure 7.

3.5. 3D-CAFIS Optimization Using the Genetic Algorithm (GA)

The GA is an evolutionary optimization algorithm that tunes system parameters by maximizing a fitness function (or minimizing a cost function) via the selection, crossover, and mutation of chromosomes comprising those parameters. Figure 8 depicts a fundamental FIS optimization diagram employing a GA. The rule base and the MFs constitute the knowledge base of the FIS. To fully define the proposed fuzzy tree, a total of 224 parameters must be tuned. These parameters include 4 for defining 2 trapezoidal MFs and 3 for defining 1 triangular MF for each of the 13 variables within the fuzzy tree, resulting in a total of 143 parameters, as well as 81 parameters corresponding to each of the rules. The computational effort required for optimizing the proposed FIS is significantly reduced through the utilization of constant rules chosen based on human intuition. Additionally, the symmetry of MFs within each FIS is maintained by considering only two tunable parameters (a and b), as shown in Figure 6. Consequently, the total number of tunable parameters is reduced to 26.
It is important to note that this optimization process is carried out before a specific mission type for a given UAV. Once optimized, the 3D-CAFIS algorithm is seamlessly integrated into the UAV system, enabling real-time CA during potentially conflicting scenarios.
This study involves the formulation of a Cost Function ( C F ) aimed at minimizing collisions, reducing path deviation during CA, and ensuring adequate separation to prevent potential collisions and preemptive collisions. The C F is defined as follows:
min C F = t = 0 t f δ 1 σ 1 ( t ) + δ 2 σ 2 ( t ) + δ 3 σ 3 ( t ) + δ 4 σ 4 ( t ) ,
where t f is the total simulation time. σ 1 , σ 2 , σ 3 , and σ 4 are the collision cost to prevent a collision, the route cost to minimize the total travel distance, the separation cost to maintain proximity to the given separation threshold ( s th ), and the time of collision cost to avoid early collisions, respectively. Correspondingly, δ i serves as the weight assigned to each respective cost, which can be tuned to achieve mission-critical performance objectives. Each of these costs at a given instant t is calculated as follows:
σ 1 ( t ) = i = 1 n coll i ,
σ 2 ( t ) = i = 1 n uavs ρ ( w / O ) i ,
σ 3 ( t ) = i = 1 n uavs s th s i ,
σ 4 ( t ) = i = 1 n coll t f t coll , i ,
where n uavs is the total number of UAVs in simulation, n coll is the total number of colliding UAVs, t coll , i is the time of collision for the i-th colliding UAV, and s th = x th , z th is the separation threshold.

4. Simulation Results and Discussion

4.1. Simulation Environment

The UAVs are considered to be operating in a cubicle space of ( 2 nm × 2 nm × 1 nm ) , moving from their arbitrary start positions at one end toward a designated waypoint at the opposite end. Such a small space is considered to evaluate the performance of the proposed algorithm in the most extreme situations. If the algorithm performs well even in severely congested airspace, the chances of collision in real life will be next to none. The UAVs are programmed to decelerate as they approach their respective waypoints. It is assumed that all UAVs are equipped with essential sensors and are operating the proposed CA algorithm. The simulations span a period ( t f ) of 130 s. The sampling time ( Δ t ) for both simulation and UAV control is set at 0.1 s. Parameters such as gains, UAV state limits, and thresholds utilized in the simulations are presented in Table 3.

4.2. 3D-CAFIS Optimization Environment

The optimization of the FIS using the GA involved a training scenario consisting of three pairwise conflicts as shown in Figure 9a. In these scenarios, the UAV pairs [ 1 , 2 ] , [ 3 , 4 ] , and [ 5 , 6 ] are each on course for a potential collision, approaching from distinct directions as they navigate toward their respective waypoints. During the training process, the range-bearing sensor Field of View (FOV) is assumed to be wider, defined by the conditions r 0.4 nm, η 90 degrees, and ϵ 90 degrees The GA parameters employed for optimizing FIS parameters are listed in Table 4.

4.3. 3D-CAFIS Optimization Results

Figure 9 shows the performance disparity among UAVs in the training scenario under three different conditions: No CA (NCA), manually defined 3D-CAFIS (MCA), and optimized 3D-CAFIS using the GA (GCA). As expected, during NCA, collisions occur, resulting in an escalating total operational cost over time with each collision. However, both MCA and GCA successfully avoided all collisions. Remarkably, the cost of operation, as calculated using Equation (31), was 15% lower with GCA than with MCA. This reduction can be attributed to GCA’s ability to maintain a greater UAV separation during CA, as evident from the significant decrease in σ 3 compared to MCA.
The inputs and outputs of FIS_H1 for both MCA and GCA are shown in Figure 10. The figure demonstrates how the optimization of the fuzzy tree altered the shape of MFs of the inputs and outputs of all the FISs to achieve the desired performance. In Figure 10a, it is evident that all the FISs were initially defined in a similar way, where the parameters of the MFs of each variable were chosen to be quite close to the center of the range of the respective variables. This was performed to maximize the values of modulation parameters to avoid potential collisions. The optimization of the fuzzy tree using GA, for the weights ( δ s) of C F tabulated in Table 3, resulted in modified MFs, as shown in Figure 10b. One can now make an inference regarding how the input variables have an impact on their respective output variables. For instance, here, it seems that the ‘Left’ and ‘Right’ MFs of Δ ψ ^ M have more influence than the ‘Center’ MF. Similarly, the ‘Go-Left’ and ‘Go-Right’ MFs of Δ χ ^ ca have some slope in them, which should induce smoother changes to χ c for a given x, Δ ψ M , and Δ χ d .
In the training scenario, there was a conflict between UAV 3 and UAV 4, as shown in Figure 11a. During NCA, both UAVs were descending towards their respective waypoints, resulting in zero vertical separation and horizontal separation dropping below the thresholds, indicating a possible collision. However, with MCA and GCA as shown in Figure 11b,c, respectively, horizontal separation remained above the threshold, preventing collisions. MCA resulted in a slight increase in vertical separation from 57.7 s to 66.6 s, which was not observed with GCA. GCA achieved most of the CA through horizontal maneuvering and concluded at 64.8 s, which was shorter than MCA. After achieving CA, both UAVs proceeded towards their respective waypoints.
Figure 12 shows how the states of the UAV 3 in the training scenario compare between the cases of MCA and GCA. In the first 5 s, the UAV changes its direction towards its waypoint as it moves away from it, leading to corrections in its ϕ and χ . The main CA phase commences at 58 s for both cases. However, with GCA, the fluctuations in states are induced due to CA ending earlier. Moreover, the UAV banks and changes flight path angles more smoothly with GCA than with MCA. In both cases, the UAV banks left and then right, which is evident from the bank angle ( ϕ ) and course angle ( χ ) plots in Figure 12. The adjustments in velocity and altitude for this conflict scenario also appear to be minimal with GCA compared to MCA. The CA for UAVs 3 and 4 in the training scenario with GCA is visualized in Figure 13. It is evident that the UAVs bank away from each other, forming an ‘S’-shaped 3D path to avoid collision. This shows that for this scenario, the CA was achieved mostly with horizontal maneuvering compared to vertical.

4.4. Optimized 3D-CAFIS Testing Results

4.4.1. Pairwise Conflict Scenario

In order to test the GCA, multiple pairwise conflict scenarios were evaluated, including one illustrated in Figure 14. This particular scenario involved UAVs 1 and 2 moving towards each other at an almost perpendicular angle. To assess the scalability of 3D-CAFIS, both UAVs had a narrower FOV, with conditions set at r 0.4 nm, η 60 degrees, and ϵ 60 degrees. These specifications were based on the EchoFlight radar sensor by Echodyne Radar. The separation plot in Figure 14a shows that when using NCA, the horizontal and vertical separation between the UAVs were below the threshold of 60 to 70 s, indicating a potential collision. However, when simulated with MCA (Figure 14b) and GCA (Figure 14c), both horizontal and vertical separation were above the thresholds at CPA. The CA was initiated at 50.9 s for both MCA and GCA but ended earlier with GCA at 61.5 s, compared to 62.6 s with MCA. Additionally, during the CA, the vertical separation was slightly greater with GCA than with MCA. Therefore, unlike the training scenario of UAVs 3 and 4 discussed above, GCA increased vertical separation in this scenario.
Figure 15 shows the comparison between the states of UAV 1 in the testing scenario for both MCA and GCA. During the first 5 s, the UAV changes its direction towards its waypoint, correcting its angle ( ϕ ) and course angle ( χ ) due to its initial divergence from the path. The main CA begins at 51 s for both cases. Similar to UAV 3 in the training scenario, here as well, the fluctuations in states are influenced by the earlier termination of CA by a second with GCA. Here too, the ϕ and χ plots show the UAV makes a noticeable horizontal maneuver, this time to the right and then to the left. Moreover, for this scenario, vertical maneuvering was also pronounced with both MCA and GCA. The airspeed ( ν air ) reduction is slightly more pronounced with GCA compared to MCA, while the altitude adjustment remains quite similar. Consequently, the added vertical separation with GCA could be associated with greater airspeed. Also, the UAV exhibits smoother banking maneuvers with GCA in comparison to MCA. This shows that the algorithm can perform a combination of horizontal and vertical maneuvers for CA based on given conflict scenarios.
The CA for the UAVs 1 and 2 in the testing scenario using GCA is visualized in Figure 16. One can see that the UAVs execute banking maneuvers to move away from each other, resulting in the formation of an ‘S’-shaped 3D path to successfully avoid collision.

4.4.2. Monte Carlo Simulations

Similarly, to evaluate the effectiveness of the algorithm over a multitude of scenarios, a Monte Carlo simulation was conducted. This simulation involved multiple UAVs operating within a confined space and moving in various directions. The UAVs were positioned at one end of the region, while their targets were at the opposite end. To begin with, the UAVs were set to operate within a range of ν = ( 40 , 80 ) knots and γ = ( 6 , 6 ) degrees, with all ϕ values set to 0 degrees. The initial χ value was determined based on the starting position of the UAVs. For instance, if the starting position was west of the confined space, χ was set between 85 and 95 degrees with respect to north, and if the starting position was east, χ was set between 265 and 275 degrees. Note that in the Monte Carlo simulation, the FOV for UAVs is assumed to be wider to evaluate the effectiveness of the algorithm when UAVs are capable of detecting almost all intruders.
The Monte Carlo simulation included 15 runs, each with randomly generated sets of 40, 60, 80, and 100 UAVs. That is, the initial states and target positions of the UAVs were randomly generated for each of the 15 runs for each set number of UAVs. This resulted in a total of 60 simulation runs performed with NCA, MCA, and GCA each, which included 4200 distinctly initialized UAVs. An example of one such simulation run is shown in Figure 17a, where 16 collisions occurred without CA, all of which are effectively prevented with GCA, as shown in Figure 17b. The simulation results are summarized in Table 5. As expected, the higher the number of UAVs in the airspace, the greater the number of collided UAVs. However, when employing MCA, the overall number of collisions reduces by an average of 98%. With GCA, this number drops even further to an average of 99.5 %. Moreover, the average total cost decreases by 15% when using GCA compared to MCA. This indicates that the proposed 3D-CAFIS is effective across diverse conflict scenarios, particularly when optimized using the GA.
It is important to note that instances of collisions still occurred when UAVs were outside of the FOV for detection and avoidance. Nevertheless, one can infer that the collisions can be effectively avoided using the proposed algorithm for most cases, as the simulations show that despite operating hundreds of UAVs in an extremely small space, the average number of collisions was near zero.

5. Conclusions

The proposed work introduces a Three-Dimensional (3D) Collision Avoidance Fuzzy Inference System (3D-CAFIS) designed for fixed-wing Unmanned Aerial Vehicles (UAVs), which is further enhanced through optimization using the Genetic Algorithm (GA). The algorithm provides a combination of horizontal and vertical resolutions to mitigate potential collisions. The 3D-CAFIS generates linguistically characterized path modulation parameters, ensuring comprehensibility and interpretability. These parameters modulate the desired path upon identifying impending collisions. The fuzzy tree used in 3D-CAFIS is optimized using the GA to ensure that the distance traveled during the mission remains minimal, even with path modulation, while ensuring that separation during CA remains above the designated threshold. In the training phase, three pairwise conflicts were considered to optimize the fuzzy tree. While the manually designed 3D-CAFIS demonstrated effective CA, the optimized 3D-CAFIS showcased superior performance, exhibiting greater efficacy in multiple UAV scenarios within confined airspace. This was evaluated through numerous pairwise conflicts and even a Monte Carlo simulation, where a multitude of randomly generated UAVs in varying numbers were simulated to operate in closed airspace multiple times. It should be noted that the algorithm is not claimed to outperform other CA algorithms. Instead, it presents an alternative approach that is characterized by its simplicity and the linguistic logic of the proposed fuzzy tree. This method can be further refined through manual tuning or optimization using a GA, and it has the flexibility to be scaled according to the performance capability of UAVs and sensor specifications. The results could be further enhanced by considering a more refined cost function that takes into account the effects of energy consumption during horizontal and vertical CA maneuvers. Therefore, future research will focus on evaluating energy effects, addressing uncertainties such as wind, and developing a modified version of the algorithm for rotary-wing UAVs, allowing for hardware verification in indoor settings.

Author Contributions

Conceptualization, S.R.; methodology, S.R. and D.K.; software, S.R.; validation, S.R.; formal analysis, S.R. and D.K.; investigation, S.R.; resources, D.K.; data curation, S.R.; writing—original draft preparation, S.R.; writing—review and editing, S.R. and D.K.; visualization, S.R.; supervision, D.K.; project administration, D.K.; funding acquisition, D.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

For convenience, the definitions of the various variables, symbols, and frames used are provided in the following abbreviations.
SymbolDefinitionUnit
A / B the A-frame with respect to the B-frame-
cthe command-
collthe colliding-
C R 3 × 3 the direction cosine matrix-
cathe collision avoidance-
C F the cost function-
dthe desired-
fthe forceN
gthe acceleration due to gravitynm/ s 2
hthe altitude of the UAVft
kthe gain
lthe look-ahead distancenm
nthe number-
optthe optimal-
p R 3 the state vectornm
rthe range of the intruder from the own-shipnm
s R 3 the separation vector from the O-frame to the I-framenm
tthe times
ththe threshold-
xthe separation vector component along m ^ 1 nm
ythe separation vector component along m ^ 2 nm
zthe separation vector component along m ^ 3 ft
wthe waypoint-
α the angle of attack of the UAVdeg
β the side-slip of the UAVdeg
γ the flight path angle of the UAVdeg
δ the weighting factor in the cost function-
ϵ the elevation of the intruder from the own-shipdeg
η the azimuth of the intruder from the own-shipdeg
θ the UAV pitch angledeg
Λ R 3 the 3-2-1 set of Euler anglesdeg
ν R 3 the velocity vector of the UAVkts
ρ R 3 the position vector of the UAVnm
σ the cost element in the cost functionnm
τ the time to the Closest Point of Approach (CPA)s
Υ the weight for path modulation in the Fuzzy Inference System (FIS)%
ϕ the UAV roll angledeg
χ the course angle of the UAVdeg
ψ the UAV yaw angledeg
ω R 3 the angular velocity vector of the UAVdeg/s
N : n ^ 1 , n ^ 2 , n ^ 3 the inertial frame
I : i ^ 1 , i ^ 2 , i ^ 3 the intruder frame
G : g ^ 1 , g ^ 2 , g ^ 3 the ground frame obtained after rotating the N-frame via the 3-2-1 set of Euler angles χ , γ , and ϕ
M : m ^ 1 , m ^ 2 , m ^ 3 the intermediate frame obtained after rotating the N-frame about n ^ 3 so that m ^ 1 is facing towards the nearest intruder while n ^ 3 = m ^ 3
O : o ^ 1 , o ^ 2 , o ^ 3 the own-ship frame obtained after rotating the N-frame via the 3-2-1 set of Euler angles ψ , θ , and ϕ
( . ) the left superscript represents the frame in which the vector in context
is expressed
( . ) the left subscript represents the frame in which the vector in context is observed
( . ) the right superscript provides additional state information
( . ) the transpose
( . ) the right subscript provides time, position, and frame information

References

  1. Cooper, T.; Smiley, J.; Porter, C.; Precourt, C. Global Fleet & MRO Market Forecast 2022–2032; Technical Report; Oliver Wyman: New York, NY, USA, 2022. [Google Scholar]
  2. Advanced Aerial Mobility Market (By Mode of Operation: Piloted, Autonomous; By End-Use: Cargo, Passenger; By Propulsion Type: Parallel Hybrid, Turboshaft, Electric, Turboelectric)—Global Industry Analysis, Size, Share, Growth, Trends, Regional Outlook, and Forecast 2022–2030. Technical Report 2110. 2022. Available online: https://www.precedenceresearch.com/advanced-aerial-mobility-market (accessed on 19 February 2023).
  3. FAA. Introduction to TCAS II: Version 7.1; Technical Report; Federal Aviation Administration (FAA): Washington, DC, USA, 2011. [Google Scholar]
  4. Manfredi, G.; Jestin, Y. An introduction to ACAS Xu and the challenges ahead. In Proceedings of the 2016 IEEE/AIAA 35th Digital Avionics Systems Conference (DASC), Sacramento, CA, USA, 25–29 September 2016; pp. 1–9. [Google Scholar] [CrossRef]
  5. Tang, J. Review: Analysis and Improvement of Traffic Alert and Collision Avoidance System. IEEE Access 2017, 5, 21419–21429. [Google Scholar] [CrossRef]
  6. Owen, M.P.; Panken, A.; Moss, R.; Alvarez, L.; Leeper, C. Acas xu: Integrated collision avoidance and detect and avoid capability for uas. In Proceedings of the 2019 IEEE/AIAA 38th Digital Avionics Systems Conference (DASC), Sacramento, CA, USA, 8–12 September 2019; pp. 1–10. [Google Scholar] [CrossRef]
  7. Alvarez, L.E.; Jessen, I.; Owen, M.P.; Silbermann, J.; Wood, P. ACAS sXu: Robust decentralized detect and avoid for small unmanned aircraft systems. In Proceedings of the 2019 IEEE/AIAA 38th Digital Avionics Systems Conference (DASC), Sacramento, CA, USA, 8–12 September 2019; pp. 1–9. [Google Scholar] [CrossRef]
  8. Kuchar, J.K.; Yang, L.C. A review of conflict detection and resolution modeling methods. IEEE Trans. Intell. Transp. Syst. 2000, 1, 179–189. [Google Scholar] [CrossRef]
  9. Albaker, B.; Rahim, N. A survey of collision avoidance approaches for unmanned aerial vehicles. In Proceedings of the 2009 International Conference for Technical Postgraduates (TECHPOS), Kuala Lumpur, Malaysia, 14–15 December 2009; pp. 1–7. [Google Scholar] [CrossRef]
  10. Huang, S.; Teo, R.S.H.; Tan, K.K. Collision avoidance of multi unmanned aerial vehicles: A review. Annu. Rev. Control 2019, 48, 147–164. [Google Scholar] [CrossRef]
  11. Trapani, A.; Erzberger, H.; Dunbar, W. Performance Analysis of a Horizontal Separation Assurance Algorithm for Short-Range Conflicts. In Proceedings of the AIAA Guidance, Navigation, and Control Conference, Chicago, IL, USA, 10–13 August 2009; p. 5746. [Google Scholar]
  12. Guo, J.; Liang, C.; Wang, K.; Sang, B.; Wu, Y. Three-Dimensional Autonomous Obstacle Avoidance Algorithm for UAV Based on Circular Arc Trajectory. Int. J. Aerosp. Eng. 2021, 2021, 1–13. [Google Scholar] [CrossRef]
  13. Yang, X.; Alvarez, L.M.; Bruggemann, T. A 3D Collision Avoidance Strategy for UAVs in a Non-Cooperative Environment. J. Intell. Robot. Syst. 2013, 70, 315–327. [Google Scholar] [CrossRef]
  14. Peng, L.; Lin, Y. A closed-form solution of horizontal maneuver to collision avoidance system for UAVs. In Proceedings of the 2010 Chinese Control and Decision Conference, Xuzhou, China, 26–28 May 2010; pp. 4416–4421. [Google Scholar] [CrossRef]
  15. van den Berg, J.; Lin, M.; Manocha, D. Reciprocal Velocity Obstacles for real-time multi-agent navigation. In Proceedings of the 2008 IEEE International Conference on Robotics and Automation, Pasadena, CA, USA, 19–23 May 2008; pp. 1928–1935. [Google Scholar] [CrossRef]
  16. Zhang, S.; Xu, T.; Cheng, H.; Liang, F. Collision Avoidance of Fixed-Wing UAVs in Dynamic Environments Based on Spline-RRT and Velocity Obstacle. In Proceedings of the 2020 International Conference on Unmanned Aircraft Systems (ICUAS), Athens, Greece, 1–4 September 2020; pp. 48–58. [Google Scholar] [CrossRef]
  17. Qu, Y.; Yi, W. Three-Dimensional Obstacle Avoidance Strategy for Fixed-Wing UAVs Based on Quaternion Method. Appl. Sci. 2022, 12, 955. [Google Scholar] [CrossRef]
  18. Khare, I.; Kandath, H.; Krishna, K.M. Predictive optimal collision avoidance for a swarm of fixed-wing aircraft in 3D space. In Proceedings of the 2022 International Conference on Unmanned Aircraft Systems (ICUAS), Dubrovnik, Croatia, 21–24 June 2022; IEEE: Piscataway, NJ, USA, 2022; pp. 1139–1145. [Google Scholar] [CrossRef]
  19. Shin, H.; Tsourdos, A.; White, B. UAS Conflict Detection and Resolution Using Differential Geometry Concepts. In Proceedings of the Sense and Avoid in UAS; John Wiley & Sons, Ltd.: Chichester, UK, 2012; pp. 175–204. [Google Scholar] [CrossRef]
  20. Chao, X.; Wu-Jie, X.; Wen-Han, D. The design of automatic air collision avoidance system based on geometric. In Proceedings of the 2018 Chinese Control And Decision Conference (CCDC), Shenyang, China, 9–11 June 2018; pp. 240–245. [Google Scholar] [CrossRef]
  21. Beary, A. Fixed Wing Aircraft Collision Avoidance Using Collision Cone Theory. Ph.D. Thesis, University of Texas, Arlington, TX, USA, August 2022. [Google Scholar]
  22. Zhu, L.; Cheng, X.; Yuan, F.G. A 3D collision avoidance strategy for UAV with physical constraints. Meas. J. Int. Meas. Confed. 2016, 77, 40–49. [Google Scholar] [CrossRef]
  23. Zhang, J.; Yan, J.; Zhang, P.; Kong, X. Collision Avoidance in Fixed-Wing UAV Formation Flight Based on a Consensus Control Algorithm. IEEE Access 2018, 6, 43672–43682. [Google Scholar] [CrossRef]
  24. Srivastava, A.; Sujit, P.B. A reformulation of collision avoidance algorithm based on artificial potential fields for fixed-wing UAVs in a dynamic environment. arXiv 2022, arXiv:2208.03117. [Google Scholar] [CrossRef]
  25. Wang, F.; Lu, Z.; Liu, Q.; Zhu, Y.; Zhang, J.; Shi, G. Research on formation collision avoidance of aircraft cooperative penetration based on improved potential field method. In Proceedings of the 2020 3rd International Conference on Unmanned Systems (ICUS), Harbin, China, 27–28 November 2020; pp. 552–557. [Google Scholar] [CrossRef]
  26. Choi, D.; Lee, K.; Kim, D. Enhanced potential field-based collision avoidance for unmanned aerial vehicles in a dynamic environment. In Proceedings of the AIAA Scitech 2020 Forum, Orlando, FL, USA, 6–10 January 2020; p. 0487. [Google Scholar] [CrossRef]
  27. Vanneste, S.; Bellekens, B.; Weyn, M. 3DVFH+: Real-time three-dimensional obstacle avoidance using an octomap. In Proceedings of the 1st International Workshop on Model-Driven Robot Software Engineering co-located with International Conference on Software Technologies: Applications and Foundations (STAF 2014), York, UK, 21 July 2014; Uwe, A., Ed.; 2014. [Google Scholar]
  28. Stepanyan, V.; Chakrabarty, A.; Ippolito, C.A. Pursuit-Evasion Game approach to Tree-Based Path Planning for Airborne Dynamic Obstacle Avoidance. In Proceedings of the AIAA Scitech 2021 Forum, Virtual, 11–21 January 2021; p. 1575. [Google Scholar] [CrossRef]
  29. Zhao, Y.; Guo, J.; Bai, C.; Zheng, H. Reinforcement Learning-Based Collision Avoidance Guidance Algorithm for Fixed-Wing UAVs. Complexity 2021, 2021, 1–12. [Google Scholar] [CrossRef]
  30. Machmudah, A.; Shanmugavel, M.; Parman, S.; Manan, T.S.A.; Dutykh, D.; Beddu, S.; Rajabi, A. Flight Trajectories Optimization of Fixed-Wing UAV by Bank-Turn Mechanism. Drones 2022, 6, 69. [Google Scholar] [CrossRef]
  31. Kosari, A.; Teshnizi, M.M. Optimal Trajectory Design for Conflict Resolution and Collision Avoidance of Flying Robots using Radau-Pseudo Spectral Approach. In Proceedings of the 2018 6th RSI International Conference on Robotics and Mechatronics (IcRoM), Tehran, Iran, 23–25 October 2018; pp. 82–87. [Google Scholar] [CrossRef]
  32. Prandini, M.; Hu, J.; Lygeros, J.; Sastry, S. A probabilistic approach to aircraft conflict detection. IEEE Trans. Intell. Transp. Syst. 2000, 1, 199–220. [Google Scholar] [CrossRef]
  33. Riviere, B.; Honig, W.; Yue, Y.; Chung, S.J. GLAS: Global-to-Local Safe Autonomy Synthesis for Multi-Robot Motion Planning With End-to-End Learning. IEEE Robot. Autom. Lett. 2020, 5, 4249–4256. [Google Scholar] [CrossRef]
  34. Berdonosov, V.; Zivotova, A.; Zhuravlev, D.; Naing, Z.H. Implementation of the speed approach for UAV collision avoidance in dynamic environment. In Proceedings of the 2018 International Multi-Conference on Industrial Engineering and Modern Technologies (FarEastCon), Vladivostok, Russia, 2–4 October 2018; pp. 1–6. [Google Scholar] [CrossRef]
  35. Zhou, D.; Wang, Z.; Bandyopadhyay, S.; Schwager, M. Fast, On-line Collision Avoidance for Dynamic Vehicles Using Buffered Voronoi Cells. IEEE Robot. Autom. Lett. 2017, 2, 1047–1054. [Google Scholar] [CrossRef]
  36. Aiyama, Y.; Ishiwatari, Y.; Seki, T. Human behavior based obstacle avoidance for human-robot cooperative transportation. In Proceedings of the Mechatronics for Safety, Security and Dependability in a New Era; Elsevier Ltd.: Amsterdam, The Netherlands, 2006; pp. 177–180. [Google Scholar] [CrossRef]
  37. Hawkins, K.; Tsiotras, P. Anticipating human collision avoidance behavior for safe robot reaction. In Proceedings of the 2018 IEEE Conference on Decision and Control (CDC), Miami Beach, FL, USA, 17–19 December 2018; pp. 6301–6306. [Google Scholar] [CrossRef]
  38. Cai, C.; Yang, C.; Zhu, Q.; Liang, Y. A fuzzy-based collision avoidance approach for multi-robot systems. In Proceedings of the 2007 IEEE International Conference on Robotics and Biomimetics (ROBIO), Sanya, China, 15–18 December 2007; pp. 1012–1017. [Google Scholar] [CrossRef]
  39. Yongjie, Y.; Yan, Z. The collision avoidance planning in multi-robot system by genetic fuzzy control algorithm. In Proceedings of the 2008 IEEE International Conference on Robotics and Biomimetics, Bangkok, Thailand, 22–25 February 2009; pp. 1020–1025. [Google Scholar] [CrossRef]
  40. Kamthan, S.; Singh, H.; Meitzler, T. Survivability: A hierarchical fuzzy logic layered model for threat management of unmanned ground vehicles. In Proceedings of the Unmanned Systems Technology XIX, SPIE, Orlando, FL, USA, 3 May 2018; Volume 10643, p. 106430W-11. [Google Scholar] [CrossRef]
  41. Najmurrokhman, A.; Kusnandar; Komarudin, U.; Sunubroto.; Sadiyoko, A.; Iskanto, T.Y. Mamdani based Fuzzy Logic Controller for A Wheeled Mobile Robot with Obstacle Avoidance Capability. In Proceedings of the 2019 International Conference on Mechatronics, Robotics and Systems Engineering (MoRSE), Bali, Indonesia, 4–6 December 2019; pp. 49–53. [Google Scholar] [CrossRef]
  42. Younas, I.; Ahmed, Z.; Mohyud-Din, S.T. A Fuzzy Based Aircraft Collision Avoidance System. In Proceedings of the 9th WSEAS International Conference on Applied Informatics and Communications, AIC’09, Moscow, Russia, 20–22 August 2009; World Scientific and Engineering Academy and Society (WSEAS): Stevens Point, WI, USA, 2009; pp. 344–350. [Google Scholar]
  43. Choi, D.; Chhabra, A.; Kim, D. Collision Avoidance of Unmanned Aerial Vehicles Using Fuzzy Inference System-Aided Enhanced Potential Field. In Proceedings of the AIAA SCITECH 2022 Forum, San Diego, CA, USA, 3–7 January 2022; p. 0272. [Google Scholar] [CrossRef]
  44. Hromatka, M. A Fuzzy Logic Approach to Collision Avoidance in Smart UAVs. Ph.D Thesis, St. John’s University, Collegeville, MN, USA, 12 August 2013. [Google Scholar]
  45. Du, B.; Liu, S. A common obstacle avoidance module based on fuzzy algorithm for Unmanned Aerial Vehicle. In Proceedings of the 2016 IEEE International Conference on Cyber Technology in Automation, Control, and Intelligent Systems (CYBER), Chengdu, China, 19–22 June 2016; pp. 245–248. [Google Scholar] [CrossRef]
  46. Wang, K.; Rong, H.; Shi, M.; Meng, Z.; Zhou, Y. Reactive Obstacle Avoidance for Unmanned Helicopter Using Fuzzy Control. In Proceedings of the 2018 37th Chinese Control Conference (CCC). Technical Committee on Control Theory, Chinese Association of Automation, Wuhan, China, 25–27 July 2018; pp. 10026–10031. [Google Scholar] [CrossRef]
  47. Cook, B.; Arnett, T.; Cohen, K. A Fuzzy Logic Approach for Separation Assurance and Collision Avoidance for Unmanned Aerial Systems. Mod. Fuzzy Control. Syst. Its Appl. 2017, 225–256. [Google Scholar] [CrossRef]
  48. Beard, R.W.; Mclain, T.W. Small Unmanned Aircraft; Princeton University Press: Princeton, NJ, USA, 2012. [Google Scholar] [CrossRef]
  49. Munoz, C.; Narkawicz, A.; Chamberlain, J. A TCAS-II Resolution Advisory Detection Algorithm. In Proceedings of the AAIA Conference Papers American Institute of Aeronautics and Astronautics, Grapevine, TX, USA, 7–10 January 2013. [Google Scholar] [CrossRef]
Figure 1. The own-ship (O) frame and the inertial (N) frame.
Figure 1. The own-ship (O) frame and the inertial (N) frame.
Electronics 12 03946 g001
Figure 2. UAV states in (a) horizontal plane and (b) vertical plane.
Figure 2. UAV states in (a) horizontal plane and (b) vertical plane.
Electronics 12 03946 g002
Figure 3. UAV system architecture integrated with the proposed CA algorithm.
Figure 3. UAV system architecture integrated with the proposed CA algorithm.
Electronics 12 03946 g003
Figure 4. Collision scenario: (a) horizontal view and (b) vertical view.
Figure 4. Collision scenario: (a) horizontal view and (b) vertical view.
Electronics 12 03946 g004
Figure 5. Proposed fuzzy tree for CA.
Figure 5. Proposed fuzzy tree for CA.
Electronics 12 03946 g005
Figure 6. Tunable parameters a and b for MFs of (a) x and (b) Δ χ ca .
Figure 6. Tunable parameters a and b for MFs of (a) x and (b) Δ χ ca .
Electronics 12 03946 g006
Figure 7. Summary of the proposed CD and CA algorithms.
Figure 7. Summary of the proposed CD and CA algorithms.
Electronics 12 03946 g007
Figure 8. GA-based FIS optimization.
Figure 8. GA-based FIS optimization.
Electronics 12 03946 g008
Figure 9. Pairwise conflict training scenario along with their associated costs over time for (a) NCA, (b) MCA, and (c) GCA.
Figure 9. Pairwise conflict training scenario along with their associated costs over time for (a) NCA, (b) MCA, and (c) GCA.
Electronics 12 03946 g009
Figure 10. MFs of inputs and output of FIS_H1 for (a) MCA and (b) GCA.
Figure 10. MFs of inputs and output of FIS_H1 for (a) MCA and (b) GCA.
Electronics 12 03946 g010
Figure 11. UAVs 3 and 4 in the training scenario along with their separation over time with (a) NCA, (b) MCA, and (c) GCA.
Figure 11. UAVs 3 and 4 in the training scenario along with their separation over time with (a) NCA, (b) MCA, and (c) GCA.
Electronics 12 03946 g011
Figure 12. States of UAV 3 in the training scenario with (a) MCA and (b) GCA; current states (Electronics 12 03946 i001), commanded states (Electronics 12 03946 i002), and waypoint position (Electronics 12 03946 i003).
Figure 12. States of UAV 3 in the training scenario with (a) MCA and (b) GCA; current states (Electronics 12 03946 i001), commanded states (Electronics 12 03946 i002), and waypoint position (Electronics 12 03946 i003).
Electronics 12 03946 g012
Figure 13. Visualization of CA of UAVs 3 and 4 in the training scenario with GCA.
Figure 13. Visualization of CA of UAVs 3 and 4 in the training scenario with GCA.
Electronics 12 03946 g013
Figure 14. Pairwise conflict testing scenario along with their separation over time with (a) NCA, (b) MCA, and (c) GCA.
Figure 14. Pairwise conflict testing scenario along with their separation over time with (a) NCA, (b) MCA, and (c) GCA.
Electronics 12 03946 g014
Figure 15. States of UAV 1 in the testing scenario with (a) MCA and (b) GCA; current states (Electronics 12 03946 i001), commanded states (Electronics 12 03946 i002), and waypoint position (Electronics 12 03946 i003).
Figure 15. States of UAV 1 in the testing scenario with (a) MCA and (b) GCA; current states (Electronics 12 03946 i001), commanded states (Electronics 12 03946 i002), and waypoint position (Electronics 12 03946 i003).
Electronics 12 03946 g015
Figure 16. Visualization of CA of UAVs 1 and 2 in the testing scenario using GCA.
Figure 16. Visualization of CA of UAVs 1 and 2 in the testing scenario using GCA.
Electronics 12 03946 g016
Figure 17. Multiple UAVs test case: (a) NCA and (b) GCA.
Figure 17. Multiple UAVs test case: (a) NCA and (b) GCA.
Electronics 12 03946 g017
Table 1. FIS tree input and output variables.
Table 1. FIS tree input and output variables.
VariablesTypeNormalization FactorSignificance
xInput 0.4 nmHorizontal proximity
zInput1562 ftVertical proximity
Δ χ d Input π Desired course
Δ ψ M Input π Course direction with respect to intruder
ν x Input140 ktsHorizontal closure rate
ν z Input6000 ft/minVertical closure rate
Δ z d * Input 42.65 ftDesired altitude
Δ ν x ca Output 1.738 ktsHorizontal speed modulation
Δ ν z ca Output 124.45 ft/minVertical speed modulation
Δ χ ca Output65 degCourse modulation
Δ h ca Output8 ftAltitude modulation
Υ x Output1Intensity of the horizontal maneuvers
Υ z Output1Intensity of the vertical maneuvers
Table 2. Normalized variable details.
Table 2. Normalized variable details.
VariablesInput forOutput ofRangeMFs
TrapezoidTriangleTrapezoid
x ^ FIS_H1, FIS_H2, FIS_V1None[0, 1]CloseNearFar
z ^ FIS_V1, FIS_V2None[ 1 , 1]AboveIn-lineBelow
Δ χ ^ d FIS_H1None[ 1 , 1]LeftCenterRight
Δ ψ ^ M FIS_H1None[ 1 , 1]LeftCenterRight
ν ^ x FIS_H2None[0, 1]SlowMediumFast
ν ^ z FIS_V2None[ 1 , 1]UpwardStraightDownward
Δ z ^ d * FIS_H1None[ 1 , 1]AboveIn-lineBelow
Δ ν ^ x ca NoneFIS_H1[ 1 , 1]Speed-DownContinueSpeed-Up
Δ ν ^ z ca NoneFIS_V1[ 1 , 1]ClimbContinueDescend
Δ χ ^ ca FIS_WFIS_H2[ 1 , 1]Go-LeftContinueGo-Right
Δ h ^ ca FIS_WFIS_V2[ 1 , 1]Go-DownContinueGo-Up
Υ x NoneFIS_W[0, 1]LowMediumHigh
Υ z NoneFIS_W[0, 1]LowMediumHigh
Table 3. Simulation parameters.
Table 3. Simulation parameters.
GainsValuesUAV StatesRangeThresholdsValues
l43 ft ν air opt [55, 65] kts τ th 25 s
k h 13 ν air [45, 70] ktsDMOD 0.4 nm
k γ 0.39 max ( λ ) 3.5 ZTHR850 ft
k ψ 1.5 γ [ 45 , 45] deg x th 197 ft
k ν air 1 ϕ [ 25 , 25] deg z th 50 ft
k ϕ 3402.97
k ϕ ˙ 116.67
T D factor 2
Table 4. GA parameters and C F parameters.
Table 4. GA parameters and C F parameters.
GA ParametersValues CF ParametersValues
Number of generations100 δ 1 300 , 000
Population size50 δ 2 10
Elitism ratio 0.05 δ 3 1000
Crossover fraction 0.8 δ 4 500
Selection algorithmTournament selection s th 0.033 nm
Crossover algorithmTwo-points crossover
Mutation algorithmAdaptive feasible
Table 5. Summary of Monte Carlo simulations.
Table 5. Summary of Monte Carlo simulations.
# of UAVsCA TypeAvg. # of Collided UAVsAvg. Total Cost (nm)
40NCA2.5 1.2051 × 10 7
MCA0.067 3.2245 × 10 5
GCA0.000 2.4129 × 10 5
60NCA6.5 2.1838 × 10 7
MCA0.067 8.7452 × 10 5
GCA0.000 7.5315 × 10 5
80NCA10.6 3.2484 × 10 7
MCA0.267 1.6881 × 10 6
GCA0.000 1.3948 × 10 6
100NCA14.8 4.3264 × 10 7
MCA0.333 2.8431 × 10 6
GCA0.267 2.4840 × 10 6
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

Rauniyar, S.; Kim, D. Genetic Fuzzy Inference System-Based Three-Dimensional Resolution Algorithm for Collision Avoidance of Fixed-Wing UAVs. Electronics 2023, 12, 3946. https://doi.org/10.3390/electronics12183946

AMA Style

Rauniyar S, Kim D. Genetic Fuzzy Inference System-Based Three-Dimensional Resolution Algorithm for Collision Avoidance of Fixed-Wing UAVs. Electronics. 2023; 12(18):3946. https://doi.org/10.3390/electronics12183946

Chicago/Turabian Style

Rauniyar, Shyam, and Donghoon Kim. 2023. "Genetic Fuzzy Inference System-Based Three-Dimensional Resolution Algorithm for Collision Avoidance of Fixed-Wing UAVs" Electronics 12, no. 18: 3946. https://doi.org/10.3390/electronics12183946

APA Style

Rauniyar, S., & Kim, D. (2023). Genetic Fuzzy Inference System-Based Three-Dimensional Resolution Algorithm for Collision Avoidance of Fixed-Wing UAVs. Electronics, 12(18), 3946. https://doi.org/10.3390/electronics12183946

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