Theoretical Study on the Dynamic Characteristics of Marine Stern Bearing Considering Cavitation and Bending Deformation Effects of the Shaft

: When the ship runs, owing to the superposition of the gravity of the shaft and resistance of water, with the increment in rotational speeds, the shaft will produce different degrees of bending deformation, which immensely reduces the power transmission efﬁciency. Based on the aforementioned problem, the present study focuses on the inﬂuences of bending deformation of the shaft with a cavitation effect on the dynamic characteristics of the stern bearing. The mixed lubrication model with bending deformation and cavitation effect is established. At present, the deﬂection curve equation is employed, the ﬁnite perturbation method is applied to calculate the dynamic coefﬁcient, and the cavitation pressure is determined by the numerical method. According to the analysis, the variation laws of equivalent stiffness and natural frequency are exhibited. It is shown that the equivalent stiffness is more affected by the speeds, especially at low speeds; There is a critical speed between 130 rpm and 150 rpm, which makes the natural frequency strike the maximum value. Finally, the research results provide a theoretical basis for the ships to avoid large vibration during navigation.


Introduction
The lubricants can be employed to form the hydrodynamic, decrease the friction and wear in bearings, intensify load-carry capacity and reduce the mechanical vibration. For the journal bearing, the lubricant forms the hydrodynamic film to support the external loads when the three conditions (the wedge gap, relative velocity and the viscous fluid) are met. In addition, the lubrication mode of journal bearings in ships is mainly water lubrication. Therefore, the present research focuses on water-lubricated bearings.
The studies on water-lubricated bearings have been implemented by plenty of scholars, and the relevant vibration performances and lubrication behaviors are acquired and analyzed. As the crucial part of the submarine or ship, the stiffness and damping properties have a crucial effect on the critical velocity and the stability of ships (as can be seen in Figure 1). For simplifying the analysis on stability, the bearing forces are linearized and adopted to be the linear function for velocity and displacement. According to the method mentioned above, there are eight dynamic coefficients, including four damping coefficients and four stiffness coefficients [1]. For these coefficients, plenty of works are implemented by scholars, and some typical conclusions and results are proposed. Sternlicht [2] expressed the bearing force vector as the linear function of the displacements and static bearing force found the equilibrium position. Sternlicht [2] carried out the relevant investigation. The force vector is adopted as the linear function of static bearing force and displacement to locate the equilibrium point. In his study, by using the Finite Difference Method (FDM), four stiffness coefficients and two direct damping coefficients were given. For the dynamic investigation. The force vector is adopted as the linear function of static bearing force and displacement to locate the equilibrium point. In his study, by using the Finite Difference Method (FDM), four stiffness coefficients and two direct damping coefficients were given. For the dynamic characteristics, a perturbation solution method of the Reynolds equation is proposed by Lund [3,4]. As noted, the method decreases the demand for numerical differentiation. The first order perturbation solution is employed in Reinhardt's study [5]. For the dynamic characteristics of the journal bearing with grooves, two perturbation methods (infinitesimal perturbation method, IFPM and finite perturbation method, FPM) are applied in Someya's investigation [6]. Nowadays, IFPM and FPM are frequently employed for analyzing dynamic characteristics [7][8][9][10]. Due to the shortcomings of conventional methods in studying the fluid-solid interaction of the bearing-rotor system, Lin [11] explored the two-way fluid-solid interaction and hydrodynamic lubrication performances with thermal effect and cavitation effect by computational fluid dynamics (CFD). In addition, the variation laws of dynamic coefficients under the misalignment are also explored [12][13][14][15]. The influences of the misalignment on the dynamic characteristics are shown in the research of Bou-Said [16] and San Andres [17]. In Qiu's research [18], for the static loads, the static and dynamic characteristics (attitude angle, friction force, moment and dynamic coefficients, etc) are analyzed under different misalignments and eccentricity ratios. For the water-lubricated bearings, the static and dynamic characteristics are investigated by Feng [19] and Hanawa [20]. Thermohydrodynamic, turbulent and inclined effects were taken into account in the revised models. Lin [21] and Liang [22] performed research on the dynamic characteristics of water-lubricated bearings with grooves. Cavitation effects were considered for the lubrication model. Research indicated that cavitation had more effects on the stiffness coefficients than the damping coefficients.
Therefore, based on the mentioned above, the studies on the influences of cavitation effect on journal bearing are abundant. However, for the coupling effect of bending deformation and cavitation effect, there is little literature. Consequently, a related investigation is carried out.
When the shaft runs, the journal deformation is caused by the force acting on it. Subsequently, the deformation has significant effects on the pressure profile and the asperity contact force under a mixed lubrication condition. Meanwhile, the hydrodynamic force and the asperity contact force jointly act on the journal and affect the deformation. In addition, when the bearing width is large, the deformation of the shaft cannot only be described by two misalignment angles (β and γ) [14]. With the exception of the influences of pressure, asperity contact force and bending deformation of the shaft, the shaft also supports the concentrated load from the propeller. In this study, the interaction between lubrication properties and deformation journals has been taken into consideration. The study applies the deflection curve equation to express the bending deformation. The finite perturbation method is selected to simplify the mathematical model. Moreover, the problem of the mixed lubrication with JFO cavitation condition on a marine stern tube bearing is presented. The effects of bending deformation of the propeller, as well as lubrication boundary conditions, on the dynamic coefficients, are analyzed. The first order perturbation solution is employed in Reinhardt's study [5]. For the dynamic characteristics of the journal bearing with grooves, two perturbation methods (infinitesimal perturbation method, IFPM and finite perturbation method, FPM) are applied in Someya's investigation [6]. Nowadays, IFPM and FPM are frequently employed for analyzing dynamic characteristics [7][8][9][10]. Due to the shortcomings of conventional methods in studying the fluid-solid interaction of the bearing-rotor system, Lin [11] explored the two-way fluid-solid interaction and hydrodynamic lubrication performances with thermal effect and cavitation effect by computational fluid dynamics (CFD). In addition, the variation laws of dynamic coefficients under the misalignment are also explored [12][13][14][15]. The influences of the misalignment on the dynamic characteristics are shown in the research of Bou-Said [16] and San Andres [17]. In Qiu's research [18], for the static loads, the static and dynamic characteristics (attitude angle, friction force, moment and dynamic coefficients, etc) are analyzed under different misalignments and eccentricity ratios. For the water-lubricated bearings, the static and dynamic characteristics are investigated by Feng [19] and Hanawa [20]. Thermohydrodynamic, turbulent and inclined effects were taken into account in the revised models. Lin [21] and Liang [22] performed research on the dynamic characteristics of water-lubricated bearings with grooves. Cavitation effects were considered for the lubrication model. Research indicated that cavitation had more effects on the stiffness coefficients than the damping coefficients.
Therefore, based on the mentioned above, the studies on the influences of cavitation effect on journal bearing are abundant. However, for the coupling effect of bending deformation and cavitation effect, there is little literature. Consequently, a related investigation is carried out.
When the shaft runs, the journal deformation is caused by the force acting on it. Subsequently, the deformation has significant effects on the pressure profile and the asperity contact force under a mixed lubrication condition. Meanwhile, the hydrodynamic force and the asperity contact force jointly act on the journal and affect the deformation. In addition, when the bearing width is large, the deformation of the shaft cannot only be described by two misalignment angles (β and γ) [14]. With the exception of the influences of pressure, asperity contact force and bending deformation of the shaft, the shaft also supports the concentrated load from the propeller. In this study, the interaction between lubrication properties and deformation journals has been taken into consideration. The study applies the deflection curve equation to express the bending deformation. The finite perturbation method is selected to simplify the mathematical model. Moreover, the problem of the mixed lubrication with JFO cavitation condition on a marine stern tube bearing is presented. The effects of bending deformation of the propeller, as well as lubrication boundary conditions, on the dynamic coefficients, are analyzed.

Characteristics of Journal Bearing
A typical lubricated cylindrical journal bearing is displayed in Figure 2. According to the references describing the lubrication properties of journal bearing, it is obvious that Lubricants 2022, 10, 242 3 of 18 the hydrodynamic pressure rises to a peak in the converging film region and drops to the atmospheric pressure level at both ends. However, in some special areas where the film thickness partially raises, the pressure descends to an ambient level on the contrary. The water evaporation is caused by the released gas from the lubricant or the ambient pressure below the vapor pressure. The rupture of water film is a crucial feature for the stability of the bearing, and the phenomenon is called lubricant cavitation. The research [23,24] explored the effects of cavitation on stability and bearing performances.

Characteristics of Journal Bearing
A typical lubricated cylindrical journal bearing is displayed in Figure 2. According to the references describing the lubrication properties of journal bearing, it is obvious that the hydrodynamic pressure rises to a peak in the converging film region and drops to the atmospheric pressure level at both ends. However, in some special areas where the film thickness partially raises, the pressure descends to an ambient level on the contrary. The water evaporation is caused by the released gas from the lubricant or the ambient pressure below the vapor pressure. The rupture of water film is a crucial feature for the stability of the bearing, and the phenomenon is called lubricant cavitation. The research [23,24] explored the effects of cavitation on stability and bearing performances. Because cavitation collapse can cause severe surface material damage, the cavitation phenomenon should be tried as a very important study direction. Moreover, the beginning and degree can change the load-carrying capacity. The cavitation has a significant effect on the stability of the bearing-rotor system and the maximum amplitude of vibration and whirl [25].
In common, the main function of hydrodynamic pressure is that it provides the wedge space and squeeze effect. The JFO model applies to the dynamic loads that the surface is at the squeeze motion of film [26]. The time variation parameters of rupture film are included in the continuity equation but the bubble dynamics are not considered. With the fluctuation of surface squeeze velocity, the shape of cavitation changes.
In the JFO model, however, since the cavitation zone is unknown during the operating condition, is hard to solve the Reynolds equation numerically. A current cavitation algorithm is developed by Elrod [27,28]. In his research, JFO model has included the Reynolds equation. For the model, the fluid bulk modulus (κ) is used to relate lubricant pressure and density. The boundary conditions at the cavitation area are automatically met by the switch function. Moreover, under cavitation area and integrity film, the feature of continuity equation is transformed from ellipse to parabola by the switch function. Then a variable (θ in [29,30]) is introduced to indicate different lubrication regions (full film region and cavitation).
For the compressible liquid, the relationship between density and pressure is shown as follows: And define the density ratio as: Because cavitation collapse can cause severe surface material damage, the cavitation phenomenon should be tried as a very important study direction. Moreover, the beginning and degree can change the load-carrying capacity. The cavitation has a significant effect on the stability of the bearing-rotor system and the maximum amplitude of vibration and whirl [25].
In common, the main function of hydrodynamic pressure is that it provides the wedge space and squeeze effect. The JFO model applies to the dynamic loads that the surface is at the squeeze motion of film [26]. The time variation parameters of rupture film are included in the continuity equation but the bubble dynamics are not considered. With the fluctuation of surface squeeze velocity, the shape of cavitation changes.
In the JFO model, however, since the cavitation zone is unknown during the operating condition, is hard to solve the Reynolds equation numerically. A current cavitation algorithm is developed by Elrod [27,28]. In his research, JFO model has included the Reynolds equation. For the model, the fluid bulk modulus (κ) is used to relate lubricant pressure and density. The boundary conditions at the cavitation area are automatically met by the switch function. Moreover, under cavitation area and integrity film, the feature of continuity equation is transformed from ellipse to parabola by the switch function. Then a variable (θ in [29,30]) is introduced to indicate different lubrication regions (full film region and cavitation).
For the compressible liquid, the relationship between density and pressure is shown as follows: And define the density ratio as: Then Equation (2) can be rewritten as: where g is defined as a switch function in [29,30] and values of g in different zones are described as: For Equation (3), the direct integration is implemented, and then the following equation can be obtained: It should be indicated that for the part of whole film, g = 1 and P > P cav , since α = ρ ρ cav > 1; and for cavitation area P = P cav because of g = 0. As noted, pressure gradient is a crucial factor for the thin film flow. Consequently, for P 1 and P 2 , they are larger than P cav , Owing to the large magnitude of κ, the slight difference in density ratio can result in a tremendous difference in pressure. In addition, for the numerical model, there are a few difficulties because of the different considerations. For the variable α, under the whole film area and cavitation area, the explanation is different [31]. In cavitation area, due to the vapor material within the cavity, the density of water is even (ρ cav ) and the wedge gap is not stuffed. Therefore, α is recognized as the small film content and the void fraction can be expressed by the (1 − α).
Based on the mass flow conversation of thin film: Under the full film area, for the mass flow rate of fluid, it can be shown as: However, for the cavitation position, ∂P ∂x = ∂P ∂z = 0, the mass flow rate is: Due to ∂P = κ ρ ∂ρ = g κ α ∂α, Equation (7) can be derived as: As noted, under the full film area (g = 1 and α = ρ/ρ cav ) and cavitation position, g = 0 with α is recognized as the small film content.
Equation (6) can be solved by the global mass conservation as: And in the cavitation zone, Equation (10) can be reduced to: This establishes a dynamic flow balance in the cavitation area. The mixed success is available in the present research by the common cavitation algorithm. However, compared with the rough technology, the common cavitation algorithm is obsolete.

Bearing Forces
There is a crucial relationship between rotational speed and journal position. The relation equation can be expressed by the bearing force.
By numerical procedure, as long as the film pressure is determined, the pressure in a horizontal and vertical direction is integrated to acquire the bearing force:

Attitude Angle
The angle between the load direction and the central line that connects the bearing and journal is the attitude angle. Under stable condition, the bearing forces are employed to balance the external loads. According to the external load W, the bearing force vector is shown by the load-carrying capacity vector: If the W acts on the vertical direction (such as the weight of the screw propeller in this study), the horizontal force is The error function is defined as: If dφ k 0 = 0, the condition of Equation (15) is met automatically. In order to calculate dφ k 0 = 0, the Newton-Raphson iteration method is adopted: The φ k 0 is considered the attitude angle φ 0 . In the numerical procedure, if Once φ 0 is converged, P i,j and water film boundaries can be determined.

Film Thickness Considering Bending Deformation of Journal
When the relative rotating speed exists between the journal and bearing, the eccentricity amount adjusts automatically until the film pressures calculated by the average Reynolds equation and the asperities contact pressures calculated by the contact model balance the external loads. From the description of lubrication properties, the film pressures and the contact pressures depend on the local film thickness. In addition, the misalignment parameters give the expression of film thickness of misaligned bearing (β and γ) [14], as Equation (19) shows: When it comes to the large width of the bearing, the deflect angle in each journal section is different as a result. After considering the bending deformation of the shaft, the film thickness can be briefly expressed as: Define e i = (υ ix + e ix ) 2 + υ iy + e iy 2 , when the shaft is bending deformation, for the ith section, the eccentricity ratio is e i ; e ix and e iy . These are the components of the eccentricity ratio without bending in horizontal (x) and vertical (y) directions (as shown in Figure 3). The deformation of this section is υ i = υ ix 2 + υ iy 2 , and in order to determine υ i accurately, the displacement superposition method is applied to calculate the bending deformation in this study (as shown in Figure 4).
sures and the contact pressures depend on the local film thickness. In addition, the misalignment parameters give the expression of film thickness of misaligned bearing ( β and γ ) [14], as Equation (19) shows: When it comes to the large width of the bearing, the deflect angle in each journal section is different as a result. After considering the bending deformation of the shaft, the film thickness can be briefly expressed as: , and in order to determine i υ accurately, the displacement superposition method is applied to calculate the bending deformation in this study (as shown in Figure 4).

Determining the Bending Deformation of Journal
Obviously, the force that acts on the shaft results in the bending deformation of the shaft. Due to the reaction relationship, the deformation will change the profile of film pressure. If the ratio of film thickness to surface roughness is less than 4 [30], the asperity contact force of the end is seriously affected. In [14], Sun used γ to describe the misaligned angle of the shaft in bearing, and it can be expressed as: However, as per the reason mentioned in Section 2.4, it should be noted that the mis-

Determining the Bending Deformation of Journal
Obviously, the force that acts on the shaft results in the bending deformation of the shaft. Due to the reaction relationship, the deformation will change the profile of film pressure. If the ratio of film thickness to surface roughness is less than 4 [30], the asperity Lubricants 2022, 10, 242 7 of 18 contact force of the end is seriously affected. In [14], Sun used γ to describe the misaligned angle of the shaft in bearing, and it can be expressed as: However, as per the reason mentioned in Section 2.4, it should be noted that the misaligned angle in each section of the shaft is various in the very wide bearing. During this study, in the axial direction, the whole bearing is divided into 20 sections. For the sake of accurate determination of the bending deformation, the differential equation of the deflection curve (given as following) is used: When it comes to the loaded stern shaft shown in Figure 3. Under sophisticated load conditions, owing to the small bending deformation of the shaft, the shaft material is recognized as elastic, based on the conclusion of [31]. The displacement superposition method can be applied in this study to solve bending deformation. From superposition theory, in the systems subject to multiple loads, resultant stress or strain is the algebraic sum of its effects when the load is applied alone. Taking a cantilever beam as an example, the solution of this beam under typically single-load conditions is shown in Table 1. The sum of these solutions can be also employed to express the deformation caused by external loads and asperity contact force.

Determining the Bending Deformation of Journal
Obviously, the force that acts on the shaft results in the bending deformation of the shaft. Due to the reaction relationship, the deformation will change the profile of film pressure. If the ratio of film thickness to surface roughness is less than 4 [30], the asperity contact force of the end is seriously affected. In [14], Sun used γ to describe the misaligned angle of the shaft in bearing, and it can be expressed as: However, as per the reason mentioned in Section 2.4, it should be noted that the misaligned angle in each section of the shaft is various in the very wide bearing. During this study, in the axial direction, the whole bearing is divided into 20 sections. For the sake of accurate determination of the bending deformation, the differential equation of the deflection curve (given as following) is used: When it comes to the loaded stern shaft shown in Figure 3. Under sophisticated load conditions, owing to the small bending deformation of the shaft, the shaft material is recognized as elastic, based on the conclusion of [31]. The displacement superposition method can be applied in this study to solve bending deformation. From superposition theory, in the systems subject to multiple loads, resultant stress or strain is the algebraic sum of its effects when the load is applied alone. Taking a cantilever beam as an example, the solution of this beam under typically single-load conditions is shown in Table 1. The sum of these solutions can be also employed to express the deformation caused by external loads and asperity contact force.

6
Pz l z EI l Concentrated load P at end As described previously, in the axis direction, the journal is divided into 20 rigid nodes, and in z direction, for the ith node, the coordinate is According to Equations (15) and (16), for all nodes of the journal, the bending deformation is displayed as: As described previously, in the axis direction, the journal is divided into 20 rigid nodes, and in z direction, for the ith node, the coordinate is Z i . For the bending deformation of the shaft, the ith section is υ ix and υ iy . For the hydrodynamic force and asperity contact force, at the ith section center, they are expressed as F ix and F iy , W ix and W iy . As noted, the subscript of x indicates the x direction, for the subscript of y, the meaning is identical. In addition, the resultant forces of both are F i and W asp , respectively.
According to Equations (15) and (16), for all nodes of the journal, the bending deformation is displayed as: when z = z n , According to the above expressions, the bending deformation at mixture load conditions, including the concentrated loads of the propeller, the hydrodynamic pressure and asperity pressures, can be determined accurately.

The Deformation of Bearing
As an external load is applied to a journal bearing, the bending deformation of the journal has happened, and as Figure 5 shows, the water film thickness is reduced at the sections close to the bearing end; as a result, the film pressure will increase rapidly at those sections. In this situation, in order to analyze the consequences of the pressure distribution and wear phenomena better, the method that the elastic deformation of bearing is included in the elastohydrodynamic model is crucial. At present, the elastic deformation is incorporated into the lubrication characteristics of dynamic bearing that analyzes the various external loads. By calculating the deformation of bearing ∆H in different conditions and adding it to the film thickness value (as shown in Equation (20)), the influences of elastic deformation on lubrication performances are explored. In this study, the equivalent model is employed to describe the deformation. From Figure 6, the whole elastic deformation of journal and bearing is expressed by H Δ , and 1 E , 2 E is Young's modulus of journal and bearing, respectively. The composite of Young's modulus of bearing and journal is E ′ . The deformation corresponds to that of the bearing surface just with composite Young's modulus. According to the equivalent model, the changing length of the object just with Young's modulus E ′ and initial length L can be expressed by H Δ . According to the above analysis, the expression [32] is obtained as: In this study, the equivalent model is employed to describe the deformation. From Figure 6, the whole elastic deformation of journal and bearing is expressed by ∆H, and E 1 , E 2 is Young's modulus of journal and bearing, respectively. The composite of Young's modulus of bearing and journal is E . The deformation corresponds to that of the bearing surface just with composite Young's modulus. According to the equivalent model, the changing length of the object just with Young's modulus E and initial length L can be expressed by ∆H. According to the above analysis, the expression [32] is obtained as: Lubricants 2022, 10, 242 9 of 18 the bearing surface just with composite Young's modulus. According to the equivalent model, the changing length of the object just with Young's modulus E ′ and initial length L can be expressed by H Δ . According to the above analysis, the expression [32] is obtained as: The normal pressure acting on the journal can be acquired by:

The Pressure in Cavitation Zone
To determine the pressure in the cavitation zone, the equations listed in reference [33] are used in this study, which is given as: The normal pressure acting on the journal can be acquired by:

The Pressure in Cavitation Zone
To determine the pressure in the cavitation zone, the equations listed in reference [33] are used in this study, which is given as: As it is summarized in reference [33], Equation (29) could give a better correlation in the transition zone. Based on those equations, the P cav can be quickly determined by ρ v , ρ l , C v , C l . Define:

Calculation of the Dynamic Coefficients
Based on Equation (12), the bearing forces are derived by: And the eight coefficients are defined as: These eight coefficients are defined as dynamic coefficients of the journal bearing, as shown in Figure 7. According to those references reviewed previously, there are two theoretical methods (IFPM and FPM) to calculate the dynamic coefficients. FPM is used in this study.

,
, , These eight coefficients are defined as dynamic coefficients of the journal bearing, as shown in Figure 7. According to those references reviewed previously, there are two theoretical methods (IFPM and FPM) to calculate the dynamic coefficients. FPM is used in this study.

Finite Perturbation Method (FPM)
FPM perturbs the journal under the equilibrium position by the small ∆x, ∆y, ∆ .
x and ∆ . y, and the force finite difference is employed to describe the force coefficient. For the calculation accuracy, the perturbation is implemented in the positive and negative directions, and the partial differential in both directions is adopted to calculate the partial derivative. At the equilibrium position (x o , y o ), the partial derivative ∂F x ∂x (force coefficient k xx ) can be solved by [7]: At the steady state, .
x, 0), (x 0 , y 0 , ∆ . x, 0), (x 0 , y 0 , 0, −∆ . y) and (x 0 , y 0 , 0, ∆ . y). If the perturbation magnitude is large, there are some errors in the linear model. This is because the bearing force is a sophisticated nonlinear function. According to the theory, there is a maximum perturbation amplitude. When the amplitude is less than the maximum value, the calculation accuracy is high and acceptable. Applying FPM, the maximum perturbation amplitude can be determined by a given accuracy.

Numerical Procedure
Based on the numerical model developed in this study, the solution procedure is shown in Figure 8. In the dynamic part, a combination of FDM and FPM methods is applied. The equation of film thickness is displayed by the nodal global coordinate of the bearing surface. For the analysis of lubrication performances, the computer program is developed. The flowchart of the numerical procedure can be seen in Figure 8. The calculation procedure is designed as follows: (1) Input the parameters of the stern bearing, initial film pressure, asperities contact forces and deformation of bearing. The assumption condition for the shaft center is given, and then the deformation of the shaft under external loads (the concentrated load of the propeller) is calculated; (2) Calculating the water film shape according to the current journal position through Equation (20) and adding the value of bearing deformation into h iθ ; (3) Obtain the hydrodynamic pressure distribution by implementing the universal cavitation algorithm; Calculate the asperities contact forces; (4) The trajectory position is acquired by the integral for the balance loads, hydrodynamic pressure and asperity contact force. In addition, the over-relaxation Newton-Raphson method is employed to accelerate convergence; (5) For the deformation of the shaft under external loads, film pressure and asperity contact force and bearing deformation, the calculation needs to be resolved until they meet the convergence criterion: (6) In order to determine the bearing force of each perturbation position, step (2)~(5) is recalculated by perturbing the position and velocity of the journal; (7) According to Equations (38)-(45), the dynamic coefficients are acquired.
In the analysis, the prediction is mainly decided by the value of the volume modulus (κ). For the accurate solution, the algebraic equations are stiff because of the actual value of the bulk modulus of water. One reasonable explanation is that the slight variation of density has a significant effect on the hydrodynamic pressure at the full film area. Therefore, for the numerical model, the round-off errors appear [7]. The artificial low bulk modulus is adopted to avoid the problem of accuracy and convergence. As noted, its magnitude is less than the physical value. The excellent results are given by the stable numerical algorithm. The numerical range is 1/100 to 1/10 of the actual physical value.

Validation
To validate the present numerical scheme, a result obtained from the present analysis is compared with the results obtained by the experimental data given by Su [34]. Table 2 lists the bearing parameters for the comparison case, and Figure 9 shows the film pressure predicted by the model developed in this study and the experimental results from [34]. As shown in this figure, compared with the results gotten from the Reynolds boundary, the film pressure predicted by the JFO boundary matches the experimental results much better. The pressure is recognized as vapor pressure in the Reynolds boundary condition. In addition, at the position of full film rupture, the circumferential gradient must be equal to zero. From the above analysis, compared with the experiment results, the error is large. Numerical research on the water-lubricated bearing with cavitation is carried out to explore the influences of propeller deformation on the behaviors of the bearing. For the following situation, the water-lubricated stern bearing (shown in Figure 10) variables used in the computation are listed in Table 3. The weight of the screw propeller is 0.21498 × 10 6 N. According to the validation case, the JFO boundary condition is much closer to the real condition. Therefore, the following analysis is only based on the JFO boundary condition. Meanwhile, the deformations of the bearing and propeller are taken into consideration. The structure of the propeller is shown in Figure 10. The shaft is divided into two parts by the flange. According to the position of the bearings, the whole shaft can be divided into five segments, as shown in Figure 10. The length of each part is listed in Table 3, together with the bearings' parameters and material properties.
To validate the present numerical scheme, a result obtained from the present analysis is compared with the results obtained by the experimental data given by Su [34]. Table 2 lists the bearing parameters for the comparison case, and Figure 9 shows the film pressure predicted by the model developed in this study and the experimental results from [34] As shown in this figure, compared with the results gotten from the Reynolds boundary the film pressure predicted by the JFO boundary matches the experimental results much better. The pressure is recognized as vapor pressure in the Reynolds boundary condition In addition, at the position of full film rupture, the circumferential gradient must be equa to zero. From the above analysis, compared with the experiment results, the error is large   Numerical research on the water-lubricated bearing with cavitation is carried out to explore the influences of propeller deformation on the behaviors of the bearing. For the following situation, the water-lubricated stern bearing (shown in Figure 10) variables used in the computation are listed in Table 3. The weight of the screw propeller is 0.21498 × 10 6 N. According to the validation case, the JFO boundary condition is much closer to the real condition. Therefore, the following analysis is only based on the JFO boundary condition. Meanwhile, the deformations of the bearing and propeller are taken into consideration. The structure of the propeller is shown in Figure 10. The shaft is divided into two parts by the flange. According to the position of the bearings, the whole shaft can be divided into five segments, as shown in Figure 10. The length of each part is listed in Table 3, together with the bearings' parameters and material properties.

The Equivalent Stiffness
The dynamic stability and maximum vibration amplitude of the bearing-rotor system are affected by the cavitation to a large extent, especially for the dynamically loaded bearing. The stability of hydrodynamic bearings is generally addressed since the early design stages of a rotor-bearing system. The water film is recognized as the spring with the equivalent stiffness: Then from [35], And this equivalent stiffness is important in the stability analysis of hydrodynamic bearings. In stability calculation, if K SS < 0, the system is absolutely unstable.
The absolute smooth surface does not exist. Generally, the surface roughness is identical to, or larger than the film thickness estimated by the smooth surface hydrodynamic lubrication theory. Based on the model developed by Patir and Cheng [36], the empirical pressure flow, shear flow factors and asperity contact factor are used in this study. The correlation between relative velocity and coefficient of friction of bearing-rotor system is constructed by Stribeck under constant load to describe the range of different operating conditions. The reference [37] has indicated the relationship between the asperity contact force and lubrication regime, and the bottom of the Stribeck curve is regarded as the transition point of mixed lubrication and hydrodynamic lubrication. Due to the concentrated load of the propeller, the bending deformation is taking place in the whole operating condition. Simultaneously, the journal deformation is affected by the hydrodynamic pressure and asperity contact force [38,39]. Furthermore, for a large width-bearing like the water-lubricated stern bearing in ships (as shown in Figure 10), in the analysis of the bearing-rotor system, the bending deformation of the shaft cannot be ignored. Figure 11 shows the Stribeck curve and the equivalent stiffness K SS in different speeds of bearing 1. From Figure 11, we can see that under hydrodynamic lubrication, the equivalent stiffness K SS decreases with speeds. As noted, the K SS descends sharply under low speeds. However, with the increase in speeds, the trend is inclined to be stable. One reasonable explanation for it is that in hydrodynamic lubrication and the hydrodynamic film forms, it is enough to support the journal away from the bearing. Finally, the solid contact between the journal and bearing disappears, and the equivalent stiffness reduces rapidly and tends to be stable.

The Natural Frequency Analysis
Based on the method developed by Zhou [40], the natural frequency of the shaft is calculated in this study. Due to the dynamic parameters and water film, supporting positions are various with speed. The natural frequency of the shaft should be different at different speeds. The natural frequencies at 90 rpm, 110 rpm, 130 rpm and 150 rpm in different frequency ratios are listed in Table 4. It is clear that the speed has an influence on the natural frequency when lubrication is taken into consideration. For natural frequency in the present, it increases first with speed and then decreases, which means there is a critical speed for the natural frequency. From the working conditions of rotational speeds, the critical speed exists between 130 rpm and 150 rpm. It can be explained by the fact that in addition to external loads and tangential forces, the shaft is also subject to large centrifugal force and gyroscopic moment. Therefore, under so many influence factors, the situation of natural frequency is caused. Compared with the results by reference [40] (the predicted model used in reference [40] only considers the bending deformation of the shaft caused by a concentrated load of screw), the results predicted in this study are smaller. That means the lubrication cannot be neglected in the natural frequency calculation [41,42].

Conclusions
The dynamic parameters of the stern bearing considering cavitation, as well as the bending deformation of the propeller, are investigated. Based on the displacement superposition method, the interaction between lubrication properties and bending de-

The Natural Frequency Analysis
Based on the method developed by Zhou [40], the natural frequency of the shaft is calculated in this study. Due to the dynamic parameters and water film, supporting positions are various with speed. The natural frequency of the shaft should be different at different speeds. The natural frequencies at 90 rpm, 110 rpm, 130 rpm and 150 rpm in different frequency ratios are listed in Table 4. It is clear that the speed has an influence on the natural frequency when lubrication is taken into consideration. For natural frequency in the present, it increases first with speed and then decreases, which means there is a critical speed for the natural frequency. From the working conditions of rotational speeds, the critical speed exists between 130 rpm and 150 rpm. It can be explained by the fact that in addition to external loads and tangential forces, the shaft is also subject to large centrifugal force and gyroscopic moment. Therefore, under so many influence factors, the situation of natural frequency is caused. Compared with the results by reference [40] (the predicted model used in reference [40] only considers the bending deformation of the shaft caused by a concentrated load of screw), the results predicted in this study are smaller. That means the lubrication cannot be neglected in the natural frequency calculation [41,42].

Conclusions
The dynamic parameters of the stern bearing considering cavitation, as well as the bending deformation of the propeller, are investigated. Based on the displacement superposition method, the interaction between lubrication properties and bending de-formation has been taken into consideration. By implementing the validated lubrication model designed in this study, the dynamic parameters of bearing systems are presented by the finite perturbation method. As an important parameter in stability analysis, the property of equivalent stiffness in the hydrodynamic lubrication regime is given. Finally, some significant conclusions are obtained, as follows: (1) For equivalent stiffness, due to the increase in hydrodynamic effect, it is more affected by the speed, especially at low speeds; (2) For natural frequency, there is a critical speed between 130 rpm and 150 rpm, which makes the natural frequency strike the maximum value because of the comprehensive influencing factors (external loads, tangential forces, large centrifugal forces and gyroscopic moment).