Next Article in Journal
A Complex Insight for Quality of Service Based on Spreading Dynamics and Multilayer Networks in a 6G Scenario
Previous Article in Journal
A Review of Mathematical Models of Elasticity Theory Based on the Methods of Iterative Factorizations and Fictitious Components
Previous Article in Special Issue
Quantification of Model Uncertainty Based on Variance and Entropy of Bernoulli Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Preventive Maintenance of k-out-of-n System with Dependent Failures

by
Vladimir Rykov
1,2,† and
Olga Kochueva
1,*,†
1
Department of Applied Mathematics and Computer Modeling, National University of Oil and Gas “Gubkin University”, 65, Leninsky Prospekt, 119991 Moscow, Russia
2
Kharkevich Institute for Information Transmission Problems of Russian Academy of Sciences (IITP RAS), 19, Building 1, B. Karetny Per., 127051 Moscow, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2023, 11(2), 422; https://doi.org/10.3390/math11020422
Submission received: 20 December 2022 / Revised: 10 January 2023 / Accepted: 11 January 2023 / Published: 13 January 2023
(This article belongs to the Special Issue Sensitivity Analysis)

Abstract

:
The paper investigates a model of a k-out-of-n system, the residual lifetime of which changes after failures of any of its components. The problem of a Preventive Maintenance (PM) organization as advice to the Decision Maker (DM) for such a system is considered. The purpose of this paper is to propose a mathematical model of the k-out-of-n system to support DM about PM. For most practical applications, it is usually possible to estimate the lifetime distribution parameters of the system components with limited accuracy (only one or two moments), which is why special attention is paid to the sensitivity analysis of the system reliability characteristics and decisions about PM to the shape of system components lifetime distributions. In the numerical examples, we consider the 3-out-of-6 model discussed in our previous works for two real systems. The novelty, significance, and features of this study consist of the following, after the failure of one of the system components, the load on all the others increases, which leads to a decrease in their residual lifetime. We model this situation with order statistics and study the quality of PM strategies with respect to the availability maximization criterion. At the same time, we are focusing on the study of the sensitivity of decision-making to the type of lifetime distribution of system components.

1. Introduction, Motivation and Examples

1.1. Introduction

Maintaining the reliability of objects, processes, and systems (especially dangerous in operation) at a high level is one of the primary tasks of their design, creation, and operation. For the solution of such a complex technical problem, it is necessary to develop adequate mathematical models and methods. Redundancy is a traditional way to improve the reliability of systems, objects, and processes. k-out-of-n models are a good configuration for describing redundant systems. A system consisting of n parallel components can be described in two ways, dual to each other, depending on the definition of the parameter k, as follows. The parameter k can represent the number of components in the system that must work for the entire system to work, referred to as a k-out-of-n: G system. On the other hand, parameter k may represent the number of components in the system that must fail to cause the failure of the entire system, referred to as a k-out-of-n: F system [1]. For our purposes, it is more convenient to consider a subset of failed system components; in this paper, we will use the second type of system description, omitting the symbol F for simplicity.
The study of the k-out-of-n system’s reliability is the subject of interest both from theoretical and practical aspects. From a theoretical point of view, it provides ample opportunities for the development of new mathematical methods. From a practical point of view, as shown by numerous publications, models of k-out-of-n systems can be applied to many real-world phenomena, including telecommunication, transmission, transportation, manufacturing, and service applications. The wide range of practical applications has led to a large number of papers devoted to the study of k-out-of-n systems. The literature on such studies is extensive (see, for example [2,3] and others [4,5,6,7,8,9]). Some applications of this model and the latest overview can be found in [10].
Engineering tasks pose ever-new formulations of problems for researchers. Thus, ref. [11] presents a problem of complex system investigation, the failures of which depend not only on the number of its failed components but also on their location in the system. Moreover, in real systems, the failure of one of its components often leads to a redistribution of the load on the survivors, which leads to a decrease in their residual lifetime. Some of the problems of system reliability studies, the failures of which depend not only on the number of failed components but also on their position in the system, have been considered in this paper.
In the present paper, we consider the situation when the failure of one of the system components leads to an increase in the load on the others, which leads to a decrease in their residual lifetime. A number of works are devoted to the problem of studying the reliability of the system under conditions of load distribution on the system after a failure of one of its components (see, for example, [12,13,14,15]). In the paper [16], the reliability function of a k-out-of-n system under the condition of a load-sharing mechanism has been studied. In this paper, we prolong this study, aiming to study the Preventive Maintenance (PM) optimization for this system with respect to its availability maximization.
Another way to maintain a high level of reliability of systems during operation is to organize their preventive maintenance. A probabilistic study of a real-world k-out-of-n system often helps to develop an optimal strategy for maintaining high system-level reliability. The idea of increasing system reliability by organizing PM is a long-standing one. A comprehensive review of PM methods can be found in the monograph of Gertsbakh [17]. Some recent developments in optimal maintenance policies can be found in [18,19,20,21,22]. Another approach to the problem of the PM optimization for linear consecutive k-out-of-n systems with dependent failures in paper [23] has been proposed.
The purpose of this paper is to propose a mathematical model to study and compare the effectiveness of different PM strategies for k-out-of-n systems based on monitoring their states as advice to DM. In this direction, in [24,25], an approach was developed to compare PM strategies according to the availability and the cost type criteria maximization for models without taking into account the impact of failures of its components on the residual life of survivors. In this paper, we take into account this factor.
Moreover, since the detailed initial information about the reliability of system components are generally lacking, and their lifetimes are usually known with an accuracy of only one or two their moments, it is of utmost importance to study the sensitivity of the system reliability indicators to the shape of their lifetime distributions. Sensitivity analysis is a topical issue attracting the attention of researchers, in particular, ref. [26] studies the role of model uncertainties in sensitivity and probability analysis of reliability and ref. [27] presents new sensitivity measures in a reliability-oriented global sensitivity analysis. Some related research can be found in a series of our papers, references to which are shown in chapter 9 of [28], as well as in [29]. Thus, special attention in this paper is paid to the sensitivity analysis of the system reliability characteristics and decision-making based on the choice of PM to the shape of its lifetime distributions.
This paper is closely related to [16,24,25], where, as well as in this one, we focus on mathematical models and methods for PM organizing for k-out-of-n systems. Since the topics and research tools in all these articles are quite close, some parts of the introduction and notations coincide with those that are used in the above-mentioned papers. However, here other aspects of the model are developed and new results are presented.
The feature of this study is that after the failure of one of the system components, the load on all the others increases, which leads to a decrease in their residual lifetime. We model this situation with order statistics and develop the procedure for comparing PM strategies with respect to the availability maximization criterion. This contributes to the novelty and significance of this paper. At the same time, we place special emphasis on investigating the sensitivity of decision-making to the type of lifetime distribution of system components.
The paper is organized as follows. In the next subsection, two examples of real-world systems for which the choice of PM strategy has to be justified are considered. In what follows, one of these examples is addressed as a demonstration of our theoretical research and proposed algorithms. Section 2 presents the problem statement and some notations. In Section 3, system behavior is modeled by the regenerative process, and the condition of strategy preference is stated. Further, in Section 4, we consider the distributions of the times to PM beginning and the system failure. An algorithm for the best PM strategy choice is also proposed in this section. Numerical experiments based on the example represented in Section 1.2 will be discussed in Section 5. In the Conclusion, we summarize the obtained results and outline further studies.

1.2. Examples

As examples of a k-out-of-n model application, we consider two real-world systems. One of them is described in [24], where an automated system for the remote monitoring of underwater sections of a gas pipeline was considered. One of the essential parts of the system is an unmanned underwater vehicle (UUV) (see Figure 1a), which carries out a visual inspection of underwater areas with the fixation of video materials. It is fitted with various equipment for monitoring, collecting, and transmitting information about the state of the object. Six engines located on different sides of the apparatus allow it to descend, float, and move along the pipe. UUV fails when three engines fail. In addition, the failure of each engine increases the load on the remaining ones, which leads to a decrease in their remaining lifetime. Thus, this system can be represented as a 3-out-of-6 model.
In papers [11,30], a similar model was used for modeling a tethered high-altitude unmanned telecommunication platform for transmitting information on the “last mile” (see Figure 1b). This system can be used in different versions: for n = 6 or for n = 8 , etc. Different ways of exploitation and failure are also possible. In our numerical experiments, we will refer to this example by modeling it with a 3-out-of-6 model.
Both examples will be used in numerical experiments to demonstrate the impact of system component failures on the decision to choose a certain PM strategy and its sensitivity to the initial information. Numerical analysis of the model is presented in Section 5.

2. The Problem Set, Notations, and Assumptions

2.1. Notations and Assumptions

We consider a k-out-of-n system that consists of n components in parallel and fails if at least k of them fail. It is assumed that the failure of any component leads to an increase in the load on others and, consequently, to a decrease in their residual lifetimes. To increase the reliability of the system, it is supposed that PM based on the observation of the state of the system is possible. We denote by L = { 1 , 2 , , k } the set of possible PM strategies, including running to the system failure for l = k . The l-th strategy ( l < k ) means that the PM begins when l components fail. The strategy l = k means that the PM is not used, and the best solution is to wait for the system to fail and start repairing it jointly with maintenance.
In this paper, we will use the following notations and main assumptions:
  • P { · } , E [ · ] , Var [ · ] —symbols of probability, expectation, and variation;
  • A i : ( i = 1 , 2 , ) is the series of component’s lifetimes, which are assumed to be independent identically distributed random variables;
  • A ( t ) = P { A i t } is their common cumulative distribution function (cdf) with their probability density function (pdf) a ( t ) = A ( t ) and expectation
    a = E [ A i ] = 0 ( 1 A i ( t ) ) d t ;
  • The times B i ( l ) ( i = 1 , 2 , ) , ( l L ) of PM and the system repair for l = k are supposed to be independent identically distributed random variables;
  • B l ( t ) = P { B i ( l ) t } are their cdf with mean values;
    b l = E [ B i ( l ) ] = 0 ( 1 B l ( t ) ) d t ;
  • It is assumed that the mean PM times b l ( l < k ) are less than the mean repair time b k , b l b k , and may or may not depend on the type of maintenance;
  • j is the system state, where j means the number of failed components;
  • E = { j = { 0 , 1 , , k } } is the set of the system states;
  • After failure of the l-th component ( l < k ) , the load on all others increases, which leads to a decrease in their residual lifetimes. It is modeled by compression of the residual components lifetime by some weight coefficient w i ( 0 w i 1 ) ;
  • In the initial time epoch, the system is absolutely reliable, i.e., all its components are in the UP states and the initial system state is j = 0 .
  • After any system repair and its PM completion, the system became “as a new one”, i.e., returns to the zero state.

2.2. The Problem Statement

The problem of comparing the PM strategies with respect to the availability of a system with independent failures was considered in paper [24]. This paper continues this study for a model where after the failure of any system component, the load on the other ones increase, resulting in a decrease in the residual lifetime, aiming to maximize the availability factor K av . , l (other quality criteria are also applicable, such as the productivity of the system and/or system service cost under different maintenance strategies (see [25]))
K av . , l = lim t 1 t { the system working time during time t under strategy l } ,
the formal definition of which is represented in (14). To achieve this goal, we define a random process J = { J ( t ) : t 0 } with the set of space E by the relation
J ( t ) = j , if at time t the system is in the state j E .
Using this process, we study the main reliability characteristics of a k-out-of-n system under the conditions when the failure of one of its components leads to an increase in the load on the others and, therefore, to a decrease in their residual lifetimes. These are:
  • Y l , ( l L ) —system working time during its regeneration cycle (see formal definition in Formula (13)).
  • R k ( t ) —system reliability function, and M k —mean time until the system fails (after the failure of k components).
    R k ( t ) = P { Y k > t } , M k = 0 R k ( t ) d t .
  • F l ( t ) —cdf of PM start times for different strategies ( l L ) , and M l —their expectations
    F l ( t ) = P { Y l t } , M l = 0 ( 1 F l ( t ) ) d t .
  • The system availability K av . , l , l L given by (3) for different PM strategies.
  • An indicator of strategy preference criterion, which will be introduced later in the next section.
Moreover, since the initial information about the system component’s lifetime is usually very limited (it is most often possible to obtain only an estimate of the first and second moments of the distribution), we concentrated on studying the sensitivity of any decision about PM strategy choice to the shape of their distributions.

3. Process J and the Main Result

3.1. Process J

For the stated problem solution, consider process J, defined by expression (4) as more detailed. First of all, we note that, by virtue of our assumptions, for any PM strategy l L , including running to the system failure for l = k , process J is the regenerative one, whose regenerative epochs S i ( l ) are the maintenance or repair completion times. For any regenerative process, the following ergodic theorem holds.
Theorem 1.
For any admissible strategy l L of a controllable (in regeneration epochs) regenerative process J = { J ( t ) , t 0 } with a finite expected regeneration period E [ T i ( l ) ] < (where T i ( l ) = S i ( l ) S i 1 ( l ) ) and any integrable function g, defined on the process set of states E, the following limit property holds
lim t 1 t 0 t g ( J ( u ) ) d u = 1 E 0 T 1 ( l ) E 0 0 T 1 ( l ) g ( J ( u ) ) d u ,
where E 0 means the expectation given initial state equals zero i = 0 .
Proof. 
For regeneration epochs S i ( l ) of process J, its renewal process will be denoted by
N ( l ) ( t ) = i 0 1 { S i ( l ) t }
and the left part of equality (7) is represented as
lim t 1 t 0 t g ( J ( u ) ) d u = = lim t N ( l ) ( t ) t 1 N ( l ) ( t ) 1 i N ( l ) ( t ) S i 1 ( l ) S i ( l ) g ( J ( u ) ) d u + 1 t S N ( l ) ( t ) ( l ) t g ( J ( u ) ) d u .
Given that, following our assumptions, the last term in this equality tends to zero, according to the limit theorem for the renewal process
lim t N ( l ) ( t ) t = 1 E 0 T 1 ( l ) ,
and the law of large numbers for independent identically distributed random variables
Y i ( l ) = S i 1 ( l ) S i ( l ) g ( J ( u ) ) d u
used under the sign of sum as follows
lim t 1 N ( l ) ( t ) 1 i N ( l ) ( t ) S i 1 ( l ) S i ( l ) g ( J ( u ) ) d u = E 0 0 T 1 ( l ) g ( J ( u ) ) d u
that ends the proof.   □

3.2. The Strategy Comparison

The theorem shows that we need to study the behavior of a process only in a separate regeneration period T i ( l ) = S i ( l ) S i 1 ( l ) (say the first one T 1 ( l ) = S 1 ( l ) ). Thus, considering function g as an indicator function of the state l , g ( j ) = 1 l ( j ) , the integral Y i ( l ) takes the value (since they are independent identically distributed random variables, we omit the lower index and place the upper index below)
Y i ( l ) Y 1 ( l ) = 0 S 1 ( l ) 1 l ( J ( u ) ) d u = inf { t : J ( t ) = l , t S 1 ( l ) } Y l ,
which equals the time of the l-th type maintenance beginning (the time for the l-th system component failure, the system lifetime for l = k ) in a separate regeneration period. We denote its mean value by M l = E [ Y l ] . Then the system availability for different PM strategies l L takes the form
K av . , l = M l E [ T l ] .
Hence, by virtue of the properties of regenerative processes, to calculate the availability K av . , l , we only need the mean value E [ T l ] of the regeneration period T l and the mean value M l = E [ Y l ] of the working time Y l in it. Since the regeneration period for any PM strategy l L is T l = Y l + B l , and it is assumed that the mean PM and repair times b l = E [ B l ] are known to a Decision Maker (DM), for the problem solution, it suffices to calculate the expectations M l = E [ Y l ] . Thus, from Theorem 1, it follows
Theorem 2.
The l-th strategy is preferable to the j-th one ( l j ) iff
b l b j < M l M j .
Proof. 
The comparison of the different PM strategies with respect to maximizing the system availability given by (14) leads to a comparison of the expressions
K av . , l = M l M l + b l
and because the l-th preventive maintenance strategy is preferable to the j-th one, if K av . , l > K av . , j from the inequality
M l M l + b l > M j M j + b j
it follows that the l-th strategy is preferable to the j-th one ( l j ) if and only if M l b j > M j b l , which in terms of dimensionless indexes, can be rewritten as (15), that ends the proof.    □

4. Distribution of PM Start and System Failure Times

4.1. Distribution of Order Statistics

It is evident that the l-th PM start time (including system failure for l = k ) is the time of the l-th component failure, in which the l-th order statistics will be denoted as X i = A ( i ) from n iid rvs A i ( i = 1 , 2 , , n ) with given cdf A ( t ) . Their distributions are very well studied (see, for example, [31]). In particular, it has been shown in this book that the joint pdf f l ( x 1 , x l ) of the l first-order statistics
X 1 X 2 X l
from a sample of n iid rv A 1 , A 2 , , A n with cdf A ( x ) and pdf a ( x ) in the domain 0 x 1 x 2 x l < equals
f l ( x 1 , x 2 , , x l ) = n ! ( n l ) ! a ( x 1 ) a ( x 2 ) a ( x l ) ( 1 A ( x l ) ) n l .
Remark 1.
From the other side, the individual l-th order statistics distribution does not depend on the positions of all its previous and equals (see also [31])
f l ( x ) = n ! ( n l ) ! ( l 1 ) ! A ( x ) l 1 ( 1 A l ) ) n l a ( x ) .
Therefore, it can be calculated by integrating the joint pdf f l ( x 1 , x l ) over all possible ( l 1 ) ! permutations of variables x i ,
f l ( x ) = 1 ( l 1 ) ! 0 x i x ( i = 1 , l 1 ¯ ) f l ( x 1 , x l 1 , x ) d x 1 d x l 1 = = n ! ( n l ) ! ( l 1 ) ! 0 x i x ( i = 1 , l 1 ¯ ) a ( x 1 ) a ( x 2 ) a ( x ) ( 1 A ( x ) ) n l = = n ! ( n l ) ! ( l 1 ) ! A ( x ) l 1 ( 1 A l ) ) n l a ( x ) .
However, if after the failure of one of the components, the load on all others changes, their residual lifetimes also change. Thus, in the next subsection, we consider the transformation of order statistics due to load changes. It will be used further to calculate the pdf of the order statistics in the case of their transformations.

4.2. Transformation of Order Statistics

For the stated problem’s solution, it is necessary to transform order statistics and calculate their distributions in accordance with the accepted rule for a component failure to influence the load on the surviving components. In accordance with this rule, after failure of the i-th component, the residual lifetimes of all survived components are reduced by multiplying by some weight coefficient w i .
Keep in mind that Y l ( l = 1 , k ¯ ) denotes the time to the l-th component failure (time to the l-th PM strategy beginning) or the system failure time for l = k in conditions of load changing after any of its components fail. We represent it in terms of its initial lifetime (order statistics X i ). It should be noted that for the considered system with dependent failure times, each subsequent l-th failure leads to compression of the intervals between successive j-th ( j > l ) and l-th failure times by w 1 w l times, Y j = Y l + w 1 w l ( X j X l ) . Thus, for the times Y l of the l-th state destination, the following generating recursive representation in terms of initial order statistics holds,
Y l + 1 = Y l + w 1 w l ( X l + 1 X l ) .
Based on this recursive relation, we prove the following explicit representation of the rvs Y l in terms of order statistics.
Lemma 1.
The time Y l to the l-th PM start (system failure time for l = k ) is a linear function of order statistics:
Y 1 = X 1 and for l = 2 , k ¯ Y l = ( 1 w 1 ) X 1 + w 1 ( 1 w 2 ) X 2 + + w 1 w l 2 ( 1 w l 1 ) X l 1 + + w 1 w l 1 X l .
Proof. 
To prove the lemma, we use mathematical induction. The initial component failure times coincide with the first-order statistics A ( 1 ) = X 1 . After the first failure at time X 1 , all other expected failure times equal
Y j = X 1 + w 1 ( X j X 1 ) = ( 1 w 1 ) X 1 + w 1 X j
and, therefore,
Y 2 = X 1 + w 1 ( X 2 X 1 ) = ( 1 w 1 ) X 1 + w 1 X 2 ,
which coincides with (23) for l = 2 and gives the possibility to begin the induction procedure. Thus, we suppose that equality (23) holds for any l < k , and check it for l + 1 . Using recursive Formula (22), representing the value Y l + 1 in terms of Y l , one can find
Y l + 1 = Y l + w 1 w l ( X l + 1 X l ) = = ( 1 w 1 ) X 1 + w 1 w l 2 ( 1 w l 1 ) X l 1 + w 1 w l 1 X l + + w 1 w l ( X l + 1 X l ) = = ( 1 w 1 ) X 1 + w 1 w l 2 ( 1 w l 1 ) X l 1 + w 1 w l 1 X l + + w 1 w l X l + 1 w 1 w l X l = = ( 1 w 1 ) X 1 + w 1 w l 1 ( 1 w l ) X l + w 1 w l X l + 1 ,
that ends the proof.    □

4.3. Distribution of the Random Variables Y l

Now we proceed to the calculation of the cdf F Y l ( t ) of the time Y l to the l-th component failure, taking into account the redistribution of the load on the components. Here we perform this using expression (23) for the time to the l-th component failure in terms of order statistics X i and Formula (19) for the joint distribution of the first l-order statistics. To simplify the representation of this cdf, let us introduce the following notations.
z 0 = t , and for ( i = 1 , l 1 ¯ ) ; z i = z i ( t ; x 1 , , x i ) = = t ( 1 w 1 ) x 1 w 1 ( 1 w 2 ) x 2 w 1 w i 1 x i 1 w 1 w i .
With these notations, the following theorem holds.
Theorem 3.
The distribution of the l-th component failure (time to the l-th PM beginning) for t 0 is
F Y l ( t ) = P { Y l t } = = n ! ( n k ) ! 0 z 0 a ( x 1 ) d x 1 x 1 z 1 a ( x 2 ) d x 2 x l 1 z l 1 a ( x l ) ( 1 A ( x l ) ) n l d x l .
Proof. 
Lemma 1 states that the time Y l of the l-th component failure is a linear function of the first l-order statistics (as shown by Formula (23)). Consequently, given the pdf f l ( x 1 , , x l ) of the l-th order statistics, we can derive the cdf F Y l ( t ) of rv Y l
F Y l ( t ) = P { Y l t } = = P { ( 1 w 1 ) X 1 + w 1 ( 1 w 2 ) X 2 + + + w 1 w l 2 ( 1 w l 1 ) X l 1 + w 1 w l 1 X l t } = = D ( t ; x 1 , , x l ) f l ( x 1 , x 2 , , x l ) d x 1 d x l ,
where the integration domain D ( t : x 1 , , x l ) is the intersection of two subsets of l-dimensional space, D ( t ; x 1 , , x l ) = D + ( l ) D w ( l ) ( t ) , where D + ( l ) is the first quadrant of the l-dimensional space: D + ( l ) = { 0 x 1 x l } and the hyperspace D w ( l ) ( t ) is
D w ( l ) ( t ) = { ( 1 w 1 ) x 1 + w 1 ( 1 w 2 ) x 2 + + w 1 w l 2 ( 1 w l 1 ) x l 1 + + w 1 w l 1 x l t } .
The multidimensional integral (29) can be transformed into an iterated one. Given 0 x 1 x 2 x l , the integration domain can be represented according to the following procedure. From inequality (30) for variable x l , it follows that
x l t w 1 x 1 w 1 ( 1 w 2 ) x 2 w 1 w i 1 x l 1 w 1 w l 1 = z l 1 ( t ; x 1 x l 1 ) .
Further, also given x l 1 x l from inequality (31), it follows that
x l 1 x l t w 1 x 1 w 1 ( 1 w 2 ) x 2 w 1 w l 1 x l 1 w 1 w l 1 .
From inequality (32), using algebraic manipulations, one can find
x l 1 t w 1 x 1 w 1 ( 1 w 2 ) x 2 w 1 w l 2 x l 2 w 1 w l 2 = = z l 2 ( t ; x 1 x l 2 ) .
In a similar way, we obtain the following inequality for the variable x 2
t ( 1 w 1 ) x 1 + w 1 ( 1 w 2 ) x 2 + w 1 w 2 x 3 ( 1 w 1 ) x 1 + w 1 ( 1 w 2 ) x 2 + w 1 w 2 x 2 = ( 1 w 1 ) x 1 + w 1 x 2 ,
from which it follows that
x 2 t ( 1 c 1 ) x 1 c 1 = z 2 ( t : x 1 ) .
Further, at last,
t ( 1 w 1 ) x 1 + w 1 x 1 = x 1 .
Consequently, 0 x 1 t . Therefore, the integration domain D ( t : x 1 , , x l ) in terms of notations (27) can be transformed as
D ( t : x 1 , , x l ) = { x i 1 x i z i ( t ; x 1 , x i 1 ) ( i = 1 , l 1 ¯ ) } .
Thus, using Formula (19) for pdf f l ( x 1 , , x l ) of the first l-order statistics and the integration domain (37), integral (29) for t 0 can be written as
F Y l ( t ) = n ! ( n l ) ! 0 t a ( x 1 ) d x 1 x 1 z 1 a ( x 2 ) d x 2 x l 1 z l 1 a ( x l ) ( 1 A ( x l ) ) n l d x l ,
which coincides with (28) and ends the proof.    □
Based on Theorem 3, the main system reliability characteristics can be obtained, such as:
(i)
its reliability function R ( t ) = 1 F Y k ( t ) ;
(ii)
its mean lifetime E [ Y k ] = 0 R ( t ) d t ;
(iii)
availability factor K av . , l for different PM strategies.

4.4. An Example: Exponential Distribution of Components Lifetimes

In a special case, when the systems’ components lifetimes A i ( i = 1 , n ¯ ) follow an exponential distribution with a parameter α ( E x p ( α )), the integral (28) can be calculated analytically. For this case, we introduce a different approach to the system’s lifetime distribution, which is derived from the memoryless property of an exponential distribution. Let T i be the time interval between the i 1 -th and i-th component’s failures, i = 1 , l 1 ¯ ( T 0 = 0 ). Then, by virtue of the memoryless property of the exponential distribution, the time to the k-th failure Y k is the sum
Y l = T 1 + T 2 + + T l ,
of l independent exponentially distributed rvs T i with mean values
E [ T 1 ] = 1 n α , E [ T i ] = w 1 w 2 w i 1 ( n i + 1 ) α = w ¯ i ( n i + 1 ) α , ( i = 2 , l ¯ ) .
This means that the parameters of their distribution are:
λ 1 = n α , λ i = ( n i + 1 ) α w 1 w 2 w i 1 = ( n i + 1 ) α w ¯ i , ( i = 2 , l ¯ ) ,
where, for purposes of simplicity, we introduce the notation:
w ¯ i = 1 , i = 1 , w 1 w i 1 , i = 2 , l ¯ .
The general case can be studied with the help of the moment-generating function (mgf) of the system lifetime, which, in this case, has the following form:
ϕ k ( s ) = E e s Y k = 1 i k E e s T i = 1 i k λ i s + λ i .
However, the cdf calculation from this expression is not a simple problem, and we limit ourselves to the simplest example of a k-out-of-n model with k = 2 and first calculate the cdf of Y 2 directly. In this case, we assume w 1 = w . Hence, following (19), the joint distribution of rv X 1 , X 2 is
f 2 ( x 1 , x 2 ) = n ! ( n 2 ) ! a ( x 1 ) a ( x 2 ) ( 1 A ( x 2 ) ) n 2 = n ! ( n 2 ) ! α 2 e α x 1 e ( n 1 ) α x 2 .
Calculate cdf F Y 2 ( t ) of the rv Y 2 = ( 1 w ) X 1 + w X 2 ,
F Y 2 ( t ) = P { ( 1 w ) X 1 + w X 2 < t } = P X 2 < t ( 1 w ) X 1 w = = n ( n 1 ) α 2 0 t e α x 1 d x 1 x 1 t ( 1 w ) x 1 w e ( n 1 ) α x 2 d x 2 = = 1 + n 1 n w ( n 1 ) e n α t n w n w ( n 1 ) e ( n 1 ) α w t ,
and, therefore, its pdf for t 0 is
f Y 2 ( t ) = n ( n 1 ) α n w ( n 1 ) e ( n 1 ) α w t e n α t .
It should be noted that this result is valid for w ( n 1 ) / n , in which case the distribution is a mixture of exponential distributions. The point w = ( n 1 ) / n is a singular point for which cdf of the rv Y 2 follows an Erlang distribution,
F Y 2 ( t ) = 1 e n α t n α t e n α t , t > 0 ,
with pdf
f Y 2 ( t ) = n 2 λ 2 t e n λ t , t > 0 .
Remark 2.
The singularity in the calculation of the cdf of the system lifetime arises because, for some special values of coefficient w i (here for w = ( n 1 ) / n ), the moment-generating function of the system’s lifetime has multiple roots that lead to changing the shape of the distribution.
Using another approach, the mgf of the system’s lifetime can be written as follows:
ϕ 2 ( s ) = n ( n 1 ) α 2 s 2 + ( 2 n 1 ) α s + n ( n 1 ) α 2 .
By decomposing this expression, one finds
ϕ 2 ( s ) = ( n 1 ) n α s + n α ( n 1 ) n α s + ( n 1 ) α ,
the inverse function, then, is calculated as
f 2 ( t ) = n ( n 1 ) α e ( n 1 ) α t e n α t ,
which coincides with the result obtained in (45) for w = 1 .
Apparently, analytic results are not to be expected in the general case. However, one can hope for numerical calculations of examples with different distributions of the random variable A i . Thus, based on Theorem 2, we propose a recursive procedure for calculating distributions F l ( t ) of random variable Y l ( l = 1 , k ¯ ) and their moments. Based on this calculation, the PM strategies comparison as advice to DM will be done. Below an Algorithm for the problem solution is proposed.

5. An Algorithm for the PM Strategies Comparison and Numerical Experiments

5.1. Algorithm

On the basis of the results of Section 4, the general procedure for the problem solution can be implemented with the following algorithm (Algorithm 1).
Algorithm 1 General algorithm for calculation of the reliability function.
         Beginning. Determine: Integers n , k , real b 0 , b l ( l = 1 , k ¯ ) ) , w l ( l = 1 , k ¯ ) , distribution A ( t ) of the system components’ lifetime, its pdf, expectation a and coefficient of variation v.
         Step 1. Calculate parameters of cdf A ( t ) in terms of its expectation a and coefficient of variation v.
         Step 2. Taking into account the joint pdf of the first l order statistics X 1 X 2 X l
f l ( x 1 , x 2 , , x l ) = n ! ( n l ) ! a ( x 1 ) a ( x 2 ) a ( x l ) ( 1 A ( x l ) ) n l ,
         calculate according to (28) the rv Y l cdf
F Y l ( t ) = n ! ( n l ) ! 0 z 0 a ( x 1 ) d x 1 x 1 z 1 a ( x 2 ) d x 2 x l 1 z l 1 a ( x l ) ( 1 A ( x l ) ) n l d x l ,
where the limits of integration are determined by the relation (27)
z 0 = t , and for ( i = 1 , l 1 ¯ ) ; z i = z i ( t ; x 1 , , x i ) = = t ( 1 w 1 ) x 1 w 1 ( 1 w 2 ) x 2 w 1 w i 1 x i 1 w 1 w i .
         Step 3. Calculate the function R l ( t ) = 1 F Y l ( t ) and the mean time to the PM strategy beginning (mean time to the l-th component failure)
M l = E [ Y l ] = 0 R l ( t ) d t ,
         and the availability factor
K av . , l = M l M l + b l .
         Step 4. Compare the PM strategies by the criterion (15)
b l b j < M l M j
and chose the best one.
         Stop.
Remark 3.
The algorithm might also be useful for other purposes, such as analysis of the sensitivity of the system reliability function and its characteristics to the shape of the lifetime distribution of the system components.

5.2. Numerical Experiments

Based on the algorithm, we consider some numerical experiments. For the numerical analysis, we use the four most popular distributions of non-negative random variables: (i) exponential ( E x p ( λ ) ), (ii) Gamma ( Γ ( λ , Θ ) ), (iii) Gnedenko–Weibull ( G W ( λ , Θ ) ), and (iv) lognormal ( L N ( μ , σ ) ) ). The corresponding rvs are denoted as E , Γ , G W , L N . All of these distributions (except for the exponential), for some of their parameters, became the so-called “heavy tailed ones”, which imposes certain restrictions on the computational procedures associated with their use.
The Gamma distribution with parameters ( λ > 0 , Θ > 0 ) has the pdf and cdf of the form
f Γ ( x ) = 1 Γ ( Θ ) λ Θ x Θ 1 e x λ , F Γ ( x ) = 1 Γ ( Θ ) λ Θ 0 x t Θ 1 e t λ d t ,
where Γ ( Θ ) is the Gamma function. Its expectation, variance, and coefficient of variation have the form
E [ Γ ] = λ Θ , Var [ Γ ] = λ 2 Θ , v Γ = 1 Θ .
Concerning the G W -distribution with parameters ( λ > 0 , Θ > 0 ) , it has the pdf and the cdf of the form
f G W ( x ) = Θ λ x λ θ 1 e ( x λ ) Θ , F G W ( t ) = 1 e t λ Θ ,
its coefficient of variation v depends only on parameter Θ
v 2 + 1 = Γ ( 1 + 2 Θ ) Γ 2 ( 1 + 1 Θ ) = 2 Θ Γ ( 2 Θ ) 1 Θ 2 Γ 2 ( 1 Θ ) = 2 Θ Γ ( 2 Θ ) Γ 2 ( ( 1 Θ ) ) .
For 0 < Θ 1 , for example Θ = 1 l it holds
v 2 + 1 = ( l + 1 ) ( 2 l 1 )
Specifically, for Θ = 0.5 , it gives v = 2 , for Θ = 1 3 it gives v 2 + 1 = 20 . Thus, the G W distribution for Θ < 1 has a heavy tail, which forces us to confine ourselves to studying the behavior of the model, in this case, by the range of values of coefficient of variation v 2 .
For Θ 1 , this distribution is a light-tailed one. In this case, putting Θ = 1 ϵ note that the behavior of its coefficient of variation depends on the behavior of the function Γ ( 2 ϵ ) ϵ Γ 2 ( ϵ ) around ϵ = 0 .
For Θ = 1 , both of the above distributions turn into an exponential.
For the L N ( μ , σ 2 ) distribution, its pdf and cdf are
f L N ( x ) = 1 x σ π exp 1 2 ln x μ σ 2 , F L N ( x ) = 1 2 1 + ln x μ σ 2 .
Its expectation, variation, and coefficient of variation are
E [ L N ] = e μ + σ 2 2 , Var [ L N ] = e σ 2 1 e μ + σ 2 2 , v [ L N ] = e σ 2 1 .
In the numerical experiments below, the parameters of the distributions were taken such that the mean operating time of the system components was equal to 1 (which, accordingly, meant the unit of time measurement), and another parameter provided the coefficient of variation in the domain 0.3 v 1.3 . Figure 2 presents a comparison of “2-strategy” and ”3-strategy” for different system components’ lifetime distributions and coefficients of variation. Solid lines present the calculations for w 1 = w 2 = 0.9 , and dashed lines demonstrate the results for w 1 = w 2 = 0.5 . Blue lines show the results for the Gamma distribution of the component’s lifetime, red is for the Gnedenko–Weibul distribution, green is for the lognormal one, and black presents the exponential distribution. For the exponential distribution, the analytical results coincide with the numerical calculations.
Let the ratio of PM time after two failures and after three failures be b 2 / b 3 = 0.75 and w 1 = w 2 = 0.9 . For the exponential distribution of the system components’ lifetime, the “3-strategy” of PM should be chosen. If the variation of system components’ lifetime is less than 0.6, the “2-strategy” is preferable to the “3-strategy” for Gamma, GW, and LN distributions. If the variation of the system components’ lifetime is greater than 0.8, the “3-strategy” is preferable to the “2-strategy” for Gamma, GW, and LN distributions. In the interval [ 0.6 , 0.8 ] , the choice of strategy depends on the lifetime distribution of system components.
We may conclude that for small values of w i and variation, the difference between Gamma and GW distributions is insignificant. For the LN distribution, the interval of variation prefers the “2-strategy” over the “3-strategy”.

5.3. An Example: Exponential Distribution

In the case of an exponential distribution based on Formula (42) for the system 2-out-of-6, the following can be written
ϕ 2 ( s ) = λ 1 s + λ 1 λ 2 s + λ 2 ,
where λ 1 = 6 α , λ 2 = 5 α w . Therefore,
ϕ 1 ( s ) = λ 1 λ 2 [ ( s + λ 1 ) ] ( s + λ 2 ) ] 2 ,
and
M 2 = ϕ 1 ( 0 ) = 1 λ 1 + 1 λ 2 = 1 6 α + w 5 α .
Analogously, for the 3-out-of-6 system it holds
ϕ 3 ( s ) = λ 1 s + λ 1 λ 2 s + λ 2 λ 3 s + λ 3 ,
where λ 1 , λ 2 are the same as before, and λ 3 = 4 α w 2 . By calculation of the derivative of mgf ϕ 3 ( s ) in points s = 0 and taking into account the values of parameters λ i , one can find
M 3 = ϕ 3 ( 0 ) = 1 6 α + w 5 α + w 2 4 α .
For the given α = 1 and w = 0.5 and w = 0.9 , these values coincide with the point in Figure 2.

6. Conclusions

The paper is devoted to the investigation of different PM strategies for the controllable k-out-of-n system based on its stated observation. It is supposed that after a failure of one of the system components, the load on all the others increases, which leads to a decrease in their residual lifetime. The novelty of the paper lies in the application of order statistics for modeling the residual lifetime of system components and developing a procedure for comparing PM strategies with respect to the availability maximization criterion, which is presented in the form of an algorithm.
A series of numerical experiments were performed for real examples modeled as a 3-out-of-6 system with exponential, Gamma, Gnedenko–Weibull, and Lognormal distributions and different coefficients of variation. The results revealed the sensitivity of the choice of PM strategy to the shape of lifetime’s distribution of the system components and the coefficient of variation. Nevertheless, it is possible to identify intervals for the coefficient of variation, where the choice of the preferable strategy is independent of the type of distribution.
The presented approach can be further expanded, involving more quality criteria, such as the productivity of the system or system service cost under different maintenance strategies, and a multi-criteria evaluation of the PM effectiveness is feasible as well. The calculations have been performed for the example of the 3-out-of-6 system, but the proposed procedure can be used for other systems. This article is theoretical in nature, and we invite persons with specific information to cooperate and further research the proposed direction.

Author Contributions

Conceptualization, methodology, formal analysis, validation, and writing—original draft preparation, V.R.; investigation, software and computer calculations, visualization, validation, and writing—original draft preparation, O.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All necessary data are contained in the article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PMPreventive Maintenance
DMDecision Maker
iidindependent and identically distributed
rvrandom variable

References

  1. Shepherd, D. k-out-of-n Systems. In Encyclopedia of Statistics in Quality and Reliability; John Wiley & Sons: Hoboken, NJ, USA, 2008; Available online: https://onlinelibrary.wiley.com/doi/pdf/10.1002/9780470061572.eqr342/ (accessed on 18 February 2022).
  2. Trivedi, K. Probability and Statistics with Reliability, Queuing and Computer Science Applications, 2nd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2016; Available online: https://onlinelibrary.wiley.com/doi/book/10.1002/9781119285441 (accessed on 18 February 2022).
  3. Chakravarthy, S.; Krishnamoorthy, A.; Ushakumari, P. A (k-out-of-n) reliability system with an unreliable server and Phase type repairs and services: The (N,T) policy. J. Appl. Math. Stoch. Anal. 2001, 14, 361–380. [Google Scholar] [CrossRef] [Green Version]
  4. Yuge, T.; Maruyama, M.; Yanagi, S. Reliability of a (k-out-of-n) system with common-cause failures using multivariate exponential distribution. Procedia Comput. Sci. 2016, 96, 968–976. [Google Scholar] [CrossRef] [Green Version]
  5. Zhang, T.; Xie, M.; Horigome, M. Availability and reliability of (k-out-of-(M+N)): Warm standby systems. Reliab. Eng. Syst. Saf. 2006, 91, 381–387. [Google Scholar] [CrossRef]
  6. Kuo, W.; Zuo, M. Optimal Reliability Modeling: Principles and Applications; John Wiley & Sons: Hoboken, NJ, USA, 2003. [Google Scholar]
  7. Wang, K.H.; Chen, W.L.; Yang, D.Y. Optimal management of the machine repair problem with working vacation: Newton’s method. Comput. Appl. Math. 2009, 233, 449–458. [Google Scholar] [CrossRef] [Green Version]
  8. Ke, J.C.; Hsu, Y.L.; Liu, T.H.; Zhang, Z.G. Computational analysis of machine repair problem with unreliable multi-repairmen. Comput. Oper. Res. 2013, 40, 848–855. [Google Scholar] [CrossRef]
  9. Wang, K.H.; Ke, J.B.; Ke, J.C. Profit analysis of the M/M/R machine repair problem with balking, reneging, and standby switching failures. Comput. Oper. Res. 2007, 34, 835–847. [Google Scholar] [CrossRef]
  10. Rykov, V.V.; Sukharev, M.G.; Itkin, V.Y. Investigations of k-out-of-n Systems Application Possibilities to Objects of Oil And Gas Industry. J. Mar. Sci. Eng. 2020, 8, 928. [Google Scholar] [CrossRef]
  11. Vishnevsky, V.M.; Kozyrev, D.V.; Rykov, V.V.; Nguyen, Z.F. Reliability Modeling of an unmanned high-altitude module of a tethered telecommunications platform. Inform. Technol. Comput. Syst. 2020, 4, 26–38. (In Russian) [Google Scholar] [CrossRef]
  12. Ying, C.; Qichao, M.; Ze, W.; Yingyi, L. Reliability analysis of k-out-of-n system with load-sharing and failure propagation effect. J. Syst. Eng. Electron. 2021, 32, 1221–1231. [Google Scholar] [CrossRef]
  13. Sharifi, M.; Taghipour, S.; Abhari, A. Inspection interval optimization for a k-out-of-n load sharing system under a hybrid mixed redundancy strategy. Reliab. Eng. Syst. Saf. 2021, 213, 107681. [Google Scholar] [CrossRef]
  14. Ahmadi, R. Reliability and maintenance modeling for a load-sharing k-out-of-n system subject to hidden failures. Comput. Ind. Eng. 2020, 150, 106894. [Google Scholar] [CrossRef]
  15. Nezakati, E.; Razmkhah, M. Reliability analysis of a load sharing k-out-of-n: F degradation system with dependent competing failures. Reliab. Eng. Syst. Saf. 2020, 203, 107076. [Google Scholar] [CrossRef]
  16. Rykov, V.; Ivanova, N.; Kozyrev, D.; Milovanova, T. On Reliability Function of a k-out-of-n System with Decreasing Residual Lifetime of Surviving Components after Their Failures. Mathematics 2022, 10, 4243. [Google Scholar] [CrossRef]
  17. Gertsbakh, I. Preventive Maintenance Models Based on the Lifetime Distribution. In Reliability Theory with Applications to Preventive Maintenance; Springer: Berlin/Heidelberg, Germany, 2000; pp. 67–106. [Google Scholar]
  18. Finkelstein, M.; Levitin, G. Preventive maintenance for homogeneous and heterogeneous systems. Appl. Stoch. Model. Bus. Ind. 2019, 35, 908–920. [Google Scholar] [CrossRef]
  19. Finkelstein, M.; Cha, J.H.; Levitin, G. On a new age-replacement policy for items with observed stochastic degradation. Qual. Reliab. Eng. Int. 2020, 36, 1132–1146. [Google Scholar] [CrossRef]
  20. Park, M.; Lee, J.; Kim, S. An optimal maintenance policy for k-out-of-n system without monitoring component. Qual. Technol. Qual. Manag. 2019, 16, 140–153. [Google Scholar] [CrossRef]
  21. Endharth, A.J.; Yun, W.Y.; Yamomoto, H. Preventive maintenance policy for linear-consequtive -k-out-of-n: F system. J. Oper. Res. Soc. Jpn. 2016, 59, 334–346. [Google Scholar]
  22. Hamdan, K.; Tavangar, M.; Asadi, M. Optimal preventive maintenance for repairable weighted k-out-of-n systems. Reliab. Eng. Syst. Safe 2021, 205, 107267. [Google Scholar] [CrossRef]
  23. Rezaei, E.; Jafary, B.; Fiondella, L. Optimal maintenance policies for linear consecutive k-out-of-n: F systems susceptible to dependent failures. Comput. Ind. Eng. 2022, 173, 108657. [Google Scholar] [CrossRef]
  24. Rykov, V.; Kochueva, O.; Farkhadov, M. Preventive Maintenance of a k-out-of-n System with Applications in Subsea Pipeline Monitoring. J. Mar. Sci. Eng. 2021, 9, 85. [Google Scholar] [CrossRef]
  25. Rykov, V.; Kochueva, O.; Rykov, Y. Preventive Maintenance of k-out-of-n System with Respect to Cost-Type Criterion. Mathematics 2021, 9, 2798. [Google Scholar] [CrossRef]
  26. Kala, Z. Quantification of Model Uncertainty Based on Variance and Entropy of Bernoulli Distribution. Mathematics 2022, 10, 3980. [Google Scholar] [CrossRef]
  27. Kala, Z. New Importance Measures Based on Failure Probability in Global Sensitivity Analysis of Reliability. Mathematics 2021, 9, 2425. [Google Scholar] [CrossRef]
  28. Rykov, V. On Reliability of Renewable Systems. In Reliability Engineering. Theory and Applications; Vonta, I., Ram, M., Eds.; CRC Press: Boca Raton, FL, USA, 2018; pp. 173–196. [Google Scholar]
  29. Rykov, V.; Kozyrev, D.; Filimonov, A.; Ivanova, N. On reliability function of a k-out-of-n system with general repair time distribution. Probab. Eng. Inf. Sci. 2020, 5, 1–18. [Google Scholar] [CrossRef]
  30. Vishnevsky, V.; Efrosinin, D.; Krishnamoorthy, A. Principles of Construction of Mobile and Stationary Tethered High-Altitude Unmanned Telecommunication Platforms of Long-Term Operation. In Communications in Computer and Information Science; Springer: Cham, Switzerland, 2018; Volume 919. [Google Scholar] [CrossRef]
  31. David, H.; Nagaraja, H. Order Statistics, 3rd ed.; John Wiley & Sons: New York, NY, USA, 2003. [Google Scholar]
Figure 1. Systems described with k-out-of-n models ((a) unmanned underwater vehicle, (b) copter).
Figure 1. Systems described with k-out-of-n models ((a) unmanned underwater vehicle, (b) copter).
Mathematics 11 00422 g001
Figure 2. Comparison of “2-strategy” and ”3-strategy” for different coefficients of variation.
Figure 2. Comparison of “2-strategy” and ”3-strategy” for different coefficients of variation.
Mathematics 11 00422 g002
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

Rykov, V.; Kochueva, O. Preventive Maintenance of k-out-of-n System with Dependent Failures. Mathematics 2023, 11, 422. https://doi.org/10.3390/math11020422

AMA Style

Rykov V, Kochueva O. Preventive Maintenance of k-out-of-n System with Dependent Failures. Mathematics. 2023; 11(2):422. https://doi.org/10.3390/math11020422

Chicago/Turabian Style

Rykov, Vladimir, and Olga Kochueva. 2023. "Preventive Maintenance of k-out-of-n System with Dependent Failures" Mathematics 11, no. 2: 422. https://doi.org/10.3390/math11020422

APA Style

Rykov, V., & Kochueva, O. (2023). Preventive Maintenance of k-out-of-n System with Dependent Failures. Mathematics, 11(2), 422. https://doi.org/10.3390/math11020422

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