Cap. 4

Properties of Dynamically
Loaded Soils
ID INTRODUCTION
It is a well-known fact that earthquake damage is strongly influenced by the
dynamic properties of local soil deposits. In addition, many problems in civil
engineering practice require the knowledge of the properties of soils subjected
to dynamic loading. These problems include the dynamic bearing capacity of
foundations, response of machine foundations subjected to cyclic loading, soil
structure interaction during the propagation of stress waves generated due to an
earthquake, and earthquake resistance of dams and embankments.
A variety of laboratory tests as well as field techniques are available, each hav
ing its own limitations as well as advantages. Some of these tests are specifically
developed for measuring properties of dynamically loaded soils whereas some
are modified versions of tests used in the domain of traditional soil mechan
ics. Some of these methods are suitable for small strain range whereas some are
suitable for large strain range. The range of strain of interest usually dictates the
type of equipment/method to be used which in turn depends on the problem to
be analyzed at hand. Some of these equipments are very specialized, expensive,
and require special training to use and interpret the results. It is worth noting
that soil behavior over a wide range of strains is nonlinear and, on unloading,
follows a different stress-strain path forming a hysteresis loop.
This chapter is devoted primarily to describing various laboratory and field
test procedures available to measure as well as estimate the soil properties using
empirical correlations subjected to dynamic loading. This chapter is divided into
three major parts:
a. Laboratory tests and results
b. Field tests and measurements
c. Empirical correlations for the shear modulus and damping ratio obtained
from field and laboratory tests. These are the two most important parame
ters needed for most design work.
Copyright 2017 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eCha 1),
Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restriction ieq11im it
4 I Properties of Dynamically Loaded Soils
Laboratory Tests and Results
lfl SHEAR STRENGTH OF SOILS UNDER RAPID

LOADING CONDITIONS
Saturated Clay
In most common soil test programs, the undrained shear strength of saturated
cohesive soils is determined by conducting unconsolidated-undrained triaxial
tests. The soil specimen for this type of test is initially subjected to a confining
pressure a 3 in a triaxial test chamber, as shown in Figure 4.la. After that an axial
stress �a is applied to the specimen (Figure 4.1 b). The axial stress �a is gradu
ally increased from zero to higher values at a constant rate of compressive strain.
The strain rate e is maintained at about 0.5% or less. The general nature of �a
versus axial strain£ diagram thus obtained is shown in Figure 4.lc. The total
major and minor principal stresses at failure can now be given as:
Major principal stress (total)= a 1u) = a 3 + �a max
Minor principal stress (total)= a3
The total stress Mohr's circle at failure is shown in Figure 4.ld. It can be
shown (see Das, 1990) that for a given saturated clayey soil, the magnitude
of �a max is practically independent of the confining pressure a 3, as shown in
Figure 4.le. The total stress Mohr's envelope for this case is parallel to the nor
mal stress axis and is referred to as the </J = 0 condition (where </J = angle of shear
ing resistance of the soil). The undrained shear strength cu is expressed as
�a max al ( f) -a3
Cu = = (4.1)
2 2
The undrained shear strength obtained by conducting tests at such low-axial
strain rates is representative of the static loading condition, or Cu = Cu(static)· Experimen
tal results have shown that the magnitude of �a max = a 1u) -a 3 gradually increases
with the increase of axial strain rate e. This conclusion can be seen from the labora
tory test results on Buckshot clay (Figure 4.2). From Figure 4.2 , it can be observed
that Cu = �a max /2 = (a 1cn -a 3 )/2 obtained between strain rates of 50% to 425%
are not too different and can be approximated to be a single value (Carroll, 1963).
This value can be referred to as the dynamic undrained shear strength, or
Cu = Cu( dynamic)
Carroll suggested that for most practical cases, one can assume that
Cu( dynamic) ::::::: 1 .S (4.2)
Cu (static )
Copyright 2017 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter{s).
Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restrictions require it.
4.2 Shear Strength of Soils Under Rapid Loading Conditions
(a) (b)
15'
I
6
II
b
<I
Axial strain,£
(c)
Total stress Mohr's

envelope(¢= 0 concept)
Cu
1 ._______.______._____._______.____
0:, Oi(f)
Normal
f-Aamax --+I stress
(d) (e)
Figure 4.1 Unconsolidated undrained triaxial tests
Sand
GRUPO A /22/04/21
Several vacuum triaxial test results on different dry sands (that is, standard
Ottawa sand, Fort Peck sand, and Camp Cooke sand) were reported by Whitman
and Healy (1963). These tests were conducted with various effective confining
200 ,--------..---------�-----------,
Moisture content= 33.5 ± 0.2%

u3 = 100 kPa
150
i= 314%/s and 426%/s
128%/s
53.6%/s
,._., 14.4%/s
0 100
I
0 4.76%/s
0.033%/s
50
OL---------'--------......________,
4 8 12
Strain, c(¾)
Figure 4.2 Unconsolidated-undrained triaxial test results on Buckshot clay

(after Carroll, 1963)
Source: Carroll, WP. (1963). "Dynamic Bearing Capacity of Soils. Vertical Displacements of
Spread Footings on Clay: Statics and Impulsive Loadings." Technical Report No. 3-599, U.S.
Army Corps of Engineers, Waterways Experiment Station, Vicksburg, Mississippi.
pressures (a3) and axial strain rates. The compressive strength Llamax determined
from these tests can be given as
(4.3)
where a3 = effective minor principle stress
ii1u) = effective major principal stress at failure
An example of the effect of axial strain rate on dry Ottawa sand is shown in
Figure 4.3. It can be seen that for a given a3 the magnitude of LlCJmax decreases
initially with the increase of the strain rate to a minimum value and increases
thereafter. From fundamentals of soil mechanics it is known that
Copyright 2017 Cengage Learning. All Rights Reserved. May not be copied. scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter{s).
4.3 Strength and Deformation Characteristics of Soils under Transient Load
All data corrected to void ratio of 0.52
r
(Todos los datos corregidos a una proporción de vacíos de 0,52)
210
a3 = 88.3kPa
I
k
'-"'
It:)
I
Range of scatter
10 140
(Rango de dispersión)
11)
rn
.....rn a3 = 33.8kPa
1
11)
1-
I
11) 70
0
0.01 0.1 1 10 100 1000
Rate of strain, &(%/s)
Figure 4.3 Strain-rate effect for dry Ottawa sand (after Whitman and Healy, 1963)
Source: Whitman, R.V., and Healy, K.A. (1963). "Shear Strength of Sands During Rapid
Loadings," Transactions, ASCE, Vol. 128, Part 1, pp. 1553-1594. With permission from ASCE.
q, = sin- 1(�1U) - �3 J (4.4)

0'1(1) + 0'3
where q, = drained soil friction angle
Based on Figure 4.3 and Eq. (4.4) it is obvious that the initial increase of the
strain rate results in a decrease of the soil friction angle. The minimum dynamic
friction angle may be given as (Vesic, 1973)
(4.5)
(obtained from static tests-

that is, small strain rate of loading )
Ill
( Obtenido a partir de pruebas estáticas, esto es, una pequeña tasa de deformación de cargado)
STRENGTH AND DEFORMATION CHARACTERISTICS

OF SOILS UNDER TRANSIENT LOAD
In many circumstances it may be necessary to know the strength and deforma
tion characteristics of soils under transient loading. A typical example of tran
sient loading is that occurring due to a blast. Figure 4.4 shows the nature of an
idealized load versus time variation for such a case. In this figure, QP is the peak
load, tL is the time of loading, and tn is the time of decay.
1----tL--------tn-----
Figure 4.4 Transient load
Casagrande and Shannon (1949) conducted some early investigations to

study the stress-deformation and strength characteristics of Manchester sand
and Cambridge clay soils. Undrained tests were conducted in three specially
devised apparatuses-one falling-beam apparatus and two pendulum-loading
apparatuses. In these specially devised pieces of equipments, the loading pattern
on soil specimens was similar to that shown in Figure 4.4. Figure 4.5a shows the
variation of stress and strain with time for an unconfined Cambridge clay spec
imen with tL = 0.02 s. Similarly, Figure 4.5b compares the nature of variation
of strain versus stress for static and transient (tL = 0.02 s) loading conditions on
unconfined Cambridge clay specimens. The unconfined compressive strength
determined in this manner with varying times of loading is shown in Figure 4.6.
Based on Figures 4.5b and 4.6, the following conclusions may be drawn.
1. qu(transient) :::::::
1.5 to 2
qu(static)
where qu = unconfined compression strength. This is consistent with the

findings of Carroll (1963) discussed in Section 4.2.
2. The modulus of deformation E as defined in Figure 4.7 is about two times as
great for transient loading as compared to that for static loading.
The nature of the stress-versus-strain plot for confined compression tests on
Manchester sand conducted by Casagrande and Shannon (1949) is as shown in
Figure 4.8. From this study it was concluded that
1. [o\n - <Y3]1ransient ::::::: 1.l and

[ (J 1(/) - G3 lstatic
2. The modulus of deformation as defined by Figure 4.7 is approximately the
same for transient and static loading conditions.
4.4 Travel-Time Test for Determination of Longitudinal and Shear Wave Velocities (u c and u5)
400 �--�-----�--�------ 8
300 6
,, ---'----- Shear failure at tL = 0.02s ,-.._

�
' 4
I \
' _,
';:; 200
''' ·a
.. ......
rJl
.....� ' r/J.

r/J.
'
'
100 ' \,
../
Stress
2
.. ..
- ..
.. -
o----....____.___.._._-_-_.._-_-_-_.__- _-________.._____.o
0 0.04 0.08 0.12 0.16 0.20 0.24
Time (s)
�:::::::::::
oiJ ...... ::r:::�==-1
(a)
0 Ir--�..� .... ..... .., ..

Static test time Shear failure --�
2 loading, 456 s ..
" /
I
I
I
I
I
I
- - _,
I
I
Transient
loading,
6 tL =0.02s
8�--�--�--�--�---�--�
40 80 120 160 200 240
Stress (kPa)
(b)
Figure 4.5 Unconfined compressive strength of Cambridge clay for varying

time of loading (after Casagrande and Shannon, 1949)
Source: Cas agrande, A., and Shannon, W.L. (1949). "Strength of Soils Under Dynamic Loads,"
Transactions, ASCE, Vol. 114, pp. 755-772. With permission from ASCE .
Ill TRAVEL-TIME TEST FOR DETERMINATION OF

LONGITUDINAL AND SHEAR WAVE VELOCITIES
(u c AND U s )
Using electronic equipment, the time tc requirement for travel of elastic waves
through a soil specimen of length L can be measured in the laboratory. For
longitudinal wave
(4.6)
260----�------------------
� 240 t-----1-----i-----------r-------------,:c-a-,,i
•
80._________._______.________._______._________,
1000 100 10 1 0.1 0.01
Time ofloading (s)
Figure 4.6 Unconfined compressive strength of Cambridge clay for varying

time of loading (after Casagrande and Shannon, 1949)
Source: Casagrande, A., and Shannon, W.L. (1949). "Strength of Soils Under Dynamic Loads,"
Transactions, ASCE, Vol. 114, pp. 755-772. With permission from ASCE.
E I
Strain
Figure 4.7 Definition of modulus of deformation, E

4.4 Travel-Time Test for Determination of Longitudinal and Shear Wave Velocities (u c and u5)
Deviator stress, a1 - a3
Soil
specimen
Figure 4.8 Confined-compression test on sand-stress- versus-strain behavior

under static and transient loading
The modulus of elasticity E can then be calculated from Eq. (3.30) as
or
L2
E = pv; = p-2- (4.7)
tc
If the soil specimen is confined laterally, then the travel time will give the
value of v'c as shown in Eq. (3.35). Thus Ve = L/t' c , and
L2
M = p- (2
(4.8)
C
where t' c = time of travel of longitudinal waves in a laterally confined specimen

M = constrained modulus
Similarly, if the travel time ts for torsional waves through a soil of length L is
determined, the velocity vs can be given as vs= L/tc , and
L2
G = pv; = p2 (4.9)
ts
Whitman and Lawrence (1963) have provided limited test results for v'c in
20--30 Ottawa sand. The schematic diagram of the apparatus used for measuring
v'c is shown in Figure 4.9a. The soil specimen was confined in 76.2 mm diam
eter Shelby tube (Figure 4.9b). Vertical load was applied by an aluminum pis
ton. In this system, a pulse was sent from one piezo-electric crystal and received
Pulse
�-� (Trigger) ,'----_--;__(P
_ _ u�ls_e _out
_ _ )�
Time
source '--------+-'
(1000 V DC)
(Trigger) Specimen
Power container
supply
Oscilloscope �---------�
(Received pulse)
(a)
------Vertical load
L--ft:=�S<=,;:::: ---Pulse in
Barium titanate ----Aluminum piston
sending transducer
l\---::�-----Lucite facing
_,___Specimen
-�---Shelby tubing (76.2 mm diameter)
-it-,--1--- -Ultrasonic waves
Barium titanate ------

receiving crystal ----Aluminum piston
--Ls�f-c::::::::-+-- Received pulse

to oscilloscope
(b)
Figure 4.9 Travel-time method: (a) schematic diagram of the laboratory

setup for measuring v'c ; (b) details of the soil specimen and container for the
laboratory setup (after Whitman and Lawrence, 1963)
Source: Whitman, R.V., and Lawrence, F.V. (1963). "Discussion on Elastic Wave Velocities
on Granular Soils," Journal of the Soil Mechanics and Foundations Division, ASCE, Vol. 89,
No. SM5, pp. 112-118. With permission from ASCE.
by a second one at the opposite end. The received signal was displayed on an
oscilloscope, which allowed measurement of t'c · It was found that the velocity v'c
increases with the increase of axial pressure.
Ill BENDER ELEMENT TEST FOR DETERMINATION

OF SHEAR WAVE VELOCITY (v 5 ) ( ASTM D 2845-08)
The bender element test is now commonly used in the laboratory to evaluate the
small-strain shear modulus (Gmax ) by measurement of the shear wave velocity
4.5 Bender Element Test for Determination of Shear Wave Velocity (u5)
(Vs ) propagation in a soil specimen. The bender element is a thin transducer that
consists of two piezoelectric plates rigidly bonded to a center shim of brass or
stainless steel plate. The bender element can be implemented in most common
laboratory soil research equipment including unconfined compression, oedome
ter, direct shear, triaxial, resonant column, centrifuges, and so forth. Two in-line
bender elements as a pair are usually installed at both ends of a soil specimen,
acting as transmitter and receiver, respectively. Figure 4.10a shows a pair of
bender elements implemented in a triaxial cell, and Figure 4.1 Ob shows a soil
specimen with a pair of bender elements.
Various waveforms such as sinusoidal and square waves can be used as an exci
tation signal. The bender element can convert between electrical voltage and mechan
ical excitation/bending motion. The wave form of the input voltage yields a bending
motion in the transmitter element, which produces a shear wave propagating through
the specimen. When the receiver element at the other end of specimen is bent by the
arrival of shear wave, an electrical signal is generated in a wave form. The transmit
ted and received wave forms can be captured and displayed by a digital oscilloscope
to determine the shear wave travel time Cts ). The shear wave velocity is calculated as
=
L
Vs -
ts
(a) (b)
Figure 4.10 Test with bender element: (a) a pair of bender elements installed
in a triaxial cell; (b) soil specimen with the pair of bender elements (Courtesy of
Geocomp Corporation, Acton, Massachusetts)
Source: Courtesy of Geocomp Corporation, Acton, Massachusetts
where L is the travel distance and can be determined as the distance between
the two tips of the two bender elements (Dyvik and Madshus, 1985). The small
strain shear modulus can be determined from elastic wave theory [similar to
Eq. (4.9)]:
G =Gmax = pvf = P( � )'
where p is the mass density of soil. The Gmax measured by bender element can
generally have good agreement with that obtained from the resonant column test
(Section 4.6). The advantages of bender element are: (a) the test procedure of
bender element is simple and efficient; (b) the bender element test is nondestruc
tive; and (c) the number of specimens required is minimal.
Ill RESONANT COLUMN TEST

The resonant column test essentially consists of a soil column that is excited to
vibrate in one of its natural modes. Once the frequency at resonance is known,
the wave velocity can easily be determined. The soil column in the resonant col
umn device can be excited longitudinally or torsionally, yielding velocities of Ve
or vs , respectively. The resonant column technique was first applied to testing
of soils in Japan by Ishimato and Iida (1937) and Iida (1938, 1940). Since then
it has been extensively used in many countries, with several modifications using
different end conditions to constrain the specimen. One of the earlier types of
resonant column devices in the United States was used by Wilson and Dietrich
(1960) for testing clay specimens.
Hardin and Richart (1963) reported the use of two types of resonant col
umn devices-one for longitudinal vibration and the other for torsional vibra
tion. The specimen were free at each end (free-free end condition). A schematic
diagram of the laboratory experimental setup is shown in Figure 4.11. The
power supply and amplifier No. 1 were used to amplify the sinusoidal output
signal of the oscillator, which had a frequency range of 5 Hz to 600,000 Hz.
The amplified signals were fed into the driver, producing the desired vibra
tions. Figure 4.12a shows the schematic diagram of the driver for torsional
oscillation. Similarly, the schematic diagram of the driver for longitudinal
vibration is shown in Figure 4.12b. These devices will give results for low-am
plitude vibration conditions. With free-free end conditions, for longitudinal
vibrations at resonance
Wn L
vc =
-- (3.53)
n'Tt
For n = l (that is, normal mode of vibration),
VC
= Wn L = 2'1tfn L = 2/, L
n
'1t '1t
4.6 Resonant Column Test
: Pickup : Amplifier
Oscillator
- - - � � ]� � � - - J No.2
Power
supply
Amplifier
No.I
--- _'::[::- - , Oscilloscope
Driver
Pressure cell
Figure 4.11 Schematic diagram of experimental setup for resonant column test
of Hardin and Richart-free-free end condition
Side view
Side view
Specimen
Specimen
Electromagnet
Rubber
Permanent magnet
Electromagnet
Bearing---1--.--r-+-�
Frame Permanent
magnet
Frame
Rubber
Top view Top view

(a ) ( b)
Figure 4.12 Drawings for steady-system vibration drivers in the resonant

column devices with free-free end conditions: (a) for torsional vibration; (b) for
longitudinal vibration (after Hardin and Richart, 1963)
Source: Hardin, B. 0., and Richard, F. E., Jr. (1963). "Elastic Wave Velocity in Granular Soils,"
Journal of the Soil Mechanics and Foundations Division, ASCE, Vol. 89, No. SMl, pp. 33-65.
With permission from ASCE.
or
or
(4.10)
Similarly, for torsional vibration, at resonance (with n = 1)

Vs =2fn L
t
or
v, = =2fn L
or
(4.11)
Once the magnitudes of E and G are known, the value of the Poisson's ratio can
be obtained as
E
(4.12) =
µ 2G -l
Hall and Richart (1963) also used two other types of resonant column devices
(one for longitudinal vibration and the other for torsional vibration). The end
conditions for these two types of devices were fixed-free-fixed at the bottom and
free at the top of the specimen. The general layouts of the laboratory setup for
this equipment were almost the same as shown in Figure 4.11, except for the fact
that the driver and the pickup were located at the top of the specimen. This is
shown in Figure 4.13. Since the driver and the pickup were located close together,
a correction circuit was introduced to correct the inductive coupling between the
driver and the pickup. The driver and pickup were attached to a common frame.
The differences in construction and arrangement of the driver and the pickup
produce either longitudinal or torsional vibration of the specimen.
fixed free
DE AQUÌ A. Derivation of Expressions for Ve and f for Use in the Fixed-Free-Type

Resonant Column Test
GRUPO/27/04/21
An equation for the circular natural frequency for the longitudinal vibration of
short rods withfixed-free end conditions was derived in Eq. (3.57) as
(2n - l)1t Ve
(J)n =-----
2 L
However, in a fixed-free type resonant column test, the driving mechanism and
also the motion-monitoring device have to be attached to the top of the speci
men (Figure 4.14), in effect changing the boundary conditions assumed in deriv
ing Eq. (3.57). So a modified equation for the circular natural frequency needs to
be derived. This can be done as follows.
Top view
Brass rods- -
--- -
Permanent magnet
Permanent magnet
Driving coil Pickup coil
4.2 mm brass frame---�

for holding coils
(a)
Side view �----Brass rod
Frame
Soft steel
Permanent magnet
(b)
Figure 4.13 Driving and measuring components for a fixed-free resonant

column device (after Hall and Richart, 1963)
Source: Hall, J.R., Jr., and Richard, F.E., Jr. (1963). "Dissipation on Elastic Wave Energy in
Granular Soils," Jour nal of the Soil Mechanics and Foundations Division, ASCE, Vol. 89,
No. SM6, pp. 27-56. With permission from ASCE.
! �
Mass=m
--
,�
Soil
column�
mass=ms �
I •
Figure 4.14 Derivation of Eq. (4.20)
Let the mass of the attachments placed on the specimen be equal to m. For
the vibration of the soil column in a natural mode,
(3.41)
and
(3.43)
At x = 0, U(x) = 0. So B 2 in Eq. (3.43) is zero. Thus
U(x) = B 1 sin( rot J (4.13)
At x = L, the inertia force of mass mis acting on the soil column, and this can
be expressed as
a2 u
F=-m (4.14)
at 2
where F = inertia force. Also, the strain
au F
(4.15)
ax AE
where A = cross-sectional area of specimen

E = modulus of elasticity.
Combining Eqs. (3.41), (4.13), and (4.15) we get
F =
au (au) (A1 sin
= OJn t + A2 cos OJn t)
AE ax ax
= � [ B, sin(
a
::x )] (Ai sin OJ nt + A2 cosw.t) (4.16)
= ( B�:· )[ cos(rot)] CAi sinw.t + A 2 cos OJ ntl
Again, combining Eqs. (3.41), (4.13), and (4.14),
F =-m �:� =-m[B, sin(rot)]( :t 2) CAi sin . t + A2 coswnt)

2
ro
OJ
= m w?; B , sin ( : x) (A, sin OJ nt + A2 cosw.t) (4.17)
,
Now, from Eqs. (4.16) and (4.17),
(4.18)
Atx = L
AE = mOJ n v, tan( OJ;,L) (4.19)
However, Ve =
�Elp; or E = v;p. Substitution of this in Eq. (4.19) gives
mnL
A Ve2 p _
- m OJn Ve tan ( )
�
Ap mn mnL
= tan( )
m Ve Ve
ALp mnL mnL

= tan( )
m Ve Ve
or
IA�y =atanal (4.20)
where r =pg= unit weight of soil

W =mg= weight of the attachments on top of the specimen
Mass=m
�
Polar moment of
mertia=Jm
--
Soil column
�
Polar moment of L
inertia= Js
X
'
I I
Figure 4.15 Derivation of Eq. (4.23)
where ls = mass polar moment of inertia of the soil specimen and Jm = mass
polar moment of inertia of the attachments with mass m.
Thus
_ mnL _ 2rcj�L
Vs------- (4.24)
a a
and
2
G = pv} = 39.48 ( /�; )p (4.25)
C. Recent Developments in Resonant Column Test Equipment

Figure 4.16 shows a recently developed fully automated resonant column assem
bly with a quasi-fixed base and free top (Werden et al., 2013). Unlike the fixed
free end condition shown in Figure 4.15, Figure 4.17 shows a schematic diagram
for the quasi-fixed base and free-top condition. This type of resonant column
device has a torque transducer attached between the specimen's bottom plate
and the fixed base. The actual torque transmitted through the specimen can be
directly measured by the torque transducer.
Figure 4.16 Photograph of a fully automated resonant column assembly

(Courtesy of Geocomp Corporation, Acton, Massachusetts).
Source: Courtesy of Geocomp Corporation, Acton, Massachusetts
SIGUE
D. Typical Laboratory Test Results from Resonant Column Tests
Most of the laboratory test results obtained from resonant column tests are for
low amplitudes of vibration. Low amplitudes of vibration mean strain ampli
tudes of the order of 10-4 or less.
Typical values of Ve and Vs with low amplitudes of vibration for No. 20-30
Ottawa sand compacted at a void ratio of about 0.55 are shown in Figures 3.13
and 3.14. These were conducted using the free-free and fixed-free types of
resonant column device developed by Hardin and Richart (1963) and Hall and
Richart (1963). Based on the results given in Figures 3.13 and 3.14, the following
general conclusions can be drawn:
1. The values of Ve and Vs in soils increase with the increase of the effective aver
age confining pressure iiO•
Legend
1. Device to apply sinusoidal torque
2 2. Soil specimen
3. Passive end platen
4. Torque transducer
5. Fixed base
Figure 4.17 Schematic diagram showing quasi-fixed base and free-top

condition for the soil specimen
2. The values of Ve and Vs for saturated soils are slightly lower than those for dry
soils. This can be accounted for by the increase of the unit weight of soil due
to the presence of water in the void spaces.
Hardin and Richart (1963) also reported the results of several resonant col
umn tests conducted in dry Ottawa sand. The shear wave velocities determined
from these tests are shown in Figure 4.18. The peak-to-peak shear strain ampli
tude for these tests was 10- 3 rad. From Figure 4.18, it may be seen that the values
of Vs are independent of the gradation, grain-size distribution, and also the rel
ative density of compaction. However, vs is dependent on the void ratio and the
effective confining pressure.
E. Shear Modulus for Large Strain Amplitudes

For solid cylindrical specimens torsionally excited by resonant column devices,
the shear strain varies from zero at the center to a maximum at the periphery, and
it is difficult to evaluate a representative strain. For that reason, hollow cylin
drical soil specimens in a resonant column device (Drnevich, Hall, and Richart,
1966, 1967) may be used to determine the shear modulus and damping at large
strain amplitudes. Figure 4.19 shows a schematic diagram of this type of appa
ratus, in which the average shearing strain in the soil specimen is not greatly
different from the maximum to the minimum. The variation of the shear mod
ulus of dense C-190 Ottawa sand with the shear strain amplitude y' is shown in
Figure 4.20. Note that the value of G decreases with r', but it decreases more
Soil
• No. 20-No. 30
o No. 80-No. 140
The peak-to-peak shear strain A 74.8% of No. 20-No. 30 and 25.2% of No. 80-No. 140
amplitude for these tests was 10-3 rad .6. No. 20-No. 140-well graded
390.-----....,.....-------,-----...-----....,....---,
360 Effective confining

pres re ( <To)
:
7
300 kPa
180 A
120 .________.________.____....________._____,
0.35 0.45 0.55 0.65 0.75
Void ratio, e
Figure 4.18 Variation of shear wave velocity with effective confining pressure
CJ0 for round-grained dry Ottawa sand (after Hardin and Richart, 1963)
Source: Hardin, B. 0., and Richard, F. E., Jr. (1963). "Elastic Wave Velocity in Granular Soils,"
Journal of the Soil Mechanics and Foundations Division, ASCE, Vol. 89, No. SMl, pp. 33-65.
With permission from ASCE.
rapidly for y' > 10-4_ This is true for all soils. The reason for this can be explained
by the use of Figure 4.21, which is shear-stress-versus-strain diagram for a soil.
The stress-strain relationships of soils are curvilinear. The shear modulus that
is experimentally determined is the secant modulus obtained by joining the
extreme points on the hysteresis loop. Note that when the amplitude of strain is
small (that is, y' = y'i; Figure 4.21), the value of G is larger compared to that for
the larger strain level (that is, r' = r'2 ).
F. Effect on Prestraining on the Shear Modulus of Soils

The effect of shear modulus of soils due to prestraining was reported by
Drnevich, Hall, and Richart (1967). These tests were conducted using C-190
Ottawa sand specimens. The specimens were first vibrated at a large amplitude
for a certain number of cycles under a constant effective confining pressure (<JO ).
Taring spring-----------____ LongLVDT
Rot. accel
Vert. accel. Rot.
LVDT
Drive coils
Top cap
... ..
..... ..��.
, ..·
...
?=
.....
Sand
.. .
•.:-,..
��
."'.
Pore pressure transducer
To back pressure� Bottom cap
Figure 4.19 Schematic diagram of hollow-specimen resonant column device

(after Drnevich, 1972).
Source: Dmevich, V.P. (1972). "Undrained Cyclic Shear of Saturated Sand," Journal of the Soil Mechanics
and Foundations Division, ASCE, Vol. 98, No. SM8, pp. 807-825. With permission from ASCE.
175.--------------------,
140
105
C!,
;:::l
;:::l
70
35
C-190 Ottawa sand
Void ratio = 0.46
0
0.1 1.0 7.0
Shear strain, r' (rad x 104)
Figure 4.20 Effect of strain amplitude on shear modulus of sand (after
Drnevich, Hall, and Richart, 1967)
Source: Dmevich, V.P., Hail, J.R., Jr., and Richart, F.E. Jr. (1967). "Effects of the Amplitude of
Vibration on the Shear Modulus of Sand," Proceedings, International Symposium on Wave Propagation
and Dynamic Properties of Earth Materials, University of New Mexico Press, pp. 189-199. 117
--r2
Shear strain, r'
Figure 4.21 Nature of variation of shear stress versus shear strain
After that the shear modulii were determined by torsionally vibrating the speci
mens at small amplitudes (shearing strain< 10- 5). Figure 4.22 shows the results
of six series of this type of test for dense sand (void ratio= 0.46). In general, the
value of G increases with increase of prestrain cycles.
Sigue
G. Determination of Internal Damping
In Section 3.15, a distinction was made between internal damping and material
damping. The internal damping of a soil specimen can be determined by reso
nant column tests.
In Chapter 2, the derivation of the expression for the logarithmic decrement
was given as
21tD
c5 = In Xn = (2 .70)
Xn+1 1- D 2
where 8 = logarithmic decrement
D = damping ratio.
The preceding equation is for the case of free vibration of a mass-spring
dashpot system. The damping ratio is given by the expression
375 .----,------------.---------,.----------,--------,
C-190 Ottawa sand
Void ratio= 0.46
a0= 56kPa
'b
-. 300
.6.
y'= 6 X 10
-4
X
�
0
C!,
,.,r 225
;::l
"O
0
a0 = 30 kPa
r;5 150 ��56kPa
-4
y'= 1.6 X 10
�
a0 = 30 kPa -4
_ _____ ______ ....._
_____ ____. y' = 1.6 X 10
_____
75 .___ ......._ ___._ _____.
10 3 10 4 10 5 10 6 107
Cycles of high amplitude torsional vibration
Figure 4.22 Effect of number of cycle of high-amplitude vibration on shear

modulus determined at low amplitude (after Drnevich, Hall, and Richart, 1967)
Source: Dmevich, V.P., Hail, J.R., Jr., and Richart, F.E. Jr. (1967). "Effects of the Amplitude of
Vibration on the Shear Modulus of Sand," Proceedings, International Symposium on Wave Propagation
and Dynamic Properties of Earth Materials, University of New Mexico Press, pp. 189-199.
- C C
D (1.47b)
- Ccr -
- 2.Jkms
where ms= mass of the soil specimen (in this case).

For soils, the value of D is small and Eq. (2. 70) can be approximated as
Xn =
u� = 1n-- 2rcD (4.26)
Xn+1
Now, combining Eqs. (1.47b) and (4.26)

rec
8= (4.27)
,Jkms
The logarithmic decrement of a soil specimen (and hence the damping ratio D)
can easily be measured by using a fixed-free type resonant column device.
The soil specimen is first set into steady-state forced vibration. The driving
power is then shut off, and the decay of the amplitude of vibration is plotted
against the corresponding number of cycles. This plots as a straight line on a
semilogarithmic graph paper, as shown in Figure 4.23. The logarithmic decre
ment can then be evaluated as
Ounoorra;ed = C) ( 1:) In (4.28)
.9 "uncorrected = n 1 n
1 X.
Xo
n
.....
(I)
Xo
.....
n Number of cycles (log scale)
Figure 4.23 Plot of the amplitude of vibration against the corresponding

number of cycles for determination of logarithmic decrement
However, in a fixed-free type of resonant column device, the driving and

the motion-monitoring equipment is placed on the top of the specimen. Hence,
for determination of the true logarithmic decrement of the soil specimen, a cor
rection to Eq. (4.28) is necessary. This has been discussed by Hall and Richart
(1963). Consider the case of longitudinal vibration of a soil column, as shown in
Figure 4.24, in which m = mass of the attachments on the top of the soil speci
men and ms = mass of the soil specimen. With the addition of mass m, Eq. (4.27)
can be modified as
1tC
8uncorrected = �k (4.29)
( ms
From Eqs. (4.27) and (4.29),
__ 8_=�m, +m =�l+-m (4.30)

8uncorrected ms ms
In order to use Eq. (4.30), it will be required to convert the mass ms into an
equivalent concentrated mass. The equivalent concentrated mass can be shown
to be equal to 0.405ms. Thus, replacing ms in Equation (4.30) by 0.405ms,
m
8 = Duncorrected 1 + --- (4.31)
0.405ms
A similar correction may be used for specimens subjected to torsional vibration,
which will be of the form
4.7 Cyclic Simple Shear Test
Mass=m
t
Mass =m8
C k
Figure 4.24 Fixed-free soil column
8 = Duncorrected
Hardin (1965) suggested a relation for 8 of dry sand in low amplitude tor
sional vibration as
8 = 9n(r')0.2(<Yo)-0.5 (4.32)
Equation (4.32) is valid for y' = 10-6 to 10-4 and <Yo = 24 kPa to 144 kPa.
PASA A LA PÁGINA 76 ♪
Ill CYCLIC SIMPLE SHEAR TEST
A cyclic simple shear test is a convenient method for determining the shear mod
ulus and damping ratio of soils. It is also a convenient device for studying the
liquefaction parameters of saturated cohesionless soils (Chapter 10). In cyclic
simple shear tests a soil specimen, usually 20-30 mm high with a side length (or
diameter) of 60-80 mm, is subjected to a vertical effective stress a v and a cyclic
shear stress r, as shown in Figure 4.25. The horizontal load necessary to deform
the specimen is measured by a load cell, and the shear deformation of the speci
men is measured by a linear variable differential transformer.
.... ··:·.: ·:·_:\·_:. :·.: ·:.· :::·.· ... :·.:·: ·::::··· :._ ... .
·.· .· . . ·.·. · ·
. .· . . .·. ·.·.·.· .
. . .·::
. ·.·
.· .. · -:·: : · :·_:: . .·::· /-:·:· './:·: .
· ·. .
·:·. · ·. :· · .: . :·.. : ·: _:.·.:·: :..-.::·· ·::· .
:
:: . ·
::.-
: .
:.
· : ·. : ·
:: . .
. .
Figure 4.25 Cyclic·:· simple shear test
The shear modulus of a soil in the cyclic simple shear test can be
determined as
amplitude of cyclic shear stress, r
G= (4.33)
amplitude of cyclic shear strain, r'
The damping ratio at a given shear strain amplitude can be obtained from the
hysteretic stress-strain properties. Referring to Figure 4.26 (also see Figure 4.21),
the damping ratio can be given as
l area of the hysteresis loop
D = __ (4_34)
21t area of traingle OAB and OA' B'
Shear stress, r
Shear strain, y'
Figure 4.26 Determination of damping ratio from hysteresis loop [Eq. (4.34)]
7.5
Uv = 200 kPa l
Rn=60%
X
c,;l 5.0
'-' Cycle 300
C!)
er.,�
;::::l
10
;::::l
s
<l'v =25kPa l
2.5
0
0.001 0.01 0.1 1.0
Cyclic shear strain, y' (%)
Figure 4.27 Shear modulus-shear strain relationship for medium dense sand
(after Silver and Seed, 1971)
Source: Silver, M.L ., and Seed, H.B. (1971). "Deformation Characteristics of Sands Under Cyclic
Loading," Journal of the Soil Mechanics and Foundation Division, ASCE, Vol. 94, No. SM8,
pp. 1081-1098. With permission from ASCE.
Figure 4.27 shows a plot of shear modulus G with cyclic shear strain y' for
two values of av (Silver and Seed, 1971) obtained from cyclic simple shear tests
on a medium dense sand (relative density, Rn= 60%). From the results of this
study, the following can be stated:
1. For a given value of y' and av , the shear modulus increases with the number
of cycles of shear stress application. Most of the increase in G takes place in
the first ten cycles, after which the rate of increase is relatively small.
2. For a given value of ifv and number of cycles of stress application, the mag
nitude of G decreases with the amplitude of shear strain y'. (Note: Similar
results are shown in Figure 4.20.)
3. For a given value of y' and number of cycles, the magnitude of G increases
with the increases of ifv.
The nature of the shear-stress-versus-shear-strain behavior of a dense sand
under cyclic loading is shown in F igure 4.28. Using the hysteresis loops of this
type and Eq. (4.34), the damping ratios obtained from a cyclic simple shear test
for a medium dense sand are shown in Figure 4.29. Note the following:
1. For a given value of ifv and amplitude of shear strain y', the damping ratio
decreases with the number of cycles. Since, in most seismic events, the number of
significant cycles is likely to be less than 20 (Chapter 7), the values determined at
5 cycles are likely to provide reasonable values for all practical purposes.
2. For a given number of cycles and av , the magnitude of D decreases with the
decrease of y'.
r(kPa) r(kPa)
5 5
Cycle 1 Cycle 2
r'(¾) r'(¾)
-0.05 0.05 -0.05 0.05
-5 -5 Rn = 80%
o'v = 25 kPa
r(kPa) r(kPa)
5 5
Cycle 10 Cycle 300
-1-------,,..._,____-+- r'(¾) r'(¾)

-0.05 0.05 -0.05 0.05
-5 -5
Figure 4.28 Stress-strain behavior of dense sand under cyclic shear (after Silver
and Seed, 1971)
Source: Silver, M.L., and Seed, H.B. (1971). "Deformation Characteristics of Sands Under Cyclic
40
Rn =60%
<Tv = 25 kPa
,...._ 30
�
r' =o.63%
·-
'-'
·-
�
20
(l)
(l) 0.088%
10
0.041%
0
1 10 100 1000
Number of cycles
Figure 4.29 Effect of number of stress cycles on hysteretic damping for

medium dense sand (after Silver and Seed, 1971)
Source: Silver, M.L., and Seed, H.B. (1971). "Deformation Characteristics of Sands Under Cyclic
Shear stress, T
Shear strain, r'
/ 1
--2r;--
Figure 4.30 Bilinear idealization of shear-stress-versus-shear-strain plots
Other parameters remaining the same (that is, Rn , number of cycles, and ampli
tude of shear strain), a vertical stress increase will decrease the damping ratio. In
many seismic analysis studies, it is convenient to represent the nonlinear shear-stress
versus-shear-strain relationship in the form of a bilinear model (also see Figure 7 .10),
as shown in Figure 4.30 (Thiers and Seed, 1968). In this figure G1 is shear modulus
up to a limiting strain of y', and G2 is the modulus for strain beyond y',.
Advantages of the Cyclic Simple Shear Test

There are several advantages in conducting cyclic simple shear tests. They are
more representative of the field conditions, since the specimens can be consoli
dated in K0 state. Solid soil specimens used in resonant column tests can provide
good results up to a shear strain amplitude of about 10- 3 %. Similarly, the hollow
samples used in resonant column studies provide results within a strain amplitude
range of 10-3 % to about 1%. However, cyclic simple shear tests can be conducted
for a wider range of strain amplitude ( that is, 10-2 % to about 5%. This range is the
general range of strain encountered in the ground motion during seismic activities.
The pore water pressure developed during the vibration of saturated soil spec
imens by a resonant column device is not usually measured. However, in cyclic
simple shear tests, the pore water pressure can be measured at the boundary (see
Section 10.10 and Figure 10.20).
BJI CYCLIC TORSIONAL SIMPLE SHEAR TEST

Another technique used to study the behavior of soils subjected to cyclic load
ing involves a torsional simple shear device. The torsional simple shear device
accommodates a "doughnut-like" specimen, as shown in Figure 4.31 (Ishibashi
and Sherif, 1974). The specimen has outside and inside radii of r 1 = 50.8 mm and
r2 = 25.4 mm. The outside and inside heights of the specimen are h 1 = 25.4 mm
and h 2 = 12.7 mm. The soil is initially subjected to a vertical effective stress av ,
an outside and inside horizontal effective stress of iih, and a cyclic shear stress of
r (Figure 4.32). When a shear stress r is applied, line AB moves to the position
of A' B' (Figure 4.32). So, the shearing strain is
, r10 , r20
YA =
-and rB =
-
h1 h2
For uniform shear strain throughout the sample,
r, -r' A - B
or
A-- -�-
, '
--A
/'
Doughnut-like specimen
�,
(a)
/-<
/
,, '
,
/ '
, '
'
,
/ '
Section at A-A
(b)
Figure 4.31 Soil specimen for torsional simple shear test

4.8 Cyclic Torsional Simple Shear Test
,____ ' 1 -------a•-I
Figure 4.32 Applied stresses on a torsional simple shear test specimen
r10
-----
r20
h, h2
So
� (4.35)
�
The following can be calculated after application of the horizontal shear

stress on the specimen.
Major effective principal stress:
(4.36a)
Intermediate effective principle stress:

(4.36b)
Minor principal effective stress:
(4.36c)
With proper design [Eq. (4.35)], a cyclic torsional shear device can apply near
uniform shear strain on the specimen. It can apply shear strains up to about 1 %.
It also eliminates any sidewall frictional stresses that are encountered in cyclic
simple shear tests.
The shear modulus of a specimen tested can be determined as
amplitude of shear stress, r
G=
amplitude of shear strain, r'
The damping ratio corresponding to a given shear strain amplitude can be deter
mined by using Figure 4.21 and Eq. (4.34).
Liquefaction studies on saturated granular soils can also be conducted by
this device along with pore water pressure measurement.
BJ CYCLIC TRIAXIAL TEST

Cyclic triaxial tests can be performed to determine the modulus of elasticity E
and the damping ratio D of soils. In these tests, in most cases, the soil specimen
is subjected to a confining pressure a O = a 3• After that, an axial cyclic stress �ad
is applied to the specimen, as shown in Fi gure 4.33. The tests conducted for the
evaluation of the modulus of elasticity and damping ratio are strain-controlled
tests. A servo-system is used to apply cycles of controlled deformation.
Figure 4.34 shows the nature of a hysteresis loop obtained from a dynamic
triaxial test. From this,
�a d
E=
e
Once the magnitude of Eis determined, the value of shear modulus can be
calculated by assuming a representative value of Poisson's ration, or
E
G=
2(1 + µ)
Again referring to Figure 4.34, the damping ratio can be calculated as
area of the hysteresis loop
D = _I_
21t area of triangle OAB and OA' B'
4.9 Cyclic Triaxial Test
CJ3 = CJo
•
� Soil specimen -
Figure 4.33 Cyclic triaxial test
Axial strain, £
A'
Figure 4.34 Determination of damping ratio from cyclic triaxial test
Stress-controlled dynamic triaxial tests are used for liquefaction studies on

saturated granular soils (see Chapter 10).
A more elaborate type of dynamic test device has also been used by several
investigators to study the cyclic stress-strain history and shear characteristics of
soils. Matsui, Ohara, and Ito (1980) used a dynamic triaxial system that could
generate sinusoidally varying axial and radial stresses.
A. Cyclic Strength of Clay

During earthquakes, the soil underlying building foundations and in structures
such as earth embankments, is subjected to a series of vibratory stress applica
tions. These vibratory stresses may induce large deformation in soil and thus
failure. In order to evaluate the strength of clay under earthquake loading con
ditions, Seed and Chan (1966) conducted a number of dynamic triaxial tests.
Lia-= amplitude of
t pulsating stress
a-1 = sustained
axial stress
Soil
0"3 � specnnen
Time
(a) (b)
T T
00
-.....
00
(!) 00
I-<
Lia-
-.....
+--' (!)
00 I-<
1
+--'
ro 00
>< ro
< <
>< Lia-
l l
Time Time
(c) (d)
Figure 4.35 Stress conditions on a soil specimen

[Note: (b) One-directional loading with symmetrical stress pulses; (c) and (d) one-directional loading with
nonsymmetrical stress pulses]
4.9 Cyclic Triaxial Test
Deviator stress (kPa)

0 100 200 300 400 500
0r-----e---.;;;;;;t;:=-----r-----�------r---�-----i
-------- r---- -l Strain induced during

normal strength test at
factor of safety= 1.5
8 3 Strain induced
4 during
c----==-=---� \
-:)l0 earthquake
c--11111111:=---�-=,5. 0 loading
_J __
I
Stress-strain curve �
for sample in normal \ Stress-strain curve
strength test \
24 for sample after
I
\ earthquake loading
Vicksburg silty clay, sample No. 8 I
Water content= 22.1%

Dry density= 16.24 kN/m3
32 Degree of saturation= 93%
Confining pressure, a-3= 100 kPa
36�----�-----�----�-----�----�
Figure 4.36 Stress-versus-strain relationship for Vicksburg silty clay under

sustained and axial pulsating stress (after Seed and Chan, 1966)
Source: Seed, H.B., and Chan, C.K. (1966). "Clay Strength Under Earthquake Loading
Conditions," Journal of the Soil Mechanics and Foundation Division, ASCE, Vol. 92, No. SM2,
Figure 4.35 shows the nature of some of the stress conditions imposed on the
soil specimens during those tests. The results of this study are very instructive
and are described in some detail in this section.
Figure 4.36 shows the results of a laboratory test on a specimen of Vicks
burg silty clay subjected to sustained and pulsating stresses. The specimen with
a degree of saturation of 93% was initially subjected to a confining pressure of
a 3 = 100 kPa and then to a conventional axial loading in undrained conditions
up to 66% of its static strength. This implies that the sustained stress a 1 - a 3 was
equal to 0.66 [a1 cn -a3 ], which corresponds to a factor of safety of 1.5. At this
time the axial deformation of the specimen was about 5%. After that, 100 tran
sient stress pulses were applied to the specimen. (Note: Loading type is similar
to that shown in Figure 4.35b). These stress pulses induced an additional axial
strain of about 11%, although the static strength was never exceeded.
Figure 4.37 shows the nature of soil deformation on three soil specimens of
San Francisco Bay mud subjected to pulsating stress levels to 100%, 80%, and
60% of normal strength (that is, static strength). For these tests, no sustained stress
was applied. (Note: Loading type is similar to that shown in Figure 4.35d.) It is
worth noting that, for each level of pulsating stress, the specimen ultimately failed.
Number of stress pulses

10 100 1000 10,000
Moisture content= 91%
No sustained load
Pulsating stress Pulsating stress Pulsating stress

equal to 100% of equal to 80% of equal to 60% of
normal strength normal strength normal strength
24 ,.__________.________.______........._________.
Figure 4.37 Deformation of San Francisco Bay mud specimens subjected to

pulsating stress (after Seed and Chan, 1966)
Figure 4.38 is a plot of the pulsating stress level (as a percent of normal
strength) versus sustained stress level (as a percent of normal strength) causing
failure of San Francisco Bay mud at various numbers of transient stress pulses.
As the number of stress pulses are increasing at the same pulsating stress level,
the sustained stress level inducing failure is decreasing. The interested readers
should refer to the original paper by Seed and Chan (1966). Similar plots could be
developed for various soils to help in the design procedure of various structures.
mJ SUMMARY OF CYCLIC TESTS

In the preceding sections, various types of laboratory test methods were pre
sented, from which the fundamental soil properties such as the shear modulus,
modulus of elasticity, and damping ratio are determined. These parameters are
used in the design and evaluation of the behavior of earthen, earth-supported,
and earth-retaining structures. As was discussed in the preceding sections, the
magnitudes of G and D are functions of the shear strain amplitude y'. Hence,
while selecting the values of G and D for a certain desi gn work, it is essential to
know the following:
a. Type of test from which the parameters can be obtained
b. Magnitude of the shear-strain amplitude at which these parameters needs
to be measured
4.10 Summary of Cyclic Tests
100
80
= Number of
20
0 ,.._____._____.___....______.___
0 20 40 60 80 100
Sustained stress
Normal strength (%)
Figure 4.38 Combinations of sustained and pulsating stress intensities causing

failure-San Francisco Bay mud (after Seed and Chan, 1966)
For example, strong ground motion and nuclear explosion can develop large
strain amplitudes whereas some sensitive equipment such as electron micro
scopes may be very sensitive to small strain amplitudes.
Figure 4.39 provides is a useful reference table for geotechnical engineers, as
it gives the amplitude of shear-strain levels, type of applicable dynamic tests, and
the area of applicability of these test results. Despite the fact that laboratory test
ing is not ideal, it will continue to be important because soil conditions can be bet
ter controlled in the laboratory. Parametric studies necessary for understanding
the soil behaviour of soils under dynamic loading conditions must be performed
in the laboratory conditions. Table 4.2 provides a comparison of the relative qual
ities (what property can be measured and what is the degree of quality of the
measured property) of various laboratory techniques for measuring dynamic
soil properties. Similarly, Table 4.3 gives a summary of the different engineering
parameters that can be measured in different dynamic or cyclic laboratory tests.
Dynamic tests in -
Dynamic
triaxial and shear
field tests apparatus
Resonance
tests
I-+- Vibration---+!
table test
Static plate---+l
load tests
--Earthquakes-
I
lQ-7 lQ-6 10-s lQ-4 10-3 10-2 10-1
Shear deformation, y(%)
Figure 4.39 Range and applicability of dynamic laboratory tests
Table 4.2 Relative Quality of Laboratory Techniques for Measuring Dynamic

Soil Propertiesa
Relative Quality of Test Results
Effect of
Shear Young's Material number of
modulus modulus damping cycles Attenuation
Resonant column Good Good Good Good
with application Fair
Ultrasonic pulse Fair Fair Poor
Cyclic triaxial Good Good Good
Cyclic simple shear Good Good Good
Cyclic torsional shear Good Good Good
aAfter Silver (1981)
Source: Silver, M.L. (1981). "Load deformation and Strength Behavior of Soild under Dynamic
Loading," State-of-the-Art Paper, Proceedings, International Conference on Recent Advances in
geotechnical Earthquake Engineering and Soil Dynamics (Ed. Shamsher Prakash), Vol. 3, pp. 873-896.
Field Test Measurements
11D REFLECTION AND REFRACTION OF ELASTIC BODY

WAVES-FUNDAMENTAL CONCEPTS
When an elastic stress wave impinges on the boundary of two layers, the wave
is reflected and refracted. As has already been discussed in Chapter 3, there are
two types of body waves-that is, compression waves (or P-waves) and shear
waves (or S-waves). In the case of P-waves, the direction of the movements of
4.11 Reflection and Refraction of Elastic Body Waves-Fundamental Concepts
Table 4.3 Parameters Measured in Dynamic or Cyclic Laboratory Testsa

Resonant column Cyclic triaxial Cyclic simple shear Torsional shear
Load Resonant frequency Axial force Horizontal force Torque
Deformation
Axial Vertical Vertical Vertical Vertical
displacement displacement displacement displacement
Shear Acceleration Not measured Horizontal Rotation
displacement
Lateral Not usually Not usually Not usually Often
measured measured measured controlled
Volumetric None for undrained tests
Volume of fluid moving into or out of the sample for drained tests
Pore water Not usually Measured at Measured at Measured at
pressure measured boundary boundary boundary
aAfter Silver (1981)
Source: Silver, M.L. (1981). "Load deformation and Strength Behavior of Soild under Dynamic Loading,"
State-of-the-Art Paper, Proceedings, International Conference on Recent Advances in geotechnical Earthquake
Engineering and Soil Dynamics (Ed. Shamsher Prakash), Vol. 3, pp. 873-896.
the particles coincides with the direction of propagation. This is shown by the
arrows in Figure 4.40a. The shear waves can be separated into two components:
a. SV-waves, in which the motion of the particles is in the plane of propaga
tion as shown by the arrows in Figure 4.40b
b. SH-waves, in which the motion of the particles is perpendicular to the
plane of propagation, as shown by a dark dot in Figure 4.40c
a C
Ray of P-wave
Ray of SV-wave
b d
(a) e (b)
Ray of SH-wave
(c)
Figure 4.40 (a) P-wave; (b) SV-wave; and (c) SH-wave

If a P-wave impinges on the boundary between two layers, as shown in

Figure 4.41a, there will be two reflected waves and two refracted waves. The
reflected waves consist of ( 1) a P-wave shown as P1 in layer 1 and (2) an SV-wave
shown as SVi in layer 1.
The refracted waves will consist of ( 1) a P-wave, shown as P2 in layer 2, and
(2) an SV-wave, shown as SV2 in layer 2. Referring to the angles in Figure 4.41a,
it can be shown that
(4.37)
and
(4.38)
SH
Layer 1 Layer 1
VP VP
1 1
Vs Vs
l l
Interface Interface
Layer 2 Layer 2
VP VP
2 2
Vs Vs
z z
SH2
(a) (b)
Interface
(c)
Figure 4.41 Reflection and refraction for (a) an incident P-ray; (b) an incident
SH-ray; and (c) an incident SV-ray
4.12 Seismic Refraction Survey (Horizontal Layering)
where
vP1 and v J>i = the velocities of the P-wave front in layers 1 and 2, respectively
vs1 and Vs2 = the velocities of the S-wave front in layers 1 and 2, respectively
If an SH -wave impinges the boundary between two layers, as shown in Fig

ure 4.41b, there will be one reflected SH-wave (shown as SH1) and one refracted
SH-wave (shown as SH2 ). For this case
/31 = /32 (4.39)
and
sin /31 _ sin /32 _ sin {33

(4.40)
Vsi Vsi Vs2
Finally, if an SV-wave impinges the boundary between two layers, as shown

in Figure 4.41c, there will be two reflected waves and two refracted waves. The
reflected waves are (1) a P-wave, shown as Pi in layer 1 and (2) an SV-wave, shown
as SVi. in layer 1. The refracted waves are (a) a P-wave, shown as Pi in layer 2, and
(b) an SV-wave, shown as SV2 in layer 2. For this case, /31 = /32 :
(4.41)
The mathematical derivations of these facts will not be shown here. For
further details the reader is referred to Kolsky (1963, pp. 24--38).
IJD SEISMIC REFRACTION SURVEY (HORIZONTAL

LAYERING)
Seismic refraction surveys are sometimes used to determine the wave propaga
tion velocities through various soil layers in the field and to obtain thicknesses
of each layer. Consider the case where there are two layers of soil, as shown in
Figure 4.42a. Let the velocities of P-waves in layers 1 and 2 be vP1 and v J>i , respec
tively, and let vPi < vPl. A is a source of impulsive energy. If seismic waves are
generated atA, the energy from the point will travel in hemispherical wave fronts.
Consider the case of P-waves, since they are the fastest. If a detecting device is
placed at point B, which is located at a small distance x from A, the P-wave that
travels through the upper medium will reach it first before any other wave. The
travel time for this first arrival may be given as
(4.42)
whereAB=x.
A B G
.... � . . -- .. ·.. . .. .. � . .. .. .. ..· .,.. . . . . : .
.....
....
.....
t'd
·i:: b
....
t'd
.....
Slope of Oa = v1
P1
1
Slope of ab = -
VP2
Distance, x
Figure 4.42 Seismic refraction survey-horizontal layering
Again, consider the first arrival time of a P-wave at a point G, which is

located at a greater distance from A. In order to understand this, one considers a
spherical P-wave front that originates at A striking the interface of the two layers.
At some point C, the refracted P-wave front in the lower medium will be such
that the tangent to the sphere will be perpendicular to the interface. In that case,
the refracted P-ray (shown as Pi in Figure 4.42a) will be parallel to the boundary
and will travel with a velocity v 1>2 . Note that because vPI < v 1>2 , this wave front will
travel faster than those described previously. From Eq. (4.38)
Since a3 = 90 °, sin a3 = 1, and
(4.43)
where ac = critical angle of incidence.

The wave front just described traveling with a velocity vP2 will create vibrat
ing stresses at the interface, and this will generate wave fronts that will spread
out into the upper medium. These P-waves will spread with a velocity of v Pi·
The spherical wave front traveling downward from D in layer 2 will have a radius
equal to DE after a time 11t. At the same time 11t, the spherical wave front trav
eling upward from point D will have a radius equal to DF. The resultant wave
front in the upper layer will follow a line EF. It can be seen from the diagram that
v PI 11t DF
= - = Slnlc (4.44)
v /J2 11t DE
So ray DFG will make an angle ic with the vertical. It can be mathematically
shown that for x greater than a critical value Xe, the P-wave that travels the path
ACDG will be the first to arrive at point G. Let the time of travel for the P-wave
along the path ACDG be equal tot. Thus, t = tA c + tcD + tD c, or
x tan i,
t=
( co: i, )( v� ) + - �� + ( co: i, )( v� )
x 2zsinic 2z
=-----+---
VP2 VP2 COS ic VPI COS ic
where x =AG .But v P2 = v PI /sin ic [from Eq. (4.44)]; thus
2-- 2zsin2 ic 2z 2- 2z 1- sin2 ic
t= + = + ( )
VP2 Vp1 COSlc V p1 COSlc V P2 Vp1 COSlc
(4.45)
X 2z
=-+-coszc
VP2 VPI
Since sin ic = vPI IvP2
(4.46)
Substituting Eq. (4.46) into Eq. (4.45), one obtains
(4.47)
If detecting instruments are placed at various distances from the source of

disturbance to obtain first arrival times and the results are plotted in graphical
form, the graph will be like that shown in Figure 4.42b. The line Oa represents
the data that follow Eq.(4.42). The slope of this line will give l/vp1 . The line ab
represents the data that follow Eq. (4.47). The slope of this line is 1/v/>}_ . Thus the
velocities of vPI and v/>}_ can now be obtained.
If line ab is projected back to x = 0, one obtains
2
2z \JfvPl. -vPl2
.=
=f
t
(vp1)(vJ>2 )
l
or
z = (ti)( )( V Pl. ) = ti VPi

Vp
1
(4.48)
2�v� - v� 2 COS ie
where is the intercept time. Hence, the thickness of layer 1 can be eas
ti
ily obtained.
The critical distance Xe (Figure 4.42b) beyond which the wave refracted at the
interface arrives at the detector before the direct wave can be obtained by equat
ing the right-hand sides of Eqs. (4.42) and (4.47):
� �
Xe =-+-----
Xe 2z�v -v
-
vP 1 V Pl. VP1 VPl.
or
Vp1VJ>l. VJ>l. +vPI

= 2z (4.49)
VJ>l. - Vp1 VJ>l. - VPI
The depth of the first layer can be calculated from Eq. (4.49) as
Xe
z=- (4.50)
2
A. Refraction Survey in a Three-Layered Soil Medium

Figure 4.43 considers the case of a refraction survey through a three-layered soil
medium. Let vPI, v/>}_ and vPJ be the P-wave velocities in layers 1, 2, and 3, respec
tively, as shown in Figure 4.43a (vPi < vPl. < vPJ ). If A in Figure 4.43a is a source
of disturbance, the P-wave traveling through layer 1 will arrive first at B, which
is located a small distance away from A. The travel time for this can be given by
Eq. (4.42) as t = x/vP1 • At a greater distance x, the first arrival will correspond
to the wave taking the path A CDE. The travel time for this case be given by
Eq. (4.47) as
A B E K
( V]
Layer 1
P1
VP !
. = sm
ic1 · -1
vP
2
C
vP
ic2 = sin-1 ( 2]
VP3
H I
(a)
.....
..........
Q)�
C
·i::
I-<
.....
.....
I-<
1i 2 Slope of Oa = -1-
vP
,
1
Slope of ab = -
VP2
Slope of be = -1
VP
3
0 Distance, x
(b)
Figure 4.43 Refraction survey in a three-layer soil
x + 2z1 v� - v�1 �
t=- -----'-----
vPl (vPI) (vPl)
where z1 = thickness of top layer.
At a still larger distance, the first arrival will correspond to the path AGHIJK.
Note that the refracted ray H-1 will travel with a velocity of v/J3 . The angle ic2 is
the critical angle for layer 3.
(4.51)
For this path (AGHIJK) the total travel time can be derived as
X __1 �
2z v 2 � /J3___P2_v2 2z v2 - v2
t = -+ __/J3_- 1 +--2 __
__p_
(4.52)
v /J3 (v /J3 )(v p1 ) ( v /J3 )(v P2 )
where z2 = thickness of layer 2.

So, if detecting instruments are placed at various distances from the source
of disturbance to obtain first arrival times, they can be plotted in a t-versus-x
graph. This graph will appear as shown in Figure 4.43b. The line Oa corresponds
to Eq. (4.42), ab corresponds to Eq. (4.47), and be corresponds to Eq. (4.52). The
slopes of Oa, ab, and be will be ljv P1 , l/v 1>2 and 1/v /J3 , respectively. The thickness
of the first layer z1 can be determined from the intercept time til in a similar man
ner, as shown in Eq. (4.48), or
The thickness of the second layer can be obtained from Eq. (4.52). Referring
to Figure 4.43b, the expression for the intercept time ti2 can be evaluated by sub
stituting x = 0 into Eq. (4. 52):
2z1 �v� - v�1 2z2 �V�3 - v�

t = ti2 = ----- + -----
( VP3 ) (V pJ ( VP3 )( VPl )
or
(4.53)
B. Refraction Survey for Multilayer Soil

In general, if there are n layers, the first arrival time at various distances from
the source of disturbance will plot as shown in Figure 4.44. There will be n seg
ments in the t-versus-x plot. The slope of the nth segment will give the value
ljv Pn (n = 1, 2, ...). More details on advanced test methods (detecting inclination
of the bedrock is briefly described later) and interpretation could be found in
several geophysics books.
The value of P-wave velocity in a natural deposit of soil will depend on sev
eral factors, such as confining pressure, moisture content, and void ratio. Some
typical values of vP are given in Table 3.1. It is worth noting that P-wave veloc
ity through saturated soils will be approximately 1500 m/sec. However, P-wave
velocity value could reach the order of few kilometers per second in case of
rocks.
.....
ro
Layer 1, vP .....I-<
1
ro
Layer 2, v
Pi
I-<
•
• Segment n
• Slope= -
1
VPn
Layer n, VPn
Distance, x
Figure 4.44 Refraction survey for multilayer soil
EXAMPLE4.1
Following are the results of a refraction survey (horizontal layering of soil).
Determine the P-wave velocities of the soil layers and their thicknesses.
Distance (m) Time of first arrival (ms)

2.5 5.5
5.0 11.1
7.5 16.1
15.0 24.0
25.0 30.8
35.0 38.2
45.0 46.1
55.0 51.3
60.0 52.8
SOLUTION
The time-distance plot is given in Figure 4.45. From the plot,
5
----=472m/s
10.6 X 10-3
10
= 1389 m/s
7.2 X 10-3
lO
= 3333 m/s
3x10-3
70----�-----------------------
60
--- - ----- - - - 10
---
- - - - - - ---
-- - ---
10
: 10.6
0 _________._____......_____._________.______._________.
0 10 20 30 40 50 60
Distance, x (m)
Figure 4.45
tn = 13.3 X 10- 3 s; ti2 = 35.6 X 10- 3 s.
From Eq. (4.48)
z 1 _ (ti 1 )(vpJ(vP1 ) _
-
(13. 3 X l0-3 )(472)(1389)
- 2Jv� -v� 2�(1389)2 - ( 472)2
= 3.39 m
From Eq. (4.53)
_1 2z1Jv�3 -v� (vp3) (vP1 )

Z2 - - [ l-1 2 - ---'----- ]
2 (vP3 )(v/Jl ) Jv� -v�
_ (2)(3.39)�(3333)2 -(472)2 ]
= _!_[35_6 X 10_ 3
2 (3333)(472)
(3333)(1389)
X
)(3333)2 -(1389)2
= _!_(0.02138) (1528) = 16.33 m

2
4.13 Refraction Survey in Soils with Inclined Layering
mJ REFRACTION SURVEY IN SOILS

WITH INCLINED LAYERING
Figure 4.46a shows two soil layers. The interface of soil layers 1 and 2 is inclined
at an angle f3 with respect to the horizontal. Let the P-wave velocities in layers 1
and 2 be v PI and v Pl' respectively ( v P1 < v Pl).
If a disturbance is created at A and a detector is placed at B, which is small
distance away from A, the detector will first receive the P-wave traveling through
layer 1. The time for its arrival may be given by
------------- X _______________,
1---------x-----------
A B E
·.,··
:: : -�- .. .- . . ' �...: ·. ·.·:.·
AF=z' ,' ,'A 1 - - - - - - - I
- I
I
z'
I
I
I
I
I
I
I
- - �I
I
A4
I
I
,' EG=z"
I
I
(a)
lu versus x
b'
� a' ....
-s
-- - ---
.....
-
..........,
Q)�
1iu
ro ro
> >
·i:: ·i::
I-< I-<
ro
....., .....,
.....� l-
I-<
�
I-<
id
td versus x
Distance
x increasing for td x increasing for lu
(b)
Figure 4.46 Refraction survey in soils with inclined layering
However, at a larger distance the first arrival will be for the P-wave following the
path ACD£-which consists of three parts. The time taken can be written as
(4.54)
Referring to Figure 4.46a,
= z'
tAC (4.55)
VPI COS le
AAt - AA1 - A2A3 - A3/4i

CD = ----------
ten= -
vP2 '\) P2
(4.56)
_ xcos/3- z'tanie - z'tanie - xsinf3tanie
Vp2
z' xsin/3
+
D A3 + A3E - ------
tnE = ---- cosie cosie
(4.57)
VP1 VP1
Substitution of Eqs. (4.55), (4.46), and (4.57) into Eq. (4.54) and simplification
yields
(4.58)
Now, if the source of disturbance is E and the detector is placed at A, the first
arrival time along the refracted ray path may be given by
2z"cos� x . ' /3)

tu =---+-sm ( le - (4.59)
VPI VPI
In the actual survey, one can have a source of disturbance such as A and
observe the first arrival time at several points to the right of A and have a source
of disturbance such as E and observe the first arrival time at several points to the
left of E. These results can be plotted in a graphical form, as shown in Figure
4.46b (time-versus-x plot). From Figure 4.46b note that the slopes of Oa and O'a'
are both l/v P1 • The slope of the branch ab will be [sin(ie + f3)]/v P1 ' as can be seen
from Eq. (4.58). Similarly, the slope of the branch of a'b' will be [sin(ie - /3)]/v /JJ
[see Eq. (4.59)]. Let
sin((+ /3)
md =
---- (4.60)
vP1
and
(4.61)
From Eq. (4.60),
(4.62)
Again, from Eq. (4.61)
ic = sin-1(m uv J}J ) + /3 (4.63)
Solving the two preceding equations,
I i = ½[sin-1(v md) + sin-1(v ,m )] I

c P1 p u (4.64)
and
/3 -- 21 [Sill
· -1 ·
(Vp md ) - Sill-1(Vp mu )]
1 1 (4.6 5)
Once ic is determined, the value of v P'2 can be obtained as
VP2 = VP
•
t
• (4.66)
Sill le
Again referring to Figure 4.46b, if the ab and a' b' branches are projected
back, they will intercept the time axes at tid and tiu, respectively. From Eqs. (4.58)
and (4.59), it can be seen that
2z' cos ic
tid =
----
VPl
or
z' = (tid )vp1

(4.67)
2 cos(
and
=
2z" cosic
fiu ----
Vpi
or
z" = ( fiu ) VPt (4.68)

2 COS ic
Since ic and vPi are known and tid and tiu can be determined from a graph, one can
obtain the values of z' and z".
EXAMPLE 4.2
Ref erring to Figure 4.46a, the results of a refraction survey are as follows. The
distance between A and E is 60 m.
Point of disturbance A Point of disturbance E

Distance Time of first Distance from Time of first
from A (m) arrival (ms) E(m) arrival (ms)
5 12.1 5 11.5
10 25.2 10 22.8
15 35.3 15 34.5
20 48.0 20 44.8
30 60.2 30 69.1
40 68.5 40 78.1
50 76.8 50 82.8
60 85.1 60 87.7
Determine
a. VPI and Vpz ,
b. z' and z", and
C. {3
SOLUTION
The time-distance records have been plotted in Figure 4.47.
a. From branch Oa,
10
vPl =---=400 m/s
25 X 10- 3
From branch O' a'
10
vPl =---= 454 m/s
22 X 10-3
The average value of v PI is 427 mis.
From the slope of branch ab,
md = 8.8 X l0-3 =0.88 X 10-
3
10
Again, from the slope of branch a'b',
m = sx10-3 =0.5 10-
u
10 x 3
50 40 30 20 10
o--------------------------1,-------------
10 20 30 40 50 60
Distance, x (m)
Figure 4.47
From Eq. (4.64),

ic = ½ [Sin - l (V PI md ) + Sin- l ( VPI mu )]
sin- 1( vPi md ) = sin- 1[( 427)(0.88 X 10- 3)] = 22.07°
sin- 1(vPi mu ) = sin- 1 [(427)(0.5 X 10- 3)] = 12.33 °
Hence
ic = ½(22.07° + 12.33 °) = 17.2°
Using Eq. (4.66)
VP -4-2 -7-
V P2 = • I• = 1444 m/s
Sln le sin (17.2)
b. From Eq. (4.67)
z' = (tid )(vpJ

2 cosic
tid = 35.9 X 10- 3 s (from Figure 4.47). So
3
z' = (35.9 X 10- )( 427) = S.03 m
2 cos(l7.2)
Again, from Eq. (4.68),

( liu ) ( VpJ
z ,, =---
2 COS ic
From Figure 4.47, tiu = 59.8 X 10- 3 s.
"(59.8 X 10- 3)( 427)
z = = 13_37 m
2cos (17.2)
c. From Eq. (4.65),
/3 =½ [sin- 1 (vp1 md )-sin- 1 (vPi mu )]=½ (22.07° -12.33° ) = 4.87°
DD REFLECTION SURVEY IN SOIL

(HORIZONTAL LAYERING)
Reflection surveys can also be conducted to obtain information about the soil
layers. Figure 4.48a shows a two-layered soil system. A is the point of distur
bance. If a recorder is placed at Cat a distance x away from A, the travel time for
the reflected P-wave can be given as
_AB+BC _ 2
t------- z
�(
+ -
x2J
(4.69)
VP1 VP1 2
where t = total travel time for the ray path ABC.

From Eq. (4.69), the thickness of layer 1 can be obtained as
(4.70)
If the travel times t for the ref lected P-waves at various distances x are
obtained, they can be plotted in a graphical form, as shown in Figure 4.48b.
Note that the time-distance curve obtained form Eq. (4.69) will be a hyperbola.
The line Oa shown in Figure 4.48b is the time-distance plot for the direct P-waves
traveling through layer 1 (compare line Oa in Figure 4.48b to the line Oa in
Figure 4.42b). The slope of this line will give l/vp1 •
If the time-distance curve obtained from the reflection data is extended back,
it will intersect the time axis at t0 • From Eq. (4.69) it can be seen that at x = 0,
2z
to =-
vP1
4.14 Reflection Survey in Soil (Horizontal Layering)
x----------------...i
C
z
I
I
�� Layer 1
B Layer 2
(a)
a
.....
[Eq. (4.69)]
0
Distance, x
(b)
x2
(c)
Figure 4.48 Reflection survey in soil-horizontal layering
or
(4.71)
With vPi and t0 known, the thickness of the top layer z can be calculated.
Another convenient way to interpret the reflection survey record is to plot a
graph of t 2 versus x2• From Eq. (4.69),
2
t =+[z 2
2 +(x) ]=+(4z 2 +x2 ) (4.72)
V
PI 2 V
Pl
This relation indicates that the plot of t 2 versus x 2 will be a straight line, as
shown in Figure 4.48c . The slope of this line give l/v�1 and the intercept on the t 2
axis will be equal to tJ. Substituting t = t 0 and x = 0 into Eq. (4.72),
4z 2
t0 -_
2
- 2
V
P1
or
t 02Vp12
z2
= (4.73)
With tJ and v�1 known, the thickness of the top layer can now be calculated.
EXAMPLE 4.3
The results of a reflection survey on a relatively flat area (shale underlain by
granite) are given here. Determine the velocity of P-waves in the shale.
Distance from point Time for first

of disturbance (m) reflection (s)
30 1.000
90 1.002
150 1.003
210 1.007
270 1.011
330 1.017
390 1.023
4.14 Reflection Survey in Soil (Horizontal Layering)
SOLUTION
Using the time-distance records, the following table can be prepared.
X x2 t ,2
(m) (m 2) (s) (.\' 2)
30 900 1.000 1.000
90 8,100 1.002 1.004
150 22,500 1.003 1.006
210 44,100 1.007 1.014
270 72,900 1.011 1.022
330 108,900 1.017 1.034
390 152,100 1.023 1.046
A plot of t 2 versus x 2 is shown in Figure 4.49. From the plot,
(�)2 79200
= 1816.6 m/s
(�t)2 0.023
1.07 ------------------------�
1.05
(�t)2 = 0.023
� 1.03
�
(Ax)2 = 79,200
1.01
0.99
0 37,000 74,000 111,000 148,000 185,000
Figure 4.49
11D REFLECTION SURVEY IN SOIL (INCLINED LAYERING)

Figure 4.50 considers the case of a reflection survey where the reflecting bound
ary is inclined at an angle f3 with respect to the horizontal. A is the point for the
source of disturbance. The reflected P-ray reaching point C will take the path
ABC. Referring to Figure 4.50,
AB+BC= A'B +BC= A'C
But
(A'C) 2 = (A'A2 ) 2 + (A2C) 2 (4.74)
A' A2 = AA' cos f3 = 2z' cos f3 (4.75)
A2C=A2A+ AC=2z' sin/3 + Xe (4.76)
Substituting Eqs. (4.75) and (4.76) into Eq. (4.74),
A'C = )(2z' cos/3)2 + (2z' sin/3 + Xe )2

= )4z'2 + x� + 4z'xe sin/3
Thus, the travel time for the reflected P-wave along the path ABC will be
=
A'C
te --
vPI
----------XE----------
-------Xe--------
Layer 2
VP2
Figure 4.50 Reflection survey in soil-inclined layering
4.15 Reflection Survey in Soil (Inclined Layering)
So
1
te = - -�4z' 2 + x� + 4z'xesin/3 (4.77)
VPI
In a similar manner, the time of arrival for the reflected P-waves received at
point E can be given as
(4.78)
Combining Eqs. (4.77) and (4.78),
v� (t}__
- fl;) XE + Xe
Sln. /3 = _r_, _
4z'(xE - Xe) 4z'
(4.79)
Now, let t = (tE + te)/2 and lit = tE - te
Substitution of the preceding relations in Eq. (4.79) gives
vp 1t(lit)
2
XE+ Xe
. /3 = -�---
sin (4.80)
2z'(xE - Xe) 4z'
If Xe is equal to zero, Eq. (4.80) will transform to
v! t(lit) _ xE
sin /3 = (4.81)
2z'xE 4z'
If Xe = 0 and /3 = 0 (that is, the reflecting layer is horizontal) then, from Eq.
(4.81),
2
A t =--
XE
o. (4.82)
2 2vPi t
If Xe = 0 and lit > x}l2v! t, the reflecting layer is sloping down in the direc
tion of positive x, as shown in Figure 4.50. If Xe= 0 and lit< x}l2v�1 t, the
reflecting layer is sloping down in the direction of negative x (that is, opposite to
that shown in Figure 4.50).
In actual practice, that point of disturbance A (Figure 4.51) is generally
placed midway between the two detectors, so xE = -xe = x. So, from Eq. (4.80)
Sln
. /3 = v� t(lit)
---=---
ri __
(4.83)
4z'x
----x----------------x----,
. ·:.·.··. . .... . · . . . . . . . . . . . · .. ·_ . ·
. . .
. .- . . . . . . . . ·. · . · . . . ·.-:-·. . . ·. .
Layer 1
VP1
I
/
I /
I
I
I
/
/
I
I
A'
Figure 4.51
Referring to Figure 4.51, AA'= 2z' = (A'C + A'E)/2. So
2z' 1 A'C A'E 1

-=- (-+-J = -(le+l
E)
= l (4.84)
V/JI 2 Vp1 Vp1 2
Combining Eqs. (4.83) and (4.84),
. V (�l)
/3 =�--
Sln p1
(4.85)
2x
EXAMPLE 4.4
Refer to Figure 4.51. Given: x = 85.5 m, le = 0.026 s, and l E = 0.038 s. Deter
mine f3 and z'. The value of v PI i.e., the velocity of the primary wave through
the top layer has been previously determined to be 410 mis.
SOLUTION
0.026 + 0.038 = _
= le + l E =
l 0 032 S
2 2
�l =l E - le = 0.038 - 0.026 = 0.012 S
4.16 Subsoil Exploration by Steady-State Vibration
From Eq. (4.85)
/3 = sin- 1 [vP1 (�t)] = sin- 1 [ ( 410 )(0.012)] = 1.650

2x (2)(85.5)
From Eq. (4.84)
or
' t v/JJ = (0.032)( 410 ) =

z = 6_56 m
2 2
11[!1 SUBSOIL EXPLORATION BY

STEADY-STATE VIBRATION
In steady-state vibration, a circular plate placed on the ground surface is vibrated
vertically by a sinusoidal loading (Figure 4.52a). This vibration will send out
Rayleigh waves (Section 3.1 3), and the vertical motion of the ground surface will
predominantly be due to these waves. This can be picked up by motion transduc
ers. The velocity of the Rayleigh waves can be given as
v, = jL (4.86)
where f = frequency of vibration of the plate and L = wavelength.
If the wavelength L can be measured, the velocity of Rayleigh waves can easily
be calculated. The wavelength is generally determined by the number of waves
occurring at a given distance x. For a given frequency Ji the wavelength can be
given as
(4.87)
where n1 = number of waves at a distance x for frequency Ji

(as shown in Figure 4.52b).
It was shown in Chapter 3 that the Rayleigh wave velocity is approximately equal
to the shear wave velocity. So
(4.88)
It was also discussed in Chapter 3 that, for all practical purposes, the Rayleigh
wave travels through the soil within a depth of one wavelength. Hence for a given
-----'-------v____,_L'\
__
1 Dynamic load
---L---
Frequency = f
__"-7.....---L
\J_C'\ __.., X
(a)
Rayleigh wave length: L1 > L2

Frequency of vibration of plate: /1 < /2
Distance, x
(b)
Figure 4.52 Subsoil exploration by steady-state vibration
frequency f, if the wavelength Lis known, the value of vs determined by the pre
ceding technique will represent the soil conditions at an average depth of L/2.
Thus for a large value off, the value of Vs is representative of soil conditions at a
smaller depth; and, for a small value off, the value of Vs obtained is representa
tive of the soil conditions at a larger depth. Figure 4.53 shows the results of wave
propagation on a stratified pavement system obtained using this technique.
11D SPECTRAL ANALYSIS OF SURFACE WAVE (SASW)

The Spectral Analysis of Surface Wave (SASW) method is a nonintrusive tech
nique developed from the steady-state vibration. The fundamentals of SASW
method involve the generation, measurement, and analysis of Rayleigh waves
at the test site. The shear wave velocity (Vs) profile of the test site can be derived
from the measured Rayleigh wave velocity.
Figure 4.52 shows a traditional configuration of equipment adopted in
SASW testing with a two-channel recording system. The surface waves are
4.17 Spectral Analysis of Surface Wave (SASW)
1000 ..--------.-------,.---.........--.........-----.
800
/
, L=0.68 m
/
/
/
/
200
o�--___._--��--�---��
120 240 360 480
Velocity, vr (mis)
Frequency
(a)
-f� -!= = -=
Velocity, vr (mis) True
o----------------
120 240 360 480 Profile
0.3 - - _Q.�5_m_ - - - - -o =
0.6
0.9
I 1.2
1.1sm 0
------ --------
0
1.5
1.8 Clay
2.1
2.4 ._______,_____.._______,_____....____.
Depth
(b)
Figure 4.53 vr as a function of frequency and depth determined by the steady

state vibration technique (after Heukelom and Foster 1960)
Source: Heukelom, W., and Foster, C.R. (1960). "Dynamic Testing of Pavements," Journal of
the Soil Mechanics and Foundations Division, ASCE, Vol. 86, No. SMl, Part 1, pp. 1-28. With
permissi on from ASCE.
generated by a dynamic vertical load (such as by hammer or bulldozer) to the

ground surface. All the equipment in SASW test is deployed on the ground sur
face, and thus no boreholes are required. The time delay between the two receiv
ers is a function of wave frequency (j) and is determined as
(f
t(f) = </J ) (4.89)
2nf
Copyright 2017 Cengage Learning. All Rights Reserved. May not be copied. scanned. or duplicated. in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter{s).
where </J(/) = phase angle for a given frequency in radians. The surface wave
phase velocity ( v,) and the surface wave length (L) are calculated as
v, = d2 (4.90)
t(f)
v
L= , (4.91)
where d2 = distance between receivers. The v, and Lare calculated for each fre
quency at various receiver spacing d2 , and the resulting v, - L relationship is the
experimental dispersion curve (for example, see Figure 4.55). A dispersion curve
is the correlation between the surface wave velocity and the wavelength ( or the
equivalent frequency).
The dispersion curves are then used to determine the shear wave velocity
(Vs) profile. The Vs profile of the test site can be determined by iteratively fitting
the experiment dispersion curve with a theoretical dispersion curve, the latter
is derived from an assumed Vs profile. In this iterative procedure, the values of
shear wave velocity and the thicknesses of each soil layer in the assumed vs
profile are updated by trial and error until a satisfactory theoretical dispersion
curve that fits the experiment dispersion curve is obtained (for example, see
Figure 4.55). The corresponding vs profile for the best-fit theoretical dispersion
curve is determined to be the Vs profile for this test site. Figure 4.56 shows the
shear wave velocity profile for the best-fit theoretical dispersion curve shown in
Figure 4.55.
The SASW method is typically used to obtain shear wave velocity profiles for
earthquake site response, liquefaction analysis, soil compaction control, pave
ment evaluation, mapping subsurface stratigraphy, etc. The shear wave velocity
profiles at greater depth require a high-energy and low-frequency wave source,
while profiles with shallower depth need a low-energy and high-frequency wave
source.
Signal
analyzer
Vertical
receiver #2
--l•---d1---1--•---dz--�•I
Figure 4.54 Standard configuration of equipment in SASW testing with a

two-channel recording system
4.17 Spectral Analysis of Surface Wave (SASW)
500
- Experimental dispersion curve
o Theoretical dispersion curve
o ,_____._................................_.....__............................_..................................
0.1 1 10 100
Wavelength, L (m)
Figure 4.55 Comparison between experimental and theoretical dispersion

curves (adapted from Stokoe et al., 1994)
Source: Based on Stokoe, K.H. II, Wright, S.G., Bay, J.A., and Roesset, J.M. (1994). "Characterization
of Geotechnical Sites by SASW Method." Technical Review: Geophysical Characterization of Sites,
/SSMFE Technical Committee 10, edited by R .D. Woods, Oxford, New Delhi, pp. 11-25.
o --------.lc===�·1:_1----.l----r--l-----.-l -----.l
0 50
Shear wave velocity, vs (mis)
100 150 200 250 300
-
1-
2-
3-
4-
5-
6-
Figure 4.56 Shear wave velocity profile determined from SASW (adapted from
Stokoe et al., 1994)
ISSMFE Technical Committee 10, edited by R .D. Woods, Oxford, New Delhi, pp. 11-25.
IJD SOIL EXPLORATION BY .. SHOOTING UP THE HOLE;'

.. SHOOTING DOWN THE HOLE," AND .. CROSS-HOLE
SHOOTING"
Shooting Up the Hole
In the technique of shooting up the hole, a hole is drilled into the ground and a
detector is placed at the ground surf ace. Charges are exploded at various depths
in the hole and the direct travel time of body waves (P or S) along the boundary
of the hole is measured. Thus the values of v P and vs of various soil layers can
be easily obtained. There is a definite advantage in this technique, since it deter
mines the shear wave velocities of various soil layers. The refraction and reflec
tion techniques give only the P-wave velocity. However, below the groundwater
table, the compression waves will travel through water. The first arrival for points
below the water table will usually be for this type, and the wave velocity will gen
erally be higher than the compression wave velocity in soils. On the other hand,
shear waves cannot travel through water, and the shear wave velocity measure
above or below the water table will be the same.
Shooting Down the Hole

Shear wave velocity determination of various soil layers by shooting down the
hole has been described by Schwarz and Musser (1972), Beeston and McEvilly
(1977), and Larkin and Taylor (1979). Figure 4.57 shows a schematic diagram for
the down-hole method of seismic waves testing as presented by Larkin and Tay
lor, which relies on measuring the time interval for SH-waves to travel between
the ground surface and the subsurface points. A bidirectional impulsive source
for the propagation of SH-waves is placed on the surface adjacent to a borehole.
A horizontal sensitive transducer is located at a depth in the borehole. The depth
of the transducer in varied throughout the length of the borehole. The shear
wave velocity can then be obtained as
(4.92)
where z = depth below the ground surf ace

t = time of travel of the shear wave from the surf ace impulsive
source to the transducer
During the process of field investigation, Larkin and Taylor (1979) deter
mined that the shear strains at depths of 3 m and 50 m were about 1 X 1o- 6 and
0.3 X 1 o- 6, respectively. In order to compare the field and laboratory values of vs ,
some undisturbed samples from various depths were collected. The shear wave
4.18 Soil Exploration by "Shooting Up the Hole," "Shooting Down the Hole
Bidirectional
surface source Time to first shear
of SH-wave wave arrival, t
.· . _ - _ . _ .
dz
V =-
s dt
Transducer
Wave path
o--+------------------- ----------
Soil layer 2
•
Depth, z
Figure 4.57 The down-hole method of seismic wave testing
velocity of various specimens at a shear strain level of 1 X 10- 6 was determined.

A comparison of the laboratory and field test results showed that, for similar
soils, the value of Vs(Iab) is considerably lower than that obtained in the field. For
the range of soil tested,
Vs( l ab) :::::: 0.25 Vs(field) + 83
Where Vs(lab) and Vs (field)are in meters per Second.

Larkin and Taylor also defined a quantity called the sample disturbance fac
tor Sn:
2
Vs field)
Sn = [ ( ] (4.93)
Vs(lab)
The average value of Sn in Larkin and Taylor's investigation varied from about
1 for Vs(field) = 140 m/s to about 4 for Vs(fiel d) = 400 m/s. This shows that small dis
turbances in the sampling could introduce large errors in the evaluation of repre
sentative shear moduli of soils.
Cross-Hole Shooting
The seismic cross-hole survey is considered by many engineers to be the most
reliable method of determining the dynamic shear modulus of soil. The tech
nique of cross-hole shooting relies on the measurement of SV-wave velocity. In
this procedure of seismic surveying, two vertical boreholes at a given distance
apart are advanced into the ground (Figure 4.58). Shear waves are generated by
a vertical impact at the bottom of one borehole. The arrival of the body wave
is recorded by a vertically sensitive transducer placed at the bottom of another
borehole at the same depth. Thus
L
Vs
=
- (4.94)
where t = travel time for the shear wave

L = length between the two boreholes.
The smallest possible borehole diameter should be used as small uncased
boreholes are more stable than larger diameter holes. Even if casing is required,
a small diameter borehole will cause less soil disturbance. Usually aluminium or
PVC casing is used instead of steel casing. Void spaces around the casing must
be filled with weak cement slurry grout. Spacing between the boreholes can be 2
Oscilloscope Trigger
Input
Vertical
velocity transducer
. . ·:
·- .·
Impulse rod
Vertical
velocity ··. . .
transducer
� � . . ..._..._........._.
'Path of body wave '"
Figure 4.58 Schematic diagram of cross-hole seismic survey technique
4.20 Cyclic Plate Load Test
4 Soil profile
Silty sand
(SM)
0.9 w=3-9%
s 1.8
.sfr
'-"
Sandy silt
Q (ML)
2.7 w= 6--14%
3.6
Silty sand
w=5-7%
w = moisture content
Figure 4.59 Shear wave velocity versus depth from cross-hole seismic survey
(redrawn after Stokoe and Woods. 1972)
ISSMFE Technical Committee 10, edited by R .D. Woods, Oxford, New Delhi, pp. 11-25.
to 3 meters. The borehole spacing at the surface can be used as L, and for deeper
boreholes (say more than 10 m in depth), inclinometers must be used to calculate L
accurately as a small error in L can lead to large differences in estimated shear wave
velocity. Figure 4.59 shows the plot of the shear wave velocity against depth for a
test site obtained from the cross-hole shooting technique of seismic surveying.
mJ CORRELATIONS FOR SHEAR WAVE VELOCITY, Vs

Several correlations between the shear wave velocity Vs and field standard pene
tration number N have been presented in the past. A few of these correlations are
given in Table 4.4. Significant differences exist among the published relations that
may be due to differences in geology along with the measurement of N and vs.
IE:J CYCLIC PLATE LOAD TEST

The cyclic field plate load test is similar to the plate bearing test conducted in
the field for evaluation of the allowable bearing capacity of soil for foundation
design purposes. The plates used for tests in the field are usually made of steel
and are 25 mm thick and 150 mm to 762 mm in diameter.
Table 4.4 Some Correlations between vs (mis) and N.

Source Correlation
Imai (1977) All soils: Vs = 91N0.37?
Sand: Vs = 80.6N°.3 3 t
Clay: Vs = 80.2N° · 292
Ohta and Goto (1978) All soils: Vs = 85.35N°· 348
Seed and Idriss ( 1981) All soils: Vs = 61.4N° ·5
Sykora and Stokoe ( 1983) Sand: Vs = 100.5N°· 29
Okamoto et al. (1989) Sand Vs = 125N° .3
Pitilakis et al. (1999) Sand: Vs = 145N°·178
Clay: = l 32N° ·271
Vs
Kiku et al. (2001) All soils: Vs = 68.3N° · 292

Jafari et al. (2002) Sand: Vs = 22N°· 77
Clay: Vs = 27N°· 73
Hasancebi and Ulusay (2006) All soils: Vs = 99N° · 309
Sand: Vs = 90.82N°.3t 9
Clay: Vs =
97.89N° ·269
Dikman (2009) All soils: Vs = 58N° · 3s
Sand: Vs = 73N° .3�
Silt: Vs = 60N° ·36
Clay Vs = 44N ° ·48
To conduct a test, a hole is excavated to the desired depth. The plate is placed
at the center of the hole, and load is applied to the plate in steps-about one
fourth to one-fifth of estimated ultimate load-by a jack. Each step load is kept
constant until the settlement becomes negligible. The final settlement is recorded
by dial gauges. Then the load is removed and the plate is allowed to rebound. At
the end of the rebounding period, the settlement of the plate is recorded. Fol
lowing that, the load on the plate is increases to reach a magnitude of the next
proposed stage of loading. The process of settlement recording is then repeated.
Figure 4.60 shows the nature of the plot of q versus settlement (s) obtained
from a cyclic plate load test. Note that
load on the plate, Q
q=
area of the plate, A
Based on field test results, the magnitude of the spring constant k [See Chapter 2,
Eq. (2.3)] and the shear modulus G of the soil can be calculated in the following
Load per unit area, q

q2
--r-- ---------
8e(2)
''
''
t -----i..----
\
''
'
__ J_ \
\
--i--------------------
I
I
Figure 4.60 Nature of load settlement diagram for a cyclic plate load test
manner. It is worth noting that in order to accurately reflect the nonlinear response
of the soil, it would be necessary to establish the similar strains between the small
scale footing and prototype footing. A number of cycles of loading of the plate
may be needed to replicate the elastic condition in the soil under footing.
Spring Constant k
1. Referring to Figure 4.60, calculate the elastic settlement [se (l), se c2), ...] for
each loading stage.
2. Plot a graph of q versus se, as shown in Figure 4.61.
3. Calculate the spring constant of the plate as
qA
�late
=
(4.95)
Se
q4
....cl:S
Q)
.....
::s q3
....
Q)
cl:S
q2
qi
Elastic rebound, se
Figure 4.61
4. T he spring constant for vertical loading for a proposed foundation can then
be extrapolated as follows (Terzaghi, 1955).
Cohesive soil:
_ kplate (foundation width]

kfoundation - (4.96)
plate width
Cohesionless soil:
foundation width+ plate width

2
_ = kplate (
kfoundation . J (4.97)
2 X plate w1dth
Shear Modulus, G
It can be shown theoretically (Barkan, 1962) that
E -1-
CZ =_q__=l.13 (4.98)
Se 1-µ 2 fl
where Cz = subgrade modulus

E = modulus of elasticity
µ = Poisson's ratio
A = area of the plate
However,
G= E
2(1 + µ)
So
2.26G(l + µ) _1_
Cz =
l-µ2 JA
or
(1-µ)Czfl
G= (4.99)
2.26
The magnitude of Cz can be obtained from the plot of q versus se (Figure 4.61).
With the known value of A and a representative value of µ, the shear modulus
can be calculated from Eq. (4.99). In nonhomogenous soils, it may be desirable
to conduct the test at different depths or one may use different plate sizes to
reflect the change in soil stiffness with depth. Again, it should be noted that
this test suffers from the same limitations as reported in traditional geotechnical
engineering practice for the design of foundations.
EXAMPLE4.5
The plot of q versus s (settlement) obtained from a cyclic plate load test is
shown in Figure 4.62. The area of the plate used for the test was 0.3 m2 •
Calculate
a. ¾late' and
b. shear modulus G (assumeµ= 0.35).
SOLUTION
a. From Figure 4.62, the following can be determined.
Load per unit area, q (k.Pa) Elastic settlement, se (mm).

75 0.53
150 1.10
225 1.50
300 2.10
100 200 300 q (kPa)
15
20
Figure 4.62
Figure 4.63 shows a plot of q versus se . From the average plot,

300
C2 = !L = = 142.86 MN/m 3
Se 0.0021
From Eq. (4.95)
kplate = qA = (142.86)(0.3) = 42.86 MN/m.

Se
b. From Eq. (4.99)
(1- µ)C2 JA (1- 0.35)(142.86)(Jo.3)

G= =
2.26 2.26
� 22.5 MPa
This method can be extremely useful in sandy soils, provided it is being
preceded by a boring program. But this test may not give good results if a
weak stratum lies below the significant depth of test plate but within the
4.21 Test Procedures for Measurement of Moduli and Damping Characteristics
Figure 4.63
significant depth of the foundation. The procedure is costly, particularly if

the ground water level is near the foundation and ground water lowering
becomes necessary.
♪ CORRELACIONES PARA EL MÒDULO DE CORTE Y RELACIÒN DE AMORTIGUAMIENTO

Correlations for Shear Modulus and Damping Ratio
Im TEST PROCEDURES FOR MEASUREMENT OF MODULI
AND DAMPING CHARACTERISTICS
For design of machine foundations subjected to vibration, calculation of ground response
during an earthquake, analysis of the stability of slopes during an earthquake, and other
dynamic analysis of soil, it is required that the shear modulus and the damping ratio of
the soil be known. The shear modulus G and the damping ratio D of soils are dependent
on several factors, such as type of soil, confining pressure, level of dynamic strain, degree
of saturation, frequency, and number of cycles of dynamic load application, magnitude of
dynamic stress, and dynamic prestrain (Hardin and Black, 1968).
precarga dinàmica
From the preceding discussions in this chapter, it is obvious that a wide vari
ety of procedures, including laboratory and field tests, can be used to obtain
the shear moduli and damping characteristics of soils. A summary of those
test conditions, range of applicability, and the parameters obtained are given in
Table 4.5. Based on these studies several correlations for estimation of G and D
have evolved during the last 30 to 40 years. Some of these correlations are sum
marized in the following sections.
In general , the shear-stress-versus-shear strain relationship for soils will be
of the nature as shown in Figure 4.64. The following can be seen from this figure:
1. The shear modulus G decreases with the increased level of shear strain.
2. At a very low strain level, the magnitude of the shear modulus is maximum
(that is, G = Gmax ),
3. The shear stress-versus-shear-strain relationship shown in Figure 4.64 can be
approximated as (Hardin and Drnevich, 1972)
r=------ r' (4.100)

1/Gmax + r' /rmax
where r = shear stress and y' = shear strain.
Table 4.5 Test Procedures for Measuring Moduli and Damping Characteristics
General Approximate Properties
procedure Test condition strain range determined
Determination of hysteretic Triaxial compression 10- 2 to 5% Modulus; damping
stress-strain relationships Simple shear 10- 2 to 5% Modulus; damping
Torsional shear 10-2 to 5% Modulus; damping
Forced vibration Longitudinal vibrations 10-4 to 10-2 % Modulus; damping
Torsional vibrations 10-4 to 10-2 % Modulus; damping
Shear vibrations-lab 10-4 to 10-2 % Modulus; damping
Shear vibration-field 10-4 to 10- 2 % Modulus
Free vibration tests Longitudinal vibrations 10- 3 to 1% Modulus; damping
Torsional vibrations 10- 3 to 1% Modulus; damping
Shear vibration-lab 10- 3 to 1% Modulus; damping
Shear vibration-field 10- 3 to 1% Modulus
Field wave velocity Comparison waves =5 X 10-4 % Modulus
measurements Shear waves =5 x10-4 % Modulus
Rayleigh waves =5 X 10-4 % Modulus
Field seismic response Measurement of motions at Modulus; damping
different levels in deposit
After Seed and Idriss, 1970
Source: Seed, H.B., and Idriss, I.M. (1970). "Soil Moduli and Damping Factors for Dynamic Response
Analysis," Report No. EERC 70-10, Earthquake Engineering Research Center, University of Calofornia,
Berkley. Reprinted by permission of the PEER Center, UC Berkeley.
4.22 Shear Modulus and Damping Ratio in Sand
- - - - - - - - - -,- - - - - - - - - - - - - - - 7 - - - - - - - - - - - - - - - - - - - - - - - - - - -
I I
/ I
/ I
I
�
I
I
'f=---- r'
1 /
I
--+-
1 r'
G
I max 'rmax
I
G /
max
I
I
Shear strain, r'

Figure 4.64 Nature of variation of shear modulus with strain
IEJ SHEAR MODULUS AND DAMPING RATIO IN SAND

MÓDULO DE CORTE Y RELACIÓN DE AMORTIGUAMIENTO EN ARENAS
Hardin and Richart (1963) reported the results of several resonant column tests
conducted in dry Ottawa sand. The shear wave velocities Vs determined from
some of these tests are shown in Figure 4.18. The peak-to-peak shear strain
amplitude for these tests was 10- 3 rad. From Figure 4.18 it may be seen that the
values of vs are independent of the gradation, grain-size distribution, and also
the relative density of compaction. However, Vs is dependent on the void ratio e
and the effective confining aO and can be expressed by the following empirical
relations:
Vs =
(19.7 - 9.06 e) ab14 for ao > 95.8 kPa (4.101)
and
vs = (11.36 - 5.35 e) a8·3 for aO < 95.8 kPa (4.102)
In Eqs. (4.101) and (4.102), the units of Vs and a 0 are meters per second and
Newtons per square meter (Pa), respectively.
Several experimental results for shear wave velocity in extremely angular
crushed quartz sands were also reported by Hardin and Richart (1963). Based on
these results, the value of vs for angular sands can be expressed by the empirical
relation
Vs = (18.43 - 6.2e)O' �4
i i (4.103)
(m/s)
Based on the shear wave velocity relations presented here, the shear modulus
of sands for low amplitudes of vibration can be given by the following relations
(Hardin and Black, 196 8):
_ 6908 (2 .17 - e)2 - 112 (round-grained)

Gmax- (4.104)
O'o
l +e
and
Gmax = 3230(2·97 - e ) 6 o112 (angular-grained)

2
(4.105)
l+e
where Gmax and iiO are in kPa.

For a soil specimen subjected to a stress condition such that a1 # a2 # a 3
(where a 1 , 62 and 63 are the major, intermediate, and minor effective principal
stresses, respectively), the average effective confining pressure is
?io = ½ (a 1 + a2 + a3) = effective octahedral stress z
This value of a0 can be used in Eqs. (4.101)- (4.105).
For field conditions at any given depth,
o\ = effective vertical stress = av

if2 = 0'3 = Koav
x
where K0 = at-rest earth pressure coefficient� 1 - sin </) y
(where </)= drained friction angle).
So
a o = ½ [iiv + 26v (1 - sin</))]

(4.106)
= av (3 - 2 sin </))
3
Several investigators (e.g., Weissman and Hart, 1961; Richart, Hall, and
Lysmer(1962); Drnevich, Hall, and Richart, 1966; Silver and Seed, 1969; Hardin
and Drnevich, 1972; Seed and Idriss, 1970; Shibata and Soelarno, 1975; and Iwa
saki, Tatsuoka, and Takagi, 1976), have reported the results of shear modulus
and damping ratio measurements using various types of test techniques. From
these test results, it appears that the shear modulus at a given strain level can be
expressed as (Seed and Idriss, 1970)
(4.107)
where G and a O are in kPa.
For low strain amplitudes(r' < 10-4 %), the preceding equation will be
(4.108)
The magnitudes of K2 cmax) vary from about 30 for loose sands to about 75 for
dense sands. Seed and Idriss (1970) recommended the following values of K2 cmax)·
Relative density,
Rv (%) K2cmax>
30 34
40 40
45 43
60 52
75 61
90 70
Hence,
G
(4.109)
Figure 4.65 shows the variation of F' with shear strain y'(¾) obtained from sev
eral studies. These values fall in a rather narrow band and, for all practical pur
poses, the average plot can be used for design and estimation purposes. Thus Eqs.
(4.104), (4.105), (4.107), (4.108), and (4.109) can be combined to estimate the
shear modulus at any required shear strain level.
Studies by Hardin and Dmevich(1972) and Seed and Idriss(1970) show that the
damping ratios for sands are affected by factors such as(a) grain-size characteristics,
(b) degree of saturation, (c) void ratio, (d) earth pressure coefficient at rest( K0 ), (e)
angle of internal friction(¢), (f ) number of stress cycles (N'), (g) level of strain, and
(h) effective confining pressure. The last two factors, however, have the major effect
on the magnitude of the damping ratio. Figure 4.66 shows a compilation of past
studies(Seed et al. 1986) to determine D. For most practical cases the average plot of
the variation of D versus r' can be used for most calculation purposes.
0.8 i-------------r-------,--------..:----1---;,......_---------t------------i
Average
F'
Range of F' -+--�
0.3
o.__________,_________....__________.___________.
10--4 10-3 10-2 10-1 1
Shear strain, y' (%)
Figure 4.65 Variation of F' with shear strain for sands (after Seed et al., 1986)
Source: SourceSeed, H.B., Wong, R.T., Idriss, I.M., and Tokimatsu, K. (1986). "Moduli and
Damping Factors for Dynamic Analyses of Cohesionless Soils," Journal of Geotechnical
Engineering, ASCE, Vol. 112, No. GT l 1, pp. 1016-1032. With permission from ASCE.
28
/
/
1-----------------------------,.-/------,,-
/
/
24
/
/
/
/
/
20
/ /
/
,...._
/
/ /
/ /
/
/
·-
/
Q
/
/
/
16
/
/
0� /
/
/
/ /
/ /
/ /
·a
/ /
s
/ /
12 / /
/
/
/ /
/
/
8
/
/
/
/
/
/
/
/
/
/
Average plot
--- (Gráfica promedio)
0 ._________...,___________._________.......__________.
10--4 10-3 10-2 10-1 1
Shear strain, y' (%)
Figure 4.66 Damping ratios for sands (after Seed et al.,1986)

Source: Seed, H.B., Wong, R.T., Idriss, I.M., and Tokimatsu, K. (1986). "Moduli and Damping
Factors for Dynamic Analyses of Cohesionless Soils," Journal of Geotechnical Engineering,
ASCE, Vol. 112, No. GT ll, pp. 1016-1032. W ith permission from ASCE.
Based on tests on dry sands using a torsional simple shear device, Sherif,
Ishibashi, and Gaddah (1977) proposed the following relationship for damping ratio.
D = 5o - 0.087 <io (73.3F - 53.3)(y')0·3 (1.01- 0.046 logN) (4.110)
38
where D = damping ratio (%)
<Jo = effective confining pressure (kPa)
r' = shear strain (%)
F = sphericity factor of the soil grains
N' = number of cycles of strain application
The sphericity factor is defined as
(4.111)
where
D}o C = (4.112)
g (D1o)(D60)
Dw , D30 , D60 = diameters, respectively, through which
10%, 30%, and 60% of the soil will pass
S'
1/f =-
s
where S' and S are, respectively, the surface area of a sphere of the same volume
as the soil particle and the actual surface area of the soil.
EXAMPLE4.6
The groundwater table in a normally consolidated sand layer is located at a depth
of 3m below the ground surface. The unit weight of sand above the groundwater
table is 15.5 kN/m3 . Below the groundwater table, the saturated unit weight of
sand is 18.5 kN/m 3 • Assuming that the void ratio and effective angle of friction
of sand below the groundwater table are 0.6 and 36° , respectively, determine the
damping ratio and the shear modulus of this sand at a depth of 7.5 m below the
ground surface if the strain is expected to be about 0.12%.
SOLUTION
From Eq. (4.106)
<J
<JO =v ( 3- 2 sin ¢)
3
<Jv = 3(15.5) + 4.5(18.5 - 9.81) = 85.61 kPa
85 61
<Jo = · [3- (2)(sin 36)] = 52.06 kPa
3
When <p is equal to 36° , Rn is about 40 to 50%. Assuming Rn :=:::: 45%, K2 (max) :::::: 43.
So, from Eq. (4.108)
or
Gmax = (218.82)(43)(52.06) 0·5
= 67,890 kPa :::::: 67.9 MPa
Referring to Figure 4.65, for r' = 0.12%, the value of F' is about 0.28. So
G = F'Gmax = (0.28)(67.9) :::::: 19 MPa
Referring to the average curve in Figure 4.66, for r' = 0.12%
»�17%
mJ CORRELATION OF Gmax OF SAND WITH STANDARD

PENETRATION RESISTANCE
The standard penetration test is used in soil-exploration programs in the United
States and other countries. In granular soils the standard penetration numbers
(N in blows/0.3 m) are widely used for the design of foundation. The standard
penetration number can be correlated (Seed et al., 1986) in the following form to
predict the maximum shear modulus:
t t (4.113)
(kPa) (kPa)
where = effective confining pressure (kPa)

CJo
N60 = N-value measured in SPT delivering 60% of the theoretical
free-fall energy to the drill rod
Equation (4.113) is very useful in predicting the variation of the maximum shear
modulus with depth for a granular soil deposit.
GRUPO A/29/04/21 GRUPO A/04/05/21
EID SHEAR MODULUS AND DAMPING RATIO FOR GRAVELS
Seed et al. (1986) provided the experimental results of several well-graded gravels.
An example of such a study on well-graded Oroville material is shown in Figure
4.67. Based on several studies of this type, Seed et al. concluded that Eqs. (4.108)
Repetido aquí por conveniencia
4.24 Shear Modulus and Damping Ratio for Gravels
Repetido aquí por conveniencia
and (4.109) can also be used to predict the variation of shear modulus with shear
strain. However, the magnitude of K2cmax) for gravels ranges between 80 to 180 (as
compared to a range of 30 to 75 for sand). Thus,
G = GmaxF' = 218.82F' K2(max){CY 0 ) 0·5

i i (4.114)
Para amplitudes bajas de deformación
unitaria se transforma en la ec(4.108) (kPa) (kPa)
The variation of F' with the level of shear strain is shown in Figure 4.68.
The equivalent damping ratio of gravelly soils determined in the laboratory
from the hysteresis loops at the fifth cycle of each strain amplitude is shown in
Figure 4.69. It can be seen that, for a given value of y', the equivalent damping
ratio increases with the increase of the relative density Rn of the gravel. Seed
et al. (1986) also observed that
a. there is not significant effect of gradation on the equivalent damping ratios
of gravelly soil, and
Q.________.,____________._________._______,
10-4 10-3 10-2 10-1 1
Cyclic shear strain, r' (%)
Figure 4.67 Shear moduli of well-graded Oroville material (after Seed et al., 1986)
ASCE, Vol. 112, No. GTl 1, pp. 1016-1032. With permission from ASCE.
1.0
0.8
plot
( gráfica promedio)
0.6
F'
0.4
0.2
0
10-4 3 10-3 3 10-2 3 10-1 3 1
Cyclic shear strain± r' (%)
Figure 4.68 Variation of F' with shear strain for gravelly soils (after Seed
et al., 1986)
ASCE, Vol. 112, No. GTll , pp. 1016-1032. W ith permission from ASCE.
b. the damping ratio is not significantly affected by the number of cycles at very
small strain amplitudes. However, it decreases to approximately three-fourths
of its original value after 60 cycles at any axial strain amplitude of+ 0.2%.
Seed et al. showed that the range and the average plot of the damping ratio D
with strain amplitude y' for gravelly soils is approximately the same as that for
sands (Figure 4.66).
24 �----------------�
-.. Rn = 100%
�
'-' 80%
.g 16 1------+-------+-------,--_65_3/c---10
f:! 15
·a�
0.1) 14
13
12
11
10
9
] 8t------+-------+,________---I
�
·3a<
�
0 .________._________._______._______.
10-4 10-3 10-2 10-1 1
Shear strain, (%) r'
Figure 4.69 Effect of relative density on damping ratio of gravelly soils (after
Seeds et al., 1986)
ASCE, Vol. 112, No. GTll , pp. 1016-1032. W ith permission from ASCE.
4.24 Shear Modulus and Damping Ratio for Gravels
Rollins et al. (1998) analyzed the results of several investigators between 1986
and 1998 which were obtained from cyclic triaxial tests and large diameter cyclic
torsional simple shear tests. Based on this analysis it was suggested that the rela
tionships for the best fit curves can be given as
G - ----------
-- 1
(4.115)
Gmax [1.2 + 16y'(l + 10- 20r')]
and
-0-9 -0·75
D(%) = 0.8 + 18(1 + 0.15y' ) (4.116)
where y' = cyclic shear strain (%)

The variations of G/Gmax and D based on Eqs. (4.115) and (4.116), along with
the standard deviation bounds, are given in Figures 4.70 and 4.71, respectively.
1.0 f"llll:��-""""'""1-:;;;::::::--------,--------,-----,
0.8
media
- - - - - - Eq. (4 .115)
Standard deviation bounds
0.2
o��-��-��-��--
10- 4 10-3 10-2 10-1
Figure 4.70 Variation of G/Gmax with y' for gravel-Eq. (4.115)

2s.-------.....-----..------...-------.
,--...
�
20
--
- - - - - Eq. (4.116)
Standard deviation bounds
�
�
.s� 15
/
/
lo.
,,
·t
�
10 /
/
/
,,
/
/
5
medio
----
0
10-4 10-3 10-2 10-1 100
Cyclic shear strain, r' (%)
Figure 4.71 Variation of damping ratio with r' for gravel-Eq. (4.116)
EXAMPLE 4.7
Consider a gravel deposit. At a certain depth below the ground surface, given:
Vertical effective stress, ii'v = 7 6 kPa

Friction angle, </J = 38 °
Relative density, Rn= 80%
Estimate the shear modulus and damping ratio at a cyclic shear strain level of
10-2 %. Use Eqs. (4.115) and (4.116).
SOLUTION
From Eq. (4.106)
o-0 = �' (3- 2sin¢) = (7 3

6
)c3- 2sin3 8° ) = 44.81 kN/m2
From Eq. (4.114)
From Figure 4.67 for y'::::::: 10- 4 %, the value of K2 ::::::: 116. Hence,
Gmax = (218.82)(116)( 44.81)0·5 = 169,915 kPa ::::::: 169.92 MPa
From Eq. (4.115)
G - --------
-- 1 l
::::::: 0.685
Gmax [1.2 + 16y'(l + 10-20r')] [1.2 + (16)(0.01)(1 + 1o-2oxO.Ol )]
Hence,
G = (169.92)(0.685) = 116.4 MPa
From Eq. (4.116)
D(%) = 0.8 + 18(1 + 0.15y'- 0·9)-0·75 = 0.8 + 18[1 + (0.15)(0.01) - 0·9 J- 0·75 = 3.89%
4.25 Shear Modulus and Damping Ratio for Clays
EID SHEAR MODULUS AND DAMPING

RATIO FOR CLAYS
Hardin and Black (1968) and Hardin (1978) proposed the following empirical
relationship for the shear modulus of clays at low amplitudes of strain, which
includes the effects of soil plasticity and the overconsolidation ratio (OCR). Or
OCRK
Gmax =625---,Jp aa o (4.11 7)
0.3 + 0.7e 2
where e = void ratio

Pa = atmospheric pressure expressed in the same units (:::::::: 100 kPa) as Gmax
K = f(plasticity index, PI )
Following are the recommended values of K for use in the preceding equation.
Plasticity index, PI (%) K

0 0
20 0.18
40 0.30
60 0.41
80 0.48
�100 0.5
For field conditions
fio = ½(fiv + 2Ko fiv ) (4.11 8)

where av = effective vertical stress
K0 = at-rest earth pressure coefficient
For normally consolidated clays (Booker and Ireland, 1965)
K0 = 0.4 + 0.007 (PI) (for O <PI< 40%) (4.119)
and
K0 = 0.68 + 0.001(PI - 40) (for 40% <PI< 80%) (4.120)
For over consolidated soils, K0 can be approximated as (Mayne and Kulhawy,
1982)
sin¢
Ko =(1- sin¢) OCR (4.121 )
1.0
0.8
0.6
�
(!JS
0.4
OCR= 1 to 15
0.2
0
10-4 10-3 10-2 10-1 10
25 Pl =0
15
20
OCR= 1 to 8 30
�
.....0 15
50
..... 100
10
200
0
10-4 10-3 10-2 10-I 10
Figure 4.72 Test results of Vucetic and Dobry (1991) - Variation of G/Gmax
andD
Vucetic and Dobry (1991) used a large database and provided the variation
of G/Gmax and D at various cyclic shear strain levels which are functions of PI
and OCR. These variations are shown in Figure 4. 72.
Correlation of Seed and Idriss

Seed and Idriss (1970) collected the experimental results for shear modulus and
damping ratio from various sources for saturated cohesive soils. Based on these
results the variation of G/cu (where cu =undrained cohesion) with shear strain
is shown in Figure 4. 73. Also, Figure 4. 74 shows the upper limit, average, and
lower limit for the damping ratio at various strain levels.
3000 .--------.----=::,-------,------,,---,------,--------.------.------,
Average
- 300
-
�l 'U
";:!.
100
10 ---�---�--�---------------�------�
10--4 10-3 10-2 10-1 1 3
Shear strain (%)
Figure 4.73 In situ shear modulus for saturated clays (after Seed and Idriss, 1970)
Source: Seed, H.B., and Idriss, I.M. (1970). "Soil Moduli and Damping Factors for Dynamic
Response Analysis," Report No. EERC 70-10, Earthquake Engineering Research Center,
University of Calofornia, Berkley. Reprinted by permission of the PEER Center, UC Berkeley.
35 ..--------,-------.----------,-------...--------,
Average
0--------------------------�-----�
10--4 10-3 10-2 10-1 1 10
Shear strain (%)
Figure 4.74 Damping ratio for saturated clays (after Seed and Idriss, 1970)
Source: Seed, H.B., and Idriss, I.M. (1970). "Soil Moduli and Damping Factors for Dynamic
Response Analysis," Report No. EERC 70-10, Earthquake Engineering Research Center,
University of Calofornia, Berkley. Reprinted by permission of the PEER Center, UC Berkeley.
EXAMPLE 4.8
A soil profile is shown in Figure 4. 75a. Calculate and plot the variation of
- ·r -
shear modulus with depth (for low amplitude of vibration) .
.·<Dry sand <

. .. . . . ··as;;; 2.65.
. ... �....
I Layer I, 3-r m e=0.7

°
¢>=30 .a.w.T.
Sand Not to
II Layer II, 1.5 m Gs= 2.65 scale
e=0.6
°
¢>=33
Normally consolidated clay
Gs= 2.78
III Layer III, 1.5 m e= 1.22
Liquid limit=48
Plastic limit =23
Rock
(a)
Figure 4. 75a
SOLUTION
At any depth z
where K0 is the coefficient of earth pressure at rest and o\ is the vertical

effective pressure. For sands,
Ko = 1- sin</)= 1- sin 30° = 0.5
In layer II,
K0 = l - sin </)= 1 - sin 33° = 0.455
For normally consolidated clays,
Ko = 0.4 + 0.007(PI) for 0<PI< 40%

= 0.4 + 0.007( 48- 23) = 0.575
Calculation of Effective Unit Weights

z = 0-3.0m
Gsrw (2.65 ) (9.81) = 15_ 2 3

r =
r eff
= =
9 kN/m
1+e l + 0. 7
dry
z = 3.0-4.5m
(Gs e) rw Gs - l
Yetr rsat - rw - rw ( J w
+
l +e l +e Y
= = =
(2.65-1)(9.81) = 10.12 kN/m3

1 +0.6
=
z = 4.5 -6.0m:
(Gs -l)rw (2.78-1)(9.81)

Yetr 7_87 kN/m3
1 +e l + 1.22
= = =
The following table can now be prepared.
Depth z - -
0'1 ii2 = <13 =Koii1 O'o e G=Gmax
(m) (kPa) (kPa) (kPa) (MPa)
0 0 0 0 0.7 0
1.5 15.29 X 1.5 = 22.94 11.47 15.29 0.7 34.348
3.0 15.29 X 3
(in layer I) =45.87 22.94 30.58 0.7 48.568
3.0 45.87 20.87 29.20 0.6 57.51 8
(in layer 11)
4.5 45.87 + 10.12 X 1.5
(in layer 11) =61.05 27.78 38.87 0.6 66.358
4.5 61.05 35.10 43.75 1.22 30.Slh
(in layer III)
61.05 + 7.87 X 1.5
6.0 =72.86 41.89 52.21 1.22 33.65b
aEq. (4.104)
hEq. (4.117). Note OCR = 1
The variation of G = Gmax with depth is plotted in Figure 4.69b.
G(MPa)
0
0 ...,..----,----r-----.---
28
84 56
1.5
4.5
6.0
(b)
Figure 4. 75b
EXAMPLE 4.9
At a given depth in a saturated clay layer, given:
Vertical effective stress, cr v = 30 kPa

Soil friction angle, ¢ = 28°
Liquid limit = 47
Plastic limit = 27
Void ratio, e = 0.92
Overconsolidation ratio = 2
Estimate the shear modulus and damping ratio at a cyclic shear strain of
0.01 %. Use Eq. (4.117) and Figure 4.72.
SOLUTION
From Eq. (4.121)
Ko = (1- sin<j>)OCR <!> = (1- sin28)(2)sin 28 = 0.734

sin
O"v 30
O"o = !( O"v + 2Ko O"v ) = (1 + 2Ko ) = (1 + 2 X 0.734) = 24.68 kPa
3 3 3
4.26 Shear Modulus and Damping Ratio for Lightly Cemented Sand
FromEq. (4.117)
OCR K
Gmax = 625
0.3 + 0.7e 2
,JPaCJo
Plasticity index = 47 - 27 = 20. Hence K = 0.18. Thus
(2)0.18
= (625)-----,j(l00)(24. 8) = 39,412 kPa � 39.4 MPa
Gmax 6
0.3 + (0.7)(0.92)2
From Figure 4.72, for y' = 0.01%
_Q_ = 0.8 4
Gmax
Or
G = (0.8 4)(39.4) � 33.1 MPa
Again, from Figure 4.72 for y' = 0.01% and PI = 20
D�4%
l!Im SHEAR MODULUS AND DAMPING RATIO

FOR LIGHTLY CEMENTED SAND
Lightly cemented sand deposits are encountered in many parts of the world. The
cementing material in the sand deposits is primarily calcium carbonate. More
recently, the results of several research projects relating to the properties of
lightly cemented sands have been published. From these studies it appears that
the behavior of lightly cemented sands can be duplicated in the laboratory by
mixing sand and Portland cement in required properties. The maximum sheer
modulus can be expressed as (Saxena, Avramidis, and Reddy, 1988)
Gmax(CS) = Gmax(S) + �Gmax(C) (4.122)
where Gmaxccs) maximum shear modulus of lightly cemented sand
=
= maximum shear modulus of sand alone
Gmax(s)
= increase of maximum shear modulus due to cementa
�Gmax(C)
tion effect
According to Saxena, Avramidis, and Reddy, the magnitudes of Gmax(S) and
�Gmax(C) can be obtained from the following empirical relationships.
428.2
Gmax(S) = ( )o.426(cr 0 )o.s14
0_3 + O. 7 e Pa
2
i i i (4.123)
(kPa) (kPa ) (kPa)
where 0,, = atmospheric pressure in the same units as Gmax(S)

- 0.515e- 0.13ee+ 0.285
/1Gmax(e) ----- 172
= (CC) 0 .88 ( <Jo )
Pa (e - 0.5168) Pa (4.124)
(for CC< 2%)
/1Gmax(e) 77 3 - 0 .698e- 0 .04ee- 0 .2

= -(CC)l.2 ( <Jo )
Pa e Pa (4.125)
(for 2% <CC<8%)
where CC = cement content (in percent ) and e = void ratio.
When using Eqs. (4.124) and (4.125), the units of Gmax(S) , Pa , and aO need to be
consistent.
The damping ratio at low strain amplitudes (r' <10- 3 %) can be expressed as
(Saxena, Avramidis, and Reddy, 1988)
Des = Ds +We (4.12 6 )
where Des = damping ratio of cemented sand (%)
Ds = damping ratio of sand alone (%)
We = increase in the damping ratio due to cementation effect
-o.Js
Ds = 0.94(� ) (4.127)
-0.36
MJc = 0.49( CC)L07 (:: ) (4.128)
where CC= cement content (in percent). The units of Pa and cr0 need to
be consistent.
EXAMPLE 4.10
If a lightly cemented sand specimen is subjected to an effective confining pres
sure of 98 kPa, estimate the value of Gmaxces), given e = 0.7 and CC= 3%.
SOLUTION
From Eq. (4.123),
428.2
Gmax(S) = ( )o.426 (a0 )o.s14
0_3 + O.7e Pa
2
Problems
Given e = 0.7, Pa= 100 kPa, and CJ'o = 98 kPa,

428.2
Gmax(S) = 2 (1OO)o.426(98)0.s14
0.3 + (0.7)(0.7)
= 65,805 kPa = 0.066 GPa
From Eq. (4.125),
�Gmax(C ) 773 - 0.698e-0.04CC -0.2
= _ ( CC)l.2 (O"0 J
Pa e Pa
or
(0 )-0 04 3)-0 2
�Gmax(C ) -773 1.2 _ 98 o.698 .1 . ( .
= (3) ( )
100 0.7 100
= (1104.3)(3.737)(0.997) = 4114.39
�Gmax(C) = 411,439 kPa = 0.411 GPa
So
Gmax (CS) = Gmax(S) + �Gmax(C )
= 0.066 + 0.411 = 0.477 GPa = 477 MPa
PROBLEMS
4.1 A uniformly graded dry sand specimen was tested in a resonant column
device. The shear wave velocity Vs determined by torsional vibration
of the specimen was 231.65 mis. The longitudinal wave velocity deter
mined by using a similar specimen was 387.40 mis. Determine each of
the following.
a. Poisson's ratio
b. Modulus of elasticity (E) and shear modulus ( G) if the void ratio
and the specific gravity of soil solids of the specimen were 0.5 and
2.65, respectively.
4.2 A clayey soil specimen was tested in a resonant column device (torsional
vibration; free-free end condition) for determination of shear modulus.
Given: length of specimen = 90 mm, diameter of specimen = 35.6 mm,
mass of specimen = 170 g, frequency at normal mode of vibration
(n = 1) = 790 Hz. Determine the shear modulus of the specimen in kPa.
4.3 The Poisson's ratio for the clay specimen described in Problem 4.2
is 0.52. If a similar specimen is vibrated longitudinally in a resonant
column device (free-free end condition), what would be its frequency at

normal mode of vibration ( n = l )?
4.4 The results of a refraction survey in terms of time of first arrival (in mil
liseconds) and distance in meters is given below in tabular form. Assum
ing that the soil layers are perfectly horizontal, determine the P-wave
velocities of the underlying soil layers and the thickness of the top layer.
Distance Time of first arrival
(m) (ms)
7.5 49.08
15.0 81.96
23.0 122.8
30.5 148.2
45.5 174.2
61.0 202.8
76.0 228.6
91.5 256.7
4.5 Repeat Problem 4.4 for the following.

Comment regarding the material encountered in the second layer.
Distance Time of first Distance Time of first

(m) arrival (ms) (m) arrival (ms)
10 19.23 100 125.82
20 38.40 150 138.72
30 57.71 200 152.61
40 76.90 250 166.81
60 115.40 300 178.31
80 120.71
4.6 Repeat Problem 4.4 with the following results. Also determine the thick-
ness of the second layer of soil encountered.
Distance Time of first Distance Time of first

(m) arrival (ms) (m) arrival (ms)
10 41.66 60 119.21
15 62.51 70 128.11
20 83.37 80 136.22
30 91.82 90 141.00
40 101.22 100 143.81
50 110.16 120 152.00
Problems
4. 7 The results of a reflection survey are given here. Determine the velocity
of P-waves in the top layer and its thickness.
Distance from Time for first arrival of

shot point (m) reflected wave (ms)
10 32.5
20 39.05
30 48.02
40 58.3
60 80.78
100 128.55
4.8 Ref er to Figure 4.46 for the results of the following refraction survey:
Distance from point of Time of Distance from Time of

disturbance, first arrival point of disturbance, first arrival
A (m) (ms) E(m) (ms)
0 0 0 0
6.0 20 6.0 20
12.0 40 12.0 40.1
18.0 60 18.0 59.8
24.5 78.2 24.5 79.7
36.5 92.8 36.5 121.0
61.0 122.2 61.0 167.2
85.5 149.8 85.5 175.1
Point E 110.0 177.9 Point A 110.0 180.2
Determine:
a. the P-wave velocities in the two layers,
b. z' and z", and
c. the angle f3.
4.9 For a reflection survey refer to Figure 4.51, in which A is the shot point. Distance
AC= AE = 180 m. The times for arrival of the first reflected wave at points C and
E are 45.0 ms and 64.1 ms, respectively. If the P-wave velocity in layer 1 is 280 mis,
determined f3 and z'.
4.10 The results of a subsoil exploration by steady-state vibration technique are given
here (Section 4.16):
Distance from the plate Number of waves Frequency of vibration

vibrated x (m) per second of the plate (Hz)
10 41.00 900
10 18.00 400
10 9.00 200
10 4.55 100
10 2.65 90
10 2.30 75
10 1.77 60
10 1.47 50
Make necessary calculations and plot the variation of the wave velocity
with depth.
4.11 Figure P4.11 shows a soil profile with the standard penetration resis
tance (N) values with depth. Using the relationships given by Imai
(1977; see Table 4.4), estimate the variation of the shear wave velocity
(Vs ) with depth.
N
. • .
. . .
1.5 0 8
3.0 0 10
Sand
4.5 0 12
7.5 o 9 Clay
9.0 8
Depth (m)
Figure P4.11
4.12 An angular-grained sand has maximum and minimum void ratios of

1.1 and 0.55, respectively. Using Eq. (4.105), determine and plot the
variation of maximum shear modulus Gmax versus relative density
(Rn = 0 -100%) for mean confining pressures of 50, 100, 150, 200 and
300 kPa.
4.13 A 20-m-thick sand layer in the field is underlain by rock. The groundwa
ter table is located at a depth of 5 m measured from the ground surface.
Problems
Determine the maximum shear modulus of this sand at a depth of 10 m

below the ground surface. Given: void ratio = 0.6, specific gravity of
soil solids = 2.68, angle of friction of sand= 36° . Assume the sand to be
round-grained. Use Eq. (4.104)
4.14 For a deposit of sand, at a certain depth in the field the effective vertical
pressure is 120 kPa. The void ratio and the relative density are 0.72 and
30° respectively. Determine the shear modulus and damping ratio for a
shear strain levels of 5 X 10-2 %. Use eqs. (4.105) and (47.109)
4.15 A remolded clay specimen was consolidated by a hydrostatic pressure
of 205 kPa. The specimen was then allowed to swell under a hydrostatic
pressure of 105 kPa. The void ratio at the end of swelling was 0.8. If
this clay is subjected to a torsional vibration in a resonant column test,
what would be its maximum shear modulus ( Gmax )? These liquid and
plastic limits of the clay are 58 and 28, respectively.
4.16 Refer to the overconsolidated soil specimen in Problem 4.15. Estimate
the shear modulus and damping ratio of the specimen at a cyclic shear
strain of 0.01%.
4.17 Ref er to F igure P4.17. Given:
H1 = 2 m Gs(l) = 2.68
H2 = 8 m Gs(2) = 2.65
H3 = 3m </>i = 35°
e1 = 0.6 </>2 = 30°
e2 = 0.7 PI of clay = 32
Estimate and plot the variation of the maximum shear modulus (Gmax)
with depth for the soil profile
Layer I
YG.W.T.
Sand
Gs(2) = 2.65
H2
e2 = 0.7
Layer II
</Ji=30
°
l
Normally consolidated clay
Gs(J) = 2.73
H3 Moisture content= 50% Layer III
Plasticity index = PI = 32
Rock
Figure P4.17
4.18 Repeat Problem 4.17 given
H1 =H2 =H3 =6 m Gs(l) = Gs(2) = 2.66

e1 = 0.88 </J1 = 28°
e2 = 0.68 </J2 = 32
°
PI of clay = 20
4.19 The unit weight of a sand deposit is 16.98 kN/m 3 at a relative density of
60%. Assume that, for this sand
</J = 30 + 0.15Rn
where </J is the drained friction angle and Rn is the relative density (in
percent). At a depth of 6.09 m below the ground surface, estimate its
shear modulus and damping ratio at a shear strain level of 0.01%. Use
the equations proposed by Seed and Idriss (1970).
4.20 The results of a standard unconsolidated undrained triaxial test on a
undisturbed saturated clay specimen are as follows:
Confining pressure = 70 kPa

Total axial stress at failure = 166.6 kPa
Using the method proposed by Seed and Idriss (1970), determine and
plot the variation of shear modulus and damping ratio with shear strain
(strain range 10-3 % to 1%).
4.21 For Example 4.10, determine the damping ratio of the cemented sand.
References
Barkan, D. D. (1962), Dynamics of Bases and Foundations, McGraw-Hill Book Company,
New York.
Beeston, H. E., and McEvilly, T. V. (1977). "Shear Wave Velocities from Down Hole Measure
ments," Journal of the International Association for Earthquake Engineering, Vol. 5, No. 2,
pp. 181-190.
Booker, E. W., and Ireland, H. 0. (1965). "Earth Pressure at Rest Related to Stress History," Cana
dian Geotechnical Journal, Vol. 2, No. 1, pp. 1-15.
Carroll, W. F. (1963), "Dynamic Bearing Capacity of Soils. Vertical Displacements of Spread
Footing on Clay: Static and Impulsive Loadings," Technical Report No. 3-599, Report 5, U.S.
Army Corps of Engineers, Waterways Experiment Station, Vicksburg, Mississippi.
Casagrande, A., and Shannon, W. L. (1949). "Strength of Soils under Dynamic Loads," Transac
tions, ASCE, Vol. 114, pp. 755-772.
Das, B. M. (1990). Principles of Geotechnical Engineering, 2nd ed., PW S-KENT, Boston.
References
Dikmen, U. (2009). Statistical Correlations of Shear Wave Velocity and Penetration Resistance for
Soils," Journal of Geophysics and Engineering, Vol. 6, No. 1, pp. 61-72.
Dmevich, V. P. (1972). "Undrained Cyclic Shear of Saturated Sand," Journal of the Soil Mechanics
and Foundations Division, ASCE, Vol. 98, No. SM8, pp. 807-825.
Dmevich, V. P., Hall, J. R., Jr., and Richart F. E., Jr. (1966). "Large Amplitude Vibration Effects
on the Shear Modulus of Sand," University of Michigan Report to Waterways Experiment
Station, Corps of Engineers, US Army Contact DA-22-079-Eng-340, October 1966.
Drnevich, V. P., Hall, J.R., Jr., and Richart, F. E., Jr. (1967). "Effects of Amplitude of Vibration
on the Shear Modulus of Sand," Proceedings, International Symposium on Wave Propagation
and Dynamic Properties of Earth Materials, Ed. G. E. Triandafilidis, University of New Mex
ico Press, Albuquerque, New Mexico, pp. 189-199.
Dyvik, R., and Madshus, C. (1985). "Lab Measurements of Gmax Using Bender Elements," Pro
ceedings, Advances in the Art of Testing Soils under Cyclic Conditions, ASCE Annual Con
vention, Detroit, Michigan, pp. 186-196.
Hall, J. R., Jr., and Richart, F. E., Jr. (1963). "Dissipation of Elastic Wave Energy in Granular
Soils," Journal of the Soil Mechanics and Foundations Division, ASCE, Vol. 89, No. SM6,
pp. 27-56.
Hardin, B. 0. (1978). "The Nature of Stress-Strain Behavior for Soils," Proceedings, Specialty
Conference on Earthquake Engineering and Soil Dynamics, ASCE, Pasadena, California,
Vol. 1, pp. 3-90.
Hardin, B. 0. (1965). "The Nature of Damping of Sands," Journal of the Soil Mechanics and Foun
dations Division, ASCE, Vol. 91, No. SMl, pp. 63-97.
Hardin, B. 0., and Black, W. L . (1968). "Vibration Modulus of Normally Consolidated Clay,"
Journal of the Soil Mechanics and Foundations Division, ASCE, Vol. 94, No. SM2, pp. 353-369.
Hardin, B. 0., and Drnevich, V. P. (1972). "Shear Modulus and Damping in Soils: Design Equa
tions and Curves," Journal of Soil Mechanics and Foundations Division, ASCE, Vol. 98, No.
SM7, pp. 667-692.
Hardin, B. 0., and Richart, F. E., Jr. (1963). "Elastic Wave Velocities in Granular Soils," Journal of
the Soil Mechanics and Foundations Division, ASCE, Vol. 89, No. SMl , pp. 33-65.
Hasancebi, N., and Ulusay, R. (2007). "Empirical Correlations between Shear Wave Velocity and
Penetration Resistance for Ground Shaking Assessments," Bulletin of Engineering Geology
and the Environment," Vol. 66, No. 2, pp. 203-213.
Heukelom, W., and Foster, C. R. (1960). "Dynamic Testing of Pavements," Journal of the Soil
Mechanics and Foundation Division, ASCE, Vol. 86, No. SMl, Part 1, pp. 1-28.
Iida, K. (1938). "The Velocity of Elastic Waves in Sand," Bulletin of the Earthquake Research Insti
tute, Tokyo Imperial University, Vol. 16, pp. 131-144.
Iida, K. (1940). "On the Elastic Properties of Soil Particularly in Relation to Its Water Con
tents," Bulletin of the Earthquake Research Institute, Tokyo Imperial University, Vol. 18,
pp. 675-690.
Imai, T. (1977). "P- and S-Wave Velocities of the Ground in Japan," Proceedings, 9th Interna
tional Conference on Soil Mechanics and Foundation Engineering, Tokyo, Japan, Vol. 2,
pp. 257-260.
Ishibashi, I., and Sherif, M. A. (1974). "Soil Liquefaction by Torsional Simple Shear Device," Jour
nal of the Geotechnical Engineering Division, ASCE, Vol. 100, No. GT8, pp. 871-888.
Ishimato, M,. and Iida, K. (1937). "Determination of Elastic Constants of Soils by Means of
Vibration Methods," Bulletin of the Earthquake Research Institute, Tokyo Imperial University,
Vol. 15, p. 67.
Iwasaki, T., Tatsuoka, F., and Takagi, Y (1976). "Dynamic Shear Deformation Properties of Sand
for Wide Strain Range," Report of the Civil Engineering Institute, No. 1085, Ministry of Con
struction, Tokyo, Japan.
Jafari, M. K., Shafiee, A., and Razmkhah, A. (2002). "Dynamic Properties of the Fine Grained
Soils in South of Tehran," Journal of Seismology and Earthquake Engineering, Vol. 4, No. 1,
pp. 25-35.
Kiku, H., Yoshida, N., Yasuda, S., Irisawa, T., Nakazawa, H., Shimizu, Y, Ansal, A., and Erkan,
A. (2001). "In-situ Penetration Tests and Soil Profiling in Adapazari, Turkey," Proceedings,
15th International Conference on Soil Mechanics and Geotechnical Engineering, TC4 Sat
ellite Conference on "Lessons Learned from Recent Strong Earthquakes," Istanbul, Turkey,
pp. 259-269.
Kolsky, H. (1963). Stress Waves in Solids, Dover Publications, Inc., New York.
Larkin, T. J., and Taylor, P. W (1979). "Comparison of Down Hole and Laboratory Shear Wave
Velocities," Canadian Geotechnical Journal, Vol. 16, No. 1, pp. 152-162.
Matsui, T., Ohara, H., and Ito, T. (1980). "Cyclic Stress-Strain History and Shear Characteris
tics of Clay," Journal of the Geotechnical Engineering Division, ASCE, Vol. 106, No. GTlO,
pp. 1101-1120.
Mayne, P. Wand Kulhawy, F.H. (1982). " K0 -OCR Relationships in Soil," Journal of the Geotech
nical Engineering Division, ASCE, Vol. 108, No, 6, pp. 851-872.
Ohta, Y, and Goto, N. (1978). "Empirical Shear Wave Velocity Equations in Terms of Charac
teristic Soil Indexes," Earthquake Engineering and Structural Dynamics, Vol. 6, pp. 167-187.
Okamoto, T., Kokusho, T., Yoshida, Y, and Kusuonoki, K. (1989). "Comparison of Surface Ver
sus Subsurface Wave Source for P-S Logging in Sand Layer," Proceedings, 44th Annual Con
ference, Japan Society of Civil Engineers, Vol. 3, pp. 996-997 (in Japanese).
Pitilakis, K., Raptakis, D., Lontzetidis, K. T., Vassilikou, T., and Jongmans, D. (1999). "Geotech
nical and Geophysical Description of Euro-Seistests Using Field and Laboratory Tests, and
Moderate Strong Ground Motions," Journal of Earthquake Engineering, Vol. 3, pp. 381-409.
Richart, F. E., Jr., Hall, J. R . Jr., and Lysmer, J. (1962). "Study of the Propagation and Dissipa
tion of 'Elastic' Wave Energy in Granular Soils," University of Florida Report to Waterways
Experiment Station, Corps of Engineers, US Army Contract DA-22-070-Eng-314.
Rollins, K. M., Evans, M. D., Diehl, N. B., and Daily, W D., III (1998). "Shear Modulus and
Damping Relationships for Gravels," Journal of Geotechnical and Geoenvironmental Engineer
ing, ASCE, Vol. 124, No. 5, pp. 396-405.
Saxena, S. K., Avramidis, A. S., and Reddy, K. R . (1988). "Dynamic Moduli and Damping
Ratios for Cemented Sands at Low Strains," Canadian Geotechnical Journal, Vol. 25, No.2,
pp. 353-368.
Schwarz, S. D., and Musser, J. M., Jr. (1972). "Various Techniques for Making In Situ Shear Wave
Velocity Measurements: A Description and Evaluation," Proceedings, International Confer
ence on Microzonation for Safer Construction, Research and Application, Seattle, Washing
ton, Vol. 2, pp. 593-608.
Seed, H. B., and Chan, C.K (1966). "Clay Strength under Earthquake Loading Conditions," Jour
nal of the Soil Mechanics and Foundations Division, ASCE, Vol. 92, No. SM2, pp. 53-78.
Seed, H.B., and Idriss, I. M. (1970). "Soil Moduli and Damping Factors for Dynamic Response
Analysis," Report No. EERC 75-29, Earthquake Engineering Research Center, University of
California, Berkeley, California.
Seed, H. B., and Idriss, I. M. (1981). "Evaluation of Liquefaction Potential of Sand Deposits
Based on Observations and Performance in Previous Earthquakes," Preprint No. 81-544, In
Situ Testing to Evaluate Liquefaction Susceptibility, Session No. 24, ASCE Annual Confer
ence, St. Louis, Missouri.
Seed, H. B., Wong, R . T., Idriss, I. M., and Tokimatsu, K. (1986). "Moduli and Damping Factors
for Dynamic Analyses of Cohesive Soils," Journal of Geotechnical Engineering, ASCE, Vol.
112, No. GT l l, pp. 1016-1032.
Sherif, M. A., Ishibashi, I., and Gaddah, A.H. (1972). "Damping Ratio for Dry Sands," Journal of
the Geotechnical Engineering Divison, ASCE, Vol. 103, No. GT7, pp. 743-756.
References
Shibata, T., and Soelarno, D. S. (1975). "Stress-Strain Characteristic of Sands under Cyclic Load
ing," Proceedings of the Japanese Society of Civil Engineers, No. 239, pp. 57-65 (in Japanese).
Silver, M. L. (1981), "Load Deformation and Strength Behavior of Soils under Loading," State-of
the-Art-Paper, Proceedings, 2nd International Conference on Recent Advances in Geotech
nical Earthquake Engineering and Soil Dynamics, St. Louis, Missouri, Vol. 3, pp. 873-896.
Silver, M. L., and Seed, H. B.(1969). "The Behavior of Sands under Seismic Loading Conditions,"
Report No. EERC 69-16, Earthquake Engineering Research Center, University of California,
Berkeley, California.
Silver, M. L., and Seed, H. B (1971). "Deformation Characteristics of Sands under Cyclic Load
ing," Journal of the Soil Mechanics and Foundations Division, ASCE, Vol. 94, No. SM8,
pp. 1081-1098.
Stokoe, K. H. II, Wright, S. G., Bay, J. A., and Roesset, J. M. (1994). "Characterization of Geo
technical Sites by SASW Method," Technical Review: Geophysical Characterization of Sites,
ISSMFE Technical Committee 10, edited by R. D. Woods, Oxford, New Delhi, pp. 11-25
Stokoe, K. H. II, and Woods, R. D. (1972), "In Situ Shear Wave Velocity by Cross-Hole Method,"
Journal of Soil Mechanics and Foundations Division, ASCE, Vol. 98, No. SM5, pp. 443-460.
Sykora, D. E., and Stokoe, K. H. (1983). "Correlations of In-Situ Measurements in Sands of Shear
Wave Velocity," Soil Dynamics and Earthquake Engineering, Vol. 20, pp. 125-136.
Terzaghi, K. (1955). "Evaluation of Coefficients of Subgrade Reaction," Geotechnique, Vol. 5,
No. 4, pp. 297-326.
Thiers, G, R., and Seed, H. B. (1968). "Cyclic Stress-Strain Characteristics of Clay," Journal of the
Soil Mechanics and Foundations Division, ASCE, Vol. 94, No. SM6, pp. 555-569.
Vesic, A. S. (1973). "Analysis of Ultimate Loads of Shallow Foundations," Journal of the Soil
Mechanics and Foundations Division, ASCE, Vol. 99, No. SMl, pp. 45-73.
Vucetic, M., and Dobry, R. (1991). "Effect of Soil Plasticity on Cyclic Response," Journal of Geo
technical Engineering, ASCE, Vol. 117, No. 1, pp. 89-107.
Weissman, G. F., and Hart, R. R. (1961). "The Damping Capacity of Some Granular Soils,"
ASTM Special Technical Publication No. 305, Symposium on Soil Dynamics, pp. 45-54.
Werden, S. K., Drnevich, V. P., Hall, J. R. Jr., Hankour, C., Conlee, C. T., and Marr, W. A. (2013).
"New Approach to Resonant Column Testing," Geotechnical Testing Journal, Vol. 36, No. 2,
pp. 1-9.
Whitman, R. V. and Healy, K. A. (1963). "Shear Strength of Sands During Rapid Loadings,"
Transactions, ASCE, Vol. 128, Part 1, pp. 1553-1594.
Whitman, R. V. and Lawrence, F. V. (1963). "Discussion on Elastic Wave Velocities in Granu
lar Soils," Journal of the Soil Mechanics and Foundations Division, ASCE, Vol. 89, No. SM5,
pp. 112-118.
Wilson, S. D. and Dietrich, R. J. (1960). "Effect of Consolidation Pressure on Elastic and Strength
Properties of Clay," Proceedings, Research Conference on Shear Strength of Cohesive Soils,
ASCE, Boulder, Colorado, pp. 419-435.

Cap. 4

Uploaded by

Copyright:

Available Formats

Cap. 4

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Cap. 4

Uploaded by

Copyright:

Available Formats

Properties of Dynamically

Laboratory Tests and Results

lfl SHEAR STRENGTH OF SOILS UNDER RAPID

Total stress Mohr's

Figure 4.1 Unconsolidated undrained triaxial tests

Moisture content= 33.5 ± 0.2%

i= 314%/s and 426%/s

Figure 4.2 Unconsolidated-undrained triaxial test results on Buckshot clay

All data corrected to void ratio of 0.52

q, = sin- 1(�1U) - �3 J (4.4)

(obtained from static tests-

STRENGTH AND DEFORMATION CHARACTERISTICS

Figure 4.4 Transient load

Casagrande and Shannon (1949) conducted some early investigations to

where qu = unconfined compression strength. This is consistent with the

1. [o\n - <Y3]1ransient ::::::: 1.l and

,, ---'----- Shear failure at tL = 0.02s ,-.._

.....� ' r/J.

0 Ir--�..� .... ..... .., ..

Figure 4.5 Unconfined compressive strength of Cambridge clay for varying

Ill TRAVEL-TIME TEST FOR DETERMINATION OF

Figure 4.6 Unconfined compressive strength of Cambridge clay for varying

Figure 4.7 Definition of modulus of deformation, E

Figure 4.8 Confined-compression test on sand-stress- versus-strain behavior

The modulus of elasticity E can then be calculated from Eq. (3.30) as

where t' c = time of travel of longitudinal waves in a laterally confined specimen

Barium titanate ------­

--Ls�f-c::::::::-+-- Received pulse

Figure 4.9 Travel-time method: (a) schematic diagram of the laboratory

Ill BENDER ELEMENT TEST FOR DETERMINATION

G =Gmax = pvf = P( � )'

Ill RESONANT COLUMN TEST

Top view Top view

Figure 4.12 Drawings for steady-system vibration drivers in the resonant

Similarly, for torsional vibration, at resonance (with n = 1)

DE AQUÌ A. Derivation of Expressions for Ve and f for Use in the Fixed-Free-Type

Driving coil Pickup coil

4.2 mm brass frame---�

Side view �----Brass rod

Figure 4.13 Driving and measuring components for a fixed-free resonant

Figure 4.14 Derivation of Eq. (4.20)

At x = 0, U(x) = 0. So B 2 in Eq. (3.43) is zero. Thus

U(x) = B 1 sin( rot J (4.13)

where A = cross-sectional area of specimen

= ( B�:· )[ cos(rot)] CAi sinw.t + A 2 cos OJ ntl

Again, combining Eqs. (3.41), (4.13), and (4.14),

F =-m �:� =-m[B, sin(rot)]( :t 2) CAi sin . t + A2 coswnt)

AE = mOJ n v, tan( OJ;,L) (4.19)

ALp mnL mnL

IA�y =atanal (4.20)

where r =pg= unit weight of soil

Figure 4.15 Derivation of Eq. (4.23)

C. Recent Developments in Resonant Column Test Equipment

Figure 4.16 Photograph of a fully automated resonant column assembly

Figure 4.17 Schematic diagram showing quasi-fixed base and free-top

E. Shear Modulus for Large Strain Amplitudes

360 Effective confining

F. Effect on Prestraining on the Shear Modulus of Soils

To back pressure� Bottom cap

Barium titanate ------

If an SH -wave impinges the boundary between two layers, as shown in Fig