Pelecanos L 2014 PHD Thesis PDF
Pelecanos L 2014 PHD Thesis PDF
Pelecanos L 2014 PHD Thesis PDF
LOIZOS PELECANOS
MEng(Hons) ACGI
ii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
iii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Declaration
The work presented in this thesis was carried out in the Department
of Civil & Environmental Engineering at Imperial College London from the
3rd of October 2009. It was jointly supervised by Dr Stavroula Kontoe and
Professor Lidija Zdravković and was funded by the Engineering and Physical
Sciences Research Council (EPSRC), UK.
This thesis is the result of my own work and any quotation from, or
description of the work of others is acknowledged herein by reference to the
sources, whether published or unpublished. This thesis is not the same as
any that I have submitted for any degree, diploma or other qualification at
any other university. No part of this thesis has been or is being concurrently
submitted for any such degree, diploma or other qualification. This document
is available online at www.imperial.ac.uk, it is less than 100,000 words long,
contains less than 300 figures and less than 20 tables.
The copyright of this thesis rests with the author and is made available
under a Creative Commons Attribution Non-Commercial No Derivatives li-
cence. Researchers are free to copy, distribute or transmit the thesis on the
condition that they attribute it, that they do not use it for commercial pur-
poses and that they do not alter, transform or build upon it. For any reuse
or redistribution, researchers must make clear to others the licence terms of
this work.
Loizos Pelecanos
iv
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Abstract
Many earth dams around the world are located in zones characterised
by moderate to high seismicity. Their seismic stability can be particularly
critical for the safety of the areas in the downstream side and therefore an in
depth understanding of their response during earthquakes is required. This
thesis describes a numerical study related to both the seismic response and
analysis of earth dams using the finite element method.
In the first part of the thesis, the effect of the upstream reservoir hy-
drodynamic pressures on the elastic seismic response of dams is explored.
Firstly, a methodology is proposed in which the reservoir domain is modelled
with finite elements focusing in particular on the accurate prediction of the
hydrodynamic pressures on the upstream dam face. Secondly, a parametric
study of dam-reservoir interaction is carried out to examine the effect of the
reservoir on the seismic response of dams.
The second part of the thesis is concerned with the nonlinear seismic
behaviour of earth dams. Firstly, a well-documented case study, the La
Villita dam in Mexico, is analysed in order to validate the numerical model
and a good agreement is obtained between the recorded and predicted data.
Subsequently, using as a reference the calibrated model, parametric studies
are performed in order to obtain a better insight into the dynamic response
and analysis of earth dams. The latter studies provide a means to assess the
effect of different modelling considerations on the seismic analysis of dams.
v
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Acknowledgements
Firstly and most importantly, I would like to thank God for everything.
Not all the people have the opportunities I had.
Secondly, I would like to express my sincere gratitude to my two research
supervisors and academic mentors, Dr Stavroula Kontoe and Professor Lidija
Zdravković. I am indebted to them, for inspiring me since my undergradu-
ate years, for providing me the opportunity to do research, for introducing
me to the scholar family of Geotechnics and for wholeheartedly welcoming
me to their research group. I need to express my appreciation to them for
everything they have taught me, both academic and non-academic matters;
for their constant availability, for the excellent computing facilities they pro-
vided, for teaching me to be professional, consistent and patient in my work
and for giving me the freedom and support to investigate my own ideas. I
am hugely grateful to them for being extremely patient with me and I need
to apologise for any inconvenience I may have caused to them.
This study was funded by a Departmental Training Grant from the En-
gineering and Physical Sciences Research Council, UK. Their generous fi-
nancial help is highly appreciated and acknowledged. This acknowledgement
is of course extended to my supervisors who made the relevant necessary
arrangements.
My PhD thesis was examined by Professors David M. Potts and Panos
Dakoulas, two eminent academics and experts in their research fields, nu-
merical analysis and dynamics of dams, respectively. I would like to express
my great appreciation to them for their keen interest in my work and all the
kind advice, feedback and directions they provided to me.
I regard it a great honour and privilege to have been part of the Geotech-
vi
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
nics Research Group at Imperial College, considering its long history and
academic reputation. The University has provided exceptional research fa-
cilities and kind access to an enormous amount of library material. I would
like to thank all the members of the staff for all their contribution and es-
pecially Professor David M. Potts for his valuable knowledge and advice on
numerical analysis and for providing access to ICFEP. I also need to thank
my MPhil/PhD transfer examiner, Dr Catherine O’Sullivan for her interest
in my work and her advice on this project and my future career. Of course,
I need to express my acknowledgements to all the previous and current re-
search students with whom we spent the last years together. Special thanks
are due to my officemate Vasilis Avgerinos for our day-to-day cooperation,
for sharing all sorts of problems and for our endless discussions. Likewise,
thanks are also due to my colleague Bo Han for sharing various thoughts
about dams and dynamic analysis. Moreover, I need to thank Miss Evan-
gelia Skiada for her assistance with a number of analyses of Chapter 5 during
her MSc studies the summer of 2012.
The last years have been very enjoyable, due to my close friends in Lon-
don. I need to express my gratefulness to my brother Fotis and my friends
Michalis Trypiniotis, Dr Yiannis Vyrides, Melinos Averkiou and Michael
Yap. Special thanks are due to my friend Mohammed Latheef, for devel-
oping tailor-made scripts for automated run of the analyses presented in
Chapter 5.
I take also the opportunity to thank all the people who offered their
assistance in any way to complete this work; assistance that I am or I am
not aware of.
Finally, I would like to express my deepest gratitude to my beloved family
and the people closest to me, from whom I had to travel away and sacrifice
the invaluable time that could have been spent with them, in order to be
educated: my wife Katerina, my respected parents Panayiotis and Irene and
my brothers Fotis and Angelos. I hope I make you happy. May God Bless
you all.
vii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Contents
Abstract v
Acknowledgements vi
Contents viii
Nomenclature xxxvii
1 INTRODUCTION 1
1.1 Background . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Scope of research . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Layout of thesis . . . . . . . . . . . . . . . . . . . . . . . . . . 3
viii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
ix
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
x
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xi
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
8 CONCLUSION 330
8.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . 330
8.2 Summary of findings . . . . . . . . . . . . . . . . . . . . . . . 331
8.2.1 Numerical modelling of the reservoir domain . . . . . . 331
8.2.2 Dynamic reservoir-dam interaction . . . . . . . . . . . 332
8.2.3 Numerical analysis of La Villita dam . . . . . . . . . . 333
8.2.4 Parametric seismic analysis of La Villita dam . . . . . 335
8.3 Recommendations for further research . . . . . . . . . . . . . 336
8.3.1 Numerical modelling of the reservoir domain . . . . . . 337
8.3.2 Dynamic reservoir-dam interaction . . . . . . . . . . . 337
8.3.3 Numerical analysis of La Villita dam . . . . . . . . . . 339
8.3.4 Parametric seismic analysis of La Villita dam . . . . . 340
8.4 Comments on the methods of analysis . . . . . . . . . . . . . . 341
8.4.1 Discussion on numerical and analytical methods . . . . 342
8.4.2 Recommendations for use of the methods . . . . . . . . 344
REFERENCES 369
xiii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
List of Figures
xiv
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xv
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xvi
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xvii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xviii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xix
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xx
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxi
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxiii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxiv
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxv
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
6.1 Aerial photo of La Villita dam (Google Earth, 2010) with ap-
purtenant facilities: (1) Reservoir, (2) Embankment, (3) Hy-
droelectric power plant, (4) Spillway, (5) Downstream Balsas
river. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 197
6.2 Recorded displacements on some points on the dam body (El-
gamal, 1992). . . . . . . . . . . . . . . . . . . . . . . . . . . . 199
6.3 Upstream-downstream cross-section of La Villita dam with
zones of different materials (After Elgamal (1992)). . . . . . . 200
6.4 Longitudinal cross-section of La Villita dam showing the geo-
logical profile of the site and the shape of the canyon (Elgamal,
1992). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 200
6.5 Plan view of La Villita dam and locations of accelerometers
(Elgamal et al., 1990). . . . . . . . . . . . . . . . . . . . . . . 202
6.6 Available bedrock records in the UD direction for EQ1-EQ3,
EQ5, EQ6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 203
6.7 Available crest records in the UD direction for EQ1-EQ6 . . . 204
6.8 Response spectra (damping, ξ = 5%) of the available acceler-
ation records in the UD direction for EQ1-EQ3 . . . . . . . . 205
6.9 Response spectra (damping, ξ = 5%) of the available acceler-
ation records in the UD direction for EQ4-EQ6 . . . . . . . . 206
6.10 The recorded UD motion on the crest of La Villita dam, during
EQ5 that shows an asymmetry in the peak values. . . . . . . . 207
xxvi
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxvii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxviii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
6.35 Contours of pore water pressure in the dam after the end of the
reservoir impoundment (positive values correspond to com-
pressive stresses). . . . . . . . . . . . . . . . . . . . . . . . . . 235
6.36 Flow net (contours of stream and potential functions) in the
clay core after the end of the reservoir impoundment. . . . . . 236
6.37 Calculated and measured crest settlement history of the dam
crest prior to the EQ events. . . . . . . . . . . . . . . . . . . . 237
6.38 End of static analysis: (a) Deformed mesh and (b) Vectors of
accumulated displacement of La Villita dam. . . . . . . . . . . 239
6.39 Location of the acceleration monitoring points of the upstream
corner of the crest (U), mid-crest (C) and downstream berm
(B) and location of application of the seismic excitation at the
rock (R) underneath the dam. . . . . . . . . . . . . . . . . . . 240
6.40 Acceleration time history at the crest of La Villita dam during
EQ2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 241
6.41 Response spectra (damping, ξ = 5%) at the crest of La Villita
dam during EQ2 . . . . . . . . . . . . . . . . . . . . . . . . . 241
6.42 Acceleration time history at the crest of La Villita dam during
EQ5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 242
6.43 Response spectra (damping, ξ = 5%) at the crest of La Villita
dam during EQ5 . . . . . . . . . . . . . . . . . . . . . . . . . 242
6.44 New calibration of the Logarithmic CNL model against the
empirical curves of Vucetic and Dobry (1991) for the clay core
- Stiffness degradation (G/G∗max − γ) . . . . . . . . . . . . . . 245
6.45 New calibration of the Logarithmic CNL model against the
empirical curves of Vucetic and Dobry (1991) for the clay core
- Damping (ξ − γ) . . . . . . . . . . . . . . . . . . . . . . . . 245
6.46 New calibration of the Logarithmic CNL model against the
empirical curves of Seed et al. (1986) for the sand filters -
Stiffness degradation (G/G∗max − γ) . . . . . . . . . . . . . . . 246
6.47 New calibration of the Logarithmic CNL model against the
empirical curves of Seed et al. (1986) for the sand filters -
Damping (ξ − γ) . . . . . . . . . . . . . . . . . . . . . . . . . 246
xxix
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxx
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxxi
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxxii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxxiii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxxiv
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxxv
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
List of Tables
xxxvi
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
xxxvii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Nomenclature
a acceleration
ao amplitude of harmonic acceleration load
am maximum value of random acceleration load
a(t) acceleration load time history
B breadth of dam slope face
c parameter of the cyclic nonlinear model
c0 cohesion
[C] damping matrix
Cp hydrodynamic pressure coefficient
d maximum element size
E Young’s modulus
Ed∗ three dimensional strain invariant
f frequency of vibration
F yield function
f1 fundamental frequency of the reservoir
Fdyn total hydrodynamic force
FH horizontal pseudo-static force
Fst total hydrostatic force
FV vertical pseudo-static force
F (ω) transfer function
FA Fourier Amplitude
g acceleration of gravity
g(θ) gradient of the yield function in the J − p0 plane
xxxviii
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
G shear modulus
Gmin , Gmax minimum and maximum values of the shear modulus
H height of the dam
I moment of inertia
Jo (·) Bessel function
J∗ three-dimensional stress invariant
JL parameter of the cyclic nonlinear model
[K] stiffness matrix
K permeability of materials
Kd bulk modulus of the dam materials
kH horizontal pseudo-static coefficient
Kw bulk modulus of water
kV vertical pseudo-static coefficient
L length
m mass
[M ] mass matrix
p0 mean effective stress
pf pore fluid pressure
Pst static water pressure
Pdyn hydrodynamic pressure
r distance of the boundary from the excitation source
Sa spectral acceleration
t time
T natural period of vibration
u displacement
u̇ velocity
ü acceleration
Un (·) shape of the nth mode of vibration
Vp p-wave velocity
Vs shear wave velocity
W width of the dam crest
x horizontal direction
y vertical direction
xxxix
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
z elevation
α parameter of the CH time-integration scheme
αf parameter of the CH time-integration scheme
αm parameter of the CH time-integration scheme
Γ(·) Gamma function
γ shear strain
γw unit weight of water
δ parameter of the CH time-integration scheme
θ Lode’s angle
λ wavelength
ν Poisson’s ratio
ξ damping
ρ mass density
ρ∞ spectral radius at infinity
σ normal stress
τ shear stress
φ angle of shearing resistance
χL parameter of the cyclic nonlinear model
ψ angle of dilation
ω circular frequency of harmonic load
ωd fundamental circular freqency of the dam
ωr fundamental circular freqency of the reservoir
Ω frequency ratio
xl
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Chapter 1
INTRODUCTION
1.1 Background
In recent years, the need for services to the global community resulted in
the construction of various types of infrastructure around the world, such as
highways, bridges and dams. A considerable number of earth and rockfill
dams are built in areas of significant seismic activity and some of them have
experienced earthquakes with catastrophic results. These dams are in need
of seismic retrofit, which includes seismic analysis, re-design and upgrade
actions. Moreover, new dams proposed to be built in such areas need to
be designed to resist expected earthquake loads according to local seismic
hazard studies and achieve a certain level of seismic performance.
Additionally, dams built in areas of insignificant or low seismic activity
have up to now been designed without any consideration of seismic loads.
Though, the design life of dams is quite large and the safety of the dam
structure is quite critical for loss of life and socio-economic losses. Conse-
quently, the need of earthquake design even for these areas of low seismicity
has been outlined recently (e.g. introduction of Eurocode 8) and it depends
on the expected ground motions based on probabilistic hazard assessments
and the design life of structures.
The performance of dams during construction, impounding and opera-
tion can be evaluated reasonably well as the methods used for static analyses
1
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
2
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
dams using numerical analysis. The La Villita earth dam in Mexico was cho-
sen as a case study and advanced nonlinear elasto-plastic finite element anal-
yses were performed. The calculated response was compared to the recorded
behaviour, exhibiting a good agreement. Finally, using the same dam, sev-
eral parametric studies were carried out to assess the effects of different mod-
elling approaches on the predicted seismic response of earth dams. Therefore,
comments were made regarding the applicability of the adopted modelling
techniques and constitutive assumptions.
3
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
4
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
interaction and (e) stiffness inhomogeneity. Previous and where available re-
sults from the literature were also discussed and compared to the outcomes
of this study. Consequently, comments were made regarding the applicability
of the adopted modelling technique and constitutive assumptions.
Finally, Chapter 8 summarizes the findings of this study and describes
the various recommendations for further investigation. At the end, some
thoughts are discussed related to the relative advantages and use of analytical
and numerical methods of analysis.
5
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Chapter 2
SEISMIC PERFORMANCE
OF DAMS
2.1 Introduction
This chapter describes issues related to the seismic performance of dams.
Information drawn from the literature is presented in five sections.
Firstly, a description of the various types of dams is presented in order
to obtain an insight into their materials and geometry and to identify any
vulnerable regions within the dam structure. This is complemented by a
description of a typical construction sequence of earth dams in order to un-
derstand the stress history within the dams and the loading conditions prior
to the seismic events. Secondly, the seismic behaviour of soil deposits is ex-
plored focusing on wave propagation, elastic ground response and nonlinear
dynamic soil behaviour. Subsequently, the behaviour of dams is investigated
and their dynamic characteristics based on both material and geometric con-
straints are discussed. This is followed by a discussion about the hydrody-
namic pressures on dams imposed by the upstream reservoir under various
conditions. Finally, building on the seismic behaviour of dams, the effects
of earthquakes on dam structures are evaluated and the consequences af-
ter a seismic event are detailed through information from post-earthquake
observations and field investigation reports.
6
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
7
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
a simple embankment of well compacted earth and they are generally very
wide (slopes 1:3) in the cross section and not very high (Fell et al., 2005).
Their major use today is for coffer dams, i.e. as temporary water retaining
structures, and they may consist of a low permeability material such as clay.
LOW PERMEABILITY
MATERIAL SUCH AS
CLAY
The second general type of earth dams, which is the most widely used,
is the zoned earth dam. These structures consist of a number of zones de-
pending on the soil used (Golzé, 1977). Generally, there are three main zones
(according to their function), starting from a central core which consists of a
low permeability material (such as clay) and its function is to minimize seep-
age of water through the dam section. Immediately after the central core
there is a thin zone of filters (which is mostly found at the downstream side
of the core). Filters consist of a coarser material, usually coarse aggregates,
similar to those used in concrete mixing, and act as a drainage path for the
water which has seeped through the core. This material is not very fine, in
order to allow some seepage, and also not very coarse so that erosion of the
core is prevented by stopping fine core material being washed through. The
final and outer zone is mainly constructed using a coarse material, which is
usually rockfill. Its function is to provide support and stability of the inner
seepage control system.
3 1. CLAY CORE
4 2. SAND FILTERS
1 3. ROCKFILL
2 4. RIP-RAP
8
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
from erosion due to waves generated in the reservoir. Figure 2.2 shows a
sketch describing the geometry and basic features of a zoned earth dam,
whereas Figure 2.3 shows the Kourris dam which is the largest zoned earth
dam in Cyprus. Zoned earth dams are also called “clay core dams” because
of the core’s major function.
Figure 2.3: The Kourris dam is the largest zoned earthfill dam in Cyprus
(Cyprus Public Works Department, 2010).
2
1 1. CONCRETE FACE
2. ROCKFILL
9
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 2.5: The Kannaviou dam, built in 2005, is the first and until today
(2012) only Concrete-Faced Rockfill Dam in Cyprus (Cyprus Public Works
Department, 2010)
10
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
SOLID
CONCRETE
Buttress dams are made from concrete or masonry and they are very sim-
11
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 2.7: The Lympia dam, built in the late 1970’s, is an example of a
small gravity dam in Cyprus (Cyprus Public Works Department, 2010)
Figure 2.8: The Stwlan dam in Wales is a buttress dam (British Dam Society,
2010)
12
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
ilar to the gravity dams. They have a watertight upstream side supported
by triangular shaped walls, called buttresses. The buttresses are spaced at
intervals on the downstream side. The main difference between gravity and
buttress dams is that they use less material because of the clear spaces be-
tween the buttresses. The buttresses act as stiffeners providing flexural sup-
port to the dam body and therefore the stability of the dam is not provided
entirely from its weight, but also from its strength. Their upstream face, is
not vertical but inclines about 25o to 45o from the vertical, so the resulting
force of the water on the upstream face inclines towards the foundation. The
Stwlan Dam in Wales shown in Figure 2.8 (British Dam Society, 2010) is
such an example. The figure shows the buttresses on the downstream side of
the dam which provide additional resistance to the dam body.
Concrete arch dams are found in very steep and narrow canyons where
space is not adequate to build a large embankment or gravity dam. These
dams are normally thick at the base and they become significantly thinner
when they reach the top of the crest. They are curved in the plan view, with
the convex side of the arch being on the upstream face (U.S. Army Corps of
Engineers, 1994). The reason for constructing dams with such geometry is to
take advantage of the concrete’s high compressive strength, as an arch loaded
on the convex side is able to carry great compressive stresses. Therefore, the
narrow canyon geometry, in which a small-span arch is built, provides great
resistance to the water load. An example of an arch dam is the Baitings dam
in Yorkshire shown in Figure 2.9 (British Dam Society, 2010). It may be
observed from the figure that the reservoir water is on the convex side of the
arch.
Roller Compacted Concrete (RCC) dams are a special type of gravity
dams built using the RCC construction method. Roller-Compacted Concrete
is a special blend of concrete that has the same ingredients as normal concrete
but in different mixing ratios (U.S. Army Corps of Engineers, 2000). Cement,
water and aggregates are the main constituents, as in the general concrete
mixing, but RCC is much drier and it is placed in a manner similar to paving.
The dam is built in successive horizontal layers resulting in a downstream
slope that makes a concrete staircase. Once a layer is placed, it can very
13
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 2.10: The Willow Creek Dam in Oregon, USA is the first Roller
Compacted Concrete dam in the world (US Army Corps of Engineers, 2010)
14
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
quickly support large loads such as the equipment to place the next layer.
After RCC is deposited on the lift surface, small dozers typically spread it
in moderately thin layers of about 0.30m.
RCC materials offer a number of advantages to the construction of con-
crete dams. Firstly, it is more economic than normal gravity dams and an
RCC dam generally requires less material than a gravity dam of the same
height. Moreover, the quick curing and strengthening of the concrete result
in rapid construction as the successive layers of concrete can be constructed
in a very short time. On the other hand however, there is a number of disad-
vantages associated with RCC dams. The main disadvantage is the seepage
control, especially between two successive horizontal layers of concrete (con-
struction joints). For this reason, particular attention needs to be paid to
watertightness and the design of the joints within the dam body, especially
in the case of extreme loading such as earthquakes. The Willow Creek Dam
in the Morrow County of Oregon, US was the first RCC dam to be built
in the world (1979-1983) and it is shown in Figure 2.10 (US Army Corps of
Engineers, 2010).
The previous sections described the various types of earth and concrete
dams and their general features. However, as this work is related to the
performance of earth dams, the rest of this review and research is focused on
the case of earth dams.
15
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
16
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
17
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
the optimum value. For this reason, the core is wetted as it is constructed,
so that it does not lose its moisture.
d. Spillway
Once the embankment is finished, the spillway is constructed in order to
accommodate overflow of the water safely. The spillway is usually placed
next to the embankment on stiff rock abutments and its crest is lower than
the embankment’s crest to allow for some freeboard. The rock abutments
on the sides of the river canyon are often strengthened (stabilized) with
anchors as the excavation of rock (especially if using blast techniques) during
construction of the spillway may affect their stability.
e. Additional facilities
Finally, after construction of the main embankment and its operational
facilities, the dam structure is complete and additional features may be added
to it for secondary purposes. In cases where a dam serves as a link between
the two sides of the river (as this is the case for some long dams), a road
might have been included in the design to provide access. Construction of
this piece of infrastructure is the final part of the construction process even
though it is not directly related to the prime function of a dam.
Besides, if the dam is used for power generation, the hydroelectric facili-
ties are also constructed nearby. They need to be close to the embankment
and in a convenient place so that there is way for incoming water from up-
stream and discharge in the downstream. Moreover, there needs to be access
to the power plant and perhaps additional space for monitoring facilities.
18
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 2.11: Normal faulting: the hanging wall has moved downwards with
respect to the footwall block (US Geological Survey, 2010)
are the result of movement of magma in the interior of the earth due to
temperature difference with the depth. This disturbance in the earth medium
propagates in a wave form in two mechanisms, compression and shear. These
mechanisms create a system of both longitudinal and transverse body waves
(also called pressure, P- and shear, S-waves respectively), as shown in Figure
2.12 (Goree, 2010). These body waves interact with each other and with
other boundaries and create more types of waves, generally called surface
waves, i.e. Rayleigh (1885) & Love (1927) waves (Kramer, 1996).
The motion of these waves causes the soil structure to deform and the soil
properties to change. Soil density tends to increase due to the relocation of
soil particles, a phenomenon called dynamic compaction (or densification). If
water is present, the strength can reduce significantly because of excess pore
19
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
20
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
conducted and then the site response problem concentrates on the response
of a soil deposit to the motion of the underlying bedrock.
The engineering significance of ground response is that the soil layer may
amplify or de-amplify the bedrock motion and even change its frequency
content. Moreover, surface topography (such as a triangular dam) and basin
effects (such as river canyons) are also important, because they affect the
dynamic response of the dam, as geometry affects the medium’s stiffness.
Several methodologies for ground response analysis were developed over time,
both analytical (Thompson, 1950; Ambraseys, 1959a; Seed and Idriss, 1969;
Gazetas, 1982b; Vrettos, 2013) and numerical (Haskell, 1953; Biot, 1963;
Kausel and Röesset, 1981; Schnabel et al., 1970; Phillips and Hashash, 2009).
Analytical methodologies adopt various simplifying assumptions and provide
valuable information about the dynamic response (such as natural periods
of vibration, amplification of the motion etc.) of soil layers (damped or
undamped) on top of the bedrock (rigid or not).
u
z uniform
visco-elastic vertically
H soil layer propagating
shear wave
rigid bedrock
umax (z, t)
F (ω) = (2.1)
umax (H, t)
where ω is the circular frequency of the applied harmonic load, umax is
21
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1
|F (ω)| = q (2.2)
cos2 (ωH/Vs ) + [ξ(ωH/Vs )]2
where Vs is the shear wave velocity of the soil, H is the height of the soil
deposit, ω is the circular frequency and ξ is damping ratio.
The frequency spectrum of the amplification (Equation 2.2) is shown
graphically in Figure 2.14 for three values of the damping ratio, ξ = 5%, 10%
and 20%. It is noticed that peak values of amplification of the response occur
at different values of the circular frequency, ωn , as described by Equation 2.3,
which correspond to the natural modes of vibration of the soil layer. It is
also shown that the higher modes of vibration are significantly affected by
the value of the damping ratio.
10
0
0 1 2 3 4 5 6
Frequency Ratio, ω/ω1 []
2n − 1 πVs
ωn = (2.3)
1 2H
22
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1 2 0.637
|F (ω1 )| = |F |max = q = = (2.4)
cos2 (π/2) + [ξ(π/2)]2 πξ ξ
πVs
ω1 = (2.5)
2H
Therefore, the first three natural periods of vibration (i.e. associated with
the first, second and third natural modes) are given by Equation 2.6.
T1
1/ω1
1
4
4H H
T2 = 2π 1/ω2 = 1/3 = 1.333 (2.6)
Vs
Vs
T3 1/ω3 1/5 0.8
23
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
The slope of the curve at each point (i.e. for each value of shear strain, γ)
is the tangential shear modulus, Gtan , whereas the ratio of the current (max-
imum values in a cycle) shear stress to strain is the secant shear modulus,
Gsec (denoted as G in Figure 2.15). The area of a hysteresis loop represents
the amount of energy dissipated during that cycle. The damping ratio, ξ
(denoted as D in Figure 2.15) is defined as the ratio of the energy dissipated
over the maximum strain energy, and it is given by Equation 2.8.
1 ∆W
ξ= (2.8)
4π W
where ∆W is the area of the hysteresis loop, denoting the energy dissi-
pated in one cycle and W is the area of the triangle with sides τa (= τmax )
and γa (= γmax ) (Figure 2.15) corresponding to the peak energy during that
cycle (Kramer, 1996).
As shown in Figure 2.15, the values of shear modulus G and damping
ξ change with shear strain γ. According to their definition, G decreases
whereas ξ increases with increasing γ. These depend also on the confining
pressure, σ00 . Figure 2.16 describes the relation between G, ξ and γ, arising
from experiments on Toyoura sand (Kokusho, 1980).
24
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 2.16: Relation between shear stiffness, G, damping, ξ and cyclic shear
strain γ, arising from experiments on Toyoura sand (Kokusho, 1980)
25
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
26
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.7 0.7
0.6 0.6
0.5 0.5
0.4 0.4
0.3 0.3
0 0
−0.5 0 0.5 1 −1 −0.5 0 0.5 1
Modal Displacement, Un(y) Modal Displacement, Ur(z)
Figure 2.19 shows the first two natural modes of vibration of a dam in
the upstream-downstream direction for a complete wedge (k 0 = 0), i.e. for
(a) n = 1 and r = 1 and (b) n = 1 and r = 2 (refer also to Figure 2.18).
However, if an infinitely wide canyon (L/H → +∞) is considered, the un-
damped circular frequency of vibration of the nth mode is given by Equation
2.11.
27
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 2.19: First two natural modes of vibration of a dam in the upstream-
downstream direction for a complete wedge (k 0 = 0) (Ambraseys, 1960a).
Vs
ω n = an (2.11)
H
Therefore, the first three natural periods of vibration (i.e. associated with
the first, second and third natural modes) are given by Equation 2.12.
T1
2.613
H
T2 = 1.138 (2.12)
Vs
T
0.7261
3
y
Un (y) = Jo an (2.13)
H
It is worth noting that according to the shear beam approach of Am-
braseys (1960a), the frequencies and modes of natural vibration do not de-
pend on the slope angle of the dam (Equations 2.12 and 2.13).
28
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1
|F (ω)| = (2.14)
Jo ( ωH
Vs
)
12
10
0
0 1 2 3 4 5 6 7
Frequency ratio, ω/ω []
1
29
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
30
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
4
f3D /f 2D
1 2 3 4 5 6 7
L/H
Figure 2.21: Comparison between natural frequencies from 2D and 3D anal-
yses of dams for different length to height, L/H ratios (Mejia and Seed,
1983).
T1 H
0.5
T1
8
0
0 1 2 3 4 5 6
L/H
31
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
DAMPING 10%
MIDCREST AMPLIFICATION FUNCTION
SEMI-CYLINDRICAL
CANYON
PLANE STRAIN
SHEAR BEAM
32
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0 1/2
G = 1000K2 (σm ) (2.15)
33
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
the 2/3 power of the depth from the crest in earthfill dams.
Figure 2.24: The inhomogeneous dam system studied by Dakoulas and Gaze-
tas (1985).
where z and H are the vertical distance from the crest of the dam and
the height of the dam respectively as shown in Figure 2.24, Gb is the shear
modulus at the base of the dam and m is exponent of the stiffness which
dictates its variation with depth.
Figure 2.25 shows the displacement modal shapes for different values of
the inhomogeneity exponent parameter m (Equation 2.16). The stiffness is
varied from constant (m=0) to linear (m=1) with depth, i.e. from zero to
large inhomogeneity. It is shown that inhomogeneity has a significant effect
on the displacement modal shapes and results in a concentration of larger
displacements at the crest of the dam.
Figure 2.26 shows the ratio of the nth natural period, Tn (m) (correspond-
ing to the nth vibration mode) of an inhomogeneous dam with inhomogeneity
exponent parameter m, over the period of a homogeneous dam (m=0). The
stiffness is varied from constant (m=0) to linear (m=1) with depth, and re-
34
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
35
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
sults are shown for two values of the truncation ratio, λ = h/H (Figure 2.24).
It is shown that inhomogeneity affects the periods corresponding to higher
modes of vibration (n = 2, 3, 4, ...) only for an untruncated dam (λ = 0, i.e.
for a complete triangle). However, in any case it does not affect significantly
the fundamental period of the dam T1 which seems to be constant for any
value of inhomogeneity parameter m and truncation ratio λ.
The effect of inhomogeneity on the variation of accelerations in a dam
was presented by Gazetas (1987) who showed that the accelerations in in-
homogeneous dams are higher close to the crest of the dam. That study
considered a 120m high uniform inhomogeneous dam under the Taft earth-
quake record scaled at 0.4g. Figure 2.27 (a) shows the variation of the peak
value of acceleration with the depth of the dam. It is clear that for highly
inhomogeneous (G ∼ z2/3 ) dams, higher accelerations are concentrated at
the top of the dam.
The previous conclusions were drawn from studies that considered linear
elastic behaviour of the dam materials. However, as reported by Gazetas
(1987), when nonlinear material behaviour is introduced, the effects of ma-
terial inhomogeneity are less pronounced as shown in Figure 2.27 (b) (Stara-
Gazetas, 1986). The figure shows the variation of the peak value of acceler-
ation in a highly nonlinear inelastic dam. That study (Stara-Gazetas, 1986)
considered a 120m uniform dam subjected to the 1952 Taft record. From the
latter figure, it is evident that the effects of material inhomogeneity diminish
and this could be due to the dissipation of energy through inelasticity.
Considering the previous studies from the literature, it may be concluded
that inhomogeneity affects the dynamic response of earth dams. This in-
fluence is mainly related to the modal shapes and the vertical variation of
accelerations. However, inhomogeneity has a smaller impact on the higher
natural periods of vibration and minimal impact on the fundamental period
of the dam. Finally, it was reported in the literature that the effects of inho-
mogeneity are less pronounced if nonlinear material behaviour is considered.
36
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
37
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
found that the amplitude of the dam response is affected by the impedance
between the dam and the canyon.
Shear beam analyses by Dakoulas and Gazetas (1987) investigated the
effect of the presence of a foundation layer on the fundamental period of a
dam. The study considered a uniform distribution of shear modulus G for
both the dam and the foundation layer. Figure 2.28 shows the ratio of the
period of a dam-layer system T f , over the period of a dam founded on rigid
1
bedrock T1 , for different dam length over height, L/H, and foundation-dam
height, h/H ratios.
38
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
(a) (b)
addition, the effect of the presence of the layer is more pronounced for longer
dams and for softer foundation layers (as compared to the stiffness of the
dam).
The flexibility of the canyon rock was also examined by Papalou and
Bielak (2001, 2004) when studying the behaviour of La Villita dam in Mexico.
The dam was modelled as a shear beam, whereas the surrounding canyon was
discretised using finite elements. It was shown that the acceleration response
of the dam reduces as the stiffness of the rock canyon reduces because of
waves radiating away. Moreover, accelerations reduce further if nonlinear
behaviour is considered. Figure 2.31 shows the acceleration time histories
predicted at the crest of the dam for different shear wave velocities of the
rock, Vr , for (a) a linear elastic and (b) a nonlinear system.
It is therefore noticed that the presence of an underlying foundation allu-
vial layer tends to soften the dam-foundation system and hence increase its
fundamental period. In addition, the presence of a dam foundation has an
effect on the dam accelerations.
39
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
(a)
(b)
(c)
40
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
(a) (b)
Figure 2.31: Accelerations predicted at the crest of La Villita dam for dif-
ferent shear wave velocities, Vr of the rock canyon for (a) a linear elastic
(Papalou and Bielak, 2001) and (b) a nonlinear system (Papalou and Bielak,
2004).
41
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
NO RESERVOIR
0 2 4 6 ω
ωd
20 0.4
1
10
1.8
ωr
ωd
(a)
2
0 2 4 6 ω
20 ωd
10
4
5
ωr
ωd
(b)
Figure 2.32: Frequency response of the accelerations at the dam crest for
different load, dam and reservoir circular frequency ratios as presented by
Chopra (1968): (a) for 0 ≤ ωr /ωd ≤ 1.8 and (b) for 2 ≤ ωr /ωd ≤ 5.
42
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 2.32 shows that the presence of the reservoir affects significantly
the dynamic response of the dam. More specifically, the amplification of
the earthquake motion in the dam depends not only on the frequency of the
loading but also on the ratio of the natural circular frequency of the reservoir
over the natural circular frequency of the dam.
The effects of dam-reservoir interaction on the dynamic response of dams
were evaluated by Hall and Chopra (1980) with the finite element method
assuming elastic dam response and considering the compressibility of water
(Hall and Chopra, 1982c) (Figure 2.33). It was shown that the presence
of the reservoir increases the fundamental period, T of the dam-reservoir
system (i.e. it becomes more flexible). This is because the mass of the water
is added to the mass, m, of the dam without though adding any flexural or
shear stiffness.
43
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
|F|
ω/ω
r
ωd
|F|
ω/ωr
ωd
Figure 2.34: Frequency spectrum of the accelerations at the dam crest for
(a) a concrete gravity dam (Hall and Chopra, 1982b) and (b) an earth dam
(Hall and Chopra, 1982a)
44
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
This is more pronounced in concrete dams (Hall and Chopra, 1982b) (and
more specifically thin concrete arch dams) which have a smaller volume (and
hence mass) than earthfill dams and usually a nearly vertical upstream face.
In contrast, the effects of reservoir-dam interaction on the dynamic response
of earth dams (mainly the fundamental period of the dam) are believed to
be less important (Hall and Chopra, 1982a). This is due to the fact that
(a) the volume (and hence the mass) of earth dams is very large compared
to the added mass from the water (see also Section 2.5) and (b) the sloped
upstream face of earthfill dams experiences smaller hydrodynamic pressures
from the reservoir than a vertical face.
Figure 2.34 (a) shows the frequency spectrum, (|F |), of the amplification
of accelerations at the dam crest for concrete dams (including damping, ξ =
5%), whereas Figure 2.34 (b) shows the corresponding spectrum for earth
dams (including damping, ξ = 10%). The horizontal axis shows the ratio
of the frequency of a harmonic load, ω over the fundamental frequency of
the reservoir, ωr . On the same axis, the values of the natural frequency
of vibration of the dam, ωd are also included. Both figures show that the
peaks in the acceleration response (which correspond to the natural modes of
vibration) occur at smaller frequencies for the case considering the reservoir
water (solid line) than the case ignoring the reservoir (dotted line). This
means that the natural frequencies of the system decrease and therefore the
system softens. Furthermore, it may be observed that the softening effect of
reservoir-dam interaction is larger in concrete dams than earth dams. More
specifically, the reduction in the fundamental frequency of the dam-reservoir
system is about 30% for concrete dams, whereas for earth dams it seems to
be less than 5%.
However, as far as the amplification is concerned, there is no clear trend
shown in these results. In Figure 2.34 (a) for concrete dams, the amplification
for the case considering the reservoir water is larger than the corresponding
amplification ignoring the water for the first natural mode, whereas it is
smaller for the higher modes of vibration. Similar observations hold for the
case of the earth dam (Figure 2.34 (b)).
It should be noted that the calculated first four natural circular frequen-
45
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
cies of the concrete dam (Hall and Chopra, 1982b) (ωd = ω1 , ω2 , ω3 and ω4 )
and the corresponding first and second calculated circular frequencies for the
earth dam (Hall and Chopra, 1982a) (ωd = ω1 and ω2 ) are also shown on
Figures 2.34 (a) & (b). It may be noticed that in the cases studied by Hall
and Chopra (1982b,a), the natural circular frequencies of the dam were not
at the proximity of the natural circular frequencies of the reservoir, ωr . It
would be therefore interesting to investigate what would be the response of
the dam in such cases where resonance is expected to occur.
Examples of dams for which the measured natural period with full reser-
voir was found to be larger than for partly full reservoir include the Pine Flat
dam in California, U. S. (Chopra and Gupta, 1978) and the Morrow Point
dam in Colorado, U. S. (Duron and Hall, 1988; Fenves et al., 1992; Tan and
Chopra, 1996). The Pine Flat dam is a 122m high concrete gravity dam
whose fundamental natural period was measured by forced vibration tests
to be 0.288 sec and 0.306 sec with the reservoir depth at 95m and 105m re-
spectively. The Morrow Point dam is a 142m high concrete arch dam whose
fundamental period was determined similarly by forced vibrations tests to
be 0.268 sec with the reservoir partially full and 0.303 sec with a full reser-
voir (Chopra, 1995). More information regarding reservoir-dam interaction
effects in embankment dams may be found in the work of Zhao et al. (1993)
and Guan and Moore (1997).
46
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
47
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
po
Figure 2.35: The dam-reservoir system considered by Westergaard (1933)
∞
ac T 2 2πt X 1 −qn nπy
v= 3
cos( ) e cos( ) (2.18)
π T n=1,3,5,... ncn 2h
where,
s
16γw h2
cn = 1− (2.19)
n2 gKT 2
nπcn x
qn = (2.20)
2h
and ac and T are the coefficient and period of the horizontal vibrations,
γw is the unit weight of the water, x is the horizontal distance from the
upstream face of the dam, y is the vertical distance from the top of the
reservoir, h is the height of the dam (Figure 2.35), K is the compressibility
of the water and t is the time. Please note that according to Westergaard
(1933), the acceleration coefficient ac is a dimensionless parameter, equal to
the peak value of acceleration, ao , divided by the acceleration of gravity, g
(=9.81 m/s2 ).
The hydrodynamic pressure distribution on the dam face is given by
Equation 2.21 and its maximum value for y = h is given by Equation 2.22.
∞
8ac γw h X 1 nπy
p(y) = 2
sin( ) (2.21)
π n=1,3,5,... ncn 2h
48
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
∞ n−1
8ac γw h X (−1) 2
pmax = (2.22)
π 2 n=1,3,5,... ncn
For incompressible fluid, i.e. the bulk modulus of water Kw → +∞,
therefore cn → 1− and hence the maximum pressure is given by:
8ac γw h 1 1
pmax = 2
(1 − 2 + 2 − ...) = 0.743ac γw h (2.23)
π 3 5
The corresponding total hydrodynamic force over the upstream face of the
dam (by integrating the hydrodynamic pressure over the vertical distance)
for incompressible reservoir water is given by Equation 2.24.
16ac γw h2 1 1
Fmax = 2
(1 + 2 + 2 − ...) = 0.543ac γw h2 (2.24)
π 3 5
This hydrodynamic pressure is considered to be constant and added to
the pseudostatic calculations for the stability of the dam. Recognising the
inertial contribution of this pressure, Westergaard (1933) suggested that the
effect of the hydrodynamic pressures could also be simulated if an additional
mass of water was added to the mass of the dam, instead of adding an external
hydrodynamic stress. This is the so-called added mass concept. The width,
b of this water mass added to the dam follows the distribution given in
Equation 2.25. Additionally, if the water mass was replaced by an equivalent
concrete mass, the width, b0 of this mass would be given by Equation 2.26.
Both distributions are shown in Figure 2.36 (Westergaard, 1933).
7q
b= hy (2.25)
8
q
b0 = 0.38 hy (2.26)
49
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
b'
h reservoir dam
0.38 h
7/8 h
50
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
∂ 2φ
Vp2 ∇2 φ = (2.29)
∂t2
The water pressure at any point p(x1 , x2 , t) (x1 is the vertical distance
from the bottom of the reservoir and x2 is the horizontal distance from the
upstream face of the dam) is given by Equation 2.30.
∞
(−1)n+1 −µn x2 µ 2 x1
p(x1 , x2 , t) = −2ac ω 2 ρh sin ωt
X
exp cos (2.30)
n=1 µ2n h h
where µn = (2n−1)π
2
, ρ is the density of water, ω is the circular frequency
of the excitation and h is the height of the reservoir.
The maximum value of the water pressure on the dam face (x2 = 0)
occurs at the bottom of the reservoir (x1 = 0) and it is given by Equation
2.31.
ẍo
pmax = 0.743ẍc ρh = 0.743 γw h (2.31)
g
where, ẍo is the maximum value of the horizontal acceleration.
Considering the three above-mentioned approaches (Westergaard, 1933;
von Kármán, 1933; Newmark and Rosenblueth, 1971) to estimate the hy-
drodynamic pressures on vertical rigid dams due to incompressible reservoir
under harmonic loading, it may be concluded that the maximum hydrody-
namic pressure occurs at the bottom of the reservoir and is given by Equation
2.32.
51
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
tions of fluid mechanics. Zangar (1952) noticed the relation of the resulting
equations with the equations describing electricity and then performed ex-
periments using an electric analog. He presented results obtained from both
the experiment and an analytical expression. The hydrodynamic pressures,
p on a dam face are given by Equation 2.33.
θo
!
C(y = 0) = 0.743 1 − o (2.35)
90
where, y is the vertical distance from the free surface, h is the height of
the reservoir and Cm is the maximum value of C. The variation of C and
Cm with the value of the slope angle, θ, are shown in Figure 2.38.
Figure 2.37 shows that the distribution of the water pressure along the
vertical dam direction, normalised with respect to its magnitude at the base
of the dam has a similar shape to the one with a vertical upstream face (θ = 0
in Figure 2.37). However, the maximum value of the hydrodynamic pressure
is not always at the bottom of the reservoir, but depending on the inclination
of the upstream face, it is located within the bottom fifth of the dam.
Following the analytical and experimental work of Zangar (1952), an an-
alytical approach was followed by Chwang and Housner (1978) and Chwang
(1978) in order to calculate the hydrodynamic pressures on dams with in-
clined upstream face. Chwang and Housner (1978) adopted the momentum-
52
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
53
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
s
h y
x θ=απ
Figure 2.40: Values of the pressure coefficient C(y) for a dam with an inclined
upstream face. The dashed line presents the results from the work of Chwang
and Housner (1978) whereas the solid line corresponds to the work of Chwang
(1978).
54
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
balance method (von Kármán, 1933) whereas Chwang (1978) used two-
dimensional potential-flow theory, which he considered to give an exact solu-
tion. In both approaches, the dam was considered to be rigid and subjected
to horizontal forces. The geometry of the system used is shown in Figure
2.39.
For both approaches, the hydrodynamic pressures may be simply obtained
using Equation 2.36.
p(y) = Cp ao ρh = Cp ac γw h (2.36)
where, γw and ρ are the unit weight and mass density of water respectively,
ao is the maximum value of the horizontal harmonic acceleration, ac is the
acceleration coefficient (=ao /g), h is the height of the reservoir and Cp is the
pressure coefficient, which can be determined from Figure 2.40. The results
shown in dashed lines come from the work based on the momentum-balance
method (Chwang and Housner, 1978), whereas the results in solid lines come
from the work based on the potential-flow theory (Chwang, 1978).
From the previous studies, it may be concluded that the hydrodynamic
pressures on a dam with an inclined upstream face are generally smaller than
those on a dam with a vertical upstream face. Moreover, the maximum value
of the hydrodynamic pressure does not occur at the base of the dam (which
is the case for a dam with a vertical upstream face) but slightly higher.
hρao 1
p(y) =
sin aπ sin βπ Γ(a)Γ(β)Γ(1 − a − β)
55
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
s
h y
x βπ θ=απ
a
y
4[sin(a + γ)π cos βπ + sin(β−)π cos aπ] Z π/2
h (tan 2θ)2β θ
dθ
sin(a + β)π 0 (tan2 θ + hy )a+β sin θ cos θ
" #
sin(a + γ)π cos(a + β)π + sin(β − γ)π h
−ρao − s(y) (2.37)
sin(a + β)π sin aπ
where, aπ, βπ and γπ are the angles of the upstream dam slope, reservoir
bottom and applied acceleration respectively, h is the height of the reservoir,
ρ is the mass density of the water, s(y) is the distance of any point on the
dam face from the bottom corner of the dam and Γ(·) is the Gamma function.
As with the work of Chwang and Housner (1978) and Chwang (1978),
the hydrodynamic pressures can be more practically obtained using Equation
2.38.
p(y) = Cp ao ρh = Cp ac γw h (2.38)
where, γw and ρ are the unit weight and mass density of water respec-
tively, ao is the maximum value of the horizontal harmonic acceleration, ac
is the acceleration coefficient (=ao /g), h is the height of the reservoir and
Cp is the pressure coefficient, which can be determined from Figures 2.42 &
2.43. In Figure 2.42, the results shown in dashed lines come from the exper-
imental work of Zangar (1952) whereas the results in solid lines come from
the theoretical work based on the potential-flow theory of Liu (1986).
Considering the study of Liu (1986), it may be concluded that the hydro-
dynamic pressures on a dam with an inclined reservoir base are even smaller
than the corresponding pressures on a dam with a horizontal reservoir base.
The maximum value of the hydrodynamic pressure does not occur at the base
of the dam, but slightly higher. The value of the hydrodynamic pressure at
56
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
57
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
the base of the dam can be extremely small, depending on the slope angle of
the reservoir base and the upstream dam face.
Fdyn (Ω) 32 +∞
X 1
= 3 q (2.40)
Fst π g n=1 (2n − 1)2 (2n − 1)2 − Ω2
58
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
3.5
2.5
g Fdyn/Fst
1.5
0.5
0
0 1 2 3 4 5 6
Frequency Ratio, Ω = ω/ω1
high for loading frequencies in the vicinity of the natural frequencies of the
reservoir, and specifically close to the first natural frequency. It should be
noted that if the water is incompressible, then the fundamental frequency of
the reservoir, ω1 becomes infinite and hence the ratio Ω = ω/ω1 → 0+ for any
loading frequency ω. The value of the hydrodynamic force corresponding to
this ratio based on Chopra (1967a) is given by Equation 2.41.
Fdyn (Ω = 0) 32 +∞
X 1
= 3 = 0.1106 (2.41)
Fst π g n=1 (2n − 1)3
The corresponding value for incompressible water (i.e. considering Equa-
tion 2.23) suggested by Westergaard (1933) is given by Equation 2.42.
0.543
Fdyn (Ω = 0) g
γw H 2
= 1 = 0.1107 (2.42)
Fst γ H2
2 w
Comparing Equations 2.41 and 2.42, it is concluded that indeed the so-
59
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
lutions of Chopra (1967a) and Westergaard (1933) agree when the water is
considered as incompressible. The work of Chopra (1967a) proved that the
general solution proposed by Westergaard (1933) (considering water com-
pressibility, i.e. Equation 2.22) is valid only for the cases where the excita-
tion frequency is smaller than the fundamental frequency, f1 (Equation 2.39)
of the reservoir. Finally, it proved that the added mass concept of Wester-
gaard (1933) is valid only when the water is considered to be incompressible,
(Kw → +∞, and therefore ω/ω1 → 0+ ).
x 4γw Vp +∞
X (−1)n−1 Z t
p (y, t) = cos (λn y) üg (τ )Jo [λn Vp (t − τ )]dτ (2.43)
πg n=1 2n − 1 0
y 4γw Vp +∞
X 1 Z t
p (y, t) = cos (λn y) v̈g (τ )sin[λn Vp (t − τ )]dτ (2.44)
πg n=1 2n − 1 0
60
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Structure
H Fluid
y
x
Figure 2.45: The fluid-structure system considered by Lee and Tsai (1991)
∂2 ∂ 2u ∂ 2u
!
EI 2 + m 2 = −müg − p(x = 0, y, t) (2.45)
∂y 2 ∂y ∂t
where, E is the Young’s modulus, I is the moment of inertia of the cross-
section about the axis of bending, m is the mass per unit length of the
structure, y is the vertical distance from the base of the dam, u is the hor-
izontal displacement of the dam structure, üg is the ground acceleration in
the upstream-downstream direction and p is the hydrodynamic pressure.
The displacement response of the dam structure, u, including the hydro-
dynamic effects can be expressed as a linear combination of modal shapes,
61
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
+∞
X
u(y, t) = φn (y)Yn (t) (2.46)
n=1
4ρVp +∞
X (−1)k+1 Z t
p(y, t) = cos λk y üg (τ )Jo [λk Vp (t − τ )]dτ
π k=1 2k − 1 0
2ρVp +∞
X +∞
X Z H Z t
+ cos λk y φn (y) cos λk ydy Ÿn (τ )Jo [λk Vp (t − τ )]dτ (2.47)
H k=1 n=1 0 o
where, Vp is the p-wave velocity in the water, H is the height of the dam
and reservoir, üg is the ground acceleration in the upstream-downstream
direction, φn (y) is the nth vibration modal shape of the dam, Jo (·) is a Bessel
function of the first kind and zero order and λk is given by Equation 2.48.
(2k − 1)π
λk = (2.48)
2H
Generally, the flexibility of the dam is assumed to result in smaller val-
ues of hydrodynamic pressures, as part of the energy is absorbed by the
deformation of the dam structure (see also the effect of the reservoir on the
dam accelerations in Section 2.4). Therefore the introduction of flexibility
seems to be beneficial as it reduces the hydrodynamic loads. Centrifuge
studies comparing small-scale models of a stiff concrete dam and a flexible
aluminium dam (Saleh and Madabhushi, 2010) have confirmed this assump-
tion, as the measured hydrodynamic pressures were found to be smaller for
the flexible case.
However, it is the Author’s opinion that the hydrodynamic pressures on
flexible dams have not been investigated (analytically, numerically or ex-
perimentally) exhaustively. Reservoir-dam interaction is considered to be
beneficial, but cases in which this interaction might be detrimental have not
been clearly identified. These could include, for example, the case of reso-
62
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
nance between dam and reservoir. Therefore, further studies are required in
order to confirm these assumptions or identify any limitations.
2.5.7 Comments
This section described several methods for calculating the hydrodynamic
pressures on dams resulting from the upstream reservoir. The hydrodynamic
pressures in the simplest case, of a rigid dam with vertical upstream face and
an incompressible reservoir, reach a maximum at the base of the dam. The
magnitude of the pressures gets smaller for inclined dams and in those cases
the maximum value occurs higher than the base of the dam.
When the compressibility of the water is considered, the magnitude of the
total hydrodynamic force on the dam depends on the ratio of the frequency
of the load to the natural frequency of the reservoir. Furthermore, if the flex-
ibility of the dam structure is considered, the hydrodynamic pressures tend
to be smaller in general because of energy absorbed by the dam. However,
the latter case still has to be further investigated. Information about the
numerical modelling of hydrodynamic pressures follows in Chapter 4.
63
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 2.46: Failure of the upstream slope of the Fategadh dam in India
during the 2001 Bhuj earthquake (Madabhushi and Haigh, 2001).
Crest settlements may also be observed and they usually occur for two
64
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
reasons. Firstly, a large magnitude slope failure (such as flow failure) on both
faces of the dam may alter the dam’s shape in such a way that it becomes
wider at the base and shorter in height. This kind of failure, firstly described
by Ambraseys (1959b) may be observed in Figure 2.47 (Newmark, 1965),
where successive shocks create a series of slope failures and severly change
the geometry of the dam. Secondly, cyclic loading causes relocation of soil
particles, and in a large intensity earthquake, dynamic compaction may occur
in the absense of slope failures. These crest settlements are important and
care should be taken in order to avoid significant loss of freeboard (which
directly affects functionality and safety of the dam) and differential crest
settlements. Both of these events become even more crucial in cases where
a road or other facilites are built on top the dam crest.
65
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
the central crest core because of tensile strains from longitudinal oscillations
or by different lateral response of adjacent material zones. Cracking is im-
portant as it may provide ways for leakage through the dam, erosion of the
internal dam material and finally local piping failures. Figure 2.48 shows an
example of cracking in dams during seismic events: wide and deep cracks on
the upstream face of Fategadh dam in India during the 2001 Bhuj earthquake
(Madabhushi and Haigh, 2001).
Figure 2.48: Cracks on Fategadh earth dam in India during the 2001 Bhuj
earthquake (Madabhushi and Haigh, 2001).
66
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
2.6.2 Liquefaction
Liquefaction is another severe effect of earthquakes on dams. The most com-
mon type of liquefaction observed in earth dams is flow liquefaction which
can be extremely dangerous for the stability of the dam, as it causes large
scale failures. Flow liquefaction occurs when the shear stress for static equi-
librium is larger than the shear strength of the liquefied soil (Kramer, 1996).
It is mainly triggered by the development of excess pore water pressures in
contractive saturated cohesionless zones of the dam (such as rockfill).
Figure 2.49: Failure of the Lower San Fernando dam due to liquefaction in the
upstream part of the rockfill (Cambridge University Engineering Department,
2012).
The most famous case of dam failure due to liquefaction is the Lower
San Fernando dam during the 1971 San Fernando earthquake (Figure 2.49).
Figure 2.50 shows a back-analysis of the liquefaction failure by Seed et al.
(1975). The hatched region in the lower part of the upstream rockfill is the
liquefied material, which caused a global collapse of the whole dam structure.
67
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 2.50: Failure and reconstructed cross section of the Lower San Fer-
nando dam during the 1971 San Fernando earthquake as presented by Seed
et al. (1975).
2.7 Summary
This chapter summarizes the main information found in the literature re-
garding the seismic behaviour of dams. It includes details about the design
and construction of dams, a description of the seismic behaviour of soil and
the dynamic response of dams, the different patterns of the reservoir hydro-
dynamic pressures and concludes with the effects of earthquakes on dams.
The important information found from the literature are summarised below:
68
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
of the soil and an increase of the damping with the induced shear
strains. Reduced soil stiffness results in longer dam periods of vibration
and increased damping results in smaller values of accelerations in the
dam.
69
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
70
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Chapter 3
SEISMIC ANALYSIS OF
DAMS
3.1 Introduction
A number of methods exist and are used today for the evaluation of the
seismic response of dams. Each method carries its own assumptions and
hence has its own limitations. It is therefore essential to be aware of the
limitations of the methods used so that useful analyses are performed.
The evolution of seismic analysis with time is interesting, as it is a valuable
indication of the advancement of science and shows how the main engineering
concerns change with time. The initial consideration of rigid body failure
mechanisms led to the development of simplified methods considering the
sliding mass as an undeformable structure. Field investigations however did
not confirm this assumption and thus more advanced methods considering
the dam as a deformable body were developed (Gazetas, 1987).
The various methods fall into several categories according to the main
objective of the analysis:
(a) assessment of the stability of the dam structure (such as the pseudo-
static method)
71
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
(c) evaluation of the dynamic response of the whole dam (such as the shear
beam and finite element methods)
72
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
kHW
kVW
W
S
N
of the seismic coefficient k for which the slope fails is called critical seismic
coefficient and is denoted as kc . Equations 3.1 and 3.2 define the horizontal,
FH , and vertical, FV , inertial forces respectively.
aH W
FH = = kH W (3.1)
g
aV W
FV = = kV W (3.2)
g
where, aH and aV are the horizontal and vertical accelerations respectively
and kH and kV are the horizontal and vertical seismic coefficients.
The basic idea of this method is that equilibrium calculations are carried
out in a similar fashion to static slope stability problems, but with the addi-
tion of the earthquake inertial loads. For the example of the horizontal load,
the objective of the analysis is either of the following:
73
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
74
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
(1960b, 1962), including the PhD research of the latter (Ambraseys, 1959b)
at Imperial College.
The most widely used method from this group is perhaps the “sliding
block analysis” method, proposed by Newmark (1965) in the fifth Rankine
Lecture. This method, known also as the Newmark’s analysis method, sug-
gests that whenever the acceleration becomes larger than the critical value, ac
(for static slope stability, i.e. FS = 1 and may be obtained using the pseudo-
static method of analysis), permanent displacements occur, according to the
difference of the applied and the critical acceleration. Consequently, by inte-
grating the acceleration time history with respect to time, the velocity time
history is obtained (Equation 3.3), and integrating this again, the displace-
ment time history is obtained (Equation 3.4).
Z t
u̇ = (a − ay )(t) dt ∀ a − ay > 0 (3.3)
0
Z tZ t
u= (a − ay )(t) dt dt ∀ a − ay > 0 (3.4)
0 0
Figure 3.2: Newmark’s sliding block method (Wilson and Keefer, 1985)
Figure 3.2 (Wilson and Keefer, 1985) illustrates this Newmark’s Sliding
Block analysis procedure. It is shown that the integrated difference of accel-
eration with respect to time yields the velocity and further, the integrated
75
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
76
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
77
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
the earthquake are mainly related to seismic wave propagation and hence
vibration mode interaction.
One of the first attempts to consider the dynamic effects of earthquakes is
the Shear beam (SB) (or shear wedge) approach, which is shown in Figure 3.3.
As stated by Gazetas (1987), the origins of this method have been attributed
to Mononobe et al. (1936), who first introduced the 1D shear beam model
for earth dams in the 1930’s.
Hatanaka (1952, 1955) showed that shearing deformations are predomi-
nant compared to the bending deformations and performed 2D shear beam
analyses in rectangular canyons. Ambraseys (1960a,b) extended the shear
beam model to account for a truncated wedge shape, rectangular canyon
and underlying elastic layer. Since then, a number of researchers extended
and used this method in various cases (Gazetas, 1987). Examples include
the work of Abdel-Ghaffar and Koh (1981) and Gazetas (1981a) for 2D lon-
gitudinal vibrations and Gazetas (1981b) for vertical vibrations.
The initial shear beam method (Mononobe et al., 1936; Hatanaka, 1952,
1955; Ambraseys, 1960a) assumes that:
• the dam deforms in simple shear and produces only horizontal displace-
ments (i.e. bending of the dam is ignored)
78
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
z
H τ
dz
cτ dz
τ +cz
2B
Figure 3.3: Shear beam approach for two-dimensional seismic analysis of
dams. After Gazetas (1987).
τ along the shear surface) on both sides of the shear beam (top and bottom)
along with the inertial force.
An example of the relevant theory for an inhomogeneous dam is given
by Dakoulas and Gazetas (1985), who provided solutions for the case of
increasing shear modulus with depth, G(z)=Gb (z/H)m , where Gb =G(z=H)
is the value of the shear modulus at the base of the dam. The resulting
governing differential equation (for a single-material inhomogeneous dam)
refers to one-dimensional wave propagation and it is described by Equation
3.5.
" #
1 ∂ ∂u
ρ(ü + üg ) = G(z)z (3.5)
z ∂z ∂z
where ρ is the density of the material, G(z) the variation of the shear
modulus with depth, z the vertical distance from the crest of the dam and u
is the horizontal displacement, as defined in Figure 3.3.
The boundary conditions are zero displacement at the bottom boundary
[u(z=H) = 0] and zero shear stress at the top boundary [free end, therefore
τ (z=0) = 0]. The horizontal displacements at any point along the vertical
direction are given by Equation 3.6. This is in fact the solution of the wave
equation in two domains, space, z and time, t.
N
X
u(z, t) = ψn (z)In (t) (3.6)
n=1
79
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
where, ψn (z) describes the mode shape of the nth natural frequency (i.e.
the variation of horizontal displacement in the vertical direction) and is given
by Equation 3.7, whereas In (t) describes the response of a single-degree-of-
freedom (SDOF) structure of the nth mode frequency subjected to the ground
excitation under study (i.e. the variation of displacement in time).
−m/2 " 1−m/2 #
z z
ψn (z) = Jq βn (3.7)
H H
where, βn is the nth root of a period relation, Jq (·) is a Bessel function of
the first kind of order q=m/(2-m), given by Equation 3.8.
∞
(−1)k
q+2k
X x
Jq (x) = (3.8)
k=0 k! Γ(q + k + 1) 2
where, Γ(·) is the gamma function, given by Equation 3.9.
Z ∞
Γ(x) = e−x xn−1 dx (3.9)
0
The amplification of the base motion at the crest of the dam is obtained
from the modulus of the Transfer Function, F (ω) given by Equation 3.10.
q
ωH
2Vs (1+iξ) 1
F (ω) =
ωH
ωH
(3.10)
Γ(q + 1)Jq Vs (1+iξ)
iaz Jq+1 Vs (1+iξ)
√
where, i = −1 is the imaginary unit, az = ρρrsVVsrs is the impedance ratio at
the soil-rock (dam-foundation) interface, ρs and ρr are the mass densities of
the dam soil and foundation rock respectively and Vs and Vsr are the average
shear wave velocities of the dam soil and foundation rock respectively.
Although an attempt is made to capture the dynamic response of dams
and analytical solutions may be obtained, obviously the shear beam method
carries its own limitations which are directly related to its assumptions.
Therefore, its main limitations (of the initial version of the shear beam, such
as Ambraseys (1960a)) are:
80
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
• 1D wave propagation
81
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
sive (in terms of time and computational power) alternative analysis tool as
compared to the more elaborate FE analysis.
The previously-mentioned approach was further extended by Elgamal
et al. (1987b) to consider the longitudinal dimension of the dam, and there-
fore was able to account for the effect of the lateral canyon constraint on the
nonlinear seismic response of earth dams. This approach was applied on the
Santa Felicia dam (Elgamal et al., 1987a) and showed that the analyses tak-
ing into account the longitudinal dimension of the dam yield different results
than the analyses that ignore the longitudal dam direction. Therefore, it
highlighted the importance of modelling the canyon effect for narrower dams
(crest length over height ratio, L/H ≈ 4.7 for Santa Felicia dam).
A two-dimensional theory was developed by Yiagos and Prevost (1991b)
which also considered the two-phase elasto-plastic seismic response of earth
dams. This theory accounted for the presence of water in the dam material
by including the pore water equations in the saturated portion of the dam
and by adopting an effective stress multi-surface elasto-plastic constitutive
model. This theory was later applied on Santa Felicia and Long Valley dams
(Yiagos and Prevost, 1991a) as an elementary attempt to investigate the ef-
fect of coupled analysis on the seismic behaviour of dams. Their computed
crest accelerations and response spectra were compared to the field measure-
ments and exhibited a fair agreement. Although the authors did not directly
compare results of single- and two-phase analyses, they commented that tak-
ing into account the pore water pressure further enhances the simulation of
the dynamic response of a dam.
The shear beam method was also extended to include three-dimensional
(3D) analysis of dams. Abdel-Ghaffar and Koh (1982) developed another
method to analyse 3D isotropic linear elastic dams whose elastic moduli
varied spatially (inhomogeneous) and which were founded on rigid canyons.
The mathematical formulation of their approach was presented in detail and
was then applied on Santa Felicia dam. The computed modes and natural
frequencies of vibration were compared to measured results from full-scale
tests and earthquake response records with reasonable agreement. However,
the authors did not compare any results from 3D and 2D analysis to identify
82
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
83
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 3.4: Geometry of the system considered by Elgamal (1992): (a) physi-
cal problem, (b) UD cross-section, (c) shear stresses acting on an soil element
and (d) geometric configuration of computational model.
the late 1950’s and early 1960’s, and it was soon widely used in studying the
response of earth dams (Clough and Chopra, 1966; Chopra, 1967a; Chopra
et al., 1960).
84
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
where [M], [C] and [K] are the mass, damping and stiffness matrices re-
spectively with entries refering to the nodal points, whereas {ü}, {u̇}, {u}
and {P (t)} are the acceleration, velocity, displacement and external loading
vectors respectively. The stiffness matrix, [K] contains information about the
constitutive behaviour of the materials (see Potts and Zdravković (1999)).
Dynamic finite element analysis involves both spatial discretisation of the
problem into small regions (the finite elements) and temporal discretisation
according to the problem’s loading history. Then, appropriate boundary
conditions (both deformation and hydraulic) are applied and advanced con-
stitutive models may be assigned to different regions and for different time
periods. The acceleration time history is applied at a part of the FE mesh
(normally at the soil-bedrock interface, depending on the stiffness of the two
materials and the impedance contrast).
The advantage of this method compared to the previous methods is the
ability to:
• Model the site more accurately (by advanced boundary conditions) and
include any nearby structures, such as spillway, drainage galleries etc.
85
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
et al. (2010) for the analysis of the Marana Capacciotti dam in Italy. A more
detailed description of the FE method along with relevant recent examples
related to dams is given in Section 3.5.2.
Figure 3.5: FE mesh used for the 2D nonlinear static and dynamic analysis
of Marana Capacciotti dam by Elia et al. (2010).
Finally, it should be noted that the FE method is actually able to fulfill all
three analysis objectives described before: (a) it can provide an assessment
of the stability of the dam structure, through the adoption of constitutive
models that consider soil plasticity, (b) it can provide a reliable calculation
of the permanent displacements of the dam body, through a proper simulation
of degradation of soil stiffness and plastic yielding, and (c) it may offer an
evaluation of the dynamic response of the whole dam through the formulation
of the equation of motion.
86
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 3.6: The finite element mesh used by Clough and Chopra (1966) (B
= 450 ft, H = 300 ft).
of two triangular dams, with slopes H:V being 1.5:1 and 3:1, at the centreline
and at the face of the dam as obtained by FE analysis and compared to the
results of the SB solution (a single result as the SB solution is unique for
any slope angle and the same at the centerline and the face of the dam,
see Section 2.4). It is shown that, the FE results differ significantly between
themselves and with respect to the SB solution (see also Section 3.6). Besides,
the fundamental period calculated with the FE method is smaller than that
of the SB method and as expected is different for different slope angles.
Figure 3.7: Shape of the first natural mode of vibration of a uniform, homo-
geneous dam as obtained by SB (a) and FE (b-e) analyses (Chopra, 1967a).
87
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
88
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 3.8: FE mesh used by Prevost et al. (1985) in the dynamic analysis
of Santa Felicia dam.
Figure 3.9: Hysteretic soil behaviour was modelled by Prevost et al. (1985)
using multi-surface plasticity.
89
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
90
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
fully coupled effective stress formulations following the Biot (1941) theory
for soil-fluid interaction (Zienkiewicz et al., 1999; Kontoe, 2006).
The first attempt to consider both the nonlinearity and the multi-phase
nature of the soil in a seismic analysis of earth dams was the work of Lacy and
Prevost (1987). That was an extention of earlier work (Prevost et al., 1985)
on Santa Felicia dam which considered the coupled dynamic equations for the
response of a two-phase soil system. The resulting maximum acceleration of
the two-phase model was found to be substantially lower than the maximum
of the one-phase model. This was attributed to the additional damping
introduced in the system from the introduction of pore fluid-soil skeleton
interaction.
Fully-coupled effective stress elasto-plastic analyses were carried out by
Muraleetharan et al. (2004) and were compared to centrifuge model tests of
two sand embankments. The constitutive model used in their code, DYSAC2
(Muraleetharan et al., 1988), was based on the bounding surface plasticity
framework (Dafalias and Popov, 1979) which allows plastic strains to occur
for stress points within the bounding surface. The numerical predictions of
accelerations, pore water pressures and deformations compared reasonably
well with the experimental results, although the peak values of accelerations
and pore pressures were underpredicted.
Moreover, Sica et al. (2008) performed coupled dynamic analyses of El
Infiernillo dam in Mexico using their in-house software GEFDYN (Aubry
et al., 1985), adopting an elasto-plastic model (Heujeux, 1985) characterised
by both isotropic and kinematic hardening. The numerical predictions were
compared to the measured response of the dam during the 14/3/1979, 19/9/1985
and 30/5/1990 Mexico seismic events. Their numerical model was calibrated
against the measured response and a good agreement was achieved between
the computed and measured accelerations and reponse spectra of the dam.
The calibration was then followed by a parametric study which aimed to
identify the influence of past loading history on the seismic response of earth
dams.
It was found that when the same seismic event was repeated after the first
event, the computed earthquake-induced crest settlements were smaller dur-
91
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
ing the second time. Moreover, it was found that the difference in the crest
settlements between two successive identical seismic events, was larger for
larger events. The authors commented that the observed difference in the
crest settlements (and therefore the influence of past loading history) was
due to strain hardening of the soil and therefore their study demonstrated
the effect of soil hardening due to past seismic history (a series of consecu-
tive earthquakes) on the computed loss of freeboard. Their conclusions were
drawn due to the adoption of the advanced elasto-plastic model that was
able to capture soil hardening. Therefore, this highlighted the advantages of
using appropriate constitutive models that take account of strain hardening.
However, it should be noted that their study did not directly compare results
from different numerical analyses that consider and ignore strain hardening
and hence it was assumed that the difference in the predicted crest settle-
ments was entirely due to the modelled strain hardening.
Fully coupled dynamic analyses of the Marana Capacciotti dam in Italy
were performed by Elia et al. (2010) employing an advanced elasto-plastic
constitutive model based on multi-surface plasticity. The model for struc-
tured soils, MSS (Kavvadas and Amorosi, 2000), uses two nested elliptic
surfaces in the stress space of the modified Cam-Clay type. This formulation
includes both isotropic and kinematic hardening and is able to describe decay
of stiffness, increase of hysteretic damping and build-up of pore water pres-
sures with shear strain. The FE mesh used is shown in Figure 3.10 and the
seismic behaviour of the dam was analysed using a number of input seismic
motions.
However, because the FE code they employed, DIANA-SWANDYNE II
(Chan, 1988, 1995) could not offer the option of dynamic absorbing bound-
aries, they had to move the lateral boundaries of the FE mesh (Figure 3.10
(a)) further away from the problem of interest (Figure 3.10 (b)). Moreover,
viscous boundaries along the left and right sides of the foundation layer were
simulated by means of two columns of elements characterised by a Rayleigh
damping equal to 25%, in order to minimise wave reflections during the seis-
mic action. This however, resulted in the adoption of a very wide mesh which
generally results in high computational cost.
92
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 3.10: The FE mesh used by Elia et al. (2010): (a) the entire FE
mesh describing the problem and (b) detail of the FE mesh close to the
embankment.
93
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
should be noted that the numerical predictions of that study were not com-
pared to any field data, as such data did not exist and therefore this raises
questions about the reliability of those results. Nevertheless, the authors
commented that their study represents a class A prediction of the dynamic
response of the dam.
Finally, it should also be noted that, as discussed earlier, the authors
attributed the prediction of significant plasticity and additional consolida-
tion settlements to their adopted advanced kinematic hardening constitutive
model. But their study did not compare results from analyses considering
and ignoring kinematic hardening, and therefore kinematic hardening cannot
be identified as the sole reason for achieving those results. If such a compar-
ison was performed, it would better highlight the benefits of such advanced
constitutive models. Therefore, further investigation is required, particularly
using available measured field data, to confirm (or perhaps question) the
results of the study of Elia et al. (2010).
94
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1
∇2 φ(x, y, t) = φ̈(x, y, t) (3.12)
Vp2
where, Vp is the p-wave (acoustic) velocity of water and φ is the velocity
potential, which is related to the pressure, p and the velocity vector, {v}
through Equations 3.13 and 3.14.
{v} = ∇φ (3.13)
p = −ρφ̇ (3.14)
∂p ṗ
= (3.15)
∂n Vp
95
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
96
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
97
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
the SB method forces the entire elastic wave energy to produce shear defor-
mations, whereas the FE method can model the actual dam’s response in
both shear and dilatational deformations.
Tsiatas and Gazetas (1982) compared the results of 2D SB and plane-
strain FE analyses on a number of idealised dam cases. This comparative
study showed that the fundamental frequency of most dams was overesti-
mated by 5% by the SB method, whereas higher percentages have been ob-
served for frequencies corresponding to higher modes (as they involve a higher
contribution of bending deformation). Although the SB displacements can
be considered relatively accurate (within 10% of those resulting from FE
predictions), crest accelerations for flexible dams appeared to be up to 50%
higher than the corresponding FE predictions, due to the so-called ’whip-
lash effect’ (Gazetas, 1987). The latter phenomenon is illustrated in Figure
3.11 which compares the peak values of acceleration in an idealised triangu-
lar dam, subjected to a real earthquake acceleration record, as these result
from the FE (dashed line) and SB (solid line) methods. It was illustrated by
Gazetas (1987) that unrealistically high values of accelerations are obtained
close to the crest of the dam from the SB method (as compared to the more
accurate FE method).
98
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
been found in the literature, the FE method seems to be the most powerful
tool, as it can overcome most of the assumptions of the other methods. In
addition, the SB method can also be used, particularly to provide analytical
relations of simplified problems in the frequency domain (such as periods
and modes of vibration) which could be used for comparison with the FE
predictions.
99
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
∂u(t)
σ(t) = ρVp (3.17)
∂t
100
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
∂v(t)
τ (t) = ρVs (3.18)
∂t
where, σ(t), τ (t) are the normal and shear stresses on the boundary, u(t)
and v(t) are the normal and tangential displacements on the boundary, Vp and
Vs are the compression and shear wave velocities of the medium respectively
and t is the time.
The Cone BC (Kellezi, 1998, 2000) was developed to deal with vibration
problems on the surface of a half-space where the generated waves propagate
in an area that increases with depth. This BC could be considered as an
extension of the Standard Viscous BC and can be physically represented by
a series of normal and tangential sets of dashpots and springs on the FE
mesh boundary. Due to the inclusion of the additional springs, the cone
BC approximates the stiffness of the unbounded soil domain and therefore it
eliminates the permanent movement which occurs with the Viscous BC at low
frequencies (Kontoe, 2006; Kontoe et al., 2009). For plane strain conditions,
the dashpots are defined by Equations 3.17 and 3.18, whereas the springs are
defined by Equations 3.19 and 3.20
Vp2
σ(t) = ρ u(t) (3.19)
2r
Vs2
τ (t) = ρ v(t) (3.20)
2r
where, r is the distance from the boundary node to the source location.
101
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
high, then the bedrock may be assumed to be rigid and hence there is no need
to be discretised. In that case, the location of the bottom boundary is at the
interface of the soil and the bedrock. On the other hand, if the impedance
constrast is not relatively high, then a part of the bedrock should be modelled
too and special techniques may be used to reduce the computational domain,
such as the Domain Reduction Method (DRM) (Bielak et al., 2003; Kontoe,
2006; Kontoe et al., 2008).
As fas as the lateral boundaries are concerned, as illustrated by Figure
3.12, they should be placed at such distance so that the response close to
the boundaries is similar to the free-field response. For this reason, 1D site
response and 2D plane strain analyses should be performed and compared
102
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
103
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
2ρ∞ − 1
αm = = 0.35 (3.22)
ρ∞ + 1
ρ∞
αf = = 0.45 (3.23)
ρ∞ + 1
1
α = (1 − αm + αf )2 = 0.3025 (3.24)
4
1
δ= − αm + αf = 0.6 (3.25)
2
In the analyses carried out in this thesis, the CNL model used was the
Logarithmic (Puzrin and Burland, 2000) model, implemented in ICFEP by
Taborda (2011). Equation 3.26 defines the backbobe curve described by the
Logarithmic model.
R #
E ∗ Gmax
"
∗
J = Ed∗ Gmax 1 − a ln (1 + d ) (3.26)
JL
104
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1 q
J = √ (σ10 − σ20 )2 + (σ20 − σ30 )2 + (σ30 − σ10 )2 (3.27)
6
2 q
Ed = √ (1 − 2 )2 + (2 − 3 )2 + (3 − 1 )2 (3.28)
6
where, σ10 , σ20 and σ30 are the principal effective stresses whereas 1 , 2 and
3 are the principal strains.
Parameters a, JL and R are model parameters described by Equations
3.29 - 3.31. It is worth to note that the maximum shear stiffness of the
soil, Gmax , is a parameter necessary for the calibration of the CNL model.
Therefore, the performance of the model depends on the spatial variation of
Gmax in a boundary value problem.
c(1 + χL ) ln (1 + χL )
R= (3.29)
χL (χL − 1)
χL − 1
a= (3.30)
χL [ln(1 + χL )]R
EdL
χL = Gmax (3.31)
JL
Figure 3.13 shows schemetically the backbone curve of the Logarithmic
model implemented in ICFEP. It is shown that the backbone curve of the
model follows a straight line after a value of the strain, EdL . The value of the
slope (Gimp ) of this straight line is determined by the value of the parameter
c. This has a direct impact on the induced damping by this model. It should
also be noted that there is an option in ICFEP to specify the minimum
value of the shear modulus, Gmin , independently of the value of parameter c
105
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
J*
Gmax
Gimp
JL
c JL
E dL E*d
Figure 3.13: Backbone curve of the Logarithmic CNL model (After Taborda
(2011)).
Damping, ξ
{
c=1
for AND
Gmin = 0
{
c<1
for OR
Gmin > 0
Shear strain, γ
Figure 3.14: Illustration of the influence of parameters c and Gmin on the
damping introduced by the Logarithmic CNL model (After Taborda (2011)).
106
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
specified. Then, the minimum value of the tangent shear modulus used by
the program will be the largest of the two, i.e. max{Gimp , Gmin } and this
value finally dictates the amount of damping introduced (see Figure 3.14).
More details about the formulation and implementation of these CNL models
in ICFEP can be found by Taborda (2011).
The previously described CNL model simulates the pre-yield elastic be-
haviour of the soil. It can be used though in conjuction with an elasto-plastic
model which can simulate the plastic yielding of the material. In this thesis,
the CNL model was coupled with a Mohr-Coulomb failure criterion. The
soil behaviour becomes plastic when the yield surface is engaged. The Mohr-
Coulomb yield surface is described by Equation 3.32 and shown in Figure
3.15.
c0
!
F =J− + p0 g(θ) = 0 (3.32)
tan φ0
where, J is the deviatoric stress (Equation 3.27), p0 is the mean effective
stress (Equation 3.33), θ is the Lode’s angle (Equation 3.34), φ0 is the angle
of shearing resistance and g(θ) defines the shape of the yield surface in the
deviatoric space (Equation 3.35).
J
g(θ)
c'
g(θ)
tanφ'
p'
Figure 3.15: Yield surface of the Mohr-Coulomb model
1
p0 = (σ10 + σ20 + σ30 ) (3.33)
3
107
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
σ 0 − σ30
" !#
1
θ = √ 2 20 −1 (3.34)
3 σ1 − σ30
sin φ0
g(θ) = φ0
(3.35)
cos θ + sin θ√sin
3
A Bωi
ξi = + (3.36)
2ωi 2
where, coefficients A and B are called Rayleigh coefficients and they are
given by Equation 3.37:
A 2ξt ω1 ω2
= (3.37)
B ω1 + ω2 1
where, ξt is the target damping ratio and ω1 and ω2 are the circular
frequencies between which the damping ratio is approximately equal to ξt as
shown in Figure 3.16. Circular frequencies ω1 and ω2 are usualy taken as the
first and third circular frequencies of the system (Zerwer et al., 2002), but
this is strongly problem dependent and should be calibrated for each case
considered (Kontoe et al., 2011).
The values A and B are used to form the damping matrix, [C], of the
dynamic FE equations (Equation 3.11) as described in Equation 3.38:
Where, [ME ] and [KE ] are the elastic mass and stiffness matrices respec-
tively, as they arise in the dynamic FE equations (Equation 3.11). More
information about Rayleigh damping in dams can be found in Woodward
and Griffiths (1996) and about geotechnical problems in general in Kontoe
et al. (2011).
108
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Rayleigh
Rayleigh damping
Damping
10
8
Rayleigh damping
[%]
7
damping
6 Target damping
[%]
ξ5
Damping
t
Rayleigh
3
Mass damping
2
Stiffness damping
1
0
0 ω1 2 3 4 ω
5 6 7
1 Frequency,
Frequency, ω ω 2
Figure 3.16: Illustration of Rayleigh damping and its mass and stiffness
components
3.8 Summary
This chapter dealt with the seismic analysis of dams. Firstly, a review of the
various available methods of seismic analysis of dams was presented explain-
ing their assumptions and describing their potentials and limitations. These
methods fall into three main categories: (a) pseudo-static methods to assess
the stability of the earth dam (i.e. the dam slope), (b) sliding block analy-
sis methods which calculate the permanent displacements and (c) methods
evaluating the dynamic response of the whole dam structure, such as the
shear beam approach. However, advanced numerical techniques, such as the
finite element method are able to satisfy all the analysis objectives: stability,
displacements and dynamic response.
Afterwards, a description of the evolution of the numerical methods of
analysis applied to dams was presented. Numerical methods were formu-
lated in order to model the nonlinear material behaviour of the dam materi-
als, which cannot be calculated analytically. Various methods were developed
over the years and these range from the simple numerical shear beam method
up to the sophisticated nonlinear coupled dynamic analysis including consid-
eration of reservoir-dam interaction effects.
Finally, taking into consideration the plethora of information found in the
109
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
110
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Chapter 4
NUMERICAL MODELLING
OF THE RESERVOIR
DOMAIN
4.1 Introduction
The various methods and approaches available for modelling the hydrody-
namic pressures on dams and dynamic reservoir-dam interaction in general
were described in Section 3.5.4. However, there are still issues that need to be
investigated so that useful and reliable numerical analyses can be performed.
These issues are related to both geometric modelling, such as absorbing
boundary conditions and material modelling such as water compressibility
and shear stiffness.
This chapter describes work carried out on modelling the hydrodynamic
pressures on dams during earthquakes and general dynamic loads using finite
elements. The reservoir water is descritised under plane-strain conditions
using displacement-based eight-noded quadrilateral solid elements, which are
the same as the elements used to discretise the dam and the foundation soil.
111
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
4.2 Background
The hydrodynamic pressures on a rigid (undeformable) dam with a vertical
upstream face under harmonic loading were investigated in the fundamental
work of Chopra (1967a) (see also Sections 2.5.4 and 2.5.5). In this analyt-
ical approach, the pressures were expressed in the frequency domain using
transfer functions for various values of loading frequency. Morever, the re-
sults of Chopra (1967a) for incompressible water agree with Westergaard
(1933). Besides, it was shown that the pressure response due to an arbitrary
load (such as earthquake) may be obtained from the unit impulse response.
Therefore, the response under a number of impulses of different frequencies
is superimposed with the aid of the convolution integral. The horizontal hy-
drodynamic pressure, px , distribution on the dam face over time, t, under a
general acceleration load, üg (t), is given by Equation 4.1.
4γw Vp +∞
X (−1)n−1 Z t
px (y, t) = cos (λn y) üg (τ )Jo [λn Vp (t − τ )]dτ (4.1)
πg n=1 2n − 1 0
where, y is the vertical distance from the bottom of the dam, t is the
time, γw is the unit weight of water, Vp is the acoustic (p-wave) velocity of
water, λn is the nth wave length, üg (t) is the ground acceleration and J0 (·) is
a Bessel Function of the first kind of order 0.
The problem under study is shown graphically in Figure 4.1. It consists
of a dam (1) which retains a large volume of water in the reservoir (2) and
rests on the ground (3) which serves as the foundation of the dam. Under
seismic or general dynamic conditions, the reservoir induces hydrodynamic
pressures on the upstream (US) dam face, A-B.
As far as efficient and economic numerical modelling of these hydrody-
namic pressures is concerned, the US reservoir is truncated at some distance
from the dam, (C-D). Likewise, the foundation soil is truncated at some
distance from the dam and reservoir, E-F-G-C. On both soil and reservoir
boundaries, special absorbing boundary conditions (BCs) need to be applied
so that reflection of outgoing waves is avoided.
112
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
A D
1 2
B C
E
3
F G
113
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
same as Equation 4.1 of Chopra (1967a). The analytical solution for the
hydrodynamic pressure time history at the base of the dam (bottom left
corner of the reservoir in Figure 4.2) is shown in Figure 4.4, whereas the
distribution of the maximum hydrodynamic pressure over the entire upstream
face is shown in Figure 4.5.
Figure 4.3: The ramp acceleration load applied to the dam-reservoir system
of Küçükarslan et al. (2005)
114
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1.2
Normalised Hydrodynamic Pressure, Pdyn/ρ α H []
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.4: The analytical solution for the hydrodynamic pressures at the
base of the dam due to the ramp acceleration using the solution of Chopra
(1967a).
0.8
0.7
Normalised Height, y/H []
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
Normalised Hydrodynamic Pressure, Pdyn/ρ α H []
Figure 4.5: The analytical solution for the distribution of the maximum
hydrodynamic pressure over the entire upstream face of the dam due to the
ramp acceleration using the solution of Chopra (1967a).
115
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 4.6: The hydrodynamic pressure time-history at the base of the dam
as obtained by Küçükarslan et al. (2005).
Figure 4.7: The distribution of the peak hydrodynamic pressure on the up-
stream face of the dam as obtained by Küçükarslan et al. (2005).
116
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
to model the reservoir domain, the stiff dam and the foundation parts. The
maximum size of the elements was 4.5m which is equal to 1/40 of the reservoir
height.
Table 4.1: Geometry of the FE meshes considered for the ramp acceleration
loading
The load was applied along the left (E-F in Figure 4.8) and the bottom (F-
G) boundary as prescribed values of acceleration in the horizontal direction.
The displacement in the vertical direction on both boundaries was prescribed
to be equal to zero. Different boundary conditions (BCs) were examined for
the upstream reservoir boundary (C-D), and these are:
• Free, i.e. zero displacement and stress in the vertical and horizontal
directions respectively
117
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
E
A D
reservoir
H
stiff material
B C
T F
G
T L
Figure 4.8: Geometry of the dam-reservoir system used for the ramp accel-
eration loading.
The water was assumed to behave as a linear elastic material with a bulk
modulus, Kw = 2.2 · 106 kPa (the exact value for water) and a nominal value
of shear modulus, Gw = 100 kPa. The reason for assigning a very small value
for the shear modulus is to avoid shear wave propagation in the water domain,
as physically the shear stresses carried by water are negligible. The dam and
foundation were modelled to behave as a “rigid” body (i.e. extremely stiff),
therefore the bulk modulus assigned is K = 108 · Kw = 2.2 · 1014 kPa. The
118
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
viscosity, CH (Equation 4.3) of the dashpots and the stiffness of the springs,
KH (Equation 4.4) were obtained from the elastic properties of the material
to which they were applied, i.e. water. No damping was assigned to the
reservoir or the dam domains.
CH = A(r)ρVp (4.3)
A(r) 2
KH = ρVp (4.4)
r
where, ρ is the mass density of water, r is the distance from the source of
the loading to the boundary and A(r) is the cross-sectional area of the cone.
The value of the distance parameter is chosen to be equal to the length
of the reservoir, i.e. r = L and A(r) was specified to be equal to the dis-
tance between the boundary nodes (see Section 3.7). More details about the
Viscous and Cone BCs as implemented in ICFEP, can be found in Kontoe
(2006).
Figure 4.10 shows the hydrodynamic pressure time-histories at the base
of the dam (Point B in Figure 4.8) for the Free BC (a) for all the five lengths
of the FE mesh along with the analytical solution. Likewise Figures 4.11,
4.12 and 4.13 show the corresponding time-histories of the Viscous (b), the
Cone (c) and the Viscous with hydrostatic pressure distribution (d) BCs
respectively. The hydrodynamic pressures are normalised with respect to
ρaH, where ρ is the mass density of water, a = 1 m/s2 is the maximum value
of the acceleration and H is the height of the reservoir.
Figures 4.10 - 4.13 show that none of the four models adequately re-
produced the hydrodynamic pressures at the base of the dam. The closest
results appear to be obtained from MESH3 (see Table 4.1), i.e. L = 1260 and
L/H = 7 with the Viscous and Cone boundary conditions. However, for all
four cases of the BCs considered, the first half cycle appears to be predicted
accurately and then depending on the distance of the lateral boundary, the
hydrodynamic pressures tend to increase with time for the shorter meshes
(L/H = 3 and 5.33) and to decrease for the longer meshes (L/H = 10 and
13.33).
119
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Numerical − L/H = 7
4 Numerical − L/H = 10
Numerical − L/H = 13.33
3
−1
−2
−3
−4
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.10: Hydrodynamic pressure at the base of the dam by using the
“Free” BC for different lengths of the reservoir
Numerical − L/H = 7
4 Numerical − L/H = 10
Numerical − L/H = 13.33
3
−1
−2
−3
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.11: Hydrodynamic pressure at the base of the dam by using the
“Viscous” BC for different lengths of the reservoir
120
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Numerical − L/H = 7
4 Numerical − L/H = 10
Numerical − L/H = 13.33
3
−1
−2
−3
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.12: Hydrodynamic pressure at the base of the dam by using the
“Cone” BC for different lengths of the reservoir
Numerical − L/H = 7
4 Numerical − L/H = 10
Numerical − L/H = 13.33
3
−1
−2
−3
−4
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.13: Hydrodynamic pressure at the base of the dam by using the
“Viscous + hydrostatic stress” BC for different lengths of the reservoir
121
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Moreover, for the Free BC (Figure 4.10), the fluctuations in the hydrody-
namic pressure tend to have a larger amplitude than the analytical solution
(dashed line). For the Viscous and Cone BCs (Figures 4.11 and 4.12) the
amplitude of these fluctuations is similar and comparable to the amplitude of
the analytical solution. In fact, the behaviour of the latter two BCs is almost
identical. Finally, for the “Viscous + hydrostatic stress” BC (Figure 4.13),
one may observe some spurious peaks in the hydrodynamic pressures which
happen at different times for the different lengths of the mesh and which
decrease with time. Furthermore, it seems that the time at which these high
peaks occur increases with the length of the mesh and this can be attributed
to waves reflected back from the upstream reservoir boundary. It seems that
the specification of a constant (hydrostatic) pressure at the boundary results
in a reflection of the waves. It is therefore concluded that both the Free and
Viscous + hydrostatic pressure BCs are not appropriate for modelling the
upstream boundary of the reservoir.
122
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
terface elements eliminated the general increase and decrease of the pressures
with time, as compared to the analytical solution.
As before, for the case of the Free BC (Figure 4.14), it is clear that the
fluctuations of the numerical results are significantly larger than the analyt-
ical solution (dashed line). A better agreement appears for the cases of the
Viscous (Figure 4.15) and Cone (Figure 4.16) BCs for which the amplitudes
of the pressures are very close to the analytical solution. Finally, the results
for the “Viscous + hydrostatic stress” BC (Figure 4.17) still show spurious
peaks in the pressure, although the values tend to get closer to the analytical
solution with time.
The effect of the mesh length for the Viscous BC is shown in more detail
in Figures 4.18 - 4.22. Clearly, the comparison of the numerical results with
the analytical solution is improved. Moreover, the numerical response starts
to deviate from the analytical one earlier (i.e. at smaller values of the time,
t) for the shorter meshes. The differences are mainly a slight increase in the
amplitude and a small shortening in the period of the fluctuations. This is
expected as any waves reflected at the boundary require less time to travel
back to the dam if the distance between the dam and the boundary is small.
Note that the results for the cases of the Cone BC plot on top of the results
of the Viscous BC.
Finally, Figure 4.23 shows the distribution of the peak hydrodynamic
pressure on the face of the dam for all 5 lengths of the mesh and how these
compare with the analytical result. It is apparent that the distributions from
all the analyses are extremely similar as there are no visible differences. This
confirms that the model used (Viscous or Cone BC with interface elements)
is adequate to describe the hydrodynamic pressures on a rigid dam with a
vertical upstream face. Again, note that the results for the cases of the Cone
BC plot on top of the results of the Viscous BC.
Considering the previous results (especially Figures 4.15 and 4.16), it may
be concluded that the best comparison with the analytical solution can be
obtained using a FE mesh with L/H = 13.33 (Figure 4.22). Of course, the
ramp acceleration may not be the most suitable load to be used to determine
an appropriate distance of the BC. However, the study carried out so far
123
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Numerical − L/H = 7
Numerical − L/H = 10
Numerical − L/H = 13.33
1
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.14: Hydrodynamic pressure at the base of the dam by using the
“Free” BC with interface elements for different lengths of the reservoir
Numerical − L/H = 7
Numerical − L/H = 10
Numerical − L/H = 13.33
1
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.15: Hydrodynamic pressure at the base of the dam by using the
“Viscous” BC with interface elements for different lengths of the reservoir
124
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Numerical − L/H = 7
Numerical − L/H = 10
Numerical − L/H = 13.33
1
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.16: Hydrodynamic pressure at the base of the dam by using the
“Cone” BC with interface elements for different lengths of the reservoir
Numerical − L/H = 10
2 Numerical − L/H = 13.33
dyn
Normalised Hydrodynamic Pressure, P
−1
−2
−3
−4
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.17: Hydrodynamic pressure at the base of the dam by using the “Vis-
cous + hydrostatic stress” BC with interface elements for different lengths
of the reservoir
125
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.18: Hydrodynamic pressure at the base of the dam by using the
“Viscous” BC with interface elements for L/H = 3
1
Normalised Hydrodynamic Pressure, P
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.19: Hydrodynamic pressure at the base of the dam by using the
“Viscous” BC with interface elements for L/H = 5.33
126
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1.2
Normalised Hydrodynamic Pressure, Pdyn/ρ α H []
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.20: Hydrodynamic pressure at the base of the dam by using the
“Viscous” BC with interface elements for L/H = 7
1.2
Normalised Hydrodynamic Pressure, Pdyn/ρ α H []
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.21: Hydrodynamic pressure at the base of the dam by using the
“Viscous” BC with interface elements for L/H = 10
127
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1.2
Normalised Hydrodynamic Pressure, Pdyn/ρ α H []
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6
Time, t [sec]
Figure 4.22: Hydrodynamic pressure at the base of the dam by using the
“Viscous” BC with interface elements for L/H = 13.33
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
Normalised Hydrodynamic Pressure, Pdyn/ρ α H []
128
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
suggests that both the Viscous and the Cone BCs can perform quite well
as absorbing BCs for the reservoir. In addition, interface elements with
appropriate values for the shear and normal stiffnesses should be introduced
at the interface between the reservoir and the solid domains.
2π πVp
ω1 = = (4.5)
T1 2H
where, T1 = 4H/Vp is the fundamental frequency of the reservoir, Vp is
the acoustic (p-wave) wave velocity of the water and H is the height of the
reservoir.
The hydrodynamic force, Fdyn is obtained from the integration of the
normal pressures on the US face of the dam. The hydrostatic force, Fst =
1
2
ρgH 2 , where ρ is the mass density of the water, g is the acceleration of
gravity and H is the height of the reservoir.
129
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
3.5
2.5
st
/F
dyn
2
gF
1.5
0.5
0
0 1 2 3 4 5 6
Frequency Ratio, Ω = ω/ω1
130
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
2.5
st
/F
dyn
2
gF
1.5
0.5
0
0 1 2 3 4 5 6
Frequency Ratio, Ω = ω/ω1
2.5
g Fdyn/Fst
1.5
0.5
0
0 1 2 3 4 5 6
Frequency Ratio, Ω = ω/ω1
131
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
(ω/ω1 = 1, 3, 5, ...) are not included as the response of the reservoir was
transient and the hydrodynamic pressures did not reach a steady state.
As it may be observed from both Figures 4.25 and 4.26, a generally good
agreement is obtained for the whole spectrum of frequencies, for both BCs.
In fact, their performance is almost identical. This suggests that both the
Viscous and the Cone BCs can be applied on the upstream boundary of the
reservoir in order to model the reservoir’s response due to a harmonic load
of a wide range of frequencies.
132
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
2.5
st
/F
dyn
2
gF
1.5
0.5
0
0 1 2 3 4 5 6
Frequency Ratio, Ω = ω/ω1 compressible
the hydrodynamic pressures in the low frequency range that occur due to
resonance between the loading and the reservoir water, or overestimate the
hydrodynamic pressures for high frequency loads.
133
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
used to model the soil domain so that wave propagation can be adequately
modelled. These studies expressed the size of the element in terms of the
wavelength, as the full wavelength has to be appropriately modelled. For
this investigation, a FE mesh of height H=90m and length L=900m (i.e.
L/H = 10) was employed. Eight-noded isoparametric quadrilateral elements
(see also Potts and Zdravković (1999)) were used as shown in Figure 4.29.
Vp
λ= (4.7)
f
The fundamental frequency of vibration of a reservoir of height H is given
by (Chopra, 1967a):
1 Vp
f1 = = (4.8)
T1 4H
134
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
4H 4H
λ= = (4.9)
f /f1 ω/ω1
Therefore, for a harmonic load with frequency such that ω/ω1 = f /f1 = 4,
the wavelength is given by:
ω 4H
λ =4 = =H (4.10)
ω1 4
A harmonic load was applied with frequency ratio, ω/ω1 = 4 and the
range of element sizes d explored is listed in Table 4.2. The reservoir was
considered compressible and the Viscous BC was applied on the upstream
reservoir boundary.
Table 4.2: Size of the finite elements considered for the harmonic acceleration
loading
135
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.9
0.8
λ/d = 5
Normalised Height, y/H []
0.7
λ/d = 1
0.6 λ/d = 2 λ/d = 4
λ/d = 3
0.5 λ/d = 4
λ/d = 1 λ/d = 3
λ/d = 5
0.4 λ/d = 2
λ/d = 6
0.3 λ/d = 8
λ/d = 10
0.2 λ/d = 12
λ/d = 15
0.1
λ/d = 20
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
Normalised Hydrodynamic Pressure, Pdyn(y)/Pdyn(0) []
are considered to be the most accurate because a large number of nodes has
been provided to model the wavelength.
As it may be observed from Figure 4.30, the peak hydrodynamic pressure
obtained using elements of size d equal or smaller than a fifth of the wave-
length, λ, deviate significantly from the most accurate prediction (λ/d = 20),
whereas the results from λ/d = 6 and 8, are reasonably close. It is therefore
suggested that for modelling dynamic problems involving water using two-
dimensional eight-noded solid elements, the elements representing the water
domain should be smaller or equal to a fifth of the wavelength.
136
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 4.31.
0.03
0.02
Acceleration, a, [g]
0.01
−0.01
−0.02
−0.03
−0.04
−0.05
0 1 2 3 4 5 6
Time, t [sec]
137
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
the corresponding values for the lengths of the reservoir are listed in Table
4.4. The maximum element size in all the analyses was 5m which is equal to
1/12 of the height of the reservoir. For the flexible dam analyses, material
properties were assigned to the dam that correspond to a shear wave velocity,
Vs = 300 m/s, which is considered to be a typical value for earth materials.
These properties are: the shear modulus, Gd = 180000 kPa, the unit weight,
γd = 20 kN/m3 and the Poisson’s ratio, νd = 0.3.
H
T
W B L
Figure 4.32: Geometry of the dam-reservoir system considered for the inves-
tigation of the effect of the slope
MESH Slope
B/H
A 0
B 1/3
C 1/2
D 1
E 2
F 3
Figures 4.33 to 4.38 show the distribution of the peak hydrodynamic pres-
sure on the upstream face of a stiff dam for all 6 slope angles and all 10 mesh
lengths considered. The hydrodynamic pressure values are normalised with
respect to ρam H, where ρ is the mass density of water, am = 0.04g is the
maximum value of the applied acceleration and H is the height of the reser-
voir. Likewise, Figures 4.39 to 4.44 show the corresponding distribution of
the peak hydrodynamic pressure for a flexible dam. The results from the
longer reservoir (i.e. the mesh corresponding to L/H = 15, according to
Table 4.4) are considered to be the most accurate because the longer the
138
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
MESH Length
L/H
1 0.5
2 1
3 2
4 3
5 4
6 5
7 8
8 10
9 12
10 15
distance of the upstream boundary, the smaller the effect on the hydrody-
namic pressures on the dam. Considering the extreme cases, the following
observations can be made:
For a stiff dam with a vertical upstream face, i.e. B/H = 0 (Figure 4.33),
the results for the peak hydrodynamic pressure on the dam for L/H = 0.5,
1, 2, 3, 4, 5 differ from the most accurate solution for L/H = 15. The
results for L/H > 5 plot very close to the most accurate solution. Therefore,
when modelling the pressures on a stiff dam with a vertical upstream face
the upstream boundary (and hence the relevant BC) should be placed at a
distance of more than 5 times the height of the reservoir.
For a stiff dam with a gentle sloped upstream face, like B/H = 3 (Figure
4.38), the results for the peak hydrodynamic pressure on the dam for L/H =
0.5, 1, 2, 3 differ from the most accurate solution for L/H = 15. The results
for L/H > 3 plot be very close to the most accurate solution. Therefore,
when modelling the pressures on a stiff dam with a sloped upstream face
the upstream boundary (and hence the relevant BC) should be placed at a
distance of more than 3 times the height of the reservoir.
For a flexible dam with a vertical upstream face, i.e. B/H = 0 (Figure
4.39), the results for the peak hydrodynamic pressure on the dam for L/H =
0.5, 1, 2 differ from the most accurate solution for L/H = 15. The results
for L/H > 2 seem to be very close to the most accurate solution. Therefore,
139
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.7 L/H = 5
L/H = 8
0.6 L/H = 10
L/H = 12
0.5 L/H = 15
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
0.7 L/H = 5
L/H = 8
0.6 L/H = 10
L/H = 12
0.5 L/H = 15
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
140
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.7 L/H = 5
L/H = 8
0.6 L/H = 10
L/H = 12
0.5 L/H = 15
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
0.7 L/H = 5
L/H = 8
0.6 L/H = 10
L/H = 12
0.5 L/H = 15
0.4
0.3
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
141
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.7 L/H = 5
L/H = 8
0.6 L/H = 10
L/H = 12
0.5 L/H = 15
0.4
0.3
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
0.7 L/H = 5
L/H = 8
0.6 L/H = 10
L/H = 12
0.5 L/H = 15
0.4
0.3
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
142
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.7 L/H = 5
L/H = 8
0.6 L/H = 10
L/H = 12
0.5 L/H = 15
0.4
0.3
0.2
0.1
0
−0.1 0 0.1 0.2 0.3 0.4 0.5 0.6
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
0.7 L/H = 5
L/H = 8
0.6 L/H = 10
L/H = 12
0.5 L/H = 15
0.4
0.3
0.2
0.1
0
−0.1 0 0.1 0.2 0.3 0.4 0.5 0.6
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
143
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.7 L/H = 5
L/H = 8
0.6 L/H = 10
L/H = 12
0.5 L/H = 15
0.4
0.3
0.2
0.1
0
−0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
0.9
0.8
Normalised Height, y/H []
0.7
0
−0.1 0 0.1 0.2 0.3 0.4 0.5 0.6
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
144
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.9
0.8
Normalised Height, y/H []
0.7
0
−0.1 0 0.1 0.2 0.3 0.4 0.5 0.6
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
0.9
0.8
Normalised Height, y/H []
0.7
0
−0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1 1.2
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
145
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
when modelling the pressures on a flexible dam with a vertical upstream face
the upstream boundary (and hence the relevant BC) should be placed at a
distance of more than 2-3 times the height of the reservoir.
For a flexible dam with a gentle sloped upstream face, like B/H = 3
(Figure 4.44), the results for the peak hydrodynamic pressure on the dam
for L/H = 0.5 and 1 differ from the most accurate solution for L/H =
15. The results for L/H > 1 seem to be very close to the most accurate
solution. Therefore, when modelling the pressures on a stiff dam with a
sloped upstream face the upstream boundary (and hence the relevant BC)
should be placed at a distance of more than 2 times the height of the reservoir.
Considering the above-mentioned observations it is concluded that for
a sloped and flexible dam, which is generally the case for earth dams, the
distance of the absorbing BC may be reduced down to twice or three times
the height of the reservoir, i.e. L/H = 2 - 3.
Combining the results for different slope angles (B/H), it may be ob-
served that for the case of the stiff dam (Figures 4.33 - 4.38), the distribution
of the hydrodynamic pressures over the height of the dam has a similar shape
for all slope angles. This may be observed from Figure 4.45 which includes
the distribution of the peak hydrodynamic pressure on the face of the stiff
dam for all values of B/H considered and for L/H = 15. Moreover, it is
shown that the value of the hydrodynamic pressures reduces with increas-
ing slope angle (B/H), and this confirms the observations of Zangar (1952)
who first investigated the hydrodynamic pressures on sloped dams (see also
Section 2.5).
As far as the case of the flexible dam is concerned (Figures 4.39 - 4.44),
the computed hydrodynamic pressures generally tend to be smaller than the
corresponding pressures for the stiff dam, but not in all the cases. This can
be observed in Figure 4.46 which compares the distribution of the peak hy-
drodynamic pressure on the face of the flexible dam for all values of B/H
considered and for L/H = 15. This is confirmed by centrifuge experiments of
model dams by Saleh and Madabhushi (2010) who observed smaller hydro-
dynamic pressures on a flexible model dam. Besides, the distributions of the
peak hydrodynamic pressures on the dam for different slope angles (B/H) do
146
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.6
0.5
0.4
0.3
0.2
0.1
0
−0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
0.6
0.5
0.4
0.3
0.2
0.1
0
−0.5 0 0.5 1
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
147
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.8
0.7
Normalised Height, y/H []
0.6
0.5
0.4
0.3
0.2
0.1
0
−0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
0.8
0.7
Normalised Height, y/H []
0.6
0.5
0.4
0.3
0.2
0.1
0
−0.5 0 0.5 1
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
148
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
not have the same shape. The negative values of the hydrodynamic pressure
imply that the deformation of the dam is dragging the reservoir and there-
fore reduces the total water (hydrostatic + hydrodynamic) pressure (which
is still positive, i.e. compressive). When comparing the values of the pres-
sures for different slope angles, in contrast to the stiff dam, the pressures do
not appear to reduce for increasing slope angles. Figure 4.47 compares the
distributions of the peak hydrodynamic pressures for both stiff and flexible
vertical dams (B/H = 0). Likewise, 4.48 compares the corresponding dis-
tributions for sloped dams (B/H = 3). It is shown that the hydrodynamic
pressures on flexible dams may be smaller or larger than those on stiff dams.
The latter observation may be attributed to the flexibility of the dam
which is a new variable in the study. It is believed that the response of
the reservoir (and hence the hydrodynamic pressures) are affected by the
response of the flexible dam, which is in return affected by the reservoir. This
is the so called reservoir-dam interaction phenomenon (see Section 2.4.5). In
order to better understand and explain the mechanisms of this phenomenon
one has to consider the natural modes of vibration of a dam (see Section 2.4)
and how these are compared to the modes of the reservoir. This of course
requires further study which is presented later in Chapter 5.
149
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
the interface elements and the Viscous BC for a stiff dam with a vertical
upstream face under the acceleration load shown in Figure 4.31 was analysed
with different values of Gw as listed in Table 4.5.
Table 4.5: Values of the Shear Modulus, Gw assigned to the reservoir water
elements
150
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.9
Gw = 10% Kw
0.8 Gw = 5% Kw
Gw = 1% Kw
Normalised Height, y/H []
0.7
Gw = 0.5% Kw
0.6
Gw = 0.1% Kw
0.5 Gw = 0.05% Kw
G = 0.01% K
0.4 w w
Gw = 0.005% Kw
0.3
Gw = 0.001% Kw
0.2 G = 0.0005% K
w w
Gw = 0.0001% Kw
0.1
G = 0.00005% K
w w
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
0.45 G = 10% K
w w
0.4 Gw = 5% Kw
Gw = 1% Kw
Normalised Height, y/H []
0.35
Gw = 0.5% Kw
0.3 Gw = 0.1% Kw
0.25 Gw = 0.05% Kw
G = 0.01% K
0.2 w w
G = 0.005% K
w w
0.15 Gw = 0.001% Kw
0.1 Gw = 0.0005% Kw
Gw = 0.0001% Kw
0.05
G = 0.00005% K
w w
0
1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
Normalised Hydrodynamic Pressure, P /ρa H []
dyn m
151
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
152
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Table 4.6: Values of the damping, ξw assigned to the reservoir water elements
CASE Damping
ξw
[%]
1 20
2 10
3 5
4 2
5 1
6 0.5
7 0.2
8 0.1
1.4
1.2
Fourier Amplitude
0.8
0.6
0.4
0.2
0
−1 0 1
10 10 10
Frequency, f [Hz]
153
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.7
ξ = 0.5%
0.6 ξ = 0.2%
ξ = 0.1%
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
0.35
ξ = 0.5%
0.3 ξ = 0.2%
ξ = 0.1%
0.25
0.2
0.15
0.1
0.05
0
1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
Normalised Hydrodynamic Pressure, Pdyn/ρamH []
154
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
20, 10, 5 and 2% differ considerably from the most accurate solution (ξw =
0.1%), as they result in smaller pressures. This is not surprising as the high
values of ξw damp significantly the response. Additionally, for ξw < 1%
the results compare very well with the accurate solution. It is therefore
recommended that the value of damping, ξw assigned to water should be
kept close to 1%.
4.6 Conclusions
This chapter describes the investigation carried out in order to numerically
model the reservoir hydrodynamic pressures on dams. The main aim was
to establish an approach/methodology to appropriately model the reservoir
with displacement-based solid elements in FE analysis. The same approach
may be used in analyses of other waterfront structures such as quay walls.
The main conclusions of this study may be summarised as follows:
• The size of the reservoir elements should be smaller than a fifth of the
acoustic (p-wave) wavelength, λ of the water.
155
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
156
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Chapter 5
DYNAMIC
RESERVOIR-DAM
INTERACTION
5.1 Introduction
The dynamic behaviour of a dam with a full reservoir is different from the
behaviour of a dam with an empty reservoir. This is because the motion of the
reservoir affects the motion of the dam and vice-versa. This phenomenon is
called dynamic reservoir-dam interaction and it could be catastrophic in cases
of resonance, i.e. when the two domains (dam and reservoir) are vibrating
in phase.
Previous studies of reservoir-dam interaction, such as Chopra (1967b,
1968) and Hall and Chopra (1980, 1982a,b,c) (see Section 2.4.5 for a more
detailed discussion) showed that the effects of reservoir-dam interaction are
more pronounced in concrete dams than in earth dams. The effects of this in-
teraction are mainly concentrated on (a) the fundamental period of vibration
of the reservoir-dam system and (b) the magnitude of the dynamic response
of the dam. It is believed that reservoir-dam interaction firstly causes the
reservoir-dam system to soften, i.e. its fundamental period elongates, and
secondly it alters the dam’s response, i.e. it amplifies or de-amplifies the seis-
157
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
mic motion and results in larger or smaller accelerations at the dam crest.
This chapter describes work carried out as part of this research project
to further investigate reservoir-dam interaction and its influence on the dy-
namic acceleration response of dams. The main aim is to examine the effects
of reservoir-dam interaction for different dynamic characteristics of the load,
the dam and the reservoir. In this study, two types of dams are considered
for this purpose. Firstly, a rectangular cantilever dam is investigated which
has a vertical upstream face (and a small mass because of its slender ge-
ometry). Secondly, a trapezoidal wide earth dam is considered which has
a sloped upstream face and a larger mass because of its wide trapezoidal
geometry. The response of both dams is investigated under harmonic and
random (seismic) acceleration loads. In all the performed analyses, both the
dam and the reservoir rest on a stiff undeformable ground and therefore the
effects of reservoir-dam-foundation interaction are not considered.
1 2
A
3
Figure 5.1: Geometry of the dam-reservoir-foundation system: (1) Dam, (2)
Reservoir, (3) Foundation.
The critical question is how the motion of each domain (dam structure
and reservoir) will be affected by the motion of the neighbouring domain, i.e.
what is the effect of reservoir-dam interaction, especially on the accelerations
158
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
of the dam crest (point B). Therefore, in this investigation, acceleration input
motions were applied at the dam base (A) and the accelerations at the crest
of the dam (B) were monitored and examined for different cases, with a full
and an empty reservoir.
E A D
Reservoir
H Flexible dam
Stiff foundation
F B C
T G H
W L
Figure 5.2: Geometry of the cantilever dam considered.
159
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Damping, ξ, of the Rayleigh type (see Section 3.7.5) was specified in the
dam with a value of 5%. The values for the two natural circular frequencies
of Rayleigh damping, ω1 and ω2 , were taken as equal to the fundamental
circular frequency of the dam, ωd , and the circular frequency of the harmonic
load, ω (Equation 5.1). The foundation was modelled to behave as rigid
(i.e. extremely stiff) and a high value of the bulk modulus was assigned,
K = 108 · Kw = 2.2 · 1014 kPa (where, Kw = 2.2 · 106 kPa is the bulk modulus
of water).
The applied boundary conditions include: zero displacements specified in
the vertical direction along the bottom boundary (G-H), whereas prescribed
values of harmonic acceleration were specified in the horizontal direction,
as given by Equation 5.1. Different values of circular frequency, ω, were
considered, for constant amplitude a0 = 1 m/s2 and for 40 cycles, in order
to reach a steady-state response.
160
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
12
Amplification, |F|
10
0
0 1 2 3 4 5 6 7 8
Frequency ratio, ω/ωd
161
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Vp
ωr = 0.25 (5.3)
H
Interface elements were placed along the interface of the dam and the
reservoir (A-B) and the interface of the reservoir and the foundation (B-
C). The values assigned for the normal and shear stiffness of the interface
elements were KN = 108 kN/m and KS = 1 kN/m respectively, as discussed
in Chapter 4.
The same boundary conditions as in the case of the dam with an empty
reservoir were applied on the bottom boundary (G-H). The horizontal bound-
ary condition applied on the upstream reservoir boundary (C-D) was the
standard viscous boundary (see Section 3.7) and the values of the viscosities
of the dashpots were obtained from the elastic properties of the material ad-
jacent to the boundary to which the dashpots were applied, i.e. water (see
also Chapter 4). Finally, zero displacements were specified in the vertical
direction on the upstream reservoir boundary.
Figure 5.4 shows the amplification of the accelerations at the dam crest,
|F | with respect to the ratio of the circular frequency of the harmonic load
to the fundamental circular frequency of the dam, ω/ωd , for various values
of the ratio of the fundamental circular frequency of the reservoir to the
fundamental circular frequency of the dam, ωr /ωd . On the same figure, the
corresponding amplification of a dam with an empty reservoir (Figure 5.3) is
plotted with a dashed line for comparison.
The amplification spectrum (|F | − ω/ωd ) follows a similar trend for all
the values of the examined ωr /ωd ratio. Namely, there are generally two peak
162
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
ωr/ωd
2
0
78
6
ω/ωd
5
4
2 3
1
0
5
0
25
20
15
10
Amplification, |F|
163
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
20
Amplification, |F|
15
10
0
0 1 2 3 4 5 6 7 8
Frequency ratio, ω/ωd
Figure 5.5: Amplification spectrum of the cantilever dam with a full (ωr /ωd =
1) and an empty reservoir.
3.5
First Second
mode mode
3
Frequency ratio, ωr/ωd
2.5
1.5
Full reservoir
1
Empty reservoir
0.5
0
0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Frequency ratio, ω/ωd
Figure 5.6: Amplification peaks for the cantilever dam with a full and an
empty reservoir.
164
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
20
Amplification, |F|
15
10
5 Full reservoir
Empty reservoir
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Frequency ratio, ωr/ωd
14 Full reservoir
Empty reservoir
12
Amplification, |F|
10
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Frequency ratio, ωr/ωd
165
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
values of the amplification as observed in Figure 5.3 for the empty reservoir
case. However, the magnitude of the peak values of the amplification and the
value of the ω/ωd ratio at which they occur is different for different values of
the ωr /ωd ratio.
Moreover, it is shown that the maximum value of amplification occurs in
the region where ω/ωd ≈ ωr /ωd ≈ 1, i.e. where ω ≈ ωd ≈ ωr . This is due
to reasonance between the harmonic load, the dam and the reservoir. There
are also large values of amplification for ω/ωd ≈ ωr /ωd (shown diagonally in
the figure), i.e. where ω ≈ ωr . The value of amplification for the latter case
(ω ≈ ωr ) can be larger than the amplification corresponding to the second
mode of vibration (i.e. close to ω/ωd ≈ 4.7). However, it seems that in some
cases (e.g. close to ω/ωd ≈ ω/ωr ≈ 3) there is a combined effect of (a) ω ≈
ωr and (b) the second natural mode of vibration of the dam.
Figure 5.5 shows a comparison of the amplification spectrum (|F | - ω/ωd )
for a dam with an empty reservoir (dashed line) and a dam with a full
reservoir (solid line), for which ωr /ωd = 1, i.e. the fundamental circular
frequency of the reservoir, ωr , is equal to the fundamental circular frequency
of the dam, ωd , namely where high amplification is expected. It may be
observed that the reservoir-dam interaction results in higher ampification for
the first natural mode of vibration, but smaller amplification for the second
natural mode. Moreover, the maximum values of amplification (for both the
first and second natural modes of vibration) occur at smaller values of the
frequency ratio ω/ωd . The harmonic load (ω) resonates with the reservoir-
dam system (RDS) at a smaller frequency than that of the fundamental
frequency of a dam, ωd , with an empty reservoir. This implies that the
fundamental frequency of vibration of the dam-reservoir system is smaller
(and hence its fundamental period is larger) than that of a dam with an
empty reservoir.
Figure 5.6 shows the value of the ratio ω/ωd for which the maximum
amplification occurs (for both of the first two natural modes), for various
values of the ratio ωr /ωd . On the same figure, the corresponding values of
ω/ωd for an empty reservoir are included for comparison. It is shown that for
both of the first two natural modes of vibration, and for all ωr /ωd ratios, the
166
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
167
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
5.3.3 Comments
For the examined case of a cantilever dam, the effects of dam-reservoir inter-
action are mainly related to changes of the fundamental period of the dam
and the amplification of its response. It is found that a dam with a full reser-
voir has a larger fundamental period, Td (i.e. the presence of the reservoir
increases the fundamental period of the dam and the latter exhibits a softer
response). This may be attributed to the added mass, m, of the reservoir
which vibrates with the dam, while it does notqprovide any additional (shear
or bending) stiffness, k to the dam (Td = 2 π m/k).
The amplification of the input acceleration can be either higher or lower
than the amplification in a dam with an empty reservoir. This depends on
the relative magnitude of the circular frequency of the load, ω, and the funda-
mental circular frequencies of the dam and reservoir, ωd and ωr respectively.
The effects of reservoir-dam interaction are more pronounced close to the
frequency ratios ω/ωd ≈ ωr /ωd ≈ 1, i.e. when there is resonance between
the load, the dam and the reservoir. Besides, high values of amplification
are observed for ω/ωd ≈ ωr /ωd , i.e. when the circular frequency of the load
is close to the fundamental circular frequency of the reservoir (ω ≈ ωr ), i.e.
when there is resonance between the load and the reservoir.
The studied case of reservoir-dam interaction effects on the dynamic re-
sponse of a cantilever dam structure could be relevant to thin slender arch
dams which have a similar cross-sectional geometry. Moreover, this study
could also provide some insight into the dynamic behaviour of canal locks,
which again have a slender geometry.
Finally, in comparing the results of this study with those of Chopra (1968)
(Figure 2.32), the following points should be noted:
• Chopra’s 3D figure shows generally only one peak value in the F −ω/ωd
curves, whereas this work (Figures 5.2 and 5.4) shows three peak values
which correspond to the first three natural modes of vibration of the
dam. This is because the work of Chopra (1968) considered a dam
deforming in its fundamental mode only.
168
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
169
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
W B L
Figure 5.9: Geometry of the triangular earth dam considered.
170
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
14
12
Amplification, |F|
10
0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5
Frequency ratio, ω/ωd
Figure 5.10: Amplification spectrum of the earth dam with an empty reser-
voir.
2π 2π πVs πVs
ωd = = = 0.765 (5.4)
Td 2.613 H H
171
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
172
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
ωr/ωd
3
2
1
0
4.5
4
3.5
d
ω/ω
2 3
2.5
0.5 1.5
1
0
8
0
18
16
14
12
10
Amplification, |F|
173
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
14
12
Amplification, |F|
10
0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5
Frequency ratio, ω/ωd
Figure 5.12: Amplification spectrum of the earth dam with a full (ωr /ωd = 1)
and an empty reservoir.
3.5
First Second Third
mode mode mode
3
Frequency ratio, ωr/ωd
2.5
1.5
Full reservoir
1
Empty reservoir
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5
Frequency ratio, ω/ωd
Figure 5.13: Amplification peaks for the earth dam with a full and an empty
reservoir.
174
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
16
14
12
Amplification, |F|
10
4 Full reservoir
Empty reservoir
2
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Frequency ratio, ωr/ωd
Figure 5.14: Amplification of acceleration, |F | for the earth dam with respect
to the ratio of the fundamental circular frequency of the reservoir over that
of the dam, ωr /ωd , for the first mode of vibration.
6
Amplification, |F|
2
Full reservoir
Empty reservoir
1
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Frequency ratio, ωr/ωd
Figure 5.15: Amplification of acceleration, |F | for the earth dam with respect
to the ratio of the fundamental circular frequency of the reservoir over that
of the dam, ωr /ωd , for the second mode of vibration.
175
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
3.5
Amplification, |F|
2.5
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Frequency ratio, ωr/ωd
Figure 5.16: Amplification of acceleration, |F | for the earth dam with respect
to the ratio of the fundamental circular frequency of the reservoir over that
of the dam, ωr /ωd , for the second mode of vibration.
for a dam with an empty reservoir (dashed line) and a dam with a full
reservoir (solid line), for which ωr /ωd = 1. In contrast to the cantilever dam
case (Figure 5.4), the interaction does not have a significant impact on the
amplification spectrum. The amplification values are very similar for the two
cases (full and empty reservoir) and their maxima occur for almost identical
values of the ω/ωd ratio. This means that in contrast to the case of the
cantilever dam, the presence of the reservoir has a minimal impact on both
the amplification and the fundamental period of the dam.
Figure 5.13 shows the value of the ω/ωd ratio for which the maximum
amplification occurs (for the first three natural modes), for various values of
the ωr /ωd ratio. In the same figure, the values of ω/ωd for an empty reservoir
are included for comparison. It is shown that generally for all three natural
modes of vibration, and for all ωr /ωd ratios, the maximum amplification
occurs very close to the value of the ω/ωd ratio of the case of the empty
reservoir, which again suggests that the presence of the reservoir does not
affect the natural periods of the dam significantly.
Finally, Figures 5.14, 5.15 and 5.16 show the maximum amplification with
176
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
respect to the ωr /ωd ratio for the first, second and third natural modes of
vibration respectively. The amplification response for an empty reservoir is
also shown on the same figure with dashed lines for comparison. Referring to
all three modes of vibration, the amplification for the case of a full reservoir
seems to be very close to the case of an empty reservoir. However, some
small differences can still be observed which show that the case of the full
reservoir generally results in slightly smaller amplification values than those
of the empty reservoir.
5.4.3 Comments
For the examined case of a trapezoidal earth dam, the effects of dam-reservoir
interaction were found to be insignificant. In contrast to the previous case
of a cantilever dam, changes of both the fundamental period of the dam and
the amplification of its response were still observed but they were very small.
This is believed to be due to (a) the sloped upstream face and (b) the large
volume of a trapezoidal earth dam. As far as the former reason is concerned,
it is long known (Zangar, 1953) that the hydrodynamic pressures induced
on a sloped upstream face are smaller than those on a vertical face (see also
Section 2.5). Moreover, regarding the second reason, the inertial effects from
the additional mass from the reservoir are small compared to the inertia of a
large earth dam. This is due to the large volume (and hence mass) of earth
dams and therefore the influence of the reservoir on the dynamic response of
the dam seems to be minimal.
A further study including trapezoidal dams with more gentle slopes could
examine whether the above-mentioned conclusions apply to other dam slopes.
That could perhaps identify for which values of the slope angle the interaction
effects become significant.
177
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
178
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.3
0.2
Acceleration, a [g]
0.1
−0.1
−0.2
−0.3
0 10 20 30 40 50 60
Time, t [sec]
0.9
0.8
Spectral Acceleration, Sa [g]
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
179
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Table 5.1: Cases of combinations of ω/ωd and ωr /ωd frequency ratios exam-
ined under seismic loading.
2.5
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Frequency ratio, ω/ωd
Figure 5.19: Cases of combinations of ω/ωd and ωr /ωd frequency ratios ex-
amined under seismic loading.
180
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.5
Acceleration, a [g]
−0.5
−1
−1.5
−2
0 10 20 30 40 50 60
Time, t [sec]
0.6
0.4
Acceleration, a [g]
0.2
−0.2
−0.4
−0.6
−0.8
0 10 20 30 40 50 60
Time, t [sec]
181
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1.5
1
Acceleration, a [g]
0.5
−0.5
−1
−1.5
−2
−2.5
0 10 20 30 40 50 60
Time, t [sec]
0.5
Acceleration, a [g]
−0.5
−1
−1.5
−2
−2.5
0 10 20 30 40 50 60
Time, t [sec]
182
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.5
Acceleration, a [g]
−0.5
−1
−1.5
0 10 20 30 40 50 60
Time, t [sec]
7
Spectral acceleration, Sa [g]
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
183
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
4
Spectral acceleration, Sa [g]
3.5
2.5
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
10
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
184
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
10
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
185
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Table 5.2: Material and geometric properties of the cantilever dam cases
considered under seismic loading.
for the dam with a full reservoir and for all 5 cases (A, B, C, D and E), the
peak value of amplification occurs at a value of the ω/ωd ratio smaller than
the corresponding for the dam with an empty reservoir (see Figures 5.25 -
5.29). This suggests that the fundamental period of the dam is larger when
reservoir-dam interaction is considered. The latter conclusion is in agreement
with the earlier observation in Section 5.3 for harmonic load, i.e. the load
resonates with the dam-reservoir system at a smaller value of natural circular
frequency, therefore the fundamental period of the dam-reservoir system (i.e.
dam with a full reservoir) is larger than that of a dam with an empty reservoir.
Moreover, the value of amplification of accelerations was larger for cases
A-D with a full reservoir than those with an empty reservoir (Figures 5.25
- 5.28). This is in agreement to the previous investigation using harmonic
loading (see Figure 5.4), where higher values of amplification were observed
for the frequency combinations of cases A-D for the dam with a full reservoir.
It can be observed from Figure 5.4, that the amplification with a full reservoir
was higher than with an empty reservoir for the following combinations of
(ω/ωd , ωr /ωd ) = (1, 1), (3, 3), (1, 3), (1/3, 1), which correspond to cases
A-D.
Besides, for case E (i.e. for ω/ωd = ωr /ωd = 1/3), the amplification of
accelerations was found to be smaller for a dam with a full reservoir than that
for a dam with an empty reservoir (Figure 5.29). This is again in agreement
with the previous investigation regarding harmonic loading (see Figure 5.4),
where the amplification for ω/ωd = ωr /ωd = 1/3 was higher in a dam with a
full reservoir than a dam with an empty reservoir.
186
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Table 5.3: Material and geometric properties of the earth dam cases consid-
ered under seismic loading.
In all the cases, the material properties of the earth dam were: Poisson’s
ratio, ν = 0.3, mass density, ρ = 2000 kg/m3 and Rayleigh damping, ξ = 5%.
The height, H and shear wave velocity, Vs of the dam were changed in each
case as listed in Table 5.3, to achieve the combinations of frequencies listed
in Table 5.1, but in all cases keeping the value of the p-wave velocity of the
reservoir constant, Vp = 1483 m/s. The reservoir water was again modelled
as a linear material with the same values for the bulk and shear modulus as
before (see Section 5.3).
Figures 5.30 - 5.34 show the acceleration time-histories at the crest of the
dam (for both full and empty reservoirs), whereas Figures 5.35 - 5.39 show
the corresponding 5% damped response spectra. It may be observed from
all figures that the observations for the cantilever dam case do not apply to
the earth dam. Although some minor differences between the full and empty
reservoir cases may be noticed in the acceleration time-histories (Figures 5.30
- 5.34), no major differences exist in the response spectra (Figures 5.35 - 5.39).
This means that both the frequency content and the value of amplification
were not significantly affected by the presence of the reservoir, for all cases
A-E, i.e. for various combinations of the frequency ratios. The previous
187
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1
Acceleration, a [g]
0.5
−0.5
−1
−1.5
0 10 20 30 40 50 60
Time, t [sec]
Figure 5.30: Acceleration time-history at the crest of the earth dam - Case
A.
0.2
Acceleration, a [g]
−0.2
−0.4
−0.6
−0.8
0 10 20 30 40 50 60
Time, t [sec]
Figure 5.31: Acceleration time-history at the crest of the earth dam - Case
B.
188
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1
Acceleration, a [g]
0.5
−0.5
−1
−1.5
0 10 20 30 40 50 60
Time, t [sec]
Figure 5.32: Acceleration time-history at the crest of the earth dam - Case
C.
0.5
Acceleration, a [g]
−0.5
−1
−1.5
0 10 20 30 40 50 60
Time, t [sec]
Figure 5.33: Acceleration time-history at the crest of the earth dam - Case
D.
189
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.2
Acceleration, a [g]
−0.2
−0.4
−0.6
−0.8
0 10 20 30 40 50 60
Time, t [sec]
Figure 5.34: Acceleration time-history at the crest of the earth dam - Case
E.
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
190
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
2.5
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
191
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
192
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
5.5.3 Comments
This section considered the effects of RDI on cantilever and earth dams under
seismic loading. A real acceleration time-history was chosen and applied on
both types of dams in a series of dynamic analyses in the time domain. The
cases were chosen in terms of the frequency ratios ω/ωd and ωr /ωd so that
a direct comparison can be made with the results of the earlier harmonic
studies. The value of ω, was related to TEQ which was called the dominant
period of the load, and it was determined as the circular frequency for which
the highest value of spectral acceleration, Sa occurs (see Figure 5.18).
However, as it may be observed from Figure 5.18, there are two peaks in
the spectral acceleration values. The second peak, at T = 0.52s is the one
chosen as dominant period, and the other one, at T = 0.27s corresponds to
a slightly smaller value of Sa . It is expected to observe resonance at both
periods (0.27s and 0.52s) if they are close to the fundamental frequencies of
the dam and the reservoir. The chosen cases A-E (Table 5.1) were based on
the peak at T = 0.52s and indeed the results for the seismic loading (Section
5.5) are in agreement to the results for the harmonic loading (Sections 5.3 and
5.4). The second peak value of Sa did not affect the results of this section,
as no resonance was expected for the cases chosen.
5.6 Conclusions
This chapter describes the investigation carried out to assess the effects of an
upstream reservoir on the dynamic behaviour of dams. The main aim was to
examine whether the presence of the reservoir alters the response of a dam
and to identify the cases in which the dynamic reservoir-dam interaction may
have detrimental effects. Both rectangular cantilever and trapezoidal homo-
geneous earth dams were considered under harmonic and seismic loading.
193
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
194
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
effects from the additional mass from the reservoir are small compared
to the inertia of a large earth dam. However, the amplification in an
earth dam with a full reservoir seems to be slightly smaller than the
corresponding amplification for a dam with an empty reservoir. This
could suggest that the presence of the reservoir may damp the dynamic
response of earth dams (through the interactive oscillation of two bodies
with different vibrational characteristics, such as fundamental period).
• The calculated values for the fundamental period for both dams (can-
tilever and trapezoidal earth) obtained from this study were found
to be slightly different than the values obtained by using analytical
relations from the literature (Euler-Bernoulli beam theory and shear
beam). This is in agreement with earlier observations of previous re-
searchers for both types of dams. The differences were attributed to the
simplifying assumptions of the analytical methods, as the FE method
(employed in this study) is able to overcome their limitations (as it
considers both shear and bending deformations).
195
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Chapter 6
NUMERICAL ANALYSIS OF
LA VILLITA DAM
6.1 Introduction
This chapter describes numerical analysis of the seismic response of a real
earth dam. The finite element method was employed to perform static and
time-domain dynamic nonlinear elasto-plastic analyses of La Villita earth
dam in Mexico. La Villita dam is considered to be a reasonably well-
documented case because of available recorded data, material properties and
because of an existing previous research that investigated its response.
The objective of the current work is to reproduce the observed behaviour
of the dam and compare the numerical predictions with the recorded re-
sponse. This was performed in order to validate the numerical approach
adopted for modelling earthfill dams in seismic conditions. In achieving this,
several issues related to numerical analysis were addressed which provided
an insight into the real behaviour of the dam and also provided valuable
feedback and information about considerations that need to be taken when
performing seismic analysis of earth dams using the finite element method.
196
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 6.1: Aerial photo of La Villita dam (Google Earth, 2010) with appur-
tenant facilities: (1) Reservoir, (2) Embankment, (3) Hydroelectric power
plant, (4) Spillway, (5) Downstream Balsas river.
The dam’s history states that between 1975 and 1985 it experienced 6
major earthquakes. Although no failure was observed, displacements of some
magnitude were recorded (Reséndiz et al., 1982; Gonzalez-Valencia, 1987;
Elgamal, 1992). It is considered to be a well documented case of seismic dam
response, because measured displacements and acceleration time histories
are available for several locations on the dam together with properties of the
materials that constitute dam’s body and foundations. Figure 6.1 shows an
aerial photo of the dam with the additional facilities.
197
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
approximately one year, during 1967. Water was impounded in the reservoir
after its construction, i.e. in early 1969 and it operated normally for 7 years
until the first major earthquake event in October 1975.
Before the construction of the dam, the river stream was diverted so that
the construction site was kept dry. Then, a series of holes were drilled and
grouted under the proposed axis of the embankment in order to create the
grout curtain whose function is to control seepage under the dam. After that,
the dam body was constructed in a number of successive thin soil layers which
were compacted upon placement.
Impoundment of the water in the reservoir was carried out over about
half a year, until the reservoir reached its highest level, at approximately
54m. Following the impoundment, the dam was in constant operation and
no rapid drawdowns were reported.
The first significant seismic event occurred in 1975, and then a series of
strong earthquakes took place in the area close to the dam in the following
ten years. Table 6.1 summarizes the details of 6 significant earthquake events
from information obtained from Sociedad Mexicana de Ingenieria Sismica
(2000). Although not all the acceleration time histories could be obtained,
the values in the table were completed using information from Elgamal et al.
(1990) and Elgamal (1992).
As it may be seen from Table 6.1, the earthquakes between 1975 and 1985
198
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
resulted in high values of PGA at the crest. This implies that they could
have caused serious damage to the dam body.
However, the dam performed well without any major damage and the con-
sequences of these events were only some permanent displacements, which left
the dam still operational. No evidence of liquefaction was observed, neither
within the embankment or nearby. Figure 6.2 shows recorded displacements
with time at several points on the dam body. Further information about
its response during the 1985 earthquakes can be found in Gonzalez-Valencia
(1987).
Figure 6.2: Recorded displacements on some points on the dam body (Elga-
mal, 1992).
199
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
through the alluvium below the dam. Figure 6.3 shows a typical transversal
(upstream-downstream) cross-section, Figure 6.4 shows a longitudinal cross-
section (Elgamal, 1992) outlining the geological profile of the site, whereas
Figure 6.5 (Elgamal et al., 1990) shows a plan view of the dam.
From the materials that make the dam, the clay and the alluvium seem to
be the less permeable and expected to behave in an undrained manner. The
rest of the materials (sand filters and rockfill shells) are reasonably coarse
and are expected to behave in a drained manner. Moreover, in-situ explosion
tests showed that the alluvium stratum underneath the dam is not considered
to be prone to liquefaction, as reported by Elgamal (1992).
A number of tests was carried out in order to characterise the materials
undertaken by the Comision Federal de Electricidad (CFE) of Mexico (Comi-
200
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
sion Federal de Electricidad, 1976, 1979, 1980, 1987) and the Universidad
Nacional Autónoma de México (UNAM) (Comision Federal de Electricidad
and Univercidad Nacional Autónoma de México, 1976). Table 6.2 provides
a summary of obtained material properties of La Villita dam, although no
detailed information was found about the types of tests carried out. Elgamal
(1992), who performed shear beam analysis of La Villita dam (see Section
6.3), provided information about the spatial variation of the shear modulus,
G of the dam and foundation according to Table 6.3.
Table 6.3: Spatial variation of the maximum shear modulus, Gmax of the
dam and foundation material, as used by Elgamal (1992).
201
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
at the right bank. The two accelerograms on the embankment are on the
downstream slope of the dam, one close to the crest and one at the down-
stream berm. Figure 6.5 (Elgamal et al., 1990) shows the location of these
accelerometers on a plan view of the dam.
202
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.05 EQ1
−0.05
0 1 2 3 4 5 6
0.05
EQ2
−0.05
0 1 2 3 4 5 6
0.02
Acceleration, a, [g]
EQ3
0.01
−0.01
−0.02
0 1 2 3 4 5 6
0.2
EQ5
0.1
−0.1
−0.2
0 10 20 30 40 50 60 70
0.04
EQ6
0.02
−0.02
−0.04
0 2 4 6 8 10 12
Time, t [sec]
203
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.5
EQ1
−0.5
0 5 10 15
0.5
EQ2
−0.5
0 2 4 6 8 10 12 14 16 18
0.5
EQ3
Acceleration, a, [g]
−0.5
0 5 10 15 20 25 30 35
0.5
EQ4
−0.5
10 5 10 15 20 25 30
EQ5
0.5
−0.5
0 10 20 30 40 50 60 70 80
0.5
EQ6
−0.5
0 10 20 30 40 50 60
Time, t [sec]
204
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1
EQ1 − Bedrock
EQ1 − Base
0.8
EQ1 − Crest
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3 3.5 4
0.7
EQ2 − Bedrock
Spectral Acceleration, Sa, [g]
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3 3.5 4
1.4
EQ3 − Bedrock
1.2 EQ3 − Base
EQ3 − Crest
1
0.8
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
205
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1.4
EQ4 − Base
1.2 EQ4 − Crest
0.8
0.6
0.4
0.2
0
20 0.5 1 1.5 2 2.5 3 3.5 4
EQ5 − Bedrock
Spectral Acceleration, Sa, [g]
EQ5 − Crest
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
0.7
EQ6 − Bedrock
0.6 EQ6 − Crest
0.5
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
206
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.4
0.2
−0.2
−0.4
0 10 20 30 40 50 60 70 80
Time, t [sec]
Figure 6.10: The recorded UD motion on the crest of La Villita dam, during
EQ5 that shows an asymmetry in the peak values.
207
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 6.11: The experimental setup used by Elgamal et al. (1990) in the pre-
diction of permanent displacements of La Villita dam: Sliding block mounted
on inclined dynamic shaker.
In their research, the crest acceleration record for the 19th of September
1985 earthquake (EQ5) was compared to the acceleration response of a block
sliding on an inclined plane. The response of the sliding block was exam-
ined using the experimental setup shown in Figure 6.11, to investigate the
observed acceleration response asymmetry. Information about the material
properties was obtained from reports of the CFE (Comision Federal de Elec-
tricidad, 1976, 1979, 1980), whereas the acceleration records were obtained
from UNAM.
A series of sliding block analyses were performed by shaking the appa-
ratus and changing the sliding slope and the friction angle until the total
displacement matched the observed one, which was about 20cm for EQ5.
The final computed displacement time history is shown in Figure 6.12, for
a sliding plane of θ = 22o . Their aim was to investigate whether a sliding
displacement of about 20cm would be associated with an asymmetric ac-
celeration response of the sliding block. And indeed, the accelerations of
208
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
the sliding block, shown in Figure 6.13 were asymmetric, i.e. the positive
(downwards sliding direction) peak accelerations had higher values than the
negative peak accelerations.
From the sliding block analyses, it was found that the strong motion
records not only indicated localised deformations but also revealed the mag-
nitude of the associated yield acceleration, namely, the value of acceleration
at which sliding occurred. The asymmetry of the positive-negative accelera-
tion peaks was attributed to a sliding block type of failure.
Following this, a three-dimensional numerical model was developed by
Elgamal (1992) based on the shear beam approach (see Section 3.4). The
numerical model was utilised in the analysis of La Villita dam, during the
15th November 1975 (EQ2) and 19th September 1985 (EQ5) earthquakes.
Initial stresses within the dam (prior to the EQs) were obtained by applying
static self weight in the whole numerical model, including the dam, thus
not modelling layered construction. The SB analysis employed a nonlinear
plasticity model (Prevost, 1977) which took account of degradation of shear
stiffness. As before, the material properties were obtained from reports of
the Comision Federal de Electricidad (1976, 1979, 1980) (CFE), whereas the
acceleration records were obtained from UNAM.
Accelerations and displacements were computed and compared to those
209
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
observed (not shown here for brevity). The calculated peak accelerations
were very similar to the recorded ones for the longitudinal and vertical direc-
tions, whereas this was not the case for the upstream-downstream direction
where the calculated values were smaller than the recorded. Furthermore,
as shown in Figure 6.14, their model predicted the dynamic characteristics
(fundamental period) reasonably well and there was a fair agreement between
the calculated and recorded response spectra. However, it was noted that the
recorded acceleration spectrum was broader in frequency content and lower
in amplitude when compared to the computed values.
The difference between the observed and predicted displacements was
attributed to the occurrence of seismically induced deformations due to lo-
calised yielding. Therefore, the calculated accelerations from the SB analy-
sis were used as input in a sliding block analysis, in order to estimate the
permanent displacements of the dam. However, the computed permanent
deformations from the sliding block analysis for EQ5 (≈ 5 cm) (Figure 6.15)
were found much smaller than the observed deformations (≈ 20 cm).
The upstream-downstream earthquake response of La Villita dam during
EQ2 and EQ5 was analysed by Succarieh et al. (1993). Firstly, a 1D nonlin-
ear elasto-plastic hysteretic inhomogeneous SB model was used to represent
210
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
211
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
212
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
213
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 6.16: The finite element mesh used by Gazetas and Uddin (1994) in
the analysis of La Villita dam with a pre-defined slip surface.
(a) (b)
Figure 6.17: Computed accelerations of two points on the crest of the dam:
(a) outside and (b) inside the failing mass, by Gazetas and Uddin (1994).
214
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
shear beam approach used by Dakoulas (1993b) (see also Section 2.4.2) for
semi-cylindrical canyons in the sense that the surrounding canyon is included
in the computational model too. Figure 6.19 (Papalou and Bielak, 2001)
presents the components of the modelling system.
The method was applied on the case of La Villita dam, but considering
its behavior as being purely elastic. The main aim of that study was to
investigate dam-canyon interaction effects. A geometric model of the dam,
the canyon and the surrounding medium was considered and the excitation
consisted of SH waves at arbitrary angles of incidence. A cross-section of
the model used in the analysis is shown in Figure 6.20 (Papalou and Bielak,
2001) and the different layers correspond to different material properties
(shear modulus, G), similar to the model of Elgamal (1992) (see also Section
6.2.3 and Table 6.3). A parametric study was undertaken to examine the
effects of a wide range of the shear wave velocity of the surrounding bedrock
on the acceleration response of the dam.
It was found that dam-canyon interaction affects greatly the seismic re-
sponse of the dam as the dynamic coupling of the dam and the surround-
ing deformable canyon resulted in smaller strains than the case with a rigid
canyon due to the energy radiated back to the canyon. Finally, a rigid canyon
yielded higher values of accelerations and therefore the authors suggested
that including dam-canyon interaction effects in seismic design of dams would
prevent unnecessary conservatism.
The procedure developed by Papalou and Bielak (2001) was later ex-
tended (Papalou and Bielak, 2004) to account for nonlinear material be-
haviour and was applied again to the case of La Villita dam. In that investi-
gation, material nonlinearity was modelled using multi-yield surface plasticity
theory (Prevost, 1977), and the seismic response in the upstream-downstream
direction was investigated. Similar to the previous paper, the material prop-
erties were taken from Elgamal (1992) and the shear wave velocity of the
rock canyon was varied in a parametric study.
Then, it was suggested that even with the inclusion of material nonlin-
earity, the effects of dam-canyon interaction were still significant. In that
case however, they were more pronounced within the interior of the dam (i.e.
215
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 6.18: Comparison of calculated (thin line) and final recorded (thick
line) permanent crest displacements of La Villita dam by Gazetas and Uddin
(1994).
Figure 6.19: The system analysed by Papalou and Bielak (2001, 2004): (a)
3D dam-canyon model, (b) cross-section of dam in the UD direction, (c)
cross-section of dam and canyon in the longitudinal direction and (d) free
body diagram of a shear wedge element.
216
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
away from the dam-canyon boundary) where shear strains were larger. Con-
sidering the results in the frequency domain, it was found that at mid-crest,
nonlinearity increased the response at high frequencies, whereas in the inte-
rior of the dam the response was significantly reduced at low and mid-range
frequencies.
The calculated peak accelerations in both of the analyses of Papalou and
Bielak (2001, 2004) (linear elastic and nonlinear elasto-plastic respectively)
were significantly smaller than the recorded values. Figure 6.21 shows a
comparison of the calculated and recorded accelerations at the midcrest.
Clearly, the elastic analyses yielded higher values of peak accelerations, which
however were still much smaller than the recorded. As an example, the
recorded peak acceleration was 1.7 times higher than that obtained from the
linear analysis for the shear wave velocity of the rock, Vsr = 4000 m/s (see
Figure 6.21).
Following the earlier explanation of Elgamal (1992), Papalou and Bielak
(2004) commented that the spikes in the crest record of La Villita dam were
due to a stick-slip displacement mechanism in the zone under the accelerom-
eter. Those spikes appeared at the end of each slip phase, as a consequence
of a sudden change in the inertial load. In addition, it was also reported
that those spikes could not be represented by the coarse mesh used in the
analyses because of their very high frequencies.
6.3.2 Comments
The case of La Villita dam is of significant research interest because of the
available data: material properties, geometry, earthquake records, and mea-
sured displacements. These reasons make it a well-documented case which
attracted the attention of various researchers. However, considering the pre-
viously reviewed information from the literature, it seems that there is still
a number of issues that need to be addressed.
Section 3.7 discussed several issues that need to be taken into account
in the dynamic analysis of earth dams. But, from the work found in the
literature, it may be noticed that those issues were not fully taken into con-
217
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 6.20: The model used by Papalou and Bielak (2001, 2004) in the anal-
ysis of La Villita dam. The different zones correspond to different material
properties (Gmax ).
218
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
sideration.
Firstly, the entire history of the dam (construction, operation and sub-
sequent earthquake events) was not simulated in order to obtain the corre-
sponding stress states just before the earthquake. This could provide a more
realistic insight into the stress states which is essential in order to appro-
priately consider soil plasticity (and hence stability and deformations of the
dam structure). Secondly, coupling of the soil skeleton and pore water (cou-
pled consolidation analyses) were not performed, which are useful in order to
capture the steady-state seepage through the dam and realistically consider
the effect of the pore water and of the reservoir.
Moreover, with the exemption of Papalou and Bielak (2004) who provided
a minimum of nine nodes per wavelength, issues related to the FE mesh dis-
cretisation and size of the elements were not strictly taken into consideration.
For example, a rather crude mesh (Figure 6.16) was used in the numerical
analyses of Gazetas and Uddin (1994). A fine mesh could allow the prop-
agation of high frequency wave components and therefore would be able to
predict better any possible high frequency peaks. Besides, the reservoir water
was not modelled and therefore hydrodynamic forces from the reservoir were
not simulated, not even as an additional external static force. Inclusion of
the reservoir could provide an insight into reservoir-dam interaction effects
and into any potential amplifications of the seismic response of dams.
Concerning the accelerations at the crest of La Villita dam, the reason
behind the record’s peculiar asymmetry was not fully understood. Most
researchers tend to agree that the spikes in the record were due to a localised
slip failure on the dam slope. This is strongly supported by the work of
Gazetas and Uddin (1994) who compared the response inside and outside a
potential failing mass. Moreover, the inability to predict such high-frequency
spikes was attributed to the geometry of the FE mesh used in the analyses
(Papalou and Bielak, 2004).
As far as the permanent displacements of the crest are concerned, earlier
work (Elgamal et al., 1990; Elgamal, 1992) suggested that a sliding block
analysis was required in order to capture their magnitudes. However, later
work (Gazetas and Uddin, 1994) provided better estimates (the calculated
219
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
displacements were 1.45 times higher than the recorded for EQ5) of the per-
manent displacements but it still requires further improvement. An elasto-
plastic analysis, taking account of soil plastic yielding could perhaps provide
a better insight into the actual seismic behaviour of La Villita dam (i.e.
without the need of a pre-defined failure surface).
220
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
60m
A B
70m
462m
Figure 6.22: 2D Finite Element mesh used for the analysis of La Villita dam
The value of the maximum shear stiffness, Gmax used by Elgamal (1992) is
listed in Table 6.3. In this work, the value of Gmax was taken to vary linearly
with depth in the dam embankment, according to Equation 6.1 (where, z is
the elevation, i.e. 0 at the base of the dam and 60m at the crest of the dam)
and shown in Figure 6.23, i.e. varying from 140000 kPa at the crest of the
221
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
dam to 260000 kPa to the base of the dam. The corresponding value of Gmax
for the foundation alluvium was assumed constant and equal to 200000 kPa.
40
Elevation, z [m]
30
20
10
0
0 50000 100000 150000 200000 250000 300000
Maximum Shear Modulus, Gmax [kPa]
Figure 6.23: Spatial variation (in the vertical direction) of shear stiffness,
Gmax . The grey solid line corresponds to the values used by Elgamal (1992)
(see also Table 6.3) and the black solid line corresponds to the values used
in this work.
The size of the elements was discussed in Section 3.7 and according to
Equation 3.21, it depends on the maximum significant frequency of the load,
fmax (≈ 10 Hz, see Figure 6.24) and the smallest considered shear wave
velocity, VSmin (≈ 265 m/s). This means that, ∆l = 1/5 ∼ 1/4 · 265/10 =
5.3 ∼ 6.625 m, and therefore, the maximum size of the element used was 6
m.
The bottom boundary of the mesh was placed at the interface of the
bedrock and the foundation alluvium, as the bedrock was considered to be
significantly stiffer than the foundation alluvium (Papalou and Bielak, 2001,
2004), and therefore it does not need to be modelled (see also Section 3.7).
222
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.012 0.14
0.12
0.01
Fourier Amplitude, FA [g]
0.004
0.04
0.002 0.02
0 0
−2 −1 0 1 2 −2 −1 0 1 2
10 10 10 10 10 10 10 10 10 10
Frequency, f [Hz] Frequency, f [Hz]
The lateral boundaries were placed at such distance from the toe of the
dam (distance A-B in Figure 6.22) so that the response close to the bound-
aries is similar to the free-field response (see Section 3.7). For this reason, 1D
column (site response) and 2D plane strain analyses (see Section 6.6) were
performed and compared to determine the appropriate distance (from the toe
of the dam) at which the lateral boundaries should be placed so that there
is no interaction with the structure studied. The required distance from the
toe of the dam was found to be 60 m.
223
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
224
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
For the calibration of the CNL model, the materials of the dam were
divided in three categories: (a) clay core, (b) sand filters and (c) rockfill and
alluvium. The clay core was considered to be a plastic material and the CNL
model was calibrated against the curves of Vucetic and Dobry (1991) for
Ip = 30%. The sand filters were considered to be of low Plasticity Index Ip ,
and therefore they were calibrated against empirical curves for such material,
which in that case, were the curves of Seed et al. (1986). Finally, the rockfill
shells and the alluvial foundation were calibrated against the curves of Rollins
et al. (1998) for rockfill materials.
Table 6.6 summarizes the results of the calibration process for the Log-
arithmic CNL and for all materials. The results, both in terms of shear
modulus (G/Gmax − γ) and damping (ξ − γ) variation are shown graphically
in Figures 6.25 - 6.30 for all the dam materials. It may be observed that
a fair match was achieved between the curves predicted by the CNL model
(solid black line, for the middle of the dam, i.e. for Gmax = 200000 kPa)
and the empirical ones from the literaure (dashed black line), for both shear
stiffness degradation and damping.
As reported earlier (Equation 6.1 and Figure 6.23), the value of the max-
imum shear stiffness, Gmax in the embankment varied linearly with depth
(140000-260000 kPa). However, one set of calibration parameters (EdL , c,
JL and Gmin ) was specified for the whole embankment domain, which corre-
sponds to a calibration for Gmax = 200000 kPa, which is the average value.
Therefore, the materials at the top and bottom of the embankment which
have smaller (140000 ∼ 200000 kPa) and larger (200000 ∼ 260000 kPa) shear
stiffness respectively, exhibit a different variation of G/Gmax −γ and ξ −γ for
225
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
the specified model parameters (EdL , c, JL and Gmin , listed in Table 6.6), as
presented in Figures 6.25 - 6.30. The solid grey line corresponds to the top of
the embankment (smallest value of Gmax = 140000 kPa) and the dashed grey
line corresponds to the bottom of the embankment (largest value of Gmax =
260000 kPa).
It should be noted that the calibration of the Logarithimc model allowed
either (a) an excellent match between the model and the empirical curves for
the variation of shear stiffness (G/Gmax − γ) but at the same time, a poor
match for the damping curves (ξ − γ) or (b) the opposite, i.e. an excellent
calibration for the damping curves (ξ − γ) and a poor calibration for the
shear stiffness curves (G/Gmax − γ). In the first case (a), i.e. if the shear
stiffness curves (G/Gmax − γ) were matched, the CNL model would predict
higher damping than the empirical curves for large values of the shear strain,
γ. That would result in an overdamped response of the dam structure. In
the second case (b), i.e. if the damping curves (ξ − γ) were matched, the
CNL would predict lower values of the shear stiffness than the empirical
curves. That would then result in a softer response of the dam, with a
longer fundamental period and perhaps smaller amplification in the higher
frequencies range. It was therefore decided to calibrate the CNL on both
G/Gmax − γ and ξ − γ curves, obtaining an optimum match for both sets of
curves.
Figure 6.31 shows a comparison between the shear stiffness degradation
(G/Gmax − γ) predicted by the Logarithmic CNL model and that used by
Elgamal (1992). It may be observed that there is a fair agreement between
the two cases. The smaller values of shear stiffness could be attributed to
the calibration against both stiffness and damping, as discussed above. The
larger difference for the case of the sand filters should not significantly affect
the dynamic characteristics of the dam, as the sand filters are just a very
small part of the whole embankment.
Figures 6.25, 6.27 and 6.29 show that the normalised shear stiffness degra-
dation, G/Gmax − γ, varies with the elevation within the dam (i.e. points
at top, middle and bottom). It is shown that the values of G at the bottom
of the dam (which are the largest values, due to the adopted spatially inho-
226
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.9
0.8
Shear Modulus, G/Gmax []
0.7
0.6
0.5
0.4
0.3
Figure 6.25: Calibration of the Logarithmic CNL model against the empirical
curves of Vucetic and Dobry (1991) for the clay core - Stiffness degradation
(G/Gmax − γ)
30
Model − Top
25 Model − Middle
Model − Bottom
Vucetic & Dobry, 1991
20
Damping, ξ [%]
15
10
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
Figure 6.26: Calibration of the Logarithmic CNL model against the empirical
curves of Vucetic and Dobry (1991) for the clay core - Damping (ξ − γ)
227
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.9
0.8
Shear Modulus, G/Gmax []
0.7
0.6
0.5
0.4
0.3
Model − Top
0.2
Model − Middle
0.1 Model − Bottom
Seed et al., 1986
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
Figure 6.27: Calibration of the Logarithmic CNL model against the empir-
ical curves of Seed et al. (1986) for the sand filters - Stiffness degradation
(G/Gmax − γ)
30
Model − Top
25 Model − Middle
Model − Bottom
Seed et al., 1986
20
Damping, ξ [%]
15
10
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
Figure 6.28: Calibration of the Logarithmic CNL model against the empirical
curves of Seed et al. (1986) for the sand filters - Damping (ξ − γ)
228
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.9
0.8
Shear Modulus, G/Gmax []
0.7
0.6
0.5
0.4
0.3
Model − Top
0.2
Model − Middle
0.1 Model − Bottom
Rollins et al., 1998
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
Figure 6.29: Calibration of the Logarithmic CNL model against the empirical
curves of Rollins et al. (1998) for the rockfill and foundation alluvium -
Stiffness degradation (G/Gmax − γ)
30
Model − Top
25 Model − Middle
Model − Bottom
Rollins et al., 1998
20
Damping, ξ [%]
15
10
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
Figure 6.30: Calibration of the Logarithmic CNL model against the empirical
curves of Rollins et al. (1998) for the rockfill and foundation alluvium -
Damping (ξ − γ)
229
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.9
0.8
Shear Modulus, G/Gmax []
0.7
0.6
0.5
0.4
0.3
Elgamal (1992)
0.2
Model (Vucetic & Dobry, 1991) − Middle
0.1 Model (Seed et al., 1986) − Middle
Model (Rollins et al., 1998) − Middle
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
300000
Model − Top
Model − Middle
Model − Bottom
250000
Shear Modulus, G [kPa]
200000
150000
100000
50000
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
230
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
300000
Model − Top
Model − Middle
Model − Bottom
250000
Shear Modulus, G [kPa]
200000
150000
100000
50000
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
300000
Model − Top
Model − Middle
Model − Bottom
250000
Shear Modulus, G [kPa]
200000
150000
100000
50000
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
231
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
mogeneous stiffness profile) degrade more with the induced shear strain, γ.
However, as shown in Figures 6.32 - 6.34, which present the absolute values
of shear modulus, G, the variation of the shear modulus within the dam,
occurs for the strains, up to γ = 0.04 ∼ 0.2 %, where it reaches the assigned
minimum stiffness value (Gmin ). Therefore, the stiffness inhomogeneity holds
for small values of the induced shear strain. Such small strains were reported
in this work (see Section 6.6).
6.5.1 Construction
After the establishment of initial stress conditions, i.e. level ground (with the
water level being at 2 m depth), the embankment was constructed. Firstly, a
concrete block which serves as a foundation of the clay core was constructed
along with the placement of the grout curtain underneath. Then the em-
bankment was built in 10 successive layers of 6 m thickness. The clay core
was considered to have suction due to its compacted nature, and a value
of 50 kPa was specified, i.e. pf = −50kP a (positive values correspond to
compressive stresses, whereas negative values correspond to tensile stresses).
This stage lasted for a year (1967) as data about the dam’s history suggested
(Table 6.4). Each layer was constructed over a number of increments in order
232
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
233
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
tion). During that stage, all the boundary conditions used before remained
the same, and only time elapsed.
Following the one-year consolidation phase, water was impounded in the
reservoir. That was modelled over a total duration of 6 months. Water level
was raised in a single layer which required 10 increments. Therefore, the wa-
ter in the reservoir was modelled as an additional external boundary stress on
the upstream face of the dam up to a height of 54 m (the higher level of the
reservoir), thus allowing 6 m for freeboard. This modelled the hydrostatic
pressure from the reservoir. Besides, an additional boundary stress was ap-
plied on the upstream riverbed alluvium equal to the maximum hydrostatic
value. At the same time of the application of the external boundary stress,
the pore pressure in the elements of the upstream rockfill and sand filters was
prescribed to be in equilibrium with the externally applied boundary stress
(i.e. hydrostatic). The approach followed to model water impoundment is
similar to that used by Kovacevic (1994) for the static analysis of Roadford
dam in the UK. A similar procedure was followed by Sica et al. (2008) for
the static and dynamic analysis of El Infiernillo dam in Mexico. They noted
that the water filling was simulated by applying a time dependent hydro-
static distribution of pore water pressure following the height of the reservoir
water and a normal stress distribution equal to the pore water pressures was
added to the wet (upstream) boundary of the core.
The deformation BCs were the same as before (i.e. during the consolida-
tion stage). The hydraulic BCs on the boundary of the mesh were also the
same except for the upstream lateral boundary of the mesh on which the pre-
scribed pore water pressures increased according to the pore pressure change
due to the water level rise. Moreover, the hydraulic BC on the upstream face
of the core was no longer precipitation. It had prescribed values of the pore
water pressure according to the elemental pore pressures prescribed within
in the upstream rockfill (i.e. a hydrostatic variation). This allowed a water
seepage through the core according to the value of the permeability, Kclay
(10−10 m/s) of the clay.
Figure 6.35 shows the pore water pressure distribution after the end of
the water impoundment, whereas Figure 6.36 shows the flownet in the core.
234
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 6.35: Contours of pore water pressure in the dam after the end of the
reservoir impoundment (positive values correspond to compressive stresses).
235
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Figure 6.36: Flow net (contours of stream and potential functions) in the
clay core after the end of the reservoir impoundment.
It may be observed from the first figure that there are pore pressures in the
upstream rockfill as a result of the water impoundment. This distribution is
linear in the vertical direction as the water is expected to penetrate quickly in
the coarse rockfill and hydrostatic conditions are established. Moreover, the
contours of pore pressure drop in the clay core which means that the pressure
reduces in the downstream side of the core, compared to the hydrostatic
values on the upstream side. There is still some suction in the upper part
of the clay core, which indicates that a part of the core is still not fully
saturated. Besides, comparing the two parts of the foundation alluvium, the
upstream part had significantly higher values of water pressure which are a
result of the reservoir water. Finally, as far as the flow net in the core is
concerned, flow lines and equipotential lines were clearly formed indicating
seepage from the upstream to the downstream side of the clay core.
236
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
6.5.3 Consolidation
After the reservoir impoundment was finished, i.e. in mid 1969, the dam
started to operate normally until late 1975, when two major seismic events
occurred, EQ1 and EQ2. The time period between the end of the reservoir
impoundment and the first two seismic events is about 6.5 years, and it was
modelled again as a pure consolidation analysis stage. For this stage of the
analysis, the BCs (both deformation and hydraulic) were the same as in the
previous consolidation stage.
6.5.4 Comments
The static part of the analysis was performed in order to obtain the correct
stress states prior to the earthquake and ensure that appropriate plasticity
is introduced in the model. This issue is discussed in more detail after the
description of the dynamic analysis, in Section 6.6.
Crest settlement of La Villita dam
0
Recorded
Calculated
Settlement, d [cm]
−5
−10
−15
1968 1969 1970 1971 1972 1973 1974 1975 1976
Time, t [year]
Figure 6.37: Calculated and measured crest settlement history of the dam
crest prior to the EQ events.
Figure 6.37 shows the calculated and recorded crest settlement history of
the dam crest prior to the EQ events which exhibit a good agreement. Figures
6.38 (a) and (b) show the deformed mesh and the vectors of accumulated
displacement respectively at the end of the static part of the analysis. It is
shown that the dam had some minor deformations resulting from the static
stresses (construction, reservoir impoundment and consolidation), but no
237
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
large movements or failure were observed. As shown in the latter figure, the
maximum value of displacement is 0.258m and this is represented by the grey
thick vector, in the core of the dam.
238
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
(a) (b)
Figure 6.38: End of static analysis: (a) Deformed mesh and (b) Vectors of
accumulated displacement of La Villita dam.
239
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
U C
B
L 420
L/H = L/Htot = = = 3.23 (6.2)
Hdam + Halluvium 60 + 70
This implies that canyon effects could be important and a 2D analysis
would be inappropriate as the real problem is stiffer than the corresponding
2D plane strain model. Such a 2D analysis would be suitable for a wide
canyon (L/H > 4), see also Mejia and Seed (1983); Dakoulas and Gazetas
(1987) and relevant discussion in Section 2.4.2. In order to overcome this
soft response of a 2D analysis, the stiffening effect of the canyon geometry
was taken into account by increasing the material stiffness of the dam. A
240
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.15
Acceleration, a [g]
0.1
0.05
−0.05
−0.1
−0.15
0 2 4 6 8 10 12 14 16 18
Time, t [sec]
Figure 6.40: Acceleration time history at the crest of La Villita dam during
EQ2
0.5
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
241
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.6
Acceleration, a [g]
0.4
0.2
−0.2
−0.4
0 10 20 30 40 50 60 70 80
Time, t [sec]
Figure 6.42: Acceleration time history at the crest of La Villita dam during
EQ5
1.2
a
0.8
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
242
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Table 6.7: New calibration parameters for the Logarithmic CNL model, used
for the dynamic analysis.
parametric study was carried out to determine the ratio of the new updated
shear stiffness over the initial stiffness, which would provide the best match
between the calculated and recorded response spectra. It was found that the
shape of the response spectrum was improved (and therefore the prediction
of the fundamental period of the dam was improved as shown later) if the
shear modulus, Gmax (z) was increased by 3.5 times.
The new shear stiffness profile, G∗max (z) was (similarly to the initial pro-
file, Gmax (z)) taken to vary linearly with depth in the dam embankment,
according to Equation 6.3 (where, z is the elevation, i.e. 0 at the base and
60m at the crest of the dam). This was varying from 420000 kPa at the crest
of the dam to 780000 kPa to the base of the dam (see also Equation 6.1 for
the initial variation of Gmax (z)). The corresponding value of G∗max for the
foundation alluvium was assumed constant and equal to 700000 kPa.
243
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
= 700000 kPa was used in the calibration process (see also Section 6.4.3).
Figures 6.54 and 6.56 show the calculated acceleration time histories at
the crest of the dam (Point C in Figure 6.39) for EQ2 and EQ5 respectively
with the new updated values of the G∗max . Likewise, Figures 6.55 and 6.57
show the corresponding response spectra for the two earthquakes respectively.
It may be observed from the former two figures that the new calculated
accelerations were found to be in better agreement with the recorded values.
Moreover, good agreement is also observed in the corresponding response
spectra (Figures 6.55 and 6.57) where the calculated spectral acceleration,
Sa values were closer to the corresponding recorded ones. The previously
observed small and high values of amplification for smaller and larger values
of the period respectively vanished and the calculated spectral accelerations
exhibit large amplifications for the smaller values of the fundamental period.
This shows that a better match of the fundamental period of the dam was
achieved by increasing the material stiffness (shear modulus, G∗max ).
The stiffening effect of narrow canyons was discussed earlier in Chapter
2 and Section 2.4.2. The studies of Mejia and Seed (1983) and Dakoulas
and Gazetas (1987) have quantified this stiffening effect and expressed it as
a function of the canyon length over height, L/H ratio. For La Villita dam,
this L/H ratio is given by Equation 6.2.
According to the study of Dakoulas and Gazetas (1987) (Figure 2.22),
the ratio of the fundamental period of vibration of a dam built in a narrow
canyon, Tn , over that of a dam built in an infinitely wide canyon, Tw , for
L/H = 3.23 and for various shapes of the canyon is given by Equation 6.4.
Tn
= 0.6 ∼ 0.95 (6.4)
Tw
In the present study, the updated value of the shear modulus was taken
as G∗ = 3.5 · G. The shear wave velocity, Vs is given by Equation 6.5 (where,
ρ is the mass density of the material).
s
G
Vs = (6.5)
ρ
244
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.9
0.8
Shear Modulus, G/Gmax []
0.7
0.6
0.5
0.4
0.3
Model − Top
0.2 Model − Middle
Model − Bottom
0.1 Vucetic & Dobry, 1991
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
Figure 6.44: New calibration of the Logarithmic CNL model against the
empirical curves of Vucetic and Dobry (1991) for the clay core - Stiffness
degradation (G/G∗max − γ)
30
Model − Top
25 Model − Middle
Model − Bottom
Vucetic & Dobry, 1991
20
Damping, ξ [%]
15
10
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
Figure 6.45: New calibration of the Logarithmic CNL model against the
empirical curves of Vucetic and Dobry (1991) for the clay core - Damping
(ξ − γ)
245
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.9
0.8
Shear Modulus, G/Gmax []
0.7
0.6
0.5
0.4
0.3
Model − Top
0.2 Model − Middle
Model − Bottom
0.1 Seed et al., 1986
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
Figure 6.46: New calibration of the Logarithmic CNL model against the em-
pirical curves of Seed et al. (1986) for the sand filters - Stiffness degradation
(G/G∗max − γ)
30
Model − Top
25 Model − Middle
Model − Bottom
Seed et al., 1986
20
Damping, ξ [%]
15
10
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
Figure 6.47: New calibration of the Logarithmic CNL model against the
empirical curves of Seed et al. (1986) for the sand filters - Damping (ξ − γ)
246
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.9
0.8
Shear Modulus, G/Gmax []
0.7
0.6
0.5
0.4
0.3
Model − Top
0.2
Model − Middle
0.1 Model − Bottom
Rollins et al., 1998
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
Figure 6.48: New calibration of the Logarithmic CNL model against the em-
pirical curves of Rollins et al. (1998) for the rockfill and foundation alluvium
- Stiffness degradation (G/G∗max − γ)
30
Model − Top
25 Model − Middle
Model − Bottom
Rollins et al., 1998
20
Damping, ξ [%]
15
10
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
Figure 6.49: New calibration of the Logarithmic CNL model against the em-
pirical curves of Rollins et al. (1998) for the rockfill and foundation alluvium
- Damping (ξ − γ)
247
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.9
0.8
Shear Modulus, G/Gmax []
0.7
0.6
0.5
0.4
0.3
Elgamal (1992)
0.2
Model (Vucetic & Dobry, 1991) − Middle
0.1 Model (Seed et al., 1986) − Middle
Model (Rollins et al., 1998) − Middle
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
1000000
Model − Top
900000 Model − Middle
Model − Bottom
800000
700000
Shear Modulus, G [kPa]
600000
500000
400000
300000
200000
100000
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
248
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1000000
Model − Top
900000 Model − Middle
Model − Bottom
800000
700000
Shear Modulus, G [kPa]
600000
500000
400000
300000
200000
100000
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
1000000
Model − Top
900000 Model − Middle
Model − Bottom
800000
700000
Shear Modulus, G [kPa]
600000
500000
400000
300000
200000
100000
0
−4 −3 −2 −1 0 1
10 10 10 10 10 10
Shear strain, γ [%]
249
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.15
0.1
Acceleration, a [g]
0.05
−0.05
−0.1
−0.15
−0.2
0 2 4 6 8 10 12 14 16 18
Time, t [sec]
Figure 6.54: Acceleration time history at the crest of La Villita dam during
EQ2 for the new updated value of G∗max (z).
0.5
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
250
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.6
Acceleration, a [g]
0.4
0.2
−0.2
−0.4
0 10 20 30 40 50 60 70 80
Time, t [sec]
Figure 6.56: Acceleration time history at the crest of La Villita dam during
EQ5 for the new updated value of G∗max (z).
1.2
0.8
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
251
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Therefore, the updated value of the shear wave velocity, Vs∗ is given by
Equation 6.6.
√
Vs∗ = 3.5Vs (6.6)
1
T∗ = √ · T = 0.54T (6.7)
3.5
Therefore, the ratio T ∗ /T = 0.54 is close to the ratio found in the liter-
ature and suggested by Dakoulas and Gazetas (1987) (Tn /Tw = 0.6 ∼ 0.95,
for various shapes of the canyon, see Figure 2.22). This observation confirms
that the calculated stiffening of the narrow canyon is in agreement with ear-
lier work from the literature. However, the stiffening observed in the present
study was found to be slightly larger, i.e. the ratio of T ∗ /T is smaller than
the ratio suggested in the literature. This difference could be attributed to
the fact that the previous theoretical studies assumed the following:
252
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
253
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.04
Acceleration, a [g]
0.02
−0.02
−0.04
−0.06
−0.08
−0.1
0 2 4 6 8 10 12 14 16
Time, t [sec]
0.25
0.2
0.15
0.1
0.05
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
254
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
along with the recorded response spectra. It is evident that the predicted
response in this work is in better agreement with the recorded one than the
response calculated by Elgamal (1992). The response spectrum of Elgamal
(1992), obtained from a 3D shear beam analysis (see Section 6.3), has a
narrower frequency content and higher amplifications at the significant fre-
quencies, whereas the broader frequency content of the spectrum from the
present study matches better the low frequency spectral ordinates. However,
it should be noted that none of the two calculated response spectra (i.e. from
this study and from Elgamal (1992)) match the spectral ordinates for higher
frequencies (smaller periods).
Those high accelerations for high frequencies in the response spectrum of
the recorded motion are believed to be originated from the high peak values of
acceleration observed in the recorded acceleration time-history at the crest of
the dam for EQ5. These high peaks are believed to be due to a localised slope
failure, as suggested by previous researchers (see also Section 6.3), which
was not predicted by the FE analysis in this study. However, if the high
frequencies are filtered from the recorded accelerations at the crest, a better
255
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
256
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.3
Acceleration, a [g]
0.2
0.1
−0.1
−0.2
−0.3
0 10 20 30 40 50 60 70 80
Time, t [sec]
Figure 6.61: Acceleration time history at the crest of La Villita dam during
EQ5. Comparison between filtered recorded and calculated response.
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
257
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
during EQ5 (around 30cm). On the latter figure, the value and orientation
of maximum displacement (which has a value of 0.036m) is represented by
a grey vector and it is located at the upstream dam slope. This means that
the deformations resulting from the earthquake events are very small.
Figure 6.65 shows the contours of stress level, S, in the dam during EQ5.
The stress level, S is defined as the ratio of the current deviatoric stress, Jc ,
over the value of the deviatoric stress at yield conditions, Jy , at the same
value of the mean efective stress, p0 . Therefore, S, takes values from 0 to 1,
and shows how close the stress state of the soil is to yielding. It is shown
that the values of stress level go up to 0.8 (contour B) at some places within
the upstream rockfill, the downstream alluvium and the downstream dam
slope. However, S is not very close to 1 and therefore the soil in the dam
is generally found not to be at the yielding stress state after the end of the
earthquake.
Figures 6.66 and 6.67 show the calculated vertical displacements at the
crest (point C in Figure 6.39) during the EQ2 and EQ5 seismic events. It may
be observed that the magnitude of the calculated settlements is extremely
small compared to the recorded values for both EQs (around 2cm and 30cm
for EQ2 and EQ5 respectively). Furthermore, it should be noted that accord-
ing to the global deformations of the dam during EQ5 (Figures 6.63 (a) and
(b)), the positive values of settlements for EQ5 are originated from the fact
that uneven displacements occur on the crest of the dam, i.e. the upstream
part of the crest experiences negative (i.e. downwards) whereas the down-
stream part of the crest experiences positive (i.e. upwards) displacements.
However, note that all crest displacements are extremely small (i.e. < 2cm).
Figures 6.68 and 6.69 show the calculated vertical displacements at the
upstream corner of the crest (point U in Figure 6.39). As discussed earlier,
according to Figures 6.63 (a) and (b), the main deformations in the dam were
concentrated in the upstream rockfill and include the upstream corner of the
crest (point U in figure 6.39). It is shown that the upstream part of the crest
experiences higher values of settlement, perhaps due to those movements of
the upstream rockfill, but still their magnitude is very small compared to the
recorded values. It should also be noted that the recorded values refer to the
258
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
(a) (b)
Figure 6.63: End of the seismic analysis (after EQ5): (a) Accumulated (for
the whole static and dynamic analysis) deformed mesh and (b) Vectors of
accumulated displacement of La Villita dam.
259
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
(a) (b)
Figure 6.64: End of the seismic analysis (after EQ5): (a) Sub-accumulated
(i.e. only for EQ5) deformed mesh and (b) Vectors of sub-accumulated dis-
pacement of La Villita dam.
260
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
BA
C
A
B
A
B
B
A
A
B
A
A
B
C
B
D
A
Figure 6.65: End of the seismic analysis (after EQ5): Stress level, S, in La
Villita dam.
261
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.01
Settlement, d [cm]
−0.01
−0.02
−0.03
−0.04
0 1 2 3 4 5 6 7
Time, t [sec]
0.2
0.15
Settlement, d [cm]
0.1
0.05
−0.05
−0.1
−0.15
−0.2
0 10 20 30 40 50 60 70
Time, t [sec]
262
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0
Settlement, d [cm]
−0.02
−0.04
−0.06
−0.08
−0.1
−0.12
0 1 2 3 4 5 6 7
Time, t [sec]
−0.2
−0.4
Settlement, d [cm]
−0.6
−0.8
−1
−1.2
−1.4
−1.6
−1.8
0 10 20 30 40 50 60 70
Time, t [sec]
Figure 6.69: Vertical displacement history at the upstream corner of the crest
(Point U in Figure 6.39) of La Villita dam during EQ5.
263
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
downstream side of the crest, where the monitoring instrument was installed.
264
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
B
A
A
B
D
C
B
A
B
C
A
D
B
C
A
B
A
C
D
B
A
E
F
Figure 6.70: End of the seismic analysis (after EQ5): Accumulated total
deviatoric strain, Ed in La Villita dam.
265
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
CC
UR DR
DA
Figure 6.71: Sketch of the dam showing the locations of the monitoring
elements: upstream rockfill (UR), clay core (CC), downstream rockfill (DR)
and downstream alluvium (DA).
Shear stress−strain response at the downstream alluvium of La Villita dam during EQ5
250
200
150
100
Shear stress, τ [kPa]
50
−50
−100
−150
−200
−250
−0.04 −0.02 0 0.02 0.04 0.06 0.08 0.1
Shear strain, γ [%]
266
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Shear stress−strain response in the clay core of La Villita dam during EQ5
150
100
Shear stress, τ [kPa]
50
−50
−100
−0.015 −0.01 −0.005 0 0.005 0.01 0.015 0.02 0.025
Shear strain, γ [%]
Figure 6.73: Shear stress-strain response in the clay core (CC) of La Villita
dam during EQ5
Shear stress−strain response in the upstream rockfill of La Villita dam during EQ5
150
100
Shear stress, τ [kPa]
50
−50
−100
−0.02 −0.01 0 0.01 0.02 0.03 0.04
Shear strain, γ [%]
267
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Shear stress−strain response in the downstream rockfill of La Villita dam during EQ5
100
50
Shear stress, τ [kPa]
−50
−100
−150
−0.04 −0.03 −0.02 −0.01 0 0.01 0.02 0.03
Shear strain, γ [%]
Shear stress−strain response at the downstream alluvium of La Villita dam during EQ2
80
60
40
Shear stress, τ [kPa]
20
−20
−40
−60
−0.015 −0.01 −0.005 0 0.005 0.01 0.015
Shear strain, γ [%]
268
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Shear strain time−history in the upstream rockfill of La Villita dam during EQ5
0.04
0.03
Shear strain, γ [%]
0.02
0.01
−0.01
−0.02
0 10 20 30 40 50 60 70
Time, t [sec]
Shear strain time−history in the downstream rockfill of La Villita dam during EQ5
0.03
0.02
0.01
Shear strain, γ [%]
−0.01
−0.02
−0.03
−0.04
0 10 20 30 40 50 60 70
Time, t [sec]
269
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
250
Start of Reservoir Impounding
Static part of analysis
200
150
100
50
0
0 50 100 150 200 250 300 350 400
Mean effective stress, p‘ [kPa]
Figure 6.79: Stress paths in the upstream rockfill (UR) of La Villita dam
200
150
100
50
0
0 100 200 300 400 500 600
Mean effective stress, p‘ [kPa]
Figure 6.80: Stress paths in the downstream rockfill (DR) of La Villita dam
270
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
271
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
in smaller values of the mean effective stress, p0 , and brings the stress path
closer to the yield surface. That is why more plasticity is introduced in the
upstream rockfill.
This seems to be a drawback of the constitutive models adopted in this
study, which introduce plasticity only when the yield surface is engaged. In
contrast, kinematic hardening models (such as Grammatikopoulou (2004))
are able to introduce plasticity much earlier, as they have a much smaller
elastic area which is more realistic. This is in agreement with earlier obser-
vations of Kontoe et al. (2011) who compared constitutive assumptions of
different levels of complexity. They commented that the prediction of plastic
strains by the constitutive assumptions used in their study was dependent
on the proximity of the initial stress state of the soil (i.e. at the beginning
of the earthquake) to the yield surface. This highlights the importance of
appropriate modelling of the construction sequence before the earthquake.
Finally, considering the results of this analysis, the constitutive assump-
tions used in this study seem to be incapable of taking account of the effects
of the previous stress history (i.e. a sequence of seismic events) and the re-
sulting development of plasticity properly. This issue is addressed in more
detail in the next chapter and particularly in Section 7.2.
6.6.3 Comments
From the dynamic analysis of La Villita dam, it was found that the three-
dimensional geometry of the surrounding canyon has an influence on the
dynamic characteristics (fundamental period and amplification of accelera-
tions) of the dam. As the current work did not involve 3D analysis, the
stiffening effect in a two-dimensional analysis was modelled by increasing
the material stiffness in order to explicitly take account of the additional
geometric stiffness from the canyon.
The dynamic response of the dam was well captured as evidenced by
the comparison of the accelerations at the crest and the resulting response
spectra. The calculated response spectrum for EQ5 at the crest of the dam
from this work was compared to that of Elgamal (1992) who performed three-
272
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
dimensional shear beam analysis of La Villita dam and was found to have a
better agreement with the recorded response spectrum.
However, the high peak values of acceleration observed in the recorded
acceleration and already attributed by previous researchers to a localised fail-
ure close to the monitoring instrument were not captured. Such a localised
failure was not predicted, from the inspection of the plots of the deformed
mesh and vectors of displacement after the earthquake. Nevertheless, when
the high frequencies (originating from the localised slip failure) from the
recorded accelerations were filtered, an excellent agreement was obtained be-
tween the calculated and recorded accelerations and response spectra. This
shows that the dynamic response of the dam was well captured and the fre-
quency content of the resulting calculated accelerograph at the crest matched
that of the recorded accelerograph.
The calculated vertical displacements at the crest of the dam were not
comparable to the recorded values. As mentioned earlier, a localised slip
failure that would result into high values of settlements was not predicted
and therefore no high values of settlements were observed. Minor deforma-
tion in the upstream rockfill showed some higher values of settlement at the
US corner of the crest, but still not large enough to be comparable to the
recorded crest settlement values. It should also be noted that both the ob-
served asymmetry in the recorded crest accelerations and the location of the
monitoring instrument suggest a localised failure on the downstream side of
the crest.
It is suggested that more research is required to be carried out in order
to further investigate the pattern of the recorded settlements. It is believed
that the failure of this study to predict the recorded displacements could be
due to:
Regarding the first possible reason, presumably some minor localised fail-
ure at the position close to the monitoring instrument may had happened
273
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
before which may had created a local discontinuity. This could have origi-
nated possibly from a previous seismic event or even from some construction
processes. Besides, it should be reminded that there is a public road on top
of the crest of La Villita dam and perhaps heavy vehicle traffic may have
caused some movements in the soil. If such a discontinuity exists, it will
form a weak zone in the dam and will be sensitive to seismic loads. As far
as modelling is concerned, the adopted procedure in this work assumed a
continuous soil profile, since no weak zones were introduced with interface
elements, such as those of Gazetas and Uddin (1994).
Regarding the second possible reason, perhaps the constitutive modelling
assumptions adopted in this work, while adequate to capture the dynamic
response of the dam (accelerations and fundamental period); were not the
most appropriate for predicting permanent movements. As observed from the
stress paths (Figures 6.79 and 6.80), plasticity (and permanent strains) were
introduced only when the yield surface was engaged, whereas the strains were
recoverable for the rest of the analysis. No pre-yield plasticity was introduced
although this is more consistent with the real behaviour of soils. Such a mod-
elling concept would introduce more strains and hence more movements (i.e.
larger settlements) should it was adopted. Moreover, kinematic hardening
and softening (i.e. change of the yield surface) would perhaps indicate some
more plasticity.
Although the yield function adopted (Mohr-Coulomb) was allowed to
change with the Lode’s angle, θ (i.e. the relative magnitude of the principal
stresses, σ1 , σ2 and σ3 , see Section 3.7.5 and Equations 3.32 - 3.35), this
change was perhaps still not sufficient. Better results could perhaps be ob-
tained should a more advanced constitutive model of the Critical State type
(Schofield and Wroth, 1968), such as the Modified Cam Clay (Roscoe and
Burland, 1968) was employed, as proved by earlier studies on embankments
(Zdravković et al., 2002).
A further improvement could also be the adoption of a kinematic harden-
ing (also known as “bubble”) model (KHM), such as that of Grammatikopoulou
(2004) which allows introduction of pre-yield plasticity and thus it pro-
vides a better prediction of deformations (Grammatikopoulou et al., 2007,
274
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
2008). This advanced model has already been applied in the static (Gram-
matikopoulou et al., 2008) and seismic analysis of tunnels (Kontoe et al.,
2009, 2011). Those studies showed that this KHM is able to provide bet-
ter predictions of strains and displacements (as compared to recorded values
from the field) than the constitutive approach used in this study, i.e. a small-
strain stiffness degradation model (such as the Logarithmic CNL (Taborda,
2011) or the Jardine et al. (1986) small-strain model) combined with a simple
failure criterion (Mohr-Coulomb or Modified Cam Clay).
Finally, a parametric analysis with different yield strengths of the rockfill
(i.e. different values of the cohesion, c0 and the angle of shearing resistance,
φ0 ) could investigate the dependency of the magnitude of displacements on
the yield strength of the rockfill.
275
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
and the alluvium (10−7 m/s) were adopted and used in the analysis. How-
ever, because of the lack of relevant data for the rockfill shoulders, these parts
of the dam were assumed to have fully drained materials during the entire
analysis.
Moreover, available field measurements were very limited. Out of the
six earthquake that were monitored between 1975-1985 (see Table 6.1), only
two input motions (EQ2 and EQ5) are useful for numerical analysis (as the
rest were incomplete). Also, for those two seismic events, the full input
bedrock motion of EQ2 was not available and also the acceleration record
at the downstream berm for EQ5 was not available at all. For this reason,
only two seismic events were used in the analysis and the remaining events
(EQ1, EQ3, EQ4 and EQ6) were not included in the analysis. An ideal back
analysis of the response of the dam would require the consideration of all
the earthquakes in order to better assess the effect of loading history on the
seismic response of earth dams.
6.8 Conclusions
This chapter describes the work carried out related to the numerical analysis
of the static and dynamic behaviour of La Villita dam. This dam is consid-
ered to be a well-documented case study because of available information,
i.e. material properties and monitored response.
The dam was analysed in both static and dynamic conditions, considering
the nonlinear elasto-plastic soil behaviour. Two-dimensional plane-strain
coupled static and dynamic finite element analyses were performed.
The (dynamic) stiffening effect of the 3D narrow canyon geometry was
taken into consideration in a 2D dynamic analysis by increasing the material
stiffness of the dam. The computed accelerations and associated response
spectra were found to be in good agreement with the corresponding avail-
able recorded values. Also, a better prediction of the response spectra was
achieved than that of Elgamal (1992).
However, the high peak values of acceleration observed in the recorded
acceleration (which were attributed to a localised failure) were not captured
276
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
as the current FE did not predict any localised failure close to the dam crest
and the calculated crest settlements were much smaller than the recorded
values. A possible pre-existing discontinuity in the downstream embankment
crest soil profile would form a weak zone which would be very sensitive to
any seismic loads. Such a discontinuity was not artificially modelled in this
work (through the introduction of interface elements in the soil profile).
The soil in the upstream rockfill experienced higher values of shear strains
and larger deformations than that in the downstream rockfill. This was due
to the reservoir impoundment during which the stress paths were directed
closer the yield surface, because of the reduction in the mean effective stress,
p0 . Finally, it is believed that an improved prediction of the deformations of
the dam may be obtained by using more advanced constitutive models (such
as kinematic hardening) able to predict plasticity more realistically.
277
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Chapter 7
PARAMETRIC SEISMIC
ANALYSIS OF LA VILLITA
DAM
7.1 Introduction
This chapter contains several parametric studies performed to assess the
effects of different modelling approaches on the predicted seismic response of
earth dams. The already studied in Chapter 6 case of La Villita dam is used
as the basis for this investigation.
The examined issues are (a) the effects of any previous seismic activ-
ity prior to the studied earthquake, (b) the relative merits of coupled and
uncoupled dynamic analysis, (c) foundation-dam interaction and modelling
approaches, (d) nonlinear reservoir-dam interaction and (e) stiffness inho-
mogeneity. Previous results from the literature are also discussed (where
available) and compared to the outcomes of this study. Finally, comments
are made regarding the applicability of the adopted modelling technique and
constitutive assumptions.
278
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
279
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.2
Acceleration, a [g]
0.1
−0.1
−0.2
−0.3
0 10 20 30 40 50 60 70
Time, t [sec]
Figure 7.1: Acceleration response at the crest of La Villita dam during EQ5
- Comparison between Cases A and B
0.2
Acceleration, a [g]
0.1
−0.1
−0.2
−0.3
0 10 20 30 40 50 60 70
Time, t [sec]
Figure 7.2: Acceleration response at the crest of La Villita dam during EQ5
- Comparison between Cases B and C
280
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
281
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.2
0.1
Displacement, d [cm]
−0.1
−0.4
0 10 20 30 40 50 60 70
Time, t [sec]
Figure 7.5: Vertical displacements at the crest of La Villita dam during EQ5
- Comparison between Cases A and B
0.4
0.3
Displacement, d [cm]
0.2
0.1
−0.2
0 10 20 30 40 50 60 70
Time, t [sec]
Figure 7.6: Vertical displacements at the crest of La Villita dam during EQ5
- Comparison between Cases B and C
282
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
283
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Shear strain time−history in the downstream alluvium of La Villita dam during EQ5
0.1
0.08
0.06
Shear strain, γ [%]
0.04
0.02
−0.04
0 10 20 30 40 50 60 70
Time, t [sec]
Shear strain time−history in the downstream alluvium of La Villita dam during EQ5
0.1
0.08
0.06
Shear strain, γ [%]
0.04
0.02
−0.06
0 10 20 30 40 50 60 70
Time, t [sec]
284
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Pore water pressure time−history in the downstream alluvium of La Villita dam during EQ5
100
B − EQ5 after EQ2
80 A − EQ5 without previous EQ
60
Pore water pressure, u [kPa]
40
20
−20
−40
−60
−80
0 10 20 30 40 50 60 70
Time, t [sec]
Pore water pressure time−history in the downstream alluvium of La Villita dam during EQ5
100
B − EQ5 after EQ2
80 C − EQ5 after EQ5
60
Pore water pressure, u [kPa]
40
20
−20
−40
−60
−80
0 10 20 30 40 50 60 70
Time, t [sec]
285
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
500
Start of Reservoir Impounding
Static part of analysis
400
300
200
Scenario AB
100
0
0 100 200 300 400 500 600 700 800 900
Mean effective stress, p‘ [kPa]
Figure 7.11: Stress paths in the downstream alluvium (DA) of La Villita dam
during EQ5 - Scenario B
500
400
300
200
100 A
Scenario B
0
0 100 200 300 400 500 600 700 800 900
Mean effective stress, p‘ [kPa]
Figure 7.12: Stress paths in the downstream alluvium (DA) of La Villita dam
during EQ5 - Scenario A
286
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
500
Start of Reservoir Impounding
Static part of analysis
400
300
200
100 Scenario C
0
0 100 200 300 400 500 600 700 800 900
Mean effective stress, p‘ [kPa]
Figure 7.13: Stress paths in the downstream alluvium (DA) of La Villita dam
during EQ5 - Scenario C
and it shows how close the stress state of the soil is from the yield condition.
Comparing all three plots, it should be mentioned that only minor differences
exist for the three scenarios considered in this study. However, one could
notice that slightly larger values of S (up to 0.9) are observed in the upstream
rockfill for Scenario C (Figure 7.14 (b)) but they are in a very localised area
around the downstream berm. This was expected as the latter scenario
involved higher values of acceleration and therefore more engagement of the
yield surface.
7.2.4 Comments
It was shown that the effects of the previous loading history were gener-
ally found to be very small in the case of La Villita dam. For the cases
considered, no difference was observed in the crest accelerations and the as-
sociated response spectra. However, minor differences were observed in the
crest settlements, the induced shear strains and the pore water pressures. It
is believed that this happened for three reasons: (a) none of the scenarios
considered a major event which would have induced some significant failure
287
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
B C
C
B
A
A
B
A
A
B
A
A
A
A
B
A
A
B
B
C
B
A
A
B
A
A
A
A
A
B
B
C
D
A
B
C
D
B
E
A
Figure 7.14: Contours of stress level, S for (a) Case A - EQ5 without any
previous earthquake, and (b) Case C - EQ5 after a large earthquake, which
is EQ5.
288
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
in the dam, prior to EQ5, (b) the nature of the materials of La Villita dam
and (c) the constitutive assumptions used in this study.
As far as the first reason is concerned, should an extremely large seismic
event (i.e. larger than EQ5) have been modelled prior to EQ5, the higher
values of plasticity would have induced large deformations and altered the
geometry of the dam significantly. In that case, the wave propagation would
have resulted in a different dynamic response of the dam structure. However,
none of the three scenarios considered such an extremely large seismic event,
and therefore the analysis of EQ5 in all three scenarios started from similar
conditions. Also, a large displacement analysis would be required to capture
the effects of a significantly altered dam geometry.
As far as the second issue is concerned, it should be noted that the clay
core constitutes of only a small part of the dam, whereas the rest consists of
sand filters and rockfill shells which were modelled (and expected to behave)
as drained materials. Therefore, only the clay core and the alluvial foun-
dation were considered as consolidating materials and were expected to be
more affected by the previous stress history. The effects of the previous stress
history would be more highlighted if all the parts of the dam were modelled
as consolidating materials (with appropriate values of permeability).
Moreover, as far as the last issue is concerned, the adopted CNL model
with a simple Mohr-Coulomb failure criterion is perhaps not adequate in
modelling rigorously the previous stress history. This was also discussed
earlier in Chapter 6, and was found inadequate to capture the recorded crest
settlements of La Villita dam. The inability of the adopted constitutive
assumptions to model pre-yield plasticity could be the main reason for the
minor differences between the three scenarios considered here. Should a
kinematic hardening model be adopted, more differences might be observed,
as plasticity would be introduced earlier in the analysis, and not only during
the limited occasions that the stress paths approached the Mohr-Coulomb
yield surface.
289
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
290
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
The investigation was concerned with both EQ2 and EQ5 seismic events.
Therefore, two separate uncoupled analyses were carried out, one for EQ2
and one for EQ5 and they both started immediately after the static analysis.
The static analysis was the same as the one used for the previous analyses,
i.e. Chapter 6 and Section 7.2. Therefore, the uncoupled analysis for EQ2
(described in this section) is compared with the coupled analysis of EQ2
described in Chapter 6. The uncoupled analysis for EQ5 (described in this
section) is comparable to the analysis of Scenario A of Section 7.2, as they
both analyse EQ5 without any previous earthquake.
291
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.05
Acceleration, a [g]
−0.05
−0.1
−0.15
−0.2
0 1 2 3 4 5 6 7
Time, t [sec]
0.5
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
Figure 7.16: Response spectra (damping, ξ = 5%) the crest of La Villita dam
during EQ2.
292
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.2
Acceleration, a [g]
0.1
−0.1
−0.2
−0.3
0 10 20 30 40 50 60 70
Time, t [sec]
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
Figure 7.18: Response spectra (damping, ξ = 5%) the crest of La Villita dam
during EQ5.
293
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
40 Coupled Analysis
Uncoupled Analysis
20
Elevation, z [m]
−20
−40
−60
−80
1 2 3 4 5 6 7 8 9 10
Displacement, d [m] −3
x 10
40
20
Elevation, z [m]
−20
−40
294
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
40
20
Elevation, z [m]
−20
−40
20
Elevation, z [m]
−20
−40
−60
−80
0.1 0.15 0.2 0.25 0.3 0.35
Acceleration, a [g]
295
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
7.3.4 Comments
The dynamic response of the dam, if coupled or uncoupled analysis is per-
formed, was found to be almost identical, as evidenced by the calculated
accelerations at the dam crest, the associated response spectra and displace-
ments. Also, the dynamic behaviour of the soil was found to be very similar,
considering the shear strains (shear-stress strain response and shear stress
time-history). However, small differences were observed in the calculated
pore water pressures, where the uncoupled analysis predicted higher fluctu-
ations of the pore pressure.
It was shown that although an uncoupled dynamic analysis may predict a
very similar seismic response of an earth dam, it may however have different
pore water pressures than a coupled dynamic analysis. The fact that the dy-
namic response of the dam is very similar could be attributed to two reasons:
(a) the consolidating materials (i.e. those for which a coupled formulation
was adopted in the coupled analysis, namely clay core and foundation al-
luvium, and which are expected to show some difference if a coupled or an
uncoupled undrained analysis is performed) consisted only a small part of the
296
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Shear stress−strain response in the clay core of La Villita dam during EQ5
200
Uncoupled analysis
150
Shear stress, τ [kPa]
100
50
−50
−0.01 −0.005 0 0.005 0.01 0.015 0.02 0.025 0.03
Shear strain, γ [%]
Figure 7.23: Shear stress-strain response in the clay core of La Villita dam
during EQ5
Shear stress−strain response in the clay core of La Villita dam during EQ5
200
Coupled analysis
150
Shear stress, τ [kPa]
100
50
−50
−100
−0.01 −0.005 0 0.005 0.01 0.015 0.02 0.025 0.03
Shear strain, γ [%]
Figure 7.24: Shear stress-strain response in the clay core of La Villita dam
during EQ5
297
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Shear strain time−history in the clay core of La Villita dam during EQ5
0.03
0.025
0.02
Shear strain, γ [%]
0.015
0.01
0.005
Figure 7.25: Shear strain time-history in the clay core of La Villita dam
during EQ5
Pore water pressure time−history in the clay core of La Villita dam during EQ5
100
Coupled analysis
50 Uncoupled analysis
Pore water pressure, u [kPa]
−50
−100
−150
−200
−250
−300
0 10 20 30 40 50 60 70
Time, t [sec]
Figure 7.26: Pore water pressure time-history in the clay core of La Villita
dam during EQ5
298
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
dam geometry, and (b) the values of the permeability of the consolidating
materials adopted in this study were not adequately large to allow signif-
icant seepage and highlight the differences of a coupled and an uncoupled
undrained analysis.
Regarding the first reason, the differences between the two analyses (cou-
pled and uncoupled) would be more visible if all the dam materials were
modelled as consolidating in the coupled analysis. In that case, appropriate
values of permeability (e.g. for rockfill) should be adopted so that seepage
and dissipation of water pressures are correctly modelled.
As far as the second reason is concerned, larger differences (at least for the
pore water pressures) would be expected if the permeability of the materials
was larger and therefore allowed for some seepage in the dam and changes in
the water pressures. Therefore, it is suggested that in order to fully examine
the differences of performing a coupled and an uncoupled analysis, a dam
consisting entirely of consolidating materials should be considered and more
permeable materials (such as rockfill shoulders) should be included in the
coupled formulation.
299
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
were: (a) free-field (FF), (b) downstream dam toe (DDT) and (c) dam base
(DB).
DB
DDT FF
Figure 7.27: Mesh of the dam showing the locations of the monitoring points:
Dam base (DB), Downstream dam toe (DDT) and Free-field (FF).
Figures 7.28 and 7.29 show the acceleration time-histories at the three
monitoring points (FF, DDT and DB) for EQ2 and EQ5 respectively, whereas
Figures 7.30 and 7.31 show the associated response spectra for the two earth-
quakes. From the figures of the accelerations, it may be observed that the
records at the FF (for both earthquakes) exhibit a higher frequency content,
whereas the records at the DB show larger periods of vibration.
This can be more clearly observed in the plots of the response spectra.
For higher periods (such as T=0.7s), the DB record has the highest spectral
ordinates whereas the FF record has smaller ordinates. At smaller values
of the fundamental period (such as T=0.2s), the opposite is observed, i.e.
the DB record has the smallest and the FF record has the largest spectral
ordinates. In both cases, i.e. at T = 0.2 and 0.7s, the DDT has spectral
ordinates between the two other records.
This means that for the DB record results in higher acceleration values at
longer fundamental periods, whereas the FF record at shorter fundamental
periods. Therefore, the DB record has a high-frequency content as compared
to the FF record which exhibits a low-frequency content. Practically, this
means that the presence of the dam results in a softer soil layer response
than the free-field soil layer response. On one hand, the presence of a large
structure (such as a massive earth dam) on top of a uniform soil layer would
provide some additional geometric stiffness to the soil layer, as the upper
300
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.05
−0.05
−0.1
0 1 2 3 4 5 6 7
Time, t [sec]
Acceleration response − EQ2
0.1
EQ2 − Dam base (DB)
Acceleration, a [g]
0.05
−0.05
−0.1
0 1 2 3 4 5 6 7
Time, t [sec]
Acceleration response − EQ2
0.1
Acceleration, a [g]
0.05
−0.05
EQ2 − Downstream dam toe (DDT)
−0.1
0 1 2 3 4 5 6 7
Time, t [sec]
301
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.1
−0.1
EQ5 − Free−field (FF)
−0.2
0 10 20 30 40 50 60 70
Time, t [sec]
Acceleration response − EQ5
0.3
EQ5 − Dam base (DB)
0.2
Acceleration, a [g]
0.1
−0.1
−0.2
0 10 20 30 40 50 60 70
Time, t [sec]
Acceleration response − EQ5
0.2
Acceleration, a [g]
0.1
−0.1
EQ5 − Downstream dam toe (DDT)
−0.2
0 10 20 30 40 50 60 70
Time, t [sec]
302
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.35
0.3
0.25
0.2
0.15
0.1
0.05
0
0 0.5 1 1.5 2 2.5 3
Period, T [sec]
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3
Period, T [sec]
303
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
boundary is no longer free to move (i.e. shear stresses being equal to zero,
τ = 0), but there are some stresses from the base of the dam. On the other
hand, an oscillator with a larger mass exhibits a softer response, i.e. it has
a larger fundamental
q
period, T as the latter is proportional to the mass, m
(T = 2π m/k, where, k is the stiffness of the oscillator). This observed
softening in the response is in agreement with earlier theoretical studies of
the shear response of a uniform elastic soil layer which was found to be stiffer
when a surface loading was applied (Ambraseys, 1959b).
304
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
For each of the two seismic events, EQ2 and EQ5, the first input motion,
“Bedrock”, was the same as the one used in Chapter 6 and shown in Figure
6.6, for both EQ2 and EQ5, i.e. the real acceleration record monitored at
the dam site. The second input motion, “Free-field” is the one calculated
from the uncoupled analysis described in Section 7.3 and discussed in Section
7.4.1. It is shown in Figures 7.28 and 7.29 for EQ2 and EQ5 respectively.
Figures 7.32 and 7.33 show the acceleration time histories at the crest
of the dam for the two scenarios considered, for EQ2 and EQ5 respectively.
Figures 7.34 and 7.35 show the associated response spectra for the two earth-
quakes, along with that from the original analysis of the dam-foundation
system discussed in Section 7.4.1.
From the acceleration time-histories at the crest of the dam (Figures 7.32
and 7.33), it is shown that the dam accelerations for the FFI (Case A) (for
both earthquakes) exhibit higher values of acceleration and higher frequency
content, whereas the dam accelerations for the BRI (Case B) show smaller
values of acceleration and larger periods of vibration.
The latter observation is also confirmed from the graphs of the response
spectra. For EQ5 (Figure 7.35), at higher values of the period (such as
T=0.7s), the BRI record (black solid line) has higher spectral ordinates
whereas the FFI record (grey solid line) has smaller ordinates. At smaller
values of the fundamental period (such as T=0.4s), the opposite is observed,
i.e. the BRI record has smaller spectral ordinates than the FFI. This was
expected as the input motion used in the BRI analysis was richer in lower
frequencies than the input motion used in the FFI analysis, whereas the lat-
ter was richer in higher frequencies than the former, as shown in the previous
section. It should however be noted however that this trend is not very clear
for EQ2.
Moreover, on the same graphs, the response spectra at the crest of the
dam, of the BRI analysis (BR input motion) (black solid line) and the original
analysis of a dam-foundation system, DFS (grey dashed lined) are presented.
It is shown that the BRI analysis with rigid foundation resulted in higher
spectral amplifications for higher frequencies, than the original analysis on
the soil layer, whereas the latter resulted in higher amplification for smaller
305
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.1
0.05
Acceleration, a [g]
−0.05
−0.1
−0.2
0 1 2 3 4 5 6 7
Time, t [sec]
0.3
0.2
Acceleration, a [g]
0.1
−0.1
−0.2
−0.3
306
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.6
0.4
Acceleration, a [g]
0.2
−0.2
−0.4
−0.6
EQ5 − Dam on rigid foundation: Bedrock Input (BRI)
−0.8
0 10 20 30 40 50 60 70
Time, t [sec]
0.5
Acceleration, a [g]
−0.5
−1
EQ5 − Dam on rigid foundation: Free−field Input (FFI)
−1.5
0 10 20 30 40 50 60 70
Time, t [sec]
307
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1.4
1.2
0.8
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3
Period, T [sec]
Figure 7.34: Response spectra (damping, ξ = 5%) under various input mo-
tions for a dam on a rigid foundation - EQ2
0
0 0.5 1 1.5 2 2.5 3
Period, T [sec]
Figure 7.35: Response spectra (damping, ξ = 5%) under various input mo-
tions for a dam on a rigid foundation - EQ5
308
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
frequencies. This was expected, as the same input record was used in the
analysis of two different systems (a dam on a rigid base and a dam-foundation
system) which had different fundamental periods and therefore amplification
of different frequencies.
7.4.3 Comments
This section examined dam-foundation interaction effects. Considering a
dam-foundation system, it was shown that the presence of a large embank-
ment dam on top of a uniform soil layer (foundation) changes the dynamic
characteristics of the layer and results in different values of amplification
and for different values of the loading frequency. It was found that smaller
amplification at the top of the soil layer occurs closer to the dam, and the
accelerations tend to contain smaller frequencies (i.e. the soil layer exhibits
a softer response).
Moreover, it was found that significantly different amplifications and fre-
quency contents are obtained if a dam is analysed on a rigid base. This was
true for a number of input motions examined, i.e. for using (a) the same
input motion (bedrock) and (b) the free-field motion (of the soil layer). This
raises issues regarding the validity of ignoring the foundation soil layer and
the use of a “de-coupled” dynamic analyses of earth dams.
Namely, considering Case A, it is incorrect to use the bedrock acceler-
ation at the base of a dam founded on a soft foundation soil layer, as this
would ignore the influence of the foundation layer, i.e. softening of the re-
sponse of the dam structure and change of the dynamic characteristics of
the earthquake (amplification and frequency content). Moreover, regarding
Case B, it is again inappropriate to consider the dam and the foundation
soil layer as separate independent domains in a de-coupled analysis, i.e. to
perform a site response analysis of the soil layer and then use the calculated
acceleration time-history at the top of the layer as input in a dynamic anal-
ysis of a dam founded on a rigid base. Therefore, it is concluded that a
dam-foundation system should be analysed together in order to capture the
effects of dam-foundation interaction appropriately and care should be taken
309
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
310
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
342m
60m
54m
70m
582m
Figure 7.36: 2D Finite Element mesh used for the analysis of La Villita dam
boundary during the dynamic analysis was the Viscous boundary condition
(see Sections 3.7 and 4.3) and the values of the viscosities were obtained from
the elastic properties of the material to which they were applied, i.e. water
(see Section 4.3). The hydraulic boundary conditions were the same as those
described in Chapter 6.
311
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.05
Acceleration, a [g]
−0.05
−0.1
−0.15
−0.2
0 1 2 3 4 5 6 7
Time, t [sec]
0.5
a
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
312
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.2
Acceleration, a [g]
0.1
−0.1
−0.2
−0.3
0 10 20 30 40 50 60 70
Time, t [sec]
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
313
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
20
Elevation, z [m]
−20
−40
−60
−80
1 2 3 4 5 6 7 8 9 10
Displacement, d [m] −3
x 10
20
Elevation, z [m]
−20
−40
−60
−80
0.13 0.135 0.14 0.145 0.15 0.155 0.16 0.165 0.17 0.175
Displacement, d [m]
314
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
20
Elevation, z [m]
−20
−40
−60
−80
0.04 0.06 0.08 0.1 0.12 0.14 0.16
Acceleration, a [g]
20
Elevation, z [m]
−20
−40
−60
−80
0.1 0.15 0.2 0.25 0.3 0.35
Acceleration, a [g]
315
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
slightly smaller values of the accelerations at the dam crest as shown for the
acceleration time-histories and the associated response spectra. Moreover,
smaller values were calculated for the horizontal accelerations and displace-
ments in the dam core for the analyses considering the reservoir. Finally,
it should be noted that the effects of reservoir-dam interaction were more
pronounced for EQ5 which had higher intensity and resulted in higher values
of accelerations in the dam structure.
It may be concluded that the presence of the reservoir “damped” the
response of the dam and resulted in slightly smaller values of accelerations.
The effects of reservoir-dam interaction were investigated in Chapter 5 for
linear elastic dams founded on a rigid foundation. A direct comparison of
this study with the results of Chapter 5 cannot be made as (a) Chapter 5
considered dams on a rigid base, whereas La Villita dam is founded on a
soil layer and (b) the fundamental period of the dam cannot be accurately
estimated as nonlinear behaviour was taken into consideration. However, the
results of this section agree qualitatively with those of Chapter 5 showing that
RDI effects are insignificant in earth dams.
316
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
400
200
−200
−400
0 1 2 3 4 5 6 7
Time, t [sec]
40
Elevation, z [m]
30
20
10
0
−5 0 5 10 15 20 25
Peak Hydrodynamic Pressure, Pdyn [kPa]
317
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1000
500
−500
−1000
−1500
0 10 20 30 40 50 60 70
Time, t [sec]
40
Elevation, z [m]
30
20
10
0
−10 0 10 20 30 40 50 60 70
Peak Hydrodynamic Pressure, Pdyn [kPa]
318
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
1
Fst = · γw · h2 = 0.5 · 9.81 · 542 = 14303kN (7.1)
2
It is shown that the maximum total hydrodynamic force for EQ5 was
larger than that for EQ2. This was expected, as EQ5 induced higher values
of acceleration in the dam structure and the reservoir domain and this is in
agreement to earlier theoretical studies (Westergaard, 1933; Zangar, 1952;
Chopra, 1968; Liu, 1986) which suggested a linear dependence of the max-
imum value of the hydrodynamic pressures on the maximum value of the
induced accelerations (see Section 2.5). Moreover, it may be observed that
the maximum value of the hydrodynamic pressure does not occur at the base
of the dam, but at some elevation. This is in agreement to the earlier obser-
vation of Zangar (1952) for dams with an inclined upstream face (see Section
2.5 and Figure 2.37).
Figures 7.46 and 7.48 show also a comparison of the calculated results
with analytical relations found in the literature. The grey line shows the
hydrodynamic pressures using the relations of Zangar (1952). The latter re-
lations were developed to estimate the peak hydrodynamic pressures on a
dam with a sloped upstream face from an incompressible reservoir (see rele-
vant discussion in Section 4.4). It is shown that the hydrodynamic pressures
calculated from this study were found to be larger than those obtained from
the relations of Zangar (1952). However, it should be noted that the study
of Zangar (1952) neglected the compressibility of the reservoir water (by as-
suming incompressible water). The significance of water compressibility was
highlighted by Chopra (1968) who showed that higher values of the hydro-
dynamic pressure may be expected in cases of resonance between the dam
and the reservoir.
319
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
• The results of Zangar (1952) were based on dams built on a rigid base,
whereas La Villita dam is built on a compliant soil layer foundation.
7.5.4 Comments
This section investigated the effects of reservoir-dam interaction on the dy-
namic response of a nonlinear earth dam built on a compliant base. This
is an extension of the work described in Chapter 5 which investigated the
effects of reservoir-dam interaction in visco-elastic dams founded on a rigid
base.
Minor differences in the dynamic response of the dam were found for
a dam with and without an upstream reservoir. Slightly smaller values of
acceleration were found for the reservoir-dam system considered and it was
concluded that the presence of the reservoir “damped” the response of the
dam. The hydrodynamic pressures calculated were found to be larger for
the larger seismic event (EQ5) and were found to be generally larger than
the values obtained using relations from the literature. This difference was
attributed to the simplifying assumptions made by the theoretical relations
found in the literature.
320
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
G(y) = linear
40 G(y) = constant
20
Elevation, z [m]
−20
−40
−60
−80
0 100000 200000 300000 400000 500000 600000 700000 800000 900000 1000000
Maximum Shear Modulus, Gmax [kPa]
Figure 7.49: Spatial variation (in the vertical direction) of maximum shear
stiffness, Gmax . The grey solid and black dashed lines correspond to a linear
and constant variation of Gmax in the dam core respectively.
321
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
a linear variation of maximum shear modulus with depth, Gmax (z). In order
to examine the effects of stiffness inhomogeneity in a non-linear analysis, an
additional analysis was carried out in which the maximum shear modulus of
the dam was taken as constant as shown in Figure 7.49. This value of Gmax
was chosen such that it is equal to the average value of the previously-used
linear variation, so that the fundamental period of vibration of the dam was
not significantly affected. It should be noted that the value of Gmax is also
used in the calibration of the CNL model (see Section 3.7.5).
7.6.3 Comments
This section investigated the effects of stiffness inhomogeneity in the dynamic
response of an earth dam. Earlier visco-elastic studies reported in the liter-
322
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.05
Acceleration, a [g]
−0.05
−0.1
−0.15
−0.2
0 1 2 3 4 5 6 7
Time, t [sec]
0.5
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
Figure 7.51: Response spectra (damping, ξ = 5%) the crest of La Villita dam
during EQ2.
323
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
0.2
Acceleration, a [g]
0.1
−0.1
−0.2
−0.3
0 10 20 30 40 50 60 70
Time, t [sec]
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Period, T [sec]
Figure 7.53: Response spectra (damping, ξ = 5%) the crest of La Villita dam
during EQ5.
324
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
20
Elevation, z [m]
−20
−40
−60
−80
1 2 3 4 5 6 7 8 9 10
Displacement, d [m] −3
x 10
20
Elevation, z [m]
−20
−40
−60
−80
0.13 0.135 0.14 0.145 0.15 0.155 0.16 0.165
Displacement, d [m]
325
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
20
Elevation, z [m]
−20
−40
−60
−80
0.04 0.06 0.08 0.1 0.12 0.14 0.16
Acceleration, a [g]
20
Elevation, z [m]
−20
−40
−60
−80
0.1 0.15 0.2 0.25 0.3 0.35
Acceleration, a [g]
326
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
7.7 Conclusions
This chapter dealt with various issues related to numerical analysis of the
seismic response of earth dams. The La Villita earth dam which was analysed
in Chapter 6, was used as an example for further investigation.
Several issues were considered related to the analysis approach followed
327
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
to analyse earth dams. Issues considered were (a) the effect of previous
earthquake history, (b) the effects of coupled and uncoupled undrained for-
mulation, (c) foundation-dam interaction, (d) reservoir-dam interaction and
(e) the effects of stiffness inhomogeneity. The main findings of this study
may be summarised as follows:
328
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
• Minor differences in the dynamic response of the dam were found for
a dam with and without an upstream reservoir (i.e. considering or
neglecting the reservoir hydrodynamic pressures). Slightly smaller val-
ues of acceleration were found for the reservoir-dam system considered
and it was concluded that the presence of the reservoir “damps” the
dynamic response of the dam.
329
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Chapter 8
CONCLUSION
8.1 Introduction
As described in the beginning of this thesis, in Chapter 1, the main aims of
this research were (a) to review the current state of the knowledge related to
the seismic response of earth dams and the available methods of analysis and
(b) investigate a real case study of seismic dam response in order to establish
and further improve the current knowledge of both seismic response and
analysis of earth dams.
The three objectives of the research, as those were described in the Intro-
duction of the thesis (Section 1.2), were presented in the previous chapters.
The first objective, which was to review the available literature and present
the main information found related to the response and analysis of earth
dams was presented in Chapters 2 and 3. The second objective, which was
to establish a methodology to model reservoir hydrodynamic pressures and
explore the effects of reservoir-dam interaction on the seismic behaviour of
dams, was presented in Chapters 4 and 5. Finally, the third objective, which
was to investigate the nonlinear behaviour of earth dams using advanced
elasto-plastic finite element analysis and assess the effects of different mod-
elling approaches, was presented in Chapters 6 and 7.
This chapter summarizes the findings of this research and recognising
its limitations, it describes the various recommendations for further inves-
330
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
tigation. At the end, some thoughts are discussed related to the relative
advantages and use of analytical and numerical methods of analysis.
• The size of the reservoir elements should be smaller than a fifth of the
acoustic (p-wave) wavelength, λ of the water.
331
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
332
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
• The calculated response spectrum for EQ5 at the crest of the dam was
compared to that of Elgamal (1992), who performed 3D shear beam
analysis, and was found to be in a better agreement with the recorded
response spectrum. Therefore, an improved prediction of the frequency
content of the acceleration was obtained in this work.
333
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
• The calculated vertical displacements at the crest of the dam were not
close to the recorded values. A localised slip failure that would result
into high values of settlements was not predicted and therefore no high
values of settlements were observed.
334
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
335
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
• Minor differences in the dynamic response of the dam were found for a
dam with and without an upstream reservoir. Slightly smaller values
of acceleration were found for the reservoir-dam system considered and
it was concluded that the presence of the reservoir damps the dynamic
response of the dam.
336
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
337
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
the early work on canyon effects by Mejia and Seed (1983) and Dakoulas
and Gazetas (1987) (see Section 2.4.2).
• The insignificant effects of RDI for earth dams were attributed (par-
tially) to the large volume of those dams, which makes the inertial
effects from the additional mass from the reservoir small compared to
the inertia of a large earth dam. Therefore, one possible direction for
further research would be to examine the amplification of accelerations
with respect to the inertial ratio (additional inertia due to the reservoir
hydrodynamic pressures over the inertia of the dam). The inertia due
to the reservoir hydrodynamic pressures is related to the magnitude of
those pressures due to the added mass effect. Therefore, examining RDI
effects with respect to the inertia ratio using numerical analyses would
be rather cumbersome, because the hydrodynamic pressures need to be
obtained for each case. It is consequently suggested that an analytical
(simplified and hence approximate) approach should be followed first
which would allow an easy analytical estimation of the hydrodynamic
pressures. It is the Author’s opinion and suggestion that the work of
Lee and Tsai (1991) (see Section 2.5.6) could be the starting point.
338
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
339
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
• Although not all (and not the whole duration of) the earthquake time
histories are available, a series of successive earthquake events should
be analysed in order to examine the effect even of the smaller previous
earthquakes. Those earthquakes may perhaps have a contribution in
the accumulation of permanent displacements of the dam.
340
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
341
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
was employed in this work, were identified as the most rigorous and reliable
methods of analysis.
342
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
343
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
344
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
References
345
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
346
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Bielak, J., Loukakis, K., Hisada, Y. and Yoshimura, C. (2003), ‘Domain Re-
duction Method for Three-Dimensional Earthquake Modelling in Localized
Regions. Part I: Theory’, Bulletin of the Seismological Society of America
93(2), 817–824.
347
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
348
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
349
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
350
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Dakoulas, P. and Gazetas, G. (2008), ‘Insight into seismic earth and water
pressures against caisson quay walls’, Géotechnique 58(2), 95–111.
Day, R. A. (1990), Finite element analysis of sheet pile retaining walls, PhD
thesis, Imperial College, University of London.
351
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Elgamal, A. W., Scott, R. F., Succarieh, M. F. and Yan, L. (1990), ‘La Villita
dam response during five earthquakes including permanent deformation’,
Journal of Geotechnical Engineering 116(10), 1443–1462.
Elia, G., Amorosi, A., Chan, C. and Kavvadas, M. J. (2010), ‘Fully coupled
dynamic analysis of an earth dam’, Géotechnique 61(7), 549–563.
352
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
353
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
354
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
355
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Idriss, I. M., J., L., Hwang, R. and Seed, H. B. (1973), QUAD-4: A computer
program for evaluating the seismic response of soil structures by variable
damping finite element procedures, Technical Report Report EERC 73-16,
Earthquake Engineering Research Center, University of California, Berke-
ley.
356
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Kokusho, T. (1980), ‘Cyclic triaxial test of dynamic soil properties for wide
strain range’, Soils and Foundations 20(2), 45–60.
Kontoe, S., Zdravković, L. and Potts, D. M. (2009), ‘An assessment of the do-
main reduction method as an advanced boundary condition and some pit-
falls in the use of conventional absorbing boundaries’, International Jour-
nal for Numerical and Analytical Methods in Geomechanics 33(3), 309–
330.
Kontoe, S., Zdravković, L., Potts, D. M. and Menkiti, C. (2008), ‘Case study
on seismic tunnel response’, Canadian Geotechnical Journal 45(12), 1743–
1764.
Kontoe, S., Zdravković, L., Potts, D. M. and Menkiti, C. (2011), ‘On the
relative merits of simple and advanced constitutive models in dynamic
analysis of tunnels’, Géotechnique 61(10), 815–829.
357
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Lysmer, J., Udaka, T., Tsai, C. F. and Seed, H. B. (1975), FLUSH: A com-
puter program for approximate 3-D analysis of soil-structure interaction
358
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
McCully, P. (1996), Silenced Rivers: The Ecology and Politics of Large Dams,
Zed Books, London.
359
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
360
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Poulos, H. G. and Davis, E. H. (1974), Elastic solutions for Soil and Rock
Mechanics, John Wiley and Sons, New York.
361
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Rollins, K., Evans, M., Diehl, N. and Daily, W. (1998), ‘Shear modulus
and damping relations for gravel.’, J. Geotechnical and Geoenvironmental
Engineering, ASCE 124(5), 396–405.
362
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Seed, H. B. and Idriss, I. M. (1970), Soil moduli and damping factors for
dynamic response analyses, Technical Report EERC 70–10, Earthquake
Engineering Research Centre, University of California, Berkeley.
Seed, H. B., Lee, K. L. and Idriss, I. M. (1969), ‘Analysis of Sheffield dam fail-
ure’, Journal of the Soil Mechanics and Foundation Engineering Division,
ASCE 95(SM6), 1453–1490.
Seed, H. B., Makdisi, F. I., Idriss, I. M. and Lee, K. L. (1975), ‘The Slides
in the San Fernando Dams During the Earthquake of February 9, 1971’,
Journal of the Geotechnical Engineering Division 101(7), 651–688.
363
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Singh, R., Roy, D. and Jain, S. K. (2005), ‘Analysis of earth dams affected
by the 2001 bhuj earthquake’, Engineering Geology 80(3-4), 282–291.
Tani, S. (2000), ‘Behavior of large fill dams during earthquake and earthquake
damage’, Soil Dynamics and Earthquake Engineering 20(1-4), 223–229.
364
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
365
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
U.S. Army Corps of Engineers (1994), Engineering and design - arch dam
design, Technical Report EM 1110-2-2201, Department of the Army.
U.S. Army Corps of Engineers (1995), Engineering and design - gravity dam
design, Technical Report EM 1110-2-2200, Department of the Army.
Von Estorff, O. and Antes, H. (1991), ‘On FEM-BEM coupling for fluid-
structure interaction analyses in the time domain’, International Journal
for Numerical Methods in Engineering 31(6), 1151–1168.
366
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
Wilson, E. L. (1975), Finite element for foundation, joints and fluid, in ‘Con-
ference on Numerical Methods in Soil and Rock mechanics’, Karlsruhe.
Wilson, E. L. and Khalvati, M. (1983), ‘Finite elements for the dynamic anal-
ysis of fluid-solid systems’, International Journal for Numerical Methods
in Engineering 19(11), 1657–1668.
367
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
368
SEISMIC RESPONSE AND ANALYSIS OF EARTH DAMS
369