2037. A fuzzy diagnosis of multi-fault state based on information fusion from multiple sensors

. This paper presents a fuzzy diagnosis for detecting and distinguishing multi-fault state, the method is constructed on the basis of possibility theory and support vector machines (SVMs) with information fusion from multiple sensors. Non-dimensional symptom parameters (NSPs) are defined to reflect the characteristics of vibration information, and principal component analysis (PCA) is used to evaluate and select sensitive NSPs of each sensor. SVMs are employed to fuse vibration information from different sensors into an effective synthetic symptom parameter (SSP) for increasing diagnostic sensitivity, then the possibility function of the SSP is used to construct a fuzzy diagnosis for fault detection and fault-type identification by possibility theory. Practical examples of diagnosis for a roller bearing used in a test bench are given to show that multi-fault states of bearing can be identified precisely by the proposed method.


Introduction
When diagnosing rotating machinery, the cases that two or more faults had been detected simultaneously often arise [1].Here, the state that two or more faults have occurred simultaneously is called "multi-fault state".The diagnosis of a multi-fault state is an increasingly active research domain [2].However, multiple faults may compensate or cover their features, and vibration information is complex, non-stationary and nonlinear, the multi-fault features are usually immersed in background noise, especially at an early stage.Consequently, fault diagnosis methods working well for single fault state often fail to detect and diagnosis multi-fault state.
Recently, a few methods have been reported to research on multi-fault detection and fault-type identification in the literature.Bachschmid et al. [3] proposed a model-based identification method which was made by a least-squares fitting approach in the frequency domain, by means of the minimization of a multi-dimensional residual between the vibrations in some measuring planes on the machine and the calculated vibrations due to the acting faults.Li et al. [4] introduced an integration of wavelet transform (WT) technique, autoregressive (AR) model and PCA for multi-fault detection.Lal and Tiwari [5] applied forced response information for detecting multi-fault in simple rotor-bearing-coupling systems.Liu et al. [6] proposed a hybrid intelligent multi-fault detection and classification method based on empirical model decomposition (EMD), distance evaluation technique and wavelet SVM with particle swarm optimization (PSO) algorithm.Jiang et al. [7] improved ensemble empirical mode decomposition (EEMD) with multi-wavelet packet for rotating machinery multi-fault diagnosis.
The above studies almost completely utilized the vibration information from single sensor to focus on signal processing.Certainly, the methods of advanced signal processing really improve to extract the multi-fault features.However, single sensor only provides unilateral vibration information, the above diagnosis methods are not enough stability.
It is a fact that multiple sensors can provide abundant vibration information from different angles or force points.Thus, a fuzzy diagnosis method based on information fusion from multiple sensors is proposed for multi-fault detection and diagnosis.In this method, six non-dimensional symptom parameters (NSPs) are defined to reflect the characteristics of vibration information from each sensor, and two sensitive NSPs that contain the most features are selected by PCA.Then SVMs are employed to fuse all sensitive NSPs from multiple sensors into an effective synthetic symptom parameter (SSP), the SSPs of two states, as the object of follow-on process, are used to construct a fuzzy diagnosis system.The rest of this paper is organized as follows: the method based on information fusion from multiple sensors is introduced in Section 2; the fuzzy diagnosis system for multi-fault detection and fault-type identification is proposed in Section 3; the experimental and practical validations are presented in Section 4; finally, discussion and conclusion are given in Section 5.

The method of information fusion from multiple sensors
In the field of condition diagnosis, the excellent symptom parameters (SPs) must be defined to distinguish intelligently states of machinery as precisely as possible, and SPs must be able to reflect sensitively the characteristics of states [1,[8][9][10][11].However, there is not an acceptable method for extracting SPs from signals measured in each state [12].In order to increase the identification sensitivity, a SSP is defined with the function of the optimal classification hyper-plane which is obtained using SVM and all sensitive NSPs from multiple sensors, as shown in Figure 1.First, six NSPs are defined to extract the features of vibration information from each sensor.Second, PCA is used to evaluate and select sensitive NSPs of each sensor.Finally, a SVM is employed to fuse all sensitive NSPs from multiple sensors into an effective SSP.

Non-dimensional symptom parameters for reflecting vibration information
Many SPs have been defined in the pattern recognition field [13,14].In this study, six NSPs in the time domain are defined to reflect vibration information: where ( = 1, 2,…, ) is digital data of vibration signal, is the number of data.̅ and are the mean value and the standard deviation: where, ( = 1, 2,…, ) is the valley data of vibration signal, is the number of data.̅ and are the mean value and the standard deviation: where, ( = 1, 2,…, ) is the valley data of vibration signal, is the number of data.̅ and are the mean value and the standard deviation.

Principal component analysis for selecting sensitive NSPs
Principal component analysis (PCA) is a multivariate statistical-analysis technique, in which a group of correlated variables are transformed into a new group of variables which are uncorrelated or orthogonal to each other [15,16].PCA can be done by eigenvalue decomposition of a data covariance matrix or the singular value decomposition of a data matrix.These decompositions are usually performed after mean centering the data for each attribute.The results of a PCA are usually discussed in terms of component scores and loadings [17].In the last few years, PCA has been applied to process fault diagnosis [18][19][20].
The principle of using the PCA for selecting sensitive SPs is illustrated on the basis of a training data set of n samples (observations) and m variables (SPs).While these data are standardized by subtracting from its mean and dividing by its standard deviation, and the standardized data is denoted by a sample matrix ( ∈ × ), it is easy to obtain the covariance matrix with the eigenvector = ( , , … , ) and the eigenvalue ( > > ⋯ > > 0).Then a new component and its cumulative contribution rate are given as follows: where, = ( , , … , ) is each row vector of , and = 1, 2,…, .In general, when ≥ 85 %, , ,…, are called as principal components which contain most of the information and the discriminatory features.Moreover, the principal components load can express the correlation degree between principal component and original variable .Therefore, the comprehensive load of can be calculated as shown in Eq. (9).It has been proved that the larger the comprehensive load , the higher the sensitivity of the corresponding will be: In this study, the comprehensive load of principal components is used as the evaluation standard of SP's sensitivity, to select two sensitive NSPs which contain most of the vibration information each sensor.

Information fusion from multiple sensors using SVM
SVM is a classifier which uses statistical learning theory to create an optimal classification hyper-plane between the two classes and ensure that the distance between the boundary and the nearest data point in each class is maximized [21].SVM can efficiently perform not only a linear classification but a non-linear classification using kernel function.In general, a classification hyper-plane can be expressed as follows: where, is the sample vector, is the dimensional number of , is a vector represented the weight coefficients of , is a classification threshold.
Recently, almost all of practical applications of SVMs have employed different kernel functions, such as polynomial kernel and radial basis function (RBF) kernel, to achieve linear classifications in high-dimensional feature spaces.However, it is difficult to choose an appropriate kernel function and determine the parameters for a given value of the regularization and kernel parameters.In order to increase the enchantment of SVM, SVM is modified to permit minimum error and relax the condition for the optimal classification hyper-plane, and then soft margin SVM is proposed for a fuzzy inference system [11,22].
In this study, the classification function of soft margin SVM is employed to achieve the optimization of fusing all sensitive NSPs selected from multiple sensors, and generate an effective SSP as the object of follow-on process.Then a SSP is defined as follows: where, , ,…, are the sensitive NSPs selected from multiple sensors by PCA, * and * are the parameters of the optimal hyper-plane by training data.

Fuzzy diagnosis system
In Section 2, the proposed SSP, as the object of follow-on process, can improve excellently the identification sensitivity of two states.Moreover, it is very difficult to generate one SSP that is employed to identify simultaneously all states.However, in the cases of the practical diagnosis, there are many different fault states for the same mechanical equipment.Then possibility theory is introduced to build a sequential diagnosis system for fault detection and fault-type recognition, as shown in Fig. 2. Therefore, this section presents specifically two contents.Firstly, the principle of fuzzy inference is defined on the basis of possibility theory.Secondly, a diagnosis system is established to perform sequentially condition identification of many different faults.

Possibility theory for fuzzy inference
Possibility theory is a mathematical theory for making a fuzzy inference in processing ambiguous information [23,24].Recently, possibility theory has been applied to condition diagnosis in rotating machinery under varying rotating speeds to process the uncertain relationship between the symptoms and fault types [11,[25][26][27].In the paper, for improving the recognition ability of ambiguous state, possibility theory is introduced to make a fuzzy inference for fault detection or fault-type recognition.

Fig. 2. Fuzzy diagnosis system based on possibility theory
For fuzzy inference, the membership function of a SP is necessary [11,[25][26][27].However, there are many methods for building the membership function of a SP.In this paper, a SSP, information fusion from multiple sensors by the proposed method in Section 2, is selected as the features for fault detection or fault-type recognition, and the possibility function of a SSP is regarded as the membership function to make a fuzzy inference.Since it is verified that each SSP follows the Weibull distribution, the possibility function of a SSP ( ( )) can be obtained from probability density functions of the SSP ( ( )) by the following formulae.Moreover, it is verified easily that the maximum value of the possibility function is 1: and can be calculated as follows: where, is the division number of the effective domain of a SSP [ ̅ − 3 , ̅ + 3 ], ̅ and are the mean and the standard deviation of the SSP, respectively., and are the parameters of the shape, scale and location of Weibull distribution.Fig. 3 displays an illustration of the possibility function and the probability density function of a SSP.Then it is easy to obtain the possibility functions of known states and diagnosis state.These known states, such as normal state and various kinds of fault states that have been diagnosed definitely, will be employed to regard the diagnosis standards as model states.Certainly, there are some faults that never occurred, and then these states are called "unknown state".Therefore, there is no way to calculate directly the possibility function of unknown state.In order to extend the diagnostic capacity and enhance the identification sensitivity, an indirect method is proposed to obtain the possibility function of unknown state, as follows: where, is the possibility functions of th known state, respectively, and is the number of all known states.Of course, it is required to define the matching approach of possibility functions between the model states and the diagnosis state for fuzzy inference.However, there are many matching methods defined in the field of fault diagnosis.In [11], the common area of the possibility functions between model state and diagnosis state was proposed as a useful matching method.However, the method must ensure a suitable range of symptom variable.Otherwise, the inference conclusion can be contorted.In [27], the possibility function value of a model state at the mean value of a diagnosis state was introduced to evaluate the degree of which the diagnosis state was the model state.The method that take the experience gained at one unit and popularize it in a whole area can ignore the non-symmetry distribution of symptom variable.Therefore, the paper proposes a novelty matching method based on Euclid close degree.Euclid close degree is one of the methods which express the proximity index of fuzzy sets [28], and can be calculated as follows: where, and are the possibility functions of th model state and diagnosis state, ̅ and are the mean and the standard deviation of symptom variable of diagnosis state, respectively.
Moreover, these close degrees are normalized by: Then the principle of selecting nearest is used to recognize the condition of diagnosis signal.The principle can use "if-then" rule as follows to explain the architecture of fuzzy inference.

Sequential diagnosis approach
Fig. 2 describes the process of the sequential diagnosis system which is proposed in the paper.If there are types of faults that have occurred frequently, along with normal state, the diagnosis system needs to operate steps.
The first step of the diagnosis system is to detect whether mechanical condition is normal.Then the previous state that has been verified as "normal" is employed to confirm the possibility function of normal-model state.In the same way, N types of faults that have occurred frequently and been known definitely are united to establish the possibility function of fault-model state.Moreover, the possibility function of unknown-model state is obtained by the Eq. ( 15).When the information of the diagnosis state is put into the system, the close degrees between three model states and diagnosis state are calculated by the Eq. ( 17), and the fuzzy inference is performed by Rule 1.If the diagnosis state is judged as "normal state", the diagnosis system automatically stops, and the diagnosis data is saved to optimize the possibility function of the future normal-model state.If the unknown state is shown, the system also stops, and we are reminded to investigate the real type of the diagnosis state.Otherwise, there is fault, and the system will come into the next step for fault-type identification.In the same way, each step only focuses to judge whether the diagnosis state is a type of fault from the second step until all fault-model states are performed sequentially.

Experimental verification
To verify the feasibility and effectiveness of the proposed fuzzy diagnosis method for detecting and distinguishing multi-fault state, the fault diagnosis of roller bearing is studied by using the real data obtained from test bench.

Experimental system
Fig. 4 shows the experimental system which is driven by a variable speed 3AC motor with speeds up to 2000 rpm.In order to research fault diagnosis of roller bearing in a real plant, it is equipped with loading setup (Model: RCS2-RA13R) which can transport the power of a fixed or varied load up to 2 T. The bearings (Type: UK204) are utilized, and the specification of the bearings is listed in Table 1.In this work, three accelerometers (Type: PCB MA352A60; Sensitivity: 10 mV/g; Frequency: 5-60000 Hz) were mounted on the three directions (horizontal, vertical and axial direction) of the bearing housing to acquire the vibration signals with a sampling frequency of 50 kHz, and the sampling time was 20 sec.As the rotating machine usually works at a constant speed and a fixed load, the vibration signals of normal state and five faults were measured at a rotation speed of 1000 rpm and 120 kg load.

Diagnosis system building
In order to fully analyze signal feature, these signals measured from different operating condition and different direction were firstly normalized using the following equation for diagnosis: where and were the mean and standard deviation of an original signal .Second, each signal was divided into 120 parts.Each part contained 8192 sampling points for high-pass filtering with 10 kHz cut-off frequency and calculation of six NSPs defined in Section 2.1.Then 120 samples of NSPs were split into two sets: 40 samples for training and the remaining 80 samples for testing.
Since six different operating conditions of roller bearing were performed in the experiment, the diagnosis system was designed to undergo five stages for fault detection and fault-type recognition, as shown in Fig. 6.Then in each stage, 40 samples of NSPs corresponding to two states were used to select two sensitive NSPs from each sensor by PCA as shown in Section 2.2.
As an example, parts of the selection results from six NSPs in horizontal direction are shown in Table 2. Similarly, higher sensitive NSPs were also selected in vertical and axial directions, as shown in Table 3.In each stage, higher sensitive NSPs selected from three directions were used as input layer of soft margin SVM, and two states that would be identified were set to become output layer and associate with labels: = 1 for the goal state and = -1 for the other state.When soft margin SVM was trained, the parameters * and * of the optimal hyper-plane were confirmed, and the definition of the effective SSP in the stage was achieved.Table 4 is the definitions of effective SSPs in each stage.
Then the possibility functions of the model states in each stage were obtained by the Eq. ( 12) and Eq.(15).Thus, the diagnosis system was built to perform sequentially fault detection and fault-type recognition.

Experimental verification
To verify the diagnostic capability of the method proposed in this study, the remaining 80 samples of NSPs from each operating condition and each direction were also split equally into four sets.Four test data sets were labeled by test-1, test-2, test-3, test-4.As an example, when the data of Test-1 were performed in the same procedure, and input into the diagnosis system, the close degrees between three model states of diagnosis system and test state were calculated in each stage, and the results of fuzzy inferences were output, as shown in Table 5.Similarly, the data of Test-2, Test-3, Test-4 were independently used to input into the diagnosis system, and the diagnosis results are displayed in Table 6.

Experimental comparison
To verify the effectiveness of synthetic symptom parameter (SSP) on the basis of vibration information fusion from multiple sensors, vibration information of single sensor was planned to fuse an object of follow-on process.Then higher sensitive NSPs from single direction as shown in Table 3 were respectively used to perform in the same procedure.The diagnosis results on the basis of vibration information from horizontal, vertical or axial direction are shown in Table 7, Table 8, Table 9, respectively.The diagnosis results from Table 6 have shown that all test states have been diagnosed correctly, and close degrees between model states and test states have remained above 0.75.However, the results from Table 7, Table 8 and Table 9 have shown that all fault states have been detected correctly, but more erroneous judges have happened for fault-type recognition, especially multi-fault states.It is very obvious that the diagnostic capability and effectiveness based on information fusion from multiple direction s are better than single direction.
The first author contributed to design and implement the methodology and wrote the manuscript.The other authors carried out the experiments and analyzed the results.

Conclusions
To improve the ambiguous relationship between multi-fault state and the vibration feature, especially at an early stage, and increase the accuracy of fault detection and fault-type recognition, a fuzzy diagnosis method for multi-fault state was presented on the basis of information fusion from multiple sensors, and the effectiveness was demonstrated experimentally.The main conclusions were obtained as follows.
1) Multiple sensors demonstrate a predictable response to complementary information when the sensors are mounted on different locations for diagnosing the same equipment or parts.Vibration information from multiple sensors can process the effective features of multi-fault state, even if at an early stage.
2) It is proven to be effective and improve the robustness of diagnosis system that soft margin SVM is used to fuse vibration information from multiple sensors.
3) The efficiency of fuzzy diagnosis proposed in the paper has been verified by applying it to a practical diagnosis for detecting and distinguishing multi-fault states of bearing.The perfect performance is attributed primarily to effective features from multiple sensors, soft margin SVM's generalization capability, SSP's high sensitivity, and possibility theory for fuzzy inference.
For future works, the method proposed in this paper will be expanded to perform in the vehicle field, the focus is how to extract fault feature under the influence of noise and complex driving conditions, and build the threshold model of fuzzy diagnosis.

Fig. 1 .
Fig. 1.Flow chart of the proposed method

Fig. 3 .
Fig. 3. Possibility function and probability density function of a SSP It is obvious that expresses the possibility function of unknown state except all known states.Then is regarded as the diagnosis standard of unknown state in the field of fuzzy inference.Of course, it is required to define the matching approach of possibility functions between the model states and the diagnosis state for fuzzy inference.However, there are many matching methods defined in the field of fault diagnosis.In[11], the common area of the possibility functions between model state and diagnosis state was proposed as a useful matching method.However, the method must ensure a suitable range of symptom variable.Otherwise, the inference conclusion can be contorted.In[27], the possibility function value of a model state at the mean value of a diagnosis state was introduced to evaluate the degree of which the diagnosis state was the model state.The method that take the experience gained at one unit and popularize it in a whole area can ignore the non-symmetry distribution of symptom variable.Therefore, the paper proposes a novelty matching method based on Euclid close degree.Euclid close degree is one of the methods which express the proximity index of fuzzy sets[28], and can be calculated as follows: Rule 1: if ( , ) = max ( , )| = 1, 2, … , , … , , then ∈ .where ( , ) is the close degree between model state and diagnosis state.FUZZY DIAGNOSIS OF MULTI-FAULT STATE BASED ON INFORMATION FUSION FROM MULTIPLE SENSORS.HONGTAO XUE, ZHONGXING LI, YONG LI, HONG JIANG, PENG CHEN possibility functions of the given model state.and express the conditions of diagnosis signal and model signal .

Fig. 4 . 5 .
Fig. 4. The rotating machine for fault diagnosis In the case of the practical diagnosis, three single-fault states, namely, the outer race defect (O), the inner race defect (I), the rolling element defect (R), and two multi-fault states, such as the

Fig. 6 .
Fig. 6.Flowchart of sequential diagnosis based on the test experiment

Table 1 .
Specification of the bearing

Table 2 .
Higher sensitive NSPs selected by PCA in horizontal direction

Table 3 .
Higher sensitive NSPs from three directions by PCA

Table 4 .
Synthetic symptom parameters (SSPs) in each stage

Table 5 .
Close degrees in each stage and diagnosis results of Test-1

Table 6 .
Diagnosis results by the proposed method

Table 7 .
Diagnosis results based of vibration information from horizontal direction

Table 8 .
Diagnosis results based of vibration information from vertical direction

Table 9 .
Diagnosis results based of vibration information from axial direction