Abstract

This paper presents analytical solutions to Fredlund and Hasan’s one-dimensional consolidation equations for unsaturated soils subjected to various cyclic loadings. Two new variables are introduced so that the governing equations for excess pore air and water pressures can be transformed to a set of conventional diffusion equations. Based on the general solutions for two introduced variables, the analytical solutions are derived for one-dimensional consolidation of unsaturated soils under trapezoidal, rectangular, triangular, and haversine cyclic loadings. It shows through the degeneration into the existing solutions for unsaturated and saturated soils that the proposed solutions are more general ones for one-dimensional consolidation of soils from unsaturated to saturated states. A comprehensive parametric study is conducted to investigate the effects of different parameters on one-dimensional consolidation of unsaturated soils under various cyclic loadings. The proposed solutions can be effectively utilized in the analysis of consolidation of unsaturated soils subjected to various cyclic loadings.

1. Introduction

In geotechnical engineering, soils beneath many structures may be subjected to complicated cyclic loadings. These cyclic loadings are induced by different vibration sources, such as filling and discharging in silos, tanks, and reservoirs; traffic loads on highways, railways, or airports; wave action for offshore structures; wind forces in high buildings; pile construction; and machine vibrations. The consolidation characteristics of soils subjected to cyclic loadings are very complex and different from those of soils under simple time-dependent loadings. Thus, consolidation behavior of subsoils under cyclic loadings has drawn considerable attention and has become the focus of academic research.

Over the past few decades, numerous studies have been conducted on consolidation behavior of saturated soils under various cyclic loadings. The existing investigation of consolidation behavior of saturated soils under cyclic loadings can be divided into two groups. One focused on the consolidation based on Terzaghi’s linear consolidation theory, while the other was based on the nonlinear consolidation theory. Preceding works on the former included the efforts by Wilson and Elgohary [1], Favaretti and Soranzo [2], Conte and Troncone [3], Razouki and Schanz [4], Razouki et al. [5], and Kim et al. [6]. Wilson and Elgohary [1] presented an analytical solution for one-dimensional consolidation of saturated soil subjected to rectangular cyclic loading, whereas Favaretti and Soranzo [2] and Conte and Troncone [3] derived solutions for different types of cyclic loadings. Razouki and Schanz [4] and Razouki et al. [5] studied one-dimensional consolidation of a clay layer subjected to haversine cyclic loading with and without rest period. Kim et al. [6] investigated one-dimensional consolidation behavior of saturated clays under time-dependent loadings including haversine cyclic loading based on the lattice Boltzmann method. The latter mainly included the work by Baligh and Levadoux [7], Xie et al. [8], Geng et al. [9], Cai et al. [10], Conte and Troncone [11], Toufigh and Ouria [12], Cheng et al. [13], Hu et al. [14], and Kim et al. [15, 16]. Baligh and Levadoux [7] and Toufigh and Ouria [12] analyzed one-dimensional nonlinear consolidation of inelastic clays under rectangular cyclic loading using the superimposing rule. Based on the nonlinear consolidation theory proposed by Davis and Raymond [17], Xie et al. [8], Conte and Troncone [11], Cheng et al. [13], and Kim et al. [15] investigated one-dimensional nonlinear consolidation of saturated soil subjected to various cyclic loadings. Hu et al. [14] proposed a numerical method for modelling one-dimensional nonlinear consolidation of multilayer soil under various cyclic loadings based on the differential quadrature method. Geng et al. [9] and Cai et al. [10] developed a semianalytical method for one-dimensional nonlinear consolidation of soil with variable compressibility and permeability under cyclic loadings. More recently, Kim et al. [16] derived analytical solutions for one-dimensional nonlinear consolidation of a saturated clay layer under cyclic loadings with the consideration of variable consolidation coefficient.

Most soils in nature are unsaturated and develop much more complex properties during the consolidation process than saturated soils. In the recent decades, a number of studies on the consolidation theory for unsaturated soil has been conducted, and a great progress has been achieved. Blight [18] derived a consolidation equation for the air phase of dry, rigid, unsaturated soil, whereas Barden [19] presented an analysis of the one-dimensional consolidation of compacted unsaturated clay. By assuming that the water and air phases are continuous, Fredlund and Hasan [20] presented an uncoupled one-dimensional consolidation theory for unsaturated soil, in which two partial differential equations for water and air phases were derived. Meanwhile, some scholars [2123] proposed coupled models for consolidation in unsaturated soil. Although coupled models are theoretically more rigorous than uncoupled ones, the uncoupled theory proposed by Fredlund and Hasan [20] is now widely accepted because it is generally simple and sufficiently accurate for practical projects [24, 25]. On the basis of the consolidation theory proposed by Fredlund and Hasan [20], many researchers have investigated consolidation of unsaturated soil under different time-dependent loadings [2643]. Among them, Feng et al. [30] presented an analytical solution for one-dimensional consolidation of unsaturated soil under sinusoidal cyclic loading. Ho and Fatahi [31] and Wang et al. [32] suggested analytical and semianalytical solutions to Fredlund and Hasan’s one-dimensional consolidation of unsaturated soils subjected to different time-dependent loadings including sinusoidal loading. On the contrary, based on the coupled consolidation model of poroelasticity, Lo et al. [44] proposed a closed-form analytical solution for one-dimensional consolidation of unsaturated soil under harmonic cyclic loading.

Despite these advances, it seems that research on consolidation of unsaturated soils subjected to different cyclic loadings is very limited. Most researchers performed investigations on consolidation behavior of unsaturated soils subjected to simple time-dependent loadings, such as constant, ramp, exponential, and a few sine loadings, but little research has been done on consolidation behavior of unsaturated soils under various cyclic loadings.

In this paper, based on the consolidation theory proposed by Fredlund and Hasan [20], analytical solutions are derived for one-dimensional consolidation of unsaturated soils under various cyclic loadings, such as trapezoidal, rectangular, triangular, and haversine cyclic loadings. The proposed solutions are degenerated into the existing solutions for unsaturated soils under constant and ramp loadings and for saturated soils under haversine cyclic loading. Then, a comprehensive parametric study is conducted to investigate the effects of different parameters on one-dimensional consolidation of unsaturated soils under various cyclic loadings.

2. Governing Equations

An unsaturated soil layer of thickness subjected to time-dependent loading is considered as shown in Figure 1. In order to obtain the analytical solution for one-dimensional consolidation of the unsaturated soil layer subjected to time-dependent loading, the main assumptions are summarized as follows:

(1)The flows of air and water phases are continuous and independent(2)The soil particles and water phase are incompressible(3)The effects of temperature change, air diffusing through water and air dissolving in the water, and the movement of vapor are ignored(4)The coefficients of permeability and volume change for both air and water phases remain constant during the consolidation process(5)The deformation takes place only along the vertical direction

The governing equations for one-dimensional consolidation of unsaturated soils subjected to time-dependent loading are as follows [19, 23]:where and are the excess pore air and pore water pressures, respectively. The consolidation parameters can be expressed as follows:where and are the coefficients of air and water volume change with respect to a change in the net normal stress , respectively; and are the coefficients of air and water volume change with respect to a change in the matric suction , respectively; and are the coefficients of permeability for air and water phases, respectively; is the gravitational acceleration; is the absolute pore air pressure, , in which is the atmospheric pressure, 101 kPa; when the pore air pressure is small or rapidly dissipates during consolidation, can be considered constant, and therefore, in this study, ; is the molecular mass of air, 0.029 kg/mol; is the universal gas constant, 8.314 J/mol/K; is the absolute temperature, (20 + 273.16)K; is the initial porosity; is the initial degree of saturation; and is the unit weight of water.

The initial excess pore air and pore water pressures can be expressed as follows:where and are the initial excess pore air and pore water pressures, respectively.

The boundary conditions for the excess pore air and pore water pressures can be given by the following equations:

Equations (1) and (2) can be presented in nondimensional form by using nondimensional parameters:where  = ultimate value of loading , ,  = coefficient of volume change with respect to a change in net normal stress , ,

The dimensionless initial and boundary conditions can be expressed as follows:where .

3. Derivation of Analytical Solutions under Cyclic Loadings

According to the simple analytical solutions presented by Zhou et al. [29], (8) and (9) can be transformed to a set of conventional diffusion equations by introducing two new variables and :where

The corresponding initial and boundary conditions for excess pore air and pore water pressures can be transformed as follows:where .

For the single drainage condition, the general solutions of (10) subjected to the solution conditions of (12)–(14) can be obtained as follows:where , , , and , in which

The solutions for the double drainage condition can be easily obtained by replacing in (16) with considering the drainage height (i.e., ).

The complete solutions of for various cyclic loadings can be calculated by substituting the different loading functions into (16).

On the basis of the definition of the average degree of consolidation, the degree of consolidation, , with respect to for different cyclic loadings can be expressed as follows:

The detailed mathematical process of obtaining the degree of consolidation, , can be found in Appendix.

3.1. Trapezoidal Cyclic Loading

Trapezoidal cyclic loading shown in Figure 2(a) can be expressed as follows:where is the period of one loading cycle; and are the loading parameters corresponding to the rate of loading increment or decrement and the rest period of loading, respectively; and is the cycle number.

Transforming (18) into nondimensional form and substituting into (16), the complete solutions of (10) can be obtained as follows:where

Substituting (19) into (17), the solutions of for trapezoidal cyclic loading can be given as follows:where

3.2. Rectangular and Triangular Cyclic Loadings

Rectangular and triangular cyclic loadings can be degenerated from trapezoidal cyclic loading. When , trapezoidal cyclic loading reduces into rectangular cyclic loading (Figure 2(b)). In addition, trapezoidal cyclic loading reduces to triangular cyclic loading (Figure 2(c)) when . Therefore, the solutions of and for rectangular and triangular cyclic loadings can be given by substituting and into (19) and (21), respectively.

3.3. Haversine Cyclic Loading

Haversine cyclic loading shown in Figure 2(d) can be expressed as follows:

Transforming (22) into nondimensional form and substituting into (16), the complete solutions of (10) can be obtained as follows:where

Substituting (23) into (17), the solutions of for haversine cyclic loading can be given as follows:

3.4. Excess Pore Pressures and Average Degree of Consolidation

Based on the solutions of , the normalized excess pore air and pore water pressures can be given by following relations:

Based on the solutions of for different cyclic loadings, the average degrees of consolidation for both air and water phases can be expressed as follows:where and are the average degrees of consolidation for air phase and water phase, respectively.

4. Degeneration of the Proposed Solution to Special Cases

The proposed solutions can be degenerated into the existing solutions for one-dimensional consolidation of both unsaturated and saturated soils subjected to different time-dependent loadings.

4.1. Constant Loading

When , and , trapezoidal cyclic loading reduces to constant loading (Figure 3(a)) expressed by . Then, the solutions of and can be given by

Equations (29) and (30) are the solutions of and for one-dimensional consolidation of unsaturated soil under constant loading and constant initial condition proposed by Zhou et al. [29].

It can be seen from Figure 4 that the degenerate solutions from this paper are consistent with the existing solutions and the proposed solutions are accurate.

4.2. Ramp Loading

When and tends to infinity, trapezoidal cyclic loading reduces to ramp loading (Figure 3(b)) given bywhere .

The solutions of and can be obtained bywhere and .

Equations (32) and (33) are the solutions of and for one-dimensional consolidation of unsaturated soil subjected to ramp loading proposed by Zhou et al. [29].

Figure 5 shows the comparison between the degenerate solutions and the existing solutions by Zhou et al. [29] for ramp loading. It can be found that the degenerate solutions from this paper show good agreement with the existing results, which proves that the proposed solutions are more general ones for one-dimensional consolidation of unsaturated soil under time-dependent loading.

4.3. Degeneration to Saturated Condition for Haversine Cyclic Loading

The consolidation coefficients for unsaturated soil can be reduced to those for saturated soil: the coefficients of water volume change ( and ) are equal to the coefficient of volume change (); the coefficients of air volume change ( and ) are zero; air pressure may be assumed to be equal to water pressure, that is, ; , , , , , and . Then, the corresponding parameters in (8) and (9) become , , , , , and . In addition, , , , , , and . It follows from (26) and (27) that . The solution for the excess pore water pressure of saturated soil with single drainage condition under haversine cyclic loading can be obtained as follows:where

Equation (34) is the solution for one-dimensional consolidation of saturated soil subjected to haversine cyclic loading proposed by Razouki et al. [5].

As observed in Figure 6, identical results are obtained from both analytical solutions, suggesting that the proposed solutions are reliable as well as more general for one-dimensional consolidation of a soil from unsaturated to saturated states.

5. Parametric Study

In this section, a comprehensive parametric study is conducted to investigate the effects of the coefficients ratio of permeability for air and water phases and the other loading parameters on consolidation behavior of unsaturated soils with single drainage condition under various cyclic loadings. The material properties adopted in this study are as follows:

The initial excess pore air and pore water pressures are assumed to be zero (i.e., and ) because of for cyclic loadings except for rectangular cyclic loading, while and (constant with depth) for rectangular cyclic loading and constant loading.

5.1. Effects of

Figures 7 and 8 represent the variation in the average degrees of consolidation for both water and air phases with time at different values of under various cyclic loadings. It can be found that the average degrees of consolidation and increase but fluctuate with the increase in , indicating that the dissipation rate of excess pore pressures increase with the fluctuation when increases. Moreover, it can be noticed that the higher values of lead to the bigger amplitudes of fluctuation in the average degrees of consolidation and .

It can be observed through the locally enlarged diagrams with different cyclic loading forms that there is similar characteristics for the dissipation process of excess pore pressures and the change in applied loading. The variation in the average degrees of consolidation with time shows the oscillation with the same frequency and the similar forms in response to the applied cyclic loadings. It can also be seen that the dissipation process of excess pore pressures does not complete after a long time but rather continues oscillating with a certain amplitude. It is due to repetitive loading and unloading stages of cyclic loadings, which results in squeezing out of pore water and air during the loading stage and then absorbing during the unloading stage.

5.2. Effects of

Figure 9 demonstrates the variation in average degrees of consolidation for both water and air phases with time under trapezoidal cyclic loading ( and ) and haversine cyclic loading with different loading parameters for . It can be seen that a smaller value of the loading parameter induces a smaller amplitude and more cycles of fluctuation in the dissipation rate of excess pore water and air pressures. In addition, the loading parameter has a significant influence on the amplitude and the cycles of fluctuation in the dissipation rate of excess pore pressures but does not affect the mean values of the dissipation rates. The loading parameter is the period of loading, and thus, it can be concluded that the longer the period of loading, the bigger the amplitude, and the less the cycles of fluctuation in the dissipation rates.

5.3. Effects of

Figure 10 presents the variation in average degrees of consolidation for both water and air phases with time under different types of cyclic loadings with and different loading parameters for . It can be seen that the smaller the loading parameter is, the more quickly the dissipation of excess pore water and air pressures develops. Since the loading parameter reflects the rate of loading increment or decrement, it can be found that the higher the rate of loading increment or decrement is, the faster the dissipation rates of excess pore pressures are. Moreover, the dissipation rates of excess pore pressures under rectangular cyclic loading are the biggest and the dissipation rates under triangular cyclic loading are the smallest. For the comparison with the results for static loading, the variation in average degrees of consolidation under constant loading is also shown in Figure 10. It can be observed that the average degrees of consolidation under cyclic loadings increase but fluctuate with a certain amplitude, whereas the average degrees of consolidation under constant loading increase monotonically with time. This phenomenon implies that consolidation of unsaturated soil under cyclic loadings does not complete after a long time, but rather continues oscillating with a certain amplitude, generating continuous volume change of the soil. Thus, neglecting the effect of cyclic loadings may result in overestimated settlement and cannot consider the continuous volume change.

5.4. Effects of

Figure 11 depicts the variation in average degrees of consolidation for both water and air phases with time under different cyclic loadings with different loading parameters for . It can be found that the dissipation rates of excess pore water and air pressures increase with the decrease in the value of loading parameter . Furthermore, a smaller value of the loading parameter results in more cycles and a smaller amplitude of fluctuation in the dissipation rates. The loading parameter is related to the rest period of loading. Thus, it is quite obvious that a longer rest period of loading induces the decrease in the dissipation rates, less cycles, and a bigger amplitude of fluctuation in the dissipation rates of excess pore pressures. As a result of that the dissipation rates of excess pore water and air pressures under the cyclic loading without rest period () are the biggest for each case. Specifically, rectangular cyclic loading without rest period reduces into constant loading, and the dissipation rates of excess pore water and air pressures are the highest and increase monotonically with time in that case, as shown in Figure 11(b).

It should be noted that the solutions presented in this paper do not consider nonlinear characteristics of soil properties such as the coefficients of permeability and volume change for both air and water phases under time-dependent loading. Even though they are limited for a linear elastic stress-strain model of unsaturated soil based on Fredlund and Hasan’s one-dimensional consolidation theory, the proposed solutions are more general ones for one-dimensional consolidation of soils from unsaturated to saturated state under time-dependent loading.

6. Conclusion

In this paper, analytical solutions are derived for one-dimensional consolidation of unsaturated soils subjected to various cyclic loadings. It shows through the degeneration into the existing solutions for both unsaturated and saturated soils that the proposed solutions are reliable and more general for one-dimensional consolidation of soils from unsaturated to saturated states under time-dependent loading. The effects of different parameters on one-dimensional consolidation of unsaturated soils under various cyclic loadings are investigated. The main conclusions can be drawn as follows:(1)The coefficients ratio of permeability for air phase and water phase, has sufficient effects on consolidation of the unsaturated soil layer under cyclic loadings. The dissipation rates of excess pore water and air pressures increase with the fluctuation when increases. Moreover, the higher values of lead to the bigger amplitudes of fluctuation in the dissipation rates.(2)Under cyclic loadings, the dissipation process of excess pore pressures does not complete after a long time, but rather continues oscillating with a certain amplitude, generating continuous volume change in the soil. Thus, neglecting the effect of cyclic loadings may result in overestimated settlement and cannot consider the continuous volume change.(3)The cyclic loading parameters such as the period of loading (), the rate of loading increment or decrement (), and the rest period of loading () have considerable impacts on consolidation of unsaturated soil. The longer the period of loading, the bigger the amplitude, and the less the cycles of fluctuation in the dissipation rates of excess pore water and air pressures. The higher the rate of loading increment or decrement is, the faster the dissipation rates of excess pore pressures are. A longer rest period of loading induces the decrease in the dissipation rates, less cycles, and a bigger amplitude of fluctuation in the dissipation rates of excess pore pressures.

Appendix

The average degrees of consolidation for air and water phases under time-dependent loading are defined as follows:

Equations (A.1) and (A.2) can be expressed in terms of based on equations (26) and (27) as follows:

Transforming (A.3) and (A.4) mathematically,

Therefore, (A.5) and (A.6) can be rewritten in the following form:where and are the degrees of consolidation with respect to and defined as follows:

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

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

Acknowledgments

This work was supported by Kim II Sung University in DPR of Korea.