U.S. Department of Transportation
Federal Highway Administration
1200 New Jersey Avenue, SE
Washington, DC 20590
202-366-4000


Skip to content
Facebook iconYouTube iconTwitter iconFlickr iconLinkedInInstagram

Federal Highway Administration Research and Technology
Coordinating, Developing, and Delivering Highway Transportation Innovations

Report
This report is an archived publication and may contain dated technical, contact, and link information
breadcrumb
Publication Number: FHWA-HRT-04-079
Date: July 2006

Seasonal Variations in The Moduli of Unbound Pavement Layers

Chapter 5: Prediction of Backcalculation Pavement Layer Moduli

INTRODUCTION

This chapter discusses the development of relationships to predict backcalculated moduli for unbound pavement layers. The models discussed address the effects of both moisture condition and stress state on the moduli of these layers in the nonfrozen condition.

The backcalculated pavement layer moduli used in this work are those from the LTPP database that remained in the working data set after the review process discussed under Evaluation of Backcalculated Layer Moduli" in Chapter 3. The stress parameters considered were computed by the author, as discussed under "STRESS PARAMETERS" in Chapter 3. Input to these computations included the backcalculated pavement layer moduli, the layer thickness, and Poisson’s ratios used in the backcalculation process, and (in the case of overburden stresses) the layer densities, daily mean layer moisture contents, and water table depth data from the LTPP database. The moisture contents considered in this work were daily mean values computed from the TDR-based moisture content data for each pavement layer (see "Moisture Parameters," Chapter 3). The sources of variation reflected in the data set included variations with FWD load (by virtue of the use of multiple FWD load levels at each test date and time), point-to point variations, within-day variations (by virtue of repeated test cycles within a day), and longer-term variations (through use of data for multiple test dates).

EXPLORATION OF RELATIONSHIPS BETWEEN BACKCALCULATED MODULI AND EXPLANATORY VARIABLES

Attempts to develop predictive models for backcalculated pavement layer moduli were preceded by examination of a series of correlation matrices exploring the relationship between the backcalculated layer moduli, stress state variables, and the TDR-monitored moisture content.

The variables considered in these matrices are: Log(E); the mean volumetric moisture content (Vm) for the layer; the applied FWD load (P); the radial load stress (sr); the bulk and octahedral shear stresses due to load and overburden (?, and t); log(?/Pa); and log(t/Pa+1). Radial overburden was computed assuming k0 = 1.0 in all cases. The stress parameters were computed for points directly beneath the center of the loaded area, and at depths of one-quarter-depth, mid-layer, and three-quarters-depth within the finite layers, and at 0.1, 0.2, and 0.3 m (3.9, 7.9, and 11.8 inches) beneath the layer interface for semi-infinite subgrade layers. Pertinent observations drawn from these correlation matrices are as follows.

The observed correlation between modulus and moisture content is summarized in Table 32.For 53 percent of the pavement layers, the correlation between modulus and moisture content is negative, indicating that increases in moisture are associated with decreases in modulus. The opposite is true for the remaining 47 percent of the pavement layers. The strength of the observed correlation, whether positive or negative, varies tremendously-from essentially zero, to 0.73 on the positive side, and 0.99 on the negative side. Both positive and negative correlations are observed for all of the soil classes represented in the data set, and all layers in the pavement structure. The prevalence of positive correlations is contrary to expectations.

Table 33 summarizes the observed correlation between modulus and applied FWD load. As was the case for the moisture-modulus correlations, the strength of the correlations between modulus and applied FWD load varies considerably, with values spanning the range from - 0.80 to 0.92. For the majority of the base layers, the expected positive correlation, implying increasing modulus with increasing FWD load, is observed. The direction of the observed correlations is less consistent for the deeper layers. The prevalence of negative correlations (implying decreasing modulus with increasing FWD load) for granular layers is contrary to expectations, and inconsistent with the findings presented in Chapter 3 that were based on paired comparisons of moduli for different load levels with all other factors constant. For this reason, it is reasonable to attribute the unexpected frequency of negative correlations to the confounding influence of variations in moisture and other factors with time.

Table 32. Summary statistics for correlation between E and mean layer moisture
  Negative Correlation Positive Correlation
Number of
Layers
Range Number of
Layers
Range
Overall 31 -0.01 to -0.99 28 0.01 to 0.73
Layer 2
(base)
10
(A-1-a,
A-1-b)
-0.03 to -0.70 12
(A-1-a, A-1-b,
A-2-4)
0.03 to 0.58
Layer 3
(subbase or
subgrade)
15
(A-1-b, A-2-
4, A-2-6,
A-3, A-4)
-0.00 to -0.99 7(A-1-a, A-2-4, A-2-6, A-3, A-4, A-6) 0.04 to 0.52
Layer 4
(subgrade)
6
(A-2-4, A-3,
A-4)
-0.21 to -0.99 9
(A-1-b, A-2-4,
A-2-6, A-3,
A-4)
0.01 to 0.73

In reviewing the correlation matrices, it was found that the backcalculated modulus is often less strongly correlated with moisture than with one or more of the stress parameters considered. This is true for all of the base layers, 82 percent of the subbase/upper subgrade (layer 3) layers, and 80 percent of the subgrade/lower subgrade (layer 4) layers. The exceptions occurred for the upper subgrade layer at section 131005 (A-4), the subgrade at section 251002 (an A-3 soil), both the subbase and subgrade layers at section 831801 (A-2-6 and A-2-4, respectively) and both the subbase and subgrade at section 241634(both A-4). For 63 percent of the layers, the backcalculated modulus is more strongly correlated with the radial load stress computed at one or more depths than it is with either moisture, or the bulk or octahedral shear stresses.

Table 33. Summary statistics with regard to observed correlation between E and applied FWD load
Layer Negative Correlation Positive Correlation
Number of
Layers
Range Number of
Layers
Range
All 25 -0.01 to -0.80 34 0.00 to 0.92
2
(base)
3
A-1-b
-0.06 to -0.15 19
A-1-a, A-1-b,
A-2-4
0.00 to 0.92
3
(subbase or
subgrade)
9
(A-1-b, A-2-4,
A-2-6, A-3,
A-4, A-6)
-0.01 to -0.66 13
(A-1-a, A-1-b,
A-2-4, A-3,
A-4)
0.00 to 0.80
4
(subgrade)
13
(A-1-b, A-2-4,
A-2-6, A-3,
A-4)
-0.03 to -0.88 2
(A-2-4, A-3)
0.10 to 0.35

The relative strength of the observed correlations between log En and the bulk and octahedral shear stress parameters is summarized in Table 34. This information is inconsistent with expectations based on laboratory Mr test results. Lab data typically indicate that bulk stress is the more important predictor of modulus for granular materials, whereas the correlations for the backcalculated moduli indicate that the octahedral shear stress is often the more important predictor for both granular and fine-grained materials. Conversely, based on laboratory data, one would expect the octahedral shear stress to be the more important predictor for the finegrained soils, but the data in Table 34 indicate that this is not always the case, as A-4 soils are among those for which the backcalculated modulus is more strongly correlated with the bulk stress than with the octahedral shear stress.

Table 34. Relative strength of correlation between log En and bulk and octahedral shear stress parameters (number of layers and soil classes represented)
Layer Bulk Stress
Correlation Is
Stronger
Octahedral Shear
Stress Correlation Is
Stronger
Relative Strength of
Correlation Varies
with Computation
Depth
2 5
(A-1-a and A-1-b)
11
(A-1-a and A-1-b)
6
(mostly ?)
(A-1-a)
3 2
(A-2-6)
12
(A-1-b, A-2-4, A-3,
A-4)
8(mostly t)
(A-1-a, A-1-b, A-2-4,
A-3, A-6)
4 4
(A-1-b, A-2-4, A-4)
11
(A-1-b, A-2-4, A-2-6,
A-3, A-4)
0

The sign (positive or negative) of the observed correlations between backcalculated modulus and the bulk and octahedral shear stress parameters is summarized in Table 35. Based on laboratory test experience, positive correlation (corresponding to increasing modulus with increasing stress) between modulus and the bulk stress is expected for all granular layers. For the backcalculated moduli, both positive and negative correlations with bulk stress occur. The relative frequency of occurrence of positive and negative correlations varies with the layer and the stress computation depth considered. Negative correlations (indicative of stresssoftening behavior) are expected for the fine-grained soils. Again, some observed correlations for the backcalculated layer moduli are consistent with this expectation, while others are not.

The correlation matrices were reviewed to identify the depth at which the computed stress parameters were most strongly correlated with the backcalculated moduli. The findings of this review are summarized in Table 36. More often than not, the correlations for the one-quarterdepth and three-quarters-depth stress computation points are stronger than those for the midlayer stresses. The frequencies with which the strongest correlations between E and bulk stress are observed at one-quarter-depth and three-quarters-depth in the layer are similar. The strongest correlations with octahedral shear stress occur most often at three-quarters-depth within the layer for the finite layers, and at 0.1-m below the layer interface for the semiinfinite subgrade layers. Between-depth differences in the strength of the observed correlation range from zero (for t and log t at section 351112, layer 4) to more than 0.60 (for log ? at section 040113, layer 2).

Table 35. Sign of observed modulus-stress correlation
Layer
(Depth)
Number and Classification of Layers with...
Positive
Correlation for
Both ? and t
Negative
Correlation for
Both ? and t
Positive
Correlation for
? and Negative
Correlation for
t
Negative
Correlation for
? and Positive
Correlation for
t
2
(1/4-depth)
13
A-1-a and
A-1-b
0 0 9
A-1-a and
A-1-b
2
(mid-layer)
5
A-1-a
0 0 17
A-1-a and
A-1-b
2
(3/4-depth)
2
A-1-a
0 0 20
A-1-a and
A-1-b
3
(1/4-depth)
17
A-1-a, A-1-b,
A-2-4, A-3, A-4
5
A-2-4, A-2-6,
A-4, A-6
0 0
3
(mid-layer)
14
A-1-a, A-1-b,
A-2-4, A-3, A-4
5
A-2-4, A-2-6,
A-4, A-6
0 3
A-1-b, A-3
3
(3/4-depth)
10
A-1-a, A-1-b,
A-2-4, A-4
5
A-2-4, A-2-6,
A-4, A-6
0 7
A-1-b, A-3, A-4
4
(1/4-depth)
3
A-2-4, A-3, A-4
10
A-1-b, A-2-4,
A-2-6, A-3,
A-4
2
A-1-b and A-4
0
4
(mid-layer)
3
A-2-4, A-3, A-4
9
A-1-b, A-2-4,
A-2-6, A-3,
A-4
3
A-1-b, A-2-4,
A-4
0
4
(3/4-depth)
3
A-2-4, A-3, A-4
9
A-1-b, A-2-4,
A-2-6, A-3, A-4
3
A-1-b, A-2-4,
A-4
0
Table 36. Distribution of strongest modulus-stress correlations with respect to stress computation depth
Stress
Parameter
Number of Layers for Which Strongest Correlation with E
Occurs at Given Depth
Finite Layers Semi-Infinite Layers
(Depth Below Interface)
1/4-
Depth
Mid-
Layer
3/4-
Depth
0.1 m 0.2 m 0.3 m
? 23 3 20 5 1 + 1tie* 7 + 1 tie
t 17 0 27 9 2 ties 3 + 1 tie
Both 15 0 15 6 1 tie 3 + 1 tie
*A "tie" occurs where the correlations for both ? and t were identical to two decimal places for two or more computation depths. If equal values occurred for only one of the two parameters, the observation was tallied with the depth for which the other stress parameter had the stronger correlation.

The relationship between variations in backcalculated modulus and moisture was also examined by looking at the ratio ?E/?Vw, where ?E and ?Vw are the change in modulus and moisture relative to a reference condition. The reference condition that was used was the earliest available modulus/moisture observation for each nominal FWD load level for a particular test section and layer. Summary values by layer type are presented in Table 37. Pertinent observations follow the table.

Table 37. Change in modulus versus change in moisture content by layer type
Layer Type Load
Level
Change in
Modulus, ?E
(kPa)
Change in
Moisture
Content, ?Vw
(%)
?E/?Vw
(kPa/%)
n
Min. Max. Min. Max. Min. Mean Max. Std.
Dev.
Base 1 -290 294 -5.9 11.4 -4590 27 5820 698 430
2 -362 428 -6.8 11.4 -1992 11 3240 295 597
3 -334 516 -4.9 11.4 -3000 6 3600 347 695
4 -359 500 -5.9 11.4 -4860 -1 7860 565 763
Subbase 1 -117 107 -13.3 5.9 -212 -11 165 41 96
2 -101 152 -12.7 5.9 -160 -7 184 49 131
3 -126 161 -8.3 5.9 -117 -7 148 39 157
4 -113 139 -6.5 9.9 -104 -7 176 37 152
Finite
Subgrade
1 -268 113 -19.4 17.4 -1055 7 2728 254 428
2 -209 157 -19.4 10.0 -2154 -16 1310 202 585
3 -230 191 -19.4 10.0 -2051 -18 1343 187 700
4 -225 121 -19.4 10.0 -1553 -7 1527 144 801
Semi-
infinite
Subgrade
1 -235 359 -6.5 20.8 -420 5 960 111 287
2 -218 347 -16.6 20.8 -3520 -23 1920 219 378
3 -106 283 -9.4 20.8 -4240 -29 800 255 448
4 -241 275 -14.7 20.8 -884 8 3520 247 492

The ratio ?E/?Vw is highly variable, with a range of -4,590 to 7,860 kPa (-666 to 1140 psi) per 1-percent change in moisture content. Both positive and negative values of the ratio are observed for 92 percent of the individual layer and load level combinations. Thus, in most cases, both increases and decreases in backcalculated modulus are associated with increases in mean layer moisture content. Exceptions to this general rule occur on a consistent basis (across all FWD load levels) for the subgrade layer at section 041024 (Arizona), the subgrade layer at section 081053 (Colorado), and the base layer at section 906405 (Saskatchewan). However, in all of these cases, the number of observations considered is small, so it is not clear that these layers are truly exceptional in this regard.

The mean values of the ratio ?E/?Vw for individual layers and load levels vary in the range of -403 to 1,089 kPa (-58 to 158 psi) per one percent change in moisture content, with the overall mean value being -4. The coefficient of variation (COV) observed for individual combinations of layer and load level varies from a low of 13 percent for the base (layer 2) at section 041024 (Arizona), to a high of more than 20,000 percent for the upper subgrade (layer 3) at section 231026 (Maine), with the overall average COV being 761 percent. In many cases, both the mean ratio and the extent of variation therein vary between load levels, but universally applicable trends are not evident.

In looking at the statistics on a layer-type basis, it may be observed that the mean ratio for the base layers decreases with increasing load on a consistent basis. However, the standard deviations associated with the mean values are so large in comparison to the means that this trend in not particularly meaningful. The COV for the base layer group increases with increasing load, from a low of 2,605 percent for the nominal 27-kN (6,070 lbf) load to a high of 111,707 percent for the 71-kN (15,961 lbf) load level. In contrast, the subbase layer values for the 40-, 53-, and 71-kN (8,992-, 11,915-, and 15,961-lbf) load levels are remarkably similar. Overall, the high degree of variability in the ratio ?E/?Vw suggests that factors other than moisture (e.g., stress conditions) may be important to explaining the observed variations.

In summary, plots of modulus versus moisture, correlation matrices, and the ratio ?E/?Vw were examined with the following findings.

  • Negative correlations between the bulk stress and E are a frequent occurrence.
  • E is often more strongly correlated with sr than with other variables.
  • The correlation between E and Vw is often relatively weak.
  • The ratio ?E/?Vw is highly variable.

Together, these observations suggest that moisture is not always the primary driver of seasonal variations in backcalculated pavement layer moduli.

MODELS OF THE FORM E/Pa=k1?/Pak2(t/Pa+1)k3

Initial efforts to develop predictive models for backcalculated pavement layer moduli focused on the model form presented as Equation 26,

Elastic modulus divided by atmospheric pressure is equal to regression constant 1 times bulk stress divided by open parenthesis atmospheric pressure close parenthesis to the regression constant 2 times open first parenthesis open second parenthesis shear stress divided by atmospheric pressure close second parenthesis plus 1 close first parenthesis to the regression constant 3 (26)

with k1 taken to be k0*10c1Vw.(Vw is the volumetric moisture content, expressed as a percent;ci and ki are regression constants; and all other variables are as previously defined.) This model form was selected to be consistent with the constitutive model to be used for unbound materials in the 2002 Guide for Design of New and Rehabilitated Pavement Structures that is currently under development. Note that it is a variant of Equation 11, in which k6 is taken to be zero, and k7 is assumed to be 1. This model will henceforth be referred to as Model 1.

Both load-induced and overburden stresses were included in the computation of the stress parameters. The horizontal overburden stresses were computed using k0 values of 0.7 and 0.5 for the base and subbase/subgrade layers, respectively. Temporal variations in overburden due to variations in mean layer moisture content were considered in the calculation (although it is likely that the impact of this is very small). The stresses considered were those computed at a point directly beneath the center of the loaded area at 0.75-depth within each finite layer, and at 0.3 m below the layer interface for each semi-infinite layer.

For the purposes of the regression modeling, Model 1 was transformed by taking the log of both sides of Equation 26, so that standard multiple regression could be used to derive the regression coefficients. The transformed model (with the moisture term explicitly included) is given in Equation 27.

Logarithm of elastic modulus divided by atmospheric pressure is equal to logarithm of regression constant K of 0 plus regression coefficient C of 1 times volumetric moisture content plus regression constant K sub 2 times log of (open first parenthesis) bulk stress divided by atmospheric pressure (close first parenthesis) plus regression constant K of 3 times log of (open second parenthesis) shear stress divided by atmospheric pressure plus 1 (close second parenthesis). (27)

The multiple regression results obtained for individual pavement layers are presented in Table 73 in Appendix E: Multiple Regression Results for Individual Pavement Layers. From a standpoint of goodness of fit, as judged by the fraction of variance explained by the model (R2) and the standard error ratio (Se/Sy) computed for log E, as well as the bias and standard error ratio for E, many (but not all) models are quite good. However, although most of the materials in question are granular, the majority of the k2 values are negative, indicating that the modulus decreases with increasing bulk stress (i.e., "stress softening"). This is inconsistent with the stress-hardening behavior typically exhibited by such materials in the laboratory and the load-hardening behavior of the pavement layers discussed previously. Most positive values for k3 are similarly inconsistent with laboratory test findings. The signs of both k2 and k3, however, are consistent with correlations discussed previously.

Several alternative approaches to computation of the stress parameters were explored in an attempt to produce models with regression coefficients that are more consistent with laboratory test results. Those alternatives are summarized in Table 38. The factors that were varied were the assumed values for K0 used in computing the horizontal overburden stresses, and whether overburden stresses were or were not included in the computation of the stress parameters. The K0 values were varied to explore the possibility that the poor agreement between the regression results obtained with backcalculated layer moduli and laboratory test results came about as a result of incorrect assumed values for K0 (the true values being unknown). The use of only the load-induced stresses for one or both of the stress parameters was pursued to explore the possibility that overburden stresses are not important to the observed variations in backcalculated pavement layer moduli. Stresses computed at threequarters-depth in the pavement layer were used in this analysis, and, k1 was taken to be constant (i.e., not a function of the volumetric moisture content). Thus, the transformed model form used in the multiple regression to obtain the model coefficients is:

Logarithm of (open first parenthesis) elastic modulus divided by atmospheric pressure (close first parenthesis) is equal to log of regression constant 1 plus regression constant 2 times log of (open second parenthesis) bulk stress divided by pressure (close second parenthesis) plus regression constant 3 times log of (open third parenthesis) shear stress divided by atmospheric pressure plus 1 (close third parenthesis). (28)

The results obtained for the alternative stress computation assumptions are presented in Table 74 through Table 80 in Appendix E: Multiple Regression Results for Individual Pavement Layers. The regression coefficients for these models are summarized by soil class in Table 39. Only those regression results meeting a minimum goodness of fit standard of R2 =0.60 (in log space) were included in the preparation of this summary. A complete evaluation of the adequacy of these regression models was not undertaken, as the goal of achieving coefficients consistent with laboratory test results was not achieved.

Table 38. Combinations of K0 and stress components considered in regression modeling
Set k0 Bulk Stress Octahedral
Shear Stress
Base Subbase/Subgrade
1 0.7 0.5 Load +
Overburden
Load +
Overburden
2 0.7 0.5 Load +
Overburden
Load Only
3 1.0 1.0 Load +
Overburden
Load +
Overburden
4 2.0 1.0 Load +
Overburden
Load +
Overburden
5 2.0 1.0 Load +
Overburden
Load Only
6 Not Pertinent Not Pertinent Load Only Load Only
7 0.7 0.5 Load Only Load +
Overburden

The results obtained with all stress data sets identified in Table 38 are qualitatively similar to those obtained in the initial effort. The majority of the k2 values are negative for both the granular and fine-grained materials, and the k3 values are typically positive. While the k1 values are similar in magnitude to the laboratory resilient modulus coefficients reported by Von Quintus and Killingsworth (see Table 9), the k2 and k3 values differ in both magnitude and sign.[5] For example, Von Quintus and Killingsworth report k1 values in the range of 250 to 2,323; k2 values in the range of -0.18 to 1.07; and k3 values in the range of -0.33 to 0.61 for base materials, while the corresponding Table 39 values for A-1-a soils vary in the ranges of 59 to 2078; -9.21 to -0.07; and -0.41 to 12.8, respectively. None of the options considered resolves the inconsistency with laboratory test results. Subsequent efforts to identify an alternative radial location for stress computation that would yield acceptable results from a standpoint of both goodness of fit and consistency with laboratory test results were similarly unsuccessful.

Table 39. Regression coefficients for alternative stress computations
  Set 1 Set 2 Set 3 Set 4 Set 5 Set 6 Set 7
k1 k2 k3 k1 k2 k3 k1 k2 k3 k1 k2 k3 k1 k2 k3 k1 k2 k3 k1 k2 k3
A-1-a
Mean 1019 -0.84 2.58 1055 -0.85 2.44 1139 -1.14 2.91 1330 -0.67 1.87 1160 -0.67 1.75 1047 -0.26 2.72 972 -0.26 2.60
Std. Dev. 449 1.60 2.95 418 1.46 2.91 409 2.34 3.00 414 0.36 1.06 495 0.37 1.04 416 0.13 2.77 419 0.12 2.12
Max 1907 -0.09 2.58 1935 -0.85 12.5 1962 -0.10 12.8 2078 -0.10 3.99 2002 -0.09 3.99 1887 -0.08 12.2 1791 -0.07 9.03
Min 59 -6.81 12.6 170 -0.09 -0.39 668 -9.21 0.83 734 -1.47 1.87 94 -1.54 -0.41 431 -0.47 0.67 247 -0.43 0.56
No.Layers 16 16 14 14 14 16 15
A-1-b
Mean 730 -0.69 3.53 913 -0.57 3.00 947 -0.91 3.56 1027 -0.72 2.93 989 -0.72 2.96 1376 -0.25 3.28 890 -0.21 2.83
Std. Dev. 536 0.70 2.64 554 0.58 2.19 500 1.01 2.64 572 0.73 2.20 524 0.73 2.19 1481 0.19 3.16 482 0.14 1.87
Max 1592 -0.13 6.37 1647 -0.13 6.21 1685 -0.06 7.79 1810 -0.12 6.52 1694 -0.12 6.52 4580 0.01 8.16 1561 -0.09 5.52
Min 232 -1.81 0.66 295 -1.69 0.63 298 -2.69 0.68 298 -2.23 0.64 298 2.23 0.64 307 -0.50 -1.35 280 0.47 1.16
No.Layers 6 7 8 7 7 7 6
A-2-4
Mean 895 -1.31 6.80 817 -1.67 8.17 842 -2.43 8.44 912 -2.42 8.42 699 -1.97 7.82 1436 -0.19 2.03 1423 -0.18 2.06
Std. Dev. 770 1.74 6.07 1045 2.15 8.74 961 3.11 9.09 1059 3.11 9.12 769 2.33 6.52
Max 1534 -0.15 13.6 1556 -0.15 14.4 1521 -0.23 14.9 1661 -0.23 14.9 1580 -0.22 14.87
Min 40 -3.31 2.00 78 -3.19 1.99 163 -4.62 2.01 163 -4.62 1.97 163 -4.62 2.01
No. Layers 3 2 2 2 3 1 1
A-3
Mean 73 -3.40 13.5 127 -3.39 12.9 236 -4.49 12.7 236 -4.49 12.7 236 -4.49 12.74 165 -1.04 13.58 110 -0.99 13.22
No. Layers 1 1 1 1 1 1 1
A-4
Mean 152 -1.21 9.70 192 -1.12 9.23 201 -1.30 9.43 201 -1.30 9.43 201 -1.30 9.43 290 -0.57 5.97 262 -0.59 6.06
Std. Dev.                         163 0.34 4.47 168 0.39 4.80
Max                         406 -0.32 9.13 381 -0.31 9.46
Min                         175 -0.81 2.81 143 -0.86 2.67
No.Layers 1 1 1 1 1 2 2

FACTORS CONTRIBUTING TO DIFFERENCES BETWEEN LABORATORY AND FIELD-BASED CONSTITUTIVE MODEL COEFFICIENTS

In comparing backcalculated pavement layer moduli with laboratory resilient moduli, it is important to keep in mind what is being compared. As discussed previously, the laboratory test result is appropriately termed a material property-a measurable characteristic of the material. In contrast, the backcalculated layer modulus is more correctly thought of as a parameter-an estimate of the "average" material characteristic for the pavement layer, which does not fully account for spatial variations in stress state (and thus, stiffness). Furthermore, the interpretation process used to obtain the estimate is based on theory that approximates, but does not match, reality. Key theoretical assumptions that are violated to one degree or another are the assumptions of homogeneity, isotropy, and linearity.

In the interpretation of laboratory resilient modulus test data we have measurements of the confining pressure and the cyclic axial load from which the major and minor principal stresses applied to the test sample are computed. For the most common laboratory situation, the radial and vertical stresses are: (1) always compressive; (2) independent of each other; (3) independent of the modulus of the material; and (4) applicable to the sample as a whole. In contrast, when looking at backcalculated layer moduli, one must compute the applicable overburden and load-induced stresses. By virtue of the theory used (linear layered-elastic theory), the computed stresses are: (1) sometimes tensile; (2) not independent of each other; (3) not independent of the modulus of the material, and (4) applicable to a particular point in the pavement layer, rather than a well-defined volume of soil. Furthermore, it is probable that the location of the point at which the computed stress parameters are most representative of thelayer as a whole varies not only from one pavement to another, but from point to point and day to day (if not hour to hour) by virtue of spatial and temporal variations in the pavement layer properties.

The stress parameters used in modeling the stress sensitivity of the backcalculated moduli are themselves a function of the backcalculated layer moduli. This has some important implications, chief among them the potential for the radial load stress to decrease as the modulus increases, with all other factors held constant. In Figure 32 and Figure 33, data for two different FWD load levels from a single FWD drop sequence are used to illustrate the impact of this possible occurrence on the relationship between backcalculated moduli and computed stress parameters.

For each set of backcalculated layer moduli, stresses were computed for the original set of backcalculated layer moduli, and for two modified data sets. In one version, the modulus for one (and only one) layer was multiplied by 1.3. In the other, the modulus for the same layer was divided by 1.3. The simulated variations (outer points in each curve) illustrate how the computed stresses vary with the modulus of one layer, when all other factors are held constant, while the actual data points (which appear as the center point in each curve) show the overall impact of FWD load level on the stresses and moduli for the layer under consideration (including load-related variations in the moduli of other layers. The stresses were computed for a point directly beneath the center of the applied load. Three different computation depths were considered for each layer, corresponding to one-quarter-depth, midlayer, and three-quarters-depth.

In Figure 32, the observed variation in computed vertical load stress is reasonably consistent with the laboratory test situation for granular materials-a direct relationship between modulus and vertical stress. However, the same cannot be said for the radial load stress, illustrated in Figure 33. For two of the three stress computation depths, there is an inverse relationship between the radial load stress computed for the original data (center point on curves) and the backcalculated modulus and FWD load. An inverse relationship between modulus and radial load stress also occurs for all three stress computation depths considered for the simulated variations in modulus. Also, note that some computed radial stresses are negative (tensile). The impact of these variations on the relationship between backcalculated modulus and the computed load-induced bulk and octahedral shear stresses is illustrated in Figures 34 and 35. In this case, the original data points all show an increase in both modulus and bulk load stress with increasing FWD load, while both direct and inverse relationships between modulus and bulk stress occur for the simulated variations. The computed octahedral shear stresses consistently increase with increasing modulus.

The trends shown in these figures do not represent universal truth. Rather, they illustrate that use of linear layered-elastic theory to backcalculate moduli and compute the associated loadinduced stresses can result in inverse relationships between backcalculated layer moduli and bulk stress, even though the expected load-hardening behavior is observed in the backcalculated moduli. This situation is compounded by the fact that the computed radial load stress may also be negative. The varying trends in the computed radial load stress, and the occurrence of negative values, manifest themselves in varying trends in the computed bulk stress. The impact of the varying radial load stress trends and tensile radial load stresses on the computed octahedral shear stress is less obvious (but no less important), by virtue of the mathematical definition of that parameter.

Figure 32. Vertical load stress versus modulus, section 040113 A-1-a base layer

The backcalculated modulus is graphed on the horizontal axis up to 400. The vertical load stress is graphed on the vertical axis from 0 to 200 kilopascals. For original data points, vertical load stress and modulus increase with increasing falling weight deflectometer (FWD) load. For simulated variations in modulus, vertical load stress increases as modulus increases in all cases. There are six approaches on the graph: one-quarter, mid-layer, and three-quarter depth for 40 kilonewtons and 71 kilonewtons. All are increasing slightly, but the depths with 71 kilonewtons are higher than the 40 kilonewtons depths

Figure 33. Radial load stress versus modulus, section 040113 A-1-a base layer

The backcalculated modulus is graphed on the horizontal axis from 0 to 400 megapascals. The radial load stress is graphed on the vertical axis from negative 50 to positive 100 kilopascals. For original data, modulus increases with increasing FWD load. Radial load stress computed at one-quarter depth increases with increasing FWD load. Radial load stress computed at mid-layer and three-quarter depth decreases with increasing load. For simulated variations in modulus, radial load stress decreases with increasing modulus in all cases. There are six approaches: one-quarter depth, mid-layer, and three-quarter depth at 40 and 71 kilonewtons. All six approaches have an inverse relationship between modulus and radial load stress.

Figure 34. Bulk load stress versus modulus, section 040113 A-1-a base layer

The backcalculated modulus is graphed on the horizontal axis from 0 to 400 megapascals. The bulk load stress is graphed on the vertical axis from 0 to 200 kilopascals. There are six approaches: one-quarter, three-quarters, and mid-layer at both 40 and 71 kilonewtons. For original data points, modulus and bulk stress both increase with increasing FWD load. For simulated variations in modulus, bulk stress increases with modulus for one-quarter depth point and 40 kilonewtons load, but decreases with increasing modulus in all other cases. All approaches are decreasing in stress as the backcalculation is increasing.

Figure 35. Load-induced octahedral shear stress versus modulus,section 040113 A-1-a base layer

The modulus is graphed on the horizontal axis from 0 to 400 megapascals. Load-induced octahedral shear stress is graphed on the vertical axis from 0 to 200 kilopascals. There are six approaches: one-quarter, mid-layer, and three-quarter depths at both 27 and 72 kilonewton loads. The modulus and computed octahedral shear stress increase with increasing load. Computed load-induced octahedral shear stress increases with increasing modulus for all simulated variations in modulus. All six approaches are increasing in stress as the modulus increases, but the approaches with 72 kilonewtons are higher than the approaches with 27 kilonewtons.

The vertical and radial stresses (sv and sr) computed for the backcalculated moduli are compared with those for two laboratory resilient modulus test protocols in Figure 36 to Figure 38. The test protocols considered are the LTPP P-46 protocol,[80] and the proposed "harmonized" protocol developed under NCHRP project 1-28A.[45] The stresses plotted in these figures were computed for points directly below the center of the loaded area, at threefourths- depth within the pavement layer, with k0 = 1.0, for the first set of backcalculated moduli for each test date and FWD load level. Plots prepared with stresses computed at onequarter- depth or mid-layer are qualitatively similar. While the stresses plotted in these figures do not represent universal truth, they are representative of the data set used in this investigation (see Table 18 and Table 19), and serve to illustrate the conditions that can exist when linear layered-elastic theory is used to backcalculate pavement layer moduli andcompute the load-induced stresses in the pavement structure.

Figure 36. Comparison of radial and vertical stress components for granular base and subbase layers

The radial stress is graphed on the horizontal axis from negative 100 to positive 150 kilopascals. The vertical stress is graphed on the vertical axis from 0 to 400 kilopascals. There are four approaches: load only, LTPP P46, harmonized, and load plus overburden. Large fractions of the computed radial load stress are negative. Both load only and load plus overburden range between radial stress of negative 100 to 25 kilopascals and a vertical stress of 0 to 200 kilopascals. Both LTPP P46 and harmonized increase in vertical stress as radial stress increases. Harmonized and LTPP are significant components for granular base and subbase layers.

The principal stresses defined in the two laboratory test protocols are similar, though those for the harmonized protocol are more extensive. While some computed stress data points for the base and subbase materials (Figure 36) fall within the bounds defined by the laboratory test protocols, many do not. Of particular significance is that a large fraction of the computed radial load stresses are negative, particularly for the base layers.

Figure 37. Comparison of radial and vertical stress components for coarse-grained subgrade layers

Graph. Comparison of radial and vertical stress components for coarse-grained subgrade layers. The radial stress is graphed on the horizontal axis from negative 50 to positive 150 kilopascals. The vertical stress is graphed on the vertical axis from 0 to 400 kilopascals. There are four approaches: LTPP P46, computed-load and overburden, harmonized, and computed-load only. The computed-load only approach has a vertical stress below 40 kilopascals. Computed-load and overburden has a vertical stress below 50 kilopascals, between radial stresses of 10 to 50 kilopascals. Both harmonized and LTPP P46 increase in vertical stress as the radial stress increases, and both are significant.

Figure 38. Comparison of radial and vertical stress components for fine-grained subgrade layers

The radial stress is graphed on the horizontal axis from negative 50 to positive 200 kilopascals. The vertical stress is graphed on the vertical axis from 0 to 200 kilopascals. There are four approaches: LTPP P46 protocol, computed-load plus overburden, harmonized, and computed-load only. Computed-load only is clustered close to the origin. Computed-load plus overburden is clustered in an ascending pattern between the vertical stress of 30 to 50 kilopascals and the radial stress of 20 to 40 kilopascals. Both harmonized and P-46 increase higher than computed.

Computed subgrade stresses illustrated in Figures 37 and 38 are more similar to those for the laboratory stresses than those for the base and subgrade layers, but many data points fall outside the envelopes defined by the laboratory test protocols. Negative computed radial stresses are again prevalent.

The corresponding bulk and octahedral shear stresses are shown in Figure 39 to Figure 41. As one would expect given the relationships observed in Figure 36 to Figure 38, the stress states computed from the backcalculated layer moduli bear little resemblance to those for either laboratory protocol. (Note that the scales used in these figures are such that the higher stress states for the harmonized protocol are not shown.) Thus, it is reasonable to conclude that the model coefficients derived for the backcalculated layer moduli differ from those obtained for laboratory resilient modulus test results because the stress states computed for the backcalculated moduli using layered-elastic theory are not directly analogous to the lab stress states. Therefore, it is impossible to derive constitutive model coefficients that are compatible with laboratory test data from the backcalculated layer moduli.

Figure 39. Comparison of bulk and octahedral shear stresses for granular base and subbase layers

The octahedral shear stress is graphed on the horizontal axis from 0 to 150 kilopascals. The bulk stress is graphed on the vertical axis from negative 300 to positive 700 kilopascals. There are four approaches: load only, LTPP P46, harmonized, and load and overburden. Several computed fall outside the test protocols. Negative computed radial stresses are prevalent. The harmonized and P-46 approaches increase in bulk stress as the shear stress increases and are significant for granular base and subbase layers.

Figure 40. Comparison of bulk and octahedral shear stresses for coarse-grained subgrade layers

The octahedral shear stress is graphed on the horizontal axis from 0 to 50 kilopascals. The bulk stress is graphed on the vertical axis from negative 50 to positive 450 kilopascals. Computed load falls outside of the test protocol. The rest of the approaches increase in bulk stress as the octahedral shear stress increases and are significant for coarse-grained subgrade layers.

Figure 41. Comparison of bulk and octahedral shear stresses for fine-grained subgrade layers

The octahedral shear stress is graphed on the horizontal axis from 0 to 50 kilopascals. The bulk stress is graphed on the vertical axis from 0 to 600 kilopascals. There are four approaches: L T P P P46 protocol, computed-load and overburden, harmonized, and computed-load only. All approaches are within the test protocol and are significant for fine-grained subgrade layers.

In light of these observations, it might be argued that the solution to the observed discrepancies is to abandon the use of linear layered-elastic theory in the interpretation of pavement deflection data. However, while nonlinear backcalculation is an essential long-term goal, abandoning backcalculation based on linear layered-elastic theory is not a practical solution at this time. Given the current state of the art, it is likely that backcalculation based on linear layered-elastic theory (with or without adaptations to consider the stress sensitivity of materials in an approximate fashion) will remain the best available technology applicable in routine practice for at least the next several years. Thus, it is appropriate to pursue models to predict the observed variations in layer moduli backcalculated using linear layered-elastic theory, even though the model coefficients obtained are not directly analogous to those derived from laboratory test results for the same or similar materials.

ALTERNATIVE MODEL FORMS

Individual Pavement Layer Models

Given the differences between the computed stresses based on linear layered elastic theory and laboratory stress states discussed in the preceding section, it is appropriate to consider possible modifications of the constitutive model form to improve the predictive capability for the backcalculated layer moduli.

The primary modification considered was the substitution of a semi-log form for the bulk stress term in the model, as shown in Equation 29 (Model 2).

Elastic modulus is equal to regression constant 1 times 10 to the open parenthesis regression constant 2 times bulk stress divided by atmospheric pressure close parenthesis times open first parenthesis open second parenthesis shear stress divided by atmospheric pressure close second parenthesis plus 1 close first parenthesis to the regression constant 3. (29)

Use of the semi-log form for the bulk stress term enables straightforward consideration of negative computed bulk stress terms. This rationale does not apply to the octahedral shear stress term, as it cannot be negative, so the log-log form was retained for the octahedral shear stress term, to maximize consistency with the original model form.

Two variations on the model form were considered: model 2A, in which k1 was taken to be a function of the volumetric moisture content, and model 2B, in which all three regression parameters were treated as being (potentially) a function of the volumetric moisture content. These variations are presented as Equations 30 and 31, respectively.

Elastic modulus is equal to regression coefficient C sub 0 times 10 to the open first parenthesis open second parenthesis regression coefficient C sub 1 times volumetric moisture content divided by 100 close second parenthesis plus open second parenthesis regression constant K sub 2 times bulk stress divided by atmospheric pressure close second parenthesis close first parenthesis times open first parenthesis open second parenthesis shear stress divided by atmospheric pressure close second parenthesis plus 1 close first parenthesis to the regression constant K of 3. (30)

Elastic modulus is equal to regression coefficient times 10 to the (open first parenthesis) regression coefficient C sub 1 times volumetric moisture content divided by 100 plus open second parenthesis) regression coefficient C sub 2 plus log of regression coefficient C sub 2 times volumetric moisture content divided by 100 (close second parenthesis) times bulk stress divided by atmospheric pressure (close first parenthesis) times (open third parenthesis) shear stress divided by atmospheric pressure plus 1 (close third parenthesis) to the (open fourth parenthesis) regression coefficient C sub 3 plus log of regression coefficient C sub 3 times volumetric moistures content divided by 100 (close fourth parenthesis). (31)

The models were fitted using multiple regression applied to the log-transformed models, as presented in Table 40. In the equations in Table 40, the c0’= log c0. The computation of the stress parameters considered in these models differed from that for Model 1 in one respect,the value of K0 used in the computation of the radial overburden stresses. For these models,K0 was taken to be 1.0 for all layers.

Table 40. Log-transformed models used in regression modeling
  Transformed Equation
2A logE =c0'+c1Vw/100+k2?/Pa+k3log(t/Pa+1)
2B log E=c0'+c1Vw/100+c2?/Pa+c2'Vw/100*?/Pa+c3log(t/Pa+1)+c3'Vw/100*log(t/Pa+1)

Regression results obtained for individual layers using 2A and 2B are presented in Table 81 through Table 82 of Appendix E: Multiple Regression Results for Individual Pavement Layers.The regression coefficients obtained for all four models are summarized in Table 41.

Table 41. Summary of regression coefficients for models 2A and 2B
Parameter 2A 2B
c0' Mean 3.479 3.740
Minimum -12.06 -38.73
Maximum 12.61 32.38
c1 Mean -1.527 -2.439
Minimum -59.78 -164.6
Maximum 3.470 237.0
c/k2 Mean -0.408 -1.048
Minimum -4.558 -28.80
Maximum 10.57 12.58
c2' Mean - 1.626
Minimum - -56.93
Maximum - 116.2
c3/k3 Mean -1.883 1.401
Minimum -160.1 -617.3
Maximum 56.53 924.0
c3' Mean - 15.17
Minimum - -5451.1
Maximum - 3256

The regression coefficients obtained for both equations appear to be rational in the context of linear layered-elastic backcalculated moduli and computed stresses for the majority of the pavement layers considered. (They would not be considered rational in the context of traditional laboratory resilient modulus test experience.) In this context, the traditional interpretation that negative k2 values are associated with stress-softening behavior is not really applicable, because the relative importance of the bulk and octahedral shear stress terms is often the reverse of what is seen in the laboratory, due to the particulars of linear layeredelastic analysis. Despite the negative k2 values, the models predict increasing modulus with increasing FWD load for the granular layers.

While the majority of the model coefficients are rational (in the context of linear layeredelastic analysis), there are several notable exceptions. For example, the c0 values for all the models for section 831801 layer 3 for both equations are very high. The goodness of fit statistics for the section 831801 layer 3 models are very good, but the sample size is very small, so the applicability of this model is limited.

The goodness of fit statistics for Models 1, 2A, and 2B are compared on a layer-by-layer basis in Table 83 of Appendix E: Multiple Regression Results for Individual Pavement Layers, and summarized in Table 42. In some instances, the number of observations considered in the regression varied between Model 1 (n1) and the other equations (n2), by virtue of the presence of negative bulk stresses in the data sets. This was the case for 11 of the 59 layers considered. In a few instances, the number of observations excluded due to negative bulk stresses approached the total number of observations available. In two cases, the exclusion of negative bulk stress observations made it impossible to complete the regression for Model 1, due to small sample sizes.

Table 42. Summary of goodness-of-fit statistics for Models 1, 2A and 2B
    1 2A 2B
Se/Sy Mean 0.58 0.55 0.52
Minimum 0.05 0.01 0.12
Maximum 1.01 1.00 1.02
Bias Mean 1 0 0
Minimum 0 0 0
Maximum 5 5 4

In most cases, providing for the possibility that all three k-coefficients vary with moisture content yields a better fit of the data (most often reflected in a reduced standard error ratio (Se/Sy)) than the assumption that only k1 is a function of moisture content. In some cases (e.g., layer 4 for section 040114), the difference is negligible (zero percent change in bias, and standard error ratio changed by less than 0.01), and in others it is more substantial (e.g., layer 4 for section 251002, where the standard error ratio was reduced from 0.96 to 0.38).

The regression results were characterized as fully acceptable if the standard error ratio, Se/Sy, computed for Ep (as opposed to log(Ep)) was less than 0.5, and the absolute value of the bias was less than 2 percent. The results were characterized as unacceptable if the standard error ratio was greater than 0.85, or the absolute bias was greater than 5 percent. Results meeting neither set of criteria were characterized as marginal. Based on these criteria, all model forms considered yielded results that were at least marginally acceptable (from a standpoint of goodness of fit) for 61 percent of the pavement layers. Of 59 pavement layers considered, the regression results were unacceptable for 2, both subgrade layers, and 17 for which all model forms yielded fully acceptable results. The regression modeling was most successful for the base layers, and least successful for the semi-infinite subgrade layers. This discrepancy in the relative rate of success is consistent with the assumption that the true modulus of some of the semi-infinite subgrade layers is relatively constant with time and FWD load, such that the observed variation is primarily due to random error in the deflection testing and backcalculation process.

Model 2 is better suited to modeling the variations in pavement layer moduli backcalculated using linear layered-elastic theory. The rationale underlying this judgment is as follows:

  1. More often than not (for 75 percent of the layers considered), the standard error ratio for one or both of Models 2A and 2B is less than or equal to that for Model 1. This observation also holds if one compares only Models 1 and 2A, though the difference between models 1 and 2A is small enough that it may be attributable to the different assumed values of K0.
  2. Models 1 and 2 are essentially equal from a standpoint of bias.
  3. Model 2 is applicable to negative computed bulk stresses, whereas Model 1 is not.That layered-elastic theory yields tensile stresses in some situations is theoretically problematic, from a practical standpoint. However, consideration of negative (tensile)stresses when working with pavement layer moduli backcalculated using linear layered elastic theory is necessary because they occur quite frequently. Exclusion of affected data points (where they occur) reduces the size of the data set on which to base the model and may result in a biased model.

Furthermore, for 58 percent of the pavement layers, the standard error ratio for Model 2B is somewhat better than for Model 2A; consideration of the potential for moisture to affect all three k coefficients thus is recommended. It must be recognized that this recommendation applies only to layer moduli backcalculated using linear layered-elastic theory. No assessment of the applicability of the model form to laboratory resilient modulus test data has been undertaken.

Soil Class Models

The regression models discussed in the preceding sections of this chapter were each derived from the data from a single pavement layer at a particular test section. In the interest of obtaining more broadly applicable models, the working data set of moduli, stress, and moisture parameters was subdivided by soil class and layer extent (i.e., whether the layer was finite or semi-infinite in depth). Multiple regression was used to fit the data for each subset to Model 2B, and several variations on Model 2 including material parameters, in addition to stress parameters and moisture content. The model forms considered for the coefficients k1, k2, and k3 are presented in Table 44. The materials data used in the models are presented in Appendix D:Materials Data Used In Development of E Predictive Models. Density and plasticity index are presented in Table 70, while gradation information is provided in Table 71. The resulting models are presented in Table 43.

It is notable that the predictive capability of the soil class models (with or without material parameters) is often (but not always) better for the semi-infinite subgrade layers than it is for the soil classes most often found in the base and subbase layers-the opposite of what was observed for the individual layer models. Factors that may contribute to the observed discrepancy are as follows.

  1. Differences in the extent to which the layer moduli vary. As noted in Chapter 4,variations in the backcalculated moduli for the base layers tends to be greater than those for the deeper layers.
  2. Differences in the extent to which the within-section variations in backcalculated moduli are attributable to differences in stress and moisture, versus other factors, including random error in the deflection measurement and backcalculation process. The regression results for the individual pavement layers indicate that variations in stress and moisture typically explain a relatively large fraction of the observed withinsection variation for the upper pavement layers, but a lesser fraction of the observed variation for the deeper layers.
  3. Differences in the extent to which between-section variations in backcalculated moduli are attributable to quantifiable material characteristics as opposed to systematic errors arising from the backcalculation process. For example, as discussed in Chapter 2, the backcalculated moduli for the upper layers of the pavement are more sensitive to small errors in layer thickness than those for deeper layers. A systematic error in backcalculated moduli arising from layer thickness errors will have little or no effect on the within-section variation, but could have a very large impact on betweensection variation.
  4. Differences in the extent to which the data for the different layers are influenced by the limitations of linear layer-elastic theory, including the prevalence of tensile radial load stresses.

Differences in the number of sections (and therefore, the array of soils) considered in the model development may also contribute to the observed differences.

The issue of the number of soils represented in the data sets on which these models are based is an important one, as the applicability of models based on data for only a few soils is limited; application of the models to the soil class in general thus is not appropriate.

Table 43. Selected soil class models for prediction of pavement layer moduli using linear layered-elastic theory
Finite Layers
A-1-a E = 10(4.270 + 0.004PR10/P200-0.490?d+0.030VwP200)*10(-0.250-1.210Vw+0.140Vw*P200)*?/Pa*(t/Pa+1)(1.490+28.56Vw-3.330Vw*P200)
A-1-b E = 10(3.671 -0.151PR10/P200-0.163?d+0.092VwP200)*10(0.217-3.343Vw+0.260Vw*P200)*?/Pa)*(t/Pa+1)(2.786+9.388Vw-0.717Vw*P200)
A-2-4 E = 103.241(1+Vw) -1.193RD-8.083M-0.249*10-0.049(1-P200)?/Pa(t/Pa+1)13.14*P200
A-2-6 E = 10(3.003-0.102PI+0.428VwP200*10(-0.010-0.546VwP200)?/Pa(t/Pa+1)(0.720+4.774VwP200)
A-3 E = 10(5.160 + 0.090PR10/P200-0.160?d-1.350VwP200)*10(-2.320-21.92Vw+5.620Vw*P200)*?/Pa*(t/Pa+1)(7.100 +374.0Vw-75.75Vw*P200)
A-4 E = 10(6.825+ -1.840PR10/P200-2.012?d+0.009VwP200)*10(0.258+0.220Vw-0.043Vw*P200)*?/Pa*(t/Pa+1)(-0.888-42.03Vw+0.863Vw*P200)
A-6 E = 10(4.070+0.560Vw)*10(-0.990-0.040Vw)*?/Pa *(t/Pa+1)(-0.700+7.260Vw)
Semi-Infinite Subgrade Layers
A-1-b E = 10(35.64 -0.190PR10/P200-15.83?d)*10(0.820Vw -0.110VwP200)*?/Pa*(t/Pa+1)(-0.980VwP200)
A-2-4 E = 10(-26.92 - 2.580PR10/P200+16.16?d+0.180VwP200)*10(2.470+13.53Vw-1.180VwP200)*?/Pa*(t/Pa+1)(-25.69 -101.0Vw+10.55VwP200)
A-3 E = 10(15.09 - 38.83PR10/P200-1.438VwP200)*10(-1.185+2.251VwP200)*?/Pa*(t/Pa+1)(1.263 -36.65Vw+1.868VwP200)
A-4 E = 10(2.850 +1.510Vw)*10(1.020-2.85Vw)*?/Pa*(t/Pa+1)(-9.380+25.44Vw)
Where: PR10/P200 = Percent retained on the No. 10 sieve/percent passing 200
P200 = Percent passing 200
Vw = Volumetric moisture content (decimal fraction)
?d = In situ dry density
RD = Relative density (?d/?optimum)
Dn= Sieve opening at which n percent of soil passes, n = 10, 20, 30,…
M = Material constant (adapted from Yang[34])
PI = Plasticity index
equation43
Table 44. Model forms considered for k1, k2, and k3 in Model 2
Model k1 k2 k3
2C 10(c0+c1PR10/P200+c1’?d+c1’’VwP200) c2 + c2’Vw+c2’’VwP200 c3+c3’Vw+c3’’VwP200
2D (coarse) 10(c0+c1D60+c1RD+c1’’VwP200) c2 + c2’Vw+c2’’VwP200 c3+c3’Vw+c3’’Vw*P200
2D (fine) 10(c0+c1PI+c1RD+c1’’’VwP200) c2 + c2’Vw+c2’’VwP200 c3+c3’Vw+c3’’Vw*P200
2E 10(c0+c1Vw+c1’P10) c2+c2’Vw+c2’’P10 c3+c3’Vw+c3’’P10+c3’’’RD
2F 10c0(1+Vw)c1RDc1’’Mc1’’’ c2(1-P200) c3P200
2G 10c0(1+Vw)c1RDc1’’Mc1’’ c2 + c2’Vw + c2’’(1-P200) c3+c3’ P200+c3’’Vw
Where: PR10/P200 = Percent retained on the No. 10 sieve/percent passing 200
P200 = Percent passing 200
Vw = Volumetric moisture content (decimal fraction)
?d = In situ dry density
RD = Relative density (?d/?optimum)
Dn= Sieve opening at which n percent of soil passes, n = 10, 20,30,...
M = Material constant (adapted from Yang[34])
PI = Plasticity index
equation

Development of more general models-i.e., models applicable to all coarse-grained or all fine-grained soils-was attempted, but was not successful.

SUMMARY

Regression models to predict backcalculated pavement layer modulus as a function of stress state, moisture content, and (for the soil class models) material parameters were developed. In developing these models, regression coefficients derived for layer moduli backcalculated using linear layered elastic theory were found to be inconsistent with those obtained for laboratory test results for the same constitutive model form and similar materials. In particular, the sign of the k2 parameter for the backcalculated moduli for granular materials was consistently negative, when positive values are expected. The reasons for this discrepancy were explored, with the finding that the stress states and the relationships between the stress parameters computed for the backcalculated moduli and those typically used in laboratory testing are dissimilar. Key factors contributing to the dissimilarity include the lack of independence between the backcalculated moduli and the (computed) stress parameters used in the model development, and the potential for tensile computed stresses.

Development of empirical models to predict layer moduli derived using linear layered-elastic theory was pursued to facilitate consideration of seasonal variations pending the development of improved backcalculation procedures.

Model 2, which is restated below, is recommended for use in modeling the stress sensitivity of layer moduli backcalculated using linear layered elastic theory. It is further recommended that Model 2B be used to consider the combined effects of stress and moisture on linear elastic backcalculated pavement layer moduli.

Elastic modulus is equal to regression constant 1 times 10 to the regression constant 2 times (open first parenthesis) bulk stress divided by atmospheric pressure (close first parenthesis) times (open second parenthesis) shear stress divided by atmospheric pressure plus 1 (close second parenthesis) to the (open third parenthesis) regression coefficient C sub 3 plus log of regression coefficient C sub 3 times volumetric moistures content divided by 100 (close third parenthesis). (32)

Elastic modulus is equal to regression coefficient times 10 to the open bracket open parenthesis regression coefficient C sub 1 times volumetric moisture content divided by 100 close parenthesis plus open parenthesis regression coefficient C sub 2 plus log of regression coefficient C sub 2 times volumetric moisture content divided by 100 close parenthesis times open parenthesis bulk stress divided by atmospheric pressure close parenthesis close bracket times open first parenthesis open second parenthesis shear stress divided by atmospheric pressure close second parenthesis plus 1 close first parenthesis to the open bracket regression coefficient C sub 3 plus open parenthesis log of regression coefficient C sub 3 times volumetric moistures content divided by 100 close parenthesis close bracket. (33)

Previous | Table of Contents | Next

 

Federal Highway Administration | 1200 New Jersey Avenue, SE | Washington, DC 20590 | 202-366-4000
Turner-Fairbank Highway Research Center | 6300 Georgetown Pike | McLean, VA | 22101