🎓 Lesson 1
D1
Getting Started with Tire–Soil Contact Pressure Distribution Modeling
It's how the weight and force of a mining vehicle's tire spreads out onto the ground beneath it — like how your foot presses into sand, but measured precisely for safety and efficiency.
🎯 Learning Objectives
- ✓ Calculate maximum contact pressure using tire load, inflation pressure, and footprint geometry
- ✓ Analyze how soil modulus and tire deflection influence pressure distribution shape (e.g., elliptical vs. uniform)
- ✓ Apply empirical and semi-empirical models (e.g., Bekker, Wong) to estimate sinkage and contact area for given soil-tire conditions
- ✓ Explain the physical significance of pressure gradient and its impact on haul road degradation
- ✓ Design appropriate tire inflation pressures for specific soil CBR values to limit peak contact pressure below allowable bearing capacity
📖 Why This Matters
In open-pit mines, 70% of operational downtime stems from haul road failures — ruts, potholes, and shoulder erosion — all rooted in excessive or uneven tire–soil contact pressure. Understanding how pressure distributes under massive off-highway trucks (e.g., 360-ton CAT 797) isn’t just academic: it directly determines road maintenance frequency, fuel consumption (up to +15% in poorly supported conditions), and tire life (a single tire costs >$50,000). Getting this right saves millions annually and improves safety by preventing loss-of-control events on degraded surfaces.
📘 Core Principles
Contact pressure distribution arises from three interacting domains: (1) Tire structure — carcass stiffness, ply orientation, and inflation pressure govern deformation under load; (2) Soil response — governed by shear strength (c, φ), modulus (Eₛ), and density, dictating whether contact behaves elastically, plastically, or viscoelastically; (3) Interface mechanics — where real contact area (often 30–60% of apparent footprint) and pressure non-uniformity (e.g., higher at edges due to ‘bulging’) emerge. Classical Hertzian theory fails here due to large deformations; instead, semi-empirical models (Bekker’s pressure-sinkage, Wong’s analytical tire model) and modern FEM-based approaches (e.g., LS-DYNA with Drucker–Prager soil) are used. A key insight: peak pressure often occurs not at the center, but 15–25% inward from the footprint edge — a consequence of radial tire carcass tension.
📐 Bekker’s Pressure–Sinkage Relationship
Bekker’s exponential model relates vertical pressure (p) to sinkage (z) and tire–soil interaction parameters. It’s foundational for estimating contact pressure distribution shape and total sinkage — essential for predicting rolling resistance and rut depth. Used in OEM simulation tools (e.g., Caterpillar’s MineStar Terrain Modeling) and haul road design guidelines.
Bekker’s Pressure–Sinkage Equation
p(z) = (k_c + k_φ z^{n−1}) z^nEmpirical relationship between vertical pressure (p) and sinkage (z), incorporating soil cohesion (k_c), shear coefficient (k_φ), and exponent (n).
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| p | Vertical contact pressure | kPa | Normal stress at the tire–soil interface |
| z | Sinkage | m | Vertical deformation of tire into soil |
| k_c | Coefficient of cohesion | kPa | Soil strength parameter related to cohesion intercept |
| k_φ | Coefficient of shear | kPa/m^(n−1) | Soil strength parameter related to internal friction |
| n | Sinkage exponent | dimensionless | Empirically derived exponent reflecting soil compressibility |
Typical Ranges:
Sandy loam (CBR 10): n = 0.9–1.0; k_c = 50–80 kPa; k_φ = 800–1200 kPa/m^(n−1)
Compacted gravel (CBR 80+): n = 1.1–1.3; k_c = 200–400 kPa; k_φ = 3000–6000 kPa/m^(n−1)
💡 Worked Example
Problem: A 63.5-ton rigid-frame haul truck applies 317.5 kN per dual-tire assembly (total axle load). Each tire has an effective contact length L = 0.85 m, width W = 0.52 m, and measured sinkage z = 0.042 m on sandy clay (n = 1.1, k_c = 120 kPa, k_φ = 1800 kPa/m^(n)). Calculate peak pressure using Bekker’s model.
1.
Step 1: Compute nominal contact area A₀ = L × W = 0.85 × 0.52 = 0.442 m²
2.
Step 2: Apply Bekker equation: p(z) = (k_c + k_φ × z^(n−1)) × z^n = (120 + 1800 × 0.042^0.1) × 0.042^1.1
3.
Step 3: Evaluate exponents: 0.042^0.1 ≈ 0.855; 0.042^1.1 ≈ 0.0347 → p ≈ (120 + 1800×0.855) × 0.0347 ≈ (120 + 1539) × 0.0347 ≈ 1659 × 0.0347 ≈ 57.6 kPa
Answer:
The calculated pressure is 57.6 kPa — well below typical allowable bearing capacity (100–150 kPa for compacted sandy clay), confirming acceptable interface performance.
🏗️ Real-World Application
At Rio Tinto’s Pilbara iron ore operation (Australia), engineers observed premature rutting on Section 7B of the East Access Road. FEM modeling revealed localized peak contact pressures of 192 kPa under loaded CAT 797 tires — exceeding the 150-kPa CBR-derived bearing capacity of the in-situ lateritic subgrade. Remediation included increasing tire inflation from 1050 kPa to 1180 kPa (reducing sinkage by 18%), adding a 150-mm crushed rock capping layer (raising Eₛ from 25 MPa to 85 MPa), and adjusting haul cycle speeds. Post-intervention monitoring showed 73% reduction in rut depth growth rate over 6 months — validating pressure distribution modeling as a predictive maintenance tool.