# A New Mathematical Inventory Model with Stochastic and Fuzzy Deterioration Rate under Inflation

Authors: Bahar Naserabadi; Abolfazl Mirzazadeh; Sara Nodoust
Journal: Chinese Journal of Engineering
 (2014)

---

## Abstract

This paper develops an inventory model for items with uncertain deterioration rate, time-dependent demand rate with nonincreasing function, and allowable shortage under fuzzy inflationary situation. The goods are not deteriorating upon reception, but the deteriorating starts after elapsing a specified time. The lead time and inflation rate are both uncertain in the model. The resultant effect of inflation and time value of money is assumed to be fuzzy in nature and also we consider lead time as a fuzzy function of order quantity. Furthermore the following different deterioration rates have been considered: for the first case we consider fuzzy deterioration rate and for the second case we assume that the deterioration rate is time dependent and follows Weibull distribution with three known parameters. Since the inflation rate, deterioration rate, and the lead time are fuzzy numbers, the objective function becomes fuzzy. Therefore the estimate of total costs for each case is derived using signed distance technique for defuzzification. The optimal replenishment policy for the model is to minimize the total present value of inventory system costs, derived for both the above mentioned policies. Numerical examples are then presented to illustrate how the proposed model is applied.


---

## Body

## 1. Introduction

Many of the physical goods undergo decay or deterioration over time. Commodities such as fruits, vegetables, and foodstuffs are subject to direct spoilage during storage period. The highly volatile liquids such as gasoline, alcohol, and turpentine undergo physical depletion over time through the process of evaporation. The electronic goods, radioactive substances, photographic film, and grain deteriorate through a gradual loss of potential or utility with passage of time. The deteriorated items repairing is a major problem in the supply chain of most of the business organizations.The first attempt to describe optimal ordering policies for deteriorating items was made by Ghare and Schrader [1]. Later, Covert and Philip [2] derived the model with variable deteriorating rate of two-parameter Weibull distribution.Inventoried goods can be broadly classified into four metacategories:(1)
obsolescence which refers to items that lose their value through time due to rapid changes of technology or the introduction of a new product by a competitor;
(2)
deterioration which refers to the damage, spoilage, dryness, vaporization, and so forth of the products;
(3)
amelioration which refers to items whose value or utility or quantity increases with time;
(4)
no obsolescence/deterioration/amelioration.If the rate of obsolescence, deterioration, or amelioration is not sufficiently low, its impact on modelling of such an inventory system cannot be ignored. There are a few papers for obsolescing and ameliorating items. Moon et al. [3] considered the ameliorating/deteriorating items on an inventory model with time-varying demand pattern. Against obsolescing and ameliorating items, the deteriorating inventory models under inflationary conditions are studied greatly. In a few of these works, deterioration rate is not constant. For example, Chen [4] proposed an inflationary model with time proportional demand and Weibull distribution for deteriorating items using dynamic programming. Sana [5] proposed a model which deals with a stochastic economic order quantity (EOQ) model over a finite time horizon where uniform demand over the replenishment period is price dependent. The selling price is assumed to be a random variable that follows a probability density function. Moreover, the article suggests a new function regarding price-dependent demand. Sana [6] also presented another EOQ model over an infinite time horizon for perishable items where demand is price dependent and partial backorder is permitted. Based on the partial backlogging and lost sale cases, the author develops the criterion for the optimal solution for the replenishment schedule and proves that the optimal ordering policy is unique. Balkhi [7] presented a production lot-size inventory model where the production, demand, and deterioration rates are known, continuous, and differentiable functions of time. Shortages are allowed, but only a fraction of the stockout is backordered, and the rest is lost. De and Sana [8] developed a research that deals with a backorder EOQ model with promotional index for fuzzy decision variables. Here, a profit function is developed, where the function itself is the function of mth power of promotional index (PI) and the order quantity; the shortage quantity and the PI are the decision variables. The demand rate is operationally related to PI variables and the model has been split into two types for the multiplication and addition operations. Lo et al. [9] introduced an integrated production-inventory model with imperfect production processes and Weibull distribution deterioration under inflation from the perspectives of both the manufacturer and the retailer. Sana [10] developed an EOQ model to determine the retailer’s optimal order quantity for similar products. It is assumed that the amount of display space is limited and the demand of the products depends on the display stock level and the initiatives of sales staff where more stock of one product makes a negative impression of the other product. Roy Chowdhury et al. [11] presented an inventory model for perishable items with stock and advertisement sensitive demand. Singh and Pattnayak [12] developed a two-warehouse inventory model for deteriorating items with linear demand under conditionally permissible delay in payment. Sivashankari and Panayappan [13] studied a production-inventory model with reworking of imperfect production, scrap, and shortages.Effect of inflation and time value of money in inventory problems is well established. The initial attempt in this direction was made by Buzacott [14]. He dealt with an EOQ model under inflation subject to different types of pricing policies. Soleimani-Amiri et al. [15] compared two methods for calculating the optimal total cost in an inventory model with time value of money and inflation for deteriorating items. Hou [16] developed an inflation model for deteriorating items with stock-dependent consumption rate and completely backordered shortages by assuming a constant length of replenishment cycles and a constant fraction of the shortage length with respect to the cycle. Chern et al. [17] proposed partial backlogging inventory lot-size models for deteriorating items with fluctuating demand under inflation.Due to increasing complexities of the world economy, it is very difficult to estimate this deference precisely. Horowitz [18] discussed a simple EOQ model with a normal distribution for the inflation rate and the firm’s cost of capital. Mirzazadeh and Sarfaraz [19] presented a multiple items inventory system with budget constraint and the uniform distribution for external inflation rate. Mirzazadeh [20] proposed a detailed comparison of the average annual cost and the discounted cost models under stochastic inflationary conditions. Mirzazadeh [21] developed an inventory model under stochastic inflationary conditions with variable probability density functions (pdfs) over the time horizon. He also assumed that the demand rate is dependent on the inflation rates. Ameli et al. [22] developed an inventory model to determine ordering policy for imperfect item with fuzzy defective percentage under fuzzy discounting and inflationary condition. Other considerations of inflation with nonconstant rate are Jana et al. [23], Gholami-Qadikolaei et al. [24], Neetu and Tomer [25], Mirzazadeh et al. [26], and so forth.The demand has a considerable influence in the decisions relating to the inventory management and production planning. Although the various formations of consumption tendency have been studied, such as constant demand, price-dependent demand (e.g., [27–29]), time-dependent demand (e.g., [30–32]), and time-and-price-dependent demand (e.g., [33, 34]), Datta and Pal [35] investigated a finite time-horizon inventory model with linear time-dependent demand rate when shortages are allowed. Here we use time-dependent rate for demand with a decreasing pattern.Wu [36] investigated an inventory model with ramp type demand rate, Weibull distributed deterioration rate, and partial backlogging. Giri et al. [37] extended the ramp type demand inventory model with a more generalized Weibull deterioration.In the literature of production planning and inventory control, most papers consider inventory systems where lead times are supposed to be equal to zero or constant. In reality, lead times are rarely constant, because unpredictable events can cause random delays. For different reasons such as machine breakdowns, transport delays, and quality problems, time of component delivery from an external supplier or processing time for the semifinished product often has an uncertain duration. Frequently, lead time fluctuations strongly degrade system performance, on the other hand, especially when we consider the effect of inflation in our inventory model. In most of the inventory models under inflation, lead times are considered to be equal to zero or constant, but we consider lead time as a fuzzy function of order quantity.In this paper we extend the inventory models for deterioration item considering the following situations.(i)
Nonconstant deterioration rate is as follows: for the first case we consider fuzzy deterioration rate and for the second case we assume that the deterioration rate is time dependent and follows Weibull distribution with three known parameters. In most of the previous inventory models with Weibull deterioration rate, the location parameter (λ) of Weibull distribution is considered to be equal to zero but we consider this parameter in our model. This means that the goods start deteriorating after elapsing a specified time (λ).
(ii)
The system is under the control of fuzzy inflation.
(iii)
We introduce the fuzzy function of order quantity for lead time.
(iv)
We consider a demand rate, which is a nonincreasing function of time.The paper is organized as follows. The notation and assumptions used in the models are given in Section2. The definitions used in the models are given in Section 3. The model with fuzzy deterioration rate is studied in Section 4.1 and the corresponding one with Weibull deterioration rate is studied in Section 4.2. The models analyses that are used to find optimum policy are given in Section 5. For each model the optimal policy is obtained. Numerical examples highlighting the results obtained are given in Section 6. The conclusion remarks are given in Section 7.


## 2. Assumption and Notation

The proposed model is developed under the following assumptions.(1)
The demand rate is a linear function of time:(1)
D
(
t
)
=
a
-
b
t
a
>
0
,
b
>
0
.
(2)
The replenishment rate is infinite.
(3)
Shortages are allowed and fully backlogged except in the last cycle.
(4)
The discount rate isr.
(5)
The inflation rate (i
~) is fuzzy in nature; therefore, the net discount rate of inflation, (R
~), is also a fuzzy number, R
~
=
i
~
-
r. The present worth of p
t is p
t
e
-
R
~
t, t
≥
0 where p
t is the value of p at time t.
(6)
Deterioration of units occurs after the specified unit of time (λ) from the order reception moment and there is no repair or replacement of deteriorated units over the period [
0
,
T
].
(7)
Lead time is a function of the total ordered quantity in a period.
(8)
L
j is the lead time for periodj:
(2)
l
j
~
=
a
~
*
w
j
,
where w
j is the total ordered quantity for period j and a
~ represents the lead time for per-unit item and it is a triangular fuzzy number (
a
1
,
a
2
,
a
3
).
(9)
θ is the on-hand inventory deterioration rate at time t and may be a fuzzy number (θ
~) or a function of time and follows a three-parameter Weibull distribution function α
β
(
t
-
λ
)
β
-
1.In addition, the following are other notations that have been used in this paper:H:
finite horizon time (time units),
n:
number of replenishment cycles during the time horizonH,
T:
length of each cycle (T
=
H
/
n),
I
i
j
(
t
):
inventory level at any timet in periodj, i
=
1,2
,
3, j
=
1,2
,
…
,
n,
C
1:
holding cost per unit per unit time at time zero ($/unit/unit time),
C
2:
shortage cost per unit per unit time at time zero ($/unit/unit time),
C
3:
purchasing cost per unit item at time zero ($),
C
4:
ordering cost per order at time zero ($/set-up).Additional notations will be introduced later.


## 3. Definition

Definition 1.
IfX is a collection of objects denoted generically by x, then a fuzzy set P
~ on X is a set of order pairs.
Ifx is continuous on X, then the fuzzy set is continuous and the range of x is called the support of fuzzy set P. Continuous fuzzy sets can have various membership function shapes such as triangular and trapezoidal. In this paper, we use triangular fuzzy numbers, special case of trapezoidal, due to it straightforward structure and computational simplicity.
IfA
~
=
(
a
1
,
a
2
,
a
3
) is a triangular fuzzy number (TFN), then μ
A
~
(
x
) is defined as follows: (a
1, a
2, and a
3 are real numbers)
(3)
μ
A
~
(
x
)
=
{
x
-
a
1
a
2
-
a
1
for
a
1
≤
x
≤
a
2
a
3
-
x
a
3
-
a
2
for
a
2
≤
x
≤
a
3
0
otherwise
.Definition 2 (α-cut of fuzzy number).
Theα-cut of a fuzzy number is a crisp set which is defined as A
~
α
=
{
x
∈
R
:
μ
A
~
(
x
)
≥
α
}. In other words the α-cut is a nonempty bounded closed interval, which is denoted by A
~
α
=
[
A
L
(
α
)
α
,
A
R
(
α
)
α
].
HereA
L
(
α
)
α and A
R
(
α
)
α are the lower and upper bounds of the closed interval and
(4)
A
L
(
α
)
α
=
inf
⁡
⁡
{
x
∈
R
:
μ
A
~
(
x
)
≥
α
}
,
A
R
(
α
)
α
=
sup
⁡
⁡
{
x
∈
R
:
μ
A
~
(
x
)
≥
α
}
.Definition 3 (signed distance).
For anya
∈
R we define d
0
(
a
,
0
)
=
a. If a
>
0 the distance from a to 0 is d
0
(
a
,
0
)
=
a. If a
<
0, the distance from a to 0 is -
d
0
(
a
,
0
)
=
-
a. This is the reason why we call d
0
(
a
,
0
) the signed distance from a to 0.
The signed distances ofA
L
(
α
) and A
R
(
α
) to 0 are d
0
(
A
L
(
α
)
,
0
)
=
A
L
(
α
) and d
0
(
A
R
(
α
)
,
0
)
=
A
R
(
α
), respectively. We define the signed distance from 0 to the interval [
A
l
(
α
)
,
A
R
(
α
)
] as
(5)
d
0
(
[
A
l
(
α
)
,
A
R
(
α
)
]
,
0
)
=
1
2
(
d
0
(
A
L
(
α
)
,
0
)
,
d
0
(
A
R
(
α
)
,
0
)
)
,
1
2
(
A
L
(
α
)
+
A
R
(
α
)
)
=
1
2
[
a
+
c
+
(
2
b
-
a
-
c
)
α
]
.
Since crisp interval [
A
L
(
α
)
,
A
R
(
α
)
] has a one-to-one correspondence with level α fuzzy interval [
A
L
(
α
)
α
,
A
R
(
α
)
α
], it is natural to define the signed distance from a fuzzy interval [
A
L
(
α
)
α
,
A
R
(
α
)
α
] to 0 as
(6)
d
(
[
A
L
(
α
)
α
,
A
R
(
α
)
α
]
,
0
)
=
d
0
(
[
A
l
(
α
)
,
A
R
(
α
)
]
,
0
)
=
1
2
[
a
+
c
+
(
2
b
-
a
-
c
)
α
]
.
This is a continuous function of α on 0
≤
α
≤
1. We can find its average value through integration. Since A
~
=
⋃
0
≤
α
≤
1
[
A
L
(
α
)
α
,
A
R
(
α
)
α
], we have the following definition.Definition 4.
IfA
~
=
(
a
,
b
,
c
) is a triangular fuzzy number, the signed distance from A
~ to 0
~ is defined as
(7)
d
(
A
~
,
0
~
)
=
∫
0
1
d
(
[
A
L
(
α
)
α
,
A
R
(
α
)
α
]
,
0
)
d
α
=
1
4
(
2
b
+
c
+
a
)
.Property 1.
By applying the definition of anα-level cut to we have the following exact arithmetic operations for continuous fuzzy set:
(8)
M
~
+
N
~
=
[
m
1
α
,
m
2
α
]
+
[
n
1
α
,
n
2
α
]
=
[
m
1
α
+
n
1
α
,
m
2
α
+
n
2
α
]
,
M
~
-
N
~
=
[
m
1
α
,
m
2
α
]
-
[
n
1
α
,
n
2
α
]
=
[
m
1
α
-
n
1
α
,
m
2
α
-
n
2
α
]
,
M
~
*
N
~
=
[
m
1
α
,
m
2
α
]
*
[
n
1
α
,
n
2
α
]
=
[
m
1
α
*
n
1
α
,
m
2
α
*
n
2
α
]
,
M
~
÷
N
~
=
[
m
1
α
,
m
2
α
]
÷
[
n
1
α
,
n
2
α
]
=
[
m
1
α
÷
n
2
α
,
m
2
α
÷
n
1
α
]
.


## 4. Model Formulation

In the development of the model, we assume that there areN cycles during the real time horizon H. The time horizon, H, is divided into n equal cycles each of length T so that T
=
H
/
n.


### 4.1. Formulation forjth (0
≤
j
≤
(
n
-
1
)) Cycle

For Case(
1
), the deterioration rate is a fuzzy number.The graphical representation of the inventory system is shown in Figure1.Figure 1
Graphical representation of the inventory system 1.Each inventory cycle except the last cycle can be divided into three parts. Company orders an amountI
1
j
(
0
) of the item at time zero. Then, the level of inventory is decreasing by only the consumption rates. At the moment of λ, the item of inventory starts deteriorating. At the moment of k
T, 0
≤
k
≤
1 inventory level leads to zero and shortages occur. And at the moment (K
T
j
-
l
j
~) of every jth (j
=
1,2
,
…
,
n) cycle, company purchases an amount Q
j units of the item. During the time interval [
k
T
j
,
T
j
] we do not have any deterioration, and therefore the shortages level linearly increases by the demand rate. Shortages are not allowed for the last replenishment cycle.Because of inflation, the purchasing cost has effect on the optimal order quantity; hence, this cost has also to be included in the analysis. At the first part ofjth cycle there is no deterioration (λ unit of time after the replenishment time). The inventory level, I
1
j
(
t
), satisfies the following differential equations:
(9)
∂
I
1
j
(
t
)
∂
t
=
-
a
+
b
t
T
(
j
-
1
)
≤
t
≤
T
(
j
-
1
)
+
λ
.
The second part of the inventory level can be described as follows:
(10)
∂
I
2
j
(
t
)
∂
t
=
-
a
+
b
t
-
θ
~
I
2
j
(
t
)
T
(
j
-
1
)
+
λ
≤
t
≤
k
T
j
.
We can solve the above differential equations by using the following equations as a boundary condition:
(11)
I
2
j
(
T
(
j
-
1
)
+
k
T
)
=
0
.
We can see
(12)
I
2
j
(
t
)
=
-
b
θ
2
-
a
-
b
t
θ
+
e
θ
(
S
-
t
)
(
a
-
b
S
θ
+
b
θ
2
)
,
I
2
j
(
T
(
j
-
1
)
+
λ
)
=
b
θ
2
-
a
+
b
[
T
(
j
-
1
)
+
λ
]
θ
+
e
θ
[
T
(
j
(
k
-
1
)
+
1
)
-
λ
]
+
a
(
T
(
j
-
1
)
+
λ
)
+
1
2
b
(
T
(
j
-
1
)
+
λ
)
2
,
I
2
j
(
T
(
j
-
1
)
+
λ
)
=
I
1
j
(
T
(
j
-
1
)
+
λ
)
,
I
1
j
(
t
)
=
b
θ
2
-
a
+
b
[
T
(
j
-
1
)
+
λ
]
θ
+
e
θ
[
T
(
j
(
k
-
1
)
+
1
)
-
λ
]
+
a
(
T
(
j
-
1
)
+
λ
)
+
1
2
b
(
T
(
j
-
1
)
+
λ
)
2
-
(
a
t
+
b
t
2
2
)
.
For the last part of jth cycle, the shortages level is governed by
(13)
∂
I
3
j
(
t
)
∂
t
=
-
(
a
+
b
t
)
,
k
T
j
≤
t
≤
j
T
,
I
3
(
T
k
j
)
=
0
,
I
3
j
(
t
)
=
a
k
T
j
+
b
k
2
T
2
j
2
-
(
a
t
+
0.5
b
t
2
)
,
k
T
j
≤
t
≤
j
T
,
I
31
(
0
)
=
0
.
The objective of the problem is minimization of the total present value of costs over the time horizon. Consider (PVCP) as the present value of purchasing cost, (PVCH) as the present value of holding cost, (PVCS) as the present value of costs of shortages, and (PVCO) as the expected present value of cost of ordering. The total present value of costs over time horizon (TPV) is
(14)
TPV
=
PVCP
+
PVCH
+
PVCS
+
PVCO
.


#### 4.1.1. Present Value of Holding Cost

The inventory is carried out over the first and second period of each cycle so the present value of holding cost forjth cycle is given by
(15)
PVCH
j
=
C
1
(
∫
T
(
j
-
1
)
T
(
j
-
1
)
+
λ
I
1
j
(
t
)
e
-
R
~
t
d
t
s
s
s
s
s
s
s
+
∫
T
(
j
-
1
)
+
λ
T
(
j
-
1
)
+
k
T
I
2
j
(
t
)
e
-
R
~
t
d
t
)
=
C
1
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
+
1
R
~
3
[
b
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
s
s
s
s
s
s
s
s
s
s
s
×
(
-
R
~
(
T
(
j
-
1
)
+
λ
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
+
R
~
(
T
(
j
-
1
)
+
λ
)
)
-
2
)
s
s
s
s
s
s
s
s
s
s
s
+
e
-
R
~
λ
(
(
R
~
(
T
(
j
-
1
)
)
×
(
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
]
+
e
-
R
~
(
T
(
j
-
1
)
+
λ
+
K
T
j
)
R
~
2
*
[
e
-
R
~
K
T
j
(
b
+
a
R
~
+
b
R
~
(
T
(
j
-
1
)
+
λ
)
)
s
s
s
s
e
-
R
~
K
T
j
+
(
b
+
a
R
~
+
b
R
~
K
T
j
)
]
+
b
R
~
θ
2
[
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
-
e
-
R
~
K
T
j
]
+
e
θ
K
T
j
[
a
+
b
K
T
j
θ
-
b
θ
2
]
*
(
e
K
T
j
(
R
~
+
θ
)
-
e
(
T
(
j
-
1
)
+
λ
)
(
R
~
+
θ
)
R
~
+
θ
)
.
So the total present value of holding cost is
(16)
PVCH
=
∑
j
=
1
n
PVCH
j
.


#### 4.1.2. Present Value of Shortage Cost

Shortage occurs during the third period of each cycle. The present value of the shortage cost for one cycle can be written as(17)
PVCS
j
=
C
2
(
∫
K
T
j
T
j
I
3
j
(
t
)
e
-
R
~
t
d
t
)
=
a
K
T
j
+
b
k
2
T
2
j
2
R
~
(
e
-
R
~
K
T
j
-
e
-
R
~
T
j
)
-
a
e
-
R
~
T
j
(
1
+
k
)
2
R
2
(
e
-
R
~
T
j
(
1
+
R
~
T
j
(
1
+
k
)
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
-
e
-
R
~
T
j
(
1
+
k
)
(
1
+
R
~
T
j
)
)
+
-
b
e
-
R
~
T
j
(
1
+
k
)
2
R
~
3
×
[
e
-
R
~
T
j
(
2
+
R
~
K
T
j
(
2
+
R
~
K
T
j
)
)
(
2
+
R
~
T
j
(
1
+
k
)
(
2
+
T
j
(
1
+
k
)
)
)
s
s
s
s
s
s
-
e
-
R
~
T
j
(
1
+
k
)
(
2
+
R
~
T
j
(
1
+
k
)
(
2
+
T
j
(
1
+
k
)
)
)
]
.
So the total present value of shortage cost is
(18)
PVCS
=
∑
j
=
1
n
-
1
PVCS
j
.


#### 4.1.3. Present Value of Purchasing Cost

The ordering quantity consists of both demand and deterioration for the relevant period excluding shortage part of the period and the amount required to satisfy the demand during the shortage period in the preceding time interval. For thejth cycle (j
=
1
,
2
,
…
,
n
-
1) the present value of purchasing cost can be formulated as follows:
(19)
PVCP
j
=
C
3
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
e
-
R
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
)
sssssss
+
I
1
j
(
T
(
j
-
1
)
)
e
-
R
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
)
)
.
Therefore, the total present value of purchasing cost is
(20)
CP
=
∑
j
=
1
m
C
3
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
s
s
s
s
s
s
s
s
s
s
+
I
1
j
(
T
(
j
-
1
)
)
)
e
-
R
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
)
.


#### 4.1.4. Present Value of Ordering Cost

Present value of sales revenue for thejth (1
≤
j
≤
N) cycle is given by
(21)
CA
=
∑
j
=
1
m
CA
j
=
C
4
∑
j
=
1
m
e
-
R
~
(
T
j
-
l
~
(
I
3
j
(
T
j
)
+
I
2
j
(
T
j
)
)
)
.


### 4.2. Formulation forjth (0
≤
j
≤
(
n
-
1
)) Cycle

For Case(
2
), time-dependent (Weibull) deterioration rate, we have the following.The graphical representation of the inventory system is shown in Figure2.Figure 2
Graphical representation of the inventory system 2.The methodology of this case is like the previous case except in this case deterioration rate follows three-parameter Weibull distributions.Three-parameter density function for shelf life product is defined as(22)
f
(
x
)
=
α
β
(
t
-
λ
)
β
-
1
e
-
α
(
t
-
λ
)
β
,
t
≥
λ
,
where α is the scale parameter (α
>
0), β is the shape parameter (β
>
0), and λ is the location parameter so the cumulative distribution function of f
(
t
) is F
(
t
)
=
e
-
α
t
β. Therefore, R
(
t
)
=
1
-
F
(
t
) and finally Weibull instantaneous rate function for the stocked item is
(23)
θ
(
t
)
=
f
(
t
)
1
-
F
(
t
)
=
α
β
(
t
-
λ
)
β
-
1
,
t
≥
λ
.
In other words when the time of deterioration of an item follows this density function, this item starts deteriorating after λ unit time after it is received to the inventory.For this case the inventory level at the first part ofjth cycle satisfies the following differential equation:
(24)
∂
I
1
j
(
t
)
∂
t
=
-
a
+
b
t
,
T
(
j
-
1
)
≤
t
≤
T
(
j
-
1
)
+
λ
,
j
=
1,2
,
…
,
m
.
In this case, the inventory levels at the second part of jth cycle satisfy the following differential equation and the equation becomes
(25)
∂
I
2
j
(
t
)
∂
t
=
-
a
+
b
t
-
α
β
(
t
-
λ
)
β
-
1
I
2
(
t
)
,
T
(
j
-
1
)
+
λ
≤
t
≤
T
(
j
-
1
)
K
,
S
=
T
(
j
-
1
)
+
K
T
,
I
2
j
(
S
)
=
0
,
I
2
j
(
t
)
=
e
-
α
(
t
-
λ
)
β
∫
t
S
(
a
-
b
u
)
e
α
(
u
-
λ
)
β
d
u
,
T
(
j
-
1
)
+
λ
≤
t
≤
T
(
j
-
1
)
+
K
T
.
So for the I
2
j
(
λ
) we have
(26)
I
2
j
(
λ
)
=
∫
λ
S
(
a
-
b
u
)
(
∑
n
=
0
∞
(
α
(
u
-
λ
)
β
)
n
n
!
)
≅
a
(
s
-
λ
)
+
α
a
(
s
-
λ
)
β
+
1
β
+
1
-
b
[
(
s
-
λ
)
2
2
+
λ
(
s
-
λ
)
s
s
s
s
s
s
s
s
+
α
(
(
s
-
λ
)
β
+
2
β
+
2
+
λ
(
s
-
λ
)
β
+
1
β
+
1
)
]
=
Q
j
,
I
1
j
(
λ
)
=
I
2
j
(
λ
)
=
Q
j
,
I
1
j
(
t
)
=
b
t
2
2
-
a
t
+
Q
j
-
b
λ
2
2
+
a
λ
,
T
(
j
-
1
)
≤
t
≤
T
(
j
-
1
)
+
λ
.
For the last part of jth cycle, the shortages level formulation is similar to the first case. The detailed analysis is given as follows.


#### 4.2.1. Present Value of Holding Cost

The present value of holding cost forjth cycle is given by
(27)
C
H
j
=
C
1
(
∫
T
(
j
-
1
)
T
(
j
-
1
)
+
λ
I
1
j
(
t
)
e
-
R
~
t
d
t
+
∫
T
(
j
-
1
)
+
λ
K
T
j
I
2
j
(
t
)
e
-
R
~
t
d
t
)
=
C
1
(
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
sssssss
+
1
R
~
3
[
b
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
-
R
~
(
T
(
j
-
1
)
+
λ
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
+
R
~
(
T
(
j
-
1
)
+
λ
)
)
-
2
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
+
e
-
R
~
λ
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
×
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
]
ssssssss
+
∫
T
(
j
-
1
)
+
λ
S
(
e
-
α
(
t
-
λ
)
β
∫
t
S
(
a
-
b
u
)
sssssssssssssssssssssss
×
∫
t
S
(
a
-
b
u
)
e
α
(
u
-
λ
)
β
d
u
)
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
e
-
R
~
t
d
t
)
.
For the simplified formulation we define
(28)
C
=
a
s
+
a
α
(
s
-
λ
)
β
+
1
β
+
1
-
b
2
S
2
-
b
α
(
s
-
λ
)
β
+
1
(
s
b
+
s
-
λ
)
(
β
+
2
)
(
β
+
1
)
.
Therefore we have
(29)
C
H
j
≅
C
1
(
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
+
1
R
~
3
s
s
s
s
s
s
s
×
[
b
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
s
s
s
s
s
s
s
s
s
s
×
(
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
2
+
R
~
(
T
(
j
-
1
)
+
λ
)
)
-
2
)
s
s
s
s
s
s
s
s
s
s
+
e
-
R
~
λ
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
]
s
s
s
s
s
s
s
+
C
(
K
T
-
λ
)
-
(
R
C
+
a
)
2
[
S
2
-
T
(
j
-
1
)
+
λ
2
]
s
s
s
s
s
s
s
-
a
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
-
(
T
(
j
-
1
)
)
β
+
2
]
(
β
+
1
)
(
β
+
2
)
ssssssss
+
(
b
+
2
R
a
)
[
S
3
-
(
T
(
j
-
1
)
+
λ
)
3
]
6
s
s
s
s
s
s
s
s
+
(
b
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
S
(
2
+
3
β
+
β
2
)
-
2
λ
-
T
(
j
-
1
)
β
+
2
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
(
T
(
j
-
1
)
+
λ
)
sssssssssssssssssssssss
b
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
×
(
2
+
3
β
+
β
2
)
-
2
λ
)
)
(
b
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
]
)
s
s
s
s
s
s
s
s
×
(
(
β
+
1
)
(
β
+
3
)
(
β
+
2
)
2
)
-
1
s
s
s
s
s
s
s
s
+
(
a
2
[
(
T
(
j
-
1
)
+
λ
)
β
+
1
(
(
T
(
j
-
1
)
+
λ
)
(
β
)
+
2
λ
+
T
(
j
-
1
)
)
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
S
β
+
λ
+
K
T
j
)
-
T
(
j
-
1
)
β
+
1
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
a
2
[
(
T
(
j
-
1
)
+
λ
)
β
+
1
(
(
T
(
j
-
1
)
+
λ
)
(
β
)
+
2
λ
+
T
(
j
-
1
)
)
)
×
(
(
T
(
j
-
1
)
+
λ
)
(
β
)
+
2
λ
+
T
(
j
-
1
)
)
]
)
s
s
s
s
s
s
s
s
×
(
(
β
+
1
)
(
β
+
2
)
)
-
1
s
s
s
s
s
s
s
s
+
α
a
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
2
(
β
+
1
)
2
s
s
s
s
s
s
s
s
-
b
a
2
(
β
+
1
)
(
β
+
3
)
(
β
+
2
)
s
s
s
s
s
s
s
s
×
[
(
s
-
λ
)
β
+
1
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
λ
S
(
β
+
1
)
+
K
2
T
2
j
2
β
(
β
+
3
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
+
2
λ
2
+
2
K
2
T
2
j
2
)
s
s
s
s
s
s
s
s
s
s
s
s
-
(
T
(
j
-
1
)
)
β
+
1
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
λ
(
T
(
j
-
1
)
+
λ
)
(
β
+
1
)
+
F
2
β
(
β
+
3
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
(
s
-
λ
)
β
+
1
+
2
λ
2
+
2
(
T
(
j
-
1
)
+
λ
)
2
)
]
s
s
s
s
s
s
s
s
-
a
c
[
(
S
-
λ
)
β
+
1
-
(
T
(
j
-
1
)
)
β
+
1
]
β
+
1
s
s
s
s
s
s
s
s
+
(
R
a
α
[
(
S
-
λ
)
β
+
2
(
S
(
β
+
2
)
+
λ
)
-
(
T
(
j
-
1
)
)
β
+
2
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
×
(
(
T
(
j
-
1
)
+
λ
)
(
β
+
2
)
+
λ
)
]
)
s
s
s
s
s
s
s
s
×
(
(
β
+
1
)
(
β
+
3
)
(
β
+
2
)
)
-
1
s
s
s
s
s
s
s
s
-
R
b
[
K
4
T
4
j
4
-
(
T
(
j
-
1
)
+
λ
)
4
]
8
)
.Present Value of Shortage Cost. The present value of shortage cost is similar to the previous case.


#### 4.2.2. The Present Value of Purchasing Cost

The total amount that we have to purchase forjth cycle is
(30)
I
1
(
j
)
(
T
j
)
+
I
3
(
j
-
1
)
(
T
j
)
.
The lead time for this amount of item is l
j
=
l
~
w
j.So the total present value of purchasing cost is(31)
PVCP
=
C
3
I
11
(
0
)
+
C
3
∑
j
=
1
n
-
1
(
I
1
(
j
)
(
T
j
)
-
I
3
(
j
-
1
)
(
T
j
)
)
×
e
-
R
~
(
T
j
-
l
~
(
(
I
1
(
j
)
(
T
j
)
-
I
3
(
j
-
1
)
(
T
j
)
)
)
.


#### 4.2.3. The Present Value of Ordering Cost

For thejth cycle we have to order at T
(
j
-
1
)
-
l
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
).So the total present value of ordering cost for thejth cycle is
(32)
PVCO
=
∑
j
=
1
m
C
A
j
=
C
4
∑
j
=
1
m
e
-
R
~
(
T
j
-
l
~
(
I
3
j
(
T
j
)
+
I
2
j
(
T
j
)
)
)
.


#### 4.2.4. Present Value of Deteriorated Item Cost

The amount of deteriorated items during[
T
(
j
-
1
)
+
λ
,
T
(
j
-
1
)
+
k
T
] is
(33)
I
2
j
(
λ
)
-
∫
λ
S
D
(
t
)
d
t
=
a
(
s
-
λ
)
+
α
a
(
s
-
λ
)
β
+
1
β
+
1
-
b
[
(
s
-
λ
)
2
2
+
λ
(
s
-
λ
)
s
s
s
s
s
s
s
s
+
α
(
(
s
-
λ
)
β
+
2
β
+
2
+
λ
(
s
-
λ
)
β
+
1
β
+
1
)
]
-
a
(
K
T
j
-
λ
)
+
b
(
k
2
T
2
j
2
-
λ
2
)
2
.
The total present value of the amount of deteriorated items cost is
(34)
PVCD
=
∑
j
=
1
m
C
D
j
=
θ
~
C
5
∑
j
=
1
m
∫
T
(
j
-
1
)
+
λ
K
T
j
(
e
-
α
(
t
-
λ
)
β
∫
t
S
(
a
-
b
u
)
e
α
(
u
-
λ
)
β
d
u
)
e
-
R
~
t
d
t
.
The total present value of costs over time horizon in this case (TPV) is
(35)
TPV
=
PVCP
+
PVCH
+
PVCS
+
PVCO
+
PVCD
.


## 4.1. Formulation forjth (0
≤
j
≤
(
n
-
1
)) Cycle

For Case(
1
), the deterioration rate is a fuzzy number.The graphical representation of the inventory system is shown in Figure1.Figure 1
Graphical representation of the inventory system 1.Each inventory cycle except the last cycle can be divided into three parts. Company orders an amountI
1
j
(
0
) of the item at time zero. Then, the level of inventory is decreasing by only the consumption rates. At the moment of λ, the item of inventory starts deteriorating. At the moment of k
T, 0
≤
k
≤
1 inventory level leads to zero and shortages occur. And at the moment (K
T
j
-
l
j
~) of every jth (j
=
1,2
,
…
,
n) cycle, company purchases an amount Q
j units of the item. During the time interval [
k
T
j
,
T
j
] we do not have any deterioration, and therefore the shortages level linearly increases by the demand rate. Shortages are not allowed for the last replenishment cycle.Because of inflation, the purchasing cost has effect on the optimal order quantity; hence, this cost has also to be included in the analysis. At the first part ofjth cycle there is no deterioration (λ unit of time after the replenishment time). The inventory level, I
1
j
(
t
), satisfies the following differential equations:
(9)
∂
I
1
j
(
t
)
∂
t
=
-
a
+
b
t
T
(
j
-
1
)
≤
t
≤
T
(
j
-
1
)
+
λ
.
The second part of the inventory level can be described as follows:
(10)
∂
I
2
j
(
t
)
∂
t
=
-
a
+
b
t
-
θ
~
I
2
j
(
t
)
T
(
j
-
1
)
+
λ
≤
t
≤
k
T
j
.
We can solve the above differential equations by using the following equations as a boundary condition:
(11)
I
2
j
(
T
(
j
-
1
)
+
k
T
)
=
0
.
We can see
(12)
I
2
j
(
t
)
=
-
b
θ
2
-
a
-
b
t
θ
+
e
θ
(
S
-
t
)
(
a
-
b
S
θ
+
b
θ
2
)
,
I
2
j
(
T
(
j
-
1
)
+
λ
)
=
b
θ
2
-
a
+
b
[
T
(
j
-
1
)
+
λ
]
θ
+
e
θ
[
T
(
j
(
k
-
1
)
+
1
)
-
λ
]
+
a
(
T
(
j
-
1
)
+
λ
)
+
1
2
b
(
T
(
j
-
1
)
+
λ
)
2
,
I
2
j
(
T
(
j
-
1
)
+
λ
)
=
I
1
j
(
T
(
j
-
1
)
+
λ
)
,
I
1
j
(
t
)
=
b
θ
2
-
a
+
b
[
T
(
j
-
1
)
+
λ
]
θ
+
e
θ
[
T
(
j
(
k
-
1
)
+
1
)
-
λ
]
+
a
(
T
(
j
-
1
)
+
λ
)
+
1
2
b
(
T
(
j
-
1
)
+
λ
)
2
-
(
a
t
+
b
t
2
2
)
.
For the last part of jth cycle, the shortages level is governed by
(13)
∂
I
3
j
(
t
)
∂
t
=
-
(
a
+
b
t
)
,
k
T
j
≤
t
≤
j
T
,
I
3
(
T
k
j
)
=
0
,
I
3
j
(
t
)
=
a
k
T
j
+
b
k
2
T
2
j
2
-
(
a
t
+
0.5
b
t
2
)
,
k
T
j
≤
t
≤
j
T
,
I
31
(
0
)
=
0
.
The objective of the problem is minimization of the total present value of costs over the time horizon. Consider (PVCP) as the present value of purchasing cost, (PVCH) as the present value of holding cost, (PVCS) as the present value of costs of shortages, and (PVCO) as the expected present value of cost of ordering. The total present value of costs over time horizon (TPV) is
(14)
TPV
=
PVCP
+
PVCH
+
PVCS
+
PVCO
.


### 4.1.1. Present Value of Holding Cost

The inventory is carried out over the first and second period of each cycle so the present value of holding cost forjth cycle is given by
(15)
PVCH
j
=
C
1
(
∫
T
(
j
-
1
)
T
(
j
-
1
)
+
λ
I
1
j
(
t
)
e
-
R
~
t
d
t
s
s
s
s
s
s
s
+
∫
T
(
j
-
1
)
+
λ
T
(
j
-
1
)
+
k
T
I
2
j
(
t
)
e
-
R
~
t
d
t
)
=
C
1
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
+
1
R
~
3
[
b
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
s
s
s
s
s
s
s
s
s
s
s
×
(
-
R
~
(
T
(
j
-
1
)
+
λ
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
+
R
~
(
T
(
j
-
1
)
+
λ
)
)
-
2
)
s
s
s
s
s
s
s
s
s
s
s
+
e
-
R
~
λ
(
(
R
~
(
T
(
j
-
1
)
)
×
(
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
]
+
e
-
R
~
(
T
(
j
-
1
)
+
λ
+
K
T
j
)
R
~
2
*
[
e
-
R
~
K
T
j
(
b
+
a
R
~
+
b
R
~
(
T
(
j
-
1
)
+
λ
)
)
s
s
s
s
e
-
R
~
K
T
j
+
(
b
+
a
R
~
+
b
R
~
K
T
j
)
]
+
b
R
~
θ
2
[
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
-
e
-
R
~
K
T
j
]
+
e
θ
K
T
j
[
a
+
b
K
T
j
θ
-
b
θ
2
]
*
(
e
K
T
j
(
R
~
+
θ
)
-
e
(
T
(
j
-
1
)
+
λ
)
(
R
~
+
θ
)
R
~
+
θ
)
.
So the total present value of holding cost is
(16)
PVCH
=
∑
j
=
1
n
PVCH
j
.


### 4.1.2. Present Value of Shortage Cost

Shortage occurs during the third period of each cycle. The present value of the shortage cost for one cycle can be written as(17)
PVCS
j
=
C
2
(
∫
K
T
j
T
j
I
3
j
(
t
)
e
-
R
~
t
d
t
)
=
a
K
T
j
+
b
k
2
T
2
j
2
R
~
(
e
-
R
~
K
T
j
-
e
-
R
~
T
j
)
-
a
e
-
R
~
T
j
(
1
+
k
)
2
R
2
(
e
-
R
~
T
j
(
1
+
R
~
T
j
(
1
+
k
)
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
-
e
-
R
~
T
j
(
1
+
k
)
(
1
+
R
~
T
j
)
)
+
-
b
e
-
R
~
T
j
(
1
+
k
)
2
R
~
3
×
[
e
-
R
~
T
j
(
2
+
R
~
K
T
j
(
2
+
R
~
K
T
j
)
)
(
2
+
R
~
T
j
(
1
+
k
)
(
2
+
T
j
(
1
+
k
)
)
)
s
s
s
s
s
s
-
e
-
R
~
T
j
(
1
+
k
)
(
2
+
R
~
T
j
(
1
+
k
)
(
2
+
T
j
(
1
+
k
)
)
)
]
.
So the total present value of shortage cost is
(18)
PVCS
=
∑
j
=
1
n
-
1
PVCS
j
.


### 4.1.3. Present Value of Purchasing Cost

The ordering quantity consists of both demand and deterioration for the relevant period excluding shortage part of the period and the amount required to satisfy the demand during the shortage period in the preceding time interval. For thejth cycle (j
=
1
,
2
,
…
,
n
-
1) the present value of purchasing cost can be formulated as follows:
(19)
PVCP
j
=
C
3
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
e
-
R
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
)
sssssss
+
I
1
j
(
T
(
j
-
1
)
)
e
-
R
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
)
)
.
Therefore, the total present value of purchasing cost is
(20)
CP
=
∑
j
=
1
m
C
3
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
s
s
s
s
s
s
s
s
s
s
+
I
1
j
(
T
(
j
-
1
)
)
)
e
-
R
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
)
.


### 4.1.4. Present Value of Ordering Cost

Present value of sales revenue for thejth (1
≤
j
≤
N) cycle is given by
(21)
CA
=
∑
j
=
1
m
CA
j
=
C
4
∑
j
=
1
m
e
-
R
~
(
T
j
-
l
~
(
I
3
j
(
T
j
)
+
I
2
j
(
T
j
)
)
)
.


## 4.1.1. Present Value of Holding Cost

The inventory is carried out over the first and second period of each cycle so the present value of holding cost forjth cycle is given by
(15)
PVCH
j
=
C
1
(
∫
T
(
j
-
1
)
T
(
j
-
1
)
+
λ
I
1
j
(
t
)
e
-
R
~
t
d
t
s
s
s
s
s
s
s
+
∫
T
(
j
-
1
)
+
λ
T
(
j
-
1
)
+
k
T
I
2
j
(
t
)
e
-
R
~
t
d
t
)
=
C
1
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
+
1
R
~
3
[
b
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
s
s
s
s
s
s
s
s
s
s
s
×
(
-
R
~
(
T
(
j
-
1
)
+
λ
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
+
R
~
(
T
(
j
-
1
)
+
λ
)
)
-
2
)
s
s
s
s
s
s
s
s
s
s
s
+
e
-
R
~
λ
(
(
R
~
(
T
(
j
-
1
)
)
×
(
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
]
+
e
-
R
~
(
T
(
j
-
1
)
+
λ
+
K
T
j
)
R
~
2
*
[
e
-
R
~
K
T
j
(
b
+
a
R
~
+
b
R
~
(
T
(
j
-
1
)
+
λ
)
)
s
s
s
s
e
-
R
~
K
T
j
+
(
b
+
a
R
~
+
b
R
~
K
T
j
)
]
+
b
R
~
θ
2
[
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
-
e
-
R
~
K
T
j
]
+
e
θ
K
T
j
[
a
+
b
K
T
j
θ
-
b
θ
2
]
*
(
e
K
T
j
(
R
~
+
θ
)
-
e
(
T
(
j
-
1
)
+
λ
)
(
R
~
+
θ
)
R
~
+
θ
)
.
So the total present value of holding cost is
(16)
PVCH
=
∑
j
=
1
n
PVCH
j
.


## 4.1.2. Present Value of Shortage Cost

Shortage occurs during the third period of each cycle. The present value of the shortage cost for one cycle can be written as(17)
PVCS
j
=
C
2
(
∫
K
T
j
T
j
I
3
j
(
t
)
e
-
R
~
t
d
t
)
=
a
K
T
j
+
b
k
2
T
2
j
2
R
~
(
e
-
R
~
K
T
j
-
e
-
R
~
T
j
)
-
a
e
-
R
~
T
j
(
1
+
k
)
2
R
2
(
e
-
R
~
T
j
(
1
+
R
~
T
j
(
1
+
k
)
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
-
e
-
R
~
T
j
(
1
+
k
)
(
1
+
R
~
T
j
)
)
+
-
b
e
-
R
~
T
j
(
1
+
k
)
2
R
~
3
×
[
e
-
R
~
T
j
(
2
+
R
~
K
T
j
(
2
+
R
~
K
T
j
)
)
(
2
+
R
~
T
j
(
1
+
k
)
(
2
+
T
j
(
1
+
k
)
)
)
s
s
s
s
s
s
-
e
-
R
~
T
j
(
1
+
k
)
(
2
+
R
~
T
j
(
1
+
k
)
(
2
+
T
j
(
1
+
k
)
)
)
]
.
So the total present value of shortage cost is
(18)
PVCS
=
∑
j
=
1
n
-
1
PVCS
j
.


## 4.1.3. Present Value of Purchasing Cost

The ordering quantity consists of both demand and deterioration for the relevant period excluding shortage part of the period and the amount required to satisfy the demand during the shortage period in the preceding time interval. For thejth cycle (j
=
1
,
2
,
…
,
n
-
1) the present value of purchasing cost can be formulated as follows:
(19)
PVCP
j
=
C
3
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
e
-
R
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
)
sssssss
+
I
1
j
(
T
(
j
-
1
)
)
e
-
R
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
)
)
.
Therefore, the total present value of purchasing cost is
(20)
CP
=
∑
j
=
1
m
C
3
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
s
s
s
s
s
s
s
s
s
s
+
I
1
j
(
T
(
j
-
1
)
)
)
e
-
R
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
)
.


## 4.1.4. Present Value of Ordering Cost

Present value of sales revenue for thejth (1
≤
j
≤
N) cycle is given by
(21)
CA
=
∑
j
=
1
m
CA
j
=
C
4
∑
j
=
1
m
e
-
R
~
(
T
j
-
l
~
(
I
3
j
(
T
j
)
+
I
2
j
(
T
j
)
)
)
.


## 4.2. Formulation forjth (0
≤
j
≤
(
n
-
1
)) Cycle

For Case(
2
), time-dependent (Weibull) deterioration rate, we have the following.The graphical representation of the inventory system is shown in Figure2.Figure 2
Graphical representation of the inventory system 2.The methodology of this case is like the previous case except in this case deterioration rate follows three-parameter Weibull distributions.Three-parameter density function for shelf life product is defined as(22)
f
(
x
)
=
α
β
(
t
-
λ
)
β
-
1
e
-
α
(
t
-
λ
)
β
,
t
≥
λ
,
where α is the scale parameter (α
>
0), β is the shape parameter (β
>
0), and λ is the location parameter so the cumulative distribution function of f
(
t
) is F
(
t
)
=
e
-
α
t
β. Therefore, R
(
t
)
=
1
-
F
(
t
) and finally Weibull instantaneous rate function for the stocked item is
(23)
θ
(
t
)
=
f
(
t
)
1
-
F
(
t
)
=
α
β
(
t
-
λ
)
β
-
1
,
t
≥
λ
.
In other words when the time of deterioration of an item follows this density function, this item starts deteriorating after λ unit time after it is received to the inventory.For this case the inventory level at the first part ofjth cycle satisfies the following differential equation:
(24)
∂
I
1
j
(
t
)
∂
t
=
-
a
+
b
t
,
T
(
j
-
1
)
≤
t
≤
T
(
j
-
1
)
+
λ
,
j
=
1,2
,
…
,
m
.
In this case, the inventory levels at the second part of jth cycle satisfy the following differential equation and the equation becomes
(25)
∂
I
2
j
(
t
)
∂
t
=
-
a
+
b
t
-
α
β
(
t
-
λ
)
β
-
1
I
2
(
t
)
,
T
(
j
-
1
)
+
λ
≤
t
≤
T
(
j
-
1
)
K
,
S
=
T
(
j
-
1
)
+
K
T
,
I
2
j
(
S
)
=
0
,
I
2
j
(
t
)
=
e
-
α
(
t
-
λ
)
β
∫
t
S
(
a
-
b
u
)
e
α
(
u
-
λ
)
β
d
u
,
T
(
j
-
1
)
+
λ
≤
t
≤
T
(
j
-
1
)
+
K
T
.
So for the I
2
j
(
λ
) we have
(26)
I
2
j
(
λ
)
=
∫
λ
S
(
a
-
b
u
)
(
∑
n
=
0
∞
(
α
(
u
-
λ
)
β
)
n
n
!
)
≅
a
(
s
-
λ
)
+
α
a
(
s
-
λ
)
β
+
1
β
+
1
-
b
[
(
s
-
λ
)
2
2
+
λ
(
s
-
λ
)
s
s
s
s
s
s
s
s
+
α
(
(
s
-
λ
)
β
+
2
β
+
2
+
λ
(
s
-
λ
)
β
+
1
β
+
1
)
]
=
Q
j
,
I
1
j
(
λ
)
=
I
2
j
(
λ
)
=
Q
j
,
I
1
j
(
t
)
=
b
t
2
2
-
a
t
+
Q
j
-
b
λ
2
2
+
a
λ
,
T
(
j
-
1
)
≤
t
≤
T
(
j
-
1
)
+
λ
.
For the last part of jth cycle, the shortages level formulation is similar to the first case. The detailed analysis is given as follows.


### 4.2.1. Present Value of Holding Cost

The present value of holding cost forjth cycle is given by
(27)
C
H
j
=
C
1
(
∫
T
(
j
-
1
)
T
(
j
-
1
)
+
λ
I
1
j
(
t
)
e
-
R
~
t
d
t
+
∫
T
(
j
-
1
)
+
λ
K
T
j
I
2
j
(
t
)
e
-
R
~
t
d
t
)
=
C
1
(
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
sssssss
+
1
R
~
3
[
b
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
-
R
~
(
T
(
j
-
1
)
+
λ
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
+
R
~
(
T
(
j
-
1
)
+
λ
)
)
-
2
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
+
e
-
R
~
λ
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
×
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
]
ssssssss
+
∫
T
(
j
-
1
)
+
λ
S
(
e
-
α
(
t
-
λ
)
β
∫
t
S
(
a
-
b
u
)
sssssssssssssssssssssss
×
∫
t
S
(
a
-
b
u
)
e
α
(
u
-
λ
)
β
d
u
)
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
e
-
R
~
t
d
t
)
.
For the simplified formulation we define
(28)
C
=
a
s
+
a
α
(
s
-
λ
)
β
+
1
β
+
1
-
b
2
S
2
-
b
α
(
s
-
λ
)
β
+
1
(
s
b
+
s
-
λ
)
(
β
+
2
)
(
β
+
1
)
.
Therefore we have
(29)
C
H
j
≅
C
1
(
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
+
1
R
~
3
s
s
s
s
s
s
s
×
[
b
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
s
s
s
s
s
s
s
s
s
s
×
(
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
2
+
R
~
(
T
(
j
-
1
)
+
λ
)
)
-
2
)
s
s
s
s
s
s
s
s
s
s
+
e
-
R
~
λ
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
]
s
s
s
s
s
s
s
+
C
(
K
T
-
λ
)
-
(
R
C
+
a
)
2
[
S
2
-
T
(
j
-
1
)
+
λ
2
]
s
s
s
s
s
s
s
-
a
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
-
(
T
(
j
-
1
)
)
β
+
2
]
(
β
+
1
)
(
β
+
2
)
ssssssss
+
(
b
+
2
R
a
)
[
S
3
-
(
T
(
j
-
1
)
+
λ
)
3
]
6
s
s
s
s
s
s
s
s
+
(
b
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
S
(
2
+
3
β
+
β
2
)
-
2
λ
-
T
(
j
-
1
)
β
+
2
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
(
T
(
j
-
1
)
+
λ
)
sssssssssssssssssssssss
b
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
×
(
2
+
3
β
+
β
2
)
-
2
λ
)
)
(
b
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
]
)
s
s
s
s
s
s
s
s
×
(
(
β
+
1
)
(
β
+
3
)
(
β
+
2
)
2
)
-
1
s
s
s
s
s
s
s
s
+
(
a
2
[
(
T
(
j
-
1
)
+
λ
)
β
+
1
(
(
T
(
j
-
1
)
+
λ
)
(
β
)
+
2
λ
+
T
(
j
-
1
)
)
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
S
β
+
λ
+
K
T
j
)
-
T
(
j
-
1
)
β
+
1
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
a
2
[
(
T
(
j
-
1
)
+
λ
)
β
+
1
(
(
T
(
j
-
1
)
+
λ
)
(
β
)
+
2
λ
+
T
(
j
-
1
)
)
)
×
(
(
T
(
j
-
1
)
+
λ
)
(
β
)
+
2
λ
+
T
(
j
-
1
)
)
]
)
s
s
s
s
s
s
s
s
×
(
(
β
+
1
)
(
β
+
2
)
)
-
1
s
s
s
s
s
s
s
s
+
α
a
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
2
(
β
+
1
)
2
s
s
s
s
s
s
s
s
-
b
a
2
(
β
+
1
)
(
β
+
3
)
(
β
+
2
)
s
s
s
s
s
s
s
s
×
[
(
s
-
λ
)
β
+
1
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
λ
S
(
β
+
1
)
+
K
2
T
2
j
2
β
(
β
+
3
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
+
2
λ
2
+
2
K
2
T
2
j
2
)
s
s
s
s
s
s
s
s
s
s
s
s
-
(
T
(
j
-
1
)
)
β
+
1
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
λ
(
T
(
j
-
1
)
+
λ
)
(
β
+
1
)
+
F
2
β
(
β
+
3
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
(
s
-
λ
)
β
+
1
+
2
λ
2
+
2
(
T
(
j
-
1
)
+
λ
)
2
)
]
s
s
s
s
s
s
s
s
-
a
c
[
(
S
-
λ
)
β
+
1
-
(
T
(
j
-
1
)
)
β
+
1
]
β
+
1
s
s
s
s
s
s
s
s
+
(
R
a
α
[
(
S
-
λ
)
β
+
2
(
S
(
β
+
2
)
+
λ
)
-
(
T
(
j
-
1
)
)
β
+
2
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
×
(
(
T
(
j
-
1
)
+
λ
)
(
β
+
2
)
+
λ
)
]
)
s
s
s
s
s
s
s
s
×
(
(
β
+
1
)
(
β
+
3
)
(
β
+
2
)
)
-
1
s
s
s
s
s
s
s
s
-
R
b
[
K
4
T
4
j
4
-
(
T
(
j
-
1
)
+
λ
)
4
]
8
)
.Present Value of Shortage Cost. The present value of shortage cost is similar to the previous case.


### 4.2.2. The Present Value of Purchasing Cost

The total amount that we have to purchase forjth cycle is
(30)
I
1
(
j
)
(
T
j
)
+
I
3
(
j
-
1
)
(
T
j
)
.
The lead time for this amount of item is l
j
=
l
~
w
j.So the total present value of purchasing cost is(31)
PVCP
=
C
3
I
11
(
0
)
+
C
3
∑
j
=
1
n
-
1
(
I
1
(
j
)
(
T
j
)
-
I
3
(
j
-
1
)
(
T
j
)
)
×
e
-
R
~
(
T
j
-
l
~
(
(
I
1
(
j
)
(
T
j
)
-
I
3
(
j
-
1
)
(
T
j
)
)
)
.


### 4.2.3. The Present Value of Ordering Cost

For thejth cycle we have to order at T
(
j
-
1
)
-
l
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
).So the total present value of ordering cost for thejth cycle is
(32)
PVCO
=
∑
j
=
1
m
C
A
j
=
C
4
∑
j
=
1
m
e
-
R
~
(
T
j
-
l
~
(
I
3
j
(
T
j
)
+
I
2
j
(
T
j
)
)
)
.


### 4.2.4. Present Value of Deteriorated Item Cost

The amount of deteriorated items during[
T
(
j
-
1
)
+
λ
,
T
(
j
-
1
)
+
k
T
] is
(33)
I
2
j
(
λ
)
-
∫
λ
S
D
(
t
)
d
t
=
a
(
s
-
λ
)
+
α
a
(
s
-
λ
)
β
+
1
β
+
1
-
b
[
(
s
-
λ
)
2
2
+
λ
(
s
-
λ
)
s
s
s
s
s
s
s
s
+
α
(
(
s
-
λ
)
β
+
2
β
+
2
+
λ
(
s
-
λ
)
β
+
1
β
+
1
)
]
-
a
(
K
T
j
-
λ
)
+
b
(
k
2
T
2
j
2
-
λ
2
)
2
.
The total present value of the amount of deteriorated items cost is
(34)
PVCD
=
∑
j
=
1
m
C
D
j
=
θ
~
C
5
∑
j
=
1
m
∫
T
(
j
-
1
)
+
λ
K
T
j
(
e
-
α
(
t
-
λ
)
β
∫
t
S
(
a
-
b
u
)
e
α
(
u
-
λ
)
β
d
u
)
e
-
R
~
t
d
t
.
The total present value of costs over time horizon in this case (TPV) is
(35)
TPV
=
PVCP
+
PVCH
+
PVCS
+
PVCO
+
PVCD
.


## 4.2.1. Present Value of Holding Cost

The present value of holding cost forjth cycle is given by
(27)
C
H
j
=
C
1
(
∫
T
(
j
-
1
)
T
(
j
-
1
)
+
λ
I
1
j
(
t
)
e
-
R
~
t
d
t
+
∫
T
(
j
-
1
)
+
λ
K
T
j
I
2
j
(
t
)
e
-
R
~
t
d
t
)
=
C
1
(
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
sssssss
+
1
R
~
3
[
b
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
-
R
~
(
T
(
j
-
1
)
+
λ
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
+
R
~
(
T
(
j
-
1
)
+
λ
)
)
-
2
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
+
e
-
R
~
λ
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
×
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
]
ssssssss
+
∫
T
(
j
-
1
)
+
λ
S
(
e
-
α
(
t
-
λ
)
β
∫
t
S
(
a
-
b
u
)
sssssssssssssssssssssss
×
∫
t
S
(
a
-
b
u
)
e
α
(
u
-
λ
)
β
d
u
)
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
e
-
R
~
t
d
t
)
.
For the simplified formulation we define
(28)
C
=
a
s
+
a
α
(
s
-
λ
)
β
+
1
β
+
1
-
b
2
S
2
-
b
α
(
s
-
λ
)
β
+
1
(
s
b
+
s
-
λ
)
(
β
+
2
)
(
β
+
1
)
.
Therefore we have
(29)
C
H
j
≅
C
1
(
[
G
e
-
R
~
T
(
j
-
1
)
R
~
(
1
-
e
-
R
~
λ
)
]
+
1
R
~
3
s
s
s
s
s
s
s
×
[
b
e
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
s
s
s
s
s
s
s
s
s
s
×
(
-
R
~
(
T
(
j
-
1
)
+
λ
)
(
2
+
R
~
(
T
(
j
-
1
)
+
λ
)
)
-
2
)
s
s
s
s
s
s
s
s
s
s
+
e
-
R
~
λ
(
(
R
~
(
T
(
j
-
1
)
)
(
2
+
R
~
(
T
(
j
-
1
)
)
)
+
2
)
)
]
s
s
s
s
s
s
s
+
C
(
K
T
-
λ
)
-
(
R
C
+
a
)
2
[
S
2
-
T
(
j
-
1
)
+
λ
2
]
s
s
s
s
s
s
s
-
a
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
-
(
T
(
j
-
1
)
)
β
+
2
]
(
β
+
1
)
(
β
+
2
)
ssssssss
+
(
b
+
2
R
a
)
[
S
3
-
(
T
(
j
-
1
)
+
λ
)
3
]
6
s
s
s
s
s
s
s
s
+
(
b
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
S
(
2
+
3
β
+
β
2
)
-
2
λ
-
T
(
j
-
1
)
β
+
2
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
(
T
(
j
-
1
)
+
λ
)
sssssssssssssssssssssss
b
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
×
(
2
+
3
β
+
β
2
)
-
2
λ
)
)
(
b
α
[
(
T
(
j
-
1
)
+
λ
)
β
+
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
]
)
s
s
s
s
s
s
s
s
×
(
(
β
+
1
)
(
β
+
3
)
(
β
+
2
)
2
)
-
1
s
s
s
s
s
s
s
s
+
(
a
2
[
(
T
(
j
-
1
)
+
λ
)
β
+
1
(
(
T
(
j
-
1
)
+
λ
)
(
β
)
+
2
λ
+
T
(
j
-
1
)
)
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
(
S
β
+
λ
+
K
T
j
)
-
T
(
j
-
1
)
β
+
1
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
a
2
[
(
T
(
j
-
1
)
+
λ
)
β
+
1
(
(
T
(
j
-
1
)
+
λ
)
(
β
)
+
2
λ
+
T
(
j
-
1
)
)
)
×
(
(
T
(
j
-
1
)
+
λ
)
(
β
)
+
2
λ
+
T
(
j
-
1
)
)
]
)
s
s
s
s
s
s
s
s
×
(
(
β
+
1
)
(
β
+
2
)
)
-
1
s
s
s
s
s
s
s
s
+
α
a
2
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
2
(
β
+
1
)
2
s
s
s
s
s
s
s
s
-
b
a
2
(
β
+
1
)
(
β
+
3
)
(
β
+
2
)
s
s
s
s
s
s
s
s
×
[
(
s
-
λ
)
β
+
1
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
λ
S
(
β
+
1
)
+
K
2
T
2
j
2
β
(
β
+
3
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
+
2
λ
2
+
2
K
2
T
2
j
2
)
s
s
s
s
s
s
s
s
s
s
s
s
-
(
T
(
j
-
1
)
)
β
+
1
s
s
s
s
s
s
s
s
s
s
s
s
×
(
2
λ
(
T
(
j
-
1
)
+
λ
)
(
β
+
1
)
+
F
2
β
(
β
+
3
)
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
(
s
-
λ
)
β
+
1
+
2
λ
2
+
2
(
T
(
j
-
1
)
+
λ
)
2
)
]
s
s
s
s
s
s
s
s
-
a
c
[
(
S
-
λ
)
β
+
1
-
(
T
(
j
-
1
)
)
β
+
1
]
β
+
1
s
s
s
s
s
s
s
s
+
(
R
a
α
[
(
S
-
λ
)
β
+
2
(
S
(
β
+
2
)
+
λ
)
-
(
T
(
j
-
1
)
)
β
+
2
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
[
(
S
-
λ
)
2
(
β
+
1
)
-
(
T
(
j
-
1
)
)
2
(
β
+
1
)
]
×
(
(
T
(
j
-
1
)
+
λ
)
(
β
+
2
)
+
λ
)
]
)
s
s
s
s
s
s
s
s
×
(
(
β
+
1
)
(
β
+
3
)
(
β
+
2
)
)
-
1
s
s
s
s
s
s
s
s
-
R
b
[
K
4
T
4
j
4
-
(
T
(
j
-
1
)
+
λ
)
4
]
8
)
.Present Value of Shortage Cost. The present value of shortage cost is similar to the previous case.


## 4.2.2. The Present Value of Purchasing Cost

The total amount that we have to purchase forjth cycle is
(30)
I
1
(
j
)
(
T
j
)
+
I
3
(
j
-
1
)
(
T
j
)
.
The lead time for this amount of item is l
j
=
l
~
w
j.So the total present value of purchasing cost is(31)
PVCP
=
C
3
I
11
(
0
)
+
C
3
∑
j
=
1
n
-
1
(
I
1
(
j
)
(
T
j
)
-
I
3
(
j
-
1
)
(
T
j
)
)
×
e
-
R
~
(
T
j
-
l
~
(
(
I
1
(
j
)
(
T
j
)
-
I
3
(
j
-
1
)
(
T
j
)
)
)
.


## 4.2.3. The Present Value of Ordering Cost

For thejth cycle we have to order at T
(
j
-
1
)
-
l
~
(
I
3
(
j
-
1
)
(
T
(
j
-
1
)
)
+
I
1
j
(
T
(
j
-
1
)
)
).So the total present value of ordering cost for thejth cycle is
(32)
PVCO
=
∑
j
=
1
m
C
A
j
=
C
4
∑
j
=
1
m
e
-
R
~
(
T
j
-
l
~
(
I
3
j
(
T
j
)
+
I
2
j
(
T
j
)
)
)
.


## 4.2.4. Present Value of Deteriorated Item Cost

The amount of deteriorated items during[
T
(
j
-
1
)
+
λ
,
T
(
j
-
1
)
+
k
T
] is
(33)
I
2
j
(
λ
)
-
∫
λ
S
D
(
t
)
d
t
=
a
(
s
-
λ
)
+
α
a
(
s
-
λ
)
β
+
1
β
+
1
-
b
[
(
s
-
λ
)
2
2
+
λ
(
s
-
λ
)
s
s
s
s
s
s
s
s
+
α
(
(
s
-
λ
)
β
+
2
β
+
2
+
λ
(
s
-
λ
)
β
+
1
β
+
1
)
]
-
a
(
K
T
j
-
λ
)
+
b
(
k
2
T
2
j
2
-
λ
2
)
2
.
The total present value of the amount of deteriorated items cost is
(34)
PVCD
=
∑
j
=
1
m
C
D
j
=
θ
~
C
5
∑
j
=
1
m
∫
T
(
j
-
1
)
+
λ
K
T
j
(
e
-
α
(
t
-
λ
)
β
∫
t
S
(
a
-
b
u
)
e
α
(
u
-
λ
)
β
d
u
)
e
-
R
~
t
d
t
.
The total present value of costs over time horizon in this case (TPV) is
(35)
TPV
=
PVCP
+
PVCH
+
PVCS
+
PVCO
+
PVCD
.


## 5. Model Analysis

The problem is to determine the optimal values ofn and k so as to minimize the total present value of inventory system costs:
(36)
TPV
=
PVCP
+
PVCH
+
PVCS
+
PVCO
+
PVCD
.
For this, the algorithm begins by setting discrete variable n
=
1 and takes the derivatives of TPV (
n
,
k
) with respect to k. Equating the derivative to zero derives the following necessary conditions of optimality ∂
T
P
V
/
∂
k
=
0.For a given value ofn derive k
* from the above equation TPV (
n
,
k
*
) is derived by substituting (
n
,
k
*
) into (18). Then, n increases by increment of one continually and TPV (
n
,
k
*
) is driven again. The above stages repeat until the minimum TPV (
n
,
k
*
) can be found.Since the inflation rate, deterioration rate, and the lead time are fuzzy numbers the objective function becomes fuzzy. We first use the singed distance method to defuzzify the objective function. Then since the derived function is complicated we use the pervious algorithm for optimizations system.


## 6. Numerical Example

The numerical example is provided to clarify how the proposed model is applied  (see Table1).Table 1
The parameter values for numerical example.
C
1
0.5
C
2
3.5
C
3
3
C
4
50
l
~
(0.000434, 0.000651, 0.000868)
θ
~
(0.086, 0.11, 0.205)
R
~
(0.5, 0.1, 0.15)
α
0.05
β
3
λ
0.028For the first case the result of the numerical example is shown in Table2.Table 2
Data for numerical example (Case  1).
N
k
TPV
1
0.30065
113917.42
2*
0.29804*
43788.11*
3
0.29315
48725.02
4
0.29307
55435.08
5
0.29104
90183.22For the second case the result of the numerical example is shown in Table3.Table 3
Data for numerical example (Case  2).
N
k
TPV
1*
0.2054*
4.6
×
10
8
*
2
0.2448
6
.
8
×
10
8
3
0.27315
7.2
×
10
8
4
0.31023
9.3
×
10
8
5
0.35002
1.2
×
10
9


## 7. Conclusion

In this paper, an inventory model with time-dependent demand studied under two different deterioration rates: (a) time-dependent (Weibull) deterioration rate and (b) fuzzy deterioration rate. The items are not deteriorating as soon as they received to the inventory but they start to deteriorate after a specified interval. The lead time is a function of the amount that we ordered; we consider both inflation rate and the lead time for preparing one item as fuzzy parameter. Total present costs of inventory system were minimized and numerical examples were proposed to illustrate findings. Different from the previous studies here we are able to consider the rate of deterioration as two kinds beside each other and compare the solutions of these two cases with each other when all of the other parameters are the same. The results show that the total cost in the second case (considering stochastic deterioration) is more than the first case (considering fuzzy deterioration).



---

*Source: 347857-2014-08-14.xml*