A Study of Static & Dynamic Modulus of Elasticity of Concrete PDF
A Study of Static & Dynamic Modulus of Elasticity of Concrete PDF
October 2008
A Study of Static and Dynamic Modulus of Elasticity of Concrete
Principal Investigator : John S. Popovics, Ph.D., P.E., FACI. Civil & Environmental Engineering,
University of Illinois, Urbana, IL
Contact: 217 244 0843 (voice); 217 265 8040 (fax); [email protected]
However, once a structure is erected the in situ elastic properties cannot be measured directly
without damaging the structure itself. Most often E is inferred from the compressive strength (fc)
of companion cylinders, rather than being measured directly, through the application of
established empirical relations. This approach often leads to overly conservative results because,
in order to meet the minimum E requirement, concrete with much higher fc is used than the
specification requires, which leads to unnecessarily high material costs. Enhanced understanding
of the relation between E and compressive strength, with respect to different types of concrete,
would improve the efficacy of the estimation of E from strength.
Non-destructive dynamic methods can be used to estimate in-place E, but the meaning of the
obtained dynamic modulus is uncertain because Ed is known to be different (higher) from that
obtained by direct static testing of a cylinder drawn from the structure. Concrete is expected to
show a nonlinear dependence between stress and strain, even at low values of deformation
caused by quasi-static tests and dynamic tests based on stress-wave propagation [Powers 1938;
Bell 1984]. Quasi-static experiments show a nonlinear dependence between stress and strain
even at infinitesimal values of deformation for a wide range of materials (metal, stone, concrete,
wood, glass, polymers, etc.). The nonlinear stress-strain relationship for concrete is well
described by a quadratic parabola [Shkolnik 1996]. The non-linear behavior provides the bases
for the conventionally accepted view why Ed is higher than E, since concrete is subjected to very
small strains in dynamic testing. [Neville 1997]. Although this argument is conceptually
satisfying, experimental test results show that this may not in fact be the the case: low-strain
static test data agree with higher strain (stressed to 40% of ultimate) static test data, as shown in
Figure 1. Therefore it is important to understand precisely how and why the E and Ed are related
to each other.
Furthermore, it is known that Ed values for a given concrete obtained by different dynamic tests
do not agree with each other [Philleo 1955]. In general, Ed obtained from pulsed wave
propagation measurements are significantly higher than those obtained from vibration resonance
measurements carried out on the same specimen. Ed varies significantly even within one type of
measurement: Ed computed from vibration resonance of prismatic beams is known to be, on
average, significantly higher than that computed from cylinders for the same concrete mixture.
The work will have significant impact on the technical community since more accurate
monitoring and in-place control of concrete E would be enabled. As such the work would likely
influence existing testing standards, concrete design specifications, concrete quality control,
inspection of structures, and development of commercial NDE equipment. The work includes
significant in-kind contributions from a certified construction materials testing firm and a
concrete ready- mix supplier.
Scope of experimental tests
Test samples
The investigating team members received from CTLGroup a test data matrix from over 200 high
strength concrete samples (150mmx300mm and 100mmx200mm cylinders).
Measured
compressive strength of these samples ranged from 24MPa to 161MPa, E from 25.3GPa to
61.7GPa, and testing age from 4 days to 730 days.
E=
V P (1+ )(1 2 )
(1 )
f
f
d = A1 2 + B1 2 + C1
f1
f1
2 f1Ro
Ed = 2(1 + d )
1
fn
where Ed = dynamic modulus of elasticity, Pa, ?d = dynamic Poissons Ratio, ? = density, kg/m3 ,
1
2
R0 = radius of the cylinder, m, f n = A2 ( d ) + B2 ( d ) + C2 and A2 , B2 , and C2 are constants based
on the dimensions of the cylinder [Subramaniam et al., 2000].
For the Loves correction method, the fundamental longitudinal resonance frequency (f) of a
cylinder with length L and diameter d is measured, and a value of Poissons ratio is assumed.
Then Ed is computed as
2 2 ( 2m 1) 2 2
)
L2
( 2m 1) 2
4 f 2 L2 (1+
E=
where m=1 for the fundamental mode and = d/2v2 for cylinders. In this work, only the
fundamental mode in cylinders was measured.
Project Findings
The project findings are based on the analysis of the experimental test data.
Dynamic test methods on control sample
The dynamic vibration methods were first applied to a homogenous, uniform material with
known properties. The tests were carried out on beam and cylinder samples made from
aluminum. The geometries of the aluminum samples replicate those of the concrete samples. The
material properties of the aluminum were verified using ultrasonic P-wave and S-wave velocity
measurements: E = 10.4x106 psi GPa and Poissons ratio = 0.33. All four vibration tests were
applied to the specimens: ASTM C215 using transverse vibrations (both beam and cylinder),
ASTM C215 using longitudinal vibrations (both beam and cylinder), Loves correction for
longitudinal vibrations (cylinders only) and the two-frequency method (cylinders only). The
results from the tests are shown in Figure 3. The ASTM C215 measurements with longitudinal
vibration consistently underestimate the actual value, while the ASTM C215 transverse
measurements overestimate the actual value. The amount of under- and over-estimation increases
when cylinders are used. The largest discrepancy with the actual value is provided by ASTM
C215 using longitudinal resonances with cylindrical specimens. The Loves correction and two
frequency methods give accurate estimates of E for Aluminum cylinders, when correct values of
Poissons ratio are assumed in the computations.
Dynamic test methods on concrete samples
The test results from the aluminum samples indicate that the various vibration tests provide
different values of E for the same sample. This was confirmed by tests carried out on high
strength concrete cylinder samples provided by CTLGroup. Figure 4 shows a comparison of Ed
values obtained by ASTM C215 using the transverse and longitudinal modes measured on the
same sample. The results confirm that the transverse vibration measurements are consistently
higher than the longitudinal measurements. This discrepancy appears to widen as the Ed value
increases. The agreement between longitudinal and transverse resonance values improves for all
values of Ed when the Loves correction in applied to the longitudinal resonances, where a
Poissons ratio of 0.25 was assumed. This is shown in Figure 5.
One issue related to the use of Loves correction for cylinders, ultrasonic pulse velocity and
ASTM C215 using the transverse mode methods is that a value of Poissons ratio must be
assumed or measured for the computation. Furthermore, it is known that both static and dynamic
Poissons ratio vary across a broad range of values for concrete, usually between 0.15 and 0.25,
and that the dynamic and static values of Poissons ratio are not strongly related to each other.
Considering this, incorrect estimates of Poissons ratio for concrete are likely when a single
value is assumed for a given concrete. The effects of this possible mis-estimation on each
method should be considered. Examination of the formulas for each method reveals that misestimates of Poissons ratio affect differently the accuracy of E estimates from each method.
ASTM C215 using longitudinal resonances and the two-frequency method for cylinders do not
rely on an estimate of Poissons ratio in the computation. Loves correction for cylinders shows
only a small disruption of E estimates owing to mis-estimation of Poissons ratio. ASTM C215
using transverse modes shows moderate disruption, while ultrasonic P-wave velocity shows large
disruption. This last point suggests that E estimates from P-wave velocity are not reliable for
concrete; this is confirmed by test results, shown in Figure 6, that show that Ed measurements
from P-wave velocity significantly over-predict those obtained from ASTM C215 longitudinal
vibration, when median values of Poissons ratio are assumed for the concrete. Interestingly, the
agreement between E measured by the two methods improves when an artificially high value of
Poissons ratio (0.28) is assumed for all concretes, as shown in Figure 7.
The behavior shown in Figures 6 and 7 is unexpected, and cannot be readily explained by
conventional theory for elastic and homogeneous materials. It is possible, however, that
composite nature of concrete has a role in this behavior. To investigate this, a set of special twophase (cement paste and coarse aggregate) concrete samples with varying proportions of the
phases were produced. Photos of the range two-phase concrete samples are shown in Figure 8.
These samples represent the extreme cases of the two-phase mixture arrangement. Any global
property of the two-phase composite should be bounded by the Voigt and Ruess (parallel and
series configurations, respectively) combinations of the properties of the individual phases
themselves. Several composite material properties were investigated across the different phase
compositions. As expected, the material mass density of the composite follows the upper bound
of the Voigt model combination of the density of the individual phases. Interestingly, the global
static and dynamic (ASTM C215 using longitudinal mode) moduli are different from each other
across the different phase compositions, as shown in Figure 9. Dynamic modulus is consistently
higher than static modulus, and it tends to follow the upper-bound Voigt mixture rule. On the
other hand, the static modulus lies in-between the upper and lower bounds. In other words, static
and dynamic moduli follow different mixture behaviors in composite materia ls such as concrete.
It is this difference that may cause dynamic moduli to be higher than static moduli in conc rete. It
should be noted that this difference between static and dynamic moduli is not observed in cement
paste samples, which are more homogeneous.
Comparison of E and Ed for high strength concrete
Several attempts have been made to correlate static (E) and dynamic (Ed) moduli for concrete.
The simplest of these empirical relations is proposed by Lydon and Balendran [Neville, 1997]
E = 0.83 Ed .
Another empirical relationship for concretes elastic moduli was proposed by Swamy and
Bandyopadhyay and is now accepted as part of British testing standard BS8110 Part2:
E = 1.25Ed 19
where both units of E and Ed are in GPa. This expression does not apply for lightweight
concretes or concrete that contains more than 500 kg of ceme nt per cubic meter of concrete
[Neville, 1997].
Shkolnik proposed relations between dynamic and static moduli for concrete, based on thermofluctation theory, which depends on strain rate of loading and the temperature [Shkolnik 1996;
Shkolnik 2005]. The vibration of atoms with average kinetic energy causes stresses on the
atomic bonds of the same magnitude as the strength of materials. Thus at atomic- molecular level
the fracture of a material is controlled by breaking of atomic bonds at thermal fluctuation. The
thermal fluctuations of energy, exceeding the average thermal energy of atoms, naturally, depend
on the time or rate of loading and predetermine the mechanical properties (modulus of elasticity,
strength) of real materials. By coupling the nonlinear stress-strain material model with the kinetic
(thermofluctuation) theory of solid strength, a physically proved expression was applied for
modulus of elasticity of concrete of a material. Considering the samples provided by the
CTLGroup, which were tested at a nominal value of stress/strength ratio of 0.24, rate of loading
0.24MPa/s, and temperature of 20 o C 9, the thermoflucutaion theory gives:
E = Ed 5864
for units of MPa. The constant at the right side of equation is determined in accordance with
ASTM C469.
For both lightweight and normal concretes, Popovics suggested a more general relationship
between the static and dynamic moduli as a function of the density of the concrete
Es = kEd1.4 1
where k = 0.23 for units of psi and ? is density, lbs/ft3 [Popovics, 1975]. Whatever the relation
between the moduli, it is thought to be unaffected by air entrainment, method of curing,
condition at test, or the type of cement used [Neville, 1997].
In Figures 10-12, the fit of the four shown equations that relate Ed to E are shown for the highstrength concrete cylinder data obtained from CTLGroup. In addition, a best fit fine to the data is
also shown. Only the fundamental longitudinal and fundamental transverse modes were provided
by CTLGroup, so only three analyses methods are possible: ASTM C215 using longitudinal
resonance for cylinders, ASTM C215 using transverse resonance for cylinders, and Loves
correction for longitudinal vibration. For ASTM C215 using transverse vibration and the Loves
correction method, a Poissons ratio of 0.25 was assumed. The predictions for Ed are shown in
Figures 10, 11 and 12, respectively. The mean absolute errors (in GPa) for each fit are also in the
figures. In general the data show a linear, or nearly linear, dependence between E and Ed across
the full range of modulus values. In all three data sets, the fitted line shows the best fit (lowest
mean absolute error) and the equation proposed by Lyndon and Balendran show the poorest fit
(highest mean absolute errors); it appears that the latter equation is inappropriate for high
strength concrete. The remaining three equations show reasonable fit, where the mean absolute
error is always below 2GPa. The best fit line exhibits absolute mean error below 1.3 GPa. It is
interesting to note that the predictions using ASTM C215 with transverse resonance show
approximately 10% higher error for the fit line, suggesting that the transverse resonance data
may exhibit a higher degree of scatter.
More complete statistical analysis of the predicted E data is shown in Tables 1 and 2. In Table 1,
column 1 has the lowest Mean E, STDEV, CV, and P-values at longitudinal and transverse
vibrations. The relative mean deviations are all negative for longitudinal vibration, but for
transverse vibration two of them are negative (columns 1 and 3), and two- positive (columns 2
and 4). The MAPE does not exceed 10% at longitudinal and transverse vibration. The column 2
illustrates the highest P-value at longitudinal vibration (P-value = 0.282), and column 4 at
transverse vibration (P-value=0.827). The relative difference between the Mean E for column 4
at transverse vibration (E=49.157 GPa), and the average value of the Mean E for columns 2 and
3, i.e. (49.777+48.774)/2= 49.276 GPa, equals 0.2%. In Table 2, the relative mean deviation
from E (C 469) decreased, and simultaneously column 2 and 4 changed the sign from negative to
positive. The values of MAPE decreased (except column 2), and The P-values increased
(columns 2, 3, 4) in comparison with corresponding data for longitudinal vibrations in Table 1.
As a consequence of accepted assumptions ( =0.25), the mean values of predicted dynamic
modulus of elasticity obtained for longitudinal vibrations (Table 2) and transverse vibrations
(Table 1) are equal within the limits of errors of measurements and calculations.
Conclusions
The following conclusions are drawn, based on the results presented and analyzed in this report:
o The dynamic and static modulus tests carried out at the University of Illinois laboratories
confirm that Ed values vary, depending on method of measurement: ultrasonic pulse
velocity (UPV ASTM C597) gives the highest predictions of Ed, and longitudinal
vibration the lowest. Of the vibration methods, longitudinal vibration of cylinders
(L/D=2) give the least accurate prediction of Ed. A modified method (Loves correction)
was introduced for longitudinal resonances for cylinders, which gives more accurate Ed
results that agree with those from the transverse method. The Loves correction was
applied to the CTLGroup data, providing a more accurate estimation of E from Ed, with
an absolute mean error of estimate of 1.10 GPa. UPV may be able to provide accurate
estimates of E, but only if excessively high values of Poissons ratio are assumed in the
calculations. More study in this area is warranted. It was found that the longitudinal and
transverse methods of ASTM C215 give different results from the same sample:
transverse provides higher dynamic modulus (E d) estimates.
o As expected, Ed is always greater than E for concrete; however, this behavior may be
caused by the composite nature of concrete, rather than non-linear behavior of concrete
exposed by varying strain levels. Stress strain curvature analyses show that the expected
non- linear behavior at higher strain levels associated with the static test is not clearly
seen. Tests on specially prepared two-phase composite samples, however, do indicate the
influenc e of composite nature on the E vs. Ed behavior.
o Several existing relations between dynamic and static modulus (E) were evaluated for the
data. The results show that the relation between Ed and E is affected by nature of Ed data:
the best relation to apply depends on how Ed was obtained. In general, though, the
relation between Ed and E appears to be linear, or nearly linear. For longitudinal
vibration, the British Standard equation (BS8110, Part 2) gives the best prediction of E,
with an average squared error of 1.36 GPa. For transverse vibration and Loves
correction of longitudinal vibration, the equation proposed by S. Popovics gives best
predictions for E: 1.41 GPa and 1.10 GPa, respectively.
o A thorough understanding of the various factors affecting Ed and E in concrete should
enable improved prediction of Ed. Once this is established, vibration measurements
should offer an effective, non-destructive, inexpensive and rugged method to estimate E
for concrete.
Acknowledgements
A portion of this work was supported through funds from the American Concrete Institute
Concrete Research Council program. The authors appreciate and thank Mike Pistilli of Prairie
Materials for providing strength and modulus data generated over a period of several years.
Technical contributions from Joni Ranchero of CTLGroup and additional assistance from Sara
Alzate, Taekeun Oh and Kristjan Kranjc are acknowledged and appreciated. The authors are
grateful to ACI Committee 363 High Strength concrete for their strong support of this work.
Cited References
Bell, J.F. (1984) Mechanics of solids: The experimental foundations of solid mechanics,
Springler-Verlag, Berlin.
Godfrey Jr., K.A. (1987) Concrete strength record jumps 36%, Civil Engineering. Volume
57(10), pp. 84-88.
Leet, K.M. and Bernal, D. (1997) Reinforced Concrete Design. 3rd Ed. McGraw-Hill
Companies, Inc. New York.
Love, A.E.H. (1944) The Mathematical Theory of Elasticity, Dover Publications Inc., New
York.
Neville, A.M. (1997) Properties of Concrete 4th Ed. John Wiley & Sons Inc., New York.
Philleo, R.E. (1955) Comparison of Results of Three Methods for Determining Young's
Modulus of Elasticity of Concrete, Journal of the American Concrete Institute, Volume
26(5), pp. 461-469.
Popovics, S. (1975). Verification of relationships between mechanical properties of concretelike materials. Materials and Structures. Vol 8. No. 45. pp. 183-191.
Powers, T.C. (1938) Measuring Youngs modulus of elasticity by means of sonic vibrations,
Proceedings of the Forty-first Annual Meeting of the American Society for Testing
Materials, Volume 38(2), pp. 460-470.
Shkolnik I.E. (2005) Effect of nonlinear response of concrete on its elastic modulus and
strength, Cement and Concrete Composites, Volume 27, pp. 747-757.
Shkolnik, I E. (1996) Evaluation of dynamic strength of concrete from results of static tests,
Journal of Engineering Mechanics, Volume 122(12), pp. 135-138.
Subramaniam, K.V., Popovics, J.S. and Shah, S.P. Determining elastic properties of concrete
using vibrational resonance frequencies of standard test cylinders, Cement, Concrete &
Aggregates (ASTM), Volume 22, Number 2, pp. 81-89, 2000.
Figures
50
Concrete
Specimens
45
40
35
Paste
Specimens
30
25
Limestone
20
River Gravel
Air-entrained
15
High Strength
unity
10
0
0
10
15
20
25
30
35
40
45
50
Figur e 1. Comparison of low-strain (initial tangent modulus) and high strain (conventional chord modulus) for static
compressive tests carried out on concrete and paste samples.
Figur e 2. Testing configuration for standard vibrational resonance tests (Naik et a. 2004)
Figur e 3. Comparison of longitudinal and transverse Ed values obtained using ASTM C215 method on aluminum
samples. The expected E value is 10.4x106 psi and is 0.33.
Figur e 4. Comparison of longitudinal and transverse Ed values obtained using ASTM C215 method on concrete
cylinders. Test data provided by CTLGroup.
Lo n g i t u d i n a l E- - Lo v e ' s Co r r e ct i o n ( p si )
11000000
10000000
9000000
8000000
7000000
6000000
5000000
4000000
4.000E+ 06
5.000E+ 06
6.000E+ 06
7.000E+ 06
8.000E+ 06
9.000E+ 06
1.000E+ 07
1.100E+ 07
Tr a n sv e r se E ( p si )
Figur e 5. Comparison of longitudinal and transverse Ed values. Longitudinal values obtained using Loves
correction (assumed v = 0.25) and transverse values obtained using ASTM C215 method on concrete cylinders. Test
data provided by CTLGroup.
60
concrete specimens
? d = 0.20
50
40
Paste Specimens
? d = 0.25
30
Limestone
River Gravel
20
Air-Entrained
High Strength
10
0
0
10
20
30
40
50
60
Figur e 6. Comparison of Ed values obtained by longitudinal resonance frequency (ASTM C215) and ultrasonic
pulse velocity (ASTM C597) tests carried out on concrete and paste cylinders. Median values of Possions ratio
were assumed for the predictions based on pulse velocity.
60
concrete specimens
? u = 0.28
50
40
Paste Specimens
? u = 0.29
30
Limestone
20
River Gravel
Air-Entrained
High Strength
10
0
0
10
20
30
40
50
60
Figur e 7. Comparison of Ed values obtained by longitudinal resonance frequency (ASTM C215) and ultrasonic
pulse velocity (ASTM C597) tests carried out on concrete and paste cylinders. Unusually high values of Possions
ratio were assumed for the predictions based on pulse velocity.
Figur e 8. Two-phase limestone and cement paste samples, with limestone content ranging from 0 to 75% by mass.
50
Line of best fit
y = 0.2272x + 19.453
R2 = 0.9802
45
40
40
Static Modulus (Chord), GPa
45
35
30
35
Upper bound (Voight)
30
25
20
Lower bound (Ruess)
20
15
10
15
0
10
20
30
40
50
60
70
80
90
100
10
20
30
% aggregate by mass
40
50
60
70
80
90
% aggregate by mass
Figur e 9. Variation of dynamic (left) and static (right) Youngs modulus as a function of aggregate content in the
two-phase composite concrete specimens. Ed determined using longitudinal resonance following the ASTM C215
procedure.
70.0
60.0
E Static(GPa)
y = 1.197x - 15.43
R = 0.967
50.0
Actual Values
40.0
30.0
20.0
30
40
50
60
70
E Dynamic(GPa)
Figur e 10. Comparison of E (ASTM C469) vs. Ed (ASTM C215) data (points) using longitudinal resonance for
concrete cylinders. Various relations between E and Ed are also shown as lines. The mean absolute error (in GPA)
for each type of fit curve is shown in the table underneath. Data obtained from CTLGroup.
100
70.0
60.0
E Static(GPa)
y = 1.137x - 13.27
R = 0.961
50.0
Actual Values
Lydon & Balendran
British Standard
40.0
Popovics
Thermofluctuation
Linear (Actual Values)
30.0
20.0
30
40
50
60
70
E Dynamic(GPa)
Lyd on &
Balen d ran
Pop ovics
1.524
1.410
4.044
Fitted line
ion
1.500
1.280
Figur e 11. Comparison of E (ASTM C469) vs. Ed (ASTM C215 using transverse resonance) data (points) for
concrete cylinders. Various relations between E and Ed are also shown as lines. The mean absolute error (in GPA)
for each type of fit curve is shown in the table underneath. Data obtained from CTLGroup.
70.0
60.0
E Static(GPa)
y = 1.173x - 15.38
R = 0.967
50.0
Actual Values
Lydon & Balendran
40.0
British Standard
Popovics
Thermofluctuation
30.0
20.0
30
40
50
60
70
E Dynamic(GPa)
Lyd on & Balen d ran
Pop ovics
Fitted lin e
4.020
1.362
1.192
1.575
1.175
Figur e 12. Comparison of E (ASTM C469) vs. Ed (Loves correction using longitudinal resonance) data (points) for
concrete cylinders. Various relations between E and Ed are also shown as lines. The mean absolute error (in GPA)
for each type of fit curve is shown in the table underneath. Data obtained from CTLGroup.
Tables
Table 1. Statistical comparison of E computed by the various relations and Ed (ASTM C215). The mean value of
E obtained by ASTM C469 is 49.313 GPa, STDEV=8.39GPa and CV=0.17.
Longitudinal vibrations
E= 0.83Ed, (MPa).
E, GPa
Mean E, GPa
STDEV, GPa
CV.
Relative mean
deviation from
E (C469)
MAPE
P-value
44.815
5.653
0.130
48.492
8.513
0.180
47.484
8.070
0.170
48.129
6.811
0.140
-0.082
-0.017
-0.036
-0.017
0.094
0.000
0.031
0.282
0.042
0.014
0.041
0.086
Tr ansver se
Mean E, GPa
STDEV, GPa
CV.
Relative mean
deviation from
E (C469)
MAPE
P-value
E=1.25 Ed - 19,
E, GPa
Data
E= Ed-5864, (MPa).
E, GPa
vibrations
1
45.688
6.006
0.130
2
49.777
9.045
0.180
3
48.774
8.648
0.180
4
49.157
7.236
0.150
-0.065
0.008
-0.011
0.003
0.079
0.000
0.033
0.557
0.030
0.485
0.033
0.827
Table 2. Statistical comparison of E computed by the various relations and Ed using Loves correction (v=0.25). The
mean value of E obtained by ASTM C469 is 49.313 GPa, STDEV=8.39GPa and CV=0.17.
Longitudinal vibrations
E=1.25 Ed - 19,
E, GPa
Data
E= 0.83Ed, (MPa).
E, GPa
E= Ed-5864, (MPa).
E, GPa
Mean E, GPa
STDEV, GPa
CV.
Relative mean
deviation from
E (C469)
MAPE
*P value
45.711
5.766
0.126
49.842
8.683
0.174
48.819
8.297
0.170
49.209
6.947
0.141
-0.063
0.010
-0.009
0.005
0.080
0.000
0.031
0.492
0.028
0.512
0.036
0.882