Decision Support System For Periodic Review Inventory Model With Backorder Price Discounts Under Controllable Lead Time
Decision Support System For Periodic Review Inventory Model With Backorder Price Discounts Under Controllable Lead Time
Decision Support System For Periodic Review Inventory Model With Backorder Price Discounts Under Controllable Lead Time
1
( ).
i
j j
j
U u
=
Thus,
the lead time crashing cost R(L) per replenishment cycle is given by
R(L) = a
i
(L
i1
L) +
1
1
( )
i
j j j
j
a U u
, for L e [L
i
, L
i1
].
4. The backorder ratio | is variable and is in proportion to the backorder price
discount offered by the supplier per unit t
x
; thus, | = |
0
t
x
/t
0
, for 0 < |
0
s 1, 0 s
t
x
s t
0
(Pan & Hsiao, 2001).
DECISION SUPPORT SYSTEM
Business and Information 2012
(Sapporo, July 3-5)
F 346
Sprague & Carson (1982) thought a complete framework of decision support
system should include the following events shown as figure 1:
Data BaseStore up the data that the decision support system requires.
Model BaseStore up the data analysis model and decision analysis model or
view these models as sub-patterns which are independent with each other.
Software SystemThe software system includes the following three subsystems.
Data Base Management Subsystem, DBMSs.
Model Base Management Subsystem, MBMSs.
Dialog Generation and Management Subsystem, DGMSs.
System Analysis
Data base Model base
MBMSs
DGMSs
DBMSs
Figure 1 DSS System Framework
Business and Information 2012
(Sapporo, July 3-5)
F 347
System Framework
Setup Cost Data
Ordering Cost
Data
Backorder Data
Crashing Cost
Data
Holding Cost
Data
Normal Demand
Cost Model
Distrbution Free
Demand Cost
Model
Uncertainty factors
Database system
Model management
system
User platform
Cost data
Model data
Figure 2 System framework
NORMAL DEMAND MODEL
The problem under study is to minimize the following total expected annual cost:
EAC(T, t
x
, k, L) = OC + HC + SC + CC
where OC stands for ordering cost, HC for holding cost, SC for shortage cost,
and CC for crashing cost. Consequently, the total expected annual cost EAC(T, t
x
, k,
L) is
EAC(T, t
x
, k, L) )] ( T ) 1 ( T
2
DT
[
T
A
k L L k h + + + + + + = o | o
) ( ] ) 1 ( [
1
0
k L T
T
x
+ + + + o t | |t
Business and Information 2012
(Sapporo, July 3-5)
F 348
(
+ +
1
1
1
) ( ) (
1
i
j
j j j i i
t T a L L a
T
for L
i
s L s L
i1
(1)
Taking the first partial derivatives of EAC(T, t
x
, k, L) with respect to T t
x
, k and
L, respectively, it follows that
2 0
0 0
2
0
EAC( , , , )
2
( )
2
x
x x
T
T L
T k L hk
T L
k
T T T L
t | o
t t | t o
t
+
( c
+
= + +
(
c +
1
1 2 2
1
1
( ) ( )
T 2 T
i
i i i j j
j
A hD
a L L a T t
=
(
+ +
(
(2)
0 0
x 0
0 0
EAC( , , , ) 2 1
= ( ) ( )
x
x
T k L
h T L k T L k
T
t | |
t | o
t t t
( c
+ + + + +
(
c
(3)
0
0
2 0
0 0
0
1
EAC( , , , )
(1- ( ))
x
x
x x
h
T k L
h T L T L k
k
T
|
t
t
t
o o
|
t t | t
t
(
| |
( |
c \ .
(
= + + u
(
c
| |
(
+ +
|
(
\ .
(4)
( ) ( ) ( )
0
1 1
0
2 2
2 0
0
0
1
EAC( , , , ) 1 1 1
4 2
1
x
x
x x
h
T k L
L T k hk L T
L T
T
|
t
t
t
o o
|
t |t
t
(
| |
( |
c \ .
(
= + + + +
(
c
| |
(
+
|
(
\ .
(5)
Setting equation (4) to zero and solving for t
x
, it follows that
0
2
x
hT t
t
+
= (6)
Setting equation (3) to zero and substituting (6) into (3) to solve for T, thus
( ) ( ) ( ) 0
4
1 2
2
1
0
0
0
0 2
=
(
+ + + |
.
|
\
|
+
|
|
.
|
\
|
+ + L R k L T A k L T
D
h
Dh T o
|
t o
t
|
(7)
Setting equation (4) to zero and substituting (6) into (4) to solve for k, then
( )
(
|
|
.
|
\
|
|
.
|
\
|
+ +
= u
4
1
4
1
)
2 2
1 (
1 ) (
0
0
2
0
0 0
0
0
|
t
t
| |
t
|
hT
T
hT
h
h
k
(8)
It can be shown from (8) that hT s t
0
, so u(k) > 0.5 holds for nonnegative safety
Business and Information 2012
(Sapporo, July 3-5)
F 349
factor k. Consequently, the value of t
x
derived in (6) will automatically be positioned
between 0 and t
0
as specified in assumption 4.
DISTRIBUTION FREE MODEL
In a distribution free model, the lead time demand assumption is relaxed to any
distribution function by only assuming that the distribution function F belongs to the
class with finite mean DL and variance o
2
L. Again, the expected total annual cost
of the model can be expressed as EAC(T, t
x
, k, L) = OC + HC + SC + CC. While the
expected ordering and expected crashing costs are of the same forms as those of the
normal demand model, the expected shortage at the end of a cycle can be expressed
by B(r) = E[X r]
+
the expected total annual cost can be presented by
0
1
1
1
DT
EAC( , , , ) T (1 ) [ ]
T 2
1
[ (1 ) ] [ ]
T
1
( ) ( )
T
x
x
i
i i j j j
j
A
T k L h k L E X r
E X r
a L L a T t
t o |
|t | t
+
+
=
(
= + + + +
(
+ +
(
+ +
(
for L
i
s L s L
i-1
(9)
Since the form of the lead time demand distribution is not known, the minimax
criterion is applied to find the least favorite distribution function inOfor each (T, t
x
, k,
L) and then find the optimal values of T
*
, t
x
*
, k
*
and L
*
that minimize the expected total
annual cost. In mathematical symbolization, the problem under investigation can be
expressed as
O e > > F L Q
max min
0 , 0
EAC(T, t
x
, k, L)
Gallego & Moon (1993) proved that the following inequality holds for any F e
O,
| | ) ) (T ( ) ) (T ( ) (T
2
1
) ( E
2 2
L D r L D r L r X + + + + s
+
o (10)
Substituting r = DL + ko L into the model, the problem can be reduced to
minimize
1
1
1
0
0
2
2 0
0 0
0
A DT 1
EAC( , , , ) T ( ) ( )
T 2 T
1
1
T ( 1 )
2
1
T
i
x i i j j j
j
x
x x
T k L h k L a L L a T t
h
L k k
t o
|
t
t
o
|
t t | t
t
=
(
(
= + + + + +
(
(
(
(
+ + +
`
(
+ +
(
)
for L
i
s L s L
i1
(11)
Business and Information 2012
(Sapporo, July 3-5)
F 350
Taking the partial derivatives of EAC(T, t
x
, k, L) with respect to T, t
x
, k and L in
each time interval [L
i
, L
i1
] in (11), we obtain
( )
2
2 2 0
0 0
2
0
1
1
2
1
EAC( , , , )
T 2 2
2
1
2
1
( ) ( )
T
x
x x
i
i i i j j
j
T k L A hD hk
T T L
T
T L
T L
k k
T
a L L a T t
t o
|
t t | t o
t
=
c
= + +
c +
+
(
+
+ +
(
(
+
(
(12)
) , , , EAC(
x
x
L k T
t
t
c
c
= ( ) k k L T
T
k L T h + +
(
+ + +
2
0 x
0
0
0
0
1
2
2
1
) ( o | t
t
|
t
|
(13)
0
0
2 2 0
0 0
0
1
EAC( , , , ) 1
(1- )
2
1
x
x
x x
h
T k L k
h T L T L
k
k
T
|
t
t
t
o o
|
t t | t
t
(
| |
( |
\ .
(
c
(
= + +
| |
c
(
+ +
|
(
\ .
+
(
(14)
( )
( )
( )
0
0
1
0
2
2
2
0
0
1
2
1
EAC( , , , ) 1
1
8
1
1 1
2
x
x
x x
h
T k L
L T k k
L
T
hk L T
T
|
t
t
t
o
|
t | t
t
o
(
| |
( |
c \ .
(
= + +
(
c
| |
(
+
|
(
\ .
+ +
(15)
Setting equation (13) to zero and solving for t
x
, it follows that
0
2
x
hT t
t
+
= (16)
Setting equation (12) to zero and substituting (16) into (12) to solve for T, it
follows that
( )
( )
( )
2 2 0
0
2 0 0
1 1
4
2 1 1 0
2 4
h
T Dh T L k k
D
A T L k k R L
|
o
t
t |
o
| |
+ +
|
\ .
(
| |
+ + + + =
| (
\ .
(17)
Business and Information 2012
(Sapporo, July 3-5)
F 351
Solving for k by setting equation (14) to zero and substituting (16) into (14), we have
( )
(
(
|
|
.
|
\
|
|
|
.
|
\
|
+ +
=
+
4
1
4
1
)
2 2
1 (
2
1
1
0
0
2
0
0 0
0
0
2
|
t
t
| |
t
|
hT
T
hT
h
h
k
k
(18)
The following algorithm can be used to find the optimal values of the order
quantity, backorder discount, reorder point and lead time.
Step 1. For i = 0, 1, 2, , n
(i) Set k
io
= 0
(ii) Substitute k
io
into (17) to evaluate T
io
.
(iii) Use T
io
to determine k
in
from (18), then find k
in
. Let k
io
= k
in
.
(iv) Repeat (ii) and (iii) until no change occurs in the values of T
i
and k
i
.
Denote these resulting solutions by T
i
and k
i
.
(v) Use T
i
and equation (16) to compute the backorder price discount t
xi
.
(vi) Use equation (1) to compute the expected total annual cost EAC(T
i
, t
xi
,
k
i
, L
i
).
Step 2. Set EAC(T
*
, t
x
*
, k
*
, L
*
) = Min{EAC(T
i
, t
xi
, k
i
, L
i
), i = 0, 1, 2, , n}
Step 3.(T
*
, t
x
*
, k
*
, L
*
) is a set of optimal solutions.
Numerical Example
Consider an inventory system with the following data: D = 600 units/year, A =
$200 per order, h = $20 per unit per year, t
0
= $150 per unit, o = 7 unit/week, and the
lead time has three components with data shown in Table 1 (Pan & Hsiao, 2001).
Apply the proposed algorithm to solve the problem for the upper bound of the
backorder ratio |
0
= 0.95, 0.80, 0.65, 0.50, 0.35 and 0.20, respectively. The resulting
optimal solutions are summarized in Table 2 and Table 3.
Table 1. Lead time data of the examples
Lead time component i 1 2 3
Normal duration U
i
(days) 20 20 16
Minimum duration u
i
(days) 6 6 9
Unit fixed crashing cost a
i
($/day)
0.4 1.2 5.0
Business and Information 2012
(Sapporo, July 3-5)
F 352
Table 2. Summary of the results for Example in Normal Demand (
i
L in weeks)
L
0
|
i 0 1 2 3
L
i
8 6 4 3
0.2
T
24.75 23.41 21.85 21.14
t
x
79.76 79.50 79.20 79.07
k
1.53 1.53 1.53 1.53
EAC(T, t
x
, k, L)
4916.69 4753.27 4596.42 4582.02
0.35
T
24.32 23.02 21.50 20.83
t
x
79.68 79.43 79.13 79.01
k
1.52 1.52 1.52 1.52
EAC(T, t
x
, k, L)
4884.67 4726.44 4576.32 4566.46
0.50
T
23.89 22.62 21.15 20.50
t
x
79.59 79.35 79.07 78.94
k
1.50 1.50 1.50 1.50
EAC(T, t
x
, k, L)
4852.73 4699.84 4556.59 4551.34
0.65
T
23.45 22.22 20.79 20.18
t
x
79.51 79.27 79.00 78.88
k
1.49 1.49 1.49 1.49
EAC(T, t
x
, k, L)
4820.93 4673.50 4537.29 4536.70
0.80
T
23.01 21.81 20.42 19.85
t
x
79.42 79.19 78.93 78.82
k
1.47 1.47 1.47 1.47
EAC(T, t
x
, k, L)
4789.29 4647.49 4518.47 4522.60
0.95
T
22.55 21.39 20.05 19.51
t
x
79.34 79.11 78.86 78.75
k
1.45 1.45 1.45 1.45
EAC(T, t
x
, k, L)
4757.88 4621.86 4500.20 4509.10
Business and Information 2012
(Sapporo, July 3-5)
F 353
Table 3 Summary of the results for Example in Distribution Free Demand (
i
L in
weeks)
L
0
|
i 0 1 2 3
L
i
8 6 4 3
0.2
T
27.32 25.80 23.98 23.09
t
x
80.25 79.96 79.61 79.44
k
2.55 2.62 2.72 2.77
EAC(T, t
x
, k, L)
5895.09 5634.61 5351.24 5250.03
0.35
T
26.84 25.35 23.58 22.72
t
x
80.16 79.88 79.53 79.37
k
2.55 2.62 2.71 2.76
EAC(T, t
x
, k, L)
5844.37 5589.60 5313.98 5218.28
0.50
T
24.82 23.48 21.91 21.20
t
x
79.77 79.52 79.21 79.08
k
2.52 2.58 2.65 2.69
EAC(T, t
x
, k, L)
5635.76 5405.93 5163.83 5091.59
0.65
T
25.85 24.43 22.76 21.97
t
x
79.97 79.70 79.38 79.23
k
2.54 2.60 2.69 2.73
EAC(T, t
x
, k, L)
5741.23 5498.49 5239.08 5154.82
0.80
T
25.34 23.96 22.34 21.59
t
x
79.87 79.61 79.30 79.15
k
2.53 2.59 2.67 2.71
EAC(T, t
x
, k, L)
5688.79 5452.39 5201.49 5123.16
0.95
T
24.82 23.48 21.91 21.20
t
x
79.77 79.52 79.21 79.08
k
2.52 2.58 2.65 2.69
EAC(T, t
x
, k, L)
5635.76 5405.93 5163.83 5091.59
Business and Information 2012
(Sapporo, July 3-5)
F 354
CONCLUSIONS
This paper extends the Pan and Hsiao model (Pan & Hsiao, 2001) to study the
impact of safety factor on the periodic review inventory model involving controllable
lead time and backorder price discount with mixture of backorder and partial lost
sales. The objective is to minimize the expected total annual cost by simultaneously
optimizing order period, backorder price discount, safety factor and lead time. There
are two models considered in the plan, one with normal demand distribution and
another with general demand distribution having both mean and variance known and
finite. Algorithms are developed to find the optimal solutions for the two models and
examples are provided to illustrate the procedures of the algorithms. We construct a
decision support system to analysis the inventory model, then find out the best choice
with period review time, backorder discount and lead time. Once we build the
decision support system, we can use the result to apply into other inventory model and
adjust the parameter to fit the other different inventory situations. The calculation of
the expected value of additional information for the distribution free problem are
demonstrated in the examples as well.
Business and Information 2012
(Sapporo, July 3-5)
F 355
REFERENCES
Ben-Daya, M. & Raouf, A. 1994. Inventory models involving lead time as a decision
variable. Journal of the Operation Research Society. 45, 579-582.
Blackburn, J. D. 1991. Time-Based Competition: The Next Battleground in American
Manufacturing. Homewood, Illinois: Business One Irwin.
Bockerstette, J. A. & Shell, R. L. 1993. Time Based Manufacturing, New York:
McGraw-Hill.
Gallego, G. & Moon, I. 1993. The distribution free newsboy problem: Review and
extensions, Journal of the Operational Research Society, 44, 825-834.
Hariga, M. A. & Ben-Daya, M. 1999. Some stochastic inventory models with
deterministic variable lead time. European Journal of Operational Research. 113,
42-51.
Liao, C. J. & Shyu, C. H. 1991. Analytical determination of lead time with normal
demand. International Journal of the Operational Research and Production
Management, 11, 7278
Moon, I. & Choi, S. 1998. A Note on lead time and distribution assumptions in
continuous review inventory models. Computers and Operations Research. 25,
1007-1012.
Naddor ,N. 1996. Inventory System. New York: John Wiley.
Ouyang, L. Y. Yen, N. C. & Wu, K. S. 1996. Mixture inventory model with
backorders and lost sales for variable lead time. Journal of the Operation
Research Society, 47, 829-832.
Ouyang, L. Y. & Wu, K. S. 1998. A minimax distribution free procedure for mixed
inventory model with variable lead time. International Journal of Production
Economics. 56, 511516.
Ouyang, L. Y. Chuang, B.R. & Lin, Y. J. 2003. Impact of backorder discounts on
periodic review inventory model. International Journal of Information and
Management Sciences. 14(3), 1-13.
Pan, J. C. & Hsiao, Y. C. 2001. Inventory models with backorder discounts and
variable lead time. International Journal of Systems Science. 32, 925929.
Silver, E. A. Pyke, D. F. & Peterson, R. 1998. Inventory Management and Production
Planning and Scheduling, 3/e. New York: John Wiley.
Sprague, R H. & Carson, E. D. 1982. Building Effective Decision Support System.
Prentice-Hall.
Tersine, R. J. 1994. Principles of Inventory and Material Management. New York:
North Holland.
Vollmann, T.E., Berry, W.L. & Whybark, D.C. 1992. Manufacturing Planning and
Control System 3/e. Chicago: Irwin.
Business and Information 2012
(Sapporo, July 3-5)
F 356