Computational Intelligence in Modeling Complex Systems and Solving Complex Problems
View this Special IssueResearch Article  Open Access
Congdong Li, Hao Guo, Ying Zhang, Shuai Deng, Yu Wang, "An Improved Differential Evolution Algorithm for a Multicommodity LocationInventory Problem with False Failure Returns", Complexity, vol. 2018, Article ID 1967398, 13 pages, 2018. https://doi.org/10.1155/2018/1967398
An Improved Differential Evolution Algorithm for a Multicommodity LocationInventory Problem with False Failure Returns
Abstract
Customer returns are a common phenomenon in many industries, and they have a significant impact on business organizations and their supply chains. False failure returns are returned products that have no functional or cosmetic defects, and they represent a large body of customer returns in practice. In this paper, we develop a mixedinteger nonlinear programming model to study a multicommodity locationinventory problem in a forwardreverse logistics network. This model minimizes the total cost in this network by considering false failure returns, and it also considers many realworld business scenarios in forward and reverse logistics flows. Moreover, we design a new heuristic approach to solve the model efficiently. Finally, numerical experiments are conducted to validate our solution approach and provide meaningful managerial insights.
1. Introduction
Reverse logistics is an emerging topic that has attracted many research attentions because of the great importance of customer returns in today’s business. Nowadays, the percentage of customer returns in total sales varies from 5% to 9% for durable goods, and it is up to 35% for fashion apparel products [1]. Since customer returns are more frequent in ecommerce, return rate in online shopping has exceeded 18% of total sales, and especially, it reaches 74% for fashion products [2]. In the retail industry, the value of customer products has exceeded $100 billion per year in the United States [3]. Although many products are returned because of quality issues, a lot of items that have no defects can still be returned for various reasons such as installation difficulties, bad user experience, or price competition. Those items represent a large body of customer returns, and they are defined as “false failure returns” [4].
Forwardreverse logistics network consist of forward flows that are directed from suppliers to customers and reverse flows of returned products from customers to suppliers. Since network design can contribute to reducing 60% the cost of a firm [5], Supply chain network design is an important strategic decision in supply chain management, and it is a big challenge to business managers. For example, many big companies such as Kodak, Xerox, and General Motors have developed their reverse logistics systems and those efforts are very successful [6]. Currently, many business decisions are made jointly to improve the efficiency of business organizations and their supply chains. Locationinventory problem (LIP) [7] is such a problem that integrates strategical facility location decisions and tactical inventory management decisions, and it can be studied by incorporating many realworld business scenarios to improve the performance of a logistics system.
In this paper, we study a multicommodity locationinventory problem (MCLIP) in a forwardreverse logistics network in which the reverse flow consists of false failure returns. This problem is formulated as a mixedinteger nonlinear programming model whose objective function is to minimize the total cost in such a network, and it incorporates many business decisions and processes such as facility locations, return processing, and order fulfillment. More specifically, the solution to this problem will determine the number and locations of hybrid distributioncollection centers (HDCCs), how to use those HDCCs to fulfill orders and collect returns, and how to replenish inventories in those HDCCs. Due to the complexity of LIPs, we also design an improved differential evolution algorithm (IDE) to solve the MCLIP efficiently. Computational results are presented to validate the solution approach, and research implications and managerial insights are also discussed.
The rest of this paper is organized as follows: In Section 2, we review related research works in the literature. In Section 3, we describe the research problem and present the mathematical model. In Section 4, we propose a new heuristic approach to solve the problem efficiently. In Section 5, we show computational results and discuss managerial insights. Section 6 concludes the paper and provides directions for future research.
2. Literature Review
The design of logistics networks involves several core components such as facility location and inventory management. LIPs have been extensively studied because of their great importance in practice, and Farahani et al. [8] provide a comprehensive review of such works in the literature. In this section, we review the research works on LIPs from the perspective of forward and reverse logistics, respectively.
In the literature, most LIPs are studied under forward logistics networks. For example, Daskin et al. [7] study a LIP with the riskpooling effect and design a Lagrangian relaxation algorithm as the solution approach. Shen et al. [9] use the same integer programming model as that in Daskin et al.[7] and convert it into a setcovering integer programming model and then solve it by column generation. Moreover, LIPs are studied by incorporating many business scenarios and processes. For example, Ozsen et al. [10] consider capacitated warehouses and the impact of multisourcing in addition to Shen et al. [9]. AmiriAref et al. [11] study multisourcing locationinventory problem with stochastic demand. Silva et al. [12] and Qu et al. [13] study the stochastic locationinventory problems under different replenishment policies. Mousavi et al. [14] study a seasonal multipleproduct location allocation and inventory control problem in a planning horizon consisting of multiple periods. Since business environments are always dynamic, demand uncertainty [15–17] and different risks [18] are also considered.
It is emerging to study LIPs in the reverse logistics. Since reverse logistics is an important part of closedloop supply chains (CLSCs), those works are usually conducted with CLSCs. For example, Diabat et al. [19] study a closedloop locationinventory problem by considering spare parts, AslNajafi et al. [20] study a dynamic closedloop locationinventory problem under disruption risk, and Kaya et al. [21] and Guo et al. [22] extend such studies by incorporating pricing and vehicle routing decisions in CLSCs, respectively. Since returned products are usually sold to secondary markets in practice, LIPs are also studied by considering CLSCs and secondary markets [23, 24].
Although LIPs have been extensively studied from the perspective of reverse logistics, those works are still very limited because the types of customer returns are rarely considered. This paper studies a LIP in a forwardreverse logistics network in which the reverse logistics flow only consists of false failure returns. To the best of our knowledge, this is the first paper that considers false failure returns in the LIP literature, and it makes a significant contribution to the literature because of the great importance of false failure returns in practice.
3. Mathematical Model
3.1. Problem Description
We consider a threeechelon supply chain network in which a manufacturing plant ships multiple commodities to customer zones via a set of facilities. Moreover, returned products will be collected from customers and shipped to the same facilities. Therefore, those facilities are hybrid distributioncollection centers (HDCCs) [25, 26], and they have many advantages such as cost saving and pollution reduction by improving the efficiency of workflows. In this network, we assume that only false failure returns will be collected. Since those returns do not have any quality issues, they will not be sent to the plant for repair, and instead, they will be resold in the market after some simple processes (e.g., inspection, repackaging, etc.) in HDCCs. Figure 1 shows the forwardreverse logistics network under study. In this network, forward flows will be unidirectional from the plant to HDCCs and then from HDCCs to customer zones, but reverse flows will only be directed from customer zones to HDCCs because false failure returns will not be sent to the plant for repair.
This paper studies a multicommodity locationinventory problem (MCLIP) in the network shown in Figure 1, and we make the following assumptions to simplify the problem: (i) customer zones are predetermined and fixed; (ii) the returned products will be the same as new items when they are resold to the market; (iii) the distances between HDCCs and customer zones are Euclidean distances. In practice, this problem can be solved to answer the following questions:(1)What are the best locations for HDCCs in this network?(2)How to assign those HDCCs to customer zones in the forward and reverse networks, respectively?(3)How to optimize inventory replenishment decisions at HDCCs?
3.2. Notations
Sets : set of commodities, indexed by ; R: set of candidate HDCC locations, indexed by r; I: set of customer zones, indexed by i.
Parameters : fixed (yearly) cost of building a HDCC at location r, for each r∈R; : fixed administrative and handling cost of placing an order for product at HDCC r, for each r∈R and ∈; : basic transportation cost per shipment of product from the plant to HDCC r, for each r∈R and ∈; : additional transportation cost per unit of product from the plant to HDCC r, for each r∈R and ∈; : shipping cost per unit of product from HDCC r to customer zone i, for each r∈R, i∈I and ∈; : repackaging/refurbishing cost per unit of returned product at HDCC r, for each r∈R and ∈; : inspection cost per unit of returned product at HDCC r, for each r∈R and ∈; : (yearly) holding cost per unit of product at HDCC r, for each r∈R and ∈; : return rate of product , for each ∈; : mean (daily) demand of product at customer zone i, for each i∈I and ∈; : variance of (daily) demand of product at customer zone i, for each i∈I and ∈; : desired percentage of market demand satisfied; : standard normal deviation such that ; L: order lead time in days at HDCCs; λ: working days per year; : mean (yearly) demand of product at HDCC r, for each r∈R and ∈.
Decision Variables : order quantity of product per order at HDCC r; X_{r} = 1 if opening one HDCC at location r, and 0 otherwise, for each r∈R; = 1 if HDCC r fulfills demand at customer zone i, and 0 otherwise, for each r∈R, i∈I and ∈; = 1 if HDCC r collects returns from customer zone i, and 0 otherwise, for each r∈R, i∈I and ∈.
3.3. The Model
In this paper, the research problem is formulated by a mixedinteger nonlinear programming model. The objective of this model is to minimize the total cost in the forwardreverse logistics network, and it comprises (i) location costs, (ii) shipping costs from HDCCs to customer zones, (iii) inventory costs including working inventory and safety stock costs, and (iv) return costs.
3.3.1. Location Cost
Obviously, the total location cost () can be written as
3.3.2. Shipping Cost from HDCCs to Customer Zones
The total shipping cost (C_{T}) from HDCCs to customer zones in the forward logistics flow can be expressed as follows:
3.3.3. Inventory Cost
Inventory strategy plays an important role in supply chain management. In this study, we assume that HDCCs order new products from the plant by using an (R, Q) policy with type I service, and the total inventory cost includes working inventory and safety stock costs. The (R, Q) model can be approximated by an economic order quantity (EOQ) model [9], and Axsater [27] shows that the maximum relative error introduced by using the EOQ model instead of the (Q, r) policy is 0.118. Therefore, we use an EOQ model to optimize the inventory strategy in this study, and the order frequency and quantity at a HDCC will be calculated by using the total demand of the customer zones served by this HDCC.
In general, the working inventory cost includes fixed cost of placing orders, shipping cost from the plant to HDCCs, and holding cost of the working inventory. The individual cost terms are formulated as follows.
(1) Order Cost. The annual fixed cost of placing orders for product from HDCC r to the plant is given by .
(2) Shipping Cost from the Plant to HDCCs. Since false failure returns will be resold to the market after they are processed in HDCCs, shipping quantities from the plant to HDCCs will equal the difference between market demands and return volumes. Therefore, the total shipping cost from the plant to HDCC r per year is given by
(3) Holding Cost of the Working Inventory. In this study, the working inventory comes from two sources which are the new products ordered from the plant and the returned products collected from customer zones. Therefore, the forward and reverse logistics flows will be integrated at HDCCs, and the average inventory of product at HDCC r during period T can be expressed as follows:where is the mean daily demand of product at HDCC . Consequently, the holding cost of the working inventory is given by
Using Eppen’s riskpooling result [28], the amount of safety stock to ensure that stockouts occur with a probability of α or less is , and the safety stock cost at HDCC r is given by .
Therefore, the total annual inventory cost (C_{I}) of product at HDCC r is given by
3.3.4. Return Cost
Returns costs are incurred by customer returns, and they include inspection cost, repackaging cost, and shipping cost from customer zones to HDCCs.
(1) Inspection Cost. Returned products will be inspected at HDCCs, and the total inspection cost is given by .
(2) Repackaging Cost. The returned products also need to be repackaged before they are resold to the market, and the total repackaging cost is given by .
(3) Shipping Cost from Customer Zones to HDCCs. The total shipping cost from customer zones to HDCCs is given by .
Therefore, the total annual return cost (C_{R}) is given by
3.3.5. Optimization Model
Given the individual cost functions above, the research problem under study can be formulated by a mixedinteger nonlinear program as follows:subject to
In this model, objective function (8) means that the total cost will be minimized. Constraint (9) means that at least one HDCC will be built. Constraint (10) means that the demand of a product in a customer zone will be fulfilled by only one HDCC. Constraint (11) means that the returns of a product from a customer zone will be collected by only one HDCC. Constraint (12) means that all types of products will be provided to a customer zone by HDCCs in the forward flow. Constraint (13) means that all types of returned products will be collected from a customer zone by HDCCs in the reverse flow. Constraints (14) and (15) mean that a customer zone can be assigned to a HDCC in the forward and reverse logistics flows, respectively, only if this HDCC is built. Constraint (16) means that a HDCC will serve at least one customer zone in the forward logistics network after it is built. Constraint (17) means that a HDCC will collect the returns of a product in the reverse flow only if it fulfills the demand of this product in the forward flow. Constraints (18)–(20) indicate that X_{r}, , and are binary variables.
Note that decision variable only appears in the objective function. Since the objective function is convex in , its optimal value can be obtained by taking the partial derivative of the objective function with respect to it, which givesTherefore, the objective function can be rewritten as follows:
4. Solution Approach
Location problems are NPhard in general [29], and LIPs can be more complicated. Evolutionary algorithms are widely used to solve LIPs, and differential evolution (DE), which was first introduced by Storn and Price [30], is popular method to solve nonlinear and nonconvex functions [31]. Although DE has a strong global search ability to solve such problems, its performance is not always guaranteed due to its weaknesses such as the poor local search ability. To obtain a more stable and effective approach, we design an improved differential evolution algorithm (IDE) to solve MCLIP. More specifically, IDE enhances the performance of DE by two ways: First, IDE introduces oppositionbased learning (OBL) [32] to generate better initial populations. Second, it uses a novel mutation factor F based on the Gaussian distribution N(0,1) to enhance the diversity of populations.
4.1. Individuals and Operators
Generally, a DE algorithm has four main operations: initialization, mutation, crossover, and selection. After an initial population is created in the initialization step, the algorithm will execute mutation, crossover, and selection operations iteratively to improve the population until a stop criterion is satisfied. IDE follows the same procedure, and it uses notations shown in Table 1.

4.1.1. Initialization
When DE is applied to solve optimization problems, an encodingdecoding scheme will be needed to convert individuals in DE to the solutions to the optimization model, and vice versa. In IDE, an individual is represented by the matrix shown in (23).
This matrix has rows and 2N columns. For the research problem under study, and N denote the numbers of commodities and customers, respectively, and the first and last N columns are related to the forward and reverse flows, respectively. More specifically, if , then the demand of product in customer zone j will be fulfilled by HDCC r in the forward flow, and if , then the returns of product from customer zone j will be collected by HDCC r in the reverse flow.
Figure 2 shows an example matrix for a MCLIP that has two commodities, ten customer zones, and five candidate HDCC locations. For product 1, HDCC 1 fulfills demands in customer zones 1, 6, and 10 and collects returns from customer zones 1, 7, and 9; HDCC 2 fulfills demands in customer zones 4 and 7 and collects returns from customer zones 2 and 6; HDCC 3 fulfills demands in customer zones 2, 3, 5, and 8 and collects returns from customer zones 4, 5, and 8; HDCC 5 fulfills the demand in customer zone 9 and collects returns from customer zones 3 and 10. For product 2, location decisions can be also explained in the same manner.
The aim of encoding is transforming the initial individual to practical solution of LIP. The formulation of decoding is as follows.
To decode an individual to a solution to the optimization problem, we can use a mechanism shown in (24) and (25).where are distinct HDCCs that provide product w in the forward flow, rand is a random variable that is uniformly distributed on , x^{L} and x^{U} are the lower and upper bounds of , and round is the rounding function. Obviously, we have x^{L} = y^{L} = 1, and x^{U} is decided by the number of candidate HDCC locations. The expression of ensures that product can be returned to a HDCC in the reverse flow only if it provides this product in the forward flow, which is reinforced by Constraint (17) in the optimization model shown in Section 3.3.5.
IDE introduce a new mechanism based on OBL to generate better initial populations. More specifically, IDE will improve an initial population that is generated randomly by checking its opposite population by the following steps.
Step 1. Generate an initial population Z^{i,0} randomly:
Step 2. Obtain its opposite population :where is the entry in row and column j of the individual in the opposite population and (i=1, 2, …, Np; j=N+1, N+2, …, 2N; =1, 2, …, ) is calculated by (25).
Step 3. Select the Np best solutions from as initial population.
4.1.2. Mutation
The mutation strategy plays an important role DE to enhance the global optimization ability and convergence speed. Usually, DE will use DE/rand/1/bin as the mutation operation that will produce a new vector by adding the weighted difference of two randomly selected vectors to the third vector in generation g [30]. In this study, we design a new mutation strategy to improve the diversity of populations, and new mutant vectors and are generated by (28) and (29).where i≠r1≠r2≠r3, r1, r2, r3 are randomly selected from , = F· and it amplifies the difference between two vectors, F is an initial mutation factor, and is independently generated from the normal distribution . In the mutation step, the mutation factor will be updated first and then used to generate new mutant vectors, which will improve local and global search abilities and generate better mutant vectors. If and in a new mutant vector exceed the upper and lower bounds, then the vector will be generated again within its range.
4.1.3. Crossover
The crossover operation mixes mutate and target vectors to generate trial vectors and increase the diversity of populations. In IDE, the trial vector is generated by (30) and (31).where rand(i) is random variable that is uniformly distributed on , and is a predefined crossover rate.
4.1.4. Feasibility Correction
A new individual that is generated by mutation or crossover operations may not be feasible by violating Constraint (17). Therefore, a feasibility correction procedure is developed to exclude infeasible individuals. In this study, a new individual will be checked after it is generated. If Constraint (17) is violated, then it will be replaced by a new individual that is randomly generated by (25).
4.1.5. Selection
After a new population is generated, the selection operation will be executed to evaluate the objective values of all trial vectors. In IDE, a trial vector f(U) will be compared with the corresponding target vector f() on their objective values by using a greedy criterion[30]. More specifically, the selection operation can be expressed as follows:
4.1.6. Stop Criterion
In DE, the mutation, crossover, and selection operations will be executed iteratively until a stop criterion is satisfied. In this study, IDE will stop until any of the criteria below is satisfied.
(1) No better solution is found in consecutive K iterations.
(2) The maximum number of iterations, i.e., G, is reached.
4.2. Algorithm Flow
In summary, IDE consists of the following steps.
Step 1. Set g = 0 and k = 0.
Step 2. Create an initial population that consists of Np individuals by applying OBL, and calculate the objective values of all individuals in the population.
Step 3. Execute mutation operation.
Step 4. Execute crossover operation.
Step 5. Perform feasibility correction.
Step 6. Execute selection operation.
Step 7. Set g = g + 1. Calculate the objective values of all individuals in the new population. If the minimal objective value is less than that in the old population, set k = 1. Otherwise, set k = k + 1.
Step 8. If k = K or g = G, stop the algorithm. Otherwise, go to Step 4.
5. Computational Results
In this section, numerical analysis is presented to validate IDE and provide managerial insights. First, we show how to obtain the optimal setting in IDE for the best possible performance. Second, the solution to an example MCLIP is presented. Third, sensitivity analysis is conducted to study the influences of the parameters in MCLIP. Last, IDE is compared with Lingo 11 and DE to validate its performance. All the experiments are implemented by Java JDK 1.7 on a Windows PC (AMD A109600P RADEON R5, 10 COMPUTE CORES 4C+6G 2.40GHz; RAM: 4.00 GB DDR; OS: Windows 10).
To get started, we use an example MCLIP that consists of eight candidate HDCC locations, eighty customer zones, and three commodities. The candidate HDCC locations and customer zones are uniformly distributed on a grid of , and the parameters are shown in Table 2.
 
Remark: denotes the uniform distribution over . 
5.1. Parameter Analysis
The performance of DE can be significantly affected by its parameters or the combinations of those parameters. In this subsection, we test the performance of IDE under different F and CR values. In this experiment, we set Np = 4S () and the algorithm was executed 30 times under each setting, and the numerical results are shown in Table 3.
 
Remark: in Table 3, the optimal solutions are obtained by Lingo 11, and “OS Ratio” means the percentage of finding optimal solutions by IDE. 
From Table 3, we can see that F and CR have a significant impact on the searching ability and convergence speed of IDE, which are indicated by OS ratio and CPU time, respectively. Obviously, IDE has the best performance when F = 0.1 and CR = 0.1, and this setting will be used in the subsequent experiments.
5.2. Illustrative Example
This subsection shows the optimal to the example MCLIP. Table 4 shows the optimal solution which indicates HDCC locations and the corresponding customer zones in the forward and reverse flows, and Table 5 shows the optimal value and individual costs in the optimal solution.


5.3. Sensitivity Analysis
The optimal solution of a MCLIP can be affected significantly by important business parameters. In this subsection, sensitivity analysis is conducted on , , , , , , , and to study their influences. More specifically, a parameter will be tested by changing its value by (30%, 30%) and fixing the other parameters and the numerical results are shown in Table 6.
 
Remark: in Table 6, “OV” means the optimal value, and “Gap” means the change of the optimal value. Note that the gap will be zero if no parameter is changed. 
From Table 6, we can see that optimal values can change significantly due to the change of those parameters. For example, the optimal value will decrease by 6.63% if decreases by 30%, and the optimal value will increase by 6.65% if increases by 30%. The influences of the parameters on individual costs are summarized as follows: When increases, C_{R} will increase, C_{I} and will decrease, but C_{L} and C_{T} will not change. When increases, C_{T} and C_{R} will increase but C_{L} and C_{I} will not change. When , , , increase, only C_{I} will increase and the other costs will not change. When or increases, only C_{I} will increase. We can see that, in general, and have the most significant influences. Moreover, the optimal value is more sensitive when , , and are relatively small. In practice, business managers may focus more on those parameters to stabilize the supply chain performance.
5.4. Performance Analysis
In this subsection, IDE is compared with Lingo 11 and DE to validate its performance in terms of solution accuracy and CPU time efficiency, and the three approaches are tested on a set of small (i.e., 20×5×2/3/4, 40×6×2/3/4), medium (i.e., 50×5/6/7×3, 70×6/7/8×3), and largesized problems (i.e., 80×8×2/3/4, 100×8/9/10×3). Since population size Np is an important parameter in IDE and it may vary according to the problem sizes, its value needs to be set properly to get the best possible performance. According to Figures 3–5 (in which “red line” represents solution accuracy in terms of the percentage of finding the optimal solutions, and “blue line” represents the mean CPU time), IDE will have the best performance when Np = 3S, 4S, 5S for small, medium, and largesized problems, respectively. Therefore, we will use these settings in the subsequent experiments.
Table 7 shows the comparison results between IDE, DE, and Lingo 11. To avoid biased results, IDE and DE were executed 30 times for each test instance. Therefore, in this table, “O.V.” means the optimal value, “O.V. (Best)” and “O.V. (Mean)” represent the best and mean optimal values in 30 runs, respectively, “OS Ratio” means the percentage of finding the optimal solutions, and “S.D.” means the standard deviation of optimal values in 30 runs.

From Table 7, we can see that IDE is the most effective and efficient method to solve MCLIPs, and its advantages over Lingo 11 and DE are summarized as follows:(1)The optimal values obtained by IDE and Lingo 11 were almost identical, but IDE has better runtimes than Lingo 11. This indicates that IDE is much more efficient than Lingo 11 with the same solution accuracy.(2)IDE and DE can achieve similar solution accuracies on small and mediumsized problems. However, the optimal values from DE are significantly larger than those from IDE on largesized problems. Particularly, DE cannot find optimal solutions for the problems whose sizes are over (100×9×3). This indicates that IDE has a great global search ability and is more stable and consistent than DE.(3)IDE has better runtimes than DE in most test instances, and it has a greater advantage on largesized problems. This indicates that IDE is more efficient than DE.
According to the comparison results, we can see that IDE is the best approach to solve MCLIPs compared with DE and Lingo 11.
6. Conclusion and Future Research
In this paper, we study a multicommodity locationinventory problem by considering false failure returns in a forwardreverse logistics network. To the best of our knowledge, this is the first work that considers false failure returns to optimize joint locationinventory decisions, and it contributes to the literature because of the great importance of false failure returns in practice.
This research can be extended in several directions: First, other types of returns such as defective items can be considered to make this study more practical. Second, the (R, Q) model is adopted to replenish inventories in this work. However, there are other popular inventory control policies in practice, and it will be interesting to explore them and compare their performances in forwardreverse logistics networks. Third, it will be valuable to study locationinventoryrouting problems by considering false failure returns.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by the National Natural Science Foundation of China under Grants nos. 71672074 and 71772075.
References
 V. D. R. Guide Jr., G. C. Souza, L. N. van Wassenhove, and J. D. Blackburn, “Time value of commercial product returns,” Management Science, vol. 52, no. 8, pp. 1200–1214, 2006. View at: Publisher Site  Google Scholar
 Z. Hua, H. Hou, and Y. Bian, “Optimal Shipping Strategy and Return Service Charge under NoReason Return Policy in Online Retailing,” IEEE Transactions on Systems, Man, and Cybernetics: Systems, vol. 47, no. 12, pp. 3189–3206, 2017. View at: Publisher Site  Google Scholar
 J. Stock, T. Speh, and H. Shear, “Many happy (product) returns,” Harvard Business Review, vol. 80, no. 7, p. 16, 2002. View at: Google Scholar
 M. Ferguson, V. D. R. Guide Jr., and G. C. Souza, “Supply chain coordination for false failure returns,” Manufacturing and Service Operations Management, vol. 8, no. 4, pp. 376–393, 2006. View at: Publisher Site  Google Scholar
 T. P. Harrison, The Practice of Supply Chain Management: Where Theory And Application Converge, Springer Science & Business Media, 2005.
 M. S. Pishvaee, F. Jolai, and J. Razmi, “A stochastic optimization model for integrated forward/reverse logistics network design,” Journal of Manufacturing Systems, vol. 28, no. 4, pp. 107–114, 2009. 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
 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
 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
 L. Ozsen, M. S. Daskin, and C. R. Coullard, “Facility location modeling and inventory management with multisourcing,” Transportation Science, vol. 43, no. 4, pp. 455–472, 2009. View at: Publisher Site  Google Scholar
 M. AmiriAref, W. Klibi, and M. Z. Babai, “The multisourcing location inventory problem with stochastic demand,” European Journal of Operational Research, vol. 266, no. 1, pp. 72–87, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 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
 Q. Hui, W. Lin, and L. Rui, “A contrastive study of the stochastic locationinventory problem with joint replenishment and independent replenishment,” Expert Systems with Applications, vol. 42, no. 4, pp. 2061–2072, 2015. View at: Publisher Site  Google Scholar
 S. M. Mousavi, N. Alikar, S. T. A. Niaki, and A. Bahreininejad, “Optimizing a location allocationinventory problem in a twoechelon supply chain network: a modified fruit fly optimization algorithm,” Computers & Industrial Engineering, vol. 87, pp. 543–560, 2015. View at: Publisher Site  Google Scholar
 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
 A. Diabat, E. Dehghani, and A. Jabbarzadeh, “Incorporating location and inventory decisions into a supply chain design problem with uncertain demands and lead times,” Journal of Manufacturing Systems, vol. 43, pp. 139–149, 2017. 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
 M. Farahani, H. Shavandi, and D. Rahmani, “A locationinventory model considering a strategy to mitigate disruption risk in supply chain by substitutable products,” Computers & Industrial Engineering, vol. 108, pp. 213–224, 2017. View at: Publisher Site  Google Scholar
 A. Diabat, T. Abdallah, and A. Henschel, “A closedloop locationinventory problem with spare parts consideration,” Computers & Operations Research, vol. 54, pp. 245–256, 2015. View at: Publisher Site  Google Scholar
 J. AslNajafi, B. Zahiri, A. BozorgiAmiri, and A. TaheriMoghaddam, “A dynamic closedloop locationinventory problem under disruption risk,” Computers & Industrial Engineering, vol. 90, pp. 414–428, 2015. View at: Publisher Site  Google Scholar
 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
 H. Guo, C. Li, Y. Zhang, C. Zhang, and Y. Wang, “A Nonlinear Integer Programming Model for Integrated Location, Inventory, and Routing Decisions in a ClosedLoop Supply Chain,” Complexity, vol. 2018, pp. 1–17, 2018. View at: Publisher Site  Google Scholar
 H. Guo, C. Li, Y. Zhang, C. Zhang, and M. Lu, “A LocationInventory Problem in a ClosedLoop Supply Chain with Secondary Market Consideration,” Sustainability , vol. 10, no. 6, p. 1891, 2018. View at: Publisher Site  Google Scholar
 H. Guo, Y. Zhang, C. Zhang, Y. Liu, and Y. Zhou, “Locationinventory decisions for closedloop supply chain management in the presence of the secondary market,” Annals of Operations Research. View at: Publisher Site  Google Scholar
 M. S. Pishvaee, R. Z. Farahani, and W. Dullaert, “A memetic algorithm for biobjective integrated forward/reverse logistics network design,” Computers & Operations Research, vol. 37, no. 6, pp. 1100–1112, 2010. View at: Publisher Site  Google Scholar
 S. M. Hatefi and F. Jolai, “Robust and reliable forwardreverse logistics network design under demand uncertainty and facility disruptions,” Applied Mathematical Modelling, vol. 38, no. 910, pp. 2630–2647, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 S. Axsäter, “Using the deterministic EOQ formula in stochastic inventory control,” Management Science, vol. 42, no. 6, pp. 830–834, 1996. View at: Publisher Site  Google Scholar
 G. D. Eppen, “Note—Effects of Centralization on Expected Costs in a MultiLocation Newsboy Problem,” Management Science, vol. 25, no. 5, pp. 498–501, 1979. View at: Publisher Site  Google Scholar
 Z.J. M. Shen, “Integrated supply chain design models: a survey and future research directions,” Journal of Industrial and Management Optimization, vol. 3, no. 1, pp. 1–27, 2007. View at: Publisher Site  Google Scholar
 R. Storn and K. Price, “Differential evolution—a simple and efficient heuristic for global optimization over continuous spaces,” Journal of Global Optimization, vol. 11, no. 4, pp. 341–359, 1997. View at: Publisher Site  Google Scholar
 B. Alatas, E. Akin, and A. Karci, “MODENAR: Multiobjective differential evolution algorithm for mining numeric association rules,” Applied Soft Computing, vol. 8, no. 1, pp. 646–656, 2008. View at: Publisher Site  Google Scholar
 R. S. Rahnamayan, H. R. Tizhoosh, and M. M. A. Salama, “Oppositionbased differential evolution,” IEEE Transactions on Evolutionary Computation, vol. 12, no. 1, pp. 64–79, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Congdong Li 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.