Advances in Civil Engineering
Volume 2019 (2019), Article ID 6042893, 14 pages
https://doi.org/10.1155/2019/6042893

Research Article

Study on the Creep Behaviors of Interactive Marine-Terrestrial Deposit Soils

Zhen Yan,¹ Gang Li,² Jinli Zhang,³ and Rui Zhang²

¹Tianjin Research Institute for Water Transport Engineering, Ministry of Transport, Tianjin 300456, China
²Shaanxi Key Laboratory of Safety and Durability of Concrete Structures, Xijing University, Xi’an, Shaanxi 710123, China
³State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian, Liaoning 116024, China

Correspondence should be addressed to Gang Li; t_bag945@126.com

Received 20 November 2018; Accepted 18 February 2019; Published 7 July 2019

Academic Editor: Giulio Dondi

Copyright © 2019 Zhen Yan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The interactive marine-terrestrial (IMT) deposit soils were formed in the complex depositional environment; their mechanical properties are different from the other deposits. The creep behaviors of Dalian clayey soils were studied according to one-dimensional creep tests and drained triaxial creep tests. Based on the creep test results, the empirical model was established to describe the one-dimensional creep behavior and triaxial creep behavior, respectively. The results showed that Dalian deposits have typical nonlinear creep behavior. With the increasing of consolidation pressure, the strain is increased, the stability time is extended, and the demarcation point between primary and secondary consolidation is more obvious. The deposits belong to medium to high secondary compressibility soil, and the secondary consolidation coefficient is decreased with the increasing of consolidation time and increased with consolidation pressure increasing. The ratio between secondary consolidation coefficient and compression index at different depths changes from 0.033 to 0.058, which conform to Mesri conclusion. Under low deviator stress, the creep processes showed the characteristic of attenuation creep and shear contraction. However, it showed the characteristic of acceleration creep, shear contraction, and shear dilatancy under damage deviator stress. The axial strain rate decreased with the increasing of creep time and increased with the deviator stress increasing, while the deviator stress has little effect on the m values. The tests results agree well with the calculation results, which showed that the creep equation is suitable for describing the creep behaviors of Dalian interactive marine-terrestrial deposits.

1. Introduction

The engineering properties of soils are related to their sedimentary environments. Interactive marine-terrestrial (IMT) deposits were formed from historical processes including transgression-regression and the reciprocating movement of coastlines [1]. Compared with typical marine deposits, IMT deposits generally have lower natural water content and void ratios but relatively better mechanical properties (compression index, friction angle, etc.) and more variable physical and mechanical properties. Chen [2] studied the consolidation properties of IMT soft deposits in the Pearl River Delta, analyzed the effects of consolidation on their basic mechanical properties, and proposed a quantitative analysis method that considers the consolidation effect. By examining the physical and mechanical properties of an IMT clay at the Yangtze River estuary, Liu and Shi [3] found that the physical property indices of soil in this region were not consistent with its mechanical property indices, analyzed the reason behind this phenomenon, and proposed a new method for evaluating the engineering properties of this soil type. Zhang et al. [4] analyzed the physical property indices of an IMT clay on the north side of the Yangtze River estuary using principal component analysis and noted unreasonableness in using the liquidity index to determine its consistency. Liu et al. [5] studied the pore pressure response of a PHC (high strength, prestressed concrete pipe piles) composite foundation in the IMT deposits on the north shore of the Yangtze River estuary and found that shallow deposits had relatively high permeability and sensitive excess pore pressure response but that deep deposits had low permeability and excess pore pressure that increased with the overlying load. Tuo et al. [6] investigated the engineering properties of Zhanjiang IMT soft clays and noted that under the influence of tectonic movement, these deposits exhibited special sedimentary characteristics due to the coupling action of complex geological stresses and had relatively high structural strength due to aggregation. Shi et al. [7] examined the mechanical properties of natural intermediate deposits on the north shore of the Yangtze River estuary and found that the physical properties of this type of soil were similar to those of soft clays, but the in situ strength and natural water content were similar to those of hard clays. They concluded that particle distribution has a relatively significant impact on the in situ strength of natural deposits.

Soil distribution exhibits typical regional characteristics. To date, relatively few studies have been conducted to examine the creep behavior of IMT deposits. In this study, the creep behavior of Dalian IMT clays was examined by one-dimensional (1D) and drained triaxial creep tests. Based on the creep test results, a creep equation that considers the primary and secondary consolidation of deposits and an improved Mesri creep equation were established, providing a reference for predicting foundation deformation and settlement.

2. Materials and Methods

2.1. Soil Samples

Test samples were collected from an undisturbed seabed at a depth of 5–30 m at a construction site. These samples were gray-grayish brown IMT and silty clays. Table 1 summarizes their physical property indices. As shown in Table 1, the marine deposits (mud and muddy clays) had relatively high water contents and void ratios and relatively small unit weights and were in a liquid-plastic state. In contrast, the IMT deposits (clays and silty clays) had relatively low water contents and void ratios and relatively large unit weights and were in a plastic state.

Table 1: Physical properties indices of clays.

Genetic factor
Soil
W(%)
γ(kN/m³)
e
(%)
IP
IL

Marine deposit
Mud
63.7
16.0
1.75
48.5
24.3
1.63

Muddy clay
47.8
17.2
1.30
40.2
19.9
1.38

Marine-terrestrial deposit
Clay
39.2
17.9
1.10
48.7
24.4
0.61

Silty clay
25.1
19.6
0.70
31.2
14.3
0.57

2.2. Test Methods

The soil samples collected at depths of 14 m (clay), 18 m (clay), and 25 m (silty clay) were subjected to a 1D creep test under drained conditions on a consolidometer (WG-1B). Repeated tests were performed for each group of samples. Each sample had a diameter of 61.8 mm and a height of 20 mm. Considering the in situ stress of the soil and the overlying backfill load, the following test load levels were applied: 50, 100, 200, 400, 800, and 1,600 kPa. After each load application, the amount of deformation was read at preset intervals until the deformation below 0.01 mm within 24 h. The temperature was controlled at 20 ± 1°C during the test process.

Drained triaxial creep tests were performed on a soil triaxial rheometer (CSS-2901TS). Each sample had a diameter of 39.1 mm and a height of 80 mm. A multistage (either four or five stages) loading procedure was adopted. Four levels of cell pressure were used in the tests, namely, 100, 200, 400, and 800 kPa. When the soil sample is failed in shear in the consolidated-drained shear test under the same cell pressure, the deviatoric stress (qf) was recorded and used to determine the load increment (qf/4 – qf/5) for each soil sample at each stage. Sample disturbance was minimized during the tests. All samples were saturated by standard increments of cell pressure and back pressure. The pore water pressure coefficient B-value was measured to confirm saturation until it reached 0.98. Each sample was successively saturated, consolidated, and subjected to multistage loading. Deformation below 0.01 within 24 h was set as the standard. The temperature was controlled at 20 ± 1°C during the test process.

3. Test Results and Analysis from One-Dimensional Creep Test

3.1. Creep Behavior

Figure 1 shows the time-history curves of strain and void ratio of the soil samples collected at depths of 14 m, 18 m, and 25 m. As shown in Figure 1, the soil strain was relatively small during the initial loading stage but gradually increased and eventually changed slowly as the consolidation time increased. The soil sample deformation was relatively large under high levels of vertical pressure (400 and 800 kPa). The longer the settling time, the more pronounced the boundary between primary and secondary consolidation. According to the laboratory experiments, the in situ void ratio e₀, in situ effective stress , preconsolidation pressure , and overconsolidation ratio (OCR) at depths of 14 m, 18 m, and 25 m are listed in Table 2. It can be concluded that the clays at depths of 14 m, 18 m, and 25 m belong to slightly overconsolidated soil. For a normal consolidated or slightly overconsolidated soil, when the load level was lower than the preconsolidation pressure, its deformation was relatively small; when the load level was higher than the preconsolidation pressure, its deformation and its rate of deformation increased, resulting in the aforementioned creep behavior. Contrary to the pattern of change in the strain, the void ratio decreased as the consolidation time increased. The slopes of the curves changed relatively significantly during the initial loading stage and remained basically unchanged during the later loading stage. The void ratio-time curves under low levels of consolidation pressure (100 and 200 kPa) were nearly straight lines, whereas those under high levels of consolidation pressure exhibited an inverse “S” shape. This is similar to the pattern of deformation of Hong Kong deposits reported by Yin [8]. Overall, the patterns of deformation of the soil samples collected at various depths (14, 18, and 25 m) were basically consistent. The results obtained using the log time method show that primary consolidation was basically completed within 100 min after each level of consolidation pressure was applied.

(a)

(b)

(c)

Figure 1: εa or e versus t at different depths. (a) 14 m, (b) 18 m, and (c) 25 m.

Table 2: Consolidation state of clays at different depths.

Depth (m)
e₀
(kPa)
(kPa)
OCR

14
1.12
100
195
1.95

18
1.10
135
210
1.56

25
0.71
200
220
1.10

To analyze the pattern of change in the strain with consolidation pressure, the creep test results for the soil samples collected at a depth of 25 m are used as an example. Figure 2 shows their isochronous stress-strain curves which depict notable nonlinear characteristics at different times. As the consolidation pressure increased, the slopes of the curves were more sharply increased. When the consolidation pressure remained constant, the strain increased with time. When the strain remained constant, the consolidation pressure decreased with time. Thus, the Dalian IMT deposits displayed typical nonlinear creep behavior which was related to the loading time. Under low consolidation pressures (50 and 100 kPa), the slopes of the isochronous curves were relatively small; under high consolidation pressures (800 and 1,600 kPa), the slopes of the isochronous curves were relatively large.

Figure 2: Isochronous curve of εa versus σ of silty clay at 25 m depth.

3.2. Secondary Consolidation Coefficient

Buisman noted that the secondary consolidation coefficient for the secondary consolidation and deformation stage could be calculated using the following equation:

where is the secondary consolidation coefficient, is the time at which 100% primary consolidation is achieved, is the calculated time for secondary consolidation, and and are the void ratios corresponding to and , respectively.

Figure 3 shows the curves of changes in the secondary consolidation coefficient of the soil samples collected at various depths (14, 18, and 25 m) with time (t₁ = 100 min), which were used to analyze the pattern of changes in the secondary consolidation coefficient with time. As shown in Figure 3, during the secondary consolidation and deformation stage, the secondary consolidation coefficient gradually decreased as time increased. Under low consolidation pressures, the slopes of the curves were relatively small, and the curves were nearly straight lines. As the consolidation pressure increased, the slopes of the curves gradually increased and the corresponding secondary consolidation coefficients also increased significantly. Overall, the soil samples collected at a depth of 18 m had a relatively large secondary consolidation coefficient, whereas the soil samples collected at a depth of 25 m had a relatively small secondary consolidation coefficient. Based on the classification of the secondary compression coefficient in Reference [9], the deposits at a depth of 18 m are soils with high secondary compressibility and those at depths of 14 and 25 m are soils with medium secondary compressibility. Thus, the Dalian IMT deposits are soils with medium to high secondary compressibility.

(a)

(b)

(c)

Figure 3: Cα versus t at different depths. (a) 14 m, (b) 18 m, and (c) 25 m.

3.3. Relationships between Secondary Consolidation Coefficient and Compression Index

Mesri and Godlewski [10] found that the ratio of the secondary consolidation coefficient to the compression index is in the range of 0.01 to 0.07 for undisturbed soils and is a constant for the same type of undisturbed soil. This conclusion has been accepted by a number of researchers [11–15]. Through calculation, it was found that the ratios of the secondary consolidation coefficient to the compression index for the Dalian IMT deposit samples collected from depths of 14, 18, and 25 m ranged from 0.033–0.058, consistent with Mesri’s conclusion.

4. One-Dimensional Creep Equations

4.1. Deriving the Creep Equations

The creep test results for the silty clay samples collected at a depth of 25 m are used as an example. Figure 4 shows the log-log plots of strain versus time based on the test results. Table 3 summarizes the fitted equations. As shown in Figure 4, there was a nearly linear strain-time relationship under each level of consolidation pressure, which satisfied the following relationship:

Figure 4: Double logarithmic curves of εa versus t.

Table 3: Fitting equation between stress and time.

Pressure (kPa)
Equation
R

50
lg (ε) = 0.11268 lg (t) − 0.48949
0.89886

100
lg (ε) = 0.07091 lg (t) + 0.37196
0.98067

200
lg (ε) = 0.05368 lg (t) + 0.53082
0.97618

400
lg (ε) = 0.04435 lg (t) + 0.74367
0.96874

800
lg (ε) = 0.03996 lg (t) + 0.88527
0.96786

1600
lg (ε) = 0.03124 lg (t) + 1.02855
0.96529

By rearranging equation (2), we have the following equation:

where is the consolidation stress, is the consolidation time, is the reference time, is the slope of the log strain versus log time, and is the strain corresponding to and .

Figure 5 shows the isochronous log-log plots of strain versus consolidation pressure for the silty clay samples collected at a depth of 25 m, showing that the strain increased as the consolidation pressure increased, with a nearly linear strain-consolidation pressure relationship at various times. The average slope of the isochronous curves was 0.63338. The strain and the consolidation stress satisfied the following relationship:

Figure 5: Double logarithmic curves of εa versus σ.

By rearranging equation (4), we have the following equation:

where is the reference stress and is the slope of the strain-consolidation stress curve.

By combining equations (3) and (5), the following creep equation is obtained:

where is the strain corresponding to the reference stress and the reference time and is the slope of the log strain versus log time.

Based on the previously mentioned analysis, the creep equation contains three parameters, namely, , m, and n. Of these, is the strain corresponding to the reference time and the reference stress , m is the slope of the log strain versus log time, and n is the slope of the strain-stress curve. All three parameters can be determined based on the creep test results. Table 4 summarizes the model parameters for various levels of consolidation stress obtained by creep test results regression analysis for the soil samples collected at depths of 14, 18, and 25 m.

Table 4: Fitting model parameters of clays at different depths.

Depth (m)
Pressure (kPa)
Reference pressure (kPa)
Reference time (min)
(%)
m
n

14
50
100
100
1.370
0.12573
0.75518

100
100
100
1.370
0.03569
0.75518

200
100
100
1.370
0.03242
0.75518

400
100
100
1.370
0.02535
0.75518

800
100
100
1.370
0.03024
0.75518

1600
100
100
1.370
0.02057
0.75518

18
50
100
100
0.955
0.11295
1.06804

100
100
100
0.955
0.09645
1.06804

200
100
100
0.955
0.08710
1.06804

400
100
100
0.955
0.07130
1.06804

800
100
100
0.955
0.04549
1.06804

1600
100
100
0.955
0.02077
1.06804

25
50
100
100
2.145
0.11268
0.63338

100
100
100
2.145
0.07091
0.63338

200
100
100
2.145
0.05368
0.63338

400
100
100
2.145
0.04435
0.63338

800
100
100
2.145
0.03996
0.63338

1600
100
100
2.145
0.03124
0.63338

4.2. Verifying the Creep Equations

To examine the equation’s reliability, calculations were performed using the creep equation based on the creep test results for soil samples collected at depths of 14, 18, and 25 m. Figure 6 shows the results which indicate a relatively large difference between the test and calculated values for the initial stage of consolidation pressure. As the consolidation time increased, the consistency between the test and calculated values gradually increased. Overall, there was a high degree of consistency between the calculated and test results. This suggests that this model is suitable for describing the creep properties of the Dalian IMT deposits. In addition, the creep equation features only a few yet highly accessible parameters and can be used to provide a reference for predicting long-term foundation deformation.

(a)

(b)

(c)

Figure 6: Test results compared with calculation results at different depths. (a) 14 m, (b) 18 m, and (c) 25 m.

5. Test Results and Analysis from Drained Triaxial Creep Test

5.1. Axial Strain and Creep Time

To analyze the drained creep properties of the IMT clays, the test data were processed using Chen’s loading method according to Reference [16]. Figure 7 shows the creep curves of the silty clay collected at depths of 26 m under cell pressures of 100, 200, 400, and 800 kPa. The creep curves exhibited notable nonlinear characteristics. During the initial stage of the deviatoric stress, the axial strain changed relatively significantly. As time increased, the strain slowly changed. The soil displayed notable attenuation creep behavior under the action of a low deviatoric stress and pronounced acceleration creep behavior under the action of the failure deviatoric stress. Overall, the failure deviatoric stress increased as the cell pressure increased. The patterns of soil deformation under the action of various cell pressures and deviatoric stresses remained basically consistent. The soil deformation basically changed slowly within 5,000 min under the action of the deviatoric stress. Based on the previously mentioned analysis, the Dalian IMT clays exhibited typical nonlinear creep behavior. The axial strains increased as the deviatoric stress increased and slowly changed with time. In addition, the IMT clays displayed attenuation and acceleration creep behavior during the entire creep process.

(a)

(b)

(c)

(d)

Figure 7: εa versus t under different cell pressures. (a) 100 kPa, (b) 200 kPa, (c) 400 kPa, and (d) 800 kPa.

5.2. Volumetric Strain and Creep Time

Figure 8 shows the volumetric strain versus time curves under cell pressures of 100, 200, 400, and 800 kPa. The volumetric strain gradually decreased as the creep time and deviatoric stress increased (except the action of the failure deviatoric stress), suggesting that the soil showed shear contraction behavior. Under the action of the failure deviatoric stress, the volumetric strain first decreased and then increased, suggesting that the soil displayed shear contraction and shear dilation behavior, which differs from the pattern of changes in the volumetric deformation of the Tianjin Binhai soft clays [17], which exhibit only shear contraction behavior.

(a)

(b)

(c)

(d)

Figure 8: εv versus t under different cell pressures. (a) 100 kPa, (b) 200 kPa, (c) 400 kPa, and (d) 800 kPa.

5.3. Deviatoric Stress and Axial Strain

To analyze the stress-strain relationship of the clay samples, the creep test results obtained under a cell pressure of 100 kPa for the silty clay collected at a depth of 26 m are used as an example. Figure 9 shows their deviatoric stress-strain isochronous curves which highlighting a notable nonlinear stress-strain relationship. As the deviatoric stress increased, the axial strain gradually increased. Under the action of a low deviatoric stress, the slopes of the isochronous curves were relatively small. As the deviatoric stress increased, the slopes of the curves gradually increased, and the corresponding axial strain also increased. Similar to the creep behavior of soft clays reported by Liu et al. [18], the isochronous curves at various times exhibited a basically consistent distribution. When the deviatoric stress remained constant, the axial strain increased with time. Based on the aforementioned analysis, the Dalian IMT clays exhibited typical nonlinear creep behavior, which became increasingly pronounced as the deviatoric stress increased.

Figure 9: Isochronal curves of εa versus q under 100 kPa cell pressure.

5.4. Axial Strain Rates and Creep Time

Figure 10 shows the curves of changes in the strain rate of the silty clay collected at a depth of 26 m with time. It illustrated that the strain rate decreased with time on the log-log plots. Under the action of the deviatoric stress, the slopes of the curves were relatively small during the initial stage (within 100 min). After the creep time exceeded 100 min, the slopes of the curves remained basically constant. This is similar to the creep behavior of soft clays reported by Lai et al. [19] and Wang and Yin [20]. In addition, the strain rate of the soil increased as the deviatoric stress increased. During the later stage of loading (after 100 min), the slopes of the curves under the action of various levels of deviatoric stress (except for the failure deviatoric stress) were basically the same and the m-value (slope) ranged from 0.77 to 0.82, suggesting that the deviatoric stress had a relatively insignificant impact on the m-value. In terms of the creep process, under the action of a low deviatoric stress (35, 75, and 120 kPa), the creep rate of the soil decreased with time, and the soil displayed attenuation creep behavior. Under the action of the failure deviatoric stress, the creep rate of the soil first decreased and then increased and the soil displayed attenuation and acceleration creep behavior. This corroborates the conclusion derived from the previous analysis.

Figure 10: Double logarithmic curves of versus t under 100 kPa cell pressure.

5.5. Axial Strain Rates and Deviatoric Stress

Figure 11 shows the axial strain rate versus deviatoric stress isochronous curves. The creep behavior of the Dalian IMT clays was similar to that of soft clays reported by Yan et al. [21]. As the deviatoric stress increased, the axial strain rate gradually increased. As the creep time increased, the axial strain rate gradually decreased. Under the action of a low deviatoric stress, the slopes of the curves were relatively small. As the deviatoric stress increased, the slopes of the curves increased, suggesting that the strain rate increased as the deviatoric stress increased. When the deviatoric stress remained constant, the strain rate decreased with creep time. When the strain rate remained constant, the deviatoric stress increased with creep time. Thus, for the Dalian IMT clays, there was a nonlinear relationship between axial strain rate and deviatoric stress, which became increasingly insignificant as creep time increased.

Figure 11: versus q under 100 kPa cell pressure.

6. Improved Mesri Creep Equations

6.1. Mesri Creep Equation

There are a relatively large number of empirical models for the creep behavior of clays, represented by the Singh–Mitchell and Mesri models [22–29]. The Singh–Mitchell model is only suitable for describing the stress-strain-time relationship within 20–80% of the deviatoric stress level. Therefore, the Mesri model is more advantageous. In this study, an attempt is made to use the Mesri creep equation to describe the stress-strain-time relationship of the Dalian IMT clays based on the triaxial creep test results. The Mesri creep equation is as follows:

where is the axial strain, is the creep time, is the initial tangent modulus, is the shear strength, is the shear stress level , is the failure stress ratio, is the reference time, and is the test constant.

When t = t₁, equation (7) can be expressed as follows:

where and are the intercept and slope of the curve, respectively, and is the slope of the curve.

The Mesri creep equation contains three parameters, namely, , Rf, and λ, all of which can be determined by the regression analysis of test result.

6.2. Establishing the Model

The creep test results obtained for the silty clay collected at a depth of 26 m and 30 m under cell pressures of 100 kPa and 200 kPa are used as examples to demonstrate how the Mesri model parameters are determined. Figure 12 shows the fitting results (reference time: 10 h) that indicate a basically linear relationship between and under cell pressures of 100 and 200 kPa, which with a correlation coefficient greater than 0.98. In addition, there was also a linear relationship between and in the log-log plot, which with a correlation coefficient greater than 0.88. Table 5 summarizes the Mesri model parameters for the silty clay samples collected at depths of 26 and 30 m under cell pressures of 100 and 200 kPa, which determined by regression analysis.

(a)

(b)

Figure 12: Fitting parameters of Mesri model under different cell pressures. (a) 100 kPa and (b) 200 kPa.

Table 5: Mesri model parameters.

Depth (m)
σ₃ (kPa)
(σ₁ − σ₃)f (kPa)
t₁ (h)

(R)f
λ

26
100
150
10
0.2333
2.7558
0.4508
0.3173

0.5000
2.7558
0.4508
0.1093

0.8000
2.7558
0.4508
0.0750

200
265
10
0.2453
2.3540
0.6136
0.3109

0.5208
2.3540
0.6136
0.1066

0.7547
2.3540
0.6136
0.0856

30
100
183
10
0.2732
4.9351
0.2289
0.1825

0.5464
4.9351
0.2289
0.0958

0.8470
4.9351
0.2289
0.0646

200
345
10
0.2667
2.8648
0.4555
0.1425

0.5420
2.8648
0.4555
0.0997

0.8116
2.8648
0.4555
0.0611

It is worth noting that the fitting results for the failure deviatoric stress are not shown in Figure 12. In view of this, the creep process under the action of the failure deviatoric stress was subjected to fitting analysis using the method of least squares. A relatively high goodness of fit was achieved when using a cubic polynomial. Therefore, based on the concept of piecewise functions, the Mesri equation was coupled with a cubic polynomial to describe the clay creep process. The Mesri equation was used for simulations in which the deviatoric stress was less than the failure deviatoric stress. A cubic polynomial was used for simulations in which the deviatoric stress was equal to or greater than the failure deviatoric stress. The model is represented by the following equation:

Table 6 summarizes the model parameters determined by fitting the creep test results of the silty clay at depths of 26 and 30 m, using a cubic polynomial.

Table 6: Polynomial model parameters.

Depth (m)
σ₃ (kPa)

a
b
c
d

26
100
1
0.0054
−0.0882
0.5448
4.1375

200
1
0.0442
−0.3451
1.0045
3.6136

30
100
1
0.0570
−0.2180
0.3130
5.4140

200
1
0.0730
−0.4540
1.3100
3.8030

6.3. Verifying the Model

The newly established model was used to calculate the creep process of the silty clay collected at depths of 26 and 30 m under cell pressures of 100 and 200 kPa. Figure 13 shows the comparison which indicates that the calculated results obtained using the Mesri model under cell pressures of 100 and 200 kPa matched the test results relatively well. The calculation results obtained using the polynomial can relatively well describe the accelerated creep process. The relative difference between the calculation and test results was small. The high degree of consistency between the test and calculation results validates the reliability of the newly established model and indicates its suitability for describing the creep behavior of Dalian IMT clay.

(a)

(b)

(c)

(d)

Figure 13: Test results compared with calculation results at different depths and cell pressures. (a) 26 m depth and 100 kPa, (b) 26 m depth and 200 kPa, (c) 30 m depth and 100 kPa, and (d) 30 m depth and 200 kPa.

7. Conclusions

The creep behavior of Dalian IMT clay was studied using 1D and triaxial drained creep tests. Based on the creep test results, a 1D creep equation in the form of a power function and a three-dimensional (3D) creep equation consisting of the Mesri equation and a cubic polynomial were established. The conclusions are summarized as follows:

(1)As the consolidation time and pressure increased, the strain of the Dalian IMT deposits increased, whereas the void ratio decreased. The clays basically entered the stage of secondary consolidation and deformation after the action of consolidation pressure at all levels over 100 min.(2)The Dalian IMT deposits exhibited typical nonlinear creep behavior. Their creep behavior became increasingly pronounced as the consolidation pressure increased. The secondary consolidation coefficient of the clays was related to the consolidation pressure. The ratio of the secondary consolidation coefficient to the compression index was consistent with the conclusion derived by Mesri.(3)The Dalian IMT clays exhibited typical nonlinear creep behavior. The clays displayed shear contraction behavior under the action of a low deviatoric stress and shear contraction and shear dilation behavior under the action of the failure deviatoric stress.(4)The axial strain rate decreased with creep time and increased with increasing deviatoric stress. The deviatoric stress had a relatively insignificant impact on the m-value. Under the action of a low deviatoric stress, the creep rate gradually decreased with creep time, and the clays displayed attenuation creep behavior. Under the action of the failure deviatoric stress, the creep rate first decreased and then increased with creep time, and the clays exhibited attenuation and acceleration creep behavior.(5)There was a relatively high degree of consistency between the test results and the calculation results, which validating the newly established model. This suggests that both the 1D and 3D creep equations are suitable for describing the creep behavior of Dalian IMT clay and can provide guidance for predicting long-term foundation deformation.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was supported by the Special Fund for Scientific Research by Xijing University (XJ18T01), Special Fund for Scientific Research by Shaanxi Provincial Education Department (18JK1199), and the Fundamental Research Funds for the Central Research Institutes (TKS170102 and TKS170108).

References

Z. B. Liu, S. Y. Liu, B. M. Zhou et al., “Geological genesis and spatial variability characteristic of the interactive marine and terrestrial soft deposit in Yangtze River estuary,” Chinese Journal of Underground Space and Engineering, vol. 10, no. 1, pp. 102–108, 2014.

X. P. Chen, “Consolidation effect of soft soil in interactive marine and terrestrial deposit,” Chinese Journal of Geotechnical Engineering, vol. 33, no. 4, pp. 520–528, 2011.

Y. H. Liu and M. L. Shi, “Evaluation research and distinguishing of cohesive soil properties in the interactive marine & terrestrial deposit,” Rock and Soil Mechanics, vol. 29, no. 2, pp. 523–528, 2008.

H. Zhang, Y. H. Liu, M. L. Shi et al., “Research of characteristics of the sedimentary clay in the alternate marine and terrigenous environment,” Rock and Soil Mechanics, vol. 26, no. 9, pp. 1491–1494, 2005.

Z. B. Liu, S. Y. Liu, F. Jing et al., “Pore pressure response of composite foundation in interactive marine and terrestrial deposit of Yangtze River estuary,” Chinese Journal of Geotechnical Engineering, vol. 37, no. S2, pp. 91–95, 2015.

Y. F. Tuo, L. W. Kong, A. G. Guo et al., “Occurrence and engineering properties of structural soft clay in Zhanjiang area,” Rock and Soil Mechanics, vol. 25, no. 12, pp. 1879–1884, 2004.

M. L. Shi, B. Zhang, and Z. S. Hong, “Mechanical properties of natural sedimentary intermediate soils,” Rock and Soil Mechanics, vol. 26, no. 11, pp. 1753–1756, 2005.

J.-H. Yin, “Non-linear creep of soils in oedometer tests,” Géotechnique, vol. 49, no. 5, pp. 699–707, 1999.

C. W. W. Ng, Q. Li, and G. B. Liu, “Characteristics of one-dimensional compressibility of Shanghai clay,” Chinese Journal of Geotechnical Engineering, vol. 33, no. 4, pp. 630–636, 2011.

G. Mesri and P. M. Godlewski, “Time and stress compressibility interrelationship,” Journal of Geotechnical Engineering, vol. 103, no. 5, pp. 417–430, 1979.

X. C. Shi, R. Wang, and Z. X. Zhang, “Study on secondary consolidation properties of Guangxi marine ooze,” Rock and Soil Mechanics, vol. 24, no. 5, pp. 863–865, 2003.

H. M. Zhang, Y. S. Xu, and Q. L. Zeng, “Deformation behavior of Shenzhen soft clay and post-construction settlement,” Chinese Journal of Geotechnical Engineering, vol. 24, no. 4, pp. 509–514, 2002.

L. C. Miao, J. H. Zhang, and Y. N. Chen, “Study on compressibility of Jiangsu marine clay,” Chinese Journal of Geotechnical Engineering, vol. 29, no. 11, pp. 1711–1714, 2007.

F. Z. Zhang, X. P. Chen, and G. Y. Huang, “Experimental investigation of consolidation behavior of the pearl delta saturated soft clay,” Rock and Soil Mechanics, vol. 24, pp. 192–194, 2003.

Q. J. Zhou and X. P. Chen, “Test study on properties of secondary consolidation of soft soil,” Rock and Soil Mechanics, vol. 27, no. 3, pp. 404–408, 2006.

P. Z. Lu, J. Zeng, and Q. Sheng, “Creep tests on soft clay and its empirical models,” Rock and Soil Mechanics, vol. 29, no. 4, pp. 1041–1045, 2008.

H. Y. Lei, J. J. Liu, Y. F. Jia et al., “Triaxial test on creep properties of coastal soft clay,” Journal of Tianjin University (Science and Technology), vol. 46, no. 5, pp. 387–392, 2013.

Y. K. Liu, Z. B. Deng, P. Cao et al., “Triaxial creep test and modified Singh-Mitchell creep model of soft clay,” Journal of Central South University (Science and Technology), vol. 43, no. 4, pp. 1440–1446, 2012.

X. L. Lai, S. M. Wang, W. M. Ye, and Y. J. Cui, “Experimental investigation on the creep behavior of an unsaturated clay,” Canadian Geotechnical Journal, vol. 51, no. 6, pp. 621–628, 2014.

L. Z. Wang and Z. Y. Yin, “Stress dilatancy of natural soft clay under an undrained creep condition,” International Journal of Geomechanics, vol. 15, no. 5, pp. 1–5, 2015.

S. W. Yan, K. J. Liu, W. Li et al., “Study of creep properties of soft clay in Tianjin Binhai new area and no-yield-surface constitutive model,” Rock and Soil Mechanics, vol. 31, no. 5, pp. 1431–1436, 2010.

M.-Q. Sun, Q. Wang, C.-C. Niu, and T. Sun, “Research on the one-dimensional rheological consolidation theory that considers secondary consolidation effect,” Journal of Computational and Theoretical Nanoscience, vol. 13, no. 2, pp. 1136–1146, 2016.

Y.-b. Zhu and H.-m. Yu, “An improved Mesri creep model for unsaturated weak intercalated soils,” Journal of Central South University, vol. 21, no. 12, pp. 4677–4681, 2014.

S.-y. Feng, L.-m. Wei, C.-y. He, and Q. He, “A computational method for post-construction settlement of high-speed railway bridge pile foundation considering soil creep effect,” Journal of Central South University, vol. 21, no. 7, pp. 2921–2927, 2014.

M.-j. Jiang, J.-d. Liu, and Z.-y. Yin, “Consolidation and creep behaviors of two typical marine clays in China,” China Ocean Engineering, vol. 28, no. 5, pp. 629–644, 2014.

M. Q. Sun, Q. Wang, Y. K. Ruan et al., “Settlement prediction of soft soil foundation based on consolidation and creep theory,” International Journal of Simulation: System, Science and Technology, vol. 16, no. 2, pp. 21–28, 2015.

A. M. Shaban and P. J. Cosentino, “Modeling long-term deformations of unbound pavement materials using the miniaturized preburemeter creep data,” Geotechnical Testing Journal, vol. 39, no. 5, pp. 813–826, 2016.

A. Singh and J. K. Mitchell, “General stress–strain–time function for soils,” Journal of the Soil Mechanics and Foundation Divisions, vol. 94, no. 1, pp. 21–46, 1968.

G. Mesri, E. Febres-Cordero, D. R. Shields, and A. Castro, “Shear stress-strain-time behaviour of clays,” Géotechnique, vol. 31, no. 4, pp. 537–552, 1981.