Issue 
Mechanics & Industry
Volume 19, Number 2, 2018



Article Number  201  
Number of page(s)  13  
DOI  https://doi.org/10.1051/meca/2018005  
Published online  03 September 2018 
Regular Article
Numerical analysis and experimental research on load carrying capacity of waterlubricated tiltingpad thrust bearings
^{1}
School of Mechanical Engineering, Shandong University of Technology,
Zibo, PR China
^{2}
State Key Laboratory of Mechanical System and Vibration, Shanghai Jiao Tong University,
Shanghai, PR China
^{*} email: gaoqi_1118@163.com
Received:
10
August
2017
Accepted:
10
January
2018
Waterlubricated bearings are expected to be widely used because of convenience, green, safe and energy saving. The purpose of this study is to investigate the load carrying property of waterlubricated tiltingpad thrust bearings. A large amount of numerical analyses are undertaken based on computational fluid dynamics and the optimization method of pivot location and the calculation method of minimum film thickness are summarized. A thrust bearing is designed according to the numerical results and is tested by experiments. The experimental results validate the numerical method and the minimum film thickness to surface roughness ratio corresponding to the change of bearing lubrication regime from mixed lubrication to hydrodynamic lubrication is obtained.
Key words: Tiltingpad thrust bearing / water lubrication / CFD / experiment
© AFM, EDP Sciences 2018
1 Introduction
In turbomachinery, hydrodynamic thrust bearing is usually used to balance axial load. It plays an important role in movement security; its failing could involve serious damage following the possible contact between the rotor and the stator. Correct choice of the bearing type and prediction of its performance are therefore of great importance for trouble free machine operation [1]. Tiltingpad thrust bearing has the advantage of excellent stability, superior durability and high load carrying capacity. Thus this study focuses on the tiltingpad thrust bearing.
Presently the conventional lubricant is oil, but it has many disadvantages such as serious waste of resources, environmental pollution caused by oil leakage, explosive hazard, large bulk and complex structure, poor flexibility and maintainability [2]. In some conditions such as nuclear pumps, hydroturbine, marine pumps and integrated motor propulsors, waterlubrication replaces oillubrication. Heretofore, researches on waterlubricated bearings have focused primarily on the friction and wear characteristics of different bearing materials under water lubrication [2–6], and the lubrication performance of waterlubricated journal bearings [7–16], only a few of them are about waterlubricated thrust bearings. Zhang et al. [17] gave a brief summary of the influence factors of bearing capacity of waterlubricated thrust bearing and they pointed out that in order to guarantee the normal work of the bearing, the bearing capacity should be improved by changing factors including material of bearing and shaft, roughness of grinding area, bearing structure, pad number, pivot location, speed and the medium temperature, etc. Zhang et al. [18] studied the load carrying property of hydrodynamic waterlubricated step thrust bearing using CFD method and proposed the design method of this kind of bearing. Nakano et al. [19] developed a prototype of the nextgeneration, highperformance microturbine system using waterlubricated journal and thrust bearings and experimental results showed that the rotor system achieved stable rotating conditions at a rated rotational speed of 51 000 rpm. Kansai and Hitachi Mitsubishi Hydro developed the waterlubricated tiltingpad thrust bearing for a vertical type hydraulic turbine generator. They studied the effect of support structure on the bearing film thickness by experiments and found that the offset support structure is much more stable than the center support structure [20].
Waterlubricated bearing has advantages of convenience, green, safe and energy saving [21], but its load carrying capacity is much smaller than that of oillubricated bearings. Works on waterlubricated bearings [4,22–24] showed that the bearings usually operate in the regime of mixed lubrication and the lubrication mode might convert to hydrodynamic when load decreases. Wang et al. [22,23] studied the critical load, considering the sudden increase in the friction coefficient as the transition of lubrication mode. In order to ensure that waterlubricated tiltingpad bearing works under hydrodynamic lubrication as far as possible, the present work first studies the design method of pivot location and the calculation method of minimum film thickness by numerical analysis, and then studies the minimum film thickness to surface roughness ratio corresponding to the change of bearing lubrication regime from mixed lubrication to hydrodynamic lubrication by experiments.
2 Tiltingpad thrust bearing model
Tiltingpad thrust bearing can be studied using taperfaced thrust bearing model. The difference is that pressure distribution on tilting pad should meet moment balance on the pivot. The supporting position of tiltingpad thrust bearing is the pressure center of taperfaced thrust bearing model. The previous researches on thrust bearing usually simplified the sector pad to rectangular pad and ignored the groove. In this work, the flow model of one sector pad and one groove of tiltingpad thrust bearing as shown in Figure 1 is studied. The thrust bearing is submersed in water. The hydrodynamic action generates dynamic pressure in water, primarily in the convergent part of the thrust pad, to counteract the load thereby separating the ring surface from the bearing surface with a thin lubricant film. Geometry parameters − inner radius R_{1}, outer radius R_{2}, number of pads n, pad ratio k, and pad supporting position all influence the load carrying capacity of tiltingpad thrust bearings. In this work, the taperfaced thrust bearing models with different radii and pad inclined angles are calculated to study the effect of geometry parameters on load carrying capacity and supporting position of tiltingpad thrust bearing. According to the design of oillubricated thrust bearings [25], R_{2} is usually 1.5–3 times of R_{1}. The number of pads n is generally 6–12. The pad ratio k, the percentage of pad area in the whole thrust surface, is typically 0.7–0.85. Pad widthtolength ratio B/L is a pad parameter determined by R_{2}/R_{1}, n and k, and the relation is (1)
It has an effect on the optimum value of pad inclined angle β and pad supporting position R_{s} and θ_{s}.
Rotational speed N, fluid viscosity µ and the minimum film thickness h_{2} determine the bearing carrying performance as well. Rotational speed and fluid viscosity are determined by work conditions. The minimum film thickness should not be less than a safety value depending on surface roughness and system vibration, or the lubrication model might change from full dynamic lubrication to boundary or mixed lubrication.
Fig. 1 Flow model of one sector pad and one groove of tiltingpad thrust bearing (not to scale). 
3 Numerical analysis
3.1 Governing equations
The flow is considered laminar, steady and incompressible, with zero gravitational and other external body forces. The bearing is fully submerged into water. As water has a low viscosity (less than 0.001 Pa.s) and a large specific heat (4200 J/(kg.K)), the temperature rise of waterlubricated bearing is small. Thus the flow is considered isothermal. To solve such flow, the following governing equations must be solved:
Mass conservation equation (2)
Momentum conservation equations (3)
When flow enters the groove, pressure might fall below the saturation water vapor pressure, and the liquid would rupture and cavitation occurs. Thus cavitation is taken into account. In the present work the CFD code ANSYS FLUENT is used. There are three available cavitation models in ANSYS FLUENT: Singhal et al. model, ZwartGerberBelamri model and Schnerr and Sauer model. The Singhal et al. model is numerically less stable and more difficult to use. The ZwartGerberBelamri and the Schnerr and Sauer models are robust and converge quickly [26]. In this case the ZwartGerberBelamri model is employed.
In cavitation, the liquidvapor mass transfer (evaporation and condensation) is governed by the vapor transport equation [26]: (4) where R_{g} and R_{c} account for the mass transfer source terms connected to the growth and collapse of the vapor bubbles, α_{v} is vapor volume fraction and ρ_{v} is vapor density. In ZwartGerberBelamri model, R_{g} and R_{c} are defined as follows[27]: (5)where F_{evap} = evaporation coefficient = 50, F_{cond} = condensation coefficient = 0.01, R_{b} = bubble radius = 10^{−6} m, α_{nuc} = nucleation site volume fraction = 5 × 10^{−4}, ρ_{l} = liquid density, p_{v} = pressure of vapor.
Supposed that there's no installation error and the load of a thrust bearing is distributed uniformly over all pads, the load and friction torque of a thrust bearing can be calculated by integrating the pressure and shear stress over the rotating wall as follows: (6)
The friction coefficient is (7)
The supporting position (R_{s}, θ_{s}) of tiltingpad thrust bearing meets the following formula (8)
To compare with results in the literature, nondimensional pressure, nondimensional load carrying capacity and nondimensional friction torque are defined as (9)
The nondimensional supporting position is defined as (10)
3.2 CFD models
The flow model is built and meshed in Gambit 2.3. The boundary condition is set as shown in Figure 2. Rotational periodic boundary condition is used to simplify the flow model and to reduce the computational cost. The operating pressure is set to 1 01 325 Pa. Since the bearing is fully submerged, the pressure at the inlet and outlet boundaries is taken as zero (gauge pressure). A noslip condition is imposed on the solid walls. The bottom wall is stationary and the upper one is assumed to be rotating at a constant rotational speed N.
In the present study, water properties at 20 °C listed in Table 1 are employed. In all cases, the groove depth h_{g} = 2 mm and pad ratio k = 0.75. For the purpose of finding relations between load carrying properties and pad parameters, the inner radius R_{1} is set to 30 mm, number of pad n is 12, rotational speed Nis 3000 rpm, R_{2} varies from 35 mm to 70;mm and δ^{*} varies from 0.4 to 3. The nondimensional inclined plane height δ^{*} varies to seek out the optimum ones. The outer radius, R_{2}, changes to find the effect of B/L on the nondimensional load F_{z}^{*}, the optimum δ^{*} and the optimum supporting position. Pad widthtolength ratio B/L can be calculated from equation (1).
Due to the existence of the thin film in hydrodynamic bearings, the dominating feature for bearing CFD models is the large aspect ratio of the grid which is 330∼4500 in this study. However, the suggested value of aspect ratio is less than 100 or 200 for normal analysis. Thus double precision calculations are employed to avoid the negative influence of the large aspect ratio according to ANSYS FLUENT user manual [26].The CFD models are meshed using hexahedron grids in Gambit 2.3. In order to obtain accurate solutions, a mesh refinement study is carried out. Since the groove area contributes little to the load carrying capacity, the groove length L_{g} is divided into 15 cells and the groove height h_{g} is divided into 20 cells for all cases. An expansion ratio of 1.4 is employed in h_{g} to avoid sudden change of mesh density near the thin film thickness. A uniform mesh is employed in other edges of the model. The influence of the mesh density of h_{1}, L, Bon load for the bearing model with R_{1} = 30 mm, R_{2} = 45 mm, n = 12, h_{2} = 10 µm, β = 0.056°, N = 3000 rpm is studied and the results are summarized in Table 2. As cell number of h_{1} increases, computed result increases. As cell number of L(B) increases, computed result decreases. Based on the result, case 8 is employed in this study. Appropriate mesh sizes are also found for other sizes of bearing models.
Fig. 2 Boundary conditions. 
Water properties at 20 °C.
Mesh refinement study.
3.3 Results and discussion
3.3.1 Pressure distribution
Figures 3a–c present the gauge pressure distribution of tiltingpad thrust bearing for different inclined angles, and D_{i} = 60 mm, D_{o} = 90 mm, h_{2} = 10 µm, N = 3000 rpm. The pressure distribution shows that the maximum pressure of the whole film is located near the outlet. Figure 4 shows the maximum pressure versus nondimensional height of inclined plane, δ^{*}. It shows that as δ increases, the pad inclined angle increases and the maximum pressure first increases and then decreases. The maximum pressure reached the maximum value at δ* = 2, but the effective load region size decreases. So there is an optimum δ* to get the maximum load carrying capacity. When the fluid flows into the groove, pressure on the stationary wall drops rapidly, which might cause cavitation if it is less than 2340 Pa. As the fluid velocity near the outer radius is larger than that near the inner radius, the pressure drop is more significant and cavitation occurs more easily. From the CFD analyses, it is found that the cavitation zone will get larger when the minimum film thickness decreases or the fluid velocity increases.
Figure 5 compares the pressure distribution of tiltingpad thrust bearing for different load and rotational speed, and D_{i} = 60 mm, D_{o} = 90 mm, R_{s}^{*} = 0.51, θ_{s}^{*} = 0.575. From Figure 5a, b, it can be seen that when the load is constant, the pressure distribution changes little with the rotational speed, but the minimum film thickness increases as the rotational speed increases. From Figure 5b, c, it can be seen that when the rotational speed is constant, the pressure increases and the minimum film thickness decreases as load increases.
Figure 6a, b show the nondimensional pressure distribution at the radius of maximum pressure and at the medium radius of tiltingpad thrust bearing for different B/L respectively, and δ^{*} = 1, h_{2} = 10 µm, N = 3000 rpm. The nondimensional pressure distribution of infinite width slider bearing in [28] is shown for comparison. Figure 6a indicates that the maximum p* at the radius of maximum pressure increases as B/L increases, and it increase sharply when B/L is less than 1.5, but increase slowly when B/L is more than 1.5. Figure 6b shows that when B/L is more than 1.5, nondimensional pressure at the medium radius is larger than that at the radius of maximum pressure, and the nondimensional pressure distribution at the medium radius for B/L = 2.04 is close to that of infinite width slider bearing. This can be explained by the inertial effect of sector pad thrust bearing. The maximum pressure moves outwards due to the inertial effect, and it reacts on the pressure buildup at the medium radius and increases p_{max}* at the medium radius.
Fig. 3 Pressure distribution for pads with different δ*, and B/L = 1.02, h_{2} = 10 µm, N = 3000 rpm. 
Fig. 4 Maximum pressure, pmax, versus δ*, and B/L = 1.02, h_{2} = 10 µm, N = 3000 rpm. 
Fig. 5 Pressure distribution for different load and rotational speed. 
Fig. 6 Nondimensional pressure distribution at (a) the radius of maximum pressure (b) the medium radius of tiltingpad thrust bearing for different B/L, and δ^{*} = 1, h_{2 }= 10 µm, N = 3000 rpm. 
3.3.2 Optimization of load carrying capacity
Figure 7 presents the nondimensional load carrying capacity F_{z}^{*}versus height of inclined plane δ and the nondimensional height of inclined plane δ^{*} for different h_{2} with B/L = 1.02. Results illustrate that although the optimum δ is different for different minimum film thicknesses, the optimum ratio of them, the nondimensional height of inclined plane δ^{*}, is always around 1.33. Figure 8 compares the optimum δ^{*}for pads with different B/L values. It can be seen that the optimum δ^{*} is 1.4 for B/L = 0.73 and the range of δ^{*} corresponding to more than 0.9 times the maximum F_{z}^{*} is [0.65, 2.6]; the optimum δ^{*} is 1.26 for B/L = 1.27 and the range of δ^{*} corresponding to more than 0.9 times the maximum F_{z}^{*} is [0.62, 2.42];the optimum δ^{*} is 1.2 for B/L = 2.04 and the range of δ^{*} corresponding to more than 0.9 times the maximum F_{z}^{*} is [0.6, 2.3].The optimum δ^{*} is summarized in Figure 9. The wider the pad, the smaller the optimum δ^{*}.The pad inclined angle can be calculated using the following equation: (11)
Fig. 7 Nondimensional load carrying capacity, F_{z}*, versus(a) height of inclined plane, δ, and (b) height of inclined plane to minimum film thickness ratio, δ*, for different minimum film thicknesses, h_{2}, and B/L = 1.02. 
Fig. 8 Nondimensional load carrying capacity, F_{z}*, versus height of inclined plane to minimum film thickness ratio, δ*, for pads with different B/L values, and h_{2} = 10 µm. 
Fig. 9 The optimum δ^{*}versus pad widthtolength ratio B/L. 
3.3.3 Optimization of pivot position
Pivot position need to be chosen appropriately when designing tiltingpad thrust bearing. For tiltingpad thrust bearing supported by ball, the pressure distribution on pad surface should meet moment balance equations on the pivot, so the pressure center is the supporting position. Figures 10 and 11 present the nondimensional radius, R_{s}^{*}, and nondimensional angle, θ_{s}^{*}, of pivot versus δ^{*} for pads with different B/L values. Figure 10 shows that R_{s}^{*} is greatly influenced by B/L and R_{s}^{*} increases as B/L increases; R_{s}^{*} is less affected by δ^{*} and R_{s}^{*}decreases as δ^{*} increases. Figure 11 shows that θ_{s}^{*} is greatly influenced by δ^{*} and θ_{s}^{*} increases as δ^{*} increases; θ_{s}^{*} is less affected by B/L and the variation of θ_{s}^{*} decreases as B/L increases. Figure 12 shows the change of pad pivot position with B/L and δ^{*} in a twodimensional coordinates.
The optimum pivot position of tiltingpad thrust bearing can be determined according to the optimum δ^{*}. The result is summarized in Figure 13. It shows that the optimum R_{s}^{*} ranges from 0.47 to 0.57 when B/L increases from 0.39 to 2.04; the optimum θ_{s}^{*} ranges from 0.567 to 0.577, and it reduces slightly with B/L when B/L is less than 1.5, while remains unchanged when B/L is more than 1.5.
Fig. 10 Nondimensional radius of pivot, R_{s}^{*}, versus δ* for pads with different B/L values, and h_{2 }= 10 µm. 
Fig. 11 Nondimensional angle of pivot, θ_{s}^{*}, versus δ*, for pads with different B/L values, and h_{2 }= 10 µm. 
Fig. 12 Change of pivot position with δ*, for pads with different B/L values, and h_{2 }= 10 µm. 
Fig. 13 Optimum position of pivot in tilting pad thrust bearing for different B/L. 
3.3.4 Check computation of the minimum film thickness
Figure 14 presents the obtained nondimensional load carrying capacity F_{z}^{*}versus pad widthtolength ratio B/L. It shows that F_{z}^{*} increases as B/L increases. The result is compared with that of rectangular pad obtained by Jakobson [29] using finite difference method. F_{z}^{*} is larger than that obtained by Jakobson. It indicates that sector pads provide greater load carrying capacity than rectangular ones and the difference increases as B/L increases. This is because as B/L increases, pressure center in a sector pad moves towards the outside radius, where the fluid velocity is much larger than that at the middle radius. Thus larger pressure and greater load are gained in a sector pad than those in a rectangular pad. It can also be seen that the slope of the load curve for small values of B/L is larger than that for large values of B/L. This is because the pressure buildup is more affected by side leakage for small B/L values.
The minimum film thickness h_{2} can be calculated from F_{z}^{*}.As illustrated in Figure 14, B/L is the main influence factor on the nondimensional load carrying capacity F_{z}^{*} for a constant δ^{*}. This means for a certain B/L value, F_{z}^{*} is definite. Then the minimum film thickness h_{2}, an important criterion of lubrication mode, can be calculated from equation (12): (12)
In order to ensure that the bearings operate under hydrodynamic lubrication, the surface roughness values of the thrust pad and the thrust ring should satisfy the following formula: (13) where R_{q,p} refers to rms (root mean square) surface roughness of the pad surface, R_{q,r} refers to rms surface roughness of the ring surface, S refers to the minimum film thickness to surface roughness ratio corresponding to the change of bearing lubrication regime from mixed lubrication to hydrodynamic lubrication, which can be obtained from experiments.
If the surface roughness values could not meet the design requirements, the bearing parameters should be corrected according to equations (12) and (13).
Fig. 14 Nondimensional load carrying capacity for different pad widthtolength ratio. 
4 Experimental test
In order to test the lubrication properties of waterlubricated tiltingpad thrust bearing and verify the theoretical calculation method, a waterlubricated bearing test rig is built, a thrust bearing is designed according to the numerical results and the friction coefficient and water film thickness of the tiltingpad thrust bearing under different load and speed are tested.
4.1 Waterlubricated bearing test rig
Figure 15 shows the waterlubricated bearing test rig. The test rig includes mechanical part and test system. Figure 16 shows schematic of the test rig. Shaft 4 is driven by servo motor and the shaft speed can be adjusted by motor driver. Shaft 4 is supported by rolling bearing 3 with high precision and its rotation accuracy reaches 2 µm. Test shaft 5 is fitted with shaft 4 through a conical surface and fixed by bolts. In order to test the friction coefficient of waterlubricated bearing directly, the test thrust bearing 8 is mounted on the ball seat 11 through a thrust ball bearing and a journal bearing component 10. The test bearing is connected to a force sensor 7 via a friction torque measuring rod 6, which can be seen more clearly in Figure 17. The error of friction torque of test bearing is mainly caused by the static friction force of component 10. When the friction torque of test bearing is less than the friction torque of component 10, the measurement is incorrect. In this case, an initial value of force sensor 7 is set artificially by tightening the line. So the measured friction torque may be a little larger. But the variation of the friction torque is accurate. The thrust bearing is loaded by pneumatic loading method. The thrust load is measured by force sensor 13.
To determine whether the waterlubricated bearing operates under hydrodynamic lubrication regime, the direct criterion is the minimum film thickness. When the minimum film thickness exceeds the minimum film thickness under mixed lubrication regime and reaches a certain value, the bearing will operated under hydrodynamic lubrication regime. The experiment aims to find the minimum film thickness to surface roughness ratio corresponding to the change of bearing lubrication regime from mixed lubrication to hydrodynamic lubrication. Four displacement sensors are used to measure the film thicknesses.
The test shaft is made of 42CrMo and it is electroplated with a layer of nickel for rust protection. The outer diameter of the test thrust plate is 150 mm and the surface roughness is Ra 0.3 µm.
The test system includes eddy current displacement sensors, force sensors, a data acquisition and analysis equipment and a computer. The performance parameters of the eddy current displacement sensor are: probe diameter 5 mm, starting position 0.5 ± 0.1 mm, range 1 mm, resolution 0.1 µm, temperature drift ≤0.05%/° C FS. Before the test, the eddy current sensor and the force sensor were calibrated by the manufacturer to ensure the measurement accuracy. As there is a problem of magnetic run out of the thrust collar when using the eddy current displacement sensor to test the change of film thickness, chebyshev lowpass filter is used to handle the signal.
Fig. 15 Waterlubricated bearing test rig. 
Fig. 16 Schematic of waterlubricated bearing test rig: 1servo motor; 2coupling; 3rolling bearings; 4shaft; 5test shaft; 6 measuring rod for friction torque; 7friction torque force sensor; 8thrust bearing; 9thrust bearing pedestal; 10thrust ball bearing and journal bearing; 11ball seat; 12guide rail; 13force sensor; 14air cylinder. 
Fig. 17 Measurement part of waterlubricated thrust bearing test rig. 
4.2 Test bearing and test method
Figure 18 shows structure of the test thrust bearing. It is composed of bushing and backing plate. The backing plate is made of stainless steel. The bearing bushing uses Polytetrafluoroethylene (PTFE) as base material and is reinforced with Polyetheretherketone (PEEK) and carbon fiber. The elastic modulus of the composite bush is about 1300 MPa. The pads of the test bearing are cut from a circular ring with a flat thrust surface, and the supporting ball surface and groove surface are processed and formed in one step. The inner diameter of the bearing pad is 80 mm, the outer diameter is 120 mm, and the surface roughness is 1.7 ∼ 2.2 µm. There are 8 pads and the water groove is designed as straight groove for easy machining. The width of water groove is 6 mm, so the pad widthtolength ratio is 0.6. The optimum nondimensional pivot position is R_{s}^{*} = 0.483, θ_{s}^{*} = 0.58 according to Figure 13, so the supporting position is designed as R_{s} = 49.65 mm, θ_{s} = 26.1°. The thrust pad is supported by a ball and limited by two pins as shown in Figure 19. Figure 20 shows the test waterlubricated thrust bearing.
As the thrust bearing is small, the displacement sensors can not be installed at the bearing bush to measure the water film thickness directly. So the sensors are installed outside the bearing bush. The mounting position of displacement sensors is shown in Figures 18 and 20. The water film thicknesses on two pads are measured and the average value of film thicknesses at the same location is used to eliminate the effect of shaft inclination. Figure 21 shows test method of the water film thickness. It can be seen that the measured data of sensor is: (14) where d_{1} is the distance from the sensor probe to the surface of bearing bush. d_{1} is constant as the sensor is fixed on the bearing backing plate. Therefore, before the test the position of sensors need to be calibrated and d_{1} is measured.
The measurement principle of bearing friction coefficient is introduced in Ch. 4.1. The friction coefficient can be calculated using the following equation: (14) where F_{f} is the measured value of force sensor 7, l is the distance from the shaft axis to the sensor measuring line, F_{z} is the load (see Fig. 17).
Fig. 18 Structure of test bearing. 
Fig. 19 Installation diagram of thrust pad. 
Fig. 20 Waterlubricated thrust bearing. 
Fig. 21 Measurement of thrust bearing film thickness. 
4.3 Results and discussion
4.3.1 Effect of load on bearing friction coefficient and film thickness
Figure 22 shows the change of friction coefficient of the waterlubricated thrust bearing with load for N = 600 rpm and 900 rpm. When the load is small, the bearing operates under hydrodynamic lubrication regime and the bearing friction is mainly the liquid shear force, which change little with load. So the friction coefficient decreases when load increases. When load continues to increases, asperity contact appears and the bearing friction includes the liquid shear force and the asperity contact friction. The larger the load, the smaller the film thickness and the greater the asperity contact frictional force. When the increase of friction is larger than the increase of load, the friction coefficient increases.
Figures 23 and 24 present the variation of water film thickness with load for N = 600 rpm and 900 rpm, respectively. The lines with marks show the mean film thicknesses at measuring point 1 and 3 or 2 and 4. The error bars present the differences between film thicknesses measured at measuring point 1 and 3 or 2 and 4. The solid line refers to the theoretical minimum film thickness. It can be seen that when the load increases, the minimum water film thickness decreases and the pad inclined angle increases. As the displacement sensor is not installed at the bearing bush, so there is deviation between measured value and theoretical value of film thickness, but the trend is the same. According to the film thickness of measuring point 1 and measuring point 2, the calculated inclined angle is about 0.01° ∼ 0.02°, which is in agreement with the theoretical value calculated by equation (11) in the order of magnitude. From the change of error bars, it can be seen that as load increases the difference between film thicknesses at the two pads decreases. This might be due to the bearing inclination. The inclination of bearing axis decreases as load increases.
From Figures 22 and 23, it can be seen that at 600 rpm the friction coefficient begins to increase when the specific pressure is 0.3 MPa, and the minimum water film thickness is 5.85 µm. From Figures 22 and 24, it can be seen that at 900 rpm the friction coefficient begins to increase when the specific pressure is 0.37 MPa, and the minimum water film thickness is 5.8 µm. The increasing friction coefficient indicates that the bearing has been in mixed lubrication regime, so the minimum film thickness corresponding to the change of bearing lubrication regime from mixed lubrication to hydrodynamic lubrication is larger than 5.8 µm.
Fig. 22 Friction coefficient versus specific pressure for N = 600 rpm and 900 rpm. 
Fig. 23 Film thickness versus specific pressure for N = 600 rpm. 
Fig. 24 Film thickness versus specific pressure for N = 900 rpm. 
4.4 Effect of rotational speed on bearing friction coefficient and film thickness
Figure 25 shows the change of friction coefficient with rotational speed for p_{s} = 0.286 MPa and 0.516 MPa. For small rotational speed, the water film thickness is small and there is asperity contact, so the friction coefficient is large. As rotational speed increases, water film thickness increases and the bearing friction coefficient decreases. When rotational speed increases to a certain value, hydrodynamic lubrication forms and the bearing friction is mainly liquid shear force, which increases with rotational speed, so the bearing friction coefficient increases when rotational speed continues to increase.
Figures 26 and 27 present the change of water film thickness with rotational speed for p_{s} = 0.286 MPa and 0.516 MPa. The lines with marks show the mean film thicknesses at measuring point 1 and 3 or 2 and 4. The error bars present the differences between film thicknesses measured at measuring point 1 and 3 or 2 and 4. The solid line shows the change of theoretical minimum film thickness. When the load is constant, the minimum water film thickness increases and the pad inclined angle decreases as rotational speed increases. The measured minimum film thicknesses are slightly larger than the theoretical values, but the trend is the same. The error bar is short for p_{s} = 0.286 MPa and 0.516 MPa and its change rule with rotational speed is not clear.
From Figures 25 and 26, it can be seen that for p_{s} = 0.286 MPa the bearing friction coefficient begins to increase when the rotational speed is about 900 rpm, and the minimum water film thickness is 7.2 µm. When p_{s} is 0.516 MPa, the friction coefficient does not show the rise when rotational speed is less than 1500 rpm, and the minimum water film thickness is 6.5 µm at N = 1500 rpm as shown in Figure 27. It indicates that the minimum water film thickness is about 7.2 µm when the bearing lubrication regime changes from mixed lubrication to hydrodynamic lubrication. Since the surface roughness of the test bearing is 1.7 ∼ 2.2 µm, the minimum film thickness to surface roughness ratio is about 3.6.
Fig. 25 Friction coefficient versus rotational speed for p_{s }= 0.286 MPa and 0.516 MPa. 
Fig. 26 Film thickness versus rotational speed for p_{s }= 0.286 MPa. 
Fig. 27 Film thickness versus rotational speed for p_{s }= 0.516 MPa. 
5 Conclusions
In this research, the hydrodynamic lubrication model of the waterlubricated tiltingpad thrust bearing with sector pads and straight radial grooves has been fully studied using CFD method. The effects of pad inclined angle, pivot position and pad widthtolength ratio on lubrication properties are discussed. Experiments are carried out to study the lubrication performance of tiltingpad thrust bearing and the correctness of theoretical method. The main conclusions are drawn as follows:

the optimum pivot position of tiltingpad thrust bearing can be chosen from Figure 14 according to B/L. The optimum R_{s}^{*} ranges from 0.47 to 0.57 when B/L increases from 0.39 to 2.04; the optimum θ_{s}^{*} ranges from 0.567 to 0.577, and it reduces slightly with B/L when B/L is less than 1.5, while remains unchanged when B/L is more than 1.5;

the nondimensional load carrying capacity of tiltingpad thrust bearing can be found in Figure 14. It varies from 0.022 to 0.125 as B/L increases from 0.39 to 2.04. Compared with rectangular pad, sector pad provides greater load carrying capacity and the difference increases as B/L increases;

the minimum film thickness should be checked using equations (12) and (13) in order to ensure that the bearings operate under hydrodynamic lubrication. Experiments showed that the minimum film thickness to surface roughness ratio corresponding to the change of bearing lubrication regime from mixed lubrication to hydrodynamic lubrication in equation (13) is about 3.6;

the measured film thicknesses are in accordance with the theoretical values, which validates the theoretical method.
The bush deformation and installation error of pads also influence the lubrication and load carrying properties. They have not been considered in this work. We'll study their effects in the future work.
Nomenclature
B/L: pad widthtolength ratio
F_{z}: load carrying capacity of a thrust bearing
F_{z}^{*}: nondimensional load carrying capacity
f: friction coefficient= M_{fr}/(F_{z}R_{m})
h_{1}, h_{2}: film thicknesses at inlet and outlet of clearance in middle radius
k: pad ratio = θ_{p}/(θ_{p}+θ_{g})
L: pad length (see Fig. 1), L = 2πR_{m}k/n
M_{fr}: friction torque on the thrust ring
M_{fr}^{*}: nondimensional friction torque= M_{fr}h_{2}(nμU_{m}LBR_{m})
p_{s}: specific pressure= F_{z}/(nBL)
p^{*}: nondimensional pressure
R_{1}, R_{m}, R_{2}: bearing inner, middle and outer radius
R_{s}^{*}: nondimensional radius of pivot= (R_{s} − R_{1})/B
U: linear velocity of the thrust collar = 2πNR, U_{m} = 2πNR_{m}
θ_{p}, θ_{g}: pad angle, groove angle
θ_{s}^{*}: nondimensional angle of pivot = θ_{s}/θ_{p}
δ: height of inclined plane = h_{1}−h_{2}
δ^{*}: nondimensional height of inclined plane = δ/h_{2}
Acknowledgements
The authors thank the editor and reviewers for their helpful comments and suggestions. The financial support from the National Natural Science Foundation of China (Grant No.11672178, Grant No.51705310 and Grant No.2017M611555) is also gratefully acknowledged.
References
 S.B. Glavatskih, Steady state performance characteristics of a tilting pad thrust bearing, J Tribol 123 (2001) 608 [CrossRef] [Google Scholar]
 J. Wang, F. Yan, Q. Xue, Tribological behavior of PTFE sliding against steel in sea water, Wear 267 (2009) 1634–1641 [CrossRef] [Google Scholar]
 J.H. Jia, J.J. Lu, H.D. Zhou, J.M. Chen, Tribological behavior of Nibased composite under distilled water lubrication, Mater Sci Eng A Struct Mater Prop Microstruct Process 381 (2004) 80–85 [CrossRef] [Google Scholar]
 W. Huang, Y. Xu, Y. Zheng, X. Wang, The tribological performance of Ti(C,N)based cermet sliding against Si3N4 in water, Wear 270 (2011) 682–687 [CrossRef] [Google Scholar]
 H. Unal, A. Mimaroglu, Friction and wear characteristics of PEEK and its composite under water lubrication, J Reinf Plast Compos 25 (2006) 1659–1667 [CrossRef] [Google Scholar]
 J. Jia, J. Chen, H. Zhou, J. Wang, H. Zhou, Friction and wear properties of bronzegraphite composite under water lubrication, Tribol Int 37 (2004) 423–429 [CrossRef] [Google Scholar]
 Y.D. Tridimas, D.R. Allanson, N.H. Woolley, D.L. Cabrera, Film pressure distribution in waterlubricated rubber journal bearings, Proc Inst Mech Eng Part J: J Eng Tribol 219 (2005) 125–132 [CrossRef] [Google Scholar]
 N. Wang, Q. Meng, P. Wang, T. Geng, X. Yuan, Experimental research on film pressure distribution of waterlubricated rubber bearing with multiaxial grooves, J Fluids Eng 135 (2013) 084501 [CrossRef] [Google Scholar]
 R. Pai, R.S. Pai, Stability of fouraxial and sixaxial grooved waterlubricated journal bearings under dynamic load, Proc Inst Mech Eng Part J: J Eng Tribol 222 (2008) 683–691 [CrossRef] [Google Scholar]
 B.C. Majumdar, R. Pai, D.J. Hargreaves, Analysis of waterlubricated journal bearings with multiple axial grooves, Proc Inst Mech Eng Part J: J Eng Tribol 218 (2004) 135–146 [CrossRef] [Google Scholar]
 B.D. Heberley, Advances in hybrid waterlubricated journal bearings for use in ocean vessels, Massachusetts Institute of Technology, Ann Arbor, USA, 2013. [Google Scholar]
 G.Y. Gao, Z.W. Yin, D. Jiang, X.L. Zhang, Numerical analysis of plain journal bearing under hydrodynamic lubrication by water, Tribol Int 75 (2014) 31–38 [CrossRef] [Google Scholar]
 W. Litwin, Marine water lubricated stern tube bearings calculations and measurements of heavy loaded bearings, in: Proceedings of the STLE/ASME International Joint Tribology Conference IJTC2008, Miami, Florida, USA, 2008. [Google Scholar]
 W. Litwin, Experimental investigation on marine main shaft bearings with reduced length to diameter ratio, in: Proceedings of the STLE/ASME 2010 International Joint Tribology Conference IJTC2010, San Francisco, California, USA 2010. [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 (2015) 153–161 [CrossRef] [Google Scholar]
 W. Litwin, C. Dymarski, Experimental research on waterlubricated marine stern tube bearings in conditions of improper lubrication and cooling causing rapid bush wear, Tribol Int 95 (2016) 449–455 [CrossRef] [Google Scholar]
 X. Zhang, X. Wang, X. Li, X. Chu, Study on influence factors of bearing capacity of waterlubricated thrust bearing and its improving methods, in Advanced Mechanical Design, 2012 [Google Scholar]
 X. Zhang, Z. Yin, D. Jiang, G. Gao, The design of hydrodynamic waterlubricated step thrust bearings using CFD method, Mech Ind 15 (2014) 197–206 [CrossRef] [EDP Sciences] [Google Scholar]
 S. Nakano, T. Kishibe, T. Inoue, H. Shiraiwa, An advanced microturbine system with waterlubricated bearings, Int J Rotating Mach (2009) 1–12 [CrossRef] [Google Scholar]
 K. Inoue, K. Deguchi, K. Okude, R. Fujimoto, Development of the waterlubricated thrust bearing of the hydraulic turbine generator, in: IOP Conference Series: Earth and Environmental Science, Vol. 15, 2012 072022 p. [CrossRef] [Google Scholar]
 R. Pai, D.J. Hargreaves, Water lubricated bearings, in: Green Tribology, Springer, Berlin Heidelberg, 2012 [Google Scholar]
 X. Wang, K. Koji, A. Koshi, A. Kohji, The effect of laser texturing of SiC surface on the critical load for the transition of water lubrication mode from hydrodynamic to mixed, Tribol Int 34 (2001) 703–711 [CrossRef] [Google Scholar]
 X.L. Wang, K. Kato, K. Adachi, K. Aizawa, Loads carrying capacity map for the surface texture design of SiC thrust bearing sliding in water, Tribol Int 36 (2003) 189–197 [CrossRef] [Google Scholar]
 D.L. Cabrera, N.H. Woolley, D.R. Allanson, Y.D. Tridimas, Film pressure distribution in waterlubricated rubber journal bearings, Proc Inst Mech Eng Part J: J Eng Tribol 219 (2005) 125–132 [CrossRef] [Google Scholar]
 X.L. Wu, Lubrication design handbook, Chemical Industry Press, Beijing, 2006 [Google Scholar]
 ANSYS, ANSYS 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 Fifth International Conference on Multiphase Flow, Yokohama, Japan, 2004 [Google Scholar]
 G. Zhou, J. Ma, Y. Quan, Fluid lubrication theory, Zhejiang University, Hangzhou, 1990 [Google Scholar]
 B. Jakobsen, L. Floberg, The rectangular plane pad bearing, Gumpert, 1958 [Google Scholar]
Cite this article as: X. Zhang, G. Gao, Z. Yin, Y. Wang, C. Gao, Numerical analysis and experimental research on load carrying capacity of waterlubricated tiltingpad thrust bearings, Mechanics & Industry 19, 201 (2018)
All Tables
All Figures
Fig. 1 Flow model of one sector pad and one groove of tiltingpad thrust bearing (not to scale). 

In the text 
Fig. 2 Boundary conditions. 

In the text 
Fig. 3 Pressure distribution for pads with different δ*, and B/L = 1.02, h_{2} = 10 µm, N = 3000 rpm. 

In the text 
Fig. 4 Maximum pressure, pmax, versus δ*, and B/L = 1.02, h_{2} = 10 µm, N = 3000 rpm. 

In the text 
Fig. 5 Pressure distribution for different load and rotational speed. 

In the text 
Fig. 6 Nondimensional pressure distribution at (a) the radius of maximum pressure (b) the medium radius of tiltingpad thrust bearing for different B/L, and δ^{*} = 1, h_{2 }= 10 µm, N = 3000 rpm. 

In the text 
Fig. 7 Nondimensional load carrying capacity, F_{z}*, versus(a) height of inclined plane, δ, and (b) height of inclined plane to minimum film thickness ratio, δ*, for different minimum film thicknesses, h_{2}, and B/L = 1.02. 

In the text 
Fig. 8 Nondimensional load carrying capacity, F_{z}*, versus height of inclined plane to minimum film thickness ratio, δ*, for pads with different B/L values, and h_{2} = 10 µm. 

In the text 
Fig. 9 The optimum δ^{*}versus pad widthtolength ratio B/L. 

In the text 
Fig. 10 Nondimensional radius of pivot, R_{s}^{*}, versus δ* for pads with different B/L values, and h_{2 }= 10 µm. 

In the text 
Fig. 11 Nondimensional angle of pivot, θ_{s}^{*}, versus δ*, for pads with different B/L values, and h_{2 }= 10 µm. 

In the text 
Fig. 12 Change of pivot position with δ*, for pads with different B/L values, and h_{2 }= 10 µm. 

In the text 
Fig. 13 Optimum position of pivot in tilting pad thrust bearing for different B/L. 

In the text 
Fig. 14 Nondimensional load carrying capacity for different pad widthtolength ratio. 

In the text 
Fig. 15 Waterlubricated bearing test rig. 

In the text 
Fig. 16 Schematic of waterlubricated bearing test rig: 1servo motor; 2coupling; 3rolling bearings; 4shaft; 5test shaft; 6 measuring rod for friction torque; 7friction torque force sensor; 8thrust bearing; 9thrust bearing pedestal; 10thrust ball bearing and journal bearing; 11ball seat; 12guide rail; 13force sensor; 14air cylinder. 

In the text 
Fig. 17 Measurement part of waterlubricated thrust bearing test rig. 

In the text 
Fig. 18 Structure of test bearing. 

In the text 
Fig. 19 Installation diagram of thrust pad. 

In the text 
Fig. 20 Waterlubricated thrust bearing. 

In the text 
Fig. 21 Measurement of thrust bearing film thickness. 

In the text 
Fig. 22 Friction coefficient versus specific pressure for N = 600 rpm and 900 rpm. 

In the text 
Fig. 23 Film thickness versus specific pressure for N = 600 rpm. 

In the text 
Fig. 24 Film thickness versus specific pressure for N = 900 rpm. 

In the text 
Fig. 25 Friction coefficient versus rotational speed for p_{s }= 0.286 MPa and 0.516 MPa. 

In the text 
Fig. 26 Film thickness versus rotational speed for p_{s }= 0.286 MPa. 

In the text 
Fig. 27 Film thickness versus rotational speed for p_{s }= 0.516 MPa. 

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.