2 Level of Repair Analysis Based On Genetic Algorithm With Tabu Search, 2010
2 Level of Repair Analysis Based On Genetic Algorithm With Tabu Search, 2010
net/publication/45534747
CITATIONS READS
3 222
3 authors, including:
All content following this page was uploaded by Mohammed Kishk on 30 May 2014.
Abstract - Genetic algorithms and their hybrid schemes have Since these support resources are very costly for systems
shown a great efficacy in solving large scale combinatorial containing thousands of components, there is always a trade-
problems in which solutions are highly time-consuming. The off between availability and maintenance costs. Ignoring these
level of repair analysis (LORA), mathematically formulised support issues in the early phase, however, may have costly
by an integer programming model (IP), is very difficult to consequences during the operation phase. This becomes in
optimize by means of traditional optimization techniques due increasingly evident with fierce market conditions. As a
to a large number of decision variables involved. In this paper, result, companies are compelled to move from fragmented
a hybridised Genetic Algorithm with Tabu Search is presented decision processes toward more integrated acquisition and
and its application to solve Level of repair analysis (LORA) operation decision in order to sustain their equipment use at
problem is investigated. The LORA, considered as an low cost and required availability. Consequently, some
important tool for strategic system maintenance decision industries become aware of the large potential for cost
making, seeks to determine the location in the repair network reductions by adopting whole life techniques in their
at which a failed component should be discarded or repaired. acquisition process. Level of repair analysis LORA is one of
The proposed algorithm is developed in order to determine the the prescribed techniques in the military and maritime
best repair decision combination. The efficacy of the industries to achieve a system design with the minimum
algorithm is investigated in the context of a case study. The whole life maintenance cost (MIL.STD.1390D).
maintenance costs of a structure of three-echelon repair and
multi-indenture is optimised under the condition that repair
decision should be taken for all system items. Typical results 1.2. Level Of Repair Analysis (LORA)
have shown that the algorithm can effectively handle a real
industrial sized case study with adequate optimisation The LORA approach was developed by military industry to
computational time. evaluate the Integrated Logistic Support (ILS) factors
contributing to the systems whole life cost. When a failure
Keywords: Level of repair analysis, maintenance occurs, failed components are removed and repaired or
optimisation, Genetic Algorithms, Tabu Search. replaced by new spare parts. In designing systems, Level of
repair analysts considers all aspects of the system design and
1. Introduction maintenance scenarios to achieve availability and cost
balanced systems. As a result, they provide essential support
1.1. Background requirements for the most effective maintenance strategy
Over the last decade considerable emphasis has been put on under predicted operational environments.
whole life costing (Tysseland, 2007; Kleyner & al., 2008;
Lindholm et al., 2004; Kishk et al., 2003). Traditional system The basic of the LORA process is the following. Level of
acquisition has taken place on the purchase costs without repair analysts have to decide for a given design which
explicit indication of operation costs. However, the decrease components to repair, which components to discard, where in
of company operational budgets has unveiled the need for the repair network to do this and finally where in the repair
techniques to forecast and to optimise ongoing cash flows network to install the required maintenance resources. Thus, a
over the system whole life. As result, a prerequisite for number of reparation locations in which systems, subsystems
effective acquisition decision is the operational parameters and components have to be repaired or discarded is set up to
such as: system availability, system reliability and the cost of satisfy maintenance requirements at minimum cost.
the required support resources. Besides, it was found that
operational readiness of expensive complex structures such as
petroleum apparatus, aircrafts, ships and military equipments
is very sensitive to the availability of spare components,
maintenance resources and manpower.
1
The Scott Sutherland School, The Robert Gordon University, Garthdee
Road, Aberdeen AB10 7QB, UK.
2
School of Engineering, The Robert Gordon University, Schoolhill, Aberdeen
AB10 1FR, UK.
Figure 1 : A multi-echelon repair network
In the literature, various models of LORA have been discarding it. Level of Repair Analysis (LORA) is an
discussed for a three echelon repair network (figure 2) and approach which examines the cost balance between repairing
multi indenture system (figure 1). These models involve a the component and discarding it. The framework of this
large number of decision variables which makes LORA approach is an iterative process ensuring the optimum
problem very difficult to optimize by means of traditional maintenance planning. However, the LORA problem as
optimization techniques. For instance, the number of all combinatorial optimization is not the most widely studied in
possible combinations (part, repair and discard decision) for a the literature (Basten & al., 2006). A limited works were
system consisting of 32 parts spread between different devoted to solve the LORA problem. Barros et al. (2001),
indentures is 6.28 × 1010 (Kumar & al., 2006). Hence, Saranga et al. (2006), Gutin et al. (2005) and Basten et al.
techniques like integer programming and branch and bound (2009) modelled LORA as Integer Programming model in
method become difficult to use. Consequently, we focus on which all repair locations at the same echelon were
the opportunity of using Genetic Algorithms which are the aggregated. Besides, they all resolved their model under
most suitable to problems involving combinatorial infinite capacity of resources. Brick et al. (2009) model
optimisations. LORA without aggregating data per echelon level for only 1
echelon and 2 indenture levels.
Component i is the parent of the component j or component algorithms to enhance their capabilities in solving
j is the child of the component i combinatorial optimisation (Zdanski & al., 2002 and
r = repair options: repair, discard or move. Hagemana & al., 2003). For instance, a GA speed is low for
i = Total number of maintenance tasks required in the the huge size population and TS relies strongly on the initial
whole life time of component i. solution. Consequently, GA and TS combination named
CFr,e,i = fixed cost related to repair option ‘r ’ at echelon e, GATS may overcome these limitations and maintain their
for component i. advantages.
VC r,e,i = variable cost related to repair option ‘r’ at echelon
e, for component i.
3.1. Genetic algorithms
Let X, the repair decision, is 1 if repair decision has been
chosen at the echelon e of the selected component i and 0 Genetic algorithms are stochastic search techniques based on
otherwise. the theory of evolution for finding the global optimum
solution. The genetic algorithm developed by Holland to
1 if repair option r at echelon e is selected for part i optimise a function F(x), where x is a vector representing
X r,e,i
0, otherwise individual solutions (Gen & al. 2000). First of all, Genetic
algorithms generate not only a single solution but a group of
Using the notation mentioned above, the total maintenance solutions, called a population. This population changes over
cost is: time, but it always keeps its initial size. The population
members are called strings or chromosomes from which a
VC it FC r,e,i X r,e,i
N 3 m
r,e,i
subset called parents is selected according to the best values
i 1 r 1 e 1 of F(x). A fitness value in Genetic algorithms is a measure of
goodness of a solution to the objective function, i.e., the
(1) fitness of an individual is directly related to its objective
Subject to function value. At any iteration, a fitness value is calculated
X r, 1, i 1 for all parts for each of the current individuals. The selection rule, called a
survivability test, exclude from the population the strings
which have the worst finesses. Second, new solutions called
(2) children (or offspring) are produced by genetic operators:
3 crossover and mutation. Together parents and new children
X move, e, i r 1
X r, e 1, i 1 are grouped in a new population which will pass again
through survival test. Thus, the population as a whole moves
iteratively towards better solutions ideally to the global
(3) optimum.
X r, e, j X r, e, i
commonality matrix (figure 4), where the column represents generated randomly. In addition, we select one of the best
parent items and in the row are child items. We start by solutions and we generate a repair decision for a component
assorting parts from the first indenture until the last but one selected at random. Again, we adjust the new changes
indenture in the column as parent items. Then, we insert parts according to the constraint Eq. (4).
from the second indenture to the last one in the commonality
matrix row. As shown, child parts 5, 6 and 7 belong to parent In our GATS algorithm, these two operators are applied for
part 3 or parent part 3 is constituted of child parts 5, 6 and 7. the individual generated by Genetic Algorithm and improved
According to this representation, whenever the parent part 3 is by tabu search.
under discard or move decision, the child parts 5, 6 and 7 will
have the same decision (constraint Eq. 4).
The procedure of TS consists of the following steps which are 3. Repeat step 2 until there is no improvement of the best
depicted in the figure (6). First, a number of neighbourhood fitness value.
solutions that can be produced from an initial solution are 4. Replace the solution by its best neighbourhood.
examined. Then, a solution with the best fitness value and it is 5. Choose two solutions to produce new chromosomes
not in the tabu list is selected from the explored using genetic operators: parent selection and crossover.
neighbourhood. This way, tabu search tries to assure that the These new solution are accepted when they verify the
method does not re-examine a solution previously generated. constraint equations 2, 3 and 4.
Finally, TS procedure iterates the previous step until no more 6. Create a new chromosomes using genetic operator:
neighbours are present (all are tabu), or when during a mutation.
predetermined number of iterations no improvements are 7. Update a tabu list of the best chromosomes.
found. 8. Repeat step 1 until there is no improvement of the best
chromosome.
4. Computational experience
3.3. Genetic & Tabu Search Algorithm Table 1: Case study data used in Saranga’s work, maintenance
costs
This approach, widely used in the literature, combines the
advantages and mitigates the disadvantages of the two
algorithms. Tabu search relies only on one solution and miss
information of a larger set of solutions, however, Genetic
Algorithms lead to lower solution quality with increasing
problem size (Zdanski & al., 2002). In this study the GATS
algorithm starts by generating N initial possible solutions
(figure 6). A tabu search, as an iterative process, is then used
for upgrading theses solutions through neighbouring
exploration. Afterwards, the flow returns to the Genetic
Algorithm which is again an iterative process. By means of
the genetic operators new off springs are produced. Then, a
tabu list of the best solutions is updated by the new off springs
according to the fitness value. The stopping criteria for the
GATS algorithm are a predefined number of consecutive
iterations attaining the same best solution is reached.
The mean steps of the algorithm are shown in figure 6 and Table 2 : Case study data used in Saranga’s work, material
described as follows: structure
1. Generate randomly a set of solutions (20 solutions)
verifying the equations 2, 3 and 4.
2. Refine each solution by the neighbourhood routine with
respect to fitness value. A neighbourhood solution is
obtained only by modifying the value of one element
from the solution under consideration to 1 or 0. Besides,
the neighbourhood solutions are not accepted until they
verify the constraint equations 2, 3 and 4. Then, a tabu
list is updated containing all the fitness values of the In this experiment, the settings are chosen as described in
solutions that have been explored. After, a new (Saranga & al., 2006) on two echelon repair network for an
neighbourhood is explored only when its fitness value aircraft engine with three-indenture structure. Different costs
does not exist to the tabu list. (CF: fixed costs and CV: variable costs) for various repair
options at different echelons for all items are listed in Table A comparison between two and three echelon network
(1). Besides, the commonality matrix shows the relationship computational time that takes GATS algorithm to come out
between first indenture parts (from 1 to 10) and the second with the optimal solution is shown in (figure 7). The
indenture parts (from 11 to 32). The optimal or near optimal computational time increases exponentially with system
solution obtained by Saranga’ work and our GATS algorithm structure size and the bigger the number of echelon is the
are similar, only part 5 got different repair decision (table 3). higher the computing time is. Thus, researchers consider
The total maintenance costs incurred are respectively three echelon repair network is enough in practice to handle
4255.274 and 4216.274. In addition, the maximum maintenance activities and to be modelled by acceptable
maintenance cost achieve by the simulation is 42 759 computational time.
representing more than ten times the optimum solution. This
witnesses the benefit of adopting LORA in maintenance
planning and support provision.