Next Article in Journal
Determining Optimal Parameters of Regular Microrelief Formed on the End Surfaces of Rotary Bodies
Next Article in Special Issue
Modeling and Optimization in Resource Sharing Systems: Application to Bike-Sharing with Unequal Demands
Previous Article in Journal
Non-Overlapping LZ77 Factorization and LZ78 Substring Compression Queries with Suffix Trees
Previous Article in Special Issue
Simulation-Optimization Approach for Multi-Period Facility Location Problems with Forecasted and Random Demands in a Last-Mile Logistics Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combining Heuristics with Simulation and Fuzzy Logic to Solve a Flexible-Size Location Routing Problem under Uncertainty

1
IN3–Computer Science Department, Universitat Oberta de Catalunya, 08018 Barcelona, Spain
2
School of Engineering, Universidad de La Sabana, Chia 250001, Colombia
3
Department of Data Analytics & Business Intelligence, Euncet Business School, 08221 Terrassa, Spain
4
International School of Economics and Administrative Sciences, Universidad de La Sabana, Chia 250001, Colombia
*
Author to whom correspondence should be addressed.
Submission received: 15 December 2020 / Revised: 21 January 2021 / Accepted: 26 January 2021 / Published: 30 January 2021
(This article belongs to the Special Issue Simulation-Optimization in Logistics, Transportation, and SCM)

Abstract

:
The location routing problem integrates both a facility location and a vehicle routing problem. Each of these problems are NP-hard in nature, which justifies the use of heuristic-based algorithms when dealing with large-scale instances that need to be solved in reasonable computing times. This paper discusses a realistic variant of the problem that considers facilities of different sizes and two types of uncertainty conditions. In particular, we assume that some customers’ demands are stochastic, while others follow a fuzzy pattern. An iterated local search metaheuristic is integrated with simulation and fuzzy logic to solve the aforementioned problem, and a series of computational experiments are run to illustrate the potential of the proposed algorithm.

1. Introduction

When designing and managing supply chains, one of the most relevant problems is the simultaneous location of distribution facilities and the routing of vehicles to deliver products to a set of geographically dispersed customers. The former is considered a strategic decision, while the latter is operational. This problem is known in the scientific literature as the location routing problem (LRP). The LRP addresses these two types of decisions in an integrated manner. From the formal view of the operational research community, the LRP is known to be NP-hard, since it can be reduced to either the facility location problem (FLP), the vehicle routing problem (VRP) or the multidepot VRP, which are all known to be NP-hard. This computational complexity means that optimal solutions are really difficult to obtain in a reasonable computational time. Thus, heuristic approaches are required to solve medium- and large-sized instances. Due to its complexity, some of the first studies tackled the problem by splitting it into the corresponding subproblems [1,2]. Nevertheless, this approach might lead to suboptimal solutions.
Due to the increase in computational power and the development of fast heuristic approaches, the LRP has been studied in an integrated way, which clearly has improved the obtained results [3]. One of the most studied versions of the LRP is the capacitated LRP, in which both depot and vehicle capacity constraints must be satisfied (the acronym LRP will henceforth refer to this version). However, all previous works consider the depot capacity as a fixed value for each location. This could not be a suitable approach when dealing with realistic problems, since it is usual that decision-makers can select the size of a facility from a discrete set of known available sizes, or even freely. For real-world problems, this set is usually associated with investment activities, such as building facilities [4], purchasing equipment [5] or qualifying workforce [6]. From an academic point of view, despite the increasing number of published works on the LRP, the consideration of flexible sizes for facilities has been rarely addressed in the literature. Nevertheless, real-life examples from both LRP [4,7,8] and non-LRP [5,6] contexts show the relevance of considering a variety of facility sizes to select from.
Traditional LRP approaches consider that parameters are deterministic or crisp, i.e., they assume that inputs are known in advance. This assumption is far from reality in many applications, such as waste collection, humanitarian logistics and urban freight distribution, where uncertainty is a key factor to consider. Despite this, the literature on the LRP addressing uncertain parameters is still scarce. In order to overcome this problem, articles employing stochastic approaches can be found in the literature. Customers’ demand is one of the most addressed stochastic parameters [9,10,11,12,13]. Other parameters might also be considered as stochastic, such as transportation costs and travel speeds [14] or logistic costs and travel distance [15]. In general, many articles addressing stochasticity in routing problems hybridize simulation models with heuristic or metaheuristic algorithms to tackle efficiently both uncertainty and NP-hardness. In many real-life situations, however, it might not be possible to accurately model all uncertainty sources as stochastic variables following a probability distribution. This might be the case, for instance, when the volume of observations is low or the available data does not have enough quality [16]. Hence, uncertainty in the LRP has also been tackled through the use of fuzzy sets. Parameters such as customers’ demands [17,18,19,20], travel times [21,22] or time windows [23] have been modeled as fuzzy in several studies. Notice that, whenever possible, modeling uncertainty as stochastic variables might allow a deeper statistical analysis of the results.
To the best of our knowledge, there are no works in the literature simultaneously addressing stochastic and fuzzy approaches to model demand uncertainty in a flexible-size LRP. This is a realistic scenario, since many companies might have historical data on trustworthy customers and not enough data on new or unreliable ones. Hence, the main contributions of this paper are two-fold: on the one hand, a new variant of the location routing problem is studied, where facility sizing decisions and hybrid fuzzy-stochastic demands are simultaneously considered. On the other hand, this paper proposes a competitive solution approach based on the hybridization of a metaheuristic algorithm with both simulation and fuzzy logic, i.e., a so-called fuzzy simheuristic, to solve the aforementioned problem. Indeed, simheuristics have been traditionally proposed to deal with stochastic issues in hard combinatorial optimization problems [24]. However, their hybridization with fuzzy logic has been rarely studied.
The remainder of this paper is organized as follows: Section 2 describes previous works on the topic. Section 3 presents a description of our addressed problem. Section 4 explains the fuzzy simheuristic approach used to solve the problem. Section 5 describes a series of computational experiments. Section 6 analyzes our obtained results. Finally, Section 7 draws some conclusions and future research perspectives.

2. Literature Review

This section presents a summary of the published manuscripts on the main topics addressed by this work. Thus, Section 2.1 outlines works related to the location routing problem in both its deterministic version and the variant including uncertain parameters. Additionally, Section 2.2 summarizes the main contributions on the field of simheuristics and fuzzy logic as methodologies to handle uncertain parameters in routing problems.

2.1. The Location Routing Problem

Perhaps the first work related to the location routing problem is the one by Maranzana [25], who analyze the influence of transportation costs on location decisions. Moreover, Salhi and Rand [1] quantify for the first time the benefits of considering routing decisions when locating facilities. They also state that solving each subproblem (location and routing) independently does not provide optimal solutions. Multiple variants of the LRP have been proposed over time. These variants depend on the characteristics of depots (capacitated or not), vehicles (capacitated or not, homogeneous or heterogeneous fleet), costs (symmetric or asymmetric) or the consideration of uncertain parameters. All capacitated variants addressed by these authors assume that depot sizes are fixed and cannot be changed.
Considering the limited computational power available at that time, the initial works on the LRP firstly solved the underlying location problem and used the obtained solution as a starting point to handle the corresponding routing problem. However, as the computational power has notably increased in recent years, the newest approaches deal with the LRP in an integrated manner [26,27]. Among the recently published works on the deterministic LRP, Escobar et al. [28] propose a granular tabu search within a VNS framework to speed up computational times without decreasing the solutions quality. A biased-randomization-based metaheuristic of two phases is developed by Quintero-Araujo et al. [27] to solve the capacitated version of the problem. Ferdi and Layeb [29] propose a GRASP with a novel technique used to create clusters around the open depots. Traditional applications of the LRP include horizontal cooperation [30], electric vehicle routing problems [31,32,33], city logistics [34], humanitarian logistics [35] or supply chain network design [36]. Moreover, most recent applications are related to environmental issues [37], cold supply chains [38] or waste management [39].
When dealing with uncertainty, most works have focused on the use of stochastic modeling. One of the utilized approaches has been the hybridization of simulation techniques with metaheuristics. For instance, Quintero-Araujo et al. [9] propose a simheuristic to solve an LRP with stochastic demands, by hybridizing Monte Carlo simulation with an iterated local search metaheuristic. A similar approach is employed by Tordecilla et al. [13], who address an LRP where the sizes of facilities to locate are also a variable to consider. Rabbani et al. [10] also propose a simheuristic approach that combines a nondominated sorting genetic algorithm-II (NSGA-II) and Monte Carlo simulation. They tackle a multiobjective multiperiod LRP in the context of the hazardous waste management industry. Both generated waste and number of people at risk are stochastic. Inventory decisions are also taken into account. Sun et al. [11] address a real-world case from an express delivery company in Shanghai. These authors tackle an LRP in which demand is stochastic and can be split for self-pickup. Then a simulation-based optimization model is proposed and two heuristics results are compared.
Other parameters are also considered to be uncertain. For instance, Herazo-Padilla et al. [14] hybridize an ant colony optimization metaheuristic with discrete-event simulation to solve an LRP in which both transportation costs and vehicle travel speeds are considered stochastic. Authors demonstrate that their proposed approach is not only efficient but is able to find statistical interactions among the different parameters. Zhang et al. [15] present an approach that hybridizes a genetic algorithm with simulation to solve a sustainable multiobjective LRP in the context of emergency logistics. The authors consider the travel distance, the demand and the cost of opening a depot as uncertain variables. Additionally, the emergence of new technologies introduces new challenges. This is the case of Zhang et al. [12], who address the problem of locating battery swap stations and routing electric vehicles with stochastic demands. This problem is solved using a hybrid approach that combines a variable neighborhood search with a binary particle swarm optimization algorithm. The problem’s complexity increases when considering the low autonomy of this type of vehicles, since route failures can frequently be present when demands are not known in an accurate manner.
Uncertainty in the LRP has been studied using either stochastic or fuzzy parameters. Table 1 shows an overview of works addressing this topic, which includes: (i) whether the uncertainty is addressed stochastically or in a fuzzy fashion; (ii) the considered uncertain parameter; (iii) the mathematical modeling approach; (iv) the approach used to solve the problem; and (v) the objective function. Analyzed works show a clear preference for considering an uncertain demand, as well as for using fuzzy chance constrained models. Given both the considered uncertainty and the combinatorial nature of the LRP, most works employ a hybrid approach combining simulation with a metaheuristic algorithm. Finally, cost minimization is the prevalent objective, although a few works also consider the minimization of risk or the minimization of the additional travel distance due to route failures. Regarding works on fuzzy parameters, Zhang et al. [17] propose a hybrid particle swarm optimization (PSO) algorithm to solve a capacitated LRP with fuzzy triangular demands (CLRP-FD). The hybrid PSO algorithm is composed of three phases including a local search method and stochastic simulation. In addition, the authors propose a chance-constrained programming model for the CLRP-FD. Zarandi et al. [21] consider a multidepot LRP with fixed depot capacity and fuzzy travel times. Mehrjerdi and Nadizadeh [18] present a fuzzy chance constrained programming model where demands are modeled as fuzzy numbers. A four-phase method called “greedy clustering” is proposed, in which both an ant colony system metaheuristic and stochastic simulation are included. Fazayeli et al. [19] propose an LRP with time windows and fuzzy demands as the delivery part of a multimodal transport network. The mixed integer mathematical fuzzy model is coded and solved using GAMS and compared to the results provided by a genetic algorithm. Nadizadeh and Kafash [20] analyze a LRP with simultaneous pick-up and delivery in the context of reverse logistics. Both types of demands (pick-up and deliveries) are fuzzy variables. A fuzzy chance constrained programming model is proposed to represent the problem, and a greedy clustering method is used to solve it.
The analyzed works show that uncertainty in the LRP has been addressed either by using stochastic or fuzzy demands but never considering both types of uncertainty at the same time—e.g., that some customers’ demands are modeled as stochastic variables while others are modeled as fuzzy values. In addition, to the best of our knowledge, there are no previous studies on the LRP with facility sizing decisions and hybrid fuzzy-stochastic demands. Only Tordecilla et al. [13] have studied a similar LRP variant, although considering all customers’ demands as stochastic. Thus, our work aims to fulfill the existing gap by considering a flexible-size LRP and two different types of uncertain parameters: stochastic and fuzzy demands.

2.2. Simheuristics and Fuzzy Logic for Vehicle Routing Problems under Uncertainty

When dealing with combinatorial optimization problems subject to uncertain parameters, one of the most recommended approaches is the combination of simulation (to handle stochasticity) with heuristic-based methods (to deal with the optimization part of the problem) [43,44]. In that sense, a simheuristic approach is a relatively new and efficient technique to tackle combinatorial optimization problems under uncertainty [24,45]. In general, a simheuristic algorithm works as follows: (i) given a stochastic problem, the random variables are transformed into their deterministic counterpart by using expected values; (ii) an approximated framework (heuristic or metaheuristic) is used to generate high-quality solutions for the transformed deterministic instance that can also be “promising” solutions for the stochastic version of the problem; (iii) these promising solutions are sent to a simulation engine in order to estimate its quality in a stochastic environment. The simulation engine, in addition, provides feedback to better guide the search used by the approximated procedure; and (iv) an improved estimation of the quality of the solutions is obtained for a subset of “elite” solutions using a longer simulation process. Different simheuristic algorithms have been presented in the literature to solve routing problems.
Stochastic demands in vehicle routing problems are addressed by Quintero-Araujo et al. [46] and Gruler et al. [47]. Moreover, stochastic demands are also studied in arc routing problems [48]. Stochastic versions of the inventory routing problem can be found in Gruler et al. [49]. Real applications like the waste collection problem with stochastic demands are analyzed in Gruler et al. [50]. Intermodal routing problems have also considered other stochastic parameters, such as capacity [51] or travel times [52,53]. Additionally, the need of using fuzzy logic in vehicle routing problems arises when there are some vague or uncertain parameters. The literature presents various works in which fuzzy logic is added, for instance, to model uncertain demands [54,55,56,57], travel times [58,59], capacity [57,59], and service times [60]. Additional aspects are also considered by these works, such as time windows [57,61,62], environmental aspects [59], multiple objectives [59], intermodal transportation [57,59] and an open VRP [63]. Additional applications of metaheuristics combined either with Monte Carlo simulation or fuzzy logic can be found in several fields, such as scheduling [64,65], controller optimization [66,67] parameter estimation [68], finance [69], facility location [70], etc.

3. Problem Description

The location routing problem is a well-known problem in which three main decisions must be made: (i) locating one or more facilities; (ii) allocating customers to open facilities without exceeding their capacity; and (iii) designing a number of routes whose aggregated customers’ demand does not exceed a vehicle capacity. Each route must start and finish at the same facility. Furthermore, we consider a location routing problem with facility sizing decisions, where the size of each open facility is also a variable to decide on. Furthermore, we also consider both stochastic and fuzzy demands. Hence, a percentage of the vehicles’ capacity is reserved as a safety stock (SS), in case the demand is higher than expected. Therefore, the main decision variables in this problem are related to the number of facilities to open, the facilities’ size and location, which customers must be allocated to each open facility, how many vehicles must be used and how to design the associated routes. This problem is NP-hard since it contains, as special cases, the capacitated vehicle routing problem, the multidepot VRP and the facility location problem, all of them known to be computationally hard. Figure 1 provides an example of a complete LRP solution. Facilities are represented by diamonds and customers by circles. Black (solid) diamonds are the open facilities, while noncolored diamonds correspond to nonopen facilities. For each open depot, a set of routes starting and finishing at the corresponding depot location is designed to serve all customers’ demands. Each route is assigned a single vehicle.
Formally speaking, the LRP can be defined on a complete, weighted, and undirected graph G ( V , E , C ) , in which V = J I is the set of nodes (comprising the subset J of potential facility locations and subset I of customers), E is the set of edges, and C is the cost matrix of traversing each edge. Delivery routes are performed by a set K of unlimited homogeneous vehicles with limited capacity. This problem also assumes that all vehicles are shared by all facilities (i.e., no depot has a specific fleet) and each edge e E satisfies the triangle inequality. The customers’ demands are uncertain and are modeled using stochastic values for a subset of customers I 1 , and fuzzy values for a subset of customers I 2 , such that I 1 I 2 = I . The variant of the LRP considered in this paper is the one in which a decision must be made about the size of the facilities to open. Hence, a set L of alternative sizes for each facility and associated fixed and variable opening costs are provided as inputs. Depots might have equal or different capacities. Each customer node must be served by exactly one vehicle that starts and finishes its route in the facility to which it has been allocated (i.e., split deliveries are not allowed). The following notation is used to describe our problem:
Parameters
s l = Available size of type l L
D i = Uncertain demand of customer i I
f j = Fixed opening cost of depot j J
o j l = Variable opening cost of depot j J with size of type l L
c e = Cost of traversing arc e E
q = Capacity of each vehicle
% S S = Safety stock percentage
Decision variables
y j l = Binary variable that indicates whether the depot j J is open with size l L or not.
x i j = Binary variable that indicates whether customer i I is assigned to the depot j J or not
w e k = Binary variable that indicates whether arc e E is used in the route performed by vehicle k K or not
The objective is to minimize the total cost ( T C ), which includes opening facilities costs ( O C ), routing costs ( R C ), and failure costs ( F C ), i.e., T C = O C + R C + F C . These parts are defined in Equations (1)–(3).
O C = j J l L ( f j + o j l ) y j l
R C = e E k K c e w e k
F C = min { c r e a c , c p r e v }
F C represents the cost incurred whenever the actual demand of a route is greater than the vehicle capacity, where c r e a c and c p r e v depend on the corrective action considered, namely:
  • A reactive strategy with a cost c r e a c , in which a vehicle must perform a round-trip to its assigned facility for a replenishment if the actual current-customer demand is higher than the vehicle’s current load.
  • A preventive strategy with a cost c p r e v , in which a vehicle must perform a detour to the facility before visiting the next customer. The decision about performing this detour depends on the type of demand of the next customer. If the demand is stochastic, the detour is carried out whenever the expected demand of the next customer is higher than the current capacity of the vehicle. Alternatively, if the demand is fuzzy, this decision depends on the comparison between the fuzzy values of both the demand of the next customer and the current capacity.
Let S V be a subset of nodes, δ + ( S ) the set of edges leaving S, δ ( S ) the set of edges entering S, and A ( S ) the set of edges with both ends in S. Hence, the location routing problem with facility sizing decisions and uncertain demands can be modeled as the following integer program:
M i n i m i z e T C
                              subject to:
k K e δ ( i ) w e k = 1 i I
i I e δ ( i ) D i w e k ( 1 % S S ) q k K
e δ + ( n ) w e k = e δ ( n ) w e k , k K , n V
e δ + ( J ) w e k 1 k K
e A ( S ) w e k | S | 1 S I , k K
e δ + ( j ) w e k + e δ ( i ) w e k 1 + x i j i I , j J , k K
j J x i j = 1 i I
i I D i x i j l L s l y j l j J
l L y j l 1 j J
y j l , x i j , w e k { 0 , 1 }
The objective function (4) minimizes the total cost. Constraint (5) ensures that each customer is served by a single route and a single vehicle. Constraint (6) guarantees that the total demand served by a vehicle in a route does not exceed its capacity. This limit is reduced by a safety stock, which is a percentage of the vehicle capacity reserved to respond more effectively to the uncertain demand. Constraint (7) guarantees the continuity of each route. Constraint (8) ensures the return of each vehicle to its starting depot. Constraint (9) guarantees the subtour elimination. Constraint (10) ensures that a customer is served by a route departing from an open depot only if this customer is allocated to this depot. Constraint (11) guarantees that a customer is assigned to only one depot. Constraint (12) ensures that the total demand served from a depot does not exceed its assigned size. Constraint (13) guarantees that only one size is assigned to an open depot. Finally, Constraint (14) determines that all decision variables are binary.

4. Solution Approach

Since the problem described in Section 3 is known for being NP-hard, the formulated mathematical model is not employed to find an optimal solution but just to provide a better understanding of the problem details Hence, we propose a fuzzy simheuristic approach [24] for minimizing the expected total cost. Traditionally, simheuristics have been used to solve optimization problems with stochastic components, such as arc routing problems with stochastic demands [48], stochastic waste collection problems [50] or team orienteering problems with stochastic travel times [71]. We have extended the simheuristic framework by including fuzzy components in order to deal with combinatorial optimization problems with uncertainty components of both stochastic and nonstochastic nature. In particular, our methodology combines an iterated local search (ILS) metaheuristic with Monte Carlo simulation and fuzzy inference systems (FIS) to deal with stochastic and fuzzy variables, respectively. As discussed in Ferone et al. [72], several metaheuristic frameworks offer a well-balanced combination of efficiency and relative simplicity and can be easily extended to a fuzzy simheuristic. In general, our approach is composed of three stages. During the first stage, a set of promising LRP solutions are generated using a constructive heuristic, which employs biased-randomization techniques [73]. In the second stage, the ILS metaheuristic tries to improve each of these promising solutions by iteratively exploring the search space and conducting a short number of simulations. Finally, in the third stage, a refinement procedure using a larger number of simulation runs is applied to these elite solutions, which allows one to obtain a more accurate estimation of the expected total cost.
Algorithm 1 outlines the main components of Stage 1. It generates quickly a ranked list of “promising” LRP solutions. The main input parameters of this heuristic are: the list of customers with both their demand and location in Cartesian coordinates, the list of facilities including their opening costs and the vehicle capacity. The algorithm procedure is as follows: initially, the minimum and maximum ( n b D e p o t s 0 and m a x N b D e p o t s , respectively) numbers of facilities required to serve the total demand are computed. Both bounds are calculated by dividing the total demand by the maximum available facility size, and the minimum available facility size, respectively, and they are rounded up to the next integer number. Then we run our algorithm for each number of facilities between n b D e p o t s 0 and m a x N b D e p o t s (line 3). Later, for each iteration of the line 4 loop, a new set of random locations are generated (line 5). This is stored in u s e d O p e n D e p o t s to avoid repeating. Next, if the available capacity of facilities in o p e n D e p o t s is enough to satisfy customers demand, customers’ allocation and routing procedures are carried out; otherwise, o p e n D e p o t s is rejected. The customers’ allocation procedure is performed by producing a new m a p (line 9) where each facility has a list of all customers sorted by savings. These savings represent the benefit of allocating each customer to the current depot instead to the best alternative facility. Then a facility in o p e n D e p o t s is selected randomly, and a biased-randomized procedure is used to allocate a customer of the list to the current depot. This procedure ends when all customers have been allocated. In the step in line 10 a VRP is solved for each subset facility-customers in the map. Finally, a feasible LRP solution is yielded and stored in the pool of solutions p o o l S o l . The algorithm ends returning a top list of complete LRP solutions, assessed in terms of opening and routing costs.
Algorithm 1 Constructive heuristic ( c u s t , d e p o t s , v e h C a p , β , i t e r m a x )
1:
u s e d O p e n D e p o t s ←∅
2:
n b D e p o t s 0 , m a x N b D e p o t s c o m p u t e D e p o t s B o u n d ( d e p o t s )
3:
for n b D e p o t s n b D e p o t s 0 to m a x N b D e p o t s do
4:
      for i t e r ← 1 to i t e r m a x do
5:
             o p e n D e p o t s d e p o t s T o O p e n ( n b D e p o t s )
6:
            if o p e n D e p o t s u s e d O p e n D e p o t s then
7:
                if c a p a c i t y ( o p e n D e p o t s ) d e m a n d ( c u s t ) then
8:
                     u s e d O p e n D e p o t s a d d ( u s e d O p e n D e p o t s , o p e n D e p o t s )
9:
                     m a p a l l o c a t e C u s t o m e r s ( o p e n D e p o t s , c u s t )
10:
                   l r p S o l C W S ( m a p , β , v e h C a p )
11:
                   p o o l S o l a d d ( p o o l S o l , l r p S o l )
12:
              end if
13:
           end if
14:
       end for
15:
end for
16:
return s o r t i n g B y C o s t ( p o o l S o l )
Algorithm 2 outlines Stages 2 and 3. During the second stage, each “promising” map generated by the constructive heuristic is processed by the simulation and the fuzzy components to estimate its safety stock (line 4). This procedure is carried out by performing a low number of runs, where a new value is assigned to each random or fuzzy element based on its probability distribution or fuzzy function, respectively. We use Monte Carlo simulation in order to estimate the stochastic variables, whilst a fuzzy inference system is used to estimate the fuzzy variables. Then, the objective function and the constraints are evaluated under the random/fuzzy generated values to compute the expected cost of each promising map. Next, the ILS metaheuristic tries to improve the set of “promising” maps by iteratively exploring the search space and conducting a second process of fuzzy/simulation runs. We start the process by perturbing the current base solution baseSol (line 8). In this phase we use two different strategies. In the first one, the algorithm randomly selects a set of customers and tries to reassign them in a random way to another facility without violating its capacity. Regarding the second strategy, the algorithm randomly exchanges the allocation of a percentage of customers among facilities. This process is dependent on the value of k, which represents the degree of exchange to be applied. This value is updated in each iteration between K m i n and K m a x , i.e., it is reset to K m i n whenever a new solution newSol outperforms the baseSol, and it is increased whenever the algorithm fails to improve the current solution until a maximum value K m a x . The strategy to be used in each iteration of the algorithm is randomly selected.
Algorithm 2 ILS-based Fuzzy Simheuristic ( i n p u t s , α , β , λ , I n c , T 0 , K m i n , K m a x , I 0 , t m a x )
1:
i n i t S o l ← genInitSol(inputs, α , β )
2:
b a s e S o l i n i t S o l
3:
b e s t S o l b a s e S o l
4:
fastSimulation(baseSol)% Fuzzy and Monte Carlo Simulation
5:
T ← T 0
6:
while (time ≤ t m a x ) do % ILS stage
7:
    k ← K m i n
8:
    perturbationSol ← perturbation(baseSol, k, α , β )
9:
    newSol ← localSearch(perturbationSol)
10:
   if (detCost(newSol) < detCost(baseSol)) then
11:
        fastSimulation(newSol) % Fuzzy and Monte Carlo simulation
12:
        if (expCost(newSol) < expCost(baseSol)) then
13:
             baseSol ← newSol
14:
             if (expCost(newSol) < expCost(bestSol)) then
15:
                 bestSol ← newSol
16:
                 insert(poolBestSol,bestSol)
17:
             end if
18:
             k ← K m i n
19:
          end if
20:
      else% SA-based acceptance criterion
21:
          temperature ← updateTemperature(detCost(newSol), detCost(baseSol), T)
22:
          if ( U (0,1) ≤ temperature) then
23:
               baseSol ← newSol
24:
               k ← K m i n
25:
          else
26:
               k ← min(k * Inc, K m a x )
27:
          end if
28:
      end if
29:
      T ← λ T
30:
end while
31:
for (sol ∈ poolBestSol) do % Refinement stage - Fuzzy and Monte Carlo simulation
32:
      longSimulation(sol)
33:
      if (expCost(sol) < expCost(bestSol)) then
34:
          bestSol ← sol
35:
      end if
36:
end for
37:
return bestSol
Afterwards, the algorithm starts a local search around the perturbed solution in order to improve it (line 9). This stage consists in a two-opt inter-route operator, which interchanges two chains of randomly selected customers between different facilities. A newSol is returned whenever no more improvements are achieved. Later, whenever the deterministic cost of the baseSol is improved (line 10), the newSol is processed by the simulation and the fuzzy components to deal with the uncertainty of the proposed problem, using a low number of runs to compute the expected cost of the solution (line 11). Notice that this procedure does not only provide estimated values to the expected cost associated with the solutions generated by our approach, but it also reports feedback to the metaheuristic search process. If the newSol is also able to improve the expected cost of the baseSol (line 12), the latter is updated. In the same way, if the expected cost of the newSol improves the cost of the best solution (bestSol) found so far (line 14), the latter is updated and added to the pool of elite solutions (line 16). This pool contains the best stochastic/fuzzy solutions found so far. The number of solutions in this pool is a known parameter that depends on the available computational time. Moreover, by limiting the size of this pool we ensure that we only keep track of the top solutions as the algorithm evolves. In order to further diversify the search, the algorithm might occasionally accept nonimproving solutions following an acceptance criterion (lines 20–28). Specifically, we have used a simulated-annealing acceptance criterion, which contains a decaying probability that is regulated by a dynamic temperature parameter (T).
Finally, a refinement procedure using a larger number of simulation runs is executed in the third stage for each elite solution (lines 31–36). Hence, a more accurate summary of output variables can be obtained. As before, both probability distributions and fuzzy functions are employed in this simulation, depending on whether the element has a stochastic or fuzzy nature. Finally, the "best" solution (or pull of best alternative solutions) is returned, considering that the decision maker might be not only interested in the average value associated with a solution but also in its variability level. Particularly, the main output variables in our experiments are: the opening and routing costs, the cost incurred whenever a route fails and the safety stock.

5. Computational Experiments

Multiple sets of instances are found in the literature to test the algorithms designed to solve the LRP [74,75,76]. Nevertheless, these sets do not consider characteristics such as parameters uncertainty and flexible facility sizes, i.e., instances must be adapted to our problem’s features. Therefore, we use Akca’s [74] instances and introduce the following modifications:
  • Traditional LRP instances consider that a single fixed size is available to assign to open depots. We extend this unit set to five alternative sizes, so that our algorithm selects one of them for each open depot. If s j is the size proposed by the original instance for each potential depot j J , and L is the set of available sizes, our approach’ alternative sizes are s j l { ( 1 2 r ) s j , ( 1 r ) s j , s j , ( 1 + r ) s j , ( 1 + 2 r ) s j } , where l L , 0.0 < r < 0.5 , and r is the range of difference between available sizes. When r = 0 , the case is the same as the traditional LRP. We consider that r = 0.25 .
  • Traditional LRP instances consider a fixed cost ( f j ) incurred whenever a depot j J is open. We keep this parameter unaltered. Additionally, we introduce a variable cost ( o j l ) depending on f j and s j l , namely: o j l = ( s j l s j ) 2 s j j f j | J | . This formula preserves o j l in the same order as f j for each depot j J . Besides, it yields negative costs whenever s j l < s j , positive costs whenever s j l > s j , and a null cost when s j l = s j . Thus our results can be compared with those found in the LRP literature.
  • An uncertain demand D i for each customer i I is considered. The demand of half of the customers is assumed to follow a log-normal probability distribution. If ϕ i is the deterministic demand in the Akca’s set, then E [ D i ] = ϕ i . In addition, three different values of variance are considered: low, medium and high, i.e., for λ { 0.05 , 0.10 , 0.20 } , V a r [ D i ] = λ ϕ i . These variability values are preserved identical to the ones used by Tordecilla et al. [13], in order to perform a suitable results comparison. The demand of the other half of the customers is considered to be fuzzy. In this case, D i can be estimated as low (DL), medium (DM) or high (DH). The demand in each of these fuzzy sets is represented by a triangular fuzzy number D i = ( d 1 i , d 2 i , d 3 i ) . If q is the vehicle total load capacity, all fuzzy demand values are expressed as a proportion of q in order to perform an appropriate comparison between the demand and the vehicle available capacity, i.e., 0 D i 1 . The membership function of these fuzzy sets are displayed in Figure 2.

Fuzzy Approach for the Demand and the Vehicle Available Capacity

When considering customers with stochastic demands, the decision about visiting the next customer in a route is made simply by comparing its expected demand with the vehicle’s current capacity. If this demand is greater, the vehicle will perform a detour to the depot for a replenishment. Nevertheless, when the next customer demand is fuzzy, the decision about serving it is made employing a preference index p i [77]. It indicates the strength of our inclination to visit the next node in a route. This index depends on both the estimated demand of the next node D i + 1 and the vehicle capacity C i that remains available after serving the customer i I . C i is expressed as a proportion of q, i.e., 0 C i 1 . It also can be treated as low (CL), medium (CM) or high (CH), and it is represented by a triangular fuzzy number C i = ( c 1 i , c 2 i , c 3 i ) . The membership function of the capacity fuzzy sets are displayed in Figure 3.
The preference index is defined between 0 and 1, i.e., 0 p i 1 . When p i = 1 , we will definitely visit the next node in a route since the vehicle available capacity can for sure meet its demand. When p i = 0 , we are sure that D i + 1 exceeds C i and the vehicle must return to the depot for a replenishment. We consider that the preference can be very low (PVL), low (PL), medium (PM), high (PH) or very high (PVH). Each of these categories is represented by a fuzzy set, whose membership function is depicted in Figure 4. Additionally, we define a set of reasoning rules (Table 2) to determine the preference to visit the next node depending on the levels of both the demand and the vehicle available capacity.
Figure 5 displays the procedure used to compute the preference index p i after serving the customer i I . This procedure is described as follows:
  • Simulate the actual demand of each customer employing a fuzzy simulation approach. Based on the works by Teodorović and Pavković [77], Sun et al. [59] and Sun [78], we follow the steps described below:
    (a)
    Generate a random demand d i between a lower bound and an upper bound. Since the objective is preserving the variability conditions similar to the stochastic demands, the lower and upper bounds are given by the expressions ϕ i 3 λ ϕ i q and ϕ i + 3 λ ϕ i q , respectively.
    (b)
    Calculate the membership degree μ ( d i ) of this demand. Notice that μ ( d i ) [ 0 , 1 ] .
    (c)
    Generate a random number ρ [ 0 , 1 ] .
    (d)
    Compare ρ and μ ( d i ) . If ρ μ ( d i ) , then assume the actual demand of the customer i as d i ; otherwise, repeat steps (a)–(d) until this condition is fulfilled.
  • Calculate the vehicle available capacity subtracting from q the sum of the simulated demand of the first m customers visited in the current route, including the customer i. Whenever the route fails and the vehicle must perform a trip to the depot for a replenishment, the counting of m starts again from 1.
  • Estimate the fuzzy demand and the fuzzy available capacity according to the categories previously defined: low, medium or high.
  • Determine the membership function of the preference index using the reasoning rules defined in Table 2.
  • Calculate a crisp preference index using the center of gravity as defuzzification method. Additional methods can be found in Klir and Yuan [79], and Opricovic and Tzeng [80].
We define a known threshold p * , such that 0 p * 1 . The computed preference index p i must be compared with p * in order to make a decision about the vehicle next destination. If p i p * , the vehicle should visit the next customer directly; otherwise, we estimate that the vehicle available capacity cannot meet the next customer demand. In this case, both preventive ( c p r e v ) and reactive ( c r e a c ) costs are calculated (see Section 3). If c p r e v < c r e a c , the vehicle should perform a detour to the depot for a preventive replenishment; otherwise, it should visit the next customer directly and react to its real demand. The lower the threshold level, the greater the inclination to unload the vehicle as much as possible before making a replenishment trip to the depot. In this case, less preventive detours are performed. Hence, the number of times that a reactive round-trip must be carried out increases. Previous tests using modified Akca’s instances yielded lower costs when p * = 0.45 .
The following parameters are used by our algorithm to run the experiments: (i) 350 iterations for map perturbations; (ii) 150 iterations for the biased-randomized savings heuristic; (iii) 150 iterations for splitting; (iv) a random value between 0.05 and 0.80 for β 1 , the parameter of the geometric distribution associated with the biased-randomized selection during the allocation map process; (v) a random value between 0.07 and 0.23 for β 2 , the parameter of the geometric distribution associated with the biased-randomized heuristic for routing; (vi) n = 100 runs for the initial simulation stage; (vii) N = 5000 runs for the intensive simulation stage; and (viii) 100 iterations to estimate the safety stock (SS), testing only discrete values between 0% and 10%. Our proposed algorithm was coded as a Java application. All experiments were executed on a standard Windows PC with a Core i 5 processor and 6 GB RAM. A total of ten different random seeds were used for each instance.

6. Results and Discussion

Table 3 shows our obtained results for 12 Akca’s instances. Five main indicators are computed: depots opening costs (OC), which is formed by both fixed and variable costs; routing costs (RC); failure costs (FC), which is incurred whenever the vehicle must perform either a detour or a round-trip to the depot; total costs (TC); and the estimated safety stock (SS) level. Four types of solutions are compared. All of them are flexible, i.e., they consider facility sizing decisions. Firstly, our best deterministic solutions are shown, i.e, there is no uncertainty in the customers’ demand and its realization is exactly as expected. In this case, a safety stock is not necessary and there are no failure costs. Secondly, we show the best stochastic solutions reported by Tordecilla et al. [13], in which the exact customers’ demand is not known. Instead, all of them follow a log-normal distribution with known mean and standard deviation. Thirdly, our best hybrid fuzzy-stochastic solutions are displayed, in which half of the customers’ demand follows a log-normal distribution, and half of the customers’ demand is considered to be fuzzy. Finally, our best fuzzy solutions are shown, in which all customers’ demand is considered to be fuzzy, due to a high level of uncertainty. Additionally, results for three levels of variability ( λ ) are shown. Clearly, our best deterministic solutions are the same regardless of the variability level, given the total absence of uncertainty.
Results in Table 3 show a slight average increase in total costs when increasing the variability level for all types of solutions, except the best deterministic solution. This growth is caused mainly by the rise in failure costs, since a greater number of detours and round-trips is expected when the demand variability level is higher. Additionally, total costs also increase when the uncertainty level is higher regardless of the variability level, i.e., the deterministic solution is the cheapest one, and the fuzzy solution is the most costly. If we compare only the average deterministic cost of each set of solutions, formed by the sum of OC and RC, we obtain values with negligible differences. Hence, the contrasts in total costs are caused mainly by failure costs. For example, for the instance Cr30x5a-3 in the low variability scenario, 1.6% of total costs are failure costs in the best stochastic solution. However, in the best fuzzy solution this percentage rises to 3.5%. Most instances show this steady growth when increasing the uncertainty level, which confirms that fuzzy scenarios have a higher uncertainty level when compared with deterministic and stochastic scenarios. Finally, the average safety stock increases when both variability and uncertainty levels rise, since more protection against uncertainty is necessary in both cases.
Results corresponding to our best deterministic solution in Table 3 were yielded assuming that the realized demand is deterministic. Hence, an additional experiment has been performed, in which this solution (called henceforth OBD) is tested in a hybrid fuzzy-stochastic environment, using 0% of safety stock protection against uncertainty. Figure 6 compares this solution’s results with our best-found hybrid fuzzy-stochastic solution (OBF) in terms of failure costs. Results for 12 Akca’s instances are depicted for each demand variability scenario. Extreme points in dashed lines indicate the average cost for each set of data. As expected, average failure costs show an increasing trend when the variability grows, regardless of the type of solution. Conversely, Figure 6 shows that OBF outperforms OBD when tested under uncertainty conditions. This fact demonstrates the quality of our fuzzy simheuristic approach, especially in scenarios where the demand variability is high.
Table 4 compares two types of hybrid fuzzy-stochastic solutions. Firstly, we show our best solution with a single facility size alternative given by the original Akca’s instances—i.e., the solution is not flexible since only one size is available to select. Secondly, we show our best flexible solution, which corresponds to our best hybrid solution in Table 3. When comparing the total costs of both types of solutions, the negative gap obtained for all instances and under all variability levels shows the advantages of considering facility sizing decisions. For example, we reach a maximum absolute gap of 7.71 % in total cost savings for a single instance. In average, both opening and routing costs decrease whenever alternative depot sizes are available. Nevertheless, each instance shows different results regarding OC and RC. The most evident case is that in which opening costs decrease. Clearly, this is a direct result of having smaller facility size alternatives. Without loss of generality, all examples below take as reference the high variability scenario. For example, the instance Cr30x5b-3 has a total demand of 1620. Both flexible and nonflexible approaches design the same routes and yield equal routing costs. Nevertheless, the nonflexible approach locates two depots of size 1000 each. Conversely, our flexible approach locates one depot of size 1000 and one depot of size 750. Hence, the nonflexible solution assigns an extra capacity that is not necessary under the problem’s current conditions.
Some instances show an opposite behavior, i.e., opening costs either increase or remain the same while routing costs decrease. For example, the nonflexible solution of the instance Cr30x5a-1 opens two depots of size 1000 each. Alternatively, the flexible solution opens one depot of size 1500 and one depot of size 500, i.e., the total capacity is equal and, given our defined costs structure, also the opening costs. However, this slight change drives a redesign of routes that decreases RC. An additional example is given by the instance Cr40x5a-2. Figure 7 depicts the best solution found by both the nonflexible approach (a) and our flexible approach (b). The solution in Figure 7a locates two depots of size 1750 each, and the solution in Figure 7b locates three depots of size 875 each. The latter case has a total capacity that is smaller than the former’s; however, opening costs are higher since the fixed cost is clearly greater when 3 facilities are open instead of 2. This new configuration decreases considerably routing costs (Table 4), which shows that considering facility sizing decisions not only reduces total costs by decreasing depots capacity but also by increasing it, since shorter routes can be designed.

Managerial Insights

From a managerial perspective, we have shown a general algorithm useful to solve a flexible-size LRP where a subset of customers provides enough information to model stochastically their demand, while the complementary subset provides scarce data. In this case, decision makers may estimate a fuzzy demand. Our algorithm is general because scenarios where the demand of all customers is deterministic, stochastic or fuzzy represent particular cases of our described problem. Hence, decision makers can employ our approach more extensively than other algorithms. We analyze these scenarios through some numerical results and assess how the level of uncertainty influences opening, routing and route-failure costs. Clearly, more precise data decrease total costs. Furthermore, we have calculated the cost of assuming a deterministic demand when the real scenario is fuzzy or stochastic. It has been shown that our hybrid approach yields less average costs, which leads to a more competitive supply chain. Additionally, we have also shown that important cost savings are generated whenever a set of facility size alternatives are analyzed by decision makers, instead of considering a single alternative—as in most LRP studies. Finally, our algorithm is able to generate detailed information about the location-allocation-routing decisions that should be made.

7. Conclusions

This work presented a location routing problem where the facility size is an additional variable, instead of a known parameter as the traditional LRP assumes. Moreover, we consider a hybrid fuzzy-stochastic setting in which some customers’ demands are fuzzy and others are stochastic. Hence, a fuzzy simheuristic approach is proposed to solve this problem cost- and time- efficiently. Initially, our algorithm selects the best size for each open facility from a set of provided alternatives. We perform an iterative procedure in which a set of location-allocation-routing configurations are assessed in terms of opening and routing costs. Then a top list of complete LRP solutions is iteratively perturbed and simulated. The perturbation stage is performed by employing an iterated local search metaheuristic. The simulation stage is carried out by running a classic Monte Carlo simulation for the stochastic demands and a fuzzy simulation for the fuzzy demands. Failure costs are introduced as an additional performance indicator. Finally, a set of elite solutions is assessed through a refinement procedure where a larger number of simulation runs is executed.
Our fuzzy simheuristic approach has been proved to be flexible enough not only to combine efficiently stochastic and fuzzy demands in a single execution but also to address less general scenarios in which demands of all customers are either deterministic or fuzzy. Our approach has also been proved to be a cost-efficient algorithm when considering uncertainty scenarios. It decreases route failure costs when compared with the best deterministic solution tested in a hybrid fuzzy-stochastic environment. The use of a safety stock policy as a protection against uncertainty has also contributed to this decrease. In order to design a time-efficient algorithm, our current approach employs stochastic and fuzzy simulation only to assess the designed routes. Hence, our algorithm results can be enhanced by introducing fuzzy-stochastic aspects from the construction stage. However, this approach might also increase computational times.
To the best of our knowledge, this is the first time that a hybrid fuzzy-stochastic LRP with facility sizing decisions is addressed. Medium-sized benchmark instances considering three demand variability levels were used. Obtained results show that introducing such flexibility decreases total costs in two mutually nonexclusive ways: firstly, yielding savings in opening costs by locating facilities of smaller size; and secondly, yielding savings in routing costs by locating facilities of higher size, which drives a routes redesign that reduces the total traveled distance. We also have demonstrated that these savings are always incurred regardless of the demand variability level.
Multiple challenges remain open for future research. Since we are considering that only routes fail when demands are higher than expected, future work can include the simulation of facility failures, which would prompt a revision of location-allocation decisions. In addition, failure costs are currently measured only by considering the distances traveled to perform round-trips and detours. Still, real-life customers might not allow a delivery delay, e.g., because a time windows constraint must be met. This delay may drive lost sales or a goodwill reduction. Hence, this type of costs can be included in the computation of failure costs. Finally, large-sized instances can be used to assess the influence of the number of nodes in our approach performance.

Author Contributions

Conceptualization, R.D.T. and A.A.J.; methodology, R.D.T., P.J.C.-M. and J.P.; software, P.J.C.-M. and J.P.; formal analysis, R.D.T.; investigation, R.D.T., P.J.C.-M., C.L.Q.-A. and J.R.M.-T.; writing—original draft preparation, all authors; writing—review and editing, R.D.T. and A.A.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been partially supported by the Spanish Ministry of Science (PID2019-111100RB-C21/AEI/10.13039/501100011033). In addition, it has received the support of the Doctoral School at the Universitat Oberta de Catalunya (Spain) and the Universidad de La Sabana (INGPhD-12-2020).

Data Availability Statement

Data are available upon reasonable request to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Salhi, S.; Rand, G.K. The effect of ignoring routes when locating depots. Eur. J. Oper. Res. 1989, 39, 150–156. [Google Scholar] [CrossRef]
  2. Nagy, G.; Salhi, S. Location-routing: Issues, models and methods. Eur. J. Oper. Res. 2007, 177, 649–672. [Google Scholar] [CrossRef] [Green Version]
  3. Prodhon, C.; Prins, C. A survey of recent research on location-routing problems. Eur. J. Oper. Res. 2014, 238, 1–17. [Google Scholar] [CrossRef]
  4. Zhou, L.; Lin, Y.; Wang, X.; Zhou, F. Model and algorithm for bilevel multisized terminal location-routing problem for the last mile delivery. Int. Trans.Oper. Res. 2019, 26, 131–156. [Google Scholar] [CrossRef]
  5. Tordecilla-Madera, R.; Polo, A.; Muñoz, D.; González-Rodríguez, L. A robust design for a Colombian dairy cooperative’s milk storage and refrigeration logistics system using binary programming. Int. J. Prod. Econ. 2017, 183, 710–720. [Google Scholar] [CrossRef]
  6. Correia, I.; Melo, T. Multi-period capacitated facility location under delayed demand satisfaction. Eur. J. Oper. Res. 2016, 255, 729–746. [Google Scholar] [CrossRef] [Green Version]
  7. Hemmelmayr, V.; Smilowitz, K.; de la Torre, L. A periodic location routing problem for collaborative recycling. IISE Trans. 2017, 49, 414–428. [Google Scholar] [CrossRef] [Green Version]
  8. Tunalıoğlu, R.; Koç, Ç.; Bektaş, T. A multiperiod location-routing problem arising in the collection of Olive Oil Mill Wastewater. J. Oper. Res. Soc. 2016, 67, 1012–1024. [Google Scholar] [CrossRef]
  9. Quintero-Araujo, C.L.; Guimarans, D.; Juan, A.A. A simheuristic algorithm for the capacitated location routing problem with stochastic demands. J. Simul. 2019, 1–18. [Google Scholar] [CrossRef]
  10. Rabbani, M.; Heidari, R.; Yazdanparast, R. A stochastic multi-period industrial hazardous waste location-routing problem: Integrating NSGA-II and Monte Carlo simulation. Eur. J. Oper. Res. 2019, 272, 945–961. [Google Scholar] [CrossRef]
  11. Sun, Z.; Yan, N.; Sun, Y.; Li, H. Location-routing optimization with split demand for customer self-pickup via data analysis and heuristics search. Asia-Pac. J. Oper. Res. 2019, 36, 1940013. [Google Scholar] [CrossRef]
  12. Zhang, S.; Chen, M.; Zhang, W. A novel location-routing problem in electric vehicle transportation with stochastic demands. J. Clean. Prod. 2019, 221, 567–581. [Google Scholar] [CrossRef]
  13. Tordecilla, R.D.; Panadero, J.; Quintero-Araujo, C.L.; Montoya-Torres, J.R.; Juan, A.A. A simheuristic algorithm for the location routing problem with facility sizing decisions and stochastic demands. In Proceedings of the 2020 Winter Simulation Conference, IEEE, Marriott, Orlando, FL, USA, 14–18 December 2020; pp. 1265–1275. [Google Scholar]
  14. Herazo-Padilla, N.; Montoya-Torres, J.R.; Nieto Isaza, S.; Alvarado-Valencia, J. Simulation-optimization approach for the stochastic location-routing problem. J. Simul. 2015, 9, 296–311. [Google Scholar] [CrossRef]
  15. Zhang, B.; Li, H.; Li, S.; Peng, J. Sustainable multi-depot emergency facilities location-routing problem with uncertain information. Appl. Math. Comput. 2018, 333, 506–520. [Google Scholar] [CrossRef]
  16. Corlu, C.G.; Panadero, J.; Onggo, S.; Juan, A.A. On the scarcity of observations when modelling random inputs and the quality of solutions to stochastic optimisation problems. In Proceedings of the 2020 Winter Simulation Conference, IEEE, Marriott, Orlando, FL, USA, 14–18 December 2020; pp. 2105–2113. [Google Scholar]
  17. Zhang, H.; Liu, F.; Ma, L.; Zhang, Z. A hybrid heuristic based on a particle swarm algorithm to solve the capacitated location-routing problem with fuzzy demands. IEEE Access 2020, 8, 153671–153691. [Google Scholar] [CrossRef]
  18. Mehrjerdi, Y.Z.; Nadizadeh, A. Using greedy clustering method to solve capacitated location-routing problem with fuzzy demands. Eur. J. Oper. Res. 2013, 229, 75–84. [Google Scholar] [CrossRef]
  19. Fazayeli, S.; Eydi, A.; Kamalabadi, I.N. Location-routing problem in multimodal transportation network with time windows and fuzzy demands: Presenting a two-part genetic algorithm. Comput. Ind. Eng. 2018, 119, 233–246. [Google Scholar] [CrossRef]
  20. Nadizadeh, A.; Kafash, B. Fuzzy capacitated location-routing problem with simultaneous pickup and delivery demands. Transp. Lett. 2019, 11, 1–19. [Google Scholar] [CrossRef]
  21. Zarandi, M.H.F.; Hemmati, A.; Davari, S. The multi-depot capacitated location-routing problem with fuzzy travel times. Expert Syst. Appl. 2011, 38, 10075–10084. [Google Scholar] [CrossRef]
  22. Zarandi, M.H.F.; Hemmati, A.; Davari, S.; Turksen, I.B. Capacitated location-routing problem with time windows under uncertainty. Knowl. Based Syst. 2013, 37, 480–489. [Google Scholar] [CrossRef] [Green Version]
  23. Ghezavati, V.; Morakabatchian, S. Application of a fuzzy service level constraint for solving a multi-objective location-routing problem for the industrial hazardous wastes. J. Intell. Fuzzy Syst. 2015, 28, 2003–2013. [Google Scholar] [CrossRef]
  24. Juan, A.A.; Kelton, W.D.; Currie, C.S.; Faulin, J. Simheuristics applications: Dealing with uncertainty in logistics, transportation, and other supply chain areas. In Proceedings of the 2018 Winter Simulation Conference, IEEE, Marriott, Orlando, FL, USA, 14–18 December 2018; pp. 3048–3059. [Google Scholar]
  25. Maranzana, F. On the location of supply points to minimize transport costs. J. Oper. Res. Soc. 1964, 15, 261–270. [Google Scholar] [CrossRef]
  26. Dai, Z.; Aqlan, F.; Gao, K.; Zhou, Y. A two-phase method for multi-echelon location-routing problems in supply chains. Expert Syst. Appl. 2019, 115, 618–634. [Google Scholar] [CrossRef]
  27. Quintero-Araujo, C.L.; Caballero-Villalobos, J.P.; Juan, A.A.; Montoya-Torres, J.R. A biased-randomized metaheuristic for the capacitated location routing problem. Int. Trans. Oper. Res. 2017, 24, 1079–1098. [Google Scholar] [CrossRef]
  28. Escobar, J.W.; Linfati, R.; Baldoquin, M.G.; Toth, P. A granular variable tabu neighborhood search for the capacitated location-routing problem. Transp. Res. Part B Methodol. 2014, 67, 344–356. [Google Scholar] [CrossRef]
  29. Ferdi, I.; Layeb, A. A GRASP algorithm based new heuristic for the capacitated location routing problem. J. Exp. Theor. Artif. Intell. 2018, 30, 369–387. [Google Scholar] [CrossRef]
  30. Quintero-Araujo, C.L.; Gruler, A.; Juan, A.A.; Faulin, J. Using horizontal cooperation concepts in integrated routing and facility-location decisions. Int. Trans. Oper. Res. 2019, 26, 551–576. [Google Scholar] [CrossRef]
  31. Hof, J.; Schneider, M.; Goeke, D. Solving the battery swap station location-routing problem with capacitated electric vehicles using an AVNS algorithm for vehicle-routing problems with intermediate stops. Transp. Res. Part B Methodol. 2017, 97, 102–112. [Google Scholar] [CrossRef]
  32. Almouhanna, A.; Quintero-Araujo, C.L.; Panadero, J.; Juan, A.A.; Khosravi, B.; Ouelhadj, D. The location routing problem using electric vehicles with constrained distance. Comput. Oper. Res. 2020, 115, 104864. [Google Scholar] [CrossRef]
  33. Theeraviriya, C.; Sirirak, W.; Praseeratasang, N. Location and routing planning considering electric vehicles with restricted distance in agriculture. World Electr. Veh. J. 2020, 11, 61. [Google Scholar] [CrossRef]
  34. Nataraj, S.; Ferone, D.; Quintero-Araujo, C.; Juan, A.; Festa, P. Consolidation centers in city logistics: A cooperative approach based on the location routing problem. Int. J. Ind. Eng. Comput. 2019, 10, 393–404. [Google Scholar] [CrossRef]
  35. Ukkusuri, S.V.; Yushimito, W.F. Location routing approach for the humanitarian prepositioning problem. Transp. Res. Record 2008, 2089, 18–25. [Google Scholar] [CrossRef]
  36. Lashine, S.H.; Fattouh, M.; Issa, A. Location/allocation and routing decisions in supply chain network design. J. Model. Manag. 2006, 1, 173–183. [Google Scholar] [CrossRef]
  37. Leng, L.; Zhao, Y.; Zhang, J.; Zhang, C. An effective approach for the multiobjective regional low-carbon location-routing problem. Int. J. Environ. Res. Public Health 2019, 16, 2064. [Google Scholar] [CrossRef] [Green Version]
  38. Wang, Z.; Leng, L.; Wang, S.; Li, G.; Zhao, Y. A hyperheuristic approach for location-routing problem of cold chain logistics considering fuel consumption. Comput. Intell. Neurosci. 2020, 2020, 8395754. [Google Scholar] [CrossRef]
  39. Rabbani, M.; Sadati, S.A.; Farrokhi-Asl, H. Incorporating location routing model and decision making techniques in industrial waste management: Application in the automotive industry. Comput. Ind. Eng. 2020, 148, 106692. [Google Scholar] [CrossRef]
  40. Ghaffari-Nasab, N.; Ahari, S.G.; Ghazanfari, M. A hybrid simulated annealing based heuristic for solving the location-routing problem with fuzzy demands. Sci. Iran. 2013, 20, 919–930. [Google Scholar]
  41. Nadizadeh, A.; Nasab, H.H. Solving the dynamic capacitated location-routing problem with fuzzy demands by hybrid heuristic algorithm. Eur. J. Oper. Res. 2014, 238, 458–470. [Google Scholar] [CrossRef]
  42. Wei, M.; Yu, L.; Li, X. Credibilistic Location-Routing Model for Hazardous Materials Transportation. Int. J. Intell. Syst. 2015, 30, 23–39. [Google Scholar] [CrossRef]
  43. Faulin, J.; Gilibert, M.; Juan, A.A.; Vilajosana, X.; Ruiz, R. SR-1: A simulation-based algorithm for the capacitated vehicle routing problem. In Proceedings of the 2008 Winter Simulation Conference, IEEE, Miami, FL, USA, 7–10 December 2008; pp. 2708–2716. [Google Scholar]
  44. Juan, A.A.; Faulin, J.; Ruiz, R.; Barrios, B.; Gilibert, M.; Vilajosana, X. Using oriented random search to provide a set of alternative solutions to the capacitated vehicle routing problem. In Operations Research and Cyber-Infrastructure; Springer: Berlin, Germany, 2009; pp. 331–345. [Google Scholar]
  45. Oliva, D.; Copado, P.; Hinojosa, S.; Panadero, J.; Riera, D.; Juan, A.A. Fuzzy simheuristics: Solving optimization problems under stochastic and uncertainty scenarios. Mathematics 2021, 1, 00005. [Google Scholar]
  46. Quintero-Araujo, C.L.; Gruler, A.; Juan, A.A.; de Armas, J.; Ramalhinho, H. Using simheuristics to promote horizontal collaboration in stochastic city logistics. Prog. Artif. Intell. 2017, 6, 275–284. [Google Scholar] [CrossRef]
  47. Gruler, A.; Panadero, J.; de Armas, J.; Moreno, J.A.; Juan, A.A. A variable neighborhood search simheuristic for the multiperiod inventory routing problem with stochastic demands. Int. Trans. Oper. Res. 2020, 27, 314–335. [Google Scholar] [CrossRef]
  48. Gonzalez-Martin, S.; Juan, A.A.; Riera, D.; Elizondo, M.G.; Ramos, J.J. A simheuristic algorithm for solving the arc routing problem with stochastic demands. J. Simul. 2018, 12, 53–66. [Google Scholar] [CrossRef]
  49. Gruler, A.; Panadero, J.; de Armas, J.; Moreno, J.A.; Juan, A.A. Combining variable neighborhood search with simulation for the inventory routing problem with stochastic demands and stock-outs. Comput. Ind. Eng. 2018, 123, 278–288. [Google Scholar] [CrossRef]
  50. Gruler, A.; Fikar, C.; Juan, A.A.; Hirsch, P.; Contreras-Bolton, C. Supporting multi-depot and stochastic waste collection management in clustered urban areas via simulation–optimization. J. Simul. 2017, 11, 11–19. [Google Scholar] [CrossRef]
  51. Uddin, M.; Huynh, N. Reliable routing of road-rail intermodal freight under uncertainty. Netw. Spat. Econ. 2019, 19, 929–952. [Google Scholar] [CrossRef]
  52. Hrušovskỳ, M.; Demir, E.; Jammernegg, W.; Van Woensel, T. Hybrid simulation and optimization approach for green intermodal transportation problem with travel time uncertainty. Flex. Serv. Manuf. J. 2018, 30, 486–516. [Google Scholar] [CrossRef] [Green Version]
  53. Zhao, Y.; Liu, R.; Zhang, X.; Whiteing, A. A chance-constrained stochastic approach to intermodal container routing problems. PLoS ONE 2018, 13, e0192275. [Google Scholar] [CrossRef]
  54. Werners, B.; Drawe, M. Capacitated vehicle routing problem with fuzzy demand. In Fuzzy Sets Based Heuristics for Optimization; Springer: Berlin, Germany, 2003; pp. 317–335. [Google Scholar]
  55. Erbao, C.; Mingyong, L. A hybrid differential evolution algorithm to vehicle routing problem with fuzzy demands. J. Comput. Appl. Math. 2009, 231, 302–310. [Google Scholar] [CrossRef] [Green Version]
  56. Xue, L.; Dai, X.X. Research on the vehicle routing problem with fuzzy demands. In Advanced Materials Research; Trans Tech Publications Ltd.: Stafa-Zurich, Switzerland, 2011; Volume 186, pp. 570–575. [Google Scholar]
  57. Sun, Y. Fuzzy approaches and simulation-based reliability modeling to solve a Road–Rail intermodal routing problem with soft delivery time windows when demand and capacity are uncertain. Int. J. Fuzzy Syst. 2020, 22, 2119–2148. [Google Scholar] [CrossRef]
  58. Zheng, Y.; Liu, B. Fuzzy vehicle routing model with credibility measure and its hybrid intelligent algorithm. Appl. Math. Comput. 2006, 176, 673–683. [Google Scholar] [CrossRef]
  59. Sun, Y.; Hrušovskỳ, M.; Zhang, C.; Lang, M. A time-dependent fuzzy programming approach for the green multimodal routing problem with rail service capacity uncertainty and road traffic congestion. Complexity 2018, 2018, 8645793. [Google Scholar] [CrossRef] [Green Version]
  60. Gupta, R.; Singh, B.; Pandey, D. Fuzzy vehicle routing problem with uncertainty in service time. Int. J. Contemp. Math. Sci. 2010, 5, 497–507. [Google Scholar]
  61. Tang, J.; Pan, Z.; Fung, R.Y.; Lau, H. Vehicle routing problem with fuzzy time windows. Fuzzy Sets Syst. 2009, 160, 683–695. [Google Scholar] [CrossRef]
  62. López-Castro, L.F.; Montoya-Torres, J.R. Vehicle routing with fuzzy time windows using a genetic algorithm. In Proceedings of the 2011 Workshop On Computational Intelligence In Production And Logistics Systems, IEEE, Paris, France, 11–15 April 2011; pp. 1–8. [Google Scholar]
  63. Cao, E.; Lai, M. The open vehicle routing problem with fuzzy demands. Expert Syst. Appl. 2010, 37, 2405–2411. [Google Scholar] [CrossRef]
  64. Gonzalez-Neira, E.M.; Ferone, D.; Hatami, S.; Juan, A.A. A biased-randomized simheuristic for the distributed assembly permutation flowshop problem with stochastic processing times. Simul. Model. Pract. Theory 2017, 79, 23–36. [Google Scholar] [CrossRef]
  65. Zarandi, M.H.F.; Asl, A.A.S.; Sotudian, S.; Castillo, O. A state of the art review of intelligent scheduling. Artif. Intell. Rev. 2020, 53, 501–593. [Google Scholar] [CrossRef]
  66. Ochoa, P.; Castillo, O.; Soria, J. High-Speed Interval Type-2 Fuzzy System for Dynamic Crossover Parameter Adaptation in Differential Evolution and Its Application to Controller Optimization. Int. J. Fuzzy Syst. 2019, 22, 414–427. [Google Scholar] [CrossRef]
  67. Bernal, E.; Castillo, O.; Soria, J.; Valdez, F. Generalized type-2 fuzzy logic in galactic swarm optimization: Design of an optimal ball and beam fuzzy controller. J. Intell. Fuzzy Syst. 2020, 39, 3545–3559. [Google Scholar] [CrossRef]
  68. Anter, A.M.; Gupta, D.; Castillo, O. A novel parameter estimation in dynamic model via fuzzy swarm intelligence and chaos theory for faults in wastewater treatment plant. Soft Comput. 2020, 24, 111–129. [Google Scholar] [CrossRef]
  69. Panadero, J.; Doering, J.; Kizys, R.; Juan, A.A.; Fito, A. A variable neighborhood search simheuristic for project portfolio selection under uncertainty. J. Heuristics 2020, 26, 353–375. [Google Scholar] [CrossRef] [Green Version]
  70. de Armas, J.; Juan, A.A.; Marquès, J.M.; Pedroso, J.P. Solving the deterministic and stochastic uncapacitated facility location problem: From a heuristic to a simheuristic. J. Oper. Res. Soc. 2017, 68, 1161–1176. [Google Scholar] [CrossRef]
  71. Panadero, J.; Juan, A.A.; Bayliss, C.; Currie, C. Maximising reward from a team of surveillance drones: A simheuristic approach to the stochastic team orienteering problem. Eur. J. Ind. Eng. 2020, 14, 485–516. [Google Scholar] [CrossRef]
  72. Ferone, D.; Gruler, A.; Festa, P.; Juan, A.A. Enhancing and extending the classical GRASP framework with biased randomisation and simulation. J. Oper. Res. Soc. 2019, 70, 1362–1375. [Google Scholar] [CrossRef]
  73. Ferrer, A.; Guimarans, D.; Ramalhinho, H.; Juan, A.A. A BRILS metaheuristic for non-smooth flow-shop problems with failure-risk costs. Expert Syst. Appl. 2016, 44, 177–186. [Google Scholar] [CrossRef] [Green Version]
  74. Akca, Z.; Berger, R.; Ralphs, T. A branch-and-price algorithm for combined location and routing problems under capacity restrictions. In Operations Research and Cyber-Infrastructure; Springer: Boston, MA, USA, 2009; pp. 309–330. [Google Scholar]
  75. Barreto, S.; Ferreira, C.; Paixao, J.; Santos, B.S. Using clustering analysis in a capacitated location-routing problem. Eur. J. Oper. Res. 2007, 179, 968–977. [Google Scholar] [CrossRef]
  76. Belenguer, J.M.; Benavent, E.; Prins, C.; Prodhon, C.; Calvo, R.W. A branch-and-cut method for the capacitated location-routing problem. Comput. Oper. Res. 2011, 38, 931–941. [Google Scholar] [CrossRef]
  77. Teodorović, D.; Pavković, G. The fuzzy set theory approach to the vehicle routing problem when demand at nodes is uncertain. Fuzzy Sets Syst. 1996, 82, 307–317. [Google Scholar] [CrossRef]
  78. Sun, Y. A Fuzzy Multi-Objective Routing Model for Managing Hazardous Materials Door-to-Door Transportation in the Road-Rail Multimodal Network With Uncertain Demand and Improved Service Level. IEEE Access 2020, 8, 172808–172828. [Google Scholar] [CrossRef]
  79. Klir, G.; Yuan, B. Fuzzy Sets and Fuzzy Logic: Theory and Applications. Possibility Theory Versus Probab. Theory 1996, 32, 207–208. [Google Scholar]
  80. Opricovic, S.; Tzeng, G.H. Defuzzification within a multicriteria decision model. Int. J. Uncertain. Fuzziness Knowl. Based Syst. 2003, 11, 635–652. [Google Scholar] [CrossRef]
Figure 1. Graphical representation of a location routing problem (LRP) solution.
Figure 1. Graphical representation of a location routing problem (LRP) solution.
Algorithms 14 00045 g001
Figure 2. Fuzzy sets for the demand of the customer i.
Figure 2. Fuzzy sets for the demand of the customer i.
Algorithms 14 00045 g002
Figure 3. Fuzzy sets for the vehicle available capacity after visiting the customer i.
Figure 3. Fuzzy sets for the vehicle available capacity after visiting the customer i.
Algorithms 14 00045 g003
Figure 4. Fuzzy sets for the preference strength to visit the customer i.
Figure 4. Fuzzy sets for the preference strength to visit the customer i.
Algorithms 14 00045 g004
Figure 5. Procedure used to compute the preference index p i .
Figure 5. Procedure used to compute the preference index p i .
Algorithms 14 00045 g005
Figure 6. Failure costs of our best deterministic and our best hybrid solutions.
Figure 6. Failure costs of our best deterministic and our best hybrid solutions.
Algorithms 14 00045 g006
Figure 7. Best-found solution by a nonflexible (a) and a flexible (b) fuzzy LRP for the instance Cr40x5a-2.
Figure 7. Best-found solution by a nonflexible (a) and a flexible (b) fuzzy LRP for the instance Cr40x5a-2.
Algorithms 14 00045 g007
Table 1. Recent works related to the location routing problem with uncertain parameters.
Table 1. Recent works related to the location routing problem with uncertain parameters.
ReferenceType of
Uncertainty
Uncertain
Parameter
Mathematical
Approach
Solving ApproachObjective Criterion
Quintero-Araujo et al. [9]StochasticDemandMixed-integer
linear programming
Simheuristic
Iterated local search
Monte Carlo simulation
Minimize cost
Rabbani et al. [10]StochasticDemand
Number of
people at risk
Mixed-integer non-
linear programming
Simheuristic
NSGA-II
Monte Carlo Simulation
Minimize cost
Minimize
environmental risk
Sun et al. [11]StochasticDemandMixed-integer
linear programming
Biogeography-based optimization
Adaptive large neighborhood search
Minimize cost
Zhang et al. [12]StochasticDemandVariable neighborhood search
Particle swarm optimization
Minimize cost
Tordecilla et al. [13]StochasticDemandSimheuristic
Iterated local search
Monte Carlo simulation
Minimize cost
Herazo-Padilla et al. [14]StochasticTransportation cost
Travel speed
Mixed-integer
linear programming
Ant colony optimization
Discrete-event simulation
Minimize cost
Zhang et al. [15]StochasticDemand
Travel distance
Depot opening cost
Mixed-integer non-
linear programming
Genetic algorithm
Uncertain simulation
Minimize travel time
Minimize emergency relief cost
Minimize CO 2 emissions
Zhang et al. [17]FuzzyDemandA fuzzy chance
constrained model
Particle swarm optimization
Variable neighborhood search
Stochastic simulation
Minimize cost
Minimize additional travel
distance due to route failures
Mehrjerdi and Nadizadeh [18]FuzzyDemandA fuzzy chance
constrained model
A greedy clustering method
Ant colony system
Stochastic simulation
Minimize cost
Fazayeli et al. [19]FuzzyDemandMixed-integer non-
linear programming
Exact approach
Genetic algorithm
Minimize cost
Nadizadeh and Kafash [20]FuzzyDemandA fuzzy chance
constrained model
A greedy clustering method
Ant colony system
Stochastic simulation
Minimize cost
Zarandi et al. [21]FuzzyTravel timeA fuzzy chance
constrained model
Simulated annealing
Fuzzy simulation
Minimize cost
Zarandi et al. [22]FuzzyDemand
Travel time
A fuzzy chance
constrained model
Simulated annealing
Fuzzy simulation
Minimize cost
Minimize additional travel
distance due to route failures
Ghezavati and Morakabatchian [23]FuzzyTime windowsMixed-integer
linear programming
Exact approachMinimize cost
Minimize risks
Ghaffari-Nasab et al. [40]FuzzyDemandA fuzzy chance
constrained model
Simulated annealing
Stochastic simulation
Minimize cost
Minimize additional travel
distance due to route failures
Nadizadeh and Nasab [41]FuzzyDemandA fuzzy chance
constrained model
A hybrid heuristic algorithm
Ant colony system
Stochastic simulation
Minimize cost
Minimize additional travel
distance due to route failures
Wei et al. [42]FuzzyTransportation cost
Number of people
that may be at risk
A fuzzy chance
constrained model
Genetic algorithm
Fuzzy simulation
Minimize cost
Minimize risks
Table 2. Reasoning rules determining the visit preference strength.
Table 2. Reasoning rules determining the visit preference strength.
DemandAvailable Capacity
CLCMCH
DLPMPHPVH
DMPLPMPH
DHPVLPLPM
Table 3. Comparative results between our flexible solutions under different uncertainty levels.
Table 3. Comparative results between our flexible solutions under different uncertainty levels.
InstanceBest Deterministic SolutionBest Stochastic Solution [13]Best Hybrid SolutionBest Fuzzy Solution
OCRCTCOCRCFCTCSSOCRCFCTCSSOCRCFCTCSS
Low variability
Cr30x5a-1200.00575.14775.14200.00575.142.37777.510%200.00575.143.31778.452%200.00575.145.86781.002%
Cr30x5a-2200.00607.28807.28200.00607.280.04807.323%200.00607.280.12807.403%200.00607.280.12807.403%
Cr30x5a-3187.50507.92695.42187.50509.2510.99707.743%187.50509.2517.48714.223%187.50509.2525.50722.253%
Cr30x5b-1225.00623.22848.22225.00623.229.37857.590%225.00623.2214.59862.810%225.00623.2222.85871.071%
Cr30x5b-2187.50625.32812.82187.50625.320.00812.822%187.50625.320.00812.822%187.50625.320.00812.822%
Cr30x5b-3187.50684.58872.08187.50684.582.25874.331%187.50684.586.35878.431%187.50684.589.50881.581%
Cr40x5a-1162.50731.84894.34162.50731.840.03894.371%162.50731.840.07894.411%162.50731.840.59894.931%
Cr40x5a-2225.00637.26862.26225.00639.020.10864.120%225.00639.020.81864.831%225.00642.020.03867.053%
Cr40x5a-3162.50752.88915.38162.50752.880.97916.350%162.50752.883.26918.640%162.50752.886.82922.211%
Cr40x5b-1162.50852.041014.54162.50852.046.901021.451%162.50852.0412.241026.781%162.50852.0420.791035.331%
Cr40x5b-2225.00690.57915.57225.00690.570.08915.651%225.00690.570.62916.181%225.00690.571.23916.791%
Cr40x5b-3175.00764.33939.33175.00772.870.07947.932%175.00772.870.29948.162%175.00772.870.35948.222%
Average191.67671.03862.70191.67672.002.76866.431.17%191.67672.004.93868.591.42%191.67672.257.80871.721.75%
Medium variability
Cr30x5a-1200.00575.14775.14200.00575.147.63782.772%200.00575.149.67784.812%200.00575.1412.91788.052%
Cr30x5a-2200.00607.28807.28200.00607.280.46807.743%200.00607.281.94809.223%200.00607.281.43808.713%
Cr30x5a-3187.50507.92695.42187.50509.2518.50715.253%187.50509.2524.10720.853%187.50509.2529.73726.483%
Cr30x5b-1225.00623.22848.22225.00623.2214.63862.850%225.00623.2218.32866.533%225.00623.2224.23872.453%
Cr30x5b-2187.50625.32812.82187.50625.320.00812.822%187.50625.320.00812.822%187.50625.320.00812.822%
Cr30x5b-3187.50684.58872.08187.50684.5810.21882.280%187.50684.5812.79884.871%187.50684.5812.88884.961%
Cr40x5a-1162.50731.84894.34162.50739.240.01901.753%162.50739.240.01901.753%162.50739.240.00901.743%
Cr40x5a-2225.00637.26862.26225.00643.523.07871.591%225.00642.020.24867.263%225.00642.020.57867.593%
Cr40x5a-3162.50752.88915.38162.50752.884.46919.851%162.50752.888.57923.951%162.50752.8811.83927.221%
Cr40x5b-1162.50852.041014.54162.50858.584.541025.622%162.50858.588.011029.092%237.50795.180.001032.684%
Cr40x5b-2225.00690.57915.57225.00690.572.06917.631%225.00690.573.77919.330%225.00690.575.80921.371%
Cr40x5b-3175.00764.33939.33175.00772.871.42949.292%175.00772.872.53950.402%175.00772.872.96950.822%
Average191.67671.03862.70191.67673.545.58870.791.67%191.67673.417.50872.572.08%197.92668.138.53874.572.33%
High variability
Cr30x5a-1200.00575.14775.14200.00575.1419.66794.802%200.00575.1419.82794.960%200.00575.1424.25799.381%
Cr30x5a-2200.00607.28807.28200.00607.740.02807.765%200.00611.410.02811.437%200.00607.740.04807.785%
Cr30x5a-3187.50507.92695.42187.50509.2527.86724.612%187.50509.2529.95726.704%187.50509.2533.41730.163%
Cr30x5b-1225.00623.22848.22225.00623.2219.99868.2110%225.00623.2220.73868.9510%225.00623.2224.86873.0810%
Cr30x5b-2187.50625.32812.82187.50625.320.10812.923%187.50625.320.20813.025%187.50625.320.15812.973%
Cr30x5b-3187.50684.58872.08187.50684.5824.93897.001%187.50684.5829.03901.115%187.50684.5834.01906.095%
Cr40x5a-1162.50731.84894.34162.50737.202.85902.552%162.50735.847.83906.171%162.50735.849.38907.711%
Cr40x5a-2225.00637.26862.26225.00642.021.79868.823%225.00642.021.48868.503%225.00642.022.25869.273%
Cr40x5a-3162.50752.88915.38162.50763.695.78931.972%162.50763.697.76933.962%162.50752.8818.65934.041%
Cr40x5b-1162.50852.041014.54237.50786.004.651028.143%237.50792.362.841032.704%237.50786.008.471031.973%
Cr40x5b-2225.00690.57915.57225.00690.579.35924.912%225.00690.5712.59928.152%225.00690.5714.96930.532%
Cr40x5b-3175.00764.33939.33175.00780.624.14959.763%175.00780.624.90960.523%175.00780.625.86961.483%
Average191.67671.03862.70197.92668.7810.09876.793.17%197.92669.5011.43878.853.83%197.92667.7714.69880.373.33%
Table 4. Comparative results between our hybrid solutions when considering facility sizing decisions.
Table 4. Comparative results between our hybrid solutions when considering facility sizing decisions.
InstanceBest Nonflexible Hybrid SolutionBest Flexible Hybrid SolutionGap
TC
OCRCFCTCSSOCRCFCTCSS
Low variability
Cr30x5a-1200.00619.513.45822.961%200.00575.143.31778.452%−5.41%
Cr30x5a-2200.00626.010.04826.051%200.00607.280.12807.403%−2.26%
Cr30x5a-3200.00507.9917.56725.552%187.50509.2517.48714.223%−1.56%
Cr30x5b-1200.00682.970.32883.292%225.00623.2214.59862.810%−2.32%
Cr30x5b-2200.00625.320.00825.322%187.50625.320.00812.822%−1.51%
Cr30x5b-3200.00684.585.95890.531%187.50684.586.35878.431%−1.36%
Cr40x5a-1200.00733.473.22936.700%162.50731.840.07894.411%−4.51%
Cr40x5a-2200.00691.4711.15902.631%225.00639.020.81864.831%−4.19%
Cr40x5a-3200.00748.649.88958.521%162.50752.883.26918.640%−4.16%
Cr40x5b-1200.00858.581.941060.532%162.50852.0412.241026.781%−3.18%
Cr40x5b-2300.00690.570.65991.222%225.00690.570.62916.181%−7.57%
Cr40x5b-3200.00780.620.07980.692%175.00772.870.29948.162%−3.32%
Average208.33687.484.52900.331.42%191.67672.004.93868.591.42%−3.45%
Medium variability
Cr30x5a-1200.00619.519.17828.680%200.00575.149.67784.812%−5.29%
Cr30x5a-2200.00626.010.60826.612%200.00607.281.94809.223%−2.10%
Cr30x5a-3200.00507.9924.30732.292%187.50509.2524.10720.853%−1.56%
Cr30x5b-1200.00681.5014.31895.801%225.00623.2218.32866.533%−3.27%
Cr30x5b-2200.00625.320.01825.332%187.50625.320.00812.822%−1.52%
Cr30x5b-3200.00684.5815.60900.181%187.50684.5812.79884.871%−1.70%
Cr40x5a-1200.00733.477.69941.171%162.50739.240.01901.753%−4.19%
Cr40x5a-2200.00700.8012.59913.393%225.00642.020.24867.263%−5.05%
Cr40x5a-3200.00748.6420.15968.790%162.50752.888.57923.951%−4.63%
Cr40x5b-1200.00863.912.321066.233%162.50858.588.011029.092%−3.48%
Cr40x5b-2300.00690.574.18994.751%225.00690.573.77919.330%−7.58%
Cr40x5b-3200.00780.620.94981.563%175.00772.872.53950.402%−3.17%
Average208.33688.589.32906.231.58%191.67673.417.50872.572.08%−3.63%
High variability
Cr30x5a-1200.00619.5120.69840.200%200.00575.1419.82794.960%−5.38%
Cr30x5a-2200.00621.455.66827.123%200.00611.410.02811.437%−1.90%
Cr30x5a-3200.00507.9930.16738.154%187.50509.2529.95726.704%−1.55%
Cr30x5b-1200.00681.5018.85900.350%225.00623.2220.73868.9510%−3.49%
Cr30x5b-2200.00625.320.14825.465%187.50625.320.20813.025%−1.51%
Cr30x5b-3200.00684.5830.23914.811%187.50684.5829.03901.115%−1.50%
Cr40x5a-1200.00737.945.78943.732%162.50735.847.83906.171%−3.98%
Cr40x5a-2200.00700.8015.98916.783%225.00642.021.48868.503%−5.27%
Cr40x5a-3200.00748.6432.89981.540%162.50763.697.76933.962%−4.85%
Cr40x5b-1200.00858.5822.531081.112%237.50792.362.841032.704%−4.48%
Cr40x5b-2300.00693.0312.661005.690%225.00690.5712.59928.152%−7.71%
Cr40x5b-3200.00772.8713.22986.092%175.00780.624.90960.523%−2.59%
Average208.33687.6817.40913.421.83%197.92669.5011.43878.853.83%−3.68%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tordecilla, R.D.; Copado-Méndez, P.J.; Panadero, J.; Quintero-Araujo, C.L.; Montoya-Torres, J.R.; Juan, A.A. Combining Heuristics with Simulation and Fuzzy Logic to Solve a Flexible-Size Location Routing Problem under Uncertainty. Algorithms 2021, 14, 45. https://0-doi-org.brum.beds.ac.uk/10.3390/a14020045

AMA Style

Tordecilla RD, Copado-Méndez PJ, Panadero J, Quintero-Araujo CL, Montoya-Torres JR, Juan AA. Combining Heuristics with Simulation and Fuzzy Logic to Solve a Flexible-Size Location Routing Problem under Uncertainty. Algorithms. 2021; 14(2):45. https://0-doi-org.brum.beds.ac.uk/10.3390/a14020045

Chicago/Turabian Style

Tordecilla, Rafael D., Pedro J. Copado-Méndez, Javier Panadero, Carlos L. Quintero-Araujo, Jairo R. Montoya-Torres, and Angel A. Juan. 2021. "Combining Heuristics with Simulation and Fuzzy Logic to Solve a Flexible-Size Location Routing Problem under Uncertainty" Algorithms 14, no. 2: 45. https://0-doi-org.brum.beds.ac.uk/10.3390/a14020045

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop