Control and Diagnostics of TV3‐117 Aircraft Engine Technical State in Flight Modes Using the Matrix Method for Calculating Dynamic Recurrent Neural Networks Serhii Vladova, Yurii Shmelova and Ruslan Yakovlieva a Kremenchuk Flight College of Kharkiv National University of Internal Affairs, vul. Peremohy, 17/6, Kremenchuk, Ukraine, 39605 Abstract In this paper, one of the methods of using dynamic recurrent neural networks for solving applied problems of control and diagnosing of aircraft engines technical state, including TV3-117, is considered. A method for coding neural networks into signal graphs is proposed, and it is shown that their adjacency matrices can be used as associative memory in step matrix algorithms for solving dynamic recurrent neural networks. It is shown that in fully connected recurrent neural networks, any neuron can be input or output, and one neuron can simultaneously be input and output. Examples of teaching by the evolutionary optimization algorithm for multiextremal problems of recurrent dynamic neural networks intended for control and diagnosing of TV3-117 aircraft engine technical state are given. The functions of activation blocks of neurons in dynamic recurrent neural networks in this work are used difference expressions of simulation models of linear dynamic links. It is shown that for identification in the time domain of transient processes in dynamic systems of the third order, satisfactory accuracy is achieved at the output of any neuron of a recurrent dynamic neural network with four neurons, while it was found that useful information about the dynamic properties of the dynamic system under study can be simultaneously obtained from the output of any neuron network. Keywords 1 Aircraft engine, neural network, signal graph adjacency matrices, synaptic weights, associative memory. 1. Introduction Flight safety of aircraft, including helicopters, is one of the key problems of aviation, largely depends on the reliability of the operation of their power plants. In turn, the reliable operation of power plants and, in particular, a gas turbine engine (GTE) (for example, TV3-117, which is part of the power plant of the Mi-8MTV helicopter and its modifications) is ensured by a whole range of measures, among which an important place is given to its diagnostics technical state. At present, the technical diagnostics of a GTE, in the broad sense of this concept, is carried out only on the ground by an engineering and technical staff, who have at their disposal the appropriate tools and methodological apparatus. But most of the failures that affect the safety of a helicopter operation arise and manifest themselves in flight, which requires an effective on-board diagnostic system that determines the technical condition of the gas turbine engine in real time directly on board the helicopter. The difficulty of carrying out diagnostics of the power plant in automatic mode in flight is associated with the high complexity of aircraft GTE, due to the multiparametry, multi-connectivity, nonlinearity of the processes occurring in them, the multi-mode application, which requires significant machine and time resources [1]. CMIS-2021: The Fourth International Workshop on Computer Modeling and Intelligent Systems, April 27, 2021, Zaporizhzhia, Ukraine EMAIL: ser26101968@gmail.com (S. Vladov); nviddil.klk@gmail.com (Yu. Shmelov); klk.univd@gmail.com (R. Yakovliev) ORCID: 0000-0001-8009-5254 (S. Vladov); 0000-0002-3942-2003 (Yu. Shmelov); 0000-0002-3788-2583 (R. Yakovliev) © 2021 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR Workshop Proceedings (CEUR-WS.org) Today engines of the 4th generation have on-board parameter control systems, which solve a wide range of tasks related to the assessment of the residual life of aircraft engines and control of their parameters. In the process of performing the task of engine life development, the problem of forecasting (short-term) the current state of the engine and trend analysis of its parameters is solved. In modern on-board system for control and diagnosing the following tasks of control of a large number of parameters of the engine [1] are solved: temperature and pressure of gases behind and in front of the turbine; rotational speed of the compressor rotor of low and high pressure; the efficiency of the oil system and engine prompting on the limit values of the parameters of the engine oil system; oil consumption to change the level in the engine oil tank; efficiency of the engine fuel system; vibration level; temperature in the engine cavities; bearing temperature; mass fuel consumption; temperature and air pressure behind the compressor; technical state of the flowing part of the engine according to its thermogas-dynamic parameters [2, 3]. However, with a wide variety of tasks, there are some difficulties: limited computing resources (amount of RAM, speed, accuracy of the results); difficulties in formalizing classical control and diagnostic algorithms and, as a consequence, their practical implementation; the need to implement the above algorithms in low-level language (assembler); difficulties in recovering information when sensors fail. Taking into account the specific features of the diagnostic object (aircraft GTE), as a complex technical system with essential nonlinear characteristics, the above difficulties can be effectively eliminated by using the mathematical apparatus of artificial neural networks, which have a number of advantages given in [4]. 2. Mathematical model of the change of aircraft engine TV3‐117 technical state To study the possible of TV3-117 aircraft engine technical states according to [5, 6], it is assumed that the engine in the helicopter flight mode can be in four states: S0 – engine idling (here we mean idling), S1 – engine works in the nominal mode (here we mean takeoff, nominal, cruising I, cruising II), S2 – the engine is working properly, readjustment is in progress, S3 – the engine is faulty (emergency mode of the engine). According to these data, it is possible to obtain a discrete Markov network based on the following possible transition states: S0 – S1, S1 – S0, S0 – S3, S1 – S3, S2 – S0, S2 – S1, S2 – S3, S3 – S2, the probability of the engine being in each of the above states must be taken into account (fig. 1) [5, 6]. P11 P22 P21 S0 S1 P12 P31 P14 P32 P24 P33 S P43 S3 2 P44 P34 Figure 1: Markov discrete network showing all kinds of changes of TV3‐117 aircraft engine technical state as a part of a helicopter power plant during flight tests [5, 6] According to fig. 1, transition probabilities matrix has the form:  P11 P12 0 P14    P P22 0 P24  P   21 . (1)  P31 P32 P33 P34     0 0 P43 P44  To find the probabilities of the stay of Markov chain in certain states as n → ∞ (final probabilities) in [5, 6] an equations system of the form is solved:  P11 1  P21 2  P31 3   1 ; P   P   P    ;  12 1 22 2 32 3 2 (2)  P  33 3  P  43 4   3 ;  P14 1  P24 2  P34 3  P44 4   4 . According to [5, 6], in matrix form, system (2) has the form π = PTπ and (PT – E)π = 0, where  P11 P21 P31 0  1 0 0 0     P P22 P32 0  0 1 0 0 PT   12 – transposed matrix; E   – identity matrix.  0 0 P33 P43   0 0 1 0      P14 P24 P34 P44  0 0 0 1 According to [5, 6], adding to this equation the normalization condition of the form: π1 + π2 + π3 + π4 = 1; (3)  P11  1 P21 P31 0    P P22  1 P32 0  the final matrix is obtained, PT  E   12 .  0 0 P33  1 P43     P14 P24 P34 P44  1 Replacing the fourth row of the matrix (PT – E) with a single row, that  P11  1 P21 P31 0     P12 P22  1 P32 0  is, P  E  T , a final system of linear algebraic equations is obtained:  0 0 P33  1 P43     1 1 1 1   P11  1 1  P21 2  P31 3  0;   P12 1   P22  1  2  P32 3  0;  (4)  P33  1  3  P43 4  0;         1.  1 2 3 4 As a result of solving (4), expressions were obtained for determining the elements of the vector of state probabilities π(π1, π2, π3, π4) (shows the probability that the engine will be in the i-th state), which will allow in the presence of values of the transition probabilities P11…P33 to obtain prognostic indicators of changes of TV3-117 aircraft engine technical state. As can be seen from system (4), to determine the probability of engine failure π4, it is enough to know the transition probabilities for the other three states, as well as the transition probabilities P34, P43, which can be obtained as a result of technical diagnostics of the helicopter. The adequacy of the developed method is confirmed in [5] by definition. The paper proposes an implementation of this method using dynamic recurrent artificial neural networks. 3. Matrix method for determining dynamic recurrent artificial neural networks In this paper, the topology (internal connections) of the neural network is proposed to be encoded by signal graphs [7, 8]. A model of a neuron of a neural network is proposed to be represented by a node in the signal graph with the assignment of a serial number 1, 2, ..., N – 1, N. The nodes are connected by directed arcs. The arcs are identified by synaptic weights wij, where i – number of the donor neuron, j – number of the acceptor neuron. A convenient descriptor of a neural network is the adjacency matrix of the signal graph W. Elements of the adjacency matrices are the storage of the associative memory of neural network for various purposes. The fig. 2 shows the signal graphs and the corresponding adjacency matrices: for the forward-directed and recurrent neural network. From the theory of graphs it follows that if the column of the adjacency matrix is zero (fig. 2, a), then the vertex of the signal graph is a sink, and in a neural network this neuron is an input one. The zero-line points to the top of the source of the signal graph, and in the neural network this neuron is the output. In a recurrent neural network (fig. 2, b) there is no formal concept of layers. Any neuron can be both input and output. Outputs 1 2 3 4 S1 1 0 w12 w13 w14  w12 w23 2  0 0 w23 0  Intputs w13 S0 S2 w14 w43 3 0 0 0 0  S3   4 0 0 w43 0  a) Outputs 1 2 3 4 w21 S1 w32 1 0 w12 w23 w12 w13 w14  2  w21 0 w23 0  Intputs w13 S0 S2 w14 w43 3 0 w32 0 w34  w41 S3 w34   4  w41 0 w43 0  b) Figure 2: Signal graphs and adjacency matrices: a – direct neural network; b – recurrent neural network The most widespread at present are direct-directional neural networks – perceptron’s [9–12], which are currently widely used in technology to solve problems of classifying the state of equipment and diagnosing events. For training forward-directed neural networks with a large number of neurons, there is a layer-by-layer backpropagation method [9–12]. For modeling and synthesis of control systems, dynamic recurrent neural networks are more attractive (fig. 2, b). Due to intranet feedbacks, to calculate each neuron in a recurrent neural network, signal values are required not only from the associated outputs of the neurons of previous layers, as in forward-directed neural networks, but also the values of signals received on this neuron itself, and others associated with him neurons. Layer-by-layer calculation of recurrent neural networks becomes impossible, in which each neuron can be both input and output at the same time. Currently, there are no formalized algorithms for calculating recurrent neural networks, including those designed to simulate dynamic systems. The known algorithms aim at adapting the backpropagation algorithm for calculating forward- directed neural networks with feedback [13, 14]. Dynamic properties (associative memory) in forward-directed neural networks are achieved by applying input signals to various layers of the neural network and increasing the total number of neurons. An analogy is seen with neural networks of existing structures of aggregate and controller control algorithms, and in the control algorithms there are real nonlinear elements (analogs of activation blocks). But they have a very small number of neurons (functional blocks). This necessitates, in contrast to multi-neural networks of scanners-classifiers, the study of dynamic neural networks from below, from a small number of neurons. And their complexity, topology and the number of neurons, must meet the requirements of the assigned tasks – high-quality implementation of various types of classical, adaptive, single-channel and hybrid multichannel, linear with switching and non-linear algorithms for regulation and logic control. Studies have shown that the listed tasks can be implemented by dynamic recurrent neural networks, for which there are currently no formalized calculation algorithms. In this article, we propose a matrix algorithm for the numerical step-by-step calculation of recurrent neural networks. In the proposed algorithm, the models of neural network neurons are represented by static equations of linear adders [15]: N s jk   xik  wij ; j = 1, 2, …, N (5) i 1 and activation block γn = φ(snk); (6) where x1n, x2n, …, xNn – signal values at the input of the n-th neuron; w1n, w2n, …, wNn – synaptic weights of channels at the input of the n-th neuron; (snk) – value of the signal at the output of the adder at the k-th step of the neural network calculation. In the processes of training or operation of a neural network, the proposed matrix algorithm requires [15]: 1) adjacency matrix of the signal graph of the neural network W filled with the values of synaptic weights; 2) vector of functions of activation blocks Ф(S): Ф  S   1  s1  ,2  s2  ,..., N 1  sN 1  , N  sN  ; (7) 3) vector of values supplied to the inputs of all neurons of the neural network of external signals in X : 0 X 0in   x10in , x20 in ,..., xNin10 , xNin0  ; (8) 4) vector of signal values at the outputs of the activation blocks: Г 0   10 ,  20 ,...,  N 10 ,  N 0  . (9) in In the process of calculating neural networks, the elements of vectors X and Г change from step to step. Suppose at the k-th step they have values X kin1 and Гk–1. Then the signals that at the kth step must be fed to the inputs of all neurons are summed up from vectors X kin1 and Гk–1 [15]: X k  X kin  Гk 1 . (10) The procedure for calculating the signal at the output of the adder of the j-th neuron of the neural network in accordance with (5) consists in multiplying the elements of the columns of the adjacency matrix W by the elements of the vector Xk, followed by summing the multiplication results with the elements of the j-th column of this matrix. It follows from the above that at the kth step, the neural network calculation procedure is reduced to three computation stages. At the first stage, using the values of the input signals X kin1 and the results of calculating the neural network at the (k – 1)-th step Гk–1, from expression (10), the elements of the vector of actual inputs of the neurons of the neural network for the k-th step Xk are calculated: X k  X kin  Гk 1 . At the second stage, the transposed adjacency matrix WT is multiplied by the vector Xk: Sk  W T  X k . (11) At the third stage, after the transformation of the vector elements by the functions of activation blocks: S k  s1k , s2 k ,..., sN 1k , sNk  ; (12) where for the k-th step the elements of the signal vector at the outputs of the neurons of the neural network are calculated: Г k   1k ,  2 k ,...,  N 1k ,  Nk  . (13) Then the algorithm returns to the first stage to calculate the values of signals at the outputs of neurons at the (k + 1)-th step of the neural network. The fig. 3 shows graphs of four fully connected recurrent neural networks as test results of the proposed matrix algorithm for calculating a fully connected dynamic recurrent neural network. The arrows in fig. 3 indicate the contacts between the variables (in this case, the possible TV3-117 aircraft engine technical states), while the signal passes only in the direction of the graph arrows. S0 S1 S1 S0 S1 S0 S2 S2 S3 a) b) c) Figure 3: Signal graphs of three fully connected recurrent neural networks 4. Training a neural network with an evolutionary algorithm based on the matrix method The topology of fully connected neural networks is interesting for maximum intelligence when learning. It is well formalized and convenient for the algorithmization of matrix transformation operations of expressions (10) – (13). In matrix form [15], these calculations are convenient because in a fully connected recurrent neural network, the number of synaptic weights is equal to the square of the number of neurons, and in the network any neuron, including the input one, can be used as an output one. For dynamic neural networks, oriented to work in regulation and control systems, an important property is their inertia. In the proposed stepwise algorithm, to calculate the elements of the vector Г at the current k-th step, the values at the outputs of all neurons at the previous (k – 1)-th step are required. Then, if the increment of the signal value at the output of the n-th neuron from step to step [15]:   kn   kn    k 1 n (14) multiply by the accepted time interval between calculation steps t  tk  tk 1 (15) then the operation of integrating the primitive of the simulation dynamic model of the neuron with N inputs will be reproduced in the neuron: N  kn    k 1n  t    k 1 z ; n = 1, 2,…, N – 1, N. (16) z 1 In this case, the step size (the time interval between calculations Δt) will determine the accuracy of the transient simulation. Then the elements of the vector Г will be the stack storage of the associative memory of the dynamic neural network. As activation blocks in the proposed algorithm, we tested simplified primitives of recurrent expressions with unit parameters of one-dimensional simulation models of a linear integral link:  k   k 1  t  sk 1 ; (17) and linear inertial (aperiodic) link  k  1  t    k 1  t  sk 1 . (18) For comparison, the algorithm was tested with the sigmoid function as an activation block: es k 1  e s k 1  k  s k 1  s k 1 . (19) e e It was found that expressions (17) and (18) have almost the same positive effect on the dynamic properties of the trained neural network, and with expression (19) it was not possible to train the neural network for a dynamic mode of operation. Obviously, the sigmoid is applicable only for neural network classifiers. The neural network was trained in the identifier mode using test graphs of the object’s response to a step effect and to harmonic oscillations of the same amplitude at three frequencies. Test signals were generated by a circuit from recurrent expressions of sequential connection of simulation models of aperiodic links (A-links) with parameters: Ta = 10 s and ka = 3 units out / units in:   t  t  yk1  1    y k 11  ka   xk 1   Ta  Ta   t   t  yk 2  1    y k 1 2  ka   yk1 (20)   Ta  Ta   yk 3  1  t   y t   k 13  ka   yk 2  Ta  Ta and a circuit of recurrent expressions of the serial connection of the A-link and the integrating (I-link) with negative feedback t     yk 1  y k 1 2  ka  T  xk 1  y k 1 2  u (21)  y  1   t   t y  ka   k 2  Ta   k 1 2 Ta  In system (21), the parameters of the A-link are Ta = 10 s and ka = 3 units out / units in, and for the I-link the neural network is trained for two options: Ti = 50 s, when the transient process is monotonic, and Ti = 3 s, when the transient is oscillatory. A probabilistic evolutionary algorithm for optimizing multiextremal problems is used as a tool for tuning a neural network. The algorithm consists of a probabilistic Monte Carlo procedure for creating a set (population) of starting points for finding an optimal solution R r1 , r2 ,..., rz 1 , rz  [16], where Z – given number, a genetic procedure for eliminating (removing) 10 % of unpromising points from a set R, an operation of replenishing a set R with new, randomly selected points and their distribution by a regular optimization algorithm over the coordinates of local extrema [17–19]. 5. Description of input data The input data for control and diagnostics of TV3-117 aircraft engine technical state are the results of modeling the transient process in it. The papers [20, 21] describe in detail the general approaches to modeling aviation gas turbine engines together with their automation in the simulator DVIG_OTLADKA. The TV3-117 aircraft engine automatic control system operates as follows: depending on the engine control joystick (ECJ) installation angle, the rotor speed is maintained with a temperature correction at the engine inlet (i.e., the following control laws are implemented (n = f(αECJ, ТН), GT = f(n)), the limiting values of the rotor speed nmax and the gas temperature in front of the compressor turbine TG* are limited. According to [22], the dynamic characteristics of an aircraft engine are determined by the moment of inertia of the compressor rotor. With the help of the structural element «Pump-regulator», the set rotor speed is maintained by changing the fuel consumption in the combustion chamber. For the simulated TV3-117 aircraft engine, the rotor moment of inertia is taken equal to 0.00045 kg ꞏ m2. As described in [23], the dynamic characteristics of various automation elements are unified in the form of coefficients of a second-order nonlinear differential equation: kx  t   TK2 y  t   TD y  t   y  t  ; (22) where k – gain; TK – oscillatory time constant, s; TD – differentiating time constant, s; x(t) – value of the function (for the structural element «Pump-regulator» – the value of the fuel consumption); y(t) – value of the controlled parameter (for the structural element «Pump-regulator» – the reduced or physical rotor speed); y  t  and y  t  – respectively the second and first time derivatives of the value of the controlled parameter. According to [20], the schedule of the transient process in the TV3-117 aircraft engine has the form shown in fig. 4. Figure 4: Fuel consumption graph in the transient process with the settings of the structural element “Pump‐regulator” TD = 0.15 s, TK = 0.25 s. 6. Results and discussion The above-considered evolutionary algorithm trained fully connected recurrent neural networks with two, three, and four neurons, which set the possible operating modes of the TV3-117 aircraft engine. For all options, the same cardinality of the set of initial points was set µ = 800, and the criterion for exiting the solution was the admissible difference between the minimum and maximum values of the optimized function ε = 10-3. The fig. 5 shows the results of the optimization program; adjacency matrices with optimal values of synaptic weights for a neural network: with two (fig. 5, a), with three (fig. 5, b) and with four neurons (fig. 5, c).  0.171 0.188 0.375 0.031     1.089 0.235 0.031   0.127 0.554 0.193 0.623     0.833 1.119 0.745 1.251   0.356 0.158   0.426 0.633 1.225       0.144 0.282 0.399   0.832 1.093     1.613 0.674 0.841 0.113  a) b) c) Figure 5: Adjacency matrices of synaptic weights for recurrent neural networks: a – with two neurons; b – with three neurons; c – with four neurons To train a neural network with four neurons (search for 16 values of synaptic weights for the matrix according to fig. 5, c), the program takes from 10 to 20 seconds. The fig. 6, 7 and 8 show graphs of the results of studying the learning processes of a neural network: with two neurons (a), with three neurons (b) and with four neurons (c). It can be seen that the identification accuracy depends on the number of neurons, and the dynamic error of the neural network depends on the rate of change of the modeled process. The dynamic error of modeling transient processes is no more than 10–4 attained by a neural network with four neurons. a) b) c) Figure 6: Graphs of learning outcomes of dynamic recurrent neural networks for system identification (20): 1 – neural network; 2 – analytical calculation a) b) c) Figure 7: Graphs of learning outcomes of dynamic recurrent neural networks for the identification of system (21) with A‐link parameters Ta = 10 s and ka = 3 units out / units in and I‐link: Ti = 50 s: 1– neural network; 2 – analytical calculation a) b) c) Figure 8: Graphs of learning outcomes of dynamic recurrent neural networks for system identification (21) with A‐link parameters Ta = 10 s and ka = 3 units out / units in and I‐link: Ti = 3 s: 1 – neural network; 2 – analytical calculation The fig. 9 and 10 show graphs of the results of studying the learning processes of dynamic recurrent neural networks with four neurons to identify the transient process (fig. 4) at the output: a – first neuron, b – second neuron, c – third neuron. a) b) c) Figure 9: Graphs of the results of studying the learning processes of a dynamic recurrent network with four neurons to identify the transient process at the output: a – first neuron, b – second neuron, c – third neuron for transient processes in system (21) with parameters Ta = 10 s and ka = 3 units out / units in and I‐link: Ti = 50 s: 1 – neural network; 2 – analytical calculation a) b) c) Figure 10: Graphs of the results of studying the learning processes of a dynamic recurrent network with four neurons to identify the transient process at the output: a – first neuron, b – second neuron, c – third neuron for transient processes in system (21) with parameters Ta = 10 s and ka = 3 units out / units in and I‐link: Ti = 3 s: 1 – neural network; 2 – analytical calculation In the work, a neural network with four neurons is tested for the identification of graphs under various stepwise disturbances. The fig. 11 shows the test results under disturbances: a – 0.5; b – 1.0; c – 2.0 units. Good agreement of the neural network with analytical calculations with a satisfactory error is seen. Figure 11: Graphs of transient processes and dynamic errors of identification of a dynamic system (21) of a recurrent neural network under stepwise disturbances, where “–––” – neural network; “•” – analytical calculation: 1 – original signal; 2, 3 – modified signal (0.5 and 2.0 units, respectively) The fig. 12 shows neural networks and graphs of transient processes obtained at the outputs of all neurons, obtained after training the neural network to identify transient processes at the output of the second neuron: 1 – output from a neuron characterizing the operating mode of the engine Si; 2 – output from neuron 2, characterizing the normal operation of the engine S1 (nominal, I cruising, II cruising) S0 S1 S2 S3 Figure 12: Signal graph of a recurrent neural network trained to identify the transient process at the output of the second neuron, signaling the normal operation of aircraft engine TV3‐117, and graphs of the corresponding transient processes at the outputs of the remaining neurons 7. Conclusions This work is devoted to the implementation of Markov discrete network showing all kinds of changes of TV3-117 aircraft engine technical state as a part of a helicopter power plant during flight tests using neural network technologies, which allows using, along with quantitative mathematical models of the TV3-117 aircraft engine, qualitative and experimental information obtained during flight tests. For modeling and implementation Markov discrete network showing all kinds of changes of TV3- 117 aircraft engine technical state as a part of a helicopter power plant during flight tests as dynamic systems in control algorithms, a dynamic recurrent neural network, encoded by signal graphs, is proposed. To calculate dynamic recurrent neural networks, a step-by-step matrix method is proposed that uses the signal graph adjacency matrix as an associative memory storage, which made it possible to obtain sufficient accuracy in approximating the thermogas-dynamic parameter signal, for example, fuel consumption in the transient process with the settings of the structural element «Pump-regulator». It was found that to identify a third-order dynamic system, it is sufficient to train a fully connected recurrent neural network with four neurons, symbolizing the vector of probabilities of TV3-117 aircraft engine technical states, in which the number of synaptic weights is equal to the square of the number of neurons, and any neuron in the network, including the input one, can be used as an output. Studies have shown that in a neural network trained for the selected (second) output neuron, active signals dynamically connected with the signal on the second neuron are observed at the outputs of its other neurons. Research prospects are the use of the obtained properties of recurrent neural networks as virtual intermediate signals in the control and diagnostics system of aircraft GTEs technical state, including TV3-117, as dynamic systems. The results of this work can be introduced into an intelligent on-board system for control and diagnosing of aircraft GTEs technical state, including TV3-117 [24]. 8. References [1] D. S. Legkonogikh, Using of neural networks technologies in diagnostics Systems of aviation powerplants, Izvestia of Samara Scientific Center of the Russian Academy of Sciences, 2012, vol. 14, no. 4 (2), pp. 639–643. [2] E. L. Ntantis, Diagnostic methods for an aircraft engine performance, Journal of Engineering Science and Technology, 2015, vol. 8, no. 4, pp. 64–72. [3] A. M. Pashayev, D. D. Askerov, C. Ardil, R. A. Sadiqov, P. S. Abdullayev, Condition monitoring system of aircraft gas turbine engine complex, International Journal of Aerospace and Mechanical Engineering, 2007, vol. 1, no. 11, pp. 689–695. [4] S. Vladov, K. Kotliarov, S. Hrybanova, O. Husarova, I. Derevyanko, S. Gvozdik, Neuro- mechanical methods of control and diagnostics of the technical state of aircraft engine TV3-117 in film regions, Visnyk of Kherson National Technical University, 2020, no. 1 (72), part 1, pp. 141–154. [5] S. Vladov, Yu. Shmelov, T. Shmelova, Modeling of the TV3-117 aircraft engine technical state as part of the helicopter power plant in the form of the Markov process of death and reproduction, in: CEUR Workshop Proceedings, 2020, vol. 2740, pp. 400–407. [6] S. I. Vladov, Ya. R. Klimova, O. I. Korniienko, Method for determining the probabilities of Mi- 8MTV helicopter engine state under flight test conditions, in: Aviation and cosmonautics: state, achievements and prospects, 2018, pp. 26–27. [7] J. Bang-Jensen, G. Gutin, Digraphs Theory, Algorithms and Applications, Springer-Verlag GmbH, 2007. [8] M. Needham, A. E. Hodler, Graph Algorithms: Practical Examples in Apache Spark and Neo4j, O'Reilly Media. 2019. [9] M. Zgurovsky, V. Sineglazov, E. Chumachenko, Artificial Intelligence Systems Based on Hybrid Neural Networks: Theory and Applications, Springer, 2021. [10] E. V. Bodyansky, O. G. Rudenko, Artificial neural networks: architectures, training, applications, Teletech, Kharkiv, 2004. [11] S. V. Zhernakov, A. T. Gilmanshin, New algorithms for on-board diagnostics of aircraft gas turbine engine based on fuzzy networks, Bulletin of USATU, 2015, vol. 19, no. 2 (68), pp. 63– 68. [12] S. V. Zhernakov, Control and diagnostics of complex failures of ACS GTE using recurrent neural networks, Aerospace Instrument-Making, 2009, no. 2, pp. 33–49. [13] V. B. Lila, E. V. Puchkov, Methodology of training recurrent artificial neural network with dynamic stack memory, Software & Systems, 2014, no. 4 (108), pp. 132–135. [14] N. I. Selvesyuk, A. S. Ostrovskiy, R. S. Aristov, A. A. Gladkikh, Object-oriented design of a neural network to automate the process of computer architecture determination in the information security problems, Science Bulletin of the NSTU, 2016, vol. 62, no. 1, pp. 133–145. [15] V. R. Sabanin, D. A. Dementyev, I. Yu. Kazmiruk, A. I. Repin, Matrix method for calculating dynamic recurrent artificial neural networks, Young scientist, 2018, no. 13 (199), pp. 1–9. [16] S. P. Li, A guided Monte Carlo method for optimization problems, International Journal of Modern Physics C, 2002, vol. 13, no. 10, pp. 1365–1374. [17] V. R. Sabanin, N. I. Smirnov, A. I. Repin, Modified genetic algorithm for optimization and control problems, Exponenta Pro. Mathematics in applications, 2004, no. 3–4, pp. 78–85. [18] A. A. Andreyev, Use of the genetic algorithms in optimization of non-linear functions, Tomsk State University Journal, 2009, vol. 14, no. 5, pp. 1036–1040. [19] A. N. Koshev, V. V. Salmin, A. A. Generalova, D. S. Bichkov, The development of genetic algorithm with adaptive mutations to determine the global extremum function of n-variables, The Eurasian Scientific Journal, 2016, vol. 8, no. 6, URL: https://naukovedenie.ru/PDF/32TVN616.pdf [20] D. A. Akhmedzyanov, I. A. Krivosheev, A. E. Kishalov, Simulation of the operation of aircraft GTE with elements of control systems, Bulletin of USATU, 2008, vol. 11, no. 2 (29), pp. 3–11. [21] D. A. Akhmedzyanov, A. E. Kishalov, Simulation of an aircraft GTE together with elements of its automation, Youth Bulletin of USATU, 2012, no. 2 (3), pp. 54–60. [22] D. A. Akhmedzyanov, A. E. Kishalov, K. V. Markina, Z. N. Bakirova, Experimental study of static and dynamic characteristics of gas turbine engines on automated laboratory units SR-30 and TJ-100 with various engine control algorithms, Young scientist, 2013, no. 3 (50), pp. 26–34. [23] D. A. Akhmedzyanov, A. E. Kishalov, An integrated approach to joint modeling and research of working processes of aircraft GTE and its automation, Young scientist, 2011, vol. 1, no. 10 (33), pp. 16–21. [24] Y. Shmelov, S. Vladov, Y. Klimova, M. Kirukhina, Expert system for identification of the technical state of the aircraft engine TV3-117 in flight modes, in: System Analysis & Intelligent Computing : IEEE First International Conference on System Analysis & Intelligent Computing (SAIC), 2018, pp. 77–82.