ToolFusion
Calculator D3

Empirical Models: Bekker, Janosi-Hanamoto, and Reece Equations Explained

These are math formulas built from real soil-and-tire experiments—not theory—to predict how deep a farm tire sinks, how much it squishes the soil sideways, and how well it grips.

⚠️ Why It Matters

1
Inaccurate sinkage prediction
↓
2
Overestimation of traction efficiency
↓
3
Excessive dynamic axle loading
↓
4
Accelerated subsoil compaction
↓
5
Reduced root zone aeration & crop yield
↓
6
Long-term loss of field productivity

📘 Definition

Empirical models—Bekker, Janosi-Hanamoto, and Reece—are semi-physical, experimentally calibrated relationships that describe the vertical and lateral pressure distribution beneath rigid or flexible wheels operating on deformable granular or cohesive soils. They express sinkage (z), shear stress (τ), and normal stress (σ) as functions of wheel geometry, load, and soil mechanical parameters (e.g., cohesion c, friction angle φ, modulus k_c, k_φ), enabling prediction of compaction depth, rut formation, and drawbar pull without requiring full continuum mechanics solvers.

🎨 Concept Diagram

Soil surfaceTire contact patchShear flow surface (Reece)Lateral compaction zoneSinkage (z)

AI-generated illustration for visual understanding

💡 Engineering Insight

Never treat Bekker’s k_c and k_φ as universal constants—they degrade rapidly with repeated trafficking. A single pass on wet soil can reduce k_c by 30–50%; always re-fit parameters after seasonal moisture shifts or cover-crop termination. The Janosi-Hanamoto model fails catastrophically when c/k_φ < 0.03: in such cases, switch to Reece’s exponential shear decay formulation with measured shear zone width.

📖 Detailed Explanation

The Bekker model (1960) was the first widely adopted empirical framework for off-road mobility. It treats soil as a nonlinear elastic-plastic continuum and expresses vertical pressure p as p = (k_c / b + k_φ) * z^n, where b is tire width and z is sinkage. This simple power law enabled early design of military and agricultural vehicles but assumes uniform pressure distribution and ignores lateral shear effects.

Janosi-Hanamoto (1961) extended this by modeling the shear stress τ along the tire-soil contact surface as τ = c + σ tan φ * (1 − e^(−j/A)), where j is shear displacement and A is a soil shear parameter. This captures the exponential buildup of shear resistance with slip—critical for predicting drawbar pull and slippage-induced compaction. However, it assumes a planar shear surface, ignoring 3D flow.

Reece (1965) introduced geometric realism: he observed that shear failure occurs along a logarithmic spiral surface, not a plane, and derived an analytical expression for the shear zone width and stress decay using plasticity theory and experimental calibration. His model couples sinkage and shear geometry explicitly—making it essential for predicting asymmetric ruts, side-slip compaction, and inter-tyre interference in modern high-clearance sprayers and harvesters operating at ≤ 0.5 m spacing.

🔄 Engineering Workflow

Step 1
Step 1: In-situ soil characterization (Cone Index, moisture, texture, bulk density)
→
Step 2
Step 2: Tire geometry & inflation measurement (rim diameter, section width, loaded radius, actual pressure)
→
Step 3
Step 3: Controlled field sinkage and drawbar pull testing at multiple loads/speeds
→
Step 4
Step 4: Regression fitting of Bekker (z vs. P), Janosi-Hanamoto (τ vs. γ), and Reece (shear zone geometry) parameters
→
Step 5
Step 5: Cross-validation using independent test runs and rut depth monitoring (LIDAR or photogrammetry)
→
Step 6
Step 6: Integration into farm machinery control logic (e.g., real-time axle load modulation or speed governors)
→
Step 7
Step 7: Annual recalibration via post-harvest soil profile assessment

📋 Decision Guide

Rock/Field Condition Recommended Design Action
Wet clay (c > 30 kPa, k_c < 2 MPa/m, moisture > 30% w.b.) Reduce axle load by ≥40%, use ultra-low-pressure tires (≤ 60 kPa), schedule operations during drying windows
Dry sandy loam (k_φ > 150 kPa/m, n ≈ 0.9, c < 5 kPa) Increase tire inflation to 120–150 kPa for reduced sinkage; accept moderate slip (8–12%) for optimal drawbar efficiency
Compacted subsoil layer (k_c > 8 MPa/m above 30 cm depth) Use Bekker-based layered model with transition depth; avoid deep tillage unless validated by cone index profiling

📊 Key Properties & Parameters

k_c (cohesion modulus)

0.5–12 MPa/m^n (n = 1.0–1.5, commonly n=1.0 for loam to clay)

Soil parameter quantifying resistance to vertical deformation due to cohesion, derived from plate sinkage tests.

⚡ Engineering Impact:

Directly governs predicted sinkage depth under static load; low k_c in wet clays leads to rapid rutting.

k_φ (friction modulus)

10–250 kPa/m (for dry sand to compact silt loam)

Soil parameter representing resistance to vertical deformation due to internal friction, obtained from circular plate tests.

⚡ Engineering Impact:

Controls load-bearing capacity at low sinkage; high k_φ enables higher axle loads without excessive rutting.

n (sinkage exponent)

0.6–1.2 (0.8 typical for medium loam; <0.7 for peat, >1.0 for dense gravel)

Empirical exponent in Bekker’s power-law sinkage equation relating pressure to sinkage (p ∝ z^n).

⚡ Engineering Impact:

Determines nonlinearity of pressure-sinkage response—low n implies progressive softening, high n indicates stiffening behavior.

c (soil cohesion)

2–50 kPa (2 kPa for saturated sand, 35+ kPa for wet clay)

Shear strength intercept in Mohr-Coulomb failure criterion, critical for lateral stress prediction in Janosi-Hanamoto.

⚡ Engineering Impact:

Dominates lateral resistance and shear failure shape—low c causes shallow, wide shear zones and poor traction retention.

📐 Key Formulas

Bekker Sinkage Equation

p = \left(\frac{k_c}{b} + k_φ\right) z^n

Predicts average vertical pressure under a rigid wheel given sinkage, width, and soil moduli.

Variables:
Symbol Name Unit Description
p average vertical pressure Pa pressure under the rigid wheel
k_c cohesive modulus Pa/m^n soil cohesive modulus
b wheel width m contact width of the rigid wheel
k_φ frictional modulus Pa/m^n soil frictional modulus
z sinkage m vertical deformation of the wheel into the soil
n sinkage exponent empirical exponent related to soil compressibility
Typical Ranges:
Dry sandy loam
40–120 kPa
Wet clay
15–45 kPa
⚠️ z/b < 0.15 for minimal permanent compaction; p > 100 kPa risks subsoil densification below 40 cm
🔧 Open Calculator

Janosi-Hanamoto Shear Stress

τ = c + σ \tan φ \left(1 - e^{-j / A}\right)

Models shear stress development with increasing slip displacement j.

Variables:
Symbol Name Unit Description
τ Shear Stress Pa Shear stress developed at the interface
c Cohesion Pa Interfacial cohesion resistance
σ Normal Stress Pa Effective normal stress acting on the interface
φ Friction Angle rad or deg Interface friction angle
j Slip Displacement m Relative displacement (slip) along the interface
A Characteristic Slip Distance m Parameter controlling the rate of shear stress development with slip
Typical Ranges:
Low-slip traction (<5%)
10–35 kPa
High-slip rutting (>15%)
40–85 kPa
⚠️ j/A > 3.0 indicates fully mobilized shear—avoid sustained operation beyond this point
🔧 Open Calculator

Reece Shear Zone Width

w = \frac{b}{2} \left[1 + \exp\left(-\frac{2c}{k_φ z}\right)\right]

Estimates lateral extent of soil shear flow adjacent to the tire footprint.

Variables:
Symbol Name Unit Description
w Reece Shear Zone Width m Lateral extent of soil shear flow adjacent to the tire footprint
b Baseline Shear Zone Width m Maximum lateral extent of shear zone at surface
c Soil Cohesion Parameter Pa Empirical parameter related to soil cohesion
k_φ Frictional Resistance Coefficient dimensionless Coefficient representing soil frictional resistance
z Depth Below Surface m Vertical distance from soil surface
Typical Ranges:
Optimal traction (z = 0.04–0.06 m)
0.12–0.22 m
Rutting threshold (z > 0.08 m)
0.25–0.45 m
⚠️ w > 0.3 m at working speed signals irreversible lateral displacement—immediate load reduction required
🔧 Open Calculator

🏭 Engineering Example

Prairie View Research Farm, University of Nebraska-Lincoln

Not applicable — soil type: Sharpsburg silty clay loam (fine, mixed, superactive, mesic Typic Argiustolls)
c
22 kPa
n
0.78
k_c
1.3 MPa/m
k_φ
48 kPa/m
bulk_density
1.32 g/cm³
moisture_content
26.4 % w.b.

🏗️ Applications

  • Precision agriculture tire selection
  • Autonomous tractor path planning with compaction constraints
  • Regenerative farming equipment certification
  • EU Soil Thematic Strategy compliance reporting

🔧 Try It: Interactive Calculator

🔧
Tire–Soil Contact Pressure Distribution Modeling Calculator
Multi-tab calculator for Tire–Soil Contact Pressur...
Open →

📋 Real Project Case

Corn Belt No-Till Field Compaction Mitigation

1,200-acre no-till corn-soy rotation in central Illinois

Challenge: Persistent surface ruts and reduced root penetration in 2022 wet season
Corn Belt No-Till Field Compaction Mitigation Persistent surface ruts Reduced root penetration (2022 wet season) Switched to 23.1R30 singles 15% lower inflation pressure + Real-time load monitoring Peak Pressure Reduction: 28% (P₁ − P₂)/P₁ × 100 Rut Depth Prediction: 1.7 cm (Measured: 1.9 cm) 20.8R42 duals High pressure → ruts 23.1R30 single Lower pressure → less compaction ~1.2 m spacing ~0.96 m footprint
Read full case study →
Precision Rice Paddy Tire Management in Vietnam High-Capacity Sprayer Tire Optimization in Western Australia Organic Vineyard Tractor Path Planning for Minimal Compaction Cold-Climate Sugar Beet Harvest Tire Selection in Minnesota

🎨 Technical Diagrams

Soil surfaceTire footprintShear zone (Reece)
Reference line (z=0)z = 0.04 mz = 0.08 mΔz = 0.04 m

📚 References

[1]
[2]
ASABE Standards D497.7: Agricultural Machinery Management Data — Tractor and Implement Performance — American Society of Agricultural and Biological Engineers (ASABE)
[3]