Next Article in Journal
Lyapunov Inequalities for Two Dimensional Fractional Boundary-Value Problems with Mixed Fractional Derivatives
Next Article in Special Issue
A Study of the Monotonic Properties of Solutions of Neutral Differential Equations and Their Applications
Previous Article in Journal
On Some Branched Continued Fraction Expansions for Horn’s Hypergeometric Function H4(a,b;c,d;z1,z2) Ratios
Previous Article in Special Issue
Further Closed Formulae of Exotic 3F2-Series
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Models of the Processes of Operation, Renewal and Degradation of a Fleet of Complex Technical Systems with Metrological Support

1
Department of “Fundamental Education”, Mytishchy Branch of University, Moscow State University of Civil Engineering, 26, Yaroslavskoye Shosse, 129337 Moscow, Russia
2
Metrology Scientific Head Center, 13 Komarova Street, Moscow Region, 141000 Mytischy, Russia
3
Department of “Standartizaion, Engineering Faculty, Moscow Polytechnic University, Metrology and Sertification”, 38, Bolshaya Semenovskaya Street, 107023 Moscow, Russia
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(3), 300; https://doi.org/10.3390/axioms12030300
Submission received: 25 January 2023 / Revised: 6 March 2023 / Accepted: 8 March 2023 / Published: 15 March 2023
(This article belongs to the Special Issue 10th Anniversary of Axioms: Mathematical Analysis)

Abstract

:
(1) Background: The aim of the study is to develop a set of models for managing a fleet of complex technical systems with metrological support, allowing the simulation and management at all the stages of the life cycle of the complex technical systems, as well as to simulate the functioning of large fleets of complex technical systems, including up to several hundred thousand samples; (2) Methods: The authors use methods of mathematical modeling, methods of the theory of Markov and semi-Markov processes, methods of optimization, methods of reliability theory, and methods of probability theory and mathematical statistics; (3) Results: an interconnected set of mathematical models for managing a fleet of complex technical systems with metrological support was developed and the applied software was developed; (4) Conclusions: The set of models presented in the article allows for the adequate simulation of all the stages of the life cycle of large complex technical systems fleets, including up to several hundreds of thousands of samples, to optimize the functioning processes of a fleet of complex technical systems, to form strategies for fleet development, and to assess the risks associated with false and undetected failures, as well as the risks associated with the degradation of complex technical systems.

1. Introduction

A considerable amount of scientific research is devoted to the problem of modeling complex technical systems (CTS) [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28]. We understand complex technical system as stationary or mobile special-purpose objects with measuring instruments (MI) installed on them, which should be metrologically maintained during long-term operation. In the last half of the century, both CTS themselves and their models have undergone a rather rapid evolution process. Starting from models with 3–5 states and going up to models with up to several hundreds and thousands of states. At the same time, the theoretical base and technical capabilities for modeling CTS with several tens of thousands and even hundreds of thousands of states have been created.
On the qualitative side, simple models allowed modeling only of the basic states of the CTS, which describe the operation processes. The models have been evolving toward a more detailed description of the operation processes (taking into account metrological support technologies, false failure states and undetected failure states), the CTS degradation processes (first degradation level, second degradation level and so on) and the CTS updating processes (by purchasing the new CTS samples, by upgrading the existing CTS samples, by developing the newest modern CTS samples). Thus, by now there is a need for models that describe all stages of the CTS life cycle. These models should allow the simulation of large CTS fleets, including up to several hundreds of thousands of STS samples. The models should make it possible to manage the process of development of such CTS fleets, taking into account the range of modern tasks to be solved by means of the CTS and the need to solve promising tasks in the future.
Let us first conduct a retrospective comparative analysis of CTS models, with a separate description of the main characteristics of each model, as well as the assumptions underlying their implementation. Let us describe the strengths and weaknesses of the models. Additionally, we will then formulate the goal of our scientific research and we will provide a statement on the problem that will be investigated in this article.

2. Scientific Literature Review

Professor L.I. Volkov [1] proposed the semi-Markov model of aircraft operation control, which has five states: workable status; periodic verifications of the operational status; recovery after the occurrence of the valid state, false failure state; the hidden failure state; the unworkable state (including the hidden failure state); and the state of periodic verifications with hidden failure.
The classical model developed by Professor E.I. Sychev [2], designed to control the process of operation of the CTS with measuring instruments (MI) installed on them to provide metrological support, in contrast to the model by Professor L.I. Volkov, already has six states. Model [2] describes the operation process more correctly. From the fourth unworkable state (including hidden failure), two states were separately highlighted: the state of undetected failure and detected failure. The model takes into account the characteristic features of the CTS with metrological support.
The model [2] assumes the identity of the recovery of the CTS after both a false failure and a valid failure. In practice, for some types of CTS, after a false failure, repeated control is carried out according to the failed technical parameter, and after a detected failure, the system is restored, for example, by adjusting or replacing the faulty element with a serviceable one. In the model [3] developed by Professor V.I. Mishchenko, which already includes seven states, the above-mentioned features and limitations have been eliminated. The model [3] takes into account the intensity of the CTS operation.
Note that the models described above do not take into account the component of maintenance efficiency, determined by the availability of spare parts and their replenishment strategy.
The further direction for the development of the models for the operation of the CTS is to take into account the possibilities of reserving the MI and the possibilities of replenishment with spare parts and tools. In [4], the model of the process for the functioning of the MI with metrological support for doubly redundant MI is proposed, which allows for the taking into account of the features of the maintenance associated with the possibility of providing spare parts, and taking into account the different strategies for replenishing spare parts, tools and accessories. In the model [4], which takes into account eight states, it is assumed: that the detection of failures by the MI occurs only during verification; there are no errors in determining the technical condition of the MI; and the MI in storage do not fail.
In [5], a new approach has been developed to assess the impact of metrological support on achieving the goals of the CTS operation: a graph with an arbitrary number of states is constructed, the edges of the graph that represent possible state transitions are attributed both probabilistic characteristics of the transitions (values of the distribution functions or simply the transition probabilities) and the costs associated with the corresponding transitions. The following states are selected: serviceable, faulty, emergency and catastrophic. The results from the study, on the influence of the volumes of metrological control for various conditions on the effectiveness of the object for its intended purpose, are presented. As a criterion of efficiency in the various solved tasks, both the readiness coefficient and the technical and economic indicator were used.
It should be noted that all the models analyzed above do not allow modeling and the taking into account of conditions corresponding to the different levels of degradation of the CTS (different levels of deterioration of the metrological reliability characteristics), leading to time and resource costs necessary both for restoring the CTS and bringing it back into working condition. Further development of the CTS operation models takes place in terms of taking into account the aging and degradation processes [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26] of the CTS (or MI installed on them) and reduction of the metrological reliability.
Thus, in [6] the model with four degradation groups is considered, having one workable state and four states corresponding to the different levels of degradation. This model describes the process of operation of the CTS, for which repair is possible with the restoration of the resource in full. In [7], a model with three degradation groups is considered, which allows for the modeling of the processes of operation, renewal and degradation of the CTS fleet. It is assumed that as a result of the repair, the resource of the CTS cannot be fully restored.
The works analyzed in this section form the basis (starting point) for the research presented in the article. This article summarizes the results of the work [5,6,7,8]: a set of models describing the processes of operation, renewal and degradation of the CTS are presented. To describe the operation process, the classical model [2] is used as it is the most adequate for the CTS class considered in the article. To describe the processes of degradation and renewal of the CTS fleet, new additions to the classical model developed by the author are presented (the model of false and undetected failures, the model of degradation and renewal of the fleet, including CTS with full and incomplete restoration of the resource during repair and metrological maintenance).

3. Statement on the Research Problem

It is necessary to develop a set of interrelated mathematical models of CTS fleet management models, allowing for the simulation and management of all stages of the CTS life cycle. The developed set of interrelated models should allow for the simulation of the functioning of large fleets of CTS, including up to several hundreds of thousands of CTS samples. The set of models shall allow for the taking into account of the degradation processes of CTS sample ageing, processes on park development due to the procurement of new samples, the modernization of existing samples and the development of new promising CTS samples. The set of interrelated models should allow for the management of the process of development of such CTS fleets, taking into account a number of modern requirements, and the need to solve promising tasks and problems in the future.

4. Materials and Methods

At first in Section 4.1.1, the results of calculating the readiness coefficient for different failure distribution laws using the classical operation model are presented. The model of false and undetected failures is described in Section 4.1.2. Section 4.1.3 describes and analyzes the models of failure and degradation of the CTS (a fan model, a drift model of the metrological characteristics and two diffusion models). In Section 4.2, the model of operation of the CTS is described, taking into account the degradation processes and the full restoration of the resource, and in Section 4.3, the model of the CTS with incomplete restoration of the resource is described.

4.1. The Classical Model

4.1.1. Construction and Study of the Classical Model for Different Laws on the Distribution of Failures of the Complex Technical System

Let us denote E i , i = 1 , 2 , n as a finite set of states in which a specific sample of the CTS can be located. The readiness coefficient of the CTS, the operation process of which is described by the semi-Markov model [2], is calculated by the formula:
K A = i = 1 n π i w i / i = 1 n π i ψ i ,
where π i is the relative fraction of the number of steps during which the CTS is in state E i , w i is the mathematical expectation of the time of operation of the CTS in state E i , and ψ i is the mathematical expectation of the time that the CTS stays in state E i .
At the same time: i = 1 n π i = 1 , ψ i = i = 1 n P i j M ( τ i j ) = i = 1 n P i j 0 τ i j d F ( τ i j ) ,
w i = ψ i for   workable   conditions   of   CTS 0 for   unworkable   conditions   of   CTS ,
where P i j are the elements of the state transition probability matrix P * = P i j * , F * ( τ i j ) is the transition probability distribution function, and M ( τ i j ) is the mathematical expectation of the transition time.
A continuously operating CTS with periodic verification of the technical condition is ready for use at that time τ if it is operational at that moment and is not under verification or repair. The results of the control are used to make a decision on the possibility of further application of the CTS. If the CTS is recognized as workable, according to the results of the verification, then it is included in the work. If the CTS is found to have failed, then its repair is carried out, as a result of which a complete restoration of its operability occurs. The transition graph is shown in Figure 1.
Possible conditions of the CTS: E 1 is workable, E 2 is unworkable (failure), E 3 is verification of the failed CTS, E 4 is recovery, E 5 is verification of a workable CTS, and E 6 is undetected failure.
The transition probability matrix has the following form:
P * = 0 F ( T K ) 0 0 1 F ( T K ) 0 0 0 1 0 0 0 0 0 0 1 β 0 β 1 0 0 0 0 0 1 α 0 0 α 0 0 0 0 1 0 0 0 ,
where F ( τ ) is the integral function of the distribution of the failure time, F ( T K ) is the probability of failure during the time between two verifications, T K is the time interval between verifications (TIBV) of the technical condition, α is the conditional probability of a false failure, and β is the conditional probability of an undetected failure.
We assume that the duration of the control (verification of the technical condition) and the duration of the restoration (repair) are deterministic values equal to t K and t B , respectively.
The system of equations for finding π i , i = 1 , 2 , , 6 has the form:
π 1 = π 4 + ( 1 α ) π 5 ,   π 2 = F ( T K ) π 1 ,   π 3 = π 2 + π 6 ,   π 4 = ( 1 β ) π 3 + α π 5 , π 5 = 1 F ( T K ) π 1 ,   π 6 = β π 3 ,   π 1 + π 2 + π 3 + π 4 + π 5 + π 6 = 1 .
The solution of the system has the form:
π 1 = 1 A ( 1 β ) π 2 = 1 A F ( T K ) ( 1 β ) π 3 = 1 A F ( T K ) π 4 = 1 A F ( T K ) + α ( 1 F ( T K ) ( 1 β ) π 5 = 1 A 1 F ( T K ) ( 1 β ) π 6 = 1 A β F ( T K ) ,
where A = 2 1 β + F ( T K ) + α 1 F ( T K ) ( 1 β ) .
The values v i , i = 1 , 2 , , 6 are equal to:
v 1 = 0 T K τ d F ( τ ) + T K 1 F ( T K ) v 2 = T K F ( T K ) 0 T K τ d F ( τ ) v 3 = t K v 4 = t B v 5 = t K v 6 = T K .
Assuming that w 1 = v 1 , w 2 = 0 , w 3 = 0 , w 4 = 0 , w 5 = 0 , w 6 = 0 , and substituting (2) and (3) into (1), we obtain the formula for calculating the CTS readiness coefficient:
K A = I + T K B B I + T K B + F ( T K ) 2 + β F ( T K ) 1 β + t k B + F ( T K ) 1 β + t B F ( T K ) + α B ,
where B = 1 F ( T K ) , I = 0 T K τ d F ( τ ) .
Next, we will conduct a study of the readiness coefficient for various laws on the distribution of the failure time. The failure time of the CTS is considered as a random variable. Analysis of statistical data has shown that the most suitable laws for describing the failure time are the exponential law, Rayleigh’s law, the Weibull distribution and the truncated normal distribution, with the appropriate choice of parameters for these distributions. The statistical function of the distribution of the failures is located inside the “curved band” covering the theoretical distribution functions.
In the case of the exponential distribution law, the expression for the readiness coefficient (4) takes the form:
K Γ = 1 e λ T K λ T K 1 β + e λ T K 1 e λ T K + λ t k ( 1 e T K ) β p 1 β + 1 + λ t B 1 e λ T K ( 1 α p ) .
For Rayleigh’s law, the integral I can be calculated numerically or using the standard Laplace function, for Weibull’s law it can be calculated numerically or using the gamma function; and for the truncated normal distribution it can be calculated numerically or using the standard Laplace function.
The calculations were carried out using the following values from the initial data: t K = 1 , t B = 1 , α = 0.1 , β = 0.1 , and λ = 0.0025 for the different values of T K . Figure 2 shows the dependences of the readiness coefficients K A on the periodicity of the verification T K , for the distribution laws described above. The maximum values of K A for the Rayleigh, normal, exponential and Weibull laws are equal to: 0.976, 0.963, 0.955, and 0.950, respectively, and reach values equal to 65, 55, 50, and 40. Note that the maximum value of the coefficient for each distribution law is reached at a single point. It can be seen that max K A is “practically insensitive” to T K . So, in a fairly wide range of changes to T K , the readiness coefficient takes values close to the maximum. In particular, when changes to T K take place in the range 25 T K 60 , the variation of K A is no more than 2–3%.
The low sensitivity of the maximum value of the readiness coefficient to the periodicity of the technical condition monitoring makes it possible to develop strategies that are “non-strict” and easy to implement in practice, for carrying out checks on the technical condition of the CTS with metrological support.

4.1.2. Development of the Classical Model: The Model of False and Undetected Failures

The probabilities of false and undetected failures [8] for the specific samples of the CTS depend on the corresponding probabilities of false and undetected failures of the individual components of the CTS (1), (2), on the configuration of the CTS using methods on the redundancy of the components, nodes and blocks of the CTS.
Let p be the actual value of the measured (controlled) parameter and ε be the measurement error. The measurement result is presented in the form r = p + ε . The general scheme of the diagnosis and decision-making based on the one-parameter method of tolerance control is shown in Figure 3.
Here δ is the tolerance for the controlled parameter, and f ( x ) and f e ( x ) are the distribution density functions of the measured parameter and the measurement error, respectively. It can be seen that the probability of making the right decision can be increased (within certain limits) by reducing the total error of the erroneous decision.
The different physical nature and, consequently, the heterogeneous range of the changes in the measured values leads to the need to introduce dimensionless standardized operational parameters for the MI. As a normalizing element, we take the mean square deviation σ x of the measured parameter x ; δ = Δ / σ x is the relative operational tolerance, where Δ is the technical tolerance; z = σ ε / σ x is the relative parametric measurement error, σ ε is the mean square deviation of the MI error.
The model is based on formulas for the conditional probabilities of false and undetected failures, respectively [8]:
α ( δ , z ) = δ δ f c u ( y ) δ y z f o ( τ ) d τ + δ y z f o ( τ ) d τ d y / δ δ f c u ( y ) d y ,
β ( δ , z ) = δ f c u ( y ) δ y z δ y z f 0 ( τ ) d τ d y + δ f c u ( y ) δ y z δ y z f 0 ( τ ) d τ d y / δ f c u ( y ) d y + δ f c u ( y ) d y ,
where f 0 ( τ ) and f c u ( y ) are the functions of the distribution densities of the measured value and the MI error, respectively.
For normally distributed measured values and MI errors, Formulas (5) and (6) take the form:
α ( δ , z ) = 1 2 π δ δ exp y 2 2 δ y z exp τ 2 2 d τ + δ y z exp τ 2 2 d τ d y / P 1 ,
β ( δ , z ) = 1 2 π δ exp y 2 2 δ y z δ y z exp τ 2 2 d τ d y + δ exp y 2 2 δ y z δ y z exp τ 2 2 d τ d y / P 2 ,
P 1 = 1 2 π δ δ exp y 2 2 d y ,   P 2 = 1 2 π δ exp y 2 2 d y + δ exp y 2 2 d y
For other distribution laws on the measured value and measurement error, the model (5), (6) were investigated in [8].
The dependences of the probabilities α ( δ , z ) and β ( δ , z ) , as well as the probability α ( δ , z ) + β ( δ , z ) of an erroneous decision δ on the tolerance value at z = 0.5 are shown in Figure 4.
Note also that the error solution function reaches its minimum at some internal point δ ( 0 ; 1 ) , as is the case with the normal distribution of the measured value and the measurement error.
The two-dimensional dependences of the probabilities of false and undetected failures on the magnitude of the dimensionless measurement error z and the dimensionless tolerance for the controlled parameter δ are shown in Figure 5a,b.

4.1.3. Development of the Classical Model: Models of Failures and Degradation of the Complex Technical System

All failure models that allow for the taking into account of the degradation processes occurring in the CTS can be conditionally divided into probabilistic, empirical, and probabilistic physical models, that includes among other things, the Markov models of degradation and failures.
In the fan model [9,10,11], also called a distribution and belonging to the category of probabilistic models, the defining parameter (DP) is represented as a linear function of time, shown in Figure 6a.
Here, t is the operating time for the failure; X is random variable of the DP; DP* is the normalized value of the DP at which the failure occurs; and f is the function of density of the distribution of the operating time for the failure.
The distribution function of the operating time up to a given level DP* is given by the distribution function [9,10,11]:
F ( t ) = Φ t μ v t ,
where Φ is the normalized normal distribution function; μ = 1 / a is the parameter of the scale of degradation, a is the mathematical expectation of the rate of change of the DP (the average rate of the degradation process), normalized to the limit value; and ν is the shape parameter (coefficient of the variation of the degradation process).
The empirical model of the “drift of metrological characteristics” [6,7] is based on the assumption on a linear law of change of the MI zero mark and an exponential law of increasing measurement error:
m 0 ( t ) = m 00 + v m t ,   σ ( t ) = σ 0 + v z a z exp ( a z t ) 1 ,
where σ 0 is the value of the initial error, v z is the average initial velocity of the error increase, a z is the parameter characterizing the acceleration of the error increase, m 00 is the initial value of the zero drift (usually assumed to be zero), and v m is the average velocity of the zero drift.
Let us now consider the Markov models of degradation and failures, widely used in applied problems. In these models, it is assumed that the degradation process can be approximated by a continuous Markov process of the diffusion type [9,10,11] and is described by a stochastic differential equation of the Ito type:
d x ( t ) = A ( t ) d t + B ( t ) d η ( t ) ,
where x ( t ) is the value of the DP; A ( t ) and B ( t ) are deterministic functions characterizing the change in the mean value and variance of the DP (drift coefficient and diffusion coefficient); and η ( t ) is a random variable of the Gaussian type.
The problem of determining the distribution of time before the first failure of the MI, in this case, is reduced to solving the problem of the first achievement of the upper limit of the DP* (see Figure 6b,c). This problem can be solved if the conditional probability density ω ( t , x ) of the process transition from one state to another is known.
For a Markov diffusion-type processes, a partial differential equation (the Fokker–Planck–Kolmogorov equation) follows from (11):
ω ( t , x ) t + A ( t ) ω ( t , x ) x B ( t ) 2 2 2 ω ( t , x ) x 2 = 0 ,
where A ( t ) and B ( t ) are the coefficients of the equation depending on the operating conditions of the MI, and the physical and chemical processes occurring in the materials from which the MI is made. To solve (12), it is necessary to set boundary conditions that depend on the type of implementation of a random process, in particular, on their monotonic nature (Figure 6b) or non-monotonic nature (Figure 6c). You also need to set the initial conditions: t = t 0 , x = x 0 .
After finding the function ω ( t 0 , x 0 ; t , x ) , satisfying the given initial conditions, the density function f ( t ) of the distribution of the time to reach the boundary DP* (the density function of the distribution of the time to failure) can be calculated by the formula [11]:
f ( t ) = t ω ( t 0 , x 0 ; t , x ) t d x
In case of one DP, Equation (12) can be integrated analytically. The distribution function for the diffusion monotone distribution (DM distribution) has the form [11]:
F ( t ) = D M ( t ; μ , ν ) = Φ t μ v μ t ,
Here, μ = 1 / a .
The distribution density function f ( t ) for (13) at μ = 0.1 is shown in Figure 7a.
The distribution function for the diffusion non-monotonic distribution (DN distribution) has the form [11]:
F ( t ) = D N ( t ; μ , ν ) = Φ t μ v μ t + exp 2 ν 2 Φ t μ v μ t .
The corresponding distribution density functions f ( t ) for (14) at ν = 0.8 are shown in Figure 7b.
The failure rates for DM distribution and DN distribution have the form:
λ D M ( t ) = ( t + μ ) exp ( t μ ) 2 2 ν 2 μ t 2 ν t 2 π μ t Φ μ t v μ t ,   λ D N ( t ) = ν t 2 π t 1 μ exp ( t μ ) 2 2 ν 2 μ t Φ μ t v μ t exp 2 ν 2 Φ μ t v μ t
Thus, distribution density functions f ( t ) , distribution functions F ( t ) , and failure rate functions λ ( t ) , are calculated using finite analytical formulas using the standard Laplace function Φ ( t ) .
In case of several DP distribution densities f ( t ) , distribution functions F ( t ) , and failure rates λ ( t ) , can only be calculated numerically.
The process of degradation of the mechanical components of the CTS, due to the irreversibility of the destruction processes (mechanical wear, fatigue straining, etc.), is considered to be a process with monotonous realizations of a random variable. DM distribution is used for CTS nodes containing electromechanical elements (relay and connector contacts, sliding electrical contacts, gears, etc.) [11].
The process of degradation of the CTS, which include integrated circuits and complex electronic devices, also has non-monotonic implementations of a random variable. Therefore, the degradation of such CTS is described by the DN distribution [11].
We will analyze the models of failures and degradation of the CTS. Degradation and failures models differ significantly from a physical point of view. In particular, the fan process assumes that its characteristics are completely determined by the initial state (the quality of the manufacturing samples of the components of the CTS), and do not depend on the mechanical, physical and chemical degradation processes occurring in the circuits and mechanisms of the components of the CTS, under the influence of external conditions and time.
The drift model of metrological characteristics [10], clearly demonstrates the departure of the zero mark of the MI and CTS with the increase in measurement error over time. The model assumes preliminary processing of statistical data in order to determine estimates of the drift parameters.
The Markov models (12), (13) are based on the use of probabilistic characteristics, the operating conditions of the CTS, as well as on the use of the physical and chemical properties of the materials. The advantage of Markov models [12,13] is that they have accurate analytical expressions for all statistical characteristics, including statistical moments. In addition, there are no analytical expressions for the statistical moments of the fan α distribution law. These moments are determined by approximate dependencies, which complicates the use of a fan distribution in practice.
The density distribution function of the DM distribution occupies an intermediate position between the, widely used in practice, normal distribution (which is symmetrical) and the more elongated distribution.
The density curves of the DN distribution have a more significant insensitivity threshold, a more positive kurtosis and are more asymmetric than the DM distribution.
The intensities of the diffusion distributions have finite limits:
lim t λ D M ( t ) = lim t λ D N ( t ) = 1 2 μ ν 2 .
Here, are some important properties of diffusion distributions for practical application:
1. Where a random variable T is described by a DM distribution of the form D M ( t ; μ , ν ) , then the random variable x = c T ( c = c o n s t ) is also described by a DM distribution of the form D M ( t ; c μ , ν ) .
2. Where a random variable T is described by a DM distribution of the form D M ( t ; μ , ν ) , then the random variable θ = 1 T is also described by a DM distribution of the form D M t ; 1 μ , ν .
3. Where a random variable T is described by a DN distribution of the form D N ( t ; μ , ν ) , then the random variable x = c T ( c = c o n s t ) is also described by a DN distribution of the form D N ( t ; c μ , ν ) .
4. The sum of n random variables obeying the distribution of the form D N ( t ; μ , ν i ) is described by the DN distribution of the form D N t ; n μ , 1 / i = 1 n ν i 2 .
5. The sum of n random variables obeying a distribution of the form D N ( t ; μ i , ν ) is described by a DN distribution of the form D N t ; i = 1 n μ i , ν / n .
6. The sum of n random variables obeying the DN distribution of the form D N ( t ; μ , ν ) is described by the DN distribution of the form D N t ; n μ , ν n .
The proof of properties 1–6 can be carried out by replacing the variables and definitions of functions (13)–(14).
Some additional properties of diffusion distributions are described in [11].
Analysis of the graphs on the distribution functions shows that distributions (9), (12), (13) have different zones of high reliability. This means that the estimation of small-level quantiles, i.e., the assignment of a gamma-percent resource, significantly depends on the selected type of failure model of the CTS.
Diffusion models can be parameterized quite simply in the presence of statistical information. For example, when parameterizing based on statistical data on the moments of failure t i , ( i = 1 , 2 , , N ) , the estimates of the parameters μ and ν calculated using the maximum likelihood method for the DM distribution have the form:
μ = 1 N i = 1 N t i ,   ν = μ 1 N i = 1 N 1 t i N i = 1 N t i 1 ,
and for DN distributions have the form:
μ = N i = 1 N 1 t i 1 + N 2 i = 1 N ( μ + t i ) 1 1 N 2 4 i = 1 N ( μ + t i ) 1 2 N i = 1 N 1 t i 1 1 N i = 1 N t i N i = 1 N 1 t i 1 ,
ν = 1 μ N i = 1 N t i + μ N i = 1 N 1 t i 2
Thus, diffusion models are more preferable (adequate), since, unlike the fan distribution and the drift model of metrological characteristics, they can be used to control the degradation and reliability of the CTS, based on the taking into account of the physical patterns implemented through time-dependent variable coefficients A ( t ) and B ( t ) in Equation (12). The task of developing models of physical processes for the purpose of constructing coefficients A ( t ) and B ( t ) is an independent scientific task and is not considered in this article.

4.2. The Model of Operation of a Complex Technical System Fleet with a Fully Recoverable Resource

The attribution of a set of CTS to one or another degradation group is carried out on the basis of structural and functional analysis of the metrological reliability indicators [2,6,7], which includes the types of failures, the consequences of the failures, as well as determination and analysis of the rational composition of the controlled parameters and an assessment of the required recovery time of the CTS. In this paper, the controlled states of the CTS will be evaluated using the probabilities of a false failure α , an undetected failure β and the time t B required for recovery after the failure is detected (the recovery time depends on the “severity” of the malfunction detected during monitoring). At low values of these estimated indicators, we will refer the CTS to the first group of degradation. As the degradation increases (as these indicators increase), we will refer the CTS to the second, third and fourth groups of degradation, respectively. Without going into the details of assigning parameters of criteria for attribution to a particular degradation group, we note that the number of degradation levels is determined by a set of types and types of CTS under consideration, their characteristic features, as well as the specific task being solved.
Figure 8a shows a graph with one fully workable state E 1 and four states corresponding to different levels of degradation (malfunction): E 2 , E 3 , E 4 and E 5 [6]. Let us distinguish the three parts in the classical model [2]: the initial operational state E 1 , the failure state E 2 and the subgraph corresponding to the control function (highlighted in Figure 8b by rectangle C2). Then, the classical model can be represented as a “serial connection” E 1 , C2 and E 2 [14,15]. The states of the subgraph: K 3 is verification of a failed MI, K 4 is the restoration of the CTS, K 5 is verification of a working MI and K 6 is the state of an undetected failure of the CTS. The probabilistic characteristics of the state transition are the same as in the classical model [2].
Note that if the control subgraph is completely removed from the graph in Figure 8b and the probabilistic characteristics are set on the edges of the remaining graph, then a simple model will be obtained that describes the operation of a small gun [6].
Figure 8a uses the notation: E1 is a fully functional state and four states corresponding to different levels of degradation of the CTS; E2 is the first group of degradation (functional state with minor deviations of the normalized metrological characteristics); E3 is the second group of degradation (a state with some deviations of the metrological characteristics, from which it is possible to return to a fully functional state with small resource costs); E4 is the third group degradation (a state from which it is possible to return to a fully functional state with costs associated with sufficiently resource-intensive maintenance); and E5 is the fourth “heavier” group of degradation. As the degradation group number increases, returning to the state E1 becomes more and more resource intensive.
Let us “attach” four metrological control systems, C2, C3, C4 and C5, between the fully functional state E 1 and the other four states, similar to the one shown in Figure 8a (“fan connection”). We will use the corresponding upper indices for the probabilistic and deterministic parameters of the model of each subsystem, describing samples of the CTS with different levels of degradation.
Then, the system of equations describing the semi-Markov stationary model will take the form [6]:
π 1 ( 1 ) = π 4 ( 1 ) + ( 1 α 1 ) π 5 ( 1 ) + π 4 ( 2 ) + ( 1 α 2 ) π 5 ( 2 ) + π 4 ( 3 ) + ( 1 α 3 ) π 5 ( 3 ) + π 4 ( 4 ) + ( 1 α 4 ) π 5 ( 4 ) π 2 ( 1 ) = γ 1 π 1 π 3 ( 1 ) = ( 1 δ 1 ) π 2 ( 1 ) + π 6 ( 1 ) π 4 ( 1 ) = ( 1 β 1 ) π 3 ( 1 ) + α 1 π 5 ( 1 ) π 5 ( 1 ) = 1 γ 1 π 1 π 6 ( 1 ) = β 1 π 3 ( 1 ) ,
π 2 ( 2 ) = γ 2 π 1 + δ 1 π 2 ( 1 ) π 3 ( 2 ) = ( 1 δ 2 ) π 2 ( 2 ) + π 6 ( 2 ) π 4 ( 2 ) = ( 1 β 2 ) π 3 ( 2 ) + α 2 π 5 ( 2 ) π 5 ( 2 ) = 1 γ 2 π 1 π 6 ( 2 ) = β 2 π 3 ( 2 ) ,   π 2 ( 3 ) = γ 3 π 1 + δ 2 π 2 ( 2 ) π 3 ( 3 ) = ( 1 δ 3 ) π 2 ( 3 ) + π 6 ( 3 ) π 4 ( 3 ) = ( 1 β 3 ) π 3 ( 3 ) + α 3 π 5 ( 3 ) π 5 ( 3 ) = 1 γ 3 π 1 π 6 ( 3 ) = β 2 π 3 ( 3 ) ,   π 2 ( 4 ) = γ 4 π 1 + δ 3 π 2 ( 3 ) π 3 ( 4 ) = π 2 ( 4 ) + π 6 ( 4 ) π 4 ( 4 ) = ( 1 β 4 ) π 3 ( 4 ) + α 4 π 5 ( 4 ) π 5 ( 4 ) = 1 γ 4 π 1 π 6 ( 4 ) = β 4 π 3 ( 4 )
Here, α i , ( i = 2 , 3 , 4 , 5 ) is the conditional probability of a false failure, β i is the conditional probability of an undetected failure and γ i = F i ( T K ) δ j , ( j = 1 , 2 , 3 ) are the probability of a transition from a state of degradation to the next, more severe, state number j + 1 .
Model (15) is a system of 21 equations. The rank of the system is 20. Exclude one of the equations (for example, the last equation of the system (15)) and add a normalization condition, as follows:
π 1 + i = 1 4 j = 2 6 π j ( i ) = 1
Then, the resulting system of linear inhomogeneous Equations (15) and (16) will have a unique solution that can be obtained using standard algorithms and methods for solving the corresponding systems [5,6].
Initial data: the total number of states is 21 and the number of degradation levels is four. As the CTS degrades, the duration of the verification and recovery time increase, and reliability decreases.
As generalized parameters characterizing the distribution of the control volumes by the degradation groups, the duration TIBV T K ( i ) for each of the four degradation groups was selected. As a result of the calculations, the dependence of the readiness coefficient on the TIBV was constructed:
K A = K A ( T K ( 1 ) , T K ( 2 ) , T K ( 3 ) , T K ( 4 ) ) ,
and the analysis of the influence of the TIBV of the different degradation groups (i = 1, 2, 3, 4) on the readiness coefficient was carried out. When constructing the dependence (17), the probabilities of false and undetected failures were set as average values for each of the degradation groups, namely: α 4 > α 3 > α 2 > α 1 , β 4 > β 3 > β 2 > β 1 .
The calculations have shown that if three arguments out of four are fixed in function (17), for example T K ( 2 ) = c ( 2 ) * , T K ( 3 ) = c ( 3 ) * , T K ( 4 ) = c ( 4 ) * , c ( i ) * = c o n s t then the dependence of function (17) on the remaining variable T K ( 1 ) will have the form shown in Figure 9. If two arguments out of four are fixed in function (17), for example T K ( 1 ) = c ( 1 ) * and T K ( 2 ) = c ( 2 ) * then the readiness coefficient K A , as functions of two variables, will be convex upwards (Figure 9). The maximum of the readiness coefficient K A is reached at a single internal point. Here and further, an asterisk in the upper index means that the corresponding value is set and fixed.
The calculations have shown that the maximum value of the readiness coefficient is achieved if the TIBV for the fourth degradation group is about 1.3 times less than for the third group.
The developed model allows us to calculate the optimal duration of the TIBV for the CTS of different degradation groups. If it is impossible to provide optimal TIBV values for some degradation groups in practice, then in (17) “possible” TIBV values should be set for these groups and local optimum TIBV durations for the remaining degradation groups should be calculated.
Note that the model of interaction of the CTS with the MI with a simplified form of technical condition control can be represented as a graph (Figure 8a), if you remove the control subgraphs C2–C5 and set the probabilistic characteristics of the state transitions on the edges of the remaining graph. Such a model, supplemented with a formula for calculating the average total resource costs S U M = S 12 p 12 + S 23 p 23 + S 34 p 34 + S 45 p 45 (where S 12 , S 23 , S 34 , S 45 are unit costs and p 12 , p 23 , p 34 , p 45 are probabilities of the state transitions), was used in [6] when calculating the technical and economic indicators of the metrological support system, when forming programs for the long-term development of the CTS fleet.
The models described in Section 4.1 and 4.2 do not allow modeling processes of CTS fleet renewal, and do not allow for the taking into account of the procurement of new CTS samples, or the modernization of existing CTS samples and the development of promising CTS samples. The model presented in Section 4.3 of the article allows modeling for all stages of the life cycle, including procurement, modernization and development of advanced CTS samples.

4.3. The Model of Operation of the Complex Technical System Fleet with a Partially Recoverable Resource

Next, we will distribute the CTS into three degradation groups [7]: the first is the start of operation of the CTS, the sample remains operational and the changes are insignificant; the second is where the operation and resource consumption of the CTS sample continue, the changes in characteristics are significant and the rate of change is average; and the third is a long-term operation, where the changes in characteristics are very significant and the rate of change is high. The number of degradation groups is determined by a set of types and the types of CTS under consideration, as well as the specific task being solved.
Figure 10 presents a graph of the operation model of the updated CTS fleet, with three degradation groups and two subgraphs modeling the process of updating the CTS fleet. The upper indices in parentheses indicate the number of the degradation group. Each degradation group will be modeled using the classical model [2], described in Section 4.1.1.
Each of the two subgraphs describing the upgrade process include three states: E 7 ( l ) ,   l = 2 ,   3 are the in-depth diagnostics of the technical condition; E 8 ( l ) ,   l = 2 ,   3 is the repair of the CTS; and E 9 ( l ) ,   l = 2 ,   3 is the purchase (or development and production) of a new similar model of the CTS. The probabilistic parameters of the main state transitions are shown in Figure 10, in Greek letters. Certain shares of the CTS, from the second ω(2) and third ω(3) degradation groups, in case of failure of the CTS are sent for in-depth diagnostics of the technical condition, in order to determine the feasibility of updating (replacing with a new model of the CTS) or continuing operation after repair. To simplify, some probabilistic characteristics are not indicated in Figure 10, but they can be easily restored, taking into account that the sum of the probabilities of the transitions from each vertex of the graph are equal to one. If one edge comes out of the vertex, then the corresponding transition probability is one, and if two edges come out of the vertex, and the probability of one transition is written on the graph, then the probability of the second transition is equal to the difference of one and the known probability of the first transition.
The model of the management of the CTS fleet is presented structurally in the form of four systems of equations:
π 1 ( 1 ) = π 4 ( 1 ) + ( 1 α ( 1 ) ) π 5 ( 1 ) + π 9 ( 2 ) + π 9 ( 3 ) π 2 ( 1 ) = γ ( 1 ) π 1 ( 1 ) π 3 ( 1 ) = ( 1 χ ( 1 ) ) π 2 ( 1 ) + π 6 ( 1 ) π 4 ( 1 ) = ( 1 β ( 1 ) ) π 3 ( 1 ) + α ( 1 ) π 5 ( 1 ) π 5 ( 1 ) = 1 γ ( 1 ) η ( 1 ) π 1 ( 1 ) π 6 ( 1 ) = β ( 1 ) π 3 ( 1 ) ,   π 1 ( 2 ) = π 4 ( 2 ) + ( 1 α ( 2 ) ) π 5 ( 2 ) + η ( 1 ) π 1 ( 1 ) + π 8 ( 2 ) π 2 ( 2 ) = γ ( 2 ) π 1 ( 2 ) + χ ( 1 ) π 2 ( 1 ) π 3 ( 2 ) = ( 1 χ ( 2 ) ) π 2 ( 2 ) + π 6 ( 2 ) π 4 ( 2 ) = ( 1 β ( 2 ) ω ( 2 ) ) π 3 ( 2 ) + α ( 2 ) π 5 ( 2 ) π 5 ( 2 ) = 1 γ ( 2 ) η ( 2 ) π 1 ( 2 ) π 6 ( 2 ) = β ( 2 ) π 3 ( 2 )
π 1 ( 3 ) = π 4 ( 3 ) + ( 1 α ( 3 ) ) π 5 ( 3 ) + η ( 2 ) π 1 ( 2 ) + π 8 ( 3 ) π 2 ( 3 ) = γ ( 3 ) π 1 ( 3 ) + χ ( 2 ) π 2 ( 2 ) π 3 ( 3 ) = π 2 ( 3 ) + π 6 ( 3 ) π 4 ( 3 ) = ( 1 β ( 3 ) ω ( 3 ) ) π 3 ( 3 ) + α ( 3 ) π 5 ( 3 ) π 5 ( 3 ) = 1 γ ( 3 ) π 1 ( 3 ) π 6 ( 3 ) = β ( 3 ) π 3 ( 3 ) ,   π 7 ( 2 ) = ω ( 2 ) π 3 ( 2 ) π 8 ( 2 ) = μ ( 2 ) π 3 ( 2 ) π 7 ( 3 ) = ω ( 3 ) π 3 ( 3 ) π 8 ( 2 ) = μ ( 3 ) π 3 ( 2 ) π 9 ( 2 ) = 1 μ ( 2 ) π 7 ( 2 ) π 9 ( 3 ) = 1 μ ( 3 ) π 7 ( 3 )
Here π i ( j ) are the stationary probabilities of finding the CTS in the corresponding states; α ( j ) ,   j = 1 ,   2 ,   3 , β ( j ) ,   j = 1 ,   2 ,   3 are the conditional probabilities of false and undetected failures, respectively; γ ( j ) = F j ( T K ) is the probability of failure during the time interval T K between the verifications; F j ( T ) is exponential distribution function; and ( χ ( j ) , η ( j ) , j = 1 , 2 ) are the probabilities of the transitions of the corresponding states from the j degradation group to the next ( j + 1 ) group. The first three systems (18) describe the processes of the CTS operation for the three degradation groups, and the fourth system (18) describes the process of updating the CTS fleet.
Model (18) is a homogeneous system of 24 linear algebraic equations. The rank of the system is 23. Exclude one of the equations (for example, the last equation of the system (18)) and add a normalization condition instead, as follows:
i = 1 6 j = 1 3 π i ( j ) + i = 7 9 j = 2 3 π i ( j ) = 1
Then, the resulting system of linear inhomogeneous algebraic equations (18) will have a unique solution [7].
The readiness coefficient K A of the fleet of the CTS is calculated by the formula [2]:
K A = j π 1 ( j ) ψ 1 ( j ) / i , j π i ( j ) ψ i ( j )
Here ψ i ( j ) is the mathematical expectation of the time (average time) of the CTS being in the corresponding states E i ( j ) (assumed to be known). In the numerator (20), summation by index j is performed for all workable states, and in the denominator (20) is the summation by both index i and index j for all states (the index i is responsible for unworkable states).
As parameters characterizing the distribution of the metrological control volumes and the quality of the metrological control by degradation groups, the duration of the TIVB T K j ,   j = 1 ,   2 ,   3 for each of the three degradation groups, the relative values of the operational tolerances for the controlled parameters δ j ,   j = 1 ,   2 ,   3 and the relative measurement errors z j ,   j = 1 ,   2 ,   3 , were selected.
As a result of the solution for system (18), the dependence of the CTS readiness coefficient for use on the above metrological parameters, organizational, technical and technical parameters is constructed:
K A = K A ( T K ( 1 ) , T K ( 2 ) , T K ( 3 ) , α ( 1 ) , α ( 2 ) , α ( 3 ) , β ( 1 ) , β ( 2 ) , β ( 3 ) , ω ( 2 ) , ω ( 3 ) ,   μ ( 2 ) , μ ( 3 ) ) ,
moreover, the functions of the conditional probabilities of false failures and undetected failures depend on the relative operational tolerance and relative measurement errors:
α ( 1 ) = α ( 1 ) ( δ ( 1 ) , z ( 1 ) ) ,   α ( 2 ) = α ( 2 ) ( δ ( 2 ) , z ( 2 ) ) ,   α ( 3 ) = α ( 3 ) ( δ ( 3 ) , z ( 3 ) ) , β ( 1 ) = β ( 1 ) ( δ ( 1 ) , z ( 1 ) ) ,   β ( 2 ) = β ( 2 ) ( δ ( 2 ) , z ( 2 ) ) ,   β ( 3 ) = β ( 3 ) ( δ ( 3 ) , z ( 3 ) )
that are calculated using Formulas (5)–(8).
In (21), the following parameters are presented: ω ( 2 ) , ω ( 3 ) are the proportion of samples sent for in-depth diagnostics of the CTS samples from the number of samples received for verification; and μ ( 2 ) , μ ( 3 ) are the parameters characterizing the process of updating the CTS fleet (so, for example, in special cases μ ( j ) = 0 , where all inoperable CTS samples are changed to new ones, and where μ ( j ) = 1 they are repaired). The parameters ω ( j ) , μ ( j ) are conditionally attributed to the organizational and technical categories. The readiness coefficient also depends on other technical parameters, for example η ( 1 ) , η ( 2 ) , χ ( 1 ) , χ ( 2 ) , which characterize the degradation process of the CTS fleet (operational parameters), and the average time ψ i ( j ) spent by the CTS sample in various states. These parameters are determined based on the processing of the available statistical information and the relevant criteria for classifying the CTS into different degradation groups.
Note that the parameters ψ 9 j (time spent in the state E 9 ( i ) ) allow you to model both the purchase and development of new samples of CTS. To simulate the procurement of new samples of CTS we have to set ψ 9 j sufficiently small, and to simulate the development of new samples, we have to set the CTS at medium and large.
Note that the constructed dependence (21), like (17), is smooth, so its extreme properties can be effectively investigated using standard gradient methods.
On the basis of solving a series of problems on the extremum of a function of several variables (21), the influence of the TIBV T K ( j ) of the CTS from different degradation groups and the relative tolerances on controlled parameters δ ( j ) on the readiness coefficient are analyzed K A .
Consider the effect of the duration of the TIBV on the readiness coefficient K A . Let us fix all the arguments (21), with the exception of three: T K i ,   i = 1 ,   2 ,   3 . If we additionally, fix any two arguments T K ( i ) out of three, for example T K 2 = C T K 2 * , T K 3 = C T K 3 * , then the dependence of function (21) on the remaining argument T K ( 1 ) will have the form given in [2]: convex upwards with a single maximum.
An asterisk in the upper index means that the corresponding value is set and fixed. If one of the three arguments is fixed in function (21) (for example, T K 3 = C T K 3 * ), then the readiness coefficient curve K A , as well as the functions of the other two arguments, T K ( 1 ) and T K ( 2 ) , will be convex upwards (Figure 11). The maximum of the readiness coefficient K A will be reached at a single internal point in the parameter plane, T K ( 1 ) × T K ( 2 ) . The dependences of the readiness coefficient K A on the frequency of the control for the first and third groups, and for the second and third groups of degradation, have a similar form.
Optimization (21) for three groups of degradation showed that the characteristic ratio of the TIBV durations is 80:45:30, thus, the higher the degradation group, the more often CTS verifications are required.
Consider the effect of the relative operating tolerances δ j ,   j = 1 ,   2 ,   3 on K A . Similarly to the above, we fix all the arguments (21), with the exception of three δ j ,   j = 1 ,   2 ,   3 . Then, the dependence of the readiness coefficient K A on relative operational tolerances is similar to its dependence on the TIBV, T K i ,   i = 1 ,   2 ,   3 .
The calculations have shown that the general form of dependence K A on two tolerances, at a fixed value of the third tolerance, has the form of surfaces shown in Figure 12. The surface of the readiness coefficient K A as a function of two arguments will be convex upwards, where the maximum is reached at a single internal point in the parameter plane, δ ( 1 ) × δ ( 2 ) , δ ( 1 ) × δ ( 3 ) or δ ( 2 ) × δ ( 3 ) . Optimization of the K A of three relative tolerances simultaneously showed that their characteristic ratio is 0.07:0.09:0.13, thus, the higher the degradation group, the greater the tolerance.
The study of the joint dependence K A on the TIBV and tolerances showed that the maximum of the function of six variables is achieved at a single internal point of a set of parameters, T K ( 1 ) × T K ( 2 ) × T K ( 3 ) × δ ( 1 ) × δ ( 2 ) × δ ( 3 ) . The optimal values of the arguments were: T K ( 1 ) = 81.59 , T K ( 2 ) = 50.46 , T K ( 3 ) = 32.55 , δ ( 1 ) = 0.072 , δ ( 2 ) = 0.0872 and δ ( 3 ) = 0.128 . At the same time, the optimal values of the probabilities of false and undetected failures were: α ( 1 ) = 0.227 , α ( 2 ) = 0.267 , α ( 3 ) = 0.347 , β ( 1 ) = 0.171 , β ( 2 ) = 0.222 and β ( 3 ) = 0.323 . The calculations have shown that with an increase in the number of the degradation group, the TIBV decreases, the tolerances for controlled parameters and the probabilities of false and undetected failures increase. The probabilities of false failures slightly exceed the corresponding probabilities of undetected failures for each degradation group.
We describe the results of a study on the stationary distribution of the CTS samples in different degradation groups, depending on the rate of degradation processes. The rate of degradation is determined using transition probabilities χ ( i ) , η ( i ) . The lower the corresponding probability, the slower the degradation processes proceed. Four variants differing in the rate of degradation were investigated: η ( 1 ) = 0.25 , χ ( 1 ) = 0.2 , η ( 2 ) = 0.35 , χ ( 2 ) = 0.3 (option 1); η ( 1 ) = 0.025 , χ ( 1 ) = 0.02 , η ( 2 ) = 0.35 , χ ( 2 ) = 0.3 (option 2); η ( 1 ) = 0.025 , χ ( 1 ) = 0.02 , η ( 2 ) = 0.035 , χ ( 2 ) = 0.03 (option 3); and η ( 1 ) = 0.025 , χ ( 1 ) = 0.02 ; η ( 2 ) = 0.035 , χ ( 2 ) = 0.003 (option 4). Note that the variants are arranged in order of decreasing degradation rate. The distribution of the proportion of working samples of the CTS by degradation groups at different values of these parameters is shown in Figure 13.
The probability of the CTS staying in the first degradation group for option 2 is about 5.5 times higher compared to option 1. At the same time, the ratio of the probability of being in the third group compared to the probability of being in the first group remains approximately the same.
The probability of the CTS staying in the first degradation group monotonically decreases, and the probability of being in the second group monotonically increases with sequential consideration of options from 2 to 4.
Next, we investigate the dependence K A on the total production capacity of the metrological units in which the MI and CTS are verified and checked. The production capacity of a metrological unit may be temporarily limited for one reason or another. The specified restriction was set in the form of an inequality, j π j ψ j C ( π ) * , where C ( π ) * is the conditional production capacity of the metrological unit, and the problem of conditional optimization was solved. In Figure 14 the dependences of the TIBV on the total conditional production capacity ζ of the metrological units and the readiness coefficient corresponding to these intervals (in percent) are presented.
If the production capacities of the metrological division do not allow for the checking of the required number of CTS, then it is possible to operate a fleet of CTS with increased TIBV. With a decrease in the production capacity of the metrological unit from 100% to 75%, the readiness coefficient K A decreases from 0.9514 to 0.8402.

5. Results

A set of interrelated mathematical models of the processes of operation, renewal and degradation of a fleet of CTS with metrological support was developed. The basis of the developed set of models consists of:
A basic model of the CTS operation;
A set of CTS operation models, having different levels of degradation (for different levels of CTS degradation a different number of system states and different variants of system maintenance are used);
A model of false failures and undetected failures;
A model of CTS fleet renewal, including such renewal methods as the purchase of new CTS samples, the modernization of existing CTS samples and the development of new promising CTS samples;
A functional dependence model of the CTS availability factor on a number of technical parameters, organizational and technical parameters, and technological parameters of the CTS belonging to different degradation groups and different methods of CTS stock renewal.
On the basis of the set of interrelated mathematical models presented in the article, the software for modeling the processes of operation, renewal and degradation of the fleet of CTS with metrological support was developed.

6. Discussions

The models developed and implemented as software allow for the parametrical optimization of the processes of CTS fleet functioning for a number of parameters, including metrological parameters, organizational and technical parameters, and technical parameters.
If in practice it is impossible to provide the optimal TIBV values or tolerances for the controlled parameters for some degradation groups, then for these groups the “possible” values of the TIBV and tolerances for the controlled parameters should be established, and the developed models should be used for calculation.
The developed set of models include a model for calculating the probabilities of false failures and undetected failures, for use in cases where the measured parameter and measurement error have a normal distribution law. The set of models developed in the article can also be used for different distribution laws of the measured parameter and measurement error: analytical defined laws, statistical defined laws or analytical and statistical distribution laws.
The constructed functional dependences of the availability factor on metrological, technical, organizational, technical and technological parameters have a smooth character, which makes it possible to effectively investigate the extreme properties of the availability factor using standard gradient methods.

7. Conclusions

Thus, the research goal has been achieved: a set of interrelated models has been developed, which solves the current need for end-to-end modeling of all the main stages of the life cycle of the CTS fleet. The developed set of models makes it possible to adequately simulate large CTS fleets, including those incorporating up to several hundreds of thousands of CTS samples.
The developed set of interrelated models allows:
Management of the process of development of CTS fleets;
Optimization of the processes of CTS fleet functioning;
Identification of problematic issues in the development of CTS fleet and the formation of strategies for CTS fleet development in the presence of various constraints;
The solving of the problem of conditional optimization in the presence of constraints on the technological parameters of the CTS fleet development (with constraints on part of the arguments of the availability factor function);
Calculation of the technological and technical–economic parameters of the CTS fleet functioning and development;
Evaluation of the risks associated with false and undetected failures, as well as the risks associated with CTS degradation;
A set of models is used in the Main Scientific Metrology Center:
To classify the designed CTS in order to establish the requirements for their metrological support;
When developing plans for medium-term and long-term development of the CTS fleet.
A set of models and software can be used by design organizations involved in the development of modern and advanced CTS with metrological support.

8. Future Works

At present, a set of models continues to develop in the direction of development and replenishment with models of CTS fleet maintenance; namely, models of workplaces for the verification of MI, taking into account the priorities of the MI samples coming for verification [27,28], as well as models of CTS fleet functioning under such modes of functioning as the mode of high readiness for use, the mode of use in extreme conditions, etc.

Author Contributions

The material is the result of the joint efforts of R.K., D.E., A.M. and T.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of Moscow State University of Civil Engineering (protocol № 1, 5 December 2022).

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Volkov, L.I. Aircraft Operation Management Moscow: Higher School. 1981. 368p. Available online: https://www.studmed.ru/volkov-l-i-upravlenie-ekspluataciey-letatelnyh-kompleksov_e28fc04269b.html (accessed on 6 March 2023).
  2. Sychev, E.I. Metrological Support of Radio-Electronic Equipment (Methods of Analysis). Moscow: RIC “Tatiana’s Day”. 1994. 277p. Available online: https://rusneb.ru/catalog/000199_000009_001691017/ (accessed on 6 March 2023).
  3. Mishchenko, V.I. Evolution of Models of the Process of Operation of Weapons and Military Equipment. Bull. Acad. Mil. Sci. 2003, 4, 200. Available online: https://militaryarticle.ru/vestnik-akademii-voennykh-nauk/2003-vavn/10460-jevoljucija-modelej-processa-jekspluatacii (accessed on 6 March 2023).
  4. Mishchenko, V.I.; Kravtsov, A.N.; Mamleev, T.F. The Semi-Markov Model of Functioning of Redundant Measuring Instruments Taking into Account the Frequency of Verification. Meas. Equip. 2021, 4, 22–27. Available online: https://link.springer.com/article/10.1007/s11018-021-01931-3 (accessed on 6 March 2023).
  5. Khayrullin, R.Z.; Popenkov, A.J. Distribution of Controlling Volumes of Metrological Support for the Objectives of Complex Organizational and Technical Systems with the Use of Semi-Markov Models. In Proceedings of the XI-th International Conference “Management of Large-Scale System Development” (MLSD), Moscow, Russia, 1–3 October 2018; p. 8551917, 5p. Available online: https://www.researchgate.net/publication/329292364_Distribution_of_controlling_volumes_of_metrological_support_for_the_objectives_of_complex_organizational_and_technical_systems_with_the_use_of_semi-Markov_models (accessed on 6 March 2023).
  6. Popenkov, A.J.; Fufaeva, O.V.; Khayrullin, R.Z. Analysis of Operation Models of Complex Technical Systems with Degrading Metrological Support. Proc. Tula State Univ. 2022, 7, 247–254. Available online: https://cyberleninka.ru/article/n/analiz-modeley-ekspluatatsii-slozhnyh-tehnicheskih-sistem-s-degradiruyuschim-metrologicheskim-obespecheniem (accessed on 6 March 2023).
  7. Khayrullin, R.Z. Optimization of the Processes of Operation and Renewal of the Park of Measuring Equipment. Meas. Equip. 2022, 8, 28–34. Available online: https://www.izmt.ru/note.php?type=TAMI_izmt&notes_id=892 (accessed on 6 March 2023).
  8. Khayrullin, R.Z.; Kornev, A.S.; Kostoglotov, A.A.; Lazarenko, S.V. Mathematical Simulation of Decision Error Functions in Tolerance Control of the Performance of Measuring Techniques. Meas. Tech. 2020, 63, 680–685. Available online: https://www.researchgate.net/publication/348442252_Mathematical_Simulation_of_Decision_Error_Functions_in_Tolerance_Control_of_the_Performance_of_Measuring_Techniques (accessed on 6 March 2023). [CrossRef]
  9. David, V.; Marie, F.; Zdeněk, V.; Jakub, G. Degradation Process and Failure Estimation of Drilling System Based on Real Data and Diffusion Process Supported by State Space Models. Measurement 2020, 164, 108076. Available online: https://www.researchgate.net/publication/342137334_Degradation_process_and_failure_estimation_of_drilling_system_based_on_real_data_and_diffusion_process_supported_by_state_space_models (accessed on 6 March 2023).
  10. Ma, J.; Fouladirad, M.; Grall, A. Flexible Wind Speed Generation Model: Markov Chain with an Embedded Diffusion Process. Energy 2018, 164, 316–328. Available online: https://www.researchgate.net/publication/327381227_Flexible_wind_speed_generation_model_Markov_chain_with_an_embedded_diffusion_process (accessed on 6 March 2023). [CrossRef]
  11. Azarskov, V.N.; Strelnikov, V.P. Reliability of Control Systems and Automation. Kiev. NAU. 2004. 164p. Available online: https://www.studmed.ru/azarskov-vn-strelnikov-vp-nadezhnost-sistem-upravleniya-i-avtomatiki_881a18fe6fa.html (accessed on 6 March 2023).
  12. Khayrullin, R.Z. Methods and Models of Calibration of Vector Network Analyzers. J. Phys. Conf. Ser. 2021, 1921, 012122. [Google Scholar] [CrossRef]
  13. Abinandhitha, R.; Sakthivel, R.; Kong, F.; Parivalla, A. Robust Non-Fragile Boundary Control for Non-Linear Parabolic PDE Systems with Semi-Markov Switching and Input Quantization. Eur. J. Control 2022, 67, 100713. Available online: https://www.researchgate.net/publication/362087099_Robust_non-fragile_boundary_control_for_non-linear_parabolic_PDE_systems_with_semi-Markov_switching_and_input_quantization (accessed on 6 March 2023). [CrossRef]
  14. Jiang, D.; Chen, T.; Xie, J.; Cui, W.; Song, B. A Mechanical System Reliability Degradation Analysis and Remaining Life Estimation Method with the Example of an Aircraft Hatch Lock Mechanism. Reliab. Eng. Syst. Saf. 2023, 230, 108922. Available online: https://www.researchgate.net/publication/365098134_A_mechanical_system_reliability_degradation_analysis_and_remaining_life_estimation_method--With_the_example_of_an_aircraft_hatch_lock_mechanism (accessed on 6 March 2023). [CrossRef]
  15. Yan, X.; Tang, G.; Wang, X. Bearing Performance Degradation Assessment Based on the Continuous-Scale Mathematical Morphological Particle and Feature Fusion. Measurement 2022, 188, 110571. Available online: https://www.researchgate.net/publication/356855692_Bearing_Performance_Degradation_Assessment_Based_on_the_Continuous-scale_Mathematical_Morphological_Particle_and_Feature_Fusion (accessed on 6 March 2023). [CrossRef]
  16. Yan, T.; Lei, Y.; Li, N.; Wang, B.; Wang, W. Degradation Modeling and Remaining Useful Life Prediction for Dependent Competing Failure Processes. Reliab. Eng. Syst. Saf. 2021, 212, 107638. Available online: https://www.researchgate.net/publication/350120821_Degradation_Modeling_and_Remaining_Useful_Life_Prediction_for_Dependent_Competing_Failure_Processes (accessed on 6 March 2023). [CrossRef]
  17. Fang, G.; Pan, R.; Hong, Y. Copula-Based Reliability Analysis of Degrading Systems with Dependent Failures. Reliab. Eng. Syst. Saf. 2020, 193, 106618. Available online: https://www.researchgate.net/publication/335646674_Copula-based_Reliability_Analysis_of_Degrading_Systems_with_Dependent_Failures (accessed on 6 March 2023). [CrossRef]
  18. Romero, E.E.; Bérenguer, C.; Martinez, J.J. Long-Term Degradation Estimation of Wind Turbine Drive-Train under a Gain-scheduling Control Strategy According to the Weather Conditions. IFAC-Pap. 2022, 55, 13–18. Available online: https://www.researchgate.net/publication/363895260_Long-term_degradation_estimation_of_wind_turbine_drive-train_under_a_gain-scheduling_control_strategy_according_to_the_weather_conditions (accessed on 6 March 2023). [CrossRef]
  19. Jantara, V.L.; Basoalto, H.; Papaelias, M. A damage mechanics approach for lifetime estimation of wind turbine gearbox materials. Int. J. Fatigue 2020, 137, 105671. [Google Scholar] [CrossRef]
  20. Komusanac, I.; Brindley, G.; Fraile, D.; Ramirez, L. Wind Energy in Europe; 2020 Statistics and the Outlook for 2021–2025; Technical Report; WindEurope Business Intelligence: Brussels, Belgium, 2021; Available online: https://www.windenergyhamburg.com/fileadmin/windenergy/2022/pdf/we22_wind-europe-stats2020.pdf (accessed on 6 March 2023).
  21. Lee, J.; Zhao, F. GWEC–Global Wind Report 2021. Technical Report. Global Wind Energy Council. 2021. Available online: https://gwec.net/global-wind-report-2021/ (accessed on 6 March 2023).
  22. Merainani, B.; Laddada, S.; Bechhoefer, E.; Chikh, M.; Benazzouz, D. An integrated methodology for estimating the remaining useful life of high-speed wind turbine shaft bearings with limited samples. Renew. Energy 2022, 182, 1141–1151. [Google Scholar] [CrossRef]
  23. Romero, E.E.; Martinez, J.J.; Berenguer, C. Degradation of a Wind-Turbine Drive-Train under Turbulent Conditions: Effect of the Control Law. In Proceedings of the 2021 5th International Conference on Control and Fault-Tolerant Systems (SysTol), Saint-Raphael, France, 29 September–1 October 2021; pp. 335–340. Available online: https://www.researchgate.net/publication/356095568_Degradation_of_a_wind-turbine_drive-train_under_turbulent_conditions_effect_of_the_control_law (accessed on 6 March 2023).
  24. Tchakoua, P.; Wamkeue, R.; Ouhrouche, M.; Slaoui, H.; Tameghe, T.; Ekemb, G. Wind Turbine Condition Monitoring: State-of-the-Art Review, New Trends, and Future Challenges. Energies 2014, 7, 2595–2630. Available online: https://www.researchgate.net/publication/273984269_Wind_Turbine_Condition_Monitoring_State-of-the-Art_Review_New_Trends_and_Future_Challenges (accessed on 6 March 2023). [CrossRef]
  25. Kampen, J.K. Reflections on and test of the metrological properties of summated rating, Likert, and other scales based on sums of ordinal variables. Measurement 2019, 137, 428–434. Available online: https://www.researchgate.net/publication/330828499_Reflections_on_and_test_of_the_metrological_properties_of_summated_rating_Likert_and_other_scales_based_on_sums_of_ordinal_variables (accessed on 6 March 2023). [CrossRef]
  26. Chen, H.B.; Zhuang, H.L. A new highly anti-interference regularization method for ill-posed problems. Vibroeng. Procedia 2017, 15, 128–133. [Google Scholar] [CrossRef]
  27. Ershov, D.S.; Khayrullin, R.Z. Workplace for Verification of Measuring Instruments as a Non-Stationary Service System. Instrum.-Mak. J. 2022, 65, 701–711. Available online: https://cyberleninka.ru/article/n/matematicheskaya-model-rabochego-mesta-poverki-sredstv-izmereniy-kak-nestatsionarnaya-sistema-obsluzhivaniya (accessed on 6 March 2023).
  28. Ershov, D.S.; Malahov, A.V.; Levina, T.A.; Khayrullin, R.Z. To Simulating the Functioning of Workplace for the Verification of Measuring Instruments. J. Phys. Conf. Ser. 2022, 2388, 012147. Available online: https://www.researchgate.net/publication/366228562_To_simulating_the_functioning_of_workplace_for_the_verification_of_measuring_instruments (accessed on 6 March 2023). [CrossRef]
Figure 1. Graph of state transitions.
Figure 1. Graph of state transitions.
Axioms 12 00300 g001
Figure 2. Dependences of the readiness coefficients on the periodicity of the control for various distribution laws: Rayleigh’s law (1), normal law (2), exponential law (3), and Weibull’s law (4).
Figure 2. Dependences of the readiness coefficients on the periodicity of the control for various distribution laws: Rayleigh’s law (1), normal law (2), exponential law (3), and Weibull’s law (4).
Axioms 12 00300 g002
Figure 3. General scheme of diagnostics and decision-making based on the tolerance control method.
Figure 3. General scheme of diagnostics and decision-making based on the tolerance control method.
Axioms 12 00300 g003
Figure 4. Dependences of the probability of an erroneous decision ( α + β ) , as well as the probabilities of false and undetected failures on the value of the reduced tolerance δ on z = 0.5 .
Figure 4. Dependences of the probability of an erroneous decision ( α + β ) , as well as the probabilities of false and undetected failures on the value of the reduced tolerance δ on z = 0.5 .
Axioms 12 00300 g004
Figure 5. Probability of false failures (a); probability of undetected failures (b).
Figure 5. Probability of false failures (a); probability of undetected failures (b).
Axioms 12 00300 g005
Figure 6. The model of a random degradation process and a scheme for the formation of a time-to-failure distribution: (a) α distribution (fan process); (b) DN distribution law; (c) DM distribution law.
Figure 6. The model of a random degradation process and a scheme for the formation of a time-to-failure distribution: (a) α distribution (fan process); (b) DN distribution law; (c) DM distribution law.
Axioms 12 00300 g006
Figure 7. Diffusion distribution functions: (a) DM distribution; (b) DN distribution.
Figure 7. Diffusion distribution functions: (a) DM distribution; (b) DN distribution.
Axioms 12 00300 g007
Figure 8. Graphs and subgraphs of the CTS operation model with full resource recovery: (a) with the control of four degradation states; (b) with the control of one degradation state (classical model).
Figure 8. Graphs and subgraphs of the CTS operation model with full resource recovery: (a) with the control of four degradation states; (b) with the control of one degradation state (classical model).
Axioms 12 00300 g008
Figure 9. Dependence of the readiness coefficient K A on the TIBV for technical conditions E 3 and E 4 .
Figure 9. Dependence of the readiness coefficient K A on the TIBV for technical conditions E 3 and E 4 .
Axioms 12 00300 g009
Figure 10. Graph of the model of operation, degradation and renewal of the CTS fleet with incomplete resource recovery.
Figure 10. Graph of the model of operation, degradation and renewal of the CTS fleet with incomplete resource recovery.
Axioms 12 00300 g010
Figure 11. Dependence of the readiness coefficient K A on the TIBV for the first and second degradation groups.
Figure 11. Dependence of the readiness coefficient K A on the TIBV for the first and second degradation groups.
Axioms 12 00300 g011
Figure 12. The dependence of the readiness coefficient K A on the relative operating tolerances for the controlled parameters: (a) on δ ( 1 ) and δ ( 2 ) ; (b) on δ ( 2 ) and δ ( 3 ) .
Figure 12. The dependence of the readiness coefficient K A on the relative operating tolerances for the controlled parameters: (a) on δ ( 1 ) and δ ( 2 ) ; (b) on δ ( 2 ) and δ ( 3 ) .
Axioms 12 00300 g012
Figure 13. Distribution of workable CTS samples by degradation levels at different values of parameters determining the rate of degradation.
Figure 13. Distribution of workable CTS samples by degradation levels at different values of parameters determining the rate of degradation.
Axioms 12 00300 g013
Figure 14. Dependence of the readiness coefficient and the TIBV on the total production capacity of metrological units.
Figure 14. Dependence of the readiness coefficient and the TIBV on the total production capacity of metrological units.
Axioms 12 00300 g014
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

Khayrullin, R.; Ershov, D.; Malahov, A.; Levina, T. Mathematical Models of the Processes of Operation, Renewal and Degradation of a Fleet of Complex Technical Systems with Metrological Support. Axioms 2023, 12, 300. https://doi.org/10.3390/axioms12030300

AMA Style

Khayrullin R, Ershov D, Malahov A, Levina T. Mathematical Models of the Processes of Operation, Renewal and Degradation of a Fleet of Complex Technical Systems with Metrological Support. Axioms. 2023; 12(3):300. https://doi.org/10.3390/axioms12030300

Chicago/Turabian Style

Khayrullin, Rustam, Denis Ershov, Alexander Malahov, and Tatyana Levina. 2023. "Mathematical Models of the Processes of Operation, Renewal and Degradation of a Fleet of Complex Technical Systems with Metrological Support" Axioms 12, no. 3: 300. https://doi.org/10.3390/axioms12030300

APA Style

Khayrullin, R., Ershov, D., Malahov, A., & Levina, T. (2023). Mathematical Models of the Processes of Operation, Renewal and Degradation of a Fleet of Complex Technical Systems with Metrological Support. Axioms, 12(3), 300. https://doi.org/10.3390/axioms12030300

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