Hala Aly Fergany and Naglaa Hassan El-Sodany
Page: 71-79 | Received 21 Sep 2022, Published online: 21 Sep 2022
Full Text Reference XML File PDF File
This study derives the probabilistic periodic review Backorders and Lost sales inventory models when the holding cost is a function of the inventory cycle. The expected total cost is composed of four components; the expected purchase cost, the expected ordering cost, the expected reviewing cost, the expected holding cost and the expected shortage cost. The objective is to minimize the expected annual total cost under a restriction on the expected annual reviewing cost when the protection interval demand follows the normal distribution. The Lagrangian multipliers are used to solve this constrained model in a closed form. Finally, some special cases are deduced and an illustrative numerical example is added with some graphs.
INTRODUCTION
Most of the literature dealing with probabilistic inventory models assumes that the demand is probabilistic since, the probability distribution of the future demand rather than the exact value of demand rate itself is known. Most of the probabilistic inventory models assume that the purchase cost of units is constant independent of the quantity ordered. Many researchers have studied unconstrained probabilistic inventory models assuming the holding cost to be constant independent of the number of periods. Hadley and Whitin (1963) and Taha (1997) have examined unconstrained probabilistic inventory models. Ben-Daya and Raouf (1994) examined unconstrained inventory model with constant units of cost, demand follows a normal distribution and the lead-time is one of the decision variables. Recently, Fergany and El-Saadani (2005) studied constrained probabilistic inventory model with continuous distributions and varying holding cost. Also, Abuo-El-Ata et al. (2003) has examined the constrained probabilistic multi-item inventory model.
Recently, El-Sodany (2011) has examined the periodic review probabilistic multi-item inventory system with zero lead time under constraint and varying holding cost. Liang et al. (2008) studied the periodic review inventory model with backorders and lost sales. Also, Lin (2008, 2010) studied the periodic review inventory models with backorders and variable lead time. The inventory model analyzed in this study assumes a periodic review system where demand is defined as a continuous random variable such that the inventory levels are reviewed at equal time intervals and orders are placed at such intervals, the quantity ordered each time depends on the available inventory level at the time of review.
One operating doctrine is that an order be placed at each review time if there have been any demands at all in the past period of time. A sufficient quantity is ordered to bring the inventory level up to level Qm where, Qm is the maximum inventory level after the arrival of the ordered quantity. The quantity ordered can vary from one review period to the next.
This operating doctrine will be called an order up to Qm policy. Inventory models which use an order up to Qm policy are referred to as (Qm, N) inventory models where, N is the time between reviews. In this study, we investigate the constrained periodic review (Qm, N) backorders and lost sales inventory models with varying holding cost.
Notations and assumptions: The following notations are adopted for developing the models:
D | = | The expected annual demand rate |
Qm | = | The maximum inventory level |
N | = | The inventory cycle |
1/N | = | The average number of cycles per year |
x | = | The continuous random variable represents the protection interval, L+ N, demand |
f (x; L+N) | = | The probability density function of the protection interval demand x in the protection interval L+N |
μ | = | The expected value protection interval demand |
σ | = | The standard deviation of the protection interval demand per year |
L | = | The lead time between the placement of an order and its receipt |
cr | = | The cost of making a review |
co | = | The order cost per cycle |
cp | = | The purchase cost of the item |
cl | = | The cost of a lost sale per cycle |
cb | = | The cost of a backorder per cycle |
ch | = | The holding unit cost per year |
Ch (N) = chNβ | = | The varying holding unit cost per year |
![]() |
= | The average on hand inventory; min on hand+max on hand/2 |
β | = | A constant real number selected to provide the benefit of estimated cost function |
Kr | = | The limitation on the expected annual review cost |
E (PC) | = | The expected annual purchase cost |
E (OC) | = | The expected annual ordering cost |
E (HC) | = | The expected annual holding cost |
E (BC) | = | The expected backorder cost |
E (LC) | = | The expected lost sales cost |
E (TC) | = | The expected annual total cost |
Min E (TC) | = | The minimum expected annual total cost |
The system is periodic review which means that the inventory levels are reviewed at equal time intervals and orders are placed at such intervals, the quantity ordered each time depends on the available inventory level at the time of review. An order is placed at each review time if there have been any demands at all in the past period of time, a sufficient quantity is ordered to bring the inventory level up to level Qm. The problem is to determine the optimal values of Qm and N which minimize the expected annual total cost. The following assumptions are made in the simple treatments for developing the mathematical model:
• | The inventory cycle N is defined as the time between the placement of two successive orders |
• | The average number of cycles per year can be written as 1/N |
• | The lead time L is constant |
• | The purchase cost cp of the item is constant independent of the quantity ordered |
• | There is never >1 order outstanding |
• | The cost cr of making a review is independent of the variables Qm and N |
• | The holding cost per unit is a varying function of the review time N. The varying holding unit cost per year takes the form: |
![]() |
• | The cost cb of a backorder is independent of the time at which the backorder exist |
Model I: Probabilistic periodic review (Qm, N) backorders inventory models under constrained reviewing cost and varying holding cost: In this model, the demands occurring when the system is out of stock are backordered until a replenishment quantity arrives. Backorders imply that when an order arrives, it is always sufficient to meet any out standing backorders. The relevant annual expected total cost is the sum of the expected purchase cost, expected holding cost, expected backorder cost, expected review cost and the expected ordering cost that is:
![]() |
Since, the purchase cost of the unit is constant independent of the quantity ordered; then the relevant annual expected total cost can be given by:
![]() |
(1) |
The expected annual reviewing cost is given by:
![]() |
(2) |
and the expected ordering cost is given by:
![]() |
(3) |
The expected annual holding cost will be found by computing the expected holding cost per period and then multiplying by the number of orders per year. The expected holding cost per period is where,
is the average on hand inventory per period given by:
![]() |
The expected holding cost per period:
![]() |
Then, the expected annual holding cost is given by:
![]() |
(4) |
The expected number of backorders incurred per year is the expected number of backorders incurred per period multiplied by the number of orders per year. Since, lead time L is constant then an order placed at time t will arrive in the system at time t+L and the next order quantity will arrive in the system at time t+L+N. After an order placed at time t, the inventory position of the system will be Qm, hence to compute the expected number of backorders occurring between t+L and t+L+N, a backorder will occur in this period of time if and only if the protection interval demand in protection interval L+N exceeds Qm. Then the expected number of backorders incurred per period is given by:
![]() |
(5) |
The expected number of backorders incurred per year is given by:
![]() |
(6) |
Then, the expected annual cost of backorders is given by:
![]() |
(7) |
Then the expected annual total cost is given by:
![]() |
(8) |
The aim is to determine the optimal values of Qm and N that minimize the expected annual total cost under the following constraint:
![]() |
In the study, we will consider Eq. 8 is a convex function. Let us write it in the following form:
![]() |
(9) |
Subject to:
![]() |
(10) |
To find the optimal values Q*m and N* which minimize Eq. 9 under Eq. 10, we use the Lagrangian multiplier technique with the Kuhn-Tucker conditions as follows: The Lagrangian function is:
![]() |
(11) |
where, λ is the Lagrangian multiplier. The optimal values Q*m and N* are found by setting the corresponding 1st partial derivatives of L(Qm, N) equal to zero as follows:
![]() |
Hence, the optimal maximum inventory level is the solution of the following equation:
![]() |
(12) |
also:
![]() |
Hence:
![]() |
(13) |
Substituting from Eq. 12 into Eq. 13 then the optimal inventory cycle is the solution of the following equation:
![]() |
(14) |
However, and N*>0 minimize the expected annual total cost (Eq. 8) since;
![]() |
(15) |
![]() |
(16) |
![]() |
(17) |
Therefore from the Hesian matrix, we have:
![]() |
(18) |
Clearly, there is no closed form solution of Eq. 12 and 14. But by the following algorithm due to Hadley and Whitin (1963), we can obtain a closed approximate solution of these equations in a finite number of iterations as follows:
Step 1: Assume that the initial value of the inventory cycle N is any constant number then from Eq. 12, we have the initial maximum inventory level Qm1.
Step 2: Substituting by Qm1 into Eq. 14 to get N1.
Step 3: Substituting by N1 into Eq. 12 to get Qm2. The procedure is to vary λ in steps 2-3 until the smallest value of λ>0 is found such that the constraint holds for the different values of β.
Step 4: Repeat steps 2 and 3 until successive values of Qm and N are sufficiently close which are the optimal values that gives the minimum annual expected total cost numerically.
Model I with normally distributed protection interval demand: Assume that the protection interval demand follows the Normal distribution with mean μ = D (L+N) and standard deviation σ/L+N. So, we can minimize the expected annual total cost mathematically as follows consider:
![]() |
(19) |
Then from Eq. 12, the optimal maximum inventory level is the solution of the following equation:
![]() |
(20) |
and also from Eq. 14, the optimal inventory cycle is the solution of the following equation:
![]() |
(21) |
and the minimum expected annual total cost is given by:
![]() |
(22) |
Clearly, there is no closed form solution of Eq. 20 and 21 and we minimize the expected annual total cost numerically using the iteration method as the previous algorithm.
Model II: Probabilistic periodic review (Qm, N) lost sales inventory models under constrained reviewing cost and varying holding cost: In the backorders case, it was assumed that all demands incurred when the system was out of stock were backordered but in the lost sales case a demand which occurs when the system is out of stock is lost forever.
We shall assume that the cost of a lost sale has the form cl that the lost sale cost will be constant independent of time since a demand which occurs when the system is out of stock in the lost sales case is lost forever.
The relevant annual expected total cost is the sum of the expected purchase cost, expected review cost, expected ordering cost, expected holding cost and the expected lost sale cost that is:
![]() |
Since, the purchase cost of the unit is constant independent of the quantity ordered then the relevant annual expected total cost can be given by:
![]() |
(23) |
The expected annual costs in the lost sales case are the same as the expected annual costs of the backorders case except the holding cost. For the lost sales case, the holding cost will be affected by the extra stock, since the demand which occurs when the system is out of stock is lost forever that there will be an increase in the holding cost.
The expected annual holding cost will be found by computing the expected holding cost per period and then multiplying by the number of orders per year. The expected holding cost per period is ch(N) . Where,
is the average on hand inventory per period given by:
![]() |
Then, the expected holding cost per period can be given by:
![]() |
The expected annual holding cost is given by:
![]() |
(24) |
The expected annual lost sales cost is given by:
![]() |
(25) |
Then the expected annual total cost is given by:
![]() |
(26) |
The aim is to determine the optimal values of Qm and N that minimize the expected annual total cost under the following constraint:
![]() |
In the study, we will consider Eq. 26 is a convex function. Let us write it in the following form:
![]() |
(27) |
Subject to:
![]() |
(28) |
To find the optimal values Q*m and N* which minimize Eq. 27 under Eq. 28, we use the Lagrangian multiplier technique with the Kuhn-Tucker conditions as follows: The Lagrangian function is:
![]() |
(29) |
The optimal values: Q*m and N* are found by setting the corresponding 1st partial derivatives of L(Qm, N) equal to zero as follows:
![]() |
Hence, the optimal maximum inventory level is the solution of the following equation:
![]() |
(30) |
also:
![]() |
Hence, the optimal inventory cycle is the solution of the following equation:
![]() |
(31) |
However, Q*m> 0 and N*>0 minimize the expected annual total cost (Eq. 27) since:
![]() |
(32) |
![]() |
(33) |
![]() |
(34) |
Therefore from the Hesian matrix, we have:
![]() |
![]() |
(35) |
Clearly, there is no closed form solution of Eq. 30 and 31 and we minimize the expected annual total cost numerically using the iteration method as the algorithm used in model I.
Model II with normally distributed protection interval demand: Assume that the protection interval demand follows the Normal distribution with mean μ = D(L+N) and standard deviation σ/L+N with the probability density function given by Eq.19. Then from Eq. 30, the optimal maximum inventory level is the solution of the following equation:
![]() |
(36) |
and from Eq. 31, the optimal inventory cycle is the solution of the following equation:
![]() |
(37) |
and the minimum annual expected total cost is given by:
![]() |
(38) |
Clearly, there is no closed form solution of Eq. 36 and 37 and we minimize the expected annual total cost numerically using the iteration method as the algorithm used in model I.
Special cases: We deduce two special cases of the models as follows:
Case 1: For model I, assume that:
![]() |
Thus from Eq. 20, the following expression for the optimal maximum inventory level can be obtained:
![]() |
Also from Eq. 21, the following expression for the optimal inventory cycle can be obtained:
![]() |
This is unconstrained probabilistic periodic review (Qm, N) backorders inventory model with protection interval demand follows the normal distribution and the holding unit cost is constant.
Case 2: For model II, assume that:
![]() |
Thus from Eq. 36, the following expression for the optimal maximum inventory level can be obtained:
![]() |
Also from Eq. 37, the following expression for the optimal inventory cycle can be obtained:
![]() |
This is unconstrained probabilistic periodic review (Qm, N) lost sales inventory model with protection interval demand follows the normal distribution and the holding unit cost is constant.
An illustrative example: A large California warehouse follows a policy of reviewing all items over fixed periods of time. It uses an order up to Qm policy. Consider one item, it carries say a particular type of tractor tire. The mean demand rate has been constant over time at the value of 600 year-1.
The warehouse orders tires directly from the manufacturer and the lead time L is nearly constant and has the value 6 months.
The demand in the protection interval L+N can be represented quite well by a normal distribution with mean 450 unit and standard deviation 25.981 unit. The warehouse uses an inventory holding cost is $3. The ordering cost is $13 and the reviewing cost is $12. The researchers will determine the optimal maximum inventory level and the optimal inventory cycle that minimize the expected annual total cost in the following cases:
• | All demands occurring when the system is out of stock are backordered and the cost of a backorder is $25, under the expected review cost limitation Kr = $44.5 |
• | All demands occurring when the system is out of stock are lost for ever and the cost of a lost sale is $25, under the expected review cost limitation Kr = $44.3 |
Solution: Using the results of the models, the optimal values and the minimum annual expected total cost for the periodic review (Qm, N) inventory models with varying holding cost can be shown in Table 1 and 2.
Table 1: | The optimal solutions and the min E (TC) for the backorders periodic review (Qm, N) inventory model when the protection interval demand follows the normal distribution |
![]() |
|
Table 2: | The optimal solutions and the min E (TC) for the lost sales periodic review (Qm, N) inventory model when the protection interval demand follows the normal distribution |
![]() |
|
![]() |
|
Fig. 1: | The values of min E (TC) for the backorders periodic review (Qm, N) inventory model when the protection interval demand follows the normal distribution at each value of β |
![]() |
|
Fig. 2: | The values of min E (TC) for the lost sales periodic review (Qm, N) inventory model when the protection interval demand follows the normal distribution at each value of β |
And the solution of the problem may be determined more readily by plotting min E (TC) against β for the two inventory shows as in Fig. 1 and 2.
MATERIALS AND METHODS
The aim of this study is to investigate probabilistic periodic review (Qm, N) backorders inventory models and probabilistic periodic review (Qm, N) lost sales inventory models under the expected reviewing cost constraint when the holding cost is a varying function of the inventory cycle.
The lagrangian multiplier and probabilistic periodic review (Qm, N) lost sales inventory models under the expected reviewing cost constraint when the holding cost is a varying function of the inventory cycle. The lagrangian multiplier approach is used to determine the optimal inventory cycle and the optimal maximum inventory level which minimize the expected annual total cost under the expected review cost constraint.
RESULTS AND DISCUSSION
The basic results of this study are: The minimum expected annual total cost of the probabilistic periodic review (Qm, N) backorders inventory model under constrained reviewing cost and varying holding cost when the protection interval demand follows the normal distribution is given by:
![]() |
The minimum expected annual total cost of the probabilistic periodic review (Qm, N) lost sales inventory model under constrained reviewing cost and varying holding cost when the protection interval demand follows the normal distribution is given by:
![]() |
At the end of this study, special cases of previously research are added. Also, a numerical illustrative example is added with some graphs by using Mathematica program.
CONCLUSION
We have evaluated the optimal inventory cycle and the optimal maximum inventory level; then we deduced the minimum annual expected total cost min E (TC) of the considered backorders periodic review inventory model and lost sales periodic review inventory model. We draw the curves min E (TC) against β for both models which indicate the value of β that gives the minimum value of the expected annual total cost of our numerical example.
Hala Aly Fergany and Naglaa Hassan El-Sodany. Probabilistic Periodic Review (Qm, N) Backorders and Lost Sales Inventory Models under Constraint and Varying Holding Cost and Normally Distributed Protection Interval Demand.
DOI: https://doi.org/10.36478/jmmstat.2011.71.79
URL: https://www.makhillpublications.co/view-article/1994-5388/jmmstat.2011.71.79