EFFECT OF SLIDING FRICTION ON SPUR AND HELICAL GEAR DYNAMICS AND VIBRO-ACOUSTICS
DISSERTATION
Presented in Partial Fulfillment of the Requirements for
the Degree Doctor of Philosophy in the Graduate School
of The Ohio State University
By
Song He, B.S., M.S.
* * * * *
The Ohio State University
2008
Dissertation Committee: Approved by Dr. Rajendra Singh, Advisor
Dr. Ahmet Kahraman ______
Dr. Ahmet Selamet Adviser
Dr. Marcelo Dapino Graduate Program in Mechanical Engineering
ABSTRACT
This study examines the salient effects of sliding friction on spur and helical gear dynamics and associated vibro-acoustic sources. First, new dynamic formulations are developed for spur and helical gear pairs based on a periodic description of the contact point and realistic mesh stiffness. Difficulty encountered in the existing discontinuous models is overcome by characterizing a smoother transition during the contact. Frictional forces and moments now appear as either excitations or periodically-varying parameters, since the frictional force changes direction at the pitch point/line. These result in a class of periodic ordinary differential equations with multiple and interacting coefficients, which characterize the effect of sliding friction in spur or helical gear dynamics.
Predictions (based on multi-degree-of-freedom analytical models) match well with a benchmark finite element/contact mechanics code and/or experimental results.
Second, new analytical solutions are constructed which provide an efficient evaluation of the frictional effect as well as a more plausible explanation of dynamic interactions in multiple directions. Both single- and multi-term harmonic balance methods are utilized to predict dynamic mesh loads, friction forces and pinion/gear displacements. Such semi-analytical solutions explain the presence of higher harmonics in gear noise and vibration due to exponential modulations of the periodic stiffness,
ii dynamic transmission error and sliding friction. This knowledge also analytically reveals the effect of the tooth profile modification in spur gears on the dynamic transmission error, under the influence of sliding friction. Further, the Floquet theory is applied to obtain closed-form solutions of the dynamic response for a helical gear pair, where the effect of sliding friction is quantified by an effective piecewise stiffness function.
Analytical predictions, under both homogeneous and forced conditions, are validated using numerical simulations. The matrix-based methodology is found to be computationally efficient while leading to a better assessment of the dynamic stability.
Third, an improved source-path-receiver vibro-acoustic model is developed to quantify the effect of sliding friction on structure-borne noise. Interfacial bearing forces are predicted for the spur gear source sub-system given two gear whine excitations (static transmission error and sliding friction). Next, a computational model of the gearbox, with embedded bearing stiffness matrices, is developed to characterize the motilities of structural paths. Radiated sound pressure is then estimated by using two numerical techniques (the Rayleigh integral method and a substitute source technique). Predicted pressures match well with measured noise data over a range of operating torques. In particular, the proposed vibro-acoustic model quantifies the contribution of sliding friction, which could be significant when the transmission error is minimized through tooth modifications.
iii
DEDICATION
Dedicated to my parents and wife
iv ACKNOWLEDGMENT
I would like to express my sincere appreciation to my advisor, Professor Rajendra
Singh, for his time, guidance, and support over the years both in my academic research and personal life. His intellectual insight and encouragement had a huge impact on my
professional growth. I would also like to thank Professor Ahmet Kahraman, Professor
Ahmet Selamet and Professor Marcelo Dapino for their services on the doctoral committee and for providing constructive suggestions.
I sincerely thank Dr. Todd Rook for providing valuable suggestions. I greatly
appreciate Dr. Rajendra Gunda for granting access to the Calyx software and for offering
insightful comments. I gracefully acknowledge the experimental work conducted by
Vivake Asnani and Fred Oswald, as well as the collaboration with Allison Lake. Dr.
Chengwu Duan is thanked for helping me in both my research and personal life. All colleagues in the Acoustics and Dynamic Laboratory are acknowledged for their encouragement and friendship. I thank Professor Goran Pavić, Professor Jean-Louis
Guyader and Corinne Lotto for their advice and kind help during my stay in INSA Lyon.
The financial support from the Army Research Office, EU’s Marie Curie
Fellowship and OSU Presidential Fellowship is gracefully appreciated.
Finally, I would like to thank my parents and my wife, Lihua, for their love and encouragement throughout my pursuit for the doctoral degree.
v VITA
October 24, 1979………………………..…..….Born – Jiangsu, China
2002………………………..………………..….B.S. Instrumentation Engineering Shanghai Jiao Tong University Shanghai, China
2004………………………..………………..….M.S. Mechanical Engineering The Ohio State University
2004 - 2006………………..………………..….Graduate Teaching and Research Associate, Mechanical Engineering The Ohio State University
2007………………………..………………..….Marie Curie Fellow (EU) National Institute of Applied Science Lyon, France
2007 - 2008………………………..……………Presidential Fellow (Graduate School) The Ohio State University
PUBLICATIONS
Research Publications
1. He, S., Gunda, R., and Singh, R., 2007, “Effect of Sliding Friction on the Dynamics of Spur Gear Pair with Realistic Time-Varying Stiffness,” Journal of Sound and Vibration, 301, pp. 927-949.
2. He, S., Gunda, R., and Singh, R., 2007, “Inclusion of Sliding Friction in Contact Dynamics Model for Helical Gears, ASME Journal of Mechanical Design, 129(1), pp. 48-57.
vi 3. He, S., Cho, S., and Singh, R., 2008, “Prediction of Dynamic Friction Forces in Spur Gears using Alternate Sliding Friction Formulations,” Journal of Sound and Vibration, 309(3-5), pp. 843-851.
FIELDS OF STUDY
Major Field: Mechanical Engineering Dynamics of Mechanical Systems Vibro-Acoustics and Noise Control
vii TABLE OF CONTENTS
Page Abstract...... ii Dedication...... iv Acknowledgment ...... v Vita...... vi
List of Tables ...... xii List of Figures...... xiii List of Symbols...... xx
Chapter 1 Introduction ...... 1 1.1 Motivation...... 1 1.2 Literature Review...... 4 1.3 Problem Formulation ...... 8 1.3.1 Key Research Issues ...... 8 1.3.2 Scope, Assumptions and Objectives...... 11 References for Chapter 1 ...... 14
Chatper 2 Spur Gear Dynamics with Sliding Friction and Realistic Mesh Stiffness ...... 18 2.1 Introduction...... 18 2.2 Problem Formulation ...... 20 2.2.1 Objectives and Assumptions...... 20 2.2.2 Timing of Key Meshing Events...... 23 2.2.3 Calculation of Realistic Time-Varying Tooth Stiffness Functions...... 24 2.3 Analytical Multi-Degree-of-Freedom Dynamic Model...... 27 2.3.1 Shaft and Bearing Stiffness Models...... 27 2.3.2 Dynamic Mesh and Friction Forces...... 28 2.3.3 MDOF Model...... 32 2.4 Analytical SDOF Torsional Model...... 33 viii 2.5 Effect of Sliding Friction in Example I...... 35 2.5.1 Validation of Example I Model using the FE/CM Code ...... 35 2.5.2 Effect of Sliding Friction ...... 40 2.5.3 MDOF System Resonances ...... 42 2.6 Effect of Sliding Friction in Example II ...... 43 2.6.1 Empirical Coefficient of Friction...... 43 2.6.2 Effect of Tip Relief on STE and k(t)...... 45 2.6.3 Phase Relationship between Normal Load and Friction Force Excitations... 50 2.6.4 Prediction of the Dynamic Responses ...... 53 2.7 Experimental Validation of Example II Models...... 56 2.8 Conclusion ...... 65
Chapter 3 Prediction of Dynamic Friction Forces Using Alternate Formulations ...... 68 3.1 Introduction...... 68 3.2 MDOF Spur Gear Model ...... 69 3.3 Spur Gear Model with Alternate Sliding Friction Models ...... 74 3.3.1 Model I: Coulomb Model ...... 74 3.3.2 Model II: Benedict and Kelley Model ...... 75 3.3.3 Model III: Formulation Suggested by Xu et al...... 76 3.3.4 Model IV: Smoothened Coulomb Model ...... 78 3.3.5 Model V: Composite Friction Model...... 78 3.4 Comparison of Sliding Friction Models ...... 80 3.5 Validation and Conclusion...... 85 References for Chapter 3 ...... 89
Chapter 4 Construction of Semi-Analytical Solutions to Spur Gear Dynamics...... 91 4.1 Introduction...... 91 4.2 Problem Formulation ...... 96 4.3 Semi-Analytical Solutions to the SDOF Spur Gear Dynamic Formulation ...... 97 4.3.1 Direct Application of Multi-Term Harmonic Balance (MHBM) ...... 99 ix 4.3.2 Semi-Analytical Solutions Based on One-Term HBM...... 100 4.3.3 Iterative MHBM Algorithm...... 103 4.4 Analysis of Sub-Harmonic Response and Dynamic Instability...... 107 4.5 Semi-Analytical Solutions to 6DOF Spur Gear Dynamic Formulation ...... 110 4.6 Conclusion ...... 116 References for Chapter 4 ...... 122
Chapter 5 Effect of Sliding Friction on the Vibro-Acoustics of Spur Gear System...... 124 5.1 Introduction...... 124 5.2 Source Sub-System Model...... 126 5.3 Structural Path with Friction Contribution ...... 130 5.3.1 Bearing and Housing Models...... 130 5.3.2 Experimental Studies and Validation of Structural Model...... 132 5.3.3 Comparison of Structural Paths in LOA and OLOA Directions ...... 138 5.4 Prediction of Noise Radiation and Contribution of Friction...... 140 5.4.1 Prediction using Rayleigh Integral Technique...... 140 5.4.2 Prediction using Substitute Source Method...... 141 5.4.3 Prediction vs. Noise Measurements...... 144 5.5 Conclusion ...... 148 References for Chapter 5 ...... 150
Chapter 6 Inclusion of Sliding Friction in Helical Gear Dynamics...... 152 6.1 Introduction...... 152 6.2 Problem Formulation ...... 153 6.3 Mesh Forces and Moments with Sliding Friction...... 155 6.4 Shaft and Bearing Models...... 169 6.5 Twelve DOF Helical Gear Pair Model ...... 173 6.6 Role of Sliding Friction Illustrated by an Example ...... 176 6.7 Conclusion ...... 183
x References for Chapter 6 ...... 184
Chapter 7 Analysis of Helical Gear Dynamics using Floquet Theory...... 186 7.1 Introduction...... 186 7.2 Linear Time-Varying Formulation ...... 187 7.3 Analytical Solutions by Floquet Theory...... 199 7.3.1 Response to Initial Conditions...... 201 7.3.2 Forced Periodic Response...... 206 7.4 Conclusion ...... 214 References for Chapter 7 ...... 216
Chapter 8 Conclusion...... 218 8.1 Summary...... 218 8.2 Contributions...... 221 8.3 Future Work...... 223 References for Chapter 8 ...... 224
Bibliography ...... 225
xi LIST OF TABLES
Table Page
2.1 Parameters of Example I: NASA-ART spur gear pair (non-unity ratio)...... 21
2.2 Parameters of Example II-A and II-B: NASA spur gear pair (unity ratio). Gear pair with the perfect involute profile is designated as II-A case and the one with tip relief is designated as II-B case ...... 22
2.3 Averaged coefficient of friction predicted over a range of operating conditions for Example II by using Benedict and Kelly’s empirical equation [2.14]...... 44
5.1 Parameters of the example case: NASA spur gear pair with unity ratio (with long tip relief)...... 128
5.2 Comparison of measured natural frequencies and finite element predictions .... 133
7.1 Relationship between Contact Zones and Contact Regions for the NASA-ART helical gear pair...... 197
xii LIST OF FIGURES
Figure Page
1.1 (a) Schematic of the spur gear contact, where LOA and OLOA represent the line- of-action direction and off line-of-action direction, respectively. (b) Directions of the sliding velocity (V), normal mesh load and friction force in spur gears...... 3
1.2 Comparison of periodic mesh stiffness functions for a spur gear pair ...... 5
1.3 MDOF helical gear dynamic model (proposed in Chapter 6) and its contact mechanics with sliding friction...... 10
1.4 Block diagram for the vibro-acoustics of a simplified geared system with two excitations at the gear mesh (as proposed in Chapter 5)...... 14
2.1 Snap shot of contact pattern (at t = 0) in the spur gear pair of Example I...... 24
2.2 Tooth mesh stiffness functions of Example I calculated by using the FE/CM code (in the “static” mode). (a) Individual and combined stiffness functions. (b) Comparison of the combined stiffness functions...... 26
2.3 Schematic of the bearing-shaft model...... 28
2.4 Normal and friction forces of analytical (MDOF) spur gear system model...... 31
2.5 Validation of the analytical (MDOF) model by using the FE/CM code (in the “dynamic” mode). Here, results for Example I are given in terms of δ ()t and its
spectral contents ∆()f with tc = 2.4 ms and fm = 416.7 Hz. Sub-figures (a-b) are for µ = 0 and (c-d) are for µ = 0.2...... 36
2.6 Validation of the analytical (MDOF) model by using the FE/CM code (in the
“dynamic” mode). Here, results for Example I are given in terms of FtpBx () and
its spectral contents FfpBx () with tc = 2.4 ms and fm = 416.7 Hz. Sub-figures (a-b) are for µ = 0 and (c-d) are for µ = 0.2...... 37
xiii 2.7 Validation of the analytical (MDOF) model by using the FE/CM code (in the
“dynamic” mode). Here, results for Example I are given in terms of FtpBy ( ) and
its spectral contents FfpBy () with tc = 2.4 ms and fm = 416.7 Hz. Sub-figures (a-b) are for µ = 0 and (c-d) are for µ = 0.2 ...... 39
2.8 Effect of µ on δ ()t based on the linear time-varying SDOF model for Example I
at Tp = 2000 lb-in. Here, tc = 1 s ...... 41
2.9 Coefficient of friction µ as a function of the roll angle for Example II, as predicted by using Benedict & Kelley’s empirical equation [2.14]. Here, oil temperature is 104 deg F and T = 500 lb-in...... 45
2.10 Mesh harmonics of the static transmission error (STE) calculated by using the FE/CM code (in the “static” mode) for Example II: (a) gear pair with perfect involute profile (II-A); (b) gear pair with tip relief (II-B) ...... 46
2.11 Tooth stiffness functions of a single mesh tooth pair for Example II: (a) gear pair with perfect involute profile (II-A); (b) gear pair with tip relief (II-B)...... 48
2.12 Combined tooth stiffness functions for Example II: (a) gear pair with perfect involute profile (II-A); (b) gear pair with tip relief (II-B) ...... 49
2.13 Dynamic loads predicted for Example II at 500 lb-in, 4875 RPM and 140 °F with tc = 0.44 ms: (a) Normal loads of gear pair with perfect involute profile (II-A); (b) normal loads of gear par with tip relief (II-B); (c) friction forces of gear pair with perfect involute profile (II-A); (d) friction forces of gear pair with tip relief (II-B)...... 51
2.14 Dynamic shaft displacements predicted for Example II at 500 lb-in, 4875 RPM
and 140 °F with tc = 0.44 ms and fm = 2275 Hz: (a)xp (t ) ; (b) X p (f ) ; (c)ytp ( )
and (d) Yfp ( ) ...... 54
2.15 Dynamic transmission error predicted for Example II at 500 lb-in, 4875 RPM and 140 °F with tc = 0.44 ms and fm = 2275 Hz: (a)δ (t ) ; (b) ∆ (f ) ...... 55
2.16 Sensors inside the NASA gearbox (for Example II-B)...... 57
2.17 Mesh harmonic amplitudes of X p as a function of the mean torque at 140 °F. All
values are normalized with respect to the amplitude of Yp at the first mesh harmonic ...... 59 xiv 2.18 Mesh harmonic amplitudes of yp as a function of the mean torque at 140 °F for
Example II-B. All values are normalized with respect to the amplitude of yp at the first mesh harmonic...... 60
2.19 Predicted dynamic transmission errors (DTE) for Example II over a range of torque at 140 °F: (a) gear pair with perfect involute profile (II-A); (b) with tip relief (II-B). All values are normalized with respect to the amplitude ofδ (II-A) at the first mesh harmonic with 100 lb-in...... 62
2.20 Mesh harmonic amplitudes of xp as a function of temperature at 500 lb-in for
Example II-B. All values are normalized with respect to the amplitude of yp at the first mesh harmonic...... 63
2.21 Mesh harmonic amplitudes of yp as a function of temperature at 500 lb-in for
Example II-B. All values are normalized with respect to the amplitude of yp at the first mesh harmonic...... 64
3.1 (a) Snap shot of contact pattern (at t = 0) in the spur gear pair; (b) MDOF spur gear pair system; here kt ( ) is in the LOA direction...... 71
3.2 (a) Comparison of Model II [3.7] and Model III [3.8] given Tp = 22.6 N-m (200
lb-in) and Ω p = 1000 RPM. (b) Averaged magnitude of the coefficient of friction
predicted as a function of speed using the composite Model V with Tp = 22.6 N-
m (200 lb-in). Here, tc is one mesh cycle...... 81
3.3 Comparison of normalized friction models. Note that curve between 0≤ tt /c < 1
is for pair # 1; and the curve between 1≤ tt /c < 2 is for pair # 0...... 83
3.4 Combined normal load and friction force time histories as predicted using
alternate friction models given Tp = 56.5 N-m (500 lb-in) andΩ p = 4875 RPM...... 84
3.5 Predicted LOA and OLOA displacements using alternate friction models given
Tp = 56.5 N-m (500 lb-in) and Ω p = 4875 RPM ...... 86
3.6 Predicted dynamic transmission error (DTE) using alternate friction models given
Tp = 56.5 N-m (500 lb-in) and Ω p = 4875 RPM ...... 87
xv 3.7 Validation of the normal load and sliding friction force predictions: (a) at Tp =
79.1 N-m (700 lb-in) and Ω p = 800 RPM; (b) at Tp = 79.1 N-m (700 lb-in) and
Ω p = 4000 RPM...... 88
4.1 Realistic mesh stiffness functions of the spur gear pair example (with tip relief)
given Tp = 550 lb-in. (b) Periodic frictional functions...... 92
4.2 (a) Normal (mesh) and friction forces of 6DOF analytical spur gear system model. (b) Snap shot of contact pattern (at t = 0) for the example spur gear pair...... 93
4.3 Semi-analytical vs. numerical solutions for the SDOF model, expressed by Eq.
(4.3), given Tp = 550 lb-in, Ω p = 500 RPM, µ = 0.04. (a) Time domain responses; (b) mesh harmonics in frequency domain...... 102
4.4 Semi-analytical vs. numerical solutions for the SDOF model as a function of pinion speed with µ = 0.04. (a) Mesh order n = 1; (b) n = 2; (c) n = 3; (d) n = 4...... 106
4.5 Normalized determinant of the sub-harmonic matrix Ku as a function of ωnNS / Ω
with µ = 0.04: (a) Tp = 100 lb-in.; (b) Tp = 550 lb-in...... 109
4.6 (a) 6DOF spur gear pair model and its subset of unity gear pair (3DOF model) used to study the natural frequency distribution; (b) Natural frequencies ΩN as a
function of the stiffness ratio KkBm / ...... 111
4.7 Semi-analytical vs. numerical solutions for the 6DOF model as a function of Ω p with µ = 0.04. (a) Mesh order n = 1; (b) n = 2; (c) n = 3; (d) n = 4...... 113
4.8 Semi-analytical vs. numerical solutions of the LOA displacement xp for the 6DOF
model as a function of Ω p with KB/km = 100, µ = 0.04 (a) Mesh order n = 1; (b) n = 2; (c) n = 3; (d) n = 4 ...... 117
4.9 Semi-analytical vs. numerical solutions of the OLOA displacement yp for the
6DOF model as a function of Ω p with KB/km = 100, µ = 0.04. (a) Mesh order n = 1; (b) n = 2; (c) n = 3; (d) n = 4...... 118
xvi 4.10 Semi-analytical vs. numerical solutions of the LOA displacement xp for the 6DOF
model as a function of Ω p with KB/km = 0.37, µ = 0.04. (a) Mesh order n = 1; (b) n = 2; (c) n = 3; (d) n = 4 ...... 120
4.11 Semi-analytical vs. numerical solutions of the OLOA displacement yp for the
6DOF model as a function of Ω p with KB/km = 0.37, µ = 0.04. (a) Mesh order n = 1; (b) n = 2; (c) n = 3; (d) n = 4...... 121
5.1 Block diagram for the vibro-acoustics of a simplified geared system with two excitations at the gear mesh...... 125
5.2 Bearing forces predicted under varying Tp given Ωp = 4875 RPM and 140 °F. (a): LOA bearing force; (b) OLOA bearing force. Key: m is the mesh frequency index ...... 129
5.3 (a) Schematic of NASA gearbox; (b) Finite element model of NASA gearbox with embedded bearing stiffness matrices ...... 131
5.4 Comparison of the gearbox mode shape at 2962 Hz: (a) modal experiment result [5.12]; (b) finite element prediction...... 134
5.5 (a) Experiment used to measure structural transfer functions; (b) comparison of transfer function magnitudes from gear mesh to the sensor on top plate...... 137
5.6 Magnitudes of the combined transfer mobilities in two directions calculated at the sensor location on the top plate...... 139
5.7 Comparison of normal surface velocity magnitudes and substitute source strength vectors under Tp = 500 lb-in and Ωp = 4875 RPM. (a) Line 1: interpolated surface velocity on top plate; (b) Line 2: simplified 2D gearbox model with 15 substitute source points; (c) Line 3: substitute source strengths in complex plane for 2D gearbox. Column 1: mesh frequency index m = 1; Column 2: m = 2; Column 3: m = 3...... 145
5.8 Normalized sound pressure (ref 1.0 Pa) predicted at the microphone 6 inch above the top plate under varying torque Tp given Ωp = 4875 RPM and 140 °F ...... 147
5.9 Overall sound pressure levels (ref: 2e-5 Pa) and their contributions predicted at the microphone 6 in above the top plate under Ωp = 4875 RPM and 140 °F. (a) Tp = 500 lb-in (optimal load for minimum transmission error); (b) Tp = 800 lb-in. 149
xvii 6.1 Schematic of the helical gear pair system...... 156
6.2 Tooth stiffness function calculated using Eq. (6.1) based on the FE/CM code [6.11]...... 157
6.3 Contact zones at the beginning of a mesh cycle. (a) In the helical gear pair; (b) contact zones within contact plane. Key: PP’ is the pitch line; AA’ is the face width W; AD is the length of contact zone Z...... 166
6.4 Predicted tooth stiffness functions...... 169
6.5 Schematic of the bearing-shaft model. Here, the shaft and bearing stiffness elements are assumed to be in series to each other. Only pure rotational or translational stiffness elements are shown. Coupling stiffness terms K , K are xθ y yθx not shown...... 170
6.6 Time and frequency domain responses of translational pinion displacements
uuuxp, yp, zp at Tp = 2000 lb-in and Ω p = 1000 RPM. All displacements are normalized with respect to 39.37 µinch (1 µm)...... 178
6.7 Time and frequency domain responses of pinion bearing forces FSB, xp , FSB, yp and
FSB, zp at Tp = 2000 lb-in and Ω p =1000 RPM. All forces are normalized with respect to 1 lb...... 179
6.8 Time and frequency domain responses of composite displacements δ xyz ,δδ , and velocity δ z at Tp = 2000 lb-in and Ω p =1000 RPM. All motions are normalized with respect to 39.37µinch (1µm) or 39.37µinch (1µm/s) ...... 181
7.1 Schematic of the helical gear pair system...... 189
7.2 Contact zones at the beginning of a mesh cycle within the contact plane. Key: PP’ is the pitch line; AA’ is the face width W; AD is the length of contact zone Z.. 191
7.3 Individual effective stiffness Ke,i(t) along the contact zone, where Tmesh is one mesh cycle...... 196
7.4 Piece-wise effective stiffness function defined in six regions within one mesh cycle with µ = 0.4 ...... 198
xviii 7.5 (a) Effective stiffness and (b) homogeneous responses predictions within two -6 mesh cycles given x0 = 2×10 in., v0 = 20 in./s at Ωp = 1000 RPM ...... 204
7.6 Predictions of damped homogeneous responses within two mesh cycles given x0 = -6 2×10 in., v0 = 20 in./s, µ = 0.2 at Ωp = 1000 RPM ...... 206
7.7 Predictions of (undamped) forced periodic responses within two mesh cycles -6 given x0 = 2×10 in., v0 = 20 in./s, Tp = 2000 lb-in, µ = 0.2 and Ωp = 1000 RPM ...... 210
-6 7.8 Steady state forced periodic responses given x0 = 2×10 in., v0 = 20 in./s, Tp = 2000 lb-in., µ = 0.1 and Ωp = 1000 RPM: (a) DTE vs. time; (b) DTE spectra... 212
7.9 Predicted mesh harmonics of (undamped) forced periodic responses as a function -6 of µ given x0 = 2×10 in., v0 = 20 in./s, Tp = 2000 lb-in and Ωp = 1000 RPM: (a) DTE; (b) slope of DTE ...... 213
7.10 Mapping of eigenvalue κ (absolute value) maxima as a function of the ratio of time-varying mesh frequency fmesh(t) to the system natural frequency fn ...... 215
xix LIST OF SYMBOLS
List of Symbols for Chapter 1
O pinion/gear center location P pitch point t time (s) V contact point speed (in./s) x line-of-action direction y off line-of-action direction z axial direction ε unloaded static transmission error (µin.)
Subscripts 1 pinion 2 gear
Abbreviations ART Advanced Rotorcraft Transmission DOF degree-of-freedom DTE dynamic transmission error EHL elasto-hydrodynamic lubrication LOA line-of-action MDOF multi-degree-of-freedom OLOA off line-of-action SDOF single degree-of-freedom STE static transmission error
List of Symbols for Chapter 2 a, b shaft distance (in) CR surface roughness constant c viscous damping (lb-s/in) E Young’s modulus (psi) F force (lb) f frequency (Hz) I area moment of inertia (in4) i, j indices of gear tooth xx J polar moment of inertia (lb-s2-in) K stiffness matrix (lb/in) K stiffness (lb/in) k tooth mesh stiffness (lb/in) L geometric length (in) N normal contact force (lb) n mesh index R tooth surface roughness (in) r radius (in) T torque (lb-in) t time (s) Ve entraining velocity (in/s) Vs sliding velocity (in/s) Wn normal load per unit length of face width (lb/in) X moment arm (in) XΓ load sharing factor x motion variable along LOA axis (in) Y frequency spectrum of motion along OLOA axis (in) y motion variable along OLOA axis (in) ∆ frequency spectrum of dynamic transmission error (in) δ dynamic transmission error (in) ε static transmission error (in) γ nominal roll angle (rad) λ base pitch (in) µ coefficient of friction 2 υ0 dynamic viscosity (lb-s/in ) θ vibratory angular displacement (rad) σ contact ratio Ω angular speed (rad/s) ζ viscous damping ratio
Subscripts 0, 1… n indices of meshing teeth/modes avg average B bearing b base c (mesh) cycle e effective f friction g gear i index of gear tooth m mesh P pitch point p pinion
xxi S shaft x LOA direction y OLOA direction z axial direction
Superscripts . first derivative with respect to time .. second derivative with respect to time
Abbreviations ART Advanced Rotorcraft Transmission DOF degree-of-freedom DTE dynamic transmission error FE/CM finite element/contact mechanics HPSTC highest point of single tooth contact LPSTC lowest point of single tooth contact LTV linear time-varying LOA line-of-action MDOF multi-degree-of-freedom OLOA off line-of-action SDOF single degree-of-freedom
Operators ceil ceiling function floor floor function mod modulus function sgn sign function
List of Symbols for Chapter 3 b empirical coefficient bH semi-width of Hertzian contact band c viscous damping (lb-s/in) E Young’s modulus (GPa) F force (lb) G dimensionless material parameter H dimensionless central film thickness J polar moment of inertia (lb-s2-in) K stiffness (lb/in) k tooth mesh stiffness (lb/in) L geometric length (in) m mass (lb⋅s2/in) N normal contact force (lb) xxii Ph maximum Hertzian pressure (GPa) r radius (in) S surface roughness (µm) SR slide-to-roll ratio T torque (lb-in) t time (s) U speed parameter V, v velocity (m/s) W load parameter wn normal load per unit length of face width (N/mm) X moment arm (in) x motion variable along LOA axis (in) y motion variable along OLOA axis (in) Z face width (mm) α pressure angle (rad) 2 ηM dynamic viscosity (N-s/mm ) Λ film parameter λ base pitch (in) µ coefficient of friction ν Poisson’s ratio θ vibratory angular displacement (rad) ρ profile radii of curvature (mm) Φ regularizing factor σ contact ratio Ω angular speed (rad/s) ζ viscous damping ratio
Subscripts 0, 1… n indices of meshing teeth avg average B bearing b base C Coulomb friction c (mesh) cycle comp composite e entraining component f friction g gear i index of gear tooth p pinion r rolling component S smoothened friction model s sliding component X Xu and Kahraman model
xxiii x LOA direction y OLOA direction
Superscripts . first derivative with respect to time .. second derivative with respect to time ’ effective value
Abbreviations DOF degree-of-freedom DTE dynamic transmission error EHL elasto-hydrodynamic lubrication HPSTC highest point of single tooth contact LPSTC lowest point of single tooth contact LTV linear time-varying LOA line-of-action MDOF multi-degree-of-freedom OLOA off line-of-action
Operators floor floor function mod modulus function sgn sign function
List of Symbols for Chapter 4
A, B harmonic balance coefficients C damping parameter (lb-s/in) c viscous damping (lb-s/in) D Fourier differentiation matrix E gear constant F harmonic balance matrix F force (lb) f frictional function i, j indices J polar moment of inertia (lb-s2-in) K harmonic balance matrix K stiffness (lb/in) k tooth mesh stiffness (lb/in) L geometric length (in) M (friction) torque (lb-in) or mass (lb⋅s2/in) m mass (lb⋅s2/in) N normal contact force (lb) or harmonic order xxiv n mesh index R base radius (in) r radius (in) S index T (normalized) period t time (s) X moment arm (in) ∆ Fourier coefficient vector δ dynamic transmission error (in) ε static transmission error (in) ϑ angle (rad) λ base pitch (in) µ coefficient of friction θ vibratory angular displacement (rad) σ contact ratio τ dimensionless time υ sub-harmonic index Ω angular speed (rad/s) ω mesh frequency (rad) ζ viscous damping ratio
Subscripts 0, 1… n indices of meshing teeth B bearing b base e effective f friction h super harmonic matrix g gear i index of gear tooth k stiffness coefficient P pitch point p pinion u sub-harmonic matrix x LOA direction y OLOA direction δ dynamic transmission error coefficient
Superscripts . first derivative with respect to time .. second derivative with respect to time ’ first derivative with respect to dimensionless time ” second derivative with respect to dimensionless time ^ differential operator xxv nominal value ~ iterative harmonic balance parameter + pseudo-inverse
Abbreviations DFT discrete Fourier transform DOF degree-of-freedom DTE dynamic transmission error FFT fast Fourier transform LOA line-of-action LTV linear time-varying MDOF multi-degree-of-freedom MHBM multi-term harmonic balance method OLOA off line-of-action SDOF single degree-of-freedom
Operators floor floor function mod modulus function sgn sign function | | matrix determinant
List of Symbols for Chapter 5 e error F force (lb) f frequency (Hz) H transfer function Hv Hankel function I identity matrix i, j indices of gear tooth J polar moment of inertia (lb-s2-in) L geometric length (in) K stiffness (lb/in) k tooth mesh stiffness (lb/in) k(ω) wave number m mass (lb⋅s2/in) N normal contact force (lb) n natural (frequency) P sound pressure (Pa) Q source strength (Pa-in2) xxvi r radius (in) S surface area (in2) T torque (lb-in) t time (s) V velocity (in/s) w weighting function X moment arm (in) x motion variable along LOA axis (in) Y mobility (in/s/lb) y motion variable along OLOA axis (in) α angle (rad) λ base pitch (in) µ coefficient of friction θ vibratory angular displacement (rad) ρ air density (lb⋅s2/in4) σ contact ratio Ω angular speed (rad/s) ω angular frequency (rad) Ξ velocity error matrix ζ viscous damping ratio
Subscripts B bearing b base e effective parameter f friction g gear i index of gear tooth m mean component P path p pinion S shaft or source x LOA direction y OLOA direction
Superscripts . first derivative with respect to time .. second derivative with respect to time ~G complex value vector -1 matrix inverse * complex conjugate Abbreviations
xxvii DOF degree-of-freedom LOA line-of-action MIMO multi-input multi-output OLOA off line-of-action RMS mean square root STE static transmission error
Operators floor floor function mod modulus function sgn sign function | | absolute value
List of Symbols for Chapter 6
E Young’s modulus (psi) e unit vector along axis F force (lb) K tooth mesh stiffness (lb/in) k tooth mesh stiffness density (lb/in2) I area moment of inertia (in4) J polar moment of inertia (lb-s2-in) L length of contact line (in) l variable along contact line (in) M moment (lb-in) m mass (lb⋅s2/in) N normal contact force (lb) N mesh index r radius (in) T torque (lb-in) Tmesh mesh period (s) t time (s) u translational motion (in) W face width (in) v velocity of contact point (in/s) x LOA coordinate of contact point (in) z axial coordinate of contact point (in) β helical angle ∆ deformation of contact point (in) µ coefficient of friction λ base pitch (in) δ dynamic transmission error (in) φ pressure angle (deg) σ contact ratio xxviii Θ (static) angular deflection (rad) θ vibratory angular displacement (rad) Ω angular speed (rad/s) ζ viscous damping ratio
Subscripts 0, 1… n indices of meshing teeth A (shaft) cross section area (in2) b base c contact point g gear h (coordinate) upper limit i index of gear tooth l (coordinate) lower limit P pitch point p pinion S shaft s sliding component V viscous component x LOA direction y OLOA direction z axial direction
Superscripts . first derivative with respect to time .. second derivative with respect to time − mean component <-1> matrix inverse T matrix transverse
Abbreviations ART Advanced Rotorcraft Transmission DOF degree-of-freedom FE/CM finite element/contact mechanics LOA line-of-action LTV linear time-varying MDOF multi-degree-of-freedom OLOA off line-of-action
Operators × cross product ceil ceiling function floor floor function xxix mod modulus function sgn sign function
List of Symbols for Chapter 7
C viscous damping coefficient (lb-s/in) e unit vector along axis F force (lb) f frequency (Hz) G state matrix H transition matrix J polar moment of inertia (lb-s2-in) K tooth mesh stiffness (lb/in) k tooth mesh stiffness (lb/in) L length of contact line (in) M moment (lb-in) m mass (lb⋅s2/in) r radius (in) T torque (lb-in) t time (s) v velocity (in/s) W face width (in) X state space response x LOA coordinate of contact point (in) Z contact zone z axial coordinate of contact point (in) β helical angle δ dynamic transmission error (in) ε static transmission error (in) Φ state transition matrix φ pressure angle (deg) γ basis solution κ eigenvalue λ base pitch (in) µ coefficient of friction Π Wronskian matrix θ vibratory angular displacement (rad) σ contact ratio τ integration variable Ω angular speed (rad/s) ζ viscous damping ratio
Subscripts
xxx 0, 1… n indices of meshing teeth c contact point e effective parameter g gear m mesh (frequency) p pinion x LOA direction y OLOA direction z axial direction
Superscripts . first derivative with respect to time .. second derivative with respect to time − time average
Abbreviations DOF degree-of-freedom DTE dynamic transmission error FE/CM finite element/contact mechanics LOA line-of-action LTV linear time-varying OLOA off line-of-action SDOF single-degree-of-freedom
Operators ceil ceiling function floor floor function LommelS1 Lommel function mod modulus function sgn sign function | | absolute value
List of Symbols for Chapter 8
µ coefficient of friction
Abbreviations DOF degree-of-freedom DTE dynamic transmission error LTV linear time-varying LOA line-of-action MDOF multi-degree-of-freedom OLOA off line-of-action xxxi SDOF single-degree-of-freedom
xxxii CHAPTER 1
INTRODUCTION
1.1 Motivation
Spur and helical gears are widely used in vehicles and mechanical devices to transmit large torques while maintaining a constant input-to-output speed ratio. One remaining challenge for modern gear engineering is the reduction of gear noise in ground and air vehicles such as heavy duty trucks and helicopters. Typically, steady state gear
(whine) noise is generated by several sources [1.1-1.2]. Virtually all of the prior researchers [1.1-1.6] have assumed that the main source is static transmission error (STE), which is defined as the derivation from the ideal (kinematic) tooth profile induced by manufacturing errors and elastic deformations. Accordingly, design engineers tend to reduce STE, via improved manufacturing processes and tooth modifications [1.7]. Yet, at higher torque loads, noise levels are still relatively high though STE might be somewhat minimal (say at the design loads). In other cases, the trend in sound pressure levels does not necessarily match the STE vs. torque curves [1.2]. Typical examples include experimental data on the Advanced Rotorcraft Transmission (ART) gears tested by
1 NASA Glenn and OSU [1.8-1.10]. These suggest that additional vibro-acoustic sources must be considered.
The relative speed between V2 and V1 of two meshing gear teeth (with centers at
O1 and O2), as depicted in Figure 1.1(a), changes direction at the pitch point P during each contact event, thus providing additional periodic excitations normal to the direction of contact, as shown in Figure 1.1(b). Certain unique characteristics of the gear tooth sliding make it a potentially dominant factor, despite the somewhat lower magnitudes of friction force. First, due to the reversal in the direction during meshing action, friction is associated with a large oscillatory component, which causes both higher magnitudes as well as higher bandwidth in dynamic responses. Furthermore, friction is more significant at higher torque and lower speeds. In reality, frictional source mechanism is associated with surface roughness, lubrication regime properties, time-varying friction forces/torques and mesh interface dynamics. These lead to interesting gear dynamic phenomena, such as super-harmonic response, unstable regimes, sub-harmonic resonance and angular modulation [1.11-1.15]. Clearly, the diverse effects of friction can only be analyzed by adopting an intra-disciplinary approach, wherein the principles of meshing kinematics, contact and tribological characteristics, dynamics and noise propagation mechanisms are integrated into a cohesive model.
2
(a)
Sliding Velocity
V2 – V1
Normal Load
Friction Force
(b) Figure 1.1 (a) Schematic of the spur gear contact, where LOA and OLOA represent the line-of-action direction and off line-of-action direction, respectively. (b) Directions of the sliding velocity (V), normal mesh load and friction force in spur gears.
3 Historically, the friction between gear teeth and its cyclic nature have been either ignored or incorporated as an equivalent viscous damping term [1.1-1.2]. Such an approach is clearly inadequate since viscous damping is essentially a passive characteristic and it cannot act as the external excitation to the governing system. Neither does it consider the dynamic effects in the off-line-of-action (OLOA) direction. Hence, there is a definite need for new or improved models that could predict the dynamic and vibro-acoustic responses of a geared system and clarify the role of sliding friction. This is the salient focus of this study.
1.2 Literature Review
In a series of recent articles, Vaishya and Singh [1.13-1.15] have provided an extensive review of prior work. They developed a spur gear pair model with sliding friction and rectangular mesh stiffness by assuming that load is equally shared among all the teeth in contact, as shown in Figure 1.2. They also solved the SDOF system equations in terms of the dynamic transmission error (DTE) by using the Floquet theory and the harmonic balance method [1.13-1.15]. While the assumption of equal load sharing yields simplified expressions and analytically tractable solutions, it may not lead to a realistic model (as shown in Figure 1.2 and then Chapter 2). Houser et al. [1.16-1.17] experimentally demonstrated that the friction forces play a pivotal role in determining the load transmitted to the bearings and housing in the OLOA direction; this effect is more pronounced at higher torque and lower speed conditions.
4 Tooth stiffness (lb/in) stiffness Tooth
Figure 1.2 Comparison of periodic mesh stiffness functions for a spur gear pair. Key: , realistic load sharing (proposed in Chapter 2); , equal load sharing assumed by Vaishya and Singh [1.13].
Velex and Cahouet [1.18] described an iterative procedure to evaluate the effects
of sliding friction, tooth shape deviations and time-varying mesh stiffness in spur and
helical gears and compared simulated bearing forces with measurements. They reported
significant oscillatory bearing forces at lower speeds that are induced by the reversal of
friction excitation with alternating tooth sliding direction. In a subsequent study, Velex and Sainsot [1.19] analytically found that the Coulomb friction should be viewed as a
non-negligible excitation source to error-less spur and helical gear pairs, especially for
translational vibrations and in the case of high contact ratio gears. However, their work
was confined to a study of excitations and the effects of tooth modifications were not
5 considered. Lundvall et al. [1.20] considered profile modifications and manufacturing
errors in a multi-degree-of-freedom (MDOF) spur gear model and examined the effect of
sliding friction on the angular dynamic motions. By utilizing a numerical method, they
reported that the profile modification has less influence on the dynamic transmission
error when frictional effects are included. However, incorporation of the time-varying
sliding friction and the realistic mesh stiffness functions into an analytical (MDOF)
formulation and their dynamic interactions remain unsolved.
In all of the work mentioned above and related literature [1.8-1.20], the sliding friction phenomenon has been typically formulated by assuming the Coulomb formulation with a constant coefficient of friction for modeling convenience. This is partially related to the difficulty associated with the measurement of friction force in a gear mesh. In reality, tribological conditions change continuously due to varying mesh properties, dynamic fluctuations and lubricant film thickness as the gears roll through a full cycle [1.21-1.26]. Thus, coefficient of friction varies instantaneously with the spatial position of each tooth and the direction of friction force changes at the pitch point.
Alternate tribological theories, such as elasto-hydrodynamic lubrication (EHL), boundary lubrication or mixed regime, have been employed to explain the sliding friction under varying operating conditions [1.21-1.23]. For instance, Benedict and Kelley [1.21] proposed an empirical dynamic friction coefficient under mixed lubrication regime based on measurements on a roller test machine. Xu et al. [1.22-1.23] recently proposed yet another friction formula that is obtained by using a non-Newtonian, thermal EHL formulation. Duan and Singh [1.27] developed a smoothened Coulomb model for dry friction in torsional dampers; it could be applied to gears to obtain a smooth transition at
6 the pitch point. Hamrock and Dawson [1.28] suggested an empirical equation to predict
the minimum film thickness for two disks in line contact. They calculated the film
parameter, which could lead to a composite, mixed lubrication model for gears. Rebbechi
et al. [1.8] have successfully used root strains to compute friction force under dynamic
conditions. Recently, Vaishya and Houser [1.9] have shown that quasi-static
measurement of friction force is possible by using the technique of digital filtering to
eliminate the dynamic effects. However, no comprehensive work could be found which
critically evaluates the existing lubrication theories in the framework of an actual gear
mesh. Also, no prior work has incorporated the time-varying coefficient of friction into
MDOF gear dynamics or examined its effect.
Sliding friction at gear teeth also manifests as a noise source, as contended by the
dynamic tests by Houser et al. [1.16]. Borner and Houser [1.17] predicted the dynamic
forces due to friction and qualitatively discussed the radiated sound from the housing.
Most studies on gearbox system dynamics [1.2] have relied on a combination of detailed
finite element, boundary element and semi-analytical methods. Van Roosmalen [1.29] formulated a gearbox model including analytical formulations for the vibration at the gears due to tooth deflections and the vibration transfer through the bearings. Lim et al.
[1.30-1.33] developed a lumped parameter model with a rigid casing and a finite element model with a flexible casing for a simple geared system. However, finite and boundary element methods often require extensive computational time. Over the last four decades, some simplified lumped parameter models have been developed though few have incorporated the torsional and translational motions in both the line-of-action (LOA) and
OLOA directions. Steyer [1.34] examined a single mesh geared system with 6 DOF. By
7 assuming the housing mass is much larger than the gears and shafts, an impedance
mismatch was created with a rigid boundary condition at the bearing location. Thus, the
internal geared system was modeled separately and analytical expressions were presented
for a unity gear pair in terms of the resulting force transmissibility curve. Kartik [1.35] developed a frequency-response based model to predict noise radiation from gearbox housings with a multi-mesh gear set. His work showed that the bearing and mesh stiffness significantly affects the sound pressure in the high frequency range while the casing stiffness controls the response in the range below 4 kHz. However, the transfer function relating the bearing forces to the equivalent force at the housing panel was based on limited experiments. Overall, the above mentioned system models fall short of providing a complete vibro-acoustic model.
1.3 Problem Formulation
1.3.1 Key Research Issues
Governing equations for gear dynamics should lead to a class of damped inhomogeneous periodic differential equations [1.36-1.38] with multiple interacting coefficients [1.13-1.20]. Although similar equations may also be found in a variety of
disciplines such as communication networks [1.36] and electrical circuits [1.39], the gear
friction problems, however, significantly differ from existing models such as the classical
Hill’s equations [1.36] in several ways. First, unlike classic friction problems in most
mechanical systems, the direction of gear friction is normal rather than in the direction of 8 nominal motions. Second, the frictional forces and moments emerge on both sides of the
governing equations as either excitations or periodically-varying parameters. Also, the periodic damping should capture not only the kinematic effects but also the energy dissipation due to sliding friction. Third, the periodic mesh stiffness is not confined to a rectangular wave assumed by Manish and Singh [1.13-1.15], or a simple sinusoid as in the Mathieu’s equation [1.36]. Instead, they should describe realistic, yet continuous,
profiles of Figure 1.2 resulting from a detailed finite element/contact mechanics analysis
[1.40]. Lastly, the stiffness and viscous damping terms incorporate combined (but phase correlated) contributions from all (yet changing) tooth pairs in contact.
Historically, such periodic differential equations are seldom investigated and limited prior research efforts, as reported in the literature review [1.1-1.20], are based largely on numerical integration and the Fast Fourier Transform algorithm. Consequently, there is a clear need for closed form analytical (say by using the Floquet theory) and semi-analytical (say by using the multi-term harmonic balance method) solutions to the dynamic responses of spur and helical gear pairs under the influence of sliding friction.
Recently, Velex and Ajmi [1.41] implemented a harmonic analysis to approximate the dynamic factors in helical gears (based on tooth loads and quasi-static transmission errors). Their work, however, does not describe the multi-dimensional system dynamics or include the frictional effect, which may lead to “multiplicative” terms as described earlier. The parametric friction force excitation may have an influence on the stability of the homogenous system.
Further, for a satisfactory understanding of dynamic behaviors of gears, a higher number of degrees-of-freedom are required for analysis, such as the MDOF helical gear
9 model of Figure 1.3. This is essential for additional phenomena like friction force, torsional-flexural coupling, shaft wobble and axial shuttling, which are yet to be fully understood. Also, to represent practical geared systems, a generalized model is required that incorporates different gear design configurations, lubrication conditions and meshing parameters. Existing solution methodology [1.29-1.35] has to be improved to compute the dynamic response of the entire gearbox, for a combined excitation of transmission error, sliding friction, mean torque and other sources. Subsequently, the relative contribution of various parameters and the resulting noise characteristics need to be understood. This requires an improved source-path-receiver model for the entire gearbox system that incorporates competing noise sources.
ε ()t
Figure 1.3 MDOF helical gear dynamic model (proposed in Chapter 6) and its contact mechanics with sliding friction.
10 1.3.2 Scope, Assumptions and Objectives
Chief goal of this research is to improve the earlier work by Vaishya and Singh
[1.13-1.15] by developing improved mathematical models and proposing new analytical solutions that will enhance our understanding of the influence of friction on gear dynamics and vibro-acoustic behavior. Many dynamic phenomena that emerge due to interactions between parametric variations (time-varying mesh stiffness and viscous damping) and sliding friction will be predicted, along with a better understanding of the relative contributions of transmission error versus sliding friction noise to the gear whine noise. The specific objectives of this study are therefore as follows:
Extend Vaishya and Singh’s work [1.13-1.15] by developing improved MDOF
dynamic system for a spur gear pair that incorporates realistic time-varying mesh
stiffness functions, accurate representations of sliding friction and load sharing
between meshing tooth pairs. (Chapter 2)
Comparatively evaluate alternate sliding friction models [1.21-1.28] and predict
the interfacial friction forces and motions in the OLOA direction. Also, validate
dynamic system models and analytical solutions by comparing predictions to
numerical solutions, the benchmark finite element/contact mechanics code as well
as measurements. (Chapters 2 and 3)
Propose a semi-analytical algorithm based on both single- and multi-term
harmonic balance methods to quickly construct frequency responses of multi-
dimensional spur gear dynamics with sliding friction. This should provide new
insights into the dynamic interactions between parametric excitations. (Chapter 4)
11 Propose a refined source-path-receiver model that characterizes the structural
paths in two directions and develop analytical tools to efficiently predict the
whine noise radiated from gearbox panels and quantify the contribution of sliding
friction to the overall whine noise. Analytical predictions of the structural transfer
function and noise radiation will be compared with measurements. (Chapter 5)
Propose a new three-dimensional formulation for helical gears to characterize the
dynamics associated with the contact plane including the reversal at the pitch line
due to sliding friction. A 12 DOF model will be developed which includes the
rotational and translational motions along the LOA, OLOA and axial directions as
well as the bearing/shaft compliances. (Chapter 6)
Develop improved closed form solutions for the linear time-varying helical gear
system in terms of the dynamic transmission error under the effect of sliding
friction by using the Floquet theory. (Chapter 7)
Scope and assumptions include the following: For the internal spur and helical gear pair sub-systems, the pinion and gear are modeled as rigid disks. The elastic deformations of the shaft and bearings are modeled using lumped elements which are connected to a rigid casing. Also, vibratory angular motions are small in comparison to the mean motion, and the mean load is assumed to be high such that the dynamic load is not sufficient to cause tooth separations [1.42]. If these assumptions are not made, the system model would have implicit non-linearities. Consequently, the position of the line of contact and relative sliding velocity depend only on the nominal angular motions; this leads to a linear time-varying system formulation. Note that different mesh stiffness schemes are assumed for spur and helical gears: For the spur gear analysis, the realistic
12 and continuous mesh stiffness is considered based on an accurate finite element/contact mechanics analysis code [1.40]; Thus, the time-varying stiffness is indeed an effective function which may also include the effect of profile modifications. For the helical gear analysis, however, only those gears with perfect involute profiles are considered and the mesh stiffness per unit length along the contact line (or stiffness density) is assumed to be constant [1.19]. This is equivalent to the equal load sharing assumption by Vishya and
Singh [1.13-1.15]. Such limitation may be further examined in future work.
For the structure-borne whine noise model of the gearbox system, a source-path- receiver model of Figure 1.4 is used. All the assumptions as mentioned above are embedded in the modeling of the internal gear pair sub-system. The unloaded static transmission error and sliding friction are considered as the two main excitations to the system; these are assumed to be most dominant in the LOA and OLOA directions, respectively. Hence, only corresponding structural paths in these two directions are considered by neglecting the moment transfer in the bearing matrices. Also, by assuming the housing mass is much larger than the gears and shafts, an impedance mismatch is created with a rigid boundary condition at the bearing location. Thus, the internal geared system could be modeled separately and its resulting force response provides force excitations to the structural paths. Finally, for the NASA gearbox used as the case study, the box plate is assumed to be the main radiator due to its relatively high mobility as well as the way the gearbox was assembled.
Finally, it is worthwhile to mention that all chapters of this thesis are written in a self-contained manner in terms of formulation, literature review, methods and results.
13 SOURCE Transmission LOA bearing forces error 6 DOF linear-time- Coupling at varying spur gear bearings Sliding pair model + shafts OLOA bearing forces friction
RECEIVER PATH
Sound Radiation Housing Housing structure pressure model velocity model
Figure 1.4 Block diagram for the vibro-acoustics of a simplified geared system with two excitations at the gear mesh (as proposed in Chapter 5).
References for Chapter 1
[1.1] Ozguven, H. N., and Houser, D. R., 1988, “Mathematical Models Used in Gear Dynamics - a Review,” Journal of Sound and Vibration, 121, pp. 383-411.
[1.2] Lim, T. C., and Singh, R., 1989, “A Review of Gear Housing Dynamics and Acoustic Literature,” NASA-Technical Memorandum, 89-C-009.
[1.3] Oswald, F. B., Seybert, A. F., Wu, T. W., and Atherton, W., 1992, “Comparison of Analysis and Experiment for Gearbox Noise,” Proceedings of the International Power Transmission and Gearing Conference, Phoenix, pp. 675-679.
[1.4] Baud, S., and Velex, P., 2002, “Static and Dynamic Tooth Loading in Spur and Helical Geared Systems-Experiments and Model Validation,” American Society of Mechanical Engineers, 124, pp. 334-346.
[1.5] Comparin R. J., and Singh, R., 1990, “An Analytical Study of Automotive Neutral Gear Rattle,” ASME Journal of Mechanical Design, 112, pp. 237-245.
[1.6] Mark, W. D., 1978, “Analysis of the Vibratory Excitation of Gear Systems: Basic Theory,” Journal of Acoustical Society of America, 63(5), pp. 1409-1430.
14
[1.7] Munro, R. G., 1990, “Optimum Profile Relief and Transmission Error in Spur Gears,” Proceedings of IMechE, Cambridge, England, 9-11 Apr., pp. 35-42.
[1.8] Rebbechi, B. and Oswald, F. B., 1991, “Dynamic Measurements of Gear Tooth Friction and Load,” NASA-Technical Memorandum, 103281.
[1.9] Vaishya, M., and Houser, D. R., 1999, “Modeling and Measurement of Sliding Friction for Gear Analysis,” American Gear Manufacturer Association Technical Paper, 99FTMS1, pp. 1-12.
[1.10] Schachinger, T., 2004, “The Effects of Isolated Transmission Error, Force Shuttling, and Frictional Excitations on Gear Noise and Vibration,” MS Thesis, The Ohio State University.
[1.11] Blankenship, G. W., and Singh, R., 1995, “A New Gear Mesh Interface Dynamic Model to Predict Multi-Dimensional Force Coupling and Excitation,” Mechanism and Machine Theory Journal, 30(1), pp. 43-57.
[1.12] Padmanabhan, C., and Singh, R., 1995, “Analysis of Periodically Excited Non- Linear Systems by a Parametric Continuation Technique,” Journal of Sound and Vibration, 184(1), pp. 35-58.
[1.13] Vaishya, M., and Singh, R., 2001, “Analysis of Periodically Varying Gear Mesh Systems with Coulomb Friction Using Floquet Theory,” Journal of Sound and Vibration, 243(3), pp. 525-545.
[1.14] Vaishya, M., and Singh, R., 2001, “Sliding Friction-Induced Non-Linearity and Parametric Effects in Gear Dynamics,” Journal of Sound and Vibration, 248(4), pp. 671- 694.
[1.15] Vaishya, M., and Singh, R., 2003, “Strategies for Modeling Friction in Gear Dynamics,” ASME Journal of Mechanical Design, 125, pp. 383-393.
[1.16] Houser, D. R., Vaishya M., and Sorenson J. D., 2001, “Vibro-Acoustic Effects of Friction in Gears: An Experimental Investigation,” SAE Paper # 2001- 01-1516.
[1.17] Borner, J., and Houser, D. R., 1996, “Friction and Bending Moments as Gear Noise Excitations,” SAE Paper # 961816.
[1.18] Velex, P., and Cahouet, V., 2000, “Experimental and Numerical Investigations on the Influence of Tooth Friction in Spur and Helical Gear Dynamics,” ASME Journal of Mechanical Design, 122(4), pp. 515-522. 15 [1.19] Velex, P., and Sainsot. P, 2002, “An Analytical Study of Tooth Friction Excitations in Spur and Helical Gears,” Mechanism and Machine Theory, 37, pp. 641- 658.
[1.20] Lundvall, O., Strömberg, N., and Klarbring, A., 2004, “A Flexible Multi-body Approach for Frictional Contact in Spur Gears,” Journal of Sound and Vibration, 278(3), pp. 479-499.
[1.21] Benedict, G. H., and Kelley B. W., 1961, “Instantaneous Coefficients of Gear Tooth Friction,” Transactions of the American Society of Lubrication Engineers, 4, pp. 59-70.
[1.22] Xu, H., Kahraman, A., Anderson, N. E., and Maddock, D. G., 2007, “Prediction of Mechanical Efficiency of Parallel-Axis Gear Pairs,” ASME Journal of Mechanical Design, 129 (1), pp. 58-68.
[1.23] Xu, H., 2005, “Development of a Generalized Mechanical Efficiency Prediction Methodology,” PhD dissertation, The Ohio State University.
[1.24] Seireg, A. A., 1998, Friction and Lubrication in Mechanical Design, Marcel Dekker, Inc., New York.
[1.25] Baranov, V. M., Kudryavtsev, E. M., and Sarychev, G. A., 1997, “Modeling of the Parameters of Acoustic Emission under Sliding Friction of Solids,” Wear, 202, pp. 125- 133.
[1.26] Drozdov, Y. N., and Gavrikov, Y. A., 1968, “Friction and Scoring under the Conditions of Simultaneous Rolling and Sliding of Bodies,” Wear, 11, pp. 291-302.
[1.27] Duan, C., and Singh, R., 2005, “Super-Harmonics in a Torsional System with Dry Friction Path Subject to Harmonic Excitation under a Mean Torque”, Journal of Sound and Vibration, 285(2005), pp. 803-834.
[1.28] Hamrock, B. J., and Dowson, D., 1977, “Isothermal Elastohydrodynamic Lubrication of Point Contacts, Part III – Fully Flooded Results,” Journal of Lubrication Technology, 99(2), pp. 264-276.
[1.29] Van Roosmalen, A., 1994, “Design Tools for Low Noise Gear Transmissions,” PhD Dissertation, Eindhoven University of Technology.
[1.30] Lim, T. C., and Singh, R., 1990, “Vibration Transmission Through Rolling Element Bearings. Part I: Bearing Stiffness Formulation,” Journal of Sound and Vibration, 139(2), pp. 179-199.
16 [1.31] Lim, T. C., and Singh, R., 1990, “Vibration Transmission Through Rolling Element Bearings. Part II: System Studies,” Journal of Sound and Vibration, 139(2), pp. 201-225.
[1.32] Lim, T. C., and Singh R., 1991, “Statistical Energy Analysis of a Gearbox with Emphasis on the Bearing Path,” Noise Control Engineering Journal, 37(2), pp. 63-69.
[1.33] Lim, T. C., and Singh, R., 1991, “Vibration Transmission Through Rolling Element Bearings. Part III: Geared Rotor System Studies,” Journal of Sound and Vibration, 151(1), pp. 31-54.
[1.34] Steyer, G., 1987, “Influence of Gear Train Dynamics on Gear Noise,” NOISE- CON 87 proceedings, pp. 53-58.
[1.35] Kartik, V., 2003, “Analytical Prediction of Load Distribution and Transmission Error for Multiple-Mesh Gear-Trains and Dynamic Studies in Gear Noise and Vibration,” MS Thesis, The Ohio State University.
[1.36] Richards, J. A., 1983, Analysis of Periodically Time-varying Systems, New York, Springer.
[1.37] Jordan, D. W., and Smith, P., 2004, Nonlinear Ordinary Differential Equations, 3rd Edition, Oxford University Press.
[1.38] Thomsen, J. J., 2003, Vibrations and Stability, 2nd Edition, Springer.
[1.39] Kenneth S. K., Jacob K. W. and Alberto S-V, 1990, “Steady-State Methods for Simulating Analog and Microwave Circuits,” Kluwer Academic Publishers, Boston.
[1.40] External2D (CALYX software), 2003, “Helical3D User’s Manual,” www.ansol.www, ANSOL Inc., Hilliard, OH.
[1.41] Velex, P., and Ajmi, M., 2007, “Dynamic Tooth Loads and Quasi-Static Transmission Errors in Helical Gears – Approximate Dynamic Factor Formulae,” Mechanism and Machine Theory Journal, 42(11), pp. 1512-1526.
[1.42] Blankenship, G. W., and Kahraman, A., 1995, “Steady State Forced Response of a Mechanical Oscillator with Combined Parametric Excitation and Clearance Type Non- linearity,” Journal of Sound and Vibration, 185(5), pp. 743-765.
17 CHAPTER 2
SPUR GEAR DYNAMICS WITH SLIDING FRICTION AND REALISTIC MESH
STIFFNESS
2.1 Introduction
In a series of recent articles, Vaishya and Singh [2.1-2.3] developed a spur gear pair model with periodic tooth stiffness variations and sliding friction based on the assumption that load is equally shared among all the teeth in contact. Using the simplified rectangular pulse shaped variation in mesh stiffness, they solved the single-degree-of- freedom (SDOF) system equations in terms of the dynamic transmission error (DTE) using the Floquet theory and the harmonic balance method [2.1-2.3]. While the assumption of equal load sharing yields simplified expressions and analytically tractable solutions, it may not lead to a realistic model. This chapter aims to overcome this deficiency by employing realistic time-varying tooth stiffness functions and the sliding friction over a range of operational conditions. New linear time-varying (LTV) formulation will be extended to include multi-degree-of-freedom (MDOF) system dynamics for a spur gear pair.
18 Vaishya and Singh [2.1-2.3] have already provided an extensive review of prior work. In addition, Houser et al. [2.4] experimentally demonstrated that the friction forces play a pivotal role in determining the load transmitted to the bearings and housing in the off-line-of-action (OLOA) direction; this effect is more pronounced at higher torque and lower speed conditions. Velex and Cahouet [2.5] described an iterative procedure to evaluate the effects of sliding friction, tooth shape deviations and time-varying mesh stiffness in spur and helical gears and compared simulated bearing forces with measurements. They reported significant oscillatory bearing forces at lower speeds that are induced by the reversal of friction excitation with alternating tooth sliding direction.
In a subsequent study, Velex and Sainsot [2.6] analytically found that the Coulomb friction should be viewed as a non-negligible excitation source to error-less spur and helical gear pairs, especially for translational vibrations and in the case of high contact ratio gears. However, their work was confined to a study of excitations and the effects of tooth modifications were not considered. Lundvall et al. [2.7] considered profile modifications and manufacturing errors in a MDOF spur gear model and examined the effect of sliding friction on the angular dynamic motions. By utilizing a numerical method, they reported that the profile modification has less influence on the dynamic transmission error when frictional effects are included. Nevertheless, two key questions remain unresolved: How to concurrently incorporate the time-varying sliding friction and the realistic mesh stiffness functions into an analytical (MDOF) formulation? How to quantify dynamic interactions between sliding friction and mesh stiffness terms especially when tip relief is provided to the gears? This chapter will address these issues.
19 2.2 Problem Formulation
2.2.1 Objectives and Assumptions
Chief objective of this chapter is to propose a new method of incorporating the sliding friction and realistic time-varying stiffness into an analytical MDOF spur gear model and to evaluate their interactions. Key assumptions are: (i) pinion and gear are modeled as rigid disks; (ii) shaft-bearings stiffness in the line-of-action (LOA) and
OLOA directions are modeled as lumped elements which are connected to a rigid casing;
(iii) vibratory angular motions are small in comparison to the mean motion; and (iv)
Coulomb friction is assumed with a constant coefficient of friction µ . If assumption (iii) is not made, the system model would have implicit non-linearities. Consequently, the position of the line of contact and relative sliding velocity depend only on the nominal angular motions.
An accurate finite element/contact mechanics (FE/CM) analysis code [2.8] will be employed, in the “static” mode, to compute the mesh stiffness at every time instant under a range of loading conditions. Here, the time-varying stiffness is calculated as an effective function which may also include the effect of profile modifications. The realistic mesh stiffness is then incorporated into the LTV spur gear model with the contributions of sliding friction. The MDOF formulation should describe both the LOA and OLOA dynamics; a simplified SDOF model will also be derived that describes the vibratory motion in the torsional direction. Proposed methods will be illustrated via two spur gear examples (designated as I and II) whose parameters are listed in Table 2.1 and
Table 2.2. The MDOF model of Example I will be validated by using the FE/CM code
20 [2.8] in the “dynamic” mode. Issues related to tip relief will be examined in Example II in the presence of sliding friction. Finally, experimental results of Example II will be used to further validate our method.
Parameter/property Pinion Gear Number of teeth 25 31 Diametral pitch, in-1 8 8 Pressure angle, deg 25 25 Outside diameter, in 3.372 4.117 Root diameter, in 2.811 3.561 Face width, in 1.250 1.250 Tooth thickness, in 0.196 0.196 Gear mass, lb⋅s2⋅in-1 6.72E-03 1.04E-02 Polar moment of inertia, lb⋅s2⋅in 8.48E-03 2.00E-02 Bearing stiffness (LOA and OLOA), lb/in 20E6 Center distance, in 3.5 Profile contact ratio 1.43 Elastic modulus, psi 30E6 Density, lb⋅s-2⋅in-4 7.30E-04 Poisson’s ratio 0.3
Table 2.1 Parameters of Example I: NASA-ART spur gear pair (non-unity ratio)
21
Parameter/property Pinion/Gear
Number of teeth 28
Diametral pitch, in-1 8
Pressure angle, deg 20
Outside diameter, in 3.738
Root diameter, in 3.139
Face width, in 0.25
Tooth thickness, in 0.191
Roll angle where the tip modification 24.5 starts (for II-B), deg
Straight tip modification (for II-B), in 7E-04
Center distance, in 3.5
Profile contact ratio 1.63
Elastic modulus, psi 30E6
Density, lb⋅s-2⋅in-4 7.30E-04
Poisson’s ratio 0.3
Range of temperatures, °F 104, 122, 140, 158, 176
Range of input torques, lb⋅in 500, 600, 700, 800, 900
Table 2.2 Parameters of Example II-A and II-B: NASA spur gear pair (unity ratio). Gear pair with the perfect involute profile is designated as II-A case and the one with tip relief is designated as II-B case
22 2.2.2 Timing of Key Meshing Events
Analytical formulations for a spur gear pair are derived via Example I (NASA-
ART spur gear pair) with parameters of Table 2.1. For a generic spur gear pair with non- integer contact ratio σ , n = ceil(σ ) meshing tooth pairs need to be considered, where the
“ceil” function rounds the σ element to the nearest integer towards a higher value.
Consequently, two meshing tooth pairs need to be modeled for Example I (σ = 1.43).
First, transitions in key meshing events within a mesh cycle need to be determined
from the undeformed gear geometry for the construction of the stiffness function. Figure
2.1 is a snapshot for Example I at the beginning of the mesh cycle (t = 0). At that time,
pair #1 (defined as the tooth pair rolling along line AC) just comes into mesh at point A
and pair #0 (defined as the tooth pair rolling along line CD) is in contact at point C,
which is the highest point of single tooth contact (HPSTC). As the gears roll, when pair
#1 approaches the lowest point of single tooth contact (LPSTC) of point B at t = tB, pair
#0 leaves contact. At t =tP , pair #1 passes through the pitch point P, and the relative
sliding velocity of the pinion with respect to the gear is reversed, resulting in a reversal of
the friction force. This should provide an impulse excitation to the system. Finally, pair
#1 goes through point C at t = tc, completing one mesh cycle (tc). These key events are
defined below, where Ω p is the nominal pinion speed, rbp is the base radius of the pinion,
length LAC is equal to one base pitch λ.
t λ tb LAB p LAP tc = , = , = . (2.1) Ω pbpr tc λ tc λ
23
Ωg
Ω p
Figure 2.1 Snap shot of contact pattern (at t = 0) in the spur gear pair of Example I.
2.2.3 Calculation of Realistic Time-Varying Tooth Stiffness Functions
The realistic time-varying stiffness functions are calculated using a FE/CM code,
External2D [2.8]. An input torque Tp is applied to the pinion rotating at Ω p , and the
mean braking torque Tg on the gear and its angular velocity Ωg obey the basic gear
kinematics. Superposed on the nominal motions are oscillatory components denoted as
θ p and θg for the pinion and gear, respectively. The normal contact forces N0(t), N1(t)
and pinion deflection θ p ()t are then computed by performing a static analysis using
FE/CM software [2.8]. The stiffness function of the ith meshing tooth pair for a generic
24 spur gear pair is given by Eq. (2.2), where the “floor” function rounds the contact ratio σ to the nearest integer towards a lower value, i.e. floor(σ ) = 1 for Example I.
Nti () kti ( )=== , i 0, 1, ... , n floor(σ ). (2.2) rtbpθ p ()
The stiffness function k(t) for a single tooth pair rolling through the entire meshing process is obtained by following the contact tooth pair for n = ceil(σ ) number
of mesh cycles. Due to the periodicity of the system, expanded stiffness function kti ( ) of the ith meshing tooth pair is calculated at any time instant t as:
kticc( )=−+ k[] ( n it ) mod( tt , ) , i = 0, 1, ... , n = floor(σ ). (2.3)
Here, “mod” is the modulus function defined as:
mod(xy , )=−⋅ x y floor( x / y ), if y ≠ 0. (2.4)
For Example I, calculated kt0 (), kt1() functions and their combined stiffness are
shown in Figure 2.2(a). Note that kt0 () and kt1() are, in fact, different portions of kt() as described in Eq. (2.3). Figure 2.2(b) compares the continuous kt() of the realistic load sharing model against the rectangular pulse shaped discontinuous kt() based on the equal load sharing formulation proposed earlier by Vaishya and Singh [2.1-2.3].
25
Figure 2.2 Tooth mesh stiffness functions of Example I calculated by using the FE/CM code (in the “static” mode). (a) Individual and combined stiffness functions. Key: , total stiffness; stiffness of pair #0; , stiffness of pair #1. (b) Comparison of the combined stiffness functions. Key: , realistic load sharing; , equal load sharing as assumed by Vaishya and Singh [2.1-2.3].
26 2.3 Analytical Multi-Degree-of-Freedom Dynamic Model
2.3.1 Shaft and Bearing Stiffness Models
Next, we develop a generic spur gear pair model with 6 DOFs including rotational
motions (θ p and θg ), LOA translations ( xp and xg ) and OLOA translations ( y p and yg ).
The governing equations are derived in the subsequent sections. First, a simplified shaft model, as shown in Figure 2.3, is developed based on the Euler’s beam theory [2.9].
Corresponding to the 6 DOFs mentioned above, only the diagonal term in the shaft stiffness matrix needs to be determined as follows, where E is the Young’s modulus,
4 Ir= π s /4 is the area moment of inertia for the shaft, and a and b are the distances from pinion/gear to the bearings.
ab+ 2 K ==KEIabab3 ⎡() −+⎤ , K = 0 . (2.5) Sx Sy ab33⎣ ⎦ Sθz
The rolling element bearings are modeled using the bearing stiffness matrix K Bm formulation (of dimension 6) as proposed by Lim and Singh [2.10, 2.11]. Assume that each shaft is supported by two identical axially pre-loaded ball bearings with a mean
axial displacement; the mean driving load Tm generates a mean radial force Fxm in the
LOA direction and a moment M ym around the OLOA direction. The time-varying friction force and torque are not included in the mean loads.
27
KBx /2
a
KBy /2 b KBx /2 x
θ z KBy /2 z y
Figure 2.3 Schematic of the bearing-shaft model.
Corresponding to the 6 DOFs considered in our spur gear model, only two
significant coefficients, KBxx and K Byy , are considered for K Bm [2.10, 2.11]. The
combined bearing-shaft stiffness ( KBx and K By in the LOA and OLOA directions) are derived by assuming that the bearing and shaft stiffness elements act in parallel.
2.3.2 Dynamic Mesh and Friction Forces
Figure 2.4 shows the mean torque and internal reaction forces acting on the pinion for Example I. For the sake of clarity, forces on the gear are not shown, which are equal in magnitude but opposite in direction to the pinion forces. Based on the Coulomb 28 friction law, the magnitude of friction force ( Ff ) is proportional to the nominal tooth
load (N) as Ff = µN where µ is constant. The direction of Ff is determined by the calculation of nominal relative sliding velocity, which results in the LTV system
formulation. Denote Xtpi () as the moment arm on the pinion for the friction force acting on the ith meshing tooth pair
Xtpi( )=+−+ L XA ( ni )λ mod( Ω p rt bp ,λσ ), i = 0, 1, ... , n = floor( ). (2.6)
The corresponding moment arm for the friction force on the gear is
Xtgi( )=+−Ω L YC iλ mod( g rt bg ,λσ ), i = 0, 1, ... , n = floor( ). (2.7)
Assume time-varying mesh (viscous) damping coefficient and relate it to kti ( ) by
22 a time-invariant damping ratio ζ m as follows, where Jepgpbggbp=+JJ/ ( Jr Jr)
ctimiie( )=⋅== 2ζ kt ( ) J , i 0, 1, ... , n floor(σ ). (2.8)
The normal forces acting on the pinion are
⎡⎤ Ntpi()== Nt gi () ktr i ()⎣⎦ bpθθε p () t − r bg g () t −+−+ p () t xt p () xt g () (2.9) ⎡⎤ ctibppbggp()⎣⎦ rθ () t−−+−= rθε () t () t x p () t x g () t , i 0, 1, ... , n = floor( σ ).
29 Here ε P (t ) is the profile error component of the static transmission error (STE), and xp(t) and xg(t) denote the translational bearing displacements of pinion and gear, respectively.
For a generic spur gear pair whose jth meshing pair passes through the pitch point within the mesh cycle, the friction forces in the ith meshing pair are derived as follows
⎧ µNtpi ( ), i= 0, 1, ... , j− 1, ⎪ ⎡⎤ Fpfi()tNt=Ω+−−=⎨ µλλ pi ()sgnmod(⎣⎦ p rtniLij bp , ) ( ) AP , , (2.10a) ⎪ ⎩⎪−=+=µNtpi ( ), i jj , 1, ... , n floor(σ ),
⎧ µNtgi ( ), i= 0, 1, ... , j− 1, ⎪ ⎡⎤ Fgfi(tNt )=Ω+−−=⎨ µλλ gi ( )sgn⎣⎦ mod( g rtniLij bg , ) ( ) AP , , (2.10b) ⎪ ⎩⎪−=+=µNtgi ( ), i jj , 1, ... , n floor(σ ).
Consequently, the friction forces for Example I of Figure 2.4 are given as:
Fpf00()tNt= µ p (), (2.10c)
⎡ ⎤ Ftpf11()=Ω−µλ Nt p ()sgnmod(⎣ p rt bp , ) L AP ⎦ , (2.10d)
Fgf00()tNt= µ g (), (2.10e)
⎡ ⎤ Ftgfg11()=Ω−µλ Nt ()sgnmod(⎣ gbgAP rt , ) L⎦ . (2.10f)
30
Figure 2.4 Normal and friction forces of analytical (MDOF) spur gear system model.
31 2.3.3 MDOF Model
The governing equations for the torsional DOFs are
nn==floor(σσ ) floor( ) J ppθ ()tT=+ p∑∑ XtFt pi () pfi () − rNt bppi (), (2.11) ii==00
nn==floor(σσ ) floor( ) J gθ g()tT=− g +∑∑ XtFt gi () gfi () + rNt bg gi () (2.12) ii==00.
The governing equations of the translational DOFs in the LOA direction are
n=floor(σ ) mxp p() t+++= 2ζ pBx K pBx mx p p () t K pBx x g () t∑ N pi () t 0 , (2.13) i=0
n=floor(σ ) mxg g() t+++= 2ζ gBx K gBx mx g g () t K gBx x g () t∑ N gi () t 0. (2.14) i=0
Here, K pBx and K gBx are the effective shaft-bearing stiffness in the LOA direction, and
ζ pBx and ζ gBx are their damping ratios. Similarly, the governing equations of the translational DOFs in the OLOA direction are
n=floor(σ ) myp p() t++−= 2ζ pB y K pBy my p p () t K pBy y p () t∑ F pfi () t 0, (2.15) i=0
n=floor(σ ) myg g() t++−= 2ζ gB y K gBy my g g () t K gBy y g () t∑ F gfi () t 0. (2.16) i=0 32 The composite DTE, which is the relative dynamic displacement of pinion and
gear along the LOA direction, is defined as
δθ()tr=−+−bp p () tr bg θ g () txtxt p () g (). (2.17)
Finally, the dynamic bearing forces are as:
FpBx()t=− K pBx xt p () − 2ζ pBx K pBx mxt p p (), (2.18a)
FpBy()t=− K pBy yt p () − 2ζ pBb K pBy myt p p (), (2.18b)
FgBx()t=− K gBx xt g () − 2ζ gBb K gBx mxt g g (), (2.18c)
FgBy()t=− K gBy yt g () − 2ζ g Bb K gBy myt g g (). (2.18d)
2.4 Analytical SDOF Torsional Model
When only the torsional DOFs of the spur gear pair are of interest, a simplified but equivalent SDOF model can be derived by assuming that the shaft-bearings stiffness is much higher than the mesh stiffness. After eliminating θp(t) and θg(t) in terms of the
DTE, δθ()tr=−bp p () tr bg θ g () t, the governing SDOF model is obtained for a generic spur
gear pair whose jth meshing pair passes through the pitch point within the mesh cycle:
33 n=floor(σ ) Jt ( )⎡⎤ cttktt ( ) ( ) ( ) ( ) eiiδδδ+++∑ ⎣⎦ i=0 ⎧⎫⎡⎤ sgn⎣⎦ mod(Ω+−−⋅pbprt ,λλ ) ( n j ) L AP n=floor(σ ) ⎪⎪ µ ⎨⎬⎡ XtJrXtJr()+ () ⎤ = ∑ pj g bp gj p bg i=0 ⎪⎪⎡⎤ctii()δδ () t+⋅ kt () () t ⎢ ⎥ ⎣⎦Jr22+ Jr ⎩⎭⎪⎪⎣⎢ pbg gbp ⎦⎥ . (2.19) T n=floor(σ ) e +++⎡⎤ct ()εε () t kt () () t 22∑ ⎣⎦ip ip Jrpbg+ Jr gbp i=0 ⎧⎫⎡⎤ sgn⎣⎦ mod(Ω+−−⋅pbprt ,λλ ) ( n j ) L AP n=floor(σ ) ⎪⎪ µ ∑ ⎨⎬⎡ XtJrXtJrpj() g bp+ gj () p bg ⎤ i=0 ⎪⎪⎡⎤ct()εε () t+⋅ kt () () t ⎢ ⎥ ⎣⎦ip ip Jr22+ Jr ⎩⎭⎪⎪⎣⎢ pbg gbp ⎦⎥
Here the effective polar moment of inertia Je is consistent with that defined in Eq. (2.8)
and the effective torque is TTJrTJrepgbpgpbg= ⋅⋅+⋅⋅. The dynamic response δ ()t is
controlled by three excitations: (i) time-varying Te , (ii) ε p ()t and its derivativeεp ()t and
(iii) sliding friction. For Example I, the governing Eq. (2.19) could be simplified as
Jteδδδ ( )++[ ctctt10 ( ) ( )] ( ) ++[ ktktt 10 ( ) ( )] ( ) + ⎡⎤XtJrXtJr()+ () µδ⎡⎤ct() () t+⋅Ω−+ kt () δ () tp11 g bp g p bg sgnmod(⎡⎤ rt , λ ) L ⎣⎦11⎢⎥22 ⎣⎦pbp AP ⎣⎦⎢⎥Jrpbg+ Jr gbp
()XtJrXtJrp00() g bp+ g () p bg µ ⎡⎤c ()ttkttδδ ()+ () () ⎣⎦0 0 22 Jrpbg+ Jr gbp T (2.20) e =+++++22[][]ct10 ( ) ct ( )εεpp ( t ) kt 10 ( ) k ( t ) ( t ) Jrpbg+ Jr gbp ⎡ XtJrXtJr()+ () ⎤ µε⎡⎤ct() () t+ kt () ε () t p11 g bp g p bg ⋅sgn⎡ mod(Ω−+rt ,λ ) L ⎤ ⎣⎦11pp⎢⎥22 ⎣ pbp AP⎦ ⎢⎥⎣ Jrpbg+ Jr gbp ⎦
()XtJrXtJrpgbpgpbg00()+ () µε⎡⎤ct ( ) ( t )+ kt ( ) ε ( t ) ⎣⎦00pp 22 Jrpbg+ Jr gbp
34 2.5 Effect of Sliding Friction in Example I
2.5.1 Validation of Example I Model using the FE/CM Code
The governing equations of either SDOF or MDOF system models are numerically integrated by using a 4th-5th order Runge-Kutta algorithm with fixed time
step. The ε p ()t and εp ()t components are neglected, i.e. no manufacturing errors other than specified profile modifications are considered. Concurrently, the dynamic responses are independently calculated by running the FE/CM code [2.8] using the Newmark method. Predicted and computed results are compared with good correlations in terms of the DTE, and LOA and OLOA forces, as shown in Figure 2.5 to Figure 2.7. Note that time domain comparisons include both transient and steady state responses but the frequency domain results report only the steady state responses. Figure 2.5 shows that the sliding friction introduces additional DTE oscillations when the contact teeth pass through the pitch point. Figure 2.6 illustrates that the sliding friction enhances the dynamic bearing forces in the LOA direction, especially at the second mesh harmonic.
This is because the moments associated with Fpfi ()t and Fgfi ()t are coupled with the moments of Npi(t) and Ngi(t).
35
(t) (in) (t) (f) (in) (t) (in) (t) (f) (in) (f)
Figure 2.5 Validation of the analytical (MDOF) model by using the FE/CM code (in the “dynamic” mode). Here, results for Example I are given in terms of δ ()t and its spectral contents ∆()f with tc = 2.4 ms and fm = 416.7 Hz. Sub-figures (a-b) are for µ = 0 and (c-d) are for µ = 0.2. Key: , Analytical (MDOF) model; , FE/CM code (in t domain); o, FE/CM code (in f domain).
36
x 103 x 103
0.2 3
2
0.1 1
0
0 0 0.5 1 1.5 2 2.5 0 2 4 6 8 10 (a)Normalized time t / t (b) c Mesh order n x 103 x 103
0.2 3
2
0.1 1
0
0 0 0.5 1 1.5 2 2.5 0 2 4 6 8 10 (c) (d) Normalized time t / t c Mesh order n
Figure 2.6 Validation of the analytical (MDOF) model by using the FE/CM code (in the
“dynamic” mode). Here, results for Example I are given in terms of FpBx ()t and its spectral contents FpBx ()f with tc = 2.4 ms and fm = 416.7 Hz. Sub-figures (a-b) are for µ = 0 and (c-d) are for µ = 0.2. Key: , Analytical (MDOF) model; , FE/CM code (in t domain); o, FE/CM code (in f domain).
37 Further, the normal loads mainly excite the vibration in the LOA direction, as illustrated by Eqs. (2.9), (2.11) and (2.13). The scales of the bearing forces of Figure
2.7(a-b) for µ = 0 case are the same as those of Figure 2.7(c-d) for the sake of comparison.
The bearing forces predicted by the MDOF model for µ = 0 case approach zero (within the numerical error range). This is consistent with the mathematical description of Eqs.
(2.15-2.16). Larger deviations at this point are observed in Figure 2.7(a-b) for the FE/CM analysis. Figure 2.7shows that the OLOA dynamics are more significantly influenced by the sliding friction when compared with the LOA results of Figure 2.6. In order to accurately predict the higher mesh harmonics, refined time steps (say more than 100 increments per mesh cycle) are needed. Consequently, the FE/CM analysis tends to generate an extremely large data file that demands significant computing time and post- processing work. Meanwhile, the lumped model allows much finer time resolution while being computationally more efficient (by at least two orders of magnitude when compared with the FE/CM). Hence, the lumped model could be effectively used to conduct parametric design studies.
38
3 x 103 x 10 0.2
5
0 (f) (lb) (f) (t) (lb) (t) 0.1 pBy pBy F F -5
-10 0 0 0.5 1 1.5 2 2.5 0 2 4 6 8 10 (a) (b) Normalized time t / tc Mesh order n
x 103 x 103 0.2 5
0 (t) (lb) (f) (lb) 0.1 pBy pBy F
F -5
-10
0 0 0.5 1 1.5 2 2.5 0 2 4 6 8 10 (c) Normalized time t / t (d) Mesh order n c
Figure 2.7 Validation of the analytical (MDOF) model by using the FE/CM code (in the
“dynamic” mode). Here, results for Example I are given in terms of FpBy ()t and its spectral contents FpBy ()f with tc = 2.4 ms and fm = 416.7 Hz. Sub-figures (a-b) are for µ = 0 and (c-d) are for µ = 0.2. Key: , Analytical (MDOF) model; , FE/CM code (in t domain); o, FE/CM code (in f domain).
39 2.5.2 Effect of Sliding Friction
Figure 2.8 shows the calculated DTE without any friction is almost identical to
the STE at a very low speed ( Ω p = 2.4 rpm). However, the sliding friction changes the
shape of the DTE curve. During the time intervaltt∈ [0,P ], the friction torque on the pinion opposes the normal load torque as shown in Figure 2.4, resulting in a higher value of the normal load that is needed to maintain the static equilibrium. Also, friction increases the peak-peak value of the DTE as compared with the STE. For the remainder
of the mesh cycle ttt∈ [Pc , ], friction torque acts in the same direction as the normal load torque. Thus a small value of normal load is sufficient to maintain the static equilibrium.
Detailed parametric studies show that the amplitude of second mesh harmonic increases with the effect of sliding friction.
40
Figure 2.8 Effect of µ on δ ()t based on the linear time-varying SDOF model for
Example I at Tp = 2000 lb-in. Here, tc = 1 s. Key: , µ = 0; , µ = 0.1.
41 2.5.3 MDOF System Resonances
6 For Example I, the nominal bearing stiffness KKKKpBx====× pBy gBx gBy 20 10
lb/in are much higher than the averaged mesh stiffness km . The couplings between the rotational and translational DOFs in the LOA direction are examined by using a simplified 3 DOF model as suggested by Kahraman and Singh [2.12]. Note that the DTE
is defined here as δ =−rrbpθθ xp bg xg , and the undamped equations of motion are
⎡⎤ ⎡⎤⎧⎫ kk− k ⎧⎫ me 00δδ⎢⎥mm m ⎧⎫ 0 ⎢⎥⎪⎪⎢⎥ ⎪⎪⎪⎪ 00mxkkKkxppmmpBxmp⎨⎬ ++−() ⎨⎬⎨⎬ = 0. (2.21) ⎢⎥⎢⎥ ⎢⎥00mx⎪⎪ ⎪⎪⎪⎪ x 0 ⎣⎦gg⎩⎭⎢⎥−−kkkK + ⎩⎭ g⎩⎭ ⎣⎦mmmgBx()
22 Here, the effective mass is defined as mJJrJrJepggppg=+/ ( ) . The eigensolutions of
Eq. (2.21) yield three natural frequencies: Two coupled transverse-torsional modes ( f1
and f3 ) and one purely transverse mode ( f2 ); numerical values are: f1 = 5,130 Hz, f2 =
8,473 Hz and f3 = 11,780 Hz. Predictions of Eq. (2.21) match well with the numerical simulations using the formulations of section 2.3 (though these results are not shown here). A comparative study verifies that one natural frequency of the MDOF model shifts away from that of the SDOF model (6,716 Hz) due to the torsional-translational coupling effects. In the OLOA direction, simulation shows that only one resonance is present at
1 f = Km/ = 9,748 Hz, which is dictated by the bearing-shaft stiffness. pBy2π pBy p
42 2.6 Effect of Sliding Friction in Example II
Next, the proposed model is applied to Example II with the parameters of Table
2.2. The chief goal is to examine the effects of tip relief and sliding friction. Further, analogous experiments were conducted at the NASA Glenn Research Center Gear Noise
Rig [2.13]. Comparisons with measurements will be given in section 2.7.
2.6.1 Empirical Coefficient of Friction
The coefficient of friction varies as the gears travel through mesh, due to constantly changing lubrication conditions between the contact teeth. An empirical equation for the prediction of the dynamic friction variable, µ , under mixed lubrication has been suggested by Benedict and Kelley [2.14] based on a curve-fit of friction measurements on a roller test machine. Rebbechi et al. [2.15] verified this formulation by measuring the dynamic friction forces on the teeth of a spur gear pair. Their measurements seem to be in good agreement with the Benedict and Kelley equation except at the meshing positions close to the pitch point. This empirical equation, when
modified to account for the average gear tooth surface roughness ( Ravg ), is
8 ⎛⎞3.17× 10XWΓ (γ ) n 44.5 µγ( )= 0.0127CRavg log10 ⎜⎟2 , CRavg = . (2.22a,b) ⎝⎠νγosVV() e () γ 44.5− Ravg
43 where CRavg is the surface roughness constant, Wn is the normal load per unit length of
face width, and υo is the dynamic viscosity of the lubricant. Here Vs (γ ) is the sliding velocity, defined as the difference in the tangential velocities of the pinion and gear, and
Ve ()γ is the entraining velocity, defined as the addition of the tangential velocities, for
roll angle γ along the LOA. Further, Ravg in our case was measured with a profilometer
using a standard method [2.13]. Lastly, X Γ (γ ) is the load sharing factor as a function of roll angle, and it was assumed based on the ideal profile of smooth meshing gears. Figure
2.9 shows µ as a function of roll angle calculated using Eq. (2.22). Since µ was assumed to be a constant earlier, an averaged value is found by taking an average over the roll angles between 19.8 and 21.8 degrees. Table 2.3 lists the µ values that were
computed at each mean torque and oil temperature for Example II (with Ravg = 0.132
µ m).
Temperature (°F) Torque (lb-in) 500 600 700 800 900 104 0.032 0.033 0.034 0.035 0.036 122 0.034 0.036 0.037 0.037 0.038 140 0.036 0.037 0.038 0.039 0.040 158 0.038 0.040 0.041 0.041 0.042 176 0.040 0.041 0.042 0.043 0.044
Table 2.3 Averaged coefficient of friction µ predicted over a range of operating conditions for Example II by using Benedict and Kelly’s empirical equation [2.14] 44
Figure 2.9 Coefficient of friction µ as a function of the roll angle for Example II, as predicted by using Benedict & Kelley’s empirical equation [2.14]. Here, oil temperature is 104 deg F and T = 500 lb-in. Key P: Pitch point at 20.85 deg.
2.6.2 Effect of Tip Relief on STE and k(t)
The STE is calculated as a function of mean torque for both the perfect involute gear pair (designated as II-A) and then one with tip relief (designated as II-B) using
FE/CM code. Figure 2.10 compares the amplitudes of STE spectra at mesh harmonics for both cases. (In this and following figures, predictions are shown as continuous lines for the sake of clarity though they are calculated only at discrete torque points.)
45 x 10-4 2.5
2
1.5
1
Transmission error (in) error Transmission 0.5
0 100 300 500 700 900 (a) Torque (lb-in) Transmission error (in) error Transmission
Figure 2.10 Mesh harmonics of the static transmission error (STE) calculated by using the FE/CM code (in the “static” mode) for Example II: (a) gear pair with perfect involute profile (II-A); (b) gear pair with tip relief (II-B). Key: , n = 1; , n = 2; , n = 3. 46 The first two mesh harmonics are most significantly affected by the tip relief and they are minimal at the “optimal” mean torque around 500 lb-in. For both the II-A and II-
B cases, typical kt() functions of a single meshing tooth over two complete mesh cycles are calculated using Eq. (2.2) for various mean torques, as shown in Figure 2.1. Note that kt() is defined as the effective stiffness since it incorporates the effect of profile modification such as the linear tip relief (II-B). Observe that although the maximum stiffness remains the same, application of the tip relief significantly changes the stiffness profile. For the perfect involute profile (II-A), steep slopes are observed in the vicinities near the single or two teeth contact regimes, and a smooth transition is observed in between these steep regimes. Also, kt() is found to be insensitive to a variation in the mean torque. However, with tip relief, an almost constant slope is found throughout the transition profile between single and two teeth contact regimes. Moreover, a smaller profile contact ratio (around 1.1 at 100 lb-in) is observed for the tip relief case when compared with around 1.6 (at all loads) for the perfect involute pair. The realistic kt() function is then incorporated into the lumped MDOF dynamic model.
Figure 2.12 shows the combined kt() with contributions of both meshing tooth pairs over two mesh cycles for Example II. Observe that the profile of case II-A is insensitive to a variation in the mean torque, but the profile of case II-B shows a minimum around 500 lb-in. Frequency domain analysis reveals that the first two mesh harmonics are most significantly affected by the linear tip modification. Overall, it is evident that significant changes take place in the STE, tooth load distribution and mesh stiffness function due to the profile modification (tip relief), which may be explained by an avoidance of the corner contact at an “optimized” mean torque. 47
Figure 2.11 Tooth stiffness functions of a single mesh tooth pair for Example II: (a) gear pair with perfect involute profile (II-A); (b) gear pair with tip relief (II-B). Key: , 100 lb-in; , 500 lb-in; , 900 lb-in. 48
Figure 2.12 Combined tooth stiffness functions for Example II: (a) gear pair with perfect involute profile (II-A); (b) gear pair with tip relief (II-B). Key: , 100 lb-in; , 500 lb-in; , 900 lb-in. 49 2.6.3 Phase Relationship between Normal Load and Friction Force Excitations
Using the 6DOF spur gear model with parameters consistent with the experimental conditions, dynamic studies are conducted for Example II. First, a mean torque of 500 lb-in is used corresponding to the “optimal” case with minimal STE.
Equations (2.13-2.16) show that the normal loads ∑ Ni and friction forces ∑ Ffi excite the LOA and OLOA dynamics, respectively. The force profile of a single tooth pair undergoing the entire meshing process is obtained by tooth pairs #0 and #1 for two continuous meshing cycles as shown in Figure 2.13(a-b) and (c-d) for II-A and II-B cases respectively. Observe that the peak-to-peak magnitude of combined pinion normal load
∑ N pi is minimized for the tip relief gear due to reduced STE at 500 lb-in. However, the combine pinion friction force ∑ Ffpi with tip relief has a higher peak-to-peak magnitude when compared with the perfect involute gear. This implies that the tip relief amplifies
∑ Ffi in the OLOA direction while minimizing ∑ Ni in the LOA direction. Such
contradictory effects are examined next using the phase relationship between N pi and
Ffpi .
50
Figure 2.13 Dynamic loads predicted for Example II at 500 lb-in, 4875 RPM and 140 °F with tc = 0.44 ms: (a) Normal loads of gear pair with perfect involute profile (II-A); (b) normal loads of gear par with tip relief (II-B); (c) friction forces of gear pair with perfect involute profile (II-A); (d) friction forces of gear pair with tip relief (II-B). Key: , combined; , tooth pair #0; , tooth pair #1. 51 At points A, B, C and D, corner contacts are observed for N pi of the perfect involute gear, corresponding to the time instants when meshing tooth pairs come into or out of contact. These introduce discontinuous points in the slope of the ∑ N pi profile.
Note that N p1 and N p2 between A and B (or C and D) are in phase with each other, which should amplify the peak-to-peak variation of ∑ N pi . For the Ffpi profile of Figure
2.13(c), an abrupt change in the direction is observed at the pitch point P in addition to
the corner contacts. Unlike N pi , the profiles of Ffp1 and Ffp2 of Figure 2.13(c) between A and B (or C and D) are out of phase with each other. This should minimize the peak-to- peak variation of ∑ Ffpi . When tip relief is applied in Figure 2.13(b), corner contacts of
N pi are reduced and smoother transitions are observed at points A, B, C and D. Unlike
the perfect involute gear, N p1 and N p2 are now out of phase with each other between A and B (or C and D), which reduces the peak-to-peak variation of ∑ N pi . However, the
profiles of Ffp1 and Ffp2 of Figure 2.13(d) are in phase with each other in the same region, which amplifies the variation of ∑ Ffpi . The out of phase relationship between N pi and
Ffpi explains why the tip relief (designed to minimize the STE) tends to increase the friction force excitations. This relationship is mathematically embedded in Eq. (2.10) and
graphically illustrated in Figure 2.13, where N p1 and N p2 are in phase while Ffp1 and
Ffp2 are out of phase. Consequently, a compromise would be needed to simultaneously address the dynamic responses in both the LOA and OLOA directions.
52 2.6.4 Prediction of the Dynamic Responses
Dynamic responses including xp ()t , ytp (), Fpbx ()t , Fpby ()t and DTE δ ()t are predicted by numerically integrating the governing equations. Predictions from both perfect and tip relief gears are compared to examine the effect of profile modification in
the presence of sliding friction. Figure 2.14 shows that the normalized xp ()t at 500 lb-in is much smaller (over 90% reduction) when the tip relief is applied. This is because that the STE is the most dominant excitation in the LOA direction and it is minimized at 500 lb-in when the tip relief is applied. An alternate explanation is that the peak-to-peak variation of ∑ N pi is minimized with the tip relief as shown in Figure 2.13.
In the OLOA direction, more significant oscillations are observed for ytp () due to increased ∑ Ffpi excitations with tip relief. Despite that the vibratory components of
∑ N pi are larger than that of ∑ Ffpi , predicted ytp () is actually higher than xp ()t . This shows the necessity of including sliding friction when other excitations such as the STE
are minimized. Note that a phase difference is present in simulated ytp () before and after the tip relief is applied. Predicted pinion bearing forces are not shown here since they depict the same features as the displacement responses of Figure 2.14.
Figure 2.15 shows the DTE predictions, as defined by Eq. (2.17), with and without the tip relief. Similarity between Figure 2.14(a-b) and Figure 2.15(a-b) suggests that the relative LOA displacement plays a dominant role in the DTE responses. However, this conclusion is somewhat case specific as the DTE results depend on the mesh stiffness, bearing stiffness, and gear geometry.
53
Figure 2.14 Dynamic shaft displacements predicted for Example II at 500 lb-in, 4875
RPM and 140 °F with tc = 0.44 ms and fm = 2275 Hz: (a) xp ()t ; (b) Xfp (); (c) ytp () and
(d) Yfp (). Key: , gear pair with perfect involute profile in t domain (II-A); o, with perfect involute profile in f domain (II-A); , with tip relief (II-B).
54
Dynamic transmission error (in)
Figure 2.15 Dynamic transmission error predicted for Example II at 500 lb-in, 4875 RPM and 140 °F with tc = 0.44 ms and fm = 2275 Hz: (a)δ ()t ; (b) ∆()f . Key: , gear pair with perfect involute profile in t domain (II-A); o, with perfect involute profile in f domain (II-A); , with tip relief (II-B). 55 2.7 Experimental Validation of Example II Models
Experiments corresponding to Example II-B were conducted at the NASA Glenn
Research Center (Gear Noise Rig) to validate the MDOF spur gear pair model and to establish the relative influence of friction force excitation on the system. Figure 2.16 shows the inside of the gearbox, where a bracket was built to hold two shaft displacement probes one inch away from the center of the gear in the LOA and OLOA directions [2.13].
The probes face a steel collar that was machined to fit around the output shaft with minimal eccentricity. Accelerometers were mounted on the bracket, so the motion of the displacement probes could be subtracted from the measurements, if necessary. A thermocouple was installed inside the gearbox to measure the temperature of the oil flinging off the gears as they enter into mesh. The thermocouple position was chosen to be consistent with Benedict and Kelley’s [2.14] experiment. A common shaft speed of
4875 rpm is used in all tests so that the first five harmonics of the gear mesh frequency
(2275, 4550, 6825, 9100, and 11375 Hz) do not excite system resonances. Data of shaft displacement in the LOA and OLOA directions are collected from the proximity sensors under oil inlet temperatures over the range of temperatures (104, 122, 140, 158, and 176
°F). At each temperature the torque is varied from 500 to 900 lb-in increments of 100 lb- in.
56 Thermocouple OLOA
LOA
Proximity Probes Bracket Accelerometer
Figure 2.16 Sensors inside the NASA gearbox (for Example II-B).
Parametric studies are conducted to examine the dynamic responses under varying operational conditions of temperature and nominal torque. Benedict and Kelly’s [2.14] friction model is used to calculate the empirical µ as given in Table 2.3 and the realistic kt()calculated using FE/CM under varying torques are incorporated into the dynamic model. Since the precise parameters of the experimental system are not known [2.13], both simulated and measured data are normalized with respect to the amplitude of their first mesh harmonic of the OLOA displacement (which is then designed as 100%). This
57 facilitates the comparison of trends and allows simulations and measurements to be viewed in the same graphs from 0 to 100%.
Figure 2.17 compares the first five mesh harmonics of the LOA displacement as a function of mean torque. (In this and other figures, predictions are shown as continuous lines for the sake of clarity though they are calculated only at discrete points like the measurements.) It was observed that overall simulation trend matches well with the experiment. Magnitudes of the first two mesh harmonics are most dominant and they have minimum values around the optimized load due to the linear tip relief. Figure 2.17 also shows predicted first two harmonics for the prefect involute gear (II-A). Compared with the tip relief gear, they increase monotonically with the mean torque and have much higher values than the tip relief gear around the “optimal” torque.
Figure 2.18 compares the first five mesh harmonics of the OLOA displacement, on a normalized basis, as a function of mean torque. The overall simulation trend again matches well with the experiment. However, unlike the LOA responses, the first harmonic of OLOA displacement grows monotonically with an increase in the mean torque. This is because the friction forces increase almost proportionally with normal loads as predicted by the Coulomb law, but the frictional contribution of each meshing tooth pair tends to be in phase with each other for the tip relief gear (II-B). Thus it should amplify the combined friction force excitation in the OLOA direction. Consequently it is not reducing the OLOA direction responses induced by the sliding friction, even though the profile modification can be efficiently used to minimize gear vibrations in the LOA direction.
58
Figure 2.17 Mesh harmonic amplitudes of X p as a function of the mean torque at 140 °F.
All values are normalized with respect to the amplitude ofYp at the first mesh harmonic. Key: , n = 1 (prediction of II-B); , n = 2 (prediction of II-B); , n = 3 (prediction of II-B); , n = 1 (prediction of II-A); , n = 2 (prediction of II-A); , n = 1 (measurement of II-B); , n = 2 (measurement of II-B); O, n = 3 (measurement of II-B).
59
Figure 2.18 Mesh harmonic amplitudes of Yp as a function of the mean torque at 140 °F for Example II-B. All values are normalized with respect to the amplitude ofYp at the first mesh harmonic. Key: , n = 1 (prediction of II-B); , n = 2 (prediction of II-B); , n = 3 (prediction of II-B); , n = 1 (measurement of II-B); , n = 2 (measurement of II-B); O, n = 3 (measurement of II-B).
60 Figure 2.19 compares first five mesh harmonics of the normalized DTE for II-A and II-B cases over a range mean torques. Observe that the DTE spectral trends are very similar to the STE spectral trends of Figure 2.10. For example, the harmonic amplitudes of the perfect involute gear grow monotonically with mean torque while the harmonic amplitudes of the tip relief gear have minimum values around the “optimal” torque. Also, the DTE spectra show a dominant second harmonic, whose magnitude is comparable to that at the first harmonic. In some cases for the tip relief gear the second harmonic becomes the most dominant component especially when the mean torque is lower than
350 lb-in.
Finally, Figure 2.20 compares the first five mesh harmonics of the normalized
LOA displacement as a function of operational temperature. The changes in temperature are converted into variation in µ of Table 2.3. Compared with the OLOA motions, both predictions and measurements in the LOA direction give almost identical results at all temperatures. Figure 2.21 shows the first five mesh harmonics of the OLOA displacement with a change in temperature. The first harmonic varies quite significantly even though the changes in µ are relatively small. Consequently, the OLOA dynamics tends to be much more sensitive to a variation in µ as compared with the LOA motions.
Measured data of Figure 2.21 show some variations due to the experimental errors [2.13].
61
Figure 2.19 Predicted dynamic transmission errors (DTE) for Example II over a range of torque at 140 °F: (a) gear pair with perfect involute profile (II-A); (b) with tip relief (II-B). All values are normalized with respect to the amplitude of δ (II-A) at the first mesh harmonic with 100 lb-in. Key: , n = 1 (II-B); , n = 2 (II-B); , n = 3 (II-B). 62
Figure 2.20 Mesh harmonic amplitudes of X p as a function of temperature at 500 lb-in for Example II-B. All values are normalized with respect to the amplitude ofYp at the first mesh harmonic. Key: , n = 1 (prediction of II-B); , n = 2 (prediction of II-B); , n = 3 (prediction of II-B); , n = 1 (measurement of II-B); , n = 2 (measurement of II-B); O, n = 3 (measurement of II-B).
63
100
80
60
40
20
0 95 115 135 155 175 Temperature (deg F)
Figure 2.21 Mesh harmonic amplitudes of Yp as a function of temperature at 500 lb-in for
Example II-B. All values are normalized with respect to the amplitude ofYp at the first mesh harmonic. Key: , n = 1 (prediction of II-B); , n = 2 (prediction of II-B); , n = 3 (prediction of II-B); , n = 1 (measurement of II-B); , n = 2 (measurement of II-B); O, n = 3 (measurement of II-B).
64 2.8 Conclusion
Chief contribution of this study is the development of a new multi-degree of freedom, linear time-varying model. This formulation overcomes the deficiency of
Vaishya and Singh’s work [2.1-2.3] by employing realistic tooth stiffness functions and the sliding friction over a range of operational conditions. Refinements include: (1) an accurate representation of tooth contact and spatial variation in tooth mesh stiffness based on a FE/CM code in the “static” mode; (2) Coulomb friction model for sliding resistance with empirical coefficient of friction as a function of operation conditions; (3) a better representation of the coupling between the LOA and OLOA directions including torsional and translational degrees of freedom. Numerical solutions of the MDOF model yield the dynamic transmission error and vibratory motions in the LOA and OLOA directions. The new model has been successfully validated first by using the FE/CM code while running in the “dynamic” mode and then by analogous experiments. Since the lumped model is more computationally efficient when compared with the FE/CM analysis, it could be quickly used to study the effect of a large number of parameters.
One of the main effects of sliding friction is the enhancement of the DTE magnitude at the second gear mesh harmonic. A key question whether the sliding friction is indeed the source of the OLOA motions and forces is then answered by our model. The bearing forces in the LOA direction are influenced by the normal tooth loads, but the sliding frictional forces primarily excite the OLOA motions. Finally, effect of the profile modification on the dynamic transmission error has been analytically examined under the influence of frictional effects. For instance, the tip relief introduces an amplification in the OLOA motions and forces due to an out of phase relationship between the normal 65 load and friction forces. This knowledge should be of significant utility to the designers.
Future modeling work should examine the effects of other profile modifications and find the conditions for minimal dynamic responses when both STE and friction excitations are simultaneously present. Also, the model could be further refined by incorporating alternate friction formulations.
References for Chapter 2
[2.1] Vaishya, M., and Singh, R., 2001, “Analysis of Periodically Varying Gear Mesh Systems with Coulomb Friction Using Floquet Theory,” Journal of Sound and Vibration, 243(3), pp. 525-545.
[2.2] Vaishya, M., and Singh, R., 2001, “Sliding Friction-Induced Non-Linearity and Parametric Effects in Gear Dynamics,” Journal of Sound and Vibration, 248(4), pp. 671- 694.
[2.3] Vaishya, M., and Singh, R., 2003, “Strategies for Modeling Friction in Gear Dynamics,” ASME Journal of Mechanical Design, 125, pp. 383-393.
[2.4] Houser, D. R., Vaishya M., and Sorenson J. D., 2001, “Vibro-Acoustic Effects of Friction in Gears: An Experimental Investigation,” SAE Paper # 2001- 01-1516.
[2.5] Velex, P., and Cahouet, V., 2000, “Experimental and Numerical Investigations on the Influence of Tooth Friction in Spur and Helical Gear Dynamics,” ASME Journal of Mechanical Design, 122(4), pp. 515-522.
[2.6] Velex, P., and Sainsot. P, 2002, “An Analytical Study of Tooth Friction Excitations in Spur and Helical Gears,” Mechanism and Machine Theory, 37, pp. 641-658.
[2.7] Lundvall, O., Strömberg, N., and Klarbring, A., 2004, “A Flexible Multi-body Approach for Frictional Contact in Spur Gears,” Journal of Sound and Vibration, 278(3), pp. 479-499.
[2.8] External2D (CALYX software), 2003, “Helical3D User’s Manual,” www.ansol.www, ANSOL Inc., Hilliard, OH.
66
[2.9] Vinayak, H., Singh, R., and Padmanabhan, C., 1995, “Linear Dynamic Analysis of Multi-Mesh Transmissions Containing External, Rigid Gears,” Journal of Sound and Vibration, 185(1), pp. 1-32.
[2.10] Lim, T. C., and Singh, R., 1990, “Vibration Transmission Through Rolling Element Bearings. Part I: Bearing Stiffness Formulation,” Journal of Sound and Vibration, 139(2), pp. 179-199.
[2.11] Lim, T. C., and Singh, R., 1990, “Vibration Transmission Through Rolling Element Bearings. Part II: System Studies,” Journal of Sound and Vibration, 139(2), pp. 201-225.
[2.12] Kahraman, A., and Singh, R., 1991, “Error Associated with a Reduced Order Linear Model of Spur Gear Pair,” Journal of Sound and Vibration, 149(3), pp. 495-498.
[2.13] Singh, R., 2005, “Dynamic Analysis of Sliding Friction in Rotorcraft Geared Systems,” Technical Report submitted to the Army Research Office, grant number DAAD19-02-1-0334.
[2.14] Benedict, G. H., and Kelley B. W., 1961, “Instantaneous Coefficients of Gear Tooth Friction,” Transactions of the American Society of Lubrication Engineers, 4, pp. 59-70.
[2.15] Rebbechi, B., Oswald, F. B., and Townsend, D. P., 1996, “Measurement of Gear Tooth Dynamic Friction,” ASME Power Transmission and Gearing Conference proceedings, DE-Vol. 88, pp. 355-363.
67 CHAPTER 3
PREDICTION OF DYNAMIC FRICTION FORCES USING ALTERNATE
FORMULATIONS
3.1 Introduction
Gear dynamic researchers [3.1-3.6] have typically modeled sliding friction phenomenon by assuming Coulomb formulation with a constant coefficient (µ) of friction
(it is designated as Model I in this chapter). In reality, tribological conditions change
continuously due to varying mesh properties and lubricant film thickness as the gears roll
through a full cycle [3.7-3.10]. Thus, µ varies instantaneously with the spatial position of
each tooth and the direction of friction force changes at the pitch point. Alternate
tribological theories, such as elasto-hydrodynamic lubrication (EHL), boundary
lubrication or mixed regime, have been employed to explain the interfacial friction in
gears [3.7-3.10]. For instance, Benedict and Kelley [3.7] proposed an empirical dynamic
friction coefficient (designated as Model II) under mixed lubrication regime based on
measurements on a roller test machine. Xu et al. [3.8, 3.9] recently proposed yet another
friction formula (designated as Model III) that is obtained by using a non-Newtonian,
thermal EHL formulation. Duan and Singh [3.11] developed a smoothened Coulomb 68 model for dry friction in torsional dampers; it could be applied to gears to obtain a smooth transition at the pitch point and we designate this as Model IV. Hamrock and
Dawson [3.10] suggested an empirical equation to predict the minimum film thickness for two disks in line contact. They calculated the film parameter Λ, which could lead to a composite, mixed lubrication model for gears (designated as Model V). Overall, no prior work has incorporated either the time-varying µ()t or Models II to V, into multi-degree- of-freedom (MDOF) gear dynamics. To overcome this void in the literature, specific objectives of this chapter are established as follows: 1. Propose an improved MDOF spur gear pair model with time-varying coefficient of friction, µ()t , given realistic mesh stiffness profiles of Chapter 2; 2. Comparatively evaluate alternate sliding friction models and predict the interfacial friction forces and motions in the off-line-of-action (OLOA) direction; and 3. Validate one particular model (III) by comparing predictions to the benchmark gear friction force measurements made by Rebbechi et al. [3.12].
3.2 MDOF Spur Gear Model
Transitions in key meshing events within a mesh cycle are determined from the undeformed gear geometry. Figure 3.1(a) is a snapshot for the example gear set (with a contact ratio σ of about 1.6) at the beginning (t = 0) of the mesh cycle (tc). At that time, pair # 1 (defined as the tooth pair rolling along line AC) just comes into mesh at point A and pair # 0 (defined as the tooth pair rolling along line CD) is in contact at point C, which is the highest point of single tooth contact (HPSTC). When pair # 1 approaches the
69 lowest point of single tooth contact (LPSTC) at point B, pair # 0 leaves contact. Further, when pair #1 passes through the pitch point P, the relative sliding velocity of the pinion with respect to the gear is reversed, resulting in a reversal of the friction force. Beyond point C, pair # 1 will be re-defined as pair # 0 and the incoming meshing tooth pair at point A will be re-defined as pair # 1, resulting in a linear-time-varying (LTV) formulation. The spur gear system model is shown in Figure 3.1(b) and key assumptions for the dynamic analysis include the following: (i) pinion and gear are rigid disks; (ii) shaft-bearings stiffness elements in the line-of-action (LOA) and OLOA directions are modeled as lumped springs which are connected to a rigid casing; (iii) vibratory angular motions are small in comparison to the kinematic motion. Overall, we obtain a LTV system formulation, as explained in Chapter 2 with a constant µ . Refinements to the
MDOF model of Figure 3.1(b) with time-varying sliding friction µ()t are proposed as
follows. The governing equations for the torsional motions θ p ()t and θg ()t are as follows:
nn==floor(σσ ) floor( ) J pθ p()tT=+ p∑∑ XtFt pi () ⋅ pfi () − rNt bp ⋅ pi (), (3.1) ii==00
nn==floor(σσ ) floor( ) J gθ g()tT=− g +∑∑ XtFt gi () ⋅ gfi () + rNt bg ⋅ gi (). (3.2) ii==00
70
(a)
−Tgg,−Ω
J g ,mg K gBy
K kt() gBx
K pBy J pp,m
x K pBx Ω pp,T θ y
(b)
Figure 3.1 (a) Snap shot of contact pattern (at t = 0) in the spur gear pair; (b) MDOF spur gear pair system; here kt() is in the LOA direction.
71 Here, the “floor” function rounds off the contact ratio σ to the nearest integer
(towards a lower value); J p and J g are the polar moments of inertia for the pinion and
gear; Tp and Tg are the external and braking torques; N pi ()t and Ngi ()t are the normal
loads defined as follows:
⎡⎤ Ntpi()== Nt gi () ktr i ()⎣⎦ bp ⋅−⋅+−+θθ p () t r bg g () t xt p () xt g () (3.3) ⎡⎤ ctibppbggp()⎣⎦ r⋅−⋅+−θ () t rθσ () t x () t x g () t , i = 0, 1, ... , n = floor( ).
where kti ( ) and cti ( ) are the time-varying realistic mesh stiffness and viscous damping
profiles; rbp and rbg are the base radii of the pinion and gear; xp ()t and xg ()t denote the
translational displacements (in the LOA direction) at the bearings. The sliding (interfacial)
th friction forces Fpfi ()t and Fgfi ()t of the i meshing pair are derived as follows; note that
five alternate µ()t models will be described later.
Fpfi()ttNt= µ () pi (), Fgfi()ttNti==µ () gi (), 0, ... , n . (3.4a,b)
th The frictional moment arms Xtpi () and Xtgi () acting on the i tooth pair are:
Xtpi( )=+−+ L XA ( ni )λλ mod( Ω p rt bp , ), i = 0, ... , n, (3.5a)
Xtgi( )=+−Ω L YC iλλ mod( g rt bg , ), i = 0, ... , n. (3.5b)
72 where “mod” is the modulus function defined as: mod(xy , )= x−⋅ y floor( x / y ), if y ≠ 0 ;
“sgn” is the sign function; Ω p and Ωg are the nominal operational speeds (in rad/s); and λ
is the base pitch. Refer to Figure 3.1(a) for length LAP . The governing equations for the
translational motions xp ()t and xg ()t in the LOA direction are:
n=floor(σ ) mxp p() t+⋅⋅++= 2ζ pBx K pBx m p x p () t K pBx x p () t∑ N pi () t 0 , (3.6) i=0
n=floor(σ ) mxg g() t+⋅⋅++= 2ζ gBx K gBx m g x g () t K gBx x g () t∑ N gi () t 0. (3.7) i=0
Here, mp and mg are the masses of the pinion and gear; K pBx and K gBx are the
effective shaft-bearing stiffness values in the LOA direction, and ζ pBx and ζ gBx are their
damping ratios. Likewise, the governing equations for the translational motions ytp ()
and ytg () in the OLOA direction are written as:
n=floor(σ ) myp p() t+⋅⋅+−= 2ζ pB y K pBy m p y p () t K pBy y p () t∑ F pfi () t 0, (3.8) i=0
n=floor(σ ) myg g() t+⋅⋅+−= 2ζ gB y K gBy m g y g () t K gBy y g () t∑ F gfi () t 0. (3.9) i=0
73 3.3 Spur Gear Model with Alternate Sliding Friction Models
Following a similar modeling strategy of Chapter 2, we obtain a LTV system formulation. Refinements to the MDOF model with time-varying sliding friction µ()t are
th proposed as follows. The sliding (interfacial) friction forces Fpfi ()t and Fgfi ()t of the i
meshing pair are
Fpfi()ttNt= µ () pi (), Fgfi()ttNti==µ () gi (), 0, ... , n . (3.10a,b)
Five alternate µ()t models are described as follows:
3.3.1 Model I: Coulomb Model
The Coulomb friction model with time-varying (periodic) coefficient of
th frictionµCi (t ) for the i meshing tooth pair is derived as follows, where µavg is the
magnitude of the time-average.
⎡⎤ µµCi(trtniLin )=⋅ avg sgn⎣⎦ mod( Ω p bp , λ ) +−− ( ) λ AP , = 0, ... , . (3.11)
74 3.3.2 Model II: Benedict and Kelley Model
The instantaneous profile radii of curvature (mm) ρ()t of ith meshing tooth are:
ρ pi(tL )=+−+ XA ( ni )λλ mod( Ω p rt bp , ) , in= 0, ... , . (3.12a)
ρρgi(tL )=− XY pi ( ti ), = 0, ... , n. (3.12b)
th The rolling (tangential) velocities vtr ( )(m/s) of i meshing tooth pair are:
Ω ρ ()t Ω ρ ()t vt()= ppi , vt( )==ggi , i 0, ... , n. (3.13a,b) rpi 1000 rgi 1000
th The sliding velocity vts ( ) and the entraining velocity vte ( ) (m/s) of i meshing tooth pair are:
vtsi()=− vtvt rpi () rgi (), vtei( )=+ vt rpi ( ) vt rgi ( ) , i = 0, ... , n. (3.14a,b)
The unit normal load (N/mm) is wTZrnp=⋅/cos( wpα ) , where α is the pressure
angle, Z is the face width (mm), Tp is torque (N-mm) and rwp is the operating pitch radius of pinion (mm). Our µ()t prediction for the ith meshing tooth pair is based on the
Benedict and Kelley model [3.7], though it is modified to incorporate a reversal in the
direction of friction force at the pitch point. Here, SSSavg=+0.5( ap ag ) is the averaged
75 surface roughness ( µm ), and ηM is the dynamic viscosity of the oil entering the gear contact.
0.0127× 1.13 ⎡⎤29700w µλλ(trtniL )=⋅ logn ⋅Ω+−− sgn⎡ mod( , ) ( ) ⎤ , Bi 10 ⎢⎥2 ⎣ pbp AP⎦ 1.13− Svtvtavg⎣⎦η M si ( ) ei ( )
in= 0, ... , . (3.15)
3.3.3 Model III: Formulation Suggested by Xu et al.
th The composite relative radius of curvature ρr (t ) (mm) of i meshing tooth pair is:
ρ pi()ttρ gi () ρri ()t = , in= 0, ... , (3.16) ρρpi()tt+ gi ()
The effective modulus of elasticity (GPa) of mating surfaces is
22 ⎡⎤11−−ν p ν g E′ =+2/⎢⎥, where E and ν are the Young’s modulus and Poisson’s ratio, ⎣⎦⎢⎥EEpg respectively. The maximum Hertzian pressure (GPa) for the ith meshing tooth pair is:
wEn ′ Pthi ()= , in= 0, ... , . (3.17) 2000πρri (t )
76 Define the dimensionless slide-to-roll ratio SR() t and oil entraining velocity
th Vte () (m/s) of i meshing tooth pair as:
2()vtsi vtei () SRi () t = , Vtei ()= , in= 0, ... , . (3.18a,b) vtei () 2
The empirical sliding friction expression (for the ith meshing tooth pair), as proposed by Xu et al. based on non-Newtonian, thermal EHL theory [3.8-3.9], is modified in our work to incorporate a reversal in the direction of the friction force at the pitch point:
fSRtP()ihiMavg(), (), tη , S b b3 bbb 2 678 ⎡ ⎤ µηλλXi()t=⋅Ω+−− e PSRtV hi i () ei () t M Rt i ()sgnmod(⎣ p rt bp , ) ( ni ) L AP ⎦ ,
− SRihi() t P ()log( t 10 η M ) Savg f( SRtPi(), hi (), tηη M , S avg) =+ b14 bSRtP i () hi ()log( t 10 M ) + be 5 + be 9 ,
in= 0, ... , . (3.19a,b)
Xu [3.9] suggested the following empirical coefficients (in consistent units) for
the above formula: b1 =−8.916465 , b2 =1.03303 , b3 = 1.036077 , b4 =−0.354068 ,
b5 = 2.812084 , b6 =−0.100601, b7 = 0.752755, b8 =−0.390958 , and b9 = 0.620305 .
77 3.3.4 Model IV: Smoothened Coulomb Model
Xu [3.9] conducted a series of friction measurements on a ball-on-disk test machine and measured the µ()t values as a function of SR; these results resemble the smoothening function reported by Duan and Singh [3.11] near the pitch point (SR = 0) especially at very low speeds (boundary lubrication conditions). By denoting the
th periodic displacement of i meshing tooth pair as xipbpAP(trtniL )=Ω mod( ,λλ ) +−− ( ) , a smoothening function could be used in place of the discontinuous Coulomb friction of
Chapter 2. The arc-tangent type function is proposed as follows though one could also use other functions [3.11]:
22µavgµσ avg µSi(txtxt )=Φ⋅+⋅ arctan[] i ( ) i ( ) , in= 0, ... , (3.20) π ⎡ 22 ⎤ π ⎣1()+Φ xi t ⎦
In the above, the regularizing factor Φ is adjusted to suit the need of smoothening requirement. A higher value of Φ corresponds to a steeper slope at the pitch point. In our work, Φ = 50 is used for a comparative study.
3.3.5 Model V: Composite Friction Model
Alternate theories (Models I to IV) seem to be applicable over specific operational conditions. This necessitates a judicious selection of an appropriate lubrication regime as indicated by the film parameter, Λ, that is defined as the ratio of minimum lubrication
78 22 film thickness and composite surface roughness RRRcomp=+ rms,, g rms p measured with a
filter cutoff wave length Lx , where Rrms is the rms gear-tooth surface roughness [3.13].
The film parameter for rotorcraft gears usually lies between 1 and 10. In the mixed lubrication regime the films are sufficiently thin to yield partial asperity contact, while in the EHL regime the lubrication film completely separates the gear surfaces. Accordingly, a composite friction model is proposed as follows:
⎧µC (t ) simplified Coulomb model, computationally efficient (Model I) ⎪ ⎪ µB (t ) 1<Λ <4, mixed lubrication, (Model II) µ()t = ⎨ (3.21) µ (t ) 4≤Λ <10, EHL lubrication, (Model III) ⎪ X ⎪ ⎩µSpp(tT ) low ΩΛ< , high , 1, boundary lubrication (Model IV)
Application of Model II, III or IV would, of course, depend on the operational and tribological conditions though Model I could be easily utilized for computationally
efficient dynamic simulations. Note that the magnitude µavg of Model I or IV should be determined separately. For instance, the averaged coefficient based on Model II was used in Chapter 2. Also, the critical Λ value between different lubrication regimes must be carefully chosen. The film thickness calculation employs the following equation developed by Hamrock and Dowson [3.10, 3.13], based on a large number of numerical solutions that predict the minimum film thickness for two disks in line contact. Here, G is the dimensionless material parameter, W is the load parameter, U is the speed parameter,
H is the dimensionless central film thickness, and bH is the semi-width of Hertzian contact band:
79
3 Htci()ρ r1 () t× 10 LX Λ=i ()t , in= 0, ... , (3.22a) Rcomp2()bt Hi
8()wtnrρ 1 s btH1()= , Gk= ηMr E, (3.22b-c) π Er
0.56 0.69 GUi () t Htci ( )= 3.06 0.10 , (3.22d) Wti ()
ηMeivt() −6 wn Uti ()=× 10 , Wti ()= . (3.22e.f) 2()Etrriρ Erriρ ()t
3.4 Comparison of Sliding Friction Models
Figure 3.2(a) compares the magnitudes of µ()t as predicted by Model II and III
for the spur gear set of Chapter 2 given Tp = 22.6 N-m (200 lb-in) and Ω p = 1000 RPM.
The LTV formulations for meshing tooth pairs # 0 and 1 result in periodic profiles for both models. Two major differences between these two models are: (1) The averaged magnitude from Model II is much higher compared with that of Model III since friction under mixed lubrication is generally higher than under EHL lubrication; and (2) while
Model III predicts nearly zero friction near the pitch point, Model II predicts the largest µ value due to the entraining velocity term in the denominator.
80 0.1
0.09
0.08
0.07
0.06 µ 0.05
0.04
0.03
0.02
0.01
0 0 0.5 1 1.5 2 Normalized time t/t c (a)
0.11
0.1
0.09
0.08
0.07 (t)
µ 0.06
0.05
0.04
0.03
0.02 0 1000 2000 3000 4000 5000 Ω (RPM) p (b)
Figure 3.2 (a) Comparison of Model II [3.7] and Model III [3.8] given Tp = 22.6 N-m
(200 lb-in) and Ω p = 1000 RPM. Key: , pair # 1 with Model II; , pair # 0 with Model II; , pair # 1 with Model III; , pair # 0 with Model III; (b) Averaged magnitude of the coefficient of friction predicted as a function of speed using the composite Model V with Tp = 22.6 N-m (200 lb-in). Here, tc is one mesh cycle.
81 As explained by Xu [3.9], three different regions could be roughly defined on a µ versus SR curve. When the sliding velocity is zero, there is no sliding friction, and only rolling friction (though very small) exists. Thus, the µ value should be almost zero at the pitch point. When the SR is increased from zero, µ first increases linearly with small values of SR. This region is defined as the linear or isothermal region. When the SR is increased slightly further, µ reaches a maximum value and then decreases as the SR value is increased beyond that point. This region is referred to as non-linear or non-Newtonian region. As the SR is increased further, the friction decreases in an almost linear fashion; this is called as the thermal region. Model II seems to be valid only in the thermal region
[3.8, 3.9]. Figure 3.2(b) shows the averaged magnitude of µavg predicted as a function of
Ω p using the composite formulation (Model V) with Tp = 22.6 N-m (200 lb-in). An abrupt change in magnitude is found around 2500 RPM corresponding to a transition from the EHL to a mixed lubrication regime. Similar results could be obtained by plotting
the composite µ()t as a function of Tp . Though our composite model could be used to predict µ()t over a large range of lubrication conditions, care must be exercised since the calculation of Λ itself is based on an empirical equation [3.10].
Figure 3.3 compares four friction models on a normalized basis. The curves
between 0≤
82 1.5
1
(t) 0.5 µ
0
-0.5 Normalized Normalized
-1
-1.5
0 0.5 1 1.5 Normalized time t/t c
Figure 3.3 Comparison of normalized friction models. Key: , Model I (Coulomb friction with discontinuity); , Model II [3.7]; , Model III [3.8]; , Model
IV (smoothened Coulomb friction). Note that curve between 0≤ tt /c < 1 is for pair # 1; and the curve between 1≤
Figure 3.4 compares the combined normal loads and friction force time histories
as predicted by four friction models given Tp = 56.5 N-m (500 lb-in) and Ω p = 4875
RPM. Note that while Figure 3.3 illustrates µ()t for each meshing tooth pair the friction forces of Figure 3.4 include the contributions from both (all) meshing tooth pairs. Though alternate friction formulations dictate the dynamic friction force profiles, they have negligible effect on the normal loads.
83
1400
1350 (N) p N
1300
0 0.5 1 1.5 2
100
50
(N) 0 fp F -50
-100
0 0.5 1 1.5 2 Normalized time t/t c
Figure 3.4 Combined normal load and friction force time histories as predicted using alternate friction models given Tp = 56.5 N-m (500 lb-in) and Ω p = 4875 RPM. Key: , Model I; , Model II; , Model III; , Model IV.
84 3.5 Validation and Conclusion
Figure 3.5 compares the predicted LOA and OLOA displacements with alternate
friction models given Tp = 56.5 N-m (500 lb-in) and Ω p = 4875 RPM. Note that the differences between predicted motions are not significant though friction formulations and friction force excitations differ. This implies that one could still employ the simplified Coulomb formulation (Model I) in place of more realistic time-varying friction models (Models II to IV). Similar trend is observed in Figure 3.6 for the dynamic
transmission errors (DTE), defined as δθ()tr=−bp p () tr bg θ g () t +−xpg()txt (). The most significant variation induced by friction formulation is at the second harmonic, which matches the results reported by Lundvall et al. [3.6].
Finally, predicted normal load and friction force time histories (with Model III) are validated using the benchmark friction measurements made by Rebbechi et al. [3.12].
Results are shown in Figure 3.7. Based on the comparison, µ is found to be about 0.004 since it was not given in the experimental study. Here, we have made the periodic LTV definitions of meshing tooth pairs # 0 and 1 to be consistent with those of measurements, where meshing tooth pairs A and B are labeled in a continuous manner. Predictions
match well with measurements at both low ( Ω p = 800 RPM) and high ( Ω p = 4000 RPM) speeds. Ongoing research focuses on the development of semi-analytical solutions given a specific µ()t model and an examination of the interactions between tooth modifications and sliding friction.
85
0.15 -59.6
0.1 m) -59.8 m) µ µ ( ( p p x x -60 0.05
-60.2 0 0 1 2 1 2 3 4 5
2 1 1 0.8 m) m)
µ 0 µ 0.6 ( ( p p y y 0.4 -1 0.2 -2 0 0 1 2 1 2 3 4 5 Normalized time t/t Mesh order n c
Figure 3.5 Predicted LOA and OLOA displacements using alternate friction models given
Tp = 56.5 N-m (500 lb-in) and Ω p = 4875 RPM. Key: in time domain: , Model I; , Model II; , Model III; , Model IV; in frequency (mesh order n) domain: + , Model I; , Model II; , , Model III; +, Model IV.
86
24
23.5 m)
µ 23
22.5 DTE ( DTE 22
21.5 0 0.5 1 1.5 2 Normalized time t/t c
0.4
0.3 m) µ 0.2
DTE ( DTE 0.1
0 1 2 3 4 5 Mesh order n
Figure 3.6 Predicted dynamic transmission error (DTE) using alternate friction models given Tp = 56.5 N-m (500 lb-in) and Ω p = 4875 RPM. Key: in time domain: , Model I; , Model II; , Model III; , Model IV; in frequency (mesh order n) domain: + , Model I; , Model II; , , Model III; +, Model IV.
87 2000
1500
(N) 1000 p N 500
0 0 0.5 1 1.5 2
5
0 (N) fp F -5
0 0.5 1 1.5 2 Normalized time t/t c (a) 2500
2000
1500 (N) p
N 1000
500
0 0 0.5 1 1.5 2
4
2
0 (N) fp
F -2
-4
-6 0 0.5 1 1.5 2 Normalized time t/t c (b)
Figure 3.7 Validation of the normal load and sliding friction force predictions: (a) at Tp =
79.1 N-m (700 lb-in) and Ω p = 800 RPM; (b) at Tp = 79.1 N-m (700 lb-in) and Ω p = 4000 RPM. Key: , prediction of tooth pair A with Model III; , prediction of tooth pair B with Model III; X, measurement of tooth pair A [3.12]; , , measurement of tooth pair B [3.12]. 88 References for Chapter 3
[3.1] Vaishya, M., and Singh, R., 2001, “Analysis of Periodically Varying Gear Mesh Systems with Coulomb Friction Using Floquet Theory,” Journal of Sound and Vibration, 243(3), pp. 525-545.
[3.2] Vaishya, M., and Singh, R., 2001, “Sliding Friction-Induced Non-Linearity and Parametric Effects in Gear Dynamics,” Journal of Sound and Vibration, 248(4), pp. 671- 694.
[3.3] Vaishya, M., and Singh, R., 2003, “Strategies for Modeling Friction in Gear Dynamics,” ASME Journal of Mechanical Design, 125, pp. 383-393.
[3.4] Velex, P., and Cahouet, V., 2000, “Experimental and Numerical Investigations on the Influence of Tooth Friction in Spur and Helical Gear Dynamics,” ASME Journal of Mechanical Design, 122(4), pp. 515-522.
[3.5] Velex, P., and Sainsot. P, 2002, “An Analytical Study of Tooth Friction Excitations in Spur and Helical Gears,” Mechanism and Machine Theory, 37, pp 641-658.
[3.6] Lundvall, O., Strömberg, N., and Klarbring, A., 2004, “A Flexible Multi-body Approach for Frictional Contact in Spur Gears,” Journal of Sound and Vibration, 278(3), pp. 479-499.
[3.7] Benedict, G. H., and Kelley B. W., 1961, “Instantaneous Coefficients of Gear Tooth Friction,” Transactions of the American Society of Lubrication Engineers, 4, pp. 59-70.
[3.8] Xu, H., Kahraman, A., Anderson, N. E., and Maddock, D. G., 2007, “Prediction of Mechanical Efficiency of Parallel-Axis Gear Pairs,” ASME Journal of Mechanical Design, 129 (1), pp. 58-68.
[3.9] Xu, H., 2005, “Development of a Generalized Mechanical Efficiency Prediction Methodology,” PhD dissertation, The Ohio State University.
[3.10] Hamrock, B. J., and Dowson, D., 1977, “Isothermal Elastohydrodynamic Lubrication of Point Contacts, Part III – Fully Flooded Results,” Journal of Lubrication Technology, 99(2), pp. 264-276.
[3.11] Duan, C., and Singh, R., 2005, “Super-Harmonics in a Torsional System with Dry Friction Path Subject to Harmonic Excitation under a Mean Torque”, Journal of Sound and Vibration, 285(2005), pp. 803-834.
[3.12] Rebbechi, B., Oswald, F. B., and Townsend, D. P., 1996, “Measurement of Gear Tooth Dynamic Friction,” ASME Power Transmission and Gearing Conference proceedings, DE-Vol. 88, pp. 355-363. 89
[3.13] AGMA Information Sheet 925-A03, 2003, “Effect of Lubrication on Gear Surface Distress.”
90 CHAPTER 4
CONSTRUCTION OF SEMI-ANALYTICAL SOLUTIONS TO SPUR GEAR
DYNAMICS
4.1 Introduction
Periodic differential equations [4.1-4.3] are usually needed to describe the gear
dynamics [4.4-4.12] since significant variations in mesh stiffness k(t) and damping c(t)
are observed, within the fundamental period tc (one mesh cycle). Additionally, dynamic
friction force Ff(t) and torque Mf(t) also undergo periodic variations, with the same period
tc, due to changes in normal mesh loads and coefficient of friction µ, as well as a reversal
in the direction of Ff(t) at the pitch point [4.4-4.6], especially in spur and helical gears.
th For the sake of illustration, typical ki(t) profiles and frictional functions fi(t) for the i meshing pair in spur gears are shown in Figure 4.1; derivations of fi(t) will be explained
later along with particulars of the example case. The fundamental nature of the linear
time-varying (LTV) system is illustrated in Figure 4.2(a); the system model is described in Chapter 2.
91 ) (lb/in) t ( i k
(a) f(t)
(b)
Figure 4.1(a) Realistic mesh stiffness functions of the spur gear pair example (with tip relief) given Tp = 550 lb-in. Key: , kt0 (); , kt1(). (b) Periodic frictional functions. Key: , f0 ()t ; , f1()t ; , f2 ()t . 92
(a)
Ωg
Ω p
(b)
Figure 4.2(a) Normal (mesh) and friction forces of 6DOF analytical spur gear system model. (b) Snap shot of contact pattern (at t = 0) for the example spur gear pair. 93 The governing single degree-of-freedom (SDOF) equation in terms of dynamic
transmission error (DTE) δθ()tr=−bp p () tr bg θ g () t is given below, where subscripts p and
g correspond to the pinion and gear, respectively; θ is the vibratory component of the
rotation; and rb is the base radius.
SS JtJ ()⎡⎤ cttktt () () () () sgnmod(⎡ rt , ) ( Sj ) L⎤ ebiδδδµλλ+++Ω+−−∑∑⎣⎦ i⎣ pbp AP⎦ ii==00 (4.1a) ⋅+⎡⎤ct ()δδ () t kt () () t⎡⎤ X () tJr + X () tJr =+ T Tt () ⎣⎦ii⎣⎦ pigbpgipbge ,
22 Jepg= JJ, Jbgbppbg=+Jr Jr , TTJrTJrepgbpgpbg=+, (4.1b-d)
Xtpi( )=+−+ L XA ( Si )λ mod( Ω p rt bp ,λσ ), i = 0, ... , S = floor( ) , (4.1e)
Xtgi( )=+−Ω L YC iλ mod( g rt bg ,λσ ), i = 0, ... , S = floor( ) . (4.1f)
Further, J is the moment of inertia; T and Ω are the nominal torque and rotation speed;
and λ is the base pitch. Tooth pairs #1 and #0 are defined as the pairs rolling along line
AC and CD in Figure 4.2(b), respectively. The jth tooth pair passes though the pitch point
P during the meshing event, and the reversal at P is characterized by the sign function
“sgn” with a constant coefficient of Coulomb friction µ [4.4]. The modulus function
(mod(x, y) = x − y⋅floor(x/y), if y ≠ 0) is used to describe the periodic friction force Ff(t) and the moment arm X(t). The “floor” function rounds off the contact ratio σ to the nearest integers towards a lower value, i.e. S = 1 for the example case. Finally, L corresponds to the geometric length in Figure 4.2(b).
94 The chief goal of this chapter is to find semi-analytical solutions to Eq. (4.1) type
periodic systems which significantly differ from the classical Hill’s equation [4.1] in
several ways. First, the periodickti ( ) is not confined to a rectangular wave assumed by
Manish and Singh [4.4-4.6], or a simple sinusoid as in the Mathieu’s equation [4.1].
Instead, Eq. (4.1) should describe realistic, yet continuous, profiles of Figure 4.1(a) resulting from a detailed finite element/contact mechanics analysis [4.7]. Hence, multiple
harmonics of kti ( ) should be considered. Second, the periodic viscous cti ( ) term should
dissipate vibratory energy due to the sliding friction besides its kinematic effect. Third,
S S the ∑δii()tc () t and ∑δii()tk () t terms of Eq. (4.1) incorporate combined (but phase i=1 i=1
correlated) contributions from all (yet changing) tooth pairs in contact. Consequently, the
relative phase between neighboring tooth pairs should play an important role in the
resulting response δ ()t . Fourth, multiplicative effects betweenkti ( ) , cti ( ) ,()Xti and
δi ()t should result in higher mesh harmonics, which poses difficulty in constructing
closed-form solutions. Lastly, Tt() of Eq. (4.1) represents the time-varying component of
the forcing function due to unloaded (manufacturing) static transmission errorε ()t . This indicates that the frictional forces and moments reside on both sides of Eq. (4.1) as either periodically-varying parameters or external excitations, thus posing further mathematical complications.
95 4.2 Problem Formulation
Sliding friction has been found as a non-negligible excitation source in spur and
helical gear dynamics by Houser et al. [4.8], Velex and Cahouet [4.9], Velex and Sainsot
[4.10], and Lundvall et al. [4.11]. Earlier, Vaishya and Singh [4.4-4.6] developed a SDOF spur gear model with rectangular k(t) and sliding friction profiles; they solved the δ(t) response by using the Floquet theory [4.4] and multi-term harmonic balance method
(MHBM) [4.5]. Their work was recently refined and extended to helical gears in our work (refer to Chapters 5 and 6) where closed-form solutions of δ(t) for a SDOF system
are derived under frictional excitations. While the equal load sharing assumption [4.4-4.6]
yields simplified expressions and analytically tractable solutions, they do not describe
realistic conditions. This particular deficiency has been partially overcame in Chapters 2
and 3 where we proposed a multi-degree-of-freedom (MDOF) model with realistic k(t) and sliding friction functions. However, we utilized numerical integration and fast
Fourier transform (FFT) analysis methods in Chapters 2 and 3 that are often computationally sensitive. Hence, a semi-analytical algorithm based on MHBM is highly desirable for quick constructing frequency responses without any loss of generality.
Recently, Velex and Ajmi [4.12] implemented a similar harmonic analysis to approximate the dynamic factors in helical gears (based on tooth loads and quasi-static transmission errors). Their work, however, does not describe the multi-dimensional system dynamics or include the frictional effect, which may lead to “multiplicative” terms as described earlier.
96 The prime objective of this chapter is thus to extend the above mentioned publications [4.4]. In particular, we intend to develop semi-analytical harmonic balance solutions to the 6DOF spur gear model of Chapter 2 with realistic k(t) and sliding friction functions. The example case used for this study is the unity ratio NASA spur gear (with tip relief); refer to Table 2.2 of Chapter 2 for its parameters. Key assumptions include: (i) the pinion and gear are rigid disks; (ii) vibratory motions are small in comparison to the nominal motion; this would lead to a linear time-varying model; (iii) Coulomb friction is assumed with a constant µ though sign is reversed at the pitch point; (iv) when the torsional component is dominant over the translational component of δ(t) for the 6DOF model of Chapter 2, the harmonic solutions of the SDOF system could be extended to predict translational responses in the line-of-action (LOA) and off-line-of-action (OLOA) directions. Note that semi-analytical method analyzes the 6DOF system as a 5DOF model as it calculates the δ(t) and not absolute angular displacements θp(t) and θp(t); all 6 motion terms are determined in the numerical method.
4.3 Semi-Analytical Solutions to the SDOF Spur Gear Dynamic Formulation
Consider the example case with only the mean load Te, i.e. T(t) = 0 including
ε ()t = 0. Equation (4.1) can be rewritten over the mesh cycle 0 ≤ tt≤ c as follows:
97 ⎡⎤ mteδδδ()++⎣⎦ cttktt00 () () () ()[ 1 ++ EEt 13] ⎡⎤L , (4.2a) ++⎡⎤ct ( )δδ ( t ) kt ( ) ( t ) 1 ++⋅−= E Et sgn( tAP ) F ⎣⎦11⎢⎥() 23 e ⎣⎦⎢⎥Ω pbpr
⎡⎤ E1 =++µλ⎣⎦()LJrLJrJXA g bp YC p bg/ b , (4.2b)
⎡⎤ E2 =++µλ⎣⎦LJrXA g bp() L YC Jr p bg/ J b , (4.2c)
E3 =Ωµ prJrJrJbp( g bp − p bg) / b , mJJeeb= / , Feeb= TJ / . (4.2d-f)
Next express Eq. (4.2) in terms of the dimensionless timeτ = tt / c , such that
222 δ ′()τδτδ==dd / ttc () andδ ′′()τδτδ==dd / ttc ():
mtctkEtEfecδτ′′()++[]00 () τδτ ′ () c ()() τδτ[ 1 ++ 130 c () τ] , (4.3a) 2 ++tccc[][]11 (τδ )′ ( τ ) tk ( τδτ ) ( ) 1 ++= Ef 2132 ( τ ) Ef ( τ ) tF ce
f0 (τ )= mod(τ ,1) , ff10(τ )=−=− sgn[ mod(ττ ,1)PP] sgn[ ( ττ ) ] , (4.3b-c)
f201(τ )=−= mod(τττττ ,1)sgn[] mod( ,1)P ff ( ) ( ) , τ PAP= L / λ . (4.3d-e)
Each periodic function, ki (τ ) , ci (τ ) , fi (τ ) or δ ()τ now has a period of T = 1; Figure
4.1(b) shows the typical f0 ()τ , f1()τ and f2 ()τ functions, which describe the periodic moment arm and sliding friction excitations for the example case.
98 4.3.1 Direct Application of Multi-Term Harmonic Balance (MHBM)
Define the Fourier series expansions of the periodic ki (τ ) and ci(t) in Eq. (4.3) up
to N mesh harmonics as follows, where ωn = 2π n (in rad/s) and n is the mesh order.
NN kAikikinn()τ =+0 ∑∑ A cos(ωτ ) + B kinn sin( ωτ ). (4.4a) nn==11
1 1 1 Akd= ()τ τ , Ak= 2cos()ω ττ d, B = 2sin()kdω ττ. (4.4b-d) ki0 ∫0 i kin∫0 in n kin∫0 in n
NN ckIAAi(τ )==++ 2ζ i ( τ ) e ci0 ∑∑ cin cos( ωτ n ) B cin sin( ωτ n ) . (4.5a) nn==11
1 1 1 Acd= ()τ τ , Ac= 2cos()ω ττ d, B = 2sin()cdω ττ. (4.5b-d) di0 ∫0 i cin∫0 in n cin∫0 in n
The fi(τ) functions could be expanded explicitly as shown below:
11N f0 ()τ =−∑ sin(ωτn ), (4.6a) 2 n=1 nπ
NN4cos⎡⎤ωτ − 1 4sin(ωτnP) ⎣⎦( nP) f1()τ =− 1 2τωτωτPn −∑∑ cos( ) + sin( n ), (4.6b) nn==11ωωnn
⎧ N ⎫ ⎪∑ ⎣⎦⎡⎤1sincoscos()−−ωτnP() ωτ nP () ωτ nP ωτ n +⎪ 14⎪ n=1 ⎪ f ()ττ=−2 + . (4.6c) 2 P 2 ⎨ N ⎬ 2 ωn ⎪ ⎪ ∑ ⎣⎦⎡⎤ωnPτωτωτωωτcos()() nP−− sin nP 0.5 n sin( n ) ⎪⎩⎭n=1 ⎪
99 Finally, assume that the periodic dynamic responseδ ()τ is of the following form:
NN δ ()τωτωτ=+AAδδ0 ∑∑nn cos( ) + B δ nn sin( ). (4.7) nn==11
Substitute Fourier series expansions of Eqs. (4.4-4.7) into Eq. (4.3) and balance
the mean and harmonic coefficients of sin(ωnτ ) and cos(ωnτ ) . This converts the linear periodic differential equation into easily solvable linear algebraic equations (as expressed
below) where Kh is a square matrix of dimension (2N+1) consisting of known
coefficients ofki (τ ) , ci (τ ) andfi (τ ) . By calculating the inverse of Kh , the 2N+1 Fourier coefficients of δ ()τ could be computed at any gear mesh harmonic (n).
⎡⎤Aδ 0 ⎡Fe ⎤ ⎢⎥A ⎢ 0 ⎥ ⎢⎥δ1 ⎢ ⎥ ⎢⎥B ⎢ 0 ⎥ ⎡⎤δ1 ⎣⎦Kh ⎢⎥= ⎢ ⎥ . (4.8) ⎢⎥... ⎢... ⎥ ⎢⎥A ⎢ 0 ⎥ ⎢⎥δ N ⎢ ⎥ ⎣⎦⎢⎥Bδ N ⎣ 0 ⎦
4.3.2 Semi-Analytical Solutions Based on One-Term HBM
Next, we construct one-term HBM [4.13] solutions to conceptually illustrate the method. Set the harmonic order N = 1 (only the fundamental mesh, in addition to the mean term) in Eqs. (4.4-4.8) and balance the harmonic terms in Eq. (4.3). This leads to a
100 Kh matrix of dimension 3. Three of its typical coefficients are given as follows and the rest could be found in a similar manner:
1 ⎛⎞tAck11 A f 21 E 3+− B k 11 B f 11 E 2(/)22 tB ck 01 E 3π +++ A k 00 A k 10 A k 11 A f 11 E 2 K = (4.9a) h11 ⎜⎟ 2 ⎝⎠++222AEtAEkckkfckfckf00 1 00 3 + AAE 10 10 2 + tAAE 10 20 3 + tBBE 11 21 3
K =++A A AEAAE + + tAAE + AAE hkk21 01 11 01 1 kfckfkf 10 11 2 10 21 3 11 10 2 (4.9b) ++tAck11 A f 20 E 3() tA ck 01 E 3 / 2
Khkfkckf31=++++++BAE 11 10 2 B 11 tBA 11 20 E 3 ABE kfckf 10 11 2 tABE 10 21 3 BE k 01 1 B k 01 ttBE (4.9c) −+cckAE 01 3 π k 00 3 2
The Fourier series coefficients of δ ()τ are then obtained by inverting Kh . Figure
4.3 shows that one-term HBM solution predicts the overall tendency (mean and first
harmonic) fairly well when compared with numerical simulations at Tp = 550 lb-in and
Ω p = 500 RPM. This confirms that the one-term HBM (and likewise the MHBM) approach coverts the periodic differential Eq. (4.3) with multiple interacting coefficients into simpler algebraic calculations that are computationally more efficient than numerical integrations and subsequent FFT analyses. Thus, the semi-analytical solution provides an effective design tool. Also, most coefficients of Eq. (4.9) show side-band effects that are introduced by k(t) (or c(t)) and the fi(t) functions.
101 x 10-4 10.4
10
(in) 9.6
9.2
0 1 2 t/t c (a) (in)
(b)
Figure 4.3 Semi-analytical vs. numerical solutions for the SDOF model, expressed by Eq.
(4.3), given Tp = 550 lb-in, Ω p = 500 RPM, µ = 0.04. (a) Time domain responses; (b) mesh harmonics in frequency domain. Key: , , numerical simulations; , , semi-analytical solutions using one-term HBM; , , semi-analytical solutions using 5-term HBM. 102 4.3.3 Iterative MHBM Algorithm
When N ≤ 5, we can utilize a symbolic software [4.14] to balance multiple
harmonic terms and calculate Kh . However, the computational cost involved with each
3 element of Kh increases by N due to the triple multiplication of periodic coefficients in
Eq. (4.3). Consequently, for higher N (say >5), a direct computation of Kh becomes inefficient and thus inadvisable. Instead, we apply a matrix-based iterative MHBM algorithm [4.15, 4.16]. First, define variables Ω and ϑ: