MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2018/4573780 4573780 Research Article Analytical Solution and Application for One-Dimensional Consolidation of Tailings Dam http://orcid.org/0000-0002-2287-5858 Liu Hai-ming 1 2 Nan Gan 1 Guo Wei 1 Yang Chun-he 2 Zhang Chao 2 Wang Shuo 1 Faculty of Civil Engineering and Mechanics Kunming University of Science and Technology Kunming Yunnan 650224 China kmust.edu.cn 2 State Key Laboratory of Geomechanics and Geotechnical Engineering Institute of Rock and Soil Mechanics Chinese Academy of Sciences Wuhan 430071 China cas.cn 2018 1932018 2018 10 11 2017 16 01 2018 12 02 2018 1932018 2018 Copyright © 2018 Hai-ming Liu 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.

The pore water pressure of tailings dam has a very great influence on the stability of tailings dam. Based on the assumption of one-dimensional consolidation and small strain, the partial differential equation of pore water pressure is deduced. The obtained differential equation can be simplified based on the parameters which are constants. According to the characteristics of the tailings dam, the pore water pressure of the tailings dam can be divided into the slope dam segment, dry beach segment, and artificial lake segment. The pore water pressure is obtained through solving the partial differential equation by separation variable method. On this basis, the dissipation and accumulation of pore water pressure of the upstream tailings dam are analyzed. The example of typical tailings is introduced to elaborate the applicability of the analytic solution. What is more, the application of pore water pressure in tailings dam is discussed. The research results have important scientific and engineering application value for the stability of tailings dam.

Chinese Academy of Sciences Z013009 National Key Research and Development Program of China 2017YFC0804601 National Natural Science Foundation of China 51764020 51741410 51234004 Natural Science Foundation of Yunnan Province 2015FB130
1. Introduction

Based on the assumption that the soil is isotropic and uniform, an external surface load is instantaneously applied and is held constant; a classical one-dimensional (1D) consolidation theory was proposed by Terzaghi . In order to analyze time-dependent loading, an analytical solution to the layered consolidation problem for a general set of boundary conditions and an arbitrary load history was presented by Schiffman and Stein . The 1D consolidation analytical solutions considering ramp loading were presented by Olson . A simple semianalytical method to solve the problem of 1D consolidation by taking into account the varied compressibility of soil under cyclic loadings was brought up by Cai et al. . A rigorous solution of the conventional Terzaghi 1D consolidation under haversine cyclic loading with any rest period was proposed by Müthing et al. , which is achieved using Fourier harmonic analysis for the periodic function representing the rate of imposition of excess pore water pressure. A semianalytical solution to 1D consolidation of viscoelastic unsaturated soils with a finite thickness under oedometric conditions and subjected to a sudden loading was put forward by Qin et al. . A semianalytical solution to 1D consolidation equation of fractional derivative Kelvin-Voigt viscoelastic saturated soils subjected to different time-dependent loadings was presented by Wang et al. . Under the condition of the increasing weight of superincumbent material and the length of the drainage path varies, a solution for the 1D consolidation of a clay layer whose thickness increases with time was proposed by Gibson . An exact analytical solution of the nonhomogeneous partial differential equation governing the conventional 1D consolidation under haversine repeated loading was derived and discussed by Razouki et al. [9, 10].

In order to analyze different boundary conditions, single drainage solutions for several specific variations of the permeability and shear modulus were given by Mahmoud and Deresiewicz . Several analytical solutions for the consolidation analysis of a soil layer with fairly general laws of variation of permeability and compressibility for both a single-drained condition and a double-drained condition were proposed by Zhu and Yin . Tang et al.  propose a closed-form solution for consolidation of three-layered soil with a vertical drain system and a method to solve a convergence solution by enlarging the precision for the entire thickness of the foundation and adding controlling precision for the overall average degree of consolidation of each soil layer. The solution of Terzaghi 1D consolidation equation with general boundary conditions is proposed by Mei and Chen , whose solution is validated by comparing it to the classical solution. Hawlader et al.  develop a new constitutive model for the compressibility behavior of soft clay sediments at low effective stress level and the model is used to solve finite strain 1D consolidation with pertinent initial and boundary conditions. A solution to the consolidation equation with boundary conditions that are cyclic with time is given by Rahalt and Vuez . A nonlinear theory of consolidation has been developed for an ideal normally consolidated soil by Davis and Raymond . A simple calculation procedure to analyze the one-dimensional response of saturated soil layers to pore pressure variations at the boundary described by a general time-dependent function is developed by Conte and Troncone .

The tailings dam is an important geotechnical structure in mining engineering. For a long time, the theory of reservoir dam is applied to tailings dam without any modification. However, there are many differences between tailings dam and reservoir dam, which lead to inaccurate calculation results of the pore water pressure. According to the classical Terzaghi consolidation theory, the analytical solution of the pore water pressure is discussed in this paper.

2. The General Equation of 1D Consolidation

According to Darcy’s law, the flow drag resistance of the z-direction is as follows:(1)Fz=-izγw=-nγwkvz.

The negative sign on the right side of (1) indicates that the drag resistance is opposite to the direction of flow velocity.

According to the mechanical equilibrium conditions, the flow drag resistance of the z-direction can also be expressed as follows:(2)Fz=pwz+γw.

Equations (1) and (2) can be combined as follows:(3)pwz+γw+nγwkvz=0.

It is assumed that the water in porosity is incompressible under the consolidation process. Thus, γw is a constant. During the process of soil consolidation, because the porosity of soil particles is continuously compressed, the permeability coefficient of soil is decreased continuously and the permeability coefficient k of soil is changed with depth, that is, k=f(z). Taking the derivative of (3), one can have(4)2pwz2+nγwzvzk=0.

According to fractional derivative rule, (4) can be expanded as follows:(5)2pwz2+nγw1kvzz-nγwvz1k2dkdz=0.

Assume that the horizontal direction consolidation of the soil can be neglected; that is, the consolidation is under confined compression condition. According to the 1D consolidation condition, (6) can be deduced(6)nvzz=εzt.

Under the condition of 1D consolidation, the axial strain εv is equal to the volume strain εv of the soil. According to Pane and Schiffman  research, (4) and (5) can be rewritten as follows:(7)2pwz2+γwkεzt=0,(8)2pwz2+γwkεzt-nγwvz1k2dkdz=0.

In the consolidation problem, it is of great significance to study the change of pore water pressure with time and position under external load. Therefore, (8) can be expanded according to pore water pressure.

Assuming that the soil profile is shown as Figure 1, the pore water pressure of a point located z under load σ0 and Δσ0 is u. According to the Terzaghi effective stress principle, the effective stress σ could be given as follows under consolidation:(9)σ=σ-pw=σ-γwh+z+u.

Schematic diagram of stress in soil.

According to Figure 1, the total stress σ of a point located z is as follows:(10)σ=σ0+Δσ+γsz+γwh.

Substituting (10) into (9) yields(11)σ=σ0+Δσ+γz-u.

Take the derivative of (11) with respect to t and consider the relationship z/t equals H/t(12)σt=tΔσ+γHt-ut.

According to the compression curve of 1D consolidation test, the volume compressibility factor mv can be defined as follows:(13)mv=11+eeσ,(14)εzt=11+eet.

Substituting (13) into (14) gives(15)εzt=mvσt.

Substituting (15) into (8), it is noted that the last term of (8) has the following relation:(16)nγwvz1k2dkdz=-Fz1kdkdz=-pwz+γw1kdkdz.

Equation (8) can be derived as(17)2pwz2+mvγwkσt+pwz+γw1kdkdz=0.

Substituting (15) into (17), (17) can be expressed as pore water pressure. Equation (18) can be obtained(18)2uz2+mvγwkΔσt+γHt-ut+1kdkdzuz=0.

Equation (18) is a general equation which reflects the 1D consolidation process. It takes into account the situation where the applied load changes with time, the soil thickness changes with time, the soil permeability changes with depth, and so on.

3. Accumulation and Dissipation of Pore Water Pressure in Tailings Dam

Tailings are mine wastes produced in the mining engineering, which are sent to tailings reservoir by pipe or flume. The tailings dam is an important part in mining engineering, which is consisted by the initial dam and fill dam. The initial dam is made by permeable rockfill generally, and the fill dam is formed by tailings. In general, the construction of a tailing dam takes many decades or even a century. According to construction method, the tailings dam can be divided into upstream tailings dam, downstream tailings dam, midline method tailings dam, and so on. The number of tailings ponds in China has reached more than 12,000 by statistics. Because the downstream of the tailings dam is residents living area or mining production area generally the social and people’s property damage a huge impact if the tailings dam is failed.

Because the upstream tailings dam has advantages of simple operation, low construction costs, the less need for coarse particles, and so on, according to statistics, 95% of the tailings dams are adopting the construction of the upstream tailings dam in China. On the other hand, the construction process of the upstream tailings dam cannot precisely control the shape of the tailings dam. The deposition structure of the tailings dam is very complex. The upstream tailings have shortcoming of long infiltration distant and poor stability. In view of the tailings slurry discharge and tailings particle deposition following the sediment mechanics, the section profile of the tailings dam has obvious regularity. The coarser the particle size is, the shorter the average distance migrates and vice versa. Generally, the deposition order of the tailings dam along the dry beach face is tailings fine sand, tailings silt, tailings sand, and tailings mud.

If particle size distribution of the tailings and the length of the dry beach face are kept constant, the interface of the tailings material should be substantially parallel to the surface of the dam slope, which is shown in Figure 2. For most of the tailings dams, because the tailings water contains more heavy metals, the initial topography and the side of tailings dam are treated as an impermeable boundary. On the other hand, the tailings dam extends very long in the direction of the reservoir area. Therefore, the consolidation of the tailings dam is a 1D problem.

Schematic diagram of the upstream tailings dam.

As can be seen from Figure 2, the closer to inside the reservoir area is, the thicker the tailings mud is and the worse the average consolidation is. In other words, the closer to inside the reservoir area is, the greater pore water pressure of the tailings dam is and the lower shear strength of the tailings material is. Generally speaking, the construction of a tailings dam needs several decades to reach the design elevation. With the continuous production of mining, the height of the tailings dam is increasing. During this process, the accumulation and dissipation of the excess pore water pressure of the tailings material experience several different stages of development.

3.1. Tailings Thickness Increases with Time

In this stage, the thickness of the tailings increases with time. The tailings material produces great compressive deformation under self-weight effect. The pore water pressure produced by the self-weight pressure at the early stage is partially dissipated. Nevertheless, the pore water pressure caused by self-weight pressure is increased with the tailings thickness at the later stage, which leads to the increase of the pore water pressure with time. Therefore, the net pore water pressure is accumulated during this process.

3.2. Load on the Tailings Layer Increases with Time

The self-deposition of the tailings mud continues to extend to interior of the tailings reservoir area in this stage. The self-deposition of the tailings layer has ended where the end of the tailings reservoir is near. It overlapped with particles coarse tailings and tailings mud. The tailings mud has to bear a growing load with time. The tailings mud continues to consolidate under the combined action of self-weight and additional load. At this stage, the pore pressure produced by the self-weight pressure dissipates a little part, and the more pore water pressure produced by additional load accumulates. The net pore water pressure tends to accumulate in this process.

For a specific part of the tailings dam, the upper boundary of the tailings dam has reached the designed elevation at a certain moment. Therefore, the additional load will remain unchanged from this moment. During this stage, the tailings mud will continue to consolidate under the combined action of self-weight and constant additional load. Thus, the pore water pressure will dissipate within the tailings dam. Obviously, the different parts of tailings mud have different consolidation stages. There is no additional load at the thickest tailings mud in the tailings reservoir until it is closed. Therefore, the tailings mud is always in the first stage and will never enter the second and third stages. From the above qualitative analysis, the pore water pressure, which is located at the slope dam segment, reaches the maximum value when the additional load just stops growing. The pore water pressure, which is located at dry beach segment and artificial lake segment, achieves the maximum value when the tailings dam reaches the maximum height. The numerical values of the pore water pressure at different locations and moments can be obtained by solving the partial differential equations of consolidation problems.

4. Analytical Solutions of the Pore Water Pressure

It is necessary to accurately calculate the accumulation and dissipation of the pore water pressure in tailings dam, of which a solution of 1D consolidation problem can be simplified. During the consolidation process, the soil parameters of tailings material (such as the change of bulk density, permeability coefficient, and consolidation coefficient) are changed with the consolidation process. The change law can be determined through many experiments. Under normal circumstances, it is difficult to obtain the analytical solution when the change law of soil parameters is considered. Therefore, the exact solution of the problem can only be depended on numerical calculation method.

As the horizontal length of the tailings reservoir is usually far greater than the vertical thickness (i.e., the horizontal length is 10 times more than the vertical thickness), the drainage consolidation effect of the horizontal direction can be ignored. Therefore, the vertical direction of drainage is needed to consider. Then, it is only a 1D consolidation problem. Based on the assumption of small-strain and constant of soil parameters, such as mν, k, Δσ/t, γ, and H/t are constant, (18) is degenerated into a constant coefficient nonhomogeneous parabolic partial differential equation(19)2uz2+1cνα+γQ-ut=0,where(20)α=Δσt,Q=Ht,cν=kmνγw.

The above three different stages are specifically discussed as follows.

4.1. The Thickness Mud of Tailings Mud Increases with Time

The overburden load of the tailings mud at this stage is zero, that is, α=0. The consolidation equation yields(21)2uz2+1cνγQ-ut=0.

If the bottom of the tailings mud is impermeable, the coordinates origin is taken as the impermeable bottom. The boundary conditions of the problem can be given as follows:(22)z=0,uz=0,(23)z=H,u=0.

The initial condition of the problem is given by(24)t=0,u=0.

According to Gibson’s  study, (25) is chosen as the solution of (21)(25)u=γh-t-1/20gξe-x-ξ2/4cνt+e-x+ξ2/4cνtdξ.

It can be verified that (25) satisfies the boundary condition (22). The choice of function g(ξ) needs to satisfy (23). Substituting (25) into (23) gives(26)12γht1/2eh2/4cνt=0gξe-x+ξ2/4cνtcoshhξ2cvtdξ.

For any function ht, the partial differential equation (26) has no numerical solution. With regard to ht=mt, the following parameters can be introduced to transform the equation(27)ξ2=τ,p=14cvt,Fτ=gττ-1/2coshkτ1/22cv.

Equation (26) can be transformed into(28)18mγcv-3/2p-3/2em2/16cv2p=0Fτe-pτdτ.

Equation (28) can be obtained through the Laplace change:(29)Fτ=γ2π1/2cv1/2sinhmτ1/2cv.

Thence(30)gξ=γξ2π1/2cv1/2tanhmξcv.

Combining the solutions of (25) and (30), the solution of (21) is obtained(31)u=γQt-γπcνt-1/2exp-z24cνt0ξtanhQξ2cνcoshzξ2cνtexp-ξ22cνtdξ.

It can be verified that (31) satisfies the initial condition of (24).

During this stage, the thickness of tailings mud not only does not increase with time but also gradually decreases with the increase of effective stress. For the sake of simplicity, assuming that the thickness of tailings mud H is constant, the basic consolidation differential (19) can be simplified as(32)cν2uz2=ut-Δσt.

The initial condition of the problem is given by(33)t=t1,u=u0z.

The boundary conditions of the problem can be given as(34)z=0,u=0z=H,uz=0.

The following three questions are called problem A, problem B, and problem C, respectively:

A :(35)ut=cν2uz2Initialconditions:t=t1,u=u0zBoundaryconditions:z=0,u=0z=H,uz=0.

B : (36)ut=cν2uz2+ΔσtInitialconditions:t=t1,u=0Boundaryconditions:z=0,u=0z=H,uz=0.

C : (37)ut=cν2uz2+ΔσtInitialconditions:t=t1,u0zBoundaryconditions:z=0,u=0z=H,uz=0.

If the solution of problem A, problem B, and problem C is uA, uB, and uC respectively, it can be directly proved through substitution method:(38)uC=uA+uB.

The solution of uA belongs to the Terzaghi classical 1D consolidation problem. The analytical solution can be obtained.(39)uA=j=1Ansinjπz2Hexp-j2π2cνt-t14H2,where(40)An=1H02Hu0zsinjπz2Hdz.

u 0 z is the solution of (31). Substituting (31) into (40) by Matlab software, An can be obtained. Substituting it into (39), the numerical solution uA can be given.

For solution to the problem uB, Schiffman  acquires the following series solutions:(41)uB=16αH2π3cν1,3,51j3sinjπz2H1-exp-j2π2cνt-t14H2.

When the tailings mud is below the groundwater level, the total stress σz of the tailings mud at the depth z can be expressed as follows:(42)σz=γz+Δσtt-t1.

The ratio of pore water pressure to total stress can be easily determined as follows:(43)Uz=uA+uBσz.

In this stage, the basic differential equation and the definite condition are problem D:

D : (44)ut=cν2uz2t=t2,u=uA+uBt=t2z=0,u=0z=H,uz=0.

The following questions can be called problem E:

E :(45)ut=cν2uz2t=t2,u=uBt=t2=u1zz=0,u=0z=H,uz=0.

If the solution of problem A is uA and the solution of problem E is uE, the solution of problem D can be easily demonstrated is uD as follows:(46)uD=uA+uE.

The definite condition of problem E and problem A is exactly the same in form, except that the function u1(z) is different from the function u0(z).

Suppose that(47)uE=j=1Bnsinjπz2Hexp-j2π2cvt-t24H2.

Then(48)Bn=1H02Hu1zsinjπz2Hdz,where(49)u1z=16αH2π3cv1,2,31j3sinjπz2H1-exp-j2π2cvt2-t14H2.

According to the orthogonal rule of the solution, the following relationship can be obtained:(50)1H02Hsinmπz2Hsinjπz2Hdz=0,mj1,m=j.

Equation (48) can be simplified through (50)(51)Bn=16αH2π3cv1j31-exp-j2π2cvt2-t14H2j=1,3,50j=2,4,6.

Substituting (51) into (47) yields(52)uE=16αH2π3cv1,3,51j31-exp-j2π2cvt2-t14H2sinjπz2Hexp-j2π2cvt-t24H2.

For the third stage, the analytical solution of pore water pressure is obtained.

Similar to (42), the ratio of pore water pressure to total stress at a certain depth of tailings mud at this stage is as follows:(53)Uz=uA+uEσz.

During this stage, σ(z) does not change with time, and it can take the following value:(54)σz=γz+αt2-t1.

5. Example

The simplified calculation section of the representative section of a tailings dam is shown in Figure 3. The length of the tailings dam is 368 m. Therefore, the problem can be simplified as 1D consolidation. The slime layer below the level is assumed for impervious boundary. The rising speed of the dam height is 1.75 m/y. The saturation density of gravel is 22.0 kN/m3, the saturation density of the tailings fine sand is 19.6 kN/m3, the saturation density of the tailings sand is 19.1 kN/m3, and the saturation density of the tailings mud is 18.8 kN/m3. The permeability coefficient of gravel is 3 × 10−3 m/s, the permeability coefficient of the tailings fine sand is 5 × 10-6 m/s, the permeability coefficient of the tailings sand is 2 × 10-6 m/s, and the permeability coefficient of the tailings mud is 1.5 × 10-8 m/s. The consolidation coefficient of gravel is 3.8 × 10-1 cm2/s, the consolidation coefficient of the tailings fine sand is 5.4 × 10-3 cm2/s, the consolidation coefficient of the tailings sand is 2.8 × 10-3 cm2/s, and the consolidation coefficient of the tailings mud is 1.2 × 10-5 cm2/s. Please determine the distribution of pore water pressure when the tailings dam reaches maximum height.

Simplified section for calculation of pore water pressure in a tailings dam.

The above problem can be solved through (31), (46), and (52). The calculation can be accomplished through the MATLAB software. Firstly, a coordinate system can be set up. Then, the corresponding coordinates are taken out of the interval 1 m. The region is determined according to the coordinates. Then, the pore water pressure can be calculated through the corresponding equation. Finally, the contour map can be drawn through a set of array. The results of the calculation are shown in Figure 4.

Distribution of pore water pressure at the moment when the dam rises to its maximum height.

6. Discussion

Due to the similarity between the tailings dam and the reservoir dam in the geotechnical structure, many scholars directly introduce the calculation methods of the reservoir dam, whose theory is relatively mature, into the tailings dam. Considering the difference on construction cycle, construction materials between the tailings dam and reservoir dam, it makes the calculation results inconsistent with the actual.

Based on the Terzaghi consolidation theory of 1D, the tailings dam is divided into the slope dam segment, the dry beach segment, and the artificial lake segment. The solutions of the pore water pressure are derived, respectively. The analysis shows that the additional load of the slope dam segment is unchanged, which can be calculated using (46). The analysis indicates that tailings thickness of the artificial lake segment increases with time during construction period, which can be calculated by (31). The situation of the dry beach section is slightly complicated. The additional load of the dry beach section which is under the slope dam segment is basically the same, which can be calculated by (52). The additional load on the other sections of the dry beach segment increases with time, which can be calculated using (31).

The theoretical derivation is based on 1D consolidation theory. It is only considering consolidation in the vertical direction. Since the horizontal scale of most tailings dams is much larger than the vertical direction, it has little influence to ignore the drainage of horizontal. From the point of engineering view, it is conservative to the stability of the tailings dam. Considering drainage of the horizontal, the pore water pressure will be lessened and the safety factor of the tailings dam will be greater.

It is assumed that the deformation of the tailings is small deformation during the consolidation process. If the actual tailings are loose relatively, the deformation of the tailings is large deformation. He et al.  studied the pore water pressure of saturated soils using the updated Lagrangian formulation of large strain method (ULM), total Lagrangian formulation of large strain method (TLM), and small strain method (SSM). The results point out that the pore water pressure of the above three methods is exactly the same when the strain is less than 2.4%. When the strain is greater than 2.4%, the pore water pressure of the ULM and SSM still keeps the same, and the pore water pressure of the TLM is obviously less than the ULM and SSM. When the strain reaches 7.5%, the pore water pressure of the ULM and SMM began to bifurcate. When the deformation of the tailings dam is examined for a relatively short time, the consolidation of prophase tailings dam has been completed under previous load. The consolidation deformation of the whole tailings dam is relatively small, and its deformation conforms to the small deformation assumption. On the other hand, it is difficult to calculate and solve partial differential equations using large deformation theory. The differences between large deformation and small deformation need to be further studied.

It is supposed that the mechanical parameters such as permeability coefficient and consolidation coefficient are constant during the consolidation process. Previous studies  have shown that the permeability coefficient of rock and soil decreases with the increase of deformation during the consolidation process. It is believed that the value of Cc/Ck decides whether it is necessary to take under consideration the effect of nonlinear property by Zhuang et al. . The results show that pore water pressure calculated depends on the ratio of loading intensity to the initial effective vertical stress subjected to time-dependent loading by Conte and Troncone . How the permeability coefficient of tailings material changes during the consolidation process is less studied by relevant scholars. How the permeability coefficient of tailings material changes during the consolidation process is less studied by relevant scholars. From the aspect of engineering view, the larger the pore water pressure of actual is, the smaller the safety factor of tailings dam is. Therefore, it is dangerous for the tailings dam. The law of permeability coefficient of tailing material during consolidation process needs to be further studied.

7. Conclusion

Tailings dam is a very important geotechnical structure of mine engineering. The calculation of pore water pressure has a great impact on the safety factor of tailings dam slope. How to accurately estimate pore water pressure is very difficult. Based on the assumption of 1D consolidation and small strain of tailings material, a general equation of the pore water pressure is proposed. According to dissipation and accumulation characteristics of the pore water pressure in the tailings dam, the tailings dam can be divided into the slope dam segment, the dry beach segment, and the artificial lake segment. The analytic solutions of the corresponding segment are obtained through solving the partial differential equation, which has great significance to the stability of the tailings dam.

Notations τ :

Shear strength of the tailings material

σ z :

Maximum increases in vertical total stress as a function of depth z

ξ :

A parameter introduced to transform an equation

α :

x :

Coordinate of x direction

p :

A parameter introduced to transform an equation

n :

Porosity of the tailings material

m , j :

Counters, 1, 2, 3,…

g ( τ ) :

A function of variable τ

g ( ξ ) :

A function introduced into solving partial differential equation

Q :

The rate of the thickness of the tailings mud, Q=H/t

F τ :

A function introduced to transform an equation

σ :

Stress of the tailings material

z :

Coordinate of z direction

u :

The pore water pressure

k :

Coefficient of coefficient

h :

The height of static water level

H :

Thickness of tailings mud

e :

Porosity ratio

Δ σ :

Increment of stress

u 1 z :

Termination pore water pressure as a function of depth z

u 0 z :

Initial pore water pressure as a function of depth z

t 2 :

Termination time

t 1 :

Initial time

c ν :

Coefficient of consolidation

U z :

Saturation

B n :

Coefficients to be determined

A n :

Coefficients to be determined

σ 0 :

Initial total stress

σ :

Effective vertical stress

ε z :

Strain in the z direction

γ w :

Bulk density of water

γ s :

Saturated bulk density of tailings dam

γ :

Buoyancy unit weight of tailings material

v z :

The actual velocity along the flow direction in the tailings dam

p w :

Water pressure

m v :

Coefficient of volume compressibility

i z :

Hydraulic gradient in the z direction

H 1 :

The total height of tailings mud plus static water level

F z :

The drag resistance force on the pore wall of a unit volume in the z direction.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This paper is supported by the Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, (Grant no. Z013009), the National Key Research and Development Program of China (Project no. 2017YFC0804601), the National Natural Science Foundation of China (Grant nos. 51764020; 51741410; 51234004), and the Natural Science Foundation of Yunnan Province (Grant no. 2015FB130). The authors would like to thank them for providing the financial support for conducting this research.

Terzaghi K. Theoretical Soil Mechanics 1943 New York, NY, USA John Wiley & Sons 10.1002/9780470172766 Schiffman R. L. Stein J. R. One-dimensional consolidation of layered systems Journal of the Soil Mechanics and Foundations Division 1970 96 4 1499 1504 Olson R. E. Consolidation under time-dependent loading Journal of the Geotechnical Engineering Division 1977 103 1 55 60 Cai Y.-Q. Geng X.-Y. Xu C.-J. Solution of one-dimensional finite-strain consolidation of soil with variable compressibility under cyclic loadings Computers & Geosciences 2007 34 1 31 40 2-s2.0-33845485857 10.1016/j.compgeo.2006.08.008 Müthing N. Razouki S. S. Datcheva M. Schanz T. Rigorous solution for 1-D consolidation of a clay layer under haversine cyclic loading with rest period SpringerPlus 2016 5 1, article no. 1987 2-s2.0-84995920160 10.1186/s40064-016-3660-9 Qin A. Sun D. Zhang J. Semi-analytical solution to one-dimensional consolidation for viscoelastic unsaturated soils Computers & Geosciences 2014 62 110 117 2-s2.0-84905393723 10.1016/j.compgeo.2014.06.014 Wang L. Sun D. Li P. Xie Y. Semi-analytical solution for one-dimensional consolidation of fractional derivative viscoelastic saturated soils Computers & Geosciences 2017 83 30 39 2-s2.0-84994012352 10.1016/j.compgeo.2016.10.020 Gibson R. E. The Progress of Consolidation in a Clay Layer Increasing in Thickness with Time Géotechnique 1958 8 4 171 182 10.1680/geot.1958.8.4.171 Razouki S. S. Schanz T. One-dimensional consolidation under haversine repeated loading with rest period Acta Geotechnica 2011 6 1 13 20 2-s2.0-79952694032 10.1007/s11440-010-0132-1 Razouki S. S. Bonnier P. Datcheva M. Schanz T. Analytical solution for 1D consolidation under haversine cyclic loading International Journal for Numerical and Analytical Methods in Geomechanics 2013 37 14 2367 2372 10.1002/nag.2188 2-s2.0-84883751636 Mahmoud M. S. Deresiewicz H. Settlement of inhomogeneous consolidating soils—I: The single‐drained layer under confined compression International Journal for Numerical and Analytical Methods in Geomechanics 1980 4 1 57 72 2-s2.0-0018911579 10.1002/nag.1610040105 Zbl0425.73095 Zhu G. Yin J.-H. Analysis and mathematical solutions for consolidation of a soil layer with depth-dependent parameters under confined compression International Journal of Geomechanics 2012 12 4 451 461 2-s2.0-84879527671 10.1061/(ASCE)GM.1943-5622.0000152 Tang X. W. Niu B. Cheng G. C. Shen H. Closed-form solution for consolidation of three-layer soil with a vertical drain system Geotextiles and Geomembranes 2013 36 81 91 10.1016/j.geotexmem.2012.12.002 2-s2.0-84872372321 Mei G.-X. Chen Q.-M. Solution of Terzaghi one-dimensional consolidation equation with general boundary conditions Journal of Central South University 2013 20 8 2239 2244 2-s2.0-84881348008 10.1007/s11771-013-1730-5 Hawlader B. C. Muhunthan B. Imai G. State-dependent constitutive model and numerical solution of self-weight consolidation Géotechnique 2008 58 2 133 141 2-s2.0-41349100271 10.1680/geot.2008.58.2.133 Rahalt M. A. Vuez A. R. Analysis of settlement and pore pressure induced by cyclic loading of silo Journal of Geotechnical and Geoenvironmental Engineering 1998 124 12 1208 1210 10.1061/(ASCE)1090-0241(1998)124:12(1208) 2-s2.0-0032459539 Davis E. H. Raymond G. P. A non-linear theory of consolidation Géotechnique 1965 15 2 161 173 10.1680/geot.1965.15.2.161 Conte E. Troncone A. Soil layer response to pore pressure variations at the boundary Géotechnique 2008 58 1 37 44 10.1680/geot.2008.58.1.37 2-s2.0-38649127568 Pane V. Schiffman R. L. A note on sedimentation and consolidation Géotechnique 1985 35 1 69 72 10.1680/geot.1985.35.1.69 Gibson R. E. A heat conduction problem involving a specified moving boundary Quarterly of Applied Mathematics 1959 16 4 426 430 10.1090/qam/103713 MR0103713 Schiffman R. L. Consolidation of soil under time-dependent loading and varying permeability 37 Proceedings of the Thirty-Seventh Annual Meeting of the Highway Research Board 1958 584 617 He K. S. Shen Z. J. Peng X. X. The comparison of large strain method using total and updated Lagrangian finite element formulation and small strain method Chinese Journal of Geotechnical Engineering 2000 22 1 30 34 Mesri G. Olson R. E. Mechanisms controlling the permeability of clays Clays and Clay Minerals 1971 19 3 151 158 2-s2.0-0001473666 10.1346/CCMN.1971.0190303 Zhuang Y.-C. Xie K.-H. Li X.-B. Nonlinear analysis of consolidation with variable compressibility and permeability Journal of Zhejiang University (Engineering Science) 2005 6 3 181 187 10.1631/jzus.2005.A0181 2-s2.0-14744286640 Conte E. Troncone A. Nonlinear consolidation of thin layers subjected to time-dependent loading Canadian Geotechnical Journal 2007 44 6 717 725 2-s2.0-34848814512 10.1139/T07-015 10.1139/t07-015