Next Article in Journal
Three-Dimensional Monitoring of Plant Structural Parameters and Chlorophyll Distribution
Next Article in Special Issue
A Mission Reliability-Driven Manufacturing System Health State Evaluation Method Based on Fusion of Operational Data
Previous Article in Journal
Experimentation Management in the Co-Created Smart-City: Incentivization and Citizen Engagement
Previous Article in Special Issue
Machine Learning-Based Sensor Data Modeling Methods for Power Transformer PHM
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differential Equation-Based Prediction Model for Early Change Detection in Transient Running Status †

1
Key Laboratory of High-efficiency and Clean Mechanical Manufacture of MOE, National Demonstration Center for Experimental Mechanical Engineering Education, School of Mechanical Engineering, Shandong University, Jinan 250061, China
2
School of Mechatronics Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
*
Author to whom correspondence should be addressed.
This paper is an extended version of paper published in Wen, X.; Lu, G.; Yan, P. Improving Structural Change Detection using a Differential Equation-based Prediction Model for Condition Monitoring of Rotating Machines. In Proceedings of the 2018 IEEE International Conference on Prognostics and Health Management (ICPHM), Seattle, WA, USA, 11–13 June 2018.
Sensors 2019, 19(2), 412; https://doi.org/10.3390/s19020412
Submission received: 23 December 2018 / Revised: 16 January 2019 / Accepted: 18 January 2019 / Published: 20 January 2019
(This article belongs to the Special Issue Sensors for Prognostics and Health Management)

Abstract

:
Early detection of changes in transient running status from sensor signals attracts increasing attention in modern industries. To achieve this end, this paper presents a new differential equation-based prediction model that can realize one-step-ahead prediction of machine status. Together with this model, an analysis of continuous monitoring of condition signal by means of a null hypothesis testing is presented to inspect/diagnose whether an abnormal status change occurs or not during successive machine operations. The detection operation is executed periodically and continuously, such that the machine running status can be monitored with an online and real-time manner. The effectiveness of the proposed method is demonstrated using three representative real-engineering applications: external loading status monitoring, bearing health status monitoring and speed condition monitoring. The method is also compared with those benchmark methods reported in the literature. From the results, the proposed method demonstrates significant improvements over others, which suggests its superiority and great potentials in real applications.

1. Introduction

Real-world systems are seldom in steady status, and they almost always operate in transient conditions that are varying over time. Detection of change(s) in running status at an early stage enables to find and locate abnormal and unexpected/undesired system behavior(s) in its successive operations. With this, corrective schedule could be made to prevent potential operation failures and ensure equipment and product reliability, safety, quality, productivity, etc. [1,2,3,4]. The problem of running status monitoring may arise in process monitoring and control applications where the system needs to make a response and/or take appropriate actions as soon as possible after a change of machine status occurs [5,6,7,8].
In industrial and manufacturing areas existing approaches of change detection are mainly data-driven methods [7,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25]. They normally use a detailed mathematical model evaluating dynamic running status, parameters of which are estimated and updated with relevant condition monitoring (CM) signals that are collected form considered machines. Once the dynamic machine behaviors are modeled, it is possible to detect potential status change(s) via analysis of changes of input parameters to the model.
Depending on the adopted models, these data-driven methods can be further classified into three categories: (1) density estimation methods [9,10,11,12], (2) rule-based or case-based methods [13,14,15,16] and (3) time-series analysis based methods [7,17,18,19,20]. The density estimation methods are developed from the availability of input/output histories of machine conditions that were collected from the field or laboratory. They are based on density estimation techniques modeling the past normal machine status, and produce change decision output based on the Bayesian probability theory for newly observed data. Typical Gaussian kernels [9] as well as relatives such as non-Gaussian process [11] and Gaussian mixture model (GMM) [12] can be employed to accomplish the density estimation. In spite of their success and popularity in simulation scenarios, they are, however, not effective when the available condition data is insufficient. The rule-based or case-based methods are expert systems that are driven by data mining, including the Dempster-Shafer (DS) theory [15], fuzzy set theory [16], etc. In these methods, there exist lots of parameter settings that need to be confirmed under human-made supervision, which would limit their practical value in real applications. In comparison, the time-series analysis methods use previous outputs regressed on to itself to provide an estimate of the current output, such that one-step-ahead prediction is allowed for residual error analysis in order to generate final detection result [19,20]. In the context of time-series analysis methods, various auto-regressive (AR) models have been proposed, going from simple linear time-series model to advanced nonlinear regression methods. For example, in [7], a linear AR model was assumed for analytical expression of those linear distributed condition signals, e.g., power consumption signal. Meanwhile, it has been widely recognized and accepted those periodic condition signals, such as vibration and sound signals, are more powerful when describing the machine system [21,22,23,24,25]. It is concluded that the statistical modelling by the way of AR models is relatively simplified in implementation, such that they can be applied to the monitoring and prognosis of industrial systems, especially when real-time and online processing is required.
In this paper, we extend the work of [26] and present a new prediction model with the purpose of monitoring and diagnosing the running status for a considered machine. The baseline of the presented method is an early change detection using a differential equation (DE) based prediction model. This model can describe dynamic properties of a system which change with time quite well and can predict the future status of the system conveniently [27,28]. Moreover, it has been also proven more powerful when discovering of sciential laws for dynamic systems [27]. The DE model has been successfully used in many fields [29,30,31,32]; however, to the best of our knowledge, it is the first time to apply it to machine running status monitoring. In this paper, together with the DE model, an analysis of continuous monitoring of CM signals by means of a null hypothesis testing is proposed to inspect/diagnose whether an abnormal running status change occurs or not in successive machine operations.
Besides of the methodological contributions put forth in the proposed method, we apply the method to three representative applications: external loading status monitoring, bearing health status monitoring, and speed condition monitoring. Three testing data sets are considered in the experiment. The first one is taken from the publicly available data set provided by the Case Western Reserve University (CWRU) [33], the second one is from the publicly available PRONOSTIA dataset provided by FEMTO-ST in France [34,35] and the last one is established based on our experimental setup. The proposed method is also compared with those reported in the literature, including autoregressive-integrated-movingaverage (ARIMA) model [18], root mean square (RMS) and kurtosis [36,37,38]. Experimental results demonstrate its significant improvements over others.
The rest of this paper are organized as follows. Section 2 describes the proposed prediction model. Section 3 gives the hypothesis testing used for decision-making. Section 4 gives the proposed framework for machine running status monitoring, followed by experiments in Section 5. Section 6 finally summarizes this paper with some concluding remarks.

2. Differential Equation-Based Prediction Model

This section first provides an overview of the proposed prediction model, and then discusses techniques carried out in the main steps.

2.1. Overview of the Proposed Model

Given a data stream composing of continuously monitoring condition signal x t , let us assume that it has been denoted by a periodic form x t X n T + v , where T is the period, v [ 1 , T ] is the phase or timing, and n is the cycle index. In this paper, we assume the condition signal is periodically stationary as made in some previous works, e.g., [25,39], thus we can employ a prior periodicity estimation proposed in [18] to confirm the value of T. Then, main steps of the proposed prediction model includes the following:
  • Model formulation, i.e., the proposed model is formulated with the considered CM signals. In this paper, a family of new time series are formed by arranging the original data at the same phase. As such, the model is formulated so as to predict next value of each phase.
  • Parameters estimation, which estimates the parameters of the model. The numerical solution method of differential equations is used to estimate the model’s parameters. The parameters of the model are constantly updated during each data prediction process.
  • Data prediction, i.e., the prediction of next data with the estimated model. The prediction value at each phase can be obtained with the model whose parameters have been estimated successfully.
The above steps are executed continuously for each observed cycles of CM data, allowing for residual error analysis that can quantify the error between the predicted values and the observed ones. The resulting error reflects the amount of temporal fluctuations in the CM data; more specifically, the higher residual error implies a higher probability that a change occurs and vice versa.
In the following, techniques used in the main steps of the model are discussed.

2.2. Proposed Model Description

2.2.1. Model Formulation

For an arrived CM data stream up to the kth cycle, for each phase v, we take i continuous cycles before k to form a new series, i.e., P v = { x ( k i ) T + v , , x ( k 1 ) T + v } , such that the DE model can be formulated for all phases.
Conventionally, the DE model relies on the following system state function f ( · ) to establish the prediction,
d x t d t = f ( x t , β , t ) ,
where x t is the value of differential equation in a series of time stamp t, β is parameter of the differential equation.
In this paper, since we use the past observed data at the same phase for prediction, the above-given new time series P v allows linear characteristics to emerge, based on which, we can obtain the prediction value, i.e., x ^ ( k ) T + v with the P v correspondingly. In the following, for each considered phase v, we re-denote the P v as the series { x t j } , j = 1 , , i for simplicity.
According to the linear assumption made already, the time-series { x t j } will satisfy,
f ( x t , β , t ) = a + b × t ,
where a ( 0 ) and b ( 0 ) are the parameters of differential equations that need to be estimated before prediction.

2.2.2. Data Prediction

Assuming that we have estimated the values of a and b, according to [40], the corresponding time response function of Equation (1) can be written as
x ^ ( k ) T + v = x t i + 1 = 1 2 b × ( t i + 1 2 t i 2 ) + a × ( t i + 1 t i ) + x ¯ t i ,
where t i = t 0 + i × h , t 0 = 0 , t i + 1 is the time stamp of the value x t i + 1 , which can also be regarded as the prediction value x ^ ( k ) T + v , and h is the distance between the two time stamps. We set h as 1 in this article. The x ¯ t i is the initial value of the response function, represented by the average value of the sequence value in order to decrease effect of errors in data collection.

2.2.3. Parameters Estimation

For the purpose of prediction of x ^ ( k ) T + v , the first step is to estimate the optimal parameters of a and b in the DE model given in Equation (2). Because the signal is usually sampled through observation, we re-define the above variables at sampling time t j ( j = 1 , , i ) with subscription j (e.g., x t j = x j ) for convenience in the following of this paper. The differential Equation (1) is calculated by simpsons method [40] by
x j + 2 x j = 1 3 h ( f j + 2 + 4 × f j + 1 + f j ) + R j + 2 ,
in which f j = f ( x j , u , t j ) , j = 1 , , i 2 and R j + 2 is local truncation error [40]. According to Equation (2), Equation (4) can be transformed into
x j + 2 x j = 2 × a + h 3 ( t j + 2 + 4 × t j + 1 + t j ) + R j + 2 ,
which is thus re-written as
x j + 2 x j 2 = a + h 6 ( t j + 2 + 4 × t j + 1 + t j ) + R j + 2 2 .
Then, Equation (6) can be re-expressed as
x j + 2 ( 0 ) = a + x j + 2 ( 1 ) × b + E j + 2 ,
where x j + 2 ( 0 ) = x j + 2 x j 2 , x j + 2 ( 1 ) = ( t j + 2 + 4 × t j + 1 + t j ) h 6 , E j + 2 = R j + 2 2 . When j = 1 i 2 E j + 2 2 has a minimum value for the observed value x j , the estimated parameters a and b can be finally obtained by
( a , b ) T = ( B T B ) 1 B T Y ,
in which
B = 1 1 x 3 ( 1 ) x i ( 1 ) T , Y = x 3 ( 0 ) x i ( 0 ) T .
The predicted value x ^ ( k ) T + v can be obtained by substituting the estimated parameters a and b into the time response Equation (3), where i = 4 is used in the following experiment.

2.3. Residual Error Analysis

The residual error q ( k ) T + v can be defined by the absolute value of the difference between prediction data x ^ ( k ) T + v and the actual monitoring data x ( k ) T + v as
q ( k ) T + v = | x ^ ( k ) T + v x ( k ) T + v | .
The cumulative value of { q ( k ) T + v } for a whole cycle is gathered, allowing for analyzing the machine condition periodically, which is calculated by
Q k = v = 1 T q ( k ) T + v ,
where Q k is the resulting value for the kth cycle.
Subsequently, we further standardize the obtained sequence { Q k } by s k = ( Q k Q ¯ ) / σ , where Q ¯ and σ are the sample mean and standard deviation of { Q k } .
The anomaly score { s k } describes the dynamic characteristics of machine condition monitoring over time. That is, when the condition of machine is stable, this implies that the condition change does not occur, and the values of { s k } will be relatively small; otherwise, the values will be large.

2.4. Simulation Validation

For the purpose of evaluation of the effectiveness of the proposed model, we applied it on synthetically generated testing data, where different levels of Gaussian White Noise are added for simulating the real engineering scenarios. Considering that the status change(s) in real machine operations can result in the change(s) of relevant CM variables in terms of amplitude, frequency, or both of them [7,18,41,42,43], the testing signals are formulated, respectively, as
x t = s i n ( ω × t ) + n o i s e , 1 t c 1 , 2 × s i n ( ω × t ) + n o i s e , c t m ,
x t = s i n ( ω × t ) + n o i s e , 1 t c 1 , s i n ( 2 ω × t ) + n o i s e , c t m ,
x t = s i n ( ω × t ) + n o i s e , 1 t c 1 , 2 × s i n ( 2 ω × t ) + n o i s e , c t m ,
where c is the change time, m is the length of data, and ω is set as π . The results are given in Figure 1, where the Figure 1a–c shows the simulated signals generated by Equations (11)–(13), respectively. For each subfigure, the SNR are set as 10, 30 and 50 dB, respectively. It can be clearly seen that the anomaly scores are stable before the status change, and then the occurrence of status change causes an abrupt increase of the anomaly score, which implies that it can be detected successfully. Here, we also found that, for amplitude change and amplitude and frequency change cases, the residuals of 30 dB and 50 dB are different from the one of the case 10 dB. The main possible reason is that, with a lower SNR, i.e., 10 dB, the noise will be more taken into account for prediction of the future value with an amplitude change; however, with a higher SNR, the noise will have less influence on the perdition process.

3. Hypothesis Testing for Decision-Making

On the basis of the resulting anomaly score s k , we test a null hypothesis using the 3 σ criterion in order to detect whether a change occurs at the current kth cycle by
H 0 : s k s ¯ k 1 < 3 σ , H 1 : s k s ¯ k 1 3 σ ,
where H 0 means that no change occurs on the kth cycle as long as s k s ¯ k 1 < 3 σ , and H 1 indicates that a change occurs when s k s ¯ k 1 3 σ (Note that the beginning nine cycles of data are used for initialization and the change detection begins with the 6th cycle after initialization). Here, s ¯ k 1 and σ are the mean and standard deviation of an assumed Gaussian distribution, respectively, and they are calculated by
s ¯ k 1 = 1 k 1 j = 1 k 1 s j ,
σ = 1 k 1 j = 1 k 1 ( s j s ¯ k 1 ) 2 .
Here, it is worth mentioning that there exist other alternatives such as Gaussian Mixed Model (GMM) [44] and other non-Gaussian assumptions [11,45] for change detection. However, since the focus of this paper is on the prediction model, we do not investigate these alternatives in this study.

4. Proposed Machine Running Status Monitoring Framework

Together with the DE model described in Section 2, an analysis of continuous monitoring of CM signals by means of a null hypothesis testing (Section 3) is proposed to inspect/diagnose whether an abnormal running status change occurs or not in successive machine operations. Total four steps are included in the framework as given below:
(1)
Collect CM data from the considered machine in a continuous manner;
(2)
Compute the prediction value using the proposed DE model;
(3)
Calculate anomaly scores based on residual error analysis at the current inspection time;
(4)
Make the change decision by testing a null hypothesis. Report an alarm to the user; otherwise, go to Step 2 to continue.
The flowchart of the proposed framework is provided in Figure 2.

5. Experimental Validation

To evaluate the effectiveness of the proposed method, we applied it to three representative industrial applications/tasks which are listed as below:
  • External loading status monitoring: External loading status is essential for condition monitoring during unsteady machine operations because a piece of equipment under operation may be exposed to a series of varying loads according to the user’s needs [46]. Moreover, the load is a critical operating condition factor which has significant impact on machine health [47]. Detection of changes in load condition makes it possible for the machine to adjust itself once a load change occurs for safety protection [48].
  • Bearing health status monitoring: As we all know, functional degeneration of machine components during the lifetime is common and unavoidable. The component degeneration will cause some undesired/unexpected consequences [48,49,50]. Based on CBM, maintenance can be scheduled in an optimal way with respect to cost, reliability, availability, or other logistic metrics of interest. Thus, automatic detection of changes in bearing health status can serve as a starting point for fault diagnosis or prediction of functional failure at an early stage.
  • Rotational speed monitoring: The rotational speed in machine operations may fluctuate due to condition variations or unsteady environments [51]. Speed condition monitoring helps to find the unexpected running behaviors for operation maintenance [52,53], thus highly desired in online process monitoring of industrial manufacturing, numerical controlled machining, ect.
The proposed method is compared with the ARIMA model proposed in [18]. In addition, we also compare it with two methods: RMS and kurtosis, which are widely used in the literature [36,37,38]. In the following experiments, we utilize the first detected alarm to determine whether the detection is successful or false. We then use the precision indicator which is defined as ratio of the number of correctly detected changes over the total number of changes to quantify the detection performance. However, for the case only including several testing data, we directly provide the detection result.

5.1. Case Study I: External Loading Status Monitoring

The testing data used in this section are the motor bearing data provided by CWRU [33]. In the experiment, the vibration data were collected from the drive end and the fan end of the motor driving mechanical system with the sampling frequency of 12 kHz [54].
The testing data for loading change detection was formed by concatenating each two load condition signals with fixed other relevant variables. We use the L L + Δ L to represent a loading status change, where L hp is the initializing load and Δ L hp is the changing load. More specifically, Δ L > 0 represents an increasing load change and Δ L < 0 represents a decreasing load change. The combination of L = { 0 , 1 , 2 , 3 } hp and Δ L = { 3 , 2 , 1 , 1 , 2 , 3 } hp (indicating different loading conditions in the employed CWRU data set) are used to simulate the loading changes. For example, when the load L begins at 3 and decreases by Δ L of 3 , the change can be expressed as 3 → 3 + ( 3 ) . In total, we have 12 combinations of simulated working load changes, as given in Table 1. In each combination, the vibration signals collected from drive end and fan end are used for analysis, respectively; thus, 24 signals are used for change detection.
Figure 3 shows an example of change detection for a working load change from 1 hp to 2 hp, where the original signal was collected from drive end with 1:10 down sampling and the change is labeled by a human instructor. From Figure 3a–d, we show the detection result by our method, ARIMA model, Kurtosis and RMS. It can be seen that our method can successfully detect the change cycle with a slight time delay. As seen in Figure 3a, the anomaly scores are relatively small and are distributed irregularly when the motor load is set around 1 hp. Then, the occurrence of loading status change causes an abrupt increase of the anomaly score which can be successfully detected through the 3 σ criterion based hypothesis testing; however, the method using the ARIMA model did not detect any alarm because the anomaly score is always in the confidence area, as shown in Figure 3b; the detected change cycle using the Kurtosis and RMS are much earlier than the actual change cycle, which belong to false detection as seen in Figure 3c,d; the main reason is that the proposed method has better power to detect weak changes than the other three methods.
Table 2 gives the comparison results. In the table, the N / 3 represents that the number of samples N in which change cycle is correctly detected in the three samples, N { 1 , 2 , 3 } . It can be clearly observed that our method achieves the best performance and the detection precision is 100%. In other words, whether the vibration signal is collected from the driven or the fan end, the proposed method can detect all the state changes precisely without generating any false alarms.

5.2. Case Study II: Bearing Health Status Monitoring

This section considers two real-engineering scenarios of bearing health status: a gradually degenerated process and a sharply degenerated process. Consequently, this case study is employed for this investigation.
The testing data set comes from the publicly available PRONOSTIA platform [34]. This system is designed to test and validate bearings fault detection, diagnostic, and prognostic approaches. The main purpose of PRONOSTIA platform is to provide real experimental data that characterize the degradation of ball bearings along their whole operational life (until their total failure). In order to reduce the bearing’s life duration, a radial force is set to the bearing’s maximum dynamic load of 4 kN and applied on the tested bearings. The force is generated by a cylinder pressure, and the pressure is delivered through a pressure regulator.
During the tests, accelerometers are fixed on the outer race of the bearing, and vibration signals are captured. The sampling rate is 25.6 kHz, the length of every sampled acceleration waveform is 2560 data points, i.e., 0.1 s, whereas recordings are repeated every 10 s. The vibration signals are transmitted into a PC for data visualization and storage through a National Instruments data acquisition (DAQ) card. Based on experience, it is considered as the end of life when the vibration level is greater than 20 g.
The experiment is carried out under three different operating conditions. In our case study, four representative bearing signals during the first operating condition are used for bearing early failure detection. The operating condition is under 1800 r/min (resolution per minute) with a 4000 N load.
Figure 4 shows detection results of the four bearings with different methods. From the original data of each Figure 4a,b, we can see that the vibration amplitude of bearing 1 and bearing 2 have gradually increasing trends, which indicates that the failure gets severe gradually. On the contrary, the vibration amplitude of bearing 3 and bearing 4 increases sharply at the end of lifetime as seen in the original data of each Figure 4c,d, which means a quick degradation processes. Table 3 gives the specific comparison results accordingly, where the number in the table represents the first early failure alarm point, and “N/A” represents false detection status change because there are too many alarm points and the first alarm point is far earlier than the early failure point. Thus, it can be clearly seen that kurtosis fails to distinguish the normal and the abnormal states, which can be verified by subfigure (3)s in all subfigures a, b, c and d in Figure 4. In addition, for bearing 3, as shown in subfigure (c), ARIMA and RMS have false alarms in a normal state. To summarize, the detection performance order is: our method > RMS > ARIMA > Kurtosis, revealing that our method has a superior ability compared to the others.

5.3. Case Study III: Speed Condition Monitoring

We used the experimental setup as shown in Figure 5 to collect testing data. In this setup, an accelerometer sensor mounted on the gearbox was used to acquire vibration signals with a sampling frequency of 1000 Hz, and then the collected signal data were transmitted to PC. During the data collection, we first set the initial speed of motor at ν , and then changed with an interval of Δ ν to simulate a speed status change (i.e., ν ν + Δ ν ). The testing values of ν and Δ ν were given as follows:
•  ν ϵ { 250 , 300 , 350 } rpm and Δ ν ϵ { 50 , 100 , 150 , 200 , 250 } rpm.
Therefore, we can obtain 15 parameter combinations with different speed condition changes as seen in Table 4. For each combination, we collected 10 data samples to form testing data set. Thus, we have a total number of 150 data sequences in our testing dataset.
Figure 6 shows an example of change detection for a speed change from 350 rpm to 400 rpm, where the change is labeled by a human instructor. From Figure 6a–d, we show the detection result by our method, ARIMA model, Kurtosis and RMS. It can be seen that our method can successfully detect the change cycle. As seen in Figure 6a, the anomaly sores show a remarkable increase at the change time such that the change was detected successfully. However, the method using the ARIMA model did not detect any alarm because the anomaly score is always in the confidence area, as shown in Figure 6b; the detected change cycle using the Kurtosis and RMS is much earlier than the actual change cycle, which is attributed to false detection, as seen in Figure 6c,d.
Table 5 gives the comparison results of the different methods under different settings of initial speed ν rpm. We can see that the proposed method also achieves a perfect detection performance on speed status change detection.

5.4. Results Summary and Discussion

To summarize, from the above experiments, we notice that the proposed differential equation-based prediction model has prior and remarkable detection performances in three different real applications, especially in weak state change. However, in the proposed method, we assume that the collected condition monitoring signal is periodically stationary. Although this assumption is reasonable and widely used in the literature [7,18,41], we have also noted that the actually collected signal in some cases may vary in periods, i.e., cycle non-stationary. With the purpose of extending the method to these cases, we would employ a preprocessing, e.g., angle re-sampling [55] and phase estimation using dynamic time warping [56,57], in order to suppress such temporal non-stationary before applying it to monitoring non-stationary condition signals.

6. Conclusions

In this paper, we have proposed a new method for monitoring and diagnosing the running condition of considered rotating machines. The base of the method is an early detection of condition changes based on the DE model that can realize one-step-ahead prediction of machine status. Together with this model, an analysis of continuously monitoring condition signal by means of sum residual error is presented to inspect/diagnose whether an abnormal condition change occurs or not based on a null hypothesis testing. The method was evaluated with three representative condition monitoring applications: external load status monitoring, bearing health status monitoring and speed condition monitoring. The method was also compared with those benchmark methods reported in the literature. The results demonstrated significant improvements of the proposed method over others, indicating its superiority and great potential in real engineering applications.
In the future, we will optimize the method to achieve a higher computational efficiency, and apply the method in the workshop for practical usages.

Author Contributions

X.W. conceived of the presented idea and developed the framework, performed the experiments and wrote the paper; G.C. contributed to the experiment; G.L. conceived the structure, provided guidance and modified the manuscript; Z.L. provided important guidance for this paper; and P.Y. gave a detailed revision. All authors have read and approved the final manuscript.

Funding

This work is supported in part by the National Natural Science Foundation of China (51775319), the National Key Research and Development Problem of China (2017YFF0105903), the Fundamental Research Funds of Shandong University (2018JC010), and the Young Scholars Program of Shandong University (YSPSDU, 2015WLJH30).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Do, P.; Voisin, A.; Levrat, E.; Iung, B. A proactive condition-based maintenance strategy with both perfect and imperfect maintenance actions. Reliab. Eng. Syst. Saf. 2015, 133, 22–32. [Google Scholar] [CrossRef]
  2. Wang, H.K.; Huang, H.Z.; Li, Y.F.; Yang, Y.J. Condition-based maintenance with scheduling threshold and maintenance threshold. IEEE Trans. Reliab. 2016, 65, 513–524. [Google Scholar] [CrossRef]
  3. Tahan, M.; Tsoutsanis, E.; Muhammad, M.; Karim, Z.A. Performance-based health monitoring, diagnostics and prognostics for condition-based maintenance of gas turbines: A review. Appl. Energy 2017, 198, 122–144. [Google Scholar] [CrossRef] [Green Version]
  4. Strączkiewicz, M.; Czop, P.; Barszcz, T. Supervised and unsupervised learning process in damage classification of rolling element bearings. Diagnostyka 2016, 17, 71–80. [Google Scholar]
  5. Wang, Z.; Bukkapatnam, S.T.; Kumara, S.R.; Kong, Z.; Katz, Z. Change detection in precision manufacturing processes under transient conditions. CIRP Ann.-Manuf. Technol. 2014, 63, 449–452. [Google Scholar] [CrossRef]
  6. Popescu, T.D. Blind separation of vibration signals and source change detection—Application to machine monitoring. Appl. Math. Model. 2010, 34, 3408–3421. [Google Scholar] [CrossRef]
  7. Reñones, A.; de Miguel, L.J.; Perán, J.R. Experimental analysis of change detection algorithms for multitooth machine tool fault detection. Mech. Syst. Signal Process. 2009, 23, 2320–2335. [Google Scholar] [CrossRef]
  8. Barszcz, T.; Randall, R.B. Application of spectral kurtosis for detection of a tooth crack in the planetary gear of a wind turbine. Mech. Syst. Signal Process. 2009, 23, 1352–1365. [Google Scholar] [CrossRef]
  9. Xiao, Y.; Wang, H.; Zhang, L.; Xu, W. Two methods of selecting Gaussian kernel parameters for one-class SVM and their application to fault detection. Knowl.-Based Syst. 2014, 59, 75–84. [Google Scholar] [CrossRef]
  10. Giantomassi, A.; Ferracuti, F.; Iarlori, S.; Ippoliti, G.; Longhi, S. Electric Motor Fault Detection and Diagnosis by Kernel Density Estimation and Kullback-Leibler Divergence Based on Stator Current Measurements. IEEE Trans. Ind. Electron. 2015, 62, 1770–1780. [Google Scholar] [CrossRef]
  11. Ge, Z.; Kruger, U.; Lamont, L.; Xie, L.; Song, Z. Fault detection in non-Gaussian vibration systems using dynamic statistical-based approaches. Mech. Syst. Signal Process. 2010, 24, 2972–2984. [Google Scholar] [CrossRef]
  12. Yan, H.C.; Zhou, J.H.; Pang, C.K. Gaussian mixture model using semisupervised learning for probabilistic fault diagnosis under new data categories. IEEE Trans. Instrum. Meas. 2017, 66, 723–733. [Google Scholar] [CrossRef]
  13. Rossi, T.M.; Braun, J.E. A statistical, rule-based fault detection and diagnostic method for vapor compression air conditioners. Hvac&R Res. 1997, 3, 19–37. [Google Scholar]
  14. Liang, X.; Wallace, S.A.; Nguyen, D. Rule-based data-driven analytics for wide-area fault detection using synchrophasor data. IEEE Trans. Ind. Appl. 2017, 53, 1789–1798. [Google Scholar] [CrossRef]
  15. Niu, G.; Yang, B.S. Dempster–Shafer regression for multi-step-ahead time-series prediction towards data-driven machinery prognosis. Mech. Syst. Signal Process. 2009, 23, 740–751. [Google Scholar] [CrossRef]
  16. Ballal, M.S.; Khan, Z.J.; Suryawanshi, H.M.; Sonolikar, R.L. Adaptive Neural Fuzzy Inference System for the Detection of Inter-Turn Insulation and Bearing Wear Faults in Induction Motor. IEEE Trans. Ind. Electron. 2007, 54, 250–258. [Google Scholar] [CrossRef]
  17. Soualhi, A.; Medjaher, K.; Zerhouni, N. Bearing health monitoring based on Hilbert-Huang transform, support vector machine, and regression. IEEE Trans. Instrum. Meas. 2015, 64, 52–62. [Google Scholar] [CrossRef]
  18. Lu, G.; Zhou, Y.; Lu, C.; Li, X. A novel framework of change-point detection for machine monitoring. Mech. Syst. Signal Process. 2017, 83, 533–548. [Google Scholar] [CrossRef]
  19. Yang, B.S.; Oh, M.S.; Tan, A.C.C. Machine condition prognosis based on regression trees and one-step-ahead prediction. Mech. Syst. Signal Process. 2008, 22, 1179–1193. [Google Scholar] [Green Version]
  20. Pham, H.T.; Yang, B.S. Estimation and forecasting of machine health condition using ARMA/GARCH model. Mech. Syst. Signal Process. 2010, 24, 546–558. [Google Scholar] [CrossRef]
  21. Randall, R.B.; Antoni, J.; Chobsaard, S. The relationship between spectral correlation and envelope analysis in the diagnostics of bearing faults and other cyclostationary machine signals. Mech. Syst. Signal Process. 2001, 15, 945–962. [Google Scholar] [CrossRef]
  22. Bouillaut, L.; Sidahmed, M. Cyclostationary approach and bilinear approach: Comparison, applications to early diagnosis for helicopter gearbox and classification method based on HOCS. Mech. Syst. Signal Process. 2001, 15, 923–943. [Google Scholar] [CrossRef]
  23. Boustany, R.; Antoni, J. A subspace method for the blind extraction of a cyclostationary source: Application to rolling element bearing diagnostics. Mech. Syst. Signal Process. 2005, 19, 1245–1259. [Google Scholar] [CrossRef]
  24. Borghesani, P.; Pennacchi, P.; Randall, R.B.; Sawalhi, N.; Ricci, R. Application of cepstrum pre-whitening for the diagnosis of bearing faults under variable speed conditions. Mech. Syst. Signal Process. 2013, 36, 370–384. [Google Scholar] [CrossRef]
  25. Antoni, J. Cyclic spectral analysis of rolling-element bearing signals: Facts and fictions. J. Sound Vib. 2007, 304, 497–529. [Google Scholar] [CrossRef]
  26. Wen, X.; Lu, G.; Yan, P. Improving Structural Change Detection using a Differential Equation-based Prediction Model for Condition Monitoring of Rotating Machines. In Proceedings of the IEEE International Conference on Prognostics and Health Management (ICPHM), Seattle, WA, USA, 11–13 June 2018; pp. 1–5. [Google Scholar]
  27. Chen, Y.; Yang, B.; Meng, Q.; Zhao, Y.; Abraham, A. Time-series forecasting using a system of ordinary differential equations. Inf. Sci. 2011, 181, 106–114. [Google Scholar] [CrossRef]
  28. Cao, H.; Kang, L.; Chen, Y.; Yu, J. Evolutionary modeling of systems of ordinary differential equations with genetic programming. Genet. Progr. Evolv. Mach. 2000, 1, 309–337. [Google Scholar] [CrossRef]
  29. Iversen, E.B.; Morales, J.M.; Møller, J.K.; Madsen, H. Short-term probabilistic forecasting of wind speed using stochastic differential equations. Int. J. Forecast. 2016, 32, 981–990. [Google Scholar] [CrossRef]
  30. Spencer, S.L.; Berryman, M.J.; Garcia, J.A.; Abbott, D. An ordinary differential equation model for the multistep transformation to cancer. J. Theor. Biol. 2004, 231, 515–524. [Google Scholar] [CrossRef] [Green Version]
  31. Zhao, Z.; Wang, J.; Zhao, J.; Su, Z. Using a grey model optimized by differential evolution algorithm to forecast the per capita annual net income of rural households in China. Omega 2012, 40, 525–532. [Google Scholar] [CrossRef]
  32. Sauhatas, A.; Bezrukovs, D. The application of stochastic differential equation models in the assessment of the economic feasibility of wind energy projects in Latvia. In Proceedings of the 57th International Scientific Conference on Power and Electrical Engineering of Riga Technical University (RTUCON), Riga, Latvia, 13–14 October 2016; pp. 1–6. [Google Scholar]
  33. Case Western Reserve University Bearing Data Center Website. Available online: http://csegroups.case.edu/bearingdatacenter/home (accessed on 20 January 2018).
  34. Nectoux, P.; Gouriveau, R.; Medjaher, K.; Ramasso, E.; Chebel-Morello, B.; Zerhouni, N.; Varnier, C. PRONOSTIA: An experimental platform for bearings accelerated degradation tests. In Proceedings of the IEEE International Conference on Prognostics and Health Management PHM’12, Denver, CO, USA, 18–21 June 2012; pp. 1–8. [Google Scholar]
  35. FEMTO Bearing Data Set. Available online: https://ti.arc.nasa.gov/tech/dash/groups/pcoe/prognostic-data-repository/ (accessed on 20 January 2018).
  36. Loutas, T.H.; Roulias, D.; Georgoulas, G. Remaining useful life estimation in rolling bearings utilizing data-driven probabilistic e-support vectors regression. IEEE Trans. Reliab. 2013, 62, 821–832. [Google Scholar] [CrossRef]
  37. Ali, J.B.; Chebel-Morello, B.; Saidi, L.; Malinowski, S.; Fnaiech, F. Accurate bearing remaining useful life prediction based on Weibull distribution and artificial neural network. Mech. Syst. Signal Process. 2015, 56, 150–172. [Google Scholar]
  38. Li, N.; Lei, Y.; Lin, J.; Ding, S.X. An improved exponential model for predicting remaining useful life of rolling element bearings. IEEE Trans. Ind. Electron. 2015, 62, 7762–7773. [Google Scholar] [CrossRef]
  39. Antoni, J.; Randall, R.B. Differential diagnosis of gear and bearing faults. J. Vib. Acoust. 2002, 124, 165–171. [Google Scholar] [CrossRef]
  40. Hu, J.W.; Tang, H.M. Numerical Methods for Differential Equations; City University: Hong Kong, China, 2003. [Google Scholar]
  41. Raad, A.; Antoni, J.; Sidahmed, M. Indicators of cyclostationarity: Theory and application to gear fault monitoring. Mech. Syst. Signal Process. 2008, 22, 574–587. [Google Scholar] [CrossRef]
  42. Figlus, T. The application of a continuous wavelet transform for diagnosing damage to the timing chain tensioner in a motorcycle engine. J. Vibroeng. 2015, 17, 1286–1294. [Google Scholar]
  43. Figlus, T.; Stańczyk, M. A method for detecting damage to rolling bearings in toothed gears of processing lines. Metalurgija 2016, 55, 75–78. [Google Scholar]
  44. Liao, L.; Lee, J. A novel method for machine performance degradation assessment based on fixed cycle features test. J. Sound Vib. 2009, 326, 894–908. [Google Scholar] [CrossRef]
  45. Tan, F.; Jiang, Z.; Bae, S.J. Generalized linear mixed models for reliability analysis of multi-copy repairable systems. IEEE Trans. Reliab. 2007, 56, 106–114. [Google Scholar] [CrossRef]
  46. Bartelmus, W.; Chaari, F.; Zimroz, R.; Haddar, M. Modelling of gearbox dynamics under time-varying nonstationary load for distributed fault detection and diagnosis. Eur. J. Mech.-A/Solids 2010, 29, 637–646. [Google Scholar] [CrossRef]
  47. Bartelmus, W.; Zimroz, R. Vibration condition monitoring of planetary gearbox under varying external load. Mech. Syst. Signal Process. 2009, 23, 246–257. [Google Scholar] [CrossRef]
  48. Zhao, F.; Tian, Z.; Bechhoefer, E.; Zeng, Y. An integrated prognostics method under time-varying operating conditions. IEEE Trans. Reliab. 2015, 64, 673–686. [Google Scholar] [CrossRef]
  49. Zhang, X.; Xu, R.; Kwan, C.; Liang, S.Y.; Xie, Q.; Haynes, L. An integrated approach to bearing fault diagnostics and prognostics. In Proceedings of the American Control Conference, Portland, OR, USA, 8–10 June 2005; pp. 2750–2755. [Google Scholar]
  50. Sawalhi, N.; Randall, R.B. Vibration response of spalled rolling element bearings: Observations, simulations and signal processing techniques to track the spall size. Mech. Syst. Signal Process. 2011, 25, 846–870. [Google Scholar] [CrossRef]
  51. McCormick, A.; Nandi, A. Cyclostationarity in rotating machine vibrations. Mech. Syst. Signal Process. 1998, 12, 225–242. [Google Scholar] [CrossRef]
  52. Combet, F.; Gelman, L. An automated methodology for performing time synchronous averaging of a gearbox signal without speed sensor. Mech. Syst. Signal Process. 2007, 21, 2590–2606. [Google Scholar] [CrossRef] [Green Version]
  53. Sawalhi, N.; Randall, R.B. Gear parameter identification in a wind turbine gearbox using vibration signals. Mech. Syst. Signal Process. 2014, 42, 368–376. [Google Scholar] [CrossRef]
  54. Smith, W.A.; Randall, R.B. Rolling element bearing diagnostics using the Case Western Reserve University data: A benchmark study. Mech. Syst. Signal Process. 2015, 64, 100–131. [Google Scholar] [CrossRef]
  55. Abboud, D.; Antoni, J.; Eltabach, M.; Sieg-Zieba, S. Angle\time cyclostationarity for the analysis of rolling element bearing vibrations. Measurement 2015, 75, 29–39. [Google Scholar] [CrossRef]
  56. Aqmar, M.R.; Fujihara, Y.; Makihara, Y.; Yagi, Y. Gait recognition by fluctuations. Comput. Vis. Image Underst. 2014, 126, 38–52. [Google Scholar] [CrossRef] [Green Version]
  57. Makihara, Y.; Trung, N.T.; Nagahara, H.; Sagawa, R.; Mukaigawa, Y.; Yagi, Y. Phase registration of a single quasi-periodic signal using self dynamic time warping. In Asian Conference on Computer Vision; Springer: Berlin/Heidelberg, Germany, 2010; pp. 667–678. [Google Scholar]
Figure 1. Results on testing data containing different changes: (a) amplitude change; (b) frequency change; (c) amplitude and frequency change.
Figure 1. Results on testing data containing different changes: (a) amplitude change; (b) frequency change; (c) amplitude and frequency change.
Sensors 19 00412 g001
Figure 2. Flowchart of the proposed framework.
Figure 2. Flowchart of the proposed framework.
Sensors 19 00412 g002
Figure 3. An example of change detection of load condition change from 1 hp to 2 hp. (a) DE model; (b) ARIMA model; (c) kurtosis; (d) RMS.
Figure 3. An example of change detection of load condition change from 1 hp to 2 hp. (a) DE model; (b) ARIMA model; (c) kurtosis; (d) RMS.
Sensors 19 00412 g003
Figure 4. Results of change detection of bearing early failure detection: (a) bearing 1; (b) bearing 2; (c) bearing 3; (d) bearing 4. In each piece of testing data: (1) the result by DE model; (2) the result by ARIMA model; (3) the result by Kurtosis, and (4) the result by RMS.
Figure 4. Results of change detection of bearing early failure detection: (a) bearing 1; (b) bearing 2; (c) bearing 3; (d) bearing 4. In each piece of testing data: (1) the result by DE model; (2) the result by ARIMA model; (3) the result by Kurtosis, and (4) the result by RMS.
Sensors 19 00412 g004aSensors 19 00412 g004b
Figure 5. Experimental setup.
Figure 5. Experimental setup.
Sensors 19 00412 g005
Figure 6. An example of change detection of speed condition change from 350 rpm to 400 rpm. (a) DE model; (b) ARIMA model; (c) Kurtosis; (d) RMS.
Figure 6. An example of change detection of speed condition change from 350 rpm to 400 rpm. (a) DE model; (b) ARIMA model; (c) Kurtosis; (d) RMS.
Sensors 19 00412 g006
Table 1. Simulated working load change of L L + Δ L hp.
Table 1. Simulated working load change of L L + Δ L hp.
L \ L + Δ L 0123
0- 0 0 + ( 1 ) 0 0 + ( 2 ) 0 0 + ( 3 )
1 1 1 + ( 1 ) - 1 1 + ( 1 ) 1 1 + ( 2 )
2 2 2 + ( 2 ) 2 2 + ( 1 ) - 2 2 + ( 1 )
3 3 3 + ( 3 ) 3 3 + ( 2 ) 3 3 + ( 1 ) -
Table 2. Results of loading change for different setting initial working loads L.
Table 2. Results of loading change for different setting initial working loads L.
Case\L0123Average
Our methodDrive3/33/33/33/3100%
Fan3/33/33/33/3
ARIMADrive3/32/33/32/387.5%
Fan3/32/33/33/3
KurtosisDrive3/31/30/30/320.8%
Fan0/30/32/30/3
RMSDrive1/31/31/31/333.3%
Fan0/30/31/33/3
Table 3. Results of bearing early failure detection.
Table 3. Results of bearing early failure detection.
Method\CaseGradual DegenerationSharp Degeneration
Bearing 1Bearing 2Bearing 3Bearing 4
Our method1441118324332190
ARIMA143512824242184
KurtosisN/AN/AN/AN/A
RMS142210834242208
Table 4. Simulated speed condition change ν ν + Δ ν rpm.
Table 4. Simulated speed condition change ν ν + Δ ν rpm.
Δ ν \ ν 250300350
50 250 250 + 50 300 300 + 50 350 350 + 50
100 250 250 + 100 300 300 + 100 350 350 + 100
150 250 250 + 150 300 300 + 150 350 350 + 150
200 250 250 + 200 300 300 + 200 350 350 + 200
250 250 250 + 250 300 300 + 250 350 350 + 250
Table 5. Results of speed condition change detection by different setting initial speeds ν rpm.
Table 5. Results of speed condition change detection by different setting initial speeds ν rpm.
250300350Average
Our method100%100%100%100%
ARIMA90%70%90%83.3%
Kurtosis72%86%82%80%
RMS100%44%80%74.6%

Share and Cite

MDPI and ACS Style

Wen, X.; Chen, G.; Lu, G.; Liu, Z.; Yan, P. Differential Equation-Based Prediction Model for Early Change Detection in Transient Running Status. Sensors 2019, 19, 412. https://doi.org/10.3390/s19020412

AMA Style

Wen X, Chen G, Lu G, Liu Z, Yan P. Differential Equation-Based Prediction Model for Early Change Detection in Transient Running Status. Sensors. 2019; 19(2):412. https://doi.org/10.3390/s19020412

Chicago/Turabian Style

Wen, Xin, Guangyuan Chen, Guoliang Lu, Zhiliang Liu, and Peng Yan. 2019. "Differential Equation-Based Prediction Model for Early Change Detection in Transient Running Status" Sensors 19, no. 2: 412. https://doi.org/10.3390/s19020412

APA Style

Wen, X., Chen, G., Lu, G., Liu, Z., & Yan, P. (2019). Differential Equation-Based Prediction Model for Early Change Detection in Transient Running Status. Sensors, 19(2), 412. https://doi.org/10.3390/s19020412

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop