Published: 30 September 2016

A novel failure analysis and diagnosis method for hydraulic-turbine generator unit

Xin Xia1
Wei Ni2
1, 2College of Automation, Huaiyin Institute of Technology, Huaian, China
Corresponding Author:
Xin Xia
Views 248
Reads 113
Downloads 1587

Abstract

The failure analysis and diagnosis of hydraulic-turbine generator unit (HGU) is important to protect the safety and stability of the electric system. But failure analysis and diagnosis for HGU is lacking of fault data, and most of the failure analysis and diagnosis methods are proposed without considering the characteristic of HGU. In this paper, a vibration dynamic modeling method for HGU is proposed by using the finite element method, and further the vibration data in different states are getting through numerical simulation. Then, the failure feature is extracted based on the nonlinear output frequency response functions (NOFRFs). Finally, a diagnosis system with Support Vector Machine (SVM) is proposed and employed for diagnosis of the HGU. The experimental results indicate that the feature which extracted from NOFRFs has a strong effect, and demonstrates that the proposed method is feasible and helpful for fault diagnosis in HGU.

1. Introduction

The hydraulic-turbine generator unit (HGU) is part of the key equipment in hydropower station. The health state of the HGU is relevant to the safety and stability of the power supply system. Fault in the HGU will not only lead to economic losses but also threaten the safety of people’s life. So, the failure analysis and diagnosis of HGU play an important role in the operation and maintenance which can reduce or eliminate the accident [1].

According to operating information of the HGU, the failure mainly performs in the form of vibration signals. Researchers have done a lot of work to find the relationship between vibration feature and the failure [1-6]. As hydropower station operating with planned maintenance, some faults are repaired in the early time without record, which makes the vibration data in fault state are difficult to acquire. So there is an urgent need to develop a new method to get faults vibration data of HGU. Some simulation works have been made to get the vibration data. Nelson [7] have proposed using the finite element method to simulate the dynamic of rotor-bearing system. Shouju Li [6] have proposed using the finite element method to identify the vibration load. But there is little research on the simulation of faults vibration data of HGU.

Failure analysis and diagnosis methods of HGU can be classified into two categories: signal analysis and modeling methods. Most signal analysis methods are evolved from the diagnosis methods of rotating machinery. Although the HGU belongs to rotating machinery, but it has its own features, such as: low speed, big weight of rotor, variety of types. The signal analysis methods cannot suit each type of HGU very well. In modeling methods, regarding the modeling approach, it is based on a reliable model of the HGU and can be used for the construction of the symptom-fault correlation, or the input-output relation [8]. Recent research has proposed many modeling methods for fault diagnosis [8-12]. Most of them are based on Volterra series and Generalized Frequency Response Function (GFRF). Volterra series is useful in modeling hydroelectric generator unit, but large number of time domain kernels are difficult for fault diagnosis, and the time domain kernels have no physical meaning. GFRF models always encounter the problem of disaster of dimensionality, and high dimensional GFRF have no intuitive feeling to people. LANG [13-15] have proposed non-linear output frequency response functions (NOFRFs) of Volterra systems. This method simplifies the nonlinear frequency domain model, and makes people have an intuitive feeling. However, their work focus on mechanism analysis and few work has done in intelligent diagnosis based on NOFRFs.

Fig. 1The fault diagnostic procedure of HGU

The fault diagnostic procedure of HGU

In this paper, a novel failure analysis and diagnosis is proposed for HGU. Firstly, in order to overcome the lacking of failure data, a finite element model of a hydroelectric generator unit is proposed, and further the vibration data in different states are acquired based on the model. Secondly, NOFRFs are built for each state, and some quantitative indicators are proposed for NOFRFs as diagnostic features. Finally, a support vector machine (SVM)[16-18] system with binary tree architecture for multi-class classification is built for fault diagnosis. The fault diagnostic procedure is shown in Fig. 1. The proposed method is employed on an experimental study which parameters are totally according to a real HGU. The results indicate that the features which extracted from NOFRFs are much more useful than other diagnosis methods. And the proposed method has a strong effect in failure analysis and diagnosis of HGU.

The rest of paper is organized as follows: in Section 2, the finite element simulation of HGU is presented. In Section 3, the feature extraction and SVM diagnostic system are presented. In Section 4, the proposed method is employed in a HGU which is simplified from a real hydroelectric generator unit, and the proposed is compare to others fault diagnosis methods. Finally, conclusions from this research as well as the advantages and limitations of the proposed method are discussed in Section 5.

2. Finite element simulation of vibration response of HGU

2.1. Dynamic model of HGU

Failure analysis and diagnosis of HGU plays an important role in the operation and maintenance which can reduce or eliminate the accident. As lacking of fault data, less progress has been made in analysis and diagnostic methods for HGU. With the development of finite element simulation technology, the vibration data of HGU in fault states can be acquired through the numerical simulation.

Fig. 2a) Model of HGU, b) Finite element model of HGU

a) Model of HGU, b) Finite element model of HGU

a)

a) Model of HGU, b) Finite element model of HGU

b)

Fig. 3The movement of elastic beam

The movement of elastic beam

Fig. 2(a) shows a HGU, and the structure is quite complex. To simplify the model, the shaft of HGU is set as N elastic beams (N+1 nodes, and the node is numbered on the figure) according to the physical structure. The following assumptions are made: the bearing force is linear; the support and foundation are rigid. Then the HGU can be described as Fig. 2(b). In this paper, the parameters of each element have been set according to a real HGU as show in Table 1.

Table 1The parameters of each element after simplified the HGU

Number
1
2
3
4
5
6
Length (mm)
332
1177
3063.5
1750
1750
1750
Outer diameter (mm)
1170
1170
8664.5
2150
2150
2150
Inner diameter (mm)
730
730
0
1707
1707
1707
Mass (Kg)
1710
6062
980000
18403
18403
18403
Mass per unit length (Kg/mm)
5.15
5.15
319.9
10.516
10.516
10.516
Number
7
8
9
10
11
Length (mm)
1750
1750
750
1750
1462.5
Outer diameter (mm)
2150
2150
2150
2150
6000
Inner diameter (mm)
1707
1707
1707
1707
0
Mass (Kg)
18403
18403
7887
18403
120000
Mass per unit length (Kg/mm)
10.516
10.516
10.516
10.516
82

For each elastic, the moving coordinate system can be set as Fig. 3.

So the movement of each elastic beam can be described as:

1
Xi=xAi yAi xAiθ yAiθ; xBi yBi xBiθ yBiθ.

The moving equation of the whole system can be described under usual assumption of time invariant physical properties, and viscous damping as follow [19, 20]:

2
MX¨+CX˙+KX=F,

where M is the physical mass symmetric mass matrix which includes translating inertial matrix MT and rotating inertial matrix MR, C is the synthetic damped matrix which includes damping matrix, coupling matrix of support and gyro force matrix. K is the stiffness matrices. X is the vector of displacement responses of elastic beams. F is the external excitation forces.

According to Nelson [6], The MT, MR, K and C are set as follow:

MT=Pl420×1560022l5400-13l0156-22l005413l00-22l4l200-13l-3l2022l004l213l00-3l2540013l15600-22l054-13l0015622l0013l-13l20022l4l20-13l00-3l2-22l004l2,
MR=Pr2120l×36003l-36003l036-3l00-36-3l00-3l4l2003l-l203l004l2-3l00-l2-3600-3l3600-3l0-363l00363l00-3l-l2003l4l203l00-l2-3l004l2,
K=EIl3×12006l-12006l012-6l00-12-6l00-6l4l2006l2l206l004l2-6l002l2-1200-6l1200-6l0-126l00126l00-6l2l2006l4l206l002l2-6l004l2,

where l, r, EI and P are the length, radius, elasticity modulus and mass per unit length of elastic beam.

C is Rayleigh damping as follow:

3
C=αM+βK,

where:

4
α=2(ξ2/ωn2-ξ1/ωn1)1/ωn22-1/ωn12β=2(ξ2ωn2-ξ1ωn1)ωn22-ωn12,

where ξ1 and ξ2 are damping coefficient that can be taken as 0.1 and 0.2, respectively. And ωn1, ωn2 are the first and second order critical speed.

It is well know that, some of the nonlinear excitation forces will cause the nonlinear vibration of hydroelectric generator unit. In this paper, we forces on four kinds of typical states of hydroelectric generator unit, as normal state, rubbing state, misalignment state and rubbing-misalignment state. The difference among the models of different states is the external excitation forces. So, the external excitation forces of each states are described as follow:

2.1.1. Normal state

In normal state, the excitation force is mainly caused by the initial unbalance mass:

5
F=Q=Meω2r,

where Q is the initial unbalance mass force, Me is initial unbalance mass, ω is the rotating speed, r is the eccentric radius.

The unbalance mass force is decompose into x and y direction:

6
Qx=Meω2rsin2πω t+φ,Qy=Meω2rcos2πω t+φ,

where φ is the initial phase.

2.1.2. Rubbing state

In rubbing state, because of the discontinuous contacts between rotors and stators, the hydroelectric generator unit present non-linearly phone vibrating phenomena. The excitation forces includes unbalance force and rubbing force:

7
F=Q+Fr.

Han Q. [20] have proposed the calculation method of partial rubbing force, as follows:

8
Frx=-krx-δ0hx-δ0,Fry=-frkrx-δ0hx-δ0,
9
h(x-δ0)=0, x<δ0,1, xδ0,

where x is radial displacement of rotor, δ0 is clearance. For detailed explanations, refer to the research of Han Q. [20].

2.1.3. Misalignment state

In misalignment state, the excitation force includes initial unbalance force and misalignment force:

10
F=Q+Fm,

where Fm is the misalignment force. K. M. Al-Hussain [21] have explained the misalignment force very detailed. In this paper, we refer to his conclusions.

2.1.4. Rubbing-misalignment state

In misalignment state, the excitation force includes unbalance force, misalignment force and rub force, as follow:

11
F=Q+Fm+Fr.

2.2. Numerical simulation

The vibration data of HGU can be get by numerical simulation from Eq. (2). According to the Newmark constant-averaged-acceleration method, Eq. (2) can be discretized into the following scheme [22]:

12
MX¨t+Δt+CX˙t+Δt+KXt+Δt=Ft+Δt,

where:

X˙t+Δt=X˙t+tt+ΔtX¨ttdt, Xt+Δt=Xt+tt+ΔtX˙ttdt,

where Δt is the step time. Substituting Eq. (12) into Eq. (2), then the derivation of acceleration can be get:

13
MX¨t+Δt=-Ctt+ΔtX¨ttdt-Ktt+ΔtX˙ttdt-CX˙t-KXt+Ft+Δt.

The initial conditions are expressed as: Xt=0=0, X˙t=0=0.

The vibration data can be get by iterative computation through Eq. (2) and Eqs. (12), (13).

The parameters of HGU and the failure parameters are set as Table 2. The initial parameters of numerical simulation are as follows: the rubbing location is at the generator rotor, and the misalignment location is at node 4 (refer to Fig. 2). The initial unbalance mass of HGU locates in runner, and the simulation step is 100 steps per round. The vibration data of node 2 of HGU are collected in different states. The wave and spectrum of simulation vibration are shown in Fig. 4 and Fig. 5.

Table 2Parameters of HGU and the failure parameters

Parameters of hydroelectric generator unit
Failure parameters
EI
Speed
Radius of eccentric mass
Rubbing coefficient
Angle of misalignment
2.1×1011 N/m2
142.9 rpm
4000 mm
6×109N/m
1o

3. Fault diagnosis with NOFRFs-SVM method

3.1. Feature extraction based on NOFRFs

HGU belongs to large rotating machinery. The vibration characteristics are similar to many rotor-bearing systems. But it has its own features as well, such as: low speed, big weight of rotor, variety of types. The traditional time-frequency domain analysis methods have difficulty in failure analysis and diagnosis. As seen from Fig. 4 and Fig. 5, the characteristic of different states cannot be described clearly in time-frequency domain. And multi fault is also hard to be explained in the traditional way. So there is an urgent demand for new viewpoints and methods for failure analysis and diagnosis of HGU. In this paper, a novel feature extraction method is proposed based on NOFRFs model for failure analysis and diagnosis.

The vibration model of HGU can be described as a time domain Volterra series, as follow:

14
yt=n=1--hnτ1,τni=1nut-τidτi,

where ut is the input vibration data, i.e. the vibration data of lower bearing. yt is the output vibration data, i.e. the vibration data of upper bearing. hnτ is called the nth order time domain kernel.

The Eq. (14) is changed with multi-dimensional Fourier Transform, and the frequency domain Volterra series is described as follow:

15
Yjω=n=1Nω1++ωn=ωHnjω1,jωni=1nUjωidσnω,

where Yjω and Ujω are the input and output spectra which are the Fourier Transforms of ut and yt. Hn is called the nth order Generalized Frequency Response Function (GFRF), and can be expressed as Eq. (16):

16
Hnjω1,jωn=--hnτ1τne-jω1τ1++ωnτni=1ndτi.

Fig. 4Time-domain vibration of HGU in each state

Time-domain vibration of HGU in each state

a) Normal state

Time-domain vibration of HGU in each state

b) Rubbing state

Time-domain vibration of HGU in each state

c) Misalignment state

Time-domain vibration of HGU in each state

d) Rubbing-misalignment state

Fig. 5Frequency domain vibration of HGU in each state

Frequency domain vibration of HGU in each state

a) Normal state

Frequency domain vibration of HGU in each state

b) Rubbing state

Frequency domain vibration of HGU in each state

c) Misalignment state

Frequency domain vibration of HGU in each state

d) Rubbing-misalignment state

The nth order NOFRFs is given by:

17
Gnjω=ω1++ωn=ωHnjω1,,jωni=1nUjωidσnωω1++ωn=ωi=1nUjωidσnω.

Substituting Eq. (16) and Eq. (17) into Eq. (15) lead to:

18
Yjω=n=1NGnjωUnjω.

When the HGU operates with fault, the change of frequency domain characteristics will react in the NOFRFs. But the NOFRFs always have a wide frequency band, which lead difficult to intelligent fault diagnosis methods. In this paper, we have proposed that use several state characteristics [23, 24] to quantify the feature of NOFRFs.

Frequency center of gravity:

19
FC=0+jωS(jω)djω0+S(jω)djω,

where jω, Sjω are frequency and power spectrum of NOFRFs.

Mean square frequency:

20
MSF=0+(jω)2S(jω)djω0+S(jω)djω.

Frequency domain variance:

21
VF=0+(jω-FC)2S(jω)djω0+S(jω)djω.

The above state characteristics of each order NOFRFs form the feature vector for failure analysis and diagnosis, as show in Table 3.

Table 3The feature vector based on NOFRFs

1st order NOFRFs
2st order NOFRFs
Nst order NOFRFs
Feature vector
FC
MSF
VF
FC
MSF
VF
FC
MSF
VF

3.2. SVM diagnosis system

Vapnik [16] introduced Support Vector Machine (SVM). It is a useful method for data classification. It is a learning system based on statistical learning theory which usually separate data into training and testing sets. And in each set, feature data have specific associated with labels. Each of these features of data can be looked upon as a dimension of a hyper-plane. SVM construct a hyper-plane that separates the hyper-space into two classes. SVM algorithm tries to achieve maximum separation between the classes. The classes are separated by a large margin minimizes the expected generalization error.

In this paper, there are four types of state for HGU. Based on the feature vector which proposed in Section 3.1, a SVM with binary tree architecture is built for multi-classification as show in Fig. 6. Firstly, the total samples are classified into two classes, one for the normal state and one for the other states. Secondly, the rest sample is classified into rubbing state and the others. And so on for the rest states. Finally, all four states can be classified by the SVM with binary tree architecture.

Fig. 6SVM diagnosis system with binary tree architecture

SVM diagnosis system with binary tree architecture

4. Experimental studies

4.1. Verification experiment of the proposed method

The HGU is a complex energy conversion system. According to the energy transfer direction, the vibration of lower guide bearing and upper guide bearing are set as input and output of the mode. In this paper, the vibration data in different states are acquired through finite element simulation which is proposed in Section 2. The vibration data of node 2 is regarded as the vibration data of upper bearing of the HGU, and the vibration data of node 9 is regarded as the vibration data of lower bearing of the HGU. For each of the four states, 20 groups of data are collected, and the sampling length is set as 4096 point/group. The sampling frequency is set as 20 Hz.

Firstly, the traditional method is used to analyze the different fault with the frequency domain spectrum of vibration data of node 2 as showed in Fig. 5. It can be seen that the rubbing state have a very small second harmonic. The misalignment state mainly performance on basic frequency and second harmonic. But for misalignment and misalignment-rubbing state, there are a lot of difficulties to distinguish between them.

To overcome the above problems, the method which is proposed in Section 3 is used for fault diagnosis. For most nonlinear system can be described as 3 order Volterra series [8], the NOFRFs model is set as 3 order system. The vibration data at node 9 and node 2 are used as input data and output data, respectively. And the NOFRFs models of different states are built. The spectrum of NOFRFs in different states is showed in Fig. 7. It can be seen that normal state mainly performance in first order NOFRFs, the second order and third order NOFRFs has few influence. In rubbing state, the first order NOFRFs focus on basic frequency and multi frequency, the second and third order NOFRFs is promiscuously. In misalignment state, the first and third order NOFRFs focus on the basic frequency, the second order NOFRFs mainly focus on the second harmonic. In misalignment-rubbing state, the rubbing changes the second order NOFRFs compare to misalignment state. So the proposed method has a strong effect to distinguish each state.

4.2. Experiment of anti-noise ability of the proposed method

Because of the complicated working condition of the HGU, the measurement data always contain noise jamming in practical engineering. So the verification test of anti-noise ability of the proposed method is necessary.

Fig. 7The NOFRFs of HGU in each state

The NOFRFs of HGU in each state

a) Normal state

The NOFRFs of HGU in each state

b) Rubbing state

The NOFRFs of HGU in each state

c) Misalignment state

The NOFRFs of HGU in each state

d) Rubbing-misalignment state

Fig. 8The NOFRFs of HGU in rubbing-misalignment state with different SNR noise

The NOFRFs of HGU in rubbing-misalignment state with different SNR noise

a) Without noise

The NOFRFs of HGU in rubbing-misalignment state with different SNR noise

b) With noise of 70 dB

The NOFRFs of HGU in rubbing-misalignment state with different SNR noise

c) With noise of 30 dB

The NOFRFs of HGU in rubbing-misalignment state with different SNR noise

d) With noise of 10 dB

In this section, the rubbing-misalignment state is used for the verification. In order to test the anti-noise ability, the white Gaussian noise of different signal-to-noise ratio are added to the measurement input-output data. The signal-to-noise ratio is set as 70 dB, 30 dB, 10 dB. The spectrum of NOFRFs in different SNR is showed in Fig. 8. It can be seen that the system’s character is performed clear even there exists noise interference.

The result indicates that the proposed method has a strong ability of restraining noise. Further pointed out that the proposed method is suitable for engineering practice.

4.3. Comparison experiment

In this section, the proposed method is compared to the wavelet analysis method and frequency domain spectrum analysis method to demonstrate the effectiveness of this method. For the proposed method, the feature is extracted with the method which proposed in Section 3.1. The 2-dimension distribution of NOFRFs features is shown in Fig. 9. And for the wavelet analysis method and frequency domain spectrum analysis method, the feature is extracted from the wavelet coefficient and frequency domain respectively.

Fig. 9The 2-dimension distribution of NOFRFs features

The 2-dimension distribution of NOFRFs features

In order to achieve the intelligent fault diagnosis and compare the proposed method to other fault diagnosis method, the SVM fault diagnosis system is proposed to verify the accuracy rate of the different methods. For the 20 groups of collected vibration data, 20 groups of feature are extracted for each method. Five groups of NOFRFs feature and traditional feature of each state are used as training data for the SVM fault diagnosis system respectively. And all groups of data of each state are used to test the fault diagnosis system. The results are showed in Table 4.

Table 4Fault diagnosis result of SVM system by different diagnosis methods

The proposed method
The wavelet analysis method
The frequency domain
spectrum analysis method
First SVM
0
0
1
Second SVM
0
1
0
Third SVM
0
1
6
Accuracy rate
100 %
97.5 %
91.25 %

The results indicate that the proposed method is useful for fault diagnosis in HGU.

5. Conclusion

Failure vibration data of HGU are difficult to acquire. And the traditional failure analysis and diagnosis methods ignore the characteristic of HGU which makes the failure analysis and diagnosis less effective. To overcome these problems, finite element method is used to simulate the vibration data of the HGU unit in different states, and a novel failure analysis and diagnosis method is proposed for HGU. The research work is conducted in three phases. Firstly, a finite element model is proposed according to the physical parameters of a real HGU, and the vibration data are acquired through numerical simulation. Secondly, a novel feature extracted method is proposed, and the feature in different state is extracted. Thirdly, a SVM system with binary tree architecture is built for fault classification. Finally, the proposed method is successfully employed on an experimental study. The results indicate that the proposed method can describe the failure characteristic better than the traditional time-frequency analysis method, and it is suitable for fault diagnosis of HGU.

Future research is to improve the computational efficiency and apply this method to practical engineering diagnostic cases.

References

  • Zhang Xiaoyuan, et al. Vibrant fault diagnosis for hydroelectric generator units with a new combination of rough sets and support vector machine. Expert Systems with Applications, Vol. 39, 2012, p. 2621-2628.
  • Dorji U., Ghomashchi R. Hydro turbine failure mechanisms: an overview. Engineering Failure Analysis, Vol. 44, 2014, p. 136-147.
  • Iliev H. Failure analysis of hydro-generator thrust bearing. Wear, Vol. 225, 1999, p. 913-917.
  • Momcilvic D., et al. Failure analysis of hydraulic turbine shaft. Engineering Failure Analysis, Vol. 20, 2012, p. 54-66.
  • Wang L., Cheung R. W., Ma Z., et al. Finite-element analysis of unbalanced magnetic pull in a large hydro-generator under practical operations. IEEE Transactions on Magnetics, Vol. 44, 2008, p. 1558-1561.
  • Li Shouju, Liu Yingxi Identification of vibration loads on hydro generator by using hybrid genetic algorithm. Acta Mechanica Sinica, Vol. 22, 2006, p. 603-610.
  • Nelson H. D., McVaugh J. M. The dynamics of rotor-bearing systems using finite elements. Journal of Engineering for Industry, Vol. 98, 1976, p. 593-600.
  • Tang Hao, et al. Fault diagnosis approach based on Volterra models. Mechanical Systems and Signal Processing, Vol. 24, 2010, p. 1099-1113.
  • Melgoza J. Jesús Rico, et al. An algebraic approach for identifying operating point dependent parameters of synchronous machines using orthogonal series expansions. IEEE Transactions on Energy Conversion, Vol. 16, 2001, p. 92-98.
  • Xiao Z., Jing X., Cheng L. Estimation of parametric convergence bounds for Volterra series expansion of nonlinear systems. Mechanical Systems and Signal Processing, Vol. 45, 2014, p. 28-48.
  • Xiao B., Gao C., Liu Z. G. Decoupling analysis on nonlinear system based on the modified generalized frequency response functions. Mechanical Systems and Signal Processing, Vol. 42, 2014, p. 283-299.
  • Cao Jianfu, et al. Fault diagnosis of complex system based on nonlinear frequency spectrum fusion. Measurement, Vol. 46, 2013, p. 125-131.
  • Lang Z. Q., Billings S. A. Energy transfer properties of non-linear systems in the frequency domain. International Journal of Control, Vol. 78, 2005, p. 345-362.
  • Lang Zi-Qiang, Billings S. A. Output frequency characteristics of nonlinear systems. International Journal of Control, Vol. 64, 1996, p. 1049-1067.
  • Lang Z. Q., Peng Z. K. A novel approach for nonlinearity detection in vibrating systems. Journal of Sound and Vibration, Vol. 314, 2008, p. 603-615.
  • Cortes Corinna, Vapnik Vladimir Support-vector networks. Machine Learning, Vol. 20, 1995, p. 273-297.
  • Lu Wenbo, et al. A gearbox fault diagnosis scheme based on near-field acoustic holography and spatial distribution features of sound field. Journal of Sound and Vibration, Vol. 332, 2013, p. 2593-2610.
  • Lu Wenbo, et al. A fault diagnosis scheme of rolling element bearing based on near-field acoustic holography and gray level co-occurrence matrix. Journal of Sound and Vibration, Vol. 331, 2012, p. 3663-3674.
  • Füllekrug U., Sinapius J. M. Identification of modal parameters, generalized and effective masses during base-driven tests. Aerospace Science and Technology, Vol. 2, 1998, p. 469-480.
  • Han Q., Zhang Z., Wen B. Periodic motions of a dual-disc rotor system with rub-impact at fixed limiter. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, Vol. 222, 2008, p. 1935-1946.
  • Al-Hussain K. M. Dynamic stability of two rigid rotors connected by a flexible coupling with angular misalignment. Journal of Sound and Vibration, Vol. 266, 2003, p. 217-234.
  • Cho J. R., Lee S. Y. Dynamic analysis of baffled fuel-storage tanks using the ALE finite element method. International Journal for Numerical Methods in Fluids, Vol. 41, 2003, p. 185-208.
  • Du R., Elbestawi Wu M. A. S. M. Automated monitoring of manufacturing processes, Part 1: Monitoring methods. Journal of Engineering for Industry, Vol. 117, 1995, p. 121-132.
  • Du R., Elbestawi M. A., Wu S. M. Automated monitoring of manufacturing processes, Part 2: Applications. Journal of Engineering for Industry, Vol. 117, 1995, p. 133-141.

About this article

Received
06 January 2016
Accepted
28 June 2016
Published
30 September 2016
SUBJECTS
Fault diagnosis based on vibration signal analysis
Keywords
hydraulic-turbine generator unit
fault diagnosis
finite element
nonlinear output frequency response functions
Acknowledgements

This work is supported by National Science Foundation for Distinguished Young Scholars of China (61203056). Platform Construction Project of Technology Innovation of Huaian China (HAP201432). We greatly thank our university for their great assistance to our study and helpful suggestions to the presentation of this paper.