Abstract
Recycling waste products is an environmentalfriendly activity that can result in manufacturing cost saving and economic efficiency improving. In the beer industry, recycling bottles can reduce manufacturing cost and the industry’s carbon footprint. This paper presents a model for a collectiondistribution center location and allocation problem in a closedloop supply chain for the beer industry under a fuzzy random environment, in which the objectives are to minimize total costs and transportation pollution. Both random and fuzzy uncertainties, for which return rate and disposal rate are considered fuzzy random variables, are jointly handled in this paper to ensure a more practical problem solution. A heuristic algorithm based on prioritybased globallocalneighbor particle swarm optimization (pbglnPSO) is applied to ensure reliable solutions for this NPhard problem. A beer company case study is given to illustrate the application of the proposed model and to demonstrate the prioritybased globallocalneighbor particle swarm optimization.
1. Introduction
Due to resource scarcity and environmental concerns, responsible companies are beginning to pay attention to the future of the planet and the global environment. Recycling used products for remanufacturing is, therefore, becoming of greater importance in supply chain management, a move that can dramatically reduce carbon emissions [1]. Closedloop supply chain (CLSC) combines the forward supply chain with a reverse supply chain to cover the whole product life cycle [2], with the manufacturing of new products and the transportation to customers via distribution centers and retailers as the forward supply chain and recycling, sorting, disposal, and remanufacturing as the reverse supply chain. In recent years, the CLSC has received a great deal of academic and business attention because of the need to be socially responsible, global environmental concerns, and government legislation [3, 4], all of which have motivated companies to pay more attention to recycling to reduce costs and lessen their carbon footprint.
Facility location and allocation problems (FLAPs) have been widely studied. Subramanian et al. [5] developed prioritybased simulated annealing to solve a CLSC network design problem, in which the distribution center (DC) and the centralized return collection center (CC) were set. Amin and Zhang [6] presented facilities location model for manufacturing and remanufacturing plants and CLSC collection centers, which included demand and return uncertainties. Subulan et al. [7] developed a CLSC network design model for the lead/acid battery industry that considered both financial and collection objectives. CLSC network design in a competitive environment with pricedependent demand was examined by Rezapour et al. [4], in which the DC and CC were separately built. Zeballos et al. [8] proposed a model for a multiperiod CLSC design and planning problem with demand uncertainty that had ten echelons in which the DC and CC were considered. Wang et al. [9] developed a granular robust model for a twostage wastetoenergy feedstock flow planning problem with uncertain capacity expansion costs. Tokhmehchi et al. [10] developed a hybrid approach to solve a closedloop supply chain location and allocation problem that considered the minimization of total cost. Vahdani and Mohammadi [11] proposed capacitated bidirectional facilities for CLSC conduct distribution, in which a multipriority queuing system was studied. As a growing number of companies are now engaging in recycling activities due to economic and environmental concerns, distribution and collection activities using the same vehicle have been found to reduce carbon emissions and transportation costs as empty loads can be avoided. In this paper, to benefit company operation and reduce construction costs, a distribution center (DC) is combined with a collection center (CC) into a collectiondistribution center (CDC). In practice, as the recycled product owners are usually at the same location as the potential new product buyer [12], a DC/CC combination has lower construction and operating expenses and can significantly reduce environmental pollution.
Ramkumar et al. [13] developed a multiechelon, multiperiod, multiproduct closedloop supply chain network model which was solved using a genetic algorithm with fixed variables. Kaya and Urek [12] presented a facility locationinventorypricing model without uncertainty to determine optimal facilities locations. Barz [14] proposed an optimization model for a twostage capacitated facility location and allocation problem with additive manufacturing, in which all variables were certain. Jindal and Sangwan [15] developed a multiobjective model for a CLSC network design problem with the economic and environmental factors being fuzzy uncertain variable and the DC and CC were separate. Ramezani et al. [16] conducted research into a CLSC network design problem that only considered of fuzzy variables. In recent years, uncertainty has attracted more research attention [17–19]. Stochastic programing, robust optimization, and fuzzy set theory have been used to present uncertainty in FLAPs [20, 21]. Wang et al. [22] used prediction sets to solve an expansion planning problem for wastetoenergy (WtE) systems facing future waste supply uncertainty. Keyvanshokooh et al. [23] proposed a novel hybrid robuststochastic programing (HRSP) approach to simultaneously model two different types of uncertainties by using stochastic scenarios for the transportation costs and polyhedral uncertainty sets for the demand and returns. However, the DC and the CC were separate and the collection disposal rate was a certain variable.
Uncertainties exist in both forward and reverse supply chains. However, the uncertainties in the reverse flow are higher than in the forward supply chain [7, 24, 25] as returned product quantity is generally seen as uncertain [23, 26]. Subjective uncertainties such as the decision maker’s choices and environmental coefficients can be dealt with using fuzziness, while objective uncertainties such as unit transportation costs, product prices, and the quantity of unusable products can be dealt with using randomness. In this paper, to reflect the study problem, the return rate and disposal rate are considered fuzzy random variables, which are concurrently handled using triangular fuzzy numbers [7]. Based on the above, a model is formulated to determine the optimum CDC number and location and the allocation strategies for the different facility types.
As facilities location and allocation problems are seen as nonconvex, nondifferentiable, strongly NPhard problems, a collection and distribution center location and allocation problem (CDCLAP) in a closedloop supply chain under a fuzzy random environment is even more complicated. Several different methods have been used to solve NPhard problems [27–29]. While particle swarm optimization (PSO) has been found to be generally effective [30–32], when the local optimal solution is found, the particle behavior in a basic PSO is directly influenced and therefore frequently falls into a local optimum [33–35]. Because of this problem, several advanced PSOs have been developed to more accurately solve supply chain management problems. Ai and Kachitvichyanukul [36] proposed a globallocalneighbor PSO which was found to be more effective. Based on this innovation, Xu et al. [33] proposed a fuzzy random simulationbased bilevel globallocalneighbor particle swarm optimization (frsbglnPSO). In this paper, a prioritybased globallocalneighbor particle swarm optimization (pbglnPSO) is applied to solve the CDCLAP.
In summary, this paper proposes a mathematical model to solve a collectiondistribution center location and allocation problem in a closedloop supply chain that considers economic and environmental factors and includes fuzzy random variables for return and disposal rates. The remainder of this paper is organized as follows. Section 2 presents the problem statement and model assumptions, after which a description of the model and its formulations is given in Section 3. The proposed hybrid solution based on the developed pbglnPSO is described in Section 4 and case study is presented in Section 5 to illustrate model formulation and the proposed method. Finally, Section 6 gives conclusions and indications for future research extensions.
2. Research Problem Statement
In this paper, a company with factories in certain locations and retailers in different customer zones is considered. The company needs to decide the locations for their integrated collection and distribution centers (CDCs), at which both a used product collection network and a new product distribution network are to be jointly [12]. As CDCs reduce construction and transportation costs because the same vehicles are used for both distribution and recycling, in this paper, only CDCs are considered.
A general illustration of the classical CDCLAP for a closedloop supply chain is shown in Figure 1, with the CLSC framework shown in Loop 1. The CLSC framework has four echelons: factories, CDCs, retailers, and disposal centers [11]. The forward supply chain begins with new production, after which the finished products are transported from the factories to the retailers via the CDCs. In the reverse supply chain, returned products are collected and transported to the CDCs, where the recycled products are inspected, consolidated, and sorted into those that are available for remanufacturing, which are sent to the factories, and those that are unsuitable for remanufacturing, which are transported to the disposal centers [23]. A CDC supplies products to multiple retailers; however, retailer demand is fulfilled by only one production site. A CDC can also handle products from different factories and dispatch returned products to multiple factories for remanufacturing.
In the CLSC examined in this paper, the retailers’ demand is estimated based on preorders. However, the return rate is considered to be a fuzzy random variable as customers may not return the used product or the product may be broken. Consequently, the availability of recycled products is unsure because of unsure transportation and carrying losses. Another fuzzy random variable considered in this paper is the returned product disposal rate, which is dependent on the inspection and consolidation at the CDC.
The assumptions for the proposed problem investigation are as follows: (1) only one product in one period is considered; (2) all alternative CDC locations have been identified; (3) recycling a used product costs less than manufacturing a new one; (4) the CDCs and factories have capacity limits [37–39]. As considering incapacitated facilities is an unrealistic assumption in many LAPs, many researchers have assigned a maximum capacity level to facilities to model more realistic decisions; (5) the factories’, retailers’, and disposal centers’ locations are known; (6) new product and returned product storage are allowed at the CDCs [21].
The initial problem is deciding which CDCs to select from the candidate sites and which allocation strategies to select to minimize total CDC costs: operating costs, transportation costs, and transportation pollution costs, while also considering flow constraints, capacity limits, and retailer demand.
3. Modelling
In this section, the mathematical formulations are given for the CDCLAP in the CLSC and the notations are given in the Notations to facilitate the problem description.
3.1. Objective Functions
Based on the variables mentioned in the Notations, the objectives are to minimize total costs and to minimize the environmental effects with the primary objective being minimizing total cost.
3.1.1. Economic Objective
In general, decision makers seek to minimize total costs, which are made up of transportation costs, fixed costs, and operating costs. The minimization objective can be described as
Equation (1) calculates the total cost, in which is the new product transport costs from the factories to the CDC, is the transport costs between the CDCs and the retailers, and is the returned product delivery costs from the CDCs to the disposal centers. The returned product transportation costs from the CDCs to disposal centers are measured as . The fixed costs for opening a new CDC are . denotes the new product variable costs. calculates the returned product operating costs.
As it is very difficult to handle objective functions with fuzzy random factors, Kruse and Meyer [40] demonstrated that the fuzzy expected value could be represented by a single fuzzy number. Based on the theory proposed by Heilpern [41], without a loss of generality, the expected value operator is used to convert the uncertain model into a deterministic model, which can then be used to transform the fuzzy random objective functions into their crisp equivalences, as shown in
Note that or above represents two expected values: the first is used to convert the fuzzy random variables into fuzzy numbers based on Kruse and Meyer’s 1987 theory, and the second is used to transform the fuzzy numbers into deterministic numbers based on Heilpern’s 1992 theory.
3.1.2. Environmental Objective
The secondary objective is to minimize the transportation carbon emissions associated with the CLSC operations, an area that has attracted significant recent research attention [42]. The following expression represents the transportation carbon emissions between the CDCs and the factories, the CDCs and the retailers, and the CDCs and the disposal centers.
is the environmental pollution caused by the transportation activities from the factories to the CDCs. is the summation of the carbon footprints for transporting products between the CDCs and retailers. is the total carbon footprint from the CDCs to the disposal centers, and is the carbon footprint from the CDCs to factories for returned products.
3.2. Constraints
Note that as each CDC has its own capacity limit, it is unable to service goods beyond capacity; therefore, the capacity limit restriction can be written as follows: is a fuzzy random variable indicating the used product return rate transported from the retailer to CDC . indicates the product quantity from CDC to retailer . indicates the product quantity transported from factory to CDC . is the capacity of CDC .
Within the capacity constraint, the factory is able to manufacture new products to meet retailer needs as well as deal with the returned products from the CDCs. is the demand from retailer and is the capacity of factory . calculates the returned product quantity transported to factory for remanufacturing.
All products for the retailers are processed through the CDCs. The recycled product quantity is always less than the product quantity transported from the factory to the CDC, which can be described as follows: is a variable indicating the new product quantity transported from factory to CDC . refers to the returned product quantity transported from retailer to CDC .
The products provided to the retailers should meet their demand. is the product quantity CDC sends to retailer . The stochastic variable is the retailer ’s demand based on order quantity.
The returned product quantity transported to the CDCs is more than the product quantity transported to disposal centers. is the expression for the returned product quantity from retailer to CDC , and is the product quantity transported from CDC to retailer .
There should be at least one CDC and there should be no more CDCs than the specified upper limit. is the upper limit for the number of CDCs, which is dependent on demand, returned product quantity, and fixed capacity constraints.
Each retailer must be serviced by only one CDC.
Since and are binary variables, the following constraints are required:
is a binary variable indicating whether a CDC is opened at point . If location is chosen to open a CDC, then ; otherwise, . is a binary variable indicating whether retailer is serviced by CDC . If , then retailer is serviced by CDC ; otherwise, .
3.3. Global Model
From the formulation above, the model for the CDCLAP with capacity, flow, and quantity constraints is developed with the aims of minimizing total costs and total transportation pollution with the primary objective of minimizing the total cost. In the CLSC, both new and returned products are considered. The product can be reproduced to save raw materials and reduce waste and pollution. In our model, all costs involved in the CDCLAP are considered as well as the influence of the transportation activity pollution. Fuzzy random theory is used to deal with the real world complex uncertainties and ensure more scientific decisions. Therefore, this CDC situation is closer to the real situation as it can deal with complicated practical problems. Finally, the global model is given:
3.4. Model Transformation
The two objectives are believed to be some similar ones to some extent. The minimization of the environmental objective requires low transportation carbon emissions which could be similar to the transportation costs of the economic objective, which are mainly affected by transportation distance. Based on previous research [33, 43, 44], the secondary environmental impact reduction objective can be transformed into a constraint to reduce the complexity of the multipleobjective model (13). By the acceptable carbon emissions level, the decision makers are concerned much more about the total cost, namely, the primary objective. Suppose that there is a maximum average environmental impact level which is acceptable to the decision maker; therefore, the secondary objective can be transformed to be a constraint as follows: Namely,
As the economic costs for the primary objective and the secondary environmental objective are transformed into a constraint, the global model can be transformed into its corresponding equivalent model as follows:
4. The Heuristic Algorithms Based on PglnPSO
Particle swarm optimization (PSO) is an evolutionary algorithm which simulates social behavior such as birds flocking and fish schooling [45]. Using a fixed population of individuals, the PSO searches the feasible zone to seek solutions, which are then updated to achieve an optimal solution. The particles [46], characterized by their position and velocity, are decided on by their flying experience, their discoveries, or the discoveries of their companions. They fly through the problem spaces following the currently optimum particles to find the best solution between the populations and the best solution for each population. Even though the PSO has been widely used to solve NPhard problems [45, 47], in the basic PSO, as the particles in the swarm are weak, they tend to cluster rapidly toward the global best particle [31]. However, the globallocalneighbor particle swarm optimization (glnPSO) developed by Ai and Kachitvichyanukul [36] has been found to improve the weakness in the basic PSO. Xu et al. [43] proposed an even more advanced globallocalneighbor particle swarm optimization with exchangeable particles (GLNPSOep). In this section, a prioritybased globallocalneighbor particle swarm optimization (pbglnPSO) is proposed to solve the CDCLAP in the CLSC.
4.1. Notations for the pbglnPSO
The basic elements of the PSO are particles, population, velocity, inertia weight, individual best, and global best. The notations needed for the pbglnPSO are shown in the Notations.
4.2. Encoding and Decoding Algorithm
The decoding process is based on the prioritybased encoding developed by Gen and Cheng and the prioritybased decoding and encoding proposed by Gen et al. [48]. The priorities for the CDCs and the retailers are equal to the total number of retailers and CDCs. At each step, the CDC (retailer) with the highest priority is selected and connected to a retailer (CDC) under a minimum transportation cost constraint. Procedure 1 shows the decoding algorithm for the prioritybased encoding and its trace table, with the prioritybased encoding being random. The CDCLAP is solved in two stages [14]. In the first stage, the CDC location is chosen and the transportation between the CDCs and retailers calculated. In the second stage, the allocations between the factories and CDCs are dealt with.

4.3. The Fitness Value Function
The CDCLAP considered in this paper has a primary objective, the minimization of total costs, and a secondary objective, the minimization of environmental pollution, with the fitness value of each particle reflecting the objective value of the total cost. As the second objective has been transformed into a constraint, when the transportation carbon emissions are beyond the acceptable level [33, 43], a penalty function is assigned according to the actual situation. If the transportation carbon emissions are within the acceptable level , the fitness value function in the pbglnPSO is as follows:
If the transportation carbon emissions are beyond the acceptable level, a penalty function is assigned and the fitness function isin which is a large enough penalty factor.
4.4. Update
Based on the notations of pbglnPSO mentioned in the Notations and the glnPSO proposed by Ai and Kachitvichyanukul [36], the inertia weight, velocity, and position are updated using the following equation:The glnPSO has been widely used in solving NPhard facilities location and allocation problems.
4.5. Overall Process of the pbglnPSO
In this paper, the glnPSO presented in Procedure 1 is used to solve a location and allocation problem. Due to uncertainties and environmental changes, a prioritybased globallocalneighbor particle swarm optimization (pbglnPSO) is proposed to solve this model. As the company pays close attention to economic costs, the environmental factor is dealt with as a constraint that has upper limits. The algorithmic details are as follows.
Step 1. Initialize particles as a swarm: (the particle is the priority).
Step 2 (constraints check). If in the feasible region, go to Step 3; otherwise, return to Step 1.
Step 3. Calculate the fitness according to the decoding algorithm in Procedure 1.
Step 4. Update the particle positions and velocities.
Step 4.1. Acquire the expected value for from the above algorithm.
Step 4.2. For , decode each particle to an installment group. Calculate the fitness value of each particle and set the position of the th particle as its personal best. The global best position is chosen from these personal best positions.
Step 4.3 (update best). For , if , .
Step 4.4 (update best). For , if , .
Step 4.5 (update best). For , among all best of neighbors around the th particle, set the personal best which has the best fitness value as .
Step 4.6 (generate best). For and , find ensuring that the FDR takes a maximum value, and set as .
Step 4.7. Update the position and the velocity of each th particle using (20) and (21).
Step 4.8. Check whether the particles are beyond the mark. If , ; otherwise, if , then .
Step 5. Based on the above calculation, replace the ranking vector using the new numbers.
Step 6. If the stopping criterion is met, stop; otherwise, and return to Step 2.
The overall process can be clearly seen in Figure 2.
5. Case Study
5.1. Case Presentation
This model is based on a beer company in a developing country that bottles beer in glass bottles. The supply chain allows customers to return empty bottles to the retailers, which are then sent to the CDCs where they are inspected, consolidated, and sorted. After processing and disinfecting, the bottles are refilled and sold again. The company is now considering the construction of several CDCs to allow for bottle recycling as producing new bottles is far more expensive than recycling used bottles.
To illustrate the validity of the model and the usefulness of the solution method, the data needed to examine the CLSC performance for the four echelons is presented here. Based on the market analysis, ten alternative CDC coordinates are suggested, which are to be assessed based on location, capacity, fixed costs, new product variable costs (NPV cost), and recycled product variable costs (RPV cost). The assessments for the ten possible CDCs are shown in Table 1. Supermarkets and restaurants are considered to be the beer retailers with flexible demand. Table 2 presents the information regarding the retailers, factories, and disposal centers. It can be seen from Table 2 that to represent 30 different retailers, to are the 4 different factories in different locations, each of which has a different capacity, and indicates the location and capacity of the disposal center. Therefore, 30 retailers, 4 factories, and 1 waste disposal center are considered in this study. The unit transportation costs and pollution are related to the distances between the facilities. The retailers’ return rates, which are fuzzy random variables, are shown in Table 3.
5.2. Maximum Generation and Population Size
With an increase in the maximum iterations, computer time may increase or the optimal result may improve. For the pbglnPSO in this case study, the maximum generation is set at and the population size is set at ; therefore, the maximum iteration is . In the test, is set from 20 to 60 with a steplength of 10, and is set from 100 to 500 with a steplength of 100, while , from which 25 different groups are obtained, and then pbglnPSO was run 20 times for each group.
Figures 3 and 4 show the test results, the computing time, and the optimal results. For the horizontal ordinate (e.g., “0–5” represents five different groups) [49], when , increases from 100 to 500 with a steplength of 100, with the remainder following the same analogy. From Figure 3, it can be seen that the maximum iterations significantly influenced the computing time. When the population size was the same, there was a positive correlation between the maximum generation and the computing time. For the optimal result, the maximum generation had an obvious influence on the result, as shown in Figure 4. As can be seen, the results were better when was from 300 to 500 and was from 40 to 60. Thus, for the above and group, the optimal results, the standard deviation of 20 optimal results, and the average computing time have been given in Table 4 for better analysis. From Table 4, it can be seen that the result was stable and reliable and that any further increase in the maximum iterations resulted in a longer computing time but no improvement in the optimal result, with the best result being when = 50 and = 400. Therefore, the suitable values for and in the pbglnPSO for this case were found to be 400 and 50.
5.3. Sensitivity Analysis on the Parameters
To find the best solution to the proposed model, a series of experiments were conducted, all of which were performed using MATLAB 7.0 on a workstation with an Intel(R) Core(TM)i7, a 2.59 GHz clock pulse with 8 GB memory, and a Windows 10 operating system. A sensitivity analysis was performed to examine the effectiveness and behavior of the proposed algorithm, as shown in Table 5. Several parameters were changed: the population size , the maximum generation , and the acceleration constants , , , and . After trying various values for the population size and maximum generations, the results were found to be better when was from 300 to 500 and was from 40 to 60. The different fitness values obtained using the pbglnPSO with the different parameters , , , and are shown in Table 5.
As can be seen from Table 5, when the parameters , , , and increased, the fitness value improved except for with the same generation and population size, with the fitness value increased from to . Therefore, when , the result was found to be optimal. For , given the same , , , and and population size, it was found that when was 400, the fitness value was better than for any other generation. Finally, for , the results improved as the population size increased and it was found to be optimal when is 50. The most effective and efficient results were gained with at 400, at 50, and .
5.4. Result Analysis
In this section, the pbglnPSO is conducted to solve the model using the above data. The parameters for the problem were set as follows: population size: = 50; maximum generation: = 400; inertia weight: and ; acceleration constant: . After running the program 20 times, the most satisfactory solution was found. Figure 6 shows the specific objective values found by the pbglnPSO in different iterations and shows the reductions in the total costs. The results are presented in Figure 5. From the calculations, at least 3 CDCs could satisfy all markets. Alternative CDC positions 4, 5, and 8 should be chosen. CDC 4 can send products to markets 2, 3, 4, 6, 7, 9, 14, 15, 25, and 30. CDC 5 can transport products for 1, 8, 11, 12, 13, 16, 17, 18, 19, 20, 21, 22, and 27 and retailers 5, 10, 23, 24, 26, 28, and 29 can be serviced by CDC 8. The total cost was found to be 8.192 million RMB, in which the fixed costs were 4.4 million RMB, the transportation costs were 3.776 million RMB, and the operating costs were 0.016 million RMB.
5.5. Algorithm Comparison
To better illustrate the effectiveness of the proposed algorithm, a brief comparison between the pbglnPSO, glnPSO, and an immune algorithm (IM) is given in this section. The glnPSO is a wellrespected evolutionary algorithm that has been successfully implemented in a variety of engineering and combinatorial problems. The IM is also being widely used to solve facilities location problems.
To establish the solution quality for the pbglnPSO, it was compared with the glnPSO and the IM. Each run time for the pbglnPSO, glnPSO, and the IM was around 10 s. The pbglnPSO, glnPSO, and IM were run 20 times using the same data. For a fair comparison between the groups, the population size was set at 50 and the maximum generation at 400. In the glnPSO, the acceleration constant was designed as and the inertia weight was and . For the IM algorithm, the crossover probability was 1 and the mutation probability was 0.1.
From Figure 7, it can be seen that the pbglnPSO outperformed both the glnPSO and the IM, and, as the glnPSO converged faster, it had a better result than the IM. This demonstrated that a better solution can be obtained using the glnPSO, and an even better result can be obtained using the pbglnPSO. The blue profile shows the convergence for the best in history for the pbglnPSO. It can be seen from Figure 7 that, as the programs ran, the results become stable for the pbglnPSO and glnPSO after about the 150th generation, while the IM became stable after the 175th generation. As shown in Figure 7, the best solution for the pbglnPSO was superior to, more stable than, and had the smallest CPU run time compared to the other algorithms (Table 6), with the IM having the highest run time.
6. Conclusion
Economic development has resulted in many environmental pollution problems, the seriousness of which has encouraged people to recycle and reuse products. To examine this problem and seek appropriate solutions, a mathematical model for a collectiondistribution center location and allocation problem in a closedloop supply chain under a fuzzy random environment was presented for the beer industry in China. For this problem, a new model was formulated, in which the decision makers sought to minimize costs and pollution under capacity and quantity constraints. To more accurately represent actual production situations, the return rate and disposal rate were considered to be fuzzy random variables. A heuristic algorithm, the pbglnPSO, was then applied to solve the problem. Based on the proposed priority, the distribution and collection activity were shown to satisfy retailer demand and reduce costs and pollution after the CDCs start operations. After calculation, the best solution was determined and the advantages of the algorithm were illustrated. The proposed model and method can be applied to the location and allocation of CDCs in the beer industry to improve supply chain management. The model was shown to assist in generating retailer demand and dealing with the returned products, which could benefit company recycling and reuse policies. At the same time, the transportation costs and pollution were reduced because of the reductions in losses from empty loads.
Notations
Sets:  Set of CDCs, 
:  Set of factories, 
:  Set of retailers, 
:  Set of disposal centers, 
:  Alternative location position for the CDCs, 
:  Known position of the factories, 
:  Known position of the retailers, 
:  Known disposal center, 
:  The upper limit of the CDCs 
:  The demand of retailer 
:  The capability of CDC 
:  The capability of factory 
:  Product quantity from factory to CDC 
:  Product quantity from CDC to retailer 
:  The product return rate from retailer 
:  The product disposal rate at CDC 
:  The fixed costs of the CDC 
:  The variable costs of the CDC for a new product unit 
:  The variable cost of the CDC triage for a returned product unit 
:  Unit transportation cost between CDC and factory 
:  Unit transportation cost between CDC and retailer 
:  Unit transportation cost between CDC and disposal center 
:  Environmental impact of transportation between CDC and factory 
:  Environmental impact of transportation between CDC and retailer 
:  Environmental impact of transportation between CDC and disposal center 
:  The environmental impact level accepted by decision makers. 
:  A binary variable indicating whether point is chosen. If point is chosen, then ; else, 
:  It indicates whether retailer is served by CDC . If is chosen, then ; else, 
:  Iteration index, 
:  Particle index, 
:  Velocity of the th particle at the th dimension in the th iteration 
:  Personal best position 
:  Local best position 
:  Personal best position acceleration constant 
:  Local best position acceleration constant 
:  Maximum position value 
:  Velocity vector of th particle 
:  Vector personal best position of th particle 
:  Vector local best position of th particle 
:  Fitness value of 
:  Dimension index, 
:  Inertia weight in th iteration 
:  Position of the th particle at the th dimension in the th iteration 
:  Global best position 
:  Nearneighbor best position 
:  Global best position acceleration constant 
:  Nearneighbor best position acceleration constant 
:  Minimum position value 
:  Position vector of th particle 
:  Vector global personal best position 
:  Uniform distributed random number within . 
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by National Natural Science Foundation of China (Grants no. 71640013 and no. 71501137), National Planning Office of Philosophy and Social Science (Grant no. 14BGL055), and System Science and Enterprise Development Research Center (Grant no. Xq16C10).