The design and reliable operation of cylindrical gears in modern mechanical systems, especially within high-performance applications such as aero-engine transmissions, face severe challenges. These environments are characterized by extreme conditions including high speed, elevated temperature, and heavy load. Under such demanding service parameters, scuffing failure emerges as a predominant failure mode alongside traditional fatigue failures. Scuffing is a sudden, severe form of adhesive wear that results in the localized welding and tearing of mating tooth surfaces, leading to catastrophic damage and compromised system reliability. Unlike fatigue, scuffing is not cycle-dependent; it can occur almost instantaneously when local operating conditions exceed a critical threshold. Therefore, accurate prediction and prevention of scuffing are paramount in the design of high-reliability gear transmissions.
The core physical mechanism behind scuffing, particularly thermal scuffing, is the excessive rise in contact temperature at the gear mesh. This temperature surge, a combination of the bulk gear temperature and instantaneous flash temperature, causes the lubricant film to rupture. Subsequently, metal-to-metal contact occurs, leading to adhesion and material transfer. Consequently, the accurate calculation of the tooth surface temperature is the foundation of any scuffing load capacity evaluation method. Internationally, two principal theoretical approaches have been developed and standardized: the Flash Temperature Method and the Integral Temperature Method.
This article conducts a comprehensive comparative study of these methodologies as embodied in three prevalent standards: GB/Z 6413.1-2003 (Flash Temperature Method), GB/Z 6413.2-2003 (Integral Temperature Method), and the older HB/Z 84.4-1984 (also based on an Integral Temperature approach). The comparison spans theoretical foundations, detailed formulaic structures, critical influencing coefficients, and practical computational outcomes for a standard gear geometry. The aim is to elucidate the differences, similarities, and specific limitations of these standards to guide engineers in their application for the design of robust cylindrical gears.
1. Theoretical Foundations of Gear Scuffing Evaluation
The evaluation of scuffing risk for cylindrical gears is fundamentally based on temperature criteria. The two main methods differ primarily in how they characterize the critical temperature and how they calculate the weighted flash temperature contribution.
1.1 The Flash Temperature Method (Blok’s Theory)
This method is rooted in the pioneering work of Blok. It postulates that scuffing initiates when the maximum instantaneous contact temperature (flash temperature) at any point along the path of contact, added to the pre-existing bulk temperature, exceeds a critical material-lubricant dependent value. The flash temperature arises from the intense, localized frictional heat generated during the brief contact of asperities or surfaces under boundary or mixed lubrication conditions. The method involves calculating the flash temperature at discrete points along the line of action. The maximum of these values is then used for the safety check. This approach is sensitive to local pressure and sliding velocity peaks, making it theoretically suitable for identifying potential scuffing initiation sites. The governing inequality for the Flash Temperature Method is:
$$S_{B} = \frac{\Theta_{S}}{\Theta_{Bmax}} \geq S_{Bmin}$$
where $\Theta_{Bmax}$ is the maximum contact temperature, $\Theta_{S}$ is the permissible scuffing temperature, and $S_{Bmin}$ is the required minimum safety factor.
1.2 The Integral Temperature Method (Winter’s Approach)
Developed by Winter and colleagues, this method introduces the concept of an integral or weighted average temperature. It argues that scuffing is better correlated with a mean value of the flash temperature over the entire engagement cycle, rather than an instantaneous peak. The average flash temperature is weighted by a factor that accounts for its greater influence on the scuffing process relative to the bulk temperature. This weighted average is then added to the bulk temperature to obtain the integral temperature. A key characteristic of this method is its lower sensitivity to localized, transient temperature spikes, as it works with averaged values. The safety condition is expressed as:
$$S_{int} = \frac{\Theta_{intS}}{\Theta_{int}} \geq S_{intmin}$$
where $\Theta_{int}$ is the calculated integral temperature, $\Theta_{intS}$ is the permissible integral temperature, and $S_{intmin}$ is the required minimum safety factor.
The following table summarizes the core theoretical and applicability aspects of the three standards under review:
| Standard | Calculation Method | Theoretical Basis | Primary Application Scope |
|---|---|---|---|
| GB/Z 6413.1-2003 | Flash Temperature Method | Blok’s Instantaneous Flash Temperature Theory | Suitable for line/band-shaped Hertzian contacts with high Péclet numbers. Sensitive to local temperature peaks. |
| GB/Z 6413.2-2003 | Integral Temperature Method | Winter’s Mean (Integral) Temperature Theory | Less sensitive to local temperature peaks. Not recommended for cases where such peaks are dominant. |
| HB/Z 84.4-1984 | Integral Temperature Method (Simplified) | Simplified derivation from Integral Temperature Theory | Developed for aviation-grade, hardened steel, involute cylindrical gears (spur, helical, double-helical). Less comprehensive in coefficient consideration. |
2. Comparative Analysis of Calculation Formulas and Coefficients
A detailed examination of the formulas for bulk temperature, flash/average temperature rise, and the resultant scoring temperature reveals the structural similarities and key divergences between the standards.
2.1 Bulk Temperature ($\Theta_{M}$)
The bulk temperature represents the temperature of the tooth body just before it enters the contact zone. All three standards calculate it as the sum of the oil/sump temperature and a portion of the average frictional heating.
- GB/Z 6413.1: $$ \Theta_M = \Theta_{oil} + 0.47 \cdot X_{s} \cdot X_{mp} \cdot \frac{\int_{\Gamma_A}^{\Gamma_E} \Theta_{fl} \, d\Gamma_y}{\Gamma_E – \Gamma_A} $$ This formulation integrates the flash temperature along the path of contact and applies both a lubricant delivery factor ($X_s$) and a contact ratio factor ($X_{mp}$).
- GB/Z 6413.2 & HB/Z 84.4: Both use a simplified approach: $$ \Theta_M = \Theta_{oil} + C_1 \cdot \Theta_{flaint} $$ where $\Theta_{flaint}$ is the mean temperature rise, and $C_1$ is an empirical constant (typically 0.7) accounting for heat dissipation. HB/Z 84.4 adds the lubricant factor $X_s$ multiplicatively outside the sum: $ \Theta_M = (\Theta_{oil} + C_1 \cdot \Theta_{flaint}) \cdot X_s $.
2.2 Flash / Mean Temperature Rise
This term quantifies the temperature increase due to friction at the contact interface.
- GB/Z 6413.1 (Flash Temperature at a point): $$ \Theta_{fl} = \mu_m \cdot X_M \cdot X_J \cdot X_G \cdot (X_{\Gamma} \cdot \omega_{Bt})^{0.75} \cdot \frac{v_t^{0.5}}{a^{0.25}} $$ It calculates the instantaneous temperature at any specific point, depending on local load ($\omega_{Bt}$), geometry ($X_G, X_J$), and thermal material properties ($X_M$).
- GB/Z 6413.2 (Mean Temperature Rise): $$ \Theta_{flaint} = X_{\epsilon} \cdot \mu_m \cdot X_M \cdot X_{BE} \cdot X_{\alpha\beta} \cdot \frac{(K_{B\gamma} \cdot \omega_{Bt})^{0.75} \cdot v_t^{0.5}}{|a|^{0.25}} \cdot \frac{X_E}{X_Q \cdot X_{Ca}} $$ This is an empirical formula for the average value. It incorporates numerous coefficients for load distribution ($K_{B\gamma}$), running-in ($X_E$), tip relief ($X_{Ca}$), and contact geometry ($X_{BE}, X_{\alpha\beta}, X_{\epsilon}$).
- HB/Z 84.4 (Simplified Mean Temperature Rise): $$ \Theta_{flaint} = \mu_m \cdot X_M \cdot X_{BE} \cdot \frac{W_t^{0.75} \cdot \tilde{v}^{0.5}}{\tilde{a}^{0.25}} \cdot \frac{X_{\epsilon}}{X_Q \cdot X_{Ca}} $$ Compared to GB/Z 6413.2, this formula omits the running-in coefficient $X_E$, the pressure angle factor $X_{\alpha\beta}$, and uses a slightly different formulation for load and geometry parameters ($W_t$, $\tilde{v}$, $\tilde{a}$). Crucially, its default or calculation path for the average friction coefficient $\mu_m$ often differs.
2.3 Calculated Tooth Surface Temperature
This is the final calculated temperature against which the permissible temperature is checked.
- GB/Z 6413.1 (Maximum Contact Temperature): $$ \Theta_{Bmax} = \Theta_M + \Theta_{flamax} $$ Uses the maximum flash temperature from the path of contact scan.
- GB/Z 6413.2 & HB/Z 84.4 (Integral Temperature): $$ \Theta_{int} = \Theta_M + C_2 \cdot \Theta_{flaint} $$ Uses the weighted mean temperature rise. The weighting factor $C_2$ (typically 1.5) accounts for the greater influence of the flash component on scuffing risk.
2.4 Permissible (Scuffing) Temperature
The allowable temperature limit is typically derived from standardized gear tests (like the FZG test).
- GB/Z 6413.1: $$ \Theta_S = 80 + (0.85 + 1.4 X_W) \cdot X_L \cdot (S_{FZG})^2 $$ Based on FZG test load stage $S_{FZG}$, weld factor $X_W$, and lubricant factor $X_L$.
- GB/Z 6413.2: $$ \Theta_{intS} = \Theta_{MT} + X_{WrelT} \cdot C_2 \cdot \Theta_{flaintT} $$ Uses test gear data directly, modified by a relative weld factor $X_{WrelT}$ which compares the actual gear material to the test gear material.
- HB/Z 84.4: $$ \Theta_{Sint} = \Theta_{MT} + C_2 \cdot X_W \cdot \Theta_{flaintT} $$ Similar in form to GB/Z 6413.2 but uses the absolute weld factor $X_W$ instead of the relative factor.
2.5 Coefficient Comparison
The standards differ significantly in the number and treatment of correction factors, which greatly impacts the calculated result. The table below provides a summary comparison:
| Coefficient Category | GB/Z 6413.1 | GB/Z 6413.2 | HB/Z 84.4 | Remarks on Differences |
|---|---|---|---|---|
| Load Application | $K_A$, $K_V$ | $K_A$, $K_V$ | $K_A$ only | HB/Z 84.4 has a narrower range for $K_A$ and omits the dynamic factor $K_V$ in its specific formulation path. |
| Load Distribution | $K_{H\beta}$ ($K_{B\beta}$), $K_{H\alpha}$ ($K_{B\alpha}$) | $K_{H\beta}$, $K_{H\alpha}$, $K_{B\gamma}$ | $K_{H\beta}$, $K_{H\alpha}$, $K_{B\gamma}$ | GB/Z 6413.1 does not explicitly use a helix load factor $K_{B\gamma}$ in its flash formula structure. |
| Geometry & Contact | $X_{mp}$, $X_G$, $X_J$, $X_{\Gamma}$ | $X_{mp}$, $X_{\alpha\beta}$, $X_{BE}$, $X_{\epsilon}$ | $X_{BE}$, $X_{\epsilon}$ | HB/Z 84.4 omits the meshing factor $X_{mp}$ and pressure angle factor $X_{\alpha\beta}$. GB/Z 6413.1 uses point-specific geometric factors. |
| Surface & Running-in | – | $X_E$, $X_R$ | – | The running-in coefficient $X_E$ and roughness coefficient $X_R$ are only considered in GB/Z 6413.2, significantly affecting the calculated friction coefficient. |
| Material & Lubricant | $X_W$, $X_L$ | $X_{WrelT}$, $X_W$, $X_L$ | $X_W$ | A key difference: GB/Z 6413.2 uses a relative weld factor $X_{WrelT}$, while the others use an absolute factor. HB/Z 84.4 omits the lubricant factor $X_L$. |
| Thermal & Lubrication | $X_M$, $X_s$ | $X_M$, $X_s$ | $X_M$, $X_s$ | All include thermal flash coefficient $X_M$. Lubrication method factor $X_s$ definitions and values may differ slightly. |
The most consequential difference often lies in the determination of the average coefficient of friction $\mu_m$. The formula in HB/Z 84.4 tends to yield higher values (often 40-60% higher) than the more comprehensive calculation in GB/Z 6413.2, which accounts for roughness ($X_R$) and running-in ($X_E$). Since $\mu_m$ directly scales the temperature rise, this leads to systematically more conservative (higher) calculated temperatures in HB/Z 84.4.
3. Case Study: Computational Results for FZG A-Type Cylindrical Gears
To quantify the differences highlighted above, a series of calculations were performed on the standard FZG A-type test gear geometry. The basic parameters are: module $m_n = 4.5$ mm, pressure angle $\alpha = 20^\circ$, pinion teeth $z_1 = 16$, gear teeth $z_2 = 24$, face width $b = 20$ mm. The material is case-hardened steel. Calculations followed the procedures of the three standards, varying key operational parameters.
3.1 Variation with Input Torque (Load Stage)
Using a constant oil temperature of $60^\circ$C and speed, the torque was varied corresponding to FZG load stages 3 through 12. The calculated tooth surface temperatures are summarized below:
| Pinion Torque (N·m) | GB/Z 6413.1 $\Theta_{Bmax}$ (°C) | GB/Z 6413.2 $\Theta_{int}$ (°C) | HB/Z 84.4 $\Theta_{int}$ (°C) |
|---|---|---|---|
| 35.3 | 81.37 | 76.78 | 106.58 |
| 60.8 | 91.58 | 87.05 | 125.66 |
| 94.1 | 103.98 | 100.14 | 149.30 |
| 135.3 | 119.75 | 116.59 | 179.65 |
| 183.4 | 137.89 | 135.48 | 215.09 |
| 239.3 | 158.70 | 157.12 | 256.24 |
| 302.0 | 181.77 | 181.08 | 302.51 |
| 372.5 | 207.44 | 207.74 | 354.40 |
| 450.1 | 235.69 | 237.06 | 411.91 |
| 534.5 | 267.16 | 268.68 | 474.46 |
Observations: The results from GB/Z 6413.1 and GB/Z 6413.2 are remarkably close across the entire load range, with differences under 5%. This indicates that for this geometry, the flash temperature maximum and the weighted integral temperature yield similar values. In contrast, HB/Z 84.4 predicts consistently and significantly higher temperatures, with the divergence increasing with load. At the highest load, its result is approximately 77% higher than the GB/Z results. This conservative bias is primarily attributable to the higher average friction coefficient $\mu_m$ used in the HB standard.
3.2 Variation with Oil Sump Temperature
Holding the torque at the FZG 12th stage level, the oil temperature was varied from $40^\circ$C to $120^\circ$C.
| Oil Temp. $\Theta_{oil}$ (°C) | GB/Z 6413.1 $\Theta_{Bmax}$ (°C) | GB/Z 6413.2 $\Theta_{int}$ (°C) | HB/Z 84.4 $\Theta_{int}$ (°C) |
|---|---|---|---|
| 40 | 240.28 | 242.72 | 390.49 |
| 50 | 253.90 | 256.38 | 434.34 |
| 60 | 267.16 | 269.68 | 476.82 |
| 70 | 280.10 | 282.66 | 517.90 |
| 80 | 292.77 | 295.36 | 557.59 |
| 90 | 305.20 | 307.82 | 595.92 |
| 100 | 317.43 | 320.08 | 632.95 |
| 110 | 329.48 | 332.15 | 668.74 |
| 120 | 341.37 | 344.06 | 703.35 |
Observations: Again, the two GB/Z standards show excellent agreement (within ~1-2°C). The temperature increase is nearly linear with oil temperature. HB/Z 84.4 maintains its significantly higher absolute temperature prediction across the entire range. The sensitivity to oil temperature (slope) is also steeper for HB/Z 84.4, amplifying the difference at higher temperatures.
3.3 Variation with Gear Module
This is a critical analysis, as the gear geometry changes fundamentally. The torque intensity (load per unit face width and tooth) was kept constant while the module was varied from 2 mm to 7 mm.
| Module $m_n$ (mm) | GB/Z 6413.1 $\Theta_{Bmax}$ (°C) | GB/Z 6413.2 $\Theta_{int}$ (°C) | HB/Z 84.4 $\Theta_{int}$ (°C) |
|---|---|---|---|
| 2.0 | 557.66 | 685.09 | 1185.38 |
| 2.5 | 451.95 | 524.39 | 917.67 |
| 3.0 | 381.97 | 423.53 | 746.63 |
| 3.5 | 332.36 | 355.06 | 628.40 |
| 4.0 | 295.58 | 306.20 | 541.95 |
| 4.5 | 267.16 | 269.68 | 476.84 |
| 5.0 | 244.64 | 241.56 | 426.18 |
| 5.5 | 227.02 | 219.96 | 386.40 |
| 6.0 | 212.63 | 202.67 | 354.55 |
| 6.5 | 200.80 | 188.67 | 328.40 |
| 7.0 | 190.56 | 176.83 | 306.16 |
Observations: This reveals the most significant divergence between the two GB/Z standards. For modules below 4 mm, GB/Z 6413.2 calculates substantially higher temperatures than GB/Z 6413.1. At a 2 mm module, the difference is about 18.6%. This is because the influencing coefficients, particularly those related to lubrication ($X_L$) and roughness ($X_R$), which affect $\mu_m$, have a more pronounced effect at smaller scales. For modules between 4 mm and 7 mm, their results converge again and are very close. HB/Z 84.4 remains the most conservative across all modules, with its results diverging most severely at the smallest module.
3.4 Variation with Rotational Speed
Speed affects the sliding velocity and thermal conditions. Calculations were performed across a wide speed range from 1,000 rpm to 10,000 rpm at the FZG 12th stage load.
| Pinion Speed (rpm) | GB/Z 6413.1 $\Theta_{Bmax}$ (°C) | GB/Z 6413.2 $\Theta_{int}$ (°C) | HB/Z 84.4 $\Theta_{int}$ (°C) |
|---|---|---|---|
| 1000 | 243.85 | 246.08 | 440.24 |
| 2000 | 289.78 | 292.58 | 510.34 |
| 3000 | 323.39 | 326.59 | 557.10 |
| 4000 | 351.36 | 354.91 | 593.27 |
| 5000 | 376.06 | 379.90 | 623.18 |
| 6000 | 398.60 | 402.72 | 648.88 |
| 7000 | 419.62 | 423.99 | 671.55 |
| 8000 | 439.51 | 444.13 | 691.90 |
| 9000 | 458.53 | 463.38 | 710.42 |
| 10000 | 476.88 | 481.95 | 727.46 |
Observations: The trend shows a sub-linear increase in temperature with speed. The results from GB/Z 6413.1 and GB/Z 6413.2 are in very close agreement across the entire speed range. HB/Z 84.4 predicts higher absolute temperatures, but the rate of increase with speed is slightly lower than for the GB/Z standards, causing the percentage difference to decrease slightly at very high speeds (though the absolute difference remains large).
4. Conclusion and Engineering Implications
This detailed comparative study of scuffing load capacity standards for cylindrical gears leads to several important conclusions for design engineers:
- Conservative Nature of HB/Z 84.4: The HB/Z 84.4-1984 standard consistently yields calculated tooth surface temperatures that are 30% to over 100% higher than those from GB/Z 6413.2 for the FZG A-type gear under a wide range of conditions. This systematic conservatism is primarily driven by its simplified calculation of the average coefficient of friction $\mu_m$, which does not account for the mitigating effects of surface roughness and running-in as the more modern GB/Z 6413.2 standard does. Designs based solely on HB/Z 84.4 may be overly conservative, potentially leading to larger, heavier gearboxes than necessary.
- Agreement Between Flash and Integral Methods: For standard gear geometries (module ~4-7 mm) under varying load, speed, and oil temperature, the results from the Flash Temperature Method (GB/Z 6413.1) and the Integral Temperature Method (GB/Z 6413.2) are in very close agreement (typically within 5%). This provides cross-validation for these two theoretically distinct approaches for common cylindrical gear designs.
- Critical Divergence at Small Modules: A significant finding is the divergence between GB/Z 6413.1 and GB/Z 6413.2 for gears with small modules (below 4 mm). In this region, GB/Z 6413.2 predicts higher temperatures, making it the more conservative of the two modern standards. This highlights the importance of geometry-dependent coefficients. Engineers designing fine-pitch, high-density cylindrical gears must be cautious and understand which standard and which set of coefficients are most applicable to their specific case.
- Parameter Sensitivity: The calculations underscore the extreme sensitivity of scuffing risk to operational parameters. For instance, increasing torque from a low load stage to the 12th stage can increase calculated contact temperature by ~250%. High rotational speeds and elevated oil temperatures also contribute substantially to the temperature rise, emphasizing the need for effective cooling and lubrication systems in high-performance applications.
In practice, for the design of modern, high-performance cylindrical gears, the GB/Z 6413 series of standards (based on ISO/TR 13989) is recommended due to their more comprehensive and updated treatment of influencing factors like surface roughness, running-in, and relative material properties. Engineers should be aware of the historical context and inherent conservatism of the older HB/Z 84.4 standard. For non-standard geometries, especially very fine-pitch gears, a careful review of the applicable coefficients and, if possible, validation through testing is strongly advised to ensure both safety and optimal design efficiency.