The European Applied Business Research Conference Rothenburg, Germany 00 Using Genetic Algorithms To Optimise Kanban-Based Production System M Al-Tahat, (E-mail: mohammadal-tahat@mailinguniboit), University of Bologna, Italy A Liverani, (E-mail: alfredoliverani@mailinguniboit), University of Bologna, Italy F Persiani, University of Bologna, Italy G Caligiana, University of Bologna, Italy Abstract Evolutionary computational methods as Genetic Algorithms GAs have proven to be a robust search technique for solving deterministic problems GA handles its search by giving out a population of solutions for an optimisation problem based on principles from evolution This paper describes the use of GA strategies in a stochastic process of a Kanban based manufacturing system, to optimise Work In Progress WIP holding cost The GA search determines the minimum WIP holding cost and the corresponding Kanban numbers and production rate required to meet customer demand rate The procedures are illustrated with an applied problem consisting of many decision variables Insights are provided on the effect of the population size, maximum generation of offspring solutions and string length on the fitness of the solutions The solutions found by the GA search indicate that the GA search provides good solutions for optimising Kanban problems Introduction In previous work a method for the purpose of determining the Kanban system parameters was introduced for a multi-product multi-stage manufacturing environment [5] Consequently we have introduced in [6] a Kanban based marketing tool, which enhances the sale department efficiency The solutions in both mentioned work were found by a traditional computing way through developing proprietary computer program commands, and were found based on the formulation approach for Kanban based production lines, which had been presented by Bhaba in terms of total production cost [] Here the production total cost will be investigated deeply; particularly the cost of WIP will be analysed As it was observed previously in Bhaba model the elements of the total production cost are: Raw material inventory holding cost, Cost of WIP inventory Tc; and finished product inventory cost The scope of this research is to analyse the second cost element cost of WIP -- and to apply Genetic Algorithms GA approach to find the best fitness for the cost of WIP inventory [4] Genetic Algorithms GAs I Rechenberg introduced idea of evolutionary computing in the 960s in his work "Evolution strategies" (Evolutions strategies in original) Other researchers then developed his idea GAs were invented by John Holland and developed by him and his students and colleagues at the university of Michigan [] This leads to Holland's book "Adoption in Natural and Artificial Systems" published in 975 In 99 John Koza has used genetic algorithm to evolve programs to perform certain tasks He called his method "genetic programming" GP LISP programs were used, because programs in this language can be expressed in the form of a "parse tree, which is the obect the GA works on [] All living organisms consist of cells In each cell there is the same set of chromosomes Chromosomes are strings of DNA and serves as a model for the whole organism A chromosome consists of genes, blocks of DNA Each gene encodes a particular protein Basically it can be said, that each gene encodes a trait, for example colour of eyes Possible settings for a trait (eg blue, brown) are called alleles Each gene has its own position in the chromosome This position is called locus Complete set of genetic material (all chromosomes) is called genome Particular set of genes in genome is called genotype The genotype constitutes, with later development after birth, the base for the organism's phenotype, it determines its physical and mental characteristics, such as eye colour, intelligence etc
The European Applied Business Research Conference Rothenburg, Germany 00 During reproduction, first of all occurs recombination (or crossover) Genes from parents form, in some way, the whole rearrange new chromosome The new created offspring can then be mutated Mutation means, that the elements of DNA can be locally changed at random These changes are mainly caused by errors in copying genes from parents The fitness of an organism is measured by success of the organism in its life GAs is stochastic search methods that mimic the behaviour of natural biological evolution Genetic algorithms operate on a population of potential solutions applying the principle of survival of the fittest to produce better and better approximations to a solution At each generation, a new set of approximations is created by the process of selecting individuals according to their level of fitness in the problem domain and of breeding them together using operators borrowed from natural genetics This process leads to the evolution of populations of individuals that are better suited to their environment than the individuals by which they were created from, ust as in natural adaptation [7] GAs model natural processes, such as selection, recombination, mutation, migration, locality, and neighbourhood Figure [], which is found in Internet websites shows the structure of a simple genetic algorithm Evolutionary algorithms work on populations of individuals instead of single solutions In this way the search is performed in a parallel manner At the beginning as indicated in the figure, population is considered made by number of individuals randomly chosen The obective function is then evaluated for these individuals The first/initial generation is produced If the optimisation criteria are not met the creation of a new generation starts Figure []: Structure of Single Population Genetic Algorithm Individuals are selected according to their fitness for the production of offspring Parents are recombined to produce offspring All offspring will be mutated with a certain probability The fitness of the offspring is then computed The offspring are inserted into the population replacing the parents, producing a new generation This cycle is performed until the optimisation criteria is satisfied Such a single population evolutionary algorithm is powerful and performs well on a broad class of problems However, better results can be obtained by introducing many populations, called subpopulations Every subpopulation evolves for a few generations isolated (like the single population evolutionary algorithm) before one or more individuals are exchanged between the subpopulations The Multi population evolutionary algorithm models the evolution of a species in a way more similar to nature than the single population evolutionary algorithm Problem Definition The following notations have been used in order to formulate the interactions of the different production cost parameters:
The European Applied Business Research Conference Rothenburg, Germany 00 The obective of this work is to use GA to optimise the total WIP holding cost in a Kanban based production line and to find the corresponding number of Kanban If the average WIP inventory between the two workstations is represented by W avg and calculated as shown below [5] then the total holding cost of WIP inventory Tc in terms of the number of Kanban k employed in Kanban stage, for a N -stage Kanban system will be given by the following [] N Tc= kw + H = When substituting the number of Kanban ( Wavg) Wavg = T N P 0 T pt = + + + k k k k W in Tc the final detailed equation for calculating the WIP holding cost in terms of avg k employed in Kanban stage, for a N -stage Kanban system will be as below: N H P0 T pt Tc = k + + + + W = T k k k k Where the values of T, T are given by: D0 ± D0 + T = ω * ω * Q R P0 ± T = 0 P + * p * D p F It is obviously shown that the total Tc is the cumulative summation of WIP cost for all Kanban stages
The European Applied Business Research Conference Rothenburg, Germany 00 Therefore WIP cost at a particular Kanban stage for a particular product will be found by: Tc H P0T pt = + + k W + + T k k k k As a result, the general fitness function for calculating the WIP cost for a product in a general Kanban stage between two adacent workstations-- Sc could be written as in the following formula This function will be utilized by genetic algorithms approach to achieve our research obectives Sc H PT 0 pt + + kw + + T k k k k = Solving of Fitness Function We will try to minimize the solutions of the above formula by genetic algorithm As shown, the cost function implies many input variables A computer program has been used to solve the best fitness of the WIP inventory cost (minimizing Sc ) Franco Persiani, Gianni Caligiana and Gian Marco Saggiani had been developed this computer program in the faculty of engineering at Bologna University System parameters used in this program are shown in Figure [] The system uses the following genetic algorithm as described by David EGoldberg [] Randomly select an initial population of popsize parent solutions with a proper number of genes ngp for all the variables npar0 of the fitness functions; therefore the chromosomes or string length chrom will be determined Popsize is recommended to be odd valued to serve the idea of elitism It score the initial selected solutions with the fitness function Sc Using roulette wheel selection method It selects two parent s solutions (chromosomes) from the population, according to their fitness, to become parents for crossover and mutating 4 After parent s selections, the system makes a step to crossover Crossover selects genes from parent chromosomes and creates a new offspring by choosing randomly single crossover point with pcross probability of crossover 5 After a crossover is performed, mutation with probability of pmutation takes place This is to prevent all solutions in population to fall into a local optimum of the solution value Mutation changes randomly the new offspring 6 The program comes back to step () until the desired number of offspring solutions is generated maxgen=max 7 All new offspring solutions Score with the fitness function Sc 8 The program returns to step () until termination criterion met 4
The European Applied Business Research Conference Rothenburg, Germany 00 Variables of fitness function Population Chromosome npar0= npar0= npar0=npar0 Size No: ngp0= ngp0= ngp0= ngp0= ngp0 ngp0= ngp0= popsize= 0 0 0 popsize= 0 0 0 0 0 popsize= 0 0 0 0 popsize= popsize popsize 0 0 0 0 Figure []: System Parameters of the Developed Computer Program ngp0= ngp0= ngp0 Starting Up the System A value for unit holding cost of WIP inventory is assumed equal to H =0 ($/unit/unit time) and a value for the cost of withdrawing Kanban between two-workstation is set equal to W = 0($/Kanban) so that, the fitness function will be as below ( + + ) ( + P0 P0 pdf P P0 P0 + pdf ) 0ω Sc = 0k + 0 + + + + + p ( D 0 D0 ωq ) k k k k R This fitness function Sc as shown is function of seven fitness variables; these are: k, p, ω, P0, DF, D0andQ R Note that the quantity of materials consumed during the production period is equal to the amount of raw materials supplied during that period of time for one-to-one conversion ratio (ie, f = D F / Q R = ) Furthermore initial values for all program system parameters are assumed The assumed initial system parameters indicated in Figure [] are described by a population of eleven chromosomes (solutions); each chromosome contains (string length) eight genes (alleles); a maximum of 500 offspring solutions (iterations) for the seven fitness variables and the corresponding fitness value is considered, A 005 probability of mutation and 06 probability of crossover is utilized to select the minimum fitness value among those solutions Limits in Figure [] are assumed for the seven variables fitness function 5
The European Applied Business Research Conference Rothenburg, Germany 00 6 kanban = variable []=k 0 drate =variable []=w prate = variable []=p 000 pratet0 = variable [4]=P0 0000 fprdemand = variable [5]=DF 000 dratet0 = variable [6]=D0 0000 rmorderedqty = variable [7]=QR W P Figure []: Variables of Fitness Function Limits The initial rough output results obtained based on the above assumptions are summarized in Figure [4] These results are ust to start up and to test the algorithm and have no physical interpretations The fitness variable values have been verified with the above corresponding fitness value using Microsoft excel and the same results were obtained Fitness The values for k,w,p,p0,df,d0 that give the minimum fitness noted that DF=QR are: {,,45,996,,,} respectively With a corresponding cost value (Sc)=000 For the following chromosome structure: {{{0,0,0,0,0,0,0,0}=k} {{0,0,0,0,0,0,0,0}=w} {{0,0,,0,0,0,0,}=p} {{,,,,,,,0}=P0} {{0,0,0,0,0,0,0,0}=DF} {{0,0,0,0,0,0,0,0}=D0}} (Off- Figure [4]: Initial Solutions Solutions Perfection Starting from the above initial solutions with the same assumptions, the model has been studied and analysed under different values of system parameters That had been done in order to get the perfect system parameters The program has been run and results have been obtained to all of the following cases: Case 0 Popsize=,maxgen=500,ngp0=8 Case 0 Popsize=5,maxgen=500,ngp0=8 Case 0 Popsize=,maxgen=500,ngp0=8 Case 04 Popsize=7,maxgen=500,ngp0=8 Case 05 Popsize=,maxgen=500,ngp0= Case 06 Popsize=,maxgen=500,ngp0=6 6
The European Applied Business Research Conference Rothenburg, Germany 00 Case 07 Popsize=,maxgen=500,ngp0=0 Case 08 Popsize=,maxgen=500,ngp0=4 Case 09 Popsize=,maxgen=400,ngp0=4 Case 0 Popsize=,maxgen=00,ngp0=4 Case Popsize=,maxgen=00,ngp0=4 Case Popsize=,maxgen=00,ngp0=4 Case Popsize=7,maxgen=00,ngp0=4 Case 4 Popsize=7,maxgen=00,ngp0=0 Case 5 Popsize=7,maxgen=00,ngp0=6 Case 6 Popsize=7,maxgen=00,ngp0=4 It has been found that the best system parameters to be used are those in case 6 where Popsize=7 maxgen=00 and ngp0=4 The corresponding solutions for these parameters based on 00 offspring solutions are summarized bellow in Figure [5] Fitness 0 0 008 006 004 00 (A): Fitness behaviour p 0 8 6 4 (B): Change of prduction rate 0 40 60 80 00 Iter No 0 40 60 80 00 Iter No k 4 (C):Change of Kanban w (D): Change of demand rate 0 8 6 0 40 60 80 00 Iter No 0 40 60 80 00 Iter No The values for k,w,p,p0,df,d0 that give the minimum fitness noted that DF=QR are: {,,5,000,,,} respectively With a corresponding cost value (Sc)=0000 For the following chromosome structure: {{{00000000000000000000}=k} {{000000000000000000}=w} {{000000000000000000000}=p} {{00000000000000000}=P0} {{000000000000000000000000000}=DF} {{0000000000000000000000}=D0}} Kanban System Design Figure [5]: Solutions for Best System Parameters Through the above sections we clearly reveal a methodology of using genetic algorithms to minimize the WIP cost by demonstrating a general solutions under specific assumptions This methodology will be used here with its corresponding assumptions to design and operate a particular Kanban-based production line given a selected production parameters values Given a values for fitness variables that give the minimum WIP cost ω, Q R andd 7 F the system will calculate the best values of the other
The European Applied Business Research Conference Rothenburg, Germany 00 The problem then is to design a Kanban production line--by finding the Kanban numbers and sizes required-- and to plan the operations of this line by finding the values of initial production and initial demand quantities-- Assume ω = 8, Q R = 500 andd F = 500 Under this situation Kanban model could be expressed by the following fitness Consequently the variables limits of fitness function will assumed to be as shown in Figure [6] ( + + 05 p) ω Sc = 0k + p + 5 p p *5000 05ω + 064 k + + + k k k + + 05 The solutions of the above problem for popsize =7,maxgen= 00, ngp0=4, 005 probability of mutation and 06 probability of crossover are shown in Figure [7] It is obviously shown in the figure that the best Kanban numbers are two with size equal to 50 unit of production for each Kanban The corresponding Machining centre for this stage of production should have a capacity of at least 9 (unit/ unit of time) This Kanban line meet a customer demand rate of 8 (units/unit time) with a minimum WIP inventory cost of 7 ($), initial production quantity at the beginning of production period of 98 units and initial demand quantity at the beginning of production period of units Results and Conclusions 6 kanban = variable []=k drate =variable []=w = 8 prate = variable []=p 000 pratet0 = variable [4]=P0 fprdemand = variable [5]=DF =500 000 dratet0 = variable [6]=D0 rmorderedqty = variable [7]=QR =500 W P Figure [6]: Variables of Fitness Function Limits Kanban technique has been used to achieve the goal of JIT philosophy in a production environment to maximize the profit and to improve the effectiveness of production line The research issue considered here is one of the main problems faced Kanban-based production lines that operate in an industry which are very susceptible to market competition Consequences obtained in this study endow with insights to the application of GA to the Kanban numbers and sizes determination problems and the utility of this technique to industry It has been shown that the application is realizable and able to create a sound results with a satisfactory quality of generated solution 8
The European Applied Business Research Conference Rothenburg, Germany 00 k p 8 8 6 6 4 5 0 5 0 5 0 Iter No 4 0 5 0 5 0 5 0 Iter No (A):Change of Kanban P0 950 900 850 800 5 0 5 0 5 0 Iter No (B):Change of production rate D0 50 5 00 75 50 5 Iter No 5 0 5 0 5 0 (C):Change of initial production (D):Change of initial demand Fitness 6 4 The values for k,w,p,p0,df,d0 that give the minimum fitness noted that DF=QR are: {,8,9,98,500,,500} respectively With a corresponding cost value (Sc)=7 for the following chromosome structure: {{{0000000000000000000000}=k} {{00000000000000000000}=w} {{0000000000000000000}=p} {{00000000000000000}=P0} {{00000000000000000000}=DF} {{0000000000000000000000000}=D0}} 0 5 0 5 0 5 0 Iter No (E): Fitness function Behavior Figure [7]: Designing and Operating of Kanban System The approach has been presented with the goal of being used between two adacent workstations in Kanban based production lines As expected, better results were obtained as the population size, offspring solutions numbers and chromosomes string length are increase till they reach a certain values equal to popsize, maxgen=max and ngp0 respectively Any further increase of these values will not recover the fitness function, but that will affect negatively the time needed to run the system It has been demonstrated how we got the best values of 7,00,4 respectively for such previous constraints 9
The European Applied Business Research Conference Rothenburg, Germany 00 Summary This paper has demonstrated the use of genetic algorithm to design and plan a Kanban based production line for two adacent workstations to attain a minimum cost of WIP between the two stations A generic production and inventory control system has been outlined An optimisation procedure was then used to determine how the system parameters would change with cost structure model within Kanban production line Considerations are on weight of WIP carrying costs The effect of each of the system parameters has been analysed Five system parameters of a production line have been described and quantified The first system of these is represented by variables of fitness function and production model, which ensures that WIP is kept to a minimum Variables of fitness function are; Kanban number, production rate, demand rate, initial production quantity, initial demand quantity, finished product demand, and raw material ordered quantity The second is population size The others are offspring solutions number, string length of chromosomes, and finally probability of mutation and crossover These ensure that the WIP costs are kept to a minimum in response to random patterns of demand and the system is robust to generate a variety of improved solutions in a proper time The fitness variables values have been verified for all types of previous solutions using Microsoft excel and the same results were obtained as shown in the Figure [8] bellow k w P0 p DF D0 QR Sc 89 58 6 58 0 7 589 4 58 5 58 996 5 0 000 5 0 8 98 9 500 500 7 Figure [8]: Verification of Results With Microsoft Excel Recommendation for Future Research The present study is limited to only two adacent Kanban-based workstations It may be worthwhile to focus future research on the design of such approach for multi stages multi-products serial production line References Bhaba R Sarker, and Chidambaram V Balan 999 Operations planning for a multi-stage Kanban system European Journal of Operational Research, Volume, Issue, Pages 84-0 Darrell Whitley December 00 An overview of evolutionary algorithms: practical issues and common pitfalls Information and software technology Volume 4, Issue 4, Pages 87-8 David EGoldberg 989 Genetic Algorithms in Search, Optimization, and machine learning Addison- Wesley 4 John D, Royce O, and John M, August 996, Using Evolution Strategies and Simulation to Optimise a Pull Production System, Journal of Materials Processing Technology, Volume 6, Issue -, Pages 47-5 5 M Al-Tahat, A Liverani, and F Persiani 00,Sept-7 Computer Supported Design Of Flexible Kanban System For A Multi-Stage Multi-Product Manufacturing First International Industrial Engineering Conference Amman-Jordan 6 M Al-Tahat, A Liverani, F Persiani 00,Sept5-7 Design of a Marketing Decision Support Method for Price Quotation In a multi-stage Single-product Kanban-Based Manufacturing Environment XII ADM International Conference - Rimini Italy 7 S M Disney, M M Naim and D R Towill December 000 Genetic algorithm optimisation of a class of inventory control systems International ournals of production economics Volume 68, Issue, Pages 59-78 0