# Impact of UAV Delivery on Sustainability and Costs under Traffic Restrictions

Authors: Yuyu Li; Wei Yang; Bo Huang
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

Compared with traditional vehicles delivery, unmanned aerial vehicle (UAV) delivery can reduce energy consumption and greenhouse gas emissions, which benefits environmental sustainability. Besides, UAVs can overcome traffic restrictions, which are the big obstacle in parcel delivery. In reality, there are two kinds of most popular traffic restrictions, vehicle-type restriction, and half-side traffic. We propose a mixed-integer (0-1 linear) green routing model with these two kinds of traffic restrictions for UAVs to exploit the environmental aspects of the use of UAVs in logistics. A genetic algorithm is proposed to efficiently solve the complex routing problem, and an experimental analysis is made to illustrate and validate our model and the algorithm. We found that, under both these two traffic restrictions, UAV delivery can accomplish deliveries that cannot be carried out or are carried out at much higher costs by vehicles only and can always effectively save costs and cut CO2 emissions, which is environmentally friendly. Furthermore, UAV delivery saves more cost and cuts more CO2 emission under the first kind of traffic restriction than that under the second.


---

## Body

## 1. Introduction

Brick-and-mortar retail stores being gradually replaced by e-commerce saves energy consumption for heating and lighting in stores and warehouses. However, the rapid development of e-commerce dramatically increases the demands on logistics or parcel delivery, which is mainly carried out by traditional vehicles. As a result, energy consumption and CO2 emission are increased [1].With the development of unmanned aerial vehicles (UAVs), they become one of the best choices improving the “last-mile” delivery of products to consumers, both economically and environmentally, due to their significant performances in saving costs and energy. Unfortunately, current UAVs have two big shortcomings: both delivery range (distance and flight time) and capacity (weight and size) are limited. Therefore, they often cannot deliver all packages themselves in a single trip, and the most reasonable way to use a UAV is to pair it with a traditional vehicle. In this “UVA-vehicle mode,” a UAV is loaded on the roof of the vehicle, and the driver guides the UAV to carry the package to the destination. At the same time of UAV delivery, the vehicle serves other customers. Although the UAV only serves part of customers, this mode of package delivery makes a best use of the advantages of UAVs and avoids their disadvantages. In other words, this mode normally delivers all packages in one signal trip with much lower time, cost, and energy consumption.Besides, in many big cities, because of heavy traffic jam or severe air pollution from vehicle exhaust, vehicle-type restriction is implemented on vehicles. Delivery vehicles are forbidden to run on some roads during some periods. Besides, during the maintenance of some roads, only one lane is available, and vehicles traveling in opposite directions have to take turns on the road with a much lower speed. This is half-side traffic, which is other most popular traffic restriction. Traffic restrictions bring big troubles to package delivery service and even some “impossible missions” to delivery by vehicles only, which can be normally overcome through UAV delivery.In order to make the most use of the advantages of this “UAV-vehicle mode” and better understand the economic and environmental impact of using UAVs on packages delivery, we should optimize the routing plan of “UVA-vehicle mode.” To this end, we propose a mixed-integer (0-1 linear) green routing model with traffic restrictions for UAVs and design a genetic algorithm (GA) to effectively solve the routing plan model. Finally, we make an experimental analysis to illustrate and validate our model and algorithm, as well as better understand the environmental benefits of delivering with UAVs.The rest of this paper is organized as follows: we review the literature on “green routing” and routing problem with UAVs in Section2; we propose a UAV-vehicle routing model with traffic restrictions and design a genetic algorithm to effectively solve the model in Sections 3 and 4, respectively; an experimental analysis is made in Section 5; and conclusions are drawn in Section 6.


## 2. Literature Review

Since the emergence of “truck dispatching problem,” vehicle routing problem (VRP) has been a hot research field [2]. And many heuristic algorithms were developed to solve VRP [3]. Our paper is mainly relevant to two topics: green vehicle routing problem and routing with UAVs, which are reviewed below.


### 2.1. Green Vehicle Routing Problem

A series of studies show that, for the logistics, transportation accounts for about 90% of the energy consumption of the whole logistics. Therefore, reducing CO2 emissions is the core of the green vehicle routing problem (GVRP). As early as 1985, Cermak and Takeda [4] outlined criteria for simulating atmospheric boundary layer and physical simulation of source characteristics in boundary layer wind tunnels and studied the air problem in urban environment. Their work led subsequent scholars to study logistics and carbon emissions.Pradenas et al. [5] studied the vehicle path problem for the energy required for each path from the perspective of vehicle distance and estimated the load and distance between customers to achieve the goal of reducing fuel consumption and carbon emission. Qian [6] developed routing and scheduling model for fleets of transport vehicles to minimize fuel emissions in the speed-dependent road network from the perspective of vehicle speed. In this paper, the route of each vehicle needs to be clarified, the speed of vehicles on different roads on their respective paths is regarded as a decision variable, and the time insertion algorithm is given for a single path. Finally, a tabu search algorithm based on column generation is proposed to solve this problem. Kwon et al. [7] built a vehicle path model that minimizes energy consumption and pollution emissions with time windows from the perspective of vehicle load. Their simulation process is based on the actual route of motor vehicles, and their method can save more than 6.9% of fuel compared with the existing method. Although these papers only studied GVRP without UAVs, we can refer to their methods of calculating delivery cost and fuel consumption in establishing our model.


### 2.2. Routing with UAVs

The “last kilometer problem” is always the bottleneck of logistics distribution due to traffic conditions, distribution personnel, customer acceptance location, and other factors. In recent years, UAVs were introduced into logistics to reduce costs and solve the “last kilometer problem.”In 2013, Jeff Bezos, CEO and founder ofAmazon.com, announced on “60 minutes” that UAVs could be used to speed up the delivery of packages to consumers. Subsequently, routing with UAVs became a hot topic in VRP field [8]. D’Andrea [9] calculated the energy consumption and the cost of high-end lithium ion batteries in the high-performance mobile delivery of UAVs based on the principle of first computability and systematically analyzed the economy of UAV technology and future applications. Zhi [10] designed a two-stage hitch-matching algorithm to serve the UAV driving equipment, which greatly reduced the logistics cost and improved the efficiency. Weng [11] made a systematic evaluation of UAV logistics from the aspects of laws and regulations, technical safety, audience preference, operation safety, and use cost. Ma et al. [12] established a model of UAV flight stability index and various factors and analyzed and studied the main factors influencing the flight stability of UAVs through Matlab.The feature of our “UVA-vehicle mode” is that UAV departures from the roof of the vehicle to deliver package; meanwhile, the vehicle serves other customers, which improves logistics efficiency and reduces energy consumption. Chiang et al. [2] proposed a UAV-vehicle model and a GA to study the economic and environmental impact of using UAVs on package delivery. Murray and Chu [13] conducted relevant research, proposing two mixed-integer linear programming formulas for unmanned delivery problems, as well as two simple but effective heuristic solutions to satisfy the scheduling of UAVs and delivery trucks. Ha et al. [14] considered a new variant of TSP-D in which the objective is to minimize operational costs including total transportation cost and one created by waste time a vehicle has to wait for the other. They formulated problem and proposed two algorithms, TSP-LS and GRASP, to solve the problem. Wang et al. [15] pose a number of questions in order to study the maximum savings that can be obtained from using drones and derived a number of worst-case results. Poikonen et al. [16] studied the UAV-vehicle delivery routing with the aim of minimizing the total operation time. Pugliese and Guerriero [17] analyzed the delivery process with drones, by taking into account the total transportation cost, under the assumption that all customers should be served within their time window.Although the above papers studied UAV-vehicle routing problems under different circumstances, they failed to research the use of UAVs and its effect on relieving one of the greatest constraints in parcel delivery and traffic restrictions [18]. In fact, there are two kinds of traffic restrictions, which allow the implementation of our UAV-vehicle mode but greatly affect the cost, energy consumption, and efficiency of parcel delivery. The first is vehicle-type restriction, mainly on delivery vehicles. For example, during the maintenance, all vehicles except public transportation were banned to run across Shimen Bridge, located in Chongqing, China. Under this circumstance, the delivery vehicle has to choose other way or use a UAV. The second is half-side traffic, under which the delivery vehicle can choose this road but may have to wait to pass through it. During its waiting, the engine is still operating and consuming gasoline, causing a higher cost and CO2 emission [19].


### 2.3. Algorithms for VRP

To solve VRP, we should design an algorithm. Generally, there are two categories of algorithms: exact algorithm and approximation algorithm [20]. The advantage of exact algorithm is that it can find the accurately optimal solution of VRP, but it is only suitable for VRP of small scale [21]. When the scale is large, it is mainly applied to obtain the initial solution for the approximation algorithm [22].Approximation algorithm can be roughly divided into two categories: heuristic algorithm, which consists of constructive heuristic algorithm and improved heuristic algorithm, as well as meta-heuristic algorithm [23]. The advantage of constructive heuristic algorithm is that it is simple and easy to understand, but the solutions found may be far from the optimal solutions. Therefore, it is no longer used for solving VRP alone but combined with the improved heuristic algorithm, being used to generate initial solutions [24]. Improved heuristic algorithm can obtain better solutions from initial solutions generated by constructive heuristic algorithm through neighborhood search. Its advantage is that the probability of obtaining the optimal solution is high, but the operation time may be very long [25].The biggest disadvantage of heuristic algorithm is that it is easy to fall into local optimal. In order to overcome it, there appeared a variety of meta-heuristic algorithms with the advantage of jumping out of the local optimal and seek the global optimal, such as genetic algorithm, simulated annealing algorithm, and tabu search algorithm [26, 27]. Among meta-heuristic algorithms, the genetic algorithm is popularly applied for solving VRP, including VRP with UAVs [28, 29]. However, the mode of UAV delivery in these studies is that the truck is just used as a mobile base, which is very different with ours. Therefore, these genetic algorithms cannot solve the key issue in our mode, that is, which customer is served by the UAV and how many and which customers should be served by the vehicle.The main contributions of this paper are as follows. Firstly, taking into account two kinds of traffic restrictions which affect parcel delivery, we propose a mixed-integer (0-1 linear) UAV-vehicle routing model with traffic restrictions to incorporate environment aspects to study the impacts of using UAVs for package delivery under these two circumstances and analyze the efficiency of UAV delivery to reduce cost and energy consumption by comparing the results under these two circumstances and that without UAVs. Secondly, we develop a genetic algorithm to effectively solve the model and analyze the impact of UAV delivery. In our GA, we make improvement and contribution in the generation of initial solution and grouping customers to determine which customer is served by the UVA and which customer(s) should be served by the vehicle during UAV delivery.


## 2.1. Green Vehicle Routing Problem

A series of studies show that, for the logistics, transportation accounts for about 90% of the energy consumption of the whole logistics. Therefore, reducing CO2 emissions is the core of the green vehicle routing problem (GVRP). As early as 1985, Cermak and Takeda [4] outlined criteria for simulating atmospheric boundary layer and physical simulation of source characteristics in boundary layer wind tunnels and studied the air problem in urban environment. Their work led subsequent scholars to study logistics and carbon emissions.Pradenas et al. [5] studied the vehicle path problem for the energy required for each path from the perspective of vehicle distance and estimated the load and distance between customers to achieve the goal of reducing fuel consumption and carbon emission. Qian [6] developed routing and scheduling model for fleets of transport vehicles to minimize fuel emissions in the speed-dependent road network from the perspective of vehicle speed. In this paper, the route of each vehicle needs to be clarified, the speed of vehicles on different roads on their respective paths is regarded as a decision variable, and the time insertion algorithm is given for a single path. Finally, a tabu search algorithm based on column generation is proposed to solve this problem. Kwon et al. [7] built a vehicle path model that minimizes energy consumption and pollution emissions with time windows from the perspective of vehicle load. Their simulation process is based on the actual route of motor vehicles, and their method can save more than 6.9% of fuel compared with the existing method. Although these papers only studied GVRP without UAVs, we can refer to their methods of calculating delivery cost and fuel consumption in establishing our model.


## 2.2. Routing with UAVs

The “last kilometer problem” is always the bottleneck of logistics distribution due to traffic conditions, distribution personnel, customer acceptance location, and other factors. In recent years, UAVs were introduced into logistics to reduce costs and solve the “last kilometer problem.”In 2013, Jeff Bezos, CEO and founder ofAmazon.com, announced on “60 minutes” that UAVs could be used to speed up the delivery of packages to consumers. Subsequently, routing with UAVs became a hot topic in VRP field [8]. D’Andrea [9] calculated the energy consumption and the cost of high-end lithium ion batteries in the high-performance mobile delivery of UAVs based on the principle of first computability and systematically analyzed the economy of UAV technology and future applications. Zhi [10] designed a two-stage hitch-matching algorithm to serve the UAV driving equipment, which greatly reduced the logistics cost and improved the efficiency. Weng [11] made a systematic evaluation of UAV logistics from the aspects of laws and regulations, technical safety, audience preference, operation safety, and use cost. Ma et al. [12] established a model of UAV flight stability index and various factors and analyzed and studied the main factors influencing the flight stability of UAVs through Matlab.The feature of our “UVA-vehicle mode” is that UAV departures from the roof of the vehicle to deliver package; meanwhile, the vehicle serves other customers, which improves logistics efficiency and reduces energy consumption. Chiang et al. [2] proposed a UAV-vehicle model and a GA to study the economic and environmental impact of using UAVs on package delivery. Murray and Chu [13] conducted relevant research, proposing two mixed-integer linear programming formulas for unmanned delivery problems, as well as two simple but effective heuristic solutions to satisfy the scheduling of UAVs and delivery trucks. Ha et al. [14] considered a new variant of TSP-D in which the objective is to minimize operational costs including total transportation cost and one created by waste time a vehicle has to wait for the other. They formulated problem and proposed two algorithms, TSP-LS and GRASP, to solve the problem. Wang et al. [15] pose a number of questions in order to study the maximum savings that can be obtained from using drones and derived a number of worst-case results. Poikonen et al. [16] studied the UAV-vehicle delivery routing with the aim of minimizing the total operation time. Pugliese and Guerriero [17] analyzed the delivery process with drones, by taking into account the total transportation cost, under the assumption that all customers should be served within their time window.Although the above papers studied UAV-vehicle routing problems under different circumstances, they failed to research the use of UAVs and its effect on relieving one of the greatest constraints in parcel delivery and traffic restrictions [18]. In fact, there are two kinds of traffic restrictions, which allow the implementation of our UAV-vehicle mode but greatly affect the cost, energy consumption, and efficiency of parcel delivery. The first is vehicle-type restriction, mainly on delivery vehicles. For example, during the maintenance, all vehicles except public transportation were banned to run across Shimen Bridge, located in Chongqing, China. Under this circumstance, the delivery vehicle has to choose other way or use a UAV. The second is half-side traffic, under which the delivery vehicle can choose this road but may have to wait to pass through it. During its waiting, the engine is still operating and consuming gasoline, causing a higher cost and CO2 emission [19].


## 2.3. Algorithms for VRP

To solve VRP, we should design an algorithm. Generally, there are two categories of algorithms: exact algorithm and approximation algorithm [20]. The advantage of exact algorithm is that it can find the accurately optimal solution of VRP, but it is only suitable for VRP of small scale [21]. When the scale is large, it is mainly applied to obtain the initial solution for the approximation algorithm [22].Approximation algorithm can be roughly divided into two categories: heuristic algorithm, which consists of constructive heuristic algorithm and improved heuristic algorithm, as well as meta-heuristic algorithm [23]. The advantage of constructive heuristic algorithm is that it is simple and easy to understand, but the solutions found may be far from the optimal solutions. Therefore, it is no longer used for solving VRP alone but combined with the improved heuristic algorithm, being used to generate initial solutions [24]. Improved heuristic algorithm can obtain better solutions from initial solutions generated by constructive heuristic algorithm through neighborhood search. Its advantage is that the probability of obtaining the optimal solution is high, but the operation time may be very long [25].The biggest disadvantage of heuristic algorithm is that it is easy to fall into local optimal. In order to overcome it, there appeared a variety of meta-heuristic algorithms with the advantage of jumping out of the local optimal and seek the global optimal, such as genetic algorithm, simulated annealing algorithm, and tabu search algorithm [26, 27]. Among meta-heuristic algorithms, the genetic algorithm is popularly applied for solving VRP, including VRP with UAVs [28, 29]. However, the mode of UAV delivery in these studies is that the truck is just used as a mobile base, which is very different with ours. Therefore, these genetic algorithms cannot solve the key issue in our mode, that is, which customer is served by the UAV and how many and which customers should be served by the vehicle.The main contributions of this paper are as follows. Firstly, taking into account two kinds of traffic restrictions which affect parcel delivery, we propose a mixed-integer (0-1 linear) UAV-vehicle routing model with traffic restrictions to incorporate environment aspects to study the impacts of using UAVs for package delivery under these two circumstances and analyze the efficiency of UAV delivery to reduce cost and energy consumption by comparing the results under these two circumstances and that without UAVs. Secondly, we develop a genetic algorithm to effectively solve the model and analyze the impact of UAV delivery. In our GA, we make improvement and contribution in the generation of initial solution and grouping customers to determine which customer is served by the UVA and which customer(s) should be served by the vehicle during UAV delivery.


## 3. Problem Description and the Model

### 3.1. Problem Description

In an urban area, a delivery company uses vehicles equipped with UAVs to deliver parcels from a warehouse to customers in its delivery area. Each customer has only one parcel waiting to be delivered. This UVA-vehicle mode is as follows. Every vehicle is equipped with a UAV. As both delivery range (distance and flight time) and capacity (weight and size) of the UAV are limited, it only delivers parts of the parcels, but the vehicle can deliver all the parcels. Each UAV may depart from its vehicle at a location of a customer or the warehouse and carry a parcel to one and only one customer. Then, the UAV returns to its vehicle to reload a parcel and recharge or swap batteries, which is instantaneous. While the UAV delivers the parcel, the vehicle carries out its delivery to one or several customers. Therefore, the UAV returns to its vehicle at a different customer location. As a matter of course, if the vehicle arrives at the customer location, where the UAV is retrieved, it waits for the UAV, and vice versa. For the sake of the safety of the UAV, if the UAV arrives before the vehicle, it hovers in the air waiting for the vehicle.In the delivery area of the warehouse, there may be traffic restrictions on the vehicles. Generally, there are two kinds of traffic restrictions, under which this UAV-vehicle mode is practicable but greatly affected. The first is vehicle-type restriction, under which delivery vehicles are forbidden to run through some paths and have to choose other way or use a UAV. The second is half-side traffic, under which only one lane of a road is available, and vehicles traveling in opposite directions can only take turns on this road. Under this circumstance, the delivery vehicle can run on this road but may have to wait for its turn with the operating engine and raise the variable cost and CO2 emission.According to the regulation on working hours, the service time of every pair of UAV and vehicle is 8 hours. Therefore, it is probable that more than one pair of UAV and vehicle are needed.There are two goals of routing the UAV-vehicle. The first is minimizing the total cost, and the alternative is minimizing the energy consumption (or CO2 emission). The total cost consists of fixed cost whenever a pair of UAV-vehicle is used, as well as variable cost, which is the function of unit route cost, travel distance, and gross weight (empty weight of the vehicle or UAV plus the payload). The total CO2 emission of a vehicle is the function of weighted average emission rate of vehicles, travel distance, and gross weight, and that of a UAV is the function of CO2 emission rate of generating per watt-hour, average energy requirement of UAV, travel distance, and gross weight.It should be noted that the two kinds of traffic restrictions may simultaneously exist, but as we analyze the impact of UAV delivery under every kind of traffic restriction, they are in our models separately.


### 3.2. The Model

According to Laporte [30], there are three ways to formulate the vehicle routing problem (VRP), including the simple set division formula of VRP first proposed by Balinski and Quandt in 1964, the commodity flow formula of Shlifer and Graves in 1979, and the two-index vehicle flow formula of Laporte and Norbert in 1983. In our study, the GVRP problem is mainly involved, and the weight of vehicles, commodities, and drones in the whole route needs to be considered. Therefore, we choose to continue to expand the commodity flow formula, which is clearly expressed and promotes the development of heuristic algorithm.The basic GVRP problem is as follows: both the UAV and the vehicle can travel back and forth on any available path. Therefore, an undirected graphR = (M, N) is given, where N is the set of edges with nonnegative routing costs, {i, j} ∈ N represents the edges from node i to node j. M = {0, 1, …, n} is the set of all nodes with a total of n customers, and the warehouse is labeled as 0. Therefore, U = M\{0} is the set of customer nodes. The distance between each node is Lij. Every customer has a demand for Pi units; for calculating energy consumption, we measure demand in units of weight.In our UAV-vehicle mode, each vehicle with a payload (weight) capacity ofQV is equipped with a UAV. Each pair of a vehicle and a UAV has a fixed cost CF, as well as a variable route cost, which is a function of range and gross weight (see next paragraph). The unit route cost CU of the UAV is expected to be much lower than that of the vehicle CV. However, the payload (capacity) QU of the UAV and the distance LU and the time TU that the UAV can fly in the air are limited. Note that the time limit may include the time the UAV has to wait for the vehicle before landing. The distance between the vehicle and the UAV may be different; in fact, another potential advantage of UAVs is that they may choose more efficient routes than vehicles. For example, vehicles need to follow the Manhattan metric, while UAVs can use Euclidean distances. Therefore, we distinguish between the distance from i to j by the vehicle, LijV, and that by the UAV, LijU. It is assumed that the UAV will serve only one customer before returning to the vehicle but can serve other customers after returning the vehicle to reload and replace batteries.If the vehicle starts from the nodei to the node j, {i, j} ∈ N, let Nij be 1; otherwise, let it be 0. Let Mij be the vehicle payload weight for edge {i, j} ∈ N. Obviously, Mij = 0, if the vehicle does not run from the node i to the node j. Finally, if the UAV leaves the vehicle at node i, serves the qj unit required by customer j, and returns to node k and lands on the roof of the vehicle, then Zijk = 1 and the flight time is Tijk; otherwise, Zijk = 0.In the delivery area of the warehouse, there separately exist two kinds of traffic restrictions on the vehicles. Under the vehicle-type restriction (labeled as “circumstance 1”), there are some edges where vehicles are forbidden to drive. LetN¯1 be the set of these edges, then, for any {i, j} ∈ N¯1, Nij = 0. Under half-side traffic (labeled as “circumstance 2”), the average speed, variable cost, and CO2 emission of the vehicle on these edges are raised by RS, RC, and RE times, respectively. Let N¯2 be the set of these edges, then RijS=RijC=RijE=1 for any i,j∉N¯2 and RijS<1, RijC>1, and RijE>1 for any i,j∈N¯2.Franzese and Davidson [31] pointed out that the increase in total vehicle weight has a certain impact on its fuel efficiency, and fuel is an important component of variable transportation costs, so they added relevant factor of fuel to GVRP. Bateman et al. [32] pointed out that the carbon footprint of transportation emissions can be estimated by transportation weight and distance traveled.On the route from node i to j, the total vehicle weight isWij, including the empty weight of the vehicle WV, plus the payload weigh Mij, and the UAV weight WU if the UAV is on the vehicle. The weight can be expressed as a nonlinear function:(1a)Wij=WVNij+Mij+WUNij1−∑h∈UZihj.The first term on the right hand of (1) is the vehicle empty weight when running from i to j, the second is the weight of packages on the vehicle from i to j, and the third is the weight of the UAV if it is on the vehicle from i to j. Another linear formula is(1b)Wij≥WVNij+Mij+WUNij−∑h∈UZihj.Although (1b) may generate negative values for edges that do not exist in the solution, these values will be 0 because of (21). To calculate the vehicle route cost, we use the unit cost, CV, time distance, LijV, and time total vehicle weight, Wij (note that CV is calculated in dollars per pound-mile and can determined by regression analysis).On the route of the UAV flying fromi to j to k, the total weights of the UAV are different. On the route of the UAV flying from i to j, the total weight of the UAV is(2a)WijU=WU+qjZijk.On the route of the UAV flying fromj to k, the total weight of the UAV is(2b)WjkU=WUZijk.To calculate the vehicle route cost, we use the unit cost,CU, time distance, LijU and LjkU, and total UAV weight, WijU and WjkU.There are two ways to analyze the impact of UAV delivery. The first is the benefit on total CO2 emission. The total CO2 emission of a vehicle can be calculated by the ways of Goodchild and Toy [33] as follows:(3a)E1V=WAER∑i∈M∑j∈MLijV×Wij,(3b)E2V=WAER∑i∈M∑j∈MRijE×LijV×Wij,where WAER is the weighted average emission rate of vehicles. (3a) is the formulation of a vehicle’s total CO2 emission under circumstance 1, and (3b) is that under circumstance 2. Obviously, in (3a), Wij=0 for any {i, j} ∈ N¯1.The total CO2 emission of UAV is calculated as(4)EU=PGFER×AER∑i∈M∑j∈M∑k∈MLijUWijU+LjkUWjkU,where PGFER is the CO2 emission rate per watt-hour (Wh) of the power generation facilities in the generation of electricity for using the UAV and AER is the average energy requirement of UAV in Wh per pound-mile.Therefore, the environmental goals of minimizing total CO2 emission under circumstances 1 and 2 are, respectively,(5a)minimizeE1V+EU,(5b)minimizeE2V+EU.The alternative way is examining its impact on the traditional objective of minimizing total cost, TC, under circumstances 1 and 2, including the fixed cost of the pair of a UAV and a vehicle and the variable route cost of the vehicle and the UAV:(6a)TC1=CF∑j∈UN0,j+CV∑i∈M∑j∈MLijV×Wij+CU∑i∈M∑j∈U∑k∈MLijUWijU+LjkUWjkU,(6b)TC2=CF∑j∈UN0,j+CV∑i∈M∑j∈MRijC×LijV×Wij+CU∑i∈M∑j∈U∑k∈MLijUWijU+LjkUWjkU.Therefore, the traditional goals of minimizing total cost under circumstances 1 and 2 are, respectively, as follows:(7a)minimizeTC1,(7b)minimizeTC2LetVijV and VijU, respectively, be the average speed of the vehicle and the UAV from node i to node j (mile/hour). The constraint set of the two goals can be expressed as follows. It should be noted that the same set of constraints are applied to the above two objectives and is solved twice:(8)∑j∈MNij+∑h∈M∑k∈MZhik=1,∀j∈U,(9)∑i∈MNij+∑h∈M∑k∈MZhjk=1,∀j∈U,(10)∑j∈UN0j=∑i∈UNi0,(11)∑j∈Mj≠iMji−∑j∈Mj≠iMij+∑h∈M∑k∈MqiZhij−∑j∈M∑k∈MqjZijk=qi,∀i∈U,(12)Nik≥∑j∈MZijk,∀i,k∈M,(13)∑j∈Mqj+WU≤QV,(14)qj∑i∈M∑k∈MZijk≤QU,∀j∈U,(15)LijU+LjkUZhik≤LU,∀i,j,k∈M,(16a)Tijk≥LijUVijU+LjkUVjkUZhjk,∀i,k∈M,∀j∈U,(16b)Tijk≥LikVVikVNik,∀i,k∈M,∀j∈U,(16c)Tijk≥LikVRikSVikVNik,∀i,k∈M,∀j∈U,(16d)Tijk≤TU+G1−Zhjk,∀i,k∈M,∀j∈U,(17)Wij≥WVNij+Mij+WUNij−∑h∈UZihj,∀i,j∈N,(18)Mij≥qjNij,∀i,j∈N,(19a)Nij=0,∀i,j∈N¯1,(19b)RijS=RijC=RijE=1,∀i,j∉N¯2,(19c)RijS<1,RijC>1,RijE>1,∀i,j∈N¯2,(20a)∑i∈M∑k∈MLikVVikVNik+∑j∈UyijkTijk−LikVVikVNik≤Tw,(20b)∑i∈M∑k∈MLikVRikSVikVNik+∑j∈UyijkTijk−LikVRikSVikVNik≤Tw,(21)Tijk≥0,∀i,k∈M,∀j∈U,(22)Wij≥0,∀i,j∈N,(23)Nij=0,1,N00=0,∀i,j∈N,(24)Zijk∈0,1,∀i,j,k∈M.Constraints (8) and (9) ensure that each customer is served by a vehicle or a UAV. Constraint (10) ensures that the numbers of vehicles leaving and entering the warehouse are the same. Constraint (11) ensures that each customer’s demand is met. Constraint (12) synchronizes vehicle and UAV so that the vehicle picks up the UAV at node k. We will relax the restrictions (12) in Section 4 to allow vehicles to make multiple stops before receiving the UAV. Constraints (13) and (14) are constraints on vehicle and UAV payloads, respectively. The distance range limit and the time range limit of UAV are, respectively, guaranteed by constraints (15) and (16d), where is a very large value. The flight time of UAV, Tijk, is the larger one in the travel time of UAV (16a) and vehicle under circumstance 1 (16b) or vehicle under circumstance 2 (16c). The constraint (17) determines the vehicle weight described in the preceding paragraph. Note that a simple lower bound (constraint (18)) can be included on the load of the vehicle at any time. Constraint (19a) ensures vehicles are forbidden to run on the edges of traffic restriction under circumstance 1, and constraints (19b) and (19c) ensure that the average speed, variable cost, and CO2 emission of the vehicle on the edges of traffic restriction under circumstance 2 are changed. (20a) and (20b) are the constraints of working hours under circumstances 1 and 2, respectively, where Tw is the upper limit of working hours, normally 8 hours. ∑j∈UyijkTijk−LikV/VikVNik is the time length of the vehicle waiting for its UAV, where yijk=1 if Tijk>LikV/VikVNik; otherwise, yijk=0. Finally, we can avoid the use of UAVs for specific routes (perhaps near airports) or for specific customers (who cannot receive such deliveries) by setting Zijk = Zhij = 0 for route {i, j} or Zijk = 0 for customer j.


## 3.1. Problem Description

In an urban area, a delivery company uses vehicles equipped with UAVs to deliver parcels from a warehouse to customers in its delivery area. Each customer has only one parcel waiting to be delivered. This UVA-vehicle mode is as follows. Every vehicle is equipped with a UAV. As both delivery range (distance and flight time) and capacity (weight and size) of the UAV are limited, it only delivers parts of the parcels, but the vehicle can deliver all the parcels. Each UAV may depart from its vehicle at a location of a customer or the warehouse and carry a parcel to one and only one customer. Then, the UAV returns to its vehicle to reload a parcel and recharge or swap batteries, which is instantaneous. While the UAV delivers the parcel, the vehicle carries out its delivery to one or several customers. Therefore, the UAV returns to its vehicle at a different customer location. As a matter of course, if the vehicle arrives at the customer location, where the UAV is retrieved, it waits for the UAV, and vice versa. For the sake of the safety of the UAV, if the UAV arrives before the vehicle, it hovers in the air waiting for the vehicle.In the delivery area of the warehouse, there may be traffic restrictions on the vehicles. Generally, there are two kinds of traffic restrictions, under which this UAV-vehicle mode is practicable but greatly affected. The first is vehicle-type restriction, under which delivery vehicles are forbidden to run through some paths and have to choose other way or use a UAV. The second is half-side traffic, under which only one lane of a road is available, and vehicles traveling in opposite directions can only take turns on this road. Under this circumstance, the delivery vehicle can run on this road but may have to wait for its turn with the operating engine and raise the variable cost and CO2 emission.According to the regulation on working hours, the service time of every pair of UAV and vehicle is 8 hours. Therefore, it is probable that more than one pair of UAV and vehicle are needed.There are two goals of routing the UAV-vehicle. The first is minimizing the total cost, and the alternative is minimizing the energy consumption (or CO2 emission). The total cost consists of fixed cost whenever a pair of UAV-vehicle is used, as well as variable cost, which is the function of unit route cost, travel distance, and gross weight (empty weight of the vehicle or UAV plus the payload). The total CO2 emission of a vehicle is the function of weighted average emission rate of vehicles, travel distance, and gross weight, and that of a UAV is the function of CO2 emission rate of generating per watt-hour, average energy requirement of UAV, travel distance, and gross weight.It should be noted that the two kinds of traffic restrictions may simultaneously exist, but as we analyze the impact of UAV delivery under every kind of traffic restriction, they are in our models separately.


## 3.2. The Model

According to Laporte [30], there are three ways to formulate the vehicle routing problem (VRP), including the simple set division formula of VRP first proposed by Balinski and Quandt in 1964, the commodity flow formula of Shlifer and Graves in 1979, and the two-index vehicle flow formula of Laporte and Norbert in 1983. In our study, the GVRP problem is mainly involved, and the weight of vehicles, commodities, and drones in the whole route needs to be considered. Therefore, we choose to continue to expand the commodity flow formula, which is clearly expressed and promotes the development of heuristic algorithm.The basic GVRP problem is as follows: both the UAV and the vehicle can travel back and forth on any available path. Therefore, an undirected graphR = (M, N) is given, where N is the set of edges with nonnegative routing costs, {i, j} ∈ N represents the edges from node i to node j. M = {0, 1, …, n} is the set of all nodes with a total of n customers, and the warehouse is labeled as 0. Therefore, U = M\{0} is the set of customer nodes. The distance between each node is Lij. Every customer has a demand for Pi units; for calculating energy consumption, we measure demand in units of weight.In our UAV-vehicle mode, each vehicle with a payload (weight) capacity ofQV is equipped with a UAV. Each pair of a vehicle and a UAV has a fixed cost CF, as well as a variable route cost, which is a function of range and gross weight (see next paragraph). The unit route cost CU of the UAV is expected to be much lower than that of the vehicle CV. However, the payload (capacity) QU of the UAV and the distance LU and the time TU that the UAV can fly in the air are limited. Note that the time limit may include the time the UAV has to wait for the vehicle before landing. The distance between the vehicle and the UAV may be different; in fact, another potential advantage of UAVs is that they may choose more efficient routes than vehicles. For example, vehicles need to follow the Manhattan metric, while UAVs can use Euclidean distances. Therefore, we distinguish between the distance from i to j by the vehicle, LijV, and that by the UAV, LijU. It is assumed that the UAV will serve only one customer before returning to the vehicle but can serve other customers after returning the vehicle to reload and replace batteries.If the vehicle starts from the nodei to the node j, {i, j} ∈ N, let Nij be 1; otherwise, let it be 0. Let Mij be the vehicle payload weight for edge {i, j} ∈ N. Obviously, Mij = 0, if the vehicle does not run from the node i to the node j. Finally, if the UAV leaves the vehicle at node i, serves the qj unit required by customer j, and returns to node k and lands on the roof of the vehicle, then Zijk = 1 and the flight time is Tijk; otherwise, Zijk = 0.In the delivery area of the warehouse, there separately exist two kinds of traffic restrictions on the vehicles. Under the vehicle-type restriction (labeled as “circumstance 1”), there are some edges where vehicles are forbidden to drive. LetN¯1 be the set of these edges, then, for any {i, j} ∈ N¯1, Nij = 0. Under half-side traffic (labeled as “circumstance 2”), the average speed, variable cost, and CO2 emission of the vehicle on these edges are raised by RS, RC, and RE times, respectively. Let N¯2 be the set of these edges, then RijS=RijC=RijE=1 for any i,j∉N¯2 and RijS<1, RijC>1, and RijE>1 for any i,j∈N¯2.Franzese and Davidson [31] pointed out that the increase in total vehicle weight has a certain impact on its fuel efficiency, and fuel is an important component of variable transportation costs, so they added relevant factor of fuel to GVRP. Bateman et al. [32] pointed out that the carbon footprint of transportation emissions can be estimated by transportation weight and distance traveled.On the route from node i to j, the total vehicle weight isWij, including the empty weight of the vehicle WV, plus the payload weigh Mij, and the UAV weight WU if the UAV is on the vehicle. The weight can be expressed as a nonlinear function:(1a)Wij=WVNij+Mij+WUNij1−∑h∈UZihj.The first term on the right hand of (1) is the vehicle empty weight when running from i to j, the second is the weight of packages on the vehicle from i to j, and the third is the weight of the UAV if it is on the vehicle from i to j. Another linear formula is(1b)Wij≥WVNij+Mij+WUNij−∑h∈UZihj.Although (1b) may generate negative values for edges that do not exist in the solution, these values will be 0 because of (21). To calculate the vehicle route cost, we use the unit cost, CV, time distance, LijV, and time total vehicle weight, Wij (note that CV is calculated in dollars per pound-mile and can determined by regression analysis).On the route of the UAV flying fromi to j to k, the total weights of the UAV are different. On the route of the UAV flying from i to j, the total weight of the UAV is(2a)WijU=WU+qjZijk.On the route of the UAV flying fromj to k, the total weight of the UAV is(2b)WjkU=WUZijk.To calculate the vehicle route cost, we use the unit cost,CU, time distance, LijU and LjkU, and total UAV weight, WijU and WjkU.There are two ways to analyze the impact of UAV delivery. The first is the benefit on total CO2 emission. The total CO2 emission of a vehicle can be calculated by the ways of Goodchild and Toy [33] as follows:(3a)E1V=WAER∑i∈M∑j∈MLijV×Wij,(3b)E2V=WAER∑i∈M∑j∈MRijE×LijV×Wij,where WAER is the weighted average emission rate of vehicles. (3a) is the formulation of a vehicle’s total CO2 emission under circumstance 1, and (3b) is that under circumstance 2. Obviously, in (3a), Wij=0 for any {i, j} ∈ N¯1.The total CO2 emission of UAV is calculated as(4)EU=PGFER×AER∑i∈M∑j∈M∑k∈MLijUWijU+LjkUWjkU,where PGFER is the CO2 emission rate per watt-hour (Wh) of the power generation facilities in the generation of electricity for using the UAV and AER is the average energy requirement of UAV in Wh per pound-mile.Therefore, the environmental goals of minimizing total CO2 emission under circumstances 1 and 2 are, respectively,(5a)minimizeE1V+EU,(5b)minimizeE2V+EU.The alternative way is examining its impact on the traditional objective of minimizing total cost, TC, under circumstances 1 and 2, including the fixed cost of the pair of a UAV and a vehicle and the variable route cost of the vehicle and the UAV:(6a)TC1=CF∑j∈UN0,j+CV∑i∈M∑j∈MLijV×Wij+CU∑i∈M∑j∈U∑k∈MLijUWijU+LjkUWjkU,(6b)TC2=CF∑j∈UN0,j+CV∑i∈M∑j∈MRijC×LijV×Wij+CU∑i∈M∑j∈U∑k∈MLijUWijU+LjkUWjkU.Therefore, the traditional goals of minimizing total cost under circumstances 1 and 2 are, respectively, as follows:(7a)minimizeTC1,(7b)minimizeTC2LetVijV and VijU, respectively, be the average speed of the vehicle and the UAV from node i to node j (mile/hour). The constraint set of the two goals can be expressed as follows. It should be noted that the same set of constraints are applied to the above two objectives and is solved twice:(8)∑j∈MNij+∑h∈M∑k∈MZhik=1,∀j∈U,(9)∑i∈MNij+∑h∈M∑k∈MZhjk=1,∀j∈U,(10)∑j∈UN0j=∑i∈UNi0,(11)∑j∈Mj≠iMji−∑j∈Mj≠iMij+∑h∈M∑k∈MqiZhij−∑j∈M∑k∈MqjZijk=qi,∀i∈U,(12)Nik≥∑j∈MZijk,∀i,k∈M,(13)∑j∈Mqj+WU≤QV,(14)qj∑i∈M∑k∈MZijk≤QU,∀j∈U,(15)LijU+LjkUZhik≤LU,∀i,j,k∈M,(16a)Tijk≥LijUVijU+LjkUVjkUZhjk,∀i,k∈M,∀j∈U,(16b)Tijk≥LikVVikVNik,∀i,k∈M,∀j∈U,(16c)Tijk≥LikVRikSVikVNik,∀i,k∈M,∀j∈U,(16d)Tijk≤TU+G1−Zhjk,∀i,k∈M,∀j∈U,(17)Wij≥WVNij+Mij+WUNij−∑h∈UZihj,∀i,j∈N,(18)Mij≥qjNij,∀i,j∈N,(19a)Nij=0,∀i,j∈N¯1,(19b)RijS=RijC=RijE=1,∀i,j∉N¯2,(19c)RijS<1,RijC>1,RijE>1,∀i,j∈N¯2,(20a)∑i∈M∑k∈MLikVVikVNik+∑j∈UyijkTijk−LikVVikVNik≤Tw,(20b)∑i∈M∑k∈MLikVRikSVikVNik+∑j∈UyijkTijk−LikVRikSVikVNik≤Tw,(21)Tijk≥0,∀i,k∈M,∀j∈U,(22)Wij≥0,∀i,j∈N,(23)Nij=0,1,N00=0,∀i,j∈N,(24)Zijk∈0,1,∀i,j,k∈M.Constraints (8) and (9) ensure that each customer is served by a vehicle or a UAV. Constraint (10) ensures that the numbers of vehicles leaving and entering the warehouse are the same. Constraint (11) ensures that each customer’s demand is met. Constraint (12) synchronizes vehicle and UAV so that the vehicle picks up the UAV at node k. We will relax the restrictions (12) in Section 4 to allow vehicles to make multiple stops before receiving the UAV. Constraints (13) and (14) are constraints on vehicle and UAV payloads, respectively. The distance range limit and the time range limit of UAV are, respectively, guaranteed by constraints (15) and (16d), where is a very large value. The flight time of UAV, Tijk, is the larger one in the travel time of UAV (16a) and vehicle under circumstance 1 (16b) or vehicle under circumstance 2 (16c). The constraint (17) determines the vehicle weight described in the preceding paragraph. Note that a simple lower bound (constraint (18)) can be included on the load of the vehicle at any time. Constraint (19a) ensures vehicles are forbidden to run on the edges of traffic restriction under circumstance 1, and constraints (19b) and (19c) ensure that the average speed, variable cost, and CO2 emission of the vehicle on the edges of traffic restriction under circumstance 2 are changed. (20a) and (20b) are the constraints of working hours under circumstances 1 and 2, respectively, where Tw is the upper limit of working hours, normally 8 hours. ∑j∈UyijkTijk−LikV/VikVNik is the time length of the vehicle waiting for its UAV, where yijk=1 if Tijk>LikV/VikVNik; otherwise, yijk=0. Finally, we can avoid the use of UAVs for specific routes (perhaps near airports) or for specific customers (who cannot receive such deliveries) by setting Zijk = Zhij = 0 for route {i, j} or Zijk = 0 for customer j.


## 4. Solution Methodology

The application of our genetic algorithm in our routing problem with UAV is described below.


### 4.1. Design

In the design of our genetic algorithm, we adapt the following service strategy. The work should be completed within one day; that is, all parcels should be delivered within legal working hours, normally 8 hours. Meanwhile, fuel consumption should be minimized, or CO2 emission should be minimized. Firstly, we can run our algorithm, where only a pair of a UAV and a vehicle is used, to calculate the total working hours. Then, the number of needed vehicles is determined by dividing the total working hours by the legal working hours. After that, the region is divided according to the number of vehicles by a line (or lines) crossing the warehouse, where the number of customers served by each vehicle and the service time of each vehicle are as close to each other as possible.In our algorithm, we first sequence the customers (please see the methods in4.1.2). The pair of a vehicle and a UAV serves customers according to the sequence. Then, we group customers randomly as follows. Generate a random positive integer x∈1,n, where n is the number of total customers served by a pair of a UAV and a vehicle. And the first x customers are put into the first group. Repeat it till all customers are grouped. As a matter of course, if there are more than one edge where vehicles are forbidden, and they are divided into different groups. If there is only one customer in the group, the vehicle serves the customer. Otherwise, the UAV serves a customer randomly selected from the group, except the last one, where the UAV is retrieved, and the vehicle serves the others. If there is an edge connecting two customers in the group, where vehicles are forbidden, we randomly select one from these two customers for the UAV to serve. Based on this strategy, we developed a genetic algorithm to solve the proposed minimization problem (fuel consumption and CO2 emissions or total cost).


#### 4.1.1. Encoding

We use natural number coding method to encode customers, warehouse, and group information. That is, we use positive integers from 1 ton to encode customers and encode the warehouse as 0. And we randomly generated positive integer(s) to encode the number of customers in a service group.


#### 4.1.2. Initial Population Generation

In our GA, the initial population is artificially optimized to reduce iteration times and improve optimization efficiency. Specifically, we equally divide the service area into subareas of an even number. Then, we arrange the service sequence of the initial individuals as follows. The pair of a vehicle and a UAV serves customers in one of the nearest subareas from the nearest customer to the farthest customer and goes to an adjacent subarea to serve customers from the farthest to the nearest. Repeat it until all customers are served, and both the UAV and the vehicle come back into the warehouse. After that, we group the customers and obtain the initial population generation.The following is an example to illustrate the service strategy, grouping, encoding, and initial population generation of our algorithm. There are 20 customers who should be satisfied. We encode them as {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20} and the warehouse as 0. As only 1 pair of a UAV and a vehicle is needed, they serve the whole delivery area. Then, we divide the whole delivery area into eight subareas and arrange the service sequence of the initial individuals as follows. The pair of a UAV and a vehicle serves customers in the first subarea from left to right and goes up to the second subarea to serve customers from right to left. Then, they go left to the third subarea to serve customers from right to left and go down to the fourth subarea to serve customers from left to right. After that, they go down to the fifth subarea to serve customers from right to left and go down to the sixth subarea to serve customers from left to right. Finally, they go right to the seventh subarea to serve customers from left to right and go up to the eighth to serve customers from right to left and go back to the warehouse. Now, we obtain the service sequence codes as {18,1,3,7,9,19,6,12,16,8,14,2,13,15,20,11,17,4,5,10} (Figure1).Figure 1
An example of partitioning service area and sequencing customers.Then, we randomly generate the group codes as {2,5,3,4,5,1}. So, the customers are divided into 6 groups with encoding result {{18,1}, {3,7,9,19,6}, {12,16,8}, {14,2,13,15}, {20,11,17,4,5}, {10}}. As there are two customers, 18 and 1, in the first group, the UAV departs from the vehicle at the location of the warehouse to deliver parcel to customer 18 and fly to customer 1. Meanwhile, the vehicle goes to serve customer 1 and retrieve the UAV there. There are 5 customers, 3, 7, 9, 19, and 6, in the second group. As customer 6 is the last one, we randomly choose 7 from 3, 7, 9, and 19 to be the customer served by the UAV. So, the UAV departs from customer 13 to serve 7 and fly to 6, while the UAV sequentially serves customers 3, 9, 19, and 6 and retrieves the UAV at 6. There are 3 customers in the third group. We randomly choose the first customer, 12, to be served by the UAV. As a result, the UAVs fly from 6 to 12 and to 8, while the vehicle runs from 6 to 16 and to 8 and picks up the UAV here. For the fourth group, the UAV serves customer 14, and the vehicle sequentially serves 2, 13, and 15 and retrieves the UAV here. For the fifth group, the UAV serves 4, and the vehicle serves 20, 11, 17, and 5 and retrieves the UAV here. It should be noticed that although there is only customer 10 in the last group, the UAV should depart from customer 5 to serve customer 10 and fly back to the warehouse, and the vehicle should directly run from customer 10 back to the warehouse. The initial population generation is shown in Figure2.Figure 2
An example of an initial population generation.


#### 4.1.3. The Algorithms

(1) Fitness Function. In Section 3, we explained the design of the model in detail, in which the objective function mentioned mainly includes the load of the vehicle and the UAV, the distance between customers and the corresponding energy consumption function and cost function. Accordingly, the fitness function can be defined as fitnessX=1/TX, where X is the combination of service sequence and grouping and TX is the target function.(2) Selection, Crossover, and Variation. In the process of optimal individual selection, the greater the fitness of an individual, the greater the probability of being selected as the parent of the next generation, and the method is similar to roulette. As long as the fitness of an individual is large enough, the same individual in the genetic algorithm can exist in different generations. In the specific selection, select the elite individual retention strategy and copy the individuals with the highest fitness to the alternative parent group of cross-matching, then cross match all the individuals in the parent group, select the individuals with the highest fitness, and repeat this process, so as to gradually eliminate the individuals with poor fitness, leaving the elite individuals. At the same time, the method of double point crossing and real-value variation is adopted.(3) Termination Conditions. After several iterations, when the output iteration results meet the following conditions, the convergence value is reached and the algorithm is terminated: (1)
There was no significant change in the fitness value after continuous iteration; in other words, the change is less than 1‰ of the fitness value(2)
The population no longer evolves and the number of iterations reaches the set value


### 4.2. Algorithm Steps

In this paper, the algorithm steps are as follows:(1)
Determine the coding mechanism, generate the initial population, and set the probability of cross variation and the maximum number of iterations.(2)
Evaluate the fitness of the current population and find the chromosome with the minimum objective function value and corresponding function value.(3)
Determine whether the termination conditions are met. Output the optimal chromosome and corresponding solutions if satisfied; otherwise, select the current population again.(4)
Repeated crossover and mutation in the newly obtained generation.(5)
Evaluate the fitness of the results and move to (3).


## 4.1. Design

In the design of our genetic algorithm, we adapt the following service strategy. The work should be completed within one day; that is, all parcels should be delivered within legal working hours, normally 8 hours. Meanwhile, fuel consumption should be minimized, or CO2 emission should be minimized. Firstly, we can run our algorithm, where only a pair of a UAV and a vehicle is used, to calculate the total working hours. Then, the number of needed vehicles is determined by dividing the total working hours by the legal working hours. After that, the region is divided according to the number of vehicles by a line (or lines) crossing the warehouse, where the number of customers served by each vehicle and the service time of each vehicle are as close to each other as possible.In our algorithm, we first sequence the customers (please see the methods in4.1.2). The pair of a vehicle and a UAV serves customers according to the sequence. Then, we group customers randomly as follows. Generate a random positive integer x∈1,n, where n is the number of total customers served by a pair of a UAV and a vehicle. And the first x customers are put into the first group. Repeat it till all customers are grouped. As a matter of course, if there are more than one edge where vehicles are forbidden, and they are divided into different groups. If there is only one customer in the group, the vehicle serves the customer. Otherwise, the UAV serves a customer randomly selected from the group, except the last one, where the UAV is retrieved, and the vehicle serves the others. If there is an edge connecting two customers in the group, where vehicles are forbidden, we randomly select one from these two customers for the UAV to serve. Based on this strategy, we developed a genetic algorithm to solve the proposed minimization problem (fuel consumption and CO2 emissions or total cost).


### 4.1.1. Encoding

We use natural number coding method to encode customers, warehouse, and group information. That is, we use positive integers from 1 ton to encode customers and encode the warehouse as 0. And we randomly generated positive integer(s) to encode the number of customers in a service group.


### 4.1.2. Initial Population Generation

In our GA, the initial population is artificially optimized to reduce iteration times and improve optimization efficiency. Specifically, we equally divide the service area into subareas of an even number. Then, we arrange the service sequence of the initial individuals as follows. The pair of a vehicle and a UAV serves customers in one of the nearest subareas from the nearest customer to the farthest customer and goes to an adjacent subarea to serve customers from the farthest to the nearest. Repeat it until all customers are served, and both the UAV and the vehicle come back into the warehouse. After that, we group the customers and obtain the initial population generation.The following is an example to illustrate the service strategy, grouping, encoding, and initial population generation of our algorithm. There are 20 customers who should be satisfied. We encode them as {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20} and the warehouse as 0. As only 1 pair of a UAV and a vehicle is needed, they serve the whole delivery area. Then, we divide the whole delivery area into eight subareas and arrange the service sequence of the initial individuals as follows. The pair of a UAV and a vehicle serves customers in the first subarea from left to right and goes up to the second subarea to serve customers from right to left. Then, they go left to the third subarea to serve customers from right to left and go down to the fourth subarea to serve customers from left to right. After that, they go down to the fifth subarea to serve customers from right to left and go down to the sixth subarea to serve customers from left to right. Finally, they go right to the seventh subarea to serve customers from left to right and go up to the eighth to serve customers from right to left and go back to the warehouse. Now, we obtain the service sequence codes as {18,1,3,7,9,19,6,12,16,8,14,2,13,15,20,11,17,4,5,10} (Figure1).Figure 1
An example of partitioning service area and sequencing customers.Then, we randomly generate the group codes as {2,5,3,4,5,1}. So, the customers are divided into 6 groups with encoding result {{18,1}, {3,7,9,19,6}, {12,16,8}, {14,2,13,15}, {20,11,17,4,5}, {10}}. As there are two customers, 18 and 1, in the first group, the UAV departs from the vehicle at the location of the warehouse to deliver parcel to customer 18 and fly to customer 1. Meanwhile, the vehicle goes to serve customer 1 and retrieve the UAV there. There are 5 customers, 3, 7, 9, 19, and 6, in the second group. As customer 6 is the last one, we randomly choose 7 from 3, 7, 9, and 19 to be the customer served by the UAV. So, the UAV departs from customer 13 to serve 7 and fly to 6, while the UAV sequentially serves customers 3, 9, 19, and 6 and retrieves the UAV at 6. There are 3 customers in the third group. We randomly choose the first customer, 12, to be served by the UAV. As a result, the UAVs fly from 6 to 12 and to 8, while the vehicle runs from 6 to 16 and to 8 and picks up the UAV here. For the fourth group, the UAV serves customer 14, and the vehicle sequentially serves 2, 13, and 15 and retrieves the UAV here. For the fifth group, the UAV serves 4, and the vehicle serves 20, 11, 17, and 5 and retrieves the UAV here. It should be noticed that although there is only customer 10 in the last group, the UAV should depart from customer 5 to serve customer 10 and fly back to the warehouse, and the vehicle should directly run from customer 10 back to the warehouse. The initial population generation is shown in Figure2.Figure 2
An example of an initial population generation.


### 4.1.3. The Algorithms

(1) Fitness Function. In Section 3, we explained the design of the model in detail, in which the objective function mentioned mainly includes the load of the vehicle and the UAV, the distance between customers and the corresponding energy consumption function and cost function. Accordingly, the fitness function can be defined as fitnessX=1/TX, where X is the combination of service sequence and grouping and TX is the target function.(2) Selection, Crossover, and Variation. In the process of optimal individual selection, the greater the fitness of an individual, the greater the probability of being selected as the parent of the next generation, and the method is similar to roulette. As long as the fitness of an individual is large enough, the same individual in the genetic algorithm can exist in different generations. In the specific selection, select the elite individual retention strategy and copy the individuals with the highest fitness to the alternative parent group of cross-matching, then cross match all the individuals in the parent group, select the individuals with the highest fitness, and repeat this process, so as to gradually eliminate the individuals with poor fitness, leaving the elite individuals. At the same time, the method of double point crossing and real-value variation is adopted.(3) Termination Conditions. After several iterations, when the output iteration results meet the following conditions, the convergence value is reached and the algorithm is terminated: (1)
There was no significant change in the fitness value after continuous iteration; in other words, the change is less than 1‰ of the fitness value(2)
The population no longer evolves and the number of iterations reaches the set value


## 4.1.1. Encoding

We use natural number coding method to encode customers, warehouse, and group information. That is, we use positive integers from 1 ton to encode customers and encode the warehouse as 0. And we randomly generated positive integer(s) to encode the number of customers in a service group.


## 4.1.2. Initial Population Generation

In our GA, the initial population is artificially optimized to reduce iteration times and improve optimization efficiency. Specifically, we equally divide the service area into subareas of an even number. Then, we arrange the service sequence of the initial individuals as follows. The pair of a vehicle and a UAV serves customers in one of the nearest subareas from the nearest customer to the farthest customer and goes to an adjacent subarea to serve customers from the farthest to the nearest. Repeat it until all customers are served, and both the UAV and the vehicle come back into the warehouse. After that, we group the customers and obtain the initial population generation.The following is an example to illustrate the service strategy, grouping, encoding, and initial population generation of our algorithm. There are 20 customers who should be satisfied. We encode them as {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20} and the warehouse as 0. As only 1 pair of a UAV and a vehicle is needed, they serve the whole delivery area. Then, we divide the whole delivery area into eight subareas and arrange the service sequence of the initial individuals as follows. The pair of a UAV and a vehicle serves customers in the first subarea from left to right and goes up to the second subarea to serve customers from right to left. Then, they go left to the third subarea to serve customers from right to left and go down to the fourth subarea to serve customers from left to right. After that, they go down to the fifth subarea to serve customers from right to left and go down to the sixth subarea to serve customers from left to right. Finally, they go right to the seventh subarea to serve customers from left to right and go up to the eighth to serve customers from right to left and go back to the warehouse. Now, we obtain the service sequence codes as {18,1,3,7,9,19,6,12,16,8,14,2,13,15,20,11,17,4,5,10} (Figure1).Figure 1
An example of partitioning service area and sequencing customers.Then, we randomly generate the group codes as {2,5,3,4,5,1}. So, the customers are divided into 6 groups with encoding result {{18,1}, {3,7,9,19,6}, {12,16,8}, {14,2,13,15}, {20,11,17,4,5}, {10}}. As there are two customers, 18 and 1, in the first group, the UAV departs from the vehicle at the location of the warehouse to deliver parcel to customer 18 and fly to customer 1. Meanwhile, the vehicle goes to serve customer 1 and retrieve the UAV there. There are 5 customers, 3, 7, 9, 19, and 6, in the second group. As customer 6 is the last one, we randomly choose 7 from 3, 7, 9, and 19 to be the customer served by the UAV. So, the UAV departs from customer 13 to serve 7 and fly to 6, while the UAV sequentially serves customers 3, 9, 19, and 6 and retrieves the UAV at 6. There are 3 customers in the third group. We randomly choose the first customer, 12, to be served by the UAV. As a result, the UAVs fly from 6 to 12 and to 8, while the vehicle runs from 6 to 16 and to 8 and picks up the UAV here. For the fourth group, the UAV serves customer 14, and the vehicle sequentially serves 2, 13, and 15 and retrieves the UAV here. For the fifth group, the UAV serves 4, and the vehicle serves 20, 11, 17, and 5 and retrieves the UAV here. It should be noticed that although there is only customer 10 in the last group, the UAV should depart from customer 5 to serve customer 10 and fly back to the warehouse, and the vehicle should directly run from customer 10 back to the warehouse. The initial population generation is shown in Figure2.Figure 2
An example of an initial population generation.


## 4.1.3. The Algorithms

(1) Fitness Function. In Section 3, we explained the design of the model in detail, in which the objective function mentioned mainly includes the load of the vehicle and the UAV, the distance between customers and the corresponding energy consumption function and cost function. Accordingly, the fitness function can be defined as fitnessX=1/TX, where X is the combination of service sequence and grouping and TX is the target function.(2) Selection, Crossover, and Variation. In the process of optimal individual selection, the greater the fitness of an individual, the greater the probability of being selected as the parent of the next generation, and the method is similar to roulette. As long as the fitness of an individual is large enough, the same individual in the genetic algorithm can exist in different generations. In the specific selection, select the elite individual retention strategy and copy the individuals with the highest fitness to the alternative parent group of cross-matching, then cross match all the individuals in the parent group, select the individuals with the highest fitness, and repeat this process, so as to gradually eliminate the individuals with poor fitness, leaving the elite individuals. At the same time, the method of double point crossing and real-value variation is adopted.(3) Termination Conditions. After several iterations, when the output iteration results meet the following conditions, the convergence value is reached and the algorithm is terminated: (1)
There was no significant change in the fitness value after continuous iteration; in other words, the change is less than 1‰ of the fitness value(2)
The population no longer evolves and the number of iterations reaches the set value


## 4.2. Algorithm Steps

In this paper, the algorithm steps are as follows:(1)
Determine the coding mechanism, generate the initial population, and set the probability of cross variation and the maximum number of iterations.(2)
Evaluate the fitness of the current population and find the chromosome with the minimum objective function value and corresponding function value.(3)
Determine whether the termination conditions are met. Output the optimal chromosome and corresponding solutions if satisfied; otherwise, select the current population again.(4)
Repeated crossover and mutation in the newly obtained generation.(5)
Evaluate the fitness of the results and move to (3).


## 5. Experimentation

In this section, we use the software Matlab R2017a to carry out the experiment and test our model and GA. The computer was used in the experiment is ASUS Rock 5 Generation, whose configuration is Inter(R) Core(TM) i7-8550u processor (1.8 GHz base, 1.99 GHz Max Turbo), 128 GB SSD, 8 GB RAM.


### 5.1. Data Parameters

The range of customer demand,qj, is 1∼10 pounds, among which 80% is no more than 5 pounds, and the position of each customer (x, y) coordinate values x and y randomly evenly generated in the interval [−10, +10]. The distance of the vehicle from customer i to j is LijV, following the Manhattan metric; the distance of the UAV from customer i to j is LijU, following the Euclidean metric; in addition, the warehouse is located at the center of the region (0,0).Referring to Chiang et al. [2], we set the other parameters as follows.The empty weight of the UAV,WU, is 55 pounds; the payload of the UAV, MU, is 5 pounds; the maximum flight distance of the UAV, LU, is 10 miles; and the average flight speed of the UAV from customer i to j, VijU, is 25 miles per hour. For roads with the second kind of traffic restriction, RijC=RijE=1.2 and RijS=0.8333.The empty weight of the vehicle (truck),WV, is 6100 pounds, its payload, MV, is 6000 pounds, and the average speed of the vehicle from customer i to j, VijV, is 25 miles per hour.The fixed cost of a pair of a vehicle and a UAV,CF, is $500, the unit route cost of a vehicle, CV, is $0.00016/pound-mile, and the unit route cost of a UAV, CU, is $0.00036364/pound-mile.Referring to Goodchild and Toy [33], the weighted average emission of vehicles, WAER, is 1.2603 Kg/pound-mile, the CO2 emission from the power generation facilities in the generation of electricity, PGFER, is 3.773 × 10−4 Kg/Wh, and the average energy requirement of UAV, AER, is 3.3333 Wh/pound-mile.


### 5.2. The Examples

Firstly, we present an example of 20 customers under the two separate kinds of traffic restrictions, in order to analyze the efficiency of UAV delivery overcoming traffic restrictions. When there are 20 customers, only one vehicle or a pair of a vehicle and a UAV is needed. Figure3(a) shows the routing of a vehicle under no traffic restriction, and Figures 3(b)–3(d) separately show the routings of a vehicle and a pair of a vehicle and a UAV under two kinds of traffic restrictions. From Figure 3(d), we can find that the routings under two kinds of traffic restrictions are the same.Figure 3
Routings under different circumstances: (a) routing of only vehicle under no traffic restriction; (b) routing of only vehicle under the first kind of traffic restriction; (c) routing of only vehicle under the second kind of traffic restriction; (d) routing of vehicle with UAV under two kinds of traffic restrictions.
(a)(b)(c)(d)The results of costs and emissions under different circumstances are as follows. When a vehicle is under no traffic restriction, the CO2 emission is 161.3032 kg, and the variable cost is $134.09. When a vehicle is under the first kind of traffic restriction, the CO2 emission is 167.4327 Kg, and the variable cost is $139.19. When a vehicle is under the second kind of traffic restriction, the CO2 emission is 166.3727 Kg, and the variable cost is $138.30. When a pair of a vehicle is under both the first and the second kinds of traffic restriction, the CO2 emissions are the same, which is 151.6250 kg, and the variable cost is the same too, which is $113.97.From Figure3(d) and the result of costs and emissions, we can find that the routing as well as cost and emission under the first kind of traffic restriction are the same as those under the second. The main reason is no matter which kind of the traffic restriction is, we always use a UAV instead of a vehicle to deliver a parcel when facing a traffic restriction.Comparing Figures3(a)–3(d), as well as the corresponding result of costs and emissions, we can obtain the following conclusions.Firstly, traffic restrictions raise the cost and emission, and first kind of traffic restriction raises them further higher. The main reason is that when a road in the routing is blocked for the delivery vehicle, it has to choose other road, while if there is only one lane available, the vehicle could still run through the road if the cost and emission are lower than those running through other roads. Therefore, the cost and emission under the second kinds of traffic restriction are higher than those under no traffic restriction but lower than those under the first kind of traffic restriction.Secondly, UAV delivery can always reduce the cost and emission no matter if there is traffic restriction and which kind of traffic restriction is. Furthermore, UAV delivery saves more cost and cuts more emission under the first kind of traffic restriction than that under the second, as the cost and emission under the first kind of traffic restriction are the highest and those under two kinds of traffic restrictions are the same.In the following section, firstly, we present the result of 200 customers under the first kind of traffic restriction. The location of the 200 customers and the weight of their demand of were randomly generated and shown in Table1. 2 cars are needed for customer service within 8 hours per day. The service times of each vehicle are 6.65 hours and 6.94 hours, respectively. The routing of the two pairs of vehicles and drones is shown in Figure 4, and the results of costs and emissions are as follows.Table 1
Location of the 200 customers and the weight of their demand.
CustomerLocationWeight of demand (pound)CustomerLocationWeight of demand (pound)1(3.5985, 3.6248)2.6464101(6.6985, −7.8311)2.46952(1.6070, 2.0104)2.0035102(−1.9531, −3.3253)4.63763(6.7836, 3.9658)4.9405103(−7.8722, −1.6863)5.39564(2.7794, 2.7118)0.4074104(5.0983, −0.4058)5.52065(5.4746, −1.4383)5.5835105(−7.3811, −9.4588)2.78916(3.7286, −5.7173)2.6218106(−9.7011, −9.1042)0.60357(−3.9446, −0.7422)0.9957107(−8.0567, 4.1457)5.72418(−6.6992, 6.3058)0.8271108(−6.8909, −9.5609)2.01719(−8.7527, −4.0516)0.7294109(−5.7863, 4.4279)6.110510(−4.4899, −5.3138)1.3197110(7.3203, 5.0068)2.323511(4.4820, −3.0529)6.0178111(8.0247, −7.8543)4.21312(1.4923, −1.8090)4.5951112(−5.0841, 1.2691)5.537813(−8.3485, 2.9984)3.7066113(3.0984, −5.5088)4.010314(3.5833, −0.3746)4.2485114(−8.4013, −2.8959)3.906615(−8.8362, 3.5809)4.8984115(−1.5125, 3.4101)4.594616(−3.4012, 2.5994)2.6331116(4.7787, 8.0351)2.548117(−3.2365, 0.3061)5.1189117(0.5036, −1.8153)4.485918(−9.3533, 3.8708)2.3671118(9.2231, −5.6816)6.126519(5.7948, 3.0688)1.6931119(−0.7709, −0.2432)0.838620(−9.3260, 2.3723)6.0346120(0.1749, 0.4201)2.517621(5.5314, −8.4922)3.9256121(6.5870, −5.7042)3.450622(−5.3389, −1.8910)4.0786122(5.1065, −4.8937)5.618423(5.0150, 1.0702)0.5409123(−6.5286, −4.1542)4.541124(5.3291, 8.8109)0.726124(4.8201, 0.3814)1.821525(5.1424, 2.8790)3.9069125(−8.0454, 2.6018)3.855326(−8.5682, −4.6079)0.2204126(5.8765, 0.1324)0.41327(2.4282, −7.1526)5.7625127(−8.9043, 2.6839)2.514728(5.3858, 3.9510)5.8148128(3.2194, 4.7627)3.514829(−3.1541, −4.4266)6.0896129(−3.9258, 0.0123)4.780830(6.9408, 0.1415)1.0124130(0.0193, 9.3367)4.643631(6.9864, 6.0812)2.2663131(−1.7297, −2.4018)5.767332(−2.2162, −7.0139)2.0512132(7.9912, 3.1555)4.044733(1.7152, 6.3001)1.4065133(−1.3450, 1.1250)2.950434(7.2730, 1.6194)0.4455134(4.6406, 2.6336)6.070435(−9.0451, −1.9021)5.6629135(−1.8672, −3.1354)3.754736(3.9235, 3.1469)5.9208136(4.4430, 7.4150)0.849137(−4.1134, 9.0178)0.7815137(−7.3532, 3.6512)0.828338(4.1405, 0.6544)5.1735138(6.8433, 5.5395)2.508139(−9.0047, 2.8633)5.2433139(8.1880, 8.9098)2.848740(3.9065, −5.4808)3.8613140(−8.6609, −7.1078)4.79841(6.8066, 6.7106)1.008141(6.5072, 1.6621)4.981542(1.4364, −6.4336)4.315142(−1.7420, 8.7680)1.447243(−8.9608, −7.4226)5.8182143(−8.7449, −4.7519)3.411144(0.6380, 1.0638)5.9632144(2.5800, 9.5467)5.510645(−5.3386, −5.5256)1.804145(−6.2282, 2.4928)3.412546(−2.7935, 0.2916)3.8737146(9.7491, 2.2369)3.690847(−3.2275, −7.6894)6.1055147(5.0762, −8.1123)3.973848(6.0844, 5.0847)6.0703148(−2.1303, 7.2874)2.636149(−5.1356, −5.9556)3.0152149(−6.2300, −6.5167)2.69950(2.6515, 7.7721)0.0072150(6.2027, −0.4723)4.616651(−0.0056, −0.1799)3.3836151(2.0571, 3.5960)5.895652(2.7594, −4.9004)5.8689152(9.7581, 5.6005)1.163253(6.6679, 7.1263)5.246153(5.2829, 1.7994)2.592554(9.8432, 3.9686)1.3284154(−2.2265, 6.9336)2.374155(−5.3313, 6.4587)1.3662155(−8.9635, 6.1876)3.386456(−7.3582, 2.8845)5.5513156(7.3399, −1.8563)2.569457(0.1812, 2.1129)3.2123157(−3.4740, −8.4158)5.95158(−7.1705, 7.0471)2.6717158(−0.1858, 1.4740)0.974259(8.5137, 5.5096)4.5286159(−0.3142, −5.9115)2.043260(−1.5415, −7.7999)3.7548160(−8.6928, 3.0577)4.506961(6.4396, 6.4118)3.0396161(−8.9804, 1.7362)3.490462(8.1804, −3.8801)5.6087162(0.8940, 8.4429)1.826263(−9.6267, 6.2938)1.5951163(−1.6117, −7.0373)2.056964(9.7571, 7.7541)5.2335164(−6.0486, −4.1884)0.806165(−7.4305, −1.9970)6.0315165(0.5796, −2.2476)4.939466(−4.1666, 9.7131)4.442166(−3.0904, −9.0519)1.678567(8.6054, −5.4922)5.3479167(8.5203, −7.3574)2.437768(−9.7305, 1.9685)3.6987168(9.3078, −5.7472)2.945869(4.0727, −8.2039)3.3089169(−1.0789, −5.0646)6.19570(−3.8506, 8.7982)3.216170(3.4054, −9.4300)0.060871(−7.1437, 7.1571)2.9886171(−8.1767, −0.5804)2.569672(8.6707, 3.3893)4.1855172(1.0544, −1.7339)4.075373(−6.2558, −1.5341)3.2457173(−8.5217, −1.5462)0.906974(−4.9621, −2.4825)1.0398174(7.3960, 5.3860)2.800875(6.5565, −1.0670)1.6786175(7.8358, −1.8235)4.369176(−9.1507, −8.5815)2.1336176(0.7854, −0.1814)2.162977(2.3058, 3.3969)3.1581177(7.3398, 4.4719)4.886278(−6.1230, 1.4946)1.635178(3.6908, 0.0740)5.874879(−6.8841, 5.0681)2.3038179(−1.7045, 5.2675)3.517780(3.1371, −3.3084)5.5688180(−4.3713, −5.1950)3.696681(1.7032, −5.9946)0.4024181(8.7083, −8.4181)2.144582(−1.5590, 5.8996)4.9225182(3.5814, 9.0920)4.384583(−5.1026, 7.0859)3.2166183(4.9116, −6.6143)2.125284(−2.3133, −0.0468)1.2746184(−9.3055, −6.6730)5.313585(−1.9291, −8.8507)5.5929185(1.1326, 9.1880)2.765386(4.6261, −7.2169)0.2662186(2.4269, 9.3005)3.812287(7.3182, 2.3801)0.2147187(−8.2876, 9.1493)4.976988(6.0759, −3.9352)5.8295188(−9.5075, 5.6753)2.537689(8.0556, −4.2223)2.5624189(8.3507, −0.2268)6.196290(−5.7516, 8.0645)5.6345190(6.5826, 3.3623)3.794491(8.1302, 8.2009)3.6932191(7.2440, −6.8963)3.564192(3.9469, 4.9619)1.1317192(−1.3492, −6.6058)4.00593(2.0032, −8.4006)1.2813193(1.2460, −5.5328)0.03694(7.8395, 7.5786)1.4329194(−1.1003, 0.0472)5.948895(−7.5101, −8.0132)0.9963195(−2.2261, −3.0986)1.103996(8.9171, 7.1257)5.1115196(2.2383, −1.0874)1.634297(−0.0308, −6.0585)2.2197(0.0418, −8.8158)5.870298(4.7966, 9.9410)4.871198(3.8139, 2.0351)4.096399(−9.4510, −7.4800)0.9984199(5.2971, −5.0442)4.8425100(6.4788, −0.9906)1.6839200(5.0138, −0.5799)4.0817Figure 4
Routing of the 2 pairs of vehicles and UAVs.The CO2 emission of the first vehicle is 209.5249 Kg, and that of the corresponding UAV is 0.1277 Kg. The CO2 emission of the second vehicle is 218.6621 kg and that of the corresponding UAV is 0.1636 kg. Therefore, the total CO2 emission from the package delivery is 428.4783 kg. If no UAV is used, the total CO2 emission will be 467.8864 kg.The variable costs of the first pair of a vehicle and a UAV are $145.52 and $2.03, respectively, and the variable costs of the second pair of a vehicle and a UAV are $152.21 and $2.60. Therefore, the total cost of serving the 200 customers is $1302.36. The total cost will be $1388.95, if no UAV is used.


### 5.3. Extended Examples

In the section, we extend our example to 300 and 400 customers.


#### 5.3.1. The Impact of UAVs on CO2 Emissions

Table2 shows the CO2 emissions without the use of UAVs, and Table 3 shows the CO2 emissions with the use of UAVs. Table 4 shows the intuitive comparison between the results without UAVs and those with UAVs.Table 2
CO2 emissions without UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)CO2 emission of each vehicle (kg)Total emission (kg)CPU time (s)20027.177.68225.9088241.9776467.886426.5830035.655.797.29178.0174182.4284229.6897590.135575.3240044.945.416.415.06155.6471170.4556201.9631159.4279687.4937156.24Table 3
CO2 emissions with UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)CO2 emission of each vehicle (kg)CO2 emission of each UAV (kg)Total emission (kg)CPU time (s)20026.656.94209.5249218.66210.12770.1636428.478329.1130027.757.89244.1831248.59420.19000.1623493.129687.1940036.316.327.08198.8123199.1274223.07310.14850.14720.1604621.4689179.21Table 4
Comparison of CO2 emissions of using UAVs and not using UAVs.
Number of customersNumber of vehiclesNumber of reduced vehiclesCO2 emission (kg)CO2 emission reductionWithout UAVWith UAVWithout UAVWith UAV(kg)(%)VehicleUAVTotal200220467.8864428.18700.2913428.478339.40819.20300321590.1355492.77730.3523493.129697.005919.67400431687.4937621.01280.4561621.468966.024810.62It can be seen from Tables2–4 that, in the case of traffic restrictions, with the increase in the number of customers, the number of vehicles and CO2 emissions are also rising. From Table 4, we can find that when there are 300 or 400 customers, 3 and 4 vehicles are needed, respectively, when no UAVs are being used, while only 2 and 3 vehicles are needed, respectively, when UAVs are being used. We also found that using UAVs can effectively reduce CO2 emission even when there are 200 customers and the number of vehicles is not reduced.


#### 5.3.2. The Impact of UAVs on Costs

Now, we analyze the impact of using UAVs on delivery costs. Table5 shows the delivery costs without UAV, Table 6 shows the delivery costs with UAV, and Table 7 shows the comparison between the costs without UAVs and that with UAVs.Table 5
Variable costs without UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)Cost per vehicle ($)Total variable cost ($)CPU time (s)20027.177.68188.32200.63388.9526.5830035.655.797.29154.96158.05198.55511.5675.3240044.945.416.415.06142.15155.39182.69145.72625.95156.24Table 6
Variable costs with UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)Cost per vehicle ($)Cost per UAV ($)Total variable cost ($)CPU time (s)20026.656.94145.52152.212.032.60302.3629.1130027.757.89173.42177.543.022.58356.5687.1940036.316.327.08131.74135.09182.262.362.342.55456.34179.21Table 7
Comparison of the costs of using UAV and not using UAV.
Number of customersNumber of reduced vehiclesNumber of reduced vehiclesNumber of reduced vehiclesPercentage of variable cost reduction by using UAV (%)Number of reduced vehiclesPercentage of total cost reduction by using UAV (%)Without UAVWith UAVWithout UAVWith UAVWithout UAVWith UAV200220388.95302.3628.641388.951302.366.65300321511.56356.5643.472011.561356.5648.28400431625.95456.3443.012625.951956.3434.23From Tables5–7, we found that the variable cost, fixed cost, and total cost rise with the increase in the number of customers. And as per our expectation, we found from Table 7 that using UAVs can always effectively reduce variable cost and total cost, as well as fixed cost when there are 300 and 400 customers because the number of required vehicles is reduced. Our expectation is that using UAVs greatly reduces the variable cost by at least 28.64% and up to 43.47%.Tables4 and 7 show that using UAVs in package delivery reduces not only the variable cost and total costs but also CO2 emission. In other words, UAV delivery is environmentally and economically beneficial.


## 5.1. Data Parameters

The range of customer demand,qj, is 1∼10 pounds, among which 80% is no more than 5 pounds, and the position of each customer (x, y) coordinate values x and y randomly evenly generated in the interval [−10, +10]. The distance of the vehicle from customer i to j is LijV, following the Manhattan metric; the distance of the UAV from customer i to j is LijU, following the Euclidean metric; in addition, the warehouse is located at the center of the region (0,0).Referring to Chiang et al. [2], we set the other parameters as follows.The empty weight of the UAV,WU, is 55 pounds; the payload of the UAV, MU, is 5 pounds; the maximum flight distance of the UAV, LU, is 10 miles; and the average flight speed of the UAV from customer i to j, VijU, is 25 miles per hour. For roads with the second kind of traffic restriction, RijC=RijE=1.2 and RijS=0.8333.The empty weight of the vehicle (truck),WV, is 6100 pounds, its payload, MV, is 6000 pounds, and the average speed of the vehicle from customer i to j, VijV, is 25 miles per hour.The fixed cost of a pair of a vehicle and a UAV,CF, is $500, the unit route cost of a vehicle, CV, is $0.00016/pound-mile, and the unit route cost of a UAV, CU, is $0.00036364/pound-mile.Referring to Goodchild and Toy [33], the weighted average emission of vehicles, WAER, is 1.2603 Kg/pound-mile, the CO2 emission from the power generation facilities in the generation of electricity, PGFER, is 3.773 × 10−4 Kg/Wh, and the average energy requirement of UAV, AER, is 3.3333 Wh/pound-mile.


## 5.2. The Examples

Firstly, we present an example of 20 customers under the two separate kinds of traffic restrictions, in order to analyze the efficiency of UAV delivery overcoming traffic restrictions. When there are 20 customers, only one vehicle or a pair of a vehicle and a UAV is needed. Figure3(a) shows the routing of a vehicle under no traffic restriction, and Figures 3(b)–3(d) separately show the routings of a vehicle and a pair of a vehicle and a UAV under two kinds of traffic restrictions. From Figure 3(d), we can find that the routings under two kinds of traffic restrictions are the same.Figure 3
Routings under different circumstances: (a) routing of only vehicle under no traffic restriction; (b) routing of only vehicle under the first kind of traffic restriction; (c) routing of only vehicle under the second kind of traffic restriction; (d) routing of vehicle with UAV under two kinds of traffic restrictions.
(a)(b)(c)(d)The results of costs and emissions under different circumstances are as follows. When a vehicle is under no traffic restriction, the CO2 emission is 161.3032 kg, and the variable cost is $134.09. When a vehicle is under the first kind of traffic restriction, the CO2 emission is 167.4327 Kg, and the variable cost is $139.19. When a vehicle is under the second kind of traffic restriction, the CO2 emission is 166.3727 Kg, and the variable cost is $138.30. When a pair of a vehicle is under both the first and the second kinds of traffic restriction, the CO2 emissions are the same, which is 151.6250 kg, and the variable cost is the same too, which is $113.97.From Figure3(d) and the result of costs and emissions, we can find that the routing as well as cost and emission under the first kind of traffic restriction are the same as those under the second. The main reason is no matter which kind of the traffic restriction is, we always use a UAV instead of a vehicle to deliver a parcel when facing a traffic restriction.Comparing Figures3(a)–3(d), as well as the corresponding result of costs and emissions, we can obtain the following conclusions.Firstly, traffic restrictions raise the cost and emission, and first kind of traffic restriction raises them further higher. The main reason is that when a road in the routing is blocked for the delivery vehicle, it has to choose other road, while if there is only one lane available, the vehicle could still run through the road if the cost and emission are lower than those running through other roads. Therefore, the cost and emission under the second kinds of traffic restriction are higher than those under no traffic restriction but lower than those under the first kind of traffic restriction.Secondly, UAV delivery can always reduce the cost and emission no matter if there is traffic restriction and which kind of traffic restriction is. Furthermore, UAV delivery saves more cost and cuts more emission under the first kind of traffic restriction than that under the second, as the cost and emission under the first kind of traffic restriction are the highest and those under two kinds of traffic restrictions are the same.In the following section, firstly, we present the result of 200 customers under the first kind of traffic restriction. The location of the 200 customers and the weight of their demand of were randomly generated and shown in Table1. 2 cars are needed for customer service within 8 hours per day. The service times of each vehicle are 6.65 hours and 6.94 hours, respectively. The routing of the two pairs of vehicles and drones is shown in Figure 4, and the results of costs and emissions are as follows.Table 1
Location of the 200 customers and the weight of their demand.
CustomerLocationWeight of demand (pound)CustomerLocationWeight of demand (pound)1(3.5985, 3.6248)2.6464101(6.6985, −7.8311)2.46952(1.6070, 2.0104)2.0035102(−1.9531, −3.3253)4.63763(6.7836, 3.9658)4.9405103(−7.8722, −1.6863)5.39564(2.7794, 2.7118)0.4074104(5.0983, −0.4058)5.52065(5.4746, −1.4383)5.5835105(−7.3811, −9.4588)2.78916(3.7286, −5.7173)2.6218106(−9.7011, −9.1042)0.60357(−3.9446, −0.7422)0.9957107(−8.0567, 4.1457)5.72418(−6.6992, 6.3058)0.8271108(−6.8909, −9.5609)2.01719(−8.7527, −4.0516)0.7294109(−5.7863, 4.4279)6.110510(−4.4899, −5.3138)1.3197110(7.3203, 5.0068)2.323511(4.4820, −3.0529)6.0178111(8.0247, −7.8543)4.21312(1.4923, −1.8090)4.5951112(−5.0841, 1.2691)5.537813(−8.3485, 2.9984)3.7066113(3.0984, −5.5088)4.010314(3.5833, −0.3746)4.2485114(−8.4013, −2.8959)3.906615(−8.8362, 3.5809)4.8984115(−1.5125, 3.4101)4.594616(−3.4012, 2.5994)2.6331116(4.7787, 8.0351)2.548117(−3.2365, 0.3061)5.1189117(0.5036, −1.8153)4.485918(−9.3533, 3.8708)2.3671118(9.2231, −5.6816)6.126519(5.7948, 3.0688)1.6931119(−0.7709, −0.2432)0.838620(−9.3260, 2.3723)6.0346120(0.1749, 0.4201)2.517621(5.5314, −8.4922)3.9256121(6.5870, −5.7042)3.450622(−5.3389, −1.8910)4.0786122(5.1065, −4.8937)5.618423(5.0150, 1.0702)0.5409123(−6.5286, −4.1542)4.541124(5.3291, 8.8109)0.726124(4.8201, 0.3814)1.821525(5.1424, 2.8790)3.9069125(−8.0454, 2.6018)3.855326(−8.5682, −4.6079)0.2204126(5.8765, 0.1324)0.41327(2.4282, −7.1526)5.7625127(−8.9043, 2.6839)2.514728(5.3858, 3.9510)5.8148128(3.2194, 4.7627)3.514829(−3.1541, −4.4266)6.0896129(−3.9258, 0.0123)4.780830(6.9408, 0.1415)1.0124130(0.0193, 9.3367)4.643631(6.9864, 6.0812)2.2663131(−1.7297, −2.4018)5.767332(−2.2162, −7.0139)2.0512132(7.9912, 3.1555)4.044733(1.7152, 6.3001)1.4065133(−1.3450, 1.1250)2.950434(7.2730, 1.6194)0.4455134(4.6406, 2.6336)6.070435(−9.0451, −1.9021)5.6629135(−1.8672, −3.1354)3.754736(3.9235, 3.1469)5.9208136(4.4430, 7.4150)0.849137(−4.1134, 9.0178)0.7815137(−7.3532, 3.6512)0.828338(4.1405, 0.6544)5.1735138(6.8433, 5.5395)2.508139(−9.0047, 2.8633)5.2433139(8.1880, 8.9098)2.848740(3.9065, −5.4808)3.8613140(−8.6609, −7.1078)4.79841(6.8066, 6.7106)1.008141(6.5072, 1.6621)4.981542(1.4364, −6.4336)4.315142(−1.7420, 8.7680)1.447243(−8.9608, −7.4226)5.8182143(−8.7449, −4.7519)3.411144(0.6380, 1.0638)5.9632144(2.5800, 9.5467)5.510645(−5.3386, −5.5256)1.804145(−6.2282, 2.4928)3.412546(−2.7935, 0.2916)3.8737146(9.7491, 2.2369)3.690847(−3.2275, −7.6894)6.1055147(5.0762, −8.1123)3.973848(6.0844, 5.0847)6.0703148(−2.1303, 7.2874)2.636149(−5.1356, −5.9556)3.0152149(−6.2300, −6.5167)2.69950(2.6515, 7.7721)0.0072150(6.2027, −0.4723)4.616651(−0.0056, −0.1799)3.3836151(2.0571, 3.5960)5.895652(2.7594, −4.9004)5.8689152(9.7581, 5.6005)1.163253(6.6679, 7.1263)5.246153(5.2829, 1.7994)2.592554(9.8432, 3.9686)1.3284154(−2.2265, 6.9336)2.374155(−5.3313, 6.4587)1.3662155(−8.9635, 6.1876)3.386456(−7.3582, 2.8845)5.5513156(7.3399, −1.8563)2.569457(0.1812, 2.1129)3.2123157(−3.4740, −8.4158)5.95158(−7.1705, 7.0471)2.6717158(−0.1858, 1.4740)0.974259(8.5137, 5.5096)4.5286159(−0.3142, −5.9115)2.043260(−1.5415, −7.7999)3.7548160(−8.6928, 3.0577)4.506961(6.4396, 6.4118)3.0396161(−8.9804, 1.7362)3.490462(8.1804, −3.8801)5.6087162(0.8940, 8.4429)1.826263(−9.6267, 6.2938)1.5951163(−1.6117, −7.0373)2.056964(9.7571, 7.7541)5.2335164(−6.0486, −4.1884)0.806165(−7.4305, −1.9970)6.0315165(0.5796, −2.2476)4.939466(−4.1666, 9.7131)4.442166(−3.0904, −9.0519)1.678567(8.6054, −5.4922)5.3479167(8.5203, −7.3574)2.437768(−9.7305, 1.9685)3.6987168(9.3078, −5.7472)2.945869(4.0727, −8.2039)3.3089169(−1.0789, −5.0646)6.19570(−3.8506, 8.7982)3.216170(3.4054, −9.4300)0.060871(−7.1437, 7.1571)2.9886171(−8.1767, −0.5804)2.569672(8.6707, 3.3893)4.1855172(1.0544, −1.7339)4.075373(−6.2558, −1.5341)3.2457173(−8.5217, −1.5462)0.906974(−4.9621, −2.4825)1.0398174(7.3960, 5.3860)2.800875(6.5565, −1.0670)1.6786175(7.8358, −1.8235)4.369176(−9.1507, −8.5815)2.1336176(0.7854, −0.1814)2.162977(2.3058, 3.3969)3.1581177(7.3398, 4.4719)4.886278(−6.1230, 1.4946)1.635178(3.6908, 0.0740)5.874879(−6.8841, 5.0681)2.3038179(−1.7045, 5.2675)3.517780(3.1371, −3.3084)5.5688180(−4.3713, −5.1950)3.696681(1.7032, −5.9946)0.4024181(8.7083, −8.4181)2.144582(−1.5590, 5.8996)4.9225182(3.5814, 9.0920)4.384583(−5.1026, 7.0859)3.2166183(4.9116, −6.6143)2.125284(−2.3133, −0.0468)1.2746184(−9.3055, −6.6730)5.313585(−1.9291, −8.8507)5.5929185(1.1326, 9.1880)2.765386(4.6261, −7.2169)0.2662186(2.4269, 9.3005)3.812287(7.3182, 2.3801)0.2147187(−8.2876, 9.1493)4.976988(6.0759, −3.9352)5.8295188(−9.5075, 5.6753)2.537689(8.0556, −4.2223)2.5624189(8.3507, −0.2268)6.196290(−5.7516, 8.0645)5.6345190(6.5826, 3.3623)3.794491(8.1302, 8.2009)3.6932191(7.2440, −6.8963)3.564192(3.9469, 4.9619)1.1317192(−1.3492, −6.6058)4.00593(2.0032, −8.4006)1.2813193(1.2460, −5.5328)0.03694(7.8395, 7.5786)1.4329194(−1.1003, 0.0472)5.948895(−7.5101, −8.0132)0.9963195(−2.2261, −3.0986)1.103996(8.9171, 7.1257)5.1115196(2.2383, −1.0874)1.634297(−0.0308, −6.0585)2.2197(0.0418, −8.8158)5.870298(4.7966, 9.9410)4.871198(3.8139, 2.0351)4.096399(−9.4510, −7.4800)0.9984199(5.2971, −5.0442)4.8425100(6.4788, −0.9906)1.6839200(5.0138, −0.5799)4.0817Figure 4
Routing of the 2 pairs of vehicles and UAVs.The CO2 emission of the first vehicle is 209.5249 Kg, and that of the corresponding UAV is 0.1277 Kg. The CO2 emission of the second vehicle is 218.6621 kg and that of the corresponding UAV is 0.1636 kg. Therefore, the total CO2 emission from the package delivery is 428.4783 kg. If no UAV is used, the total CO2 emission will be 467.8864 kg.The variable costs of the first pair of a vehicle and a UAV are $145.52 and $2.03, respectively, and the variable costs of the second pair of a vehicle and a UAV are $152.21 and $2.60. Therefore, the total cost of serving the 200 customers is $1302.36. The total cost will be $1388.95, if no UAV is used.


## 5.3. Extended Examples

In the section, we extend our example to 300 and 400 customers.


### 5.3.1. The Impact of UAVs on CO2 Emissions

Table2 shows the CO2 emissions without the use of UAVs, and Table 3 shows the CO2 emissions with the use of UAVs. Table 4 shows the intuitive comparison between the results without UAVs and those with UAVs.Table 2
CO2 emissions without UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)CO2 emission of each vehicle (kg)Total emission (kg)CPU time (s)20027.177.68225.9088241.9776467.886426.5830035.655.797.29178.0174182.4284229.6897590.135575.3240044.945.416.415.06155.6471170.4556201.9631159.4279687.4937156.24Table 3
CO2 emissions with UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)CO2 emission of each vehicle (kg)CO2 emission of each UAV (kg)Total emission (kg)CPU time (s)20026.656.94209.5249218.66210.12770.1636428.478329.1130027.757.89244.1831248.59420.19000.1623493.129687.1940036.316.327.08198.8123199.1274223.07310.14850.14720.1604621.4689179.21Table 4
Comparison of CO2 emissions of using UAVs and not using UAVs.
Number of customersNumber of vehiclesNumber of reduced vehiclesCO2 emission (kg)CO2 emission reductionWithout UAVWith UAVWithout UAVWith UAV(kg)(%)VehicleUAVTotal200220467.8864428.18700.2913428.478339.40819.20300321590.1355492.77730.3523493.129697.005919.67400431687.4937621.01280.4561621.468966.024810.62It can be seen from Tables2–4 that, in the case of traffic restrictions, with the increase in the number of customers, the number of vehicles and CO2 emissions are also rising. From Table 4, we can find that when there are 300 or 400 customers, 3 and 4 vehicles are needed, respectively, when no UAVs are being used, while only 2 and 3 vehicles are needed, respectively, when UAVs are being used. We also found that using UAVs can effectively reduce CO2 emission even when there are 200 customers and the number of vehicles is not reduced.


### 5.3.2. The Impact of UAVs on Costs

Now, we analyze the impact of using UAVs on delivery costs. Table5 shows the delivery costs without UAV, Table 6 shows the delivery costs with UAV, and Table 7 shows the comparison between the costs without UAVs and that with UAVs.Table 5
Variable costs without UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)Cost per vehicle ($)Total variable cost ($)CPU time (s)20027.177.68188.32200.63388.9526.5830035.655.797.29154.96158.05198.55511.5675.3240044.945.416.415.06142.15155.39182.69145.72625.95156.24Table 6
Variable costs with UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)Cost per vehicle ($)Cost per UAV ($)Total variable cost ($)CPU time (s)20026.656.94145.52152.212.032.60302.3629.1130027.757.89173.42177.543.022.58356.5687.1940036.316.327.08131.74135.09182.262.362.342.55456.34179.21Table 7
Comparison of the costs of using UAV and not using UAV.
Number of customersNumber of reduced vehiclesNumber of reduced vehiclesNumber of reduced vehiclesPercentage of variable cost reduction by using UAV (%)Number of reduced vehiclesPercentage of total cost reduction by using UAV (%)Without UAVWith UAVWithout UAVWith UAVWithout UAVWith UAV200220388.95302.3628.641388.951302.366.65300321511.56356.5643.472011.561356.5648.28400431625.95456.3443.012625.951956.3434.23From Tables5–7, we found that the variable cost, fixed cost, and total cost rise with the increase in the number of customers. And as per our expectation, we found from Table 7 that using UAVs can always effectively reduce variable cost and total cost, as well as fixed cost when there are 300 and 400 customers because the number of required vehicles is reduced. Our expectation is that using UAVs greatly reduces the variable cost by at least 28.64% and up to 43.47%.Tables4 and 7 show that using UAVs in package delivery reduces not only the variable cost and total costs but also CO2 emission. In other words, UAV delivery is environmentally and economically beneficial.


## 5.3.1. The Impact of UAVs on CO2 Emissions

Table2 shows the CO2 emissions without the use of UAVs, and Table 3 shows the CO2 emissions with the use of UAVs. Table 4 shows the intuitive comparison between the results without UAVs and those with UAVs.Table 2
CO2 emissions without UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)CO2 emission of each vehicle (kg)Total emission (kg)CPU time (s)20027.177.68225.9088241.9776467.886426.5830035.655.797.29178.0174182.4284229.6897590.135575.3240044.945.416.415.06155.6471170.4556201.9631159.4279687.4937156.24Table 3
CO2 emissions with UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)CO2 emission of each vehicle (kg)CO2 emission of each UAV (kg)Total emission (kg)CPU time (s)20026.656.94209.5249218.66210.12770.1636428.478329.1130027.757.89244.1831248.59420.19000.1623493.129687.1940036.316.327.08198.8123199.1274223.07310.14850.14720.1604621.4689179.21Table 4
Comparison of CO2 emissions of using UAVs and not using UAVs.
Number of customersNumber of vehiclesNumber of reduced vehiclesCO2 emission (kg)CO2 emission reductionWithout UAVWith UAVWithout UAVWith UAV(kg)(%)VehicleUAVTotal200220467.8864428.18700.2913428.478339.40819.20300321590.1355492.77730.3523493.129697.005919.67400431687.4937621.01280.4561621.468966.024810.62It can be seen from Tables2–4 that, in the case of traffic restrictions, with the increase in the number of customers, the number of vehicles and CO2 emissions are also rising. From Table 4, we can find that when there are 300 or 400 customers, 3 and 4 vehicles are needed, respectively, when no UAVs are being used, while only 2 and 3 vehicles are needed, respectively, when UAVs are being used. We also found that using UAVs can effectively reduce CO2 emission even when there are 200 customers and the number of vehicles is not reduced.


## 5.3.2. The Impact of UAVs on Costs

Now, we analyze the impact of using UAVs on delivery costs. Table5 shows the delivery costs without UAV, Table 6 shows the delivery costs with UAV, and Table 7 shows the comparison between the costs without UAVs and that with UAVs.Table 5
Variable costs without UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)Cost per vehicle ($)Total variable cost ($)CPU time (s)20027.177.68188.32200.63388.9526.5830035.655.797.29154.96158.05198.55511.5675.3240044.945.416.415.06142.15155.39182.69145.72625.95156.24Table 6
Variable costs with UAVs.
Number of customersNumber of vehiclesService time of each vehicle (hours)Cost per vehicle ($)Cost per UAV ($)Total variable cost ($)CPU time (s)20026.656.94145.52152.212.032.60302.3629.1130027.757.89173.42177.543.022.58356.5687.1940036.316.327.08131.74135.09182.262.362.342.55456.34179.21Table 7
Comparison of the costs of using UAV and not using UAV.
Number of customersNumber of reduced vehiclesNumber of reduced vehiclesNumber of reduced vehiclesPercentage of variable cost reduction by using UAV (%)Number of reduced vehiclesPercentage of total cost reduction by using UAV (%)Without UAVWith UAVWithout UAVWith UAVWithout UAVWith UAV200220388.95302.3628.641388.951302.366.65300321511.56356.5643.472011.561356.5648.28400431625.95456.3443.012625.951956.3434.23From Tables5–7, we found that the variable cost, fixed cost, and total cost rise with the increase in the number of customers. And as per our expectation, we found from Table 7 that using UAVs can always effectively reduce variable cost and total cost, as well as fixed cost when there are 300 and 400 customers because the number of required vehicles is reduced. Our expectation is that using UAVs greatly reduces the variable cost by at least 28.64% and up to 43.47%.Tables4 and 7 show that using UAVs in package delivery reduces not only the variable cost and total costs but also CO2 emission. In other words, UAV delivery is environmentally and economically beneficial.


## 6. Discussion

In this paper, we classified traffic restrictions into two circumstances and studied the environmental and economic impact of UAV delivery under traffic restriction, and we found that UAV delivery can effectively save energy and cost. In other words, UAV delivery economically benefits not only logistics firms but also the sustainability of environmental development, which is supported by Chiang et al. [2]. We also found that, under the two kinds of traffic restrictions, UAV delivery can accomplish deliveries that cannot be carried out or are carried out at much higher costs by vehicles only, and it is more effective under the first kind of traffic restriction, which has not been studied.We would like to point out that, with the development of UAV technology, UAV delivery is more and more popularly accepted as the best resolution of “last-mile” delivery problem. Amazon, UPS, Walmart, Google, JD.com, and Alibaba started to apply UAVs in package delivery. In this paper, we proposed managers an effective model to create competitiveness by optimally coordinating vehicles and UAVs. Shifting small package delivery from trucks to UAVs greatly reduces energy consumption and costs in package delivery. As a matter of course, to realize the environmental and economic benefits of UAV delivery, firms should carefully plan and control the routing and cooperation of vehicles and UAVs. Operational decisions can effectively reduce delivery costs and CO2 emission, as well as improve the environment [34].Absolutely, the current power technology of UAVs limits the range and payload capacity of UAVs. At present, the battery used in UAVs is lithium ion battery with an energy density of around 300 Wh/kg, while the fluorinated battery can reach 2585 Wh/kg, which can greatly improve the performance of UAVs [35]. This means future developments in emerging energy sources will ease current restrictions on drone use and may provide additional cost savings that will allow UAV-vehicle delivery of last-mile packages to be further improved.


## 7. Conclusions

In this paper, we classified traffic restrictions into two circumstances and proposed a mixed-integer (0-1 linear) UAV-vehicle routing model with these two kinds of traffic restrictions to exploit the environmental aspects of the use of UAVs in logistics. A genetic algorithm was proposed to solve the model, and an experimental analysis was made to illustrate and validate our model and the algorithm. We found that, under the two kinds of traffic restrictions, delivering with UAVs can accomplish deliveries that cannot be carried out or are carried out at much higher prices by vehicles only and can effectively save costs and reduce CO2 emissions, which is environmentally friendly. Furthermore, UAV delivery saves more cost and cuts more CO2 emission under the first kind of traffic restriction than that under the second.



---

*Source: 9437605-2020-08-08.xml*