Research Article  Open Access
Dynamic Model of Spur Gear Pair with Modulation Internal Excitation
Abstract
In the actual measurements, vibration and noise spectrum of gear pair often exhibits sidebands around the gear mesh harmonic orders. In this study, a nonlinear timevarying dynamic model of spur gear pair was established to predict the modulation sidebands caused by the AMFM modulation internal excitation. Here, backlash, modulation timevarying mesh stiffness, and modulation transmission error are considered. Then the undamped natural mode was studied. Numerical simulation was made to reveal the dynamic characteristic of a spur gear under modulation condition. The internal excitation was shown to exhibit obvious modulation sideband because of the modulation timevarying mesh stiffness and modulation transmission error. The RungeKutta method was used to solve the equations for analyzing the dynamic characteristics with the effect of modulation internal excitation. The result revealed that the response under modulation excitation exhibited obvious modulation sideband. The response under nonmodulation condition was also calculated for comparison. In addition, an experiment was done to verify the prediction of the modulation sidebands. The calculated result was consistent with the experimental result.
1. Introduction
The first systematic efforts to analyze the gear dynamics were in the 1920s and early 1930s [1]. Many researchers have made a great contribution to the mathematical models and dynamic behaviors of gear transmission systems. Lumpedparameter models are widely used in modeling the gear pair transmission systems. Kahraman and Singh [2] assumed the shafts and bearings to be rigid and a lumpedparameter model including timevarying mesh stiffness, transmission error, and backlash was developed and solved with Harmonic Balance Method (HBM). Considering the coupled transversetorsional motion of a spur gear pair, a model of 3 DOF (degrees of freedom) was established by deducing the freedom of the dynamic model [3] and the dynamic equations were solved with HBM and RungeKutta method. Then a series of studies were done to reveal the nonlinear behavior of gearbearing system with backlash and timevarying mesh stiffness by Kahraman et al. [4, 5]. Cai and Hayashi [6] developed a 2 DOF model including timevarying mesh stiffness and transmission error and the calculated result with this 2 DOF linear model was consistent with the experiment result well. In [7], Cai proposed a mathematical method to simulate the timevarying mesh stiffness of helical gear.
Shen et al. [8] solved the singleDOF model of a spur gear pair considering the timevarying mesh stiffness, backlash, and static transmission error with Incremental Harmonic Balance Method (IHBM) and studied the effect of damping coefficient and excitation amplitude. Zhang et al. [9, 10] and Ma et al. [11] developed a model of multishaft helical gear considering geometric eccentricity. Cui et al. [12] studied the effect of revolution speed, backlash of gear pair, radial clearance of bearings, and shaft stiffness on the nonlinear dynamic behavior by a finite element/contact model. Li and Kahraman [13] proposed a tribodynamics model for spur gear pair which couples a mixed elastohydrodynamic lubrication model with a transversetorsional dynamic model. Liu et al. [14] focused on the nonlinear behavior of a spur gear with backlash by a 2 DOF lumpedparameter model, solved the dynamic equations with RungeKutta method, and expounded the influence of revolution speed, backlash, and mesh damping coefficient on the dynamic response. In [15], a dynamic model of a spur gear pair is employed to investigate the influence of gear tooth indexing errors on the dynamic response. In [16], philosophy of mathematical phenomenological mapping has been applied to the nonlinear dynamics of spur gears and radial ball bearings. Reference [17] describes a nonlinear dynamic model for the study of the vibration signals generated by gear transmissions. Liu et al. [18] examine the dynamic response of the spur gear pair system and interactions between bearing clearance and the backlash. In [19], a multidegreesoffreedom model is presented mathematically to characterize the static transmission error based on fractal theory, investigate the relative sliding friction using an EHLbased friction coefficient formula, and detail the timevarying friction coefficient suitable for HCR gear.
In most studies above, the excitation is usually given under nonmodulation condition. But, in the actual measurements, vibration and noise spectrums of gear pairs often exhibit sidebands around the gear mesh harmonic orders. Blankenship and Singh [20] revealed the modulation sidebands of a spur gear pair by a singledegreeoffreedom model. In [21], the AMFM modulation characteristic of gear box vibration under nonstationary operating conditions was studied.
To the authors’ knowledge, the internal excitation of the spur gear pair caused by the timevarying meshing stiffness and transmission error was usually regarded as nonmodulated signals. In this study, a 6 DOF lumpedparameter model of a spur gear pair considering bearing stiffness, timevarying mesh stiffness, backlash, and transmission errors is established in Section 2.1. The modulation internal excitation that caused the modulation sidebands is discussed in Section 2.2. Numerical simulation is made both under the nonmodulation condition and under modulation condition using RungeKutta method in Section 3. Then the calculated results under two conditions are compared. At the end of this study, an experiment is proposed to verify the modulation sidebands predicted by the model which is proposed in Section 2.1.
2. Dynamic Model of Spur Gear Pair
2.1. LumpedParameter Model
Here, some assumptions are presented to simplify the dynamic model: (a) friction forces at the gear mesh point are neglected, (b) the bearings and shafts that support the gears are modeled by equivalent elements with viscous damping and linear spring, and (c) the mesh damping is modeled by equivalent element with viscous damping. The simplified mechanical model to predict the modulation sidebands of a spur gear pair is established in Figure 1.
As shown in Figure 1, / and / are the translational displacements of the driving/driven gear, / is the torsional displacement of the driving/driven gear, / is the base circle radii of the driving/driven gear, andα is the pressure angle.
The relative gear mesh displacement is defined as where is the transmission error. The backlash nonlinearity function, , is defined as where is the halfbacklash value.
Define the dynamic gear mesh force aswhere is the mesh damping of the spur gear pair system and is the timevarying mesh stiffness.
Equations of multidegreesoffreedom motion of spur gear pair with lumpedparameter method are written as follows:where is the mass of driving/driven gear, is the rotational inertia of driving/driven gear, is the damping of driving/driven gear, is the bearing support stiffness, and is the input/output torque of the spur gear pair system.
2.2. Excitation of the Spur Gear Pair System
The input torque and output torque are assumed to be constant; namely, and =. The highfrequency internal excitation of the timevarying mesh stiffness and transmission error caused by transmission error is considered in the formulations as shown in (5)–(8c).
In the nonmodulation form, the timevarying mesh stiffness and the static transmission error e(t) are usually written as follows:
In (5), is the mesh frequency, is the mean mesh stiffness value, is the dimensionless jth harmonic amplitude coefficient, and is the dimensionless jth phase angle. Meanwhile, in (6), is the mean error displacement, is the amplitude coefficient of the jth harmonic terms, and is the dimensionless th phase angle.
If the system was to be exercised with the internal excitation defined as (5) and (6), the dynamic mesh force would observe little or no modulation sidebands and as a result the nonlinear response will observe little or no modulation sidebands.
The modulation timevarying mesh stiffness is rewritten as (7a), (7b), and (7c).
Here, are the dimensionless th modulation coefficients of amplitude/frequency, which are defined as follows:
In (7b), is the modulation amplitude coefficient of the th harmonic term of driving/driven gear and is the corresponding phase angle. For (7c), is the modulation frequency coefficient of the jth harmonic term of driving/driven gear and is the corresponding phase angle.
Like (7a), (7b), and (7c), the static transmission error e(t) is rewritten as where and are the modulation amplitude coefficient and modulation frequency coefficient of th harmonic term of the driving/driven gear and , , , and are the corresponding phase angles.
3. Numerical Simulation
Numerical simulations are carried out to investigate the modulation sidebands of the spur gear pair and a pair of spur gears are taken as an example and the system parameters of gear transmission are shown in Table 1. The RungeKutta method is used for solving the dynamic equations given in (4) with the modulation transmission errors given in (7a), (7b), and (7c) and modulation timevarying mesh stiffness given in (8a), (8b), and (8c).

Neglecting gear mesh damping and bearing support damping and considering the timevarying mesh stiffness as a constant , with the values of basic design parameters defined in Table 1, the six undamped natural modes of the lumpedparameter model of spur gear pair are determined as rigidbody mode at 0 Hz, gear pair mode at 1149 Hz, translational mode of driven gear at 1293 Hz, gear pair mode at 1604 Hz, translational mode of driving gear at 1799 Hz, and gear pair mode at 5043 Hz, and the corresponding natural mode types are shown in Figure 2.
(a)
(b)
(c)
(d)
(e)
(f)
For further calculation, the mean mesh stiffness and mean transmission error are defined as shown in Table 1. Considering only the fundamental wave of and , the corresponding modulation parameters in (7a)–(8c) are defined as , , , , , , , , , , , , , , , , , and .
The dynamic equations are solved where the input revolution speed is used as a control parameter. The relationship between amplitude of displacement along axis of driving gear and input revolution speed is shown in Figure 3.
(a)
(b)
The calculated undamped resonance frequency is mostly distributed in 1.2–1.8 kHz, and the corresponding input revolution speed of driving gear is distributed in 1000–2000 r/min. It can be observed that two peaks appear at 1250 r/min and 1750 r/min both in Figures 3(a) and 3(b). As shown in Figure 3, the resonant frequency bands under modulation condition are broader than corresponding ones under nonmodulation condition.
For better clarity, the response of the spur gear pair under offresonance conditions at the input revolution speed of 1500 r/min is further expounded in Section 3.1 and the response under nearresonance conditions at the input revolution speed of 1800 r/min is further expounded in Section 3.2.
3.1. Response under OffResonance Condition
Define the internal excitation as . The transmission error and timevarying mesh stiffness are predicted using (5)(6) under nonmodulation condition and using (7a)–(8c) under modulation condition. The two kinds of internal excitation are calculated and the results are shown in Figure 4. Here, axis presents frequency and axis presents the relative internal excitation defined as . Here, is the mean value of .
(a)
(b)
While the input revolution speed is 1500 r/min, only two peaks appear at 1375 Hz and 2750 Hz in the internal excitation spectrum under nonmodulation condition as shown in Figure 4(a). For the modulation condition, as shown in Figure 4(b), the modulation sidebands at Hz and Hz are obvious. In addition, two peaks at input revolution frequency of Hz and output revolution frequency of Hz are obvious.
To study the nonlinear dynamic behaviors of spur gear pair under offresonance condition at revolution speed of 1500 r/min, the internal excitation under nonmodulation condition as shown in Figure 4(a) and that under modulation condition as shown in Figure 4(b) are applied to the system. The response of the input bearing along axis under nonmodulation condition is shown in Figure 5(a) and the response of the input bearing along axis under modulation condition is shown in Figure 5(b). Under the nonmodulation condition, four peaks appear in the speed spectrum at the mesh harmonic order of Hz, Hz, Hz, and Hz and the peak at Hz is the max one. The sidebands disappear in the speed spectrum as shown in Figure 5(b). Under the modulation condition, sidebands appear around the mesh harmonic order of Hz, Hz, Hz, and Hz.
(a)
(b)
3.2. Response under NearResonance Condition
For revealing the modulation response of the spur gear pair under nearresonance condition, the input revolution speed of 1800 r/min is applied to the gearbearing system. The internal excitation under nonmodulation form and that under modulation form are shown in Figure 6.
(a)
(b)
At the revolution speed of 1800 r/min, the spur gear pair is under nearresonance condition as shown in Figures 6(a) and 6(b). Both internal excitation under nonmodulation condition and that under modulation condition at revolution speed of 1800 r/min are similar to the corresponding ones at revolution speed of 1500 r/min. For the nonmodulation condition, only two peaks appear in the spectrum at Hz and Hz as shown in Figure 6(a) and the sidebands disappear. For the modulation condition, sidebands appear around the mesh harmonic order of Hz and Hz. In addition, two more peaks appear in the spectrum at input revolution frequency of Hz and output revolution frequency of Hz.
To reveal the nonlinear dynamic response of the spur gear pair under nearresonance conditions, the internal excitation shown in Figure 6(a) and that shown in Figure 6(b) are applied to the system. The dynamic response of driving gear bearing along axis under nonmodulation condition is shown in Figure 7(a) and the one under modulation condition is shown in Figure 7(b).
(a)
(b)
For the nonmodulation condition, the frequency components under nearresonance condition are more complex than the ones under offresonance condition. As shown in Figure 7(a), the peaks appear at mesh harmonic orders of Hz, Hz, Hz, Hz, and Hz and the main peak appears at Hz. In addition, some of the energy of the gear mesh harmonic order at frequency of Hz, Hz, Hz, and Hz is distributed to small sidebands. For the modulation condition, four sidebands around mesh harmonic order at frequency of Hz, Hz, Hz, and Hz appear in the speed spectrum of the driving gear bearing along axis as shown in Figure 7(b) and the main peak appears at frequency of Hz. The peak value appearing at mesh frequency under modulation condition at revolution speed of 1800 r/min is much larger than the corresponding one under nonmodulation condition. The sideband around gear mesh harmonic order at frequency of Hz is not symmetric to mesh frequency.
4. Experimental Validation
An experiment is done to verify the prediction of the modulation sidebands. The parameters of the spur gear pair in the experiment are shown in Table 1. Figure 8(a) shows the test gear box of a spur gear pair and Figure 8(b) shows the spur gear pair with the parameters given in Table 1. As shown in Figure 8(a), use a frequency converter to control the driving motor and a magnetic particle brake is used as the loader. The rated power of the drive motor is 0.55 kW. The motor input current frequency ranges from 0 to 60 Hz. Motor synchronous speed ranges from 0 to 1800 r/min. The maximum braking current of the magnetic brake is 0.5 A and the maximum braking torque is 5 N·m. A revolution speed transducer is used for measuring the revolution speed of driving gear and an acceleration sensor is mounted on the gear box of the input bearing.
(a)
(b)
Since the maximum synchronous speed of the motor is only 1800 r/min, the experiment only verifies the system dynamics response at 1500 r/min. The braking current is 0.4 A. The sampling frequency of this experiment is 65536 Hz and the sampling time is 1 s. The measured speed spectrum of input bearing at revolution speed of 1498 r/min is shown in Figure 9.
As shown in Figure 5(b), the sidebands at frequency of , 2, 3, and 4 can be resolved in the calculation results. However, as shown in Figure 9, the sidebands of the measured spectrum only appear at frequency of and 2. There are two significant peaks at the meshing frequency and its 2x multiplier. The maximum peak appears at the meshing frequency of 1373 Hz. There are eight significant peaks near the meshing frequency, at 1330 Hz, 1355 Hz, 1348 Hz, 1366 Hz, 1380 Hz, 1391 Hz, 1398 Hz, and 1416 Hz, respectively.
As shown in Figure 5(b), there is a significant peak at the meshing frequency of 1375 Hz in the calculation result. There are also eight significant peaks near the meshing frequency, at 1332 Hz, 1350 Hz, 1357 Hz, 1368 Hz, 1382 Hz, 1393 Hz, 1400 Hz, and 1418 Hz, respectively. At the frequency of 2, both the calculated and the measured values show a significant peak. However, at the frequency of 2, neither the calculated spectrum nor the measured spectrum shows significant symmetry. The measured peak value is slightly larger than the calculated one. The theoretical value and experimental values are basically consistent.
5. Conclusion
In this study, a 6 DOF dynamic model is established to predict modulation sidebands of a spur gear pair due to the modulation internal excitation. The spectrum of internal excitation is shown to contain sidebands because of the modulation timevarying mesh stiffness and modulation transmission errors caused by manufacturing errors and so forth. The sidebands of modulation internal excitation spectrum mainly appear around the mesh harmonic order at frequency of and 2. Then the modulation internal excitation is applied to the spur gear pair and the dynamic equations are solved by RungeKutta method. As a result, the dynamic responses contain sidebands. At the end, the measured spectrums of spur gear pair acceleration are provided for comparisons to verify the dynamic responses predicted by the 6 DOF model proposed in this paper. The calculated results are basically in conformity with the measured results.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This project is supported by the National Natural Science Foundation of China (no. 51675350), the PhD Startup Foundation of Liaoning Institute of Science and Technology (no. 1604B03), and the University Science Research Project of Liaoning Education Department (no. L2016lky02).
References
 H. N. Özgüven and D. R. Houser, “Mathematical models used in gear dynamics—a review,” Journal of Sound and Vibration, vol. 121, no. 3, pp. 383–411, 1988. View at: Publisher Site  Google Scholar
 A. Kahraman and R. Singh, “Nonlinear dynamics of a spur gear pair,” Journal of Sound and Vibration, vol. 142, no. 1, pp. 49–75, 1990. View at: Publisher Site  Google Scholar
 A. Kahraman and R. Singh, “Nonlinear dynamics of a geared rotorbearing system with multiple clearances,” Journal of Sound and Vibration, vol. 144, no. 3, pp. 469–506, 1991. View at: Publisher Site  Google Scholar
 A. Kahraman and R. Singh, “Interactions between timevarying mesh stiffness and clearance nonlinearities in a geared system,” Journal of Sound and Vibration, vol. 146, no. 1, pp. 135–156, 1991. View at: Publisher Site  Google Scholar
 G. W. Blankenship and A. Kahraman, “Steady state forced response of a mechanical oscillator with combined parametric excitation and clearance type nonlinearity,” Journal of Sound and Vibration, vol. 185, no. 5, pp. 743–765, 1995. View at: Publisher Site  Google Scholar
 Y. Cai and T. Hayashi, “Linear approximated equation of vibration of a pair of spur gears (theory and experiment),” Journal of Mechanical Design, Transactions Of the ASME, vol. 116, no. 2, pp. 558–564, 1994. View at: Publisher Site  Google Scholar
 Y. Cai, “Simulation on the rotational vibration of helical gears in consideration of the tooth separation phenomenon (a new stiffness function of helical involute tooth pair),” Journal of Mechanical Design, vol. 117, no. 3, pp. 460–469, 1995. View at: Publisher Site  Google Scholar
 Y. Shen, S. Yang, and X. Liu, “Nonlinear dynamics of a spur gear pair with timevarying stiffness and backlash based on incremental harmonic balance method,” International Journal of Mechanical Sciences, vol. 48, no. 11, pp. 1256–1263, 2006. View at: Publisher Site  Google Scholar
 Y. Zhang, Q. Wang, H. Ma, J. Huang, and C. Zhao, “Dynamic analysis of threedimensional helical geared rotor system with geometric eccentricity,” Journal of Mechanical Science and Technology, vol. 27, no. 11, pp. 3231–3242, 2013. View at: Publisher Site  Google Scholar
 Y.M. Zhang, J. Yang, P. Hu, and Q.B. Wang, “Meshing characteristics analysis of spur gear pair considering modification coefficient,” Journal of Northeastern University (Natural Science), vol. 34, no. 9, pp. 1287–1291, 2013. View at: Google Scholar
 H. Ma, L. Zhu, Q. Wang, and Y. Zhang, “Modal coupling characteristic analysis of a helical gear rotor system with parallel shafts,” Proceedings of the Chinese Society of Electrical Engineering, vol. 32, no. 29, pp. 131–136, 2012. View at: Google Scholar
 L. Cui, X.G. Song, and J.R. Zheng, “Nonlinear dynamic analysis for a geared flexible rotor system considering multiclearance,” Journal of Vibration and Shock, vol. 32, no. 8, pp. 171–198, 2013 (Chinese). View at: Google Scholar
 S. Li and A. Kahraman, “A tribodynamic model of a spur gear pair,” Journal of Sound and Vibration, vol. 332, no. 20, pp. 4963–4978, 2013. View at: Publisher Site  Google Scholar
 J. Liu, S. Wang, S. Zhou et al., “Nonlinear behavior of a spur gear pair transmission system with backlash,” Journal of Vibroengineering, vol. 16, no. 8, pp. 3850–3861, 2014. View at: Google Scholar
 M. Inalpolat, M. Handschuh, and A. Kahraman, “Influence of indexing errors on dynamic response of spur gear pairs,” Mechanical Systems and Signal Processing, vol. 6061, pp. 391–405, 2015. View at: Publisher Site  Google Scholar
 I. Atanasovska, “The mathematical phenomenological mapping in nonlinear dynamics of spur gear pair and radial ball bearing due to the variable stiffness,” International Journal of NonLinear Mechanics, vol. 73, pp. 114–120, 2015. View at: Publisher Site  Google Scholar
 A. FernándezdelRincón, M. Iglesias, A. deJuan, A. DiezIbarbia, P. García, and F. Viadero, “Gear transmission dynamics: effects of index and run out errors,” Applied Acoustics, vol. 108, pp. 63–83, 2016. View at: Publisher Site  Google Scholar
 Z. Liu, Z. Liu, J. Zhao, and G. Zhang, “Study on interactions between tooth backlash and journal bearing clearance nonlinearity in spur gear pair system,” Mechanism and Machine Theory, vol. 107, pp. 229–245, 2017. View at: Publisher Site  Google Scholar
 K. Huang, Y. Xiong, T. Wang, and Q. Chen, “Research on the dynamic response of highcontactratio spur gears influenced by surface roughness under EHL condition,” Applied Surface Science, vol. 392, pp. 8–18, 2017. View at: Publisher Site  Google Scholar
 G. W. Blankenship and R. Singh, “Analytical solution for modulation sidebands associated with a class of mechanical oscillators,” Journal of Sound and Vibration, vol. 179, no. 1, pp. 13–36, 1995. View at: Publisher Site  Google Scholar
 F. Chaari, W. Bartelmus, R. Zimroz, T. Fakhfakh, and M. Haddar, “Gearbox vibration signal amplitude and frequency modulation,” Shock and Vibration, vol. 19, no. 4, pp. 635–652, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Zhong Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.