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
Back to Publication List        
Publication Number:  FHWA-HRT-15-074     Date:  September 2016
Publication Number: FHWA-HRT-15-074
Date: September 2016

 

Pavement Structural Evaluation at the Network Level: Final Report

 

CHAPTER 7. 3D-MOVE CALIBRATION

7.1 Introduction

This chapter presents the results of the 3D-Move software calibration using data collected from installed sensors at the MnROAD facility during the September 2013 TSDD field trials. The objective of this calibration was to enable the use of the 3D-Move software in the development of methodologies and/or procedures for incorporating TSDD measurements into network-level PMS applications. This software was chosen because it is ideally suited to evaluate and compare pavement responses measured with TSDDs. It estimates the dynamic pavement responses at any given point within the pavement structure using a continuum-based finite layer approach.

The dynamic pavement surface deflections measured by TSDDs are affected by many factors, which include pavement layer characteristics (e.g., thickness, stiffness, and damping properties), vehicle loading (e.g., tire configuration and contact stress distribution), and vehicle speed. Since rate-dependent (viscoelastic) material properties can be accommodated by 3D-Move, it is an ideal tool to study pavement response as a function of vehicle speed through the direct use of the frequency sweep test data (dynamic modulus and damping) of AC mixture.

The applicability and versatility of the 3D-Move approach has been reported in literature.(35,31) Measured pavement responses with the existing and newly installed geophones and accelerometers during MnROAD trials provided a great opportunity to validate the applicability of 3D-Move due to loading from TSDDs. This chapter focuses on details relating to the 3D-Move calibration effort.

Purpose and Overview of Calibration Effort

One of the major tasks of this project was to simulate pavement surface deflections using numerical models with a focus on understanding the parameters that affect the TSDD measurements. Those parameters include changes in TSDD vehicle speed, pavement layer properties (e.g., age and moisture), and vehicle loading (e.g., tire configuration, load, and inflation pressure). To accomplish this, calibration of the analytical tool that simulates TSDDs was required (i.e., it is critical to verify that the analytical tool and the corresponding software are used correctly to assess pavement responses and that they provide reliable results). Toward this end, calibration of 3D-Move relative to TSDDs was undertaken first, and the calibrated software was subsequently used in the simulations of pavement responses.

As noted earlier in this report, during the September 2013 TSDD field trials, data were collected from the embedded geophones and accelerometers (which were installed as a part of this project) as well as also by the two TSDDs. The measured velocity time histories from the geophones along with the computed displacement time histories (referred to as “project measurements”) are important since they are the basic measurements made by the TSDDs and subsequently used for pavement condition evaluations. In addition, the MnROAD facility had existing sensors (subsequently referred as MnROAD sensors) capable of measuring earth pressures and pavement strain responses under vehicle loading (i.e., dynamic). Those MnROAD sensor data were included in the project database (see section 5.7). The resulting data were used to create a comprehensive database of independently measured pavement responses (i.e., stresses, strains, and displacements) from different sources, including MnROAD sensors, project sensors, and TSDDs.

The material properties that are appropriate for the conditions (e.g., temperature and moisture) that existed at the time of TSDD trials are critical inputs to 3D-Move. Some pavement layer properties are sensitive to temperature changes, while others are affected by moisture changes. Consequently, some of the TSDD trials were conducted at different times of the day (e.g., early morning and late afternoon). A study of the effects of moisture changes was not possible within the project schedule. Data from FWD testing that was conducted prior to the TSDD field trials were available. An iterative procedure that involved selecting numerous case scenarios where 3D-Move computed responses for a wide range of inputs were compared with those actually measured at the MnROAD facility. The main variables relative to the 3D-Move inputs were the pavement layer configuration and material properties. These material properties must be reasonable and consistent with available data (e.g., FWD measurements). It is important that the appropriate case scenarios should bracket the measured responses in a consistent manner covering all independent measurements (i.e., stresses, strains and displacements).

7.2 3D-Move Program Overview and Calibration

The 3D-Move program was developed at UNR to evaluate pavement response. The program can account for important pavement response factors such as the moving traffic-induced complex 3D contact stress distributions (normal and shear) of any shape, vehicle speed, and viscoelastic material characterization for the pavement layers. The finite-layer approach used in the formulation of the program treats each pavement layer as a continuum and uses the Fourier transform technique; therefore, it can handle complex surface loadings such as multiple loads and non-uniform tire pavement contact stress distribution. This approach is suitable to analyze tire imprints of any shape, including those generated by wide-base tires.(32,33,2) The finite-layer method is much more computationally efficient than the moving load models based on the finite element method.(51,52) This is because often times, the pavements are horizontally layered, and pavement responses are required only at a few selected locations. For such problems, the finite layer approach of 3D-Move is ideally suited. The frequency domain solutions along with complex modulus formulation adopted in 3D-Move enables a direct use of the frequency sweep test data of AC mixture in the analysis.

A number of field calibrations (e.g., Penn State University test track, MnROAD, and UNR Off-Road Vehicle study) that compared a variety of independently-measured pavement responses (i.e., stresses, strains, and displacements) with those computed by 3D-Move have been reported in the literature.(35,31) ViscoRoute, which is a semi-analytical dynamic multi-layer model, is currently the only other available software capable of handling loaded areas with uniform contact vertical pressure.(10) It is an ideal program that can be compared with 3D-Move in terms of incorporating viscoelastic properties in its formulation. Hajj et al. compared computed maximum transverse strains by 3D-Move and ViscoRoute at the bottom of thin and thick AC layers.(36) The computed responses were found to be within 6 percent in the 60 cases considered. These verification studies demonstrate the validity, applicability, and versatility of the 3D-Move approach. Since details of 3D-Move are available in the literature, they are not addressed in this report. (See references 35, 31–33, and 2.) However, details concerning the application of 3D-Move to this project (e.g., inputs and response outputs) are presented next.

7.3 Materials Characterization of MnROAD Accuracy Cells

MnROAD Accuracy Cells

Cells 3 and 19 in the MnROAD mainline and cell 34 in the MnROAD LVR were selected to evaluate the accuracy of the TSDDs and hence are referred to as the “accuracy cells.” All three cells consisted of AC pavements. As described in chapter 5, each cell was instrumented with four geophones and one accelerometer during the project. These cells were also previously instrumented during construction by MnROAD with sensors such as strain gauges, pressure cells, LVDTs, moisture gauges, TCs, and tipping buckets. In the calibration process, only load-induced measured response values of strains, stresses, and displacements from the referenced sensors at the three accuracy cells were compared with those computed with 3D-Move.

Material Characterization

The characterization of the material properties for the three accuracy cells is an important step in the calibration process.

The required 3D-Move inputs for the surface layer (AC considered viscoelastic) were as follows:

  1. Thickness.

  2. Unit weight.

  3. Dynamic modulus and damping as a function of temperature and frequency.

  4. Poisson's ratio.

The required 3D-Move inputs for unbound layers (base and subgrade, which were considered LE) were as follows:

  1. Thickness.

  2. Unit weight.

  3. Resilient modulus.

  4. Poisson's ratio.

MnROAD maintains a database containing laboratory and field-testing results on soils, aggregates, asphalt mixtures, asphalt binders, concrete mixtures, and other materials. The database also contains cell-specific information, including layer thickness, type of layers, and cross sections of cells at time of construction and subsequent treatments. Figure 124 shows the pavement structure for the three MnROAD accuracy cells. As shown in table 29, FWD testing was done a few days prior to the TSDD field trials, which facilitated the characterization of existing material properties for the pavement layers in the accuracy cells.

Figure 124. Illustration. Pavement structure of MnROAD accuracy cells. This illustration shows the pavement structure for the three MnROAD accuracy cells. The accuracy cells from left to right are cell 34 in the low-volume road (LVR) and cells 19 and 3 in the mainline. Cell 34 has 4 inches (101.6 mm) of asphalt concrete (AC), 12 inches (304.8 mm) of class 6 clay as the unbounded aggregate base, and clay as the subgrade. Cell 19 has 5 inches (127 mm) of AC, a 31-inch (787.4-mm) base which includes 12 inches (304.8 mm) of class 5 clay, 12 inches (304.8 mm) of class 3 clay, and 7 inches (177.8 mm) of granular blend, followed by clay as the subgrade. Cell 3 has 3 inches (76.2 mm) of AC, a 43-inch (1,092.2-mm) base which includes 6 inches (152.4) of full-depth reclamation (FDR) with engineered emulsion, 2 inches (50.8 mm) of FDR, 2 inches (50.8 mm) of class 5 clay, and 33 inches (838.2 mm) of class 3 clay, followed by clay as the subgrade.

1 inch = 25.4 mm

Figure 124. Illustration. Pavement structure of MnROAD accuracy cells.

Table 29. Date and time of FWD and TSDD testing.
Device Accuracy Cell Date and Time
Cell 3 Cell 19 Cell 34
FWD 9/16/2013 at 2 p.m. 9/20/2013 at 2 p.m. 9/5/2013 at 3 p.m.
TSD 9/26/2013 at 12:30 p.m. 9/26/2013 at 10 a.m. 9/24/2013 at 3 p.m.
RWD 9/23/2013 at 2:30 p.m. 9/23/2013 at 9 a.m. 9/23/2013 at 5 p.m.

Table 30 shows the backcalculated layer moduli for the accuracy cells, which were estimated using the widely used program MODULUS.(53) Other outputs of the program include standard deviation of the backcalculated moduli computed from multiple FWD drop data. These backcalculated moduli may be considered appropriate for the site at the time of FWD testing. The AC layer moduli shown in table 30 had to be adjusted to actual temperature at the time of the TSDD testing for the purposes of the 3D-Move simulations. The procedure used for temperature adjustment is detailed in the subsection, Pavement Material Properties at the Time of TSDD Trials, later in this chapter.

Table 30. Backcalculated moduli of pavement layers for accuracy cells.
Cell Material Thickness (inches) Average Modulus (ksi) Standard Deviation (ksi) COV (Percent)
3 AC 3 554 34 14.0
Base 43 68.8 13.6 19.8
Subgrade 122.4 17.7 2.2 12.3
19 AC 5 301 65 22.0
Base 31 32 5.8 18.0
Subgrade 18.1 6.1 0.6 10.2
34 AC 4 299 67 22.0
Base 12 15.7 3.1 19.9
Subgrade 46.3 8.5 0.9 10.2

1 inch = 25.4 mm

1 ksi = 6.89 MPa

The properties of the unbound layers (base and subgrade) were generally unaffected by temperature (as long as there was no freezing), but moisture could play a significant role. However, the FWD and TSDD testing were performed within a short period of each other, and a careful review of climate at the site revealed that no significant changes in moisture (i.e., no rain or snow) occurred. Therefore, the FWD backcalculated moduli for the unbound layers at the three MnROAD accuracy cells were used without adjustments in the 3D-Move runs.

Evaluation of AC Pavement Temperature During FWD Testing

Accurate estimates of the AC layer temperatures at the time of FWD testing were necessary since AC stiffness is significantly affected by temperature. Accordingly, adjustments to the AC stiffness to reflect actual temperatures at the time of TSDD testing were needed. A TC tree device was used for the MnROAD accuracy cells to measure temperature within the pavement layers. A typical MnROAD TC tree is 6 ft (1.83 m) long. Sensor depths were selected during construction to provide a temperature profile within the AC, base, and subbase layers. Table 31 shows the TC sensor depths from the surface at the MnROAD accuracy cells, while table 32 shows the AC layer temperatures recorded by the TC device at the three cells.

Table 31. Depths of TC sensors in accuracy cells.
Cell 3 Cell 19 Cell 34
Sensor Depth (inches) Sensor Depth (inches) Sensor Depth (inches)
101 0.5 101 0.5 201 0.5
102 1.5 102 1.5 202 1.6
103 2.5 103 2.5 203 2.5
104 3.5 104 3.5 204 3.5
105 5 105 4.5 205 7
106 7 106 6 206 9
107 9 107 9 207 11
108 11 108 12 208 13
109 13 109 15 209 15
110 15 110 18 210 20
111 24 111 24 211 24
112 28 112 30 212 30
113 36 113 36 213 36
114 48 114 48 214 48
115 60 115 60 215 60
116 72 116 72 216 72

1 inch = 25.4 mm

Note: Bolded cells indicate the sensors were located within the AC layer.

Table 32. Measured temperature within AC at the time of FWD test.
Cell Day Hour Temperature at Location from Surface (ºF)
Sensor 101 Sensor 102 Sensor 103 Sensor 104 Sensor 105 Sensor 201 Sensor 202 Sensor 203 Sensor 204
3 9/16/13 14 126 93 90 84 N/A N/A N/A N/A N/A
19 9/20/13 14 84 ** 135* 73 72 N/A N/A N/A N/A
34 9/5/13 15 N/A N/A N/A N/A N/A ** 108 ** **

1 °F = 1.8° C + 32

*Possible incorrect data.

**Missing data.

N/A = Not applicable.

As shown in table 32, there are missing and possibly incorrect data from the MnROAD TC devices. Although AC temperature depth coverage seemed adequate for cell 3, which had the lowest thickness of 3 inches (76.2 mm), coverage in the other cells were not. For example, in cell 19, there are four data points, but one of them (sensor 103) was unreliable. For cell 34, only one data point is available. Accordingly, an alternate robust and defensible approach was adapted to estimate temperatures within the AC layers during FWD testing and also during the TSDD field trials. Towards this end, the BELLS equation was selected for a number of reasons, including its extensive calibration using data from the Long-Term Pavement Performance (LTPP) Program.(40,54)

The BELLS equation was used to predict temperatures within the AC layer for cells 19 and 34 at locations other than those for which field measurements were available. The BELLS equation is given by the following:

Figure 125. Equation. BELLS equation. T subscript d equals 0.95 plus 0.892 times IR plus the quantity of open bracket logarithm open parenthesis d closed parenthesis minus 1.25, end quantity, closed bracket times the quantity of open bracket -0.448 times IR plus 0.621 times open parenthesis 1-day closed parenthesis plus 1.83 times sine times the quantity of open parenthesis hr subscript 18 minus 15.5, end quantity, end quantity, closed parenthesis closed bracket plus 0.042 times IR times sine times the quantity of open parenthesis hr subscript 18 minus 13.5, end quantity closed parenthesis.

Figure 125. Equation. BELLS equation.

Where:

Td= Pavement temperature at depth d.

IR = Pavement surface temperature.

d = Depth at which pavement temperature is to be predicted.

1-day = Average air temperature the day before testing.

hr18= Time of day, in a 24-h clock system, calculated using an 18-h AC and temperature rise-and-fall time cycle.

Figure 126 shows the variation of air temperature in Monticello, MN, which is close to the MnROAD facility, during September 2013, which is when the FWD testing and TSDD field trials under consideration took place.

Figure 126. Graph. September 2013 temperature variations in Monticello, MN. This graph shows the variation of air temperature in Monticello, MN. The y-axis shows temperature from 10 to 110 °F (-12.2 to 43.3 °C), and the x-axis shows days of the month of September from 1 to 31. Two data points are shown: actual low (squares) and actual hi (circles). Actual low varies from 40 to 60 °F (4.4 to 15.5 °C), and actual hi varies from 65 to 90 °F (18.3 to 32.2 °C). Actual low starts at 50 °F (10 °C), decreases to 40 °F (4.4 °C) on September 3, increases to 60 °F (15.5 °C) on September 6, remains almost same until September 10, then decreases to 40 °F (4.4 °C) on September 13, has small changes until September 16, then increases to 60 °F (15.5 °C) on September 18, then decreases to 40 °F (4.4 °C) on September 21, constantly increases to 60 °F (15.5 °C) on September 27, and decreases to 40 °F (4.4 °C) on September 29, and finally increases to 50 °F (10 °C) on September 30. Actual high has a same trend with an offset of about 20 °F (11.1 °C) except on September 8 when actual high and actual low are close to each other. Actual low and actual high temperatures were 58 and 73 °F (14.44 and 22.78 °C) on September 19 and 53 and 85 °F (11.6 and 29.4 °C) on September 4, respectively.

1 °F = 1.8 °C + 32

Figure 126. Graph. September 2013 temperature variations in Monticello, MN.

Figure 126 was used to find the average air temperature of the day before FWD testing, which was 66 °F (18.9 °C) on September 19, 2013, and 69 °F (20.6 °C) on September 4, 2013, for cells 19 and 34, respectively. Appropriate surface temperatures were determined based on the data available from the top-most sensor (sensor 101 in cell 19 and sensor 202 in cell 34) of the TC device and matching the prediction by the BELLS equation. Having estimated the required surface temperatures, the BELLS equation was used to predict the variation of temperature with depth within the AC layer at both cells in question. Table 33 and table 34 show the estimated temperatures within AC layer for cells 19 and 34 at the time of FWD testing.

Table 33. Temperature within AC from BELLS equation in cell 19 during FWD test.
Matched Surface Temperature (°F) Time Depth (inches) 1-Day Temperature (°F) Temperature from BELLS Equation (°F)
87.3 2 p.m. 0.5 65.5 84.2
2 p.m. 1.5 65.5 82.4
2 p.m. 2.5 65.5 80.6
2 p.m. 3.5 65.5 78.8
2 p.m. 4.5 65.5 78.8

1 °F = 1.8 °C + 32

1 inch = 25.4 mm

Table 34. Temperature within AC from BELLS equation in cell 34 during FWD test.
Matched Surface Temperature (°F) Time Depth (inches) 1-Day Temperature (°F) Temperature from BELLS Equation (°F)
119.3 (48.5) 3 p.m. 0.5 68.9 116.6
3 p.m. 1.6 68.9 107.6
3 p.m. 2.5 68.9 104
3 p.m. 3.5 68.9 102.2

1 °F = 1.8 °C + 32

1 inch = 25.4 mm

The data presented in table 32 to table 34 were then used to estimate representative average AC layer temperatures at the time of FWD testing for cells 3, 19, and 34, which resulted in temperatures of 99, 81, and 108 °F (37.2, 27.2, and 42.2 °C), respectively.

Pavement Material Properties at the Time of TSDD Trials

The 3D-Move program requires viscoelastic characterization (frequency dependent) of the AC layers to evaluate pavement responses as a function of vehicle speed. Since the TSDD field trials were conducted under different temperature conditions, a master curve approach was adopted so that the AC layer modulus versus frequency relationship could be determined at any given temperature.

The FWD backcalculated AC layer modulus was appropriate for the temperature at the time of testing and a loading frequency of about 30 Hz.(55) Using this AC layer modulus and frequency as an anchoring point, the AC master curve (i.e., modulus versus frequency) could be established.

The first step in the development of the AC master curves was determining the representative AC layer temperatures during TSDD testing. The same procedure used to determine the AC layer temperatures during FWD testing was used to estimate the average AC layer temperatures during TSDD testing. Table 35 summarizes the average AC layer temperatures at the time of FWD and TSDD testing.

Table 35. Average temperatures within AC layer.
Cell Temperature at Time of FWD Testing (°F) Temperature at Time of TSD Testing (°F) Temperature at Time of RWD Testing (°F)
3 99 91 99
19 81 68 63
34 108 91 90

1 °F = 1.8 °C + 32

Using the average AC layer temperature, damping was estimated using dynamic modulus test data available in the MnROAD database. In those tests, the phase angels of the AC layer at several temperatures and frequencies were measured. In its complex form, the dynamic modulus is given by the following:

Figure 127. Equation. Complex dynamic modulus. E star equals E prime times the quantity of open parenthesis 1 plus 2 times i times zeta subscript AC, end quantity, closed parenthesis.

Figure 127. Equation. Complex dynamic modulus.

Where:

E* = Complex dynamic modulus.

E' = Elastic (or storage) modulus.

i = Unit imaginary number.

ζAC= Measure of internal damping of AC.

For viscoelastic layers, E* can be presented as the sum of the real and imaginary components, which is given by the following:

Figure 128. Equation. Modified modulus equation. E star equals E prime plus i times E double prime.

Figure 128. Equation. Modified modulus equation.

Where:

E" = Damping (or loss) modulus.

The equations may be re-written as follows:

Figure129. Equation. Rewritten modulus equation. Zeta subscript AC equals E double prime divided by the quantity of open parenthesis 2 times E prime, end quantity, closed parenthesis equals 0.5 times tan times open parenthesis phi closed parenthesis.

Figure 129. Equation. Rewritten modulus equation.

Where:

φ= Phase angle.

Table 36 and table 37 show the values of damping associated with the TSDD field trials. The procedure used to obtain these values consisted of fitting a best curve through the available AC dynamic modulus test (phase angle measured at given frequencies) results and then interpolating the phase angles for the temperature corresponding with the TSDD field trials.

Table 36. Phase angle and damping for accuracy cells for TSD field trials.
Cell 3 (91 °F) Cell 19 (68 °F) Cell 34 (91 °F)
Frequency (Hz) Phase Angle (Degree) Damping (Percent) Frequency (Hz) Phase Angle (Degree) Damping (Percent) Frequency (Hz) Phase Angle (Degree) Damping (Percent)
0.1 36.3 36.8 0.1 34.2 33.9 0.01 29.0 27.8
0.5 36.7 37.3 0.5 26.4 24.8 0.1 32.0 31.2
1 36.2 36.6 1 23.7 21.9 1 31.6 30.7
5 33.2 32.7 5 19.1 17.3 25 27.6 26.1
10 34.4 34.2 10 16.9 15.1 N/A N/A N/A
25 31.6 30.7 25 13.5 11.9 N/A N/A N/A

1 °F = 1.8 °C + 32

N/A = Not available.

Table 37. Phase angle and damping for accuracy cells for RWD field trials.
Cell 3 (99 °F) Cell 19 (63 °F) Cell 34 (90 °F)
Frequency (Hz) Phase Angle (Degree) Damping (Percent) Frequency (Hz) Phase Angle (Degree) Damping (Percent) Frequency (Hz) Phase Angle (Degree) Damping (Percent)
0.1 36.8 37.4 0.1 32.7 32.1 0.01 29.7 28.5
0.5 36.6 37.2 0.5 22.6 20.8 0.1 33.0 32.4
1 36.4 36.9 1 19.4 17.6 1 32.5 31.9
5 34.7 34.6 5 15.4 13.7 25 28.0 26.5
10 36.7 37.2 10 12.4 10.9 N/A N/A N/A
25 34.6 34.5 25 9.4 8.3 N/A N/A N/A

1 °F = 1.8 °C + 32

N/A = Not available.

Undamaged Modulus of AC Layer:

The procedure used to develop the appropriate AC master curves took into account the following considerations:

One of the most comprehensive AC mixture stiffness models is Witczak's dynamic modulus predictive equation. This model predicts modulus as a function of temperature and frequency based on volumetric AC mix design information. The Witczak-Andrei AC dynamic modulus equation is given by the following:(56)

Figure 130. Equation. Witczak equation. Logarithm E star equals -1.25 plus 0.029 times rho subscript 200 minus 0.0018 times open parenthesis rho subscript 200 closed parenthesis raised to the power of 2 minus 0.0028 times rho subscript 4 minus 0.058 times V subscript a minus 0.822 times V subscript beff divided by the quantity of V subscript beff plus V subscript a, end quantity, plus the quantity of 3.872 minus 0.0021 times rho subscript 4 plus 0.003958 times open parenthesis rho subscript 38 closed parenthesis minus 0.000017 times open parenthesis rho subscript 38 closed parenthesis raised to the power of 2 plus 0.0055 times rho subscript 34, end quantity, divided by quantity of 1 plus e raised to the power of the quantity of open parenthesis -0.603313 minus 0.313351 times logarithm open parenthesis f closed parenthesis minus 0.393532 times logarithm open parenthesis eta, end quantity, end quantity closed parenthesis closed parenthesis.

Figure 130. Equation. Witczak equation.

Where:

ρ200 = Percent passing #200 (0.075 mm) sieve.

ρ4 = Cumulative percent retained on #4 (4.76 mm) sieve.

ρ38 = Cumulative percent retained on 0.375-inch (9.5 mm) sieve.

ρ34 = Cumulative percent retained on 0.75-inch (19 mm) sieve.

Va = Air void, percent by volume.

Vbeff = Effective binder content, percent by volume.

e = Euler's number.

f = Loading frequency, Hz.

η = Viscosity (cP).

The Witczak equation requires input such as gradation and viscosity. The gradation data for the three MnROAD accuracy cells were available from the MnROAD database, while viscosity values were estimated as a function of temperature based on the regression intercept, A, and the regression slope of viscosity temperature susceptibility (VTS) as shown in figure 131.(58)

Figure 131. Equation. Temperature-viscosity relationship. Logarithm logarithm of eta equals A plus VTS times logarithm T subscript R.

Figure 131. Equation. Temperature-viscosity relationship.

Where:

A = Regression intercept.

TR = Temperature at which the viscosity is estimated (Rankine).

AASHTO T315-10 gives guidelines for calculating A and VTS.(59) The dynamic shear rheometer test results are available in the MnROAD databases for the three cells under consideration. The performance grade (PG) for these cells, which were estimated from available data, are PG 70-16, PG 64-22, and PG 58-34 for cells 3, 19, and 34, respectively. Accordingly, the calculated A and VTS values were 10.641 and -3.548 for cell 3, 10.98 and -3.68 for cell 19, and 10.149 and -3.359 for cell 34.

Having determined the A and VTS values, the equation given in figure 130 was used to estimate the undamaged AC modulus as a function of temperature and frequency. The next step entailed the derivation of the existing AC modulus at various frequencies and at the AC layer temperature corresponding to the time of the TSDD testing based on the FWD backcalculated layer moduli. This step is explained next.

Obtaining Existing AC Modulus for Use with 3D-Move:

The existing AC moduli as a function of frequency was estimated by using the backcalculated AC layer moduli as anchor points and shifting the Witczak-derived AC modulus relationships as shown in figure 132. In this figure, ENDT and tr are the AC modulus and reduced time from the non-destructive testing (NDT), respectively. EPRED is the AC modulus from the prediction equation and δ is the minimum AC modulus.

Figure 132. Graph. Development of existing AC modulus master curves from undamaged AC moduli. This illustration shows the development of existing asphalt concrete (AC) modulus master curves from undamaged AC moduli. The y-axis shows modulus in logarithmic scale, and the x-axis shows reduced time in logarithmic scale. The graph includes two curves. The first curve is the undamaged master curve, which has a high value at low x values and decreases to the lowest value but with smooth variations in both the low and high x values. This curve can be produced using the Witczak prediction equation. The second curve is the damaged master curve, which has the same trend as the undamaged master curve but smaller values in all x values. The damaged master curve reaches undamaged master curve at high x values. The difference between the highest and lowest values is specified by alpha in the undamaged master curve and alpha prime in the damaged master curve. The lowest value of the two curves is specified by delta. At one specific time or t subscript r, the non-destructive test (NDT) is done and corresponding modulus on damaged master curve and undamaged master curve are E subscript NDT and E subscript PRED or predicted modulus.

Figure 132. Graph. Development of existing AC modulus master curves from undamaged AC moduli.

The AC existing modulus master curve was obtained by applying the following equation to E* computed from the original master curve:

Figure 133. Equation. Existing modulus equation. E star subscript dam equals 10 raised to the power of delta plus the quantity of E star minus 10 raised to the power of delta, end quantity, divided by the quantity of 1 plus e raised to the power of the quantity of -0.3 plus 5 times logarithm open parenthesis d subscript AC closed parenthesis, end quantity, end quantity.

Figure 133. Equation. Existing modulus equation.

Where:

E*dam = Existing modulus.

δ= Regression parameter (from master curve), it is also the minimum AC modulus.

dAC = Fatigue damage in AC layer.

As stated earlier, the AC moduli obtained from the backcalculation of FWD deflection data reflect the temperatures at the time of FWD testing as well as for a frequency of about 30 Hz. Figure 132 illustrates the procedure used to derive the AC existing modulus master curves from the undamaged AC moduli.

Accordingly, using Witczak's equation and taking into consideration the fatigue damage factor, a reasonable AC modulus master curve can be generated. Moreover, this approach is similar to that proposed by Seo et al.(60)

AC Master Curves During TSDD Field Trials:

AC modulus master curves for the three MnROAD accuracy cells were developed using the procedure described in the previous subsection. Figure 134 through figure 136 show the resulting master curves for the temperatures associated with the TSDDs field tests. These curves appear realistic, with smooth variation in both the low and high frequencies. Hence, they were used as input to 3D-Move.

Figure 134. Graph. Master curves for cell 3 in TSDD field trials. This graph shows the resulting master curves for cell 3 for the temperatures associated with the traffic speed deflection device (TSDD) field tests. The y-axis shows logarithm E star from 2 to 7 psi (13.78 to 48.23 kPa), and the x-axis shows logarithm f from -20 to 15. Two lines are shown: Traffic Speed Deflectometer with a temperature of 91 °F (32.78 °C) and Rolling Wheel Deflectometer with a temperature of 99 °F (37.22 °C). Both curves begin at (-15, 3) and end around (13, 6.75). These curves show smooth variation in both the low and high frequencies.

1 psi = 6.89 kPa

1 °F = 1.8 °C + 32

Figure 134. Graph. Master curves for cell 3 in TSDD field trials.

Figure 135. Graph. Master curves for cell 19 in TSDD field trials. This graph shows the resulting master curves for cell 19 the temperatures associated with the traffic speed deflection device (TSDD) field tests. The y-axis shows logarithm E star from 2 to 7 psi (13.78 to 48.23 kPa), and the x-axis shows logarithm f from -20 to 15. Two lines are shown: Traffic Speed Deflectometer with a temperature of 68 °F (20 °C) and Rolling Wheel Deflectometer with a temperature of 63 °F (17.22 °C). Both curves begin at (-15, 3) and end around (13, 6.4). These curves show smooth variation in both the low and high frequencies.

1 psi = 6.89 kPa

1 °F = 1.8 °C + 32

Figure 135. Graph. Master curves for cell 19 in TSDD field trials.

Figure 136. Graph. Master curves for cell 34 in TSDD field trials. This graph shows the resulting master curves for cell 34 for the temperatures associated with the traffic speed deflection device (TSDD) field tests. The y-axis shows logarithm E star from 2 to 7 psi (13.78 to 48.23 kPa), and the x-axis shows logarithm f from -20 to 15. Two lines are shown: Traffic Speed Deflectometer with a temperature of 91 °F (32.78 °C) and Rolling Wheel Deflectometer with a temperature of 90 °F (32.22 °C). Both curves begin at (-15, 3) and end around (13, 6.5). These curves show smooth variation in both the low and high frequencies.

1 psi = 6.89 kPa

1 °F = 1.8 °C + 32

Figure 136. Graph. Master curves for cell 34 in TSDD field trials.

Summary

This section detailed the procedure that was used to obtain the pavement layer material properties, including viscoelastic characterization of the AC layers. The procedure relied on three data sources for estimating material properties: (1) laboratory testing, which included inspection, sampling, and testing data for various pavement materials used at MnROAD; (2) MnROAD cell specifications, which included layer thickness and type of layers; and (3) FWD deflection test results, which were done a few days before the TSDD field trials. Accordingly, it was possible to estimate the existing material layer properties for the three MnROAD accuracy cells. For unbound materials, layer moduli backcalculated from the FWD deflection measurements (see table 30) were used without modification. Table 38 presents the values used in 3D-Move for the unbound layers properties.

Table 38. Unbound layer properties used in 3D-Move.
Cell Layer Modulus (ksi) Thickness (inches) Poisson's Ratio Unit Weight (pcf) Damping Ratio (Percent)
3 Base 68.8 43 0.35 105 5
Subgrade 17.7 122.4 0.4 105 5
19 Base 32 31 0.35 105 5
Subgrade 6.1 18.1 0.4 105 5
34 Base 15.7 12 0.35 110 5
Subgrade 8.5 46.3 0.4 110 5

1 ksi = 6.89 MPa

1 inch = 25.4 mm

1 pcf = 16 Kg/m3

Table 39 provides the AC layer moduli corresponding to the temperature at the time of the TSDD field trials at frequency of 30 Hz, where standard deviation was obtained from MODULUS backcalculation shown in table 30.

Table 39. AC modulus at the time of TSDD trials for frequency of 30 Hz.
Parameter Device Cell 3 Cell 19 Cell 34
Average AC temperatures, °F RWD 99 63 90
TSD 91 68 91
Existing modulus at f = 30 Hz (mean), ksi RWD 554 550.1 522
TSD 716.8 468.5 506.5
Existing modulus at f = 30 Hz (mean minus standard deviation), ksi RWD 520 431.2 404.9
TSD 672.8 367.2 392.9

1 °F = 1.8° C + 32

1 ksi = 6.89 MPa

7.4 Loading Characteristics of TSDDs

Characterization of the applied load is an important input for analyzing pavement structures. During the TSDD field trials, the loads carried by the axles were determined using a static scale owned and operated by MnROAD. Details concerning the axle configuration and contact pressure distribution of the TSD and RWD are presented next.

Because the weights of the devices were measured using a static scale, the 3-D Move results should be used with caution when simulating moving vehicles. Figure 137 illustrates the variation in the rear trailing axle load for a five-axle truck semitrailer traveling 25–30 mi/h (4.25–48.3 km/h). The dynamic load could vary by as much as 33 percent of the static load. This suggests that the variation of axle load, which has a direct influence on the computed deflection response, should be addressed.

Figure 137. Graph. Variation of five-axle truck semitrailer. This graph illustrates the variation in the rear trailing axle load (right and left wheel) for a five-axle truck semitrailer traveling at 24.85 to 31.07 mi/h (40 to 50 km/h). The y-axis shows dynamic wheel load from 0 to 60 kN (0 to 13,636 lbf), and the x-axis shows time from 100 to 220 s. The right wheel trend is represented by circles over a solid line, and the left wheel trend is represented by a solid line. The average of the dynamic load of the right and left wheel in one particular time is close to static load (F subscript stat), which is 39 kN (8,863 lbf). The variations of the dynamic load (F subscript dyn) of the right wheel shows that it often has a smaller value than static load and reaches the minimum value of 25 kN (5,682 lbf) or 64 percent of the static load in some time. The variation of the dynamic load of left wheel shows that it often has higher value than static load and reaches the maximum value of 52 kN (11,818 lbf) or 133 percent of the static load in some time.

1 mi/h = 1.61 km/h

1 lbf = 0.0044 kN

Figure 137. Graph. Variation of five-axle truck semitrailer.(61)

Axle Configurations and Contact Pressure Distribution

Figure 138 shows TSD axle configuration and axle loads, while figure 139 shows RWD axle configuration and axle loads. Since the deflection measuring sensors for the two TSDDs were mounted close to the dual tire of the rear trailing axles (shown in red box in the figures), the 3D-Move comparisons in this calibration effort focused on the responses generated by those axles. Because data on the pavement-tire contact pressure distribution were not available, the rear axles were modeled as dual circular loads with uniform contact pressure in the 3-D Move analyses. Based on the information presented in figure 137, the static loads measured with the two TSDDs were increased by 25 percent in two 3D-Move scenarios to explore the impact of axle load variations caused by the dynamic interaction with pavements.

Figure 138. Illustration. TSD axle configuration and loads. This illustration shows Traffic Speed Deflectometer (TSD) axle configuration and loads. Tire pressure and load for rear axle are 116 psi (800 kPa) and 5,575 lb (2,528.7 kg) per tire, respectively. The tire span of dual tires is 13.5 inches (342.9 mm) in the axle. Based on information from the manufacturer, the tire size code is 275/70 R22.5. The dual tire at the right rear end is enclosed in a red box. The distance between the outside of the tires in the rear axle is 95 inches (2,413 mm). The distance between the outside of the tires in the front axle is 89 inches (2,260.6 mm). The front axle load is 5,450 lb (2,472.1 kg) per tire. The load for the first driver axle and second driver axle are 2,437.5 and 2,525 lb (1,105.6 and 1,145.3 kg) per tire, respectively. The distances between the rear and second driver and first driver axles and front axels are 554, 252, and 200 inches (14,071.6, 6,400.8, and 5,080 mm), respectively.

1 inch = 25.4 mm

1 T = 0.907 Mg

1 lb = 0.454 kg

1 psi = 6.89 kPa

Figure 138. Illustration. TSD axle configuration and loads.

Figure 139. Illustration. RWD axle configuration and loads. This figure shows Rolling Wheel Deflectometer (RWD) axle configuration and loads. Tire pressure and load for the rear axle are 100 psi (690 kPa) and 4,750 lb (2,154.5 kg) per tire, respectively. The dual tire at the rear right side is enclosed in a red box. The tire span of dual tires is 14.5 inches (368.3 mm) in the axle. The distance between centers of dual tires in the rear axle is 74 inches (1,879.6 mm). The distance between the centers of the tires in the steering axle is 82 inches (2,082.8 mm). The steering axle load is 4,650 lb (2,109.2 kg) per tire. The load for the driver axles is 2,900 lb (1,315.4 kg) tire. The distance between the rear axle and the first driver axle is 473 inches (12,014.2 mm), the distance between the first driver axle and the second driver axle is 51 inches (1,295.4), and the distance between the second driver axle and the steering axle is 204 inches (5,181.6 mm).

1 inch = 25.4 mm

1 T = 0.907 Mg

1 lb = 0.454 kg

1 psi = 6.89 kPa

Figure 139. Illustration. RWD axle configuration and loads.

7.5 3D-Move Calibration using Project Sensor Data

Project Sensors

As discussed in chapter 5, four geophones and one accelerometer were installed near the pavement surface at each of the MnROAD accuracy cells. Since the measurement sensors for both TSDDs were located along the midline of the rear axle tires, two geophones (GEO 1 and GEO 3) were located along this midline plane, as shown in figure 140. To estimate the lateral wheel location relative to installed sensors, two video cameras were also used to record each pass of the TSDDs. A laser trigger device was also installed across the travel lane (transverse plane) for timing purposes. The location of the rear wheels relative to the location of the maximum displacement (i.e., lag distance) could be determined by superimposing the data from the laser device and measured displacement bowls. As noted in chapter 5, the geophone data were filtered and analyzed to arrive at the final surface deflection data. The TSDD field trials were conducted at various vehicle speeds, as summarized in table 40. Figure 141 provides a sketch of the pavement sections and average FWD backcalculated layer moduli and the standard deviation of cells 3, 19, and 34. These pavement structures represent pavements with AC thickness varying from 3–5 inches (76.2–127 mm) and a base thickness varying from 12–43 inches (304.8– 1,092.2 mm).

Figure 140. Illustration. Configuration and spacing of embedded project sensors. This illustration shows four geophones installed at each of the MnROAD accuracy cells. The direction of the vehicle pass is aligned from left to right, which is the same as the sequence of geophones labeled 1 to 4. Two geophones (GEO 1 and GEO 3) are placed along the midline of the inner and outer wheel, which is the plane of measurement for both traffic speed deflection device. The others geophones (GEO 2 and GEO 4) are offset by 6 inches (152.4 mm) on either side of them. The longitudinal distance between geophones is 1 ft (0.305 m). A laser beam is shown as being perpendicular to GEO 3.

1 inch = 25.4 mm

1 ft = 0.305m

Figure 140. Illustration. Configuration and spacing of embedded project sensors.

Table 40. TSDD accuracy runs at MnROAD.

Cell TSDD Velocity (mi/h)
3 TSD 30 and 45
RWD 30, 45, and 60
19 TSD 30, 45, and 60
RWD 30, 45, and 60
34 TSD 30 and 45
RWD 30 and 45

1 mi/h = 1.61 km/h

Figure 141. Illustration. Sketch of pavement structures used for cells 3, 19, and 34. This illustration provides a sketch of the pavement sections and falling weight deflectometer (FWD) backcalculated layer moduli of cells 3, 19, and 34. These pavement structures represent pavements with asphalt concrete (AC) thicknesses of 3, 4, and 5 inches (76.2, 101.6, and 
127 mm) for cells 3, 34, and 19, respectively, base thicknesses of 12, 31, and 43 inches (304.8, 787.4, and 1,092.2 mm) for cells 34, 19, and 3, respectively, and subgrade thicknesses of 18.1, 46.3, and 122.4 inches (460, 1,176, and 3,109 mm) for cells 19, 34, and 3, respectively. AC FWD modulus for cells 3, 19, and 34 are 554 ksi (3,820 MPa) with a standard deviation of 34 ksi (234 MPa), 301 ksi (2,075 MPa) with a standard deviation of 65 ksi (448.2 MPa), and 299 ksi (2,062 MPa) with a standard deviation of 67 ksi (462 MPa), respectively. The base modulus for cells 3, 19, and 34 are 68.8 ksi (474.4 MPa) with a standard deviation of 8.5 ksi (58.6 MPa), 
32 ksi (220.6 MPa) with a standard deviation of 5.8 ksi (40 MPa), and 15.7 ksi (108.2 MPa) with a standard deviation of 3.1 ksi (21.3 MPa), respectively. Subgrade modulus for cells 3, 19, and 34 are 17.7 ksi (122 MPa) with a standard deviation of 2.9 ksi (20 MPa), 6.1 ksi (42 MPa) with a standard deviation of 0.6 ksi (4.1 MPa), and 8.5 ksi (58.6 MPa) with a standard deviation of 0.9 ksi (6.2 MPa), respectively.

1 ksi = 6.89 MPa

1 inch = 25.4 mm

σ = Standard deviation.

Figure 141. Illustration. Sketch of pavement structures used for cells 3, 19, and 34.

Case Scenarios for 3D-Move to Bracket the Project Sensor Results

Initial attempts revealed that surface displacements computed from the 3D-Move analyses with mean tire loads and material properties did not match well with those given by the project sensors. The 3D-Move displacement basins were wider and of a smaller magnitude for all three accuracy cells and passes. Many possible factors might have contributed to this. They include the following:

Accordingly, after an extensive trial and error process, acceptable bracketing of the measured deflection responses was obtained with the following three 3D-Move case scenarios:

Case X was only used for cell 19 to bracket the measured values. As is shown later in this chapter, the 3D-Move predictions based on the above three scenarios bracketed well the measured responses (stresses, strains, and displacements) collected by the project and MnROAD sensors as well as with those measured by TSDDs.

3D-Move Results for Accuracy Runs

The results from the analyses with 3D-Move using the case scenarios described in the previous section were initially compared with the measured values from the project sensors. TSDD measurements were also included in the comparisons. It should also be noted that TSD provided surface velocities at six points ahead of the wheels, while RWD measured displacement at two points.

Figure 142 illustrates the comparison of 3D-Move results and measured deflections for cell 34 based on testing done by the TSD at a vehicle velocity of 30 mi/h (48.3 km/h). Since the main focus is the comparison of the deflection bowls (shape and maximum value), the plots were shifted to align so that the maximum displacement locations coincided. Since GEO 1 and GEO 3 are located along a plane parallel to the vehicle direction, they were expected to have similar deflection bowls. Therefore, the responses from these two sensors are shown in the figure. The variation between the deflection bowls given by these two sensors may be viewed as a measure of the overall variability in the measurements made by the project sensors and possibly any spatial variability between 2 ft (0.61 m) of distance between geophones (figure 140). The lower and upper bound of the project sensor data were determined by treating GEO 1 and GEO 3 data as independent datasets in figure 142 and figure 143. In all cases, 3D-Move adequately captured the maximum and shape of measured displacements. Case X1 provided the closest deflection basin to that measured. Figure 143 compares predicted and measured TSD surface velocities (maximum from field trials). While all 3D-Move case results had similar shapes as the measured ones, the case X1 predictions were closer to both project sensors and TSD measurements.

Figure 142. Graph. 3D-Move predictions and measured deflections for cell 34 in TSD trials (device velocity = 30 mi/h (48 km/h)). This graph illustrates the comparison of 3D-Move predicted (3D-Move/case 1 and 3D-Move/case X1) and measured deflections (GEO1/GEO 3/ lower bound and GEO1/GEO 3/upper bound) for cell 34 based on testing done by the Traffic Speed Deflectometer (TSD) at a vehicle velocity of 30 mi/h (48 km/h). The y-axis shows deflection from -39 to 1 mil (-0.975 to 0.025 mm), and the x-axis shows distance from -10 to 10 ft (-3.05 to 3.05 m). In all cases, 3D-Move adequately captured the bell shape of measured displacements. 3D-Move/case X1 provides the closest deflection to that measured deflection in GEO1/GEO 3/upper bound which, is 36 mil (0.91 mm). The maximum deflection in 3D-Move/case 1 is 27 mil (0.686 mm). The maximum deflection in GEO1/GEO 3/lower bound is 32 mil (0.813 mm).

1 mil = 0.025 mm

1 ft = 0.305 m

Figure 142. Graph. 3D-Move predictions and measured deflections for cell 34 in TSD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 143. Graph. 3D-Move predictions and measured velocities for cell 34 in TSD trials 
(device velocity = 30 mi/h (48 km/h)). This graph illustrates the comparison of 3D-Move predicted (3D-Move/case 1 and 3D-Move/case X1) and measured velocities (GEO1/GEO3/upper bound, GEO1/GEO3/lower bound, and Traffic Speed Deflectometer (TSD)) for cell 34 based on testing done by the TSD at a vehicle speed of 30 mi/h (48 km/h). The y-axis shows velocity from -1,500 to 1,500 mil/s (-38.1 to 38.1 mm/s), and the x-axis shows distance from -10 to 10 ft (-3.05 to 3.05 m). In both cases (case 1 and case X1), 3D-Move adequately captures the shape of measured displacements. Also, the velocity decreases from 0 mil/s (0 mm/s) to the minimum value, then increases sharply, and then decreases to 0 mil/s (0 mm/s). The minimum velocity ranges from about -600 to -1,000 mil/s (-15.24 to -25.4 mm/s) for all cases. The 3D-Move/case X1 provides the closest deflection velocity to TSD measurements and project sensors. The maximum velocity ranges from about 700 to 1,100 mil/s (17.8 to 27.9 mm/s).

1 mil/s = 0.0254 mm/s

1 ft = 0.305 m

Figure 143. Graph. 3D-Move predictions and measured velocities for cell 34 in TSD trials (device velocity = 30 mi/h (48.3 km/h)).

To compare the shapes of the computed and measured basins, normalized deflection basins were used. This was done by dividing each deflection basin by its maximum value. Figure 144, which shows normalized basins, reveals that the differences in shape among the computed and measured deflection basins were small. A comparison between the normalized basins can be achieved by using the widths of the deflection basins at three levels as shown in figure 145. Figure 146 shows the comparison of the pulse widths at various levels of normalized deflection. This figure shows that 3D-Move produced comparable normalized deflection basins. Similar plots were developed for all possible combinations of TSDD speeds and accuracy cells, and they showed the same trends as shown in the example figures (i.e., figure 144 and figure 146 through figure 150) included within this chapter. The data used to generate all of these figures may be obtained directly from FHWA.

Figure 144. Graph. Normalized basins of 3D-Move predictions and measured deflections for cell 34 in TSD trials (device velocity = 30 mi/h (48 km/h)). This graph shows normalized basins of 3D-Move predicted (3D-Move/case 1 and 3D-Move/case X1) and measured deflections (GEO1/GEO3/upper bound and GEO1/GEO3/lower bound) for cell 34 based on testing done by the Traffic Speed Deflectometer (TSD) at a vehicle speed of 30 mi/h (48 km/h). The y-axis shows the normalized deflection from -1.1 to 0.1 mil/mil (-1.1 to 0.1 mm/mm), and the x-axis shows distance from -10 to 10 ft (-3.05 to 3.05 m). The predicted and measured deflection shapes agree closely. This graph reveals that the differences in shape among the computed and measured deflection basins are small at the center but vary slightly as distance increases from the center.

1 mil/mil = 1 mm/mm

1 ft = 0.305 m

Figure 144. Graph. Normalized basins of 3D-Move predictions and measured deflections for cell 34 in TSD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 145. Graph. Definition of percentage of maximum displacement. This graph defines three levels (25, 50, and 75 percent) of maximum displacement in the normalized basins which were used for comparisons of width of normalized basins between 3D-Move and project 
sensors in these levels. The y-axis shows normalized deflection from -1.1 to 0.1 mil/mil (-1.1 to 
0.1 mm/mm), and the x-axis shows distance from -10 to 10 ft (-3.05 to 3.05 m). The widths of the normalized deflection basin, which is a bell-shaped curve, at the three levels are shown. These levels are defined based on normalized deflection at 75, 50, and 25 percent of maximum displacement.

1 mil/mil = 1 mm/mm

1 ft = 0.305 m

Figure 145. Graph. Definition of percentage of maximum displacement.

Figure 146. Graph. Comparison of normalized deflection basins at various levels for cell 34 in TSD trials (device velocity = 30 mi/h (48 km/h)). This bar graph shows the comparison of the pulse widths of normalized basins of 3D-Move predicted (case 1 and case X1) and measured deflections (GEO1 and GEO3) at three levels of maximum displacement percentage for cell 34 based on testing done by the Traffic Speed Deflectometer at a vehicle speed of 30 mi/h (48 km/h). The y-axis shows pulse width from 0 to 4.5 ft (0 to 1.37 m), and the x-axis shows percentage of maximum displacement which are 75, 50, and 25 percent. A similar trend is visible at all three levels, and the 3D-Move (case 1 and case X1) produces close normalized deflection basins compared to the measured ones (GEO1 and GEO3). Maximum measured pulse width varies from 1.5 to 3.6 ft (0.457 to 1.1 m), and maximum computed pulse width varies from 1.4 to 4.2 ft (0.43 to 1.28 m). Measured and computed pulse widths at 75 percent of maximum displacement are closer than those at 50 and 25 percent of maximum displacement.

1 ft = 0.305 m

Figure 146. Graph. Comparison of normalized deflection basins at various levels for cell 34 in TSD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 147 shows the comparison of 3D-Move and measured deflection basins for cell 19 during TSD testing at 30 mi/h (48.3 km/h). Unlike cell 34, case X yielded closer results to those measured. Figure 148 compares the prediction and measured surface velocities. Figure 149 and figure 150 present the comparison of the normalized deflection basins. Again, a close match can be observed between the 3D-Move predictions and those measured.

Figure 147. Graph. 3D-Move predictions and measured deflections for cell 19 in TSD trials (device velocity = 30 mi/h (48 km/h)). This graph illustrates the comparison of 3D-Move predicted (3D-Move/case 1, 3D-Move/case X, and 3D-Move/case X1) and measured deflections (GEO1/GEO 3/lower bound and GEO1/GEO 3/upper bound) for cell 19 based on testing done by the Traffic Speed Deflectometer (TSD) at a vehicle velocity of 30 mi/h (48 km/h). The y-axis shows deflection from -25 to 5 mil (-0.635 to 0.127 mm), and the x-axis shows distance from  -10 to 10 ft (-3.05 to 3.05m). In all cases, 3D-Move adequately captures the bell shape of measured displacements. 3D-Move/case X provides the closest deflection to that measured, which is 21 mil (0.533 mm). The maximum deflection in 3D-Move/case 1 and case X1 are  13 and 18 mil (0.33 and 0.457 mm), respectively.

1 mil = 0.0254 mm

1 ft = 0.305 m

Figure 147. Graph. 3D-Move predictions and measured deflections for cell 19 in TSD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 148. Graph. 3D-Move predictions and measured velocities for cell 19 in TSD trials (device velocity = 30 mi/h (48 km/h)). This graph illustrates the comparison of 3D-Move predicted (3D-Move /case 1, 3D-Move/case X, and 3D-Move/case X1) and measured velocities (GEO1/GEO3/upper bound, GEO1/GEO3/lower bound, and Traffic Speed Deflectometer (TSD)) for cell 19 based on testing done by the TSD at a vehicle speed of 30 mi/h (48 km/h). The y-axis shows velocity from -800 to 1,000 mil/s (-20.32 to 25.4 mm/s), and the x-axis shows distance from -6 to 6 ft (-1.83 to 1.83 m). In all cases, 3D-Move adequately captures the shape of measured displacements. Additionally in all cases, the velocity decreases from 0 mil/s (0 mm/s) to a minimum value, then increases sharply, and then decreases to 0 mil/s (0 mm/s). The minimum velocity ranges from about -400 to -650 mil/s (-10.16 to 16.5 mm/s) for all cases. The 3D-Move/case X provides the closest deflection velocity to TSD measurements and project sensors. The maximum velocity ranges from about 400 to 800 mil/s (10.16 to 20.32 mm/s).

1 mil/s = 0.0254 mm/s

1 ft = 0.305 m

Figure 148. Graph. 3D-Move predictions and measured velocities for cell 19 in TSD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 149. Graph. Normalized basins of 3D-Move predictions and measured deflections for cell 19 in TSD trials (device velocity = 30 mi/h (48 km/h)). This graph shows normalized basins of 3D-Move predicted (3D-Move/case 1, 3D-Move/case X, and 3D-Move/case X1) and measured deflections (GEO1/GEO3/upper bound and GEO1/GEO3/lower bound) for cell 19 based on testing done by the Traffic Speed Deflectometer (TSD) at a vehicle speed of 30 mi/h (48 km/h). The y-axis shows normalized deflection from -1.2 to 0.2 mil/mil (-1.2 to 0.2 mm/mm), and the x-axis shows distance from -10 to 10 ft (-3.05 to 3.05m). The predicted and measured deflection shapes agree closely. This graph reveals that the differences in shape among the computed and measured deflection basins are small at the center but vary slightly as distance increases from the center.

1 mil /mil = 1 mm/mm

1 ft = 0.305 m

Figure 149. Graph. Normalized basins of 3D-Move predictions and measured deflections for cell 19 in TSD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 150. Graph. Comparison of normalized deflection basins at various levels for cell 19 in TSD trials (device velocity = 30 mi/h (48 km/h)). This bar graph shows the comparison of the pulse widths of normalized basins of 3D-Move predicted (case 1, case X, and case X1) and measured deflections (GEO 1 and GEO 3) at three levels of maximum displacement for cell 19 based on testing done by the Traffic Speed Deflectometer (TSD) at a vehicle speed of 30 mi/h (48 km/h). The y-axis shows the pulse width from 0 to 5 ft (0 to 1.524 m), and the x-axis shows percentage of maximum displacement, which are (from left to right) 75, 50, and 25 percent. This graph shows a similar trend at all three levels, and the predicted deflections produce close normalized deflection basins compared to the measured ones. Maximum measured pulse width varies from 1.2 to 3.6 ft (0.366 to 1.1 m), and maximum computed pulse width varies from 1.2 to 4.5 ft (0.366 to 1.37 m). Measured and computed pulse widths at 75 percent of maximum displacement are closer than those at 50 and 25 percent of maximum displacement.

1 ft = 0.305 m

Figure 150. Graph. Comparison of normalized deflection basins at various levels for cell 19 in TSD trials (device velocity = 30 mi/h (48.3 km/h)).

The 3D-Move results were also compared with those from the RWD device. Again as a representative plot, figure 151 shows 3D-Move and measured deflection bowls as well as the two measured deflection values from RWD. 3D-Move results fall within the geophone measured responses, and the RWD measurements were close to the predicted results. The pulse widths of the normalized deflection basins shown in figure 152 and figure 153 were also similar.

Figure 151. Graph. 3D-Move predictions and measured deflections for cell 34 for RWD trials (device velocity = 30 mi/h (48 km/h)). This graph illustrates the comparison of 3D-Move predicted (3D-Move/case 1 and 3D-Move/case X1) and measured deflections (GEO 1/GEO 3/ lower bound and GEO 1/ GEO 3/upper bound) and Rolling Wheel Deflectometer (RWD) measurements (RWD/D0/upper, RWD/D0/bottom, RWD/D15/upper, and RWD/D15/bottom) for cell 34 based on testing done by the RWD at a vehicle velocity of 30 mi/h (48 km/h). The y-axis shows deflection from -35 to 5 mil (-0.889 to 0.127 mm), and the x-axis shows distance from -10 to 10 ft (-3.05 to 3.05 m). In all cases, 3D-Move adequately captures the bell shape of measured displacements. Maximum displacements from 3D-Move/case 1 and 3D-Move/case X1, which are 22 and 30 mil (0.56 to 0.762 mm), respectively, cover the ranges of measured displacements (i.e., GEO 1/GEO 3/lower bound and GEO 1/GEO 3/upper bound) which are 22 and 25 mil (0.559 to 0.635 mm), respectively. RWD/D0 is about 25 mil (0.635 mm), and RWD/D15 is 15 mil (0.381) for both upper and bottom measurements.

1 mil = 0.0254 mm

1 ft = 0.305 m

Figure 151. Graph. 3D-Move predictions and measured deflections for cell 34 for RWD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 152. Graph. Normalized basins of 3D-Move predictions and measured deflections for cell 34 in RWD trials (device velocity = 30 mi/h (48 km/h)). This graph shows normalized basins of 3D-Move predicted (3D-Move/case 1 and 3D-Move/case X1) and measured deflections (GEO 1/GEO 3/upper bound and GEO 1/GEO 3/lower bound) for cell 34 based on testing done by the Rolling Wheel Deflectometer (RWD) at a vehicle speed of 30 mi/h (48 km/h). The y-axis shows normalized deflection from -1.2 to 0.2 mil/mil (-1.2 to 0.2 mm/mm), and the x-axis shows distance from -10 to 10 ft (-3.05 to 3.05 m). The predicted and measured deflection shapes agree closely. This graph reveals that the differences in shape among the computed and measured deflection basins are small at the center but vary slightly as distance increases from the center.

1 mil/mil = 1 mm/mm

1 ft = 0.305 m

Figure 152. Graph. Normalized basins of 3D-Move predictions and measured deflections for cell 34 in RWD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 153. Graph. Comparison of normalized deflection basins at various levels for cell 34 in RWD trials (device velocity = 30 mi/h (48 km/h)). This bar graph shows the comparison of the pulse widths of normalized basins of 3D-Move predicted (case 1 and case X1) and measured deflections (GEO 1 and GEO 3) at three levels of maximum displacement percentage for cell 34 based on testing done by the Rolling Wheel Deflectometer (RWD) at a vehicle speed of 30 mi/h (48 km/h). The y-axis shows pulse width from 0 to 5 ft (0 to 1.52 m), and the x-axis shows percentage of maximum displacement, which are (from left to right) 75, 50, and 25 percent. This graph shows a similar trend at all three levels, and the 3D-Move (cases 1 and X1) produce close normalized deflection basins compared to the measured ones (GEO 1 and GEO 3). Maximum measured pulse width varies from 1.4 to 4 ft (0.427 to 1.22 m), and maximum computed pulse width varies from 1.4 to 4.3 ft (0.427 to 1.31 m). Measured and computed Pulse widths at 75 percent of maximum displacement are closer than those at 50 and 25 percent of maximum displacement.

1 ft = 0.305 m

Figure 153. Graph. Comparison of normalized deflection basins at various levels for cell 34 in RWD trials (device velocity = 30 mi/h (48.3 km/h)).

Summary

The 3D-Move analyses were conducted using inputs derived based on the following considerations:

Several 3D-Move analyses were undertaken in an attempt to bracket the measured surface deflection data (peak and basin shape). In the process, due consideration was given to the selection of 3D-Move inputs (material properties and loading) so that they were rationally determined and the adjustments were defensible.

It was shown that 3D-Move results fell within the measured time histories from project geophones. Figure 154 shows the comparison of maximum displacements computed by 3D-Move and those measured by the project sensors for all accuracy runs made with the TSD and RWD. In this comparison, case X1 was used for cells 3 and 34, and case X was used for cell 19. When project sensor measurements were plotted, the largest displacement given by either GEO 1 or GEO 3 (from all field trials) was selected. The figure shows a good match between computed and measured maximum displacements.

Figure 154. Graph. 3D-Move computed maximum displacement versus measured displacement for all cells and vehicle velocities during RWD and TSD trials. This graph shows the comparison of maximum displacements computed by 3D-Move and those measured by the project sensors for all accuracy runs made with the Traffic Speed Deflectometer (TSD) and Rolling Wheel Deflectometer (RWD). The y-axis shows computed maximum displacement/3D-Move from 0 to 40 mil (0 to 1.02 mm), and the x-axis shows measured maximum displacement/project sensors from 0 to 40 mil (0 to 1.02 mm). The graph shows a good match between computed and measured maximum displacements. The constant is -1.47 mil (-0.037 mm), the slope is 1.1087, the R square value is 0.939, and the standard error of estimate is 2.26 mil (0.057 mm). The measured maximum displacement corresponding to three cells are grouped with in the box with cell 3 values ranging between 10 and 15 mil (0.254 and 0.381 mm), cell 19 ranging between 17 and 22 mil (0.43 and 0.56 mm), and cell 34 ranging between 24 and 36 mil (0.61 and 0.91 mm).

1 mil = 0.0254 mm

Figure 154. Graph. 3D-Move computed maximum displacement versus measured displacement for all cells and vehicle velocities during RWD and TSD trials.

In addition to maximum displacement, it was also important to have a good representation of the deflection bowls since TSDDs use deflections at various points of a bowl to predict pavement layer condition. Figure 155 shows the comparison of the shapes of normalized deflection bowls that gave the maximum deflection with 3D-Move analyses and sensor measurements (GEO 1 and GEO 3). Three levels of the normalized deflection bowls (75, 50, and 25 percent of maximum displacement) are reported. In this figure, wmax is the normalized maximum deflection. Based on this information, it was concluded that 3D-Move can predict the displacement bowl similar to those measured by geophones.

Figure 155. Graph. 3D-Move computed pulse width versus measured pulse for all cells and velocities in RWD and TSD trials. This graph shows the comparison of the shapes of normalized deflection bowls that gave the maximum deflection with 3D-Move analyses and sensor measurements for all cells and velocities in the Rolling Wheel Deflectometer (RWD) and Traffic Speed Deflectometer (TSD) trials. The y-axis shows computed pulse width for 3D-Move from 0.2 to 5.2 ft (0.061 to 1.58 m), and the x-axis shows the measured pulse width for the project sensors from 0.2 to 6.2 ft (0.061 to 1.9 m). Three levels of the normalized deflection bowls (75, 50, and 25 percent of maximum displacement, W subscript max) are shown. 3D-Move can capture the displacement bowl well. The constant is -0.2951 ft (0.09 m), the slope is 0.9843, the R square value is 0.8734, and the standard error of estimate is 0.43 ft (0.13 m). The measured pulse width corresponding to three cells are grouped with in the box with cell 3 values ranging between 1.2 and 1.6 ft (0.366 and 0.488 m), cell 19 ranging between 2.3 and 3 ft (0.7 and 0.91 m), and cell 34 ranging between 3.1 and 5.1 ft (0.945 and 1.55 mm).

1 ft = 0.305 m

Figure 155. Graph. 3D-Move computed pulse width versus measured pulse for all cells and velocities in RWD and TSD trials.

7.6 3D-Move Calibration using MnROAD Sensor Data

MnROAD Sensors

This subsection provides a comparison between the 3D-Move results and the measured MnROAD sensor data. Unlike the earlier sensor comparisons, which focused on surface deflections, the comparisons presented in this subsection consider stresses and strains at various interior pavement locations. Since the load-induced stresses and strains are critical inputs to pavement performance predictions, this effort is important in the validation of the applicability of the 3D-Move for pavement response predictions to be used in identifying the most promising indices from TSDD measurements that best relate to pavement structure.

The existing instrumentations at the MnROAD facility include normal strain responses in the longitudinal direction and in some cases in the transverse direction as well as vertical pressure histories in unbound base and subgrade layers. The measurements from these sensors can be directly compared with those computed by 3D-Move. However, it is important to recognize that measurements made by these sensors reflect an average value over the entire dimensions of the sensors (e.g., pressure cell measurements reflect an average Earth pressure over a diameter of 9 inches (228.6 mm)). Table 41 presents the sensors that were monitored during the accuracy runs.

Table 41. Available sensor data in cells 3, 19, and 34 (used in 3D-Move calibration).
Sensor ID Cell Number Facility Sensor Type Latitude (Degree)  Longitude (Degree) Depth (inches)
03LE101 3 ML Longitudinal SG 45.27331792 -93.73110608 9.3
03LE102 3 ML Longitudinal SG 45.27332065 -93.73111185 9.4
03LE103 3 ML Longitudinal SG 45.27332372 -93.73111822 9
03PG103 3 ML Pressure cell 45.27329061 -93.73105228 14
19LE101 19 ML Longitudinal SG 45.25624056 -93.69966629 4.9
19LE104 19 ML Longitudinal SG 45.2562641 -93.69970936 4.7
19PG101 19 ML Pressure cell 45.25621817 -93.69962245 18
19PG103 19 ML Pressure cell 45.25625157 -93.69968785 18.2
34LE202 34 LVR Longitudinal SG N/A N/A 3.5
34LE203 34 LVR Longitudinal SG N/A N/A 3.5
34PG201 34 LVR Pressure cell N/A N/A 15.9
34PG202 34 LVR Pressure cell N/A N/A 15.9
34PG203 34 LVR Pressure cell N/A N/A 15.9

1 inch = 25.4 mm

ML = Main line.

SG = Strain gauge.

N/A = Not available.

3D-Move Results for Accuracy Runs

The MnROAD pressure cells measured the average dynamic pressure applied to the entire plate, while the 3D-Move responses were computed point by point (i.e., at any given single location). Therefore, an average estimate of the computed pressure over the entire plate diameter was used for comparison purposes. Since the location of the wheel during the field trials were not known precisely, it was decided to compute 3D-Move responses at many transverse locations as an attempt to bracket the possible variation in the responses as a result of wheel wander. Figure 156 shows the response points considered along the transverse direction in 3D-Move runs. These points are located at 2-inch (50.8-mm) intervals with enough coverage so that the role of the wheel wander can be investigated. The difference between the average pressure on the plate and pressure at the center of the plate was less than 5 percent.

Figure 156. Illustration. Response points for normal pressure in 3D-Move. This illustration shows the response points considered along the transverse direction in 3D-Move runs. The figure includes the inner and outer wheels. The first response point is located on the midline between the tires, and the other points are located at 2-inch (50.8-mm) intervals with enough coverage so that the role of the wheel wander can be investigated.

1 inch = 25.4 mm

Figure 156. Illustration. Response points for normal pressure in 3D-Move.

The average MnROAD strain responses were computed over the length of the strain gauges, which were 6 inches (152.4 mm) long. Figure 157 shows the locations used in 3D-Move to determine the role of wheel wander in the comparison of strain responses.

Figure 157. Illustration. Response points for strain gauges in 3D-Move. This illustration shows the locations used in 3D-Move to determine the role of wheel wander in the comparison of strain responses. It includes the inner and outer wheels. The first response point is located on the midline between tires, and the other points are located at 2-inch (50.8-mm) intervals with enough coverage so that the role of the wheel wander can be investigated.

1 inch = 25.4 mm

Figure 157. Illustration. Response points for strain gauges in 3D-Move.

In the 3D-Move versus project sensor comparison, the following same three scenarios for analyses were found to be appropriate:

These cases were used to compare the MnROAD sensor measurements with the 3D-Move responses.

It became apparent from the start that the earth pressures predicted by 3D-Move were substantially greater than those measured by the MnROAD sensors. The differences were larger than what could be realistically attributed to the wheel wander. To reconcile this difference, other sources where earth pressures had been predicted and measured were reviewed. As a first step, the static LE software WESLEA results were evaluated for the TSDDs loading. Unlike 3D-Move, WESLEA cannot provide the response time history nor can it provide the results as a function of vehicle speed. The results of WESLEA were compared only with the responses from the slowest TSDD speed, which was 30 mi/h (48.3 km/h) for both the RWD and TSD. The 3D-Move and WESLEA pressures were compared mainly to determine whether the measured earth pressures were reasonable.

Table 39, which was presented earlier, contains the AC moduli used in the WESLEA analyses. The existing AC moduli at a frequency of 30 Hz were adjusted to the temperature at the time of the TSDDs trials. The FWD backcalculated values were directly used for unbound layers.

Table 42 and figure 158 show the 3D-Move and WESLEA computed vertical normal pressures along various transverse offsets for the TSD pass along cell 34. The values for each case are close to one another. Minor differences could be attributed to issues related to the viscoelastic characterization of the AC layer and to vehicle speed in 3D-Move. As seen in figure 158, the values measured by the MnROAD sensors were significantly lower (a factor of about 4 to 6). Conversely, the normalized predicted and measured earth pressures shown in figure 159 agree closely, leading to questioning of the correctness of the measurements. These sensors were installed during construction about 5 years before the TSDDs trials, and hence it is possible that the calibration was off or they are damaged. It should also be noted that the earth pressures measured in cell 19 in both TSDDs had substantial negative values, which were unrealistic. This information, coupled with the close match with 3D-Move comparisons with WESLEA, led to the conclusion that the MnROAD earth pressure sensor data were suspect.

Table 42. 3D-Move versus WESLEA predictions at various offsets for cell 34 during TSD trials (device velocity = 30 mi/h (48.3 km/h)).
Offset Distance (inches) Pressure From 3D-Move Case 1 (psi) Pressure From WESLEA Analyses Case 1 (psi) Pressure From 3D-Move Case X1 (psi) Pressure from WESLEA Analyses Case X1 (psi)
0 7.7 7.2 10.6 10.2
2 7.5 7.2 10.4 10.1
4 7.2 7 10 9.9
6 6.8 6.8 9.4 9.5
8 6.3 6.4 8.6 8.9
10 5.7 5.9 7.8 8.2
12 5 5.4 6.8 7.4
14 4.4 4.8 5.9 6.6
16 3.7 4.2 5 5.7
18 3.2 3.6 4.2 4.9
20 2.7 3.1 3.5 4.1
22 2.2 2.6 2.8 3.5
24 1.9 2.2 2.3 2.9

1 inch = 25.4 mm

1 psi = 6.89 kPa

Figure 158. Graph. 3D-Move versus WESLEA predictions and MnROAD pressure cell for cell 34 in TSD trials (device velocity = 30 mi/h). This scatter plot shows the 3D-Move (case 1 and case X1), Waterways Experiment Station Linear Elastic Analysis (WESLEA) (case 1 and 
case X1), and MnROAD (pressure gauge (PG) 1/max pressure, PG2/max pressure, and PG3/max pressure) computed vertical normal pressures along various transverse offsets for the Traffic Speed Deflectometer (TSD) pass along cell 34 at a vehicle speed of 30 mi/h (48.3 km/h). The y-axis shows maximum of pressure from 0 to 12 psi (0 to 82.7 kPa), and the x-axis shows distance of center pressure cell to the midline between tires from 0 to 30 inches (0 to 762 mm). The 3D-Move and WESLEA predicted maximum pressure values for each case are relatively close to one another. The maximum pressure ranges from 2 to 7.5 psi (13.8 to 51.7 kPa) for case 1 and from 3 to 10.5 psi (20.7 to 72.3 kPa) for case X1 for both 3D-Move and WESLEA. The average maximum pressure value measured by the MnROAD sensors is about 1 psi (6.89 kPa), which is lower than the 3D-Move and WESLEA predicted pressure values.

1 psi = 6.89 kPa

1 inch = 25.4 mm

Figure 158. Graph. 3D-Move versus WESLEA predictions and MnROAD pressure cell for cell 34 in TSD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 159. Graph. Normalized basins of 3D-Move predictions and MnROAD pressure cell measurement for cell 34 in TSD trials (device velocity = 30 mi/h (48.3 km/h)). This scatter plot shows the normalized 3D-Move predicted (3D-Move/case 1 and 3D-Move/case X1) and MnROAD measured maximum pressure in cell 34 for the Traffic Speed Deflectometer (TSD) trials at a vehicle speed of 30 mi/h (48.3 km/h). The y-axis shows normalized pressure from -0.2 to 1.2 psi/psi (-0.2 to 1.2 kPa/kPa), and the x-axis shows distance from -10 to 10 ft (-3.05 to 3.05 m). The normalized predicted and measured pressures agree closely for all cases.

1 psi/psi = 1 kPa/kPa

1 ft = 0.305 m

Figure 159. Graph. Normalized basins of 3D-Move predictions and MnROAD pressure cell measurement for cell 34 in TSD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 160 compares the measured and predicted longitudinal strains in cell 34 for the TSD runs. Unlike the case of the earth pressures, the calculated longitudinal strains from 3D-Move matched well with the measured data from the MnROAD strain gauges. Table 43 and figure 161 through figure 163 show the comparison of the measured and calculated vertical pressure, normalized vertical pressure, and longitudinal strains for the RWD runs for cell 3. Similar to the TSD findings, the values for longitudinal strain were predicted closely by 3D-Move.

Similar plots were developed for all possible combinations of TSDD speeds and accuracy cells, and they showed the same trends as shown in the example figures included with in this chapter. The data used to generate all of these figures/tables may be obtained directly from FHWA.

Figure 160. Graph. 3D-Move versus MnROAD strain gauge measurement for cell 34 in TSD trials (device velocity = 30 mi/h (48.3 km/h)). This scatter plot compares the MnROAD measured (strain gauge (SG)/Max) and 3D-Move predicted (3D-Move/case 1 and 3D-Move/ case X1) longitudinal strains in cell 34 for the Traffic Speed Deflectometer (TSD) runs at a vehicle speed of 30 mi/h (48.3 km/h). The y-axis shows longitudinal strain from -200 to 500 microstrain, and the x-axis shows distance from -10 to 10 ft (-3.05 to 3.05 m). The calculated longitudinal strains from 3D-Move match well with the measured data from the MnROAD strain gauges. All cases follow a similar shape, and 3D-Move/case X1 provides the maximum longitudinal strain value closest to the measured longitudinal strain. In all curves, strains are equal to 0 microstrain until -4 ft (-1.22 m) and then decrease to minimum values which are -100 microstrain for the 3D-Move/case X1 and MnROAD sensors and -50 microstrain for 3D-Move/case 1. They then sharply increase to maximum values at 0 ft (0 m) and 460 microstrain for the 3D-Move/case X1 and MnROAD sensors and 320 microstrain for 3D-Move/case 1. They sharply decrease to minimum values, which are -100 microstrain for 3D-Move/case X1, -50 microstrain for 3D-Move/case 1, and -20 microstrain for the MnROAD sensor. Finally they all increase to 0 microstrain at 4 ft (1.22 m).

1 ft = 0.305 m

Figure 160. Graph. 3D-Move versus MnROAD strain gauge measurement for cell 34 in TSD trials (device velocity = 30 mi/h (48.3 km/h)).

Table 43. 3D-Move versus WESLEA predictions in various offsets for cell 3 RWD trials (device velocity = 30 mi/h (48.3 km/h)).
Offset Distance (inches) Pressure From 3D-Move Case 1 (psi) Pressure From WESLEA Analyses Case 1 (psi) Pressure From 3D-Move Case X1 (psi) Pressure From WESLEA Analyses Case X1 (psi)
0 10.5 10.7 12.9 13.2
2 10.4 10.7 12.8 13.2
4 10.2 10.7 12.5 13.1
6 9.7 10.4 11.8 12.7
8 8.8 9.8 10.7 12
10 7.7 8.9 9.2 10.8
12 6.3 7.6 7.6 9.4
14 5 6.2 6.1 7.7
16 3.8 4.9 4.6 6.1
18 2.8 3.7 3.4 4.7
20 2 2.7 2.5 3.5
22 1.4 2 1.8 2.6
24 1 1.4 1.3 1.8

1 inch = 25.4 mm

1 psi = 6.89 kPa

sensor. Finally they all increase to 0 microstrain at 4 ft (1.22 m). Figure 161. Graph. 3D-Move versus WESLEA predictions and MnROAD pressure cell for cell 3 in RWD trials (device velocity = 30 mi/h (48.3 km/h)). This scatter plot shows the 3D-Move (case 1 and case X1), Waterways Experiment Station Linear Elastic Analysis (WESLEA) (case 1 and case X1), and MnROAD (pressure gauge (PG) 1/max pressure) computed vertical normal pressures along various transverse offsets for the Rolling Wheel Deflectometer (RWD) pass along cell 3 at a vehicle speed of 30 mi/h (48.3 km/h). The y-axis shows maximum of pressure from 0 to 14 psi (0 to 96.5 kPa), and the x-axis shows distance of center of pressure cell to midline between tires form 0 to 30 inches (0 to 762 mm). The 3D-Move and WESLEA predicted maximum pressure values for each case are relatively close to one another, where the maximum pressure ranges from 1.5 to 10.5 psi (10.3 to 72.35 kPa) for case 1 and from 2 to 13 psi for case X1. The average maximum pressure value measured by the MnROAD sensor is about 2 psi (13.78 kPa), which is lower than the 3D-Move and WESLEA predicted pressure values.

1 psi = 6.89 kPa

1 inch = 25.4 mm

Figure 161. Graph. 3D-Move versus WESLEA predictions and MnROAD pressure cell for cell 3 in RWD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 162. Graph. Normalized basins of 3D-Move predictions and MnROAD pressure cell measurement for cell 3 in RWD trials (device velocity = 30 mi/h (48.3 km/h)). This scatter plot shows the normalized 3D-Move predicted (3D-Move/case 1 and 3D-Move/case X1) and MnROAD measured maximum pressure in cell 3 for the Rolling Wheel Deflectometer (RWD) trials at a vehicle speed of 30 mi/h (48.3 km/h). The y-axis shows normalized pressure from -0.2 to 1.2 psi/psi (-0.2 to 1.2 kPa/kPa), and the x-axis shows distance from -10 to 10 ft (-3.05 to 3.05 m). The normalized predicted and measured pressures agree closely for all lines.

1 psi/psi = 1 kPa/kPa

1 ft = 0.305 m

Figure 162. Graph. Normalized basins of 3D-Move predictions and MnROAD pressure cell measurement for cell 3 in RWD trials (device velocity = 30 mi/h (48.3 km/h)).

Figure 163. Graph. 3D-Move versus MnROAD strain gauge measurement for cell 3 in RWD trials (device velocity = 30 mi/h (48.3 km/h)). This scatter plot compares the MnROAD measured strain gauge (SG)/Max and 3D-Move predicted (3D-Move/case 1 and 3D-Move/ case X1) longitudinal strains in cell 3 for the Rolling Wheel Deflectometer (RWD) runs at a vehicle speed of 30 mi/h (48.3 km/h). The y-axis shows longitudinal strain from -100 to 200 microstrain, and the x-axis shows distance from -10 to 10 ft (-3.05 to 3.05 m). The calculated longitudinal strains from 3D-Move/case 1 and 3D-Move/ case X1 match well with the measured data from the MnROAD strain gauge. All cases follow a similar shape. In all curves, strains are equal to 0 microstrain until -6 ft (-1.83 m) and then decrease to minimum values which are -25 microstrain for 3D-Move/case 1 and the MnROAD sensor and -45 microstrain for 3D-Move/case X1. Values then sharply increase to maximum values at 0 ft, (0 m) which are 150 microstrain for 3D-Move/case X1, 130 microstrain for the MnROAD sensor, and 
100 microstrain for 3D-Move/case 1. Values then sharply decrease to minimum values, which are -45 microstrain for 3D-Move/case X1, and -25 microstrain for 3D-Move/case 1, and -20 microstrain for MnROAD sensor. Values then increase to 0 microstrain at 6 ft (1.83 m) for all lines.

1 ft = 0.305 m

Figure 163. Graph. 3D-Move versus MnROAD strain gauge measurement from cell 3 in RWD trials (device velocity = 30 mi/h (48.3 km/h)).

Summary

The vertical pressures measured by MnROAD sensors were significantly less by a factor of 4 to 6 than those computed by 3D-Move. The possible transverse wheel wander could not explain the differences. However, the shapes of the pressure histories, when normalized, closely matched with those given by 3D-Move. Figure 164 shows the comparison of 3D-Move prediction of normal pressures and those computed from WESLEA for all accuracy passes and cells. In this graph, only the 3D-Move results computed with the lowest vehicle speed are shown. The 3D-Move normal pressure values were consistent and reasonable with WESLEA predictions, unlike those measured with MnROAD sensors.

The maximum longitudinal strains from the MnROAD measurements and 3D-Move predictions were similar (see figure 165). Accordingly, it can be concluded that 3D-Move is capable of predicting field measured surface displacement histories and interior pavement responses (stresses and strains) and could be used to evaluate pavement responses under TSDD loadings.

Figure 164. Graph. Computed normal pressure for 3D-Move versus WESLEA. This graph shows a comparison of 3D-Move prediction of normal pressures and those computed from Waterways Experiment Station Linear Elastic Analysis (WESLEA) for all accuracy passes and cells. The y-axis is labeled WESLEA prediction of normal pressure at various distances from 0 to 18 psi (0 to 124 kPa), and the x-axis shows 3D-Move prediction of normal pressure at various distances from 0 to 18 psi (0 to 124 kPa). The 3D-Move normal pressure values are consistent and reasonable with WESLEA values. The slope of the best fitted line is y equals 1.0034 times x, and the R square value is 0.994.

1 psi = 6.89 kPa

Figure 164. Graph. Computed normal pressure for 3D-Move versus WESLEA.

Figure 165. Graph. Maximum longitudinal strains from MnROAD sensors and 3D-Move computations. This bar graph shows a comparison of computed (3D-Move/upper bound and 3D-Move/lower bound) and measured MNROAD longitudinal strains in Traffic Speed Deflectometer (TSD) and Rolling Wheel Deflectometer (RWD) passes. The y-axis shows longitudinal strain from 0 to 500 microstrain, and the x-axis shows TSD and RWD trials in the accuracy cells at various speeds. The maximum longitudinal strains from the MnROAD measurements and 3D-Move/upper bound predictions agree well. The minimum and maximum measured longitudinal strain are 80 and 460 microstrain, and the minimum and maximum computed longitudinal strain from 3D-Move/upper bound are 145 and 460 microstrain.

Figure 165. Graph. Maximum longitudinal strains from MnROAD sensors and 3D-Move computations.

 

 

 

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