Exact solutions for one‐dimensional consolidation of single‐layer unsaturated soil

Based on the Fredlund consolidation theory of unsaturated soil, exact solutions of the governing equations for one‐dimensional consolidation of single‐layer unsaturated soil are presented, in which the water permeability and air transmission are assumed to be constants. The general solution of two coupled homogeneous governing equations is first obtained. This general solution is expressed in terms of two functions psi₁ and ψ₂, where ψ₁ and ψ₂, respectively, satisfy two second‐order partial differential equations, which are in the same form. Using the method of separation of variables, the two partial differential equations are solved and exact solutions for three typical homogeneous boundary conditions are obtained. To obtain exact solutions of nonhomogeneous governing equations with three typical nonhomogeneous boundary conditions, the nonhomogeneous boundary conditions are first transformed into homogeneous boundary conditions. Then according to the method of undetermined coefficients and exact solutions of homogenous governing equations, the series form exact solutions are put forward. The validity of the proposed exact solutions is verified against other analytical solutions in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Analytical solution for one‐dimensional consolidation of unsaturated soils using eigenfunction expansion method

This paper introduces an exact analytical solution for governing flow equations for one‐dimensional consolidation in unsaturated soil stratum using the techniques of eigenfunction expansion and Laplace transformation. The homogeneous boundary conditions adopted in this study are as follows: (i) a one‐way drainage system of homogenous soils, in which the top surface is considered as permeable to air and water, whereas the base is an impervious bedrock; and (ii) a two‐way drainage system where both soil ends allow free dissipation of pore‐air and pore‐water pressures. In addition, the analytical development adopts initial conditions capturing both uniform and linear distributions of the initial excess pore pressures within the soil stratum. Eigenfunctions and eigenvalues are parts of the general solution and can be obtained based on the proposed boundary conditions. Besides, the Laplace transform method is adopted to solve the first‐order differential equations. Once equations with transformed domain are all obtained, the final solutions, which are proposed to be functions of time and depth, can be achieved by taking an inverse Laplace transform. To verify the proposed solution, two worked examples are provided to present the consolidation characteristics of unsaturated soils based on the proposed method. The validation of the recent results against other existing analytical solutions is graphically demonstrated. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical Laplace-Fourier transform inversion technique for layered-soil consolidation problems: I. Fundamental solutions and validation

A numericl method for solving consolidation problems of layered soils is developed. Starting from the governing differential equations for the coupled poro-elastic medium, the governing partial differential equations are reduced to ordinary differential equations by means of the appropriate displacement functions and Laplace-Fourier transformation. Once the fundamental solution in the transformed domain has been found, the solution in the physical domain is obtained by numerically inverting the transformations. A series of soil consolidation problems have been solved and validated against existing solutions in order to compare the feasibility and the accuracy of the present technique.     

A closed form analytical solution for two‐dimensional plane strain consolidation of unsaturated soil stratum

This paper discusses the excess pore‐air and pore‐water pressure dissipations and the average degree of consolidation in the 2D plane strain consolidation of an unsaturated soil stratum using eigenfunction expansion and Laplace transformation techniques. In this study, the application of a constant external loading on a soil surface is assumed to immediately generate uniformly or linearly distributed initial excess pore pressures. The general solutions consisting of eigenfunctions and eigenvalues are first proposed. The Laplace transform is then applied to convert the time variable t in partial differential equations into the Laplace complex argument s. Once the domain is obtained, a simplified set of equations with variable s can be achieved. The final analytical solutions can be computed by taking a Laplace inverse. The proposed equations predict the two‐dimensional consolidation behaviour of an unsaturated soil stratum capturing the uniformly and linearly distributed initial excess pore pressures. This study investigates the effects of isotropic and anisotropic permeability conditions on variations of excess pore pressures and the average degree of consolidation. Additionally, isochrones of excess pore pressures along vertical and horizontal directions are presented. It is found that the initial distribution of pore pressures, varying with depth, results in considerable effects on the pore‐water pressure dissipation rate whilst it has insignificant effects on the excess pore‐air pressure dissipation rate. Copyright © 2015 John Wiley & Sons, Ltd.     

A quasistatic analysis of a plate on consolidating layered soils by analytical layer‐element/finite element method coupling

In this paper, a coupling method between finite element and analytical layer‐elements is utilized to analyze the time‐dependent behavior of a plate of any shape and finite rigidity resting on layered saturated soils. Based on the integral transform techniques together with the aid of an order reduction method, an analytical layer‐element solution is derived from the governing equations for three‐dimensional Biot consolidation with respect to a Cartesian coordinate system and then extended to be the fundamental solution for the layered saturated soil under a point load. The Mindlin plate is modeled by eight‐noded isoparametric elements. The governing equations of the interaction between soil and plate in the Laplace‐Fourier transformed domain are deduced by referring to the coupling theory of FEM/BEM, and the final solution is obtained by applying numerical inversion. Numerical examples concerned with the time‐dependent response of a plate are performed to demonstrate the influence of soil and plate properties on the interaction process. Copyright © 2014 John Wiley & Sons, Ltd.     

竖向简谐荷载下二维层状饱和地基的解析层元解

基于Biot固结理论,运用解析层元方法求解竖向简谐荷载作用下二维层状饱和地基的动力响应问题。从直角坐标平面应变问题控制方程出发,通过Fourier-Laplace变换将偏微分方程组转化为常微分方程组,求解得到单层饱和地基的解析层元。结合层间连续条件和边界条件,组装得到多层饱和地基的总刚度矩阵方程,进而求得变换域内的解。借助Fourier-Laplace逆变换的数值积分方法,获得平面应变动力问题在物理域内的解,编制了相应的计算程序,其计算结果与已有文献结果吻合较好。通过算例分析了荷载圆频率、荷载作用深度及地基成层性对地基竖向位移的影响。计算结果表明：随荷载圆频率的增大,地基竖向位移先增加后减小;地基竖向位移在荷载作用点处呈现波峰,且受表层土性的影响较大。     

Response of a tunnel deeply embedded in a viscoelastic medium

This paper presents a three‐dimensional energy‐based solution for the time‐dependent response of a deeply embedded and unsupported semi‐infinite tunnel of circular cross‐section. The tunnel is taken to be excavated quasi‐instantaneously from an infinite rock body that initially exhibits an isotropic stress state and that is made up of a homogeneous, isotropic and viscoelastic material. The viscoelastic behaviour is modelled by means of Burger's model, and the rock is taken to behave volumetrically linear elastic and to exhibit exclusively deviatoric creep. This viscoelastic problem is transformed into the Laplace domain, where it represents a quasi‐elastic problem. The displacement fields in the new solution are taken to be the products of independent functions that vary in the radial and longitudinal directions. The differential equations governing the displacements of the system and appropriate boundary conditions are obtained using the principle of minimum potential energy. The solutions for these governing equations in the Laplace domain are then obtained analytically and numerically using a one‐dimensional finite difference technique. The results are then transformed back into the time domain using an efficient numerical scheme. The accuracy of the new solution is comparable with that of a finite element analysis but requires much less computation effort. Copyright © 2011 John Wiley & Sons, Ltd.     

Analytical solutions of a finite two‐dimensional fluid‐saturated poroelastic medium with compressible constituents

An analytical solution is proposed for transient flow and deformation coupling of a fluid‐saturated poroelastic medium within a finite two‐dimensional (2‐D) rectangular domain. In this study, the porous medium is assumed to be isotropic, homogeneous, and compressible. In addition, the point sink can be located at an arbitrary position in the porous medium. The fluid–solid interaction in porous media is governed by the general Biot's consolidation theory. The method of integral transforms is applied in the analytical formulation of closed‐form solutions. The proposed analytical solution is then verified against both exact and numerical results. The analytical solution is first simplified and validated by comparison with an existing exact solution for the uncoupled problem. Then, a case study for pumping from a confined aquifer is performed. The consistency between the numerical solution and the analytical solution confirms the accuracy and reliability of the analytical solution presented in this paper. The proposed analytical solution can help us to obtain in‐depth insights into time‐dependent mechanical behavior due to fluid withdrawal within finite 2‐D porous media. Moreover, it can also be of great significance to calibrate numerical solutions in plane strain poroelasticity and to formulate relevant industry norms and standards. Copyright © 2014 John Wiley & Sons, Ltd.     

多层横观各向同性地基轴对称固结的传递矩阵解

从横观各向同性地基轴对称Biot固结的基本方程出发,通过关于t的Laplace变换和关于r的Hankel变换,得到关于z的一阶常微分方程组。然后,对变换域内的基本未知量进行线性化处理,建立了变换域内的基本状态变量在z = 0处和任意深度处z的显式关系。利用传递矩阵法,结合层间连续性条件和边界条件,得到了多层横观各向同性地基的Biot固结轴对称问题的解答。该解答能避免随着层数增加而需要求解大型方程组的困难,明显地提高了计算效率。     

Time‐dependent behaviour of a rigid foundation on a transversely isotropic soil layer

An analytical solution is presented in this paper to study the time‐dependent settlement behaviour of a rigid foundation resting on a transversely isotropic saturated soil layer. The governing equations for a transversely isotropic saturated soil, within Biot's poroelasticity framework, are solved by means of Laplace and Hankel transforms. The problem is subsequently formulated in the Laplace transform domain in terms of a set of dual integral equations that are further reduced to a Fredholm integral equation of the second kind and solved numerically. The developed analytical solution is validated via comparison with the existing analytical solution for an isotropic saturated soil case, and adopted as a benchmark to examine the sensitivities of the mesh refinement and the locations of truncation boundaries in the finite element simulations using ABAQUS. Particular attention is paid to the influences of the degree of soil anisotropy, boundary drainage condition, and the soil layer thickness on the consolidation settlement and contact stress of the rigid foundation. Copyright © 2009 John Wiley & Sons, Ltd.     

Consolidation analysis of saturated multi‐layered soils with anisotropic permeability caused by a point sink

This paper presents the analytical layer‐element method to analyze the consolidation of saturated multi‐layered soils caused by a point sink by considering the anisotropy of permeability. Starting from the governing equations of the problem, the solutions of displacements and stresses for a single soil layer are obtained in the Laplace–Hankel transformed domain. Then, the analytical layer‐element method is utilized to further derive the solutions for the saturated multi‐layered soils in the transformed domain by combining with the boundary conditions of the soil system and continuity conditions between adjacent layers. The actual solutions in the physical domain can be acquired by the inversion of Laplace–Hankel transform. Numerical results are carried out to show the accuracy and stability of the proposed method and evaluate the influence of sink depth and anisotropic permeability on excess pore pressure and surface settlement. Copyright © 2011 John Wiley & Sons, Ltd.     

Semi‐analytical solution to one‐dimensional consolidation for unsaturated soils with semi‐permeable drainage boundary under time‐dependent loading

This paper presents semi‐analytical solutions to Fredlund and Hasan's one‐dimensional consolidation of unsaturated soils with semi‐permeable drainage boundary under time‐dependent loadings. Two variables are introduced to transform two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved by the Laplace transform. The pore‐water pressure, pore‐air pressure and settlement are obtained in the Laplace domain. Crump's method is adopted to perform the inverse Laplace transform in order to obtain semi‐analytical solutions in time domain. It is shown that the present solutions are more general and have a good agreement with the existing solutions from literatures. Furthermore, the current solutions can also be degenerated into conventional solutions to one‐dimensional consolidation of unsaturated soils with homogeneous boundaries. Finally, several numerical examples are provided to illustrate consolidation behavior of unsaturated soils under four types of time‐dependent loadings, including instantaneous loading, ramp loading, exponential loading and sinusoidal loading. Parametric studies are illustrated by variations of pore‐air pressure, pore‐water pressure and settlement at different values of the ratio of air–water permeability coefficient, depth and loading parameters. Copyright © 2017 John Wiley & Sons, Ltd.     

Numerical solution for consolidation and desiccation of soft soils

The consolidation and desiccation behaviour of soft soils can be described by two time‐dependent non‐linear partial differential equations using the finite strain theory. Analytical solutions do not exist for these governing equations. In this paper, we develop efficient numerical methods and software for finding the numerical solutions. We introduce a semi‐implicit time integration scheme, and show numerically that our method converges. In addition, the numerical solution matches well with the experimental result. A boundary refinement method is also developed to improve the convergence and stability for the case of Neumann type boundary conditions. Interface governing equations are derived to maintain the continuity of consolidation and desiccation processes. This is useful because the soil column can undergo desiccation on top and consolidation on the bottom simultaneously. The numerical algorithms has been implemented into a computer program and the results have been verified with centrifuge test results conducted in our laboratory. Copyright © 2001 John Wiley & Sons, Ltd.     

Response of multilayered transversely isotropic medium due to axisymmetric loads

A novel procedure associated with the precise integration method (PIM) and the technique of dual vector is proposed to effectively calculate the magnitude and distribution of deformations in a homogeneous multilayered transversely isotropic medium. The planes of transverse isotropy are assumed to be parallel to the horizontal surface of the soil system. The linearly elastic medium is subjected to four types of vertically acting axisymmetric loads prescribed either at the external surface or in the interior of the soil medium. There are no limits for the thicknesses and number of soil layers to be considered. By virtue of the governing equations of motion and the constitutive equations of the transversely isotropic elastic body, and based on the Hankel integral transform and a dual vector formulation in a cylindrical coordinate system, the partial differential motion equations can be converted into first‐order ordinary differential matrix equations. Applying the approach of PIM, it is convenient to obtain the solutions of ordinary differential matrix equations for the continuously homogeneous multilayered transversely isotropic elastic soil in the transformed domain. The PIM is a highly accurate algorithm to solve the sets of first‐order ordinary differential equations, which can ensure to achieve any desired accuracy of the solutions. What is more, all calculations are based on the standard method with the corresponding algebraic operations. Computational efforts can be reduced to a great extent. Finally, numerical examples are provided to illustrate the accuracy and effectiveness of the proposed approach. Some more cases are analyzed to evaluate the influences of the elastic parameters of the transversely isotropic media on the load‐displacement responses. Copyright © 2015 John Wiley & Sons, Ltd.     

Semi-analytical solutions to one-dimensional consolidation for unsaturated soils with symmetric semi-permeable drainage boundary

This paper presents semi-analytical solutions to Fredlund and Hasan’s one-dimensional consolidation for unsaturated soils under symmetric semi-permeable drainage boundary conditions. Two variables are introduced to transform two coupled governing equations of pore-air and pore-water pressures into an equivalent set of partial differential equations, which are easily solved by the Laplace transform. Then, the pore-air and pore-water pressures, and soil settlement are obtained in the Laplace domain. Crump’s method is adopted to perform the inverse Laplace transform in order to obtain semi-analytical solutions in time domain. It is shown that the present solution is more applicable to various types of drainage boundary conditions, and in a good agreement with existing solutions from the literature. Furthermore, several numerical examples are provided to investigate the consolidation behavior of an unsaturated single-layer soil with traditional drainage boundary (single or double), and single-sided and double-sided semi-permeable drainage boundaries. Finally, it illustrates the changes in pore-air and pore-water pressures and soil settlement with time at different values of symmetric semi-permeable drainage boundary conditions parameters. In addition, parametric studies are conducted by the variations of pore-air and pore-water pressures at different ratios of air-water permeability coefficient and the depth.     

Axisymmetric thermal consolidation of multilayered porous thermoelastic media due to a heat source

This paper presents an analytical layer element solution to axisymmetric thermal consolidation of multilayered porous thermoelastic media containing a deep buried heat source. By applying the Laplace–Hankel transform to the state variables involved in the basic governing equations of porous thermoelasticity, the analytical layer elements that describe the relationship between the transformed generalized stresses and displacements of a finite layer and a half‐space are derived. The global stiffness matrix equation is obtained by assembling the interrelated layer elements, and the real solutions in the physical domain are achieved by numerical inversion of the Laplace–Hankel transform after obtaining the solutions in the transformed domain. Finally, numerical calculations are performed to demonstrate the accuracy of this method and to investigate the influence of heat source's types, layering, and the porous thermoelastic material parameters on thermal consolidation behavior. Copyright © 2015 John Wiley & Sons, Ltd.     

Reflection–transmission matrix method for the consolidation of a multilayered saturated soil

The consolidation of the layered saturated soil is an important issue in civil engineering and has been investigated extensively during the past decades. In this study, based on the Biot's theory, the reflection–transmission matrix (RTM) method for treating the layered saturated soil under axisymmetric consolidation is developed. To decouple the governing equations of the Biot's theory, the McNamee displacement functions are introduced, and the general solution for the saturated soil is obtained using the Laplace and Hankel transforms. In order to develop the RTM method for the layered saturated soil, based on the obtained general solution, the static wave vector corresponding to the state vector of the saturated soil and the transform matrix relating the aforementioned two vectors are defined. Also, the transfer matrices corresponding to the two vectors are introduced, and the representations of the RTMs for the static wave vector of the saturated soil are presented. As the state vector, static wave vector, and the transform matrix relating the two vectors are all defined in the global coordinate system, the RTMs obtained in this study thus have a reasonable physical meaning. By using the RTMs for the layered saturated soil, the solutions for the layered saturated soil subjected to external sources are derived. Comparison of results due to the proposed RTM method with some existing results and results due to the transfer matrix method validates the developed RTM method. Some numerical results are obtained based on the proposed RTM method for the layered saturated soil. Copyright © 2016 John Wiley & Sons, Ltd.     

A simple analytical solution to one‐dimensional consolidation for unsaturated soils

This paper presents a simple analytical solution to Fredlund and Hasan's one‐dimensional (1‐D) consolidation theory for unsaturated soils. The coefficients of permeability and volume change for unsaturated soils are assumed to remain constant throughout the consolidation process. The mathematical expression of the present solution is much simpler compared with the previous available solutions in the literature. Two new variables are introduced to transform the two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved with standard mathematical formulas. It is shown that the present analytical solution can be degenerated into that of Terzaghi consolidation for fully saturated condition. The analytical solutions to 1‐D consolidation of an unsaturated soil subjected to instantaneous loading, ramp loading, and exponential loading, for different drainage conditions and initial pore pressure conditions, are summarized in tables for ease of use by practical engineers. In the case studies, the analytical results show good agreement with the available analytical solution in the literature. The consolidation behaviors of unsaturated soils are investigated. The average degree of consolidation at different loading patterns and drainage conditions is presented. The pore‐water pressure isochrones for two different drainage conditions and three initial pore pressure distributions are presented and discussed. Copyright © 2013 John Wiley & Sons, Ltd.     

半透水边界下分数阶黏弹性饱和土一维固结特性分析

为描述饱和土体的流变特性,引入分数阶导数Kelvin-Voigt黏弹性模型,采用解析方法对半透水边界下的分数阶黏弹性饱和土一维固结特性进行了研究。分别对骤加恒载下饱和土一维固结微分方程和分数阶Kelvin-Voigt黏弹性本构方程进行Laplace变换,并联立求解得到了双边半透水边界条件下分数阶黏弹性饱和土在Laplace变换域内的解析表达式。通过Crump方法实现Laplace数值反演,得到时间域内的半解析解。将所得到的解分别退化为分数阶黏弹性饱和土一维固结半解析解和双边半透水黏弹性饱和土一维固结半解析解,结果与已有文献半解析解相同,验证了提出的双边半透水边界条件下分数阶黏弹性饱和土一维固结解的可靠性。通过算例考察了半透水边界条件和分数阶黏弹性饱和土参数对一维固结特性的影响。研究表明,双边半透水边界下分数阶黏弹性饱和土一维固结发展过程与半透水边界条件、分数阶次和黏滞系数有关,且土体的压缩模量对饱和土一维固结最终沉降量有显著影响。