Next Article in Journal
Corrosion Resistance and Nano-Mechanical Properties of a Superhydrophobic Surface
Previous Article in Journal
Contact Load and Elastohydrodynamic Lubrication Analysis of Eccentric Bearings in RV Reducer Considering the Effects of Roller Profile Modification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Response and Lubrication Performance of Spur Gear Pair Under Time-Varying Rotation Speeds

1
School of Mechanical and Equipment Engineering, Hebei University of Engineering, Handan 056038, China
2
School of Mechanical and Aerospace Engineering, Nanyang Technological University, Singapore 639798, Singapore
3
State Key Laboratory of Reliability and Intelligence of Electrical Equipment, Hebei University of Technology, Tianjin 300401, China
4
College of Mechanical Engineering, Tianjin University of Technology, Tianjin 300401, China
*
Authors to whom correspondence should be addressed.
Lubricants 2025, 13(1), 15; https://doi.org/10.3390/lubricants13010015
Submission received: 4 December 2024 / Revised: 20 December 2024 / Accepted: 2 January 2025 / Published: 3 January 2025
(This article belongs to the Special Issue Surface Machining and Tribology)

Abstract

:
The rotation speed directly influences the vibration and lubrication behaviors of gear pairs, but studying the effects of time-varying rotation speeds during their operation poses substantial challenges. The present work proposed an approach to analyzing the dynamic response and lubrication performance of spur gear pairs under time-varying rotation speeds. A single-degree-of-freedom torsional dynamics model was established to capture the vibration responses and meshing forces of a gear pair, with the meshing stiffness modulated by the time-varying rotation speed. Additionally, a transient elastohydrodynamic lubrication model of the gear system was proposed to obtain the pressure pro-file and film shape, incorporating the effects of time-varying rotation speeds. Three types of time-varying rotation speeds were investigated: acceleration, deceleration, and oscillation. The results reveal that the time-varying rotation speed induces chaotic motion of the gear system, resulting in significant changes in the dynamic meshing force, entrainment velocity, and curvature radius of the gear pair compared to those in constant-speed scenarios. The lubrication performance under time-varying rotation speeds also shows diverse dynamic characteristics, highlighting significant differences from that observed under a constant rotation speed. These insights contribute to a more comprehensive understanding of gear dynamics under realistic operating conditions, enhancing gears’ performance and reliability in practical applications.

1. Introduction

The gear system has many distinguished merits, such as high transmission efficiency, high transmission accuracy, and high reliability. Consequently, gear transmissions are used widely in mechanical equipment. Investigating gears has attracted a large number of researchers, and research on gear vibration and lubrication has always been important. The gear pair’s vibration response is extremely complex, and it directly affects the smoothness and reliability of gear operation. The dynamic load caused by vibrations affects the formation of a lubricating oil film. Lubrication failure can increase the friction and wear and even result in tooth surface failure. Nevertheless, the dynamic response and lubrication behaviors of certain gear systems with time-varying rotation speeds, such as wind turbine and robot reducer gears, is still not known clearly.
A gear transmission system is a typical nonlinear dynamics system, and the earliest gear dynamics analyses can be traced back to the 1920s and early 1930s [1]. After nearly 100 years of development, several methods and models have been developed for studying the vibration responses of gear pairs. Kahraman et al. [2] established a single-degree-of-freedom model that accounted for the backlash and static transmission error, but the time-dependent characteristics of the meshing stiffness were neglected. Chaotic and subharmonic resonances were observed in their study. Blankenship and Kahraman [3] studied the complex behavior of a vibration system with clearance through analytical and numerical methods, and a gear dynamics test rig was designed to verify the solutions. Kahraman and Blankenship [4] observed the nonlinear behavior of a gear pair considering the clearance and a periodically time-varying mesh stiffness function through experiments. Cheng et al. [5] established a dynamics model of a gear system which could consider the lubrication. Their conclusion revealed that the surface roughness played an important role in the meshing stiffness, damping, and gear performance. Wu et al. [6] proposed a gear dynamics analysis approach that accounted for the effect of tooth cracks. Specifically, the influence of tooth cracks on the engaged stiffness of the gear pair was considered. Yi et al. [7] proposed a nonlinear dynamics model of a gear system that accounted for the time-varying pressure angle and the backlash caused by bearing deformation. They found that the system might approach a chaotic state early. Shahmasoorian et al. [8,9] undertook a tremendous amount of work on chaos estimation for and control over nonlinear systems. Shi et al. [10] developed a gear dynamics model that included the effect of the tooth separation and the backside mesh. Cheng et al. [11] studied the impact of teeth wear on the dynamic characteristics of gears. They found that the wear of the teeth changed the meshing stiffness and meshing position and further impacted the dynamic response of the gear pair. Cao et al. [12] investigated the effect of the load on the meshing stiffness and developed a nonlinear dynamics model that accounted for the force-dependent meshing stiffness, backlash, and profile deviation. Hasnijeh et al. [13] studied the vibration behaviors of a spur gear system considering stochastic external excitations. From the above discussion, the time-varying mesh stiffness, backlash, static transmission error, and external load have been the focuses of gear dynamics research. The effects of these external and parametric excitations on the vibration responses of gear pairs have been studied systematically. However, the studies discussed above neglected the variations in the rotation speed, which play a crucial role in the meshing stiffness. The meshing stiffness is associated with the engagement position of the gear joint, and it is a periodic function of time when the rotation speed of the gears is constant. For this condition, the time for the gear pair to rotate by one tooth is fixed, and it is called the meshing cycle. However, in practice, the rotation speed of the gear does not always remain constant, such as in industrial robots and machine tools; the speed of a gear pair varies over time. The meshing stiffness is no longer periodic, with a constant period as a result. The meshing stiffness has to be converted from the space domain into the time domain based on the rotation speed. Khabou et al. [14] studied the dynamic behavior of gear sets powered by electric and four-cylinder diesel engines. The variations in the rotation speed and torque were considered in their work. However, the influence of a time-varying rotation speed was not investigated.
Elastohydrodynamic lubrication (EHL) is another significant field of gear research. The analysis approach to investigations of gear EHL has developed maturely. The current EHL model, which can account for factors including the temperature effect, surface roughness, non-Newtonian behavior of the lubricant, and transient effects, meets the engineering practice requirements. Li and Kahraman [15,16,17] studied the effect of transient and dynamic behaviors on the lubrication performance of gear transmission. The surface roughness was considered in their works. As with surface roughness, the manufacturing errors is significant to the lubrication performance of gear systems. Hjelm and Wahlström [18] established a thermal elastohydrodynamic lubrication model which could take into account the manufacturing errors. The results revealed that manufacturing errors had a great influence on the lubrication behavior of the gear system. Xiao et al. [19] studied normal and tangential oil film damping of gears using a transient thermal EHL model. The non-Newtonian effect of the lubricant was also considered in their work. Jamali et al. [20] proposed a transient EHL model for point contact to consider the longitudinal crowning of the spur gear teeth. They found that the film thickness close to the tooth edges obviously increased. Liu and Gong [21] established a thermal EHL model for helical gears to investigate the influences of profile modifications, the input torque, and the rotation speed on the flash temperature. To study the lubrication characteristics of spur gears under transient impacts and vibration conditions, Huang et al. [22] proposed an EHL model that accounted for thermal and squeeze effects. The results revealed that the film temperature and pressure increased remarkably while the film thickness decreased. Yin et al. [23] studied the lubrication performance of double involute gears at the graded position of the tooth waist. They found the thermal effect had a negligible influence on the lubrication characteristics at this position. In summary, gear EHL has already been studied thoroughly. As is well known, the rotation speed, which affects the curvature radius, entrainment velocity, and contact force, is essential to the lubrication performance of gear pairs. The rotation speed of the gear pair has usually been assumed to be constant in previous studies. As discussed earlier, the rotation speed is not always constant. However, an EHL model and an analysis method for gear pairs with time-varying rotation speeds are still absent, and thus the lubrication characteristics of gear systems with time-varying rotation speeds remain unclear.
This work proposed an approach to studying the transient response and lubrication performance of spur gear pairs under time-varying rotation speeds. First, a nonlinear dynamics model of the gear system was established, incorporating the time-varying meshing stiffness, static transmission error, and backlash. The meshing stiffness model, based on the Ishikawa method, was developed by considering the variation in the rotation speed. Additionally, a transient EHL model for the gear pair was formulated to account for the dynamic changes in lubrication caused by variations in speed. Three types of rotation speeds of acceleration, deceleration, and oscillation were analyzed to assess their impacts on the gear dynamics and lubrication behavior. This study provides insights into the influence of time-varying rotation speeds on both the transient responses and lubrication characteristics of a spur gear pair, advancing our understanding of the gear performance under realistic operating conditions.

2. Materials and Methods

2.1. Time-Varying Stiffness

Several approaches are usually employed to obtain the meshing stiffness along the meshing line, and the Ishikawa approach was used in this work. As can be seen in Figure 1, the gear tooth is simplified as a combination of a rectangle and a trapezoid according to the Ishikawa method. The calculation process is shown in detail in Refs. [24,25]. The meshing stiffness under different rotational angles was ultimately obtained as follows:
k ( θ ) = F n δ
where k denotes the meshing stiffness, which is a function of the rotational angle. An involute spur gear pair is employed in this paper, and its parameters are shown in Table 1; the meshing stiffness is obtained and is shown in Figure 2.
For simplicity, the meshing stiffness is translated into a square wave. The maximum stiffness is the mean of the double-tooth meshing region, and the minimum is the mean of the single-tooth area. The meshing stiffness and the simplified stiffness are depicted in Figure 2. The maximum and minimum values of the square wave stiffness are defined as follows:
k max = 1 θ d 0 θ d k ( θ ) d θ
k min = 1 θ s θ d θ d + θ s k ( θ ) d θ
where kmax and kmin represent the maximum and minimum square wave stiffness; θd and θs are the angular displacement in double-tooth and single-tooth meshing, and they can be defined as follow:
θ s = θ t ( 2 ε )
θ d = θ t ( ε 1 )
where θt is the meshing angle of one tooth, which defines θt = 2π/z, and ε is the contact ratio.
Three types of time-varying rotations are employed in this paper. The first and second types of rotation speed are acceleration and deceleration, respectively.
θ ˙ = θ ¨ t
where θ ˙ is the rotation speed, and θ ¨ denotes the acceleration of the rotation speed; if θ ¨ is negative, this corresponds to deceleration.
The third type is a harmonic function of time, which can be expressed as follows:
θ ˙ = θ ˙ c + A sin ( ω r t + φ r )       ω r = 2 π n 1 / 60
where θ ˙ c denotes the constant component of the rotation speed, A represents the amplitude of varying components, ωr is the angular frequency, and φr denotes the phase angle. If the initial rotation angle is zero, the rotation angle can be calculated as follows, respectively:
θ = θ ˙ t + θ ¨ t 2 2
θ = θ ˙ c t + A ω r 1 cos ( ω r t + φ r )
Based on Equations (5) and (6), the engaged stiffness is a time function.
k ( t ) = k max , i + 1 ε θ t < θ < i θ t k min , e l s e
where i is the cycle number. The gear meshing stiffness model can consider time-varying rotation speeds. The meshing stiffness can easily transform from the space domain to the time domain regardless of the form of the time-varying speed.

2.2. The Dynamics Model of the Gear Pair

A gear pair is a strong nonlinear dynamics system, and several typical nonlinear phenomena, such as chaos and jump discontinuities, exist in gear dynamic systems. These nonlinear phenomena mainly originate from the time-variant mesh stiffness and backlash. The rotation speed plays an important role in the time-varying meshing stiffness. A torsional dynamic model is developed to investigate the influence of the time-dependent speed on the dynamic behavior of the gear system. As illustrated in Figure 3, the pinion and gear are considered as two rigid disks with corresponding masses. The contact of the tooth is replaced by a massless spring and massless damping. The kinetics model can include the backlash and static transmission error. According to Newton’s law, a dynamics equations can be obtained [26]:
I 1 θ ¨ 1 + R b 1 c R b 1 θ ˙ 1 R b 2 θ ˙ 2 e ˙ + R b 1 k ( t ) R b 1 θ 1 R b 2 θ 2 e = T 1 I 2 θ ¨ 2 + R b 2 c R b 1 θ ˙ 1 R b 2 θ ˙ 2 e ˙ R b 2 k ( t ) R b 1 θ 1 R b 2 θ 2 e = T 2
where Ii is the rotary inertia, θi is the torsional angular displacement, Rbi denotes the base cycle radius of the gear, and i = 1, 2 are the pinion and the gear, respectively. T1 and T2 are the input and output torque, and e denotes the transmission error due to manufacturing. k(t) represents the mesh stiffness, which is time-varying. c denotes constant viscous damping:
c = 2 ξ m e k m
where ξ represents the damping ratio, me denotes the equivalent mass of the gear system, and km is the average mesh stiffness.
Let x = Rb1θ1Rb2θ2e; Equation (11) can be rewritten as follows [26]:
m e x ¨ + c x ˙ + k ( t ) g ( x ) = f m e d 2 e d t 2
in which
m e = I 1 I 2 I 1 R b 2 2 + I 2 R b 1 2
g ( x ) = x b , x > b 0 , b x b x + b , x < b
f = T 1 R b 1 = T 2 R b 2
where g(x) is the nonlinear displacement function. b is the backlash. The static transmission error is assumed to be sinusoidal [26]:
e = e r sin ( ω e t + φ e )
where er represents the amplitude; ωe denotes the angular frequency; φe represents the phase angle.
Based on the dimensionless variables and parameters below, X = x/b,   ω n = 2 ξ k m / m e , K = k/km, F = f/(bkm), T = n, Ωe = ωe/ωn, F = er/b. The dimensionless form of Equation (14) can be written as follows:
X ¨ + 2 ξ X ˙ + K G ( X ) = F + F a Ω e 2 sin ( Ω e T + φ e )
G ( X ) = X 1 , X > 1 0 ,         1 X 1 X + 1 , X < 1
The time-varying rotation speed mainly affects the meshing stiffness. The meshing stiffness is no longer a function of time with a fixed period. The advantage of this model is that it can be used to tackle dynamic response problems for a gear pair with any form of time-varying rotation speed.
The meshing force on the tooth can be calculated according to the dynamic response and the meshing stiffness:
f t = k i g ( x ) + c x ˙

2.3. Geometry and Kinematics Analysis

According to EHL theory, the contact of the gear teeth can be treated as the contact between a plane and a cylinder. The radius of the cylinder defines the equivalent radius which can be obtained using the geometry analysis. The curvature radii of the two gear teeth at the contact point can be expressed as follows:
R 1 = R b 1 tan ( θ 0 + θ m i )
R 2 = s 2 ( R 1 s 1 )
where θ0, s1, s2 are the geometry parameters displayed in Figure 4. θmi is the rotation angle of the ith tooth, which changes from zero to θi. θi is the rotation angle of the ith tooth from the approach point to the recession point. The curvature radii can be transformed into the time domain with the known gear rotational speed. With the curvature radii of the contact teeth, the equivalent radius can be obtained [18]:
R = R 1 R 2 R 1 + R 2
The transient tangential velocity of the teeth at the contact point can be expressed as follows [27]:
u 1 = θ ˙ 1 R 1
u 2 = θ ˙ 2 R 2
where θ ˙ 1 and θ ˙ 2 are the rotation speed of the gear and the pinion, respectively. The entrainment velocity can be written as follows:
u = u 1 + u 2 2
where u1 and u2 are the rolling speeds of the teeth at the contact point, respectively.

2.4. The Transient EHL Model

A lubricant film forms in the contact region of the teeth. The pressure of the lubricant film is controlled by the Reynolds equation. When the thermal effect is neglected, the line contact Reynolds equation for Newtonian fluid reads as follows [27]:
x ρ h 3 η p x = 12 u x ρ h + t ρ h
where x represents the coordinate along the rolling direction; p donates the hydrodynamic pressure; h represents the film’s thickness; ρ and η represent the density and viscosity of the lubricant, respectively; and t represents the time.
The boundary and cavitation conditions should be known when solving the Reynolds equation, which define the following:
x = x 0 , p = 0 x = x e , p = p x x = x e = 0
where x0 and xe are, respectively, the inlet and outlet locations of the lubricant.
The film thickness equation reads as follows [27]:
h ( x ) = h 0 + x 2 2 R 4 π E x 0 x e p ( s ) ln x s d s
where h0 donates the rigid body displacement, and E denotes the composite Young’s modulus, which can be expressed as follows [18]:
1 E = 1 2 1 υ 1 2 E 1 + 1 υ 2 2 E 2
The equilibrium equation reads as follows:
w = x 0 x e p ( x ) d x
The Dowson–Higginson equation and the Roelands equation are employed to describe the relation of density and viscosity to pressure, respectively [28].
ρ = ρ 0 ( 1 + 0.6 p 1 + 1.7 p )
η = η 0 exp ( ln η 0 + 9.67 ) [ 1 + ( 1 + 5 . 1 × 10 - 9 p ) z
where ρ0 and η0 represent the density and viscosity of the lubricant in ambient conditions; Z denotes the viscosity–pressure index and is set to 0.68 in this paper [29,30].
A lack of dimensions can reduce the number of parameters and simplify the control equations. Hence, the following parameters are introduced to obtain a dimensionless form of the control equations.
X = x b 0 , H = h R 0 b 0 2 , P = p P H , η ¯ = η η 0 , ρ ¯ = ρ ρ 0
T = t u s b 0 , P H = E b 0 4 R 0 , b 0 = 8 w 0 R 0 π E , U = u u s , ε = ρ ¯ H 3 η ¯ λ , λ = 12 η 0 u s R 0 2 P H b 0 3
where b0 represents the half-width of the Hertz contact region at the approach point of the teeth; R0 is the equivalent radius at the approach point; PH donates the maximum Hertz pressure of the approach point of contact of the teeth; us is the rolling speed of the approach point of contact of the teeth when the gear rotation speed corresponds to deceleration or sinusoidal variation, and the rolling speed of the recession point when the gear rotation speed corresponds to acceleration; w0 represents the force on a unit width of the gear tooth at the approach point.
The dimensionless forms of the Reynolds equation, the film thickness equation, and the equilibrium equation are defined as follows:
X ε P X = U ρ ¯ H X + ρ ¯ H T
H = H 00 + X 2 R 0 2 R R 0 π R X 0 X e P ( S ) ln X S d S

2.5. The Numerical Procedure

Varying rotation speeds mainly impact the mesh stiffness of the gear system. Obtaining the mesh stiffness is the first step of a dynamics and lubrication analysis. The meshing stiffness along the line of action is calculated using the Ishikawa approach. The meshing stiffness is transformed into the time domain according to the time-varying rotation speed. A dynamics model of the gear system considering time-varying rotation speeds is proposed. The fourth-order Runge–Kutta method is employed for the solution of the dynamics equations. Then, the meshing force can be obtained according to the transient response of the gear pair. The curvature radii and the entrainment velocity are obtained through a geometry and kinematics analysis of the gear system. The transient EHL model is established to obtain the film thickness profile and the pressure distribution. The meshing time of each tooth is not equal. The number of discrete points is also unequal in each meshing cycle period. The inlet and outlet positions of the lubricant are x0 = −4 and xe = 1.5, and 512 nodes exist in the X direction. The pressure and the film thickness from the previous time step are employed to initialize the current time step. The pressure and film thickness at the approach point are the solutions for the steady state. The pressure is Hertz pressure when the speed of the gear is equal to zero. The Gauss–Seidel iterative method is employed to obtain the pressure in light-load regions, and the Jacobi dipole iterative method in is used high-load regions. The iteration ends and proceeds to the next time step when the convergence criteria for the pressure and load are both satisfied. The convergence criteria are defined as follows:
p i n + 1 p i n p i n + 1 ε p
x 0 x e p ( x ) d x π / 2 π / 2 ε w
where εp and εw are the thresholds for the relative error in the pressure and load; we take εp = 10−5 and εw = 10−3 in this paper. The equation was solved in this article using MATLAB R2021a.

3. Results and Discussions

Three types of time-varying rotation speeds are investigated in this paper. The first type of speed accelerates from zero to 200 rad/s at constant acceleration. The second decelerates from 200 rad/s to zero at constant deceleration. The third is the sum of a constant speed and a sinusoidal wave. The amplitude of the sinusoidal wave is 0.5 times of the constant speed, and the frequency is the same as the pinion’s rotation.
The parameters of the gear pair and the lubricant are summarized in Table 1. The types of time-varying rotation speeds and corresponding time-varying meshing stiffnesses are shown in Figure 5, Figure 6 and Figure 7. The meshing stiffness is obtained using the model proposed in Section 2.1. The gears’ rotation speed has a significant effect on the meshing stiffness. No periodicity exists in the meshing stiffness when the rotation speed corresponds to acceleration or deceleration. As depicted in Figure 5b, the meshing stiffness is simplified as a square wave; the lower edge of the square wave is a single-tooth mesh region, and the higher edge is a double-tooth mesh region. The tooth mesh cycle decreases as the rotation speed increases. Figure 7b demonstrates the meshing stiffness when the gear’s rotation speed is sinusoidal. The meshing stiffness shows regular changes in density. The red lines in Figure 5, Figure 6 and Figure 7 denote the meshing cycles selected for studying the following lubrication performance of the gear pair.

3.1. Dynamic Response of the Gear System with Time-Varying Rotation Speeds

The transient displacement of the gear pair with different types of rotation speed is shown in Figure 8. The dynamic response is obtained by solving the dynamics model for the gear pair. As can be seen, the response is more complex than the response of a constant-rotation-speed gear pair. No evident period is observed in the dynamic response. Severe fluctuation appears when the teeth enter into meshing and then gradually disappears. Fluctuation also appears when the number of teeth engaged changes. These fluctuations are the transient response of the gear pair, and the frequency is close to the natural frequency of the system. As shown in Figure 8a,b, the transient response is followed by a steady response. The fluctuation in the steady response regions is caused by the static transmission error, which is an internal excitation with a high frequency. The backlash is taken as 10−4 m in this paper, so the teeth are always stay in contact, and no impact occurs. The other fluctuation is caused by time-varying meshing stiffness, and its frequency is not constant. The frequency increases with the rotation speed.
The Fourier spectra of the dynamic transmission error shown in Figure 9a–c correspond to acceleration, deceleration, and a sinusoidal rotation speed, respectively. As can be seen, rich frequency content exists in the amplitude-frequency spectrum The conclusion can be drawn that the motions of the three rotational speeds are all chaotic. The amplitude–frequency spectra for acceleration motion and deceleration motion are essentially identical. The amplitude–frequency spectrum for a sinusoidal rotation speed is rougher. The gear system is more likely to enter into chaotic motion when the rotation speed is time-varying.

3.2. Effect of Time-Varying Rotation Speeds on Gear Lubrication

Compared to the transmitted torque, the rotation speed plays a more significant role in the formation of the film of gear lubrication. The transmitted torque affects the dynamic meshing force, while the rotation speed also affects the curvature radius and the entrainment velocity. The curvature radius and the entrainment velocity in the third, fourth, and fifth stages of the acceleration process are illustrated in Figure 10. The curvature radius, entrainment velocity, and meshing force in each tooth meshing cycle are different when the rotation speed is time-varying. Figure 10a shows the curvature radii in different tooth meshing cycles. The values for the curvature radius from the approach point to the recession point are equivalent. However, the meshing time for each tooth is different. For the acceleration process, the meshing time is decreased. A trend can also be seen in Figure 10b for the entrainment velocity in different tooth meshing cycles. Unlike the curvature radius, the value of the entrainment velocity is different in each tooth meshing cycle. During the acceleration process, the entrainment velocity increases with time. This is because the entrainment velocity is positively correlated with the rotation speed. Figure 11 depicts the engagement force of the third to the fifth tooth. Fluctuation is also observed when the number of teeth engaged changes. The transient response disappears, and a steady-state response is entered. The frequency of the steady-state response is the same as the static transmission error. The engagement force in the single-tooth mesh region is larger than that in the double-tooth mesh regions. The contact area of the meshing teeth is related to the meshing force. So, the contact area for each tooth is different due to the meshing force on each tooth varying when the rotation speed is time-varying. The curvature radius, entrainment velocity, and engagement force on the tooth in the deceleration process and with a sinusoidal rotation speed also demonstrate different characteristics compared with those at a constant rotation speed. The Dowson and Higginson’s (D-H) formula for the minimum film thickness is employed to validate the proposed model. The 15th tooth’s meshing cycles at a sinusoidal rotation speed are selected. As shown in Figure 12, the minimum film thickness obtained using the present model is slightly smaller than that with the D-H formula. This is because the D-H formula neglects the transient effects.
Considering the complex calculation of the EHL simulation, six tooth meshing cycles of each type of time-varying speed are selected. For the acceleration and deceleration processes, the 3rd, 13th, 23rd, 33rd, 43rd, and 43rd tooth meshing cycles are chosen. For a sinusoidal rotation speed, the 3rd, 9th, 15th, 21st, 27th, and 33rd tooth meshing cycles are selected. The selected meshing cycles are highlighted in red in Figure 5, Figure 6 and Figure 7. The dimensionless pressure and the film thickness distribution in the third engagement cycle of the acceleration process are depicted in Figure 13. The typical characters of EHL line contact, such as the quadratic pressure spike and the neck of the film thickness, are displayed at each time point. The pressure profile and the film distribution exhibit fluctuations caused by the meshing force along the time axis. The trend is the same for the meshing force. The pressure in the engaged single-tooth region is larger, and the film thickness is thinner than that in the double-tooth meshing regions.
The dimensionless minimum and central film thicknesses under three types of rotation speed are depicted in Figure 14, Figure 15 and Figure 16. Figure 14 demonstrates the minimum and central film thicknesses in different meshing cycles of the acceleration process. The minimum and central film thickness both increase from mesh-in to mesh-out. The film thickness in different meshing cycles increases with an increasing rotation speed. Hence, the wear of the gear pair is more serious at the starting stage. Fluctuations exist in both the minimum and central film thicknesses, and the frequency is identical to the meshing force. As discussed previously, the meshing force on the tooth has two dominant frequencies: the higher frequency is the natural frequency of the system, and the lower frequency is caused by the static transmission error. The two frequencies are constant, and the tooth meshing cycle in the acceleration process decreases. So, the number of fluctuations decreases as the rotation speed increases. The range of fluctuations in the minimum film thickness is greater than that for the central film thickness. The high-frequency component is less obvious in the central film thickness. The film thickness in the engaged single-tooth region is thinner than that in the double-tooth meshing regions. The minimum film thickness suddenly increases when single-tooth meshing changes to double-tooth meshing. The minimum and central film thicknesses in different meshing cycles of the deceleration process are depicted in Figure 15. The film thickness decreases with a decreasing rotation speed. The film thicknesses in different meshing cycles increase from mesh-in to mesh-out. Nevertheless, they decrease as the rotation speed decreases.
Figure 16 depicts the minimum and central film thicknesses in different meshing cycles when the rotation speed changes in sin form. The film thickness in different meshing cycles changes periodically with the rotation speed. The range of variation in the film thickness increases with rotation speed.

4. Conclusions

The dynamic response and lubrication performance of gear pairs at time-varying rotation speeds were investigated in the present work. Three time-varying rotation speed types were selected: acceleration, deceleration, and oscillation. The meshing stiffness, dynamic response, and lubrication characteristics of an involute spur gear pair were obtained. Some of the findings of this study are shown as follows:
  • The time-varying rotation speeds affect the meshing stiffness due to the variation in the rotation time per tooth. The type of time-varying rotation plays a significant role in the shape of the meshing stiffness of the gear pair.
  • The time-varying rotation speeds induce chaotic motion of the gear pair. The vibration is higher when the frequency of the rotation speed is lower, and the amplitude of the rotation speed primarily influences the fluctuation range of the dynamic transmission error.
  • The time-varying rotation speed influences the lubrication characteristics through changes in the curvature radius, entrainment velocity, and meshing force. The film thickness of each tooth is different, which makes it difficult to predict the lubricant and wear performance of the gear pair.
This study highlights the complex relationship between speed variations and dynamic behavior and lubrication performance, offering valuable insights for optimizing gear performance in variable-speed applications. The rotation speed of a given gear system not only exhibits variability but also randomness, such as for the gear pairs used in wind turbines. The models and method proposed in the paper will support dynamics and lubrication studies in the future.

Author Contributions

Conceptualization: J.P. and Y.T. (Yourui Tao). Methodology, J.P. Software: J.P. and H.H. Validation: Y.T. (Yuanyuan Tian) and M.W. Formal analysis: M.W. Investigation: Z.G. Resources: Y.T. (Yuanyuan Tian). Data curation: Y.T. (Yourui Tao). Writing—original draft preparation: J.P. Writing—review and editing: H.H. Visualization: Y.T. (Yuanyuan Tian). Supervision: Y.T. (Yourui Tao). Project administration: Y.T. (Yuanyuan Tian). Funding acquisition: Y.T. (Yourui Tao). All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (U23A6017, 52075146).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Özgüven, H.N.; Houser, D.R. Dynamic Analysis of High Speed Gears by Using Loaded Static Transmission Error. J. Sound Vib. 1988, 125, 71–83. [Google Scholar] [CrossRef]
  2. Kahraman, A.; Singh, R. Non-linear Dynamics of a Spur Gear Pair. J. Sound Vib. 1990, 142, 49–75. [Google Scholar] [CrossRef]
  3. Kahraman, A.; Blankenship, W. Interactions between commensurate parametric and forcing excit-tions in a system with clearance. J. Sound Vib. 1996, 194, 317–336. [Google Scholar] [CrossRef]
  4. Kahraman, A.; Blankeship, G.W. Experiments on nonlinear dynamic behavior of an oscillator with clearance and periodically time-varying parameters. J. Appl. Mech. 1997, 64, 217–226. [Google Scholar] [CrossRef]
  5. Cheng, G.; Ma, J.; Li, J.; Sun, K.; Wang, K.; Wang, Y. Study on the dynamic characteristics of gears considering surface topography in a mixed lubrication state. Lubricants 2024, 12, 7. [Google Scholar] [CrossRef]
  6. Wu, S.; Zuo, M.; Parey, A. Simulation of spur gear dynamics and estimation of fault growth. J. Sound Vib. 2008, 317, 608–624. [Google Scholar] [CrossRef]
  7. Yi, Y.; Huang, K.; Xiong, Y.; Sang, M. Nonlinear dynamic modelling and analysis for a spur gear system with time-varying pressure angle and gear backlash. Mech. Syst. Signal Process. 2019, 132, 18–34. [Google Scholar] [CrossRef]
  8. Gharagozloo, M.; Shahmansoorian, A. Chaos control in gear transmission system using GPC and SMC controllers. J. Appl. Comput. Mech. 2022, 8, 545–556. [Google Scholar]
  9. Shams, Z.; Shahmansoorian, A. Fault estimation based on observer for chaotic Lorenz system with bifurcation problem. Trans. Inst. Meas. Control 2020, 42, 576–585. [Google Scholar] [CrossRef]
  10. Shi, J.; Goua, X.; Zhu, L. Modeling and analysis of a spur gear pair considering multi-state mesh with time-varying parameters and backlash. Mech. Mach. Theory 2019, 134, 582–603. [Google Scholar] [CrossRef]
  11. Geng, Z.; Chen, M.; Wang, J.; Xia, Y.; Kong, Y.; Xiao, K. Nonlinear Dynamic analysis of a spur gear pair system with wear considering the meshing position. Lubricants 2024, 12, 25. [Google Scholar] [CrossRef]
  12. Cao, Z.; Chen, Z.; Jiang, H. Nonlinear dynamics of a spur gear pair with force-dependent mesh stiffness. Nonlinear Dyn. 2020, 99, 1227–1241. [Google Scholar] [CrossRef]
  13. Hasnijeh, S.G.; Poursina, M.; Leira, B.J.; Karimpour, H.; Chai, W. Stochastic dynamics of a nonlinear time-varying spur gear model using an adaptive time-stepping path integration method. J. Sound Vib. 2019, 447, 170–185. [Google Scholar] [CrossRef]
  14. Khabou, M.T.; Bouchaala, N.; Chaari, F.; Fakhfakh, T.; Haddar, M. Study of a spur gear dynamic behavior in transient regime. Mech. Syst. Signal Process. 2011, 125, 3089–3101. [Google Scholar] [CrossRef]
  15. Li, S.; Kahraman, A. A transient mixed elastohydrodynamic lubrication model for spur gear pairs. ASME J. Tribol. 2010, 132, 011501. [Google Scholar] [CrossRef]
  16. Li, S.; Kahraman, A. Influence of dynamic behaviour on elastohydrodynamic lubrication of spur gears. Proc. I. Mech. E. Part J J. Eng. Tribol. 2011, 225, 740–753. [Google Scholar] [CrossRef]
  17. Li, S.; Kahraman, A. prediction of spur gear mechanical power losses using a transient elastohydrodynamic lubrication model. Tribol. Trans. 2010, 53, 554–563. [Google Scholar] [CrossRef]
  18. Hjelm, R.; Wahlström, J. Influence of manufacturing error tolerances on thermal EHL behavior of gears. Lubricants 2022, 10, 323. [Google Scholar] [CrossRef]
  19. Xiao, Z.; Li, Z.; Shi, X.; Zhou, C. Oil film damping analysis in non-Newtonian transient thermal elastohydrodynamic lubrication for gear transmission. J. Appl. Mech.-Trans. ASME 2018, 85, 035001. [Google Scholar] [CrossRef]
  20. Jamali, H.U.; Aljibori, H.S.S.; Jweeg, M.J.; Abdullah, O.I.; Ruggiero, A. An investigation on the teeth crowning effects on the transient EHL performance of large-scale wind turbine spur gears. Lubricants 2023, 11, 462. [Google Scholar] [CrossRef]
  21. Liu, Y.; Gong, J. Study on TEHL flash temperature of helical gear pair considering profile modification. Ind. Lubr. Tribol. 2023, 75, 17–26. [Google Scholar] [CrossRef]
  22. Huang, X.B.; Yang, B.T.; Wang, Y.Q. Influences of transient impact and vibration on the lubrication performance of spur gears. Proc. I Mech. E. Part J J. Eng. Tribol. 2021, 235, 274–289. [Google Scholar] [CrossRef]
  23. Yin, Z.; Fan, Z.; Wang, M. Thermal elastohydrodynamic lubrication characteristics of double involute gears at the graded position of tooth waist. Tribol. Int. 2020, 144, 106028. [Google Scholar] [CrossRef]
  24. Shi, J.; Ma, X.; Xu, C.; Zang, S. Meshing stiffness analysis of gear using the Ishikawa method. Appl. Mech. Mater. 2013, 401–403, 203–206. [Google Scholar] [CrossRef]
  25. Wang, J.; He, G.; Zhang, J.; Zhao, Y.; Yao, Y. Nonlinear dynamics analysis of the spur gear system for railway locomotive. Mech. Syst. Signal Process. 2017, 85, 41–55. [Google Scholar] [CrossRef]
  26. Pei, J.; Han, X.; Tao, Y.; Feng, S. Lubrication reliability analysis of spur gear systems based on random dynamics. Tribol. Int. 2021, 153, 106606. [Google Scholar] [CrossRef]
  27. Larsson, R. Transient non-Newtonian elastohydrodynamic lubrication analysis of an involute spur gear. Wear 1997, 207, 67–73. [Google Scholar] [CrossRef]
  28. Masjedi, M.; Khonsari, M.M. Film thickness and asperity load formulas for line-contact elastohydrodynamic lubrication with provision for surface roughness. J. Tribol. 2012, 134, 011503. [Google Scholar] [CrossRef]
  29. Gu, C.; Meng, X.; Wang, S.; Ding, X. Research on mixed lubrication problems of the non-gaussian rough textured surface with the influence of stochastic roughness in consideration. J. Tribol. 2019, 141, 121501. [Google Scholar] [CrossRef]
  30. Pei, J.; Han, X.; Tao, Y.; Feng, S. Mixed elastohydrodynamic lubrication analysis of line contact with non-gaussian surface roughness. Tribol. Int. 2020, 151, 106449. [Google Scholar] [CrossRef]
Figure 1. Approximate tooth profile using the Ishikawa method.
Figure 1. Approximate tooth profile using the Ishikawa method.
Lubricants 13 00015 g001
Figure 2. Meshing stiffness.
Figure 2. Meshing stiffness.
Lubricants 13 00015 g002
Figure 3. Dynamics model of gear system.
Figure 3. Dynamics model of gear system.
Lubricants 13 00015 g003
Figure 4. Geometric size of gear system.
Figure 4. Geometric size of gear system.
Lubricants 13 00015 g004
Figure 5. Acceleration and corresponding meshing stiffness: (a) rotational speed, (b) meshing stiffness.
Figure 5. Acceleration and corresponding meshing stiffness: (a) rotational speed, (b) meshing stiffness.
Lubricants 13 00015 g005
Figure 6. Deceleration and corresponding meshing stiffness: (a) rotational speed, (b) meshing stiffness.
Figure 6. Deceleration and corresponding meshing stiffness: (a) rotational speed, (b) meshing stiffness.
Lubricants 13 00015 g006
Figure 7. Sinusoidal rotational speed and corresponding meshing stiffness: (a) rotational speed, (b) meshing stiffness.
Figure 7. Sinusoidal rotational speed and corresponding meshing stiffness: (a) rotational speed, (b) meshing stiffness.
Lubricants 13 00015 g007
Figure 8. Dynamic response of gear pair with time-varying rotation speeds: (a) acceleration, (b) deceleration, and (c) sin form.
Figure 8. Dynamic response of gear pair with time-varying rotation speeds: (a) acceleration, (b) deceleration, and (c) sin form.
Lubricants 13 00015 g008
Figure 9. Fourier spectra of gear pair with time-varying rotation speeds: (a) acceleration, (b) deceleration, and (c) sin form.
Figure 9. Fourier spectra of gear pair with time-varying rotation speeds: (a) acceleration, (b) deceleration, and (c) sin form.
Lubricants 13 00015 g009
Figure 10. Curvature radius and entrainment velocity in different tooth meshing cycles in the acceleration process: (a) curvature radius and (b) entrainment velocity.
Figure 10. Curvature radius and entrainment velocity in different tooth meshing cycles in the acceleration process: (a) curvature radius and (b) entrainment velocity.
Lubricants 13 00015 g010
Figure 11. Meshing force in different meshing cycles of the acceleration process.
Figure 11. Meshing force in different meshing cycles of the acceleration process.
Lubricants 13 00015 g011
Figure 12. The minimum film thickness obtained using the D-H formula and the present model.
Figure 12. The minimum film thickness obtained using the D-H formula and the present model.
Lubricants 13 00015 g012
Figure 13. The pressure and film of the 3rd tooth from mesh-in to mesh-out in the acceleration process: (a) pressure profile and (b) film distribution.
Figure 13. The pressure and film of the 3rd tooth from mesh-in to mesh-out in the acceleration process: (a) pressure profile and (b) film distribution.
Lubricants 13 00015 g013
Figure 14. The minimum and central film thicknesses in different meshing cycles of the acceleration process (af) six different meshing cycles.
Figure 14. The minimum and central film thicknesses in different meshing cycles of the acceleration process (af) six different meshing cycles.
Lubricants 13 00015 g014
Figure 15. The minimum and central film thicknesses in different meshing cycles of the deceleration process (af) six different meshing cycles.
Figure 15. The minimum and central film thicknesses in different meshing cycles of the deceleration process (af) six different meshing cycles.
Lubricants 13 00015 g015
Figure 16. The minimum and central film thicknesses in different meshing cycles of the sin form: (af) six different meshing cycles.
Figure 16. The minimum and central film thicknesses in different meshing cycles of the sin form: (af) six different meshing cycles.
Lubricants 13 00015 g016
Table 1. Gear data.
Table 1. Gear data.
Symbol, UnitGear Pair
Number of teeth, z1:z235:58
Pressure angle, °20
Module, mm3
Face width, mm20
Rotational inertia, kg∙mm21860/13,590
Constant component of rotation speed, rad/min2000
Center distance, mm139.5
Modulus of elasticity, GPa206
Damping ratio0.05
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pei, J.; Tian, Y.; Hou, H.; Tao, Y.; Wu, M.; Guan, Z. Dynamic Response and Lubrication Performance of Spur Gear Pair Under Time-Varying Rotation Speeds. Lubricants 2025, 13, 15. https://doi.org/10.3390/lubricants13010015

AMA Style

Pei J, Tian Y, Hou H, Tao Y, Wu M, Guan Z. Dynamic Response and Lubrication Performance of Spur Gear Pair Under Time-Varying Rotation Speeds. Lubricants. 2025; 13(1):15. https://doi.org/10.3390/lubricants13010015

Chicago/Turabian Style

Pei, Jiaxing, Yuanyuan Tian, Hongjuan Hou, Yourui Tao, Miaojie Wu, and Zhigang Guan. 2025. "Dynamic Response and Lubrication Performance of Spur Gear Pair Under Time-Varying Rotation Speeds" Lubricants 13, no. 1: 15. https://doi.org/10.3390/lubricants13010015

APA Style

Pei, J., Tian, Y., Hou, H., Tao, Y., Wu, M., & Guan, Z. (2025). Dynamic Response and Lubrication Performance of Spur Gear Pair Under Time-Varying Rotation Speeds. Lubricants, 13(1), 15. https://doi.org/10.3390/lubricants13010015

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop