变荷载下层状非饱和土地基全耦合固结特性研究

非饱和土是地球表层土体一种常见的存在形式。与经典的饱和土Biot固结理论相比,非饱和土固结理论还亟待发展。基于Fredlund非饱和土双应力变量固结理论,放弃传统理论固结过程中骨架总应力不变的假设,推导出变荷载作用下非饱和土全耦合轴对称固结理论的控制方程组。通过Laplace-Hankel变换处理变量t和r,所得的控制方程被处理为常微分方程组;扩展的精细积分法将被进一步运用求解方程组以得到层状非饱和土地基在变换域内的固结解答,而最终的解答将通过Laplace-Hankel逆变换技术求得。将非饱和土地基退化为饱和土地基,与现有文献结果进行对比,验证所提方法结果的可靠性;最后,提供3个数值算例,以讨论加载时间T0、孔隙水关于净应力的体积变化系数m1w以及土体分层性对非饱和土固结特性的影响。     

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.     

非饱和土热湿耦合传输问题的数值解——理论及一维问题的解

本文介绍了求解非饱和土壤中热量和水分耦合传输问题的一种数值方法——积分有限差分方法(IFDM)。基于菲利普-迪弗瑞斯(Philip-De Vries)多孔介质热湿耦合流动模型,采用积分有限差分方法,编制了求解热湿传输问题的计算程序 HM1,该程序可用于求解各类工程中遇到的多孔介质一维传热传湿问题。文中还给出了计算结果与实测结果相比较的实例。     

Numerical analysis of seepage–deformation in unsaturated soils

A coupled elastic–plastic finite element analysis based on simplified consolidation theory for unsaturated soils is used to investigate the coupling processes of water infiltration and deformation. By introducing a reduced suction and an elastic–plastic constitutive equation for the soil skeleton, the simplified consolidation theory for unsaturated soils is incorporated into an in-house finite element code. Using the proposed numerical method, the generation of pore water pressure and development of deformation can be simulated under evaporation or rainfall infiltration conditions. Through a parametric study and comparison with the test results, the proposed method is found to describe well the characteristics during water evaporation/infiltration into unsaturated soils. Finally, an unsaturated soil slope with water infiltration is analyzed in detail to investigate the development of the displacement and generation of pore water pressure.     

非饱和土二维固结简化计算的研究

针对较高饱和度的非饱和土二维固结问题展开研究。对于较高饱和度的非饱和土,如饱和度 时,将孔隙中气、水近似地看作可压缩的气、水混合流体后,非饱和土近似为土骨架和混合流体的二相土。考虑混合流体的压缩性,建立混合流体的连续方程。联立平衡方程和混合流体的连续方程,求出应力-应变和混合流体压力;再建立水连续方程求解水压力,继而求出气压力、吸力等。算例表明：加载和消散过程中,混合流体与水压力变化基本一致,气压的作用并不大;地基变形过程与高速公路填筑过程中地基变形发生规律一致。说明该简化方法是合理的,并促进了非饱和土固结变形计算走向实用化。     

Analytical solution for the 1D consolidation of unsaturated multi-layered soil

The 1D consolidation of unsaturated multi-layered soil is studied based on the theory proposed by Fredlund and Hasan, and an analytical solution for a typical boundary condition is obtained by assuming all material parameters remain constant during consolidation. In the derivation of the analytical solution, the eigenfunction and eigenvalue for the multi-layered problem are first derived through the transfer matrix method. Then, by using the method of undetermined coefficients and the orthogonal relation of the eigenfunction, the analytical solution is obtained. The present method is applicable to various types of boundary conditions. Finally, numerical examples are provided to investigate the consolidation behavior of unsaturated multi-layered soil.     

Gibson地基的一维固结解

针对单层Gibson地基模型,运用Laplace变换,求得了单层Gibson地基的一维固结问题,得到了频域内的通解。通过Laplace逆变换,即可计算单层Gibson地基在任意荷载下的一维固结特性。此外,结合工程实例,对解进行了探讨,揭示了单层Gibson地基的固结特性。     

非饱和红层路堤沉降理论与应用研究

路堤本体主要是非饱和土,其变形由3部分组成,即瞬时压密产生的变形、主固结产生的变形和次固结产生的变形,因此其理论也应反映这3部分变形。以单轴压缩红层蠕变试验为基础,提出基于非饱和土单变量理论的红层非饱和土路堤沉降计算理论;通过算例对遂渝客运专线遂宁段某工点实际填筑过程的进行了数值模拟,并对其工后沉降进行了预则,结果表明,模拟结果与现场监测数据预测的结果相符,说明该理论是正确的?     

非饱和土中热−湿−盐耦合作用的稳态分析

以多孔介质理论为基础,研究了稳态条件下非饱和土中温度?水分?盐分多场耦合问题。考虑非饱和土的孔隙被液态水、溶解的盐分、水蒸气和干燥气体等填充,在质量和能量守恒的基础上获得了非饱和土中水分、气体、盐分的质量守恒方程以及能量守恒方程。考虑一维稳态问题,选取温度、孔隙气压、孔隙水压和盐溶液浓度以及它们的导数作为状态变量,得到了问题的状态方程组。在给定的边界条件下,采用打靶法求解了该强耦合的非线性变系数微分方程组,通过与已有的试验结果相比较,验证了模型的有效性。基于数值算例,参数分析了含水率、温度边界、孔隙率等条件对非饱和土中温度场、水分场和盐分场分布的影响规律。     

膨胀土路基变形的数值模拟

在非饱和土固结理论的基础上,分析了降雨入渗条件下全填筑和全开挖膨胀土路基变形的特征,获取了浸水后膨胀土路基内塑性区的发展过程,得出了膨胀土路基一般发生浅层滑动的结论,这种滑动与通常饱和土的边坡失稳的原因是不同的。     

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

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

分段循环荷载作用下二维非饱和土固结特性分析

基于Dakshanamurthy和Fredlund提出的二维非饱和土固结理论,利用Fourier正弦级数展开、Laplace变换,分别给出了分段循环荷载作用下二维非饱和土固结问题的超孔隙气压力、超孔隙水压力和沉降的半解析解,并应用退化法验证了本文所得半解析解的正确性。然后,结合3种具体的荷载形式,分析了分段循环荷载作用下气相与液相渗透系数之比(k_a/k_w)、水平方向与竖直方向渗透系数之比(k_x/k_z)和荷载特征参数(a)对二维非饱和土固结特性的影响。结果表明：k_a/k_w和k_x/k_z的增大均会加速固结沉降进程;荷载特征参数越大,沉降发展越早,沉降值越小;二维非饱和土固结特性受分段循环荷载作用影响明显。因此,在实际施工过程中改变施工速度和设置径向排水装置可有效控制二维条件下非饱和土体的固结过程,该研究成果可为非饱和土地基的设计和施工提供重要理论依据。     

Semi-analytical solution to one-dimensional consolidation for viscoelastic unsaturated soils

This paper presents a semi-analytical solution to one-dimensional consolidation of viscoelastic unsaturated soils with a finite thickness under oedometric conditions and subjected to a sudden loading. The solution is obtained by using Lee’s correspondence principle based on the semi-analytical solution to one-dimensional consolidation of elastic unsaturated soils. The boundary contains the top surface permeable to water and air and the bottom impermeable to water and air. A typical example is given to show the evolution of excess pore-air and pore-water pressures as well as the total degree of consolidation of the soil layer with time for different ratios of air–water permeability coefficient, elastic modulus and viscoelastic coefficient. The one-dimensional consolidation behavior of viscoelastic unsaturated soil is discussed according to the semi-analytical solution. These results contribute to a better understanding of the consolidation behavior of viscoelastic unsaturated soils.     

Influence of gas generation in electro‐osmosis consolidation

This paper presents a numerical model for the elasto‐plastic electro‐osmosis consolidation of unsaturated clays experiencing large strains, by considering electro‐osmosis and hydro‐mechanical flows in a deformable multiphase porous medium. The coupled governing equations involving the pore water flow, pore gas flow, electric flow and mechanical deformation in unsaturated clays are derived within the framework of averaging theory and solved numerically using finite elements. The displacements of the solid phase, the pressure of the water phase, the pressure of the gas phase and the electric potential are taken as the primary unknowns in the proposed model. The nonlinear variation of transport parameters during electro‐osmosis consolidation are incorporated into the model using empirical expressions that strongly depend on the degree of water saturation, whereas the Barcelona Basic Model is employed to simulate the elasto‐plastic mechanical behaviour of unsaturated clays. The accuracy of the proposed model is evaluated by validating it against two well‐known numerical examples, involving electro‐osmosis and unsaturated soil behaviour respectively. Two further examples are then investigated to study the capability of the computational algorithm in modelling multiphase flow in electro‐osmosis consolidation. Finally, the effects of gas generation at the anode, the deformation characteristics, the degree of saturation and the time dependent evolution of the excess pore pressure are discussed. Copyright © 2016 John Wiley & Sons, Ltd.     

考虑渗透系数变化的非饱和土固结性状分析

基于Fredlund非饱和土一维固结理论,研究了有限厚度的表面透水透气、底面不透水不透气的线弹性和黏弹性非饱和土地基在加荷随时间指数性变化时的一维固结特性。分别得到了两类地基在固结过程中同时考虑液相、气相渗透系数非线性变化和仅考虑液相渗透系数变化两种情况下的半解析解答。利用典型算例进行计算,分析了不同情况下两类地基中超孔隙水、气压力消散以及地基固结度随时间的变化规律,并与不考虑渗透系数变化时的半解析解计算结果进行了对比。结果发现：固结过程中渗透系数呈非线性变化;只考虑液相渗透系数变化时,超孔隙气压力的消散变化不大,超孔隙水压力的消散加快;气相渗透系数变化对超孔隙气的消散产生明显影响,对超孔隙水压力消散影响不大。同时考虑液相和气相渗透系数变化时,土体中超孔隙水、气压力的消散均有明显变化,土体固结速度也相应加快;分析结果对非饱和土固结的进一步研究具有重要意义。     

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.     

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.     

粘土干湿循环中裂缝演变过程的数值模拟

在非饱和土简化固结理论基础上,对粘土在干湿过程中,表面裂缝从形成、发展到闭合的发展过程进行了数值模拟,模拟结果说明,本文采用的基于非饱和简化固结理论的计算方法是可行、实用的.     

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.     

不同饱和度土层分界面上剪切波的反射与透射

基于多相孔隙介质弹性理论,给出了非饱和土中不同弹性波的传播方程。根据分界面上的边界条件,建立了各势函数波幅值之间的关系式,讨论了入射剪切波在不同饱和度土层分界面上的反射与透射问题。在无限空间非饱和土体中存在3种压缩波和1种剪切波,因此,当剪切波传播到不同饱和度的非饱和土层分界面上将分别在上、下土层激发产生4种反射波和4种透射波。推导出不同反射波和透射波的振幅比例系数和能量比例系数的理论表达式,并且在此基础上进行数值分析。在数值算例中分别研究了各反射波与透射波的能量比例系数（即能量反射率和能量透射率）受入射频率、入射角度以及上、下土层土体饱和度变化的影响情况。计算结果表明：各能量反射率和能量透射率不仅与入射角和入射频率有关,而且其受上、下土层饱和度变化的影响也同样不能忽视。