Abstract
In order to investigate the nonlinear characteristics of gear transmission system under the action of external and internal excitations, a dynamic model of a spur gear pair was established involving the backlash, damping, transmission error and the meshing stiffness. Based on the incremental harmonic balance method (IHBM), the general forms of the periodic solution with arbitrary precision are deduced. The vibration response obtained by IHBM compare very well with the results obtained by NewMark method, which verifies the accuracy and electiveness of the analytical methodology (IHBM) and provide information on the dynamic characteristic of spur gear. The simulation results revealed that several types of steadystate periodic solution are identified and determined by employing the IHBM. Due to the effect of backlash, the nonlinear characteristics of jump discontinuity phenomena and multiple stable solutions coexist and the meshing impact phenomenon are obvious. In addition, the influences of the system damping, transmission error and excitation amplitude on the amplitude frequency characteristic are illustrated by a series of diagrams. The results implicate that increasing the external excitation amplitude and system damping can effectively decrease the system resonant amplitude and control the nonlinear vibration response of the gear system and the effect of hardening spring behavior becomes weaker. The transmission error excitation amplitude variation also tends to worse the degree of nonlinearity. Therefore, it presents some useful information to reduce the vibration and noise of gear system.
1. Introduction
It is well known that the gear mechanisms have the advantages of wide range power, high torque to weight ratios, high reliability and stationarity, high transmission efficiency, and accurate drive ratio. Therefore, in the manufacture of industrial rotational machinery, gear transmission systems have been widely used in various applications such as wind turbines, ships, automobiles, and aircrafts. However, the special geometrical characteristics of the gear systems affect the dynamics and vibrational behavior in a significant way. The dynamic analysis of gear system is essential in describing noise and vibration characteristics. Due to the spur gear meshing, the system has different vibration characteristics to compare with simplified rotor system. Therefore, it is important to establish an accurate model for dynamic characteristics of the gear transmission system. In recent year, many researches have been reported on the dynamic models to analyze the dynamic characteristics of gear systems. As a consequence, research in the area of dynamics of mechanical system involving gear mechanisms has been intensified. A number of models were proposed in past to describe the dynamic behavior of gear system [14]. Ozguven [5] has comprehensively reviewed the mathematical models and dynamic phenomena in gear dynamics, where the focus was mainly on linear case. Chen [6] developed a nonlinear dynamic model of spur gear pair system with consideration of timevarying mesh stiffness and backlash, which was solved using the numerical method. Ma [7] established a general total degree of freedom dynamic model of a helical gear pair by considering the static transmission error and geometric eccentricities of the gears, and the system was simulated by using beam element and the gear pair by using mesh stiffness matrix and damping matrix.
With the development of nonlinear dynamics, some special issues for nonlinear of gear system will arise from numerical method. Therefore, special measures and small step size are needed which make parametric studies extremely expensive. The IHBM is a powerful approach, which is able to deal with strongly nonlinear system to any desired accuracy, and it is ideally suited to large range parametric studies [8]. A. Kaharman and R. Singh [9] studied the nonlinear dynamic characteristics of spur gear system with the harmonic balance method (HBM). In order to further study the nonlinear characteristics of gear transmission system, A. Kaharman [10] deduced the nonlinear dynamic equation of spur gear rotor system and the influences of various parameters were studied, which solved the nonlinear dynamic response of spur gear system with HBM and RungeKutta. Due to the effect of timevarying mesh stiffness, A. Kaharman [1112] sequentially analyzed the influence of the gear rotorbearing system with the backlash and the timevarying meshing stiffness, and carried on a detailed analysis to the nonlinear system. Zhang [13] performed dynamic analysis on gear system with the effect of timevarying mesh stiffness and friction between tooth pairs by HBM. In addition, the results obtained by HBM method compare very well with the results obtained by numerical integration method. Raghothama [14] investigated the periodic motions of a nonlinear geared rotorbearing system by the IHBM. A single degree of freedom model was proposed which included nonlinearities associated with periodic timevarying stiffness, backlash and viscous damping between a spur gear pair. Several key issues such as nonlinear modal interactions and differences between internal static transmission error excitation and external torque excitation were discussed by Li [15]. Shen [1619] comprehensively studied the periodic response of a single degree of freedom spur gear model with/without static transmission error by using IHBM, and demonstrated the effectiveness of increasing the damping coefficient ratio or excitation amplitude as a means of controlling the dynamics of the system. Yang [20] established a rightangle gear pair nonlinear timevarying dynamic model by considering both backlash and asymmetric mesh effects and the dynamic response was obtained using an enhanced multiterm harmonic balance method. It can be seen from foregoing literature that there are numerous studies on the topic, but the shafts and bearings are considered as rigid and actually vibration forms are negligible in a spur gear pait system. Due to the reduction of vibration and noise is still one of the issues in spur gear design, it is essential to accurately represent the different components of the gear system, so the relevant parameters ruling system’s dynamic behavior could be identified.
In this paper, considering the complex nonlinear vibration characteristic, an analytical model considered in the following sections proposes a spur gear pair that belongs to a partial of wind turbine gearbox, which can study the vibration characteristics of gear system systematically. And the present work investigates the complex, nonlinear dynamic behaviors of the system by discussing the influence of the system damping, the static transmission error and the excitation amplitude parameters on the periodic motions are illustrated by a series of response diagrams.
2. Dynamic model and equations of motion
2.1. Dynamic model of a spur gear
A simplified twodegreeoffreedom purely torsional generalized lumped parameter model of the spur gear transmission system investigated in the present study is shown in Fig. 1. Due to the effects of backlash, relevant parametric excitation and forced excitation, the gear system is a strongly nonlinear system. The shafts and bearings are considered as corresponding stiffness and damping. In addition, the centers of both gears are not allowed to move laterally. The mesh coupling between the driving and driven gears is represented using a set of mesh stiffness and a constant mesh damping elements acting along in the direction of the gear mesh line of action. In addition, backlash nonlinearity and static transmission error are also included.
In order to keep the model formulations at a manageable level, it is assumed that the meshing line direction is $y$direction in coordinate system in Fig. 1. ${O}_{1}$ and ${O}_{2}$ are the center of driving and driven gears. The torsional angular displacement of gear is assumed to result from a constant angular velocity term ${\omega}_{i}t$ ($i=$ 1, 2) plus a small variation displacement ${\theta}_{i}\left(t\right)$ due to vibrations originating from the flexibility of the mating gear teeth. Therefore, the angle displacements ${\phi}_{i}\left(t\right)$ of the driving and driven gears can be expressed by the following equations:
where ${\omega}_{1}$ and ${\omega}_{2}$ are the constant angular velocity components of the driving and driven gears.
Fig. 1Two degrees of freedom torsional dynamic model of spur gear system
2.2. Equations of motion
According to the previous geometrical relationships, as shown in Fig. 1, a displacement vector $\mathbf{X}$ of a spur gear pair can be defined from the pressure line coordinate system. So the generalized deformation vector can be expressed as follows:
From the proposed concept, taking into the mesh torque, input/output torque. In addition, the displacement force relation at the bearings is taken as linear in this paper. The kinetic energy $T$, the potential energy $U$ and the dissipation function $R$ are established. Utilizing the Lagrange equation, the equation of twodegreeoffreedom torsional motion of the nonlinear gear system can be expressed as:
${J}_{2}{\ddot{\theta}}_{2}+{c}_{t2}{\dot{\theta}}_{2}{r}_{b2}{c}_{m}\left({r}_{b1}{\dot{\theta}}_{1}{r}_{b2}{\dot{\theta}}_{2}\dot{e}\right)+{k}_{t2}{\theta}_{2}{r}_{b2}{k}_{m}f\left({r}_{b1}{\theta}_{1}{r}_{b2}{\theta}_{2}e\right)={T}_{l},$
where, ${J}_{1}$ and ${J}_{2}$ are the polar mass moment of inertias; the number of overdots represents the order of differentiation with respect to time $t$; ${\theta}_{i}$ ($i=$ 1, 2) is the torsional displacement of the driving/driven gear; ${r}_{b1}$ and ${r}_{b2}$ represent the base radii of gears, ${k}_{m}$ and ${c}_{m}$ are the mesh stiffness and mesh damping; ${k}_{ti}$ and ${c}_{ti}$ indicate the stiffness and damping of bearing; ${T}_{d}={T}_{dm}+{T}_{dr}\mathrm{s}\mathrm{i}\mathrm{n}({\omega}_{1}t+{\theta}_{d})$, ${T}_{l}={T}_{lm}+{T}_{lr}\mathrm{s}\mathrm{i}\mathrm{n}({\omega}_{2}t+{\theta}_{l})$ are the torque applied on driving and driven gears; $f\left(x\right)$ is a nonlinear displacement function of the dynamic transmission error $\delta $ and static transmission error $e$, which represents geometrical errors of the teeth profile and spacing. The dynamic transmission error $\delta $ and the new define variables $x$ can be expressed as follow:
while the nonlinear displacement function $f\left(x\right)$ is defined as:
where, ${e}_{m}$ and ${e}_{r}$ stand for the mean and fluctuation; ${\omega}_{e}$ is the gear meshing frequency, and ${\omega}_{e}=2\pi {n}_{1}{z}_{1}/60=2\pi {n}_{2}{z}_{2}/60$, ${\phi}_{e}$ is the initial phase of transmission error; ${n}_{1}$ and ${n}_{2}$ are rotational speed of gears; integer ${z}_{1}$ and ${z}_{2}$ represent the teeth number of driving and driven gears.
The backlash nonlinearity function $f\left(x\right)$ and illustration of impact is shown in Fig. 2, which presents three possible impacts, such as no impact, singlesided impact and doublesided impact.
Fig. 2Nonlinear displacement function
For simplicity, to derive the dimensionless equation of gear transmission system, the following transformation parameters are applied:
$\stackrel{~}{e}=\frac{1}{{b}_{c}}\left({e}_{m}+{e}_{0}{e}_{r}sin\left(\frac{{\omega}_{e}}{{\omega}_{0}}t+{\phi}_{e}\right)\right),$
where, ${b}_{c}$ is a characteristic length, ${e}_{0}$ is the coefficient of error fluctuation amplitude.
Therefore, nondimensionalized form of Eq. (3) is obtained by considering the foregoing manipulation:
${k}_{11}=\frac{{k}_{t1}}{{J}_{1}{\omega}_{0}^{2}},{k}_{22}=\frac{{k}_{t2}}{{J}_{2}{\omega}_{0}^{2}},$
${F}_{L2}=\frac{{r}_{b2}^{2}{k}_{m}}{{J}_{2}{b}_{c}{\omega}_{0}^{2}}\left({T}_{lm}+{f}_{{}_{0}}{T}_{lr}\mathrm{s}\mathrm{i}\mathrm{n}\left(\frac{{\omega}_{2}}{{\omega}_{0}}\tau +{\phi}_{l}\right)\right)\frac{{r}_{b2}^{2}{c}_{m}}{{J}_{2}{b}_{c}{\omega}_{0}^{2}}{e}_{0}{e}_{r}\mathrm{s}\mathrm{i}\mathrm{n}\left(\frac{{\omega}_{e}}{{\omega}_{0}}\tau +{\phi}_{e}\right),$
where, ${f}_{0}$ is the coefficient of excitationr fluctuation amplitude, ${\stackrel{~}{c}}_{m}={c}_{0}{c}_{m}$ and ${c}_{0}$ is the system damping coefficient.
3. Incremental harmonic balance method for MDOF gear system
For the multiple degrees of freedom (MDOF) gear transmission system, the generalized nonlinear equations can be written as follows:
where, $\mathbf{M}$, $\mathbf{C}$, $\mathbf{K}$ are mass, linear viscous damping and linear stiffness matrices of gear system, respectively; $\mathbf{f}(\dot{\mathbf{X}},\mathbf{X})$ is a nonlinear force; $\mathbf{X}={\left[{x}_{1},{x}_{2},\dots ,{x}_{i}\right]}^{T}$ is the degrees of freedom vector of the gear system; $i$ represents the number of the degrees of freedom; $\mathbf{P}={\left[{p}_{1},{p}_{2},\dots ,{p}_{i}\right]}^{T}$ is the external excitation of the gear system.
In this paper, efforts are devoted to the application of the IHBM to a multidegree of freedom (MDOF) nonlinear gear system. So $\tau =\omega t$, and the Eq. (9) can be expressed as follows:
where, $\omega $ is the excitation frequency (minimum rotational frequency) of gear system; the dots denote derivatives with respect to the dimensionless time $\tau $.
Using the NetwonRaphson procedure of the IHBM. Letting ${x}_{i0}$, $\omega $ denote a state of vibration in hand, the neighboring state can be expressed by adding the corresponding increments to them as follows:
where, ${\mathbf{X}}_{0}={\left[{x}_{10},{x}_{20},\dots ,{x}_{i0}\right]}^{T}$; $\mathrm{\Delta}{\mathbf{X}}_{0}={\left[\mathrm{\Delta}{x}_{10},\mathrm{}\mathrm{\Delta}{x}_{20},\dots ,\mathrm{\Delta}{x}_{i0}\right]}^{T}$.
If the gear system concerned is rotational frequency, meshing frequency, combinatorial frequency and so on. The steady state response can be expressed as follows [12]:
${x}_{i0}\left(\tau \right)={a}_{i00}+\sum _{n=1}^{N}\left({a}_{in0}\mathrm{cos}\left(n\tau \right)+{b}_{in0}\mathrm{sin}\left(n\tau \right)\right)={\mathbf{C}}_{s}{\mathbf{A}}_{i0},$
$\mathrm{\Delta}{x}_{i}\left(\tau \right)=\mathrm{\Delta}{a}_{i0}+\sum _{n=1}^{N}\left(\mathrm{\Delta}{a}_{in}\mathrm{cos}\left(n\tau \right)+\mathrm{\Delta}{b}_{in}\mathrm{sin}\left(n\tau \right)\right)={\mathbf{C}}_{s}\mathrm{\Delta}{\mathbf{A}}_{i},$
where, $N$ is the sum of $N$ order harmonic terms:
Substituting Eq. (12) into Eq. (11), the vector of the gear system can be expressed as follows:
where:
Substituting Eq. (13) into Eq. (10) and neglecting small terms of high order, one obtains the following equation:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\mathbf{P}\left(\tau \right){\mathbf{f}}_{n}\left({\mathbf{X}}_{0},{\dot{\mathbf{X}}}_{0}\right)\left({\omega}_{0}^{2}\mathbf{M}{\ddot{\mathbf{X}}}_{0}+{\omega}_{0}\mathbf{C}{\dot{\mathbf{X}}}_{0}+\mathbf{K}{\mathbf{X}}_{0}\right)\left(2{\omega}_{0}\mathbf{M}{\ddot{\mathbf{X}}}_{0}+\mathbf{C}{\dot{\mathbf{X}}}_{0}\right)\mathrm{\Delta}\omega ,$
where, ${\mathbf{C}}_{n}=\partial {\mathbf{f}}_{n}(\dot{\mathbf{X}},\mathbf{X})/\partial \dot{\mathbf{X}}$, ${\mathbf{K}}_{n}=\partial {\mathbf{f}}_{n}(\dot{\mathbf{X}},\mathbf{X})/\partial \mathbf{X}$, ${\mathbf{C}}_{n}$ and ${\mathbf{K}}_{n}$ are the Jacobian matrixes.
Using the Galerkin procedure gives:
$={\int}_{0}^{2\pi}\delta \left(\mathbf{X}{)}^{T}\right[\mathbf{P}\left(\tau \right){\mathbf{f}}_{n}({\mathbf{X}}_{0},{\dot{\mathbf{X}}}_{0})({\omega}_{0}^{2}\mathbf{M}{\ddot{\mathbf{X}}}_{0}+{\omega}_{0}\mathbf{C}{\dot{\mathbf{X}}}_{0}+\mathbf{K}{\mathbf{X}}_{0})$
$(2{\omega}_{0}\mathbf{M}{\ddot{\mathbf{X}}}_{0}+\mathbf{C}{\dot{\mathbf{X}}}_{0})\mathrm{\Delta}\omega ]d\tau .$
Substituting Eq. (13) into Eq. (15) and simplified by:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}={\int}_{0}^{2\pi}\left(\delta \mathbf{A}{)}^{T}{\mathbf{S}}^{T}\right[\mathbf{P}\left(\tau \right){\mathbf{f}}_{n}({\mathbf{X}}_{0},{\dot{\mathbf{X}}}_{0})({\omega}_{0}^{2}\mathbf{M}\ddot{\mathbf{S}}{\mathbf{A}}_{0}+{\omega}_{0}\mathbf{C}\dot{\mathbf{S}}{\mathbf{A}}_{0}+\mathbf{K}\mathbf{S}{\mathbf{A}}_{0})$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}(2{\omega}_{0}\mathbf{M}\ddot{\mathbf{S}}{\mathbf{A}}_{0}+\mathbf{C}\dot{\mathbf{S}}{\mathbf{A}}_{0})\mathrm{\Delta}\omega ]d\tau .$
According to the variational characteristics, the Eq. (16) can be expressed as follow:
$={\int}_{0}^{2\pi}{\mathbf{S}}^{T}\left({\omega}_{0}^{2}\mathbf{M}\ddot{\mathbf{S}}+{\omega}_{0}\mathbf{C}\dot{\mathbf{S}}+\mathbf{K}\mathbf{S}\right)d\tau \cdot {\mathbf{A}}_{0}+{\int}_{0}^{2\pi}{\mathbf{S}}^{T}\left[\mathbf{P}\left(\tau \right){\mathbf{f}}_{n}\left({\mathbf{X}}_{0},{\dot{\mathbf{X}}}_{0}\right)\right]d\tau $
${\int}_{0}^{2\pi}{\mathbf{S}}^{T}\left(2{\omega}_{0}\mathbf{M}\ddot{\mathbf{S}}+\mathbf{C}\dot{\mathbf{S}}\right)d\tau \cdot {\mathbf{A}}_{0}\mathrm{\Delta}\omega .$
Therefore, a set linear equations in terms of $\mathrm{\Delta}\mathbf{A}$ and $\mathrm{\Delta}\omega $ can be obtained readily:
where:
${\mathbf{R}}_{m1}={\int}_{0}^{2\pi}{\mathbf{S}}^{T}\left({\omega}_{0}^{2}\mathbf{M}\ddot{\mathbf{S}}+{\omega}_{0}\mathbf{C}\dot{\mathbf{S}}+\mathbf{K}\mathbf{S}\right)d\tau ,$
${\mathbf{R}}_{m2}={\int}_{0}^{2\pi}{\mathbf{S}}^{T}\left[\mathbf{P}\right(\tau ){\mathbf{f}}_{n}({\mathbf{X}}_{0},{\dot{\mathbf{X}}}_{0}\left)\right]d\tau ,{\mathbf{R}}_{m3}={\int}_{0}^{2\pi}{\mathbf{S}}^{T}\left(2{\omega}_{0}\mathbf{M}\ddot{\mathbf{S}}+\mathbf{C}\dot{\mathbf{S}}\right)d\tau .$
4. The validation of analytical methods
From the previous conclusion and analysis, it can be seen that the spur gear pair is a complicated system with the strong nonlinearity, time variance and complicated working environment. Therefore, it is necessary to give a detailed analysis of the gear system. The dynamic behaviors of system are investigated by IHBM. In this paper, Table 1 summarizes the geometrical and physical parameters of a spur gear pair.
Table 1System parameters of gear transmission
Case study  Driving/driven gear  Case study  Driving/driven gear 
Number of teeth ${z}_{1}/{z}_{2}$  20  Meshing stiffness ${k}_{m}$  5.0×10^{8} N/m 
Module $m$  10 mm  Meshing damping ${c}_{m}$  8.0×10^{2} N/(m/s) 
Moment of inertia ${J}_{1}/{J}_{2}$  0.05 kg.m^{2}  Error mean ${e}_{m}$  2.0×10^{5} m 
Radius ${r}_{b1}/{r}_{b2}$  0.1 m  Error fluctuation ${e}_{r}$  3.0×10^{5} m 
Torsional stiffness ${k}_{t1}/{k}_{t2}$  4.0×10^{6} N.m/rad  Torque mean ${T}_{dm}/{T}_{lm}$  100 N/m 
Torsional damping ${c}_{t1}/{c}_{t2}$  5.0×10^{2} N/(rad/s)  Torque fluctuation ${T}_{dr}/{T}_{lr}$  300 N/m 
In order to demonstrate the effectiveness and accuracy of the IHBM, the nonlinear dynamic response characteristics of the spur gear pair are obtained by adopting the IHBM and numerical integration NewMark method (NMM). The timedomain waveforms of the twodegreeoffreedom spur gear system in this study are shown in Fig. 3, where the circle denotes the steady state calculation results by the IHBM and the solid line denotes stable solutions by the NMM, respectively.
Fig. 3Timedomain waveform of the spur gear system
a)
b)
As shown in Fig. 3, it is indicated that the result by IHBM agrees with the NMM’s, which verifies the precision of the IHBM presented in this paper to have high enough. In addition, it takes 300 cycles and 78.84 s, when the gear system closes to steady state response of convergence by NMM. However, the IHBM only needs 13.51 s, which also demonstrate the effectiveness of the IHBM presented in this paper to have a significant advantage.
5. Nonlinear dynamic response of gear system
In this section, the influence of the backlash, the transmission error and the external excitation amplitude parameters on the periodic motions are investigated to detailedly understand the effects of key parameters on the gear dynamical response by a series of response diagrams.
5.1. The jump discontinuity phenomena and multiple stable solutions
The nonlinear dynamic characteristics of the gear system are detailedly obtained by the IHBM. Due to the driving gear and driven gear have the same geometrical parameters, one of the gear vibration response characteristics is given in this paper. The amplitude frequency characteristic curves of gear transmission system are shown in Fig. 4, where the circle denotes the amplitude frequency curve with considering the backlash ($b\ne $ 0) and the solid line denotes amplitude frequency curve without considering the backlash ($b=$ 0) by the IHBM. It can be seen from Fig. 4 that the gear transmission system is a strongly nonlinear system with the influence of backlash factor. There are obviously typical nonlinear characteristics between approximately linear system ($b=$0) and nonlinear system ($b\ne $ 0), and the gear system presents hardening spring behavior. As a consequence, the amplitude frequency curve in the solution branches, multiple stable solutions may coexist, making possible the appearance of classical jump discontinuity phenomena with the frequency sweeping. In addition, it appears resonance characteristics closeby the natural frequency and the response amplitude increases slightly.
Fig. 4Jump discontinuity phenomena and multiple stable solutions
5.2. Effect of system damping coefficient
System damping coefficient ${c}_{0}$ is an important parameter that it affects the frequency response substantially. To illustrate this point, in this section, it keeps other parameters constant and only change the system damping coefficient ${c}_{0}\text{.}$ The effect of system damping coefficient on bifurcation diagram and amplitude frequency curves are shown in Fig. 5. It can be seen from the bifurcation diagram that the dynamic behavior of the spur gear system performs chaotic at low values of the system damping coefficient ${c}_{0}<$0.96. At higher values of the dimensionless system damping coefficient, i.e. ${c}_{0}>$0.95, that the motion state of the system mainly shows quasiperiodic motion state. But when the dimensionless system damping coefficient ${c}_{0}>$2.0, the system has a change trend to 1Tperiodic motion. It can be noted from Fig. 5(b) that the gear system shows the rich nonlinear behaviors including obvious jump discontinuities phenomena and multiple stable solutions respectively. When the system damping coefficient ${c}_{0}$ is small, it can be found the amplitude of response increase and two jump discontinuities are seen coupled with singlesided impact and doublesided impact. In addition, the system exhibits hardening spring behavior due to contact. As the system damping coefficient ${c}_{0}$ increase, the nonlinear characteristic of jump discontinuity phenomena is disappearing and further increase of ${c}_{0}$ will decrease the response amplitude and multiple stable solutions withdraw gradually, which makes the vibration response behavior more linearly. In addition, the large amplitude of the primary resonance frequency is decreasing and the effect of hardening spring becomes weaker. Therefore, the effect of backlash nonlinearity can decrease with the increasing system damping coefficient ${c}_{0}$.
Fig. 5Dynamics response of system damping coefficient
a)
b)
5.3. Effect of static transmission error
The transmission error is a displacement excitation at the mesh point due to the geometrical error of the teeth profile and spacing, which not only affects the response amplitude of the system, but also directly affects the degree of nonlinearity. To examine the effect of transmission error on the frequency response, in this section, the bifurcation diagram with the changing transmission error coefficient ${e}_{0}$ is shown in Fig. 6(a). And three different values of the transmission error coefficient ${e}_{0}$ are investigated and the results are shown in Fig. 6(b).
Fig. 6Dynamics response of system with changing meshing error
a)
b)
It can be seen from Fig. 6(a) that the gear system exhibits 1Tperiodic motion at low values of the dimensionless transmission error ${e}_{0}$, i.e., ${e}_{0}<$0.48. At higher values of the transmission error coefficient ${e}_{0}$, i.e. 0.48$<{e}_{0}<$0.79, that the motion state of the system mainly shows quasiperiodic motion state. However, as ${e}_{0}$ is increased, the $n$Tperiodic motion is replaced by nonperiodic motion. As shown in Fig. 6(b), when the transmission error coefficient ${e}_{0}$ is larger, the system clearly performs hardening spring behavior, multiple stable solutions and jump discontinuity phenomena when the gear pair undergoes the transition from singlesided impact to doublesided impact behaviors with the frequency sweping down. However, it can be seen that the jump discontinuity phenomena, tooth impact behavior and multiple stable solutions are gradually disappearing with the decreasing of the transmission error coefficient ${e}_{0}$. In addition, the response amplitude also decreases. Hence, it should reduce the gear transmission error between gears in manufacturing and design phase of gear, which could decrease vibration and noise of gear transmission system.
5.4. Effect of external excitation amplitude
In this analysis, only the effect of the external excitation amplitude is analyzed and the other parameters remain the same. The bifurcation diagram and the amplitude frequency curves are shown in Fig. 7. Fig. 7(a) presents the bifurcation diagram of the spur gear system using the dimensionless external excitation amplitude ${f}_{0}$ as a bifurcation parameter. A transition from periodic motion to the chaotic dynamic are seen when ${f}_{0}$ is increase from 0 to 2.68. It can be observed from Fig. 7(a) that the gear system exhibits a quasiperiodic response at low values of the external excitation amplitude, i.e., ${f}_{0}<$0.5. However, as ${f}_{0}$ is increase from 0.5 to 2.68, the quasiperiodic response is replaced by chaotic motion. As the control parameter ${f}_{0}$ is further increased, i.e. ${f}_{0}>$2.68, the chaotic motion transits to 1Tperiodic motion. Fig. 7(b) presents the results for the spur system with three different values of the external excitation amplitude parameter ${f}_{0}$. As the ${f}_{0}$ increases, significant changes in the amplitude frequency curves are observed. When the ${f}_{0}$ is relatively small values (i.e. light load), the gear system gradually transition between no impact, singlesided impact and doublesided impact with the frequency sweeping respectively. In addition, the amplitude frequency curves exhibit strongly hardening spring behavior and the response amplitude smaller. As the amplitude of external excitation parameter ${f}_{0}$ increases (i.e. heavy load), the singlesided impact and doublesided impact gradually withdraw. Namely, the jump discontinuity phenomena and multiple stable solutions also disappear gradually. When the amplitude of excitation parameter ${f}_{0}$ increases to a certain value, the characteristics of gear system response approximate to a linear system and there is no jump discontinuity phenomena seen in the case at all, but the response amplitude increases obviously. Therefore, it is concluded that increasing the amplitude of external excitation can effectively control the nonlinear dynamic response characteristics of the gear transmission system. However, the external excitation amplitude should be limited in a certain scope because of its dual character.
Fig. 7Dynamics response of system with changing excitation amplitude
a)
b)
6. Conclusions
The nonlinear dynamic characteristic of a spur gear pair is studied considering parametrically excited and piecewise linear of two degrees of freedom system that includes the backlash nonlinearity, static transmission error and external excitation. The amplitude frequency characteristic analysis is carried out by the IHBM under a complex dynamical condition. The system shows nonlinear resonances with hardening behavior, jump and multiplicity. Based on the results of the parametric study presented in the previous section, the conclusions can be summarized as follows:
1) Due to the effects of the backlash and the internal and external excitations, the gear system exists obviously typical nonlinear characteristics with hardening spring behavior, jump discontinuity phenomena, such as singlesided impact, doublesided impact, and multiple stable solutions and then the gear pair system can also exhibit more complicated dynamic response.
2) The effect of the system damping parameter is verified that an increase in its value causes a decrease in the amplitude of the solution and the jump discontinuities phenomena and multiple stable solutions withdraw gradually. In addition, the large amplitude of the primary resonance frequency is decreasing and the hardening spring behavior will fade away gradually. Error excitation amplitude variation also tends to worse the degree of nonlinearity, i.e. the nonlinear characteristics are more and more obvious by increasing the error excitation amplitude and the response amplitude increases slightly. The effect of hardening spring behavior gradually strengthens. The unpredictable nonlinear characteristics of jump and multiplicity shoule be avoided in the actual engineering, so it try to increase the system damping system and decrease the installation error and manufacture error during the design phase of gear system.
3) The excitation amplitude has different level effects on amplitude frequency characteristics of the system. The jump discontinuity phenomena and multiple stable solutions of the system become very obvious with smaller excitation amplitude. As the increase of excitation amplitude, the doublesided impact and singlesided impact are disappearing gradually. i.e. nonlinear characteristic becomes weaker, but the response amplitude increases obviously. Therefore, the excitation amplitude should be limited in a certain scope because of its dual character and then preferably control of the dynamic characteristic of gear transmission system.
4) Increasing load, systen dmping and decreasing the error can reduce the amplitude of vibration, the increased complexity of the dynamical behaviour with a spur gear suggests further studies. In particular, a global optimization of gear should be developed in order to suggest the best geometry of spur gears versus vibration reduction.
References

Li R. F., Wang J. J. Gear System Dynamics: Vibration, Impact, Noise. Beijing, Science Press, 1997, p. 13, (in Chinese).

Wojnarowski J., Onishchenko V. Tooth wear effects on spur gear dynamics. Mechanism and Machine Theory, Vol. 38, Issue 2, 2003, p. 161178.

Kuang J. H., Lin A. D. The effect of tooth wear on the vibration spectrum of a spur gear pair. Journal of Vibration and Acoustics, Vol. 123, Issue 3, 2001, p. 311317.

Wang J. J., Li R. F., Peng X.H. Survey of nonlinear vibration of gear transmission system. Applied Mechanics Reviews, Vol. 56, Issue 3, 2003, p. 309329.

Ozguven H. N., Houser D. R. Mathematical models used in gear dynamicsa review. Journal of Sound and Vibration, Vol. 121, Issue 3, 1988, p. 383411.

Chen H. T., Wu X. L., Qin D. T., et al. Dynamic characteristics of gear transmission system subjected to random internal and external excitation. China Mechanical Engineering, Vol. 24, Issue 4, 2013, p. 233537, (in Chinese).

Ma H., Wang Q. B., Huang J. Vibration response analysis of gear coupled rotor system considering geometric eccentric effect of helical gears. Journal of Aerospace Power, Vol. 28, 2013, p. 1624, (in Chinese).

Zhou J. X., Zhang L. Incremental harmonic balance method for predicting amplitudes of a multid.o.f. nonlinear wheel shimmy system with combined Coulomb and quadratic damping. Journal of Sound and Vibration, Vol. 279, Issue 12, 2005, p. 403416.

Kahraman A., Singh R. Nonlinear dynamic of a spur gear pair. Journal of Sound and Vibration, Vol. 142, Issue 1, 1990, p. 4975.

Kahraman A., Singh R. Nonlinear dynamic of geared rotorbearing system with multiple clearances. Journal of Sound and Vibration, Vol. 144, Issue 3, 1991, p. 469506.

Kahraman A., Singh R. Interactions between timevarying mesh stiffness and clearance nonlinear in a geared system. Journal of Sound and Vibration, Vol. 146, Issue 1, 1991, p. 135156.

Blankenship G. W., Kahraman A. Steady state forced response of a mechanical oscillator with combined parametric excitation and clearance type nonlinearity. Journal of Sound and Vibration, Vol. 185, Issue 5, 1995, p. 743756.

Zhang S. H., Li Y. P., Qiu D. M. Dynamic analysis of a rotorbearing system geared by means of the harmonic balance method. Chinese Journal of Mechanical Engineering, Vol. 36, Issue 7, 2000, p. 1822, (in Chinese).

Raghothama A., Narayana S. Bifurcation and chaos in geared rotor bearing system by incremental harmonic balance method. Journal of Sound and Vibration, Vol. 226, Issue 3, 1999, p. 469492.

Li Y. G., Chen T. N., Wang X. P., et al. Nonlinear dynamics of spur gear pair under external periodic excitation. Journal of Xi’An Jiao Tong University, Vol. 24, Issue 4, 2013, p. 233537, (in Chinese).

Yang S. P., Shen Y. J., Liu X. D. Nonlinear dynamics of gear system based on incremental harmonic balance method. Journal of Vibration and Shock, Vol. 24, Issue 3, 2005, p. 4043, (in Chinese).

Shen Y. J., Yang S. P., Liu X. D. 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, Issue 8, 2006, p. 12561263.

Shen Y. J., Yang S. P, Pan C. Z., et al. Nonlinear dynamics of a spur gear pair with timevarying stiffness and backlash. Journal of Low Frequency Noise, Vibration and Active Control, Vol. 23, Issue 3, 2004, p. 178187, (in Chinese).

Shen Y. J., Yang S. P, Pan C. Z., etc. Nonlinear dynamics of a spur gear pair to parametric and external excitations. Journal of Beijing Jiaotong University, Vol. 29, Issue 1, 2005, p. 6973, (in Chinese).

Yang J. Y., Peng T., Lim T. C. An enhanced multiterm harmonic balance solution for nonlinear periodone dynamic motions in rightangle gear pairs. Nonlinear Dynamics, Vol. 67, Issue 3, 2012, p. 10531065.
About this article
The project was supported by China Natural Science Funds (No. 51305276) and the Fundamental Research Funds for the Central Universities (No. N120403004).