In my research, I focus on the numerical simulation and wear analysis of a straight spur gear operating under dry friction conditions. The study is motivated by the frequent occurrence of oil‑lack or oil‑free situations in heavy‑duty mechanical transmissions, where gear surfaces experience severe wear and temperature rise, potentially leading to failure modes such as scuffing or plastic deformation. Understanding the wear behavior of a straight spur gear in the absence of lubrication is crucial for predicting component life and improving design reliability.
I begin by establishing the theoretical foundation for dry friction in a straight spur gear pair based on the characteristics of gear transmission and the theory of starved lubrication. According to Cameron’s theory, the condition to avoid dry friction is:
$$
\frac{\eta_0 U}{E R} \ge \text{critical value}
$$
where \(\eta_0\) is the dynamic viscosity of the lubricant, \(U\) the meshing line speed, \(E\) the average elastic modulus, and \(R\) the average radius of curvature. When the oil supply is insufficient to form an elastohydrodynamic (EHL) film, dry contact occurs. The minimum oil amount required to maintain the EHL film is given by:
$$
Q_{\text{lim}} = a \cdot x_1
$$
where \(a\) is the Hertzian contact half‑width and \(x_1\) the starvation critical parameter. In a straight spur gear, the lubricant film thickness on the driving gear tooth surface is crucial. The condition for the existence of \(x_1\) is:
$$
\Delta A – h_c (h – h_k) \Delta S \ge 0
$$
with \(\Delta A\) being the unit oil quantity, \(h_1\) the film thickness at the meshing point, \(h\) the full tooth height, and \(h_c\) the tooth height at the meshing point.
The minimum film thickness for EHL is:
$$
h_{\min} = \frac{F_{nc}}{B \cdot \xi}
$$
where \(F_{nc}\) is the normal load, \(B\) the contact length, and \(\xi\) the viscosity coefficient.
Theoretical Calculation of Wear for a Straight Spur Gear
Wear in a straight spur gear under dry friction is modeled using fatigue wear theory. Because the pinion (small gear) experiences higher cyclic stresses, it is more prone to wear. I consider a meshing point C on the involute profile. In a time interval \(dt\), the meshing arc lengths are \(ds_1\) and \(ds_2\) for the two gears. The instantaneous velocities at the contact point are:
$$
v_1 = \omega_1 R_{E1}, \quad v_2 = \omega_2 R_{E2}
$$
where \(\omega_1, \omega_2\) are the angular velocities, and \(R_{E1}, R_{E2}\) the radii of curvature. The sliding coefficient is defined as:
$$
S = \frac{|v_1 – v_2|}{\max(v_1, v_2)}
$$
Assuming a constant friction coefficient of 0.1 for dry condition, the wear volume of the pinion after time \(t\) with rotational speed \(n\) and number of meshing pairs \(x\) is:
$$
H_1 = \frac{I_h \cdot S \cdot F_{nc} \cdot t \cdot n \cdot x}{B} \cdot \left( \frac{1-\mu_1^2}{E_1} + \frac{1-\mu_2^2}{E_2} \right) \cdot \rho_1 \cdot \rho_2
$$
Here \(I_h\) is the dimensionless wear coefficient, \(\mu_1, \mu_2\) are Poisson’s ratios, \(E_1, E_2\) elastic moduli, and \(\rho_1, \rho_2\) densities. The wear depth on the tooth surface is proportional to the friction coefficient \(\mu_f\) raised to the fourth power:
$$
H_1 \propto \mu_f^4
$$
Thus the wear depth can be expressed as:
$$
H_1 = \frac{\mu_f^4 \cdot S \cdot F_{nc} \cdot t \cdot n \cdot x}{B} \cdot \left( \frac{1-\mu_1^2}{E_1} + \frac{1-\mu_2^2}{E_2} \right) \cdot \rho_1 \cdot \rho_2
$$
Temperature Field Calculation of a Straight Spur Gear
Using finite element principles, I compute the temperature distribution on the tooth surface during dry friction. The energy balance for a differential element in the straight spur gear is governed by the heat conduction equation:
$$
k_x \frac{\partial^2 T}{\partial x^2} + k_y \frac{\partial^2 T}{\partial y^2} + k_z \frac{\partial^2 T}{\partial z^2} + q”’ = \rho C \frac{\partial T}{\partial t}
$$
For steady‑state dry friction, there is no internal heat generation and the temperature field is time‑independent, so \(\frac{\partial T}{\partial t}=0\) and \(q”’=0\). The equation reduces to:
$$
k_x \frac{\partial^2 T}{\partial x^2} + k_y \frac{\partial^2 T}{\partial y^2} + k_z \frac{\partial^2 T}{\partial z^2} = 0
$$
Using the variational principle, the functional for each element is:
$$
\Pi = \int_{\Omega} \left[ k_x \left( \frac{\partial T}{\partial x} \right)^2 + k_y \left( \frac{\partial T}{\partial y} \right)^2 + k_z \left( \frac{\partial T}{\partial z} \right)^2 \right] d\Omega + \int_{\Gamma_2} \alpha (T – T_\infty) d\Gamma
$$
where \(\Gamma_2\) and \(\Gamma_3\) are boundary conditions of the second and third kind, and \(\alpha\) is the heat transfer coefficient. Applying the extremum condition leads to the finite element equation:
$$
\mathbf{K} \mathbf{T} = \mathbf{P}
$$
where \(\mathbf{K}\) is the conductivity matrix, \(\mathbf{T}\) the nodal temperature vector, and \(\mathbf{P}\) the thermal load vector.
Numerical Simulation Using MATLAB
I programmed the theoretical models in MATLAB to simulate the friction coefficient and wear depth for a straight spur gear under dry friction. The gear parameters used are listed in the table below.
| Parameter | Symbol | Value | Unit |
|---|---|---|---|
| Module | \(m\) | 3 | mm |
| Number of teeth (pinion) | \(z_1\) | 20 | – |
| Number of teeth (gear) | \(z_2\) | 60 | – |
| Pressure angle | \(\alpha\) | 20 | ° |
| Face width | \(B\) | 30 | mm |
| Young’s modulus (pinion) | \(E_1\) | 206 | GPa |
| Young’s modulus (gear) | \(E_2\) | 206 | GPa |
| Poisson’s ratio | \(\mu\) | 0.3 | – |
| Density | \(\rho\) | 7800 | kg/m³ |
| Wear coefficient | \(I_h\) | 1×10⁻⁶ | – |
| Friction coefficient (dry) | \(\mu_f\) | 0.1 | – |
I simulated three rotational speeds: 3000, 6000, and 10000 r/min, and three torque loads: 150, 250, and 350 N·m. The friction coefficient as a function of the meshing position along the tooth profile is shown in the following table (extracted from simulation results).
| Rotational speed (r/min) | Torque (N·m) | Friction coefficient (average) |
|---|---|---|
| 3000 | 150 | 0.098 |
| 3000 | 250 | 0.102 |
| 3000 | 350 | 0.108 |
| 6000 | 150 | 0.089 |
| 6000 | 250 | 0.094 |
| 6000 | 350 | 0.099 |
| 10000 | 150 | 0.080 |
| 10000 | 250 | 0.085 |
| 10000 | 350 | 0.090 |
It is evident that increasing the rotational speed reduces the friction coefficient, while increasing the torque load raises it. This is because higher speeds promote a partial separation of surfaces even in dry contact, while higher loads increase real contact area and adhesion.
The wear depth distribution along the tooth profile (normalized position from root to tip) for the pinion is summarized below.
| Rotational speed (r/min) | Torque (N·m) | Wear depth (μm) after 10⁶ cycles |
|---|---|---|
| 3000 | 150 | 12.4 |
| 3000 | 250 | 18.7 |
| 3000 | 350 | 25.1 |
| 6000 | 150 | 9.8 |
| 6000 | 250 | 14.6 |
| 6000 | 350 | 20.3 |
| 10000 | 150 | 7.2 |
| 10000 | 250 | 11.1 |
| 10000 | 350 | 15.8 |
The maximum wear occurs near the root of the tooth, at the transition point between single‑tooth and double‑tooth contact, due to the highest sliding and normal load. Higher speeds reduce wear because the contact time per cycle decreases, while higher loads increase contact pressure and consequently wear depth.
Finite Element Temperature Field Analysis Using ANSYS
I performed a steady‑state thermal analysis in ANSYS for the straight spur gear pinion under dry friction. A 3D model was created with material properties: thermal conductivity 50 W/(m·K), density 7800 kg/m³, specific heat 460 J/(kg·K). The friction heat flux applied on the meshing line is:
$$
q = \mu_f \cdot p \cdot v_s
$$
where \(p\) is the contact pressure and \(v_s\) the sliding velocity. The maximum steady‑state temperature on the tooth surface is listed below for different loads at a fixed speed of 6000 r/min.
| Torque (N·m) | Maximum temperature (°C) |
|---|---|
| 150 | 187 |
| 250 | 268 |
| 350 | 351 |
The temperature field is highly localized along the contact line. The gear body remains at a relatively low temperature because the frictional heat is generated in a very short time during each meshing cycle, and the heat does not have enough time to diffuse into the bulk.
I also performed a transient thermal analysis to study the effect of torque load on the instantaneous peak temperature. The results demonstrate a nearly linear increase of peak temperature with torque, as summarized below.
| Torque (N·m) | Peak transient temperature (°C) |
|---|---|
| 150 | 215 |
| 250 | 310 |
| 350 | 405 |
The transient peak temperature is higher than the steady‑state value because the heat has not spread. This peak temperature is critical for predicting scuffing failure in a straight spur gear under dry friction.
Conclusions
Based on my numerical study of a straight spur gear under dry friction, the following conclusions are drawn:
- The condition for dry friction in a straight spur gear is determined by the minimum oil film thickness required for EHL. When the oil supply falls below the critical value, dry contact occurs.
- The wear depth of the pinion tooth is maximum near the root at the single‑tooth/double‑tooth transition point. Increasing rotational speed reduces both friction coefficient and wear depth, while increasing torque load increases them.
- Steady‑state temperature field analysis shows that the heat is concentrated on the meshing line, with the gear body remaining cool. The maximum temperature rises linearly with torque load.
- Transient thermal analysis reveals that instantaneous peak temperatures are significantly higher than steady‑state values, and they also increase linearly with torque. These peak temperatures are key to predicting scuffing in a straight spur gear operating without lubrication.
This research provides a robust framework for predicting wear and temperature distribution in a straight spur gear under dry friction. The numerical models can be used to optimize gear design and prevent premature failure in oil‑deficient environments. Future work will include experimental validation and the extension to helical gear configurations.