# Full-Vector Signal Acquisition and Information Fusion for the Fault Prediction

Authors: Lei Chen; Jie Han; Wenping Lei; Yongxiang Cui; Zhenhong Guan
Journal: International Journal of Rotating Machinery
 (2016)

---

## Abstract

Fault prediction is the key technology of the predictive maintenance. Currently, researches on fault prediction are mainly focused on the evaluation of the intensities of the failure and the remaining life of the machine. There is lack of methods on the prediction of fault locations and fault characters. To satisfy the requirement of the prediction of the fault characters, the data acquisition and fusion strategies were studied. Firstly, the traditional vibration measurement mechanism and its disadvantages were presented. Then, the full-vector data acquisition and fusion model were proposed. After that, the sampling procedure and information fusion algorithm were analyzed. At last, the fault prediction method based on full-vector spectrum was proposed. The methodology is that of Dr. Bently and Dr. Muszynska. On the basis of this methodology, the application study has been carried out. The uncertainty of the spectrum structure can be eliminated by the designed data acquisition and fusion method. The reliability of the diagnosis on fault character was improved. The study on full-vector data acquisition system laid the technical foundation for the prediction and diagnosis research of the fault characters.


---

## Body

## 1. Introduction

As the key technology of the predictive maintenance [1], fault prediction is utilized to predict the operating condition of the machine and give an early warning of the failure which may occur in the future. Researchers have conducted lots of work on the aspects of machine fault diagnosis [2–6], operation status prediction [7–10], and feature extraction and information fusion [11–13], and the application of fault prediction has been proved to be favorable to instruct production dispatch and equipment maintenance. However, existing researches on mechanical fault prediction mainly focus on the forecasting of single specific fault trend and evaluating equipment’s residual life. There are few studies and reports about the aspect of fault character, failure type, and failure position [14, 15].To analyze the fault source and predict the fault character, failure type, and failure position of the equipment, one of the feasible means is to predict the trend of frequency spectrum. The failure type which may occur can be determined by analyzing the predicted spectrum structure [16, 17]. To ensure the reliability of the prediction, the acquired signal should more faithfully and comprehensively reflect the operation status of the machine. At the same interface of rotating machine, the signal spectrum of different vibration will be obtained because of the different installation direction of the transducers. It is obvious that the information collected by a single transducer is not comprehensive [16]. Including the Chinese Power Utility Industry, 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. Bently et al. [18]. Dr. Bently’s classic work methodically and accurately explains in great detail the use of diagnostic equipment, theory, and applications. On the basis of these methodologies, the application has been carried out in this study. At the same test position, spectrum structure with one installation direction may be different from that of another direction. An error diagnosis may be made by analyzing such a signal spectrum. In order to make a comprehensive prediction and feasible diagnosis, making information fusion becomes necessary. The fusion signal should more faithfully and comprehensively reflect the machine operating condition. The fusion signal’s spectrum is not concerned with the installation direction of the transducer. The failure type can be predicted by analyzing such a signal’s spectrum trend. Therefore, this paper will take the rotating machine as the research object. Firstly, the data acquisition and signal fusion model for the fault character prediction will be analyzed. Then, the vibration signal acquisition system based on full-vector spectrum will be designed. Lastly, the signal fusion algorithm will be further studied.


## 2. Vibration Mechanism and Signal Monitoring

Generally, the rotating machine has complex structure. For the convenience of analysis, the mechanical model of the machine rotor was simplified as a disc with massm and an elastic shaft with no mass, and the shaft ends are supported by a rigid bearing and a bearing base [19–21], as shown in Figure 1.Figure 1
Jeffcott rotor.Assume that the two support points of the rotor are in the same horizontal line and the disc is located in the center of the two fulcrums of the rotor. When the rotor is rotating, it generates dynamic torsion with diameterr due to the centrifugal inertia force. There are two kinds of rotations at this time: (1) rotation of the disc around the axis AO′B and (2) rotation of the plane of the AO′B and the AOB of the bearing line around the AB axis [16, 22].Assuming the mass of the disc ism, the elastic recovery force of the rotor is F:(1)F=-kr,wherek is the rotor stiffness coefficient, r=OO′. The differential equation of motion of the disc is(2)mx¨=Fx=-kx,my¨=Fy=-ky,Makeω2=km,x¨+ω2x=0,y¨+ω2y=0.Its solution can be written as(3)x=Xcos⁡ωt+φxy=Ycos⁡ωt+φy,where X stands for the amplitude of the motion of the disc center in the direction of x and Y stands for the amplitude of the motion of the disc center in the direction of y. φx and φy stand for the phase angles of the corresponding direction and t stands for the time. The rotor vibration of x- and y-directions of a stationary rotating machinery is a combination of motion orbits for some harmonic waves wi (i=1,2,…).Generally, the vibration of the rotating machine is measured by noncontact eddy current transducers. The transducers are installed with theX-Y layout or V installation as is shown in Figure 2. The measured vibration signal of x or y transducer is a combination of motion orbits for some harmonic waves wi (i=1,2,…). For a certain harmonic i, suppose(4)Xci=Xicos⁡∅xi,Xsi=Xisin⁡∅xi,Yci=Yicos⁡∅yi,Ysi=Yisin⁡∅yi,where Xi stands for the amplitude of vibration signal of the harmonic i 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 the direction x of harmonic i when t=0; Yci stands for the amplitude of vibration in the 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.Figure 2
Layout of transducers.Movement equation ofx- and y-direction can be expressed as (5)x=∑i=1∞Xicos⁡wit+∅xi=Re⁡∑i=1∞X-iejwit=Re⁡∑i=1∞xci+jxsiejwit,y=∑i=1∞Yicos⁡wit+∅yi=Re⁡∑i=1∞Y-iejwit=Re⁡∑i=1∞yci+jysiejwit.In the above formula,X-i and Y-i represent the complex amplitudes of the vibration signals in the directions x, y of the harmonic i (i=1,2,3,…); the trajectories of each harmonic are an ellipse; that is, (6)Yci2+Ysi2xi2+Xci2+Xsi2yi2-2XciYsi+XsiYcixiyi=XsiYsi-XciYci2.Suppose there exists a vibration signal with three harmonics. The angular frequencies of the harmonics arew3, w2, and w1. The harmonics have the same length of the major axis A. The angles between the elliptical major axis and the x-direction are π/2, 3π/4, and 0, as shown in Figure 3. Figure 3 indicates that when measuring the same signal with transducer in different directions, different frequency spectrums will be obtained. That is to say, whether the transducer mounts in x-direction or in y-direction, the acquired signal in one signal transducer is one-sided for fault diagnosis. The reliability of conclusion concerning fault diagnosis and fault character prediction with such signals cannot be guaranteed. In order to obtain the reliable fault diagnosis and prediction, it is required to overcome the one-sidedness of the vibration information obtained by a single transducer. One-sided information can be overcome by key-phasor [18]. Information fusion of the transducers in data layer is a feasible mean. In this study, full-vector signal is fused in the process of vibration data acquisition, and the full-vector spectrum will be extracted for the purpose of fault diagnosis and prediction.Figure 3
Spectrum structure error of transducers in different directions.


## 3. Full-Vector Data Acquisition and Information Fusion

### 3.1. Full-Vector Data Acquisition

As is shown in Figure4, transducers of x- and y-direction sampled integer periodic vibration signal synchronously. The vibration acquired by the transducers of x- and y-direction is fused based on full-vector spectrum in the data preprocessing stage, and a set of full-vector signals is obtained [18]. In order to obtain the full-vector phase spectrum of the vibration, a key phase trigger signal is needed. A key phase transducer is used to provide reference position of vibration phase and measure rotational speed. The data sampling process is controlled by a field-programmable gate array (FPGA). Application running on the central processing unit (CPU) completes the full-vector fusion operation and obtains feature vectors for predicting and diagnosing.Figure 4
Structure of full-vector data acquisition.


### 3.2. Control Logic

The sampling process was controlled by a FPGA. The control logic is shown in Figure5.Figure 5
Procedure of full-vector acquisition.FPGA receives sampling parameters and instructions such as sampling frequency, sampling length, and sampling start from the CPU and converts them into the sampling control logic of analog-to-digital (AD) converter. As for the rotating machine, due to the importance of speed and phase for the fault diagnosis, the sampling is triggered by an external key. When the falling edge of the trigger signal comes, integrated periodic signal was acquired synchronously [23]. As is shown in Figure 6, the key phase transducer generates a level jump in every turn of the rotor. The interval between the two jumps is the rotation period. After the sampling instruction is issued, current speed is calculated in the FPGA and sampling rate is set according to the speed. AD conversion starts when the first pulse falling edge comes after setting the sampling parameters. Continuous 2N sampling points are acquired for each channel. The signals of x- and y-channel are acquired synchronously. Each sampling signal contains W cycles and each cycle contains 2M sampling points. Then, the sampling length 2N=W×2M. The synthetic signal of X and Y sensor is a motion ellipse. The major axis of the ellipse, called main vector of the full-vector spectrum, 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. The vibration signal can be considered as superposition of a series of harmonic signals. The main vibration vectors of each harmonic can be obtained by data fusion, and then the full-vector spectrum is obtained. The obtained full spectrum is good for fault diagnosis and prediction.Figure 6
Diagram of synchronous integrated period sampling.


### 3.3. Fusion Algorithm

According to the theory of full-vector spectrum, the main vector of vibration harmonic can be simplified as a simple calculation by FFT. Assume thatxi and yi are the signals acquired by the transducers x and y. Here, i=0,1,2,…,N-1, where N is the sampling points. A complex sequence zi can be constituted by 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, (7) can be obtained. Hence,(7)RLk=12NZk+ZN-kRSk=12NZk-ZN-kk=0,1,2,…,N2-1,whereRLk and RSk are the main and assistant vibration vectors. Since the synthetic signal of X and Y transducer at each harmonic is a motion ellipse, the major axis of the ellipse reflects the most significant vibration. RLk is half major axis of the ellipse at the frequency point k. RSk is half minor axis of the ellipse at the frequency point k. RLk is the maximum amplitude of the vibration in the motion plane. On the same rotor section, RLk and RSk are not affected by the radial mounting position of the two transducers installed vertically with each other. The full-vector spectrum, spectrum of RLk, is the synthesis of two homologous transducers’ information. Its spectral structure can reflect the operation state of the machine comprehensively. Fusion calculation in the data acquisition system can be carried out according to Figure 7.Figure 7
Full-vector fusion calculation.


## 3.1. Full-Vector Data Acquisition

As is shown in Figure4, transducers of x- and y-direction sampled integer periodic vibration signal synchronously. The vibration acquired by the transducers of x- and y-direction is fused based on full-vector spectrum in the data preprocessing stage, and a set of full-vector signals is obtained [18]. In order to obtain the full-vector phase spectrum of the vibration, a key phase trigger signal is needed. A key phase transducer is used to provide reference position of vibration phase and measure rotational speed. The data sampling process is controlled by a field-programmable gate array (FPGA). Application running on the central processing unit (CPU) completes the full-vector fusion operation and obtains feature vectors for predicting and diagnosing.Figure 4
Structure of full-vector data acquisition.


## 3.2. Control Logic

The sampling process was controlled by a FPGA. The control logic is shown in Figure5.Figure 5
Procedure of full-vector acquisition.FPGA receives sampling parameters and instructions such as sampling frequency, sampling length, and sampling start from the CPU and converts them into the sampling control logic of analog-to-digital (AD) converter. As for the rotating machine, due to the importance of speed and phase for the fault diagnosis, the sampling is triggered by an external key. When the falling edge of the trigger signal comes, integrated periodic signal was acquired synchronously [23]. As is shown in Figure 6, the key phase transducer generates a level jump in every turn of the rotor. The interval between the two jumps is the rotation period. After the sampling instruction is issued, current speed is calculated in the FPGA and sampling rate is set according to the speed. AD conversion starts when the first pulse falling edge comes after setting the sampling parameters. Continuous 2N sampling points are acquired for each channel. The signals of x- and y-channel are acquired synchronously. Each sampling signal contains W cycles and each cycle contains 2M sampling points. Then, the sampling length 2N=W×2M. The synthetic signal of X and Y sensor is a motion ellipse. The major axis of the ellipse, called main vector of the full-vector spectrum, 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. The vibration signal can be considered as superposition of a series of harmonic signals. The main vibration vectors of each harmonic can be obtained by data fusion, and then the full-vector spectrum is obtained. The obtained full spectrum is good for fault diagnosis and prediction.Figure 6
Diagram of synchronous integrated period sampling.


## 3.3. Fusion Algorithm

According to the theory of full-vector spectrum, the main vector of vibration harmonic can be simplified as a simple calculation by FFT. Assume thatxi and yi are the signals acquired by the transducers x and y. Here, i=0,1,2,…,N-1, where N is the sampling points. A complex sequence zi can be constituted by 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, (7) can be obtained. Hence,(7)RLk=12NZk+ZN-kRSk=12NZk-ZN-kk=0,1,2,…,N2-1,whereRLk and RSk are the main and assistant vibration vectors. Since the synthetic signal of X and Y transducer at each harmonic is a motion ellipse, the major axis of the ellipse reflects the most significant vibration. RLk is half major axis of the ellipse at the frequency point k. RSk is half minor axis of the ellipse at the frequency point k. RLk is the maximum amplitude of the vibration in the motion plane. On the same rotor section, RLk and RSk are not affected by the radial mounting position of the two transducers installed vertically with each other. The full-vector spectrum, spectrum of RLk, is the synthesis of two homologous transducers’ information. Its spectral structure can reflect the operation state of the machine comprehensively. Fusion calculation in the data acquisition system can be carried out according to Figure 7.Figure 7
Full-vector fusion calculation.


## 4. Experiment and Discussion

### 4.1. Experiment

To verify the correctness of the data acquisition system and fusion theory, an online monitoring system was installed in a thermal power plant. TheBently Nevada 3300 XL series of sensor systems were used in the system. And theeM3000 series data acquisition system was provided byZhengzhou Expert P.D.E. Co., Ltd. Arrangement of measuring points on the feed water pump of steam turbine is shown in Figure 8. Two eddy current transducers are installed with V layout at the position of input shaft bearing of the turbine. Vibration signals were acquired by synchronous integrated period sampling. The sampling frequency is 2000 Sa/s, sampling length is 1024 points, and the rotor speed is 2930 r/min. The signals of x- and y-direction of the same cross section of an axis of the feed pump for a turbine are shown in Figure 9. The amplitude of first-order frequency, the amplitude of second-order frequency, and possible conclusions are shown in Table 1. Although the signal vibration amplitudes are all within the normal limits, there are obvious differences between the x and y signals in both time and frequency domain. According to the signal 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 signal in the y-direction, it seems that the vibration relates to the rotating frequency; the fault may be related to the unbalance. Different fault conclusions may be drawn on the machine by different signals. The synthetic signal of x and y transducer at each harmonic is a motion ellipse. 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 [18]. It works as a virtual transducer in the plane to detect the maximum vibration of the interface at each harmonic. So its spectral structure can reflect the operation state of the machine comprehensively.Table 1
Difference of thex-direction and y-direction and fusion signals.
Signal
Amp. of 1st order (μm)
Amp. of 2nd order (μm)
Possible conclusions
x-direction
2.61
3.19
Fault with double frequency
y-direction
5.42
2.41
Fault with rotating frequency
Fusion signal
5.52
3.97
Fault with rotating frequencyFigure 8
Machine structure and point position.Figure 9
(a) Vibration wave ofx-direction; (b) spectrum of x-direction; (c) vibration wave of y-direction; (d) spectrum of y-direction; (e) full-vector spectrum.
(a)
(b)
(c)
(d)
(e)R
L
k is half major axis of the ellipse at the frequency point k. The spectrum of RLk, called full-vector spectrum, is the synthesis of two homologous transducers’ information and is not affected by the radial mounting position of the two transducers installed vertically with each other. Comparing with the spectrum in one direction only, its spectral structure can reflect the operation state of the machine comprehensively. Autoregressive (AR) model can be used for time series vibration analysis. Figure 10 is a series of predicted full-vector spectrums based on time series signals. The order of the AR model is 13. The signals were acquired from a steam turbine. High pressure rotor speed of the turbine is 3000 r/min; sampling frequency is 1600 Hz; sampling points are 1024. The AR model was constructed by 67 sets of vibration data. Since change trend of vibration amplitude of each harmonic is different, the AR models were constructed for the time series of each harmonic. The spectrum predicted model is an AR model array. Since the fault characteristic is concerned with the spectrum structure, the characteristic and position of machine fault which may occur in the future can be determined by analysis of the predicted full-vector spectrum. By this means, one may now determine a potential fault such as cracked rotor, bad bearing, and seal instability and make a reasonable arrangement of equipment maintenance.Figure 10
Results of full-vector spectrum prediction.


### 4.2. Discussion

The main vector sequenceRLk (k=0,1,2,…,N/2-1) reflects the distribution of the vibration energy along the frequency. Comparing with the spectrum in one direction only, its spectral structure can reflect the operation state of the machine comprehensively. At the same test position, spectrum structure with one installation direction may be different from that of another direction. On the same rotor section, the full-vector spectra’s structure is not affected by the radial mounting position of the two transducers installed vertically with each other. The full-vector spectrum is the true reflection of the vibration characteristics. The characteristic and position of machine fault which may occur in the future can be determined by analysis of the predicted full-vector spectrum.


## 4.1. Experiment

To verify the correctness of the data acquisition system and fusion theory, an online monitoring system was installed in a thermal power plant. TheBently Nevada 3300 XL series of sensor systems were used in the system. And theeM3000 series data acquisition system was provided byZhengzhou Expert P.D.E. Co., Ltd. Arrangement of measuring points on the feed water pump of steam turbine is shown in Figure 8. Two eddy current transducers are installed with V layout at the position of input shaft bearing of the turbine. Vibration signals were acquired by synchronous integrated period sampling. The sampling frequency is 2000 Sa/s, sampling length is 1024 points, and the rotor speed is 2930 r/min. The signals of x- and y-direction of the same cross section of an axis of the feed pump for a turbine are shown in Figure 9. The amplitude of first-order frequency, the amplitude of second-order frequency, and possible conclusions are shown in Table 1. Although the signal vibration amplitudes are all within the normal limits, there are obvious differences between the x and y signals in both time and frequency domain. According to the signal 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 signal in the y-direction, it seems that the vibration relates to the rotating frequency; the fault may be related to the unbalance. Different fault conclusions may be drawn on the machine by different signals. The synthetic signal of x and y transducer at each harmonic is a motion ellipse. 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 [18]. It works as a virtual transducer in the plane to detect the maximum vibration of the interface at each harmonic. So its spectral structure can reflect the operation state of the machine comprehensively.Table 1
Difference of thex-direction and y-direction and fusion signals.
Signal
Amp. of 1st order (μm)
Amp. of 2nd order (μm)
Possible conclusions
x-direction
2.61
3.19
Fault with double frequency
y-direction
5.42
2.41
Fault with rotating frequency
Fusion signal
5.52
3.97
Fault with rotating frequencyFigure 8
Machine structure and point position.Figure 9
(a) Vibration wave ofx-direction; (b) spectrum of x-direction; (c) vibration wave of y-direction; (d) spectrum of y-direction; (e) full-vector spectrum.
(a)
(b)
(c)
(d)
(e)R
L
k is half major axis of the ellipse at the frequency point k. The spectrum of RLk, called full-vector spectrum, is the synthesis of two homologous transducers’ information and is not affected by the radial mounting position of the two transducers installed vertically with each other. Comparing with the spectrum in one direction only, its spectral structure can reflect the operation state of the machine comprehensively. Autoregressive (AR) model can be used for time series vibration analysis. Figure 10 is a series of predicted full-vector spectrums based on time series signals. The order of the AR model is 13. The signals were acquired from a steam turbine. High pressure rotor speed of the turbine is 3000 r/min; sampling frequency is 1600 Hz; sampling points are 1024. The AR model was constructed by 67 sets of vibration data. Since change trend of vibration amplitude of each harmonic is different, the AR models were constructed for the time series of each harmonic. The spectrum predicted model is an AR model array. Since the fault characteristic is concerned with the spectrum structure, the characteristic and position of machine fault which may occur in the future can be determined by analysis of the predicted full-vector spectrum. By this means, one may now determine a potential fault such as cracked rotor, bad bearing, and seal instability and make a reasonable arrangement of equipment maintenance.Figure 10
Results of full-vector spectrum prediction.


## 4.2. Discussion

The main vector sequenceRLk (k=0,1,2,…,N/2-1) reflects the distribution of the vibration energy along the frequency. Comparing with the spectrum in one direction only, its spectral structure can reflect the operation state of the machine comprehensively. At the same test position, spectrum structure with one installation direction may be different from that of another direction. On the same rotor section, the full-vector spectra’s structure is not affected by the radial mounting position of the two transducers installed vertically with each other. The full-vector spectrum is the true reflection of the vibration characteristics. The characteristic and position of machine fault which may occur in the future can be determined by analysis of the predicted full-vector spectrum.


## 5. Conclusions

The machine running condition can be predicted by means of fault prediction. Usually, the fault prediction uses the data-driven model. Traditionally, the signals are acquired by one or more transducers. 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. Bently et al. [18]. However, researches on mechanical fault prediction, currently, are still mainly focused on the forecasting of single specific fault trend and evaluating equipment’s residual life. In this study, the application has been carried out on the basis of the methodologies proposed by Dr. Bently and Dr. Muszynska. The spectrum based on one transducer has disadvantages of uncertainty and one-sidedness. The proposed full-vector data acquisition and information fusion method for fault prediction effectively overcame the above problems. In theory, the data acquisition system, on the one hand, can eliminate the uncertainty of the spectrum structure caused by a single vibration signal; on the other hand, data fusion in the data acquisition layer can also provide a feasible method for the establishment of the full-vector prediction model. The study on full-vector data acquisition system lays the technical foundation for the further research of the prediction and diagnosis of fault characters.



---

*Source: 5980802-2016-03-20.xml*