# Analysis of an Automated Vehicle Routing Problem in Logistics considering Path Interruption

Authors: Yong Zhang; Lei Shi; Jing Chen; Xuefeng Li
Journal: Journal of Advanced Transportation
 (2017)

---

## Abstract

The application of automated vehicles in logistics can efficiently reduce the cost of logistics and reduce the potential risks in the last mile. Considering the path restriction in the initial stage of the application of automated vehicles in logistics, the conventional model for a vehicle routing problem (VRP) is modified. Thus, the automated vehicle routing problem with time windows (AVRPTW) model considering path interruption is established. Additionally, an improved particle swarm optimisation (PSO) algorithm is designed to solve this problem. Finally, a case study is undertaken to test the validity of the model and the algorithm. Four automated vehicles are designated to execute all delivery tasks required by 25 stores. Capacities of all of the automated vehicles are almost fully utilised. It is of considerable significance for the promotion of automated vehicles in last-mile situations to develop such research into real problems arising in the initial period.


---

## Body

## 1. Introduction

Due to fierce competition in the global market, an intelligent, efficient logistics system has become one of the most critical assets needed for an enterprise to stand out from others. The rapid development of the urbanisation and car ownership has produced big impacts on the transport systems, leading to challenges to many transport topics/areas, including demand management [1–3], traffic signalization [4, 5], traffic safety [6, 7], and public transport [8–11]. Modern logistics industry asserts a strict claim to high efficiency.Compared to conventional modes of transport, automated vehicles are proved to offer a huge improvement in service efficiency and quality. In dense urban cities, labour costs are high. Manpower is one of the significant elements in any logistics system, which undoubtedly leads to the high level of inputs into the logistics system. Additionally, fatigue driving and bad driving behaviours cause potential risks for the transportation safety, but automated vehicles are free from these worries. To be concluded, the application of automated vehicles in logistics can efficiently cut down the operation costs, improve the efficiency of freightage, and eliminate the potential risks caused by human beings. In this regard, the methodology of automated vehicle routing planning in logistics is proposed. Implementing automated equipment in the last mile can largely cut down the cost of the logistics and render the logistics system “smart” or “intelligent,” which thus meets consumer expectations; however, unmanned equipment is still on the test bed. Defects and instabilities in the technology limit the scale of the application of automated vehicles in the initial period. According to a report of public views on automated vehicles from the American Automobile Association (AAA), 75% of interviewees hold pessimistic views on unmanned technology: 54% of drivers maintain that automated vehicles on real roads must increase the potential risk of a traffic accident and intensify public anxiety. Thus, selectively opening or prohibiting sections of road is an essential stage in the promotion of the use of automated vehicles.The conventional VRP model is modified by considering a traffic interruption. The partial blocking of roads is considered as a constraint in the proposed model. This work focuses on the automated vehicle routing problem with time windows (AVRPTW) model of single distribution centre-mass client mode. The purpose of the distribution centre is to plan the driving circuit and satisfy the demands of all of the consumers: the objective function is to minimise the operational cost. Thus, we examine the similarity between the conventional VRP and the AVRPTW. The greatest difference is that AVRP model is faced with stricter constraint conditions, such as path interruption. Additionally, cost structure in VRP and AVRP is different: cost in the VRP model consists of dispatching and human costs. However, these costs are simplified in the AVRPTW model. Based on the fundamental particle swarm optimisation (PSO) algorithm, an improved PSO algorithm is designed by modifying the inertia weights and learning factors to improve the performance thereof. The AVRPTW model and the improved PSO algorithm are tested in a case study: their validity and practicability are certified. Additionally, the application of automated vehicles in logistics is proved to offer a huge improvement in efficiency and other benefits.The paper is organised as follows: Section2 summarises the literature review, in Section 3, the AVRPTW model considering a traffic interruption is formulated, Section 4 introduces the improved PSO algorithm, the AVRPTW model and improved PSO algorithm are applied in a case study in Section 5, and Section 6 draws the conclusions.


## 2. Literature Review

The application of unmanned vehicle technology and vehicle networking technology can effectively reduce the severity of traffic congestion, traffic accidents, and environmental pollution. Some scholars applied this technology to some fields of transportation. For example, unmanned vehicle technology was applied in intersection control to ensure the safety and efficiency thereof [12, 13]; however, as for last-mile logistics distribution, transportation, and other areas, the application of unmanned vehicles has not yet started.Jozefowiez et al. [14] analysed a multiobjective path optimisation problem by comparing the definition of each problem, objectives, and algorithms [14]. Baldacci et al. [15] reviewed recent research into the exact algorithms for the VRP with capacity and time window constraints [15]. Nagy and Salhi [16] proposed an improved heuristic algorithm to solve the pick-up and delivery VRP with multiple depots [16].In Europe, as a source of air pollution, transportation activities in the logistics sector activities account for 24% of all European emissions. In recent decades, more environmental problems related to economic activities have been transferred to the logistics field. Based on this background, a multiobjective vehicle routing problem was proposed [17]. Michallet et al. [18] studied a periodic VRP with time constraints in high value goods transportation. Archetti et al. [19] considered enterprises that not only used company drivers and delivery vehicles, but also allowed others willing to participate in one-way distribution runs, close to their delivery destination, so as to reduce the cost of distribution [19].Considering a multivehicle scheduling problem, an example of a distribution enterprise was studied by Coelho et al. [20]. The problem consisted of a number of service constraints for specific customer requirements, considering fixed vehicle, and transportation costs. An adaptive search procedure with a local search and an improved process using a variable neighbourhood were designed. Coelho et al. [20] used this algorithm to derive an optimal transportation plan and reduced the cost thereof [20]. The path planning of an unmanned aerial vehicle (UCAV) was a complex global optimisation problem, which could be used to search for the optimal path on the field of combat. Zhang et al. [21] proposed a new heuristic algorithm to solve the wolf UCAV 2D path planning problem. Then, UCAV could avoid the threat by connecting the two-dimensional coordinates of the node, achieving the lowest cost while finding a safe path. The simulation results showed that the proposed method considered the quality, speed, and ultimate stability in path planning [21]. An energy-efficient path planning problem was developed: the experimental results showed that the method could save energy effectively [22]. In the last-mile logistics considering the use of electric cars as a means of conveyance, Klumpp [23] studied the logistics cost and undertook a sustainability evaluation [23]. A mixed integer linear programming optimisation algorithm was proposed to solve the VRP, to validate the algorithm. A periodic vehicle routing problem (PVRP) was used as a case study by Cacchiani et al. [24]. The results verified the effectiveness of the proposed algorithm and were the best known solutions [24]. The hybrid particle swarm optimisation algorithm was used to solve the conventional VRP [25].VRP is an NP-hard problem, and a heuristic algorithm is usually needed to solve this kind of problem. Chen et al. [26] used system dynamics to modify the particle and established the adaptive particle swarm algorithm to solve the set of multiproduct material transport path optimisation problems. To minimise the total cost, corresponding constraints were set and this model algorithm was used to solve various problem instances, which further verified the effectiveness of the algorithm [26]. In order to generate an effective solution way to fix the VRPSPD, which was known as an NP-hard problem, Goksal et al. [27] proposed an approach which was a solution based on particle swarm optimisation algorithm. At the same time, an annealing strategy was used to maintain the diversity of population. The effectiveness of the algorithm was verified by experiments, and the optimal solution was obtained by using the improved ant colony algorithm [27]. Combined with a real-life application, Goel and Gruhn [28] studied the rich VRP. GVRP was a solution in the highly restricted search space, which was suggested to improve the iterative search process by changing the neighbourhood structure [28]. Although 2L-VRPB was often applied, few researchers studied this method. Based on this fact, Liu [29] proposed a hybrid metaheuristic algorithm to solve the problem and verified the effectiveness of the algorithm experimentally [29]. Cinar et al. [30] designed a finite time, two-phase, cumulative algorithm to solve the VRP and proposed a fast, easy way to implement their CumVRP-LD algorithm. The algorithm could change in response to changes in the actual situation; the high response speed was necessary to develop the system [30].Yu et al. [31] proposed an improved ant colony algorithm to solve the VRP. The effectiveness of the algorithm was verified by the data of fourteen datum points [31]. Pisinger and Ropke [32] proposed a general heuristic algorithm to solve the following five kinds of VRP: VRP with time windows, capacity-limited VRP, multidepot VRP, site-dependent VRP, and open VRP [32]. Christian [33] established a relatively simple, but effective, hybrid genetic algorithm. In terms of the average cost, the performance of the hybrid genetic algorithm was better than that of most heuristic algorithms. It was proved that this algorithm was the best solution in 20 large cases studies [33]. Schyns [34] used the ant colony algorithm to solve the VRP including mixed batch arrival and dynamic capacity with time windows [34]. The ant colony algorithm is a metaheuristic algorithm, which was used to solve the VRP. Bell and McMullen [35] compared the ant colony algorithm with the TS search algorithm and genetic algorithm and found that the ant colony algorithm required less computation time [35]. Gribkovskaia et al. [36] proposed a mixed integer programming model (SVRPPD), which was compared with the classical algorithm and the improved heuristic algorithm. The results showed that the best solution was often generated by the heuristic algorithm [36]. Kumar et al. [37] used a self-learning particle swarm optimisation (PSO) algorithm to solve the multiobjective model of production and pollution, and the total operating cost was minimised [37].Previous studies are limited with regard to their application in practice. To sum up, contributions of this paper are two twofold: (i) AVRPTW model considering a traffic interruption is built for the consideration of the policy restriction when applying automated vehicles in last-mile situations in the initial stage; (ii) an improved PSO algorithm is developed for the AVRPTW model with path interruption.


## 3. Mathematical Model for AVRPTW

### 3.1. Problem Statement

The automated vehicle routing problem (AVRP) model can be described as follows, and the model is shown in Figure1: the total number of the consumers in one certain area is N; the distribution centre has to distribute a certain amount of goods to certain consumers; presently, this distribution centre is preparing to put in several automated vehicles to carry out the assignment. In the initial stage of the application of the automated vehicle, considering the restrictions imposed by the transportation policy, some sections among certain distribution paths can be seen as a break-off point, which is shown as the red line in Figure 1. Meanwhile, when vehicles arrive, consumers have to sign for the goods. Thus, the AVRP model can be identified as an automated vehicle routing problem with time windows (AVRPTW). Time windows are defined as ei,li. When the vehicle arrives before ei, then a vehicle has to wait, which is the cause of the penalty cost. When the vehicle arrives after li, this distribution service is regarded as a failure. The distribution centre has to plan a certain path for each automated vehicle and all the consumer demands must be satisfied: the objective is to minimise the total cost of the operation.Figure 1
AVRPTW model diagram.Considering practical circumstances, several assumptions are made: each automated vehicle is assumed to have the same capacity. Additionally, the total load in one journey cannot surpass this value. Each automated vehicle starts from the distribution centre and must return there after the delivery. Besides, each consumer is only served by one automated vehicle. Considering policy restrictions, some paths among certain consumers are break-off points.In addition, this model is defined as a VRP model with hybrid time windows. That is to say, when the vehicle arrives beforeei, the penalty cost and waiting time have a linear relationship. When the vehicle arrives after li, the penalty cost is seen as infinite. Considering the real operational process and the cost structure of a distribution centre, the total operating cost in this model is composed of three parts: fixed costs, transportation costs, and the penalty cost.The notation in this paper is summarised in Notations.


### 3.2. AVRPTW Model with Path Interruption

(1)
M
i
n
a
∑
k
=
1
K
∑
j
=
1
n
X
0
j
k
+
∑
k
=
1
K
∑
j
=
0
n
∑
i
=
0
n
c
i
j
X
i
j
k
+
∑
k
=
1
K
∑
i
=
1
n
α
max
⁡
0
,
s
i
-
T
i
k
+
M
max
⁡
0
,
T
i
k
-
l
i
,
(2)
s
.
t
.
∑
k
=
1
K
∑
i
=
0
n
X
i
j
k
=
1
j
=
1,2
,
…
,
n
,
(3)
∑
j
=
1
n
∑
i
=
1
n
q
i
X
i
j
k
≤
Q
k
=
1,2
,
…
,
K
,
(4)
∑
k
=
1
K
∑
i
=
0
n
X
i
0
k
≤
K
k
=
1,2
,
…
,
K
,
(5)
∑
i
=
0
n
X
i
0
k
≤
1
k
=
1,2
,
…
,
K
,
(6)
∑
i
=
0
n
X
i
j
k
=
∑
i
=
0
n
X
j
i
k
j
=
1,2
,
…
,
n
k
=
1,2
,
…
,
K
,
(7)
T
j
k
≥
T
i
k
+
f
i
+
t
i
j
+
max
⁡
0
,
s
i
-
T
i
k
i
=
1,2
,
…
,
n
j
=
0,1
,
2
,
…
,
n
k
=
1,2
,
…
,
K
,
(8)
∑
i
∈
S
∑
j
∈
S
X
i
j
k
≤
S
-
1
∀
S
⊆
I
k
=
1,2
,
…
,
K
,
(9)
X
i
j
k
=
0,1
,
(10)
X
g
h
k
=
0
g
,
h
∈
I
k
=
1,2
,
…
,
K
,
(11)
T
0
k
=
0
k
=
1,2
,
…
,
K
,
(12)
T
i
k
≥
0
i
=
1,2
,
…
,
n
k
=
1,2
,
…
,
K
.The objective function in (1) aims to minimise the fixed costs, transport costs, and penalty cost. Equation (2) implies that each consumer is only served by one automated vehicle. Equation (3) denotes the capacity limitation of each automated vehicle. Equation (4) requires that the total number of automated vehicles cannot surpass the sum of the vehicles that the distribution centre has at its disposal. Equation (5) implies that each path can only pass by the distribution centre once or not at all. Equation (6) shows that the number of vehicles entering a location equals the number of vehicles leaving it. Equation (7) defines the time windows. Equation (8) is the branch elimination constraint. Equation (9) defines the decision variables as binary values. Equation (10) denotes the interrupted paths. Equation (11) denotes the time when the automated vehicles leave the distribution centre. Equation (12) denotes that the time is nonnegative.The conventional VRP model with hybrid windows is known as a combinational optimisation problem, which is an NP-hard problem and cannot be solved by simple algorithms. Clearly, solving the AVRP problem is also difficult; thus to solve the AVRP model, existing exact algorithms cannot efficiently obtain the solution. Therefore, heuristic algorithms are efficient ways of solving this NP-hard problem. The particle swarm optimisation algorithm offers excellent performance with regard to solving the various VRP models. Therefore, the PSO algorithm is adopted here to solve the AVRPTW model.


## 3.1. Problem Statement

The automated vehicle routing problem (AVRP) model can be described as follows, and the model is shown in Figure1: the total number of the consumers in one certain area is N; the distribution centre has to distribute a certain amount of goods to certain consumers; presently, this distribution centre is preparing to put in several automated vehicles to carry out the assignment. In the initial stage of the application of the automated vehicle, considering the restrictions imposed by the transportation policy, some sections among certain distribution paths can be seen as a break-off point, which is shown as the red line in Figure 1. Meanwhile, when vehicles arrive, consumers have to sign for the goods. Thus, the AVRP model can be identified as an automated vehicle routing problem with time windows (AVRPTW). Time windows are defined as ei,li. When the vehicle arrives before ei, then a vehicle has to wait, which is the cause of the penalty cost. When the vehicle arrives after li, this distribution service is regarded as a failure. The distribution centre has to plan a certain path for each automated vehicle and all the consumer demands must be satisfied: the objective is to minimise the total cost of the operation.Figure 1
AVRPTW model diagram.Considering practical circumstances, several assumptions are made: each automated vehicle is assumed to have the same capacity. Additionally, the total load in one journey cannot surpass this value. Each automated vehicle starts from the distribution centre and must return there after the delivery. Besides, each consumer is only served by one automated vehicle. Considering policy restrictions, some paths among certain consumers are break-off points.In addition, this model is defined as a VRP model with hybrid time windows. That is to say, when the vehicle arrives beforeei, the penalty cost and waiting time have a linear relationship. When the vehicle arrives after li, the penalty cost is seen as infinite. Considering the real operational process and the cost structure of a distribution centre, the total operating cost in this model is composed of three parts: fixed costs, transportation costs, and the penalty cost.The notation in this paper is summarised in Notations.


## 3.2. AVRPTW Model with Path Interruption

(1)
M
i
n
a
∑
k
=
1
K
∑
j
=
1
n
X
0
j
k
+
∑
k
=
1
K
∑
j
=
0
n
∑
i
=
0
n
c
i
j
X
i
j
k
+
∑
k
=
1
K
∑
i
=
1
n
α
max
⁡
0
,
s
i
-
T
i
k
+
M
max
⁡
0
,
T
i
k
-
l
i
,
(2)
s
.
t
.
∑
k
=
1
K
∑
i
=
0
n
X
i
j
k
=
1
j
=
1,2
,
…
,
n
,
(3)
∑
j
=
1
n
∑
i
=
1
n
q
i
X
i
j
k
≤
Q
k
=
1,2
,
…
,
K
,
(4)
∑
k
=
1
K
∑
i
=
0
n
X
i
0
k
≤
K
k
=
1,2
,
…
,
K
,
(5)
∑
i
=
0
n
X
i
0
k
≤
1
k
=
1,2
,
…
,
K
,
(6)
∑
i
=
0
n
X
i
j
k
=
∑
i
=
0
n
X
j
i
k
j
=
1,2
,
…
,
n
k
=
1,2
,
…
,
K
,
(7)
T
j
k
≥
T
i
k
+
f
i
+
t
i
j
+
max
⁡
0
,
s
i
-
T
i
k
i
=
1,2
,
…
,
n
j
=
0,1
,
2
,
…
,
n
k
=
1,2
,
…
,
K
,
(8)
∑
i
∈
S
∑
j
∈
S
X
i
j
k
≤
S
-
1
∀
S
⊆
I
k
=
1,2
,
…
,
K
,
(9)
X
i
j
k
=
0,1
,
(10)
X
g
h
k
=
0
g
,
h
∈
I
k
=
1,2
,
…
,
K
,
(11)
T
0
k
=
0
k
=
1,2
,
…
,
K
,
(12)
T
i
k
≥
0
i
=
1,2
,
…
,
n
k
=
1,2
,
…
,
K
.The objective function in (1) aims to minimise the fixed costs, transport costs, and penalty cost. Equation (2) implies that each consumer is only served by one automated vehicle. Equation (3) denotes the capacity limitation of each automated vehicle. Equation (4) requires that the total number of automated vehicles cannot surpass the sum of the vehicles that the distribution centre has at its disposal. Equation (5) implies that each path can only pass by the distribution centre once or not at all. Equation (6) shows that the number of vehicles entering a location equals the number of vehicles leaving it. Equation (7) defines the time windows. Equation (8) is the branch elimination constraint. Equation (9) defines the decision variables as binary values. Equation (10) denotes the interrupted paths. Equation (11) denotes the time when the automated vehicles leave the distribution centre. Equation (12) denotes that the time is nonnegative.The conventional VRP model with hybrid windows is known as a combinational optimisation problem, which is an NP-hard problem and cannot be solved by simple algorithms. Clearly, solving the AVRP problem is also difficult; thus to solve the AVRP model, existing exact algorithms cannot efficiently obtain the solution. Therefore, heuristic algorithms are efficient ways of solving this NP-hard problem. The particle swarm optimisation algorithm offers excellent performance with regard to solving the various VRP models. Therefore, the PSO algorithm is adopted here to solve the AVRPTW model.


## 4. Solution Algorithm

### 4.1. The Particle Swarm Optimisation (PSO) Algorithm

The PSO algorithm is an evolutionary algorithm, which originates from research into bird flock preying behaviour. The PSO algorithm is similar to a genetic algorithm, both of which are based on iteration.In the PSO algorithm, a particle swarm is initialised, which is also identified as the initial solution. Then the particle follows the best particle through iteration. For theith particle in the dth dimension, the position vector and the velocity vector are defined as follows: Xi=xi,1xi,2⋯xi,d and Vi=vi,1vi,2⋯vi,d. During each iteration, each particle refreshes itself by following the two types of best solutions: (i) the optimal solution for each individual and (ii) the optimal solutions of the population. After finding these two optimal solutions, the particle refreshes its position and velocity according to(13)vi,jt+1=ωvi,jt+c1r1pi,j-xi,jt+c2r2pg,j-xi,jt,xi,jt+1=xi,jt+vi,jt+1j=1,2,…,d.The specific steps are as follows.Step 1.
Random initialisation of the position and velocity of all particles.Step 2.
Calculation of the adaptive value for each particle and then storing the optimal solution of the present individual inpbest, while storing the optimal solution of the population ingbest.Step 3.
Refreshment of the position vector and velocity vector.Step 4.
Adjustment of the optimal solution of each individual.Step 5.
Adjustment of the optimal solution of the population.Step 6.
Execution of a stop-check as to whether, or not, the number of iterations is greater than the preset end value.Step 7.
Output and end.


### 4.2. The Improvement of the PSO Algorithm

The adjustable parameters in the PSO algorithm include learning factorsc1, c2 and an inertia weight w. Inertia weight is used to control the impact of histories of velocities on the current velocity. The particle adjusts its trajectory based on information about its previous best performance and the best performance of its neighbours. Learning factors, c1 and c2, decide how the interaction between the particle itself and its neighbour affects the adjustment of both position and velocity. Thus, they are used to control the convergence behaviour of the PSO. And the function y=cos⁡x is applied to adjust the inertia weight and learning factors.(1
) The Improvement of Inertia Weight. Apply y=cos⁡x in the formula for the inertia weight, and the improved expression for inertia weight is given by(14)W=wmax-wmax-wmincos⁡tT-1π,where wmax and wmin denote the maximum and the minimum inertia weight, respectively: these two values are defined as 0.9 and 0.4 according to the last experienced value. The range of t/T is 0,1.(2
) The Improvement of Learning Factors. According to the improvement of inertia weight, the calculation of learning factors is given by(15)c1=c1max-c1max-c1mincos⁡tT-1π,c1=c2.


### 4.3. The PSO Algorithm for AVPRTW Model

Applying the PSO algorithm to solve the AVRP model, it is vital to realise that each particle should be in correspondence with its coding. Two spatial vectors are designed to represent the solution:k denotes the number of the vehicle for each path, while r denotes the precedence order of vehicle arrival at the consumers. These two spatial vectors are Xk and Xr. To illustrate the coding, an example is given; assuming that there are 10 consumers in this area and the distribution centre is equipped with three vehicles, we assume an optimal solution as follows:(16)Xk1332233111Xr2321241413The optimal solution denotes that all of the three vehicles are assigned to execute the distribution task. The first vehicle serves consumers 1, 8, 9, and 10. The second vehicle serves consumers 4 and 5, and consumers 2, 3, 6, and 7 are served by the third vehicle. According toXr, the distribution order of the first vehicle is 0→9→1→10→8→0, that of the second vehicle is 0→4→5→0, and that of the third vehicle is 0→7→3→2→6→0.The specific steps are as follows.Step 1 (initialisation of the particle swam).
Firstly, several parameters of the PSO algorithm have to be set, including the population, the maximum number of iterations, and the spatial dimension.Xk and Xr are, respectively, confined to [1,K] and [1,N]. The corresponding velocity vectors Vk and Vr are, respectively, confined to [-(K-1),(K-1)] and [-(N-1),(N-1)].Step 2.
According to (14) and (15), the position and velocity are refreshed. The refreshed values of position and velocity are substituted into the fitness function and the optimal value of each individual and the population are initialised.Step 3.
Compare the adaptive value with its optimal value, and refreshpbest.Step 4.
Calculate the refreshed optimal value of the population, compare this value to the presentgbest, and refresh the value ofgbest.Step 5.
Judge whether, or not, the maximum number of iterations has been reached, and if so, output the optimal solution (otherwise, return to Step2).


## 4.1. The Particle Swarm Optimisation (PSO) Algorithm

The PSO algorithm is an evolutionary algorithm, which originates from research into bird flock preying behaviour. The PSO algorithm is similar to a genetic algorithm, both of which are based on iteration.In the PSO algorithm, a particle swarm is initialised, which is also identified as the initial solution. Then the particle follows the best particle through iteration. For theith particle in the dth dimension, the position vector and the velocity vector are defined as follows: Xi=xi,1xi,2⋯xi,d and Vi=vi,1vi,2⋯vi,d. During each iteration, each particle refreshes itself by following the two types of best solutions: (i) the optimal solution for each individual and (ii) the optimal solutions of the population. After finding these two optimal solutions, the particle refreshes its position and velocity according to(13)vi,jt+1=ωvi,jt+c1r1pi,j-xi,jt+c2r2pg,j-xi,jt,xi,jt+1=xi,jt+vi,jt+1j=1,2,…,d.The specific steps are as follows.Step 1.
Random initialisation of the position and velocity of all particles.Step 2.
Calculation of the adaptive value for each particle and then storing the optimal solution of the present individual inpbest, while storing the optimal solution of the population ingbest.Step 3.
Refreshment of the position vector and velocity vector.Step 4.
Adjustment of the optimal solution of each individual.Step 5.
Adjustment of the optimal solution of the population.Step 6.
Execution of a stop-check as to whether, or not, the number of iterations is greater than the preset end value.Step 7.
Output and end.


## 4.2. The Improvement of the PSO Algorithm

The adjustable parameters in the PSO algorithm include learning factorsc1, c2 and an inertia weight w. Inertia weight is used to control the impact of histories of velocities on the current velocity. The particle adjusts its trajectory based on information about its previous best performance and the best performance of its neighbours. Learning factors, c1 and c2, decide how the interaction between the particle itself and its neighbour affects the adjustment of both position and velocity. Thus, they are used to control the convergence behaviour of the PSO. And the function y=cos⁡x is applied to adjust the inertia weight and learning factors.(1
) The Improvement of Inertia Weight. Apply y=cos⁡x in the formula for the inertia weight, and the improved expression for inertia weight is given by(14)W=wmax-wmax-wmincos⁡tT-1π,where wmax and wmin denote the maximum and the minimum inertia weight, respectively: these two values are defined as 0.9 and 0.4 according to the last experienced value. The range of t/T is 0,1.(2
) The Improvement of Learning Factors. According to the improvement of inertia weight, the calculation of learning factors is given by(15)c1=c1max-c1max-c1mincos⁡tT-1π,c1=c2.


## 4.3. The PSO Algorithm for AVPRTW Model

Applying the PSO algorithm to solve the AVRP model, it is vital to realise that each particle should be in correspondence with its coding. Two spatial vectors are designed to represent the solution:k denotes the number of the vehicle for each path, while r denotes the precedence order of vehicle arrival at the consumers. These two spatial vectors are Xk and Xr. To illustrate the coding, an example is given; assuming that there are 10 consumers in this area and the distribution centre is equipped with three vehicles, we assume an optimal solution as follows:(16)Xk1332233111Xr2321241413The optimal solution denotes that all of the three vehicles are assigned to execute the distribution task. The first vehicle serves consumers 1, 8, 9, and 10. The second vehicle serves consumers 4 and 5, and consumers 2, 3, 6, and 7 are served by the third vehicle. According toXr, the distribution order of the first vehicle is 0→9→1→10→8→0, that of the second vehicle is 0→4→5→0, and that of the third vehicle is 0→7→3→2→6→0.The specific steps are as follows.Step 1 (initialisation of the particle swam).
Firstly, several parameters of the PSO algorithm have to be set, including the population, the maximum number of iterations, and the spatial dimension.Xk and Xr are, respectively, confined to [1,K] and [1,N]. The corresponding velocity vectors Vk and Vr are, respectively, confined to [-(K-1),(K-1)] and [-(N-1),(N-1)].Step 2.
According to (14) and (15), the position and velocity are refreshed. The refreshed values of position and velocity are substituted into the fitness function and the optimal value of each individual and the population are initialised.Step 3.
Compare the adaptive value with its optimal value, and refreshpbest.Step 4.
Calculate the refreshed optimal value of the population, compare this value to the presentgbest, and refresh the value ofgbest.Step 5.
Judge whether, or not, the maximum number of iterations has been reached, and if so, output the optimal solution (otherwise, return to Step2).


## 5. Case Study

### 5.1. Background Information

One certain chain supermarket has one distribution centre (DC) in this area and there are 25 stores. The distribution centre adopts a self-distribution strategy and has to schedule deliveries and plan the path of each vehicle, so that all demands from the 25 stores are satisfied. The geographical position of the distribution and all stores are shown in Figure2 and the position coordinates are simplified and specific data are given in Table 1.Table 1
Position coordinates and time windows.
Number
Coordinate
Earliest
Latest
Demand
Service time
DC
(
35,35
)
—
—
—
—
1
(
41,49
)
161
171
10
10
2
(
35,17
)
50
60
7
10
3
(
55,45
)
116
126
13
10
4
(
55,20
)
149
159
19
10
5
(
15,30
)
34
44
26
10
6
(
25,30
)
99
109
3
10
7
(
20,50
)
81
91
5
10
8
(
10,43
)
95
105
9
10
9
(
55,60
)
97
107
16
10
10
(
30,60
)
124
134
16
10
11
(
20,65
)
67
77
12
10
12
(
50,35
)
63
73
19
10
13
(
30,25
)
159
169
23
10
14
(
15,10
)
32
42
20
10
15
(
30,5
)
61
71
8
10
16
(
10,20
)
75
85
19
10
17
(
5,30
)
157
167
2
10
18
(
20,40
)
87
97
12
10
19
(
15,60
)
76
86
17
10
20
(
45,65
)
126
136
9
10
21
(
45,20
)
62
72
11
10
22
(
45,10
)
97
107
18
10
23
(
55,5
)
68
78
29
10
24
(
65,35
)
153
163
3
10
25
(
65,20
)
172
182
6
10Figure 2
Locations of the distribution centre and consumers.In the practical delivery operation, considering the employee cost and risk cost caused by fatigue while driving, the enterprise is preparing to assign several automated vehicles in these real services. Due to anxiety from the public, prohibited sections are used in the initial stage and break-off sections are marked by the red dotted line in Figure2. Each store should arrange some people to be on duty within a specific time frame, the time windows, and the service time of each store is given in Table 1.


### 5.2. Analysis of the Case Study

Based on the AVRPTW model, the improved PSO algorithm is used to solve the problem. To simplify the calculation, some parameters of the model are assumed when carrying out the experimental test. The specific parameters of the model are as listed in Table2. With the help of MATLAB™ 2015b, this problem is solved and the solutions are given in Figure 3 and Table 3.Table 2
Model parameters.
Parameter
Value
Fixed cost of each automated vehicle
200
Index of the penalty function
1.5
Capacity
100
Transport cost per unit hauled
8
Average speed of an automated vehicle
40Table 3
Paths and load factors.
Number
Path
Load factor (%)
1
0
→
24
→
17
→
23
→
18
→
21
→
4
→
0
76
2
0
→
11
→
1
→
7
→
5
→
13
→
0
76
3
0
→
20
→
9
→
10
→
3
→
12
→
2
→
0
80
4
0
→
16
→
14
→
15
→
19
→
8
→
6
→
22
→
25
→
0
100Figure 3
Convergence of the solution algorithm.


### 5.3. Analysis and Discussion

According to the solution, the distribution centre has to allocate four automated vehicles to execute the task. All of the stores are included in these planned paths. Based on Figure3, after the 145th iteration, the algorithm converges and the optimal solution is given. The minimum total cost is 11,328. From Figure 3, it is seen that the convergence performance of the PSO algorithm is good, and its antijamming capability is realised.Considering four paths are interrupted, a total four automated vehicles are assigned to serve all the stores, and the load factor of each one is 76%, 76%, 80%, and 100%, respectively. The capacities of the automated vehicles are almost fully utilised. The results show that the transportation resources available have been used to full advantage.


## 5.1. Background Information

One certain chain supermarket has one distribution centre (DC) in this area and there are 25 stores. The distribution centre adopts a self-distribution strategy and has to schedule deliveries and plan the path of each vehicle, so that all demands from the 25 stores are satisfied. The geographical position of the distribution and all stores are shown in Figure2 and the position coordinates are simplified and specific data are given in Table 1.Table 1
Position coordinates and time windows.
Number
Coordinate
Earliest
Latest
Demand
Service time
DC
(
35,35
)
—
—
—
—
1
(
41,49
)
161
171
10
10
2
(
35,17
)
50
60
7
10
3
(
55,45
)
116
126
13
10
4
(
55,20
)
149
159
19
10
5
(
15,30
)
34
44
26
10
6
(
25,30
)
99
109
3
10
7
(
20,50
)
81
91
5
10
8
(
10,43
)
95
105
9
10
9
(
55,60
)
97
107
16
10
10
(
30,60
)
124
134
16
10
11
(
20,65
)
67
77
12
10
12
(
50,35
)
63
73
19
10
13
(
30,25
)
159
169
23
10
14
(
15,10
)
32
42
20
10
15
(
30,5
)
61
71
8
10
16
(
10,20
)
75
85
19
10
17
(
5,30
)
157
167
2
10
18
(
20,40
)
87
97
12
10
19
(
15,60
)
76
86
17
10
20
(
45,65
)
126
136
9
10
21
(
45,20
)
62
72
11
10
22
(
45,10
)
97
107
18
10
23
(
55,5
)
68
78
29
10
24
(
65,35
)
153
163
3
10
25
(
65,20
)
172
182
6
10Figure 2
Locations of the distribution centre and consumers.In the practical delivery operation, considering the employee cost and risk cost caused by fatigue while driving, the enterprise is preparing to assign several automated vehicles in these real services. Due to anxiety from the public, prohibited sections are used in the initial stage and break-off sections are marked by the red dotted line in Figure2. Each store should arrange some people to be on duty within a specific time frame, the time windows, and the service time of each store is given in Table 1.


## 5.2. Analysis of the Case Study

Based on the AVRPTW model, the improved PSO algorithm is used to solve the problem. To simplify the calculation, some parameters of the model are assumed when carrying out the experimental test. The specific parameters of the model are as listed in Table2. With the help of MATLAB™ 2015b, this problem is solved and the solutions are given in Figure 3 and Table 3.Table 2
Model parameters.
Parameter
Value
Fixed cost of each automated vehicle
200
Index of the penalty function
1.5
Capacity
100
Transport cost per unit hauled
8
Average speed of an automated vehicle
40Table 3
Paths and load factors.
Number
Path
Load factor (%)
1
0
→
24
→
17
→
23
→
18
→
21
→
4
→
0
76
2
0
→
11
→
1
→
7
→
5
→
13
→
0
76
3
0
→
20
→
9
→
10
→
3
→
12
→
2
→
0
80
4
0
→
16
→
14
→
15
→
19
→
8
→
6
→
22
→
25
→
0
100Figure 3
Convergence of the solution algorithm.


## 5.3. Analysis and Discussion

According to the solution, the distribution centre has to allocate four automated vehicles to execute the task. All of the stores are included in these planned paths. Based on Figure3, after the 145th iteration, the algorithm converges and the optimal solution is given. The minimum total cost is 11,328. From Figure 3, it is seen that the convergence performance of the PSO algorithm is good, and its antijamming capability is realised.Considering four paths are interrupted, a total four automated vehicles are assigned to serve all the stores, and the load factor of each one is 76%, 76%, 80%, and 100%, respectively. The capacities of the automated vehicles are almost fully utilised. The results show that the transportation resources available have been used to full advantage.


## 6. Conclusion

Considering the partial sections interruption and the time windows in the initial stage of the application of automated vehicles in last-mile mode, the conventional VRP model is modified. The AVRPTW model with traffic interruptions is established. Additionally, an improved PSO algorithm is designed by changing the inertia weight and learning factors. Finally, a case study has been analysed and the validity and practicability of the mathematical model and algorithm have been verified.Due to the immaturity of unmanned vehicle technology, the promotion and popularisation of automated vehicles in real applications are likely to be faced with strict restrictions in the initial stage. First, it is inevitable that automated vehicles must exert a negative influence on the public and general transport environment. Second, the potential risks caused by traffic congestion and the unregulated behaviour of other drivers are likely to affect the promotion of such unmanned technology and unmanned vehicles in general. Accordingly, a policy that partially opens, or restricts, road sections for automated vehicles can both ensure the safety of the automated vehicles and reduce anxiety about unmanned vehicle technology from the public. Thus, it is of practical significance to study automated vehicle route problems.Apart from the aforementioned problems, there may also be other issues that must be taken into consideration when applying unmanned vehicles in the logistics sector, such as the influence that automated vehicles have on ordinary vehicles and the safety of automated vehicles. These issues should be addressed in future research.



---

*Source: 1624328-2017-08-07.xml*