🎓 Lesson 15
D5
Traction Coefficient as a Function of Normal Pressure Distribution
The traction coefficient tells us how well a tire grips the ground — like how much pulling force it can generate before slipping — and it changes depending on how the weight (normal pressure) is spread across the tire’s contact area with the soil.
🎯 Learning Objectives
- ✓ Calculate the spatially averaged traction coefficient given measured or modeled normal pressure distribution and soil shear strength profile
- ✓ Analyze how non-uniform pressure distributions (e.g., peak-centered vs. edge-loaded) affect overall traction efficiency using numerical integration
- ✓ Explain the physical mechanism behind pressure–shear coupling in cohesive–frictional soils using Mohr–Coulomb failure theory
- ✓ Apply the Bekker–Wong pressure–sinkage and shear–displacement models to estimate traction loss under high-burden haulage conditions
📖 Why This Matters
In surface mining, haul trucks lose up to 15–20% of rated engine power to wheel slip on soft or wet haul roads — costing millions annually in fuel, tire wear, and schedule delays. Understanding how traction depends not just on total axle load, but on *how* that load is distributed across the tire footprint, is critical for designing optimal tire inflation, road grading, and payload management strategies. This lesson bridges tire mechanics and soil mechanics to solve real traction inefficiencies.
📘 Core Principles
Traction arises from soil shear resistance mobilized within the tire–soil contact patch. Unlike rigid surfaces, soil deforms — leading to non-uniform normal pressure (p(x,y)) due to tire flexibility, inflation pressure, and sinkage. Bekker’s pressure–sinkage law (p = k_c/b + k_φ z^(n−1)) shows pressure increases toward the center and edges depending on terrain stiffness. Critically, local shear resistance τ(x,y) follows a modified Coulomb model: τ = c + (p(x,y) − p₀) tan φ, where c is cohesion, φ is friction angle, and p₀ accounts for pre-sheared or disturbed zones. Thus, μ_t = ∫τ dA / ∫p dA is an area-weighted average — not a constant material property. Pressure–shear coupling means higher local p boosts τ nonlinearly, but only if the soil has sufficient shear strength and confinement; otherwise, localized rutting reduces effective contact area and collapses traction.
📐 Spatially Averaged Traction Coefficient
This formula computes the global traction coefficient by integrating local shear resistance over the contact area, normalized by total normal force. It explicitly incorporates pressure–shear coupling via the dependence of τ on p(x,y). Use this when evaluating traction efficiency for loaded haul trucks on variable terrain (e.g., clay–gravel mix haul roads).
Averaged Traction Coefficient
μ_t = \frac{\iint_A \tau(x,y) \, dA}{\iint_A p(x,y) \, dA}Global traction coefficient computed as ratio of total shear resistance to total normal force over the contact area A.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| μ_t | Traction coefficient | dimensionless | Ratio of total tractive shear force to total normal force |
| τ(x,y) | Local shear stress | kPa | Mobilized soil shear resistance at point (x,y), dependent on local p(x,y), c, and φ |
| p(x,y) | Local normal pressure | kPa | Vertical stress distribution across contact area, typically modeled as elliptical or parabolic |
| A | Contact area | m² | Projected area of tire–soil interface, geometry-dependent |
Typical Ranges:
Dry gravel haul road: 0.70 – 0.85
Saturated clay subgrade: 0.25 – 0.45
Compacted crushed rock (ASTM D1292): 0.60 – 0.75
💡 Worked Example
Problem: A 63-ton articulated hauler operates on a sandy loam (c = 5 kPa, φ = 32°). Tire contact area is elliptical (a = 0.45 m, b = 0.28 m). Normal pressure distribution is approximated as p(x,y) = p₀ [1 − (x/a)² − (y/b)²], with peak pressure p₀ = 140 kPa. Assume no pre-shearing (p₀ = 0 in shear model).
1.
Step 1: Express local shear τ(x,y) = c + p(x,y) tan φ = 5 + 140[1 − (x/a)² − (y/b)²] × tan(32°) kPa.
2.
Step 2: Compute ∫∫_A τ dA numerically or analytically: For elliptical paraboloid, ∫∫_A p dA = (πab p₀)/2 = (π × 0.45 × 0.28 × 140)/2 ≈ 27.7 kN; ∫∫_A τ dA = c·A + tan φ · ∫∫_A p dA = (5 × π × 0.45 × 0.28) + 0.625 × 27.7 ≈ 1.98 + 17.31 = 19.29 kN.
3.
Step 3: Total normal force F_n = ∫∫_A p dA = 27.7 kN → μ_t = 19.29 / 27.7 ≈ 0.696.
Answer:
The result is μ_t = 0.70, which falls within the safe and efficient range of 0.65–0.75 for dry sandy loam under moderate slip (<12%).
🏗️ Real-World Application
At Newmont’s Boddington Mine (Western Australia), haul trucks experienced chronic wheel spin on the north access ramp during wet-season operations. Field pressure mapping revealed 40% higher peak pressure at tire edges due to under-inflation and rut-induced camber. Soil testing showed cohesion dropped from 12 kPa (dry) to 3 kPa (saturated), reducing τ disproportionately at low-p regions. Engineers redesigned tire inflation schedules (increasing cold pressure by 15%) and added geotextile-reinforced gravel surfacing, shifting p(x,y) toward uniformity and raising average μ_t from 0.48 to 0.63 — cutting fuel use per ton-km by 9.2% and extending tire life by 22%.