Mathematical Problems in Engineering

Volume 2017, Article ID 8763101, 8 pages

https://doi.org/10.1155/2017/8763101

## Determining Optimal Replacement Policy with an Availability Constraint via Genetic Algorithms

^{1}School of Management, Lanzhou University, Lanzhou, Gansu 730000, China^{2}College of Economics and Management, Lanzhou Institute of Technology, Lanzhou 730050, China

Correspondence should be addressed to Shengliang Zong; moc.liamg@30lhsgnoz

Received 22 April 2016; Revised 23 November 2016; Accepted 15 December 2016; Published 18 January 2017

Academic Editor: Marco Mussetta

Copyright © 2017 Shengliang Zong 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.

#### Abstract

We develop a model and a genetic algorithm for determining an optimal replacement policy for power equipment subject to Poisson shocks. If the time interval of two consecutive shocks is less than a threshold value, the failed equipment can be repaired. We assume that the operating time after repair is stochastically nonincreasing and the repair time is exponentially distributed with a geometric increasing mean. Our objective is to minimize the expected average cost under an availability requirement. Based on this average cost function, we propose the genetic algorithm to locate the optimal replacement policy* N* to minimize the average cost rate. The results show that the GA is effective and efficient in finding the optimal solutions. The availability of equipment has significance effect on the optimal replacement policy. Many practical systems fit the model developed in this paper.

#### 1. Introduction

Power equipment such as a wind power generation system is operating in a dynamic and uncertain environment. Such equipment may fail due to random events such as sand/rain storms, ice rain, heavy snow, and extreme temperatures. The service interruptions caused by the equipment failure will affect the company’s performance in terms of revenue and customer satisfaction. For example, during January 10 to February 5, 2008, an unexpected strong East Asian winter monsoon occurred with extremely low temperatures, blizzard conditions, and freezing rain in southern China. Due to its unusual persistence and intensity, this extreme event caused extensive damage and general disruption over southern China. Most of the infrastructure damages were attributed to freezing rain, which led to a huge amount of broken power equipment, electric power transmission lines, and chaotic traffic conditions for southern China. The disastrous heavy snow caused direct economic losses of 53.8 billion RMB [1, 2].

In practice, the power equipment may be interrupted by random shocks, which represent severe weather conditions such as thunder storm and freezing rain or operational contingencies such as short circuit. When a shock arrives, it will cause a random amount of damage to the system. When the accumulated amount of shock damage exceeds a specified threshold, the system will fail. Therefore, it is necessary to consider the maintenance problem under random shocks. Chen et al. [3] investigated the impact of the voltage imbalance on the reliability of power distribution systems. Kiureghian [4] considered the problem of distribution of electric substation equipment suffering from earthquake.

Recently, shock models were utilized to model the operating time. Conceptually, a system fails due to the shock effect on the system. While most of existing shock models were based on the accumulated or extreme damage causing a system failure. Li [5] first introduced the *δ*-shock model to avoid measuring amount of damage which may not be easy in many situations. Although the shock’s damage amount is difficult to measure, its arrival time can be accurately recorded. Therefore, the *δ*-shock model focuses mainly on the frequency of shocks rather than their accumulated damage amount of shocks. In a *δ*-shock model, a shock causes a system failure which also called a deadly shock if the time interval between two successive shocks is smaller than a critical value *δ*. The threshold *δ* is usually a constant. In this paper, we adopt a general *δ*-shock model by letting *δ* be an exponentially distributed random variable with parameter varying with number of repairs [6]. The power equipment after repair can be more fragile and more prone to fail again. As a result, the threshold of a deadly shock is increasing with* n*, which is the number of repairs.

There are extensive studies on the maintenance problems for the system with operating and repair stages. This is mainly because some classical assumptions are not realistic in modeling the real systems. These assumptions include repaired system becoming “as good as new” and a dialed system being replaced by a new one immediately. Barlow and Proschan [7] introduced an imperfect repair model, where the repair is prefect with probability and minimal with probability . We also consider the nonzero repair and replacement times in contrast to a classical assumption that repairs are instantaneous. In most practical situations, to reflect the aging process of the system, the consecutive repair times are assumed to become longer and longer till the system is replaced with a new one according to some replacement rule. The geometric processes (GP) method, with its practical advantages in characterizing the whole operation process of the components in a system [8], gains worldwide scholar’s attention in recent years, and many researches have been carried out based on GP in different application domains. Lam [9] first introduced GP to study the maintenance problem for the deteriorating system. Finkelstein [10] generalized Lam’s work by scale transformation. Most of the studies based on GP denote a good applicability in the related fields, but their assumption that the property of the component after premaintenance is as good as new is not fully compatible with the electrical power equipment power [11].

Ensuring system function and life are two very important maintenance goals [12]. Li et al. [13] designed an age-related dynamical model to create a timely maintenance schedule that extends the life of the equipment and improves its efficiency. Preventive maintenance also needs to consider how to reduce costs [14]. Castanier et al. [15] proposed a parametric maintenance decision framework to coordinate inspection and replacement of the two components and to minimize the long-run maintenance cost of the system. Lin and Huang [16] used a nonhomogeneous Poisson process with a power law failure intensity to describe the deterioration of a repairable system and to minimize the expected total cost per unit time to determine the optimal nonperiodic maintenance schedule. As we all know, the goal of the maintenance decision is not only to minimize costs. We adopt the replacement policy by satisfying certain requirements of equipment availability. Farag and Al-Baiyat [17] consider two major problems of power system design optimization. One is to maximize the reliability under the cost constraint and the other is to minimize the cost under the reliability constraint. Levitin and Lisnianski [18] studied the optimization problem of the power system containing a bridge substructure in which redundant elements are included in order to provide a desired level of reliability. There are more research works on power system structure optimization under reliability constraints, for example, Lisnianski et al. [19] and Bevrani et al. [20].

Existing research goal is to minimize the total cost or the long-run maintenance cost of the system, and the objective function mainly is linear. In this paper, the objective function is an integer nonlinear programming, and constraints are an integer nonlinear inequality. In order to try to solve this kind of problem in an appropriate way, a method based on genetic algorithms (GA) is proposed to explore the optimality of this problem. GA is heuristic search techniques, which are full of promise to deal with complex optimization problems [21]. GA can also solve optimization problems possessing arbitrary degrees of nonlinearity, discontinuity, and stochasticity. Levitin and Lisnianski [22] used GA to minimize the sum of costs of system modernization actions over the study period while satisfying reliability constraints at each stage. Abdulwhab et al. [23] used GA optimization technique to maximize the overall system reliability for a specified future time period in which a number of generating units are to be removed from service for preventive maintenance.

In this paper, we propose a *δ*-shock model, which will be helpful for practitioners to find the appropriate machine maintenance and replacement policies, for electrical power equipment. We assume that the shocks will arrive according to a renewal process with an exponential distributed interval time of shocks and that the threshold of a deadly shock forms a geometrical increase. Furthermore, the distribution function of the threshold follows exponential distribution. The assumptions and analysis of the model are introduced in Section 2. Section 3 presents an optimization model which is to minimize the long-run average cost per unit time subject to some availability constraints. The GA is introduced in Section 4. Finally, Section 5 provides a case study using our approach.

#### 2. Model and Analysis

##### 2.1. Model Assumptions

To construct the *δ*-shock model, we first introduce the definition of GP as follows.

*Definition 1. *A stochastic process is called a geometric increasing (or decreasing) process, if there is a real number , such that forms a renewal process. is called the ratio of the GP.

Moreover, let and . Then, we haveTherefore, , , and are three parameters for a GP.

We then propose a repair model for electrical power equipment under availability constraints based on the following assumptions.

(I) At the beginning, a new system is installed. Whenever the system fails, it will be immediately either repaired or replaced by an identical system.

(II) The system in operation is subject to environment shocks which arrive according to a Poisson process with rate or . Let be the time interval between the ()th and the th shocks. Then is assumed to be exponentially distributed with the distribution function denoted by . Let be another exponentially distributed random variable associated with . We assume that the sequence forms an increasing geometric process with . Then, has cumulative distribution function , where is the cumulative distribution function of which is a brand new system’s threshold of the shock model. Assume that . follows a -shock model if the system fails at th shock which satisfies , and then the life time or equivalently the operating time is the sum of all until the one satisfying the above condition. Further, we assume that is independent of .

(III) Let be the operating time between the ()th and th repairs. is a sequence of stochastically decreasing random variables.

(IV) Let denote the repair time after the th failure. Then the distribution function of is , where is a constant. Let denote the expected repair time after the th failure. Assume that .

(V) and , are independent.

(VI) The reward and cost structure is assumed as follows: the reward rate when the system is operating is . The repair cost rate of the system is . The replacement cost consists of a fixed replacement cost and variable cost , where is the variable cost rate and is the replacement time with a mean of .

(VII) The -replacement policy is adopted. Under such a policy, the system will be replaced with a new and identical one at the th failure.

(VIII) The equipment must meet the minimum availability requirements. Let represent the availability of the equipment. In other words, the proportion of the nonproductive time, for example, repair or replacement, should not be too high.

##### 2.2. Maintenance Cost

Since , forms an increasing GP with ratio , we have

Thus, the maintenance cost function is shown as

##### 2.3. Operating Revenue

Let be the expected operating time of the system between the ()th and th repairs, and let be the interval time between the ()th and th shock following the ()th repair. Then the number of shocks experienced until the next failure can be written as

We have

Here follows a geometric distribution. That iswhere (a shock following the ()th repair is a deadly shock) and

Therefore, we have . As is a stopping time with respect to the random sequence , which are independent identically distributed random variables. According to Wald’s equation [24], we have

Based on Assumption (II), and the fact that and are exponential distributions, we have

Furthermore, as and are independent and have the marginal exponential distributions with means of and , therefore, we obtain

Then

Consequently,

Then, the operating revenue of the system, denoted by OR, is given by

##### 2.4. Replacement Cost

Based on Assumption (VI), the replacement cost function, denoted by RC, can be written as

##### 2.5. Long-Run Average Cost

Let denote the long-run average cost per unit time of the system. According to renewal reward theorem [24]

Let be the expected length of a renewal cycle of the system under replacement policy . Thus, we have

Let TC be the total expected cost of the renewal cycle of the system under replacement policy . Then,

Consequently,

#### 3. Optimization

##### 3.1. System Availability

Using (13) and (17), the availability function of the system can be written as

##### 3.2. Minimizing the Average Cost

The optimal policy can be found by minimizing the cost of each unit time subject to a minimum level of system availability. The optimization problem is

Thus, we can see that the optimization problem is an integer nonlinear programming problem. The analytical solution is not effective to solve this optimization problem, because the objective function is a nonconvex and extremely complex function. Several classical computational techniques, for example, branch and bound technique, cutting planes technique, implicit enumeration, and out approximation, which are reasonably efficient, have been proposed in literature for solving integer nonlinear programming problems [25–27]. These techniques are applicable to a particular class of problem. In the case of nonconvex problems these techniques may cut-off the global optimal.

#### 4. Genetic Algorithms

In the last two decades, a number of stochastic algorithms have been proposed for solving the integer nonlinear programming problems. Genetic algorithms have been proved to be a useful tool in solving global optimization problems [28]. GA is general purpose population based stochastic search techniques which mimic the principles of natural selection and genetics. It can be implemented easily with the computer programming. In particular, GA is more effective to solve the large scale, real life, discrete, and continuous optimization problems without making unrealistic assumptions and approximations. Recently, GA has been proposed to solve integer programming problems by Soltani and Shafiei [29] and Wong et al. [30]. Applying GA to solve nonlinear mixed integer optimization problems was reported by Pal et al. [31] and Deep et al. [32].

Our objective function is an integer nonlinear programming and constraints are an integer nonlinear inequality. For solving such a problem, we have developed a real coded GA for integer variables.

It is generally accepted that any GA consists of the following basic components.

*Step 1 (set the parameters). *The parameters of the GA in our problem include population size, probability of crossover (pc), probability of mutation (pm), and stopping conditions. The population size should be neither too large nor too small. We set the parameters as follows:

Population size = 10, pc = 0.90, pm = 0.10, stall generations as the stopping conditions.

*Step 2 (generate initial populations). *An initial population of size is randomly generated. Any general procedure to get them can be applied. Nevertheless, in the implementation, we prefer to build the initial population with feasible chromosomes. This device guarantees that, upon termination, the method will provide a bilevel feasible solution. A positive integer is randomly generated according to the nature of the problem. It follows the next law (min (1) + (max (1) − min (1)) rand (Size, 1)).

*Step 3 (calculate the fitness values). *Each chromosome in the population is evaluated by creating the estimated value for each coefficient in (19). The fitness value of each chromosome is then obtained from (20). Let, be the fitness value for each chromosome and the total fitness value of the populations is , where = .

*Step 4 (perform selection). *The roulette wheel selection mechanism is the selection strategy chosen in our paper. Each time, a single chromosome is selected from the current generation to create a new generation. The selection process is as follows. Each time, a random number is generated from the range . Let be the weight of the th chromosome in the whole populations, where . If , then this chromosome is selected. Otherwise, the th chromosome is selected if . This selection process is continued until the new population has been created. The chromosome which has the higher fitness value tends to be picked out to the next generation.

*Step 5 (perform crossover). *Crossover is the key to the power of genetic algorithms. The purpose of the crossover is to generate offspring which endorses them with the common character of their parents. Calvete et al. [33] proposed two types of crossover methods. One is called “variable-to-variable” crossover and the other is called “basis-to-basis” crossover. Lin [34] used the one-point crossover method, which randomly selects one cut-point and exchanges the right parts of two randomly selected parents in the population to generate offspring. The crossover method of this paper uses also the one-point method.

*Step 6 (perform mutation). *Mutation is an operator that acts on various chromosomes by randomly exchanging one of its indices. Mutation resets a selected position in a chromosome to a randomly generated real number in . The mutation rate is usually a very small value; on average, 0.3% of the total populations undergo mutation [34]. The mutation method used here is the midvalue method, where the mutation of chromosome took place in the middle of the population. The new chromosome contains the indices of the new basis.

*Step 7 (return to Step 3 until the stopping condition is met). *There are many stopping conditions for GA, including generations, time limit, fitness limit, stall generations, and stall time. The programs will stop until one of the stopping conditions is met. The stopping condition used here is stall generations.

#### 5. Empirical Results

In this section we analyze an example using the proposed GA. All experiments were performed in MATLAB on a PC with an Intel Core i7 2.80 GHz and 4 G RAM. The data is simulation data which are used to further evaluate the proposed approaches. The simulation data are generated from the analytical curves. Table 1 shows the value of the parameters.