Abstract

The bent-axis piston pump is the core component of electrohydrostatic actuators (EHA) in aerospace applications, and its wear of key friction interfaces is greatly related to the healthy operation of pumps. The leakage of the piston ring-cylinder bore interface (PRCB), as the important part of the return oil flow of the pump house that commonly assesses the wear of key friction interfaces in piston pumps, is changed with the rotation speed. Thus, the wear of key friction interfaces is usually inaccurate by using the leakage of PRCB. In order to obtain the relationship between the PRCB leakage and the rotation speed, an elastohydrodynamic lubrication model is proposed. First, the proposed model includes a minimum film thickness model of PRCB to analyze the dynamic change of oil film of PRCB when subject to the elastohydrodynamic lubrication. After that, a mathematical model of PRCB is induced by combining the minimum film thickness model with the flow equation, which helps produce the effects of the oil film on the leakage of PRCB. The proposed model is verified by numerical simulation and experiment. The results show that the leakage of PRCB has a negative effect on the return oil flow of the pump case in the range of rotation speed of 700–1300 r/min and discharge pressure of 10–20 MPa. Furthermore, the leakage of PRCB is proportional to the rotation speed, but the return oil flow of the pump case is decreased. The effects of rotation speed are enhanced under the high discharge pressure conditions.

1. Introduction

Bent-axis piston pumps are widely used in electrohydrostatic actuators (EHA) in aerospace applications due to the advantage of miniaturization and lightweight through high speed and high pressure. Therefore, the bent-axis piston pump greatly affects the healthy operation of EHA [13].

Performance degradation prediction and fault diagnosis are of great significance to ensure the safe and stable operation of mechanical equipment and reduce maintenance cost [4, 5]. In hydraulic pumps, the wear of key friction pairs is a major cause of the degradation of hydraulic pump performance [69]. As the wear increases with pump working time, the lubricating interface gap increases and the lubricating interface leakage increases, which eventually leads to the increase of the return oil flow rate of the pump case [10, 11]. Hence, the return oil flow rate is often used to evaluate the wear state of the key friction pairs in the piston pump [1214]. In engineering scenarios, however, the lubricating interface leakage is inevitably changed with the varying working conditions, such as the rotation speed, which affects the accurate evaluation of wear states.

The bent-axis piston pump is shown in Figure 1. The cylinder block moves synchronously with the drive shaft through the bevel gears, and the piston reciprocates in the cylinder bore. The piston ring with the spherical cross-section is usually installed on the end of the piston to reduce the leakage between the piston and the cylinder bore. There are two lubricating interfaces, that is, the valve plate-cylinder block interface (VPCB) and the piston ring-cylinder bore interface (PRCB). As one of the two lubricating interface leakages, PRCB leakage directly affects the return oil flow rate in bent-axis piston pump.

The lubrication interface gap, as the main parameter affecting PRCB leakage, is influenced by multifactor interactions. Pelosi et al. developed the fully coupled multibody dynamics model to capture the complex fluid-structure interaction phenomena affecting the nonisothermal fluid film conditions and analyzed the influence of the squeeze film effect and elastic deformation on the piston-cylinder interface [15]. In order to improve energy efficiency of piston pump, Song et al. presented a numerical method to analyze influence of temperature on lubricant film characteristics of the piston-cylinder interface [16]. Qian et al. presented a nonisothermal fluid-structure mathematical model considering the piston eccentricity, elastic deformations, and the fluid physical properties to analyze the piston-cylinder bore interface leakage and the fluid pressure distribution [17]. Zhang et al. established the mathematical model of the oil film thickness by using the cosine theorem in the cross section of the piston and studied the piston-cylinder interface leakage and pressure field under the condition of ultrahigh pressure [18]. Nie et al. established a parameterized elastohydrodynamic lubrication model and studied the influence of the elastohydrodynamic behavior, viscosity temperature effect, and deep-sea environmental pressure on piston-cylinder interface leakage in seawater hydraulic axial piston pump [19]. Cho et al. analyzed the film thickness, the pressure distribution, and the friction force of piston ring-cylinder bore interface under elastic deformation of plunger ring in bent-axis piston pump with the tapered piston [20]. Kumar established the theoretical model of the piston ring-piston bore pair in the bent-axis piston pump with rectangular cross-section of piston ring and analyzed the leakage flow characteristics between piston ring and cylinder bore under different working conditions through numerical simulation and test [21]. Manhartsgruber analyzed the dynamics of the leakage of the bent-axis piston pump with the tapered piston and revealed the influence of the piston rotation inside the cylinder bore on the pressure fluctuation under low speed and low pressure condition by experiment [22]. Zhang et al. analyzed the effect of centrifugal force and reciprocation inertial force on the leakage between piston and cylinder in the bent-axis piston pump under high speed and show the leakage of piston-cylinder pair increases obviously at high speed [23].

From the literature review, the previous studies mentioned mainly focused on the influence of the piston eccentricity, elastic deformation, and thermal effect on the piston-cylinder bore interface leakage in the swash-plate piston pump. Some researchers studied the effects of the elastic deformation and micromotion of piston ring on the oil film thickness, pressure field, and friction force of PRCB in the bent-axis piston pump under certain rotation speed conditions. However, the variable speed and wide speed are the typical working conditions of bent-axis piston pump in EHA, and the influence of rotation speed on PRCB leakage considering elastohydrodynamic lubrication effect has not been fully elucidated.

Rotation speed and its fluctuation are very meaningful monitoring parameters, and the leakage analysis of the friction pairs is also related to the rotation speed. Our previous research has proved that the operating state information of the swash-plate piston machine can be extracted from the speed fluctuation signal [2427]. In this study, an elastohydrodynamic lubrication model is proposed to discuss the influence of rotation speed on the leakage of PRCB. The proposed model includes a minimum oil film model of PRCB to analyze the dynamic change of minimum oil film of PRCB under the elastohydrodynamic lubrication. Combining the minimum oil film model with the flow equation, a mathematical model of PRCB is induced to help produce the effects of the oil film on the leakage of PRCB. The accuracy of the proposed model is verified by simulation and experiment, which provides guidance on the piston pump wear diagnosis based on the leakage. The layout of this paper is organized as follows. A theoretical model of PRBC is established and explained in detail in Section 2. The simulation results are discussed in Section 3. The experimental results are presented in Section 4, followed by the conclusions in Section 5.

2. Elastic Deformation and Leakage Mechanism of PRBC

2.1. Elastohydrodynamic Lubrication Equations

In order to reduce the leakage of PRCB, a ring groove is usually opened on the spherical surface of the large end of the piston, and a piston ring is installed in the ring groove. The boundary conditions between the piston ring and the cylinder bore are different because of the change of the angle position of the piston in the suction and discharge region, as shown in Figure 2.

Because of the elastic deformation of the contact surface in the piston ring at relative motion, the gap of PRCB is changed, which affects the leakage. Therefore, it is necessary to access elastohydrodynamic lubrication (EHL) to analyze the leakage between the piston ring and the cylinder bore. The analysis model between the piston ring and the cylinder bore is shown in Figure 3. V represents the sliding velocity of the piston and the piston ring relative to the cylinder bore, R is the radius of the piston ring section, W is the unit load on the piston ring. Figures 3(a) and 3(b) show the pressure distribution of lubricating oil film between the piston ring and the cylinder bore, and its typical characteristics are the primary peak pressure and the secondary peak pressure. P1 represents the initial pressure of the oil film, and it equals the pressure in the displacement chamber, and P2 represents the pressure at the end of the film, and it equals the pressure in the pump case.

The assumptions for the EHL analysis between the piston ring and the cylinder bore are listed as follows:(1)The lubricating pressure remains unchanged along the direction of film thickness(2)The lubricating oil film remains unchanged along the circumferential direction of piston bore(3)The fluid inertia force is ignored(4)The compressibility and viscosity temperature characteristics of the lubricating oil film are ignored(5)The only concern is the velocity gradient along the piston ring direction of motion

The Reynolds equation for one-dimensional isothermal elastohydrodynamic lubrication is shown as follows [28]:where is oil film density, is the oil film thickness, is the pressure distribution of lubricating oil film, p=p1 at x=xa, p=p2, ∂p/∂x= 0 at x=xb, us is the average oil film velocity, us=v/2, and η is the viscosity oil film.

The velocity of the piston relative to the cylinder bore can be expressed as follows:where is the nominal pitch circle radius of piston balls on the driving flange, is the angular velocity of the driving shaft, is the angular position of piston relative to the top dead point, γ is the angle between the axis of the cylinder and axis of the driving shaft.

The film thickness is expressed as follows [29]:where , μ1, E is the equivalent Young’s modulus, μ1 and μ2 are Poisson's ratio of the piston ring and the cylinder, E1 and E2 are Young’s modulus of the piston ring and the cylinder, h0 is the film thickness at x=0, and C is the integral calculus constant.

The viscosity coefficient and density coefficient of lubricating oil are expressed as follows [30]:

The density coefficient of lubricating oil is expressed as follows [31]:where η0 is the viscosity of oil under atmospheric pressure, p0 and z are the pressure-viscosity coefficients, and ρ0 is the density of oil under atmospheric pressure.

The unit load on piston ring is expressed as follows:where is the starting point of lubricating film in the analysis and is the ending point.

2.2. PRCB Leakage Model

The PRCB leakage is composed of two parts. One is the Poiseuille flow by the pressure difference between the displacement chamber and the pump case, and the other is the Couette flow by the piston ring relative to the cylinder bore motion and the hydraulic oil viscosity. As the most important parameter of PRCB, the minimum oil film varies dynamically with the piston velocity due to the elastohydrodynamic lubrication effect, which will change the proportion of the Poiseuille flow and the Couette flow in the gap and then affect the leakage flow of PRCB. The dynamic model of the leakage flow of PRCB is established by combining the minimum oil film of PRCB with the flow equation, which is the theoretical basis for the subsequent study of the leakage flow of the plunger pair under the variable speed condition.

The instantaneous leakage flow rate between the piston ring and the cylinder bore is expressed as follows [23]:where is the diameter of piston, is the gap between piston ring and cylinder bore, is the pressure difference between displacement chamber and pump case, is the oil viscosity, is the width of piston ring, is the eccentricity of the piston ring, and is the velocity of motion of piston ring in the cylinder bore.

The average leakage flow rate of PRCB is expressed aswhere is number of the pistons and is the time during one rotating speed of the bent-axis piston pump.

3. Simulation Results

3.1. Simulation Parameter Setting

The finite difference method is used to solve the theoretical model. In order to improve the calculation accuracy, the central difference scheme is used to discretize the Reynolds equation. The flow chart of the numerical calculation program for the elastohydrodynamic lubrication problem of PRCB is shown in Figure 4.

According to the given working parameters and the initial values of oil film thickness, distribution pressure, viscosity, density and other parameters, the Reynolds equation is calculated to obtain the new oil film distribution pressure, and then the elastic deformation, viscosity, density, and other parameters of the piston ring are calculated under the obtained pressure.

The Gauss–Seidel iterative method is used to solve the algebraic equations of Reynolds equation discretization, and the converging condition is expressed as follows:where is the relaxation factor and α=0.8 is chosen to obtain a fast convergence, and are the pressure values at steps k and k + 1, and is the temporary value at step k + 1.

The whole process is repeated until the error tolerance of the Reynolds equation at each node is less than the relative precision requirement ε = 1 × 10−5; the error tolerance of the Reynolds equation ε is expressed as follows:

MATLAB platform is used to calculate the leakage flow rate between piston ring and cylinder bore. In order to compare the characteristics of the leakage with different operating conditions, the parameters of simulation calculation process are set in Table 1.

3.2. Analysis of Simulation Results
3.2.1. The Pressure Distribution of PRCB

The pressure distribution of PRCB in the discharge region is shown in Figure 5.

As shown in Figure 5, the pressure of oil film changes in the lubrication zone at a certain angle position, and there is the maximum pressure value (the local maximum pressure). The local maximum pressure varies with the angular position of the discharge region, and there is also the maximum pressure (the global maximum pressure). As a whole, the local maximum pressure in the discharge region decreases first and then increases with the increase of angle position, and the global maximum pressure is 146 MPa at angle position  = 4.7°, while the angle position is near the top dead point 0° and the bottom dead point 180°; the local maximum pressure changes dynamically again. The reason is that the piston velocity is small near the top and bottom dead points, and the local maximum pressure is determined by the primary peak pressure, as shown in Figure 3(b). However, the local maximum pressure is determined by the secondary peak pressure in other region, and the secondary peak pressure decreases with the increase of the piston velocity.

The local maximum pressure distribution of PRBC in the discharge region at different speeds is shown in Figure 6.

As shown in Figure 6, when the speed increases from 700 r/min to 1300 r/min with a step of 200 r/min, the global maximum pressure is basically not affected by the change of speed, which is kept at 146 MPa, but the angular position of the global maximum pressure decreases from 4.7° to 2.5°.

The local maximum pressure distribution of PRBC in the discharge region at different discharge pressure is indicated in Figure 7.

As indicated in Figure 7, when the discharge pressure increases from 10 MPa to 20 MPa at n = 700 r/min, the local maximum pressure becomes larger as a whole, the global maximum pressure increases from 146 MPa to 205 MPa, and the angular position of the global maximum pressure increases from 4.7° to 15.3°.

The rotation speed has a weak effect on the value of the global maximum pressure, which has a significant effect on its the angular position. The discharge pressure has important influence on the value and the angle position of the global maximum pressure. The reason is that the global maximum pressure of PRBC in the discharge area is determined by the secondary peak pressure of the lubricating oil film. The secondary peak pressure is positively correlated with the discharge pressure. The angle position of the secondary peak pressure is affected by the rotation speed and the discharge pressure. The angular position of the secondary peak is positively correlated with the discharge pressure, and it is negatively correlated with the rotation speed.

3.2.2. The Minimum Thickness of PRCB

The minimum thickness of PRCB at different speed is displayed in Figure 8.

As displayed in Figure 8, the minimum film thickness increases first and then decreases with the increase of angle position in the discharge region and the suction region. However, the minimum film thickness of PRCB in the discharge region is obviously larger than that in the suction region as a whole. When the speed increases from 700 r/min to 1300 r/min with a step of 100 r/min, the minimum thickness increases in the discharge region and the suction region, and its increase in the discharge region is greater.

The effect of discharge pressure on the minimum thickness of PRCB is shown in Figure 9.

As shown in Figure 9, the minimum thickness of PRBC in the discharge region and the suction region increases first and then decreases with the increase of angular position, but the former is larger than the latter as a whole. When the discharge pressure increases gradually from 10 MPa to 20 MPa with a step of 5 MPa, the minimum oil film increases obviously in the oil drainage region, but only slightly in the suction region.

The reason is that the minimum thickness is determined by the velocity and the load of the piston ring. When the rotation speed and the discharge pressure are constant, the velocity of piston ring increases first and then decreases from the top dead point (TDC) to the bottom dead point (BDC) in the drainage region and is the same change trend from BDC to TDC in the suction region, so the change law of minimum thickness is the same as the velocity of piston ring. The discharge pressure mainly affects the load of the piston ring in the discharge region, and the discharge pressure increases, and the load of the piston ring becomes larger, which leads to the increase of the elastic deformation of the piston ring, so the minimum thickness in the discharge region is greater than that in the suction region.

The mean value of minimum thickness of PRCB is shown in Figure 10.

As shown in Figure 10, it is found that the mean value of the minimum thickness in the discharge region and the suction region increases with the increasing rotation speed. When the discharge pressure is 20 MPa, with the rotation speed increasing from 700 r/min to 1300 r/min, the mean value of the minimum thickness changes the most, and the variation in the mean value in the discharge region is 2.87 um, while the variation in the mean value in the suction region is 0.28 um.

The speed has obvious influence on the mean value of minimum thickness in the discharge region and the suction region, and the mean increases with the increase of the speed. The discharge pressure has obvious influence on the mean in the discharge area, and the mean value increases with the increase of discharge pressure, but the discharge pressure has a weak effect on the mean in the suction area. The mean of minimum thickness in the discharge area is obviously higher than that in the suction area.

As the speed increases from 700 r/min, the variation of the mean minimum thickness increases with the increase of the variation of the speed in the discharge region and suction region, especially high discharge pressure, as described in Figure 11.

In this picture, the variation of mean minimum thickness in the discharge region is greater. Moreover, the variation in the discharge region is obviously higher than that in the suction region. The maximum of variation in the discharge region is 2.87 um, and the maximum of variation in the suction region is 0.28 um, and the former is 10 times the latter at p = 20 MPa.

3.2.3. The Leakage Flow Rate of PRCB per Piston

The leakage of PRCB is caused by the Poiseuille flow and the Couette flow. The leakage from the displacement chamber to the pump case is positive, whereas the leakage is negative. When the piston is in the suction region, the Poiseuille flow is not considered because the pressure difference between the displacement chamber and the pump case is very small.

The effect of rotation speed on the instantaneous flow of PRCB per piston is presented in Figure 12.

As presented in Figure 12, the Poiseuille flow first increases and then decreases with the increase of angular position in the discharge region, and the Poiseuille flow in the discharge pressure increases with the increase of speed, as shown in Figure 12(a). The Couette flow is negative in the discharge pressure region and is positive in the suction region. The value of the Couette flow first increases and then decreases with the increase of angular position, but the value of the Couette flow in the discharge region is significantly larger. When the speed increases from 700 r/min to 1300 r/min, the leakage flow in the discharge region and suction region increases, but the variation of the former is larger, as shown in Figure 12(b). The total flow rate of PRCB is negative in the discharge area and positive in the suction area. The value of the total leakage increases with the increase of rotation speed, but the value in the discharge area is obviously larger than that in the suction region as shown in Figure 12(c).

The effect of the discharge pressure on the instantaneous flow rate of PRCB per piston is shown in Figure 13.

As shown in Figure 13, the Poiseuille flow first increases and then decreases with the increase of angular position in the discharge region and increases with the increase of speed as a whole, as shown in Figure 13(a). The Couette flow is negative in the discharge pressure region and is positive in the suction region. The value of the Couette flow first increases and then decreases with the increase of angular position, but the value of the Couette flow in the discharge region is significantly larger. When the speed increases from 700 r/min to 1300 r/min, the leakage flow in the discharge region and suction region increases, but the variation of the former is larger, as shown in Figure 13(b). The total flow rate is negative in the discharge region and is positive in the suction region. The total flow rate increases with the increase of the speed and the discharge pressure, but the value in the discharge region is obviously larger than that in the suction region, as shown in Figure 13(c).

The reason is that the Poiseuille flow is determined by the minimum thickness and the pressure difference between the displacement chamber and the pump case. The pressure difference and the minimum thickness increase with the increase of discharge pressure. Once the discharge pressure is determined, the pressure difference will basically remain unchanged. Therefore, the change trend of the Poiseuille flow consistent with the change trend of the minimum thickness in the discharge region under different the discharge pressure. The Couette flow is determined by the velocity of the piston ring and the minimum thickness between the piston ring and the cylinder, and the positive and negative of the Couette flow are determined by the direction of piston velocity. The leakage flow is caused by the Poiseuille flow and the Couette flow, and the two flows oppose in the discharge pressure region and support each other depending on in the suction region. According to Figures 12 and 13, the Couette flow plays a leading role in leakage flow of PRBC, especially the Couette flow in the discharge region.

The average leakage flow rate of PRCB per piston at different speeds is shown in Figure 14.

As shown in Figure 14, it is found that the average leakage flow is negative, which indicates the flow rate flows from the pump case to the displacement chamber. When the speed is constant, the average leakage flow rate increases with the increase of the discharge pressure. When the discharge pressure is constant, the average leakage flow rate increases with the increase of the speed. If the discharge pressure is 20 MPa, the mean increased from 27.9 ml/min to 75.6 ml/min with the increase of the speed from 700 r/min to 1300 r/min, which increased by 171 percent.

4. Experimental Results

4.1. Test Rig Description

The leakage flow rate of PRCB is difficult to measure directly. According to the previous works [21, 23], the leakage flow rate of PRCB is an important part of the return oil flow rate. Therefore, the return oil flow rate is measured to investigate the effects of rotation speed on the leakage flow of PRCB. The bent-axis piston pump experiments are implemented on a test rig, as shown in Figure 15.

As shown in Figure 15, the test rig consists of a bent-axis piston pump, a motor, a speed sensor, a pressure senor, a temperature senor, and a relief valve and so on. The bent-axis piston pump is continuously adjusted from 0 to 2000 r/min by the motor. The return oil flow rate is measured by the flow sensor which is installed on the outlet of the pump case in the bent-axis piston pump. The discharge pressure is controlled by adjusting the relief valve. The signals obtained from the main sensors are collected using NI acquisition card and LabVIEW software. The detailed descriptions of the main parameters of the test rig are shown in Table 2.

4.2. Data Processing

The data acquisition diagram of return oil flow rate in the bent-axis piston pump is presented in Figure 16. The sampling frequency was 1 Hz, and the time interval was 35 s under different working conditions. The empirical mode decomposition (EMD) was used to deal with the return oil flow rate under p = 20 MPa n = 1300 r/min to obtain its intrinsic mode functions (IMF), as shown in Figure 16(a). According to the noise reduction and smoothing of the collected the return oil flow rate data, the return oil flow rate diagram under different working conditions is shown in Figure 16(b). It can be seen from this figure that the overall trend of the return oil flow rate decreases with the increase of rotation speed.

4.3. Results Analysis

The simulated flow rate of PRCB and the measured the return oil flow rate in the bent-axis piston pump is shown in Figure 17.

In Figure 17, it is obvious that the simulated total flow rate of PRCB is negative, and its value increases with the increase of speed and discharge pressure, but the measured the return oil flow rate decreases with the increase of speed and increases with the increase of discharge pressure. The reason is that the measured the return oil flow rate is composed of the leakage of VPCB and the leakage of PRCB. The leakage of VPCB is the main part of the measured the return oil flow rate, and its values increases with the increase of the discharge pressure and is less affected by rotation speed. But the flow rate of PRCB is negative, which leads to the decrease of the measured return oil flow rate.

The variations of the simulated flow of PRCB and the measured return flow rate relative to the flow rate under the same discharge pressure and 700 r/min is shown in Figure 18.

In this picture, it is obvious that the variations of the simulated flow of PRCB increase with the increase of rotation speed, indicating that the simulated flow of PRCB increases with the increase of speed. The variation of the measured return flow rate is negative, and it decreases with the increase of rotation speed, indicating that the measured return flow rate decreases with the increase of rotation speed. Moreover, the trend of the two variations with rotation speed is basically the same at the different discharge pressure. However, there is a certain difference between the two; its maximum is 44 ml/min at n = 1300 r/min and p = 20 MPa. The main reason is that the measured result is by the leakage flow of PRCB but also by the leakage flow of VPCB.

5. Conclusions

The model is suitable for analyzing the dynamic change of the minimum oil film of PRCB under the action of elastohydrodynamic lubrication. Combining the minimum oil film model with the flow equation, a mathematical model of PRCB is introduced to help produce the effect of the oil film on the leakage of PRCB. The effect of rotation speed on the leakage flow of PRCB is analyzed through simulation and experiment, and the following conclusions are obtained.(1)The lubrication of PRCB in the discharge of oil segment belongs to the elastohydrodynamic lubrication, and the angle position of the maximum oil film pressure is close to the top dead point 0° with the increase of rotation speed, but its maximum pressure value is basically not affected by the change of the rotation speed.(2)The minimum thickness of PRCB increases with the increase of rotation speed, and its value in the discharge of oil segment is greater than that in the suction of oil segment, and the former varies more with the rotation speed.(3)The average flow rate of PRCB is from the pump case to the displacement chamber, which reduces the return oil flow rate in the pump case. The return oil flow rate decreases with the increase of rotation speed, and the effect of rotation speed on the average flow rate of PRCB and the return oil flow rate is also remarkable under the high discharge pressure.

The return oil flow rate in the bent-axis piston pump should be the synthesis of the two friction pairs, which are related to the rotation speed. In this paper, the relationship between the elastohydrodynamic lubrication characteristics of PRCB and the rotation speed is analyzed. The relationship between leakage of the valve plate-cylinder block interface and rotation speed in the bent-axis piston pump will be the focus of further work.

Data Availability

The data that support the findings of this study are available on request from the corresponding author. These data were derived from the test rig built by the authors.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by the National Natural Science Foundation of China (no. 51675399), Special Research Project of Education Department of Shaanxi Province (no. 18JK0408), and the Science Foundation of Xi’an Aeronautical University (no. 2020KY0223).