非饱和土拟二维平面应变固结问题的解析计算方法

基于Fredlund非饱和土一维固结理论,建立了二维平面应变条件下的固结方程组,并得到了单层非饱和土平面应变条件下的解析解。首先,基于相关理论,假设体变系数和渗透系数都为常量,同时考虑到瞬时加载条件下,沿着土体深度方向上产生均匀或者线性分布的初始超孔隙压力,建立了二阶二元偏微分方程组。求解时,引入函数方法来降低方程的阶数,然后通过分离变量法获得方程的通解。在此基础上,结合一个针对单面排水条件下二维平面应变问题的案例,通过与数值解对比验证本文方法的正确性。并采用本文方法计算获得二维平面下超孔隙水、气压力沿垂直和水平方向消散的等时线,通过计算对比,分析了不同线性分布情况下,初始超孔隙压力对固结消散过程的影响。研究结果表明：初始超孔隙压力的不同分布对超孔隙气压力消散影响几乎可以忽略,而对超孔隙水压力的消散影响更大。     

基于COMSOL Multiphysics的非饱和裂隙土降雨入渗特性研究

基于COMSOL Multiphysics软件对非饱和裂隙土降雨入渗特性进行数值模拟研究。通过将裂隙和基质分别离散成有限单元,建立了能充分模拟土中裂隙流、基质流以及裂隙-基质流量交换的离散裂隙-孔隙介质模型。结合"空气单元"的概念,对裂隙土的上边界进行模拟。该方法不仅能描述降雨初期雨水沿裂隙优先入渗的现象,还能描述当降雨量大于裂隙土入渗量时雨水沿地表流走的现象。通过对地表以下2 m深度内低渗含裂隙土体进行模拟,分析了裂隙的几何特征、基质的水力特性、前期水分条件以及降雨强度对非饱和裂隙土降雨入渗过程的影响。结果表明,在非饱和裂隙土中,存在两个主要的渗流过程:一是水沿裂隙优先流动;二是水不断从裂隙吸入基质中,基质吸收水的作用抑制了裂隙中优势流的发展。与裂隙的几何特征相比,基质的水力特性对非饱和裂隙土渗流的影响较大。增大基质的饱和渗透系数可能使由裂隙流主导的渗流过程转变为由基质流主导的渗流过程,而基质的非饱和特性与裂隙土的初始含水率改变了土体的储水能力,从而加速或延缓了降雨入渗至某一深度的时间。降雨强度对土体入渗速率和入渗量均有影响,当超过裂隙土的入渗能力时,多余积水沿地表流走,断面入渗率随...     

非饱和土拟二维平面应变固结问题的解析计算方法

基于Fredlund非饱和土一维固结理论,建立了二维平面应变条件下的固结方程组,并得到了单层非饱和土平面应变条件下的解析解。基于相关理论,假设体变系数和渗透系数都为常量,同时考虑到瞬时加载条件下,沿着土体深度方向上产生均匀或者线性分布的初始超孔隙压力,建立了二阶二元偏微分方程组。求解时,引入函数方法来降低方程的阶数,然后通过分离变量法获得方程的通解。在此基础上,结合一个针对单面排水条件下二维平面应变问题案例,通过与数值解对比,验证了所提方法的正确性。并采用所提方法计算获得了二维平面下超孔隙水压力、气压力沿垂直和水平方向消散的等时线,通过计算对比,分析了不同线性分布情况下,初始超孔隙压力对固结消散过程的影响。研究结果表明：初始超孔隙压力的不同分布对超孔隙气压力消散的影响几乎可以忽略,而对超孔隙水压力消散的影响更大。     

成层非饱和土渗流的耦合解析解

成层土在工程中很常见,研究降雨过程中成层非饱和土的渗流-变形耦合对非饱和土土力学的发展具有重要的意义。由流体质量守恒,Darcy定律和Lloret等的非饱和土本构模型可得成层非饱和土渗流-变形耦合的控制方程。采用Gardner的非饱和土的渗透系数公式以及Boltzman模型,基于Laplace变换得到耦合方程的解析解。解析及其参数分析表明,渗流和变形耦合是具有时间效应的。与吸力变化相关的土的模量F,对成层土的孔隙水压力分布有明显影响。两层土的F差异越大,孔隙水压力消散得越慢,耦合效应越不显著。增大表层土的F值有利于降低耦合效应。成层土饱和体积含水率变化对吸力变化产生有限的影响     

真空预压加固区硬壳层的水分运移解析解

本文根据非饱和土体的渗流理论,建立水分运移模型,研究真空预压在稳态下的水分运移机制,给出了孔隙水压力分布的控制方程及其解析解。根据工程实测数据验证解析解的结果,并对不同参数进行分析。结果表明：稳态下真空预压水分运移模型计算的结果和现场实测数据基本吻合。扩散通量Q和降饱和系数?对上层非饱和区域的形成位置、强度及厚度的影响较大。扩散通量Q越大,土体表层的基质吸力越大,孔隙水压力变化趋势就越大;降饱和系数?在0.3~1.0 kPa-1内对非饱和土体部分基本没有影响。     

混合非齐次边界下非饱和土一维固结解析解

非饱和土的固结研究对道路工程、软土地基处理工程等具有十分重要的意义。基于Fredlund和Hasan提出的非饱和土一维固结理论,给出了土体内孔隙水压力和孔隙气压力变化的控制方程,并给出了单层非饱和土的初始条件与一类随时间变化的混合非齐次边界条件,构成了非饱和土一维固结的定解问题。通过非齐次边界条件齐次化和特征函数展开法,得到了土体内孔隙水压力和孔隙气压力消散的精确时域解析解。最后,通过对比验证了解析方法的正确性,并分析了边界条件指数变化对非饱和土体内孔隙水压力和孔隙气压力的消散以及土体沉降的影响。结果表明,边界处孔压或孔压梯度随时间的指数变化对非饱和土固结过程有明显影响。     

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

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

福建德化霞碧滑坡渗流稳定性模拟分析

降雨入渗导致土体孔隙水压力的增大造成土体有效抗剪强度的降低,进而发生滑坡。本文建立滑坡饱和-非饱和渗流有限元模型,考虑孔隙水压力随时间和空间任意不规则分布,计算特定降雨量、降雨强度、降雨历时不同滞后时段基于渗流场的滑坡稳定性。     

非饱和土入渗的数值模拟

入渗引起土中孔隙水压力的变化,是影响非饱和土性状的最重要的原因之一。因此,正确地分析和模拟非饱和土入渗过程具有重要意义。通过对一维和二维降雨入渗问题的研究,展示了湿润锋面、地下水位面及孔隙水压力分布随时间的变化,使得对入渗问题的认识更加直观和清楚。     

考虑土体结构性损伤的沉桩挤土效应解析分析

根据土塑性力学基本原理,采用简化的线性软化模型并考虑土的结构性损伤影响,推导出沉桩时产生的挤压应力的空间解析解,并在此基础上分析了土体具有不同的损伤程度时沉桩产生的挤压应力和极限扩张压力的空间分布规律。同时,采用Henkel孔隙水压力公式,进一步分析孔隙水压力的发展规律,该解析解可以弥补经典圆孔扩张理论在沉桩应用中的不足,能更好地揭示沉桩时的挤土机理。     

非饱和土一维固结混合流体方法的解析分析

对采用混合可压缩流体方法分析非饱和土一维固结问题的固结方程进行了求解,在得到的解析解的基础上,对影响非饱和土一维固结的因素进行了分析。分析结果表明,在采用混合流体方法计算非饱和土一维固结的孔隙水压力时,所用公式与计算饱和土一维固结的太沙基理论公式基本相同,不同之处在于引入Bishop有效应力系数来体现孔隙气对孔隙水的影响。而在非饱和土孔隙气压的计算公式中除了体现孔隙水对孔隙气的影响参数以外,还有体现孔隙气体的可压缩性对固结影响的参数。在所有影响因素中,影响非饱和土一维固结最重要的因素是孔隙流体的渗流路径。     

降雨诱发边坡破坏数值模拟两个关键问题的解决方法

研究降雨因素对边坡稳定性的影响,最有效的方法是室内外试验和数值模拟方法,最为重要的是设置与实际情况接近的降雨入渗边界与动态流量边界下初始孔隙水压力分布问题。为有效地解决这两方面的问题,首先将现有的降雨数值边界施加方法所得模拟结果与室内试验结果进行了对比分析,指出了传统方法的不足之处,并在此基础上提出了新的降雨数值边界条件,即在土体表面施加一层很薄的"空气单元",表述的是一个动态的边界条件。降雨入渗边界施加在"空气单元"表面,而不直接施加在土体单元表面,应用该方法完成的入渗土柱数值模拟结果更加接近现场和试验测得的降雨入渗率,为有效模拟降雨诱发边坡破坏的边界条件提供准确、真实可行的方法。同时,由于相对渗透系数和负孔隙水压力均为饱和度的函数,通过设置流体的抗拉强度来赋予节点负孔隙水压力模拟基质吸力,在数值模拟中每一个计算步根据节点饱和度自动更新负孔隙水压力和相对渗透系数,可实现FLAC软件中非饱和土的有效渗流分析。此外,为解决非饱和区土体施加负孔隙水压力所带来的问题,提出数值模拟过程中施加动态的流量边界,通过模拟自然界中地下水的补给,很好地解决使用FLAC进行数值模拟时动态流量边界条件下初始孔隙水压力的分布问题。     

复杂水力路径下非饱和流动的数值分析

多孔介质中的流动问题,与孔隙介质的特性,含水量状态以及含水量的变化历史密切相关。基于毛细循环滞回理论模型,考虑含水量变化历史对土水特征关系的影响,在开发的U-DYSAC2有限元程序中进行了相应的数值实施。在试验给定的初边值条件下进行了非饱和渗流模拟分析,并将模拟结果与实测数据比较,表明在压力边界条件反复变化下,考虑滞回效应能获得更接近实测的结果,证实该模型在模拟各种循环变化条件下非饱和土渗流初边值问题的适用性与必要性。对入渗重分布反复变化条件下非饱和土柱流动的数值模拟表明,考虑滞回与不考虑滞回条件下,含水量、孔隙水压力和湿峰的迁移的预测在入渗后的重分布过程差异较大。考虑滞回效应时,土柱上部的脱湿速率、下部的吸湿速率比不考虑滞回时要低。从而证实了非饱和多孔介质中的土水状态依赖于含水量变化,而且强烈依赖于土体的水力路径变化。因此,循环边界条件变化下,毛细滞回效应在非饱和渗流模拟中的影响显著,必须加以考虑。     

降雨入渗条件下非饱和膨胀土边坡原位监测

为了对降雨诱发的非饱和膨胀土边坡失稳的机理有较深入的了解,在湖北枣阳选取了一个11 m高的典型的非饱和膨胀土挖方边坡进行人工降雨模拟试验和原位综合监测。监测成果表明：降雨入渗造成2 m深度以内土层中孔隙水压力和含水量大幅度增加,致使膨胀土体的抗剪强度由于有效应力的减少及土体吸水膨胀软化而降低;同时,降雨入渗造成土体中水平应力与竖向应力比显著增加,并接近理论的极限状态应力比,以致软化的土体有可能沿着裂隙面发生局部被动破坏,此破裂面在一定条件下（如持续降雨条件下）可能会逐渐扩展,最后发展成为膨胀土中常见的渐进式滑坡。     

乐坝火车站后山边坡降雨作用下渗流数值分析

研究拉张裂缝充水和降雨共同条件下边坡渗流特征对边坡稳定性分析具有关键作用。基于饱和—非饱和土渗流理论及边坡饱和—非饱和土渗透特性,设计4种降雨类型并对案例边坡有限元渗流数值计算分析,分析表明:(1)在稳态渗流下,边坡上部负孔隙水压力值最大,边坡下部值最小,渗流速度在边坡上部、中上部、中下部值增加缓慢,在下部值骤增;(2)在瞬态渗流下,考虑拉张裂缝已充满水,坡体孔隙水压力由正值变为负值,起初边坡上部正孔隙水压力最小,下部正孔隙水压力最大,随时间推移,上部负孔隙水压力最大,下部负孔隙水压力最小,渗流速度边坡上部最小,边坡下部最大。因此,坡体不同位置渗流速度和孔隙水压力值存在差异,这对边坡治理极其重要。     

降雨诱发边坡破坏数值模拟两个关键问题的解决方法

研究降雨因素对边坡稳定性的影响,最有效的方法是室内外实验和数值模拟方法,最为重要的是设置与实际情况接近的降雨入渗边界与动态流量边界下初始孔隙水压力分布问题。为有效地解决这两方面的问题,首先将现有的降雨数值边界施加方法所得模拟结果与室内试验结果进行了对比分析,指出了传统方法的不足之处,并在此基础上提出了新的降雨数值边界条件,即在土体表面施加一层很薄的“空气单元”,表述的是一个动态的边界条件。降雨入渗边界施加在“空气单元”表面,而不直接施加在土体单元表面,应用该方法完成的入渗土柱数值模拟结果更加接近现场和试验测得的降雨入渗率,为有效模拟降雨诱发边坡破坏的边界条件提供准确、真实可行的方法。同时,由于相对渗透系数和负孔隙水压力均为饱和度的函数,通过设置流体的抗拉强度来赋予节点负孔隙水压力模拟基质吸力,在数值模拟中每一个计算步根据节点饱和度自动更新负孔隙水压力和相对渗透系数,可实现FLAC软件中非饱和土的有效渗流分析。此外,为解决非饱和区土体施加负孔隙水压力所带来的问题,本文提出数值模拟过程中施加动态的流量边界,通过模拟自然界中地下水的补给,很好地解决使用FLAC进行数值模拟时动态流量边界条件下初始孔隙水压力的分布问题。     

降雨诱发边坡破坏数值模拟两个关键问题的解决方法

研究降雨因素对边坡稳定性的影响,最有效的方法是室内外实验和数值模拟方法,最为重要的是设置与实际情况接近的降雨入渗边界与动态流量边界下初始孔隙水压力分布问题。为有效地解决这两方面的问题,首先将现有的降雨数值边界施加方法所得模拟结果与室内试验结果进行了对比分析,指出了传统方法的不足之处,并在此基础上提出了新的降雨数值边界条件,即在土体表面施加一层很薄的“空气单元”,表述的是一个动态的边界条件。降雨入渗边界施加在“空气单元”表面,而不直接施加在土体单元表面,应用该方法完成的入渗土柱数值模拟结果更加接近现场和试验测得的降雨入渗率,为有效模拟降雨诱发边坡破坏的边界条件提供准确、真实可行的方法。同时,由于相对渗透系数和负孔隙水压力均为饱和度的函数,通过设置流体的抗拉强度来赋予节点负孔隙水压力模拟基质吸力,在数值模拟中每一个计算步根据节点饱和度自动更新负孔隙水压力和相对渗透系数,可实现FLAC软件中非饱和土的有效渗流分析。此外,为解决非饱和区土体施加负孔隙水压力所带来的问题,本文提出数值模拟过程中施加动态的流量边界,通过模拟自然界中地下水的补给,很好地解决使用FLAC进行数值模拟时动态流量边界条件下初始孔隙水压力的分布问题。     

高压环境下流体非饱和渗流对含水合物沉积物的影响

沉积物与流体流动的性质是影响水合物形成和聚集的两个重要因素,为研究水合物在沉积地层中的赋存机制必须探明高压环境下含水合物沉积物在非饱和渗流条件下的相互影响关系。以逸度差为水合物反应驱动力,反应动力学常数为Arrhenius类型,建立了包括非饱和流体流动-沉积物特征-水合物形成动力学耦合的二维模型,从理论上研究了多孔介质内流体与沉积物参数如含水率、去饱和系数、水力分布和水合物饱和度等在孔隙内的相互影响规律。结果表明,在设定的条件下,随着反应的进行孔隙水压力随时间逐渐大,在相同条件下水合物饱和度与温度增加导致孔隙水压力变大,其中水合物饱和度的影响较小,而沉积物基质吸力、去饱和系数与本征动力学常数则与孔隙水压力成反向变化,其中本征动力学常数的影响较大。     

三峡库区降雨型滑坡入渗特征及变形机制——基于一维和二维模型试验研究

强降雨诱发古滑坡的浅层变形是三峡库区最为严重的地质灾害,因此,探索暴雨作用下滑坡土体的入渗规律及其浅层变形机制具有重要的研究意义。以三峡库区的降雨型滑坡为研究对象,统计了降雨型滑坡土体的渗透特性分布特征,考虑强降雨作用,分别开展一维土柱入渗试验和二维滑坡模型试验,研究了不同降雨强度条件下滑坡土体的入渗特征以及对应的滑坡浅层变形机制。降雨入渗试验结果表明：降雨入渗土体的速度取决于降雨强度与土体渗透系数的相对大小,当降雨强度小于或等于土体渗透系数时,入渗能力随降雨强度增大而增大。当降雨强度大于土体渗透系数时,入渗能力反而随之减小。模型试验结果表明：强降雨入渗时,会造成地表土体暂态饱和,表层土体以下非饱和区内的气体被暂态封闭并受到表层孔隙水压力作用而压缩,导致孔隙气压力随降雨入渗迅速增大。对于三峡库区的降雨型滑坡,短时的暴雨会产生瞬态饱和区。封闭气体,这也是影响强降雨入渗能力的主要原因;同时封闭气体传递水压使得浅层土体的孔隙水压力骤增,这也是许多滑坡发生浅层变形和破坏的主要原因。     

一维稳态入渗时非饱和土基质吸力分布

雨水的入渗和蒸发对土体的其质吸力分布有着重要的影响作用。利用水相连续性条件、广义Darcy定律及Gardner关于非饱和土导水系数和基质吸力的经验公式,建立并求解出非饱和土一维稳态渗流情况下基质吸力沿深度变化的关系表达式。根据所解得的解析解,详细分析了在单一土层情况下参数α、入渗强度q和qlks对基质吸力分布的影响。