Nonlinear Dynamics of Driveline Systems with Hypoid Gear Pair
A dissertation submitted to the
Division of Research and Advanced Studies
of the University of Cincinnati
in partial fulfillment of the requirements
for the degree of
DOCTOR OF PHILOSOPHY
in the Program of School of Dynamic Systems
of the College of Engineering and Applied Science
April 2012
by
Junyi Yang
M. S. Southeast University, Nanjing, P. R. China, 2007
B. S. Southeast University, Nanjing, P. R. China, 2004
Academic Committee Chair: Dr. Teik C. Lim
Members: Dr. Dong Qian
Dr. David F. Thompson
Dr. Murali M. Sundaram
ABSTRACT
This dissertation research focuses on evaluating the nonlinear dynamics of driveline systems employed in motor vehicles with emphasis on characterizing the excitations and response of right-angle, precision hypoid-type geared rotor structure. The main work and contribution of this dissertation is divided into three sections. Firstly, the development of an asymmetric and nonlinear gear mesh coupling model will be discussed. Secondly, the enhancement of the multi-term harmonic balance method (HBM) is presented. Thirdly and as the final topic, the development of new dynamic models capable of evaluating the dynamic coupling characteristics between the gear mesh and other driveline structures will be addressed.
A new asymmetric and nonlinear mesh model will be proposed that considers backlash, and the fact that the tooth surfaces of the convex and concave sides are different. The proposed mesh model will then be fed into a dynamic model of the right-angle gear pair to formulate the dimensionless equation of motion of the dynamic model. The multi-term HBM will be enhanced to simulate the right-angle gear dynamics by solving the resultant dimensionless equation of motion. The accuracy of the enhanced HBM solution will be verified by comparison of its results to the more computationally intensive direct numerical integration calculations. The stability of both the primary and sub-harmonic solutions predicted by applying multi-term HBM will be analyzed using the Floquent Theory. In addition, the stability analysis of the multi-term HBM solutions will be proposed as an approximate approach for locating the existence of sub- harmonic and chaotic motions.
In this dissertation research, a new methodology to evaluate the dynamic interaction between the nonlinear hypoid gear mesh mechanism and the time-varying characteristics of the rolling element bearings will also be developed. The time-varying mesh parameters will be
II
obtained by synthesizing a 3-dimensional loaded tooth contact analysis (TCA) results. The time- varying stiffness matrix approach will be used to represent the dynamic characteristics of the rolling element bearings.
An overall nonlinear dynamic model of the hypoid gear box considering elastic housing
structure will be developed as well. A lumped parameter model of the flexible housing will be
extracted form an appropriate set of frequency response functions through modal parameter identification method. In order to obtain the rotational coordinates, a rigid body interpolation of
the translational responses at the bearing locations on the housing structure will be applied. The
reduced model will be then coupled with the hypoid gear-shaft-bearing assembly model by
applying a proposed dynamic coupling procedure.
Finally, a hypoid geared rotor system model considering the propeller shaft flexibility
will be established. The propeller shaft bending flexibility will be modeled as lumped parameter
model through using the component mode synthesis (CMS). The torsional flexibility of propeller
shaft will be simplified as a torsional spring connecting the inertia of moment of engine and
pinion. Physically, the pinion input shaft is driven by the propeller shaft through a universal
joint, which will be modeled as a flexible simple supported boundary condition as well as
fluctuating rotation speed and torque excitation.
III
IV
ACKNOWLEDGEMENTS
I would like to thank Prof. Teik C. Lim, who is serving as my academic advisor and the chair of my academic committee, for his great instructions and support throughout my graduate study. I would also like to thank Dr. Dong Qian, Dr. David F. Thompson and Dr. Murali M.
Sundaram for serving as my doctoral academic committee members.
I wish to thank all my colleagues at the Vibro-Acoustic and Sound Quality Research
Laboratory in University of Cincinnati for their friendship. I would like to express my special gratitude to Dr. Tao Peng for his valuable academic suggestions.
Finally, I would like to thank my parents and my wife Si Chen for their support and patience during my graduate study.
V
TABLE OF CONTENTS
ABSTRACT II
ACKNOWLEDGEMENTS ...... V
TABLE OF CONTENTS ...... VI
LIST OF TABLES ...... X
LIST OF FIGURES ...... XI
LIST OF SYMBOLS ...... XVIII
Chapter 1. Introduction ...... 1
1.1 Literature Review ...... 1
1.2 Scope and Objectives ...... 6
1.3 Organization ...... 8
Chapter 2. A Review of Loaded Tooth Contact Analysis Approaches ...... 12
2.1 Introduction ...... 12
2.2 Gear Body Flexible Deflection ...... 14
2.3 Gear Tooth Flexible Deflection ...... 17
2.3.1 Cantilever Beam ...... 17
2.3.2 Cantilever Plate ...... 18
2.3.3 Shell ...... 21
2.3.4 FEM Method ...... 23
2.3.5 Finite Strip Method...... 26
2.3.6 Finite Prism Method ...... 27
2.3.7 Combined Surface Integral and Finite Element Method ...... 27
2.4 Conclusion ...... 29 VI
Chapter 3. An Enhanced Multi-term Harmonic Balance Solution for Non-linear Period-one
Dynamic Motions in Right-angle Gear Pairs ...... 31
3.1 Introduction ...... 31
3.2 Mathematical Model ...... 34
3.3 Period-one Dynamics ...... 42
3.4 Parametric Studies ...... 46
3.4.1 Numerical Validation ...... 47
3.4.2 Numerical Analysis ...... 50
3.5 Conclusion ...... 58
Chapter 4. An Enhanced Multi-term Harmonic Balance Solution for Non-linear Period-
Dynamic Motions in Right-angle Gear Pairs ...... 59
4.1 Introduction ...... 59
4.2 Period- Sub-harmonic Response ...... 62
4.3 Results and Discussion ...... 65
4.3.1 Comparison of HBM and numerical integration results ...... 66
4.3.2 Effect of Static Load ...... 75
4.3.3 Effect of Static Transmission Error Excitation ...... 80
4.3.4 Sub-harmonic and Chaotic Motions ...... 83
4.4 Conclusion ...... 84
Chapter 5. Dynamics of Coupled Nonlinear Hypoid Gear Mesh and Time-varying Bearing
Stiffness Systems ...... 86
5.1 Introduction ...... 86
5.2 Nonlinear Geared Rotor System Model ...... 88
VII
5.2.1 Basic Assumptions ...... 88
5.2.2 Mesh Model ...... 89
5.2.3 Shaft-bearing Assembly Model ...... 92
5.2.4 Time-varying Bearing Stiffness ...... 93
5.2.5 Formulation ...... 95
5.3 Case Study ...... 97
5.3.1 Mesh and Bearing Interaction ...... 99
5.3.2 Backlash and Bearing Interaction ...... 100
5.4 Summary ...... 104
Chapter 6. Nonlinear Dynamic Simulation of Hypoid Gearbox with Elastic Housing ...... 105
6.1 Introduction ...... 105
6.2 Geared Rotor System Dynamic Model ...... 107
6.3 Lumped Parameter Model of Housing ...... 109
6.4 Coupled Dynamic Model ...... 115
6.5 Case Study ...... 116
6.5.1 Methodology Validation ...... 118
6.5.2 Parametric Study ...... 120
6.6 Conclusion ...... 125
Chapter 7. Propeller Shaft Bending and Effect on Gear Dynamics ...... 126
7.1 Introduction ...... 126
7.2 Mathematical Model ...... 128
7.2.1 Simplified Propeller Shaft Model ...... 128
7.2.2 Geared Rotor System Model ...... 130
VIII
7.2.3 Coupled Boundary Condition ...... 133
7.3 Numerical Results and Discussion ...... 134
7.4 Conclusion ...... 141
Chapter 8. Conclusions and Proposed Future Studies ...... 143
8.1 Conclusions ...... 143
8.2 Proposed Future Studies ...... 145
BIBLIOGRAPHY ...... 147
IX
LIST OF TABLES
Table 3.1, Dimensionless dynamic parameters for a typical automotive hypoid gear pair ...... 47
Table 3.2, Physical parameters of a real application hypoid gear pair ...... 48
Table 4.1, Dimensionless dynamic parameters for a light load automotive hypoid gear pair ...... 66
Table 5.1, Design parameters: (a) Gear data, (b) System data, (c) Data of bearing on pinion shaft,
(d) Data of bearing on gear shaft ...... 97
Table 6.1, Design parameter: (a) Gear data, (b) System data, (c) Data of bearing on pinion shaft,
(d) Data of bearing on gear shaft ...... 117
Table 6.2, Stiffness and damping of the analytical model ...... 119
Table 7.1, Geometry dimension of propeller shaft ...... 134
Table 7.2, Design parameter: (a) Gear data, (b) System data, (c) Data of bearing on pinion shaft,
(d) Data of bearing on gear shaft ...... 134
X
LIST OF FIGURES
Fig. 2.1, Linear normal stress distribution (Weber, 1949) ...... 14
Fig. 2.2, Schematic of test rig (O’Donnell, 1960) ...... 15
Fig. 2.3, Finite element model of gear tooth (Stegemiller and Houser, 1993) ...... 17
Fig. 2.4, Cantilever beam model of gear tooth (Weber, 1949) ...... 18
Fig. 2.5, Cantilever plate model of gear tooth (MacGregor, 1936) ...... 18
Fig. 2.6, Moment imaging: (a) plate of infinite length, (b) plate of finite length (Wellauer et al.,
1960) ...... 19
Fig. 2.7, Tapered cantilever plate model of gear tooth (Yau, et al., 1994) ...... 20
Fig. 2.8, Annual sector plate model of gear tooth (Vaidyanathan, et al., 1994) ...... 21
Fig. 2.9, Thick cylinder circular shell model of gear tooth (Vaidyanathan, et al. 1993) ...... 22
Fig. 2.10, Discretization of instant contact line (Wilcox, 1981) ...... 24
Fig. 2.11, Finite mesh of gear tooth (Gosselin, et al., 1995)...... 26
Fig. 2.12, Finite strip model of gear tooth (Gagnon, et al., 1996) ...... 26
Fig. 2.13, Finite prism element (Vijayakar and etc., 1987) ...... 28
Fig. 2.14, Two-dimensional formulation of finite strip element (Vijayakar and etc., 1989) ...... 28
Fig. 3.1, (a) Two degrees of freedom torsional vibration model of a hypoid gear pair, (b) Pinion
and gear coordinate systems ...... 35
Fig. 3.2, Asymmetric nature of hypoid gear mesh stiffness: (a) Drive side (b) Coast side ...... 38
Fig. 3.3, Comparison of multi-term HBM and numerical integration: (a) RMS of the dynamic
displacement. (b) Mean value of the dynamic displacement. Stable steady-state
solution by multi-term HBM; Unstable steady-state solution by multi-terms HBM; ○
Solutions by numerical integration ...... 49
XI
Fig. 3.4, Effect of the directional rotation radii variation on the dynamic response: (a) RMS of
the dynamic displacement. (b) Primary resonance of the dynamic displacement.
Stable steady-state solution by multi-term HBM; Unstable steady-state solution by
multi-term HBM ...... 52
Fig. 3.5, Effect of directional rotation radii asymmetric nonlinearity on the dynamic response: (a)
RMS of the dynamic displacement. (b) Primary resonance of the dynamic displacement.
Stable steady-state solution by multi-term HBM; Unstable steady-state
solution by multi-term HBM ...... 53
Fig. 3.6, Effect of the mesh stiffness variation on the dynamic response: (a) RMS of the dynamic
displacement. (b) Primary resonance of the dynamic displacement. Stable steady-
state solution by multi-term HBM; Unstable steady-state solution by multi-term
HBM ...... 55
Fig. 3.7, Effect of the mesh stiffness asymmetric nonlinearity on the dynamic response: (a)
RMS of the dynamic displacement. (b) Primary resonance of the dynamic displacement.
Stable steady-state solution by multi-term HBM; Unstable steady-state
solution by multi-term HBM ...... 57
Fig. 4.1, RMS of dynamic transmission error. Stable steady solution by multi-terms
HBM; Unstable steady solution by multi-terms HBM; ○ Solutions by Numerical
Integration ...... 67
Fig. 4.2, Bifurcation diagram of baseline case (frequency sweep up) ...... 68
Fig. 4.3, a) Period 1 HBM solutions Stable steady solution; Unstable steady
solution; Δ Compared HBM solution; b)FFT of NI result; c) Poincare plot of NI result; d)
Trajectory of NI result ...... 69
XII
Fig. 4.4, a) Period 2 HBM solutions Stable steady solution; Unstable steady
solution; Δ Compared HBM solution; b)FFT of NI result; c) Poincare plot of NI result; d)
Trajectory of NI result ...... 71
Fig. 4.5, a) Period 3 HBM solutions Stable steady solution; Unstable steady
solution; Δ Compared HBM solution; b)FFT of NI result; c) Poincare plot of NI result; d)
Trajectory of NI result ...... 72
Fig. 4.6, a) Period 2 HBM solutions Stable steady solution; Unstable steady
solution; Δ Compared HBM solution; b)FFT of NI result; c) Poincare plot of NI result; d)
Trajectory of NI result ...... 73
Fig. 4.7, a) Period 3 HBM solutions Stable steady solution; Unstable steady
solution; Δ Compared HBM solution; b)FFT of NI result; c) Poincare plot of NI result; d)
Trajectory of NI result ...... 74 ~ ~ Fig. 4.8, a) Period 2 HBM solutions (Tp 0. 3 ); b) Period 2 HBM solutions (Tp 0. 5 ); c)
~ ~ Bifurcation diagram (Tp 0. 3 ); d) Bifurcation diagram (Tp 0. 5 ); Stable steady
solution; Unstable steady solution ...... 77
Fig. 4.9, a) Trajectory on phase plane (~ 1.33 ); b) Poincare plot (~ 1.33 ); c) Trajectory on
phase plane (~ 1.45 ); d) Poincare plot (~ 1.45 ); e) Trajectory on phase plane
(~ 1.51 ); f) Poincare plot (~ 1.51 );g) Trajectory on phase plane (~ 1.69 ); h) Poincare
plot (~ 1.69 ); i) Trajectory on phase plane (~ 1.95 ); j) Poincare plot (~ 1.95 ) ...... 79 ~ ~ Fig. 4.10, a) Period 2 HBM solutions (e3 0. 225 ); b) Period 2 HBM solutions (e3 0. 115 );
~ ~ c) Bifurcation diagram (e3 0. 225 ); d) Bifurcation diagram (e3 0. 115 ); Stable
steady solution; Unstable steady solution ...... 81
XIII
Fig. 4.11, a) Trajectory on phase plane (~ 1.1); b) Poincare plot (~ 1.1); c) Trajectory on
phase plane (~ 1.38 ); d) Poincare plot (~ 1.38 ); e) Trajectory on phase plane
(~ 1.91 ); f) Poincare plot (~ 1.91 ) ...... 83
Fig. 4.12, a) Existence spaces of sub-harmonic motions; b) Existence spaces of chaotic motions;
Ο start point; + end point ...... 84
Fig. 4.13, a) Existence spaces of sub-harmonic motions; b) Existence spaces of chaotic motions;
Ο start point; + end point ...... 84
Fig. 5.1, Schematic of a 14-DOF nonlinear dynamic model of hypoid geared rotor system...... 89
Fig. 5.2, Loaded tooth contact analysis model: (a) Gear pair geometry, and (b) Contact cells on
engaging tooth surface...... 90
Fig. 5.3, Shaft bearing assembly: (a) Physical structure, (b) Beam finite element representation.
...... 92
Fig. 5.4, Rolling element bearing kinematics and the corresponding coordinate systems...... 93
Fig. 5.5, Shaft-bearing assembly: (a) Pinion, (b) Gear...... 97
Fig. 5.6, Time-varying bearing stiffness ...... 98
Fig. 5.7, Time-varying meshes stiffness ...... 99
Fig. 5.8, Dynamic mesh force comparison: time-invariant bearing stiffness, time-
varying bearing stiffness...... 100
Fig. 5.9, Comparison of bearing axial reaction force: time-invariant bearing stiffness;
time-varying bearing stiffness...... 100
Fig. 5.10, Dynamic mesh force: frequency sweep up; frequency sweep down.
...... 101
XIV
Fig. 5.11, Dynamic mesh force comparison: time-invariant bearing stiffness;
time-varying bearing stiffness ...... 102
Fig. 5. 12, Comparison of bearing radial reaction force: time-invariant bearing stiffness;
time-varying bearing stiffness ...... 102
Fig. 5.13, Dynamic mesh force comparison: with backlash; without backlash. 103
Fig. 5. 14, Comparison of bearing reaction moment: with backlash; without
backlash...... 103
Fig. 6.1, Dynamic Model of hypoid geared rotor assembly ...... 107
Fig. 6.2, a) Housing continuous structure; b) Excitation and response points ...... 110
Fig. 6.3, Frequency response function of housing structure ...... 111
Fig. 6.4, Poles estimation from FRFs ...... 112
Fig. 6.5, Plot of Enhanced frequency response function; and Single-DOF
frequency response function fit ...... 113
Fig. 6.6, Illustration of 2-D rigid body interpolation ...... 113
Fig. 6.7, Shaft-bearing assembly: (a) Pinion, (b) Gear ...... 117
Fig. 6.8, Lumped parameter model representing a housing structure ...... 118
Fig. 6.9, Comparison of dynamic mesh force: Full coupling; Partial coupling;
No coupling...... 120
Fig. 6.10, Acceleration of mass block 1: Full coupling; Partial coupling ...... 120
Fig. 6.11, Comparison of dynamic mesh force: 10-modes coupling; 5-modes
coupling; No coupling...... 121
Fig. 6.12, Comparison of surface acceleration: 10-modes coupling; 5-modes
coupling ...... 122
XV
Fig. 6.13, Comparison of dynamic mesh force: No external excitation; Excitation
exerted on point 1; Excitation exerted on point 2...... 123
Fig. 6.14, FFT of dynamic mesh force with external excitation on point 1...... 124
Fig. 6.15, FFT of dynamic mesh force with external excitation on point 2...... 124
Fig. 6.16, Comparison of surface acceleration: No external excitation; Excitation
exerted on point 1; Excitation exerted on point 2 ...... 125
Fig. 7.1, Simplified propeller shaft (a) beam with different cross-section, (b) Cross-section
dimension definition ...... 128
Fig. 7.2, (a) Static system with unit forces applied at interface, (b) Dynamic system with fixed
interface boundary condition ...... 130
Fig. 7.3, Lumped parameter model of hypoid geared rotor system ...... 130
Fig. 7.4, Schematic of propeller shaft and pinion coupling ...... 133
Fig. 7.5, Shaft-bearing assembly: (a) Pinion, (b) Gear ...... 134
Fig. 7.6, Comparison of dynamic mesh force: Rigid Propeller Shaft; Propeller
Shaft Bending (L=0.1); Propeller Shaft Bending (L=0.15) ...... 136
Fig. 7.7, Modal shape of the linear model ...... 136
Fig. 7.8, Comparison of reaction forces on pinion bearing: Rigid Propeller Shaft;
Propeller Shaft Bending (L=0.1) (a) radial force, (b) axial force, (c) moment ...... 138
Fig. 7.9, Comparison of reaction forces on gear bearing: Rigid Propeller Shaft;
Propeller Shaft Bending (L=0.1) (a) radial force, (b) axial force, (c) moment ...... 139
Fig. 7.10, Effect of coupling stiffness: (a) comparison of dynamic mesh force, (b) affected
resonance Baseline Coupling Stiffness ( ); Low Coupling
Stiffness ( ); High Coupling Stiffness ( )...... 140
XVI
Fig. 7.11, Effect of coupling damping: (a) comparison of dynamic mesh force, (b) affected
resonance Baseline Coupling Damping ( ); High Coupling
Damping ( ); High Coupling Damping ( ) ...... 140
Fig. 7.12, Comparison of dynamic mesh force: Baseline Modal Damping Ratio (