CC BY
32
THE REVERSE LOGISTICS MODEL WITH REUSING OF COMPONENTS OF SERIES
SYSTEM PRODUCT
M. Plewa, A. Jodejko-Pietruczuk
Wroclaw University of Technology, Wroclaw, Poland e-mail: [email protected]
ABSTRACT
The main goal of this paper is to create the reverse logistics model that uses reliability theory to describe reusability of product parts with assumption that recovered components are used in production process but they aren't as good as new ones. The model allows to estimate the potential profits of the reusing policy in a production and gives the base to optimize some of the process parameters: the threshold work time of returns or the warranty period for products containing reused elements.
1 INTRODUCTION
Until recently logistics systems supported only processes carried out in classical material flow from producer to final user. Recently it has been a remarkable growth of interest in optimizing logistics processes that supports recapturing value from used goods. The process of planning, implementing, and controlling the efficient, cost effective flow of raw materials, in-process inventory, finished goods and related information from the point of consumption to the point of origin for the purpose of recapturing value or proper disposal is called reverse logistics. Reverse logistics has become one of the logicians' key areas of interest. It enjoys ever-increasing interest of many industrial branches. Nowadays a growing number of companies realize the meaning of that field of logistics. Reuse of products or product parts can bring direct advantages to the company because it reduces costs associated with acquiring new components by using recycled materials or recovered components instead of expensive raw material. Literature survey that has been done around the theme of the reverse logistics area, allowed to set out this article aims and objectives.
In the reverse supply chain the issue of: timings, quantities and conditions of returned and reusable components make production planning difficult (Murayama, T. et al. 2006). The majority of models assume that demand for new products and returns quantity are independent Poisson random variables (Plewa, M. & Jodejko-Pietruczuk, A. 2011) and very few use the reliability theory to estimate the number of reusable products. Murayama et al. (Murayama, T. & Shu, L. H. 2001, Murayama, T. et al. 2005, Murayama, T. et al. 2006) propose the method to predict the number of quantities of returned products and reusable components at each time period by using series system reliability models. The condition aspect of returns in a reverse logistic system is usually omitted by using assumption that all the returns are reusable and usually "as good as new" (Kiesmuller, G.P. 2003). Only few models use the reliability theory to diversify returned elements' reusability, but they don't give any guidelines for the way to optimize the threshold value of components' residual life (e.g. (Murayama, T. & Shu, L. H. 2001, Murayama, T. et al. 2005, Murayama, T. et al. 2006)).
In literature the field of reverse logistics is usually subdivided into three areas: inventory control, production, recovery and distribution planning. The model presented in this paper is a
inventory model and hence the literature survey refers to this area. Basic inventory model in reverse logistics is built on following assumptions (Fleischmann, M. et al. 1997):
- the manufacturer meets the demand for the final products;
- the manufacturer receives and stores the products returned from the final user;
- demand for the new products can be ful-filled by the production or by the recovery of returned products;
- recovered products are as good as new ones;
- the goal of the model is to minimize the total costs.
Inventory models in reverse logistics can be divided into deterministic and stochastic models.
Deterministic models presented in the literature are mainly the modifications of the EOQ model (e.g. (Dobos, I. & Richter, K. 2000, Dobos, I. 2002, Mabini, M.C. et al. 1992, Richter, K. 1994, Richter, K. 1997, Richter, K. 1996a, b, Schrady, D.A. 1967, Teunter, R.H. 2001)). There are also some models based on dynamic programming which are the extensions of the classical Wagner Whitin model (e.g. (Kleber, R. et al. 2002, Konstantaras, I. & Papachristos, S. 2007, Richter, K. & Sombrutzki, M. 2000, Richter, K. & Weber, J. 2001)).
In the class of stochastic models there are two model groups:
- models, in which demand for the new product is s a consequence of the return;
- models, in which the demand for new products does not depend on the number of returns, but the number of returns may depend on the previous demand.
In first group there are models of repair systems. These are closed systems in which the number of elements remains constant. The main goal of such a system is to keep sufficient number of spare parts to provide the required level of availability of the technical system (e.g. (Muckstadt, J.A. 1973, Sherbrooke, C.C. 1971, Sherbrooke, C.C. 1968)). The second group of stochastic inventory models can be divided into continuous (e.g. (Fleischmann, M. et al. 2002, Fleischmann, M. & Kuik, R. 2003, Heyman, D.P. 1997, Korugan, A. & Gupta, S. M. 1998, Muckstadt, J.A. & Isaac. M.H. 1981, Van der Laan, E.A. et al. 1996, Van der Laan, E.A. et al. 1996, Van der Laan, E.A. & Salomon, M. 1997]) and periodic review models (e.g. (Inderfurth, K, Inderfurth, K. & van der Laan, E. 2001, Kelle, P. & Silver, E.A. 1989, Simpson, V.P. 1978)). A more detailed description of inventory models in reverse logistics can be found in (Plewa, M. & Jodejko-Pietruczuk, A. 2011).
Literature review allows to summarize the current state of knowledge and to define the main shortages of existing logistics models that deal with the reverse logistics problem. The majority of models assume that demand for new products and returns quantity are independent Poisson random variables (Plewa, M. & Jodejko-Pietruczuk, A. 2011). Few authors examine the relationship between the demand, and the number of returns but there are no inventory models in reverse logistics that use reliability theory to assess the number of returns and very few use the reliability theory to estimate the number of reusable products. Murayama et al. (Murayama, T. & Shu, L. H. 2001, Murayama, T. et al. 2005, Murayama, T. et al. 2006) propose the method to predict the number of quantities of returned products and reusable components at each time period by using series system reliability models. The condition aspect of returns in a reverse logistic system is usually omitted by using assumption that all the returns are reusable and usually "as good as new" (Kiesmuller, G.P. 2003). Only few models use the reliability theory to diversify returned elements' reusability, but they don't give any guidelines for the way to optimize the threshold value of components' residual life (e.g. (Murayama, T. & Shu, L. H. 2001, Murayama, T. et al. 2005, Murayama, T. et al. 2006)). Most of created models assume single component product.
Main goal of this paper is to create the reverse logistics inventory model that uses the reliability theory to describe reusability of product parts with assumption that recovered components are used in a production process but they aren't as good as new ones. The model allows to estimate the potential profits of the reusing policy in production and inventory management. It
gives the base to optimize some of the process parameters: the threshold work time of returns, the warranty period for products containing reused elements.
2 THE MODEL OF THE REUSING POLICY
The model that is presented in the paper is based on the following assumptions:
- A company produces the object composed of two elements (A and B). The product fails when one of components fails - series reliability structure.
- A failure of each component occurs independently on other components' failures.
- If the product fails during the warranty period, it is returned to the manufacturer and he has to pay some penalty cost (e.g. the cost of a new product).
- The products are returned as soon as their lives are ended and reusable B components are stored in a stock until new production batch running, when they may be reused.
- The component B of the product may be reused in a new production, if it was not the cause of a product failure and its total work time up to this moment is not greater than some acceptable - threshold time T (Fig. 1).
- Neither failed elements B can be reused in a new production (not repairable) nor any A element. All A components are new in a new production.
- Demand for the products is determined and fixed.
- New products are manufactured and sold periodically in established moments.
The process of reusing of the component B, dependently on its threshold age T, is shown in (Fig. 1).
Figure 1. Process of element reusing in a new production.
Despite the fact that companies realize the potential of reusing products, the question "is it worth to do it", is not so simple to answer. Within main reasons for products reusing are: difficulties with raw material supplying, high cost of utilization of returned and damaged products or lower cost of reusing of products' components than buying new ones. The objective of the presented model is to estimate profitability of using returned and recovered elements in a new production, in the case when they are not as good as new.
According to the assumptions, before every production beginning, the manufacturer has to make the decision: which of returned elements B should be used in the new production. The usage of recovered components decreases production costs but also increases the risk that additional costs occur because of larger amount of returns during the warranty period.
The objective is to find the threshold work time T for returned element that equalizes potential cost and profits of the reusing policy:
e(cwa (tw,t)) e(cwn (tw , t))- cb cr cm
(1)
E(cwa (tw , t)) -
RB (Tw + T )Ra (Tw )
. RB(T) .
c
E(Cwn (Tw , T )) - [l _ Rb (Tw )Ra (Tw )lCa
(2)
(3)
where CWO = the cost of warranty services if an "old" element is used in a new production; TW = the warranty period of the product; T = the threshold age of the element B, after that the further exploitation isn't continued; CWN = the cost of warranty services if a "new" element is used in production; CB = the purchase cost of a new element B; CR = the total cost of all activities of: decomposition, cleaning, preparing of the returned B element to reusing in a production; Co = penalty cost resulting from a product failure during warranty period (e.g. the total cost of production of a new object ); RA(t) = reliability of the element A in t moment; RB(t) = reliability of the element B in t moment; n = production batch percent of reusable B elements, that return during warranty period; E(x) = the expected value of variable x.
The left side of the equation 1 specifies the increase in expected costs of product warranty services (during a single warranty period) caused by the reusing in a production element that is not as good as new. This part of the expression depends on: both elements' reliability (A and B), the length of warranty period and the length of the acceptable total work time T of the returned B element. The increase in expected cost of warranty services is calculated for the case, when all reused elements are in the age of T. The real age of reused elements is equal or lower than T and this way the left side of the expression 1 estimates the maximum possible growth in warranty cost when "old" elements are reused in the production. The direction of changes in the expected cost of reusing elements is not so obvious. In special cases (low value of time T and short warranty period) it can happen that reused product is more reliable than a new one.
The right side of the equation 1 determines the potential cost savings resulting from using cheaper, recovered components instead of new elements.
On the basis of this expression the threshold value T of total work time of the element B can be found, above which reusing the component is not economical. An analytical solution of the equations can't be achieved for the majority of probability distributions describing components' time to failure, but the value of T can be easily found by applying numerical calculations for given vector of T.
The practical application of the proposed model is limited because of the mentioned simplification that all reused elements are in the age of T. According to model assumptions, the demand for the products is determined and fixed. It means that a new production is usually mix of
new and reused elements, dependently on their accessibility. The threshold work time T (for which savings from components' reusing are equal losses) can be specified on the base of Equation 1, but practical questions require more precise data: how many reusable elements will return before the new production batch, how many new components must be kept in a stock or what is the expected profit/cost when mix of new-old elements is used in the production?
To answer this questions, the model should consider the percent of returns used in the production. The number of returns that can return between two moments of a production beginning and may be reused depends on the number of the products that were sold earlier, the length of the period between two consecutive production batch, the length of the warranty period and threshold work time T. The number of returns (calculated as a percent of the production batch size) may be estimated as follows (Plewa, M. & Jodejko-Pietruczuk, A. 2011):
«(ta , tb , TW >T) = ni + n2 + n3 (4)
J'^ 0)D(t)F(ta -1, TM )dt Jmax(ta-tm ) (5)
2 =f T D (t )F (ta -t' tb -t )dt (6)
n3 = [ibD(t)F(°, tb -t)dt ^a (7)
tm = min(tw,t) (8)
F (ti, t2 ) = F (t 2)-F (ti) (9)
F (t ) = (1- Ra (t ))RB (t) (1°)
F (t ) = (1-Ra (t ))RB (t + T-Tm ) (11)
where n = production batch percent of reusable B elements, that return during warranty period; TW = the warranty period of the product; T = the threshold age of the element B, after that the further exploitation isn't continued; TM = minimal value of threshold time T and warranty period TW; D(t) -size of the production sold in t moment; F(ti,t2) = the increase of product unreliability (caused by A component when B element is still working) in period between ti and t2 moments; RA(t) = reliability of the element A in t moment; RB(t) = reliability of the element B in t moment; min, max (ti,t2) = minimum/maximum value of variables ti and t2.
0 t
0
T
-l J
M
T
T
W
11
a i
i
reusable returns
1-1
global tb time
t -1 T
La 1 1M
local time
Figure 2. The period when reusable components may be returned to a producer according to
equation 5.
ni =
0
Figure 3. The period when reusable components may be returned to a producer according to
equation 6.
I-----f
0
t t
a
Figure 4. The period when reusable components may be returned to a producer according to
equation 7.
The expression 4 allows to calculate the percent of B elements that can be returned to a manufacturer between two consecutive production batch and may be reused for the case when t2 - t1 < TM. Value TM is the minimal value of warranty period and threshold time lengths. It limits the possibility of returns: component B may be returned only during the product warranty period and reused only if it is not older then T. The expression may estimate maximal or minimal number of reusable B elements dependently on the form of the reliability function of B component. If the formula 10 is used, the maximum number of reusable elements is estimated according to the assumption that all products were as good as new when they were sold. The minimal number of reusable returns may be calculated when the expression 11 is considered in the equation 4 because it assumes that all components B in new products are in the maximal allowed age of T - TM. Other
cases (when t2 - t1 > TM) are presented in detail in literature (Plewa, M. & Jodejko-Pietruczuk, A. 2011).
Real values of costs and savings coming from the reusing policy is proportional to the variable n:
E Co Tw , T)) - E(Cwn {Tw ))]n = C - CR )n - CM (12)
where CWO = the cost of warranty services if "old" element is used in a new production; Tw = the warranty period of the product; T = the threshold age of the element B, after that the further exploitation isn't continued; CWN = the cost of warranty services if "new" element is used in a production; n = production batch percent of reusable B elements, that return during warranty period; Cb = the purchase cost of a new element B; CR = the total cost of all activities of preparing of returned B element to reusing in a production; CO = penalty cost resulting from a product failure during warranty period; E(x) = the expected value of variable x.
The cost of elements' identification CM is independent on the number of returns, because if a company decides to apply the reusing policy, all elements used in a production have to be identified.
3 RESEARCH RESULTS
Some example of analytical model results are presented in figures 5-17 (for process parameters: CO = 5, CB = 1, CR = CM = 0, RA = RB = exp(-(t/100))2.
Figure 5 presents the number of reusable returns assuming that all new products sold earlier:
- were as good as new (maximal number),
- include reused B elements (minimal number)
Figure 5. Maximal and minimal production batch percent of reusable B elements, that may be used
in new production.
The minimal and maximal number of reusable elements are similar for low values of threshold time T and warranty period Tw. When an acceptable time T exceeds the average life time of the element B, the minimal number of reusable elements decreases very fast - second reusing of "old" elements is little probable.
Figure 6. Increase of unreliability costs and savings resulting from the reusing policy according to
equation 1.
Figure 7. Increase of unreliability costs and savings resulting from the reusing policy according to equation 12 with maximal possible number of reusable returns.
Figure 8. Increase of unreliability costs and savings resulting from the reusing policy according to equation 12 with minimal possible number of reusable returns.
The results in Figures 6-8 are obtained for the same process parameters, but are different because of the influence of returned elements amount, used in a new production batch. Savings
coming from element reusing in Figure 6 (without considering the number of returns) are constant for all analysed values of T and TW, but they have irregular shape in Figure 7,8. Considering two alternative values of the total work time T, sometimes it is more economical to get lower unit profit but coming from greater quantity of reused elements. Potential savings and costs are proportional to the number of reusable elements (Fig. 5). Their real values will take place between plates presented in Figure 9 and it is seen that incorrect estimation of T time (too long) causes great costs of the reusing policy. Irregular shape of reusing policy benefits and costs shows that the plane may have more than one local maximum or minimum. The research conducted for various process parameters haven't given any general rules where the optimal value of work time T may be found.
The estimation of threshold and optimal work time T for various warranty periods proves that they are the same for the minimal and maximal number of returns (Fig. 10,11).
maximal number of returns ♦
minima! number of returns + +
Figure 9. Summary profit resulting from reusing of elements B in new production.
Figure 10. Optimum value of T (savings resulting from the reusing policy get maximum).
Figure 11. Threshold value of T (reusing is not economical above this value).
The optimum work time T, after that the further exploitation isn't continued, rises when a warranty period gets longer, but does not exceed the average lifetime of the component B.
The shape of threshold values of T time may also be interesting (Fig. 11). For very short warranty periods costs are usually equal to savings and are very close to zero. It is the effect of very low changes of reliability of the product in very short time. In such cases reusing even "very old" (high values of threshold time T) component is profitable.
In order to assess the influence of the reusing process parameters on its cost results, the sensitivity analysis of the proposed model was conducted. The parameters that was concerned as the most meaningful were tested and the results obtained during the research are shown in Figures 12 -17.
When a returned component that was working some time in the past is used in the product offered as new for a customer, the element's intensity of failures is one of the determinant of possible reusing. For this reason the impact of the reusable component failure rate on production process costs and profit was tested.
Figure 12. Threshold value of T for various failure intensity of the component B: RB = exp(-
(t/100))D.
- <-
❖ ❖ O ❖
-a -
+
20 40 60 "BÛ 100120140160 Tw
Figure 13. Optimal value of T for various failure intensity of the component B: RB = exp(-(t/100))D.
Exponentially distributed time to failure has a constant failure rate and reusing of such elements always is profitable (Fig. 12,14). The optimal value of work time T, for which the profit is the highest, gets maximum possible value for all tested warranty periods (Fig. 13). The sudden decrease in optimal T value for warranty periods two times longer than co-component's (A) expected time to failure is the effect of lower number of possible returns of the product. Only in the case when A element fails during a warranty period and is returned to the manufacturer, reusing of B components is possible. The optimum T takes value of double A lifetime.
Figure 14. Summary profit resulting from reusing of elements B in new production for RB = exp(-(t/100))1.
summary profit fora = 1
Figure 15. Summary profit resulting from reusing of elements B in new production for RB = exp(-(t/100))2.
Figure 16. Summary profit resulting from reusing of elements B in new production for RB = exp(-(t/100))3.
Figures 15,16 present cost results of the reusing policy for the growing failure rate of the reusable component. The faster growth of this parameter shifts the threshold and the optimal T to lower values - reusing is less profitable (Fig. 12,13).
The second most meaningful parameter of the reusing process is the relationship between possible costs of warranty services and savings coming from the lower number of elements in the production process. The research results are presented in Figures 17-18.
Co/CB = 5 Co/Cb =10
i$a o 20 40 60 80^00120«0«0
Figure 17. Summary profit resulting from reusing of elements B for cO = 5 and cO = 10, cB = 1, cR =
cm = 0.
Figure 18. Threshold value of T for various relationship between possible costs ratios.
The results of the research are quite obvious - higher cost of a product failure during warranty period causes lower profitability of reusing policy.
4 CONCLUSIONS
The model presented in this paper is the continuation of wide researches in the reverse logistics area. The majority of models deal with single - element system or with the assumption that reused elements are as good as new. The proposed model develops the previous ones by releasing both assumptions and gives the base to determine some of reusing policy parameters such as: the threshold work time of returned element that can be used again, the warranty period for the product containing elements which have some history of work, the size of new elements' stock necessary to fulfil production planes. The model is presented and tested for two-element, series system but it is very simple to be developed it to the case of x-element, series system. From the point of view of
CO/CB = 5 ♦ CO/CB =10 +
O ^ :::::: a a
6 20 40 60 "B0 100120140160 Tw
possible product returns during a warranty period this assumption is usually real because only a failure of a one product component allows to reuse others.
Although the majority of presented expressions is difficult to solve analytically, their analytical form enables easy implementation to a numerical search of optimal process parameters.
Practical application of the model is possible in companies which are interested in the reusing policy because the decision must be supported by the calculation of possible cost and benefits. There are many possible extensions of the presented model by taking into account:
- parameters associated with the process of collection and transportation of returns. Particularly important may be the consideration of different return sources;
- longer than a period, time of recovery;
- random times of recovery;
- random lead time for external orders;
- ability to deliver returns in batches;
- non-deterministic demand for the final product;
- complex reliability structure of the technical objects.
5 REFERENCES
Dobos, I. & Richter, K. 2000. The integer EOQ repair and waste disposal model - Further analysis.
Central European Journal of Operations Research 8: 173-194. Dobos, I. 2002. The generalization of Schrady's model: a model with repair, Working Paper No. 7, Department of Business Economics, Budapest University of Economics and Public Administration.
Fleischmann, M. & Bloemhof-Ruwaard, J.M. & Dekker, R. & van der Laan, E. & van Nunen, J.A.E.E. & van Wassenhove, L.N. 1997. Quantitative models for reverse logistics: A review. European Journal of Operational Research 103: 1-17. Fleischmann, M. & Dekker, R. & Kuik, R. 2002. Controlling inventories with stochastic item
returns: A basic model. European Journal of Operational Research 138: 63-75. Fleischmann, M. & Kuik, R. 2003. On optimal inventory control with independent stochastic item
returns. European Journal of Operational Research 151: 25-37. Heyman, D.P. 1997. Optimal disposal policies for a single-item inventory system with return.,
Naval Research Logistics Quarterly 24: 385-405. Inderfurth, K., Simple optimal replenishment and disposal policies for a product recovery system
with leadtimes. OR Spektrum 19: 111-122. Inderfurth, K. & van der Laan, E. 2001. Leadtime effects and policy improvement for stochastic inventory control with remanufacturing., International Journal of Production Economics 71: 381390.
Kelle, P. & Silver, E.A. 1989. Purchasing Policy of New Containers Considering the Random
Returns of Previously Issued Container. IIE Transactions 21 (4): 349-354. Kiesmuller, G.P. 2003. A new approach for controlling a hybrid stochastic manufacturing/remanufacturing system with inventories and different leadtimes. European Journal of Operational Research 147: 62-71. Kleber, R. & Minner, S. & Kiesmuller, G. 2002. A continuous time inventory model for a product
recovery system with options. International Journal of Production Economics 79: 121-141. Konstantaras, I. & Papachristos, S. 2007. Optimal policy and holding cost stability regions in a periodic review inventory system with manufacturing and remanufacturing options. European Journal of Operational Research 178: 433-448. Korugan, A. & Gupta, S. M. 1998. A Multi-Echelon Inventory System with Returns. Computers & Industrial Engineering 53 (1-2): 145-148.
Mabini, M.C. & Pintelon, L.M. & Gelders, L.F. 1992. EOQ type formulations for controlling repairable inventories. International Journal of Production Economics 28: 21-33.
Muckstadt, J.A. 1973. A model for a multi-item, multi-echelon, multi-indenture inventory system. Management Science 20 (4): 472-481.
Muckstadt, J.A. & Isaac. M.H. 1981. An analysis of single item inventory systems with returns. Naval Research Logistics Quarterly 28: 237-254.
Murayama, T. & Shu, L. H. 2001. Treatment of Reliability for Reuse and Remanufacture, Proceedings of the 2nd International Symposium on Environmentally Conscious Design and Inverse Manufacturing (EcoDesign'01), Tokyo, Japan.
Murayama, T. & Yoda, M. & Eguchi, T. & Oba, F. 2005. Adaptive Production Planning by Information Sharing for Reverse supply Chain, Proceedings of the 4th International Symposium on Environmentally Conscious Design and Inverse Manufacturing (EcoDesign'05), Tokyo, Japan.
Murayama, T. & Yoda, M. & Eguchi, T. & Oba F. 2006. Production Planning and Simulation for Reverse Supply Chain, Japan Society Mechanical Engineering International Journal, Series C, Vol. 49, No. 2.
Plewa, M. & Jodejko-Pietruczuk, A. 2011. The reverse logistics model of single-component product recovery. European Safety and Reliability Conference, Troyes, France.
Richter, K. 1994. An EOQ repair and waste disposal model. Proceedings of the Eight International Working Seminar on Production Economics 3: 83-91. Ingls/Innsbruck.
Richter, K. 1997. Pure and mixed strategies for the EOQ repair and waste disposal problem. OR Spectrum 19: 123-129.
Richter, K. & Sombrutzki, M. 2000. Remanufacturing planning for the reverse Wagrner/Whitin models. European Journal of Operational Research 121 (2): 304-315.
Richter, K. 1996a. The EOQ repair and waste disposal model with variable setup numbers. European Journal of Operational Research 95: 313-324.
Richter, K. 1996b. The extended EOQ repair and waste disposal model. International Journal of Production Economics 45: 443-447.
Richter, K. & Weber, J. 2001. The reverse Wagner/Whitin model with variable manufacturing and remanufacturing cost. International Journal of Production Economics 71: 447-456.
Schrady, D.A. 1967. A deterministic inventory model for repairable items. Naval Research Logistics Quarterly 14 (3): 391-398.
Sherbrooke, C.C. 1971. An Evaluator for the Number of Operationally Ready Aircraft in a Multilevel Supply System. Operations Research 19: 618- 635.
Sherbrooke, C.C. 1968. METRIC: A Multi-echelon Technique for Recoverable Item Control. Operations Research 16: 122-141.
Simpson, V.P. 1978. Optimum solution structure for a repairable inventory problem. Operations Research 26 (2): 270-281.
Teunter, R.H. 2001. Economic ordering quantities for recoverable item inventory system. Naval Research Logistics 48: 484-495.
Van der Laan, E.A. & Dekker, R. & Salomon, M. 1996. Product remanufacturing and disposal: A numerical comparison of alternative control strategies. International Journal of Production Economics 45: 489-498.
Van der Laan, E.A. & Dekker, R. & Salomon, M. & Ridder, A. 1996. An (S,Q) inventory model with remanufacturing and disposal. International Journal of Production Economics 46-47: 339350.
Van der Laan, E.A. & Salomon, M. 1997. Production planning and inventory control with
remanufacturing and disposal. International Journal of Operational Research 102: 264-278.
