Next Article in Journal
Optimal Scheduling of Industrial Task-Continuous Load Management for Smart Power Utilization
Previous Article in Journal
The Reduction of Vertical Interchannel Crosstalk: The Analysis of Localisation Thresholds for Natural Sound Sources
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Engine Fault Diagnosis Method Based on Multi-Sensor Data Fusion

School of Electronics and Information, Northwestern Polytechnical University, Xi’an 710072, Shaanxi, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2017, 7(3), 280; https://doi.org/10.3390/app7030280
Submission received: 29 January 2017 / Revised: 1 March 2017 / Accepted: 9 March 2017 / Published: 14 March 2017

Abstract

:
Fault diagnosis is an important research direction in modern industry. In this paper, a new fault diagnosis method based on multi-sensor data fusion is proposed, in which the Dempster–Shafer (D–S) evidence theory is employed to model the uncertainty. Firstly, Gaussian types of fault models and test models are established by observations of sensors. After the models are determined, the intersection area between test model and fault models is transformed into a set of BPAs (basic probability assignments), and a weighted average combination method is used to combine the obtained BPAs. Finally, through some given decision making rules, diagnostic results can be obtained. The proposed method in this paper is tested by the Iris data set and actual measurement data of the motor rotor, which verifies the effectiveness of the proposed method.

1. Introduction

Engine fault diagnosis, one of the aspects of fault diagnosis [1,2], which is useful for engine design [3,4], can make the maintenance personnel quickly grasp the engine fault type and make the corresponding maintenance measures. At present, there are many methods for fault diagnosis. For example, the methods based on expert systems [5,6,7], which utilize the experience accumulated by domain experts in long-term practice. In these kinds of methods, a computer program is designed to simulate the reasoning and decision-making process of human experts to diagnose faults. Other fault diagnosis methods, such as machine learning [8,9] and signal processing [10,11,12], where there is a comparison of frequency analysis and SVM (support vector machine) classifiers in [11,12], are also widely used in practical applications [13,14]. There are also fault diagnosis methods in multi-sensor networks. As another efficient method to fault diagnosis, the residual-based unsupervised fault diagnosis method addresses uncertainty in the form of error bars around the soft computing models extracted from data [15,16,17]. Multi-sensor data fusion fault diagnosis [18,19], which is widely used in sensor design [20,21], as a data-driven approach, has been paid more and more attention. The advantages of multi-sensor data fusion for fault diagnosis can be concluded as follows:
  • Compared with single source data, multi-source information fusion or multi-sensor information fusion can improve the accuracy and promptness of fault diagnosis.
  • Multi-sensor data fusion method has the ability to deal with the uncertainty of data in fault diagnosis so as to enhance the credibility of diagnostic results.
A number of theories have been developed on multi-sensor data fusion, such as the Dempster–Shafer (D–S) evidence theory [22,23,24,25,26], fuzzy set theory [27], rough sets theory [28,29], Z numbers [30,31], and D numbers [32,33,34,35]. The D–S evidence theory, first proposed by Dempster [22] and then developed by Shafer [23], is able to deal with uncertain information without a prior probability. It has been applied to wide areas, such as decision making [36,37,38], pattern recognition [39,40,41] and risk assessment [42,43,44,45]. In the D–S evidence theory, mass function (or basic probability assignment (BPA)), belief function and plausibility function are defined to express uncertain information. Many methods can be used to integrate multi-source information [46,47,48], while Dempster’s combination rule, which is effective at addressing uncertainty information, can be used to fuse multi-source information in uncertain environments to make a decision.
Much research have been carried out in fault diagnosis based on multi-sensor data fusion [31,49,50,51,52]. Generally, the overall process of fault diagnosis based on multi-sensor data fusion is as follows. Firstly, prior to diagnosis, fault models are established using features extracted from known faulty equipment. Secondly, test models are established using the extracted features from a machine to be detected. Thirdly, BPAs are obtained by computing the BPA determination algorithm using the obtained fault models and test models. Finally, final BPA is obtained by using a combination rule to fuse all of the obtained BPAs. Decision making is carried out through the obtained final BPA. In this process, the key issues are the reasonable establishment of models and the reasonable determination of BPA. However, in previous studies, they are not well solved, especially considering the uncertainty. For example, in [31], the authors take ordinates of intersection points of functions as BPA, while it is not enough to express the uncertain information in fault models due to inadequate consideration on model variance information. Xu [53] utilized the relationship between test data and the Gaussian distribution model generated by training data to determine BPA, while using single measurement data as test data may fail in addressing uncertainty when the fluctuation of sensor data is large. Hégarat-Mascle et al. [54] and Salzenstein et al. [55] modeled the knowledge provided by each information source using probability density functions, while the mass value associated with each compound hypothesis is obtained by subtracting the mass values of the involved individual hypotheses. In this paper, a new fault diagnosis method based on multi-sensor data fusion is proposed. In the step of modeling, for both fault models and test models, considering the Gaussian type interference of mechanical noise and electromagnetic waves in sensor working environments, the Gaussian type membership function is established as the model. The proposed modeling method matches measured data of sensor well. In the step of BPA determination, the intersection area between model curves under the identical attribute is transformed into a set of BPA. The proposed BPA determination method, which fully utilize the variance information of models, can better express uncertain information. In the step of BPA combination and decision making, a weighted average approach for evidence combination [56] is used to combine the obtained BPAs and multiple criteria are established for decision making. The proposed method establishes reasonable fault models and test models based on the complex working environment of sensors and takes into account the variance information of the model, which can better express the uncertain information in the system and improve the reliability of the fault diagnosis.
The paper is organized as follows: in Section 2, the preliminaries are briefly introduced. In Section 3, a new fault diagnosis method based on multi-sensor data fusion is proposed. The proposed method is used for Iris data set classification and engine fault diagnosis in Section 4. Conclusions are made in Section 5.

2. Preliminaries

2.1. Dempster–Shafer Evidence Theory

The D–S evidence theory [22,23], as introduced by Dempster and then developed by Shafer, has emerged from their works on statistical inference and uncertain reasoning. This theory is widely used in decision making [57,58,59,60,61] and information fusion [62,63,64,65]. In this part, a few concepts about D–S evidence theory are given.
Let Θ be a finite nonempty set of mutually exclusive alternatives and Θ is called the frame of discernment
Θ = θ 1 , θ 2 , θ i , , θ N .
The power set of Θ is indicated by 2 Θ , namely:
2 Θ = , θ 1 , θ 2 . . . , θ N , θ 1 , θ 2 , . . . , θ 1 , . . . , θ i , . . . , θ 1 , . . . , θ N .
For a frame of discernment Θ, a basic probability assignment (BPA) is a mapping m from 2 Θ to [0,1], formally defined by
m : 2 Θ 0 , 1 ,
which satisfies the following condition:
m = 0 ,
A 2 Θ m A = 1 .
The mass m A represents how strongly the evidence support A. A is called a focal element when m A > 0 A 2 Θ .
Two BPAs m 1 and m 2 can be combined with Dempster’s combination rule as follows:
m A = 0 , A = , 1 1 K B C = A m 1 B m 2 C , A ,
K = B C = m 1 B m 2 C ,
where K is a normalization constant, reflecting the conflict between the two BPAs m 1 and m 2 .
For the sake of understanding, an illustrative example is given as follows to show how to combine two BPAs according to Equations (6) and (7).
Example 1.
Assuming m 1 ( { a } , { a , b } , { c } ) = { 0 . 6 , 0 . 3 , 0 . 1 } , m 2 ( { a } , { b } , { c } ) = { 0 . 7 , 0 . 2 , 0 . 1 } , computational method of K and m ( { a } , { b } , { c } ) is as follows:
K = m 1 ( { a } ) × m 2 ( { b } ) + m 1 ( { a } ) × m 2 ( { c } ) + m 1 ( { a , b } ) × m 2 ( { c } ) + m 1 ( { c } ) × m 2 ( { a } ) + m 1 ( { c } ) × m 2 ( { b } )     = 0 . 6 × 0 . 2 + 0 . 6 × 0 . 1 + 0 . 3 × 0 . 1 + 0 . 1 × 0 . 7 + 0 . 1 × 0 . 2 = 0 . 3 m ( { a } ) = m 1 ( a ) × m 2 ( a ) + m 1 ( a , b ) × m 2 ( a ) 1 K = 0 . 6 × 0 . 7 + 0 . 3 × 0 . 7 1 0 . 3 = 0 . 9 m ( { b } ) = m 1 ( a , b ) × m 2 ( b ) 1 K = 0 . 3 × 0 . 2 1 0 . 3 = 0 . 09 m ( { c } ) = m 1 ( c ) × m 2 ( c ) 1 K = 0 . 1 × 0 . 1 1 0 . 3 = 0 . 01 .

2.2. Evidence Distance

Jousselme [66] proposed the concept of evidence distance to measure the distance of BPAs. Let m 1 and m 2 be two BPAs on the frame of discernment Θ, containing N mutually exclusive and exhaustive hypotheses. The distance between m 1 and m 2 is defined as
d B P A m 1 , m 2 = 1 2 m 1 m 2 T D ̲ ̲ m 1 m 2 ,
where m 1 and m 2 are the BPAs according to Equations (4) and (5) and D ̲ ̲ is a 2 N × 2 N matrix whose elements are D A , B = A B / A B ( A , B 2 Θ ) .

2.3. Weighted Average Approach for Evidence Combination

In the D–S evidence theory, Dempster’s combination rule is used to fuse certain multi-source information to make a decision. However, if there is a high level of conflict between evidence, the combination rule will fail [67]. Many methods were proposed to solve this problem [66,68,69,70]. Among them, Murphy [71] proposed a simple average approach to solve the problem. However, in Murphy’s approach, all bodies of evidence are equally important. Hence, a weighted average approach [56] was proposed to improve Murphy’s approach. The weighted average approach is given as follows.
Suppose the distance between two BPAs m i and m j is calculated and denoted as d B P A ( m i , m j ) . The similarity measure S i m i j between m i and m j is defined as:
S i m m i , m j = 1 d B P A m i , m j .
If n BPAs need to be combined, the support degree of BPA m i is defined as:
S u p m i = j = 1 , j i . n S i m m i , m j ,
Then, the credibility degree C r d i of the BPA m i is calculated as follows:
C r d i = S u p m i i = 1 n S u p m i .
The credibility degree is actually a weight that shows the relative importance of the collected evidence. According to C r d i , the weighted average of the evidence MAE is obtained:
m ¯ = M A E m = i = 1 n C r d i × m i .
Finally, the weighted average evidence m ¯ is combined n-1 times by using Dempster’s combination rule:
m f = m ¯ m ¯ . . . m ¯ n 1 t i m e s .
The obtained m f is the combination result of the n BPAs.

3. The Proposed Method for Fault Diagnosis

In this paper, a new multi-sensor data fusion method is proposed for fault diagnosis, and the detailed description of the proposed method is shown in Figure 1. One or more attribute features, which can be used for fault diagnosis, can be extracted from each sensor signal. Firstly, in step 1 of Figure 1, Gaussian types of fault models are established using the extracted attribute features from the signal of sensors in each known type of fault equipment. In step 2, Gaussian types of test models are generated as well using the attribute features extracted from sensor signals of the equipment. In step 3, the produced fault models and test models are matched to produce the BPAs. Finally, in step 4, the weighted average approach is used to combine the produced BPAs. Through the decision making rules, the diagnostic results that the equipment is working properly or which type of fault the equipment has is diagnosed. If the equipment has a type of fault, corresponding intervention measures will be implemented.

3.1. Fault Model Modeling

In this section, we give the method to determine the fault model. Due to the change of time and environment of the sensor measurements, the measured data have a certain degree of fuzziness. The membership function can be used to represent the fault features extracted from the observed data.
Suppose X is a sample space of a type of attribute feature variable extracted from the sensor—for example, the amplitude of fundamental harmonic in the frequency spectrum of sensor signal, the membership function of fault model can be modeled as:
μ F x : X 0 , 1 , x X ,
where μ F represents the membership function of fault F under the corresponding attribute.
Two main factors are considered in the determination of membership function: the working performance of the sensor itself such as the linearity of sensor measurements and various kinds of interferences in the working process of the sensor such as the interference of mechanical noise and electromagnetic waves. If only the second factor is considered, the probability density function of the measured value of the same physical quantity is generally regarded as a form of Gaussian distribution. Gaussian distribution possesses the following advantages [72]: first, the Gaussian distribution is manageable, that is, a large number of data that is related to Gaussian distribution can be put through it to obtain a exact form. Second, if the error can be seen as superposition of many independent random variables. According to the central limit theorem, the error is supposed to have the form of Gaussian distribution. Third, Gaussian distribution has an extremely wide range of practical background, and a lot of the probability distribution of random variables in production and scientific experiments can be approximately described by Gaussian distribution. We take the Gaussian type of membership function as a fault model based on the above several points, and the details of membership function determination are presented as follows:
  • Assuming n groups of data are observed, each group consists of m consecutive observations in the same time interval Δ t . In order to make the generated fault model representative, n should be more than 3 and 30, respectively, and n groups of data should be observed at different times.
  • Calculate the mean value μ and standard deviation σ of the n groups’ observations. Here, observation i in group j is denoted as x i , j . The mean value μ and standard deviation σ can be obtained as:
    X ¯ = i = 1 m j = 1 n x i , j n m ,
    σ = i = 1 m j = 1 n ( x i , j X ¯ ) 2 n m 1 .
  • The Gaussian type of membership function is obtained as:
    μ F x = exp ( x X ¯ ) 2 2 σ 2 .

3.2. Test Model Modeling

In fault diagnosis, considering the influences of the interferences in the working environment of the sensor, it is difficult to judge the fault class precisely from a single measurement. When multiple measurements on the equipment in the same time interval Δ t are obtained, the Gaussian type membership function can be obtained using Equations (15)–(17), which is denoted as:
μ T x : X 0 , 1 , x X .
μ T is taken as a test model. Compared with the single measurement and mean value of multiple measurements, it is more comprehensive and objective to reflect the change of the value of equipment fault characteristic parameters.

3.3. BPA Determination

The determination of BPA is a key step in the application of D–S evidence theory. In the application of data fusion based on D–S evidence theory, the BPAs generated from different attributes will be combined to make a decision. While how to determine BPA from multi-sensor data fusion system is still an open issue, much research has been carried out [53,73]. The BPA determination method in [53] utilized the relationship between test data and Gaussian distribution model generated by training data. However, the method of determining BPA with single measurement values can not cover the fuzziness of sensors in complex working environments. Taking this fact into consideration, a BPA determination method based on multiple measurements value is proposed in this section.
Assuming the equipment in Figure 1 may have M faults occur, a total of N (N >L) attributes are extracted from L sensors. The frame of discernment should be Θ = F 1 , F 2 , F i , , F M . For a certain attribute k (k = 1, 2, ...,N), the test model under attribute k is put into its corresponding M fault models under attribute k to determine B P A k .
Let μ F i ( x ) , X ¯ i , σ i be the function curve, mean value, standard deviation of generated model of fault F i (i = 1, 2, 3, ..., M) under attribute k, respectively. Let μ T ( x ) , X ¯ T , σ T be the function curve, mean value and standard deviation of generated model using observation under attribute k, respectively. The B P A k calculation procedures are as the following: the proposition support degree calculation method is given in Section 3.3.1, Section 3.3.2 and Section 3.3.3. B P A k is obtained by normalizing all obtained support degrees to propositions, which is shown in Section 3.3.4.

3.3.1. Support Degree Calculation of Singleton Propositions

In this part, the support degree of the test model to proposition { F i } denoted as S u p ( { F i } ) (i = 1, 2, ..., M) is obtained.
Taking into account the property of the model curve, firstly let us calculate the integral interval:
I n M i n = min ( X ¯ i + 3 σ i , X ¯ i 3 σ i , X ¯ T + 3 σ T , X ¯ T 3 σ T ) ,
I n M a x = max ( X ¯ i + 3 σ i , X ¯ i 3 σ i , X ¯ T + 3 σ T , X ¯ T 3 σ T ) .
Then, the intersection area between the test model curve μ T ( x ) and the fault model curve μ F i ( x ) is calculated:
A r e a ( { F i } ) = I n M i n I n M a x min ( μ F i ( x ) , μ T ( x ) ) d x .
For example, A r e a ( { F 1 } ) is as shown in Figure 2a if M = 3 .
Finally, calculate the support degree S u p ( { F i } ) of the test model to proposition { F i } :
S u p ( { F i } ) = A r e a ( { F i } ) A r e a T ,
where the A r e a T is the area under the test model curve μ T x , A r e a T = 2 π σ T .

3.3.2. Support Degree Calculation of Propositions with Two Elements

In this part, the support degree of the test model to compound proposition { F i , F j } denoted as S u p ( { F i , F j } ) ( i , j = 1, 2, ..., M , i < j) is obtained.
Considering the nature of the model curve, firstly let us calculate the integral interval:
I n M i n = min ( X ¯ i + 3 σ i , X ¯ i 3 σ i , X ¯ j + 3 σ j , X ¯ j 3 σ j , X ¯ T + 3 σ T , X ¯ T 3 σ T ) ,
I n M a x = max ( X ¯ i + 3 σ i , X ¯ i 3 σ i , X ¯ j + 3 σ j , X ¯ j 3 σ j , X ¯ T + 3 σ T , X ¯ T 3 σ T ) .
Then, the intersection area between the test model curve μ T ( x ) and the fault model curves μ F i ( x ) , μ F j ( x ) is calculated:
A r e a ( { F i , F j } ) = I n M i n I n M a x min ( μ F i ( x ) , μ F j ( x ) , μ T ( x ) ) d x .
For example, A r e a ( { F 1 , F 2 } ) is as shown in Figure 2b if M = 3 .
Finally, the support degree S u p ( { F i , F j } ) of the test model to compound proposition { F i , F j } is obtained:
S u p ( { F i , F j } ) = A r e a ( { F i , F j } ) A r e a T ,
where the A r e a T is the area under the test model curve μ T x , A r e a T = 2 π σ T .

3.3.3. Support Degree Calculation of Propositions with Multiple Elements

In this part, the support degree of the test model to compound proposition { F n 1 , F n 2 , . . . , F n i } denoted as S u p ( { F n 1 , F n 2 , . . . , F n i } ) ( n 1 = 1 , 2 , . . . , M , n 2 = 1 , 2 , . . . , M , . . . , n i = 1 , 2 , . . . , M , n 1 < n 2 < . . . < n i ) is obtained.
Taking into account the nature of the model curve, firstly let us calculate the integral interval:
I n M i n = min ( min n j ( X ¯ n j + 3 σ n j , X ¯ n j 3 σ n j ) , X ¯ T + 3 σ T , X ¯ T 3 σ T ) ( n j = n 1 , n 2 , . . . , n i ) ,
I n M a x = max ( max n j ( X ¯ n j + 3 σ n j , X ¯ n j 3 σ n j ) , X ¯ T + 3 σ T , X ¯ T 3 σ T ) ( n j = n 1 , n 2 , . . . , n i ) .
Then, the intersection area between the test model curve μ T ( x ) and the fault model curves μ F n j ( x ) is obtained as:
A r e a ( { F n 1 , F n 2 , . . . , F n i } ) = I n M i n I n M a x min ( min n j ( μ F n j ( x ) ) , μ T ( x ) ) d x ( n j = n 1 , n 2 , . . . , n i ) .
For example, A r e a ( { F 1 , F 2 , F 3 } ) is as shown in Figure 2c if M = 3 .
Finally, the support degree S u p ( { F n 1 , F n 2 , . . . , F n i } ) of the test model to compound proposition { F n 1 , F n 2 , . . . , F n i } is obtained:
S u p ( { F n 1 , F n 2 , . . . , F n i } ) = A r e a ( { F n 1 , F n 2 , . . . , F n i } ) A r e a T ,
where the A r e a T is the area under the test model curve μ T x , A r e a T = 2 π σ T .

3.3.4. Normalization

The sum of all obtained support degrees is denoted as S u p S .
S u p S = S u p ( { F 1 } ) + S u p ( { F 2 } ) + . . . + S u p ( { F M } ) + S u p ( { F 1 , F 2 } ) + . . . + S u p ( { F 1 , . . . , F i } ) + . . . + S u p ( { F 1 , . . . , F M } ) .
The obtained support degrees are normalized to obtain B P A k under attribute k:
m F i = S u p F i S u p S ( i = 1 , 2 , . . . , M ) m F i , F j = S u p F i , F j S u p S ( i , j = 1 , 2 , . . . , M , i < j ) . . . m F 1 , F 2 , . . . , F M = S u p F 1 , F 2 , . . . , F M S u p S .

3.4. BPA Combination and Decision Making

After B P A 1 , B P A 2 , . . . , B P A N are obtained through the method in Section 3.3, these N BPAs can be combined using Equations (8)–(12) to obtain final BPA. Decision making rules based on BPA are required because it is a set of BPA ultimately obtained through the fault diagnosis method based on multi-sensor data fusion. Hence, the decision making rules based on the obtained final BPA in this paper are as follows. If F k is the final diagnostic fault class of the equipment, the following decision making rules should be simultaneously met:
  • F k should have the maximum belief in final BPA that is m ( { F k } ) = max i ( m ( { F i } ) ) (i = 1, 2, ..., M) and m ( { F k } ) should exceed a threshold k 1 . k 1 is the minimum belief of F k . k 1 should be a slightly larger value to ensure the high credibility of diagnostic result F k so as to make sure the correctness of diagnostic result. Here, k 1 is designated as 0.5. The purpose of this rule is to make the reliability of the diagnostic result F k sufficiently large.
  • The BPA of dual set proposition and multi-subset proposition, such as m ( { F 1 , F 2 } ) , m ( { F 1 , F 3 , F 6 } ) , m ( { Θ } ) , should be less than a certain threshold k 2 . k 2 is the maximum belief of the uncertain components. If larger than k 2 , it is considered that the uncertainty of diagnostic result is too large, which may lead to the wrong diagnostic result. k 2 should be a smaller value so that the uncertainty of the diagnostic result is small. Here, k 2 is designated as 0.3. The purpose of this rule is to make the uncertainty of the diagnostic result smaller.
  • The difference between m ( { F k } ) and m ( { F l } ) should be larger than a certain threshold k 3 that is m ( { F k } ) m ( { F l } ) > k 3 , where m ( { F l } ) = max i ( m ( { F i } ) ) (i = 1, 2, ..., M , i k ). k 3 is minimum difference between the maximum and the second largest in m ( { F i } ) ) (i = 1, 2, ..., M). If larger than k 3 , it is considered that the diagnostic result can be clearly distinguished from the credibility of the other fault types, which can avoid false diagnostic results. Hence, to avoid false diagnostic results, k 3 can not be too small. Here, k 3 is designated as 0.2 because k 3 that is too large may cause a missing alarm. The purpose of this rule is to make the difference between the maximum and the second largest in m ( { F i } ) ) larger, which means that m ( { F k } ) is relatively larger in m ( { F i } ) ) (i = 1, 2, ..., M).
If such F k can not be found, it is considered that the decision making can not be made by the obtained final BPA. The above decision making rules are based on the same fundamental principle with traditional decision making methods based on BPA [74].

4. Applications

In this section, two examples are given to demonstrate the effectiveness of the proposed method. In the first example of Iris data set classification, which is similar to fault diagnosis, multiple instances are selected as the test set to simulate the selection of the test sample in fault diagnosis, and the identification accuracy is high. In the second application, the fault diagnosis experiment of a motor rotor is carried out to demonstrate the validity in engine fault diagnosis.

4.1. Iris Data Set Classification

The Iris data set [75], which is perhaps the best known database to be found in pattern recognition literature, involves classification of three species of the Iris flowers, Iris setosa (S), Iris versicolour (E), and Iris virginica (V) on the basis of four numeric attributes of the Iris flower: sepal length ( S L ) in cm, sepal width ( S W ) in cm, petal length ( P L ) in cm, and petal width ( P W ) in cm. The frame of discernment should be Θ = { S , E , V } .
In the Iris data set, there are 50 instances for each of three species. The data are obtained from the UCI (University of California Irvine) Machine Learning Repository [76] of machine learning databases. Among 50 instances of each species, 20 instances are randomly selected as the training set; thus, there are a total of 12 training models (Table 1 and Figure 3). Ten instances are randomly selected in the rest of 30 instances under each attribute of E (Iris versicolour) as the test set; thus, there are a total of four test models (Table 1 and Figure 3). Each of the four attributes is regarded as an information source. The final identification result should be E (Iris versicolour).
First, calculate the support degrees, and the results are shown in Table 2. Then, the BPA of each attribute is obtained through the normalization of support degrees (Table 3). After that, using the combination method mentioned in Section 2.3, the final BPA is obtained (Table 3). Finally, through the decision making rules mentioned in Section 3.4, the final identification result is E (Iris versicolour).
Samples from S, E, and V were tested as test sets 10,000 times, respectively, and the accuracy of the method on Iris data set classification is presented in Table 4 using the decision method mentioned in Section 3.4. The overall recognition rate of the three categories is 99.98%.

4.2. Motor Rotor Fault Diagnosis

In this section, the proposed method is applied in practical engine fault diagnosis, and experimental data is from actual measurement data of the motor rotor [77]. Feature j under fault i is denoted as F i j , where i = 1 , 2 , 3 , respectively, on behalf of the three kinds of faults: rotor imbalance, rotor misalignment and the base loose. In addition, j = 1 , 2 , 3 , 4 , respectively, on behalf of the four features: the amplitude of fundamental harmonic ( 1 X ), the amplitude of second harmonic ( 2 X ), and the amplitude of the third harmonic ( 3 X ) in the frequency spectrum of the rotor vibration acceleration sensor signal and shifting of rotor vibration ( S h i f t i n g ). For example, F 11 means feature 1 X of fault rotor imbalance. The frame of discernment should be Θ = { F 1 , F 2 , F 3 } . There are five sets of measured data for specific i and j and each set of measured data is composed of forty continuous measurements in sixteen seconds.
Step 1: Establish fault models μ F i j . For specific i and j, four sets of feature data are randomly selected as fault samples to generate fault model μ F i j ; thus, a total of 12 fault models can be established for all i and j. For example, ( X ¯ μ F i 4 , σ μ F i 4 ) is shown in Table 5 and the function curve of μ F i 4 is shown in Figure 4a.
Step 2: Generate test models μ T j . Assuming the left four sets of feature data in F 2 j are selected as the test sample, which can be used to generate test models μ T j , the final diagnostic result should be F 2 . For example, ( X ¯ μ T 4 , σ μ T 4 ) is shown in Table 5 and the function curve of μ T 4 is shown in Figure 4b.
Step 3: Generate B P A j . After 12 fault models μ F i j and four test models μ T j are obtained, the B P A j can be generated by matching the fault models μ F i j with the test model μ T j . Taking test model F 24 as an example, a detailed introduction to the BPA generation process of B P A 4 is given as follows with graphical assistant instructions:
In Figure 5, first, the integral interval is calculated:
  • For A r e a ( { F 1 } ) , the integral interval is (3.0740,6.9470).
  • For A r e a ( { F 2 } ) , the integral interval is (3.0740,6.9470).
  • For A r e a ( { F 1 , F 2 } ) , the integral interval is (3.0740,6.9470).
Then, the intersection area between fault models and test model, which is represented by the shadow in Figure 5, is obtained as:
A r e a ( { F 1 } ) = 0 . 5686 , A r e a ( { F 1 } ) = 0 . 7205 , A r e a ( { F 1 , F 2 } ) = 0 . 4666 .
Third, the calculated support degrees are as follows:
S u p F 1 = 0 . 3514 , S u p F 2 = 0 . 4453 , S u p F 1 , F 2 = 0 . 2884 .
Finally, normalizing the obtained support degrees, B P A 4 (BPA under feature 4 ( S h i f t i n g )) is obtained as:
m F 1 = 0 . 3238 , m F 2 = 0 . 4104 , m F 1 , F 2 = 0 . 2658 .
Use the same method to determine B P A i under feature i (i =1, 2, 3). The results are shown in Table 6.
Step 4: Combine the obtained BPAs and give the result of diagnosis. Using the combination method mentioned in Section 3.4, the final BPA is obtained (Table 6). Through the decision making rules mentioned in Section 3.4, the final diagnostic result is fault F 2 (rotor misalignment).
Through five experiments, the experiment results of the proposed fault diagnosis method is presented in Table 7 using decision making rules mentioned in Section 3.4. The mean value and standard deviation of models in each case are as follows:
For case 1: ( X ¯ μ F 11 , σ μ F 11 ) = ( 0 . 1614 , 0 . 0108 ) , ( X ¯ μ F 12 , σ μ F 12 ) = ( 0 . 1508 , 0 . 0171 ) , ( X ¯ μ F 13 , σ μ F 13 ) = ( 0 . 1131 , 0 . 0071 ) , ( X ¯ μ F 14 , σ μ F 14 ) = ( 4 . 3241 , 0 . 3240 ) , ( X ¯ μ F 21 , σ μ F 21 ) = ( 0 . 1794 , 0 . 0125 ) , ( X ¯ μ F 22 , σ μ F 22 ) = ( 0 . 3300 , 0 . 0110 ) , ( X ¯ μ F 23 , σ μ F 23 ) = ( 0 . 2332 , 0 . 0299 ) , ( X ¯ μ F 24 , σ μ F 24 ) = ( 4 . 6414 , 0 . 3087 ) , ( X ¯ μ F 31 , σ μ F 31 ) = ( 0 . 3301 , 0 . 0075 ) , ( X ¯ μ F 32 , σ μ F 32 ) = ( 0 . 3422 , 0 . 0150 ) , ( X ¯ μ F 33 , σ μ F 33 ) = ( 0 . 1367 , 0 . 0119 ) , ( X ¯ μ F 34 , σ μ F 34 ) = ( 9 . 8220 , 0 . 1010 ) , ( X ¯ μ T 1 , σ μ T 1 ) = ( 0 . 1616 , 0 . 0006 ) , ( X ¯ μ T 2 , σ μ T 2 ) = ( 0 . 1428 , 0 . 0025 ) , ( X ¯ μ T 3 , σ μ T 3 ) = ( 0 . 1095 , 0 . 0021 ) , ( X ¯ μ T 4 , σ μ T 4 ) = ( 4 . 3319 , 0 . 0347 ) .
For case 2: ( X ¯ μ F 11 , σ μ F 11 ) = ( 0 . 1619 , 0 . 0101 ) , ( X ¯ μ F 12 , σ μ F 12 ) = ( 0 . 1481 , 0 . 0164 ) , ( X ¯ μ F 13 , σ μ F 13 ) = ( 0 . 1131 , 0 . 0071 ) , ( X ¯ μ F 14 , σ μ F 14 ) = ( 4 . 3035 , 0 . 3204 ) , ( X ¯ μ F 21 , σ μ F 21 ) = ( 0 . 1848 , 0 . 0109 ) , ( X ¯ μ F 22 , σ μ F 22 ) = ( 0 . 3320 , 0 . 0106 ) , ( X ¯ μ F 23 , σ μ F 23 ) = ( 0 . 2492 , 0 . 0342 ) , ( X ¯ μ F 24 , σ μ F 24 ) = ( 4 . 7851 , 0 . 4010 ) , ( X ¯ μ F 31 , σ μ F 31 ) = ( 0 . 3271 , 0 . 0050 ) , ( X ¯ μ F 32 , σ μ F 32 ) = ( 0 . 3446 , 0 . 0141 ) , ( X ¯ μ F 33 , σ μ F 33 ) = ( 0 . 1367 , 0 . 0113 ) , ( X ¯ μ F 34 , σ μ F 34 ) = ( 9 . 8182 , 0 . 1090 ) , ( X ¯ μ T 1 , σ μ T 1 ) = ( 0 . 1696 , 0 . 0096 ) , ( X ¯ μ T 2 , σ μ T 2 ) = ( 0 . 3187 , 0 . 0041 ) , ( X ¯ μ T 3 , σ μ T 3 ) = ( 0 . 2131 , 0 . 0053 ) , ( X ¯ μ T 4 , σ μ T 4 ) = ( 4 . 4358 , 0 . 4015 ) .
For case 3: ( X ¯ μ F 11 , σ μ F 11 ) = ( 0 . 1619 , 0 . 0107 ) , ( X ¯ μ F 12 , σ μ F 12 ) = ( 0 . 1481 , 0 . 0164 ) , ( X ¯ μ F 13 , σ μ F 13 ) = ( 0 . 1114 , 0 . 0050 ) , ( X ¯ μ F 14 , σ μ F 14 ) = ( 4 . 3241 , 0 . 3240 ) , ( X ¯ μ F 21 , σ μ F 21 ) = ( 0 . 1790 , 0 . 0101 ) , ( X ¯ μ F 22 , σ μ F 22 ) = ( 0 . 3320 , 0 . 0106 ) , ( X ¯ μ F 23 , σ μ F 23 ) = ( 0 . 2492 , 0 . 0342 ) , ( X ¯ μ F 24 , σ μ F 24 ) = ( 4 . 6484 , 0 . 4402 ) , ( X ¯ μ F 31 , σ μ F 31 ) = ( 0 . 3299 , 0 . 0079 ) , ( X ¯ μ F 32 , σ μ F 32 ) = ( 0 . 3446 , 0 . 0141 ) , ( X ¯ μ F 33 , σ μ F 33 ) = ( 0 . 1367 , 0 . 0113 ) , ( X ¯ μ F 34 , σ μ F 34 ) = ( 9 . 8220 , 0 . 1010 ) , ( X ¯ μ T 1 , σ μ T 1 ) = ( 0 . 3275 , 0 . 0023 ) , ( X ¯ μ T 2 , σ μ T 2 ) = ( 0 . 3409 , 0 . 0135 ) , ( X ¯ μ T 3 , σ μ T 3 ) = ( 0 . 1341 , 0 . 0113 ) , ( X ¯ μ T 4 , σ μ T 4 ) = ( 9 . 7652 , 0 . 0953 ) .
For case 4: ( X ¯ μ F 11 , σ μ F 11 ) = ( 0 . 1619 , 0 . 0101 ) , ( X ¯ μ F 12 , σ μ F 12 ) = ( 0 . 1488 , 0 . 0173 ) , ( X ¯ μ F 13 , σ μ F 13 ) = ( 0 . 1139 , 0 . 0062 ) , ( X ¯ μ F 14 , σ μ F 14 ) = ( 4 . 3414 , 0 . 0745 ) , ( X ¯ μ F 21 , σ μ F 21 ) = ( 0 . 1790 , 0 . 0101 ) , ( X ¯ μ F 22 , σ μ F 22 ) = ( 0 . 3279 , 0 . 0115 ) , ( X ¯ μ F 23 , σ μ F 23 ) = ( 0 . 2492 , 0 . 0342 ) , ( X ¯ μ F 24 , σ μ F 24 ) = ( 4 . 6414 , 0 . 3087 ) , ( X ¯ μ F 31 , σ μ F 31 ) = ( 0 . 3271 , 0 . 0050 ) , ( X ¯ μ F 32 , σ μ F 32 ) = ( 0 . 3422 , 0 . 0150 ) , ( X ¯ μ F 33 , σ μ F 33 ) = ( 0 . 1367 , 0 . 0119 ) , ( X ¯ μ F 34 , σ μ F 34 ) = ( 9 . 7948 , 0 . 0888 ) , ( X ¯ μ T 1 , σ μ T 1 ) = ( 0 . 3387 , 0 . 0071 ) , ( X ¯ μ T 2 , σ μ T 2 ) = ( 0 . 3503 , 0 . 0060 ) , ( X ¯ μ T 3 , σ μ T 3 ) = ( 0 . 1341 , 0 . 0080 ) , ( X ¯ μ T 4 , σ μ T 4 ) = ( 9 . 8739 , 0 . 1267 ) .
For case 5: ( X ¯ μ F 11 , σ μ F 11 ) = ( 0 . 1607 , 0 . 0106 ) , ( X ¯ μ F 12 , σ μ F 12 ) = ( 0 . 1486 , 0 . 0076 ) , ( X ¯ μ F 13 , σ μ F 13 ) = ( 0 . 1114 , 0 . 0050 ) , ( X ¯ μ F 14 , σ μ F 14 ) = ( 4 . 3286 , 0 . 3241 ) , ( X ¯ μ F 21 , σ μ F 21 ) = ( 0 . 1790 , 0 . 0101 ) , ( X ¯ μ F 22 , σ μ F 22 ) = ( 0 . 3300 , 0 . 0110 ) , ( X ¯ μ F 23 , σ μ F 23 ) = ( 0 . 2461 , 0 . 0362 ) , ( X ¯ μ F 24 , σ μ F 24 ) = ( 4 . 7533 , 0 . 4653 ) , ( X ¯ μ F 31 , σ μ F 31 ) = ( 0 . 3294 , 0 . 0079 ) , ( X ¯ μ F 32 , σ μ F 32 ) = ( 0 . 3445 , 0 . 0118 ) , ( X ¯ μ F 33 , σ μ F 33 ) = ( 0 . 1370 , 0 . 0122 ) , ( X ¯ μ F 34 , σ μ F 34 ) = ( 9 . 8169 , 0 . 1019 ) , ( X ¯ μ T 1 , σ μ T 1 ) = ( 0 . 1932 , 0 . 0138 ) , ( X ¯ μ T 2 , σ μ T 2 ) = ( 0 . 3266 , 0 . 0108 ) , ( X ¯ μ T 3 , σ μ T 3 ) = ( 0 . 2255 , 0 . 0135 ) , ( X ¯ μ T 4 , σ μ T 4 ) = ( 4 . 5631 , 0 . 0678 ) .

5. Conclusions

This paper presented a new fault diagnosis method based on multi-sensor data fusion. The Gaussian type of membership function is used to construct a fault model and a test model due to its advantages in processing sensor data. First, the distribution of the measurements of the sensor, which is working in a complex environment, is generally regarded as Gaussian distribution if the influence of the performance of the sensor is ignored. Second, the Gaussian type of membership function is tractable in dealing with a large number of data involving Gaussian distribution. The overall diagnostic processes are as follows: first, obtain the Gaussian type of fault models through a certain amount of fault samples. Secondly, through the measurements of sensors in equipment, Gaussian types of test models are determined. Thirdly, the intersection area between test models and fault models with identical attributes is transformed into a set of BPA. Finally, a weighted average approach is used to combine the obtained BPAs, and final diagnostic result is obtained through the proposed multiple decision making rules. On the whole, the proposed fault diagnosis method has the following advantages: the Gaussian type of fault model and test model determination method can fit the measurement data well due to the complex working environment of sensors in fault diagnosis. The basic probability assignment determination method based on the intersection area between the test model and fault models can better express uncertain information because this method takes into account the variance information of the models. Finally, the effectiveness of the proposed method is verified by some illustrative examples.

Acknowledgments

The work is partially supported by the National Natural Science Foundation of China (Grant No. 61671384), the Natural Science Basic Research Plan in Shaanxi Province of China (Program No. 2016JM6018), the Aviation Science Foundation (Program No. 20165553036), the Shanghai Aerospace Science and Technology Innovation Fund (Program No. SAST2016083), and the Seed Foundation of Innovation and Creation for Graduate Students in Northwestern Polytechnical University (Program No. Z2016122).

Author Contributions

Wen Jiang and Weiwei Hu conceived and designed the experiments; Wen Jiang performed the experiments; Chunhe Xie and Weiwei Hu analyzed the data; Chunhe Xie contributed analysis tools; and Weiwei Hu wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, J.; Patton, R.J. Robust Model-Based Fault Diagnosis for Dynamic Systems; Springer: New York, NY, USA, 1999. [Google Scholar]
  2. Isermann, R. Fault-Diagnosis Systems. An Introduction from Fault Detection to Fault Tolerance; Springer: New York, NY, USA, 2006; Volume 28, pp. 195–197. [Google Scholar]
  3. Al-Busaidi, W.; Pilidis, P. Modelling of the non-reactive deposits impact on centrifugal compressor aerothermo dynamic performance. Eng. Fail. Anal. 2016, 60, 57–85. [Google Scholar] [CrossRef]
  4. Xue, R.; Hu, C.; Sethi, V.; Nikolaidis, T.; Pilidis, P. Effect of steam addition on gas turbine combustor design and performance. Appl. Therm. Eng. 2016, 104, 249–257. [Google Scholar] [CrossRef]
  5. Ramesh, T.S.; Shum, S.K.; Davis, J.F. A structured framework for efficient problem solving in diagnostic expert systems. Comput. Chem. Eng. 1988, 12, 891–902. [Google Scholar] [CrossRef]
  6. Wu, J.D.; Wang, Y.H.; Bai, M.R. Development of an expert system for fault diagnosis in scooter engine platform using fuzzy-logic inference. Expert Syst. Appl. 2007, 33, 1063–1075. [Google Scholar] [CrossRef]
  7. Yang, J.B.; Liu, J.; Wang, J.; Sii, H.S.; Wang, H.W. Belief rule-base inference methodology using the evidential reasoning Approach-RIMER. IEEE Trans. Syst. Man Cybern. A Syst. Hum. 2006, 36, 266–285. [Google Scholar] [CrossRef]
  8. Tayarani-Bathaie, S.S.; Vanini, Z.S.; Khorasani, K. Dynamic neural network-based fault diagnosis of gas turbine engines. Neurocomputing 2014, 125, 153–165. [Google Scholar] [CrossRef]
  9. Tamilselvan, P.; Wang, P. Failure diagnosis using deep belief learning based health state classification. Reliab. Eng. Syst. Saf. 2013, 115, 124–135. [Google Scholar] [CrossRef]
  10. Yan, R.; Gao, R.X.; Chen, X. Wavelets for fault diagnosis of rotary machines: A review with applications. Signal Process. 2014, 96, 1–15. [Google Scholar] [CrossRef]
  11. Pichler, K.; Lughofer, E.; Pichler, M.; Buchegger, T.; Klement, E.P.; Huschenbett, M. Detecting cracks in reciprocating compressor valves using pattern recognition in the pV diagram. Pattern Anal. Appl. 2015, 18, 461–472. [Google Scholar] [CrossRef]
  12. Pichler, K.; Lughofer, E.; Pichler, M.; Buchegger, T.; Klement, E.P.; Huschenbett, M. Fault detection in reciprocating compressor valves under varying load conditions. Mech. Syst. Signal Process. 2016, 70–71, 104–119. [Google Scholar] [CrossRef]
  13. Fernando, H.; Surgenor, B. An unsupervised artificial neural network versus a rule-based approach for fault detection and identification in an automated assembly machine. Robot. Comput. Integr. Manuf. 2015, 43, 79–88. [Google Scholar] [CrossRef]
  14. Delgado-Arredondo, P.A.; Morinigo-Sotelo, D.; Osornio-Rios, R.A.; Avina-Cervantes, J.G.; Rostro-Gonzalez, H.; Romero-Troncoso, R.D.J. Methodology for fault detection in induction motors via sound and vibration signals. Mech. Syst. Signal Process. 2016, 83, 568–589. [Google Scholar] [CrossRef]
  15. Serdio, F.; Lughofer, E.; Pichler, K.; Buchegger, T.; Pichler, M.; Efendic, H. Fault detection in multi-sensor networks based on multivariate time-series models and orthogonal transformations. Inf. Fusion 2014, 20, 272–291. [Google Scholar] [CrossRef]
  16. Serdio, F.; Lughofer, E.; Zavoianu, A.C.; Pichler, K.; Pichler, M.; Buchegger, T.; Efendic, H. Improved Fault Detection employing Hybrid Memetic Fuzzy Modeling and Adaptive Filters. Appl. Soft Comput. 2017, 51, 60–82. [Google Scholar] [CrossRef]
  17. Serdio, F.; Lughofer, E.; Pichler, K.; Pichler, M.; Buchegger, T.; Efendic, H. Fuzzy Fault Isolation using Gradient Information and Quality Criteria from System Identification Models. Inf. Sci. 2015, 316, 18–39. [Google Scholar] [CrossRef]
  18. Mitchell, H.B. Multi-Sensor Data Fusion; Artech House: Norwood, MA, USA, 1990; pp. 585–610. [Google Scholar]
  19. Crowley, J.L. Principles and Techniques for Sensor Data Fusion; Springer: Berlin/Heidelberg, Germany, 1993; pp. 5–27. [Google Scholar]
  20. Alexandridis, A. Evolving RBF neural networks for adaptive soft-sensor design. Int. J. Neural Syst. 2013, 23, 1350029. [Google Scholar] [CrossRef] [PubMed]
  21. Graziani, S.; Pagano, F.; Xibilia, M.G. Soft sensor for a Propylene Splitter with seasonal variations. In Proceedings of the 2010 IEEE Instrumentation and Measurement Technology Conference (I2MTC), Austin, TX, USA, 3–6 May 2010; pp. 273–278.
  22. Dempster, A.P. Upper and lower probabilities induced by a multivalued mapping. Ann. Math. Stat. 1967, 38, 325–339. [Google Scholar] [CrossRef]
  23. Shafer, G. A mathematical theory of evidence. Technometrics 1976, 20, 242. [Google Scholar]
  24. Deng, X.; Xiao, F.; Deng, Y. An improved distance-based total uncertainty measure in belief function theory. Appl. Intell. 2017, in press. [Google Scholar] [CrossRef]
  25. Deng, X.; Han, D.; Dezert, J.; Deng, Y.; Shyr, Y. Evidence combination from an evolutionary game theory perspective. IEEE Trans. Cybern. 2016, 46, 2070–2082. [Google Scholar] [CrossRef] [PubMed]
  26. Jiang, W.; Zhan, J. A modified combination rule in generalized evidence theory. Appl. Intell. 2017. [Google Scholar] [CrossRef]
  27. Zadeh, L.A. Fuzzy sets. Inf. Control 1965, 8, 338–353. [Google Scholar] [CrossRef]
  28. Walczak, B.; Massart, D.L. Rough sets theory. Chemom. Intell. Lab. Syst. 1999, 47, 1–16. [Google Scholar] [CrossRef]
  29. Shen, L.; Tay, F.E.H.; Qu, L.; Shen, Y. Fault diagnosis using Rough Sets Theory. Comput. Ind. 2000, 43, 61–72. [Google Scholar] [CrossRef]
  30. Zadeh, L.A. A Note on Z-numbers. Inf. Sci. 2011, 181, 2923–2932. [Google Scholar] [CrossRef]
  31. Jiang, W.; Xie, C.; Zhuang, M.; Shou, Y.; Tang, Y. Sensor Data Fusion with Z-Numbers and Its Application in Fault Diagnosis. Sensors 2016, 16, 1509. [Google Scholar] [CrossRef] [PubMed]
  32. Zhou, X.; Shi, Y.; Deng, X.; Deng, Y. D-DEMATEL: A new method to identify critical success factors in emergency management. Saf. Sci. 2017, 91, 93–104. [Google Scholar] [CrossRef]
  33. Zhou, X.; Deng, X.; Deng, Y.; Mahadevan, S. Dependence assessment in human reliability analysis based on D numbers and AHP. Nucl. Eng. Des. 2017, 313, 243–252. [Google Scholar] [CrossRef]
  34. Deng, X.; Hu, Y.; Deng, Y.; Mahadevan, S. Supplier selection using AHP methodology extended by D numbers. Expert Syst. Appl. 2014, 41, 156–167. [Google Scholar] [CrossRef]
  35. Deng, X.; Lu, X.; Chan, F.T.; Sadiq, R.; Mahadevan, S.; Deng, Y. D-CFPR: D numbers extended consistent fuzzy preference relations. Knowl. Based Syst. 2015, 73, 61–68. [Google Scholar] [CrossRef]
  36. Yao, J.; Wu, C.; Xie, X.; Qian, K. A new method of information decision-making based on D–S evidence theory. In Proceedings of the IEEE International Conference on Systems Man and Cybernetics, Istanbul, Turkey, 10–13 October 2010; pp. 1804–1811.
  37. Li, Z.; Wen, G.; Xie, N. An approach to fuzzy soft sets in decision making based on grey relational analysis and Dempster–Shafer theory of evidence. Artif. Intell. Med. 2015, 64, 161–171. [Google Scholar] [CrossRef] [PubMed]
  38. Xie, N.; Wen, G.; Li, Z. A method for fuzzy soft sets in decision making based on grey relational analysis and d-s theory of evidence: Application to medical diagnosis. Comput. Math. Methods Med. 2014, 2014, 581316. [Google Scholar] [CrossRef] [PubMed]
  39. Li, L.L.; Li, Z.G.; Zhao, Q.M.; Zhang, H.J. An Algorithm of General Closeness Degree Based on Evidence Theory and Its Applications in Fuzzy Pattern Recognition. In Proceedings of the 9th International Conference on Control, Automation, Robotics and Vision, Singapore, 5–8 December 2006; pp. 1–5.
  40. Jiang, H.Q.; Gao, J.M.; Liang, Z.M.; Gao, Z.Y.; Zhang, X.W. Pattern recognition of die casting process based on D–S evidence theory. Jisuanji Jicheng Zhizao Xitong/Comput. Integr. Manuf. Syst. Cims 2015, 21, 1343–1349. [Google Scholar]
  41. Denœux, T. Application of evidence theory to k -NN pattern classification. Mach. Intell. Pattern Recognit. 1994, 16, 13–24. [Google Scholar]
  42. Mu, C.P.; Li, X.J.; Huang, H.K.; Tian, S.F. Online Risk Assessment of Intrusion Scenarios Using D–S Evidence Theory. In Proceedings of the 13th European Symposium on Research in Computer Security: Computer Security, Málaga, Spain, 6–8 October 2008; pp. 35–48.
  43. Lu, S.; Zhang, J.; Hao, S.; Luo, L. Security Risk Assessment Model Based on AHP/D–S Evidence Theory. In Proceedings of the 2009 International Forum on Information Technology and Applications, Chengdu, China, 15–17 May 2009; pp. 530–534.
  44. Zhang, G.; Duan, M.; Zhang, J.; Chen, D.; Yang, J. Power system risk assessment based on the evidence theory and utility theory. Autom. Electr. Power Syst. 2009, 33, 1–7. [Google Scholar]
  45. Dutta, P. Uncertainty Modeling in Risk Assessment Based on Dempster–Shafer Theory of Evidence with Generalized Fuzzy Focal Elements. Fuzzy Inf. Eng. 2015, 2, 15–30. [Google Scholar] [CrossRef]
  46. Jiang, W.; Wei, B.; Zhan, J.; Xie, C.; Zhou, D. A visibility graph power averaging aggregation operator: A methodology based on network analysis. Comput. Ind. Eng. 2016, 101, 260–268. [Google Scholar] [CrossRef]
  47. Jiang, W.; Wei, B.; Tang, Y.; Zhou, D. Ordered visibility graph average aggregation operator: An application in produced water management. Chaos 2017, 27, 023117. [Google Scholar] [CrossRef] [PubMed]
  48. Song, M.; Jiang, W.; Xie, C.; Zhou, D. A new interval numbers power average operator in multiple attribute decision making. Int. J. Intell. Syst. 2016. [Google Scholar] [CrossRef]
  49. Cheng, G.; Chen, X.H.; Shan, X.L.; Liu, H.G.; Zhou, C.F. A new method of gear fault diagnosis in strong noise based on multi-sensor information fusion. J. Vib. Control 2014, 22. [Google Scholar] [CrossRef]
  50. Li, Y.; Yang, T.; Liu, J.; Fu, N.; Wu, G. A fault diagnosis method by multi sensor fusion for spacecraft control system sensors. In Proceedings of the IEEE International Conference on Mechatronics and Automation, Harbin, China, 7–10 August 2016; pp. 748–753.
  51. Hang, J.; Zhang, J.; Cheng, M. Fault diagnosis of wind turbine based on multisensors information fusion technology. IET Renew. Power Gener. 2014, 8, 289–298. [Google Scholar] [CrossRef]
  52. Song, M.; Jiang, W. Engine fault diagnosis based on sensor data fusion using evidence theory. Adv. Mech. Eng. 2016, 8, 1–9. [Google Scholar] [CrossRef]
  53. Xu, P.; Deng, Y.; Su, X.; Mahadevan, S. A new method to determine basic probability assignment from training data. Knowl. Based Syst. 2013, 46, 69–80. [Google Scholar] [CrossRef]
  54. Hégarat-Mascle, S.L.; Bloch, I.; Vidal-Madjar, D. Application of Dempster–Shafer evidence theory to unsupervised classification in multisource remote sensing. IEEE Trans. Geosci. Remote Sens. 1997, 35, 1018–1031. [Google Scholar] [CrossRef]
  55. Salzenstein, F.; Boudraa, A.O. Unsupervised multisensor data fusion approach. In Proceedings of the Symposium on Signal Processing and ITS Applications, Sixth International, Kuala Lumpur, Malaysia, 13–16 August 2001; Volume 1, pp. 152–155.
  56. Yong, D.; WenKang, S.; ZhenFu, Z.; Qi, L. Combining belief functions based on distance of evidence. Decis. Support Syst. 2004, 38, 489–493. [Google Scholar] [CrossRef]
  57. Sirbiladze, G.; Khutsishvili, I.; Badagadze, O.; Kapanadze, M. More Precise Decision-Making Methodology in the Temporalized Body of Evidence. Application in the Information Technology Management. Int. J. Inf. Technol. Decis. Mak. 2016, 15, 1–34. [Google Scholar] [CrossRef]
  58. Bauer, M. Approximation algorithms and decision making in the Dempster–Shafer theory of evidence—An empirical study. Int. J. Approx. Reason. 1997, 17, 217–237. [Google Scholar] [CrossRef]
  59. Merigó, J.M.; Casanovas, M.; Martínez, L. Linguistic aggregation operators for linguistic decision making with Dempster–Shafer theory of evidence. Int. J. Uncertain. Fuzziness Knowl. Based Syst. 2010, 18, 287–304. [Google Scholar] [CrossRef]
  60. Baba, V.V.; Hakemzadeh, F. Toward a theory of evidence based decision making. Manag. Decis. 2012, 50, 832–867. [Google Scholar] [CrossRef]
  61. Peng, L.I.; Liu, S.F. Intuitionistic fuzzy decision-making methods based on grey incidence analysis and D–S theory of evidence. Grey Syst. 2012, 2, 54–62. [Google Scholar]
  62. Li, F.J.; Qian, Y.H.; Wang, J.T.; Liang, J.Y. Multigranulation information fusion: A dempster-shafer evidence theory based clustering ensemble method. In Proceedings of the 2015 International Conference on Machine Learning and Cybernetics, Guangzhou, China, 12–15 July 2015; pp. 58–63.
  63. Leung, Y.; Ji, N.N.; Ma, J.H. An integrated information fusion approach based on the theory of evidence and group decision-making. Inf. Fusion 2013, 14, 410–422. [Google Scholar] [CrossRef]
  64. Zhang, X.; Mahadevan, S.; Deng, X. Reliability analysis with linguistic data: An evidential network approach. Reliab. Eng. Syst. Saf. 2017, 162, 111–121. [Google Scholar] [CrossRef]
  65. Sun, R.; Huang, H.Z.; Miao, Q. Improved information fusion approach based on D–S evidence theory. J. Mech. Sci. Technol. 2008, 22, 2417–2425. [Google Scholar] [CrossRef]
  66. Jousselme, A.L.; Grenier, D.; Bossé, É. A new distance between two bodies of evidence. Inf. Fusion 2001, 2, 91–101. [Google Scholar] [CrossRef]
  67. Zadeh, L.A. A simple view of the Dempster–Shafer theory of evidence and its implication for the rule of combination. Ai Mag. 1986, 7, 85–90. [Google Scholar]
  68. Yager, R.R. On the Dempster–Shafer framework and new combination rules. Inf. Sci. Int. J. 1987, 41, 93–137. [Google Scholar] [CrossRef]
  69. Dubois, D.; Prade, H. Representation And Combination Of Uncertainty With Belief Functions And Possibility Measures. Comput. Intell. 2007, 4, 244–264. [Google Scholar] [CrossRef]
  70. Lefevre, E.; Colot, O.; Vannoorenberghe, P. Belief function combination and conflict management. Inf. Fusion 2002, 3, 149–162. [Google Scholar] [CrossRef]
  71. Murphy, C.K. Combining belief functions when evidence conflicts. Decis. Support Syst. 2000, 29, 1–9. [Google Scholar] [CrossRef]
  72. Casella, G.; Berger, R.L. Statistical Inference, 2nd ed.; Duxbury Press: Pacific Grove, CA, USA, 2002. [Google Scholar]
  73. Jiang, W.; Zhang, A.; Yang, Q. Fuzzy Approach to Construct Basic Probability Assignment and Its Application in Multi-Sensor Data Fusion Systems. Chin. J. Sens. Actuators 2008, 21, 1717–1720. [Google Scholar]
  74. Pan, Q.; Cheng, Y.; Liang, Y.; Yang, F.; Wang, X. Multi-Source Information Fusion Theory and Its Application; Tsinghua University Press: Beijing, China, 2013. [Google Scholar]
  75. Fisher, R.A. The Use of Multiple Measurements in Taxonomic Problems. Ann. Hum. Genet. 1936, 7, 179–188. [Google Scholar] [CrossRef]
  76. Lichman, M. UCI Machine Learning Repository; University of California, School of Information and Computer Science: Irvine, CA, USA, 2013; Available online: http://archive.ics.uci.edu/ml (accessed on 12 March 2017).
  77. Wen, C.; Xu, X. Theories and Applications in Multi-Source Uncertain Information Fusion—Fault Diagnosis and Reliability Evaluation; Beijing Science Press: Beijing, China, 2012. [Google Scholar]
Figure 1. The procedure of fault diagnosis with the proposed method.
Figure 1. The procedure of fault diagnosis with the proposed method.
Applsci 07 00280 g001
Figure 2. Graphical introduction of intersection area if M = 3 . (a), (b), (c) shows the A r e a { F 1 } , A r e a { F 2 , F 3 } , A r e a { F 1 , F 2 , F 3 } , respectively.
Figure 2. Graphical introduction of intersection area if M = 3 . (a), (b), (c) shows the A r e a { F 1 } , A r e a { F 2 , F 3 } , A r e a { F 1 , F 2 , F 3 } , respectively.
Applsci 07 00280 g002
Figure 3. Training models and test models of four attributes. (a), (b), (c), (d) depicts training models and test model under attribute S L , S W , P L , P W , respectively.
Figure 3. Training models and test models of four attributes. (a), (b), (c), (d) depicts training models and test model under attribute S L , S W , P L , P W , respectively.
Applsci 07 00280 g003
Figure 4. Generated fault models μ F i 4 (i = 1, 2, 3) and test model μ T 4 . (a) shows the generated fault models μ F i 4 (i = 1, 2, 3); (b) shows the generated test model μ T 4 .
Figure 4. Generated fault models μ F i 4 (i = 1, 2, 3) and test model μ T 4 . (a) shows the generated fault models μ F i 4 (i = 1, 2, 3); (b) shows the generated test model μ T 4 .
Applsci 07 00280 g004
Figure 5. Intersection area of fault models and test model under feature 4 ( S h i f t i n g ). (a), (b), (c) shows the A r e a { F 1 } , A r e a { F 2 } , A r e a { F 1 , F 2 } , respectively.
Figure 5. Intersection area of fault models and test model under feature 4 ( S h i f t i n g ). (a), (b), (c) shows the A r e a { F 1 } , A r e a { F 2 } , A r e a { F 1 , F 2 } , respectively.
Applsci 07 00280 g005
Table 1. Mean value and standard deviation of training models ( S , E , V ) and test models.
Table 1. Mean value and standard deviation of training models ( S , E , V ) and test models.
Category ( X ¯ S L , σ S L ) ( X ¯ S W , σ S W ) ( X ¯ P L , σ P L ) ( X ¯ P W , σ P W )
S(5.0050, 0.3203)(3.4000, 0.3061)(1.4450, 0.1356)(0.2450, 0.0887)
E(5.8600, 0.5510)(2.7250, 0.3567)(4.1300, 0.5497)(1.2650, 0.2277)
V(6.6250, 0.5466)(3.0000, 0.2847)(5.4400, 0.4627)(1.9350, 0.2815)
TestModel(6.1800, 0.4492)(2.7500, 0.3375)(4.4100, 0.3635)(1.3400, 0.1647)
Table 2. Support degrees.
Table 2. Support degrees.
Category S u p S S u p E S u p V S u p S , E S u p S , V S u p E , V S u p S , E , V
S L 0.10860.83240.72540.10860.05950.59340.0595
S W 0.29820.99190.63360.29740.29800.62890.2974
P L 2 × 10 9 0.91960.24161.9 × 10 9 2 × 10 11 0.23852 × 10 11
P W 1.2 × 10 5 0.96600.24701.1 × 10 5 5.5 × 10 6 0.24715.5 × 10 6
S L : sepal length; S W : sepal width; P L : petal length; P W : petal width.
Table 3. BPA (basic probability assignment) under each attribute and final BPA.
Table 3. BPA (basic probability assignment) under each attribute and final BPA.
Category m S m E m V m S , E m S , V m E , V m S , E , V
B P A S L 0.04370.33460.29160.04370.02390.23850.0239
B P A S W 0.08650.28790.18390.08630.08650.18250.0863
B P A P L 1.4 × 10 9 0.65700.17261.3 × 10 9 1.4 × 10 11 0.17041.4 × 10 11
B P A P W 8.2 × 10 6 0.66160.16928.2 × 10 6 3.8 × 10 6 0.16923.8 × 10 6
Final BPA4.9 × 10 4 0.87980.11303.3 × 10 5 2.2 × 10 6 0.00661.5 × 10 6
Table 4. The identification accuracy.
Table 4. The identification accuracy.
CategoryTest TimesIdentification Accuracy
S10,000 100 %
E10,000 99 . 95 %
V10,000 99 . 99 %
Table 5. Mean value and standard deviation of μ F i 4 (i = 1, 2, 3) and μ T 4 .
Table 5. Mean value and standard deviation of μ F i 4 (i = 1, 2, 3) and μ T 4 .
Category of ModelsMean Value ( X ¯ ) Standard Deviation ( σ )
μ F 14 4.32410.3240
μ F 24 4.64140.3087
μ F 34 9.82200.1010
μ T 4 5.01050.6455
Table 6. B P A i (i = 1, 2, 3) and final BPA.
Table 6. B P A i (i = 1, 2, 3) and final BPA.
Category m F 1 m F 2 m F 3 m F 1 , F 2 m F 1 , F 3 m F 2 , F 3 m F 1 , F 2 , F 3
B P A 1 0.07730.845200.0775000
B P A 2 00.40520.2974000.29740
B P A 3 00.99950.0002000.00030
Final BPA0.00200.99730.000560.000100.000040
Table 7. Five experiment results of motor rotor fault diagnosis.
Table 7. Five experiment results of motor rotor fault diagnosis.
CaseReal CategoryCombined BPADiagnostic Category
1 F 1 m ( F 1 ) = 0 . 9791 , m ( F 2 ) = 0 . 0199 , m ( F 1 , F 2 ) = 0 . 0010 F 1
2 F 2 m ( F 1 ) = 0 . 0681 , m ( F 2 ) = 0 . 9302 , m ( F 3 ) = 0 . 0008 , m ( F 1 , F 2 ) = 0 . 0009 F 2
3 F 3 m ( F 2 ) = 0 . 0002 , m ( F 3 ) = 0 . 9998 F 3
4 F 3 m ( F 2 ) = 0 . 0001 , m ( F 3 ) = 0 . 9999 F 3
5 F 2 m ( F 1 ) = 0 . 0124 , m ( F 2 ) = 0 . 9863 , m ( F 3 ) = 0 . 0005 , m ( F 1 , F 2 ) = 0 . 0008 F 2

Share and Cite

MDPI and ACS Style

Jiang, W.; Hu, W.; Xie, C. A New Engine Fault Diagnosis Method Based on Multi-Sensor Data Fusion. Appl. Sci. 2017, 7, 280. https://doi.org/10.3390/app7030280

AMA Style

Jiang W, Hu W, Xie C. A New Engine Fault Diagnosis Method Based on Multi-Sensor Data Fusion. Applied Sciences. 2017; 7(3):280. https://doi.org/10.3390/app7030280

Chicago/Turabian Style

Jiang, Wen, Weiwei Hu, and Chunhe Xie. 2017. "A New Engine Fault Diagnosis Method Based on Multi-Sensor Data Fusion" Applied Sciences 7, no. 3: 280. https://doi.org/10.3390/app7030280

APA Style

Jiang, W., Hu, W., & Xie, C. (2017). A New Engine Fault Diagnosis Method Based on Multi-Sensor Data Fusion. Applied Sciences, 7(3), 280. https://doi.org/10.3390/app7030280

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

Article Metrics

Back to TopTop