Apiems2004 34.7 PDF
Apiems2004 34.7 PDF
Apiems2004 34.7 PDF
Using the former definition, one has
( ) ( )
y
xt
e t F 1 t R
= =
( )
( )
( )
1 y
xt
xt 1 y
xyt
e
e xyt
t R
t f
t Z
y
y
= = =
( )
1 y
xyt t Z
= will be used in the model development in this paper.
When , deteriorating rate increases with time; 1 y >
When , deteriorating rate decreases with time; and 1 y <
When , deteriorating rate is constant; the two-parameter Weibull distribution reduces to
the exponential distribution
1 y =
3. MODEL DEVELOPMENT
The manufacturers inventory system in Fig. 1 can be divided into two independent phases
depicted by to . This methodology reduces the complexity in our problem derivation
and analysis. Each phase has its own time unit, , which starts from the beginning of the
phase, . During time period, there is an inventory buildup and so there is deterioration
becomes effective. It increases to its maximum, , at
1
T
2
T
i
t
i
T
1
T
m
I
1 1
T t = . There is no production
during time period and the inventory level decreases due to demand and deterioration.
The inventory level becomes zero at
2
T
2 2
T t = .
[See APPENDIX Figure 1]
The manufacturers inventory system can be represented by the following differential
equations:
( )
( ) (
1 1
1
1
1
1 1
t I t d p
dt
t dI
=
)
1 1
T t 0 (1)
( )
( )
2 2
1
2
2
2 2
t I t d
dt
t dI
=
2 2
T t 0 (2)
34.7.5
Proceedings of the Fifth Asia Pacific Industrial Engineering and Management Systems Conference 2004
The first-order differential equations can be solved by using the boundary conditions
. From Spiegel1960, one has ( ) 0 0 I =
1 2 2
( ) 0 T I =
( )
( )
1
1
t
t
0
u
1 1
e
du e d p
t I
= ,
1 1
T t 0 (3)
( )
2
2
t
t
0
u
mv
2 2
e
du e d I
t I
= ,
2 2
T t 0 (4)
For small value, the approximate solution is derived by neglecting the second and
higher term of . From Fig. 1,
( ) 0 T I
2 2
= ; hence from equation (4),
du e d I
2
T
0
u
m
=
du
! m
u
d
2
T
0
0 m
m m
=
=
( )
=
+
+
=
0 m
1 m
2
m
1 m ! m
T
d
+
+
+
1
T
T d
1
2
2
(5)
From Fig. 1, ( ) ( )
m 2 1 1
I 0 I T I = = , one has
( )
2
T
0
u
1
T
1
T
0
u
du de
e
du e d p
Expand the exponential and the second or higher-order terms in are neglected, the
approximate solution is
( )
+
+ =
+
+ +
1
T
T d
1
T
T d p
1
2
2
1
1
1
(6)
The method used in Misra (1975), when 1 << where ( ) 1 T
1
1
+
+
is neglected,
results in the following approximate value for
1
T
34.7.6
Proceedings of the Fifth Asia Pacific Industrial Engineering and Management Systems Conference 2004
+
+
+
1
T
T
d p
d
T
1
2
2 1
(7)
From the fact , one has
2 1
T T T + =
+
+
1
T d
p
d p
T
T
2 2
(8)
The change in the retailers inventory level against time is depicted in Fig. 2.
[See APPENDIX Figure 2]
The retailers inventory system can be represented by the following differential equation:
( )
( ) t I t d
dt
t dI
r
1 r
=
,
n
T
t 0 (9)
Using the boundary condition, ( )
r r
I 0 I = , the first-order differential equation can solved..
From Spiegel (1960), one has
( )
t
t
0
u
r
r
e
du e d I
t I
= ,
n
T
t 0 (10)
For very small value, the approximate solution is derived by neglecting the second and
higher term of . From Fig. 2, 0
n
T
I
r
=
du
! m
u
d
n
T
0
0 m
m m
=
=
( )
=
+
+
=
0 m
1 m
m
1 m ! m
n
T
d
34.7.7
Proceedings of the Fifth Asia Pacific Industrial Engineering and Management Systems Conference 2004
+
+
1
n
T
n
T
d
1
(11)
For very small value, the approximate solution is derived by neglecting the second and
higher term of . The retailers total cost function per unit time,
r
TC , is the sum of the
ordering cost, the incoming quality cost, the holding cost and deteriorating. One has
r
TC
( ) ( ) ( )
+ + + =
d
n
T
0 I
T
n
C dt t I
T
n
C nK C
T
1
r r
n
T
0
r r 2 r 0 r 1
( )
+ + =
d
n
T
du e d
T
n
C dt
e
du e d du e d
T
n
C nK C
T
1
n
T
0
u
r
n
T
0
t
t
0
u
n
T
0
u
r 2 r 0 r 1
( )
( )
( )
+ + =
=
+
+
d
n
T
du e d
T
n
C
dt
! m
t
1 m ! m
t
n
T
d
T
n
C nK C
T
1
n
T
0
u
r
n
T
0
0 m
m
0 m
1 m
1 m
m
r 2 r 0 r 1
( )
( )( )
+ +
+ +
+ +
1
n
T
d
T
n
C
2 1
n
T
2
n
T
d
T
n
C nK C
T
1
1
r
2 2
r 2 r 0 r 1
(12)
The manufacturers total cost function per unit time, , is the sum of the setup cost,
the incoming quality cost, the holding cost and deteriorating. One has
m
TC
m
TC
( ) ( ) ( ) ( )
( )
+ + + =
d
n
T
0 I n dT pT
T
C
dt t I n dt t I dt t I
T
C
nK C
T
1
r 1
m
n
T
0
r 2
T
0
2 2 1
T
0
1 1
m 2
m 0 m 1
2 1
34.7.8
Proceedings of the Fifth Asia Pacific Industrial Engineering and Management Systems Conference 2004
( ) ( )
( )
( )
( )
( )
( )
( )
( )
+
+ + =
=
+
+
=
+ +
=
+
d
n
T
du e d n dT pT
T
C
dt
! m
t
1 m ! m
t
n
T
d n dt
! m
t
1 m ! m
t T
d
dt
! m
t
1 m ! m
t
d p
T
C
nK C
T
1
n
T
0
u
1
m
0 m
n
T
0
0 m
1 m
1 m
m
2
0 m
2
T
0
0 m
1 m
2
1 m
2
m
1
T
0
0 m
1
0 m
1 m
1
m
m 2
m 0 m 1
2
1
( ) ( )
( )( ) ( )( )
( )( )
+ +
+ +
+ +
+ +
+ +
+ +
+ +
1
n
T
d
n dT pT
T
C
2 1
n
T
2
n
T
d n
T
C
2 1
T
2
T
d
T
C
2 1
T
2
T
d p
T
C
nK C
T
1
1
1
v
2 2
v 2
2
2
2
2 v 2
2
1
2
1 v 2
v 0 v 1
(13)
The integrated total cost function TC for the manufacturer and the retailers is the sum of
and
m
TC
r
TC . From equations (7), (8), (12), and (13), the integrated total cost per unit time
can be written as a function of and n, one has
2
T
( )
m r 2
TC TC n , T TC + =
( )
( )( )
+ +
+ + =
+ +
1
n
T
d
T
n
C
2 1
n
T
2
n
T
d
T
n
C nK C
T
1
1
r
2 2
r 2 r 0 r 1
( ) ( )
( )( ) ( )( )
( )( )
+ +
+ +
+ +
+ +
+ + +
+ +
+ +
1
n
T
d
n dT pT
T
C
2 1
n
T
2
n
T
d n
T
C
2 1
T
2
T
d
T
C
2 1
T
2
T
d p
T
C
nK C
T
1
1
1
m
2 2
m 2
2
2
2
2 m 2
2
1
2
1 m 2
m 0 m 1
(14)
34.7.9
Proceedings of the Fifth Asia Pacific Industrial Engineering and Management Systems Conference 2004
If the retailers storage space limited to a maximum inventory size of W, the problem can
be stated as
Minimize: ( )
v d 2
TC TC n , T TC + =
Subject to: W fI g
r
=
The optimal inventory levels can be derived by using Lagrange multipliers, one has
( ) ( ) ( ) W fI n , T TC n , , T H
r 2 2
+ = (15)
4. SOLUTION PROCEDURE
The following heuristic technique (Yang and Wee, 2003) is used to derive the optimal n and
2
T
Step 1Since the number of delivery per order, n, is a discrete variable. Start by choosing a
discrete variable n, where n is any integer number equal or greater than 1.
Step 2Take the partial derivatives of ( ) n , , T H
2
with respect to and
2
T , and equating
the results to zero. The necessary conditions for optimality are
( )
0
T
n , , T H
2
2
=
and
( )
0
n , , T H
2
=
The simultaneous equations above can be solved for and
2
T .
Step 3Using found at step 2, substitute ( ( ) n T
2
( ) n T
2
,n) into equation (14) and find the
. ( ) ( ) n , n T TC
2
Step 4Repeat steps 2 and 3 for all possible n values until the minimum is
found. The
( ) ( ) n , n T TC
2
( ) ( )
* *
2
n , n T and ( ) ( )
* *
2
n , n T TC values constitute the optimal solution
and satisfy the following condition
( ) ( ) ( ) ( )
* *
2
* *
2
n , n T TC 1 n , 1 n T TC , and ( ) ( ) ( ) ( ) 1 n , 1 n T TC n , n T TC
* *
2
* *
2
+ +
Step 5Derive the and
*
T
1
*
T from equations (7) and (8), and the production quantity,
can be found.
1
pT
Step 6Derive the delivery quantity from equation (11).
If the objective function is convex, the following sufficient conditions must be satisfied:
34.7.10
Proceedings of the Fifth Asia Pacific Industrial Engineering and Management Systems Conference 2004
0
T
H H
T
H
2
2
2
2
2
2
2
>
, and 0
T
H
2
2
2
>
, 0
H
2
2
>
Since the total cost is a very complicated function due to high-power expression of the
exponential function, it is not possible to show analytically the validity of the above sufficient
conditions, a search procedure is used instead.
5. NUMERICAL EXAMPLE
The following parameters are assumed. The manufacturers production rate p is
2,000,000(unit/year), the demand rate d is 500,000(unit/year). The ordering cost
000 , 2 C
r 1
= ($/order), the production set-up cost 000 , 100 C
m 1
= ($/cycle), the
quality-control cost ($/delivery), the transportation charge
($/delivery), the carrying cost for the retailer and for the manufacturer are
500 K
r 0
=
000 , 1 K
m 0
=
60 C
r 2
= ($/unit/year) and ($/unit/year) respectively, the deteriorated cost for the
retailer and for the manufacturer are
40 C
m 2
=
600 C
r
= ($/unit) and 400 C
m
= ($/unit). The scale and
shape parameters of the deterioration rate are 1 . 0 = and 2 = . The retailers maximum
limited storage space for each unit is W2,000 feet
3
, the retailers storage requirement for
each unit is f2 feet
3
.
By applying the solution procedure in Section 4, the results are derived as shown in Table 1.
Table 2 is the comparison of results for the above special conditions.
[See APPENDIX Table 1]
[See APPENDIX Table 2]
From Table 1, if all the conditions and constraints are satisfied, optimal solution can be
derived. In this example, the integrated optimal solution which minimizes TC is found when
the number of delivery per order, n, is 59. The optimal , and T value are 0.0295,
0.0885, and 0.01180 year, respectively. The optimal production quantity of the manufacturer
is 59,000 units, and the optimal delivery quantity of the retailer is 1,000 units. The optimal
integrated total cost per year is $2,500.708. The concept of lot splitting comes from the idea
of single order with multiple deliveries. It is analyzed from the perspective of how material
flows in the system. This analysis highlights the value of the lot splitting. The number of
deliveries is inversely proportional to the lot size.
1
T
2
T
Table 2 shows the results comparisons for the above special conditions. From the above
34.7.11
Proceedings of the Fifth Asia Pacific Industrial Engineering and Management Systems Conference 2004
analysis, the following conclusions can be made
1. When the deterioration of item is not considered (i.e. 0 = ) when there is not
deterioration, the number of delivery decreases and the integrated total cost per year is
$8602 higher than the case when the deterioration is considered.
2. When the rate of deterioration is exponential rather than Weibull (i.e. 1 = ), the
integrated total cost per year is $128572 lower.
3. When the limited retailer storage space is not considered (i.e. 0 = ), the number of
delivery decreases and the integrated total cost per year is $585805 lower.
6. SUMMARY
Multi-lot size is one of the important policies of a successful enterprise. The study is
particularly useful for the inventory systems where the manufacturers and their retailers form
a strategic alliance with a mutually beneficial objective. To make it acceptable to both parties,
the integrated policy should offer some kind of profit sharing policy. The profit sharing policy
can be in the form of advanced payment and quantity discounts. As a result of this policy, both
the manufacturer and the retailer will benefit in the long run. The application of profit sharing
policy is an area worthy of studying for future research.
This study develops an optimal joint cost from the perspectives of both the manufacturer
and the retailer. The integrated two-stage inventory model with Weibull distribution
deteriorating items is assumed to have a constant demand rate and a limited retailer storage
space. A computer code is developed to derive the optimal solution, and a numerical example
is given to validate the results of the inventory system.
ACKNOWLEDGEMENTS
The authors wish to thank the referees for their constructive comments and suggestions on an
earlier version of the paper. This study is partially supported by the National Science Research
Council of the ROC.
REFERENCES
Abad, P. L., (2003), Optimal pricing and lot-sizing under conditions of perishability, finite
production and partial backordering and lost sale, European Journal of Operational
Research, 2003; 144: 677-685.
Conver, R. P., Philip, G. C. (1973), An EOQ model for items with Weibull distribution
deterioration, AIIE Transactions, 5, 323-326.
See APIEMS web site for reformating
Ghare, P. M., Schrader, G. F., (1963), A model for exponential decaying inventory, Journal
34.7.12
Proceedings of the Fifth Asia Pacific Industrial Engineering and Management Systems Conference 2004
of Industrial Engineering, 14, 238-243.
Goyal, S. K., (1976), An integrated inventory model for a single-supplier single-customer
problem, International Journal of Production Research, 15, 107-111.
Goyal, S. K., and Giri, B. C., (2003), The production-inventory problem of a product with
time varying demand, production and deterioration rates, European Journal of
Operational Research, 147, 549-557.
Goyal, S. K., Nebebe, F., (2000), Determination of economic production-shipment policy for
a single-vendor-single buyer system, European Journal of Production Research, 121,
175-178.
Hill, R. M., (1997), The single-vendor single-buyer integrated production-inventory model
with a generalized policy, European Journal of Production Research, 97, 493-499.
Misra, R. B. (1975), Optimal production lot size model for a system with deteriorating
inventory, International Journal of Production Research, 15, 495-505.
Misra, R. B. (1979). A note on optimal inventory management under inflation, Naval Logist
Quarterly, 26, 161-165.
Papachristos, S., and Skouri K., (2003), An inventory model with deteriorating items,
quantity discount, pricing and time-dependent partial backlogging, International Journal
of Production Economics, 83, 247-256.
Rau, H., Wu, M. Y., Wee, H. M., (2003), Integrated inventory model for deteriorating items
under a multi-echelon supply chain environment, International Journal of Production
Economics, 86, 155-168.
Sana, S. and Goyal, S. K. and Chaudhuri, K. S., (2004), A production-inventory model for a
deteriorating item with trended demand and shortages, European Journal of Operational
Research, 157, 357-371.
Skouri, K., and Papachristos, S., (2003), Optimal stopping and restarting production items
for an EOQ model with deteriorating items and time-dependent partial backlogging,
International Journal of Production Economics, 81-82, 525-531.
Spiegel, M. R., (1960), Applied Differential Equations, Prentice Hall, New York.
Wang, S. P., (2002), An inventory replenishment policy for deteriorating items with shortage
and partial backlogging, Computers and Operations Research, 29, 2043-2051.
Wee, H. M., (1997), A replenishment policy for item with a price-dependent demand and a
varying rate of deterioration, Production Planning and Control, 8.5, 494-499.
Wee, H. M., J . F. J ong, (1998), An integrated multi-lot-size production inventory model for
deteriorating items, Management System, 5.1, 97-114.
Wee, H. M., (1998), Optimal buyer-seller discount pricing and ordering policy for
deteriorating items, The Engineering Economist, 43.2, 151-168.
Wee, H. M., Law, S. T., (1999), Economic production lot size for deteriorating items taking
account of the time-value of money, Computers and Operations Research, 26, 545-558.
Wee, H. M., Law, S. T., (2001), Replenishment and pricing policy for deteriorating items
34.7.13
Proceedings of the Fifth Asia Pacific Industrial Engineering and Management Systems Conference 2004
taking account of the time-value of money, International Journal of Production
Economics, 71, 213-220.
Woo, Y. Y., Hsu, S. L., Wu, S., (2001), An integrated inventory model for a single vendor
and multiple buyers with ordering cost reduction, International Journal of Production
Economics, 73, 203-215.
Yang, P. C., Wee, H. M., (2002), A single-vendor and multiple-buyers production-inventory
policy for deteriorating item Computers and Operations Research, 143, 570-581.
Yang, P. C., Wee, H. M., (2003), An integrated multi-lot-size production inventory model for
deteriorating item Computers and Operations Research, 30, 671-682.
APPENDIX
( ) t I
m
m
I
0 t
1
t
2
t
1
T
2
T
Figure 1 Graphical representation for the manufacturers production system
34.7.14
Proceedings of the Fifth Asia Pacific Industrial Engineering and Management Systems Conference 2004
( ) t I
r
r
I
0
n
T
n
T 2
n
T 3
. t
Figure 2 Graphical representation for the buyers inventory system
Table 1 The numerical results
The time intervalyear
n
1
T
4
10
2
T
4
10
T
4
10
d
TC
( )
3
10
v
TC
( )
3
10
TC
( )
3
10
1 25862.487 5 15 20 1280.402 50494.983 51775.023
2 13105.000 10 30 40 7800.040 25509.967 26290.007
3 8847.518 15 45 60 613.373 17191.602 17804.975
20 1496.233 100 300 400 330.040 3278.645 3608.683
58 381.071 290 870 1160 297.281 2203.558 2500.839
59
*
366.232 295 885 1180 296.989 2203.719
2500.708
*
*
708 . 2500
60 351.638 300 900 1200 296.707 2204.357 2501.064
The optimal n value that minimizes the integrated TC
34.7.15
Proceedings of the Fifth Asia Pacific Industrial Engineering and Management Systems Conference 2004
Table 2 Comparison of results for special conditions
The time intervalyear
Special
conditions
No. of
delivery
n
1
T
4
10
2
T
4
10
T
4
10
d
TC
( )
3
10
v
TC
( )
3
10
TC
( )
3
10
Our example 59 295 885 1180 296.989 2203.719 2500.708
When 0 = 58 290 870 1160 297.241 2212.069 2509.310
When 1 =
61 305 914 1219 326.416 2045.720 2372.136
When 0 = 7 294 883 1177 301.791 1613.112 1914.903
34.7.16