Transformation between polar and rectangular coordinates of stiffness and dampness parameters in hydrodynamic journal bearings

The stiffness and dampness parameters of journal bearings are required in rectangular coordinates for analyzing the stability boundary and threshold speed of oil film bearings. On solving the Reynolds equation, the oil film force is always obtained in polar coordinates; thus, the stiffness and dampness parameters can be easily obtained in polar coordinates. Therefore, the transformation between the polar and rectangular coordinates of journal bearing stiffness and dampness parameters is discussed in this study.


Introduction
In modern industry, rotating parts of engineering equipment are supported by journal bearings in the vertical direction. The machine characteristics are significantly dependent on the performances of journal bearings. Due to the oil whip effect in a rotating hydrodynamic journal bearing, self-excited vibration occurs in the oil film, which increases with an increase in the rotation speed [1][2][3][4]. As a result of the self-excited vibration, threshold speed and stability boundary exist for the rotating bearings. When the rotating speed is larger than the threshold speed, the vibration results in large orbiting amplitudes of the journal and leads to the contact between the journal and bearing, causing bearing failure. Further, the bearing is stable at the journal bearing center in the stability boundary. Khonsari and Chang [5] analyzed the nonlinear stability of journal bearings, and obtained the stability boundary by tracking the journal center trajectory. To easily derive the linear threshold speed, Huang et al. [6] chose polar coordinates instead of Cartesian coordinates to signify the state vector. Lin et al. [7,8] discussed the threshold speed and stability boundary of hydrodynamic journal bearings lubricated using non-Newtonian fluids. Kushare and Sharma [9] dealt with the threshold speed of two lobe symmetric hole entry worn hybrid journal bearing by theoretically considering the non-Newtonian behavior of the lubricants.
The threshold speed and stability boundary of the hydrodynamic journal bearing are determined based on the stiffness and damping coefficients. Lund and Thomsen [10] proposed a method to study these coefficients. The Reynolds equation is derived with respect to Cartesian coordinates and stiffness and damping coefficients are obtained by integrating the new equations. Based on this method, Jang and Kim [11] studied the dynamic characteristics of journal bearings with five degrees of freedom (DOF) and derived the perturbation equations. Crooijmans et al. [12] discussed the self-excited vibration of hydrodynamic journal bearings, and interpreted the dynamic performance

Subscripts and superscripts
gz gradient on z direction s stability state ε component on eccentric direction ϕ component perpendicular to eccentric direction * non-dimensional parameter of rotor-bearing systems. Wang and Khonsari [13] proposed a method to calculate the dynamic coefficients to obtain an improved physical meaning. It is known that oil stiffness and damping characteristics are important for linear threshold speed and stability boundary of hydrodynamic bearings.
The coordinate system applied in the general Reynolds equation is polar coordinates. Thus, the oil film force is always obtained in polar coordinates and the stiffness and dampness parameters can be easily obtained in polar coordinates. While discussing linear threshold speed and stability boundary of hydrodynamic bearings, the stiffness and dampness parameters in rectangular coordinates are required. As per our understanding, the coordinate system applied in the previous study on linear threshold speed and stability boundary was either polar coordinates or rectangular coordinates for the complete calculation, which com-plicated the analysis. To address this complication, the transformation between polar and rectangular coordinates of journal bearing stiffness and dampness parameters is discussed in this study. The stiffness and dampness parameters are first calculated in polar coordinates, and then, converted into rectangular coordinates. Further, the stiffness and dampness parameters in rectangular coordinates can be applied to calculate linear threshold speed and stability boundary of hydrodynamic bearings. Figure 1 shows the journal bearing structure. The journal rotates in an anticlockwise direction with an angular speed ω and a radius R. The radial clearance is expressed as C, eccentricity is expressed as e, angular  The non-dimensional parameters are as follows:

Transformation between polar and rectangular coordinates
Threshold speed is an important parameter for the stability of the oil film bearing, and the non-dimensional linear threshold speed of the journal bearing  * s is as follows [5,8] where * ij K and * ij C (i, j = X, Y) are the non-dimensional stiffness and dampness parameters of the bearing in rectangular coordinates, respectively.
Although the results of oil film forces in journal bearings have always been expressed in polar coordinates, in this study, the transformation between polar and rectangular coordinates of journal bearing stiffness and dampness parameters is discussed as follows: The non-dimensional oil film forces are expended in the X and Y directions as a first-order Taylor series approximation, as shown in Eq. (3).
where K r and C r are the non-dimensional stiffness and dampness matrixes in rectangular coordinates, respectively. The expressions for K r and C r are as follows: In Fig. 1, the non-dimensional resultant forces in the X and Y directions are as follows: where S is the Sommerfeld number of the bearing. At the equilibration position, the expression of S is where the subscript s denotes the equilibrium position. The differential of Eqs. (5) and (6) is calculated as where A is a transformation matrix as follows: The differential of non-dimensional oil film forces 204 in the ε and ϕ directions are where K p and C p are the non-dimensional stiffness and dampness matrixes in rectangular coordinates, respectively. The expressions of K p and C p are By substituting Eq. (10) in Eq. (8), we get where K p0 is the non-dimensional stiffness correction matrix as follows: The coordinate transformations between (X * , Y * ) and (ε, ϕ) at the equilibrium position are When Eqs. (14) and (15)  Equations (17) and (18) are the general transformation of the hydrodynamic journal bearing stiffness and dampness coefficients between polar and rectangular coordinates, respectively. The transformation validity for short journal bearings is discussed in Section 3.

Verification and discussion
When the short-bearing approximate is applied with half-Sommerfeld boundary conditions, the nondimensional oil film forces of the bearing are as follows [6,13].