# An Elastic-Viscoplastic Model for Time-Dependent Behavior of Soft Clay and Its Application on Rheological Consolidation

Authors: Jinzhu Li; Wenjun Wang; Yijun Zhu; Haofeng Xu; Xinyu Xie
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

To describe the time-dependent behavior of soft clay, this paper extended one-dimensional Nishihara model to three-dimensional stress state based on the framework of Perzyna’s overstress theory and modified cam-clay model. The yield criterion of modified cam-clay model was used to describe the plastic properties of soft clay, and the overstress theory was used to describe the strain rate effect. Triaxial rheological tests were carried out on Ningbo soft clay and the rheological characteristics were studied. Based on laboratory results, the parameters of proposed model were determined by curve fitting, which show that this model is suitable for the rheological characteristics of Ningbo soft clay. The analysis of parameters shows that, the value of parameters changes slightly with different deviatoric stress when the confining pressure was constant, but changes notably with the increase of confining pressure. A user material subroutine of the proposed constitutive mode was coded on the platform of the FEM software ABAQUS and verified by triaxial compression of soil column. A plain strain problem was computed to analyze the rheological consolidation properties of soft clay, in which the rheological effect and the finite strain effect were considered.


---

## Body

## 1. Introduction

According to the classic consolidation theory presented by Terzaghi, when soft clay is consolidated under external loads, the excessive pore pressure dissipates and the effective stress increases, while the deformation grows until the excessive pore pressure dissipates completely. This theory has laid the foundation of soft soil mechanics and even became the main theoretical basis of the computation of soft soil foundation. However, in the former works, Buisman first discovered the long term deformation of soils after the excessive pore pressure was dissipated in experiments, which is also reported in the works of Zeevaert [1], Leonards and Ramiah [2], and also in lots of engineering practices. This long term deformation which cannot be explained by classic theories is usually called rheology, which contains creep, stress relaxation, long-term strength, elastic aftereffect, hysteresis effect, and so on. Because of strict requirement of post-construction settlement, researchers are paying increasing attention on the rheology of soft clay. To propose a suitable constitutive model to describe the rheology of soft clay and to determine the model parameters accurately by laboratory tests, these are the key points of the topic. Since 1960s, a number of constitutive models were proposed for time-dependent behavior of soft clay. Commonly used rheological constitutive models can be divided into three categories: empirical models, component models, and microscopic models, in which the first two categories were most widely used. Empirical models such as Singh-Mitchell model [3], Mesri et al. model [4, 5] are mostly one dimensional models, which were often used to calculate the rheological deformation in one dimensional problem but were limited in complex stress state. Component models are intuitive in physical meaning, which can describe numerous rheological models by different combinations of components and can be easily extended to multidimensional problems to apply in numeral computation for complex stress state, such as the work of Adachi and Oka [6], Fodil et al. [7], Yin et al. [8, 9], Hinchberger and Rowe [10], Hinchberger and Qu [11]. Based on the framework of Perzyna’s overstress theory [12] and modified cam-clay model, this paper proposed a three-dimensional rheological model. The model was verified by laboratory triaxial rheological tests of Ningbo soft clay, and the model parameters were determined by curve fitting. On the platform of ABAQUS, a user material subroutine of the model was coded and verified.


## 2. Model for Soft Clay

As one of commonly used component models, Nishihara model was connected by Hooke body, Kelvin body, and Bingham body, which can describe elastic deformation, creep, elastic aftereffect, and viscous flow comprehensively, such as Figure1. When a constant external pressure σ was given, the stress and strain of the model can be written as follows:(1a)
σ
=
σ
e
=
σ
v
e
=
σ
v
p
,
(1b)
ε
=
ε
e
+
ε
v
e
+
ε
v
p
.Figure 1
Sketch of Nishihara model.In the above equation,σ
e, σ
v
e, σ
v
p denote the stress of Hooke body, Kelvin body, and Bingham body, respectively; ε, ε
e, ε
v
e, ε
v
p denote the total strain, the strain of Hooke body, Kelvin body, and Bingham body, respectively.Substituting the stress-strain relationship of Kelvin body and Bingham body into (1a) and (1b), the constitutive equation of one-dimensional Nishihara model can be obtained as follows:
(2)
ε
=
σ
E
0
+
σ
E
1
(
1
-
e
-
(
E
1
/
η
1
)
t
)
+
〈
σ
-
σ
s
〉
η
2
t
,
where E
0 denotes the elastic modulus of Hooke body, E
1 and η
1 denote the elastic modulus and viscosity coefficient of Kelvin body, η
2 and σ
s denote the viscosity coefficient and yield stress of Bingham body, 〈
σ
-
σ
s
〉 denotes a switch function, which can be defined as follows:
(3)
〈
f
〉
=
{
0
(
f
≤
0
)
f
(
f
>
0
)
.Equation (2) contains the first two stages of creep, which are often called primary creep stage and steady creep stage. In fact, for complex stress state, (2) cannot be used to calculate the deformation of soil, it is necessary to extend (2) to three dimensional, so some hypotheses can be presented as follows:(1)
the soil material is isotropic;
(2)
the volume deformation is just caused by spherical stress, and it is unrelated to shearing stress;
(3)
the Poisson’s ratio of soil is constant and does not change with the stress or time.Based on the above hypotheses, the stress-strain relationship of Hooke body is given as(4)
ε
i
j
e
=
I
1
9
K
δ
i
j
+
s
i
j
2
G
0
,
where I
1 is the first stress invariant and K and G
0 denote the bulk modulus and shear modulus of Hooke body.For Kelvin body, the stress-strain relationship is given as(5)
ε
˙
i
j
v
e
=
s
i
j
2
H
1
+
G
1
H
1
ε
i
j
v
e
,
where G
1 and H
1 denote the shear modulus and 3D viscosity coefficient of Kelvin body, in which H
1
=
2
(
1
+
μ
)
/
η
1, μ is the Poisson’s ratio of soil.For Bingham body, the stress-strain relationship is given as(6)
ε
˙
i
j
v
p
=
1
H
2
〈
ϕ
(
f
)
〉
∂
Q
∂
σ
i
j
.In (6), H
2 denotes 3D viscosity coefficient of Bingham body, and 〈
ϕ
(
f
)
〉 is a switch function to judge whether plastic yield is occurring and whether the amplitude of the plastic yield occurred, Q is the plastic potential function, when the associated flow rule is adopted, Q
=
F, where F is the yield function. ϕ
(
f
) can be defined as follows:
(7)
ϕ
(
f
)
=
(
f
-
f
0
f
0
)
N
,
where f is the current yield function, f
0 is the initial yield function, and N is a constant derived from experiments; according to the work of Wei [13], the value of N can be assigned to 1.0 approximately for soft clay.According to the work of Adachi and Oka [6] and Yin et al. [8, 9], we define a static yield criterion f
0, which represents a reference yield surface for the material, such as Figure 2. Its initial shape depends on the consolidation pressure p
c
s. The expansion of the static yield surface, which describes the hardening of the material, is expressed by the variation of the consolidation pressure due to the inelastic volumetric strain ε
v
v
p as follows:
(8)
p
c
s
=
p
0
·
exp
⁡
(
1
+
e
0
λ
-
κ
ε
v
v
p
)
.Figure 2
Yield criterion of model.In (8), p
0 is the intercept of initial static yield surface in the p
′ axis, for soft clay, it can be assigned to the biggest consolidation stress in history; λ is the slope of normal consolidation curve in e
~
ln
⁡
p
′ plane; κ is the slope of recompression curve in e
~
ln
⁡
p
′ plane.A dynamic yield criterionF
d is defined to describe the current state of stress as follows:
(9)
F
d
=
q
2
M
2
+
p
′
(
p
′
-
p
c
d
)
=
0
.From (8) and (9), the f and f
0 in (7) can be written as
(10)
f
=
p
c
d
=
p
′
+
q
2
M
2
p
′
,
f
0
=
p
c
s
=
p
0
·
exp
⁡
(
1
+
e
0
λ
-
κ
ε
v
)
.So, we can write (7) as follows:
(11)
ϕ
(
f
)
=
p
c
d
p
c
s
-
1
.Combining (6)~(11), (6) can be written as follows:
(12)
ε
˙
i
j
v
p
=
1
H
2
〈
ϕ
(
f
)
〉
(
3
s
i
j
′
M
2
+
(
2
p
′
-
p
c
d
)
δ
i
j
3
)
.From (4), (5), and (12), when the stress is constant, the total strain of the presented model can be written as follows:
(13)
ε
i
j
=
I
1
9
K
δ
i
j
+
s
i
j
′
2
G
0
+
s
i
j
′
2
G
1
(
1
-
e
-
(
E
1
/
η
1
)
t
)
+
1
H
2
〈
ϕ
(
f
)
〉
(
3
s
i
j
′
M
2
+
(
2
p
′
-
p
c
d
)
δ
i
j
3
)
t
.


## 3. Experiments and Computation

### 3.1. Test Program

In order to verify the presented model, laboratory triaxial rheological tests were performed to observe the mechanical behavior of soft clay under long-term loading. The soil investigated in the experiments is Ningbo soft clay, which is a kind of problematic soil for low strength, high compressibility, and time-dependent behavior and deposits in Hangzhou Bay, China. The basic properties of the soil were summarized in Table1. Both moisture content and moist unit weight of this material are significantly less than those of typical natural sedimentary deposits.Table 1
Property of testing material.
Property
Value
Embedded depth (m)
3.5~4.5
Specific gravity
2.67
Moist unit weight (kN/m3)
17.3
Water content (%)
40.4
Void ratio
1.08
Liquid limit (%)
42.5
Plastic limit (%)
24.9
Plasticity index
17.6
Slope of the critical state line
0.898
Slope of compression curve
0.102
Slope of the recompression curve
0.015The experimental equipment was refitted on the platform of strain controlled triaxial apparatus, the original strain loading method was changed to stress loading method, while the confining pressure system, the back pressure loading system, and the measurement system were reserved. The soil was cut into replicate specimens with a diameter of 39.8 mm and a height of 80.0 mm and placed in the triaxial test apparatus. Both soil specimens were consolidated 24 hours under a confining pressure of 100 kPa and 200 kPa, respectively; then, keep the confining pressure as constant and apply the deviatoric stress increment until the specimens were damaged or the total strain exceeds 15%. In the whole process, the free drainage conditions were kept. Table2 showed the different loading rates applied in the experiments. In order to make sure that the creep deformation can be fully developed, each load increment lasted no less than 7 days until the deformation of the specimen is lower than 0.01 mm within 24 h.Table 2
Loading scheme of triaxial rheological tests.
Specimen number
Confining pressure (kPa)
Loading scheme (kPa)
1
100
30-60-90-12-150-180-210-240
2
200
40-80-120-160-200-240-280-320


### 3.2. Experimental Results

The experimental results were presented to observe the rheological behavior of the soft clay. It is considered that deformation occurs during the whole process, but only the deformation that occurs after the dissipation of the excessive pore water pressure can be regarded as creep only. The full shearing process stress-strain-time curve of each specimen is illustrated in Figure3, and the mutistage strain-time curves of each specimen are illustrated in Figure 4. The figures show that the creep process under triaxial loading exhibits attenuation characteristics when the deviatoric stress is small, which contains two stages of creep: primary stage and steady stage; while when the deviatoric stress is large, it may exhibit accelerated characteristics, which contains three stages of creep: primary stage, steady stage and accelerated stage, such as the last load increment of specimen 2.Full-process press-strain-time curve.
(a)
Specimen 1
(b)
Specimen 2Mutistage strain-time curves.
(a)
Specimen 1
(b)
Specimen 2


### 3.3. Computation and Curve Fitting

For triaxial rheological tests,I
1
=
σ
1
+
2
σ
3
=
3
p
′, s
11
′
=
2
(
σ
1
-
σ
3
)
/
3
=
2
q
/
3, so (13) can be written as follows:
(14)
ε
11
=
P
1
+
P
2
(
1
-
e
-
P
3
t
)
+
P
4
t
,
where, P
1
=
(
p
/
3
K
)
+
(
q
/
3
G
0
), P
2
=
(
q
/
3
G
1
), P
3
=
(
G
1
/
H
1
), P
4
=
(
1
/
H
2
)
〈
ϕ
(
f
)
〉
(
(
2
q
/
M
2
)
+
(
1
/
3
)
(
2
p
′
-
p
c
d
)
).Use (14) to fit the experimental data, the results were showed in Figures 5 and 6, in which the parameters SSE and R
2 denote the residual sum of squares and squares of correlation coefficient, respectively. From Figures 5 and 6, we can see that the model fits testing data of primary creep stage and steady creep stage very well, while for the accelerated creep stage, it needs further study. According to the value of P
1, P
2, P
3, and P
4, the value of the model parameters E
0, E
1, η
1, and η
2 can be determined easily, such as in Tables 3 and 4.Table 3
Model parameters of specimen 1 (σ
3
=
100 kPa).
σ
1
-
σ
3 (kPa)
E
0 (MPa)
E
1 (MPa)
η
1 (MPa·s)
η
2 (GPa·s)
30
2.49
18.56
501.60
16073
60
2.51
10.69
189.93
17950
90
2.49
11.69
181.98
21297
120
2.47
11.58
145.44
24253
150
2.36
10.49
104.97
25049
180
2.36
10.18
98.91
22046
210
2.35
9.45
120.32
26240Table 4
Model parameters of specimen 2 (σ
3
=
200 kPa).
σ
1
-
σ
3 (kPa)
E
0 (MPa)
E
1 (MPa)
η
1 (MPa·s)
η
2 (GPa·s)
40
10.75
18.04
138.81
12223
80
9.03
18.70
258.62
11747
120
8.34
16.94
299.47
12591
160
7.63
16.92
300.57
12704
200
6.60
17.29
261.18
11732
240
6.04
18.31
274.81
12823
280
5.08
16.02
281.65
11589Fitting curves of specimen 1 (σ
3 = 100 kPa).
(a)
q
=
30 kPa
(b)
q
=
60 kPa
(c)
q
=
90 kPa
(d)
q
=
120 kPa
(e)
q
=
150 kPa
(f)
q
=
180 kPa
(g)
q
=
210 kPaFitting curves of specimen 2 (σ
3 = 200 kPa).
(a)
q
=
40 kPa
(b)
q
=
80 kPa
(c)
q
=
120 kPa
(d)
q
=
160 kPa
(e)
q
=
200 kPa
(f)
q
=
240 kPa
(g)
q
=
280 kPa


### 3.4. Discussion

Table3 shows that the model parameters of specimen 1 (σ
3 = 100 kPa) change slightly with the change of deviatoric stress except the first stage; Table 4 shows that the model parameters of specimen 2 (σ
3 = 200 kPa) change slightly with the change of deviatoric stress except E
0. Compared with the model parameters of specimen 1 and specimen 2, it can be found that both E
0 and E
1 increase notably with the increase of σ
3, which means that the elastic strain reduces when σ
3 increases, and η
1 increases notably while η
2 is reduced notably with the increase of σ
3, which means that the viscoelastic strain rate and viscoplastic strain rate change significantly when σ
3 increases. For the lack of more experiment data, the further variation of model parameters was unable to find out, which needed further study.Based on the parameters in Tables3 and 4, the components of elastic strain, viscoelastic strain, and viscoplastic strain at any time can be calculated with (13). Tables 5 and 6 show ε
11
e, ε
11
v
e, and ε
11
v
p of all loading levels when t = 7 d; it can be found that the elastic strain accounted for the major part of the total strain, but the viscoelastic strain and viscoplastic strain increase with the increase of σ
3, which takes about 15%~30% for each loading increment of the total and cannot be neglected.Table 5
Percentage of strain components of specimen 1 (σ
3
=
100 kPa, t
=
7 d).
σ
1
-
σ
3 (kPa)
ε
11
e (%)
ε
11
ve (%)
ε
11
vp (%)
30
92.06
5.56
2.38
60
85.30
12.02
2.68
90
84.83
12.21
2.96
120
82.98
12.72
4.30
150
81.51
13.76
4.73
180
81.09
14.48
4.43
210
80.41
15.73
3.85Table 6
Percentage of strain components of specimen 2 (σ
3
=
200 kPa, t
=
7 d).
σ
1
-
σ
3 (kPa)
ε
11
e (%)
ε
11
ve (%)
ε
11
vp (%)
40
77.03
16.53
6.45
80
67.74
16.83
15.43
120
66.77
19.71
13.53
160
67.13
19.81
13.06
200
69.51
18.37
12.12
240
70.96
16.86
12.17
280
73.97
17.40
8.63Figure7 shows several typical curves of the development of computed viscoelastic strain and viscoplastic strain of specimen 1 and specimen 2. It can be found that the viscoelastic strain increases quickly in the early stage and then reaches steady in a short time (about 30 hours), while the viscoplastic strain increases the whole process, which indicates that the long-term deformation after consolidation of soft clay foundation is mainly caused by the viscoplasticity rather than viscoelasticityDevelopment of computed viscous strain (VES: viscoelastic strain; VPS: viscoplastic strain).
(a)
Specimen 1
(b)
Specimen 2


## 3.1. Test Program

In order to verify the presented model, laboratory triaxial rheological tests were performed to observe the mechanical behavior of soft clay under long-term loading. The soil investigated in the experiments is Ningbo soft clay, which is a kind of problematic soil for low strength, high compressibility, and time-dependent behavior and deposits in Hangzhou Bay, China. The basic properties of the soil were summarized in Table1. Both moisture content and moist unit weight of this material are significantly less than those of typical natural sedimentary deposits.Table 1
Property of testing material.
Property
Value
Embedded depth (m)
3.5~4.5
Specific gravity
2.67
Moist unit weight (kN/m3)
17.3
Water content (%)
40.4
Void ratio
1.08
Liquid limit (%)
42.5
Plastic limit (%)
24.9
Plasticity index
17.6
Slope of the critical state line
0.898
Slope of compression curve
0.102
Slope of the recompression curve
0.015The experimental equipment was refitted on the platform of strain controlled triaxial apparatus, the original strain loading method was changed to stress loading method, while the confining pressure system, the back pressure loading system, and the measurement system were reserved. The soil was cut into replicate specimens with a diameter of 39.8 mm and a height of 80.0 mm and placed in the triaxial test apparatus. Both soil specimens were consolidated 24 hours under a confining pressure of 100 kPa and 200 kPa, respectively; then, keep the confining pressure as constant and apply the deviatoric stress increment until the specimens were damaged or the total strain exceeds 15%. In the whole process, the free drainage conditions were kept. Table2 showed the different loading rates applied in the experiments. In order to make sure that the creep deformation can be fully developed, each load increment lasted no less than 7 days until the deformation of the specimen is lower than 0.01 mm within 24 h.Table 2
Loading scheme of triaxial rheological tests.
Specimen number
Confining pressure (kPa)
Loading scheme (kPa)
1
100
30-60-90-12-150-180-210-240
2
200
40-80-120-160-200-240-280-320


## 3.2. Experimental Results

The experimental results were presented to observe the rheological behavior of the soft clay. It is considered that deformation occurs during the whole process, but only the deformation that occurs after the dissipation of the excessive pore water pressure can be regarded as creep only. The full shearing process stress-strain-time curve of each specimen is illustrated in Figure3, and the mutistage strain-time curves of each specimen are illustrated in Figure 4. The figures show that the creep process under triaxial loading exhibits attenuation characteristics when the deviatoric stress is small, which contains two stages of creep: primary stage and steady stage; while when the deviatoric stress is large, it may exhibit accelerated characteristics, which contains three stages of creep: primary stage, steady stage and accelerated stage, such as the last load increment of specimen 2.Full-process press-strain-time curve.
(a)
Specimen 1
(b)
Specimen 2Mutistage strain-time curves.
(a)
Specimen 1
(b)
Specimen 2


## 3.3. Computation and Curve Fitting

For triaxial rheological tests,I
1
=
σ
1
+
2
σ
3
=
3
p
′, s
11
′
=
2
(
σ
1
-
σ
3
)
/
3
=
2
q
/
3, so (13) can be written as follows:
(14)
ε
11
=
P
1
+
P
2
(
1
-
e
-
P
3
t
)
+
P
4
t
,
where, P
1
=
(
p
/
3
K
)
+
(
q
/
3
G
0
), P
2
=
(
q
/
3
G
1
), P
3
=
(
G
1
/
H
1
), P
4
=
(
1
/
H
2
)
〈
ϕ
(
f
)
〉
(
(
2
q
/
M
2
)
+
(
1
/
3
)
(
2
p
′
-
p
c
d
)
).Use (14) to fit the experimental data, the results were showed in Figures 5 and 6, in which the parameters SSE and R
2 denote the residual sum of squares and squares of correlation coefficient, respectively. From Figures 5 and 6, we can see that the model fits testing data of primary creep stage and steady creep stage very well, while for the accelerated creep stage, it needs further study. According to the value of P
1, P
2, P
3, and P
4, the value of the model parameters E
0, E
1, η
1, and η
2 can be determined easily, such as in Tables 3 and 4.Table 3
Model parameters of specimen 1 (σ
3
=
100 kPa).
σ
1
-
σ
3 (kPa)
E
0 (MPa)
E
1 (MPa)
η
1 (MPa·s)
η
2 (GPa·s)
30
2.49
18.56
501.60
16073
60
2.51
10.69
189.93
17950
90
2.49
11.69
181.98
21297
120
2.47
11.58
145.44
24253
150
2.36
10.49
104.97
25049
180
2.36
10.18
98.91
22046
210
2.35
9.45
120.32
26240Table 4
Model parameters of specimen 2 (σ
3
=
200 kPa).
σ
1
-
σ
3 (kPa)
E
0 (MPa)
E
1 (MPa)
η
1 (MPa·s)
η
2 (GPa·s)
40
10.75
18.04
138.81
12223
80
9.03
18.70
258.62
11747
120
8.34
16.94
299.47
12591
160
7.63
16.92
300.57
12704
200
6.60
17.29
261.18
11732
240
6.04
18.31
274.81
12823
280
5.08
16.02
281.65
11589Fitting curves of specimen 1 (σ
3 = 100 kPa).
(a)
q
=
30 kPa
(b)
q
=
60 kPa
(c)
q
=
90 kPa
(d)
q
=
120 kPa
(e)
q
=
150 kPa
(f)
q
=
180 kPa
(g)
q
=
210 kPaFitting curves of specimen 2 (σ
3 = 200 kPa).
(a)
q
=
40 kPa
(b)
q
=
80 kPa
(c)
q
=
120 kPa
(d)
q
=
160 kPa
(e)
q
=
200 kPa
(f)
q
=
240 kPa
(g)
q
=
280 kPa


## 3.4. Discussion

Table3 shows that the model parameters of specimen 1 (σ
3 = 100 kPa) change slightly with the change of deviatoric stress except the first stage; Table 4 shows that the model parameters of specimen 2 (σ
3 = 200 kPa) change slightly with the change of deviatoric stress except E
0. Compared with the model parameters of specimen 1 and specimen 2, it can be found that both E
0 and E
1 increase notably with the increase of σ
3, which means that the elastic strain reduces when σ
3 increases, and η
1 increases notably while η
2 is reduced notably with the increase of σ
3, which means that the viscoelastic strain rate and viscoplastic strain rate change significantly when σ
3 increases. For the lack of more experiment data, the further variation of model parameters was unable to find out, which needed further study.Based on the parameters in Tables3 and 4, the components of elastic strain, viscoelastic strain, and viscoplastic strain at any time can be calculated with (13). Tables 5 and 6 show ε
11
e, ε
11
v
e, and ε
11
v
p of all loading levels when t = 7 d; it can be found that the elastic strain accounted for the major part of the total strain, but the viscoelastic strain and viscoplastic strain increase with the increase of σ
3, which takes about 15%~30% for each loading increment of the total and cannot be neglected.Table 5
Percentage of strain components of specimen 1 (σ
3
=
100 kPa, t
=
7 d).
σ
1
-
σ
3 (kPa)
ε
11
e (%)
ε
11
ve (%)
ε
11
vp (%)
30
92.06
5.56
2.38
60
85.30
12.02
2.68
90
84.83
12.21
2.96
120
82.98
12.72
4.30
150
81.51
13.76
4.73
180
81.09
14.48
4.43
210
80.41
15.73
3.85Table 6
Percentage of strain components of specimen 2 (σ
3
=
200 kPa, t
=
7 d).
σ
1
-
σ
3 (kPa)
ε
11
e (%)
ε
11
ve (%)
ε
11
vp (%)
40
77.03
16.53
6.45
80
67.74
16.83
15.43
120
66.77
19.71
13.53
160
67.13
19.81
13.06
200
69.51
18.37
12.12
240
70.96
16.86
12.17
280
73.97
17.40
8.63Figure7 shows several typical curves of the development of computed viscoelastic strain and viscoplastic strain of specimen 1 and specimen 2. It can be found that the viscoelastic strain increases quickly in the early stage and then reaches steady in a short time (about 30 hours), while the viscoplastic strain increases the whole process, which indicates that the long-term deformation after consolidation of soft clay foundation is mainly caused by the viscoplasticity rather than viscoelasticityDevelopment of computed viscous strain (VES: viscoelastic strain; VPS: viscoplastic strain).
(a)
Specimen 1
(b)
Specimen 2


## 4. Derivation and Verification of UMAT

### 4.1. Derivation of UMAT

ABAQUS is a FEM software which has powerful computing capabilities of strongly nonlinear problems. A user material subroutine (UMAT) is coded for the proposed model on the platform of ABAQUS to study the elastic-viscoplastic consolidation properties of soft clay. According to ABAQUS, the stiffness coefficient matrix (Jacobian matrix) and the equation to update the stress and strain should be defined in the UAMT. The Jacobian matrixD is defined as follows:
(15)
D
=
∂
Δ
σ
∂
Δ
ε
,
where Δ
σ and Δ
ε denote the stress increment matrix and strain increment matrix, respectively.For the proposed model, (4)~(6) can be written as follows in 2D or 3D problems:(16a)
ε
n
e
=
μ
E
0
A
σ
n
,
(16b)
ε
˙
n
v
e
=
1
η
1
A
σ
n
-
G
1
η
1
ε
n
v
e
,
(16c)
ε
˙
n
v
p
=
1
η
2
〈
ϕ
(
F
)
〉
∂
F
∂
σ
.In (16a), (16b), and (16c) σ
n, ε
n, ε
n
v
e, and ε
n
v
p denote the Voigt matrix form of stress and strain component tensors; A denotes the flexibility matrix, for plane strain problem and 3D problems, and it can be written as follows:(17a)
A
=
[
1
-
μ
-
μ
0
-
μ
1
-
μ
0
-
μ
-
μ
1
0
0
0
0
2
(
1
+
μ
)
]
,
(17b)
A
=
[
1
-
μ
-
μ
0
0
0
-
μ
1
-
μ
0
0
0
-
μ
-
μ
1
0
0
0
0
0
0
2
(
1
+
μ
)
0
0
0
0
0
0
2
(
1
+
μ
)
0
0
0
0
0
0
2
(
1
+
μ
)
]
.Suppose that the time step of thenth increment is Δ
t
n, so the viscoelastic strain increment in the increment step can be written as follows:
(18)
Δ
ε
n
v
e
=
Δ
t
n
[
(
1
-
Θ
)
ε
˙
n
v
e
+
Θ
ε
˙
n
+
1
v
e
]
,
where Θ is differential coefficient, 0
≤
Θ
≤
1. Suppose that B
n
=
1
/
(
1
+
Θ
Δ
t
n
E
1
/
η
1
), C
n
=
B
n
Θ
Δ
t
n
A
/
η
1; (18) can be written as follows:
(19)
Δ
ε
n
v
e
=
B
n
ε
˙
n
v
e
Δ
t
n
+
C
n
Δ
σ
n
.The viscoplastic strain increment in the increment step can be written as follows:(20)
Δ
ε
n
v
p
=
Δ
t
n
[
(
1
-
Θ
)
ε
˙
n
v
p
+
Θ
ε
˙
n
+
1
v
p
]
.In order to calculate the viscoplastic strain rate at the end of thenth increment, apply the Taylor series expansion and ignore the second order trace; then we can obtain the following equation:
(21)
ε
˙
n
+
1
v
p
=
ε
˙
n
v
p
+
∂
ε
˙
n
v
p
∂
σ
∂
σ
∂
t
Δ
t
n
.Suppose thatH
n
=
∂
ε
˙
n
v
p
/
∂
σ; (21) can be written as
(22)
ε
˙
n
+
1
v
p
=
ε
˙
n
v
p
+
H
n
Δ
σ
n
.Substituting (22) into (21), we can obtain the following equation:
(23)
Δ
ε
n
v
p
=
Δ
t
n
ε
˙
n
v
p
+
Θ
Δ
t
n
H
n
Δ
σ
n
.From (16c), H
n can be written as
(24)
H
n
=
1
η
2
[
〈
ϕ
(
F
)
〉
∂
a
T
∂
σ
+
∂
〈
ϕ
(
F
)
〉
∂
F
a
a
T
]
,
where a
=
∂
F
/
∂
σ. According to (9), a can be written as follows:
(25)
a
=
∂
F
∂
σ
=
C
1
a
1
+
C
2
a
2
,
where C
1
=
2
q
/
M
2, C
2
=
2
p
-
p
0, a
1
=
∂
q
/
∂
σ, a
2
=
∂
p
/
∂
σ.From (25), the following equation can be derived:
(26)
∂
a
T
∂
σ
=
∂
C
1
∂
q
a
1
a
1
T
+
C
1
∂
a
1
T
∂
σ
+
∂
C
2
∂
p
a
2
a
2
T
,
where (
∂
a
1
T
/
∂
σ
)
=
(
3
/
2
q
)
M
1
-
(
9
/
4
q
3
)
M
2, a
1
a
1
T
=
(
9
/
4
q
2
)
M
2, a
2
a
2
T
=
(
1
/
9
)
M
3. M
1, M
2, and M
3 are matrices written as follows:
(27)
M
1
=
[
2
/
3
-
1
/
3
-
1
/
3
0
0
0
-
1
/
3
2
/
3
-
1
/
3
0
0
0
-
1
/
3
-
1
/
3
2
/
3
0
0
0
0
0
0
2
0
0
0
0
0
0
2
0
0
0
0
0
0
2
]
,
M
2
=
[
S
11
2
S
11
S
22
S
11
S
33
2
S
11
S
12
2
S
11
S
23
2
S
11
S
13
S
11
S
22
S
22
2
S
22
S
33
2
S
22
S
12
2
S
22
S
23
2
S
22
S
13
S
11
S
33
S
22
S
33
S
33
2
2
S
33
S
12
2
S
33
S
23
2
S
33
S
13
2
S
11
S
12
2
S
22
S
12
2
S
33
S
12
4
S
12
2
4
S
12
S
23
4
S
12
S
13
2
S
11
S
23
2
S
22
S
23
2
S
33
S
23
4
S
12
S
23
4
S
23
2
4
S
23
S
13
2
S
11
S
13
2
S
22
S
13
2
S
33
S
13
4
S
12
S
13
4
S
23
S
13
4
S
13
2
]
,
M
3
=
[
1
1
1
0
0
0
1
1
1
0
0
0
1
1
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
]
.Combining the above equations, the value ofH
n can be derived. Substituting H
n into (22), the viscoplastic strain increment Δ
ε
n
v
p is obtained.For current increment,Δ
σ
n is calculated by (28) as follows:
(28)
Δ
σ
n
=
D
e
Δ
ε
n
e
.Substituting (19) and (23) into (28), we can obtain the following equation:
(29)
Δ
σ
n
=
D
(
Δ
ε
n
-
B
Δ
t
n
ε
˙
n
v
e
-
Δ
t
n
ε
˙
n
v
p
)
.Equation (29) is the relationship between stress increment and strain increment, in which matrix D is defined as follows:
(30)
D
=
[
(
D
e
)
-
1
+
C
n
+
Θ
Δ
t
n
H
n
]
-
1
.Based on the above equations, a user material subroutine was coded for the proposed constitutive model on the platform of ABAQUS.


### 4.2. Verification of UMAT

In order to verify the validity of the UMAT, we study a soil column of triaxial compression such as Figure8, the size of the column is 0.1 m × 0.1 m × 0.2 m. Supposing that σ
1 = 100 kPa, σ
2
=
σ
3
=
50 kPa, and the model parameters are given in Table 7. The column was cut into one element, the element type is C3D8, and the consolidation effect was not considered here. Embedding the UMAT coded in this paper on the platform of ABAQUS, we calculated the triaxial compression property of the column; the results were summarized in Figure 9~Figure 10.Table 7
Material parameters of verification model.
E
0 (MPa)
E
1 (MPa)
η
1 (GPa·s)
η
2 (GPa·s)
P
0 (kPa)
e
0
μ
M
λ
κ
1.0
2.0
1000
1000
50
1.08
0.35
0.898
0.102
0.015Figure 8
Sketch of verification model.Triaxial compression property of soil column.
(a)
(b)
(c)
(d)Sensitivity of model parameters.
(a)
(b)
(c)From Figure9(a), we can see the deformation of the column exhibits apparent rheology, the rheological deformation takes about 37% of the whole deformation. Figure 9(b) shows that the static yield surface gradually expanded with the time until it reached steady state, when it coincided with the dynamic yield surface. From Figures 9(c) and 9(d), we can see that the viscoelastic strain and the viscoplastic strain both gradually increased with time until they reached their steady state; accordingly, the viscoelastic strain rate and the viscoplastic strain rate were gradually deduced to nearly zero. The above results fit the property of Nishihara model and the hardening law of modified Cambridge model very well, which proved the correctness of the UMAT.In the proposed model,η
1, η
2, and p
0 are the most important parameters; Figure 10 showed the sensitivity of these parameters. From Figures 10(a) and 10(b), we can see the final viscous deformation is a constant regardless of the different values of η
1, η
2, but the time when it reaches steady varies greatly, so, the accuracy of η
1 and η
2 is directly related to the accuracy of computed displacement in the rheological process. Figure 10(c) shows that the value of p
0 has a great influence on the viscoplastic strain: the smaller p
0 is, the more rapidly the viscoplastic strain develops and the larger the final viscoplastic strain occurs, and so, the accuracy of p
0 is not only directly related to the accuracy of computed displacement in the rheological process, but also affected the final displacement, which needs to be paid more attention.


## 4.1. Derivation of UMAT

ABAQUS is a FEM software which has powerful computing capabilities of strongly nonlinear problems. A user material subroutine (UMAT) is coded for the proposed model on the platform of ABAQUS to study the elastic-viscoplastic consolidation properties of soft clay. According to ABAQUS, the stiffness coefficient matrix (Jacobian matrix) and the equation to update the stress and strain should be defined in the UAMT. The Jacobian matrixD is defined as follows:
(15)
D
=
∂
Δ
σ
∂
Δ
ε
,
where Δ
σ and Δ
ε denote the stress increment matrix and strain increment matrix, respectively.For the proposed model, (4)~(6) can be written as follows in 2D or 3D problems:(16a)
ε
n
e
=
μ
E
0
A
σ
n
,
(16b)
ε
˙
n
v
e
=
1
η
1
A
σ
n
-
G
1
η
1
ε
n
v
e
,
(16c)
ε
˙
n
v
p
=
1
η
2
〈
ϕ
(
F
)
〉
∂
F
∂
σ
.In (16a), (16b), and (16c) σ
n, ε
n, ε
n
v
e, and ε
n
v
p denote the Voigt matrix form of stress and strain component tensors; A denotes the flexibility matrix, for plane strain problem and 3D problems, and it can be written as follows:(17a)
A
=
[
1
-
μ
-
μ
0
-
μ
1
-
μ
0
-
μ
-
μ
1
0
0
0
0
2
(
1
+
μ
)
]
,
(17b)
A
=
[
1
-
μ
-
μ
0
0
0
-
μ
1
-
μ
0
0
0
-
μ
-
μ
1
0
0
0
0
0
0
2
(
1
+
μ
)
0
0
0
0
0
0
2
(
1
+
μ
)
0
0
0
0
0
0
2
(
1
+
μ
)
]
.Suppose that the time step of thenth increment is Δ
t
n, so the viscoelastic strain increment in the increment step can be written as follows:
(18)
Δ
ε
n
v
e
=
Δ
t
n
[
(
1
-
Θ
)
ε
˙
n
v
e
+
Θ
ε
˙
n
+
1
v
e
]
,
where Θ is differential coefficient, 0
≤
Θ
≤
1. Suppose that B
n
=
1
/
(
1
+
Θ
Δ
t
n
E
1
/
η
1
), C
n
=
B
n
Θ
Δ
t
n
A
/
η
1; (18) can be written as follows:
(19)
Δ
ε
n
v
e
=
B
n
ε
˙
n
v
e
Δ
t
n
+
C
n
Δ
σ
n
.The viscoplastic strain increment in the increment step can be written as follows:(20)
Δ
ε
n
v
p
=
Δ
t
n
[
(
1
-
Θ
)
ε
˙
n
v
p
+
Θ
ε
˙
n
+
1
v
p
]
.In order to calculate the viscoplastic strain rate at the end of thenth increment, apply the Taylor series expansion and ignore the second order trace; then we can obtain the following equation:
(21)
ε
˙
n
+
1
v
p
=
ε
˙
n
v
p
+
∂
ε
˙
n
v
p
∂
σ
∂
σ
∂
t
Δ
t
n
.Suppose thatH
n
=
∂
ε
˙
n
v
p
/
∂
σ; (21) can be written as
(22)
ε
˙
n
+
1
v
p
=
ε
˙
n
v
p
+
H
n
Δ
σ
n
.Substituting (22) into (21), we can obtain the following equation:
(23)
Δ
ε
n
v
p
=
Δ
t
n
ε
˙
n
v
p
+
Θ
Δ
t
n
H
n
Δ
σ
n
.From (16c), H
n can be written as
(24)
H
n
=
1
η
2
[
〈
ϕ
(
F
)
〉
∂
a
T
∂
σ
+
∂
〈
ϕ
(
F
)
〉
∂
F
a
a
T
]
,
where a
=
∂
F
/
∂
σ. According to (9), a can be written as follows:
(25)
a
=
∂
F
∂
σ
=
C
1
a
1
+
C
2
a
2
,
where C
1
=
2
q
/
M
2, C
2
=
2
p
-
p
0, a
1
=
∂
q
/
∂
σ, a
2
=
∂
p
/
∂
σ.From (25), the following equation can be derived:
(26)
∂
a
T
∂
σ
=
∂
C
1
∂
q
a
1
a
1
T
+
C
1
∂
a
1
T
∂
σ
+
∂
C
2
∂
p
a
2
a
2
T
,
where (
∂
a
1
T
/
∂
σ
)
=
(
3
/
2
q
)
M
1
-
(
9
/
4
q
3
)
M
2, a
1
a
1
T
=
(
9
/
4
q
2
)
M
2, a
2
a
2
T
=
(
1
/
9
)
M
3. M
1, M
2, and M
3 are matrices written as follows:
(27)
M
1
=
[
2
/
3
-
1
/
3
-
1
/
3
0
0
0
-
1
/
3
2
/
3
-
1
/
3
0
0
0
-
1
/
3
-
1
/
3
2
/
3
0
0
0
0
0
0
2
0
0
0
0
0
0
2
0
0
0
0
0
0
2
]
,
M
2
=
[
S
11
2
S
11
S
22
S
11
S
33
2
S
11
S
12
2
S
11
S
23
2
S
11
S
13
S
11
S
22
S
22
2
S
22
S
33
2
S
22
S
12
2
S
22
S
23
2
S
22
S
13
S
11
S
33
S
22
S
33
S
33
2
2
S
33
S
12
2
S
33
S
23
2
S
33
S
13
2
S
11
S
12
2
S
22
S
12
2
S
33
S
12
4
S
12
2
4
S
12
S
23
4
S
12
S
13
2
S
11
S
23
2
S
22
S
23
2
S
33
S
23
4
S
12
S
23
4
S
23
2
4
S
23
S
13
2
S
11
S
13
2
S
22
S
13
2
S
33
S
13
4
S
12
S
13
4
S
23
S
13
4
S
13
2
]
,
M
3
=
[
1
1
1
0
0
0
1
1
1
0
0
0
1
1
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
]
.Combining the above equations, the value ofH
n can be derived. Substituting H
n into (22), the viscoplastic strain increment Δ
ε
n
v
p is obtained.For current increment,Δ
σ
n is calculated by (28) as follows:
(28)
Δ
σ
n
=
D
e
Δ
ε
n
e
.Substituting (19) and (23) into (28), we can obtain the following equation:
(29)
Δ
σ
n
=
D
(
Δ
ε
n
-
B
Δ
t
n
ε
˙
n
v
e
-
Δ
t
n
ε
˙
n
v
p
)
.Equation (29) is the relationship between stress increment and strain increment, in which matrix D is defined as follows:
(30)
D
=
[
(
D
e
)
-
1
+
C
n
+
Θ
Δ
t
n
H
n
]
-
1
.Based on the above equations, a user material subroutine was coded for the proposed constitutive model on the platform of ABAQUS.


## 4.2. Verification of UMAT

In order to verify the validity of the UMAT, we study a soil column of triaxial compression such as Figure8, the size of the column is 0.1 m × 0.1 m × 0.2 m. Supposing that σ
1 = 100 kPa, σ
2
=
σ
3
=
50 kPa, and the model parameters are given in Table 7. The column was cut into one element, the element type is C3D8, and the consolidation effect was not considered here. Embedding the UMAT coded in this paper on the platform of ABAQUS, we calculated the triaxial compression property of the column; the results were summarized in Figure 9~Figure 10.Table 7
Material parameters of verification model.
E
0 (MPa)
E
1 (MPa)
η
1 (GPa·s)
η
2 (GPa·s)
P
0 (kPa)
e
0
μ
M
λ
κ
1.0
2.0
1000
1000
50
1.08
0.35
0.898
0.102
0.015Figure 8
Sketch of verification model.Triaxial compression property of soil column.
(a)
(b)
(c)
(d)Sensitivity of model parameters.
(a)
(b)
(c)From Figure9(a), we can see the deformation of the column exhibits apparent rheology, the rheological deformation takes about 37% of the whole deformation. Figure 9(b) shows that the static yield surface gradually expanded with the time until it reached steady state, when it coincided with the dynamic yield surface. From Figures 9(c) and 9(d), we can see that the viscoelastic strain and the viscoplastic strain both gradually increased with time until they reached their steady state; accordingly, the viscoelastic strain rate and the viscoplastic strain rate were gradually deduced to nearly zero. The above results fit the property of Nishihara model and the hardening law of modified Cambridge model very well, which proved the correctness of the UMAT.In the proposed model,η
1, η
2, and p
0 are the most important parameters; Figure 10 showed the sensitivity of these parameters. From Figures 10(a) and 10(b), we can see the final viscous deformation is a constant regardless of the different values of η
1, η
2, but the time when it reaches steady varies greatly, so, the accuracy of η
1 and η
2 is directly related to the accuracy of computed displacement in the rheological process. Figure 10(c) shows that the value of p
0 has a great influence on the viscoplastic strain: the smaller p
0 is, the more rapidly the viscoplastic strain develops and the larger the final viscoplastic strain occurs, and so, the accuracy of p
0 is not only directly related to the accuracy of computed displacement in the rheological process, but also affected the final displacement, which needs to be paid more attention.


## 5. Analysis on Rheological Consolidation of Soft Clay

### 5.1. Computed Model

Consolidation problems can be classified as small strain theory and finite strain theory according to the difference of geometric assumption. Small strain theory assumed that the strain is little, but the consolidation deformation is usually very large for deep soft clay foundation. Take dredger fill for example, the deformation may exceed 80%, for which the small strain theory is no longer suitable [14, 15]. As the pioneer researchers studied on finite strain consolidation theory, Mikasa [16] and Gibson et al. [17] derived their one-dimensional finite strain consolidation equations, respectively, in 1960s, and many scholars developed further research on the topic in the following decades. Because of strong nonlinearity, most of existing researches on finite strain consolidation adopted some simplified assumptions without considering rheological characteristics. By means of ABAQUS, this paper studied the rheological consolidation of soft clay considering finite strain effect.The computed model is plane strain problem as shown in Figure11. The soil foundation is homogeneous, anisotropic, saturated, and normally consolidated soft clay with a thickness of 10.0 m, and only the top boundary is permeable. In addition, the semilogarithmic empirical formula proposed by Tavenas et al. [18] was adopted to consider the variation of permeability in the process of consolidation as follows:
(31)
k
v
=
k
v
0
·
1
0
(
e
-
e
0
)
/
C
k
,
where k
v
0 is the initial permeability coefficients, and k
v is the permeability coefficient in the consolidation process; e
0 is the initial void ratio, and e is the void ratio in the consolidation process; C
k is the permeability index. According to the above work, the model parameters are summarized in Table 8 approximately.Table 8
Material parameters of computed model.
E
0 (MPa)
E
1 (MPa)
η
1 (MPa·s)
η
2 (GPa·s)
γ
sat
⁡ (kN·m−3)
k
v
0 (m·s−1)
P
0 (kPa)
M
λ
κ
e
0
μ
C
k
2.5
11.5
150
24000
17.7
4.5
E
-
9
50
0.898
0.102
0.015
1.08
0.35
0.5Figure 11
Calculation model and mesh generation.


### 5.2. Consolidation Property of Soft Clay

From Figure12 we can see the differences of the consolidation process between finite strain and small strain. Figure 12(a) shows that for any location the excess pore pressure of finite strain is less than that of small strain at the same time, and the difference is increasing with the increase of time in beginning period. Figure 12(b) shows that the displacement of finite strain is less than that of small strain and the difference increased significantly if the load is increasing. For example, when q = 200 kPa, the difference may take about 20% while it is nearly 0 when q = 50 kPa. So, it is essential to consider the error caused by small strain if the load is large and the compressibility of the soil is large.Results comparison of large strain and small strain (SS: small train; FS: finite strain).
(a)
(b)From Figure13 we can see the influence of rheology on the consolidation process. Figure 13(a) shows that the dissipation of excess pore pressure is significantly slower if the rheological effect is considered; for example, when t = 1000 d, the excess pore pressure is nearly 0 if the rheological effect is not considered, but it is still larger than 20 kPa if the rheological effect is considered. Figure 13(a) shows that the displacement is influenced by rheological effect too; for example, in the early stage (t ≤ 100 d), the displacement is nearly the same for two cases, but the final displacement may have a difference of more than 35%, which indicates that the displacement in the early stage is mainly caused by consolidation, while the displacement in the late stage is mainly caused by rheology.Effects of rheology on calculated results. (NR: rheological effect not considered; R: rheological effect considered).
(a)
(b)


## 5.1. Computed Model

Consolidation problems can be classified as small strain theory and finite strain theory according to the difference of geometric assumption. Small strain theory assumed that the strain is little, but the consolidation deformation is usually very large for deep soft clay foundation. Take dredger fill for example, the deformation may exceed 80%, for which the small strain theory is no longer suitable [14, 15]. As the pioneer researchers studied on finite strain consolidation theory, Mikasa [16] and Gibson et al. [17] derived their one-dimensional finite strain consolidation equations, respectively, in 1960s, and many scholars developed further research on the topic in the following decades. Because of strong nonlinearity, most of existing researches on finite strain consolidation adopted some simplified assumptions without considering rheological characteristics. By means of ABAQUS, this paper studied the rheological consolidation of soft clay considering finite strain effect.The computed model is plane strain problem as shown in Figure11. The soil foundation is homogeneous, anisotropic, saturated, and normally consolidated soft clay with a thickness of 10.0 m, and only the top boundary is permeable. In addition, the semilogarithmic empirical formula proposed by Tavenas et al. [18] was adopted to consider the variation of permeability in the process of consolidation as follows:
(31)
k
v
=
k
v
0
·
1
0
(
e
-
e
0
)
/
C
k
,
where k
v
0 is the initial permeability coefficients, and k
v is the permeability coefficient in the consolidation process; e
0 is the initial void ratio, and e is the void ratio in the consolidation process; C
k is the permeability index. According to the above work, the model parameters are summarized in Table 8 approximately.Table 8
Material parameters of computed model.
E
0 (MPa)
E
1 (MPa)
η
1 (MPa·s)
η
2 (GPa·s)
γ
sat
⁡ (kN·m−3)
k
v
0 (m·s−1)
P
0 (kPa)
M
λ
κ
e
0
μ
C
k
2.5
11.5
150
24000
17.7
4.5
E
-
9
50
0.898
0.102
0.015
1.08
0.35
0.5Figure 11
Calculation model and mesh generation.


## 5.2. Consolidation Property of Soft Clay

From Figure12 we can see the differences of the consolidation process between finite strain and small strain. Figure 12(a) shows that for any location the excess pore pressure of finite strain is less than that of small strain at the same time, and the difference is increasing with the increase of time in beginning period. Figure 12(b) shows that the displacement of finite strain is less than that of small strain and the difference increased significantly if the load is increasing. For example, when q = 200 kPa, the difference may take about 20% while it is nearly 0 when q = 50 kPa. So, it is essential to consider the error caused by small strain if the load is large and the compressibility of the soil is large.Results comparison of large strain and small strain (SS: small train; FS: finite strain).
(a)
(b)From Figure13 we can see the influence of rheology on the consolidation process. Figure 13(a) shows that the dissipation of excess pore pressure is significantly slower if the rheological effect is considered; for example, when t = 1000 d, the excess pore pressure is nearly 0 if the rheological effect is not considered, but it is still larger than 20 kPa if the rheological effect is considered. Figure 13(a) shows that the displacement is influenced by rheological effect too; for example, in the early stage (t ≤ 100 d), the displacement is nearly the same for two cases, but the final displacement may have a difference of more than 35%, which indicates that the displacement in the early stage is mainly caused by consolidation, while the displacement in the late stage is mainly caused by rheology.Effects of rheology on calculated results. (NR: rheological effect not considered; R: rheological effect considered).
(a)
(b)


## 6. Conclusions

(1)
An elastic-viscoplastic model developed to describe the time-dependent behavior of normally consolidated soft clay was presented in the paper on the framework of modified cam-clay elastic plastic model and Perzyna’s overstress viscoplastic theory.
(2)
A series of laboratory triaxial rheological tests of Ningbo soft clay with different confining pressure and deviatoric stress are performed; the results show that the presented constitutive model was suitable for the rheological characteristic of Ningbo soft clay, and it is achievable to determine the parameters of presented model by curve fitting.
(3)
The analysis on the model parameters shows that the value of all parameters is related to confining pressure and deviatoric stress. It may need further study on the relationship between parameters and the stress level so as to use the model conveniently.
(4)
In the proposed model, the viscoelastic viscosity coefficientη
1, the viscoplastic viscosity coefficient η
2, and the initial yield surface p
0 are the most important parameters. η
1 and η
2 have important influence on the time of viscous strain developed, and p
0 have important influence on both the viscoplastic strain rate and the final viscoplastic strain.
(5)
The existence of rheology will significantly slow down the dissipation of the excess pore pressure. Considering the rheological effect, the total displacement is larger, and the displacement in the early stage is mainly caused by consolidation, while in the late stage it is mainly caused by rheology.
(6)
For consolidation problems of high compressibility soft clay foundation, great error will be caused by small strain assumption, so finite strain effect should be considered. When finite strain effect is considered, the consolidation process is developed more rapidly, and the computed consolidation displacement will be less than that of small strain.



---

*Source: 587412-2014-06-05.xml*