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.
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:
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;
deterioration which refers to the damage, spoilage, dryness, vaporization, and so forth of the products;
amelioration which refers to items whose value or utility or quantity increases with time;
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.
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 (λ).
The system is under the control of fuzzy inflation.
We introduce the fuzzy function of order quantity for lead time.
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 Section 2. 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.
The demand rate is a linear function of time:
(1)D(t)=a-bta>0,b>0.
The replenishment rate is infinite.
Shortages are allowed and fully backlogged except in the last cycle.
The discount rate is r.
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 pt is pte-R~t, t≥0 where pt is the value of p at time t.
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].
Lead time is a function of the total ordered quantity in a period.
Lj is the lead time for periodj:
(2)lj~=a~*wj,
where wj 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 (a1,a2,a3).
θ 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:
finite horizon time (time units),
number of replenishment cycles during the time horizon H,
length of each cycle (T=H/n),
inventory level at any time t in periodj, i=1,2,3, j=1,2,…,n,
holding cost per unit per unit time at time zero ($/unit/unit time),
shortage cost per unit per unit time at time zero ($/unit/unit time),
purchasing cost per unit item at time zero ($),
ordering cost per order at time zero ($/set-up).
Additional notations will be introduced later.
3. DefinitionDefinition 1.
If X is a collection of objects denoted generically by x, then a fuzzy set P~ on X is a set of order pairs.
If x 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.
If A~=(a1,a2,a3) is a triangular fuzzy number (TFN), then μA~(x) is defined as follows: (a1, a2, and a3 are real numbers)
(3)μA~(x)={x-a1a2-a1fora1≤x≤a2a3-xa3-a2fora2≤x≤a30otherwise.
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~α=[AL(α)α,AR(α)α].
Here AL(α)α and AR(α)α are the lower and upper bounds of the closed interval and
(4)AL(α)α=inf{x∈R:μA~(x)≥α},AR(α)α=sup{x∈R:μA~(x)≥α}.
Definition 3 (signed distance).
For any a∈R we define d0(a,0)=a. If a>0 the distance from a to 0 is d0(a,0)=a. If a<0, the distance from a to 0 is -d0(a,0)=-a. This is the reason why we call d0(a,0) the signed distance from a to 0.
The signed distances of AL(α) and AR(α) to 0 are d0(AL(α),0)=AL(α) and d0(AR(α),0)=AR(α), respectively. We define the signed distance from 0 to the interval [Al(α),AR(α)] as
(5)d0([Al(α),AR(α)],0)=12(d0(AL(α),0),d0(AR(α),0)),12(AL(α)+AR(α))=12[a+c+(2b-a-c)α].
Since crisp interval [AL(α),AR(α)] has a one-to-one correspondence with level α fuzzy interval [AL(α)α,AR(α)α], it is natural to define the signed distance from a fuzzy interval [AL(α)α,AR(α)α] to 0 as
(6)d([AL(α)α,AR(α)α],0)=d0([Al(α),AR(α)],0)=12[a+c+(2b-a-c)α].
This is a continuous function of α on 0≤α≤1. We can find its average value through integration. Since A~=⋃0≤α≤1[AL(α)α,AR(α)α], we have the following definition.
Definition 4.
If A~=(a,b,c) is a triangular fuzzy number, the signed distance from A~ to 0~ is defined as
(7)d(A~,0~)=∫01d([AL(α)α,AR(α)α],0)dα=14(2b+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~=[m1α,m2α]+[n1α,n2α]=[m1α+n1α,m2α+n2α],M~-N~=[m1α,m2α]-[n1α,n2α]=[m1α-n1α,m2α-n2α],M~*N~=[m1α,m2α]*[n1α,n2α]=[m1α*n1α,m2α*n2α],M~÷N~=[m1α,m2α]÷[n1α,n2α]=[m1α÷n2α,m2α÷n1α].
4. Model Formulation
In the development of the model, we assume that there are N 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 for jth (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 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 amount I1j(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 kT, 0≤k≤1 inventory level leads to zero and shortages occur. And at the moment (KTj-lj~) of every jth (j=1,2,…,n) cycle, company purchases an amount Qj units of the item. During the time interval [kTj,Tj] 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 of jth cycle there is no deterioration (λ unit of time after the replenishment time). The inventory level, I1j(t), satisfies the following differential equations:
(9)∂I1j(t)∂t=-a+btT(j-1)≤t≤T(j-1)+λ.
The second part of the inventory level can be described as follows:
(10)∂I2j(t)∂t=-a+bt-θ~I2j(t)T(j-1)+λ≤t≤kTj.
We can solve the above differential equations by using the following equations as a boundary condition:
(11)I2j(T(j-1)+kT)=0.
We can see
(12)I2j(t)=-bθ2-a-btθ+eθ(S-t)(a-bSθ+bθ2),I2j(T(j-1)+λ)=bθ2-a+b[T(j-1)+λ]θ+eθ[T(j(k-1)+1)-λ]+a(T(j-1)+λ)+12b(T(j-1)+λ)2,I2j(T(j-1)+λ)=I1j(T(j-1)+λ),I1j(t)=bθ2-a+b[T(j-1)+λ]θ+eθ[T(j(k-1)+1)-λ]+a(T(j-1)+λ)+12b(T(j-1)+λ)2-(at+bt22).
For the last part of jth cycle, the shortages level is governed by
(13)∂I3j(t)∂t=-(a+bt),kTj≤t≤jT,I3(Tkj)=0,I3j(t)=akTj+bk2T2j2-(at+0.5bt2),kTj≤t≤jT,I31(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 for jth cycle is given by
(15)PVCHj=C1(∫T(j-1)T(j-1)+λI1j(t)e-R~tdtsssssss+∫T(j-1)+λT(j-1)+kTI2j(t)e-R~tdt)=C1[Ge-R~T(j-1)R~(1-e-R~λ)]+1R~3[be-R~(T(j-1)+λ)((R~(T(j-1))(2+R~(T(j-1)))+2))sssssssssss×(-R~(T(j-1)+λ)sssssssssssssss×(2+R~(T(j-1)+λ))-2)sssssssssss+e-R~λ((R~(T(j-1))×((2+R~(T(j-1)))+2)sssssssssssssssssss×(2+R~(T(j-1)))+2))e-R~(T(j-1)+λ)]+e-R~(T(j-1)+λ+KTj)R~2*[e-R~KTj(b+aR~+bR~(T(j-1)+λ))sssse-R~KTj+(b+aR~+bR~KTj)]+bR~θ2[e-R~(T(j-1)+λ)-e-R~KTj]+eθKTj[a+bKTjθ-bθ2]*(eKTj(R~+θ)-e(T(j-1)+λ)(R~+θ)R~+θ).
So the total present value of holding cost is
(16)PVCH=∑j=1nPVCHj.
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)PVCSj=C2(∫KTjTjI3j(t)e-R~tdt)=aKTj+bk2T2j2R~(e-R~KTj-e-R~Tj)-ae-R~Tj(1+k)2R2(e-R~Tj(1+R~Tj(1+k))ssssssssssssssssss-e-R~Tj(1+k)(1+R~Tj))+-be-R~Tj(1+k)2R~3×[e-R~Tj(2+R~KTj(2+R~KTj))(2+R~Tj(1+k)(2+Tj(1+k)))ssssss-e-R~Tj(1+k)(2+R~Tj(1+k)(2+Tj(1+k)))].
So the total present value of shortage cost is
(18)PVCS=∑j=1n-1PVCSj.
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 the jth cycle (j=1,2,…,n-1) the present value of purchasing cost can be formulated as follows:
(19)PVCPj=C3(I3(j-1)(T(j-1))+I1j(T(j-1))e-R~(I3(j-1)(T(j-1))+I1j(T(j-1)))sssssss+I1j(T(j-1))e-R~(I3(j-1)(T(j-1))+I1j(T(j-1)))).
Therefore, the total present value of purchasing cost is
(20)CP=∑j=1mC3(I3(j-1)(T(j-1))ssssssssss+I1j(T(j-1)))e-R~(I3(j-1)(T(j-1))+I1j(T(j-1))).
4.1.4. Present Value of Ordering Cost
Present value of sales revenue for the jth (1≤j≤N) cycle is given by
(21)CA=∑j=1mCAj=C4∑j=1me-R~(Tj-l~(I3j(Tj)+I2j(Tj))).
4.2. Formulation for jth (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 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-λ)β-1e-α(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 of jth cycle satisfies the following differential equation:
(24)∂I1j(t)∂t=-a+bt,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)∂I2j(t)∂t=-a+bt-αβ(t-λ)β-1I2(t),T(j-1)+λ≤t≤T(j-1)K,S=T(j-1)+KT,I2j(S)=0,I2j(t)=e-α(t-λ)β∫tS(a-bu)eα(u-λ)βdu,T(j-1)+λ≤t≤T(j-1)+KT.
So for the I2j(λ) we have
(26)I2j(λ)=∫λS(a-bu)(∑n=0∞(α(u-λ)β)nn!)≅a(s-λ)+αa(s-λ)β+1β+1-b[(s-λ)22+λ(s-λ)ssssssss+α((s-λ)β+2β+2+λ(s-λ)β+1β+1)]=Qj,I1j(λ)=I2j(λ)=Qj,I1j(t)=bt22-at+Qj-bλ22+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 for jth cycle is given by
(27)CHj=C1(∫T(j-1)T(j-1)+λI1j(t)e-R~tdt+∫T(j-1)+λKTjI2j(t)e-R~tdt)=C1([Ge-R~T(j-1)R~(1-e-R~λ)]sssssss+1R~3[be-R~(T(j-1)+λ)((R~(T(j-1))(2+R~(T(j-1)))+2))ssssssssssssss×(-R~(T(j-1)+λ)ssssssssssssssssss×(2+R~(T(j-1)+λ))-2)ssssssssssssss+e-R~λ((R~(T(j-1))(2+R~(T(j-1)))+2)ssssssssssssssssssssse-R~(T(j-1)+λ)×(2+R~(T(j-1)))+2))]ssssssss+∫T(j-1)+λS(e-α(t-λ)β∫tS(a-bu)sssssssssssssssssssssss×∫tS(a-bu)eα(u-λ)βdu)[Ge-R~T(j-1)R~(1-e-R~λ)]e-R~tdt).
For the simplified formulation we define
(28)C=as+aα(s-λ)β+1β+1-b2S2-bα(s-λ)β+1(sb+s-λ)(β+2)(β+1).
Therefore we have
(29)CHj≅C1([Ge-R~T(j-1)R~(1-e-R~λ)]+1R~3sssssss×[be-R~(T(j-1)+λ)((R~(T(j-1))(2+R~(T(j-1)))+2))ssssssssss×(-R~(T(j-1)+λ)(2+R~(T(j-1)+λ))-2)ssssssssss+e-R~λ((R~(T(j-1))(2+R~(T(j-1)))+2))]sssssss+C(KT-λ)-(RC+a)2[S2-T(j-1)+λ2]sssssss-aα[(T(j-1)+λ)β+2-(T(j-1))β+2](β+1)(β+2)ssssssss+(b+2Ra)[S3-(T(j-1)+λ)3]6ssssssss+(bα[(T(j-1)+λ)β+2[(S-λ)2(β+1)-(T(j-1))2(β+1)]ssssssssssssssss×(S(2+3β+β2)-2λ-T(j-1)β+2ssssssssssssssssssss×((T(j-1)+λ)sssssssssssssssssssssssbα[(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)]])ssssssss×((β+1)(β+3)(β+2)2)-1ssssssss+(a2[(T(j-1)+λ)β+1((T(j-1)+λ)(β)+2λ+T(j-1)))ssssssssssssssss×(Sβ+λ+KTj)-T(j-1)β+1sssssssssssssssa2[(T(j-1)+λ)β+1((T(j-1)+λ)(β)+2λ+T(j-1)))×((T(j-1)+λ)(β)+2λ+T(j-1))])ssssssss×((β+1)(β+2))-1ssssssss+αa2[(S-λ)2(β+1)-(T(j-1))2(β+1)]2(β+1)2ssssssss-ba2(β+1)(β+3)(β+2)ssssssss×[(s-λ)β+1ssssssssssss×(2λS(β+1)+K2T2j2β(β+3)ssssssssssssssss+2λ2+2K2T2j2)ssssssssssss-(T(j-1))β+1ssssssssssss×(2λ(T(j-1)+λ)(β+1)+F2β(β+3)sssssssssssssss(s-λ)β+1+2λ2+2(T(j-1)+λ)2)]ssssssss-ac[(S-λ)β+1-(T(j-1))β+1]β+1ssssssss+(Raα[(S-λ)β+2(S(β+2)+λ)-(T(j-1))β+2sssssssssssssssss[(S-λ)2(β+1)-(T(j-1))2(β+1)]×((T(j-1)+λ)(β+2)+λ)])ssssssss×((β+1)(β+3)(β+2))-1ssssssss-Rb[K4T4j4-(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 for jth cycle is
(30)I1(j)(Tj)+I3(j-1)(Tj).
The lead time for this amount of item is lj=l~wj.
So the total present value of purchasing cost is
(31)PVCP=C3I11(0)+C3∑j=1n-1(I1(j)(Tj)-I3(j-1)(Tj))×e-R~(Tj-l~((I1(j)(Tj)-I3(j-1)(Tj))).
4.2.3. The Present Value of Ordering Cost
For the jth cycle we have to order at T(j-1)-l~(I3(j-1)(T(j-1))+I1j(T(j-1))).
So the total present value of ordering cost for the jth cycle is
(32)PVCO=∑j=1mCAj=C4∑j=1me-R~(Tj-l~(I3j(Tj)+I2j(Tj))).
4.2.4. Present Value of Deteriorated Item Cost
The amount of deteriorated items during [T(j-1)+λ,T(j-1)+kT] is
(33)I2j(λ)-∫λSD(t)dt=a(s-λ)+αa(s-λ)β+1β+1-b[(s-λ)22+λ(s-λ)ssssssss+α((s-λ)β+2β+2+λ(s-λ)β+1β+1)]-a(KTj-λ)+b(k2T2j2-λ2)2.
The total present value of the amount of deteriorated items cost is
(34)PVCD=∑j=1mCDj=θ~C5∑j=1m∫T(j-1)+λKTj(e-α(t-λ)β∫tS(a-bu)eα(u-λ)βdu)e-R~tdt.
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 of n 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 ∂TPV/∂k=0.
For a given value of n 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 Table 1).
The parameter values for numerical example.
C1
0.5
C2
3.5
C3
3
C4
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.028
For the first case the result of the numerical example is shown in 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.22
For the second case the result of the numerical example is shown in Table 3.
Data for numerical example (Case 2).
N
k
TPV
1*
0.2054*
4.6×108*
2
0.2448
6.8×108
3
0.27315
7.2×108
4
0.31023
9.3×108
5
0.35002
1.2×109
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).
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
GhareP. M.SchraderG. H.A model for exponentially decaying inventory system196321449460CovertR. P.PhilipG. C.An EOQ model for items with Weibull distribution deterioration19735432332610.1080/056955573089749182-s2.0-0015726898MoonI.GiriB. C.KoB.Economic order quantity models for ameliorating/deteriorating items under inflation and time discounting2005162377378510.1016/j.ejor.2003.09.025MR2104064ZBL1067.900042-s2.0-9544222622ChenJ.-M.An inventory model for deteriorating items with time-proportional demand and shortages under inflation and time discounting1998551213010.1016/S0925-5273(98)00011-52-s2.0-0032496519SanaS. S.The stochastic EOQ model with random sales price2011218223924810.1016/j.amc.2011.05.015MR28204862-s2.0-79960838435SanaS. S.Optimal selling price and lotsize with time varying deterioration and partial backlogging2010217118519410.1016/j.amc.2010.05.040MR2672576ZBL1231.900552-s2.0-77955422425BalkhiZ. T.An optimal solution of a general lot size inventory model with deteriorated and imperfect products, taking into account inflation and time value of money2004352879610.1080/002077203100016570452-s2.0-2142662250DeS. K.SanaS. S.Fuzzy order quantity inventory model with fuzzy shortage quantity and fuzzy promotional index201331135135810.1016/j.econmod.2012.11.0462-s2.0-84871947116LoS.-T.WeeH.-M.HuangW.-C.An integrated production-inventory model with imperfect production processes and Weibull distribution deterioration under inflation2007106124826010.1016/j.ijpe.2006.06.0092-s2.0-33845508611SanaS. S.The EOQ model—a dynamical system2012218178736874910.1016/j.amc.2012.02.032MR29213612-s2.0-84859432600Roy ChowdhuryR.GhoshS. K.ChaudhuriK. S.An inventory model for perishable items with stock and advertisement sensitive demand20149316917710.1080/17509653.2013.867097SinghT.PattnayakH.Two-warehouse inventory model for deteriorating items with linear demand under conditionally permissible delay in payment201492104113SivashankariC. K.PanayappanS.Production inventory model with reworking of imperfect production, scrap and shortages201491920BuzacottJ. A.Economic order quantity with inflation197526355355810.1057/jors.1975.1132-s2.0-0016557640Soleimani-AmiriM.NodoustS.MirzazadehA.MohammadiM.The average annual cost and discounted cost mathematical modeling methods in the inflationary inventory systems—a comparison analysis2014113141HouK.An inventory model for deteriorating items with stock-dependent consumption rate and shortages under inflation and time discounting2006168246347410.1016/j.ejor.2004.05.011MR2170249ZBL1101.900062-s2.0-24944521899ChernM.YangH.TengJ. L.Partial backlogging inventory lot-size models for deteriorating items with fluctuating demand under inflation2008191112714110.1016/j.ejor.2007.03.053MR24325652-s2.0-43049141197HorowitzI.EOQ and inflation uncertainty200065221722410.1016/S0925-5273(99)00034-12-s2.0-0033907763MirzazadehA.SarfarazA. R.Constrained multiple items optimal order policy under stochastic inflationary conditionsProceedings of the 2nd Annual International Conference on Industrial Engineering Application and Practice1997San Diego, Calif, USA725730MirzazadehA.A comparison of the mathematical modeling methods in the inventory systems under uncertain conditions2011361316142MirzazadehA.Inventory management under stochastic conditions with multiple objectives20132162610.12677/AIRR.2013.21003AmeliM.MirzazadehA.ShiraziM. A.Economic order quantity model with imperfect items under fuzzy inflationary conditions20116329430310.3923/tasr.2011.294.303JanaD. K.DasB.RoyT. K.A partial backlogging inventory model for deteriorating item under fuzzy inflation and discounting over random planning horizon: a fuzzy genetic algorithm approach2013201313973125MR308556510.1155/2013/973125Gholami-QadikolaeiA.MirzazadehA.Tavakkoli-MoghaddamR.A stochastic multiobjective multiconstraint inventory model under inflationary condition and different inspection scenarios2013227710571074NeetuA. K.TomerA.Deteriorating inventory model under variable inflation when supplier credits linked to order quantity20123812411263MirzazadehA.Seyyed EsfahaniM. M.Fatemi GhomiS. M. T.An inventory model under certain inflationary conditions, finite production rate and inflation-dependent demand rate for deteriorating items with shortages2009401213110.1080/00207720802088264MR24899742-s2.0-60749107006AbadP. L.Optimal price and order size for a reseller under partial backordering2001281536510.1016/S0305-0548(99)00086-62-s2.0-0033791466AbadP. L.Optimal pricing and lot-sizing under conditions of perishability and partial backordering19964281093110410.1287/mnsc.42.8.1093ZBL0879.900692-s2.0-0001021905HouK. L.LinL. C.Optimal pricing and ordering policies for deteriorating items with multivariate demand under trade credit and inflation2013503404417MR3105537SarkarB.SanaS. S.ChaudhuriK.An imperfect production process for time varying demand with inflation and time value of money: an EMQ model20113811135431354810.1016/j.eswa.2011.04.0442-s2.0-79959940448MousaviS. M.HajipourV.NiakiS. T.AlikarN.Optimizing multi-item multi-period inventory control system with discounted cash flow and inflation: two calibrated meta-heuristic algorithms20133742241225610.1016/j.apm.2012.05.019MR30023142-s2.0-84870240357SinghS. R.GuptaV.GuptaP.Three stage supply chain model with two warehouse, imperfect production, variable demand rate and inflation201341819210.5267/j.ijiec.2012.010.0052-s2.0-84871631460GhoreishiM.Arshsadi KhamsehA.MirzazadehA.Joint optimal pricing and inventory control for deteriorating items under inflation and customer returns20132013770908310.1155/2013/709083HsiehT.DyeC. Y.Pricing and lot-sizing policies for deteriorating items with partial backlogging under inflation201037107234724210.1016/j.eswa.2010.04.0042-s2.0-78049353790DattaT. K.PalA. K.Effects of inflation and time-value of money on an inventory model with linear time-dependent demand rate and shortages199152332633310.1016/0377-2217(91)90167-T2-s2.0-0026173567WuK. S.An EOQ inventory model for items with Weibull distribution deterioration, ramp type demand rate and partial backlogging200112878779310.1080/095372801100518192-s2.0-0035782366GiriB. C.JalanA. K.ChaudhuriK. S.Economic order quantity model with Weibull deterioration distribution, shortage and ramp-type demand200334423724310.1080/0020772131000158500MR20091752-s2.0-0142186093