Dilatancy For Cohesionless Soils DAFALIAS PDF
Dilatancy For Cohesionless Soils DAFALIAS PDF
Dilatancy For Cohesionless Soils DAFALIAS PDF
Gc!U/(!Chllicjl{(!
INTRODUCTION
= d(17 , C)
(I)
449
450
Ll AND DAFALlAS
parameters for a single sand and thus does not have a good
control over changes in the material state during loading.
Attempts have been made in recent years to tackle this issue
from the perspective of dilatancy. With the concept of critical
state as basis, Been & Jefferies (1985) introduced a scalar
quantity Vi called the state parameter, which measures the
difference between the current and critical void ratios at the
same p' . Kabilamany & Ishihara (1990) provided experimental
evidence showing that d + 'q increases as shear deformation
increases. Manzari & Dafalias (1997) presented a sand model in
which a linear dependence of the phase transformation, or
dilatancy stress ratio (the stress ratio at which the response
changes from contractive to dilative) on Vi was introduced. Li
(1997) investigated the response of sand at the ultimate stress
ratio and explicitly pointed out that the dilatancy d is not
related only to the stress ratio but is also a function of plastic
volumetric strain. More recently, Wan & Guo (1998) proposed a
model with its dilatancy modified from Rowe's stress-dilatancy
equation . The modified dilatancy equation includes the density
dependence with the critical void ratio as a reference.
Cubrinovski & Ishihara (1998) also showed a dilatancy relationship that depends on the material state represented by cumulative plastic shear strain. Li et 01 . (1999) introduced a statedependent dilatancy .into an existing hypoplasticity sand model
(Wang et al., 1990), resulting in a successful simulation of the
responses of Toyoura sand to both drained and undrained
triaxial loading over a wide range of densities and pressures.
In the present paper a number of issues on this subject are
discussed, starting from some microscopic analytical considerations and ending with the presentation of a simple macroscopic
constitutive framework and modelling, the simulative capability
of which is shown by comparison with experimental data of the
response of a sand under various initial state and loading
conditions.
0;
tan(fJ
+ (3)
0;
= -02 = tan(q')". . + (3) tanf3
(3)
= 8 sin f3 cos f3 _ 1
(4)
:rr
0;
2sinf3
--,
= tan(fJ + (3) I + 2 cos f3
02
(5)
and
e=
(6)
3:rr
(2)
(7)
02(\
+ dV / VEt}
tanf3
i. \
1\ =60
L~. --II~
L, --II~
L
I
L .
L,
451
I.. .
L?
L2
.Lt.
L2
Lz
L1
L1
10 r - -- - - -- - -- - - - -- -- - - - . -- - - - ,
O~;
............
+.. 0'6
~, 0-4
't;
--
02
F'ael--in9 A
- - Packing B
<!
0:
'cl
(1 ,261
Equation Ci!
2
dVi V~;
452
Ll AND DAFALIAS
1GOO
'.000
r-------------------,,-----,
e = 083:3
Dr = 3'7 '9%
1400
1200
1500
0:;
0:;
2;: 1000
0-
~
0
ro;;
800
2;:
')0/0//0/
'0'
ui
<fl
Q)
~ 1000
0
ro;;
Q)
<ll
600
...........
400
500
Dense state, d < 0
200
o .""
o
200
400
600
800
1000
1200
1500
Effecti ve mean normal stre ss. p' : kPa
(a )
(b )
2000
Fig. 4. Variation in dilatancy with material state (data from Verdugo & Ishihara (1996 . Undrained response of a sand with (a) different
densities and (b) the same density but under different confining pressures
4000 , - - - - - - - - - - - - - - - - - - - - - - --,
~ ~IOOO
Va ri eci pll ase
tidflsior lTIatiC! n
Ul
~ "O()(1
tn
L.
C>
axis. Eventually, when the sand becomes too loose, the phase
transformation phenomenon totally disappears. The assumption
that d is uniquely related to 1J again contradicts the observation .
Furthermore, undrained tests on dense sand often show that
the q- pI stress path eventually converges with a line of more
or less constant 1J = M towards an ultimate state (Figs 5 and
6). At the ultimate state, both the stresses and the plastic
volumetric strain stop changi ng, as shown in Fig. 6 (see also
Figs 10- 12). Since the stresses do not change, neither does the
elastic volumetric strain . By definition , this is a critical state
where dpl = dq = dE" = 0 while dEq =f. O. Along the approximately '7 = M path and before the critical state is reached, the
fact that Ee tends towards a constant value implies that the
dilatancy d tends towards a zero va lue . If d were a unique
function of 'I , however, along thi s path d wou ld be esse nti ally a
constant, which means that, as shear proceeds, pI would increase contin uou sly as a result of the undrained constraint of
zero total vo lum etric change, and the criti cal state wou ld never
be reached. One m ay argue that along this path the stress ratio
'I only approaches M asymptotica ll y, and correspondingly the
evo lu tion of d from a no n-zero value towards zero would be a
result of the tiny deviations in '7 frolll M . However, this
argum en t on ly facilitates a mathematical description that barely
make s equation ( I) not violate the concept of a critical state.
Cons id ering the un certaint ies involved in soil testing, one would
not be able to physically identify and quantify sllch tiny
deviations , if any, in a mean ingful manner. On the other hand,
wit h the hypothesis that d and '] are not uniquely re lated. such
an argume nt becomes unnecessary.
The above observations lead to the concl usion that a sa nd
model with its dilatancy fo ll owing equatio n (I) works well on ly
when the change in the material internal state is minor.
(ii
;;
2;
1000
GOO
1000
1500
200\)
2500
Efieclive mean narrn,]1 sir ' 5S, p ' kP,]
3000
(8)
453
0)
dL P
Ultimate/critica l slate
.. .
(d = 0)
....
... . .. .. .. .......... - ... ... .... .. .. .... .. .. .. ...... .. ...... .. - .. _. _._. -... -. -- .._. _. _. - .:.;
-
"-;
Ultimate/critical
states
(d;i 0)
(b)
(8)
Fig. 6. Illustration of the dilative shear on the failure surface: (a) stress path; (b) stress-strain response
M, e
= ee, Q, C) = 0
(9)
ClllTen l sl<l te ;0
Q.
Secondly, it is possible for sand to hav e a so-called ' phase
transformation state' at which d = 0 but 17 i- M and e i- ee, as
discussed earlier. Analytically thi s means that the equation
Il' > 0
(contractive)
u
>
( 10)
can be lI sed to specify th e combination of 'I , e , and Q that
defin es a phase transformation state . Co nverse ly, one may use
an ([ priori experimental knowledge of phase transformation
states to specify via equation (10), together with equation (9),
app ropri ate forms of the dilatancy function in equation (8).
(pip ):
Ll AND DAFALIAS
454
= d(rJ , l/J, C)
(12)
_~
3G
I
KMd
_ _
11 =114-
( 13)
I)
( 14)
r = q -lIP' = 0
(16)
(Of ,
I
of )
L =-,dp + -::----dq
Kp up
()q
dq - '7 dp'
Kp
P'd17
Kp
(17)
d f~
{ dE~
} _
L{
1 } _ { p'
d
Therefore. for L
df q
_
c
- d Ep
( 18 )
> 0:
p _
+ d Eq
dl}1Kp
d p ' dlJ l Kp
dq
p' elI7 _
3G +
( 19)
dp'
d
Kp
d e\'= dE \, + d E~ =- +d d f ~=-dq+
( I
d17)
~--
1\
Kp
dp '
(20)
Equ ation s (19) and (20) establish the relationship between the
stress and strain increments. They can be inverted by a straightforwa rd algebraic manipulation , and expressed in a matri x form
as:
[(
30G
0)
h(L)
( 9C
Kp + 3C - K1}d 3KCd
(21 )
G = Go
(2.97 - e)2
I+ e
V P Pa
G _2_(1_+_v_)
3(1 - 2v)
(22)
dEq
(3~ +
q
- =
(a)
For the hardening and softening responses before and after the
peaks , Kp is positive and negative as 17 < Mb and 17 > MO,
respectively. It can be seen that conditions (a) and (b) are
dq =
C~ + hG(q~ _ q) dq
(I I I
- - - - - I n ( 1 - r)
3G hG hGr .
Eq
(24)
;J
(26)
(23)
455
)-1
(27)
where r = q I qr. Equation (27) can be converted into a normalized modulus reduction curve (secant shear modulus normalized
to its maximum value G max versus shear strain) with h as a
parameter. Fig. 8 shows a family of such curves together with
the curve based on the hyperbolic stress-strain relationship
ql Eq = 3G ma x /( I + EIEqr) , where G max and Eqr are two material
constants. It can be seen that the Kp function used here allows
more flexibility in calibrating the shear stress-strain response
than does using the hyperbolic stress-strain response.
12 ,------,---
--,---
--,---,----,-----,.----,------,
11 :: 6
(5
~ (Hi
-+-- h :::3
~ 0 4 - --- h::: 1 2
::;
~ 02 - ... 11 =06
.......
h :~ 0 3
0 03
0 1
03
10
"
LI AND DAFALIAS
456
Table I. Model parameters calibntted for Toyoura sand
Critical state
parameters
Elastic
parameters
Gn = 125
v = 005
Hardening
parameters
Dilatancy
parameters
125
do
er = 0934
Xe = 00 19
III
= 088
hi
= 3.5
3 15
h2 = 305
n
= 11
,; = 07
test data for the critical stress ratio and the critical state line in
the e- p' plane.
The parameter 111 can be determined by equation (15) at a
phase transformation state, at which d = O. Hence,
I
Md
m=-In-
(28)
ljJd
where l/Jd and Md are the values of 1/) and 17 at the phase
transformation state, measured from drained or undrained test
results .
The parameter 11 can be determined by equation (24) at a
drained peak stress state, at which Kp = O. Hence,
I
n =-In l/) b
Mb
(29)
where l/Jb and Mb are the values of l/J and 17 at the drained peak
stress state, measured from test results.
Next, consider the drained tests. Ignoring the smal1 elastic
deformations,
dEy
~ dE~
dE~
dEq
d= do (e
1111/' _
!l..)
M
(30)
hGo{
+ e)( 1- m71
dp'
Kd
2Go( 1 + v)
dq
dEq
I I)
(-+
3G
Kp -
-I
1] Kd
dq
~- -
dE~ -
K -17Kd
(33)
(31 )
in which the parameter (J is either equal to 1/ 3 (for conventional tests) or to zero (for constant p' tests). As al l the material
constants in the brackets have been predetermined, the COl1lbined parameter hG o ca n be calibrated independently based on
the experimental q - ELI curves. It may be found during calibration that the quantity hGo varies with density. Fitting these
values of hGn into equation (25) yie lds the constants hi and h2
(aftcr Go has been determincd).
Now let us turn to undrained (constant volume) tests . For
dE" = 0, equation (20) yields :
value of 2Go( I + 1')/ 3(1 - 2v) , and could be negative. A micromechanics study (Chang & Misra, 1990) has shown that the
Poisson 's ratio of an assembly of particles is predominantly
control1ed by the ratio of the shear stiffness to the normal
stiffness at particle contacts. The value of the Poisson's ratio of
the assembly could be much lower than that for the particle
material itself. If this stiffness ratio is high, v could be negative.
Even though a negative v is theoretically justifiable and affects
nothing but the volumetric strain at extremely low strain ievel
(pure elastic range), if it is encountered and disliked, as an
alternative one may pick a v value first and then calculate Go.
As this alternative approach does not guarantee the accuracy of
the elastic shear response, it should be used only when the
shear stiffness in the elastic range is unimportant or when
accurate values of G are unavailable. Once Go has been
determined, h, and therefore hI and h2, can finally be found
from equations (25) and (31).
Last, but not least, one can obtain the undrained deviatoric
stress-strain response by substituting dp' = K dE~ = - K dE~
= -Kd dE~ into equation (32) and accounting for dq = 3G dE~.
The relationship is as follows:
(32)
As al1 the material co nstants in the brackets have been predetermined, the combined parameter 2Go( I + v)/3( I - 2v) becomes the only means at this stage of adjusting the undrained
p' - q responses of the model. By matching these responses with
their experimental counterparts, the value of 2G o( I + v) /
3( I - 2v) can be determined.
Finally, one needs to separate the parameter Go from hand
v . If shear stiffness at small strains is important , Go should be
determined by independent small strain tests, such as resonant
column tests or bendcr element tests, through fitting the test
data into equation (22). However, if Go is high, the value of
Poisson 's ratio ]I is reduced based on the already calibrated
457
2500 r---------------------------------------------.
OY5 r----------------------------------------------,
0-90
2000
ro
2z
'C7
u;
~ 085
(f,)
,9
Ui
ro
'D
g 080
'~
Expllrirn ental
075
500
1000
1500
2000
2500
3000
25
30
2500
co 2000
2z
'0-
~ 1500
e:
5000r---------------------------------------------~
co 4000
---::::::::-;:::-.:.-...
Ui
- - Model simulation
ro
':;;
Q)
2z
'0~
3000
~
E
2000
co
10
Q)
20
15
Axial strain , 1'1: a::n
(b)
:;;
1000
1000
500
2000
2500
3000
1200
5000 r---------------------------------------------~
- - Model simulation
co
ro
4000
2z
~
.u;.
E
2000
cr
2z
.. ~ .....--.. - .--.
u; JOO()
./
.--
~
' " Pu'
..Plj' : i
1000
:.:1- SOO
e:
00 kPa
"th
)00
(5
1000 kPa
n,;
:;
.~
400
Q;
200
10
15
20
2[;
2500
30
:lOOO
(8)
(tJ)
12ClO
1000
ro
c:
C'"
2000 kPa
800
UJ
if;
7T:
.9
\"1:.;;:
CJ
[)
10
15
P.Xi 81 ,;trail1 , I'
20
30
. ~:i,
(Il)
LI AND DAFALIAS
458
3'.in
CONCLUSIONS
- - IvlocJ el sirnulil lion
300
Tes t results
ru
~ 250
.&
~ 200
150
g
in
0:;
100
Q)
50
08
095
09
085
Void ratio.
10
(a )
350
300
m
22 250
0-
m
sQ)
0
10
15
20
25
30
A xi al strain , {' 1: %,
(b)
1600
Model sim ulation
ACKNOWLEDGEMENTS
Test results
~ 1 200
.;,L
"c
t~
BOO
1]
>
(!)
-100
(34)
where L; is the load per rod in direction i, 1)" is the angle of friction
between the surfaces of the rods in contact, and (J is the deviation of the
tangent at the contact points from the direction I. The size of a basic cell
in the packing is characteri zed by
Of)
09
() 85
Void
095
"li'o, e
la)
[35a)
I, = 4rsin/j
1600
and
(35b)
I" = 4reos/i
Therefore, the ratio of the major principal st ress to the minor principal
stress is
~ 1200
0-
~:
rr,
'>J)
800
(36)
m
;;
E"
-=
r.J
-dii/ v-f,
E,
400
E,
=-
( 1+
dV)
VEl
i"l,
')
=~ = -tan - II
(37)
I,I?
Therefore, the ratio of the work done per unit volume by the major
principal stress to the work done on the minor principal stress is
0
10
15
A xial strain
(t))
20
I,
25
30
0'
: \;
aiEl
0;
OlE?
tan(<j.l,,+13)
tan 13
(38)
test
I:'
= 1,1" - 2m'?'
2:71'~
::r
(39)
459
iGOO
Moclel sim ula ti on
300
Tes t re sults
Test results
~ 1200
t:i-
Q'
~ 250
~ 200
~
0
rn',>
vi
<1)
(jj
150
800
rn
'>
Q)
100
Q)
400
50
0
08
0,88
084
0,96
092
Void ralio,
10
084
08
088
092
096
Void ratio, e
1600
350
300
~ 250
-.-..-.
t:i-
t:i-
~ 200
~
(jj
e!
(jj
rn
',>
Q]
~ 1200
ru
150
800
e Q =0'810
rn
',>
100
Q)
400
50
0
0
20
10
30
10
Fig. 15. Simulations of drained tests with h varying with eo instead of e: (a) pu
Packing B (Fig. 2)
At each contact point, sliding takes place when the ratio of thc loads
LI
'
(3'
Ld2 = tan(if>J1 + , )
(40)
I I = 4rsin/~
and
(4Ib)
Thereforc, the ratio of the major principal stress to the minor principal
stress is
0;
2L I/ I
J 2sin/3
L212 = tan(r/J p +/) 1+ 2cos(J
(42)
h
(
di; )
~= - I+~
2sin2/J
i2/1
=- - (I
= il/2
+ 2 cos (3)cos/J
(43)
Thercfore, the ratio of the work done per unit voluille by the major
principal strcss to the work done on the minor principal stress is equal to
tJ; EI
tJ2E2
0;
02(1 +dil / vEI)
tan(if>f l + /-1)
tan/3
REFERENCES
Becn, K. & Jefferies, M.
Geofechnique 35, No , 2,
Bolton, M. D. (\ 986). The
niqlle 36. No. I , 65 - 78.
Chang. C. S. & Misra. A.
30
E 1: ~/O
(b)
(a)
01 =
20
Axial strain.
(44)
1=
= 100 kPa,
It
hI - "zeo; (b) Po
= 500 kPa,
hI - "zen
460
LI AN D DAFALIAS
idealized 'wet' clay. Pmc-. EII/: COllI Sail Mech. FOllnd Ellgllg.
Wiesbaden 1,47 - 54.
Roscoe, K. H., Schofield, A. N. & Wroth , C. P. (1 958). On the yielding
of soi ls. Geatechnique 8, No. I , 22 - 53.
Rowe, P. W. (1962). The stress - dilatancy relation for stnti c equilibrium
of an assembly of particles in contact. Pmc. R. Soc. , LOlldon, SeJ: A
269, 500-527.
Ri chart , F. E. Jr., Hall , J. R . & Woods. R. D. (1970). Vibrations of" soils
and foundations. Englewood Cliffs, NJ: Prentice-Hall.
Taylor, D. W. ( 1948). Fllndam el1lals of soil mechanics. New York:
Wiley.
Vardoulakis, I. & Sulem. J. (1995). BijiJrcatian all a lysis in geomeehanics. Glasgow : Blackie Academic & Professional.