Complex third-order cumulant has different definition forms. Different forms have different coupling properties, and the generated complex cumulants slices contain different coupling information of signals. In experiments, using the different definitions, the same coupling method is applied to both specific fault signals and normal signals. Furthermore, complex third-order cumulant slices spectrum is defined, and it is used to analyse the coupling features of normal signals and fault signals. Experiments indicate that the detection accuracy rate on the same fault is not the same with the different coupling method, thus, it provides an alternative method to diagnose the specific fault.
Using higher-order cumulant can suppress the gaussian background noise automatically (colored or white), because of this, the unify measurement of higher order cumulant has been brought to the people’s attention and is becoming increasingly a very useful tool of signal processing. In general, the mechanical vibration signals are nonlinear, non-gaussian, the noise in the environment can be approximated as gaussian noise, and if vibration signals are analyzed by the high-order cumulant, it is easier to extract characteristic information related to the workpiece. Therefore, high order cumulant has been widely used in mechanical vibration, fault diagnosis and other fields [
If { x ( n ) } is a zero mean k-order stationary random process, then the k-order cumulative quantity c k x ( τ 1 , τ 2 , ⋯ , τ k − 1 ) of this process is defined as the k-order cumulant of random variable { x ( n ) , x ( n + τ 1 ) , ⋯ , x ( n + τ k − 1 ) } , namely:
c k x ( τ 1 , τ 2 , ⋯ , τ k − 1 ) = c u m { x ( n ) , x ( n + τ 1 ) , x ( n + τ 2 ) , ⋯ , x ( n + τ k − 1 ) } (1)
The k-order moment m k x ( τ 1 , τ 2 , ⋯ , τ k − 1 ) of this process is defined as the k-order joint moment of random variable { x ( n ) , x ( n + τ 1 ) , ⋯ , x ( n + τ k − 1 ) } , namely:
m k x ( τ 1 , τ 2 , ⋯ , τ k − 1 ) = m o m { x ( n ) , x ( n + τ 1 ) , ⋯ , x ( n + τ k − 1 ) } (2)
Here, mom() represents the joint moment, and the third-order cumulant is
c 3 x ( τ 1 , τ 2 ) = E { x ( n ) x ( n + τ 1 ) x ( n + τ 2 ) } (3)
According to literature (7), in Equation (3), x ( n ) is the plural signal and is defined as follows:
Define one:
c 3 x ( τ 1 , τ 2 ) = E { x ( n ) x ( n + τ 1 ) x ( n + τ 2 ) } (4)
Define two:
c 3 x ( τ 1 , τ 2 ) = E { x ∗ ( n ) x ( n + τ 1 ) x ( n + τ 2 ) } (5)
Define three:
c 3 x ( τ 1 , τ 2 ) = E { x ∗ ( n ) x ( n + τ 1 ) x ∗ ( n + τ 2 ) } (6)
where x ∗ ( n ) is the conjugate complex number of x ( n ) , the complex signal in this paper is obtained by Hilbert transformation of the original signal collected. Literature (6), points out that the way to define a, can use the plural form of a harmonic of the quadratic phase coupling signals (such as type (7), the ϕ 3 = ϕ 2 + ϕ 1 , harmonic component ω 3 is by the harmonic component of ω 1 and ω 2 by quadratic phase coupling, and ω 3 = ω 1 + ω 2 ) and third order cumulant to 0, the way to define the second and third order cumulant, as shown in the type (8) to define the three ways of third order cumulant as shown in type (9).
x ( n ) = ∑ i = 1 3 A i exp ( ω i n ) + ϕ i (7)
A 1 A 2 A 3 [ exp [ j ( ω 1 τ 1 + ω 2 τ 2 ) ] + exp [ j ( ω 2 τ 1 − ω 1 τ 2 ) ] ] (8)
A 1 A 2 A 3 [ exp [ j ( ω 3 τ 1 − ω 1 τ 2 ) ] + exp [ j ( ω 3 τ 1 − ω 2 τ 2 ) ] ] (9)
Let τ 2 = 0 in Equations ((8) and (9)) respectively, then Equation (8) becomes:
A 1 A 2 A 3 [ exp [ j ( ω 1 τ 1 ) ] + exp [ j ( ω 2 τ 1 ) ] ] (10)
then Equation (9) becomes:
A 1 A 2 A 3 [ exp [ j ( ω 3 τ 1 ) ] + exp [ j ( ω 3 τ 1 ) ] ] (11)
Can be seen from the type (3), when τ 2 = 0 for the section, is the third order cumulants slices. Formulas (10) and (11) show that when the signal can be expressed in the form of complex harmonics, the third-order cumulant slices obtained by definition 2 and definition 3 contain only the previous-coupled signals and the after-coupled signals, respectively.
The pressure-reducing valve studied in this experiment is the pilot pressure-reducing valve. When a foreign body oil press reducing valve in and out of the mouth, the different pressure will affect the normal operation of the valve’s pressure, in order to get the pressure in operation under the fault state of the signal, this paper tests of artificial set the fault as follows:
Fault 1: add Φ 3 mm core to the relief valve outlet
Fault 2: Φ 3 mm iron core was added at the inlet and outlet of the pressure relief valve.
The working failure condition of pressure relief valve can be approximately simulated by experiments.
In this paper, LabVIEW software and pci-6014 data acquisition card and an acceleration sensor are used to collect vibration signals of pressure relief valve in normal and fault state successively. In each measurement, the oil pressure is divided into five pressure levels from 1 MPa to 5 MPa. The sampling frequency is 250 Hz, the reading frequency is 125 Hz, and the sampling process took about 2 minutes. In this experiment, a total of 54 data sets of 18 groups were collected in the normal working state and two failure states of the pressure relief valve respectively. The number of data used in the experiment is 1536.
Due to the system in the process of testing the influence of the external and internal factors inevitable in the process of the output with unwanted noise, In this paper, the median method is used to pre-process the collected vibration signal and eliminate the interference noise mixed in the signal. The filtered signal is shown in
Fault recognition, this paper in order to will be measured in the experiments of normal state and failure state of the 36 sets of data, respectively, according to the type (4)-(6), and τ 2 = 0 , calculate the amount of plural third-order cumulant slices and take its modulus value, by contrast, let − 8 ≤ τ 1 ≤ 8 and the calculated results are input support vector machine (SVM) to diagnose faults, A more satisfactory result can be obtained. Among them, the complex third-order cumulant slices calculated by definition 1 are shown in
The same LSSVM is established in this paper. When trained LSSVM, respectively the normal state and fault state code is 1 and 1, due to the normal data and fault data are obtained under the five different oil pressure, training respectively in table a normal data and fault data of different hydraulic conditions corresponding 10 groups, 5 sets of data to set up by the least squares support vector
Normal state | 1891,222.1180 | 2,055,712.5807 | 1,950,438.9033 |
---|---|---|---|
724,726.0347 | 549,599.1331 | 525,556.0427 | |
448,594.3487 | 537,601.5387 | 579,072.0499 | |
Failure state | 3,935,434.0292 | 2,980,882.4356 | 3,161,106.0132 |
7,032,668.3428 | 6,825,709.3659 | 5,024,387.5584 | |
9,151,062.1560 | 10,907,623.2398 | 1,291,835.9094 |
machine (SVM) training, training after recycling the support vector machine (SVM) is used to identify the fault, use the network after training for simulating the remnant of 26 sets of data to identify, as shown in
In order to have an intuitive understanding of the properties of third-order cumulant slices in fault diagnosis based on the above 3 definitions, this paper defines three-order cumulant slice spectra as follows:
D i a g ( ω ) = ∑ τ 1 = − ∞ ∞ c 3 x ( τ 1 , 0 ) e − j ω τ 1 (12)
In this experiment, 108 sets of third-order cumulant slices obtained from 36 sets of data in normal state and fault 1 state according to the 3 definitions mentioned above were calculated, and then the complex third-order cumulant slice spectrum according to formula (12) were calculated. In this paper, when the oil pressure is 1 MPa, 3 MPa and 5 MPa respectively in normal state and failure state, a set of third-order cumulant slice spectra of the two states is selected and shown in Figures 3-8.
In the figure, the X-axis and Y-axis represent the frequency in HZ, and the vertical axis represents the normalized amplitude, without dimension. From
Different definitions of complex third-order cumulants determine that they contain different coupling information, which is bound to be reflected in the complex third-order cumulant slices. In this paper, the different coupling characteristics are analyzed by fault diagnosis experiments. It is concluded that the ability of complex third-order cumulant slices with different coupling properties to maintain the characteristics of the same kind signals is different. It also provides an alternative method for fault diagnosis. The fault diagnosis method provided in this paper can be realized through the single-chip computer system, in which a certain amount of data can be stored first, and then the fault diagnosis can be made according to the data collected in real time.
This paper is supported by Research Foundation of the Nanchang Normal University for Doctors (NSBSJJ2018014), National Natural Science Foundation of China (61562063), Science and Technology Project of Jiangxi Provincial Education Department (GJJ171113).
The authors declare no conflicts of interest regarding the publication of this paper.
Wu, W.B. (2018) The Application of Cumulants Slices in Fault Diagnosis. International Journal of Modern Nonlinear Theory and Application, 7, 89-96. https://doi.org/10.4236/ijmnta.2018.73007