Dynamic Parameter Identification of Hydrodynamic Bearing-Rotor System

A newmethod calledmodal parameter genetic time domain identification was employed to study the characteristics of the bearingrotor system.Amultifrequency signal decomposition technology to identify themain components of themeasured signal and reject the image mode produced by noise has been used. The firstand second-order 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.Themethod has been effective albeit with few testing points and operation cases.The derivation of oil-film 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 firstand second-order 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.


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][4][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.

Shock and Vibration
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 oil-film 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 oil-film coefficients [6][7][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 oil-film coefficients and presented the characteristics under various operational conditions.They had pointed out that linear oil-film dynamic coefficients are sensitive to excitation amplitude.Of course, their linear oil-film 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 fluid-film-bearing parameters from larger journal orbital motion.Nonlinear coefficients are defined by third-order 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 oil-film 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 oil-film 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 rotor-bearing-seal test rigs, and theoretical procedures to extract the relevant bearing and seal dynamics parameters and estimate-associated parameter uncertainties.They had concluded that the synchronous unbalance response, which can easily be obtained from the run-down/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 second-order natural frequencies matched quite well those obtained from modal parameter identification.The conclusion is presented in Section 4.

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 as where   is damping,   is vibration frequency without damping, while   is vibration frequency with damping, and The coordinates of transformation of displacement vector {} are given by Substituting (3) into (1), Equation ( 4) could be solved taking the initial conditions into account: where   and   are constants depending on the initial conditions: 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.

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 ∑  =1    (  +  ) 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 from where  = 0, 1, . . .,  − 1;   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  2    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 V +1 , V +2 , . . ., V  are noise Eigenvectors, which have the same Eigenvalues  2  ; V 1 , V 2 , . . ., V  are signal Eigenvector whose Eigenvalues are  1 +  2  ,  2 +  2  , . . .,   +  2  .It could be seen that white noise has an influence on sine signal Eigenvalues.Because the signal vectors   ( = 1, 2, . . ., ) and noise Eigenvectors V +1 , V +2 , . . ., V  are orthogonal [20], their linear combination is also orthogonal: Let   = 1 ( = +1, +2, . . ., ).Now the function   () could be defined similarly to the power spectrum function: The function-independent 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 V  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  = 1, 2, 3, . .., 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.

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 run-down 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 two-order natural frequencies are 3.9 Hz and 9.7 Hz, and the first two-order damping ratios are 2.23% and 1.56%.

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 oil-film force: where {} is the node displacement vector; {ℎ} is the mass unbalanced force vector; {} is the bearing oil-film force matrix.
The response of the rotor system could be expressed by the superposition of multiple sine and cosine components: where   ( = 1, 2, . . ., ) 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 as Arbitrary journal bearing node displacement is given by where  = 1, 2, . . .,  denotes the journal bearing number, the total number of journal bearings being .
The derivate { Q  ()} could be expressed as The corresponding bearing oil-film force is given by Let The oil-film force could be expressed as where Then, where Let Now {   } and {   } could be simplified as Substituting ( 27) and ( 21) into (20) and decomposing cos(  ) and sin(  ) components, Let Finally, the following equation could be obtained: The arbitrary node displacement of the shaft system could be expressed by the unbalanced force and oil-film force.The structural parameters of the shaft system are known; that is, matrix  is known.The oil-film force is a function of journal bearing stiffness and damping coefficients.Hence, the genetic algorithm could be utilized to identify the oil-film 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 second-order natural frequencies are 4.23 Hz and 9.37 Hz.The calculated first-and second-order 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

Conclusions
(1) 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.
(2) The oil-film 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 oil-film 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 second-order 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.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.

Figure 1 :Figure 2 :
Figure 1: Vibration response of the shaft system.

Table 1 :
Parameters of the generator set.

Table 2 :
Identification results of dynamic characteristic coefficients.