Research Article  Open Access
Zhiqiang Song, Yunhe Liu, "Dynamic Parameter Identification of Hydrodynamic BearingRotor System", Shock and Vibration, vol. 2015, Article ID 959568, 7 pages, 2015. https://doi.org/10.1155/2015/959568
Dynamic Parameter Identification of Hydrodynamic BearingRotor System
Abstract
A new method called modal parameter genetic time domain identification was employed to study the characteristics of the bearingrotor system. A multifrequency signal decomposition technology to identify the main components of the measured signal and reject the image mode produced by noise has been used. The first and secondorder natural frequency and damping ratios of the shaft system are identified. Furthermore, because of the deficiency of the traditional least square method, a new genetic identification method to identify the bearing dynamic characteristic parameters has been proposed. The method has been effective albeit with few testing points and operation cases. The derivation of oilfilm dynamic coefficients could also provide a basis for shaft system natural vibration characteristic and vibration response analysis. Using the identified dynamic coefficients as the supporting condition, the shaft system modal characteristics were studied. The calculated first and secondorder natural frequencies match quite well those obtained from the modal parameter identification. It was proved that the modal parameter and physical parameter identification methods utilized in this paper are reasonable.
1. Introduction
Recently, the finite element methods have been widely used in the study on vibration characteristics of rotating machines such as the water turbine generator set shaft system [1, 2]. However, due to the uncertainty of the boundary condition or dynamic parameters, the simulation results were sometimes far from practical. The system identification theory deals with the vibration problem depending on the theory analysis and prototype testing. It could obtain the boundary condition or dynamic physical parameters according to the measured structure response. The identification could provide the rationality of boundary condition or dynamic parameters for further response characteristic analysis of the rotating machine.
To date, system identification can generally be classified into two categories: (a) identification of modal parameters of the system and (b) identification of physical parameters of the system. Modal parameters reflect the natural vibration characteristics of the system better than physical parameters. The number of modal parameters is also fewer than that of physical parameters. Hence, the modal parameter identification is the basis of system identification. Over the decades, modal parameter identification has become well accepted and widely used in the study of vibration stability of rotating machines. However, the traditional method of frequency domain identification must require a significant number of tests to obtain the frequency response functions. Obtaining the frequency response function may be timeconsuming and costly; the identification would not be carried out without the function. In the recent past, it has been observed that time domain identification has played a significant role in the rapid resolution for modal parameters [3–5].
Based on the linear structure vibration mode superposition theory, the time domain identification method was proposed to obtain the modal parameters from the measured structure response signals. This method could identify the modal parameters of large structures that work only under environment excitation. It did not need system input information. It had no influence on the normal operation of the structure or equipment. Hence, this method is very suitable for modal identification of continuously operating equipment such as the water turbine generator set shaft system.
Compared with modal parameter identification, the physical parameters were much more direct and visual. The physical parameters obtained by identification are the basis of structural characteristic analysis and vibration response resolutions. The physical parameters that need to be identified for the generator set shaft system are just the oilfilm dynamic characteristic coefficients of the journal bearings. Rotating machines are supported by bearings, which play a vital role in determining the behavior of rotating systems under the action of dynamic loads. One of the most important factors governing the vibration characteristics of rotating machinery is the parameters of bearing dynamics. Although the importance of rotor support dynamic stiffness is generally well recognized by the design engineer, it is often the case that theoretical models available for predicting it are insufficiently accurate or are accurate only in very specific cases. Moreover, the stiffness and damping characteristic coefficients are greatly dependent on many physical and mechanical parameters such as lubricant temperature, bearing clearance and load, journal speed, and machine misalignment in the system, and these are difficult to obtain accurately under test conditions. Uncertainties about machine parameters can render the results obtained with the best theoretical methods inaccurate in the study of the behavior of fluidfilm journal bearings. Owing to this, it may be very useful to determine the unbalance and bearing dynamic stiffness by means of identification methods based on experimental data and machine models.
Many researchers have worked on the identification of oilfilm coefficients [6–8]. Yao et al. [9] have utilized the rotating machinery analysis to obtain the critical speed and oilfilm stiffness. Because the critical speed of water electric units is higher than the rated speed or working speed, this method is not suitable for water electric units. Lees and Friswell [10] have presented a method to evaluate the state of unbalance and bearing coefficients of a rotating machine using the measured pedestal vibration. Subsequently, Edwards et al. [11] have presented experimental verification of the method to evaluate the unbalance and bearing coefficients of a rotating machine. Tiwari and Chakravarthy [12] have proposed an identification algorithm for simultaneous estimation of residual unbalances and bearing dynamic parameters. The algorithm uses standard condensation to overcome the difficulty with too much response freedom that needs to be measured. The identification parameters matched quite well the parameters assumed for the simulation of responses. Zhao et al. [13] have proposed an experimental method to identify linear oilfilm coefficients and presented the characteristics under various operational conditions. They had pointed out that linear oilfilm dynamic coefficients are sensitive to excitation amplitude. Of course, their linear oilfilm model would be invalid if the condition of small perturbation is not satisfied. Meruane and Pascual [14] have proposed a framework for numerical identification of nonlinear fluidfilmbearing parameters from larger journal orbital motion. Nonlinear coefficients are defined by thirdorder Taylor expansion of bearing reaction forces and are evaluated through least mean square in time domain technique. Lou and Bi [15, 16] have also utilized the least square method to obtain oilfilm coefficients. This method needs data from different running cases to satisfy the request about the equation number. However, if the operation cases change greatly, the oilfilm coefficients will vary greatly, which may lead to the premise that the identification formula is not satisfied any more. Recently, Tiwari et al. [17, 18] have provided a review of the identification procedures applied to bearing and seal models, the existing experimental techniques for acquiring the measurement data from the rotorbearingseal test rigs, and theoretical procedures to extract the relevant bearing and seal dynamics parameters and estimateassociated parameter uncertainties. They had concluded that the synchronous unbalance response, which can easily be obtained from the rundown/up of large turbo machines, should be exploited more for the identification of bearing dynamics parameters along with the estimation of residual unbalance.
The organization of the rest of the paper is as follows: theory of the modal parameter identification method is described in Section 2; a newly developed multifrequency signal decomposition technology has been used to identify the main components of measured signal and reject the image mode produced by noise. Section 3 presents a new genetic identification method for the journal bearing dynamic characteristic parameters; using the dynamic coefficients obtained, the shaft system modal characteristics can be then studied using finite element method; the calculated first and secondorder natural frequencies matched quite well those obtained from modal parameter identification. The conclusion is presented in Section 4.
2. Modal Parameter Genetic Time Domain Identification
2.1. Algorithm for Genetic Time Domain Identification
For a system with freedoms, the equation of free motion is given as Let , where is called the free response amplitude matrix. Hence, the Eigenvalue could be expressed aswhere is damping, is vibration frequency without damping, while is vibration frequency with damping, and .
The coordinates of transformation of displacement vector are given bySubstituting (3) into (1),Equation (4) could be solved taking the initial conditions into account:where and are constants depending on the initial conditions:Substituting (5) into (3),where and .
The precision request could be satisfied well as long as the modal displacement responses of first few orders are taken.
The modal parameter genetic identification method is based on the system identification and genetic optimization theory [19]. This method utilizes the free vibration response measurement signals of the shaft system. The chromosomes are used to express the modal parameters such as natural frequencies and damping ratios. The parameters that made the objective function optimal are just the structural modal parameters.
The objective function could be expressed by an equation that includes the measured response and calculated response: where is the measured response.
is the calculated response ( is the modal order number).
2.2. Multifrequency Signal Process
Vibration signals are always composed of many frequency components. The noise would produce the image mode with poor identification accuracy. Using the signal decomposition technology, a new multifrequency identification algorithm to identify the main components of a multifrequency signal and reject the image mode produced by noise has been proposed.
A random process could be expressed as follows: where is the signal composed by sine waves, the signal determined order is to make sure of the number of sine waves , and is additive white noise.
The autocorrelation function of signal could be calculated fromwhere ; is the square variance of the random signal. Now, define the signal vectors where is the frequency of sine signals.
The autocorrelation matrix of signal is expressed by the summation of signal autocorrelation matrix and noise autocorrelation matrix: where is the noise autocorrelation matrix and is an dimensional unit matrix.
could be expressed by its Eigenvalues and Eigenvectors:
Generally, the rank of signal autocorrelation matrix (i.e., ) is less than that of its dimension number . Hence, matrix has Eigenvalues that are nonzero:
also could be expressed by its Eigenvalues and Eigenvectors. Because all its Eigenvalues are 1, any vectors could be regarded as its Eigenvectors : Substituting (14) and (15) into (12), where are noise Eigenvectors, which have the same Eigenvalues ; are signal Eigenvector whose Eigenvalues are . It could be seen that white noise has an influence on sine signal Eigenvalues. Because the signal vectors and noise Eigenvectors are orthogonal [20], their linear combination is also orthogonal:Let . SoLet . Now the function could be defined similarly to the power spectrum function:
The functionindependent variable is just a sine frequency of the signal when the function has a wave peak. It seems that the identification frequency could make the function tend to infinity. In fact, the Eigenvectors would be obtained from the approximate autocorrelation matrix. Because of the estimate error, the function only has a finite peak.
Now is the order number of the signal that needs to be determined. Assuming , find out whether the function has a peak in the possible frequency domain. If the function has a peak, the corresponding is just the sine wave number in the signal. Thus, the vibration modal order number of the structure could be obtained. For the hydroelectric unit shaft system, the natural frequency region is about 0~150 Hz and may be 1~5. Hence, the computation cost is not very high.
2.3. Modal Parameter Identification of the Shaft System
Figure 1 shows the free vibration response signal of a factual unit shaft at the turbine guide bearing location at rundown operating condition. The new multifrequency identification algorithm was used to identify the main components of a multifrequency signal and reject the image mode produced by noise. The result is shown in Figure 2. It can be seen that the signal includes two orders of the vibration modals. The modal parameters were identified using the genetic time domain identification method. The identification results show that the first twoorder natural frequencies are 3.9 Hz and 9.7 Hz, and the first twoorder damping ratios are 2.23% and 1.56%.
3. Bearing Dynamics Parameter Identification
The identification of bearing dynamics parameters is helpful to understand the shaft system support conditions at normal operation cases. This may provide the basis for the study of shaft system modal characteristics and vibration responses.
The dynamic differential equation could be expressed as follows considering the mass unbalanced force and bearing oilfilm force: where is the node displacement vector; is the mass unbalanced force vector; is the bearing oilfilm force matrix.
The response of the rotor system could be expressed by the superposition of multiple sine and cosine components:where is the ratio of a sine or cosine component to the rated rotating speed, which could be only a plus integer or rational fraction; is the total number of humorous wave components. Now, the mass unbalanced force could be written asArbitrary journal bearing node displacement is given by where denotes the journal bearing number, the total number of journal bearings being .
The derivate could be expressed asThe corresponding bearing oilfilm force is given by LetThe oilfilm force could be expressed aswhereThen,whereLetNow and could be simplified asSubstituting (27) and (21) into (20) and decomposing and components,LetFinally, the following equation could be obtained:
The arbitrary node displacement of the shaft system could be expressed by the unbalanced force and oilfilm force. The structural parameters of the shaft system are known; that is, matrix is known. The oilfilm force is a function of journal bearing stiffness and damping coefficients. Hence, the genetic algorithm could be utilized to identify the oilfilm dynamic characteristic parameters.
The genetic algorithm evaluates the chromosomes only according to the objective function. The objective of optimization is to find a group of parameters that make the square variances between the measured response and calculated response minimum:where is the number of the harmonics.
There are 28 parameters that need to be identified, including the journal bearing stiffness and damping coefficients and residual unbalances. The individual chromosome string could be very long and the colony size very large because of too many parameters. Few methods are available to reduce the huge calculation costs. In one method, each parameter’s search range could be set smaller according to experience and is so classified to ensure that the parameters obtained are available.
The main structural parameters for the practical water turbine generator set are given in Table 1. The identification results are given in Table 2.


The journal guide bearing stiffness coefficients are identified and the supporting boundary condition of the generator set shaft system is found. Now, the modal characteristic of shaft system could be calculated by finite element modal analysis. The simulation results show that the first and secondorder natural frequencies are 4.23 Hz and 9.37 Hz. The calculated first and secondorder natural frequencies matched quite well those obtained from modal parameter identification. This proves that the modal parameter and physical parameter identification methods utilized in this paper are reasonable. Further, it can be confirmed that the unit would be steady because the firstorder natural frequency is far from the rated frequency and runaway frequency.
4. Conclusions
Modal parameter genetic time domain identification method proposed in this paper is suitable for larger structures or equipment that works only under environment excitation. The method has no influence on the normal operation of structures or equipment. The method makes the modal parameter identification of the running unit shaft system convenient. A new multifrequency signal decomposition technology to identify the main components of the measured signal and reject the image mode produced by noise has been used. This identification method offers a new study approach on the modal vibration characteristics of the rotating machine shaft system.
The oilfilm dynamic parameter genetic identification method proposed in this paper avoids the deficiency of traditional least square method, such as the need for many case responses, which is often unrealistic because oilfilm dynamic parameters change greatly in different cases. Using the dynamic coefficients identified, the shaft system modal characteristics could then be studied by finite element method. The calculated first and secondorder natural frequencies matched quite well those obtained from modal parameter identification. It was proved that the modal parameter and physical parameter identification methods proposed in this paper are reasonable. The result of identification, that is, stiffness coefficients, could also provide a basis for the study of shaft system vibration response. Otherwise, because of the complex of the vibration sources loads, the work about the feed calculation of shaft system vibration response and contrast with the prototype test dates is being processed.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research was supported by the Natural Science Foundation of China (51479165 and 51179154) and Program 2013KCT15 for Shaanxi Provincial Key Innovative Research Team. Sincere gratitude is extended to the editor and anonymous reviewers for their professional comments and corrections, which greatly improved the presentation of the paper.
References
 Z. Wang, J. Yu, Y. Fang, X. Wen, J. Cao, and Q. Shi, “Characteristic analysis of rotor dynamics of large hydraulic generating unit,” Journal of Hydroelectric Engineering, vol. 24, no. 4, pp. 62–66, 2005. View at: Google Scholar
 X. An, J. Zhou, L. Liu, and X. Xiang, “Lateral vibration characteristic analysis of the hydro generator set,” Lubrication Engineering, vol. 33, no. 2, pp. 40–44, 2008. View at: Google Scholar
 X. Xu, H. Hua, and Z. Chen, “Review of modal identification method based on ambient excitation,” Journal of Vibration and Shock, vol. 21, no. 3, pp. 1–5, 2002. View at: Google Scholar
 S. Li and J. Lian, “Genetic algorithm for hydraulic power house modal parameters identification,” Journal of Tianjin University, vol. 42, no. 1, pp. 11–16, 2009. View at: Google Scholar
 H.D. Zhang, Y. Zhou, and J.J. Lian, “Method based on ARMA model for modal parameter identification of a power house,” Journal of Vibration and Shock, vol. 26, no. 5, pp. 115–118, 2007. View at: Google Scholar
 P. Pennacchi, N. Bachschmid, A. Vania, G. A. Zanetta, and L. Gregori, “Use of modal representation for the supporting structure in modelbased fault identification of large rotating machinery: part 1—theoretical remarks,” Mechanical Systems and Signal Processing, vol. 20, no. 3, pp. 662–681, 2006. View at: Publisher Site  Google Scholar
 P. Mohanty and D. J. Rixen, “A modified Ibrahim time domain algorithm for operational modal analysis including harmonic excitation,” Journal of Sound and Vibration, vol. 275, no. 12, pp. 375–390, 2004. View at: Publisher Site  Google Scholar
 Y. Guo, Y. Chen, J. Wang, and X. Yuan, “Principle and measure point research of bearing oilfilm force identification,” Lubrication Engineering, vol. 32, no. 8, pp. 39–41, 2007. View at: Google Scholar
 D. Yao, W. Huang, and J. Zhou, “Identification of the journal bearing stiffness coefficients,” Dynamic Engineering, vol. 25, no. 4, pp. 483–486, 2005. View at: Google Scholar
 A. W. Lees and M. I. Friswell, “The evaluation of rotor imbalance in flexibly mounted machines,” Journal of Sound and Vibration, vol. 208, no. 5, pp. 671–683, 1997. View at: Publisher Site  Google Scholar
 S. Edwards, A. W. Lees, and M. I. Friswell, “Experimental identification of excitation and support parameters of a flexible rotorbearingsfoundation system from a single rundown,” Journal of Sound and Vibration, vol. 232, no. 5, pp. 963–992, 2000. View at: Publisher Site  Google Scholar
 R. Tiwari and V. Chakravarthy, “Simultaneous identification of residual unbalances and bearing dynamic parameters from impulse responses of rotorbearing systems,” Mechanical Systems and Signal Processing, vol. 20, no. 7, pp. 1590–1614, 2006. View at: Publisher Site  Google Scholar
 S. X. Zhao, H. Zhou, G. Meng, and J. Zhu, “Experimental identification of linear oilfilm coefficients using leastmeansquare method in time domain,” Journal of Sound and Vibration, vol. 287, no. 45, pp. 809–825, 2005. View at: Publisher Site  Google Scholar
 V. Meruane and R. Pascual, “Identification of nonlinear dynamic coefficients in plain journal bearings,” Tribology International, vol. 41, no. 8, pp. 743–754, 2008. View at: Publisher Site  Google Scholar
 X. Lou, Unbalance Identification Method Study for the Operating Rotor, Zhejiang University, Zhejiang, China, 1998.
 S. Bi, W. Jing, and N. Cui, “Recursive damping least square estimation of dynamic characteristic for journal bearings,” Journal of Vibration Engineering, vol. 11, no. 4, pp. 410–415, 1998. View at: Google Scholar
 R. Tiwari, A. W. Lees, and M. I. Friswell, “Identification of dynamic bearing parameters: a review,” Shock & Vibration Digest, vol. 36, no. 2, pp. 99–124, 2004. View at: Publisher Site  Google Scholar
 R. Tiwari, S. Manikandan, and S. K. Dwivedy, “A review of the experimental estimation of the rotor dynamic parameters of seals,” Shock & Vibration Digest, vol. 37, no. 4, pp. 261–284, 2005. View at: Publisher Site  Google Scholar
 X. Wang and L. Cao, Theory, Application and Software Implementation of Genetic algorithm, Xi'an Jiaotong University Press, Xi'an, China, 2002.
 G. Lu, X. Peng, and L. Zhang, Random Signal Process, Xi Dian University Press, Xi'an, China, 2002.
Copyright
Copyright © 2015 Zhiqiang Song and Yunhe Liu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.