Dynamic Soil Properties
Dynamic Soil Properties
Faculty of Engineering
Department of Civil Engineering (TW15)
Laboratory of Soil Mechanics
+36 9 264 57 23
+36 9 264 57 17
Fax +36 9 264 57 49
http://terzaghi.ugent.be
Preface
I wish to thank all persons who have contributed to this thesis with their
work, their ideas and helpful comments. First of all, I would like to thank my
supervisor Prof. Wim Haegeman for providing the testing equipment and for
the extensive assistance throughout this research.
I am also very grateful to Prof. Geert Degrande from the K.U. Leuven for
his numerous valuable suggestions and support.
Furthermore I wish to thank the co-workers of the STWW-project Traffic
induced vibrations in buildings from the K.U. Leuven, especially Lincy Pyl,
Dr. Janusz Kogut, Serge Jacobs and Kathleen Geraedts.
The help of Prof. Mia Loccufier, Peter Buffel, Wouter Ost and the staff of
the Laboratory of Soil Mechanics in Ghent is very much acknowledged. Special
thanks deserves also Lou Areias for shearing his knowledge on the SCPT, for
performing of field tests and for the visual wave velocity interpretation at the
sites in Retie, Waremme and Lincent.
The research presented in this thesis is embedded in the STWW-project
IWT000152 Traffic induced vibrations in buildings. The financial support of
the Ministry of the Flemish Community for this project is likewise gratefully
acknowledged.
Lutz Karl
January 2005
Contents
Preface
xi
1 Introduction
1.1 Purpose and scope . . . . . . . . . . . . . . . . . . . . . . . . .
1.2 Outline of the thesis . . . . . . . . . . . . . . . . . . . . . . . .
3
3
4
7
7
11
19
19
21
22
23
23
23
24
28
28
29
30
31
32
33
33
34
vii
. . .
tests
. . .
. . .
. . .
. . .
. . .
. . .
. . .
. . .
. . .
. . .
. . .
. . .
. . .
. . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
viii
II
CONTENTS
35
.
.
.
.
.
.
.
.
.
.
.
37
37
38
41
41
41
44
45
47
47
48
49
53
53
53
57
60
60
61
62
62
62
63
64
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
67
67
67
71
71
72
76
77
77
79
80
81
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
85
85
85
87
87
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
CONTENTS
7.4
7.5
III
ix
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
97
.
97
.
97
.
99
.
99
.
99
. 100
. 100
. 101
88
88
90
90
91
91
92
93
94
105
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
107
107
108
108
111
117
119
121
121
125
125
128
131
142
151
152
156
160
160
160
166
169
CONTENTS
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
173
173
174
174
176
178
178
179
179
182
186
186
204
210
213
213
213
215
216
235
IV
239
Appendix
tests
. . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . .
241
242
243
244
249
G
Gmax , G0
H
h
IFFT
m, m/s, m/s2 ;
m2
m;
m, m/s, m/s2
Pa; kg/s;
m/s
, %
m
, %
Pa; J
Pa
Nm2
Hz
Hz
Pa
Pa
; m
m
xii
n , n+1
t
x
m4 ; kg/(m2 s)
, %
kg m2
Pa
kg/s2 ; 1/m;
; m/s
m2
m
kg
kg/m
m, m/s, m/s2
Pa
N/m2
s m2 ; ...;
m
m, m/s, m/s2
s
s
m
m/s
m/s
m
m
s/m;
m
1/m
, %; N/m3
m
s
m
Pa s
m
, %
kg/m3
Pa
Pa
xiii
Pa
Pa; s
; o
1/s
1/s
Subscripts
B
k
kin
max
min
p
pot
R
r
s
T
Superscripts
Other symbols
xiv
Part I
Introduction and
background
Chapter 1
Introduction
1.1
The evaluation of a number of civil engineering problems involve the transmission of waves through soil, such as seismic response under earthquake loading,
foundation response under dynamic loading and vibrations in buildings, induced by various sources. Important sources are industrial activities (looms,
printing presses), building activities (pile driving, installation of sheet piles and
demolition of buildings), rail and road traffic.
Mechanical waves can be divided in body and surface waves. Body waves
can exist in an ideal full space or they travel in a region that is not affected
by a free surface. P- (primary, compression, longitudinal) waves and s- (secondary, shear, transverse) waves are types of body waves. The particle motion
of p-waves is in propagation direction; the particle motion of s-waves is perpendicular to the direction of propagation. Surface waves may only exist at
the surface or the boundary, separating media of different properties. Rayleigh
(vertically polarized) and Love (horizontally polarized) waves are examples of
surface waves.
The prediction of elastic wave propagation from the source, through the
soil to the receiver can be performed by means of numerical models. Effects of
construction and renovation of traffic infrastructure or other vibration sources
can already be evaluated in the phase of planning. Unreasonable nuisances to
residents or damage to buildings can be avoided by design of a proper vibration
isolation system. Experimental validation has shown that the estimated soil
characteristics crucially affect the vibration prediction accuracy in the free field
and in buildings.
The most important characteristic parameters are the velocity of s- and pwaves and the material damping ratios of these body waves. Material damping
in the soil represents energy dissipation caused by friction between solid particles in the skeleton and by the relative motion between the soil skeleton and the
pore fluid. Material damping must be distinguished from geometrical or radial
3
CHAPTER 1. INTRODUCTION
1.2
Chapter 9 is dedicated to the SCPT. The selected test equipment and developed data acquisition tools are described. After a summary on methods to
determine the s-wave velocity, the chapter focuses on possibilities to measure
the damping ratio based on time recods gathered by a SCPT. The results of
an evaluation of the spectral ratio slope and the attenuation coefficient method
by means of a numerical simulation of the test are given. Finally the obtained
results on the five testing sites are summarized.
Chapter 10 deals with the methods based on the bender element test. Procedures to measure the arrival time of the s-wave are described at first. Then
methods to determine attenuation parameters are given. A resonant method is
applied on samples from two of the testing sites and the results are discussed.
Chapter 11 concludes the main findings of the thesis and gives suggestions
for further research.
The obtained insights in the accuracy of the testing techniques concerning
the G-modulus and the damping ratio are used to compare the methods under
consideration of their costs. The results of the comparison can be found in the
appendix and can be considered as a conclusion of this work.
CHAPTER 1. INTRODUCTION
Chapter 2
A harmonic loading of a soil element as depicted in Figure 2.1 leads to a stressstrain diagram as shown in Figure 2.2. It is a typical outcome of a cyclic
torsion or a cyclic triaxial test. A hysteresis loop accrues for each period of the
oscillation. The shear modulus is found as the slope of the line that connects
the point of origin and the inversion point between loading and unloading. This
modulus is also called secant modulus or equivalent modulus. It decreases with
an increase of load and shear amplitude. Therefore the stress-strain relationship
is principally non-linear. A schematic sketch of the shear modulus in function
of the shear strain is given in Figure 2.3.
The first loading curve, sometimes called the backbone curve, connects
the load inversion points of oscillation periods with different levels of shear
strain and has a hyperbolic shape. The slope in the origin point to this curve
corresponds to the dynamic shear modulus Gmax or G0 . It is also called tangent
modulus. Gmax is the shear modulus in the low strain range, usually assumed
at values below the linear elastic threshold strain of about = 5 103 %. It
is in general expected that Gmax does not change in the low-strain range.
Vibrations caused by seismic in situ tests, traffic, construction works, weak
earthquakes or even blastings usually have shear strain amplitudes below
5 103 %. This opens the opportunity to model the wave propagation with a
comparatively simple linear elastic constitutive model with Gmax as its essential
parameter.
Gmax in Pa can be calculated as:
Gmax = Vs2
(2.1)
with the soil density () in kg/m3 and the shear wave velocity (Vs ) in m/s.
This essential relation is used to obtain Gmax by means of various seismic
geophysical tests providing values of Vs .
7
Damping ratio, D
Shear modulus, G
D
Shear strain,
Figure 2.3: Shear modulus and damping ratio in function of shear strain.
Laboratory tests have shown that soil stiffness is influenced mainly by cyclic
strain amplitude, void ratio, mean effective stress, plasticity index, overconsolidation ratio and number of loading cycles. The secant shear modulus decreases
with increasing shear strain amplitude. Gmax is, according to the definition,
strain independent. The parameters which can influence Gmax are summarized
in Table 2.1 taken from Dobry & Vucetic (1987). In function of the type of soil
some of these parameters might be irrelevant.
Increasing parameter
mean effective stress, 0
void ratio, e0
geological age, tg
cementation, c
overconsolidation ratio, OCR
plasticity index, Ip
Strain rate,
Gmax
increases with 0
decreases with e0
increases with tg
increases with c
increases with OCR
increases with Ip if OCR > 1;
stays approximately constant
if OCR = 1
no effect for non-plastic soils;
increases with for the case of
plastic materials
decreases after N cycles of large
cyclic strain amplitudes,
but recovers with time in clays;
increases with N for sand
(2.2)
with the atmospheric pressure (pa 100kPa) in the same units as Gmax and
0 , an empirical exponent k related to the plasticity index Ip
k
Ip0.72
0.5
50
(2.3)
the void ratio function Fe , varying somewhat in different studies, and a stress
10
(2.973 e0 )2
1 + e0
(2.4)
this has later been improved to cover a wider range of void ratios in Hardin
(1978):
1
0.3 + 0.7e20
Fe =
(2.5)
Fe =
(2.6)
e1.3
0
The stress exponent is often taken as n=0.5 but can be calculated for individual soils from the results of laboratory tests at different effective confining
pressures.
Other empirical relationships are proposed for different soil types. Some of
them can be found in Kramer (1996) together with correlation functions based
on parameters obtained by conventional in-situ tests as CPT (cone penetration test), SPT (standard penetration test), DMT (dilatometer test) and PMT
(pressuremeter test).
The shear modulus G at higher shear strain amplitudes can be obtained if
a hyperbolic non-linear constitutive model Duncan & Chang (1970) is used. It
shows the relation between G and Gmax :
G=
Gmax
1 + /r
(2.7)
r =
max
Gmax
(2.8)
with
max represents the maximum shear stress before failure occurs. is the shear
strain related to the calculated G. Details can be found for instance in Studer
& Koller (1997).
Ishibashi & Zhang (1993) describe the modulus reduction G/Gmax with
as a function of mean effective stress and plasticity index only:
G
Gmax
(2.9)
with
m(, Ip ) m0 = 0.272 1 . . .
. . . tanh ln
0.000556
0.4
1.3
e0.0145 Ip
(2.10)
11
K(, Ip ) = 0.5 1 . . .
. . . + tanh ln
0.000102 + n(Ip )
0.492
0.0
for Ip = 0
(sandy soils)
6 1.404
3.37
10
I
for
0
< Ip 15
5 1.115
2.7
10
I
for
I
p > 70
p
(2.11)
(2.12)
2.2
(2.13)
= k u0 sin t + c u0 cos t
(2.14)
The energy dissipated during the oscillation can be obtained using the relation du = u dt as:
ED =
c u du =
c u u dt
(2.15)
12
Q(t)
Aloop, vis =E
D, vis =
0
c u u dt = u20 c
(2.16)
1 2
ku
2 0
(2.17)
whereas AAOB represents the area between the points A, O, B in Figure 2.2.
The ratio of dissipated energy and strain energy ED, vis /ES gives:
ED, vis
2 c
=
ES
k
and by using the relationsfor the natural frequency 0 =
damping cc = 2m0 = 2 k m:
ED, vis
c
= 4
= 4
ES
c c 0
(2.18)
k/m and the critical
(2.19)
c
ED, vis
Aloop, vis
=
=
cc
4 ES
4 AAOB
(2.20)
is a constant for the viscous damped system. But the dissipated energy per
loading cycle ED, vis is proportional to the loading frequency as can be seen in
Equation (2.16). Since damping in soil is in general assumed to be frequency
independent, that means ED, vis is not a function of frequency, the viscous
damped system has to be adapted to meet this requirement. The adapted
system is called a system with hysteretic or rate independent damping.
13
(2.21)
ceq
k
k
=
eq =
=
=
(2.22)
cc
2
2
m
2 mk
eq is a function of and therefore also function of the loading frequency .
The loss factor can be expressed as energy ratio. The dissipated energy
of the hysteretic system can be given as
ED, hys = u20 ceq
(2.23)
(2.24)
(2.25)
ED, hys
2 ES
(2.26)
The afterwards used material damping ratio D can be derived from the loss
factor by
D=
(2.27)
ED, hys
Aloop, hys
=
4 ES
4 AAOB
(2.28)
Experiments show that some energy is dissipated even at very low strain
levels, so the damping ratio is never zero. Above the threshold strain, the
width of the hysteresis loop increases with increasing cyclic strain amplitude,
indicating an increasing damping ratio with increasing strain amplitude.
The concept of the equivalent damping coefficient ceq works well in the
frequency domain, but if the equation of motion should be used in the time
domain it cannot be applied.
14
In this case hysteretic damping can be approximated by an equivalent hysteretic model. This is a purely viscous system with a constant c matched to
a ceq at a certain frequency. Preferably the natural frequency of the system
is used ( = 1). In this way the amplification function of the viscous damped
system shows a reasonable agreement with the hysteretic system, at least in
the important range of the natural frequency.
Isotropic visco-hysteretic elastic material model The further treatment
of the wave attenuation needs a constitutive equation for a continuum. Basis is
an isotropic visco-hysteretic elastic material model as described in Molenkamp
& Smith (1980).
The relation between stress, strain and strain rate is given by
ij = 2Geij + 3K
m ij
+ 2G e ij + 3K
m ij
(2.29)
in which G is the shear modulus, K the bulk modulus, G the viscous shear
modulus, K the viscous bulk modulus, m the isotropic strain component, eij
the deviatoric strain tensor component and ij the Kronecker delta.
The equation can be divided into the isotropic stress component m and
the deviatoric component sij
ij = sij + m ij
(2.30)
m = 3K
(2.31)
with
m
+ 3K m
and
sij = 2Geij + 2G e ij
(2.32)
For cyclic deformations the Equations (2.31) and (2.32) can be written as
m = 3K m ei(t+m )
(2.33)
(2.34)
and
(2.35)
K = K + iK
(2.36)
and
15
G
2G
(2.37)
K
2K
(2.38)
for dilatation.
Since D should be frequency independent, G and K have to be inversely
proportional to . This concept is already used in connection with the introduction of the ceq in the SDOF system. G and K can be expressed by
converting of the Equations (2.37) and (2.38) as
G =
2Ds G
K =
2Dk K
with leads introduced in the Equations (2.35) and (2.36) to the complex moduli
G = G + 2iDs G
(2.39)
K = K + 2iDk K
(2.40)
and
Vs =
(2.41)
Vs
(2.42)
16
kk
k2 k 2
(2.43)
which can be simplified for a small attenuation (k << k) and writing the attenuation coefficient as = k to
Ds =
=
k
k
(2.44)
Ds
2 f Ds
=
Vs
Vs
(2.45)
is obtained.
Other measures of energy dissipation are the quality factor Q and the specific damping capacity related to the damping ratio D and the loss factor
by:
Q=
1
,
2D
= 2 D,
= 2D
(2.46)
(2.47)
(2.48)
1
1 2
(2.49)
Vs
Vp
(2.50)
17
The increase of the damping ratio with the shear strain can be estimated
under the same assumptions used for Equation (2.8):
D
/r
=
Dmax
1 + /r
(2.51)
Damping ratio, D
decreases with 0 ;
effect decreases with increasing Ip
decreases with e0
decreases with tg
may decrease with c
not affected
decreases with Ip
increases with
stays constant or
may increase with
not significant for
moderate and N
1.3
1 + e0.0145Ip
2
. . . 1.547
G
+1
Gmax
0.586
G
Gmax
...
(2.52)
If Equation (2.9) is used to calculate the modulus reduction, the damping ratio
can be described as a function of mean effective stress, shear strain amplitude
and plasticity index.
18
Chapter 3
Methods to determine G
and D
This chapter gives an overview of the standard soil tests for the determination
of the dynamic parameters with emphasis on the resonant column test and the
free torsion pendulum test since both tests are performed in this research to
obtain reference values for Gmax and D on two of the testing sites (Waremme
and Sint-Katelijne-Waver).
It is useful to divide the testing procedures for dynamic soil parameters in
tests working under low strain conditions so the deformations can be assumed
as elastic and tests under high strain conditions with non negligible plastic
deformations. Some of the high strain tests are able to observe the dynamic
soil behavior to the range of failure.
An overview of the relevant shear strain amplitudes in different engineering
applications and test methods is given in Figure 3.1.
3.1
Laboratory tests
20
104
103
101
102
10
seismic surv ey
v ibrations caused by traffic,
blastings and weak earthquakes
machine foundations
strong earthquakes,
farfield of explosions
nearfield of explosions
In situ Tests
Laboratory Tests
ultrasonic tests
resonant column tests
(cylindrical specimen)
resonant column
tests (hollowcylinder)
cyclic torsional shear
tests (hollowcylinder)
cyclic triaxial test
shaking table
104
103
102
101
10
Figure 3.1: Overview of possible shear strain amplitudes, Studer & Koller
(1997).
21
common test, the cyclic triaxial test, at first provides values for the dynamic
elasticity modulus and only indirectly the G-modulus.
Besides the direct testing of soil specimens, also small-scale physical models
can be subjected to a cyclic loading. These tests are performed on shaking
tables or, for models whose stress dependency has to match that of the fullscale problem, more commonly in centrifuges.
3.1.1
22
Top cap
(receiv er)
Bender
elements
Sample
Bottom cap
(sender)
of transmitting and receiving both s- and p-waves. Such elements are already
commercially available.
The dynamic shear modulus Gmax and the dynamic elasticity modulus Emax
can be calculated out of the s- and p-wave velocity using Equation (2.1) and
Emax = Vp2
(1 + )(1 2)
1
(3.1)
3.1.2
The test device consists of the standard triaxial testing equipment extended
with a cyclic axial loading unit. In some cases, the cell pressure is also applied cyclically. Isotropic or anisotropic initial stress conditions are possible. A
sketch of the device in given in Figure 3.3. Bedding errors and system compliance effects generally limit the measurements to shear strains greater than
102 %, although local strain devices can produce accurate measurements at
strain levels as small as 104 %.
3.1.3
23
The cyclic simple shear test device, as shown in Figure 3.2 is most commonly
used for liquefaction testing. A short cylindrical specimen is restrained against
lateral expansion. By applying cyclic horizontal shear stresses to the top or
bottom of the specimen, the test specimen is deformed in much the same way
as an element of soil subjected to vertically propagating s-waves. Simple shear
devices that control the vertical and horizontal stresses independently are able
to impose stresses other than those corresponding to K0 conditions.
3.1.4
The cyclic torsional shear test works with a torsional loading of a cylindrical
soil specimen. The equipment looks like a conventional triaxial device except
for the added cyclic loading apparatus. Isotropic and anisotropic initial stresses
are possible. The test is most commonly used to measure stiffness and damping
characteristics over a wide range of strain levels. Torsional testing of soil specimens produce shear strains that range from zero along the axis of the specimen
to a maximum value at the outer edge. To increase the radial uniformity of
shear strains, testing devices for hollow cylinder specimen are used.
3.1.5
24
J 0, m 0
J, m
Sample
Figure 3.6: SDOF system assumed for the behavior of sample and apparatus.
J
n h
n h
=
tan
J0
Vs
Vs
(3.2)
m
n h
n h
=
tan
m0
Vl
Vl
(3.3)
J, J0 are the mass polar moments of inertia of respectively the sample and
the top cap; m, m0 are the mass of respectively the sample and the top cap; h
is the sample height and n the circular frequency of the system at resonance.
The G-modulus can be calculated out of Vs by means of Equation (2.1). The
constrained elasticity modulus follows from the longitudinal wave velocity Vl :
E = Vl2
(3.4)
Alternatively to a harmonic sinusoidal excitation, the fundamental natural frequency can also be obtained by a single pulse excitation. The frequency of the
developing free oscillation corresponds to the natural frequency of the system.
The damping ratio can be determined by the logarithmic decrement method
as shown in the description of the free torsion pendulum technique later.
The measured response curve, in the case of harmonic excitation, can be
analyzed using the half-power bandwidth method or the circle-fit method. Both
system identification techniques will be presented later on to obtain the damping ratio by means of bender elements. Finally, if the applied dynamic force
is quantitatively known by a careful calibration of the apparatus, the amplification factor at resonance can be used to obtain the damping ratio as well.
The devices used in this research are of the Drnevich type using a Hardin
oscillator to apply a torque to the top of the sample. A dynamic axial loading
is not possible in such a device. Both isotropic and anisotropic stress stages
can be imposed.
3.1.6
The free torsion pendulum test, sometimes also called Zeevaert test after its inventor, is performed on a sample from the site in Waremme and Sint-KatelijneWaver. A schematic sketch of the device at the soil mechanics laboratory at
Ghent University is given in Figure 3.9. For earlier publications on the device
can be referred to Storrer et al. (1986) and Van Impe (1977).
25
Hanger
Dial gauge
Confinement
pressure
Air
Hardin
oscillator
Oil
Sample
Cell
Water
Horizontal arm
Counter weights
Sealing bus
Confinement
pressure
Contactless
displacement sensor
Torsion shaft
Cell
Sample
Air
Water
Figure 3.9: Test set-up of the free torsion pendulum test at Ghent University.
26
0.2
n+1
0
0.2
0.4
0
Time [s]
2h JS
a4
(3.5)
27
In this equation S is the undamped circular frequency of the equipment including the soil sample, h the height of the sample, a the radius of the soil sample
and JS the polar moment of inertia of all oscillating parts of the apparatus
(Ja ) and of the soil sample (Jp ). In the framework of this research Ja has been
carefully recalculated using the mass and the dimensions of all oscillating parts,
as shown in Figure 3.11. A value of Ja = 3.26158 kgm2 has been found for the
device in Ghent. The undamped natural circular frequency S can be calculated using the damped frequency D , obtained out of the measured vibration
period, and the fraction of critical damping :
S =
D
1 2
(3.6)
n
1
log
2
n+1
(3.7)
D = = for = 1
(3.8)
2
The level of shear strain can be estimated by:
=
2 n,n+1 a
3
rh
(3.9)
n,n+1 is the mean oscillation amplitude between the two amplitudes n and
n+1 used to calculate the damping ratio, r is the distance between the measuring point on the beam and the center of rotation. The factor 2/3 indicates
that strain is calculated for a cylindrical face on 2/3 of the sample radius.
The initial shear strain applied in this test is, in dependence of the force
of the hammer blow, in general higher than 103 %. Typical values for small
hammer blows are 102 %. However, the amplitudes decrease after a couple
of free oscillations below the linear elastic limit. A MATLAB algoritmn was
developed to pick all usable peaks of the recorded damped beam oscillation,
instead of using only the first peaks. In this way the change of damping ratio
and G-modulus could be computed during the whole decay process. Because
the strain amplitude and frequency change from cycle to cycle, G-modulus
and damping ratio can be plotted versus shear strain. The resulting curves
of several hammer blows, usually 20, are averaged to improve the accuracy of
the obtained values. The G-modulus and the damping ratio could finally be
obtained in a range between 104 % and 101 % of strain. Because of the
28
Sample
JS=J p+J a
d=2a
3.2
In situ tests
The in situ tests can also be divided in tests belonging to the small-strain
range and others belonging to the high strain range. The various seismic tests,
using an artificial vibration source and vibration sensors are classified into the
first category. To the second category belong the conventional static tests like
cone and standard penetration test, pressuremeter and dilatometer tests. They
provide an indirect way to obtain dynamic moduli by means of correlation functions. Dynamically loaded steel plates open the possibility to gather stiffness
information of the first decimeters close to the surface in the low- and in the
high-strain range.
The seismic tests focus in general on the determination of the velocity of
p- and s-waves. If the generated wave is measured at several sensors also a
calculation of the damping ratio is possible. Especially for the cross-hole test
attempts are known from different authors.
3.2.1
The seismic reflection test measures wave propagation velocities and thickness
of superficial layers. The method follows the principle of echo-sounding and
radar. The test is performed by producing an impulsive disturbance at the
source, S, and measuring the arrival time at the receiver, R, located at a certain
distance from the source as shown in Figure 3.12. Dependent on the used source
29
the test can analyze p- and s-waves. However, because the generation of high
energy s-waves is difficult, the separation of the s-wave from the first arriving
p-waves might fail.
R
direct wave
H
2 ic
v1
v2
II
Figure 3.12: Ray paths in the seismic reflection test.
Some of the wave energy will follow the direct path from source to receiver.
Another part will travel downwards until it is reflected at the boundary of an
underlaying layer. The wave velocity of the superficial layer can be calculated
from the arrival time of the direct wave.
Because the angle of incidence at the layer interface has to be equal to the
reflexion angle and using the wave velocity already known, the thickness of the
superficial layer can be calculated.
If two or more receivers are used, a possible inclination of the layer interface
can be theoretically estimated. However, the practical realization fails in the
most cases. The properties of deeper layers may be evaluated using reflections
from deeper interfaces.
The method is limited to situations where the arrival times of the direct
and the reflected wave are sufficiently different. This means for instance that
the method is especially confident for deep layers and less for shallow layers. If
used, also arrival times from waves reflected at several layer interfaces have to
be distinguishable.
Additional information on the reflection test can be found in Kramer (1996).
3.2.2
30
(3.10)
At the critical angle ic , that means short before a total reflection appears,
the refracted ray travels in layer two horizontally, parallel to the boundary.
This ray will send continuously parts of its energy back to layer one. At the
passage of the layer boundary refraction reoccurs and a head wave traveling
towards the surface develops. The refraction angle in layer one is the same
as the critical angle of incidence ic . Because the wave velocity in layer two is
higher than in layer one the head wave arrives from a certain distance from the
source xc on the surface before the component taking the shorter direct path
through layer one with the lower wave speed.
The wave velocity in layer two, v2 , can be calculated based on the velocity
in the superficial layer, v1 , and the arrival time of the head wave at at least
two distances from the source or graphically from the slope in the time arrival
diagram as seen in Figure 3.14. Because the length and the shape of the travel
path of the refracted wave is known if v1 and v2 are obtained, the thickness H
of the superficial layer might be calculated.
The method is also applicable for inclined layer interfaces and multi-layered
stratifications, if the wave propagation velocity increases with the layer depth.
Closer details for the practical application can be found in Kramer (1996)
or in the geophysical literature for instance in Udas (1999).
3.2.3
31
x
t
R
head wave
direct wave
v1
v1
I
ic
II
ic
head wave
1/v2
direct wave
1/v1
v2>v 1
xc
(refracted wav e)
A cross power spectrum between the two signals is calculated. The unwrapped
phase of this spectrum is used to calculate an experimental dispersion curve of
the Rayleigh wave velocity.
R1
d1
R2
d2
3.2.4
The seismic cross-hole test measures the p- and s-wave velocities between
boreholes. At least two boreholes are necessary. The first for installing a
seismic source. This might be a mechanical or an explosive source. If the focus
is on the survey of the s-wave velocity the preference is on a mechanical source
able to produce s-wave impulses with reversible polarity.
32
In the second hole a receiver is installed at the same depth as the source in
the first borehole. The measurement is triggered at the source and the arrival
time of s- respectively p-wave is obtained by visual interpretation of the signal
from the receiver borehole. By testing at various depths, a velocity profile can
be drawn. A sketch of the test set-up is given in Figure 3.16.
Because the trigger time measurement is potentially inaccurate it is desirable to use more than two boreholes, that means more than one receiver
point. The wave velocity is then calculated from the difference in the arrival
times at the receiver holes. This has the additional advantage that instead of
a visual interpretation, cross correlation can be used. Typical distances between the boreholes are 5 to 12 m for layered soils and up to 30 m for nearly
homogeneous sites.
Mok et al. (1988) describe the application of the attenuation coefficient
method for the determination of the damping ratio based on a cross-hole test.
3.2.5
The source in the down-hole test is located at the surface close to a borehole
with an installed receiver at a certain depth. The generated waves travels nearly
vertically from source to receiver. In the up-hole test the source is situated in
the borehole and the measurement is done at the surface. The set-up is shown
in Figure 3.16.
The down-hole test is more commonly used than the up-hole test because
it is more convenient to place and adjust a seismic source at the surface than
in a borehole.
The analysis of the arrival times is done as in the cross-hole test by visual
interpretation or if more than one receiver at different distances from the source
are used also cross correlation can be applied.
R
a)
R1
R2
b)
c)
Figure 3.16: a) Seismic cross-hole test, b) seismic up-hole test, c) seismic downhole test.
3.2.6
33
The seismic cone penetration test (SCPT) can be seen as a special version of
a down-hole test with the receivers (geophones or accelerometers) installed in
the tip of a cone pushed into the ground by a conventional cone penetration
equipment (CPT truck, Figures 3.17 and 3.18). Since no borehole is necessary
the test is much less expensive than a down-hole test.
P- and s-wave sources are placed at the surface beside the penetration point
of the cone. The sources consist of steel beams or plates which are hit by a
hammer horizontally (s-wave source) or vertically (p-wave source). The cone
with the receivers is pushed stepwise into the ground. Usual intervals are 0.5
or 1.0 m. At each step the source generates a seismic pulse recorded by the
cone receivers. The determination of the p- or s-wave arrival can be performed
visually. If two receivers in a certain distance are installed in a cone (dual cone)
the travel time between these receivers can be calculated by cross correlation.
The travel time leads directly to the wave velocity using the direct wave travel
path from the source to the cone.
More details are discussed later when a method to obtain damping ratio
out of SCPT data is presented in Chapter 9.
3.2.7
Geotomography
34
3.2.8
High-strain tests
Various in-situ tests working in the high-strain range are in use. They provide
soil stiffness parameters in this strain range either directly like for instance
dilatometer test (DMT) and pressuremeter test (PMT) or by means of correlation like cone penetration test (CPT) and standard penetration test (SPT).
Furthermore also correlations to the parameters in the low strain range are
developed.
The elasticity modulus of superficial layers can also be obtained by measuring the settlement of plates loaded statically or dynamically. Dependent on
the type of dynamic loading device the high strain range can be covered as well
as the low strain range. The falling weight device is an example in the low
strain range. The water cannon developed by the ETH in Zurich (Studer &
Koller (1997)) works for instance in the high-strain range.
In the case of cohesive soils and rocks dynamic stiffness parameters are
obtained by means of a free oscillation test of a laterally free part after releasing
an applied lateral force. Dependent on the initial deformation the values are
valid for low- or high-strain conditions.
Part II
Characterization of the
testing sites
35
Chapter 4
Introduction
The site is situated on a field next to the property of the architect P. Mertens, Molsebaan 43 in Retie and around the house of the architect itself. The
following tests were performed or previous data were available on the site in
Retie:
Date
01/04/1936
01/08/1977
28/01/2000
08/11/2000
13/12/2000
10/04/2001
12/05/2003
Boring
3 SASW set-ups (two from KUL
and one from UGent, concluding
in an inversion calculation from
KUL and one from UGent)
2 SCPT with a single geophone
cone
2 SCPT with a single geophone
cone
1 SCPT with a dual accelerometer
cone
Abbr.
186.0
CPT1
CPT2
CPT3
B1
SASW KUL
SASW UGent
7.4
3.2
9.0
18.0
SCPT1
SCPT2
SCPT3
SCPT4
SCPT5
10.8
12.3
12.5
12.5
12.3
Depth
[m]
68.0
38
The position of each boring or penetration test can be found in Figure 4.1.
A detailed description of the testing procedure of the SASW tests and the
SCPTs in combination with the visual interpretation of the recorded signals
can be found in the report Areias & Haegeman (2001). This report gives all
numerical values of wave velocities from SCPT1 to SCPT4 and the results of
inversion approaches with different numbers of layers on the SASW test of K.U.
Leuven.
4.2
A drilling was performed until a depth of 18 m. The soil consists of sand and fine
sand over the whole drilling. The top tertiary layer is sand from the formation
of Mol. Data of deep drillings until 68 m in close proximity to the testing site
and even deeper borings from the region of Retie are available from the Flemish
Subsoil Data Base. They confirm that the sand reaches from the surface to a
depth of at least 186 m. The stratification found by the performed boring B1
is given in Table 4.2. The soil description is not detailed enough to distinguish
between the shallow layer of quaternary deposits and the deeper tertiary layers.
The borderline is estimated at a depth of 5 m, considering the results of the
CPT, the geological map and other drillings in the neighborhood. The profile
of a drilling in Retie, Lageweg 19 is described in Table 4.3. It includes the
geological stratification of the tertiary deposits to a depth of 190 m. It can be
used to estimate the soil structure at Molsebaan 43 for great depths.
From
To
Color
[m]
[m]
0.00
0.50 dark brown
0.50
1.50 light brown
1.50
2.50 beige
2.50
3.50 beige
3.50
4.00 beige
4.00
5.00 beige
5.00
7.00 beige, green
7.00
9.00 beige, green
9.00 14.00 beige, green
14.00 15.00 green
15.00 16.00 green
16.00 17.00 green
17.00 18.00 green
Lithographic Stratification: 0
mation of Mol.
Main Component
Admixtures
topsoil
sand
sand
fine sand
fine sand
fine sand
sand
fine sand
sand
sand
sand
silt
sand
to 5 m quaternary
silt
silt
gravel
silt
silt
silt
sand
deposits, 5 to 18 m for-
39
40
From
[m]
0.00
6.00
10.00
34.00
38.00
Color
yellow
yellow
white
green
gray, white
Main
Component
fine sand
fine sand
fine sand
fine sand
sand
Admixtures
quartz
partly lignite
quartz, partly lignite
clay, partly much glauconite
much clay, glauconite, partly
quartz
46.00
58.00 gray, white unknown
glauconite, partly clay
58.00 174.00 gray, green fine sand
glauconite
174.00 182.00 gray, green fine sand
glauconite, lime
182.00 186.00 black
fine sand
much glauconite
Lithographic Stratification: 0 to 6 m quaternary deposits, 6 to 26 m formation of Mol, 26 to 46 m formation of Kasterlee, 46 to 102 m formation
of Diest, 102 to 158 m member of Dessel (formation of Diest), 158 to 190 m
formation of Berchem.
Table 4.3: Results of boring Lageweg 19.
However due to the non plastic behavior of the material, undisturbed sampling was impossible. As a consequence laboratory tests on undisturbed material
could not be performed and the density of the material could not be obtained.
The physical characteristics of Mol sand have been widely studied through
former static and dynamic tests at Ghent University and the Flemish Geotechnical Institute. Some results are given in the following section.
Sand of Mol
The sand is geologically referred to a Tertiary-Pliocene deposit at Mol in the
north-east of Belgium. Mol sands are nearly pure quartz sands. A typical
composition is 96 % of quartz mineral and 4 % of mica and traces of other
minerals. The physical characteristics and the curve of grain size distribution
are given in Table 4.4, Table 4.5 and Figure 4.2. The data were collected by
Yoon (1991) and can be used as reference for Mol sand.
The majority of the grains of Mol sand falls in the fraction of fine sand.
The given uniformity coefficient and the degree of curvature classify the sand
as poorly graded (SP) following ASTM D-2487. The maximal and minimal
possible void ratio found by laboratory compaction experiments allows to calculate a range of the dry and saturated unit weight. The dry unit weight is
between 13.55 and 16.39 kN/m3 . The saturated value reaches from 18.25 to
20.02 kN/m3 and can be taken as an assumption of the unit weight in situ
below the groundwater table. The density above the groundwater can not be
given because of the missing value of the saturation degree. But it has to be
situated between the range for the dry and the saturated density.
V Gravel
41
II Silt
I Clay
10
Sieve residue [%]
20
30
40
50
60
70
80
90
100
10
0.2
0.06
0.1
Grain size [mm]
0.002
0.01
0.001
Grain diameter
[mm]
<0.06
0.06-0.2
0.2-2.0
>2.0
Classification
Clay to Silt
fine Sand
medium to coarse Sand
Gravel and Stones
Mass
[%]
0.6
63.4
36.0
0.0
4.3
(CPT)
Three CPTs have been performed by Geologica N.V. around the building
no. 43. A 50 kN CPT mobile apparatus has been used. The device is not
automotive. It is maneuvered by hand and needs to be anchored in the soil in
order to archive the required reaction force for a CPT sounding.
The resulting profiles of the cone resistance qc for each single CPT and an
averaged profile are given in Figure 4.3. The cone resistance of the tests clearly
indicates a weaker layer at the depth of 4 to 5 m with a thickness of 0.5 and
1 m. Such a soft layer is not seen in the profile of the drilling. Figure 4.4
shows the undrained angle of internal friction calculated from the averaged
cone resistance. The value is approximately 30 at depths below the weak
intermediate layer and 34 above.
The groundwater table measured in the holes of the CPTs is 1.15 m below
the surface.
4.4
4.4.1
(SCPT)
Test description
SCPT1 to SCPT4 were performed using a 200 kN CPT truck fitted with additional tracks. The tracks can be lowered to support the conventional wheel
drive. Areas of difficult access can be reached by this means. For SCPT5 a
200 kN automotive remote-controlled track vehicle was available that is trans-
42
Unit
mm
kN/m3
kN/m3
kN/m3
kN/m3
Value
0.195
1.6
1.02
2.65
0.918
0.585
13.553
16.387
18.25
20.02
CPT 1
0
CPT 3
CPT 13
0
Mean q
Depth [m]
10
q [MN/m]
c
20
10
q [MN/m]
c
20
10
q [MN/m]
c
20
Range of
standard
deviation
10
q [MN/m]
20
ported on the whole on top of a truck in the public road traffic. The CPT
truck and the track vehicle provide a sufficient dead load therefore anchoring
in the ground is not necessary.
Test depths went to 12.5 m and were limited by the mechanical strength of
the cone and the extension rods or by the resolution of the vibration sensors.
The maximum safe load on the cone at penetration refusal is 100 kN.
A cone equipped with a triaxial geophone and a cross-section of 10 cm2 has
been applied in SCPT1 to SCPT4. A cone with two triaxial accelerometers
43
CPT 1-3
0
Mean
Range of
standard deviation
Depth [m]
9
10
15
20
25
30
35
40
[]
44
4.4.2
The results are presented in graphical form for shear wave velocity Vs and compression wave velocity Vp in the Figures 4.5 and 4.6. Vp values are available
for SCPT1 and SCPT2 only. Both Vs and Vp generally show a trend of increasing velocity with depth. Because of less arbitrary of the cross correlation
and the smaller scatter in the profile it can be assumed that SCPT5 gives the
most reliable values. At larger depths the compression wave velocities come
close to the velocity of compression waves in the present groundwater, about
1450 m/s.
SCPT 14
SCPT 12
0
No. 1
No. 2
No. 1
No. 2
No. 3
No. 4
6
Depth [m]
Depth [m]
10
10
12
12
14
200
400
600
14
500
1000
1500
2000
SCPT1 to SCPT4 belong to the first SCPTs performed by Ghent University. The adjustment of the seismic source was not always optimal. Because
with SCPT5 a reliable profile is available a re-analysis of SCPT1 to SCPT4 with
the cross correlation technique is not considered and also not very promising.
45
SCPT 5
0
Depth [m]
8
Mechanical Hammer
Sledgehammer
10
12
50
100
150
200
250
300
350
400
450
500
4.4.3
The gathered time records of SCPT5 have been used for an analysis of the
damping ratio D of the shear wave component by means of the spectral ratio
slope method. The technique is described in Section 9.4.2. The results of the
damping ratio given in Figure 4.7 show major variation. At a few depths even
negative values are found. It is not reasonable to fit a line to obtain a closed
profile. The shallow layer until about 6 m depth shows a mean damping ratio
of 4 %, below this layer the value exceeds 10 % following the increasing trend of
the shear wave velocity. Under the assumption that the reached shear strain is
below 103 %, the values in comparison with values found in literature appear
to be quite high. Figure 4.8 gives example results from laboratory tests on
sands collected by Seed & Idriss (1970).
To provide alternative and perhaps more reliable values of the damping
ratio the correlation function of Ishibashi & Zhang (1993) is applied (Equation
(2.52)). The correlation function has been developed from a statistical analysis
of results of a large number of laboratory tests, mostly resonant column tests,
published in literature. It allows to estimate the damping ratio of a large range
of soils, from sands to moderately over-consolidated clays. The mean effective
stress
0 , the plasticity index Ip and the shear strain are necessary input
values. A profile of the damping ratio for the site of Retie is calculated with
the assumption of the soil parameters in Table 4.6 and = 103 %.
46
Depth [m]
Depth [m]
Mechanical Hammer
Sledgehammer
10
10
12
10
20
30
40
12
Figure 4.7: Values of damping ratio found from SCPT5 (left) and profile obtained by Ishibashis equation (right).
Figure 4.8: Damping ratio for sand, mean values and scattering range from
Seed & Idriss (1970).
Unit weight:
Angle of internal
friction :
K0 :
Plasticity index Ip :
47
Unit
kN/m3
kN/m3
Value
18.0
10.0
34.0
30.0
0.441
0.500
0.0
4.5
4.5.1
(SASW)
Test description
48
Depth [m]
Depth [m]
SASW UGent
0
100
200
Swave velocity [m/s]
300
100
200
Swave velocity [m/s]
300
Figure 4.9: SASW performed by Ghent University and K.U. Leuven, S-wave
velocity, Retie.
4.5.2
Test results
A model with 3 layers on a half-space has been used in the inversion calculation of K.U. Leuven and Ghent University. The density was chosen as
= 2000 kg/m3 . The coefficient of Poisson was assumed as = 0.48 from
K.U. Leuven and as = 0.30 from Ghent University. It should be noted that
the inversion calculation is very less sensitive regarding the chosen density and
the Poisson ratio. Deviating assumptions as in this case have nearly no influence on the results. A variation of the number of layers in the initial profile of
49
the inversion was done by K.U. Leuven. However the variation in the obtained
shear wave velocities was minor.
Even though the site has a homogeneous sand profile, the shear wave profiles from K.U. Leuven and Ghent University show a weaker layer with lower
propagation velocity. The profile from Ghent University has the soft layer between 1.2 and 2 m depth and the profile from K.U. Leuven at a depth between
about 0.46 and 0.77 m. The calculated velocity at depths below 3.70 m is with
250 and 253 m/s for both tests similar.
4.6
Figure 4.10 summarizes the geotechnical findings for the site in Retie. It gives
a simplified profile of the performed drilling B1 and the shear wave velocity.
The latter profile is based on the results of SCPT5 for depths below 1.50 m.
The velocity of the shallow layer is taken from both SASW tests. The soft layer
indicated by the CPTs and SASW tests is because of its uncertain depth and
thickness not reflected in the simplified shear wave velocity profile. Furthermore
it could not be confirmed by the boring and the SCPTs. The averaged CPT
profiles for tip resistance qc and internal friction are repeated. Table 4.7 gives
the numerical values.
50
Averaged CPT-Profile
Depth
qc
[m] [MN/m2 ]
[ ]
0.20
0.60
26.70
0.40
1.53
29.96
0.60
2.90
31.57
0.80
3.47
31.55
1.00
5.70
33.37
1.20
8.07
34.83
1.40
9.33
34.97
1.60
10.60
35.31
1.80
10.63
34.90
2.00
10.77
34.64
2.20
10.00
33.95
2.40
10.50
33.92
2.60
10.07
33.37
2.80
9.50
32.65
3.00
11.00
33.53
3.20
10.50
33.03
3.40
9.60
32.27
3.60
10.10
32.33
3.80
7.15
29.13
4.00
4.80
27.23
4.20
5.20
27.53
4.40
3.75
24.80
4.60
3.00
22.32
4.80
3.20
20.23
5.00
9.60
30.58
5.20
8.50
29.67
5.40
7.60
28.62
5.60
9.65
30.09
5.80
9.75
29.98
6.00
8.20
28.58
6.20
8.00
28.19
6.40
7.60
27.48
6.60
8.35
27.87
6.80
8.20
27.58
7.00
9.90
29.00
7.20
11.35
29.68
7.40
15.75
31.42
7.60
15.20
31.03
7.80
15.50
31.02
8.00
11.00
28.80
8.20
11.20
28.77
8.40
11.00
28.48
8.60
12.00
28.95
8.80
13.50
29.63
Depth
[m]
1.34
1.86
2.38
2.89
3.41
3.93
4.44
4.96
5.48
5.99
6.51
7.03
7.54
8.06
8.58
9.09
9.61
10.12
10.64
11.16
11.67
Layer
1
2
3
4
Layer
1
2
3
4
51
Profile of
Vs[MH]
[m/s]
232.53
236.93
213.81
201.65
200.16
211.03
243.12
240.19
223.52
216.02
256.70
304.54
344.96
356.10
362.70
360.40
339.19
314.46
378.21
302.68
353.76
SCPT5
Vs[SH]
[m/s]
225.53
230.33
228.68
218.30
206.83
210.14
226.73
222.77
212.77
226.58
260.62
282.66
331.56
318.50
311.04
362.95
333.69
324.38
329.72
388.34
389.94
SASW UGent
Thickness
[m]
1.2
0.3
0.5
Vs
[m/s]
220
80
120
250
SASW KUL
Thickness
[m]
0.46
0.31
2.93
Vs
[m/s]
143
46
159
253
Vs[MEAN]
[m/s]
229.03
233.63
221.24
209.97
203.49
210.58
234.92
231.48
218.15
221.30
258.66
293.60
338.26
337.30
336.87
361.68
336.44
319.42
353.96
345.51
371.85
Table 4.7: Numerical values of the averaged CPT, SCPT5 and SASW results.
52
Chapter 5
Introduction
The testing site is located on a field between the Rue de la Bruy`ere in Lincent
and the high speed railway track Brussels-Cologne. A very elaborated soil
testing campaign is performed in preparation of the construction works to the
new build high speed track in 1993 and 1994. Information on the results of in
situ and laboratory tests in close proximity to the testing site is provided by
TucRail. They delivered conventional soil parameters of the site. SASW tests
and SCPTs are performed by the Ghent University and K.U. Leuven. Table
5.1 gives an overview of all tests performed at this site.
A site plan is given in Figure 5.1. It shows the position of each in situ
test. For the testing procedure of the SASW test and the visual interpretation
of the SCPT can be referred to the report Areias & Haegeman (2001). The
report gives the numerical values of wave velocities from SCPT1 to SCPT3.
The SASW test of K.U. Leuven is described in Pyl & Degrande (2001a).
5.2
The drilling with the closest proximity to the locations of the SCPTs, SASW
tests and train vibration measurements is boring B1108. It reaches a depth of
15 m. A simplified description of the boring profile is given in Table 5.2.
A top layer with a thickness of 1.2 m is followed by a layer of fine sand
reaching to a depth of 3.2 m. Between 3.2 and 7.5 m a sequence of very stiff
layers of arenite embedded in clay has been found. Arenite is a sediment of the
residue of a sand stone transported from another place. The fixations between
the individual grains can reach from intact to very loose. None of the performed
CPTs or SCPTs could pass this layer because the maximum pushing force was
reached.
Below the arenite layers is a layer of clay from 7.5 to 8.5 m depth, followed
by fine sand from 8.5 to 10.0 m. Between 10 and 15 m thin layers of fine sand
53
54
01/02/1939
24/06/1991
18/08/1993
06/07/1993
09/04/2001
01/06/2001
3 SASW set-ups
15/05/2003
55
Abbr.
Depth
[m]
34.0
74.0
B1108
CPT1107
CPT1108
CPT1109
SCPT1
SCPT2
SCPT3
SASW KUL1
SASW KUL2
SASW UGent
SCPT4
SCPT5
15.0
4.0
5.5
7.2
5.0
5.5
7.5
5.7
7.2
56
From
[m]
0.00
1.20
2.00
2.50
To
[m]
1.20
2.00
2.50
3.20
brown
gray, brown to olive-green
beige to gray beige
gray to dark gray
Main
Component
silt
fine sand
fine sand
fine sand
3.20
7.50
a)
Sequence of:
dark gray to green
arenite
b)
clay
8.00
8.50
10.00
10.50
11.00
11.50
12.50
a)
clay
clay
fine sand
clay
fine sand
clay
sand
clay
7.50
8.00
8.50
10.00
10.50
11.00
11.50
b)
12.50
13.00
13.00
15.00
a)
Color
clay
Admixtures
gravel
clay, glauconite
clay, glauconite
silt, clay,
glauconite
fine sand,
glauconite
fine sand, silt,
glauconite
silt
silt, glauconite
glauconite
glauconite
glauconite
glauconite
clay
fine sand,
glauconite
b) dark gray to dark green
fine sand
glauconite
Lithographic Stratification: 0 to 1 m quaternary deposits, 1 to 2.5 m member of Halen/Lincent (formation of Hannut), 2.5 to 13 m member of Waterschei (formation of Hannut), 13 to 15 m member of Gelinden/Orp (formation of Heers).
Table 5.2: Results of boring B1108 at the Lincent site.
surface. Sample 2 consists of the clay material in which the hard arenite layers
are embedded. Sample 3 is finally taken from an arenite layer. The grain size
distribution curve of sample 3 shows a significant amount of gravel.
The density is obtained for the clay sample only. The wet unit weight is
15.9 kN/m3 and the dry unit weight is 11.1 kN/m3 .
57
From
To
Color
Main
Admixtures
[m]
[m]
Component
0.00
1.00 brown
back filling
1.00 13.00 tawniness
clay
stones
13.00 45.00 gray
clay
petrified clay layers
45.00 50.00 light gray
silt
50.00 58.00 gray, green sand
stones
58.00 60.00 gray, green sand
glauconite
60.00 74.00 white
porous chalk Lithographic Stratification: 0 to 1 m quaternary deposits, 1 to 45 m formation of Hannut and formation of Tienen, 45 to 50 m member of Gelinden
(formation of Heers), member of Orp (formation of Heers), 60 to 74 m
cretaceous formations.
Table 5.3: Results of the boring in Waasmont.
From
To
Color Main
Admixtures
[m]
[m]
Component
0.00
4.00 brown silt
glauconite, clay, sand
4.00
9.00 gray
clay
sand, glauconite
9.00 17.00 gray
silt
17.00 21.00 gray
silt
sand, mica, sporadic glauconite
21.00 31.00 gray
silt
much sand
31.00 34.00 white
chalk
Lithographic Stratification: 0 to 4 m quaternary deposits, 4 to 9 m member
of Lincent (formation of Hannut), 9 to 31 m member of Gelinden (formation
of Heers), 31 to 34 m chalk (cretaceous formations).
Table 5.4: Results of the boring in Lincent, Rue de Li`ege.
Groundwater table
The borehole of drilling B1108 was converted to a monitoring station for the
groundwater table. Form August 1993 to July 1997 the water level is measured.
It changes in a range between 6.0 m to 12.2 m below the surface, the mean
value is 10.4 m.
5.3
(CPT)
58
Depth
Grain size distribution
Frac- Grain
tion
diameter
I
<0.002 mm
(clay)
II
0.002-0.06 mm
(silt)
III
0.06-0.2 mm
(fine sand)
IV
0.2-2.0 mm
(medium/
coarse sand)
V
2.0-20.0 mm
(fine gravel)
VI
>20.0 mm
(coarse gravel/
stones)
Specific gravity, GS
Water content, w
Void ratio, e
Degree of Saturation, Sr
Unit weight,
Dry unit weight, d
Saturated unit weight, r
Liquid limit, wL
Plastic limit, wP
Plasticity index, IP
Classification
Unit
m
Sample 1
1.0-1.4
Sample 2
3.0-3.4
Sample 3
5.0-5.4
16
34
18
57
26
15
24
32
15
13
32
%
%
%
kN/m3
kN/m3
kN/m3
%
%
%
0
2.633
12.87
0
2.594
42.99
1.2864
86.7
15.9
11.1
16.7
71.1
39.8
31.3
OH
organic
clay
7
2.514
24.50
31.6
21.7
9.9
CL
silty
clay
59.0
38.3
20.7
GC
gravel-sandclay mixture
VI
V fine Gravel
59
II Silt
I Clay
20
30
40
50
60
70
80
Depth 1 (1.01.4 m)
10
Depth 2 (3.03.4 m)
Depth 3 (5.05.4 m)
90
20
10
0.2
0.1
Grain size [mm]
0.01
0.06
0.002
CPT 1107
CPT 1109
CPT 1108
CPT 11071109
0
Mean q
Range of
standard deviation
Depth [m]
100
50
q [MN/m]
c
100
50
q [MN/m]
c
100
50
100
q [MN/m]
c
50
100
q [MN/m]
c
0.001
60
Depth [m]
10
15
20
25
30
35
40
45
[]
5.4
5.4.1
(SCPT)
SCPT1 to SCPT3 are performed using the 200 kN CPT truck fitted with
additional tracks. SCPT4 and SCPT5 are done by the 200 kN automotive
remote-controlled track vehicle. The tests penetrate to depths between 5.0 and
7.2 m and are limited again by the mechanical strength of the cone and the
extension rods as soon as the hard arenite layer is reached.
The seismic single cone with the triaxial geophone array is used in SCPT1
to SCPT3. The dual cone with two triaxial accelerometers in 1 m distance
is used for SCPT4 and SCPT5. The de-coupled seismic source provides the
61
needed shear wave excitation in most cases. Only for SCPT2 the coupled
seismic source is in use. Some closer details on cone and source can be found
in Section 4.4.1.
SCPT 13
SCPT 45
0
No. 1
No. 2
No. 3
3
Depth [m]
Depth [m]
No. 4
No. 5
100
200
300
400
Swave velocity [m/s]
500
100
200
300
Swave velocity [m/s]
400
5.4.2
Figure 5.5 presents the profiles of shear wave velocity from all performed
SCPTs. The plot on the left side shows the results of SCPT1 to SCPT3,
analysed using the visual interpretation of the first shear wave arrival. The
right side gives the profiles of SCPT4 and SCPT5. These are calculated by
means of cross correlation between the two simultaneously recorded signals of
the dual cone in one of the horizontal axes. Each of the lines in the right figure
is the average of the two independently obtained velocity profiles using a sledgehammer on one side of the beam and the mechanical swing hammer on the
other side. The results of the shear wave velocity show a very good agreement
with both types of hammers.
As already seen on the site of Retie, the profiles obtained by visual interpretation show a much higher scatter and variation than those obtained by
the cross correlation technique. Therefore also in this case it is assumed that
SCPT4 and SCPT5 deliver more reliable values on the soil profile than the
other SCPTs.
The curves of SCPT4 and SCPT5 show a quite linear increase of the shear
wave velocity from about 160 m/s at 1 m depth to about 310 m/s at 6.5 m
62
5.4.3
Depth [m]
5
6
Damping Ratio [%]
10
5.5
5.5.1
(SASW)
63
the excitation weight is placed close to the rail track as shown in the site plan.
In the first setup the accelerometers are installed at 3, 4, 6, 8, 12, 16, 24, 32
and 48 m distance from the source; setup two uses accelerometers at distances
of 8, 12, 16, 20, 24, 32, 40, 48, 64, 80, 96, 128, 160 and 198 m.
Detailed information on the SASW tests of K.U. Leuven are given in Pyl &
Degrande (2001a). Inversion approaches with different numbers of layers are
commented in this report. The final inversion calculation for the first setup uses
a model with 4 layers on a halfspace; for setup two three layers on a halfspace
are used. The coefficient of Poisson = 0.33 and the density = 1800 kg/m 3
are kept constant in the inversion procedure.
SASW performed by Ghent University - SASW UGent
A SASW test is performed in the direct neighbourhood to the SASW measurement line of K.U. Leuven. The drop weight of the Ghent University and a
hand-held hammer is used to generate the surface waves. The testing procedure
is the same as at the Retie site.
The inversion calculation shows a profile of two layers on the halfspace.
The density and Poissons ratio are chosen = 2000 kg/m3 and = 0.3, respectively.
5.5.2
Test results
SASW KUL1/KUL2
SASW UGent
0
3
Depth [m]
Depth [m]
Setup 1 (KUL1)
Setup 2 (KUL2)
100
200
300
Swave velocity [m/s]
400
100
200
300
Swave velocity [m/s]
400
Figure 5.7: SASW performed by Ghent University and K.U. Leuven, S-wave
velocity, Lincent.
64
The shear wave velocity profiles are given in Figure 5.7. The left side shows
the profile obtained by the Ghent University and the right side the two profiles
of K.U. Leuven.
The results of all inversion calculations show an increasing stiffness with
depth. At 3.0 respectively 3.5 m depth a stiffness increase is observed. This is
probably caused by the arenite layer.
5.6
Figure 5.8 summarizes the gathered geotechnical data for the site in Lincent.
It gives a simplified profile of drilling B1108. The layers are labeled based on
there main components. Names in brackets give information on the components
besides this main component. They are partly questionable. Which means a
clayey sand could also be a silty sand, a sandy clay, a silty clay.
The profile of the shear wave velocity is summarized out of the results of
SCPT4 and SCPT5 for depths below 1 m. It is assumed that both tests deliver
the most reliable values. They show a very uniform increase of the velocity
with depth. Furthermore the SASW profiles plot in the same range. The swave velocity of the very shallow layer is taken as a mean value from all three
performed SASW tests.
Table 5.6 presents selected numerical values for CPT, SCPT and SASW
tests.
Mean
Groundwater
10.40 m
1.00
(Silty Clay)
15.00
13.00
(Silty Clay)
Sequence of
Fine Sand / Clay
10.00
Fine Sand
8.50
Clay
7.50
(Sandy Clay
to Silty Clay)
Sequence of
Arenite / Clay
3.20
(Clayey Sand)
Fine Sand
1.20
Silt
0.00 m
Quaternary
Deposits
Formation of Hannut
(Tertiary Deposit)
Formation of Heers
(Tertiary Deposit)
7
0
20 40
phi []
310
160
140
Vs [m/s]
0 200 400
S-wave velocity [m/s]
Simplified S-Wave
Velocity Profile
Site Lincent
0 20 40 60
qc [MN/m]
Internal
0 Friction
Tip
0 Resistance
Depth [m]
Depth [m]
Averaged CPT-Profiles
Depth [m]
Simplified Stratification
of Drilling B1108
7.20
1.70
0.00 m
66
Averaged CPT-Profile
Depth
qc
[m] [MN/m2 ]
[ ]
0.20
12.20
40.00
0.40
9.75
40.00
0.60
9.37
37.80
0.80
6.50
34.10
1.00
9.03
34.67
1.20
10.90
34.97
1.40
8.33
34.40
1.60
4.17
29.47
1.80
1.37
20.80
2.00
2.53
23.53
2.20
13.53
30.93
2.40
14.37
32.37
2.60
1.30
17.13
2.80
2.33
20.87
3.00
3.83
24.10
3.20
40.77
36.07
3.40
47.07
37.50
3.60
10.10
28.20
3.80
6.27
23.70
4.00
33.30
33.53
4.20
48.00
36.35
4.40
43.00
35.85
4.60
19.60
28.45
4.80
8.50
28.10
5.00
42.80
35.00
5.20
58.00
37.55
5.40
15.10
31.05
5.60
11.00
28.95
5.80
8.15
26.15
6.00
23.95
29.20
6.20
44.15
32.50
6.40
52.50
35.75
6.60
60.00
36.80
6.80
44.00
35.30
Profiles of SCPTs
SCPT4
SCPT5
Depth
Vs
Vs
[m]
[m/s]
[m/s]
0.83
421.31
173.21
1.34
174.63
166.65
1.86
166.23
190.10
2.38
176.08
203.97
2.89
216.80
201.14
3.41
223.96
221.10
3.93
231.00
251.71
4.44
231.92
274.00
4.96
253.00
279.62
5.48
272.43
5.99
280.41
6.51
324.83
SASW UGent
Thickness
[m]
1
1.0
2
2.5
3
Vs
[m/s]
110
165
300
SASW KUL1
Thickness
[m]
0.19
2.53
0.22
3.33
Vs
[m/s]
108
158
176
248
249
SASW KUL2
Thickness
[m]
3.02
0.10
0.56
Vs
[m/s]
150
280
278
282
Layer
Layer
1
2
3
4
5
Layer
1
2
3
4
Table 5.6: Numerical values of the averaged CPT, SCPT4, SCPT5 and SASW
results.
Chapter 6
Introduction
The site is situated in the garden of the house Rue de la Costale 55 in Waremme.
The property borders directly on the high speed railway track Brussels-Cologne.
Data from the elaborated soil testing campaign performed in preparation of
the construction works to the new build high speed track in 1993 and 1994
are provided by TucRail. They delivered conventional soil parameters of the
site. SASW tests, SCPTs and a boring are performed by the Ghent University
and K.U. Leuven. Undisturbed and also disturbed material obtained by the
boring are tested by means of bender element, resonant column and free torsion
pendulum test to determine Gmax and D. Table 6.1 gives an overview of all
considered tests performed at this site. A site plan is given in Figure 6.1.
For the testing procedure of the SASW test and the visual interpretation
of the SCPT can be referred to the report Areias & Haegeman (2001). The
report gives the numerical values of wave velocities from SCPT1 and SCPT2.
The SASW test of K.U. Leuven is described in Pyl & Degrande (2001b).
6.2
The boring B1 is driven on the meadow on a direct line between the house
Rue de la Costale 55 and the train track. The description of the boring profile
is given in Table 6.2. It shows beneath a shallow top layer of 1 m thickness
consisting of topsoil mixed with gravel a quite homogeneous layer of a clayey
silt. The color of the silt layer changes from rust brown over gray to black. It
is resting on a cretaceous formation of chalk. The transition zone consists of a
mixture of silt with field-stones and starts at about 12 m depth.
By means of the boring B3064 driven in 1993 on the other side of the train
track and a boring, not closer described, available from the Flemish Subsoil
Database, the first 7.5 m from the surface is assigned geologically to a quaternary deposit. The rest of the silt layer until the beginning of the chalk is
67
68
01/08/1936
20/09/1993
27/02/2001
18/08/1993
19/08/1993
31/05/2001
01/06/2001
06/04/2001
13/08/2002
13/08/2002
18/09/2002
69
Abbr.
Depth
[m]
52.0
B3064
20.0
B1
12.5
CPT3063
12.0
CPT3064
12.2
SASW KUL1
SASW KUL2
SASW UGent1
SASW UGent2
SCPT1
SCPT2
SCPT3
12.5
12.0
11.4
SCPT4
11.4
SCPT5
11.4
70
From
To
Color
Main Component Admixtures
[m]
[m]
0.00
1.00 dark brown gravel
1.00
7.50 rust brown silt
clay
7.50
8.25 gray brown silt
clay
8.25 11.50 gray
silt
clay
11.50 12.00 gray
silt
field-stone
12.00 12.50 black
silt
field-stone
12.50 12.75 black
field-stone
chalk
Lithographic Stratification: 0 to 7.5 m quaternary deposits, 7.5 to 12.5 m
member of Orb (formation of Heers), 12.5 to 12.75 m chalk (cretaceous
formation).
Table 6.2: Results of boring B1 at the Waremme site.
the dry unit weight is around 16.71 kN/m3 . More details and the plasticity
parameters are given in Table 6.5.
0
VI
V fine Gravel
II Silt
I Clay
20
30
40
50
60
70
80
10
Depth 1 (3.03.4 m)
Depth 2 (5.05.4 m)
90
100
20
10
0.2
0.1
Grain size [mm]
0.06
0.01
0.002
0.001
To
[m]
1.40
Color
1.40
3.40
4.50
7.40
3.40
4.50
7.40
8.00
8.00
9.40
9.40
11.40
11.40
12.50
brown ocher
brown, yellow ocher
yellow ocher
olive-green to
yellow ocher
dark green, grayly,
partly yellow ocher
dark green, grayly,
rust brown
dark gray, greenish
12.50
14.00
dark green
14.00
15.50
light gray
brown
71
Main
Component
silt
Admixtures
silt
silt
silt
silt
silt
fine sand
silt
silt
sand, small
angular flint
stones
fine sand
angular flint
stones (1-3 cm)
chalk
sand, small
angular flint
stones (2-5 cm)
15.50 16.00 whitish beige
chalk
sand
16.00 19.50 light gray, beige
chalk
sand, small
angular flint
stones (1-5 cm)
19.50 20.00 whitish beige
chalk
sand
Lithographic Stratification: 0 to 8.0 m quaternary deposits, 8.0 to 12.5 m
formation of Heers, 12.5 to 20.00 m chalk (cretaceous formation).
Table 6.3: Results of boring B3064 at the Waremme site.
6.3
6.3.1
In situ tests
Cone penetration test (CPT)
Two CPT profiles provided by TucRail are taken into consideration. CPT3064
was performed at the same place as B3064. CPT3063 was situated at about
85 m distance from the test site in the direction of Brussels. The profiles of the
cone resistance qc are given in the Figures 6.3 and 6.4. The average of both
is shown in Figure 6.5. CPT3063 is shifted about 1.5 m in depth to take the
different levels of the surface between both CPTs into account so that they
can be directly compared.
The penetration depth of both tests is limited by the stiff transition layer
between the clayey silt and the chalk bedrock at 12 respectively 13 m.
72
Figure 6.6 gives the profile of the correlated undrained angle of internal
friction .
6.3.2
Depth
Grain size distribution
Fraction Grain diameter
I
<0.002 mm
(clay)
II
0.002-0.06 mm
(silt)
III
0.06-0.2 mm
(fine sand)
IV
0.2-2.0 mm
(medium/coarse sand)
V
2.0-20.0 mm
(fine gravel)
VI
>20.0 mm
(coarse gravel/stones)
Specific gravity, GS
Water content, w
Void ratio, e
Degree of Saturation, Sr
Unit weight,
Dry unit weight, d
Saturated unit weight, r
Liquid limit, wL
Plastic limit, wP
Plasticity index, IP
Classification
73
Unit
m
Sample 1
3.0-3.4
Sample 2
5.0-5.4
20.1
12.3
69.0
76.4
10.6
9.3
0.3
1.1
0.9
%
%
%
kN/m3
kN/m3
kN/m3
%
%
%
0
2.641
17.9
0.5764
83.6
19.74
16.74
20.34
34.6
19.9
14.7
CL
silty
clay
31.9
3
0
2.644
22.1
0.5956
100.0
20.24
16.58
20.23
29.5
23.3
6.2
CL/ML
silty clay to
clayey silt
39.2
5
kN/m2
Depth [m]
Depth [m]
Depth [m]
74
10
10
10
12
12
12
14
10
qc [MN/m]
20
14
10
qc [MN/m]
20
14
10
qc [MN/m]
20
75
Depth [m]
10
12
14
10
20
[]
30
40
50
SCPT 45
0
Depth [m]
Depth [m]
SCPT 12
0
10
10
12
12
No. 4
No. 5
No. 1
No. 2
14
200
400
Swave velocity [m/s]
600
14
200
400
Swave velocity [m/s]
600
76
Depth [m]
10
12
Mechanical Hammer
Sledgehammer
14
2.5
5
7.5
Damping Ratio [%]
10
6.3.3
77
6.4
6.4.1
Two samples of boring B1 were installed in the triaxial bender element test
device in a way that the s-wave velocity in the vertical direction can be determined. Two other were installed turned by 90o . By this means the horizontal
s-wave velocity can be measured. The sample was consolidated under isotropic
stress conditions. The s-wave velocity is measured at the end of the consolidation by a visual interpretation of the arrival signals. A back pressure of 100 kPa
was used in al tests.
Table 6.6 summarizes the results and gives also calculated values of the
maximum shear modulus Gmax . The necessary density are actual values at the
end of the consolidation step. Figure 6.12 and Figure 6.13 show the results in
graphical form.
It is obvious that the differences in the results of horizontal and vertical
propagation are small if samples from a similar depth are compared. This
means that no significant structural anisotropy is found.
The development of Gmax during consolidation and further details can be
found in the technical report Haegeman (2002).
78
Setup 2
0
3
Depth [m]
Depth [m]
Setup 1
0
UGent
KUL
UGent
KUL
8
100
200
Swave velocity [m/s]
300
horiz.
vert.
Sample /
Depth
B1/2
2.0-2.5 m
B1/10
10.0-10.5 m
B1/3
3.0-3.5 m
B1/11
11.0-11.5 m
Vs [m/s]
Gmax [MPa]
Vs [m/s]
Gmax [MPa]
Vs [m/s]
Gmax [MPa]
Vs [m/s]
Gmax [MPa]
25
120.1
26.7
110.4
27.3
-
100
200
Swave velocity [m/s]
300
50
143.2
40.9
168.2
57.1
138.9
43.8
173.1
55.5
1,2,3
100
182.5
66.9
198.2
79.1
175.1
70.6
211.7
82.7
[kPa]
200
221.6
99.4
239.8
116.6
217.7
111.0
245.5
111.1
400
279.4
159.6
294.3
177.0
293.1
159.1
600
315.6
204.8
-
Table 6.6: Results of isotropic bender element tests with vertical and horizontal
propagation direction.
79
250
350
200
300
Gmax [MPa]
Vs [m/s]
250
200
150
50
100
200
300
400
c [kPa]
500
B 1/2, vertical
B 1/10, vertical
B 1/3, horizontal
B 1/11, horizontal
50
600
0
0
700
6.4.2
100
B 1/2, vertical
B 1/10, vertical
B 1/3, horizontal
B 1/11, horizontal
100
0
0
150
100
200
300
400
c [kPa]
500
600
700
A sample from the boring B1 from a depth of 12.0 to 12.5 m was tested in
the free torsion pendulum device of the Ghent University. The sample consists
mainly of silt and is taken from the transition area between the silt layer and
the chalk bedrock.
The dimensions of the sample were: 20 cm height, 10 cm diameter. The
isotropic stress stages have been chosen to 50, 100, 200 and 400 kPa. The
modulus at the reference shear strain of 105 is 35.1, 62.1, 100.6 and 152.9 MPa
at 50, 100, 200 and 400 kPa confining stress respectively.
The obtained results of the 200 kPa consolidation stage show an exceptional
behavior. The G-modulus increases extremely from = 105 in the direction of
the lower shear strain. The damping ratio shows a peak at = 3 106 . Even
though the test was repeated several times with the same exceptional result,
both phenomena are considered to be caused by the testing device. Therefore
the given values for the 200 kPa stage should be trusted only above a shear
strain of 105 for the G-modulus and above 104 for the damping ratio.
The damping ratio at a reference shear strain was found in the range of 1.0
to 3.0 %, neglecting the 200 kPa stage. The smaller values correspond to the
higher confining stress stages and the larger ones to smaller stresses.
Shear modulus and damping ratio in function of the applied shear strain can
be found in Figure 6.14. The same parameters in dependence of the confining
stress at a shear strain of 105 are given in Figure 6.15.
80
G-Modulus
180
50 kPa
100 kPa
200 kPa
400 kPa
160
G [MPa]
140
120
100
80
60
40
20
0 -6
10
10
-5
10
-4
10
-3
Damping Ratio
0.15
50 kPa
100 kPa
200 kPa
400 kPa
D [-]
0.1
0.05
0 -6
10
10
-5
10
-4
10
-3
Figure 6.14: G-modulus and damping ratio versus applied shear strain obtained
by the Free Torsion Pendulum Test, sampling depth: 12.0 to 12.5 m, boring
B1, Waremme.
6.4.3
81
Dampig Ratio
160
0.05
0.045
140
0.04
120
0.035
0.03
D []
G [MPa]
100
80
0.025
0.02
60
0.015
40
0.01
20
0
0.005
0
100
200
300
[kPa]
400
100
200
300
[kPa]
400
Figure 6.15: G-Modulus and damping ratio versus isotropic confining stress,
shear strain 105 , boring B1, Waremme.
6.5
Figure 6.16 summarizes the gathered geotechnical data for the site in Waremme.
It gives a simplified profile of drilling B1. The layers are labeled based on their
main components.
The profile of the shear wave velocity is concluded out of the results of
SCPT4 and SCPT5 for depths deeper than 1 m. It is assumed that both tests
deliver the most reliable values. They show a uniform increase of the velocity
with depth, which is considered by the assumption of a linear increase from
115 m/s at 1 m depth to 265 m/s at 11 m depth. Furthermore the results
obtained by the SASW tests fall in the same range. The s-wave velocity of
the very shallow layer is taken from the SASW tests on setup 1. Its value is
155 m/s.
Table 6.7 presents selected values for qc and of the averaged CPT-profile;
Table 6.8 gives the s-wave velocities determined by SCPT4, SCPT5 and the
four SASW tests.
Groundwater
6.72 m
Chalk
Site Waremme
12.75
12.50
Silt/Stone Mixture
11.50
(Clayey Silt
to Silty Clay)
Silt
1.00
Gravel
0.00 m
Simplified Stratification
of Drilling B1
12.50
7.50
Quaternary
Deposits
Creataceous Formation of Heers
(Tertiary Deposit)
formation
Depth [m]
14
12
10
10
15
q c [MN/m]
20
14
12
[]
20
40
14
12
10
10
200
400
S-wave velocity [m/s]
Simplified S-Wave
Velocity Profile
Internal
Friction
Depth [m]
Depth [m]
Averaged CPT-Profiles
Tip
Resistance
265
115
155
Vs [m/s]
0.00 m
11.00
1.00
82
CHAPTER 6. TEST SITE WAREMME
Depth
[m]
0.20
0.40
0.60
0.80
1.00
1.20
1.40
1.60
1.80
2.00
2.20
2.40
2.60
2.80
3.00
3.20
3.40
3.60
3.80
4.00
4.20
4.40
4.60
4.80
5.00
5.20
5.40
5.60
5.80
6.00
6.20
Averaged CPT-Profile
qc
Depth
qc
[MN/m2 ]
[ ]
[m] [MN/m2 ]
4.69 41.10
6.40
1.68
5.47 38.70
6.60
1.49
5.08 36.50
6.80
1.49
2.77 32.20
7.00
1.30
3.54 32.30
7.20
1.49
2.77 30.00
7.40
1.50
3.17 29.90
7.60
1.88
2.00 31.65
7.80
1.89
2.59 30.85
8.00
2.08
2.97 30.05
8.20
4.20
2.40 27.05
8.40
4.97
1.83 24.20
8.60
5.36
1.25 20.05
8.80
6.13
0.87 14.70
9.00
5.36
1.25 19.25
9.20
6.33
1.25 18.10
9.40
5.95
1.65 20.30
9.60
4.41
1.26 16.75
9.80
4.22
1.46 17.15
10.00
5.18
1.27 16.15
10.20
4.41
1.27 15.55
10.40
6.54
1.47 16.35
10.60
5.96
1.47 15.95
10.80
7.31
1.28 13.55
11.00
7.31
1.47 14.20
11.20
6.93
1.66 15.95
11.40
7.32
1.67 16.05
11.60
5.97
2.05 17.60
11.80
4.24
1.87 15.95
12.00
28.68
1.87 14.30
12.20
30.99
1.87 16.05
83
[ ]
16.50
14.35
11.25
12.45
12.10
11.90
13.60
13.95
14.55
20.95
22.30
22.40
23.85
22.50
23.45
22.70
20.35
19.85
21.35
19.90
22.60
21.05
20.80
20.65
21.80
22.30
21.45
18.20
25.65
26.75
84
Profiles of SCPTs
SCPT4 SCPT5
Depth
Vs
Vs
[m]
[m/s]
[m/s]
0.7
52.55
138.88
1.20
65.75
193.76
1.70
92.25
263.57
2.20
116.95
200.76
2.70
117.45
154.75
3.20
135.10
190.88
3.70
145.55
170.25
4.20
141.55
171.41
4.70
145.15
176.11
5.20
158.95
183.72
5.70
178.45
189.64
6.20
186.55
198.07
6.70
191.75
199.98
7.20
193.60
199.84
7.70
204.80
201.18
8.20
219.25
212.40
8.70
227.20
219.79
9.20
251.90
244.33
9.70
273.95
259.74
10.20
264.20
267.79
10.70
264.90
252.49
216
SASW KUL Setup 2
Layer Thickness
Vs
[m] [m/s]
1
0.27
145
2
0.70
179
3
179
Chapter 7
Test site
Sint-Katelijne-Waver
7.1
Introduction
The test location (I) is situated on a field on the property of the De Nayer
Institute in Sint-Katlijne-Waver. Two borings were driven and two SCPTs
were performed. In the framework of a research project on Soil Displacement
Screw Piles a previous elaborate soil testing campaign has been performed on
a research site at a distance of about 150 m. The surface of this site, later
on labelled as II, is about 2.3 m below the surface of site I. The tests have
been performed at the end of the 1990th. They deliver beside conventional soil
parameters also information on shear and compression wave velocity obtained
by SASW and seismic refraction tests. A detailed summary on this results can
be found in Menge (2001), but the most essential findings are also repeated in
the following. Table 7.1 gives an overview of the in situ tests which results are
needed later.
A sketch of the site is given in Figure 7.1. It shows the position of the in
situ test on location I.
7.2
The stratification of the site is quite simple. It consists of a disturbed top layer
of clayed sand with stones and gravel, followed by a homogeneous layer of clay
which extends to great depths. The interface between the sand and the clay is
for testing site I at about 3 m depth and for site II at 0.7 m depth.
The clay belongs geologically to the tertiary formation of Boom, the sand
of the top layer is a quaternary deposit.
85
86
Date
Site
Abbr.
14/05/2003
B1I
Depth
[m]
11.0
I
I
I
II
B2I
SCPT1
SCPT2
B1II
11.0
12.9
12.9
14.0
II
II
II
II
CPT-E
DMT
SASW1-5
SRT
14.0
13.0
14/05/2003
13/05/2003
14/05/2003
3.60
Building of the
De Nayer Institute
Entrance
14.00
12.10
7.20
3.10
Parking Space
5.00
SCPT 2
5.20
4.00
Access Road
Tree
Tree
B2
SCPT 1
Tree
Tree
B1
1.60 1.80
[m]
Groundwater table
The groundwater table is found at site I at depths between 0.5 and 1 m depth.
If a horizontal groundwater table between site I and II is assumed the table at
site II can be estimated 2.8 to 3.3 m deep.
87
Unit
kN/m3
kN/m3
%
%
%
%
%
a
15.7
19.7
25.3
(101.8)
72.2
25.4
46.8
Sample
b
c
15.3
15.5
19.3
19.5
26.0
25.8
99.2 (100.9)
65.1
75.9
25.3
26.2
39.9
49.7
d
15.4
19.4
26.1
(101.0)
71.8
26.0
45.8
Mean
Value
15.5
19.5
25.8
100.0
71.6
25.7
45.6
VI
V fine Gravel
II Silt
I Clay
20
30
40
50
60
70
80
10
90
100
20
10
0.2
0.1
Grain size [mm]
0.06
0.01
0.002
0.001
7.3
7.3.1
In situ tests
Marchetti dilatometer test (DMT)
Altogether four DMTs were performed on site II. A profile of the overconsolidation ratio OCR and another for the K0 -value were correlated using the
Depth [m]
Depth [m]
88
10
10
12
12
14
20
40
60
Overconsolidation ratio
14
80
Figure 7.3: Overconsolidation ratio derived from dilatometer test results at site II.
horizontal stress index KD , obtained directly from the test. Both are shown
in Figure 7.3 and Figure 7.4. Menge (2001) gives details on the correlation
equation.
The OCR-value in the boom clay is in the range between 70 at shallow
depths and 15 at 14 m depth. K0 decreases in the same zone from 3 to 1.8.
Because the top layer at site I is 2.3 m thicker than at site II it is not possible
to convert the values directly to site I. However, it can be assumed that site I
would show similar results.
7.3.2
CPTs have been done during the research project on Soil Displacement Screw
Piles using electrical CPT-cones. The averaged profile of the tip resistance q c
of 6 of the CPTs is drawn in Figure 7.5. A high reliability of the CPT results
is concluded from the small range of standard deviation. Some discontinuities
in the averaged profile can be explained by the presence of a few boulders at
several depths. Also in this case the different thick top layers at site I and II
have to be kept in mind if the profile should be used also for site I.
7.3.3
Test setup
The two SCPTs on site I are performed using the 200 kN automotive remotecontrolled track vehicle. The cone was driven to a final depth of about 13 m in
89
Depth [m]
10
12
14
2
3
q [MN/m]
Figure 7.5: Averaged cone resistance qc and range of standard deviation, site II.
SCPT 12
SCPT 2
0
Mechanical Hammer
Sledgehammer
Depth [m]
Depth [m]
No. 1
No. 2
100
200
Swave velocity [m/s]
300
4
6
8
Damping Ratio [%]
10
90
7.3.4
5 measurement lines were tested by the SASW method on site II. The resulting
profiles are presented in Figure 7.8, the numerical values are given in Table
7.3. The inversion calculation was performed assuming three layers on the
halfspace, for SASW 2 and 5 four layers of the halfspace were preferred because
the measured dispersion curve could be better approximated and because a fill
material was used at both measurement lines to replace the first 30 cm of the
top soil material. The very similar profiles from the five SASW tests show an
increasing stiffness with depth. Within the clay layer the shear wave velocity
increases from about 170 to 470 m/s. In comparison with the SCPTs are these
noticeable higher values.
7.3.5
3 measurement lines on site II are also tested by seismic refraction. Figure 7.9
shows the obtained profiles. A s-wave velocity of VS = 145 m/s was found in
the top layer until a depth of 2 m. This corresponds to a G0 of 40 MPa. The
p-wave velocity gives values below 1000 m/s.
Table 7.4 gives the obtained velocity values for the clay layer beneath 2 m
depth. The s-wave velocity averages to about 350 m/s; the p-wave velocity
0
No. 1
No. 2
No. 3
No. 4
No. 5
200
400
600
Swave velocity [m/s]
1
2
3
4
5
No. 1
d
Vs
[m] [m/s]
0.8
90
1.0
180
4.1
380
450
800
Profile A
Profile B
Profile C
Depth [m]
Depth [m]
91
No. 2
d
Vs
[m] [m/s]
0.3
170
0.5
100
1.1
160
4.3
340
500
200
400
600
Swave velocity [m/s]
800
470
No. 4
d
Vs
[m] [m/s]
0.6
90
1.2
170
4.2
350
470
No. 5
d
Vs
[m] [m/s]
0.3
200
0.5
110
1.1
160
4.4
300
700
7.4
7.4.1
The undisturbed samples from site II were tested in a triaxial cell under isotropic confining stress. The results of the analysis of the arrival of the shear
92
Vp
[m/s]
1667
1457
1669
Vs
[m/s]
320
348
385
[-]
0.48
0.47
0.47
Edyn
MPa
576
675
840
G0
MPa
195
231
286
Table 7.4: Seismic refraction test results for the clay material beneath a depth
of 2 m, site II.
160
300
140
250
120
200
VS [m/s]
G0 [MPa]
100
80
150
60
100
40
50
20
0
0
100
200
300
c [kPa]
400
500
0
0
600
Figure 7.10: Results of bender element tests, shear modulus G0 versus isotropic consolidation stress c ,
site II.
100
200
300
c [kPa]
400
500
600
Figure 7.11: Results of bender element tests, s-wave velocity VS versus isotropic consolidation stress c ,
site II.
wave are drawn in Figure 7.10 and Figure 7.11. Both diagrams show the same
results once in terms of shear modulus and the other in terms of shear wave
velocity. They can be linked if a density of 1988 kg/m3 is used.
7.4.2
A sample taken from the boring segment between 4.0 and 4.5 m depth of
boring B2I was installed in the free torsion pendulum test device of the Ghent
University. The dimensions of the sample were: 20 cm height, 10 cm diameter.
The isotropic stress stages have been chosen to 50, 100, 200 and 400 kPa. The
damping ratio at a reference shear strain of 105 was found in the range of 3.0 to
4.0 %. The smaller values correspond as expected to the higher confining stress
stages and the larger ones to smaller stresses. The modulus at the reference
shear strain is 35.0, 42.1, 53.4 and 69.2 MPa at 50, 100, 200 and 400 kPa
confining stress respectively.
Shear modulus and damping ratio in function of the applied shear strain can
be found in Figure 7.12. The same parameters in dependence of the confining
93
GModulus
90
80
70
G [MPa]
60
50
40
30
20
10
50 kPa
100 kPa
200 kPa
400 kPa
0
6
10
10
Strain []
10
10
Damping Ratio
0.08
0.07
0.06
D []
0.05
0.04
0.03
0.02
0.01
0
6
10
50 kPa
100 kPa
200 kPa
400 kPa
5
10
10
10
Strain []
Figure 7.12: G-Modulus and damping ratio versus applied shear strain obtained
by the Free Torsion Pendulum Test, sampling depth: 4.0 to 4.5 m, site I.
stress at a shear strain of 105 are given in Figure 7.13. The found increase
of the damping ratio from a shear strain of 2 106 to lower strains is not
expected. It is assumed that this is caused by limitations in the test equipment
to measure very small movements of the excitation beam.
7.4.3
A Boom clay sample has been sent to the private laboratory GEOLABO in
France to perform a resonant column test. The isotropic stress stages have
been again chosen to 50, 100, 200 and 400 kPa. Shear modulus and damping
ratio were measured during consolidation at a strain level of 105 . The values
of the shear modulus increase with the confining stress from 27.9 MPa at 50 kPa
isotropic stress to 32.3 MPa at 100 kPa, to 42.0 MPa at 200 kPa and finally
to 55.9 MPa at 400 kPa. The damping ratio is for all stress stage between 3.5
and 4.0 %. Both can be seen in Figure 7.13.
The damping ratio agrees very well with the findings from the free torsion
pendulum test. The shear modulus obtained by the free torsion pendulum test
is in all cases 20 to 30 % higher than the results from the resonant column test.
94
Damping Ratio
80
0.06
70
0.05
60
0.04
D []
G [MPa]
50
40
0.03
30
0.02
20
0.01
10
100
200
300
[kPa]
c
400
100
200
300
[kPa]
400
Figure 7.13: G-Modulus and damping ratio versus isotropic confining stress,
shear strain 105 , sampling depth: 4.0 to 4.5 m, site I.
7.5
Figure 7.14 summarizes the geotechnical test results for the site in SintKatelijne-Waver. The simplified profile of the shear wave velocity at site I
is based on both SCPTs. But there are no SCPT data available for depths
smaller than 0.8 m. The velocity for this shallow layer is assumed out of the
SASW data and the seismic refraction tests at site II. The soil material close
to the surface has to be similar at site I and II to justify this assumption.
Therefore the reliability of the shear wave velocity values in the first meter is
limited. Table 7.5 presents numerical values for the SCPTs.
Groundwater,
max. 2.80 m
min. 3.30 m
0.00 m
8.50
3.00
0.80
Reference level:
180
150
210
Vs [m/s]
100
210
Site Sint-Katelijne-Waver
11.00
Clay
3.00
Clayey Sand
Quaternary
Deposits
Formation of Boom
0.00 m
Depth [m]
Simplified S-Wave
Velocity Profile at Site I
surface site I
16.30 m
Groundwater
as at site I
3.00 m
2.30 m
Clay
Clayey
Sand
surface site II
14.00 m
0.70 m
0.00 m
Simplified Stratification
at Site II
Quaternary
Deposits
Formation of Boom
Simplified Stratification
at Site I
14
12
10
2
3
q [MN/m]
c
Averaged CPT-Profile
at Site II
96
Profiles of SCPTs
SCPT1 SCPT2
Depth
Vs
Vs
[m]
[m/s]
[m/s]
0.83
204.27
1.34
142.06
172.30
1.86
234.22
181.08
2.38
244.37
208.20
2.89
255.15
166.95
3.41
173.68
140.83
3.93
145.93
151.18
4.44
151.52
153.12
4.96
155.63
151.81
5.48
157.93
160.79
5.99
162.15
163.29
6.51
164.29
174.00
7.03
166.16
169.56
7.54
164.36
165.63
8.06
179.59
8.58
184.12
Chapter 8
Introduction
The test location in Ghent is situated close to a storage and garage building
of the Flemish Geotechnical Institute on the campus of the Ghent University
in Zwijnaarde. A SCPT is performed on this site. Boring and CPT data
are additionally available. They have been gathered in preparation of the
construction of the buildings of the geotechnical institute in the late 1960th.
A lay-out plan of the testing site is given in Figure 8.1 and the overview of
all used in situ test results can be found in Table 8.1.
Date
Abbr.
15/07/1968
B6
Depth
[m]
14.0
B21
5.0
19/07/1968
19/07/1968
25/11/2003
CPT21
SCPT1
14.75
17.5
8.2
The boring B21 in close proximity to the performed SCPT reached only to a
depth of 5 m. Therefore it is combined with the deeper, but also more far,
boring B6. The stratification of the combined profile at the place of boring
97
98
Parking Space
B6
Parking Space
CPT 21
B 21
SCPT 1
Storage and
Garage Building
B21 is given in Table 8.2 and in the site overview in Figure 8.5. Because of the
extrapolation, the accuracy of the layer interface-levels at depths deeper than
5 m might be rather low.
However, the site has a 0.75 m thick top layer of a silty sand followed by
fine sand to a depth of about 8.25 m. The sand is mixed from this depth on
with an increasing amount of silt. The material is classified as silty sand to
sandy silt.
All described soil layers are deposits originated in the Quaternary. The first
tertiary layer at the site is referred to the formation of Tilt. It is deeper than
the final depth of the evaluated boring. Its presence is seen in boring profiles
available in the Flemish Geotechnical Database on other locations inside the
university campus.
Groundwater table
The groundwater table is found at the boring B21 at 1.85 m depth.
99
Main
Admixtures
Component
fine sand
silt
fine sand
fine sand
silt
silt
fine sand
to 13.5 m quaternary deposits.
Table 8.2: Combined results of boring B21 and B6 for the location of B21.
8.3
A CPT sounding has been performed at the place of boring B21 at one corner
of the garage. The profile of the tip resistance of the cone qc is given in Figure
8.2.
0
Depth [m]
10
12
14
16
6
8 10
q [MN/m]
12
14
16
8.4
8.4.1
The SCPT was performed using a 200 kN CPT-tuck without track support.
Because of a very good signal quality the test was driven to a final depth of
about 17.5 m.
100
8.4.2
The shear wave velocity profiles obtained by cross correlation from both excitation directions are very consistently. They are given in Figure 8.3. A simplified
profile can be found in the summary in Figure 8.5.
8.4.3
The data of the SCPT are analyzed to obtain the damping ratio using the
spectral ratio slope method. The results are presented in Figure 8.4. The
values are mainly in the range between 0 and 10 % at all depths. The variation
is qualitativelly similar to the other test sites. For closer details on the variation
in the results can be refered to Section 4.4.3.
SCPT 1
SCPT 1
Depth [m]
Depth [m]
10
Left side
Right side
12
12
Left side
Right side
14
14
16
16
18
10
100
200
300
Swave velocity [m/s]
400
18
10
20
Damping Ratio [%]
30
8.5
101
Figure 8.5 summarizes the gathered geotechnical data for the site in Ghent.
The simplified s-wave velocity profile is only derived from the SCPT in Figure
8.3 because there are no other measurements of wave velocity performed at this
site. Table 8.3 presents numerical values for the SCPT.
(Silty Sand)
Fine Sand
13.50
(Sandy Silt)
Silt
10.80
8.25
(Sand)
Fine Sand
Depth [m]
Quaternary
Deposits
16
14
12
10
Groundwater
1.85 m
Fine Sand
0.00 m
8
10
q [MN/m]
c
12
14
Averaged CPT-Profile
at CPT21
16
18
16
14
12
10
100
200
300
400
S-wave velocity [m/s]
250
280
160
180
Vs [m/s]
Depth [m]
Simplified Stratification
at Drilling B21
17.00
7.80
3.00
0.00 m
102
CHAPTER 8. TEST SITE GHENT
103
Profiles of SCPTs
SCPT1
Depth Vs, Left Vs, Right
[m]
[m/s]
[m/s]
1.34 228.02
170.23
1.86 198.58
192.48
2.38 179.91
173.01
2.89 156.71
164.08
3.41 172.16
176.04
3.93 182.39
181.69
4.44 203.71
195.25
4.96 214.44
220.84
5.48 231.68
227.22
5.99 242.09
260.00
6.51 269.30
272.48
7.03 265.46
286.80
7.54 287.33
284.05
8.06 261.24
261.19
8.58 251.82
258.82
9.09 247.12
250.81
9.61 238.09
245.00
10.13 246.17
242.83
10.64 242.77
239.01
11.16 248.12
249.57
11.67 256.85
251.50
12.19 228.94
238.74
12.71 243.60
248.02
13.22 257.02
256.12
13.74 261.63
268.89
14.26 259.69
279.15
14.77 267.13
270.41
15.29 256.53
260.21
15.81 245.56
260.25
16.32 252.19
264.30
16.84 230.26
241.43
Table 8.3: Numerical values of Vs from SCPT1; excitation from the left and
from the right side.
104
Part III
105
Chapter 9
The seismic cone penetration test (SCPT) is developed from the cone penetration test (CPT) using electrical cones. The capabilities of an electrical cone are
extended by adding geophones or accelerometers to the usual transducers for
tip resistance, sleeve friction and pore pressure. Also cones that are limited to
the vibration sensors are in use. These cones are called dummy cones. Waves
generated at the surface travel to the cone pushed into the ground. The test is
similar to a down-hole test except for the borehole, necessary in the down-hole
test to install seismic sensors. The main purpose of the method is to obtain the
velocity of shear and compression waves leading to values for Gmax and Emax .
Objectives and methodology
Aim of this part of the research is to obtain at first confidence in the determination of wave velocity profiles. Afterwards the possibilities to use hidden
information in the time histories to calculate the damping ratio of the soil
should be investigated. The incorporation of the damping ratio measurement
would mean a significant improvement of the SCPT method. This is especially
true considering the fact that the few other approaches (e.g. SASW with extended inversion, cross-hole test) to measure material damping in-situ are still
under research and not state-of-the-art.
In order to achieve the aims the following topics are treated and described
in this chapter:
- Description of the test equipment
- Selection of suitable sensors for the installation in the seismic cone
- Summary of available direct and indirect time methods to measure wave
velocity, including example calculations
107
108
9.2
An overview of a typical seismic cone set-up is given in Figure 9.1; the equipment described in detail in the following paragraphs is shown in Figure 9.2.
Lateral Offset
Mechanical
Hammer
Static
Load
Sledge
hammer
tT
LT
tB
Depth
Beam
LB
Triaxial
Accelerometers
Seismic Cone
Penetrometer
9.2.1
Seismic source
The seismic source consists of a steel beam loaded by the weight of the CPT
truck and different types of hammers. A horizontal hammer blow at one end
of the beam produces a shear wave dominated pulse and a vertical blow on top
of the beam a compression wave dominated pulse.
109
6 channels
6 channels
X1, Y1,
Z1
3 channels
X1, Y1,
Z1
Geophones
X2, Y2,
Z2
Hamlin accelerometers or
geophones
Termination box
X2, Y2,
Z2
X1, Y1,
Z1
PC with data
acquisition card
Trigger
channel
Hammer and beam
(shearwav e source)
Single
seismic cone
Applied load
provided by the CPT truck
a)
Source beam
b)
Impact anvil
Applied load
provided by the CPT truck
Source beam
c)
Applied load
provided by the CPT truck
Impact anvil
Wooden beams
110
The SCPTs performed in the framework of this research use in most of the
cases a source beam with a system of rollers on top, under the loading weight of
the CPT truck. It is developed by Areias et al. (1999) and called a de-coupled
source. This increases the horizontal mobility of the beam resulting in a more
efficient transfer of the impact energy of the hammer to the actually generated
shear waves. Comparative studies between the de-coupled and the coupled
source without the roller system are found in Areias et al. (2004) showing the
positive effect of the rollers.
According to the same author there is an optimum vertical stress at the
boundary between beam and soil. Consequently, it is always attempted to
adjust the vertical force on the beam, applied by hydraulic cylinders, to its
optimum value. This is checked by comparing time records of additional hammer blows performed in the beginning of each SCPT with the cone at a fixed
shallow depth.
In the SCPT test on the university campus in Ghent two wooden plates are
successfully used instead of the steel beam. Their size suits more the demands
of an easy transport than the steel beam. Also in this case rollers under the
loading piston are placed. Figure 9.3 shows principle sketches of the source
beams.
The impact energy is given by a sledge hammer or a mechanical swing
hammer with an adjustable falling height. With the latter is aimed to generate
repeatable shear pulses of the same intensity. However, comparisons of repeated blows with the cone at a certain depth show that the assumption that the
mechanical swing hammer provides a repeatable source is questionable. Differences in the amplitudes and shapes of the generated shear pulse are observed,
probably caused by a slight backwards movement of the hammer support and
beam from blow to blow. A fixing of the hammer frame with metal pins on
the ground could not improve the result. Furthermore it can be assumed that
the load on the source beam changes with time due to a slowly sinking of the
CPT truck additionally causing a changing shear amplitude.
The acceleration of the head of the mechanical swing hammer is measured
on the occasions of the SCPT5 in Retie and SCPT1 in Sint-Katelijne-Waver.
A small shock accelerometer is installed on the backside of the hammer head.
The weight of the whole moving arm of the hammer is approximately 25 kg.
Because of the unknown weight distribution in the hammer arm, the impact
force was not calculated.
Four selected hammer impacts at the site in Retie are shown in Figure 9.4.
Figure 9.5 provides the associated amplitude spectra. Even so the spectra are
partly different, it can be concluded that the main energy of the hammer impact
is located between 200 and 600 Hz. It has to be remarked that the measured
time histories characterize only the impact force applied to the source beam.
The energy distribution of the dynamic force transfered to the ground might
be different.
Another point concerns the measurement itself: the signal is sampled with
2000 Hz, more than sufficient for the ground vibrations, but the frequency
111
content of the impact pulse approaches the Nyquist frequency of 1000 Hz.
Because there is no anti-aliasing filter applied to the hammer channel, it is
not completely impossible that existing signal parts at frequencies higher than
1000 Hz are wrapped into the high frequency parts of the spectrum just below
1000 Hz. However, this disturbance is probably limited to the higher frequency
range.
4
x 10
8
No. 295
4
2
0
2
4
6
8
0.48
0.5
0.52
Time [s]
0.54
4
2
0
2
4
6
8
0.48
0.56
8
No. 321
0.5
0.52
Time [s]
0.54
0.56
x 10
4
2
0
2
4
6
8
0.48
No. 309
x 10
0.5
0.52
Time [s]
0.54
0.56
x 10
No. 325
6
4
2
0
2
4
6
8
0.48
0.5
0.52
Time [s]
0.54
0.56
Figure 9.4: Acceleration of the mechanical swing hammer during the SCPT
No. 5 in Retie, selected impacts.
9.2.2
Seismic cones
112
80
70
60
No. 325
50
40
No. 309
30
20
No. 295
10
0
0
200
400
600
Frequency [Hz]
800
1000
Figure 9.5: Amplitude spectra of the selected hammer impacts from Figure
9.4.
Figure 9.6.
In the frame of this research both cones are used. The seismic sensors
in the dual cone are replaced two times in the run of the project. The first
change replaced the originally installed geophones by VTI Hamlin accelerometers because of insufficient frequency response characteristics of the geophones
as described later on. These accelerometers where later replaced by accelerometers from Br
uel & Kjr to improve the signal-to-noise ratio leading to larger
testing depths.
In the following the characteristics of the sensors used are described.
GeoSpace geophones
The triaxial array consists of three single geophones GSC-14-L3 produced by
GeoSpace. The geophones are designed and built to maintain performance
characteristics even after being subjected to high shock forces. The main specifications of the sensors given by the manufacturer are:
Natural frequency: 28 Hz 5 Hz
Moving mass: 2.15 g
Dimensions (diameter/height): 1.7 cm / 1.8 cm
Open circuit damping: 18 %
Intrinsic voltage sensitivity: 0.114 V/(cm/s), amplified by an internal
amplifier to 27 V/(cm/s)
35.682
(A=10 cm)
43.702
(A=15 cm)
47.873
(A=18 cm)
Triaxial Sensor
(top-part)
35.682
(A=10 cm)
Triaxial Sensor
(bottom-part)
113
60
Cone-tip
43.702
(A=15 cm)
189
1033
392
260
720
370
1482
400
Conetip
[mm]
51
230
133
Dual Cone
(1 m configuration)
35.682
(A=10 cm)
Triaxial
Sensor
Friction
transducer
60
Load cell/Pore
pressure transducer
446
630
Single Cone
Because the suspension of the geophones in the housing of the cones might
have an influence on the natural frequency and to validate the given sensitivity,
a dynamical calibration of the single cone has been performed in the Laboratory
for Mechanical Engineering at Ghent University. The testing set-up to calibrate
the vertically orientated geophone is given in Figure 9.7 and the horizontally
orientated sensors in Figure 9.8.
In both set-ups the cone is excited by a mechanical shaker using a homogeneous sine function and the time records of the output signal of the geophones
are stored. The movement of the cone is additionally measured by a well calibrated ICP accelerometer mounted at the surface of the cone in close proximity
to the geophones. At each frequency of the sinusoidal excitation the amplitude
of the signals from the geophone output and from the accelerometer are measured. After integrating the acceleration amplitude into velocity, the ratio to the
velocity amplitude of the geophone provides one point in the response curve,
describing the sensitivity dependency on the frequency. The whole procedure is
repeated with a changed excitation frequency until the desired frequency range
from 1 to 500 Hz is covered.
In an early stage a LVDT, an eddy current proximity transducer and a laser
interferometer are also used to provide reliable reference values for the cone
movement. However, the results show that the surface accelerometer delivered
the clearest signals easily to analyze. Figure 9.9 shows the obtained sensitivity
curves, Figure 9.10 the general calibration curve of the manufacturer. The
response curves of the horizontal geophones fall together in Figure 9.9, therefore
only one line is drawn.
114
Triaxial geophone
(inside)
Single Cone
Accelerometer
Accelerometer
Single cone
Shaker
Roller bearing
(rotation centre)
Leaf spring
115
6
Vertical Geophone
Horizontal Geophones
Sensitivity [V / (mm/s)]
50
100
150
200
250
300
350
400
450
500
Frequency [Hz]
Figure 9.9: Calibration curves of the seismic single cone equipped with geophones GeoSpace GSC-14-L3.
OUTPUT [V/(IN/SEC)]
10
10
GeoSpace
GSC14L3
10
10
10
10
FREQUENCY [Hz]
Figure 9.10: Calibration curve of the geophones GeoSpace GSC-14-L3 as provided by the manufacturer.
116
Dual cone
Reference
accelerometer
Roller bearing
(rotation centre)
Reference
accelerometer
Leaf spring
Triaxial
accelerometer (inside)
Triaxial
accelerometer (inside)
Shaker
Figure 9.11: Calibration of the dual cone with installed Hamlin accelerometers,
horizontal.
117
1.4
chosen characteristic sensitivity: 1.116 V/g
(exact at 23 and 50 Hz )
1.2
Sensitivity [V/g]
0.8
0.6
0.4
0.2
0
0
20
40
60
80
100
120
140
Frequency [Hz]
9.2.3
118
and applies an analogue bandpass filter to the signals with the edge frequencies
1 and 350 Hz (high-pass 1 Hz, 6 dB/oct; low-pass 350 Hz, 12 dB/oct). This
conditioner module of the seismic cone is connected to the termination box of
the data acquisition card inside the PC by conventional BNC-cables. The first
channel is used to trigger the measurement. The hammer and the beam are
connected to this port. A voltage of 5 V is applied to the channel inside the
termination box. When hammer and beam get in contact, the two poles are
bypassed over a resister and the voltage on the channel drops. This can be
used to generate a software trigger in the acquisition program.
In the case that the Br
uel & Kjr accelerometers are used the conditioner
unit is to be replaced by an ICP or Isotron source. At the test site in Ghent
a filter and amplifier system KEMO VBF35 provided by the K.U. Leuven is
used. The lowpass cut-off frequency is set to 1000 Hz. Alternatively the 9
channel source and amplifier of Ghent University ENDEVCO 2792B without
filter can be applied.
SCPTs using the single cone are performed with a sampling frequency of
4000 Hz; vibration events measured with the dual cone are sampled at 2000
Hz. So the Nyquist frequency is 2000, respectively 1000 Hz, much higher than
the 350 Hz cut-off frequency of the bandpass filter, or equal to the 1000 Hz cutoff frequency of the KEMO VBF35 system. Aliasing effects are consequently
impossible or in case of the KEMO amplifier sufficiently reduced. On the site
in Ghent at some depths measurements are repeated without analog low-pass
filter. The time histories show a similar quality as those with filter except for
some additional noise.
Main specifications of the data acquisition card manufactured by Microstar
Laboratories are:
Type: DAP 4000a/112, PCI bus
Number of analog inputs: 16
119
9.3
120
elastic threshold strain the modulus calculated with Equation (2.1) can be
assumed as Gmax or G0 .
The elasticity modulus Emax or E0 can be calculated out of the compression
wave velocity Vp using Equation (9.1) for an infinite elastic medium.
Emax = Vp2
(1 + )(1 2)
1
(9.1)
Vp
Vs
Vp
Vs
2
2
(9.2)
121
Stewart (1992) has shown that the consideration of a bended path has nearly
no influence on the calculated wave velocity and is therefore not taken into
account.
The main focus to calculate the wave velocity is the determination of the
travel time between two measurement points. The methods to achieve this can
be divided into direct and indirect time methods.
The signals of the horizontal sensors in the direction of the hammer impact
(y-axis) are used for the s-wave velocity determination. The vertical sensors
(z-axis) are used for the p-wave velocity determination.
9.3.1
Direct time methods are based on the visual interpretation of the time histories.
Characteristic points in the signals are identified and their shifts in time with
increasing depth are used to calculate the wave velocity.
The most obvious way is to determine the time of the first arrival of the por s-wave, but also the first peak or the first trough might be used. Another
possibility is to draw two signals from hammer impacts from both sides of the
source beam on top of each other. The shear wave component of the two signals
will be polarized. This can be seen in the example time histories in Figure 9.14.
Components not related to the generated s-wave will show the same polarity in
both signals. In this way the s-wave component can be identified more easily.
The reference point in time is the first crossing point of the two signals close
to the time of first arrival. The method is therefore called cross-over method.
The direct time methods represent the most straightforward way to obtain
the wave velocity. However, there are cases where the identification of the chosen reference point is difficult and requires some experience from the operator.
Beyond that is an automation of the first arrival or peak picking fault-prone.
Therefore the application of the methods is in general time consuming.
9.3.2
The most popular indirect methods are the cross correlation method, still performed in the time domain, and the method based on the phase of the cross
power spectrum as a frequency domain technique.
It is found in this research that the cross correlation in the time domain
works most effectively and reliably. In addition it provides the opportunity for
automation. Therefore the applied technique is explained in more detail using
example time histories.
Cross correlation method
The time shift, t = , between the signals of the bottom and the top accelerometer in the cone, sB and sT , is obtained at the maximum of the cross
122
0.4
Hammer blow left
Hammer blow right
0.2
0.2
0.4
0.5
0.52
0.54
0.56
0.58
Time [s]
0.6
0.62
0.64
0.66
0.4
Hammer blow left
Hammer blow right
0.2
0.2
0.4
0.5
0.52
0.54
0.56
0.58
Time [s]
0.6
0.62
0.64
0.66
Figure 9.14: Example time histories from the site in Ghent, cone tip at a fixed
depth, hammer blows from left and right.
rBT ( ) =
sB (t) sT (t + ) dt
(9.3)
t=
The difference in the travel path length divided by the time shift gives the
shear wave velocity.
Because the sampling frequency is limited, the accuracy of the time shift
t, obtained by the cross correlation is also limited. In most cases the sampling
frequency is set to 2000 Hz, corresponding to a sampling interval of 0.5 ms.
The calculated time shifts are consequently multiples of 0.5 ms. This leads, in
consideration of the small travel path length between the two receivers of about
1 m, to an insufficient accurate velocity calculation. The obvious solution for
this problem is to increase the sampling frequency. This is done in reality by
changing the settings of the analogue/digital converter in the data acquisition
unit or, more conveniently, by interpolating between the measured data points.
The latter technique allows to increase the sampling frequency artificially to
virtually any value.
123
The interpolation is based on the sampling theorem of the digital signal processing. The theorem implies that a band limited signal is completely described
if it is sampled with at least the double of the maximum signal frequency. This
is ensured by the selection of a sufficient high sampling frequency. Therefore
it is possible to reconstruct the signal between the sampling intervals unambiguously. With other words, it is possible to interpolate. Kammeyer & Kroschel
(1998) provide an equation for this interpolation for an arbitrary time t:
s(t) =
k=
s(k/fS )
sin[(fS t k)]
(fS t k)
(9.4)
with the signal s(t) originally sampled as s(k/fS ) at times k/fS . fS is the original
sampling frequency and k the sample index.
The signals are processed by another algorithm more straightforward to
implement. The source can be found in Weinstein & Clifford (1979) and is
realized in a MATLAB script. The main steps are the following:
1. expanding of the input data vector by inserting zeros between the original
data values
2. design of a special symmetric FIR (finite impulse response) filter that
allows the original data to pass through unchanged and interpolates in
between
3. application of the filter to the input data vector to produce the interpolated output vector
The time records obtained from SCPTs are interpolated to a 50 times
higher sampling rate. This means, for the case of fS = 2 kHz, an artificial
sampling frequency of 100 kHz, corresponding to a sampling interval of 0.01 ms.
This leads consequently to a finer time shift and also wave velocity resolution.
Figure 9.15 gives an example time history from the testing site in Ghent.
The actual acquired data points are indicated as dots, the line in between is
interpolated and further on used for the cross correlation.
Figure 9.16 shows the cross correlation of the left side time records given
in Figure 9.14. The maximum of the curve appears at a t of 4.23 ms. The
travel-path difference is 1.0304 m, giving a s-wave velocity of 243.60 m/s.
Phase of the cross power spectrum
The phase of the cross power spectrum can be used to analyze the variation of
the velocity with frequency. The technique is the basis of the SASW method.
After the calculation of the discrete Fourier transformation of both signals
sT (t) and sB (t) the cross power spectrum is obtained by Equation (9.7). SB (f)
is the complex conjugate spectrum of the bottom signal. The discrete Fourier
transformation is abbreviated for simplicity as FFT (Fast Fourier Transformation) even though the FFT is only a certain but most common numerical
implementation of the discrete Fourier transformation.
124
0.4
Acquired data points
Interpolated time history
0.2
0.2
0.4
0.54
0.55
0.56
0.57
Time [s]
0.58
0.59
0.6
Crosscorrelation [m/s]
x 10
t = 4.23 ms
4
2
0
2
4
0.1
0.05
0
Time shift [s]
0.05
0.1
0.15
Figure 9.16: The cross correlation of the signals from Figure 9.14; picked maximum at t = 4.23 ms.
125
ST (f ) = F F T (sT (t))
(9.5)
SB (f ) = F F T (sB (t))
(9.6)
RT B = ST (f ) SB (f )
(9.7)
Im(RT B )
Re(RT B )
(9.8)
The phase shift is wrapped, that means the values are limited by and .
Before the time shift t(f) between the two signal is calculated using Equation
(9.9), it has to be unwrapped.
t(f ) =
(f )
2 f
(9.9)
Finally the wave velocity V(t) is obtained in consideration of the travel path
difference x.
V (t) =
x
t(f )
(9.10)
Figure 9.17 shows the application on signals from the site in Ghent as used
in Figure 9.14. The velocity in the range of 50 to 150 Hz is nearly constant,
confirming the assumption that the s-wave is not dispersive, and is in good
agreement with the result of the cross correlation in the time domain. Outside
the frequency range of 50 to 150 Hz, the velocity fluctuates. This is in consideration of the frequency content of the signals also reasonable. Most of the
energy is concentrated in the range of 50 to 150 Hz. The calculation of the
coherence function would probably confirm this finding.
9.4
9.4.1
This paragraph will introduce the equations describing body waves and discuss
factors affecting the amplitude of these waves. This leads to a method to
quantify the wave attenuation from seismic measurements.
First consider a simple sine wave with no attenuation traveling along a
string with wavelength , velocity c and source amplitude A0 . t marks the
time and x the distance of the current observation point from the source, then:
A++ (x, t) = A0 sin
2
(x ct)
(9.11)
126
Amplitude [m/s]
2.5
x 10
2
1.5
1
0.5
0
50
10
300
Unwrapped phase
Wrapped phase
5
0
100
150
200
250
Wrapped / unwrapped phase of the crosspower spectrum
50
100
500
150
Swave velocity
200
250
300
250
300
400
300
200
Result of time domain crosscorrelation, V = 243.60 m/s
100
0
50
100
150
Frequency [Hz]
200
Figure 9.17: Analysis of the phase of the cross power spectrum; amplitude
spectra of the signals, phase shift and resulting s-wave velocity; travel-path
difference between the two signals is 1.0304 m.
127
(9.12)
(9.13)
For spherical waves in a homogeneous medium, neglecting near field terms, the
amplitude decays inversely with the distance x:
A+ (x, t) =
A0 i(kxt)
e
x
(9.14)
This attenuation is also called geometrical damping. Factors for other types
of waves can be found in Table 9.1.
Wave type
Body waves at depth
Body waves at the surface
Rayleigh waves
Table 9.1: Geometrical damping; the amplitude at a distance x from the source
is proportional to the tabulated values (Studer & Koller (1997)).
In order to introduce material damping caused by a not perfectly elastic
medium, the wave number is allowed to be complex: k = k + i. Equation
(9.14) becomes:
A+ (x, t) =
A0 x ik x it
e
e
e
x
(9.15)
A0 x
e
x
(9.16)
(9.17)
128
ln
= ln
x1
x2
(x2 x1 )
(9.18)
A1 x 1
A2 x 2
x2 x 1
(9.19)
2 f D
c
(9.20)
(9.21)
If the dual cone is used and the technique is applied to s-waves with velocity
Vs , Equation (9.21) changes to:
AT x T
Vs ln AB xB
Ds =
2f xB xT
(9.22)
xT and xB represent the travel path length from the source to the top, respectively the bottom sensor array; AT is the amplitude of the frequency spectrum
of the signal of the top sensor, AB is the corresponding amplitude of the bottom
sensor.
The time records of a dual cone SCPT given in Figure 9.14 are used to
give an example of the described theory. Figure 9.18 shows, based on Equation (9.19), the calculated attenuation coefficient versus the frequency. The
highest signal energy is in the frequency range between 90 and 130 Hz. This
can be seen in Figure 9.17. The curve of shows, neglecting some scatter, a
quite linear slope at the discussed frequencies.
This slope appears in the diagram of the damping ratio Ds given in Figure
9.19 as a region of an only slightly changing Ds . The mean value between 90
and 130 Hz is in this case 2.7 %.
9.4.2
The attenuation coefficient method can be adapted to make the results independent from the selection of a factor for the geometrical attenuation that does
not depend on frequency.
129
0.5
0.4
Attenuation Coefficient, [m ]
0.3
0.2
0.1
0
0.1
0.2
0.3
0.4
0.5
0
50
100
150
200
Frequency [Hz]
250
300
Figure 9.18: Attenuation coefficient for the example signals from Figure 9.14.
20
15
10
5
0
5
10
15
20
0
50
100
150
200
Frequency [Hz]
250
300
130
(9.23)
d ln
df
= z(x2 x1 )
(9.24)
It can be noted that the factor to correct the geometrical damping from
Equation (9.14) is eliminated by differentiating. Any geometric term affecting
the amplitude that does not depend on frequency will also be eliminated. This
will include the transmission and reflection corrections on layer interfaces if the
velocities are independent of frequency, which is the case if the frequency range
for the analysis is properly selected. If x1 is held constant and x2 (or simply
x) is varied, differentiating with respect to x gives:
2 ln
A2
A1
f x
= z
(9.25)
Techniques which are based on the Equations (9.24) and (9.25) are called
spectral ratio slope methods (SRS) and are described among others in Udas
(1999), Buttkus (1991), Kudo & Shima (1981) and Stewart (1992). z is linked
to the damping ratio D and the geophysical properties quality factor Q and
dissipation factor Q1 :
D=
1
Q1
zc
=
=
2Q
2
2
(9.26)
(9.27)
131
0.5
0.4
0.3
0.2
0.1
0
0.1
0.2
0.3
0.4
0.5
0
50
100
150
200
Frequency [Hz]
250
300
Figure 9.20: Logarithmic amplitude ratio of the signals from Figure 9.14, fitted
slope for the range of highest signal energy, corresponding damping ratio.
9.4.3
132
a1
Free ground
surface
Accelerometers
1
2
z1
a3
z2
a4
z3
z3
z4
a2
a5
a6
ai1
ai
ai+1
i1
i
i+1
zi1
zi
(z, t) =
u
(, t) d
(9.28)
=0
= i1 (t) +
(9.29)
(9.30)
in which subscripts i and (i-1/2) refer to levels zi (of the ith accelerometer)
and (zi1 + zi )/2 (halfway between accelerometers i and (i-1)) respectively,
i = (zi , t), u
i = u
(zi , t) and zi is the distance between the accelerometers.
The corresponding shear strains at levels zi and (zi1 + zi )/2 may be expressed as:
i (t)
1
zi1 + zi
(ui+1 ui )
. . . +(ui ui1 )
ui ui1
zi1
i = 2, 3, . . .
i1/2 (t)
133
zi1
zi
zi
zi1
(9.31)
(9.32)
fmax =
c
min
(9.33)
134
F F T (a2 )
A2
=
F F T (a1 )
A1
(9.34)
[m]
1.033
1.033
1.033
1.033
i1
m=1
Index SCPT
subsetup (n)
n=1
Dual cone
sensors
A1, 2
A1,3
Surface
sensor
H 2/3 =
Depth
Index (m)
An , m = FFT ( a n , m )
Pm = FFT ( p m )
A2,3
a1,3
A2, 4
H 3/ 4 =
a1, 2
H 4/5 =
a 2, 4
a 2,3
A3, 4
A3,5
H i 1 / i =
Ai 2,i 1
Ai 2,i
a 3, 5
a 3, 4
i2
ai 2 , i
ai 2,i 1
= IFFT A1,3
k =4
H k 1 / k
pi = IFFT ( Pi 1 H i 1 / i )
= ...
pi 1 = IFFT (Pi 2 H i 2 / i 1 )
Free ground
surface
= IFFT (A1,3 H 3 / 4 H 4 / 5 )
p5 = IFFT (P4 H 4 / 5 )
= IFFT (A1,3 H 3 / 4 )
p 4 = IFFT (P3 H 3 / 4 )
p3 = a1,3
p 2 = a1, 2
p1 = a1,1
136
137
G and D
for the selected loop
Figure 9.23: Application of the hysteresis method to data from the Waremme
site, SCPT3.
138
Figure 9.24:
SCPT3.
Figure 9.25:
SCPT3.
139
Figure 9.26: Normalized shear strain histories versus depth, Waremme, SCPT3.
140
100
Hysteresis-Loop II
Hysteresis-Loop I
150
150
Loop II
t = 0.523 ... 0.544 s,
G = 63.436 MPa
100 Vs = 176 m/s (with = 2040 kg/m),
D = 27.57 %
Loop I
t = 0.505 ... 0.524 s,
G = 87.127 MPa
Vs = 207 m/s (with = 2040 kg/m),
50
D = 25.75 %
0
Loop I
-150
-2
0
t = 0.544 s
t = 0.523 s
-50
-50
-100
50
t = 0.505 s
-1
t = 0.524 s
0
Shear strain, [-]
Loop II
-100
2
x 10
-6
-150
-2
-1
0
Shear strain, [-]
2
x 10
-6
141
Shear strain
x 10
Shear strain, []
Loop I
1
1
Loop II
2
0.45
0.475
0.5
0.525
0.55
0.575
0.6
Time, t [s]
Shear stress
200
Shear stress, [Pa]
Loop I
100
100
Loop II
200
0.45
0.475
0.5
0.525
0.55
0.575
0.6
Time, t [s]
Figure 9.30: Shear strain and shear stress time histories at a depth of 3.71 m
after sledgehammer excitation, Waremme, SCPT3.
142
gnal drifts. In this case a rectangular high-pass filter (FIR) is used. This filter
provides a unit gain factor over the desired frequency pass band and introduces
no phase shift in the filtering procedure. Also an IIR filter with smoother filter
flanks designed after Ormsby (1961) is tried. Both filters are only efficient in
removing the disturbing signal drifts if all components below about 10 Hz are
significantly reduced or removed. This means of course a loss of important
signal information.
Stress-strain calculation The receiver spacing between the dual cone sensors is not optimum. A low-pass filter has to be applied to remove frequency
components above 60 Hz. Unfortunately, the main wave energy is with 50 to
150 Hz in this frequency range.
9.5
Preliminary trials of the SRS method on SCPT data show a wide result scattering concerning the damping ratio. Because other in situ testing techniques
were not available it was not clear if this variation is due to real changes in the
soil properties with depth or caused by the algorithm of analyzing.
Therefore the seismic cone test is simulated using the program SPECTRAL
providing artificial input signals for the analyzing procedures.
The algorithm of analyzing applied on the measured signals at the testing
sites is also used to analyze the simulated signals. Essentially three artificial
soil profiles as given in Figure 9.31 are evaluated, a half-space, a one layer and
a two layer profile. Damping ratio, shear wave and compression wave velocity
are kept constant within one layer.
The geometrical conditions for a typical SCPT are also used for the simulated tests as given in Figure 9.32. In particular a horizontal off-set of the source
of one meter and a depth distance between the simulated signal records of 0.5
m are chosen.
The excitation hammer blow is simulated by a point load on the surface.
In the cases of profile a) and c) the direction of the force is once horizontal
and once vertical. For profile b) only the horizontal excitation is used. The
behavior of the excitation force follows a squared sine function and is given in
Figure 9.33.
The calculation of the p-wave velocity is performed using the vertical particle acceleration caused by the simulated vertical impact. The s-wave velocity
is obtained based on the horizontal particle acceleration due to the horizontal
impact force.
The damping ratio is calculated for the shear wave only, therefore the horizontal acceleration in connection with the horizontal impact is used.
The Figures 9.34, 9.35 and 9.36 show the time histories used for the s-wave
velocity and the damping ratio calculation. The direct s- and p-wave compo-
+/ 0.00 m
+/ 0.00 m
VS = 107.5 m/s
VP = 215.0 m/s
D = 1.0 %
143
+/ 0.00 m
VS = 107.5 m/s
VP = 215.0 m/s
D = 1.0 %
3.00 m
7.00 m
VS = 107.5 m/s
VP = 215.0 m/s
D = 3.0 %
II
VS = 215.0 m/s
VP = 430.0 m/s
D = 2.3 %
III
VS = 107.5 m/s
VP = 215.0 m/s
D = 1.0 %
10.00 m
II
15.00 m
VS = 152.1 m/s
VP = 304.2 m/s
D = 2.3 %
15.00 m
a) HalfSpace
15.00 m
b) 1 Layer Profile
c) 2 Layer Profile
Excitation
Axis
Cone
Axis
1m
v ertical
horizontal
+/ 0.00 m
Source
0.5 m
Calculation
Points
0.5 m
15.00 m
144
Normalised force []
15
20
25
t [ms]
30
35
40
Normalised force []
200
400
600
800
1000
f [Hz]
1200
1400
1600
1800
2000
Figure 9.33: Time history and frequency content of the simulated normalized
impact force.
nent is indicated. There is an unexpected component in profile b) between pand s-wave arrival assumed to be caused by an internal calculation instability
in the SPECTRAL program algorithm. Besides this, reflections on the layer
interfaces can be very well seen. Also nicely seen, at least in profile c), is
that the refraction of the s-wave component at a layer boundary generates an
additional wave component traveling with p-wave velocity.
The complex structure of the time records of the 2 layer profile c) emphasizes
the need to implement a possibility to select parts of the signals and not to
use the whole signal for the velocity and damping analysis. It is expected
that at least the common treatment of wave components of different speed
will cause significant errors in the calculation of the damping ratio as the wave
velocity is an essential part of the SRS method and the ACM. The wave velocity
calculation based on cross correlation should be not affected as long as the
amplitudes of the component of interest are higher than the others.
The analysis of the simulated signals is performed through a couple of MATLAB script files. The algorithm is outlined in Figure 9.37. The SRS method
and the ACM are applied to signals from depths with one meter distance,
meter-by-meter approach. The SRS method is additionally tried on two signals from a depth of 5 and 12 m. The slope of the spectral logarithmic ratio
is obtained by a linear curve fit in a frequency range of 300 Hz always centered at the frequency of the maximum spectral amplitude. The damping ratio
145
a) HalfSpace
Depth [m]
10
12
14
16
0
0.05
0.1
0.15
0.2
t [s]
Figure 9.34: Simulated time histories for profile a) Half-Space, horizontal excitation.
146
b) 1 Layer Profile
Layer I
Depth [m]
6
S
Layer II
10
12
16
0
14
0.05
0.1
0.15
0.2
t [s]
Figure 9.35: Simulated time histories for profile b) one layer on a half-space,
horizontal excitation.
147
c) 2 Layer Profile
Layer I
Depth [m]
Layer II
10
Layer III
12
14
16
0
0.02
0.04
0.06
0.08
0.1
0.12
t [s]
Figure 9.36: Simulated time histories for profile c) two layers on a half-space,
horizontal excitation.
148
determined by the ACM is an average value for the same frequency range. To
study the influence of the different length of a cutting window, different shapes
of window functions and the meaning of additional noise to the results of the
half-space of profile a), some parameters are changed which will be explained
further on.
Input data from
SPECTRAL program
Signal adaption
Length adjustment of the windowed
signal part (E)
Adding of noise (N)
Applying of Hanning (H) or tapert (T) window
AdaptFFT.m
Meterbymeter approach
Cros s Correlation
Cros s Correlation
autoSCPTFromWinSig.m
autoSCPTFromWinSig3.m
Vs v s. Depth
autoSCPTFromWinSig2.m
Vs for the selected
depth area
Spectral Ratio
Slope Method (SRS)
dualmanu.m
Attenuation Coefficient
Method (ACM)
dualmanu3.m
Spectral Ratio
Slope Method (SRS)
dualmanu2.m
D v s. Depth
D v s. Depth
Figure 9.37: Flow chart for the calculations on the simulated SCPTs.
Window-length
The signal from the half-space profile a) is analyzed using three different window lengths. First the complete signal without window is used, secondly a
signal of a length of the first 1.5 cycles of the s-wave component is used which
is virtually the whole s-wave component and finally a signal containing only
the first cycle of the s-wave component is analyzed. Figure 9.38 illustrates the
149
shape of the signal after the window is applied. The windowing is done manually by selecting the crossing points with the neutral position. In this way it
can happen that the two proceeded windowed signals have a slightly different
length. To investigate influences of this fact all calculations are done additionally with signals of exactly the same length. For this case it is accepted that
one of the signals contains more than one respectively 1.5 cycles. Results of
exact equal signal length are indicated with E in the summary Table 9.2.
1.5 Cycles
1 Cycle
Figure 9.38: Schematic view of the results after windowing of 1.5 cycles (1.5
periods) and one cycle (1 period) of the signal.
Window-function
Hanning window Additionally to the simple rectangular shaped window, a
tapered window and a Hanning, also called Hann window, is applied. The shape
of the window function is given in Figure 9.39. The Hanning window follows
a cosine function and is defined by Equation (9.35). n describes the length of
the window in number of samples; k is the sample index. The Hanning window
is a common window type and for instance used by Kudo & Shima (1981) for
damping measurements on a down-hole setup with two boreholes.
w(k + 1) = 0.5 1 cos 2
k
n1
, k = 0, . . . , n 1
(9.35)
150
Time domain
Frequency domain
50
Hanning
T apered
Rectangular
Magnitude (dB)
Amplitude
0.8
0.6
0.4
-50
0.2
5
Time (ms)
10
15
-100
0
0.2
0.4
0.6
0.8
Frequency (kHz)
Figure 9.39: Hanning, tapered and rectangular window with a typical length
of 1.5 cycles in time and frequency domain.
151
9.5.1
2.5
Depth [m]
5
Vp
7.5
Vs
10
12.5
15
50
100
150
200
Wave velocity [m/s]
250
300
152
2.5
2.5
Vs
5
Depth [m]
Depth [m]
7.5
7.5
V
10
10
12.5
12.5
15
50
100
150
200
Wave velocity [m/s]
250
15
300
Figure 9.41: Calculated s-wave velocity, profile b), one layer on a halfspace.
100
200
300
Wave velocity [m/s]
400
500
Figure 9.42: Calculated wave velocity, profile c), two layers on a halfspace.
and shape of the applied window as long as the wave component of interest
has a higher amplitude than all others. The input velocities are nearly exactly
found except for the shallow depths smaller than 2.5 m in case of the profile
a). This is probably due to the fact that the difference in travel path length
between two measurement points close to the surface is rather small, so small
deviations in the time arrival determination have an unequal higher disturbing
influence in the resulting velocity.
9.5.2
The input parameter for the damping ratio is Ds = Dp = 1% over the whole
depth of profile a). Table 9.2 and Figure 9.44 summarize the results concerning
the damping ratio for the half-space profile a). The Figures 9.45 to 9.51 show
the results for selected parameter configurations. Table 9.2 gives the value of
the averaged damping ratio over all n obtained points of the profile between 5.5
and as a measure for scattering, the mean absolute deviation
and 11.5 m, D,
sD0 from the target value D0 = 1%. The latter is defined as:
sD0 =
1
n
n
i=1
|Di D0 | .
(9.36)
153
In the case where the signals at 5 and 12 m are compared the single calculated value D is given. sD0 is then simply the deviation from the target input
value.
As an example Figure 9.43 shows for the calculation point at 10 m depth the
influence of the length of the rectangular window to the shape of the amplitude
spectra of the signals and the logarithmic ratio between both.
Normalized
Amplitude []
0.5
Normalized
Amplitude []
200
400
600
800
1000
1200
1400
1600
1800
2000
Spectrum at depth 9.5 m, AT
200
400
600
800
1000
1200
1400
1600
1800
2000
Logarithmic Ratio, ln(A /A )
0.5
1.5
Logarithmic
Ratio []
1
0.5
0
200
400
600
800
1000
1200
Frequency [Hz]
1400
1600
1800
2000
154
cN
1 cycle
cET
cEH
cT
cH
cE
c
bN
1.5 cycles
bET
bEH
bT
bH
bE
Full signal
aN
Selected-Depth, SRS
aT
Metre-by-Metre, ACM
aH
Metre-by-Metre, SRS
a
0.000
0.200
0.400
0.600
0.800
1.000
Figure 9.44: Mean absolute deviation of the damping ratio profiles obtained
from the simulated SCPTs.
Metre-by-Metre
Method
SRS
Id-Code
Tapered Window
Hanning Window
ACM
D
[%]
sD0
[%]
D
[%]
sD0
[%]
D
[%]
sD0
[%]
a
aH
aT
aN
1.003
1.002
1.003
1.025
0.009
0.005
0.009
0.225
0.975
0.505
0.975
0.968
0.025
0.495
0.025
0.111
1.002
1.002
1.002
1.063
0.002
0.002
0.002
0.063
b
bE
bH
bT
bEH
bET
bN
1.137
0.987
0.869
1.360
0.682
0.910
1.139
0.147
0.032
0.148
0.398
0.318
0.092
0.432
1.895
0.978
1.644
1.842
1.037
0.983
1.896
0.945
0.024
0.665
0.888
0.044
0.022
0.985
1.029
1.108
0.815
1.089
0.718
1.174
1.078
0.029
0.108
0.185
0.089
0.282
0.174
0.078
c
cE
cH
cT
cEH
cET
cN
0.580
0.469
0.360
0.496
0.271
0.329
0.584
0.566
0.645
0.640
0.593
0.729
0.671
0.586
1.281
1.173
0.993
1.034
1.031
1.108
1.273
0.476
0.179
0.581
0.414
0.466
0.395
0.483
0.570
0.774
0.480
0.561
0.569
0.761
0.638
0.430
0.226
0.520
0.439
0.431
0.239
0.362
x
x
x
x
x
x
x
x
x
x
x
x
Selected-Depth
Area Method
Noise
c) 1 cycle
b) 1.5 cycles)
a) Full signal
Window length
SRS
x
x
155
x
x
156
9.5.3
This parameter study has shown that overall the influence of the type of window
is rather small. It seems useful to choose an exact equal window length for both
signals. However, the main conclusion is, that the length of the applied window
has the most significant influence on the result and its scatter. As shorter the
window is chosen as larger the deviation and as larger the scatter too.
157
2.5
Depth [m]
7.5
10
12.5
SRS a
ACM a
15
0.5
0.5
1
1.5
2
Damping ratio [%]
2.5
Figure 9.45: Damping ratio, profile a), based on the full signal.
158
2.5
2.5
Depth [m]
Depth [m]
7.5
10
12.5
0.5
1
1.5
2
Damping ratio [%]
2.5
15
0.5
0.5
1
1.5
2
Damping ratio [%]
2.5
2.5
2.5
Depth [m]
Depth [m]
7.5
SRS c
SRS cE
SRS cH
SRS cT
10
7.5
ACM c
ACM cE
ACM cH
ACM cT
10
12.5
15
0.5
ACM b
ACM bE
ACM bH
ACM bT
10
SRS b
SRS bE
SRS bH
SRS bT
12.5
15
0.5
7.5
12.5
0.5
1
1.5
2
Damping ratio [%]
2.5
15
0.5
0.5
1
1.5
2
Damping ratio [%]
2.5
159
2.5
2.5
Depth [m]
Depth [m]
SRS a
SRS aN
SRS b
SRS bN
SRS c
SRS cN
7.5
10
10
12.5
15
0.5
7.5
ACM a
ACM aN
ACM b
ACM bN
ACM c
ACM cN
12.5
0.5
1
1.5
2
Damping ratio [%]
2.5
15
0.5
0.5
2.5
Depth [m]
2.5
7.5
D = 2.3 %
10
D = 2.3 %
7.5
10
SRS
SRS E
ACM
ACM E
12.5
12.5
5
Damping ratio [%]
SRS
SRS E
ACM
ACM E
D = 3.0 %
D = 1.0 %
Depth [m]
2.5
15
5
1
1.5
2
Damping ratio [%]
10
15
15
5
D = 1.0 %
5
Damping ratio [%]
10
15
160
9.6
Details on the practical test setup on the five research sites can be found in the
Chapters 4 to 8. They include also results for the wave velocities and damping
ratio profiles. The SCPTs are divided in three measurement campaigns: campaign I) includes tests with the geophone equipped single cone, campaign II)
are tests with the dual cone and the Hamlin accelerometers and campaign III)
covers a SCPT with the dual cone and the Br
uel & Kjr accelerometers.
Site
Retie
Lincent
Waremme
Sint-Katelijne-Waver
Ghent
I
SCPT1 to SCPT4
SCPT1 to SCPT3
SCPT1 to SCPT2
Campaign
II
SCPT5
SCPT4 to SCPT5
SCPT3 to SCPT5
SCPT1 to SCPT2
III
SCPT1
9.6.1
Campaign I (geophones)
The first tests were performed using the single cone. The arrival time of the
s-wave component and exceptionally the p-wave velocity is visually determined.
Because the resonant frequency of the geophones lays in the frequency range
of interest the test data are not used for a damping calculation. It would also
be necessary to calculate the ratio of the spectra of time records related to
different hammer blows, because of only one sensor array in the cone. Because
of the limited reproducibility of the mechanical swing hammer this method is
less promising than measurements with two sensors at the same time
9.6.2
The Hamlin accelerometers installed in the dual cone are used for these measurements.
Algorithm of data processing
The analysis of the test data follows in general the flow chart in Figure 9.54.
The SRS method coupled with the meter-by-meter approach is applied according to the algorithm used for the simulated SCPTs. The attenuation
coefficient method is not tried.
Because the real measurements contain noise and other influences additional
steps for the data preparation are added.
161
Windowing / Filtering
manual / semiautomatic signalpartselection
(rectangular window)
poss. lowpass filtering (200 Hz cutoff frequency)
srsdatanew.m
Meterbymeter approach
CrossCorrelation
autoSCPTFromWinSig.m
autoSCPTFromWinSig3.m
Vs vs. Depth
Spectral Ratio
Slope Method (SRS)
dualmanu.m
D vs. Depth
Figure 9.54: Flow chart for the calculations on the data from the SCPT testing
sites, SRS method.
162
Depth [m]
Depth [m]
After the conversion of the ASCII data, gathered with the LabVIEWprogram, taking into account the sensors sensitivity and possibly applied gain
factors by the conditioner equipment, all signals are visually checked. Usually
5 repetition of excitation hammer blows at each cone depth are recorded. This
opens the possibility to skip obviously disturbed hammer blows that have been
overlooked during the measurement.
The remaining time records are averaged in the time domain to improve
the signal to noise ratio. A trial with averaging in the frequency domain led to
virtually the same averaged signals.
However, even after averaging of up to 5 signals the obtained time records
have still remarkable noise components. Figure 9.57 shows an averaged time
record from the site in Sint-Katelijne-Waver. Figure 9.55 shows the averaged
signals according to their depth.
8
10
8
10
12
12
14
14
16
16
18
0.45
0.5
0.55
t [s]
0.6
0.65
0.7
18
0.45
0.5
0.55
t [s]
0.6
0.65
0.7
The high noise level suggests the application of a low-pass filter. In this
way it is more straightforward to find reasonable points for the beginning and
the end of the rectangular window. A corner frequency of 200 Hz is the best
compromise between smoothing the signal and keeping its frequency content.
The filtering is done before the signals are windowed.
Afterwards the main cycle, following the same definition given in Figure
9.38 (1 cycle), of every signal is cut by a rectangular window.
Because most of the signal energy is situated in the frequency range between
163
0.1
0.05
0.05
0.1
0.55
0.555
0.56
0.565
0.57
0.575
Time [s]
0.58
0.585
0.59
0.595
0.6
0.59
0.595
0.6
0.1
0.05
0.05
0.1
0.55
0.555
0.56
0.565
0.57
0.575
Time [s]
0.58
0.585
Figure 9.57: Example of averaged time histories from the top accelerometer,
y-direction, Depth 11.67 m, Sint-Katelijne-Waver (SCPT2) before and after
filtering.
164
75 and 130 Hz, the slope of the logarithmic spectral ratio is obtained by fitting
a straight line in this range.
The s-wave velocity calculated by the cross correlation technique is found
to be independent of the application of the low-pass filter. Also if the signal is
not windowed the velocity profile stays the same.
The results of s-wave velocity and damping ratio are given in the Figures
9.58 to 9.65. The same results are already included in the previous chapters
on the test sites Retie, Lincent, Waremme and Sint-Katelijne-Waver, Chapter
4 to 7.
Retie The site in Retie is characterized by a sand profile. The soil changes at
about 5 to 6 m depth from a fine sand to the coarser sand of the formation of
Mol. This boundary layer is confirmed by a change in the shear wave velocity
from averaged 215 to 335 m/s.
The results of the damping ratio show major variation. The shallow layer
until about 6 m depth has a mean damping ratio of 4 %, below this layer the
value exceeds 10 %. This means the layer interface fine sand to sand is also
identifiable in the damping ratio measurements.
Although comparative damping values found in literature, based on laboratory tests and a shear strain of smaller than 103 %, differ a lot, a damping
ratio of more than 10 % is not expected. Seed & Idriss (1970) designate a D
of 1 to 4 % as a resonable value for the low strain range in sands.
Retie
Retie
Depth [m]
Depth [m]
Mechanical Hammer
Sledge Hammer
10
10
12
12
50
100
150
200
250
300
350
400
14
0
10
20
30
40
165
Lincent
0
Mechanical Hammer
Sledge Hammer
Depth [m]
Depth [m]
100
150
200
250
300
350
400
10
166
Waremme
0
Depth [m]
Depth [m]
10
10
12
12
14
16
0
50
100
150
200
250
300
Mechanical Hammer
Sledge Hammer
Resonant Column Test
Free Torsion Pendulum Test
14
350
16
0
10
9.6.3
The SCPT on the campus of the Ghent University is done with the Br
uel &
Kjr accelerometers, providing a much higher resolution than the Hamlin sensors and so an improved signal-to-noise ratio. Figure 9.56 shows the normalized
time records of the test in Ghent. The difference to the time histories obtained
with the Hamlin sensors on the site in Sint-Katelijne-Waver given on the left
side is quite obvious. While the test in Sint-Katelijne-Waver had to be stopped
at a depth of 13 m because of a bad signal quality, the test in Ghent provides
until the final depth of about 18 m signals with very small noise components.
Because of this good signal quality no low-pass filter has to be applied. Regarding to the experiences with small window lengths on the simulated signals,
Depth [m]
Depth [m]
Sint KatelijneWaver
0
5
6
5
6
9
Mechanical Hammer (SCPT)
Sledge Hammer (SCPT)
Resonant Column Test
Free Torsion Pendulum Test
10
11
0
167
50
100
150
200
250
10
300
11
0
10
different window lengths are tried. The procedure is the same as for campaign
II given in Figure 9.54. The averaging of up to 5 signals is done in the same way
as in campaign II. The slope of the logarithmic spectral ratio is obtained by
fitting a straight line in a range of 150 Hz on both sides of the amplitude peak
of the frequency spectrum. This is the same range as used in the simulated
signal analysis. Because the frequency of maximum energy is in general below
150 Hz, the application of the rule leads practically to a range from 0 to about
250 Hz.
The results of the s-wave velocity given in Figure 9.68 (identical with Figure
8.3) are again independent from the applied window. The damping ratio in
Figure 9.69 (identical with Figure 8.4) is calculated with the SRS method
with meter-by-meter approach and no window is used, so averaged signals are
analyzed without further modifications on the full signals.
For comparative reasons the SRS method is also applied with rectangular
windowed signals of one cycle length and on the full signal. Additionally a
rectangular window is tried with 0 from the beginning of the time record to the
first appearance of the shear wave component and 1 for the rest of the signal.
This window aims to exclude components arriving before the shear wave. Such
fast components, possibly p-waves, with remarkable amplitudes are found at
some depths.
As an example Figure 9.67 shows the amplitude spectra used for the dam-
168
ping calculation at the depth of 12.71 m. Figure 9.66 gives the resulting profiles
of the damping ratio obtained after applying the three window lengths.
Ghent, SCPT 1
0
2.5
Depth [m]
7.5
10
12.5
15
17.5
Full signal
2nd part
1 cycle
20
22.5
10
0
10
Damping ratio [%]
20
30
Figure 9.66: Influence of the length of the rectangular window to the damping
profile on the Ghent site.
The profile based on the full signal agrees essentially with the profile based
on the second part of the signal, at least for greater depths. The profile of the
one-cycle-window shows some deviation.
The site in Ghent, the only one covered by campaign III, has an extended
layer of fine sand. From about 8 m depth to greater depths it changes stepwise
to a silt. The shear wave velocity increases with depth in the superficial fine
sand. In the underlaying slightly coarser material it is almost constantly 250
m/s.
The damping ratio is again unexpectedly high. Even though the majority of
the calculated values is below 10 % some of them fall also in the range between
10 and 20 %. The damping ratio seems to be not depth dependent.
It is remarkable that the damping ratio at both sites dominated by sand,
Retie and Ghent, has unexpected high values, not agreeing with values given
by other authors. However, the results at the sites characterized by silt and
clay give reasonable result ranges.
Beside this it cannot be hidden that the result scattering is not satisfying.
Logarithmic Ratio []
Normalized Amplitude []
Normalized Amplitude []
169
Spectrum at depth 13.22 m, AB
Full signal
2nd part
1 cycle
0.5
50
100
150
200
250
300
Spectrum at depth 12.19 m, AT
50
100
150
200
250
300
Logarithmic Ratio, ln(AT/AB)
0.5
0
1
0.5
0
0.5
50
100
150
Frequency [Hz]
200
250
300
9.7
170
Ghent
Left side
Right side
Depth [m]
Depth [m]
10
10
12
12
14
14
16
16
18
0
50
100
150
200
250
300
350
400
18
0
Left side
Right side
10
15
20
25
30
171
Caparison tests using the resonant column and the free torsion pendulum
tests on samples from Waremme and Sint-Katelijne-Waver show a good agreement with the SCPTs concerning the s-wave velocity. The range of damping
ratios obtained by the SRS method matches with the results from the laboratory tests.
An outline of the main findings and conclusions:
Accelerometers have because of their flat response in the low frequency
range in contrast to miniature geophones advantages if accurate vibration
measurements are needed. The finally installed Br
uel & Kjr accelerometers are highly suitable for this purpose.
The averaging of time histories is an important tool to improve the signal
to noise ratio.
Vs can be accurately determined by the cross correlation technique.
It is questionable if a mechanical swing hammer with an adjustable falling height can be considered as a highly repeatable source. An increasing
settlement of the source beam during a SCPT changes the coupling stress
beam-soil and therewith the energy transfer characteristic from the hammer to the soil.
An artificial increase of the sampling frequency to 100 kHz improves
the Vs determination. By this way the minimum sampling frequency
of 20 kHz, requested by the guideline for SCPT of the ISSMGE TC 10
(2004), can be fulfilled even if the physical sampling frequency is far
below.
The application of a window function to a time record should be restricted. The s-wave velocity determination on real measured signals needed
no windowing in the case of the testing sites. The application of the SRS
method to calculate the damping ratio needs a windowing but it should
be as long as possible. It is useful to choose exactly the same length for
both compared time histories.
A window function other than the rectangular window brings no further
improvement.
172
Chapter 10
174
10.2
10.2.1
Bender elements
Bender elements consist of two thin piezoceramic plates rigidly bonded to a central metallic plate. Two thin conductive layers, electrodes, are glued externally
to the bender.
The polarization of the ceramic material in each plate and the electrical
connections are such that when a driving voltage is applied to the element, one
plate elongates and the other shortens. The net result is a bending displacement. On the other hand, when an element is forced to bend, an electrical
signal can be measured through the wires leading to the element.
There are two possible versions, differing only in the electrical connection
of the two polarized plates: a series or a parallel connection (Figure 10.1). As
the series version gives a higher output for a given distortion, it is better to
use it as a receiver. The parallel version has an additional electrode between
the piezoceramic plates and is better used as transmitter because it gives the
largest distortion for a given input.
A transmitter and receiver element can be placed in various soil testing
devices, such as oedometer, simple shear test device or, as in this research, in
the bottom and top cap of a conventional triaxial device.
Although a combination of a parallel element as sender and a serial element
as receiver would improve the signal quality, in the frame of this research only
175
b)
a)
Figure 10.1: (a) Series and (b) parallel connected piezoceramic element (Dyvik
& Madshus (1985)).
series connected elements are used. This gives the opportunity to change the
pulse direction in the sample without the element configuration influencing the
receiving signal. It means the received signals should have the same shape
independent of the sending direction.
The electromechanical properties of the series connected bender elements
are given in Table 10.1.
Material:
Dimensions,
length:
width:
thickness:
Resonant frequency,
unmounted:
capacitance:
Maximum drive voltages,
low frequency:
at resonance:
used range 0.1 - 15 kHz:
Sample penetration,
bottom cap:
top cap:
Lead Zirconate
Titanate (PZT-5HN)
20 mm
10 mm
0.5 mm
610 Hz
6.5 nF
50 V
15 V
20 V
1.6 mm
2.4 mm
176
10.2.2
The equipment to perform pulse bender element tests is given in Figure 10.2.
The setup for a test with a continuous excitation with a sample of a small
height, as it has been used in the later explained resonant method, is depicted
in Figure 10.3.
Signal Analyzer
Top cap
Source
Bender
elements
Signal Analyzer
Ch.1 Ch.2
Sample
Amplifier
IN OUT TRIGG.
Source
Top cap
Bottom cap
Var.: trigger on
step pulse
Bender
elements
Ch.1 Ch.2
Amplifier
IN OUT TRIGG.
Sample
Bottom cap
Var.: trigger on
sending signal
Signal generation
A HP 3562A Dynamic Signal Analyzer provides a continuous sine output signal.
The frequency can be step-less chosen between 0 and 100 kHz. Because the
amplitude of the output is limited to 5 V, an additional power amplifier is
necessary.
Alternative signal generation Instead of the signal analyzer, not able to
generate a sine pulse directly, a sound card of a computer is also tried as signal
source with success. It provides a very cost-effective way to produce signals of
arbitrary characteristics, including sine pulses. The sound card is steered by a
MATLAB script but also a control by LabVIEW can be implemented. Because
the output amplitude of the sound card is also limited to about 5 V a subsequent
amplification is unavoidable. However, to obtain the results presented in the
frame of this research the sound card method is not used, it is stuck to the
signal analyzer as source.
Signal amplification
The driving signal is amplified by a specially developed power amplifier to an
amplitude of 20 V (40 V peak-to-peak amplitude). The amplifier is also able
to cut single sine pulses from the continuous sine to provide a pulse source.
Simultaneously to the sine pulse it also generates a function of rectangular
177
shape with the same length as the circle duration of the sine pulse. This signal
is intended to trigger the data acquisition, because it is supposed to be perfectly
in phase with the sine pulse. It is shown later that the analyzer for the data
acquisition also can directly trigger on the sine pulse making the trigger pulse
unnecessary and eliminating the incorporation of triggering errors caused by
observed small deviations of the starting times of the rectangular trigger pulse
and the sine pulse.
If a continuous signal with amplified amplitudes is desired the pulse cutting
function of the device can be switched off. The signal analyzer and the amplifier
are shown in Figure 10.4.
Signal acquisition
The receiver signal is sampled by the same dynamic analyzer used already for
the generation of the continuous sine function. It is able to capture signals at
two channels. The first channel is in general connected in parallel to the signal
that is driving the sender element and is used to trigger the measurement.
The second channel is connected to the output of the receiver element. If the
rectangular pulse of the amplifier is used as trigger signal, it is connected to
the first channel instead of the driving signal. The collected data is transferred
to a computer by means of a floppy disk.
A function of the analyzer to average repetitions of the receiver signal is
used to reduce the signal noise during pulse tests. Usually 20 single pulses are
averaged and summarized to one signal for the further interpretation.
To calculate a response function between the output and the input signal
of the bender elements, the swept sine mode is used. If the dynamic analyzer
178
is set to this mode and a frequency span with a step width is given, the device
changes the source frequency automatically, calculates the ratio between source
and receiver amplitudes and draws the response curve. Such a test takes about
1.5 minutes.
Triaxial cell
The test triaxial cell, subjecting the sample to an isotropic confining water
pressure, is shown in Figure 10.5. It is a conventional triaxial cell connected to
a cell- and back-pressure supply panel and has a graduated burette to measure
the volume change. Bender elements are mounted at the center of the bottom
and top cap.
A diagram of the bender element mounted in the base of the cell is presented
in Figure 10.6.
4
5
3
9
7
Figure 10.6: Bottom cap with bender element and compression wave element;
1) compression wave transducer, 2) isolation support, 3) epoxy, 4) bender element (shear wave transducer), 5) latex, 6) wire leads, 7) bottom cap, 8) porous
stone, 9) drainage.
The mass of a standard metal upper cap (about 600 g) is several times larger
than the mass of the soil specimen. Therefore it was necessary in application of
the later described resonant method, to redesign the metal cap and to replace
it by a new one, made of plastic material with a total mass of about 51 g. For
the test in resonant mode the plunger is disconnected from the top cap to allow
its free movement in all directions.
It was very useful to connect the electrical ground of all components of the
measuring system to the metal parts of the cell housing.
10.3
10.3.1
The smallest distance between the bender elements is assumed to be the length
of the wave travel path needed for the calculation of the wave velocity. This
179
is the distance between the tip of the transmitter element and the tip of the
receiver element, also called the tip-to-tip distance. This assumption is based
on the work of Dyvik & Madshus (1985), confirmed by studies using different
element penetration depths by Brignoli et al. (1996), for instance.
10.3.2
Although pulse tests with bender elements started with a step function, as seen
in Dyvik & Madshus (1985), a sine pulse with a shape of a single sine cycle
is now most commonly used. It opens the opportunity to study the influence
of the frequency on the signals. Occasionally also driving signals consisting
of several sine cycles and deformed sine cycles, as suggested by Jovicic et al.
(1996), are used to simplify the identification of the s-wave arrival. For special
signal analyzes also continuous sine signals are applied.
10.3.3
Figure 10.7 shows a typical receiver signal obtained from a bender test using a
single sine pulse of 4.5 kHz for the silt material of the testing site in Waremme.
As it will be explained later the arrival time corresponds with the first inflexion
at point D. The example time record shows beside the first arrival also clearly
a second and a third wave arrival. However, in the majority of bender element
measurements it is not possible to identify such multiple arrivals.
20
Input
Signal
0.25
10
A
C
0
A
0
0.5
0.5
10
0.5
20
0.5
1.5
2.5
Time [ms]
3.5
4.5
Figure 10.7: Bender element record on silt from the Waremme site, depth:
11.0 -11.5 m, 200 kPa isotropic effective stress, pulse frequency f = 4.5 kHz,
Vs = 229 m/s, arrival time at point D: 371 s, travel path length 8.5 cm, input
and output signal.
In the following, methods for the determination of the travel time of the
shear wave through the sample are summarized. Because the s-wave velocity
is not changing with frequency, the simple visual interpretation of the receiver
signal is still the most commonly used technique. Sanches-Salinero et al. (1986)
developed an analytical solution for the time record at a monitoring point
that would result from a transverse sine pulse of a point source within an
180
infinite isotropic medium. Jovicic et al. (1996) made finite element calculations
confirming the analytical calculation. The essential findings are the definition
of a phenomenon called near field effect, explained in detailed later in this
section, which is the reason for a deflection, as labeled in Figure 10.7, direct in
front of the arrival point. Viggiani & Atkinson (1995) and Jovicic et al. (1996)
concluded from their laboratory work and from the numerical study that the
most consistent results are obtained by assuming that the shear wave arrival
is marked by this first inflexion labeled as D. Therefore in the framework of
this research this point is chosen in all cases, if it is identifiable. Additionally
the frequency of the sine pulse is varied to check if the selected point is not
moving with frequency. So the test is repeated at 3, 6, 9 and 15 kHz sine pulse
frequency.
First direct arrival in the output signal
The travel time can be estimated as the time between the start of the input
voltage pulse to the sending bender (point O) and the first deflection in the
output signal from the receiver bender (between point C and D). Depending
on the bender element signals, the identification of the first deflection point
can be straightforward or difficult. The arrival point of the shear wave front
can be masked by faster traveling p-wave components, the described near field
phenomenon and by an inductive coupling between the sender and receiver
channels.
Travel time between characteristic points of input and output signals
The travel time of an impulse wave between two points in space may be taken
as the time between characteristic points in the signals recorded at these two
points. Two measured receiver signals are needed for this approach.
Because in the case of a bender element test only one actual receiver signal
is available, it has to be assumed that the input signal to the transmitting
element has the same shape as a fictive receiver signal at the source. The input
signal is then used as the first receiver signal.
This assumption is questionable because of two reasons. The first is the capability of the bender elements to follow the driving signal. Studies of different
authors as for instance by Jovicic et al. (1996), applying self-monitoring elements, show that the ability of the bender elements to follow the driving signal
decreases in the very high frequency range. This is no problem for single and
multiple excitation sine pulses since the commonly used frequencies are below
this range. But if rectangular step function are used, such high frequencies are
reached.
The second problem is related to the transfer characteristic element soil
and the attenuation in the soil itself. Even though the transfer characteristic
mechanism is not well understood it can be seen from practical tests that
both phenomena cause a change of the frequency content of the signal. The
compared signals, i.e. driving signal and one receiver signal, show therefore
181
a different shape also in the time domain. This again makes the picking of
comparable points at the two signals often erroneous. The method is therefore
not used in this work.
It has to be remarked that this doubt is also valid for other methods comparing a driving and a receiver signal, especially the cross correlation of the
two signal types in the time domain.
The most commonly used characteristic points are the first peak (point A
and A), first trough point (B and B) or zero crossings of the input and output
signals.
Travel time out of cross correlation of input and output signal
Another way for the determination of the travel time between sender and receiver element is the calculation of the cross correlation between sender and
receiver signal. The equation for the cross correlation is given, analogue to the
technique used in the SCPT, by:
rRS ( ) =
sR (t) sS (t + ) dt
(10.1)
t=
where sS is the sender signal, sR the receiver signal and t the time. rRS will
reach a maximum value for the time shift that equals the travel time of the
impulse from source to the receiver.
If reflected or refracted waves are strong and the cross correlation fails, a
windowing of the receiver signal can be useful.
Travel time out of the phase of the cross power spectrum of input
and output signal
Instead of the cross correlation in the time domain also cross correlation in the
frequency domain might be applied. This method leads, again very analogue to
the described technique for the SCPT, using the wrapped and unwrapped phase
of the cross-power spectrum to the wave velocity as function of the frequency.
Because the shear wave velocity is in general not frequency dependent, that
means not dispersive, an average value of the velocity over the frequency range
with a reasonable high signal energy, can be calculated. Viggiani & Atkinson
(1995) describe the steps involved in this method and apply it to bender element tests. Their results show that the travel times calculated using the cross
correlation and the cross-power-spectrum method are almost identical.
Travel time using multiple arrivals in the output signal
Output signals show in some cases a clear second arrival or even multiple arrivals of the input wave. The second arrival is the input wave after it reflects
from the receiver cap, travels back to the transmitter cap where it reflects again
and then returns to the receiver cap a second time. Arrivals of a higher order
are caused by successive refections on receiver and top cap. Assuming plane
182
wave propagation, the time between the multiple arrivals in the output signal
is equal to multiples of the double of the travel time from cap to cap.
The travel time may be determined using either characteristic peaks or
the cross correlation method. Possible characteristic points in the example
time record are A-A-A, B-B-B or D-D-D. The cross correlation of
the output signal with itself, also called an auto-correlation, is in principle
sufficient to determine the travel time. If the different arrivals are not very
well separated it is useful to create, by windowing, for each wave arrival a
dummy signal, where all parts outside the time window containing the arrival,
are set to zero. The dummy signals are than cross-correlated which each other.
Arulnathan et al. (1998) found this method helpful in their studies.
The advantage using multiple arrivals in opposite to the comparison of
sender and receiver signal is the avoidance of the assumption that the electrical
input signal is transfered perfectly into a mechanical oscillation of the bender
tip. That means for instance that the method is also open for a step pulse
excitation.
However, multiple arrivals are not very often observed in the bender element
output signal. Therefore the field of application is quite restricted.
Travel time out of the phase shift of continuous cycled bender tests
To remove uncertainties in the identification of the travel time, recent innovations have led to continuous wave cycles being used. The technique is described
by Blewett et al. (1999) and a FE interpretation is given in Hardy et al. (2003).
In this method a continuously cycled sine voltage is applied to the transmitter element at a low frequency. The output from the receiver is displayed
together with the driving voltage. The frequency of the input voltage is increased until the input and the output signals come into phase. If the initial
frequency was low enough, under this conditions exactly one complete wavelength is created between the bender element tips. The frequency is again
increased until a couple of full wavelengths is obtained. Knowing the frequency
of the input voltage, the number of full wavelengths and the tip-to-tip distance,
the shear wave velocity can be calculated. In the case that it is uncertain whether input and output are fully in phase or not, the cross correlation can help
to determine the phase shift and to correct the result.
10.3.4
There are effects which make the identification of the s-wave arrival difficult, if
not impossible. These are shortly described now. It is often possible to reduce
their influence by adapting the measurement configuration.
P-wave components
Bender elements generate besides the transversely polarized shear waves also
a certain amount of compression waves, traveling with p-wave velocity. This
183
velocity is faster than the s-wave and can reach, if the sample is fully saturated, the compression wave velocity of the pore water. Therefore these wave
components arrive before the s-wave and might overlap with the actual s-wave
arrival.
Brignoli et al. (1996) show a collection of typical receiver signal shapes,
including partly significant p-wave components. It is noticed that this kind of
disturbances are more frequent in measurements with shear-plate transducers
and rather seldom seen if bender elements are used.
Near field effects
This phenomenon was first treated by Sanches-Salinero et al. (1986) by means of
an analytical model. A transverse sine pulse of a point source within an infinite
isotropic elastic medium is modeled and the time record at a monitoring point
at some distance is calculated. It was found that the resulting wave consists
of three coupled wave components with different characteristics. All three
components represent transverse motion, but they propagate with different
velocities; the first two (1 and 2 ) with the velocity of a shear wave and the
third (3 ) with the velocity of a compression wave. For the three components
the attenuation arising from the geometrical damping occurs at different rates,
2 and 3 attenuating an order of magnitude faster than the 1 . 1 represents
the far field pure shear wave. 2 and 3 are the near field components, because
they have only significant amplitudes in close proximity to the source.
The initial polarity of 3 is opposite to 2 and more important also to
1 . This agrees very well with the observations in typical time records of
experimental data as shown in Figure 10.7. The deflection caused by the near
field component is always opposite to the first deflection of the s-wave. The
arrival of the near field component 3 is almost identical with the p-wave
arrival. This is confirmed experimentally by Brignoli et al. (1996) for dry
specimens. In saturated soils the arrival of the near field component does,
following Brignoli et al. (1996), not show a strong initial excursion from the zero
axis so the arrival time cannot be clearly determined. Nevertheless, the near
field component arrival gives much lower velocities than p-waves in the fluidsoil system, measured with compression transducers. Therefore it is assumed
that the near field component travels with the speed of the p-wave of the dry
specimen.
Even so the analytical solution is based on a single sine pulse the results
are transferable to other excitation shapes such as step pulses.
Brignoli et al. (1996) conclude from the findings of Sanches-Salinero et al.
(1986) that it is in general useless to reverse the electrical impulse to the transmitter element for the purpose of reversing only the shear wave arrival (1 )
since also 2 and 3 reverse.
In order to separate the near field from the far field Sanches-Salinero et al.
(1986) introduced the parameter Rd . The amplitude of the near field component in relation to those of the first s-wave arrival becomes less significant with
184
d
df
=
Vs
(10.2)
where d is the travel distance, is the wavelength and f is the frequency of the
input signal. Because of the increasing deviation of the frequency of the received s-wave arrival from the input pulse frequency with increasing frequency,
Brignoli et al. (1996) use the frequency of the first sine cycle of the received
s-wave arrival.
While Sanches-Salinero et al. (1986) suggested Rd 2, Brignoli et al. (1996)
obtain results without near field component from Rd > 4 to 5 only. Viggiani &
Atkinson (1995) specify that the near field effect may mask the arrival of the
s-wave when Rd = 0.25 to 4. Arroyo (2001) calculates a lower limit for Rd of
0.625, based on an allowed error in the Vs determination of less than 5 %. The
minimum frequency in function of the wave velocity and the Rd value for a
sample in a typical triaxial cell with 11.7 cm tip-to-tip distance can be selected
using Figure 10.9. Figure 10.8 gives the same diagram for a sample with 2.6 cm
height, later used to obtaine the damping ratio.
30
7 6
5
Rd=3
25
Rd=2
Rd=4
20
R =1
f mean [kHz]
10
Rd = d /
5
R =0.25
d
0
0
100
200
300
Vs [m/s]
400
500
600
Figure 10.8: Mean frequency of the received signal fmean in function of Vs . and
Rd for a sample with 26 mm distance between the bender tips.
If signal interpretation is difficult due to the near field phenomenon, an
increase of the excitation frequency will make the arrival point identification
easier. A high excitation frequency has also the advantage that the received
185
30
10
7
6
25
5
R =d/
d
20
R =4
f mean [kHz]
15
Rd=3
Rd=2
10
Rd=1
Rd=0.25
0
0
100
200
300
Vs [m/s]
400
500
600
Figure 10.9: Diagram like Figure 10.8 for 117 mm tip-to-tip distance.
186
10.4
10.4.1
Resonant method
The work of Brocanelli & Rinaldi (1998) describes a method to measure the
damping ratio and shear wave velocity, using bender elements while working
in the frequency domain.
The idea of the method is to bring a short cylindrical soil sample in a shear
movement. For that reason the bottom bender element is excited with a steady
sine signal of constant voltage and the amplitude is measured at the receiver
element. To make this value independent from the source amplitude it is normalized by this amplitude. This process is repeated at different frequencies
until the whole spectrum of sample response is defined. The damping ratio is
estimated at the points of the curve around the natural frequency of the shear
mode. For this purpose different techniques are available such as the half-power
bandwidth and the more general circle-fit method, the latter is also using the
phase components of the resonant curves. Preliminary tests on a clay sample
showed that it could be difficult to find the correct peak in the response spectrum corresponding with the shear mode. So it is useful to have knowledge of
the frequency range were this mode dominates.
The main tasks to apply this method are: a) finding a suitable way to
pinpoint the shear mode frequency and b) the calculation of the damping value
from the response spectrum at that frequency.
Natural frequency of the shear movement
Analytical formulations A cylindrical soil sample placed in the triaxial
cell can be considered in a simplified assumption as a fixed-free beam with the
mass of the top cap attached at the free end. Brocanelli & Rinaldi (1998) have
shown by means of a finite element analysis of a 3 cm high sample that, with
a relatively light mass at the top, the sample deforms in perfect shear in the
first mode and with an increasing mass at the top, the shape of the deformed
specimen tends to be flexural. For a pure shear deformation an expression
for the angular natural frequency of the first shear mode (s ) can be derived
from the wave equation for a shear wave propagating in a rod and the force
equilibrium at the boundaries:
k mT
s L
s L
=
tan
m
Vs
Vs
(10.3)
where m is the mass of the top cap, mT is the mass of the sample, L is the
sample length, k is the shear factor and Vs is the shear wave velocity. The latter
is connected to the shear modulus G and the unit mass by the expression
already mentioned:
Vs =
(10.4)
187
A detailed formulation can be found for instance in Graff (1991) and Brocanelli & Rinaldi (1998). For long samples or large masses of the top cap, the
first resonant mode of the beam will be predominately of the bending. For
the case of a pure bending deformation Cascante et al. (1998) obtained the
following equation from the Rayleigh approximation:
f2 =
3 EI
L3
33
140 mT
h
+ m 1 + 3L
+
9
4
h 2
L
(10.5)
E is the modulus of elasticity, I the area moment of inertia and h the distance
between the end of the sample and the centroid of the mass of the top cap.
The basic concept in the Rayleigh method is the principle of conservation of
energy. To apply the Rayleigh procedure, it is necessary to assume the shape of
the system in its fundamental mode of vibration. This assumption of a shape
function effectively reduces the system to a SDOF system. Thus the frequency
of vibration can be found by equating the maximum strain energy developed
during the motion to the maximum kinetic energy. Detailed information concerning Rayleighs method can be found in Clough & Penzien (1993).
Because the real behavior of the sample is always influenced by shearing
and bending it could be difficult to find criteria to decide what expression is
more suitable. Therefore the Rayleigh expression was extended by a term for
the horizontal shear displacement. This led to the following equation:
2 = 42
60 EI L2 + 5 GA k 2
mT + m
(10.6)
with
= 198L5 + 231L3 + 70L 2
and
The derivation of Equation (10.6) can be found at the end of this section.
The actual shape of the mode is supposed as the superposition of a bending
and a shearing displacement. depends on the ratio between these two parts.
The unit of is m2 . A is the area of the cross section. The resonant frequency
can be calculated at those where becomes a minimum. This follows from
the consideration that any shape other than the true vibration shape would
require the action of additional external constraints to maintain equilibrium.
These extra constrains would stiffen the system, adding to its potential strain
energy, and thus cause an increase in the computed frequency. Consequently,
the true vibration shape will yield the lowest frequency obtainable by Rayleighs
188
method. The solution of the extreme value problem is possible but leads to a
very large expression and is therefore not given here.
Another approach to calculate the natural frequency of the soil sample is
to use the thick beam theory, the so-called Timoshenko beam theory. The
Timoshenko beam includes the bending effects, and shear and rotary inertia
effects. A suitable solution of the system of differential equations of a cantilever
beam with a mass elastically mounted at the free end can be found in Rossit
& Laura (2001). Simplified to the case of a rigid mounted mass, this solution
is given by the following determinant expression:
0=
+ L
sin L
cos L
L
mr =
m
mT
L2
2
L
L
sinh L
L cosh L
mr
cos L
L sin L
L2
2 mr
cosh L
L sinh L
(10.7)
with
F = 2
2 =
2 2 1
L4
= 2 + 2 L2
=
H+
T2 L3 mT
EI
H = (1 + )
I
AL2
2(1 + )
k
2
L2
= 2 2 L2
H 2 4F
2
H +
H 2 4F
m
mT, A, I, , k
(10.8)
189
(10.9)
(10.10)
The mode shape for the shearing displacement yS is assumed to be linear. This
is described by:
yS = x for x L
(10.11)
yS = L for x > L
(10.12)
and
(10.13)
(10.14)
and
The maximum internal potential energy Epot is computed by taking into consideration the internal energy of the bending and shearing component of the
deformed specimen:
Epot =
1
EI
2
L
0
1
(y )2 dx + GAk
2
(y )2 dx
(10.15)
= 6EI2 L3 . . .
12 2 5
1
L + 2L3 + 2 L
5
2
. . . + GAk
(10.16)
The maximum kinetic energy Ekin is computed by presuming harmonic oscillation with frequency . The component for the specimen is
Ekin, mT =
1 2
A
2
L
0
y 2 dx
(10.17)
190
Ekin, mT = mT 2
(10.18)
The kinetic energy for a concentrated mass added at a distance h above the
specimen can be evaluated using Equation (10.14) at x = h + L:
Ekin, m =
1
2
m yh+L
2
2
(10.19)
and
Ekin, m
1
m 2 (42 L6 + 122 L5 h + 4L4 . . .
2
. . . + 6L3 h + 92 L4 h2 + 2 L2 )
(10.20)
Finally, the circular resonant frequency for the flexural mode is computed by
equating the maximum internal potential energy Epot and the maximum kinetic
energy Ekin, mT + Ekin, m :
2 = 42
60 EI L2 + 5 GA k 2
mT + m
(10.21)
with
= 198L5 + 231L3 + 70L 2
and
Top Cap
Sam ple
mT
x
E ,G=
G=0
m
mT
E ,G=
E =0
a)
b)
191
(10.22)
where u = ur er + u e + uz ez is the displacement vector in cylindrical coordinates r, , z and and k are the circular frequency and the wavenumber,
respectively. The circumferential order n is an integer that describes the symmetry of the waveform around the z-axis. The modes with n=0 (axisymmetric)
can be decoupled in two families: torsional (T(0,m)), with azimuthal displacement only, and longitudinal (L(0,m)), with both axial and radial displacements.
The modes without axial symmetry are called flexural (F(n,m), with n 1) and
exhibit all three displacements coupled. The three mode families are shown in
Figure 10.12.
The circumferential order n specifies for n > 0 the number of cycles around
the z axis. The index m is used to sort the modes for a given family, in ascending
phase speed. It gives the number of counterrotating sections along the radius.
192
Torsional
ur
uz
Longitudinal
ur
uz
Flexural
Figure 10.12: Torsional, longitudinal and flexural modes in a cylindrical waveguide (Seco et al. (2002)).
(10.23)
where Df is a 6x6 matrix. The cylinder parameters are the outer radius b, the
p-wave velocity in a rod c0 (c0 = E/) and the Poisson ratio . For given
cylinder parameters and circumferential order n, the result is a relationship
between the frequency and the wavenumber, which can be solved to find a
finite number of propagation modes, with real wavenumbers, and an infinite
number of evanescent modes, which have complex wavenumbers and do not
propagate. This leads to the conclusion that all modes except the fundamental
modes of each mode family are only present above a cut-off frequency, different
for each mode.
Because the practical computation is done using a set of MATLAB scripts
developed by Fernando Seco and because the Pochhammer-Chree approach
plays only a secondary roll in this research, the reader is referred to Eringen
& Suhubi (1975), Thurston (1992) and Thurston (1978) for closer details in
the mathematical theory. The latter article provides also some insights in the
effects appearing if the waveguide is not of infinite length.
Calculation of cut-off frequencies and dispersion curves for the
silt from the Waremme site A computation based on the properties of
the tested reconstituted silt sample S1 from the Waremme site is done using
a set of MATLAB scripts (PCDISP) provided by Fernando Seco. Information on these tools can be found for instance in Seco et al. (2002). The input parameters for the program are a cylinder radius b = 2.5 cm, a density
= 1984 kg/m3 , a rod p-wave velocity Vp = c0 = 348.4 m/s and a Poisson ratio = 0.29. The parameters correspond to a s-wave velocity measured in the
specimen of Vs = 216.9 m/s.
193
Objective in the calculation are the dispersion curves of the modes below
10 kHz and the cut-off frequencies of every mode. Flexural modes are calculated
for n = 1 and n = 2 only.
Figure 10.13 shows the obtained frequency dependent phase velocities of
every identified mode and Table 10.2 gives the related cut-off frequencies.
Torsional
Type
fcut
[kHz]
T(0,1)
0.0
T(0,2) 7.0916
Modes
Longitudinal
Type
fcut
[kHz]
L(0,1)
0.0
L(0,2) 5.2910
L(0,3) 5.3678
L(0,4) 9.6875
Transversal
Type
fcut
[kHz]
F(1,1)
0.0
F(1,2) 2.5424
F(1,3) 3.8733
F(1,4) 7.3620
F(1,5) 8.7852
F(1,6) 9.7086
F(2,1) 3.2417
F(2,2) 4.2175
F(2,3) 6.0311
F(2,4) 9.2602
Table 10.2: Calculated cut-off frequencies, fcut , for the example silt material
S1.
Even so Figure 10.13 is quantitatively valid only for the used input values,
the qualitative findings are valid also for other soil properties. Some essential
findings are therefore pointed out in the following:
- The fundamental torsional mode T(0,1) is the only mode showing no
dispersive behavior. It travels with the velocity of a shear wave in a
rod, equal to the shear wave velocity in an infinite medium. This is the
typical mode used in a resonant column device with torsional excitation
and agrees with the elementary theory of torsional vibration in a rod.
- Higher torsional modes start asymptotic at a certain cut-off frequency
and at high frequencies they approach asymptotic the s-wave velocity.
- The basic longitudinal mode L(0,1) agrees at low frequencies with the
velocity of a compression wave (c0 = Vp (rod)) in the thin rod model.
However, at high frequencies the velocity approaches the Rayleigh-wave
velocity Vr .
- Higher longitudinal rods have a cut-off frequency and approach at infinite
high frequencies the s-wave velocity Vs .
0
0
100
200
300
400
500
600
700
800
900
1000
V (half space)
r
F(1,1)
F(2,1)
T(0,1)
L(0,1)
V (infinite medium/rod)
s
c0 =Vp (rod)
V (infinite medium)
p
F(1,2)
L(0,2)
F(2,2)
5
6
Frequency [kHz]
F(1,3)
F(2,3)
T(0,2)
F(1,4)
L(0,3)
F(1,4)
F(2,4)
10
L
(0,4)
194
CHAPTER 10. BENDER ELEMENT TECHNIQUE
195
Vs
d
(10.24)
196
z y
x
Figure 10.14: Discretization of sample and
top-cap.
elasticity is chosen for the behavior of the soil sample and the cap. The latter
is modeled with the actual properties of the plastic material, i.e. a mass of 51
g, an elasticity modulus of 3400 MPa, a height of 19 mm and a diameter of
50 mm. The frequency extraction procedure of ABAQUS is used to calculate
the natural frequencies and mode shapes. This procedure uses the subspace
iteration eigenvalue extraction technique. The modeling of the bender elements
itself and the excitation is not necessary for this purpose.
The first five modes together with a verbal description are given in Table
10.3. This table is valid for a tested reconstituted silt sample S1 at 100 kPa
confining pressure, as described in Section 10.5.2, and for the later defined
added mass coefficient CA = 0.
From the FEM-calculation it is seen that the lowest mode is indeed a flexural
mode dominated by shear deformation. The frequency of the second mode
seems to be sufficiently higher so not interacting with the first mode.
Mode 1
Mode 2
Mode 3
Mode 4
Mode 5
875 Hz
flexural
1301 Hz
torsional
2171 Hz
longitudinal
2354 Hz
flexural
4251 Hz
torsional,
2nd order
Table 10.3: Natural frequencies found by finite element calculation for sample
S1, 100 kPa and CA = 0.
197
6(1 + )
7 + 6
(10.25)
198
decreases with the distance from a fixed bottom of the sample. The k values
for the different sections are again averaged taking into account its height of
influence, hinflu. .
Section
1
2
3
4
5
6
7
8
9
10
z
[mm]
0.40
0.80
1.60
5.66
9.72
13.77
17.83
21.89
25.94
30.00
hinflu.
[mm]
0.40
0.80
2.43
4.06
4.06
4.06
4.06
4.06
4.06
2.03
k
[-]
1.0487
1.0254
0.9551
0.8552
0.7824
0.7556
0.7609
0.7781
0.7642
0.7412
Average:
=
k
0.8039
d2
4
(10.26)
Cross-section
199
Center-line
150
150
200
20
100
Isolines
labelled in N/m
30
30
y
x
0
10
250
250
200
150
10
y [mm]
z = 0, k = 1.0487
200
150
20
20
10
0
10
x [mm]
20
300
200
100
0
30
30
20
10
0
x [mm]
10
20
30
20
0
150
140 90
130
1
120 1
Isolines
labelled in N/m
30
30
120
160
12 1150
110 13400
0
20
110
10090
10090
120 1
2
160150 0
160
90
10
130
10400
15 16
y
x
150
0
14
130
110
90
10
y [mm]
140130 1110
1020
150
0
160
0
11010130
90 140
20
10
0
10
x [mm]
20
300
200
100
0
30
30
20
10
0
x [mm]
10
20
30
40
1
8000
60
16
0
20
160
40
12
0
y [mm]
160
14
121000 0
8
10
y
x
120
140
40
80
100
10
60
80
100
140 120
20
Isolines
labelled in N/m
30
30
20
10
0
10
x [mm]
20
30
300
200
100
0
30
20
10
0
x [mm]
10
20
30
Table 10.5: Shear stress distribution over the sample cross-section and at the
x-axis at z=0/5.66/13.77 mm, k-value for each height.
200
Cross-section
Center-line
60
120
100
80
80
120
100
80
80
60
100
120
10
40
y
x
0
40
y [mm]
100
60
10
120
20
20
Isolines
labelled in N/m
30
30
20
10
0
10
x [mm]
20
300
200
100
0
30
30
20
10
0
x [mm]
10
20
30
100
60
100
80
20
60
80
10
20
y
x
100
80
40
80
60
40
60
10
y [mm]
20
20
Isolines
labelled in N/m
30
30
20
10
120
0
10
x [mm]
20
30
300
200
100
0
30
20
10
0
x [mm]
10
20
30
Table 10.6: Shear stress distribution over the sample cross-section and at the
x-axis at z=21.89/30.00 mm, k-value for each height.
201
1
1 + (d/l)2
(10.27)
Equation (10.27) is given for the case of a moving cylinder with two free ends.
Even though the considered system of sample and top cap is fixed at one side
it can be assumed that this equation gives a reasonable approximation. The
length l is assumed to be the total length of sample and top cap.
So resonance frequency calculations, including the effects of the confining
water, can be easily performed using a virtual mass composed of the actual
body mass and the above given added mass.
Viscous material damping
The basis for the analysis of the frequency response of the soil sample is the
identification of different modes of vibration at resonance. The damping ratio
D is calculated at these points of the response spectrum in the neighbourhood
of a resonance peak.
Amplitude of the response spectrum The most common method of measuring damping uses the relative width of the response spectrum. Using the
variables indicated on the curve in Figure 10.16, the logarithmic decrement
and the damping ratio D, can be calculated from (Richart et al. (1970)):
2D
f22 f12
=
=
2
2
fm
1 D2
A2
A2max A2
1 2D2
1 D2
(10.28)
When D is small, the last term and (1 D2 )0.5 can be taken as equal to
1.0. Further simplification is possible, if A is chosen equal to Amax /20.5 :
D=
f22 f12
2
4fm
(10.29)
202
Figure 10.16: Resonant curve with variables for half-power bandwidth method.
20
22 12
2 tan 22 + 1 tan 21
(10.30)
with:
- 0 , angular frequency corresponding to the maximum sweep angular velocity;
- 1 , 2 , angular frequencies;
- 1 , 2 , angles at both sides of 0 .
A circle is fitted to the points of the response curve close to the resonant
frequency to find the center. Knowing this point makes it possible to determine
the necessary angles .
Multi-degree of freedom methods The frequency response function from
a linear system is defined as the ratio between the frequency domain representations of the signal leaving the system Y() and the input signal X(). The
function can be described further as a ratio of two polynomials, as shown in
Equation (10.31).
H() =
Y ()
=
X()
m
k=0 ak
n
k=0 bk
(i)k
(i)k
(10.31)
203
Im
Re
c
1
1
0
at Resonance
am
bn
m
k=1 (i i0,k )
n
k=1 (i i,k )
a m im
b n in
m
k=1 ( 0,k )
n
k=1 ( ,k )
(10.32)
For a model with N-degrees-of-freedom the frequency response function contains N-pole pairs. The order of nominator, m, and denominator, n, polynomial has not to be different. Iglesias (2000) shows that both can be chosen as
n = m = 2N. This leads to the following simplified equation:
a2N
H() =
b2N
2N
k=1
0,k
,k
(10.33)
|,k |
2
Dk =
Re(,k )
|,k |
(10.34)
(10.35)
The practical difficulties in applying this approach are essentially the calculation of appropriate curve fitting parameters ak and bk . Fortunately the problem is almost identical with the design of a digital infinite-impulse-response
204
filter (IIR-filter) based on a given filter frequency response function. The obtained parameters ak and bk are in this case called filter coefficients and could
be directly used to filter a signal in the time domain. For the filter design a
wide number of software packages are available.
In this research a MATLAB script developed by Iglesias (2000) for the
parameter estimation based on an article of Richardson & Formenti (1982) is
used.
10.4.2
1
|F F T (r1st )|
ln
2L |F F T (r2nd )|
(10.36)
Pulse frequency:
3000 Hz
Pulse frequency:
4500 Hz
Pulse frequency:
6000 Hz
1st arrival
0.8 0.5
0.5
2nd arrival
3rd arrival
1.5
t in s
205
2.5
3.5
3
x 10
|F F T (r1st )|
1
ln
4L |F F T (r3rd )|
(10.37)
z Vs
2
(10.38)
The results are summarized in Table 10.7. The calculated damping ratios are
between 6.3 and 9.3 %.
Amplitude ratio at the frequencies with the highest energy
Because, in opposite to a SCPT, the frequency content of the BE-excitation
can be varied easily by changing the length, respectively the frequency, of the
206
10
8
6
z = 2.0352 10
DS = 7.45 %
s/m,
2
0
0
1000
5000
6000
10
8
6
3
z = 2.3415 10
D = 8.57 %
s/m,
2
1st to 2nd arrival
1st to 3rd arrival
0
0
2000
4000
6000
Frequency [Hz]
8000
20
15
10
z = 2.1814 10
D = 7.99 %
s/m,
0
0
2000
4000
6000
Frequency [Hz]
8000
Parameter
Unit
3000 Hz
fmean
S
z
Ds
[Hz]
[m1 ]
[s/m]
[%]
3074
5.13575
2.0352 103
7.45
2922
3.89520
2.0041 103
7.34
4500 Hz
fmean
S
z
Ds
[Hz]
[m1 ]
[s/m]
[%]
4427
7.25377
1.7205 103
6.30
3830
5.91632
2.3415 103
8.57
6000 Hz
fmean
S
z
Ds
[Hz]
[m1 ]
[s/m]
[%]
4953
8.18255
2.5276 103
9.25
4166
6.96369
2.1814 103
7.99
207
Table 10.7: Parameters of the spectral ratio slope method applied to bender
element time histories showing multiple arrival; data based on the signals of
Figure 10.18; fmean : mean dominant frequency, S : attenuation coefficient, z:
rise of S with frequency, Ds : damping ratio.
sinusoidal pulse, the spectral ratio slope method can also be applied in another
way.
The attenuation coefficient at the mean dominant frequency, fmean , for each
excitation frequency is drawn versus the frequency in a separate diagram, as
in Figure 10.22. It is then again possible to fit a line to these points and to
calculate the damping ratio based on the slope z of this line.
A damping ratio of 5.9 % is obtained if first and second arrival are compared.
First and third arrival lead to a damping ratio of 8.9 %. The deviation of these
values to the mean values of Table 10.7 for each considered travel path is slightly
higher than the deviation of the single values in the table from its mean values.
A particular reason for that cannot be given yet.
Correction for the boundaries
Until now it is assumed that during the refections no loss of energy appears.
The SRS method would be able to compensate a frequency independent change
of amplitudes caused by a reflection but no change depending on the frequency.
A frequency independent amplitude decay appears in terms of an S -f diagram
as a parallel upward shift of the curve. Only a frequency dependent decay
changes the slope z and as a consequence Ds .
The loss at the caps can be considered by the introduction of a reflection
208
10
8
6
z = 2.4163 103 s/m,
DS = 8.85 %
4
2
0
2500
3000
5000
5500
factor c() in the Equations (10.36) and (10.37). This means for the first and
the second arrival:
=
=
F F T (r1st ) c()2
1
ln
2L
|F F T (r2nd )|
|F F T (r1st )|
1
1
ln
+
ln c()2
2L |F F T (r2nd )| 2L
(10.39)
=
=
F F T (r1st ) c()4
1
ln
4L
|F F T (r3rd )|
1
|F F T (r1st )|
1
ln
+
ln c()4
4L |F F T (r3rd )| 4L
(10.40)
The exponent of c() indicates the number of reflections the signal has to
undergo. Two reflections are considered for the travel path between 1st and
2nd arrival and four reflections between 1st to 3rd arrival.
The simplest way to model the boundary effects is to assume at the interface
soil to metal cap a two media boundary as shown in Figure 10.23. The factor c
can be expressed based on a constant cross section of both media. The equation
follows Clough & Penzien (1993) and Fratta & Santamarina (1996):
c=
ISample
ICap 1
I
1 + Sample
ICap
(10.41)
209
s1
Loading
Piston
s2 = IFFT(FFT(s1) c())
Sample
c1
Top/Bottom
Cap
Base Plate
The impedances I can be calculated by the material density and the shear
wave velocity Vs :
I = Vs
(10.42)
According to the equation, the reflection factor c is not depending on frequency and therefore has no influence on the damping calculation using the
spectral ratio slope technique. For the used testing device and the tested sample |c| = 0.9586. This value is based on the material parameters in Table 10.8.
Parameter
Density,
S-wave velocity, Vs
Impedance, I
Reflection coefficient, |c|
Unit
[kg/m3 ]
[m/s]
[kg/(m2 s)]
[-]
Sample
Caps
2342
8009
230
3180
5.3866 105 2.5470 107
0.9586
Table 10.8: Calculation of the reflection coefficient c for the case of the silt
sample Waremme B1/3.
Fratta & Santamarina (1996) provide also an equation for a three media
boundary. The base plate of the triaxial cell and the loading piston, as shown in
in Figure 10.24, could be considered as the third medium. In this way the effects
of reflection and refraction on the interface cap to base plate respectively piston
could be taken into account. However, there are some difficulties in determining
a suitable impedance for the third medium. Fratta & Santamarina (1996)
found, based on their sample and apparatus parameters, a reflection coefficient
c() depending on frequency in a range below 1000 Hz. At higher frequencies
c() is nearly constant.
210
10.4.3
The direct application of the spectral ratio slope method or alternatively the attenuation coefficient method on the setup of two bender elements in a triaxial
cell sample is not possible because there is no sufficient information on the
transfer characteristics soil - BE - output voltage, respectively driving voltage
- BE - soil. This transfer characteristic might be expressed as a transfer function (f), potentially highly frequency dependent. (f) incorporates effects of
different resonance frequencies of the elements, different penetration depths of
the elements and the coupling of the element tip and the specimen.
However, there is a theoretical approach to eliminate the influence of (f)
if two measurements with opposite wave propagation directions are combined.
The realization of this idea needs two receiver elements and two transmitter
elements installed on one sample.
55
10
10
7
35
48
35
76
Figure 10.25: Frame with mounted sender (S1/S2) and receiver elements
(R1/R2), unit: mm, Wang et al. (2004).
211
A1
A2
x1
x2
x2 x 1
Setup A
Setup B
Receiver 1
FFT(r1A)
Sender 1
FFT(r1A)1
BE
z1A
BE
Receiver 1
FFT(r1B)
BE
Sender 1
FFT(r1B)1
z2A
Receiver 2
FFT(r2A)
BE
BE
z1B
Sender 2
BE
FFT(r2A)2
Receiver 2
FFT(r2B)
BE
Soil Specimen
FFT(r2B)2
Sender 2
z2B
BE
Soil Specimen
F F T (r1A ) 1
F F T (r2A ) 2
z1A
z2A
z2B
z1B
(10.43)
z2A z1A
F F T (r2B ) 2
F F T (r1B ) 1
(10.44)
z1B z2B
F F T (r1A )
F F T (r2A )
+ ln
1
2
+ ln
z1A
z2A
(10.45)
(f ) (z1B z2B ) = ln
F F T (r2B )
F F T (r1B )
+ ln
2
1
+ ln
z2B
z1B
(10.46)
and
212
F F T (r1A )
F F T (r2A )
F F T (r2B )
F F T (r1B )
z1A
z2A
z2B
z1B
(10.47)
213
10.5
Measurements of damping ratio on samples taken at the test sites based on the resonant method
10.5.1
Introduction
The bender element test is performed on samples from the site in Waremme
and Sint-Katelijne-Waver. The specimens are subjected to different isotropic
stress levels by an external water pressure. The back pressure applied to the
sample is set to 100 kPa. The external isotropic pressure is therefore chosen
always 100 kPa above the value of the target effective stress. The absorption of
water by the sample during consolidation is measured and taken into account
for the calculation of the sample mass. The shear wave velocity is determined
by the analysis of the wave arrival at each stress level. The properties of the
free top cap are: height 19 mm, mass 50.74 g and diameter 50 mm.
The damping ratio is calculated based on half-power bandwidth and circlefit method.
10.5.2
Site in Waremme
Altogether three samples taken from the drilling B1 in Waremme are tested.
Specimen S1 is a reconstituted sample; S2 and S3 are undisturbed samples. All
samples have a length of 29 to 30 mm. The 300 kPa consolidation stage of the
sample S3 is repeated with air as confining medium to study its influence on
the response spectrum.
The Figures 10.27 to 10.29 present the measured response spectra, the Tables 10.9 and 10.10 the predicted and measured natural frequencies including
the physical parameters of the samples.
The theoretical estimation of the natural frequencies of the first mode is
done with an added mass coefficient of CA = 0.48, which corresponds to a
length to diameter ratio of 0.96, and as well with no added mass (CA = 0) to
allow comparisons. The coefficient is calculated by means of Equation (10.27).
As expected, the frequencies including the added mass component are lower
than the values without the additional mass. For the same CA , the results from
the Rayleigh expression show a good agreement with the FEM-values. The
frequencies calculated by the Timoshenko equation are always higher.
The test results show that the natural frequencies and the shear wave velocities increase with an increasing cell pressure. The amplitude ratios of the
observed peaks decrease with increasing stress level.
The predicted frequencies are in all cases higher than the measured ones.
The closest agreement to the experimental results is achieved at CA = 0.48 by
the Rayleigh equation and the FEM. A frequency comparison for the case of
sample S2 can be seen in Figure 10.30
The damping ratio versus effective confinement stress is shown in Figure
10.31. D is in the range of 7.2 to 9.9 % for S1, between 4.7 and 7.3 % for S2
214
and between 2.3 and 7.8 % for S3. The deviation between the values from the
half-power method and the circle-fit method are in nearly all cases less than
1 %. A dependence of the damping ratio on the confining stress could be not
noticed in the stress range applied. The damping values are high in comparison
with the resonant column and free torsion pendulm tests.
Example: half-power bandwidth and circle-fit method
Half-power and circle-fit method are used to obtain the damping ratio. The
Figures 10.32 shows an example of the half-power bandwidth technique applied
on the frequency response curve of the undisturbed silt sample S2 at 100 kPa
confining pressure. The calculation of the damping ratio D is based on Equation
(10.29).
Figure 10.33 gives insides in the circle-fit technique based on the same response spectrum. The following algorithm is used:
Drawing of the Nyquist plot from the experimental frequency response
curve.
Manual selection of some points of the response curve close to the resonance peak. These data points are indicated in the Nyquist plot of Figure
10.33 as data points used for curve fitting.
Based on these selected points a circle is fitted.
The center point of this circle is the origin of the angles . The data
point in the Nyquist plot with the largest distance to its adjoining points
is the point of the natural frequency.
Based of the selected range of data points for the curve fitting many
combinations of - and -values are possible. They are all evaluated by
Equation (10.30) and the damping ratio is obtained.
All D-values are finally averaged and this average value provides the result
of the test.
Influence of the confining medium
For the evaluation of the influence of the water pressure to the response curves
in the experimental bender element test, the manner of applying the confining
stress is changed. The triaxial cell is connected to a source of controllable air
pressure and the 300 kPa consolidation stage of sample S3 is repeated using
air pressure instead of the usual water pressure.
The response curve of this stage, given in Figure 10.34, has a similar shape
as this found using water pressure. The resonant peak of the first mode is
shifted towards higher frequencies but not as much as predicted by analytical
and FEM calculations.
215
10.5.3
Site in Sint-Katelijne-Waver
Two Boom clay samples from the drilling B1I in Sint-Katelijne-Waver are tested
by the bender element resonant technique. The length of sample C4 is about
30 mm and sample C5 is 60 mm long. The additional sample C5 is selected to
prove if the deviated equations for the resonance frequency can properly predict
216
10.6
217
and flexible.
In the following methods to calculate the damping ratio are treated. An
approach to analyze multiple refections by the spectral ratio slope method is
applied on a sample of the site in Waremme. Damping ratios calculated by
this technique for the silt material are in the same range as obtained by the
BE resonant test.
A self correcting method is described to eliminate the influence of the transfer characteristic of the bender elements. This enables the application of the
spectral ratio slope method.
The determination of the damping ratio focuses on the bender element
resonant method. At first the behavior of the sample-top-cap system is modeled
analytically, using the Timoshenko beam theory and Rayleighs approximation,
and numerically, using the FE program ABAQUS. This provides information
on the natural frequency of the first flexural mode to be used in the modal test.
The influence of the confinement medium water is considered by an added mass
component and experimentally validated.
The first flexural mode is identified in experimental response curves by
means of the natural frequency prediction. The response curves are analyzed
using circle-fit and half-power method to obtain the damping ratio. Undisturbed samples from Waremme and Sint-Katelijne-Wave are tested.
Two basic problems appear looking at the test results. The first is the
deviation of the calculated and the measured resonant frequencies. The assumption of a beam with top-mass at its one end, even if shear effects are taken
into account, might be a relatively crude way to describe the very short, about
3 cm long, sample and the top cap. The finite element model can partly overcome this geometrical problem, but also in this case, influences of the rubber
membrane and the not fully rigid mounting of the bottom cap are not taken
into account. Furthermore the consideration of the water in the cell by an
additional mass might not be sufficiently correct. If the viscosity of the water
is taken into account, a frictional (viscous) drag force, proportional to the velocity of the moving cylinder, has to be introduced (Wilson (1984)). The used
coefficient CA for the added mass depends, besides the cylinder geometry, also
on the Reynolds number and the cylinder roughness. A more detailed investigation could help to find a more exact theoretical solution. On the other hand
the prediction of the natural frequency is only used to find the right resonant
peak among several, so the calculations are already sufficient for a successful
selection.
The second problem concerns the difference between the results of the BEresonant test and the other laboratory tests. Possible explanations for this
phenomenon can be external factors as the supply cables to the top-cap, the
rubber membrane and the confining medium affecting slightly the oscillation.
Another factor could be a possible frequency dependence of the damping ratio.
Remarks on the frequency dependence of the damping ratio A possible reason of the high damping values obtained by the bender element resonant
218
technique might be the different range of working frequency of the tests. Cyclic
tests are usually performed at frequencies below 2 Hz; free torsion pendulum
test show commonly natural frequencies of sample and device below 10 Hz.
Even resonant column tests devices reach, dependent on the sample geometry
seldom resonance frequencies higher than 250 Hz, in the most cases much lower.
In contrast to that, samples subjected to a BE-resonant test have natural frequencies, crucial for the damping determination, in the range between 400 to
900 Hz.
A frequency dependence of the damping ratio could possibly explain the
deviation in the results. Unfortunately, this question is not completely solved
yet. Dobry & Vucetic (1987) write in their summary, given in table Table 2.2,
on influence factors of the damping ratio, that with an increasing strain rate,
i.e. frequency, the damping ratio stays or may increase.
Mok et al. (1988) suggested based on the work of others that, for dry
rocks and dry soils, D is independent of frequency. On the other hand, watersaturated rocks, sands and silts would show a definite dependency of D on frequency. For saturated clays the authors conclude a relatively small frequency
dependency.
Several authors compare results from cyclic torsional shear tests and resonant column tests. Zavoral & Campanella (1994) are doing this for a marine
silt to clay deposit and found no frequency dependence. The cyclic torsional
apparatus was driven at frequencies between 0.01 and 1 Hz; the resonance was
reached in the range of 10 to 30 Hz. Bolton & Wilson (1989) come to the same
finding for a dry sand. The cyclic testing frequency was fixed at 0.001 Hz,
the authors call it a pseudo-static test; the resonant conditions were obtained
between 45 and 95 Hz.
Lin et al. (1996), tested two types of sand at various frequencies in the range
from almost static to 20 Hz using again a cyclic torsional device. They found
an increase of damping ratio from 2.6 to 6.4 % for the first sand and between
1.9 to 6.2 % for the other sand. These tests were performed at a shear strain of
4 103 %. Remarkably no significant change in the G-modulus was observed.
A few in situ tests give also some indication of the influence of the frequency
on the damping ratio. Hall & Bodare (2000) describe a cross-hole test performed in layers of mud and clay where the attenuation coefficient method is used
to analyze the attenuation. In the range between 20 and 30 Hz the damping
ratio is found to be frequency independent. Mok et al. (1988) came to the same
result for frequencies between 170 and 300 Hz and a site dominated by a clay.
Kudo & Shima (1981) confirm an almost frequency independent behavior for
an even larger range from 5 to 100 Hz based on a down-hole set-up in sand.
Mok et al. (1988) investigated the damping of p-waves in dry sand based on
a miniaturized cross-hole set-up and the attenuation coefficient method. They
found in the range of 2200 to 3000 Hz an almost constant damping ratio.
A mathematical model including a frequency dependent damping can be
obtained by Biots theory concerning the interaction of the pore water and the
soil fabric. If the soil is assumed to be a saturated two-phase medium the
219
damping ratio of the shear wave, following Biots theory, has to be frequency
dependent. The damping ratio increases form the static case with the frequency
until a local maximum is reached after which the ratio decreases again. Figure
10.42, taken from Stoll (1979), shows the logarithmic decrement as function of
frequency according to Biots theory for the parameters given in Table 10.13.
Beside self-explanatory values, a pore-size parameter (a) and a structure constant (c) are influencing the shape of the curves. The pore-size parameter
(a) depends on the size and shape of the pores and can be estimated by the
mean grain diameter. However, the structure constant (c) is an experimentally
determined parameter not easy to quantify.
Curve No. 1 shows the results of the calculation for a less permeable material. The logarithmic decrement stays constant in the whole observed frequency
range and is equal to the decrement of the frame. All other parameter sets with
higher permeability lead to the described frequency dependent curves.
Further information on Biots theory can be found in the basic papers Biot
(1956a) and Biot (1956b).
As a conclusion there is evidence that comparing free torsion pendulum
tests, resonant column tests and BE tests on the level of damping ratio can
give different results due to frequency dependence of the ratio.
220
x 10
582 Hz
100 kPa
200 kPa
300 kPa
0.8
708 Hz
0.6
807 Hz
0.4
0.2
0
100
200
300
400
500
600
700
800
900
1000
Frequency [Hz]
x 10
636 Hz
100 kPa
200 kPa
300 kPa
1.5
766 Hz
0.5
820 Hz
0
100
200
300
400
500
600
700
800
900
1000
Frequency [Hz]
221
1.2
x 10
100 kPa
200 kPa
300 kPa
631 Hz
1
0.8
735 Hz
0.6
811 Hz
0.4
0.2
0
100
200
300
400
500
600
700
Frequency [Hz]
800
900
1000
1200
Timoshenko Equation
1000
FEM
800 Rayleigh
Equation
Experiment
600
400
200
0
100
200
300
222
10
HalfPower Method
CircleFit Method
S1
9
S1
8
S1
S3
7
S3
S2
S2
Sampling Depths:
S1, S2 > 4.04.5 m,
S3 > 5.05.5 m
S2
5
4
3
2
1
0
0
100
200
300
Stress, 1, 2, 3 [kPa]
400
500
Figure 10.31: Damping ratios obtained from the samples S1, S2 and S3 in
comparison with resonant column test and free torsion pendulum test on other
samples of the same material (reference shear strain 103 %).
2.5
x 10
f = 636 Hz
m
max
A = Amax / 20.5
1.5
=(f f ) / (2 f ) = 6.2 %
HPM
f1 = 593 Hz
0.5
f = 672 Hz
2
0
100
200
300
400
500
600
700
800
900
1000
Frequency [Hz]
223
x 10
Nyquist plot of
the respose spectrum
Imaginary Component
DCFM = 5.9 %
0
fitted circle
5
center point
10
15
20
1.5
resonance
frequency
0.5
0.5
Real Component
1.5
5
x 10
Figure 10.33: Circle-fit technique on the Nyquist plot of the response spectrum
of Figure 10.32.
x 10
811 Hz
3
930 Hz
2
0
0
500
1000
Frequency [Hz]
1500
2000
224
Waremme
100 kPa
33 kPa
200 Pa
50 kPa
25 kPa
400 kPa
62 kPa
400 kPa
100 kPa
200 kPa
77 kPa
50 kPa
Depth [m]
14
0
50
100
600 kPa
400 kPa
100 kPa
50 kPa
100 kPa
150
200 kPa
141 kPa
12
200 kPa
126 kPa
10
200
250
300
350
225
Waremme
Mechanical Hammer
Sledge Hammer
2
62 kPa
62 kPa
Depth [m]
77 kPa
BEresonant method
10
141 kPa
12
14
0
50 kPa
100 kPa
10
226
x 10
397 Hz
50 kPa
100 kPa
200 kPa
400 kPa
460 Hz
6
5
609 Hz
4
3
2
735 Hz
1
0
100
200
300
400
500
600
700
Frequency [Hz]
800
900
1000
x 10
50 kPa
100 kPa
200 kPa
400 kPa
0
100
200
300
400
500
600
700
Frequency [Hz]
800
900
1000
227
16
HalfPower Method
CircleFit Method
14
12
Sampling Depth:
C4 > 4.04.5 m
10
8
6
4
2
0
0
100
200
300
Stress, 1, 2, 3 [kPa]
400
500
228
SintKatelijneWaver
400 kPa
200 kPa
50 kPa
100 kPa
400 kPa
100 kPa
121 kPa
200 kPa
BEresonant method
Depth [m]
121 kPa
9
0
50
100
150
200
250
300
229
SintKatelijneWaver
confining stress stages used in laboratory tests
0
Mechanical Hammer
Sledge Hammer
1
00
00
,2
200 kPa
50 kPa
400 kPa
BEresonant method
Free Torsion
Pendulum Test
100 kPa
Depth [m]
,1
50
a
kP
4
0,
9
0
10
15
230
Sample Description:
Sample index
Site
Description
Sampling
Depth
[m]
S1
Waremme
Silt
reconst. sample
4.0-4.5
S2
Waremme
Silt
undist. sample
4.0-4.5
29.0
50.0
25.0
101.53
1783.10
29.5
50.0
25.5
119.49
2062.90
100
11.4
113.0
216.9
200
10.9
112.5
224.5
300
10.2
111.8
246.3
100
-3.2
116.3
197.7
200
-5.6
113.9
246.1
300
-7.0
112.5
271.3
Resonant frequencies:
CA = 0
Rayleigh Eq., fRay
[Hz]
Timosh. Eq., fTBT
[Hz]
FEM, fFEM
[Hz]
888
962
875
918
995
904
1005
1090
989
799
864
787
988
1070
972
1084
1175
1066
CA = 0 .48
Rayleigh Eq., fRay
Timosh. Eq., fTBT
FEM, fFEM
[Hz]
[Hz]
[Hz]
773
843
763
799
871
788
875
954
862
696
757
686
860
937
847
943
1028
929
[Hz]
582
708
807
636
766
820
Damping ratio:
Half-Power, DHPM
Circle-fit, DCFM
[%]
[%]
7.9
7.2
9.1
9.9
8.0
8.9
6.2
5.9
7.3
4.7
5.4
4.9
Table 10.9: Test data for sample S1 (reconstituted silt) and S2 (undisturbed
silt).
Sample Description:
Sample index
Site
Description
Sampling
Depth
Initial test parameters:
Length, L
Diameters, d
Distance BE tips
Mass, mT
Unit mass,
231
[m]
S3
Waremme
Silt
undist. sample
5.0-5.5
[mm]
[mm]
[mm]
[g]
[kg/m3 ]
30.0
50.0
26.0
116.54
1978.44
Consolidation parameters:
Effective stress, peff
[kPa]
100
200
300
Water absorption
Effective mass, mT,eff
S-wave velocity, Vs
[g]
[g]
[m/s]
1.3
117.8
195.1
0.7
117.2
230.6
0.4
116.9
255.0
300
air
0.4
116.9
257.0
Resonant frequencies:
CA = 0
Rayleigh Eq., fRay
Timosh. Eq., fTBT
[Hz]
[Hz]
774
837
913
988
1009
1092
1017
1100
CA = 0 .48
Rayleigh Eq., fRay
Timosh. Eq., fTBT
[Hz]
[Hz]
674
766
796
904
879
999
886
1007
[Hz]
631
735
811
930
Damping ratio:
Half-Power, DHPM
Circle-fit, DCFM
[%]
[%]
7.8
7.3
6.8
6.2
(10.0)
2.3
(17.5)
2.1
232
Sample Description:
Sample index
Site
Description
Sampling
Depth
[m]
C4
Sint-Katelijne-Waver
Boom-Clay
undist. sample
4.0-4.5
29.8
50.0
25.8
114.60
1961.86
Consolidation parameters:
Effective stress, peff
[kPa]
Water absorption
[g]
Effective mass, mT,eff
[g]
S-wave velocity, Vs
[m/s]
50
5.7
120.3
157.0
100
8.6
123.2
178.2
200
8.5
123.1
253.5
400
8.1
122.7
299.6
Resonant frequencies:
CA = 0 .48
Rayleigh Eq., fRay
[Hz]
553
634
902
1064
[Hz]
397
460
609
735
Damping ratio:
Half-Power, DHPM
Circle-fit, DCFM
[%]
[%]
15.7
10.7
12.1
7.3
14.9
12.8
13.5
13.4
Sample Description:
Sample index
Site
Description
Sampling
Depth
[m]
233
C5
Sint-Katelijne-Waver
Boom-Clay
undist. sample
4.0-4.5
59.5
50.0
55.5
221.62
1885.65
Consolidation parameters:
Effective stress, peff
[kPa]
Water absorption
[g]
Effective mass, mT,eff
[g]
S-wave velocity, Vs
[m/s]
50
8.1
229.7
175.4
100
7.3
228.9
205.9
200
5.5
227.1
232.9
400
2.8
224.4
290.0
Resonant frequencies:
CA = 0 .71
Rayleigh Eq., fRay
[Hz]
255
299
337
418
[Hz]
Damping ratio:
Half-Power, DHPM
Circle-fit, DCFM
[%]
[%]
234
Log. Decrement []
Curve No.
1
2
3
4
5
101
102
1
102
Frequency [Hz]
10
103
104
Material property
Mass density of grains, r
Mass density of fluid, f
Viscosity of fluid,
Shear modulus of frame, G
Structure constant, c
Intrinsic permeability, ki
Unit
kg/m3
kg/m3
Pa s
MPa
m2
Permeability, k = ki f g/
with g = 9.81m/s2
Pore-size parameter, a
m/s
Porosity,
Log. decrement of frame
for shearing motion,
Damping ratio, D = /(2)
2.5
1014
2.45
107
3.8
106
0.49
0.08
1.27
No. 1
0.03
0.48
0.03
0.48
No. 5
5
1011
4.91
104
2.8
105
0.36
0.02
0.32
Table 10.13: Material properties of the Biot theory example, Stoll (1979).
Chapter 11
236
237
certain mode shapes. Because the first flexural mode is used to obtain the
damping ratio, its frequency is predicted by means of two analytical formulations (Timoshenko beam theory, Rayleighs method) and a FE calculation
(ABAQUS). The influence of the confining medium is taken into account by
the introduction of an added mass component. A reasonable approximation of
measured and predicted natural frequencies is obtained.
Results based on bender element resonant method are gathered on samples
of a silt (5 10 %) and a clay (10 16 %) material. These results are repeatable but higher than the results from resonant column and free torsion pendulum
tests (1 4 %) performed for comparative purposes. There is some evidence
that this deviation is caused by the higher frequency used in the BE-test in
comparison with the two other laboratory tests.
Another interpretation technique uses the property of shear wave pulses reflected on the end caps of the triaxial device. The SRS method can be applied
based on the spectral ratio between reflections of different order. Damping
ratios calculated in this approach are between 7 and 9 % for the silt material.
This is the same range as obtained by the BE resonant test. Unfortunately multiple reflections could be not identified in samples of the clay material probably
because of the higher attenuation in the clay.
A SRS based method to eliminate the influence of the transfer characteristic
of the bender elements by means of a combined analysis of different travel paths
is discussed.
Practical relevance The SCPT and BE test methods for the determination
of the damping ratio need still a further refinement until they can be recommended for the practical application. Especially the scattering seen in the ACM
and SRS method should be reduced. However, by taking in consideration that
there are nearly no other, at least no cost-equivalent, field tests available to
determine the damping ratio, it seems obvious to spend additional effort in
this technique.
Reasons for the high damping values of the BE-resonant method should
also be closer investigated.
238
Part IV
Appendix
239
Appendix A
242
A.1
Laboratory tests
243
time, the method is able to provide results of high accuracy and reproducibility.
Since the technique is in use for several years a lot of experience is available
and its application is widely spread.
The comparability to the resonant and the cyclic tests suffers sometimes
under the unknown level of shear strain reached. But it is assumed that the
shear strain is always lower during a bender element test than in the other
laboratory tests.
Damping ratio determination In comparison to the time arrival techniques, the determination of the damping ratio by the described resonant method
is not a commercial testing technique. It is only used in a research environment.
In the frame of this research the BE-resonant method and, to a smaller
extend, also the multiple reflection approach is applied. Both techniques lead
to significantly (1 to 10 %) higher damping ratios than the resonant techniques
(resonant column and free torsion pendulum test). The results vary also with
the used analyzing technique, in case of the BE-resonant method, and the used
excitation frequency, in case of the multiple reflection approach. This variation
is in the range of +/ 0.25 to 2.5 % damping ratio. Based on these results the
accuracy of this techniques is considered as comparatively small.
The costs of the tests are similar to the time arrival method. It has to be
ensured only that the top cap of the triaxial apparatus is light in weight and
can move freely.
Ultrasonic tests / shaking table
These tests use ultrasonic acoustic waves to determine the p-wave velocity only.
Since they are based on the same time arrival techniques as the bender element
test, the accuracy of this method is considered as high. Application to obtain
the damping ratio of soil materials are not known.
Shaking table set-ups are always very specifically orientated to a certain
problem. They are therefore rather expensive and the accuracy depends on the
chosen configuration.
A.2
Field tests
The principal advantage of field tests is that they avoid the need to obtain
undisturbed samples by drilling. This reduces on one side the testing costs
and on the other side effects caused by sample disturbances. Furthermore field
tests can be performed rapidly in comparison to laboratory tests where it needs
some consolidation time to adjust the in situ stress conditions.
Surface methods
Surface techniques need no costly equipment to penetrate the ground and are
therefore the most inexpensive in situ testing techniques. Seismic reflection
and refraction tests but also the SASW tests belong to this group.
244
The reflection and refraction test provide, dependent on the actual soil stratification, wave velocities of medium to high accuracy. Approaches to determine
also the damping ratio are not common.
The accuracy of the SASW technique depends on the complexity of the soil
layering but can be considered as medium to high. Experiences concerning the
accuracy of the obtained damping ratio, if a combined Vs and D inversion is
used, are very limited in literature, but promising.
Penetrating methods
The SCPT belongs to this group . Its price is situated between a purely surface
technique and a borehole method. The accuracy concerning the s-wave velocity
has shown to be very high coequal to a cross- or down-hole test.
However, the determination of the damping ratio by spectral ratio slope
or attenuation coefficient method is much less accurate and produces quite
scattered profiles.
Borehole methods
These are the classical seismic methods to obtain Vs and Vp , the cross-, downand up-hole methods. The determination of the wave arrival is no questionable
issue anymore. Therefore the accuracy concerning the stiffness parameters is
high.
The determination of the damping ratio is also based on spectral ratio slope
or attenuation coefficient method. Published results indicate a slightly better
accuracy than found with SCPT. But also here the experiences are limited.
The borehole tests are quite expensive since the sender and / or the receiver
are to be lowered in a cased borehole. In case of a cross-hole test at least two
holes are necessary, for a damping ratio determination even three holes are
needed.
A.3
245
Currency of
determination of
small
large
Applicability at
G / Gmax /
+
+
++
++
++
+++
++++
+
+++
+
++++
++
++
++
++
++
+
+
?
+++
++
+++
+
+ (?)
+
+
+
++
+
+
-
Costs
+++
++
+
+
+
+
+++
++
+++
+
++++
+++
++
Accuracy / Reliability
G
D
strains
strains
V s / Vp
Laboratory tests
Resonant column test (cylinder / hollow-cylinder)
f
f
g/g
l/m
+++
Free torsion pendulum test
i
i
m
m
++
BE - time arrival method
f
g
+++
BE - resonant method
r
r
g
+
BE - multiple reflection approach
r
r
g
+++
BE - combination of different travel paths
r
r
g
+++
Cyclic triaxial test
f
f
l
m
+++
Cyclic simple shear test
i
i
l
m
++
Cyclic torsional shear test
f
f
g
m
+++
i
g
+++
Ultrasonic tests (only Vp )
Shaking table
r
r
g
++
In situ tests
f
g
++
SASW - Vs profile
SASW - combined inversion for Vs and D profile
r
r
g
++
SCPT - time arrival
f
g
+++
SCPT - spectral ratio slope / attenuation coefficient method
r
g
SCPT - hysteresis approach
r
r
g
+
Down- and up-hole tests
i
r
g
+++
Cross-hole tests
i
r
g
+++
Reflection-/refraction tests
f
g
++
Free oscillation tests, enforced oscillations
r
r
g
m
+
Dynamically loaded plates
f
l
g
+
In situ shear tests
r
r
l
g
+
Geotomography
i
g
++
Dilatometer test (DMT)
i
g
++
Pressuremeter test (PMT)
i
g
++
Standard penetration test (SPT) with correlation
f
l
+
Cone penetration test (CPT) with correlation
f
l
+
Legend:
Currency: f frequently, i infrequent, r under research conditions
Shear strain level: small strain < 103 %, large strain > 103 %
Applicability: g good, m medium, l less, - not suitable
Accuracy/Costs: + low, ++ medium, +++ high, ++++ very high, - not applicable, ? unknown
Testing method
246
APPENDIX A. ACCURACY AND COSTS OF SOIL TESTS
Testing method
Laboratory tests
Resonant column test
Laboratory
Costs
(excl.
VAT)
a)
b)
c
c
2300 e
1850 e
650 e
a)
250 e
b)
a)
2100 e
b)
550 e
Field tests
SCPT
Cross-hole test
SASW test
247
a)
b)
a)
b)
n
n
680 e
900 e
100 to
250 e
60 e
profile to 12 m depth, seismic registrations every meter, application of a pand a s-wave source, p- and s-wave velocity determination, including one source
borehole and two receiver holes
same as a)
3000 e
6000 e
one profile
one day of site work and one day of signal
processing and reporting
200 e
1500 to
1800 e
248
Bibliography
Airey, D. W. & Wood, D. M., 1987. An evaluation of direct simple shear tests
on clay. Geotechnique, 37(1):2535.
Areias, L. & Haegeman, W., 2001. In Situ Tests Retie-Waremme-Lincent.
Interim Report IV.1.16.3 (19/09/2001), Ghent University.
Areias, L., Van Impe, W. F. & Haegeman, W., 1999. Variation of shear wave
energy with coupling stress in the SCPT method. European Journal of Environmental and Engineering Geophysics, 4:8795.
Areias, L., Haegeman, W. & Van Impe, W. F., 2004. New de-coupled shear
wave source for the SCPT test. In Proceedings of the 2nd International
Conference on Geotechnical Site Characterization, volume 1, pages 643645.
Porto.
Arroyo, M., 2001. Pulse tests in soil samples. Ph.D. thesis, University of
Bristol.
Arulnathan, R., Ross, W. B. & Riemer, M. F., 1998. Analysis of Bender
Element Tests. Geotechnical Testing Journal, 21(2):120131.
Best, R., 1991. Digitale Mewertverarbeitung. Oldenbourg, Munich, Vienna.
Biot, M., 1956a. Theory of propagation of elastic waves in a fluid saturated
porous solid. I. Low frequency range. Journal of the Acoustical Society of
America, 28(2):168178.
Biot, M., 1956b. Theory of propagation of elastic waves in a fluid saturated
porous solid. II. High frequency range. Journal of the Acoustical Society of
America, 28(2):179191.
Blewett, J., Blewett, I. J. & Woodward, P. K., 1999. Measurement of shearwave velocity using phase-sensitive detection techniques. Canadian Geotechnical Journal, 36:934939.
Bolton, M. D. & Wilson, M. R., 1989. An experimental and theoretical comparison between static and dynamic torsional soil tests. Geotechnique, 39(4):585
599.
249
250
BIBLIOGRAPHY
BIBLIOGRAPHY
251
252
BIBLIOGRAPHY
BIBLIOGRAPHY
253
254
BIBLIOGRAPHY
Richart, F. E., Hall, J. R. & Woods, R. D., 1970. Vibration of soils and
foundations. Prentice-Hall, Englewood Cliffs, N. J.
Rossit, C. A. & Laura, P. A. A., 2001. Transverse, normal modes of vibration
of a cantilever Timoshenko beam with a mass elastically mounted at the free
end. Journal of the Acoustical Society of America, 6(110):28372840.
Sanches-Salinero, I., Roesset, J. M. & Stokoe, K. H. II, 1986. Analytical studies
of body wave propagation and attenuation. Report GR 86-15, University of
Texas, Austin.
Seco, F., Martn, J. M., Jimenez, A., Pos, J. L., Calder
on, L. & Ceres, R.,
2002. PCDISP: A tool for the simulation of wave propagation in cylindrical
waveguides. In 9th International Congress on Sound and Vibration. Orlando,
Florida.
Seed, H. B. & Idriss, I. M., 1970. Soil moduli and damping for dynamic response
analysis. Report EERC 70-10, UCB.
Stewart, W. P., 1992. Insitu measurements of dynamic soil properties with
emphasis on damping. Ph.D. thesis, University of British Columbia, Department of Civil Engineering, Vancouver.
Stoll, R. D., 1979. Experimental studies of attenuation in sediments. Journal
of the Acoustical Society of America, 66(4):11521160.
Storrer, O., Van Den Broeck, M. & Van Impe, W. F., 1986. A numerical model
analysing free torsion pendulum results. In 2nd International Symposium
on Numerical Models in Geomechanics, pages 265273. Ghent University,
Ghent.
Studer, J. A. & Koller, M. G., 1997. Bodendynamik; Grundlagen, Kennziffern,
Probleme. Springer, Berlin, 2 edition.
Thurston, R. N., 1978. Elastic waves in rods and clad rods. Journal of the
Acoustical Society of America, 64(1):137.
Thurston, R. N., 1992. Elastic Waves in Rods and optical Fibers. Journal of
Sound and Vibration, 159(3):441467.
Udas, A., 1999. Principles of Seismology. Cambridge University Press, Cambridge, 1 edition.
Van Impe, W., 1977. De vrije wringingsproef: een belangrijke bijdrage tot
de studie van het gedrag van funderingen onder invloed van grondtrillingen.
Tijdschrift der Openbare Werken van Belgie, pages 118.
Viggiani, G. & Atkinson, J. H., 1995. Interpretation of Bender Element Tests.
Geotechnique, 45(1):149154.
BIBLIOGRAPHY
255
Wang, Y. H., Yan, W. M. & Lo, K. F., 2004. Laboratory and in-situ measurement of attenuation in soil. In Viana da Fonseca & Mayne, editors, 2nd
Conference on Geotechnical and Geophysical Site Characterization, volume 2,
pages 18831889. Porto.
Weinstein & Clifford, 1979. Programs for Digital Signal Processing. IEEE
Press, New York.
Wilson, J. F., 1984. Dynamics of offshore structures. John Wiley & Sons, New
York.
Winkler, K. & Nur, A., 1979. Pore Fluids and Seismic Attenuation in Rocks.
Geophysical Research Letters, 6(1):14.
Yoon, Y.-W., 1991. Static and dynamic behavior of crushable and non-crushable
sands. Ph.D. thesis, Ghent University, Laboratory of Soil Mechanics.
Yuan, D. & Nazarian, S., 1993. Automated Surface Wave Method: Inversion
Technique. Journal of Geotechnical Engineering, 119(7):11121126.
Zavoral, D. Z. & Campanella, R. G., 1994. Frequency effects on damping/modulus of cohesive soil. In R. J. Ebelhar, V. P. Drnevich & B. L.
Kutter, editors, Dynamic Geotechnical Testing II, ASTM STP 1213, pages
191201. ASTM, Philadelphia.
Zeghal, M. & Elgamal, A.-W., 1994. Analysis of site liquefaction using earthquake records. Journal of Geotechnical Engineering, 120(6):9961017.
Zeghal, M., Elgamal, A.-W., Tang, H. T. & Stepp, J. C., 1995. Lotung Downhole Array. II: Evaluation of Soil Nonlinear Properties. Journal of Geotechnical Engineering, 121(4):363377.