#### Abstract

We review the problem for the design of supply chains called Capacitated Fixed Cost Facility Location Problem with Transportation Choices (CFCLP-TC). The problem is based on a production network of two echelons with multiple plants, a set of potential distribution centers, and customers. The problem is formulated as an optimization model with two objective functions based on time and cost. This paper proposes three changes to the original model to compare the sets of efficient solutions and the computational time required to obtain them. The main contribution of this paper is to extend the existing literature by incorporating approaches for the supply of product to customers through multiple sources, the direct flow between plants and customers, without this necessarily implying removing the distribution centers, and the product flow between distribution centers. From these approaches, we generate mathematical programming models and propose to solve through the epsilon-constraint approach for generating Pareto fronts and thus compare each of these approaches with the original model. The models are implemented in GAMS and solved with CPLEX.

#### 1. Introduction

Supply Chain Management (SCM) is the process of planning, implementing, and controlling the operation of the supply chain efficiently. SCM spans all movements and storage of raw materials, work-in-process inventory, and finished goods from the point of origin to the point of consumption [1]. Part of the planning processes in SCM aim at finding the best possible supply chain configuration so that all operations can be performed in an efficient way. The Capacitated Fixed Cost Facility Location Problem with Transportation Choices (CFCLP-TC) proposed by [2, 3] is a combinatorial optimization problem for supply chain design. It is an extension of the CFLP as a biobjective mixed-integer program. It is based on a two-echelon system for the distribution of one product in a single time period with two objectives: to minimize cost and to minimize the time of transportation from plants to customers. This approach considers several alternatives to transport the product from one facility to the other in each echelon of the network. The criterion of cost is an aggregate function of variable cost and fixed cost. At difference from similar works in the literature, the aim here is to provide the decision maker with a set of nondominated alternatives to allow her deciding. Some qualitative information only known by the decision maker may motivate the selection of one of those alternatives.

This paper presents three innovative approaches for the design of a supply chain for biobjective problems. In the first approach, it is proposed that customers can be supplied by more than one distribution center; in the second approach it is proposed that customers can be supplied directly by plants, but the distribution centers can be used. In the third approach proposed here a product flow may exist between the distribution centers. For the study of this problem and the proposed variations, instances of different sizes were used. The results are compared on the basis of three metrics. The first metric is called , proposed by [4], and the second and third metrics are called and , proposed by [3] to compare two biobjective Pareto fronts.

This paper is divided into five sections as follows: the first section presents a general introduction to the work. The second section presents the literature review to define the opportunity area that the work will fill. The third section presents an overview of the problem and its variations. The fourth section describes the computational experiment, the metrics used to evaluate the Pareto fronts, and the results of the computational implementation. Finally, in the fifth section the conclusions and proposals to carry out further work are exposed.

#### 2. Literature Review

Historically, researchers have focused on the design of distribution systems [5], without considering the supply chain as a whole. Typically, discrete location models were proposed to include additional features. Reference [6] reviews some important mixed-integer formulations for productions-distribution systems. However, those models had limited scope and could not cope with realistic supply chain structures. Reference [5] proposes the inclusion of relevant features for the Supply Chain Management (SCM) in the facility location models that gradually began to be considered. These include(a)subsets of the products (customer specifications),(b)upper and lower limits on shipments of a product to a particular plant,(c)specifications of the product weights for performance measures in the distribution centers,(d)piecewise linear approximation for nonlinear costs,(e)the ability to locate plants and distribution centers,(f)inclusion of capacity constraints of the products in plants,(g)conversion of raw materials in activities of one or two levels,(h)additional distribution and production levels.

Reference [7] suggested the inclusion of additional elements in facility location models such as the inclusion of new objectives (maximum return of investment) and decisions regarding the selection of equipment for new installations. In discrete location problems, selecting sites to establish new facilities in is restricted to a finite set of places available for the location. The simplest example of this approach is where* p* facilities must be selected to minimize the total distance (weighted) or the costs to supply customer demand.

This is a classic problem called the* p*-median problem, which has been extensively studied by [8–11]. This problem assumes that all candidate sites are equivalent in terms of installation cost for a new facility. When this is not the case, the objective function is extended with a term for the fixed cost of location, and as a result, the number of facilities to be open is an endogenous decision. This new approach is known in the literature as the uncapacitated facility location problem (UFLP). There are many references to this problem like in [12–14]. In both cases, the* p*-median and UFLP, each client is assigned to an open facility that minimizes the allocation cost. One of the most important extensions to the UFLP is the capacitated facility location problem (CFLP) in which exogenous values are considered to maximize the demand that can be supplied for each potential site. In this case, the closest assignment property is no longer valid as proposed by [15–17].

The models above have several common characteristics as follows:(a)single period of planning horizon,(b)deterministic parameters (demand, costs),(c)single product,(d)a single type of facility,(e)decisions of location-allocation.

Clearly, these models are insufficient to handle a realistic facility location scenario. Therefore, many extensions of this basic problem have been proposed and widely studied [18]. A crucial issue in many practical problems of localization is to consider the existence of different types of facilities [19], each playing a specific role (production or storage), and a natural flow of the material (i.e., a hierarchy) between them [20]. Each set of facilities of the same type is usually denoted by a level or echelon in the hierarchy of facilities.

Part of the planning processes in SCM aim at finding the best possible supply chain configuration so that all operations can be performed in an efficient way. The coming back to the capacitated facility location problem (CFLP) is a well-known combinatorial optimization problem. It consists in deciding which facilities to open from a given potential set, and how to assign customers to those facilities. The objective is minimizing total fixed and shipping cost. Applications of the CFLP include location and distribution planning, lot sizing in production planning, and telecommunication network design as mentioned by [21]. Numerous heuristics and exact algorithms for the CFLP have been proposed in the literature. Heuristic solution methods as well as approximation algorithms were proposed by [22–24]. Tabu Search methods for the related* p*-median problem and the CFLP with single source were developed [25, 26]. Exact solution methods based on the Benders decomposition algorithm are considered [27]. Polyhedral results for the CFLP have been obtained by [28]. Reference [29] uses these results in a branch and cut algorithm for the CFLP.

Moreover, several variants of the CFLP have been investigated. Reference [30] formulated a stochastic integer linear programming model for the CFLP with stochastic demands. A branch and cut approach was applied to find the optimal solution of the problem. Reference [31] formulated the mathematical model of the two-echelon single-source CFLP and considered six Lagrangian relaxation based approaches for the solution. In the recent years, many metaheuristic approaches have been applied to combinatorial optimization problems successfully, such as Simulated Annealing (SA), Genetic Algorithms (GAs), Tabu Search (TS), and Ant Colony Optimization (ACO). Some recent work in this field includes those presented by [32] in which they use an MAX-MIN ant system approach for the design of a supply chain. Reference [33] presented a memetic algorithm for a multistage supply chain problem. Reference [34] proposes a simulated annealing algorithm for an allocation problem. Reference [35] presents a hybrid approach using an artificial bee algorithm (BA) with mixed integer programming (MIP) applied to a large-scale CFLP; BA is applied for the purpose of solving the location problem, and the MIP is applied for the purpose of finding the optimal mathematical problem.

The biobjective location problems are extensions of classic locations problems. These problems are biobjective median, knapsack, quadratic, covering, unconstrained, location-allocation, hub, hierarchical, competitive, network, and undesirable and semidesirable location problems. Considering capacities in location problems, there are capacitated and uncapacitated problems in the literature. For instance, [36] has considered an uncapacitated facility location problem with two maxisum objectives (net profit and profitability of investment) and modeled it as parametric integer program with fractional and linear objectives. Reference [37] has modeled a supply network as a biobjective uncapacitated facility location problem with minisum and maxisum objectives (cost and coverage). In contrast, [38] developed an extension of the capacitated model to deal with locating maternity facilities with minisum objectives (distance traveled and load imbalance). Reference [39] has used a different bicriteria approach to the single hub location/allocation problem. This approach has two objectives; the first has a minisum form (cost), while the second objective (processing time) has two alternative forms.

The Capacitated Fixed Cost Facility Location Problem with Transportation Choices (CFCLP-TC) proposed by [2] is an extension of the CFLP with a biobjective mixed-integer program approach (cost and time). This approach considers several alternatives to transport the product from one facility to the other in each echelon of the network. At difference from similar works in the literature, the aim here is to provide the decision maker with a set of nondominated alternatives to allow him deciding. Some qualitative information only known by the decision maker may motivate the selection of one of those alternatives.

#### 3. Problem Description

The Capacitated Fixed Cost Facility Location Problem with Transportation Choices (CFCLP-TC) proposed by [2] is based on a two-echelon system for the distribution of one product in a single time period. In the first step, the product is sent from manufacturing plants to distribution centers .

The second step corresponds to the flow of product from distribution centers to customers . In this problem, the number and location of plants and customers are known a priori. This includes a further decision on the selection of channels of transportation between facilities using a bi-objective approach that simultaneously minimizes the time of transportation of product from plants to customers and the combined costs of locating facilities and transportation. This solution approach builds a set of alternative nondominated solutions for the decision maker. This problem has a set of possible locations for the opening of distribution centers and their number is not defined. Each candidate site has a fixed cost for opening a facility, and each site has limited capacity. Manufacturing plants have limited capacity and their production is sent from each plant to the distribution centers.

An important feature of the problem is to consider various alternatives for transportation of product from one facility to another in each step of the network. Each option represents a type of service with associated costs and time parameters. The existence of third-party logistics companies (3PL) causes that different transportation services are available in the market. The alternatives are generated by the supply of different companies, the availability of different types of services (urgent or regular), and the use of different modes of transportation (truck, train, plane, ship, or intermodal). Commonly, these differences involve an inverse correlation between time and cost; a faster service is more expensive. The outline of the distribution network is displayed in Figure 1.

##### 3.1. Approach That Allows Supplying from Multiple Distribution Centers

In the original model it is established as a restriction that each customer is served by a single source (distribution center ). At this point, we ask what would happen if the delivery to customers from multiple sources is allowed? The main idea of ****this variation is to allow customers to be supplied in some cases by more than one source (distribution center ).

The implementation of this variation to the original model will be evaluated over the objective functions (1) and (2) to determine which has a better Pareto front. Also we study the behavior of the time required to obtain results compared with that required in the original model.

*Sets*Set of plants Set of potential distribution centers Set of customers Set of arcs between nodes and Set of arcs between nodes and .

*Parameters*Cost of transporting one unit of product from plant to distribution center using the arc; Cost of transporting one unit of product from distribution center* j* to customer using the arc ; Time for transporting any quantity of product from plant to the distribution center using arc ; Time for transporting any quantity of product from distribution center to customer using arc ; Capacity of plant ; Capacity of distribution center ; Demand of customer ; Fixed cost for opening distribution center ; .

*Decision Variables*Quantity transported from plan to distribution center using arc ; Quantity transported from distribution center to customer using arc ; Binary variable equal to 1 if distribution center is open and equal to 0 otherwise; Binary variable equal to 1 if arc is used to transport product from plant to distribution center and equal to 0 otherwise; , , Binary variable equal to 1 if arc is used to transport product from distribution center to customer and equal to 0 otherwise; , *, *.

*Auxiliary Variables*A variable that computes the longest time that takes sending product from any plant to any customerLongest time in the first echelon of the supply chain for active distribution center ; that is, Longest time in the second echelon of the supply chain for active distribution center ; that is, .

*Model 1. *Consider

In this model, the objective function (1) minimizes the transportation costs and the cost of opening the distribution centers. The objective function (2) minimizes the time of transportation. Restriction (3) calculates the longest total time between and . The restrictions (4) and (5) allow calculating the transportation time between and and then from to . The restriction (6) allows the satisfaction of the demands of each client. The restriction (7) impedes to exceed the capacity of each plant . The restriction (8) impedes to exceed the capacity of the distribution centers (). The restriction (9) allows the balance of flow between and . Restrictions (10) and (11) provide that the transportation of the material can only be done through a single arc. Restrictions (12) and (13) provide that an arc will be inactive if there is no flow through it. Equations (14) and (15) state the shipment of product only through active arcs. Restriction (16) provides that the distribution centers having no product flow through them must be closed. Restrictions (17) and (18) provide a definition of the variables in the model.

##### 3.2. Approach Allowing the Direct Flow between Plants and Customers

This variation suggests that some plants supply customers in a direct way, that is, without necessarily passing through the distribution centers . The main reason is because in some cases customers may be closer to the plants than the distribution centers . This could make the flow more efficient in terms of cost and time without considering distribution centers . As in the previous variation, we want to evaluate the efficiency of this approach through the construction of a new Pareto front and the performance evaluation of the computational time required. Figure 2 shows the outline of this proposal.

The mathematical model for the variation that allows the direct flow between plants and costumers in some cases is as follows.

*Sets*Set of arcs* l* between nodes and *: **, *

*Parameters*Cost of sending one unit of product from plant to customer using the arc *: **, **, *Time to transport any quantity of products of plant to customer using the arc *; **, **, **. *

*Decision Variables*Quantity transported from plant to customer using the arc *; *Binary variable equal to 1 if the arc is used to transport product from plant to customer and equals 0 otherwise and *; **, **, **. *

*Auxiliary Variables*The maximum time it takes to ship the to.

*Model 2. *Consider

In this new formulation (19) replaces (1) in the original formulation as the objective function that searches for the best cost. Equation (20) is added to the original model and this calculates the time between the plants and customers . Equation (21) assigns the longest time between plants and customers to variable . Equation (22) replaces (6) in the original model and is aimed at satisfying the customer demand . Equation (23) replaces (7) in the original model and this requires that the transported amount from to does not exceed the capacity of the plant . Equation (24) is added to the model, and this ensures that customers can be supplied only by a single source. Equation (25) is added to the original model and states that the transportation of the material can only be done through a single arc. Equation (26) is added to the original model to establish that an arc is inactive whether its flow is zero. Equation (27) is added to the original model and states that the product flow will be made only through active arcs. Equations (28) and (29) replace (18) and (19) in the original model and set the domain of the variables.

##### 3.3. Approach That Allows Flow between Distribution Centers

This variation suggests the possibility of exchange of goods between distribution centers . For some alternatives this variation would compete directly with the flow of plants , distribution centers , plants with the alternative flow of plants , distribution centers , distribution centers , and plants . The main idea is that for some cases it is cheaper and faster to send a product from a distribution center to other distribution centers and then send it to customers . Figure 3 shows the outline of this proposal.

The mathematical model allowing flow between distribution centers, -, is as follows.

*Sets*Set of arcs between nodes :Set of arcs between nodes : .

*Parameters*Cost of sending one unit of product between distribution centers using :*, *Cost of sending one unit of product between distribution centers using :*, *Time to transport any quantity of products between distribution centers using the *, * Time to transport any quantity of products between distribution centers using the *; **, **. *

*Decision Variables*Quantity transported from to using *; **, *Quantity transported from to using *; **, *Binary variable equal to 1 if the arc is used to transport product between distribution centers and equals 0 otherwise; *, **, *Binary variable equal to 1 if the arc is used to transport product between distribution centers and equals 0 otherwise*; **, **, *

*Auxiliary Variables*Maximum time it takes to ship the product using the alternative , , , , ; Alternatively time , A very large positive valueBinary variable equal to 1 if the arc is used to transport product from to and equals 0 otherwiseBinary variable equal to 1 if the arc is used to transport product from to and equals 0 otherwise *. *

*Model 3. *Consider

In this new formulation (30) replaces (1) in the original model as an objective function that looks for the best possible cost. Equation (31) is added to the original model and calculates the longest time considering the flow -. Equations (32) and (33) are added to the original set as a condition if-then to determine the longest time in the flow of product from to , considering the flow between and . Equation (34) is added to the original model and this limits multiple routing between distribution centers . Equations (35) and (36) are added to the original model and complete the if-then condition to calculate the longest time in the flow of product from to considering the flow of -. Equation (37) replaces (8) in the original model and states that the amount transported from to does not exceed the capacity of distribution center by considering the flow between distribution centers -. Equation (38) replaces (9) in the original model and provides the flow balance in the distribution centers. Equations (39) and (40) are added to the original model and these provide that the transportation of the material can only be done through a single arc. Equation (41) is added to the original model and states that the exchange of product between and does not create a cycle. Equations (42) and (43) provide that an arc will be inactive if there is no flow through it. Equations (44) and (45) are added to the original model and establish that the distribution centers , having no product flow through them, will be closed. Equations (46) and (47) replace (18) and (19) in the original model and set the domain of the variables in the model.

#### 4. Computational Experiment

For the process of computational experiments, five sets of instances of each size were used as shown in Table 1. The encoding of the instance sizes is as follows. The first index indicates the number of plants , the second index indicates the number of distribution centers , the third index indicates the number of customers , and finally the fourth index indicates the number of arcs between nodes in each echelon. In each size, 5 instances were tested.

##### 4.1. Instance Design

The instances used in the experimental process were initially proposed by [2, 3]. These instances were used to test the original model and the variation that allows multiple sources of supply for customers.

For the variation that permits direct supply from the plants to customers , these instances were extended to accommodate this approach. The times for the arcs of the alternative - were generated at random on the basis of a normal distribution, with values ranging from 5 to 50. The aim of these values is that they are competitive with those of the original instance so that the model can choose to select sometimes the direct flow of - and in other cases to select the flow --. We consider that time and cost for the alternative - are negatively correlated. The unit cost of transportation in the flow - is a floating point variable calculated as follows:

This factor of “7” was determined after some experiments with some instances in such a way that the different solutions along the Pareto front included direct paths - and indirect paths --. This was done to impede that the model would prefer a certain type of path.

To test the model that allows the exchange of product between distribution centers, we intended to use the original data extended by adding new times and costs for the arcs that allow the exchange between distribution centers . With these data the model did not select arcs -, because in fact the solutions had few open distribution centers . To promote the opening of more distribution centers the fixed opening cost should be reduced in the new instances. With this change, we expected that some arcs - could be used in the solutions. Therefore, the extension to the original instance is as follows: the fixed cost of distribution center was generated as a random variable with integer values between 10,000 and 15,000.

These values will force the model to have more open distribution centers so that the product flow between them is more likely. From this point, the instances were expanded, creating additional time for flows - at random with a range of 1 to 5. The costs of the new alternatives are generated using the following relationship: This factor of “1/10” was determined after some experiments with some instances in such a way that the different solutions along the Pareto front included paths and paths. This was done to impede that the model would prefer a certain type of path.

To perform the computational experiment we used a computer equipped with the following features: Workstation with Intel (R) Xenon (R) CPU X5550 2.67 GHz with 12 GB of RAM and 64-bit operating system (Windows 7). The implementation of the models was performed in GAMS 12 and solved with CPLEX 23.6.2.

##### 4.2. Metrics

To make the comparison of the new Pareto fronts with the original ones, the metric proposed by [4] was used. Additionally, we registered the average number of Pareto-optimal solutions in each front. To calculate the , consider that and are the sets of Pareto-optimal solutions obtained from each model and is the union of the sets of Pareto-optimal solutions such that it includes only nondominated solutions* Y'*s*.* The ratio of Pareto-optimal solutions in that are not dominated by any other solutions in is calculated as follows:
where means that the solution is dominated by solution . The higher the ratio is, the better the solution set is. Similarly, we used the metrics proposed by [3] called and . These were developed to give practical meaning to the comparison of sets point by point. The discretization of objective and the number of objectives allow proceeding as follows for a pair of sets and .

Let and be the objective functions of the problem and and the sets of nondominated solutions to be compared. By discretization of function we can construct a set , such that its elements are those values of that exist in and ,

Then computes an average rate deviation of the objective function for each value of , which is in the set ,

The metric indicates the quality of a set compared to another. The following relationship can be established:

It is important to establish that an important parameter is the computational time required to solve each model for an instance.

##### 4.3. Results

Table 2 shows the comparison results for the solution of the original model and the first variation (without single source constraint) of the CFCLP-TC problem. In this evaluation five instances are presented with metrics , , and and the processing time in seconds of each one.

Table 2 shows that in all cases D_{avg} is less than 1; this indicates a superior quality of the Pareto fronts of the variation compared to the original model. D_{min} in all cases indicates a small difference between both fronts, since the values are close to 1. The values of in the original model are on average 64% for instance 5-5-5-2, 69% for instances 5-5-5-5, and 0% for instances 5-10-10-2, 5-10-15-2, and 5-10-20-2. This indicates that in all cases the variation that allows the supply of product to customers without single source constraint has always better Pareto fronts compared with those obtained in the original model.

The processing time for the variation increased on average by 9% for instances 5-5-5-2, increased on average 627% for instances 5-5-5-5, and increased by 76% on average for instances 5-10-10-2 compared to the original model. However for instances 5-10-15-2 time decreased by 16% and by 500% for instances 5-10-20-2.

In Table 3 the results of the comparison between the original model and the variation that allows the direct flow between the plants and customers are shown. It should be noted that, in the majority of the cases, is less than 1; this indicates that the Pareto fronts of the variation are better compared to the original model. in all cases indicates the smallest difference comparing both fronts and provides a measure of the difference of the fronts compared.

In relation to it is observed that the variation of the model presents values of 1 in all cases, compared with the original model. And the values of in the original model are below the proposed variation on average 53% for instance 5-5-5-2, 73% for instance 5-5-5-5, 64% for the instances 5-10-10-2, 26% for instance 5-10-15-2, and 18% for instance 5-10-20-2. This indicates that in all cases the direct flow variation between and has always better Pareto fronts as compared with those obtained in the original model.

Concerning the processing time, comparing the original model with the variation that allows the flow from to , we have the following: for instance 5-5-5-2 time is increased on average by 1.6% and for instance 5-5-5-5 increased on average by 6.2% over the original model, for instance 5-10-10-2 a decrease is observed on average by 3.4%, and for instances 5-10-15-2 and 5-10-20-2 processing time decreased by 22% and increased by 16%, respectively.

Table 4 shows the results of the variation that allows the flow between distribution centers -. is greater than 1 in all cases; this indicates a lower quality of the Pareto fronts of this variation compared with the original model.

The values in the original model are on average of 1 with respect to the variation; this indicates that the original model presents always better Pareto fronts. With respect to processing time, for instance 5-5-5-2 time is increased by average of 460%. For instance 5-5-5-5 time is increased on average by 560% over the original model.

For the analysis of this variation only the first two groups of instances were solved, because the processing time required for evaluating the variation was very long. Solving some instances of the group 5-5-5-5 required up to 34 hours to get results, however for instances 5-10-10-2 after 143 hours of processing the solver did not yield a result, so we decided not testing for the rest of the instances.

#### 5. Conclusions

Reference [1] defines the management of the supply chain as the process of planning, implementation, and operational control of the supply chain in an efficient manner. This aspect is defined in the context of tactical decisions that allow for more efficiency in the full cycle of manufacture. The work developed in this research explores an area that has not been sufficiently analyzed and incorporated into mathematical models of supply chain design with selection of distribution channels according to [40, 41].

The CFCLP-TC proposed by [2, 3] incorporates in a novel way the selection of transportation alternatives in the context of a (two-echelon) problem considering plants, distribution centers, and customers. However, the changes proposed in this paper can bring theoretical models to real applications, as they consider situations that could occur in real context.

The first variation that allows the flow to customers from multiple sources (DC) states that the proposed approach generally results in better cost compared with the original proposal where customers are limited to obtaining the product from a single source (DC). This enables better Pareto fronts in the proposed alternative. Although the results obtained suggest that it is cheaper to allow customers to be supplied by more than a distribution center, this increases the complexity of cross-docking, and in accordance with [40] the increase of the level of complexity of the supply chain will negatively affect the performance of a manufacturing plant. However, a greater variety of products may create economies of a scale enough to reduce this effect.

The second variation allows in some cases the product flows from plants directly to customers without necessarily passing through the distribution centers. This variation generally obtained better costs compared to the original proposal. The proposal can be justified when customers are geographically closer to the plants than the distribution centers.

When planning the configuration of the supply chain, it is not possible to accurately predict in the future where they will reside geographically, and it is clear that it is more costly to relocate distribution centers based on the geographical configuration of new customers. Therefore, an alternative that would optimize the cost and time of the supply chain is based on the ability to send directly from plants to customers. The results obtained allowed us to determine the best Pareto fronts with this approach; however, a problem is the time required to solve the instances, since the increase of time is significant as we try to solve the larger instances that resembled real problems. The proposal explores a configuration of the new supply chain that has not been considered in the literature.

The third variation allowing in some cases the exchange of product between distribution centers generally gets worse costs compared with the original model. This proposal is considered in the overall context of managing the supply chain, where a distribution center can supply to another distribution center if the times to get the product to the customer and the costs associated with this alternative are lower. An example is the exchange of product between car dealers to reduce lead times, since placing an order to assembly plant takes time and cost, and these impact the service level.

The results determined that this approach provides Pareto fronts worse than the original model and therefore they have higher costs. For this approach, processing times are greatly increased as the size of the instances grows; this situation is not rare in practice. But one aspect that the model does not consider is the level of customer service, since for the example above it is difficult to predict with accuracy the products requested by the customers, and when this happens, you must deliver the product the customer demands in the shortest possible time and cost.

The proposed approaches are novel and allow mathematical modeling to configure supply chains to situations that occur most commonly in practice and allow closing the gap between theory and practice, contributing to the state of the art in this context.

It is clear that the new approaches produced an increase in the required processing time. The findings and conclusions presented are based only on instances that were tested, which in a context of real implementation are of small size. It is possible that the results of the Pareto front for larger instances may be different; therefore, it is important to determine these fronts, and however the impossibility of doing so with exact methods requires us to try to get them with heuristics and metaheuristics. In the revised literature the Lagrangian relaxation method is widely used to solve similar problems [42], and this may be applied to the problem addressed in this work. Another approach is the use of genetic algorithms proposed by [43] for a similar problem found in the literature. Evolutionary algorithms like Nondominated Sorting Genetic Algorithm (NSGA-II) and Strength Pareto Evolutionary Algorithm (SPEA-II) are widely used in multiobjective problems as in [17, 44, 45].

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors appreciate the support of the administration of Polytechnic University of Tulancingo, for the equipment and facilities provided for the development of this work, as well as the comments from the Faculty of Applied Mathematics and Approximation in the same university.