<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-title-group>
        <journal-title>Journal of Computational Physics</journal-title>
      </journal-title-group>
      <issn pub-type="ppub">1613-0073</issn>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.35546/kntu2078-4481.2021.4.7</article-id>
      <title-group>
        <article-title>Classifier of Helicopters Turboshaft Engines Operational Status</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Serhii Vladov</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Ruslan Yakovliev</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Oleksandr Hubachov</string-name>
          <email>oleksandrgubachov@gmail.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Juliia Rud</string-name>
          <email>juliarud25@gmail.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Workshop</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Kremenchuk Flight College of Kharkiv National University of Internal Affairs</institution>
          ,
          <addr-line>Peremohy street 17/6</addr-line>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Kremenchuk</institution>
          ,
          <addr-line>39605</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2022</year>
      </pub-date>
      <volume>206</volume>
      <issue>4</issue>
      <fpage>52</fpage>
      <lpage>63</lpage>
      <abstract>
        <p>The work is devoted to the development of an on-board neural network classifier of helicopters turboshaft engines operational status. Proposed neural network classifier was developed based on an ensemble of neural networks, which consists of radial basis networks (RBF), perceptron, Kohonen network and hybrid neural network. The main task of the developed on-board neural network classifier is to determine defects in helicopters turboshaft engines units (air inlet section, compressor, combustion chamber, compressor turbine, free turbine, exhaust unit). It has been proven that to obtain the best result, it is advisable to use the following algorithms for training neural networks: a modified gradient algorithm for training an RBF network, a backpropagation method for training a multilayer perceptron, a competition method between neurons for training a Kohonen neural network, a hybrid algorithm - for training a hybrid neural network. The results of testing the developed on-board neural network classifier showed the ability to determine a compressor defect, a compressor turbine defect, and simultaneous compressor and compressor turbine defects. The effectiveness of the developed on-board neural network classifier for recognizing defects in helicopters turboshaft engines has been proven. A comparative assessment of the effectiveness of the developed on-board neural network classifier and existing methods for parametric diagnostics of the operational status of complex dynamic objects was carried out. The results of the studies showed that the developed on-board neural network classifier can identify defects in helicopters turboshaft engines components with an accuracy of up to 99.8%. Helicopters turboshaft engines, neural network classifier, RBF network, multilayer perceptron, Kohonen neural network, hybrid neural network, thermogas-dynamic parameters, training, 0000-0002-0328-5895 (J. Rud) Proceedings ITTAP'2023: 3rd International Workshop on Information Technologies: Theoretical and Applied Problems, November 22-24, Proceedings</p>
      </abstract>
      <kwd-group>
        <kwd>operational status</kwd>
        <kwd>diagnostics</kwd>
        <kwd>defects</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>Helicopters</p>
    </sec>
    <sec id="sec-2">
      <title>1. Introduction</title>
      <p>During the operation of aviation equipment such as aircraft and helicopters, a primary objective is
to diagnose the parameters of their aircraft engines. The total count of monitored (diagnosed)
parameters can extend to 500 or more. The existing methods and techniques for diagnosing gas turbine
engines (GTE) demand substantial enhancements, particularly as new generations of aircraft GTE
necessitate advanced intelligent computer diagnostic technologies founded on the principles of expert
systems (ES), neural networks (NN), fuzzy logic (FL), and genetic algorithms (GA). These technologies
should be capable of incorporating the accumulated experience from prior work in this domain and
devising (generalizing) innovative methods and techniques for further exploration. This imperative
applies equally to aircraft gas turbine engines with a free turbine, known as turboshaft engines (TE),
ORCID:
EMAIL:
(S.
️©</p>
      <p>2020 Copyright for this paper by its authors.
CEUR</p>
      <p>
        ceur-ws.org
which form integral components of helicopter power plants. Among the intricate tasks that significantly
enhance the efficiency of GTE diagnostics and elements of automatic control systems (ACS), there is
a need to address several issues aimed at overcoming obstacles in identifying the operational status of
GTE. These challenges are interconnected [
        <xref ref-type="bibr" rid="ref1 ref2">1, 2</xref>
        ]:
      </p>
      <p>– in the event of primary information sensor malfunctions, there is a risk of generating false alarms
within the gas turbine engine (GTE) control system, leading to a considerable decline in failure
identification reliability. To ensure the accurate operation of the GTE operational status monitoring
system, it becomes essential to distinguish (classify) deviations stemming from alterations in power
plant characteristics from deviations in measured parameters linked to sensor malfunctions. In essence,
this involves concurrently identifying the engine status, the parameters of its gas flow duct, and the
measurement system, all while simultaneously identifying the program regulation;
– complications arise in distinguishing failures of engine components and sensor malfunctions
during engine failures, including its subsystems. This challenge is particularly evident when facing
minor deviations in gas-dynamic parameters, such as those occurring when individual turbine blades
burn out, which are comparable to random errors in the measuring channels;</p>
      <p>– challenges arise in the automatic acquisition and extraction of essential, reliable information
during each flight, encompassing both steady and transient operating modes. This involves obtaining
independent measurements following each engine transition to a new steady-state mode and ensuring a
distinct separation between transient and steady modes. Addressing these challenges is imperative for
enhancing the reliability of monitoring and diagnosing the operational status of the engine and the
elements of the ACS during flight, especially when dealing with slight deviations from the anticipated
standard of the measured parameters.</p>
      <p>
        In such circumstances, the application of neural network technologies holds significant promise. A
review of research in the domain of gas turbine engine operational status diagnostics utilizing neural
networks [
        <xref ref-type="bibr" rid="ref3 ref4">3, 4</xref>
        ] indicates that ongoing work is underway. However, due to various reasons such as
secrecy and the narrow specialization of the tasks at hand, most publications lack engineering
methodologies, as well as theoretical and practical guidance for addressing such issues.
      </p>
      <p>Therefore, the objective of this study is to formulate methods and techniques for the comprehensive
diagnostics of helicopters TE operational status during flight modes utilizing neural network
technologies.</p>
    </sec>
    <sec id="sec-3">
      <title>2. Related works</title>
      <p>
        A novel and promising domain in the realm of automatic control for intricate dynamic systems,
operational status diagnostics, and predictive tasks involves the utilization of intelligent control systems
founded on artificial neural networks [
        <xref ref-type="bibr" rid="ref5 ref6">5, 6</xref>
        ].
      </p>
      <p>
        However, prevailing approaches to employing intelligent diagnostic methods are constrained by the
specificity of the tasks, the underdeveloped theory pertaining to the use of neural networks in gas turbine
engine (GTE) diagnostics, the absence of universal and formalized approaches, and the inherent
imperfections within neural network methods themselves [
        <xref ref-type="bibr" rid="ref7 ref8">7, 8</xref>
        ].
      </p>
      <p>
        Investigations into the creation of automated systems for diagnosing the operational status of
complex dynamic objects [
        <xref ref-type="bibr" rid="ref10 ref9">9, 10</xref>
        ], including aircraft GTE [
        <xref ref-type="bibr" rid="ref11 ref12">11, 12</xref>
        ], reveal the inadequacy of relying
solely on one of the known diagnostic methods. This is because none of these methods is universally
applicable and entirely reliable. Consequently, monitoring and diagnostics systems built on the
foundation of a single classifier will likely fall short of meeting the escalating requirements for gas
turbine engine diagnosis. To enhance the efficiency of on-board technologies for GTE operational status
diagnostics, several directions are identified. The primary focus should be on intellectualizing
information processing processes through the application of neural network methods [
        <xref ref-type="bibr" rid="ref13 ref14">13, 14</xref>
        ], which
can enhance the quality of on-board algorithms for GTE operational status diagnostics.
      </p>
      <p>In this context, the development of neural network methods for diagnosing the operational status
and identifying potential defects in helicopters' turboshaft engine units during flight modes remains
pertinent.</p>
    </sec>
    <sec id="sec-4">
      <title>3. Proposed technique 3.1.</title>
    </sec>
    <sec id="sec-5">
      <title>Problem statement</title>
      <p>Let an N-dimensional feature space be given, each point of which can be represented by an
Ndimensional vector X = X1, …, XN. Let's divide this space into Q regions, which correspond to one class
or another. Let the training set {X, В} = (X1, B1), (X2, B2), ..., (XN, BN) be given, where Xi – point in the
feature space, Bi – label of the class to which this point belongs. The job of the classifier is to indicate
for each new point X that is not included in the training sample, under conditions of partial or complete
uncertainty, which class this point belongs to, using the training sample {X, В} for this. Let the state of
a complex dynamic object (CDO) (helicopter TE) X at each discrete time t be described by an
Ndimensional vector X t = ( X1t , X 2t ,..., X Nt ) of variables satisfying N equations X t+1 = F ( X t ,Qt ) , k = 1,
2, …, M, where Qt = (q1t ,q2t ,...,qmt ) – reference (defect-free) vector the state of the CDO (helicopter
aircraft TE). Changes in CDO operational status at any time can be described by the equation
Y t = H  X t , where Y t = ( y1t , y2t ,..., yRt ) – operational status vector of real output parameters, H –
transformation matrix. It is required to determine the diagnostic state vector of CDO (helicopter TE),
minimizing the root-mean-square error between the reference (desired) Y and the real YRt outputs.</p>
      <p>The organization of complex diagnostics for helicopters TE is notably intricate due to the need for
objective and comprehensive information about their operational status. This requires incorporating a
substantial number of diverse physical quantities (parameters) into the diagnostic procedure to capture the
behavior of various subsystems. The transition of a gas turbine engine from one state to another is marked
by discernible changes in the controlled and diagnosed parameters. Fig. 1 provides a geometric
representation of complex diagnostics for helicopters' turboshaft engines, which can be articulated as
follows: based on a limited number of measurements of the object being diagnosed, an optimal decision
must be made regarding its classification into one of several classes, specifically: S1 – the area of
serviceable states; S2 – the area of critical states; S3 – the area of faulty states; P1, P2, P3 – classes (areas)
of uncertain states.</p>
      <p>p2</p>
      <p>S2
x3
S1
p3</p>
      <p>S3
xN
p1</p>
      <p>x2</p>
      <p>
        Geometrically, the operational status of helicopters TE can be visualized as an N-dimensional vector
(XN) (fig. 1), where the spatial coordinates represent N input parameters of the engine (X1, X2, …, XN).
The position of this state vector in space corresponds to a specific level of engine performance, and the
establishment of standards involves creating separating hypersurfaces within this space. These
hypersurfaces act as boundaries between different classes, determined by a decision rule that dictates
their construction. Consequently, decision-making revolves around assigning the diagnosed object to a
specific class. The input parameters X1, X2, …, XN refer to the thermogas-dynamic parameters of the
working process of helicopters TE (Table 1). These parameters are either recorded on board the
helicopter or computed using a comprehensive mathematical model of an aircraft engine with a free
turbine [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]. The transition from the physical parameters of the engine to the specified values (and vice
versa) follows a developed methodology [
        <xref ref-type="bibr" rid="ref16">16, 17</xref>
        ].
      </p>
      <p>gas generator rotor r.p.m., nTC
air flow through the compressor, Gair
air pressure behind the compressor, PC
air temperature behind the compressor, TC
total gas pressure behind the combustion chamber, PG
gas temperature in front of the compressor turbine, TG</p>
      <p>fuel consumption, GT
total gas pressure behind the compressor turbine, PTC
gas temperature behind the compressor turbine, TTC
total gas pressure behind the free turbine, PFT</p>
      <p>free turbine rotor speed, nFT
total gas pressure behind the exhaust unit, POUT
gas temperature behind the exhaust unit, TOUT</p>
      <p>Determination
calculated analytically</p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
calculated analytically
      </p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
registered on board
      </p>
      <p>the helicopter
calculated analytically</p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
calculated analytically
      </p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
calculated analytically
      </p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
calculated analytically
      </p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
registered on board
      </p>
      <p>the helicopter
calculated analytically</p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
calculated analytically
      </p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
calculated analytically
      </p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
calculated analytically
      </p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
calculated analytically
      </p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
registered on board
      </p>
      <p>the helicopter
calculated analytically</p>
      <p>
        according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
calculated analytically
according to [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]
Free turbine
      </p>
      <p>gas temperature behind the free turbine, TFT</p>
      <p>The output diagnostic parameters of helicopters TE are: degree of increase in the total pressure in
the compressor  C* , compressor efficiency ηC, mechanical compressor efficiency ηMC, recovery factor
of the total gas pressure in the combustion chamber σCC, combustion chamber cross-sectional area FCC,
compressor turbine efficiency ηTC, compressor turbine operation ATC, degree of reduction of the total
gas pressure in the compressor turbine  T*C , total pressure reduction ratio in the free turbine  *FT , power
efficiency of a free turbine ηΣFT, total pressure reduction ratio in the free turbine and exhaust unit  *FT
, total gas pressure recovery factor in the exhaust unit σEU.
3.2.</p>
    </sec>
    <sec id="sec-6">
      <title>Neural network classifier development</title>
      <p>In addressing the intricate challenges of diagnosing helicopters' turboshaft engines, hybrid
ensembles of neural networks [18] can serve effectively as dynamic repositories of expert knowledge.
In comparison with traditional neural networks, these ensembles offer additional practical advantages,</p>
      <p>To address the issue at hand, a neural network classifier was devised (fig. 2), utilizing a composite
ensemble of neural networks. This ensemble incorporates radial basis networks (RBF), perceptron,
Kohonen network, and a hybrid neural network.
including: decomposition of complex dynamic objects (CDO) into simpler entities or subsystems;
enhanced adaptability to changing external conditions, positioning them within the class of adaptive
and self-adjusting systems; optimization of the neural ensemble's structure to cater to specific diagnostic
tasks; superior speed and accuracy compared to classical fully connected networks; improved
approximation of piecewise continuous functions by the neural ensemble [19, 20].</p>
      <p>To adapt the recognition of defects by a neural ensemble, as part of the training sample, we
distinguish five generalized classes of engine status (table 2): S0 – serviceable (reference) status
corresponding to the vector R = [0; 0; 0]; S1 – compressor defect corresponding to the vector R = [0; 1;
0]; S2 – combustion chamber defect corresponding to the vector R = [0; 1; 1]; S3 – compressor turbine
defect corresponding to the vector R = [1; 0; 0]; S4 – free turbine defect corresponding to the vector R
= [1; 1; 0]; S5 – exhaust unit defect corresponding to the vector R = [1; 0; 1].
... ... ... ...</p>
      <p>S</p>
      <p>Opting for the radial basis function (RBF) architecture as the recognition neural network, as opposed
to the perceptron neural network, is more advantageous. This preference stems from the fact that the
determination of weight coefficients in the RBF network is accomplished more swiftly and accurately
compared to adjusting the parameters of the perceptron. This efficiency arises because the utilization of
gradient methods for parameter adjustment in the perceptron often leads to the attainment of local minima.
3.3.</p>
    </sec>
    <sec id="sec-7">
      <title>RBF network algorithm training</title>
      <p>
        The structure of the RBF network entails a two-layer design, where the first layer executes a
predetermined non-linear transformation without engaging in parameter tuning. This process maps the
input space to a new space. In this context, considering the 16 input parameters of the helicopters TE
(table 1), an optimal configuration, in terms of parameter decomposition, would involve six RBF
architecture neural networks. These networks would align with the number of parameters at the input
(state vector) based on the engine node (ranging from 2 to 4, depending on the node), and three output
parameters in line with the binary classification of statuses (table 2). The training algorithm for the RBF
neural network employs the modified gradient training of RBF networks, as developed in [
        <xref ref-type="bibr" rid="ref16">16</xref>
        ], with the
block diagram illustrated in fig. 3, where n – number of parts in the first layer; x1, x2, ..., xn – input signals;
m – number of elements in the second layer; ci1, ci2, ..., cin – coordinates of the center of the i-th element;
σi – width of the radial function of the i-th element; θi – output signal of the i-th element; wi – weight
coefficient of the initial connection of the i-th element; y – output signal of the RBF neural network.
      </p>
      <sec id="sec-7-1">
        <title>Start 1</title>
      </sec>
      <sec id="sec-7-2">
        <title>Receiving input signals</title>
        <p>x1, x2, xn
and output signal y
2
Calling the procedure for
adding and removing
elements
i = 1...m, m – number of
elements</p>
      </sec>
      <sec id="sec-7-3">
        <title>Calculation output value</title>
        <p>of і-th element θі</p>
        <p>θі &gt; θinput
6</p>
        <p>Yes
Calculation Δwi, Δσi, ρi,</p>
        <p>wi = wi + Δwi
ρi &gt; ρгр і Δσi &gt; 0</p>
        <p>Yes
Δσi = Δσi kcond
σi = σi + Δσi
10
j = 1...n, n – number of</p>
        <p>input signals</p>
        <p>Calculation Δcij
12 ρi &gt; ρlim і
No sign(Δcij) = sign(cdj – cij)</p>
        <p>Yes
Δcij = Δcij kcond
cij = cij + Δcij
3
4
8
9
No 5</p>
        <p>No 7</p>
        <p>
          The Gaussian function governs the output signal of each component within the radial basis function
neural network [
          <xref ref-type="bibr" rid="ref16">16, 21</xref>
          ]:
        </p>
        <p> i = e . (1)</p>
        <p>The output signal of the radial basis function neural network is computed by aggregating the signals
of its elements through a weighted sum:
n
(xj −cij )2
− j=1</p>
        <p>2 i2
m
y =  wi  i . (2)</p>
        <p>i=1</p>
        <p>A gradient-based algorithm is employed to train the radial basis function (RBF) neural network by
minimizing the network error objective function. This algorithm computes the adjustments for each
element, including changes in the weight factor (wi), the element width (Δσi), and the element center
coordinates (cij).</p>
        <p>Through experiments, certain drawbacks of the conventional gradient algorithm for RBF neural
network training were identified:</p>
        <p>1. The RBF neural network training algorithm lacks specific guidelines for the initial assignment of
network elements and their parameters, as well as rules for adjusting the element count during training.
The uniform distribution of elements across the working area may not always be optimal, and there
could be instances where the initially specified number of elements proves insufficient for achieving
the required training quality.</p>
        <p>2. Throughout the training process, the parameters of all network elements undergo changes,
resulting in an escalation of computational costs as the element count increases.</p>
        <p>3. The RBF network faces challenges in reaching a stable state during training when elements with
closely positioned center coordinates (cij) and radial function widths (σi) exist. The occurrence of such
situations is highly dependent on the chosen number of elements and their initial parameters. The
degradation of training quality stems from the assumption in the gradient algorithm that the output value
of the RBF neural network at any given point in the working area is primarily influenced by only one
element. When multiple elements are present in a specific region of the working area, adjusting their
parameters according to the gradient algorithm does not consistently lead to a reduction in training error.</p>
        <p>To identify instances where the parameters of certain elements converge, the notion of the mutual
intersection coefficient of elements has been introduced. To compute this coefficient for a specific
element within the RBF network, it is essential to identify a second element whose center is in closer
proximity to the center of the analyzed element. The mutual intersection coefficient is determined by
summing the initial value of the current element at the center of the second element and the initial value
of the second element at the center of the current element:
n n
(cij −cdj )2 (cij −cdj )2
− j=1 − j=1
 i = e + e ; (3)
where i – number of the element for which the value of the coefficient of mutual intersection is
calculated; d – element number, the center of which is located closer to the center of the element with
number i, which is determined according to the expression:
2 i2</p>
        <p>2 d2
n 2
d = arg min (cij − ckj ) . (4)</p>
        <p>k j=1</p>
        <p>The coefficient of mutual intersection resides within the range (0; 2). It achieves its maximum value
when the centers of the examined elements coincide. As the coefficient of mutual intersection surpasses
1.95, it becomes imperative to cap the maximum value at 1.95 to attain optimal training quality for the
RBF network.</p>
        <p>
          To address the deficiencies of the conventional gradient algorithm in training the RBF network for
helicopter TE identification tasks, this study employs a modified gradient algorithm [
          <xref ref-type="bibr" rid="ref16">16, 22</xref>
          ], as
illustrated in Fig. 3. Blocks absent in the classic algorithm are denoted by asterisks. The key distinctions
from the classical algorithm are outlined as follows:
        </p>
        <p>1. Introducing regulations for altering the RBF network structure during training (block 2). At the
onset of neural network training, the RBF network is devoid of elements. New elements are introduced
as necessary, and unused elements are eliminated.</p>
        <p>2. A reduction in the computational costs associated with each training cycle is achieved by
modifying the parameters not for all elements, as in the classic algorithm, but exclusively for elements
whose initial value at the specified point surpasses the value of θzm (blocks 4 and 5).</p>
        <p>3. The likelihood of a scenario where the parameters of certain elements closely align is mitigated.
To achieve this, the calculated values Δcij and Δσi are diminished if the coefficient of mutual intersection
among elements surpasses the threshold value ρgr, set at 1.95 (blocks 7, 8, 12, 13).</p>
        <p>Altering the configuration of the RBF network through the addition or removal of elements results
in a modification of the RBF network's output value solely in the proximity of the added or removed
element’s center. This effect is localized and does not impact the entire working area, unlike the
alteration of a multilayer perceptron's structure. Consequently, the addition and removal of elements in
the RBF network can be executed during the training process without the necessity of restarting the
training process from the beginning.
3.4.</p>
      </sec>
    </sec>
    <sec id="sec-8">
      <title>Perceptron neural network algorithm training</title>
      <p>Within the hybrid ensemble, the perceptron neural network serves as a focal field, amalgamating the
outputs of six RBF-architecture neural networks, each with three outputs. In the experimental phase of
this neural network, following the contrasting process, its architecture manifested as a six-layer
perceptron: the initial input layer contained 18 neurons, succeeded by a second hidden layer with 15
neurons, a third hidden layer with 12 neurons, a fourth hidden layer with 9 neurons, a fifth hidden layer
with 6 neurons, and a final output layer comprising three neurons. The perceptron neural network
contributes to the enhancement of defect recognition quality by "fine-tuning" its weight coefficients.
General equations describing the operation of the perceptron: for the input layer (k = 1) – U1 = X (input
n2
vector); for the first hidden layer (k = 2) – U 2j = f  wi1j Ui1 ; … ; for the sixth output layer (k = 6) –
i=1
n6
Yj = U 6j = f  wi6j Ui6 . Table 3 shows a comparative analysis of the perceptron training results, on the
i=1
basis of which the error backpropagation algorithm is selected as the perceptron training algorithm.</p>
    </sec>
    <sec id="sec-9">
      <title>Kohonen neural network algorithm training</title>
      <p>The Kohonen neural network, functioning as a classifier with three inputs and six outputs, exhibits
a high level of precision in categorizing (identifying) the operational status of helicopter TE. This
includes accommodating the partial or complete uncertainty of its parameters. Simultaneously, the
Kohonen neural network serves the purpose of initial sorting and clustering of incoming values,
contributing to the structuring of the initial data for the hybrid neural network. As per [23], the training
foundation for the Kohonen neural network involves a competitive process among neurons. In this
scenario, the Kohonen network features three inputs and six outputs (R1...R6, corresponding to the
number of generalized state classes). The weight coefficients of synaptic connections for each i-th
neuron in the output layer of the Kohonen neural network collectively form a vector
wi = ( wi1 ; wi2 ;...; wi10 )T at i = 1, 2, ..., n. When the Kohonen neural network is activated by the input
vector ΔY, the neuron whose weights are the least different from the corresponding components of the
input vector wins the competition, i.e., for the winning neuron wp, the relation is holds:
d(ΔY, wp) = mind(ΔY, wi), 1 ≤ i ≤ n.</p>
      <p>In equation (5), d(ΔY, wi) represents the distance (as per the chosen metric) between the vectors ΔY
and w = (w1, w2,…, wn)T, where n – number of outputs in the neural network’s output layer (in this
instance, n = 6). The victorious neuron establishes a topological neighborhood Sp(k) around itself,
characterized by a specific energy. All neurons situated within this neighborhood undergo adaptation,
wherein their vectors of weight coefficients undergo changes in the direction of the vector ΔY, following
a prescribed rule:</p>
      <p>wi (k +1) = wi (k ) +i (k )(Y − wi (k ));
where ηi(k) – communication coefficient of the first neuron from the neighborhood Sp(k) at the k-th
moment of time. The training coefficient diminishes as the separation between the i-th neuron and the
victor expands, and the weights of neurons beyond the confines of the Sp(k) neighborhood remain
unchanged.</p>
      <p>The objective of training the Kohonen neural network through neuron competition is to arrange the
neurons (determine the values of their weight coefficients) in a manner that minimizes the expected
distortion value. This distortion is gauged by the approximation error of the input vector ΔY relative to
the weight coefficients of the winning neuron. In the context of L input vectors (ΔY)j, where j = 1, 2, ...,
L, and utilizing the Euclidean metric, this error can be formulated as:</p>
      <p>1 L 2
E =  (Yi ) j − wp ( j ) ; (7)</p>
      <p>L j=1
where wp(j) – weighting coefficients of the winning neuron when the vector network is presented (ΔY)j.</p>
      <p>The results of the training process of the Kohonen neural network (after 700...800 training cycles) are
presented in table 4.
(5)
(6)</p>
    </sec>
    <sec id="sec-10">
      <title>Hybrid neural network algorithm training</title>
      <p>The hybrid neural network enables the assessment of the membership level of the input set of
indicator values to a predetermined class, indicative of the operational status of helicopter TE.</p>
      <p>As previously mentioned, the conclusive step in diagnostics involves determining the type of failure
in helicopter TE based on the analysis of the numerical vector R. Consequently, the graphical
representation of the task (see Fig. 1) undergoes modification to the format depicted in Fig. 4. In this
representation, the vertices of the cube align with the centers of clusters (reference engine states) as
outlined in table 2, that is, S0 – center of the cluster (precedent) corresponding to the serviceable
(reference engine status); S1 – center of the cluster corresponding to a defect in the compressor; S2 –
cluster center corresponding to a defect in the combustion chamber; S3 – center of the cluster
corresponding to a defect in the compressor turbine; S4 – cluster center corresponding to a defect in a
free turbine; S5 – cluster center corresponding to a defect in the exhaust unit. The actual engine status
vector S can take on a value at any point inside the given cube S = (R1, R2, R3)T, 0 ≤ Ri ≤ 1.</p>
      <p>S2
1</p>
      <p>S1</p>
      <p>R3</p>
      <p>The assessment of the operational status of helicopter TE follows the "nearest neighbor" rule,
wherein the engine is categorized into the class that includes its closest neighbor (or the majority of its
closest neighbors). The decision rule, governing the decision-making process (diagnosis), is formulated
as follows:</p>
      <p>S → Sp, if d(S, Sp) → min;
(8)
where d – distance to the center of the nearest (p-th) cluster (precedent), which is calculated, for
example, using the Euclidean metric.</p>
      <p>
        The hybrid training algorithm devised for the hybrid neural network relies on the backpropagation
algorithm as its foundation. Over a specified number of epochs, the neural network undergoes training
employing a modified backpropagation algorithm. Simultaneously, the residual for neurons in the
output layer is computed in a manner consistent with the backpropagation approach. The inconsistency
of the hidden layers is derived from the inconsistencies across all variations of the preceding layer, as
described in [
        <xref ref-type="bibr" rid="ref16">16, 24</xref>
        ]:
 in−1 = f ( Sin )    wln,k ,i  li,k . (9)
l k
      </p>
      <p>Utilizing the residual outlined, akin to the backpropagation algorithm, the weights are incremented
in the direction opposite to the gradient. However, in this instance, the gradient takes the form of a
matrix rather than a vector. Consequently, the adjustments to the weights are made to minimize errors
across all directions (fig. 5). The objective is to take steps that effectively diminish the disparities
between the neural network outputs and the target values at various points.
where x – offset increment.</p>
      <p>When crossing over neurons during training steps using a genetic algorithm, the function is used as

a fitness function f ( x) = 1 − , where ε – accumulated error of the neuron of the output layer (in this
m
case, after the training step using the genetic algorithm, it is reset and takes on a zero value), m = o ∙ l,
where l – number of training sets, o – number of training epochs.</p>
      <p>Block diagram of the developed hybrid training algorithm for hybrid neural networks is shown in
fig. 6, information of the blocks of which is given in table 5.</p>
    </sec>
    <sec id="sec-11">
      <title>4. Experiment</title>
      <p>The examination and initial processing of the input data were conducted by the authors group and
comprehensively detailed in [25, 26]. The input parameters for the mathematical model of helicopter
TE encompass atmospheric values (h – flight altitude, TN – temperature, PN – pressure, ρ – air density).
The parameters recorded aboard the helicopter (nTC – gas generator rotor r.p.m., nFT – free turbine rotor
speed, TG – gas temperature in front of the compressor turbine) were adjusted to absolute values
according to the gas-dynamic similarity theory developed by Professor Valery Avgustinovich (refer to
Number
1
2
3
4
5
6
7
8
9
10
…
256
table 6). In this study, we posit the constancy of atmospheric parameters (h – flight altitude, TN –
temperature, PN – pressure, ρ – air density) [25, 26].</p>
      <p>The assessment of the homogeneity of the training and test samples is a crucial consideration. In this
regard, the Fisher-Pearson criterion χ2 with degrees of freedom r – k –1 is employed [25, 26]:
 2 = min =r1 i  mi n−pni(pi () ) ; (12)
where θ – maximum likelihood estimate found from the frequencies m1, …, mr; n – number of elements
in the sample; pi(θ) – signifies the probabilities of elementary outcomes up to a certain indeterminate
k-dimensional parameter θ.</p>
      <p>The conclusive step in statistical data processing involves their normalization, a procedure that can
be carried out using the following expression:
yi =</p>
      <p>yi − yi min ;
yi max − yi min
where y i – dimensionless quantity in the range [0; 1]; yimin and yimax – minimum and maximum values
of the yi variable.</p>
      <p>The previously mentioned χ2 statistics, under the stated assumptions, enables the testing of the
hypothesis regarding the representativeness of sample variances and the covariance of factors within
the statistical model. The realm of accepting the hypothesis is denoted by  2   n−m, , with α
representing the significance level of the criterion. The outcomes of the computations based on equation
(12) are provided in table 7 [25, 26].</p>
      <p>To verify the representativeness of the training and test samples, an initial data cluster analysis was
conducted (refer to table 7), resulting in the identification of eight classes (fig. 7, a). Subsequent to the
randomization procedure, the effective training (control) and test samples were chosen in a 2:1 ratio,
specifically 67 % and 33 %, respectively. The clustering process applied to the training (fig. 7, b) and
test samples reveals that, akin to the original sample, each of them comprises eight classes. The
distances between the clusters align nearly perfectly across all examined samples, indicating the
representativeness of both the training and test samples [25, 26].</p>
      <p>To form the training and test subsets in the work, cross-validation was used [27] to estimate the
values of the parameters of TV3-117 engine, the results obtained in [27] of which are shown in fig. 8.</p>
    </sec>
    <sec id="sec-12">
      <title>5. Results</title>
      <p>The outcomes of assessing the unknown condition of helicopter TE, demonstrated through the
TV3117 aircraft engine (integral to the power plant of the Mi-8MTV helicopter), using the developed neural
network classifier, are outlined in table 6. Similar to the earlier discussed examples, the performance of
the developed neural network classifier was examined on a test sample (the third and fourth rows of
table 6) in noise-free settings (the first and third rows of table 6), as well as under the influence of the
additive component of white noise (σ = 0.01; M = 0), considering a 1% and 3% alteration in compressor
efficiency (second and fourth rows of table 8). The final line (table 8) corresponds to a dual defect,
involving simultaneous changes in the efficiency of both the compressor and compressor turbine.
operational status deviation (defect size) of the compressor unit; δx3 – operational status deviation
(defect size) of the combustion chamber unit; δx4 – operational status deviation (defect size) of the
compressor turbine unit; δx5 – operational status deviation (defect size) of the free turbine unit; δx6 –
operational status deviation (defect size) of the exhaust unit (see table 1).</p>
      <p>
        Fig. 9 shows the results of the efficiency of the neural network classifier for recognizing defects in
helicopter TE (using the example of the TV3-117 engine): GT – fuel consumption; TC – air temperature
behind the compressor (calculated using [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]); TG – gas temperature in front of the compressor turbine;
nTC, nFT – gas generator rotor r.p.m. and free turbine rotors, respectively; TTC, TFT – gas temperature
behind the compressor turbine and free turbine, respectively (calculated using [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]).
      </p>
      <p>Wherein:
fig. 9, a: curve 1 – GT; curve 2 – TC; curve 3 – TG;
fig. 9, b: curve 1 – TC, TG; curve 2 – GT, TC; curve 3 – GT, TG;
fig. 9, c: curve 1 – GT, TC, TG;
fig. 9, d: curve 1 – TTC; curve 2 – GT; curve 3 – nFT; curve 4 – TFT; curve 5 – nTC;
fig. 9, e: curve 1 – GT, nTC, nFT, TTC; curve 2 – GT, nTC, TTC, TFT; curve 3 – GT, nFT, TTC, TFT; curve 4
– GT, nTC, nFT, TFT; curve 5 – nTC, nFT, TTC, TFT;
fig. 9, f: curve 1 – GT, nTC, nFT, TTC, TFT.
d e
Figure 9: Results of the effectiveness of the neural network classifier for recognizing defects in
helicopter turboshaft engines</p>
      <p>It was found (fig. 7) that as the depth of node defects δxi increases in the training interval 0 ... 10 %,
the recognition efficiency first increases (the value of ∆Σ decreases), reaching its maximum (the
minimum value of ∆Σ), and then decreases. This feature of the dependence ∆Σ = f(δxi) can be used for
practical purposes, for example, the boundaries of the training interval should be selected from the
condition of obtaining the greatest efficiency in recognizing the defective state of the engine (for an
engine compressor unit this may be the boundary value of the reduced efficiency of the compressor at
which it is required engine flushing, etc.).</p>
    </sec>
    <sec id="sec-13">
      <title>6. Discussions</title>
      <p>The results of a comparative analysis of the accuracy and error in diagnosing defects in helicopters
TE units (on the example of defects in the compressor and compressor turbine) are given in table 9.</p>
      <p>Currently, several methods are known for parametric diagnostics of the GTE operational status, which
can be divided into methods A, B, C, D and E [29]: A – method of diagnostic matrices; B – method
based on solving a system of normal equations; C – method based on nonlinear optimization of a criterion
characterizing the state of the engine; D – method of adjustment using a square objective function; E –
method of adjustment using a modular objective function. The results of comparing the effectiveness of
the neural network approach with methods A, B, C, D and E show (fig. 10) that the advantage of the
developed neural network classifier over other methods increases as information about the controlled
engine parameters decreases. Wherein: fig. 10, a: diagnostics by parameter TC; fig. 10, b: diagnostics by
parameter TG; fig. 10, c: diagnostics by parameter GT; fig. 10, d: diagnostics by parameter GT and TG; fig.
10, e: diagnostics by parameter GT and TC; fig. 11, f: diagnostics by parameter TC and TG.
d e f
Figure 10: Comparative assessment of the effectiveness of methods for diagnosing the condition using
the developed neural network classifier and methods A, B, C, D and E</p>
      <p>Since the effectiveness of the considered methods for diagnostics the operational status of helicopter
TE varies in different situations, it is obvious that the combined method is optimal.</p>
    </sec>
    <sec id="sec-14">
      <title>7. Conclusions</title>
      <p>The enhancement of the neural network diagnostics method involves the utilization of an ensemble
comprising six radial-basis neural networks, a perceptron, a Kohonen neural network, and a hybrid
neural network. This approach, applied to experimental data recorded during helicopter operation or
data obtained through a mathematical model, proves effective in addressing the challenge of diagnosing
the operational status of helicopter engines during flight.</p>
      <p>Unlike diagnostic methods reliant on calculating the thermogas-dynamic parameters of helicopters
turboshaft engines using nonlinear element-by-element engine models, the neural network diagnostics
method is refined by training the neural network based on a small training sample. The quality of the
resulting neural network model is then evaluated on a meticulously organized test sample.</p>
      <p>The hybrid neural network training algorithm is enhanced by combining the error backpropagation
algorithm with a genetic algorithm. This involves altering the range of initialization for neuron weights
and biases in the input, hidden, and output layers. This refinement enables the training of the hybrid
neural network to determine the degree of membership of indicator values at their inputs to a specified
class characterizing the operational status of helicopters turboshaft engines.</p>
      <p>The merit of incorporating the Kohonen network in a neural network classifier (part of a neural
networks ensemble) for diagnosing helicopters turboshaft engines operational status in flight modes lies
in its capacity for automatic classification (clustering) without expert instructions. This involves
processing a training sample composed of real or calculated (reference) data across various engine
operating modes.</p>
      <p>For decision-making regarding the location and nature of defects in helicopters turboshaft engines,
the estimation of neuroclassifier outputs using the nearest neighbor rule with precedents (codes of
reference condition) can be employed. The metric value (distance to the nearest precedent) provides
insight into defect intensity or multiplicity (i.e., the number of simultaneously manifesting defects).</p>
      <p>The developed neural network classifier, initially designed for helicopters turboshaft engines, has
the potential for diagnosing other engine types, such as turbojet and turboprop engines used in aircraft
power plants. However, when extending the use of the classifier to other aircraft engines, considerations
must be made for their structural blocks (e.g., fan, low-pressure compressor, high-pressure compressor,
etc.) and, consequently, the applied thermogas-dynamic parameters.</p>
    </sec>
    <sec id="sec-15">
      <title>8. References</title>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>E. L.</given-names>
            <surname>Ntantis</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Botsaris</surname>
          </string-name>
          ,
          <article-title>Diagnostic methods for an aircraft engine performance</article-title>
          ,
          <source>Journal of Engineering Science and Technology</source>
          , vol.
          <volume>8</volume>
          , no.
          <issue>4</issue>
          (
          <year>2015</year>
          )
          <fpage>64</fpage>
          -
          <lpage>72</lpage>
          . doi:
          <volume>10</volume>
          .25103/jestr.084.10
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>B.</given-names>
            <surname>Li</surname>
          </string-name>
          .
          <string-name>
            <given-names>Y.-P.</given-names>
            <surname>Zhao</surname>
          </string-name>
          ,
          <article-title>Group reduced kernel extreme learning machine for fault diagnosis of aircraft engine</article-title>
          ,
          <source>Engineering Applications of Artificial Intelligence</source>
          , vol.
          <volume>96</volume>
          (
          <year>2020</year>
          )
          <article-title>103968</article-title>
          . doi:
          <volume>10</volume>
          .1016/j.engappai.
          <year>2020</year>
          .
          <volume>103968</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>H.</given-names>
            <surname>Hanachi</surname>
          </string-name>
          , J. Liu,
          <string-name>
            <surname>C.</surname>
          </string-name>
          <article-title>Mechefske, Multi-mode diagnosis of a gas turbine engine using an adaptive neuro-fuzzy system</article-title>
          ,
          <source>Chinese Journal of Aeronautics</source>
          , vol.
          <volume>31</volume>
          , issue
          <volume>1</volume>
          (
          <year>2018</year>
          )
          <fpage>1</fpage>
          -
          <lpage>9</lpage>
          . doi:
          <volume>10</volume>
          .1016/j.cja.
          <year>2017</year>
          .
          <volume>11</volume>
          .017
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>S.</given-names>
            <surname>Pang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Q.</given-names>
            <surname>Li</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.</given-names>
            <surname>Ni</surname>
          </string-name>
          ,
          <article-title>Improved nonlinear MPC for aircraft gas turbine engine based on semialternative optimization strategy</article-title>
          ,
          <source>Aerospace Science and Technology</source>
          , vol.
          <volume>118</volume>
          (
          <year>2021</year>
          )
          <article-title>106983</article-title>
          . doi:
          <volume>10</volume>
          .1016/j.ast.
          <year>2021</year>
          .106983
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>S.</given-names>
            <surname>Herzog</surname>
          </string-name>
          ,
          <string-name>
            <given-names>C.</given-names>
            <surname>Tetzlaff</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Worgotter</surname>
          </string-name>
          ,
          <article-title>Evolving artificial neural networks with feedback</article-title>
          ,
          <source>Neural Networks</source>
          , vol.
          <volume>123</volume>
          (
          <year>2020</year>
          )
          <fpage>153</fpage>
          -
          <lpage>162</lpage>
          . doi:
          <volume>10</volume>
          .1016/j.neunet.
          <year>2019</year>
          .
          <volume>12</volume>
          .004
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>M.</given-names>
            <surname>Lungu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Lungu</surname>
          </string-name>
          ,
          <article-title>Automatic control of aircraft lateral-directional motion during landing using neural networks and radio-technical subsystems</article-title>
          ,
          <source>Neurocomputing</source>
          , vol.
          <volume>171</volume>
          (
          <year>2016</year>
          )
          <fpage>471</fpage>
          -
          <lpage>481</lpage>
          . doi:
          <volume>10</volume>
          .1016/j.neucom.
          <year>2015</year>
          .
          <volume>06</volume>
          .084
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>F. D.</given-names>
            <surname>Amare</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S. I.</given-names>
            <surname>Gilani</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B. T.</given-names>
            <surname>Aklilu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Mojahid</surname>
          </string-name>
          ,
          <article-title>Two-shaft stationary gas turbine engine gas path diagnostics using fuzzy logic</article-title>
          ,
          <source>Journal of Mechanical Science and Technology</source>
          , vol.
          <volume>31</volume>
          (
          <year>2017</year>
          )
          <fpage>5593</fpage>
          -
          <lpage>5602</lpage>
          . doi:
          <volume>10</volume>
          .1007/s12206-017-1053-9
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>J.</given-names>
            <surname>Zeng</surname>
          </string-name>
          , Y. Cheng,
          <article-title>An Ensemble Learning-Based Remaining Useful Life Prediction Method for Aircraft Turbine Engine, IFAC-PapersOnLine</article-title>
          , vol.
          <volume>53</volume>
          , issue
          <volume>3</volume>
          (
          <year>2020</year>
          )
          <fpage>48</fpage>
          -
          <lpage>53</lpage>
          . doi:
          <volume>10</volume>
          .1016/j.ifacol.
          <year>2020</year>
          .
          <volume>11</volume>
          .009
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>D.</given-names>
            <surname>El-Masri</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Petrillo</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.-G.</given-names>
            <surname>Gueheneuc</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Hamou-Lhadj</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Bouziane</surname>
          </string-name>
          ,
          <article-title>A systematic literature review on automated log abstraction techniques</article-title>
          ,
          <source>Information and Software Technology</source>
          , vol.
          <volume>122</volume>
          (
          <year>2020</year>
          )
          <article-title>106276</article-title>
          . doi:
          <volume>10</volume>
          .1016/j.infsof.
          <year>2020</year>
          .106276
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>J.</given-names>
            <surname>Friederich</surname>
          </string-name>
          , G. Lugaresi,
          <string-name>
            <given-names>S.</given-names>
            <surname>Lazarova-Molnar</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Matta</surname>
          </string-name>
          ,
          <article-title>Process Mining for Dynamic Modeling of Smart Manufacturing Systems: Data Requirements</article-title>
          ,
          <string-name>
            <surname>Procedia</surname>
            <given-names>CIRP</given-names>
          </string-name>
          , vol.
          <volume>107</volume>
          (
          <year>2022</year>
          )
          <fpage>546</fpage>
          -
          <lpage>551</lpage>
          . doi:
          <volume>10</volume>
          .1016/j.procir.
          <year>2022</year>
          .
          <volume>05</volume>
          .023
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <given-names>S.</given-names>
            <surname>Kiakojoori</surname>
          </string-name>
          ,
          <string-name>
            <given-names>K.</given-names>
            <surname>Khorasani</surname>
          </string-name>
          ,
          <article-title>Dynamic neural networks for gas turbine engine degradation prediction, health monitoring and prognosis</article-title>
          ,
          <source>Neural Computing &amp; Applications</source>
          , vol.
          <volume>27</volume>
          , no.
          <issue>8</issue>
          (
          <year>2016</year>
          )
          <fpage>2151</fpage>
          -
          <lpage>2192</lpage>
          . doi:
          <volume>10</volume>
          .1007/s00521-015-1990-0
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [12]
          <string-name>
            <given-names>Y.</given-names>
            <surname>Shen</surname>
          </string-name>
          ,
          <string-name>
            <given-names>K.</given-names>
            <surname>Khorasani</surname>
          </string-name>
          ,
          <article-title>Hybrid multi-mode machine learning-based fault diagnosis strategies with application to aircraft gas turbine engines</article-title>
          ,
          <source>Neural Networks</source>
          , vol.
          <volume>130</volume>
          (
          <year>2020</year>
          )
          <fpage>126</fpage>
          -
          <lpage>142</lpage>
          . doi:
          <volume>10</volume>
          .1016/j.neunet.
          <year>2020</year>
          .
          <volume>07</volume>
          .001
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          [13]
          <string-name>
            <given-names>B.</given-names>
            <surname>Li</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.-P.</given-names>
            <surname>Zhao</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.-B.</given-names>
            <surname>Chen</surname>
          </string-name>
          ,
          <article-title>Unilateral alignment transfer neural network for fault diagnosis of aircraft engine</article-title>
          ,
          <source>Aerospace Science and Technology</source>
          , vol.
          <volume>118</volume>
          (
          <year>2021</year>
          )
          <article-title>107031</article-title>
          . doi:
          <volume>10</volume>
          .1016/j.ast.
          <year>2021</year>
          .107031
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          [14]
          <string-name>
            <given-names>Z.</given-names>
            <surname>Wu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Li</surname>
          </string-name>
          ,
          <article-title>A Framework of Dynamic Data Driven Digital Twin for Complex Engineering Products: the Example of Aircraft Engine Health Management, Procedia Manufacturing</article-title>
          , vol.
          <volume>55</volume>
          (
          <year>2021</year>
          )
          <fpage>139</fpage>
          -
          <lpage>146</lpage>
          . doi:
          <volume>10</volume>
          .1016/j.promfg.
          <year>2021</year>
          .
          <volume>10</volume>
          .020
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          [15]
          <string-name>
            <given-names>O.</given-names>
            <surname>Avrunin</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Vladov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Petchenko</surname>
          </string-name>
          ,
          <string-name>
            <given-names>V.</given-names>
            <surname>Semenets</surname>
          </string-name>
          ,
          <string-name>
            <given-names>V.</given-names>
            <surname>Tatarinov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>H.</given-names>
            <surname>Telnova</surname>
          </string-name>
          ,
          <string-name>
            <given-names>V.</given-names>
            <surname>Filatov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Shmelov</surname>
          </string-name>
          and
          <string-name>
            <given-names>N.</given-names>
            <surname>Shushlyapina</surname>
          </string-name>
          .
          <article-title>Intelligent automation systems</article-title>
          , Kremenchuk, Novabook,
          <year>2021</year>
          . 322 p. doi:
          <volume>10</volume>
          .30837/
          <fpage>978</fpage>
          -617-639-347-4
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          [16]
          <string-name>
            <given-names>S.</given-names>
            <surname>Vladov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>I.</given-names>
            <surname>Dieriabina</surname>
          </string-name>
          ,
          <string-name>
            <given-names>O.</given-names>
            <surname>Husarova</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Pylypenko</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Ponomarenko</surname>
          </string-name>
          <article-title>, Multi-mode model identification of helicopters aircraft engines in flight modes using a modified gradient algorithms</article-title>
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>