[Retracted] Combined Location-Inventory Optimization of Deteriorating Products Supply Chain Based on CQMIP under Stochastic Environment

Wang, Yong Ming; Yin, Hong Li

doi:https://doi.org/10.1155/2018/8210254

Mathematical Problems in Engineering

On this page

Abstract Introduction Literature Review Conclusions Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article Retraction

!

This article has been Retracted. To view the article details, please click the ‘Retraction’ tab above.

Research Article | Open Access

Volume 2018 | Article ID 8210254 | https://doi.org/10.1155/2018/8210254

[Retracted] Combined Location-Inventory Optimization of Deteriorating Products Supply Chain Based on CQMIP under Stochastic Environment

Yong Ming Wang¹and Hong Li Yin²

Academic Editor: Mauro Gaggero

Received02 Sept 2017

Revised13 Dec 2017

Accepted03 Jan 2018

Published05 Feb 2018

Abstract

The design and optimization of combined location-inventory model for deteriorating products are a main focus in supply chain management. There were many combined location-inventory design models in this field, but these models are under the assumptions of adequate capacity facilities, invariable lead time, unique product, and uncorrelated retailer’s demands. These assumptions have a big gap in the practical situation. In this paper, we design a combined location-inventory model for deteriorating products under capacitated facilities, stochastic lead time, multiple products, and correlated retailers’ stochastic demands assumptions. These constraints are near to actual supply chain circumstance. The problem is modeled as conic quadratic mix-integer programming (CQMIP) to minimize the total expected cost. We explain how to formulate these problems as conic quadratic mixed-integer problems, and in order to obtain better computational results we use extended cover cuts. Simultaneously we compare our method with the previous Lagrange methods; the result is that the new CQMIP method can get better solution.

1. Introduction

Design and optimization of the combined location-inventory for deteriorating products are extremely important. Deterioration is defined as “the damage, spoilage, dryness, etc. that results in a decrease of usefulness of the original one” [1]. Deteriorating products are classified as two kinds. The products that have a maximum usable lifetime are seen as “perishable products,” like foodstuffs, green vegetables, human blood, and photographic film. And the products that have no shelf-life at all are seen as “decaying products,” like alcohol, gasoline, and radioactive substances [2].

Deteriorating products about a third of sales at retailers, combined location-inventory design and optimization of the deteriorating products become a attention point. These problems are more difficult than combined location-inventory design and optimization of products with an infinite lifetime. We will consider minimizing the total expected cost under satisfying customer demand. In many cases, the holding cost and transportation cost per unit or per unit time for deteriorating products are not fixed. The cost of each unit deteriorating product is alterable; it depends on time, inventory level, storage location, and so forth. Therefore, the combined location-inventory system’s costs of deteriorating products are more complex. The influence of deteriorating products on the environment is a real and big problem; more and more academic researchers are paying attention to these deteriorating products supply chain.

There are numerous challenges and opportunities for deteriorating products supply chain. However, few rigorous and general supply chain models have been developed. Therefore, it is worthwhile and necessary to research deteriorating products supply chain. There are impacts on the economy and positive effects on health. At present, few works consider combined location-inventory in deteriorating products supply chain design and optimization. The majority of studies on deteriorating products are limited to inventory level or location level, rarely involved in combined location-inventory aspect. Hence, this paper considers a combined deteriorating supply chain design and optimization problem; we design a combined location-inventory model for deteriorating products. In our model, in order to make the model more realistic we consider capacitated facilities, stochastic lead time, multiple products, and correlated retailer’s demands.

In this paper, we propose a novel method based on conic quadratic mixed-integer programming, which is more flexible and global. We use optimization software CPLEX to solve this problem. The experimental results show that, compared to Lagrangian relaxation, however, our method has many advantages. And this research work provides an effective reference for managers of deteriorating products supply chain.

The remainder of this paper is organized as follows. In Section 2, we review the main researches in combined location-inventory models of supply chain. The notation and modeling assumptions are described in Section 3. In Section 4, a CQMIP model is presented for the combined location-inventory of deteriorating products. Numerical results together with related analyzes are presented in Section 5. Section 6 concludes this paper and briefly discusses future research directions.

2. Literature Review

Supply chain design and optimization are critical and difficult. It involves not only facility location but also facility inventory. Most of existing supply chain design and optimization models separately consider location decision and inventory decision, with no regard for the interaction between decision of location and decision of inventory. In order to reach efficient cost savings, the major cost components of location and inventory should be considered in combination, rather than separately. Thus, the idea of considering inventory costs while making facility location decisions of a supply chain has been developed during the past decade by different researchers. These models concurrently consider decisions at both the strategic location decisions level and the tactical inventory decisions level. Snyder et al. [3] proposed a stochastic location-inventory model with uncertainty. Shen and Qi [4] proposed a combined location-inventory routing decisions model. Ozsen et al. [5], directed at limited capacities for facilities, proposed an integrated location-inventory model. In these works, because of inventory consideration, the objective function is nonlinear. These models were complex, usually with integer decision variables, nonlinear costs, and constraints. So they proposed many heuristic algorithms for different special cases. Shen et al. [6] proposed a column generation method for this problem. Daskin et al. [7] and Shen and Daskin [8] all used Lagrangian relaxation method to solve this class of problems. Miranda and Garrido [9] considered facility location and inventory control decisions and presented an integrated model; the solution method involves Lagrangian relaxation and subgradient method. Erlebacher and Meller [10] proposed a two-stage heuristic algorithm to estimate the number of DCs. Dasci and Verter [11] introduced concave technology selection cost into the objective function of a multiproduct location model. Asl-Najafi et al. [12] considered a dynamic closed-loop location-inventory problem and designed a hybrid metaheuristic algorithm based on multiobjective particle swarm optimization (MOPSO) and nondominated sorting genetic algorithm-II (NSGA-II) to solve the problem. Zhang and Unnikrishnan [13] presented a location-inventory model with uncertain demands, which is based on integer nonlinear programming and can be transformed to conic quadratic mixed-integer programming, and used CPLEX software solve this problem. Diabat et al. [14] presented a joint location-inventory model based on uncertain demands and lead times, the use of which can determine not only the location and number of distribution centers but also the size, and adopted a hybrid algorithm to solve the presented model; this algorithm arose from simulated annealing and direct search. Vahdani et al. [15] considered a model based on mixed-integer nonlinear programming for a location-inventory problem, which assumed that the demands of retailers are correlated, and presented two metaheuristic algorithms, including genetic algorithm and simulated annealing to solve this problem.

Because the feature of the deteriorating products, there is a strict time restrict about deteriorating products. The product will lose its value rapidly, which can lead to extra waste and costs and loss potential market share. With increasing time-sensitivity for deteriorating products, timely delivery is becoming a strategy. A great number of tactics for time competition have been proposed; several methods were utilized for time-dependent supply chain design and optimization. Rich and Hines [16] proposed mathematical modeling, which focuses on time reduction for a single specific process. Guide Jr. et al. [17] discussed profit maximization in reverse supply chains for deteriorating return. Blackburn and Scudder [18] developed a cost minimization model for one specific deteriorating product supply chain, concerning the declining value of the product over time. Wang et al. [19] discussed inventory policies for a deteriorating product but based on a simple supply chain system consisting of one producer, one distributor, and only one retailer. Arntzen et al. [20] proposed mixed-integer linear programming model to minimize the weighted combination of total cost and activity days, but the unit processing activity days are fixed, with no regard for the facility capacities.

In the existing literature on deteriorating products, for describing either the decrease in quantity or the degeneration in quality, exponential time decay has been utilized. Regarding the decrease in quantity, Nahmias [21] and Yu and Nagurney [22] discussed perishable inventory optimization. Tijskens and Polderdijk [23] and Akkerman et al. [24] discussed the degeneration in quality but assumed that all the products deteriorate at the same rate.

3. Model Assumptions and Notations

The structure of the integrated supply chain is as in Figure 1. It is a three-level supply chain under the single-sourcing assumption; it has a central manufacturer, a set of distribution centers (DCs) that store and distribute supplies, and the retailers.

We used parameters and notations for the CQMIP model as follows:

Sets : set of deteriorating products, indexed by : set of retailers, indexed by : set of potential distribution center locations, indexed by

Demands : the daily demand at retailer for deteriorating product : standard deviation of daily demand at retailer for deteriorating product : variance-covariance matrix of retailers daily demand about deteriorating product

Costs : fixed cost of DC located at retailer site per unit time : cost per unit to transport deteriorating product between site and : fixed cost of placing an order deteriorating product at DC : fixed cost of placing an order deteriorating product at central manufacturer : per unit cost to transport deteriorating product from the central manufacturer to DC : fixed cost of per unit transport deteriorating product from the central manufacturer to DC : unit inventory holding cost per unit of deteriorating product per unit time : fixed disposing cost of deteriorating product after it becomes putrid

Weights : weight factor associated with the transportation cost for the deteriorating product : weight factor associated with the inventory cost for the deteriorating product

Other Parameters : days worked per unit time for the deteriorating product : service level of deteriorating product : standard normal deviation associated with : lead time in days for deteriorating product at DC : standard deviation of lead time in days for deteriorating product at DC : time of transportation of deteriorating product from the central manufacturer to DC : time of transportation of deteriorating product from site to site : maximum usable lifetime of deteriorating product : maximum inventory capacity of DC : quantity of deteriorating product at DC : decay rate of deteriorating product

Decision Variables : if demand for deteriorating product of retailer is assigned to DC at retailer site : if demand for deteriorating product of retailer is not assigned to DC at retailer site : if DC is located at retailer site : if DC is not located at retailer site : the reorder number of deteriorating product at DC .

We also consider the following assumptions:(1)Each retailer is supplied from unique DC.(2)Each DC has nonlinear capacity constraints.(3)Deteriorating products are transported directly from DCs to retailers.(4)Objective function is nonlinear.(5)Although the demand of retailer for deteriorating product is stochastic, there is a correlation between different retailers’ demands.(6)Lead times of deteriorating products are independent and orders not crossed.(7)The lead time of deteriorating product is stochastic.(8)Every deteriorating product has a lifetime; product value decreases over time.(9)The environmental conditions associated with each deteriorating product are fixed, such as temperature and transportation.

4. CQMIP Model

At least one of the variables is unknown in our model; we assume they follow some probability distribution. Stochastic supply chain optimization and design problems were usually described as integer or mixed-integer nonlinear optimization problems. In our paper, we will propose a combined location-inventory model based on conic quadratic mixed-integer programming. In this section, the conic quadratic mixed-integer programming (CQMIP) optimization problem is described first. Then, we use conic quadratic mixed-integer programming model of complex deteriorating product supply chain optimization and design problems.

4.1. CQMIP Optimization Problem

Now we refer to Ben-Tal and Nemirovski [25] and Alizadeh and Goldfarb [26] who define conic quadratic mixed-integer programming (CQMIP) optimization problem:where , is the Euclidean norm, and all the parameters are rational numbers. Let , ; we can obtain convex quadratic constraint.

4.2. CQMIP Model

In this section, we are going to consider the multiple deteriorating products combined location-inventory model with inventory capacity constraint, correlated retailers’ stochastic demands, stochastic lead times, and maximum usable lifetime.

Based on the above assumptions, the multiple deteriorating products combined location-inventory model is indicated in problem as follows:where is the cost of transportation product to retail site . , , is the assignment decision matrix for the th DC for the deteriorating product . is variance-covariance matrix of daily demand at retailers to deteriorating product . And is a survival probability of deteriorating product . The quantity surviving of deteriorating product at the end of the time interval follows a binomial distribution with parameters and . Hence, express the expected surviving quantity. The objective of the model is to minimize total expected costs of location, inventory, transportation, and loss of deteriorating product value. The first objective term is fixed cost of DC located at retailer site per unit time. The second objective term is the sum of transportation costs, inventory costs, and the loss of deteriorating product value. Constraint (2) defines the capacity of each DC. Constraint (3) ensures that each deteriorating product can be delivered in usable lifetime period. Constraint (4) ensures that each retailer with demand for deteriorating product is supplied only by one DC. Constraint (5) illustrates that DC can provide service for retailers only when potential DC becomes a real DC. Constraint (6) defines the range of the reorder number of deteriorating product at DC .

By introducing auxiliary variables , for each and , is defined by constraint (14), and is defined by constraint (15). Now we, based on the above formulation, give conic quadratic mixed-integer programming model (CQMIP):where .

Variable is used to substitute the term , and we can prove that problem is equivalent to CQMIP.

4.3. Extended Cover Cuts

We use optimization software CPLEX to solve CQMIP problem. In order to get a better solution we strengthen the formulations (CQMIP) with cover inequalities, utilize the left side of inequality (14) to substitute , and then we receive the nonlinear 0-1 constraint knapsack result:

For inequality (17), define a set function that is from set to set :

Using submodularity of set function, we can gain the extended cover cuts for the submodular knapsack set :

Because the validity of the cover inequality depends on not only the submodular knapsack set but also the cover set, define a cover set as follow: , which is called a cover set for if .

Narayanan [27] shows that for cover set , cover inequality is valid for . We will utilize this cover inequality in our computations.

5. Numerical Experiments and Sensitivity Analysis

In this section, we make numerical analysis to computational results about CQMIP formulations of the combined location-inventory problems in the deteriorating products supply chain and compare our computational results with the earlier Lagrangian relaxation method.

There are three data sets in our experiment: 18-node, 80-node, and 150-node date set. In Table 1, we give the results obtained by running CQMIP, compared with Lagrangian relaxation method. The results indicate that our method performs quite well in this experiment, reaching optimum solution faster than Lagrangian relaxation method.

In our method, we add extended cover inequalities; in Table 2, we give the results obtained by adding extended cover inequalities. We can find that, with capacity decrease, the problem becomes more difficult, but, adding extended cover inequalities, we can find that the solution time is reduced and the number of nodes is also reduced.

Now we examine the influence of correlated retailers’ stochastic demands on the combined location-inventory problems in the deteriorating products supply chain. We utilize 18-node data set to illustrate the effect of retailers’ stochastic demand correlation. In Figure 2, the retailers’ demands are uncorrelated, and the links show retailer assignments in the optimal solution. Four DCs are opened in Beijing, Shanghai, Guangzhou, and Chengdu, and the total expected cost is 180090. Now we add correlation to the demands of the retailers, to see how correlation will change the solution. We assume the correlation between retailers is up to 80%; in this case the total expected cost is 186595. The current solution is no longer optimal, but the inventory level and subsequent costs are lower.

We consider positive correlation between retailers; Figure 3 shows the total expected cost is a function of the retailers’ demands correlation; the total expected cost increases monotonically. Moreover, Figure 3 shows the number of opened DCs is a function of the retailers’ demands correlation; from this we can conclude that the retailers’ positive demands correlation reduces the risk of location.

In our model we consider lead time is stochastic; now let the standard deviation of the lead time be 0; Figure 4 shows this particular case. Based on Figure 4, in order to explain lead time variability, we add a third axis. The influence of retailers’ correlation demands and stochastic lead times on total cost is captured by Figure 4; we can find that when retailers’ stochastic demands are correlated and/or stochastic lead time standard deviation increases, the total cost of supply chain will be increased.

Now let the number of opened DCs increase to 18, but keep lead time standard deviation invariability; from Figure 5 we can conclude that with the retailers’ stochastic demands correlation increase, there is an increase in the number of opened DCs. We also conclude that the number of opened DCs is a function of retailers’ stochastic demands correlation and lead time standard deviation.

Now we examine the influence of and on the number of opened distribution centers. We utilize 150-node data set to illustrate the influence; in our experiment the values of are 0.001, 0.003, 0.012, 0.060, and 0.360, and the values of are 0.100, 0.300, 0.900, 2.700, and 4.100. From Figures 6 and 7 we can see that when the value of increases, the number of opened distribution centers increased, and when the value of increases, the number of opened distribution centers also increased. Because increases equally to the corresponding transportation cost increase, increases equally to the corresponding inventory cost increase, so the fixed-location cost has fewer impacts on the cost of the objective function; this results in more distribution centers being opened.

From Figures 8 and 9 we can see that when increases, the sum of total expected cost decreases, and when increases, the sum of total expected cost increases. Figures 8 and 9 show the total expected cost as a function of and . Taking into account inventory decision and location decision at the same time, it will bring a larger total expected cost reduction.

6. Conclusions

In this paper, we propose a new method to model and optimize combined location-inventory deteriorating products supply chain problem. We think capacity of facilities are limited, lead time is random, the products are varied, and retailers’ random demands are correlated. These assumptions are more practical and near to actual supply chain circumstance. We use conic quadratic mixed-integer programming to model integrated location and inventory problem of deteriorating products supply chain, and we used optimization software CPLEX to solve this problem. Not only does this new method obtain better computational solution than previous Lagrange methods, but also our model is useful to more general problems. Our method can be used for inventory management problems, location management problems, and even transportation management problems. Conic quadratic mixed-integer programming also can be used for service supply chain.

The research work provides an effective reference for managers of deteriorating products supply chain. Considering the correlation of the demand of the retailers, it can effectively reduce the inventory and reduce the cost of excessive inventory. It can better control the location of the distribution center and reduce the risks caused by improper location selection. With stochastic lead time standard deviation increase, the total cost of supply chain will increase. So supply chain managers need to control the stochastic lead time standard deviation in the right range. The number of opened DCs is a function of retailers’ stochastic demands correlation and lead time standard deviation, so supply chain managers can effectively control the number of opened DCs in virtue of retailers’ stochastic demands correlation and lead time standard deviation. The above measures make the supply chain achieve the minimum expected total cost.

In our model, we only consider combined location-inventory condition. In subsequent studies, we will continue studying combined location-inventory model with vehicle routing, applying conic integer programming approach to model combined production and transportation case.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work is supported by Natural Science Foundation of China (no. 71640028, no. 71362030, no. 71262029) and Natural Science Foundation of Yunnan (no. 2013FZ048).

References

K.-S. Wu, L.-Y. Ouyang, and C.-T. Yang, “An optimal replenishment policy for non-instantaneous deteriorating items with stock-dependent demand and partial backlogging,” International Journal of Production Economics, vol. 101, no. 2, pp. 369–384, 2006.
View at: Publisher Site | Google Scholar
S. K. Goyal and B. C. Giri, “Recent trends in modeling of deteriorating inventory,” European Journal of Operational Research, vol. 134, no. 1, pp. 1–16, 2001.
View at: Publisher Site | Google Scholar | MathSciNet
L. V. Snyder, M. S. Daskin, and C.-P. Teo, “The stochastic location model with risk pooling,” European Journal of Operational Research, vol. 179, no. 3, pp. 1221–1238, 2007.
View at: Publisher Site | Google Scholar
Z. J. Shen and L. Qi, “Incorporating inventory and routing costs in strategic location models,” European Journal of Operational Research, vol. 179, no. 2, pp. 372–389, 2007.
View at: Publisher Site | Google Scholar
L. Ozsen, C. R. Coullard, and M. S. Daskin, “Capacitated warehouse location model with risk pooling,” Naval Research Logistics (NRL), vol. 55, no. 4, pp. 295–312, 2008.
View at: Publisher Site | Google Scholar | MathSciNet
Z.-J. M. Shen, C. R. Coullard, and M. S. Daskin, “A joint location-inventory model,” Transportation Science, vol. 37, no. 1, pp. 40–55, 2003.
View at: Publisher Site | Google Scholar
M. S. Daskin, C. R. Coullard, and Z. . Shen, “An inventory-location model: formulation, solution algorithm and computational results,” Annals of Operations Research, vol. 110, pp. 83–106, 2002.
View at: Publisher Site | Google Scholar | MathSciNet
Z.-J. M. Shen and M. S. Daskin, “Trade-offs between customer service and cost in integrated supply chain design,” Manufacturing and Service Operations Management, vol. 7, no. 3, pp. 188–207, 2005.
View at: Publisher Site | Google Scholar
P. A. Miranda and R. A. Garrido, “Incorporating inventory control decisions into a strategic distribution network design model with stochastic demand,” Transportation Research Part E: Logistics and Transportation Review, vol. 40, no. 3, pp. 183–207, 2004.
View at: Publisher Site | Google Scholar
S. J. Erlebacher and R. D. Meller, “The interaction of location and inventory in designing distribution systems,” IIE Transactions, vol. 32, no. 2, pp. 155–166, 2000.
View at: Google Scholar
A. Dasci and V. Verter, “A continuous model for production-distribution system design,” European Journal of Operational Research, vol. 129, no. 2, pp. 287–298, 2001.
View at: Publisher Site | Google Scholar | MathSciNet
J. Asl-Najafi, B. Zahiri, A. Bozorgi-Amiri, and A. Taheri-Moghaddam, “A dynamic closed-loop location-inventory problem under disruption risk,” Computers & Industrial Engineering, vol. 90, pp. 414–428, 2015.
View at: Publisher Site | Google Scholar
Z.-H. Zhang and A. Unnikrishnan, “A coordinated location-inventory problem in closed-loop supply chain,” Transportation Research Part B: Methodological, vol. 89, pp. 127–148, 2016.
View at: Publisher Site | Google Scholar
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
B. Vahdani, M. Soltani, M. Yazdani, and S. Meysam Mousavi, “A three level joint location-inventory problem with correlated demand, shortages and periodic review system: Robust meta-heuristics,” Computers & Industrial Engineering, vol. 109, no. 7, pp. 113–129, 2017.
View at: Publisher Site | Google Scholar
N. Rich and P. Hines, “Supply-chain management and time-based competition: The role of the supplier association,” International Journal of Physical Distribution & Logistics Management, vol. 27, no. 3-4, pp. 210–225, 1997.
View at: Publisher Site | Google Scholar
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
J. Blackburn and G. Scudder, “Supply chain strategies for perishable products: the case of fresh produce,” Production Engineering Research and Development, vol. 18, no. 2, pp. 129–137, 2009.
View at: Publisher Site | Google Scholar
K.-J. Wang, Y. S. Lin, and J. C. P. Yu, “Optimizing inventory policy for products with time-sensitive deteriorating rates in a multi-echelon supply chain,” International Journal of Production Economics, vol. 130, no. 1, pp. 66–76, 2011.
View at: Publisher Site | Google Scholar
B. C. Arntzen, G. G. Brown, T. P. Harrison, and L. L. Trafton, “Global Supply Chain Management at Digital Equipment Corporation,” Interfaces, vol. 25, no. 1, pp. 69–93, 1995.
View at: Publisher Site | Google Scholar
S. Nahmias, “Perishable inventory theory: a review,” Operations Research, vol. 30, no. 4, pp. 680–708, 1982.
View at: Publisher Site | Google Scholar
M. Yu and A. Nagurney, “Competitive food supply chain networks with application to fresh produce,” European Journal of Operational Research, vol. 224, no. 2, pp. 273–282, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
L. M. M. Tijskens and J. J. Polderdijk, “A generic model for keeping quality of vegetable produce during storage and distribution,” Agricultural Systems, vol. 51, no. 4, pp. 431–452, 1996.
View at: Publisher Site | Google Scholar
R. Akkerman, P. Farahani, and M. Grunow, “Quality, safety and sustainability in food distribution: A review of quantitative operations management approaches and challenges,” OR Spectrum, vol. 32, no. 4, pp. 863–904, 2010.
View at: Publisher Site | Google Scholar
A. Ben-Tal and A. Nemirovski, Lectures on Modern Convex Optimization: Analysis, Algorithms, and Engineering Applications, MOS-SIAM Series on Optimization, Society for Industrial and Applied Mathematics (SIAM), Pennsylvania, Pa, USA, 2001.
View at: Publisher Site | MathSciNet
F. Alizadeh and D. Goldfarb, “Second-order cone programming,” Mathematical Programming, vol. 95, no. 1, pp. 3–51, 2003.
View at: Publisher Site | Google Scholar
V. Narayanan, Lifting for Conic Mixed-Integer Programming, Springer-Verlag Inc, New York, NY, USA, 2011.

Copyright

Copyright © 2018 Yong Ming Wang and Hong Li Yin. 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.

PDF Download Citation

Download other formats

Order printed copies