A Continuous Hyperplasticity Model For Sands Under Cyclic Loading
A Continuous Hyperplasticity Model For Sands Under Cyclic Loading
A Continuous Hyperplasticity Model For Sands Under Cyclic Loading
+ +
+ =
N
n
n
q
n n
q
n
p
H q p
N
G
q
K
p
g
1
2
) ( ) ( ) ( ) (
2 2
2
3
'
1
6 2
'
(5)
( )
= =
= =
N
n
n
q
N
n
n
p
N
n
M
N
d
N
d
1
) (
1
'
1 1
& (6)
where
) (n
H is the hardening modulus related to the
th
n mechanism and M is the value that the stress
ratio ' / p q attains at critical state conditions. The
normalisation term N / 1 in (5) and (6) makes the
formulation independent on the number of surfaces.
In the previous model (Houlsby, 1992) the energy
function used was the Helmholtz free energy func-
tion, while the dissipation function was formulated
in terms of strains. However, the two energy func-
tions are linked by the Legendre transformation
q p
q p g f + + = and either g and f can be used
for describing the behaviour of material.
The definition of appropriate constrains enables
the introduction of dilation as well as anisotropy into
the model (Houlsby, 1992):
( )
0
) ( ) ( ) (
= + + = a C
n
q d
n
q c
n
p
n
d
& & & (7)
0
1
1
) (
1
) (
=
=
= =
N
n
n
q
N
n
n
q a
a A
N
a C & & & (8)
The first constraint specifies that dilation is made up
from isotropic (compressive) and anisotropic (dila-
tive) parts given by functions
c
and
d
respec-
tively. The second constraint specifies the evolution
of the anisotropy parameter a, which varies between
1 + and 1 for positive and negative shearing. The
rate of evolution of anisotropy is determined by the
constant A. With dissipation specified, it is possible
to obtain the yield function
) (n
y associated with each
set of plastic strains through the degenerate special
case of the Legendre transformation of
) (n
d , which
is homogeneous of degree 1 in the rates:
0
) ( ) ( ) ( ) ( ) ( ) ( ) (
= + =
n n
q
n
q
n
p
n
p
n n
d y & & (9)
where by definition:
) (
1
) (
) (
) (
) (
) (
) (
1
) (
) (
) (
) (
) (
n
q
a
a
N
i
n
q
i
d i
d
n
q
n
q
n
p
a
a
N
i
n
p
i
d i
d
n
p
n
p
C
C
d
C
C
d
=
=
& & &
& & &
(10)
ln v
ln p'
ln p'
ref
ln B
ln
ln D
ln v
compression dilation
ln B
ln
ln D
M
n
N
(1- a sgn(q ))
. (n)
. (n)
max
n
N
a sgn(q )
where
) (n
d
and
a
are Lagrangean multipliers.
From the constraints (7) and (8) it follows that
( )
p
n
d
= and 0 =
a
. The generalised stresses are:
) ( ) (
) (
) (
) (
) (
3
n
q
n
n
q
n
q
n
p
n
p
H q
g
p
g
=
=
=
=
(11)
and Zieglers orthogonality condition
) ( ) ( n n
=
leads to the yield function in terms of stress:
( ) ' sgn 3
) ( ) ( ) (
p a
N
n
M H q
n
q d c
n
q
n
+ + = & (12)
The yield surfaces exhibit kinematic hardening,
given by the term
n
q
n
H
) (
3 , where the expression for
the variation of the hardening modulus is:
b
n
N
n
h H
= 1
) (
(13)
with h and b being parameters of the model. To in-
troduce the difference between compression and ex-
tension, the critical stress ratio M is given by:
( ) ( ) ( ) ( )
) (
sgn
2
1
n
q e c e c
M M M M M + + = & (14)
where
c ec e
M r M = , and
ec
r is the ratio between the
critical stress ratios in extension
e
M and compres-
sion
c
M . As in the previous model (Houlsby, 1992),
the values of
c
and
d
depend on the state of the
soil, defined by the distance between the current
specific volume and the critical state line, which is
assumed to be linear in a bi-logarithmic plot:
= +
ln ln
ln ln
ln ln
ln ln
max
D
v
N
n
D B
D v
N
n
M
d c
c
(15)
where
max
is the maximum rate of dilation while
&
&
&
&
(16)
where the first term refers to compression and the
second to dilation. It is worth noting that for =
v
the second term is always zero while the first one is
null only when the term ) sgn( 1
) (n
q
a vanishes.
Figure 1 shows the graphical interpretation of
contractive and dilative terms in (16).
Figure 1. Density constants and graphical interpretation of
compression and dilation rules (equation 16).
3 EXAMPLE ANALYSES
We illustrate the model by example analyses of ide-
alised tests. Although we do not compare these here
with specific data sets, the patterns of behaviour cor-
respond to those that are well-established empiri-
cally. The example calculations are carried out using
the parameter values given in Table 1.
Figure 2 shows a set of drained constant mean
pressure tests on sands with different initial densi-
ties. As the index of density
v increases (looser
samples) the strength reduces and the samples
change from being strongly dilative to contractive.
Although not apparent in Figure 2, the denser sands
show a mild peak in the stress-strain response.
Figure 3 shows the results for drained cycling
over a constant range of strain. The upper plots are
for a loose sample, which exhibits an accumulation
of compressive strain. The resulting densification
causes a slight increase in stiffness of the response to
the cycles. The lower plots show the equivalent for a
dense sand. This time the sand dilates during the
K Bulk modulus 20 MPa
G Shear modulus 10 MPa
H Hardening modulus constant 2 MPa
B Hardening modulus constant 2
A Rate of anisotropy development 100
B Specific volume at loosest state at
ref
p p = 1.9
Specific volume at critical state at
ref
p p = 1.8
D Specific volume at densest state at
ref
p p = 1.6
M
c
Critical stress ratio in compression 1.2
max
Maximum rate of dilation 0.4
r
ec
Ratio between critical stress ratios in extension
and compression
0.8
Table 1: Example parameters for model
0 2 4 6 8 10
deviatoric strain
q
(%)
0
40
80
120
160
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
0 2 4 6 8 10
deviatoric strain
q
(%)
-2
-1
0
1
v
o
l
u
m
e
t
r
i
c
s
t
r
a
i
n
p
(
%
)
v
Figure 2. Response of the model in drained compression tests
on sands with different initial densities.
-2 0 2
deviatoric strain
q
(%)
-120
-80
-40
0
40
80
120
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
-2 0 2
deviatoric strain
q
(%)
0
2
4
6
8
v
o
l
u
m
e
t
r
i
c
s
t
r
a
i
n
p
(
%
)
-2 0 2
deviatoric strain
q
(%)
-150
-100
-50
0
50
100
150
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
-2 0 2
deviatoric strain
q
(%)
-6
-4
-2
0
2
v
o
l
u
m
e
t
r
i
c
s
t
r
a
i
n
p
(
%
)
p' = 100 kPa
v
= 1.85
p' = 100 kPa
v
= 1.65
p' = 100 kPa
v
= 1.85
p' = 100 kPa
v
= 1.65
Figure 3. Effect of relative density on drained cyclic tests.
cycling, resulting in a slight reduction of stress in the
cycles as the material loosens.
Stress-strain curves and effective stress paths for
undrained monotonic tests on samples of different
densities are shown in Figure 4. The loose sands
show a reduction in effective stress, whilst dense
sands show a strong increase in effective stress.
Sands of medium density show a slight reduction of
mean effective stress before the increase. This
pattern of behaviour is well known.
Undrained cycling over a constant stress ampli-
tude is shown in Figure 5 for two densities of sand.
A loose sand (upper plots) shows an initial reduction
in effective stress, after which the q p , plot settles
into a characteristic butterfly shape, and the strain
amplitude also becomes constant. There is a large
amount of hysteresis in each cycle. The dense sand
shows a similar pattern, but the butterfly plot is
narrower, and the response both stiffer and with less
hysteresis.
Figure 6 shows undrained cycling on a dense
sand at a higher stress range. Although a stable
butterfly pattern is developed, note that this time
the dilation during each cycle means that there is a
net increase rather than a decrease of mean effective
stress.
Finally, figure 7 shows the effect of the strain
amplitude in on the stress-paths of a loose sample in
constant strain amplitude undrained cycling. For the
larger amplitude ( % 2 =
q
) stabilization of stress
is achieved while for the other ( % 1 =
q
) cyclic
liquefaction is obtained after just two cycles.
4 CONCLUSIONS
A model for the cyclic behaviour of sand under
triaxial conditions has been presented. The model
successfully describes typical trends of behaviour
for undrained and drained cycling, including typical
variation of volumetric behaviour for sands of
different densities.
REFERENCES
Houlsby, G.T. (1992) "Interpretation of Dilation as a
Kinematic Constraint", Proceedings of the Workshop on
Modern Approaches to Plasticity, Horton, Greece, J une 12-
16, ISBN 0-444-89970-7, pp 19-38
Puzrin, A.M. and Houlsby, G.T. (2001a) "A Thermomechani-
cal Framework for Rate-Independent Dissipative Materials
with Internal Functions", Int. J our. of Plasticity, Vol. 17, pp
1147-1165
Puzrin, A.M. and Houlsby, G.T. (2001b) "Fundamentals of
Kinematic Hardening Hyperplasticity", Int. J our. of Solids
and Structures, Vol. 38, No. 21, May, pp 3771-3794
0 2 4 6 8
deviatoric strain
q
(%)
0
100
200
300
m
e
a
n
s
t
r
e
s
s
p
'
(
k
P
a
)
0 100 200 300
mean stress p' (kPa)
0
100
200
300
400
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
0 2 4 6 8
deviatoric strain
q
(%)
0
100
200
300
400
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
v
Figure 4. Response of the model in undrained compression
tests on sands with different initial densities.
-4 -3 -2 -1 0 1 2
deviatoric strain
q
(%)
-80
-40
0
40
80
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
0 20 40 60 80 100
mean stress p' (kPa)
-80
-40
0
40
80
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
-4 -2 0 2
deviatoric strain
q
(%)
-80
-40
0
40
80
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
0 20 40 60 80 100
mean stress p' (kPa)
-80
-40
0
40
80
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
-2 -1 0 1 2
deviatoric strain
q
(%)
-120
-80
-40
0
40
80
120
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
80 100 120 140 160
mean stress p' (kPa)
-120
-80
-40
0
40
80
120
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
0 20 40 60 80 100
mean stress p' (kPa)
-80
-40
0
40
80
120
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
0 20 40 60 80 100
mean stress p' (kPa)
-80
-40
0
40
80
120
d
e
v
i
a
t
o
r
i
c
s
t
r
e
s
s
q
(
k
P
a
)
q = 2%
q = 1%
Figure 5. Simulation of stress controlled undrained cyclic tests
for loose and dense sand.
Figure 6. Large stress controlled undrained cyclic tests for
dense sand showing increase of mean effective stress.
Figure 7. Effect of cyclic strain amplitude on the behaviour of
loose samples