Research Article  Open Access
Yanling Feng, RenQian Zhang, Guozhu Jia, "Vehicle Routing Problems with Fuel Consumption and Stochastic Travel Speeds", Mathematical Problems in Engineering, vol. 2017, Article ID 6329203, 16 pages, 2017. https://doi.org/10.1155/2017/6329203
Vehicle Routing Problems with Fuel Consumption and Stochastic Travel Speeds
Abstract
Conventional vehicle routing problems (VRP) always assume that the vehicle travel speed is fixed or timedependent on arcs. However, due to the uncertainty of weather, traffic conditions, and other random factors, it is not appropriate to set travel speeds to fixed constants in advance. Consequently, we propose a mathematic model for calculating expected fuel consumption and fixed vehicle cost where average speed is assumed to obey normal distribution on each arc which is more realistic than the existing model. For smallscaled problems, we make a linear transformation and solve them by existing solver CPLEX, while, for largescaled problems, an improved simulated annealing (ISA) algorithm is constructed. Finally, instances from real road networks of England are performed with the ISA algorithm. Computational results show that our ISA algorithm performs well in a reasonable amount of time. We also find that when taking stochastic speeds into consideration, the fuel consumption is always larger than that with fixed speed model.
1. Introduction
The vehicle routing problem (VRP) is one of the most important and studied combinatorial optimization problems [1]. It has gained great attentions from many researchers, especially in distribution and logistics management fields. More than fifty years have elapsed since Dantzig and Ramser [2] firstly introduced the problem in in 1959. They described a realworld application concerning the delivery of gasoline to service stations and proposed the first mathematical programming formulation and algorithmic approach. Following this seminal paper, hundreds of models and algorithms were proposed for the optimal and approximate solution of the different versions of VRP. Then in order to be better aligned with the realworld applications, many different VRP versions have been studied. The most common version is the CVRP [3], where each vehicle has limited load capacity. The VRP with Time Window (VRPTW) [4–6] aims to find optimal route sets with minimum total travel cost while serving each customer within specified time window. Other extended versions include VRP with backhauls [7], VRP with pickups and deliveries [8], and the multidepot VRP [9].
The deterministic VRP cannot cover all the situations in reality while considering stochastic VRP components. Consequently, Stochastic VRP are developed. For example, Ritzinger et al. [10] made a detailed review of the stochastic VRP. Mehrjerdi [11] combined chance constrained programming and multiple objective programming to obtain satisfactory solutions. Taş et al. [12, 13], Ehmke et al. [14], and Laporte et al. [15] put forward different heuristic methods for VPR with stochastic travel time and time windows. Marinakis et al. [16] developed a particle swarm algorithm for VRP with stochastic demands. For more stochastic demand results, one may refer to [17–19].
Based on the NP hardness of VRP, multiple heuristic algorithms have been put forward to solve this problem. Xiao and Konak [20] present simulating annealing algorithm to solve the green vehicle routing and scheduling problem with hierarchical objectives and weighted tardiness. Kondekar et al. [21] provide a mapreduce based hybrid genetic solution for solving largescale vehicle routing problems in dynamic network with fluctuant link travel time. Neural network is also applied to solve stochastic multiconstraint problems with different timescales; see Zhang et al. [22] and MeyerBäse et al. [23].
Demir et al. [24] proposed an adaptive large neighborhood search algorithm (ALNS) to minimize the fuel consumption and the driving time with Pareto optimality. Garaix et al. [25] proposed a column generation algorithm for a dialaride problem.
Due to the uncertain traffic factors, such as unexpected workload or bad weather, stochastic VRP has attracted more and more attentions. Cao et al. [26] proposed a partial Lagrange multiplier method. Ishigaki [27] considered a dynamic collection plan with stochastic demand and apply their search algorithm to an actual trash collection problem. Novel applications of the wellestablished problem can also be found in optical flow and vehicle systems [28–30].
In most literatures, travel speeds are assumed to be fixed or timedependent (see, e.g., [31, 32]). However, in practice, due to the uncertainty of weather, traffic conditions, and other random factors, it is not appropriate to set travel speeds to be fixed constants in advance. The interest in stochastic VRP in this paper is motivated by both its practical relevance and its considerable difficulty: large VRP instances may be solved to optimality only in particular cases. Therefore, we study vehicle routing problems with stochastic travel speeds. Moreover, as each arc has limited speed and other random factors, the average speed of the same type of vehicles on the same arc approximately obeys normal distribution.
Figure 1 is an illustration of this VPR with stochastic average travel speed. Assume that a logistics company owns a fleet of trucks and one truck delivers goods for customers 1, 2, and 3. On the first day, the truck travels with an average speed of 25 m/s on arc. On the second day, it rains heavily when it delivers goods for customer 1. The truck has to get over the poor road conditions caused by the heavy rain; as a result, the average speed it travels on arc becomes 15 m/s. Several days later, a traffic accident occurred on arc, so it may travel with an average speed of 10 m/s. Practically, during a relatively long period of time, the average speed on arc is not fixed but with slight fluctuation. From this point of view, the average speed is a stochastic variable. Here we assume that the average speed on each arc follows normal distribution.
On the other hand, environmental issues have become worldwide problems. In fact, fuel consumption of traffic vehicles is a great contributor to CO_{2 }emissions which is the main culprit of global warming [33, 34]. As a result, we take the sum of the expected fuel consumption and total vehicle cost as the objective to evaluate different routes. It has been studied that the fuel consumed of a vehicle traveling along a route depends on many factors, which include distance, load, speed, road conditions, and vehicle types. Models can be found in Bektaş and Laporte [34] and Xiao et al. [35].
Table 1 summarizes and compares the formulation of VRP in the stream of uncertain VRP and green VRP.

Two main contributions of this paper are described as follows:(1)Based on the fuel consumption model introduced by Bektaş and Laporte [34], the fuel consumption is a nonlinear function of travel speed, distance, and vehicle load. Due to the stochasticity of the travel speeds, we extend the model with stochastic speeds considered.(2)An ISA (simulated annealing algorithm) is presented to solve largescaled VRP problems.
The remainder of this paper is as follows. In Section 2, the description and formulation of this model is provided. In Section 3, we made a linearization of this problem. Then Section 4 introduces the improved simulated annealing algorithm with memory for optimizing the routing plan. In Section 5, the computational experiments are performed. Finally, Section 6 presents the conclusions and managerial insights.
2. Problem Statement and Model Formulation
Generally speaking, logistics companies tend to make budget decisions for a planning period; for example, they figure out the amount of the fuel consumption, the cost of keeping, and maintaining a fleet of vehicles in advance. Thus, we try to develop a fuel consumption model and algorithm to obtain more precise budgets. We use a digraph to describe the vehicle routing problems with fuel consumption and stochastic speeds (VRPFSV). Let a complete connected digraph be the logistics network with a node set and an arc set . Arc represents the path from node to node . The depot is denoted by node 0 and the other nodes in represent customers with nonnegative demand . The depot owns enough homogeneous vehicles with limited capacity , so the total demands of customers assigned to the same vehicle must be less than or equal to . Each route is finished by only one vehicle and each customer is served only once. The average travel speed on arc is stochastic. After serving all the assigned customers, each vehicle has to return to the depot.
2.1. Assumptions and Notations
In this section, the assumptions and notations used in the model formulation are listed. First, we assume the following:(1)Each demand must be satisfied and each customer is served only once.(2)The depot owns enough homogenous vehicles.(3)Each vehicle must departure from the depot and after having served its customers it must return to the depot.(4)Time window constraints are not considered here.
Then the notations used in the VRPFSV formulation are listed as follows: : total number of customers. : index of nodes (the depot is represented by 0). : fixed cost of using a vehicle. : the capacity of a vehicle. : set of nodes including the depot, . : set of arcs formed by all pairs of nodes, . : distance between node and node . : stochastic average travel speed on arc.
Two decision variables are as follows: : binary variable indicating whether arc is traveled. If arc is traveled by a vehicle, then ; otherwise, . : load carried from node to node .
2.2. Formulation of the Model of VRPFSV
In this section, a mixed integer linear programming model is developed to construct the VRPFSV. The VRPFSV is an extension of the classical VRP and assumes that average travel speed on arc is a nonnegative random variable. Furthermore, the average travel speed on arc is assumed to obey normal distribution with mean and variance .
The fuel consumption model applied in this paper is based on the comprehensive model formed by Barth and Boriboonsomsin [33]. It is given by where and are constants decided by different fuel properties, , and are vehicle characteristics related constants. is a constant associated with the road characteristics and acceleration, , where is acceleration and is the road angle inclination. The reference values of all parameters are given in Table 2. Assuming that acceleration and road inclination are zeros, then the fuel consumption on arc can be rewritten as where , and the price of the fuel . Equation (2) establishes a good correlation between fuel consumption and travel speed, weight, and distance.
 
Source: Demir et al., 2012 [24]. 
In our research, the service cost includes two components: the expected fuel consumption and the fixed vehicle cost.
Consider arc; if it is traveled by a vehicle with average speed of , then it has an expected fuel consumption; according to (2), the expected fuel consumption is given in
Consequently, the mathematical formulation of VRPFSV can be expressed as follows:
In this optimization model, the two variable sets are and . Equation (4) is the objective function which consists of the expected fuel consumption and the total vehicle cost. indicates a departure of a newly used vehicle from the depot; as any used vehicle will trigger a fixed cost, the sum of represents the total vehicle cost. Constraint (5) represents that the vehicle visits each node (except the depot) only once. Equation (6) represents the conservation of flow constraints. Equation (7) ensures that the demand of each customer must be met and it also eliminates subtours. Equation (8) guarantees that the vehicle load cannot exceed the vehicle capacity . Equation (9) indicates that is a binary variable. Equation (10) is the nonnegative constraint of vehicle load.
3. Linearization and Solution by CEPLEX
Since the VRP is NPhard and this VRPFSV model includes stochastic travel speed, it is at least NPhard. In this section, we transform the original model into a linear one, and then small instances can be solved easily by CEPLEX 12.6.2.
This VRPFSV model is nonlinear due to the existence item in the objective. In fact, as the Hessian matrix of the objective is nonpositive, the VRPFSV model is even not convex. Fortunately, we manage to make the linearization and obtain the linear model which is equivalent to the original one.
For any given arc, consider its expected fuel consumption:
In fact, the total fuel consumption can be reduced as follows:
Property 1. In the VRPFSV, (11) is equivalent to (12), subjecting to constraint (8).
Proof of Property 1. In order to prove this property, we only need to prove that the product of and is equivalent to . In fact, is a binary variable; if , then , else if , due to constraints (8) and (10), we get . Consequently, . Thus, is equivalent to which indicates that (11) is equivalent to (12).
According to (12), objective (4) is converted to the equivalent linear form:LetAs a result, objective (4) can be transformed into , which is a linear objective function about and now. Moreover, constraints (5)–(10) are linear about the decision variables, so we established a linear mixed integer programming model for VRPFSV.
Then the following intractable problem is how to calculate and , as it is unable to get the closed form of the integral in (14). In the following section, we use numerical integration to obtain and store the value of and on every arc in advance.
Now all the expressions in the model are linear, so we got the coefficient value in the objective. It is possible to optimally solve problem VRPFSV by the existing solver CPLEX (version 12.6.2) for small sized problem. The linearization plays a significant role in finding optimal solution and reducing computation time.
4. ISA Algorithm for LargeScaled VRPFSV
For largescaled instances, the existing solver will go expired, as the computation process will cost a very large amount of time. We develop an algorithm based on SA algorithm. This ISA algorithm includes four parts, construction of the initial solution, generating neighborhood solutions, local search, and replacement of the current best solution. And theoretically, the SA algorithm is guaranteed to converge to the global optimization solution with probability one [43].
4.1. Construction of Initial Solution
In the string model, the unique depot is encoded as 0 and the customers are encoded as a series of positive numbers. The string vector which represents the initial route is formatted as and the length of is .
Then make feasible for VRPFSV; that is, assign customers to appropriate number of vehicles. Usually, vehicle load contributes a lot to the fuel consumption. For example, if a vehicle carries a heavy load starting from the depot, it manages to serve more customers. On the contrary, if the initial load is small, more vehicles will be needed because of the capacity limit. According to vehicle capacity, string is transformed to feasible solution with the following procedures. First, start from the second character of , and then accumulate the demands of nodes successively. If the cumulative demands exceed the vehicle capacity, start a new vehicle; that is, insert a 0 into string before the current character and then reset the cumulative demands to 0. Repeat the above steps until the end of . Take as an example, assuming that is the initial string. As we do not know how many zeros are needed to be inserted into initially, we set the length of VRPX to the largest likely length 21. Because the cumulative demands of 3, 4, and 1 are less than the capacity and the cumulative demands of 3, 4, 1, and 7 exceed the capacity, 0 is inserted between 1 and 7. Continue checking until the end of ; then we can obtain , which represents ten customers being visited by four vehicles.
4.2. Generate Neighborhood Solutions
Exchange rule: three commonly used exchange rules swap, relocation, and 2opt (see Figure 2) are used to configure the neighborhood solutions. After exchanging, we get a neighborhood solution, and then make it feasible following steps of Section 4.1. Finally, the neighborhood solutions are obtained.
(a) Swap
(b) Relocation
(c) 2opt
4.3. Local Search and Update the Current Best Solution
The objective function value is the criterion to evaluate which solution is better. If the total cost becomes smaller, the newly generated solution is accepted. Otherwise, the solution is accepted with certain probability. As poor solutions maybe accepted at certain probability, so the best of the output cannot be guaranteed. Consequently, we use a memory array in the algorithm to record the best solution. Only when a better solution appears, will be updated. After continuous improvement, we can get the best solution in all searched neighborhoods.
Furthermore, it is not so easy to calculate the vehicle load on each arc, so once a route is generated, the load will be calculated in a reverse order. More details are listed in in Appendix B.
This algorithm adopts adaptive cooling processing, and the temperature cooling coefficient can be set to , where is the accepted solution numbers of the current temperature and loop is the total loop numbers. In this way, if the number of accepted solutions is small, that is, the current solution is close to optimal solution, the coefficient may lead to smaller search scope which may cost less time. On the other hand, if the number of accepted solutions is big, a relatively big cooling coefficient will help expand the search scope.
The main steps of this improved simulated annealing algorithm are as follows.
Step 1. Initialize parameters. Input customers with demand and distance matrix , vehicle capacity , fixed single vehicle cost , speed distribution parameters of each arc and , the end temperature , and the number of inner loops loop.
Step 2. Generate initial string according to the method introduced in Section 4.1.
Step 3. Generate initial temperature. Randomly change to its neighbor for 1000 times; the maximum objective deviation is selected as the initial temperature.
Step 4. Whether the inner loop number is reached, if it is, go to Step 8. Otherwise, generate new solutions and transform it into feasible solution.
Step 5. Calculate the deviation between the current cost and the former one and accept or reject it according to Metropolis rule.
Step 6. Update the best so far solution.
Step 7. Whether the improvement is less than 0.01, if it is, go to Step 9.
Step 8. Whether the end temperature is reached, if so, go to Step 9. Otherwise, drop the temperature and go to Step 4.
Step 9. Stop and output the best solution. The pseudocode of the improved annealing algorithm is listed in Appendix A.
In Step 3, we use selfadaptive method to obtain the best initial temperature. As too high initial temperature will increase the computational time and too low initial temperature will trap in local optimal solution. Here, we generate the maximum temperature adaptively according to the problem; that is, randomly change the incumbent solution for 1000 times. The maximum cost deviation between two neighboring solution is selected as the initial temperature.
Moreover, the temperature drops adaptively, under certain temperature; if the number of solutions being accepted is large, then the temperature drops very quickly. Else, if the number of accepted solutions is small, the temperature drops relatively slowly.
In Step 5, let and be the objective value of the th iteration and the best objective of the accepted routes, and . and denote, respectively, the minimum cost route among the accepted ones and the th route. If , then replace with . Else, replace with if .
The algorithm terminates on two conditions: the current temperature is below or equal to the end temperature , or the improvement of the best solution is less than 0.01.
Theory 2. In the VRPFSV, the time complexity of the proposed SA algorithm is .
Proof. In the outer temperature loop, the number of temperatures in the cooling procedure is ; in the internal loop, the number of loops is loop, and each loop has a calculation of cost function with time complexity . Thus, the total time complexity becomes . As , , , and are all given constants, Theory 2 is proved.
5. Computational Experiments of the ISA Algorithm
In this section, the results of computational applications of ISA algorithm are presented.
First, we introduce the PollutionRouting Problem Instance Library (PRPLIB) (http://www.apollo.management.soton.ac.uk/prplib.htm). This library contains nine differentscale groups of instances of the PollutionRouting Problem (PRP). Each group consists of 20 different instances. These instances are based on real distances collected from UK cities. The first number in the file name after UK shows the number of nodes contained in the instance. The second is the order number of the instance within the group. The problem consists of nodes, and the format of each file includes data about number of customers, vehicle capacity, city name and demand, distance matrix, and minimum and maximum speed level.
As there are no suitable benchmark problems of this model, we modified the benchmark problems of the PRPLIB to test our algorithm. Suppose that is the average travel speed on arc and it obeys normal distribution with mean parameter and variance parameter .
In these examples, speed limit of all arcs is 5 m/s~25 m/s. For each arc, the integer mean parameters are generated randomly from interval [5, 25] and set . At the same time, the other parameters stay unchanged. The ISA algorithm was implemented in MATLAB and executed on an Intel 2.0 GHz processor with 1.59 G of RAM.
5.1. Quality and Efficiency of the ISA Algorithm
Tables 2 and 3 compared results of the optimal solutions found by CPLEX, SA, and ISA algorithm when is small. The SA and ISA algorithm are performed 10 times for each instance. The average and best objective value of the 10 runs are listed in the tables, respectively. We label the metrics associated with each objective as follows: the first column is the instance ID, and the second and third columns are the optimal objective value and computation time of solutions solved by the solver CPLEX. AvgObj is the average total cost found by our SA or ISA algorithm in 10 runs. BestObj is the best objective of the 10 runs. Dev. is the relative deviations between the best objective obtained in 10 runs and the optimal cost solved by CPLEX or SA; that is, . AvgCPU time is the average run time (CPU) of the SA or ISA algorithm.

5.1.1. Parameter Settings
As different parameter settings may influence the performance of ISA heuristic, we test different parameter values in advance. Furthermore, the initial temperature and temperature drop factor are generated adaptively; we only need to determine the appropriate values of and . We change one parameter at a time to observe how that parameter affects the solution. It can be noticed from Figure 3 that using a larger or a smaller may improve the solution quality a little but increase the computation time significantly.
By trading off the effectiveness and the efficiency of the algorithm, we found the values of and are appropriate parameters combination.
As analyzed above, we set , , , and . , , , are set to the same value as mentioned in Section 2.2; that is, , , , and . In order to compare the results obtained by SA and ISA, we set the same parameter with the same value both in SA and in ISA. Thus, the parameters in SA is as follows: , , and .
In Table 3, the CPLEX is always superior to SA and ISA in both objective and computation time, and the largest relative deviation between ISA and CPLEX is 5.88%. Table 4 depicts in most cases that the ISA manage to obtain the optimal solution and the computation time turns out to be less than CPLEX. Bold numbers indicate that, in four instances, the ISA algorithm finds the optimal solution. In all cases, ISA consumes much less time than classic SA while obtaining the same solution.

It is shown in Table 5 that when , the largest deviation of the ISA and CPLEX is 1.46%, and the ISA obtains the near optimal solution in less time than both CPLEX and ISA. For three instances the ISA managed to get the optimal solution.

In Table 6, results from the ISA outperform the solver CPLEX. As can be seen, the biggest deviation between ISA and CPLEX is 5.25% and the longest running time of ISA is less than 60 seconds, while solver CPLEX costs more than 8 hours. On the other hand, ISA costs less time and obtain better solutions than SA. These results highlight the high quality and time advantage of the proposed SA algorithm for largescaled problems.

For largescaled problems, if we set the time limit of CPLEX solver to 12 hours, all the computing process terminated without any solution. Thus, Tables 7 and 8 only perform the ISA and SA results. As can be seen, the longest running time of ISA is less than 100 seconds, which illustrates that the proposed algorithm performs well especially for largescaled problems.


5.1.2. The Effect of Fixed Cost and Variable Cost
As our interest lies in the robustness of the final approach of the problem, we consider different fixed cost and variable cost values based on average over ten runs for each instance.
Take as an example. First let , and variable values are tested over the ranges . The total cost and number of vehicles are listed in Table 9.

From Table 9, we notice that the total cost increase as the fixed cost increases. The number of vehicles increases as the vehicle cost increases, but the increase is not obvious. It only changes slightly when the difference between and is large enough. The reason is that our heuristic algorithm always tends to accept the solution with the minimum number of vehicles with a given depots visiting sequence.
5.2. Compared Results between Models with and without Stochastic Speed
In order to observe the effects brought by the stochastic speeds, we focus on the compared results when . In the stochastic model, we set the means of average speed on all arcs to the given value in the first column of Table 10. In the fixed speed model, the speed of each arc is set to the given value, that is, the corresponding mean in the stochastic model. Table 10 lists the objectives obtained by the ISA algorithm.

As can be seen from Table 10, over the speed interval [5, 25], cost of model with fixed speed declines first; when it reaches 15 m/s, the minimum cost is obtained; then the cost begins to increase until the end of the speed interval.
In fact, fuel consumption takes up a large proportion in the total cost, so the trend of it (Figure 4) is in line with the trend of the total cost. Generally, costs obtained by the stochastic speed model is larger than that with fixed speed over the speed interval (the bold numbers in Table 10 are relatively high cost of the stochastic speed model), but when the speed is close to the end points of the interval, model with fixed speed will obtain smaller cost. The reason of this is that the integral interval will be ineffective once speeds exceed the two interval ends.
6. Conclusion
A new variant of the VRP with fuel consumption is presented, where the travel speed is considered to be stochastic. To solve this problem, an improved simulated annealing algorithm is proposed and presented in this paper. The proposed ISA is effective for solving the fuel VRP especially for largescale problems. Experimental results showed that when the value of the expected speed is not close to the end points of speed limit, total travel cost of the stochastic model is always larger than that of the fixed speed model. In order to make reliable logistics routing decisions, managers should take the stochasticity of speed into consideration.
Several extensions are possible for further research. One worth mentioning extension is when average travel speed follows other forms of distribution, and another one would be the time window version or traffic congestion version of this problem.
Appendix
A. The Pseudocode of the ISA Algorithm
See Algorithm 1.

B. Calculation of the Cost Function (VRPS)
See Algorithm 2.

Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The research of RenQian Zhang is supported by NSFC under 71571006 and 71271010, the research of Guozhu Jia is supported by NSFC under 71372007, and the research of Yanling Feng is supported by China Scholarship Council under 201606020056.
References
 J. K. Lenstra and A. H. G. R. Kan, “Complexity of vehicle routing and scheduling problems,” Networks, vol. 11, no. 2, pp. 221–227, 1981. View at: Publisher Site  Google Scholar
 G. B. Dantzig and J. H. Ramser, “The truck dispatching problem,” Management Science, vol. 6, no. 2, pp. 80–91, 1959. View at: Publisher Site  Google Scholar  MathSciNet
 P. Toth and D. Vigo, “Models, relaxations and exact approaches for the capacitated vehicle routing problem,” Discrete Applied Mathematics, vol. 123, no. 1–3, pp. 487–512, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 S. Kassem and M. Chen, “Solving reverse logistics vehicle routing problems with time windows,” International Journal of Advanced Manufacturing Technology, vol. 68, no. 14, pp. 57–68, 2013. View at: Publisher Site  Google Scholar
 R. Spliet and G. Desaulniers, “The discrete time window assignment vehicle routing problem,” European Journal of Operational Research, vol. 244, no. 2, pp. 379–391, 2015. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. Baldacci, A. Mingozzi, and R. Roberti, “Recent exact algorithms for solving the vehicle routing problem under capacity and time window constraints,” European Journal of Operational Research, vol. 218, no. 1, pp. 1–6, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. TavakkoliMoghaddam, A. R. Saremi, and M. S. Ziaee, “A memetic algorithm for a vehicle routing problem with backhauls,” Applied Mathematics and Computation, vol. 181, no. 2, pp. 1049–1060, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Hoff, I. Gribkovskaia, G. Laporte, and A. Løkketangen, “Lasso solution strategies for the vehicle routing problem with pickups and deliveries,” European Journal of Operational Research, vol. 192, no. 3, pp. 755–766, 2009. View at: Publisher Site  Google Scholar
 S. Allahyari, M. Salari, and D. Vigo, “A hybrid metaheuristic algorithm for the multidepot covering tour vehicle routing problem,” European Journal of Operational Research, vol. 242, no. 3, pp. 756–768, 2015. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 U. Ritzinger, J. Puchinger, and R. F. Hartl, “A survey on dynamic and stochastic vehicle routing problems,” International Journal of Production Research, vol. 54, no. 1, pp. 215–231, 2016. View at: Publisher Site  Google Scholar
 Y. Z. Mehrjerdi, “A multiple objective stochastic approach to vehicle routing problem,” The International Journal of Advanced Manufacturing Technology, vol. 74, no. 5–8, pp. 1149–1158, 2014. View at: Publisher Site  Google Scholar
 D. Taş, M. Gendreau, N. Dellaert, T. van Woensel, and A. G. de Kok, “Vehicle routing with soft time windows and stochastic travel times: a column generation and branchandprice solution approach,” European Journal of Operational Research, vol. 236, no. 3, pp. 789–799, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 D. Taş, N. Dellaert, T. van Woensel, and T. de Kok, “The timedependent vehicle routing problem with soft time windows and stochastic travel times,” Transportation Research Part C: Emerging Technologies, vol. 48, pp. 66–83, 2014. View at: Google Scholar
 J. F. Ehmke, A. M. Campbell, and T. L. Urban, “Ensuring service levels in routing problems with time windows and stochastic travel times,” European Journal of Operational Research, vol. 240, no. 2, pp. 539–550, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 G. Laporte, F. Louveaux, and H. Mercure, “Vehicle routing problem with stochastic travel times,” Transportation Science, vol. 26, no. 3, pp. 161–170, 1992. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y. Marinakis, G.R. Iordanidou, and M. Marinaki, “Particle Swarm Optimization for the vehicle routing problem with stochastic demands,” Applied Soft Computing Journal, vol. 13, no. 4, pp. 1693–1704, 2013. View at: Publisher Site  Google Scholar
 M. Dror, G. Laporte, and P. Trudeau, “Vehicle routing with stochastic demands: properties and solution frameworks,” Operations Research Society of America. Transportation Science Section. Transportation Science, vol. 23, no. 3, pp. 166–176, 1989. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. J. Bertsimas, “A vehicle routing problem with stochastic demand,” Operations Research, vol. 40, no. 3, pp. 574–585, 1992. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Zhang, W. H. Lam, and B. Y. Chen, “Ontime delivery probabilistic models for the vehicle routing problem with stochastic demands and time windows,” European Journal of Operational Research, vol. 249, no. 1, pp. 144–154, 2016. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Xiao and A. Konak, “A simulating annealing algorithm to solve the green vehicle routing & scheduling problem with hierarchical objectives and weighted tardiness,” Applied Soft Computing Journal, vol. 34, pp. 372–388, 2015. View at: Publisher Site  Google Scholar
 R. Kondekar, A. Gupta, G. Saluja, R. Maru, A. Rokde, and P. Deshpande, “A MapReduce based hybrid genetic algorithm using island approach for solving time dependent vehicle routing problem,” in Proceedings of the International Conference on Computer and Information Science (ICCIS '12), pp. 263–269, IEEE, Kuala Lumpur, Malaysia, June 2012. View at: Publisher Site  Google Scholar
 J. Zhang, X. Zhao, and X. He, “A minimum resource neural network framework for solving multiconstraint shortest path problems,” IEEE Transactions on Neural Networks and Learning Systems, vol. 25, no. 8, pp. 1566–1582, 2014. View at: Publisher Site  Google Scholar
 A. MeyerBäse, G. Botella, and L. RybarskaRusinek, “Stochastic stability analysis of competitive neural networks with different timescales,” Neurocomputing, vol. 118, pp. 115–118, 2013. View at: Publisher Site  Google Scholar
 E. Demir, T. Bektaş, and G. Laporte, “An adaptive large neighborhood search heuristic for the pollutionrouting problem,” European Journal of Operational Research, vol. 223, no. 2, pp. 346–359, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 T. Garaix, C. Artigues, D. Feillet, and D. Josselin, “Optimization of occupancy rate in dialaride problems via linear fractional column generation,” Computers & Operations Research, vol. 38, no. 10, pp. 1435–1442, 2011. View at: Publisher Site  Google Scholar
 Z. Cao, H. Guo, J. Zhang, D. Niyato, and U. Fastenrath, “Improving the efficiency of stochastic vehicle routing: a partial lagrange multiplier method,” IEEE Transactions on Vehicular Technology, vol. 65, no. 6, pp. 3993–4005, 2015. View at: Publisher Site  Google Scholar
 A. Ishigaki, “An application to stochastic vehiclerouting problem in a waste collection,” in Proceedings of the 5th IIAI International Congress on Advanced Applied Informatics (IIAIAAI '16), pp. 1151–1156, Kumamoto, Japan, July 2016. View at: Publisher Site  Google Scholar
 L. Schulze and A. Wüllner, “The approach of automated guided vehicle systems,” in Proceedings of the IEEE International Conference on Service Operations and Logistics, and Informatics (SOLI '06), pp. 522–527, IEEE, Shanghai, China, June 2006. View at: Publisher Site  Google Scholar
 G. Botella, A. García, M. RodríguezÁlvarez, E. Ros, U. MeyerBaese, and M. C. Molina, “Robust bioinspired architecture for opticalflow computation,” IEEE Transactions on Very Large Scale Integration (VLSI) Systems, vol. 18, no. 4, pp. 616–629, 2010. View at: Publisher Site  Google Scholar
 S. Mafrica, A. Servel, and F. Ruffier, “Towards an automatic parking system using bioinspired 1D optical flow sensors,” in Proceedings of the IEEE International Conference on Vehicular Electronics and Safety (ICVES '15), pp. 96–103, Yokohama, Japan, November 2015. View at: Publisher Site  Google Scholar
 Y. Kuo, “Using simulated annealing to minimize fuel consumption for the timedependent vehicle routing problem,” Computers & Industrial Engineering, vol. 59, no. 1, pp. 157–165, 2010. View at: Publisher Site  Google Scholar
 J. Qian and R. Eglese, “Finding least fuel emission paths in a network with timevarying speeds,” Networks, vol. 63, no. 1, pp. 96–106, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 M. Barth and K. Boriboonsomsin, “Energy and emissions impacts of a freewaybased dynamic ecodriving system,” Transportation Research Part D: Transport and Environment, vol. 14, no. 6, pp. 400–410, 2009. View at: Publisher Site  Google Scholar
 T. Bektaş and G. Laporte, “The pollutionrouting problem,” Transportation Research Part B: Methodological, vol. 45, no. 8, pp. 1232–1250, 2011. View at: Publisher Site  Google Scholar
 Y. Xiao, Q. Zhao, I. Kaku, and Y. Xu, “Development of a fuel consumption optimization model for the capacitated vehicle routing problem,” Computers and Operations Research, vol. 39, no. 7, pp. 1419–1431, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Figliozzi, “Vehicle routing problem for emissions minimization,” Transportation Research Record, no. 2197, pp. 1–7, 2010. View at: Publisher Site  Google Scholar
 F. Errico, G. Desaulniers, M. Gendreau, W. Rei, and L.M. Rousseau, “A priori optimization with recourse for the vehicle routing problem with hard time windows and stochastic service times,” European Journal of Operational Research, vol. 249, no. 1, pp. 55–66, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 M. Marinaki and Y. Marinakis, “A glowworm swarm optimization algorithm for the vehicle routing problem with stochastic demands,” Expert Systems with Applications, vol. 46, pp. 145–163, 2016. View at: Publisher Site  Google Scholar
 B. Sarasola, K. F. Doerner, V. Schmid, and E. Alba, “Variable neighborhood search for the stochastic and dynamic vehicle routing problem,” Annals of Operations Research, vol. 236, no. 2, pp. 425–461, 2016. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 O. Jabali, T. van Woensel, and A. G. de Kok, “Analysis of travel times and CO_{2} emissions in timedependent vehicle routing,” Production and Operations Management, vol. 21, no. 6, pp. 1060–1074, 2012. View at: Publisher Site  Google Scholar
 A. Franceschetti, D. Honhon, T. Van Woensel, T. Bektaş, and G. Laporte, “The timedependent pollutionrouting problem,” Transportation Research Part B: Methodological, vol. 56, pp. 265–293, 2013. View at: Publisher Site  Google Scholar
 Y. Xiao and K. Abdullah, “A genetic algorithm with exact dynamic programming for the green vehicle routing & scheduling problem,” Journal of Cleaner Production, 2016. View at: Publisher Site  Google Scholar
 V. Granville, M. Krivánek, and J.P. Rasson, “Simulated annealing: a proof of convergence,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 16, no. 6, pp. 652–656, 1994. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Yanling Feng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.