# Prediction Model of Vibration Feature for Equipment Maintenance Based on Full Vector Spectrum

Authors: Lei Chen; Jie Han; Wenping Lei; ZhenHong Guan; Yajuan Gao
Journal: Shock and Vibration
 (2017)

---

## Abstract

Establishing a prediction model is a key step for the implementation of prognostic and health management. The prediction model can be used to forecast the change trend of the characteristics of the vibration signal and analyze the potential failure in the future. Taking the vibration of power plant steam turbine as an example, the full vector fusion and fault prediction were studied. Due to the fact that the evaluation of the machine fault with only one transducer may result in a fault judgement with partiality, an information fusion method based on the theory of full vector spectrum was adopted to extract the vibration feature. An autoregressive prediction model was established. The collected vibration signals with pairing channels were fused. The time sequence of the fused vectors and spectrums were used to build the prediction model. The amplitude of main vector of rotating frequency and spectrum order structure were analyzed and predicted. The uncertainty of the spectrum structure can be eliminated by the information fusion. The reliability of the fault prediction was improved. The study on vibration prediction model system laid a technical foundation for the fault prognostic research.


---

## Body

## 1. Introduction

By the techniques of prediction, equipment status deterioration laws can be drawn out from the machine historical state parameters and thus potential faults can be predicted. The maintenance plan can be arranged reasonably to ensure the safe operation of the equipment [1]. Prediction modeling and vibration data feature extraction are the key predictive maintenance technologies. In recent years, scholars have conducted lots of work on the aspect of equipment fault prediction, and many prediction methods are proposed [2–6]. These methods can be roughly divided into three categories: knowledge based, model based, and data based. Some expert systems and fuzzy logic [7–9] are the method based on knowledge. This method can make use of expert knowledge and experience, and it does not need to be very accurate model. But the acquisition of the knowledge is quite difficult. The method of filter prediction and the method of modeling based on the fault mechanism are the methods based on model [10–12]. When the mathematical model of the diagnosis object is accurate, this kind of method can get the accurate result of fault prediction. However, it is difficult to establish accurate mathematical model for the complex dynamic system. So this method is very limited and expensive. The prediction methods of autoregressive forecast, grey theory, neural network, and so forth are the methods based on data [4–6, 13–16]. By mining the implicit information from the industrial field data, this type of method has a wide range of engineering application. The prediction method based on data has a wide range and low cost. It is the most practical method and has become the research focus and development trend in the field of fault diagnosis and prediction. In this paper, an autoregressive forecasting model is established based on historical vibration data.When a fault prediction model was established, fault signature was obtained from vibration transducers. Although multiple transducers were used in the monitoring system, fault signature was still extracted from one of them. Based on the change of eigenvalues or characteristic vector, possible failures could be predicted. According to the theory of full spectrum and holographic and full vector spectrum [1, 17–19], the vibration track under a vibration harmonic is an ellipse. The vibration obtained by a transducer is the projection of the vibration trace in the transducer’s installation direction. The measurement results have certain limitation and one-sidedness, and even in some cases different conclusions will be drawn according to the transducers installed in different directions. So in this study signal vector was extracted based on the theory of full vector spectrum in the process of modeling to overcome the partial impact. The equipment vibration vector prediction model was established to predict degradation trend and analyze the possibility of equipment failure.


## 2. Structure of Full Vector Prediction

The structure of the full vector fault prediction and diagnosis system is shown in Figure1. The Arabic numerals 1~14 is the measuring points. In the same interface, two transducers are installed in X-Y or vertical layout. For example, in the A-A section, two transducers are installed in X and Y direction to acquire the rotor vibration. The data acquisition system transforms the analog vibration into digital data after signal conditioning, filter, and analog to digital conversion. The equipment running status can be evaluated and possible failure can be analyzed according to the result of prediction. When the rotor is running, the orbit of shaft center on each harmonic is an ellipse. Under each harmonic, the signal obtained by the transducer is the projection of the vibration trace in the corresponding direction. The overall signal measured by a transducer is the superposition of projections of all measurable frequencies. So signals of a single transducer cannot reflect the vibration of the equipment comprehensively and information fusion is necessary. In this study, the difference of spectrums is addressed between signals obtained by a signal transducer and signals after information fusion. The time sequence vibration vector is used as input parameters and a data driven AR prediction model is built with fixed order. The prediction model’s coefficient was updated with the changing of the vibration signal to improve the model with the ability of self-correcting.Figure 1
Diagram of prediction system based on full vector spectrum.


## 3. Feature Extraction and Prediction

Since, in the same measurement section, vibration obtained by a single transducer has certain one-sidedness, information fusion is needed to extract the vibration vector. The vibration vector is the fused signal of the two transducers installed in different directions. It can reflect the actual circumstances of the machine vibration. It is a feasible way to judge the trend of equipment operating status by vibration vector.


### 3.1. Extract of Vibration Vector

Assuming that the two transducers ofx and y direction were installed in section A-A of Figure 1 to measure the rotor vibration, vibrations acquired by the two transducers of a steady rotating machinery are a combination of motion orbits for some harmonic waves wi(i=1,2,…). For a certain harmonic i, assume the following definitions: Xi stands for the amplitude of vibration signal in the direction x; Yi stands for the amplitude of vibration signal of the harmonic i in the direction y; ∅xi stands for the phase angle of the harmonic i in the direction x; ∅yi stands for the phase angle of the harmonic i in the direction y; Xci stands for the amplitude of vibration in direction x of harmonic i when t=0; Yci stands for the amplitude of vibration in direction y of harmonic i when t=0; Xsi and Ysi have no physical meaning and were used for the geometry calculation of harmonic track. Then there are the following equations:(1)Xci=Xicos⁡∅xiXsi=Xisin⁡∅xiYci=Yicos⁡∅yiYsi=Yisin⁡∅yi.The motion equations inX and Y direction can be generally expressed as(2)x=∑i=1∞Xicos⁡ωit+∅xi=Re∑i=1∞X-iejωit=Re∑i=1∞xci+jxsiejωity=∑i=1∞Yicos⁡ωit+∅yi=Re∑i=1∞Y-iejωit=Re∑i=1∞yci+jysiejωit.After some trigonometric function calculations, (3) can be reached. The equation indicates that the vibration trajectory is an ellipse in every harmonic i(i=1,2,3,…).(3)Yci2+Ysi2xi2+Xci2+Xsi2yi2-2XciYsi+XsiYcixiyi=XsiYsi-XciYci2.As shown in Figure2, RL is the half major axis of the ellipse. RS is the half minor axis. Elliptic trajectory in kinematics can be regarded as the synthesis of two circular motions of the same frequency and opposite direction [1, 17, 18]. There exists the following equation.(4)z-=x-p+x-r.Figure 2
Ellipse orbit.In Figure2, z- represents a point on the rotor trajectory at a certain time. x-p is the forward precession circle vector with phase angle ∅p. x-r is the backward precession circle vector with phase angle ∅r.If the ellipse was drawn in the same coordinate system, the running state of the rotor can be understood directly [18]. Industry wide of a single transducer to assess dynamics has been replaced by use of two or more transducers and key-phasor as fully described in detail in the classic works by Dr. Donald Bently [17]. Dr. Bently’s classic work methodically and accurately explains in great detail the use of diagnostic equipment, theory, and applications. Theory of full spectrum proposed by Dr. Bently considers that the radius of the forward circle and backward circle which form the harmonic ellipse can be displayed, respectively, according to the frequency axis of positive and negative. By comparing the amplitudes of the harmonic under positive and negative frequencies, the direction of precession of the ellipse and vibration intensity can be determined.The maximum displacement of the rotor is at the vertex point of the major axis of the ellipse. The major axis of the ellipse is the biggest displacement of the vibration. It can stand for the size of the rotor vibration. Here the ellipse’s long half axis was called the main vibration vector, and the short half axis was called the assistant vibration vector.According to the theory of full vector spectrum, the main vector of vibration harmonic can be simplified as a simple calculation by FFT [1, 2]. Assuming that xi and yi are the signals acquired by the transducers x and y. Here i=0,1,2,…,N-1, N are the sampling points. A complex sequence zi can be constituted by the xi and yi; that is, zi=xi+jyi. Zk can be obtained by Fourier transform of zi data series; here k=1,2,…,N/2-1. After a series of calculations, (5) can be obtained.(5)RLk=12NZk+ZN-kRSk=12NZk-ZN-kk=0,1,2,…,N2-1.In (5), RLk and RSk are the main and assistant vibration vectors. Its spectral structure can reflect the operation state of the machine comprehensively.


### 3.2. Autoregressive Prediction Model

Prediction method of vibration vector was based on autoregressive (AR) model. Autoregressive model is a common form of time series [20]. It is a special case of Autoregressive Moving Average (ARMA) model.For a stationary time sequence{xt}, if the value of xt is concerned not only with the data xt-1,xt-2,…,xt-n(n=1,2,…,n) but also with the disturbance αt-1,αt-2,αt-m(m=1,2,…,m), the general structure of ARMA (n, m) according to the idea of multiple linear regression is as follows:(6)xt=∑i=1nφixt-i-∑j=1mθjαt-j+αt.Heren stands for the order of AR(n); m stands for the order of MA(m); φi(i=1,2,…,n) stands for the parameters of autoregressive; θj(j=1,2,…,m) stands for the parameters of moving average; αt stands for a sequence of independent and identically distributed random variables. The mean of αt:Eαt=0; the variance of αt:Varαt=σε2>0.Whenθj=0, the ARMA(n, m) model has no moving average. At this circumstance, the model was called autoregressive model with order n, that is, AR(n). Its structure is as follows:(7)xt=∑i=1nφixt-i+αt.The parameter estimation of AR model is linear estimation, which is fast and robust, and it is easy to be used in practical engineering. ARMA(n, m) model is equivalent to the AR(n) model in the case of large order, so the fault prediction can be carried out using AR(n) model.The establishment of the model is shown in Figure3. The first step is to identify the model, that is, to determine whether the vibration vector data sequence can be used for the AR prediction. The basic condition of the AR model is that the time sequence data must meet the requirement of stability. Whether the data is stable can be determined by the changes of self-correlation coefficient ρj with j [21, 22]. If ρj approach rapidly zero with the increase of j, the time series data meets the requirement of stability. On the contrary, it does not meet the requirement of stability. For the nonstationary data, advanced stationary processing is needed.Figure 3
Process of prediction model establishment.The order of the model is an important factor that affects the accuracy of the prediction results. Once the order number is selected, the prediction model structure is determined. At present, the commonly used order selection criteria were information criterion [23–25] and the empirical criteria given by the field experts [26]. In the literature [27], some empirical criterion for selecting the model order under different data number of time series was proposed. The selection of the order also needs to combine the experimental analysis with the current data.For the AR(n) model, the parameters needing estimation were the mean μ, variance σα2, and the coefficients φ=[φ1,φ2,…,φp].Mean value was estimated by(8)μ=X-=∑i=1nXin.Self-covariance was estimated by(9)γk^=1n∑i=1n-kZi-Z-Zi+k-Z-.For the sample{Xt~}, if the self-covariance γk is known, the solution can be obtained by the Yule-Waller estimation equation. That is, (10)γ1γ2⋮γn=γ0γ1⋯γn-1γ1γ2⋯γn-2⋮⋮⋮γn-1γn-2⋯γ0φ1φ2⋮φn.Replaceγk of the formula (10) with γk^; the parameter estimates φ^=φ^1,φ^2,…,φ^n will be obtained. After determining the order of the AR model and the coefficients of each order, the prediction model is established. Statistical value test can be used to determine whether the model can be used for prediction.


### 3.3. Application of Prediction Model

The AR prediction model is established according to the historical data. In the process of using, as the data is updated, the number of order of the AR model remains unchanged, but the coefficients of each order are revised. A relational database based on Microsoft SQL Server was built in the monitoring system for the purpose of facilitating the trend analysis and fault prediction. There are real-time data tables and historical data tables in the database. Every two seconds the acquisition system collects a set of data for real-time monitoring. The prediction system reads the data from the historical data table and sends it into the established AR model for the trend analysis and numerical prediction. Historical data is divided into seven levels. The first level stores the data of the last 8 hours; the second level stores the data of the last day; the third level stores the data of the last week; the fourth level stores the data of the last month; the fifth level stores the data of the last three months; the sixth level stores the data for the last six months; and the seventh level stores the data for the last year. Historical data older than one year is no longer used for fault prediction. When the number of predicted steps was determined, the time span of prediction results was different with the different time data series.


## 3.1. Extract of Vibration Vector

Assuming that the two transducers ofx and y direction were installed in section A-A of Figure 1 to measure the rotor vibration, vibrations acquired by the two transducers of a steady rotating machinery are a combination of motion orbits for some harmonic waves wi(i=1,2,…). For a certain harmonic i, assume the following definitions: Xi stands for the amplitude of vibration signal in the direction x; Yi stands for the amplitude of vibration signal of the harmonic i in the direction y; ∅xi stands for the phase angle of the harmonic i in the direction x; ∅yi stands for the phase angle of the harmonic i in the direction y; Xci stands for the amplitude of vibration in direction x of harmonic i when t=0; Yci stands for the amplitude of vibration in direction y of harmonic i when t=0; Xsi and Ysi have no physical meaning and were used for the geometry calculation of harmonic track. Then there are the following equations:(1)Xci=Xicos⁡∅xiXsi=Xisin⁡∅xiYci=Yicos⁡∅yiYsi=Yisin⁡∅yi.The motion equations inX and Y direction can be generally expressed as(2)x=∑i=1∞Xicos⁡ωit+∅xi=Re∑i=1∞X-iejωit=Re∑i=1∞xci+jxsiejωity=∑i=1∞Yicos⁡ωit+∅yi=Re∑i=1∞Y-iejωit=Re∑i=1∞yci+jysiejωit.After some trigonometric function calculations, (3) can be reached. The equation indicates that the vibration trajectory is an ellipse in every harmonic i(i=1,2,3,…).(3)Yci2+Ysi2xi2+Xci2+Xsi2yi2-2XciYsi+XsiYcixiyi=XsiYsi-XciYci2.As shown in Figure2, RL is the half major axis of the ellipse. RS is the half minor axis. Elliptic trajectory in kinematics can be regarded as the synthesis of two circular motions of the same frequency and opposite direction [1, 17, 18]. There exists the following equation.(4)z-=x-p+x-r.Figure 2
Ellipse orbit.In Figure2, z- represents a point on the rotor trajectory at a certain time. x-p is the forward precession circle vector with phase angle ∅p. x-r is the backward precession circle vector with phase angle ∅r.If the ellipse was drawn in the same coordinate system, the running state of the rotor can be understood directly [18]. Industry wide of a single transducer to assess dynamics has been replaced by use of two or more transducers and key-phasor as fully described in detail in the classic works by Dr. Donald Bently [17]. Dr. Bently’s classic work methodically and accurately explains in great detail the use of diagnostic equipment, theory, and applications. Theory of full spectrum proposed by Dr. Bently considers that the radius of the forward circle and backward circle which form the harmonic ellipse can be displayed, respectively, according to the frequency axis of positive and negative. By comparing the amplitudes of the harmonic under positive and negative frequencies, the direction of precession of the ellipse and vibration intensity can be determined.The maximum displacement of the rotor is at the vertex point of the major axis of the ellipse. The major axis of the ellipse is the biggest displacement of the vibration. It can stand for the size of the rotor vibration. Here the ellipse’s long half axis was called the main vibration vector, and the short half axis was called the assistant vibration vector.According to the theory of full vector spectrum, the main vector of vibration harmonic can be simplified as a simple calculation by FFT [1, 2]. Assuming that xi and yi are the signals acquired by the transducers x and y. Here i=0,1,2,…,N-1, N are the sampling points. A complex sequence zi can be constituted by the xi and yi; that is, zi=xi+jyi. Zk can be obtained by Fourier transform of zi data series; here k=1,2,…,N/2-1. After a series of calculations, (5) can be obtained.(5)RLk=12NZk+ZN-kRSk=12NZk-ZN-kk=0,1,2,…,N2-1.In (5), RLk and RSk are the main and assistant vibration vectors. Its spectral structure can reflect the operation state of the machine comprehensively.


## 3.2. Autoregressive Prediction Model

Prediction method of vibration vector was based on autoregressive (AR) model. Autoregressive model is a common form of time series [20]. It is a special case of Autoregressive Moving Average (ARMA) model.For a stationary time sequence{xt}, if the value of xt is concerned not only with the data xt-1,xt-2,…,xt-n(n=1,2,…,n) but also with the disturbance αt-1,αt-2,αt-m(m=1,2,…,m), the general structure of ARMA (n, m) according to the idea of multiple linear regression is as follows:(6)xt=∑i=1nφixt-i-∑j=1mθjαt-j+αt.Heren stands for the order of AR(n); m stands for the order of MA(m); φi(i=1,2,…,n) stands for the parameters of autoregressive; θj(j=1,2,…,m) stands for the parameters of moving average; αt stands for a sequence of independent and identically distributed random variables. The mean of αt:Eαt=0; the variance of αt:Varαt=σε2>0.Whenθj=0, the ARMA(n, m) model has no moving average. At this circumstance, the model was called autoregressive model with order n, that is, AR(n). Its structure is as follows:(7)xt=∑i=1nφixt-i+αt.The parameter estimation of AR model is linear estimation, which is fast and robust, and it is easy to be used in practical engineering. ARMA(n, m) model is equivalent to the AR(n) model in the case of large order, so the fault prediction can be carried out using AR(n) model.The establishment of the model is shown in Figure3. The first step is to identify the model, that is, to determine whether the vibration vector data sequence can be used for the AR prediction. The basic condition of the AR model is that the time sequence data must meet the requirement of stability. Whether the data is stable can be determined by the changes of self-correlation coefficient ρj with j [21, 22]. If ρj approach rapidly zero with the increase of j, the time series data meets the requirement of stability. On the contrary, it does not meet the requirement of stability. For the nonstationary data, advanced stationary processing is needed.Figure 3
Process of prediction model establishment.The order of the model is an important factor that affects the accuracy of the prediction results. Once the order number is selected, the prediction model structure is determined. At present, the commonly used order selection criteria were information criterion [23–25] and the empirical criteria given by the field experts [26]. In the literature [27], some empirical criterion for selecting the model order under different data number of time series was proposed. The selection of the order also needs to combine the experimental analysis with the current data.For the AR(n) model, the parameters needing estimation were the mean μ, variance σα2, and the coefficients φ=[φ1,φ2,…,φp].Mean value was estimated by(8)μ=X-=∑i=1nXin.Self-covariance was estimated by(9)γk^=1n∑i=1n-kZi-Z-Zi+k-Z-.For the sample{Xt~}, if the self-covariance γk is known, the solution can be obtained by the Yule-Waller estimation equation. That is, (10)γ1γ2⋮γn=γ0γ1⋯γn-1γ1γ2⋯γn-2⋮⋮⋮γn-1γn-2⋯γ0φ1φ2⋮φn.Replaceγk of the formula (10) with γk^; the parameter estimates φ^=φ^1,φ^2,…,φ^n will be obtained. After determining the order of the AR model and the coefficients of each order, the prediction model is established. Statistical value test can be used to determine whether the model can be used for prediction.


## 3.3. Application of Prediction Model

The AR prediction model is established according to the historical data. In the process of using, as the data is updated, the number of order of the AR model remains unchanged, but the coefficients of each order are revised. A relational database based on Microsoft SQL Server was built in the monitoring system for the purpose of facilitating the trend analysis and fault prediction. There are real-time data tables and historical data tables in the database. Every two seconds the acquisition system collects a set of data for real-time monitoring. The prediction system reads the data from the historical data table and sends it into the established AR model for the trend analysis and numerical prediction. Historical data is divided into seven levels. The first level stores the data of the last 8 hours; the second level stores the data of the last day; the third level stores the data of the last week; the fourth level stores the data of the last month; the fifth level stores the data of the last three months; the sixth level stores the data for the last six months; and the seventh level stores the data for the last year. Historical data older than one year is no longer used for fault prediction. When the number of predicted steps was determined, the time span of prediction results was different with the different time data series.


## 4. Test and Experiment

### 4.1. Comparison of Full Vector Spectrum and Traditional Spectrum

In order to verify the correctness of the method for data fusion and prediction, an online monitoring system was installed in a power plant in China, and the vibration data were collected. The steam turbine of the power plant drives the feed water pump to supply water for the boiler. The rotor structure of the steam turbine feed pump group is composed of a steam turbine rotor, a feed water pump shaft, and a flexible coupling which is used for connecting the rotor and the feed pump shaft. The steam turbine feed pump on-site photo is shown in Figure4. Figure 5 is the internal structure. Vibration measuring points are shown in Figure 6. In the figure, S1~S10 are eddy current transducers. S1~S8 were used to measure the rotor radial displacement; S9 was used to measure the axial position of the rotor; S10 was used to measure the rotating speed and provide the key phase signal. S5 and S6 are eddy current transducers installed perpendicular to each other, which are used to monitor the radial vibration of the rotor at the output end of the turbine. S5 was installed in the x direction. S6 was installed in the y direction.Figure 4
Photo of the steam turbine feed pump.Figure 5
Internal structure of the steam turbine feed pump.Figure 6
Layout of vibration measuring point.The sampling frequency is 1600 Sa/s; sampling length is 1024 points; the rotor speed is 3000 r/min. The vibration waveforms of S5 and S6 are shown in Figure7. In time domain, the waveforms indicate that the overall vibration amplitude of X and Y direction is of considerable size. In frequency domain as shown in Figure 8, the frequency structure of x and y direction is quite different. According to the spectrum of x direction, it seems that the vibration amplitude mainly relates to the double frequency, and perhaps there exists fault of misalignment; however, after observing the spectrum in y direction, it seems that the vibration relates to the rotating frequency. Therefore, the fault may be related to the unbalance. Different diagnosis conclusions may be drawn on the machine by signals of different directions on the same rotor interface. Making a diagnosis conclusion with signals of any one direction will be one-sided to some extent.Figure 7
Vibration wave ofx and y direction of the same cross section.
(a)
Wave ofx-direction
(b)
Wave ofy-directionFigure 8
Spectrum ofx and y direction of the same cross section.
(a)
Spectrum ofx-direction
(b)
Spectrum ofy-directionFigure9 is a fusion signal spectrum of x and y direction. Figure 9(a) is full vector spectrum; Figure 9(b) is the shaft orbit of the signal at different order. The synthetic signal of X and Y transducer at each harmonic is a motion ellipse, as shown in Figure 9(b). The major axis of the ellipse reflects the most significant vibration. No matter where the transducer mounts, the main vibration vector is the maximum amplitude of the vibration in the motion plane. On the same rotor section, the main vibration vector is not affected by the radial mounting position of the two transducers installed vertically with each other. The full vector spectrum is the synthesis of two homologous transducers’ information. It works as a virtual transducer in the plane to detect the maximum vibration of the interface at each frequency. So its spectral structure can reflect the operation state of the machine comprehensively. The vibration amplitude distribution of the full vector spectrum is in accordance with shaft orbit; Figure 9 shows that the vibration amplitude mainly relates to the double frequency, and the actual situation is also the case.Figure 9
Full vector spectrum and shaft orbit.
(a)
Full vector spectrum
(b)
Shaft orbit


### 4.2. Prediction of Full Vector Amplitude

The fault severity can be determined by analyzing the vibration amplitude. Full vector data fusion was carried out based on the collected vibration signals to extract the full vector amplitude of rotating frequency. The AR(4) data prediction model was established. The time interval is 2 hours. The method is suitable for the long-term monitoring and numerical prediction of stable operating machines. Figure 10 indicates that the predicted value is almost the same as the measured vibration value. The established AR prediction model of full vector information fusion could predict the change of vibration effectively. By using full vector amplitude, the one-sidedness of one single transducer can be overcome.Figure 10
Comparison between predicted and measured values.


### 4.3. Prediction of Full Vector Spectrum

The fault position and property can be determined in accordance with the machine structure by analysis the spectrum structure [28]. To predict the fault property, spectrum structure prediction is a feasible method.Obtaining signal with one single transducer, whether the transducer mounts inx direction or in y direction, the acquired signal is one-sided for fault prediction. The reliability of conclusion concerning fault diagnosis and fault character prediction with such signals cannot be guaranteed. Date fusion must be carried out to enhance the reliability of spectrum. The data fusion is based on the theory of full vector spectrum, and the main vector was used for the prediction of spectrum structure. The main vector sequence reflects the distribution of the vibration energy along the frequency. Compared with the spectrum in one direction only, the full vector spectral structure can reflect the operation state of the machine comprehensively. The result of spectrum structure prediction with AR(13) is shown in Figure 11. The potential fault property can be determined by analyzing the spectrum structure in the prediction domain and reasonable maintenance can be arranged.Figure 11
Full vector spectrum prediction with AR(13).


### 4.4. Equipment Maintenance Using Predicted Feature

Vibration analysis is one of the important methods for fault diagnosis of rotating machinery. By collecting the vibration signal and extracting the eigenvector, the established vibration analysis and prediction system can be used as an important support for equipment health management. By analyzing and predicting the vibration feature, it is possible to evaluate the operating status of the equipment, predict deterioration trends, and identify potential failures that may occur.The established AR prediction model can predict the changes of vibration amplitude and full vector spectrum. By predicting the change of the amplitude, it is possible to determine the deterioration tendency of the equipment and to predict how long the alarm will occur. When the alarm is predicted to occur in the near future, the full vector spectrum prediction model is called to check the amplitude changes at different characteristic frequencies. The characteristic frequency is related to the corresponding equipment failure; for example, rotating frequency is related to the failure of dynamic unbalance. By examining the changes of the amplitude of the characteristic frequency in the spectrum, it is possible to judge out the possible fault to occur. The vibration feature prediction can provide a decision-making basis for the equipment management and maintenance department.


### 4.5. Discussions

In the machine condition monitoring system, the AR model is an effective prediction method. It can be used to predict the change value of the vibration. The structure of the prediction model can be determined by the order number of the AR model. The order of the model was fixed. The coefficients of each order can be determined by the calculation of the sample data. In the actual condition monitoring system, the sampling data was changed as time goes on. The coefficients of each order can be modified to fit in with the latest running state by calculating the sample data in the latest period of time.For the data feature extraction, since the full vector amplitude can reflect the maximum vibration amplitude of the equipment, the prediction model was constructed based on the full vector amplitude. The AR amplitude prediction model can forecast the change trend of the equipment effectively. By analyzing the structure characteristics and the change of full vector spectrum, the fault property and failure position can be predicted. Compared with signals with a single transducer, the full vector spectrum is more effective to reflect the fault condition of the equipment.


## 4.1. Comparison of Full Vector Spectrum and Traditional Spectrum

In order to verify the correctness of the method for data fusion and prediction, an online monitoring system was installed in a power plant in China, and the vibration data were collected. The steam turbine of the power plant drives the feed water pump to supply water for the boiler. The rotor structure of the steam turbine feed pump group is composed of a steam turbine rotor, a feed water pump shaft, and a flexible coupling which is used for connecting the rotor and the feed pump shaft. The steam turbine feed pump on-site photo is shown in Figure4. Figure 5 is the internal structure. Vibration measuring points are shown in Figure 6. In the figure, S1~S10 are eddy current transducers. S1~S8 were used to measure the rotor radial displacement; S9 was used to measure the axial position of the rotor; S10 was used to measure the rotating speed and provide the key phase signal. S5 and S6 are eddy current transducers installed perpendicular to each other, which are used to monitor the radial vibration of the rotor at the output end of the turbine. S5 was installed in the x direction. S6 was installed in the y direction.Figure 4
Photo of the steam turbine feed pump.Figure 5
Internal structure of the steam turbine feed pump.Figure 6
Layout of vibration measuring point.The sampling frequency is 1600 Sa/s; sampling length is 1024 points; the rotor speed is 3000 r/min. The vibration waveforms of S5 and S6 are shown in Figure7. In time domain, the waveforms indicate that the overall vibration amplitude of X and Y direction is of considerable size. In frequency domain as shown in Figure 8, the frequency structure of x and y direction is quite different. According to the spectrum of x direction, it seems that the vibration amplitude mainly relates to the double frequency, and perhaps there exists fault of misalignment; however, after observing the spectrum in y direction, it seems that the vibration relates to the rotating frequency. Therefore, the fault may be related to the unbalance. Different diagnosis conclusions may be drawn on the machine by signals of different directions on the same rotor interface. Making a diagnosis conclusion with signals of any one direction will be one-sided to some extent.Figure 7
Vibration wave ofx and y direction of the same cross section.
(a)
Wave ofx-direction
(b)
Wave ofy-directionFigure 8
Spectrum ofx and y direction of the same cross section.
(a)
Spectrum ofx-direction
(b)
Spectrum ofy-directionFigure9 is a fusion signal spectrum of x and y direction. Figure 9(a) is full vector spectrum; Figure 9(b) is the shaft orbit of the signal at different order. The synthetic signal of X and Y transducer at each harmonic is a motion ellipse, as shown in Figure 9(b). The major axis of the ellipse reflects the most significant vibration. No matter where the transducer mounts, the main vibration vector is the maximum amplitude of the vibration in the motion plane. On the same rotor section, the main vibration vector is not affected by the radial mounting position of the two transducers installed vertically with each other. The full vector spectrum is the synthesis of two homologous transducers’ information. It works as a virtual transducer in the plane to detect the maximum vibration of the interface at each frequency. So its spectral structure can reflect the operation state of the machine comprehensively. The vibration amplitude distribution of the full vector spectrum is in accordance with shaft orbit; Figure 9 shows that the vibration amplitude mainly relates to the double frequency, and the actual situation is also the case.Figure 9
Full vector spectrum and shaft orbit.
(a)
Full vector spectrum
(b)
Shaft orbit


## 4.2. Prediction of Full Vector Amplitude

The fault severity can be determined by analyzing the vibration amplitude. Full vector data fusion was carried out based on the collected vibration signals to extract the full vector amplitude of rotating frequency. The AR(4) data prediction model was established. The time interval is 2 hours. The method is suitable for the long-term monitoring and numerical prediction of stable operating machines. Figure 10 indicates that the predicted value is almost the same as the measured vibration value. The established AR prediction model of full vector information fusion could predict the change of vibration effectively. By using full vector amplitude, the one-sidedness of one single transducer can be overcome.Figure 10
Comparison between predicted and measured values.


## 4.3. Prediction of Full Vector Spectrum

The fault position and property can be determined in accordance with the machine structure by analysis the spectrum structure [28]. To predict the fault property, spectrum structure prediction is a feasible method.Obtaining signal with one single transducer, whether the transducer mounts inx direction or in y direction, the acquired signal is one-sided for fault prediction. The reliability of conclusion concerning fault diagnosis and fault character prediction with such signals cannot be guaranteed. Date fusion must be carried out to enhance the reliability of spectrum. The data fusion is based on the theory of full vector spectrum, and the main vector was used for the prediction of spectrum structure. The main vector sequence reflects the distribution of the vibration energy along the frequency. Compared with the spectrum in one direction only, the full vector spectral structure can reflect the operation state of the machine comprehensively. The result of spectrum structure prediction with AR(13) is shown in Figure 11. The potential fault property can be determined by analyzing the spectrum structure in the prediction domain and reasonable maintenance can be arranged.Figure 11
Full vector spectrum prediction with AR(13).


## 4.4. Equipment Maintenance Using Predicted Feature

Vibration analysis is one of the important methods for fault diagnosis of rotating machinery. By collecting the vibration signal and extracting the eigenvector, the established vibration analysis and prediction system can be used as an important support for equipment health management. By analyzing and predicting the vibration feature, it is possible to evaluate the operating status of the equipment, predict deterioration trends, and identify potential failures that may occur.The established AR prediction model can predict the changes of vibration amplitude and full vector spectrum. By predicting the change of the amplitude, it is possible to determine the deterioration tendency of the equipment and to predict how long the alarm will occur. When the alarm is predicted to occur in the near future, the full vector spectrum prediction model is called to check the amplitude changes at different characteristic frequencies. The characteristic frequency is related to the corresponding equipment failure; for example, rotating frequency is related to the failure of dynamic unbalance. By examining the changes of the amplitude of the characteristic frequency in the spectrum, it is possible to judge out the possible fault to occur. The vibration feature prediction can provide a decision-making basis for the equipment management and maintenance department.


## 4.5. Discussions

In the machine condition monitoring system, the AR model is an effective prediction method. It can be used to predict the change value of the vibration. The structure of the prediction model can be determined by the order number of the AR model. The order of the model was fixed. The coefficients of each order can be determined by the calculation of the sample data. In the actual condition monitoring system, the sampling data was changed as time goes on. The coefficients of each order can be modified to fit in with the latest running state by calculating the sample data in the latest period of time.For the data feature extraction, since the full vector amplitude can reflect the maximum vibration amplitude of the equipment, the prediction model was constructed based on the full vector amplitude. The AR amplitude prediction model can forecast the change trend of the equipment effectively. By analyzing the structure characteristics and the change of full vector spectrum, the fault property and failure position can be predicted. Compared with signals with a single transducer, the full vector spectrum is more effective to reflect the fault condition of the equipment.


## 5. Conclusions

Signal feature extraction and prediction model establishment were two important aspects in the process of the fault prediction. The full vector fusion and fault prediction were studied with the turbine as the research object. The construction of the AR prediction model can effectively solve the problem of vibration data prediction. Since the signal of a single transducer has shortage of one-sidedness for the machine fault diagnosis, the information fusion was accomplished based on the theory of full vector spectrum. The full vector amplitude was extracted to reflect the strength of the equipment. The full vector spectrum was figured out to analyze the fault property. The amplitude and spectrum of fused signal can overcome the shortage of one-sidedness.The established AR prediction model based on the time series full vector amplitude can determine the change of vibration intensity effectively. The established AR spectrum prediction model based on the time series full vector spectrum can predict the change of the structure of the vibration spectrum effectively. According to the failure mechanism of the machine equipment, it is possible to determine the nature and location of the fault by analyzing the predicted frequency spectrum. The study on signal feature extraction based on the theory of full vector spectrum and the construction of prediction model based on the theory of AR model lays a technical foundation for the further research of prognostic and health management.



---

*Source: 6103947-2017-03-12.xml*