Research Article  Open Access
Optimizing the Joint Replenishment and Delivery Scheduling Problem under Fuzzy Environment Using Inverse Weight Fuzzy Nonlinear Programming Method
Abstract
In reality, decisionmakers are always in front of imprecise and vague operational conditions. We propose a practical multiobjective joint replenishment and delivery scheduling (JRD) model with deterministic demand and fuzzy cost. This model minimizes the total cost defuzzified by the signed distance method and maximizes the credibility that the total cost does not exceed the budget level. Then, an inverse weight fuzzy nonlinear programming (IWFNLP) method is adopted to formulate the proposed model. This method embeds the idea of inverse weights into the MaxMin fuzzy model. Thirdly, the fuzzy simulation approach and differential evolution algorithm (DE) are utilized to solve this problem. Results show that solutions derived from the IWFNLP method satisfy the decisionmakerās desirable achievement level of the cost objective and credibility objective. It is an effective decision tool since it can really reflect the relative importance of each fuzzy component. Our study also shows that the improved DE outperforms DE with a faster convergence speed.
1. Introduction
As a multiitem inventory problem, the joint replenishment problem (JRP) has been widely applied to lots of sizing problems in manufacturing applications (Hsu [1]; Goyal [2]; Wang et al. [3, 4]). Beside the possible quantity discount, gathering several items on a single order to reduce the total of these fixed ordering costs is pretty reasonable (Kaspi and Rosenblatt [5]). The studies of the JRPs can be divided into two categories: (a) heuristics for the classic JRP under constant demand and special applications (Olsen [6]; Lee and Yao [7]; Porras and Dekker [8]) and (b) the JRPs under dynamic and/or stochastic demand. The extensive literature review is available in Khouja and Goyal [9] and Narayanan et al. [10].
Typically speaking, if the warehouse is assumed as the central of a supply chain for all the new JRPs, two extensions of JRP should be noted; one extension is in the supplying end and the other extension is in the selling end. For both of two extensions, the delivery considerations should be considered. Here we call them the joint replenishment and delivery scheduling (JRD) problems. In fact, most corporations with global purchasing have realized that considerable cost savings can be achieved by a JRD policy (Sindhuchao et al. [11]). The relatively scarce literature on the JRDs can be classified into two categories: (a) JRDs with deterministic demand, Chan et al. [12] addressed scheduling issues for multibuyer joint replenishments and Cha et al. [13] studied a JRD model of the onewarehouse and retailer system and developed a genetic algorithm (GA) and an improved heuristic named RAND to solve this problem and Moon et al. [14] provided joint replenishment and consolidated freight delivery policies with deterministic demand for a third party warehouse; (b) JRDs with stochastic demand, Wang et al. [15] proposed an effective and efficient differential evolution algorithm (DE) for the integrated stochastic JRD. Similar studies can be found in [16ā18]. A limitation common in all the JRDs mentioned above is that all the key factors are assumed to be certain.
In reality, it is more realistic to handle imprecise values using the fuzzy theory for the JRD. In fact, decisionmakers are always in front of imprecise and vague operational conditions (Pishvaee and Torabi [19]; Zeng et al. [20]). Uncertainties have been tackled in a lot of ways and fuzzy set theory has a long history for managing inventory (Chiang et al. [21]; Kacprzyk and Stanieski [22]; Wang et al. [23, 24]; Wee et al. [25]). Generally, the imprecision may originate from two aspects in the JRD modeling: (1) the imprecise specification of objectives, for example, a decisionmaker may have to face vague goals such as āthis total cost should be around $20000ā; (2) the imprecise specification of related parameters, an important task involved in the fuzzy JRD model is to predict parametersā values such as the holding cost and major ordering cost. However, due to the nonavailability of sufficient and precise input data, the precise predicted values cannot be obtained easily; just an āapproximateā value may be ascertained, while fuzzy numbers can efficiently model the imprecise values.
At present, there are two common approaches to handle fuzzy parameters. (1) Defuzzification. As a favorite approach in many inventory modeling for its simplicity, defuzzification can easily transfer fuzziness to be explicit without complex analysis (Roy et al. [26]). Handfield et al. [27] defuzzified a (Q, r) model based on fuzzyset representations of various sources of uncertainly in the supply chain. (2) Fuzzy dependentchance programming (DCP). Fuzzy DCP maximizes the credibility of event such that the total cost in the planning periods does not exceed a certain budget level (Liu [28]). Wang et al. [29] studied a novel JRP model based on DCP with fuzzy minor replenishment cost and fuzzy inventory holding cost. Similar paper can be found in Peng and Zhou [30]. For the JRP decision, managers also hope that the total cost does not exceed the budget level, especially when enterprises have much cash flow pressure. In this situation, a fuzzy DCP model is always regarded as a good choice.
However, no study has simultaneously considered total cost and credibility of the total cost does not exceed the budget level as performance criterions for the JRD with fuzzy cost. This unexplored area is important and interesting since it integrates into a single model two main decision criterions: total cost and credibility. The aim of this paper is to develop a practical multiobjective JRD (MJRD) model with deterministic demand and fuzzy cost firstly. Moreover, an inverse weight fuzzy nonlinear programming (IWFNLP) adopted by [25] is applied to the MJRD to satisfy the decisionmakerās desirable achievement level of each fuzzy component. At last, novel hybrid algorithms are provided using the fuzzy simulation approach (FSA) and DE for this model handled by IWFNLP method. Results of examples show solutions derived from the IWFNLP method satisfy the decisionmakerās desirable achievement level of the cost and credibility objective. It is an effective decision tool since it can really reflect the relative importance of each fuzzy component.
The rest of this paper is organized as follows. The multiobjective JRD model is given in Section 2. In Section 3, the fuzzy simulation for credibility calculation and DE for searching for an optimum solution are introduced. Section 4 is numerical examples and analysis. Section 5 contains conclusions and future research directions.
2. Mathematical Model and Analysis
2.1. Assumptions and Formulation of Fuzzy JRD Model
(1) Assumptions and Notations. Lu and Posner [31] addressed the problem of supplying multiple retail outlets with constant and continuous demand from a single warehouse. Each retailer sells one product, and other retailers sell the same product. All demands must be met without shortages or backlogging. Orders placed by a retailer will result in demand generated at the warehouse. The objective is to minimize the average total cost. This one warehouse, nretailer JRD problem has received considerable attention from researchers [13, 14]. In this study, we focus on the JRD model under fuzzy environment. The following assumptions are similar with Lu and Posner [31] and Cha et al. [13]. (1) All parameters including demand rates and costs (except warehouseās major ordering cost) are known and constant; (2) replenishment is instantaneous; (3) replenishment lead time is constant.
As reported in Wang et al. [23], an interesting real problem is when human originated data like ordering cost are not precisely known but subjectively estimated or linguistically expressed because of the lack of the accurate history data. Suppose that the major ordering cost (S) is not precisely known. For example, let the linguistic estimates of S be as follows. āThe major ordering cost is about 40 dollars per order, but not less than 20 dollars per order nor more than 60 dollars per order.ā In reality, it is very hard to obtain the precise cost. So fuzzy variables are also utilized to handle the JRD problem under uncertainty. In this study, warehouseās major ordering cost is treated as a fuzzy number. Our work introduces fuzziness into the JRD which makes it become more practical. Due to the JRDās difficult mathematical properties, we just suppose the major ordering cost as a fuzzy number to make the study become more possible because two different objectives are considered simultaneously.
As presented in Cha et al. [13], the warehouse will deliver to an individual retailer after it replenishes goods jointly from the suppliers by taking into account the demand of each item at the retailer. The warehouse that has a lot of associated retailers can obtain significant cost savings (logistics cost and ordering cost) by replenishing jointly. The reduction effects will be higher if materials are imported from overseas with the high procurement cost. Referring to [13], the JRD model in this study can be described in Figure 1. Note the warehouseās major ordering cost is a fuzzy number.
In order to discuss the JRD problem, the following notations are defined:ā: number of items;ā: index of item, ;ā: demand rate of item ;ā: warehouseās major ordering cost, a triangle fuzzy number;ā: warehouseās minor ordering cost of item ;ā: warehouseās inventory holding cost of item per unit time;ā: warehouseās outbound transportation cost of item ;ā: retailerās inventory holding cost of item per unit per time;ā: warehouseās basic cycle time (decision variable);ā: integer number that decides the outbound schedule of item (decision variable);ā: integer number that decides the replenishment schedule of item (decision variable), it means the replenishment cycle of item is .
(2) Formulation with Fuzzy Major Ordering Cost. Similar to [13], we also discuss different types of productions and assume that all productions can be stored at the warehouse. Further, we assume that item is stored and sold only by retailer . The warehouse replenishes item at every integer multiple () of the basic cycle time () and delivers it to retailer .
Procedures to find optimal policies are very difficult. There are no known good approaches for solving this problem in time polynomial in the number of retailers (Lu and Posner [31]). Most scholars have concentrated on developing good heuristics for special policies. In this study, a popular stationary policy used in [13] is also adopted that a warehouse delivers item to retailer at the same time interval with a fixed order quantity.
The total cost (TC) is composed of the sum of the major ordering cost, minor ordering cost, inventory holding cost, and outbound transportation cost of a warehouse as well as the total of the inventory holding costs of retailers. According to the above definitions, the total relevant fuzzy cost per unit time to be minimized is given by
2.2. Proposed Fuzzy Multiobjective JRD (MJRD) Model and Analysis
2.2.1. Defuzzified Total Cost by Signed Distance
(1) Defuzzified TC. The signed distance method is simple and easy to handle, and this is why the extension principle and centroid method are not applied to this fuzzy JRD model. It is difficult for the extension principle and centroid to obtain the estimated TC in fuzzy sense. Moreover, Chiang et al. [21] found there was not a significant difference between signed distance method and the extension principle or centroid method. So the signed distance method is used to obtain the TC in fuzzy sense.
From (1), we can obtain the following equations:
From (2), the defuzzified TC to be minimized using signed distance can be obtained as
2.2.2. Credibility Object Based on DCP
Generally speaking, DCP is related to maximizing some chance functions events in an uncertain environment (Liu [28]). As a widely used type of stochastic programming, DCP has been extended to the area of fuzzy programming (R. Wang and L. Wang [32]; Wang et al. [33]). In practice, the goal of JRD policy cannot be confirmed exactly due to inevitable uncertainty. Hence, a realistic approach for decision makers (DMs) may be to maximize the credibility of achieving the optimization goals. Sometimes, DMs are not concerned with minimizing the total cost but hope that TC does not exceed the budget level TC_, especially under much cash flow pressure. In this situation, a natural idea is to maximize the credibility of the event such that the total cost is less than or equal toTC_. So an objective function can be written as
2.2.3. Proposed Fuzzy Multiobjective JRD (MJRD) Model
Due to the uncertainty of the decisionmaking, it is quite natural to assume two main goals: (1) nonrigid total cost goal; (2) credibility goal to assure the safety of cash flow. These goals represent different attitudes of managers for handling the inevitable uncertainty. In reality, it is not surprising that managers have their own opinions on the goal under uncertainty. For the fuzzy objectives, the MJRD model can be described as where the wavy bar āā denotes the fuzziness of the two objectives goal.
To describe the fuzzy objectives in (13), the decisionmaker may subjectively define the acceptable intervalāā for the defuzzification objective and for the credibility objective, where Ā· denotes .
Alternatively, the following four cases can be solved to confirm the values of , , , and , respectively.(1)The first one minimizes and the result is set as .(2)The second one minimizes and the result is set as .(3)The third one maximizes , and the result is set as .(4)The forth one maximizes , and the result is set as .
Then, the membership function is used to describe the attainable degree of the two objectives. In this study, the linear membership functions , are used, where
The pictorial representations of these membership functions are given in Figure 2.
(a)
(b)
2.3. Mathematical Analysis Using IWFNLP Method
2.3.1. Formulation Using Traditional Fuzzy Additive Goal Programming
In order to specify imprecise aspiration levels of the goals under fuzzy environment, Narasimhan [34] had firstly developed fuzzy goal programming (FGP) using membership functions. Chen and Tsai [35] reformulated the fuzzy additive goal programming (FAGP) by incorporating different important and preemptive priorities of fuzzy goals. In contrast to other methods, the FAGP allows managers to determine a desirable achievement degree for each fuzzy goal to reflect explicitly the relative importance of these goals. In order to better understand the following proposed MJRD model using IWFNLP, we give the formulation using traditional FAGP as follows: where and are weights for the defuzzification objective and credibility objective, respectively, and āā. These weights reflect the relative importance of each goal in the decision model.
2.3.2. Formulation Using IWFNLP Approach
(1) Rationales of Using IWFNLP. Although the FAGP can be applied to handle the fuzzy decision problem, Wee et al. [25] found that the ratio of the achievement levels (, ) obtained by this approach cannot really reflect the ratio of the weights. To overcome the shortages of the FAGP method, a novel IWFNLP method was adopted. This method embeds the idea of inverse weights into the MaxMin fuzzy model. Results of a lot of experiments show IWFNLP method can give solutions in which the ratio of the fuzzy components achievement levels is as close to the ratio of the assigned weights as possible. So the IWFNLP is also adapted to handle the MJRD. The effectiveness of this method is further verified in Section 4.
(2) Formulation. For the weights āā, onetoone transformation from the weight to the inverse weight is defined as
In this study, . So the fuzzy model is formulated as where , .
According to the definition of inverse weight, (12) is equivalent to
3. The Proposed Algorithm for Fuzzy MJRD Using FSA and DE
3.1. Outlet of the Proposed Algorithm
(1) The FSA Is Utilized to Solve the Fuzzy DCP. Usually, it is hard to obtain credibility value with an analytical method. Wang et al. [33] designed a FSA to calculate the credibility value for fixed decision variables. The advantage of the FSA is that it can calculate the credibility whether the membership functions of the fuzzy numbers are simple or not and whether several fuzzy numbers are involved or not. In fact, analytical methods usually can also be utilized for the same target. However, it can work only when member functions are simple, for example, linear functions and triangular functions. Moreover, analytical methods do not work well when too many fuzzy numbers are involved. But the FSA does not subject to such a restraint. So the FSA is adopted in this study.
(2) DE Is Adopted to Find the Solution Quickly. When the objective functions to be optimized are multimodal or the search spaces are particularly irregular, intelligent algorithms should be designed to solve the fuzzy DCP models. Moreover, optimization algorithms need to be highly robust in order to avoid getting stuck at the local optimal solution. Among these algorithms, genetic algorithm (GA) has been proved to be effective for the DCP models (Ke and Liu [36]). However, the GA displays inherent difficulties in performing local search for some applications because of the difficulty of the selection of the suitable probability of crossover and mutation. So it is important to find a novel algorithm to handle the fuzzy DCP more effectively.
DE is one of the best evolution algorithms (EAs) for solving nonlinear, nondifferentiable, and multimodal optimization problems (Storn and Price [37]). Due to its simple structure, easy implementation, quick convergence, and robustness, DE has been applied to a variety of fields [38ā43]. Wang et al. [4] found that the basic DE is a good candidate for the similar JRPs. But, the effectiveness of DEs for the fuzzy JRD should be studied further because of the difficult mathematical properties.
DE has a good performance in convergence speed, but the faster convergence may cause the diversity of population to descend quickly during the solution process [42]. A good tradeoff between convergence and diversity should be designed. Zou et al. [44] developed an improved DE (IDE) with a modified mutation factor according to the objective function values in mutation steps and adjusted crossover factor in terms of iteration number in crossover step. The IDE can not only diversify candidate solutions, but also increase the convergence rate. So the IDE is utilized in case of DE having unsatisfactory performance.
3.2. Fuzzy Simulation Approach for Solving DCPBased Model
Usually, it is difficult to compute with an analytic method. Instead, we use a FSA to obtain the value of for a fixed . According to [33], a classic fuzzy simulation approach can be adopted as follows.
Step 1. Set , , and . Set the maximum iteration times .
Step 2. Uniformly generate a sequence from such that Pos where is a sufficiently small number. Thus, a real vector can be obtained.
Step 3. Calculate and .
Step 4. If and , set ; otherwise, if and , set .
Step 5. Set and return to Step 2 until the maximum iteration times is reached.
Step 6. Return , where denotes Pos and denotes Pos or 1Nec.
3.3. Basic DE and Improved DE (IDE)
When the objective functions to be optimized are multimodal or the search spaces are particularly irregular, algorithms need to be highly robust in order to avoid getting stuck at local optimal solution. The advantage of DE is just to obtain the global optimal solution fairly. In the following, we discuss two DEs.
3.3.1. Basic DE
Basic DE consists of three evolution operators: mutation, crossover, and selection. In mutation operator, DE uses the differences between randomly selected individuals to generate a trial individual. Then crossover operator is used to produce one offspring which is only accepted if it improves on the fitness of parent individual. The process of choosing individuals is called selection. A brief description of the DE algorithm is as follows.
Step 1 (initialization and representation)
Initialization. The initial population is created by assigning random values which lies inside the feasible bound of the decision variable. Each individual is generated by where is the number of individuals; is the dimension of each individual; and represent the low and upper bound of the th decision parameter respectively; is a uniformly distributed random number in range [0,1].
Representation. Now, we will demonstrate how our chromosomes can be decoded to a feasible solution and how each chromosome of population is evaluated. For , Figure 3 shows what the chromosomes is composed of.
Assuming is the lower bound of th gene and is the upper bound of th gene, it is obviously that is 1 for and and is 0 for . It is difficult to give the maximum values for decision variables of the fuzzy MJRD model because and influence each other. We consider that RAND algorithm mentioned in Section 3.4 can obtain the optimum solution for the crisp JRD model. So we set as values of three times the optimal upper bounds obtained by the RAND_JRD for this fuzzy JRD model following the similar experience of Cha et al. [13]; Wang et al. [15].
Step 2 (mutation). For each target individual āā () of generation , a mutant individual is generated according to with randomly chosen integer indexes . Note that indexes have to be different from each other and from the running index. is called mutation factor between [0, 1] which controls the amplification of the differential variation .
Step 3 (crossover). The basic DE crossover operator implements a discrete recombination of the trial individual and the parent individual to produce offspring . The crossover is implemented as follows: where refers to the th element of the individual . and are similarly defined. is the th evaluation of a uniform random number generator between [0,1]. is a randomly chosen index from which ensures that gets at least one parameter from . Otherwise, no new parent individual would be produced and the population would not alter. CR is the crossover or recombination rate between [0, 1] which has to be determined by the user.
Step 4 (selection). The selection in DE is deterministic and simple. The evaluation function of an offspring is onetoone competition in the DE. It means the resulting trial individual will only replace the original if it has a lower objective function value. Otherwise, the parent will remain in the next generation if ,āāand then where is obtained by (13). This is the same for all variants of the DE. Although the selection pressure is only one, the best individual of the next generation will be at least as fit as the best individual of the current generation.
Step 5 (stop and output results). When stopping criterion is met, output the optimal results; otherwise, repeat Step 2āStep 4. In this study, stopping criterion is met when maximum number of iteration () is reached.
3.3.2. The IDE
The key difference between the IDE and basic DE is in the way of adjusting scale mutation factor and crossover rate CR. The method modifies mutation factor according to the objective function values of all candidate solutions in mutation step and adjusts crossover rate CR in terms of iteration number in crossover step. Both modified operators can not only diversify candidate solutions, but also increase the convergence of algorithm. In short, the IDE and DE are different in two aspects.
(1) For mutation of DE, it is set to a fixed value for all candidate solutions over all iterations. That is to say, all candidate solutions have the same magnification factor of the differential variation . Here, if the objective function is minimizing the total cost, an adaptive scale factor for the th candidate solution is stated as follows:
If the objective function is maximizing the credibility, an adaptive scale factor for the th candidate solution is stated as follows: where belongs to a uniform distribution in the ranges [0, 1]; represents the objective value of the th solution; and represent the minimal and maximal objective function values of all candidate solutions, respectively; represents the average objective function value of all solutions.
(2) For crossover rate CR of DE, it is set to a fixed value for any dimension of any candidate solution over all iterations. In other words, any dimension of any candidate solution has the same crossover rate, which does not change with the evolution process. Here, a dynamic crossover is adopted for all candidate solutions at iteration , and it is stated at follows: where and are current iteration number and maximal iteration number, respectively; and represent minimal crossover rate and maximal crossover, respectively.
3.4. Modified Algorithm for the Defuzzified JRD Model
(1) Optimal T. For a given set of and , the optimal basic cycle time can be easily obtained as shown in (23) from the first order derivative of the defuzzified TC function since the TC function is convex in :
(2) Optimal and . According to experiences of [13ā15, 17], for a given set of and , the optimality condition of can be derived from the following two conditions:
Therefore, the optimality condition of is defined as follows:
Similarly, the optimality condition of is defined as follows:
(3) RAND_JRD: A Modified Algorithm for the Defuzzified JRD. Cha et al. [13] proposed a modified RAND algorithm for the crisp JRD model to minimize the TC. Making corresponding changes in , an algorithm named RAND_JRD can be used to find the optimal solution. The procedures are as follows:
Step 1. Compute = and = āāfor each .
Step 2. Divide the range into different equally spaced values of . The value of is decided by the decisionmaker. Set .
Step 3. Set and . Put āā and .
Step 4. Set .
Step 5. For a given value of āā and a given set of āā , find the optimal values of using (25). Put .
Step 6. For a given value of āā and a given set of , find the optimal values of using (26). Put .
Step 7. For a given set of s and s, find the optimal value of using (23). Put .
Step 8. If , go to Step 4. Otherwise, put ,āā,ā andā. Compute for this .
Step 9. If , go to Step 3.
Otherwise, stop and select with the minimum .
3.5. Flowchart of Proposed Algorithm for Fuzzy MJRD Model
Main flowchart of the proposed algorithm is shown in Figure 4. In the adjustment, if or , then , where .
4. Numerical Examples
4.1. Experiment 1
Basic Experiment 1 is given to compare the DE and IDE. According to the recommendation of the inventor of DE [37] and the similar experiences of [4, 15, 17], the factors setting for DE used in all experiments are listed as follows: the population , the mutation factor , the crossover factor , and the maximum generation is set to 500. Factors of IDE are set as follows: ,ā ā. The data and results of this experiment are shown in Tables 1 and 2. Convergent curves of DE and IDE are given in Figure 5.


(a) Convergent curve of DE
(b) Convergent curve of IDE
We can calculate , , , and .
Then, DE and IDE are performed 50 times and the results are reported in Table 3.

From Table 2, the radios of the achievement levels gained by DE and IDE can be computed as follows:
Results in Tables 2 and 3 show (1) the ratios of the achievement levels and the weights for the fuzzy objectives are nearly equivalent using IWFNLP and DE/IDE; (2) IDE outperforms DE with a rapid convergence speed.
4.2. Experiment 2 with Different Major Ordering Cost
Khouja and Goyal [9] indicated that the value of major ordering cost affected the performance of algorithms. So three examples were designed to test the IDE. The data are same for , , , , , , and as shown in Table 2. But the values of are different as shown in Table 4. The results are listed in Table 5. Figure 6 shows the convergent processes.


(a) Experiment 2a
(b) Experiment 2b
(c) Experiment 2c
From Table 5 and Figure 6, we can conclude that the IDE can also find satisfactory solutions with different values of major ordering cost. The ratios of the achievement levels and the weights for the fuzzy objectives are still nearly equivalent using IWFNLP method and IDE.
5. Conclusions and Future Research
This paper is an interdisciplinary research of the fuzzy inventory model and intelligent optimization algorithm. Due to the inevitable uncertainty, it is quite natural for decisionmaker to assume two main goals: (1) nonrigid total cost goal; (2) credibility goal to assure the safety of cash flow. We developed a practical JRD model under uncertainty and provided an effective algorithm for this model. The main contributions are as follows.
(1) Actually, there are lots of papers discussed inventory and risk management issues [45]. However, to our best knowledge, the JRD model under fuzzy costs that simultaneously minimizes the total cost and maximizes the credibility to assure the safety of cash flow is nonexistence. Our work provides a useful approach for the joint replenishment and delivery scheduling under uncertainty.
(2) The formulation of the proposed MJRD model is handled by the IWFNLP which can make the ratios of the achievement levels of objectives and the weights for the fuzzy objectives are nearly equivalent. The IWFNLP method gives the solution that satisfies the decisionmakerās desirable achievement level of the total cost objective and credibility objective. It is an effective decision tool to ensure a decisionmakerās expectation is achieved.
(3) Hybrid intelligent algorithms are designed to solve the proposed JRD handled by the IWFNLP method using the FSA and DE/IDE. Results of numerical examples show the IDE can find satisfactory solutions faster than DE.
Other intelligent algorithms, such as geneticsimulated annealing algorithm [46] and quantum evolution algorithm [47], also show good performances to solve complex optimization problems. In the future, we will design hybrid algorithms by taking the advantages of the above algorithms to handle more complex fuzzy JRD problems. Moreover, another future research direction is to optimize the fuzzy JRD problem under different operational risks because of the inevitable uncertainty under supply chain environment for enterprises in the network economic era [48ā50].
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research is partially supported by National Natural Science Foundation of China (nos. 71371080 and 71131004), Humanities and Social Sciences Foundation of Chinese Ministry of Education (no. 11YJC630275), and Fundamental Research Funds for the Central Universities (HUST: 2014QN201).
References
 S.L. Hsu, āOptimal joint replenishment decisions for a central factory with multiple satellite factories,ā Expert Systems with Applications, vol. 36, no. 2, pp. 2494ā2502, 2009. View at: Publisher Site  Google Scholar
 S. K. Goyal, āDetermination of optimum packaging frequency of items jointly replenished,ā Management Science, vol. 21, no. 4, pp. 436ā443, 1974. View at: Publisher Site  Google Scholar
 L. Wang, J. He, and Y. Zeng, āA differential evolution algorithm for joint replenishment problem using direct grouping and its application,ā Expert Systems, vol. 29, no. 5, pp. 429ā441, 2012. View at: Publisher Site  Google Scholar
 L. Wang, J. He, D. Wu, and Y. Zeng, āA novel differential evolution algorithm for joint replenishment problem under interdependence and its application,ā International Journal of Production Economics, vol. 135, no. 1, pp. 190ā198, 2012. View at: Publisher Site  Google Scholar
 M. Kaspi and M. J. Rosenblatt, āOn the economic ordering quantity for jointly replenished items,ā International Journal of Production Research, vol. 29, no. 1, pp. 107ā114, 1991. View at: Publisher Site  Google Scholar
 A. L. Olsen, āAn evolutionary algorithm to solve the joint replenishment problem using direct grouping,ā Computers and Industrial Engineering, vol. 48, no. 2, pp. 223ā235, 2005. View at: Publisher Site  Google Scholar
 F. C. Lee and M. J. Yao, āA global optimum search algorithm for the joint replenishment problem under poweroftwo policy,ā Computers and Operations Research, vol. 30, no. 9, pp. 1319ā1333, 2003. View at: Publisher Site  Google Scholar
 E. Porras and R. Dekker, āA solution method for the joint replenishment problem with correction factor,ā International Journal of Production Economics, vol. 113, no. 2, pp. 834ā851, 2008. View at: Publisher Site  Google Scholar
 M. Khouja and S. Goyal, āA review of the joint replenishment problem literature: 1989–2005,ā European Journal of Operational Research, vol. 186, no. 1, pp. 1ā16, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 A. Narayanan, P. Robinson, and F. Sahin, āCoordinated deterministic dynamic demand lotsizing problem: a review of models and algorithms,ā Omega, vol. 37, no. 1, pp. 3ā15, 2009. View at: Publisher Site  Google Scholar
 S. Sindhuchao, H. E. Romeijn, E. Akçali, and R. Boondiskulchok, āAn integrated inventoryrouting system for multiitem joint replenishment with limited vehicle capacity,ā Journal of Global Optimization, vol. 32, no. 1, pp. 93ā118, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 C. K. Chan, L. Y. Li, C. T. Ng, B. K. Cheung, and A. Langevin, āScheduling of multibuyer joint replenishments,ā International Journal of Production Economics, vol. 102, no. 1, pp. 132ā142, 2006. View at: Publisher Site  Google Scholar
 B. C. Cha, I. K. Moon, and J. H. Park, āThe joint replenishment and delivery scheduling of the onewarehouse, nretailer system,ā Transportation Research E: Logistics and Transportation Review, vol. 44, no. 5, pp. 720ā730, 2008. View at: Publisher Site  Google Scholar
 I. K. Moon, B. C. Cha, and C. U. Lee, āThe joint replenishment and freight consolidation of a warehouse in a supply chain,ā International Journal of Production Economics, vol. 133, no. 1, pp. 344ā350, 2011. View at: Publisher Site  Google Scholar
 L. Wang, C.X. Dun, W.J. Bi, and Y.R. Zeng, āAn effective and efficient differential evolution algorithm for the integrated stochastic joint replenishment and delivery model,ā KnowledgeBased Systems, vol. 36, pp. 104ā114, 2012. View at: Publisher Site  Google Scholar
 W. W. Qu, J. H. Bookbinder, and P. Iyogun, āAn integrated inventorytransportation system with modified periodic policy for multiple products,ā European Journal of Operational Research, vol. 115, no. 2, pp. 254ā269, 1999. View at: Publisher Site  Google Scholar
 L. Wang, H. Qu, S. Liu, and C. X. Dun, āModeling and optimization of the multiobjective stochastic joint replenishment and delivery problem under supply chain environment,ā The Scientific World Journal, vol. 2013, Article ID 916057, 11 pages, 2013. View at: Publisher Site  Google Scholar
 L. G. Cui, L. Wang, and J. Deng, āRFID technology investment evaluation model for the stochastic joint replenishment and delivery problem,ā Expert Systems with Applications, vol. 41, no. 4, pp. 1792ā1805, 2014. View at: Google Scholar
 M. S. Pishvaee and S. A. Torabi, āA possibilistic programming approach for closedloop supply chain network design under uncertainty,ā Fuzzy Sets and Systems, vol. 161, no. 20, pp. 2668ā2683, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Zeng, L. Wang, and J. He, āA novel approach for evaluating control criticality of spare parts using fuzzy comprehensive evaluation and GRA,ā International Journal of Fuzzy Systems, vol. 14, no. 3, pp. 392ā401, 2012. View at: Google Scholar
 J. Chiang, J. S. Yao, and H. M. Lee, āFuzzy inventory with backorder defuzzification by signed distance method,ā Journal of Information Science and Engineering, vol. 21, no. 4, pp. 673ā694, 2005. View at: Google Scholar
 J. Kacprzyk and P. Stanieski, āLongterm inventory policymaking through fuzzy decisionmaking models,ā Fuzzy Sets and Systems, vol. 8, no. 2, pp. 117ā132, 1982. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 L. Wang, H. Qu, Y. H. Li, and J. He, āModeling and optimization of stochastic joint replenishment and delivery scheduling problem with uncertain costs,ā Discrete Dynamics in Nature and Society, vol. 2013, Article ID 657465, 12 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 L. Wang, C. Dun, C. Lee, Q. Fu, and Y. Zeng, āModel and algorithm for fuzzy joint replenishment and delivery scheduling without explicit membership function,ā International Journal of Advanced Manufacturing Technology, vol. 66, no. 9–12, pp. 1907ā1920, 2013. View at: Publisher Site  Google Scholar
 H. Wee, C. Lo, and P. Hsu, āA multiobjective joint replenishment inventory model of deteriorated items in a fuzzy environment,ā European Journal of Operational Research, vol. 197, no. 2, pp. 620ā631, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Roy, S. Kar, and M. Maiti, āA deteriorating multiitem inventory model with fuzzy costs and resources based on two different defuzzification techniques,ā Applied Mathematical Modelling, vol. 32, no. 2, pp. 208ā223, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. Handfield, D. Warsing, and X. M. Wu, ā(Q, r) inventory policies in a fuzzy uncertain supply chain environment,ā European Journal of Operational Research, vol. 197, no. 2, pp. 609ā619, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 B. D. Liu, āDependentchance programming in fuzzy environments,ā Fuzzy Sets and Systems, vol. 109, no. 1, pp. 97ā106, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 L. Wang, Q.L. Fu, C.G. Lee, and Y.R. Zeng, āModel and algorithm of fuzzy joint replenishment problem under credibility measure on fuzzy goal,ā KnowledgeBased Systems, vol. 39, pp. 57ā66, 2013. View at: Publisher Site  Google Scholar
 H. Peng and H. Zhou, āA fuzzydependent chance multiobjective programming for water resources planning of a coastal city under fuzzy environment,ā Water and Environment Journal, vol. 25, no. 1, pp. 40ā54, 2011. View at: Publisher Site  Google Scholar
 L. Lu and M. E. Posner, āApproximation procedures for the onewarehouse multiretailer system,ā Management Science, vol. 40, no. 10, pp. 1305ā1316, 1994. View at: Publisher Site  Google Scholar
 R. Wang and L. Wang, āMaximum cut in fuzzy nature: models and algorithms,ā Journal of Computational and Applied Mathematics, vol. 234, no. 1, pp. 240ā252, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 X. B. Wang, W. S. Tang, and R. Q. Zhao, āFuzzy economic order quantity inventory models without backordering,ā Tsinghua Science and Technology, vol. 12, no. 1, pp. 91ā96, 2007. View at: Publisher Site  Google Scholar
 R. Narasimhan, āGoal programming in a fuzzy environment,ā Decision Science, vol. 11, no. 2, pp. 325ā336, 1980. View at: Google Scholar
 L. H. Chen and F. C. Tsai, āFuzzy goal programming with different importance and priorities,ā European Journal of Operational Research, vol. 133, no. 3, pp. 548ā556, 2001. View at: Publisher Site  Google Scholar
 H. Ke and B. Liu, āFuzzy project scheduling problem and its hybrid intelligent algorithm,ā Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 34, no. 2, pp. 301ā308, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 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  MathSciNet
 L. Wang, Q. Fu, and Y. Zeng, āContinuous review inventory models with a mixture of backorders and lost sales under fuzzy demand and different decision situations,ā Expert Systems with Applications, vol. 39, no. 4, pp. 4181ā4189, 2012. View at: Publisher Site  Google Scholar
 H. Qu, L. Wang, and Y. R. Zeng, āModeling and optimization for the joint replenishment and delivery problem with heterogeneous items,ā KnowledgeBased Systems, vol. 54, pp. 207ā215, 2013. View at: Google Scholar
 L. Wang, H. Qu, T. Chen, and F. P. Yan, āAn effective hybrid selfadapting differential evolution algorithm for the joint replenishment and locationinventory problem in a threelevel supply chain,ā The Scientific World Journal, vol. 2013, Article ID 270249, 11 pages, 2013. View at: Publisher Site  Google Scholar
 F. AlObeidat, N. Belacel, J. A. Carretero, and P. Mahanti, āDifferential Evolution for learning the classification method PROAFTN,ā KnowledgeBased Systems, vol. 23, no. 5, pp. 418ā426, 2010. View at: Publisher Site  Google Scholar
 L. Wang, H. Qu, S. Liu, and C. Chen, āOptimizing the joint replenishment and channel coordination problem under supply chain environment using a simple and effective differential evolution algorithm,ā Discrete Dynamics in Nature and Society, vol. 2014, Article ID 709856, 12 pages, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 Q. K. Pan and L. Wang, āA novel differential evolution algorithm for noidle permutation flowshop scheduling problems,ā European Journal of Industrial Engineering, vol. 2, no. 3, pp. 279ā297, 2008. View at: Publisher Site  Google Scholar
 D. Zou, H. Liu, L. Gao, and S. Li, āAn improved differential evolution algorithm for the task assignment problem,ā Engineering Applications of Artificial Intelligence, vol. 24, no. 4, pp. 616ā624, 2011. View at: Publisher Site  Google Scholar
 Y. M. Gong, Global Operations Strategy: Fundamentals and Practice, Springer, 2013.
 Y. H. Li, H. Guo, L. Wang, and J. Fu, āA hybrid geneticsimulated annealing algorithm for the locationinventoryrouting problem considering returns under esupply chain environment,ā The Scientific World Journal, vol. 2013, Article ID 125893, 10 pages, 2013. View at: Publisher Site  Google Scholar
 L. Cui, L. Wang, J. Deng, and J. Zhang, āA new improved quantum evolution algorithm with local search procedure for capacitated vehicle routing problem,ā Mathematical Problems in Engineering, vol. 2013, Article ID 159495, 17 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 S. Liu and L. Wang, āUnderstanding the impact of risks on performance in internal and outsourced information technology projects: the role of strategic importance,ā International Journal of Project Management, 2014. View at: Publisher Site  Google Scholar
 H. Soleimani, M. SeyyedEsfahani, and G. Kannan, āIncorporating risk measures in closedloop supply chain network design,ā International Journal of Production Research, vol. 52, no. 6, pp. 1843ā1867, 2014. View at: Google Scholar
 S. Liu and L. Wang, āUser liaisons perspective on behavior and outcome control in IT projects: role of IT experience, behavior observability, and outcome measurability,ā Management Decision, vol. 52, no. 6, pp. 1148ā1173, 2014. View at: Google Scholar
Copyright
Copyright © 2014 YuRong Zeng 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.