Advances in Modelling, Monitoring, and Control for Complex Industrial Systems
View this Special IssueResearch Article  Open Access
Aaron Guerrero Campanur, Elias OlivaresBenitez, Pablo A. Miranda, Rodolfo Eleazar PerezLoaiza, Jose Humberto AblanedoRosas, "Design of a Logistics Nonlinear System for a Complex, Multiechelon, Supply Chain Network with Uncertain Demands", Complexity, vol. 2018, Article ID 4139601, 16 pages, 2018. https://doi.org/10.1155/2018/4139601
Design of a Logistics Nonlinear System for a Complex, Multiechelon, Supply Chain Network with Uncertain Demands
Abstract
Industrial systems, such as logistics and supply chain networks, are complex systems because they comprise a big number of interconnected actors and significant nonlinear and stochastic features. This paper analyzes a distribution network design problem for a fourechelon supply chain. The problem is represented as an inventorylocation model with uncertain demand and a continuous review inventory policy. The decision variables include location at the intermediate levels and product flows between echelons. The related safety and cyclic inventory levels can be computed from these decision variables. The problem is formulated as a mixed integer nonlinear programming model to find the optimal design of the distribution network. A linearization of the nonlinear model based on a piecewise linear approximation is proposed. The objective function and nonlinear constraints are reformulated as linear formulations, transforming the original nonlinear problem into a mixed integer linear programming model. The proposed approach was tested in 50 instances to compare the nonlinear and linear formulations. The results prove that the proposed linearization outperforms the nonlinear formulation achieving convergence to a better local optimum with shorter computational time. This method provides flexibility to the decisionmaker allowing the analysis of scenarios in a shorter time.
1. Introduction and Context of the Problem
Supply chain design is a critical strategy for achieving competitiveness in current global business environments. The supply chain consists of all functions involved, directly and indirectly, in meeting customer needs. Decisionmaking in supply chain management is classified into three hierarchical levels: strategic (longterm), tactical (mediumterm), and operational (shortterm). Frequently, these decisions are analyzed and solved independently at each planning stage, which generates an overall suboptimal solution when compared to solutions from comprehensive models. Therefore, supply chain management demands the use of new strategies and technologies to meet the current challenges of economic globalization [1], especially the highly dynamic behavior of customer demand. This characteristic imposes a greater challenge in developing and designing responsive and efficient supply chain networks. The mathematical formulation in the operational planning stage is highly critical and must model the relevant characteristics of the supply chain.
Facility location problems, which are typically used to design distribution networks, involve determining the sites where to install resources, as well as the assignment of potential customers to those resources. This family of problems typically assumes a linear cost function and a set of deterministic customer demand. These assumptions avoid modelling interactions between facility location and inventory management decisions. According to these facts, the inventorylocation research literature is aimed to integrate the demand uncertainty and the riskpooling effect into supply chain network design [2].
This paper presents a mixed integer nonlinear programming formulation aimed to optimize the distribution network design. The model includes decisions on plants and warehouses location, transportation and assignment between facilities (plantswarehouses, warehousesretailers), supplier selection, and implicit inventory decisions. The model is named a FourEchelon InventoryLocation Problem with Supplier Selection (4EILPSS). The mathematical formulation models a fourechelon supply chain, with a set of potential sites for allocation of plants and warehouses, as shown in Figure 1. Each warehouse and plant has a limited capacity and manages a continuous review inventory policy to satisfy the demand of retailers and warehouses, respectively.
The objective is to determine the flow of material in each echelon and the location of plants and warehouses. Inventory management decisions are simultaneously modeled with the design of the distribution network. The problem of locating facilities is commonly used as a base for the design of supply chains. However, the effect of risk integration by inventory management decisions is not directly modeled by the classic facility location problem (FLP). Additionally, supplier selection is a relevant issue because it affects many important aspects of the supply chain performance. In some real cases, for example in the context of the automotive industry, it is common to consider the involvement of suppliers in the supply chain planning stage. Supplier selection is a complex process that involves identifying potential suppliers, requesting information, establishing contract terms, negotiation, and evaluation. Several criteria are considered for supplier selection: quality, delivery time, price, service, and flexibility, among others. In the problem studied, we are assuming that the available suppliers were preselected considering several of the criteria described before, and the last decision was made based on the outcome of the optimization model that considers cost, lead time, and capacity. This is a common practice in industry, especially automotive industry, which is the motivation of this work, where an exclusive group of suppliers has the appropriate certifications by Original Equipment Manufacturers to participate in bids and requests.
Additionally, in the recent years, there has been an offshore moving of automotive manufacturing plants from industrialized countries to emerging economies. This happened in China, India, Mexico, Brazil, Thailand, and Turkey, among others [3]. This required the design of new supply chains to meet the demand, including new plants and warehouses.
The remainder of this paper is organized as follows. Section 2 reviews the relevant literature in supply chain design, with special emphasis on risk integration that results from the integration of an inventory management policy. Section 3 describes the features of the problem and presents two mathematical formulations: a nonlinear mixed integer programming model and a linearization of the aforementioned model using a piecewise linear function to obtain a linear mixed integer programming model. Section 4 presents the experimentation and results for some instances. Lastly, Section 5 presents the conclusions and final remarks.
2. Literature Review
Inventorylocation problems (ILPs) have been extensively studied over the past two decades. This section presents a review of the main contributions reported in the literature, showing trends in formulations and solution approaches. Although ILPs are extensions of Facility Location Problems (FLPs), this paper does not discuss the related FLP literature. Interested readers may consult Hamacher and Drezner [4], Daskin [5], Nickel and Puerto [6], Melo et al. [7], and Eiselt and Marianov [8, 9].
Jayaraman [10] studied a basic ILP, based on a mixed integer programming model for locating warehouses and satisfying customer’s demand for different products. The inventory modeling included the EOQ (Economic Order Quantity) model with fixed lot size and deterministic demand. Nozick and Turnquist [11] investigated an ILP, which included fixed warehouse location costs, and transportation and inventory costs. The problem was solved using a combination of a greedy and improvement algorithm. Erlebacher and Meller [12] presented an ILP with stochastic demands for a twostage supply chain network, along with some heuristic procedures to solve a variety of instances. Another twoechelon ILP was studied by Nozick and Turnquist [13]. The model minimizes the inventory and unfulfilled demand costs. They used a combination of iterative greedy and improvement algorithms as a solution strategy. Nozick and Turnquist [14] analyzed a distribution system and the tradeoff between facility costs, inventory costs, transportation costs, and customer responsiveness. They built an efficient frontier for assessing the different solution configurations.
Daskin et al. [15], Shen et al. [16], and Miranda and Garrido [17] presented similar versions of ILPs, proposing mixed nonlinear integerprogramming formulations. The model considers a continuous control inventory policy with stochastic demand, the safety stock cost is determined by a chanceconstrained approach, and the ordering decisions are based on the EOQ model with order quantities as decision variables. Daskin et al. [15] and Miranda and Garrido [17] solved the model using Lagrangian relaxation, while Shen et al. [16] reformulated the model as a setcovering problem, and developed a columngeneration algorithm. Shu et al. [18] studied the column generation approach proposed by Shen et al. [16], showing that the pricing problem gives rise to a new class of the submodular function minimization problem. They discussed how the general pricing problem could be solved by exploiting certain special structures.
Snyder et al. [19] presented a stochastic programmingbased ILP including some random parameters defined by discrete scenarios. The location model explicitly handles the economies of scale and riskpooling effects that result from consolidating inventory sites. They used a Lagrangian relaxationbased algorithm. Shu et al. [20] proposed a twostage stochastic model to address an integrated location, and twoechelon, inventory network design under uncertainty. They defined the problem as a twostage nonlinear discrete optimization problem. The first stage decides which warehouses to open and the second one decides the warehouseretailer assignments and the twoechelon inventory replenishment. They used a set covering formulation transformed into a pricing problem, which was solved using column generation.
ILPs with inventory capacity constraints were analyzed by Miranda and Garrido [21, 22] and Ozsen et al. [23]. Ozsen et al. [23] introduced a 100% service level constraint for inventory capacity to a previous ILP presented by Daskin et al. [15], while Miranda and Garrido [21, 22] discussed how the inventory capacity constraint can be controlled by a previously defined service level (usually less than 100%), based on a chanceconstrained formulation. Ozsen et al. [23] and Miranda and Garrido [21, 22] developed a Lagrangian relaxationbased solution approach. Miranda and Cabrera [24] and Cabrera et al. [25] introduced an ILP with periodic inventory control review, where inventory capacity constraints are modeled in a stochastic manner.
Notice that all previous works are focused on optimizing location and inventory costs and decisions are taken at a single stage of the supply chain (i.e., warehouses), along with transport or assignment decisions (plantwarehouses and/or warehousesretailers). In the works described forward, some additional considerations of costs and operations management decisions are integrated into similar ILP formulations.
Kang and Kim [26] investigated an ILP for warehouse location, integrating inventory costs and assignment decisions at the warehouse and retailers. They formulated a nonlinear mixed integer programming model and developed a Lagrangian relaxationbased heuristic.
Silva and Gao [27] introduced a joint replenishment ILP. They proposed a Greedy Randomized Adaptive Search Procedure to solve the problem. Tancrez et al. [28] studied an ILP in threelevel supply chain networks. The decision set included distribution centers location, flows allocation, and shipment sizes. They proposed a nonlinear formulation that decomposes into a closedform equation and a linear programming model when the distribution center flows are fixed. They developed an iterative heuristic that estimates the distribution center flows a priori, solves the linear program, and then improves the distribution center flows estimations. Similarly, Diabat et al. [29] studied a simplified multiechelon ILP (with location decisions at a single level), where the problem was formulated as a nonlinear mixed integer program and was solved using a Lagrangian relaxation approach.
Shahabi et al. [30] developed a model to coordinate facility location and inventory control in a fourechelon supply chain with hubs, which helps in reducing transportation costs by consolidating shipments. There were three decisions: warehouses and hubs locations, assignment of suppliers and retailers, and inventory control decisions and costs at a single stage (warehouses). A mixed integer nonlinear programming formulation was first presented and then transformed into a compact conic mixed integer programming formulation. Commercial solvers were used to solve the problem.
Shen and Daskin [31], Gaur and Ravindran [32], Miranda and Garrido [33], Mak and Shen [34], and Escalona et al. [35] are a variety of works in which customer or system service levels are integrated and analyzed explicitly based on previous ILP formulations. Genetic Algorithms, Lagrangian relaxation, and twostep greedybased heuristics are some of the approaches employed to solve the analyzed ILPs.
Atamtürk et al. [36] introduced a conic integer programming approach to reformulate a variety of ILPs. These reformulations, for different cases of nonlinearities, can be solved directly using standard optimization software packages without the need for designing specific algorithms. They concluded that their approach led to similar or better computational times than other approaches reported in the literature, and they claimed the applicability of their approach to model more general, related problems.
Kaya and Urek [37] proposed a model for a closed loop supply chain that integrates location, inventory and pricing decisions. The location decision is applied to the level of collectiondistribution centers. Inventory is modeled without uncertainty in demand but affected by a pricedemand function. The mixed integer nonlinear programming model is solved using commercial optimization software for small instances, and several heuristics for large instances. The heuristics include a Piecewise linearization, and hybrids of Simulated Annealing, Tabu Search, Genetic Algorithms, and Variable Neighborhood Search.
Schuster Puga and Tancrez [38] analyze a supply chain network with three levels. The flow of product is from a central plant to distribution centers to retailers. Location decisions are applied to the distribution centers. The model developed considers inventory decisions, and transportation costs follow a full truckload approach. They developed an iterative heuristic to solve the continuous nonlinear programming model. Also, they developed a conic quadratic mixedinteger program to transform the original formulation. They solved instances with up to 1000 retailers and 1000 distribution centers.
Escalona et al. [39] studied a supply chain network with three levels and inventorylocation decisions. They add a stratification of customers according to classes based on the service level required. The mixed integer nonlinear programming model is transformed into a conic quadratic mixedinteger problem. Their main objective is to analyze the effect of different inventory policies on the configuration of the supply chain network.
TapiaUbeda et al. [40] studied an ILP for a supply chain consisting of two levels, with warehouses shipping product to customers. They developed a mixed integer nonlinear programming model to be solved using a Benders Decomposition algorithm. The nonlinear nature of the original model is transferred to the subproblems, but for the instances solved the subproblems can be solved to optimality without a big effort.
It must be noticed that most of the previous ILP formulations integrate location decisions at one or two stages (warehouse location, plant location and supplier selection, warehouse and consolidating center location, etc.). Inventory costs and decisions are modeled at one or two levels (warehouses, warehouses, and plants, warehouses and retailers, etc.), but none of these previous work simultaneously integrates location and inventory management decisions at two stages in the supply chain. Therefore, this research paper fills an important gap in the field of supply chain network design.
To the best of our knowledge, only a previous work by ourselves has considered the design of such a complex supply chain as in this research, with some differences in the model and more important, in the solution method. Perez Loaiza et al. [41] studied a supply chain of four levels with location and inventory decisions and a biobjective approach. The mixed integer nonlinear programming model is solved using an evolutionary algorithm. Pareto fronts are compared with respect to solutions obtained with commercial optimization software.
The following part of the literature review discusses linearization strategies for nonlinear problems. You and Grossmann [42] addressed an ILP from the chemical industry, the design of a multiechelon supply chain and inventory system in the presence of uncertain customer demands. A mixed integer nonlinear program modeled the transportation, inventory, and network structure of a multiechelon supply chain. The model had a nonconvex objective function. They reformulated the problem as a separable concave program. A spatial decomposition algorithm based on the integration of Lagrangian relaxation and piecewise linear approximation was proposed to solve the model. The transformation assures all the constraints are linear and the only nonlinear terms are univariate concave terms in the objective function. Petridis [43] addressed a multiproduct, multiechelon supply chain network with demand uncertainty. Decisions about the selection of facilities and their capacity are made. Information about the flows of products transferred and the safety stock at each distribution center was derived. The lead time of an order to a customer is computed, using the probabilities of overstocking and understocking. The problem was formulated as a single period mixed integer nonlinear programming problem. Linearization techniques for selected nonlinear terms of the models were explored in order to reduce the computational effort for solving the model. The linearization was achieved by rewriting one of the nonlinear constraints that had a product of a continuous and a binary variable. The objective function was not linearized.
According to the discussed literature review, this paper contributes with an inventorylocation model, under continuous review policy, to optimize a fourechelon supply chain network, encompassing multiple suppliers, plants, warehouses, and retailers. The problem is formulated as a mixed integer nonlinear programming model, considering a single item, with stochastic demands across de supply chain network. Decision variables include warehouse and plant location, plantwarehouse assignments, retailer assignment to warehouses, and shipments from suppliers to plants. This ILP integrates inventory costs and decisions at plants and warehouses. As it can be observed from the related literature, the proposed network design model of a fourechelon supply chain has a nonlinear formulation that arises when inventory aspects are integrated into facility location problems, yielding to a nonlinear mixed integer programming model.
Nonlinear mixed integer programming models are extremely hard to solve to optimality, particularly when commercial optimization solvers are employed. They combine the mixed integer programming (NPHard) nature with nonlinear and nonconvex components. In many cases, researchers and practitioners face two major challenges: local optimum solutions and prohibitive computational times for tackling practical problems. Moreover, these challenges increase when more stages and decision variables are added to the model, such as in this research article. In order to reduce the computational effort, and as a strategy to develop a heuristic framework to tackle this complex problem, a linearization of the model is proposed based on a piecewise linear approximation of the objective function, and a linear reformulation of nonlinear constraints, yielding to a mixed integer linear programming model. The approach is developed following the approach proposed by Diabat and Theodorou [44], who applied this technique in a smaller supply chain network with less supply chain stages and decisions.
3. Problem Description and Modeling
3.1. Continuous Review Based Inventory Control Policy
A continuous review based inventory control policy is assumed at plants and warehouses with the aim of dealing with uncertain and stochastic demands (warehouses deal with retailer demands and plants deal with warehouse demands or orders). Figure 2 shows the evolution of the physical inventory level when this inventory control policy is assumed. When the inventory level reaches reorder point , an order of units of the product is requested. The units of product are received after a fixed and known lead time LT. In other words, RP is the critical inventory level that generates a new order. In the case of warehouses, they generate an order to plants, while in the case of plants they generate a production order or batch. In both cases, a widely accepted approach states that the reorder point at each location, RP, must be set in order to ensure that served demand arising during lead time LT will not be greater than this reorder point, at least with a fixed and known probability. This probability is also known as inventory service level.
The existing inventory level just before the order arrives at each location is known as the safety stock, while the inventory level that is observed over the safety stock is known as cyclic or working inventory.
Warehouse demand depends on customers assigned to each warehouse , and plant demand depends on warehouses assigned to each plant .
To calculate the working inventory cost (WIC), the total inventory cost is estimated bywhere UC is the purchasing cost, d is the demand, q is the order quantity, OC is the ordering cost, and HC is the inventory holding cost. The first derivative of (1) with respect to is obtained, equal to zero, and solved for . The optimal order quantity is obtained as Then, the working inventory costs WIC, the last component of (1), is obtained using from (2) inAssuming an uncertain demand with normal distribution, the reorder point RP can be calculated according tosuch thatwhere ZT represents the value of the Standard Normal Distribution (0,1) that accumulates a probability of 1 (service level). This value of ZT in (5) determines the safety stock level at each facility. The variance of the demand is given by u, and LT represents the lead time.
RP is the amount of inventory that should be available in each plant and warehouse to absorb fluctuations in demand from retailers during the lead time. Safety stock is the carried out inventory to prevent stock outs. It describes the balance between the cost of holding inventory and profit foregone as a result of outages. In this case, it was not considered the cost of not having stock, only the safety stock cost (SSC) in each plant or warehouse. SSC is given by A similar explanation of these terms can be found in Daskin et al. [15] and Shen et al. [16].
3.2. Description of the Model Formulation
In this section, the mathematical formulation of the problem studied in this paper is presented. Section 3.2.1 discusses the mixed integer nonlinear model for distribution network design, 4EILPSS. The linearization process for the constraints is explained in this section too. Section 3.2.2 illustrates the linearization process of both safety stock and working inventory costs using an approximation based on piecewise linear functions; thereby a mixed integer linear model is obtained (4EILPSSL). The motivation to linearize the original model is based on the characteristics of commercial optimization software. In most of the cases, the standard methods used to solve nonlinear models do not guarantee to find a global optimum solution, and when this is possible the computational time may be too long and unreasonable.
3.2.1. Mixed Integer Nonlinear Model
Notation and Definitions
Sets : Set of suppliers, . : Set of potential plant locations, . : Set of potential warehouse locations, . : Set of retailers, .
Parameters : Production capacity of supplier . : Maximum capacity of warehouse . : Maximum capacity of plant . : Average demand for each retailer . : Fixed cost of opening a plant at site . : Fixed cost of opening a warehouse at site . : Inventory holding cost per unit in plant . : Inventory holding cost per unit in warehouse . : Deterministic delivery lead time to plant . : Deterministic delivery lead time to warehouse . : Ordering cost at plant . : Ordering cost at warehouse . : Transportation cost per unit from supplier to plant . : Transportation cost per unit from plant to warehouse . : Transportation cost per unit from warehouse to retailer . : Purchase cost per unit at supplier . : Demand variance at retailer . : Values of ZT for a given service level, assuming a Normal Distribution. : Working inventory cost at plant . : Safety stock at plant . : Working inventory cost at warehouse . : Safety stock at warehouse .
Decision Variables : 1 if warehouse is opened; 0 otherwise. : 1 if plant is opened; 0 otherwise. : 1 if plant serves warehouse ; 0 otherwise. : 1 if warehouse serves retailer ; 0 otherwise. : Number of units of product produced by supplier and shipped to plant .
Auxiliary Variables : Average demand assigned to plant . : Demand variance assigned to plant . : Average demand assigned to warehouse . : Demand variance assigned to warehouse . : Product of ordering cost, inventory holding cost and demand assigned to plant . : Product of lead time and demand variance assigned to plant : Product of ordering cost, inventory holding cost and demand assigned to warehouse . : Product of lead time and demand variance assigned to warehouse .
Thus, the mathematical model for the problem is formulated as The objective function of the optimization model, 4EILPSS, is to minimize the total cost given by the cost components in (7). These components represent the fixed cost if plant and warehouse are open, the cost of safety stock and working inventory at each plant and warehouse, the purchase and transportation cost to serve the average demand shipped from supplier to plant , and the transportation cost from plant to warehouse and from there to retailer . Equation (8) ensures that each retailer is served exactly by one warehouse. Constraints in (9) guarantee that each opened warehouse is served only by one plant. These constraints follow the single source assumption used frequently in the literature and practice. Equation (10) indicates that suppliers should ship at least one unit to open plants. Equations (11) and (12) compute the transferred demand average and variance from retailers to warehouses. Similarly, nonlinear constraints in (13) and (14) compute the demand average and its variance transferred from warehouses to plants. Equation (15) meets the demand from suppliers to plants. Equations (16) to (19) calculate auxiliary variables used into the objective function. Their use is not essential at this point, but these equations are helpful later to address the nonlinearity of the model. Inequalities in (20), (21), and (22) ensure that the capacities of the warehouses, plants, and suppliers are not exceeded, respectively, for open facilities. Finally, (23) and (24) show the nature of the variables.
A reformulation of the constraints in (13) and (14) is presented forward. In these constraints, the product of an integer auxiliary variable and a binary decision variable occurs, for instance and in (13). Combining (13) and (11), the first can be reformulated asNonlinear constraints are produced by the product of two binary decision variables, i.e., . The nonlinearity in (25) is removed when are added, with the introduction of a binary variable , such that if and , then variable , and otherwise may take the value of zero [45].
If is substituted by in (25), the result will be the linear constraintsinstead of those in (13). The constraints in (29) involve the product of a parameter and a binary decision variable . Using the same procedure, constraints in (14) are reformulated and presented asOnce a solution for the model is obtained, the order quantity and reorder point at each open facility can be calculated. Using the results for variables as demands for plants, with their corresponding values for the ordering cost and holding cost, the order quantity at each plant can be calculated using (2). Using the results for variables as demands and as demand variances for plants, with their corresponding values for the lead time and service level factor, the reorder point at each plant can be calculated using (5). The same can be done for warehouses using the values of as demands and as demand variances.
3.2.2. Linear Approximation of the Objective Function
This section presents a linear approximation approach for the safety stock costs SSC and working inventory costs WIC. WIC in plant is calculated as , as derived from (3) and (16). The approximation is based on a piecewise linear function (PLF) defined by two or more straightlines, presented in different sections of the domain of the function as shown in Figure 3.
Considering a partition of in sections, the function is approximated linearly by the expression , where and , for sections. Here, is the value of at the upper bound of section , is the value of at the upper bound of section , is the slope along the section , is equal to 1 if section of the piecewise linear function of is used or 0 otherwise. Also, are the units of in any point along section .
The linear approximation of WIC at each plant is given by Equation (31) represents the linear approximation of the WIC at plant . Equation (32) ensures that the average demand in plant is considered for the linear approximation. Equation (33) establishes that only a section must be used in the linear approximation. Equation (34) guarantees that the linear approximation is carried out only in open plants. Finally, (35) describes the domain of the variables. The same analysis is considered to generate the linear approximation of SSC in each plant and WIC and SSC in each warehouse.
An alternative objective function Z1’ is defined. Z1’ is equivalent to Z1 but considering the linear approximation of WIC and SSC introduced before. The new optimization model, with objective function Z1’ and all the linearized constraints, is named 4EILPSSL.
4. Results and Discussion
In this section, five instance sizes with 10 cases each one are presented. The code of the instance size indicates the number of suppliers followed by the number of potential plants, followed by the number of potential warehouses, and the number of retailers. For example, instance size 53510 indicates 5 suppliers, 3 potential plants, 5 potential warehouses, and 10 retailers.
A commercial license of LINGO 14 was used to solve every case, with the nonlinear and linear models. The nonlinear model corresponds to the formulation described by (7) to (24). The linear model corresponds to the linearization of the objective function and the constraints applied to the original nonlinear model. A service level of 95% was considered for all the cases. Tables 1–6 show the parameters values of the base case, which were varied to construct the different cases. The case base has 10 potential suppliers, 3 potential plants, 5 potential warehouses, and 25 potential retailers. For example, to generate one instance of the size 53510, 5 out of the 10 available suppliers in the base case were selected randomly, the 3 plants are used, the 5 warehouses are used, and 10 out of the 25 retailers were selected randomly from the base case. To generate one instance of the size 73520, 7 out of the 10 available suppliers in the base case were selected randomly, the 3 plants are used, the 5 warehouses are used, and 20 out of the 25 retailers were selected randomly from the base case.






The uncertainty in demand is represented by the variance values shown in Table 3. In all the instances of the linear model, four intervals were used for the linearization, following the results from Diabat and Theodorou [44].
Table 7 shows the computational results of applying the nonlinear and the linear model with a service level of 95% (ZT_{(0.95)} = 1.64). Table 7 shows the objective function value Z1, the saved cost from the comparison, the % saving with respect to the solution of the nonlinear model, the run elapsed time for the solution with the linear model, and the cases where the solution was obtained only with the linear model. The comparison denotes the significant benefit (% Saving) that can be obtained by using the linear model instead of the nonlinear model. In this study, the average elapsed run time from the solutions of the linear model was set as a time limit on each group of instances when solving the nonlinear model. In the last column in Table 7, the cases where the nonlinear model does not converge (in the time limit) are shown as well. In order to make a fair comparison between both models, the total cost Z2 was obtained by using the values of the solution obtained by the linear model that optimizes Z1’ but evaluated in the original objective function Z1. Let X’ be the optimal solution of 4EILPSSL, then Z2 = Z1 (X’).

Table 8 shows the structure of the supply chain according to the results obtained with the models. In many cases, the structures are very similar, but the assignments between customers and warehouses, warehouses and plants, and flows from suppliers produced the differences in costs.

Figure 4 shows the comparison of the average total cost of the nonlinear and the linear models, where a solution was found. In all cases, the solutions of the linear model are cheaper than the solutions from the nonlinear model. The range of the average percent savings (% saving) is from 0.51% to 5.11%. The results of the nonlinear model may be better using a more efficient nonlinear optimizer, especially for the small instances, but it is expected that the linear model will be more efficient for larger instances.
The comparison of the number of solutions obtained between the linear and nonlinear model is presented in Figure 5. The nonlinear model does not converge in all cases with a similar run time than the linear model. The range of performance of the nonlinear model goes from 50% to 100% among the instances. Using the linear model an optimal solution was obtained in all cases. For these instance sizes, one can be confident in obtaining a good solution with the linear model, which is not the case with the nonlinear model.
5. Conclusions and Final Remarks
This research introduced an inventorylocation model, 4EILPSS, for the network design of a fourechelon supply chain, considering the location of warehouses and plants, transportation cost, inventory costs, and supplier selection. The inventory policy and network design decisions are tactical and strategic decisions that should be analyzed concurrently. This is attained using the aforementioned model. This model can be applied to design complex and long supply chains, like in the automotive industry because of the offshore movement of production facilities from industrialized countries to emerging economies. In many cases, there is infrastructure available in the country, and it is easy to adapt the model to that situation fixing locations for the available facilities and allowing the model to decide for the new optimum flows and locations of new facilities.
The model proposed becomes nonlinear because of the application of the continuous review inventory policy. This feature added to the computational complexity inherited from the Facility Location Problem makes the problem hard to solve. The mixed integer nonlinear model was reformulated using a piecewise linear function and a reformulation of nonlinear constraints in order to generate a mixed integer linear model. Several instances were solved using the linear and nonlinear models, comparing the results in terms of cost savings, number of solutions obtained, and computational runtime. It was observed that solving the nonlinear model with the commercial software generated local optimal solutions, and in several cases the software either did not converge or resulted in infeasible solutions. The total cost of the original objective function is lower, in most of the cases, with the solution obtained by the linear model than the solution obtained by the nonlinear model.
The main contribution of the paper is twofold. Firstly, a network optimization model for a particular four level supply chain, where facility location and inventory management decisions are considered in two stages (plants and warehouses), and transportation or assignment decisions are modeled in three stages of the supply chain, i.e., suppliersplants, plantswarehouses, and warehousescustomers. Secondly, a linearization approach of the proposed nonlinear model, in order to facilitate its applicability and its effective and efficient implementation, improves both solutions quality and computational times.
The proposed linear model proves to be useful for decision makers interested in analyzing and designing supply chain networks with the structure of a 4EILPSS. The problem can be analyzed and solved in a short period of time and with significant confidence in the solution quality.
The experiments conducted allowed us to understand the complexity of the problem. For this study, LINGO was used to solve the linear and the nonlinear models, but future work may take advantage of using more efficient commercial software to solve larger instances. As observed, using commercial optimizers may not be the best alternative when it comes to solve very large instances. Therefore, immediate future research should consider the development of other heuristics aimed to tackle large instances (Lagrangian relaxation, Bender decompositions, Cutting planes, and Metaheuristics). Another direction is to extend the pricing approach proposed by Shu et al. [18] to solve a similar problem to the long supply chain proposed in this paper.
In addition, a deeper study may involve the analysis of the relations between parameters like costs, capacities, lead times, and demand variability, to understand their impact in the supply chain network configuration.
Finally, it would be interesting to extend the model to analyze multiobjective situations, more common when designing supply chain networks for real situations.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
Acknowledgments
This research was supported by Programa para el Mejoramiento del Profesorado PROMEP (Grant ITESUR001 Charter no. PROMEP/103.5/10/5709) and Universidad Panamericana (Grant no. UPCI2018INGGDL06).
References
 M. Mourits and J. J. Evers, “Distribution network design: an integrated planning support framework,” Logistics Information Management, vol. 9, pp. 45–54, 1995. View at: Google Scholar
 R. Z. Farahani, H. Rashidi Bajgan, B. Fahimnia, and M. Kaviani, “Locationinventory problem in supply chains: A modelling review,” International Journal of Production Research, vol. 53, no. 12, pp. 3769–3788, 2015. View at: Publisher Site  Google Scholar
 OICA, International Organization of Motor Vehicle Manufacturers, Production statistics, 2017, http://www.oica.net/category/productionstatistics/2017statistics/.
 H. Hamacher and Z. Drezner, Facility Location: Applications and Theory, Springer, Berlin, Germany, 2002.
 M. S. Daskin, Network and Discrete Location: Models, Algorithms, and Applications, WileyInterscience, New York, NY, USA, 1st edition, 1995. View at: Publisher Site  MathSciNet
 S. Nickel and J. Puerto, Location Theory: A Unified Approach, Springer, New York, NY, USA, 2005.
 M. T. Melo, S. Nickel, and F. SaldanhadaGama, “Facility location and supply chain management—a review,” European Journal of Operational Research, vol. 196, no. 2, pp. 401–412, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 H. A. Eiselt and V. Marianov, Foundations of Location Analysis, Springer, New York, NY, USA, 2011.
 H. A. Eiselt and V. Marianov, Applications of Location Analysis, Springer, New York, NY, USA, 2015.
 V. Jayaraman, “Transportation, facility location and inventory issues in distribution network design: An investigation,” International Journal of Operations and Production Management, vol. 18, no. 5, pp. 471–494, 1998. View at: Publisher Site  Google Scholar
 L. K. Nozick and M. A. Turnquist, “Integrating inventory impacts into a fixedcharge model for locating distribution centers,” Transportation Research Part E: Logistics and Transportation Review, vol. 34, no. 3, pp. 173–186, 1998. View at: Publisher Site  Google Scholar
 S. J. Erlebacher and R. D. Meller, “The interaction of location and inventory in designing distribution systems,” IIE Transactions, vol. 32, no. 2, pp. 155–166, 2000. View at: Google Scholar
 L. K. Nozick and M. A. Turnquist, “A twoechelon inventory allocation and distribution center location analysis,” Transportation Research Part E: Logistics and Transportation Review, vol. 37, no. 6, pp. 425–441, 2001. View at: Publisher Site  Google Scholar
 L. K. Nozick and M. A. Turnquist, “Inventory, transportation, service quality and the location of distribution centers,” European Journal of Operational Research, vol. 129, no. 2, pp. 362–371, 2001. View at: Publisher Site  Google Scholar
 M. S. Daskin, C. R. Coullard, and Z. J. M. Shen, “An inventorylocation model: formulation, solution algorithm and computational results,” Annals of Operations Research, vol. 110, pp. 83–106, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 Z.J. M. Shen, C. R. Coullard, and M. S. Daskin, “A joint locationinventory model,” Transportation Science, vol. 37, no. 1, pp. 40–55, 2003. View at: Publisher Site  Google Scholar
 P. A. Miranda and R. A. Garrido, “Incorporating inventory control decisions into a strategic distribution network design model with stochastic demand,” Transportation Research Part E: Logistics and Transportation Review, vol. 40, no. 3, pp. 183–207, 2004. View at: Publisher Site  Google Scholar
 J. Shu, C.P. Teo, and Z. M. Shen, “Stochastic transportationinventory network design problem,” Operations Research, vol. 53, no. 1, pp. 48–60, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 L. V. Snyder, M. S. Daskin, and C.P. Teo, “The stochastic location model with risk pooling,” European Journal of Operational Research, vol. 179, no. 3, pp. 1221–1238, 2007. View at: Publisher Site  Google Scholar
 J. Shu, Q. Ma, and S. Li, “Integrated location and twoechelon inventory network design under uncertainty,” Annals of Operations Research, vol. 181, pp. 233–247, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 P. A. Miranda and R. A. Garrido, “A simultaneous inventory control and facility location model with stochastic capacity constraints,” Networks and Spatial Economics, vol. 6, no. 1, pp. 39–53, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 P. A. Miranda and R. A. Garrido, “Valid inequalities for Lagrangian relaxation in an inventory location problem with stochastic capacity,” Transportation Research Part E: Logistics and Transportation Review, vol. 44, no. 1, pp. 47–65, 2008. View at: Publisher Site  Google Scholar
 L. Ozsen, C. R. Coullard, and M. S. Daskin, “Capacitated warehouse location model with risk pooling,” Naval Research Logistics (NRL), vol. 55, no. 4, pp. 295–312, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 P. A. Miranda and G. Cabrera, “Inventory location problem with stochastic capacity constraints under periodic review (R, s, S),” in Proceedings of the International Conference on Industrial Logistics: Logistics and Sustainability, (ICIL '10), pp. 289–296, 2010. View at: Google Scholar
 G. Cabrera, P. A. Miranda, E. Cabrera et al., “Solving a novel inventory location model with stochastic constraints and inventory control policy,” Mathematical Problems in Engineering, vol. 2013, Article ID 670528, 12 pages, 2013. View at: Publisher Site  Google Scholar
 J.H. Kang and Y.D. Kim, “Inventory control in a twolevel supply chain with risk pooling effect,” International Journal of Production Economics, vol. 135, no. 1, pp. 116–124, 2012. View at: Publisher Site  Google Scholar
 F. Silva and L. Gao, “A joint replenishment inventorylocation model,” Networks and Spatial Economics, vol. 13, no. 1, pp. 107–122, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 J.S. Tancrez, J.C. Lange, and P. Semal, “A locationinventory model for large threelevel supply chains,” Transportation Research Part E: Logistics and Transportation Review, vol. 48, no. 2, pp. 485–502, 2012. View at: Publisher Site  Google Scholar
 A. Diabat, J.P. Richard, and C. . Codrington, “A Lagrangian relaxation approach to simultaneous strategic and tactical planning in supply chain design,” Annals of Operations Research, vol. 203, pp. 55–80, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 M. Shahabi, S. Akbarinasaji, A. Unnikrishnan, and R. James, “Integrated inventory control and facility location decisions in a multiechelon supply chain network with hubs,” Networks and Spatial Economics, vol. 13, pp. 497–514, 2013. View at: Publisher Site  Google Scholar
 Z.J. M. Shen and M. S. Daskin, “Tradeoffs between customer service and cost in integrated supply chain design,” Manufacturing and Service Operations Management, vol. 7, no. 3, pp. 188–207, 2005. View at: Publisher Site  Google Scholar
 S. Gaur and A. R. Ravindran, “A bicriteria model for the inventory aggregation problem under risk pooling,” Computers & Industrial Engineering, vol. 51, no. 3, pp. 482–501, 2006. View at: Publisher Site  Google Scholar
 P. A. Miranda and R. A. Garrido, “Inventory servicelevel optimization within distribution network design problem,” International Journal of Production Economics, vol. 122, no. 1, pp. 276–285, 2009. View at: Publisher Site  Google Scholar
 H.Y. Mak and Z. M. Shen, “A twoechelon inventorylocation problem with service considerations,” Naval Research Logistics (NRL), vol. 56, no. 8, pp. 730–744, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 P. Escalona, F. Ordóñez, and V. Marianov, “Joint locationinventory problem with differentiated service levels using critical level policy,” Transportation Research Part E: Logistics and Transportation Review, vol. 83, pp. 141–157, 2015. View at: Publisher Site  Google Scholar
 A. Atamtürk, G. Berenguer, and Z.J. Shen, “A conic integer programming approach to stochastic joint locationinventory problems,” Operations Research, vol. 60, no. 2, pp. 366–381, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 O. Kaya and B. Urek, “A mixed integer nonlinear programming model and heuristic solutions for location, inventory and pricing decisions in a closed loop supply chain,” Computers & Operations Research, vol. 65, pp. 93–103, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 M. a. Schuster Puga and J.S. Tancrez, “A heuristic algorithm for solving large locationinventory problems with demand uncertainty,” European Journal of Operational Research, vol. 259, no. 2, pp. 413–423, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 P. Escalona, V. Marianov, F. Ordóñez, and R. Stegmaier, “On the effect of inventory policies on distribution network design with several demand classes,” Transportation Research Part E: Logistics and Transportation Review, vol. 111, pp. 229–240, 2018. View at: Publisher Site  Google Scholar
 F. J. TapiaUbeda, P. A. Miranda, and M. Macchi, “A generalized Benders decomposition based algorithm for an inventory location problem with stochastic inventory capacity constraints,” European Journal of Operational Research, vol. 267, no. 3, pp. 806–817, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 R. E. Perez Loaiza, E. OlivaresBenitez, P. A. Miranda Gonzalez, A. Guerrero Campanur, and J. L. Martinez Flores, “Supply chain network design with efficiency, location, and inventory policy using a multiobjective evolutionary algorithm,” International Transactions in Operational Research, vol. 24, no. 12, pp. 251–275, 2017. View at: Publisher Site  Google Scholar
 F. You and I. E. Grossmann, “Integrated multiechelon supply chain design with inventories under uncertainty: MINLP models, computational strategies,” AIChE Journal, vol. 56, no. 2, pp. 419–440, 2010. View at: Google Scholar
 K. Petridis, “Optimal design of multiechelon supply chain networks under normally distributed demand,” Annals of Operations Research, vol. 227, no. 1, pp. 63–91, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 A. Diabat and E. Theodorou, “A locationinventory supply chain problem: Reformulation and piecewise linearization,” Computers & Industrial Engineering, vol. 90, pp. 381–389, 2015. View at: Publisher Site  Google Scholar
 A. Nyberg, I. E. Grossmann, and T. Westerlund, “An efficient reformulation of the multiechelon stochastic inventory system with uncertain demands,” AIChE Journal, vol. 59, no. 1, pp. 23–28, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Aaron Guerrero Campanur et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.