Semi-analytical solution to one-dimensional consolidation in unsaturated soils

In this paper, a series of semi-analytical solutions to one-dimensional consolidation in unsaturated soils are obtained. The air governing equation by Fredlund for unsaturated soils consolidation is simplified. By applying the Laplace transform and the Cayley-Hamilton theorem to the simplified governing equations of water and air, Darcy＇s law, and Fick＇s law, the transfer function between the state vectors at top and at any depth is then constructed. Finally, by the boundary conditions, the excess pore-water pressure, the excess pore-air pressure, and the soil settlement are obtained under several kinds of boundary conditions with the large-area uniform instantaneous loading. By the Crump method, the inverse Laplace transform is performed, and the semi-analytical solutions to the excess pore-water pressure, the excess pore-air pressure, and the soils settlement are obtained in the time domain. In the case of one surface which is permeable to air and water, comparisons between the semi-analytical solutions and the analytical solutions indicate that the semi-analytical solutions are correct. In the case of one surface which is permeable to air but impermeable to water, comparisons between the semi-analytical solutions and the results of the finite difference method are made, indicating that the semi-analytical solution is also correct.     

CONSOLIDATION THEORY OF UNSATURATED SOIL BASED ON THE THEORY OF MIXTURE(Ⅱ)

The present paper uses the mathematics model for consolidation of unsaturatedsoil developed in ref.[1]to solve boundary value problems.The analytical solutionsfor one-dimensional consolidation problem are gained by making use of Laplacetransform and finite Fourier transform.The displacement and the pore water pressureas well as the pore gas pressure are found from governing equations simultaneously.The theoretical formulae of coefficient and degree of consolidation are also given inthe paper.With the help of the method of Galerkin Weighted Residuals,the finiteelement equations for two-dimensional consolidation problem are derived.A FORTRANprogram named CSU8 using8-node isoparameter element is designed.A plane strainconsolidation problem is solved using the program,and some distinguishing features onconsolidation of unsaturated soil and certain peculiarities on numerical analysis arerevealed.These achievements make it convenient to apply the theory proposed by theauthor in engineering practice.     

A new analytical solution to axisymmetric Blot's consolidation of a finite soil layer

A new analytical method is presented to study the axisymmetric Blot's consolidation of a finite soil layer. Starting from the governing equations of axisymmetric Blot's consolidation, and based on the property of Laplace transform, the relation of basic variables for a point of a finite soil layer is established between the ground surface (z= 0) and the depth z in the Laplace and Hankel transform domains. Combined with the boundary conditions of the finite soil layer, the analytical solution of any point in the transform domain can be obtained. The actual solution in the physical domain can be obtained by inverse Laplace and Hankel transforms. A numerical analysis for the axisymmetric consolidation of a finite soil layer is carried out.     

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

The semi-analytical solutions to Fredlund and Hasan's one-dimensional(1 D)consolidation for unsaturated soils with a semi-permeable drainage boundary are presented. Two 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(PDFs), which are easily solved by the Laplace transform method. Then, the pore-water pressure, pore-air pressure, and soil settlement are obtained in the Laplace domain. The Crump method is adopted to perform the inverse Laplace transform in order to obtain the semi-analytical solutions in the time domain. It is shown that the proposed solutions are more applicable to various types of boundary conditions and agree well with the existing solutions from the literature. Several numerical examples are provided to investigate the consolidation behavior of an unsaturated single-layer soil with single, double, mixed, and semi-permeable drainage boundaries. The changes in the pore-air and pore-water pressures and the soil settlement with the time factor at different values of the semi-permeable drainage boundary parameters are illustrated. In addition, parametric studies are conducted on the pore-air and pore-water pressures at different ratios(the air permeability coefficient to the water permeability coefficient) and depths.     

A new analytical solution to axisymmetric Biot＇s consolidation of a finite soil layer

A new analytical method is presented to study the axisymmetric Biot＇s consolidation of a finite soil layer. Starting from the governing equations of axisymmetric Blot＇s consolidation, and based on the property of Laplace transform, the relation of basic variables for a point of a finite soil layer is established between the ground surface （z= 0） and the depth z in the Laplace and Hankel transform domains. Combined with the boundary conditions of the finite soil layer, the analytical solution of any point in the transform domain can be obtained. The actual solution in the physical domain can be obtained by inverse Laplace and Hankel transforms. A numerical analysis for the axisymmetric consolidation of a finite soil layer is carried out.     

A new method for solving Biot＇s consolidation of a finite soil layer in the cylindrical coordinate system

A new method is developed to solve Biot＇s consolidation of a finite soil layer in the cylindrical coordinate system. Based on the governing equations of Biot＇s consolidation and the technique of Laplace transform, Fourier expansions and Hankel transform with respect to time t, coordinate θ and coordinate r, respectively, a relationship of displacements, stresses, excess pore water pressure and flux is established between the ground surface （z = 0） and an arbitrary depth z in the Laplace and Hankel transform domain. By referring to proper boundary conditions of the finite soil layer, the solutions for displacements, stresses, excess pore water pressure and flux of any point in the transform domain can be obtained. The actual solutions in the physical domain can be acquired by inverting the Laplace and the Hankel transforms.     

Plane strain consolidation of soil layer with anisotropic permeability

This paper presents an alternative analytical technique to study a plane strain consolidation of a poroelastic soil by taking into account the anisotropy of permeability. From the governing equations of a saturated poroelastic soil, the relationship of basic variables for a point of a soil layer is established between the ground surface （z=0） and the depth z in the Laplace-Fourier transform domain. Combined with the boundary conditions, an exact solution is derived for plane strain Biot＇s consolidation of a finite soil layer with anisotropic permeability in the transform domain. Numerical inversions of the Laplace transform and the Fourier transform are adopted to obtain the actual solution in the physical domain. Numerical results of plane strain Biot＇s consolidation for a single soil layer show that the anisotropic of permeability has a great influence on the consolidation behavior of the soils.     

Dynamic analysis of functionally graded plates using the hybrid Fourier-Laplace transform under thermomechanical loading

In this study, the analytical solution is presented for dynamic response of a simply supported functionally graded rectangular plate subjected to a lateral thermomechanical loading. The first-order and third-order shear deformation theories and the hybrid Fourier-Laplace transform method are used. The material properties of the plate, except Poisson’s ratio, are assumed to be graded in the thickness direction according to a power-law distribution in terms of the volume fractions of the constituents. The plate is subjected to a heat flux on the bottom surface and convection on the upper surface. A third-order polynomial temperature profile is considered across the plate thickness with unknown constants. The constants are obtained by substituting the profile into the energy equation and applying the Galerkin method. The obtained temperature profile is considered along with the equations of motion. The governing partial differential equations are solved using the finite Fourier transformation method. Using the Laplace transform, the unknown variables are obtained in the Laplace domain. Applying the analytical Laplace inverse method, the solution in the time domain is derived. The computed results for static, free vibration, and dynamic problems are presented for different power law indices for a plate with simply supported boundary conditions. The results are validated with the known data reported in the literature. Furthermore, the results calculated by the analytical Laplace inversion method are compared with those obtained by the numerical Newmark method.     

条形荷载下梯度非均匀非饱和土的动力响应分析

基于多孔介质混合物理论, 建立了梯度非均匀非饱和土地基模型, 研究了条形荷载作用下梯度非均匀非饱和土地基的动力响应问题. 通过傅里叶积分变换和Helmholtz矢量分解原理, 获得频域内非饱和土地基动力响应问题的通解, 结合回传射线矩阵法和边界条件, 求解获得了非均匀非饱和土层中位移、应力以及孔隙压力的计算列式. 假设沿深度方向梯度非均匀非饱和土的物理力学性质按幂函数连续变化, 通过数值傅里叶逆变换得到了非均匀非饱和土地基中的应力、位移以及孔隙压力等物理量的数值解, 分析讨论了土体非均匀性对非饱和土介质动力响应的影响规律. 结果表明: 土体非均匀性显著改变了非饱和土中竖向位移、正应力和孔隙压力在其深度方向上的振动模态, 其中孔隙气压在其深度方向的振动频率随着梯度因子的增加而不断增大, 波峰值不断靠近地表处附近; 竖向位移随着梯度因子的增大不断减小; 正应力和孔隙水压随着梯度因子的增大先增大后减小, 并且土体非均匀程度越高, 正应力与孔隙水压的幅值越大.      

Two-dimensional transient thermo-hydro-mechanical fundamental solutions of multiphase porous media in frequency and time domains

In this paper, the closed form two-dimensional fundamental solutions for a non-isothermal unsaturated deformable porous medium have been derived for a symmetric polar domain in both Laplace transform and time domains. The governing differential equations of the non-isothermal unsaturated soil consist of equilibrium, moisture, air and heat transfer equations including the suction effect, temperature effect and dissolved air in water. The derived fundamental solution has been verified mathematically by comparison with the previously presented corresponding fundamental solutions in three limiting cases including the steady-state thermo-hydro-mechanical, steady-state hydro-mechanical and elastostatic fundamental solutions. Also these 2D kernel functions are tested in comparison with a finite element method (FEM).     

Exact solutions of multi-term fractional difusion-wave equations with Robin type boundary conditions

General exact solutions in terms of wavelet expansion are obtained for multiterm time-fractional difusion-wave equations with Robin type boundary conditions. By proposing a new method of integral transform for solving boundary value problems, such fractional partial diferential equations are converted into time-fractional ordinary diferential equations, which are further reduced to algebraic equations by using the Laplace transform. Then, with a wavelet-based exact formula of Laplace inversion, the resulting exact solutions in the Laplace transform domain are reversed to the time-space domain. Three examples of wave-difusion problems are given to validate the proposed analytical method.     

General semi-analytical solutions to one-dimensional consolidation for unsaturated soils

This paper presents general semi-analytical solutions to Fredlund and Hasan's one-dimensional(1D) consolidation equations for unsaturated soils subject to different initial conditions, homogeneous boundaries and time-dependent loadings. Two variables are introduced to transform the two-coupled governing equations of pore-water and poreair pressures into an equivalent set of partial differential equations(PDEs), which are solved with the Laplace transform method. The pore-water and pore-air pressures and settlement are obtained in the Laplace transform domain. The Crump's method is used to perform inverse Laplace transform to obtain the solutions in the time domain. The present solutions are more general in practical applications and show good agreement with the previous solutions in the literature.     

Exact solutions of multi-term fractional diffusion-wave equations with Robin type boundary conditions

General exact solutions in terms of wavelet expansion are obtained for multi- term time-fractional diffusion-wave equations with Robin type boundary conditions. By proposing a new method of integral transform for solving boundary value problems, such fractional partial differential equations are converted into time-fractional ordinary differ- ential equations, which are further reduced to algebraic equations by using the Laplace transform. Then, with a wavelet-based exact formula of Laplace inversion, the resulting exact solutions in the Laplace transform domain are reversed to the time-space domain. Three examples of wave-diffusion problems are given to validate the proposed analytical method.     

Analysis of One-Dimensional Advection–Diffusion Model with Variable Coefficients Describing Solute Transport in a Porous medium

This work presents the analytical solution and temporal moments of one-dimensional advection–diffusion model with variable coefficients. Two case studies along with the two different sets of boundary conditions are considered at the inlet and outlet of the domain. In the first case, a time-dependent solute dispersion in the homogeneous domain along uniform flow is taken into account, whereas in the second case, due to inhomogeneity of domain, velocity is taken spatially dependent and the dispersion is assumed proportional to the square of the velocity. The Laplace transform is used to obtain the analytical solutions. The analytical temporal moments are derived from the Laplace domain solutions. To verify the correctness of the analytical solutions, a high-resolution second-order finite volume scheme is applied. Different case studies are considered and discussed. Both analytical and numerical results are in good agreement with each other.     

A Quasilinear Based Procedure for Saturated/Unsaturated Water Movement in Soils

The quasilinear form of Richards equation for one-dimensional unsaturated flow in soils can be readily solved for a wide variety of conditions. However, it cannot explain saturated/unsaturated flow and the constant diffusivity assumption, used to linearise the transient quasilinear equation, can introduce significant error. This paper presents a quasi-analytical solution to transient saturated/unsaturated flow based on the quasilinear equation, with saturated flow explained by a transformed Darcy's equation. The procedure presented is based on the modified finite analytic method. With this approach, the problem domain is divided into elements, with the element equations being solutions to a constant coefficient form of the governing partial differential equation. While the element equations are based on a constant diffusivity assumption, transient diffusivity behaviour is incorporated by time stepping. Profile heterogeneity can be incorporated into the procedure by allowing flow properties to vary from element to element. Two procedures are presented for the temporal solution; a Laplace transform procedure and a finite difference scheme. An advantage of the Laplace transform procedure is the ability to incorporate transient boundary condition behaviour directly into the analytical solutions. The scheme is shown to work well for two different flow problems, for three soil types. The technique presented can yield results of high accuracy if the spatial discretisation is sufficient, or alternatively can produce approximate solutions with low computational overheads by using large sized elements. Error was shown to be stable, linearly related to element size.     

Effect of quadratic pressure gradient term on a one-dimensional moving boundary problem based on modified Darcy’s law

A relatively high formation pressure gradient can exist in seepage flow in low-permeable porous media with a threshold pressure gradient, and a significant error may then be caused in the model computation by neglecting the quadratic pressure gradient term in the governing equations. Based on these concerns, in consideration of the quadratic pressure gradient term, a basic moving boundary model is constructed for a one-dimensional seepage flow problem with a threshold pressure gradient. Owing to a strong nonlinearity and the existing moving boundary in the mathematical model, a corresponding numerical solution method is presented. First, a spatial coordinate transformation method is adopted in order to transform the system of partial differential equations with moving boundary conditions into a closed system with fixed boundary conditions; then the solution can be stably numerically obtained by a fully implicit finite-difference method. The validity of the numerical method is verified by a published exact analytical solution. Furthermore, to compare with Darcy’s flow problem, the exact analytical solution for the case of Darcy’s flow considering the quadratic pressure gradient term is also derived by an inverse Laplace transform. A comparison of these model solutions leads to the conclusion that such moving boundary problems must incorporate the quadratic pressure gradient term in their governing equations; the sensitive effects of the quadratic pressure gradient term tend to diminish, with the dimensionless threshold pressure gradient increasing for the one-dimensional problem.     

Transient Heat and Moisture Flow Around Heat Source Buried in an Unsaturated Half Space

This article presents solutions for the transient heat and moisture transport due to both disk heat source and cylindrical heat source buried in an unsaturated half space. The solutions are presented in Hankel–Laplace transform domain and in dimensionless style. Coupled effect of thermally driven moisture transport is especially investigated because of its importance to alter the flow field in low-permeability medium. Parametric study has been performed to assess the effects of five independent dimensionless parameters on flow field. The stability and accuracy of the present solutions are demonstrated from the comparison between the results obtained from these solutions and those by using a well-established finite element code CODE_BRIGHT. Despite the simplified assumptions required in order to obtain analytical solutions in Hankel–Laplace transform domain, the results incorporate the main mechanisms involved in the coupled thermo-hydraulic (T-H) problem, and they may be eventually used for validation purposes.     

Linear viscoelastic analysis of a semi-infinite porous medium

This paper considers the problem of a semi-infinite, isotropic, linear viscoelastic half-plane containing multiple, non-overlapping circular holes. The sizes and the locations of the holes are arbitrary. Constant or time dependent far-field stress acts parallel to the boundary of the half-plane and the boundaries of the holes are subjected to uniform pressure. Three types of loading conditions are assumed at the boundary of the half-plane: a point force, a force uniformly distributed over a segment, a force uniformly distributed over the whole boundary of the half-plane. The solution of the problem is based on the use of the correspondence principle. The direct boundary integral method is applied to obtain the governing equation in the Laplace domain. The unknown transformed displacements on the boundaries of the holes are approximated by a truncated complex Fourier series. A system of linear equations is obtained by using a Taylor series expansion. The viscoelastic stresses and displacements at any point of the half-plane are found by using the viscoelastic analogs of Kolosov–Muskhelishvili’s potentials. The solution in the time domain is obtained by the application of the inverse Laplace transform. All the operations of space integration, the Laplace transform and its inversion are performed analytically. The method described in the paper allows one to adopt a variety of viscoelastic models. For the sake of illustration only one model in which the material responds as the standard solid in shear and elastically in bulk is considered. The accuracy and efficiency of the method are demonstrated by the comparison of selected results with the solutions obtained by using finite element software ANSYS.     

The Flow Problem of Fluids Flow in a Fractal Reservoir with Double Porosity

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｆｌｏｗｔｈｅｏｒｙａｎｄｉｔｓａｐｐｌｉｃａｔｉｏｎｏｆｆｌｕｉｄｓｆｌｏｗｉｎａｆｒａｃｔａｌｒｅｓｅｒｖｏｉｒｈａｖｅｃｏｎｔｉｎｕａｌｌｙｇｏｎｅｄｅｅｐｉｎｔｏｓｔｕｄｙｓｉｎｃｅＣｈａｎｇａｎｄＹｏｒｔｓｏｓ[1]ｂｕｉｌｔｔｈｅｆｌｏｗｍｏｄｅｌｏｆｆｌｕｉｄｔｈｒｏｕｇｈａｆｒａｃｔａｌｒｅｓｅｒｖｏｉｒ.ＴＯＮＧＤｅｎｇ＿ｋｅ[2 ]ｐｒｅｓｅｎｔｅｄｔｈｅｅｘａｃｔｓｏｌｕｔｉｏｎａｎｄｉｔｓｐｒｅｓｓｕｒｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｆｏｒｔｈｅｖａ…     

Plane strain consolidation of a multi-layered soil with anisotrpic permeability and compressible constituents

对多层地基的平面应变固结问题进行了研究,并同时考虑了土体的渗透各向异性和孔隙 流体的可压缩性. 从平面应变Biot固结的控制方程出发,对时间t, 坐标z和x进行 Laplace和Fourier变换,建立了地基表面(z=0)和任意深度z处的基本量 在Laplace-Fourier变换域内的传递矩阵关系. 利用传递矩阵 法,结合土层连续条件和边界条件,并应用Laplace-Fourier逆变换技术,推导出渗透各向 异性可压缩多层地基平面应变固结的理论解. 基于该解,编制了计算程序,并进行了 数值计算. 讨论了土体的渗透各向异性、孔隙流体的可压缩性以及地基的分层特性对地基固 结的影响,分析结果表明：土体的渗透各向异性、孔隙流体的可压缩性,以及地基的分层特 性对地基的固结行为有着重要的影响.