Your Equation

Library / Terramechanics and off-road mobility / TER-36

Terramechanics and off-road mobility

Bekker moduli from two plate widths

Cohesive modulus kc from the pressure–sinkage coefficients B1 and B2 fitted at two plate widths b1 and b2 with a common exponent, using B = kc/b + kφ.

Ref
TER-36
Rev
v1
Field
Terramechanics and off-road mobility
Result
kc [Pa m^(1−n)]

Plates of two widths give pressure–sinkage coefficients B1 and B2 in p = Bzⁿ. With Bekker’s form B = kc/b + kφ, subtracting the two relations gives kc = (B1 − B2)/(1/b1 − 1/b2), and kφ = B1 − kc/b1. The cohesive modulus kc captures the effect of plate width, and the frictional modulus kφ the part independent of it.

Cohesive modulus kc from the pressure–sinkage coefficients B1 and B2 fitted at two plate widths b1 and b2 with a common exponent, using B = kc/b + kφ. The frictional modulus kφ = B1 − kc/b1 is in the derived rows. Keep the width convention (rectangular width, circular radius) and units with the data.

Calculator

Solve for
Pa m^(1−n)
m
0.020.19
m
0.21
Pa/mⁿ
1,0001e7
Pa/mⁿ
1,0001e7
  • Frictional modulus, kφ = B1 − kc/b1100,000 Pa/m^(n)
  • Coefficient predicted for a 0.3 m plate, kc/0.3 + kφ166,700 Pa/mⁿ
  • Cohesive share of the first plate’s coefficient, kc/(b1 B1)66.67 %

What drives the result

Each variable is moved across its typical range while the others stay at the values above. The longest bar is the key driver.

Cohesive modulus (Pa m^(1−n)) current value 20,000 Pa m^(1−n)
B1Fitted coefficient, plate 1Key driver
−39,8001,960,000
B2Fitted coefficient, plate 2
−1,940,00059,800
b1Width of plate 1
2,220380,000
b2Width of plate 2
11,100
  • The cohesive modulus is 20,000 and the frictional modulus 100,000, so the cohesive term makes up 66.7 % of the narrow plate’s coefficient.
  • The two-plate method subtracts two similar numbers, so a 5 % error in B1 or B2 can change kc by far more than 5 %; use more widths when you can.
Chart data as tables

Ranked variables

VariableTypical rangeLowest resultHighest resultSwing (Pa m^(1−n))Local stepChange in result
B1 Fitted coefficient, plate 11,000 to 1e7 Pa/mⁿ−39,8001,960,0002,000,000+10%+30 %
B2 Fitted coefficient, plate 21,000 to 1e7 Pa/mⁿ−1,940,00059,8002,000,000+10%−20 %
b1 Width of plate 10.02 to 0.19 m2,222380,000377,800+10%+22.2 %
b2 Width of plate 20.2 to 1 m11,11020,0008,889+10%−8.33 %

Sweep of Fitted coefficient, plate 2

Fitted coefficient, plate 2 (Pa/mⁿ)Cohesive modulus (Pa m^(1−n))
1,00059,800
3,16259,370
10,00058,000
31,62053,680
100,00040,000
316,200−3,246
1,000,000−140,000
3,162,000−572,500
1e7−1,940,000

Sweep one variable

Cohesive modulus (Pa m^(1−n))Hover, touch, or use the arrow keys on the chart
−2M−1.5M−1M−500k0 100010k100k1M10M Fitted coefficient, plate 2 (Pa/mⁿ), log scale now 20,000

Variables and how each one contributes

kc

Cohesive modulus

Result · Pa m^(1−n)

What the equation gives you. Everything else feeds it.

b1

Width of plate 1

m · typical 0.02 to 0.19

The width of the narrower plate, in m. A narrower first plate widens the 1/b spacing and gives a better conditioned estimate. Raising it by 10 % (to 0.11) moves the result by +22.2 %, from 20,000 to 24,440 Pa m^(1−n).

b2

Width of plate 2

m · typical 0.2 to 1

The width of the wider plate, in m. A wider second plate does the same, so well-separated widths are better. Raising it by 10 % (to 0.22) moves the result by −8.33 %, from 20,000 to 18,330 Pa m^(1−n).

B1

Fitted coefficient, plate 1

Pa/mⁿ · typical 1,000 to 1e7

The fitted coefficient from the first plate, in Pa/mⁿ. Raises kc in direct proportion to the difference B1 − B2. Raising it by 10 % (to 330,000) moves the result by +30 %, from 20,000 to 26,000 Pa m^(1−n).

B2

Fitted coefficient, plate 2

Pa/mⁿ · typical 1,000 to 1e7

The fitted coefficient from the second plate, in Pa/mⁿ. Lowers kc in direct proportion to the same difference. Raising it by 10 % (to 220,000) moves the result by −20 %, from 20,000 to 16,000 Pa m^(1−n).

Worked example with the values above

  1. Width term: 1/b1 − 1/b2 = 10 − 5 = 5 1/m.
  2. Coefficient difference: B1 − B2 = 100,000.
  3. kc = 100,000 / 5 = 20,000.

Assumptions and limits

  • The same soil state, plate or ring geometry and loading rate for each point.
  • A power-law pressure–sinkage law, a monotonic exponential shear curve and a linear Mohr–Coulomb envelope in the tested range.

Common mistakes

  • Mixing plate widths, or using width where the radius is needed.
  • Using test points from different soil states or different loading branches.
  • Using a shear point at or beyond the plateau, where the deformation length cannot be found.

Statistics and uncertainty

Soil strength and stiffness parameters from two-point methods carry the full measurement scatter of two points, which is typically 10 % to 30 %, and give no estimate of the fit uncertainty. Use a least-squares fit when more tests are available.

Where it is used

  • Fitting Bekker and Janosi shear parameters from bevameter and ring-shear tests.
  • Getting first-estimate strength parameters from a pair of simple field tests.

Maths you need

  • Bekker relation: p = (kc/b + kφ) zⁿ
  • Janosi–Hanamoto curve: τ = τmax(1 − e^(−j/K))
  • Mohr–Coulomb strength: τ = c + σ tan φ

Related formulas

Practise this formula