^{1}

^{1}

^{1}

This study attempts to solve vehicle routing problem with time window (VRPTW). The study first identifies the real problems and suggests some recommendations on the issues. The technique used in this study is Genetic Algorithm (GA) and initialization applied is random population method. The objective of the study is to assign a number of vehicles to routes that connect customers and depot such that the overall distance travelled is minimized and the delivery operations are completed within the time windows requested by the customers. The analysis reveals that the problems experienced in vehicle routing with time window can be solved by GA and retrieved for optimal solutions. After a thorough study on VRPTW, it is highly recommended that a company should implement the optimal routes derived from the study to increase the efficiency and accuracy of delivery with time insertion.

In daily distribution business activities, no company is free from vehicle routing problem especially in terms of searching the minimum distance to deliver the goods. The travelled distance related to people and goods is very significant. It is motivated by the continuously increasing business complexity experienced today [

In today’s world, everything is desired by the customers just at one click of a button and the requirement of products becomes so fast. Therefore, delivering products faster and in efficient ways have become very important. The objective of delivering products at a minimum cost (as highlighted earlier, in this study, we consider travelled distance as the cost) involves several other factors. The factors are routes taken, capacity of the vehicles, requirement and constraints sets by the customers and many more. To achieve the objective, several models have been introduced in the previous years.

Vehicle Routing Problem with Time Windows (VRPTW) involves several customers (nodes) to be served by a fleet of vehicles within time interval and each vehicle has its limited capacity. The vehicles deliver the products from the depot to the customers and the trip ends at same depot. The size of the delivered products is identical and consumes same amount of capacity of each vehicle. In addition, the delivery locations are unique and any activities of feeding customers at a node other than main depot are not allowed. The service level efficiency depends on the time required by the customers. Thus, the aim of the study is to minimize the total travel distance of the fleet vehicle every time it travels to deliver the items requested by each customer. Travelling at the minimum route can reduce cost and also increase the efficiency of service time [

Vehicle Routing Problem (VRP) is described as the designation of least cost routes from a central depot to a set of geographically dispersed points with various demands [

Since year 1999, the methods for solving VRPTW have been developed and the most efficient method is two-phased hybrid algorithms. Two-phase hybrid algorithms are divided into two stages which are the minimization of routes number and minimization of travel cost. Basically, the two-phase hybrid algorithm uses the local search procedures to minimize the routes and then minimize the total cost of vehicle routing. Gehring and Homberger [

Genetic Algorithm (GA) has been popular due to its contribution in obtaining good solutions for complicated optimization problems in a reasonable amount of time. In particular, heuristic search strategies based on GA have been explored in recent years in order to improve the solutions in VRPTW problems. GA deals with a population of solutions subjected to reproductive processes. Those processes occur under the direct influence of each individual’s relative fitness. The fitness of each individual of a population is defined by the quality of its solution with respect to the objective under consideration. Such a setup contributes to increasing the average fitness of a population in each subsequent generation. After a certain number of generations, we are likely to get a solution of the quality we need. The applicability of GA is spread over a wide variety of fields, such as single to multi-objective optimization, engineering, bin packing, multidimensional numerical optimization, pattern recognition and routing. However, GA has obtained very limited success in the area of vehicle routing, due to their high computational cost [

In order to ensure that the data is suitable for the comparison purpose, a model has been adapted. Together with the model, several assumptions have been identified to describe the model [

1) There is a demand required by each and every customer in the network.

2) All of the capacities of the vehicles used to transport the goods are identical and limited.

3) The customers have to be served by the fixed number of vehicles.

4) Vehicles have to leave the depot carrying the amount of goods equal to the amount it must deliver to the customers.

5) Time window is given in two values described as lower bound (start time of service at visited customer by respected vehicle) and upper bound (finish time of service at visited customer by respected vehicle).

6) Vehicle is not supposed to begin service at a customer once it crosses the latest time limit.

7) If the vehicle reaches before the earliest time limit, the respected route is rejected and search for new route is analyzed.

8) Service time is allotted to each vehicle for loading and unloading goods to each customer in the network.

9) The total route time of each vehicle is calculated as the sum of travel time (proportional to travel distance), waiting time and service time.

10) The maximum route travelled by a vehicle must be higher than the total route travel distance.

The following are the notations used to describe the model:

The VRPTW model is adopted from [

Subject to

The objective function of VRPTW model is to minimize the distance travel. Constraint 1) ensures that each customer is visited by a vehicle while Constraint 2) ensures that same vehicle arrives and departs from each customer it serves. Constraint 3) is to ensure that at most k vehicles are used and Constraint 4) indicates the flow equations for dynamic demands. Constraint 5) is the flow equation for delivery demands and Constraint 6) is to make sure that demands will only transported in the exact route. Constraint 7) is time window constraint and Constraint 8) defines the nature of the decision variables.

There are many studies in solving VRPTW using GA that have been done. However, in this study the application does not only concentrate on the trial solution but also generating the initial solution. In doing so, random population initialization (GA/R) is used to generate candidate solutions. In addition, the fitness function of this study is the summation of the traveling distance of all the routes as discussed in [

The modified Genetic Algorithm in VRPTW consists of three phases namely: 1) initialization of population and chromosome representation; 2) selection and; 3) crossover and mutation. All three phases are discussed in detail as follows.

It is the first stage where GA starts with a set of chromosomes. Initial population is randomly generated in the case where the time taken for converging to the solution is longer. The initial population size used is 100. It defines the size of population at each generation in a GA. Hence 100 sets of random integer N is created where N is the number of customers in the problem. The solution space may not be adequately explored by the algorithm, if the population size is too small. Increasing the population size enables the GA to search more points thereby obtain a better result. However, the larger the population size, the longer the GA takes to compute each generation [

In this study, chromosome representation is in the form of network configuration given by an integer string of length N. A gene in a given chromosome is the point assigned to a customer. The sequence of genes in chromosome string is the order of visitations of the customers. An example of a chromosome resulting in a solution for the network is given in

Selection function is used to select the prospective parents based on evaluation of the fitness function. Then, the selected parent chromosomes are recombined by the crossover operator for creating new population. This study uses a tournament selection in which two identical copies of the population are kept. For each generation, the adjacent chromosomes in one copy of the population are compared pair by pair, and then the chromosome with greater fitness value is selected. Tournament selection is applied because genetically fitter chromosomes are given priority in this mechanism [

Single point crossover operator that swaps the whole portion in

To prevent the premature convergence by involving the random changes of the segments of the parent string,

this study uses swapping mutation where the selection is random for the two segments in the string and then swaps them.

Iteration process will run for hundred times before terminated.

The analysis is divided into two phases; the first phase analysis is on VRP without time windows and second phase of the analysis is on VRPTW. In phase one, two sets of data comprising small and medium sizes are used to validate the performance of the proposed algorithm (random population initialization (GA/R)) in solving VRP. In phase two, the same sets of data are used to validate the performance of the proposed algorithm (random population initialization (GA/R)) in solving VRPTW. The comparison is made with other type of initializations, GA/NN and GA/GRASP. Below are the conditions of each parameter:

1) Population size (Popsize) defined as the number of chromosomes in each generation (10, 20, 40 and 60).

2) Maximum generation (Maxgen) which is a termination criterion that sets the maximum number of chromosome populations generated before the top scoring chromosome return as the search answer (100).

3) Crossover probability is the probability that a pair of chromosomes cross (0.95).

4) Mutation probability is the probability that a gene on a chromosome randomly mutate (0.01).

In this benchmarking phase, the results indicate that the performance of GA/R is better compared to other competing heuristics. Moreover, GA/R is also the fastest in computing the results. By observing the shortest distance in all three algorithms and different problem sizes, GA/R provides excellent results by producing the smallest standard deviation. For GA/R, the smallest computational time of 2.49 occurs when population size used is 10.

The validation phase proved that the proposed method is capable to be applied for small and medium size data. In this phase, the analysis is carried out on two types of data with the insertion of time windows.

Population Size | Shortest Distance (km) | Standard Deviation | Average Distance (km) | Average Computational Time (seconds) | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|

GA/ NN | GA/ GRASP | GA/R | GA/ NN | GA/ GRASP | GA/R | GA/ NN | GA/ GRASP | GA/R | GA/ NN | GA/ GRASP | GA/R | |

10 | 790.00 | 790.00 | 730.00 | 34.35 | 34.35 | 6.58 | 801.23 | 801.23 | 759.50 | 10.50 | 2.40 | 2.49 |

20 | 810.00 | 790.00 | 730.00 | 36.76 | 34.35 | 4.48 | 824.93 | 801.23 | 753.00 | 10.50 | 2.08 | 2.66 |

40 | 795.00 | 765.00 | 730.00 | 32.69 | 30.31 | 4.14 | 801.18 | 784.60 | 741.00 | 16.88 | 1.44 | 2.94 |

60 | 790.00 | 765.00 | 730.00 | 34.35 | 30.31 | 3.15 | 801.18 | 784.60 | 736.00 | 23.45 | 1.01 | 3.51 |

Population Size | Shortest Distance (km) | Standard Deviation | Average Distance (km) | Average Computational Time (seconds) | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|

GA/NN | GA/ GRASP | GA/R | GA/ NN | GA/ GRASP | GA/R | GA/NN | GA/ GRASP | GA/R | GA/ NN | GA/ GRASP | GA/R | |

10 | 1925.00 | 1925.00 | 1655.00 | 37.91 | 40.28 | 22.03 | 1931.98 | 1935.88 | 1755.00 | 1.94 | 1.55 | 2.45 |

20 | 1815.00 | 1925.00 | 1660.00 | 70.31 | 42.56 | 17.39 | 1837.35 | 1940.75 | 1801.50 | 6.50 | 3.21 | 2.52 |

40 | 1815.00 | 1815.00 | 1645.00 | 70.31 | 74.44 | 15.07 | 1837.35 | 1839.30 | 1758.50 | 15.83 | 2.36 | 3.76 |

60 | 1780.00 | 1805.00 | 1645.00 | 106.23 | 70.85 | 15.89 | 1816.38 | 1835.90 | 1645.00 | 31.16 | 2.11 | 4.65 |

distance is higher (795 km). For the standard deviation, mixed results are found between GA/R and GA/N in which, when only population sizes are 10 and 40, GA/R outperforms GA/NN. In term of average distance GA/R also did well as it has the smallest average distance contributed by population size 40 (783 km) compared to its competitors. However, for bigger population size, the initialization based on GA/R needs to be improved.

From the analysis, GA/R has performed well where it generates smaller distance as well as computational time. The best route obtained from GA/R initialization can be used and applied to the vehicle routing problem either with or without time window. The analysis has proved that the VRPTW can be solved by the GA/R and this may helps organizations to consider of using the routes in order to reduce distance travel.

Population Size | Shortest Distance (km) | Standard Deviation | Average Distance (km) | Average Computational time (seconds) | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|

GA/ NN | GA/ GRASP | GA/R | GA/ NN | GA/ GRASP | GA/R | GA/NN | GA/ GRASP | GA/R | GA/NN | GA/ GRASP | GA/R | |

10 | 795.00 | 795.00 | 780.00 | 22.98 | 22.98 | 13.34 | 800.63 | 800.63 | 830.00 | 3.26 | 2.48 | 3.42 |

20 | 795.00 | 795.00 | 765.00 | 6.50 | 22.98 | 6.85 | 795.65 | 800.63 | 809.00 | 4.55 | 5.10 | 4.86 |

40 | 795.00 | 795.00 | 765.00 | 11.08 | 22.98 | 4.14 | 795.93 | 800.63 | 783.00 | 4.72 | 8.41 | 7.63 |

60 | 795.00 | 795.00 | 765.00 | 11.08 | 22.98 | 25.33 | 795.93 | 800.63 | 866.50 | 5.14 | 10.19 | 7.49 |

Population Size | Shortest Distance (km) | Standard Deviation | Average Distance (km) | Average Computational Time (seconds) | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|

GA/NN | GA/ GRASP | GA/R | GA/NN | GA/ GRASP | GA/R | GA/NN | GA/ GRASP | GA/R | GA/NN | GA/ GRASP | GA/R | |

10 | 1395.00 | 1355.00 | 1705.00 | 64.22 | 40.21 | 15.45 | 1431.93 | 1385.65 | 1813.00 | 189.80 | 190.28 | 105.31 |

20 | 1395.00 | 1360.00 | 1660.00 | 64.22 | 46.52 | 17.41 | 1431.93 | 1373.48 | 1748.00 | 189.80 | 210.75 | 110.82 |

40 | 1320.00 | 1355.00 | 1600.00 | 68.51 | 51.99 | 19.01 | 1359.83 | 1399.65 | 1764.50 | 168.00 | 109.30 | 104.53 |

60 | 1345.00 | 1340.00 | 1645.00 | 63.93 | 46.55 | 17.07 | 1379.53 | 1384.58 | 1755.00 | 168.41 | 108.56 | 156.82 |

Problem Size | Distance (km) | Route |
---|---|---|

Small (8 Customers and 3 Vehicles) | 765.00 | Vehicle 1: 0-6-7-4-0 Vehicle 2: 0-3-5-1-0 Vehicle 3: 0-8-2-0 |

Medium (25 Customers and 5 Vehicles) | 1600.00 | Vehicle 1: 0-1-2-7-4-5-6-3-0 Vehicle 2:0-8-9-10-11-12-13-0 Vehicle 3:0-14-15-16-17-18-0 Vehicle 4:0-19-20-21-22-23-0 Vehicle 5: 0-24-25-0 |

Large (70 Customers and 8 Vehicles) | 9225.00 | Vehicle 1: 0-20-43-40-4-57-1-2-53-66-24-0 Vehicle 2: 0-47-22-19-49-35-25-21-6-0 Vehicle 3: 0-61-51-11-13-23-27-17-0 Vehicle 4: 0-36-58-46-18-16-14-64-12-0 Vehicle 5: 0-26-5-63-31-10-59-44-45-41-38-0 Vehicle 6: 0-52-28-39-42-33-8-69-7-48-0 Vehicle 7:0-65-68-70-62-29-67-54-32-34-37-50-0 Vehicle 8: 0-55-56-9-3-60-30-15-0 |

Population Size | ||||
---|---|---|---|---|

10 | 20 | 40 | 60 | |

Travelled Distance (km) | 37.50 | 36.40 | 36.40 | 37.50 |

Route | 0-2-8-4-5-1-3-7-6-0 | 0-2-8-5-4-1-3-7-6-0 | 0-2-8-5-4-1-3-7-6-0 | 0-8-4-5-2-1-3-7-6-0 |

Depot-Pandamaran-Crystal Crown-Jalan Kem-Jalan Limbongan-Kg Idaman-Bandar Botanik-Tmn Pendamar Indah 2-Tmn Teluk Gedung Indah-Depot | Depot-Pandamaran-Crystal Crown-Jalan Limbongan-Jalan Kem-Kg Idaman-Bandar Botanik-Tmn Pendamar Indah 2-Tmn Teluk Gedung Indah-Depot | Depot-Pandamaran-Crystal Crown-Jalan Limbongan-Jalan Kem-Kg Idaman-Bandar Botanik-Tmn Pendamar Indah 2-Tmn Teluk Gedung Indah-Depot | Depot-Crystal Crown- Jalan Kem-Jalan Limbongan-Pandamaran- Kg Idaman-Bandar Botanik-Tmn Pendamar Indah 2-Tmn Teluk Gedung Indah-Depot | |

Standard Deviation | 2.56 | 1.74 | 2.45 | 1.24 |

Computational Time (second) | 2.15 | 2.04 | 2.62 | 2.30 |

Nowadays, logistics and transportation are needed for daily delivery in order to deliver the items demanded by customers. In daily delivery, they practically experience many types of the vehicle routing problems. Par- ticularly, the problem is based on the distribution time which does not meet the time that is required by customers. The problem that considers the distribution time is specified as the vehicle routing problem within time interval or VRPTW and is designed to solve and increase the company’s delivery efficiency in future. This is a case study of a district in Selangor, Malaysia in which the company is supplying ice cubes for several cus- tomers within the area. We capture and analyze the problem in details and report the solutions.

Based on the results in

We would like to acknowledge Ministry of Higher Education Malaysia and Research Management Institute (RMI) of Universiti Teknologi MARA, Malaysia, for funding this research through Grant No: 600-RMI/FRGS 5/3 (9/2013).

Nor Edayu Abdul Ghani,S. Sarifah Radiah Shariff,Siti Meriam Zahari, (2016) An Alternative Algorithm for Vehicle Routing Problem with Time Windows for Daily Deliveries. Advances in Pure Mathematics,06,342-350. doi: 10.4236/apm.2016.65025