Revised Lewis Bending Stress Capacity Model

Background: During operation, cylindrical gearset experiences tangential, radial, and axial (helical gears only) force components that induce bending, compressive, and shear stresses at the root area of the gear tooth. Accurate estimation of the effective bending stress at the gear root is a challenge. Lewis was the first person who attempted estimating the root bending stress of spur gears with some reasonable accuracy. Various gear standards and codes in use today are modifications and improvements of the Lewis model. Objective: This research aims at revising the Lewis model by making adjustments for dynamic loads, shear stresses, axial bending stress for helical gears, and stress concentration factor that is independent on the moment arm of tangential or axial force component. Methods: An analytical approach is used in formulating a modified formula for the root bending stress in cylindrical gears starting with the original Lewis model. Intermediate expressions are developed in the process and works from many previous authors are reviewed and summarized. The new model developed is used to estimate the root bending stress in four example gearsets of 0 to 41.41 helix angle and the results are compared with those of AGMA (American Gear Manufacturers Association) formula. Results: Analysis from the examples shows that neglecting the radial compressive stress over-estimated the root bending stress by 5.27% on average. When shear stresses are ignored, the root bending stress is under-estimated by 7.49% on average. It is important, therefore, to account for both compressive and shear stresses in cylindrical gear root bending stress. When the root bending stress estimates from the revised Lewis model were compared with AGMA results, deviations in the range of -4.86% to 26.61% were observed. The stress estimates from the revised Lewis formulae were mostly higher than those of AGMA. Conclusion: The new root bending stress model uses stress concentration factors (normal and shear) that are independent of the point of load application on the gear tooth. This decoupling of stress concentration factor from the load moment arm distinguishes the new model from AGMA formula and brings bending stress analysis in gear design in line with classical bending stress analysis of straight and curved beams. The model can be used for both normal contact ratio and high contact ratio cylindrical gears.


INTRODUCTION
A gear is a toothed disk used to transmit power and motion when mounted on a rotating shaft. Cylindrical gears have teeth-angle. A spur gear may be treated as a helical gear with a zerohelix angle [3]. The helix allows more gradual engagement of meshing teeth, which results in a smoother transfer of load from the driving teeth to the driven teeth. Therefore, helical gears are increasingly being used because of their relatively smooth and quiet operation, large load carrying capacity, and higher operating speed [4,5].
The resistance of gear teeth to failure in bending is called beam strength and gear failure in bending fatigue is a common problem [6 -8]. In 1892, Lewis modeled a gear tooth as a short cantilever beam on a rigid base with the transmitted load applied near the tip of the gear tooth [9, 10]. The maximum tensile stress occurs at the root radius on the loaded or active side of the gear tooth. Due to the repeated loading of a gear tooth, this region becomes the preferential site for the initiation of a fatigue crack. Fatigue failure is due to crack formation and propagation induced by repeated loading. A crack normally initiates at a discontinuity where there is a cyclic maximum stress [9, 10]. Cracks grow along planes normal to the maximum tensile stress [9] and when the growth becomes unstable, brittle fracture rapidly follows. Thru-hardened gears most often fail in bending fatigue due to a crack initiated at the surface in the root area. Because the surface hardness of casehardened gears is higher than the core value, the bending fatigue strength of the gear root surface can be higher than that of the core. Case-hardened gears generally fail in fatigue at the boundary of case-core hardness, except when there is a sharp stress raiser at the surface [11].
Case-hardening processes, such as induction hardening, flame hardening, carburizing, nitriding, and some other processes like shot-peening, hammering, and cold rolling, can increase the fatigue resistance of gears because they generate surface compressive residual stresses [12 -16]. The surface stresses tend to weaken the effects of stress raisers reducing fatigue damage. Finishing processes, such as grinding or skiving, tend to develop tensile residual stresses, thereby reducing the benefits of compressive residual stress in casehardening processes, but some compressive surface residual stresses usually persist.

MOTIVATIONS AND OBJECTIVE OF THE STUDY
The tangential load in spur gears induces both bending and direct shear stresses on the gear tooth, while the radial load produces only compressive stress. In helical gears, an axial load exists and it introduces both bending and shear stresses, in addition to those experienced by spur gears. The bending stress capacity models of American Gear Manufacturers Association (AGMA), International Standardization Organization (ISO), and Japanese Industrial Standards (JIS) account for the bending stress from the tangential force, compressive stress from the radial load and the effect of stress concentration, but not explicitly, for the shear stresses which are considered negligible. Accounting for shear stresses will improve the estimation of gear root bending stress.
The AGMA model is based on the load point being at the worst position, which corresponds to the highest point of single tooth contact (HPSTC) in a gear mesh for normal contact ratio gears. In normal contact ratio gear meshes, the contact ratio is less than 2. Gear meshes with a contact ratio of at least 2 are described as high contact ratio gears and multiple gear pairs are always in contact during operation. In normal contact ratio gears, one pair of gear teeth carries the transmitted load at HPSTC in the AGMA model, but near the tip of the gear, two pairs carry the transmitted load, indicating that load sharing occurs. The main issue in high contact ratio gears is the identification of the point of contact that results in the maximum bending stress for each tooth and whether that point occurs simultaneously for all the teeth in the mesh [17]. Hence the concept of "HPSTC" seems inapplicable when multiple pairs of teeth are in contact in a mesh. AGMA, therefore, has two bending stress capacity models for cylindrical gears, based on contact ratio. Also, the computation of the J-value in the AGMA formula is evaluated for HPSTC, which is complicated, especially for addendum corrected gearsets.
The AGMA bending stress concentration factor model is based on the work of Dolan and Broghamer, who experimented on photoelastic plastic gear materials. The stress concentration factor is related to the bending moment arm, which makes it to be dependent on the point of load application on the gear tooth. In classical stress concentration analysis of straight and curved beams, the location of the bending load on the beam does not affect the stress concentration factor. In these cases, the stress concentration factor is dependent on the local geometry and the notch sensitivity of the component material. This appears to be an obvious weakness of the AGMA bending capacity models. Another issue is that tests on metal gears suggest that the stress concentration factor obtained is not the same as that for photoelastic plastic gears [17]. The actual stress concentration value depends on the presence of deep tool marks and scratches (surface roughness), surface hardness, surface residual compressive stress, and load cycles [17,18]. Stress concentration arises from changes in the cross-sectional dimensions of components that are under stress [19].
As mentioned earlier, contact of two gear pairs occurs even in normal contact ratio meshes near the gear tips. Hence, in general, multiple pairs of gear teeth contacts exist in all gear meshes, but it is perhaps impossible to locate HPSTC for high contact ratio gear meshes. Because load sharing occurs near the gear tip in all gear meshes and the Lewis model assumes a load point near the gear tip, the model appears to hold promise in estimating gear bending stress for both normal and high contact ratio gear meshes. The Lewis model is also attractive due to the computational simplicity of the bending stress form factor. The objective of this study is to modify Lewis bending stress model such that it accommodates stress concentration factor that is independent of load application point and accounts for shear stresses in the root bending stress of cylindrical gears.

FORCES IN CYLINDRICAL GEARS
The power transmitted by a gearset generates torque and force loads that must be safely carried by the gears. Fig. (1) shows the forces generated during power transmission in helical gears. The torque and tangential force on the pinion tooth are: (1) Equation (1) has two entries and should be interpreted as Eq. (1a) and Eq. (1b) from left to right. The same rule should be applied to other equations of similar nature. The subscript 1 in Eq. (1) applies to the pinion in a gearset. Similar equations can be deduced for the gear by using Eq.2 as a subscript.
The radial and axial forces in Fig. (1) are given in Eq. (2).

(2)
The normal forces in the gear mesh

(3)
The influence of the helix angle in gear performance can be summarized by two parameters: the transverse pressure angle, and the base helix angle ѱ b . The relationships of these parameters with the helix angle are given in Eq. (4).

(4)
At low helix angles, ѱ ≤ 20 o [20], the base helix angle is approximately equal to the helix angle. At high values of helix angle, a significant difference arises between the base helix angle and the helix angle.

SERVICE LOAD FACTOR (K s )
Experience shows that the forces acting on equipment in service are generally higher than the rated or nominal values in gear drives. Practically, the design or service load is often estimated by multiplying the rated load with a service load factor, which is used to account for load increases during normal operations of gearsets. It is a load magnification factor and in gear design, it may be evaluated by a multiplicative model as was done in a previous study [21]:

(5)
Readers are referred to Osakue [21], AGMA documents [22], Osakue and Anetor [23] for the selection of K a and evaluation of K v K m and K r . However, a brief summary of the parameters is given in Appendix B.

LEWIS BENDING CAPACITY MODEL
Gear failure in bending fatigue is one of the common modes of failure. The bending capacity of gear teeth was first calculated to a reasonable degree of accuracy by Wilfred Lewis He considered a gear tooth as a cantilever beam on a rigid base with a maximum parabola inscribed inside the tooth shape. A parabolic shaped cross-sectional beam develops constant bending stress at the surface. For a gear loaded in bending, the maximum tensile stress occurs at the root radius on the loaded or active side of the gear tooth.
The assumptions made or implied in the derivation of the Lewis formula are [24, 25]: 1. The applied load is assumed to be static. In practice, it is dynamic and depends on several factors, therefore this assumption is not conservative. The internal dynamic overload factor component in the service load factor in Eq. (5) was introduced in gear design to remedy this problem.
2. The tangential or transmitted load is applied to the tip of a single tooth. This assumption is conservative because when gears mesh at the tip, more than one pair is usually in engagement. This means load sharing is ignored.
3. The transmitted load is distributed uniformly across the full face width. This is a non-conservative assumption and can be instrumental in gear failures involving wide teeth and misaligned or deflecting shafts. The mesh or mounting overload factor component of the service load factor in Eq. (5) is used to account for gear mesh misalignment, which leads to non-uniform load distribution over the gear face width. 4. The radial force component is negligible. This is a conservative assumption because it produces a compressive stress that subtracts from the tensile stress at the point of maximum stress. 5. Forces due to tooth sliding friction are negligible. This is reasonable since the sliding friction coefficient in properly lubricated gearsets is small. However, considering mesh friction improves the estimation of the bending stress at the gear root.
6. Shear stress from the transmitted load is ignored because it is considered negligible. This is a non-conservative assumption, though it seems reasonable. But considering them makes the bending stress estimate more realistic.
7. Stress concentration in the tooth fillet is negligible. Stress concentration effects were unknown in the time of Lewis but are now known to be important [19]. Therefore, stress concentration can no longer be ignored.
The Lewis formula may be expressed as: Eq. (6) is evaluated separately for the pinion and gear in gear design. The pinion is usually more vulnerable to bending stress failure, being of a smaller root tooth thickness. Eq. (7a) gives the expression for the modified Lewis bending stress form factor when the radial compressive stress is considered [26,27]. Thus if the radial compressive force is neglected as assumed by Lewis, the Eq. (7b) is obtained.
It is clear that Eq. (7a) will always give higher values compared to Eq. (7b) because of the compressive stress considered in Eq. (7a) resulting in lower bending stresses. This explains why the value of Y for modern gear standards that account for the compressive radial force is slightly higher than Y / . Note that Y or Y / values for the pinion and gear can be estimated from a single curve. This makes its use relatively easier than the AGMA-J values, which are evaluated at HPSTC and leads to the use of two separate curves, one for the pinion and the other for the gear. The AGMA-J value incorporates a stress concentration factor.

MODIFICATIONS OF LEWIS FORMULA
In general, when multiple teeth pairs are in contact in a mesh, then load shearing should be accounted for as in Eq. (8):

(8)
By incorporating the load sharing parameter in Eq. (8), the bending stress in both normal and high contact ratio gears can be estimated using a single expression.

Tangential and Axial Bending Stresses
For helical gears, Eq. (8a) neglects the bending stress from the axial force and also the shear stresses from the transverse and axial forces. To account for these stresses, a fictitious spur gear is considered that has an axial load component in addition to a tangential load component, as depicted in Fig. (2). The gear tooth is modeled as a cantilever beam loaded in two planes in the  1) The influence of the helix angle in actual helical gears will be considered later.

Revised Lewis Bending Stress Capacity Model
The Open Mechanical Engineering Journal, 2020, Volume 14 5 Fig. (2). Cantilever models of gear tooth in bending.

Combined Bending and Compressive Stresses
Considering the left diagram of Fig. 3, where a beam of a rectangular cross-section is loaded in two planes. The two perpendicular forces produce bending moments about two axes: the vertical force about the horizontal axis and the horizontal force about the vertical axis at the centroid of the beam.
The distribution of the bending stresses from the forces at the cross-section is shown to the right of the loaded beam. At point A, the two bending stresses have positive values and they add up to give the maximum tensile stress. At point D, the two bending stresses have negative values and they added to give the maximum compressive stress. At point B or C, one bending stress is positive while the other is negative, therefore, the resultant stress in magnitude will be less than that for A or D. Since tensile stresses are responsible for fatigue failure, point A is the critical point in this figure. Applied to a gear tooth, one of the corner points at the root will experience additive tensile stresses.
Eq. (13b) gives the estimate of the bending root stress from the bending influence of the tangential and axial forces in combination with the compressive radial force. But the shear stresses from both the tangential and axial forces are not accounted for. This is done in the next section.

Direct Shear Stresses
The transverse and axial forces also induce direct shear stresses on the equivalent spur gear. The direct shear stress from the transverse load is: The direct shear stress from the axial load is:

(15)
The tangential and axial shear stresses act on perpendicular planes, so the resultant shear stress is given in Eq. (16a). Eq.   Fig. (3). Bending stress distribution in rectangular section.

Equivalent Root Tensile Stress
The equivalent tensile stress at the gear root may be based on the distortion energy theory or maximum shear stress theory depending on whether the material is ductile or brittle [12, 28]. Most gears are made from ductile materials, so the equivalent tensile stress at the tooth root may be estimated by applying the distortion energy theory. For a plane stress situation, the equivalent tensile stress based on the distortion energy theory is:

REVISED LEWIS FORMULA
The dedendum circle is generally connected with the involute profile of a gear tooth with a fillet. This introduces a geometric discontinuity at the gear root resulting in stress concentration. The locations of size changes or discontinuity are sites of stress concentration where the maximum stress values are considerably higher than the nominal or average values. The ratio of the maximum stress to the nominal stress is generally called the stress concentration factor. Stress concentration was not known in the days of Lewis; but it is very important today because experimental and simulated results indicate that it can significantly increase local stresses [19]. Therefore, a stress concentration factor should be included in the Lewis formula for it to be more realistic. Since there are both normal and shear stresses, Eq. (17) may be modified as: where:

(22)
In Appendix A, a method for estimating λ a , using the rack tooth profile is described. It allows a single moment arm factor to be used for a gear profile standard. This is a new concept in gear design.

(23)
Eq. (19) is true for a fictitious spur gear with tangential and axial loads, which are considered as static forces. It needs to be corrected for dynamic load effects. Also, to apply it to an actual helical gear, it must account for the influence of the helix angle.

Accounting for Helix Angle of Helical Gears
The transfer of power between gears in a mesh takes place by means of contact between the active teeth. In helical gears, contact of active teeth occurs in the normal plane. Therefore, the design analysis of helical gears should be preferably done on the normal plane. The normal plane of a helical gear intercepts the pitch cylinder so that the gear tooth profile generated in it has the same properties as the actual helical gear [29]. The normal plane, therefore, may be used to define an equivalent spur gear for a helical gear. The diameter of the equivalent spur gear depends on the base helix angle because it gives an accurate estimate of the radius of curvature of the equivalent spur gear on the normal plane of contact [28]. In helical gears, the normal bending force, acts perpendicular to the tooth similar to the transverse force F t in spur gears. The analogous bending stress from this force may be expressed using normal plane parameters as was done in a previous work

Corrections for Dynamic Load, Mesh Misalignment, and Effective Contact Length
Eq. (26) makes a correction for the helix angle, but the loading is still considered static and ignores friction load effects. Introducing the service load factor into the Lewis model will account for dynamic load and mesh misalignment effects.
Surface roughness seems to have a great influence on the actual contact area when two bodies are in relative sliding motion [31]. Factors, such as thermal gradient, centrifugal forces, work hardening, residual stresses [12], etc., can distort pinion or gear shape and lead to teeth mismatch so that full contact does not occur over the nominal face width of meshing gears. In general, the effective face width factor is assumed to account for surface roughness, surface treatment quality, and miscellaneous effects that make contact over the full nominal face width of a gear impossible. AGMA [22] suggests a value of 0.95 for the effective face width factor of helical gears; this value is adopted here for spur gears also.
Taking account of the service load factor,K s and the effective facewidth factor λ e ; the revised Lewis capacity model for the gear root bending stress of cylindrical gears is: The parameters T, b, d and Y / n in equations (25b) to (30) apply to either the pinion or gear in a mesh. Separate equations for the pinion or gear may be obtained by appending subscript 1 or 2, respectively, to these parameters. It is important to note that the Lewis form factor Y / n is applicable only to standard gears; that is, gears with no profile shift or modification.

Stress Concentration Factors
Stress concentration gives rise to unusually high stresses near a stress raiser, while much lower values exist on the remainder of the cross-section. Practically, some factors tend to limit stress concentration effects, and these include local plastic deformation, residual stresses, notch radius, component size, and load type. To account for these influences, the ideas of notch sensitivity and effective stress concentration were introduced [19]. Notch sensitivity is a material's response to stress concentration, and it is assessed by the notch sensitivity factor that has values between 0 for complete insensitivity and 1 for perfect sensitivity.
Due to the notch sensitivity of materials, the effective stress concentration factor is less than the theoretical value. Eq. (31) defines the relationship between these parameters.

(31)
According to Zahavi [32], the theoretical bending stress concentration factor in (steel) gears is in the range from 1.5 to 1.7. For steel gear materials, it is reasonable to expect normal notch sensitivity values in the range of 0.80 for normalized materials to 0.95 for quenched and tempered materials. Applying Eq. (31) to Zahavi's findings [32], the effective bending or normal stress concentrator factor would be in the range of 1.4 to 1.67 for steel gears. A bending stress concentration factor of 1.4 is used in a Japanese standard [33] for steel gear materials. MITCal [34] used a bending stress concentration factor of 1.5 in some gear design examples. According to Dobrovolsky et al. [11], the normal stress concentration factor is approximately 1.1 to 1.7 for steel gears and 1.2 for cast iron gears. Therefore, the estimates of the effective normal stress concentration factor in the range from 1.4 to 1.67 seem reasonable for steel gears.
The effective shear stress concentration is taken as 1.7 to 2.2 for flexible spline or harmonic gear teeth generated with pinion cutter and 1.6 to 2.0 when the teeth are hobbed [35]. Harmonic gear teeth are produced with hobbing and shaping methods, similar to those of cylindrical gears. Based on the above information and using engineering judgment, Table 1 is suggested for use in gear bending stress estimation during initial sizing. The values of the stress concentration factors in the table have no relationship with the point of load application on a gear tooth. The entries for case-hardened steel gears are for surface root bending stress and are informed by the beneficial effects of compressive surface residual stresses they generate. For subsurface root stress, values for quenched and tempered steels may be used. When values of stress concentration factors in Table 3 are used in Eq. (27), Eq. (29), and Eq. (30), the stress concentration factor is independent of the point of load application on the gear tooth. This eliminates the weaknesses highlighted about the AGMA model, whose stress concentration factor depends on the point of load application on the gear tooth.

DESIGN EXAMPLES
The new design formulas presented in the previous sections are applied in four design examples, taken from the indicated references. The problem statements in the design cases have been paraphrased and the design parameters have been converted to metric units where necessary by the authors. The equations presented were coded in Microsoft Excel for computational efficiency. The goal is to estimate the root bending stresses using the new formulae and make comparisons with those from AGMA. AGMA standards are perhaps the most popular gear standards in use and have a good reputation amongst gear designers and manufacturers.

Design Problems
The four of the design problems are considered below. The solutions to the problems are available from the stated references and comparisons will be made with the estimates from the formulae presented in the sections above. The problems cover a wide range of helix angles, which span from 0 to 41.41. Note that the 0 helical gear is actually a spur gear.

Example 1
A gearset of steel and ductile cast iron transmits a torque of 1694.8 nm at the pinion at 406 rpm. The gearset basic size data are: z 1 = z 2 = 20, d 1 = d 2 = 127 mm, b = 25.4 mm. The gearset has a normal pressure angle of 20° and a helix angle of 0°. What is the expected root bending stress? This example was presented in a previous study [21].

Example 2
A helical steel gearset for a milling machine drive is to transmit 48.5 kW from an electric motor with a pinion speed of 3450 rpm and a gear speed of 1100 rpm. The gearset has a normal pressure angle of 20 o and a helix angle of 15 o . The pinion has 24 teeth, gear has 75 teeth, the normal module is 2.17 mm and the face width is 57.15 mm. Determine the root bending stress on the pinion [6, pp. 461 -462].

Example 3
A 17-tooth helical steel pinion with a right-hand helix angle of 30 o rotates at 1800 rpm when transmitting 3 kW to a 52-tooth helical steel gear. The gearset has a normal pressure angle of 20 o , normal module of 2.54 mm, and a face width of 38.1 mm. Determine the root bending stress on the pinion [12 pp. 771 -773].

Example 4
A 15-tooth helical steel pinion with a right-hand helix angle of 41.41 o rotates at 2500 rpm when transmitting 3.75 kW to a 24-tooth helical steel gear. The gearset has a normal pressure angle of 20 o , normal module of 2.54 mm, and a face width of 29.21 mm. Determine the root bending stress on the pinion [3, pp. 658 -662]. Table 2 summarizes the basic gearset dimensions and load data, for example, 1 to 4. The service load factors were evaluated and used for root bending stress estimation. Table 3 gives the results for the problems for the revised Lewis model. The entries in the "A" column are obtained from Eq. (27) and Eq. (30). The entries in the "B" column are obtained from Eq. 29a and Eq. (30). The entries in the "C" column are obtained from Eq. 29b) and Eq. (30). Table 4 shows the AGMA bending stress solutions being compared with the revised Lewis model. The stress concentration values used in the bending stress estimations were taken from Table 3.

DISCUSSION
The examples cover a wide range of helix angles from 0 o to 41.41 o . This range essentially covers values in the common practice of cylindrical gear design and manufacture. The influence of the helix angle is directly indicated in equation (30) through the base helix angle. However, the stress correction factor k t also indicates a secondary influence of the helix as revealed in (Eq. (27). The solution of Eq. (30) is unique when compared with current gear design standards. For instance, a different chart is not required for the helix angle factor, as is done in the AGMA approach. Table 3 gives the results for the new revised Lewis model. Column"A" entries are the new accepted solutions that take account of both the radial load and shear stresses. Column "B" entries ignore the radial load but account for shear stresses. Column "C" entries account for the radial load but ignore shear stresses. The "B-DEV" column shows the percentage deviation of the "B" column from the "A" column. This column shows that the root bending stress is over-estimated by 4.66 to 5.73%, with an average of 5.27%. The "C-DEV" column shows the percentage deviation of the "C" column from the "A" column. This column shows that the root bending stress is underestimated by 6.67 to 8.16%, with an average of 7.49%. Clearly the influence of shear stresses is important if a more accurate model is desired in root bending stress estimation. Table 4 shows the bending stresses for comparison purposes. The AGMA values are given in column 3 and column 4 gives values for the revised Lewis formula. Column 5 shows the percentage deviations of the new solutions form the AGMA values in the range of -4.83% to 26.61%. It should be noted that getting accurate values of AGMA-J factor for helical gears is somewhat challenging. The new model solutions are thus, somewhat conservative compared with AGMA results but indicate good comparison.

CONCLUSION
A new revised Lewis bending stress capacity model is formulated and provides a single expression (Eq. (30)) for both spur and helical gears and normal and high contact ratio cylindrical gears. It accounts for both radial compressive stress and shear stresses in cylindrical gearsets. Also, the helix angle is expressly taken into account in Eq. (30) and Eq. (27). Note that a spur gear is a special case of a helical gear when the helix angle is 0 o . Eq. (27) and Eq. (30) eliminate the use of multiple charts or tables in estimating root bending stress of thecylindrical gears. The equations use stress concentration factors that are independent of the point of load application on the gear tooth. This brings bending stress analysis in gear design in agreement with classical bending stress analysis of straight and curved beams. Also, the need to identify HPSTC is eliminated and the estimation of Lewis form factor is easier. Load shearing is uniformly applied to both spur and helical gears through the parameter . The use of the rack profile for a gear profile standard in the determination of the moment arm is a new concept introduced in this study. It allows a single moment arm factor to be used for all modules in a gear profile standard. These features make the new model unique from those of AGMA and ISO.
Four design examples of bending stress evaluations with helix angles from 0 to 41.41 were considered. The torque range is 14 to 1695 nm and the estimated root bending stress range is from 27 to 882 MPa. Though higher torques may occur in practice, the stress of 882 MPa is too high and likely unacceptable, since most steel gears would not be able to sustain such a stress value for long. The study results from the new model show that ignoring the radial load led to an overestimation of 5.27% in the root bending on the average. When shear stresses are ignored, an under-estimation of 7.49% in the root bending on the average was obtained. Accounting for shear stresses is, thus, very important. The root bending stress estimates from the new model for four gearsets were compared with the AGMA results. The differences between the two estimates are in the range of -4.83 to 26.61%. Because gear design is complicated, the design procedures are usually not precise [36]. Hopefully, this study should help improve the accuracy of root bending stress estimation for cylindrical gears.

ETHICS APPROVAL AND CONSENT TO PARTI-CIPATE
Not applicable.

HUMAN AND ANIMAL RIGHTS
Not applicable.

CONSENT FOR PUBLICATION
Not applicable.

AVAILABILITY OF DATA AND MATERIALS
The authors confirm that the data supporting the findings of this study are available within the article.

FUNDING
No grant money from private or public sources or any other funds were used in the research work presented.

CONFLICT OF INTEREST
The author declares no conflict of interest, financial or otherwise.

ACKNOWLEDGEMENTS
This study was supported with research facilities from the College of Science, Engineering and Technology (COSET) and the University Faculty Development unit of Texas Southern University.

APPENDIX A: BENDING MOMENT ARM FACTOR
A diagram of the bending stress form factor curves [33] for different addendum correction factors converge to one point for a rack tooth. This indicates that the bending stress form factor for a rack is independent of the amount of addendum correction applied to the gear tooth. Now the rack tooth profile envelops all possible involute gear tooth shapes for a tooth standard and hence has a fixed tooth width at the root where the maximum tensile stress occurs. These reasons suggest that the bending moment arm factor λ a may be estimated reliably by considering the basic rack profile of a gear tooth system. Fig. (A1) shows the basic rack profile for 20 and 25 involute gear tooth standards. The fillet radius factor of 0.35 is a popular value in AGMA recommendations [13]. The root thickness is defined at the intersection of the fillet radius and the straight flank line of the rack tooth. A direct measurement of the root thickness which is equal to K for a module of 1 mm (Eq. 22a), can be made from these diagrams. Note that the fillet is tangential to both the straight flank line and the root or deddendum circle horizontal line shown as hidden line. Eq. (A1) is obtained from Eq. (22b).

(A1)
The value of Y / for a rack tooth is required in Eq. (A1). In Table A1, column 2 shows Y / values for some popular involute gear tooth standards. The values of K in the same table are obtained from Fig. (A1). The forth column in the table gives the values of λ a for 20° and 25° involute gear tooth standards based on Eq. (A1).

APPENDIX B: SERVICE LOAD FACTOR ESTIMATION
The service load factor accounts for the fact the forces acting on equipment in service are generally higher than the rated or nominal values in gear drives. It is a load magnification factor and may be estimated [11,12,21,22,23,33,37] as:

(A2) Application Overload Factor, K a
The application overload factor accounts for load increases caused by the power source device and the driven or load device on gear drives. It is indicative of the influences of the accelerations and decelerations of external masses connected to the gear drive. Table A2 is an example of values for this factor for 8 to10 hrs of daily operation of a gear drive.

Internal Overload Factor, K v
The internal overload factor, K V accounts for load excitations caused by non-conjugate action, backlash, profile error, pitch error, dynamic imbalance, etc. of meshing gears The internal overlad factor may be estimated using Eq. (A3a) for spur gears and Eq. (A3b) for helical gears.

(A4)
Parameters a1 and a2 depend on AGMA gear quality numbers that range from 0 to 12, and are similar to ISO quality numbers. Gear hobbing can produce 7 to 10 gear quality numbers, shaved gears may have 6 to 8 quality numbers, and ground gear can have 2 to 7 quality numbers. Careful production practice can improve on quality numbers. Lower numbers represent higher gear quality. For AGMA/ISO gear quality numbers in the range of 6 ≤ q n ≤ 12:

(A5)
An estimate of K v above 1.5 should be considered as probably unacceptable [23]. For uni-directional loading and as a guide for gear tooth profile quality selection, commercial quality gears may have 1.25, < K v <1.5 premium quality gears may have 1.15 ≤ K v <1.5, and precision quality gears may have K v <1.15. For high speed applications, especially those above 20 m/s, methods that account for gear material properties, mass and inertia of the gears, and actual tooth profile errors should be used to estimate K v .

Mesh Overload Factor, K m
The mesh or mounting overload factor K m takes care of non-uniform load distribution along the tooth facewidth. The AGMA method presented here is recommended for normal, relatively stiff gear assemblies. It is limited to straddle mounted gears, gear aspect ratio of λ b ≤ 2.0, gear face width, b ≤ 1000 mm (40"), and assumes full mesh contact over the face of narrowest gear.
The mesh overload factor may be estimated using Eq. (A6a) for spur gears and Eq. (A6b) for helical gears.

(A6)
Based on AGMA recommendation [12, 22]: (A7) C mn accounts for misalignment of gear pitch cylinders due to non-elastic deformations, principally influenced by the quality of gear mounting. It is evaluated as: Table A3 gives the values of the coefficients A 1 , A 2 , and A 3 . C me accounts for misalignment of gear pitch cylinders due to elastic deformation of shaft, housing, etc. It is a function of the gear aspect ratio and gear face width. Eq. (9) and Eq. (A10) are for λ b and C me , respectively.

(A9)
Based on the range of values for the gear facewidth, then:  Eq. (A10) assumes λ b = 0.5 when λ b < 0.5. C mp makes correction for the location of a gearset on a shaft relative to the mid-span. It's value is based on λ P and Fig. (A2) indicates how it is evaluated. That is:

(A11)
Other coefficients in Eq. (A7) are given in Table A4 .  From [ 37 ], the frictional load factor for enclosed gear drives may be approximated as:

(A17)
Eq. (A17) is useful during design verification or validation when the gearset is already sized. During initial sizing of gears, the gear dimensions are unknown, hence V t , V s and ℓ m cannot be evaluated. Therefore, a value of 1.1 is suggested for K in enclosed gear drives and 1. 15 is suggested for open gear drives during initial sizing of cylindrical gears.