General Log-Linear Weibull Model Combining Vibration and Temperature Characteristics for Remaining Useful Life Prediction of Rolling Element Bearings

In industrial applications, the vibration and temperature measurements of rolling element bearings are known as two popular condition monitoring methods. *e previously published method for remaining useful life (RUL) prediction has been limited to using the vibration signal. However, a single signal source cannot fully reflect the degradation trend of bearings, influencing the RUL prediction precision. In this paper, a novel general log-linear Weibull (GLL-Weibull) model is constructed by considering vibration and temperature conditionmonitoring signals to estimate themodel parameters. During the feature extraction stage, the relative root mean square (RRMS) is derived from the monitored vibration signal, and the relative temperature trend value is extracted from the monitored temperature signal to eliminate individual differences in bearings and random signal fluctuations. *en, a fuzzy operator is introduced to describe the degree of an “overheated bearing” and “excessive bearing vibrations.” During the RUL prediction stage, both the extracted vibration and temperature characteristics are used to create the GLL-Weibull model. *e best parameters are attained by employing the maximum likelihood estimation approach. *e algorithm performance is checked with criteria like the root mean square error (RMSE) and the mean absolute percentage error (MAPE). *e effectiveness and superiority of the presented approach are validated by two real-life prognosis cases. According to the experimental results, the presented approach provides superior prediction precision and lower computational cost than other approaches for bearings under constant or variable operating conditions.


Introduction
Rolling element bearings are known as important parts in various machines, and bearing failure can lead to considerable losses in production and human casualties. erefore, a higher reliability and readiness state for bearings is urgently required, especially for complex and high-maintenance-cost machines such as wind turbines [1]. In Industry 4.0, industrial cyberphysical systems (ICPSs) have attracted increasing attention both in industry and academia [2][3][4]. With the development of new intelligent sensors, it is possible to replace the condition monitoring maintenance of large and complex machines with preventive maintenance. To ensure machinery safety, prognostics and health management, including the prediction of performance degradation, have drawn much attention. Prognostics and health management can be utilized to conduct timely maintenance and attain timely RUL predictions and intelligent accident prevention [5,6]. e effective prediction of RUL must overcome two fundamental deficiencies: developing a useful model and choosing the parameters that reflect the efficiency degradation procedure. Statistical model-based methods estimate the machinery RUL by developing empirical knowledgebased statistical models and employing a probability density function (PDF) to present the observation-dependent prediction results. In the mentioned methods, RUL prediction models are established by fitting the acquired measurements into stochastic coefficient or stochastic process models by using a probabilistic approach without utilizing any physical concepts. erefore, statistical model-based methods can be utilized to model the degradation procedure uncertainty and its impact on RUL prediction. Among these approaches, the proportional hazards model (PHM) was first presented by Cox [7] in 1972. In the mentioned model, it was assumed that the system hazard rate consists of two multiplicative functions, including baseline hazard and covariate functions. PHM has no special requirements for data distribution and can analyze failure data and suspension data at the same time. ey have been extensively studied and applied in the area of reliability assessment and RUL predicting of mechanical parts. In the literature, Kundu et al. [8] presented a Weibull accelerated failure time regression model for bearings working under multiple operating conditions for the URL prediction. Regattieri et al. [9] presented a novel generalized log-linear Weibull (GLL-Weibull) distributionbased method to estimate the gear life. Kundu et al. [10] combined the clustering output with a change point identification technique to develop RUL prediction models using a general log-linear Weibull (GLL-Weibull) distribution.
Most of the current research is limited to the statistical analysis of a single signal source. Using signal processing techniques, the condition monitoring information is obtained as the input of the prediction model. Vibration signals have been extensively utilized for bearing fault detection and prognosis because they are easy to collect and sensitive to the degradation performance of rolling bearings [11][12][13][14]. For instance, Wang et al. [15] extracted time-frequency domain degradation features to construct an enhanced Wiener process model to predict the RUL by utilizing the vibration information of the bearing and predicting the rolling element bearings RUL. However, the vibration signal easily interferes with rolling element bearings working in the presence of variable load and speed situations, and the prediction accuracy is reduced due to the poor signal-tonoise ratio. erefore, other signals should be used. Recently, the temperature signal has attracted the attention of some scholars. For example, Hu et al. [16] used temperature characteristics to realize the RUL prediction of a wind turbine generator rear bearing. Although the temperature signal has good anti-interference performance, its linear change trend cannot fully demonstrate the nonlinear degradation procedure of the bearing, and the available degradation information can be limited.
In light of the above problems, we propose a novel GLL-Weibull model for the RUL prediction of bearings that considers both vibration and temperature characteristics by taking full advantage of the high sensitivity of vibration characteristics and the powerful anti-interference capability of temperature characteristics (temperature can be utilized to detect the bearings decay caused by the thermal inertia). e extracted features can contain more effective degradation information; thus, the performance of the model is greatly improved. According to the experimental results, the presented approach provides a superior prediction precision and lower computational cost than other approaches for rolling bearings under constant or variable operating conditions. e rest of the current study is organized as follows: Section 2 first presents the GLL-Weibull model theory. en, the methods for multicovariate selection are presented, which include the extraction of vibration features from vibration signals and the extraction of temperature trends from temperature signals. Next, the procedure for the RUL estimation is described in detail in Section 3. e performance of the constructed model is evaluated through a PRONOSTIA data set and a real wind turbine high-speed shaft (HSS) bearing data set for 50 days in Section 4. In the end, conclusions are provided in Section 5.

GLL-Weibull Prediction Model
e evolution law of bearing failure probability basically follows the classic "bathtub curve" with service time shown in Figure 1, which can be divided into three state stages: early failure stage, occasional failure stage, and fatigue failure stage. e early failure stage is mainly caused by the initial design, manufacture, and transportation of the equipment and defects resulting from the installation. Before service, each bearing is subjected to a series of industrial tests according to predefined rules, which can effectively suppress the probability of early failure for the bearing in operation. e occasional failure stage and fatigue failure stage occur during the normal bearing service life and are affected by the combined actions of thermal, mechanical, and chemical stresses as well as environmental factors. During these two stages, defects are formed in the bearing, or fatigue gradually deteriorates the bearing, which leads to an increase in the failure rate.

Fundamental eory.
In accordance with the proportional hazards model (PHM) proposed by Cox, the system failure rate depends on its operation time and the operation covariates. Mathematically, the widely accepted cumulative density function (CDF) of the 2-parameter Weibull distribution can be defined as follows: where t denotes the service time of the wind turbine bearing, i.e., the primary ageing parameter, η describes the characteristic life parameter that indicates when the failure probability is equal to 63%, and β represents the shape parameter that is used to describe different failure stages. e shape parameter defines failure dispersion. β < 1: e failure component risk reduces by increasing the ageing parameter value, which belongs to the early failure stage. β � 1: e failure component risk is not changed while increasing the ageing parameter value, which belongs to the occasional failure stage.

2
Shock and Vibration β > 1: e failure component risk grows by increasing the ageing parameter value, which belongs to the fatigue failure stage.
Although the failure probability will change with time, various stimulus parameters will also accelerate the expansion of failure rate improvement. e mentioned parameters include the exerted load, operating velocity, pressure, vibration, temperature, torque, etc. Here, these parameters are called covariates. Covariates can be categorized into external and internal covariates. Factors like speed, load, and pressure are named as external covariates. Measurement data from various sensors are used to monitor health status, including vibration, temperature, and torque signals.
e outputs of these sensors are termed internal covariates. For bearing problems, when external operating conditions (load, speed) change, the vibration and temperature feature indicators will also be affected. In other words, when external covariates are uncertain or not easy to identify analytically and experimentally, the utilization of various measurable parameters (i.e., vibration and temperature signal) can guarantee the reduction of various uncertainties that originate from environmental conditions, loading levels, and measurement equipment. erefore, this paper introduces the internal covariates into the Weibull distribution model to establish a general log-linear Weibull (GLL-Weibull) model, which can also be called the Weibull accelerated failure model. e CDF for the GLL-Weibull model is given by the following: where the intercept is denoted by a 0 , the i−th covariate (i � 1, 2, . . . , k), is described by x i , and a i represents the influence coefficient of the covariates, which can also be called the regression coefficient.
Letting η ′ � exp(a 0 + k i�1 a i x i ), the CDF for the GLL-Weibull model can be rewritten as follows: e hazard rate function for the GLL-Weibull model is as follows: According to the principle of reliability analysis [17], the CDF reliability is also referred to as unreliability. Consider that the distribution reliability function can be easily obtained by subtracting the CDF from one. us, the reliability function of the GLL-Weibull model can be written as follows: Similarly, F(t) is defined as the probability that failure happens at time t or the cumulative distribution function of the failure distribution. Accordingly, the density function of the failure probability is described as follows: where i is the failure time, j is the suspension time, n describes the number of failure samples, and m indicates the number of suspension samples. Usually, the log-likelihood function is more tractable than the likelihood function, and the simplified form of the log-likelihood function for the GLL-Weibull model can be expressed as follows: Equation (8) can also be written as follows: By setting the first partial derivatives of β, a 0 , a in equation (9) equal to zero and then adopting the Newton-Raphson [18] iterative method for numerical calculation, the optimal estimation of each parameter can be obtained as β, a 0 , a. e GLL-Weibull model is a completely parameterized model that combines historical condition monitoring data to estimate the RUL at any time. In the presented method, multiple covariates (through measured vibration signals and temperature signals) have been considered to enhance the modeling accuracy and prediction precision.

Feature Extraction Method.
In condition monitoring, vibration and temperature features are two important indicators that can reflect the dynamic behavior of bearing degradation. Vibration monitoring is a useful tool for the primary prediction and identification of the bearing failures. On the other hand, bearing temperature can provide information on the upcoming deterioration procedure from extreme mechanical friction.

Vibration Feature Extraction.
e root mean square (RMS) describes the effective amount (magnitude) of the vibration signal. It can be obtained by calculating the square root of the mean of the sum of the squares of the signal samples using the following equation: where x(i) denotes the ith sample of the primary signal, while the number of samples is represented by n. Usually, the vibration signals' RMS values are utilized to verify the whole status of the mechanical parts. e reason is that the whole vibration extent generally grows with the service life due to deterioration of the parts like gearbox and bearing. RMS is used as a widely accepted characteristic index for bearing fault diagnosis.
However, even under similar operating situations, the RMS values of various bearings in a stable working stage differ from the others. erefore, the relative root mean square (RRMS) is proposed as a measurement in this paper, which can be computed by the following equation: where x base represents the baseline value. Note that the average mean square root within the relatively flat occasional failure stage of a bearing can be taken as x base . Single differences in bearings do not influence the RRMS as a simple index. e RRMS is sensitive to the primary damage and continuously increases with the damage growth. en, the smoothing method with 7 points 3 times is adopted to reduce the impact of stochastic vibration, which can be described in the following equation: where x is the signal before smoothing, x MA is the signal after smoothing, and k � 1, 2, . . . , N denotes the number of signal samples.

Temperature Trend Analysis.
Temperature characteristics can reflect the degradation process of wind turbine bearings caused by many external uncertain factors. To extract the temperature trend of the degradation period from the temperature monitoring data, (t 1 , r 1 ), (t 2 , r 2 ), . . ., (t m , r m ) denoting the bearing speed time series and (t 1 , c 1 ), (t 2 , c 2 ), . . ., (t m , c m ) denoting the bearing temperature time series are initially obtained. en, the relative temperature data (t 1 , c 1 ′ ), (t 2 , c 2 ′ ), . . ., (t m , c m ′ ) under different speed operating conditions can be calculated as follows: Furthermore, the moving average method is used to eliminate random fluctuations, and the temperature trend data of the bearing degradation period is obtained, which is expressed in the following equation: where x represents the temperature trend data at the kth time instant and p stands for the number of adjacent data points calculating the mean value. Based on equation (14), the time series of the temperature trend can be described as (t 1 , x 1 ), (t 2 , x 2 ), . . ., (t m , x m ), which will gradually increase with the deterioration of the bearing.

Fuzzification Processing.
During practical testing, although the vibration and temperature trend values can reflect the evolution of the bearing from the degradation state to final failure, the sensitivity of the two signals to different types of faults is inconsistent. For example, the vibration and temperature signals are susceptible to position and coupling faults, respectively. Each signal will grow faster under the corresponding sensitive faults, so the change rate of the two signals may be different under the same fault.
By using fuzzy logic, it is possible to process the indications from temperature measurements in parallel with indications from vibration measurements, which can be accomplished such that the signals from the temperature and vibration recording will be processed with a disjunctiveprobabilistic fuzzy operator. In a fuzzy set, for each variable t, a fuzzy operator is defined, which may have values in the range between 0 and 1. Expressions of fuzzy operators for vibration and temperature signals are given by the set of equations (15) and (16), respectively: where Vib and Tem are the fuzzy operators for vibration and temperature indicators, respectively, v(t) and x(t) are the indicator values corresponding to each variable t, the vibration limitations are v a and v b , and the temperature limitations are θ a and θ b . rough fuzzification processing, real data are transformed into fuzzy operators. In this paper, the fuzzy subsets of vibration and temperature correspond to three various scenarios, as shown in Figure 2.
(i) Normal state: All the parameters have under control amounts, while the bearing status is accurate. e fuzzy operator is smaller than 0.5. (ii) Early warning: Bearings have begun to degrade, and defects will continue to expand if preventive maintenance is not available. e fuzzy operator ranges from 0.5 to 0.8. (iii) Near failure: Recognition and immediate intervention are required to set proper values for the status parameters. In other words, the area belongs to the fuzzy set named "excessive vibrations" or "overheated bearing." e fuzzy operator is greater than 0.8.
Based on the vibration and temperature fuzzy operators proposed in this paper, the reliability function at time t for the bearing issue with vibration and temperature coordination is described as follows: where the intercept is represented by a 0 , β describes the distribution shape parameter and a 1 and a 2 are the regression coefficients for vibration and temperature fuzzy operators, respectively. en, the bearing RUL is determined using the following conditional reliability: where t 0 denotes the present time, t indicates the bearing failure time, while R(t | t 0 ) represents the conditional reliability, which is a component or system that has survived to the given time t 0 and works without failure until time t. At the current time, the bearing reliability can be obtained from equation (17). e time to failure of the bearing can be approximately obtained by extrapolating the present condition monitoring to the critical conditional reliability presented in (18) in which its value is determined by the user, and the time t is the RUL of the bearing at t 0 .

Performance
where y i and y i describe the real and predicted RUL values, respectively. n denotes the number of data points, and i is the predicted point sequence number. Smaller values of these two criteria indicate higher prediction.

Remaining Useful Life Prediction
Procedure for Rolling Bearings e prognostic procedure for determining the RUL is shown in Figure 3. e proposed method can be divided into offline and online stages. During the offline stage, by analyzing the data of historically failed bearings, the RRMS and the relative temperature value are derived from the raw vibration and temperature signals as some features, respectively. e smoothing process is performed so that the obtained features better reflect bearing degradation. en, according to the allowable variation range of vibration and temperature, the two feature values are fuzzified, and the fuzzy operators are used as two internal covariates to estimate the values of unknown parameters of the GLL-Weibull model. e model performance is evaluated using specific performance criteria. During the online stage, the GLL-Weibull model coefficients attained from model training are employed to estimate the RUL. For a service bearing, the RUL estimation is performed based on the current time. en, the maintenance decision is made by the user.  [19]. e experimental equipment is shown in Figure 4.

Experimental Validation and Analysis
is platform provides practical data describing the normal degradation process of rolling bearings from run to failure under different working situations. Two accelerometers and a thermocouple are installed to pick up the horizontal and vertical vibration signals and the temperature signal. e vibration and temperature signals are acquired every 10 s and 60 s at sampling rates of 25.6 kHz and 10 Hz, respectively. e experiments are finished when the amplitude of the monitoring vibration data becomes greater than 20 g. Hence, the vibration data from run to failure is chosen as a threshold for model training and validation.
e observed time to failure information for different bearings under specific working conditions is presented in Table 1. Some studies have shown that some bearings in the data set have an abnormal degradation trend due to improper installation or operation [20][21][22]. Based on overall vibration and temperature history, bearings 4, 5, 6, and 7 from operating condition (1) (OC1), bearing 4 from OC2, and bearing 3 from OC3 have been considered during initial modeling.

Experimental Results.
e presented prediction approach includes two parts. e first one processed the vibration and temperature signals of the bearing's degradation status, including feature extraction, smoothing, and fuzzification, to observe the trend behavior of the two signals.   Shock and Vibration First, according to the content in Section 2, the RRMS amounts of the target bearing vibration signal and the relative temperature amounts of the target bearing temperature signal are calculated. At the same time, the smoothing process is performed, as small fluctuation regions in the feature can influence the stability of the prognosis model. Figure 5 shows the corresponding feature after the smoothing process for bearing 4 from OC1.

Shock and Vibration
In addition, the extracted feature must appropriately correlate with the age and health of the bearing. e mentioned correlation level can also affect the further stages and even the ultimate prognosis results. e fuzzy operator indicates the feature correlation to time and varies among zero and one and tells us when "overheated bearing" and "excessive bearing vibrations" appear. Figure 6 shows the corresponding fuzzy operator for 6 bearings from different operation conditions. As can be seen from Figure 6, the vibration characteristics of all bearings change slightly in the early stage and sharply increase in the later stage. It is complicated to detect the early degradation point of the bearing using the vibration characteristics alone. A single vibration signal cannot fully reflect the degradation trend of the bearing. erefore, other signals should be used. In this paper, temperature characteristics are introduced to assist the diagnosis of early bearing degradation. Rapid temperature change or continuous high temperature can cause bearing thermal strain, which accelerates bearing failure. erefore, the point when the rising rate of the temperature fuzzy operator suddenly accelerates (i.e., bearings 1-4 and 1-6) or the temperature fuzzy operator fluctuates near 0.8 (i.e., bearings 1-5, 1-7, 2-4, 3-3) is selected as the starting point of bearing degradation to prepare for the next RUL prediction. e second part of the presented approach describes the GLL-Weibull model establishment and RUL prediction. Six bearing failure records are utilized for model training and validation, as previously illustrated in Section 2.1. e leaveone-out cross-validation technique is utilized to check the model efficiency. Five bearing failure records are employed to estimate the model parameters. e model's effectiveness is evaluated with the sixth bearing, which was not assumed within estimating model parameters. For example, bearings 1-5, 1-6, 1-7, 2-4, and 3-3 are utilized to estimate the model parameters for the RUL prediction of bearing 1-1. Now, the same procedure is performed for remaining bearings (see Table 2). Table 3 presents the corresponding values of the model parameters attained from various training data sets. e calculated values for the bearing shape coefficients belong to [1.67, 1.68]. e intercept amount is in a range from 6.02 to 6.15, the coefficient for the vibration fuzzy operator lies in a range from 0.063 to 0.099, and the coefficient range for the temperature fuzzy operator lies between 0.0024 and 0.0057. Considering several training data sets, a small change in the model coefficient values indicates model consistency. Generally, when the bearings approach the failure status (i.e., failure state degradation or the wear increase), both the vibration and temperature fuzzy operators increase. ese coefficients usually describe the impact of covariates on the failure procedure. e coefficient values of the vibration and temperature fuzzy operators are positive. It can be  concluded that an increase in the amount of vibration and fuzzy temperature operators leads to degradation or wear in the bearing. Besides, the corresponding value for the vibration fuzzy operator is higher compared with the fuzzy temperature operator, demonstrating a powerful impact of vibration on the bearing life. e RUL predicting results for the GLL-Weibull model are presented in Figure 7. e mentioned results correspond to the 95% critical failure probability, i.e., (1 − R(t + t 0 )). e user can select various values for the critical failure probability, and in the current state, the optimal RUL predicting results are achieved with the 95% failure probability. Twenty sets of prediction points were obtained at each 5% time range from the beginning of the degradation phase to its end.
Moreover, α − λ performance [23] is utilized to verify the predicting performance of the presented approach, where the parameter α is adopted to draw upper and lower bounds for error in the RUL estimation procedure, and the parameter λ describes the relative time distance of a determined point within working from the end of life. In the current work, the efficiency of the presented approach with α � 20% is verified. As seen in Figure 7, the prediction accuracy can be guaranteed around α � 20% for almost all instances (the experimental results are achieved by computing the mean values). Accordingly, after training the GLL-Weibull by historical data, the RUL value can also be predicted for different values of the load and speed and similar bearing types.

Comparison Analysis.
To further evaluate the efficiency and advantage of the presented approach, it is compared with previous approaches for the RUL estimating of bearings, like a linear regression model (LRM) [24], a support vector machine (SVM) [25], a Kalman filter (KF) [26], and a convolutional neural network (CNN) [27]. ( e similar data set is utilized in these approaches-the PRO-NOSTIA data set.) Table 4 shows the performance metrics of the different methodologies. In Table 4, it is observed that the LRM has the largest prediction error because the LRM is not suitable for describing nonlinear and nonstationary processes. An LRM trained on all the samples of the entire history may not be able to model the life patterns well. e SVM provides better prediction for RUL with large data sets, but improper parameters significantly affect the prediction results due to the lack of any standard approach for choosing the kernel function. e present and future failure states can be predicted through the KF. However, a system dynamics knowledge is required to construct straightforward model relations that depend on domain knowledge and a large data set. Although CNN has a higher and longer prediction horizon, finding the optimal network framework is complicated. Moreover, the neural network is unstable and provides various results for various numbers of runs. e proposed method can obtain more accurate RUL results with fewer condition monitoring histories from each operating condition, and the computation time is low.

Testing Data Description.
To validate the efficiency of the presented RUL predicting approach, a high-speed shaft bearing of a 2 MW wind turbine at a wind farm in Northeast China is taken as another example for analysis. e running speed range of the bearing is 1200 rpm∼1800 rpm. In March 2019, a wind turbine stopped working due to a high temperature exceeding the maximum permissible limit of 95°C. e wind turbine gearbox suffered from excessive vibration and a continuous rise in temperature. e root cause of the failure was the high-speed shaft (HSS) bearing in the  Training  Training  Training  Training  Training  Training  Test  Training  Training  Training  Training  Training  Training  Test  Training  Training  Training  Training  Training  Training  Test  Training  Training  Training  Training  Training  Training  Test  Training  Training  Training  Training  Training Training Test  gearbox, and the vibration data were collected for 50 consecutive days by an acceleration sensor mounted on the horizontal direction of the bearing seat. e sampling rate of the bearing vibration signal was 97 656 Hz. e time domain signal of the bearing life is presented in Figure 8. At the same time, viewing the temperature and speed data recorded in the supervisory control and data acquisition (SCADA) system (data sampled in 10 min intervals), it was found that the bearing temperature showed a total temperature trend rising sharply as a fluctuation until failure (see Figure 9).
According to the above data, the RUL prediction process shown in Section 3 is used to verify and analyze the HSS bearing of the wind turbine. Figures 10 and 11 show the corresponding features after the smoothing process and the fuzzy operator for the HSS bearing from run to failure, respectively. As seen from Figures 10 and 11, the vibration feature is sensitive to early bearing faults but has larger fluctuations. On the other hand, the temperature feature demonstrates a steady upward trend with bearing degradation, but the detection of weak faults is not effective. In summary, the vibration signal combined with the temperature signal can fully reflect the dynamic degradation of the HSS bearing. e proposed approach is an effective and reliable technique for monitoring the health status and predicting the RUL of bearings.

Results and Discussion.
To estimate the GLL-Weibull model parameters, the 50-day monitoring data are divided into 7 monitoring periods. e time of every monitoring period is 168 hours (h), including 1008 sample points (10 min per sample interval). According to the discussion in Section 2, the four parameters of each monitoring period are estimated using equation (9), as shown in Table 5. e probability density function (PDF) distribution for the RUL of the HSS bearing is obtained from the estimated parameters presented in Table 5. After 168 h monitoring interval, the PDF distribution of the RUL is determined (see Figure 12). According to Figure 12, the maximum PDF value Fis attained at the corresponding time of the predicted RUL amount, and the RUL time of the HSS bearing after the initial monitoring time range is 2192 h for about 91 days. Likewise, a similar procedure is employed to attain the RUL conditions under the remaining monitoring periods. e predicted and real RUL amounts of each monitoring period for the HHS bearing are illustrated in Figure 13, and their two error types are also listed as shown in Table 6. Figure 13 demonstrates that the efficiency of the HSS bearing degenerates considerably as the monitoring period interval increases. Narrowing the PDF width makes the RUL predicted results closer to their corresponding real ones. Figure 13 shows that the probability density distribution changes from wide to narrow with the deterioration of the bearing and the increase of the monitoring data, and the predicting result of the RUL gradually converges to the real one. In Table 6, two different errors are calculated with equations (19) and (20). e degradation procedure of the HSS bearing includes initial, middle, and final phases. It can be observed that the two errors at the initial stage are larger than those at the middle and final phases. e errors decrease as the monitoring periods increase, and there are some small fluctuations at the middle stage, as shown in Figure 14.
e reasons are analyzed in accordance with Table 6, Figures 13 and 14. At the initial stage of RUL prediction, there are fewer data points available and more uncertainties in the parameter estimation. Compared with the real RUL values, the predicted ones are large. us, the two errors are also large. At the middle stage, the vibration and temperature fuzzy operator enters a fluctuating state with the development of the fault. e parameters vary over a small range, so the prediction results also fluctuate. At the final stage, the parameters of the model gradually converge as the monitoring data increase, and the predicted results are more accurate.
In summary, this study demonstrates that the real-time RUL of HSS bearings can be attained by the presented approach using the recorded data of vibration and temperature.

Conclusions
e present paper proposed a novel general log-linear Weibull (GLL-Weibull) model combining vibration and temperature characteristics for bearing RUL prediction. According to the attained results, the upcoming results can be derived.
(1) e GLL-Weibull model is a multicovariate-based reliability model that can be used for the reliable estimation of mechanical components. During model development, the effects of vibration and temperature monitoring data are considered. (2) By extracting the RRMS of the vibration signal, the influence of individual differences on different degradation stages can be eliminated. (3) By extracting the relative temperature feature, it can be observed that the relative temperature increases the trend of the bearing from run to failure.    (4) Due to asynchronous problems in vibration and temperature characteristics, the fuzzy operator is a reliable tool for describing the degree of faulty bearing degradation. (5) e effectiveness and superiority of the presented approach are evaluated using two different bearing data sets. Due to its robustness, this approach can be applied to bearings under constant or variable working situations.
Although the presented approach is evaluated with bearings examined in different cases, its usefulness and progress under strong interference working situations will be verified in the future. In addition, a mixed RUL prediction model after bearing maintenance will be explored in our next research project.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare no conflicts of interest.