水位波动下非饱和心墙土坝体的渗流和稳定性

基于Richards方程,对坝体的饱和-非饱和渗流场进行了模拟,再根据饱和-非饱和渗流场和非饱和土抗剪强度公式,对坝体的稳定性进行了分析。结果表明,水位骤降过程中坝体的安全系数通常呈现先缓慢增加后迅速减小的变化过程,分析坝体失稳时塑性区和位移场发现,水位下降的初期,坝体左侧坡体的安全系数要低于坝体右侧坡体,但水位下降到一定程度,右侧坡体的安全系数迅速减小,并先于左侧坡体失稳;采用有限元强度折减法用于多坡面边坡稳定分析时,只能获得最小安全系数的包络线;心墙具有隔水防渗的作用,对水位变化渗流具有阻尼作用。     

非饱和地震滑坡堆积体降雨破坏水-力耦合行为试验

降雨滑坡的水-力耦合过程和机制研究是开展准确预测预报的基础和前提。选取典型滑坡模型,采用足尺人工降雨模型试验模拟了典型砾石土滑坡堆积体破坏过程,并采用Brooks-Corey(BC)和Van Genuchten(VG)模型分别建立了降雨入渗吸湿过程的土-水特征曲线,再通过一维非饱和无限边坡稳定性分析模型对滑坡开展应力状态和稳定性分析,揭示了斜坡破坏过程是优先流和基质流的双渗流场共同作用的结果,而VG模型适用于吸湿条件下宽级配砾石土优先流的土-水特征曲线的重构;稳定性计算和试验过程的对比分析表明,试验中斜坡的变形破坏过程中观测到的含水率、基质吸力和地表倾斜角度的变化过程与计算的吸力应力、稳定系数的演进过程有显著的对应关系。通过试验还揭示了非饱和滑坡堆积体渗流潜蚀现象,定量分析了细颗粒迁移程度。基于物理模型试验,揭示了降雨激发条件下滑坡堆积体形成优先流作用的土-水特征曲线和变形破坏的水力耦合过程,可为地震滑坡堆积体非饱和失稳的预测预警提供参考和基础。     

非饱和黄土土-水特征曲线与渗透系数Childs & Collis-Geroge模型预测

我国黄土厚度大,黄土地区地下水位深,降雨量少,黄土多处于非饱和状态,土-水特征曲线是非饱和黄土应力状态、强度及渗透性研究的基础。以陇东高原马兰黄土为试验对象,采用张力计法测定原状土样的土-水特征曲线,采用三种理论模型对试验数据拟合,其中Gardner模型简单、参数少,但Fredlund & Xing模型拟合效果最好。基于已测得的土-水特征曲线,采用Childs & Collis-Geroge预测非饱和渗透系数的模型,计算得到非饱和黄土渗透系数与基质吸力或含水率的关系,发现黄土从饱和到非饱和,其渗透系数急剧降低,黄土非饱和渗透系数与基质吸力或体积含水率的关系均可用指数函数表示。本文对典型黄土土-水特征曲线的研究及渗透性的预测为黄土工程问题,如降雨入渗的边坡稳定性评价、非饱和地基湿陷变形计算等提供理论基础。     

非饱和堤岸的渗流特征及其稳定性研究

通过试验确定了非饱和黏土层的土-水特征曲线和强度参数,通过理论分析建立了非饱和堤岸渗流-应力耦合模型。在此基础上,利用耦合计算程序分析了非饱和土堤岸在河水位变动时的非稳定渗流场特征,结合强度折减有限元法分析了河水位反复升降后非饱和堤岸稳定性的变化。结果表明：水位快速上升时,堤脚渗透流速下降后又逐步上升并趋向稳定,河水位快速上升使堤岸边坡的稳定性降低,随着渗流场中孔隙水压力的调整,堤岸边坡的安全系数又有所回升;河水位骤降时,黏性土层饱和度变化相对滞后,水位骤降加大了堤岸的渗透流速,边坡的稳定性迅速降低,水位下降约120 h后堤岸边坡进入较危险时段;非饱和土的基质吸力提高了堤岸边坡整体稳定性,河水位反复升降降低了堤岸边坡整体稳定性,且河水位越低,水位反复升降对堤岸稳定性产生的影响越大     

非饱和土土水特征曲线(SWCC)测试与预测

非饱和土土水特征曲线(SWCC)表示了土中含水量与吸力之间的关系。文章介绍了6种常用方法,各有其适用范围。体积压力板仪可量测最大基质吸力值为1500kPa的干燥曲线和浸湿曲线;超过1500kPa时,可用盐溶液法进行量测;Tem-ple仪可量测基质吸力达100kPa的干燥曲线;滤纸法可用于测量土体的基质吸力与总吸力;Dew-point电位计可用于量测土样总吸力变化,尤其适合渗透吸力的量测;TDR探头适合于量测小于300kPa的基质吸力。用GDS非饱和土三轴仪可以进行SWCC测试,测试范围主要取决于陶土板的进气值。用准确的数学模型对测得的含水量、吸力数据进行拟合,对于预测非饱和土力学性质、渗透系数、抗剪强度及分析边坡稳定性有重要意义。由于准确测试SWCC难度较大,并且测试影响因素较多,所以根据土体孔隙大小分布和颗粒大小分布情况预测SWCC,也是一种较好的方法。     

土水特征曲线的概率分析

土水特征曲线(SWCC)反映非饱和土基质吸力与含水率的关系,是非饱和土研究中重要的关系曲线。为进行与非饱和土有关工程的可靠度分析,需要对土水特征曲线进行概率分析。文章基于理论分析及数值模拟试验,分析比较了求解曲线拟合参数的直接回归法及数理统计法,探讨了试验点数量及位置对计算结果的影响。研究了求解响应量置信带的直接回归法,以及由自变量的置信区间求解响应量置信带的间接方法,指出了考虑变量间相关性的重要性,提出了由相关非正态变量求解响应量置信带的方法。通过土水特征试验,对合肥非饱和土的SWCC概率特征进行了分析,求得了SWCC拟合参数的均值、变异系数及相关系数,绘出了SWCC响应量的置信范围与基质吸力的关系。     

甘肃永靖黑方台地区黄土滑坡土的力学性质

与地下水位以下的饱和黄土力学性质相比,包气带黄土力学性质具有不同的表征参数和测试方法,其对黄土斜坡渗流场和稳定性的影响不容忽视.以甘肃黑方台地区黄土为例,分别采用张力计和轴平移技术量测了黄土的部分区段土-水特征曲线,依据Gardner、Van Genuchten和Fredlund and Xing 3种经验公式对轴平移测试数据进行拟合,获取了完整的土-水特征曲线.基于GDS多功能三轴仪,通过控制基质吸力的三轴剪切试验,获取了考虑基质吸力的黄土强度参数.结果显示:张力计和轴平移技术量测的部分区段土-水特征曲线具有相近的曲线形态,但张力计曲线位于轴平移技术曲线的上方;当基质吸力由100kPa减小至20kPa时,黄土内摩擦角减小约1°,粘聚力下降近22kPa;拟合所得的黑方台地区黄土的基质吸力摩擦角为15.6°.试验为进行黑方台黄土斜坡的非饱和渗流和稳定性分析提供了参数.     

基于VG模型确定土水特征曲线基本参数

土水特征曲线基本参数(进气值、残余吸力值和反弯点的斜率等)在非饱和土的强度理论、渗流理论以及体变特性等方面都是非常重要的参数,常常通过传统画图法来确定,其精确度难以保证。以洛川标准剖面4层原状黄土为研究对象,基于VG模型提出单峰和双峰SWCC基本参数的确定方法;采用滤纸法测得SWCC数据点,基于VG模型进行最优化拟合,获取拟合参数,利用单峰和双峰SWCC数据点验证该方法的可行性;进行自然状态下的水分蒸发试验,根据质量含水率与蒸发时间的关系确定蒸发残余饱和度Srzf,依据Sr1、Sr2(Sr1和Sr2分别由确定残余状态的两种方法得到)与Srzf的相对误差值,比较了两种确定残余状态的方法。结果表明:对于单峰和双峰SWCC曲线,该方法都可以得到合理的基本参数并有效地确定其残余状态。     

降雨条件下土坡饱和-非饱和渗流及稳定性分析

基于非饱和土力学理论,利用有限元法对降雨条件下边坡的饱和-非饱和渗流及稳定性进行了探讨。数值分析表明,雨水入渗使边坡非饱和区土体的基质吸力减小,是导致边坡稳定性降低的主要因素。进而考察了基质吸力的不同取值及裂隙的不同位置对边坡稳定性的影响,结果表明：在降雨初期基质吸力取值越大,得到的边坡安全系数越高,但经过长时间降雨即边坡土体基本达到饱和后,基质吸力的取值不同对边坡安全系数基本上没有影响;裂隙的存在为雨水入渗提供了良好的通道,因而导致边坡稳定性相对于无裂隙状态来说有所降低。降雨初期下游裂隙使边坡稳定性降低较快,但随着降雨历时的增加,即下游土体逐渐变为饱和状态后,边坡的稳定性变化则主要由上游裂隙控制。     

密度和干湿循环对黄土土-水特征曲线的影响

土-水特征曲线在非饱和土力学中扮演着重要的角色。影响非饱和土土-水特征曲线的因素很多,如初始密度、干湿循环、荷载等。利用GCTS土-水特征曲线仪,主要测定了不同初始干密度的重塑非饱和黄土试样的减湿土-水特征曲线,选取合理的土-水特征曲线模型进行数据拟合,通过分析模型拟合参数来研究初始干密度对非饱和黄土土-水特征曲线的影响;并对给定初始干密度的试样进行了2次干湿循环试验,测定相应的减-增湿土-水特征曲线,通过对比干湿循环曲线上体积含水率的差异,提出了“滞回度”的概念,初步研究了黄土土-水特征曲线在干湿循环下的滞回特性规律。这些成果将为研究应力作用下的土-水特征曲线模型以及考虑干湿循环的力-水耦合本构模型发挥重要作用     

Influence of hydraulic characteristics on stability of unsaturated slope under transient seepage conditions

This paper examines the current procedure for determining the soil-water characteristic curve (SWCC) model with a particular focus on its application to slope stability analysis under transient unsaturated seepage conditions. A series of laboratory experiments was performed to determine the SWCC of different soils, ranging from high plasticity clay to silty sand, found across the Korean Peninsula. The experimental results were utilized to identify the suitable SWCC model for each soil type based on the fitting criterion. Also, this paper developed a numerical framework for infinite slope stability analysis under transient unsaturated seepage conditions. The significant advantage of the proposed framework, from the practical viewpoint, is to directly predict the timing of failure and potential failure plane based on rainfall recording. The effect of choice of SWCC models on predictability in stability analysis was evaluated by adopting the present framework along with the identified SWCC models. Furthermore, a case study of landslides after a 3-month rainfall in Pohang, Korea, was revisited to assess the performance of the proposed framework. The obtained results demonstrate the significant role of SWCC model on the results of slope stability analysis. The analysis using the SWCC model satisfying the fitting criterion could still not capture the real behavior of unsaturated soil. The comprehensive transient analysis is strongly suggested as a complementary means to the current fitting criterion for determining the suitable SWCC model for stability analysis under transient seepage conditions.     

Water characteristic curve of soil with bimodal grain-size distribution

Soil–water characteristic curve (SWCC) is the most fundamental and important soil property in unsaturated soil mechanics. It has been used for analyzing slope stability due to the infiltration of rainfall into slopes and water flow in unsaturated embankments. Generally, SWCC is obtained by laboratory tests. However high cost, long duration and difficulty of the tests impede the application of unsaturated soil mechanics to practical design or analysis. Therefore, several equations have been developed to predict the SWCC using grain-size distribution (GSD) curve. However, most of the equations were limited to soils with unimodal characteristics and the parameters of the equations are not related to the physical properties of the soil. In this paper, an equation to predict SWCC for soils with bimodal characteristics is proposed. The parameters of the proposed equation are related to the physical properties of soil and the variables of SWCC closely. The proposed equation is evaluated with data from the literature and laboratory tests carried out in this study. In addition, the computer codes for the computation of the predicted bimodal SWCC are presented.     

考虑饱和-非饱和渗流的土坡极限分析

给出了均质土坝的下游坝坡的安全系数的计算方法。均质土坝的下游坝坡有可能会通过滑裂面发生破坏。土坝中的土坡通常处于非饱和状态。非饱和土坡的安全系数计算需要考虑吸力对抗剪强度的贡献以及土坡中的非饱和渗流。给出处于饱和-非饱和渗流状态下的土坡的安全系数,有助于评价均质土坝的安全系数。下游坝坡的安全系数计算方法有：极限平衡法、上限解法和下限解法,该算法适用于非饱和土坡,且是在饱和土坡安全系数的计算方法上修正得到的。算例中非饱和土坡安全系数的计算考虑了吸力对抗剪强度的贡献。考虑非饱和渗流理论的土坡安全系数计算方法通常更加接近现场实际情况,并且对于同一坝坡,考虑非饱和渗流计算出的土坡安全系数要比饱和渗流理论计算出的安全系数大。     

降雨入渗条件下非饱和土边坡稳定分析

针对降雨入渗土坡的稳定问题,建立一个考虑水力渗透系数特征曲线、土-水特征曲线以及修正的Mohr-Coulomb破坏准则的非饱和土流固耦合有限元计算模型,进行雨水入渗下非饱和土边坡渗流场和应力场耦合的数值模拟,得到非饱和土边坡变形与应力的若干重要规律。研究成果为降雨入渗条件下非饱和土边坡的稳定分析提供了基础。     

土-水特征曲线的试验研究及曲线拟合

土-水特征曲线是非饱和土力学研究中的一条重要曲线,其拟合参数对于分析非饱和土的各种性质具有重要意义。基于渗析法及滤纸法试验,测定了合肥市某工地非饱和膨胀土的土-水特征。以van Genuchten模型为基础,采用科学计算语言Matlab中的Lsqcurvefit函数对试验结果进行了曲线拟合,分析了土-水特征曲线拟合参数的多解性及其影响因素。拟合结果表明：渗析法试验得到的吸力范围偏低,拟合参数多解性较为明显;综合2种试验结果,可以得到较大吸力范围内的土-水特征曲线,且参数多解性较小;参与土-水特征曲线拟合的试验数据点的吸力范围应足够大,试验数据点应能够代表整条曲线的形状,试验点不必太过于密集,否则土-水特征曲线拟合参数的多解性较为明显;另一方面,只要试验数据点能代表整条曲线的形状,试验数据点的个数不必太多。这一研究结果对提高土-水特征曲线拟合参数的精度、提高试验效率有重要指导意义。     

Effect of soil hydraulic properties model on slope stability analysis based on strength reduction method

Soil hydraulic properties models which have been proposed were derived based on the empirical fitting curve such as Brooks-Corey model (BC) and Van Genuchten model (VG), or based on soil pore radius distribution such as Lognormal model (LN). Each model has different accuracy for predicting soil moisture distribution. In the analysis of rainfall-induced slope failure, the soil hydraulic properties model was needed to describe the physical phenomena of behavior characteristic of water in unsaturated soil. As moisture content has an effect on soil strength, it is vital to select the suitable soil hydraulic properties model for predicting Factor of Safety (FOS) especially in forecasting landslide hazard. In this study, a numerical model of seepage finite element analysis using BC, VG, and LN model were used and compared in order to analyze the soil moisture distribution, water movement phenomenon, and slope stability characteristic in unsaturated soil slope based on the strength reduction method (SRM). The results showed that the variations of the parameters predicting the moisture content of soil leads to differences of FOS in some cases. The parametric study showed that for the unsaturated soil condition, BC model has the greatest FOS value than the other model, while VG model has the lowest. On the other hand, the FOS of all models have the same result for the saturated condition. Other than that, it was found that the increasing of ESP value in the surface layer has significant effect in the sub-surface layer.     

Strength reduction FEM in stability analysis of soil slopes subjected to transient unsaturated seepage

The instability of soil slopes induced by the fluctuation of water level or rainfall has received much attention in the literature recently and the failure mechanism is closely related to the change of matric suction of unsaturated soils. Such a change is basically induced by unsaturated transient seepage. The traditional approach for estimating the stability of slopes subjected to unsaturated transient seepage is based on the limit equilibrium method (LEM). The limit equilibrium approach is limited by assumptions about analysis method and failure mechanism. In order to overcome those limitations, a finite element method with shear strength reduction technique (SSRFEM) has been successfully applied to the slope stability analysis in absence of seepage. The main objective of this paper is to extend the use of the strength reduction FEM to include the effects of unsaturated transient seepage and some primary numerical results concerning the stability of an earth dam under rapid drawdown are presented. Emphasis has been given to comparison of the safety factors obtained by LEM and SSRFEM. Topics for the further research in this area are also suggested.     

土-水特征曲线及其相关性研究

土-水特征曲线（SWCC）描述了土体吸力与含水率之间的关系,是非饱和土力学研究的重要内容。从理论和试验两个方面研究了土-水特征曲线和渗透曲线的相互关系。在理论上,应用自然比例法则研究了土-水特征曲线和渗透曲线的特征,从物理和力学意义两个方面探讨了土-水特征曲线及其导数、渗透曲线、颗粒分布曲线之间的特征,给出了土-水特征曲线和渗透曲线两指数函数表达式幂的取值范围和相互关系。试验研究结果表明：土体的土-水特征曲线受到土体的物质成分、塑性指数、孔隙结构、应力状态等多种因素影响,土体持水特性研究对工程建设具有重要的指导意义。     

水位下降条件下黏性土边坡稳定性的图表法

库（河）岸边坡由于坡体迎水面水位下降经常造成坡体稳定性降低,建立快速评估水位变化条件下的边坡稳定性分析方法具有重要的工程实用价值。基于库（河）水位下降过程及坡体内非稳定渗流条件,通过对ABAQUS程序进行二次开发,发展了考虑地表水-地下水联动作用下的黏性土边坡强度折减有限元分析方法;在分析边坡土性参数、相对渗透比值、边界条件对边坡稳定性影响的基础上,建立了考虑土性参数、相对渗透比值、库水位下降比、坡角等一体化的相对稳定安全系数综合图表表示方法。该方法能够简便、快捷地查出实际涉水边坡在不同工况下的稳定安全系数及设计边坡坡比,可作为现有图表法的有益补充。