Impact of travel time uncertainties on the solution cost of a two-echelon vehicle routing problem with synchronization

Two-echelon vehicle routing problems which contain synchronization between vehicles can be deeply impacted by time uncertainty, because one vehicle’s delay can propagate to other vehicles. In this paper, we evaluate the deterministic solution of such a problem based on simulated travel time scenarios. The information obtained by simulation is incorporated in the optimization procedure iteratively. Computational results show that the degree of synchronization in an instance is directly correlated with the potential improvements by reoptimization. We present findings on the number of travel time scenarios required to obtain a representative picture of the stochastic solutions. In addition, we demonstrate that time dependent travel times can be aggregated on a city-wide level and linearized as a function of free flow times without major loss of reliability.


Introduction
Emission-free, agile vehicles are increasingly used for freight deliveries within cities and especially within historic city centers.This development is triggered by the ongoing urbanization as well as the negative effects of climate change caused by greenhouse gas emissions.Also negative health effects in cities caused by other transport-related emissions (WHO 2016) contribute to this trend.However, innovative logistic concepts are required to enable the efficient use of these vehicles within a distribution network.Since the loading capacity and the operating distance of these vehicles are lower, certain limitations have to be taken into consideration.Either consolidation centers near the city center are required, where those vehicles can reload goods (Quak et al. 2014), or a reorganized distribution scheme with synchronization between different types of vehicles can be used to allow for reloading.
The latter approach does not require additional storage facilities, but deals with the problem of reloading agile emission-free vehicles by synchronized meetings with possibly conventional vehicles of larger capacity and operating distance.The conventional vehicles are used for transporting goods from a depot to customers located outside the city center as well as to supply the agile emission-free vehicles at so-called satellites.Satellites are easily accessible spaces which do not offer any storage possibilities.Therefore, temporal and spatial synchronization between reloading vehicles is required.One main characteristic of a solution to such a problem is its inherent interdependence, which makes it sensitive to uncertainties.
The aim of this paper is to investigate the impact of stochastic travel times on the cost of a deterministic solution of a two-echelon vehicle routing problem (2eVRP) with spatial and temporal synchronization.Additionally, we explore how the deterministic solution can be improved with respect to stochasticity by exploiting the available stochastic information.Therefore, we build on an existing solution algorithm for the deterministic problem (Anderluh et al. 2017) and add additional buffer time at time-critical nodes based on information we gain by simulating travel time scenarios.
The contribution of this paper is twofold.On the one hand it provides first insights into the quality of a deterministic solution to a complex VRP with synchronization constraints when applied to a real-life stochastic setting.On the other hand ways of improving the deterministic solution by exploiting stochastic information are evaluated.
The remainder of the paper is organized as follows: Sect. 2 reviews the literature related to stochasticity in routing problems, whereas Sect. 3 describes the problem at hand and Sect. 4 focuses on the solution procedure in use.In Sect. 5 computational results are presented and Sect.6 concludes the paper.

Literature review
An early work on stochastic vehicle routing can be found in Stewart and Golden (1983).In the 1990s Bertsimas and van Ryzin (1991), Bertsimas (1992) and Bastian and Rinnooy Kan (1992) continued dealing with this class of problems.Gendreau et al. (1996) provided a first literature review on the topic.A recent literature review by Braekers et al. (2016) about VRPs shows that only 58 out of 277 articles published between 2009 and 2015 focus on the stochastic aspects.The majority of these papers-34 articles-deal with stochasticity in demand and dynamic requests.Some of these papers deal with specific solution methods like a branch-cut-and-price algorithm (Gauvin et al. 2014) or a self-adaptive robust search (Chen et al. 2012).Others focus on additional constraints like time windows (Lei et al. 2011), pickup and delivery (Minis and Tatarakis 2011), weighted related cost (Luo et al. 2016) or maximum tour duration (Erera et al. 2010).Juan et al. (2011) generate different scenarios for the VRP with stochastic demand to solve the resulting deterministic problems efficiently with existing solvers.Only 14 articles included in the review by Braekers et al. (2016) deal with stochastic or unknown travel times, only 10 papers focus on time-dependent travel times and none considers stochastic as well as time-dependent travel times at once.
Papers addressing stochastic and dynamic vehicle routing problems of the last three decades have been recently surveyed by Psaraftis et al. (2016), who consider 117 papers in total in their review.10 papers deal with travel times as dynamic element and 8 papers take travel times as stochastic element into account.Only one out of these papers considers travel times as stochastic and dynamic element simultaneously.Based on Pillac et al. (2013), the difference between stochastic and dynamic problems can be defined as such: Dynamic problems refer to information evolution-i.e. when do we get the information.Stochastic problems refer to information quality-i.e.how reliable is the information at hand.
VRPs solved with stochastic travel times typically take into account time window restrictions as shown in Taş et al. (2014) who solve such a problem with an iterated local search.In some cases, in addition to stochastic travel times and time window restrictions, also stochastic service times are considered (Li et al. 2010;Zhang et al. 2013).Ehmke et al. (2015) focus additionally on the service level of the solution.Besides, Zhang et al. (2012) solve a VRP with simultaneous pickup and delivery while taking into account stochastic travel times, whereas Chen et al. (2014) focus on an arc routing problem with stochastic travel and service times.
The two-echelon capacitated VRP has already been investigated in detail as shown by Cuda et al. (2015) in their survey as well as by Cattaruzza et al. (2017).Synchronization aspects in a deterministic 2eVRP are tackled by Grangier et al. (2016).Stochastic information in a 2eVRP is only taken into account by Tavakkoli-Moghaddam and Raziei (2016) who consider fuzzy demand for multiple products in multiple periods in a two-echelon bi-objective distribution network.To our knowledge there is no paper dealing with stochastic travel times in a 2eVRP.
When dealing with stochastic travel times, one important issue is the probability distribution of such travel times.This distribution can either be determined based on historic empirical data (Aron et al. 2014;Rahmani et al. 2015) or by using 1 3 Impact of travel time uncertainties on the solution cost of… traffic-related factors like the reduction in traffic capacities because of incidents as well as daily variation in travel activity (Boyles et al. 2010).Although different probability distributions have already been tested in vehicle routing (Gomez et al. 2016;Aron et al. 2014), in the majority of stochastic VRP models a lognormal distribution is assumed as a good approximation (Lecluyse et al. 2009).
In contrast to stochastic travel times, time-dependent travel times can also be taken into account.Gendreau et al. (2017) just recently edit a special section on time-dependent routing in which Montero et al. (2017) focus on a time-dependent TSP with time windows while Mancini (2017) deals with a combined construction heuristic for a time-dependent VRP which includes service times and maximum tour duration as additional constraints.Time-dependent travel times are also considered in green VRPs (Soysal and Cimen 2017) as well as green drone routing problems (Coelho et al. 2017).In addition, time-dependent travel times can also be taken into account for point-to-point evacuation routing and scheduling (Osman and Ram 2017).
Furthermore, Huang et al. (2017) deal with path flexibility in a time-dependent VRP while to our knowledge Soysal et al. (2015) are the only ones who consider time-dependent travel times in a 2eVRP.Wang et al. (2016) evaluate the value of deterministic solutions in service network design.They show that parts of deterministic service network design solutions can be used effectively in a stochastic setting.To our knowledge the basic idea of their work has not been applied to a 2eVRP with spatial and temporal synchronization yet.
Therefore, in this paper we focus on the impact of stochastic travel times on the solution cost of a 2eVRP with synchronization and customer deliveries on both echelons.Our aim is to gain insights into the impact of these uncertainties to the deterministic solution of the problem based on fixed as well as time-dependent travel times and to evaluate methods of reoptimizing the deterministic solution based on scenario simulation of lognormally distributed travel times.

Problem description
In this paper we consider a 2eVRP with spatial and temporal synchronization, and customer deliveries on both echelons.Customers are divided into two groups.Customers located in the predefined city center have to be served by cargo bikes, whereas all other customers have to be served by vans.Vans as first echelon vehicles start their routes at the van depot at the outskirts of the city (see red large rectangle in Fig. 1), which is also the storage facility for the goods.
On the second echelon, cargo bikes start their routes at the bike depot inside the city center (see green small rectangle in Fig. 1).This bike depot is only used for parking and has no storage facilities.Thus, a cargo bike requires at least one temporally and spatially synchronized meeting with a van at one of the predefined synchronization points, the set of satellites (see blue triangles in Fig. 1), right after starting from the bike depot.Because of a cargo bike's limited capacity, reloading may also be required, once the cargo bike runs out of cargo during its tour (see cargo bike trips (dotted lines in Fig. 1) starting from the bike depot (green small rectangle), meeting a van at a satellite (blue triangle) to get the first load, then serving customers (black dots) and visiting again a satellite for loading the next part of goods).
The overall objective of this optimization problem is to minimize total cost for transporting all ordered goods to the customers either by van or by bike.This total cost consists of distance-based vehicle cost, time-based personnel cost and fixed cost for each vehicle in use.In addition to a maximum allowed route duration, synchronized meetings between vans and cargo bikes at satellites are modeled in a way that any vehicle involved in a synchronization must not wait longer at a satellite than a predefined fraction of the maximum allowed route duration.Capacity constraints for vans on the first echelon must not be exceeded during each van route, whereas capacity constraints for bikes on the second echelon become effective on any route segment between two satellite visits, as well as the last part of the bike route from a satellite to the bike depot.For the complete mathematical model the interested reader is referred to Anderluh et al. (2017).The multiple synchronizations between vans and cargo bikes in the described city distribution scheme, create a high degree of interdependence between the routes.A cargo bike can meet with a van to get the first part of goods, deliver it and then meet with another van to get the next load (see upper part in Fig. 2).
A delay of one vehicle naturally has an impact on the vehicle itself (see Van 1 in the lower part of Fig. 2) but-because of its synchronized meeting with another vehicle-it can also cause delays at that second vehicle (see Cargo bike in the lower part of Fig. 2).Furthermore, the delay can also propagate to other vehicles which have a synchronized meeting with any of the already affected vehicles (see Van 2 in the lower part of Fig. 2).In the end all vehicles displayed in Fig. 2 are affected by only one delay.Therefore, it is of interest how such a deterministic solution behaves under more realistic conditions characterized by uncertainties.In this context, we explicitly investigate the impact of stochastic travel times on the deterministic solution, because travel time deviations are constantly caused by different amounts of traffic, congestion, and/or traffic jams.
In this paper we rely on findings by Lecluyse et al. (2009) for the probability distributions of travel times and assume a lognormal distribution for the travel times of vans.For extensive information on determining appropriate travel time distributions we refer the interested reader to Gomez et al. (2016) and Aron et al. (2014).
Travel times of bikes are assumed as rather unaffected by congestion and hardly uncertain as shown for example by Jensen et al. (2010).

Solution procedure
For the deterministic problem described in Sect. 3 a two-stage GRASP with path relinking is used (see Pseudocode 1, line 1), where cargo bike routes and van routes are built sequentially.After building the cargo bike routes, improving them by local search with a 2opt-, a move-, and a swap-operator and inserting necessary satellites for synchronization purposes, the van routes are built accordingly, based on the demand and synchronization information of the already fixed cargo bike routes.Path relinking then enables finding additional promising solutions by converting one solution out of a pool of solutions into the pre-selected so-called guiding solution (Anderluh et al. 2017).
Then we apply a Monte Carlo simulation to generate N travel time scenarios for the deterministic solution based on lognormally distributed travel times.This yields different costs C s stoch for each scenario s.Based on these simulations we then calculate the average cost of stochastic travel times Cstoch = 1 N ∑ N s=1 C s stoch of the solution (see Pseudocode 1, lines 2-3) to estimate the performance of the deterministic solution under uncertain travel times.Hence, the higher the average cost of stochastic travel times Cstoch in proportion to the cost of the deterministic solution C det , the worse the performance of the deterministic solution under travel time uncertainties.Therefore, the focus of this work is to find a good way of reducing Cstoch without increasing C det significantly.

3
Impact of travel time uncertainties on the solution cost of… After the average cost of stochastic travel times Cstoch is determined, we mark in all simulated travel time scenarios infeasible paths with respect to the maximum route duration t max , which is assumed as hard constraint.Such a path is called critical path.
For each node i we evaluate the critical value i | 0 ≤ i ≤ 1 , which represents the fraction of scenarios in which node i is part of a critical path (see Pseudocode 1, line 4).Then we assign an additional buffer time b i to the given service time i of a critical node, i.e. a node i where  i > 0 (see Pseudocode 1, lines 5-7).Then this new deterministic problem is reoptimized (see Pseudocode 1, lines 10-30).
The reoptimization phase of the deterministic solution now requires additional local search operators which respect all synchronization constraints of the existing deterministic solution-that is the required spatial and temporal synchronization of vehicles at satellites and the related maximum allowed waiting time for vehicles at satellites traced through the complete solution: • A Move*-operator, which moves one customer to a position in one of the routes of the same type with respect to the synchronization constraints.• A Swap*-operator, which exchanges two customers in routes of the same type with respect to the synchronization constraints.• A MoveVanSat*-operator, which inserts available satellites in van routes with respect to the synchronization constraints.This includes also changing the satellites in the respective bike routes.• A DeleteSat*-operator, which removes unnecessary satellites from the solution.
Conducted pretests show that the most preferable order of applying these operators is: In this way, the Swap*-operator which searches a larger neighborhood than the Move*-operator, is applied second.The DeleteSat*-operator is used to remove unnecessary satellites from the solution after applying the former two operators.Then the MoveVanSat*-operator is applied to check if satellites can be shifted in the solution.
The re-evaluation of the critical paths is done regularly after each ŷ th accepted local search step (see Pseudocode 1, lines 13-28).This reassessment of critical paths is based on simulated travel time scenarios and yields new i values (see Pseudocode 1, line 21).Old i -values are halved and stored (see Pseudocode 1, lines 17-19).Then the new values for i are summed up with the old ones and new values for b i are calculated (see Pseudocode 1, lines 22-25).This is done to incorporate a kind of memory of the critical nodes during the reoptimization process.Finally, we re-evaluate the average cost of stochastic travel times of the new found solution C * stoch without buffer times and compare it with the respective average cost Cstoch calculated before (see Pseudocode 1, lines 31-34).
Table 1 provides an overview of all parameters used in this paper.

Computational results
The algorithm was coded in C/C++.All tests were performed under Linux Ubuntu 16.04 LTS on a Virtual Machine using 2 processors and 2GB memory on a host Intel(R) Core(TM) i5-3320 M CPU @ 2.60GHz 4GB RAM and running on a single thread.
Because of the different scale of the objective function values of our test instances, we focus in our tests on the relative change in average cost of stochastic travel times of the reoptimized solution C * stoch compared to the respective cost 1 3 Impact of travel time uncertainties on the solution cost of… of the original deterministic solution Cstoch (see Sect. 4 for details).Therefore, we look at to assess the performance of the reoptimized solutions.
The instances used and data sets generated during the current study are available from the corresponding author on request.
The remainder of the section starts with a detailed description of the test instances used, followed by results of pretests and parameter tuning.Then the computational results regarding the degree of synchronization of an instance are discussed.Eventually, the use of time-dependent travel times is investigated.

Test instances
For our computational tests we use three types of instances.Six adapted Solomon instances (one out of each group of Solomon instances) are used for pretests and parameter tuning.The adaption of the Solomon instances by Grangier et al. (2016) is extended by using the Solomon depot as bike depot and by adding additional vehicle information to reflect our two assumed vehicle types.
In addition, 12 synthetic test instances which cover three different layouts and two different instance sizes are generated.The three different layouts are based on the idea of the Solomon instances with randomly located customers (r), clustered customers (c), and a combination of the former two customer location schemes (rc).All customers are located in an area representing a city like Vienna.The instances   2017)).
To cover also time-dependent travel times a real-life instance of the city of Vienna is used.This instance includes 100 randomly selected pharmacies in Vienna and the van depot is the location of a pharmacy wholesaler in Vienna.18 potential satellites are assumed at appropriate locations along the outer ring street.The bike depot is assumed at an appropriate location in the city center (for further details see Anderluh et al. (2017)).We expand this instance by time-dependent travel times for the first echelon vehicles (vans) based on data described as follows.Leodolter et al. (2015) use historical floating car data (FCD) collected from 3500 taxis in the region of Vienna.The taxis are equipped with an on-board-unit providing global positioning system (GPS) data of the position of the taxis every 30 to 60 s.The GPS data is analyzed using the real-world FCD system FLEET (Toplak et al. 2011) and projected on the road network graph in order to obtain vehicle speed measurements.
The goal of the work of Leodolter et al. (2015) is to estimate the variation of travel speed in the course of a typical day.The authors present a method using a linear regression model based on static map parameters to estimate the diurnal variation of travel times for vehicles in urban areas.This enables the estimation of realistic travel times from static map coefficients without the need for further data collection.Data is extracted from the FLEET system, containing daily time series of 96 15-min intervals for 4 different speed profiles on each link of the road network of Vienna (Toplak et al. 2011).

Pretests and parameter settings
To figure out an appropriate number of scenarios for each evaluation step, we conduct pretests with an increasing number of scenarios.In these tests a sample size of N = 60 yields already an acceptably low deviation in the simulated travel time scenarios, which we express in the additional average cost of stochastic travel times ( Cstoch − C det ) in cost units.Furthermore, this sample size also guarantees good results in a reasonable amount of computational time (see Fig. 3).
Three different methods for allocating buffer times are tested.The total amount of buffer time B, which is a fraction of the maximum route duration t max , is assigned to each node i: (1) Proportionally based on i .
(2) As equal share to any node i with  i > 0.5.
(3) As a combination of the former two ways. 2 / 3 of B are based on (1), the rest is based on (2).
The reassignment of B is done either only once before starting the reoptimization or after every ŷ th accepted improvement step in the local search of the reoptimization with ŷ ∈ (5000, 1000, 500, 100, 10).

Impact of travel time uncertainties on the solution cost of…
The different parameter settings were tested with all pretest instances and an average of five runs was considered for each instance.Table 2 shows the ranking of all combinations with respect to change in average cost of stochastic travel times, where 1.000 as the maximum value represents the best parameter setting.
The tests show that using type (1) for the assignment of buffer time B and reoptimizing the critical values i after every 1000th acceptance of a local search step yield the best results.
Nevertheless, the impact of the assignment type of B is rather small.Furthermore, the assessment of different reoptimization modes depicted in Table 2 show that calculating the critical values i after every 1000th acceptance step in the local search yields the best results for each assignment type of B. Doing the reoptimization only once at the beginning performs worst for assignment types (1) and (3), while assignment type (2) gives the worst result while reoptimizing after every 100th acceptance step.In all cases reoptimizing after every 5000th, 500th, 100th or 10th acceptance step in the local search does not improve the results, although no clear pattern can be distinguished.In addition, computational times for all settings show that ŷ = 1000 yields good results in a small amount of com- putational time (see Fig. 4).
To determine an appropriate value for B, which we assume as fixed value, tests for different percentages of t max (5, 10, 15, 20 and 25) over all test instances have

Impact of the degree of synchronization
The level of interdependence in a solution depends on several drivers like the number of customers in the city center, the demand of these customers, and the specific location of these customers.These factors impact the total load of cargo bikes, and such the required number of synchronized meetings with vans.One of the main influencing factors is the degree of synchronization in a problem instance.This factor is strongly correlated with Q b z -the capacity of the second echelon vehicles (cargo bikes) which reflects a certain synchronization level z.We performed our computational tests with a capacity varying from a minimum required capacity represents the highest demand of all customers supplied by a second echelon vehicle in the instance, as every customer has to be served by exactly one visit.Q b max is assumed as half of the first echelon vehicle's capacity.Between these two bounds the capacity A major insight from our computational results is that the amount of the decrease in average cost of stochastic travel times clearly depends on the degree of synchronization which is reflected by Q b z .Lower capacity and hence the necessity of more synchronizations (see dashed line in Fig. 6) for a feasible solution of the problem is related to a stronger decrease in average cost of stochastic travel times as depicted by the solid line in Fig. 6.A low level of bike capacity also implies a higher average number of successors (see dash-dotted line in Fig. 6), which reflects the average number of nodes in an instance that follows directly (in the same route) or indirectly (in a synchronized route) each node of the instance.These findings hold for all our test instances.
Figure 6 also shows that for low bike capacity levels 1 and 2-which reflect a high degree of synchronization-the average cost of stochastic travel times can be reduced by the reoptimization step by around 4% on average over all test instances.From bike capacity levels 3-6 this reduction is around 2% and for bike capacity levels 7-9 this reduction stabilizes around 1%.
Therefore, our results show that the more interdependent the instance is, the more the performance of the deterministic solution under travel time uncertainties can be improved by reoptimization.
On the other hand the increase in deterministic cost by reoptimization is rather constant at around 1% on average for all bike capacity levels (see dotted line in Fig. 6).
Detailed results for each instance and bike levels are depicted in Tables 3, 4 and 5 in the "Appendix".These results show that, especially for the clustered test instances (n*_c*, and 'vienna'), the average cost of stochastic travel times can be reduced by up to around 9%, whereas instances with randomly distributed customers (n*_r*) in most cases perform worse in our setting.
In contrast to that the size of the instance does not show a significant impact in our tests.

Use of time-dependent Travel Times
For the Vienna test instance we test also the influence of time-dependent travel times in the reoptimization process.Therefore, we calculate solutions for starting times of first echelon vehicles from 7:00 to 13:00 in 1-h steps.We compare these results with the solution calculated based on free flow times.Additionally, we use the idea of linear regression (Leodolter et al. 2015) based on time-dependent travel times in the provided 1-h slots from 7:00 am to 7:00 pm on typical working days.Comparing free flow times for all arcs in the Vienna instance with the average time based on the time-dependent times over the working day is depicted in Fig. 7. Therefore, recalculating travel times based on free flow times with the above mentioned linear regression seems to be a viable estimate.
Travel times for the second echelon vehicles (cargo bikes) are assumed as steady over the working day as bikes are usually not delayed by accidents or traffic jams.This assumption is also supported by Jensen et al. (2010), who show in their survey that cyclists travel at nearly the same speed during the day.Furthermore, Fig. 8 shows the comparable behavior when using time-dependent travel times as well as travel times gained by applying linear regression.On the other hand the use of free flow travel times causes proportionally higher additional cost of stochastic travel times.
Based on these results we conclude that, at least for our realistic test instance of the city of Vienna, using adapted free flow times performs similarly to using timedependent travel times.This is of importance because deriving time-dependent travel times based on floating car data for the complete road network requires significantly more effort than determining free flow travel times from publicly available sources as for example openstreetmap.org,surveying time-dependent data only for some representative arcs and deriving the linear regression based on these observations.

Conclusion
In this paper we analyze the performance of an innovative city distribution scheme for a two-echelon vehicle routing problem with temporal and spatial synchronization between vans and cargo bikes in a realistic setting.Travel time uncertainties are taken into account.
We present a reoptimization method for the deterministic solution of the problem based on stochastic information to improve the quality of the overall solution.Computational results show that our simple iterative reoptimization of the deterministic solution based on information gained by evaluating simulated travel time scenarios can improve the behavior of the solution under real-life circumstances significantly.The evaluation of all test instances shows that the performance of our method depends on the degree of synchronization-that is the number of synchronizations between firstechelon and second-echelon vehicles-in the solution.The higher the degree of synchronization the more can be gained by applying the reoptimization step.We expect this to generalize to other vehicle routing problems containing synchronization constraints.
Using time-dependent travel times has turned out to yield no better results than transformed free flow travel times, where the transformation is based on a simple linear regression.These results indicate the good performance of the latter approach, which is of importance especially for decision makers, because such travel times are easier to obtain than time-dependent ones.
Additional stochastic information concerning for example uncertain service times as well as dynamic requests are one direction for future research.
To cover the main advantages of the proposed cargo bikes-especially their lack of local emissions and of noise-future research should also take into account the evaluation of emissions, like greenhouse gases, particulate matters and nitrogen oxides, as well as the evaluation of other external factors, like noise, side-effects of congestion and accidents within the optimization process.

3
Impact of travel time uncertainties on the solution cost of…

Fig. 1
Fig. 1 City distribution scheme with vans and cargo bikes

Fig. 2
Fig. 2 Interdependence between routes because of synchronization (upper part) and potential impact of a delay (lower part) Capacity of a second echelon vehicle at capacity level z Q v Capacity of a first echelon vehicle C det Cost of original solution based on deterministic travel times C s stoch Cost of original solution based on travel times generated in scenario s Cstoch Average cost of the original solution evaluated averaged over the N generated scenarios C * det Cost of reoptimized solution based on deterministic travel times C * stoch Average cost of the reoptimized solution evaluated averaged over the N generated scenarios ŷ Number of accepted LS-steps after which the recalculation of i is applied include 100/125 (n100, n125) customers and 10 satellites.The satellites are located around the assumed city center approximately in the middle of the instance area.The values for the vehicles used reflect realistic values for speed, capacity as well as costs of cargo bikes and 3.5-ton vans.Service time and demand of customers are random numbers in an appropriate range (for further details see Anderluh et al. (

Fig. 3
Fig. 3 Impact of sample size of scenario simulation on deviation in additional average cost of stochastic travel times (in cost units) and average cpu-time (in s)Table 2 Ranking of buffer time assignment type and reoptimization mode with respect to change in average cost of stochastic travel times

Fig. 4
Fig. 4 Average computational time for recalculation of critical values after ŷ ∈ (10, 100, 500, 5000) local search steps and only once at the beginning

Fig. 5
Fig. 5 Trade-off between increase in deterministic cost and decrease in average cost of stochastic travel times depending on the amount of B as fraction of t max

Fig. 7
Fig. 7 Linear regression of average time-dependent travel times to free flow times

Fig. 8
Fig. 8 Notion of additional average cost of stochastic travel times on cost of deterministic original and reoptimized solution calculated for different travel time basis

Table 1
Overview of parameters and variables

Table 3
Detailed results for instance 'Vienna'

Table 5
Detailed results for instances 'n125*' InstQ z Impact of travel time uncertainties on the solution cost of… InstQ z