Rolling Contact Fatigue & Pitting Calculator Back
Contact Fatigue Analysis

Rolling Contact Fatigue & Pitting Resistance Calculator

Compute Hertz contact pressure, subsurface shear stress, and L10 fatigue life in real time. Essential for bearing and gear pitting resistance design in CAE workflows.

Contact Type
Geometry
Radius R₁ (mm)20
Radius R₂ (mm)30
Material
Elastic Modulus E₁ (GPa)210
Elastic Modulus E₂ (GPa)210
Poisson's Ratio ν0.30
Loading
Normal Force F (kN)10.0

Color map: subsurface shear stress field | White dashed line: depth of τ_max

Contact Radius a (µm)
Peak Pressure p₀ (GPa)
Max Shear Stress τ_max (GPa)
Depth of τ_max (µm)
Elastic Approach δ (nm)
Relative L10 Index

Hertz Contact (Point)

$$a = \left(\frac{3FR^*}{4E^*}\right)^{1/3}, \quad p_0 = \frac{3F}{2\pi a^2}$$ $$\tau_\text{max}\approx 0.31\,p_0 \quad \text{at depth}\approx 0.48\,a$$

$\frac{1}{R^*}=\frac{1}{R_1}+\frac{1}{R_2}$,   $\frac{1}{E^*}=\frac{1-\nu_1^2}{E_1}+\frac{1-\nu_2^2}{E_2}$

What is Rolling Contact Fatigue?

🧑‍🎓
What exactly is "pitting" on gears and bearings? I see it mentioned in the tool's title.
🎓
Basically, pitting is a type of surface failure. When two curved surfaces, like a gear tooth and a pinion, roll under high load, tiny cracks can form just below the surface. These cracks grow and eventually cause small chunks of material to flake out, leaving pits. In the simulator, the "Normal Force F" slider directly controls the load that drives this process.
🧑‍🎓
Wait, really? The cracks start *below* the surface? So the contact looks fine until it suddenly fails?
🎓
Exactly! The maximum shear stress, which initiates the cracks, occurs at a subsurface point. For instance, in a ball bearing, the highest stress isn't right where the ball touches the raceway. Try adjusting the "Radius R₁" and "R₂" in the simulator. Making one surface flatter (larger radius) changes the stress depth and the contact patch size instantly.
🧑‍🎓
So how do engineers predict when pitting will happen? Is that what the "L10 fatigue life" is?
🎓
Yes! L10 life is the number of stress cycles where 90% of a group of identical components are expected to survive without pitting. It's calculated from the subsurface shear stress. A common case is designing a gearbox for a wind turbine—you need a 20-year life. In the tool, watch how increasing the "Elastic Modulus E" (making the material stiffer) reduces the contact area and increases pressure, which can drastically shorten the calculated L10 life.

Physical Model & Key Equations

The core of the analysis is Hertzian contact theory, which models the stress when two elastic, curved bodies are pressed together. The first step is to find the size of the contact area and the maximum pressure at its center.

$$a = \left(\frac{3FR^*}{4E^*}\right)^{1/3}, \quad p_0 = \frac{3F}{2\pi a^2}$$

Here, a is the radius of the circular contact area, F is the normal force, R* is the effective radius (from R₁ and R₂), and E* is the effective elastic modulus (from E₁, E₂, and ν). p₀ is the maximum Hertz contact pressure at the center.

Pitting fatigue is initiated by the oscillating subsurface shear stress. The maximum value of this shear stress and its depth below the surface are critical for life prediction.

$$\tau_\text{max}\approx 0.31\,p_0 \quad \text{at depth} \approx 0.48\,a$$

The factor 0.31 is derived from elasticity theory. This τ_max is the key driver for crack nucleation. The L10 fatigue life is then typically calculated using a power-law relationship (like the Lundberg-Palmgren model) where life is inversely proportional to a high power (e.g., cube or 9th power) of this shear stress.

Real-World Applications

Automotive Transmissions: Gear pitting is a primary failure mode. Engineers use this exact calculation to select gear steel hardness and heat treatment, ensuring the contact stress stays below the material's endurance limit for the vehicle's target lifespan.

Wind Turbine Main Bearings: These massive bearings support the rotor. They experience variable loads from wind gusts. Calculating Hertz pressure and L10 life is essential for maintenance scheduling and preventing catastrophic failures that require a costly crane-assisted replacement.

Railway Wheels and Rails: The contact between a steel wheel and the rail is a classic Hertzian contact problem. Managing contact stress is crucial to prevent rolling contact fatigue (RCF) defects like "squats" or "head checks," which can lead to rail breaks.

Precision Ball Screws and Linear Guides: In CNC machines and robotics, these components transmit motion with minimal friction. Their rated life is directly based on rolling contact fatigue calculations, balancing load, preload, and material to achieve the required precision and longevity.

Common Misconceptions and Points to Note

Let's go over some typical pitfalls you might encounter first with this type of calculation. The first is the misconception that a larger contact half-width 'a' is always safer. While the contact pressure p0 does decrease, the depth z_max at which the maximum shear stress τ_max occurs also becomes deeper. For instance, in a case-hardened gear, if the hardened layer depth is shallower than this z_max, the point of highest stress ends up in the softer core material, which can actually cause premature failure. If you use the tool to increase F and observe how z_max moves, you should get a clear picture of this relationship.

Next, watch out for errors in inputting material parameters. Poisson's ratio ν, in particular, is often casually entered as around 0.3, but it's part of the formula for calculating the equivalent Young's modulus E*: $1/E^* = (1-\nu_1^2)/E_1 + (1-\nu_2^2)/E_2$. Even changing ν from 0.25 to 0.33 can alter E* by several percent, which has an even larger impact on the L10 life. In practice, always verify the values from the material certificate.

Finally, don't forget that this calculation provides a baseline for an "idealized state". Actual bearings and gears are subject to countless factors not included in the calculation, such as lubricant effects, surface roughness, residual stresses, and assembly errors. For example, even if the tool calculates a life of 10,000 hours, poor lubrication can easily reduce it to less than a tenth of that. Treat these simulation results as a "benchmark for comparison" and aim for designs with a solid safety factor.

Related Engineering Fields

Hertzian contact and fatigue life calculations actually pop up in a much wider range of fields than you might think. The first that comes to mind is Tribology. Rolling contact surfaces always involve a lubricant film. This film pressure alters the contact pressure distribution, which can either reduce the maximum shear stress value or conversely promote surface pitting. The contact pressure you find with this tool becomes a crucial input condition for Elastohydrodynamic Lubrication (EHL) analysis.

Another is its close relationship with Material Surface Engineering. As briefly mentioned earlier, surface modification techniques like shot peening or carburizing specifically target the depth where maximum shear stress occurs (z_max≈0.48a). By introducing compressive residual stresses at the surface, they counteract the operational shear stress, dramatically extending fatigue life. This calculation is the first step in determining the optimal hardened layer depth.

Broadening the view further, the same equations are at work in the field of Biomechanics. Artificial joints, especially the ball-and-socket contact in hip joints, are essentially Hertzian contact. Here, evaluation goes beyond just the durability of metals or ceramics to assess whether contact pressure might cause bone resorption or dissolution. Fields that seem unrelated to mechanical design at first glance are actually connected by the same physical principles.

For Further Learning

Once you're comfortable with this tool's calculations, learning the "why" behind them will really open up new horizons. Start by grasping the concepts of stress tensors and principal shear stress. The τ_max displayed by the tool is the largest shear stress value among countless points within the contact body. This value is calculated from the three principal stresses (σ1, σ2, σ3) as $ \tau_{max} = ( \sigma_1 - \sigma_3 ) / 2 $. In a Hertzian contact field, it's theoretically derived that this occurs at a depth of 0.48a.

Next, I recommend looking into alternative approaches to fatigue life calculation. The method used in this tool is a relatively simple one based on maximum contact pressure p0 (the basis of Lundberg-Palmgren theory). For more detailed evaluation, you can use material-specific "stress-life (S-N) curves" and consider fluctuating loads with cumulative damage rules (like Miner's rule), or use methods based on "fracture mechanics" to calculate the number of cycles from an initial crack to failure. For example, Paris' law: $ da/dN = C(\Delta K)^m $.

Ultimately, you can move on to the step of making this calculation multidimensional. This tool models 2D (plane strain) cylindrical contact. Actual ball bearings involve 3D point contact. In that case, the contact area becomes elliptical, and the depth of maximum shear stress also changes (from about 0.48a to about 0.78a). Furthermore, handling "rolling contact" where the load position moves, as in gears, requires tracking a time-varying stress field. Fully understanding this simple model first is the best stepping stone to tackling those more complex problems.