Research Article  Open Access
Lisong Wang, Xiaoliang Zhang, Pingyu Deng, Jiexiang Kang, Zhongjie Gao, Liang Liu, "An EnergyBalanced Path Planning Algorithm for Multiple Ferrying UAVs Based on GA", International Journal of Aerospace Engineering, vol. 2020, Article ID 3516149, 15 pages, 2020. https://doi.org/10.1155/2020/3516149
An EnergyBalanced Path Planning Algorithm for Multiple Ferrying UAVs Based on GA
Abstract
When performing a search and rescue mission, unmanned aerial vehicles (UAVs) should continuously search targets above the mission area. In order to transfer the search and rescue information quickly and efficiently, two types of UAVs, the ferrying UAVs and the searching UAVs, are used to complete the mission. Obviously, this application scenario requires an efficient path planning method for ferrying UAVs. The existing path planning methods for ferrying UAVs usually focus on shortening the path length and ignore the different initial energy of ferrying UAVs. However, the following problem does exist: if the ferrying UAV with less initial energy is assigned a longer path, meaning that the ferrying UAV with less initial energy will ferry messages for more searching UAVs. When the lowerinitialenergy ferrying UAV is running out of energy, more searching UAVs will no longer deliver messages successfully. Therefore, the mismatch between the planned path length and the initial energy will eventually result in a lower global message delivery ratio. To solve this problem, we propose a new concept energyfactor for a ferrying UAV and use the variance of all ferrying UAVsâ€™ energyfactor to measure the balance between the planned path length and the initial energy. Further, we model the energybalanced pathplanning problem for multiple ferrying UAVs, which actually is a multiobject optimization problem of minimizing the planned path length and minimizing the variance of all ferrying UAVsâ€™ energyfactor. Based on the genetic algorithm, we design and implement an energybalanced path planning algorithm (EMTSPA) for multiple ferrying UAVs to solve this multiobject optimization problem. Experimental results show that EMTSPA effectively increases the global message delivery ratio and decreases the global message delay.
1. Introduction
In recent years, multiUAVs have been used in many domains such as industry, agriculture, military, and disaster relief. Since there are many restrictions on network communication based on the infrastructure, such as terrestrial relay stations and satellites, flying ad hoc networks (FANETs) are generally used to realize communication between multiple UAVs [1â€“3]. Delaytolerant networking (DTN), which is famous of its â€śstoragecarryforwardâ€ť feature, can solve network interruption problems well [4, 5].
In the UAV search and rescue mission, we need to search the mission area continuously. Inspired by DTN mechanism, two types of UAVs, searching UAVs and ferrying UAVs, are used in the mission. Searching UAVs take picture or shot video while ferrying UAVs transfer messages to the ground station [6]. Based on the message routing algorithm, messages are first stored in searching UAVs, and then, messages are transferred to ferrying UAVs at appropriate time; finally, messages are delivered to the ground station.
Intelligent path planning methods and great routing protocols can effectively reduce packet loss and message delay [7]. Regarding the path planning problem of the ferrying node, there have been many studies in recent years. In [8], Message Ferrying (MF) mechanism is first proposed. The paper pointed out that the path planning problem of a ferrying node is a traveling salesman problem (TSP). In addition, the author designed a path for one ferrying node, while ordinary nodes are stationary. In [9], MF mechanism was extended. According to the mobility of ferrying nodes and ordinary nodes, MF mechanism is divided into MF initiated by ordinary nodes and MF initiated by ferrying nodes. Literature [10] designed the route for multiple ferrying nodes in DTN, and it points that there are three ways to interact between ferrying nodes: no interaction, interaction through ferrying nodes, and interaction through ordinary nodes. In [11], moving ordinary nodes are clustered into many clusters according to geographic location, and then, one ferrying node was designed to pass through centers of these clusters along a TSP route. In [12], for the case of moving ordinary nodes, the algorithm firstly selected appropriate points on the paths of the moving ordinary nodes, and then calculated the appropriate encounter time by strengthening the learning algorithm, so that ferrying node could accept ordinary nodesâ€™ messages in time. In [13], ferrying UAVs fly in a pingpong mode along waypoints between the mission area and the ground station. In [14], for mission with noflying zone, paths of ferrying UAVs are designed with determined and online path planning methods. The determined method generates picture points and then generates a route according to picture points, while the online method generates route first, and then generates picture points according to the route. However, since the online path planning method is inefficient in most cases, it is better to adopt the determined path planning method which is performed before flight [14].
Some literatures have further simulated the motion trajectory of a mobile robot. The simulation of the flight trajectory in [15] not only stays in simple kinematic equations but also considers the systemâ€™s dynamic constraints and corresponding dynamic constraints. This study considers the optimization of minimizing motion time and minimizing height and distance, and uses an enhanced particle swarm algorithm to solve the trajectory planning problem. In addition, the path planning algorithm integrates a trajectory tracking algorithm which utilizes the Lyapunovbased constrained backstepping approach and command filters. Reference [16] introduced a novel method to solve the optimal control problem with free initial conditions and verified it in the scenario of minimizing the flight time of the lowthrust orbital transfers. This method uses two evolutionary optimization methods: genetic algorithmparticle swarm optimization and imperial competition algorithm and three orthogonal equations in Hibert space. The above two references both propose corresponding algorithms for the optimal control problem which is to find a control law for a given dynamical system over a period of time. However, this problem in this study only considers the initial energyfactor of the UAV without considering other constraints and control equations of the UAV, which can be regarded as an optimization problem.
In summary, the path planning problem of multiferrying UAVs can be solved by two steps. The first step is to get the set of ferrying UAVsâ€™ waypoints. The second step is to solve the multiple traveling salesman problem (MTSP) according to these waypoints. But the following problem does exist: if the ferrying UAV with less initial energy is assigned a longer path, meaning that the ferrying UAV with less initial energy will ferry messages for more searching UAVs. When the lowerinitialenergy ferrying UAV is running out of energy, more searching UAVs will no longer deliver messages successfully. This unbalanced allocation will result in more searching UAVsâ€™ message being unable to be transmitted to the ground station when the ferrying UAV is energyexhausted, eventually resulting in a lower global message delivery ratio. Based on the above facts, we first put forward the concept of energyfactor, the variance of which reflects the balance between the length of the planning path and the initial energy. Then, we model the energybalanced path planning problem which actually is a multiobjective optimization problem of minimizing the planned path length and minimizing the variance of energyfactor. Finally, an energybalanced path planning algorithm based on the classic genetic algorithm (GA) is proposed to solve this multiobjective optimization problem.
The rest of the paper is organized as follows. In Section 2, we introduce a typical search and rescue scenario and model the network from aspects of mobility models, message routing algorithms, and energy consumption. In Section 3, we model the energybalanced path planning problem for multiferrying UAVs. In Section 4, we introduce the energybalanced path planning algorithm for multiferrying UAVs (EMTSPA) based on GA. In Section 5, we present a comparison of EMTSPA with other path planning algorithms. We conclude in Section 6.
2. Application Scene and Network Model
In this section, we will first introduce a typical search and rescue scene. Next, we will model the flying ad hoc network from aspects such as mobility model, message routing model, and energy consumption.
2.1. The Scenario of Application
This article takes fourrotor drones as an example for discussion. The search and rescue scene is shown in Figure 1. The searching UAVs hover in a fixed area and take pictures. The ferrying UAVs fly in a circular way according to the planned path. In addition, each ferrying UAV ferry messages for at least one searching UAV and messages of one searching UAV should be ferried by only one ferrying UAV. Based on DTN mechanism, messages are first stored in searching UAVs, then transferred to ferrying UAVs at appropriate time, and finally delivered to the ground station.
2.2. Mobility Model
References [17, 18] discuss the mobility model of UAV in threedimensional space, while references [19â€“21] discuss it in a twodimensional plane. In order to facilitate discussion, this paper discusses the mobility model of UAV in twodimensional plane and ignores the collision between drones.
The motion of UAV can be represented by two kinds of mobility models. One is the famous GaussMarkov mobility model, which reflects the correlation between the current state and future state of UAV. The other one reflects the character that UAV flies along the planned path in the search and rescue mission. We call it MapRoute mobility model in this article. (i)GaussMarkov mobility model
The mobility model of UAV is linear and deterministic [13], and the future position of UAV can be estimated by current position, direction, and speed. Therefore, the mobility model of UAV can be represented by the GaussMarkov mobility model [22]. The GaussMarkov equation is shown as follows: where and are the new speed and direction of UAV at interval ; , where , is the tuning parameter used to vary the randomness; and are constants representing the mean value of speed and direction as and and are random variables from a Gaussian distribution. (ii)MapRoute mobility model
When carrying out search and rescue missions, UAVs usually fly in a straight line between the adjacent waypoints of the planned path. This mobility model is called MapRoute mobility model. In this model, the UAV has a target waypoint at any time during the flight. After reaching one target point, UAV stops time () and then continues to move to the next target point.
Suppose the interval time is , and are positions at and . and are speed and direction at . indicates the horizontal axis, and is the vector from the previous waypoint to the next waypoint . represents the angle between the horizontal axis and the vector . So, in the MapRoute mobility model, the position equation of UAV can be expressed as:
2.3. Message Routing Algorithm
In FANET, the communication connection is intermittent due to the fast flight speed of the drone. Therefore, traditional mobile ad hoc network routing algorithms such as AODV, DSR are not suitable for UAV communication [1]. The delay tolerant network (DTN) mechanism is a good solution to communication interruption problems.
Literature [23â€“25] proposed messages routing algorithms based on location information. The main idea of this type of algorithm is as follows: if there is a node in the communication range closer to the destination than the current stored node, messages will be transferred to the closer node. Otherwise, messages will be stored at the current node. However, this type of algorithm has obvious drawbacks. Since nodes may meet for multiple times, the node may receive messages that were just transferred from itself to other nodes. The phenomenon that messages are transferred back and forth is called pingpong effect, which consumes unnecessary energy [13].
In [13], the GPS information is transmitted through a lowthroughput network (IEEE 802.15.4); then, , , and based on geographic location prediction are proposed. Based on the Gaussâ€“Markov mobility model, messages are transferred to the nodes closer to the ground station in the future time. By this way, it decreases the pingpong effect to some extent. However, if the distance between the search area and the ground station exceeds the coverage range of the lowthroughput network, there will be two problems. First, the message routing algorithms cannot completely eliminate the pingpong effect. Second, since GPS data cannot be broadcast between all nodes through the lowthroughput network, the routing algorithm is not feasible in the longdistance search and rescue mission.
Based on the literature [8], we design a servicebased message routing algorithm. The main idea of the routing algorithm is as follows. Messages of searching UAV can only be transmitted to the nearest ferrying UAV or saved in the searching UAV. Messages can be transmitted to the ground station when the ferrying UAVS return back to the ground station. This routing algorithm ensures that messages can be reliably delivered to the ground station even without a lowthroughput network. In addition, this routing algorithm completely eliminates the PingPong effect.
The hop count is the number of hops a message passes until it reaches the destination. This metric allows to reflect the efficiency of the forwarding algorithm. The pingpong effect increases the average number of message hops. Obviously, the servicebased routing algorithm designed in this paper makes the hop count of messages to be 2.
2.4. Energy Consumption Model
Energy is mainly consumed in communication and movement during the flight. Literature [17, 26] pointed out that energy used for communication is much smaller than energy used for movement, so we ignore the communicationrelated energy here.
According to the droneâ€™s MapRoute motion model, we know that energy consumption of UAVs can be divided into two parts: energy consumed in the linear motion and energy consumed in adjusting direction. Let represent the energy consumed in the linear motion phase, while represent the energy consumed in adjusting direction. Then, the total consumed energy can be expressed as follows.
Reference [27] proposed a quadrotor droneâ€™s energy consumption equation in a linear flight at a constant speed. Let denotes the speed, denotes the flight distance, and denotes the speedrelated energy consumption ratio. If is constant, will be unchanged, and energy consumption equation can be expressed as follows.
Generally, in the longdistance mission of UAV, the energy consumed in the linear motion is much larger than the energy consumed in adjusting direction, so we ignore here. According to formulas (4) and (5), if there is a fixed value Î¸ (Î¸â€‰>0), the energy consumption of a drone at a constant speed can be expressed as:
Based on the above analysis, we simply consider that the energy consumption of UAVs in a constantspeed flight is directly proportional to the flight distance.
3. Modeling EnergyBalanced Path Planning Problem for Multiferrying UAVs
Message delivery efficiency and reliability are critical to an UAV search and rescue mission. Therefore, the goal of path planning is not only to shorten the length of path but also to increase the efficiency and reliability of message transmission.
In the following, we first give the evaluation metrics of message delivery, and then model the energybalanced path planning problem for multiferrying UAVs. According to analysis, the problem is actually a multiobjective optimization problem of minimizing length and minimizing the variance of energyfactor. Finally, it is proved that solving this multiobjective optimization problem can increase the message delivery ratio and decrease the message delay under certain conditions.
3.1. Metrics
The message delivery effect is generally evaluated by message delivery ratio, message delay, etc. [13]. However, for a complete search and rescue mission, we insist that it is better to consider the global delivery effect of messages in the network. So, before giving the definition of energybalanced path planning problem, we introduce some important global metrics and its definition firstly.
Definition 1 (global message delivery ratio (GMDR)). The delivery ratio is defined as the fraction of messages that have been successfully delivered to the destination out of the messages that have been generated. It is an indicator for a node. Global message delivery ratio (GMDR) is the average message delivery ratio of all nodes. Global message delivery ratio measures the reliability of message transmission in the global network.
Let denotes the number of searching UAVs, denotes the number of messages generated by searching UAV and finally delivered to the destination successfully, and denotes the number of messages generated by searching UAV . So, GMDR can be expressed as:
Definition 2 (global message delay (GMD)). Message delay is the sum of the communication delay occurred on each hop a message traverses to reach the destination. Global message delay is the average message delay for all messages arriving at the ground station successfully. Global message delay measures the efficiency of message transmission in a global network.
Let denotes the number of messages arriving at ground station successfully, denotes the delivery time of message , denotes the time when message is stored in the searching UAV, and denotes the time when message is relayed by the ferrying UAV. So, GMD can be expressed as:
3.2. Definition of EnergyBalanced Path Planning Problem for Multiferrying UAVs
Before giving the definition of energybalanced path planning problem for multiferrying UAVs, we will give three base definitions firstly: PPL of multiferrying UAVs, EF of ferrying UAV, and variance of energyfactors. Finally, we abstract objective functions and restrictions from the problem.
Definition 3 (planned path length (PPL) of multiferrying UAVs). Planned path length of multiferrying UAVs is the sum of all ferrying UAVâ€™s PPL.
MTSP based on assignment is usually represented with a twoindexes integer linear programming formula [28, 29]. Unlike this, we define PPL of multiferrying UAV with threeindexes.
Suppose there is a graph, where is the set of waypoints of multiferrying UAVs (vertices in the standard MTSP), is the set of all arcs (edges). There is a set of ferrying UAVs (traveling salesmen in the standard MTSP). represents the distance between point s and point t. is a matrix of cost (distance) about set A.
Firstly, we define a binary variable as shown in formula (8). If ferrying UAV is planned to fly from point to point , the value of is 1. Otherwise, the value of is 0.
Then, the PPL of ferrying UAV is
Therefore, PPL of multiferrying UAVs is
Definition 4 (energyfactor (EF) of ferrying UAV). Energyfactor of ferrying UAV is the ratio of a ferrying UAVâ€™s PPL to its initial energy.
Assuming that represents the PPL of the ferrying UAV , represents the initial energy of the ferrying UAV , and then, the EF of the ferrying UAV j can be expressed as:
Definition 5 (variance of energyfactors (VEF)). Variance of energyfactors is the variance of all ferrying UAVsâ€™ EF. It reflects the balance between the PPL and the initial energy of ferrying UAVs in the global network.
Let donates the number of all ferries, donates the average value of all EFs, and donates the VEF. Then, and can be expressed as:
Definition 6 (energybalanced path planning problem for multiferrying UAVs). It is assumed that there is a set N of waypoints (including the ground station). All ferrying UAVs depart from the ground station and return to the ground station. Each ferrying UAV transfers messages for at least one searching UAVs. Messages of each searching UAV are transferred by only one ferrying UAV. The initial energy of ferrying UAVs is different. is the PPL of multiferrying UAVs, and is the VEF. The problem is to find an optimal path for multiferrying UAVs, which can cover all points in N and minimize and .
Obviously, the problem is a multiobjective optimization problem, and the optimization model can be expressed as follows.
Objective function:
Restrictions:
Formulas (16) and (17) indicate that all ferrying UAVs depart from the ground station (point 1 indicates the ground station) and return to the ground station. Formulas (18) and (19) indicate that there is only one ferrying UAV passing through one waypoint. Formula (20) indicates that the ferrying UAV should pass through at least two points and not more than points.
3.3. Basic Properties of the Model
Based on the above definition and abstraction of energybalanced path planning problem, we can find the flowing properties and theorems of the model which provide a basis for designing algorithms and experiments. It should be emphasized that the following properties and theorems are based on the premise that each drone is flying at a constant speed and all dronesâ€™ speed is the same. For the convenience of follow discussion, we first give unified descriptions of the symbols and definitions in Table 1.
Property 7. If , then, .
Proof. According to equation (6), we can further express GMDR with mission time T. Assuming that the number of messages generated by the ferrying UAV during is a constant value (if ), GMDR can be expressed as follows.

Otherwise, if a ferrying UAVâ€™s energy is exhausted, messages of related searching UAVs will no longer be transmitted to the ground station. Therefore, the summation term in equation (22) satisfies:
Then, the maximum GMDR can be expressed as:
Before time , the number of messages generated by the ferrying UAV and delivered successfully to the ground station equals multiply by , and it also equals the number of messages generated by the searching UAV during a lap of the ferrying UAV () multiply by the number of flight laps (). Therefore, the integral term in equation (24) can be converted as follows:
In addition, messages of several searching UAVs may be ferried by one ferrying UAV. Let the number of searching UAVs assigned to the ferrying UAVj be . Then, the maximum GMDR can be further derived from equations (24) and (25).
In summary, Property 7 is proved.
Property 8. If , then, .
Proof. If searching UAVs are stationary, the message delay of message equals the sum of the time that message waits for ferrying UAV and the time that messages are transferred by ferrying UAV [8]. Therefore, the average message delay of the searching UAV () is the time that the ferrying UAV flies for one lap.
If searching UAVs fly in the fixed area, we should consider the delay error caused by the movement of searching UAVs. Therefore, the average message delay of the searching UAV can be expressed as equations (28).
If the PPL of the ferrying UAV () is much larger than the radius (r) of the active area of searching UAVs, will be ignored. Moreover, since one ferrying UAV can ferry messages for multiple searching UAVs, GMD can be expressed as equations (29).
Therefore, Property 8 is proved.
Theorem 9. Suppose there are a PPL set L, an initial energy set E, and a positive number set X, Y. There is a onetoone mapping relationship between L and X, and E and Y. Besides, , , , , , . If Ď‰ () satisfies , then, takes the maximum when is the minimum.
Proof. In general, the more searching UAVs that the ferrying UAV travels, the longer PPL will be. In particular, we assume that the PPL of ferrying UAVs is proportional to the number of searching UAVs which ferrying UAV travels. In other world, we assume there is that satisfies .
It is assumed that PPL of ferrying UAVs constitute vector and initial energy constitute vector . Based on Property 7, we can deduce formula (30).
Since, the smaller is, the greater the balance between the PPL and the initial energy of ferrying UAV is. From formula (13), we know that longer as well as higher or shorter as well as lower can result in better balance. Obviously, when and , is larger than other combinationsâ€™ dot product. In other world, takes the maximum when is the minimum. So, Theorem 9 is proved.
Theorem 10. If , then, there is that satisfies .
Proof. When the number of searching UAVs () whose messages are transferred by the ferrying UAV equals constant (if ), we can reduce following equation (21) from Property (2). So Theorem 10 is proved.
4. EnergyBalanced Path Planning Algorithm for Multiferrying UAVs
Before path planning, we need to get all the waypoints of UAVs. The waypoints can be obtained mainly by clustering, division, and other algorithms [11, 14]. The path planning algorithm of this paper is based on the existing waypoints.
Because solving the shortest path is a part of energybalanced path planning problem for multiferrying UAVs, and the standard MTSP problem is an NPComplete problem. So, energybalanced path planning problem for multiferrying UAVs proposed in this paper is also an NPcomplete problem. Therefore, we can solve it with the heuristic algorithms [30â€“32].
Energybalanced path planning algorithm for multiferrying UAVs (EMTSPA) proposed in this paper is based on the genetic algorithm (GA). Standard genetic algorithms usually consist of these steps: chromosome coding, initialization, selection, evolution, and finally getting the best individual. EMTSPA has made a new design from the fitness function and so on. The flow chart of the algorithm shows in Figure 2.
4.1. Chromosome Coding and Initialization
In the search and rescue mission, each ferrying UAV departs from the ground station. Therefore, the first point of path should be the ground station in the chromosome coding. To better explain the algorithm, we assume that there are 4 ferrying UAVs (ferry 14), 9 waypoints (node 19), and 1 ground station (node 0).
The chromosome coding operation is divided into the following steps. (1) Build a base array of waypoints 19; (2) â€śshuffleâ€ť the base array; (3) generate array â€ślengthsâ€ť which contains the number of waypoints that each ferrying UAV travels (except the ground station); (4) take points from the shuffled base array according to array â€ślengths.â€ť The chromosome coding process is shown in Figure 3.
The above operation generates an individual in the first generation. Initial operation repeats the chromosome coding populationsize times. After initial operation, the firstgeneration population is generated.
4.2. Fitness Function
There are two general methods to solve the multiobjective optimization problem: weightbased method and Paretobased method [33â€“35]. From Theorems 9 and 10, we know that reducing the PPL and the VEF can improve the performance of GMDR and GMD under certain conditions. We solve the multiobjective optimization problem with weightbased methods which combine multiple indicators with weight value. The fitness value is used to measure the individual. So, the fitness function can be expressed as formula (32) where and are weights, is the PPL of multiferrying UAVs, and is the VEF.
We found from formula (32) that the smaller and are, the larger fitness is. Conversely, when or is larger, fitness is smaller. In the verification below, and are set to 10000.
4.3. Genetic Operation
(i)Selection
The purpose of the selection operation is to select the appropriate individuals as parents of the individual in the next generation. EMTSPA uses a twostep selection elite strategy. (1) Copy the individual with the highest fitness from previous generation directly to next generation. (2) Randomly select two subsets from the previous generation with the same size, and then pick out two best individual from two subsets as parents of an individual in the next generation. (ii)Crossover
After selecting the parents, crossover operation is performed to generate a new individual in next generation. We use Figure 4 to illustrate the main process of crossover operation. (1) Pick a continuous segment [1, 4, 5] from array â€śFather.baseâ€ť. (2) Insert the continuous fragment into appropriate place (â€śstartâ€ť to â€śendâ€ť) of array â€śChild.baseâ€ť. (3) Insert other points except [1, 4, 5] of array â€śMother.baseâ€ť into array â€śChild.baseâ€ť orderly; (4) finally, perform chromosome coding based on array â€śChild.baseâ€ť with method described in Section 4.1. (iii)Mutation
To avoid early convergence in the genetic algorithms, individuals need to be mutated. When the generated random value is less than the setting parameter, the individual needs to be mutated. The mutation operation contains the following steps. (1) Select some ferrying UAVs (ferry 2 and ferry 4) from all ferrying UAVs. (2) Select a randomlength segment from base arrays of selected ferrying UAVs and reverse the sequence of segment. As shown in Figure 5, segment [1, 9] is mutated to [1, 9], and segment [3, 6, 7] is mutated to [3, 6, 7]. The mutation operation is shown in Figure 5.
5. Experiments and Analysis
5.1. Setting of Parameters
We experiment in the ONE simulator. The ONE simulator is a popular opportunistic network environment simulator. It can simulate message transmission in different mobility models and different message routing protocol and can provide realtime visual simulation and result output [36, 37].
In search and rescue missions, there are two types of flight modes for ferrying UAVs. (1) Along the path of TSP, multiple ferrying UAVs fly one after another. The following algorithm TSPA belongs to this type. (2) Along the path of MTSP, each ferrying UAV flies cyclically. The following algorithms MTSPA, B_MTSP, and EMTSPA belong to this type.
We mainly compare the performance of the following four algorithms in terms of the PPL of multiferrying UAVs and communication metrics. The four algorithms are as follows: (i)The general algorithm for standard MTSP (MTSPA). It does not consider other factors and selects a random number of waypoints for a ferrying UAV [6](ii)The algorithm for MTSP considering pure balance (BMTSPA). It considers that all ferrying UAVs are the same and aims to make the length of all path completely the same [38](iii)The algorithm for standard TSP (TSPA). It allows multiple ferrying UAVs to fly at equal intervals one after the other [39](iv)The algorithm proposed in this paper for MTSP considering the balance between the initial energy and the PPL of multiferrying UAVs (EMTSPA)
From the perspective of algorithm design, EMTSPA has the following advantages over other algorithms. EMTSPA avoids too long PPL in TSPA; it supports the case where the initial energy of ferrying UAVs is different, while BMTSPA does not support; it avoids complete randomness in MTSPA. Below, we discuss further with the experimental results.
The experimental scene is shown in Figure 6. Suppose there are 4 ferrying UAVs (f10f13), 9 searching UAVs (u1u9), and 1 ground station (g0). Searching UAVs keep flying cyclically along the Zshaped path shown in Figure 6, and centers of the active area are c1c9. Ferrying UAVs depart from the ground station (g0), pass through waypoints (c1c9) along the planned path, and keep flying round by round. Suppose the communication range of the ferrying UAV is not smaller than the range of active area where the searching UAVs flies. In other words, when a ferrying UAV passes through the center of the active area, messages of the searching UAV can be sent to the ferrying UAV. Besides, the ground station remains stationary.
We take searching UAVs flying within a certain range as an example to explain the path planning algorithm. In this article, it is assumed that when ferry UAVs fly through the center of the searching UAVs flight area, ferrying UAVs could receive messages stored in the searching UAVs. Therefore, the set of all points that the multiferrying UAVs need to pass is the all center points of each searching UAVâ€™s flight area. However, the searching UAVs maybe not move in a small range in practical applications. When the searching UAV has a large flight area, we can first divide the area into a number of small subareas using a square with a communication radius of the drone. In this case, the points where the multiferrying UAVs pass are the centers of these subareas. Because this is another category of problems, and in order to concisely grasp the key of the energy equalization algorithm proposed in this paper, we assume the searching UAVs fly in a â€śZâ€ť flight in a square area to explain the algorithm.
Other test parameters are shown in Table 2.

5.2. Results and Analysis
5.2.1. Convergence Behaviour of EMTSPA
Because the above four algorithms with different fitness functions are all based on genetic algorithms, EMTSPA proposed in this study is used as an example to verify the convergence of the genetic algorithm. When the population size is set to 20, one convergence behaviour of fitness function is shown in Figure 7.
5.2.2. Average Planned Path Length
Due to the randomness of the genetic algorithm, in order to compare the performance of the above four algorithms, this paper conducts 50 rounds of tests on the four algorithms to test the average path planning length and energyfactor variance of each algorithm and compares the corresponding average and minimum value.
According to Theorem 10, generally speaking, the smaller the average length, the smaller the global average message delay. Therefore, the average PPL is an important comparison index in path planning. The average PPL of each algorithm in 50 experiments is shown in Figure 8. Moreover, Figure 9 shows the average and minimum values of the average PPL obtained by each algorithm.
It can be seen from Figure 8 that the genetic algorithm has randomness, but the average PPL obtained by the algorithm is stable in a certain range. This is because these four algorithms all aim at least to shorten the distance. Therefore, it is reasonable to compare the performance of each algorithm with the average and minimum values. As shown in Figure 9, the average PPL of EMTSPA, MTSPA, and BMTSPA are similar while the average PPL of TSPA is the largest. Specifically, in terms of the average value of average PPL, EMTSPA is 33.3% better than TSPA. In terms of the minimum value, EMTSPA is 36.7% better than TSPA.
5.2.3. Variance of EnergyFactors
The variance of energyfactors reflects the balance between the PPL and the initial energy of ferrying UAVs. It can be seen from Theorem 9 that the smaller the VEF is, the larger the GMDR is under certain conditions. The experimental results of the energyfactor variance in 50 experiments are shown in Figure 10. The average and minimum values are shown in Figure 11.
It can be seen from Figure 10 that VEFs of TSPA and MTSPA fluctuate greatly, while the VEFs of BMTSPA and EMTSPA are stable. This is because TSPA and MTSPA only consider shortening PPL, EMTSPA considers optimizing VEF, and BMTSPA considers optimizing the variance of the planned length of each aircraft. When PPLs of each aircraft are the same in each round of experiments, VEF must be equal.
A comparison of the VEF is shown in Figure 8. EMTSPA is far superior to other algorithms in the performance of VEF. Since the planned path of each ferrying UAV by TSPA are the same, the VEF of TSPA is much larger than other algorithms when the initial energy is different. Although the VEF of MTSPA is smaller than BMTSPA, MTSPA has a stronger randomness. Specifically, in terms of the average value of VEF, EMTSPA is 99.1% better than TSPA, 97.6% better than MTSPA, and 97.3% better than BMTSPA.
5.2.4. Global Message Delivery Ratio
This study aims to improve the message delivery rate and reduce the message delay by designing a better multiferrying drone path. Therefore, it is necessary to test the planned path in a simulation environment. Considering the randomness of the genetic algorithm, we have run several times to get the path and find the average and minimum values of average PPL and VEF. In this series of experiments, suitable paths whose average PPL and VEF are close to the above results (Figures 9 and 11) are chosen.
This series of experiments attempts to discover the effects of the paths designed by four algorithms on GMDR and GMD before the droneâ€™s energy is exhausted. We have taken six different mission time (200â€‰s, 400â€‰s, 600â€‰s, 800â€‰s, 1000s, 1200â€‰s) to test GMDR. The experimental results of GMDR are shown in Figure 12. (The unit of yaxis is %, and the unit of xaxis is second).
It can be seen from Figure 12: (i)GMDR of EMTSPA is slightly lower than GMDR of BMTSPA when test time is short. We are convinced as follows. In EMTSPA, the ferrying UAV with longer PPL is responsible for more messageferrying task. Due to the short test time, the ferrying UAV with longer PPL fails to deliver messages to the ground station in time. It results that GMDR of EMTSPA is slightly lower than GMDR of BMTSPA. However, GMDR of EMTSPA is still greater than GMDR of TSPA and MTSPA(ii)Most of the time, GMDR of EMTSPA is greater than other algorithms. As test time increases, some ferrying UAVs may be energyexhausted, then, these energyexhausted ferrying UAVs can no longer ferry messages. While in EMTSPA, the ferrying UAVs with more energy are assigned longer PPL. In other words, ferrying UAVs with more energy can ferry more messages and ferry messages in longer time. Undoubtedly, this will increase GMDR. The experimental results have verified this(iii)From the overall performance of GMDR, EMTSPA performed best, BMTSPA performed second, and MTSPA and TSPA performed poorly. The comparison is shown as Table 3. (The â€ś+â€ť means that EMTSPAâ€™s performance is better while the â€śâ€ť means that EMTSPAâ€™s performance is worse. The default unit is %. For example, +925.9 in Table 3 means that EMTSPAâ€™s GMDR is 925.9% of TSPAâ€™s GMDR)

5.2.5. Global Message Delay
When test time is different, the result of the global message delay in different test time is shown in Figure 13. As can be seen from Figure 13, GMD of TSPA is much larger than other algorithms, and GMD of EMTSPA is close to BMTPSAâ€™s. The PPL of MTSPA is the smallest, and GMD of MTSPA is the smallest. It can be found that GMD has a positive correlation with the PPL. The results also verify the correctness of Theorem 10 in Section 3. Specifically, EMTSPAâ€™s GMD is only 50% of TSPAâ€™s.
6. Conclusions
In this paper, we considered the path planning problem in search and rescue scenarios using UAV networks. In particular, we have noticed such a problem that the existing path planning methods for ferrying UAVs usually focus on shortening the path length and ignore the different initial energy of ferrying UAVs, and with ignoring this problem, the next problem will arise: lowenergy ferry UAV is assigned a longer path, and this lowenergy ferry UAV may be no longer deliver messages successfully. To solve the mismatch between the planned path length and the initial energy, we propose a new concept named energyfactor for the ferrying UAV, and use the variance of all ferrying UAVsâ€™ energyfactor to measure the balance between the planned path length and the initial energy. Based on multiobject optimization method, we propose a novel method named EMTSPA to deal with this problem. We have completed the modeling, algorithm design, analysis, and experimental verification; the results show that EMTSPA has a significant improvement in global message delivery ratio and global message delay compared to other algorithms. However, because the algorithm in this study uses weight coefficients to superimpose two optimization goals, the choice of weight coefficients is subjective. The algorithm can be further optimized based on Pareto solution which finally get a set of solutions.
Data Availability
If need any information about the paper, please contact the authors.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
References
 Ä°. Bekmezci, O. K. Sahingoz, and Ĺž. Temel, â€śFlying adhoc networks (FANETs): a survey,â€ť Ad Hoc Networks, vol. 11, no. 3, pp. 1254â€“1270, 2013. View at: Publisher Site  Google Scholar
 S. Hayat, E. Yanmaz, and R. Muzaffar, â€śSurvey on unmanned aerial vehicle networks for civil applications: a communications viewpoint,â€ť IEEE Communications Surveys & Tutorials, vol. 18, no. 4, pp. 2624â€“2661, 2016. View at: Publisher Site  Google Scholar
 E. W. Frew and T. X. Brown, â€śAirborne communication networks for small unmanned aircraft systems,â€ť Proceedings of the IEEE, vol. 96, no. 12, pp. 2008â€“2027, 2008. View at: Publisher Site  Google Scholar
 K. Fall, â€śA delaytolerant network architecture for challenged internets,â€ť in Proceedings of the 2003 conference on Applications, technologies, architectures, and protocols for computer communications  SIGCOMM '03, pp. 27â€“34, Karlsruhe Germany, August 2003. View at: Google Scholar
 S. Jain, K. Fall, and R. Patra, â€śRouting in a delay tolerant network,â€ť ACM SIGCOMM Computer Communication Review, vol. 34, no. 4, p. 145, 2004. View at: Publisher Site  Google Scholar
 M. M. B. Tariq, M. Ammar, and E. Zegura, â€śMessage ferry route design for sparse ad hoc networks with mobile nodes,â€ť in ACM International Symposium on Mobile Ad Hoc NETWORKING and Computing, pp. 37â€“48, Florence Italy, 2006. View at: Google Scholar
 V. Soundaryah and D. Shanthi, â€śEfficient path planning based on hybrid routing protocols over a VANET,â€ť International Journal of Applied Engineering Research, vol. 10, no. 20, pp. 16654â€“16659, 2017. View at: Google Scholar
 W. Zhao and M. H. Ammar, â€śMessage ferrying: proactive routing in highlypartitioned wireless ad hoc networks,â€ť in The Ninth IEEE Workshop on Future Trends of Distributed Computing Systems, 2003. FTDCS 2003. Proceedings, pp. 308â€“314, San Juan, Puerto Rico, USA, 2003. View at: Google Scholar
 W. Zhao, â€śA message ferrying approach for data delivery in sparse mobile ad hoc networks,â€ť in Proceedings of the 5th ACM international symposium on Mobile ad hoc networking and computing  MobiHoc '04, pp. 187â€“198, Roppongi Hills, Tokyo Japan, 2004. View at: Google Scholar
 W. Zhao, M. Ammar, and E. Zegura, â€śControlling the mobility of multiple data transport ferries in a delaytolerant network,â€ť in Proceedings IEEE 24th Annual Joint Conference of the IEEE Computer and Communications Societies, vol. 2, pp. 1407â€“1418, Miami, FL, USA, 2005. View at: Google Scholar
 H. Miura, D. Nishi, N. Matsuda, and H. Taki, â€śMessage ferry route design based on clustering for sparse ad hoc networks,â€ť in KnowledgeBased and Intelligent Information and Engineering Systems, pp. 637â€“644, SpringerVerlag, 2010. View at: Google Scholar
 K. Scott, R. Dai, and J. Zhang, â€śOnlinerelayingbased image communication in unmanned aerial vehicle networks,â€ť in 2017 IEEE International Conference on Communications (ICC), pp. 1â€“6, Paris, France, 2017. View at: Google Scholar
 M. Asadpour, K. A. Hummel, D. Giustiniano, and S. Draskovic, â€śRoute or carry: motiondriven packet forwarding in micro aerial vehicle networks,â€ť IEEE Transactions on Mobile Computing, vol. 16, no. 3, pp. 843â€“856, 2017. View at: Publisher Site  Google Scholar
 E. Yanmaz, R. Kuschnig, M. Quaritsch, C. Bettstetter, and B. Rinner, â€śOn path planning stratiogies for networked unmanned aerial vehicles,â€ť in 2011 IEEE Conference on Computer Communications Workshops (INFOCOM WKSHPS), pp. 212â€“216, Shanghai, China, 2011. View at: Google Scholar
 J. Karimi and S. H. Pourtakdoust, â€śIntegrated motion planning and trajectory control system for unmanned air vehicles,â€ť Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 227, no. 1, pp. 3â€“18, 2012. View at: Publisher Site  Google Scholar
 I. Shafieenejad, A. B. Novinzadeh, and V. R. Molazadeh, â€śIntroducing a novel algorithm for minimumtime lowthrust orbital transfers with free initial condition,â€ť Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 229, no. 2, pp. 333â€“351, 2014. View at: Publisher Site  Google Scholar
 Y. Zeng and R. Zhang, â€śEnergyefficient UAV communication with trajectory optimization,â€ť IEEE Transactions on Wireless Communications, vol. 16, no. 6, pp. 3747â€“3760, 2017. View at: Publisher Site  Google Scholar
 Y. Zeng, J. Xu, and R. Zhang, â€śEnergy minimization for wireless communication with rotarywing UAV,â€ť IEEE Transactions on Wireless Communications, vol. 18, no. 4, pp. 2329â€“2345, 2019. View at: Publisher Site  Google Scholar
 A. Konak, G. E. Buchert, and J. Juro, â€śA flockingbased approach to maintain connectivity in mobile wireless ad hoc networks,â€ť Applied Soft Computing, vol. 13, no. 2, pp. 1284â€“1291, 2013. View at: Publisher Site  Google Scholar
 Q. Wu, Y. Zeng, and R. Zhang, â€śJoint trajectory and communication design for multiUAV enabled wireless networks,â€ť IEEE Transactions on Wireless Communications, vol. 17, no. 3, pp. 2109â€“2121, 2018. View at: Publisher Site  Google Scholar
 Y. Zeng, X. Xu, and R. Zhang, â€śTrajectory design for completion time minimization in UAVenabled multicasting,â€ť IEEE Transactions on Wireless Communications, vol. 17, no. 4, pp. 2233â€“2246, 2018. View at: Publisher Site  Google Scholar
 T. Camp, J. Boleng, and V. Davies, â€śA survey of mobility models for ad hoc network research,â€ť Wireless Communications & Mobile Computing, vol. 2, no. 5, pp. 483â€“502, 2010. View at: Google Scholar
 E. Kuiper and S. NadjmTehrani, â€śGeographical routing in intermittently connected ad hoc networks,â€ť in 22nd International Conference on Advanced Information Networking and Applications  Workshops (aina workshops 2008), pp. 1690â€“1695, Okinawa, Japan, 2008. View at: Google Scholar
 R. Shirani, M. StHilaire, T. Kunz, Y. Zhou, J. Li, and L. Lamont, â€śThe performance of greedy geographic forwarding in unmanned aeronautical adhoc networks,â€ť in 2011 Ninth Annual Communication Networks and Services Research Conference, pp. 161â€“166, Ottawa, ON, Canada, 2011. View at: Google Scholar
 M. Y. Arafat and S. Moh, â€śLocationaided delay tolerant routing protocol in UAV networks for postdisaster operation,â€ť IEEE Access, vol. 6, pp. 59891â€“59906, 2018. View at: Publisher Site  Google Scholar
 C. Desset, B. Debaillie, V. Giannini et al., â€śFlexible power modeling of LTE base stations,â€ť in 2012 IEEE Wireless Communications and Networking Conference (WCNC), pp. 2858â€“2862, Paris, France, April 2012. View at: Google Scholar
 C. Di Franco and G. Buttazzo, â€śEnergyaware coverage path planning of UAVs,â€ť in 2015 IEEE International Conference on Autonomous Robot Systems and Competitions, pp. 111â€“117, Vila Real, Portugal, 2015. View at: Google Scholar
 T. Bektas, â€śThe multiple traveling salesman problem: an overview of formulations and solution procedures,â€ť Omega, vol. 34, no. 3, pp. 209â€“219, 2006. View at: Publisher Site  Google Scholar
 A. KirĂˇly and J. Abonyi, â€śA Novel Approach to Solve Multiple Traveling Salesmen Problem by Genetic Algorithm,â€ť in Computational Intelligence in Engineering, pp. 141â€“151, Springer, Berlin, Heidelberg, 2010. View at: Google Scholar
 A. Singh, â€śA Review on Algorithms Used to Solve Multiple Travelling Salesman Problem,â€ť vol. 3, pp. 598â€“603, 2016. View at: Google Scholar
 J. Gupta and S. Lal, â€śComparison of some Approximate Algorithms Proposed for Traveling Salesmen and Graph Partitioning Problems,â€ť in Proceedings of 3rd International Conference on Internet of Things and Connected Technologies (ICIoTCT), Jaipur, India, March 2018. View at: Google Scholar
 J. Henrio, T. Deligne, T. Nakashima, and T. Watanabe, â€śRoute planning for multiple surveillance autonomous drones using a discrete firefly algorithm and a Bayesian optimization method,â€ť Artificial Life and Robotics, vol. 24, no. 1, pp. 100â€“105, 2019. View at: Publisher Site  Google Scholar
 A. Konak, D. W. Coit, and A. E. Smith, â€śMultiobjective optimization using genetic algorithms: a tutorial,â€ť Reliability Engineering & System Safety, vol. 91, no. 9, pp. 992â€“1007, 2006. View at: Publisher Site  Google Scholar
 J. Li, H. Sang, Y. Han, C. G. Wang, and K. Z. Gao, â€śEfficient multiobjective optimization algorithm for hybrid flow shop scheduling problems with setup energy consumptions,â€ť Journal of Cleaner Production, vol. 181, pp. 584â€“598, 2018. View at: Publisher Site  Google Scholar
 S. ZapotecasMartinez, C. A. Coello Coello, H. E. Aguirre, and K. Tanaka, â€śA review of features and limitations of existing scalable multiobjective test suites,â€ť IEEE Transactions on Evolutionary Computation, vol. 23, no. 1, pp. 130â€“142, 2019. View at: Publisher Site  Google Scholar
 A. KerĂ¤nen, J. Ott, and T. KĂ¤rkkĂ¤inen, â€śThe one simulator for DTN protocol evaluation,â€ť in Proceedings of the 2nd international conference on simulation tools and techniques. ICST (Institute for Computer Sciences, SocialInformatics and Telecommunications Engineering), p. 55, Rome Italy, 2009. View at: Google Scholar
 M. S. Desta, E. HyytiĂ¤, A. KerĂ¤nen, T. KĂ¤rkkĂ¤inen, and J. Ott, â€śEvaluating (Geo) content sharing with the ONE simulator,â€ť in Proceedings of the 11th ACM international symposium on Mobility management and wireless access  MobiWac '13, pp. 37â€“40, Barcelona Spain, 2013. View at: Google Scholar
 M. Zhang and W. Cai, â€śUnderwater targets tracking path planning based on task cooperation of multiple AUVs,â€ť Chinese Journal of Sensors and Actuators, vol. 31, no. 7, pp. 1101â€“1107, 2018. View at: Publisher Site  Google Scholar
 B. Fraser, R. Hunjet, and A. Coyle, â€śA Swarm Intelligent Approach to Data Ferrying in Sparse Disconnected NetworksA swarm intelligent approach to data ferrying in sparse disconnected networks,â€ť International Journal of Parallel, Emergent and Distributed Systems, pp. 1â€“22, 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Lisong Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.