=Paper= {{Paper |id=Vol-3664/paper5 |storemode=property |title=Monitoring State of Marine Plain Bearings Based on Exponential Degradation Model |pdfUrl=https://ceur-ws.org/Vol-3664/paper5.pdf |volume=Vol-3664 |authors=Artem Sharko,Oleksandr Sharko,Dmitro Stepanchikov,Artem Yanenko |dblpUrl=https://dblp.org/rec/conf/colins/SharkoSSY24 }} ==Monitoring State of Marine Plain Bearings Based on Exponential Degradation Model== https://ceur-ws.org/Vol-3664/paper5.pdf
                         Monitoring State of Marine Plain Bearings Based on
                         Exponential Degradation Model
                         Artem Sharko1, Oleksandr Sharko2, Dmitro Stepanchikov3and Artem Yanenko2
                         1 Technical University of Liberec, Studentská 1402/2, Liberec, 46117, Czech Republic
                         2 Kherson State Maritime Academy, 20, Ushakov ave., Kherson, 73000, Ukraine
                         3 Kherson National Technical University, 24 Beryslavske Shose, Kherson, 73008, Ukraine



                                         Abstract
                                         Some aspects of mathematical, algorithmic, and software apparatuses for modelling ship power
                                         equipment condition monitoring are considered. A predictive model of technical state change by
                                         parameters of generated vibration is developed. The principal components method is used to reduce
                                         the dimensionality of input information. The results of digitalisation and spectral analysis of
                                         vibroacoustic characteristics generated during the operation of dynamic equipment of ship sliding
                                         bearings are presented. The sharp increase of the first principal component compared with the rest
                                         and its monotonous increase at degradation of material properties are established. It is shown that
                                         developing predictive models of technical state change by vibration parameters is a major step in the
                                         transition to qualitatively new forms of maintenance and repair necessary for equipment safety.

                                         Keywords
                                                                                                                                              1
                                         monitoring, exponential model, degradation, prediction, residual life, principal component method.


                         1. Introduction
                         Ship sliding bearings with antifriction layers made of tin and lead-based alloys are used on ship low-
                         and medium-speed diesel engines, turbines, and shafts. By their design and direct purpose, they
                         support rotating parts and are the most important units of power equipment. Nowadays, methods of
                         vibration signal processing in frequency and frequency-time domains are used to estimate the
                         technical condition of plain bearings. The most widespread are spectral methods of vibration
                         diagnostics on the basis of the Fourier transform, where each defect of the equipment corresponds to a
                         set of discrete frequencies. Because defects of plain bearings generate both low-frequency and high-
                         frequency vibration, their diagnostics are carried out by joint analysis of the spectrum of low-
                         frequency vibration and envelope spectrum of the vibration signal. Monitoring of technical conditions
                         and diagnostics of arising defects is the basis for ensuring the high reliability of equipment elements.
                         According to the totality of defective frequencies detected in the spectrum a conclusion is made about
                         the defect and the degree of its development. The disadvantage of spectral methods is the low
                         efficiency of non-stationary signal processing due to the blurring of spectra. The causes of defects are
                         the variability of the shaft rotation speed when changing the operating mode, different degrees of
                         wear of elements, etc.
                             Physical processes occurring in sliding bearings depend on the design's peculiarities and the
                         correlation of many external and internal factors determining bearing operating conditions. The
                         operation of plain bearings is based on the principles of sliding friction, the realisation of which
                         requires constant lubrication control. The inner ring is made of an antifriction material to
                         ensure a low coefficient of friction, and the outer ring is made of high-strength material. The
                         advantages of plain bearings are their resistance to radial loads and vibrations, their

                         COLINS-2024: 8th International Conference on Computational Linguistics and Intelligent Systems, April 12–13, 2024,
                         Lviv, Ukraine
                            sharko_artem@ukr.net (A. Sharko); mvsharko@gmail.com (O. Sharko); dmitro_step75@ukr.net (D. Stepanchikov);
                         yanenko9494@gmail.com (A. Yanenko)
                            0000-0002-6350-7189 (A. Sharko);0000-0001-9025-7990 (O. Sharko); 0000-0001-5027-2213 (D. Stepanchikov);
                         0009-0004-7992-8369 (A. Yanenko)
                                  © 2024 Copyright for this paper by its authors.
                                  Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0).




CEUR
                  ceur-ws.org
Workshop      ISSN 1613-0073
Proceedings
separability, which simplifies installation and maintenance, and their lower operating noise.
The disadvantage is the high wear of contacting surfaces. The main reason for the occurrence of
defects in plain bearings is abrasive wear of mating surfaces of shaft journal and liners, as a
result of getting on the contact surface of various particles, being in lubricants and load
differences, leading to semi-dry friction at starting and stopping of the bearing. Defects in the
bearing shells indicate a malfunction of the diesel engine. If the bearing shell has been in service
for a long time, various defects will appear on its surfaces. Under the weight of the rotor and
other static loads in the plain bearing, oscillating forces of kinematic and impulse character can
occur caused by the action of friction forces. Besides, these forces in plain bearings can be
observed as vibrations with a rotation frequency of 0,42...0,48, excited by auto-oscillations of
the rotor on the oil wedge. The short-term appearance of cavitation zones causes shocks of
hydrodynamic origin.
    The time interval until the first signs of material fatigue appear depends on the bearing
speed and load. The actual durability of plain bearings is lower than the calculated durability
due to high loads, insufficient lubrication, incorrect choice of lubricant, and installation errors.
    The relevance of the work is to investigate the possibility and develop a methodology to
observe the development of defects in a plain bearing in real-time.
    The exponential model of material degradation is chosen as the hypothesis used in the
analysis of monitoring of vibration signals received during the operation of sliding bearings.
    The experience of creating models of degradation of the technical state of dynamic
equipment based on the analysis of mechanical vibrations is described in [1]. In [2], the results
of vibration characteristics are considered in detail, and the classification of defects is given,
based on which specific vibration analysis methods are justified.
    In [3], the statistical analysis of failures and the causes of their occurrence is described.
Peculiarities of monitoring transport technology elements are devoted to [4-6]. It is shown that
the timely detection of the origin and development of defects in bearings, compressors, and
turbines of high-speed technological equipment is the basis for the safety of transport devices.
Frequency-time analysis of the steady-state helpful life of bearings is presented in [7]. An
improved strength degradation model for fatigue life prediction considering material
characteristics is presented in [8]. The model has high prediction accuracy, but its
implementation needs to be improved to account for the influence of associated factors. The
prediction of the remaining service life of a self-lubricating liner of a high-frequency plain
bearing is presented in [9]. However, it should be noted that the problem of anomaly evaluation
in trends still needs to be addressed. Life prediction of self-lubricating spherical plain bearings
based on a failure physics model and accelerated degradation tests is presented in [10]. The
accuracy of the forecast is determined by the number of diagnostic measurements, which
restrains the broad application of the method.
    The information support system for monitoring the development of equipment degradation
state is presented in [11-13]. It is shown that the organisation of monitoring observation of the
state of the operating equipment is a means of increasing its reliability and failure-free
operation. Some scientific interest in solving the problem of identifying the state of shipboard
equipment is represented by [14-16]; intellectual support of forecasts of approaching the state
of destruction and estimation of residual life is presented in [17-19].
    The given review of literature sources shows that, despite the variety of scientific directions
to improve the efficiency of equipment monitoring in the process of its operation, the interest in
publications in this area is particularly noticeable. The purpose of the present work is to
monitoring the condition of ship sliding bearings on the basis of the exponential degradation
model.

2. Materials and Methods
Methods of analysing vibrodiagnostic signals with an exponential degradation model were applied to
monitoring state of marine plain bearings.
   The exponential degradation model forecasts residual life based on prior information about
equipment operation conditions and the latest diagnostic measurements. The model allows for
observing the degradation trend in real time and updating its priori parameters, when new
information is received using data collection, post-processing, ranking, and merging essential
functions with a selection of forecasting curves. This realises the principle of observing the
continuous wear of the equipment material during its lifetime.
   Materials for analysing vibrodiagnostic signals ship sliding bearings of mark ТРL-88 were
used. Their parameters are given in Table 1.

Table 1
Parameters of bearing ТРL-88
Maximum shaft Bearing          inner Bearing     outer Weight, kg               Anti-friction
speed, rpm         diameter, mm      diameter, mm                               coating
12000              114.073           164.672           7.7                      CuPb15Sn

   The arrangement of the plain bearings on the motor is shown in Fig. 1.




Figure 1: Placing of plain bearings on the engine: 1 - plain bearing, 2 - compressor, 3 – turbine

   Measurements are carried out on the bearing unit housing, namely in its lower part, because
here the loads on the unit are maximum [20]. The signals from the sensors can be digitised and
recorded for trend analysis. An accelerometer is used to record the vibration levels. A vibration
signal of 6 s duration was obtained daily for 50 consecutive days. A bearing failure occurred and
caused the bearing to fail.

3. Methodology
The exponential degradation model at its inception assumes that the failure rate function
failures λ(t) is constant during the entire life of the bearing.
                                          λ (t ) = const > 0                           (1)
    Probability of its failure-free operation
                                          P (t , λ ) = exp(− λt )                      (2)
    Accordingly, the probability failures of bearing failures
                              Q(t , λ ) = 1 − P(t , λ ) = 1 − exp(− λt )               (3)
    In the framework of such a model, aging and wear processes are absent. This model is used
to approximate the unknown failure rate at the stage of normal operation under the assumption
that the failure rate remains constant.
    The failure rate during damage accumulation is described by a Weibull distribution
                             λ (t ,α ) = α (λ , t ) , t > 0,α > 0, λ > 0               (4)
                                                   α −1


   where α - intensity index.
    The probability of failure-free operation for time t is equal to
                                    Q(t ) = 1 − exp(− (λt ) )                               (5)
                                                            α


    If α>1 the failure rate function increases monotonically, which describes the processes of
wear and ageing of equipment.
    Maintenance based on interactive condition monitoring, fault detection and resource
forecasting of power equipment requires a whole set of information characteristics: amplitude,
time, frequency, etc. The method of principal components allows to reduce their dimensionality
without significant loss of input information. Based on linear algebra and mathematical
statistics, the method of principal components, which works with continuous data streams,
allows us to avoid multidimensionality and select only the main characteristics of objects. It is
designed to divide the matrix of initial data into two parts: meaningful and noise. The principal
component method approximates an n-dimensional cloud of observations to an n-dimensional
ellipsoid, the semi-axes of which will be the principal components. The exponential model
allows real-time prediction of residual life based on statistical processing of data and their
integration. In this paper, the method of principal components is applied to the data written in
the form of a matrix of numbers of dimensionality n×m. In preparing the data, the sample is
centred so that the mean of the features is zero by replacing the finite set of points with lines
and planes. For a finite set of vector sets х1,х2,.., хn∈Rn where Rn – linear manifolds, we need to
find such values of the set of linear combinations of Sk⊂Rn, so that the sum of squares of
deviations xi from Sk would be minimal.
                                     m
                                       dist 2 ( x , S ) → min                               (6)
                                   ∑i =1
                                               i   k


    where dist (xi, Sk) is the Euclidean distance from a point to a line, k = 0,1,2…n-1.
    Linear varieties are represented by a set of principal components. When describing a
random variable, the mathematical expectation representing the centre of gravity of this
variable and dispersion representing its dimensions in the form of spread are used. To describe
a multidimensional special quantity at occurrence of vibration signals, in addition to the
mathematical expectation Е(х) and dispersion of its projections on the axis, we used the concept
of covariance matrix, the elements of which are correlations of the features хi and xj.
                                           [           (      )]
                             cov( xi , x j ) = E (xi − E xi ) ⋅ x j − E x j               (7)
    In the covariance matrix, the eigenvalues are used to estimate the contribution of the
principal components to the overall variability of the process.
    Before proceeding to the calculation of principal components, it is necessary to standardise
the data to zero mean and unit variance. The algorithm of the principal component method is an
orthogonal linear transformation that maps the data from the original feature space into a new
space of linear dimension. By projecting the principal components of the axis, a new basis is
formed, dimensionality is reduced, and the greatest amount of information is retained. This
method is incorporated into analytical platforms and used based on data preprocessing. Having
centred the input data sample, we shift it linearly so that the mean values of the features are
zero. To describe the shape of the random vector, we need a covariance matrix that has i,j
elements are a correlation of attributes хi and xj. The covariance matrix is a generalisation of the
variance of a multivariate random variable.
    The calculation of principal components is reduced to the calculation of eigenvectors and
eigenvalues of the covariance matrix. In the matrix on the diagonal will be the variance of the
features, and in the remaining cells - the covariance of pairs of these features.
    Modelling of the exponential probability distribution of the diagnostic parameters is
conveniently performed using the logarithmic subharmonic function to preserve the scale. A
random variable obeys the logarithmic normal distribution θ, if its logarithm has a normal
distribution.
    The formula for the density of the log-normal distribution of a random variable has the
following form
                                       1 1         1  ln x − m  2                        (8)
                          P( x ) =          + exp −             
                                     σ 2p x        2 σ          
   σ,m – distribution parameters.
   The logarithmic normal distribution of a random variable is asymmetric when compared to
the Gaussian normal distribution of a random variable β.
   The exponential degradation model is defined as
                                                           σ2                    (9)
                              h (t ) = φ + θ exp βt + e −    
                                                           2 
   where h(t) – bearing state indicator as a function of time; φ – tipping constant; θ and β are
random parameters that determine the slope of the model, where θ is a logarithmic normal
distribution, and β is Gaussian distributed. At each time step t distribution θ and β is updated to
a posteriori based on the latest observation of h(t). Parameter e is a Gaussian white noise, σ 2 –
variance.
   Member        exponentially should force the expectation h(t) fulfil the condition

                                                                                         (10)
   In information theory, Gaussian white noise is an abstract mathematical model of a
stationary random process with constant spectral density at all frequencies, which is a mixture
of sounds reproduced simultaneously at all frequencies. Gaussian white noise is completely
uncorrelated, i.e. any instantaneous noise value is uncorrelated with the previous one.
   Degradation models extrapolate past behaviour to predict future behaviour. This type of
calculation is suitable when determining the steady-state useful life of complex equipment. The
degradation profile of the test component is then used to statistically calculate the remaining
time until the diagnostic indicator reaches some level of degradation consistent with regulatory
guidelines.

4. Experiment, results and discussion
Fig.2 shows in one scale the whole set of received vibration signals in the order of their occurrence
and fixation in the course of 50-day operation life of the plain bearing, which gives a visual
representation of the change of peak amplitudes.




Figure 2: Sequence of received vibration signals in the order of their occurrence and fixation during
the operation of the plain bearing

    As can be seen from the consideration of Fig.2, the vibration signals in the time domain show
a tendency to increase the impulsiveness of the signal.
    Statistical characteristics of vibration signals in time and frequency domains were
determined during post-processing. In the time domain, 11 statistical characteristics were
determined: Mean, Std, Skewness, Kurtosis, Peak2Peak, RMS, CrestFactor, ShapeFactor,
ImpulseFactor, MarginFactor, Energy. In the frequency domain 4 statistical characteristics were
determined: spectral mean (SKMean), spectral standard deviation (SKStd), spectral asymmetry
(SKSkewness) and spectral excess (SKKurtosis). All the above statistical characteristics of
vibration signals can serve as potential indicators of bearing condition degradation (Fig.2).
   The extracted statistical characteristics of vibration signals are related to noise. Noise is
detrimental for residual life prediction. In addition, one of the most important properties of the
function, its monotonicity, is unstable to noise. Therefore, a filtering and smoothing procedure
is applied to the extracted statistical features. The comparison of the original and filtered
characteristics of vibration signals is presented in Fig.3.




Figure 3: Evolution of dimensionless statistical characteristics of vibration signals in the process of
sliding bearing operation (blue line - initial values, red - smoothed values)
   The main criterion of belonging of the extracted statistical characteristics of vibration signals
in the frequency and time domain to their further analysis is the numerical value of
monotonicity. For quantitative estimation of the belonging of statistical characteristics to their
further analysis it is proposed to use the formula
                                                                                             (11)

     where n – number of measurement points, in our case n = 50. m – number of controlled
samples, in our case m = 1, xij – i-th characteristic measured on the j-th sample, diff(xji) = xji(t) –
 xji(t – 1). The results of calculation of monotonicity of statistical characteristics of vibration
signals in frequency and time domain by formula (11) are presented in Fig.4.




Figure 4: Monotonicity of statistical characteristics of vibration signals arising in the process of
operation of the sliding bearing

    Characteristics with a monotonicity score greater than 0.3 are selected to combine them in
the subsequent analysis using the principal component method. As the monotonicity
calculations have shown, the following values greater than 0.3 were obtained: standard
deviation (Std), full range of oscillations (Peak2Peak), RMS value (RMS), energy (Energy),
marginal factor (MarginFactor), excess (Kurtosis), shape factor (ShapeFactor), standard spectral
deviation (SKStd).
    In order to reduce the dimensionality of the analysed quantities and to combine the features
of vibration signals, the PCA principal component analysis method was applied. Using the
computer mathematics system Matlab 2018b, the first PCA principal component1 was
calculated when processing the statistical characteristics of vibration signals with monotonicity
greater than 0.3. Its distribution depending on the measurement day, which correlates with the
successive approach of the plain bearing state to failure, is shown in Fig. 5.
    The graph in Fig.5 shows that the first principal component increases as the bearing
approaches failure. Thus, the first principal component is a promising combined indicator of the
condition of a plain bearing.
    In practice, the data of the whole life cycle of a plain bearing is not enough to develop a
prognostic algorithm, but it can be considered as data for training a system for monitoring the
condition of marine plain bearings during their operation. Therefore, the data collected in the
first 20 days (40% of the life cycle) are considered as training data. In the subsequent
operations of ranking and combining the importance of features, only the training data are used.
Figure 5: The first principal component for successively recorded vibration signals during the
operation of a plain bearing and its approaching to the failure state

    In this paper, the Echronic Degradation Model algorithm implemented in the Matlab 2018b
computer mathematics system is used to model the exponential degradation process and
estimate the Remaining Useful Life (RUL) of a bearing. The degradation models estimate RUL by
predicting when the monitored signal crosses a predetermined threshold value.
    The selection of the condition indicator threshold value may be based on chronological
records of bearing condition or some subject matter knowledge. Since there is no chronological
data in the study set, the last value of the health indicator is selected as the threshold value h(t).
The threshold value is selected on the basis of smoothed chronological data.
    As follows from the graphs Fig.2,5 the study of statistical characteristics of vibration signals
reveals predominantly ascending trend with increasing operation time of the sliding bearing. A
priori of slope parameters varies within the limits (E(θ) = 1, Var(θ) = 106, E(β) = 1,
Var(β) = 106), where 106 – scale factor of the distribution. The model relies mainly on observed
data. Based on the relationship E[h(0)] = φ + E(θ), intersection φ is set to -1 so that the model
also starts with a 0.
    The relationship between the change in state index and the change in noise can be obtained
as follows:
                                                                                              (12)
    Here it is assumed that the standard deviation of noise causes a 10% deviation in the
performance indicator when it approaches a threshold value.
    The exponential degradation model also provides functionality to assess the significance of
the slope. Once a significant slope of the performance indicator is detected, the model will forget
previous observations and restart the estimation based on the original a priori values.
    Fig.6 shows the operation of the exponential model of the sliding bearing state degradation
on the example of the vibration signals considered above. It is clearly seen how the prediction
changes depending on the slope of the state indicator curve. The threshold value of the state
indicator is chosen equal to h(t) = 36 on the 20th day of monitoring, as the first twenty days
were used to train the model.
Figure 6: Estimation of the residual life of a plain bearing based on the exponential model of state
degradation (red line - threshold value of the state indicator, blue line - exponential degradation
model, yellow line - current value of the state indicator, dotted line - confidence interval limits)

    As can be seen from Fig.6, if the data for the first twenty days of operation of the plain bearing are
used for training of the exponential model, the exponential model gives exceeding of the threshold
value on the fiftieth day. However, it should be noted that the boundary of the confidence interval,
shown by the dotted line in Fig.6, crosses the threshold already on the fortieth day of the bearing
operation. Therefore, it can be stated that after forty days of operation, further operation of the bearing
is unsafe.
    The chronological data on the technical state of the object may include structural and functional
schemes of subsystems interaction, control algorithms, data on modes, characteristics of control and
diagnostics systems, methods of operation and technical support. The disadvantage of such an
assessment is the lack of uncertainty accounting, which is characterised by the actual results of the
diagnostics device operation, as well as the assessment of the influence of external influences.
    For predictive performance analysis the following is used α-λ graph shown in Fig. 7, where the
threshold is set to α=20%. The probability that the estimated RUL is between the boundary of the true
RUL is calculated as a measure of model performance:
                                                                                                    (13)
    where r(t) - estimated RUL during t, r*(t) - true RUL during t, Θ(t) - are the estimated
parameters of the model during t.




Figure 7: α-λ graph used for prognostic analysis of the performance of a plain bearing
   Since the preset a priori values do not reflect the true a priori value, the model usually needs
several time steps to adjust to the correct parameter distribution. Emergence, recognition, and
evolution are statistical characteristics of prediction that require a large amount of input information.
In these states, the initial position is uncertain and therefore the statistical characteristics cannot be
accurately calculated. Improving the accuracy of prediction requires a constant replenishment of
current information at each RUL prediction step, estimated by the corresponding probabilities. The
prediction becomes more accurate as more data points become available.

    Conclusions
       1. A framework for conducting turbine plain bearing condition monitoring based on
          vibration signal measurements is proposed. A set of statistical characteristics of
          vibration signals in time and frequency domains has been identified, and their
          statistical processing and analysis has been performed.
       2. On the basis of the performed calculations of the sliding bearings state monitoring
          using the method of principal components a new information diagnostic parameter
          in the form of the first principal component is established, which allows to combine a
          set of statistical characteristics of vibration signals at reduction of the noise
          component.
       3. The exponential degradation model used in the paper, taking into account
          probabilistic uncertainty estimates and confidence estimation of predictions, has
          shown that for cases in which a normal distribution is used it is sufficient to use a
          confidence interval in the form of a propagating band around the predicted point and
          to estimate the RUL based on a predetermined failure threshold.
       4. It is found that forecasts made at an early stage have access to less information on
          damage dynamics and require forecasting to more distant horizons.

    References
[1] L.-L. Cui, O. Jin, X. Wang Rolling bearing life prediction method based on improved
    similarity theory. Zhendong Gongcheng Xuebao. Journal of Vibration Engineering. (2023)
    36(3), с. 854-860 DOI: 10.16385/j.cnki.issn.1004-4523.2023.03.028
[2] L. Liangxing, M. Guozheng, S. Jianfang, H. Cuihong, Y. Qingsong, S. Fenghua, W. Haidou
    Remaining Useful Life Prediction Method of Coated Spherical Plain Bearing Based on VMD-
    EEMD-LSTM. Jixie Gongcheng Xuebao/Journal of Mechanical Engineering. (2023) 59(9), с.
    125-136 DOI: 10.3901/JME.2023.09.125
[3] B.A. Abass, S.Y. Ahmed, Z.H. Kadhim Analysis and Optimization of Nanolubricated Journal
    Bearing under Thermoelasto-Hydrodynamic Lubrication Considering Cavitation Effect.
    Tribology in Industry. (2023) 45(4), с. 618-633 DOI: 10.24874/ti.1441.01.23.04
[4] M. Qiu, K. Niu, J. Li, Y. Li, Y. Xu Residual life prediction of rolling bearing under multiple
    stresses. Hangkong Dongli Xuebao/Journal of Aerospace Power. (2022) 37(5), с. 980-988
    DOI: 10.13224/j.cnki.jasp.20210246
[5] C. Zhang, S. Li, J. Wen, S. Fu Remaining Useful Life Prediction Method Based on the
    Improved Holt Double Exponential Model. 2022 5th International Conference on Pattern
    Recognition       and       Artificial  Intelligence,    PRAI    (2022)     с.   47-52     DOI:
    10.1109/PRAI55851.2022.9904057
[6] Y. Xu, M. Qiu, J. Li, L. Liu, K. Niu Remaining useful life prediction method of rolling bearing
    based on SKF-KF-Bayes. Zhendong yu Chongji/Journal of Vibration and Shock. (2021)
    40(19), с. 26-31 and 40 . DOI:10.13465/j.cnki.jvs.2021.19.004
[7] D. Huang, G. Yu, J. Zhang, J. Tang, J. Su An Accurate Prediction Algorithm of RUL for
    Bearings: Time-Frequency Analysis Based on MRCNN. Wireless Communications and
    Mobile Computing. (2022), 2222802. DOI:10.1155/2022/2222802
[8] M. Zhang, G.Hu, X.Liu, X.Yang An improved strength degradation model for fatigue life
     prediction considering material characteristics . Journal of the Brazilian Society of
     Mechanical Sciences and Engineering. Journal of the Brazilian Society of Mechanical
     Sciences and Engineering. (2021) 43(5),275. DOI:10.1007/s40430-021-02997-4
[9] X. Hao, S. Wang, M. Chen, D. Pan, Remaining Useful Life Prediction of High-Frequency Swing
     Self-Lubricating      Liner.          Shock      and       Vibration.     (2021),     8843374.
     DOI:10.1155/2021/8843374
[10] Y. Wang, X. Fang, C. Zhang, X. Chen, J.Lu Lifetime prediction of self-lubricating spherical
     plain bearings based on physics-of-failure model and accelerated degradation test |
     [Prognozowanie czasu pracy samosmarujących łożysk ślizgowych w oparciu o model fizyki
     uszkodzeń oraz przyspieszone badania degradacji]. Eksploatacja i Niezawodnosc. (2016)
     18(4), с. 528-538. DOI: 10.17531/ein.2016.4.7
[11] A. Abuelo, J. Brester, K. Starken, L. Neuder. Technical note: Comparative evaluation of 3
     methods for the quantification of nonesterified fatty acids in bovine plasma sampled
     prepartum. Journal of Dairy Science. 2020. Vol. 103 (3), P. 2711-2717
[12] T. Noguerol, R. Barousse, T. Amrhein, J. Royuela-Del-val, P. Montesinos, A. Luna Optimizing
     diffusion-tensor imaging acquisition for spinal cord assessment. Physical basis and
     technical adjustments Radiographics. 2020. 40 (2), Р. 403-427
[13] V. Levin, A. Yahya, Adaptive management of technical condition of power transformers.
     International Journal of Electrical and Computer Engineering. 2020. 10 (4), Р. 3862-3868
[14] S. Zinchenko, V. Kobets, O. Tovstokoryi, K. Kyrychenko, P.Nosov, I. Popovych. Control of the
     Pivot Point Position of a Conventional Single-Screw Vessel, CEUR-WS.org, Vol.3513, P.130-
     140, 2023 (ICST-2023). https://ceur-ws.org/Vol-3513/paper11.pdf
[15] S. Zinchenko, K. Kyrychenko, O. Grosheva, P. Nosov, I. Popovych, P. Mamenko Automatic
     reset of kinetic energy in case of inevitable collision of ships, IEEE Xplore, p.496-500, 2023
     13th International Conference on Advanced Computer Information Technologies (ACIT),
     Wrocław, Poland, 2023. DOI: 10.1109/ACIT58437.2023.10275545.
[16] S. Zinchenko, V. Kobets, O. Tovstokoryi, P. Nosov and I. Popovych. Intelligent System
     Control of the Vessel Executive Devices Redundant Structure //CEUR Workshop
     Proceedings, Vol-3403, pp. 582-594, 2023. https://ceur-ws.org/Vol-3403/
[17] N. Petrushenko, M. Sharko, O. Gonchar, N. Vasylenko, K. Vorobyova, I. Zakryzhevska
     Information Support of Intelligent Decision Support Systems for Managing Complex
     Organizational and Technical Objects Based on Markov Chains CEUR Workshop
     Proceedings, 2022, 3171, рр. 986-998. https://ceur-ws.org/Vol-3171/ (https://ceur-
     ws.org/Vol-3171/paper71.pdf)
[18] M. Sharkо, O. Gonchar, M. Tkach, A. Polishchuk, N. Vasylenko, M. Mosin, N. Petrushenko
     Intellectual Information Technologies of the Resources Management in Conditions of
     Unstable External Environment. Book Chapter. Lecture Notes on Data Engineering and
     Communication Technologies, 2022, 77, 1158, рр. 519-533. DOI:10.1007/978-3-030-
     82014-5
[19] C. Surace, A. Bovsunovsky The use of frequency ratios to diagnose structural damage in
     varying environmental conditions. Mechanical Systems and Signal Processing. 2020.
     Vol.136, 106523
[20] P. Louda, V. Marasanov, O. Sharko, D. Stepanchikov, A. Sharko The Theory of Similarity and
     Analysis of Dimensions for Determining the State of Operation of Structures Under Difficult
     Loading Conditions. Materials 2022, 15(3), 1191. DOI:10.3390/ma15031191.
[21] P. Louda, O. Sharko, D. Stepanchikov, A. Sharko Experimental and Theoretical Study of
     Plastic Deformation of Epoxy Coatings on Metal Substrates Using the Acoustic Emission
     Method. Materials 2022, 15(11), 3791. DOI: 10.3390/ma15113791.