Establishment of Approximate Analytical Model of Oil Film Force for Finite Length Tilting Pad Journal Bearings

Abstract

Tilting pad bearings offer unique dynamic stability enabling successful deployment of high-speed rotating machinery. The model of dynamic stiffness, damping, and added mass coefficients is often used for rotordynamic analyses, and this method does not suffice to describe the dynamic behaviour due to the nonlinear effects of oil film force under larger shaft vibration or vertical rotor conditions. The objective of this paper is to present a nonlinear oil force model for finite length tilting pad journal bearings. An approximate analytic oil film force model was established by analysing the dynamic characteristic of oil film of a single pad journal bearing using variable separation method under the dynamic π oil film boundary condition. And an oil film force model of a four-tilting-pad journal bearing was established by using the pad assembly technique and considering pad tilting angle. The validity of the model established was proved by analyzing the distribution of oil film pressure and the locus of journal centre for tilting pad journal bearings and by comparing the model established in this paper with the model established using finite difference method.

Yongliang Wang 1 and Yu Gao 1 and Ying Cui 1 and Zhansheng Liu 2

Academic Editor:Paolo Pennacchi

1, Marine Engineering College, Dalian Maritime University, Dalian, Liaoning 116026, China
2, School of Energy Science and Engineering, Harbin Institute of Technology, Harbin, Heilongjiang 150001, China

Received 26 January 2015; Accepted 29 April 2015; 20 May 2015

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.

1. Introduction

Tilting pad journal bearings are widely used in large steam turbine-generator, heavy duty gas turbines, vertically oriented machines, and so forth, due to their excellent dynamic stability compared with other oil film bearing designs [1-3]. Dynamic performance of tilting pad bearing is one of the important parameters affecting the dynamic characteristics of high-speed rotor supported by tilting pad journal bearings.

In system with horizontal rotors, when the journal dynamic displacements about the equilibrium position are far smaller than the bearing clearance, the bearing reaction forces can be expressed as a stiffness, damping, and added mass coefficients model [4, 5], but lager amplitude journal motions mean the linear model has limits to apply. And in vertical rotors with tilting pad bearing, there is no equilibrium position; for example, in vertical hydropower units and pumps, the position of rotor determined by the generator's properties, the flow conditions in the turbine, the mass unbalance, and the linear oil film force models cannot perform simulations of a vertical rotor's behavior [2]. It is therefore of great significance to calculate the nonlinear oil film force of a tilting pad journal bearing for rotordynamic analysis under severe unbalance conditions or large amplitude journal motions [6].

"Lund's pad assembly method" presented by Lund in 1964 is a commonly known procedure for determining the dynamic characteristics of tilting pad journal bearing [7]; many design and analysis works of tilting pad journal bearing have been done in the last 60 years based on "Lund's pad assembly method" [8, 9]. And the nonlinear oil film force of a tilting pad bearing can be calculated at present using numerical method [10, 11], data base method [12, 13], or nonlinear oil film force formula, and so forth. The nonlinear formula is more efficient and convenient compared with the other methods. Okabe and Cavalca [14, 15] presented an analytical model of a tilting pad journal bearing with the short bearing assumption and reported that the analytical model was faster than numerical method. The short bearing model has limits in performing dynamic characteristics of tilting pad journal bearing for length-diameter ratio larger than 0.5.

The aim of the present work and its main contribution is related to the establishment of an analytical oil film force model of finite length tilting pad journal bearings. An approximate analytic model of a single pad is established through variable separation method under the dynamic [figure omitted; refer to PDF] oil film boundary condition. And an approximate analytical model of a 4-pad finite length tilting pad journal bearings is presented by using the pad assembly method and by considering pad tilting angle.

2. Equation for Pad Tilting Angle

During the calculation of oil film force for tilting pad bearing under given disturbance displacements and velocities of journal, if tilting angle can be quickly obtained, the distribution of oil film thickness can be obtained, and consequently the dynamic characteristics of oil film force of the tilting pad are obtained efficiently.

As shown in Figure 1(a), a pad is installed on supporting block, the journal is directly above the pad, and the gap between the pad and the journal is full of lubricant. Assume that the journal is motionless at the beginning and that [figure omitted; refer to PDF] as the center of pad arc, [figure omitted; refer to PDF] as the center of journal, and [figure omitted; refer to PDF] as the pivot of pad are on the same straight line. The distance between the center of journal and the center of pad, [figure omitted; refer to PDF] , is defined as original eccentricity [figure omitted; refer to PDF] , and the original eccentricity ratio is defined as [figure omitted; refer to PDF] , where [figure omitted; refer to PDF] is bearing radius clearance.

Figure 1: Schematic diagram of pad-swing angle. (a) With journal motionless. (b) With journal rotating at velocity [figure omitted; refer to PDF] .

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

When the journal starts to rotate clockwise at angular velocity [figure omitted; refer to PDF] , assume [figure omitted; refer to PDF] is motionless while the journal is rotating. According to fluid lubrication theory, an oil film is formed and pressure is generated under the shear stress of fluid film, as the lubricant is moved between the journal and the pad. After the pad and the neck journal get stable, the tilting pad will rotate through angle [figure omitted; refer to PDF] about the pivot point. As shown in Figure 1(b), angle [figure omitted; refer to PDF] is the pad tilting angle.

The relationship among pad tilting angle [figure omitted; refer to PDF] , original eccentricity [figure omitted; refer to PDF] , and tilting bearings clearance ratio [figure omitted; refer to PDF] , with a wrap angle of 80° and a pad pivot offset of 50%, can be expressed as shown below as reported in [16]: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] .

3. Analysis of Oil Film Characteristics of a Single Tilting Pad

3.1. Oil Film Force Model for Cylindrical Journal Bearing

Figure 2 illustrates the structure of oil film in cylindrical journal bearing, where [figure omitted; refer to PDF] is the radius, [figure omitted; refer to PDF] is bearing center, and [figure omitted; refer to PDF] is journal center. The journal rotates at velocity [figure omitted; refer to PDF] in the bearing and [figure omitted; refer to PDF] is the eccentricity of journal in bearing. The [figure omitted; refer to PDF] Cartesian coordinate is established as shown in Figure 2. The journal displacement and velocity are setting as [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] , respectively. [figure omitted; refer to PDF] is circumferential coordinate, while the start is [figure omitted; refer to PDF] -axis positive and the direction is anticlockwise and [figure omitted; refer to PDF] is the oil film thickness.

Figure 2: Structure of oil film in cylindrical journal bearing.

[figure omitted; refer to PDF]

According to fluid lubrication theory, the nondimensional Reynolds equation of oil film pressure [figure omitted; refer to PDF] in a cylindrical journal bearings can be expressed as [figure omitted; refer to PDF]

The corresponding nondimensional boundary conditions can be expressed as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the nondimensional oil film pressure; [figure omitted; refer to PDF] is radius clearance of bearing; [figure omitted; refer to PDF] is the viscosity of lubrication oil; [figure omitted; refer to PDF] is the diameter of bearing; [figure omitted; refer to PDF] is the length of bearing; [figure omitted; refer to PDF] is the length-diameter ratio; [figure omitted; refer to PDF] is the axial coordinate; [figure omitted; refer to PDF] is the nondimensional length-diameter ratio; and [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] are nondimensional displacements and velocities of disturbance, respectively. [figure omitted; refer to PDF] is the thickness of the nondimensional oil film; [figure omitted; refer to PDF] is the nondimensional time; and [figure omitted; refer to PDF] is the time.

The inlet and outlet circumferential boundary conditions for oil film under the dynamic [figure omitted; refer to PDF] boundary conditions can be expressed as shown below: [figure omitted; refer to PDF]

The approximate analytic model of nondimensional oil film pressure can be established as shown below using variable separation method under the dynamic [figure omitted; refer to PDF] boundary conditions as reported in [17]: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] Consider [figure omitted; refer to PDF] , [figure omitted; refer to PDF] . The expression of [figure omitted; refer to PDF] is shown in the Appendix.

The expression of oil film force can be established by integrating the oil film force along the oil film domain as shown below: [figure omitted; refer to PDF]

3.2. Model of Single Tilting Pad Oil Film Force in a Local Coordinate System

As shown in Figure 3, the disturbance displacement and velocity of a neck journal in a local coordinate system are ( [figure omitted; refer to PDF] ) and ( [figure omitted; refer to PDF] ).

Figure 3: Tilting pad oil film in local coordinates system.

[figure omitted; refer to PDF]

The expression of oil film pressure in a local coordinate system can be established by analyzing the oil film characteristics of a tilting pad in a local coordinate system and by determining the circumferential boundary conditions according to the pad parameters and the dynamic [figure omitted; refer to PDF] oil film boundary condition for a cylindrical bearing. The dynamic [figure omitted; refer to PDF] oil film boundary condition in a local coordinate system can be expressed as [figure omitted; refer to PDF] and from formula (4) [figure omitted; refer to PDF]

For a tilting pad, the actual oil film is the area [ [figure omitted; refer to PDF] projected on the pad. The tilting pad bearing can be divided into five areas in the circumferential direction according to its geometric structure. When [figure omitted; refer to PDF] is in these areas, the expressions of oil film boundaries for both upper and lower pads are as follows:

: If [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] ;

: if [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] ;

: if [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] ;

: if [figure omitted; refer to PDF] , the oil film pressure is zero;

: if [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] .

The expression of oil film pressure in a tilting pad under local coordinate is as follows: [figure omitted; refer to PDF]

And the dimensionless oil film force of a tilting pad under local coordinate is as follows: [figure omitted; refer to PDF]

4. Relations of Oil Film Force in Different Coordinates Systems

In Figure 4(a), the center of [figure omitted; refer to PDF] coordinate system is on the center of bearings [figure omitted; refer to PDF] , [figure omitted; refer to PDF] is the center of journal, the radius of journal is [figure omitted; refer to PDF] , the radius clearance between the pad and journal is [figure omitted; refer to PDF] , and the arc radius of the tilting pad is [figure omitted; refer to PDF] .

Figure 4: Tilting pad bearing and rotor journal.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

As shown in Figure 4(b), the pivot point of the pad is determined by the pivot angle [figure omitted; refer to PDF] , if [figure omitted; refer to PDF] is the center of pad arc, and the preload coefficient of the bearings is [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] . If the pad rotates through angle [figure omitted; refer to PDF] , the pad arc center moves from [figure omitted; refer to PDF] to [figure omitted; refer to PDF] ; the journal center [figure omitted; refer to PDF] is on the line connecting the pivot point and arc center [figure omitted; refer to PDF] , as shown in Figure 4(c); [figure omitted; refer to PDF] is the original eccentricity [figure omitted; refer to PDF] mentioned above. When the journal rotates at velocity [figure omitted; refer to PDF] , the pad will consequentially tilt. As shown in Figure 4(d), the pad tilting angle is [figure omitted; refer to PDF] and the pad arc center moves from [figure omitted; refer to PDF] to [figure omitted; refer to PDF] .

By setting the coordinates of NO, [figure omitted; refer to PDF] pad pivot in a global coordinate system as [figure omitted; refer to PDF] , their specific expression is as follows: [figure omitted; refer to PDF]

Assume dimensionless parameter [figure omitted; refer to PDF] : [figure omitted; refer to PDF]

Angle [figure omitted; refer to PDF] can be calculated using the triangular relations. As shown in Figure 4(c), according to geometric relationship [figure omitted; refer to PDF] , [figure omitted; refer to PDF]

By defining counterclockwise is positive [figure omitted; refer to PDF]

The original eccentricity can be expressed as [figure omitted; refer to PDF]

After pad original eccentricity [figure omitted; refer to PDF] is obtained, pad tilting angle [figure omitted; refer to PDF] can be calculated using (1) with a pad clearance ratio [figure omitted; refer to PDF] .

The corresponding local coordinate system is defined in Figure 5 which is an enlarged view of Figure 4(d), where coordinate [figure omitted; refer to PDF] is a global coordinate system with its origin at the bearing center [figure omitted; refer to PDF] and [figure omitted; refer to PDF] as its basis; [figure omitted; refer to PDF] is a local coordinate system fixed on the pad fulcrum and bearing center [figure omitted; refer to PDF] , with [figure omitted; refer to PDF] as its basis; and [figure omitted; refer to PDF] is a moving coordinate system with its original at the pad center [figure omitted; refer to PDF] and [figure omitted; refer to PDF] as its basis.

Figure 5: Relationship of coordinates systems.

[figure omitted; refer to PDF]

Using the algorithm for space vector, the following expression can be established in accordance with the geometric relationship shown in Figure 5: [figure omitted; refer to PDF] where [figure omitted; refer to PDF]

Similarly, the transitions of velocity are [figure omitted; refer to PDF]

And the relation of dimensionless oil film force under global coordinate system and local coordinate system is as follows: [figure omitted; refer to PDF]

The dimensionless forms of (17), (19), and (20) are as follows: [figure omitted; refer to PDF]

When disturbance displacement [figure omitted; refer to PDF] and velocity [figure omitted; refer to PDF] of a journal in bearings are known, disturbance displacement [figure omitted; refer to PDF] and velocity [figure omitted; refer to PDF] of the journal in a local coordinate system can be calculated using (17) and (19). Nondimensional oil film force [figure omitted; refer to PDF] is calculated in the local coordinate system of the pad. Finally, it can be transformed to using (23) into the oil film force [figure omitted; refer to PDF] in a global coordinates system.

5. Approximate Analysis Model of Oil Film Force for Tilting Pad Bearing

The nonlinear oil film force of a four-tilting-pad bearing can be calculated through the process shown in Figure 6 when the parameters of bearing, disturbance displacement, and disturbance velocity of a journal in the bearing are identified.

Figure 6: Process of solving oil film force of tilting pad bearing.

[figure omitted; refer to PDF]

6. Analysis of Oil Film Characteristics of Tilting Pad Bearing

The dimensionless pressure distributions of tilting pad journal bearing are got by evaluating (10) for different length-diameter ratio ( [figure omitted; refer to PDF] ), shown in Figure 7, when [figure omitted; refer to PDF] [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , attitude angle of 45°, [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] .

Figure 7: Dimensionless pressure of four-tilting-pad bearing in different [figure omitted; refer to PDF] : (a) 0.1, (b) 0.5, (c) 1, and (d) 5.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

With these parameters, the axis center lies in first quadrant, and the original eccentricities of the pads in second quadrant and fourth quadrant are the same; the oil film pressure distribution of pads #1 and #3 in second and fourth quadrant should agree. The original eccentricity in third quadrant is relatively larger than those in other quadrants, and so the oil film pressure of pad #2 is relatively smaller. Similarly, the oil film pressure of pad #4 in the fourth quadrant is the largest.

As shown in Figure 7, the dimensionless pressure of oil film increases as the length-diameter ratio increases. The oil film pressure distribution in axial direction exhibits as parabola form when the length-diameter ratio is smaller; the axial distribution of oil film pressure is generally considered consistent when the length-diameter ration reaches 5.

Figure 8 is the displacement of journal center within bearing clearance for different operating conditions, where, [figure omitted; refer to PDF] the length-diameter ratio of tilting pad bearing is 0.5; [figure omitted; refer to PDF] the preload coefficient is 0.5; and, [figure omitted; refer to PDF] , the clearance ratio is 3[per thousand].

Figure 8: Locus of journal center for tilting pad journal bearing.

(a) Static equilibrium position

[figure omitted; refer to PDF]

(b) Partial enlarged view

[figure omitted; refer to PDF]

The static equilibrium position curve of a tilting pad journal bearing is approximately a line segment along the direction of bearing center, which is quite different from the balanced semicircle of circular-pad journal bearing in an equilibrium position. The static equilibrium position for each state is closer to [figure omitted; refer to PDF] -axis, because the pad of a journal bearing tilts as the journal center changes with weak horizontal and vertical couplings. As shown in Figure 8(b), the static balance curve is not exactly a straight line but a line with weak fluctuation, which means the oil film forces in horizontal and vertical directions exhibit a weak coupling relationship.

In order to validate the accuracy of the present oil film force model, comparisons have been made with the finite difference method model, and the boundary conditions are consistent in the models; oil film forces for different models with different eccentricities were investigated. Based upon Figure 9, the total oil film force can be derived as follows: [figure omitted; refer to PDF]

Figure 9: Dimensionless force versus original eccentricity for different [figure omitted; refer to PDF] : (a) 0.01, (b) 0.5, (c) 1, and (d) 5.

(a) [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

Figure 9 shows the comparison of the total dimensionless oil film force versus original eccentricity from finite difference model with the current approach for different length-diameter ratio, where eccentricity is 0.01-0.49, the pad preload coefficient [figure omitted; refer to PDF] , bearing clearance ratio [figure omitted; refer to PDF] , pad pivot offset is 50%, pad arc angle is 80°, attitude angle is 45°, and length-diameter ratios [figure omitted; refer to PDF] , 0.5, 1, and 5.

As shown in Figure 9, the variation of approximate analysis model proposed in the paper is the same as that of the finite difference model for different length-diameter ratios, which means the model proposed in this paper is correct and reliable.

7. Conclusions

The oil film characteristics of a single tilting pad are analyzed in this paper, and an oil film force model is established for a single pad under local coordinates, the relations between local and global coordinates are discussed, and a nonlinear oil film force model is proposed for a four-tilting-pad journal bearing.

The model proposed in this paper was proved to be correct and accurate by analyzing the distribution of oil film pressure and by comparing the model proposed in this paper with the model established using finite difference method.

Acknowledgment

This work is funded by the Fundamental Research Funds for the Central Universities (Grant no. 3132015026).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] T. Dimond, A. Younan, P. Allaire, "A review of tilting pad bearing theory," International Journal of Rotating Machinery , vol. 2011, 2011.

[2] M. Nasselqvist, R. Gustavsson, J.-O. Aidanpaa, "Experimental and numerical simulation of unbalance response in vertical test rig with tilting-pad bearings," International Journal of Rotating Machinery , vol. 2014, 2014.

[3] D. P. Tschoepe, D. W. Childs, "Measurements versus predictions for the static and dynamic characteristics of a four-pad, rocker-pivot, tilting-pad journal bearing," Journal of Engineering for Gas Turbines and Power , vol. 136, no. 5, 2014.

[4] L. S. Andres, Y. J. Tao, "The role of pivot stiffness on the dynamic force coefficients of tilting pad journal bearings," Journal of Engineering for Gas Turbines and Power , vol. 135, no. 11, 2013.

[5] D. W. Childs, J. Harris, "Static performance characteristics and rotordynamic coefficients for a four-pad ball-in-socket tilting pad journal bearing," Journal of Engineering for Gas Turbines and Power , vol. 131, no. 6, 2009.

[6] J. M. Cao, T. Dimond, P. Allaire, "Nonlinear modeling of tilting-pad bearings with application to a flexible rotor analysis," in Proceedings of the ASME Design Engineering Technical Conference, of Paper no. DETC2013-13712, Portland, Ore, USA, August 2013.

[7] J. W. Lund, "Spring and damping coefficients for the tilting-pad journal bearing," ASLE Transactions , vol. 7, pp. 342-352, 1964.

[8] J. C. Nicholas, E. J. Gunter, P. E. Allaire, "Stiffness and damping coefficients for the five-pad tilting-pad bearing," ASLE Transactions , vol. 22, no. 2, pp. 113-124, 1979.

[9] J. C. Nicholas, "Lund's tilting pad journal bearing pad assembly method," Transactions of the ASME, Journal of Vibration and Acoustics , vol. 125, no. 4, pp. 448-454, 2003.

[10] S. Pagano, E. Rocca, M. Russo, R. Russo, "Dynamic behaviour of tilting-pad journal bearings," Proceedings of the Institution of Mechanical Engineers J: Journal of Engineering Tribology , vol. 209, no. 4, pp. 275-285, 1995.

[11] R. Brancati, E. Rocca, R. Russo, "Non-linear stability analysis of a rigid rotor on tilting pad journal bearing," The American Society of Lubrication Engineers , vol. 29, no. 7, pp. 571-578, 1996.

[12] Y. Hou, T. Lai, S. Chen, B. Ma, J. Liu, "Numerical analysis on the static performance of tilting pad-journal gas bearing in subsystems," Tribology International , vol. 61, pp. 70-79, 2013.

[13] W. Wang, Z. M. Wang, "Nonlinear oil film force database," Journal of Shanghai University of Technology , vol. 14, no. 4, pp. 299-305, 1993.

[14] E. P. Okabe, K. L. Cavalca, "Rotordynamic analysis of systems with a non-linear model of tilting pad bearings," in Proceedings of the 7th International Conference on Rotor Dynamics (IFToMM '06), of Paper No. 222, Technische Universitat Wien, Vienna, Austria, 2006.

[15] E. P. Okabe, K. L. Cavalca, "Rotordynamic analysis of systems with a non-linear model of tilting pad bearings including turbulence effects," Nonlinear Dynamics , vol. 57, no. 4, pp. 481-495, 2009.

[16] Y. L. Wang, Z. S. Liu, "Study on tilting characteristic of pads in tilting pad bearing," in Proceedings of the International Conference on Electronic and Mechanical Engineering and Information Technology (EMEIT '11), vol. 9, pp. 4653-4657, Harbin, China, August 2011.

[17] Y. L. Wang, Z. S. Liu, W. J. Kang, J. J. Yan, "Approximate analytical model for fluid film force of finite length plain journal bearing," Proceedings of the Institution of Mechanical Engineers Part C: Journal of Mechanical Engineering Science , vol. 226, no. 5, pp. 1345-1355, 2012.

Appendix

Appendix

Consider [figure omitted; refer to PDF]