Helical Gears, Effects of Tooth Deviations and Tooth Modifications On Load Sharing and Fillet Stresses
Helical Gears, Effects of Tooth Deviations and Tooth Modifications On Load Sharing and Fillet Stresses
Helical Gears, Effects of Tooth Deviations and Tooth Modifications On Load Sharing and Fillet Stresses
net/publication/275382387
CITATIONS READS
26 419
3 authors:
L. Cloutier
Laval University
39 PUBLICATIONS 365 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Claude Gosselin on 11 August 2016.
1 Introduction the reference section 关2–9兴. One or more of these aspects, under
specific meshing conditions, are explored in these studies. For
With the industrial revolution, early gears made of wood were
instance, Wagaj et al. 关5兴 present a parametric study of the effects
transformed into modern gears, offering remarkable tooth shapes. of profile modifications. The authors illustrate the influence of a
Nowadays, gear design is largely standardized. Organizations perfectly assembled gear set over one mesh cycle.
such as AGMA or ISO supply standards that establish accuracy The main objective of this paper is to illustrate the need for the
classification, rating formulas, etc. Common gear practice also careful study of variation in design parameters, especially the po-
offers common design techniques that are largely employed, but tential consequences, using specific case studies. The intention
not always fully or explicitly normalized, and which refer to past here is not to portray the effects over the complete meshing period
experience as a guide. Tooth tip relief is a good example. In this of a tooth pair, which could be cumbersome or could average
context, the gear design process could appear to be a straightfor- some local upshots. Instead, it is preferred to select a particular
ward and easy task, with no deficient result likely to occur. In position that includes all possible contact types at the same time.
reality, however, the guidelines provided by the standards should From our point of view, this will offer a clearer understanding of
be combined with complementary experience and analysis, other- the phenomenon and be more condusive to general conclusions.
wise poor gear sets may result. The gear designer should be able The express model 关10兴 is used for all simulations. This model
to predict the performances of a gear pair from its accuracy grade combines the finite strip method with a pseudo 3-D model of the
and assembly quality. tooth base solved with the finite differences method to establish
The ANSI/AGMA ISO 1328-1 Standard for Cylindrical Gears tooth displacement and fillet stresses. This representation is
关1兴 sets 13 gear accuracy grades with allowable helix slope and coupled to a cell discretization of the contact area in order to
form deviations. While these inspection items are not individually determine the contact pressure distribution.
subject to mandatory tolerances, they are recognized as having a
significant influence on the performance of a gear pair. A com-
parative investigation is conducted on the effects of helix slope 2 The Express Model Brief Description
and form deviation tolerances, as specified by grades 5 and 7 for The load distribution between and along meshing teeth is ruled
this standard. by the stiffness of the contacting tooth pairs. Equations 共1兲 and 共2兲
The assembly precision also has a strong influence on gear express the displacement compatibility condition between mesh-
behavior. To reduce the effects of assembly deviation such as ing tooth pairs, and the torque equilibrium between the loaded
misalignment, flank crowning and tip relief are often recom- surfaces and the external torque, respectively.
mended 关2,3兴. The resulting consequences of such tooth modifi- ជi + ⌺
ជc + ⌺
ជe + ⌺
ជg = ⌺
ជa
cations are studied when applied on a misaligned gear set. The ⌺ kc kc kc kc 共1兲
冉冕 冊
investigation is conducted from the perspective of static load shar-
ing and its influence on surface contact and fillet bending stress.
Various papers dealing with assembly errors, profile modifica- 兺D
p共兲储rជ储共nជ · ជt 兲d =T 共2兲
tions, and errors have been published, many of which are listed in Lj D
Equation 共1兲 is the contact condition for a point with coordinates
Contributed by the Power Transmission and Gearing Committee of ASME for
kc in the plane of contact—when Eq. 共1兲 is satisfied, point kc
publication in the JOURNAL OF MECHANICAL DESIGN. Manuscript received November 8, belongs to the loaded surfaces. This equation regroups five angu-
2004; final manuscript received May 26, 2005. Review conducted by Teik C. Lim. lar terms:
444 / Vol. 128, MARCH 2006 Copyright © 2006 by ASME Transactions of the ASME
ជ a:
⌺ ⫽ the angular approach of the mating teeth;
ជ c:
⌺ ⫽ the contact deformation;
ជ e:
⌺ ⫽ the tooth bending and shearing deflection;
ជ i:
⌺ ⫽ the initial profile separation;
ជ
⌺ g: ⫽ the surface deformation due to frictional traction;
Fig. 1 Displacement-stress model of a helical gear tooth Fig. 2 Pressure cells in contact plane
Equation 共2兲 joins the individual torque produced by the load Over the segments the 2-D governing equations are solved with
distribution p共兲 along each of the D tooth pairs in mesh and the the finite difference method 共FDM兲. The main advantage of the
external torque 共T兲. Terms nជ , 储rជ储 and ជt are the unit normal to the FDM is the mesh where irregular grids are very easy to create.
contact plane, the radius of the meshing position, and a unit vector Since in the FDM, nodes are not attached to elements, there is no
in a plane perpendicular to the axis of rotation, respectively. need to impose common node sharing as with FEM meshing.
The simultaneous solution of the preceding equations leads to Thus, the mesh density can be refined rapidly in regions of stress
the load distribution along each meshing surface. Over the years, concentration.
numerous models have been presented 关2,3,6–9,11–21兴. Due to Contact analysis is a major concern in this kind of simulation.
the 3D integration implied and the nonlinear behavior of the con- The rolling-sliding load on meshing tooth surfaces is composed of
tact zones, most of the calculations were based on analytical rep- normal pressure and shear traction. An accurate and economic
resentations, often including important approximations 关2,6–12兴 way to handle such a problem is the numerical algorithm based on
or on different numerical approaches such as the finite element a discretization of the contact plane in pressure shear cells—
method 共FEM兲 关3,9,13–21兴. The express model was developed in influence functions for uniform pressure and shear tractions acting
order to eliminate the loss of precision associated with common on a rectangular zone are introduced in an iterative calculation
simplifications, while reducing the computer time usually attrib- scheme. The contact cells easily cope with any surface irregular-
uted to numerical simulation of contact gear problems. Equations ity, without increasing time consumption. The basis of the method
共1兲 and 共2兲 can be solved concurrently for up to five tooth pairs. is given in Hartnett’s work 关23兴. The approach has been applied to
The finite strip method 共FSM兲 introduced earlier in the gear the gear contact problem by a few authors 关15,17,18,24兴. Figure 2
domain 关22兴 establishes the tooth displacement from its neutral schematically illustrates the procedure. Directions and run
surface represented by a clamped 2-D strip. This approach con- along and across the contact line, respectively. This contact model
siderably reduces the node number while ensuring a precision in was also validated in a previous paper 关10兴 by comparison to
bending and shearing displacement within 5% compared to a 3-D Hertz’s prediction: for a spur gear set, the agreement on the maxi-
FEM analysis. mum contact pressure was higher than 95%.
The FSM clamps the tooth at the root radius. Therefore, to work Figure 3 graphically presents the express model.
out the complete tooth displacement, the strips are combined to a
tooth base discretization. This base model calculates the displace- 3 Helix Slope and Form Deviations: Effects of Toler-
ments of the tooth support as well as the fillet stresses. Figure 1
ance Margins for ANSI/AGMA ISO 1328-1 Grades 5
illustrates the resulting representation of a helical tooth used to
ជ e兲. and 7 on Contact and Fillet Stresses
calculate the complete tooth deflection 共⌺
Figure 1 shows that the tooth base is divided into a series of 3.1 Helix Slope Deviation „fh…. Real-life gear teeth must
segments of equal width. This strategy is comparable to the one comply with standardized tolerance or deviation margins that
employed in other approaches such as the models of Elkholy 关11兴 specify the grade to which the gears belong. The deviations are
or Smith 关9兴, and thus offers similar speed gain. Nevertheless, the normally measured in the direction of the transverse base tangent
loss of physical continuity usually associated with such a tech- and apply to an evaluation length that can be shorter than the tooth
nique is eliminated in the model—the junction between the base face width. In this study, the evaluation length is equal to the face
and the strips that are lengthwise continuous links the response of width.
the base segments. Since at any position along the tooth, there is a A 200 N · m torque is applied to the helical pinion of Table 1.
direct relation between load distribution and displacement, the The studied meshing position is shown schematically on the pin-
discrete load acting on each tooth base segment is calculated from ion member in Fig. 4, where two tooth pairs mesh simultaneously.
the FSM displacement at the junction line. Reference 关10兴 com- Table 2 gives the contact line endpoint coordinates. Figure 4 also
pares and validates the model with 3-D FEM analysis. presents the tooth meshing sequence: tooth pair 0 is the main
Helical gear
pinion/gear
共R , X3兲
共mm兲 Endpoint I Endpoint II
Fig. 3 Flowchart of the express model Table 3 Displacement-stress and contact cell model
parameters
meshing tooth pair, while tooth pair ⫺1 follows. Value Pinion Wheel
The theoretical contact lines observed on the pinion teeth ex-
tend, for tooth pair ⫺1, from the fillet radius to a radius just above D⬘ 30.00 mm 30.00 mm
M 15.00 mm 15.00 mm
the pitch circle and, for tooth pair 0, from a radius above the pitch
l⬘ 14.09 mm 13.88 mm
circle to the tip circle 共Fig. 4兲. Table 3 presents the parameter
 11.70 deg 15.11 deg
values introduced in the express model 关10兴. Finite 0.23 mm 0.23 mm
The measurement of helix slope deviation is obtained between differences
two helix traces that intersect the mean helix trace at the endpoints increment size
of the evaluation length 关1兴. Table 4 lists the tolerance values for No. of finite strips 14 14
grades 5 and 7 from their tooth numbers. 共5 nodes per
strip-radial
For a given grade, the combined deviation of a tooth pair is a direction兲
function of pinion and gear tooth accuracy. In order to limit the Contact plane 18 cells along, 11 cells across contact line
combined pinion and gear deviation range, individual deviations
serve here for the purpose of analysis. Table 5 lists all the condi-
tions modeled with helix slope deviation. For comparison pur-
poses, the results of the unmodified theoretical gear set are also
presented.
The results obtained for the operating conditions of Table 5 and
the perfect gear pair are shown in Figs. 6–16.
The helix slope deviation of a tooth, if not compensated by an
inverse deviation in the mating tooth, results in a shift of the
contact zone along the theoretical contact line. In case A, the
contact zone on tooth pair ⫺1 shifts toward endpoint I, Fig. 1, and
toward point II for tooth pair 0. For case B, the contact zone
moves toward point II for tooth pair ⫺1 and toward point I for
tooth pair 0. In short, for both cases, the deviations will result in a
load reduction of the tooth pair where the deviation is negative,
and in a load increase for the tooth pair with a positive deviation.
The first series of results, Figs. 6–10, shows the pressure dis-
tribution in the contact plane of the active tooth pairs.
For the theoretical gear pair 共Fig. 6, no deviation兲, pressure Fig. 6 Pressure distributions on contact plane—no helix slope
peaks related to tip edge contact can be observed. For tooth pair deviation
⫺1, the pressure elevation is clearly visible; for pair 0, the phe-
nomenon is of lesser importance. The real contact pressure value
at tooth edge cannot be predicted exactly, since the model is based The assumption made in this study is that a similar convergence
on a contact cell discretization 关10兴 where the solution gives an could be reached using the cells in edge contact simulation. In
average pressure for each contact cell, and thus depends on cell fact, even if a convergence demonstration is not presented, it is
size and the position of the pressure center relative to the tooth assumed that the 80 m width cells 共ten divisions along or the
edge. On the other hand, Bhushan 关26兴 suggested that, when the 0.8 mm width of meshed contact plane兲 employed here produce a
contact cells are small enough, the predicted pressure on contact- sufficient precision level in the flank limit vicinity. To the authors’
ing asperities of rough surfaces converge to the real amplitude. knowledge, the current literature does not offer experimental
works dealing with such a particular contact type that are suitable
to validate this constant pressure cells approach.
The effects of helix slope deviation on contact pressure are
observed in Figs. 7–10. For the studied meshing position, case A
Pair ⫺1 Pair 0
Case A Grade 5 − ⫹
Grade 7 − ⫹
Case B Grade 5 ⫹ −
Grade 7 ⫹ − Fig. 7 Pressure distributions on contact plane—case A, grade
5
entire contact line. On the other hand, when the first contact point
⫺1. is at the tip radius, such as tooth pair 0 共Case A兲, the perturbation
The analysis of Table 6 and Fig. 11 reveals the importance of is localized. In others words, the effect of deviations depends on
the load increment and the influence of the helix slope deviation, the nature of the contacting surfaces.
depending on the position of the contact line. With a first contact As Table 6 clearly shows, load elevation for grade 7 is almost
point moved to a position close to the pitch radius, for tooth pair twice that of grade 5, and the deviations double from grade 5 to 7.
⫺1 共case A兲, the helix slope deviation influence spreads over the The results in Figs. 12–16 present the maximum principal tooth
Fig. 13 Maximum principal pinion fillet stress distributions— Fig. 15 Maximum principal pinion tooth fillet stress
case A, grade 5 distributions—case B, grade 5
Theoretical Grade 7
value Grade 5 Grade 7
共N/mm兲 共N/mm兲 共N/mm兲 Grade 5
Theoretical Grade 7
value Grade 5 Grade 7
共MPa兲 共MPa兲 共MPa兲 Grade 5
Theoretical Grade 7
value Grade 5 Grade 7
共N/mm兲 共N/mm兲 共N/mm兲 Grade 5
clearly visible. In addition, in Figs. 19共b兲 and 20共b兲 共grade 7兲, the
elastic deformation fails to produce contact in some areas of the
tooth where load distribution falls to zero.
In the case of grade 5, the contact areas where the contact
pressure is reduced sustain an average load of 128 N / mm, equiva-
lent to 21.5% of the mean peak value of 594 N / mm and 37.0% of
the theoretical 共unmodified teeth兲 average value. Moreover, from
Table 9, one can see that the maximum load value nearly doubles
when the grade changes from 5 to 7, where the deviations are
doubled.
Helix form deviations also have a visible effect on fillet stress,
as shown in Figs. 21 and 22 where the maximum fillet stress
Fig. 19 Pressure distributions in the contact plane—grade 7 follows the evolution of helix form deviations. However, the high-
est values remain under acceptable levels, while the intensity
on tooth pair 0, while tooth pair ⫺1 remains unaffected. reached by the contact pressures 共Figs. 18 and 19兲 was significant.
Table 9 lists, for tooth pair 0, the maximum load increase for Table 10 lists the increases in maximum fillet stress. These
the modified teeth; the mean value is calculated along the contact results are averaged values calculated at the peak location on the
line of the unmodified tooth. deviated tooth flanks and at the corresponding position on the
The importance of load increase under helix form deviation is theoretical unmodified pinion tooth 关Fig. 11共b兲兴. The values in
Table 10 show a much weaker influence of the helix form devia-
tions on fillet stress than on contact pressure, since fillet stress
distribution is a response to flank loading, which depends on the
shape, amplitude, and position of the helix form deviation. Con-
sequently, helix form errors will affect contact surface pressure
Fig. 20 Load distributions on simultaneous contact lines „un- Fig. 21 Maximum principal pinion fillet stress distributions—
modified teeth and grades 5 and 7 with undulated profiles… grade 5, deviated tooth
Fig. 26 Pressure distributions on the contact plane—crowned, Fig. 28 Pressure distributions on the contact plane—crowned
misaligned gear set and tip relieved, misaligned gear set
Crowned
Theoretical Tip and tip
value Crowned relieved relieved
共N/mm兲 共N/mm兲 共N/mm兲 共N/mm兲
Crowned
Theoretical Tip and tip
Value Crowned relived relived
共MPa兲 共MPa兲 共MPa兲 共MPa兲
5 Conclusion
The first part of this paper illustrates the influence of helix slope
and form deviations. Longitudinal flank crowning and radial tip
relief modifications were studied.
Helix slope deviation and misalignment of a gear set have com-
parable consequences: the first contact point is driven to an end-
point of the theoretical line of contact, which results in localized
overloading. As a consequence, contact pressure is increased lo-
cally, which generates a corresponding increase in bending stress
Fig. 33 Maximum principal pinion fillet stress distributions— in the fillet zone.
tip relieved, misaligned gear set Flank crowning and tip relief modifications are then introduced;
results show that they should always be considered in the design
process and should be used in concert, since partial correction
pressure. For example, the fillet stress reduction obtained is on the may cause load increase rather than decrease: when used sepa-
order of 9%, which is an improvement. It is also worth emphasiz- rately on a misaligned cylindrical gear set, flank crowning, while
ing that the stress levels always stayed within a normal range, reducing contact stress at tooth ends, will amplify the overload at
whereas the contact pressure amplitudes attained dangerous inten- tip edge. Conversely, tip relief will overload tooth ends and dis-
sities. charge the tip edge.
While helix slope deviations and misalignment errors have
similar consequences, helix slope deviations are potentially more
dangerous since opposite errors on meshing teeth may occur. This
condition results in load transfer from one tooth pair to another.
For example, for a grade 5 gear pair with helix slope deviation,
our results show that load increase is more than 70% with fillet
stress increases of more than 40%. A comparative analysis of
grades 5 and 7 reveals values virtually proportional to the ampli-
tude of the error.
The helix form deviation, modeled as a sine wave superim-
posed on a theoretical tooth surface, presents the same direct re-
lationship between the amplitude of the error and the load state of
the teeth.
Therefore, it can be concluded that reducing or increasing the
quality class of a gear pair will proportionally influence the gear
performance as well as its life.
The analysis presented here also indicates, for all cases, that
contact perturbations led to very high or even excessive contact
pressure, while the corresponding fillet stress remained moderate.
Comparable observations have also been made by other investi-
gators 关5兴. Thus, we can assume that any contact disturbance will
primarily damage the contacting surfaces and might have no sig-
nificant effect on the bending life.
Nomenclature
D ⫽ number of tooth pairs in mesh
D⬘ ⫽ width of flexible support under the tooth
Lj ⫽ length of contact line j
M ⫽ height of flexible support under the tooth
T ⫽ applied torque
Fig. 34 Maximum principal pinion fillet stress distributions— X ⫽ addendum modification coefficient
crowned and tip-relieved, misaligned gear set l⬘ ⫽ height of the tooth base model