Determining The Shaper Cut Helical Gear Fillet Profile: George Lian Management Summary
Determining The Shaper Cut Helical Gear Fillet Profile: George Lian Management Summary
Management Summary
This article describes a root fillet form calculating method for a helical gear generated with a shaper cutter. The shaper cutter con-
sidered has an involute main profile and an elliptical cutter edge in the transverse plane. Since the fillet profile cannot be determined
with closed-form equations, a Newton’s approximation method was used in the calculation procedure. The article also explores the
feasibility of using a shaper tool algorithm for approximating a hobbed fillet form. Finally, the article discusses some of the applications
of fillet-form calculation procedures, such as form diameter (start of involute) calculation and finishing stock analysis.
Introduction
Analytical methods for determining the gear fillet pro- generated external or internal helical gear can be calculated.
file (trochoid) have been well documented. Khiralla (Ref. 1) A shaper tool generating a gear can be visualized as a gear
described methods for calculating the fillet profile of hobbed set meshing with zero backlash. The algorithm in this article is
and shaped spur gears. Colbourne (Ref. 2) provided equations based on a shaper tool in tight mesh with a semi-finished helical
for calculating the trochoid of both involute and non-involute gear. The semi-finished gear geometry was used for calculation
gears generated by rack or shaper tools. The MAAG Gear because the shaper tool used as the semi-finishing tool is usually
Handbook (Ref. 3) also provided equations for calculating the one that generates the trochoid. However, if the shaper cutter
trochoids generated with rack-type tools that have circular tool is the finishing tool, the algorithm presented will also work by
tips. Vijayakar, et al. (Ref. 4) presented a method of determining letting the finishing stock equal zero. The trochoid of a spur gear
spur gear tooth profiles using an arbitrary rack. The above men- can also be calculated by letting the helix angle equal zero.
tioned are only samples of many published works. However, the The shaper tool used in this algorithm may have a different
method for determining the trochoid of a helical gear generated reference normal pressure angle than that of the gear. A neces-
with a shaper tool is not widely published. This article presents sary condition for a shaper tool to generate the correct involute
an intuitive algorithm where the fillet profile of a shaper-tool- profile on a gear is that both the tool and the gear must have
ww w . p o w e r t r a ns m i s s i o n . c o m • w w w . g e a r t e chnology.com • G E A R T E C H N O L O G Y • S E P T E M B E R / O C T O B E R 2 006 57
� inv is the involute function of an angle
inv α = tan α – α
Standard transverse pressure angle of tool, φ0 Base circular thickness of semi-finished gear, sbr (in.)
tan φn0 sr
φ0 = arctan ( ) (1) sbr = 2rbr ( ± inv φr ) (14)
cos ψ0 2r
Standard reference pitch radius of tool, r0 (in.) Center of Tool Tip on a Shaper Tool
n0 A shaper tool for gear semi-finishing usually has protuber-
r0 = 2P (2) ance. It generates undercut on a gear, so that the finishing tool
nd0cosψ0 only needs to machine the involute profile of the gear. To obtain
Base radius of tool, rb0 (in.) the designed amount of protuberance on a shaper tool, the tool
tip is made tangent to the involute profile that is temporarily
rb0 = r0cosφ0 (3) formed by increasing the shaper tooth thickness to include the
protuberance (see Fig. 3). The tangent point, common to the
Reference transverse circular thickness of tool, S0 (in.) tool tip and the involute profile, will be referred to as the profile
sn0 tangent point, P0 . When the temporarily formed involute profile
s0 = cosψ (4) is removed, the shaper tool will have the designed amount of
0 protuberance.
Transverse base pitch of tool, pb0 (in.) The shaper tool tip is also made tangent to the outside diam-
2πrrb0 eter of the tool (see Fig. 4) so that the transition from the outside
pb0 = n (5) diameter to the tool tip will be smooth. The common tangent
0 point on the shaper tool tip and the outside diameter of the tool
Normal base pitch of tool, pnb0 (in.) will be referred to as the end tangent point, E0.
πcosφn0 The following are the required data for calculating the cen-
pnb0 = (6) ter of the shaper tool tip:
Pnd0
Base helix angle of tool, ψb0 da0 is the outside diameter, tool (in.)
pnb0 sb0 is the base circular thickness, tool (in.)
ψb0 = arccos( ) (7) ρ0 is the tool tip radius (in.)
pb0
δ0 is the protuberance (in.)
Base circular thickness of tool, sb0 (in.) ψb0 is the base helix angle, tool
s0 ψ0 is the reference helix angle, tool
sb0 = 2rb0 ( + invφ0) (8)
2r0
The base circular thickness of the involute profile, formed
where
cos θPn0 0
P0 = S0 + ( ρ0 , ρ0 sin θPn0 ) (17)
cos ψ0
En0
where
En0
θPn0 is the auxiliary angle that locates P0. The angle is
Figure 4—End of tool tip (with helix angle exaggerated).
measured in the normal plane, clockwise from the
horizontal axis of the tool tip. θPn0 will usually have a The tool radius to end tangent point, rE0 (in.), is
negative value.
rE0 = || E0 || (24)
The tool radius to profile tangent point, rP0 (in.), is
The following are conditions for the tool tip to position properly
rP0 = || P0 || (18) on a shaper tool tooth:
The transverse pressure angle, φP0, at P0 is
r 1) The profile tangent point, P0, on the tool tip must also be
φP0 = arccos( rb0 ) (19)
P0 a point on the involute profile that includes the protuberance,
thus
The tangent angle, αP0, at P0 (the derivation of Equation 20 is
given in Appendix A) is π
αP0 + φP0 – ζP0 – 2 = 0 (25)
–cos ψ0
αP0 = arctan ( ) (20) 2) The angle, ζP0, subtended by one half of the transverse cir-
tanθPn0
tanθ
cular thickness of the involute curve (include the tool protuber-
The angle between the y-axis and the radius to the profile tan- ance) at P0, must equal the angle formed by the y-axis and the
gent point, ζP0, is line connecting the center of the tool to P0.
sb0_pr
ζP0 = – inv
invφP0 (21) xP0
2rb0 ζP0 – arcsin ( r ) = 0 (26)
P0
ww w . p o w e r t r a n s m i s s i o n . c o m • w w w . g e a r t e chnology.com • G E A R T E C H N O L O G Y • S E P T E M B E R / O C T O B E R 2 006 59
∂ffi fi(X + ∆X
Xj) – fi(X)
Equations 25–28 must all be satisfied for the tool tip to be ≈ (36)
∂xj δδxxj
correctly positioned on a shaper tool tooth. The variables to be
determined are rS0, λS0, θPn0 and θEn0. Since the systems of the where
equations are transcendental and cannot be solved directly, the i is the ith row of the Jacobian matrix
Newton’s method is used to calculate the roots for Equations j is the jth column of the Jacobian matrix
25–28. ∆XXj is a vector with its jth element equal to the jth element
Solving the System of Non-linear Equations of the current Newton’s step, δX, and all remaining
for Center of Tool Tip elements equal 0
For simplicity, rewrite Equations 25–28 as generic vector
equations in the form For each iteration, the sum of the absolute values of the
functions (errors) is calculated.
F(X) = 0 (29)
4
ERR(X1) = Σ |f (X1)| (37)
i=1 i
where
F(X) = (f1(X), f2(X), f3(X), f4(X))T The Newton’s iteration procedure is terminated when the error
= (Eq. 25, Eq. 26, Eq. 27, Eq. 28)T (30) (see Eq. 37) becomes smaller than a predetermined tolerance, or
when a predetermined number of iterations has been reached.
0 = (0,0,0,0)T (31) The Newton’s iteration procedure is described below:
X = (x1, x2, x3, x4)
T
= (rS0, λS0, θPn0, θEn0)T (32) 1) Select a set of initial guess values for the new root, X1. The
following are the suggested values:
The Newton’s iteration equation (Ref. 6) is written as da0
x1(rS0) = – ρ0
2
X1 = X + δX
δX (33) x2(λS0) = 0.0175
x3(θPn0) = –φn0
where δX satisfies the following system of linear equations x4(θEn0) = –1.4835
J • δX = – F(X) (34) 2) Select the initial Newton’s steps, δX. The following values
work satisfactorily:
where
δX = (0.01, 0.01, 0.01, 0.01)T
X1 is the vector of the new roots for the next iteration
X is the vector of current roots
δX is the vector of Newton’s steps for the next iteration 3) Evaluate the system of non-linear equations (see Eq. 30) at
J is the Jacobian matrix the new root, F(X1).
4) Calculate the error ERR(X1) (see Eq. 37).
where 5) The iteration is terminated, if ERR(X1) ≤ 10–10 or if a pre-
determined number of iterations (30 should be sufficient)
∂∂ff1
∂x1
∂x
∂f1
∂x2
…
∂f1
∂x4
have been reached. Otherwise, continue with the next
step.
6) Save the new roots as the current roots, so that a new set of
roots can be calculated
∂∂ff2 ∂ff2
… (35)
…
7) Calculate the Jacobian matrix, column by column, starting
…
…
ww w . p o w e r t r a n s m i s s i o n . c o m • w w w . g e a r t e chnology.com • G E A R T E C H N O L O G Y • S E P T E M B E R / O C T O B E R 2 006 61
�
r X0 rg0 rS0
rX0
rg0
rS0
0
G G
G0
rg '
X0
'
S0
rg TX0
Normal
X0 Normal (after rotation)
S0 X0
Figure 6—The arbitrary point X0 and the normal after rotating the shaper tool Figure 7—The arbitrary point X0 on the tool tip and its conjugate point TX0
for an angle –ξ0. on the trochoid.
Semi-finished
involute profile
Trochoid
P0
Finished involute
Invo 3 -
Pro ute
E0 profile
file
e
(SOI)
2- Conjugate pt.
ne l
Zo iptica of P0
p
Zone 1 - Ell ol Ti
Cutter OD To Conjugate pt.
of E0
x M(–ξ0)X0
X0' = M(–ξ (50)
ξ0 = arcsin( rG0 ) (49)
G0
where
Note: The angle ξ0 may be positive or negative. If G0 is
ϕ
cosϕ ϕ
–sinϕ
M(ϕ)) =
on the left side of the y-axis, xG0 (see Eq. 49) will be negative, M(ϕ (51)
and so will ξ0. On the other hand, if G0 is on the right side of sinϕ
sinϕ cosϕ
cosϕ
the y-axis, ξ0 will have a positive value. M(ϕ) is a rotation matrix. When multiplied to a vector, the
To find the conjugate point of X0, the shaper tool is rotated vector would be rotated an angle ϕ about the origin (0,0). If
from its reference position (see Fig. 6) by an angle, –ξ0. The ϕ > 0, the rotation is counterclockwise. Otherwise, the rotation
arbitrary point X0 will rotate to a new position, X0' is clockwise.
ww w . p o w e r t r a n s m i s s i o n . c o m • w w w . g e a r t e chnology.com • G E A R T E C H N O L O G Y • S E P T E M B E R / O C T O B E R 2 006 63
� the hob to be 0.001901". For a 10,000-tooth shaper cutter,
the difference decreased to merely 0.000011". The difference
between the shaper-generated and the hobbed SOI’s followed a
similar trend. For the 100-tooth shaper cutter, the difference was
SOI-HOB HOB
0.010209" and for 10,000-tooth shaper cutter, 0.000117".
SOI-10000T 10000T
The trochoid curves plotted in Figure 10 show the shaper-
SOI-1000T 1000T
SOI-100T
generated trochoid converging to that of the hobbed one when
100T
the number of teeth in the shaper tool is large (e.g. 10,000).
Applications for the Shaper Tool Algorithm
The shaper tool algorithm can be used in computer-aided
0.0025
gear design and gear tooth modeling as shown in Figure 11. The
algorithm is also useful for calculating the trochoid geometry
for finite element or boundary element analysis. The following
Figure 10—Pinion trochoid (Ex. 3.1.5-AGMA 918-A93) generated with a hob sections describe applications of the shaper tool algorithm in
and shaper cutters.
form diameter calculation and gear finishing stock analysis.
Calculating Form Diameter
�
The form diameter or the start of involute (SOI) of a
finished gear is the gear diameter where the trochoid joins or
intersects the involute profile. When the two curves intersect,
two intersection points may appear to exist. The intersection
point that is closer to the tip diameter of the gear is the SOI.
The other “intersection” point is an artificial one, as the involute
curve has already been truncated at the SOI. When calculating
the SOI of a gear by iteration, it is important to make sure that
the algorithm converges to the SOI. Plotting the trochoid and
the involute profile will provide a visual verification that the
iteration process converges correctly (see Fig. 12).
Figure 11—A 23-tooth internal spur gear model.
The SOI can be calculated by comparing the polar angles of
� a trochoid point and an involute profile point, εtro and εinv respec-
Fillet profile tively, on the same gear diameter (see Fig. 12) (Ref. 9). When
P'tro
the two polar angles become equal, the trochoid and involute
Pinv Involute profile points will coincide, and the gear diameter at the intersection
'
Pinv
point is the SOI. If the two polar angles are unequal, compare
Form diameter (SOI)
Ptro the polar angles for a new set of points at slightly larger or
smaller gear diameter than the current one. Repeat the process
until the two polar angles become equal.
Table 2 compares the calculated SOI’s of the selected
numerical examples in AGMA 918-A93 (Ref. 8) using the
tro
i nv
Table 3—Comparison of Form Diameters of a Pinion Generated with Normal Lead and Short Lead Hobs.
Oper. Cntr. Dist. 14.500 in.
Pinion Hob A Pinion Hob B
Description Unit
Gear Pinion (normal lead) (short lead)
Therefore, hob (A) does not provide the required grinding stock on a class of shaper tool that has an involute main profile and
while keeping the SOI below the required SAP. In order to push elliptical tool tip in the transverse plane. However, the algorithm
the SOI closer to the root diameter, a short lead hob (B) was will also work for a shaper tool with other tool tip geometries,
designed. The hob had a 14.5° reference normal pressure angle. provided the coordinates and the normal of the tool tip profile
The calculated SOI based on the short lead hob (B) was 4.4550", are known.
smaller than the SAP. Hob B provided the required finishing The shaper tool algorithm can also approximate the tro-
stock with satisfactory SOI (see Fig. 13). choid generated with a rack-type tool if the number of shaper
Conclusions tool teeth is large. The numerical examples showed that a
A method for determining the shaper-tool-generated fillet trochoid curve generated with a 10,000-tooth shaper tool can
profile (trochoid) was presented. The method is applicable to approximate that generated with a hob with small error.
both external and internal helical gears. The algorithm is based The algorithm presented in this article does not require the
ww w . p o w e r t r a n s m i s s i o n . c o m • w w w . g e a r t e chnology.com • G E A R T E C H N O L O G Y • S E P T E M B E R / O C T O B E R 2 006 65
2. Colbourne, J.R. The Geometry of Involute Gears, New
York, Springer-Verlag, 1987.
3. MAAG Gear Company Ltd. MAAG gear handbook. Zurich,
SOI-Hob A
SAP Switzerland, 1990.
(Normal lead)
SOI-Hob B
4. Vijayakar, S.M., B. Sarkar and D.R. Houser. Gear Tooth
(Short lead)
Profile Determination from Arbitrary Rack Geometry.
AGMA Technical Paper No. 87 FTM 4, 1987.
TROCHOID-Hob A
5. American Gear Manufacturers Association, AGMA 908-
(Normal lead)
B89: Geometry Factors for Determining the Pitting
Resistance and Bending Strength of Spur, Helical and
Herringbone Gear Teeth (AGMA Information Sheet),
TROCHOID-Hob B
(Short lead)
American Gear Manufacturers Association, Alexandria,
VA, 1989.
6. Press, W., S. Teukolsky, W. Vetterling, and B. Flannery.
Figure 13—Trochoid curves generated with a normal lead and a short-lead
hob.
Numerical Recipes in Fortran 77, The Art of Scientific
Computing, Second Edition. New York, University of
tool and the gear to have equal reference normal pressure angle. Cambridge, 1986, 1992.
Consequently, a trochoid generated with a non-standard cutter 7. Cheney, W., and D. Kincaid. Numerical Mathematics and
such as a short lead hob can also be calculated. Computing. Belmont, CA, Wadsworth, 1985, 1980.
Examples for the form diameter (SOI) calculation and the 8. American Gear Manufacturers Association, AGMA 918-
finishing stock analysis were provided using the shaper tool A93 : A Summary of Numerical Examples Demonstrating
algorithm presented. the Procedures for Calculating Geometry Factors for Spur
A computer program was developed using the algorithm and Helical Gears (AGMA Information Sheet), American
described in this article. The calculated form diameters (SOI’s) Gear Manufacturers Association, Alexandria, VA, 1993.
for both external and internal gears compare well to those 9. American Gear Manufacturers Association. Gear Rating
calculated with other gear software. An internal spur gear was Suite User’s Manual, American Gear Manufacturers
used to verify the shaper tool algorithm. Association, Alexandria, VA, 2003.
Printed with permission of the copyright holder, the Ameri- Differentiating Equation A.1 and Equation A.2 with respect to
can Gear Manufacturers Association, 500 Montgomery the auxiliary angle θXn0 we get
Street, Suite 350, Alexandria, Virginia 22314-1560. sinθXn0
Statements presented in this paper are those of the author dxX0 = –ρ0 dθXn0 (A.3)
cosψ0
and may not represent the position or opinion of the
American Gear Manufacturers Association. dyX0 = ρ0 cosθXn0 dθXn0 (A.4)
Taking an arc tangent on both sides of Equation A.6 completes View "A-A" Normal
the derivation for Equation 20. Transverse plane
A A
ψ0
–cosψ
αP0 = arctan ( ) (A.7) 0
tanθPn0
The normal at the given arbitrary point on a shaper tool tip is
Cylinder of
perpendicular to the tangent. Therefore, the slope of the normal, radius 0
sb0 is the transverse base circular thickness, tool (in.); Substituting Equation A.9 and Equation A.10 into Equation
rb0 is the base radius, tool (in.); A.8 and dividing both sides of the new equation by 2rrg0, we get
sbr is the transverse base circular thickness, semi-finished pg0 sb0 rgr sbr rgr
gear (in.). If shaping is the finishing operation, the = – invφg + ± invφg (A.12)
2rg0 2rb0 rg0 2rbr rg0
base circular thickness for the finished gear should be
using the following established relationships
used; and
rbr is the base radius, semi-finished gear (in.). pg0 pb0
= (A.13)
2rg0 2rb0
The sum of the transverse circular thickness of the tool and
rgr rbr
the gear equals the circular pitch at the generating pitch circle. = (A.14)
rg0 rb0
pg0 = sg0 + sgr (A.9) where
ww w . p o w e r t r a n s m i s s i o n . c o m • w w w . g e a r t e chnology.com • G E A R T E C H N O L O G Y • S E P T E M B E R / O C T O B E R 2 006 67