Skip to content

Deflection Based Nondestructive Pavement Analyses

Introduction

This topic will be oriented toward techniques which can be used to interpret nondestructive testing, NDT, data from the Falling Weight Deflectometer, FWD. There are all kinds of NDT data which can be collected on or about pavements but concentration is placed on measured surface deflections. This article is separated into three broad themes: (1) an introduction, (2) indices for project analysis, and (3) fundamentals of backcalculation. The introduction will start with an overview of some of the deflection basic parameters that have been used in the past; although, some of the NDT equipment noted is rarely used today (such as the Dynaflect). This is followed by a discussion of some easy to use models to predict layer moduli that were developed by the University of Washington for the Washington State DOT (WSDOT). A short discussion on AASHTO developed models follows along with some of the work done in South Africa and companies such as the Shell Oil Co. The second theme will cover basic parameters which can be quite useful in analyzing existing pavement structures. The last theme focuses on the basics associated with backcalculation of layer moduli. These basics do not rely on specific software but do largely involve layered elastic analysis.

Other Nondestructive Tools

The FWD is one of several tools that evaluate a pavement in place without removing material. Others commonly used are:

  • Seismic pavement analyzer. Pneumatic hammers of different sizes strike the pavement and produce vibrations of different frequencies, which sensors at fixed distances from the hammers record. From the response at each frequency the analyzer estimates layer thickness and stiffness and can detect cracks or delaminations beginning to form.
  • Ground penetrating radar (GPR). Radar pulses sent into the pavement are reflected back to an antenna, and variations in the returned signal are used to estimate individual layer thicknesses and the air void content of asphalt layers.
  • Dynamic cone penetrometer (DCP). A weight is repeatedly dropped to drive a hardened cone into unbound base or subgrade, or weakly stabilized material, and the rate of penetration indicates shear strength. It is regarded as nondestructive because the small hole does not affect the material's ability to support the pavement.
  • Nuclear density gauges and electrical (non-nuclear) density gauges. These measure the in-place density of asphalt layers with little or no disturbance of the mat.

Nondestructive tools are particularly useful for finding variability that may indicate weak areas and for checking whether load-bearing capacity is structurally adequate, though some structural properties still require laboratory testing of cores or samples.

Deflection Basin Parameters

Over the years numerous techniques have been developed to analyze deflection data from various kinds of pavement deflection equipment. A fairly complete summary of deflection basin parameters was provided by Horak at the Sixth International Conference Structural Design of Asphalt Pavement [Horak, 19871] and is shown in Table 1.

Table 1. Summary of Deflection Basin Parameters [modified from Horak, 1987
ParameterFormulaMeasuring Device
Maximum deflectionD0Benkelman Beam, Lacroux deflectometer, FWD
Radius of curvatureR = r2/2 D0(D0/Dr – 1); r = 5″Curvaturemeter
SpreadabilityS = + D1 + D2 + D3/5100/D0; D1 … D3 spaced 12″ apartDynaflect
AreaA = 6[1 + 2 (D1/D0) + 2 (D2/ D0) + D3/D0); with sensors located at 0, 1, 2, 3 ftFWD
Shape factorsF1 = (D0 -D2) / D1 and F2 = (D1 -D3) / D2FWD
Surface curvature indexSCI = D0 – Dr, where r = 12″ or r = 20″Benkelman Beam, Road Rater, FWD
Base curvature indexBCI = D24″ – D36″Road Rater
Base damage indexBDI = D12″ – D24″Road Rater
Deflection ratioQr = Dr/D0, where Dr ≈ D0/2FWD
Bending indexBI = D/a, where a = Deflection basinBenkelman Beam
Slope of deflectionSD = tan-1 (D0 – Dr )/r where r = 24 in.Benkelman Beam

All of these parameters tend to focus on four major areas:

  • Plate or center load deflection which represents the total defection of the pavement. This was obviously the first deflection parameter which came with the Benkelman Beam. It has been used for many years as the primary input for several overlay design procedures.
  • The slope or deflection differences close to the load such as Radius of Curvature (R), Shape Factor (F₁), and Surface Curvature Index (SCI). These parameters tend to reflect the relative stiffness of the bound or upper regions of the pavement section.
  • The slope or deflection differences in the middle of the basin about 11.8 in. (300 mm) to 35.4 in. (900 mm) from the center of the load. These parameters tend to reflect the relative stiffness of the base or lower regions of the pavement section.
  • The deflections toward the end of the basin. Deflections in this region relate quite well to the stiffness of the subgrade below the pavement surfacing.

The subsequent parameters to be presented were developed to provide a means of obtaining the resilient modulus values of the surfacing layers more easily or quickly than full backcalculation. In general, the success of these regression equations to predict the resilient modulus of the surfacing layers has been limited. There is a clear consensus, however, that the deflections measured beyond the primary effects of the load stress bulb relate quite well to the resilient modulus of the subgrade.

WSDOT Equations

Several researchers have developed regression equations to predict the subgrade modulus, ESGE_{SG}, from plate load and deflections measured at distances from about 24 in. (600 mm) to 48 in. (1200 mm) from the center of the plate load. Newcomb developed such regression equations to predict ESGE_{SG} as part of an overall effort to develop a mechanistic empirical overlay design procedure for WSDOT [Newcombe, 1986 2]. For two layer cases, the subgrade modulus, can be estimated from:

ESGE_{SG} = -466 + 0.00762 (P/D₃), where 3 = 3 feet from middle of the load plate (Eq. 1)

ESGE_{SG} = -198 + 0.00577 (P/D₄), where 4 = 4 feet from middle of the load plate (Eq. 2)

ESGE_{SG} = -371 + 0.00671 (2P/(D₃ + D₄)), (Eq. 3)

and for three layer cases:

ESGE_{SG} = -530 + 0.00877 (P/D₃), (Eq. 4

ESGE_{SG} = -111 + 0.00577 (P/D₄), (Eq. 5)

ESGE_{SG} = -346 + 0.00676 (2P/(D₃ + D₄)) (Eq. 6)

The R² ≈ 99% for all equations and the sample sizes were 180 (two layer case) and 1,620 (three layer case).

These equations were developed from generated data (ELSYM5) with the following levels of input data: (Tables 2 and 3)

Table 2. Two Layer Cases
Load, P, lb (kN)Surface Thickness, hAC, in. (mm)Surface Modulus, EAC, psi (MPa)Subgrade Modulus, ESG, psi (MPa)
5,000 (22)2 (50)2,000,000 (13800)50,000 (345)
10,000 (44)6 (150)500,000 (3450)30,000 (207)
15,000 (67)12 (300)100,000 (690)10,000 (69)
18 (450)5,000 (35)
2,500 (17)
Table 3. Three Layer Cases
Load, P, lb (kN)Surface Thickness, hAC, in. (mm)Base Thickness, hB, in (mm)Surface Modulus, EAC, psi (MPa)Base Modulus, EB, psi (MPa)Subgrade Modulus, ESG, psi (MPa)
5,000 (22)2 (50)4 (100)2,000,000 (13800)100,000 (690)50,000 (345)
10,000 (44)6 (150)10 (250)500,000 (3450)50,000 (345)30,000 (207)
15,000 (67)12 (300)18 (450)100,000 (690)30,000 (207)10,000 (69)
10,000 (69)5,000 (35)
2,500 (17)

(assumed that load applied on a 11.8 in. (300 mm) diameter load plate)

From this generated data (no rigid base), regression equations were also developed for estimating the surface modulus (AC) for a two layer case (for example a "full-depth" pavement):

log EACE_{AC} = -0.53740 - 0.95144 log₁₀ ESGE_{SG} -1.2118(hACh_{AC})0.33{}^{0.33} + 1.78046 log₁₀ PA₁/D₀²) (Eq. 7)

where R² = 0.83

For a three layer case, equations were developed for both EACE_{AC} and EBE_{B} as follows:

If both EACE_{AC} and EBE_{B} unknown:

log EACE_{AC} = -4.13464 + 0.25726 (5.9/hACh_{AC}) + 0.92874(5.9/hBh_{B})0.5{}^{0.5} - 0.69727 (hACh_{AC}/hBh_{B})0.5{}^{0.5} - 0.96687 log₁₀ ESGE_{SG} + 1.88298 log₁₀ (PA₁/D₀²) (Eq. 8)

where R² = 0.78.

log EBE_{B} = 0.50634 + 0.03474 (5.9/hACh_{AC}) + 0.12541(5.9/hBh_{B})0.5{}^{0.5} - 0.09416 (hACh_{AC}/hBh_{B})0.5{}^{0.5} + 0.51386 log ESGE_{SG} + 0.25424 log₁₀ (PA₁/D₀²) (Eq. 9)

where R² = 0.70.

The following variables were used in the equations shown above:

P = applied load (lbs) on a 11.8 in. (300 mm) plate,

hACh_{AC} = surface course thickness (in.),

hBh_{B} = base course thickness (in.),

EACE_{AC} = surface course modulus (psi),

EBE_{B} = base course modulus (psi),

ESGE_{SG} = subgrade modulus (psi),

D₀ = deflection under center of applied load (in.),

D0.67D_{0.67} = deflection at 8 in. (0.67 ft) from center of applied load (in.),

D₁ = deflection at 1 ft from center of applied load (in.),

D₂= deflection at 2 ft from center of applied load (in.),

D₃ = deflection at 3 ft from center of applied load (in.),

D₄ = deflection at 4 ft from center of applied load (in.), and

A₁ = approximate area under deflection basin out to 3 ft = 2 + D0.67D_{0.67}) + (D0.67D_{0.67} + D₁) + 3(D₁ + D₂) + 3(D₂ + D₃) = 4D₀ + 6D0.676D_{0.67} + 8D₁ + 12D₂ + 6D₃

AASHTO Equations

Witczak presented a similar regression equation in Part III of the 1986 AASHTO Guide for Design of Pavement Structures [AASHTO, 19863]. That equation is in the form:

ESGE_{SG} = (P)(SfS_{f})/(DrD_{r})(r) (Eq. 10)

where

P = plate load (lb),

SfS_{f} = subgrade modulus prediction factor,

= 0.2686 for µ = 0.50

= 0.2792 for µ = 0.45

= 0.2892 for µ = 0.40

= 0.2874 for µ = 0.35

= 0.2969 for µ = 0.30

DrD_{r} = pavement surface deflection measured at r distance, from the load, and

r = distance from the load to DrD_{r} (in.).

Using an SfS_{f} value of 0.2892 where the Poisson's ratio is 0.40, the equation reduces to the following equations for deflections measured at 2 ft (610 mm), 3 ft (914 mm), and 4 ft. (1,219 mm).

ESGE_{SG} = 0.01205(P/D₂) (Eq. 11)

ESGE_{SG} = 0.00803(P/D₃) (Eq. 12)

ESGE_{SG} = 0.00603(P/D₄) (Eq. 13)

In the AASHTO Guide, detailed procedures are provided to insure that the deflections used to determine ESGE_{SG} are outside the pressure bulb from the test load. To most accurately represent the subgrade stiffness; however, the deflections closest to the load without being directly effected by the pressure bulb should be used. Stiff underlying layers will have the greatest effect on deflections furthermost from the test load. For example in cases where total surfacing depths are around 11.8 in. (300 mm), deflections taken around 23.6 in. (600 mm) should be used to determine ESGE_{SG}.

A NCHRP study [Darter et al., 19914] was used to revise Part III of the AASHTO Pavement Guide [AASHTO, 19863], it is recommended that the following equation be used to solve for subgrade modulus:

MRM_{R} = P(1 - µ2)/(π)(DrD_{r})(r) (Eq. 14)

where

MRM_{R} = backcalculated subgrade resilient modulus (psi),

P = applied load (lbs),

DrD_{r} = pavement surface deflection a distance r from the center of the load plate (inches), and

r = distance from center of load plate to DrD_{r} (inches).

Using a Poisson's ratio of 0.40, Equation 14 reduces to

MRM_{R} = 0.01114 (P/D₂) (Eq. 15)

MRM_{R} = 0.00743 (P/D₃) (Eq. 16)

MRM_{R} = 0.00557 (P/D₄) (Eq. 17)

for sensor spacings of 2 ft (610 mm), 3 ft (914 mm), and 4 ft (1219 mm). If a Poisson's ratio of 0.45 is used instead for the same sensor spacings, the equations become:

MRM_{R} = 0.01058(P/D₂) (Eq. 18)

MRM_{R} = 0.00705 (P/D₃) (Eq. 19)

MRM_{R} = 0.00529 (P/D₄) (Eq. 20)

Darter et al. [19914] recommends that the deflection used for subgrade modulus determination should be taken at a distance at least 0.7 times r/aea_{e} where r is the radial distance to the deflection sensor and aea_{e} is the radial dimension of the applied stress bulb at the subgrade "surface." The aea_{e} dimension can be determined from the following:

aea_{e} =

Darter et al. [1991<pi-cite data-n=
Figure 1.

where

aea_{e} = radius of stress bulb at the subgrade-pavement interface,

a = NDT load plate radius (inches),

D = total thickness of pavement layers (inches)

EPE_{P} = effective pavement modulus (psi), and

MRM_{R} = backcalculated subgrade resilient modulus.

For "thin" pavements (such as Section A, Figure 3), aea_{e} ≈ 15 in.; "medium" to "thick" pavements (Sections B and C, Figure 3), aea_{e} ≈ 26 to 33 in. Thus, the minimum r is usually 24 to 36 in. (recall r > 0.7 (aea_{e})).

South African Equations

A similar relation between deflections taken at 2000 mm (78.7 in.) was shown by Horak [1987]1:

log₁₀ ESGE_{SG} = 9,727 - 0.989 log₁₀ δ₂₀₀₀

where

ESGE_{SG} = subgrade elastic modulus (Pa), (same as ESGE_{SG}), and

δ₂₀₀₀ = deflection at a distance of 2000 mm from the (point) of loading (µm).

Additional Modulus Estimates

Subgrade Soils

  • Table 4 is used to show typical subgrade moduli as listed by Chou et al. [19895].
  • CBR Correlation- A classic correlation that has been widely used is:

ESGE_{SG} (psi) = 1,500 (CBR)

or

ESGE_{SG} (MPa) = 10 (CBR)

where

CBR = California Bearing Ratio

Note: Limited to fine-grained soils with a CBR of 10 or less.

Table 4. Typical Values of Subgrade Moduli [after Chou et al., 1989
Climate Condition
DryWet – No FreezeWet - Freeze
FrozenUnfrozen
Materialpsi (MPa)psi (MPa)psi (MPa)psi (MPa)
Clay15,000 (103)6,000 (41)6,000 (41)50,000 (345)
Silt15,000 (103)10,000 (69)5,000 (34)50,000 (345)
Silty or Clayey Sand20,000 (138)10,000 (69)5,000 (34)50,000 (345)
Sand25,000 (172)25,000 (172)25,000 (172)50,000 (345)
Silty or Clayey Gravel40,000 (276)30,000 (207)20,000 (138)50,000 (345)
Gravel50,000 (345)50,000 (345)40,000 (276)50,000 (345)

Granular Materials Shell Method [Chou et al., 1989 5; Smith and Witczak, 1981 6]

The granular base modulus is a function of the subgrade modulus and base course thickness.

EBE_{B} = 0.2(25.4h₂)0.45{}^{0.45}ESGE_{SG}

where

EBE_{B} = base modulus (psi),

h₂ = base thickness (in.)

ESGE_{SG} = subgrade modulus (psi)

Example

If

ESGE_{SG} = 15,000 psi and the base thickness = 10 in., then

EBE_{B} = 0.2(25.4 (10))0.45{}^{0.45}(15,000) ≈ 36,000 psi (250 MPa)

Stabilized Base Materials

From Chou et al. [22] 5:

  • Asphalt stabilized materials: use information available for asphalt concrete.
  • Lime-stabilized materials: a modulus range of 100,000 to 1,000,000 psi (690 to 6,895 MPa) is suggested.
  • Cement stabilized materials: a range of 300,000 to 4,000,000 psi (2,070 to 27,600 MPa) is suggested.

Asphalt Concrete

From Chou et al. [19895]:

  • Absolute minimum and maximum values for AC are 50,000 to 2,000,000 psi (345 to 13,790 MPa).
  • To estimate the lowest and highest AC moduli, use expected high and low temperatures along with the Asphalt Institute equation [ASTM, 19847] which is of the form:

| E* | = f (P₂₀₀, f, VvV_{v}, η_(70˚F), T, PacP_{ac})

where

| E* | = dynamic modulus (stiffness of asphalt concrete),

P₂₀₀ = percent aggregate passing No. 200 sieve,

f = frequency of loading,

VvV_{v} = percent air voids,

η_(70˚F) = original absolute viscosity used in mix at 70F (21.1C)

T = temperature, and

PacP_{ac} = asphalt content, by weight of mix.

Chou et al. [19895] use default values of

P₂₀₀ = 6%

f = 25 Hz

VvV_{v} = 7%

η_(70˚F) = 10⁶ poises

PacP_{ac} = 5%

Indices For Project Analysis

Introduction

Over the years it has found that the use of selected indices and algorithms provide a "good picture" of the relative conditions found throughout a project. This picture is very useful in performing backcalculation and may at times be used by themselves on projects with large variations in surfacing layers.

WSDOT is currently using deflections measured at the center of the test load combined with Area values and ESGE_{SG} computed from deflections measured at 24 in. (610 mm) presented in a linear plot to provide a visual picture of the conditions found along the length of any project (as illustrated in Figure 1).

WSDOT Non-destructive Pavement Testing
Figure 2. WSDOT Non-destructive Pavement Testing

Area Parameter

The Area value represents the normalized area of a slice taken through any deflection basin between the center of the test load and 3 ft (914 mm). By normalized, it is meant that the area of the slice is divided by the deflection measured at the center of the test load, D₀. Thus the Area value is the length of one side of a rectangle where the other side of the rectangle is D₀. The actual area of the rectangle is equal to the area of the slice of the deflection basin between 0 and 3 ft (914 mm).

The original Area equation is:

A = 6(D₀ + 2D₁ + 2D₂ + D₃)/D₀

where

D₀ = surface deflection at center of test load,

D₁ = surface deflection at 1 ft (305 mm),

D₂ = surface deflection at 2 ft (610 mm), and

D₃ = surface deflection at 3 ft (914 mm).

The approximate metric equivalent of this equation is:

A = 150(D₀ + 2D₃₀₀ + 2D₆₀₀ + D₉₀₀)/D₀

where

D₀ = deflection at center of test load,

D₃₀₀ = deflection at 300 mm,

D₆₀₀ = deflection at 600 mm, and

D₉₀₀ = deflection at 900 mm.

Figure 2 shows the development of the normalized area for the Area value using the Trapezoidal Rule to estimate area under a curve.

Computing an Area Parameter
Figure 3. Computing an Area Parameter

The basic Trapezoidal Rule is:

K = h (¹/₂y₀ + y₁ + y₂ + ¹/₂y₃)

where

K = any planar area,

y₀ = initial chord,

y₁, y₂ = immediate chords,

y₃ = last chord, and

h = common distance between chords.

Thus, to estimate the area of a deflected basin using D₀, D₁, D₂, and D₃, and h = 6 in., then:

K = 6 (D₀ + 2D₁ + 2D₂ + D₃)

Further, normalize by dividing by D₀:

Area = (6/D₀)*(D₀ + 2D₁ + 2D₂ + D₃)

or

Area = 6 (1 +2D₁/D₀ + 2D₂/D₀ + D₃/D₀)

Thus, since we normalized by D₀, the Area Parameter's unit of measure is inches (or mm) not in² or mm² as one might expect.

The maximum value for Area is 36.0 (915 mm) and occurs when all four deflection measurements are equal (not likely to occur) as follows:

If, D₀ = D₁ = D₂ = D₃

Then, Area = 6(1 + 2 + 2 + 1) = 36.0 in.

For all four deflection measurements to be equal (or nearly equal) would indicate an extremely stiff pavement system (like portland cement concrete slabs or thick, full-depth asphalt concrete.)

The minimum Area value should be no less than 11.1 in. (280 mm).

This value can be calculated for a one-layer system which is analogous to testing (or deflecting) the top of the subgrade (i.e., no pavement structure). Using appropriate equations, the ratios of

D₁/D₀ , D₂/D₀, D₃/D₀

always result in 0.26, 0.125, and 0.083, respectively. Putting these ratios in the Area equation results in:

Area = 6(1+ 2(0.26) + 2(0.125) + 0.083) = 11.1 in.

Further, this value of Area suggests that the elastic moduli of any pavement system would all be equal (e.g., E₁ = E₂ = E₃ = …). This is highly unlikely for actual, in-service pavement structures. Low area values suggest that the pavement structure is not much different from the underlying subgrade material (this is not always a bad thing if the subgrade is extremely stiff — which doesn't occur very often).

Typical Area values were computed for pavement Sections A, B, and C (Figure 3) and are shown in Table 5 (along with the calculated surface deflections (D₀, D₁, D₂, D₃)). Table 6 provides a general guide in the use of Area values obtained from FWD pavement surface deflections. For Figure 3, the initial layer moduli are HMA = 500,000 psi, crushed stone base = 25,000 psi, and fine-grained subgrade = 7,500 psi.

"Typical" Pavement Sections
Figure 4. "Typical" Pavement Sections
Table 5. Estimates of Area from Pavement Sections - Cases A, B, and C
Pavement CasesPavement Surface Deflections, inchesArea
D⁰D₁D₂D₃in.(mm)
Standard Pavement (HMA @ 500,000 psi, base course @ 25,000 psi, and subgrade @ 7,500 psi)
Section A (thin)0.0480.0260.0140.00917.1(434)
Section B (med)0.0270.0200.0140.01023.3(592)
Section C (thick)0.0180.0150.0120.00927.0(686)
Stabilize Subgrade (upper 6 in. of subgrade increased from 7,500 to 50,000 psi)
Section A (thin)0.0360.0200.0130.00918.5(470)
Section B (med)0.0230.0170.0120.00923.5(597)
Section C (thick)0.0160.0130.0110.00927.4(696)
Asphalt Treated Base (increase base course modulus from 25,000 to 500,000 psi)
Section A (thin)0.0210.0180.0130.01026.6(676)
Section B (med)0.0140.0120.0100.00928.7(729)
Section C (thick)0.0120.0110.0090.00830.0(762)
Moisture Sensitivity (decrease HMA modulus from 500,000 to 200,000 psi)
Section A (thin)0.0530.0260.0140.00916.1(409)
Section B (med)0.0330.0220.0140.00920.7(526)
Section C (thick)0.0240.0180.0130.01024.0(610)
Table 6. Trends of D₀ and Area Values
FWD Based ParameterGeneralized Conclusions*
AreaMaximum Surface Deflection (D₀)
LowLowWeak structure, strong subgrade
LowHighWeak structure, weak subgrade
HighLowStrong structure, strong subgrade
HighHighStrong structure, weak subgrade

Naturally, exceptions can occur.

Table 7. Typical Area Values
PavementArea
in.(mm)
PCCP24-33(610-840)
“Sound” PCC*29-32(740-810)
BST flexible pavement (relatively thin structure)15-17(380-430)
Thick ACP (≥ 4.2 in. HMA)21-30(530-760)
Weak BST12-15(300-380)
Thin ACP (< 4.2 in. HMA)16-21(410-530)

from 1993 AASHTO Guide, p. III-117

Fundamentals of Backcalculation

Introduction

A method which has received extensive development work is the "backcalculation" method. This is essentially a mechanistic evaluation of pavement surface deflection basins generated by various pavement deflection devices. Backcalculation is one where measured and calculated surface deflection basins are matched (to within some tolerable error) and the associated layer moduli required to achieve that match are determined. The backcalculation process is usually iterative and normally done with software. An illustration of the backcalculation process is shown in Figure 4.

The Backcalculation Process - Matching Measured and Calculated Deflection Basins
Figure 5. The Backcalculation Process - Matching Measured and Calculated Deflection Basins

Typical Flowchart

A basic flowchart which represents the fundamental elements in all known backcalculation programs is shown as Figure 5. This flowchart was patterned after one shown by Lytton (19898). Briefly, these elements include:

  1. Measured deflections - Includes the measured pavement surface deflections and associated distances from the load.
  2. Layer thicknesses and loads - Includes all layer thicknesses and load levels for a specific test location.
  3. Seed moduli - The seed moduli are the initial moduli used in the computer program to calculate surface deflections. These moduli are usually estimated from user experience or various equations.
  4. Deflection calculation - Layered elastic computer programs are generally used to calculate a deflection basin.
  5. Error check - This element simply compares the measured and calculated basins. There are various error measures which can be used to make such comparisons (more on this in a subsequent paragraph in this section).
  6. Search for new moduli - Various methods have been employed within the various backcalculation programs to converge on a set of layer moduli which produces an acceptable error between the measured and calculated deflection basins.
  7. Controls on the range of moduli - In some of the backcalculation programs, a range (minimum and maximum) of moduli are selected or calculated to prevent program convergence to unreasonable moduli levels (could be too high or low).
Common Elements of Backcalculation Programs [modified after Lytton, 1989
Figure 6. Common Elements of Backcalculation Programs [modified after Lytton, 1989

Measure of Convergence

Measure of Deflection Basin Convergence

In backcalculating layer moduli, the measure of how well the calculated deflection basin matches (or converges to) the measured deflection basin was previously described as the "error check." This is also referred to as the "goodness of fit" or "convergence error."

The primary measure of convergence will be Root Mean Square (RMS) or Root Mean Square Error (RMSE) and is used in the EVERCALC program. Both terms are taken to be the same.

RMS (%)=(1nd∑i=1n(dci−dmidmi)2)(100)\mathrm{RMS}\ (\%) = \left(\sqrt{\frac{1}{n_d}\sum_{i=1}^{n}\left(\frac{d_{ci} - d_{mi}}{d_{mi}}\right)^2}\right)(100)

where

RMS = root mean square error,

dcid_{ci} = calculated pavement surface deflection at sensor i,

dmid_{mi} = measured pavement surface deflection at sensor i, and

ndn_{d} = number of deflection sensors used in the backcalculation process.

Table 8. Example RMS Calculation
n_dDeflections (mils)
MeasuredCalculated
1 (0″)5.074.90
2 (8″)4.323.94
3 (12″)3.673.50
4 (18″)2.993.06
5 (24″)2.402.62
6 (36″)1.691.86
7 (60″)1.010.95
RMS (%)=[17(4.90−5.071.01)2+⋯+(0.95−1.011.01)2]=6.9 % (very high)\begin{aligned} \mathrm{RMS}\ (\%) &= \left[\frac{1}{7}\left(\frac{4.90 - 5.07}{1.01}\right)^2 + \dots + \left(\frac{0.95 - 1.01}{1.01}\right)^2\right] \\ &= 6.9\ \%\ \text{(very high)} \end{aligned}

Generally an adequate range of RMS is 1 to 2 percent.

Depth of Stiff Layer

Research provided at least two approaches for estimating the depth to stiff layer [Rohde and Scullion, 19909; Hossain and Zaniewski, 199110]. The approach used by Rohde and Scullion [19909] will be summarized below. There are two reasons for this selection: (1) initial verification of the validity of the approach is documented, and (2) the approach is used in MODULUS and EVERCALC — backcalculation programs widely used in the U.S.

(a) Basic Assumptions and Description

A fundamental assumption is that the measured pavement surface deflection is a result of deformation of the various materials in the applied stress zone; therefore, the measured surface deflection at any distance from the load plate is the direct result of the deflection below a specific depth in the pavement structure (which is determined by the stress zone). This is to say that only that portion of the pavement structure which is stressed contributes to the measured surface deflections. Further, no surface deflection will occur beyond the offset (measured from the load plate) which corresponds to the intercept of the applied stress zone and the stiff layer (the stiff layer modulus being 100 times larger than the subgrade modulus). Thus, the method for estimating the depth to stiff layer assumes that the depth at which zero deflection occurs (presumably due to a stiff layer) is related to the offset at which a zero surface deflection occurs. This is illustrated in Figure 6 where the surface deflection DcD_{c} is zero.

Illustration of Zero Deflection Due to a Stiff Layer
Figure 7. Illustration of Zero Deflection Due to a Stiff Layer

An estimate of the depth at which zero deflection occurs can be obtained from a plot of measured surface deflections and the inverse of the corresponding offsets (1/r). This is illustrated in Figure 7. The middle portion of the plot is linear with either end curved due to nonlinearities associated with the upper layers and the subgrade. The zero surface deflection is estimated by extending the linear portion of the D vs. 1/r plot to a D = 0, the 1/r intercept being designated as r₀. Due to various pavement section-specific factors, the depth to stiff layer cannot be directly estimated from r₀ — additional factors must be considered. To do this, regression equations were developed based on BISAR computer program generated data for the following factors and associated values:

Load = P = 9000 lbs (40 kN) (only load level considered)

Moduli ratios:

E₁/ESGE_{SG} = 10, 30, 100

E₂/ESGE_{SG} = 0.3, 1.0, 3.0, 10.0

ErigidE_{rigid}/ESGE_{SG} = 100

Thickness levels:

T₁ = 1, 3, 5, and 10 in. (25, 75, 125, 250 mm)

T₂ = 6, 10, and 15 in. (150, 250, 375 mm)

B = 5, 10, 15, 20, 25, 30, and 50 ft (1.5, 3.0, 4.5, 6.0, 7.5, 9.0, 15.0 m)

where

EiE_{i} = elastic modulus of layer i,

TiT_{i} = thickness of layer i,

B = depth of the rigid (stiff) layer measured from the pavement surface (ft).

Plot of Inverse of Deflection Offset vs. Measured Deflection (Mahoney et al., 1991)
Figure 8. Plot of Inverse of Deflection Offset vs. Measured Deflection (Mahoney et al., 1991)

The resulting regression equations follow in (b).

(b) Regression Equations

Four separate equations were developed for various HMA layer thicknesses. The dependent variable is 1/B and the independent variables are r₀ (and powers of r₀) and various deflection basin shape factors such as SCI, BCI, and BDI.

(i) HMA less than 2 in. (50 mm) thick

1/B = 0.0362 - 0.3242 (r₀) + 10.2717 (r₀²) - 23.6609 (r₀³) - 0.0037 (BCI)

R² = 0.98

(ii) HMA 2 to 4 in. (50 to 100 mm) thick

1/B = 0.0065 + 0.1652 (r₀) + 5.4290 (r₀²) - 11.0026 (r₀³) - 0.0004 (BDI)

R² = 0.98

(iii) HMA 4 to 6 in. (100 to 150 mm) thick

1/B = 0.0413 + 0.9929 (r₀) - 0.0012 (SCI) + 0.0063 (BDI) - 0.0778 (BCI)

R² = 0.94

(iv) HMA greater than 6 in. (150 mm) thick

1/B = 0.0409 + 0.5669 (r₀) + 3.0137 (r₀²) + 0.0033 (BDI) - 0.0665 log (BCI)

R² = 0.97

where

r₀ = 1/r intercept (extrapolation of the steepest section of the D vs. 1/r plot) in units of 1/ft

SCI = D₀ - D_(12") (D₀ - D305mmD_{305mm}), Surface Curvature Index,

BDI = D_(12") - D_(24") (D₃₀₅ - D610mmD_{610mm}), Base Damage Index,

BCI = D_(24") - D_(36") (D₆₁₀ - D914mmD_{914mm}), Base Curvature Index,

DiD_{i} = surface deflections (mils) normalized to a 9,000 lb (40 kN) load at an offset i in feet.

(c) Example

Use typical deflection data to estimate B (depth to stiff layer). The drillers log suggests a stiff layer might be encountered at a depth of 198 in. (5.0 m).

(i) First, calculate normalized deflections (9,000 lb (40 kN) basis).

Table 9. Normalized Deflections
Load LevelDeflections (mils)
D₀D₈D₁₂D₁₈D24D36D60
9,512 lb5.074.323.672.992.401.691.01
9,000 lb (normalized)4.764.043.442.802.261.590.95
6,534 lb3.282.692.331.881.561.090.68

(ii) Second, estimate r₀. Plot DrD_{r} vs. 1/r (refer to Figure 7) :

Table 10. Table of Values D_r vs. 1/r
D_r (mils)r1/r(1/ft)
4.760″—
4.048″1.50
3.4412″1.00
2.8018″0.67
2.2624″0.50
1.5936″0.33
0.9560″0.20

where all DrD_{r} normalized to 9,000 lb (40 kN)

(iii) Third, use regression equation in (b)(iv) (for HMA = 7.65 in. (194 mm)) to calculate B:

1/B = 0.0409 + 0.5669 (r₀) + 3.0137 (r₀²) + 0.0033 (BDI) - 0.0665 log (BCI)

where

r₀ = 1/r intercept (refer to Figure 8)

≈ 0 (used steepest part of deflection basin which is for Sensors at 36 and 60 inches)

BDI = D_(12") - D_(24") = 3.44 - 2.26 = 1.18 mils

BCI = D_(24") - D_(36") = 2.26 - 1.59 = 0.67 mils

1/B = 0.0409 + 0.5669 (0) + 3.0137 (0²) + 0.0033 (1.18) - 0.0665 log (0.67)

= 0.0564

1/B = 10.0564 = 17.7 feet (213 inches or 5.4 m)

This values agrees fairly well with "expected" stiff layer conditions at 16.5 feet (5.0 m) as indicated by the drillers log even though the 1/r value equaled zero.

Plot of Inverse of Deflection Offset vs. Measured Deflection
Figure 9. Plot of Inverse of Deflection Offset vs. Measured Deflection

Example of Depth to Stiff Layer Estimates - Road Z-675 (Sweden)

(a) Overview:

This road located in south central Sweden is used to illustrate calculated and measured depths to stiff layers (the stiff layer apparently being rock for the specific road).

(b) Measurement of Measured Depth:

The depth to stiff layer was measured using borings (steel drill) and a mechanical hammer. The hammer was used to drive the drill to "refusal." Thus, the measured depths could be to bedrock, a large stone, or hard till (glacially deposited material). Further, the measured depths were obtained independently of the FWD deflection data (time difference of several years).

(c) Deflection Measurements:

A KUAB 50 FWD was used to obtain the deflection basins. All basins were obtained within ± 5 m (16 ft) of a specific borehole. The deflection sensor locations were set at 0, 200, 300, 450, 600, 900, and 1200 mm (0, 7.9, 11.8, 17.7, 23.6, 35.4, and 47.2 in.) from the center of the load plate.

(d) Calculations:

The equations described in Section 3.3.4 were used to calculate the depth to stiff layer. Since the process requires a 9000 lb (40 kN) load and 1 ft (305 mm) deflection sensor spacings, the measured deflections were adjusted linearly according to the ratio of the actual load to a 9000 lb (40 kN) load.

(e) Results:

The results of this comparison are shown in Figure 9. Given all the uncertainties concerning the measured depths, the agreement is quite good.

Plot of Measured and Calculated Depths to Stiff Layer (Road Z-675, Sweden)
Figure 10. Plot of Measured and Calculated Depths to Stiff Layer (Road Z-675, Sweden)

Verification of Backcalculation Results

There have been and, undoubtedly, will be a number of attempts to verify that backcalculated layer moduli are "reasonable." One way to do this is to compare measured, in-situ strains to calculated (theoretical) strains based on layered elastic analyses.

Two studies have done this. The first study summarized (Winters, 199311) was conducted at the PACCAR Technical Center located at Mt. Vernon, Washington. Horizontal strains were measured on asphalt concrete instrumented cores (inserted into the original asphalt concrete surface material). The HMA layer was approximately 5.5 in. thick over a 13.0 in. thick base course. The WSDOT FWD was used to apply a known load to induce the strain response (measured strain). The calculated strains were determined by use of the EVERCALC program (which used FWD deflection basins to backcalculate layer moduli which, in turn, were used to calculate the corresponding strains.)

Figures 10 and 11 show the measured and calculated strains in the two horizontal directions at the bottom of the AC layer (longitudinal and transverse). This specific testing was conducted during February 1993. In general, the agreement is good.

Measured vs. Calculated Strain for Axial Core Bottom Longitudinal Gauges - February 1993 FWD Testing (from Winters, 1993)
Figure 11. Measured vs. Calculated Strain for Axial Core Bottom Longitudinal Gauges - February 1993 FWD Testing (from Winters, 1993)
Measured vs. Calculated Strain for Axial Core Bottom Transverse Gauges - February 1993 FWD Testing (from Winters, 1993)
Figure 12. Measured vs. Calculated Strain for Axial Core Bottom Transverse Gauges - February 1993 FWD Testing (from Winters, 1993)

The second study was reported by Lenngren [199012]. That work was conducted for RST Sweden and used backcalculated layer moduli (from a modified version of the EVERCALC program) to estimate strains at the bottom of two thicknesses of HMA (3.1 and 5.9 in. thick). The pavement structures were located in Finland. Tables 10 and 11 are used as summaries. The backcalculated strains are generally within ± 5 percent of the measured strains.

Table 11. Backcalculated and Measured Tensile Strains - 3.1 in. (80 mm) HMA Section (1990)
Time of Day (a.m. or p.m.)Tensile Strain Bottom of HMA (x 10-⁶)
Backcalculated*Measured% Difference
a.m.119123-3
a.m.119122-2
a.m.7464.9+14
a.m.6064.7-8
p.m.284292-3
p.m.284283~0
p.m.167159+5
p.m.167158+6
p.m.8784.8+2
p.m.8184.2-4
  • Backcalculation process used sensors @ D₀, D₃₀₀, D₆₀₀, D₉₀₀, and D₁₂₀₀
Table 12. Backcalculated and Measured Tensile Strains —5.9 in. (150 mm) AC Section (1990)
Time of Day (a.m. or p.m.)Tensile Strain Bottom of AC (x 10-⁶)
Backcalculated*Measured% Difference
a.m.6669.5-6
a.m.7169.0+3
a.m.6868.7-1
a.m.3834.7+9
a.m.127130-2
a.m.119130-8
p.m.178185-4
p.m.182183-1
p.m.10495.9+8
p.m.5148.0+6
p.m.5648.5+14

Teaching this? Figures and worked examples are free to use with credit. Course reading lists

Cite this page

APA

Pavement Interactive. (2026, November 1). Deflection based nondestructive pavement analyses. HeadLight. Retrieved [the date you read it], from https://pavementinteractive.org/reference-desk/pavement-management/pavement-evaluation/deflection-based-nondestructive-pavement-analyses/

ASCE

Pavement Interactive. (2026). “Deflection Based Nondestructive Pavement Analyses.” HeadLight, <https://pavementinteractive.org/reference-desk/pavement-management/pavement-evaluation/deflection-based-nondestructive-pavement-analyses/> ([the date you read it]).

BibTeX

@misc{pavint-deflection-based-nondestructive-pavement-analyses,
  author = {{Pavement Interactive}},
  title = {Deflection Based Nondestructive Pavement Analyses},
  year = {2026},
  publisher = {HeadLight},
  howpublished = {\url{https://pavementinteractive.org/reference-desk/pavement-management/pavement-evaluation/deflection-based-nondestructive-pavement-analyses/}},
  note = {Accessed [date]}
}
This address will not change. The citation records the date you read it.
Sources (12)
  1. Horak (1987). "The Use of Surface Deflection Basin Measurements in the Mechanistic Analysis of Flexible Pavements," Proceedings, Vol. 1, Sixth International Conference Structural Design of Asphalt Pavements, University of Michigan, Ann Arbor, Michigan, USA, 1987. ↩ ↩
  2. Newcombe (1986). Newcomb, D. E., "Development and Evaluation of Regression Methods to Interpret Dynamic Pavement Deflections," Ph.D. Dissertation, Department of Civil Engineering, University of Washington, Seattle, Washington, 1986. ↩
  3. AASHTO (1986). AASHTO Guide for Design of Pavement Structures, American Association of State Highway and Transportation Officials, Washington, D.C., 1986. ↩ ↩
  4. Darter et al. (1991). Darter, M.I., Elliott, R.P., and Hall, K.T., "Revision of AASHTO Pavement Overlay Design Procedure," Project 20-7/39, National Cooperative Highway Research Program, Transportation Research Board, Washington, D.C., September 1991. ↩ ↩
  5. Chou et al. (1989). Chou, Y. J., Uzan, J., and Lytton, R. L., "Backcalculation of Layer Moduli from Nondestructive Pavement Deflection Data Using the Expert System Approach," Nondestructive Testing of Pavements and Backcalculation of Moduli, ASTM STP 1026, American Society for Testing and Materials, Philadelphia, 1989, pp. 341 - 354. ↩ ↩ ↩ ↩ ↩
  6. Smith and Witczak (1981). Smith. B. E. and Witczak, M. W., "Equivalent Granular Base Moduli: Prediction," Journal of Transportation Engineering, Vol. 107, TE6, American Society of Civil Engineering, November 1981. ↩
  7. ASTM (1984). "Standard Method of Indirect Tension Test for Resilient Modulus of Bituminous Mixtures," D4123-82, Volume 04.03, American Society for Testing and Materials, Philadelphia, PA, 1984. ↩
  8. Lytton (1989). Lytton, R.L., "Backcalculation of Pavement Layer Properties," Nondestructive Testing of Pavements and Backcalculation of Moduli, ASTM STP 1026, American Society for Testing and Materials, Philadelphia, 1989, pp. 7-38. ↩
  9. Rohde and Scullion (1990). Rohde, G.T., and Scullion, T., "MODULUS 4.0: Expansion and Validation of the MODULUS Backcalculation System," Research Report No. 1123-3, Texas Transportation Institute, Texas A&M University System, College Station, Texas, November 1990. ↩ ↩
  10. Hossain and Zaniewski (1991). Hossain, A.S.M., and Zaniewski, J.P., "Detection and Determination of Depth of Rigid Bottom in Backcalculation of Layer Moduli from Falling Weight Deflectometer Data," Transportation Research Record No. 1293, Transportation Research Board, Washington, D.C., 1991. ↩
  11. Winters (1993). "The PACCAR Pavement Test Section — Instrumentation and Validation," Master's thesis, Department of Civil Engineering, University of Washington, Seattle, Washington, March 1993. ↩
  12. (1990). "Relating Bearing Capacity to Pavement Condition," Ph.D. Dissertation, Royal Institute of Technology, Department of Highway Engineering, Stockholm, Sweden, 1990. ↩
Report an error

Reports go to the editors, who check them against the cited sources. Corrections are noted at the top of the article.

Only if you want a reply. It is not published.