One-Dimensional Consolidation of Double-Layered Foundation with Depth-Dependent Initial Excess Pore Pressure and Additional Stress

A model for one-dimensional consolidation of a double-layered foundation considering the depth-dependent initial excess pore pressure and additional stress and time-dependent loading under different drainage conditions was presented in this study and its general analytical solution was deduced. The consolidation solutions of several special cases of single-drained and doubledrained conditions under an instantaneous loading and a single-level uniform loading were derived. Then, the average degree of consolidation of the double-layered foundation defined by settlement was gained and verified. Finally, the effects of the initial excess pore pressure distributions, depth-dependent additional stress, and loading modes on the consolidation of the soft foundation with an upper crust with different drainage conditions were revealed. The results show that the distributions of initial excess pore pressure and additional stress with depth and loading rates have a significant influence on the consolidation process of the soft foundation with an upper crust. This influence is larger with the single-drained condition than that with the double-drained condition. Comparing the consolidation rate with a uniform initial pore pressure and additional stress, their decreasing distribution with depth quickens the consolidation at the former and middle stages. Moreover, the larger the loading rate is, the quicker the consolidation of the soft foundation with an upper crust is.


Introduction
Soft soil foundations are widely distributed in coastal and inland lake areas.Because of some natural actions or longterm engineering practices, a kind of natural or artificial hard crust with a thickness of a few meters on the soft substratum forms.The consolidation behavior of this kind of doublelayered foundation with an upper crust is quite different from that of a homogeneous foundation.
Since Gray [1] first discussed the consolidation of two contiguous layers, many research methods have been used by researchers to evaluate the effects of different factors on the consolidation behavior of a double-layered soil.In the analytical method aspect, Schiffmann and Stein [2] developed a solution of a layered foundation according to a series of complex boundary conditions and stress history.Garlanger [3] and Mesri and Rokhsar [4] analyzed the onedimensional nonlinear consolidation for a double-layered overconsolidated soil foundation.Subsequently, Zhu and Yin [5] presented a model of the double-layered foundation to account for the depth-dependent vertical stress and timedependent loading.Xie et al. [6] gained a fully explicit analytical solution for one-dimensional consolidation of twolayered soils with partially drained boundaries.Tang and Onitsuka [7] gave an analytical solution for consolidation of double-layered foundation with vertical drains under quasiequal strain condition.Xie et al. [8] derived an analytical solution for one-dimensional nonlinear consolidation of doublelayered soil considering a decreasing permeability.Miao et al. [9] proposed an analytical solution for consolidation 2 Advances in Materials Science and Engineering of a double-layered foundation partially penetrated by deep mixed columns considering one-side or two-side vertical drainage.And, then, Miao et al. [10] studied the analytical solution for one-dimensional consolidation of the doublelayered foundation with a multilevel loading condition.Zhu and Yin [11] derived several analytical solutions to analyze the consolidation of a soil layer with fairly general laws of variation of permeability and compressibility for different drainage conditions.In the finite element method aspect, Desai and Saxena [12] studied the consolidation behaviour of layered anisotropic foundations.Pyrah [13] evaluated the effects of the permeability and coefficient of volume compressibility of soil on the consolidation of layered foundations.Huang and Griffiths [14] solved the coupled (settlement and excess pore pressure), uncoupled (excess pore pressure only), and classical Terzaghi equation.In addition, Chen et al. [15] introduced a differential quadrature method to analyze the one-dimensional nonlinear consolidation of multilayered soil under partially drained boundaries and arbitrary loading condition.Fox et al. [16] built a numerical model for onedimensional large strain consolidation of layered soils with different material properties, which is able to account for time-dependent loading and boundary conditions.Kim and Mission [17] and Sadiku [18] presented solutions to the consolidation of multilayered soil based on the finite differences method.
For a double-layered foundation, if the loading, which is always time-dependent, applied on the foundation surface is large enough, such as a high embankment or a high building, its influence region inside the foundation is deep and the vertical additional stress may vary with depth.At the same time, there may exist different distributions of initial excess pore pressure for a practical project.Therefore, a consolidation model of the double-layered foundation which takes account of the depth-dependent initial excess pore pressure and additional stress, time-dependent loading needs to be built.Based on the preceding discussions, the current researchers are not able to consider the above three factors in one consolidation model.Though Zhu and Yin [5] considered the depth-dependent vertical stress and timedependent loading, their model assumed the initial excess pore pressure to be zero.They also did not analyze the effect of the depth-dependent vertical stress on the consolidation behavior of the double-layered foundation.At the same time, the numerical model built by Fox et al. [16] is able to account for time-dependent loading.However, the solution of this model should be gained using the numeric method and is an approximate solution.Therefore, this study aims at building a one-dimensional (1D) consolidation model for a doublelayered foundation to consider the depth-dependent initial pore pressure and additional stress and time-dependent loading and analyzing the effect of them on the consolidation behavior of the soft foundation with an upper crust.
The paper is organized as follows.The forthcoming section introduces the mathematical model used in this study.Subsequently, the general analytical consolidation solutions of the double-layered foundation considering the depthdependent initial excess pore pressure and additional stress, time-dependent loading under different drainage conditions

Pervious or impervious
are derived.The next section deduces the consolidation solutions of some several cases of a bilinear initial excess pore pressure and additional stress with depth under an instantaneous loading and a single-level uniform loading.After that, the average degree of consolidation defined by settlement of the special cases is verified and the effects of the depth-dependent initial excess pore pressure and additional stress as well as loading rate on the consolidation of the soft foundation with an upper crust under different drainage conditions are discussed.The final section summarizes the major findings of this study.

Mathematical Modelling
A double-layered soil profile model is presented in Figure 1.The upper layer may be a natural or an artificial hard crust.The lower layer is a native soft soil layer.There are 2 kinds of drainage conditions in this model: (1) a single-drained condition, that is, only the top of upper layer is drained; (2) a double-drained condition, that is, both top and bottom of layers are freely drained.The soil properties of the th layer are the coefficient of consolidation  V , the coefficient of permeability  V , and the compression modulus   .The compressible stratum has a total thickness of .A timedependent loading, (), is applied on the foundation surface, as shown in Figure 2.
When the time-dependent loading, (), as shown in Figure 2 is applied on the foundation surface, the resulting additional stress along the depth inside the foundation, (, ), can be expressed as follows: where  0 is the initial loading;  is duration and  is the vertical coordinate; () and () are the additional stress coefficient and initial excess pore pressure, respectively; () is the loading rate and () = /.For a 1D consolidation problem, the assumptions in Terzaghi's [19] consolidation theory are retained except for q 0 0 t c t q u q(t) the depth-dependent initial pore pressure and additional stress and time-dependent loading.Then, let (, ) = /; the partial differential equation for 1D consolidation of soils by vertical drainage is given as follows: where  is the excess pore pressure.Two commonly encountered drainage conditions are studied in this study: The former and the latter represent the drainage conditions (2) and (1) above, respectively.The initial and continuity boundary conditions are given in the following forms: where   and   () are the excess pore pressure and initial excess pore pressure in soil layer  at depth , respectively ( = 1, 2).Equation ( 4) is the initial condition.Equation ( 5) is the continuity condition of the two layers.

General Analytical Consolidation Solutions
The general consolidation solutions to consider depthdependent initial excess pore pressure and additional stress, time-dependent loading, and different drainage conditions are deduced as follows.

Single-Drained Condition.
The following dimensionless parameters were defined to simplify the expression for consolidation: where  V is the compression coefficient of the layers ( = 1, 2).
According to the derivation processes of   and (13), the coefficient   can be given as follows: where Comparing ( 21) with (22), it is found that the expressions of   and   are similar.However,   () of ( 21) and   () of ( 22) represent the additional stresses coefficient and initial excess pore pressure in soil layer  at depth  ( = 1, 2).
So, for the consolidation problem of a double-layered foundation with the pervious top surface of the upper layer, the 1D consolidation solution is gained as follows: (23)

Double-Drained Condition.
According to the similar principle, for the consolidation problem of a double-layered foundation with a double-drained condition, the consolidation solution is where   and  are similar to (15)∼( 17);   is determined by ;

Special Cases
Based on the above solutions, results for two kinds of loading conditions, (1) instantaneous loading and (2) single-level uniform loading, are given.The bilinear distributions with depth of initial excess pore pressure and additional stress are considered.

Solution with an Instantaneous
Loading.For the instantaneous loading,   = 0. So, (23)∼(24) can be simplified as follows: Assume the initial excess pore pressure has a bilinear distribution with depth; that is, where  1 ,  2 , and  3 are the initial excess pore pressure at the top of the upper layer, interface between the upper and lower layers, and bottom of the lower layer, respectively.Substituting (27) into ( 22) and (25), the expression of   is determined.
For the single-drained condition, For the double-drained condition, Then, the average excess pore pressure inside the upper and lower soils is derived as follows: So, the average consolidation degree inside the foundation defined by settlement is Advances in Materials Science and Engineering 7 where   is the compression of the two soil layers at time  and  ∞ is the final compression when the excess pore pressure is zero;   is the average initial excess pore pressure in the layers ( = 1, 2).From (31), it is evident that   only relates to  and  when the parameters of , , and  are constant.

Solution with a Single-Level Uniform
Loading.If a singlelevel uniform loading (dash line in Figure 2) is applied on the foundation surface, the following relationships can be gained: where   is the time when the loading becomes a constant value,   .
When the self-weight consolidation is not considered, the initial excess pore pressure inside the foundation is equal to zero; that is,   () = 0; then   = 0. So, the consolidation solutions of the double-layered foundation with the depthdependent additional stress are determined as follows.
When 0 <  ≤   , When  ≥   , where According to the research results by Zhang et al. [20], the distribution of the depth-dependent additional stress can be simplified to a bilinear line; that is, where  1 ,  2 , and  3 are the additional stress at the top of the upper layer, the interface between the upper and lower layers, and the bottom of the lower layer, respectively, as shown in Figure 1.Substituting (35) into ( 21) and ( 25), the consolidation solutions with a single-level loading are determined as follows.
For the single-drained condition, For the double-drained condition, (37) So, the average degree of consolidation inside doublelayered foundation defined by settlement is where   is the average coefficient of additional stress in the layers, respectively ( = 1, 2).

Model Verification. A calculation program was developed
to analyze the consolidation of the soft foundation with an upper curst in terms of (26)∼(38).In order to verify the consolidation model above, the degree of consolidation of a double-layered foundation with a uniform additional stress with depth under a single-level uniform loading was calculated as shown in Figure 3; that is,  =  = 1,  = 1,  = 5,  = 1, and   =  V1   /ℎ 2 1 = 0.2, 1, 2, 5, and 10.In Figure 3, SDC and DDC denote the single-drained condition and double-drained condition.It is obvious from Figure 3 that the consolidation rate relates to the loading rates (the values of   ) and drainage conditions, which is larger with the increasing loading rate (i.e., with a smaller value of   ) and under the double-drained condition.Moreover, the degree of consolidation of a double-layered foundation with a uniform initial excess pore pressure with depth under an instantaneous loading was also obtained, as shown in Figure 4.For this calculation,  =  = 1,   = 0,  = 1,  = 0.1, 2, 5, and 10,  = 10 for the single drainage condition, and  = 1 for the double drainage condition.It can be seen from Figure 4 that the degree of consolidation relates to both the ratios of the modulus of the upper and lower layers (the values of ) and the drainage conditions.The consolidation rate becomes smaller with the increasing , which is larger under the double-drained condition than that under the single-drained condition.In addition, compared with the analytical results of Xie et al. [21] and Xie et al. [22], it can be found that the consolidation curves under the single and double drainage conditions in Figures 3 and 4 are perfectly similar to Figures 13 and 3 from Xie et al. [21] and Figures 9 and 2 from Xie et al. [22], respectively.
In addition, Pyrah [13] discussed the influence of soil parameters on the consolidation behaviour of a doublelayered system through four idealized soil profiles using the finite element method.Zhu and Yin [5] gained the analytical solutions of the first three soil profiles and found that they are similar to the results of Pyrah [13].In order to further verify the consolidation model in this study, the degree of consolidation of the first three soil profiles (namely,  =  = 10, 0.1, 1) was calculated with a uniform additional stress ( 1 :  2 :  3 = 1: 1 : 1).The result with a depthdependent additional stress ( 1 :  2 :  3 = 1 : 0.4 : 0), which decreases along the depth, also was gained to investigate the effects of the distribution of the additional stress.For all the calculations,  = 1 and   = 0. Figure 5 presents the calculation results with the corresponding curves of Zhu and Yin [5].1 and 3 denote the consolidation curves with the uniform and depth-dependent additional stress.2 denotes the curves of (1), (2), and (3) with   = 0 in Figure 10 from Zhu and Yin [5].It is evident from Figure 5 that 1 curves are similar to 2 curves, which indicates that the results in this study agree well with the analytical solutions of Zhu and Yin [5] and the finite element results of Pyrah [13].Meanwhile, the decreasing additional stress with depth ( 1 :  2 :  3 = 1 : 0.4 : 0) quickens the consolidation process of the doublelayered foundation.For the last comparison, the degree of consolidation was calculated for the single layer defined by Case 1 in Example 3 from literature [16], as shown in Figure 6.The analytical solution was gained by the model in this study and the numerical solution was calculated by the curve with the bilinear additional stress in Figure 12 and the parameters in Table 4 from literature [16].It can be concluded from Figure 6 that the analytical solution matches the numerical results well except a minor difference at the middle stage.
Therefore, the consolidation model in this study is rational by several comparisons above with other analytical and numerical results, which can be further used to analyze the effects of distributions of initial excess pore pressure and  depth-dependent additional stress on the consolidation of the double-layered foundation.

Effect of Distributions of Initial Excess
Pore Pressure.For the soft foundation with an upper crust, the upper layer always forms due to the soil sedimentation or some longterm engineering practices; it has the same ingredient and smaller compressibility.Then,  will be smaller than 1 while  will be larger than 1.So, assume  = 0.2,  = 5, and  = 2 in the following analyses.In order to investigate the effect of the initial pore pressure on the consolidation of the double-layered foundation, three distributions of the initial pore pressure with depth were selected, namely,  1 :  2 :  3 = 1 : 1 : 1, 0 : 2 : 5, and 5 : 2 : 0, which represent the uniform, gradually increasing, and decreasing distributions with depth, respectively.
Figure 7 shows the curves of degree of consolidation under different drainage conditions.It is evident from Figure 7 that the initial excess pore pressure distributions have a similar influence on the consolidation of the doublelayered foundation under different drainage conditions.Comparing with the consolidation curve of a uniform initial excess pore pressure with depth ( 1 :  2 :  3 = 1 : 1 : 1), an increasing initial excess pore pressure with depth ( 1 :  2 :  3 = 0 : 2 : 5) slows down the consolidation process while the consolidation rate increases with a decreasing initial excess pore pressure with depth ( 1 :  2 :  3 = 5: 2 : 0).At the same time, the maximum differences of degree of consolidation between the uniform and increasing and decreasing initial excess pore pressure with depth under the single-drained condition (Curve 1 and Curve 2 in Figure 7(a)) are −9.9% and 20.8%, respectively, while the corresponding differences under the double-drained condition (Curve 1 and Curve 2 in Figure 7(b)) are −1.8% and 6.5%.Therefore, the initial excess pore pressure has a greater influence on the consolidation rate under the single-drained condition than that with the double-drained condition.For a practical project, the variation of initial excess pore pressure with depth should be considered, especially for the project with the single-drained condition.

Effect of Depth-Dependent Additional Stress.
Typically,  3 is not equal to zero and 10% of the gravity stress is always used in settlement calculations.In order to compare the effect of different additional stress, two kinds of distributions were selected, namely,  1 :  2 :  3 = 1 : 1 : 1 and 1 : 0.4 : 0. In order to investigate the effect of loading rate on the consolidation of the soft foundation with an upper crust, different values of   (i.e.,   ) were used.
The consolidation curves of different additional stress with a pervious top of the upper layer are shown in Figure 8.It is evident from Figure 8 that the additional stress distributions have a significant influence on the consolidation of the double-layered foundation.The consolidation with the decreasing additional stress with depth ( 1 :  2 :  3 = 1 : 0.4 : 0) is quicker than that with the uniform additional stress ( 1 :  2 :  3 = 1 : 1 : 1) for different   , especially at the middle stage.When   = 0.1, 0.5, 1, and 5, the maximum differences of degree of consolidation, which appear nearly at the middle stage, are 21%, 21%, 20%, and 17%, respectively.They are large enough to be a crucial factor to the success or failure of a road or building project.In addition, the differences of the degree of consolidation for different   indicate that the loading rate also affects the consolidation process.The larger the loading rate is (i.e., with a smaller value of   ), the greater the effect of the additional stress on the consolidation of the soft foundation with an upper crust is and the quicker the consolidation process is.  Figure 9 gives the consolidation curves of different additional stresses with the double-drained conditions.It is clear from Figure 9 that the additional stress also affects the consolidation of the soft foundation with an upper crust under a double-drained condition and the maximum differences of the consolidation degree are 6%, 6%, 5%, and 2% for different   .Comparing with the curves in Figure 8, it indicates that the effect of additional stress on the consolidation degree with the single-drained condition is much greater than that with the double-drained condition.At the same time, the effect of additional stress on the degree of consolidation with a double-drained condition also gradually weakens with the increase of   .

Conclusions
(1) The 1D consolidation model of the double-layered foundation proposed in this study can account for the depth-dependent initial excess pore pressure and additional stress and time-dependent loading under different drainage conditions, which can be utilized to comprehensively investigate the consolidation behavior of the double-layered foundation.The general solutions of the model under different drainage conditions were gained.
(2) The consolidation solutions of special cases, which consider two loading modes, that is, instantaneous

1 Figure 5 :Figure 6 :
Figure 5: Degree of consolidation for different soil profiles.

Figure 7 :
Figure 7: Consolidation curves with different drainage conditions.

Figure 8 :
Figure 8: Consolidation curves with a single-drained condition.