ESTIMATION OF THE PROBABILITY OF ELIMINATING CRITICAL FAILURES OF DISTANTLY CONTROLLED
OBJECTS
Alexander Dokukin, Mikhail Lomakin
All-Russian Research Institute for Civil Defense and Emergency Situations of the Ministry of Emergency Situations of Russia, Moscow, Russia [email protected] [email protected]
Abstract
The problem of determining the probability of timely elimination of critical failures under conditions of incomplete data represented by small samples of hardware dwell time in failure state and critical hardware time is considered. This probability is the probability of stochastic dominance of the critical time of the hardware over the time of the hardware in the state of failure, i.e. the probability that the random time of the hardware in the state of failure is less than the critical time of the hardware. Critical time is the time during which the elimination of the occurred failures does not lead to the occurrence of an unrecoverable failure of the object or to a significant reduction in the efficiency of its functioning. The problem of determining the probability of timely elimination of critical failures is formulated as a problem of finding guaranteed (lower and upper) probability estimates on the set of available initial data, represented by small samples of the time of stay of the equipment in the state of failure and critical time of the equipment. The guaranteed estimates of the probability of timely elimination of critical failures are found using the results of solving the problem of Markov moments.
Keywords: probability, critical failures, distantly controlled objects
I. Introduction
Remotely Operated Controlled Objects (ROOs) are complex and expensive objects designed to solve important national economic, scientific and applied tasks. RCEs include artificial satellites, spacecraft, probes, pipelines and oil pipelines. In the process of their target functioning they are exposed to a large number of external factors. [1] This impact in most cases is the cause of failures of their hardware. Hardware failures lead to different consequences for DKO. In the case of some hardware failures, the recovery of its performance is automatic with the built-in system of control, diagnostics and recovery, in the case of other failures for their detection, diagnosis and recovery requires the involvement of ground control complex facilities, including ground control and management facilities [1]. Among the failures of the second group there are such failures, the untimely elimination of which leads (may lead) to the occurrence of unrecoverable hardware failure, and further to the failure of RSC or to a significant reduction in the efficiency of its operation. Let us call such hardware failures critical failures; they include failures of the power
153
supply system, failures of the satellite navigation system, failures of the orientation and stabilization system, etc.
Let at a random moment of time there is a hardware failure, which after some random time % manifests itself, as a rule, is detected during the telemetry control session following the failure. The detected failure is diagnosed and eliminated within the period of time y.
Currently, only the task of technical condition monitoring is solved on board some ROVs. The task of diagnostics and technical condition management in most cases is realized remotely. [2].
The random dwell time of the hardware in a failure state a will consist of the time of failure occurrence Z, and the time of diagnostics and restoration of the equipment operability y. The value of time a will be determined by the ratio:
a = { + y, (1)
where a - is the random dwell time of the equipment in the state of failure, % - random time of hardware failure manifestation, y - random time of diagnostics and restoration of equipment operability. Let t] - a random time during which the elimination of critical failures does not lead to an unrecoverable hardware failure or to a significant reduction in the efficiency of its operation. Let us call this time the critical time.
The task is to determine the probability of timely elimination of DKO equipment failures - the probability that the random dwell time of the equipment in the state of failure will be less than the critical time, i.e. to determine the probability ofP = P < ]), the probability of stochastic dominance of the critical time of the hardware over the time of the hardware stay in the state of failure.
At known functions of distribution of time of stay of the equipment in a state of failure and critical time this problem is solved. Ratios for determining the probability P = P < ]) with respect to the basic distribution functions of random variables a, t] are given, for example, in [3]. In a real situation, the distribution functions of the time of the equipment stay in the state of failure and the critical time are not known. The present paper is devoted to the consideration of this case.
II. Main results
Let av = (ov1,ov2,... ,avnfT)Rn(T - a sample of values of the time of the equipment staying in the failure state, tjv = (rjv,, >jv2,..., rjvnY)Rnrt - sample of values of critical time of the equipment® The components of the sample avi > 0 are independent equally distributed values from some unknown distribution H(t), ]vi > 0 are independent equally distributed values from some unknown distribution G(t). The samples av, ]v are finite samples of small volume, for which it is impossible to recover the original distributions F(t) and G(t).
Let us transform the expression for the probability of timely elimination of critical failures. We have
P = P(a < t]) = P (a < l) = P(aA <1)=P(J3 < 1). (2)
whereA = 1;3 = aA. v
Let's find on the basis of samples av, ]v sample initial moments (hereinafter - moments) of random variables a, A, ft:
na nv .
= av} ;^A=—Y(—) ;j = 1k;k>0. (3)
i=l i=l
Random variables a, A are independent quantities, therefore, the random variables a},A} are also independent. The mathematical expectation of the product of random variables a} ,A} is equal to the product of mathematical expectations of these quantities
M(a}A})=M)M(A};j = 1k;k>0. (4)
Then the moments of the random variable /3 will be defined by the expression:
Vjii = VjeVjX ;j = l,k;k>0. (5)
Let us define the set F0 - the set of distribution functions F(t) of the random variable ft with moments equal to the moments defined by the relation (8)
Fo = [F(t) ■ f™ t'dF(t) = Hrf, j= lj;k> 0}. (6)
The problem of determining the probability of timely elimination of critical failures of DKO equipment due to incompleteness of initial data represented by small samples of the time of equipment stay in the state of failure and critical time of the equipment should be formulated as a problem of determining the guaranteed (lower and upper) estimates. It is necessary to determine the lower and upper (guaranteed) estimates of probability P = P(ft < 1) on the set of distributions F0 , i.e. to find
= < 1):= < 1). (7)
Ratio P = P(p <1) can be written in the form:
P = P(ft<1)=f dF(t). (8)
■>0
The solution of the latter problem, defined by relation (10), is obtained in [5-10] and it is summarized as follows.
Highest (lowest) value of the integral
J(F) = f c(t)dF(t). (9)
■>0
at F(t)-rFo is achieved on the only stepwise distribution F(t)e which among the growth points ti-r try ...,t v has a point-»-;
at odd k the number of growth points v of the distribution function F(t) is determined by the relation v = ( k + 3) 2tt# and to = 0 < t < tu < ...< t < v •«- •«at even k the number of growth points v of the distribution function F(t) is determined by the relation v = ky 2 + ly with 0 <ti<ti < ...< t<u yy
numbers p ,y tjyj = 1 2y*y*y*y*y*y v satisfy the system of equations:
Vip=^t}iPi;j = 1,k;k>0. (10)
i=i
The function c(t) must have a non-negative k+1th derivative.
Let us consider for certainty the problem of determining the lower estimate of probability P = P(P <1) on the set of distribution functions F0 .
Let two moments of a random variable be defined ft - (further, where it does not cause
discrepancies, the index ft is omitted) then using the above result we find the lower estimate of probability P = P(ft < 1)which is defined by the following expression [6]:
(P1 - 1)2
^^-rf.^-^»- (11)
In the case when three moments of a random variable are defined ft - we find the
lower estimate of the probability P = P(ft < 1)which is determined by the following expression [6]:
Vf13 + Vs1 - 2^.12 - V3V1 + vl „ ^ ^-^+,31-2,212-,1>"1- (12)
In the last ratio, the symbol "1" (t=1 - growth point of the distribution function F(t)) is used to confirm that the numerator and denominator dimensions are the same.
In the case when more than three moments of a random variable are defined ft the estimate of the probability P = P(ft < 1) should be calculated numerically.
Let us consider the case when four moments of the random variable are defined / -In this case, the number of growth points of the distribution function F(t) is equal to three, the system of equations (10) will have the form:
Pl+P2+P3 = 1, Pltl + P2t2 + Psts =Vl. Plt^+P2t^+Pst32 = V2,
Pltl+P2tl + PstSS =
(13)
Vs,
Pltl+P2t2 + Pst3 = V4.
Since the distribution function F(t) is a left continuous function, the growth point t=1 can be either t2 or tsthen we should consider two systems of equations
Pl + P2 + Ps = 1, Pltl + P2T + Psts =Vl,
Pltl+P2T2+Pstl = V2, Pltl+P2T3 +Pst3 = Vs, Pltt+P2ji + Pst3 = V4,
(14)
Pl + P2 + Ps = 1,
Pltl + P2t2 + PS* = Vl, Pltl+P2t2+PsT2 = V2, Pltl+P2tl + PS?3 = Vs, ■Pltt+P2t2+Ps*4 = V4.
(15)
To solve the systems of equations (14), (15) we used Microsoft Excel package (item) "Solution Search" sub-item "Search for the solution of nonlinear problems by OPG method1 ". In this case, the first equation of the systems of equations (17), (18) was considered as a target function of the following form
CF = Pl + P2+Ps- 1
the value of which should be equal to zero, the other equations were considered as constraints in the optimization problem.
The results of numerical solution of the problem of determining the lower estimate of the probability of P = P(J3 < 1) - of the probability of timely elimination of failures of DKO equipment on the set of distribution functions Fo at different values of the first four moments are given in Table 1.
Table 1: Lower estimate of the probability of timely elimination of failures space craft hardware
V2 Vs V4 P 1 n
0,200 0,080 0,048 0,038 0,959
0,300 0,180 0,162 0,194 0,832
0,400 0,320 0,384 0,614 0,648
0,500 0,500 0,750 1,500 0,579
0,600 0,720 1,296 3,110 0,488
0,700 0,980 2,058 5,762 0,429
0,800 1,280 3,072 9,830 0,330
1 The OPG method is the generalized reduced gradient method.
156
III. Conclusion
The problem of determining the probability of timely elimination of critical failures of DKO equipment in the conditions of incomplete data represented by small samples of the time of equipment stay in the state of failure and critical time of the equipment is considered. This problem is formulated as the problem of determining guaranteed (lower and upper) estimates of the probability that the random dwell time of the equipment in the failure state is less than the random critical time of the equipment, i.e., the probability of stochastic dominance of the critical time of the equipment over the dwell time of the equipment in the failure state. To determine the probability of timely elimination of critical failures of DKO equipment, the results of solving the problem of Markov moments are used, with their help we obtained relations for determining the lower estimate of this probability at two and three moments of a random variable equal to the ratio of the time of stay of the equipment in the state of failure to the critical time of the equipment. At four moments of the last random variable, the lower estimates of the probability of timely elimination of critical failures of DKO equipment were determined numerically using Microsoft Excel.
References
[1] Novikov L.S. Radiation Effects on Spacecraft Materials. - M.: Universitetskaya kniga, 2010. -192 c.
[2] Baranovsky, A.M.; Privalov, A.E. System of control and diagnostics of the onboard equipment of a small spacecraft (in Russian) // Izv. VUZov. Instrument engineering. 2009. T.52, № 4. C.51-56.
[3] Ostreikovsky V.A. Multifactor reliability tests. - M.: Energia, 1978. -152 c.
[4] Gnedenko B.V. Probability Theory Course. - Moscow: Fizmatgiz, 1988. - 406 c.
[5] Lomakin M® And® Guaranteed estimates of probability of failure-free operation in the class of distributions with fixed moments®® Izvestia AS USSR. Automatics and Telemechanics. 1991. NO. 1® FROM® 154-161 ®
[6] Lomakin M., Buryi A., Dokukin A., Niyazova J. Strekha A., Balvanovich A.. Estimation of quality indicators based on sequential measurements analysis // International Journal for Quality Research. 2020. No. 1. pp. 823-834.
[7] Lomakin, M.I.; Sukhov, A.V.; Dokukin, A.V.; Niyazova, Yu.M. Estimation of the spacecraft reliability indicators in the conditions of incomplete data // Space Research. 2021. T. 59. № 3. C. 235-239.
[8] Lomakin M.I., Niyazova Y.M., Dokukin A.V., Zlydnev M.I., Garin A.V. Quality assessment of business processes in the conditions of incomplete data // Welding production. 2022. № 4. C. 52 - 58.
[9] Crane M.G., Nudelman A.A. The problem of Markov moments and extremal problems (Ideas and problems of P.L. Chebyshev and A.A. Markov and their further development). -Moscow: Nauka, 1973 . - 551 c.
[10] Akhiezer N.I. Classical problem of moments. - Moscow: Fizmatgiz, 1961. - 310 c.