Fault diagnosis technology plays a vital role in the variety of critical engineering applications. Fuzzy approach is widely employed to cope with decision-making problems because it is in the simplest and most used form. This paper proposed a new similarity measure of generalized trapezoidal fuzzy numbers used for fault diagnosis. The presented similarity measure combines concepts of the geometric distance, the center of gravity point, the perimeter, and the area of the generalized trapezoidal fuzzy numbers for calculating the degree of similarity between generalized trapezoidal fuzzy numbers. This method is proposed to deal with both standardized and nonstandardized generalized trapezoidal fuzzy numbers. Some properties of the proposed similarity measure have been proved, and 12 sets of generalized fuzzy numbers have been used to compare the calculation results of the proposed similarity measures with the existing similarity measures. Comparison results indicate that the proposed similarity measure can overcome the drawbacks of existing similarity measures. Finally, a fault diagnosis experiment is carried out in laboratory based on multifunctional flexible rotor experiment bench. Experimental results demonstrate that the proposed similarity measure is more effective than other methods in terms of rotor fault diagnosis.
1. Introduction
Fault diagnosis technology has gained attention from researchers for its various applications [1–3]. To date, a large number of valuable approaches have been proposed for dealing with fault analysis issues, such as fuzzy theories [4, 5], expert system [6], wavelet analysis [7, 8], data fusion [9, 10], and neural network [11, 12]. Particularly, fuzzy approach is most successfully applied in fault diagnosis because it is in the simplest and most used form [13]. Using the fuzzy theories during the decision-making process, it is necessary to develop suitable technology to measure fuzzy similarity. Ranking the fuzzy numbers plays an important role in decision-making uncertain environment. Many methods [4, 14–20] have been presented for handing fuzzy multiple attributes decision-making problems. But they are only designed for standardized fuzzy numbers, that is,A~=(a1,a2,a3,a4;wA~) and 0≤a1≤a2≤a3≤a4≤1. However, nonstandardized fuzzy numbers are very common in reality. Wen et al. [21] proposed a new method for ranking nonstandardized fuzzy numbers, that is, A~=(a1,a2,a3,a4;wA~) and 0≤a1≤a2≤a3≤a4. The method combines the concepts of exponential distance, the perimeter, and the area of the generalized trapezoidal fuzzy numbers for calculating the degree of similarity between nonstandardized generalized trapezoidal fuzzy numbers. But we find that the result of Wen and Zhou’s method is unreasonable. Even though Wen and Zhou use exponential distance e-|xA~*-xB~*| to solve the unreasonable negative result caused by geometric distance, the problem caused by exponential distance e-|xA~*-xB~*| is obvious: when |xA~*-xB~*|>1, we always obtain unreasonable result S(A~,B~)<0.37 though coincidence of the area is large and the perimeter and area of generalized trapezoidal fuzzy numbers are almost same.
Based on Wen and Zhou’s method, this paper proposed a new similarity measure combines concepts of the geometric distance, the center of gravity point, the perimeter and the area of generalized trapezoidal fuzzy numbers for calculating the degree of similarity between generalized trapezoidal fuzzy numbers. Some properties of the proposed similarity measure have been proved, and 12 sets of generalized fuzzy numbers have been used to compare the calculation results of the proposed similarity measures with the existing similarity measures. Comparison results indicate that the proposed similarity measure can overcome the drawbacks of the existing similarity measures. Finally, a fault diagnosis experiment is carried out in laboratory based on multifunctional flexible rotor experiment bench to demonstrate the proposed similarity measure. Based on the proposed similarity measure, we present a new fuzzy fault diagnosis algorithm for dealing with rotor fault diagnosis problem, where the values of the evaluating items are represented by generalized trapezoidal fuzzy numbers. The proposed method provides us a useful way for handling fault diagnosis problem, and experimental results demonstrate that the proposed similarity measure is more effective than other methods in terms of rotor fault diagnosis.
The rest of this paper is organized as follows. In Section 2, we briefly review basic concepts of generalized fuzzy numbers. In Section 3, we briefly review some existing similarity measures of fuzzy numbers. In Section 4, we present a new similarity measure between nonstandardized generalized trapezoidal fuzzy numbers. In Section 5, we make a comparison of the calculation results of the proposed similarity measure with existing similarity measures. In Section 6, we apply the proposed similarity measure to deal with rotor fault diagnosis problem. The conclusions are discussed in Section 7.
2. Basic Concepts of Generalized Fuzzy Numbers
In this section, we review basic concepts of generalized fuzzy numbers from Chen [22]. Chen represented generalized trapezoidal fuzzy number A~=(a1,a2,a3,a4,wA~), where a1, a2, a3 and a4 denote real values, 0≤a1≤a2≤a3≤a4≤1 and 0<wA~≤1, as shown in Figure 1. The generalized fuzzy number A~ is a fuzzy subset of the real line R. The membership function μA~ of a generalized fuzzy number A~ satisfies the following conditions:
μA~ is a continuous mapping from the universe of discourse R to the closed interval in [0,1],
μA~(x)=0, where -∞<x≤a1,
μA~(x) is monotone increasing in [a1,a2],
μA~(x)=wA~, where a2≤x≤a3,
μA~(x) is monotone decreasing in [a3,a4],
μA~(x)=0, where a4<x≤∞.
A generalized trapezoidal fuzzy number.
If wA~=1, then the generalized fuzzy number A~ is a normal trapezoidal fuzzy number, denoted as A~=(a1,a2,a3,a4). If a1=a2 and a3=a4, then the generalized fuzzy number A~ is called a crisp interval. If a1<a2=a3<a4, then A~ is normal triangular fuzzy number. If a1<a2<a3<a4, then A~ is a generalized trapezoidal fuzzy number. If a1=a2=a3=a4 and wA~=1, then A~ is a real number.
3. Existing Similarity Measure between Fuzzy Numbers
In this section, some works have to be done to review some existing similarity measures between fuzzy numbers from Chen [14], Lee [15], S.-J. Chen and S.-M. Chen [4], Wei and Chen [20], and Wen et al. [21]. Consider two trapezoidal fuzzy numbers A~ and B~, where A~=(a1,a2,a3,a4,wA~) and B~=(b1,b2,b3,b4,wB~).
3.1. Chen’s Similarity Measures between Fuzzy Numbers
Chen [14] presented a distance-based similarity measure between two trapezoidal fuzzy numbers A~ and B~, where A~=(a1,a2,a3,a4) and B~=(b1,b2,b3,b4). The degree of similarity S(A~,B~) between A~ and B~ is calculated as follows:
(1)S(A~,B~)=1-∑i=14|ai-bi|4,
where S(A~,B~)∈[0,1], if A~ and B~ are triangular fuzzy numbers, where A~=(a1,a2,a3) and B~=(b1,b2,b3). The degree of similarity S(A~,B~) between A~ and B~ is calculated as follows:
(2)S(A~,B~)=1-∑i=13|ai-bi|3.
The larger the value of S(A~,B~), the more similarity between the fuzzy numbers A~ and B~.
3.2. Lee’s Similarity Measures between Fuzzy Numbers
Lee [15] presented a similarity measure between two trapezoidal fuzzy numbers A~ and B~, where A~=(a1,a2,a3,a4) and B~=(b1,b2,b3,b4). The degree of similarity S(A~,B~) between A~ and B~ is calculated as follows:
(3)S(A~,B~)=1-∥A~-B~∥lp∥U∥×4-1/p,
where
(4)∥A~-B~∥lp=(∑i=14(|ai-bi|)p)1/p,∥U∥=max(U)-min(U).
3.3. Chen and Chen’s Similarity Measures between Fuzzy Numbers
Chen’s and Lee’s similarity measures cannot deal with generalized fuzzy numbers. In order to solve this problem, S.-J. Chen and S.-M. Chen [4] presented a similarity measure between two generalized trapezoidal fuzzy numbers A~=(a1,a2,a3,a4;wA~) and B~=(b1,b2,b3,b4;wB~) by using the simple center of gravity method (SCGM). The degree of similarity S(A~,B~) between A~ and B~ is calculated as follows:
(5)S(A~,B~)=[1-∑i=14|ai-bi|4]×(1-|xA~*-xB~*|)B(SA~,SB~)×min(yA~*,yB~*)max(yA~*,yB~*),
where S(A~,B~)∈[0,1], and B(SA~,SB~) is defined as follows:
(6)B(SA~,SB~)={1,ifSA~+SB~>0,0,ifSA~+SB~=0,
where SA~ and SB~ are defined as follows:
(7)SA~=a4-a1,SB~=b4-b1.
Then, the center of gravity (COG) point (xA~*,yA~*) and (xB~*,yB~*) of generalized trapezoidal fuzzy numbers A~ and B~ can be calculated as follows:
(8)yA~*={wA~×((a3-a2)/(a4-a1)+2)6,ifa4≠a1,0<wA~≤1wA~2,ifa4=a1,0<wA~≤1xA~*=yA~*(a3+a2)+(a4+a1)(wA~-yA~*)2wA~.
The values of xB~* and yB~* can be calculated as the same.
3.4. Wei and Chen’s Similarity Measures between Fuzzy Numbers
Wei and Chen [20] presented a new method by combining concepts of the geometric distance, the perimeter, and the height of generalized fuzzy numbers to calculate the degree of similarity S(A~,B~) between A~ and B~, where A~=(a1,a2,a3,a4;wA~), B~=(b1,b2,b3,b4;wB~), 0≤a1≤a2≤a3≤a4≤1, and 0≤b1≤b2≤b3≤b4≤1. The degree of similarity S(A~,B~) between A~ and B~ is calculated as follows:
(9)S(A~,B~)=[1-∑i=14|ai-bi|4]×min(P(A~),P(B~))+min(wA~,wB~)max(P(A~),P(B~))+max(wA~,wB~),
where S(A~,B~)∈[0,1], and the perimeters P(A~) and P(B~) are defined as follows:
(10)P(A~)=(a1-a2)2+wA~2+(a3-a4)2+wA~2+(a3-a2)+(a4-a1),P(B~)=(b1-b2)2+wB~2+(b3-b4)2+wB~2+(b3-b2)+(b4-b1).
3.5. Wen and Zhou’s Similarity Measures between Fuzzy Numbers
All of former four similarity measures are only designed for standardized fuzzy numbers, that is, A~=(a1,a2,a3,a4;wA~) and 0≤a1≤a2≤a3≤a4≤1. They cannot be used to calculate the degree of similarity between nonstandardized fuzzy numbers (i.e., A~=(a1,a2,a3,a4;wA~) and 0≤a1≤a2≤a3≤a4). This is because they all use geometric distance ∑i=14|ai-bi|/4 to measure the degree of similarity S(A~,B~), which may lead to unreasonable negative result easily, that is, S(A~,B~)<0. However, nonstandardized fuzzy numbers are very common in reality. Therefore, it is necessary to put forward a new similarity measure. In order to deal with both standardized and nonstandardized fuzzy numbers, Wen et al. [21] put forward a new method combines concepts of the exponential distance, the perimeter, and the area of the generalized trapezoidal fuzzy numbers for calculating the degree of similarity S(A~,B~) between A~ and B~, where A~=(a1,a2,a3,a4;wA~), B~=(b1,b2,b3,b4;wB~), 0≤a1≤a2≤a3≤a4, and 0≤b1≤b2≤b3≤b4. The degree of similarity S(A~,B~) between A~ and B~ is calculated as follows:
(11)S(A~,B~)=e-|xA~*-xB~*|×min(P(A~),P(B~))+min(a(A~),a(B~))max(P(A~),P(B~))+max(a(A~),a(B~)),
where S(A~,B~)∈[0,1], the center of gravity point (xA~*,yA~*) is calculated as (8), and the perimeters P(A~) and the area a(A~) of generalized trapezoidal fuzzy number A~ are calculated as follows:
(12)P(A~)={(a1-a2)2+wA~2+(a3-a4)2+wA~2+(a3-a2)+(a4-a1),ifa4-a1≠0wA~,ifa4-a1=0a(A~)=12×[(a2-a1)+(a4-a3)]×wA~+(a3-a2)×wA~.
The center of gravity (COG) point (xB~*,yB~*), the perimeter P(B~), and the area a(B~) of generalized trapezoidal fuzzy number B~ can be calculated as the same.
Even though Wen and Zhou use exponential distance e-|xA~*-xB~*| to solve the unreasonable negative result caused by geometric distance ∑i=14|ai-bi|/4, the problem caused by exponential distance e-|xA~*-xB~*| is obvious: the effect of exponential distance on S(A~,B~) is more than the perimeter and the area of the generalized trapezoidal fuzzy numbers. When |xA~*-xB~*|>1, we always obtain the unreasonable result S(A~,B~)<0.37, though the coincidence of the area is large and the perimeter and area of generalized trapezoidal fuzzy numbers are almost the same. Therefore, it is necessary to put forward a new similarity measure to solve the existing defects of Wen and Zhou’s similarity measure.
4. A New Similarity Measure between Generalized Trapezoidal Fuzzy Numbers
In this section, a new method is presented to calculate the degree of similarity between generalized fuzzy numbers based on Wen and Zhou’s method. The proposed method combines concepts of the geometric distance, the center of gravity point, the perimeter, and the area of the generalized trapezoidal fuzzy numbers for calculating the degree of similarity S(A~,B~) between A~ and B~, where A~=(a1,a2,a3,a4;wA~), B~=(b1,b2,b3,b4;wB~), 0≤a1≤a2≤a3≤a4, 0≤b1≤b2≤b3≤b4, 0≤wA~≤1, and 0≤wB~≤1. The proposed method not only can deal with standardized fuzzy numbers but also can deal with nonstandardized fuzzy numbers. In order to deal with nonstandardized fuzzy numbers, normalization process is used to transform it into standardized fuzzy numbers which helps to avoid obtaining unreasonable negative result S(A~,B~)<0. Normalization process is shown as follows:
(13)a1′=a1-min(a1,b1)max(a4,b4)-min(a1,b1),a2′=a2-min(a1,b1)max(a4,b4)-min(a1,b1),a3′=a3-min(a1,b1)max(a4,b4)-min(a1,b1),a4′=a4-min(a1,b1)max(a4,b4)-min(a1,b1).
Then, the new standardized fuzzy number A′~ can be got, where A′~=(a1′,a2′,a3′,a4′;wA~) and 0≤a1′≤a2′≤a3′≤a4′≤1. The new standardized fuzzy number B′~ can be got in the same way, where B′~=(b1′,b2′,b3′,b4′;wB~) and 0≤b1′≤b2′≤b3′≤b4′≤1. The degree of similarity S(A~,B~) between A~ and B~ is calculated as follows:
(14)S(A~,B~)=S(A′~,B′~)=[1-∑i=14|ai′-bi′|4](1-|xA′~*-xB′~*|)×min(P(A′~),P(B′~))+min(a(A′~),a(B′~))max(P(A′~),P(B′~))+max(a(A′~),a(B′~)).
The normalization process is unworkable when the generalized trapezoidal fuzzy numbers A~ and B~ are real numbers at the same time, where A~=(a1,a2,a3,a4;wA~), B~=(b1,b2,b3,b4;wB~), a1=a2=a3=a4, b1=b2=b3=b4, 0≤wA~≤1, and 0≤wB~≤1. In order to deal with this problem, the geometric distance should be transformed into exponential distance. The degree of similarity S(A~,B~) between A~ and B~ is calculated as follows:
(15)S(A~,B~)=e-|(∑i=14|ai-bi|/4)×(1-|xA~*-xB~*|)|×min(P(A~),P(B~))+min(a(A~),a(B~))max(P(A~),P(B~))+max(a(A~),a(B~)),
where S(A~,B~)∈[0,1], and the larger the value of S(A~,B~), the more similarity between the fuzzy numbers A~ and B~.
4.1. Properties of the New Similarity Measure
We can verify the following properties of the new similarity measure.
Property 1.
S(A~,B~)∈[0,1].
Proof.
When the generalized trapezoidal fuzzy numbers A~ and B~ are not real numbers at the same time, and because(1-(∑i=14|ai′-bi′|/4))∈[0,1], (1-|xA′~*-xB′~*|)∈[0,1], min(P(A′~),P(B′~))≤max(P(A′~),P(B′~)), and min(a(A′~),a(B′~))≤max(a(A′~),a(B′~)), we can obtain S(A~,B~)∈[0,1].
When the generalized trapezoidal fuzzy numbers A~ and B~ are real numbers at the same time, and because e-|(∑i=14|ai-bi|/4)×(1-|xA~*-xB~*|)|∈[0,1], min(P(A~),P(B~))≤max(P(A~),P(B~)), and min(a(A~),a(B~))≤max(a(A~),a(B~)), we can obtain S(A~,B~)∈[0,1].
The proof is completed.
Property 2.
Two generalized trapezoidal fuzzy numbers A~ and B~ are identical if and only if S(A~,B~)=1.
Proof.
When the generalized trapezoidal fuzzy numbers A~ and B~ are not real numbers at the same time, if A~ and B~ are identical, then a1=b1, a2=b2, a3=b3, a4=b4 and wA~=wB~. Thus, a1′=b1′, a2′=b2′, a3′=b3′, a4′=b4′, min(P(A′~),P(B′~))=max(P(A′~),P(B′~)), min(a(A′~),a(B′~))=max(a(A′~),a(B′~)), and xA′~*=xB′~*. The degree of similarity S(A~,B~) between A~ and B~ is calculated as follows:
(16)S(A~,B~)=S(A′~,B′~)=[1-∑i=14|ai′-bi′|4](1-|xA′~*-xB′~*|)×min(P(A′~),P(B′~))+min(a(A′~),a(B′~))max(P(A′~),P(B′~))+max(a(A′~),a(B′~))=(1-0)1×1=1.
When the generalized trapezoidal fuzzy numbers A~ and B~ are real numbers at the same time, if A~ and B~ are identical, thus, a1=a2=a3=a4=b1=b2=b3=b4, min(P(A~),P(B~))=max(P(A~),P(B~)), min(a(A~),a(B~))=max(a(A~),a(B~)), and xA~*=xB~*. The degree of similarity S(A~,B~) between A~ and B~ is calculated as follows:
(17)S(A~,B~)=e-|(∑i=14|ai-bi|/4)×(1-|xA~*-xB~*|)|×min(P(A~),P(B~))+min(a(A~),a(B~))max(P(A~),P(B~))+max(a(A~),a(B~))=e-|0×1|=1.
When the generalized trapezoidal fuzzy numbers A~ and B~ are not real numbers at the same time, if S(A~,B~)=S(A′~,B′~)=1, then
(18)S(A~,B~)=S(A′~,B′~)=[1-∑i=14|ai′-bi′|4](1-|xA′~*-xB′~*|)×min(P(A′~),P(B′~))+min(a(A′~),a(B′~))max(P(A′~),P(B′~))+max(a(A′~),a(B′~))=1.
It implies that a1′=b1′, a2′=b2′, a3′=b3′, a4′=b4′, min(P(A′~),P(B′~))=max(P(A′~),P(B′~)), min(a(A′~),a(B′~))=max(a(A′~),a(B′~)), and xA′~*=xB′~*. Therefore, a1=b1,a2=b2,a3=b3,a4=b4, and wA~=wB~, the generalized trapezoidal fuzzy numbers A~ and B~ are identical.
When the generalized trapezoidal fuzzy numbers A~ and B~ are real numbers at the same time, if S(A~,B~)=S(A′~,B′~)=1, then
(19)S(A~,B~)=e-|(∑i=14|ai-bi|/4)×(1-|xA~*-xB~*|)|×min(P(A~),P(B~))+min(a(A~),a(B~))max(P(A~),P(B~))+max(a(A~),a(B~))=1.
It implies that a1=b1, a2=b2, a3=b3, a4=b4, min(P(A~),P(B~))=max(P(A~),P(B~)), min(a(A~),a(B~))=max(a(A~),a(B~)), and xA~*=xB~*. Therefore, a1=b1,a2=b2,a3=b3,a4=b4, and wA~=wB~, the generalized trapezoidal fuzzy numbers A~ and B~ are identical.
The proof is completed.
Property 3.
S(A~,B~)=S(B~,A~).
Proof.
When the generalized trapezoidal fuzzy numbers A~ and B~ are not real numbers at the same time,
(20)S(A~,B~)=S(A′~,B′~)=[1-∑i=14|ai′-bi′|4](1-|xA′~*-xB′~*|)×min(P(A′~),P(B′~))+min(a(A′~),a(B′~))max(P(A′~),P(B′~))+max(a(A′~),a(B′~)),S(A~,B~)=S(B′~,A′~)=[1-∑i=14|bi′-ai′|4](1-|xB′~*-xA′~*|)×min(P(B′~),P(A′~))+min(a(B′~),a(A′~))max(P(B′~),P(A′~))+max(a(B′~),a(A′~)).
Because ∑i=14|ai′-bi′|=∑i=14|bi′-ai′|, based on (12), we can obtain that min(P(A′~),P(B′~))=min(P(B′~),P(A′~)), max(P(A′~),P(B′~))=max(P(B′~),P(A′~)), min(a(A′~),a(B′~))=min(a(B′~),a(A′~)), max(a(A′~),a(B′~))=max(a(B′~),a(A′~)), and |xA′~*-xB′~*|=|xB′~*-xA′~*|. Therefore, S(A~,B~)=S(B~,A~).
When the generalized trapezoidal fuzzy numbers A~ and B~ are real numbers at the same time,
(21)S(A~,B~)=e-|(∑i=14|ai-bi|/4)×(1-|xA~*-xB~*|)|×min(P(A~),P(B~))+min(a(A~),a(B~))max(P(A~),P(B~))+max(a(A~),a(B~)),S(B,~A~)=e-|(∑i=14|bi-ai|/4)×(1-|xB~*-xA~*|)|×min(P(B~),P(A~))+min(a(B~),a(A~))max(P(B~),P(A~))+max(a(B~),a(A~)).
Because ∑i=14|ai-bi|=∑i=14|bi-ai|, based on (12), we can obtain that min(P(A~),P(B~))=min(P(B~),P(A~)), max(P(A~),P(B~))=max(P(B~),P(A~)), min(a(A~),a(B~))=min(a(B~),a(A~)), max(a(A~),a(B~))=max(a(B~),a(A~)), and |xA~*-xB~*|=|xB~*-xA~*|. Therefore, S(A~,B~)=S(B,~A~).
The proof is completed.
5. The Comparison of the Similarity Measures
In this section, the proposed similarity measure is compared with five existing methods (Chen’s method, Lee’s method, Chen and Chen’s method, Wei and Chen’s method, and Wen and Zhou’s method) using 12 sets of generalized fuzzy numbers shown in Figure 2. Set 1 to set 8 are 8 sets of standardized generalized fuzzy numbers, and set 9 to set 12 are 4 sets of nonstandardized generalized fuzzy numbers. Table 1 compares the calculation results of all six similarity measures.
The comparison of the calculation results of the proposed similarity measure and the existing similarity measures.
Sets
Chen’s method
Lee’s method
Chen and Chen’s method
Wei and Chen’s method
Wen and Zhou’s method
The proposed method
Set 1
0.7000
0.4965
0.4200
0.6820
0.7018
0.6567
Set 2
0.7000
0.5000
0.4900
0.7000
0.7408
0.7071
Set 3
1
1
0.8000
0.8248
0.8324
0.8763
Set 4
0.7000
0.3675
0.4900
0.6222
0.5875
0.4743
Set 5
0.7000
0.2500
0.4900
0.7000
0.7408
0.7071
Set 6
1
1
0.7000
0.7209
0.7291
0.8218
Set 7
0.9500
0.6464
0.9048
0.6215
0.6061
0.2810
Set 8
1
#
0.5000
0.5000
0.5000
0.5000
Set 9
−2
0.5000
4
−2
0.0498
0.7071
Set 10
−4
0.9180
22.6667
−3.3481
0.0011
0.8736
Set 11
−4
0.9412
22.6667
−3.5314
0.0011
0.9079
Set 12
−1
0
1
−1
0.1353
0.1353
Note: “#” means that the similarity measure cannot calculate the degree of similarity between two generalized fuzzy numbers.
“bold” means incorrect results.
12 sets of generalized fuzzy numbers.
Some drawbacks of the existing similarity measures are shown in Table 1, which are described as follows.
Set 1 and set 2 of Figure 2 are clearly two different sets of generalized fuzzy numbers. The fuzzy numbers A~ and B~ in set 2 are more similar than set 1. However, according to Table 1, Chen’s method obtains the same degree of similarity S(A~,B~) for set 1 and set 2, which is an incorrect result.
Set 3 in Figure 2 clearly indicates that the generalized fuzzy numbers A~ and B~ are not the same, that is, S(A~,B~)≠1. However, according to Table 1, Chen’s Method and Lee’s method obtain the incorrect results, that is, S(A~,B~)=1.
Set 4 and set 5 of Figure 2 are clearly two different sets of generalized fuzzy numbers. The fuzzy numbers A~ and B~ in set 5 are more similar than set 4. However, according to Table 1, Chen’s method and Chen and Chen’s method obtain the same degree of similarity S(A~,B~) for set 4 and set 5, which are incorrect results. Furthermore, Lee’s method obtains an incorrect result, in which the degree of similarity S(A~,B~) for set 4 is higher than that of set 5.
Set 6 in Figure 2 clearly indicates that the generalized fuzzy numbers A~ and B~ are not the same, that is, S(A~,B~)≠1. However, according to Table 1, Chen’s method and Lee’s method obtain the incorrect results, that is, S(A~,B~)=1.
Set 4 and set 5 of Figure 2 are clearly two different sets of generalized fuzzy numbers. The fuzzy numbers A~ and B~ in set 5 are more similar than set 4. However, according to Table 1, Chen and Chen’s method obtains an incorrect result, in which the degree of similarity S(A~,B~) for set 7 is higher than that of set 6.
Set 8 in Figure 2 clearly indicates that the degree of similarity between two real values cannot be calculated by Lee’s method because the denominator ∥U∥=max(U)-min(U)=0, which results in the incorrect result S(A~,B~)=∞. Furthermore, set 8 in Figure 2 clearly indicates that the generalized fuzzy numbers A~ and B~ are not the same, that is, S(A~,B~)≠1. However, Chen’s method obtains the incorrect results, that is, S(A~,B~)=1.
Set 9 in Figure 2 clearly indicates that Chen’s method, Chen and Chen’s method and Wei and Chen’s method cannot deal with nonstandardized fuzzy numbers, which all give the incorrect result, that is, S(A~,B~)∉[0,1]. Moreover, set 9 has the same degree of similarity S(A~,B~) as set 2, because the geometric distance, the center of gravity point, the perimeter, and the area of the generalized trapezoidal fuzzy numbers of the generalized trapezoidal fuzzy numbers in set 9 should be equal to those in set 2 after normalization process. However, Wen and Zhou’s Method obtains an incorrect result, in which the degree of similarity S(A~,B~) for set 9 is particularly lower than that of set 2.
Set 10 and set 11 of Figure 2 clearly indicates that Chen’s method, Chen and Chen’s Method and Wei and Chen’s method cannot deal with nonstandardized fuzzy numbers, which all give the incorrect result, that is, S(A~,B~)∉[0,1]. Moreover, set 10 and set 11 of Figure 2 are clearly two different sets of nonstandardized generalized fuzzy numbers. Nonstandardized generalized fuzzy numbers A~ and B~ in set 11 are more similar than set 10. However, according to Table 1, Wen and Zhou’s method obtains the same degrees of similarity S(A~,B~) for set 10 and set 11. Finally, set 10 and set 11 in Figure 2 clearly indicates that nonstandardized fuzzy numbers A~ and B~ are approximately similar, that is, S(A~,B~)>0.5. But Wen and Zhou’s method obtains the incorrect results, that is, S(A~,B~)≪0.5.
Set 12 in Figure 2 clearly indicates that Chen’s method, Lee’s method, and Wei and Chen’s method cannot deal with nonstandardized fuzzy numbers, which all give the incorrect result, that is, S(A~,B~)∉[0,1]. Moreover, set 12 in Figure 2 clearly indicates that nonstandardized generalized fuzzy numbers A~ and B~ are not same, that is, S(A~,B~)≠1. However, Chen and Chen’s method obtains the incorrect results, that is, S(A~,B~)=1.
From Table 1 and Figure 2, it clearly indicates that the proposed similarity measure can overcome the drawbacks of the existing similarity measures.
6. Rotor Fault Diagnosis Based on the Proposed Similarity Measure
In this section, the proposed similarity measure is used to deal with the rotor fault diagnosis problem. The major cause of machine vibration is rotor unbalance and misalignment, which leads to additional dynamic load and to accelerate machine deterioration [23]. The proposed similarity measure is a very helpful method in diagnosing the rotor unbalance and misalignment to avoid any failures or damages that may arise.
The proposed similarity measure will be verified by fault diagnosis experiments on multifunctional flexible rotor experiment bench. The vibration acceleration and vibration displacement sensors are, respectively, installed on supporting rotor bases to collect vibration signals by data collector. Experimental procedures involve following main steps.
6.1. Determine the Rotor Fault Diagnosis Set
Rotor unbalance and misalignment are two major rotor fault features. We determine the rotor fault diagnosis set as F={X,Y}, where X=rotor unbalance and Y=rotor misalignment.
6.2. Structure the Fault Template of Rotor Fault Feature
Rotor speed is 1500 rpm, and we choose 1X, 2X, and 3X amplitude frequency response (fundamental frequency is 25 Hz) and average amplitude of vibration displacement as four types of rotor fault feature parameters. In order to structure the fault template of rotor fault feature, we adopt experimental statistics method to simulate rotor fault operation mode and collect five sets of data for each type of rotor fault feature parameter. Every set of data consists of 40 observations. Let Fij = [fj,1,fj,2,…,fj,40] indicate the 40 observations of the j set of data within 16s, where F={X,Y} represent two rotor fault features, respectively, i=1,2,3,4 represent four types of rotor fault feature parameters, respectively, and j=1,2,3,4,5 represent the set number of the collected data. The arithmetic average Mj and mean deviation σj of the j set of data can be got as follows:
(22)Mj=(fj,1+fj,2+⋯+fj,40)40,σj=((fj,1-Mj)2+(fj,2-Mj)2+⋯+(fj,40-Mj)2)40.
Then, the generalized trapezoidal fuzzy number F~ij=(f1,f2,f3,f4,wF~) can be constructed based on these collected data as follows:
(23)f1=min(fj,k),f2=Mj-σj,f3=Mj+σj,f4=max(fj,k),wF~=1.
We can obtain 40 generalized trapezoidal fuzzy numbers F~ij, where F~ij={X~ij,Y~ij}, i=1, 2, 3, 4 and j=1, 2, 3, 4, 5. The result of 40 generalized trapezoidal fuzzy numbers F~ij is shown in Table 2.
The result of 40 generalized trapezoidal fuzzy numbers F~ij.
Fault feature parameters
Rotor unbalance
Rotor misalignment
Generalized trapezoidal fuzzy number
Generalized trapezoidal fuzzy number
1X
X~11
(0.1485, 0.1517, 0.1770, 0.2006; 1)
Y~11
(0.1567, 0.1600, 0.1793, 0.1834; 1)
X~12
(0.1518, 0.1523, 0.1669, 0.1820; 1)
Y~12
(0.1644, 0.1698, 0.1787, 0.1809; 1)
X~13
(0.1629, 0.1635, 0.1653, 0.1656; 1)
Y~13
(0.1625, 0.1794, 0.2069, 0.2038; 1)
X~14
(0.1606, 0.1611, 0.1622, 0.1630; 1)
Y~14
(0.1841, 0.1878, 0.1948, 0.1958; 1)
X~15
(0.1576, 0.1588, 0.1607, 0.1615; 1)
Y~15
(0.1736, 0.1773, 0.1835, 0.1857; 1)
2X
X~21
(0.1210, 0.1426, 0.1651, 0.1665; 1)
Y~21
(0.3071, 0.3158, 0.3375, 0.3466; 1)
X~22
(0.1343, 0.1457, 0.1560, 0.1563; 1)
Y~22
(0.3079, 0.3196, 0.3361, 0.3388; 1)
X~23
(0.1420, 0.1444, 0.1491, 0.1534; 1)
Y~23
(0.3095, 0.3269, 0.3499, 0.3507; 1)
X~24
(0.1438, 0.1452, 0.1486, 0.1512; 1)
Y~24
(0.3215, 0.3287, 0.3413, 0.3510; 1)
X~25
(0.1382, 0.1402, 0.1453, 0.1488; 1)
Y~25
(0.3132, 0.3145, 0.3228, 0.3291; 1)
3X
X~31
(0.0899, 0.1065, 0.1261, 0.1296; 1)
Y~31
(0.2338, 0.2522, 0.3022, 0.3218; 1)
X~32
(0.1044, 0.1074, 0.1116, 0.1131; 1)
Y~32
(0.2485, 0.2631, 0.2820, 0.2919; 1)
X~33
(0.0990, 0.1026, 0.1101, 0.1150; 1)
Y~33
(0.1865, 0.1878, 0.2556, 0.2991; 1)
X~34
(0.1059, 0.1095, 0.1138, 0.1158; 1)
Y~34
(0.2022, 0.2078, 0.2185, 0.2250; 1)
X~35
(0.1137, 0.1166, 0.1199, 0.1211; 1)
Y~35
(0.2031, 0.2120, 0.2390, 0.2595; 1)
Displacement
X~41
(4.1820, 4.1933, 4.4181, 4.6660; 1)
Y~41
(4.7450, 4.7943, 5.1707, 5.3920; 1)
X~42
(4.3740, 4.3917, 4.4369, 4.4650; 1)
Y~42
(4.3910, 4.4618, 4.7070, 4.8610; 1)
X~43
(4.3190, 4.3452, 4.3799, 4.3890; 1)
Y~43
(4.0940, 4.0343, 4.8373, 5.2340; 1)
X~44
(4.2800, 4.2889, 4.3387, 4.4000; 1)
Y~44
(4.6650, 4.7705, 5.0505, 5.1510; 1)
X~45
(4.2870, 4.2962, 4.3662, 4.4620; 1)
Y~45
(4.4510, 4.4953, 4.6309, 4.6770; 1)
The generalized trapezoidal fuzzy number for each fault feature parameter can be got as follows:
(24)F~i=∑j=15F~ij5.
The result of generalized trapezoidal fuzzy numbers for the fault template is shown in Table 3.
The fault template of rotor fault feature.
Fault feature
Fault feature parameters
Generalized trapezoidal fuzzy numbers
Rotor unbalance X
1X
X~1= (0.1563, 0.1575, 0.1664, 0.1745; 1)
2X
X~2= (0.1359, 0.1436, 0.1528, 0.1552; 1)
3X
X~3= (0.1026, 0.1085, 0.1163, 0.1189; 1)
Displacement
X~4= (4.2884, 4.3031, 4.3880, 4.4764; 1)
Rotor misalignment Y
1X
Y~1= (0.1683, 0.1749, 0.1886, 0.1899; 1)
2X
Y~2= (0.3118, 0.3211, 0.3375, 0.3432; 1)
3X
Y~3= (0.2148, 0.2246, 0.2595, 0.2795; 1)
Displacement
Y~4= (4.4692, 4.5112, 4.8793, 5.0630; 1)
6.3. Structure the Diagnosing Rotor Fault Feature
We let multifunctional flexible rotor experiment bench bed run in rotor misalignment mode Y and collect data for each parameter of rotor fault feature. The generalized trapezoidal fuzzy numbers A~i of diagnosing rotor fault feature can be obtained through the collected data. The result is shown in Table 4.
The diagnosing rotor fault feature.
Fault feature parameters
Generalized trapezoidal fuzzy numbers
1X
A~1= (0.1515, 0.1550, 0.1743, 0.1784; 1)
2X
A~2= (0.3221, 0.3308, 0.3525, 0.3611; 1)
3X
A~3= (0.1338, 0.1522, 0.2021, 0.2214; 1)
Displacement
A~4= (4.2450, 4.2943, 4.6707, 4.8910; 1)
6.4. Calculate the Degree of Similarity between Fuzzy Numbers A~i and Each Fault Feature Parameter
In the following, Wen and Zhou’s and the proposed methods are used to deal with the rotor fault diagnosis problem, respectively. Based on (11) and (14), the degree of similarity between generalized trapezoidal fuzzy numbers A~i and each fault feature parameter shown in Table 3 can be calculated. The calculating rotor fault diagnosis results of Wen and Zhou’s and the proposed methods are shown in Tables 5 and 6, respectively.
Rotor fault diagnosis result of Wen and Zhou’s method.
Fault feature parameters
The degree of similarity
Rotor fault diagnosis result
Rotor unbalance X
Rotor misalignment Y
1X
S(A~1,X~1)=0.9853
S(A~1,Y~1)=0.9769
Rotor unbalance
2X
S(A~2,X~2)=0.8038
S(A~2,Y~2)=0.9776
Rotor misalignment
3X
S(A~3,X~3)=0.8638
S(A~3,Y~3)=0.9105
Rotor misalignment
Displacement
S(A~4,X~4)=0.5756
S(A~4,Y~4)=0.7930
Rotor misalignment
Rotor fault diagnosis result of the proposed method.
Fault feature parameters
The degree of similarity
Rotor fault diagnosis result
Rotor unbalance X
Rotor misalignment Y
1X
S(A~1,X~1)=0.6404
S(A~1,Y~1)=0.6512
Rotor misalignment
2X
S(A~2,X~2)=0.6961
S(A~2,Y~2)=0.7145
Rotor misalignment
3X
S(A~3,X~3)=0.4282
S(A~3,Y~3)=0.6350
Rotor misalignment
Displacement
S(A~4,X~4)=0.4617
S(A~4,Y~4)=0.7797
Rotor misalignment
From Table 5, we obtain S(A~1,X~1)>S(A~1,Y~1), S(A~2,X~2)<S(A~2,Y~2), S(A~3,X~3)<S(A~3,Y~3), and S(A~4,X~4)<S(A~4,Y~4). It means that the rotor fault diagnosis result for 1X fault feature parameter is rotor unbalance and the result for the other three types of rotor fault feature parameters are rotor misalignment. Thus, it is unable to determine fault feature from the rotor fault diagnosis result of Wen and Zhou’s method. Comparing with the result of Wen and Zhou’s method, the rotor fault diagnosis results of the proposed method for four types of rotor fault feature parameters are all rotor misalignment. The rotor fault diagnosis result of the proposed method coincides with the fault feature of multifunctional flexible rotor experiment bench. Therefore, experimental results demonstrate that the proposed similarity measure is more effective than Wen and Zhou’s methods in terms of rotor fault diagnosis.
7. Conclusions
This paper presents a new similarity measure used for fault diagnosis technology, which can deal with both standardized and nonstandardized generalized trapezoidal fuzzy numbers. It combines concepts of the geometric distance, the center of gravity point, the perimeter, and the area of generalized trapezoidal fuzzy numbers for calculating the degree of similarity S(A~,B~) between A~ and B~. Some properties of the proposed similarity measure have been proved, and 12 sets of generalized fuzzy numbers have been used to compare the calculation results of the proposed similarity measures with the existing similarity measures. From Table 1 and Figure 2, it clearly indicates that the proposed similarity method can overcome the drawbacks of existing similarity measures. Finally, a rotor fault diagnosis example has been given to show the proposed similarity measure.
Acknowledgment
The authors would like to thank the anonymous reviewers for critical comments and suggestions leading to an improved presentation.
GhainaniA. T.ZinA. A. M.IsmailN. M.Fuzzy timing petri net for fault diagnosis in power system201220121271719510.1155/2012/717195RobertsC.DassanayakeH. P. B.LehrasabN.GoodmanC. J.Distributed quantitative and qualitative fault diagnosis: railway junction case study20021044194292-s2.0-003653206210.1016/S0967-0661(01)00159-9UbarR.KostinS.RaikJ.Embedded fault diagnosis in digital systems with BIST2008325-62792872-s2.0-5094910799910.1016/j.micpro.2008.03.006ChenS.-J.ChenS.-M.Fuzzy risk analysis based on similarity measures of generalized fuzzy numbers200311145562-s2.0-003733146510.1109/TFUZZ.2002.806316EslamloueyanR.Designing a hierarchical neural network based on fuzzy clustering for fault diagnosis of the Tennessee-Eastman process20111114071415LuoH.YangS. L.HuX. J.HuX. X.Agent oriented intelligent fault diagnosis system using evidence theory20123925242531DuZ.JinX.YangY.Fault diagnosis for temperature, flow rate and pressure sensors in VAV systems using wavelet neural network2009869162416312-s2.0-6344910793010.1016/j.apenergy.2009.01.015MuralidharanV.SugumaranV.Feature extraction using wavelets and classification through decision tree algorithm for fault diagnosis of mono-block centrifugal pump201346353359ZhangJ.Improved on-line process fault diagnosis through information fusion in multiple neural networks200630558571GrbovicM.LiW. C.XuP.UsadiA. K.SongL. M.VuceticS.Decentralized fault detection and diagnosis via sparse PCA based decomposition and Maximum Entropy decision fusion201222738750BarakatM.DruauxF.LefebvreD.KhalilM.MustaphaO.Self adaptive growing neural network classifier for faults detection and diagnosis20117438653876WeerasingheM.Barry GommJ.WilliamsD.Neural networks for fault diagnosis of a nuclear fuel processing plant at different operating points1998622812892-s2.0-034490624710.1016/S0967-0661(97)00003-8XuX. B.FengH. S.WangZ.WenC. L.An information fusion method of fault diagnosis based on interval basic probability assignment20112022552602-s2.0-79958121685ChenS.-J.New methods for subjective mental workload assessment and fuzzy risk analysis19962754494722-s2.0-0030231260LeeH. S.Optimal aggregation method for fuzzy opinions of group decision3Proceedings of the IEEE International Conference on Systems, Man, and Cybernetics (SMC '99)October 19993143192-s2.0-17544405065ChenS.-J.ChenS.-M.A new method to measure the similarity between fuzzy numbersProceedings of the 10th IEEE International Conference on Fuzzy Systems2011Melbourne, Australia208214ChenJ. H.ChenS.-M.A new method for ranking generalized fuzzy numbers for handling fuzzy risk analysis problemsProceedings of the 9th Joint Conference on Information Sciences (JCIS '06)October 2006Kaohsiung, Taiwan119611992-s2.0-3384770980210.2991/jcis.2006.80ChenS.-J.ChenS.-M.Fuzzy risk analysis based on the ranking of generalized trapezoidal fuzzy numbers20072611112-s2.0-3384653857010.1007/s10489-006-0003-5WeiS.-H.ChenS.-M.A new similarity measure between generalized fuzzy numbersProceedings of the 3rd International Conference on Soft Computing and Intelligent System and 7th International Symposium on Advanced Intelligent System2006Tokyo, Japan315320WeiS.-H.ChenS.-M.A new approach for fuzzy risk analysis based on similarity measures of generalized fuzzy numbers20093615895982-s2.0-5384912923410.1016/j.eswa.2007.09.033WenC. L.ZhouZ.XuX. B.A new similarity measure between generalized trapezoidal fuzzy numbers and its application to fault diagnosis2011393162-s2.0-79959574212ChenS.-HRanking generalized fuzzy number with graded mean integration2Proceedings of the 8th international Fuzzy Systems Association World Congress1999Taipei, Taiwan899902ArebiL.GuF.BallA.Rotor Misalignment detection using a wireless sensor and a shaft encoderProceedings of the Computing and Engineering Annual Researchers' Conference (CEARC '10)2010Huddersfield, UKUniversity of Huddersfield613