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
Publication Number: FHWA-RD-02-088
Date: May 2003

Evaluation of Joint and Crack Load Transfer Final Report

Previous | Table of Contents | Next

Load Transfer Efficiency

When a traffic load is applied near a joint in a PCC pavement, both loaded and unloaded slabs deflect because a portion of the load applied to the loaded slab is transferred to the unloaded slab.  As a result, deflections and stresses in the loaded slab may be significantly less than if, instead of a joint with another slab, there was a free edge.  The magnitude of reduction in stresses and deflections by a joint compared to a free edge depends on the joint's LTE.

Traditionally, LTE at the joint is determined based on the ratio of the maximum deflection at the joint of the loaded slab and the deflection of the unloaded slab measured right across the joint from the maximum deflection.  Two equations for the deflection LTE are used most often:

Equation 1 Load transfer efficiency index equals the maximum deflection at the joint of the loaded slab divided by the corresponding deflection at the joint of the unloaded slab times 100 percent. (1)

and

Equation 2 Load transfer efficiency index equals 2 times the corresponding deflection at the joint of the unloaded slab divided by the sum total of the maximum deflection at the joint of the loaded slab plus the corresponding deflection at the joint of the unloaded slab, times 100 percent. (2)

where:

-          dl = the maximum deflection at the joint of the loaded slab.

-          du = the corresponding deflection at the joint of the unloaded slab.

-          LTE and LTE* = load transfer efficiency indexes. 

If a joint exhibits a poor ability to transfer load, then the deflection of the unloaded slab is much less than the deflection at the joint of the loaded slab and both LTE indexes have values close to 0.  If joint's load transfer ability is very good, then the deflections at the both sides of the joint are equal and both indexes have values close to 100 percent.  Moreover, these two indices are related by the following equation:

Equation 3 Load transfer efficiency index equals 2 times the sum total of 1 minus the sum total of 1 divided by 1 plus the load transfer efficiency divided by 100, times 100 percent. (3)

Therefore, these indexes are equivalent, and if one of them is known, the other can be determined.  In this study, we will define the deflection LTE using the index from equation 1 because it is much more widely used.

The definitions of the LTE described above are based on deflections.  A joint LTE based on stress can be defined as follows:

Equation 4 Load transfer efficiency in stress equals the corresponding stress at the joint of the unloaded slab divided by the maximum stress at the joint of the loaded slab times 100 percent. (4)

where:

-          the maximum stress at the joint of the loaded slab = the maximum stress at the joint of the loaded slab.

-          the corresponding stress at the joint of the unloaded slab = the corresponding stress at the joint of the unloaded slab.

-          LTEsigma subscript = load transfer efficiency in stress.

The stress-based LTE indicates the degree of stress reduction at the joint of the loaded slab caused by the presence of the unloaded slab.  Studies have indicated that there is no one-to-one relationship between stress-based and deflection-based LTE indexes.  Because it is difficult to measure stresses in a concrete slab and stress-based LTE is much more affected by geometry of the applied load than deflection LTE, the deflection-based LTE is commonly used to measure load transfer in concrete pavements.  If deflection LTE is known, stress reduction due to load transfer for any given load configuration may be calculated using a finite element model.

Joint LTE depends on many factors, including the following:

Load transfer between the slabs occurs through aggregate particles of the fractured surface below the saw cut at a joint, through steel dowels (if they exist), and through the base and subgrade.  LTE may vary throughout the day and year because of variation in PCC temperature.  When temperature decreases, a joint opens wider, which decreases contact between two slabs and also may decrease LTE, especially if no dowels exist.  Also, PCC slab curling may change the contact between the slab and the underlying layer and affect measured load-induced deflections.

Mechanistic Modeling of LTE

Mechanistic modeling of the load transfer mechanism is a complex problem.  Frieberg (1940) was one of the first researchers who attempted to tackle this problem.  The introduction of a finite element method for analysis of JCP (Tabatabie and Barenberg [1980]) gave a significant boost to understanding load transfer mechanisms.  However, although many comprehensive finite element models have been developed (Scarpas et al. [1994], Guo et al. [1995], Parson et al. [1997], Shoukry and William [1998], Davids et al. [1998], Khazanovich et al. [2001]), work on the development of a comprehensive, practical, and reliable model for joints of rigid pavements is far from complete. 

Ioannides and Korovesis (1990, 1992) made an important step forward in PCC joint analysis. They identified the following nondimensional parameters governing joint behavior:

Nondoweled pavements:

 Equation 5 Nondimensional stiffness of nondoweled joints equals shear stiffness of a unit length of an aggregate interlock divided by subgrade K-value times portland cement concrete slab radius of relative stiffness. (5)

Doweled pavements:

Equation 6 Nondimensional stiffness of doweled joint equals shear stiffness of a single dowel divided by dowel spacing times subgrade K-value times portland cement concrete slab radius of relative stiffness. (6)

where:

-          AGG* and D* are nondimensional stiffnesses of nondoweled and doweled joints.

-          AGG = a shear stiffness of a unit length of an aggregate interlock.

-          D = a shear stiffness of a single dowel (including dowel-PCC interaction).

-          Radius of relative stiffness = a PCC slab radius of relative stiffness.

-          k = subgrade k-value.

-          s = dowel spacing.

Using the finite element program ILLI-SLAB, Ioannides and Korovesis also identified a unique relationship between these parameters and LTE (see figure 1).

Figure 1. Load transfer efficiency versus nondimensional joint stiffness. The figure is a line graph. The nondimensional joint stiffness (AGG*), of 0.01 to 1000, is graphed on the horizontal axis. The load transfer efficiency, percent, is graphed on the vertical axis. The graph starts at the origin of point (0, 0.01) and increases in an S-shaped curve to the approximate points (98, 500). As the load transfer efficiency increases, so does nondimensional joint stiffness.

Figure 1. LTE versus nondimensional joint stiffness.

The following assumptions were made in derivation of these relationships:

The relationships developed by Ioannides and Korovesis form a basis for back calculation of joint aggregate interlock stiffness of nondoweled joints or dowel shear stiffness of doweled joints if their LTEs are known.  In both cases, however, the backcalculated stiffness overestimates real aggregate interlock stiffness or dowel stiffness because the entire joint stiffness is attributed to a single (although perhaps prevailing) component.  The cause of this limitation is the inability of the ILLI-SLAB model to distinguish between load transfer mechanisms.  This limitation does not cause a significant problem because the backcalculated joint stiffness provides sufficient information for accurate joint modeling in the forward analysis.  Moreover, if the addition information is available, individual components can be determined more accurately.

For example, if AGGtot is the backcalculated stiffness of a doweled joint and D is the known dowel stiffness, then the "true" aggregate interlock factor, AGG0, for this joint can be determined from the following relationship:

Equation 7 The true aggregate interlock factor equals the backcalculated stiffness of a doweled joint minus the sum total of the known dowel stiffness times dowel spacing. (7)

where:

-          s = dowel spacing.

In this study, Ioannides and Korovesis's relationship was further investigated and an efficient backcalculation procedure for joint stiffness determination was developed.

Joint Stiffness Backcalculation

The relationship identified by Ioannides and Korovesis was further elaborated by Crovetti (1994), Zollinger et al. (1999), and Ioannides et al. (1996).  Crovetti and Zollinger developed regression models for that relationship, whereas Ioannides and Hammons developed a neural network prediction model. 

Crovetti proposed the following relationship between nondimensional joint stiffness and LTE:

Equation 8 Load transfer efficiency index equals 100 percent divided by the sum total of 1 plus the total of 1.2 times the total of the total stiffness divided by the coefficient of subgrade reaction (K-value) times the portland cement concrete slab radius of relative stiffness, all to the -0.849 power.                                                           (8)

where:

-          AGGtot = the total stiffness.

-          Radius of relative stiffness = the PCC slab radius of relative stiffness.

-          k = a coefficient of subgrade reaction (k-value).

Zollinger's model for this relationship has the following form:

Equation 9 Zollinger's model.  Load transfer efficiency index equals 100 percent divided by the sum total of 1 plus the total of 10 multiplied by the sum total of 0.214 minus 0.183 times the equivalent radius of the applied load divided by the portland cement concrete slab radius of relative stiffness minus the log rhythm of the total joint stiffness divided by the subgrade K-value times the portland cement concrete slab radius of relative stiffness, divided by 1.18.                                                       (9)

where

-          AGGtot = the total joint stiffness.

-          Radius of relative stiffness = the PCC slab radius of relative stiffness.

-          k = the subgrade k-value.

-          a = the equivalent radius of the applied load.

The neural network model developed by Ioannides and Hammons cannot be expressed as a simple equation, but rather takes the form of a computer program that relates LTE and joint stiffness, radius of relative stiffness, subgrade k-value, and load geometry.  Although this model is more accurate than the regression models, it is not publicly available and was not evaluated in this study.

To evaluate the Crovetti and Zollinger models, a factorial of 375 runs was performed to simulate FWD testing at the PCC joint.  ISLAB2000 (Khazanovich et al. [2000]) is a completely rewritten version of ILLI-SLAB that retains all the positive features of ILLI-SLAB but is more convenient to use and is free from many unnecessary limitations (including limitation on the number of nodes in a finite element model).

A finite element model developed in this study has four slabs in the longitudinal direction and three slabs in the transverse direction.  The system was assumed to be symmetrical with respect to the longitudinal axis, so only half of the system was modeled.  Because the focus of this analysis was on the deflections near the loaded joint along the centerline, a much finer finite element mesh was used along the centerline and loaded joint, as shown in figure 2.  The LTE of longitudinal joints was selected to be equal to 70 percent.  Also, due to symmetry, only half of an FWD load was applied.  The slab radius of relative stiffness was varied from 508 to

2,032 mm (80 inches), and the nondimensional transverse joint stiffness was varied from 0.1 to 278, which resulted in LTE from 8 to 99 percent.

Figure 2. Finite element model for FWD loading simulation. The figure is a rectangle divided into four squares. In each square is a pattern of grids. The falling weight deflectometer load resides on the lower part of the rectangle and slightly to the left of the center horizontal line. The figure shows that the system is not symmetrical.

Figure 2.  Finite element model for FWD loading simulation.

Figure 3 presents comparisons of the LTEs calculated from the ISLAB2000 results and LTEs from Crovetti's and Zollinger's equations.  Although Crovetti's equation is simpler than Zollinger's, it better corresponds with ISLAB2000 LTEs.  Because the Crovetti model was adopted for use in the 2002 design procedure (National Cooperative Highway Research Program Project 1-37A) and it compares well with the finite element analysis, it was selected for use in this study.

Figure 3. Comparison of LTE calculated from ISLAB2000 results with predictions using Crovetti's and Zollinger's models. The ISLAB 2000 load transfer efficiency, percent, is graphed on the horizontal axis. The predicted load transfer efficiency, percent, is graphed on the vertical axis. A straight line starts at the origin, increasing in a 45-degree angle ends at the point (100, 100). The figure has two sites, which are Crovetti's and Zollinger's equation. Both sites are scattered close together along this line.

Figure 3.  Comparison of LTE calculated from ISLAB2000 results with predictions using Crovetti's and Zollinger's models.

"True" versus "Measured" LTE

As stated above, the joint LTE is a ratio of the maximum deflection at the joint of the unloaded slab to the deflection of the loaded slab measured directly across the joint from the maximum deflection.  However, measurement of such deflections in the field may be quite cumbersome.  In the LTPP program, joint deflection testing is conducted by placing the load plate tangential to the edge of the joint.  The loaded slab joint deflections are measured under the center of the load plate (152 mm [6 inches] away from the joint).  The deflection of the unloaded slab is also measured at some distance (152 mm [6 inches]) from the joints.  This raises an issue about the possibility of error as a result of differences in deflections directly at the joint and measured deflections 152 mm (6 inches) away from the joint because of slab bending.  Some experts advocate the need to use a correction ("bending") factor to adjust the measured deflections.

To investigate this problem, the results of the 375 ISLAB2000 runs were analyzed.  The LTEs calculated from the deflections located exactly at the joints ("true" LTEs) were compared to the ratios of the deflections located 152 mm (6 inches) away from the joint ("measured" LTEs similar to the FWD procedure).  Figure 4 presents comparisons of true and measured LTEs.

Figure 4. Comparison of true and measured LTEs. True load transfer efficiency, percent, is graphed on the horizontal axis. Measured load transfer efficiency, percent, is graphed on the vertical axis. The figure is a line graph increasing in a straight line at a 45-degree angle. There are plots scattered closely above the slope. Measured load transfer efficiency is close to true load transfer efficiency.

Figure 4.  Comparison of true and measured LTEs.

In most cases, measured LTE is close to true LTE.  The exceptions are the cases with very low PCC slab stiffness (radius of relative stiffness is less than 750 mm (30 inches).  For those cases, measured LTE overestimates true LTE.  To address this discrepancy, the research team attempted to develop bending correction factors. 

PCC slab bending depends primarily on the radius of relative stiffness.  Since stiffer slabs require less bending correction than slabs with low radii of relative stiffness, the following functional form was proposed for the correction factor:

Equation 10 Bending correction factor equals the true load transfer efficiency divided by the measured load transfer efficiency, which equals 1 plus the sum total of the model parameter 1 divided by the radius of the relative stiffness, plus the sum total of the model parameter 2 divided by the radius of relative stiffness squared.                                                                     (10)

where:

-          B = the bending correction factor.

-          LTEtrue = the true LTE.

-          LTEmes = the measured LTE.

-          Radius of relative stiffness = the radius of relative stiffness.

-          a1 and a2 = model parameters.

A regression analysis was performed to determine coefficients a1 and a2 and the following expression for the bending factor was obtained:

Equation 11 Bending correction factor equals 1 minus the sum total of 6.285 divided by the radius of relative stiffness in millimeters, plus the sum total of 36661 divided by the radius of relative stiffness in millimeters squared.                                                                       (11)

R2  = 0.9998

SEE = 0.000495

where:

-          Radius of relative stiffness = a radius of relative stiffness in mm.

Although the correction factor has excellent statistics, its practical applicability is quite limited because the radius of relative stiffness may be unknown.  To avoid this limitation, a correction factor based on sensor deflections located on the leave slab was used to correct measured LTE.  The AREA parameter was used for this purpose.  This parameter combines the effect of several measured deflections in the basin and is defined as follows:

Equation 12 AREA parameter equals the sum total of 1 divided by 2 times the measure deflections at initial sensor position, times the total quantity of measured deflection at initial sensor position times the distance in millimeters between the first and second sensor plus the sum total of the summation from the second sensor position to the second from the last sensor position of the measured deflections at all sensor positions times multiplied by the difference in distance in millimeters of the corresponding adjacent sensor positions, plus the measured deflection at the sensor position times multiplied by the difference of the distance of the last sensor position and the next-to-the-last sensor position.           (12)

where:

-          Wi = measured deflections (i = 0, n).

-          n    = number of FWD sensors minus 1.

-          ri    = distances between the center of the load plate and sensors in mm.

The AREA parameter has been used extensively to analyze concrete pavement deflection basins since 1980.  Ioannides et al. (1989) identified the unique relationship between AREA and the radius of relative stiffness.  The AREA parameter is not truly an area, but rather has dimensions of length, since it is normalized with respect to one of the measured deflections in order to remove the effects of load magnitude.  For a given number and configuration of deflection sensors, AREA may be computed using the trapezoidal rule.  It was found that the correction factor depends on the AREA parameter and the magnitude of the measured LTE itself.  The following relationship between the true and measured LTE was proposed:

Equation 13 Predicted loaded transfer efficiency equals the sum total of the regression coefficient subscript 1 times the measured loaded transfer efficiency, plus the sum total of the regression coefficient subscript 2 times the sum total of the measured loaded transfer efficiency divided by AREA parameter, plus the sum total of the regression coefficient subscript 3 times the sum total of the measured loaded transfer efficiency divided by the AREA parameter, squared,  plus the sum total of the regression coefficient subscript 4 times the sum total of the measured loaded transfer efficiency divided by the AREA parameter, cubed.    (13)

where:

-          K1, K2, K3, and K4 = regression coefficients, depending on the sensor configuration used in AREA calculation.

Using the results from the 375 ISLAB2000 runs, the regression analysis for determining K1, K2, K3, and K4 was performed for the following sensor configurations:

Table 1 presents the determined coefficients along with basic statistics.

Table 1.  Regression coefficients for bending correction factors.

Test

Sensor Configuration

K1

K2

K3

K4

R2

SEE

Approach

O5

0.929710

  23.61

   -45.73

 1171.28

0.99729

0.108298

Approach

O4

0.790073

116.74

 -166.34

  856.58

0.99724

0.231601

Leave

C6

0.924255

  92.96

 -508.52

5302.13

0.99727

0.163864

Leave

C5

0.806827

155.35

 -491.51

3472.81

0.99715

0.366376

Leave

O5

0.923371

  72.03

 -182.53

2143.95

0.99729

0.114590

Leave

O4

0.764173

129.91

-209.29

1105.75

0.99720

0.296859

To investigate the applicability of the bending factor developed in the study, bending factors were calculated for more than 600,000 FWD basins from the LTPP database. Figures 5 and 6 present the distribution of these factors for the FWD basins from approach slab test calculated using the O5 sensor configuration and for the basins from the leave slab test calculated using the C6 sensor configuration.  In most cases, the correction was less than 4 percent.

During the bending factor testing using data from the LTPP database, it was found that significant discrepancies existed in measured LTEs for the same joints measured with a load plate placed on the approach and the leave side of the joint.  The correction factors presented above did not significantly reduce these discrepancies.  An additional investigation was conducted in an attempt to resolve this issue.

The finite element model used in the development of the correction factors did not account for differences in pavement responses as a result of the location of the loaded plate.  Also, the PCC slabs were assumed to be in full contact with the subgrade.  In reality, the application of a large number of a heavy axle loads moving in one direction may cause the formation of permanent voids under the leave side of the joint.  An assumption was made that those voids are responsible for the discrepancy between leave and approach test results.  In addition, PCC slab curling causes slabs to separate, creating temporary voids.  This effect also was analyzed in this study.

Figure 5. Distribution of bending correction factors for approach slab testing based on O5 sensor configuration.  The bending correction factor, from 0.98 to 1.08, is graphed on the horizontal axis. The cumulative percentage of measurements is graphed on the vertical axis. The graph increases in an S-shaped curve beginning at the lowest percentage at 0.98 bending correction factor. The line starts to increase in a steep incline at 2 percent and 0.99 bending factor to 92 percent and 1.03 bending factor. The line then tapers off to a horizontal line.

Figure 5.  Distribution of bending correction factors for approach slab testing based on O5 sensor configuration.

A factorial of 504 finite element runs with different temperature gradients through slab thickness, joint stiffness, and subgrade stiffness was performed.  In all cases, a 1.52-m (5-ft)-long permanent void occupying the entire leave slab width was modeled.  In half of the cases, the load was placed at the approach side of the joint, and in the other half it was placed at the leave side.  The load transfer values from the deflections induced by the loads placed at the approach side of the joint, LTE1, were compared with the corresponding load transfer values calculated for the same systems but the load placed on the leave side on the joint, LTE2.  Figure 7 presents a comparison of those LTEs.  As shown in the figure, the presence of a void significantly skewed computed LTE values.  The figure shows that LTEs computed from the leave slab test should be lower than those computed from the approach slab test.  Moreover, this effect may be much more pronounced than the effect of slab bending. 

Figure 6. Distribution of bending correction factors for leave slab testing based on C6 sensor configuration. The figure is a line graph. Bending correction factor, of 0.98 to 1.08, is graphed on the horizontal axis. Cumulative percentage of measurements is graphed on the vertical axis. The graph increases in an S-shaped curve beginning at the origin (0, 0.98) and rising in a steep angle. The line tapers off at (86, 1.03) and continues to increase but at a sloping horizontal angle, never reaching 100.

Figure 6.  Distribution of bending correction factors for leave slab testing based on C6 sensor configuration.

Figure 7. Comparison of loaded transfer efficiencies predicted using ISLAB2000 for the approach slab test (LTE subscript 1) and the leave slab test (LTE subscript 2) for pavements with voids. The ISLAB 2000 for the approach slab, percent, is graphed on the horizontal axis. The leave slab test for pavements with voids, percent, is graphed on the vertical axis. The figure has a linear slope beginning at the origin and increasing at a 45-degree angle. There are scatter plots below the slop also increasing. As the approach slab increases, so does the leave slab. The figure shows that the leave slab test should be lower than the approach slab.

Figure 7.  Comparison of LTEs predicted using ISLAB2000 for the approach slab test (LTE1) and the leave slab test (LTE2) for pavements with voids.

This analysis agrees with the analysis of LTPP data that discrepancies may exist between the LTEs computed from leave and approach tests; however, in the vast majority of cases analyzed, an opposite trend was observed.  As will be discussed in chapter 3, LTEs computed from the leave slab test usually were higher than LTEs computed from the approach slab test.  At this stage, no mechanistic explanation of this phenomenon was found.  More research is needed to address this issue.  Meanwhile, considering that the correction factors presented above did not change the calculated LTE values significantly, and that void/curling correction should be incorporated into this factor, a decision was made not to use the bending factor in this study.

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