# Modified Liu-Carter Compression Model for Natural Clays with Various Initial Water Contents

Authors: Sen Qian; Jian Shi; Jian-wen Ding
Journal: Advances in Civil Engineering
 (2016)

---

## Abstract

The initial water content has a significant effect on the compression behaviour of reconstituted clays. This effect has to be considered in the Liu-Carter model to ensure the addition voids ratio only related to soil structure. A modified Liu-Carter compression model is proposed by introducing the empirical equations for reconstituted clays at different initial water contents into the Liu-Carter model. The proposed model is verified against the experimental results from the literature. The simulations by the proposed method are also compared with that by old method where the influence of initial water content is not considered. The results show that the predicted virgin compression curves of natural clays are similar, but the values ofb and Δ
e
y may be very different.


---

## Body

## 1. Introduction

Natural clays generally sustain a higher void ratio than the corresponding reconstituted clays at the same stress level in the compression curves [1–4]. The void ratio difference (Δ
e) between natural clays and reconstituted clays at the same effective stress often decreases with increasing consolidation stress when the applied stress exceeds the consolidation yield stress [4–7]. This void ratio difference Δ
e, named as the additional voids ratio, is commonly attributed to the effect of soil structure that is developed during the depositional and the postdepositional processes [1–3]. The decrease of Δ
e with stress level is owing to the progressive destructuring [5, 6, 8].Liu and Carter [5, 6] proposed a compression model to quantify the effect of destructuring on the additional voids ratio Δ
e. They assumed that Δ
e is only controlled by soil structure and inversely proportional to the current effective stress during virgin compression. This model is termed as Liu-Carter compression model in this study. The formulation of this model is simple and the involved parameters have the clear physical meanings. As a consequence, it has been widely and successfully employed for constitutive modelling of structured clays [9–11].Since the compression line of reconstituted clay is not routinely measured [3], the reference line of determining Δ
e is generally estimated from the intrinsic compression line (ICL) [6, 8, 9]. Note that the determination of ICL is based on the oedometer test data on reconstituted clays at an initial water content equal to 1.25 times the liquid limits [3]. It has been well documented that the initial water contents of natural clays vary within a wide spectrum of liquid limits [8, 12, 13]. In addition, several researchers have illustrated that the initial water content significantly affects the compression behaviour of reconstituted clays [13–15]. Hong et al. [13] reported that for a given reconstituted clay the virgin compression line at a higher initial water content lies above that at a lower initial water content, as shown in Figure 1. Thus when the initial water content of the reconstituted clay is not close to the natural water content of the natural clay, the additional voids ratio Δ
e will depend on not only soil structure but also the difference of initial water content. To ensure that Δ
e is directly related to soil structure, the estimated reference line is therefore necessary to consider the influence of the initial water content.Figure 1
Typical compression curves of reconstituted clays at different initial water contents (test results from [13]).In this study, a modified Liu-Carter compression model is proposed by incorporating the empirical equations proposed by Zeng et al. [15], covering a wide range of initial water contents for reconstituted clays. This modified model can take into account the effect of initial state on the compression behaviour of reconstituted clays as a reference line. The proposed model is subsequently verified against 16 oedometer tests on 8 different natural clays available from the literature.


## 2. Modified Liu-Carter Compression Model

Taking the compression curve of reconstituted clay as the reference line, Liu and Carter [5, 6] suggested that the void ratio of a natural clay, e, can be expressed in terms of the void ratio of corresponding reconstituted clay e
∗ and the additional voids ratio Δ
e at the same stress level, as shown in Figure 2. Through the quantification of Δ
e at a given consolidation stress, the one-dimensional compression behaviour of natural clays can be evaluated as(1)e=e∗+Δeyσvy′σv′bforσv′≥σvy′,where e and e
∗ represent the void ratio of natural clay and reconstituted clay, respectively; σ
v
′ is the effective vertical stress; σ
v
y
′ is the effective vertical yield stress at which the natural clay begins to yield; Δ
e
y is the addition voids ratio at σ
v
y
′, owing to the original soil structure (Figure 2). b is destructuring index that is a constant parameter for quantifying the rate of destructuring. The more rapidly the destructuring develops, the higher the value of b is. Its value depends on soil type and soil structure [9].Figure 2
Idealization of the compression behaviour of natural clay and reconstituted clays at different initial water contents (after Bian et al. [16]).In the absence of the laboratory compression curve for corresponding reconstituted clay, ICL proposed by Burland [3] is often used as the reference line [6, 8, 9]. ICL is obtained by using the void index (I
v) to correlate the compression curves of reconstituted clays with initial water contents of about 1.25 times the liquid limits. Based on the definition of I
v, the virgin compression curve of reconstituted clays can be expressed as [3](2)e∗=e100∗+Cc∗×Iv,where e
100
∗ is the void ratio of reconstituted clays at the effective vertical stress of 100 kPa; C
c
∗ is termed the intrinsic compression index. In I
v-log
⁡
σ
v
′ space, ICL can be described by the following equation:(3)Iv=2.45-1.285log⁡σv′+0.015log⁡σv′3.Burland [3] also showed that components e
100
∗ and C
c
∗ in (2) can be estimated by the following two empirical equations:(4)e100∗=0.109+0.679eL-0.089eL2+0.016eL3Cc∗=0.256eL-0.04,where e
L is the void ratio at the liquid limit, and it can be evaluated directly using the equation e
L
=
G
s
×
w
L (G
s is the density of soil particle; w
L is the liquid limit). Therefore, when e
L ranges from 0.6 to 4.5 and the initial void ratio e
0
=
1.25
e
L, the virgin compression curve of reconstituted clays can be estimated by G
s and w
L.However, the initial water contents of natural clays are not always around 1.25 times the liquid limits [8, 12, 13], and the intrinsic compression curves or intrinsic compression parameters for reconstituted clays may be prominently affected by initial water contents (Figure 1) [13–15]. Hong et al. [13] illustrated that the virgin compression lines of various reconstituted clays with the initial water content at 0.7–2.0 times the liquid limits can also be normalized by the void index I
v. The normalized curve can be represented by an extended intrinsic compression line (EICL), which is expressed as [13](5)Iv=3.0-1.87log⁡σv′+0.179log⁡σv′2.Based on this intrinsic concept, new empirical correlations of e
100
∗ and C
c
∗ were proposed by Zeng et al. [15]:(6)e100∗=0.357+0.171e0+0.223eLCc∗=-0.069+0.109e0+0.152eL.It should be noted that (6) are based on the test data of various reconstituted clays with e
L in the range between 0.76 and 2.7 and e
0
/
e
L in the range between 0.7 and 2.2. Zeng et al. [15] also recommended that for the reconstituted clays within an extended range of e
L from 0.66 to 5.72 the intrinsic parameters can be estimated by following relationships with a possible reduction of accuracy:(7)e100∗=0.223+0.261e0+0.282eL-0.018e02-0.05eL2+0.015eL3Cc∗=-0.064+0.153e0+0.11eL-0.006e02.Knowing three parameters G
s, w
L, and e
0, the compression line of reconstituted clays could be simulated by combining (2) and (5). e
100
∗ and C
c
∗ can be evaluated by (6) or (7). Using this simulated compression curve as the new reference line in the Liu-Carter model, a modified Liu-Carter compression model is proposed:(8)e=e100∗+Cc∗×Iv+Δeyσvy′σv′bforσv′≥σvy′,where I
v, e
100
∗, and C
c
∗ are determined by corresponding equations (5)–(7). It is showed that this modified Liu-Carter model has six basic parameters G
s, w
L, e
0, σ
v
y
′, Δ
e
y, and b. The first three parameters are physical parameters, and the last three parameters can be obtained from the laboratory consolidation tests on natural clays. The virgin compression behaviour of natural clays with different initial states could be described using this modified Liu-Carter compression model.


## 3. Validation of Modified Liu-Carter Model

To validate the modified Liu-Carter compression model, the results of 16 oedometer tests on 8 types of natural clays with different initial water contents are compiled from the literature, as shown in Table1. Values of all parameters needed for simulation are also summarized in Table 1. The void ratio at the liquid limit e
L can be simply obtained from G
s and w
L. The value of σ
v
y
′ is identified directly from the measured compression curves using the method of Liu and Carter [6], as shown in Figure 2. The additional voids ratio at yield stress Δ
e
y is determined by comparing the measured compression curve with the predicted reference line suggested by Zeng et al. [15]. The destructuring index b is evaluated by best fitting equation (8) with the measured compression data of natural clays.Table 1
Details of soil samples and model parameters.
Soils
G
s (g⋅cm−3)
w
L (%)
e
0
e
0
/
e
L
σ
v
y
′ (kPa)
Δ
e
y
b
Reference
Singapore clay
2.74
84
1.59
0.69
218
0.50
0.54
Low et al. [8]
Gulf of Guinea clay
2.65
118
3.28
1.05
98
1.17
0.55
Hattab et al. [17]
Lianyungang clay
2.73
53
1.79
1.24
69
0.61
0.60
Deng et al. [18]
Jangyu clay
2.71
60
2.35
1.45
60
0.97
0.44
Chung et al. [19]
Undisturbed Ariake clay
2.66
127
4.30
1.27
38
1.15
0.30
Liu and Hong [20]
Disturbed Ariake clay
2.66
127
4.30
1.27
29
1.27
0.26
Liu and Hong [20]
Shanghai clay
2.70
45
1.15
0.94
160
0.31
0.24
Li et al. [21]
Batiscan clay
2.67
43
2.13
1.85
117
1.15
1.05
Leroueil et al. [22]
Wenzhou clay
Bian [23]
5 m
2.70
65
1.91
1.05
61
0.60
0.53
7.5 m
2.73
70
1.94
1.01
58
0.52
0.41
10 m
2.71
69
1.90
1.02
80
0.58
0.53
12.5 m
2.71
73
2.01
1.02
92
0.74
0.56
15 m
2.70
76
1.84
0.90
100
0.53
0.85
17.5 m
2.71
74
1.73
0.87
148
0.54
0.42
20 m
2.71
65
1.60
0.91
138
0.55
0.57
22.5 m
2.73
67
1.42
0.78
156
0.37
0.79As listed in Table1, the samples are investigated over a wide spectrum of liquid limits (43–127%) and initial state e
0
/
e
L (0.69–1.85). And it can be found that except Gulf of Guinea clay and Ariake clay, e
L value is within 0.76 to 2.7. According to the suggestion of Zeng et al. [15], e
100
∗ and C
c
∗ in (8) for Gulf of Guinea clay and Ariake clay should be estimated by using (7), while for other clays equations in (6) should be employed. The comparisons between the simulations and the measurements of the virgin compression curves are shown in Figures 3–7. The initial yield point of each sample is represented by a solid circle. The experimental data and simulated compression curves for natural clays are indicated by open dots and thick solid lines, respectively. The predicted reference lines are showed as broken lines. It should be pointed out that the compression behaviour of natural clays prior to yielding was not investigated here.Figure 3
Comparison between simulated and experimental results for Singapore clay and Gulf of Guinea clay based on test data reported by Low et al. [8] and Hattab et al. [17], respectively.Figure 4
Comparison between simulated and experimental results for Lianyungang clay and Jangyu clay based on test data reported by Deng et al. [18] and Chung et al. [19], respectively.Figure 5
Comparison between simulated and experimental results for undisturbed and disturbed Ariake clay based on test data reported by Liu and Hong [20].Figure 6
Comparison between simulated and experimental results for Shanghai clay and Batiscan clay based on test data reported by Li et al. [21] and Leroueil et al. [22], respectively.Figure 7
Comparison between simulated and experimental results for Wenzhou clay based on test data reported by Bian [23]: (a) 5 m and 12.5 m; (b) 7.5 m and 15 m; (c) 10 m and 17.5 m; (d) 20 m and 22.5 m.
(a)
(b)
(c)
(d)Figures3–7 show that all the simulated virgin compression curves agree well with the experimental data of the natural clays at different initial states. It demonstrates that the modified Liu-Carter compression model is capable of quantifying the influence of soil structure on the virgin compression of natural clays.Furthermore, for the purpose of model application, some discussions on the key parametersΔ
e
y and b are necessary. The addition voids ratio Δ
e
y owing to the influence of original soil structure should be closely correlated with the stress sensitivity S
σ
y (=
σ
v
y
′
/
σ
y
e
∗), where σ
y
e
∗ is the equivalent stress on the intrinsic compression line at the void ratio of natural clay under σ
v
y
′ [8]. Low et al. [8] suggested that if an experimental compression curve is unavailable, S
σ
0 (=
σ
v
y
′
/
σ
v
0
e
∗) instead of S
σ
y can be used to estimate the value of Δ
e
y, where σ
v
0
e
∗ is the equivalent stress on the intrinsic compression line at the initial void ratio. But the ratio between Δ
e
y and ln
⁡
S
σ
0 is not always a constant for different clays [8] (Figure 8). Similarly, a different parameter Δ
I
v
y (=
Δ
e
y
/
C
c
∗) could be used here, which is the additional void index at σ
v
y
′. Figure 9 shows the relationship between Δ
I
v
y and ln
⁡
S
σ
0 for the soils in Table 1, and regression equation (9) is obtained with R
2
=
0.80. Hence, when the experimental data is unavailable, the parameter Δ
e
y could be estimated via (10):(9)ΔIvy=0.52ln⁡Sσ0(10)Δey=0.52Cc∗ln⁡σvy′σv0e∗.Currently, it is not possible to perform a systematic and detailed study on the prediction of b due to the limited amount of experimental data [6, 8]. But some changes in the b value can be generalized. The b value for these natural clays is between 0.24 and 1.05, which is comparable with that reported by Liu and Carter [6] for natural soft clays (ranging from 0.38 to 1.32). As seen in Table 1, the b value of undisturbed Ariake clay is higher than that of disturbed Ariake clay (Figure 5). This implies that the rate of destructuring can be reduced by disturbance. Liu and Carter [6] hypothesized that for clay samples of a given mineralogy and with similar geological stress history but at different depths, b may be mainly dependent on the liquidity index. Similarly, for natural Wenzhou clay at different depths taken from the same location, b tends to decrease with the increase of the value of e
0
/
e
L, as shown in Figure 10.Figure 8
The relationship betweenΔ
e
y and ln
⁡
S
σ
0 for the clays in Table 1.Figure 9
The relationship betweenΔ
I
v
y and ln
⁡
S
σ
0 for the clays in Table 1.Figure 10
The values ofb versus e
0
/
e
L for Wenzhou clay at different depths.


## 4. Discussion on the Effect of the Initial State

In order to investigate the influence of initial state on the simulations, the predicted curves and parametersb and Δ
e
y in the modified Liu-Carter model are compared with the ones determined by (8) but with I
v, e
100
∗, and C
c
∗ evaluated via (3)-(4), where the influence of initial water contents is not considered. Singapore clay (e
0
/
e
L
=
0.69), Lianyungang clay (e
0
/
e
L
=
1.24), and Batiscan clay (e
0
/
e
L
=
1.85) are chosen to be analyzed in this part because that they cover the whole range of initial water contents listed in Table 1.The comparison results of predicted compression curves are summarized in Figure11. The predicted compression curves of natural clays and reconstituted clays without considering the initial state are represented by thin solid lines and thin broken lines, respectively. The corresponding ones for the modified Liu-Carter model are represented by thick solid lines and thick broken lines. It can be seen from Figure 11 that the predicted curves of natural clays without considering the initial state are almost identical to that considering the initial state, and they both agree well with the experimental data. It means that the initial states may have no noticeable effect on the predicted compression curves of natural clays. But as seen from Figures 11(a) and 11(c), the predicted compression lines of reconstituted Singapore clay and reconstituted Batiscan clay (e
0
/
e
L is not around 1.25) differ greatly in locations. It is expected that the different locations of the reference line would cause great difference in the value of b and Δ
e
y.Figure 11
The influence of initial state on simulated curves: (a) Singapore clay,e
0
/
e
L
=
0.69; (b) Lianyungang clay, e
0
/
e
L
=
1.24; (c) Batiscan clay, e
0
/
e
L
=
1.85.
(a)
(b)
(c)Figures12(a) and 12(b) show the comparison results of the key parameters b and Δ
e
y. For Lianyungang clay with e
0 equal to 1.24
e
L, the values of b and Δ
e
y obtained from two different methods are almost the same with a little difference less than 5 percent. The reason is that the reference lines estimated nearly coincide (Figure 11(b)). However, for Singapore clay with e
0 equal to 0.69
e
L, compared to the corresponding values for the modified Liu-Carter model, b and Δ
e
y without considering the initial state are 67 percent higher and 48 percent lower, respectively. For Batiscan clay with e
0 equal to 1.85
e
L, b and Δ
e
y without considering the initial state are 11 percent lower and 14 percent higher than the ones for the modified Liu-Carter model. This observation indicates that when e
0
/
e
L is not around 1.25, using the estimated reference line without considering the initial state may lead to a noticeable error in the determination of b and Δ
e
y. In other words, it is essential for parameters b and Δ
e
y of different natural clays to consider the influence of initial state. Thus, the quantification of destructuring by the modified Liu-Carter model can be successful for various natural clays at a wide range of initial water content.Figure 12
The influence of initial state on the model parameters: (a)b versus e
0
/
e
L; (b) Δ
e
y versus e
0
/
e
L.
(a)
(b)


## 5. Conclusions

Based on the modelling, validation, and discussion presented in this study, the following conclusions can be drawn:(1)
By introducing the empirical equations for reconstituted clays proposed by Zeng et al. [15] to Liu-Carter compression model, a modified Liu-Carter model is proposed to consider the effect of initial state. The mathematical formulation is e
=
e
100
∗
+
C
c
∗
×
I
v
+
Δ
e
y
σ
v
y
′
/
σ
v
′
b, for σ
v
′
≥
σ
v
y
′, where I
v
=
3.0
-
1.87
log
⁡
σ
v
′
+
0.179
(
log
⁡
σ
v
′
)
2, and e
100
∗ and C
c
∗ are determined by three physical parameters G
s, w
L, and e
0.
(2)
The numerous comparisons of simulations and experimental results show that the modified Liu-Carter model could quantify successfully the virgin compression behaviour of natural clays over a wide spectrum of initial water content.
(3)
Compared to the corresponding values ofb and Δ
e
y for the modified Liu-Carter model, the values of b and Δ
e
y without considering the influence of initial water contents will be higher and lower, respectively, when e
0
/
e
L is less than 1.25. On the contrary, the values of b and Δ
e
y will be lower and higher, respectively, when e
0
/
e
L is higher than 1.25.



---

*Source: 1691605-2016-06-20.xml*