Advances in Acoustics and Vibration

Volume 2011, Article ID 583678, 13 pages

http://dx.doi.org/10.1155/2011/583678

## On the Boundary between Nonlinear Jump Phenomenon and Linear Response of Hypoid Gear Dynamics

^{1}Advanced Components Systems Division, Caterpillar Inc., Peoria, IL 61630, USA^{2}Mechanical Engineering Program, School of Dynamic Systems, University of Cincinnati, 598 Rhodes Hall, P.O. Box 210072, Cincinnati, OH 45221, USA

Received 27 October 2010; Revised 10 March 2011; Accepted 17 March 2011

Academic Editor: Mohammad Tawfik

Copyright © 2011 Jun Wang and Teik C. Lim. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A nonlinear time-varying (NLTV) dynamic model of a hypoid gear pair system with time-dependent mesh point, line-of-action vector, mesh stiffness, mesh damping, and backlash nonlinearity is formulated to analyze the transitional phase between nonlinear jump phenomenon and linear response. It is found that the classical jump discontinuity will occur if the dynamic mesh force exceeds the mean value of tooth mesh force. On the other hand, the propensity for the gear response to jump disappears when the dynamic mesh force is lower than the mean mesh force. Furthermore, the dynamic analysis is able to distinguish the specific tooth impact types from analyzing the behaviors of the dynamic mesh force. The proposed theory is general and also applicable to high-speed spur, helical and spiral bevel gears even though those types of gears are not the primary focus of this paper.

#### 1. Introduction

Gear dynamics have been studied intensively as evident from the discussions in [1–5]. In one of the earlier studies, Comparin and Singh investigated the nonlinear frequency response of an impact pair and located the transition conditions for no-impact, single-sided impact, and double-sided impact cases [1]. Later, Kahraman, and Singh examined the nonlinear dynamics of a spur gear pair [2] as well as a geared rotor-bearing system [3] and also studied the interaction between mesh stiffness and clearance nonlinearities [4]. More recently, Cheng and Lim studied the vibratory response of a hypoid geared rotor system with nonlinear time-varying mesh characteristics [5]. From those studies, it is now well known that nonlinear phenomena like jump discontinuity frequently occurs for lightly loaded gear pairs. As the drive load is progressively increased, the response becomes more linear until the jump response completely disappears. Understanding the conditions ahead of time when the geared system will exhibit nonlinear behaviors is highly desirable since the nonlinear dynamic analysis is typically very time consuming. However, to the best of our knowledge this transitional phase has never been studied extensively in the past or fully understood. Hence, the goal of this study is to seek out the conditions dictating the transitional boundary between nonlinear and linear responses. In other word, we would like to know precisely when jump phenomenon stops in advance as drive load goes from light to heavy.

To assist in the above-mentioned analysis, a nonlinear time-varying dynamic model with time-dependent mesh position, line-of-action mesh vectors, mesh stiffness, mesh damping, and backlash nonlinearity is proposed similar to the one employed by the authors in earlier studies [6, 7]. Based on the model, different levels of torques and damping ratios are applied so that the transition condition between linear and nonlinear responses can be located. The analysis clearly shows that when the dynamic mesh force is greater than mean mesh force, the jump phenomena will occur. On the other hand, if the dynamic mesh force is smaller in value than the mean mesh force, the jump behavior tends to disappear making the response appearing to be quite linear.

#### 2. Mathematical Model

The two degrees-of-freedom (DOF) lumped parameter torsional vibration model of a hypoid gear pair as shown in Figure 1(a), which was also applied in earlier studies by the authors [6, 7], is again assumed for the current study. This simplified representation is chosen over more complex system to allow for the study to focus on the behavior of the boundary between nonlinear jump phenomenon and linear dynamic response. The pinion and gear bodies each are allowed to rotate only. The mesh coupling is represented by a pair of mesh damping and stiffness elements. Using either the Newton’s or Lagrangian’s methods, the equations of motion can be easily derived as where and are the mass moments of inertias of pinion and gear, and are torque loads a on the pinion and gear, and are mesh stiffness and mesh damping coefficients, is unloaded kinematic transmission error, and is the gear backlash. The dynamic transmission error can be written as , while the directional rotation radius can be expressed as ( for pinion and gear resp.). Also, is the unit vector of the line of action along the mesh force direction in the coordinate system as shown in Figure 1(b), is the position vector of mesh point, and is the unit vector along the rotating axes of pinion and gear.

By assuming , (1a) and (1b) can be combined and simplified into the following form:where . Since mesh damping is generally time varying, proportional mesh damping is assumed here, which is more realistic than simply applying constant damping. Thus, the damping model can be expressed as where , .

The time-varying mesh parameters in (2a), (2b), and (2c) can be described as a superposition of a mean value and a sinusoidal component as shown below: The dynamic mesh force can be computed from From the torque and directional rotation radii, the mean value of the mesh force can be calculated from

By further assuming dimensionless parameters , , , , , , and , the dimensionless form of (2a) can be obtained aswhere , , and . The dimensionless mesh parameters are described as Subsequently the dimensionless dynamic mesh force can be derived as and the dimensionless mean value mesh force can be calculated from

The numerical integration method applying an explicit Runge-Kutta (4) and (5) formulation is employed to compute the response from (8a) since there is no analytical method available. The computed time domain response is then applied to calculate the dynamic mesh force using (10). From the predicted time history response function, the frequency spectrum can be obtained by taking the standard deviation (STD) and mean values of time series data at each operating frequency. From here onwards, the STD value is considered as the amplitude of dynamic mesh force.

#### 3. Dynamic Analysis

Using the proposed dynamic model, a typical rear axle hypoid gear set is analyzed. The basic system parameters of this rear axle unit are given in Table 1, which were also applied in an earlier numerical study by the authors [6]. Tooth contact analysis is first conducted to obtain the necessary mesh parameters for five different loads as shown in Table 2. The procedure for computing mesh parameters obtained from output results of a tooth contact analysis [8] is described in an earlier publication by the authors [9]. The extracted mesh parameters are readily employed in the dynamic model. For brevity, details of the process will not be repeated here.

##### 3.1. Frequency Domain Analysis

The dynamic mesh force response calculated using the above procedure for different mean loads is shown in Figure 2. Here mesh damping ratio of is used. As expected, when the torque load increases, the mean value of mesh force increases too. In fact, the mean value of mesh force response is very close to the theoretical static mesh force derived from the input torque load and is almost frequency invariant compared to the dynamic mesh force that varies significantly with frequency. Also, the dynamic mesh force appears to decrease with increasing input torque. A more in-depth analysis of the dynamic mesh force results are discussed next.

At , several jump discontinuities and very rich single and double-sided tooth impacts can be observed at the primary resonance and higher frequencies. As is increased to 1.37, only two jumps are seen in which one is a softening type (lower frequency) and the other is a hardening one (upper frequency). At the same time, less tooth impacts occur around the primary resonance. For the case of only a softening type jump, and very few single and double-sided tooth impact points appear near the primary resonance. At there is no jump response observed and only a single-sided tooth impact is found. As is increased further to 3.27, the dynamic response over the entire frequency range shown is completely void of jump phenomenon and tooth impacts.

From the results, an interesting finding is that when the dynamic mesh force is greater than the mean value of mesh force or the theoretical static mesh force, jump discontinuity can be seen. On the other hand, when the dynamic mesh force is lower than the mean value of mesh force or the theoretical static mesh force, virtually no jump discontinuity is seen. This is because the mean value of mesh force is able to keep the meshing teeth engaged all the time when its magnitude is higher than the dynamic counterpart that tends to cause the meshing teeth to lose contact. The results of mean value of mesh force response for different levels of input torque load are shown as a closed-up view in Figure 3. Again, as seen in these results, as the drive torque load is increased, the mean mesh force level increases. However, its variation decreases, yielding a more linear characteristic.

The effect of damping is analyzed next. Figure 4 shows the dynamic mesh force response for different mesh damping ratios. From these plots, it can be seen that as is increased, the dynamic mesh force tends to decrease, while mean value does not change at all. This is consistent with the results depicted in (11). At the same time, occurrence of jump discontinuity and tooth impacts lessen. Also, the mean value of mesh force is very close to the theoretical static mesh force and changes very little. Similarly, the results again show that as dynamic mesh force becomes larger than the mean value of mesh force, jump response appears, while no evidence of jump response is seen when the dynamic mesh force is lower than the mean value.

The effect of mesh damping on the mean value of mesh force is shown in Figure 5. Note that these plots represent a closed-up view of the ones shown in Figure 4. Here, it can be seen that as mesh damping ratio is increased, the mean value changes very little and its variation tends to decrease. This former result is unlike the previously seen effect of drive torque load, but the latter observation is the same as the trend of drive torque load effect.

##### 3.2. Time Domain Analysis

Due to the design of backlash in the hypoid gear set, there are three classes of tooth meshing cases that are no impact, single-sided impact, and double-sided impact as illustrated in Figure 6. The zone between −1 and 1 is backlash zone. So there are three types of no impact (, , and ), two types of single-sided impact ( and ), and one type of double-sided impact (*D*). These tooth impact behaviors were first reported by Comparin and Singh in 1989 [1] for parallel axis gears. For nonparallel axis gears, there are no known prior studies. The current study attempts to fill this gap using the hypoid gear design as the example case. In the ensuing discussion, let us consider the predicted dynamic mesh force spectrum for as shown in Figure 7.

For the case of as shown in Figure 7, the time history plots of dynamic response for increasing mesh frequencies from 0.9 to 1.1 are shown in a set of plots in Figure 8. In the plot given by Figure 8(a) for the case of , the time history response shown is of type where no tooth impact response occurs. The corresponding time history plot of the dynamic mesh force at this frequency is shown in Figure 8(f). From this figure, it is seen that the trough value of dynamic mesh force is above zero. To explain this observation, consider (10) and (8b) that show the dynamic mesh force being dominated by . Since is greater than zero for type no impact response, the dynamic mesh force must therefore be greater than zero. Figures 8(b) and 8(g) for shows the time history response of type *D* double-sided tooth impact behavior. It shows that the trough value of dynamic mesh force reaches below zero because is less than zero for . Figure 8(c) shows the time history of type single-sided tooth impact response at , with the corresponding dynamic mesh force shown in Figure 8(h). Here, the trough value of dynamic mesh force is equal to zero because for . This is found to be similar for the time history plots of type single-sided tooth impact response and dynamic mesh force at , as shown in Figures 8(d) and 8(i), respectively. At as shown in Figures 8(e) and 8(j), the response goes back to type no tooth impact behavior.

Time history plots of dynamic response for decreasing mesh frequencies from 1.1 to 0.85 are shown in Figure 9. In this set of plots, Figure 9(a) shows time history of type no tooth impact response at and Figure 9(g) shows the corresponding time history of dynamic mesh force in which the trough value is above zero. As is decreased to 1.05, the gear response becomes type single-sided tooth impact as depicted by the time history response in Figures 9(b) and 9(h) in which the trough value is equal to zero. Similar type single-sided tooth impact response occurs at as shown in Figures 9(c) and 9(i). At and 0.9, double-sided tooth impacts occur and the trough value of dynamic mesh force reaches below zero as shown in Figures 9(d), 9(j), 9(e), and 9(k), respectively. Note that the negative mesh force implies a reversal of its vector along the line of action due to the tooth in question back-colliding with the preceding tooth of the mating gear. As is further decreased to 0.85 given by Figures 9(f) and 9(l), type no tooth impact response appears and the trough value of dynamic mesh force is well above zero.

From the above analysis, it is shown that for type no tooth impact case, the trough value of dynamic mesh force is greater than zero; for type single-sided tooth impact case, the trough value of dynamic mesh force is equal to zero; and finally for *D* type double-sided tooth impact case, the trough value of dynamic mesh force reaches below zero. Furthermore, it can be concluded that for type no tooth impact case (where no mesh contact at all), the dynamic mesh force will always be equal to zero; for type no tooth impact case (where the tooth-tooth engagement happens between the current and preceding tooth of the mating gear), the peak value of dynamic mesh force will be less than zero implying a mesh vector in the opposite direction from the case; for type single-sided tooth impact case, the peak value of dynamic mesh force will be equal to zero since the current tooth in question never makes contact with the tooth on the mating gear it is designed to engage with.

#### 4. Concluding Remarks

A nonlinear time-varying torsional vibration model of a hypoid gear pair system with time-dependent mesh point, line-of-action vector, mesh stiffness, mesh damping, and backlash nonlinearity is formulated to study the condition dictating the boundary between nonlinear jump phenomena and linear responses, that is, when jump discontinuity begins to disappear as mean load is increased from light to heavy levels. An interesting finding is that when the dynamic mesh force exceeds the mean value of mesh force, jump discontinuity will appear, but when the dynamic mesh force is lower than the mean value of mesh force, jump response is less likely to occur.

Based on the time domain analysis results, it is shown that for type no tooth impact case, the trough value of dynamic mesh force is greater than zero; for type no tooth impact case, the dynamic mesh force is always equal to zero; for type no tooth impact case, the peak value of dynamic mesh force is less than zero; for type single-sided tooth impact case, the trough value of dynamic mesh force is equal to zero; for type single-sided tooth impact case, the peak value of dynamic mesh force is equal to zero; for *D* type double-sided tooth impact case, the dynamic mesh force reaches both above and below zero. From these observations, it is clear that different tooth impact types can be distinguished by the dynamic mesh force response behavior.

#### Nomenclature

: | Gear backlash |

: | Mesh damping coefficient |

: | Unloaded kinematic transmission error |

: | Nonlinear displacement function |

, : | Mass moments of inertias of pinion and gear |

: | Unit vector along pinion or gear rotating axis |

, : | Mesh stiffness |

: | Equivalent mass |

: | Unit normal vector of mesh point |

: | Difference between dynamic and kinematic transmission errors |

: | Position vector of mesh point |

: | Coordinate system for dynamic formulation |

: | Time |

, : | Mean value of pinion and gear torque loads |

: | Dynamic transmission error |

: | Directional rotation radius |

: | Excitation frequency |

: | Mesh damping ratio |

: | Dimensionless torque ratio |

: | Phase angle. |

*Subscripts*

: | Label for pinion () and gear () |

: | Fundamental harmonic |

: | Mean |

: | Natural. |

*Superscripts*

~: | Dimensionless quantities |

: | Vector quantities |

′: | Derivative with respect to time. |

#### References

- R. J. Comparin and R. Singh, “Non-linear frequency response characteristics of an impact pair,”
*Journal of Sound and Vibration*, vol. 134, no. 2, pp. 259–290, 1989. View at Google Scholar · View at Scopus - A. Kahraman and R. Singh, “Non-linear dynamics of a spur gear pair,”
*Journal of Sound and Vibration*, vol. 142, no. 1, pp. 49–75, 1990. View at Google Scholar · View at Scopus - A. Kahraman and R. Singh, “Non-linear dynamics of a geared rotor-bearing system with multiple clearances,”
*Journal of Sound and Vibration*, vol. 144, no. 3, pp. 469–506, 1991. View at Google Scholar · View at Scopus - A. Kahraman and R. Singh, “Interactions between time-varying mesh stiffness and clearance non-linearities in a geared system,”
*Journal of Sound and Vibration*, vol. 146, no. 1, pp. 135–156, 1991. View at Google Scholar · View at Scopus - Y. Cheng and T. C. Lim, “Dynamics of hypoid gear transmission with nonlinear time-varying mesh characteristics,”
*Journal of Mechanical Design, Transactions of the ASME*, vol. 125, no. 2, pp. 373–382, 2003. View at Publisher · View at Google Scholar · View at Scopus - J. Wang, T. C. Lim, and M. Li, “Dynamics of a hypoid gear pair considering the effects of time-varying mesh parameters and backlash nonlinearity,”
*Journal of Sound and Vibration*, vol. 308, no. 1-2, pp. 302–329, 2007. View at Publisher · View at Google Scholar · View at Scopus - J. Wang and T. C. Lim, “Effect of tooth mesh stiffness asymmetric nonlinearity for drive and coast sides on hypoid gear dynamics,”
*Journal of Sound and Vibration*, vol. 319, no. 3-5, pp. 885–903, 2009. View at Publisher · View at Google Scholar · View at Scopus - S. Vijayakar,
*Contact Analysis Program Package: Calyx*, Advanced Numerical Solutions, Hilliard, Ohio, 2003. - T. C. Lim and J. Wang, “Effects of assembly errors on hypoid gear mesh and dynamic response,” in
*Proceedings of the ASME Power Transmission and Gearing Conference (DETC '05)*, Long Beach, Calif, USA, 2005.