Next Article in Journal
Single-Machine Maintenance Activity Scheduling with Convex Resource Constraints and Learning Effects
Next Article in Special Issue
Gibbs Distribution and the Repairman Problem
Previous Article in Journal
On the Equilibrium in a Queuing System with Retrials and Strategic Arrivals
Previous Article in Special Issue
Event-Triggered Extended Dissipativity Fuzzy Filter Design for Nonlinear Markovian Switching Systems against Deception Attacks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Redundancy Allocation of Components with Time-Dependent Failure Rates

1
Center of Research on Risk and Crises (CRC), MINES-Paris, PSL University, 75272 Paris, France
2
Energy Department, Politecnico di Milano, 20133 Milan, Italy
3
Department of Industrial Engineering, University of Gonabad, Gonabad 96919-57678, Iran
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(16), 3534; https://doi.org/10.3390/math11163534
Submission received: 3 July 2023 / Revised: 11 August 2023 / Accepted: 12 August 2023 / Published: 16 August 2023

Abstract

:
The Redundancy Allocation Problem (RAP) is well-known in the field of reliability optimization. In this paper, RAP is investigated assuming that the distribution of the time to failure of the components has the form of an Erlang distribution with a time-dependent rate parameter and considering that the choice of redundancy for each subsystem can be none, active, standby or mixed. A genetic algorithm is used to solve the problem of optimal allocation. To analyze the effect of the time dependence, some numerical examples are worked out. Then, a case study of RAP from the literature is analyzed. The obtained results show that time dependence of the failure time distribution parameters can lead to significant differences in the optimal redundancy allocation.

1. Introduction

The reliability of a system can be improved at design by allocating redundant components to the subsystems. In a redundant system, when a working component fails, the failed component function is taken over by a redundant one which can provide the same function. The problem to do this optimally is called Redundancy Allocation Problem (RAP), which amounts to optimally selecting the redundant components to be allocated to the subsystems. The choice of the redundancy allocation strategy is driven by objectives, and/or subject to constraints, such as reliability, cost, volume and weight, depending on the system application. Then, RAP can be defined as the problem of selecting the number and type of the redundant components to allocate to the system, so to optimize a specific objective, such as system reliability, under given constraints, such as total system cost and weight.
The two redundancy strategies commonly considered in RAPs are “active” and “standby”, although one called “mixed” has been introduced relatively recently [1]. In the active redundancy strategy, all components start to work at time zero and a redundant component might be found already failed before it is called to operation. In the standby redundancy strategy, one component is active and the others (at least one component) are standby. The standby redundancy strategy comes in the three types: cold, warm and hot [2]. In cold-standby, redundant components do not fail before they are used. In warm-standby, redundant components are subject to failure even when they are dormant, but at a rate lower than that of the main component which is working. Finally, in hot-standby, redundant components are subject to failure even while dormant, at a rate equal to that of the main components. Without loss of generality of the analysis, among the standby strategies, the present paper considers the cold-standby redundancy strategy, only for simplicity of illustration. In the mixed redundancy strategy, more than one component can be active while some others (one or more) are standby. The mixed redundancy strategy tries to combine the advantages of both active and standby strategies. Indeed, it is a general form of the active and standby strategies. The mixed redundancy strategy can be used in most practical cases, but the active and standby redundancy strategies are also used in certain cases, due to their simple forms. In the standby and mixed redundancy strategies, a switching system has to be used to activate a standby redundant component.
Some researchers considered RAP with the same redundancy strategy for all subsystems in the system [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34]. Following Coit [35], some researchers have turned to consider the redundancy strategy for each subsystem, either active or standby, as a decision variable that can be optimally determined by an optimization model [36,37,38,39,40,41,42,43]. The mixed redundancy strategy has also been considered as an alternative for the subsystems [44,45,46,47,48,49,50,51,52].
Other aspects come realistically into play in a RAP; these include the objective function, the states of the components and that of the entire system, the type of the components to be allocated to each subsystem, the reparability of the components, and so on. For example, RAPs may be considered with a single objective function, e.g., aiming to maximizing system reliability [1,3,4,5,6,7,8,9,10,11,27,28,30,31,35,36,37,39,40,41,44,46,47,48,53,54,55], minimizing total system cost [12,13,14,15] or maximizing/minimizing other parameters [16,29]. Alternatively, RAPs may consider multi-objective functions [18,19,20,21,22,23,24,25,32,34,38,42,45,50,56,57,58,59].
Regarding the components, they might be represented as binary or multi-state. In the binary state representation, the components can only be totally healthy or completely failed [1,3,4,5,6,7,8,9,10,11,12,13,17,18,19,20,21,22,23,25,26,27,28,29,30,31,32,33,34,35,40,41,43,44,45,46,47,48,53,55]; in the multi-state, the components might have other states, intermediate between these two [14,15,16,24,60,61,62,63,64,65]. The type of the components in the subsystems can be characterized from different viewpoints. For instance, allocation of either identical [1,3,4,5,6,7,12,18,19,20,35,36,40,45,46,53,54] or non-identical components to each subsystem might be allowed [8,9,10,11,13,14,15,16,17,21,22,25,27,28,33,37,41,44,48,49,62]. Moreover, the components characteristics, including their cost, weight and performance function, might be deterministic or uncertain functions [66]. From another aspect, the components can be repairable [14,15,20,24,26,30,56] or non-repairable [1,3,4,5,6,7,9,11,13,16,17,18,19,21,22,27,28,29,31,35,36,37,40,41,43,44,45,46,47,48,53,54,55]. In this regard, a comprehensive categorization of reliability optimization problems is provided in [67].
This paper considers the RAP for a system with identical binary components in each subsystem. The redundancy strategy for each subsystem is taken as the decision variable that can be active, standby or mixed, and the objective function requires system reliability to be maximized. Also, the distribution of the time to failure of the components is assumed to have the form of an Erlang distribution, whose rate parameter (λ(t)) is not constant in time. To the best of our knowledge, no RAP with these novel realistic properties has ever been investigated in the literature.
The failure of the components allocated to a system may be due to degradation or due to occurrence of some shocks. In the case where component failure is due to shock occurrence, the Erlang distribution can be considered suitable for modeling the stochastic time to failure of the components. If the distribution of the time to failure of a component is considered as an Erlang distribution with a shape parameter of k, it means that the failure of the component is due to the occurrence of some shocks, so that the component will fail as soon as the kth shock occurs. So, the problem is of practical interest, for example for electronic components which often fail due to some shocks. In the Erlang distribution, the rate parameter indicates the rate of occurrence of the shocks. The rate parameter is usually assumed to be constant, to describe the failure of components subject to shocks occurring at a constant rate. But the occurrence rate of the shocks may not be constant in time and, for this reason, in this paper the rate parameter is considered time-dependent, λ(t).
The main contribution of this paper is to propose a model in which the rate parameter of the time to failure distribution of the components is time-dependent and the redundancy allocation strategy for the subsystems can be mixed or other redundancy strategies.
The rest of the paper is structured as follows. In Section 2, the mathematical formulation of the RAP is described in detail. A solution method based on a genetic algorithm is proposed in Section 3. A numerical example is elaborated in Section 4 with and without time dependence, to demonstrate the effect of the time dependence of the rate parameter. Then, a case study taken from literature is considered. Finally, the conclusions and the suggestions for future work are given in Section 5.

2. Materials and Methods

2.1. Problem Formulation

A single-objective RAP is considered, in which the distribution of the time to failure of any component has the form of an Erlang distribution. The main specificities are the time-dependent parameter λ(t) in the time to failure distribution function and the redundancy strategy of any subsystem that can select among active, standby, mixed or no redundancy strategies. Figure 1 depicts the structure of a system composed of a series of five subsystems.
A number of redundant components can be allocated to each subsystem, each of which may be active or standby. In Figure 1, the active components are shown in green and the standby ones are shown in grey. The purpose of allocating the redundant components is to increase the reliability of the subsystems and consequently to increase the reliability of the system. But, the existence of some physical and/or economical restrictions makes it impossible to allocate any desired number of components to the subsystems. Therefore, the goal is to determine the best choices including determining the type and number of active and standby components allocated to each subsystem so that the reliability of the system is maximized, under given constraints.
The number of active and standby components allocated to each subsystem determines the redundancy strategy. If no redundant component is allocated to a subsystem, the subsystem has no redundancy strategy. The redundancy strategy of a subsystem is active if all of the components allocated to the subsystem are active. If one of the components allocated to a subsystem is active and the rest are standby, the redundancy strategy is standby, and if more than one allocated component is active and one or more components are standby, the redundancy strategy is mixed. Since the number of active and standby components allocated to each subsystem is a decision variable, the redundancy strategy is also a decision variable that is determined during the optimization process.
In Figure 1, the redundancy strategy for each subsystem can be determined according to the number of active and standby components allocated to the subsystem. Accordingly, the redundancy strategy for subsystems 1 and 4 is standby, that for subsystems 2 and 5 is mixed, and the one for subsystem 3 is active. Before addressing the mathematical model, the assumptions considered and the notations used are presented below.

2.1.1. Notation

i index   for   subsystems   ( i = 1 , 2 , , s )
s number of subsystems
n A , i number of active components allocated to subsystem i
n S , i number of standby components allocated to subsystem i
n A ,   n S vectors   of   [ n A , i ]   and   [ n S , i ] ,   ( n A = n A , 1 , n A , 2 , , n A , s ,   n S = n S , 1 , n S , 2 , , n S , s )
n i number   of   components   allocated   to   subsystem   i ,   ( n i = n A , i + n S , i )
n m a x , i maximum   value   allowed   for   n i
A R S i redundancy   strategy   assigned   for   subsystem   i ,   A R S i N R , A c t i v e , S t a n d b y , M i x e d
r i ( t ) reliability at time t of the component (all identical) allocated to subsystem i
R s y s ( t ; n A , n S ) system   reliability   at   time   t   due   to   the   vectors   n A   and   n S
R i , A R S i ( t ) reliability   of   subsystem   i   at   time   t   with   redundancy   strategy A R S i
k i shape parameter of the Erlang distribution for the component allocated to subsystem i
λ i ( t ) rate parameter of the Erlang distribution for the component allocated to subsystem i
c i , w i cost and weight of the component allocated to subsystem i
c s w i t c h , i , w s w i t c h , i cost and weight of the switching system used in subsystem i (if present)
W , C maximum allowed amount for the system weight and cost
ρ i ( t ) switching system reliability at time t for subsystem i (if present)
f i ( j ) ( t ) pdf of the jth failure time of standby components allocated to subsystem i at time t
f i M a x , n A , i ( t ) pdf of time to failure of the last active component allocated to subsystem i at time t
f i ( t ) pdf of the time to failure of the component allocated to subsystem i at time t
F i ( t ) cdf of the time to failure of the component allocated to subsystem i at time t
λ a b average rate of occurrence of events between times a and b ( λ a b = a b λ ( t ) d t )
N ( t ) number of events occurred until time t.

2.1.2. Assumptions

As previously mentioned, the objective is to maximize the reliability of a system under a given number of constraints. The following assumptions are made:
  • The components and the entire system are binary, i.e., they can be either totally healthy or completely failed.
  • Components failures are independent events that individually cause no damage to the system.
  • The redundancy strategy for each subsystem is a decision variable that can be selected to be active, standby, mixed or none.
  • The distribution of the time to failure of the components has the form of an Erlang distribution with a time-dependent parameter λ(t).
  • The components are non-repairable and there is no preventive maintenance.
  • The components of each subsystem are identical, i.e., mixing of component types is not allowed.

2.1.3. Mathematical Model

The mathematical model for the considered problem can be formulated as follows:
max   R s y s ( t ; n A , n S )
s . t .
i = 1 s c i . n i + i S t a n d b y o r   M i x e d c s w i t c h , i C
i = 1 s w i . n i + i S t a n d b y o r   M i x e d w s w i t c h , i W
1 n i n m a x , i
Equation (1) expresses the objective of maximizing the system reliability by determining the type and the number of active and standby components allocated to each subsystem of the system. Constraints (2) and (3) express the maximum cost and weight for the entire system, respectively, whereas constraint set (4) considers the maximum number of components that can be allocated to each subsystem. The reliability of a system composed of a series of subsystems can be calculated as in Equation (5):
R s y s ( t ; n A , n S ) = i N R R i , N R ( t ) i A c t i v e R i , A c t i v e ( t ) i S t a n d b y R i , S t a n d b y ( t ) i M i x e d R i , M i x e d ( t )
where NR represents subsystems with no redundancy strategy, whereas Active, Standby and Mixed represent those with active, standby and mixed redundancy strategies, respectively. As previously mentioned, the redundancy strategy for each subsystem is determined according to the number of active and standby components allocated. Therefore, the number of active and standby components allocated to each subsystem determines in which part of Equation (5) the reliability calculation of the subsystem should be placed. After determining the redundancy strategy for each subsystem, the reliability of the subsystems at time t is calculated with one of the following equations:
R i , N R ( t ) = r i ( t )
R i , A c t i v e ( t ) = 1 1 r i ( t ) n i
R i , S t a n d b y ( t ) = r i ( t ) + j = 1 n i 1 0 t ρ i ( u ) f i ( j ) ( u ) r i ( t u ) d u
R i , M i x e d ( t ) = 1 1 r i ( t ) n A , i + 0 t ρ i ( u ) f i M a x , n A , i ( u ) r i ( t u ) d u + j = 1 n S , i 1 0 t v t f i M a x , n A , i ( v ) ρ i ( u ) f i ( j ) ( u v ) r i ( t u ) d u d v
In these equations, ri(t) is the reliability of the components used in subsystem i at time t, ni is the number of components allocated to subsystem i, ρi(t) is the reliability of the switching system at time t, f i ( j ) ( t ) is the pdf of the time to failure of the jth component in subsystem i in standby and mixed redundancy strategies, whereas nA,i and nS,i represent the numbers of active and standby components assigned to subsystem i, respectively, in the case a mixed redundancy strategy is used (in this case nis = nA,i + nS,i), and f i M a x , n A i ( t ) is the pdf of the time to failure of the last active component in subsystem i in the mixed redundancy strategy.
In Equation (6), the reliability of subsystem i with no redundancy strategy is equal to the reliability of the only component allocated. In Equation (7), the reliability of subsystem i at time t with an active redundancy strategy is equal to the probability that at least one of the allocated components can remain healthy until time t. In Equation (8), the reliability of subsystem i at time t with standby redundancy strategy consists of two parts. The first part is the reliability of the active component allocated and the second one is equal to the sum of the probabilities that the jth standby component fails at time u (a time before t) and the next standby component is activated by the switching system at this time and can be healthy until time t. In Equation (9), the reliability of subsystem i at time t with a mixed redundancy strategy consist of three parts. The first part is equal to the probability that at least one of the active components remains healthy until time t. The second part is equal to the probability that the last active component fails at time u (a time before t) and the first standby component is activated by the switching system at this time and can remain healthy until time t. The third part is equal to the sum of the probabilities that the last active component fails at time v (a time before t) and, then, the standby components will be activated one by one after the failure of the previous one, until the jth standby component fails at time u (a time between v and t), and the next standby component is activated by the switching system at this time and can remain healthy until time t.
In Equation (9), f i M a x , n A , i ( t ) can be calculated using Equation (10):
f i M a x , n A , i ( t ) = n A i F i ( t ) n A , i 1 f i ( t )
For example, if n A , i = 3 f i M a x , n A , i ( t ) can be calculated as in Equation (11):
f i M a x , 3 ( t ) = f i ( t ) F i ( t ) F i ( t ) + F i ( t ) f i ( t ) F i ( t ) + F i ( t ) F i ( t ) f i ( t ) = 3 F i ( t ) 3 1 f i ( t )
If the time to failure of a component in subsystem i follows an Erlang distribution, then its reliability at time t is given by Equation (12):
r i ( t ) = t f i ( t ) d t = j = 0 k i 1 e λ i t ( λ i t ) j j !
When the time to failure of a generic component follows an Erlang distribution with shape parameter k, this component is, then, exposed to shocks that follow a Poisson distribution, and fails at the occurrence of the kth shock. The reliability of the component at time t is, therefore, the probability that less than k shocks occur until time t.
Considering that the times to failure of the components follow Erlang distribution, the reliability of subsystem i with no redundancy strategy and with an active redundancy strategy is obtained as Equations (13) and (14), respectively:
R i , N R ( t ) = j = 0 k i 1 e λ i t ( λ i t ) j j !
R i , A c t i v e ( t ) = 1 1 j = 0 k i 1 e λ i t ( λ i t ) j j ! n i
If the redundancy strategy of subsystem i is standby, considering that  ρ i ( t )  is a non-increasing function, we have:
R i , S t a n d b y ( t ) r i ( t ) + j = 1 n i 1 ρ i ( t ) 0 t f i ( j ) ( u ) r i ( t u ) d u = j = 0 k i 1 e λ i t ( λ i t ) j j ! + ρ i ( t ) j = k i n i k i 1 e λ i t ( λ i t ) j j !
The reliability formulation for a subsystem with a mixed redundancy strategy is slightly more complicated than the others. In Equation (9), we can rewrite ri(t) and Fi(t) as in Equation (12). Also, we know that for a component with Erlang time to failure distribution of shape parameter ki, fi(t)dt is the probability that ki − 1 shocks occur before time t and the kith shock occurs at time t. So, we can write fi(t)dt so as to explicitly represent that ki − 1 shocks occur before time t and one shock occurs between time t and t + dt:
f i ( t ) d t = e λ i t ( λ i t ) k i 1 ( k i 1 ) ! . e λ i d t ( λ i d t ) 1 1 !
For f i ( j ) ( t ) , the pdf of the time to failure of the jth component in subsystem i, we have:
f i ( j ) ( t ) d t = e λ i t ( λ i t ) j . k i 1 ( j . k i 1 ) ! . e λ i d t ( λ i d t ) 1 1 !
This means that f i ( j ) ( t ) d t  is the probability that j.ki − 1 shocks occur before time t and one shock occurs between time t and t + dt. So, if the redundancy strategy of subsystem i is mixed, the reliability of subsystem i, whose components times to failure follow an Erlang distribution, can be calculated as in Equation (18):
R i , M i x e d ( t ) = 1 1 j = 0 k i 1 e λ i t ( λ i t ) j j ! n A , i + 0 t ρ i ( u ) n A , i 1 j = 0 k i 1 e λ i u ( λ i u ) j j ! n A , i 1 e λ i u ( λ i u ) k i 1 ( k i 1 ) ! . e λ i d u ( λ i ) 1 1 ! j = 0 k i 1 e λ i ( t u ) ( λ i ( t u ) ) j j ! d u + j = 1 n S , i 1 0 t v t ρ i ( u ) n A , i 1 j = 0 k i 1 e λ i v ( λ i v ) j j ! n A , i 1 e λ i v ( λ i v ) k i 1 ( k i 1 ) ! . e λ i d v ( λ i ) 1 1 ! . e λ i ( u v ) ( λ i ( u v ) ) j . k i 1 ( j . k i 1 ) ! . e λ i d u ( λ i ) 1 1 ! j = 0 k i 1 e λ i ( t u ) ( λ i ( t u ) ) j j ! d v d u
In this paper, it is assumed that the rate parameter is not constant, but a function of time, i.e., we have λ(t) instead of λ. So, the components will be exposed to shocks with a variable rate over time. For a nonhomogeneous Poisson process with rate λ(t), we have:
P N ( b ) N ( a ) = j = e λ a b ( λ a b ) j j !
where
λ a b = a b λ ( t ) d t
and N(t) is the number of events occurring prior to time t. Therefore, the reliability at time t for a component exposed to the shocks that follow a nonhomogeneous Poisson distribution, i.e., the probability that less than k shocks occur until time t, can be calculated as in Equation (21):
r ( t ) = P N ( t ) N ( 0 ) < k = j = 0 k 1 e λ 0 t ( λ 0 t ) j j !
Now, if the rate parameter is time-dependent (λ(t)), by replacing λ i t by λ i , 0 t in Equations (13)–(15) and (18), the reliability of the subsystems with no redundancy strategy, and of those with active, standby and mixed redundancy strategies can be written as Equations (22)–(25), respectively:
R i , N R ( t ) = j = 0 k i 1 e λ i , 0 t ( λ i , 0 t ) j j !
R i , A c t i v e ( t ) = 1 1 j = 0 k i 1 e λ i , 0 t ( λ i , 0 t ) j j ! n i
R i , S t a n d b y ( t ) = j = 0 k i 1 e λ i , 0 t ( λ i , 0 t ) j j ! + ρ i ( t ) j = k i n i k i 1 e λ i , 0 t ( λ i , 0 t ) j j !
R i , M i x e d ( t ) = 1 1 j = 0 k i 1 e λ i , 0 t ( λ i , 0 t ) j j ! n A , i + 0 t ρ i ( u ) n A , i 1 j = 0 k i 1 e λ i , 0 u ( λ i , 0 u ) j j ! n A , i 1 e λ i , 0 u ( λ i , 0 u ) k i 1 ( k i 1 ) ! . e λ i , u u + d u ( λ i ( u ) ) 1 1 ! j = 0 k i 1 e λ i , u t ( λ i , u t ) j j ! d u + j = 1 n S , i 1 0 t v t ρ i ( u ) n A , i 1 j = 0 k i 1 e λ i , 0 v ( λ i , 0 v ) j j ! n A , i 1 e λ i , 0 v ( λ i , 0 v ) k i 1 ( k i 1 ) ! . e λ i , v v + d v ( λ i ( v ) ) 1 1 ! . e λ i , v u ( λ i , v u ) j . k i 1 ( j . k i 1 ) ! . e λ i , u u + d u ( λ i ( u ) ) 1 1 ! j = 0 k i 1 e λ i , u t ( λ i , u t ) j j ! d v d u

2.2. RAP Solution Method

RAP is a problem of NP-hard class [68], and a meta-heuristic algorithm is here developed to search for the optimal solution. The genetic algorithm (GA) is a well-known meta-heuristic method for solving combinatorial optimization problems [69] and, indeed, it has been successfully employed for solving optimization problems in various fields [1,4,5,6,8,24,30,38,41,42,43,44,45,46,47,56]. In the following, the deployment of the algorithm for the solution of the proposed RAP is briefly described.

2.2.1. Solution Encoding (Chromosome)

Solution encoding in the GA used is a 3 × s matrix, in which s represents the number of the subsystems in the system. For each subsystem, the first row is the type of selected component and the second and third rows give the numbers of allocated active and standby components, respectively. There is no need to add a row for the redundancy strategy, as it may be determined using the second and third rows: if the number of standby components is equal to zero, then the redundancy strategy is active (or no redundancy strategy, in case the number of active components is equal to one); if the number of active components is equal to one and there are some standby components, then the redundancy strategy is standby; finally, if the number of active components is larger than one and there are also some standby components, then the redundancy strategy is mixed. Figure 2 shows an example for the solution encoding with s = 6 subsystems. In the Figure, the active components are shown in green and the standby ones are shown in grey. As it can be seen in the Figure, the redundancy strategy for subsystems 1 and 4 is standby, that for subsystems 2 and 6 is mixed, and that for subsystems 3 and 5 is active. Figure 3 represents the structure of the system.

2.2.2. Initial Population

The initial population is composed of the nPop solution matrix, whose entries are generated randomly. The solution matrices should be generated in such a way that in each column (subsystem), the first element is an integer number between 1 and the number of components types available for the subsystem, and the second and third elements are integer numbers so that their summation is less than or equal to the maximum number of components allowed for the subsystem.

2.2.3. Fitness Function

The fitness function is defined as the system reliability minus the penalty for constraints violation. Indeed, if a solution goes beyond the constraints, a large enough amount of penalty is deducted from the fitness function to guarantee that the final optimal solution is feasible. The amount of the penalty is equal to the product of the total amount of violations of the constraints and a fixed number. This penalty also allows the algorithm to search in the infeasible space for maintaining appropriate diversity in the search.

2.2.4. Selection

At each iteration, new offspring are generated by implementation of the crossover and mutation operators. The roulette wheel selection method is used to select the parents [8,40]. In the roulette wheel selection method, selection of the parents occurs based on their fitness function value: in other words, a solution is selected with probability proportional to its fitness function. After the selection process, the crossover and mutation operators are performed.

2.2.5. Crossover Operator

The crossover operator activates at a prespecified rate, rc. To explore a large variety of solutions, use has been made of different types of single-point, double-point, uniform and max–min crossover operators [1,4,44], each with a predetermined probability.

2.2.6. Mutation Operator

The mutation operator is implemented at a prespecified rate, rm. The main purpose of using this operator is to not get caught in the local optimal solutions. In this paper, use has been made of the random and max–min mutations [1,4,44], each with a predetermined probability.
To create the next generation, the offspring generated by implementation of the crossover and mutation operators are combined with the current population and the nPop best ones are selected. These best solutions, then, form the next-generation population.

2.2.7. Stopping Criteria

The GA is terminated after a prespecified number of iterations, MaxIt.

3. Results

As a numerical example, the new optimization model proposed in this paper is applied to three simple systems composed of series subsystems, for each of which some component choices are available for selection. The first system consists of 6 subsystems and the second one consists of 15 subsystems. The objective is to maximize system reliability at a mission time t = 100 (in arbitrary units). The third example is a case study introduced in [6]. The system is a pharmaceutical plant containing 10 series subsystems as follows:
1. Weighting machine, 2. Sifter machine, 3. Mass machine,
4. Granulator, 5. Fluid bed dryer, 6. Octagonal blender,
7. Rotary compression machine,8. Coating machine, 9. Air compressor,
10. Strip packing machine.
In this work, we extend this case study so that the rate parameters for the time to failure of the components are time-dependent. The mission time is equal to t = 1000.
The rate parameter of the Erlang distribution is assumed to be in the form of a bathtub curve [70,71,72,73,74], as shown in Figure 4. As can be seen, the rate parameter first decreases over time, representing infant mortality, i.e., the process for which defective items fail early and the occurrence rate of the shocks decreases over time as the defective items are weeded out of the population. Then, the rate parameter remains constant for a period of time, as fully random shocks are the cause of failure, independent of time. Finally, the rate parameter increases, due to the aging process, for which the components are more likely to fail as time goes on.
We consider the mathematical form of the rate parameter as in Equation (26) to describe all three parts, decreasing, constant and increasing:
λ ( t ) = c λ 0 ( λ 0 t ) α 1
In the above equation, c is a positive constant parameter. Different values of the parameter α allow representing the different parts of the rate parameter: a value of α < 1 indicates a rate parameter decreasing over time; α = 1 indicates a constant rate parameter over time; α > 1 indicates a rate parameter increasing over time. Then, the rate parameter for the components in Figure 4 can be expressed by Equation (27), in which α1 is a value less than one and α2 is a value larger than one:
λ ( t ) = c 1 λ 0 ( λ 0 t ) α 1 1 ;   0   t     t 1 λ 0 ;   t 1   t     t 2 c 2 λ 0 ( λ 0 t ) α 2 1 ;   t t 2
In this equation, the values of c1 and c2 should be determined so that the rate function is continuous at points t1 and t2. For this purpose, their values should be considered as given in Equation (28):
λ ( t ) = ( λ 0 t 1 ) 1 α 1 λ 0 ( λ 0 t ) α 1 1 = λ 0 ( t t 1 ) α 1 1 ;   0   t     t 1 λ 0 ;   t 1   t     t 2 ( λ 0 t 2 ) 1 α 2 λ 0 ( λ 0 t ) α 2 1 = λ 0 ( t t 2 ) α 2 1 ;   t t 2
According to Equation (28), λ 0 t can be calculated as Equation (29) to be used in Equations (22)–(25).
λ 0 t = 0 t λ ( t ) d t = 0 t λ 0 ( u t 1 ) α 1 1 d u = λ 0 t 1 α 1 ( t t 1 ) α 1 ;   0   t     t 1 0 t 1 λ 0 ( u t 1 ) α 1 1 d u + t 1 t λ 0 d u = λ 0 t 1 α 1 + λ 0 ( t t 1 ) ;   t 1   t     t 2 0 t 1 λ 0 ( u t 1 ) α 1 1 d u + t 1 t 2 λ 0 d u + t 2 t λ 0 ( u t 2 ) α 2 1 d u = λ 0 t 1 α 1 + λ 0 ( t 2 t 1 ) + λ 0 t 2 α 2 ( t t 2 ) α 2 1 ;   t t 2
The parameters α1, α2, t1 and t2 must be determined for each component type allocated in each subsystem. The components characteristics, including cost, weight and rate parameters, for the numerical examples are presented in Table A1, Table A2 and Table A3 in Appendix A. The values of C, W, reliability of the switching system at mission time and the maximum number of components that can be allocated to each subsystem are provided in Table A4 (in Appendix A) for the examples considered.
In order to tune the parameters of the GA, the Taguchi method has been used. For this purpose, three levels have been considered for population size, crossover rate and mutation rate.
The results obtained by the implementation of the Taguchi method are shown in Figure 5. The considered levels, as well as the best values obtained by this method for population size, crossover rate and mutation rate are shown in Table 1. The number of iterations is set to be MaxIt = 100 for all the examples.
As the GA is a stochastic algorithm, 20 trials of the algorithm are performed for each example and the best solution found amongst them is selected as the optimal solution. The trend of the algorithm towards the optimal solution in the 20 trials is shown in Figure A1, Figure A2 and Figure A3 in Appendix A. In order to show the robustness of the algorithm, the values of the minimum, maximum, average, standard deviation and coefficient of variation values of the reliability, weight and cost of the system in the optimal solutions obtained by different trials are reported in Table 2. The values of the reliability, weight and cost of the system in the optimal solutions obtained by different trials are reported in Table A5, Table A6 and Table A7 in Appendix A.
To the best of our knowledge, there is no benchmark problem to compare with the RAP presented in this paper. So, to show the effect of time dependence of the rate parameter, the above-mentioned examples have also been considered with a constant rate parameter (i.e., useful lifetime), of the same value for all components. The best solutions obtained for the examples by considering a constant rate parameter (λ) and by considering a time-dependent rate parameter (λ(t)) (i.e., also initial early failure and final wear-out parts of the lifetime) are given in Table 3, Table 4 and Table 5.

4. Discussion

In this paper, a formulation for RAP is proposed in which the redundancy strategy for each subsystem is taken as the decision variable and the distribution of the time to failure of the components is assumed to have the form of an Erlang distribution, whose rate parameter (λ(t)) is not constant in time. Then, the genetic algorithm is proposed to search for the optimal solution.
As a numerical example, the proposed model was applied to three simple systems composed of series subsystems. As the GA is a stochastic algorithm, 20 trials of the algorithm were performed for each example. The trend of the algorithm towards the optimal solution was shown in Figure A1, Figure A2 and Figure A3. Due to the random nature of the GA, the trend of the algorithm towards the optimal solution varies in different trials. But, it can be seen that the differences are small, which shows the robustness of the algorithm in finding the optimal solution. It can also be seen that, in all three numerical examples, the different trials have converged to a similar solution, with only a slight difference, which shows the convergence of the algorithm to a same optimal solution.
Also, the results obtained by different trials including the values of the minimum, maximum, average, standard deviation and coefficient of variation of the reliability, the weight and cost of the system in the optimal solutions were reported in Table 2. As can be seen, the values of the coefficient of variation of the reliability, cost and weight of the system are small. This shows that the optimal solutions obtained by different trials are only slightly different from each other, which confirms the robustness of the algorithm in finding the optimal solution and its convergence to the optimal solution.
Since we did not find any benchmark problem to compare with the RAP presented in this paper, to show the effect of time dependence of the rate parameter, the examined examples have also been considered with a constant rate parameter, of the same value for all components. By comparing the best solution obtained for the problems with a constant rate parameter and with a time-dependent rate parameter (given in Table 3, Table 4 and Table 5), it can be seen that the structure of the obtained solutions, the type of components allocated to the subsystems and the redundancy strategy of the subsystems may be different.
For example, in the first numerical example, the type of the component allocated to the first subsystem is different in the cases of time-dependent rate parameter and constant rate parameter. In the case with time dependence, component type 3 is allocated to this subsystem, and in the case without time dependence, component type 2 is allocated to this subsystem. In this subsystem, the number of standby components is also different. In the case with time dependence, three standby components are allocated, and in the case without time dependence, one standby component is allocated to this subsystem. As another example, in the first numerical example, the redundancy strategy of the third subsystem is different in the two considered cases. In the case with time dependence, the redundancy strategy is active and, therefore, no switching system is required, but in the case without time dependence, the redundancy strategy is standby and, therefore, a switching system is required. These differences can also be seen in other cases.
In general, it can be seen that the structure of the optimal solution is different in the two cases. In other words, the optimal solution obtained in the case with time dependence is not optimal in the case without time dependence and vice versa. This issue is important in considering assumptions in designing a system. Because, if the rate parameter is variable over time, by considering it as constant, the real optimal solution cannot be found and, therefore, the system designed based on the considered assumptions is no longer optimal and its reliability in practice will be lower than what is theoretically expected. Therefore, to reach a real optimal solution, correct assumptions should be made.
In addition to the difference in the structure of the subsystems, there are also differences in the values of reliability, weight and cost in the obtained optimal solutions at the system level. For the cases studied, the obtained results show that, when the rate parameters are considered as time-dependent, the system reliability is lower. This is expected, because when the time-dependent rate parameters are considered for the components, the failure rates in the early and wear-out parts of the lifetime are higher than the failure rates in the useful part of the lifetime. Therefore, the probability of failure of the components is larger and, as a result, the reliability of the system is smaller than with constant failure rates. This decrease in reliability of the system in the case with time dependence is even accompanied by an increase in the cost of the system in some cases. For example, in the first numerical example, in the case with time dependence, the system reliability is 0.9702 and the system cost is 45, but in the case without time dependence, the system reliability is 0.9877 and the system cost is 37. This means that if the rate parameter is time-dependent, the system cost is higher than the case of a constant rate parameter but the system reliability is lower.
All of these confirm that it is important to consider the actual time-dependence of the rate parameter according to reality.

5. Conclusions

A redundancy allocation problem with multiple redundancy types has been investigated, in which the distribution of time to failure of the components has the form of an Erlang distribution with time-dependent parameter λ(t). The mathematical model has been formulated and a genetic algorithm has been developed to find the optimal solution.
To demonstrate the effect of time dependence of the time to failure distribution, numerical examples have been worked out with and without time dependence of the rate parameter. In the cases studied, it has been shown that taking the time dependence of the parameter λ(t) into account changes the reliability of the components and as a result, both the optimal redundancy structure of the subsystems of the system and the values of its reliability, cost and weight are different from those of the case with constant rate parameter.
In this paper, only identical components have been assumed to be allocated to each subsystem. However, mixing non-identical components might lead to improved system reliability; so, it is suggested that future studies examine the effect of mixing non-identical components. Another aspect to consider is the uncertainty in the characteristics and the time to failure distribution parameters of the components. Considering multi-state, repairable components, other pdfs for the times to failure of the components and/or preventive maintenance is another interesting area for future investigation.
Since the model presented in each research needs to be compared with the previous ones, it can be useful to define benchmark problems that can show the impact of the assumptions of each model well and be used in future research. Therefore, the definition of suitable benchmark problems to be used in the research in this field can also be a suggestion for future research.

Author Contributions

Conceptualization, E.Z. and H.G.; Investigation, E.Z. and H.G.; Methodology, H.G.; Software, H.G.; Supervision, E.Z.; Validation, E.Z. and H.G.; Writing—original draft, H.G.; Writing—review and editing, E.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

In this of article, no new data were created.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Component data for the first numerical example.
Table A1. Component data for the first numerical example.
Component Type 1 (j = 1)Component Type 2 (j = 2)Component Type 3 (j = 3)Component Type 4 (j = 4)
iλ0ijα1ijα2ijt1ijt2ijkijcijwijλ0ijα1ijα2ijt1ijt2ijkijcijwijλ0ijα1ijα2ijt1ijt2ijkijcijwijλ0ijα1ijα2ijt1ijt2ijkijcijwij
10.0520.3310906240.0070.1212902240.0490.42101004320.0810.331080434
20.0170.52151204250.1100.44181004350.1240.335756260.0460.15575443
30.0430.4481105340.0560.63201106450.0280.7410953350.0040.5415100224
40.0260.4312804460.0010.358801380.0040.6312851340.0090.4220120236
50.0230.656703240.0770.8215854260.1330.337905250.1220.731785745
60.0110.779953450.0830.2314955540.0350.22151103260.0430.331490458
Table A2. Component data for the second numerical example.
Table A2. Component data for the second numerical example.
SubsystemComponent Typeλ0ijα1ijα2ijt1ijt2ijkijcijwij
i = 110.0140.552070342
20.0340.46780543
30.0140.2410110236
40.0520.22785754
50.0480.331685826
60.0130.43675234
70.0210.145100342
80.0450.331980755
90.0230.351880458
100.0710.75590724
i = 210.0270.538100353
20.0170.4416105233
30.1050.261180824
40.0190.831980326
50.0290.731080334
60.0490.541980645
70.0530.8311105526
80.0460.151490758
90.0170.83680255
100.0230.761080343
i = 310.0710.651085828
20.0180.33795228
30.0220.631585252
40.0340.621885424
50.0210.521270258
60.0190.261670454
70.0710.641595837
80.0510.852085642
90.0630.6511105756
100.0180.421070242
i = 410.0210.866105227
20.0140.141490336
30.0250.121395523
40.0320.2420105525
50.0470.331175432
60.0170.141080327
70.0490.14575456
80.0510.351785838
90.0210.762075327
100.0430.252095746
i = 510.0730.531085853
20.0480.551085632
30.0220.246105323
40.0550.831185638
50.0540.261195857
60.0170.3310100228
70.0320.65695432
80.0630.361290944
90.0510.761090648
100.0730.547105922
i = 610.0160.821995226
20.0270.13990538
30.0120.252090223
40.0090.26870247
50.0510.752070757
60.0290.5419110443
70.0270.848105337
80.0170.641575247
90.0210.741585234
100.0310.16585435
i = 710.0320.56970647
20.0220.229100336
30.0530.461390745
40.0170.721495255
50.0290.3613105433
60.0340.141390748
70.0370.641670546
80.0510.362075838
90.0240.741695224
100.0150.422090235
i = 810.0180.451085233
20.0230.742075336
30.0120.278110224
40.0240.242080532
50.020.241190353
60.0580.661995725
70.0610.64670824
80.0240.342085326
90.0680.7312105732
100.0730.5514105827
i = 910.0490.72885557
20.0680.661185827
30.0240.661570435
40.0190.231090237
50.0240.761080355
60.0490.731375635
70.0350.2320110655
80.0250.83780336
90.0490.84575628
100.0390.221490545
i = 1010.0280.152090742
20.0480.42590542
30.0190.2220105334
40.0150.241270234
50.0180.34685255
60.0190.861495223
70.0580.361580838
80.0160.146110244
90.0510.561170926
100.0650.3510100823
i = 1110.0330.651480433
20.0270.32585337
30.0350.441590528
40.0090.1210110258
50.0380.261785623
60.0350.3320110533
70.0290.161175623
80.0220.631070323
90.0580.641890728
100.0620.231190822
i = 1210.0220.6414105256
20.0410.251885726
30.0370.56695434
40.0230.566105234
50.0320.661890424
60.0590.73680753
70.0380.332080654
80.0540.442090827
90.0570.745100523
100.0150.52885238
i = 1310.0470.831195527
20.0190.531275255
30.0530.331080748
40.0390.621170453
50.0190.851785232
60.0480.421190626
70.0230.13790348
80.0340.14680537
90.0720.331185928
100.0510.822085543
i = 1410.0190.531675255
20.0380.62785457
30.0450.351985756
40.0330.26995427
50.0450.4312105642
60.0240.23980345
70.0230.838110324
80.0480.331490733
90.0240.831485358
100.0410.22670647
i = 1510.0210.841475242
20.0410.471190553
30.0430.53670548
40.0250.6416110333
50.0120.221895256
60.0390.84775454
70.0470.35990526
80.0130.261495226
90.0150.641890244
100.0440.261675823
Table A3. Component data for the third example (pharmaceutical plant).
Table A3. Component data for the third example (pharmaceutical plant).
Component Type 1Component Type 2Component Type 3Component Type 4
iλ0ij × 105α1ijα2ijt1ijt2ijkijcijwijλ0ij × 105α1ijα2ijt1ijt2ijkijcijwijλ0ij × 105α1ijα2ijt1ijt2ijkijcijwijλ0ij × 105α1ijα2ijt1ijt2ijkijcijwij
11280.5 5 200880 3 16 261760.8 2 50 840 3 20 192760.6 3 140 860 5 18 176070.9 4 190 1060 8 21 24
21950.7 6 120790 4 14 20780.3 4 90 870 2 21 23920.7 2 100 1000 2 13 264300.1 4 80 1010 9 19 18
34190.6 3 2201050 7 20 244730.8 5 150 1090 7 15 173680.2 6 80 720 8 16 264550.8 2 240 1010 7 17 27
44800.5 4 90850 7 19 223880.7 6 120 990 6 17 243940.9 7 80 750 8 20 171590.1 8 90 780 5 19 21
54950.7 6 80820 8 14 191730.5 7 70 990 3 18 22970.5 7 150 900 2 19 19880.9 5 190 840 2 18 24
61010.2 6 150850 3 15 231060.1 2 210 940 4 13 261940.2 4 70 980 4 20 19680.2 6 210 930 2 14 22
7680.3 7 180830 2 14 262450.5 7 160 800 5 15 182410.2 8 240 760 8 21 234550.9 3 100 820 7 14 22
81940.2 4 2201070 5 18 231890.2 5 220 920 5 14 262050.3 4 70 870 4 16 201900.3 4 90 1100 4 20 25
92650.8 6 50780 5 21 274150.8 5 240 830 7 15 192930.7 3 210 1000 5 16 264290.3 5 250 820 9 20 18
10990.1 2 1701050 4 14 185970.6 4 90 840 9 15 252130.7 6 150 910 4 21 234520.2 4 120 990 8 18 25
Table A4. Values of C, W, nmax,i and ρi(t) for the examples.
Table A4. Values of C, W, nmax,i and ρi(t) for the examples.
CWnmax,iρi (t)
First example507060.99
Second example31040060.99
Third example48051960.99
Figure A1. Trend of GA towards the optimal solution in 20 different trials; first example.
Figure A1. Trend of GA towards the optimal solution in 20 different trials; first example.
Mathematics 11 03534 g0a1
Figure A2. Trend of GA towards the optimal solution in 20 different trials; second example.
Figure A2. Trend of GA towards the optimal solution in 20 different trials; second example.
Mathematics 11 03534 g0a2
Figure A3. Trend of GA towards the optimal solution in 20 different trials; third example.
Figure A3. Trend of GA towards the optimal solution in 20 different trials; third example.
Mathematics 11 03534 g0a3
Table A5. Details of the optimal solutions obtained by different trials; first example.
Table A5. Details of the optimal solutions obtained by different trials; first example.
System Characteristics
Iteration NumberReliabilityCostWeight
10.96124568
20.95413770
30.92844570
40.95783968
50.96924869
60.97024570
70.96194966
80.96624668
90.96954570
100.96525070
110.96734968
120.96694968
130.93763968
140.96894869
150.9344766
160.96243868
170.95494765
180.95774570
190.93154669
200.95624668
Table A6. Details of the optimal solutions obtained by different trials; second example.
Table A6. Details of the optimal solutions obtained by different trials; second example.
System Characteristics
Iteration NumberReliabilityCostWeight
10.9626200265
20.9506306367
30.9614259302
40.9557251380
50.9425259352
60.9584250337
70.9397271372
80.9574303319
90.9584249339
100.9351242325
110.9337235311
120.9384237298
130.9565219393
140.9643248333
150.9488236285
160.9707272282
170.946240323
180.9524277384
190.9493302384
200.9628260343
Table A7. Details of the optimal solutions obtained by different trials; third example.
Table A7. Details of the optimal solutions obtained by different trials; third example.
System Characteristics
Iteration NumberReliabilityCostWeight
10.9866458519
20.9865448518
30.9881466517
40.9846478518
50.9896454514
60.9894453516
70.9889469513
80.9881469519
90.9886464509
100.9894453516
110.9889469513
120.9869471512
130.9864429518
140.9896454514
150.9868451519
160.9891468510
170.9869424509
180.9896467517
190.9896461519
200.9894453516

References

  1. Ardakan, M.A.; Hamadani, A.Z. Reliability optimization of series–parallel systems with mixed redundancy strategy in subsystems. Reliab. Eng. Syst. Saf. 2014, 130, 132–139. [Google Scholar] [CrossRef]
  2. Ebeling, C.E. An Introduction to Reliability and Maintainability Engineering; Tata McGraw-Hill Education: New York, NY, USA, 2004. [Google Scholar]
  3. Hsieh, T.J. Hierarchical redundancy allocation for multi-level reliability systems employing a bacterial-inspired evolutionary algorithm. Inf. Sci. 2014, 288, 174–193. [Google Scholar] [CrossRef]
  4. Abouei Ardakan, M.; Mirzaei, Z.; Zeinal Hamadani, A.; Elsayed, E.A. Reliability Optimization by Considering Time-Dependent Reliability for Components. Qual. Reliab. Eng. Int. 2017, 33, 1641–1654. [Google Scholar] [CrossRef]
  5. Gupta, R.K.; Bhunia, A.K.; Roy, D. A GA based penalty function technique for solving constrained redundancy allocation problem of series system with interval valued reliability of components. J. Comput. Appl. Math. 2009, 232, 275–284. [Google Scholar] [CrossRef]
  6. Mellal, M.A.; Zio, E. System reliability-redundancy allocation by evolutionary computation. In Proceedings of the 2017 2nd International Conference on System Reliability and Safety (ICSRS) 2017, Milan, Italy, 20–22 December 2017; pp. 15–19. [Google Scholar]
  7. Coit, D.W. Cold-standby redundancy optimization for nonrepairable systems. Iie Trans. 2001, 33, 471–478. [Google Scholar] [CrossRef]
  8. Zhao, J.; Zeng, S.; Guo, J.; Yang, C. Redundancy allocation with non-identical component and uncertainty. In Proceedings of the 2015 Annual Reliability and Maintainability Symposium (RAMS) 2015, Palm Harbor, FL, USA, 26–29 January 2015; pp. 1–6. [Google Scholar]
  9. Liang, Y.C.; Smith, A.E. An ant colony optimization algorithm for the redundancy allocation problem (RAP). IEEE Trans. Reliab. 2004, 53, 417–423. [Google Scholar] [CrossRef]
  10. Onishi, J.; Kimura, S.; James, R.J.; Nakagawa, Y. Solving the redundancy allocation problem with a mix of components using the improved surrogate constraint method. IEEE Trans. Reliab. 2007, 56, 94–101. [Google Scholar] [CrossRef]
  11. Caserta, M.; Voß, S. An exact algorithm for the reliability redundancy allocation problem. Eur. J. Oper. Res. 2015, 244, 110–116. [Google Scholar] [CrossRef]
  12. Ramirez-Marquez, J.E.; Coit, D.W. A heuristic for solving the redundancy allocation problem for multi-state series-parallel systems. Reliab. Eng. Syst. Saf. 2004, 83, 341–349. [Google Scholar] [CrossRef]
  13. Sharma, V.K.; Agarwal, M.; Sen, K. Reliability evaluation and optimal design in heterogeneous multi-state series-parallel systems. Inf. Sci. 2011, 181, 362–378. [Google Scholar] [CrossRef]
  14. Lai, C.M.; Yeh, W.C. Two-stage simplified swarm optimization for the redundancy allocation problem in a multi-state bridge system. Reliab. Eng. Syst. Saf. 2016, 156, 148–158. [Google Scholar] [CrossRef]
  15. Li, Y.; Zio, E. A quantum-inspired evolutionary approach for non-homogeneous redundancy allocation in series-parallel multi-state systems. In Proceedings of the 2014 10th International Conference on Reliability, Maintainability and Safety (ICRMS) 2014, Guangzhou, China, 6–8 August 2014; pp. 526–532. [Google Scholar]
  16. Wang, Y.; Li, L. A PSO algorithm for constrained redundancy allocation in multi-state systems with bridge topology. Comput. Ind. Eng. 2014, 68, 13–22. [Google Scholar] [CrossRef]
  17. Hsieh, T.J.; Yeh, W.C. Penalty guided bees search for redundancy allocation problems with a mix of components in series–parallel systems. Comput. Oper. Res. 2012, 39, 2688–2704. [Google Scholar] [CrossRef]
  18. Garg, H.; Rani, M.; Sharma, S.P.; Vishwakarma, Y. Bi-objective optimization of the reliability-redundancy allocation problem for series-parallel system. J. Manuf. Syst. 2014, 33, 335–347. [Google Scholar] [CrossRef]
  19. Zhang, E.; Chen, Q. Multi-objective reliability redundancy allocation in an interval environment using particle swarm optimization. Reliab. Eng. Syst. Saf. 2016, 145, 83–92. [Google Scholar] [CrossRef]
  20. Guo, J.; Wang, Z.; Zheng, M.; Wang, Y. Uncertain multiobjective redundancy allocation problem of repairable systems based on artificial bee colony algorithm. Chin. J. Aeronaut. 2014, 27, 1477–1487. [Google Scholar] [CrossRef]
  21. Kulturel-Konak, S.; Coit, D.W. Determination of Pruned Pareto sets for the multi-objective system redundancy allocation problem. In Proceedings of the 2007 IEEE Symposium on Computational Intelligence in Multi-Criteria Decision-Making 2007, Honolulu, HI, USA, 1–5 April 2007; pp. 390–394. [Google Scholar]
  22. Mariano, C.H.; Kuri-Morales, A.F. Complex componential approach for redundancy allocation problem solved by simulation-optimization framework. J. Intell. Manuf. 2014, 25, 661–680. [Google Scholar] [CrossRef]
  23. Cao, D.; Murat, A.; Chinnam, R.B. Efficient exact optimization of multi-objective redundancy allocation problems in series-parallel systems. Reliab. Eng. Syst. Saf. 2013, 111, 154–163. [Google Scholar] [CrossRef]
  24. Sun, M.X.; Li, Y.F.; Zio, E. On the optimal redundancy allocation for multi-state series–parallel systems under epistemic uncertainty. Reliab. Eng. Syst. Saf. 2019, 192, 106019. [Google Scholar] [CrossRef]
  25. Valaei, M.R.; Behnamian, J. Allocation and sequencing in 1-out-of-N heterogeneous cold-standby systems: Multi-objective harmony search with dynamic parameters tuning. Reliab. Eng. Syst. Saf. 2017, 157, 78–86. [Google Scholar] [CrossRef]
  26. Bei, X.; Chatwattanasiri, N.; Coit, D.W.; Zhu, X. Combined redundancy allocation and maintenance planning using a two-stage stochastic programming model for multiple component systems. IEEE Trans. Reliab. 2017, 66, 950–962. [Google Scholar] [CrossRef]
  27. Dobani, E.R.; Ardakan, M.A.; Davari-Ardakani, H.; Juybari, M.N. RRAP-CM: A new reliability-redundancy allocation problem with heterogeneous components. Reliab. Eng. Syst. Saf. 2019, 191, 106563. [Google Scholar] [CrossRef]
  28. Hsieh, T.J. Component mixing with a cold standby strategy for the redundancy allocation problem. Reliab. Eng. Syst. Saf. 2021, 206, 107290. [Google Scholar] [CrossRef]
  29. Bhattacharyee, N.; Kumar, N.; Mahato, S.K.; Bhunia, A.K. Development of a blended particle swarm optimization to optimize mission design life of a series–parallel reliable system with time dependent component reliabilities in imprecise environments. Soft Comput. 2021, 25, 11745–11761. [Google Scholar] [CrossRef]
  30. Wang, W.; Wu, Z.; Xiong, J.; Xu, Y. Redundancy optimization of cold-standby systems under periodic inspection and maintenance. Reliab. Eng. Syst. Saf. 2018, 180, 394–402. [Google Scholar] [CrossRef]
  31. Mellal, M.A.; Zio, E. System reliability-redundancy optimization with cold-standby strategy by an enhanced nest cuckoo optimization algorithm. Reliab. Eng. Syst. Saf. 2020, 201, 106973. [Google Scholar] [CrossRef]
  32. Yeh, W.C.; Su, Y.Z.; Gao, X.Z.; Hu, C.F.; Wang, J.; Huang, C.L. Simplified swarm optimization for bi-objection active reliability redundancy allocation problems. Appl. Soft Comput. 2021, 106, 107321. [Google Scholar] [CrossRef]
  33. Yeh, W.C. BAT-based algorithm for finding all Pareto solutions of the series-parallel redundancy allocation problem with mixed components. Reliab. Eng. Syst. Saf. 2022, 228, 108795. [Google Scholar] [CrossRef]
  34. Kundu, P. A multi-objective reliability-redundancy allocation problem with active redundancy and interval type-2 fuzzy parameters. Oper. Res. 2021, 21, 2433–2458. [Google Scholar] [CrossRef]
  35. Coit, D.W. Maximization of system reliability with a choice of redundancy strategies. IIE Trans. 2003, 35, 535–543. [Google Scholar] [CrossRef]
  36. Kong, X.; Gao, L.; Ouyang, H.; Li, S. Solving the redundancy allocation problem with multiple strategy choices using a new simplified particle swarm optimization. Reliab. Eng. Syst. Saf. 2015, 144, 147–158. [Google Scholar] [CrossRef]
  37. Tavakkoli-Moghaddam, R.; Safari, J. A new mathematical model for a redundancy allocation problem with mixing components redundant and choice of redundancy strategies. Appl. Math. Sci. 2007, 45, 2221–2230. [Google Scholar]
  38. Qiu, X.; Ali, S.; Yue, T.; Zhang, L. Reliability-redundancy-location allocation with maximum reliability and minimum cost using search techniques. Inf. Softw. Technol. 2017, 82, 36–54. [Google Scholar] [CrossRef]
  39. Kim, H. Optimal reliability design of a system with k-out-of-n subsystems considering redundancy strategies. Reliab. Eng. Syst. Saf. 2017, 167, 572–582. [Google Scholar] [CrossRef]
  40. Kim, H.; Kim, P. Reliability–redundancy allocation problem considering optimal redundancy strategy using parallel genetic algorithm. Reliab. Eng. Syst. Saf. 2017, 159, 153–160. [Google Scholar] [CrossRef]
  41. Kim, H. Parallel genetic algorithm with a knowledge base for a redundancy allocation problem considering the sequence of heterogeneous components. Expert Syst. Appl. 2018, 113, 328–338. [Google Scholar] [CrossRef]
  42. Wang, W.; Lin, M.; Fu, Y.; Luo, X.; Chen, H. Multi-objective optimization of reliability-redundancy allocation problem for multi-type production systems considering redundancy strategies. Reliab. Eng. Syst. Saf. 2020, 193, 106681. [Google Scholar] [CrossRef]
  43. Li, X.Y.; Li, Y.F.; Huang, H.Z. Redundancy allocation problem of phased-mission system with non-exponential components and mixed redundancy strategy. Reliab. Eng. Syst. Saf. 2020, 199, 106903. [Google Scholar] [CrossRef]
  44. Gholinezhad, H.; Hamadani, A.Z. A new model for the redundancy allocation problem with component mixing and mixed redundancy strategy. Reliab. Eng. Syst. Saf. 2017, 164, 66–73. [Google Scholar] [CrossRef]
  45. Ardakan, M.A.; Hamadani, A.Z.; Alinaghian, M. Optimizing bi-objective redundancy allocation problem with a mixed redundancy strategy. ISA transactions 2015, 55, 116–128. [Google Scholar] [CrossRef] [PubMed]
  46. Peiravi, A.; Karbasian, M.; Ardakan, M.A.; Coit, D.W. Reliability optimization of series-parallel systems with K-mixed redundancy strategy. Reliab. Eng. Syst. Saf. 2019, 183, 17–28. [Google Scholar] [CrossRef]
  47. Peiravi, A.; Ardakan, M.A.; Zio, E. A new Markov-based model for reliability optimization problems with mixed redundancy strategy. Reliab. Eng. Syst. Saf. 2020, 201, 106987. [Google Scholar] [CrossRef]
  48. Ouyang, Z.; Liu, Y.; Ruan, S.J.; Jiang, T. An improved particle swarm optimization algorithm for reliability-redundancy allocation problem with mixed redundancy strategy and heterogeneous components. Reliab. Eng. Syst. Saf. 2019, 181, 62–74. [Google Scholar] [CrossRef]
  49. Reihaneh, M.; Ardakan, M.A.; Eskandarpour, M. An exact algorithm for the redundancy allocation problem with heterogeneous components under the mixed redundancy strategy. Eur. J. Oper. Res. 2022, 297, 1112–1125. [Google Scholar] [CrossRef]
  50. Azizi, S.; Mohammadi, M. Strategy selection for multi-objective redundancy allocation problem in a k-out-of-n system considering the mean time to failure. Opsearch 2023, 60, 1021–1044. [Google Scholar] [CrossRef]
  51. Kim, H. Parallel Genetic Algorithm with Knowledge Archives for the Redundancy Allocation Problem in a Mixed Redundant System. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=4283392 (accessed on 22 November 2022).
  52. Zhang, J.; Lv, H.; Hou, J. A novel general model for RAP and RRAP optimization of k-out-of-n: G systems with mixed redundancy strategy. Reliab. Eng. Syst. Saf. 2023, 229, 108843. [Google Scholar] [CrossRef]
  53. Yeh, W.C. Solving cold-standby reliability redundancy allocation problems using a new swarm intelligence algorithm. Appl. Soft Comput. 2019, 83, 105582. [Google Scholar] [CrossRef]
  54. Yeh, W.C. A novel boundary swarm optimization method for reliability redundancy allocation problems. Reliab. Eng. Syst. Saf. 2019, 192, 106060. [Google Scholar] [CrossRef]
  55. Huang, X.; Coolen, F.P.; Coolen-Maturi, T. A heuristic survival signature based approach for reliability-redundancy allocation. Reliab. Eng. Syst. Saf. 2019, 185, 511–517. [Google Scholar] [CrossRef]
  56. Lins, I.D.; Droguett, E.L. Redundancy allocation problems considering systems with imperfect repairs using multi-objective genetic algorithms and discrete event simulation. Simul. Model. Pract. Theory 2011, 19, 362–381. [Google Scholar] [CrossRef]
  57. Marseguerra, M.; Zio, E.; Podofillini, L.; Coit, D.W. Optimal design of reliable network systems in presence of uncertainty. IEEE Trans. Reliab. 2005, 54, 243–253. [Google Scholar] [CrossRef]
  58. Monalisa, P.; Satchidananda, D.; Kumar, J.A. Multi-objective artificial bee colony algorithm in redundancy allocation problem. Int. J. Adv. Intell. Paradig. 2023, 25, 24–50. [Google Scholar] [CrossRef]
  59. Maji, A.; Duary, A.; Bhunia, A.K.; Mondal, S.K. A Redundancy Allocation Problem for a Series-Parallel System with Multiple Choice Technologies Considering Fuzzy Sense with Ambiguity and Vagueness; Springer: Berlin/Heidelberg, Germany, 2022. [Google Scholar] [CrossRef]
  60. Yeh, C.T.; Fiondella, L. Optimal redundancy allocation to maximize multi-state computer network reliability subject to correlated failures. Reliab. Eng. Syst. Saf. 2017, 166, 138–150. [Google Scholar] [CrossRef]
  61. Li, Y.F.; Zhang, H. The methods for exactly solving redundancy allocation optimization for multi-state series–parallel systems. Reliab. Eng. Syst. Saf. 2022, 221, 108340. [Google Scholar] [CrossRef]
  62. Xu, Y.; Pi, D.; Yang, S.; Chen, Y. A novel discrete bat algorithm for heterogeneous redundancy allocation of multi-state systems subject to probabilistic common-cause failure. Reliab. Eng. Syst. Saf. 2021, 208, 107338. [Google Scholar] [CrossRef]
  63. Li, J.; Wang, G.; Zhou, H.; Chen, H. Redundancy allocation optimization for multi-state system with hierarchical performance requirements. Proc. Inst. Mech. Eng. Part O J. Risk Reliab. 2022. [Google Scholar] [CrossRef]
  64. Blokus, A. Multistate System Reliability with Dependencies; Academic Press: Cambridge, MA, USA, 2020. [Google Scholar]
  65. Yingkui, G.; Jing, L. Multi-state system reliability: A new and systematic review. Procedia Eng. 2012, 29, 531–536. [Google Scholar] [CrossRef]
  66. Paramanik, R.; Mahato, S.K.; Bhattacharyee, N. Optimization of Redundancy Allocation Problem Using Quantum Particle Swarm Optimization Algorithm Under Uncertain Environment. In Advances in Reliability, Failure and Risk Analysis; Springer Nature: Singapore, 2023; pp. 177–197. [Google Scholar]
  67. Coit, D.W.; Zio, E. The evolution of system reliability optimization. Reliab. Eng. Syst. Saf. 2019, 192, 106259. [Google Scholar] [CrossRef]
  68. Chern, M.S. On the computational complexity of reliability redundancy allocation in a series system. Oper. Res. Lett. 1992, 11, 309–315. [Google Scholar] [CrossRef]
  69. Holland, J.H. Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence; MIT Press: Cambridge, MA, USA, 1992. [Google Scholar]
  70. Scheidegger, A.; Leitao, J.P.; Scholten, L. Statistical failure models for water distribution pipes–A review from a unified perspective. Water Res. 2015, 83, 237–247. [Google Scholar] [CrossRef]
  71. Romaniuk, M.; Hryniewicz, O. Estimation of maintenance costs of a pipeline for a U-shaped hazard rate function in the imprecise setting. Eksploat. I Niezawodn. 2020, 22, 352–362. [Google Scholar] [CrossRef]
  72. Al-nefaie, A.H.; Ragab, I.E. A Novel Lifetime Model with A Bathtub-Shaped Hazard Rate: Properties & Applications. J. Appl. Sci. Eng. 2022, 26, 413–421. [Google Scholar]
  73. Bai, B.; Li, Z.; Wu, Q.; Zhou, C.; Zhang, J. Fault data screening and failure rate prediction framework-based bathtub curve on industrial robots. Ind. Robot. Int. J. Robot. Res. Appl. 2020, 47, 867–880. [Google Scholar] [CrossRef]
  74. Ahsan, S.; Lemma, T.A.; Gebremariam, M.A. Reliability analysis of gas turbine engine by means of bathtub-shaped failure rate distribution. Process Saf. Prog. 2020, 39, e12115. [Google Scholar] [CrossRef]
Figure 1. Structure of a system of redundant subsystems.
Figure 1. Structure of a system of redundant subsystems.
Mathematics 11 03534 g001
Figure 2. Chromosome representation.
Figure 2. Chromosome representation.
Mathematics 11 03534 g002
Figure 3. Structure of the system corresponding to the solution in Figure 2.
Figure 3. Structure of the system corresponding to the solution in Figure 2.
Mathematics 11 03534 g003
Figure 4. Rate parameter (bathtub curve).
Figure 4. Rate parameter (bathtub curve).
Mathematics 11 03534 g004
Figure 5. The results obtained by the implementation of Taguchi method: (a) first example; (b) second example; (c) third example.
Figure 5. The results obtained by the implementation of Taguchi method: (a) first example; (b) second example; (c) third example.
Mathematics 11 03534 g005aMathematics 11 03534 g005b
Table 1. The results obtained by the implementation of Taguchi method.
Table 1. The results obtained by the implementation of Taguchi method.
LevelnPoprcrm
11000.20.2
22000.30.3
33000.40.4
Best levelExample 12000.40.2
Example 23000.30.3
Example 32000.40.3
Table 2. The results obtained by different trials.
Table 2. The results obtained by different trials.
System Characteristics
First ExampleReliabilityCostWeight
Minimum0.92843765
Maximum0.97025070
First exampleAverage0.957145.1568.4
Standard deviation0.01313.78521.4283
Coefficient of variation0.01370.08380.0209
Minimum0.9337200265
Maximum0.9707306393
Second exampleAverage0.9522255.8334.7
Standard deviation0.010226.564336.5392
Coefficient of variation0.01070.10380.1092
Minimum0.9846424509
Maximum0.9896478519
Third exampleAverage0.9882457.95515.3
Standard deviation0.001413.18133.2879
Coefficient of variation0.00150.02880.0064
Table 3. Optimal solution with time dependence and without time dependence; first example.
Table 3. Optimal solution with time dependence and without time dependence; first example.
SubsystemWith Time DependenceWithout Time Dependence
zinAinSiRed. StrategyReliabilityzinAinSiRed. StrategyReliability
1323Mixed0.9938221Mixed0.9994
2111Standby0.9977111Standby0.9987
3420Active0.9939411Standby0.9986
4211Standby0.9909421Mixed0.9981
5122Mixed0.9959112Standby0.9953
6111Standby0.9977111Standby0.9980
System reliability
System cost
System weight
0.9702
45
70
0.9877
37
70
Table 4. Optimal solution with time dependence and without time dependence; second example.
Table 4. Optimal solution with time dependence and without time dependence; second example.
SubsystemWith Time DependenceWithout Time Dependence
zinAinSiRed. StrategyReliabilityzinAinSiRed. StrategyReliability
1131Mixed0.99997751Mixed0.99788
2221Mixed0.99835942Mixed0.99938
3322Mixed0.99999311Standby0.99977
4933Mixed0.99912623Mixed0.99995
5923Mixed0.99428322Mixed0.99876
6822Mixed0.99916132Mixed0.99888
71032Mixed0.99968433Mixed0.99755
8523Mixed0.99717334Mixed0.99764
9823Mixed0.99796741Mixed0.99554
10633Mixed0.99567533Mixed0.99686
11123Mixed0.99761122Mixed0.99853
12113Standby0.999791032Mixed0.99915
13223Mixed0.99877223Mixed0.99921
14923Mixed0.99365833Mixed0.99493
15523Mixed0.99914523Mixed0.99921
System reliability
System cost
System weight
0.9707
272
282
0.9736
254
335
Table 5. Optimal solution with time dependence and without time dependence; third example.
Table 5. Optimal solution with time dependence and without time dependence; third example.
SubsystemWith Time DependenceWithout Time Dependence
zinAinSiRed. StrategyReliabilityzinAinSiRed. StrategyReliability
1321Mixed0.99956321Mixed0.99980
2121Mixed0.99928421Mixed0.99999
3221Mixed0.99950221Mixed0.99966
4321Mixed0.99957411Standby0.99976
5111Standby0.99740121Mixed0.99987
6321Mixed0.99918321Mixed0.99982
7221Mixed0.99923320Active0.99999
8321Mixed0.99924121Mixed0.99998
9211Mixed0.99755421Mixed0.99999
10121Mixed0.99908111Standby0.99981
System reliability
System cost
System weight
0.9896
454
514
0.9987
480
517
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zio, E.; Gholinezhad, H. Redundancy Allocation of Components with Time-Dependent Failure Rates. Mathematics 2023, 11, 3534. https://doi.org/10.3390/math11163534

AMA Style

Zio E, Gholinezhad H. Redundancy Allocation of Components with Time-Dependent Failure Rates. Mathematics. 2023; 11(16):3534. https://doi.org/10.3390/math11163534

Chicago/Turabian Style

Zio, Enrico, and Hadi Gholinezhad. 2023. "Redundancy Allocation of Components with Time-Dependent Failure Rates" Mathematics 11, no. 16: 3534. https://doi.org/10.3390/math11163534

APA Style

Zio, E., & Gholinezhad, H. (2023). Redundancy Allocation of Components with Time-Dependent Failure Rates. Mathematics, 11(16), 3534. https://doi.org/10.3390/math11163534

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop