0
Research Papers

# Nonlinear Vibration of Gears With Tooth Surface ModificationsOPEN ACCESS

[+] Author and Article Information
Tugan Eritenel

Mem. ASME
Department of Mechanical Engineering,
The Ohio State University,
Columbus, OH 43210
e-mail: eritenel.1@osu.edu

Robert G. Parker

Research Professor
Fellow ASME
The Ohio State University,
Columbus, OH 43210;
Distinguished Professor Chair and Executive Dean,
University of Michigan–Shanghai Jiao Tong University Joint Institute,
Shanghai Jiao Tong University,
Shanghai 200240, China
e-mail: parker.242@osu.edu

1Corresponding author.

Contributed by the Design Engineering Division of ASME for publication in the Journal of Vibration and Acoustics. Manuscript received June 30, 2011; final manuscript received May 6, 2012; published online June 18, 2013. Assoc. Editor: Philip Bayly.

J. Vib. Acoust 135(5), 051005 (Jun 18, 2013) (11 pages) Paper No: VIB-11-1149; doi: 10.1115/1.4023913 History: Received June 30, 2011; Revised May 06, 2012

## Abstract

This work provides an analytical solution for the nonlinear vibration of gear pairs that exhibit partial and total contact loss. Partial contact loss is where parts of contact lines lose contact although other parts remain in contact. The gear tooth surface modifications admit an arbitrary combination of profile and lead modifications. Modifications are a source of partial contact loss. The analysis also applies for total contact loss. Unlike models in the literature that are excited by static transmission error or time-varying mesh stiffness, the excitation and the nonlinearity are not a priori specified. Instead, the force-deflection function of the gear pair is provided by an independent source, such as a finite element model or Hertz contact formula. The manipulation of the single-degree-of-freedom oscillator equation of motion yields the excitation and the nonlinearity that arise from Fourier and Taylor series expansions of the force-deflection function. These expansions capture the essential contact behavior that includes tooth profile and lead modifications as well as the bending and shear flexibility of the gear teeth and gear blanks. The method of multiple scales gives the steady-state dynamic response in terms of a frequency-amplitude relation. Comparisons with gear vibration experiments and simulations from the literature that include spur and helical gears with tooth profile and lead modifications verify the method.

<>

## Introduction

Vibration of gear pairs can be considered under the general category of contact vibrations, which also includes vibration of bearings, splines, linkages, and other mechanical connections. This work focuses on nonlinear rotational vibration of gears near resonance speeds. Although gear resonance can be avoided in some constant speed applications, many gears operate over a wide range of speeds where some gear resonances cannot be avoided. In addition, design constraints may prevent changes in the system to move the natural frequencies away from the gear mesh frequency.

In gear applications, the involute tooth surface shape is usually modified slightly. Lead modifications are used to improve misalignment tolerance [1], and profile modifications are used to avoid undesirable edge or corner contact [2]. With or without such tooth modifications, gear vibrations exhibit softening nonlinearity [3-9] near resonance as a result of partial or total contact loss. Partial contact loss is where portions of nominally contacting surfaces lose contact (and other nominally noncontacting portions potentially gain contact) as a result of dynamic motions, profile and lead modifications, or misalignment. Partial contact loss includes reductions in both the profile and face contact ratios. Modifications are the major source of partial contact loss [10]. Although gear systems that exhibit total contact loss have been analyzed, only numerical solutions for the vibration of gears that exhibit partial contact loss are available in the literature; for example Refs. [4,11-15]. This work applies a perturbation method that yields a closed-form solution for nonlinear gear contact vibrations. The method admits spur and helical gears with tooth profile and lead modifications. The solution is compact and exposes the effects of design parameters on the dynamic response. Although harmonic balance can yield closed-form analytical solutions in some nonlinear systems, refined frequency and time discretizations are needed when applying the method to gear contact problems. As a result, the process becomes numerical without expressions showing parameter dependencies of the dynamic response.

The equations of motion for gear vibrations include time-dependent parameters such as periodic variation in mesh stiffness. In many cases, lumped-parameter gear models excited by static transmission error or time-varying mesh stiffness give satisfactory results [6,11,16-21]. A single spring is used to model the gear mesh interface in these works. Velex and Ajmi [21] examine the validity of approximating transmission error as the excitation source. Liu and Parker [22] explore the conditions under which the aforementioned approximations work.

Harmonic balance [23-26] and perturbation methods [27-29] are used when the gear vibration is approximated using a lumped-parameter model. The mesh stiffness formulation in Refs. [24-28] include the periodic fluctuation due to gear tooth engagement/disengagement, while Ref. [23] uses a constant mesh stiffness. The sole source of nonlinearity in the works listed above is total contact loss, a piecewise nonlinearity in which the gear mesh ceases to transmit any force. Gear tooth surface modifications are not included in these works. The modeling in Refs. [22,30] for the dynamic response of multimesh gears differs from the literature listed above because it considers tooth profile modifications and contact loss at each of the individual meshing gear teeth rather than the gear mesh as a whole. In Refs. [22,30], such contact loss is due to linear tooth profile modifications, which is the only type of modification they consider.

Detailed tooth contact models allow dynamic analysis for specified tooth profile and lead modification [3,4,9,31,32]. Such models predict partial contact loss arising from arbitrary tooth surface modifications [4,9,13,15,33,34]. The contact algorithms that allow for partial contact loss, however, are prohibitively complex for analytical methods, and nonlinear dynamic response can only be obtained numerically. Numerical methods give only restricted generalizations about dynamic behavior because they are limited to selected parametric studies. There are no studies in the literature that provide closed-form solutions for the nonlinear vibration of gears when partial contact loss is present.

This work gives closed-form solutions for nonlinear rotational gear vibrations near primary resonance, that is, when the excitation frequency is close to the natural frequency. The solution can accommodate spur and helical gears, and it includes the nonlinear behavior due to partial contact loss and admits arbitrary modifications of the gear tooth surface including profile and lead modifications. The use of a general force-deflection function as an input frees the dynamic equation of motion from the physical gear modeling assumptions and can apply to other contact problems. For gears, the force-deflection function captures the phenomena of partial contact loss, parameter time-dependence, tooth profile and lead modifications, changing number of teeth in contact, changing radius of curvature at contact, and the elastic deformation of the gear blank. With Taylor and Fourier series expansions of this force-deflection function, the equation of motion takes a form that enables use of the method of multiple scales.

## Mathematical Model

###### Equation of Motion.

This section considers a single degree-of-freedom equation of motion. The single degree-of-freedom model includes fluctuating mesh stiffness and mesh force nonlinearity. The nonlinear and periodically varying mesh force is represented by a general, time-dependent force-deflection curve that can account for partial contact loss, gear misalignment, tooth surface modifications, varying number of teeth in contact, torque dependence, and total contact loss. Mathematical manipulations of the tooth mesh force-deflection term using Fourier and Taylor series expansion puts the equation of motion in a form suitable for perturbation analysis. These manipulations and the closed-form perturbation approximations for a general mesh force-deflection relationship are what distinguish this work from past studies. The basic gear modeling has been proposed in past research. It is assumed that the force-deflection function and the applied load vary periodically with time. No further assumptions are made at this point.

The equation of motion isDisplay Formula

(1)$mx··+cx·+C(x,t)f(x,t)=F(t)$

where $F(t)$ is the periodic applied load with period $PF$, $f(x,t)$ is the nonlinear, time-dependent, periodic force-deflection function with period $Pf$, and $C(x,t)$ is the separation function given byDisplay Formula

(2)$C(x,t)={0if x

where $g(t)$ is the unloaded ($F(t)=0$), quasi-static ($x·=x··=0$) transmission error one would measure if the gears are rotated through a gear mesh cycle with zero torque. It satisfies $f(g(t),t)=0$.

This single degree-of-freedom mechanical oscillator approximates the rotational vibration of gear pairs with the following assumptions. The bearings are assumed to be rigid, allowing a model with only rotational degrees of freedom. This assumption is justified when the natural frequencies of vibration modes with high bearing deformations are sufficiently separated from the natural frequencies of vibration modes with primarily rotational oscillations of the gears.

The force-deflection function $f(x,t)$ is general and there are no restrictions on the type of the elastic deformation that can be included. In gears, the effects include, but are not limited to, Hertz contact deformations at the gear teeth, tooth bending, shear, and gear blank deformations. Because the number of teeth in contact changes with rotation of gears, the force-deflection function is periodic with tooth pass frequency. An important distinction must be made between misaligned gears in the original assembly and dynamic tilting motions that occur during vibrations. We use the term misalignment to refer to a mounting error such that the tooth load is unevenly distributed across the facewidth in its assembled state. Misalignment affects the mesh stiffness and would be reflected in the general force-deflection function, which can be calculated by finite element analysis. Therefore the analysis considers the effects of misalignment. The dynamic tilting motions that would require additional degrees-of-freedom are neglected.

Equation (2) restricts contact loss to single-sided impacts. Elimination of the rigid body motion [19,35] lumps the rotational inertias of both the pinion ($Ip$) and gear ($Ig$) into $m=IpIg/(Iprg2+Igrp2)$ in Eq. (1), where $rp$ and $rg$ are the pinion and gear base radii, and $θp$ and $θg$ are the rotational deflections of the pinion and gear. In this case, the load is constant and $F(t)=F=Tp/rp$ represents the mesh force due to a constant applied torque $Tp$. Gear vibrations are lightly damped. Damping is reported to be about 2% in Ref. [5], so approximation with viscous damping $c$ is adopted. Assuming well-lubricated gear contact, the tooth friction forces are assumed small with no significant impact on the vibration. With this formulation, $x(t)=θprp+θgrg$ is the dynamic transmission error in Eq. (1), and $xs=g(t)$ is the unloaded static transmission error.

The force-deflection function $f(x,t)$ is the crucial quantity in this work. For a variety of systems, the force-deflection function can be obtained from experiments, finite element analysis, and analytical contact models. In gears, it is routinely calculated from computational models and readily measured. The force-deflection function, independent of the physical system in question and the method used to obtain it, can be represented by a Taylor series around $xm$ up to the $nth$ order byDisplay Formula

(3)$f(x,t)=∑i=0n1i!∂if(x,t)∂xi|x=xm(x-xm)i=∑i=0nBi(t)(x-xm)i$

where $xm$ is the deflection induced by the mean applied load. Averaging the force-deflection function and the applied force over their shortest common period $P$, $xm$ is found fromDisplay Formula

(4)$∫0PF(t)dt=∫0Pf(xm,t)dt$

The Fourier series expansion of $Bi(t)$ in Eq. (3) isDisplay Formula

(5)$Bi(t)=βi,0+∑r=1pβi,rcos(rξt-Φi,r), i=0,1,…,n$

where $ξ=2π/Pf$ is the frequency of $f(x,t)$. The Fourier expansion of $F(t)$ isDisplay Formula

(6)$F(t)=F0+∑r=1pFrcos(rζt-Rr)$

where $ζ=2π/PF$ is the frequency of $F(t)$. Substitution of $y=x-xm$, Eqs. (5), and (6) into Eqs. (1) and (2) givesDisplay Formula

(7)$my··+cy·+C(y,t)∑i=0n[βi,0+∑r=1pβi,rcos(rξt-Φi,r)]yi =F0+∑r=1pFrcos(rζt-Rr)$
Display Formula
(8)$C(y,t)={0if y<-xm+g(t)1if y≥-xm+g(t)$

Substitution of $β0,0=F0$, a result of Eqs. (4), (5), and (6) into Eq. (7) eliminates the mean load $F0$ under contact, givingDisplay Formula

(9)$my··+cy·+C(y,t)∑i=1n[βi,0+∑r=1pβi,rcos(rξt-Φi,r)]yi =∑r=1pFrcos(rζt-Rr)-C(y,t)∑r=1pβ0,rcos(rξt-Φ0,r) -F0[C(y,t)-1]$

so the mean load $F0$ acts only when contact is lost.

Fourier expansion of $g(t)$ givesDisplay Formula

(10)$g(t)=g0+∑r=1pgrcos(rξt-Ψr)$

The definition $d=xm-g0$ and substitution of $T=ωnt$, $ωn=β1,0/m$, and $y=ud$ give the nondimensional forms of Eqs. (8) and (9) asDisplay Formula

(11)$u"+2μu'+C(u,T)∑i=1n[αi,0+∑r=1pαi,rcos(rΩfT-φi,r)]ui =∑r=1pqrcos(rΩFT-ρr) -C(u,T)[α0,1cosΩfT+∑r=2pα0,rcos(rΩfT-φ0,r)] -α0,0[C(u,t)-1](11)$
Display Formula
(12)$C(u,T)={0if u<-1+∑r=1percos(rΩfT-ψr)1if u≥-1+∑r=1percos(rΩfT-ψr)$

where $()'=d/dT$, the nondimensional quantities areDisplay Formula

(13)$Ωf=ξωn, ΩF=ζωn, μ=c2mωn, qr=Frdβ1,0(r≠0), αi,r=βi,rdi-1β1,0, er=grd(13)$

and the phase angles relative to $Φ0,1$ are $(i=0,1,…,n$ and $r=0,1,…,p)$Display Formula

(14)$φi,r=Φi,r-Φ0,1, ρr=Rr-Φ0,1, ψr=Ψr-Φ0,1$

###### Nondimensional Parameters.

The physical interpretation of the important nondimensional parameters are listed below and summarized in Table 1.

1. (1)The mean load $α0,0$, linear mean stiffness $α1,0$, and nonlinear stiffnesses $αi,0$$(i>1)$ come from Taylor expansion of the force-deflection function by Eq. (3). In gears, $α0,0$ represents the constant mesh force, $α1,0$ the linear mesh stiffness, and $αi,0$ ($i>1$) the stiffness nonlinearities that include partial contact loss, nonlinear Hertz contact and others.
2. (2)The excitation harmonics from the force-deflection function, $α0,r$, the harmonics of linear stiffness, $α1,r$, the harmonics of stiffness nonlinearity $αi,r$ ($i>1$) $r>0$ and associated phase angles $φi,r$ come from Fourier expansion of the Taylor coefficients $Bi(t)$ by Eq. (5). They exist in gears due to periodic engagement/disengagement of the gear teeth. The quantity $α0,r$ excites the dynamics. It is close to what is called the “static transmission error excitation” for unmodified gears. The $α1,r$ represents the periodic change in the linear mesh stiffness. The periodicity of the stiffness nonlinearity, that is, changes in the strength of nonlinearity as the gear teeth engage/disengage, is contained in $αi,r$, $(i>1)$.
3. (3)The harmonics of the external excitation $qr$ and associated phase angles $ρr$ come from Fourier expansion of the applied force $F(t)$ by Eq. (6). When the applied torque is constant, $qr=0$, and the first harmonic of the fluctuating torque is reflected in $q1$ that, if exists, further excites the system.
4. (4)The harmonics of the unloaded deflection $er$ and associated phase angles $ψr$ come from Fourier expansion of the unloaded deflection $g(t)$ by Eq. (10) and represent the time-dependent unloaded transmission error.

###### Partial Contact Loss.

Partial contact loss occurs when portions of nominal gear contact lines lose contact while the other parts are still in contact [9,10]. This contrasts with total contact loss, where the gear mesh ceases to transmit any force. Gear vibrations combined with tooth profile and lead modifications cause partial contact loss. Partial contact loss depends heavily on applied torque. The mesh stiffness of modified gears, which depends on the total length of the contacting lines, changes with applied torque. Figure 1 demonstrates this dependence by comparing the dimensional mesh stiffness of unmodified and modified gears. The modified gears depend heavily on torque; the unmodified gears do not. This dependence of stiffness on torque is a source of nonlinearity, and it results from partial contact loss. The parameters $αi,r$ and $φi,r$ for $i>1$ incorporate partial contact loss, and any other nonlinearity present in the force-deflection function, into the equation of motion. Static finite element analysis at many points in a tooth mesh cycle is one effective way to compute the force-deflection function considering partial contact loss and Hertz contact nonlinearity. Commercial software can do this accurately.

###### Independence of Total Contact Loss From Applied Torque.

The analysis helps explain the observation in Refs. [6,22,30,37,38] that, for steady torque systems that do not involve gear rattle, increasing the applied torque does not reduce the occurrence of contact loss near primary resonance in unmodified gears; that is, if contact loss occurs at one torque, it will occur at any other torque. Once contact loss occurs, however, the extent of contact loss and the nature of the nonlinear response will depend on applied torque.

To explain, Fig. 2 shows dimensional, mean-removed excitations from an unmodified and a modified gear pair. The dimensional mean removed excitation $B0(t)-F0$ and the mean deflection $xm$ for unmodified gears are almost directly proportional to the applied torque. The modified gears show no such proportionality. The first harmonic of the nondimensional excitation $α0,1$ is the dominant driver of vibrations near primary resonance. This quantity from the unmodified and modified gear pairs is shown in Fig. 3. Changing the applied torque does not change the first harmonic of the excitation for the unmodified gear pair much, but it dramatically changes this quantity for the modified gear pair. In other words, the first harmonic of the excitation is sensitive to torque for modified gears, but insensitive to torque for unmodified gears. Because the first harmonic of the nondimensional excitation $α0,1$ from unmodified gears is similar for any applied torque, the applied torque has no effect upon the occurrence of contact loss. This mathematical explanation translates into a physical explanation as follows. In unmodified gears, a higher torque nearly linearly increases both: (a) the excitation and so the vibration amplitude, and (b) the threshold of vibration amplitude needed to trigger contact loss. In summary

1. (1)For unmodified gears, the occurrence of contact loss is independent of the applied torque.
2. (2)For modified gears, the occurrence of contact loss depends on the applied torque.
3. (3)Once contact loss starts to occur, the vibrations depend on applied torque. The amplitude of applied torque affects the duration of separation in a mesh period, and the vibration amplitude.

## Analysis Method

###### Dynamic Response Near Primary Resonance.

This section presents approximate analytical solutions of Eqs. (11) and (12) near the primary resonance region using the method of multiple scales. To unify the analysis, we consider that parametric and external excitations are both present and their frequencies are equal and close to the natural frequency, $Ω=Ωf=ΩF≈1$. The reason for setting $Ωf=ΩF$ is to keep the length of the closed-form solution to a minimum and to allow for periodic external torque excitation near resonant speed; thus there is no reason to consider $Ωf≠ΩF$. Following Ref. [39] to expand $u$, $T$, and the forcing frequency $Ω$ near the primary resonance givesDisplay Formula

(15)$u(T,ɛ)=u0(T0,T1,T2)+ɛu1(T0,T1,T2)+ɛ2u2(T0,T1,T2)Tn=ɛnT, Ω=1+σ=1+ɛσ∧$

The small parameter $ɛ$ is a detuning parameter such that the excitation frequency is close to the natural frequency. The separation function in Eq. (12) is rewritten as $C=1+H$ whereDisplay Formula

(16)$H(u,t)=12sgn[u+1-∑r=1percos(rT-ψr)]-12$

This allows $p$ contact losses per period. Fourier expansion of $H(u,t)$ is admissible because, like $u(t)$, it has period $2π$, givingDisplay Formula

(17)$H=Δ0+∑r=1NΔrcos(rT0-νr)$

The quantities $Δ0$ and $Δr$ depend on $u(t)$ and are not yet known. Arbitrarily, many harmonics can be used to represent the separation function $H$, but $N$ must not be less than the number of contact losses per period, $N≥p$. These steps turn contact loss into a form manageable by the method of multiple scales [22,30,40,41].

The $O(ɛ)$ parameters areDisplay Formula

(18)$α0,1=ɛα∧0,1, α2,0=ɛα∧2,0, α3,0=ɛα∧3,0, q1=ɛq∧1, μ=ɛμ∧,σ=ɛσ∧, Δ0=ɛΔ∧0(18)$

These orderings reflect the assumption that the nonlinearities, forcing, and damping are small. The mean value $Δ0$ of the separation function $H$ is ordered by assuming that the duration of contact loss is small compared to the period of vibration. Furthermore, all harmonics are assumed one order smaller than the corresponding mean values. This gives the $O(ɛ2)$ parameters asDisplay Formula

(19)$α2,1=ɛ2α∧2,1, Δ1=ɛ2Δ∧1, Δ2=ɛ2Δ∧2$

We consider up to the third order polynomial approximation of the nonlinear force-deflection function, that is, $αi,r≈0$ for $i≥4$, and the first harmonic of the periodically varying force-deflection function, that is, $αi,r≈0$ for $r≥2$. Parametric instability due to the first harmonic of the linear stiffness variation $α1,1$ is possible, but this needs to be treated separately. The third and higher harmonics of the separation function ($Δr$ for $r≥3$) and excitations $α0,r$ and $qr$ for $r≥2$ do not contribute to the solution, so they are excluded from subsequent equations for brevity.

Substitution of Eq. (15) into Eq. (11) and combining like orders of $ɛ$ gives the perturbation equationsDisplay Formula

(20)$D02u0+u0=0$
Display Formula
(21)$D02u1+u1=-2D0D1u0-α∧2,0u02-α∧3,0u03-2μ∧D0u0 -Δ∧0(u0+α0,0)-α∧0,12(eiT0+e-iT0) + q∧12[ei(T0-ρ1)+e-i(T0-ρ1)]$
Display Formula
(22)$D02u2+u2=-2μ∧(D1u0+D0u1)-2D0D1u1-2D0D2u0 -D12u0-(u1+α∧2,0u02+α∧3,0u03)Δ∧0-2α∧2,0u0u1 -3α∧3,0u02u1-α∧0,1Δ∧02(eiT0+e-iT0) -(α0,0+u0)Δ∧12[ei(T0-ν1)+e-i(T0-ν1)] -(α0,0+u0)Δ∧22[ei(2T0-ν2)+e-i(2T0-ν2)] -α∧2,1u022[ei(T0-φ2,1)+e-i(T0-φ2,1)](22)$

where $Di=∂/∂Ti$. The general solution of Eq. (20) isDisplay Formula

(23)$u0=A(T1,T2)eiT0+A¯(T1,T2)e-iT0$

Substitution of Eq. (23) into the first order Eq. (21) yields secular terms that must be eliminated for $u0$ to be periodic, giving the conditionDisplay Formula

(24)$2iD1A+A[2iμ∧+Δ∧0+3α∧3,0A2A¯]+α∧0,12eiσ∧T1-q∧12ei(σ∧T1-ρ1)=0$

Transforming to polar coordinates by setting $A=1/2aeiβ$, $β=σ∧T1-γ$, with $a,γ∈R$ gives the real and imaginary parts of Eq. (24) asDisplay Formula

(25)$aD1γ=aσ∧-12aΔ∧0-38α∧3,0a3-12α∧0,1cosγ+12q∧1cos(γ-ρ1)D1a=-μ∧a-α∧0,12sinγ+12q∧1sin(γ-ρ1)$

###### First Order Perturbation.

If the process is terminated at the first order perturbation, only Eq. (25) is considered. The steady state solutions of Eq. (23) result from $D1a=0$, $D1γ=0$ in Eq. (25). This gives, after use of Eq. (18), the frequency-amplitude relation cast in terms of the original parameters independent of $ɛ$ asDisplay Formula

(26)$aσ-12aΔ0-38α3,0a3-12α0,1cosγ+12q1cos(γ-ρ1)=0-μa-12α0,1sinγ+12q1sin(γ-ρ1)=0$

These two equations can be combined by eliminating $γ$, givingDisplay Formula

(27)$μ2a2+(aσ-12Δ0a-38α3,0a3)2-α0,12+q124+α0,1q12cosρ1=0$

Many parameters ($α0,0$, $α2,0$, $α2,1$, $φ2,1$, $Δ1$, $Δ2$, $ν1$, and $ν2$) do not appear in Eq. (26) and only contribute through the second order perturbation. After calculation of the general solution of Eq. (21), the approximate solution is constructed from $u=u0+ɛu1+O(ɛ2)$, givingDisplay Formula

(28)$u= acos(T0-γ)+α2,0a22[13cos(2T0-2γ)-1]+α3,0a332cos(3T0-3γ)-α0,0Δ0+O(ɛ2)$

###### Second Order Perturbation.

Proceeding to the second order perturbation, substitution of Eqs. (28) into (22) yields secular terms that need to be eliminated for $u1$ to have a periodic solution, giving the conditionDisplay Formula

(29)$D12A+2iD2A+2μD1A-103α∧2,02AA¯+38α∧3,02A3A¯2+Δ∧0(α∧0,12eiσ∧T1-2Aα0,0α∧2,0+3A2A¯α∧3,0)+α0,0Δ∧12ei(σT1-ν1)+A2α∧2,12e-i(σ∧T1-φ2,1)+A¯Δ∧22ei(2σ∧T1-ν2)+AA¯α∧2,1ei(σ∧T1-φ2,1)=0(29)$

Substitution of $D1A$ from Eq. (24) into Eq. (29) and use of the polar coordinates $a(T1,T2)$ and $γ(T1,T2)$ give the real and imaginary parts of Eq. (29) asDisplay Formula

(30)$aD2γ=[α∧0,1cosγ-q∧1cos(γ-ρ1)](332α∧3,0a2+14σ∧)-38Δ∧0α∧0,1cosγ+14μ∧α∧0,1sinγ -14μ∧q∧1sin(γ-ρ1)-12α0,0Δ∧1cos(ν1-γ)-14Δ∧2acos(ν2-2γ)-38α∧2,1a2cos(Φ2,1-γ) +18Δ02a-316α∧3,0Δ∧0a3+12μ∧2a+15256a5α∧3,02+aα∧2,0α0,0Δ∧0+512a3α∧2,02D2a=[α∧0,1sinγ-q∧1sin(γ-ρ1)](932α∧3,0a2+14σ∧)-38Δ∧0α∧0,1sinγ-14μ∧α∧0,1cosγ +14μ∧q∧1cos(γ-ρ1)+38μ∧α∧3,0a3+12α0,0Δ∧1sin(ν1-γ) +14Δ∧2asin(ν2-2γ)+18α∧2,1a2sin(Φ2,1-γ)$

Steady-state vibration requires $D1a=D1γ=0$ in Eq. (25) and $D2a=D2γ=0$ in Eq. (30). Combination of the real and imaginary parts of Eqs. (25) and (30) according to the reconstitution method [41,42], substitution of $D1a=D1γ=0$ into Eq. (25) and $D2a=D2γ=0$ into Eq. (30), and use of Eqs. (18) and (19) in that sequence give the frequency-amplitude relationsDisplay Formula

(31)$aσ-Δ02a-3α3,0a38-α0,12cosγ+q12cos(γ-ρ1)+{14μα0,1sinγ-14μq1sin(γ-ρ1)+[α0,1cosγ-q1cos(γ-ρ1)](332α3,0a2+14σ) -38Δ0α0,1cosγ-12α0,0Δ1cos(ν1-γ)-14Δ2acos(ν2-2γ)-38α2,1a2cos(Φ2,1-γ)+18Δ02a-316α3,0Δ0a3+12μ2a+15256a5α3,02 +aα2,0α0,0Δ0+512a3α2,02}=0-μa-12α0,1sinγ+12q1sin(γ-ρ1)+{-14μα0,1cosγ+14μq1cos(γ-ρ1) +[α0,1sinγ-q1sin(γ-ρ1)](932α3,0a2+14σ)-38Δ0α0,1sinγ+38μα3,0a3+12α0,0Δ1sin(ν1-γ)+14Δ2asin(ν2-2γ) +18α2,1a2sin(Φ2,1-γ)}=0(31)$

from which $a$ and $γ$ can be solved. The expressions outside of the brackets are the first order perturbation solution, and the expressions inside the brackets are the second order corrections.

###### Separation.

The quantities $Δr$ and $νr$ in Eq. (17) are found using the separation instants $Ti$. These separation instants are solved by substitution of $u$ from Eq. (28) into Eq. (16) for $H=-1$ (i.e., $C=0$). For arbitrarily many contact losses per period $p$, the solutions $Ti$, $i=1,…,2p+1$, mark the instants when contact is lost or re-established. Let $Hi$ denote whether there is contact ($Hi=0$) or not ($Hi=-1$) when $Ti≤T. The Fourier coefficients and phase angles of the separation function $H$ in Eq. (17) are thenDisplay Formula

(32)$Δ0=12π∑i=12p(Ti+1-Ti)Hi$
Display Formula
(33)$Δr=Δr,a2+Δr,b2, νr=tan-1(Δr,b,Δr,a),r=1,…,NΔr,a=1rπ∑i=12p(sinrTi+1-sinrTi)Hi,Δr,b=-1rπ∑i=12p(cosrTi+1-cosrTi)Hi, r=1,…,N$

The perturbation solution considers a maximum of $N=2$ because harmonics of $H$ higher than the second, i.e., $Δr$ for $r>2$ in Eq. (17), do not yield secular terms. This limits contact loss to twice-per-period because $N≥p$. Even when the contact loss is limited to twice-per-period, finding the separation instants $Ti$ requires solving Eq. (16) for $p=2$. Finding a closed-form solution is not feasible unless: (1) contact loss occurs once per period ($p=1$), and (2) the mean change in solution $u$ due to contact loss does not affect when contact loss starts and ends ($α0,0Δ0=0$). Adopting these assumptions gives the closed-form expressionsDisplay Formula

(34)$Δ0=-1+τπ, Δr=2rπsinrτνr=rtan-1(-asinγ+e1sinψ1,acosγ-e1cosψ1),r=1,2$
Display Formula
(35)$τ=cos-1(-1+α2,0a/2(acosγ-e1cosψ1)2+(asinγ-e1sinψ1)2)$

The numerical solution for $Δ0$, $Δr$, and $νr$ using Eqs. (32) and (33) is computationally demanding. To save computation when analyzing a system with the second order perturbation solution Eq. (31), we use a coarse frequency resolution across the frequency range of interest to confirm that contact loss occurs once per period using numerical solutions of Eqs. (32) and (33). Once confirmed, the analytical expressions in Eqs. (34) and (35) are used. The analytical expressions in Eqs. (34) and (35) and numerical solution of Eqs. (32) and (33) give almost identical results for the cases analyzed in this work, so the presented results use the analytical expressions. This comparison justifies the assumption of only one instance of contact loss per mesh period for the current system. For a set of parameters with even stronger nonlinearity at resonance, the perturbation solution assumptions must be re-evaluated to confirm their suitability.

The vibration amplitude and frequency at the onset of total contact loss can be found from the foregoing analytical solution. Just before the onset of total contact loss, the mean value of the separation function and its harmonics are zero and $Δ0=Δr=0$ in Eq. (34). This gives $τ=π$. Substitution of $τ=π$ into Eq. (35) gives the amplitude at the onset of total contact loss asDisplay Formula

(36)$a=4e1cos(γ-ψ1)-2α2,0±2e12α2,02+4-4e12sin2(γ-ψ1)-4e1α2,0cos(γ-ψ1)4-α2,02$

Substitution of $a$ from Eq. (36) and $Δ0=Δr=0$ into Eq. (31) gives the frequency $σ$ and the phase $γ$ at the onset of total contact loss. The closed-form solution is lengthy using the second order perturbation solution, but a compact expression can be obtained using the first order perturbation solution. The first order does not consider quadratic nonlinearity and unloaded fluctuations (they do not yield secular terms), so $α2,0=e1=0$. Subsequently, Eq. (36) reveals that $a=1$ initiates total contact loss. Substitution of $a=1$ and $Δ0=0$ into Eq. (27) gives the two frequencies where total contact loss starts asDisplay Formula

(37)$σ=±α0,12+q124-μ2-α0,1q12cosρ1+38α3,0$

###### Key Points from the Analytical Solution

1. (1)The first order perturbation gives the frequency amplitude relation by Eq. (27). It considers only the cubic nonlinearity $α3,0$ and the mean reduction of stiffness due to total contact loss $Δ0$.
2. (2)The second order perturbation solution gives the frequency amplitude relation by Eq. (31). It considers quadratic nonlinearities ($α2,0$, $α2,1$, $φ2,1$), and cubic nonlinearity ($α3,0$). The mean value of the separation ($Δ0$) and its first two harmonics ($Δ1$, $Δ2$) with associated phase angles ($ν1$, $ν2$) contribute. These quantities account for total contact loss.
3. (3)The mean value ($Δ0$) and the first two harmonics ($Δ1$, $Δ2$) of the separation function and their phase angles ($ν1$, $ν2$) are found from the perturbation equations numerically from Eqs. (32) and (33) or analytically from Eqs. (34), (35) when contact loss is limited to once per period ($p=1$). Both harmonics are crucial for an accurate response as the upcoming results will show.
4. (4)When the expression inside $cos-1$ in Eq. (35) is not in the range $(-1,1)$, total contact loss does not occur. The amplitude and frequency at onset of contact loss is analytically given by Eq. (37) from the first order perturbation. The second order approximation of the amplitude at onset of contact loss is found from Eq. (36) and the frequencies can be solved by substitution of that amplitude into Eq. (31).

###### Stability.

The solution of Eq. (31) can give one or three solutions for a given frequency. When more than one solution is possible, only the stable solutions would be observed in practice or calculated by numerical integration. The stability of a solution is found by linearizing Eqs. (25) and (30) about a stationary point $(a,γ)$. This yields the matrix formDisplay Formula

(38)$ddT(aγ)=J(aγ)$

where the matrix $J$ denotes the Jacobian. The solution is unstable if any eigenvalues of $J$ have a positive real part. The algebraic expression for $J$ is prohibitively long and provides no practical insight. For these reasons, the stability is found numerically from the Jacobian matrix $J$ at the stationary points.

## Results

In this section, dynamic response of the equation of motion Eq. (11) obtained using the perturbation solution in Eq. (28) is compared with gear vibration experiments from the literature, numerical integration, and the linearized solution. The first order perturbation solution is found from Eq. (27), and the second order perturbation solution is found from Eq. (31).

###### Key Regions in Dynamic Response and Comparison With Numerical Solution.

Independent of the physical system analyzed, there are four possible distinct regions in the nonlinear response. To investigate the characteristics of these four regions, time domain numerical integration and the second order perturbation solution of an example system are shown in Fig. 4. The linear time-invariant solution (only linear stiffness and external excitation exist, $α1,0,α0,1≠0$) is plotted to provide a comparison:

###### Region 1.

The linear and nonlinear solutions give almost identical responses. The response is relatively small. There is no contact loss, and nonlinear terms negligibly affect the response.

###### Region 2.

The nonlinear solution deviates significantly from the linear solution as the vibration amplitude grows. This region does not exist, i.e., becomes identical to region 1, if $α2,0=0$ and $α3,0=0$. Although the contact is maintained, this region represents partial contact loss in physical systems, where the nonlinearities are due to a reduction in instantaneous dynamic stiffness. The nonlinear terms $α2,0$, $α3,0$, and the time-dependence $α2,1$, $φ2,1$ dominate the response.

###### Region 3.

Total contact loss, manifest from the kink in the response at $Ω=0.97$, starts in this region. Here $C=0$. The vibration amplitude $a$ and the frequency of excitation $Ω$ at which total contact loss starts is given by Eq. (36). It depends on the quadratic nonlinearity, the fluctuation in the unloaded deflection $e1$, and the relative phase angle $γ-ψ1$.

###### Region 4.

Total contact loss occurs. The solution is dominated by the mean reduction in stiffness $Δ0$ due to separation, the first two harmonics of the separation function $Δ1$, $Δ2$, their phase angles $ν1$, $ν2$, and the mean load $α0,0$. The agreement between the numerical and perturbation solution is excellent, justifying the assumptions in the perturbation solution. The nonlinearities are stronger if the periods of separation become long, which would degrade the agreement somewhat as shown in Ref. [37] for a simpler tooth mesh contact model in planetary gears.

###### Computational Results.

Computational dynamic response of a helical gear pair is obtained using the contact algorithm in Refs. [9,10]. That model has been successfully compared against experiments and a specialized finite element/contact mechanics model of gear dynamics. Some elastic behaviors, such as corner contact, radius of curvature effects, buttressing effects, and the coupling between the deflections of adjacent teeth, are neglected in this contact algorithm.

The gear parameters, which are from [9], are listed in Table 2. The nondimensional parameters are given in Table 3.

The dynamic response of the helical gear pair with profile and lead modifications and the force-deflection function shown in Fig. 5 is plotted in Fig. 6. The computational dynamic response and perturbation solution agree throughout most of the operating frequencies. There is a clear difference between the linear and the nonlinear response. Region 2 of Fig. 4 identifies this type of response. Physically, this difference is attributed to partial contact loss, where some portions of the nominal contact lines separate while the gear mesh as a whole is still engaged. The perturbation solution predicts the nonlinear dynamics arising from partial contact loss.

The investigation of practical gear vibration raises the question of whether all the nonlinear terms ($α2,0$, $α3,0$, $α2,1$, $φ2,1$) are necessary to obtain an accurate response or not. The importance of these terms are assessed by investigating the perturbation solution when these are set to zero one at a time, as shown in Fig. 7. The response becomes inaccurate when any one of them is set to zero; all nonlinear terms are necessary to obtain an accurate response.

###### Experimental Results.

In this section, experimental measurements of rotational vibration of unity ratio gear pairs with profile and lead modifications [5] are compared with the second order perturbation solution. The experimental measurements are available for various tip relief starting roll angles and for different applied torques.

The gear parameters from Ref. [5] are listed in Table 4. The nondimensional parameters are given in Table 3. The gears in this test rig are sufficiently isolated from the vibrations of the supporting structure to justify use of a single degree-of-freedom oscillator to model the dynamics [24].

The force-deflection function would ideally come from the experimental setup by slowly rotating the gears through a mesh period and measuring transmission error at various applied torques. Because these data are not given, the force-deflection function is obtained from finite element analysis [43] of the gears. When the tip relief starts at roll angle 20.9 deg, the force-deflection function calculated by finite element analysis is shown in Fig. 8 for various torques. Even if the finite element analysis could precisely replicate the elastic behavior of the experimental setup, the crucial gear tooth profile and lead modifications are specified within a manufacturing tolerance ($μm$). These errors can cause a mismatch between the experimental data and perturbation solution.

The experimental measurements and the second order perturbation solution for three gear pairs are compared in Fig. 9 (where the applied torque is 85 Nm and the 10$μm$ profile modification starts at 20.9 deg), Fig. 10 (where the applied torque is 170 Nm and the 10$μm$ profile modification starts at 20.9 deg), and Fig. 11 (where the applied torque is 170 Nm and the profile 10$μm$ modification starts at 23.6 deg). Table 3 lists the nondimensional parameters used in this analysis. All gears have 5$μm$ lead crown modification.

The perturbation solutions agree accurately with the experimental results in all three cases. Damping is estimated to be between 2% and 5% as shown in Table 3. This agrees with the reported damping of 2% in Ref. [5] and yields perturbation results that closely match the experiments, as shown in the figures. The agreement with the experimental measurements justifies the use of the single-degree-of-freedom rotational oscillator model and the exclusion of friction forces for gear vibrations. These justifications, however, are specific to the experimental setup because compliant bearings or poorly lubricated gears with high friction may require different physical assumptions.

The onset of total contact loss is manifest by the kinks in the perturbation solution curves. Even in regions without total contact loss, the linear response differs from the experimental measurements and perturbation solutions. Region 2 of Fig. 4 characterizes this type of response. The difference is attributed to partial contact loss and occurs due to the profile and lead modifications. The good agreement provides experimental evidence that the perturbation solution predicts the nonlinear dynamic response due to partial contact loss.

## Conclusions

This work derives approximate, closed-form analytical solutions for the nonlinear vibration of gear pairs with profile and lead modifications described by a specified nonlinear, periodic force-deflection curve. The solution includes reductions in profile and face contact ratios as a result of vibration combined with modifications. Misalignments in the gear assembly are considered by their effects on the force-deflection curve. The solution includes nonlinear effects of partial contact loss, that is, reductions in profile and face contact ratios as a result of vibration combined with tooth modifications. The method of multiple scales provides the analytical perturbation solutions in the primary resonance region.

1. (1)The dynamic excitation in gear pairs is mathematically shown to come from the periodicity of the force-deflection function. This periodicity is due to the periodic engagement/disengagement of the gear teeth resulting in linear and quadratic mesh stiffness fluctuations that drive the gear vibration. The occurrence of total contact loss in unmodified gears is shown to be independent of applied torque because the vibration amplitude increases with larger torques, so does the threshold of vibrations needed to trigger total contact loss. This is not true for modified gears; the occurrence of contact loss depends on applied torque in modified gears. Once contact loss starts to occur, however, the amplitude of nonlinear vibration does depend on applied torque.
2. (2)The method of multiple scales gives the approximate, closed-form analytical solutions, which provide physical insight on the effect of design parameters on nonlinear vibration. Comparisons with gear vibration experiments and numerical integration verify the analytical solution by perturbation. It is observed that the second order perturbation solution is significantly more accurate than the first order perturbation solution, indicating the nonlinearity is strong in the physical experiments.
3. (3)Key regions in the nonlinear response are:
1. (a)Linear region: vibration amplitude is relatively low, and the system behaves linearly.
2. (b)Partial contact loss region: The mean stiffness drops due to partial contact loss from reduced profile and face contact ratios as a result of vibration; nevertheless, contact is maintained. During partial contact loss, nominal contact lines lose contact although the gear teeth remain engaged. Partial contact loss is most prominent in systems with tooth profile and lead modifications. Quadratic and cubic nonlinearities capture the partial contact loss nonlinearity in modified gears; quadratic and cubic nonlinearities drive the nonlinear response.
3. (c)Onset of total contact loss region: Vibrations reach the threshold to cause total contact loss. The onset frequency and amplitude is independent of applied torque in unmodified gears. The amplitude and the frequency at the onset of total contact loss is analytically given and is consistent with experiments.
4. (d)Total contact loss region: The contact separates fully and the mean load brings the system back into contact. Total contact loss in gears occurs when the gear teeth disengage completely. When contact loss occurs, vibration amplitude depends on applied torque. Contact loss is captured by the mean reduction in mesh stiffness and its two harmonics. For contact loss, the perturbation solution considers one separation per period, and excludes further period doubling and chaos. The validity of this assumption is borne out by the correlation with experimentally measured vibration.

## References

Chang, S., Houser, D. R., and Harianto, J., 2003, “Tooth Flank Corrections of Wide Face Width Helical Gears That Account for Shaft Deflections,” Proceedings of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering, Chicago, IL, September 2–6, ASME Paper No. DETC2003/PTG-48072.
Litvin, F. L., Fuentes, A., Gonzalez-Perez, I., Carnevali, L., and Kawasaki, K., 2003, “Modified Involute Helical Gears: Computerized Design, Simulation of Meshing, and Stress Analysis,” Tech. Report No. ARL-CR-514.
Umezawa, K., Sato, T., and Ishikawa, J., 1984, “Simulation of Rotational Vibration of Spur Gears,” Bull. JSME, 27(223), pp. 102–109.
Velex, P., and Maatar, M., 1996, “A Mathematical Model for Analyzing the Influence of Shape Deviations and Mounting Errors on Gear Dynamic Behavior,” J. Sound Vib., 191(5), pp. 629–660.
Kahraman, A., and Blankenship, G. W., 1999, “Effect of Involute Tip Relief on Dynamic Response of Spur Gear Pairs,” J. Mech. Des., 121(5), pp. 313–315.
Parker, R. G., Vijayakar, S. M., and Imajo, T., 2000, “Non-Linear Dynamic Response of a Spur Gear Pair: Modelling and Experimental Comparisons,” J. Sound Vib., 237(3), pp. 435–455.
Ambarisha, V. K., and Parker, R. G., 2007, “Nonlinear Dynamics of Planetary Gears Using Analytical and Finite Element Models,” J. Sound Vib., 302(3), pp. 577–595.
Eritenel, T., and Parker, R. G., 2011, “Nonlinear Vibration of Gear Pairs With Tooth Surface Modifications at Primary Resonance Using a Perturbation Method,” Proceedings of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering, Washington, DC, August 28–31, ASME Paper No. DETC2011/MSNDC-48689.
Eritenel, T., and Parker, R. G., 2012, “Three-Dimensional Nonlinear Vibration of Gear Pairs,” J. Sound Vib., 331, pp. 3628–3648.
Eritenel, T., and Parker, R. G., 2012, “An Investigation of Tooth Mesh Nonlinearity and Partial Contact Loss in Gear Pairs Using a Lumped-Parameter Model,” Mech. Mach. Theory, 56, pp. 28–51.
Gregory, R. W., Harris, S. L., and Munro, R. G., 1963, “Dynamic Behavior of Spur Gears,” Proc. Inst. Mech. Eng., 178(8), pp. 261–266.
Andersson, A., and Vedmar, L., 2003, “A Dynamic Model to Determine Vibrations in Involute Helical Gears,” J. Sound Vib., 260(2), pp. 195–212.
Baud, S., and Velex, P., 2002, “Static and Dynamic Tooth Loading in Spur and Helical Geared Systems-Experiments and Code Validation,” ASME J. Mech. Des., 124(2), pp. 334–346.
Kubur, M., Kahraman, A., Zini, D. M., and Kienzle, K., 2004, “Dynamic Analysis of a Multi-Shaft Helical Gear Transmission by Finite Elements: Model and Experiment,” ASME J. Vib. Acoust., 126(3), pp. 398–406.
Ajmi, M., and Velex, P., 2005, “A Model for Simulating the Quasi-Static and Dynamic Behaviour of Solid Wide-Faced Spur and Helical Gears,” Mech. Mach. Theory, 40(2), pp. 173–190.
Harris, S. L., 1958, “Dynamic Loads on the Teeth of Spur Gears,” IMechE Appl. Mech., 172, pp. 87–100.
Singh, R., Houser, D. R., and Kahraman, A., 1990, “Non-Linear Dynamic Analysis of Geared Systems,” Tech. Report No. NASA-CR-180495.
Özgüven, H. N., and Houser, D. R., 1988, “Mathematical-Models Used in Gear Dynamics—A Review,” J. Sound Vib., 121(3), pp. 383–411.
Özgüven, H. N., and Houser, D. R., 1988, “Dynamic Analysis of High Speed Gears by Using Loaded Static Transmission Error,” J. Sound Vib., 125(1), pp. 71–83.
Lee, C., Lin, H.-H., Oswald, F. B., and Townsend, D. P., 1991, “Influence of Linear Profile Modification and Loading Conditions on the Dynamic Tooth Load and Stress of High-Contact-Ratio Spur Gears,” ASME J. Mech. Des., 113(4), pp. 473–480.
Velex, P., and Ajmi, M., 2006, “On the Modelling of Excitations in Geared Systems by Transmission Errors,” J. Sound Vib., 290(3–5), pp. 882–909.
Liu, G., and Parker, R. G., 2008, “Dynamic Modeling and Analysis of Tooth Profile Modification for Multimesh Gear Vibration,” ASME J. Mech. Des., 130(2), p. 121402.
Kahraman, A., and Singh, R., 1990, “Non-Linear Dynamics of a Spur Gear Pair,” J. Sound Vib., 142(1), pp. 49–75.
Blankenship, G. W., and Kahraman, A., 1995, “Steady State Force Response of a Mechanical Oscillator With Combined Parametric Excitation and Clearance Type Non-Linearity,” J. Sound Vib., 185(5), pp. 743–765.
Padmanabhan, C., and Singh, R., 1996, “Analysis of Periodically Forced Nonlinear Hill's Oscillator With Application to a Geared System,” J. Acoust. Soc. Am., 99(1), pp. 324–334.
Kahraman, A., and Blankenship, G. W., 1996, “Interactions Between Commensurate Parametric and Forcing Excitations in a System With Clearance,” J. Sound Vib., 194(3), pp. 317–336.
Kahraman, A., and Singh, R., 1991, “Interactions Between Time-Varying Mesh Stiffness and Clearance Non-Linearities in a Geared System,” J. Sound Vib., 146(1), pp. 135–156.
Theodossiades, S., and Natsiavas, S., 2000, “Non-Linear Dynamics of Gear-Pair Systems With Periodic Stiffness and Backlash,” J. Sound Vib., 229(2), pp. 287–310.
Natsiavas, S., Theodossiades, S., and Goudas, I., 2000, “Dynamic Analysis of Piecewise Linear Oscillators With Time Periodic Coefficients,” Int. J. Nonlinear Mech., 35(1), pp. 53–68.
Liu, G., and Parker, R. G., 2008, “Nonlinear Dynamics of Idler Gear Systems,” Nonlinear Dyn., 53(4), pp. 345–367.
Kubo, A., 1978, “Stress Condition, Vibration Exciting Force and Contact Pattern of Helical Gears With Manufacturing and Alignment Error,” ASME J. Mech. Des., 100, pp. 77–84.
Matsumura, S., Umezawa, K., and Houjoh, H., 1996, “Rotational Vibration of a Helical Gear Pair Having Tooth Surface Deviation During Transmission of Light Load,” JSME Int. J., Ser. C, 39(3), pp. 614–620.
Velex, P., Maatar, M., and Raclot, J. P., 1997, “Some Numerical Methods for the Simulation of Geared Transmission Dynamic Behavior Formulation and Assessment,” ASME J. Mech. Des., 119(2), pp. 292–298.
Raclot, J. P., and Velex, P., 1999, “Simulation of the Dynamic Behaviour of Single and Multi-Stage Geared Systems With Shape Deviations and Mounting Errors by Using a Spectral Method,” J. Sound Vib., 220(5), pp. 861–903.
Kubo, A., Kiyono, S., and Makoto, F., 1986, “On Analysis and Prediction of Machine Vibration Caused by Gear Meshing: 1st Report, Nature of Gear Vibration and the Total Vibrational Excitation,” Bull. JSME, 29(258), pp. 4424–4429.
Kahraman, A., and Blankenship, G. W., 1999, “Effect of Involute Contact Ratio on Spur Gear Dynamics,” ASME J. Mech. Des., 121(1), pp. 112–118.
Bahk, C.-J., and Parker, R. G., 2011, “Analytical Solution for the Nonlinear Dynamics of Planetary Gears,” J. Comput. Nonlinear Dyn., 6(2), p. 021007.
Blankenship, G. W., and Kahraman, A., 1996, “Gear Dynamics Experiments Part 1: Characterization of Forced Response,” Proceedings of the ASME Design Engineering Technical Conferences, Irvine, CA, August 18–22.
Nayfeh, A. H., and Mook, D. T., 1979, Nonlinear Oscillations, Wiley, New York.
Nayfeh, A. H., 1993, Problems in Perturbation, John Wiley & Sons, New York.
Zhu, F., and Parker, R. G., 2006, “Perturbation Analysis of a Clearance-Type Nonlinear System,” J. Sound Vib., 292(3–5), pp. 969–979.
Luongo, A., and Paolone, A., 1999, “On the Reconstitution Problem in the Multiple Time-Scale Method,” Nonlinear Dyn., 19(2), pp. 133–156.
Vijayakar, S. M., 1991, “A Combined Surface Integral and Finite-Element Solution for a Three-Dimensional Contact Problem,” Int. J. Numer. Methods Eng., 31(3), pp. 525–545.
View article in PDF format.

## References

Chang, S., Houser, D. R., and Harianto, J., 2003, “Tooth Flank Corrections of Wide Face Width Helical Gears That Account for Shaft Deflections,” Proceedings of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering, Chicago, IL, September 2–6, ASME Paper No. DETC2003/PTG-48072.
Litvin, F. L., Fuentes, A., Gonzalez-Perez, I., Carnevali, L., and Kawasaki, K., 2003, “Modified Involute Helical Gears: Computerized Design, Simulation of Meshing, and Stress Analysis,” Tech. Report No. ARL-CR-514.
Umezawa, K., Sato, T., and Ishikawa, J., 1984, “Simulation of Rotational Vibration of Spur Gears,” Bull. JSME, 27(223), pp. 102–109.
Velex, P., and Maatar, M., 1996, “A Mathematical Model for Analyzing the Influence of Shape Deviations and Mounting Errors on Gear Dynamic Behavior,” J. Sound Vib., 191(5), pp. 629–660.
Kahraman, A., and Blankenship, G. W., 1999, “Effect of Involute Tip Relief on Dynamic Response of Spur Gear Pairs,” J. Mech. Des., 121(5), pp. 313–315.
Parker, R. G., Vijayakar, S. M., and Imajo, T., 2000, “Non-Linear Dynamic Response of a Spur Gear Pair: Modelling and Experimental Comparisons,” J. Sound Vib., 237(3), pp. 435–455.
Ambarisha, V. K., and Parker, R. G., 2007, “Nonlinear Dynamics of Planetary Gears Using Analytical and Finite Element Models,” J. Sound Vib., 302(3), pp. 577–595.
Eritenel, T., and Parker, R. G., 2011, “Nonlinear Vibration of Gear Pairs With Tooth Surface Modifications at Primary Resonance Using a Perturbation Method,” Proceedings of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering, Washington, DC, August 28–31, ASME Paper No. DETC2011/MSNDC-48689.
Eritenel, T., and Parker, R. G., 2012, “Three-Dimensional Nonlinear Vibration of Gear Pairs,” J. Sound Vib., 331, pp. 3628–3648.
Eritenel, T., and Parker, R. G., 2012, “An Investigation of Tooth Mesh Nonlinearity and Partial Contact Loss in Gear Pairs Using a Lumped-Parameter Model,” Mech. Mach. Theory, 56, pp. 28–51.
Gregory, R. W., Harris, S. L., and Munro, R. G., 1963, “Dynamic Behavior of Spur Gears,” Proc. Inst. Mech. Eng., 178(8), pp. 261–266.
Andersson, A., and Vedmar, L., 2003, “A Dynamic Model to Determine Vibrations in Involute Helical Gears,” J. Sound Vib., 260(2), pp. 195–212.
Baud, S., and Velex, P., 2002, “Static and Dynamic Tooth Loading in Spur and Helical Geared Systems-Experiments and Code Validation,” ASME J. Mech. Des., 124(2), pp. 334–346.
Kubur, M., Kahraman, A., Zini, D. M., and Kienzle, K., 2004, “Dynamic Analysis of a Multi-Shaft Helical Gear Transmission by Finite Elements: Model and Experiment,” ASME J. Vib. Acoust., 126(3), pp. 398–406.
Ajmi, M., and Velex, P., 2005, “A Model for Simulating the Quasi-Static and Dynamic Behaviour of Solid Wide-Faced Spur and Helical Gears,” Mech. Mach. Theory, 40(2), pp. 173–190.
Harris, S. L., 1958, “Dynamic Loads on the Teeth of Spur Gears,” IMechE Appl. Mech., 172, pp. 87–100.
Singh, R., Houser, D. R., and Kahraman, A., 1990, “Non-Linear Dynamic Analysis of Geared Systems,” Tech. Report No. NASA-CR-180495.
Özgüven, H. N., and Houser, D. R., 1988, “Mathematical-Models Used in Gear Dynamics—A Review,” J. Sound Vib., 121(3), pp. 383–411.
Özgüven, H. N., and Houser, D. R., 1988, “Dynamic Analysis of High Speed Gears by Using Loaded Static Transmission Error,” J. Sound Vib., 125(1), pp. 71–83.
Lee, C., Lin, H.-H., Oswald, F. B., and Townsend, D. P., 1991, “Influence of Linear Profile Modification and Loading Conditions on the Dynamic Tooth Load and Stress of High-Contact-Ratio Spur Gears,” ASME J. Mech. Des., 113(4), pp. 473–480.
Velex, P., and Ajmi, M., 2006, “On the Modelling of Excitations in Geared Systems by Transmission Errors,” J. Sound Vib., 290(3–5), pp. 882–909.
Liu, G., and Parker, R. G., 2008, “Dynamic Modeling and Analysis of Tooth Profile Modification for Multimesh Gear Vibration,” ASME J. Mech. Des., 130(2), p. 121402.
Kahraman, A., and Singh, R., 1990, “Non-Linear Dynamics of a Spur Gear Pair,” J. Sound Vib., 142(1), pp. 49–75.
Blankenship, G. W., and Kahraman, A., 1995, “Steady State Force Response of a Mechanical Oscillator With Combined Parametric Excitation and Clearance Type Non-Linearity,” J. Sound Vib., 185(5), pp. 743–765.
Padmanabhan, C., and Singh, R., 1996, “Analysis of Periodically Forced Nonlinear Hill's Oscillator With Application to a Geared System,” J. Acoust. Soc. Am., 99(1), pp. 324–334.
Kahraman, A., and Blankenship, G. W., 1996, “Interactions Between Commensurate Parametric and Forcing Excitations in a System With Clearance,” J. Sound Vib., 194(3), pp. 317–336.
Kahraman, A., and Singh, R., 1991, “Interactions Between Time-Varying Mesh Stiffness and Clearance Non-Linearities in a Geared System,” J. Sound Vib., 146(1), pp. 135–156.
Theodossiades, S., and Natsiavas, S., 2000, “Non-Linear Dynamics of Gear-Pair Systems With Periodic Stiffness and Backlash,” J. Sound Vib., 229(2), pp. 287–310.
Natsiavas, S., Theodossiades, S., and Goudas, I., 2000, “Dynamic Analysis of Piecewise Linear Oscillators With Time Periodic Coefficients,” Int. J. Nonlinear Mech., 35(1), pp. 53–68.
Liu, G., and Parker, R. G., 2008, “Nonlinear Dynamics of Idler Gear Systems,” Nonlinear Dyn., 53(4), pp. 345–367.
Kubo, A., 1978, “Stress Condition, Vibration Exciting Force and Contact Pattern of Helical Gears With Manufacturing and Alignment Error,” ASME J. Mech. Des., 100, pp. 77–84.
Matsumura, S., Umezawa, K., and Houjoh, H., 1996, “Rotational Vibration of a Helical Gear Pair Having Tooth Surface Deviation During Transmission of Light Load,” JSME Int. J., Ser. C, 39(3), pp. 614–620.
Velex, P., Maatar, M., and Raclot, J. P., 1997, “Some Numerical Methods for the Simulation of Geared Transmission Dynamic Behavior Formulation and Assessment,” ASME J. Mech. Des., 119(2), pp. 292–298.
Raclot, J. P., and Velex, P., 1999, “Simulation of the Dynamic Behaviour of Single and Multi-Stage Geared Systems With Shape Deviations and Mounting Errors by Using a Spectral Method,” J. Sound Vib., 220(5), pp. 861–903.
Kubo, A., Kiyono, S., and Makoto, F., 1986, “On Analysis and Prediction of Machine Vibration Caused by Gear Meshing: 1st Report, Nature of Gear Vibration and the Total Vibrational Excitation,” Bull. JSME, 29(258), pp. 4424–4429.
Kahraman, A., and Blankenship, G. W., 1999, “Effect of Involute Contact Ratio on Spur Gear Dynamics,” ASME J. Mech. Des., 121(1), pp. 112–118.
Bahk, C.-J., and Parker, R. G., 2011, “Analytical Solution for the Nonlinear Dynamics of Planetary Gears,” J. Comput. Nonlinear Dyn., 6(2), p. 021007.
Blankenship, G. W., and Kahraman, A., 1996, “Gear Dynamics Experiments Part 1: Characterization of Forced Response,” Proceedings of the ASME Design Engineering Technical Conferences, Irvine, CA, August 18–22.
Nayfeh, A. H., and Mook, D. T., 1979, Nonlinear Oscillations, Wiley, New York.
Nayfeh, A. H., 1993, Problems in Perturbation, John Wiley & Sons, New York.
Zhu, F., and Parker, R. G., 2006, “Perturbation Analysis of a Clearance-Type Nonlinear System,” J. Sound Vib., 292(3–5), pp. 969–979.
Luongo, A., and Paolone, A., 1999, “On the Reconstitution Problem in the Multiple Time-Scale Method,” Nonlinear Dyn., 19(2), pp. 133–156.
Vijayakar, S. M., 1991, “A Combined Surface Integral and Finite-Element Solution for a Three-Dimensional Contact Problem,” Int. J. Numer. Methods Eng., 31(3), pp. 525–545.

## Figures

Fig. 1

Dimensional gear mesh stiffness at 50 Nm (dashed-dotted line), 150 Nm (dashed line), and 250 Nm (solid line) using finite element analysis. (a) Unmodified spur gear pair from [36] (ICR = 1.37). (b) Modified spur gear pair from [5] (tip relief starts at 22.2 deg).

Fig. 2

Dimensional mean-removed excitation levels from the gears in Fig. 1 at 50 Nm (dashed-dotted line), 150 Nm (dashed line), and 250 Nm (solid line) using finite element analysis. (a) Unmodified spur gear pair. (b) Modified spur gear pair.

Fig. 3

The first harmonic of nondimensional excitation from the unmodified and modified gears with varying applied torque. The gears are the same as in Fig. 1.

Fig. 4

Dynamic response of equation of motion (Eqs. (11)–(12)) using second order perturbation and numerical integration. Parameters are: α0,1 = 2/75, α0,0 = 2/3, α2,0 = 1/4, α3,0 = -1/24, e1 = 0.1, ψ1 = 0, α2,1 = 0, μ = 0.0082. Important regions of solutions are marked with (1) through (4). (Stable perturbation solution: solid line, unstable perturbation solution: dashed line, numerical integration: circles, linear solution: dashed-dotted line.)

Fig. 5

Force-deflection function f(x,t) of a helical gear pair obtained by the method of Ref. [9]. Tooth surface modifications are: quadratic profile tip and root crown of 10μm and lead crown of 10μm on the gear.

Fig. 6

Primary resonance of a helical gear pair at 200 Nm using second order perturbation, numerical integration, and the linearized model. The numerical integration results are obtained using the method of Ref. [9] (stable perturbation solution: solid line, computational data by numerical integration: circles, linear solution: dashed-dotted line).

Fig. 7

Perturbation solution when α2,0, α3,0, α2,1, and φ2,1 are set to zero one at a time are shown. The nonzero parameters are the same as in Fig. 6. The circles denote numerical integration with no nonzero parameters.

Fig. 8

Force-deflection function f(x,t) obtained by finite element analysis of a spur gear pair with increasing applied torque. Gear data from Ref. [5], linear tip relief of 10μm starts at 20.9 deg.

Fig. 9

Primary resonance of a spur gear pair (tip relief start at 20.9 deg) at 85 Nm using second order perturbation and linearized model. Experimental data are from Fig. 3(c) of Ref. [5] (stable second order perturbation: solid line, unstable second order perturbation solution: dashed line, experimental measurement: circles, linear solution: dashed-dotted line).

Fig. 10

Primary resonance of a spur gear pair (tip relief start at 20.9 deg) at 170 Nm using second order perturbation and linearized model. Experimental data are from Fig. 3(b) of Ref. [5] (stable second order perturbation: solid line, unstable second order perturbation solution: dashed line, experimental measurement: circles, linear solution: dashed-dotted line).

Fig. 11

Primary resonance of a spur gear pair (tip relief start at 23.6 deg) at 170 Nm using second order perturbation and linearized model. Experimental data are from Fig. 3(b) of Ref. [5] (stable perturbation solution: solid line, unstable perturbation solution: dashed line, experimental measurement: circles, linear solution: dashed-dotted line).

## Tables

Table 1 Summary of nondimensional parameters and their physical interpretation for gear vibrations
Table 2 Physical parameters of the helical gears analyzed in Fig. 6
Table 3 Parameters of the analyzed gears in Figs. 6, 9, 10, and 11
Table 4 Physical parameters of the spur gears analyzed in Figs. 9, 10, and 11

## Discussions

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections