Vehicle routing problem explained

The vehicle routing problem (VRP) is a combinatorial optimization and integer programming problem which asks "What is the optimal set of routes for a fleet of vehicles to traverse in order to deliver to a given set of customers?" It generalises the travelling salesman problem (TSP). It first appeared in a paper by George Dantzig and John Ramser in 1959,[1] in which the first algorithmic approach was written and was applied to petrol deliveries. Often, the context is that of delivering goods located at a central depot to customers who have placed orders for such goods. The objective of the VRP is to minimize the total route cost. In 1964, Clarke and Wright improved on Dantzig and Ramser's approach using an effective greedy algorithm called the savings algorithm.

Determining the optimal solution to VRP is NP-hard,[2] so the size of problems that can be optimally solved using mathematical programming or combinatorial optimization may be limited. Therefore, commercial solvers tend to use heuristics due to the size and frequency of real world VRPs they need to solve.

VRP has many direct applications in industry. Vendors of VRP routing tools often claim that they can offer cost savings of 5%–30%.[3]

Setting up the problem

The VRP concerns the service of a delivery company. How things are delivered from one or more depots which has a given set of home vehicles and operated by a set of drivers who can move on a given road network to a set of customers. It asks for a determination of a set of routes, S, (one route for each vehicle that must start and finish at its own depot) such that all customers' requirements and operational constraints are satisfied and the global transportation cost is minimized. This cost may be monetary, distance or otherwise.

The road network can be described using a graph where the arcs are roads and vertices are junctions between them. The arcs may be directed or undirected due to the possible presence of one way streets or different costs in each direction. Each arc has an associated cost which is generally its length or travel time which may be dependent on vehicle type.

To know the global cost of each route, the travel cost and the travel time between each customer and the depot must be known. To do this our original graph is transformed into one where the vertices are the customers and depot, and the arcs are the roads between them. The cost on each arc is the lowest cost between the two points on the original road network. This is easy to do as shortest path problems are relatively easy to solve. This transforms the sparse original graph into a complete graph. For each pair of vertices i and j, there exists an arc (i,j) of the complete graph whose cost is written as

Cij

and is defined to be the cost of shortest path from i to j. The travel time

tij

is the sum of the travel times of the arcs on the shortest path from i to j on the original road graph.

Sometimes it is impossible to satisfy all of a customer's demands and in such cases solvers may reduce some customers' demands or leave some customers unserved. To deal with these situations a priority variable for each customer can be introduced or associated penalties for the partial or lack of service for each customer given [2]

The objective function of a VRP can be very different depending on the particular application of the result but a few of the more common objectives are:[2]

VRP variants

Several variations and specializations of the vehicle routing problem exist:

Several software vendors have built software products to solve various VRP problems. Numerous articles are available for more detail on their research and results.

Although VRP is related to the Job Shop Scheduling Problem, the two problems are typically solved using different techniques.[9]

Exact solution methods

There are three main different approaches to modelling the VRP

  1. Vehicle flow formulations—this uses integer variables associated with each arc that count the number of times that the edge is traversed by a vehicle. It is generally used for basic VRPs. This is good for cases where the solution cost can be expressed as the sum of any costs associated with the arcs. However it can't be used to handle many practical applications.
  2. Commodity flow formulations—additional integer variables are associated with the arcs or edges which represent the flow of commodities along the paths travelled by the vehicles. This has only recently been used to find an exact solution.
  3. Set partitioning problem—These have an exponential number of binary variables which are each associated with a different feasible circuit. The VRP is then instead formulated as a set partitioning problem which asks what is the collection of circuits with minimum cost that satisfy the VRP constraints. This allows for very general route costs.

Vehicle flow formulations

The formulation of the TSP by Dantzig, Fulkerson and Johnson was extended to create the two index vehicle flow formulations for the VRP

min\sumi\in\sumjcijxij

subject to

In this formulation

cij

represents the cost of going from node

i

to node

j

,

xij

is a binary variable that has value

1

if the arc going from

i

to

j

is considered as part of the solution and

0

otherwise,

K

is the number of available vehicles and

r(S)

corresponds to the minimum number of vehicles needed to serve set

S

. We are also assuming that

0

is the depot node.

Constraints and state that exactly one arc enters and exactly one leaves each vertex associated with a customer, respectively. Constraints and say that the number of vehicles leaving the depot is the same as the number entering. Constraints are the capacity cut constraints, which impose that the routes must be connected and that the demand on each route must not exceed the vehicle capacity. Finally, constraints are the integrality constraints.

One arbitrary constraint among the

2|V|

constraints is actually implied by the remaining

2|V|-1

ones so it can be removed. Each cut defined by a customer set

S

is crossed by a number of arcs not smaller than (minimum number of vehicles needed to serve set

S

).

An alternative formulation may be obtained by transforming the capacity cut constraints into generalised subtour elimination constraints (GSECs).

\sumi\in\sumj\inxij\leq|S|-r(S)

which imposes that at least arcs leave each customer set

S

.

GCECs and CCCs have an exponential number of constraints so it is practically impossible to solve the linear relaxation. A possible way to solve this is to consider a limited subset of these constraints and add the rest if needed. Identification of the needed constraints is done via a separation procedure. Efficient exact separation methods for such constraints (based on mixed integer programming) have been developed.[10]

A different method again is to use a family of constraints which have a polynomial cardinality which are known as the MTZ constraints, they were first proposed for the TSP [11] and subsequently extended by Christofides, Mingozzi and Toth.[12]

uj-ui\geqdj-C(1-xij)~~~~~~\foralli,j\inV\backslash\{0\},ij~~~~s.t.di+dj\leqC

0\lequi\leqC-di~~~~~~\foralli\inV\backslash\{0\}

where

ui,~i\inV\backslash\{0\}

is an additional continuous variable which represents the load left in the vehicle after visiting customer

i

and

di

is the demand of customer

i

. These impose both the connectivity and the capacity requirements. When

xij=0

constraint then

i

is not binding' since

ui\leqC

and

uj\geqdj

whereas

xij=1

they impose that

uj\gequi+dj

.

These have been used extensively to model the basic VRP (CVRP) and the VRPB. However, their power is limited to these simple problems. They can only be used when the cost of the solution can be expressed as the sum of the costs of the arc costs. We cannot also know which vehicle traverses each arc. Hence we cannot use this for more complex models where the cost and or feasibility is dependent on the order of the customers or the vehicles used.

Manual versus automatic optimum routing

There are many methods to solve vehicle routing problems manually. For example, optimum routing is a big efficiency issue for forklifts in large warehouses. Some of the manual methods to decide upon the most efficient route are: Largest gap, S-shape, Aisle-by-aisle, Combined and Combined + While Combined + method is the most complex, thus the hardest to be used by lift truck operators, it is the most efficient routing method. Still the percentage difference between the manual optimum routing method and the real optimum route was on average 13%.[13] [14]

Metaheuristic

Due to the difficulty of solving to optimality large-scale instances of vehicle routing problems, a significant research effort has been dedicated to metaheuristics such as Genetic algorithms, Tabu search, Simulated annealing and Adaptive Large Neighborhood Search (ALNS). Some of the most recent and efficient metaheuristics for vehicle routing problems reach solutions within 0.5% or 1% of the optimum for problem instances counting hundreds or thousands of delivery points.[15] These methods are also more robust in the sense that they can be more easily adapted to deal with a variety of side constraints. As such, the application of metaheuristic techniques is often preferred for large-scale applications with complicating constraints and decision sets.

See also

Further reading

Notes and References

  1. Dantzig. George Bernard. Ramser, John Hubert . The Truck Dispatching Problem. Management Science. October 1959. 6. 1. 80–91. 10.1287/mnsc.6.1.80.
  2. Book: Toth, P. . Vigo, D. . The Vehicle Routing Problem . 9 . Monographs on Discrete Mathematics and Applications . 2002 . Society for Industrial and Applied Mathematics . Philadelphia . 0-89871-579-2.
  3. Book: Geometric Modelling, Numerical Simulation, and Optimization:: Applied Mathematics at SINTEF . 2007 . Springer Verlag . 978-3-540-68783-2 . Geir Hasle . Berlin . 397–398 . Knut-Andreas Lie . Ewald Quak.
  4. Chao . I-Ming . Golden . Bruce L . Wasil . Edward A . The Team Orienteering Problem . European Journal of Operational Research . 1996 . 88 . 3 . 464–474 . 10.1016/0377-2217(94)00289-4.
  5. Book: Archetti, C.. Sperenza, G.. Vigo, D. . Toth, P. . Vigo, D.. 2014. Vehicle routing problems with profits. Vehicle Routing: Problems, Methods, and Applications . Second . 10.1137/1.9781611973594.ch10. 273–297.
  6. Hammami. Farouk . Rekik . Monia . Coelho . Leandro C. . A hybrid adaptive large neighborhood search heuristic for the team orienteering problem . Computers & Operations Research . 2020 . 123 . 105034 . 10.1016/j.cor.2020.105034. 221134904 .
  7. Ekici . Ali . Özener . Okan Örsan . Kuyzu . Gültekin . Cyclic Delivery Schedules for an Inventory Routing Problem . Transportation Science . November 2015 . 49 . 4 . 817–829 . 10.1287/trsc.2014.0538.
  8. Mahmud . Nafix . Haque . Md. Mokammel . Solving Multiple Depot Vehicle Routing Problem (MDVRP) using Genetic Algorithm . 2019 International Conference on Electrical, Computer and Communication Engineering (ECCE) . February 2019 . 10.1109/ECACE.2019.8679429.
  9. Beck, J.C. . Prosser, P. . Selensky, E. . 2003 . Vehicle routing and job shop scheduling: What's the difference? . Proceedings of the 13th International Conference on Artificial Intelligence Planning and Scheduling .
  10. Pavlikov . K. . Petersen . N. C. . Sørensen . J. L. . 2023 . Exact Separation of the Rounded Capacity Inequalities for the CapacitatedVehicle Routing Problem . Networks . 83 . 197–209 . 10.1002/net.22183 . 263321558 . Networks. free .
  11. Miller. C. E.. Tucker. E. W.. Zemlin. R. A.. Integer Programming Formulations and Travelling Salesman Problems. J. ACM. 1960. 7. 326–329. 10.1145/321043.321046. 2984845. free.
  12. Book: Christofides. N.. Mingozzi. A.. Toth. P.. The Vehicle Routing Problem. 1979. Wiley. Chichester, UK. 315–338.
  13. Web site: Why Is Manual Warehouse Optimum Routing So Inefficient? . Locatible.com . 2016-09-26 . 2016-09-26.
  14. Web site: Roodbergen . Kees Jan . Routing methods for warehouses with multiple cross aisles . roodbergen.com . 2001 . 2016-09-26.
  15. Vidal T, Crainic TG, Gendreau M, Prins C. 2014. A unified solution framework for multi-attribute vehicle routing problems. European Journal of Operational Research. 234. 3. 658–673. 10.1016/j.ejor.2013.09.045. 21037953 .