- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Advances in Operations Research
Volume 2013 (2013), Article ID 973125, 13 pages
A Partial Backlogging Inventory Model for Deteriorating Item under Fuzzy Inflation and Discounting over Random Planning Horizon: A Fuzzy Genetic Algorithm Approach
1Department of Applied Sciences, Haldia Institute of Technology, Haldia, Purba Medinipur 721657, India
2Department of Mathematics, Sidho-Kanho-Birsha University, Purulia 723101, India
3Department of Mathematics, Bengal Engineering and Science University, Howrah 711103, India
Received 5 December 2012; Accepted 2 June 2013
Academic Editor: Koji Shingyochi
Copyright © 2013 Dipak Kumar Jana et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
An inventory model for deteriorating item is considered in a random planning horizon under inflation and time value money. The model is described in two different environments: random and fuzzy random. The proposed model allows stock-dependent consumption rate and shortages with partial backlogging. In the fuzzy stochastic model, possibility chance constraints are used for defuzzification of imprecise expected total profit. Finally, genetic algorithm (GA) and fuzzy simulation-based genetic algorithm (FSGA) are used to make decisions for the above inventory models. The models are illustrated with some numerical data. Sensitivity analysis on expected profit function is also presented. Scope and Purpose. The traditional inventory model considers the ideal case in which depletion of inventory is caused by a constant demand rate. However, to keep sales higher, the inventory level would need to remain high. Of course, this would also result in higher holding or procurement cost. Also, in many real situations, during a longer-shortage period some of the customers may refuse the management. For instance, for fashionable commodities and high-tech products with short product life cycle, the willingness for a customer to wait for backlogging is diminishing with the length of the waiting time. Most of the classical inventory models did not take into account the effects of inflation and time value of money. But in the past, the economic situation of most of the countries has changed to such an extent due to large-scale inflation and consequent sharp decline in the purchasing power of money. So, it has not been possible to ignore the effects of inflation and time value of money any more. The purpose of this paper is to maximize the expected profit in the random planning horizon.
In the past few decades, many researches have studied an inventory model with constant demand or dynamic demand (cf. M. K. Maiti and M. Maiti , Taleizadeh et al. , Jana et al. , and others). Moreover, in a competitive situation attractive display of units in the showroom is an important factor. Levin et al.  noted that at times the presence of inventory has a motivational effect on the people around it. It is a common belief that large piles of goods displayed in a supermarket will lead the customer to buy more. Thus, many business people use showrooms and the attractive display of units in the showroom to influence the customers. Roy et al.  and Maiti  have developed an inventory model with stock-dependent demand.
In most of the earlier inventory models, lifetime of an item is assumed to be infinite while it is in storage. But in reality, many physical goods deteriorate due to dryness, spoilage, vaporization, and so forth and are damaged due to hoarding longer than their normal storage period. The deterioration also depends on preserving facilities and environmental conditions in warehouses/storage. So, due to deterioration effect, a certain fraction of the items either damaged or decayed are not in perfect condition to satisfy the future demand of customers for good items. Deterioration for such items is continuous and constant or time-dependent and/or dependent on the on-hand inventory. A number of research papers have already been published on the above type of items by Roy et al.  and others.
Moreover, the effects of inflation and time value of money are vital in practical situation, especially in the developing countries with large-scale inflation. Therefore, the effect of inflation and time value of money cannot be ignored in real situations. To relax the assumption of non inflationary effects on costs, Buzacott  and Misra  simultaneously developed an EOQ model with a constant inflation rate for all associated costs. Bierman and Thomas  then proposed an EOQ model under inflation that also incorporated the discount rate. Misra  then extended the EOQ model with different inflation rates for various associated costs. Recently, Chern et al.  proposed partial backlogging inventory lot-size models for deteriorating items with fluctuating demand under inflation. Maity and Maiti  have developed a multiobjective optimal inventory control problem for deteriorating multi-items under fuzzy inflation and discounting. Yang et al.  proposed an inventory model under inflation for deteriorating items with stock-dependent consumption rate and partial backlogging shortages.
Use of GA in complex decision making problem is well established Michalewicz . A simple GA starts with a set of potential solutions (called initial population) of the decision making problem under consideration. Individual solutions are called chromosome. Crossover and mutation operations happen among the potential solutions with some probability and respectively, to get a new set of solutions and it continues until terminating conditions are encountered. Behavior and performance of a GA is directly affected by the interaction between the parameters, that is, selection process of chromosomes for mating pool, () and so forth. Poor parameter settings usually lead to several problems such as premature convergence. Extensive research work has been made to improve the performance of GA for single/multiobjective continuous/discrete optimization problems during the last two decades. Michalewicz  proposed a genetic algorithm, named contractive mapping genetic algorithm (CMGA), where movement from old population to new population takes place only when average fitness of new population is better than the old one and proved the asymptotic convergence of the algorithm by Banach fixed-point theorem. Bessaou and Siarry  proposed a GA where initially more than one population of solutions are generated. Genetic operations are done on every population a finite number of times to find a promising zone of optimum solution. Finally, a population of solutions is generated in this zone and genetic operations are done on this population a finite number of times to get a final solution. Last and Eyal  developed a GA with varying population size, where chromosomes are classified into young, middle-aged, and old according to their age and lifetime. Genotype diversity and phenotype diversity of the final population are obtained to measure the performance of the GA. Pezzella et al.  developed a GA for the flexible Job-shop scheduling problem, which integrates different strategies for generating the initial population, selecting the individuals for reproduction and reproducing new individuals. In this research paper, an EPQ model of an item is developed in a random planning horizon; that is, lifetime of the product is assumed as random in nature and it follows an exponential distribution with known mean. Unit production cost decreases in each production cycle due to learning effects of the workers on production. Similarly setup cost in each cycle is partly constant and partly decreases in each cycle due to learning effects of the employees. Model is formulated to maximize the expected profit from the whole planning horizon. Following Last and Eyal , a GA with varying population size is implemented where chromosomes are classified into young, middle-aged and old according to their age and lifetime. In this GA, crossover probability is a function of parents' age type (young, middle-aged, old, etc.) and is obtained using a fuzzy rule base and fuzzy possibility theory Dubois and Prade . It is an improved GA where a subset of better children is included with the parent population for next generation, and size of this subset is a percentage of the size of its parent set. This GA is used to make optimal decision for the above production inventory model. Performance of the proposed GA for solving the model is compared with that of basic GA and CMGA. The model is illustrated with some numerical data. Sensitivity analysis on expected profit function is also presented. Recently, many papers have been developed in GA (cf. Narmatha Banu and Devaraj , Kar et al. ).
Due to fuzzy inflation rate, the objective function is fuzzy in nature, and then following Liu and Iwamura , M. K. Maiti and M. Maiti , the said objective function is converted to a crisp objective function. Some research papers have been already published considering imprecise planning horizon (cf. M. K. Maiti and M. Maiti , Roy et al. [5, 22], and others). Till now, none has developed inventory models with both random planning horizon and imprecise effect due to inflation and discounting and stock-dependent demand. In this paper, a partial backlogging inventory model with stock-dependent demand for deterioration item has been developed under imprecise inflation rate over a random planning horizon. Here, the planning horizon is stochastic in nature and follows exponential distribution. The expected profit is maximized using a FSGA with roulette wheel selection, arithmetic crossover, and random mutation (Michalewicz, ). The models are illustrated with some numerical data. Sensitivity analysis on expected profit is presented.
2. Assumptions and Notations
The mathematical model of the inventory replenishment problem is based on the following assumptions. Demand rate is assumed to depend on the existing stock level. The time horizon (a random variable) is finite. The time horizon completely accommodates first cycles and ends during th cycle. Lead time is negligible. Replenishment rate is infinite but replenishment size is finite. Shortages are allowed. Unsatisfied demand is partially backlogged. The fraction of shortages backordered is a differentiable and decreasing function of time , denoted by , where is the waiting time up to the next replenishment, and with . Note that if (or 0) for all , then shortages are completely backlogged (or lost). For deteriorating items, a constant fraction of the on-hand inventory deteriorates per unit of time and there is no repair or replacement of the deteriorated inventory during the planning period.
For convenience, the following notations are used throughout the entire paper. : on-hand inventory of a cycle at time , .: the demand rate, where : number of fully accommodated cycles to be made during the prescribed time horizon. : duration of a complete cycle. : total ordered quantity in th cycle. : selling price per unit. : is setup cost in th cycle, ; ; . Here , , and are so chosen to best fit the setup cost function. It is also noted that, set-up cost decreases with the number of cycle due to the learning effect.: the external variable purchasing cost per unit. : the inventory holding cost per unit per unit time. : the backlogging cost per unit per unit time. : the deterioration rate per unit per unit time. : the discount rate. : the inflation rate, which is varied by the social economical situations. is the discount rate minus the inflation rate.: the backlogging rate which is a decreasing function of the waiting time , we here assume that , where , and is the waiting time. : holding cost in the th cycle. : holding cost in the last cycle. : sales revenue in the th cycle. : sales revenue in the last cycle. : purchasing cost in the th cycle. : purchasing cost in the last cycle. : total profit after completing fully accommodated cycles. : total profit for the last cycles. : expected total profit after completing fully accommodated cycles. : expected total profit from the last cycle. : expected total profit from the planning horizon. : total time horizon (a random variable) and is the finite time horizon. Here, it is assumed that the planning horizon is a random variable and follows exponential distribution with probability density function (p.d.f) as follows: Here, is the parameter of the distribution. : the amount of backorders at time . : order quantity at time .
3. Mathematical Formulation
In the development of the model, we assume that there are full cycles during the random time horizon and the planning horizon ends within th cycle, that is, within and . The th replenishment is made at time . The quantity received at is used partly to meet the accumulated backorders in the previous cycle. The inventory at gradually reduces to zero at (cf. Figure 1). For the last cycle, some amount may be left after the end of planning horizon. This amount is sold at a reduced price in a lot.
3.1. Formulation for th Cycle
The differential equation describing the inventory level in the interval is given by where and .
Subject to the conditions that at .
The solution of the differential equation (3) is given by During the time interval , the demand rate and backlogged rate . Hence, the amount of backorders is governed by the following differential equation: subject to the conditions that at . The solutions of the differential equation (5) is given by From (4) and (6), we have the order quantity at in the th replenishment cycle: Present value of holding cost of the inventory for the th cycle is given by Present value of ordering cost for the th cycle is given by Present value of purchasing cost for the th cycle is given by Present value of sales revenue for the th cycle is given by Then, the present value of the total profit from the full cycle is as follows:
Now, So, Since the planning horizon has a p.d.f , the present value of expected total profit from complete cycles is given by
3.2. Formulation for Last Cycle
The differential equations describing the inventory level in the interval are given by where , subject to the conditions that The solutions of the differential equations (15) and (16) are given by From (17), we have the order quantity at in the th replenishment cycle: Present value of ordering cost is given by Present value of purchasing cost is given by In last cycle, for simplicity we consider two cases only depending upon the cycle length. Let be the real value corresponding to the random variable .
Case 1 (). Present value of holding cost of the inventory for the last cycle (cf. Figure 2) is given by Present value of sales revenue is given by
Case 2 (). Present value of holding cost of the inventory for the last cycle (Figure 3) is given by Present value of sales revenue is given by So, expected holding cost for the last cycle is given by Expected sales revenue from the last cycle is given by Expected reduced selling price from the last cycle is given by Expected setup cost from the last cycle is given by So, expected total profit from last cycle is given by
3.3. Total Profit from the System
Now, total expected profit from the complete time horizon is given by
4. Models in Different Environments
4.1. Stochastic Model (Model 1)
When the resultant effect of inflation and discounting is crisp in nature, then the present problem is to determine and so as to
4.1.1. Fuzzy Stochastic Model (Model 2)
In the real world, deterioration () and rate of inflation () are imprecise in nature, that is, vaguely defined in some situations. So we take , as fuzzy number, that is, as and . Then, due to this assumption, our objective function becomes . Since optimization of a fuzzy objective is not well defined, so instead of one can optimize its equivalent optimistic or pessimistic return of the objective function. When decision maker likes to optimize the optimistic equivalent of , then the problem reduces to the determination of , so as to Following Liu and Iwamura  and others, it can be defuzzified and rewritten as (for details see Das et al. ).
5. Solution Procedure
To solve the stochastic model 1, GA is used. The basic technique to deal with problems (35) is to convert the possibility constraint to its deterministic equivalent. However, the procedure is usually very hard and successful in some particular cases. Following Liu and Iwamura , M. K. Maiti and M. Maiti , here two simulation algorithms are proposed to determine in (35), for a feasible .
To determine , for feasible variables, roughly find a point from fuzzy number , which approximately minimizes . Let this value be and set . (For simplicity one can take .) Then, is randomly generated in -cut set of and let value of for and if replace with . This step is repeated a finite number of times and final value is taken as value of . This phenomenon is used to develop Algorithm 1.
Now roughly find a point from fuzzy number , which approximately minimizes . Let this value be (for simplicity one can take also) and a positive number. Set and if , then increase with . Again check and it continues until . At this stage, decrease the value of and again try to improve . When becomes sufficiently, small, then we stop, and the final value of is taken as value of . Using this criterion, required algorithm is developed as shown in Algorithm 2. In this algorithm, the variable is used to store initial assumed value of and is used to store the value of in each iteration.
Therefore for feasible value of the variables, we determine using the above algorithms and to optimize we use GA. Since fuzzy simulation algorithm is used to determine in the algorithm, this GA is named a fuzzy simulation based genetic algorithm (FSGA). This algorithm is named FSGA when fuzzy simulation process is used to determine objective function value.
5.1. Fuzzy Simulation Based Single Objective Genetic Algorithm (FSGA)
In natural genesis, we know that chromosomes are the main carriers of the hereditary information from parents to offsprings and that genes, which present hereditary factors, are lined up in chromosomes. At the time of reproduction, crossover and mutation take place among the chromosomes of parents. In this way, hereditary factors of parents are mixed up and carried over to their offsprings. Darwinian principle states that only the fittest animals can survive in nature. So a pair of fittest parents normally reproduce better offspring.
The above-mentioned phenomenon is followed to create a genetic algorithm for an optimization problem. Here potential solutions of the problem are analogous with the chromosomes and chromosome of better offspring with the better solution of the problem. Crossover and mutation occur among a set of potential solutions and obtained a new set of solutions and it continues until terminating conditions are encountered. Michalewicz  proposed a genetic algorithm named the contractive mapping genetic algorithm (CMGA) and proved the asymptotic convergence of the algorithm by the Banach fixed-point theorem. In CMGA, movement from an old population to a new population takes place only when the average fitness of a new population is better than the old one. This algorithm is modified with the help of a fuzzy simulation process to solve the model in some cases. The algorithm is named FSGA and this is presented below. In this algorithm, , are probabilities of the crossover and the probability of mutation, respectively, is the iteration counter, and is the population of potential solutions for iteration . The initialize function initializes the population at the time of initialization. The evaluate function evaluates the fitness of each member of , and at this stage an objective function value due to each solution is evaluated via the fuzzy simulation process (using Algorithm 1).
5.1.1. FSGA Algorithm
See Algorithm 3.
5.1.2. FSGA Procedures
Representation: An “-dimensional real vector” is used to represent a solution, where represent decision variables of the problem.Initialization: such solutions are randomly generated by random number generator such that each satisfies the constraints of the problem. Constraints of the problem are satisfied using Algorithm 1. This solution set is taken as initial population . Also set , , and .Fitness value: value of the objective function due to the solution is taken as fitness of . Let it be . Objective function is evaluated via fuzzy simulation process (using Algorithm 1 or Algorithm 2).Selection process for mating pool: the following steps are followed for this purpose. Find total fitness of the population . Calculate the probability of selection of each solution by the formula . Calculate the cumulative probability for each solution by the formula . Generate a random number “” from the range . If , then select otherwise select where . Repeat step, (iv) and (v) times to select solutions from old population. Clearly, one solution may be selected more than once. Selected solution set is denoted by in the proposed FSGA algorithm. CrossoverSelection for crossover: for each solution of generate a random number from the range . If , then the solution is taken for crossover, where is the probability of crossover. Crossover process: crossover took place on the selected solutions. For each pair of coupled solutions , , a random number is generated from the range , and their offsprings and are obtained by the formula: MutationSelection for mutation: for each solution of generate a random number from the range . If , then the solution is taken for mutation, where is the probability of mutation. Mutation process: to mutate a solution select a random integer in the range . Then, replace by randomly generated value within the boundary of th component of .
6. Numerical Illustration
To solve the stochastic model (Model 1), genetic algorithm (GA) (Section 5.1) is used and fuzzy stochastic model (Model 2) is solved by fuzzy simulation based genetic algorithm (FSGA) (Section 6). The corresponding parameters in GA and FSGA are POPSIZE = 50, PCROS = 0.2, PMUTE = 0.2, and MAXGEN = 50. A real-number presentation is used here. In this representation, each chromosome X is a string of (here, ) number of genes, these represent decision variables. For each chromosome X, every gene, which represents the independent variables (here, ), is randomly generated between its boundaries until it is feasible. In this problem, arithmetic crossover and random mutation are applied to generate new offsprings.
To illustrate the models, we consider the following numerical data: , , , , , , , , , , , and ; that is, in appropriate units.
6.1. Stochastic Model 1
The optimal values of along with maximum expected total profit have been calculated for different values of and ; results are displayed in Table 1.
It is observed that as and increase, expected profit increases due to the increase of demand henceforth selling amount. Moreover for increasing values of and , as length of business periods decreases, average expected profit increases. All these observations agree with reality.
6.2. Fuzzy Stochastic Model 2
Here, the resultant inflationary effect and deterioration rate are considered as a triangular fuzzy number; that is, and for , , and all other data remain the same as in stochastic model. The maximum optimistic returns have been calculated for different levels of optimistic.
6.3. Sensitivity Analysis
A sensitivity analysis is performed for stochastic model with respect to different resultant inflationary effect () for crisp inflation, and results are presented in Table 4. It is observed that as increases profit decreases which agrees with reality.
A sensitivity analysis is performed for the maximum expected total profit with respect to the different values of parameter for stochastic model and is presented in Table 5. It is observed that as decreases profit increases. This happens because as decreases expected time horizon increases which increases the total expected profit.
Results due to different values of confidence levels for Model 2 is calculated and depicted in Figure 2. In both cases, as expected, profit decreases with the increase of confidence levels. The graphical representation of possibility threshold versus expected profit is depicted in Figure 4.
7. Concluding Remarks
In this paper, a realistic stock-dependent inventory model with shortages and partial backlogging has been formulated with fuzzy deterioration, inflation rate, and setup cost with learning effect in random planning horizon. Until now, no inventory model has been formulated in such consideration, that is, stock-dependent demand, shortages with partial backlogging, effect of inflation, learning effect, deterioration, random planning horizon, and so forth, in both stochastic and fuzzy stochastic environments. The proposed models are optimized via soft computing methods: GA and FSGA. The present concept can be extended to production planning model, multi-item production planning model, optimal control problem for multiproduct manufacturing also, and so forth, which may be areas of future research.
The authors are very grateful to the respected editors of this journal for improvement of the paper. This work is supported by the University Grants Commission, India, vide F.PSW-103/10-11 (ERO) dated 20.10.10.
- M. K. Maiti and M. Maiti, “Fuzzy inventory model with two warehouses under possibility constraints,” Fuzzy Sets and Systems, vol. 157, no. 1, pp. 52–73, 2006.
- A. A. Taleizadeh, S. T. A. Niaki, and R. Nikousokhan, “Constraint multiproduct joint-replenishment inventory control problem using uncertain programming,” Applied Soft Computing Journal, vol. 11, no. 8, pp. 5143–5154, 2011.
- D. K. Jana, K. Maity, and T. K. Roy, “A bi-fuzzy approach to a production-recycling-disposal inventory problem with environment pollution cost via genetic algorithm,” International Journal of Computer Applications, vol. 61, pp. 1–10, 2013.
- R. I. Levin, C. P. Mcaughlim, P. R. Lamone, and J. F. Kottas, Production Management/ Operations Management: (Contemporary Policy for Managing Operating System), McGraw-Hill, New York, NY, USA, 1972.
- A. Roy, M. K. Maiti, S. Kar, and M. Maiti, “Two storage inventory model with fuzzy deterioration over a random planning horizon,” Mathematical and Computer Modelling, vol. 46, no. 11-12, pp. 1419–1433, 2007.
- M. K. Maiti, “A fuzzy genetic algorithm with varying population size to solve an inventory model with credit-linked promotional demand in an imprecise planning horizon,” European Journal of Operational Research, vol. 213, no. 1, pp. 96–106, 2011.
- J. A. Buzacott, “Economic order quantities with inflation,” Operational Research Quarterly, vol. 26, no. 3, pp. 553–558, 1975.
- R. B. Misra, “A study of inflationary effects on inventory systems,” Logistic Spectrum, vol. 9, pp. 260–268, 1975.
- H. Bierman and J. Thomas, “Inventory decisions under inflationary conditions,” Decision Sciences, vol. 8, pp. 151–155, 1997.
- R. B. Misra, “Note on optimal inventory management under inflation,” Naval Research Logistics Quarterly, vol. 26, no. 1, pp. 161–165, 1979.
- M.-S. Chern, H.-L. Yang, J.-T. Teng, and S. Papachristos, “Partial backlogging inventory lot-size models for deteriorating items with fluctuating demand under inflation,” European Journal of Operational Research, vol. 191, no. 1, pp. 127–141, 2008.
- K. Maity and M. Maiti, “A numerical approach to a multi-objective optimal inventory control problem for deteriorating multi-items under fuzzy inflation and discounting,” Computers & Mathematics with Applications, vol. 55, no. 8, pp. 1794–1807, 2008.
- H. L. Yang, J. T. Teng, and M. S. Chern, “An inventory model under inflation for deteriorating items with stock-dependent consumption rate and partial backlogging shortages,” International Journal of Production Economics, vol. 123, pp. 8–19, 2010.
- Z. Michalewicz, Genetic Algorithms + Data Structures = Evolution Programs, Springer, Berlin, Germany, 1992.
- M. Bessaou and P. Siarry, “A genetic algorithm with real-value coding to optimize multimodal continuous functions,” Structural and Multidisciplinary Optimization, vol. 23, no. 1, pp. 63–74, 2001.
- M. Last and S. Eyal, “A fuzzy-based lifetime extension of genetic algorithms,” Fuzzy Sets and Systems, vol. 149, no. 1, pp. 131–147, 2005.
- F. Pezzella, G. Morganti, and G. Ciaschetti, “A genetic algorithm for the flexible job-shop scheduling problem,” Computers and Operations Research, vol. 35, no. 10, pp. 3202–3212, 2008.
- D. Dubois and H. Prade, Fuzzy Sets and Systems, Theory and Applications, vol. 144, Academic Press, New York, NY, USA, 1980.
- R. Narmatha Banu and D. Devaraj, “Multi-objective GA with fuzzy decision making for security enhancement in power system,” Applied Soft Computing, vol. 12, pp. 2756–2764, 2012.
- S. Kar, D. Das, and A. Roy, “A production-inventory model for a deteriorating item incorporating learning effect using genetic algorithm,” Advances in Operations Research, vol. 2010, Article ID 146042, 26 pages, 2010.
- B. Liu and K. Iwamura, “Chance constrained programming with fuzzy parameters,” Fuzzy Sets and Systems, vol. 94, no. 2, pp. 227–237, 1998.
- A. Roy, M. K. Maiti, S. Kar, and M. Maiti, “An inventory model for a deteriorating item with displayed stock dependent demand under fuzzy inflation and time discounting over a random planning horizon,” Applied Mathematical Modelling, vol. 33, no. 2, pp. 744–759, 2009.
- B. Das, K. Maity, and M. Maiti, “A two warehouse supply-chain model under possibility/ necessity/credibility measures,” Mathematical and Computer Modelling, vol. 46, no. 3-4, pp. 398–409, 2007.