Abstract

Creep is a common phenomenon for soft clays. The paper focuses on investigating the influence of temperature on the time-dependent stress-strain evolution. For this purpose, the temperature-dependent creep behavior for the soft clay has been investigated based on experimental observations. A thermally related equation is proposed to bridge the thermal creep coefficient with temperature. By incorporating the equation to a selected one-dimensional (1D) elastic viscoplastic (EVP) model, a thermal creep-based EVP model was developed which takes into account the influence of temperature on creep. Simulations of oedometer tests on reconstituted clay are made through coupled consolidation analysis. The bonding effect of the soil structure on compressive behavior for intact clay is studied. By incorporating the influence of the soil structure, the thermal creep EVP model is extended for intact clay. Experimental predictions for thermal creep oedometer tests are simulated at different temperatures and compared to that obtained from reconstituted clay. The results show that the influence of temperature on the creep behavior for intact clay is significant, and the model, this paper proposed, can successfully reproduce the thermal creep behavior of the soft clay under the 1D loading condition.

1. Introduction

It is well known that soft clays exhibit time-dependent behavior due to their viscosity. The long-term settlement of these clays after the dissipation of excess pore water pressure, which is sometimes called creep deformation, has been an important issue in geotechnical engineering. The creep behavior of the soft clay has been investigated experimentally in [15]. Based on that, some practical models have been developed [69].

Due to the deposition effect, interparticle bonds are usually formed in soft clays referred as the soil structure. When suffering loading for the soft clay, a significant progressive loss of bonds will happen. By comparing the compression curves in (void ratio versus vertical effective stress in the log scale) for intact and reconstituted samples, the large differences observed are induced by the bond elimination. In addition, the bond elimination effect on the creep behavior of soft clays has been studied in [10, 11].

Soft clays are also subjected to the action of heat under many circumstances, for example, the nuclear waste isolation, heat energy storage, and geothermal development. Studies show that the creep behavior of soft clays is strongly related not only to the bonding structure but also to temperature. The creep behavior of these clays will be changed accordingly [1216]. However, the effect of temperature on the creep behavior was somehow showed, but not well documented. It will be nice if there exists a direct way for modeling the thermal creep behavior.

For this purpose, we focus on the 1D behavior which can bring fundamental features for more mechanical behavior. Firstly, the temperature dependency behavior of creep is studied based on the experimental observations. Then, a thermally related EVP model for reconstituted clay is proposed by incorporating the effect of temperature. Furthermore, the bonding elimination of the structure is incorporated into the thermal-based model. Finally, the prediction ability of the proposed model is shown by simulations at different conditions.

2. Temperature Dependency of Creep of Soil

2.1. Experimental Observations

The consolidation process of soil is usually divided into primary consolidation and secondary consolidation, and the boundary point is whether the excess pore pressure completely dissipates. Similarly, thermal secondary consolidation occurs for thermal consolidation after the dissipation of excess pore pressure [13]. Figure 1 shows an idealized thermal consolidation test result. Actually, thermal secondary consolidation occurs throughout the thermal consolidation process, and the induced deformation only relates to time and temperature. Thus, thermal creep used in this paper is more suitable for the process. The slope of the linear portion of the thermal consolidation curve is the thermal creep coefficient , given in void ratio per log cycle.

A number of studies demonstrate that the temperature affects significantly the thermal creep coefficient. For example, Figure 2 presents the evolution of the thermal creep coefficient with increasing temperature for intact Pacific illite conducted by Houston et al. [13]. The values range from about 0.01 at 40°C to about 0.06 at 200°C. The increase in the rate of thermal creep deformation at the elevated temperatures was quite significant. Additionally, the creep-temperature tests on peats conducted by Fox and Edil [12] also show that creep dominates the consolidation process and temperature influences the creep rate significantly. Thus, it can be concluded that it is necessary to account for the effect of temperature on creep for soft intact clay.

2.2. Proposed Thermal Creep Equation

Based on the experiments, different expressions were given to describe the temperature-dependent behavior of the thermal creep coefficient. The relationship between thermal creep coefficient and temperature can be linearly functioned by Houston et al. [13]. That is,where is the thermal creep coefficient under temperature and is the thermal relate parameter, which can be obtained by correlating the results as shown in Figure 2. Equation (1) describes increasing linearly with temperature straightforwardly. However, we can observe that will be zero when the temperature decreases to zero and will be negative if the temperature still decreases and is unreasonable.

Based on the experimental results on peat, the following equation is used by Fox and Edil [12] to predict the value of , due to a change in soil temperature:where is the reference thermal creep coefficient under the reference temperature , thermal relate parameter B is equal to 0.25 0.02/°C for the peat in [12], and the parameter is independent of vertical effective stress and the magnitude of temperature change. Equation (2) describes that increases nonlinearly with temperature. However, due to its mathematic structure, Equation (2) will deduce a very large number of at a higher temperature, far beyond the normal range.

In this study, we propose a new equation to describe the nonlinear increase of the thermal creep coefficient with temperature, which can be written aswhere is the thermal relate parameter which can be correlated with the experiment results. This equation can overcome the deficiency exposed by Equations (1) and (2). Figure 3 plots the comparisons of the three equations describing the relationship between thermal creep coefficient and temperature. The point = 20°C and (20) = 0.006 is adopted as the reference. As indicated, three equations give different shape curves. Overall, the curve plotted from Equation (3) is more practical than the others.

3. EVP Model considering the Temperature Effect on Creep

3.1. Adopted One-Dimensional Elastic Viscoplastic Model

First, a time-dependent stress-strain model without the thermal effect needs to be selected as a base [6]. The adopted model is presented briefly in this part.

Following the classic elastic viscoplastic approach, the total strain rate contains the elastic and viscoplastic strains rate, that is,where represents the total strain rate and the superscripts “” and “” denote the elastic and viscoplastic parts, respectively. The elastic strain rate is expressed aswhere is the initial void ratio, is the effective vertical stress, and is the slope of recompression lines in space.

For the viscoplastic strain rate, a one-dimensional formulation proposed by Kutter and Sathialingam [17] was adopted based on the creep coefficient :where denotes the slope of the normal compression line in space; is the reference time, and it equals to the duration of each load increment in the oedometer test; and is the reference stress corresponding to the incremental time and increases with the development of the viscoplastic strain according towhere denotes the preconsolidation pressure.

The above relationships have been suggested in [1820] and validated in [6, 10].

3.2. Incorporation of the Thermal Effect

The effects of temperature on the stress-strain behavior of clay have been observed in the laboratory [12, 13], which are helpful to discuss the thermal effect on the parameters in the above constitutive model. Increase and decrease in temperature may produce changes in the bonding of clay particles and the viscosity of absorbed water. These changes alternately produce more or fewer changes in compressibility. However, the experiments conducted in [15, 21, 22] show that the variations of and with temperature are negligible. Furthermore, a volume change due to elastic expansion of the clay particle will occur during the drainage of thermal consolidation. Considering that the strain under a constant effective stress remains small [14, 23, 24] and the emphasis of this paper, the parameters , , and will remain constant when clays suffering a change of temperature.

The thermal creep coefficient can be incorporated into the basic EVP model directly. By substituting the parameter in Equation (6) by Equation (3), the 1D thermal-related viscoplastic strain rate changes to

Actually, the preconsolidation pressure also varies with temperature. The thermal-dependent behavior of has been studied from the oedometer tests or isotropic compression tests with variable temperatures. All of the results indicate that will decrease with an elevated temperature [2529]. Based on the data collected from literature, Figure 4 plots the relationship between the normalized preconsolidation pressure and temperature. The regression analyses show that it is also reasonable to assume a linear relationship between and , firstly proposed by [28]. Thus, the relationship between the preconsolidation pressure and temperature can be fitted bywhere is a thermal parameter and and are the preconsolidation pressures at temperature and the reference temperature , respectively. The present model has no elastic limit, which is different from Perzyna’s overstress method [30].

3.3. Coupled Consolidation Analysis

In the following numerical simulations of the conventional oedometer tests, soil-water coupling analysis will be performed. Darcy’s law was adopted for the consolidation process:where is the vertical depth, is the excess pore pressure, and is the hydraulic conductivity. Actually, is also influenced by temperature. A decrease in pore water viscosity will happen with increasing temperature, which will result in an increase in the permeability of the soil. Thus, increasing temperature will speed up the consolidation process. Considering the emphasis of this paper, only the void ratio-dependent behavior of k is considered, and according to experimental results, the evolution of can be expressed aswhere the initial hydraulic conductivity is corresponding to and the permeability coefficient can be easily measured from the oedometer test results by plotting .

The above equations were implemented in finite element software PLAXIS Version 8, but in 1D, finite element simulations can be established for modeling the primary consolidation and the creep process in an oedometer test. Details of coupled consolidation and creep analysis can be found in [19, 31] and are not repeated here.

3.4. Simulated One-Dimensional Behavior

In order to validate the proposed thermal-dependent EVP model, numerical simulations for the assumed oedometer creep tests were performed at three temperatures ( = 20°C, 50°C, and 80°C). The results are shown in Figure 5, and the parameters adopted for these simulations are listed in Figure 5(a). The simulated temperature behavior agrees with the common experimental phenomena on unstructured clay, as expected by the model’s principle. For example, the simulated relationship between the preconsolidation pressure and temperature corresponds to the input value of .

Take the simulated test on = 80°C; for example, the compression behavior for each load increment is presented in Figure 5(b). The simulated thermal creep coefficient = 0.013 agrees well with that obtained from Equation (3) with = 20°C and = 0.0065.

4. Thermal-Dependent EVP Model for Intact Clay

4.1. Bonding Effect on Compression Behavior

During the oedometer tests, the difference of the compression curves obtained on intact and reconstituted clay is caused by bonding elimination as shown in Figure 6(a) for Wenzhou clay [32]. The structures between soil particles for intact clay will be eliminated gradually during compressing. The tests conducted under variable temperatures show that the shape of the compression curves does not change with temperature, for instance, tests on intact Berthierville clay [25] and Linköping clay [28]. Thus, we assume that the process of bonding elimination is thermal independent but only relates to the strain level. Figure 6(b) presents the schematic plot of the stress-strain curve at an arbitrary temperature for soft intact clay. For a given viscoplastic strain level , the bond elimination results in the current stress reaching point for intact clay. Corresponding to the same viscoplastic strain, we define an intrinsic stress on the reconstituted sample. We assume that the difference between the current stress and intrinsic stress is due to the existing of the soil structure, based on which a bonding ratio can be proposed. Thus, the current stress during straining can be expressed as

Initially, the bonding ratio . Following the increasing of strain, the bonds are broken gradually and decreases from its initial value ultimately towards zero when the bonds are completely eliminated as plotting in Figure 6(a). According to the definition, bonding ratio and the corresponding viscoplastic strain during compression is measured and plotted in Figure 6(c) for Wenzhou clay. Based on the results, we propose the following relationship to express the attenuation of bonding ratio:where the parameter controls the rate of bonding elimination ( = 8.0 for Wenzhou clay in Figure 6(c)). Actually, the intrinsic stress in Equation (12) can be regarded as the reference stress as indicated in Equation (8), and the bonding ratio can be regarded as the scaling parameter. Thus, the present model is then composed of Equations (8), (9), and (12). Combining with the elastic strain rate in Equation (5), the stress-strain curve for a given temperature can be obtained.

4.2. Model Parameters

The present model that combined creep and temperature involves a number of parameters which can be divided into four groups:(a)Parameters related to compressibility: initial void ratio (), the intrinsic slope of the compression line (), and the slope of the recompression line (). The values of and can be measured from the oedometer tests on the reconstituted and intact samples, respectively. As the thermal expansion and contraction are neglected in this paper, can be measured initially at the reference temperature.(b)Parameter related to bonding elimination: the initial bonding ratio () and the parameter . The value of can be measured from the oedometer tests on the intact and reconstituted samples conducted at the same temperature. It needs to point out that, for , high-quality intact samples are needed. The parameter representing the bonding elimination rate can be derived from Equations (12) and (13):where is the volumetric viscoplastic strain corresponding to (Figure 6(b)) at an arbitrary temperature. Thus, can be calculated by selecting a point on the postyield curve.(c)Parameters of creep: reference thermal creep coefficient () and thermal parameter (). can be measured directly from the oedometer tests on reconstituted at a reference temperature (). can be obtained by correlating with temperature.

Parameters related to preconsolidation pressure: the reference preconsolidation pressure () at the reference temperature and the thermal parameter (). can be obtained at the intersection of the compression curves for the reconstituted and intact samples, as shown in Figure 6(b). can be obtained directly from the oedometer tests on the reconstituted or intact samples at different temperatures. Wang et al. [33] investigated the value of for seven clays and summarized that varies from 0.125 to 0.194. Furthermore, can be obtained by the empirical correlation of liquid limit () expressed as

Thus, can be obtained by correlating without carrying out the temperature-controlled tests.

Thus, ten model parameters (, , , , , , ,, , and ) are required for the model, and all of the parameters can be determined straightforwardly from the temperature-controlled oedometer test.

5. Calibration and Validation

5.1. Predictions on Wenzhou Clay

The thermal creep coefficient is also influenced by the bonding elimination process during straining. In order to evaluate the model ability to reproduce the behavior of soft intact clay, simulations were performed for the thermal oedometer tests with different temperatures. The input ordinary parameters are the same as those given in Table 1. For intact clay, two parameters = 3.5 and = 8 were used. Three temperatures with = 20°C, 50°C, and 80°C are adopted in the simulations.

Figure 7 shows typical results of the thermal oedometer tests in natural intact clay simulated by the new model:(1)In Figure 7(a), the compression curves in different temperatures are similar to the oedometer test on intact Wenzhou clay (Figure 6(a)). The preconsolidation pressure is influenced significantly by the temperature as presented in Figure 5 on reconstituted clay.(2)In Figure 7(b), the typical curves of strain versus time were reproduced for natural intact clay, where the thermal creep coefficient can be obtained at the end of each load; here, only the simulation results for = 80°C are presented.(3)In Figure 7(c), the measured at three temperatures were plotted versus the normalized applied stress by the preconsolidation pressure. Here, the preconsolidation pressure is different for the three tests. The value for intact clay increases rapidly with the applied stress, and when the stress reaches the preconsolidation pressure, reaches a peak value and then decreases. The difference between these curves is due to the combined influence of temperature and bonding elimination.(4)In Figure 7(d), the measured for intact and reconstituted clay are compared for the test at =80°C. The difference between the two curves is due to the existence of bonding and its elimination.

5.2. Predictions on Utby Clay

Li et al. [34] presented a set of the long-term oedometer tests on Utby clay (6 m depth) with different temperatures for both intact and reconstituted samples. The clay properties are as follows: liquid limit = 55%, = 1.977, and = 55 kPa for the intact sample at the reference room temperature = 20°. The other parameters for the model were correlated with thermal oedometer tests results. Figure 8(a) shows the compression tests on intact and reconstituted samples at = 20°. The parameter = 0.18 is determined based on the results of the reconstituted sample. The parameter =0.003 is correlated from the recompression curve. Bonding ratio = 54 is obtained based on the method described above. The value of = 15 was determined by selected a stress-strain point in Figure 8(a) and using Equation (14). Figure 8(b) presents the variation of the thermal creep coefficient with temperature. Adopting at = 20° as a reference, the predicted values of for =5° agree well with the experimental results with = 0.31. Furthermore, an average value of = 0.02 is used based on the experimental results at = 20°. With Equation (15), = 0.15 is adopted in the simulation. To simulate the long-term creep test, the permeability of soil was taken as 4 × 10−6 m/h estimated from the consolidation curves of vertical strains versus time. The value of was equal to , as suggested by Tavenas et al. [35] based on observations on soft marine clays. All of the parameters are summarized in Table 1.

Figure 9 shows the comparison between experimental and simulation results. The predicted compression behavior at different temperatures shows good agreement with the experimental results for the values of the preconsolidation pressures and for the shape of the compression curves (Figure 9(a)). Furthermore, the predicted thermal creep coefficient also agrees well with the experimental results. For the values of vertical stress equal to the preconsolidation pressure, the thermal creep coefficient reaches to the maximum value (Figure 9(b)). It can be concluded that it is necessary to account for the coupling effect of temperature and destructuration for accurate predictions of the thermal compression behavior of the soft clay.

5.3. Predictions on Tokyo Bay Clay

Tsuchida et al. [36] conducted the oedometer tests on Tokyo bay clays which are consolidated at room temperature (25°) and at high temperature (75°), respectively. The sample which is suffered consolidating at a high temperature and cooled after the completion consolidation will behave like the lightly aged clay. Adopting the sample consolidated at room temperature as a reference, the properties are as follows: liquid limit = 78% and = 2.4. From the experimental results on the sample consolidated at room temperature (Figure 10(a)), the parameter = 0.36, =0.05, and = 79 kPa can be obtained. In addition, the reference = 0031 is averaged for the vertical stress larger than preconsolidation pressure (Figure 10(b)). Also, the parameter = 0.17 can be calculated by Equation (15). = 0.2 is obtained by the increase of preconsolidation pressure due to the cooling of the sample after high temperature consolidation. Similarly, = 7 is correlated based the compression curve, and the consolidation-related parameters = 1 × 10−7 m/h, = 1.2. All of the parameters are collected in Table 1.

Figure 11(a) shows that the predicted compressibility of Tokyo bay clay consolidated at a high temperature agrees well with the experiment. Furthermore, the predicted thermal creep coefficient has the same shape with experiment. The predicted values are a little smaller than the experimental results (Figure 11(b)). Overall speaking, the model can reproduce well the thermal creep behavior for the soft clay.

6. Conclusions

The temperature-dependent behavior of creep for soft intact clay has been investigated based on experimental observations from experimental results. A thermally related equation is proposed to bridge the thermal creep coefficient with temperature. By incorporating the equation to a selected one-dimensional EVP model, a thermal creep-based EVP model is developed taking into account the temperature dependency of creep. The determination of the model parameters is straightforward. Numerical simulations have been conducted to examine the predictive ability of the model for the soft clay.

Experimental predictions have carried out the thermal oedometer tests at different temperatures. The bonding elimination effect on the evolution of the thermal creep coefficient has been highlighted by comparing predictions with and without considering bonding elimination.

The results demonstrate that the proposed model can well reproduce temperature-dependent creep behavior of soft intact clay under the one-dimensional loading condition. Future work will be done to extend the proposed model to three-dimensional general stress space.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interests regarding the publication of this article.

Acknowledgments

We acknowledge with gratitude the financial support provided by the Fundamental Research Funds for the Central Universities (2015QNA64).