Diss Benz
Diss Benz
Thomas Benz
Thomas Benz
IGS
vorgelegt von
T HOMAS B ENZ
aus Tubingen
Hauptberichter:
Mitberichter:
Mitteilung 55
Geotechnik
des Instituts fur
Universitat Stuttgart, Germany, 2007
Editor:
Prof. Dr.-Ing. P. A. Vermeer
c
Thomas
Benz
Geotechnik
Institut fur
Universitat Stuttgart
Pfaffenwaldring 35
70569 Stuttgart
All rights reserved. No part of this publication may be reproduced, stored in a
retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, scanning or otherwise, without the permission in
writing of the author.
Keywords: Small-strain stiffness, constitutive soil models
Printed by e.kurz + co, Stuttgart, Germany, 2007
ISBN 978-3-921837-55-9
(D93 - Dissertation, Universitat Stuttgart)
Preface
In the seventies when I began to do research it was realised that soil behaviour was
highly non-linear, but at that time we still underestimated the extent of it. At that
time I should have had the book by Tsytovich (1973)1 , who taught soil mechanics at
the Moscow Civil Engineering Institute. In 1976 this book was translated into English
and I read it some ten years later. In this book Professor Tsytovich introduces a structural soil strength for stress increments as induced by external loads. Up to its structural
strength soil behaves extremely stiff as there is no rearrangement of particles. Once stress
increments exceed this threshold value, soils show the nowadays well-known stressdependent stiffness, as for instance expressed by hyperbolic and/or exponential rules.
In practical analysis of foundation settlements the high initial small-strain stiffness has
always been taken into account by the introduction of a so-called limit depth; below this
depths strains were simply assumed to be negligibly small. Tsytovich does not use the
name limit depth, but calls it depth of the active compression zone. In finite element
calculations I always accounted for this limit depth by using relatively shallow meshes.
From a mathematical point of view the subsoil should obviously be modelled by infinite
boundaries, but because of the small-strain stiffness only relatively shallow finite layer
had to be discretised. On using the new constitutive model by Dr. Thomas Benz this
dilemma on the depth of finite element meshes is non-existent. Boundaries can now
be chosen so for away that they do not influence the computed stress distribution and
the constitutive model ensures that the deep part of the mesh behaves virtually incompressible. The new HS-Small model is obviously not only relevant for the analysis of
foundations, but also for the analysis of settlements due to tunnelling or deep excavations, as shown in this dissertation study.
The focus of this study is on the small-strain stiffness of non-cemented soils, but the
model may also be applied to somewhat bonded soils. In recent years the behaviour of
sensitive and/or structured clays has got much attention from researches and it would
seem that their findings can also be modelled by the new HS-small model. In order to
do so one may have to choose relatively large values for the strain range parameter 0.7 ,
as used in the HS-Small model.
This thesis on a new computational model for the behaviour of soils at small strains
is a scholarly piece of work, covering a very wide scope of material. The new model
is developed, and calibrated, and in my opinion extremely useful for engineering problems. I am also indebted to the Federal Waterways Engineering and Research Institute
in Karlsruhe who have sponsored this study. Within this institute Dr. Radu Schwab has
advised Thomas Benz and it has been a real pleasure for me to work with them.
Pieter A. Vermeer
Stuttgart, March 2007
N. Tsytovich. Soil Mechanics. Mir Publishers, Moscow, 1976. (translation from the 1973 Russian edition)
Acknowledgments
This Dissertation would not have been completed without the support and the patience
of Professor Pieter Vermeer and Radu Schwab. I thank Pieter for giving me the opportunity to join his research team and for his extensive help. Likewise, I thank Radu for the
unlimited help, advice and support he provided at the Federal Waterways Engineering
and Research Institute (BAW). This Dissertation emanated from a research project that
was initiated and funded by the BAW.
I thank all colleagues at the BAW in Karlsruhe and at the Institute of Geotechnical
Engineering in Stuttgart for their friendly help and for the pleasant years. In particular,
I am grateful to the following persons (in random order) for the exchange of thoughts,
critical reviews, beta testing, proof reading, etc: Markus Herten, Marcin Cudny, Markus
Contents
Introduction
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
9
10
12
12
12
13
14
14
15
15
15
15
16
16
19
21
23
23
24
25
28
30
31
31
34
34
35
35
38
Contents
4
ii
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
their use in
. . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
39
39
41
42
44
46
47
49
49
50
51
53
55
56
58
.
.
.
.
.
.
.
.
.
.
.
59
59
59
62
63
65
66
66
66
69
72
72
.
.
.
.
.
.
.
.
.
.
.
.
75
75
77
82
83
85
86
91
94
97
98
100
102
Contents
8
Lock
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
105
105
115
118
118
120
121
123
126
130
130
134
136
139
10 Conclusions
147
Bibliography
150
167
173
177
D Material data
179
iii
Abstract
This thesis is concerned with the stiffness of soils at small strains and constitutive models that can be applied to simulate this. The strain range in which soils can be considered
truly elastic, i.e. where they recover from applied straining almost completely, is very
small. With increasing strain amplitude, soil stiffness decays non-linearly: Plotting soil
stiffness against log(strain) yields characteristic S-shaped stiffness reduction curves. At
the minimum strain which can be reliably measured in classical laboratory tests, i.e. triaxial tests and oedometer tests without special instrumentation, soil stiffness is often
less than half its initial value. If this non-linear variation of soil stiffness at small strains
is considered in the analysis of soil-structure interaction, analysis results improve considerably: The width and shape of settlement troughs are more accurately modeled,
excavation heave is reduced to a more realistic value, etc. The importance of smallstrain stiffness in engineering practice is therefore well recognized. As the few existing
small-strain stiffness models are either research orientated, and/or applicable to specific
loading paths only, the objective of this thesis is to develop a small-strain stiffness model
for the engineering community.
The newly developed Small-Strain Overlay model is based on the Hardin-Drnevich
model [52]. Similar to the Hardin-Drnevich model, the Small-Strain Overlay model has
two input parameters. Both input parameters have a clear physical meaning. For strains
below the limit of classical laboratory testing, the Small-Strain Overlay model calculates
the actual strain dependent stiffness. This stiffness can subsequently be deployed in
many existing elastoplastic constitutive models. One of them is the well known Hardening Soil (HS) model as implemented in the finite element code P LAXIS V8 [21]. In this
thesis, the Small-Strain Overlay model is combined with the HS model. The resulting
HS-Small model is additionally enhanced by a Matsuoka-Nakai failure criterion. The
HS-Small model is validated in element tests and various 2D and 3D boundary value
problems. In the validation sequence, all effects that are commonly attributed to smallstrain stiffness in soil-structure interaction can be recognized.
Although the thesis main concern is with the development of a small-strain stiffness
model, it also provides much information and data that may be helpful in using it. This
information and data comprise an introduction to available small-strain stiffness testing
methods, many empirical small-strain stiffness correlations and sample model applications.
After a short introduction to the topic in Chapter 1 and clarification of important terminology in Chapter 2, the thesis is organized thematically in three parts:
Abstract
Tests and correlations - Chapter 3 Firstly, small-strain stiffness is discussed from
an experimental point of view. The most important small-strain stiffness testing
methods are introduced. Different parameters that affect small-strain stiffness are
identified. The influence of the most significant of these parameters, regarding
small-strain stiffness, is quantified in a collection of empirical correlations. In the
absence of experimental data, these correlations may prove very helpful in using
the new model.
Model building and formulation - Chapter 4 to 7 Secondly, the new Small-Strain
Overlay model and its combination with the HS model is formulated. This main
part of the thesis also contains a critical review of existing small-strain models.
Important numerical aspects, for example the local integration of the incrementally
formulated HS-Small equations, are presented at the end of the model building and
formulation section.
Code verification and model validation - Chapter 8 to 9 Finally, the HS-Small
model is validated in element tests, several 2D boundary value problems, and a
state-of-the-art 3D boundary value problem. Initialization of the HS-Small model
and a new possibility for verifying the proper extent of boundary conditions are
discussed within the validation section as well.
Results from the study presented in this thesis are concluded in Chapter 10.
vi
Zusammenfassung
Die vorliegende Arbeit beschaftigt sich mit der Steifigkeit von Geomaterialien unter
kleinen Dehnungen und deren Beschreibung in Stoffmodellen. Der Dehnungsbereich, in
Entlastung der ursprungliche Dehnungszustand nahezu wieder einstellt - ist sehr klein.
Auerhalb dieses sehr kleinen Bereiches zeigen Versuche eine nicht lineare Steifigkeits
abnahme mit zunehmender Dehnung. Aufgetragen uber
den Logarithmus der aufge die Abnahmefunktion der Steifigkeit eine charakterbrachten Dehnung, ergibt sich fur
istische S-formige
Kurve. Bereits bei der kleinsten Dehnung, die noch zuverlassig in
klassischen Laborversuchen (d.h. Dreiaxiale Versuche und Oedometer Versuche ohne
spezielle Instrumentierung) gemessen werden kann, hat sich die Steifigkeit, bezogen
auf ihren Ausgangswert, oftmals schon um mehr als die Halfte reduziert. Bei Berucksichtigung dieses nicht-linearen Verhaltens in Boden-Bauwerk Interaktionsproblemen
sich die Genauigkeit von berechneten Setzungen und Verschiebungen erheblich:
erhoht
Die laterale Ausdehnung und Form von Setzungsmulden wird erheblich besser abgebildet, Hebungen in Baugrubensohlen realistischer berechnet, etc. Die Bedeutung von
nicht-linearem Bodenverhalten bei kleinen Dehnungen wird daher auch in der praktischen geotechnischen Entwurfs- und Bemessungsarbeit allgemein als hoch eingeschatzt.
vii
Zusammenfassung
Obwohl das primare Ziel der vorliegenden Arbeit die Entwicklung eines Stoffmodells
ist, enthalt diese auch zahlreiche Informationen und Daten, die sich bei der Nutzung
Moglichkeiten
der Modellinitialisierung diskutiert sowie eine Bewertungsmethode
die Randabstande von Berechnungsausschnitten vorgeschlagen.
fur
Die Schlufolgerungen aus der vorliegenden Arbeit sind in Kapitel 10 zusammengefat.
viii
Chapter 1
Introduction
One of the major problems in ground engineering in the 1970s and earlier was the
apparent difference between the stiffness of soils measured in laboratory tests and those
back-calculated from observations of ground movements (e.g. Cole & Burland [29], St
John [84], Wroth [195], Burland [23]). These differences have now largely been reconciled
through the understanding of the principal features of soil stiffness and, in particular,
the very important influence of non-linearity. This is one of the major achievements of
geotechnical engineering research over the past 30 years. (Atkinson [7])
In particular the non-linear influence of strain on soil stiffness has been extensively investigated over the past decades. The maximum strain at which soils exhibit almost fully
recoverable behavior is found to be very small. The very small-strain stiffness associated
with this strain range, i.e. shear strains s 1 106 , is believed to be a fundamental
property of all types of geotechnical materials including clays, silts, sands, gravels, and
rocks (Tatsuoka et al., 2001) under static and dynamic loading (Burland [24]) and for
drained and undrained loading conditions (Lo Presti et al. [140]). With increasing strain,
soil stiffness decays non-linearly. On a logarithmic scale, stiffness reduction curves exhibit a characteristic S-shape, see Figure 1.1.
The smallest shear strain that can be reliably measured in conventional soil testing,
e.g. triaxial or oedometer tests without special instrumentation, is s 1 103 . By
definition Atkinson [7] terms strains smaller than the limit of classical laboratory testing
(s < 1 103 ), small strains. Strains, s > 1 103 are termed large or larger strains. The
limit of classical laboratory testing coincides at the same time with characteristic shear
strains that can be measured near geotechnical structures (Figure 1.1). However, the soil
stiffness that should be used in the analysis of geotechnical structures is not the one that
relates to these final strains. Instead, very small-strain soil stiffness and its non-linear
dependency on strain amplitude should be properly taken into account in all analysis
that strive for reliable predictions of displacements. The Rankine lectures by Simpson
[168] and Atkinson [7], or the Bjerrum Memorial lecture by Burland [24] are just a few
occasions where this has been highlighted. As yet, small-strain stiffness has not been
widely implemented in engineering practice. Considering numerical analysis, this may
be due to a lack of capable, yet user-friendly constitutive models. The main objective
of this thesis is to provide such a capable and sufficiently simple small-strain stiffness
model for engineering practice.
The user-friendliness of a constitutive model largely depends on the input parameters. They should be limited in their number, easy to understand in their physical meaning, and easy to quantify based on test data or experience. The level of sophistication
Chapter 1 Introduction
Retaining walls
Foundations
Tunnels
Very
small
strains
Small strains
Larger strains
0
-6
10
10
-5
10
-4
10
-3
10
-2
10
-1
Dynamic methods
Local gauges
Figure 1.1: Characteristic stiffness-strain behavior of soil with typical strain ranges for
laboratory tests and structures (after Atkinson & Sallfors [9] and Mair [109])
of a model is tied to the extent to which it is able to reproduce the experimentally observed functional relationships. The premise of fulfilling the objective defined above is
therefore a thorough description of experimental observations. It is only through this
that the models parameters and mechanisms can be decided upon, or in the words of
Ernst Mach: The description of functional relationships is an explanation in itself. The
methodology applied in this thesis is mostly inductive reasoning. Figure 1.2 explains the
terminology used to describe the different stages within the model building and validation process. The model is formulated incrementally. Therefore the code verification and
model validation ar regarded as separate processes.
Although the main objective of this thesis is the development of a simple and capable small-strain stiffness constitutive model, it also acts as a compilation of information
and data that might be helpful in using it. This information and data comprise empirical correlations that can be useful in quantifying small-strain stiffness, an introduction
to available testing methods and sample applications in 2- and 3-dimensional finite element analysis.
The outline of this thesis is as follows:
Chapter 2 introduces some definitions and conventions used throughout the thesis.
Chapter 3 concentrates on the experimental aspects and the quantitative description of small-strain stiffness. After introducing the experimental concepts of laboratory and in-situ testing, the influence of various parameters on small-strain stiffness is studied. Readily available test data and correlations from literature are
presented at the end of this chapter.
Experiment
Model
Validation
Simulation
Numerical Model
Abstraction
Programming
Model
Building
Model
Code Verification
Figure 1.2: Terminology used in the model building and validation process.
Chapter 1 Introduction
different problem categories is quantified. More general issues, as for example
possible initialization procedures for the HS-Small model, are discussed at the end
of this chapter.
Sulfeld
is currently being reconstructed right next to the high-speed railway link
between Hannover and Berlin. A state of the art 3D finite element analysis of the
deep excavation is conducted in order to predict the construction activitys influence on the railway tracks. It is shown how the HS-Small model could significantly
enhance the reliability of this analysis.
Finally, Chapter 10 presents the main conclusions of this study.
Chapter 2
Terminology and definitions
This thesis is mainly concerned with pre-failure deformation characteristics of soils.
Continuum mechanics is generally employed as the analysis tool. In soil mechanics
however, there is not a clear consensus in how to define some stress and strain properties in this framework (e.g. sign convention). This chapter briefly introduces the stress
and strain definitions, notations, and invariants used in the following chapters.
Tensor notation Tensorial quantities are generally expressed in indicial notation.
The order of a tensor is indicated by the number of unrepeated (free) subscripts.
Whenever a subscript appears exactly twice in a product, that subscript will take
on the values 1, 2, 3 successively, and the resulting terms are summed (Einsteins
summation convention). The Kronecker delta ij takes the value 1 or 0 if i = j or
i 6= j respectively. In the rare occasions where indicial notation is not used in this
thesis, first-order tensors (vectors) are denoted by lowercase Latin or Greek bold
letters, second-order tensors are denoted by uppercase Latin or Greek bold letters,
and fourth-order tensors are denoted by Calligraphic letters.
Stress and strain Infinitesimal deformation theory is applied. Cauchy stress is
related to linearized infinitesimal strain. Eigenvalues of stress and strain tensors
(principal stresses and strains) are denoted by one subscript only, e.g. i and i
with i = 1, 2, 3. Without loss of generality, let
1 2 3 .
(2.1)
The Roscoe stress invariants p (mean stress) and q (deviatoric stress), are defined
as:
r
ii
3
1
1
p=
and
q=
(ij ij kk )(ij ij kk ),
(2.2)
3
2
3
3
In triaxial compression with 1 2 = 3 , the Roscoe invariants simplify to
1
(axial + 2lateral )
3
q = (axial lateral ).
p =
(2.3)
(2.4)
In analogy to the stress invariants, volumetric strain v and shear strain s are defined as:
r
1
1
3
v = ii
and
s =
(ij ij kk )(ij ij kk ),
(2.5)
2
3
3
ETan(A)
E0
E0
ESec(A)
1
E0
eaxial
Figure 2.1: Definition of secant and tangent moduli in triaxial stress-strain space.
which simplify to
v = axial + 2lateral
s = (axial lateral ),
(2.6)
(2.7)
laboratory testing, which is according to Atkinson [7] s 1 103 (see Figure 1.1).
In this thesis however, the limit of classical laboratory testing is taken as the shear
strain where the stiffness modulus reduction curve equals the unloading-reloading
stiffness observed in classical laboratory testing.
Small-strain stiffness Small-strain stiffness is the stiffness of soils at small shear
strains. The maximum soil stiffness found in stiffness reduction curves is in this
thesis referred to as initial, or maximum soil stiffness. Initial shear and Youngs
moduli are denoted as G0 and E0 respectively. Often, initial soil stiffness is also
called very small-strain stiffness, as it is commonly associated to the very-smallstrain range as defined above. In this thesis, small-strain stiffness is always discussed along with the respective stiffness reduction curve in the range of very
small and small strains. The term small-strain stiffness model consequently refers
in this thesis to a constitutive soil model that is applicable to very small and small
strains; Small-strain stiffness measurements include stiffness measurements at very
small and small strains; etc.
Chapter 3
Experimental evidence for small-strain stiffness
From dynamic response analysis, it has been found that most soils have curvilinear
stress-strain relationships as shown in Figure 3.1. The shear modulus is usually expressed as the secant modulus determined by the extreme points on the hysteresis loop
while the damping factor is proportional to the area inside the hysteresis loop. It is
readily apparent, that each of this properties will depend on the magnitude of the strain
for which the hysteresis loop is determined and thus, both shear moduli and damping
factors must be determined as functions of the induced strain in a soil specimen or soil
deposit. (Seed & Idriss [162])
In soil dynamics, small-strain stiffness has been a well known phenomena for a long
time. The above conclusions by Seed & Idris for example were drawn more than 35
years ago. In static analysis however, the findings from soil dynamics have long been
disregarded, or even considered not to be applicable. Seemingly differences between
static and dynamic soil stiffness have been attributed to the nature of loading (e.g. inertia
forces and strain rate effects) rather than to the magnitude of applied strain, which is
generally small in dynamic conditions (earthquakes excluded).
Nowadays, it is commonly accepted that inertia forces and strain rate have little influStress
G1
e1
G2
1
e2
Strain
Figure 3.1: Stress-strain hysteresis loops presented by Seed & Idris [162]. The secant
shear modulus of a closed stress-strain loop decays monotonically with its
strain amplitude: G2 < G1 for 2 > 1 .
10
+ 2G
and
(3.1)
(3.2)
Travel direction
Travel direction
Figure 3.2: P-wave (top) and S-wave (bottom) particle motion: The particle motion in Pwaves is longitudinal whereas it is transverse in S-waves. The particle motion
vector in the plane perpendicular to the direction of propagation, is referred
to as polarization (only S-waves).
where vp is the propagation velocity of pressure, or primary (P-) waves, vs is the propagation velocity of shear or secondary (S-) waves, and , and G are Lames constants (G
is also termed the shear modulus). The expression of Lames constants as a function of
Youngs modulus and Poisons ratio is generally more convenient for engineers:
E
2(1 + )
E
=
.
(1 + )(1 2)
G =
(3.3)
(3.4)
11
12
horizontal
benders
axial
benders
25mm
Figure 3.3: Bender elements mounted in a top plate and to a radial belt. Photograph by
Pennington et al. [136]
Originally, bender elements could transmit and receive S-waves only. Providing a
samples density is known, the knowledge of shear wave velocity is sufficient to calculate its shear modulus G0 from Equation 3.2. Unfortunately, S-waves are slower than
P-waves (vs 12 vp ) so that they always have noisy first arrivals (in this context, reflected
P-waves and surface waves are considered noise). Less valuable is the isolated knowledge of P-wave velocity as this is governed by both elastic constants. Ideal is the knowledge of both, S and P wave velocities, which can be simultaneously measured using
newly developed bender-extender elements (Lings & Greening [104]).
Compared to local strain transducers, the obvious disadvantage of bender, and bender-extender elements is their restriction to the very small-strain range, e.g. s < 1
106 . The geophysical nature of the indirect bender element stiffness measurements can
sometimes be troublesome too: Due to the low signal to noise ratio, the time at which
the first wave arrives is usually subject to interpretation. Sample preparation on the
other hand is simple. Commercially offered bender element mounts for top cap and
base pedestal along with user-friendly software promote bender elements as a feasible
option for routine lab testing. Figure 3.3 shows the integration of bender elements into a
top cap and into a radial belt.
3.1.1.3 Resonant column and torsional shear
Resonant column and torsional shear devices can load soil samples not only triaxially
but also torsionally. The basic difference between resonant column and torsional shear
testing is the frequency and amplitude of loading.
Torsional shear tests are static, or quasi-static cyclic tests where an axially confined
13
14
15
16
Figure 3.4: Small-strain stiffness decay as a function of (a) effective friction angle 0 , (b)
vertical effective stress v0 , (c) void ratio e, and (d) K0 [162].
17
Table 3.1: Parameters affecting the stiffness of soils at small-strains (modified after
Hardin & Drnevich [53]).
Parameter
Importance to a
G0
Clean
Cohesive
sands
soils
V
V
V
V
V
V
V
R
L
V
V
R
R
R
R
L
L
L
L
0.7
Clean
Cohesive
sands
soils
V
V
V
V
R
V
V
R
L
R
R
V
V
R
R
L
L
R
R
Strain amplitude
Confining stress
Void ratio
Plasticity index (PI)
Overconsolidation ratio
Diagenesis
Strain history
Strain rate
Effective material strength
Grain Characteristics
(size,shape,gradation)
Degree of saturation
R
V
L
L
Dilatancy
R
R
R
R
a
V means Very Important, L means Less Important, and R means
Relatively Unimportant
18
1.4
1.2
0.8
0.4
0.0
400
800
1200
1600
Figure 3.5: Isotropic compression test interrupted by small-strain cycles (after Lade &
Abelev [98]).
is much less than in primary loading: Figure 3.5 shows a recent test result by Lade &
Abelev [98]. In their test, Lade & Abelev compared the stiffness in isotropic loading
and unloading to that obtained when interrupting the continuous loading process by
small load cycles. In primary loading, soil stiffness within the load cycles is found to be
higher than the soil stiffness in continuous loading. Load cycling during unloading on
the other hand did not show any significant stiffness increase. This is also in agreement
with the findings of Zdravkovic & Jardine [202]: The secant and tangent bulk moduli
curves developed during swelling fall relatively gently with strain, remaining far above
the tangent compression value.
The remaining part of this chapter is hence focused on the reduction of soil stiffness
with shear strain. The terms strain and shear strain are sometimes used as synonyms.
Small-strain stiffness in volumetric loading is again discussed in the next chapter.
(3.5)
which has not yet been superseded. Hardin & Richard themselves used the power law
exponent m = 0.5 for both, cohesive and non-cohesive soils. Today, their exponent is
widely confirmed for non-cohesive soils: All recent correlations use exponents in the
range of 0.40 m 0.55. For cohesive soils, the exponent m = 0.5 is controversial.
Many researchers confirmed it, others found exponents as high as m = 1.0.
At this point, it has to be considered that confining stress influences the degradation
of small-strain stiffness as well: Damping decreases with increasing confining stress.
19
1.0
1.0
0.9
0.9
Exponent m
Exponent m
0.8
0.7
0.6
0.5
0.8
0.7
0.6
10
20
30
40
50
60
0.5
20
40
60
80
100
120
140
160
180
Figure 3.6: The power law exponent m as a function of plasticity index (PI) and liquid
limit wL (after Viggiani & Attkinson [187], and Hicher [57]).
Correlating secant stiffnesses at greater strain amplitudes (s > 1 106 ) yields therefore higher m values than strictly using low strain bender element measurements (e.g.
s < 1 106 ). Additionally, it has to be taken into account that Hardin & Richard, like
many others, use a void ratio term in their relationship (see Equation 3.8), which also
relates to confining stress. Considering that non-cohesive soils are typically less compressible than cohesive soils, the scatter in the exponent m for cohesive soils can readily
be explained. Without taking into account void ratio in their relationship, Viggiani &
Attkinson [187] compiled the exponents m for different clays at very small strains as a
function of plasticity index; Hicher [57] compiled them as a function of Liquid Limit.
Both charts are shown in Figure 3.6. Small-strain stiffness data for a Kaolin clay with
PI 20 and LL 43 acquired by Rammah et al. [146] even correlated well, with m = 1.0.
This suggests that whenever void ratio is assumed constant in the corresponding relationship, the small-strain power law exponents m are equal (sands) or only a little less
(clays) than the exponents generally obtained in large strain oedometer and triaxial tests,
e.g. m = 0.5 for sands and m = 0.7 1.0 for clays.
Normalized modulus reduction curves, as for example the ones shown in Figure 3.8,
indicate that the threshold shear strain 0.7 is also mean stress dependent. Only a few relationships between confining stress and reduction of small-strain stiffness are proposed
in literature though. That by Ishibashi & Zhang [70] varies the exponent m in Equation
3.5 non-linearly with strain, as proposed earlier by Iwasaki et al. [76], [75]. A general recommendation for practical application can hardly be given though from their analysis.
Darendeli & Stokoe [30] developed normalized modulus reduction curves from an extensive study on over 100 undisturbed specimens from depths of 3 to 263 m. From their
data, the threshold shear strain 0.7 correlates very well to confining pressure p0 where:
0 m
p
,
(3.6)
0.7 = (0.7 )ref
pref
and pref = 100 kPa is a reference pressure, (0.7 )ref is the threshold shear strain at
20
2.5
m=
0.6
5
m=
0.5
0
Iwasaki et al.
m
=0
.3
5
Kallioglou et al.
2
1.5
1
g0.7(p'):=(g0.7)ref*(p'/pref)
0.5
0
10
100
1000
Figure 3.7: Correlation between confining stress p0 and shear strain 0.7 .
p0 = pref , and m
= 0.35. In Figure 3.7, the result from Darendeli & Stokoe is presented
together with other test data on non-plastic soils, i.e. data by Iwasaki et al. [76] that
was later used by Ishibashi & Zhang [70] for their correlation. From this representation
it can be concluded that a) the power law works reasonably well also for the threshold
shear strain and b) the power law exponent for non-plastic soils under moderate confining pressures is typically 0.35 < m
< 0.65. According to Stokoe et al. [89], Equation 3.6
can be used in combination with m
= 0.35 in plastic soils, as well. Other test data does
not always support this statement. Biarez & Hicher [16], for example, find the threshold
shear strain 0.7 of Kaolinite (PI = 30) largely unaffected by confining stresses in the range
of p0 = 100 300 kPa.
v = a(b e)p0 2 .
(3.7)
From this linear dependency, Hardin & Richart derived their well known formula for
void ratio dependency of G0 to:
(2.17 e)2
1+e
(2.97 e)2
1+e
G0
G0
and
(3.8)
(3.9)
21
1.1
1.0
1.0
0.9
G / G0 [-]
G / G0 [-]
0.9
0.8
ID=0.95
ID=0.85
0.7
ID=0.70
ID=0.63
0.6
ID=0.50
0.5 all tests:
p'=80 kPa
0.4
10 -6
10 -5
0.8
0.7
0.6
10 -4
10 -3
0.5
10 -6
p'=400 kPa
p'=200 kPa
p'=100 kPa
p'= 50 kPa
all tests:
ID=0.63-0.66
10 -5
10 -4
10 -3
Figure 3.8: Influence of void ratio e (left) and confining stress p0 (right) on the decay of
small-strain stiffness (after Wichtmann & Triantafyllidis [193]). The influence
eemin
of void ratio is expressed by the density index ID = emax
.
emin
for round-grained sands (e < 0.80), and angular-grained sands (e > 0.60) respectively.
Hardin & Black [50], [51] later indicated that Equation 3.9 also correlates reasonably
well for clays with low surface activity. For clays with higher surface activity, the basic
structure of relationship 3.9 is often maintained, only the coefficient 2.97 is replaced by a
somewhat increased one (see Table 3.3).
Other relationships between void ratio and initial stiffness found in literature are typically of the form:
G0 ex
(3.10)
where the exponent x is quantified for instance as:
x = 0.8 (for sand - Fioravante [40])
x = 1.0 (for sand and clay - Biarez & Hicher [16])
x = 1.3 (for cemented sands, fine-grained soils - Lo Presti [141])
1.1 x 1.5 (for various clays - Lo Presti & Jamiolkowski [139]).
The influence of void ratio on the threshold shear strain 0.7 is apparently very limited in non-cohesive soils. This can be observed for example in normalizing the results
of Figure 3.4.c or in the more recent results presented by Wichtmann & Triantafyllidis
[193] (Figure 3.8). Damping in cohesive soils on the other hand is linked to void ratio.
In cohesive soils, void ratio also correlates to plasticity index (PI), as higher plasticity is
generally a prerequisite for a more open soil structure and thus a higher void ratio. Yet,
void ratio is also related to the confining stress and overconsolidation ratio. Therefore,
the influence of void ratio on the threshold shear strain 0.7 is not discussed in this section. Instead the effects of PI, OCR, and p0 on 0.7 are discussed in the following sections.
22
1.1
1.0
1.0
0.9
0.9
G / G0 [-]
1.1
0.8
PI
=1
5
0.8
PI
=5
PI
=1
00
0.7
=3
0
0.6
=0
PI
PI=0
0.7
PI=1-7
PI=8-20
0.6
PI=21-40
PI=41-44
0.5 47 soils tesed:
sV'=24-1700 kPa
0.4
10 -6
10 -5
PI=200
PI
G / G0 [-]
0.5
10 -4
10 -3
0.4
10 -6
OCR=1-15
10 -5
10 -4
10 -3
Figure 3.9: Influence of plasticity index (PI) on stiffness reduction: Left database for soils
with different PI; Right: PI-chart by Vucetic & Dobry (after Hsu & Vucetic
[64], and Vucetic & Dobry [188] respectively).
(3.11)
where they define OCR as the ratio of maximum past vertical effective stress to the cur0
(3.12)
23
(3.13)
Again the empirical parameter k increases with clay plasticity. For clays with 10 < PI <
40, Atkinson & Viggiani found 0.20 < k < 0.25. For practical applications this seems
a relatively small variation, such that it is sometimes proposed to neglect the effect of
OCR on G0 completely, e.g. in Lo Presti & Jamiolkowski [139]. However, overconsolidation cannot be neglected in the inherent and stress-induced anisotropy of small-strain
stiffness, which is discussed later in Section 3.2.7.
In non-plastic soils, overconsolidation has only a small effect on the normalized modulus reduction curve. In plastic soils, overconsolidation markablly increases the threshold
shear strain 0.7 . The increase of 0.7 with overconsolidation OCR can after Stokoe et al.
[89] be approximated as:
0.7 = (0.7 )ref + 5 106 PI(OCR)0.3 ,
(3.14)
where (0.7 )ref is here a reference threshold shear strain for a non-overconsolidated, nonplastic soil, e.g. (0.7 )P I=0,OCR=1 1e4. For overconsolidated soils, the findings by
Stokoe et al. [89] and Vucetic & Dobry [188] are thus also in better agreement for soils
with higher plasticity (Note that the PI-chart by Vucetic & Dobry is independent of overconsolidation).
t
G0 (t)
= 1 + NG,1 log
(3.15)
G0 (tp )
tp
24
(3.16)
where tp is the time required to reach the end of primary consolidation, t is any time
(t > tp ), G0 (tp ) is the maximum stiffness at time tp , and NG,x is an empirical material
factor. Lo Presti et al. [140] proposed the following relationship:
NG,1 = C0.5 ,
(3.17)
(3.18)
25
2
Holocene sand
(thin-wall sampling)
In-situ frezzing
method
1 Holocene clay
Soft rock
(block)
Pl
e
(tu isto
Soft rock
be ce
sa ne
(rotary core)
m sa
pl
n
in d
g)
0
10
100
0.10
0.25
1.50
2.00
150
300
Remolded
cemented
sandy soils
450
Range from
ROSRINE
study
600
750
General trend
900
1050
0.0
1000
0.5
G0,Field [MPa]
1.0
1.5
Figure 3.10: Differences in the ratio of labratory-to-field stiffness. Left: Data from
Japanese case studies (after Toki et al. [180]). Right: Results from the US
American ROSRINE study (after Stokoe & Santamarina [173]).
200
first increase of g
180
re-increase of g
160
140
reduction of g
plateau
increase of
gS > gprestarin
prestraining
amplitude
-4
gprestarin=10
120
10 -6
10 -5
10 -4
10 -3
Figure 3.11: Influence of cyclic prestraining on small-strain stiffness (after Wichtmann &
Triantafyllidis [193]).
Although experimental data often confirms Masings rules, the threshold shear strain
0.7 in unloading-reloading increases not always to the extent described by Masing. An
example presented in Wichtmann & Triantafyllidis [193] is shown in Figure 3.11. Here,
a sample was cyclic loaded to investigate the effects of cyclic prestraining. In this particular example, the modulus G0 slightly decayed in cyclic loading. In other samples,
which are not shown here, it also slightly increased so that no correlation between the
number of applied prestraining cycles and G0 is found. Part one of Masings rule on the
other hand is confirmed. Normalizing all load cycles for their individual maximum also
confirms the second part.
26
400
200
S
VX,v
-400
-200
200
200
-200
VX,h
VX
VX,v
-400
400
VX
-200
200
VX,h
VX
400
-200
VX,h
-400
sV=100 kPa; sH=100 kPa
VX,v
-400
sV=150 kPa; sH=75 kPa
Figure 3.12: Polar plot of wave speeds as a function of their inclination ( = 0 := vertically travelling waves, = 90 := horizontally travelling waves). Isotropic
loading yields almost isotropic wave speeds (left); axial unloading to a stress
ratio of 2 decreases axial stiffness (right). Dotted lines give the results from
a cross anisotropic model after Bellotti et al. [14].
27
gS,1
GHigh,1
aG,1
1
GHigh,2
1
GLow,1
aG,2
GLow,2
0
gS,Low
gS,High
-5
gS
.4
<1
-5
e
2.5
-5
e
9.7
<
e
6.0
g S<
-4
< gS
.2e
<1
10
20
30
40
50
Figure 3.13: The strain-rate shear modulus parameter and its dependence on strain amplitude, and plasticity Index (PI). After Vucetic et al. [190].
(GS )
.
log
(3.19)
In this way, the strain rate shear modulus parameter G defines the slope of the stress
versus strain rate line in a semi-logarithmic format. In plotting the parameter G versus
strain rate and amplitude (Figure 3.13), the following correlation becomes clear: The
strain rate effect increases with strain rate and soil plasticity, but decreases with shear
strain amplitude. Vucetic et al. [190] first quantified the strain rate effect as a function of
plasticity index (PI) and liquid limit wL but with significant uncertainties.
From Figure 3.13 and 3.14 it can be concluded that small-strain stiffness is only a little
biased by common straining rates. This is in agreement with findings from dynamic
tests, where inertia effects have to be considered on top of strain rate effects. Stokoe et
al. [172] conclude from a recent study that for excitation frequencies changing from 1 to
about 100 Hz, G0 increases by about 5% to 30%, with the effect generally increasing with
increasing PI. Many other researchers found similar trends. Georgiannou et al. [44], for
example, report a 10% increase of G0 in different clays and marl due to dynamic loading.
In sands, only a very small or no increase is found (Hicher [57]). The Japanese working
hypothesis that soil parameters for earthquake analysis can be derived from static tests
is therefore validated from a practical point of view.
28
10
Hardrock core
Concrete
Mortar
Ultrasonic
Resonant-column
4
E0 [MPa]
10
10
10
NC Kaolin (CU)
Vallerica clay
10
-5
10
-4
10
-3
10
-2
10
-1
10
10
10
10
10
Figure 3.14: Strain rate effect on the Youngs modulus E0 after Tatsuoka [177].
29
(3.20)
Anderson & Woods [5] determined the dimensionless soil parameter IR for many clays
within the band 300 < IR < 1800. Considering that the undrained shear strength of
most clays decays to cu 1.7 kPa when they approach the liquid limit [199], their shear
modulus can be as low as 500 kPa. Such low shear moduli at the liquid limit have, for
30
s
p0
pref
(3.23)
can be used to estimate E0 and G0 in various soils. Here again, e is void ratio, p0 is the
mean effective stress in kPa, and pref = 100 kPa is the reference pressure equal to the
atmospheric pressure. The relationship by Biarez & Hicher was proposed for all soils
with wl < 50%, the one by Hardin & Black was derived for undisturbed clayey soils and
crushed sands.
31
Table 3.2: Relations for the shear modulus G0 of CLEAN SANDS and GRAVELS.
Soil tested
Kenya carbonate sand
D50
[mm]
UC
[-]
A
[-]
f (e)
[-]
k
[-]
m
[-]
Ref.
0.13
1.86
101-129
e0.8
0.45-0.52
[40]
(2.17e)2
0.41-0.51
[61]
0.50-0.57
[28]
0.50
[88]
0.50
[88]
0.50
[88]
0.16
1.46
71-87
0.19
1.56
84-104
0.20
1.10
80
0.20
1.70
62
0.20
1.10
62
1+e
(2.17e)2
1+e
(2.17e)2
1+e
(2.17e)2
1+e
(2.17e)2
1+e
0.22
1.35
72
e1.3
0.45
[141]
0.27
1.67
72-81
0.50-0.52
[96]
0.27
1.28
64-69
0.55-0.56
[96]
0.31
1.94
80
0.47
[60]
0.32
2.80
48
0.50
[88]
0.50
1.33
61-64
0.44-0.53
[61]
0.54
1.50
71
0.43
[141]
Silica sand
0.55
1.80
275
(2.17e)2
1+e
(2.17e)2
1+e
(2.17e)2
1+e
(2.17e)2
1+e
(2.17e)2
1+e
(2.27e)2
1+e
(1.46e)2
1+e
0.42
[193]
0.55
1.66
79-90
e0.8
0.43-0.48
[40]
(2.17e)2
0.44-0.53
[61]
0.50
[54]
0.62
1.11
82-130
0.72
1.20
69
1+e
(2.17e)2
1+e
0.75
4.40
71
e1.3
0.62
[141]
Decomposed granite
1.30
75
45
0.88
[86]
1.73
1.33
53-94
0.45-0.51
[61]
Chiba gravel
7.90
10
76
(2.17e)2
1+e
(2.17e)2
1+e
0.50
[122]
32
A
[-]
f (e)
[-]
k
[-]
m
[-]
Ref.
74
e1.27
NA
0.46
[139]
10-40
44
e1.11
NA
0.58
[139]
15-34
50
e1.33
NA
0.40
[139]
24
40
0.202
0.65 [187]
27
44
...
0.85
[147]
Kaolin clay
35
45
(2.97e)2
1+e
NA
0.50
[110]
23-46
50
e1.43
NA
0.44
[139]
41
13
0.252
0.76 [187]
44
52
e1.30
NA
0.50
[139]
45-75
64
e1.52
NA
0.40
[139]
4.5
(4.40e)2
NA
0.50
[110]
PI
[%]
Bentonite
60
1+e
Table 3.4: Proposed relationships for the shear modulus G0 of entire soil groups.
Soil tested
emin
[-]
emax
[-]
A
[-]
f (e)
[-]
m
[-]
Ref.
0.5
1.1
57
(2.17e)2
1+e
0.40 [74]1
0.4
1.8
59
0.50 [16]2
0.6
1.5
33
0.6
1.5
16
Loess
1.4
4.0
1.4
1
e
(2.97e)2
1+e
(2.97e)2
1+e
(7.32e)2
1+e
0.50 [51]3
0.50 [90]4
0.60 [92]5
33
200
Kim et al.
150
Kokusho et al.
Sand &
Gravel
100
Kenya
sand
Clay
50
0
0.5
1.5
400
Kim et al.
Kokusho et al.
300
200
100
2.5
0
0.5
1.5
2.5
Figure 3.15: Shear modulus G0 versus void ratio e, after Table 3.2, 3.3, and 3.4.
qt0.695
e1.13
(3.24)
where qt is the corrected tip resistance in units of MPa. The correlation by Mayne & Rix
[116] and that by Lunne et al. [108] are shown in Figure 3.16.
Since G0 is related to CPT tip resistance, which can be correlated to the undrained
shear strength, G0 should also be directly related to the undrained shear strength. Larsson & Mulabdic [100] proved this for Scandinavian clays by formulating:
A
G0 =
+ B cu
(3.25)
PI
where PI is the plasticity index, cu is undrained shear strength in kPa, and A, B are
empirical soil parameters. For Scandinavian clays they found a reasonable fit for A =
208, and B = 250.
34
Gravelly
sands
100
Sand
nte
me
nce
gu
n
You
Sand mixtures
10
oils
ds
1000
Silt mixtures
100
.
e=0
10
Multiple
Regression
n=481
r=0.901
.
e=1
.0
e=2 0
.
e=3
Clays
G0=49.4 qt
0.695
-1.130
1
1.0
10.0
100.0
0.1
1.0
10.0
Figure 3.16: Correlation between CPT results and very small-strain stiffness after Lunne
et al. [108] (left) and Mayne & Rix [116] (right).
Most correlations between standard penetration test (SPT) results and very smallstrain stiffness take the form:
G0 = A((N1 )60 )B
(3.26)
which was first proposed in Ohsaki & Iwasaki [132]. In Equation 3.26, (N1 )60 is the number of blows per 0.3 m, and A and B are again empirical soil parameters. The relationship
is generally limited to a maximum blow count of (N1 )60 = 30. For cohesive soils, Imai &
Tonouchi [69] found a best fit to test data for A = 15.56 and B = 0.68 (G0 in MPa units).
3.3.2.2 Estimating G0 from conventional tests - Chart by Alpan
Alpan [3] published the chart shown in Figure 3.17 which, as he stated, relates static to
dynamic soil stiffness (see axis labels Es , Ed ). Actually, this chart relates stiffness from
conventional laboratory tests to very small-strain stiffness: The dynamic modulus is the
initial, or very small-strain modulus E0 , and the static modulus is the apparent elastic
Youngs modulus in conventional soil testing, e.g. at axial strains of a 1 103 in
triaxial testing. For a soil with known Youngs modulus in triaxial unloading-reloading,
the Alpan chart can thus provide an estimate for its very small-strain modulus E0 .
35
Figure 3.17: Correlation between very small-strain stiffness and stiffness at larger strains
from conventional laboratory tests after Alpan (10 kg/cm2 1 MPa).
directly, others use more complex fitting procedures. Only the former are discussed in
this section. The latter are discussed in detail in Chapter 5.
From test data, sufficient agreement is found that the stress-strain curve for small
strains can be adequately described by a simple hyperbolic law. An analogy to the hyperbolic law for larger strains by Kondner [94] (see Section 7) was proposed by Hardin
& Drnevich [52]:
G
1
.
=
(3.27)
G0
1+
r
max
G0
(3.28)
with max being the shear stress at failure. Essentially, Equation 3.28 relates large (failure)
strain to small-strain properties, which astonishingly often works well. More straightforward, and less prone to error is probably the use of a smaller threshold shear strain.
Stokoe et al. [89] for example, suggest to use the shear strain r = 0.5 at which the shear
modulus G0 is reduced to its half. Additionally they modify Equation 3.27 by introducing an additional exponent = 0.92:
1
G
.
=
G0
1 + r
36
(3.29)
0
Secant modulus
G/G0 [-]
1.0
0.8
0.6
0.4
Hardin &
Drnevich
0.2
0.0 -3
1e
1e
-2
1e
-1
1e
1e
1e
1e
0.8
0.6
0.4
Santos &
Correia
0.2
0.0 -3
1e
-2
1e
-1
1e
1e
1e
1e
1e
Figure 3.18: The original Hardin & Drnevich [52] relationship (left) and its modification
by Santos & Correia [155] (right) compared to test data.
Using the threshold shear strain a = 0.7 , Santos & Correia [155] have put forward the
following modified Hardin & Drnevich relationship:
G
1
.
=
G0
1 + a a
(3.30)
The constant a in Equation 3.30 is problematic though: For a 6= 73 the relationship is inconsistent as the shear modulus G is not reduced exactly to 0.7G0 at shear strain = 0.7 .
This for example is the case for the best fit a = 0.385, Santos & Correia [155] found from
a correlation using many test results. The modification by Santos & Correia combined
with a constant of a = 73 on the other hand is identical to the Hardin-Drnevich relationship. Looking at both reduction curves (Figure 3.18), their deviation is < 3% and thus
negligible for practical purposes anyway. The approach by Santos & Correia gives a
slight upshift of the reduction curve which is indeed desirable in comparison to the not
truly hyperbolic test data collected by Santos & Correia [155] shown in Figure 3.18).
In the following, mostly the modified Hardin-Drnevich relationship (3.30) with a =
0.385 is used. For the original Hardin-Drnevich model, the reference shear strain a =
0.7
is then taken at the shear strain where G decayed to 0.722G0 . For the Santos & Correia
upshift, a = 0.7 where G decayed to 0.7G0 is to be selected. The exponent proposed
by Stokoe et al. is not adopted as it may sometimes flatten the shear modulus reduction
curve too much in the small-strain range.
In Equation 3.28 the threshold shear strain r is related to the shear stress at failure. In
the Mohr-Coulomb failure criterion [123] the shear stress at failure can be expressed as
follows:
1
(3.31)
max = (2c(1 + cos 2) + (10 + 30 ) sin 2) .
4
of normally
Remembering that a = amax /G0 , the threshold shear strains 0.7 , and 0.7
37
0.385
(2c(1 + cos 2) + 10 (1 + K0 ) sin 2) ,
4G0
(3.32)
0.7 =
3
(2c(1 + cos 2) + 10 (1 + K0 ) sin 2) ,
28G0
(3.33)
and
3.3.4 Summary
The stiffness of soils at very small strains is affected by a number of parameters. Its
correct determination is not a trivial task: Test results on disturbed samples may give
very different results than in-situ tests. In-situ testing methods use geophysical methods,
which are not yet common in geotechnical engineering.
However, it has been shown that many of the factors which affect the shear modulus
G0 and its reduction can be regarded as being of minor importance. Neglecting these
factors leads to the simple relationships shown in Table 3.4. and illustrated in Figure
3.15. Empirical relationships for estimating the shear modulus reduction with strain
are available as well. From a practical point of view, small-strain stiffness can thus be
incorporated into routine design at minimum cost. The use of laboratory and/or in-situ
testing data is still preferably though.
If no in-situ test data are available, laboratory small-strain stiffness data can be extrapolated to the field. In the case of common tube or rotary core sampled cemented soils,
the actual soil stiffness in the field G0,F ield should be expected to be considerably higher
than that measured in the laboratory G0,Lab (Figure 3.10). The modulus reduction curve
found in the laboratory can be extrapolated to the field by assuming [173]:
G,F ield =
G0,F ield
G,Lab .
G0,Lab
(3.34)
A numerical model that allows for the incorporation of small-strain stiffness into routine design is developed in the following chapters. The input parameters of this model
will be exactly the ones quantified above: The shear modulus G0 and the threshold shear
.
strain 0.7 or 0.7
38
Chapter 4
Small-strain stiffness at the soil particle level
The most important parameters that affect small-strain stiffness have been identified in
the previous chapter. In summary, these are void ratio, confining stress, cementation
and strain amplitude. The observed influence of these parameters on material stiffness
are now interpreted at the soil particle level. However, it is not the aim of this chapter
to derive a complete small-strain stiffness model from micromechanical considerations.
It is aimed rather at providing a thorough understanding of the small-strain stiffness
phenomenon.
The material, which we commonly refer to as soil is usually composed of air, water,
and solid particles, which may vary in size, shape, and mineralogic composition. The
mechanical properties of soil depend directly on the interaction of these phases with
each other and with externally applied potentials such as stress and temperature. Some
mechanical properties of soil, for example suction, cannot be understood without explicitly considering multi-phase physics. Others, that are purely determined by the net
forces between particles, for example small-strain stiffness, can be sufficiently explained
by the micromechanical considerations presented in this chapter.
39
Figure 4.1: Scanning electron microscopy images of sand Particles: Single-grain sand
fabric (left); Close-up of amorphous silica and quartz bonds in a sandstone (right). Photographs by IWF, Universitat Marburg, and Tore Stendahl,
NTNU, Trondheim, Norway.
Clay matrix
Silt or
Sand
(a)
(b)
(c)
Figure 4.2: Particle arrangements after Mitchell [120]: (a) single-grain sand or silt fabric,
(b) multi-grain clay fabric, (c) mixed fabric.
In contrast to sand, clay has typically a multi-grain fabric. If not an artefact from a
preexisting rock, multi-grain clay fabrics develop due to high surface forces and chemical activity. Assemblages of clay particle groups or clusters can be arranged within clay
fabric in several ways, i.e. edge-to-face or face-to-face. The particle arrangement within
the clusters themselves depends mainly on the clays mineralogy. Possible inter-particle
and intercluster forces in clay, in addition to the aforementioned ones, include electrostatic forces in the form of ionic bonds, hydrogen bonds, and hydration forces, as well as
capillary stresses.
Although sand and clay structures are different in many ways, their qualitative smallstrain behavior is not (see Chapter 3). Both soil types are modeled as granular assemblies
here, which represents the level of single-grain sand fabric, and multi-grain clay fabric.
Figure 4.2 gives an illustration of the above introduced sand and clay fabrics and the
more common case of mixed fabric.
40
t
th
c -contact
branch vector li
x
Figure 4.3: The cth inter-particle contact in global and local coordinates.
M
1 X c c c
=
l k l
V c=1 i jl k
(4.1)
where V is the representative volume of the granular assembly, M is the total contact
c
is the
number in the volume V, lic is the branch vector as introduced in Figure 4.3, and kjl
contact stiffness tensor in the global coordinate system. In a local coordinate system the
contact stiffness tensor can be written in the form:
c
kjl
= knc ncj ncl + ksc scj scl + ktc tcj tcl
(4.2)
where knc , ksc , and ktc are the local contact stiffnesses along directions nc , sc , and tc respectively. Several homogenization techniques have been proposed in literature to rewrite
Equation 4.1 in an integral form. However, all these homogenization techniques involve
the definition of density functions for directional contact frequency, branch lengths, and
the definition of local stiffness properties. Postulating these density functions is not an
easy task for a non-ideal particle assembly: The perturbation function of branch vector
lengths for example should reflect particle gradation, shape and preferred particle orientation. As an alternative to explicitly introducing these density functions, numerical
discrete element codes can be employed. At the cost of calculation time, these codes
keep track of all inter-particle contacts.
Either way, using homogenization techniques or discrete element codes, deriving the
global material stiffness from Equation 4.1 is an elaborate task. Unfortunately, not even
the constitutive law for the inter-particle contacts is known exactly. Often the HertzMindlin theory (see Section 4.5) is used in combination with correction factors, bonding,
viscous components etc. Relating local inter-particle parameters to global soil parameters from laboratory or field testing is also difficult. Largely due to these reasons, micromechanical models are not yet used in engineering practice. Nevertheless, they can
41
42
a1
b1
a2
b2
Figure 4.4: Disturbance of inter-particle forces in a biaxial test. A bonded discrete element model is first consolidated (a) and then biaxially sheared (b). From
hydrostatic loading equally distributed chains of contact forces (a1) become
concentrated during deviatoric loading (b1). At the same time, most of the
initially sticking contacts (a2) slide so that the bonds break up (b2).
43
2GS
a,
1 S
(4.3)
where GS , and S are the elastic constants of the sphere material, and a is the radius of
the circular contact area, which is calculated as:
s
3N R(1 S )
a= 3
.
(4.4)
8GS
The normal contact force N , the contact radius a, and the sphere radius R are illustrated
in Figure 4.5. Contact shear stress distributions caused by different amounts of slip xT
inbetween the spheres are shown in Figure 4.5. The partial slip introduced in Figure 4.5
is a refinement of the previously discussed simple Coulomb-type frictional law. Energy
dissipation and stiffness reduction start with the onset of partial slip. Johnson [85] gives
as the resulting tangential contact stiffness:
s
2(1 S ) 3
T
ks = kt = kn
1
,
(4.5)
2 S
N tan S
where S is the inter-particle friction angle, and T is the shear force in the inter-particle
contact. Both normal and shear stiffness therefore increase with the normal contact force
to the power of m = 31 .
The empirical power law by Hardin & Richard [54] that relates initial stiffness to confining pressure (Equation 3.5) can therefore be readily understood from a micromechanical point of view. The power law exponent observed in experiments is generally somewhat higher though. More sophisticated, and for soil particles probably more adequate,
contact theories, yield these higher exponents. Rough surface contacts for example yield
a power law exponent of m 0.5 0.6 [197]. Meta-stable cone to sphere contacts yield a
power law exponent of m = 0.5 [45], too. Higher power law exponents for large strains
can be explained by slippage between the particles.
Local contact stiffness is not an explicit function of void ratio as can be seen in Equation
4.3, and 4.5. Therefore, the effect of void ratio is studied here by comparing two different
cubical assemblies: simple cubic (SC) packing and face-centered cubic (FCC) packing of
44
xT/2
xT/2
T
T
N
T
T
T
xT/2
xT/2
gross slip
contact
stress
partial slip
xT/2
T
T
b
a
xT/2
xT/2
b
Gross
slippage
No slippage
Zone of
slippage
Zone of
no slippage
2R
Unit cell SC
2 2R
Figure 4.6: Simple cubic packing (SC) versus face-centered cubic packing (FCC).
45
(k)F CC
=2 2
(k)SC
2 2
13
= 2.
(4.6)
Assuming Poisons ratio to be constant, the empirical relationship given by Hardin &
Richart [54] presented in Equation 3.8, suggests a bulk stiffness increase of 65% in the
FCC packing. Compared to the empirical value, the above deduced stiffness increase of
100% is too high. The difference is probably due to the simplification of real soil structure
in SC and FCC assemblies. Nevertheless, the stiffness increase is of the same magnitude.
The number of contacts per particle and particle bracing can therefore be assumed to be
the main factors, which lead to higher stiffness in real soils.
32
GS
3(1 S ) N
2
(kn )af ter =
(CC + 1) 3 1 +
(4.9)
1 S
2
4R2 GS
where kn is the local normal contact stiffness as introduced in Equation 4.2. In both cases,
cementation is causing a stiffness increase independent of confining pressure. With in-
46
R
N
Soil
particle
TN
a
R+t
a0
Cementation
Figure 4.7: Cemented inter-particle contact of two spheres after Fernandez & Santamarina [39].
4.2.4 Summary
The influences of parameters that appear in almost all very small-strain stiffness relationship (see Table 3.2, 3.3, and 3.4) were studied at the inter-particle contact level. Most
previously presented experimental findings were qualitatively reproduced. In order to
obtain a quantitative reliable constitutive response of a micromechanical model, a number of conditions have to be fulfilled: First, the inter-particle contact needs to be defined
adequately, second the soil grains and its structure have to be reasonably well abstracted
and finally the problem of homogenization or calculation time must be solved. For the
moment, a small-strain stiffness model aimed for practical use can hardly be formulated
from a particle contact level. The following chapters therefore concentrate on continuum
models. Existing continuum small-strain stiffness models are discussed in Chapter 5. A
new small-strain model is developed in Chapter 6.
47
Chapter 5
Existing small-strain stiffness constitutive
models
The first small-strain stiffness constitutive models that were intended for static soil analysis date from the late 1970s. This chapter gives an overview of the developments in
the field of small-strain stiffness since then. It begins with a historical review followed
by a discussion of the most well known models and modeling directions. Some of the
models that will be introduced in this chapter should be classified rather as relationships
because they lack a constitutive response in non-monotonic loading patterns. Nevertheless they are included in this chapter, since they have been extensively used in the past
despite this shortcoming.
The models and modeling directions covered in more detail are:
The Simpson Brick model
Models known from soil dynamics
The Jardine model
Multi surface models
Intergranular Strain
The order of the above list is chosen according to model complexity rather than according the date they were proposed. Only continuum models are considered in this
chapter. Although micromechanical concepts might be very helpful in the model abstraction process, they cannot yet be used in the analysis of boundary value problems.
49
50
log (strain)
Figure 5.1: Simpson brick model with four bricks (left). Stiffness decay after load reversal
is a discontinuous approximation of the real stiffness-strain curve (right).
the model by Lansivaara [105], its basic idea of history tracking has been widely used.
Doing so, the number of bricks are often reduced to one only.
Bilinear:
Ramberg-Osgood:
Hardin-Drnevich:
G0
for < y
y + G0 ( y ) for y
=
1 +
G0
y
G0
=
1 + r
=
(5.1)
(5.2)
(5.3)
where G0 is the initial shear modulus, , and 1 are material constants, and y , y ,
and r are user defined threshold values for shear- strain and stress respectively. Rewriting the Ramberg-Osgood, and Hardin-Drnevich models in terms of normalized secant
51
t
G0
1
t
G0
G0
ty
G0
tf
ty
gy
gy
gr
Figure 5.2: Skeleton curves in -space of a) the bilinear model, b) the RambergOsgood model, and c) the Hardin-Drnevich model.
stiffness,
G
G0
(Ramberg-Osgood)
(5.4)
G
1
=
G0
1 + r
(Hardin-Drnevich).
(5.5)
The Ramberg-Osgood defines stiffness decay as a function of stress whereas the HardinDrnevich model defines it as a function of strain. Other than that, the two formulations are very similar. The hysteresis effect described in Masings second rule (see Section 3.2.7) can be easily incorporated in both models using the expression proposed in
Hashiguchi [55]:
0 R
0 R = f
(5.6)
L
with shape Factor L = 1 in primary loading and L = 2 in unloading and reloading.
In Equation 5.6, f () represents the monotonic stress-strain relationship = f (). The
strain and stress levels 0 and 0 denote the last load reversal point in stress-strain space,
where R , and R specify the actual strain and stress respectively. Figure 5.2 shows examples of hysteresis loops derived by the bilinear, Ramberg-Osgood and Hardin-Drnevich
models.
The Ramberg-Osgood model has the probably unintended quality that it can be easily used in combination with stress space based elastoplastic models, as shown in Papadimitriou et al. [134]. For the coupling of the Ramberg-Osgood formulation to their
elastoplastic model, Papadimitriou et al. defined two scalar valued functions for the
52
(e elr ) : (e elr )
(5.7)
where r and e are deviatoric stress and strain ratio tensors respectively and superscript
() denotes tensors that are associated to the last memorized load reversal point. A load
reversal is assumed whenever the scalar valued function e looses its monotony, e.g.
e e < 0. In terms of a one-brick Simpson model in strain space, a load reversal is
thus assumed whenever the string between the man and the brick is slacking. Then the
position of the brick is memorized and stiffness is calculated as a function of the distance
to this memorized point in stress space, e.g. as the resulting stiffness of the RambergOsgood formulation 5.4.
Small-strain stiffness models similar to the one published by Papadimitriou et al. can
be found in Hueckel & Nova [65] or Pestana & Whittle [137]. The model of Hueckel &
Nova mainly differs from the one by Papadimitriou et al. in the stress memory. The
model by Hueckel & Nova memorizes not only one, but several load reversal points
and does not explicitly use the Ramberg-Osgood formulation. Likewise, the hysteretic
equations of the MIT-S1 model described in Pestana & Whittle [137] as:
(
1
for unloading
2G/K
1+r
=
(5.8)
1
for reloading
2G0 /K0
1+p r
lr
53
2B
ea= C
log10 e a
Projected minimum
Figure 5.3: Jardine model - secant stiffness defined as a trigonometric function of strain
(projected curve).
(5.10)
Today, this type of non-linear elastic stress-strain law is commonly known as the Jardine
model. Its empirical model constants A, B, C, , and can be determined from (smallstrain) triaxial test data. The secant Youngs modulus Eu is normalized with undrained
shear strength cu . Due to the Jardine models trigonometric nature, it is mandatory to
specify the exact strain range in which the non-linear stress-strain law is to be applied.
When exceeding the upper or lower limit of the strain range specified, stiffness is set as
constant (Figure 5.3).
The normalized tangent Youngs modulus Eut corresponding to the secant Youngs
modulus Eu can be derived by differentiating and rearranging Equation 5.10:
1
BI
Eut
= A + B cos(I )
cu
2.303
sin(I ),
(5.11)
where I = log10 Ca . Since most numerical procedures make use of the tangent modulus Eut instead of the secant modulus Eu , Equation 5.11 represents the more commonly
used form of the Jardine model. Jardine et al. [79] also mention in their original work
the possibility of expanding the proposed nonlinear stress-strain law to bulk and shear
moduli. Potts & Zdravkovic [138] propose the following notation for such an expansion:
54
c1
3G
p0
2
= C1 + C2 cos(c1 X c2 ) C2 c1 c2 X2.303
sin(c1 X c2 )
K
p0
4
= C4 + C5 cos(c3 X c4 ) C5 c3 c4 X2.303
sin(c3 X c4 ) with
1
2 2 [(1 2 )2 (2 3 )2 (1 3 )2 ]
2d
= log10 3C3 = log10
,
3C3
X
Y
c1
and
(5.12)
where G is the shear modulus, K is the bulk modulus, p0 is the mean effective stress,
i are the principal strains, and the number of material parameters (C1 C6 , c1 c4 ) is
increased to 10. In its expanded version, the Jardine model implicitly varies Poissons
ratio as a function of strain. When selecting the models parameters the user needs to
be aware of this in order to keep Poissons ratio within its physical limits. Since the
Jardine model is often used in Great Britain, model parameters for many British clays
are available in the literature.
fn
= f (
ij )n rn f (q?,n )
(
ij )n = ij (ij )n
with
(5.13)
where rn is the ratio of the nth yield surface to the outermost, or distinct yield surface,
(ij )n gives the internal backstress of the nth yield surface, and q?,n is a suitable set of
scalar internal variables.
One of the most important features of multi-surface models is their ability to capture
recent stress history in their internal variables (mainly backstress). Small-strain stiffness
models based on the multi-surface concept associate small-strain stiffness with the innermost yield surface. By defining an appropriate set of hardening laws for all yield surfaces, such multi-surface models can give realistic predictions over a wide strain range.
A special subclass of kinematic multi-surface small-strain stiffness models is the class
of bubble models (Al-Tabbaa [1], [2]). Bubble models have only one small kinematic yield
surface (bubble), with elastic stiffness equal to the soils small-strain stiffness. Increasing
the number of kinematic yield surfaces to two, Stallebrass & Taylor [171] extended the
bubble model by Al-Tabbaa & Wood [2] to the 3-SKH model. Although the 3-SKH model
might be considered sufficiently complex, the number of yield surfaces could also theoretically be increased to infinite.
55
q
D
40
40
20
0
-20
20
A
20
B
40
C
60
80
100
20
40
60
80
100
-20
-40
-40
Figure 5.4: Position of 10 yield surfaces after stress path 0ACBDB. The inner surfaces
have been dragged by the stress point. The outermost yield surface has never
been reached because of volumetric hardening: Plastic volumetric strains on
inner yield surfaces increase the size of all yield surfaces.
An illustration of the multi-surface concept is given in Figure 5.4. Here the evolution of
10 yield surfaces for a given load path is shown, using a discrete form of Houlsbys [62]
continuous hyperplastic formulation. Note that this formulation abandons the concept
of nested and non-overlapping yield surfaces, which goes back originally to Prevost
[142].
||||
R
(5.14)
where R is a material constant associated with the very small-strain range. Evolution of
Intergranular Strain is controlled by the loading direction; loading history is memorized
as Intergranular Strain. Again, a one-brick like loading-unloading criterion is used:
= (1 ) : D f or : D > 0
(5.15)
D
f or : D 0
where D is the stretch tensor and R is a material parameter that controls the evolution
of Intergranular Strain. At the same time, R controls the nonlinearity of the resulting
stiffness-strain curve. Numerically integrating Equation 5.15 over time gives continuous
Intergranular Strain with 0 1. For = 0 stiffness is increased in all loading directions; for = 1 stiffness is a function of loading direction. Minimum stiffness is obtained
56
Figure 5.5: Interpolated stiffness M as a function of angular change in the loading path
according to Equation 5.16 with L = 1, N = 0, mR = 5, and mT = 2. Response
curves are drawn at Intergranular Strains of = 0 to 1 in steps of 0.1 for
= 6.0 on the left and = 2.0 on the right respectively.
for monotonic loading, neutral loading gives a somewhat increased stiffness defined by
the multiplier mT whereas reversed loading gives maximum stiffness by applying the
multiplier mR . Loading directions in between are interpolated according to:
M = [ mT + (1 )mR ]L +
(1 mT )L : + N f or
(mR mT )L :
f or
: D > 0
: D 0
(5.16)
where M is the interpolated stiffness and is the direction of Intergranular Strain defined as = /||||. Material parameters mR , mT and affect the interpolation, where
purely increases the impact of on the interpolation and so together with R controls
the non-linear stiffness decay. The hypoplastic operators L and N represent the linear
and the nonlinear terms in the original hypoplastic equation respectively [93]. Assuming a straight monotonic loading history in a plane problem, the directional dependency
of M on the loading direction can be visualized as shown in Figure 5.5. For increasing
Intergranular Strain the directional dependency grows to its maximum at = 1.
The high number of parameters and also the model abstraction might pose a serious
problem for routine application of Intergranular Strain. Engineers working in practise
may simply not understand the meaning of the parameters they should provide. Different aspects of the practical usability of all the models introduced so far are discussed in
more detail in the next section.
57
58
Chapter 6
The Small-Strain Overlay model
All isotropic hardening, and most kinematic hardening elastoplastic models use yield
loci of finite size. Strategies to implement small-strain stiffness in just one very small
elastic domain (bubble models), or in many yield loci (multi-surface plasticity), have
been discussed in the previous chapter. However, models that are used in engineering practice today, rarely have more than one elastic domain. The Cam-Clay, or the
Hardening Soil model are just two examples of such models. It is impossible for very
small-strain stiffness to be associated with the entire elastic domain in these models as
they would respond too stiff in larger unloading-reloading cycles as a consequence.
Therefore, their elastic stiffness is generally taken as secant stiffness of larger strain
unloading-reloading loops. Small-strain stiffness is neglected. Here, the Small-Strain
Overlay model which will be developed in this chapter, comes into play.
The Small-Strain Overlay model is, similar to some of the models discussed in the previous chapter, a paraelastic approach. The overlay model can be used in combination
with elastoplastic formulations: If it detects small or very small-strain amplitudes, it increases the stiffness of the elastoplastic model accordingly. For larger strain amplitudes,
the overlay model is inactive, but still monitors loading history.
In order to guarantee the new models compatibility with many existing constitutive
models, it was decided to keep its output isotropic. Its objective multi-axial formulation
though accounts for strain induced anisotropy. In this chapter, the small-strain stiffness
overlay model is first formulated, and secondly validated.
59
(m)
e kl (m) kl Sl = 0
(6.1)
a set of orthogonal Eigenvectors Slm for e kl is found. Unlike the transformed deviatoric
strain increment
e kl = Skm e mn Snl ,
(6.2)
the transformed strain history
H kl = Skm Hmn Snl ,
(6.3)
1
1 + u((2) H22 )( H22 + 1 1)
= H22
+1
1
(3)
1
+
u(
H
)(
= H33
H
+
1
1)
33
33
+1
T 11 =
T 22
T 33
60
1
H11 +1
(6.4)
e1
e1
e2
e2
e3
e3
e1
e1
e2
e2
e3
e3
Figure 6.1: Strain history interpreted as 2nd order surface in strain space. From left to
right: (a) hydrostatic loading, (b) triaxial extension, (c) simple shear, (d) simple shear after triaxial extension.
61
u(x) =
0 f or x < 0
1 f or x 0
(6.5)
(6.6)
From the updated strain history a scalar valued shear strain measure Hist is defined
next:
ke H k
Hist = 3 km ml
(6.7)
ke kl k
where k...k denotes the Hilbert-Schmidt norm kAk = aij aij , which can be proven to
conserve objectivity. In a geometric context, Equation 6.7 is the projection of the strain
history onto the actual loading direction. The shear strain measure Hist can then subsequently be used to define an isotropic material stiffness as shown in the next section. In
an incremental form, the model formulation is summarized in Box 6.1.
G0
,
1 + a Hist
0.7
(6.8)
where G is the actual secant shear modulus, G0 is the initial shear modulus at very smallstrains, and Hist is the monotonic shear strain defined above. The parameter a and
the threshold shear strain 0.7 is discussed in more detail in Section 3.3.3. In numerical
applications, the secant modulus given in Equation 6.8 must be converted to a tangent
modulus:
2
0.7
.
(6.9)
G = G0
0.7 + aHist
62
Gi+1 =
i+1
Hist
G0
i
Hist
1
i+1
Hist
0.385 i+1
+ 0.7Hist
i
Hist
0.385 i
+ 0.7Hist
(6.10)
where i and i + 1 denote quantities of the previous and the actual calculation step respectively.
A complete numerical implementation of the Small-Strain Overlay model that integrates over the stiffness reduction curve is shown in Appendix A. Here, a lower cut-off
of the models stiffness reduction curve is introduced. The cut-off should be activated in
combination with hardening plasticity formulations as otherwise, material stiffness for
larger strains would become too small (see also Section 7.3.1).
Bulk modulus
Assuming Poissons ratio to be constant (see Section 3.3) the bulk modulus can be calculated from the actual shear modulus as:
K=G
2(1 + )
.
3(1 2)
(6.11)
63
1.0
t/ty
0.8
t/ty
C
1
0.8
A
0.6
0.6
3
0.4
0.4
0.2
0.2
g/ gy
-8.0
-6.0
-4.0
-2.0
2.0
4.0
6.0
8.0
g/ gy
-8.0
-6.0
-4.0
-2.0
2.0
4
-0.2
-0.4
-0.6
-0.6
-1.0
6.0
8.0
-0.2
-0.4
-0.8
4.0
-0.8
0
-1.0
Figure 6.2: Hysteresis loops in symmetric (a) and irregular (b) loading according to the
extended Massing rules (after Pyke [143]).
All four rules together are also referred to as the extended Masing rules. A discussion
of several possible extensions replacing the fourth of the above rules, which goes back
to Rosenblueth & Herrera [152], is given in Pyke [143]. Jennings [81], [82] for example
proposes that the upper and lower bounds of the stress-strain relationship could be given
by the reloading and unloading curves corresponding to the previous minimum and
maximum strain amplitudes. Richart [149] on the other hand, simply limits the shear
stress to a value no greater than the shear strength.
The Masing rule extensions by Rosenblueth & Herrera, Jennings and Richart are illustrated in Figure 6.2 for a one-dimensional simple shear problem. If point 2 is the previous minimum strain amplitude and curve 3 is the greatest previous reloading curve,
on reloading from point 4, the solutions suggested by Rosenblueth & Herrera and that
proposed by Jennings would follow path A. If, however, the dashed curve 0-1 is the
greatest previous reloading curve, the solution used by Jennings follows stress path B.
The solution used by Richart always follows path C.
The Small-Strain Overlay model adopts the first Masing rule by assuming G0 = constant in Equation 6.8. The second rule is obeyed by assuming a shape factor of 2 in the
relationship given by Hashiguchi [55] (see Equation 5.6) i.e. (0.7 )ur = 20.7 where (0.7 )ur
is the threshold shear strain in unloading and reloading. The third and fourth rule request an internal memory of previous changes in the stress-strain path. For all of the
above introduced Masing rule extensions, the memory i.a. requires loci in stress space.
64
65
66
i
(m)
H ikl = Skm Hmn
Snl where (ei+1
kl )Sl
kl
=0
1
(1)
T 11 = H11
1
+
u(
H
)(
H
+
1
1)
11
11
+1
1
(2)
H
+
1
1)
T 22 = H22
1
+
u(
H
)(
22
22
+1
1
T 33 = H33
1 + u((3) H33 )( H33 + 1 1)
+1
T ij
= 0 for i 6= j
(
0 f or x < 0
where u(x) =
1 f or x 0
d) Update strain history:
i
i
i
i+1
H i+1
kl = T km (H mn + mn )T nl + ekl kl
G0
i+1
i
Hist
Hist
i+1
Hist
0.385 i+1
1+ 2 Hist
0.7
i
Hist
0.385 i
1+ 2 Hist
0.7
67
Esec [kN/m2]
reloading
4E+5
primary loading
2E+5
0
-5
10
-4
10
-3
10
e1[-]
Esec [kN/m2]
q [kN/m2]
100
4E+5
50
0
2E+5
-50
Calculated
Rivera/Bard (initial)
-100
Rivera/Bard (reloading)
0
-5
10
-4
10
-3
10
e1[-]
-4
-5x10
-4
-3x10
0x10
-4
3x10
-4
5x10
e1[-]
Figure 6.4: Standard drained triaxial test by Rivera & Bard (after Biarez & Hicher [16])
and simulation with the Small-Strain Overlay model (G0 = 185MPa, = 0.20,
0.7 = 2e4 ).
68
q [kPa]
DOX
AOX
150
100
COX
30
BOX
50
Calculated (HS)
20
0
0
100
O
200
A
300
400
p' [kPa]
-50
10
-100
0
10-4
10-3
10-2
e1[-]
-150
Figure 6.5: Triaxial test by Richardson [148] and its simulation with the Small-Strain
4
Overlay model (Gref
0 = 42 MPa, = 0.20, 0.7 = 5.5e ) Left: path dependent
shear stiffness. Right: loading paths in q-p space.
model response of the Small-Strain Overlay model is compared to that of the Intergranular Strain concept.
Both, the model proposed here and the Hypoplastic model with Intergranular Strain
show a similar behavior. The least stiffness is obtained during monotonic loading; highest stiffness is obtained for a full load reversal, which in the case of biaxial loading is
a reversal in both loading directions. Reversing only one biaxial loading direction and
keeping the other constant, results in a somewhat increased stiffness which lies between
the minimum and the maximum stiffness. Here the Small-Strain Overlay model clearly
shows its ability to capture deviatoric strain history. Loading to the side yields less stiffness than reversed loading.
69
2
e2
Figure 6.6: Stiffness-strain decay in a biaxial test on sand. Comparison of the tangent
stiffness generated by the Small-Strain Overlay model to that, generated by a
Hypoplastic model with Intergranular Strain range.
which is reached by isotropic consolidation to 100 kPa and application of 100 kPa vertical stress afterwards. The magnitude of the then applied stress increments is 10 kPa.
Experimentally tested stress increment directions are shown in Figure 6.7. Only the intermediate stress state is considered here.
Doanhs strain response envelopes are derived in the triaxial plane. In the triaxial
plane however, the unloading-reloading criterion of the Small-Strain Overlay model reduces to a one brick criterion. A sharp stiffness jump is obtained whenever the deviatoric
strain increment reverses its direction. In order to avoid sharp stiffness jumps around
hydrostatic loading, an alternative to the above model formulation is explored in Figure
6.7. Assuming Poissons ratio not to be constant, but varying from = 0.05 0.25 in the
case of the sand tested, smooth stiffness transitions could be obtained: The bulk modulus is assumed constant as long as Poissons ratio is within an admissible range, i.e.
min . In the experiment, the bulk modulus is therefore unaffected by the non-linear
small-strain shear stiffness decrease. However, as such large variations in Poissons ratio
within the small strain range is not supported by experimental evidence, it is not further
considered in this thesis.
In Figure 6.7 experimentally derived strain responses are compared to these derived
by the Small-Strain Overlay model. Strain responses after load reversals are modeled
reasonably well except for those which involve more soil plasticity. However, a combination of the Overlay model with an elastoplastic model can also reproduce monotonic
and almost monotonic strain responses.
70
6x10
15
-4
125
(a)
90
90
10
61
125
4x10-4
145
135
145
35
61
2x10-4
0
180
193
-5
180
e1 [-]
s1 [kN]
(b)
135
0x10
215
35
230
215
230
-15
-15
-2x10
305
243
-10
-4
270
-4
-10
-5
10
-4x10
-4
-6x10
15
-4
-4x10
0x10
-4
2x10
-4
4x10
e3 2 [-]
s3 2 [kN]
-4
6x10
-4
-2x10
125
Plot (b):
(c)
135
Test by Doanh
linear elastic
90
small-strain overlay
-4
4x10
145
61
-4
2x10
e1 [-]
180
215
35
230
0x10
Plot (c):
-4
-2x10
Test by Doanh
elastoplastic
elastoplastic with small-strain overlay
-4
-4x10
-4
-6x10
-4
-4x10
-2x10
-4
0x10
-4
2x10
4x10
-4
e3 2 [-]
Figure 6.7: Strain response envelopes experimentally derived by Doanh [32] and their
simulation with the Small-Strain Overlay model: (a) strain increment directions, (b) Small-Strain Overlay model, (c) Small-Strain Overlay combined
with hardening plasticity.
71
200
100
-100
-200
Calculated
Rivera/Bard (Monotonic)
-300
-3
-3
-3x10 -2x10
10
-3
0x10
-3
10
-3
2x10
-3
3x10 e1[-]
Figure 6.8: Strain accumulation in closed cycles due to mean stress dependent elastic
stiffness.
72
where denotes specific entropy, % mass density, temperature, q heat absorption, and
h is the heat flux vector.
Constitutive equations have then not only to be defined for stress, but also for specific
entropy, internal energy, and heat flux. This, however, is not within the scope of the
Small-Strain Overlay model, which simulates near elastic material behavior with a relatively small change in specific entropy compared to the overall entropy change during
plastic loading of the background model. Nevertheless, a rigorously formulated constitutive model should obey thermodynamic principles. The model has not, however, been
demonstrated to obey such principles, although its constitutive parts have been chosen
to represent soil behavior realistically.
An example of the potential problems due to the non-rigorous thermodynamic formulation of the Small-Strain Overlay model is given next. Considered are isothermal
processes, starting from a state of equilibrium only. For such processes, the theorem on
isothermal cyclic processes hold (e.g. Truesdell & Noll [181]: In every isothermal process
starting from a state of equilibrium, the total stress work around a closed path is non-negative.
Partly as a consequence of not having implemented the entire set of Masing rules in
the Small-Strain Overlay model, the irregular cyclic loading situation shown on the left
hand side of Figure 6.9 may violate the theorem on isothermal cyclic loading processes.
The same loading situation, however, with less negative stress work release during unloading does not violate the above theorem (Figure 6.9 right hand side). The unloading
stiffness in the second case is calculated as follows:
Let t0 be the time at the onset of loading. For any time t1 > t0 the total stress work
associated with the overlay model is given by:
Z
t1
W (t1 ) =
(6.16)
t0
(6.17)
The shear modulus Gel (t1 ) of the iso-energetic linear elastic material can thus be explicitly
determined. In order to avoid negative stress work in a closed path simulated by the
73
q [kPa]
250
250
200
200
150
150
100
100
50
50
-4
-2x10
0x10
-4
2x10
-4
4x10
-4
6x10
-4
8x10
e1[-]
-4
-2x10
0x10
-4
2x10
-4
4x10
-4
6x10
-4
8x10
e1[-]
overlay model, its stiffness Gol (t) in unloading should at all times t > t0 be greater than
or equal to that of the corresponding linear elastic material:
(t) < 0.
Gol (t) Gel (t) t > t0 with W
(6.18)
In this way, it is ensured that the total negative stress work during unloading is not
less than that of a linear elastic material with identical internally stored strain energy at
the onset of unloading. Furthermore, as the positive stress work of the materials before
unloading are equal, the total stress work in a closed path cannot be negative.
74
Chapter 7
HS-Small, a small-strain extension of the
Hardening Soil model
In this chapter, the Small-Strain Overlay model is combined with an existing doublehardening elastoplastic model. As the main focus of this thesis is with the development
of a model for routine design, the combination is made with the Hardening Soil (HS)
model as implemented in the finite element code P LAXIS V8. The small-strain enhanced
version of the HS model resulting from the combination is referred to as the HS-Small
model. Compared to the HS model, the HS-Small model additionally offers a MatsuokaNakai [113] [114] failure criterion (optional) and a refined flow rule (mandatory). Following a brief introduction to the basic equations of infinitesimal elastoplasticity, both
models constitutive equations are developed within this chapter.
Because of the availability of two yield criteria, it is generally distinguished between
the HS-Small(MC) model with Mohr-Coulomb failure criterion, and the HS-Small(MN)
with Matsuoka-Nakai failure criterion. The terminology HS-Small is used whenever the
actual type of yield criterion is not of interest for a particular discussion. The term HSSmall is thus a collective term for HS-Small(MC) and HS-Small(MN).
Within this chapter, all stresses are effective although not explicitly indicated by a
prime. In triaxial conditions, 1 = axial , denotes axial strain and 2 = 3 = lateral , denotes
lateral strain. In the following, it is generally assumed axial /lateral 1. The reader can
easily verify all triaxial equations to hold likewise in triaxial extension (axial /lateral < 1)
by substituting 1 with 3 and vice versa.
(7.1)
(7.2)
(7.3)
(7.4)
where ij , eij , and pij denote the total, elastic, and plastic strain tensor respectively, ij
is the Cauchy stress tensor, q is a suitable set of scalar or tensorial internal material
75
(7.5)
(7.6)
(7.7)
f
D dkl
ij ijkl
f
g
D
+
ij ijkl kl
with h =
f
h .
q
(7.9)
If all derivatives in Equation 7.9 are constant or linear functions in stress and the internal
variables, the plastic multiplier may be found analytically. However, commonly used
yield functions and plastic potentials do not satisfy these requirements. Calculation of
the plastic multiplier and at the same time the integration of the constitutive equations
therefore requires an efficient iterative scheme.
Many constitutive laws can be formulated within the framework of elastoplasticity
presented above. These laws differ only in their yield function, plastic potential, and
evolution laws for the internal variables. When functions of strain are applied the above
formulation can also be recast in strain space in terms of { ij , pij }, by eliminating ij
with the help of Equation 7.1 and Equation 7.2.
The implications of large strains are not considered explicitly in this thesis. However,
it should be noted that conserving objectivity, for example by the consideration of a
Jaumann term in the stress-rate, will not alter the following model description.
76
q
qa
asymptote
qa
qf
asymptote
1
E50
E50
qf
qa
Eur
1
e1
e1
Figure 7.1: Hyperbolic stress-strain law by Kondner & Zelasko (left) and its modification
after Duncan & Chang (right).
q
2 sin
q50
qa
with qa =
(3 + c cot ) and 50 =
=
qa q
1 sin
E50
2E50
(7.10)
Duncan & Chang [36] based their hypoelastic model on the above formulation by Kondner & Zelasko, additionally introducing the deviatoric measure qf in the form:
1 = 50
q
2 sin
qf
f or q < qf =
(3 + c cot ) and qa =
qf
qa q
1 sin
Rf
(7.11)
The conceptual difference in the formulations by Kondner & Zelasko and that by Duncan
& Chang is illustrated in Figure 7.1. Extending the hypoelastic Duncan-Chang model to
an elastoplastic formulation, Schanz [157] proposed the following yield function:
fs =
qa
q
2q
ps
E50 qa q Eur
(7.12)
where ps is an internal material variable for the accumulated plastic deviatoric strain,
q = 1 3 is defined for triaxial loading, and qa is the asymptotic deviatoric strain
77
re
1200
om
b fa
ilu
1000
300
600
400
200
p
on
g =c
st.
200
100
0
-Co
ul
800
Mo
hr
400
q (kN/m)
q (kN/m)
q (kN/m)
500
100
200
300
400
500
p (kN/m)
200
400
600
800
p (kN/m)
p (kN/m)
Figure 7.2: Experimentally derived shear strain contours versus HS yield loci. Left:
Shear strain contours of loose Fuji River sand [71]. Middle: Shear strain contours of dense Fuji River sand [178]. Right: Yield loci of the HS model for
stress independent stiffness (m=0).
(7.13)
q
1
qa
q
q50 q
q
fs =
=
= 50
.
Eur 2
2E50 qa q
E50 qa q
qa q
(7.14)
In accordance with the basic idea of Vermeer [185], the defined yield loci are lines
of constant plastic shear strain in p-q space. The striking similarity between experimentally derived shear strain contours and the yield loci defined by Equation 7.12 is depicted
in Figure 7.2. All the shear strain contours shown in Figure 7.2 (left and middle) were
experimentally derived by Tatsuoka & Ishihara [178], and Ishihara et al. [71]. Their
triaxial testing procedure was as follows: They first sheared a given sample at mean
stress p1 to shear stress q1 . After reducing shear stress (unloading) and increasing mean
stress to p2 , they sheared the sample again until it yielded at shear stress q2 . This procedure was repeated many times so that the regression line through the n-tuple of points
{(p1 , q1 ), (p2 , q2 ), , (pn , qn )}, called the shear strain contour, marks a line of constant
shear strain in p-q space.
When translating Equation 7.12 to principal stress space, the resulting cone-type yield
78
2qa (1 2 )
Ei qa (1 2 )
2qa (1 3 )
Ei qa (1 3 )
2(1 2 )
Eur
2(1 3 )
Eur
ps , and
ps .
(7.15)
The transition from E50 to Ei is made because of a second yield surface that will be introduced later in this section. The second yield surface will affect material stiffness such
that, the meaning of Ei in the full HS model is not as closely related to the hyperbolic
model by Kondner & Zelasko as the one of E50 .
Associated plasticity is an unrealistic assumption for most geotechnical materials.
Thus, an additional plastic potential of the form:
s
g12
=
s
g13 =
(1 2 )
2
(1 3 )
2
1 +2
2
1 +3
2
sin m , and
sin m
(7.16)
is introduced. The mobilized dilatancy angle m in the above equation is defined according to Rowes stress dilatancy theory [153]:
sin m =
sin m sin cs
1 sin m sin cs
(7.17)
where, cs is the critical state friction angle and the mobilized friction angle m is calculated as follows:
1 3
sin m =
.
(7.18)
1 + 3 + 2 c cot
The stiffness moduli Ei and Eur are scaled for their stress dependency with an Ohde
[131] or Janbu [77] power law:
cot m
) , and
Ei = Eiref ( pref3 +c
+c cot
ref 3 +c cot m
Eur = Eur ( pref +c cot )
(7.19)
where Eiref and Eurref are the material stiffness moduli at the reference pressure pref , and m
is the exponent of the power law. In Equation 7.19 the minor principal stress 3 is used
as an indicator of the actual stress state in the material instead of the mean stress p = 3ii .
The cone-type yield loci of the HS model mainly accounts for plastic deviatoric strains,
or shear hardening. A second cap-type yield surface is introduced next. The cap-type
yield surface accounts for plastic volumetric strains, or volumetric hardening. Volumetric hardening corrects too stiff primary oedometric or isotropic loading, obtained in pure
shear hardening models. The HS models cap-type yield surface is defined as follows:
q2
p2 p2p
(7.20)
2
where p = 3ii is the mean stress, is an internal material constant, controlling the steepness of the cap in p-q space as shown in Figure 7.4, pp is an internal material variable for
pre-consolidation stress, and q is a special stress measure, defined as:
fc =
q = 1 + ( 1 1)2 1 3 with =
3 sin
.
3 + sin
(7.21)
79
s1
q
s2
s3
s2
s3
Figure 7.3: Yield surfaces of the HS model for cohesionless soil. Left: Cap and cone-type
yield surfaces in principal stress space with the cone being in its ultimate
Mohr-Coulomb failure criterion position. Right: p-q slice through the yield
surfaces.
The definition of the special stress measure q is necessary to adopt the cap-type yield
surfaces deviatoric shape to this of the cone-type yield surface as shown in Figure 7.3.
The plastic potential for the cap-type yield surface is chosen equal to its yield surface
(g c = f c ), so that plastic strain on the cap-type yield surface is associated, in contrast to
plastic strain, on the cone-type yield surface.
A total of two state variables are used within the HS models yield and potential functions: The plastic deviatoric strain ps , and the pre-consolidation stress pp . Their evolution laws are defined as follows:
d ps = ds h ps
dpp
= dc hpp
g
g
g
)=1
with h ps = (
1
2
3
m
3 +c cot
with hpp = 2H pref +c cot
p
(7.22)
where m represents the power law exponent, and H relates plastic volumetric strain
pv = p1 + p2 + p3 to pre-consolidation stress pp as follows:
m
3 + c cot
dpp = H
dpv .
(7.23)
pref + c cot
In decomposing volumetric strain in elastic and plastic contributions, H can be rewritten
as a function of the bulk stiffness in unloading-reloading, Ks and the bulk stiffness in
primary loading, Kc as:
1
Ks Kc
= Ks
Ks .
(7.24)
H=
Ks Kc
1
Kc
where due to the assumption of isotropic elasticity, the elastic bulk stiffness Ks relates to
Eurref as follows:
Eurref
Ks =
.
(7.25)
3(1 2)
80
en
str
h
gt
ite
cr
elastic region
c cotj
ria
pp
Figure 7.4: Evolution of the HS modelss cap and cone-type yield surfaces in p-q space.
The model parameter H can therefore be determined by the bulk stiffness ratio Ks /Kc .
As the physical meaning of the latter is more evident, it is often used to quantify H. An
illustration of yield surface evolution in p q space is given in Figure 7.4. Note that by
the above definition of d ps = ds , it is always assured that for zero volumetric strain,
d ps = 21 as required in Equation 7.13.
In order to formally complete the model description in the framework of infinitesimal
elastoplasticity, the elastic stiffness tensor Dijkl must be defined. The HS model assumes
isotropic elasticity inside the yield loci, so that Dijkl is assembled from the user defined
stiffness Eur and Poisons ratio, ur as follows:
Dijkl =
Eur
((1 2ur )ik jl + ur jk jl ).
(1 + ur )(1 2ur )
(7.26)
81
Eur
((1
(1+ur )(12ur )
2ur )ik jl + ur jk jl )
2. Yield functions:
(1 2 )
s
a
f12
= 2q
Ei qa (1 2 )
s
f13
=
fc
2(1 2 )
Eur
2(1 3 )
Eur
ps
2qa (1 3 )
Ei qa (1 3 )
q2
2
2
p pp where
2
1 + ( 1 1)2 1 3
ps and
with =
3sin
3+sin
3. Plastic potentials:
2)
s
g12
= (1
2
s
g13
=
gc
1 +2
sin m
2
(1 3 )
1 +3
2 sin m and
2
q2
2
p p2p where
2
1 + ( 1 1)2 1 3 with
3sin
3+sin
sin m sin cs
1sin m sin cs
sin sin
1sin sin
5. Hardening laws:
d ps = ds h ps
dpp
= dc hpp
Ks Kc
Ks Kc
0, where
g
g
g
with h ps = (
)=1
1
2
3
m
3 +c cot
p, where
with hpp = 2H pref +c cot
and Ks
ref
Eur
3(12)
82
s1
q
s2
s2
s3
s3
Figure 7.5: Yield surfaces of the HS-Small model for cohesionless soil. Left: Cap- and
cone-type yield surfaces in principal stress space with the cone being in its
ultimate Matsuoka-Nakai failure criterion position. Right: p-q slice through
the yield surfaces.
E
2(1+)
E
3(12)
(7.27)
hold. Two additional material parameters are needed to control the stress and strain history dependent stiffness of the HS-Small model. These are the initial shear modulus Gref
0
defined for the reference pressure pref and the shear strain 0.7 , at which the shear modulus has decayed to 70 percent of its initial value. The shear modulus G0 is calculated
from:
m
3 + c cot
ref
G0 = G0
,
(7.28)
pref + c cot
where m is the power law exponent that scales also the other stiffness parameters of the
HS-Small model (see Equation 7.19). The threshold shear strain 0.7 is taken independent
of mean stress.
A lower cut-off in the hyperbolic small-strain stiffness reduction curve is introduced at
the shear strain c where tangent stiffness is reduced to the unloading-reloading stiffness
Gur in larger strain cycles. The unloading-reloading shear modulus Gur relates to the HS
model parameter Eur as follows:
Gur =
Eur
.
2(1 + ur )
(7.29)
83
g0.7
HS-Small
Hardin & Drnevich
Gur
gc
Figure 7.6: Cut-off in the tangent stiffness degradation curve as used in the HS-Small
model.
The shear strain c is then obtained from the modified Hardin-Drnevich relationship,
which for the tangent stiffness Gur yields:
2
ac
0.7
1
Gur = G0
= G0
(7.30)
2 +
c
1 a 0.7
0.7 + ac
c
0.7 1 a 0.7
and hence:
0.7
c =
a
!
G0
1 .
Gur
(7.31)
The resulting small-strain stiffness behavior of the HS-Small model is illustrated in Figure 7.6. Its isotropic elasticity tensor is calculated as:
Dijkl
2G
=
((1 2ur )ik jl + ur jk jl ) with G = G0
1 2ur
0.7
0.7 + aHist
2
.
(7.32)
The stiffness tensor given by Equation 7.32 is not a truly elastic, but a paraelastic one
(Hueckel & Nova [65]). Nevertheless, the terminology elastic stiffness is maintained for
the use of the above defined stiffness in the HS-Small model as a reminder that it applies
to all stress states, within as well as on the yield surface, similar to the elastic stiffness
known from traditional elastoplastic models.
By default, the HS model is hardened under initial loading. Without altering its hardening laws, plastic straining will reduce the HS-Small models small-strain stiffness. In
order to preserve the stiffness defined by the material parameters G0 and 0.7 , the HSSmall hardening laws are rewritten as:
d ps = ds hi h ps
dpp = dc hi hpp
84
with h ps
with hpp
(7.33)
hi := Gm
(7.34)
where Gm is the stiffness multiplier defined in Section 6.1.3. Note, that in combination
with the HS model, the threshold shear strain 0.7 in the Small-Strain Overlay model
is always doubled for reloading. If a primary loading condition exists, this will be accounted for by plastic hardening of the yield locus. For further details on the distinction
between initial and reloading, it is referred to Section 6.1.3.
sin m sin cs
0.
1 sin m sin cs
(7.35)
Unfortunately, the dilatancy cut off introduced in Equation 7.17 typically yields too little
plastic volumetric contraction. Recently, several formulations have been proposed to
resolve this problem.
Wehnert [191] proposed a non-zero lower cut-off value, similar to Equation 7.35, for
the allowable contraction in loose sands:
sin m =
sin m sin cs
0.07 .
1 sin m sin cs
(7.36)
Even though the method by Wehnert is efficient, the lower cut-off value effectively represents an additional material parameter. Soreide [169] proposed to scale Rowes equation
as a function of mobilized friction:
sin m sin cs
sin m
sin m =
.
(7.37)
1 sin m sin cs
sin
From a quantitative point of view, the methods by Wehnert, and Soreide seem to be very
much alike as shown in Figure 7.7 and again in the model validation (Chapter 8). Soreide
additionally claims numerical benefits from his formulation: The likelihood of stress
return towards the apex is less when using big step sizes in the return mapping (see
Section 7.4 for further details on return mapping). However, having well defined initial
stresses and a suitable return mapping scheme this advantage should not be overrated.
In contrast to the Rowe trouble shooting strategies by Wehnert and Soreide, the modified flow rule of the HS-Small model makes use of critical state soil mechanics. Li
& Dafalias [103] proposed to treat dilatancy as a state-dependent quantity within the
85
d0
(M em ),
M
(7.38)
where d0 , and m are two material parameters, M is the critical stress ratio, = pq is the
actual stress ratio, and is the state parameter defined by Been & Jefferies [13]. The state
parameter is a scalar quantity that measures the difference between the actual and
the critical void ratio at the same mean stress p. Equation 7.38 reduces to the Cam-Clay
dilatancy (Roscoe & Schofield [151]) by specifying m = 0, and d0 = M .
Using the original formulation by Li & Dafalias in the HS-Small model would require
four additional material input parameters. These are the parameters d0 , and m introduced in Equation 7.38, and two parameters which define the critical state line in e ln p
space. The critical state line can for example, be defined by an initial void ratio e0 and
its slope . Unfortunately, the parameters d0 , and m can be obtained in numerical fitting
procedures only. Altogether, the input to the original approach by Li & Dafalias seems
too demanding for a model that is meant to be used in practical engineering applications. Therefore, the formulation by Li & Dafalias is simplified by a) calculating the state
parameter, based on constant void ratio, and b) providing default parameters for m,
and d0 . With these simplifications, Equation 7.38 can be rewritten as:
1
where
(7.39)
pcs
sin cs (1 sin m )
=
p
M sin m (1 sin cs )
(7.40)
1
, and d0 = M
. These default parameters are in good agreement to the
and m = 15
10
ones presented in Li & Dafalias [103]. As the default parameters cannot be guaranteed to
work in dilatant conditions, they are used for contraction only. Whenever dilatant material behavior is detected, the HS-Small model falls back to the original Rowe equation.
Figure 7.7 compares the mobilized dilatancy angles of all flow rules discussed above.
The comparison of mobilized dilatancy angles is supplemented by test data in Chapter 8.
However, it can already be noted that not one of the discussed approaches is considered
fully satisfactory. In future work, the HS-Small model could perhaps fully adopt the
approach by Li & Dafalias, therefore extending its use to dilatant behavior in a void
ratio dependent formulation (see Chapter 10).
86
j = 35, y = 10
10
ym []
ym []
j = 35, y = 0
10
-10
-20
Rowe
-10
-20
Rowe
Wehnert
Soreide
-30
Wehnert
Soreide
-30
HS-Small
0
10
20
jm []
30
HS-Small
0
10
20
30
jm []
Figure 7.7: Mobilized dilatancy angle versus mobilized friction angle. Results from the
original Rowe equation, its modification by Wehnert, and Soreide, and the
modified Li & Dafalias formulation.
between triaxial compression and extension. Figure 7.8, for example, shows experimental results for reconstituted Osaka alluvial clay. Similar findings for sand and sandstone
have been for example reported by Lade [97], Ramamurthy & Rawat [144], Sutherland
& Mesdary [176], and Yamada & Ishihara [196].
Matsuoka & Nakai [113], [114] (MN) proposed a failure criterion that is in better agreement with the experimental surface shown in Figure 7.8. They propose the concept of
a Spatial Mobilized Plane (SMP), which defines the plane of maximum spatial, averaged
particle mobilization in principal stress space. The SMP is geometrically constructed
by deriving the mobilized (Mohr-Coulomb) friction angles for each principal stress pair
separately (Figure 7.9, left) and sketching the respective mobilized planes in principal
stress space (Figure 7.9, right). Matsuoka & Nakai derive their failure criterion by limiting the averaged ratio of spatial normal stress to averaged spatial shear stress on this
plane. Their failure stress ratio can be expressed as a simple function of the first, second,
and third stress invariant, I1 , I2 , andI3 as shown in Equation 7.41. With the SMP concept, the Matsuoka-Nakai criterion automatically retains the well established material
strength of the Mohr-Coulomb criterion in triaxial compression and extension.
The likewise well-known failure criterion by Lade & Duncan [99] (Lade) appears compared to the Mohr-Coulomb criterion and the Matsuoka-Nakai criterion rather optimistic in plane strain conditions and triaxial extension. Therefore, it has not (unlike
the Matsuoka-Nakai criterion), been considered as an option for the HS-Small model,
though it is mentioned here for the sake of completeness. It should also be mentioned
that using bifurcation analysis, progressive failures would most likely correct for the
87
s1
Experimental surface
Mohr-Coulomb surface
(j = 33.7, c = 2 kPa)
s3
s2
Legend:
sc = 100 kPa
sc = 200 kPa
sc = 300 kPa
Figure 7.8: Test results by Shibata & Karube [163] on reconstituted Osaka alluvial clay.
sI
s1
SMP
Mobilized
Planes
fmob23
fmob13
45+
fmob12
s3
s2
s1
45+
fmob12
2
s3
45+
s2
fmob13
2
sIII
fmob23
2
sII
Figure 7.9: The SMP concept of Matsuoka & Nakai. Left: Three mobilized planes where
the maximum shear stress to normal stress ratio is reached for the respective
principal stresses. Right: SMP in principal stress space.
88
I13
c1
I3
I1 I2
c2
I3
=0
with
c1 =
=0
with
c2 =
(3+sin )3
(1sin )(1+sin )2
9sin2
1+sin2
I1 = ii
I2 = 12 (ij ij ii jj )
I3 = 61 (ii jj kk + 2ij jk ki 3ij ji kk ).
(7.41)
(7.42)
(7.43)
The constants c1 and c2 in Equation 7.41 are defined so that both failure criteria are
identical to the Mohr-Coulomb criterion in triaxial compression.
The HS-Small model optionally uses a Matsuoka-Nakai yield surface and criterion instead of the original HS models Mohr-Coulomb criterion. In the HS-Small(MN) model
the non-associated Mohr-Coulomb potential is substituted by a Drucker-Prager [35] potential. The cap-type yield function remains unchanged, only its potential is adjusted to
the new deviatoric flow of the cone-type yield surface.
The true triaxial experiments by Yamada & Ishihara [196] shown in Figure 7.10 not
only support the failure hypothesis by Matsuoka & Nakai but at the same time proof
that the Drucker-Prager potentials radial deviatoric flow direction is a reasonable assumption for low mobilized friction angles. For higher mobilized friction angles, flow
directions increasingly tend towards an associated deviatoric flow, the solid curve in
Figure 7.10 represent the Matsuoka-Nakai failure criterion ( = 30 ). However, in the
example presented in Figure 7.10, the error in assuming radial flow for high mobilized
friction angles is of the same magnitude or less than that made in assuming deviatoric
associated flow. A more detailed discussion of the test results by Yamada & Ishihara is
89
Direction of strain
increment
Q=-15
Failure stress
Q=-30
Q=-15
Q=0
Q=15
s3
s2
Q=15
Q=-30 Q=-15
Q=0
Figure 7.10: Experimentally derived shear strain increment vectors projected on the octahedral plane (after Nakai & Matsuoka [128]).
HS-Small(MN) model or its return mapping. As the Lode dependent formulation of the
original Matsuoka-Nakai yield criterion is bulky, several other simplified Lode dependent yield criterion were also proposed in the literature.
The Lode angle [106] in this thesis is defined as:
=
1
3
arcsin( 3
3J3
3/2 )
2J2
(7.44)
where J2 and J3 are the second and third deviatoric stress invariants. These are defined
for the deviatoric stress tensor, ij = ij 31 ii analog to Equation 7.42, substituting by
. In principal stress space, J2 and J3 can be expressed as:
J2 =
J3 =
1
((1 2 )2 + (2 3 )2 + (3 1 )2 )
6
1
(21 2 3 )(22 3 1 )(23
27
1 2 ).
(7.45)
90
=
with
and
1
3
2 2 +1 cos
(
1
6
arccos(1 +
= 3 16
3sin
3+sin
arccos(1
27 2 (1)2
sin2 (3))
2( 2 +1)
2
(1)2
+ 27
sin2 (3))
2
2( +1)
for 0
for > 0
(7.46)
Unfortunately, the derivatives of the LMN dependence are bulky compared to those
of other Lode dependent strength criteria proposed in the literature. Here the criteria
by Argyris [6], Papadimitriou & Bouckovalas (PB) [133], and Wang [102] shall be mentioned:
2
(7.47)
(1 + ) (1 ) sin(3)
4
1
PB: () =
((1 + ) + (1 ) sin(3)) (7.48)
(1 + ) (1 ) sin(3) 2
p
(1 + 2 )2 + 4(1 2 ) sin(3) (1 + 2 )
Wang: () =
(7.49)
2(1 ) sin(3)
Argyris: () =
Results of all the above strength criteria are plotted in Figures 7.11, 7.12, and 7.13 for
triaxial compression friction angles of 20 , 30 , and 40 respectively. Obviously the
less complex formulations 7.47 to 7.49 tradeoff convexity at high friction angles for their
simplicity. The Matsuoka-Nakai (MN), and Lade yield surfaces stay convex even at high
friction angles and hence, fulfill Druckers convexity postulate [33].
Due to the Drucker-Prager plastic potential of the HS-Small(MN) model, the Lode
angle is constant during return mapping (see Section 7.4). This favors the use of the
LMN dependence in the HS-Small(MN) model. The less complex and very similar criterions by Wang or Papadimitriou & Bouckovalas would seem more appropriate in case
the condition
= 0, does not hold while higher order derivatives are needed in the
i
algorithmic setting.
1 3
1 + 3 + 2c cot
and
qf =
q = sin m (1 + 3 + 2c cot )
2 sin
(3 + c cot )
1 sin
(7.50)
(7.51)
91
1.2
30
s1
MN
Lade
0.8
Argyris
28
PB
0.4
Wang
j []
26
0
24
-0.4
s3
s2
22
-0.8
-1.2
-1.2
20
-0.8
-0.4
0.4
0.8
1.2
0.2
0.4
0.6
0.8
b = (s2-s3)/(s2-s1) [-]
deviatoric plane (p = 1)
Figure 7.11: Different failure criteria for a friction angle of = 20 in triaxial compression. Left: Deviatoric shape of the yield criteria. Right: Friction angle as a
function of the intermediate stress parameter b.
1.2
40
s1
MN
Lade
0.8
Argyris
38
PB
0.4
Wang
j []
36
0
34
-0.4
s3
s2
-0.8
-1.2
-1.2
32
30
-0.8
-0.4
0.4
deviatoric plane (p = 1)
0.8
1.2
0.2
0.4
0.6
0.8
b = (s2-s3)/(s2-s1) [-]
Figure 7.12: Different failure criteria for a friction angle of = 30 in triaxial compression. Left: Deviatoric shape of the yield criteria. Right: Friction angle as a
function of the intermediate stress parameter b.
92
50
s1
MN
Lade
0.8
Argyris
48
PB
0.4
Wang
j []
46
0
44
-0.4
s3
s2
42
-0.8
-1.2
-1.2
40
-0.8
-0.4
0.4
0.8
1.2
0.2
0.4
0.6
0.8
b = (s2-s3)/(s2-s1) [-]
deviatoric plane (p = 1)
Figure 7.13: Different failure criteria for a friction angle of = 40 in triaxial compression. Left: Deviatoric shape of the yield criteria. Right: Friction angle as a
function of the intermediate stress parameter b.
and thus:
q
= Rf
qa
1 sin
sin
sin m
1 sin m
(7.52)
(7.53)
which reduces to s = 1 3 in triaxial conditions, one obtains s = 23 1 for zero volumetric strain. Adopting the original HS models yield function (Equation 7.12) to shear
strain s and considering Equation 7.52 results in:
1sin m
sin m
3 q
3 q
sps
(7.54)
fs =
1sin
1sin
m
2 Ei
2 Eur
R
sin m
sin
where sin m is the mobilized friction angle in triaxial compression. In mobilized friction,
the Matsuoka-Nakai yield criterion defined in Equation 7.41, can be written as:
sin2 m
9
1
I1 I2
I3
I1 I2
I3
(7.55)
This definition yields deviatoric isolines of mobilized friction that are similar to the shape
of the Matsuoka-Nakai yield criterion. Alternatively, Equation 7.55 can be expressed in
93
3q
.
6(p + c cot ) + q
(7.56)
6 sin m
3 sin m
(7.57)
Due to the cone-type yield surfaces new shape, the cap-type yield surface is also to
be reformulated. The specialstress measure defined in Equation 7.21 is now replaced
by the Roscoe invariant q = 3J2 . At the same time the Lode angle dependency of the
cone-type yield surface is translated to the cap by scaling its steepness:
q2
f =
p2 p2p .
2
()
c
(7.58)
q2
p2 p2p
2
()
(7.59)
where = (T rial ). The caps deviatoric plastic flow direction is thus consistent with
the cones radial Drucker-Prager potential.
Although the HS-Small model should respond slightly different under general (plastic) loading conditions, its triaxial behavior should be equal to the HS model. The HSSmall evolution laws are therefore adapted to the objective shear strain measure as follows:
q
g 2
g
g 2
g
g 2
1
(( g
) + (
) + (
) ) = 32
dsps = ds hi hsps with hsps =
2 1
2
2
3
3
1
m
cot
dpp = dc hi hpp
with hpp = 2H pref3 +c
p
+c cot
(7.60)
The evolution laws finally round up the formulation of the HS-Small model, which is
summarized for the HS-Small(MN) model in Box 7.2. A complete list of model parameters can be found in Table 7.1. The integration of the constitutive equations is discussed
in the next section.
94
2
G
0.7
for s < c
0 0.7 +aHist
G
=
E
ur
for s c
2(1+ur )
2. Yield functions:
fs
= 32 Eqi
m
( 1sin
sin m )
(
)Rf ( 1sin
sin )
q2
2
2
= ()
2 p pp where
fc
sin2 m =
1sin m
sin m
3 q
2 Eur
sps
I1 I2
I3
I1 I2
1 I3
3. Plastic potentials:
6 sin m
g s = q (p + c cot ) 3sin
m
q2
()
2
gc =
p2 p2p where
= (T rial )
4. Mobilized
(dilatancy (modified Rowe):
sin m sin cs
for sin m sin cs 0
1sin m sin cs
sin m =
1
1
ln(p
/p)
cs
(M e 15
) for sin m sin cs < 0
10
sin cs =
sin sin
1sin sin
5. Hardening laws:
dsps = ds hsps
dpp
H
= dc hpp
=
Ks Kc
Ks Kc
and
pcs
p
sin cs (1sin m )
M sin m (1sin cs )
with hsps =
with hpp
and Ks
3
2
= 2H
=
3 +c cot
pref +c cot
ref
Eur
3(12)
7. Lode Matsuoka-Nakai
dependency:
3
1
() = 22 +1
with
cos
(
2 (1)2
1
arccos(1 + 27
sin2 (3))
for 0
6
2( 2 +1)
() =
2
2
(1)
16 arccos(1 + 27
sin2 (3)) for > 0
3
2( 2 +1)
3sin
3+sin
95
96
Symbol
ref
E50
ref
Eoed
Eurref
m
c
ur
pref
K0nc
Rf
Tension
E0ref
0.7
Eiref
Ks /Kc
ps
sps
pp
Hij
Unit
kN
m2
kN
m2
kN
m2
]
[ kN
m2
[ ]
[ ]
[ ]
kN
m2
[ ]
[ ]
kN
m2
kN
m2
[ ]
kN
m2
[ ]
[ ]
[ ]
[ ]
kN
m2
[ ]
HS
HS-Small
X
X
X
X
X
X
X
X
X
X
X
X
-
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
-
X
X
X
f
n+1
nij =
ij n+1
g
n+1
mij =
(7.61)
ij n+1
f
n+1
= q
n+1
where again the asterisk denotes n indices or vanishes for scalar functions.
n+1
ij = Dijkl n+1kl n+1 pkl
(7.62)
n+1 p
n p
n+
kl = kl + mij
(7.63)
n+1
n
n+
q = q + h
(7.64)
Fn+1 = 0
(7.65)
where
n+
h = h (1 ) nij + n+1ij , (1 ) nq + n+1q
n+
and
(7.66)
(7.67)
and the superscripts n and n + 1 denote two subsequent steps in the integration procedure and 0 1, defines the time when the flow rule is evaluated. For = 0 the flow
rule is evaluated at the initial stress or, if this is elastic, at the intersection with the yield
97
n+1
cross
n+a
s
n+1
n+1
n+1
n+a
n+a
sTrial
n
m
f
s33
s22
Figure 7.14: Deviatoric return mapping using the generalized midpoint algorithm with
different values. The chosen yield surface (Matsuoka-Nakai) and plastic
potential (Drucker-Prager) correspond to the HS-Small(MN) model.
ij =n ij + Dijkl
98
n+1
dkl .
(7.68)
mkl
mkl
n+1
dmn + Dijkl
dq = 0.
(7.71)
rij + dij + Dijkl
mkl + Dijkl
mn n+1
q n+1
By algebraic manipulation the above equation can be solved to dij :
mkl
ijmn =
im nj + Dijkl
and
mn n+1
mkl
n+1
kl =
mkl +
h
q n+1
(7.72)
(7.73)
(7.74)
(7.75)
which is derived by a first order Taylor expansion around the final stress state n+1f = 0.
Solving the above equation for d, the solution for the iterative update of the plastic
multiplier is readily found as:
n+1
d =
(7.76)
where again the conventions defined in Equation 7.61 have been used.
In order to successfully apply the closest point projection as an iterative procedure, a
suitable starting point is necessary. A good guess for the starting point can typically be
found by using the normal to the elastic trial point [83]:
T rial
(7.77)
and assuming that any stress relaxation from the trial stress is due to plastic straining
and thus:
(7.78)
dT rialmn = Dmnpq d(0) T rialmpq .
99
T rial
T rialn
T rialm
nm Dnmpq
pq
T rial h
(7.79)
(0)
T rial
mpq .
(7.80)
The algorithmic setting for the closest point projection algorithmic can now be summarized as follows (Box 7.3):
f
f
dn+1f =
dq
d
+
ij
q
ij n+1
(7.81)
n+1
= n+1nij dij +n+1 dq = 0.
The unknown stress increment dij is found by linearization or differentiation of the
implicit Backward Euler equation (7.69) to:
mkl
mkl
n+1
dij = Dijkl dkl dDijkl
mkl Dijkl
dmn Dijkl
dq (7.82)
mn n+1
q n+1
and isolation through simple algebraic manipulations:
dmn = ijmn Dijkl (dkl dkl )
(7.83)
where and are defined in Equation 7.73 and 7.74 respectively. Finally substituting
7.83 into 7.81 gives an expression for the consistent elastoplastic tangent:
ep
dmn
dpq = Dpqmn
ep
Dpqmn
= (rspq Drsmn )
100
with
(rspq Drskl )kl n+1nij (rsij Drsmn )
n+1n (
n+1 h
ot
rsot Drspq )pq +
(7.84)
(7.85)
ij =n ij + Dijkl
T rial (0)
n+1
dkl
(0)
f (T rialij ,T rial q (0) )
d(0) =
T rialn
(k)
(k)
(k)
mkl
mn
n+1
, and
(k)
mkl
q
n+1
(k+1)
(k)
n+1f r
n+1n
mn
ij ijmn
(k+1)
(k)
d
= d + Dijkl kl ijmn n+1nmn n+1 h
d) Update stress and state variables
(k)
n+1 (k+1)
mn =T rial mn d(k+1) Dmnpq n+1mpq
n+1 (k+1)
q
=n q + d(k+1) h
e) SET k = k + 1
END DO
101
corner region
n+1 s
n ij
n+1
msij
sTrial
ij
n+1 c n+1
n ij = mcij
n+1
sij
s
c cotj
sij
c
f =g
elastic region pp
Figure 7.15: Corner problem at the intersection of cone- and cap-type yield surface.
When using an implicit return mapping scheme as the one presented in the previous
section, it is only a small effort to calculate the consistent elastoplastic tangent. However, the main drawback associated with the implicit closest point projection algorithm
and the calculation of the consistent tangent is their need for the plastic potentials second order gradients. This may prove exceedingly laborious for more complex plastic
potentials, but is straightforward in the HS-Small model.
102
q
corner region
f
msij
n+1
n+1
mtij
apex gray
region
sTrial
ij
f =g
n
sij
n+1
sij
apex c cotj st
elastic region pp
Figure 7.16: Apex gray region and tension cut-off illustrated in a section with 2 = 0.
If not smoothly curved, the apex is another singularity in the yield surface. The apex
gray region is defined by the gradient to the cone-type potential surface as shown in Figure 7.16. If the apex corresponds to an admissible tensile stress, the trial stress is returned
to the apex. If the apex point violates the user defined maximum allowable tensile stress,
a two vector return mapping scheme to the respective tension cut-off planes is evoked.
The HS-Small model uses a tension cut-off criterion based on minimum principal stress,
which implies three (fixed) orthogonal tension cut-off planes in principal stress space:
fit = Tension i ,
(7.87)
103
Chapter 8
Validation and verification of the HS-Small
model
The HS-Small model is validated by comparing numerically simulated experimental
data to real (measured) data. Model validation is divided into two main parts: Element
tests and boundary value problems.
The element tests section strictly speaking also contains partly code verification. Iteratively derived stiffness decay curves are checked against the closed form model by
Hardin & Drnevich in the version of Santos & Correia. This is necessary due to the hardening plasticitys impact on total elastoplastic stiffness. For clarification of the differences
implied by the terms validation and verification see also Figure 1.2.
All element tests are calculated on the Gau-Point level. All boundary value problems
are analyzed using the finite element code P LAXIS V 8.2 [21]. As the HS-Small model was
implemented in the calculation kernel of this code, standard interface definitions and
calculation settings could be used. The help and cooperation of Plaxis B.V, especially the
help of Dr. Paul Bonnier, during model implementation is greatly acknowledged here.
All calculations are performed with triangular elements (6 nodes, 3 integration points)
and default settings in the iteration controls of the finite element code P LAXIS.
Drome,
Rhone-Alpes,
France. Hostun sand results from the specification of its sieving
process (Flavigny et al. [41]) uniformly graded. Its coefficient of uniformity is U =
D60
= 1.70, where its medium grain size is D50 = 0.35mm. The minimum and maximum
D10
void ratio are experimentally derived as emin = 0.63 and emax = 1.00. The experimental
database distinguishes between dense and loose Hostun sand, based on its initial void
ratio e. For initial void ratios in the range of e = 0.63...0.68 the sand is classified as
dense, whereas in the range of e = 0.85...0.92, it is classified as loose. The corresponding
105
106
s1[kN/m2]
400
300
s1[kN/m2]
200
HS-Small
HS (original)
150
Experiments
200
100
100
50
0
0.000
0.004
0.008
0.012
e1[-]
0
0.000
0.005
0.010
0.015
0.020
e1[-]
Figure 8.1: Oedometer test on dense (left) and loose (right) Hostun sand.
Strain Overlay) yields too low small-strain stiffness. Here, plasticity is over-dominating
the initial elastoplastic stiffness.
It can be noted that there is a small numerical discrepancy between the material input
parameter 0.7 and the computed shear strain (HS-Small), at which small-strain stiffness
is decreased to 70% of its initial value. For practical applications this discrepancy can be
easily neglected since it is probably well below the experimental precision and repeatability. Alternatively, the internal model parameter hi (see Section 7.3.1), could be fitted in
an iterative procedure by including it with the internal parameter conversion discussed
in Chapter 7. In either case, a perfect match between the HS-Small decay curves and the
analytically derived curves is obtained for zero plastic straining. Then, the elastoplastic
model degenerates to the Small-Strain Overlay model as presented in Chapter 6.
Altogether, the simulated triaxial tests match the test data presented, reasonably well.
It turns out that the modified flow rule has a vanishing effect on the volumetric behavior
in drained tests. This will be proven to be different for undrained tests later in this
section.
The implementation of the Small-Strain Overlay model into the HS-Small model does
not overly affect the overall stress-strain curve. There is no over-pronounced shift or
bump in the stress-strain curve. Combining the Small-Strain Overlay model with a bubble model, or any other model that has an initial elastic range would probably not give
such smooth results.
In oedometer and triaxial tests the intermediate stress parameter b always equals 0.
Both failure criteria used in the HS-Small model, Mohr-Coulomb (MC) and MatsuokaNakai (MN), are coincident in this case. For this reason, there has so far been no distinction made towards the HS-Small models failure criteria in the validation process.
In biaxial tests (plane strain), however, the non-zero intermediate stress parameter increases the material strength estimate of the Matsuoka-Nakai criterion compared to that
of the Mohr-Coulomb criterion.
Figures 8.5, and 8.6 show biaxial test results for dense and loose Hostun sand respectively. As expected, the HS-Small model with the Matsuoka-Nakai yield criterion
107
s3 = 100 kPa CD
s1/s3
GSecant [kN/m2]
eVol[-]
-0.20
-0.16
80000
-0.12
2
-0.08
40000
-0.04
0
0.00
0.00
0.02
0.04
0.06
0.08
s3 = 300 kPa CD
s1/s3
e1[-]
0.0001
0.01
e1-e3[-]
0.01
e1-e3[-]
0.01
e1-e3[-]
GSecant [kN/m2]
eVol[-]
-0.20
0.001
160000
-0.16
120000
-0.12
80000
-0.08
-0.04
0
0.00
40000
0.00
0.02
0.04
0.06
0.08
s3 = 600 kPa CD
s1/s3
e1[-]
0.0001
0.001
GSecant [kN/m2]
eVol[-]
-0.20
-0.16
200000
-0.12
2
-0.08
100000
-0.04
0
0.00
0.00
0.02
0.04
0.06
0.08
e1[-]
0.0001
Experiment
HS (original)
HS-Small
0.001
HS + Small-Strain Overlay
Figure 8.2: Drained triaxial tests on dense Hostun sand at confining pressures of 100,
300, and 600 kPa. Left: Stress-strain data. Right: Shear modulus reduction.
108
s3 = 100 kPa CD
s1/s3
GSecant [kN/m2]
eVol[-]
-0.20
60000
-0.16
-0.12
40000
-0.08
20000
-0.04
0
0.00
0.00
0.04
0.08
s3 = 300 kPa CD
s1/s3
e1[-]
0.0001
0.001
0.01
e1-e3[-]
0.01
e1-e3[-]
0.01
e1-e3[-]
GSecant [kN/m2]
eVol[-]
-0.20
120000
-0.16
-0.12
80000
-0.08
40000
-0.04
0.00
0.00
0.04
0.08
s3 = 600 kPa CD
s1/s3
0
0.0001
e1[-]
0.001
GSecant [kN/m2]
eVol[-]
-0.20
-0.16
200000
-0.12
2
-0.08
100000
-0.04
0
0.00
0.00
0.04
0.08
0
0.0001
e1[-]
Experiment
HS (original)
HS-Small
0.001
HS + Small-Strain Overlay
Figure 8.3: Drained triaxial tests on loose Hostun sand at confining pressures of 100, 300,
and 600 kPa. Left: Stress-strain data. Right: Shear modulus reduction.
109
s3 = 100 kPa CD
s1/s3
GSecant [kN/m2]
eVol[-]
-0.20
30000
-0.16
-0.12
20000
-0.08
10000
-0.04
0
0.00
0.00
0.05
0.10
0.15
0.20
s3 = 300 kPa CD
s1/s3
e1[-]
0.0001
0.001
0.01
e1-e3[-]
0.01
e1-e3[-]
0.01
e1-e3[-]
GSecant [kN/m2]
eVol[-]
-0.20
60000
-0.16
-0.12
40000
-0.08
20000
-0.04
0
0.00
0.00
0.05
0.10
0.15
0.20
s3 = 600 kPa CD
s1/s3
e1[-]
0.0001
0.001
GSecant [kN/m2]
eVol[-]
-0.20
120000
-0.16
-0.12
80000
-0.08
40000
-0.04
0
0.00
0.00
0.05
0.10
0.15
0.20
e1[-]
0.0001
Experiment
HS (original)
HS-Small
0.001
HS + Small-Strain Overlay
Figure 8.4: Drained triaxial tests on reconstituted Kaolinite clay at confining pressures of
100, 300, and 600 kPa. Left: Stress-strain data. Right: Shear modulus reduction.
110
s3 = 100 kPa CD
s1/s3
eVol[-]
-0.20
s3 = 200 kPa CD
s1/s3
eVol[-]
-0.20
-0.16
4
-0.16
4
-0.12
-0.08
-0.12
-0.08
-0.04
0
0.00
0.00
0.02
e1[-]
0.04
s3 = 400 kPa CD
s1/s3
-0.04
0
0.00
eVol[-]
-0.20
0.00
0.02
0.04
0.06
e1[-]
s3 = 800 kPa CD
s1/s3
eVol[-]
-0.20
-0.16
4
-0.16
4
-0.12
-0.08
-0.12
-0.08
-0.04
0
0.00
0.00
0.02
0.04
0.06
Experiment
0.08
e1[-]
-0.04
0
0.00
HS-Small(MN)
0.04
0.08
0.00
0.12 e1[-]
HS-Small(MC)
Figure 8.5: Drained biaxial tests on dense Hostun sand at confining pressures of 100, 200,
400, and 800 kPa.
111
s3 = 100 kPa CD
s1/s3
eVol[-]
-0.20
s3 = 200 kPa CD
s1/s3
eVol[-]
-0.20
-0.16
3
-0.16
3
-0.12
2
-0.12
2
-0.08
1
0
0.00
0.02
0.04
0.06
s3 = 400 kPa CD
s1/s3
0.08
-0.04
0.00
0
0.00
e1[-]
eVol[-]
-0.20
-0.08
-0.04
0.00
0.04
0.08
e1[-]
s3 = 800 kPa CD
s1/s3
eVol[-]
-0.20
-0.16
3
-0.16
3
-0.12
2
-0.12
2
-0.08
-0.08
-0.04
-0.04
0.00
0.00
0.00
0.04
0.08
Experiment
0.12 e1[-]
0.00
HS-Small(MN)
0.04
0.08
0.12 e1[-]
HS-Small(MC)
Figure 8.6: Drained biaxial tests on loose Hostun sand at confining pressures of 100, 200,
400, and 800 kPa.
112
q [kN/m2]
100
4E+5
50
0
2E+5
-50
Calculated
Rivera/Bard (initial)
-100
Rivera/Bard (reloading)
0
-5
10
-4
10
-3
10
e1[-]
-4
-5x10
-4
-3x10
0x10
-4
3x10
-4
5x10
e1[-]
Figure 8.7: Hysteresis loop in a drained triaxial test on dense Hostun sand.
(HS-Small(MN)) always predicts higher peak strengths than the model with the MohrCoulomb yield criterion (HS-Small(MC)). The test data fully validates the MatsuokaNakai approach for dense sand whereas for loose sand the strength is slightly overpredicted. Here it seems that bifurcations set in even before the limit state according to
the Matsuoka-Nakai criterion is reached.
The element tests presented in Figures 8.7, 8.8, and 8.9 recapitulate some of the tests
already presented and introduced in Chapter 6. Now, the HS-Small and the original HS
model are evaluated in these tests instead of the pure Small-Strain Overlay model.
The initial stiffness in Figure 8.7 is slightly over-predicted since the standard material
set for dense Hostun sand given in Table D.1 is used. The calculated hysteresis loop is not
closed anymore because of too much unloading plasticity. Here, as for other cyclic phenomena as well, kinematic hardening formulations should be used rather than purely
isotropic hardening models. Nevertheless, the hysteresis loop is simulated not too badly
by the (isotropic hardening) HS-Small model.
Figures 8.8, and 8.9 show numerical simulations of the tests by Richardson [148] on
reconstituted London clay. The London clay could be reasonably modeled with the data
set defined for Kaolin clay by slightly reducing the initial stiffness to Gref
0 = 26 MPa and
ref
ref
ref
= 7 MPa).
increasing the unloading-reloading stiffness to Eur = 21 MPa (E50 = Eoed
Again, there is no difference between the Matsuoka-Nakai and the Mohr-Coulomb flavor of the HS-Small model as b = 0. The improvement achieved by using the HS-Small
model is obvious. However, a pronounced distinction between stress paths AOX and
COX is missing. The HS-Small model tracks the difference in these stress paths by means
of different overconsolidation pressures pp only. It does not track the second main difference between these stress paths, which is the materials void ratio at the start of devi-
113
GTangent [MPa]
40
q [kPa]
DOX
AOX
150
100
COX
30
BOX
50
Calculated (HS)
20
0
0
100
O
200
A
300
400
p' [kPa]
-50
10
B
-100
0
10-4
10-3
10-2
e1[-]
-150
GTangent [MPa]
40
q [kPa]
DOX
150
100
AOX
COX
30
BOX
50
Calculated
(HS-Small)
20
0
0
100
O
200
A
300
400
-50
10
-100
0
10-4
10-3
10-2
e1[-]
-150
114
p' [kPa]
115
q [kN/m2]
(s1-s3)/2 [kN/m2]
s3 = 200 kPa CU
2000
300
1500
200
1000
100
500
0
0.00
0.02
q [kN/m2]
0.04
0.06
e1[-]
200
300
400
400
1000
200
(s1+s3)/2 [kN/m2]
0
0.02
q [kN/m2]
0.04
0.06
500
(s1-s3)/2 [kN/m2]
s3 = 300 kPa CU
2000
0
0.00
(s1+s3)/2 [kN/m2]
e1[-]
400
600
800
(s1-s3)/2 [kN/m2]
s3 = 600 kPa CU
3000
800
2000
400
1000
0
0.00
(s1+s3)/2 [kN/m2]
0
0.02
0.04
Experiment
HS (original)
0.06
e1[-]
HS-Small
800
1200
1600
Figure 8.10: Undrained triaxial tests on dense Hostun sand at confining pressures of 200,
300, and 600 kPa. Left: Stress-strain data. Right: s-t diagram.
116
(s1-s3)/2 [kN/m2]
s3 = 200 kPa CU
160
200
120
80
100
40
0
0.00
0
0.02
0.04
0.06
e1[-]
100
200
(s1+s3)/2 [kN/m2]
q [kN/m2]
(s1-s3)/2 [kN/m2]
s3 = 300 kPa CU
160
300
120
200
80
100
40
0
0.00
0
0.02
0.04
0.06
e1[-]
100
200
300
(s1+s3)/2 [kN/m2]
q [kN/m2]
(s1-s3)/2 [kN/m2]
s3 = 600 kPa CU
400
600
300
400
200
200
0
0.00
100
0
0.02
0.04
0.06
e1[-]
200
400
600
(s1+s3)/2 [kN/m2]
Experiment
HS (original)
HS-Small
Figure 8.11: Undrained triaxial tests on loose Hostun sand at confining pressures of 200,
300, and 600 kPa. Left: Stress-strain data. Right: s-t diagram.
117
Figure 8.12: Steinhaldenfeld NATM tunnel before concreting the inner liner.
8.2.1 Tunnels
Schmidt [160] and Peck [135] were the first to propose that green field settlement troughs
above tunnels can be approximated by a (Gaussian) distribution function of the form:
y2
s(y) = smax exp( 2 )
2i
(8.1)
where smax is the maximum vertical settlement above the tunnel axis, y is the horizontal
distance to the tunnel axis, and i is the horizontal distance of the settlement troughs
inflection point from the tunnel axis.
The empirical relationship expressed in Equation 8.1 holds very well, independent of
the applied tunnelling method. Numerical methods often fail in predicting the correct
shape of settlement troughs. Especially, if the applied soil model does not account for
small-strain stiffness, numerical methods commonly predict too flat and wide settlement
troughs. This can be observed for instance in the following two case studies, where the
HS-Small model yields significantly improved results compared to the HS model.
8.2.1.1 Steinhaldenfeld NATM tunnel
Figure 8.12 shows Steinhaldenfeld tunnel, a tunnel near Stuttgart that is constructed using the new Austrian tunnelling method (NATM). The soil layers together with the FE
mesh used in the plane strain analysis of Steinhaldenfeld tunnel are shown in Figure
8.13. The lower soil layers (Keuper marl and Limestone) are overconslidated. The overconsolidation ratio is assumed to be 2.0. The material behavior is drained in all the soil
layers. A complete list of the soil layers material parameters used in the calculation is
given in Table D.2.
The weak marl around the tunnel excavation was improved by rock bolts during the
NATM construction process. The strength increase due to rock bolting is modeled by
118
10.00
6.20 m
20.00
Lacustrine Limestone
0.00
10.00
20.00
30.00
40.00
introducing additional cohesion of 25 kPa in the bolted area. The tunnels linear elastic
lining plate elements have a stiffness of EA = 3.75 106 , a flexural rigidity of EI =
1.95 104 , a Poissons ratio of = 0.2, and a unit weight of C = 25 kN
.
m3
The meshed problem area below the tunnel is relatively small since the Limestone is
much stiffer than the overlaying marl. For the mesh boundaries it is generally proposed
here to require the stiffness multiplier:
Gm 0.8
G0
Gur
(8.2)
in the meshs outermost stress points (at boundaries). For a definition of the stiffness
multiplier Gm refer to Section 7.3.1.
For the numerical analysis of Steinhaldenfeld tunnel, the stress reduction method is
adopted. Using this method, the tunnel excavation remains unsupported by the lining
till, an artificial support pressure p0 is reached, where p0 is the initial rock mass pressure before excavation. Then, the tunnel lining is activated and the support pressure is
removed completely. The factor 0 < < 1 is called the load reduction factor.
Figure 8.14 shows results of the numerical analysis. The settlement trough calculated
by the HS-Small model is clearly steeper and more accurate compared to that calculated
by the original HS model. The remaining discrepancy to the measured settlement trough
can possibly be explained by the inadequateness of 2D analysis for 3D problems. Also
the small-strain stiffness parameters could be underestimated because they are based
solely on literature data and correlations. Set 2 in Figure 8.14 therefore indicates the
changes in settlement from a 40 percent increase of the initial stiffness, displayed in Table
D.2 (Set 1). The failure criterion (Mohr-Coulomb or Matsuoka-Nakai) has almost no
influence on the shape of the settlement trough. However, it has an influence on the load
reduction factor. This is chosen so that in all calculations the settlement directly above
119
10
20
30
-0.01
-0.01
Measurement
-0.02
HS (original)
10
20
30
40
Measurement
-0.02
HS (original)
Settlement [m]
Settlement [m]
Figure 8.14: Green field settlement trough due to construction of Steinhaldenfeld NATM
tunnel. Left: HS-Small(MC) prediction. Right: HS-Small(MN) prediction.
the tunnel axis agrees with the measurement. It is found that the load reduction factor is
highest for the HS model = 0.36, followed by the HS-Small(MC) model = 0.32, and
the HS-Small(MN) model = 0.25. Unfortunately the load reduction factor cannot be
experimentally derived. The difference in the calculated factors is however reasonable:
The Matsuoka-Nakai criterion yields the highest strength and stiffness in the given plane
strain conditions; its load reduction factor should therefore be least.
8.2.1.2 Heinenoord Tunnel
The second Heinenoord tunnel, south of Rotterdam is a shield tunnel project in largely
sandy soils. Figure 8.15 shows the finite element mesh, used in the following plane strain
analysis. Figure 8.15 depicts also the different soil layers, for which Table D.3 provides
HS material parameters based on the work by Bakker [10]. Drained material behavior
is assumed in all soil layers. The lining of the second Heinennord tunnel is modeled by
linear elastic plate elements with EA = 3.75 106 , EI = 1.95 104 , = 0.2, and = 25 kN
.
m3
Several numerical simulation methods for the analysis of shield tunnels can be found
120
GW
Fill
10.00
Sand
8.30 m
20.00
Clay
0.00
10.00
20.00
30.00
40.00
50.00
Small(MN) analysis for its highest strength and stiffness. The highest roof pressure of
130 kPa is obtained in the analysis using the HS model, in between is the HS-Small(MC)
analysis with a pressure of 129 kPa. All pressures seem very reasonable as they are of
the same size as the overburden pressure.
The Inclinometer measurements (Figure 8.16) at a distance of 6, 10, and 16 meters from
the tunnel axis, as well as the surface settlement troughs (Figure 8.17) underline again the
HS-Small models improved precision in displacement analysis. The HS-Small model is
not perfect though. The Inclinometer measurement at a 6 meter distance from the tunnel
axis show a somewhat different characteristic to the calculated one. Nevertheless, the
decay of displacement with distance from the tunnel axis is modeled nicely. This is
entirely due to the new small-strain formulation. The choice of failure criterion, besides
its influence on the total material stiffness (see roof pressure), has a negligible effect on
the calculated strain distribution.
Due to missing experimental data, the two tunnel examples presented above are not
suited to validate the HS-Small model completely. However, it can be concluded that the
HS-Small model can improve displacement analysis considerably by causing a diffuse
8.2.2 Excavations
Simpson introduced his brick model in a case study of the British Library deep excavation. The House of Parliament underground car park case study is another deep excavation study, that is often cited along with small-strain stiffness. The tilt of Big Ben clock
tower right next to this excavation illustrates nicely the effects of advanced (small-strain
121
0.01
-0.01
0.01
-0.01
-10
0.01
0
-10
-10
Experiment
HS
HS-Small (MC)
-20
HS-Small (MN)
-20
-20
b
Depth [m]
Depth [m]
Depth [m]
Figure 8.16: Results from Inclinometer measurements next to the second Heinenoord
tunnel. Horizontal displacements are displayed at a distance of (a) 6 m,
(b) 10 m, and (c) 16 m from the tunnel axis.
10
20
30
-0.01
Settlement [m]
20
30
40
-0.01
Measurement
-0.02
10
HS (original)
Measurement
-0.02
HS (original)
Figure 8.17: Green field settlement trough due to construction of the second Heinenoord
tunnel. Left: HS-Small(MC) prediction. Right: HS-Small(MN) prediction.
122
30.00
excavation step 1
Sand (Layer 1)
-3.00 (GW)
-4.80
excavation step 2
-9.30
excavation step 3
-14.35
excavation step 4
-16.80
27
27
20
Sand (Layer 2)
40
27
8.0
0
A1
8.0
8.0
-17.90 (GW)
A2
A3
hydraulic barrier
60
-30.00
-32.00
80 diaphragm wall
Anchor Prestress
80
A1
A2
A3
Sand (Layer 3)
100
-20
20
40
60
80
100
120
768 kN
945 kN
980 kN
Distance
2.30 m
1.35 m
1.35 m
Length l
15 cm
15 cm
15 cm
19.80 m
23.30 m
23.80 m
Figure 8.18: Excavation in Berlin sand: 2D mesh and geometry detail (right).
stiffness) models. These can predict the total excavation heave more correctly and thus
do not suggest that the clock tower tilts away from the pit, which less refined models
may do. Beyond the constitutive model however, there are many more things to consider in the analysis of deep excavations, for example; the interface between soil and
retaining structure, initial stress conditions, structural elements, time effects, and pore
water pressures. Some of these issues should always be thoroughly considered in excavation analysis (interface, initial stresses, structural elements), others are less critical in
less difficult soil conditions. Two examples with less critical soil conditions are chosen
for the validation process of the HS-Small model: An excavation in Berlin sand and an
excavation in Rupel boom clay.
8.2.2.1 Excavation in Berlin sand
The working group 1.6 Numerical methods in Geotechnics of the German Geotechnical
Society (DGGT) has organized several comparative finite element studies (benchmarks).
One of these benchmark examples is the installation of a triple anchored deep excavation
wall in Berlin sand. The reference solution by Schweiger [161] is used here as the starting
point for the next validation example: Both, the mesh shown in Figure 8.18, and the
soil parameters given in Table D.3 are taken from this reference solution. However, the
bottom soil layer defined by Schweiger could be omitted in the analysis when using the
HS-Small model. In the reference solution this layers only purpose is the simulation of
small-strain stiffness due to a lack of small-strain stiffness constitutive models back then.
Figure 8.19, and 8.20 show results from the finite element calculation using the HSSmall(MC), and HS-Small(MN) model respectively. The original reference solution is
calculated with the HS model (Figure 8.19, Table D.4). A second reference solution utilizing the HS model with Matsuoka-Nakai failure criterion, but no small-strain stiffness
123
20
40
0
0
-0.004
-0.008
-10
-0.012
-0.016
-20
Settlement [m]
Measurement
HS (Reference)
-30
HS-Small (MC)
Depth below surface [m]
Figure 8.19: Excavation in Berlin sand. HS-Small(MC) model predictions versus measured displacements after the final excavation step. Left: Surface settlement
trough. Right: Lateral wall deflection.
formulation (HS(MN)) is shown in Figure 8.20. From the comparison of HS-Small results
with the reference solutions, the impact, which a) the small-strain stiffness, and b) the
failure criteria have can now be derived.
The small-strain stiffness formulation accumulates more settlements right next to the
wall, whereas the settlement trough is smaller. The triple anchored retaining wall is generally deflected less when using the two HS-Small models. The deflection calculated
by the HS-Small(MC) almost fits that measured. The Matsuoka-Nakai criterion in the
HS-Small(MN) model gives less displacements compared to the Mohr-Coulomb criterion. In comparison to the measured wall deflection, the HS-Small(MN) model even
gives lower displacements. However, this might be related to the problems original reference solution. The chosen parameters are most likely partly back-analyzed from HS
model calculations. In this case the Mohr-Coulomb strength parameters already include
an increase for plane strain conditions. Additionally increasing the strength by using the
HS-Small(MN) model may thus overdo things.
Calculated excavation heaves at the end of excavation are shown in Figure 8.21. The
heave, which is due to elastic unloading, is roughly halved when using the HS-Small
models, independent of the failure criterion chosen. From Figure 8.21 it can be further
concluded that the meshed problem area is not quiet sufficient in the example since the
heave is not asymptotically decaying to zero. Figure 8.22 verifies in this respect the
criterion proposed in Equation 8.2. The monotonic stiffness multiplier Gm at the end of
excavation is well below 3 at the lower mesh boundary.
However, defining the mesh boundary in a depth of 2-3 times the excavation width as
124
20
40
60
-0.004
-0.008
-10
-0.012
-0.016
-20
Settlement [m]
Measurement
-30
HS (MN)
HS-Small (MN)
Figure 8.20: Excavation in Berlin sand. HS-Small(MN) model predictions versus measured displacements after the final excavation step. Left: Surface settlement
trough. Right: Lateral wall deflection.
Section A-A
20.00 m
0.01
0.02
0.03
0.04
-20
-40
-60
-80
HS (Reference)
HS (MN)
HS-Small (MC)
-100
HS-Small (MN)
Figure 8.21: Excavation in Berlin sand: Vertical displacements in the excavation pit at a
distance of 10 m from the retaining wall.
125
3.0
20
G m [-]
40
2.0
60
80
1.0
100
-20
20
40
60
80
100
120
Figure 8.22: Excavation in Berlin sand: Factorized increase of unloading/ reloading stiffness at the end of excavation.
proposed by the AK 1.6, and seen in the reference solution by Schweiger, seems reasonable, too. The error in the excavation heave is relatively small. A mesh according to the
requirements of the criterion defined in Equation 8.2 is used in the next example of an
excavation in Ruple clay.
8.2.2.2 Excavation in Ruple clay (Offenbach)
Figure 8.23 shows the geometry of an excavation in Offenbach. Similar to the well known
subsoil conditions of the neighboring town Frankfurt, the predominant soil here is here
a highly overconsolidated clay. Different to Frankfurt clay, the Rupel clay is a oligocene
clay, with no solid rock intrusions. Besides the excavations geometry and measurement
data, Krajewski [95] also published a material parameter set for Rupel clay that can be
used together with the HS model. The published material parameters have been obtained by back-analysis of lab experiments and a large hopper, built at a neighboring
chemical plant. Due to the axisymmetric conditions of all back-analyzed experiments, it
can be assumed that the parameter set given is not adjusted for plane strain conditions
(as it might have been the case in the Berlin sand excavation). Again, there is no smallstrain stiffness testing data available. The parameter set given in Table D.5 is therefore
supplemented by correlated small-strain stiffness data. Figure 8.24 shows the 2D plane
strain mesh used in the numerical analysis of the excavation.
The retaining wall in the example is a bored pile wall with pile diameter r = 0.90m
and axial distance d = 1.30m. The struts are IPBs 340 and IBPs 360 with a clearance of
2.45 m, 1.95 m, and 1.80 m in the top, second, and third strut layer respectively.
The Rupel boom clay has an overconsolidation ratio (OCR) well above 2. Using the
relationship
(8.3)
K0oc = (1 sin )(OCR)sin
126
p =30 kN/m
12.00
0.00
p = 39 kN/m
-2.50
Fill
Rupel boom clay
-2.80 (GW)
IPB 340
-7.60
IPB 360
-13.00
IPB 360
Bored pile wall
Pile diameter: 90 cm
Axial distance: 1.20 m
-19.30
-22.70
Inclinometer
-30.00
Figure 8.23: Excavation in Rupel clay: Geometry of the Offenbach excavation for a railway tunnel.
0
20
40
60
Rupel boom clay
80
-100
-80
-60
-40
-20
20
40
60
80
100
Figure 8.24: Excavation in Rupel clay: 2D mesh used in the plane strain analysis. Neighboring buildings to the excavation are modeled by a line load at the top
ground surface.
127
-0.02
-5
-0.04
0.04
0.08
-10
-0.06
-15
Settlement [m]
Measurement
HS (original)
-20
HS-Small (MC)
HS-Small (MN)
Figure 8.25: Excavation in Rupel clay. HS-Small model predictions versus measured displacements after the final excavation step. Left: Surface settlement trough.
Right: Lateral deflection of the left retaining wall.
presented by Mayne & Kulhawy [115], the K0 value for generating the initial stress field
is calculated as 0.95, for an assumed OCR of 2.8. Because of stress relaxation, no initial
elastic clay behavior is assumed. Consolidation is considered only to the extend that no
excess pore pressures built up behind the wall.
Figure 8.25 shows the calculated results in comparison to the available measurements
at the end of excavation. The original HS model overall gives too much surface settlements accompanied by too much lateral deformations in the retaining structure. The
HS-Small(MC) model improves results significantly, the HS-Small(MN) model gives almost perfect results. The settlements to the right hand side of the excavation are still
too big, probably due to a slightly overestimated surcharge at the ground surface. Analogous to the previous example, the settlements right next to the wall increase slightly
when using the HS-Small model. This is again related to the reduced excavation heave,
which is shown in Figure 8.26. Excavation heave is calculated to be lowest for the HSSmall(MN) model and highest for the HS model.
Figure 8.27 shows the elastic material stiffness at the end of excavation. Here, stiffness
is plotted as a factorized increase of the unloading-reloading stiffness (see Equation 6.13).
Stiffness increases due to strain reversals are not considered in this representation. In
agreement with the criterion defined in Equation 8.2, all calculated displacements (see
Figure 8.25, and 8.26) are negligible in areas where Gm > 0.8 GGur0 .
In conclusion, the HS-Small model(s) performed very well in both excavation examples, although there has been no small-strain stiffness test data available at all. All small-
128
0.04
0.08
0.12
0.16
-20
-40
-60
HS (Reference)
HS (MN)
-80
HS-Small (MC)
HS-Small (MN)
Figure 8.26: Excavation in Rupel clay: Vertical displacements in the excavation pits symmetry.
4.0
0
20
G m [-]
3.0
40
2.0
60
80
1.0
-100
-80
-60
-40
-20
20
40
60
80
100
Figure 8.27: Excavation in Rupel clay: Factorized increase of unloading-reloading stiffness at the end of excavation.
129
130
5.00
10.00
15.00
0.00
5.00
10.00
15.00
20.00
Figure 8.28: Load test on a spread footing on sand. Left: Texas A&M Universitys National Geotechnical Experimental Site. Right: 2D axisymmetric mesh used
in the finite element analysis.
available CPT and Pressuremeter field tests it can be derived that the sandy site is overconsolidated. The amount of overconsolidation is here back-analyzed from the load settlement curves. Yet, the stiffness in initial loading and in unloading-reloading cycles is
still underestimated. Increasing the initial stiffness derived from laboratory test data by
40% yields the results presented in Figure 8.30.c. Now, the unloading-reloading stiffness
is simulated almost perfectly. However, the overall load settlement curve overestimates
soil stiffness slightly, which may be due to the axisymmetric simplification of the given
3D problem.
As shown in Chapter 3, stiffness discrepancies in the magnitude found in the A&M
footing example is not very uncommon. Nevertheless, it should be noted that underestimating small-strain stiffness by using laboratory data or correlations is always better
than not accounting for small-strain stiffness at all. This, for example, can also be concluded from Figure 8.30: The HS-Small model gives more realistic simulation results
than the original HS model regardless of initial stress and stiffness is used.
From a practical point of view, the most important difference in the load settlement
curves is probably their initial slope (assuming that ultimate load design can also be
handled by other methods). In the given example, service load settlements are overestimated by more than 80% when not using the HS-Small model.
The Texas A&M spread footing example is a nice demonstration of the benefits of
small-strain stiffness on the one hand. On the other hand, it demonstrates the sensitivity
of numerical calculations to initial stress conditions. As there is also a link between
initial stress conditions and a materials strain history, the initialization of strain history
is briefly discussed in the next section.
131
s3 = 34.5 kPa CD
s1-s3
GSecant [kN/m2]
eVol[-]
-0.03
80
80000
60000
-0.02
40000
40
-0.01
0.04
0.08
20000
0.12
0.00
e1[-]
s3 = 138 kPa CD
s1-s3
0.0001
0.01
e1-e3[-]
GSecant [kN/m2]
eVol[-]
400
0.001
-0.03
120000
-0.02
80000
-0.01
40000
300
200
100
0.04
0.08
0.12
0.00
e1[-]
s3 = 345 kPa CD
s1-s3
0.0001
800
0.01
e1-e3[-]
GSecant [kN/m2]
eVol[-]
-0.03
0.001
200000
150000
-0.02
100000
400
-0.01
0
0.00
50000
0.00
0.04
0.08
0.12
e1[-]
0
0.0001
0.001
0.01
e1-e3[-]
HS (original)
Experiment
HS-Small
Figure 8.29: Triaxial tests on sand samples from the A&M test site (0.6 m sampling
depth). Small-strain stiffness parameters are obtained by averaging the results from several resonant column tests as shown on the right hand side.
132
4000
8000
4000
8000
0
Load [kN]
-0.04
-0.04
-0.08
-0.08
-0.12
-0.12
1a
Load [kN]
4000
1b
Settlement [m]
8000
4000
Settlement [m]
8000
Load [kN]
-0.04
-0.04
-0.08
-0.08
-0.08
-0.12
-0.12
-0.12
2b
Settlement [m]
Experiment
8000
Load [kN]
Load [kN]
-0.04
2a
4000
Settlement [m]
HS (original)
2c
Settlement [m]
HS-Small
Figure 8.30: Calculated and experimentally derived load-settlement curves for the A&M
3x3 m footing on sand (HS - top row; HS-Small - bottom row): a) NC sand,
b) OC sand with laboratory small-strain stiffness, c) OC sand with increased
(in-situ) small-strain stiffness.
133
0.00
0.00
110.00
110.00
98.00
5.00
98.00
5.00
86.00
86.00
74.00
74.00
10.00
10.00
62.00
62.00
50.00
15.00
50.00
15.00
G [kN/m]
0.00
5.00
10.00
15.00
20.00
G [kN/m]
0.00
5.00
10.00
15.00
20.00
Figure 8.31: Initial elastic stiffness of the HS-Small model obtained by different initialization procedures in the A&M footing example. Left: No strain history
reset between removing the surcharge used to model the materials geologic
history and construction of the footing. Right: Full strain history reset by
applying a small reverse load step before introducing the footing.
134
A = 4500 kN (reloading)
0.00
0.00
110.00
110.00
98.00
5.00
98.00
5.00
86.00
86.00
74.00
74.00
10.00
10.00
62.00
62.00
50.00
15.00
50.00
15.00
G [kN/m]
0.00
5.00
10.00
15.00
20.00
G [kN/m]
0.00
A = 8000 kN
5.00
10.00
15.00
20.00
A = 12000 kN
0.00
0.00
110.00
110.00
98.00
5.00
98.00
5.00
86.00
86.00
74.00
74.00
10.00
10.00
62.00
62.00
50.00
15.00
50.00
15.00
G [kN/m]
0.00
5.00
10.00
15.00
20.00
G [kN/m]
0.00
5.00
10.00
15.00
20.00
An overburden surcharge will therefore trigger secondary loading although the strain
history is reset. If that is not desired, the monotonic strain history can be manually
reset or an initial stress procedure as offered for example in the finite element program
P LAXIS, can be used.
In discussing the strain historys reset, loading-unloading-reloading cycles as seen in
the A&M footing example, are also worth a closer look. Strain history is fully reset in 180
degree load reversals. Due to isotropic hardening, the accumulated strains are usually
less in reloading than in primary loading. Strain history is hence more salient in virgin
loading than in reloading. This can be observed for example in the decay of elastic
material stiffness shown in Figure 8.32.
135
136
s1[kN/m2]
s1[kN/m2]
800
800
600
600
400
400
200
200
0
0.000
0.004
0.008
Experiment
0.012
ev[-]
0
0.000
HS (original)
0.004
0.008
0.012
ev[-]
HS-Small
Figure 8.33: Isotropic compression test on dense (left) and loose (right) Hostun sand.
137
Chapter 9
3D case study Sulfeld
Lock
The Mittellandkanal (MLK) links the rivers Rhine, Ems, Weser, and Elbe to Berlin and to
the channels peak altitude of 65 m to 56 m height above sea level. The existing Sulfeld
Lock was constructed between 1934-37 and was opened in 1938. The locks southern
chamber is currently being rebuilt in order to adapt it to raise ship traffic and generations
of new ships. The northern lock chamber including its six economizing basins is to be
kept fully operational during the 4 year construction period. The high-speed railway
link between Hanover and Berlin lies close to the construction site and also needs to be
kept operational at all times during construction. Figures 9.1 and 9.2 show the Sulfeld
lock and the high-speed railway link passing nearby.
Investigation of the site shows that overconsolidated glacial deposits overlie Lias bedrock. The glacial deposits are composed of clayey silt, sand, and glacial till (Figure 9.3).
The cohesive layers (silt and glacial till) were characterized in laboratory tests. A two-
139
Chapter 9 3D case study Sulfeld
Lock
step model calibration procedure was employed: First, all test results were statistically
averaged before secondly, selected tests were back analyzed in numerical element tests.
Properties of the sand layer were estimated from in-situ heavy dynamic penetration test
results. The obtained material parameters were finally calibrated in a 2D finite element
model by comparing calculated displacements of the existing lock during operation to
measured ones. In the calibration procedure, the confidence in the correlated smallstrain stiffness parameters could be specifically increased. An overview of the final material sets used in the 3D analysis is given in Table D.7.
The excavation for the new lock chamber is fully housed by vertical retaining walls.
The selected type of retaining wall is closely related to the type and presence of adjacent
structures. A tied back diaphragm wall is built in sections, close to structures that are
sensitive to settlements. The most sensitive part of the excavation pit next to the railway
link is strutted. A cut-off wall is built in all sections where settlements are tolerable. If
required by design, the cut-off walls are reinforced by sheet piles.
The HS-Small model is deployed here in a 3D finite element analysis of the interaction
between the excavation pit and railway bridge abutments. As solely the Mohr-Coulomb
yield criterion is used in the following, the term HS-Small refers to HS-Small(MC) within
this chapter. Figure 9.3 shows the geometry, structural features, and excavation stages
modeled within the finite element code ABAQUS. The structural features include two
abutments on pile foundations, back-anchored and strutted retaining walls, sheet pile
walls, and a floor slab in the strutted section. Embedded elements are used to model the
floor slab, sheet pile walls and piles. All other structural features and the subsoil are dis-
140
Ship channel
Abutments with piled found.
80 m
Excavation
+57 mNN
Fill
Glacial till
27 m
2.0 m (step 1)
Sand
Clayey silt
4.0 m (step 2)
+46 mNN
5.0 m (step 3)
17.7 m
70 m
Figure 9.3: Geometry and Geology of the Sueldfeld excavation in the 3D FE model. Lias
bedrock marks the lower boundary of the model.
cretized by full integrated volume elements. Four excavation stages are introduced: One
before wall construction (pre-excavation), and three after wall construction (excavation
step 1 to 3). Construction of the retaining walls is not modeled accurately. Instead, they
are placed as a whole. The phreatic level is adjusted several times in the calculation due
to groundwater discharge in the bedrock or groundwater lowering in the fill. After the
final excavation step, the groundwater table is lowered due to works near the pumping
station. This change in phreatic level is the final load step considered in the analysis. A
list of all calculation stages is given in Table 9.1. As pore water measurements indicate
a relatively fast dissipation of excess pore pressures in the vicinity of the excavation pit,
drained conditions are assumed.
The Sulfeld
analysis main focus is on displacements of the railway bridge abutments
and their backfill. Due the rather unprecise numerical retaining wall placement method,
measured and calculated abutment displacements are compared for excavation stages 1
141
Chapter 9 3D case study Sulfeld
Lock
Action
1
2
3
4
5
6
7
8
9
10
Initialization
Pre-excavation (Excavation to +57 m NN)
Installation cut-off wall
Phreatic level at +53 m NN (local vacuum)
Installation diaphragm wall
Phreatic level at +57 m NN
Excavation to +55 m NN (excavation step 1)
Excavation to +51 m NN (excavation step 2)
Excavation to +46 m NN (excavation step 3)
Phreatic level at +53 m NN (channel section only)
Bottom slab
Excavation step 3
8
4
0
-4
-8
06.2005
07.2005
08.2005
09.2005
Figure 9.4: Heave measured in the excavation pit during excavation step 3.
to 3 and the final change in phreatic level only. The well equipped site allows for various
other comparisons, as well. The excavation heave for example, is presented in Figure 9.4.
From extensometer readings in the excavation pit, the final heave during excavation step
3 can be quantified as 23 mm. In the analysis, heave 4 mm is obtained when using
the HS-Small model. The HS analysis over-predicts the excavation heave by more than
8 mm. The overall vertical displacements from the two analysis are compared in Figure
9.5. The measured settlement (extensometer) of the abutment closest to the excavation
pit is shown in Figure 9.6. From geodetic measurements it can be additionally derived
that the abutments are not tilting. Exactly the same behavior is found in the HS-Small
analysis: Almost no tilt of the abutments occurs with a total settlement of 2 mm. The
HS calculation on the other hand, shows settlements of up to 10 mm in combination with
tilt.
Figure 9.8 shows a comparison of measured and calculated horizontal displacements
of both bridge abutments. Measured displacements were separated in displacements
due to pre-excavation and wall construction (not shown) and displacements due to ex-
142
uz [cm]
6.5
4.5
2.5
Figure 9.5: Vertical displacements due to excavation (step 1 to 3) and groundwater lowering.
settlement [mm]
Pre-excavation
(wall constr. etc.)
Excavation Lowering
step 1-3
GW
0
1
2
3
4
01.2004
01.2005
01.2006
143
Chapter 9 3D case study Sulfeld
Lock
a.) Vertical displacement uz - HS
uz [mm]
0.0
-5.0
-10.0
-15.0
Figure 9.7: Vertical displacements of the abutments due to excavation and groundwater
lowering.
cavation step 1 to 3 and final groundwater lowering. The latter can be directly compared
to the deformed abutments drawn to scale for both, HS and HS-Small analysis. The overall trend in the HS-Small calculation is in better agreement with the measured displacements compared to those of the HS model. Right next to the abutment front walls, the
calculated displacements resemble reasonably well the measured displacements. However, the wing walls do not deform extensively enough, which could be partly due to the
linear elastic concrete material model chosen in the calculation.
Figure 9.9 finally depicts the factorized increase of the unloading-reloading stiffness
at several calculation stages. For illustration purposes, the factorized stiffness increase is
shown in the clayey silt only. It is only towards the end of excavation, that small-strain
stiffness disappears with greater distance to the excavation pit. Remarkably is also the
effectiveness of the pile foundation. Within the reinforced part of the soil, almost no
stiffness reduction can be observed.
144
5 mm
4 mm
11 mm
5 mm 4 mm
18 mm
8 mm
5 mm
4 mm
11 mm
5 mm 4 mm
18 mm
8 mm
145
Chapter 9 3D case study Sulfeld
Lock
Gm [-]
6.0
3.5
1.0
Legend:
Construction stages
1: Pre-excavation
2: Excavation step 1
3: Excavation step 2
4: Excavation step 3
5: Lowering GW
146
Chapter 10
Conclusions
The strain range in which soils can be considered truly elastic, i.e where they recover
from applied straining almost completely, is very small. With increasing strain amplitude, soil stiffness decays non-linearly: Plotting soil stiffness against log(strain) yields
characteristic S-shaped stiffness reduction curves. Hardin & Drnevich [52] first described
these small-strain stiffness reduction curves with a simple hyperbolic formula. Since
then, many additional experimental findings have been reported in the literature. These
findings concern mainly the qualitative and quantitative influence of various parameters on small-strain stiffness and stiffness reduction curves. Void ratio and mean stress,
for example, are two parameters that strongly affect small-strain stiffness. The influence
of many other parameters, such as cementation, grain size distribution, plasticity index
etc. is likewise appreciated, but is less straightforward to quantify. Recently, inherent
and stress induced anisotropic features of small-strain stiffness have became important
themes in research. At the same time, engineers working in practice have only very
limited possibilities to include small-strain stiffness in design. There is a lack of simple
and capable small-strain stiffness models. Those small-strain stiffness models available
(Chapter 5) are either research orientated, and/or applicable to specific loading paths
only. The objective of this thesis was therefore, to develop a capable small-strain stiffness model suitable for the engineering community.
The Small-Strain Overlay model proposed in this thesis generalizes the hyperbolic
Hardin-Drnevich model to multi-axial strain space. Although it takes into account strain
induced anisotropy, the overlay model renders small-strain stiffness isotropic. It can
therefore be easily combined with many existing elastoplastic constitutive models that
are based on isotropic elasticity. In fact, these are the kind of models most commonly
used in engineering practice today. One of these is the Hardening Soil (HS) model, as
implemented in the finite element code P LAXIS V8 [21]. The combination of the HS
model with the small-strain overly model introduced in this thesis, is named the HSSmall model. The HS-Small model was validated in a single stress point and various
boundary value problems. In these, all effects that are commonly attributed to smallstrain stiffness in soil-structure interaction could be recognized: The width and shape
of settlement troughs, for example, are more precisely modeled and excavation heave
is reduced to a realistic value. The overall reliability of numerical displacement analysis is considerably increased. Analysis results are also less sensitive to the choice of
proper boundary conditions. Large meshes no longer cause extensive accumulation of
displacements, because marginally strained mesh parts are very stiff. Conserved smallstrain stiffness at the meshs boundaries is an indicator for its proper extent.
147
Chapter 10 Conclusions
The HS-Small model with the Mohr-Coulomb failure criterion is recommended for
all applications where the HS model could also be used. The HS-Small model with the
Matsuoka-Nakai failure criterion seems to perform a little better than its Mohr-Coulomb
counterpart. Its increased strength and stiffness should however be further validated
before it is applied in design. If, for a certain design situation the original HS model is
already at hand, it may be beneficial to also stay with the HS-Small(MC) model: Plane
strain calibrated strength and stiffness parameters yield overly high material stiffness
and strength using the HS-Small(MN) model.
Application of the HS-Small model is straightforward. Its Small-Strain Overlay extension requires only two additional parameters. These are the very small-strain stiffness
G0 and the threshold shear strain 0.7 at which stiffness is reduced to 0.7G0 . In case no
small-strain stiffness in-situ or laboratory test data are available, these two parameters
can be correlated with various other soil or test parameters. In the validation process
of the HS-Small model, the correlations summarized in Chapter 3 of this thesis could
be successfully applied. Although these correlations bear uncertainties, it seems advisable to use them if the only alternative is to neglect small-strain stiffness. Neglecting
small-strain stiffness certainly produces higher analysis errors than using correlations
with slight uncertainties. Most correlations are at the lower end of in-situ determined
small-strain stiffness anyway.
The adaptation of G0 and 0.7 to parameters that affect them within the Small-Strain
Overlay model could be improved in future work. In the present implementation, no
void ratio dependency is assumed. This is in clear contradiction to the findings presented in Chapter 3. However, as most models that are currently used in engineering
practice assume a constant void ratio, this concept is also adopted in the present implementation. The main drawback of void ratio independent constitutive models is that
their material parameters have to be adjusted to the actual void ratio via manual user
input. Mean stress dependency of small-strain stiffness is taken into account using the
same power law as for large strains. In this way, the number of additional material parameters is kept at a minimum. For clays, the large strain power law exponent might
be slightly too high, though. The threshold shear strain on the other hand, is not explicitly corrected for mean stress in the present implementation. The latter compensates
for overly high power law exponents. Future work on the Small-Strain Overlay model
should also include further validation of the strain history mapping presented in Chapter 6. This can presently not be undertaken due to a lack of experimental data.
Future work on the HS-Small model might include the formulation of a cap-type yield
loci similar to that of the Cam-Clay model. This would allow for degenerating the HSSmall model to a full Cam-Clay model. Kinematic hardening should be introduced to
open up the HS-Small model to cyclic loading applications as well. Within a preliminary
implementation, the possibilities of introducing bounding surface plasticity to the HSSmall model are shown in [186]. The resulting model is illustrated in Figure 10.1.
In summary, the HS-Small model is a suitable tool to incorporate small-strain stiffness
into routine design. Testing procedures and correlations for determining its small-strain
parameters are described in Chapter 3. From experience so far, no situations are known
where the HS-Small model gives less reliable results than the often validated HS model.
148
Figure 10.1: Bounding surface model with cap in principal strain space. Only the highlighted cap portion is active. Right: View inside the model onto yield-,
dilatancy-, critical-, and bounding surface.
In contrary, analysis results are throughout, in better agreement with experimental findings when the HS-Small model is deployed. Combining the Small-Strain Overlay model
with other existing elastoplastic models should turn out to be equally beneficial.
149
Bibliography
[1] A. Al-Tabbaa. Permeability and stress strain response of speswhite kaolin. PhD thesis,
University of Cambridge, 1987.
[2] A. Al-Tabbaa and D. Muir Wood. An experimentally based bubble model for
clay. In S. Pietruszczak and G. N. Pande, editors, Numerical Models in Geomechanics:
NUMOG III, pages 9199. London: Elsevier Applied Science, 1989.
[3] I. Alpan. The geotechnical properties of soils. Earth-Science Reviews, 6:549, 1970.
[4] D.G. Anderson and K.H. Stokoe. Shear modulus: A time-dependent soil property,
chapter Dyn. Geotech. Test., pages 6689. Number STP 654. ASTM Spec. Tech.
Publ., Philadelphia, 1978.
[5] D.G. Anderson and R.D. Woods. Time-dependent increase in shear modulus of
clay. Proc. ASCE: Journal of the Geotechnical Engineering Division, 102(GT5):525537,
1976.
[6] J.H. Argyris, G Faust, J Szimmat, P Warnke, and K William. Recent developments
in finite element analyses of prestressed concrete reactor vessels. Nuclear Engineering and Design, 28:4275, 1974.
[7] J. H. Atkinson. Non-linear soil stiffness in routine design. Geotechnique, 50(5):
487508, 2000.
[8] J.H. Atkinson and J.A. Little. Undrained triaxial strength and stress-strain characteristics of a glacial till soil. Can. Geotech. J., 25(3):428439, 1988.
[9] J.H. Atkinson and G. Sallfors. Experimental determination of soil properties. In
Proc. 10th ECSMFE, volume 3, pages 915956, Florence, 1991. General Report to
Session 1.
[10] J.K. Bakker. Soil retaining structures - development of models for structural analysis.
PhD thesis, Technical Univesity of Delft, 2000.
[11] G. Baldi, R. Bellotti, V.N. Ghionna, M. Jamiolkowski, and D.C.F. Lo Presti. Modulus of sands from cpts and dmts. In Proc. 12th Int. Conf. SMFE, volume 1, pages
165170, Rio de Janeiro, 1989.
[12] J.P. Bardet. Lode dependences for isotropic pressure-sensitive elastoplastic materials. Transactions of the ASME, 57(9):498506, 1990.
151
Bibliography
[13] K. Been and M.G. Jefferies. A state parameter for sands. Geotechnique, 35(2):99112,
1985.
[14] R. Bellotti, M. Jamiolkowski, D.C.F. Lo Presti, and D.A. ONeill. Anisotropy of
small strain stiffness in ticino sand. Geotechnique, 46(1):115131, 1996.
[15] T. Belytschko, W.K. Liu, and B. Moran. Nonlinear Finite Elements for Continua and
Structures. John Wiley & Sons Ltd., Chichester, 2000.
[16] J. Biarez and P.-Y. Hicher. Elemantary Mechanics of Soil Behaviour. Balkema, 1994.
[17] L. Bjerrum. Problems of soil mechanics and construction on soft clays. In Proc. 8th
Int. Conf. SMFE, volume 3, pages 111159, Moscow, 1973.
[18] P.G. Bonnier. Implementational aspects of constitutive modelling. SCMEP Workshop No 1 at NTNU, August 2000.
[19] J.-L. Briaud and R. Gibbens. Large scale load tests and data base of spread footings on sand. Final Report FHWA-RD-97-068, U.S. Department of Transportation.
Federal Highway Administration, 11 1997.
[20] E.G.M. Brignoli, M. Gotti, and K.H. Stokoe. Measurement of shear waves in laboratory specimens by means of piezoelectric transducers. Geotechnical Testing Journal,
19(4):384397, 1996.
[21] R.B.J. Brinkgreve, editor. P LAXIS, 2D Version 8. AA. Balkema, 2002.
[22] B.B. Broms and A.O. Casbarian. Effects of rotation of the principal stress axes and
of the intermediate principal stress on the shear strength. In Proc. 6th Int. Conf.
SMFE, volume 1, pages 179183, 1965.
[23] J.B. Burland. Contribution to discussion on session 4. In Proc. 7th Eur. Conf. Soil
Mech., volume 4, page 137, Brighton, 1979. SMFE.
[24] J.B. Burland. Small is beautiful - the stiffness of soils at small strains, Ninth Laurits
Bjerrum memorial lecture. Canadian Geotechnical Journal, 26:499516, 1989.
[25] J.B. Burland, B. Simpson, and H.D. St. John. Movements around excavations in
London clay. In Proc. 7th Eur. Conf. Soil Mech., volume 1, pages 1330, Brighton,
1979. SMFE.
[26] T.K. Caughey. Sinusoidal excitation of a system with bilinear hysteresis. Journal of
Applied Mechanics, ASME, 27(4):640643, 1960.
[27] C.S. Chang. Micromechanical modelling of constitutive relations for granular material. In M. Satake and J.T. Jenkins, editors, Micromechanics of Granular Materials,
pages 271278, Amsterdam, The Netherlands, 1978. Elsevier.
152
Bibliography
[28] S.K. Chaudhary, J. Kuwano, and Y. Hayano. Measurement of quasi-elastic stiffness parameters of dense Toyoura sand in hollow cylinder apparatus and triaxial
apparatus with bender elements. Geotechnical Testing Journal, 27(1):113, 2004.
[29] K.W. Cole and J.B. Burland. Observation of retaining wall movements associated
with a large excavation. In Proc. 5th ECSMFE, volume 1, page 445453, Madrid,
1972.
[30] M.B. Darendeli and K.H. Stokoe. Development of a new family of normalized
modulus reduction and material damping curves. Engrg. Rpt. GD01-1, University
of Texas, Austin, Texas, 2001.
[31] J. Desrues, P.A. Vermeer, and B. Zweschper. Database for tests on Hostun RF sand.
Institutsbericht 13, Universitat Stuttgart, 2000.
[32] T. Doanh. Constitutive Modelling of Granular Materials, chapter Strain resonse envelope: a complementary tool for evaluating hypoplastic constitutive equations,
pages 375396. Berlin Heidelberg: Springer, 2000.
[33] D.C. Drucker. On uniqueness in the theory of plasticity. Quart. Appl. Math., XIV:
3542, 1956.
[34] D.C. Drucker. Introduction to Mechanics of Deformable Solids. McGraw-Hill, 1967.
[35] D.C. Drucker and W. Prager. Soil mechanics and plastic analysis or limit design.
Quart. Appl. Math., 10(2):157165, 1952.
[36] J.M. Duncan and C.-Y. Chang. Nonlinear analysis of stress and strain in soil. Proc.
ASCE: Journal of the Soil Mechanics and Foundation Division, 96:16291653, 1970.
[37] R. Dyvik and C. Madshus. Lab measurements of Gmax using bender elements. In
Advances in the Art of Testing Soils Under Cyclic Conditions, pages 186196., Detroit,
Michigan, 1985. ASCE. ASCE convention.
[38] W. Ehlers. A single-surface yield function for geomaterials. Arch. Appl. Mech., 65:
246259, 1995.
[39] A.L. Fernandez and J.C. Santamarina. Effect of cementation on the small-strain
parameters of sands. Can. Geotech. J., 38:191199, 2001.
[40] V. Fioravante. Anisotropy of small-strain stiffness of Ticino and Kenya sands from
seismic wave propagation measured in triaxial testing. Soils and Foundations, 40(4):
129142, 2000.
[41] E. Flavigny, J. Desrues, and B. Palayer. Le sable dHostun RF. Rev. Franc. Geotechn.,
53:6770, 1990.
153
Bibliography
[42] D. Forel, T. Benz, and W.D. Pennington. Seismic Data Processing with Seismic Un*x.
Number 12 in Course Notes. Society of Exploration Geophysicists, Houston, TX,
2005.
Normung (DIN). DIN EN ISO 14688. Berlin: Beuth Verlag,
[43] Deutsches Institut fur
11 2004. Geotechnical investigation and testing. Identification and classification of
soil.
[44] V.N. Georgiannou, S. Rampello, and F. Silvestri. Static and dynamic measurements
of undrained stiffness on natural over consolidated clays. In Proc. 10th ECSMFE,
volume 1, pages 9195, Florence, 1991.
[45] J.D. Goddard. Nonlinear elasticity and pressure-dependent wave speeds in granular media. Proceedings of the Royal Society London, (430):105131, 1990.
[46] A. Grammatikopoulou. Development, implementation and application of kinematic
hardening models for overconsolidated clays. PhD thesis, University of London (Imperial College), 2004.
[47] G. Gudehus. A comparison of some constitutive laws for soils under radially symmetric loading and unloading. In Proc. of 3rd Num. Meth. in Geomechanics, volume 4,
pages 13091323. A.A. Balkema, 1979.
[48] M.J. Gunn. The prediction of surface settlement profiles due to tunnelling. In
G.T. Houlsby and A.N. Schofield, editors, Predictive Soil Mechanics, pages 304316.
Thomas Telford, 1993. Proc. Wroth Memorial Symposium.
[49] B.O. Hardin. The nature of stress-strain behaviour of soils. In Proc. Earthquake
Engineering and Soil Dynamics, volume 1, pages 390, Pasadena, CA, 1978. ASCE,
New York. State-of-the-art report.
[50] B.O. Hardin and W.L. Black. Vibration modulus of normally consolidated clays.
Proc. ASCE: Journal of the Soil Mechanics and Foundations Division, 94(SM2):353369,
1968.
[51] B.O. Hardin and W.L. Black. Closure to vibration modulus of normally consolidated clays. Proc. ASCE: Journal of the Soil Mechanics and Foundations Division, 95
(SM6):15311537, 1969.
[52] B.O. Hardin and V.P. Drnevich. Shear modulus and damping in soils: Design
equations and curves. Proc. ASCE: Journal of the Soil Mechanics and Foundations
Division, 98(SM7):667692, 1972.
[53] B.O. Hardin and V.P. Drnevich.
Shear modulus and damping in soils:
Measurement and parameter effects. Proc. ASCE: Journal of the Soil Mechanics and
Foundations Division, 98(SM6):603624, 1972.
154
Bibliography
[54] B.O. Hardin and F.E. Richart Jr. Elastic wave velocities in granular soils. Proc.
ASCE: Journal of the Soil Mechanics and Foundations Division, 89(SM1):3365, 1963.
[55] K. Hashiguchi. Fundamental requirements and formulations of elastoplastic constitutive equations with tangential plasticity. Int. J. Plasticity, 9(5):525549, 1993.
[56] P. Hepton. Shear wave velocity measurements during penetration testing. In Penetration Testing in the UK, pages 275283, Birmingham, UK, 1989. Thomas Telford,
London. Proceedings of the Geotechnology Conference on Penetration Testing in
the United Kingdom (PTUK).
[57] P.-Y. Hicher. Elastic properties of soils. Journal of Geotechnical Engineering, 122(8):
641648, 1996.
[58] R. Hill. The Mathematical Theory of Plasticity. Oxford University Press, London,
1950.
[59] U. Holzlohner.
Dynamische Bodenkennwerte - Meergebnisse und Zusammenhange. Bautechnik, 65(9):306312, 1988.
[60] E. Hoque and F. Tatsuoka. Kinematic elasticity of a granular material. In Proc.
GeoEng2000, Melbourne, 2000. GEE0085.PDF.
[61] E. Hoque and F. Tatsuoka. Effects of stress ratio on small-strain stiffness during
triaxial shearing. Geotechnique, 54(7):429439, 2004.
[62] G.T. Houlsby. Critical state models and small-strain stiffness. In Developments in
Theoretical Geomechanics, pages 295312. Balkema, 2000. Proc. of the Booker Memorial Symposium (Sydney).
[63] G.T. Houlsby and C.P. Wroth. The variation of shear modulus of a clay with pressure and overconsolidation ratio. Soils and Foundations, 31(3):138134, 1991.
[64] C-C. Hsu and M. Vucetic. Dynamic and cyclic behavior of soils over a wide range
of shear strains in ngi-type simple shear testing device. UCLA Reserarch Report
ENG-02-228, Civil and Environmental Engineering Department, University of California, Los Angeles, 2002.
[65] T. Hueckel and R. Nova. Some hysteresis effects of the behavior of geological
media. International Journal of Solids and Structures, 15(8):625642, 1979.
[66] T.J.R. Hughes. The Finite Element Method; Linear Static and Dynamic Finite Element
Analysis. Prentice-Hall, 1987.
[67] T.J.R. Hughes and R.L. Taylor. Unconditionally stable algorithms for quasi-static
elasto/viscoplastic finite element analysis. Computers and Structures, 8:169173,
1978.
155
Bibliography
[68] I.M. Idriss and H.B. Seed. Seismic response of horizontal soil layers. Proc. ASCE:
Journal of the Soil Mechanics and Foundation Division, 94(SM4):10031031, 1968.
[69] T. Imai and K. Tonouchi. Correlation of N-nalue with S-wave velocity and shear
modulus. In Proc. 2nd European Symposium on Penetration Testing, pages 6672,
Amsterdam, 1982.
[70] I. Ishibashi and X. Zhang. Unified dynamic shear moduli and damping ratios of
sand and clay. Soils and Foundations, 33(1):182191, 1993.
[71] K. Ishihara, F. Tatsuoka, and S. Yasuda. Undrained deformation and liquefaction
of sand under cyclic stresses. Soils and Foundations, 15(1):2944, 1975.
[72] J.N. Israelachvili. Intermolecular and surface forces. Academic Press, 1992.
[73] W.D. Iwan. On a class of models for the yielding behaviour of continuous and
composite systems. Journal of Applied Mechanics, ASME, 34(3):612617, 1967.
[74] T. Iwasaki and F. Tatsuoka. Effects of grain size and grading on dynamic shear
moduli of sands. Soils and Foundations, 17(3):1935, 1977.
[75] T. Iwasaki, F. Tatsuoka, and Y. Takagi. Hysteretic damping of sands under cyclic
loading and its relation to shear modulus. Soils and Foundations, 18(2):2540, 1978.
[76] T. Iwasaki, F. Tatsuoka, and Y. Takagi. Shear moduli of sands under cyclic torsional
shear loading. Soils and Foundations, 18(1):3956, 1978.
[77] N. Janbu. Soil compressibility as determined by oedometer and triaxial tests. In
Proc. 3rd ECSMFE, volume 1, pages 1925, Wiesbaden, 1963.
[78] R.J. Jardine, A.B. Fourie, J. Maswoswe, and J.B. Burland. Field and labaratory
measurement of soil stiffness. In Proc. 11th Int. Conf. SMFE, volume 2, pages 511
514, San Fransisco, 1985.
[79] R.J. Jardine, D.M. Potts, A.B. Fourie, and J.B. Burland. Studies of the influence of
non-linear stress-strain characteristics in soil-structure interaction. Geotechnique,
36(3):377396, 1986.
[80] R.J. Jardine, M.J. Symes, and J.B. Burland. The measurement of soil stiffness in the
triaxial apparatus. Geotechnique, 34(3):323340, 1984.
[81] P.C. Jennings. Periodic response of a general yielding structure. Proc. ASCE: Journal
of the Engineering Mechanics Division, 90(EM2):131166, 1964.
[82] P.C. Jennings. Earthquake response of a yielding structure. Proc. ASCE: Journal of
the Engineering Mechanics Division, 91(EM4):4167, 1965.
[83] B. Jeremic and S. Sture. Implicit integrations in elastoplastic geotechnics. Mechanics
of Cohesive-Frictional Materials, 2:165183, 1997.
156
Bibliography
[84] H.D. St John. Field and theoretical studies of behaviour of ground around deep excavations
in London Clay. PhD thesis, University of Cambridge, 1975.
[85] K.L. Johnson. Contact Mechanics. Cambridge University Press, 1985.
[86] V. Jovicic and M.R. Coop. Stiffness of coarse-grained soils at small strains.
Geotechnique, 47(3):545561, 1997.
[87] F.E. Richart Jr., R.D. Woods, and J.R. Hall Jr. Vibrations of Soils and Foundations.
Civil Engineering and Engineering Mechanics Series. Prentice-Hall, Inc., Englewood Cliffs, N.J., 1970.
[88] T.T. Kallioglou, A. Papadopoulou, and K. Pitilakis. Shear modulus and damping
of natural sands. In H. Di Benedetto, T. Doanh, H. Geoffroy, and C. Sauzeat, editors, Deformation Characteristics of Geomaterials, pages 401407. A.A. Balkema, 2003.
Proc. IS-Lyon03.
[89] Stokoe K.H., M.B. Darendeli, R.B. Gilbert1, F.-Y. Menq, and W.K. Choi. Development of a new family of normalized modulus reduction and material damping
curves. In International Workshop on Uncertainties in Nonlinear Soil Properties and
their Impact on Modeling Dynamic Soil Response, UC Berkeley, CA, 2004.
[90] T.C. Kim and M. Novak. Dynamic properties of some cohesive soils of ontario.
Canadian Geotech. Journal, 18(3):371389, 1981.
[91] W.T. Koiter. General theorems for elastic-plastic solids. In I.N. Sneddon and R. Hill,
editors, Progress in Solid Mechanics, volume 1, pages 165221, North-Holland, Amsterdam, 1960.
[92] T. Kokusho, Y. Yoshida, and E. Yasuyuki. Dynamic properties of soft clay for wide
strain range. Soils and Foundations, 22(4):118, 1982.
[93] D. Kolymbas. An outline of hypoplasticity. Archive of Applied Mechanics, 61:143
151, 1991.
[94] R.L. Kondner and J.S. Zelasko. A hyperbolic stress-strain formulation for sands.
In 2nd Pan. Am. Conf. Soil Mech. Found. Eng., volume 1, pages 289394, Brazil, 1963.
[95] W. Krajewski and O. Reul. Erfahrungen zur numerischen Berechnung von
157
Bibliography
[98] P.V. Lade and V. Abelev. Characterization of cross-anisotropic soil deposits from
isotropic compression tests. Soils and Foundations, 45(5):89102, 2005.
[99] P.V. Lade and J.M. Duncan. Elasto-plastic stress-strain theory for cohesionless soil.
Proc. ASCE: Journal of the Geotechnical Engineering Division, 101(GT10):10371053,
1975.
[100] R. Larsson and M. Mulabdic. Piezocone tests in clay. Technical Report 42, Swedish
158
Bibliography
[113] H. Matsuoka. Stress-strain relationships of sands based on the mobilized plane.
Soils and Foundations, 14(2):4761, 1974.
[114] H. Matsuoka and T. Nakai. A new failue criterion for soils in three dimensional
stresses. In IUTAM Conference on Deformation and Failure of Granular Materials,
pages 253263, Delft, 1982.
[115] P.W. Mayne and F.H. Kulhawy. K0 -OCR relationships in soil. Journal of Geotechnical
Engineering, 108(6):851872, 1982.
[116] P.W. Mayne and G.J. Rix. Gmax qc relationships for clays. Geotechnical Testing
Journal, 16(1):5460, 1993.
[117] G. Mesri, T.W. Feng, and J.M. Benak. Postdensification penetration resistance of
clean sands. Journal of Geotechnical Engineering, ASCE, 116(7):10951115, 1990.
[118] R.D. Mindlin. Compliance of elastic bodies in contact. Journal of Applied Mechanics,
ASME, 16:259268, 1949.
[119] R.D. Mindlin and H. Deresiewicz. Elastic spheres in contact under varying oblique
forces. Journal of Applied Mechanics, ASME, 20(3):327344, 1953.
[120] J.K. Mitchell. Fundamentals of Soil Behavior. John Wiley & Sons, 2 edition, 1993.
159
Bibliography
[128] T. Nakai and H. Matsuoka. Deformations of soils in three-dimensional stresses. In
IUTAM Conference on Deformation and Failure of Granular Materials, pages 275285,
Delft, 1982.
[129] A. Niemunis and I. Herle. Hypoplastic model for cohesionless soils with elastic
strain range. Mechanics of Cohesive-Frictional Materials, 2(4):279299, 1997.
[130] R.L. Nigbor and T. Imai. The suspension p-s velocity logging method. In Proc. 13th
Int. Conf. SMFE, volume 1, pages 5761, New Delhi, India, 1994. TC10.
160
Bibliography
[143] R. Pyke. Nonlinear soil models for irregular cyclic loadings. Proc. ASCE: Journal of
the Geotechnical Engineering Division, 105(GT6):715726, 1979.
[144] T. Ramamurthy and P.C. Rawat. Shear strength of sand under general stress system. In Proc. 8th Int. Conf. SMFE, volume 1, pages 339342, Moscow, 1973.
[145] W. Ramberg and W. R. Osgood. Description of stress-strain curve by three parameters. Technical Note 902, National Advisory Committee for Aeronautics, Washington, DC, 1943.
[146] K.I. Rammah, D.V. Val, and A.M. Puzrin. Effects of ageing on small-strain stiffness
of overconsolidated clays. Geotechnique, 54(5):319322, 2004.
[147] S. Rampello, G.M.B. Viggiani, and A. Amorosi. Small-strain stiffness of reconstituted clay compressed along constant triaxial effective stress ratio paths.
Geotechnique, 47(3):475489, 1997.
[148] D. Richardson. Investigations of threshold effects in soil deformation. PhD thesis, City
University, London, 1988.
[149] F.E. Jr. Richart. Some effects of dynamic soil properties on soil structure interaction.
Proc. ASCE: Journal of the Geotechnical Engineering Division, 101(GT12):11931240,
1975.
[150] P.K. Robertson, R.G. Campanella, D. Gillespie, and A. Rice. Seismic cpt to measure in-situ shear wave velocity. Journal of Geotechnical Engineering, 112(8):791803,
1986.
[151] K.H. Roscoe and A.N. Schofield. Mechanical behaviour of an idealised wet clay.
In Proc. 2nd European Conf. SMFE, volume 1, pages 4754, Wiesbaden, 1963.
[152] E. Rosenblueth and I. Herrera. On a kind of hysteretic damping. Proc. ASCE:
Journal of the Engineering Mechanics Division, 90(EM4):3747, 1964.
[153] P.W. Rowe. The stress-dilatancy relation for static equilibrium of an assembly of
particles in contact. In Proc. of the Royal Society of London. Series A, Mathematical and
Physical Sciences, volume 269, pages 500527, 1962.
[154] K. Runesson and A. Samuelsson. Numeta 85 Numerical Methods in Engineering,
Theory and Applications, chapter Aspects on numerical techniques in small deformation plasticity, pages 337347. A.A. Balkema, 1985.
[155] J.A. Santos and A.G. Correia. Reference threshold shear strain of soil. Its application to obtain a unique strain-dependent shear modulus curve for soil. In 15th Int.
Conf. SMGE, volume 1, pages 267270, Istanbul, 2001. A.A. Balkema.
[156] S.A. Savidis. Empfehlungen des Arbeitskreises Baugrunddynamik. Grundbauinstitut der Technischen Universitat Berlin, 2002.
161
Bibliography
[157] T. Schanz. Zur Modellierung des mechanischen Verhaltens von Reibungs Geotechnik 45, Universitat Stuttgart, 1998.
materialien. Mitt. Inst. fur
[158] T. Schanz, P.A. Vermeer, and P.G. Bonnier. Beyond 2000 in Computational Geotechnics,
chapter Formulation and verification of the Hardening-Soil Model, pages 281290.
Balkema, Rotterdam, 1999.
[159] J.H. Schmertmann. The mechanical aging of soils. Journal of Geotechnical Engineering, 117(9):12881330, 1991.
[160] B.F. Schmidt. Settlements and ground movements associated with tunnelling in soils.
PhD thesis, University of Illinois, Urbana, 1969.
[161] H.F. Schweiger. Results from numerical benchmark exercises in geotechnics. In
P. Mestat, editor, Numge 2002. 5th European Conference Numerical Methods in Geotechnical Engineering, volume 1, pages 305314, Paris, 2002. LCPC.
[162] H.B. Seed and I.M. Idriss. Soil moduli and damping factors for dynamic response
analysis. Report 70-10, EERC, Berkeley, CA, 1970.
[163] T. Shibata and D. Karube. Influence of the variation of the intermediate principal
stress on mechanical properties of normally consolidated clay. In Proc. 6th Int. Conf.
SMFE, volume 1, pages 359363, Paris, 1985.
[164] S. Shibuya, T. Mitachi, S. Yamashita, and H. Tanaka. Effect of sample disturbance
on Gmax of soils - a case study. In S. Shibuya T. Mitachi and S. Miura, editors,
Prefailure Deformations of Geomaterials, volume 2, pages 7782. Balkema, Rotterdam,
1995.
[165] D.J. Shirley and L.D. Hampton. Shear-wave measurements in laboratory sediments. The Journal of the Acoustical Society of America, 63(2):607613, 1978.
[166] J.C. Simo and T.J.R. Hughes. Computational Inelasticity, volume 7 of Interdisciplinary
Applied Mathematics. Springer, 1998.
[167] J.C. Simo and R.L. Taylor. Consistent tangent operators for rate-independent
elastoplasticity. Computer Methods in Applied Mechanics and Engineering, 48:101
118, 1985.
[168] B. Simpson. Retaining structures: displacement and design, the 32nd Rankine
lecture. Geotechnique, 42(4):541576, 1992.
[169] O.K. Soreide. Mixed hardening models for frictional soils. PhD thesis, Norwegian
University of Science and Technology (NTNU), Trondheim, 1990.
[170] S.A. Stallebrass. Modelling the effect of recent stress history on the deformation of overconsolidated soils. PhD thesis, City University, London, 1990.
162
Bibliography
[171] S.E. Stallebrass and R.N. Taylor. The development and evaluation of a constitutive model for the prediction of ground movements in overconsolidated clay.
Geotechnique, 47(2):235253, 1997.
[172] K.H. Stokoe, M.B. Darendeli, R.D. Andrus, and L.T. Brown. Dynamic soil properties: laboratory, field and correlation studies. In Proc. 2nd Int. Conf. on Earthquake
Geotech. Eng., volume 3, pages 811845. A.A. Balkema, 1999.
[173] K.H. Stokoe and J.C. Santamarina. Seismic-wave-based testing in geotechnical engineering. In GeoEng 2000: An International Conference on Geotechnical and Geological
Engineering, volume 1, pages 14901536, Melbourne, Austtralia, 2000. Technomic
Publishing Company.
[174] K.H. Stokoe, S.G. Wright, A.B. James, and M.R. Jose. Geophysical characterization of
sites, chapter Characterization of geotechnical sites by SASW method, pages 1525.
Balkema, A.A., Rotterdam, 1994.
[175] K.H.I. Stokoe and R.D. Woods. In situ shear wave velocity by cross-hole method.
Proc. ASCE: Journal of the Soil Mechanics and Foundations Division, 98(SM5):443460,
1972.
[176] H.B. Sutherland and M.S. Mesdary. The influence of the intermediate principal
stress on the strength of sand. In Proc. 7th Int. Conf. SMFE, volume 1, pages 391
399, Mexico City, 1969.
[177] F. Tatsuoka. Impacts on geotechnical engineering of several recent findings from
laboratory stress-strain tests on geomaterials. Department of Civil Engineering
Mechanics, Columbia University, 2000. The 2000 Burmister Lecture.
[178] F. Tatsuoka and K. Ishihara. Yielding of sand in triaxial compression. Soils and
Foundations, 14(2):6376, 1974.
[179] K. Terzaghi, R.B. Peck, and G. Mesri. Soil Mechanics in Engineering Practice. WileyInterscience; 3 edition (January 1, 1996), 3 edition, 1996.
[180] S. Toki, S. Shibuya, and S. Yamashita. Standardization of laboratory test methods to determine the cyclic deformation properties of geomaterials in Japan. In
S. Shibuya T. Mitachi and S. Miura, editors, Pre-failure Deformations of Geomaterials,
volume 2, pages 741784. Balkema, Rotterdam, 1995.
[181] C. Truesdell and W. Noll. The Non-Linear Field Theories of Mechanics, volume III of
Encyclopedia of Physics. Berlin, Heidelberg: Springer, 1965.
[182] T.W. and R.V. Whitman. Soil Mechanics. John Wiley & Sons, New York, NY, 1969.
[183] P.A. Vermeer. A double hardening model for sand. Geotechnique, 28(4):413433,
1978.
163
Bibliography
[184] P.A. Vermeer. A modified initial strain method for plasticity problems. In Proc. 3th
Int. Conf. Numer. Meth. Geomech., pages 377387. Balkema, Rotterdam, 1979.
[185] P.A. Vermeer. Formulation and analysis of sand deformation problems. PhD thesis,
Delft University of Science and Technology, 1980.
[186] P.A. Vermeer, T. Benz, and R. Schwab. Modern Trends in Geomechanics, chapter Two
Elastoplastic Models for Small and Large Strains and Their Use in Engineering
Practise, pages 175190. Springer, 2006.
[187] G. Viggiani and J.H. Atkinson. Stiffness of fine grained soils at very small strains.
Geotechnique, 45(2):249265, 1995.
[188] M. Vucetic and R. Dobry. Effect of soil plasticity on cyclic response. Journal of
Geotechnical Engineering, ASCE, 117(1):89107, 1991.
[189] M. Vucetic and K. Tabata. Influence of soil type on the effect of strain rate on
small-strain cyclic shear modulus. Soils and Foundations, 43(5):161173, 2003.
[190] M. Vucetic, K. Tabata, and L. Matesic. Effect of average straining rate on shear
modulus at small cyclic strains. In H. Di Benedetto, T. Doanh, H. Geoffroy, and
C. Sauzeat, editors, Deformation Characteristics of Geomaterials, pages 321328. A.A.
Balkema, 2003. Proc IS-Lyon03.
[191] M. Wehnert. Ein Beitrag zur drainierten und undrainierten Analyse in der Geotechnik.
PhD thesis, Universitat Stuttgart, 2006.
[192] R.V. Whitman. The behaviour of soils under transient loadings. In Proc. 4th Int.
Conf. SMFE, volume 1, pages 207210, London, 1957.
[193] T. Wichtmann and T. Triantafyllidis. Influence of a cyclic and dynamic loading
history on dynamic properties of dry sand, Part i: cyclic and dynamic torsional
prestraining. Soil Dynamics and Earthquake Engineering, 24(2):127147, 2004.
[194] R.D. Woods. Measurement of dynamic soil properties. In Proc. Earthquake Engineering and Soil Dynamics, volume 1, pages 9178, Pasadena, CA, 1978. ASCE, New
York.
[195] C.P. Wroth. In-situ measurement of initial stresses and deformation characteristics.
In Proc. Geot. Engng Div. Speciality Conf. on In-situ Measurement of Soil Properties,
pages 181230, Rayleigh, NC, 1975. ASCE.
[196] Y. Yamada and K. Ishihara. Anisotropic deformation characteristic of sand under
three dimensional stress conditions. Soils and Foundations, 19(1):7994, 1979.
[197] S. Yimsiri and K. Soga. Effect of surface roughness on small-strain modulus: Micromechanics view. In M. Jamiolkowski, R. Lancellotta, and D. Lo Presti, editors,
Pre-failure Deformation Characteristics of Geomaterials, volume 1, pages 597602. A.A.
Balkema, Rotterdam, 1999.
164
Bibliography
[198] R.N. Yong and R.D. Japp. A flow law for clays in dynamic compression. In
Proc. Int. Symp. Wave Propagation and Dynam. Prop. of Earth Materials, Albuquerque,
N.M., 1967. University of New Mexico Press.
[199] M.S. Youssef, A.H. El Ramli, and M. El Demery. Relationships between shear
strength, consolidation, liquid limit, and plastic limit for remoulded clays. In Proc.
6th Int. Conf. SMFE, volume 1, pages 126129, Montreal, 1965.
[200] P. Yu and Jr. F.E. Richart. Stress ratio effects on shear modulus of dry sands. Journal
of Geotechnical Engineering, 110(3):331345, 1984.
[201] C. Zambelli. Modellazione constitutiva del comportamento meccanico in campo dinamico
di material granulari. Tesi di laurea, Politecnico di Milano, 2002.
[202] L. Zdravkovic and R.J. Jardine. Some anisotropic stiffness characteristics of a silt
under general stress conditions. Geotechnique, 47(3):407437, 1997.
165
Appendix A
Small-Strain Overly Fortran Code
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
!+
! 23456789012345678901234567890123456789012345678901234567890123456789|
!+
! Small S t r a i n O v e r l a y
|
!+
! P u r p o s e : Add s m a l l s t r a i n s t i f f n e s s t o e l a s t o p l a s t i c m o d e l s
|
!
|
!
I /O Type
|
! Eps0
I R ( 6 ) : t o t a l train at s t a r t of increment
|
! dEps
I R ( 6 ) : s t r a i n increment
|
! UsrSta0 I R ( ) : s t a t e v a r i a b l e s a t s t a r t o f i n c ; u s e s # 2 t o # 1 0 |
! U s r S t a 0 O R ( ) : s t a t e v a r i a b l e s a t end o f i n c ; u s e s # 2 t o # 1 0
|
! Sig0
I R ( 6 ) : s t r e s s at s t a r t of increment
|
! G
I R
: e l a s t i c modulus o f b a c k g r o u n d m o d e l ( l a r g e r s t r . ) |
! GFac
I R
: r a t i o G0 /G w h e r e G0 i s t h e s m a l l s t r a i n modulus |
! Gam07
I R
: 70% threshold shear strain
|
! GAct
O R
: A c t u a l e l a s t i c modulus t o b e u s e d f o r c a l c u l a t i o n |
!
|
! The s t a t e v a r i a b l e s a r e o r g a n i z e d a s f o l l o w s
|
! UsrSta ( 1 ) : not used
|
! UsrSta ( 2 ) : not used
|
! U s r S t a ( 3 ) : S t r a i n h i s t o r y component 1
|
! U s r S t a ( 4 ) : S t r a i n h i s t o r y component 2
|
! U s r S t a ( 5 ) : S t r a i n h i s t o r y component 3
|
! U s r S t a ( 6 ) : S t r a i n h i s t o r y component 4
|
! U s r S t a ( 7 ) : S t r a i n h i s t o r y component 5
|
! U s r S t a ( 8 ) : S t r a i n h i s t o r y component 6
|
! U s r S t a ( 9 ) : A c t u a l s t i f f n e s s GAct f o r o u t p u t e t c . / r e d u n d a n t
|
! U s r S t a ( 1 0 ) : Monoton ic s t i f f n e s s ; r e l o a d i n g f o r GAct>G UsrSt ( 1 0 )
|
!
|
! Needs : S u b r o u t i n e P r n S i g ( dEpsD , xN1e , xN2e , xN3e , dEpsD1 , dEpsD2 , dEpsD3 ) |
!
t h a t f i n d s t h e E i g e n s y s t e m ( xN1e , xN2e , xN3e ) and E i g e n v a l u e s |
!
( dEpsD1 , dEpsD2 , dEpsD3 ) o f t h e s t r a i n i n c r e m e n t dEpsD
|
!+
! Author : Thomas Benz < t h o m a s . benz@wewinumge . de >
|
!+
!
SUBROUTINE S m a l l S t r a i n ( Eps0 , dEps , UsrSta , UsrSta0 , Sig0 ,
167
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
168
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
2
3
4
141
142
143
144
145
146
T r a f o ( 1 , 3 ) = xN3e ( 1 )
T r a f o ( 2 , 3 ) = xN3e ( 2 )
T r a f o ( 3 , 3 ) = xN3e ( 3 )
T r a n s f o r m s t a r i n h i s t o r y i n p r i n c i p a l s y s t e m o f dEpsD
DO i = 1 , 3
DO j = 1 , 3
X1 = 0 d0
DO k = 1 , 3
X2 = 0 d0
DO l = 1 , 3
X2 = X2 + H i s t 1 ( k , l ) T r a f o ( l , j )
END DO
X1 = X1 + X2 T r a f o ( k , i )
END DO
H i s t 2 ( i , j ) = X1
END DO
END DO
Elongation of uncorrected history
DO i = 1 , 3
Chi ( 1 , i ) = H i s t 2 ( i , i ) 1 0 0 d0
END DO
Build transformation matrix
nCountXT = 0
DO i = 1 , 3
xT ( i ) = 1 d0
IF ( c h i ( 1 , i ) dEpsD ( i ) . LT . 0 d0 ) THEN
xT ( i ) = 1 0 d0/DSQRT( c h i ( 1 , i ) + 1 0 0 d0 )
nCountXT = nCountXT + 1
END IF
END DO
C a l c u l a t e new h i s t o r y
DO i = 1 , 3
DO j = 1 , 3
H i s t 2 ( i , j ) = H i s t 2 ( i , j ) xT ( i ) xT ( j )
END DO
Chi ( 2 , i )
= H i s t 2 ( i , i ) 1 0 0 d0
H i s t 2 ( i , i ) = H i s t 2 ( i , i ) + dEpsD ( i )
Chi ( 3 , i )
= H i s t 2 ( i , i ) 1 0 0 d0
END DO
Calculate history strain ( iso )
dEps2 = ( dEpsD ( 1 ) 2 + dEpsD ( 2 ) 2 + dEpsD ( 3 ) 2 )
DO i = 1 , 2
GamHis ( i ) = DSQRT( 3 d0
( ( dEpsD ( 1 ) Chi ( i + 1 , 1 ) ) 2 +
( dEpsD ( 2 ) Chi ( i + 1 , 2 ) ) 2 +
( dEpsD ( 3 ) Chi ( i + 1 , 3 ) ) 2 )/ dEps2 )
END DO
From s t r a i n h i s t o r y t o s t i f f n e s s
Fac0 = GamHis(2) GamHis ( 1 )
Fac1 = Gam07/ 0 . 3 8 5 d0 ! u s e 0 . 4 2 9 f o r HD
DO i = 1 , 2 ! i =1:= l o a d i n g ; i = 2 : = r e l o a d i n g
IF ( i . EQ . 2 ) Fac1 = 2 d0 Fac1
169
147
148
149
150
151
152
153
154
155
156
157
158
2
3
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
!
!
188
189
190
191
192
193
194
195
196
197
198
IF ( ( UsrSta0 ( 1 0 ) . NE . 1 ) . OR . ( i . EQ . 2 ) ) THEN
GamOff = ( DSQRT( GFac)1d0 ) Fac1
IF ( Fac0 . LE . 0 ) THEN
GBct ( i ) = GGFac ( 1 / ( 1 +GamHis ( 1 ) / Fac1 )
GamHis ( 1 ) / Fac1 /(1+GamHis ( 1 ) / Fac1 ) 2 )
ELSE
IF ( GamHis ( 1 ) . LT . GamOff .AND. GamHis ( 2 ) . LT . GamOff ) THEN
GBct ( i ) = GGFac ( GamHis ( 2 ) / ( 1 d0+GamHis ( 2 ) / Fac1 )
GamHis ( 1 ) / ( 1 d0+GamHis ( 1 ) / Fac1 ) ) / Fac0
ELSE
IF ( GamHis ( 1 ) . LT . GamOff .AND. GamHis ( 2 ) . GT . GamOff ) THEN
GBct ( i ) = ( GGFac ( GamOff/(1 d0+GamOff/Fac1 )
GamHis ( 1 ) / ( 1 d0+GamHis ( 1 ) / Fac1 ) ) +
G ( GamHis(2) GamOff ) ) / Fac0
ELSE
GBct ( i ) = G
END IF
END IF
END IF
END IF
END DO
Update S t a t e V a r w i t h new h i s t o r y ( B a c k T r a f o )
DO i = 1 , 3
DO j = 1 , 3
X1 = 0 d0
DO k = 1 , 3
X2 = 0 d0
DO l = 1 , 3
X2 = X2 + H i s t 2 ( k , l ) T r a f o ( j , l )
END DO
X1 = X1 + X2 T r a f o ( i , k )
END DO
H i s t 1 ( i , j ) = X1
END DO
END DO
UsrSta ( 3 ) = H i s t 1 ( 1 , 1 )
UsrSta ( 4 ) = H i s t 1 ( 2 , 2 )
UsrSta ( 5 ) = H i s t 1 ( 3 , 3 )
UsrSta ( 6 ) = H i s t 1 ( 1 , 2 ) + H i s t 1 ( 2 , 1 )
UsrSta ( 7 ) = H i s t 1 ( 2 , 3 ) + H i s t 1 ( 3 , 2 )
UsrSta ( 8 ) = H i s t 1 ( 1 , 3 ) + H i s t 1 ( 3 , 1 )
L o a d i n g / r e l o a d i n g i n c o m b i n a t i o n w i t h an e l a s t o p l a s t i c m o d e l :
Update o f U s r S t a ( 1 0 )
IF ( UsrSta0 ( 1 0 ) .NE . 1 . ) UsrSta ( 1 0 ) =MIN( GBct ( 1 ) /G, UsrSta0 ( 1 0 ) )
Copy s t i f f n e s s ( a l w a y s r e l o a d i n g i n c o m b i n a t i o n w i t h p l a s t i c i t y )
GAct = GBct ( 2 )
END IF
Plateau avoids numerical problems in i s o t r o p i c loading
I f ( GamHis ( 2 ) . Lt . 1 d 6 ) GAct = GGFac
!
!
! E i t h e r RETURN h e r e o r c h e c k s t r e s s work . . .
170
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
!
Do i =1 ,6
Eps ( i ) = Eps0 ( i ) + dEps ( i )
End Do
EpsV0 =( Eps0 ( 1 ) + Eps0 ( 2 ) + Eps0 ( 3 ) )
EpsD0 = ( ( Eps0 (1) Eps0 ( 2 ) ) 2
2
+( Eps0 (2) Eps0 ( 3 ) ) 2
3
+( Eps0 (3) Eps0 ( 1 ) ) 2 ) / 6 d0
4
+( Eps0 ( 4 ) / 2 d0 ) 2 + ( Eps0 ( 5 ) / 2 d0 ) 2 + ( Eps0 ( 6 ) / 2 d0 ) 2
IF ( DABS( EpsV0 ) . Lt . 1 d 1 2 .OR.
2
DABS( EpsD0 ) . Lt . 1 d12
) THEN
! No Check o r a l t e r n a t i v e c o d e
ELSE
SigV0 =( S i g 0 ( 1 ) + S i g 0 ( 2 ) + S i g 0 ( 3 ) ) / 3 d0
SigD0 = ( ( S i g 0 (1) S i g 0 ( 2 ) ) 2
2
+( S i g 0 (2) S i g 0 ( 3 ) ) 2
3
+( S i g 0 (3) S i g 0 ( 1 ) ) 2 ) / 6 d0
4
+( S i g 0 ( 4 ) ) 2 + ( S i g 0 ( 5 ) ) 2 + ( S i g 0 ( 6 ) ) 2
xxK=SigV0/EpsV0
xxG=DSQRT( SigD0/EpsD0 )/2 d0
xxF1=xxK + 4 . d0xxG / 3 . d0
xxF2=xxK2.d0xxG / 3 . d0
DEl=0
DO i =1 ,3
DO j =1 ,3
DEl ( i , j ) = xxF2
END DO
DEl ( i , i ) = xxF1
DEl ( i +3 , i + 3 ) = xxG
END DO
xWElasticD2 = DOT PRODUCT(MATMUL( DEl , Eps ) , Eps )
2
DOT PRODUCT(MATMUL( DEl , Eps0 ) , Eps0 )
END IF
IF ( xWElasticD2 . LT . 0 d0 ) THEN ! U n l o a d i n g
GAct=MAX( GAct , xxG )
! Correction
END IF
! For output e t c .
UsrSta ( 9 ) = GAct
RETURN
END SUBROUTINE S m a l l S t r a i n
171
Appendix B
Return mapping in Fortran code
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
!+
! 23456789012345678901234567890123456789012345678901234567890123456789|
!+
! A cutting plane algotithm a f t e r Boris Jeremic
|
!+
! P u r p o s e : I m p l i c i t r e t u r n mapping
|
!
|
!
I /O Type
|
! D
I R ( 6 , 6 ) : e l a s t i c s t i f f n e s s matrix
|
! Sig0
I R(6)
: s t r e s s at s t a r t of increment
|
! SigE
I R(6)
: elastic trial stress
|
! dEps
I R(6)
: s t r a i n increment
|
! UsrSta0 I R ( )
: s t a t e v a r i a b l e s at s t a r t of increment
|
! UsrSta O R ( )
: s t a t e v a r i a b l e s a t end o f i n c r e m e n t
|
! Sig
O R(6)
: s t r e s s a t end o f i n c r e m e n t
|
! R0 , Q0
I R
: i n t e r n a l v a r i a b l e s at s t a r t of increment
|
! R ,Q
O R
: i n t e r n a l v a r i a b l e s a t end o f i n c r e m e n t
|
! ...
I R or I : model s p e c i f i c p a r a m e t e r s
|
!
|
! Needs : D e r i v a t i v e s o f t h e y i e l d f u n c t i o n ( s ) and p l a s t i c p o t e n t i a l ( s )
! Surface1 ( yield surface ) :
dfdSig :=n
|
! S u r f a c e 1 ( p l a s t i c p o t e n t i a l ) : d g d S i g : =m , d 2 g d S i g 2 : =m2
|
! Surface1 ( int . varaibles ) :
R
|
! Surface2 ( yield surface ) :
dfdSig :=u
|
! S u r f a c e 2 ( p l a s t i c p o t e n t i a l ) : d g d S i g : =w , d 2 g d S i g 2 : = w2
|
! Surface2 ( int . variables ) :
Q
|
! Assumes : dmdR = 0 and dwdQ= 0
|
!+
! Thomas Benz < t h o m a s . benz@wewinumge . de >
|
!+
31
32
33
34
!+
! Case 1 : S i n g l e s u r f a c e r e t u r n
|
!+
35
36
! 1. Initialize
37
38
39
40
41
42
i t =0
itmax =20
xTol = 1 d9
Unity = 0 . 0 d0
Unity ( 1 , 1 ) = 1 . 0 d0
173
43
44
45
46
! 2 . Starting Point
47
48
S i g = SigE
xF1 = f 1 ( Sig , R . . . )
xF1Tol = xTol DAbs ( xF1 )
CALL C a l c D e r i v a t i v e s 1 ( Sig , R , n , n2 ,m, . . . )
xLI1 = xF /(DOT PRODUCT( n ,MATMUL(D,m) ) ( dF1dRdR ) )
S i g = SigExLI1 MATMUL(D,m)
R
= R0 + dR xLI1
xF1 = f 1 ( Sig , R . . . )
49
50
51
52
53
54
55
56
57
! 3 . E u l e r Backward . . . ( I t e r a t i o n )
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
RETURN
END
80
81
82
83
84
85
!+
! C a s e 2 : Two s u r f a c e r e t u r n
|
!+
86
! 1. Initialize
87
88
i t =0
itmax =20
xTol = 1 . 0 d9
Unity = 0 . 0 d0
Unity ( 1 , 1 ) = 1 . 0 d0
Unity ( 2 , 2 ) = 1 . 0 d0
89
90
91
92
93
94
174
95
96
97
Unity ( 3 , 3 ) = 1 . 0 d0
R = R0
Q = Q0
98
99
! 2 . Starting Point
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
S i g = SigE
xF1 = f 1 ( Sig , R . . . )
xF2 = f 2 ( Sig ,Q . . . )
xF1Tol = xTol DAbs ( xF1 )
xF2Tol = xTol DAbs ( xF2 )
CALL C a l c D e r i v a t i v e s 1 ( Sig , R , n , n2 ,m, . . . )
CALL C a l c D e r i v a t i v e s 2 ( Sig , Q, u , u2 ,w , . . . )
xOmega11 = dF1dRdR + DOT PRODUCT( n ,MATMUL(D,m) )
xOmega12 = DOT PRODUCT( n ,MATMUL(D,w) )
xOmega21 = DOT PRODUCT( u ,MATMUL(D,m) )
xOmega22 = dF2dQdQ + DOT PRODUCT( u ,MATMUL(D,w) )
xLI1
= ( xF1
xOmega22 xF2
xOmega12)/
2
( xOmega11 xOmega22 xOmega12 xOmega21 )
xLI2
= ( xF2
xOmega11 xF1
xOmega21)/
2
( xOmega11 xOmega22 xOmega12 xOmega21 )
S i g = SigE ( xLI1 MATMUL(D,m) + xLI2 MATMUL(D,w) )
R = R0 + dR xLI1
Q = Q0 + dQ xLI2
xF1 = f 1 ( Sig , R . . . )
xF2 = f 2 ( Sig ,Q . . . )
121
122
! 3 . E u l e r Backward . . . ( I t e r a t i o n )
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
175
147
+dxLI2 MATMUL(D,w) , Ti )
R = R0 + dR xLI1
Q = Q0 + dQ xLI2
xF1 = f 1 ( Sig , R . . . )
xF2 = f 2 ( Sig ,Q . . . )
i t = i t +1
END DO
148
149
150
151
152
153
154
155
156
157
158
159
160
RETURN
END
161
162
176
Appendix C
Two surface return strategy
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
!+
! 23456789012345678901234567890123456789012345678901234567890123456789|
!+
! HS g r a y c o r n e r s t r a t e g y o r when t o a c t i v a t e c o n e , cap , o r b o t h
|
!+
! After P .G . Bonnier ( 2 0 0 0 )
|
!+
!
1 . I n i t ( Small s t r a i n s t i f f n e s s , D i l a t a n c y , e t c . )
!
2. Trial stress
!
3 . Check y i e l d s u r f a c e s and r e t u r n s t r e s s
IF ( Cone y i e l d s u r f a c e i s v i o l a t e d ) THEN
CALL Cone Return ! R e t u r n and h a r d e n c o n e
IF ( F a i l u r e c r i t e r i a i s v i o l a t e d ) THEN
CALL C o n e F a i l u r e ! R e t u r n t o u l t i m a t e c o n e
IF ( Tension c r i t e r i a i s exceeded ) THEN
CALL Tension CutOff ! Do a 3 s u r f a c e t e n s i o n r e t u r n
END IF
END IF
IF ( Cap y i e l d s u r f a c e i s v i o l a t e d ) THEN
CALL Cap Return ! R e t u r n and h a r d e n c a p
IF ( Cone y i e l d s u r f a c e i s v i o l a t e d ) THEN
CALL Mixed Return ! R e t u r n and h a r d e n c o n e and c a p
IF ( F a i l u r e c r i t e r i a i s v i o l a t e d ) THEN
CALL C a p F a i l u r e ! R e t u r n t o c a p and u l t i m a t e c o n e
END IF
END IF
END IF
ELSE ! No f r i c t i o n h a r d e n i n g ( y e t ) : c h e c k Cap o n l y
IF ( Cap y i e l d s u r f a c e i s v i o l a t e d ) THEN
CALL Cap Return ! R e t u r n and h a r d e n c a p
IF ( Cone y i e l d s u r f a c e i s v i o l a t e d ) THEN
CALL Mixed Return ! R e t u r n and h a r d e n c o n e and c a p
IF ( F a i l u r e c r i t e r i a i s v i o l a t e d ) THEN
CALL C a p F a i l u r e ! R e t u r n t o c a p and u l t i m a t e c o n e
END IF
END IF
END IF
END IF
END IF
!
4 . Check Apex ; Update s t r e s s and s t a t e v a r i a b l e s ; E x i t
177
Appendix D
Material data
Unit
ur
pref
K0nc
Rf
Tension
E0ref
0.7
II. Internal parameters
Eiref
KS /KC
kN
m2
kN
m2
kN
m2
]
[ kN
m2
[ ]
[ ]
[ ]
kN
m2
[ ]
[ ]
kN
m2
kN
m2
[ ]
kN
m2
[ ]
[ ]
Sand(D) Sand(L)
30000
30000
90000
0.55
0.00
42.0
16.0
0.25
100
0.40
0.90
0.00
270000
0.0002
65488
1.47
1.84
Clay
12000
2150
16000
1050
60000 11500
0.75
0.80
0.00
0.00
34.0
20.0
0.0
0.0
0.25
0.20
100
100
0.44
0.66
0.90
0.90
0.00
0.00
168000 80000
0.0001 0.0002
23800
1.56
2.01
14050
0.78
4.76
179
Unit
ur
pref
K0nc
Rf
Tension
E0ref
0.7
II. Internal parameters
Eiref
KS /KC
kN
m3
kN
m3
kN
m2
kN
m2
kN
180
m2
]
[ kN
m2
[ ]
[ ]
[ ]
kN
m2
[ ]
[ ]
kN
m2
kN
m2
[ ]
kN
m2
[ ]
[ ]
Fill
Keuper (up)
Keuper (lo)
Limestone
20.0
10000
10000
30000
0.50
10.00
25.0
0.0
0.20
100
0.58
0.90
0.00
120000
0.0002
24.0
33000
33000
100000
0.40
25.00
25.0
0.0
0.20
100
0.58
0.90
0.00
400000
0.0002
23.0
16000
16000
48000
0.40
25.00
25.0
0.0
0.20
100
0.58
0.90
0.00
192000
0.0002
23.0
190000
190000
575000
0.30
200.00
35.00
0.0
0.20
100
0.43
0.90
0.00
2400000
0.00005
24206
0.80
1.62
81216
0.83
1.72
39353
0.83
1.71
415649
1.52
1.95
Unit
ur
pref
K0nc
Rf
Tension
E0ref
0.7
II. Internal parameters
Eiref
KS /KC
kN
m3
kN
m3
kN
m2
kN
m2
kN
m2
]
[ kN
m2
[ ]
[ ]
[ ]
kN
m2
[ ]
[ ]
kN
m2
kN
m2
[ ]
kN
m2
[ ]
[ ]
Fill
Sand
Clay
17.5
17.5
20.0
20.0
14000
35000
12000
14000
35000
7000
42000 105000
35000
0.50
0.50
0.90
3.00
0.00
7.00
27.0
35.0
31.0
0.0
0.0
0.0
0.20
0.20
0.20
100
100
100
0.55
0.43
0.48
0.90
0.90
0.90
0.00
0.00
0.00
126000 240000 210000
0.0002 0.0002 0.0002
32786
0.88
1.62
72888
1.43
1.64
34582
0.93
1.93
181
182
Parameter (Symbol)
Unit
ur
pref
K0nc
Rf
Tension
E0ref
0.7
II. Internal parameters
Eiref
KS /KC
kN
m3
kN
m3
kN
m2
kN
m2
kN
m2
]
[ kN
m2
[ ]
[ ]
[ ]
kN
m2
[ ]
[ ]
kN
m2
kN
m2
[ ]
kN
m2
[ ]
[ ]
Sand (L1)
Sand (L2)
Sand (L3)
19.0
20.0
45000
45000
180000
0.55
1.00
35.0
5.0
0.20
100
0.43
0.90
0.00
405000
0.0002
19.0
20.0
75000
75000
300000
0.55
1.00
38.0
6.0
0.20
100
0.38
0.90
0.00
675000
0.0002
19.0
20.0
105000
105000
315000
0.55
1.00
38.0
6.0
0.20
100
0.38
0.90
0.00
na
na
96662
1.48
2.15
154447
1.87
2.07
208642
1.88
1.59
Unit
ur
pref
K0nc
Rf
Tension
E0ref
0.7
II. Internal parameters
Eiref
KS /KC
kN
m3
kN
m3
kN
m2
kN
m2
kN
m2
]
[ kN
m2
[ ]
[ ]
[ ]
kN
m2
[ ]
[ ]
kN
m2
kN
m2
[ ]
kN
m2
[ ]
[ ]
Fill
Rupel Clay
18.0
20.0
20000
20000
60000
0.50
1.00
32.0
2.0
0.20
100
0.43
0.90
0.00
180000
0.00015
19.0
20.0
15000
15000
50000
0.70
25.00
20.0
0.0
0.15
100
0.66
0.90
0.00
207000
0.00025
43421
1.16
1.63
39992
5.31
1.36
183
184
Parameter (Symbol)
Unit
ur
pref
K0nc
Rf
Tension
E0ref
0.7
II. Internal parameters
Eiref
KS /KC
kN
m3
kN
m3
kN
m2
kN
m2
kN
m2
]
[ kN
m2
[ ]
[ ]
[ ]
kN
m2
[ ]
[ ]
kN
m2
kN
m2
[ ]
kN
m2
[ ]
[ ]
18.0
20.0
16000
22000
90000
0.50
0.00
37.0
1.5
0.20
100
0.41
0.90
0.00
390000
0.0004
30785
5.08
2.38
30785
5.08
2.38
Table D.7: Sulfeld
excavation - HS-Small parameters
Parameter (Symbol)
Unit
ur
pref
K0nc
Rf
Tension
E0ref
0.7
II. Internal parameters
Eiref
KS /KC
kN
m2
kN
m2
kN
m2
kN
m2
kN
m2
]
[ kN
m2
[ ]
[ ]
[ ]
kN
m2
[ ]
[ ]
kN
m2
kN
m2
[ ]
kN
m2
[ ]
[ ]
18.5
21.5
8500
6150
25750
0.70
6.00
28.0
6.0
0.29
100
0.80
0.90
0.00
154500
0.0003
18.5
21.5
8500
6000
23000
0.90
30.00
32.0
10.0
0.30
100
0.80
0.90
0.00
138000
0.0003
82037
1.62
1.65
27783
0.90
2.78
20715
1.10
2.22
185
Curriculum Vit
Projektingenieur Bruckenbau
Wasserbau, Karlsruhe, Referat Grundbau
2002-2006 Bundesanstalt fur
Wissenschaftlicher Angestellter
Geotechnik
2002-2006 Universitat Stuttgart, Institut fur
Wissenschaftlicher Mitarbeiter
20.09.2006 Promotion
187
Nr. 01
Thamm, B. R.
(1974)
Nr. 02
Gumann, P.
(1975)
Nr. 03
Feeser, V.
(1975)
Nr. 04
Du Thin, K.
(1976)
Nr. 05
Smoltczyk, U. /
Pertschi, O. /
Hilmer, K.
(1976)
Nr. 06
Hilmer, K.
(1976)
Nr. 07
Laumans, Q.
(1977)
Nr. 08
Lachler, W.
(1977)
Nr. 09
Spotka, H.
(1977)
uberdr
ucke
eines normalverdichteten wassergesattigten Tones.
e 5,11
Einheitliche Berechnung von Grundbruch und
Boschungsbruch.
e 2,56
die
Die Bedeutung des Kalziumkarbonats fur
Stutzwand
bei Sand.
Vergriffen
189
Geotechnik
Mitteilungen des Instituts fur
Nr. 10
Schad, H.
(1979)
Nr. 11
Ulrich, G.
(1980)
Gumann, P.
Nr. 12
Salden, D.
(1980)
Nr. 13
Seeger, H.
(1980)
Nr. 14
Schmidt, H.H.
(1981)
Nr. 15
Smoltczyk, U. /
Schweikert, O.
(1981)
Nr. 16
Malcharek, K. /
Smoltczyk, U.
(1981)
Nr. 17
Gruhle, H.D.
(1981)
Nr. 18
Kobler, W.
(1982)
Nr. 19
Lutz, W.
(1983)
Nr. 20
Smoltczyk, U.
(1983)
190
Boden
Stutzw
ande bei nachgiebigem Baugrund.
e 12,78
Vorstudie uber
bauliche Alternativen fur
Durchgangsstraen in Siedlungen.
e 6,14
die BerechVergleich nationaler Richtlinien fur
nung von Fundamenten.
e 7,67
Das Verhalten des Baugrundes unter Einwirkung vertikal gezogener Ankerplatten als
raumliches Problem des Erdwiderstandes.
Vergriffen
Untersuchungen uber
Boschungsund Grundbruch bei begrenzten Lastflachen.
e 12,78
Tragfahigkeit des geschlitzten Baugrunds
neben Linienlasten.
e 12,78
Studienunterlagen
Bodenmechanik
und
Grundbau; uberarbeitete
Ausgabe 1983.
e 20,45
Geotechnik
Mitteilungen des Instituts fur
Nr. 21
Schweikert, O.
(1984)
Nr. 22
Vogt, N.
(1984)
Nr. 23
Buchmaier, R.
(1985)
Nr. 24
Schad, H.
(1985)
Smoltczyk, U. /
Schad, H. /
Zoller, P.
Nr. 25
Gumann, P.
(1986)
Nr. 26
Steinmann, B.
(1985)
Nr. 27
Lee, S.D.
(1987)
Nr. 28
Kolb, H.
(1988)
Nr. 29
Ochmann, H.
(1988)
Nr. 30
Breinlinger, F.
(1989)
Nr. 31
Smoltczyk, U. /
Breinlinger, F. /
Schad, H. /
Wittlinger, M.
(1989)
Moglichkeiten
der Boschungssicherung
bei
kleinen Baugruben.
e 12,78
191
Geotechnik
Mitteilungen des Instituts fur
Nr. 32
Gumann, P. /
Schanz, T. /
Smoltczyk, U. /
Willand, E.
(1990)
Nr. 33
Gruhle, H.D.
(1990)
Nr. 34
Henne, J.
(1995)
Nr. 35
Wittlinger, M.
(1994)
Nr. 36
Schad, H.
(1992)
Nr. 37
Belz, I.
(1992)
Nr. 38
Ma, J.
(1994)
Nr. 39
Smoltczyk, U.
(1994)
Nr. 40
Rilling, B.
(1994)
Nr. 41
Vermeer, P.A.
(1996)
Nr. 42
Vermeer, P.A.
(1997)
Nr. 43
Brinkmann, C.
(1998)
192
Vergriffen
Stutzscheiben
stabilisierten Boschungen.
e 15,34
Sonderheft: 25 Jahre Lehre und Forschung in
der Geotechnik.
e 15,34
Untersuchungen zur Grenztragfahigkeit bindi
ger Schuttstoffe
am Beispiel von Lolehm.
e 17,90
Deponiebau und Geotechnik.
e 17,90
Baugruben in Locker- und Festgestein.
e 17,90
Untersuchungen zum Verhalten von Dich
tungsuberg
angen im Staudammbau.
e 17,90
Geotechnik
Mitteilungen des Instituts fur
Nr. 44
Fiechter-Scharr,
I.
(1998)
Nr. 45
Schanz, T.
(1998)
Nr. 46
Akinrogunde,
A.E.
(1999)
Nr. 47
Vogt-Breyer, C.
(1999)
Nr. 48
Vermeer, P.A.
(1999)
Nr. 49
Marcher, T.
(2002)
Nr. 50
Marcher, T.
(2003)
Nr. 51
Ruse, N.M.
(2004)
Nr. 52
Beutinger, P.H.
(2005)
Nr. 53
Wehnert, M.
(2006)
Nr. 54
Moller,
S.C.
(2006)
Nr. 55
Benz, T.
(2007)
193