Abstract
Based on regular analysis of the spectrum, PQ value, precision moisture content and kinematic viscosity of a certain sliding bearing lubricating oil, this paper used grey theory to study the correlation degree within these four kinds of analysis indicators. It is hoped that it can achieve mutual verification and judgment of the state evaluation of sliding bearings. It can provides reliable basis for studying the internal operating rules of sliding bearings. It can also provide scientific basis for the management and maintenance of sliding bearings.
1. Introduction
Sliding bearing was a complex frictional system with uncertainty [1-2]. Generally speaking, there were many theoretical methods for solving uncertainty problems, such as grey theory [3], probability and statistics theory [4], fuzzy theory [5] and rough set theory [6]. They were applied in different fields. Grey theory had distinct advantages in solving problems with unclear connotations among them [7]. It could use simple methods to study complex systems and obtain reliable analysis conclusions. In recent years it has been widely used in the fields of ships, railways, aviation and aerospace. Overall, the main application methods of grey system theory currently include grey correlation analysis [8], grey forecasting prediction [9], grey decision-making [10] and grey control.
Correlation degree was called as a measure of the magnitude of the correlation between two systems or two factors [11]. It described the relative changes in factors during the development process of the system. It was considered that the correlation was high if the relative changes between the two were basically the same during the system development process. Correlation analysis was a part of grey system theory which used grey correlation to describe the strength, magnitude and order of relationship between sequences. The greater the correlation was, the stronger the correlation was between the change factor and the reference factor. This reflected the greater the degree of mutual influence among various factors in the system.
The sliding bearings were mainly used to support the weight of each shaft and coupling, as well as the radial load generated during the rotation of the rotation of the shaft system, to ensure the normal operation of the shaft system. A sliding bearing was mainly consisted by shell components, bearing components, oil filling components, and connecting components. The used lubricating oil was Great Wall TSE68# turbine oil.
This article used grey theory to calculate the behavioral characteristics of various groups of lubricating oil in sliding bearings. The correlation degree was calculated using the measurement data of wear element content, PQ value, precision moisture and 40 ℃ kinematic viscosity of lubricating oil as the mother sequence and the measurement data of other items as the sub sequence. The purpose was to deeply study the relationship between different monitoring parameters in sliding bearing lubricating oil and provide scientific basis for the management and maintenance of sliding bearings.
2. Materials and methods
2.1. Acquisition of experimental data
Tracking and monitoring were carried out by using a sliding bearing lubricating oil as the monitoring object. Samples were taken and analyzed every 300 hours of operation. Each oil sample was tested for spectrum, PQ value, precision moisture, and 40℃ kinematic viscosity. Six main wear elements such as Fe, Cr, Pb, Cu, Sn, and Ni were mainly analyzed in the spectral analysis (unit: μg/g). The used instrument was a M-type oil emission spectrometer produced by American Hyperspectral Company. The PQ value was obtained using a TTL-M ferromagnetic abrasive particle tester produced by Shenzhen Yateks Optical Electronic Technology Co., Ltd. The precision moisture content (unit: μg/g) was tested using DL32 precision moisture tester produced by Mettler Toledo Technology Co., Ltd. The kinematic viscosity (unit: mm2/s) was tested using S-FLOW850 automatic viscometer produced by Omnitek in the Netherlands. The results were shown in Table 1.
Table 1The data of oil monitoring
S/N | Fe | Cr | Pb | Cu | Sn | Ni | PQ | Moisture | Viscosity |
1 | 8.9 | 0.1 | 3.8 | 2.3 | 28 | 0.3 | |||
2 | 3.6 | 0.2 | 2.6 | 1.7 | 16.2 | 0 | |||
3 | 0.7 | 0.1 | 1.6 | 2.5 | 26.9 | 0 | 29 | 1588 | 66 |
4 | 7.2 | 0.2 | 0 | 2.5 | 43.7 | 1.2 | 42 | 163 | 67 |
5 | 1.5 | 0.1 | 0.5 | 2.1 | 26.3 | 0 | 32 | 283 | 64 |
6 | 8.0 | 0 | 0 | 1.4 | 22.1 | 0.1 | 35.5 | 536 | 68 |
7 | 3.5 | 0 | 0.2 | 1 | 7 | 0 | 23.5 | 328 | 61 |
8 | 10.4 | 0 | 0 | 1.9 | 26 | 0.7 | 30 | 367 | 68 |
9 | 6.6 | 0 | 0.4 | 1.8 | 17.9 | 0.3 | 45 | 717 | 65 |
10 | 5.0 | 0 | 0 | 0 | 15.7 | 0 | 45 | 159 | 65 |
11 | 9.7 | 0 | 0.9 | 1.6 | 22.7 | 0.3 | 30.5 | 92 | 65 |
12 | 9.4 | 0.1 | 0.1 | 1.3 | 22.8 | 0.4 | 44.5 | 115 | 64 |
13 | 7.8 | 0.5 | 1.6 | 1.2 | 20.3 | 0.1 | 31 | 948 | 64 |
14 | 7.2 | 0.1 | 0.3 | 0.9 | 16.5 | 0 | 18 | 384 | 64 |
15 | 8.8 | 0.5 | 2.0 | 1.2 | 19.5 | 0 | 36.5 | 641 | 66 |
16 | 10.7 | 0.1 | 10.1 | 10.2 | 19.2 | 0 | 37 | 421 | 66 |
17 | 10.3 | 0 | 1.1 | 6 | 19 | 0.5 | 26.5 | 200 | 67 |
Fig. 1Variation of wear element content in different metals

The changing trends of the six main wear elements, Fe, Cr, Pb, Cu, Sn, and Ni were showed in Fig. 1. It could be seen from Fig. 1 that the Fe element content was higher in the 8th, 16th and 17th tests. The Pb and Cu element content both were higher in the 16th test. The Sn element content was higher in the 4th test. The measured values of Cr and Ni element content were both smaller in each test.
The change of trend chart of PQ value was showed in Fig. 2. The change of the trend chart of precision moisture content was showed in Fig. 3. The change of the trend chart of dynamic viscosity test was showed in Fig. 4. There were many results which could be seen from Figs. 2-4. The PQ value was relatively high at the 9th and 10th test which were both 45. The moisture content was the highest at the 3rd test which was 1588 μg/g. The results of the kinematic viscosity test at 40 ℃ did not change significantly with a minimum measurement of 61 mm2/s and a maximum measurement of 68 mm2/s.
Fig. 2The chart of PQ value trend

Fig. 3The chart of precision moisture content trend

Fig. 4The chart of kinematic viscosity trend at 40 ℃

2.2. Calculation method for correlation degree
Correlation analysis generally includes following calculations and steps: (a) transformation of raw data, (b) calculation of correlation coefficients, (c) find the degree of closure, (d) arrange the closure sequence.
There are time series:
That is , ,..., (1, 2,..., ). is the length of each sequence, which is the number of data. These m sequences represent m factors (variables). Set another time series (1, 2,..., ), which is called the mother sequence. The time series mentioned above are called sub sequences.
2.2.1. Data transformation: eliminate dimensions (or units) from the original data which were necessary and convert it into a comparable data sequence
A. Mean transformation: First, calculate the average of each sequence separately, and then remove the corresponding original data in each column from the average to obtain a new data column, which is the mean sequence. Each value is greater than 0, and most of them are close to 1. The sequence curves are interrelated intersection.
B. Initial value transformation: Use the first data of the same sequence to remove the subsequent raw data. A new multiple sequence is obtained, which is the initial value sequence. Each value is greater than 0 and the sequence has a common starting point.
C. Standardization transformation: First, calculate the mean and standard deviation of each sequence separately. Then subtract the mean from each original data and divide by the standard deviation. The new data sequence obtained in this way is the standardized sequence. The sequence mean is 0 and the variance is 1.
2.2.2. Calculate correlation coefficient
The mother sequence after data transformation is denoted as . The child sequence is denoted as . Therefore, the correlation coefficient between the mother sequence and the child sequence can be calculated by Eq. (1) at time :
where: represents the absolute difference between two comparison sequences at time , i.e. . and respectively represent the maximum and minimum absolute differences of all comparison sequences at each time step. Because of the comparison sequences intersect, is generally taken. is called the resolution coefficient, which means to weaken the distortion caused by the maximum absolute difference value being too large. It also can improve the significance of the difference between correlation coefficients. and it is generally taken as 0.1-0.5.
The correlation coefficient reflects the degree of closeness between two compared sequences at a certain moment. For example at moment, 1, while at moment, the correlation coefficient is at its minimum value. Therefore, the range of correlation coefficient is 01.
2.2.3. Find correlation degree
From the above, it can be seen that correlation analysis is actually a geometric comparison of time series data. If two sequences coincide at each time point, that means the correlation number is equal to 1, then the correlation between the two sequences must also be equal to 1. On the other hand, the two comparison sequences cannot be perpendicular at any time, so the correlation coefficient is greater than 0 and the correlation degree is also greater than 0. Therefore, the correlation between the two sequences is calculated as the average of the correlation coefficients at each time point of the two comparison sequences. That is Eq. (2):
where represents the degree of association between subsoquence and mother sequence 0. represents the length of the comparison sequence (the number of data).
2.2.4. Sorting sequence
By arranging the correlation between m subsequences and the same mother sequence in order of magnitude, the maximum correlation factor can be determined with the parent sequence.
3. Results and analysis
3.1. Calculation of correlation degree of lubricating oil indicators for sliding bearings
The experimental data of sliding bearing lubricating oil were averaged firstly. The resolution coefficient is selected as to calculate the correlation between indicators in various row in the Table 2. Fe element data were used as the mother sequence. Cr, Pb, Cu, Sn, Ni, PQ values, moisture and kinematic viscosity (40 ℃) were calculated and used as sub sequences. The correlation calculation results were shown in Table 2. The second row shows the correlation calculated using Cr as the mother sequence and Pb, Cu, Sn, Ni, PQ values, moisture, and 40C kinematic viscosity as sub sequences, and so on.
Table 2The correlation degree of lubricating oil indicators for sliding bearings
Correlation degree | Cr | Pb | Cu | Sn | Ni | PQ | Moisture | Viscosity |
Fe | 0.6977 | 0.6824 | 0.7377 | 0.8640 | 0.7168 | 0.7831 | 0.6492 | 0.4723 |
Cr | 0.8722 | 0.7856 | 0.7200 | 0.8102 | 0.6020 | 0.8781 | 0.3837 | |
Pb | 0.8383 | 0.7036 | 0.7797 | 0.5958 | 0.8815 | 0.3941 | ||
Cu | 0.7460 | 0.7803 | 0.6393 | 0.7723 | 0.4331 | |||
Sn | 0.7417 | 0.7575 | 0.6726 | 0.4563 | ||||
Ni | 0.6295 | 0.8199 | 0.4213 | |||||
PQ | 0.6179 | 0.5581 | ||||||
Moisture | 0.3958 |
3.2. Analysis
The data in column 8 of Table 2 are all greater than those in column 9 of the corresponding row. It was indicated that the concentration of elements such as Fe, Cr, Pb, Cu, Sn, Ni was more closely related to moisture than to kinematic viscosity in the lubricating oil of sliding bearings. This indicated that the moisture content in the lubricating oil of sliding bearings was closely related to wear. Controlling the moisture content was an important measure to prevent wear of sliding bearings.
The correlation between Fe and Sn, PQ values was high. The correlation between Sn and PQ values was high too in the tested sliding bearing lubricating oil. These indicated that they increased simultaneously during bearing wear. Fe was the substrate element of the rotating shaft. Sn was the substrate element of the bearing shell. PQ value was the quantity of ferromagnetic abrasive particles such as Fe and the friction. The rub between the shaft and the bearings could cause significant wear on the substrate which was consistent with objective reality.
The correlation between PQ value and kinematic viscosity was the highest among all wear indicators and kinematic viscosity. It was indicated that the amount of ferromagnetic abrasive particles had a significant impact on kinematic viscosity in the lubricating oil.
4. Conclusions
In the daily monitoring of sliding bearing lubricating oil, the following conclusions can be drawn: (1) Controlling the moisture content in the shaft lubricating oil was an important measure to prevent sliding bearing wear; (2) Special attention should be paid to the growth of wear indicators such as Fe, Sn, and PQ values. If there were abnormalities, timely measures should be taken to prevent excessive wear of sliding bearings; (3) We need to strengthen the monitoring of PQ values, and take timely measures to control any abnormalities to prevent accelerated deterioration of lubricating oil.
References
-
P. Qin, B. Yan, and R. Y. Song, “The condition monitoring of contacting fault of plain bearings,” China Mechanical Engineering, Vol. 13, No. 8, pp. 689–691, 2002.
-
A. R. Chen, “Fault Diagnosis technology for slide bearing based on Fourier Transform,” Bearing, No. 6, pp. 56–60, 2009.
-
Z. Y. Wang, Z. Li, and J. H. Wu, “On grey choice of the best measurement plan,” The Journal of Grey System, Vol. 9, No. 4, pp. 371–378, 1997.
-
S. Swaminathan and C. Smidts, “The event sequence diagram framework for dynamic probabilistic risk assessment,” Reliability Engineering and System Safety, Vol. 63, No. 1, pp. 73–90, 1999.
-
S. Al-Sharhan, F. Karray, W. Gueaieb, and O. Basir, “Fuzzy entropy: a brief survey,” in 10th Annual IEEE Conference on Fuzzy Systems, Vol. 2, pp. 1135–1139, Oct. 2024, https://doi.org/10.1109/fuzz.2001.1008855
-
Z. Pawlak, “Rough set theory and its application to data analysis,” Cybernetics and Systems, Vol. 29, No. 7, pp. 661–688, 1998.
-
Z. Y. Wang, P. Qin, and Y. S. Gao, “Grey evaluation of measurement uncertainty,” The Journal of Grey System, Vol. 11, No. 4, pp. 347–352, 1999.
-
T. L. Tien, “A research on the deterministic grey dynamic model with multiple inputs DGDMMI (1, 1, 1),” Applied Mathematics and Computation, Vol. 139, pp. 401–416, 2003.
-
D. Q. Truong and K. K. Ahn, “An accurate signal estimator using a novel smart adaptive grey model SAGM (1, 1),” Expert System with Applications, Vol. 39, No. 9, pp. 7611–7620, 2012.
-
Y. Yang and D. Xue, “Continous fractional-order grey model and electricity prediction research based on the observation error feedback,” Energy, Vol. 115, pp. 722–733, 2016.
-
C.-H. Lin, C.-H. Wu, and P.-Z. Huang, “Grey clustering analysis for incipient fault diagnosis in oil-immersed transformers,” Expert Systems with Applications, Vol. 36, No. 2, pp. 1371–1379, Mar. 2009, https://doi.org/10.1016/j.eswa.2007.11.019
About this article
The authors have not disclosed any funding.
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
The authors declare that they have no conflict of interest.