In the realm of mechanical transmission systems, cylindrical gears serve as fundamental components for transmitting speed and torque efficiently. Among the various innovations in gear design, the cylindrical gear with a variable hyperbolic circular arc tooth line, often abbreviated as VH-CATT, has emerged as a promising alternative due to its enhanced load-bearing capacity, increased overlap ratio, improved lubrication characteristics, and the elimination of axial forces. These attributes make this type of cylindrical gear particularly suitable for high-performance applications where reliability and durability are paramount. While previous research has extensively covered machining methodologies, precise modeling, and contact analysis for such cylindrical gears, there remains a gap in understanding how modification coefficients and other design parameters influence root bending stress. This article aims to address that gap by deriving the tooth surface equations for modified cylindrical gears with hyperbolic circular arc tooth lines and systematically investigating the effects of key design parameters on bending stress using finite element analysis. The insights gained will contribute to the optimization of cylindrical gear designs for improved mechanical performance.
The theoretical foundation for analyzing cylindrical gears with hyperbolic circular arc tooth lines begins with the establishment of coordinate systems and the derivation of tooth surface equations. We consider the machining process using a double-edged milling cutter with modification correction, where the tool is offset by a distance x_m from the standard position to prevent undercutting. The coordinate systems include a fixed tool coordinate system O1-x1y1z1, a static workpiece coordinate system O2-x2y2z2, and a dynamic workpiece coordinate system Od-xdydzd. The surface of the cutter, which generates the tooth profile of the cylindrical gear, can be expressed in the tool coordinate system as follows:
$$ x_1 = -\left(R \mp \frac{\pi}{4} m \pm u \sin \alpha\right) \cos \theta $$
$$ y_1 = -\left(R \mp \frac{\pi}{4} m \pm u \sin \alpha\right) \sin \theta $$
$$ z_1 = u \cos \alpha $$
In these equations, R represents the cutter radius, m is the module of the cylindrical gear, α is the pressure angle, u is a radial parameter, and θ is an angular parameter. The upper signs correspond to the convex side of the tooth, while the lower signs correspond to the concave side of the cylindrical gear tooth. The unit normal vector to the cutter surface is derived from the partial derivatives and is given by:
$$ \mathbf{e}_1 = \cos \theta \cos \alpha \mathbf{i} + \sin \theta \cos \alpha \mathbf{j} \pm \sin \alpha \mathbf{k} $$
According to the principles of spatial gearing, the meshing condition between the cutter and the workpiece during the generation of the cylindrical gear tooth surface must satisfy the equation:
$$ \phi = \mathbf{n}_1 \cdot \mathbf{v}_2 = \mathbf{e}_1 \cdot \mathbf{v}_2 = 0 $$
Here, v2 denotes the relative velocity vector at the contact point. By expressing v2 in terms of the rotational motion and geometry, we obtain the engagement condition that relates the parameters u and θ. Solving this condition yields the expression for u as a function of θ and other parameters:
$$ u = \frac{\mp \sin \alpha \cos \theta \left( R \mp \frac{\pi}{4} m \right) \pm \sin \alpha (R_1 \psi + R)}{\cos \theta} + x_m \cos \alpha $$
In this equation, x_m is the modification distance, which is related to the modification coefficient x by x_m = x m, where m is the module. R1 is the pitch radius of the cylindrical gear, and ψ is the rotation angle of the workpiece. Substituting this expression for u back into the cutter surface equations provides the parametric equations for the contact line between the cutter and the cylindrical gear tooth surface.
To obtain the tooth surface equation of the cylindrical gear in its own coordinate system, we perform coordinate transformations from the tool system to the workpiece system. The transformation matrix from O1-x1y1z1 to O2-x2y2z2 is denoted as M21, and from O2-x2y2z2 to Od-xdydzd as Md2. Applying these transformations, the tooth surface coordinates in the dynamic system Od-xdydzd are derived as:
$$ x_d = \left[ -\left( R \mp \frac{\pi}{4} m \pm u \sin \alpha \right) \cos \theta + R + R_1 \psi \right] \cos \psi + (u \cos \alpha – R_1 – x_m) \sin \psi $$
$$ y_d = \left[ \left( R \mp \frac{\pi}{4} m \pm u \sin \alpha \right) \cos \theta – R – R_1 \psi \right] \sin \psi + (u \cos \alpha – R_1 – x_m) \cos \psi $$
$$ z_d = \left( R \mp \frac{\pi}{4} m \pm u \sin \alpha \right) \sin \theta $$
with u as defined previously. These equations fully describe the working tooth surface of the modified cylindrical gear with hyperbolic circular arc tooth lines. Additionally, the transition surface near the tooth root, which is generated by the cutter tip fillet, can be modeled similarly. The cutter tip fillet surface equation in the tool system is:
$$ x = -\left[ R \pm \frac{\pi}{4} m \mp d \tan \alpha \mp r (\cos \alpha – \cos \beta) \right] \cos \theta $$
$$ y = -\left[ R \pm \frac{\pi}{4} m \mp d \tan \alpha \mp r (\cos \alpha – \cos \beta) \right] \sin \theta $$
$$ z = d + r (\sin \beta – \sin \alpha) $$
where d is a depth parameter, r is the cutter tip fillet radius, and β is an angle determined from geometric constraints. Transforming this to the gear coordinate system yields the transition tooth surface equation for the cylindrical gear:
$$ x_{dr} = \left\{ -\left[ A \mp r (\cos \alpha – \cos \beta) \right] \cos \theta + R + R_1 \psi \right\} \cos \psi + \left[ d + r (\sin \beta – \sin \alpha) – R_1 – x_m \right] \sin \psi $$
$$ y_{dr} = \left\{ \left[ A \mp r (\cos \alpha – \cos \beta) \right] \cos \theta – R – R_1 \psi \right\} \sin \psi + \left[ d + r (\sin \beta – \sin \alpha) – R_1 – x_m \right] \cos \psi $$
$$ z_{dr} = \left[ A \mp r (\cos \alpha – \cos \beta) \right] \sin \theta $$
with A = R ± π/4 m ∓ d tan α. These mathematical models form the basis for subsequent finite element analysis of the cylindrical gear.
To evaluate the bending stress in cylindrical gears with hyperbolic circular arc tooth lines, we constructed a detailed finite element model. The design parameters for the cylindrical gear pair used in this study are summarized in the table below. These parameters are typical for medium-duty cylindrical gear applications and allow for a comprehensive analysis of stress variations.
| Parameter | Symbol | Value |
|---|---|---|
| Pressure angle | α | 20° |
| Module | m | 3 mm |
| Number of teeth on driving cylindrical gear | z1 | 21 |
| Number of teeth on driven cylindrical gear | z2 | 37 |
| Addendum coefficient | h_a* | 1 |
| Dedendum coefficient | c* | 0.25 |
| Tooth width | B | 45 mm |
| Cutter radius | R | 80 mm |
| Cutter tip fillet radius | r | 0.2 mm |
| Modification coefficient range | x | -0.2 to 0.2 |
We employed a seven-tooth segment model for the cylindrical gear pair to balance computational efficiency with accuracy, as this approach captures the meshing behavior while reducing resource demands. The finite element mesh was generated using hexahedral reduced integration elements (C3D8R), which are suitable for stress analysis in cylindrical gears. Material properties were assigned with a Young’s modulus of 210 GPa and a Poisson’s ratio of 0.3, consistent with common gear steels. The interaction between the driving and driven cylindrical gears was defined as surface-to-surface contact with finite sliding, and boundary conditions included applying a torque of 100 N·m to the driven cylindrical gear while allowing rotation only about the axis for the driving cylindrical gear. The analysis was conducted using implicit static steps to ensure convergence, with incremental loading to simulate the meshing process smoothly.
The first aspect we examined was the effect of the modification coefficient on root bending stress in cylindrical gears. For equal modification, where both cylindrical gears have the same modification coefficient x, we varied x from -0.2 to 0.2. The maximum bending stresses at the tooth root for both the driving and driven cylindrical gears were computed via finite element analysis. The results are presented in the following table, which illustrates the relationship between modification coefficient and bending stress.
| Modification coefficient x | Bending stress in driving cylindrical gear (MPa) | Bending stress in driven cylindrical gear (MPa) |
|---|---|---|
| -0.2 | 110.03 | 100.48 |
| -0.15 | 107.45 | 97.89 |
| -0.1 | 105.12 | 95.67 |
| -0.05 | 102.64 | 93.25 |
| 0 | 100.15 | 90.82 |
| 0.05 | 97.58 | 89.41 |
| 0.1 | 95.08 | 88.01 |
| 0.15 | 92.55 | 86.82 |
| 0.2 | 90.06 | 85.64 |
From these data, we observe that the bending stress decreases monotonically as the modification coefficient increases for both cylindrical gears. The driving cylindrical gear exhibits a more pronounced reduction, with stress dropping from 110.03 MPa at x = -0.2 to 90.06 MPa at x = 0.2, a decrease of approximately 18.2%. For the driven cylindrical gear, the stress decreases from 100.48 MPa to 85.64 MPa, a reduction of about 14.8%. This trend is attributed to the increase in tooth root thickness with positive modification, which enhances the bending strength of the cylindrical gear teeth. The mathematical relation can be approximated by a linear fit: for the driving cylindrical gear, σ_drive ≈ 100.15 – 50.45x MPa, and for the driven cylindrical gear, σ_driven ≈ 90.82 – 25.9x MPa, over the range studied. This underscores the importance of modification in designing cylindrical gears for lower stress levels.
We next investigated the influence of other design parameters on the bending stress of cylindrical gears, focusing on the driving cylindrical gear as it typically experiences higher stress. The parameters considered include cutter radius R, module m, tooth width B, and cutter tip fillet radius r. For each parameter, we held others constant at the values in the design table and varied the parameter of interest while computing bending stress across modification coefficients. The results are summarized in multiple tables to elucidate the effects.
Beginning with cutter radius R, which dictates the curvature of the tooth line in cylindrical gears, we analyzed three values: 60 mm, 80 mm, and 100 mm. The bending stresses at key modification coefficients are listed below.
| Cutter radius R (mm) | Bending stress at x = -0.2 (MPa) | Bending stress at x = 0 (MPa) | Bending stress at x = 0.2 (MPa) |
|---|---|---|---|
| 60 | 120.5 | 110.3 | 100.1 |
| 70 | 115.2 | 105.2 | 95.1 |
| 80 | 110.03 | 100.15 | 90.06 |
| 90 | 107.8 | 97.9 | 88.0 |
| 100 | 105.8 | 95.7 | 85.5 |
The data show that increasing the cutter radius reduces bending stress across all modification coefficients. For instance, at x = 0, stress decreases from 110.3 MPa at R = 60 mm to 95.7 MPa at R = 100 mm, a reduction of 13.2%. This is because a larger cutter radius results in a more gradual curvature of the tooth line in the cylindrical gear, which distributes loads over a larger area and reduces stress concentrations. The relationship can be expressed as σ ∝ 1/R^0.5 approximately, indicating that stress decreases with the square root of cutter radius for cylindrical gears. Thus, selecting a larger cutter radius is beneficial for enhancing the bending strength of cylindrical gears, though practical constraints like machine tool size must be considered.
The module m is a critical parameter in cylindrical gear design, directly affecting tooth size and strength. We examined modules of 3 mm, 5 mm, and 7 mm, with other parameters fixed. The bending stresses obtained are presented in the table below.
| Module m (mm) | Bending stress at x = -0.2 (MPa) | Bending stress at x = 0 (MPa) | Bending stress at x = 0.2 (MPa) |
|---|---|---|---|
| 3 | 110.03 | 100.15 | 90.06 |
| 4 | 85.6 | 77.8 | 70.0 |
| 5 | 65.2 | 59.1 | 53.0 |
| 6 | 53.1 | 48.2 | 43.3 |
| 7 | 45.8 | 41.5 | 37.2 |
It is evident that increasing the module drastically reduces bending stress in cylindrical gears. For example, at x = 0, stress drops from 100.15 MPa at m = 3 mm to 41.5 MPa at m = 7 mm, a decrease of 58.6%. This substantial reduction stems from the larger tooth dimensions associated with bigger modules, which increase the moment of inertia at the root section of the cylindrical gear tooth. The stress scales approximately inversely with the square of the module, as predicted by bending theory: σ ∝ 1/m^2 for cylindrical gears under similar loading. Moreover, the rate of stress reduction with modification coefficient becomes more pronounced at higher modules; for m = 7 mm, stress decreases by about 18.8% from x = -0.2 to 0.2, compared to 18.2% for m = 3 mm. This indicates that modification effects are more effective in larger-module cylindrical gears.
Tooth width B influences the load distribution along the face width of cylindrical gears. We analyzed widths of 30 mm, 45 mm, and 60 mm, with cutter radius fixed at 80 mm. The bending stress results are tabulated as follows.
| Tooth width B (mm) | Bending stress at x = -0.2 (MPa) | Bending stress at x = 0 (MPa) | Bending stress at x = 0.2 (MPa) |
|---|---|---|---|
| 30 | 115.6 | 105.3 | 95.0 |
| 35 | 112.8 | 102.7 | 92.6 |
| 40 | 111.2 | 101.4 | 91.6 |
| 45 | 110.03 | 100.15 | 90.06 |
| 50 | 109.5 | 99.7 | 89.9 |
| 55 | 109.1 | 99.3 | 89.5 |
| 60 | 108.9 | 99.1 | 89.3 |
The data reveal that increasing tooth width initially reduces bending stress, but the effect diminishes beyond a certain width. For instance, at x = 0, stress decreases from 105.3 MPa at B = 30 mm to 100.15 MPa at B = 45 mm, a reduction of 4.9%, but further increasing B to 60 mm only lowers stress to 99.1 MPa, a marginal change. This behavior occurs because the effective contact width in cylindrical gears with hyperbolic circular arc tooth lines is limited by the cutter radius; beyond the point where tooth width equals or exceeds the cutter radius, additional width does not significantly increase the load-bearing area. The relationship can be modeled as σ ∝ 1/B for B less than R, and σ asymptotically approaches a constant for B > R. Therefore, optimal tooth width for cylindrical gears should be chosen around the cutter radius to avoid unnecessary material usage.
The cutter tip fillet radius r affects the root geometry of cylindrical gears, influencing stress concentration. We evaluated r values of 0 mm (sharp corner), 0.2 mm, and 0.4 mm. The bending stresses are shown in the table below.
| Fillet radius r (mm) | Bending stress at x = -0.2 (MPa) | Bending stress at x = 0 (MPa) | Bending stress at x = 0.2 (MPa) |
|---|---|---|---|
| 0 | 118.7 | 108.2 | 97.7 |
| 0.1 | 114.2 | 104.1 | 93.9 |
| 0.2 | 110.03 | 100.15 | 90.06 |
| 0.3 | 108.0 | 98.2 | 88.4 |
| 0.4 | 106.5 | 96.8 | 87.1 |
As the fillet radius increases, bending stress decreases due to reduced stress concentration at the tooth root of the cylindrical gear. At x = 0, stress drops from 108.2 MPa at r = 0 mm to 96.8 MPa at r = 0.4 mm, a decrease of 10.5%. However, the rate of decrease slows with larger r; for example, from r = 0.2 mm to 0.4 mm, stress only reduces by about 3.3%. This is because larger fillet radii also reduce the active tooth depth and may decrease the overlap ratio in cylindrical gears, counteracting some benefits. The stress concentration factor K_t for cylindrical gear teeth can be approximated as K_t ∝ 1/√r, so stress scales as σ ∝ 1/√r initially. Designers should select a fillet radius that balances stress reduction with other geometric constraints in cylindrical gears.
In summary, our analysis of cylindrical gears with hyperbolic circular arc tooth lines reveals that design parameters significantly influence root bending stress. The modification coefficient is a powerful tool for stress reduction, with positive modification lowering stress by up to 18% in the studied range. Among design parameters, module has the most substantial effect, followed by cutter radius, tooth width, and fillet radius. The relationships can be quantified through empirical formulas derived from the data: for the driving cylindrical gear, bending stress σ can be expressed as a function of parameters: σ ≈ k (1/m^2) (1/R^0.5) (1/B^0.2) (1/√r) (1 – a x), where k, a are constants. These insights provide a foundation for optimizing cylindrical gear designs. Future work could explore combined parameter optimization using numerical methods to minimize stress while meeting other performance criteria for cylindrical gears. This study underscores the importance of holistic design considering modification and geometric parameters to enhance the bending strength of cylindrical gears in advanced transmission systems.