A Fluid–Solid–Heat Coupling Analysis for Water-Lubricated Rubber Stern Bearing Considering the Deﬂection of Propeller Shaft

: The novelty of the present study is that it investigates the effect of propeller shaft deﬂection, caused by the propeller self-weight and interfacial mixed forces, on the ﬂuid–solid–heat (FSH) coupling performance of water lubricated rubber stern bearing (WLRSB). In the FSH coupling model, the generalized average Reynolds equation and the Kogut–Etsion asperity contact model are used to determine the hydrodynamic and the elastic–plastic contact behaviors of WLRSB. In the thermal analysis, the journal, water ﬁlm, and rubber bushing are considered as an integrated system (JWR system) using the Euler method. To prove the correctness of the developed model, the predicted results are veriﬁed by comparisons with the experimental results given in the literature. In addition, to assess the effect of the force-driven deﬂection during FSH simulation, comparisons of the FSH predictions between the aligned journal case and the deﬂected journal case are carried out. The results indicate that, especially under a heavy load condition, the deﬂection of the stern shaft should be incorporated into the lubrication gap between the journal–rubber interface during the analysis of FSH performance of the JWR system.


Introduction
Water-lubricated rubber stern bearings, which are widely used in marine propulsion systems, have been attracting increasing attention owing to their excellent performance [1], such as low temperature rise, environmental friendliness, convenient assembly, etc. Figure 1 depicts a typical stern bearing lubricated with water, which exhibits a lower rate of viscous dissipation and a better cooling capacity compared with an oil-lubricated journal bearing [2]. However, the viscosity of water is extremely lower than that of mineral oil, implying that the hydrodynamic lifting generated by a water film is typically not able to support the external load [3][4][5], and therefore, asperity contact occurs at the journal-rubber interface [1,6], which affects the mixed friction behaviors by the generated frictional heat and thermal expansion. Therefore, the fluid field (from the hydrodynamic water film), solid field (from the asperity contact and elastic-thermal deformation), and thermal field (from the mixed friction heat) coexist (FSH coupling effect) at the journal-rubber interface. Consequently, a multi-field coupling numerical analysis of the water-lubricated stern bearings is necessary to investigate their reliability and failure mechanism. However, comprehensive theoretical works related to this topic are still limited. Studies have been published to investigate the fluid-solid coupling performances (or mixed friction performances) of stern bearings [7][8][9][10][11][12][13]. He et al. [8] and Lv et al. [12] presented a mixed lubrication model for a misaligned marine stern tube bearing, in which the effects of cavitation, bend deformation of the propeller shaft, and the wall slip of thin film have been considered. Kraker et al. [9] investigated the Stribeck curves of waterlubricated bearings experimentally under a wide range of operating conditions. To improve the computation efficiency of the fluid-solid coupling model of journal bearings, a fast computing technology was developed by Han et al. [11] and Chan et al. [10]. However, the studies mentioned above ignored the thermal effects caused by the asperity contact behavior between the journal and the bush.
The thermal effects caused by fluid film shearing significantly affect the performance of a stern bearing and have been studied numerically and experimentally by scholars. The experimental results of the oil film temperature distribution presented by Majumdar et al. [14], Bouyer et al. [15], and Pierre et al. [16] demonstrated that the maximum temperature of oil film occurred near the maximum fluid pressure location, and that oil film temperature distribution was highly dependent on operating conditions such as the journal speed, external load, journal misalignment, etc. It is noteworthy that a hydrodynamic bearing system can be considered as a heat conduction coupling system during thermal analysis. Related studies were experimentally performed by Mitsui et al. [17], Fitzgerald et al. [18], and Wang et al. [19]. In the field of theoretical modeling, thermohydrodynamic analyses for journal bearings were performed by Yadav et al. [20], Sahu et al. [21], Sun et al. [22], Mishra et al. [23], Lin et al. [24], Hoon et al. [25], Pu et al. [26], and Li et al. [27]. For the thermal study of water-lubricated bearings, Feng et al. [28] numerically investigated the thermohydrodynamic performance of a water-lubricated bearing considering misaligned and turbulent effects. In their study, the asperity contact effect was absent and the misaligned effect was considered based on a predetermined misalignment value. More recently, Wodtke and Litwin [29] performed an interesting study to investigate, numerically and experimentally, the effects of forced axial flow on the thermal behaviors of multi-groove water-lubricated bearings, in which the journal was treated as aligned. These theoretical studies primarily focused on the thermal performance of the fluid film and assumed that the heat was only generated by the viscous shear of the fluid film. However, during operation, the asperity contact friction heat caused by the contact effect at the journal bearing interface affects the performance significantly. The aforementioned studies fail to evaluate the thermal behavior originating from the asperity contact effect. To predict the thermal performance of WLRSBs more accurately, in the present study, the heat generated by the contact effect has been integrated in the thermal analysis. Subsequently, a general transfer equation based on the Euler method was used for simulating the thermal behavior of the journal, water film, and rubber bushing (JWR) system.
It is noteworthy that misalignment caused by the deflection of the propeller shaft has been observed to be a familiar phenomenon in a stern bearing system because WLRSBs exhibit a larger length-to-diameter ratio compared with general journal bearings. Recently, the effects of journal misalignment on the performance of the journal bearing interface have been investigated by many scholars (He T  In these studies, He T et al. [8] and Sun J et al. [22] calculated the misalignment using the external load acting on the journal and concluded that the tribology performance of the bearing was sensitive to misalignment of the journal, especially under a large eccentricity ratio. Mallya et al. analyzed the effects of shaft misalignment on the stability of journal bearings, and the thermal effects due to viscous shearing of oil film were also included in Mallya's study. Han et al. [11], Lv et al. [12], Litwin et al. [30], and Mallya et al. [31] investigated the tribology performance of a misaligned bearing numerically under given misalignment parameters, i.e., degree of misalignment and misalignment direction angle. In the current study, in order to obtain the misalignment of the propeller shaft more accurately, the external force applied to the end face of the WLRSB and the hydrodynamic and contact forces were combined to determine the deflection of the propeller shaft under the FSH coupling conditions. The purpose of this study was to investigate the FSH coupling performance of a WLRSB when the force-driven deflected shaft is taken into account. In order to evaluate the deflection of the shaft more accurately, the segmented mixed forces along the axial direction (i.e., hydrodynamic and contact forces) serve as the input information in the deformation analysis of the shaft, and the superposition method is applied to attain the deflection of the propeller shaft. Subsequently, the FSH coupling model is coupled with the deflection model to identify the deflected shaft effect on the FSH predictions. The proposed model is validated by comparing it with published experimental results. Therefore, this study can provide a more practical model to predict the FSH coupling performance for WLRSBs.

Governing Equation
The water film pressure of a WLRSB considering surface roughness effects and 3D viscosity distribution can be determined by the 3D generalized average Reynolds equation, which can be expressed as: where φ θ and φ z [4] are the pressure flow factor in θ and z directions, respectively; φ c [32] and φ s [4] are the contact factor and shear factor, respectively; p h is the water film pressure, ρ W is the density of water, and η is the viscosity of water. Furthermore, ω is the speed of the journal, h is the average film thickness, and σ is the composite standard deviation of the rubber surface. In Equation (1), the parameters I R , J R , and G R can be determined by: where y indicates the radial coordinate along the cross-film direction. The profile of the WLRSB is shown in Figure 2. As shown, the lubrication clearance can be evaluated as the sum of the average geometric gap considering the deflection of the propeller shaft, h 0 , the groove depth, h g , and the elastic/thermal deformation of the shaft and the bush, δ J and δ R . Consequently, the water film thickness, h(θ, z, ∆T), is given by the following: where z represents the axial direction shown in Figure 2, ∆T is the temperature rise at the coordinate of (θ, z), and h 0 can be calculated as: where ε i is the eccentricity ratio of the ith cross-section of the bearings, and C is the radial clearance of the bearing.  In the force analysis of the propeller shaft, a simplification has been made so that the propeller shaft in the bearing hole is treated as a cantilever beam, which can be supported by two aspects. The first one is that the length-to-diameter ratio of the WLRSB is greatly larger than that of the general journal bearing, the second one is that the external load, i.e., the self-weight of the propeller, locates at the end face of the water-lubricated stern bearing. Additionally, a similar simplification was made in the study performed by He et al. [8]. As shown, the bearing is discretized using n girds in the axial direction, the shaft is simplified as a cantilever beam, and the end face of the bearing that is close to the propeller is assumed to be fixed. The force acting on the ith section of the shaft is defined as F i , which is equal to the vector sum of F iz in the z direction and F ix in the x direction. It is noteworthy that both  Figure 3. According to the superposition principle, d ix and d iy can be calculated using the following formula: where M is the self-weight of the propeller, L is the width of the bearing, E J is the elastic modulus of the shaft, I is the inertia moment of the shaft section, and I = πD 4 /64, where D is the diameter of the shaft.
The eccentricity, e i , between the point O i and the point O b , as shown in Figure 3, can be calculated as: where the distance d i can be calculated by Further, e s is the eccentricity of the end face of the shaft. Subsequently, the ε i in Equation (4) can be written as: The average geometric gap, h 0 , can be obtained by substituting Equations (5)-(7) into Equation (4). When asperity contact occurs, the total pressure, p, acting on any nodal (θ j , z k ) of the stern bearing inner face consists of two parts: the water film pressure, p h , and the contact pressure, p c . The latter can be determined using the Kogut-Etsion elastic-plastic contact model [33], which is given by: where the superscript "*" denotes dimensionless parameters, H B is the hardness of the bush, whose unit is Pa; β and D are the curvature radius and the density of the asperity, respectively; and I c is the integral operator, which is calculated by: where φ * (y * ) can be expressed as: where σ s is the surface roughness of the shaft, and the value of σ s /σ is set to 1.0 in the present simulation. In Equation (8), the dimensionless parameters ω * c and K can be calculated by the following expressions: where E c is the composite elastic modulus and ν B is the Poisson ratio of the bush. The deformation of the bearing, δ R , can be expressed as the sum of the elastic deformation, δ RE , and the thermal deformation, δ RT . The elastic deformation of the journal, δ JE , can be ignored because its elastic modulus is significantly larger than that of the rubber bushing. Additionally, the thermal deformation of the journal, δ JT , is assumed to be a linear thermal expansion so as to simplify the analysis. Consequently, the deformation equations of the journal and the rubber bushing are given as follows: where the node (θ ξ , z η , r ζ ) is shown in Figure 4; ∆T(θ ξ , z η , r ζ ) is the temperature rise of the rubber bushing at the nodal of (θ ξ , z η , r ζ ), and G RE (θ j , y k , θ ξ , y η ) and G RT (θ j , z k , θ ξ , z η , r ζ ) represent the influence functions of the elastic and thermal deformations, respectively. The calculated method of the former is presented by Chan et al. [10], and that of the latter can be referred to in Ref. [34]. In the thermal analysis, the journal, water film, and rubber bushing are treated as a heat conduction coupling system (JWR system for simplicity), as shown in Figure 4, using the Euler method. Based on this, the thermal behavior of the water film, propeller shaft, and rubber bushing can be described by a general transfer equation in polar coordinates, which is given by: where ρ is the density of the rubber bushing (ρ R ), journal (ρ J ), or lubricant (ρ W ); C P is the specific heat of the bushing (C PR ), journal (C PJ ), or water (C PW ); k r (k rR , k rJ , k rW ), k θ (k θR , k θJ , k θW ), and k z (k zR , k zJ , k zW ) denote the thermal conductivities in r, θ, and z directions, respectively; v r , v θ , and v z denote the velocity at which the rotation speed corresponds to the journal and the flow velocity corresponds to the water in r, θ, and z directions, respectively. One study [29] pointed out that the thermal behaviors of water-lubricated bearings depend significantly on the intensity of the enforced axial flow. In the present study, we assume the axial flow is sufficiently low, and therefore, the effect of the axial flow on the thermal behavior of the water-lubricated bearing is neglected in the current simulation.
Subsequently, the flow velocity in the θ direction, v Wθ , and in the z direction, v Wz , is given as follows: v Wθ = u + uφ s where u is the linear velocity of the journal, u = ωR 1 . Furthermore, the parameter z m can be evaluated by The heat source in the JWR system is composed of the viscous shear heat of the water film, Φ W , and the friction heat, Φ C , generated from the asperity contact occurring at the journal-rubber interface, which can be described by the following formulas: (15) where µ c denotes the friction coefficient of the asperity contact, and its value was assumed to be 0.15 in this study. Subsequently, the heat source, Φ, can be obtained as follows: In Equation (16), dV is the differential volume for the water film, and its corresponding area at the journal-rubber interface is represented as dA. Consequently, the thermal properties of the JWR system can be determined by coupling Equations (13)- (16). Similar to the evaluation of the thermal deformation, the influence function method was also employed to obtain the temperature distribution of the JWR system. The corresponding influence functions can be referred to in Ref. [34].
The effect of the temperature rise on the water viscosity can be evaluated by: where η(θ, z, r) and T W (θ, z, r) are the 3D average viscosity and 3D average temperature of the water film along the radial direction, respectively. It should be noted that the unit of the temperature is K in Equation (17).
Once the hydrodynamic and asperity contact pressures are obtained, the external force acting on the propeller shaft can be calculated as follows: where F x and F y represent the force, including hydrodynamic and contact forces, in the x and y directions, respectively. Subsequently, the load capacity can be described by which is equal to the external force M. The friction force, F f , considering the water film shearing in the circumferential direction and the asperity contact at the journal-rubber interface can be estimated by: where ω is the rotation speed of the propeller shaft. Subsequently, the average friction coefficient, µ, that combines the effects of viscous shearing and asperity contact is given by:

Boundary Condition
In this study, the following fluid pressure boundary condition is applied to solve the average Reynolds equation: where θ 0 is the interface between the pressured region and the cavitation region.
For simulating the elastic deformation, the outer surface of the bearing bush, which connects with the bearing house shown in Figure 1, is assumed to be fixed, as shown in Figure 5a, while all of the other boundaries in the calculation of elastic deformation are considered as elastically free. The thermal boundary conditions for the JWR system, shown in Figure 5b, can be classified into two categories. (1) Internal boundary conditions: B.C.1 and B.C.2 represent the bearing-water and waterjournal interfaces; which can be evaluated by the following Equation (22).
It should be mentioned that the journal is treated as isothermal along the circumferential direction.
(2) External boundary conditions: B.C.3 through B.C.7 are the external boundaries for the system and they must be specified. In the present thermal model, the convective boundary is applied to describe the thermal behavior at B.C.3-B.C.7, which is given by: where T c and T ∞ represent the water inlet and ambient temperatures, respectively, and h h denotes the convective heat transfer coefficient.

Numerical Algorithm
The discrete form of the average Reynolds equation (Equation (1)) can be obtained by finite difference method and the obtained linear equations at each unknown pressure nodes are given by: The successive overrelaxation (SOR) method was used for solving Equation (24), which is described by: /A j,k , and λ is the relaxation factor-its value is assumed to be 1.5 in the current simulation.
The flow chart used to simulate the fluid-solid-heat coupled performance of the JWR system is illustrated in Figure 6, which includes a procedure for the calculation of fluid pressure, contact pressure, thermal-elastic deformation, temperature distribution, and shaft deflection. In the present study, a calculation procedure based on Fortran is developed to determine the FSH coupling performance of the water-lubricated rubber bearing considering propeller shaft deflection. After numerical testing at a computer with a 1.61 GHz CPU and 16 GB RAM, it was observed that the evaluation time for the bearing properties is around 30 min. Furthermore, it is worth mentioning that for the statistical hydrodynamic and contact models, only the average Reynolds equation was calculated for the hydrodynamic lubrication regime; however, both the average Reynolds equation and asperity contact model were calculated for the mixed friction regime.
The pressure, including the water film pressure and contact pressure, is controlled by the error given in Equation (26). The temperature and attitude angle are controlled by the error given in Equations (28) and (29), respectively.
where T max = (T Jmax +T Wmax +T Rmax )/3, with T Jmax , T Wmax , and T Rmax being the maximum temperature of the journal, water film, and rubber bushing, respectively.
The attitude angle can be adjusted by the force component in x and y directions, which is given by:

Model Validation
In order to validate the FSH coupling model, the journal bearing parameters given by Wang et al. [19] were used to obtain the temperature distributions under different external loads. As shown in Figure 7, a comparison of the temperature distribution between the present model and the experimental data provided by Wang et al. [19] was performed. A good agreement can be found between the predicted results and the experiment data. Verification of the thermal distribution of the bearing-lubricant-journal system was carried out by comparison with Mitsui's experiment results (Mitsui et al. [17]), as shown in Figure 8, and the geometry, material, and lubricant parameters used in the presented model are consistent with those from Mitsui's experiment. It was observed that the maximum temperature predicted by the present model agrees well with that tested by Mitsui, and the distribution trends of the temperature obtained by the developed model agree well with those given in Mitsui's experiment. Therefore, the comparisons of the temperature distribution support the correctness of the developed thermal model.  The proposed model was further validated by comparing the predicted load capacity with the experiment results from Hirs et al. [35]. Similarly, the geometry, material, and lubricant parameters used in the simulations were all the same as Hirs' parameters. As shown in Figure 9, the load capacities under various eccentricity ratios agree well with those given in Hirs' experiment [35]. Consequently, the comparisons above demonstrate the correctness of the present model, and therefore, it can be used to analyze the FSH coupling performance of the WLRSB. In this study, the basic parameters, including the dimensions, materials, thermal parameters, and surface texture of the JWR system, used in the proposed model are listed in Table 1.  Figure 10 shows the distribution of the water film pressure, contact pressure, elastic deformation, and thermal deformation under the deflected propeller shaft. Figure 10a shows that the deflected shaft causes asymmetric hydrodynamic pressure distribution. By comparing the results plotted in Figure 10a, it was found that the external force is mainly carried by the asperity contact in the current simulated case, which demonstrates that the asperity contact plays a crucial role in the actual operating condition of the WLRSB. The results shown in Figure 10 indicate that the maximum fluid pressure, contact pressure, elastic deformation, and thermal deformation are mainly located on the water outlet side when the deflection of the propeller shaft is considered. This is because the deflected shaft decreases the lubrication gap close to the water inlet, thereby resulting in an increase in hydrodynamic pressure, asperity contact, and frictional heat. Figure 10c shows that negative elastic deformation occurs near the maximum elastic deformation; this is due to the fact that the rubber can be treated as an incompressible material owing to its Poisson's ratio being near to 0.5. Moreover, Figure 10d indicates that the thermal deformation of the groove region is less than that of the pressured region; this is because the water flow in the groove removes a large amount of heat, thus reducing the thermal expansion in such a region. In general, it can be concluded that strong hydrodynamic and contact effects are yielded at the outlet water side of the JWR system when deflection of the propeller shaft is considered.  Figure 11 compares the circumferential, axial, and radial temperature distributions yielded by the deflected shaft case and the aligned shaft case. Figure 11a illustrates the circumferential temperature distribution of the JWR system at the section containing the maximum water film temperature node. As shown, the temperature of the groove region is obviously smaller than that of the pressured region because the water flow in the groove region can facilitate heat dissipation. Furthermore, the distribution trend simulated by the deflected shaft case is similar to that simulated by the aligned shaft case; however, the maximum temperature predicted by the deflected shaft case is obviously smaller than that predicted by the aligned shaft case. More specifically, the maximum temperature of the misaligned shaft case is about 11.5 • C higher than that of the aligned shaft case.

Thermal Analysis of the JWR System
The axial temperature distribution of the water film is shown in Figure 11b, and the radial temperature distribution of the JWR system at the maximum temperature section is shown in Figure 11c. As shown in Figure 11b, compared with the aligned shaft case, the maximum temperature zone generated by the misaligned shaft case is closer to the water outlet side. The results demonstrate that inaccurate FSH predictions may be yielded if the deflection of the propeller shaft is ignored. As shown in Figure 11c, the shaft exhibits a lower temperature compared with the water film and the rubber bushing owing to the higher thermal conductivity coefficient (Table 1) compared with the other two. It is noteworthy that, as shown in Figure 11d, the temperature of the misaligned shaft case is lower than that of the aligned shaft case for z/L in the range from 0 to 0.58 and that the maximum water temperature is located nearby the water outlet for both aligned and misaligned shaft cases. These discrepant results between the two cases prove that the deformation of the propeller shaft, resulting from the propeller self-weight, the hydrodynamic force of water film, and contact force of the journal-rubber interface, cannot be neglected in the thermal analysis of WLRSBs.
The results shown in Figure 12 are the temperature distributions in the circumferential direction for three different axial sections: the maximum temperature section, the middle section (z/L = 0.5), and the end face of the bearing (z/L = 1). As shown for the case of the z/L = 0.5 section, the maximum water film temperature of the aligned shaft case is slightly higher than that of the misaligned shaft case. However, Figure 12b,c show that for the cases of the maximum temperature section and the z/L = 1.0 section, the water film temperature yielded by the misaligned shaft case is obviously higher than that yielded by the aligned shaft case. These results indicate that the thermal performance of the WLRSB depends highly on the fluid-solid coupling effects of the outlet water side when the deflection of the shaft has been considered.  Figure 11. Comparisons of the thermal distribution in the circumferential section, axial direction, and radial section obtained with a misaligned shaft and an aligned shaft: (a-c) are the temperature distribution in the circumferential section, axial section, and radial section, respectively: (d) is the temperature distribution at the maximum water temperature section in the axial direction.

Comparisons of FSH Coupling Performance
In this section, the FSH coupling performances of the JWR system, including the contact load, maximum temperature, and maximum thermal-elastic deformation, obtained with the aligned shaft case and the misaligned shaft case under different external forces are studied comparatively. In order to make the sensitive studies clearer for readers, under various parameters, a dimensionless process for pressure (hydrodynamic and contact) and deformation (elastic and thermal) is performed in the figures of the current section. The dimensionless process of pressure and load can be expressed as: where N is the revolutions of the journal per second. The dimensionless process of deformation can be expressed as: The comparisons of the contact load and the friction coefficient between the two cases are shown in Figure 13a. The contact load yielded by the misaligned shaft case was almost same as that from the aligned shaft case when the external load was less than 3 kN. However, when the external load increases from 3 to 8 kN, the contact load generated by the misaligned shaft case becomes higher than that of the aligned shaft case; the reason for this is that the deflected shaft effect is gradually enhanced with the increasing external load. Additionally, it can be found from Figure 13a that the friction coefficient for both cases decreased when the load capacity increased from 1 to 3 kN, which conforms to the elastohydrodynamic lubrication theory [4]. Figure 13b,c show comparisons of the maximum temperature and thermal-elastic deformation, respectively. The maximum temperatures of the misaligned shaft case and the aligned shaft case were almost the same when the external load was less than 4 kN, implying that the deflection of the propeller shaft almost does not affect the thermal performance of the WLRSB under a relatively light load condition. Meanwhile, the maximum water film temperature of the misaligned shaft case became higher than that of aligned shaft case when the load capacity was larger than 4 kN. Furthermore, Figure 14a,b can be used to explain the relationship of the maximum water film temperature between the two shaft cases. As shown in Figure 14b, for the misaligned shaft case, the contact pressure occurred mainly in the outlet water side, and its value increased sharply with the increase in load capacity. It can also be found from Figure 14b that for the load capacity in the range from 4 to 10 kN, the load capacity yielded by the misaligned shaft case was obviously larger than that yielded by the aligned shaft case (see Figure 14a). The results indicate that the misaligned shaft generates more friction heat, Φ C , compared with the aligned shaft. However, the misaligned shaft provides a smaller average temperature of journal than that of aligned shaft case when the load capacity is larger than 4 kN. This may be due to the fact that the aligned shaft yields a larger contact area under heavy load conditions. As expected, the misaligned shaft gives a larger elastic-thermal deformation than the aligned shaft, as shown in Figure 13c.  Additionally, it can be found from Figure 14a,b that "Edge contact" occurs for both shaft cases, and this phenomenon has also been observed from Litwin's experiment (Litwin et al. [6]). Figure 14c illustrates the mechanism of the "Edge contact" shown in Figure 14a,b. As shown for the aligned shaft case, the deformation in the central area is larger than that on both sides of the WLRSB due to the "water bag" caused by hydrodynamic effects (Zhou et al. [36]), implying that the contact effect will occur on the inlet and outlet sides of the WLRSB. However, the "Edge contact" for the aligned shaft case disappears when the load capacity is larger than 4 kN, as shown in Figure 14a, which may be due to thermal expansion leading to a contact effect in the central area. As shown in Figure 14c, for the misaligned shaft case, a "water bag" may occur on the divergent side (z/L: 0-0.5) owing to this region providing a relatively large clearance to contain hydrodynamic water. The results shown in Figure 14b also indicate that for the misaligned shaft case, the "Edge contact" increased with the increase in external load.

Conclusions
A fluid-solid-heat coupling model for WLRSBs considering the deflection of the propeller shaft was built, and the proposed model was used to investigate the fluid-solid-heat performance of the JWR system. Based on the Euler method, a general transfer equation in polar coordinates was employed to describe the heat convection in a water film and the heat conduction in rubber bushing and a journal. The segmented mixed forces, including the hydrodynamic and contact forces, and the superposition method were bridged to evaluate the deflection of the stern shaft. Additionally, the FSH coupling performances obtained with two propeller shaft cases, i.e., aligned shaft case and misaligned shaft case, were compared systematically. The conclusions are as follows: (1) The temperature distribution and the load capacity obtained by the proposed model agree well with the experimental data given in the literature, thus validating that the presented model can be used to predict the FSH coupling performances of hydrodynamic bearings. (2) The asymmetric tribological performance distributions along the axial direction occur due to the propeller shaft deflection, including the hydrodynamic force, asperity contact, and temperature. The maximum asperity contact and temperature rise of the misaligned shaft are obviously larger than those of the aligned shaft. Furthermore, the temperature yielded by the water outlet is larger than that yielded by the water inlet for both the misaligned and aligned shaft cases. (3) Owing to the WLRSB having a larger length-to-diameter ratio compared with the general journal bearing, the misalignment caused by the deflection of the propeller shaft affects the FSH performance of the JWR system significantly, especially under a relatively heavy load condition. Furthermore, the aligned shaft case was found to be a proper substitute for a misaligned shaft case when the load capacity is smaller than 4 kN for the current simulation.
This study develops a comprehensive steady FSH coupling model for a water-lubricated rubber bearing operated at the stern shaft system, in which shaft deflection is caused by the self-weight of the propeller that is integrated. The present mathematical model conducts toward a guideline to help design appropriate parameters for increasing the reliability of water-lubricated rubber bearings. In addition, future works are needed to integrate rotor dynamic characteristics into the present model and to study the interaction between the FSH coupling performance and the high inertia components of a propeller in specific applications.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.