Next Article in Journal
Speed of Convergence of Time Euler Schemes for a Stochastic 2D Boussinesq Model
Next Article in Special Issue
Preface to the Special Issue on Probability and Stochastic Processes with Applications to Communications, Systems and Networks
Previous Article in Journal
Improved G-Optimal Designs for Small Exact Response Surface Scenarios: Fast and Efficient Generation via Particle Swarm Optimization
Previous Article in Special Issue
Analytic and Computational Analysis of GI/Ma,b/c Queueing System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Reliability Function of a k-out-of-n System with Decreasing Residual Lifetime of Surviving Components after Their Failures

1
Department of Applied Probability and Informatics, Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya Str., 117198 Moscow, Russia
2
Department of Applied Mathematics and Computer Modelling, National University of Oil and Gas “Gubkin University”, 65 Leninsky Prospect, 119991 Moscow, Russia
3
Institute for Transmission Information Problems (Named after A.A. Kharkevich) RAS, Bolshoy Karetny, 19, GSP-4, 127051 Moscow, Russia
4
V.A. Trapeznikov Institute of Control Sciences of Russian Academy of Sciences, 65 Profsoyuznaya Str., 117997 Moscow, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(22), 4243; https://doi.org/10.3390/math10224243
Submission received: 25 October 2022 / Revised: 10 November 2022 / Accepted: 11 November 2022 / Published: 13 November 2022

Abstract

:
We consider the reliability function of a k-out-of-n system under conditions that failures of its components lead to an increase in the load on the remaining ones and, consequently, to a change in their residual lifetimes. Development of models able to take into account that failures of a system’s components lead to a decrease in the residual lifetime of the surviving ones is of crucial significance in the system reliability enhancement tasks. This paper proposes a novel approach based on the application of order statistics of the system’s components lifetime to model this situation. An algorithm for calculation of the system reliability function and two moments of its uptime has been developed. Numerical study includes sensitivity analysis for special cases of the considered model based on two real-world systems. The results obtained show the sensitivity of system’s reliability characteristics to the shape of lifetime distribution, as well as to the value of its coefficient of variation at a fixed mean.

1. Introduction and Motivation

Ensuring the reliability of systems, objects, and processes is one of the main goals in their creation and further operation. Redundancy serves this aim, and a k-out-of- n : F model is a very popular configuration for it. This is a model of a system that consists of n components in parallel that fails when at least k of them fail. Hereinafter, we will use this notation omitting the symbol “F”.
Due to the wide range of practical applications of k-out-of-n systems, many papers have been devoted to their study. The bibliography on the related topics is extensive (see Trivedi [1], Chakravarthy et al. [2] and the bibliography therein). For a brief overview of further investigations, see, for example, [3] by Rykov et al. An overview of recent publications on k-out-of-n multi-state systems can be found in [4]. Furthermore, the  k-out-of-n systems with several types of failure have been considered in [5,6]. In the 1980s in [7,8] for the investigation of heterogeneous systems, Ushakov proposed the method of Universal Generating functions. At present, it has become a very popular technique and has been used in different applications (see, for example, a monograph by Levitin [9] and the bibliography therein). Recently, in [10], Kala proposed new sensitivity measures for the system’s reliability function based on the entropy of its structural function. Engineering applications of this model to the study of real-world systems can be found in [11] for the reliability study of some structures in the oil and gas industry. In [12], the model is used for reliability analysis of a remote monitoring system of underwater sections of gas pipeline, and in [13] for the reliability study of a rotary-wing flight module of a high-altitude telecommunications platform.
Another interesting line of research within the framework of the problem of system reliability enhancement is the prediction of the remaining useful life (RUL), which is an indispensable indicator to measure the degradation process of system components. In [14], a novel adaptive approach based on Kalman filter and expectation maximum with Rauch–Tung–Striebel was proposed to solve the problem of the RUL prediction of lithium-ion battery which is critical for the normal operation of electric vehicles [15].
A data-driven RUL prediction approach based on deep learning was proposed in [16] and verified by two real-world datasets—the aircraft engines dataset and the actual milling machine dataset.
Recently, an interesting approach of stress–strength reliability characteristic study was proposed (see, especially, [17,18,19]). It is interesting to study this index with respect to its sensitivity to both stress and strength distribution. We do not touch this approach here, but it will be in our plans in future.
In paper [13], a wide range of issues was posed for the study of systems whose failure depend not only on the number of failed components, but also on their location in the system. Moreover, it is also very important to take into account that failures of system components lead to the increase in the load on the remaining ones. A simple load-sharing model, in which the lifetime is exponentially distributed and the load from the failed components is distributed proportionally among the survivors, is considered in [20] through the example of a 2-out-of-3 system. A load-sharing k-out-of- n : G system with identical components and arbitrary distribution of lifetime under the equal load-sharing rule in the context of semi-Markov embedded processes was studied in [21].
The study of a k-out-of-n model in which failed components do not affect the residual lifetime of surviving components, using order statistics, is considered in [22]. On the other hand, the increase in the load on working components after the stop of functioning of the failed ones can lead to the decrease in their residual lifetime. Such a problem has been studied in our previous papers [12,23]. In addition, in [24] this problem was modeled by the changing in components’ hazard rate function.
The application of order statistics to the study of k-out-of-n models is not new [1,25]. Previously, in [26], the so-called sequential order statistics (which is some extension of ordinary order statistics) were considered for the study of a k-out-of-n system, in which a failure of any component can affect other components, so that their basic failure rate is corrected in relation to the number of previous failures. A similar model of the impact of a component’s failure on the functioning of the survived ones has been developed for example in [27,28], where it was supposed that the failure of any component influences the others, so that their failure rate is adjusted with respect to the number of preceding failures.
However, the problem of system failure, associated with a change in the residual component lifetime, depending on the increase in load after the failure of any component, has not yet been solved. Thus, the present article is devoted to the solution of this problem. The novelty of this investigation consists of the following:
we perform the reliability study of a k-out-of-n system, whose component failures change residual lifetime of the other components;
in the current paper, despite the fact that order statistics have already been applied to the study of k-out-of-n system reliability characteristics, we propose a novel application of order statistics to study of the lifetimes of components and the whole system.
The paper is organized as follows. In the next section, the problem is set up, the main notations and some practical examples of k-out-of-n models are given. Then, in Section 3 the necessary preliminaries are introduced and in Section 4 the general procedure for the solution of the stated problem is proposed. The numerical study of different scenarios for the investigation of a 2-out-of-6 system is made in Section 5. In conclusion, directions for further research are outlined.

2. State of the Problem: Notations and Examples

2.1. Problem Setting

Usually, real-world redundant systems are constructed based on the same type of components. Thus, we consider a k-out-of-n system that consists of n identical components in parallel and fails if at least k of them fail. At that point, it is supposed that the failure of an i-th component for i < k leads to the increase of load on the others and therefore to the decrease of their residual lifetimes. It is modeled by multiplying the residual lifetime of the surviving components by some weighting factor c i < 1 , ( i = 1 , k 1 ¯ ) . We will consider the system operation up to its first failure.
In the present paper, the main reliability characteristics of such a system are studied, namely:
time T to the first failure of the system,
reliability function R ( t ) = P { T > t } of the system,
its two first moments,
high confidence quantiles;
sensitivity analysis of the system’s reliability function to the shapes of its components’ lifetime distribution.

2.2. Notations: Assumptions

To study the system, introduce the following notations:
  • P { · } , E [ · ] are symbols of probability and expectation;
  • A i : ( i = 1 , n ¯ ) is the series of components’ lifetimes, which are supposed to be independent identically distributed (iid) random variables (rv);
  • A ( t ) = P { A i t } is their common cumulative distribution function (cdf);
  • j is the system state, which means the number of failed components;
  • E = { j = { 0 , 1 , , k } } is the set of the system states.
Under the set of states E, define a stochastic process J = { J ( t ) : t 0 } by the expression
J ( t ) = j , if in time t the system is in the state j E
and denote by T and R ( t ) time to the first system failure and the reliability function, respectively,
T = inf { t : J ( t ) = k } , R ( t ) = P { T > t } .

2.3. Examples

As mentioned in the Introduction, k-out-of-n models have a wide sphere of applications (see [1] and others), including the study of energy (see [11,12]), and  telecommunication [13] problems. Let us focus on two examples of applying the k-out-of-n model. In the numerical analysis, we will use these examples for the special case of n = 6 , k = 2 .

2.3.1. A Flight Module of a Tethered High-Altitude Telecommunication Platform

As an application example of the proposed k-out-of-n model, consider the model of a multi-copter flight module, which is part of the tethered high-altitude telecommunications platform [13]. The main area of its application is solving problems related to the long-term operation (tens of hours) without lowering the unmanned flight module to the ground. Therefore, unlike autonomous Unmanned Aerial Vehicles (UAVs) reliability parameters are of crucial importance for the tethered UAV-based high-altitude platforms.
A multi-rotor UAV is a system consisting of n rotors arranged uniformly in a circle and pairwise symmetrically with respect to the center of the circle [29]. The multi-copter may malfunction due to the failure of the propeller engines. There are various modifications of multi-rotor UAVs. The most common architectures are quad-, hexa-, and octocopters. The higher the redundancy ratio, the higher the reliability of the system. Therefore, in practice, flight modules with 6 or 8 rotors are most often used. In this example, we consider a hexacopter as a hot standby system consisting of n = 6 components (rotors) that work and fail independently of each other (see Figure 1).
If the location of the failing components is not taken into account, this system fails when k = 2 out of 6 rotors fail.For practical use, various reliability characteristics of such a system, including those considered in the general model, are of interest.

2.3.2. An Automated System for Remote Monitoring of a Sub-Sea Pipeline

As another application example of the k-out-of-n model, we consider an automated system for remote monitoring of a sub-sea pipeline. This system has been considered in [12], where its description has been given in details. One of the main parts of this system is an Unmanned Underwater Vehicle (UUV), the structure of which is illustrated in Figure 2.
The UUV consist of 6 motors, indicated by numbers 1–6, which allow it to rise, fall and move in various directions, including along the pipeline. The UUV is equipped with various devices, indicated by numbers 7–15, for receiving and transmitting information about the state of the pipe. In paper [12] the reliability function of this model has been studied in two scenarios:
(1)
in the case, when the system’s failure depends only on the number of its failed components. At that point, it is assumed that the device can perform its functions until at least 3 of its engines fail;
(2)
in the case, when the system’s failure depends also on the position of the failed components in the system. At that point, the UUV can perform its functions as long as at least two engines located on opposite sides, or any three engines are operational. Therefore, it could be considered to be a combination of 3 + 1 -out-of- 6 : F and 5-out-of- 6 : F systems. For such a system, the  special notation such as ( 5 , 3 + 1 ) -out-of- 6 : F system was used.
However, the influence of the number of failed components on the residual lifetimes of the survived ones was not taken into account earlier. In the current paper, this model has been studied under the condition that failed components reduce the residual lifetime of surviving system’s components.

3. Distribution of the System’s Time to Failure

3.1. Preliminaries

It is evident that if a k-out-of-n system’s failure depends only on several its failed components, it coincides with the k-th order statistic from n iid rv A i ( i = 1 , n ¯ ) with a given cdf A ( t ) . For simplicity, further we will denote order statistics A ( 1 ) A ( k ) A ( n ) of iid rv A i ( i = 1 , n ¯ ) by X i , i.e., X i = A ( i ) and X 1 X k X n . Distributions of order statistics are well studied (see, for example, [30]), where it was shown that the joint probability density function (pdf) f n ( x 1 , x n ) of all order statistics X 1 X 2 X n from n iid rv A 1 , A 2 , , A n with a given pdf a ( x ) has the following form:
f n ( x 1 , x 2 , , x n ) = n ! a ( x 1 ) a ( x 2 ) a ( x n ) ( x 1 x 2 x n ) .
By integration of this pdf with respect to last n k variables one can simply find the joint pdf f k ( x 1 , x k ) of the first k order statistics X 1 X 2 X k from the n iid rv A i ( i = 1 , n ¯ ) in the domain x 1 x 2 x k in the form
f k ( x 1 , x 2 , , x k ) = n ! ( n k ) ! a ( x 1 ) a ( x 2 ) a ( x k ) ( 1 A ( x k ) ) n k .
However, if a failure of one of the system’s components leads to the change in the residual lifetimes of all survived components, then their distributions are also changed.

3.2. Transformation of Order Statistics

Following the proposed model of the influence of components’ failures on the residual lifetime of survivors, they are reduced by multiplying by some constant c i depending on the number of failed components. Denote by Y i ( i = 1 , k ¯ ) the time of an i-th component failure under the conditions of increasing the load on survived components. To  simplify the representation of these values in terms of order statistics X 1 X 2 X n , we introduce the following notations,
C 1 = ( 1 c 1 ) , C 2 = c 1 ( 1 c 2 ) , , C k 1 = c 1 c k 2 ( 1 c k 1 ) , C k = c 1 c k 1 .
In these notations, the following theorem holds.
Theorem 1. 
The time to the considered system failure Y k is a linear function of order statistics of the following form:
Y k = C 1 X 1 + C 2 X 2 + + C k 1 X k 1 + C k X k .
Proof. 
To calculate the time of the system failure, we slightly expand the problem statement and calculate the successive time moments Y i ( i = 1 , k ¯ ) of failures of the system’s components under conditions of increasing load on the surviving components. To do that, we use a recursive procedure and denote by X i ( j ) the expected time moment of the i-th failure after the failure of the j-th component ( i > j ) .
Thus, to start the induction, we have X i ( 0 ) = X i . After the first failure of a component at time Y 1 = X 1 ( 0 ) = X 1 all residual lifetimes of surviving components that equal X i X 1 for i > 1 decrease by a factor of c 1 , and therefore the expected failure times X i ( 1 ) for i > 1 take the form
X i ( 1 ) = X 1 ( 0 ) + c 1 ( X i ( 0 ) X 1 ( 0 ) ) = ( 1 c 1 ) X 1 + c 1 X i , i = 2 , k ¯ .
Therefore Y 2 = X 2 ( 1 ) = ( 1 c 1 ) X 1 + c 1 X 2 .
Similarly, after the j-th failure at time Y j = X j ( j 1 ) , the residual lifetimes X i ( j 1 ) X j ( j 1 ) of all surviving components for all i > j decrease by a factor c j , 0 < c j < 1 ( j = 1 , k ¯ ) and the expected failure times of components take the following form:
X i ( j ) = X i ( j 1 ) , i j , X i ( j ) = X j ( j 1 ) + c j ( X i ( j 1 ) X j ( j 1 ) ) = ( 1 c j ) X j ( j 1 ) + c j X i ( j 1 ) , i > j .
Thus, the expected failure times of the system components Y j ( j = 1 , k ¯ ) under conditions of load redistribution equal to Y j = X j ( j 1 ) ( j = 1 , k ¯ ) . Expressing X i ( j ) in terms of the original order statistics, we obtain the following expression for i > j :
X i ( j ) = ( 1 c 1 ) X 1 + c 1 ( 1 c 2 ) X 2 + c 1 c 2 ( 1 c 3 ) X 3 + + c 1 c j 1 ( 1 c j ) X j + c 1 c j X i = l = 1 j c 1 c l 1 ( 1 c l ) X l + c 1 c j X i .
Supposing that the last expression is true for a given j check it for all i > j :
X i ( j + 1 ) = ( 1 c 1 ) X 1 + c 1 ( 1 c 2 ) X 2 + c 1 c 2 ( 1 c 3 ) X 3 + + c 1 · · c j 1 ( 1 c j ) X j + c 1 · · c j ( 1 c j + 1 ) X j + 1 + c 1 · · c j c j + 1 X i = l = 1 j c 1 c l 1 ( 1 c l ) X l + c 1 c j ( 1 c j + 1 ) X j + 1 + c 1 c j + 1 X i = l = 1 j + 1 c 1 c l 1 ( 1 c l ) X l + c 1 c j + 1 X i .
Hence, by the principle of mathematical induction, the equality (3) holds for any j. In terms of the original order statistics X i ( i = 1 , k ¯ ) , we obtain for all j = 1 , k ¯ :
Y j = X j ( j 1 ) = ( 1 c 1 ) X 1 + c 1 ( 1 c 2 ) X 2 + + c 1 c j 2 ( 1 c j 1 ) X j 1 + c 1 c j 1 X j ,
which, using the notation introduced earlier, leads to (2) for j = k , which completes the proof.    □

3.3. Distribution of the System Failure Time

Now move on to the calculation of the cdf F Y k ( y ) of the system’s time to failure Y k under the condition of redistribution of the load on the components. We will do that by taking into account expression (2) for the time of the system failure in terms of order statistics X i and using Formula (1) for the joint distribution of the first k order statistics.
To simplify the representation of this cdf we introduce the following notations,
z 0 = y , z i = z i ( y ; x 1 , , x i ) = y C 1 x 1 + C 2 x 2 C i x i C i + 1 ( i = 1 , k 1 ¯ ) .
With these notations the following theorem holds.
Theorem 2. 
The distribution of the system’s time to failure for y > 0 is
F Y k ( y ) = P { Y k < y } = n ! ( n k ) ! 0 z 0 a ( x 1 ) d x 1 x 1 z 1 a ( x 2 ) d x 2 x k 1 z k 1 a ( x k ) ( 1 A ( x k ) ) n k d x k .
Proof. 
According to Theorem 1 (see Formula (2)) the time Y k of the system failure is the linear function of the first k order statistics
Y k = C 1 X 1 + C 2 X 2 + + C k 1 X k 1 + C k X k .
Therefore, for cdf F Y k ( y ) of rv Y k in terms of pdf f k ( x 1 , , x k ) of the first k order statistics we obtain
F Y k ( y ) = P { Y k < y } = P { C 1 X 1 + C 2 X 2 + + C k 1 X k 1 + C k X k < y } = D ( x 1 , , x k ; y ) f k ( x 1 , x 2 , , x k ) d x 1 d x k ,
where the integration domain is
D ( x 1 , , x k ; y ) = { 0 x 1 x k , C 1 x 1 + C 2 x 2 + C 3 x 3 + + C k 1 x k 1 + C k x k y } .
Let us represent this multidimensional integral as an iterated one. Taking into account that x 1 x 2 x k , the integration domain can be transformed in the following way. For the last variable x k from the inequality
C 1 x 1 + C 2 x 2 + C 3 x 3 + + C k 1 x k 1 + C k x k y ,
it follows that
x k y C 1 x 1 C 2 x 2 C 3 x 3 C k 1 x k 1 c 1 c k 1 = z k 1 ( y ; x 1 x k 1 ) .
Furthermore, taking into account that x k 1 x k , from the last inequality, it follows that
x k 1 x k y C 1 x 1 C 2 x 2 C 3 x 3 C k 1 x k 1 c 1 c k 1 .
From this inequality with the simple algebra one can find
x k 1 y C 1 x 1 C 2 x 2 C k 2 x k 2 c 1 c k 2 = z k 2 ( y ; x 1 , x k 2 ) .
Following in the same way we obtain for variable x 2 the inequality
y ( 1 c 1 ) x 1 + c 1 ( 1 c 2 ) x 2 + c 1 c 2 x 3 ( 1 c 1 ) x 1 + c 1 ( 1 c 2 ) x 2 + c 1 c 2 x 2 = ( 1 c 1 ) x 1 + c 1 x 2 ,
from which it follows that
x 2 y ( 1 c 1 ) x 1 c 1 ,
and, at last,
y ( 1 c 1 ) x 1 + c 1 x 1 = x 1 .
It means that 0 x 1 y . This argumentation shows that the integration domain D ( x 1 , , x k ; y ) in terms of notations (5) can be represented as
D ( x 1 , , x k ; y ) = { x i 1 x i z i ( y ; x 1 , x i 1 ) ( i = 1 , k ¯ ) } .
Thus, using Formula (1) for pdf f k ( x 1 , , x k ) for the first k order statistics and the above form of the integration domain, we can rewrite integral (7) for y 0 as
F Y k ( y ) = n ! ( n k ) ! 0 y a ( x 1 ) d x 1 x 1 z 1 a ( x 2 ) d x 2 x k 1 z k 1 a ( x k ) ( 1 A ( x k ) ) n k d x k ,
that ends the proof.    □
As a consequence of the theorem, the main system reliability characteristics can be calculated.
Remark 1. 
Based on the distribution of the system’s time to failure, any other system’s reliability characteristics can be calculated, such as:
its reliability function R ( y ) = 1 F Y ( y ) ;
its mean lifetime E [ T ] = 0 R ( t ) d t ;
its lifetime variation var [ T ] .

3.4. A Special Case: Exponential Distribution

In a special case, when the system components’ A i ( i = 1 , n ¯ ) lifetimes have exponential ( E x p ) distribution with a parameter α the integral (6) can be calculated analytically, but the calculations are rather cumbersome. We show it for the given value of k = 2 . But for exponential distribution of the system components’ lifetime, we propose another approach for the system lifetime distribution. It is based on the memoryless property of any exponentially distributed rv.
Denote by T i the time interval between i 1 -th and i-th components failures, i = 1 , k 1 ¯ ( T 0 = 0 ). Then due to the memoryless property of the exponential distribution the time to the k-th failure Y k is the sum
Y k = T 1 + T 2 + + T k ,
of k independent exponentially distributed rv T i with parameters
λ 1 = n α , λ i = c 1 c 2 c i 1 ( n i + 1 ) α = c ¯ i ( n i + 1 ) α , i = 2 , k ¯ ,
where for simplicity additional notations are used:
c ¯ i = 1 , i = 1 , c 1 c i 1 , i = 2 , k ¯ .
The moment generating function (mgf) of the system’s lifetime 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 c ¯ i λ i s + c ¯ i λ i .
To apply the above theorem and the proposed approach, let us consider the simplest example of a k-out-of-n model with k = 2 . In this case, suppose c 1 = c . Thus, according to (1) 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 ( y ) of the rv Y 2 = ( 1 c ) X ( 1 ) + c X ( 2 ) ,
F Y 2 ( y ) = P { ( 1 c ) X ( 1 ) + c X ( 2 ) < y } = P X ( 2 ) < y ( 1 c ) X ( 1 ) c = n ( n 1 ) α 2 0 y e α x 1 d x 1 x 1 y ( 1 c ) x 1 c e ( n 1 ) α x 2 d x 2 = 1 + n 1 n c ( n 1 ) e n α y n c n c ( n 1 ) e ( n 1 ) α c y ,
and therefore its pdf for y 0 is
f Y 2 ( y ) = n ( n 1 ) α n c ( n 1 ) e ( n 1 ) α c y e n α y .
Please note that this result holds for c ( n 1 ) / n and in this case the distribution is a mixture of exponential distributions. The point c = ( n 1 ) / n is a singular point for which cdf of the rv Y 2 is the Erlang distribution,
F Y 2 ( y ) = 1 e n α y n α y e n α y , y > 0 ,
with pdf
p Y 2 ( y ) = n 2 λ 2 y e n λ y , y > 0 .
Remark 2. 
The singularity in the calculation of the cdf of the system’s lifetime arises because for some special values of the coefficient c i (here for c = ( n 1 ) / n ) the moment generating function of the system’s lifetime has multiple roots that leads to changing of the shape of distribution.
With the help of another approach one can find mgf of the system’s lifetime in the following form:
ϕ 2 ( s ) = n ( n 1 ) α 2 s 2 + ( 2 n 1 ) α s + n ( n 1 ) α 2 .
By expanding this expression into simple fractions, we find
ϕ 2 ( s ) = n ( n 1 ) α s + n α n ( n 1 ) α s + ( n 1 ) α ,
then, by calculating the inverse function, we obtain
f 2 ( y ) = n ( n 1 ) α e ( n 1 ) α y e n α y ,
which is the same as the result above for c = 1 .
The analytical calculations of the reliability characteristics are not always possible. Nevertheless, their numerical analysis in the wide domain of initial data is possible. Therefore, in the next section a procedure for the numerical calculation of different reliability characteristics of the considered system will be proposed. Furthermore, in Section 5 this procedure will be used for the numerical analysis of the model with some examples.

4. The General Calculation Procedure of the System Reliability Characteristics and Numerical Experiments

Based on the results of the previous section, the general procedure for the problem solution can be implemented with the help of the following algorithm (Algorithm 1).
Algorithm 1: General algorithm for calculation of reliability function
Beginning. Determine: Integers n , k , real c i ( i = 1 , k ¯ ) , distribution A ( t ) of the system components’ lifetime and its pdf.
Step 1. Taking into account that the system’s failure moment Y k according to Formula (2) equals
Y k = C 1 X 1 + C 2 X 2 + + C k 1 X k 1 + C k X k ,
calculate the following,
C i = 1 c i , i = 1 , c 1 c i 1 ( 1 c i ) i = 2 , k 1 ¯ , c 1 c k 1 i = k .
Step 2. Taking into account that according to Formula (1), the joint distribution density of first k order statistics X 1 X 2 X k holds
f X 1 X 2 X k ( x 1 , x 2 , , x k ) = n ! ( n k ) ! a ( x 1 ) a ( x 2 ) a ( x k ) ( 1 A ( x k ) ) n k ,
with which following to (6) calculate the reliability function
R ( y ) = 1 F Y k ( y ) = 1 n ! ( n k ) ! 0 y a ( x 1 ) d x 1 x 1 z 1 a ( x 2 ) d x 2 x k 1 z k 1 a ( x k ) ( 1 A ( x k ) ) n k d x k ,
where the limits of integration are determined by the relation (5)
z 0 = y , z i = z i ( y ; x 1 , , x i ) = y C 1 x 1 + C 2 x 2 C i x i c 1 c 2 c i ( i = 1 , k 1 ¯ ) .
Find the values of the constants c i (singular points at which the denominator of the cdf F Y k ( y ) turns into 0) for which the cdf changes its appearance.
Step 3. From the system reliability function R ( y ) , calculate
– mean time to the system failure
μ T = E [ Y k ] = 0 R ( y ) d y ;
– its variance
σ T 2 = Var [ Y k ] = 0 ( y μ T ) 2 f ( y ) d y , where f ( y ) = d d y F Y k ( y ) ,
and coefficient of variation
v = σ μ .
Stop.
Remark 3. 
The algorithm can also be used to solve other different problems, for example, to analyze the sensitivity of the system’s reliability function and its characteristics to the shape of the lifetime distribution of the system’s components.
Furthermore, the algorithm will be applied to some examples.

5. Numerical Experiments: 2-out-of-6 System

According to Algorithm 1, we calculate the reliability function of a 2-out-of-6 system. Since such a system fails due to the failure of two components, we have only one constant that defines the decreasing residual lifetime of surviving components. Therefore, hereafter, we suppose c 1 = c . Consider the Gnedenko–Weibull ( G W ) distribution as the lifetime distribution of the system’s components, A ( t ) G W θ , a Γ ( 1 + θ 1 ) , with the corresponding cdf
A ( t ) = 1 exp t Γ ( 1 + θ 1 ) a θ , t > 0 ,
where
  • a is a fixed mean components’ lifetime,
  • θ is the shape parameter of G W distribution calculated based on the preset value of the coefficient of variation,
  • v = σ a = a 1 · Γ ( 1 + 2 · θ 1 ) Γ ( 1 + θ 1 ) 2 1 is the coefficient of variation,
  • σ is the standard deviation.
Additionally, consider the Erlang ( E r l ) distribution, A ( t ) E r l l , θ with pdf
a ( y ) = θ l Γ ( l ) y l 1 e θ y , y > 0 .
In this case, the distribution’s parameters can be represented via the corresponding mean a and coefficient of variation v as follows,
l = v 2 , θ = ( a v 2 ) 1 .
For numerical experiments, we consider the reliability function and its characteristics of a 2-out-of-6 system for given distributions with a fixed mean a and different values of v. Thus, we can analyze the influence of the coefficient of variation of the repair time on the reliability characteristics of the system. In other words, investigate its sensitivity.
Suppose that the mean lifetime of the component a = 1 . If θ = 1 , G W and E r l distributions transform into the exponential one with the mean time a and the coefficient of variation v = 1 . In this case, its reliability function is
R ( t ) = 5 e 6 t 6 c · e 5 t c 5 6 c .
From Formula (8) it is clear that c = 5 6 leads to changing of the shape of distribution.
Since calculating the coefficient θ for G W through the value of v is quite difficult, we define the parameter θ so that v 0.5 . Moreover, if θ of G W takes non-integer values, it is not always possible to obtain a closed-form reliability function R ( t ) according to Algorithm 1 (the integrand takes a complex form). Therefore, define θ = 2 , then, the coefficient of variation v = 0.5227 . For E r l distribution, suppose that v = 0.5 , which leads to θ = 4 .
Suppose c = 0.1 ; 0.5 ; 1 . Figure 3 illustrates the reliability function of the 2-out-of-6 system for different distributions, as well as c and v. Here, solid line means v = 1 and reliability function (8), dashed one is for G W with v = 0.5227 and dash-dotted is for E r l with v = 0.5 . The legend of the figure denotes the color of line for different c.
The figure shows that the higher reliability coincides with the lower value of v. The case c = 1 means the absence of load from the failed components to the surviving ones, thus this case corresponds to the highest reliability for different v compared to the values c < 1 . Moreover, dependence of the reliability function curve on the shape of lifetime distribution is observed. On a small interval y, the system reliability as A E r l is higher than as A G W for each c, despite close values v. This may indicate the sensitivity of the reliability function not only to the shape of the lifetime distribution, but also to the corresponding value of the coefficient of variation.
According to the algorithm, we calculate other reliability characteristics of the 2-out-of-6 system (Table 1 and Table 2). These characteristics correspond to the system’s reliability behavior, shown in Figure 3. The lower value of v leads to the higher value of the system lifetime expectation E [ Y 2 ] , and the lower value of c leads to the lower value of E [ Y 2 ] . Moreover, as v 0.5 the relative error between the considered distributions is 14.11 % for c = 0.1 , 7.98 % for c = 0.5 and 3.86 % for c = 1 .
To distinguish coefficients of variation of the components and the whole system, denote them as v c o m p and v s y s , respectively. Thus, Table 2 shows the following. With a decrease in c, the coefficient of variation of the system v s y s grows and tends to the value of the coefficient of variation of each system component v c o m p . The increasing v c o m p leads to the increasing v s y s for all distributions and c. Thus, the coefficient of variation of the system v s y s confirms that as c tends to 0 and v c o m p tends to 1, variability with respect to the average lifetime of the system E [ Y 2 ] grows.
Table 1 showed that with increasing c and v 0.5 , the mean system lifetime E [ Y 2 ] is very close. However, Figure 3 shows that over a small interval y with these c and v, the reliability of the system has significant differences. This leads to the study the quantiles of the system reliability. This measure shows how long the system will be reliable with a fixed probability. The quantiles q γ = R 1 ( γ ) of the reliability function are presented in Figure 4, Figure 5 and Figure 6. In all cases, red bullets correspond to γ = 0.99 , whereas black bullets correspond to γ = 0.9 .
All the values for quantiles γ = 0.999 ; 0.99 ; 0.9 are presented in Table 3 for different distributions. The values in the table show that for the presented quantiles q γ , the shape of the lifetime distribution of the system’s components as well as its coefficient of variation play a critical role on the system’s reliability. Therefore, for example, as c = 0.1 and A E r l a given reliability level 0.9 will last about 8 times longer than for c = 0.1 and A E x p . At that for q 0.999 , the difference for similar case is almost 40 times. As the coefficient c increases, this difference decreases for all values of the quantiles and lifetime distributions of the components. As c = 1 this difference is reduced by about two times. Thus, even as c = 1 , which defines no changing in components’ residual lifetimes, the influence not only of the lifetime distribution of the components but also its coefficient of variation on the reliability of the system is huge. This once again confirms the sensitivity of the reliability characteristics of the k-out-of-n system to the shape of the lifetime distribution and the coefficient of variation of system’s components.

6. Conclusions and the Further Investigations

The reliability function of a new k-out-of- n : F model is investigated, under the new assumptions that the failures of its components lead to the increase in the load on the remaining ones and, consequently, to the change in their residual lifetimes. To model the situation, we proposed a novel approach based on the transformation of the order statistics of the system components’ lifetimes, which is the new field of application of order statistics. An algorithm for calculation of the system’s reliability function and its moments has been developed. Numerical experiments for the special case of the considered model based on the real-world systems have been carried out. The experiments show an essential sensitivity of the model reliability function and its moments to the shapes of the lifetime distributions of the system’s components and their coefficient of variation.
Furthermore, it is proposed we extend this approach to the investigation of stationary characteristics of the model and consider its preventive maintenance, aiming to improve its reliability characteristics.

Author Contributions

Conceptualization, writing—original draft preparation, supervision, project administration, V.R.; validation, investigation, visualization, N.I.; writing—review and editing, methodology, D.K.; data curation, software, T.M. All authors have read and agreed to the published version of the manuscript.

Funding

This paper has been supported by the RUDN University Strategic Academic Leadership Program (recipients V.R, supervision and problem setting, N.I., visualization, D.K. writing—review and editing, T.M., analytic results). This paper has been partially funded by RFBR according to the research projects No.20-01-00575A (recipients V.R., conceptualization, and N.I., formal analysis) and RSF according to the research projects No.22-49-02023 (recipient N.I., validation, D.K. review and analytic results).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors express their gratitude to the Referees for the valuable suggestions, which improved the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
iidindependent and identically distributed
rvrandom variable
cdfcumulative distribution function
pdfprobability density function
UAVUnmanned Aerial Vehicle
UUVUnmanned Underwater Vehicle
mgfmoment generating function
E x p Exponential distribution
G W Gnedenko–Weibull distribution
E r l Erlang distribution

References

  1. Trivedi, K.S. Probability and Statistics with Reliability, Queuing and Computer Science Applications, 2nd ed.; John Wiley & Sons: New York, NY, USA, 2016. [Google Scholar] [CrossRef]
  2. Chakravarthy, S.R.; Krishnamoorthy, A.; Ushakumari, P.V. 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]
  3. 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, 35, 885–902. [Google Scholar] [CrossRef]
  4. Pascual-Ortigosa, P.; Sáenz-de-Cabezón, E. Algebraic Analysis of Variants of Multi-State k-out-of-n Systems. Mathematics 2021, 9, 2042. [Google Scholar] [CrossRef]
  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. Gertsbakh, I.; Shpungin, Y. Reliability Of Heterogeneous ((k, r)-out-of-(n, m)) System. RTA 2016, 3, 8–10. [Google Scholar]
  7. Ushakov, I. A universal generating function. Sov. J. Comput. Syst. Sci. 1986, 24, 37–49. [Google Scholar]
  8. Ushakov, I. Optimal standby problem and a universal generating function. Sov. J. Comput. Syst. Sci. 1987, 25, 61–73. [Google Scholar]
  9. Levitin, G. The Universal Generating Function in Reliability Analysis and Optimization; Springer Series in Reliability Engineering; Springer: London, UK, 2005. [Google Scholar] [CrossRef]
  10. Kala, Z. New Importance Measures Based on Failure Probability in Global Sensitivity Analysis of Reliability. Mathematics 2021, 9, 2425. [Google Scholar] [CrossRef]
  11. Rykov, V.; Sukharev, M.; Itkin, V. Investigations of the Potential Application of k-out-of-n Systems in Oil and Gas Industry Objects. J. Mar. Sci. Eng. 2020, 8, 928. [Google Scholar] [CrossRef]
  12. 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]
  13. Vishnevsky, V.M.; Kozyrev, D.V.; Rykov, V.V.; Nguyen, D.P. Reliability modeling of an unmanned high-altitude module of a tethered telecommunication platform. Inf. Technol. Comput. Syst. 2020, 4, 26–36. [Google Scholar] [CrossRef]
  14. Zhang, J.; Jiang, Y.; Li, X.; Huo, M.; Luo, H.; Yin, S. An adaptive remaining useful life prediction approach for single battery with unlabeled small sample data and parameter uncertainty. Reliab. Eng. Syst. Saf. 2022, 222, 108357. [Google Scholar] [CrossRef]
  15. Zhang, J.; Jiang, Y.; Li, X.; Luo, H.; Yin, S.; Kaynak, O. Remaining Useful Life Prediction of Lithium-Ion Battery with Adaptive Noise Estimation and Capacity Regeneration Detection. IEEE/ASME Trans. Mechatron. 2022, 1–12. [Google Scholar] [CrossRef]
  16. Zhang, J.; Jiang, Y.; Wu, S.; Li, X.; Luo, H.; Yin, S. Prediction of remaining useful life based on bidirectional gated recurrent unit with temporal self-attention mechanism. Reliab. Eng. Syst. Saf. 2022, 221, 108297. [Google Scholar] [CrossRef]
  17. Eryilmaz, S. Phase type stress-strength models with reliability applications. Commun. Stat.—Simul. Comput. 2018, 47, 954–963. [Google Scholar] [CrossRef]
  18. Bai, X.; Shi, Y.; Liu, Y.; Liu, B. Reliability estimation of stress-strength model using finite mixture distributions under progressively interval censoring. J. Comput. Appl. Math. 2019, 348, 509–524. [Google Scholar] [CrossRef]
  19. Zhang, L.; Xu, A.; An, L.; Li, M. Bayesian inference of system reliability for multicomponent stress-strength model under Marshall-Olkin Weibull distribution. Systems 2022, 10, 196. [Google Scholar] [CrossRef]
  20. Tang, Y.; Zhang, J. New model for load-sharing k-out-of-n : G system with different components. J. Syst. Eng. Electron. 2008, 19, 748–751, 842. [Google Scholar] [CrossRef]
  21. Hellmich, M. Semi-Markov embeddable reliability structures and applications to load-sharing k-out-of-n system. Int. J. Reliab. Qual. Saf. Eng. 2013, 20, 1350007. [Google Scholar] [CrossRef]
  22. Bairamov, I.; Arnold, B.C. On the residual lifelengths of the remaining components in an nk + 1 out of n system. Stat. Probab. Lett. 2008, 78, 945–952. [Google Scholar] [CrossRef]
  23. Nguyen, D.P.; Kozyrev, D.V. Reliability Analysis of a Multirotor Flight Module of a High-altitude Telecommunications Platform Operating in a Random Environment. In Proceedings of the 2020 International Conference Engineering and Telecommunication (En&T), Dolgoprudny, Russia, 25–26 November 2020; pp. 1–5. [Google Scholar] [CrossRef]
  24. Rykov, V.; Ivanova, N.; Kochetkova, I. Reliability Analysis of a Load-Sharing k-out-of-n System Due to Its Components’ Failure. Mathematics 2022, 10, 2457. [Google Scholar] [CrossRef]
  25. Katzur, A.; Kamps, U. Order statistics with memory: A model with reliability applications. J. Appl. Probab. 2016, 53, 974–988. [Google Scholar] [CrossRef]
  26. Cramer, E.; Kamps, U. Sequential order statistics and k-out-of-n systems with sequentially adjusted failure rates. Ann. Inst. Stat. Math. 1996, 48, 535–549. [Google Scholar] [CrossRef]
  27. Navarro, J.; Marco, B. Coherent Systems Based on Sequential Order Statistics. Nav. Res. Logist. 2011, 58, 123–135. [Google Scholar] [CrossRef]
  28. Sutar, S.; Naik-Nimbalkar, U.V. A load share model for non-identical components of a k-out-of-m system. Appl. Math. Model. 2019, 72, 486–498. [Google Scholar] [CrossRef]
  29. Kozyrev, D.V.; Phuong, N.D.; Houankpo, N.G.K.; Sokolov, A. Reliability evaluation of a hexacopter-based flight module of a tethered unmanned high-altitude platform. Commun. Comput. Inf. Sci. 2019, 1141, 646–656. [Google Scholar] [CrossRef]
  30. David, H.A.; Nagaraja, H.N. Order Statistics, 3rd ed.; John Wiley & Sons: New York, NY, USA, 2003. [Google Scholar] [CrossRef]
Figure 1. An unmanned hexacopter flight module of a tethered high-altitude telecommunications platform.
Figure 1. An unmanned hexacopter flight module of a tethered high-altitude telecommunications platform.
Mathematics 10 04243 g001
Figure 2. An unmanned multi-functional underwater vehicle.
Figure 2. An unmanned multi-functional underwater vehicle.
Mathematics 10 04243 g002
Figure 3. Reliability function of a 2-out-of-6 system.
Figure 3. Reliability function of a 2-out-of-6 system.
Mathematics 10 04243 g003
Figure 4. Reliability function with v = 1 and quantiles ( A E x p ). Red and black bullets are the points of intersection of the reliability function curves with fixed reliability levels of 0.99 and 0.9, respectively.
Figure 4. Reliability function with v = 1 and quantiles ( A E x p ). Red and black bullets are the points of intersection of the reliability function curves with fixed reliability levels of 0.99 and 0.9, respectively.
Mathematics 10 04243 g004
Figure 5. Reliability function with v = 0.5227 and quantiles ( A G W ). Red and black bullets are the points of intersection of the reliability function curves with fixed reliability levels of 0.99 and 0.9, respectively.
Figure 5. Reliability function with v = 0.5227 and quantiles ( A G W ). Red and black bullets are the points of intersection of the reliability function curves with fixed reliability levels of 0.99 and 0.9, respectively.
Mathematics 10 04243 g005
Figure 6. Reliability function with v = 0.5 and quantiles ( A E r l ). Red and black bullets are the points of intersection of the reliability function curves with fixed reliability levels of 0.99 and 0.9, respectively.
Figure 6. Reliability function with v = 0.5 and quantiles ( A E r l ). Red and black bullets are the points of intersection of the reliability function curves with fixed reliability levels of 0.99 and 0.9, respectively.
Mathematics 10 04243 g006
Table 1. E [ Y 2 ] of a 2-out-of-6 system.
Table 1. E [ Y 2 ] of a 2-out-of-6 system.
c = 0.1 c = 0.5 c = 1
v = 0.5 ( A E r l ) 0.4925 0.5670 0.6668
v = 0.5227 ( A G W ) 0.4316 0.5251 0.6420
v = 1 ( A E x p ) 0.1867 0.2667 0.3667
Table 2. v s y s of a 2-out-of-6 system.
Table 2. v s y s of a 2-out-of-6 system.
c = 0.1 c = 0.5 c = 1
v c o m p = 0.5 ( A E r l ) 0.3802 0.3049 0.2986
v c o m p = 0.5227 ( A G W ) 0.4813 0.3854 0.3641
v c o m p = 1 ( A E x p ) 0.8993 0.7289 0.7100
Table 3. Quantiles of reliability function q γ .
Table 3. Quantiles of reliability function q γ .
c = 0.1 c = 0.5 c = 1
q 0.999 A E x p 0.0026 0.0059 0.0083
A G W 0.0419 0.0804 0.1027
A E r l 0.1019 0.1641 0.1945
q 0.99 A E x p 0.0088 0.0192 0.0271
A G W 0.0804 0.1458 0.1859
A E r l 0.1576 0.2345 0.2784
q 0.9 A E x p 0.0344 0.0691 0.0972
A G W 0.1813 0.2784 0.3517
A E r l 0.271 0.3574 0.424
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/math10224243

AMA Style

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(22):4243. https://doi.org/10.3390/math10224243

Chicago/Turabian Style

Rykov, Vladimir, Nika Ivanova, Dmitry Kozyrev, and Tatyana Milovanova. 2022. "On Reliability Function of a k-out-of-n System with Decreasing Residual Lifetime of Surviving Components after Their Failures" Mathematics 10, no. 22: 4243. https://doi.org/10.3390/math10224243

APA Style

Rykov, V., Ivanova, N., Kozyrev, D., & Milovanova, T. (2022). On Reliability Function of a k-out-of-n System with Decreasing Residual Lifetime of Surviving Components after Their Failures. Mathematics, 10(22), 4243. https://doi.org/10.3390/math10224243

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