Abstract

This paper studies a coordinated system for multidrone single-truck distribution, where a truck delivers goods to a group of customers along a closed ground path with the help of a number of drones. For each delivery, the truck departs from the distribution centre with drones and all goods needed and returns back to the centre after fulfilling the delivery tasks. That is, the truck assigns these delivery tasks to several of its drones, each of which is responsible for sending goods to a different subgroup of customers in the empty air space. This study provides a new mixed-integer programming model of the routing problem with this distribution system based on urban road network. Meanwhile, a hybrid genetic algorithm and a hybrid particle swarm algorithm are designed. Experimental results show that the performance of the hybrid algorithms is better than that of the corresponding basic algorithms.

1. Introduction

Drones or UAVs (Unmanned Aerial Vehicle) could be used to transport goods such as packages, food, and medicine. Compared with traditional delivery with trucks, delivery with drones has several advantages. Firstly, a drone is able to be unmanned without well-paid drivers so that human and operational costs can be reduced. Secondly, a drone has the ability of flying at a high speed in the empty airspace to avoid ground road congestion and then improve distribution efficiency effectually. Thirdly, drones are especially suitable for the areas or situations in which land transport is not much convenient, such as remote mountains and islands and transport for emergency supplies.

However, there are also some obvious drawbacks for drone delivery. The low load capacity of a drone limits the weight or size of the goods it carries. In addition, since drones are powered by batteries, they have a very short flight radius compared to ordinary fuel trucks. Both of these disadvantages require the drone to travel frequently between central warehouse and customer demand point, which means that, after each delivery, the drone needs to return back to the warehouse to take the next goods, as shown in Figure 1. However, although a normal fuel truck can travel longer distances and carry more goods, it may travel slowly and delivery efficiency may be affected easily due to the large weight itself and road traffic jams. In addition, truck drivers have to operate during the whole driving route so that high salary will be paid to drivers, which means high human cost. In contrast, drones are able to save this kind of cost by delivering automatically without being operated during the whole time. Besides, electric cost of drones is much cheaper than fuel cost of trucks, so truck expense is more in every mile, especially when customer points are very scattered.

Figure 1 

A distribution system that only uses a drone. (i) The drone delivers to Customer 1 from the distribution centre and returns to it. (ii) The drone delivers to Customer 2 from the distribution centre and then returns back. (iii) The drone delivers to another customer from the distribution centre and then returns back.

A solution to the limited flight distance and load capacity of the drone is to allow the drone to collaborate with other types of delivery tools to complete the delivery tasks. Figure 2 shows a collaborative distribution system that uses a drone and a truck. The truck loads the goods for a group of customers at the distribution centre in advance, and it also needs to carry a drone. After leaving from the distribution centre, the truck is ought to move along a closed path on the ground. When the truck travels to a location near a customer demand point, it sends the drone carrying goods in the direction of this customer, while the truck continues to move on its intended path. The drone carries the goods along a straight line to the customer point, then returns back to the truck in another straight line after delivering, gets ready to pick-up new goods, and delivers them to the next customer. For every customer demand point, the drone performs the same actions as mentioned above so that the goods are able to be delivered to all the customers. At last, the truck returns to its original distribution centre.

Figure 2 

A collaborative distribution system that uses a drone and a truck. (i) The truck is loading goods from the distribution centre. (ii) When the truck approaches a customer, it sends a drone to deliver goods to the customer, and the truck continues to move. (iii) The drone returns to the truck and prepares the next customer's delivery.

At present, researchers have also begun to study the optimization of the collaborative distribution system where trucks and drones collaborate to complete delivery tasks, rather than simply relying on single delivery tools. Murray and Raj [1] studied a new extension of the “flying sidekick traveling salesman problem,” where a truck works with multiple various drones to deliver packages [1]. They described the problem as a mixed-integer linear programming and proposed a method by using three subproblem sequences. Zhen et al. [2] modelled and solved a nonlinear MIP model for different instance sizes with two efficient algorithms, which are SOCP-based algorithm and dynamic linearization algorithm to handle nonconvex constraints and the model under different scales of route networks, respectively [2]. Bai et al. [3] studied the efficient routing problem for precedence-constrained package delivery for heterogeneous vehicles, where a truck and a microdrone are utilized to deliver packages to dispersed customers [3]. Trucks are restricted to the street network, while microdrones are available for the last-mile delivery with the restriction of their capacity and flight range. The goal of the problem is to minimize the delivery time in the priority constraint. The authors combined topological sorting technology and proposed several task distribution algorithms. Poikonen and Golden considered the k-multivisit drone routing problem, where every drone can take off from a truck and deliver one or more packages to customers [4]. Each drone can return to a truck to replace a battery or charge, take a new package, and then fly again to another customer demand point. de Freitas and Penna studied a flying sidekick traveling salesman problem, where the drone takes off from the truck, delivers the package to the customer, and returns to the truck in the third place [5]. When the drone flies, the truck delivers packages to other customers as long as the drone has enough batteries to hover in the air. Zhen imposed both a stochastic programming formulation that can cope with arbitrary probability distribution of operation time deviation and a robust formulation that is applicable to situation in which limited information about probability distributions is available [6]. Ha et al. studied the Traveling Salesman Problem with Drone (TSP-D) and proposed a hybrid genetic search with dynamic population management and adaptive diversity control based on a split algorithm to solve the completion time for the truck and drone [7]. Results show that the proposed algorithm outperforms existing methods in terms of solution quality and improves best known solutions found in the literature. Zhen et al. developed a column generation solution approach on a set partitioning-based reformulation of the original model [8]. Results show that this simple but practical solution approach can optimally solve the daily berth planning problem instances with reasonable and acceptable for the real-world applications. Liu presented a mixed-integer programming (MIP) model for on-demand meal delivery using drones and proposed an optimization-driven, progressive algorithm for online fleet dispatch operations [9]. Schermer et al. proposed an extension of the VRPD that is called the Vehicle Routing Problem with Drones and En Route Operations (VRPDERO), where drones may not only be launched and retrieved at vertices but also on some discrete points are located on each arc [10]. They also studied the Vehicle Routing Problem with Drones (VRPD), where given a fleet of trucks, each truck carries a given number of drones, and the objective consists in designing feasible routes and drone operations such that all customers are served and minimal make span is achieved [11]. Jeong et al. extended the previous vehicle routing models to the hybrid delivery systems by considering two important practical issues: the effect of parcel weight on drone energy consumption and restricted flying areas [12]. Zhen considered a combination of probabilistic and physics-based models for truck interruptions [13]. It enables to exactly evaluate link travel time, which then acts as the basis for proposing a mixed-integer programming model that minimizes the total expected travel time of moving containers. Wang et al. studied the traveling salesman problem with drone (TSP-D), which aims to find the coordinated routes of a drone and a truck to serve a list of customers. In practice, managers sometimes intend to attain a compromise between operational cost and completion time. Therefore, this article addresses a biobjective TSP-D considering both objectives [14]. Sacramento et al. published a paper where a mathematical model is formulated, defining a problem similar to the Flying Sidekick Traveling Salesman Problem, but for the capacitated multiple-truck case with time limit constraints and minimizing cost as objective function [15]. Karak and Abdelghany presented a mathematical formulation and efficient solution methodology for the hybrid vehicle-drone routing problem (HVDRP) for pick-up and delivery services. The problem is formulated as a mixed-integer program, which minimizes the vehicle and drone routing cost to serve all customers [16]. Javadi et al. extended the traveling Repairman Problem (TRP) by assuming a single truck which can stop at customer locations and launch drones multiple times for each stop location to serve customers. They also developed an efficient hybrid Tabu Search-Simulated Annealing algorithm to solve the problem [17]. Moshref-Javadi et al. presented a mathematical formulation and a heuristic solution approach for the optimal planning of delivery routes. They formulated the problem as a Mixed-Integer Linear Programming model and proposed an efficient Truck and Drone Routing Algorithm (TDRA) [18]. Gonzalez-R et al. proposed an iterated greedy heuristic based on the iterative process of destruction and reconstruction of solutions. This process is orchestrated by a global optimization scheme using a simulated annealing algorithm. The obtained results are quite promising even for large-size scenarios [19]. Chang and Lee focused on finding an effective delivery route for trucks carrying drones and proposed a new approach on a nonlinear programming model to find shift-weights that move the centres of clusters to make wider drone-delivery areas along shorter truck route after initial K-means clustering and TSP (Traveling Salesman Problem) modelling [20]. Carlsson and Song proved that the improvement in efficiency is related to the square root of the ratio of the speeds of the truck and the UAV by combining a theoretical analysis in the Euclidean plane with real-time numerical simulations on a road network [21]. Xia et al. studied how drones can be scheduled to monitor the sailing vessels in ECAs. They also modelled the dynamics of each sailing vessel using a real-time location function which allows to approximately represent the problem on a time-expanded network. Furthermore, they developed a Lagrangian relaxation-based method to obtain near-optimal solutions [22].

The contributions to this article are as follows. (1) This paper delved into the issue of collaborative distribution route planning for drones and trucks based on urban road networks. Although there is a lot of research on co-distribution in recent years, research on co-distribution path planning based on urban road networks is not enough, and our paper filled this gap. (2) This problem involves many decision variables, including determining the routes of trucks and multiple drones and distributing tasks between multiple drones. A number of constraints are also involved in the issue, such as the closeness of the route of the truck, the cargo capacity of the truck, the cruising ability of a drone, and the sequential relationship between all launch sites and landing sites for one drone. (3) For megacities such as Shanghai, London, and Paris, the urban road network is particularly complex and large with a particularly large number of customers. It is apparent that, with the increasing size of the problem and the growing number of decision variables, performance of the algorithm gets worse. Thus, we designed a hybrid algorithm that considers the efficiency and quality of the solution and achieves better computational results at a faster rate.

2. Problem

In this section, we describe a coordinated system for multidrone single-truck distribution and then illustrate the routing problem of the coordinated system on urban road network.

2.1. Multidrone Single-Truck Distribution System

A coordinated system for multidrone single-truck distribution means that a truck works with multiple drones to deliver to a set of determined customers. That is, the truck assigns these delivery tasks to several of its drones, each of which is responsible for sending goods to a different subgroup of customers in the empty air space. In particular, a truck loads goods and a number of drones from the distribution centre at first. After leaving the distribution centre, the truck travels along a closed ground path. When approaching a customer demand point, the truck releases one of the remaining drones on it and continues to move forward on its intended path, while the drone carries goods to the customer point along a straight line and returns to the truck along another straight line after delivery. When the truck approaches the next customer demand point, the truck releases another drone from the remaining drones, the truck goes on to move forward, and the drone executes the same operation as mentioned above before returning to the truck. At this point, some of previously released drones may still be on delivery to customers and have not yet returned to the truck. If there are no extra drones on the truck, the truck has to wait for a drone to come back and then launch it to another customer point.

2.2. Collaborative Distribution System on Urban Road Network

A city road network refers to a network of all roads in the city. It can be expressed as , where is a collection of all road intersections and is a collection of the whole roads. In , the number of intersections is , and , where the intersection of number 0 represents the distribution centre. The coordinate of Intersection is

A directed line from Intersection to Intersection is defined as follows:

Figure 3 is an example of an urban road network. This network is composed of 32 intersections, where the distribution centre is Intersection 0, and the other intersections are numbered as 1, 2, …, 30, and 31. The distribution centre 0 requires to deliver goods to six customers, a, b, c, d, e, and f, respectively. A truck will start from the distribution centre to deliver those six customers in turn and then return to the distribution centre. There are many options for this closed ground route driven by the truck. Here are two options:

Figure 3 

Example of an urban road network.

Figure 4 shows the route of the truck and the flight route of multiple drones, where the black thick arrow represents the route of the truck, and the red and green light arrows represent the flight routes of Drone 1 and Drone 2, respectively. In each delivery, Drone 1 moves in a red line to a customer and then returns back to the truck in another red line. The process is the same as Drone 2 while the colour of line is green. It is obvious to see that this collaborative distribution system has improved the delivery efficiency to a great extent.

Figure 4 

The movement route of a truck and drones.

The following is the procedure of collaborative distribution in Figure 4. The truck starts from the distribution centre and moves along a closed ground route indicated with the black thick arrow. First of all, it travels along the line 0⟶11⟶10⟶9. When the truck is gradually close to intersection 9, distance between the truck and Customer a decreases. The truck launches Drone 1 along with the red arrow to Customer a. After sending goods to Customer a, Drone 1 goes back along another red arrow to the truck which is moving forward and preparing the next delivery at the same time. The truck continues to drive along line 9⟶16⟶17. When traveling to Intersection 16, the truck is getting closer to Customer c, and it launches Drone 1 to deliver to Customer c once again. Similarly, the drone moves along the red arrow, performs the same set of actions as mentioned above, and finally returns to the truck which is moving forward after Customer c's delivery is completed. In the next stage, the truck continues to drive along the line 17⟶22⟶23, and when approaching Intersection 22, the truck sends Drone 1 to Customer d along the red arrow. Noted that Drone 1 does not return to the truck when the truck is close to Intersection 23, so Drone 2 is sent to Customer f along the green arrow. Then, the truck continues to drive along the line 23⟶24⟶18⟶0, and when it is near to Intersection 24, it sends Drone 2 along the green arrow to Customer e. Since Drone 2 has not returned to the truck when it is approaching to Intersection 18, the truck needs to send Drone 1 instead along the red arrow to Customer b. At this moment, all customers have been delivered and the truck could return to the distribution centre.

2.3. Notational Conventions

Before discussing the mathematical model, it is necessary to illustrate some mathematical symbols involved. Here are some of the main mathematical symbols. : we use Euclidean distance to denote travel distance according to Carlsson [21] and Agatz [23], for any two points and : : the Euclidean plane. : the number of customers to visit in the Euclidean plane. : the collection of the customers that need to be served, . : customer number, . : the coordinates of the position of Customer , . : the amount of demand for Customer . : the total number of drones on a truck. : the collection of all the drones, . : drone number, . : load weight of the truck. : closed transport route of the truck. : urban road network . : the total number of the road intersections. : a collection of all road intersections, . : a directed line from Intersection to Intersection , , called a directed road. : a binary indicating whether a directed line segment is present in the urban road network . Note that there is no directed line segment from any intersection to itself, that is to say, : : a collection of all the directed roads, . : The coordinates of the road intersections, . : speed of the truck. : speed of each drone, . : the battery life of the drone , determining its maximum endurance.

2.4. Decision Variables

The decision variables for the routing problem of collaborative distribution system are listed as follows.  binary equals one if the directed line segment belongs to a transport closed route , while u, v = 0, 1, …, n; otherwise, equals zero:  binary equals one if Customer is delivered by Drone ; otherwise, equals zero.  “Launch site” of the drone released from the truck to Customer .  “Landing place” of the drone which returns onto the truck from Customer .

Note that each Customer corresponds to a Launching Point and a Landing Point for the drone. The drone starts from Launching Site , sends goods to Customer along a straight line, and then returns from Customer along another straight line to Landing Site on the truck. The details are shown in Figure 5, where the thick black line indicates the route of the truck, the blue dotted lines indicate the flight path of the first drone, the pink dotted lines indicate the flight path of the second drone, and the green one indicates the flight path of the third drone. and are, respectively, the launching point and landing point of the first drone for delivery of the customer point ; and are that of customer point ; and are the launching point and landing point of the second drone for delivery of the customer point ; and are the launching point and landing point of the third drone for delivery of customer point .

Figure 5 

Launching sites and landing sites of drones.

2.5. Routing Problem of the Collaborative Distribution System

The routing problem of the coordinated system for multidrone single-truck distribution can be defined as follows. Given a city road network , the distribution centre is numbered as 0 while each customer is numbered as ; the location coordinates and demand quantity of all the customers are known. In addition, some assumptions are mentioned as follows to complete our problem:(i)Only one truck is provided, but is equipped with a certain number of drones(ii)A customer demand point can only be delivered by only one drone, while a drone is able to deliver to none or more customer demand points(iii)The weight of cargo for every customer is no more than the capacity of a drone(iv)All drones have the same speed(v)Each drone can only transport only one good each time(vi)The truck is just a platform for drones to launch and land, while itself does not visit customer demand points(vii)Drones could land on or fly off the truck no matter if the truck is motionless or locomotor(viii)The number of drones on each truck is fixed and the load weight of the truck excludes the total weight of drones on it, which means the weight of drones is negligible

Our problem aims to find out two sorts of routes in the shortest distribution time: one is the closed transport route in which starting and terminal points are both distribution centre in the whole delivery of the truck, and the other is flight path of each drone on the truck, which means that launching and landing points of all drones are ought to be decided in this problem.

3. Mathematical Model

Some of relevant intermediate variables will be illustrated before recommending the mathematical model.

3.1. Intermediate Variables Related to the Closed Route of the Truck

We define a symbol that indicates the mileage of the truck from the distribution centre to Intersection . Obviously, .

If the intersection is not on the transport closed route , then there is no intersection that makes . Set when .

If the intersection is on the transport closed route , then there is only one intersection that makes ; in other words, . Because the length of the directed line segment isso

In a word,

3.2. Intermediate Variables Related to Launching and Landing Sites of Drones

In order to simplify the mathematical model, the following decision variables are utilized directly to substitute launching site and landing site . : a binary indicating whether the launching point is located on the directed line segment : : a decimal indicating dividing the distance between and its start point by the length of line segment , so : : an integer indicating the line segment of in which landing point belongs to : a decimal indicating dividing the distance between and its start point by the length of line segment , so :

Depending on the decision variables above, coordinates of the launching point and landing point could be calculated as follows: : the x-coordinate of the launching point:  : the y-coordinate of the launching point:  : the x-coordinate of the landing point:  : the y-coordinate of the landing point:

Define a symbol that indicates the mileage of the truck from the distribution centre to launching point :

That is,

Define a symbol that indicates the mileage of the truck from the distribution centre to landing point :

That is,

The flight mileage of the drone when delivering to Customer can also be calculated as follows:

3.3. Mathematical Model and Its Description

subject to

Here, , , , , , and are decision variables with and . In particular, , , , or is an integer for 0 or 1, while or is a decimal between 0 and 1, that is, and .

The objective function (19) minimizes the total delivery time. We do not simply add the truck time and drone time to obtain the whole completion time during delivery tasks. This is because these two vehicles need to be synchronized at the drone’s landing point, which means that the truck have to stop and wait for the drone if it arrives to the landing point early, while the drone is also ought to wait until the truck comes if the drone arrives to the landing point early.

In function (19), the first term represents the total time which the truck spends to wait for the drone or the drone spends to wait for the truck in all pairs of launching point and landing point . In this term, shows the time which the truck takes from a launch point to the corresponding landing point, and indicates the total time that the drone takes to leave the truck from the launching point, reach the customer demand point, and then return to the truck at the landing point. The second term represents the total time of the truck on the closed loop when the drone and the truck are fully synchronized for each delivery. In other words, flight time of the drone is completely equal to the time of the truck for each pair of take-off points and landing points.

Constraint (20) indicates that there is no directed line segment from the intersection to itself. Constraint (21) says if , that is to say , then . Otherwise, if , then or . Constraint (22)shows that the number of input lines of any Intersection cannot exceed 1. Constraint (23) ensures that the quantity of input lines of any Intersection is 1 at most. These two constraints are used to guarantee that any line segment is passed by the truck only once. Constraint (24) means that an output or input line segment of an intersection is allocated to at most one input or output line segment or not. These constraints enable a closed loop of the truck, and there is no need for it to travel through all intersections in the road network. Instead, the truck just requires to go through part of the intersections, starting from and returning back to the distribution centre. Constraint (25) indicates that the total distance of the truck driving from the distribution centre to the distribution centre is 0. Constraint (26) represents the total distance of the truck driving from the distribution centre to the intersection , . Constraint (27) says if , that is to say , then . Otherwise, if , then or . Constraint (28) says if , that is to say , then . Otherwise, if , then or . Constraint (29) indicates that a launching point can only be on one directed line segment. Constraint (30) indicates that a landing point can only be on one directed line segment. Constraints (31) and (32) represent the x-coordinate and y-coordinate of the launching point, respectively, while Constraints (33) and (34) are the x-coordinate and y-coordinate of the landing point, respectively. Constraints (35) and (36) show the mileage of the truck from distribution centre to launching point and to landing point separately. Constraint (37) imposes that, for the drone of each customer, the launching point corresponded on the truck trajectory is earlier than the landing point. Constraint (38) guarantees that a customer demand point is allocated to only one drone, while a drone could deliver to none or some customer demand points. Constraint (39) states the relationship between endurance mileage of a drone and its battery life. Constraints (40) and (41) ensure that for any Customer i and Customer j which are delivered by the same drone, both launching point and landing point , are either earlier than launching point or later than landing point , which means that and are not able to locate between and on the truck trajectory. Constraint (42) imposes that the gross weight of goods loading on the truck is no more than the maximum carrying weight .

Since our problem is a new extension of the Vehicle Routing Problem (VRP) which has been shown to be NP-hard, the problem in this paper is also NP-hard and the model is nonlinear. In addition, CPLEX is always used to solve small-scale vehicle routing problem, but not to large-scale vehicle routing problem due to the fact of low solution efficiency. Therefore, we utilize heuristic algorithm in this paper for large-scale problem.

4. Algorithm

With the increasing number of decision variables, problem scale, and search space of genetic algorithm (GA) or particle swarm optimization (PSO) for a complex problem, performances of algorithm is impacted negatively due to the slow convergence rate of algorithm and large amount of iterations.

In order to improve the convergence speed, this paper designs a hybrid genetic algorithm (h-GA) and a hybrid particle swarm optimization (h-PSO). The basic idea of these two hybrid algorithms is to decompose a problem into a master problem and a child problem and just encode decision variables in the master problem as a chromosome or a particle instead of encoding them all. Decision variables of the child problem are solved by a nonlinear programming algorithm or a heuristic method. This way greatly reduces the search space of genetic algorithm or particle swarm algorithm and then enhances the convergence rate to a great extent.

4.1. Hybrid Genetic Algorithm

Hybrid genetic algorithm breaks down this problem into a master problem and a child problem. The master problem is to device the best closed transport route, which the truck along will start from and return to the distribution centre after delivery tasks are completed.

According to the closed route solved by the master problem, the child problem is utilized for every customer to plan a best flight route, which the drone flies along from the truck to the customer and returns back to the moving truck after present delivering. As the demand and coordinate location of each customer is known, this child problem is equivalent to finding out the best launching point and landing point of a drone for each customer on the known closed transport route of the truck.

4.1.1. Genetic Algorithm of the Master Problem

In this section, the master problem is solved by genetic algorithm. Figure 6 shows the flowchart of this algorithm.

Figure 6 

Genetic algorithm of the master problem.

The main steps of the genetic algorithm are encoding and decoding of chromosomes and the fitness function. Please refer to relevant documents for details of other steps.

(1) Encoding and Decoding of Chromosomes of the Master Problem. Because the master problem is simply solving a best closed transport line of the truck, the only decision variable needs to be encoded into chromosomes is , where . Here, we adopt natural number encoding for . On the contrary, other decision variables such as , , , , and do not need to be encoded into chromosomes, while will be solved by the child problem. As shown in Figure 7, genetic chromosomes just contain , which is either 0 or 1. In this way, search space of the genetic algorithm is drastically reduced, and decision variables in the master problem can be solved in short time. For instance, if the closed loop of the truck is asthen, in the chromosomes, the values of genes , , , , , , , , , and are 1, and the value of the rest genes are 0.

Figure 7 

Chromosome of the master problem.

For a chromosome, if the line segments of its related genes with value 1 can form a closed loop from distribution centre and back to it, the chromosome corresponds to a feasible solution of the master problem.

(2) The Fitness Function of the Master Problem. The fitness function shows whether the solution to the master problem is good or not, and the degree is depended on the following factors. The first factor is to see if this solution to the master problem satisfies the constraints; the second one is to see if the solution to its related child problem meets the constraints; and the third one is to see whether the target value of its child problem is good or not. Specifically, the solution to the master problem requires to be a closure of the path of the truck.

 Then, the best flight route of the drone for each customer and target value of child problem are figured out based on the solution of the master problem.

The solution in child problem needs to consider the following conditions: (1) sequence of departure points; (2) landing points of the drone; (3) maximum endurance of the drone; (4) total weight of demand goods, which is at most the maximum load of the truck.

(3) Random Generation of Initial Populations. Each gene is randomly generated within its range of values so as to get a random chromosome. Repeat the above operation until a certain number of chromosomes are produced.

(4) Replication of Chromosomes and Selection of the Elite. Copy a certain number of chromosomes into the next generation with a roulette method. And for each generation of chromosomes, choose the best fitness individuals depending on a certain selection rate to replace the lowest one for the next generation.

(5) Crossover Operation. In the last generation, select two individuals randomly as fathers, randomly identify an intersection, and then perform the crossover operation.

(6) Mutation Operation. Randomly select mutation points in the chromosome, and randomly generate a value for the genes at these mutation points.

(7) Construct New Population. Production of individuals in a new population comes in three approaches: (1) produce new individuals by crossover operation; (2) generate them by mutation operation; (3) select 10% of the optimal fitness individuals of the original population.

(8) Termination Condition. If the number of iterations in the algorithm reaches the specified maximum, terminate the algorithm and output the optimal solutions. Otherwise, return to Step 4 and continue to perform the genetic operations above.

4.1.2. Heuristic Method of the Child Problem

This child problem enables to figure out the best take-off point and landing point for each customer demand point on a known closed transport line of the truck. Since the child problem is nonlinear programming, it can be solved by a heuristic method. Figure 8 shows the details of the heuristic.(1)Determine the order of assigned customer points. First, the child problem needs to determine the order of allocated customer points which the drone deliveries. If the total number of customers is n, the amount of all possible delivery orders is . The calculated amount will be enormous if all possible delivery orders are solved with the child problem. Thus, a construction method is used here to get a better delivery order in a shorter time. Given the closed transport route of the truck, solution of this construction method is as follows. For each customer demand point , calculate point which is the closest to on Route . Then, according to the sequence of these points on Route , determine the order of assigned customer points. Specifically, if is earlier than on Route , then corresponding demand point in the delivery order is arranged before . On the contrary, if is later than , is arranged after .(2)Given the truck route and delivery order as above, plan the drone route When the truck closed route and customer delivery order are solved, this child problem is transformed into a convex optimization problem for the take-off point and landing point of the drone and can also be solved by a structured method as follows. When solving customer’s delivery order following the above step, point which is the closest to each point on Route is calculated. The corresponding decision variables for the delivery of a drone to the customer are and . The truck and a drone are ought to avoid waiting for each other at a landing point in order to guarantee time effect. In other words, the best case is a drone and the truck arrives at the landing point at the same time. Thus, the total delivery time of the drone is equal to that of the truck: Assuming that take-off point and landing point are symmetric about the line segment from customer point to point , the flying distance of the drone from to is equal to that from to , and the driving distance of the truck from to is equal to that from to . So, an approximation can be derived: Assume that the line segment from customer point to point is perpendicular to the line segment from to . The approximate solution of the take-off point and landing point can be figured out according to the Pythagorean theorem, respectively.(3)Solve decision variable according to the order of delivery and the drone route. After the delivery order and drone route are solved by the above steps 1 and 2, decision variables can be solved in order by the sequence loop traversal method as follows. For each drone k, choose each unvisited customer demand point in turn according to the solved delivery order. If a customer point is the ith one for Drone k to visit successfully, then set . Otherwise, if launching point of the ith customer point is before the landing point of the customer point, set and turn to consider the next customer point until all points are delivered.