A Practical Approach For Modeling A Bevel Gear: Brendan Bijonowski
A Practical Approach For Modeling A Bevel Gear: Brendan Bijonowski
A Practical Approach For Modeling A Bevel Gear: Brendan Bijonowski
Printed with permission of the copyright holder, the American Gear Manufacturers Association, 1001 N. Fairfax Street, Fifth Floor, Alexandria, VA 22314-1587. Statements
presented in this paper are those of the author(s) and may not represent the position or opinion of the American Gear Manufacturers Association.
Basic Generation
The majority of all generated spiral bevel gears are manufac-
tured in one of two processes, i.e. — face milling or face hob-
bing; both manufacturing methods have advantages and disad-
Figure 3 Generating triangle in pitch plane for face milling.
vantages. For the purposes of this method, a brief understand-
ing of the generating method utilized during these processes is
necessary to realize the 3D model.
Face milling. The face milling manufacturing method
employs a circular, cup-shaped cutting tool moving in a timed
relationship with the workpiece to roll through the gear blank
and generate an individual slot. The cutter is then withdrawn,
the work is indexed to the location of the next slot, and the pro- Since the lengthwise shape of a face hobbed part creates an
cess repeated. See Figure 3. extended, epicycloidal shape, a little trigonometry is necessary
Figure 3 displays three instances of the cutter as it passes to calculate the discrete points along the path created by the cut-
through the workpiece; these points are at the toe, mean, and ter.
heel of the crown gear. The path of the cutter sweeps a circu- The crown gear tooth count,
lar arc in the lengthwise direction, with the same radius as the (12)
radius of the cutter, Rc. The axis of the crown gear, Cg, is known N
Nc =
as the machine center. The local coordinate system for the pitch sin Γ
plane is located at the machine center. Xpp correlates to X and The lead angle of the cutter,
Ypp lies along the line describing Ao in Figure 2. (13)
The cutter axis, Cc. This location can be found by:
V = Rc cos ψm
(1)
-1
Rc Nc [
ν = sin Am Ns cos ψm ]
(2) The first auxiliary angle,
H = Am − Rc sin ψm (14)
(3)
Cc = (V, H) λ = π − ψ+ ν
2
The sign of the vertical term, V, may either be positive or neg-
ative, depending on the hand of spiral to be modeled. Figure 3
shows a right-hand member (use negative value for V when cal-
culating left-hand members). The cutter sweep angle, θr, is nec-
essary to determine how much rotation is used during the cut-
ting process. When calculating points to model the gear tooth,
these will be the endpoints of the working portion of the cutter
path.
(4)
S = √V2 + H2
(5)
[2 2 2
] [ 2 2 2
θr = cos-1 S + Rc − Ao − cos-1 S + Rc − Ai
2S Rc 2S Rc ] (6)
-1
[
2 2
2S Rc
2
] [ 2 2
2S Rc
2
θro = cos S + Rc − Ao − cos-1 S + Rc − Am ] Figure 4 Generating triangle in pitch plane for face hobbing.
(7)
θri = θr − θro
The center distance from crown gear center to cutter (Radial),
(15)
Discrete points can be calculated along the cutter sweep on
S1 = √Am + Rc − 2Am Rc cos λ
2 2
the pitch plane. Using the general parametric formulas for a cir-
cle, the solution for a point, p, along the cutter path is: The second auxiliary angle,
(16)
(8)
px(t) = V + Rc cos t
(9)
[
η = cos Am cos ψm (Nc + Ns)
-1
S1Nc ]
py(t) = V + Rc sin t
The second roll angle,
(17)
Where t is a parameter that has the following ranges: φ2 = π − η1 − ψm − λ
For left-hand members,
(10) The auxiliary roll angle,
π + ψm − θri ≤ t ≤ π + ψm − θro (18)
Where, Therefore, all the local cutter positions can be wrapped and
(23)
Ai ≤ Aμ ≤ Ao transformed into the global Cartesian coordinate system.
The first roll angle as a function of cone distance, Calculating Local Cutter Coordinate System
(24) Since AGMA 929 effectively calculates the normal circular tooth
ρ thicknesses at a specific spot along the cutter path, the next step
φp1 (Aμ) = φp2 (Aμ) ρ2
1
is to determine the correct orientation of the normal plane. A
The local Cartesian coordinates describing a point, p, along complete coordinate system will be oriented to have the xnyn
the cutter sweep path can be calculated. plane normal to the cutter path, with the origin at each global
(25) coordinate of the cutter path calculated previously; this coordi-
px = −Rc sin (φp1 + φp2) + (ρ1 + ρ2) sin φp1
(26) nate system will be called Cn.
py = −Rc cos (φp1 + φp2) + (ρ1 + ρ2) cos φp1 Tangent axis, zn. The tangent axis is defined by a vector that is
tangent to the cutter sweep path at the location of the wrapped
Wrapping the Pitch Plane point. There are a couple of options for calculating this tan-
The pitch plane is a tangent plane to the lateral surface of a right gent vector. One could calculate the first derivative of the cut-
cone. This cone, as described earlier, is the “pitch cone.” The ter sweep path formulas for both face milling and face hobbing
diagrams take into account the set-up of the machine, the cut- so that the slope anywhere along that path can be predicted.
ter size, and the motion of the cutter. Rotation of the workpiece Once that is accomplished a vector can be constructed in three
also needs to be accounted for; this is accomplished by wrap- dimensions to describe this tangent axis.
ping the pitch plane around the pitch cone. This transforms the Since many assumptions are made throughout this method,
local coordinates calculated for the cutter sweep path into global the simplest method for calculating an approximate tangent
coordinates. These points in the global coordinates will define vector is using finite difference; the finite difference method is
the center of a tooth slot. The following formulas will transform described in the following equation:
(32)
a point, p, from the local XppYpp plane to the global Cartesian
coordinate system, CG.
The first step is to determine the location of the point in the
f '(x) = lim Δx—>0 ( f (x + Δx) − f (x)
Δx )
global Z axis direction, It is difficult to determine the correct size of Δx that will
(27) approximate the tangent close enough for this method. Since
zp = √px + py cos Γ
2 2
the cutter sweep path is smooth and continuous for all locations
Figure 5 Normal view to pitch plane displaying wrap point. Figure 6 Definition of a wrapped point’s local coordinate system.
along the arc, it is possible to quantify the error of our approxi- the gear terminology should be replaced with pinion terminol-
mation using a Taylor expansion. This is beyond the scope of ogy in the formulas if the pinion is the part being modeled.
this document, but is mentioned here for further exploration. Spiral angle for face milling with respect to cone distance,
(42)
The practical approach to applying the finite difference meth- 2 2
od is to calculate a neighboring wrapped point along the cutter -1 2Am Rc sin ψm − Am + Aμ
ψ (Aμ) = sin
2Aμ Rc
sweep path for each wrapped calculation point. A vector can
then be defined by passing through both sweep path points. As Generating angle for face hobbing with respect to cone dis-
the points become closer and closer the vector connecting the tance,
(43)
two points approaches a tangent line.
Normal axis — yn. The normal vector is defined as a vector Aμ2 + S12 − Rc2
-1
q (Aμ) = cos
2Aμ S1
radiating perpendicularly from the surface of the right cone at
the wrapped point location (Fig. 6). The normal vector at wrap Spiral angle for face hobbing with respect to cone distance,
(44)
point p = (px, py, pz) can be calculated by the following equation:
(33) -1 Aμ − Q cos q (Aμ)
ψ (Aμ) = tan
yn = (2px cos2 Γ, 2py cos2 Γ, − 2pz sin2 Γ) Q sin q (Aμ)
A brief derivation will follow to explain the calculation of the Slot width with respect to cone distance for the gear member,
normal axis. The general function of a right cone, (45)
C (x, y, z) = (x2 + y2) cos2 Γ − z2 sin2 Γ
(34)
[
WG (Aμ) = W'e 1 −
Aμ cos ψ (Aμ)
Am cos ψm ]
A
+ μ [tmP cos ψ (Aμ)] −
Am
A gradient, u, is always normal to a function when,
= uxex + uyey + uzez
(35) [ Aμ cos ψ (Aμ)
Am cos ψm G ]
A cos ψ (Aμ)
b (tan Φ1 − tan Φ2) + μ B
Ao cos ψ (Ao) cos
+
u (Am − Aμ)(tan Φ1 − tan Φ2) tan ξ
Where, Slot width with respect to cone distance for the pinion member,
< >
(36) (46)
ex = 1, 0. 0
A cos ψ (Aμ)
< >
(37) WP (Aμ) = μ p − Σbμ (tan Φ1 − tan Φ2) − WG (Aμ) +
ey = 1, 0. 0 Am cos ψm n
< >
(38)
ez = 1, 0. 0 Aμ cos ψ (Aμ) B
Ao cos ψ (Ao) cos
Let,
(39) Normal pitch radius with respect to cone distance (substitute
u (x, y, z) = x2 cos2 Γ + y2 cos2 Γ − z2 sin2 Γ proper pitch diameter and pitch angle for member to model),
Therefore, (47)
(40) D Aμ
u <
= 2x cos2 Γ, 2y cos2 Γ − 2z sin2 Γ > rN (Aμ) =
2 cos Γ cos2 ψm Ao
The cross vector will complete the definition of the local Normal base radius with respect to cone distance, concave,
(48)
Cartesian coordinate system. This vector is calculated by taking
rbN1(Aμ) = rN (Aμ) cos Φ1
the cross product of the tangent vector and the normal vector.
(41) Normal base radius with respect to cone distance, convex,
xn = zn yn (49)
rbN2(Aμ) = rN (Aμ) cos Φ2
After calculating all the vector directions for the local coordi-
nate system, all three of the vectors should be normalized. Figure 7 displays the base radii graphically. The value of rbN1
and rbN2 are shown as the same value because the pressure angle
Normal Circular Tooth Thickness Calculations of the concave and convex flanks are the same. This is typical
Generalizing AGMA 929-A06. As previously discussed AGMA but not mandatory for bevel gears without a hypoid offset.
929-A06 utilizes the technique of converting the spiral bevel Normal pinion circular tooth thickness at pitch line with
gear tooth to a virtual spur gear tooth to calculate the top lands respect to cone distance,
(50)
at the toe, mean, and heel; the equations presented in AGMA
929-A06 are unique for each of these points of interest. General A cos ψ (Aμ) B
tNP (Aμ) = bG (Aμ)(tan Φ1 − tan Φ2) − WG (Aμ) − μ
Ao cos ψ (Ao) cos
equations can be derived from AGMA 929-A06 so that the
tooth thickness can be calculated anywhere along the profile
and lengthwise direction. The generalized formulas for the con-
version to a virtual spur gear are presented here (Fig. 7). Also,
a handful of helpful formulas for calculating some spiral gear
tooth geometry with respect to cone distance, Aμ, are developed.
The cone distance range variable must correlate to the cone dis-
tances first chosen in the cutter path sections.
The notation in this section utilizes the terminology as if the
member being modeled is the gear. Unless specifically specified
Figure 7 Definition of the geometry of the virtual spur gear in the
normal plane.
Where
rwN (Aμ) = rN (Aμ) + w
[ t (A )
ywsP (Aμ) = rwN (Aμ) 1 − cos wNP μ
2rwN (Aμ) ]
w is a range variable for calculating the working radius. The Gear shift factor at working distance with respect to cone dis-
valid range for w is, tance,
(53) (61)
-b (Aμ) ≤ w ≤ a (Aμ)
coordinate system, Cn, is constructed that passes through a cal- are also broken into 10 segments. The part being modeled has a
culated tooth flank point, the transformation can be calculated: short addendum and long dedendum. Figure 10 shows the cal-
(67)
VG = CT VP
culated points and solid model sectioned through the normal
plane (Fig. 9) of a tooth. The calculated points go beyond the
Where, root line of the model because the calculated points were gen-
VG is a vector that passes through the exact same flank point, erated to the base circle radius of the virtual spur gear. When
but is described relative to the global coordinate system. doing the comparison the last four points of each profile will be
omitted, as they fall below the root fillet tangency point.
Results A linear measurement normal to the T900 surface to each cal-
Now that the method is complete, the next logical step is to culated point is used to measure the deviation between the cal-
determine just how accurate this model is when compared to culated points and the surface. The normal tooth thicknesses
the theoretical geometry of a spiral bevel gear. This method is at the toe and heel will be omitted in the comparison, as these
strictly approximate; many caveats have been discussed in the points are beyond the bounds of the T900 surface. The abso-
previous sections. To accomplish a comparison between this lute values of the measured deviations for each flank are given
method and a theoretical part there must be a reliable standard (Tables 2 and 3).
by which to compare it. Comparing the calculated points from The red values are measured deviations that exceed one one-
this method with a physical, cut part would introduce poten- thousandth of an inch. Overall, the results depict a model that
tial variations from the manufacturing process. For this reason very closely approximates the flanks predicted by the Gleason
alone it was decided to compare the calculated points to a differ- Works T900 software. The average deviation for each normal
ent, yet trusted, mathematical model. Gleason Works has devel- cross-section in the lengthwise direction is displayed in Table 4.
oped a commercial software package — T900 — that generates
an accurate point cloud describing the geometry of the tooth Conclusion
flank and fillet condition of bevel gears and pinions. While the The purpose of this presentation is to describe a procedure for
purpose of this software is far greater than just the generation calculating a very close approximation of the geometry of a spi-
of a point cloud, the other functionality is beyond the scope of ral bevel gear tooth.
this document. The point cloud produced from T900 can be This method is built upon the techniques and formulas
imported into a CAD package to assist in the generation of the described in AGMA 929-A06.
The second portion of the document compared the results of
Table 1 Basic spiral geometry of model this method with the results from a proven Gleason software
Pinion Gear package.
Number of teeth 11 57
Face width 2.0 2.0
The results are extremely close.
Diametral pitch 4.0 Many models have been generated since this original case
Shaft angle Non 90° study, and the subsequent models correlate well with the geom-
etry predicted by Gleason software.
exact geometry produced by the machine settings
for a particular design. The point cloud represents
the standard as to which this method is compared.
The member chosen for this analysis is a left-hand
spiral bevel gear generated by face mill completing.
Some of the basic geometry has been provided in
Table 1, but the exact details of the geometry are rela-
tively unimportant. An angular set was chosen as this
method has no shaft angle restrictions.
The output from T900 is a point cloud that
describes the flanks of a spiral bevel or hypoid gear
set. A gear tooth is divided into a 13- (profile) by-10
(length) grid for each flank; therefore 260 unique
points describe each tooth. For this analysis these Figure 9 Calculated points overlaid on a tooth modeled using Gleason T900.
discrete points were bridged together using Siemens
NX 8.5 CAD software. The point cloud was con-
nected in the profile and lengthwise direction with
curves generated by a cubic, polynomial regression.
Once the lattice of curves is generated the curves
are used as ribs to create a bi-cubic surface. Figure
9 shows the tooth surface generated from the T900
point cloud in gray.
The face width of the part is broken into 10 seg-
ments for the calculation; addendum and dedendum
Figure 10 Calculated points in the normal plane.