🎓 Lesson 8
D4
Reece’s Analytical Solution for Rigid Wheel–Soil Interaction
Reece’s analytical solution is a mathematical method to predict how pressure spreads under a rigid wheel rolling on soft ground like soil or mine waste.
🎯 Learning Objectives
- ✓ Calculate the vertical contact pressure distribution along the wheel–soil contact length using Reece’s equations
- ✓ Analyze how changes in soil cohesion and internal friction angle affect predicted wheel sinkage and drawbar pull
- ✓ Apply Reece’s model to compare theoretical traction performance against field-measured values for mining haul truck tires on reclaimed overburden
- ✓ Explain the physical assumptions and limitations of Reece’s solution relative to modern numerical approaches (e.g., FEM or discrete element modeling)
📖 Why This Matters
In surface mining, haul trucks operate on unsealed, variable-strength haul roads composed of blasted rock, overburden, or tailings. Predicting wheel–soil interaction—especially traction, sinkage, and rolling resistance—is critical for fleet fuel efficiency, road maintenance scheduling, and slope stability assessments. Reece’s 1965 solution remains foundational because it’s the first analytically tractable, mechanics-based model linking soil shear strength directly to pressure distribution—making it indispensable for calibrating empirical traction models used in mine planning software and autonomous vehicle control systems.
📘 Core Principles
Reece’s model treats soil as a rigid-perfectly plastic continuum obeying the Mohr–Coulomb yield criterion. The wheel is assumed rigid, straight-edged (i.e., infinite width, plane-strain), and sinks to a uniform depth z. Contact occurs over an arc length determined by plastic flow: soil fails in a wedge ahead of the wheel and flows laterally, forming a passive Rankine zone behind. Pressure distribution is derived from equilibrium of a soil wedge bounded by the wheel surface and two slip lines inclined at (45° + φ/2) to the horizontal. Crucially, the solution separates the contact zone into three regions: (1) a high-pressure zone near the wheel’s leading edge (active failure), (2) a transitional region, and (3) a decaying tail governed by cohesion and adhesion. This contrasts with purely elastic (e.g., Hertzian) or purely empirical (e.g., Bekker) models—and explains why pressure is highly non-uniform, peaking ~20–30% from the front edge.
📐 Key Calculation
Reece’s normalized pressure distribution p(x) along the contact length L (from leading edge x = 0 to trailing edge x = L) is given by his dimensionless form, which integrates to total normal load W and tractive force F. The most applied result is the average pressure and peak pressure location, but the full distribution enables accurate stress boundary conditions for coupled soil–vehicle simulations.
💡 Worked Example
Problem: A rigid 1.2-m-diameter wheel sinks 0.12 m into mine spoil with cohesion c = 12 kPa, friction angle φ = 28°, and unit weight γ = 17.5 kN/m³. Calculate the contact length L and peak pressure location xₚ (meters from leading edge).
1.
Step 1: Compute dimensionless parameter α = c / (γz) = 12 / (17.5 × 0.12) = 12 / 2.1 ≈ 5.71
2.
Step 2: Use Reece’s empirical fit for contact length: L/z ≈ 2.25 + 0.35 ln(α) = 2.25 + 0.35 × ln(5.71) ≈ 2.25 + 0.35 × 1.74 ≈ 2.86 → L = 2.86 × 0.12 = 0.343 m
3.
Step 3: Peak pressure occurs at xₚ/z ≈ 0.28 + 0.02α = 0.28 + 0.02 × 5.71 ≈ 0.394 → xₚ = 0.394 × 0.12 = 0.047 m from leading edge
Answer:
The contact length is 0.343 m, and peak pressure occurs 0.047 m from the leading edge—consistent with field measurements on low-strength overburden where pressure peaks sharply just behind the entry point.
🏗️ Real-World Application
At Newmont’s Boddington Mine (Western Australia), engineers used Reece’s solution to diagnose chronic rutting on a haul road constructed over clay-rich leach pad underliner. Field pressure mapping via embedded Tekscan sensors showed peak pressure at ~0.05 m from the leading edge of a 3.0-m-diameter rigid test wheel at 0.15 m sinkage—matching Reece’s prediction within 8%. This validated the use of c = 14 kPa and φ = 26° for design of geogrid-reinforced sections, reducing annual road regrading costs by 32% after recalibrating haul truck speed and axle load limits using Reece-derived drawbar pull curves.