1. Introduction
Compared with belts, chains and other transmission devices, gear drives are the most widely used in a variety of engines, such as those in automobiles, ships and spacecraft, due to their high efficiency, compact structure and smooth operation [
1]. As a major problem in gear drives, lubrication exerts a significant influence on gear drive performance [
2]. Good lubrication can improve the life of gears, reduce noise, and increase meshing efficiency [
3].
A review of the published literature shows that the field of gear lubrication has two major research focuses. The first focuses on the state of the lubricating oil film on the meshing tooth surfaces, including the model analysis of EHL, in order to establish the nonlinear relationships among oil film thickness, pressure and density [
4]; numerical methods for solving the EHL [
5]; and measurement of and experimentation on the contact area [
6]. The second group focuses on the use of commercial computational fluid dynamics (CFD) software to simulate the flow state of lubricating oil during the gear meshing process and visualise the flow field of lubricating oil [
7].
Gears operate under heavy loads and high pair meshing, and contact pressure on tooth surfaces is extremely high. The characteristics of lubricating oil change, and meshing tooth surfaces deform, when under high pressure. The gears are in a state of EHL [
4]. The state of the lubricating oil film between the contacting tooth surfaces in EHL is described using the Reynolds equation. This equation reasonably explains the formation mechanism of oil film; this mechanism is the theoretical basis of EHL [
8]. The pressure–viscosity effect of lubricating oil and the elastic deformation of contact surfaces are considered in EHL. Thus, the EHL problem is described using a series of nonlinear coupled equations that are extremely difficult to solve. The earliest method used to solve EHL was a direct iteration method that produced considerable errors in the actual results. The inverse, multigrid [
9] and finite element [
10] methods were subsequently introduced, and their results can accurately describe the shape of the lubricating oil film and the pressure distribution in linear contact EHL.
Recent studies have used the Reynolds equation to solve the general lubrication problem in non-Newtonian fluids by regarding lubricating oil flow as the superposition of Poiseuille and Couette flows. The velocity of lubricating oil flow in this equation was solved using the characteristics of Poiseuille and Couette flows, and the other characteristics of lubricating oil film were solved on the basis of Newtonian fluids [
11]. The traditional Integral Approach, the Single Domain Full-System Approach and the Double Domain Full-System Approach have been used to solve the EHL problem [
12]. The comparisons of velocity and oil film thickness between this method and the conventional approach established the feasibility and validity of the method and illustrated its simplicity of calculation. The influence of the surface roughness of gear teeth on the lubricating oil film was significant [
13]. Thus, the modified Reynolds equation in Patir and Cheng [
14] and the surface deformation and statistical elasto-plastic asperity micro-contact model in Zhao et al. [
15] were solved using a systematic approach under steady-state and isothermal operating conditions. The central oil film thickness formula, minimum oil film thickness and asperity-to-total load ratio were obtained via the regression of simulation data [
13]. To examine the influence of roughness on EHL, the rough surface was regarded as a rough plane surface with uniformly distributed cylindrical asperities that were evenly truncated on top, and the inlet of the contact area was analysed. Several valuable conclusions, including the reduction of load capacity in the oil film and the increase in friction coefficient due to roughness, particularly at low speeds and under heavy loads, were drawn [
16]. The load transmitted by a gear tooth during the meshing process may fluctuate significantly due to time-varying meshing stiffness and torque fluctuation. The influence of dynamic load on the minimum oil film thickness, pressure distribution and temperature distribution were analysed on the basis of the quasi-static tooth force solved via finite element contact analysis [
17]. The results showed that a dynamic load will reduce the thickness of the oil film, increase its pressure and temperature and destroy the EHL of the contact surface. A tribo-dynamic model for spur gear pairs combining dynamic loading with EHL theory was proposed in 2016 [
18]; this model can predict the contact force and velocity of dynamic teeth. A method that combined the Runge Kutta and multigrid methods was used to solve the aforementioned tribo-dynamic model, and the results were compared with those of the model under quasi-static conditions. Starved lubrication conditions will occur when an incomplete oil film is formed between contact surfaces. To investigate the influences of starved lubrication on minimum film thickness, maximum temperature and friction coefficient, a thermally starved EHL model was proposed [
19]. This previous research found that the minimum oil film thickness was thinner, and the friction coefficient and maximum temperature were higher, in starved lubrication than those in EHL.
CFD is a new branch of fluid mechanics that has been widely used in gear lubrication. Liu et al. [
20] built a CFD model based on the finite volume method to investigate the lubricating oil distribution and the churning power losses of a gearbox. The results were obtained by considering the 3D finite volume simulation model as a two-phase flow, and the results were compared with the measurements of churning losses and high-speed camera recordings. The comparisons exhibited good agreement in predicting lubricating oil distribution and churning power losses. Hu et al. [
21] presented a CFD model combined with turbulence to investigate the lubricating oil distribution and churning power losses of a moving gearbox. The influences of gear rotation, immersion depth, gearbox motion frequency and amplitude on the results were determined.
On the basis of the preceding statement, current studies have explored effective methods for establishing a lubrication analysis model that can realistically and accurately describe the contact state among tooth, the characteristics of oil film and the distribution of lubricating oil. However, these studies on gear lubrication are based on the condition that the lubricating oil supply and the lubrication state of the gear are known. A method that can analyse the lubricating oil demand of gears under specific operating conditions could help significantly in the design of gear lubrication systems.
In the current work, new calculation and analysis models of the lubricating oil demand of spur gears are proposed by considering lubricating oil demand as the superposition of oil for lubrication and cooling. This study aims to understand variations in oil for lubrication and cooling during a mesh cycle under different parameters and working conditions. The calculation models of lubricating oil demand are established in accordance with the meshing of spur gears and EHL theories. The influences of modulus, number of teeth, input speed and input torque on oil for lubrication and cooling are simulated. This study provides guidance for the appropriate design of gear lubrication systems.
3. Lubricating Oil Demand during Meshing Process
The curvature radius and load of a tooth surface are continuously changing during the gear meshing process. Thus, the change in lubricating oil demand during the meshing process should be analysed. The parameters of gears, lubricating oil and working conditions are provided in
Table 1.
The solution for lubricating oil demand can be attributed to solving complicated nonlinear equations. A compound iterative method is used to obtain the demand of lubricating oil.
The central film thickness and entrainment velocity of each meshing tooth are compared during the meshing process, as shown in
Figure 2. The two vertical coordinates shown in
Figure 2 represent the central film thickness (μm) and entrainment velocity (m/s), respectively. The abscissa denotes the coordinates on the meshing line.
Figure 3 shows the oil volume for the lubrication of each meshing tooth during the meshing process and the total oil volume for the lubrication of a gear.
Meshing teeth 1 and 2 mesh simultaneously in the first two tooth meshing area. The entrainment velocity of meshing teeth 1 and 2 in the first two tooth meshing area increases continuously with gear meshing. An oil film easily forms at high entrainment velocities. Thus, the oil film thickness of the two meshing teeth in the meshing area of the first two teeth increases with gear meshing. A positive correlation exists between oil film thickness and oil volume for lubrication. The oil volume for the lubrication of meshing teeth increases along the meshing line, reaching the maximum value at the end of the meshing area of the first two teeth. The total oil volume for the lubrication of a gear is the sum of the oil volume for the lubrication of meshing teeth 1 and 2.
As the gear rotates, meshing tooth 2 moves gradually away from the meshing. Thus, the gear is in the single tooth meshing area. The sudden increase in load of meshing tooth 1 reduces the oil film thickness. Therefore, the oil volume for the lubrication of meshing tooth 1 plummets at the intersection of the meshing areas of the first two teeth and the single tooth. The total oil volume for the lubrication of the gear is the same as the oil volume for the lubrication of meshing tooth 1. Oil film thickness and oil volume for the lubrication of meshing tooth 1 increase linearly in the single tooth meshing area.
Meshing tooth 3 begins to mesh at the end of the single tooth meshing area, and the gear is located in the second two tooth meshing area. Entrainment velocity, oil film thickness and oil volume for the lubrication of meshing tooth 1 are larger than those for meshing tooth 3, and they increase continuously with gear meshing. The oil volume for the lubrication of two meshing teeth reaches the maximum value at the end of the second two tooth meshing area. The total oil volume for the lubrication of the gear is the sum of the oil volume for the lubrication of meshing tooth 1 and 3.
The meshing process involves the continuous and repeated movement of a meshing point on a meshing line. In a meshing circle, two equal maxima for the oil volume for lubrication exist at the ends of the first two tooth meshing and the second two tooth meshing areas.
The analysis of the second section shows that the relative slide velocity and friction coefficient play key roles in the calculation of the oil volume for cooling. Therefore, the relative slide velocity and friction coefficient are compared during the meshing process, as shown in
Figure 4. The two vertical coordinates that represent the relative slide velocity and friction coefficient are established in this figure. The abscissa denotes the coordinates along the meshing line. The area selected by the wireframe is magnified 10 times for clear visualization.
Figure 5 illustrates the oil volume for the cooling of each teeth and the total oil volume for the cooling of a gear. The image is partially enlarged.
Figure 4 shows that the friction coefficient at the pitch point is equal to zero and that of the tooth surface increases as the distance between the meshing and pitch points increases, because the friction force between meshing tooth surfaces is the same as the shear force of the oil film in the EHL of a smooth surface. The friction coefficient increases as the relative slide velocity increases and oil film thickness decreases. The friction coefficient in the beginning of the meshing process is greater than that at the end of meshing process because the oil film thickness and relative slide speed at the end of the meshing process help reduce friction.
The friction coefficient of meshing tooth 1 is at its maximum value at the beginning of the first two tooth meshing area, resulting in the maximum value of oil volume for cooling, as shown in
Figure 5. This finding is attributed to the relative slide velocity being at its maximum value and the oil film thickness being at its minimum value. In the first two tooth meshing area, the total oil volume for the cooling of the gear is the sum of the oil volume for cooling of meshing tooth 1 and 2. The relative slide velocity, oil film thickness and oil volume for the cooling of meshing tooth 1 decrease with gear rotation, and those of meshing tooth 2 increase with gear rotation, albeit with a small increment, and the total oil for cooling continuously decreases.
When the gear is in the single tooth meshing area, meshing tooth 1 meshes separately. The friction coefficient and oil volume for the cooling of meshing tooth 1 increase suddenly at the beginning of the single tooth meshing area due to the decrease in oil film thickness. The oil volume for the cooling of the gear is equal to that of meshing tooth 1, but less than that of two tooth meshing area.
Meshing tooth 3 starts to mesh at beginning of the second two tooth meshing area. The maximum friction coefficient and oil volume for the cooling of meshing tooth 3 occur at this point. The total oil volume for the cooling of the gear is the sum of the oil volume for the cooling of meshing teeth 1 and 3. As the gear meshes, the relative slide velocity, oil film thickness and oil volume for the cooling of meshing tooth 3 decrease, but those of meshing tooth 1 increase. The total oil for the cooling of the gear still decreases because of the small increment of oil for the cooling of meshing tooth 1. Oil volume for the cooling of the gear has two maxima, which are located at the beginning of the first two tooth meshing and the second two tooth meshing areas.
4. Effect of Gear Parameters
Under the same working conditions, various parameters of the gear pair, such as number of teeth, modulus and transmission ratio, will lead to different contact states between tooth surfaces that will further affect the oil volumes for the lubrication and cooling of the gear.
4.1. Effect of Modulus
Modulus is a basic parameter of the gear that determines its size and exerts a significant impact on transmission performance. Moduli are set to 3, 4, 5, 6 and 8, and the other parameters listed in
Table 1 remain the same. Several lubrication simulations are conducted to analyse the influence of modulus on lubricating oil demand.
The variation in oil volume for the lubrication of the gear under different moduli during the meshing process is obtained, as shown in
Figure 6a.
This figure indicates that the length of the actual meshing line, the width of the first two teeth area, the width of the single tooth area and the width of the second two tooth area increase as the modulus of the gear increases. In addition, the oil volume for the lubrication of the gear increases with increasing modulus of the gear because the diameter of the gear increases with increasing modulus. The entrainment velocity and oil film thickness at the same rotational speed increase as the modulus of the gear increases. Notably, the gaps of the oil volume for the lubrication between the two tooth meshing and the single tooth meshing areas increase as the modulus of the gear increases.
The preceding analysis demonstrates that the friction coefficient is the primary influencing factor of the oil volume for cooling. The variations in the friction coefficient at different moduli with gear rotation are obtained, as shown in
Figure 6b, including a partially enlarged image. The friction coefficient decreases as the modulus increases, because the radius of the tooth surface curvature increases as the modulus of the gear increases and the friction coefficient is inversely proportional to the radius of the curvature.
A series of simulations is performed to examine the influence of different moduli of the gear on the oil volume for cooling, and an enlarged image of the single meshing area is shown in
Figure 6c. The oil volume for cooling decreases as the modulus of the gear increases due to the friction coefficient.
4.2. Effect of Number of Teeth
Number of teeth is an important parameter of gears that significantly affects the lubrication state between tooth surfaces during the meshing process and the performance of a gear transmission device. This section focuses on the influence of different numbers of teeth on the oil volumes for lubrication and cooling. The numbers of teeth in the pinion are set to 20, 23, 25 and 30. The other gear parameters listed in
Table 1 remain the same. A series of calculations is used to analyse the influence of different numbers of teeth on the oil volume for lubrication during the meshing process, as shown in
Figure 7a.
This figure indicates that the width of the single tooth meshing area decreases as the number of teeth increases, and the oil volume for the lubrication of the gear increases as the number of teeth increases because the size of the gear increases as the number of teeth increases. Consequently, the entrainment velocity and oil film thickness of all the meshing teeth increase at the same rotational speed. Notably, the difference in the oil volume for cooling increases as the number of teeth increases.
To compare the oil volume for cooling at different numbers of teeth, the variations in the friction coefficient that play a crucial role in the oil volume for cooling are obtained first, as shown in
Figure 7b, which includes a partially enlarged image.
Figure 7b shows that the friction coefficient decreases as the number of teeth increases because oil film thickness increases as the number of teeth increases. A series of simulations is performed to examine the influence of different numbers of teeth, as shown in
Figure 7c, and the enlarged image of the single meshing area is presented.
Figure 7c shows the decrease in the oil volume for cooling as the number of teeth increases. This behaviour is affected by the friction coefficient.
4.3. Effect of Transmission Ratio
In this section, the transmission ratio of a gear is used as the variable to examine the influence of transmission ratio on the oil volume for lubrication and cooling. The transmission ratios are set to 1, 2.95, 4.7 and 9.3. The other gear parameters are identical to those described in
Section 3. The variations in the oil volume for lubrication during the meshing process at different transmission ratios are obtained using a series of simulations, as shown in
Figure 8a.
This figure indicates that the oil volumes for lubrication in the first two tooth meshing and the second two tooth meshing areas are invariable, and their values are identical when the transmission ratio of the gear is 1. The oil volume for lubrication in the single tooth meshing area is half that in the two tooth meshing area, because the entrainment velocity and oil film thickness of each meshing tooth during the meshing process are invariable. Moreover, the distance between the start of the meshing and pitch point increases as the transmission ratio increases and the width of the single meshing area decreases. The rangeability of the oil volume for lubrication along the mesh line (i.e., the slope of the curve of the oil volume for lubrication) increases as the transmission ratio increases. The oil volume for lubrication at the beginning of meshing decreases and the oil volume for lubrication at the end of meshing increases as the transmission ratio increases, because a similar variation trend occurs in entrainment velocity and oil film thickness at different transmission ratios.
A series of simulations used to examine the influence of different transmission ratios on the friction coefficient during the meshing process and the variations in the oil volume for cooling at different transmission ratios during the meshing process are shown in
Figure 8b,c, respectively. These figures contain partially enlarged images.
Figure 8b indicates that the friction coefficient curve is symmetric around the pitch point when the transmission ratio of the gear is 1 and the friction coefficients at the start and end of meshing are consistent. The friction coefficient during the meshing process decreases as the transmission ratio increases. However, the friction coefficient at the start of meshing increases as the transmission ratio increases.
Figure 8c illustrates that the oil volume for cooling decreases as the transmission ratio increases. However, the oil volumes for cooling at the start of meshing and at the intersection of the single meshing teeth and the second meshing teeth areas initially decrease and then increase as the transmission ratio increases because two teeth are meshed simultaneously at these meshing points. One of the meshing teeth is found at start of meshing, and the friction coefficient and oil volume for cooling of this meshing tooth increase as the transmission ratio increases. However, the friction coefficient and oil volume for cooling of the other meshing tooth decrease as the transmission ratio increases. When the increase in the transmission ratio is small, the increase in the oil volume for cooling of the meshing tooth at the start of meshing is smaller than the decrease in the oil volume for cooling of the other meshing tooth. Thus, the total oil volume for cooling of the gear decreases. Moreover, when the transmission ratio increases considerably, the increase in the oil volume for cooling of the meshing tooth at the start of meshing is greater than the decrease in the oil volume for cooling of the other meshing tooth. Hence, the total oil volume for cooling of the gear increases.
5. Effect of Working Conditions
The contact and oil film states between the tooth surfaces under different working conditions, including the input speed and input torque, exhibit significant differences. This situation results in different oil volumes for lubrication and cooling. Therefore, the influence of working conditions on the oil volumes for lubrication and cooling during the meshing process should be examined.
An analysis method similar to that described in
Section 4 is used to investigate the influence of input speed on lubricating oil demand. The input speeds are set to 4000, 8000, 10,000 and 12,000 rpm. The parameters of the gear and input torque are identical to those described in
Section 3. The variations in the oil volume for lubrication during the meshing process at different input speeds are obtained through a series of simulations, as shown in
Figure 9a.
This figure indicates that the input speed of the gear has no effect on the meshing characteristics of the gear. The oil volume for lubrication increases as the input speed increases because the entrainment velocity and oil film thickness of each meshing tooth during the meshing process increases as the input speed increases. A series of simulations is performed to examine the influence of different input speeds on the friction coefficient during the meshing process, as shown in
Figure 9b, and partially enlarged images are obtained concurrently.
Oil film thickness and relative slide velocity increase as the input speed increases. The friction coefficient decreases as oil film thickness increases, but it increases as relative slide velocity increases.
Figure 9b illustrates that the friction coefficient decreases slightly as the input speed increases because of the considerable influence of oil film thickness on the friction coefficient. The oil volume for cooling at different input speeds during the meshing process is compared, as shown in
Figure 9c with a partially enlarged image.
This figure indicates that the oil volume for cooling increases as the input speed increases because the relative slide velocity exerts a significant influence on the oil volume for cooling although the friction coefficient decreases as input speed increases. The oil for cooling increases as the relative slide velocity increases.
To compare lubricating oil demand at different input torques during the meshing process, the input torques are set to 40, 60, 80 and 100 Nm. The parameters of the gear and input speed are the same as those described in
Section 3. The variations in the oil volume for lubrication during the meshing process at different input torques are obtained using a series of simulations, as shown in
Figure 10a.
This figure illustrates that the oil volume for lubrication decreases as the input torque increases because oil film thickness decreases as the input torque increases. The thin oil film leads to an increase in the friction coefficient as the input torque increases, as shown in
Figure 10b. A series of simulation models is used to examine the influence of different input torques on the oil volume for cooling during the meshing process, as shown in
Figure 10c.
The oil volume for cooling increases as the input torque increases due to the friction coefficient, as shown in
Figure 10c.