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.
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.
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))).
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 systemCovertR. P.PhilipG. C.An EOQ model for items with Weibull distribution deteriorationMoonI.GiriB. C.KoB.Economic order quantity models for ameliorating/deteriorating items under inflation and time discountingChenJ.-M.An inventory model for deteriorating items with time-proportional demand and shortages under inflation and time discountingSanaS. S.The stochastic EOQ model with random sales priceSanaS. S.Optimal selling price and lotsize with time varying deterioration and partial backloggingBalkhiZ. T.An optimal solution of a general lot size inventory model with deteriorated and imperfect products, taking into account inflation and time value of moneyDeS. K.SanaS. S.Fuzzy order quantity inventory model with fuzzy shortage quantity and fuzzy promotional indexLoS.-T.WeeH.-M.HuangW.-C.An integrated production-inventory model with imperfect production processes and Weibull distribution deterioration under inflationSanaS. S.The EOQ model—a dynamical systemRoy ChowdhuryR.GhoshS. K.ChaudhuriK. S.An inventory model for perishable items with stock and advertisement sensitive demandSinghT.PattnayakH.Two-warehouse inventory model for deteriorating items with linear demand under conditionally permissible delay in paymentSivashankariC. K.PanayappanS.Production inventory model with reworking of imperfect production, scrap and shortagesBuzacottJ. A.Economic order quantity with inflationSoleimani-AmiriM.NodoustS.MirzazadehA.MohammadiM.The average annual cost and discounted cost mathematical modeling methods in the inflationary inventory systems—a comparison analysisHouK.An inventory model for deteriorating items with stock-dependent consumption rate and shortages under inflation and time discountingChernM.YangH.TengJ. L.Partial backlogging inventory lot-size models for deteriorating items with fluctuating demand under inflationHorowitzI.EOQ and inflation uncertaintyMirzazadehA.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 conditionsMirzazadehA.Inventory management under stochastic conditions with multiple objectivesAmeliM.MirzazadehA.ShiraziM. A.Economic order quantity model with imperfect items under fuzzy inflationary conditionsJanaD. 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 approachGholami-QadikolaeiA.MirzazadehA.Tavakkoli-MoghaddamR.A stochastic multiobjective multiconstraint inventory model under inflationary condition and different inspection scenariosNeetuA. K.TomerA.Deteriorating inventory model under variable inflation when supplier credits linked to order quantityMirzazadehA.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 shortagesAbadP. L.Optimal price and order size for a reseller under partial backorderingAbadP. L.Optimal pricing and lot-sizing under conditions of perishability and partial backorderingHouK. L.LinL. C.Optimal pricing and ordering policies for deteriorating items with multivariate demand under trade credit and inflationSarkarB.SanaS. S.ChaudhuriK.An imperfect production process for time varying demand with inflation and time value of money: an EMQ modelMousaviS. M.HajipourV.NiakiS. T.AlikarN.Optimizing multi-item multi-period inventory control system with discounted cash flow and inflation: two calibrated meta-heuristic algorithmsSinghS. R.GuptaV.GuptaP.Three stage supply chain model with two warehouse, imperfect production, variable demand rate and inflationGhoreishiM.Arshsadi KhamsehA.MirzazadehA.Joint optimal pricing and inventory control for deteriorating items under inflation and customer returnsHsiehT.DyeC. Y.Pricing and lot-sizing policies for deteriorating items with partial backlogging under inflationDattaT. K.PalA. K.Effects of inflation and time-value of money on an inventory model with linear time-dependent demand rate and shortagesWuK. S.An EOQ inventory model for items with Weibull distribution deterioration, ramp type demand rate and partial backloggingGiriB. C.JalanA. K.ChaudhuriK. S.Economic order quantity model with Weibull deterioration distribution, shortage and ramp-type demand