Design of a Logistics Nonlinear System for a Complex, Multiechelon, Supply Chain Network with Uncertain Demands

1 ITSUruapan, Tecnologico Nacional de Mexico, Carr. Uruapan-Carapan 5555, Uruapan, Michoacan 60015, Mexico 2Facultad de Ingenieria, Universidad Panamericana, Prolongacion Calzada Circunvalacion Poniente 49, Zapopan, Jalisco 45010, Mexico 3Department of Engineering Sciences, Universidad Andres Bello, Quillota 980, Viña del Mar, 2531015, Chile 4ITApizaco, Tecnologico Nacional de Mexico, Av. Instituto Tecnológico S/N, Apizaco, Tlaxcala 90300, Mexico 5College of Business Administration, University of Texas at El Paso, 500 W. University Avenue, El Paso, TX 79968, USA


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. Decision-making in supply chain management is classified into three hierarchical levels: strategic (long-term), tactical (medium-term), and operational (short-term). 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 inventory-location research literature is aimed to integrate the demand uncertainty and the risk-pooling 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 (plants-warehouses, warehouses-retailers), supplier selection, and implicit inventory decisions. The model is named a Four-Echelon Inventory-Location Problem with Supplier Selection (4EILP-SS). The mathematical formulation models a four-echelon 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.

Literature Review
Inventory-location 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 two-stage supply chain network, along with some heuristic procedures to solve a variety of instances. Another two-echelon 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 trade-off 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 integer-programming formulations. The model considers a continuous control inventory policy with stochastic demand, the safety stock cost is determined by a chance-constrained 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 set-covering problem, and developed a column-generation 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 risk-pooling effects that result from consolidating inventory sites. They used a Lagrangian relaxation-based algorithm. Shu et al. [20] proposed a two-stage stochastic model to address an integrated location, and two-echelon, 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 warehouse-retailer assignments and the two-echelon 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 chance-constrained formulation. Ozsen et al. [23] and Miranda and Garrido [21,22] developed a Lagrangian relaxation-based 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 (plant-warehouses and/or warehouses-retailers). 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 relaxation-based 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 three-level supply chain networks. The decision set included distribution centers location, flows allocation, and shipment sizes. They proposed a nonlinear formulation that decomposes into a closed-form 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 four-echelon 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 two-step greedy-based 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 collection-distribution 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 inventory-location 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 mixed-integer problem. Their main objective is to analyze the effect of different inventory policies on the configuration of the supply chain network.
Tapia-Ubeda 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 bi-objective 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 non-convex 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 inventory-location model, under continuous review policy, to optimize a four-echelon 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, plant-warehouse 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 four-echelon 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 (NP-Hard) nature with nonlinear and non-convex 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.

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 by where 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 * = √ 2 ⋅ ⋅ .
Then, the working inventory costs WIC, the last component of (1), is obtained using * from (2) in Assuming an uncertain demand with normal distribution, the reorder point RP can be calculated according to such that where 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].

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, 4EILP-SS. 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 (4EILP-SS-L). 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. : Ordering cost at plant .

Mixed Integer Nonlinear
: 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 .
The objective function of the optimization model, 4EILP-SS, 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 8 Complexity 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 as Nonlinear constraints are produced by the product of two binary decision variables, i.e., ⋅ . The nonlinearity in (25) are added, with the introduction of a binary variable V , such that if = 1 and = 1, then variable V = 1, and otherwise V may take the value of zero [45].
If ⋅ is substituted by V in (25), the result will be the linear constraints instead of those in (13). The constraints in (29) involve the product of a parameter and a binary decision variable V . Using the same procedure, constraints in (14) are reformulated and presented as 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.

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 straight-lines, 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 −1 ⋅ + ⋅ , where 1 < 2 < 3 . . . < and 1 > 2 > 3 . . . > , for = 1, . . . 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, = − −1 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 4EILP-SS-L.

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 5-3-5-10 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    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 5-3-5-10, 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 7-3-5-20, 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 4EILP-SS-L, 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.

Conclusions and Final Remarks
This research introduced an inventory-location model, 4EILP-SS, for the network design of a four-echelon 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 two-fold. 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, plants-warehouses, and warehouses-customers. 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 4EILP-SS. 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.