Issue 
Mechanics & Industry
Volume 24, 2023



Article Number  3  
Number of page(s)  15  
DOI  https://doi.org/10.1051/meca/2022029  
Published online  12 January 2023 
Regular Article
Study on elastohydrodynamic lubrication performance of doublelayer composite waterlubricated bearings
^{1}
School of Mechanical Engineering, Shandong University of Technology, Zibo, China
^{2}
Shandong Provincial Key Laboratory of Precision Manufacturing and Nontraditional Machining, Shandong University of Technology, Zibo, China
^{3}
Marine Design and Research Institute of China, Shanghai, China
^{*} email: zhangxiulli@163.com
Received:
12
August
2022
Accepted:
12
December
2022
Doublelayer composite waterlubricated bearing is a new type of waterlubricated bearing which can integrate the good damping performance of low elastic underlayer bush and good tribological performance of plastic layer bush. This paper analyzes its elastohydrodynamic lubrication performance by fluid–structure interaction (FSI) method, and studies the effects of eccentricity ratio, rotational speed, elastic modulus distribution and thickness distribution of bearing bush on its lubrication performance. Results show that the lubrication performance of doublelayer bearing is more like that of plastic bearing. As rotational speed and eccentricity ratio increase, the maximum water film pressure, the load carrying capacity and the maximum bush deformation increase significantly. As the elastic modulus of the low elastic underlayer bush decreases, the total bush deformation increases significantly, but the load carrying capacity decreases slightly. The bush thickness distribution influences the deformation distribution of both low elastic underlayer bush and plastic layer bush, but have little impact on the total bush deformation and bearing lubrication performance.
Key words: Doublelayer bush / water lubrication / journal bearing / elastohydrodynamic lubrication / fluid–structure interaction
© Z. Yang et al., Published by EDP Sciences 2023
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Waterlubricated bearings have the characteristics of environmental protection and energy conservation, and are now widely used in ship propulsion system, nuclear main pump, water turbine and other hydraulic machinery. At present, stern tube waterlubricated bearings mainly use rubber as bush material, because rubber has good damping performance and can weaken the adverse effect of journal misalignment and vibration. But rubber bearings have large bush deformation and low hydrodynamic load carrying capacity, and the large friction coefficient usually causes high power consumption and noise [1,2]. In order to improve the bearing load carrying capacity and reduce power consumption, some turbine and marine pump waterlubricated bearings use plastic such as polytetrafluoroethylene (PTFE) and ultrahigh molecular weight polyethylene (UHMWPE) as bush material [3–6]. While plastics have low friction coefficient, but worse compliance and damping properties than rubber. In recent years, some research institutions and researchers have studied doublelayer composite waterlubricated bearings. Duramax Marine has developed a doublelayer stern tube bearing ROMOR I with nitrile rubber (NBR) as one layer bush and UHMWPE as the other, which has been applied in US surface ships and submarines [7]. Litwin [8–10] studied the lubrication performance of a threelayer waterlubricated bearing (bronzeNBRPTFE). Results show that the bearing works under hydrodynamic lubrication state for most conditions, and its running resistance is less than the NBR bearing.
As a new kind of waterlubricated bearing, the doublelayer composite bearing has not been studied extensively. But many scholars have studied the lubrication characteristics of waterlubricating rubber bearings or plastic bearings. Cabrera et al. [11] studied the film pressure distribution in waterlubricated rubber journal bearings by experiments, and found that the bearings operate in the regime of mixed lubrication. Due to the elastic deformation of rubber, the film pressure profiles are obviously different from those of rigid bearings, and peak pressures are greatly reduced. Kraker et al. [12] calculated the Stribeck curves for waterlubricated journal bearings by a mixed elastohydrodynamic lubrication model, and found that within a certain range, the smaller the radial clearance, surface roughness, external load and elastic modulus of the bearing, the smaller the rotational speed at the mixed lubrication (ML) to elastohydrodynamic lubrication (EHL) transition. Liao et al. [13] studied the effects of different lubricants (pure water, saline and sandy water) on the lubricating properties of waterlubricated rubber bearings. Results showed that the water film pressure of the bearing generally increases and stabilizes with increasing rotational speed and eccentricity, but sandy water shows a greater increase compared with saline and pure water. The deformation and the stress of the carrying surface shows the same trend. Wang et al. [14] studied the lubrication characteristics of waterlubricated bearing by CFD and FSI method. Results showed that the bearing bush deformation reduces the water film pressure and increases the film thickness, which reduces the bush stress concentration. Lao et al. [15] modified UHMWPE by adding Polyimide (PI) and studied its tribological performance in artificial seawater. They found that the modified material has the lowest friction coefficient and wear height than Thordon, Tenmat, and PEEK. Zhou et al. [16,17] studied the friction coefficient, water film pressure and shaft orbits of multiple groove rubber alloy bearing, and compared the Stribeck curve of straight groove and spiral groove bearing. Results showed that groove and bearing bush deformation have a great influence on the pressure distribution of bearing, and spiral groove bearing has better lubrication performance than straight groove bearing. Zhao and Rui [18] studied the hydrodynamic characteristics of a waterlubricated rubber journal bearing with a large aspect ratio by finite difference numerical method and experiments. They found that the liquid film was not continuous at low speeds. Wang et al. [19] studied the tribological properties of the rubber bearings under different loads and velocities. The experimental results showed that the friction coefficient decreases with increasing velocity, the friction coefficient increases at first and then decreases with increasing load. Huang et al. [20] studied the dynamic performance of a waterlubrication rubber bearing with rubber damping layer, and found that the damping layer can reduce bearing vibration, especially under high rotational speed and high load conditions. Xiang et al. [21] explored nonlinear dynamic mixed friction behaviors of Thorlon waterlubricated bearings.
Although the doublelayer composite waterlubricated bearing has been confirmed to have superior properties, due to the large difference in the elastic modulus of low elastic underlayer bush and plastic layer bush, the bush deformation of this kind of bearing is quite different from those of waterlubricated rubber bearing and plastic bearing, and its lubrication performance is also different. Therefore, this paper built a fluid–structure interaction model of doublelayer composite waterlubricated bearing to analyze its elastohydrodynamic lubrication performance, and studied the influence of eccentricity ratio, rotational speed, elastic modulus distribution and thickness distribution of the doublelayer bush on the lubrication performance of this kind of bearing, so as to provide reference for the design of such bearings.
2 Bearing structure and parameters
Figure 1 shows the structure of the doublelayer composite waterlubricated bearing. The low elastic underlayer bush is between the plastic layer bush and stainlesssteel sleeve. There is radial clearance, c, between bearing and journal. Under external load, when the journal rotates, the journal center O_{j} is away from the bearing center O_{b}, forming wedgeshaped lubrication film and generating dynamic pressure and load carrying capacity to support the external load. Parameters affecting the bearing load carrying capacity include journal rotational speed N, bearing inner diameter D, bearing length L, plastic bush layer thickness t_{1}, plastic layer bush elastic modulus E_{1} and Poisson's ratio ν_{1}, low elastic underlayer bush thickness t_{2}, low elastic underlayer bush elastic modulus E_{2} and Poisson's ratio ν_{2}, radial clearance c, journal eccentricity ratio ε (=e/c), viscosity of lubricant µ, etc. eθz is the coordinate system.
Fig. 1 Waterlubricated doublelayer bearing. 
3 Simulation analysis method
3.1 Governing equations
The basic equations for solving the fluid model are mass conservation equation, momentum conservation equation and energy conservation equation. As the temperature rise of waterlubricated bearing is very small, the temperature change is ignored and the energy conservation equation is not considered in the calculation, ⊗ the fluid inside bearings are taken as incompressible.
The mass conservation equation (the continuous equation) is expressed as:
where ρ is the fluid density, v is the fluid velocity vector, .
The momentum conservation equation (NavierStokes equation) is given as:
where p is the static pressure, τ is the stress tensor, ρg and F are the gravitational force and external body force. Stress tensor τ is defined as:
where µ is the fluid viscosity.
When the pressure is reduced to a certain value, the water will evaporate, so the cavitation model needs to be coupled with the above equations to calculate the pressure distribution. This paper uses the ZwartGerberBelamri cavitation model in ANSYS FLUENT. The mass transfer of liquid and gas (evaporation and condensation) is governed by the vapor transport equation [22]:
where α_{v} is the vapor volume fraction, ρ_{v} is the vapor density, R_{g} and R_{c} are the mass transfer source terms connected to the growth and collapse of the vapor bubbles. ZwartGerberBelamri cavitation model assumes that all vapor bubbles in the system have the same size, and R_{g} and R_{c} are defined as [23]
where F_{evap} is the evaporation coefficient, F_{cond} is the condensation coefficient, R_{b} is the vapor bubble radius, α_{nuc} is the nucleation site volume fraction, ρ_{l} is the liquid density, p_{v} is the vapor pressure. In this work, F_{evap} = 50, F_{cond} = 0.01, R_{b} = 10^{−6} m, α_{nuc} = 5 × 10^{−4}.
The load carrying capacity can be calculated by integrating the pressure over the journal surface as follows:
The total load carrying capacity of bearing is defined as:
The attitude angle is
The motion of solid domain is governed by Newton's second law, which can be described as:
where ρ_{s} is the solid density, a_{s} is the local acceleration vector of the solid domain, σ_{s} is the solid stress tensor, and f_{s} is the body force vector.
The displacement d, stress τ, heat flow q, temperature T and other variables of the fluid and solid models are equal on the fluid–solid interface. Ignoring the temperature change, the equilibrium equation can be expressed as:
where index f refers to fluid, index s refers to solid.
Twoway FSI analyses were carried out in this work. Figure 2 shows the flowchart of calculation process. The fluid control equations and solid control equations are solved in different solvers respectively, and the calculation results of fluid domain and solid domain are exchanged through the fluid–structure coupling interface. In this work, the water film pressure distribution was solved in FLUENT firstly, and then transferred to Transient Structure to solve the bearing bush deformation. Then the bush deformation distribution was transferred to FLUENT to update the film thickness distribution and calculate the water film pressure. The time step size was 0.02s and the number of time steps was 50. The maximum number of iterations in FLUENT was 20, and that of the coupling interaction was 5. When the calculation is converged or reaches the number of iterations at this time step, it goes to the next time step. The lift forces are monitored and when the forces remain unchanged, the calculation will be terminated to reduce the processing time.
Fig. 2 Flowchart of calculation process. 
3.2 Models
Figure 3a shows the model of composite bush and its boundary condition. The inner surface of the plastic bush is the FSI interface, and the outer surface of the low elastic underlayer bush is fixed. The outer surface of the plastic bush and the inner surface of the low elastic underlayer bush are bonded.
Figure 3b shows the water film model and its boundary condition. In this work the Reynolds number of the fluid is less than 1000, so the viscosity model was set to a laminar model. To consider cavitation, the multiphase flow model was used. The first phase was water liquid, and the second phase was water vapor. The cavitation model is the ZwartGerberBelamari model, and the cavitation pressure is 2340 Pa. The operating pressure is set to 101,325 Pa. Because the bearing is fully submerged, two sides of the water film model are set as “pressure inlet” and “pressure outlet” with gauge pressure as zero Pa. The outer surface of the water film is the FSI interface, and the inner surface is a rotating wall with a rotational speed equal to that of the journal. Figure 4 shows the mesh model of bush and water film. All grids are hexahedron grids. The sensitivity of mesh density to results was checked for both the fluid and solid domain. To obtain accurate results, the mesh distributions in the axial, circumferential and radial directions are 60 × 360 × 8 for the water film model, and 60 × 90 × 12 for the bush model.
In the FSI calculation, for models with small bush elastic modulus, negative grids often occur when updating the fluid model mesh, causing calculation error. Thus, in this work the elastic modulus of the low elastic underlayer bush was set as more than 200 MPa. To study the bearings' lubrication performance, the model parameters were set as: D = 60 mm, L/D = 1, c = 30 µm, E_{1 }= 1000 MPa, E_{2 }= 200–1000 MPa, ν_{1 }= 0.4, ν_{2 }= 0.47, t_{1} + t_{2} =10 mm. To compare the lubrication performance of doublelayer bearing and those of low elastic material bearing and plastic bearing, the following parameters were used: ε = 0.6, N = 1500 rpm, t_{1 }= t_{2 }= 5 mm, E_{1 }= 1000 MPa, E_{2 }= 400 MPa. To study the influence of eccentricity ratio on the lubrication performance of doublelayer bearing, ε ranged from 0.2 to 0.6, and N = 1500 rpm, t_{1 }= t_{2 }= 5 mm, E_{2 }= 400 MPa. To study the influence of the journal rotational speed on the lubrication performance of this kind of bearing, N ranged from 500 rpm to 2500 rpm, and ε = 0.6, t_{1 }= t_{2 }= 5 mm, E_{2 }= 400 MPa. To study the influence of elastic modulus of low elastic underlayer bush, E_{2} ranged from 200 to 1000 MPa, and ε = 0.6, N = 1500 rpm, t_{1 }= t_{2 }= 5 mm. To study the influence of the bearing bush thickness distribution, the low elastic underlayer bush thickness t_{2} ranged from 1 to 5 mm, and t_{1} ranged from 9 to 5 mm, and ε = 0.6, N = 1500 rpm, E_{2 }= 400 MPa.
Fig. 3 Model and boundary conditions of (a) Bush and (b) Water film. 
Fig. 4 Mesh model of (a) Bush and (b) Water film. 
4 Results and discussion
4.1 Lubrication performance
Figures 5 and 6 compare the water film pressure distribution and the bush radial deformation distribution of doublelayer bearing (t_{1 }= t_{2 }= 5 mm, E_{1 }= 1000 MPa, E_{2 }= 400 MPa, ν_{1 }= 0.4, ν_{2 }= 0.47), singlelayer plastic bearing (t_{1 }= 10 mm, t_{2 }= 0, E_{1 }= 1000 MPa, ν_{1 }= 0.4) and singlelayer low elastic material bearing (t_{1 }= 0, t_{2 }= 10 mm, E_{2 }= 400 MPa, ν_{2 }= 0.47) for ε = 0.6, N = 1500 rpm. Figure 7 shows the total bush deformation and the water film thickness in the midplane of the three kinds of bearings. Table 1 summarizes the lubrication performance parameters of the three kinds of bearings.
It can be seen that the water film pressure contours have positive pressure zone and cavitation zone, and the bush deformation is positive at the positive pressure zone, and negative at the cavitation zone. The maximum bush deformation occurs at the position with the maximum water film pressure. The maximum bush deformation of doublelayer bearing, plastic bearing and low elastic material bearing is 2.56 µm, 1.84 µm and 4.06 µm respectively, which means that the maximum bush deformation of doublelayer bearing is 0.72 µm more than that of plastic bearing, but 1.5 µm less than that of low elastic material bearing. The minimum water film thickness of low elastic material bearing is larger than those of doublelayer bearing and plastic bearing. The maximum water film pressure of doublelayer bearing, plastic bearing and low elastic material bearing is 0.39 MPa, 0.415 MPa and 0.34 MPa respectively, and the load carrying capacity is 835 N, 854 N, 788 N respectively. The load carrying capacity of doublelayer bearing decreases by 2.22% as compared with that of plastic bearing, but increases by 5.96% as compared with that of low elastic material bearing. The friction coefficient of doublelayer bearing, plastic bearing and low elastic material bearing is 0.00272, 0.00269 and 0.00280 respectively, which means that the friction coefficient of doublelayer bearing is lower than that of low elastic material bearing and close to plastic bearing. It can be found that the lubrication performance of doublelayer bearing is more like that of plastic bearing.
Fig. 5 The water film pressure distribution contour of (a) Doublelayer bearing, (b) Plastic bearing, and (c) Low elastic material bearing for ε = 0.6, N = 1500 rpm. 
Fig. 6 Bush radial deformation contours of (a) Doublelayer bearing, (b) Plastic bearing, and (c) Low elastic material bearing. 
Fig. 7 (a) Bush total deformation distribution, (b) Water film thickness distribution in the midplane of bearing. 
The lubrication performance of the bearings.
4.2 Effect of eccentricity ratio
Figures 8a–8c shows the radial deformation distribution curves of doublelayer bush, low elastic underlayer bush, and plastic layer bush in the midplane of bearing for different eccentricity ratio, and N = 1500 rpm, t_{1 }= t_{2 }= 5 mm, E_{1 }= 1000 MPa, E_{2 }= 400 MPa. The plastic layer deformation is equal to the total bush deformation minus the low elastic underlayer deformation. It can be seen that the trend of the deformation curves is similar, and the bush deformation increases significantly as the eccentricity ratio increases. The maximum deformation appears at around θ = 153°. Figure 8d shows the water film thickness distribution for different eccentricity ratio. It can be seen that the minimum water film thickness decreases as the eccentricity ratio rises.
Figures 9a–9e shows the variation of load carrying capacity, attitude angle, maximum water film pressure, friction coefficient and the maximum bearing bush deformation with eccentricity ratio for different rotational speed. It can be seen that the load carrying capacity, maximum water film pressure and maximum bush deformation increase as the eccentricity ratio rises. The load carrying capacity increases from 431 N to 1255 N as eccentricity ratio moves up from 0.2 to 0.6 for N = 2500 rpm. Here it increases about 2.9 times. While without considering the bush deformation, the load carrying capacity increases from 437 N to 1408 N, about 3.2 times. The deformation of bush has a great impact on the lubrication performance of the bearing especially for large eccentricity ratios. For N = 2500 rpm, the maximum water film pressure increases from 0.114 MPa to 0.588 MPa, and the maximum deformation of bearing bush increases from 0.75 µm to 3.89 µm. The attitude angle decreases as the eccentricity ratio rises, from 81.9° to 52.6°, and the friction coefficient of the bearing decreases from 0.0071 to 0.0029 as the eccentricity ratio rises.
Fig. 8 (a) Bush total deformation distribution, (b) Deformation distribution of low elastic bush, (c) Deformation distribution of plastic layer bush, (d) Water film thickness distribution for different eccentricity ratio (N = 1500 rpm, t_{1 }= t_{2 }= 5 mm, E_{1 }= 1000 MPa, E_{2 }= 400 MPa, ν_{1 }= 0.4, ν_{2 }= 0.47). 
Fig. 9 Effect of eccentricity ratio on the lubrication performance of doublelayer bearing (a) Load carrying capacity, (b) Attitude angle, (c) Maximum water film pressure, (d) Friction coefficient, (e) Maximum bearing bush deformation, for t1 = t2 = 5 mm, E1 = 1000 MPa, E2 = 400 MPa, n1= 0.4, n2 = 0.47. 
4.3 Effect of Rotational speed
Figures 10a–10c shows the radial deformation distribution curves of doublelayer bush, low elastic underlayer bush, and plastic layer bush in the midplane of bearing for different rotational speed. It can be seen that the maximum bush deformation increases significantly as the rotational speed rises. Figure 10d shows the water film thickness for different rotational speed. With the increase of the rotational speed, the water film thickness in the wedge convergence zone increases, and that in the divergent region changes little.
Figures 11a–11e shows the variation of load carrying capacity, attitude angle, maximum water film pressure, friction coefficient and the maximum bearing bush deformation with rotational speed for different eccentricity ratio. It can be seen that the load carrying capacity, maximum water film pressure and maximum bush deformation increase as the rotational speed rises. When without considering the bush deformation, the load carrying capacity and the maximum pressure increase almost linearly with the rotational speed, but after considering the bush deformation, the increasing rate slows down especially for larger rotational speeds. The load carrying capacity increases from 343 N to 1255 N, about 3.7 times, as rotational speed moves up from 500 rpm to 2500 rpm for ε = 0.6. The water film maximum pressure increases from 0.15 MPa to 0.59 MPa, about 3.9 times. The maximum bush deformation increases from 0.98 µm to 3.89 µm. The attitude angle decreases as the rotational speed rises, from 79.8° to 52.6°. The friction coefficient increases from 0.0024 to 0.0029 for ε = 0.6.
Fig. 10 (a) Bush total deformation distribution, (b) Deformation distribution of low elastic underlayer bush, (c) Deformation distribution of plastic layer bush, (d) Water film thickness distribution for different rotational speed (ε = 0.6, t_{1 }= t_{2 }= 5 mm, E_{1 }= 1000 MPa, E_{2 }= 400 MPa). 
Fig. 11 Effect of rotational speed on the lubrication performance of doublelayer bearing (a) Load carrying capacity, (b) Attitude angle, (c) Maximum water film pressure, (d) Friction coefficient, (e) Maximum bearing bush deformation, for t1 = t2= 5 mm, E1 = 1000 MPa, E2 = 400 MPa, n1 = 0.4, n2 = 0.47. 
4.4 Effect of bush elastic modulus
Figures 12a–12c shows the radial deformation distribution of doublelayer bush, low elastic underlayer bush, and plastic layer bush in the midplane of bearing for different low elastic underlayer bush elastic modulus. It can be seen that the bush deformation decreases as the low elastic underlayer bush elastic modulus increases, and the maximum deformation still appears at around θ = 153°. With the increase of the elastic modulus of low elastic underlayer, the deformation curve of low elastic underlayer changes significantly, and the deformation curve of plastic layer bush changes relatively small, which indicates that the change of the total deformation mostly comes from the change of the low elastic underlayer bush deformation. Figure 12d shows the water film thickness distribution for different low elastic underlayer bush elastic modulus. As the low elastic underlayer bush elastic modulus increases from 200 MPa to 1000MPa, the water film thickness in wedge convergence zone decreases, and the maximum reduction is 1.7 µm.
Figures 13a–13e shows the variation of load carrying capacity, attitude angle, maximum water film pressure, friction coefficient and the maximum bearing bush deformation with the low elastic underlayer bush elastic modulus for different eccentricity ratio. Low elastic underlayer bush elastic modulus has a significant impact on the maximum bush deformation especially for large eccentricity ratios, but has little impact on the other lubrication performance parameters. This also indicates that the lubrication performance of doublelayer bearing is more similar to that of plastic bearing. For ε = 0.6, as low elastic underlayer bush elastic modulus decreases from 1000 MPa to 200 MPa, the maximum bush deformation increases from 1.8 µm to 3.5 µm, almost doubles. But the load carrying capacity decreases from 854 N to 810 N, only dropping 5.2%. The water film maximum pressure decreases from 0.41 MPa to 0.36 MPa, the attitude angle decreases slowly from 61.5° to 58.2°, and the friction coefficient changes little.
Fig. 12 (a) Bush total deformation distribution, (b) Deformation distribution of low elastic underlayer bush, (c) Deformation distribution of plastic layer bush, (d) Water film thickness distribution for different low elastic underlayer bush elastic modulus (ε = 0.6, N = 1500 rpm, t_{1 }= t_{2 }= 5 mm, E_{1 }= 1000 MPa, ν_{1 }= 0.4, ν_{2 }= 0.47). 
Fig. 13 Effect of low elastic underlayer bush elastic modulus on the lubrication performance of doublelayer bearing (a) Load carrying capacity, (b) Attitude angle, (c) Maximum water film pressure, (d) Friction coefficient, (e) Maximum bearing bush deformation, for N = 1500 rpm, t1 = t2 = 5 mm, E1 = 1000 MPa, n1 = 0.4, n2 = 0.47. 
4.5 Effect of bush thickness distribution
Figures 14a–14c shows the radial deformation distribution curves of doublelayer bush, low elastic underlayer bush, and plastic layer bush in the midplane of bearing for different low elastic underlayer bush thickness. The total bush thickness is 10 mm, as the low elastic underlayer bush thickness increases from 1 mm to 5 mm, the plastic layer bush thickness decreases from 9 mm to 5 mm. With the increase of low elastic underlayer bush thickness, the deformation of low elastic underlayer bush increases and the deformation of plastic layer bearing bush decreases, but the total bush deformation changes little. Figure 14d shows the water film thickness distribution for different low elastic underlayer bush thickness. It can be seen that the low elastic underlayer bush thickness has little impact on the water film thickness.
Figure 15a–15e shows the variation of load carrying capacity, attitude angle, maximum water film pressure, friction coefficient and the maximum bush deformation with the low elastic underlayer bush thickness for different eccentricity ratio. It can be seen that when the low elastic underlayer bush thickness increases from 1 mm to 5 mm, the lubrication performance parameters change little. For ε = 0.6, the load carrying capacity decreases from 841 N to 836 N, only with a 0.6% drop. The attitude angle decreases from 60.3° to 59.9°, the maximum water film pressure decreases from 0.397 MPa to 0.390 MPa, the maximum bush deformation increases from 2.46 µm to 2.6 µm, and the friction coefficient almost remains the same.
Fig. 14 (a) Bush total deformation distribution, (b) Deformation distribution of low elastic underlayer bush, (c) Deformation distribution of plastic layer bush, (d) Water film thickness distribution for different low elastic underlayer bush thickness (ε = 0.6, N = 1500 rpm, E_{1 }= 1000 MPa, E_{2 }= 400 MPa, ν_{1 }= 0.4, ν_{2 }= 0.47). 
Fig. 15 Effect of low elastic underlayer bush thickness on the lubrication performance of doublelayer bearing (a) Load carrying capacity, (b) Attitude angle, (c) Maximum water film pressure, (d) Friction coefficient, (e) Maximum bearing bush deformation, for N = 1500 rpm, E1 = 1000 MPa, E2 = 400 MPa, n1 = 0.4, n2 = 0.47. 
5 Conclusions
This work studies the elastohydrodynamic lubrication performance of doublelayer composite waterlubricated bearings using FSI method, and analyzes the influence of rotational speed, eccentricity ratio, low elastic underlayer bush elastic modulus and bush thickness distribution on the lubrication performance of this kind of bearing. The following conclusions are drawn:
The load carrying capacity of doublelayer composite waterlubricated bearing is a little lower than that of plastic bearing and much higher than that of low elastic material bearing. The friction coefficient of doublelayer bearing is lower than that of low elastic material bearing and close to plastic bearing. The elastohydrodynamic lubrication performance of doublelayer bearing is more like that of plastic bearing.
Eccentricity ratio and journal rotational speed have significant impacts on the lubrication performance of the bearing. As the eccentricity ratio and journal rotational speed increases, the maximum pressure, the load carrying capacity and the maximum bush deformation increase significantly.
Low elastic underlayer bush elastic modulus has a significant impact on the maximum bush deformation especially for large eccentricity ratios, but has little impact on the other lubrication performance. For ε = 0.6, as low elastic underlayer bush elastic modulus decreases from 1000 MPa to 200 MPa, the maximum bush deformation almost doubles, but the load carrying capacity drops only 5.2%. The change of the total deformation mainly comes from the low elastic underlayer bush deformation, which indicates that the doublelayer bearing has better load carrying capacity than low elastic material bearing.
Bush thickness distribution influences the deformation of both low elastic underlayer bush and plastic bush, but have little impact on the total bush deformation. The load carrying capacity, the attitude angle, the maximum water film pressure and the friction coefficient change little as the low elastic underlayer bush thickness increases from 1 mm to 5 mm.
Nomenclature
p_{max}: Maximum water film pressure
t_{1}: Thickness of plastic bush
ν_{1}: Poisson's ratio of plastic bush
E_{1}: Elastic modulus of plastic bush
E_{2}: Elastic modulus of low elastic underlayer bush
h_{min}: Minimum film thickness
N: Rotational speed of journal
S: Maximum deformation of bearing bush
t_{2}: Thickness of low elastic underlayer bush
ν_{2}: Poisson's ratio of low elastic underlayer bush
θ: Bearing angle, or angular coordinate
Funding Information
The financial supports from National Natural Science Foundation of China (Grant No. 51905317) and Shandong Provincial Key Laboratory of Precision Manufacturing and Nontraditional Machining are gratefully acknowledged.
References
 W. Litwin, Waterlubricated bearings of ship propeller shafts − problems, experimental tests and theoretical investigations, Pol. Maritime Res. 16, 41–49 (2009) [Google Scholar]
 H. Wang, Z. Liu, L. Zou et al., Influence of both friction and wear on the vibration of marine water lubricated rubber bearing, Wear 376–377, 920–930 (2017) [Google Scholar]
 E. Piątkowska, Influence of solid particle contamination on the wear process in water lubricated marine strut bearings with NBR and PTFE bushes, Pol. Maritime Res. 28, 167–178 (2022) [Google Scholar]
 Z. Yan, X. Zhou, H. Qin et al., Study on tribological and vibration performance of a new UHMWPE/graphite/NBR water lubricated bearing material, Wear 332–333, 872–878 (2015) [Google Scholar]
 A.B. Ali, M.A. Samad, N. Merah, UHMWPE hybrid nanocomposites for improved tribological performance under dry and waterlubricated sliding conditions, Tribol. Lett. 65, 102 (2017) [Google Scholar]
 Z. Yang, Z. Guo, Z. Yang et al., Study on tribological properties of a novel composite by filling microcapsules into UHMWPE matrix for water lubrication, Tribol. Int. 153, 106629 (2021) [Google Scholar]
 Duramax Marine, https://www.duramaxmarine.com/advancedbearings/romoridovetailstavebearings.htm (accessed 17 February 2022) [Google Scholar]
 M. Wodtke, W. Litwin, Waterlubricated stern tube bearing − experimental and theoretical investigations of thermal effects, Tribol. Int. 153, 106608 (2021) [Google Scholar]
 W. Litwin, Properties comparison of rubber and three layer PTFENBRbronze water lubricated bearings with lubricating grooves along entire bush circumference based on experimental tests, Tribol. Int. 90, 404–411 (2015) [Google Scholar]
 W. Litwin, Experimental research on water lubricated three layer sliding bearing with lubrication grooves in the upper part of the bush and its comparison with a rubber bearing, Tribol. Int. 82, 153–161 (2015) [Google Scholar]
 D. Cabrera, N. Woolley, D. Allanson et al., Film pressure distribution in waterlubricated rubber journal bearings, Proc. Inst. Mech. Eng. J 219, 125–132 (2005) [Google Scholar]
 A. Kraker, R. Ostayen, D. Rixen, Calculation of Stribeck curves for (water) lubricated journal bearings, Tribol. Int. 40, 459–469 (2007) [Google Scholar]
 Y. Liao, Y. Zhou, Analysis on the lubrication characteristics of a few aqueous media that affect water lubricated rubber alloy bearing, J. Mech. Sci. Technol. 31, 4771–4779 (2017) [Google Scholar]
 Y. Wang, Z. Yin, D. Jiang, G. Gao, X. Zhang, Study of the lubrication performance of waterlubricated journal bearings with CFD and FSI method, Ind. Lubric. Tribol. 68, 341–348 (2016) [CrossRef] [Google Scholar]
 X.S. Lao, C.H. Dai, Y. Liu et al., Performance study on modified UHMWPE composites water lubricated bearing material, Mater. Sci. Forum 947, 167–172 (2019) [Google Scholar]
 G. Zhou, J. Wang, Y. Han et al., An experimental study on film pressure circumferential distribution of waterlubricated rubber bearings with multiple grooves, Tribol. Trans. 60, 385–391 (2017) [Google Scholar]
 G. Zhou, X. Mi, J. Wang et al., Experimental comparison between the Stribeck curves of water lubricated rubber bearing with straight and spiral grooves, Ind. Lubric. Tribol. 70, 1326–1330 (2018) [CrossRef] [Google Scholar]
 Z. Zhao, Z. Rui, Theoretical and experimental analysis of a waterlubricated rubber journal bearing with a large aspect ratio, Ind. Lubric. Tribol. 72, 797–803 (2020) [CrossRef] [MathSciNet] [Google Scholar]
 Y. Wang, X. Shi, L. Zhang, Experimental and numerical study on waterlubricated rubber bearings, Ind. Lubric. Tribol. 66, 282–288 (2019) [Google Scholar]
 L. Huang, Y. Jin, Z. Liu, L. Zou, The effect of damping layer on the dynamic performance of waterlubricated rubber bearings, Noise Vibr. Control 36, 32–37 (2016) [Google Scholar]
 J.L. Cai, G. Xiang, S. Li, Mathematical modeling for nonlinear dynamic mixed friction behaviors of novel coupled bearing lubricated with low viscosity fluid, Phys. Fluids 34 (2022) [Google Scholar]
 ANSYS, ANSYS F FLUENT, version 14.0: user manual (ANSYS Inc., Canonsburg, USA, 2011) [Google Scholar]
 P.J. Zwart, A.G. Gerber, T. Belamri, A twophase flow model for predicting cavitation dynamics, in Proceedings of the fifth international conference on multiphase flow, Yokohama, Japan (2004) [Google Scholar]
Cite this article as: Z. Yang, X. Zhang, Y. Wang, Y. Zhao, Study on elastohydrodynamic lubrication performance of doublelayer composite waterlubricated bearings, Mechanics & Industry 24, 3 (2023)
All Tables
All Figures
Fig. 1 Waterlubricated doublelayer bearing. 

In the text 
Fig. 2 Flowchart of calculation process. 

In the text 
Fig. 3 Model and boundary conditions of (a) Bush and (b) Water film. 

In the text 
Fig. 4 Mesh model of (a) Bush and (b) Water film. 

In the text 
Fig. 5 The water film pressure distribution contour of (a) Doublelayer bearing, (b) Plastic bearing, and (c) Low elastic material bearing for ε = 0.6, N = 1500 rpm. 

In the text 
Fig. 6 Bush radial deformation contours of (a) Doublelayer bearing, (b) Plastic bearing, and (c) Low elastic material bearing. 

In the text 
Fig. 7 (a) Bush total deformation distribution, (b) Water film thickness distribution in the midplane of bearing. 

In the text 
Fig. 8 (a) Bush total deformation distribution, (b) Deformation distribution of low elastic bush, (c) Deformation distribution of plastic layer bush, (d) Water film thickness distribution for different eccentricity ratio (N = 1500 rpm, t_{1 }= t_{2 }= 5 mm, E_{1 }= 1000 MPa, E_{2 }= 400 MPa, ν_{1 }= 0.4, ν_{2 }= 0.47). 

In the text 
Fig. 9 Effect of eccentricity ratio on the lubrication performance of doublelayer bearing (a) Load carrying capacity, (b) Attitude angle, (c) Maximum water film pressure, (d) Friction coefficient, (e) Maximum bearing bush deformation, for t1 = t2 = 5 mm, E1 = 1000 MPa, E2 = 400 MPa, n1= 0.4, n2 = 0.47. 

In the text 
Fig. 10 (a) Bush total deformation distribution, (b) Deformation distribution of low elastic underlayer bush, (c) Deformation distribution of plastic layer bush, (d) Water film thickness distribution for different rotational speed (ε = 0.6, t_{1 }= t_{2 }= 5 mm, E_{1 }= 1000 MPa, E_{2 }= 400 MPa). 

In the text 
Fig. 11 Effect of rotational speed on the lubrication performance of doublelayer bearing (a) Load carrying capacity, (b) Attitude angle, (c) Maximum water film pressure, (d) Friction coefficient, (e) Maximum bearing bush deformation, for t1 = t2= 5 mm, E1 = 1000 MPa, E2 = 400 MPa, n1 = 0.4, n2 = 0.47. 

In the text 
Fig. 12 (a) Bush total deformation distribution, (b) Deformation distribution of low elastic underlayer bush, (c) Deformation distribution of plastic layer bush, (d) Water film thickness distribution for different low elastic underlayer bush elastic modulus (ε = 0.6, N = 1500 rpm, t_{1 }= t_{2 }= 5 mm, E_{1 }= 1000 MPa, ν_{1 }= 0.4, ν_{2 }= 0.47). 

In the text 
Fig. 13 Effect of low elastic underlayer bush elastic modulus on the lubrication performance of doublelayer bearing (a) Load carrying capacity, (b) Attitude angle, (c) Maximum water film pressure, (d) Friction coefficient, (e) Maximum bearing bush deformation, for N = 1500 rpm, t1 = t2 = 5 mm, E1 = 1000 MPa, n1 = 0.4, n2 = 0.47. 

In the text 
Fig. 14 (a) Bush total deformation distribution, (b) Deformation distribution of low elastic underlayer bush, (c) Deformation distribution of plastic layer bush, (d) Water film thickness distribution for different low elastic underlayer bush thickness (ε = 0.6, N = 1500 rpm, E_{1 }= 1000 MPa, E_{2 }= 400 MPa, ν_{1 }= 0.4, ν_{2 }= 0.47). 

In the text 
Fig. 15 Effect of low elastic underlayer bush thickness on the lubrication performance of doublelayer bearing (a) Load carrying capacity, (b) Attitude angle, (c) Maximum water film pressure, (d) Friction coefficient, (e) Maximum bearing bush deformation, for N = 1500 rpm, E1 = 1000 MPa, E2 = 400 MPa, n1 = 0.4, n2 = 0.47. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.