Novak 1976
Novak 1976
Novak 1976
Frtculty of Engineering Science, University of Western Ontario, London, Ont., Crtnndrc N6A 3K7
Received November 12, 1975
Accepted June 14, 1976
Dynamic experiments with small pile foundations were conducted in the field. Individual piles
and a group of piles were tested. The experimental data were compared with theoretical predic-
tions made on the basis of a theory described previously by Novak. It was found that the theory
predicts the general characteristics of the response very well. Resonant amplitudes and natural
frequencies can be predicted reasonably well when a shear modulus of the soil derived from a
static test of a pile is used and when a correction for the effect of pile group (interaction) is made.
Des essais dynamiques sur des fondations de pieux modeles ont CtC rCalisCs sur le terrain. Des
pieux isolCs et un groupe de pieux ont Cte soumis B essais. Les resultats experimentaux ont it15
compares aux previsions theoriques faites a partir d'une thiorie proposCe par Novak. On a
constate que la thCorie permet une trbs bonne prevision des caractkristiques gCnerales de la
rkponse. Les amplitudes de rksonnance et les frequences naturelles peuvent ktre prevues
raisonnablement bien lorsque le module de cisaillement du sol diduit de i'essai statique d'un pieu
est utilisi et lorsqu'on applique une correction pour tenir compte des effets de groupe.
[ ~ r a d u i t - p alar revue]
Can. Geotech. J., 13,372 (1976)
For personal use only.
be adjusted for large amplitudes by means of depends on the slenderness ratio l/ro, where
equivalent linearization. 1 = length of the pile and yo = radius of the
pile, and the stiffness of the soil. The condi-
Experiments tions of the site suggested a maximum pile
Can. Geotech. J. Downloaded from www.nrcresearchpress.com by UNIVERSITY OF ARIZONA LIBRARY on 09/06/13
The Site and Soil Properties diameter of about 4.5 in. (1 1.5 cm). With
A suitable site was found on the campus of respect to the ease of installation standard
the University of Western Ontario, close to black steel pipes were chosen with diameters
the Thames River. It is a layer of fine silty of 3.5 and 2.4 in. (9 and 6 cm). The following
sand approximately 6 to 7 ft (1.8 to 2.1 m) test situations were examined: ( a ) a single
thick underlain by gravel and till. 3.5 in. (9 cm) diameter pile, ( b ) a single 2.4
The soil properties were determined by in in. ( 6 cm) diameter pile, (c) a group of four
situ and laboratory tests. The in sifir density 2.4 in. (6 cm) diameter piles, and ( d ) a group
was determined by means of the Rubber Bal- of four 2.4 in. (6 cm) diameter piles with soil
loon Method in an excavation. Verifications packed beneath the pile cap. The properties of
of the soil profile and soil samples were ob- the piles used are given in Table 2.
tained through the use of a hand augered test The piles were hydraulically driven into the
hole. soil until they reached the solid layer and by
The shear wave velocity of the soil was dropping a 140 Ib (62 kg) hammer 30 in. (76
established by means of the steady state vibra- cm). The piles used in the group were cut
tion technique. The ground motion was excited short and did not reach the solid layer.
by the piles and the results were clearer and T o examine the effect of the pile length,
more consistent than is usual with soil alone. several stiffness and damping parameters were
For personal use only.
The depth of the hard layer was evident from calculated from formulas given in Novak
the variation of the shear wave velocity with (1974a) as a function of slenderness. These
frequency and was also verified by a seismic parameters are shown in Fig. 1. It can be seen
test. The soil properties useful in dynamic that for the slenderness used, the stiffness and
analysis are givcn in Table 1. More details on
the shear wave velocities measured are given in
the Appendix.
Test Piles
The size of the piles used was limited,
among other considerations, by the force and
frequency range of the oscillator (exciter)
available and by the requirement that thc mo-
tion of the pile tips be kept as small as pos-
sible. The satisfaction of the latter condition is
facilitated by the theoretical finding that it is
only the upper part of the pile that actively
participates in the motion and determines the
pile stiffness and damping. This active length
,I 1, 1 ,
TABLE
1. Properties of soil (fine silty sand)
TABLE
2. Properties of piles
TABLE
3. Properties of test footings
Height of Height of
Mass moment
Can. Geotech. J. Downloaded from www.nrcresearchpress.com by UNIVERSITY OF ARIZONA LIBRARY on 09/06/13
(6 cm) diameter pile were established in two FIG.3. Static displacement of pile head vs. horizon-
directions differing by 90". tal load.
From these static tests an apparent shear
modulus of the soil can be established by com- static test is significantly lower than that estab-
paring the measured displacements with those lished from the wave propagation test. In other
predicted theoretically. Such a comparison is words, the theory significantly under-estimates
facilitated by fair linearity of the relationships the static response if the shear modulus estab-
shown in Fig. 3 for displacements smaller than lished from wave propagation is used. This
about 0.1 in. (2.5 mm). difference seems to be due to two main factors
Using Fig. 3 of the theoretical static solution which are the separation of the pile from the
in Poulos (1972) and assuming Poisson's ratio soil in the region of tension and the diminishing
v = 0.25, an iterative procedure yields the of the shear modulus with decreasing depth.
apparent shear modulus and wave velocity Both of these factors, indicated in Fig. 4,
given in Table 1 and denoted Vst. depend on confining pressure which is very
The shear wave velocity obtained from the small in the upper region of the soil layer. It
CAN. GEOTECH J. VOL. 13, 1976
vibrations are
smaller than that obtained from wave propaga- If the vibrating system were perfectly linear
tion which reaches the deeper regions of soil. all the dimensionless curves shown in each
The difference between V , and V,, would figure should be identical. The amount by
probably be much smaller if V,,. were deter- which they differ, indicates the degree of non-
mined for the surface layer of the soil or if the linearity. The nonlinearity is demonstrated by
confining pressure were greater as it is in a real shifting of the resonant frequency and variation
pile foundation. in the resonant amplitudes. The nonlinearity,
The effect of pile separation indicated in though distinct, is not very strong. It could be
Fig. 4a, depends on the magnitude of the analyzed using a technique described by Novak
applied load and manifests itself by the non- (1971 ). However, for the purpose of this
linearity of the response. However, no strong study it is sufficient to use average natural
nonlinearity was observed in the experiments frequencies and resonant amplitudes. These
and this seems to indicate that the reduction are listed in Table 4.
of the shear modulus towards the surface was In this table, subscript 0 refers to the natural
a more important factor than the separation of frequency for vertical vibration, subscripts 1
the pile. In any case, it may be assumed that and 2 refer to the first or second resonance
the shear modulus obtained from the static test with horizontal excitation. The magnitudes not
incorporates in an approximate way the re- given either lay beyond the frequency range
duced stiffness of the surface layer of the soil of the equipment or were not obtained reliably
and the lack of bond between the soil and the because of malfunctioning of the oscillator.
pile. The group of four piles was tested for the
horizontal excitation in two directions. The
Dynamic Tests force acted in the direction of either the shorter
A few typical sets of response curves mea-
or the larger distance x , between the piles. An
sured with vertical or horizontal excitation are
example of the horizontal response of the pile
plotted in Figs. 5 to 8. Each figure shows a
group is shown in Fig. 7. As expected, mar-
set of response curves obtained with different
kedly different response was observed for the
intensity of excitation (excitation moment
two horizontal directions tested. The difference
m,e)l under otherwise unchanged conditions.
in resonant frequencies f l and f 2 can be seen
'The excitation moment is given in the figures as from comparison of the cases 3 and 4 in Table
m,ge, lb-in. 4. The vertical amplitude is not shown because
NOVAK AND GRIGG
Can. Geotech. J. Downloaded from www.nrcresearchpress.com by UNIVERSITY OF ARIZONA LIBRARY on 09/06/13
FREQUENCY (CPS)
A
SYMBOL LB-IN.
For personal use only.
FREQUENCY (CPS)
FIG.5. Experimental vertical response of vertically loaded pile ( 1 pile 3.5 in. diameter, mass
of footing = 86.41 slugs). ( a ) dimensional, ( b ) dimensionless.
the resonant peak lay above the frequency in Novak (1974a, 1976) with footing, pile and
range of the oscillator and could not be soil data corresponding to the experiments and
reached. given above. For horizontal excitation two
Finally, the group of four piles was tested shear wave velocities were used: VTV,derived
with cohesive soil compacted into the space from wave propagation, and Vst, obtained from
beneath the pile cap so that full contact was the static tests of the pile.
secured between the test footing and the soil. For vertical excitation two tip conditions
The response curves measured under these were assumed: A fixed tip and a relaxed tip.
otherwise identical conditions are shown in The degree of the relaxation of the tip can be
Fig. 8. The resonant frequencies measured are described by the shear wave velocity of the
given in Table 4. (The excitation load is not soil layer lying below the tip, V,. In addition
given in Fig. 8 because it consisted of both to the relaxation of the tip a reduction of the
forces and moments in this case). shear wave velocity appeared desirable but this
reduction was much smaller than that needed
Comparison of Theory With Experiments for the horizontal excitation. Good results were
Theoretical Response obtained with V, = 2 to 3 V , and V, = 1 to
Theoretical response curves were computed 2 v,.
by means of the general expressions derived An example of the theoretical response to
CAN. GEOTECH. J. VOL. 13, 1976
0 40 SYMBOL LB-IN.
fl>
Can. Geotech. J. Downloaded from www.nrcresearchpress.com by UNIVERSITY OF ARIZONA LIBRARY on 09/06/13
3
a
30
W
0
3
k
_I
n
2 20
(0
ln
W
_I
Z
0
y lo
W
E
0
0
0 10 20 30 40 50 60
FREOUENCY (CPS)
SYMBOL LB-IN
For personal use only.
o 0.855
A 0.598
0 0.427
o 0.171
FREQUENCY (CPS)
FIG.6. Experimental response of horizontaIIy loaded pile (1 pile 2.4 in. diameter, mass of
footing = 3 1.06 slugs). ( a ) horizontal component (sliding), ( b ) rocking component.
vertical as well as horizontal loads is plotted ing constants of individual piles assumed to
in Fig. 9. Comparison of the individual the- act in isolation. However, if the distance be-
oretical and experimental response curves is tween the piles is small, the stiffness and damp-
shown in Figs. 10 to 12. The theoretical nat- ing of each pile in a group are reduced because
ural frequencies and resonant amplitudes are the displacement of one pile contributes also
compared with the experimenta1 values in to the displacements of the other piles. No
Table 4. dynamic solution to this interaction of piles is
available at present but some insight can be
Eflect of Pile Group obtained from the approximate static solution
The total stiffness and damping constants of Poulos (1971). The static solution may
of the pile group used in the analysis were be used for guidance because the behaviour of
established as a sum of the stiffness and damp- a pile group should be similar in both static
NOVAK AND GRIGG
fl SYMBOL L8-IN
2.558
Can. Geotech. J. Downloaded from www.nrcresearchpress.com by UNIVERSITY OF ARIZONA LIBRARY on 09/06/13
0 10 20 30 40 50 60
FREQUENCY (CPS)
For personal use only.
FREQUENCY (CPS )
FIG.7. Experimental response of pile group to horizontal load (4 piles 2.4 in. diameter, mass
of footing = 103.3 slugs). ( a ) horizontal component (sliding), (b) rocking component.
and dynamic cases. The reason for this sim- where a , = the interaction factor describing the
ilarity is that damping is the imaginary part of contribution of the rth pile to the displacement
the complex dynamic stiffness and the dimen- of the reference pile; al = 1 and other factors
sionless frequency a, = r, w / V , , controlling the are smaller than unity and decrease with the
dynamic response, is very low. distance between the piles.
The relationship between the displacement Using [3] the effect of pile interaction upon
u~ of a group of n piles and the displacement the response of the pile group as shown in
ul of a single pile carrying the same load as a Fig. 11 can be examined. The first resonance
pile in the group, can be written as peak, which is of primary importance, is dom-
inated by the first vibration mode and consists
mainly of horizontal translation as can be seen
from Figs. 7a and 7b. Interaction factors CY, for
Can. Geotech. J. Downloaded from www.nrcresearchpress.com by UNIVERSITY OF ARIZONA LIBRARY on 09/06/13
TABLE
4. Natural frequencies and resonant amplitudes
Dimensionless amplitude
Natural frequency
Piles Footing C P (Hz)
~ Vert Horizontal Rocking
diameter m
Case (in) 4 Values fo fi fi Aw A, 1 Auz A* I Aw
For personal use only.
Experimental
Theoretical V,
1 1 pile Theoretical V,,
2.4 Experimental
Theoretical V,
Theoretical V,,
Experimental 4.36
Theoretical V, 8.17
1 pile Theoretical V,, 4 . OO*
2 3.5 Experimental -
Theoretical V , 4.92
Theoretical V,, 9.20
4 piles Experimental
3 2.4 Theoretical V .
x, = 9 Theoretical V,,
4 piles Experimental
4 2.4 Theoretical V,"
x, = 13.5 Theoretical V,,
4 piles
2.4 Experimental
5 x, = 13.5 Theoretical V,"
+ soil Theoretical V,,
*Calculated with relaxed tip and V,, < V, < V,.
NOVAK A N D GRIGG 38 1
Can. Geotech. J. Downloaded from www.nrcresearchpress.com by UNIVERSITY OF ARIZONA LIBRARY on 09/06/13
f" =m=
14'7 12.3 cps
FIG. 8. Influence of soil beneath pile cap upon amplitude of sliding is approximately
experimental response of pile group ( 4 piles 2.4 in.
diameter, mass of footing = 103.3 slugs).
3
0
'Z
26
I O ~
SLIDING
1 \,VERTICAL
where A = the first resonant amplitude ob-
tained without the interaction effect. Equation
6 yields
AG = 9.13 J145= 11.0
2 The first resonant peak approximately adjusted
for the effect of the pile interaction is shown
in dotted lines in Fig. 1 1.
In other vibration modes or with free-head
piles the interaction effect can be expected to
be less pronounced.
The above approach based on summation of
FREQUENCY (CPS) a,. is approximate because it implies that all
FIG.9. Theoretical response of pile group to verti- piles carry the same load.
cal load and horizontal load ( 4 piles 2.4 in. diameter, Discussion
mass of footing = 103.3 slugs, V , = V , , = 147 ft/s,
fixed tip). A number of observations can be made on
the basis of the comparison of the theoretical
sliding fixed-head piles can be obtained from with the experimental results.
Fig. 3 in Poulos (1971). For the conditions of It can be seen that the theory predicts the
the experiment a value of a,. = 1.45 was general character of the response very well.
found using shear wave velocity V,,. The piled foundations respond to changes in
Stiffness of the pile group kG is inversely stiffness, number, and distance between piles
proportional to displacement (flexibility) and and in the footing mass in the manner that the
hence k , = k/C a,.. Therefore, the adjusted theory predicts. In particular, the following
natural frequency of the group is approx- trends predicted theoretically agree with the
imately experiments:
CAN. GEOTECH. J. VOL. 13, 1976
Can. Geotech. J. Downloaded from www.nrcresearchpress.com by UNIVERSITY OF ARIZONA LIBRARY on 09/06/13
For personal use only.
FIG.10. Comparison of theoretical and experimental response ( 1 pile 2.4 in. diameter, mass
of footing = 31.06 slugs; horizontal excitation). (a) sliding component, ( b ) rocking component.
For the single piles subjected to coupled and this in turn increases the natural frequen-
motion, the first resonant peak is larger than cies and reduces the resonant amplitudes.
the second for both sliding and rocking com- The quantitative comparison of the theoret-
ponents. ical predictions with experimental data can be
For the pile group subjected to coupled summarized as follows.
motion, the first resonant peak is dominated by The theory significantly over-estimates the
sliding and the second peak by rocking. The natural frequencies and, in some cases, also the
larger the distance from the center line of the resonant amplitudes if the shear wave velocity
footing to the pile, the greater are the natural derived from wave propagation in deeper
frequencies and the smaller the amplitudes. layers of soil is used. This over-estimation is
This implies that the vertical stiffness should primarily due to the reduction of the shear
not be neglected when considering horizontal modulus with decreasing depth. This reduction
loads as is often done. of the shear modulus can be accounted for,
The effect of soil beneath the pile cap is to approximately, by using the apparent shear
increase both the stiffness and the damping modulus derived from the static test in the
NOVAK AND GRIGG
Can. Geotech. J. Downloaded from www.nrcresearchpress.com by UNIVERSITY OF ARIZONA LIBRARY on 09/06/13
0 j I
0 1",?.3 20 30 40 50 60
FREOUENCY ICPSI
FIG.11. Comparison of theoretical and experimental response (group of 4 piles 2.4 in. diam-
eter, mass of footing = 103.3 slugs, horizontal excitation). Sliding component.
of soil a few pile diameters deep effectively Canada. The assistance of John Howell and
participates in the motion. In the full scale, Gary Lusk is gratefully acknowledged.
this effective layer is deeper because the pile
diameter is greater and the confining pressure KUHLEMEYER, R. L. 1976. Static and dynamic laterally
Can. Geotech. J. Downloaded from www.nrcresearchpress.com by UNIVERSITY OF ARIZONA LIBRARY on 09/06/13
pp. 341-360.
depends very much on the shear modulus of RICHART, F. E., HALL,J. R., and WOODS,R. D. 1970.
soil used in the analysis. With a shear modulus Vibrations of soils and foundations. Prentice-Hall, Inc.,
characteristic of the deeper layers of soil res- New York, N.Y.
onant frequencies and amplitudes can be sig-
nificantly over-estimated. Notations
Better results can be obtained with lower A = area of pile cross section
values of the shear modulus characteristic of = din~ensionlessamplitude of foot-
Ai,j
the uppermost layer of soil. For small displace- ing in direction i and mode j
ments and horizontal excitation, a suitable AG = dimensionless amplitude adjusted
value of this shear modulus can be derived for the effect of pile interaction
from a static test of the pile. This approach a, = r,w/V, = dimensionless frequency
also approximately accounts for the separation c = coefficient of damping
of the pile from the soil. CG = coefficient of damping adjusted
With the shear modulus from the static test, for the effect of pile interaction
natural frequencies are predicted very well and D = damping ratio of footing
the prediction of resonant amplitudes is in DG = damping ratio of footing adjusted
most cases also improved for horizontal ex- for the effect of pile interaction
citation. For pile groups, a further improve- e = eccentricity of rotating mass
ment can be achieved by the correction for the f~ = natural frequency adjusted for the
effect of pile interaction. This correction can effect of pile interaction
be based on the static analysis by Poulos f = natural frequency
(1968, 1971). G = shear modulus of soil
For vertical excitation, the pile tip should I = moment of inertia of pile cross
be considered as relaxed unless it rests on bed- section
rock or the pile is extremely long. Also, a = mass moment of inertia of footing
smaller reduction of the shear modulus appears k = stiffness constant
sufficient for the vertical response. k~ = stiffness constant adjusted for the
effect of pile interaction
Acknowledgements I = length of pile
This study was supported by a research m = mass of footing
grant from the National Research Council of me = unbalanced rotating mass
NOVAK AND GRIGG 385
The velocity of shear waves in the soil was the upper foot or two of the soil. The introduc-
established by means of the steady state vibra- tion of values V,, attempts to describe the soil
tion technique. (The method is described in stiffness in the uppermost layer.