非饱和土力学吸力与.pptx
刘 艳北 京 交 通 大 学 土 建 学 院 岩 土 工 程 系E m a i l：非饱和土力学理论简介第1页/共75页三.吸力与土水特征曲线1.吸力的概念吸力的概念2.土水特征曲线土水特征曲线3.滞后现象滞后现象4.SWCC模型模型5.SWCC影响因素影响因素第2页/共75页1.吸力20世纪初，土中吸力的概念已在土壤物理学中得到发展。土的吸力理论主要是同土-水-植物相关连而发展起来的。1956年召开的力学研讨会上，Aitchison的一篇论文从热动力学角度对土的吸力及其组成下来定义。这些定义已在岩土工程学中被广泛接受。土中吸力反映土中水的自由能状态1.吸力的概念第3页/共75页自由能什么是自由能？freeenergy1.吸力的概念 在热力学当中，自由能指的是在某一个热力学过程中，系统减少的内能中可以转化为对外做功的部分，它衡量的是：在一个特定的热力学过程中，系统可对外输出的“有用能量”。通常有 Helmholtz自由能和Gibbs自由能lHelmholtz自由能是等温过程中系统能对外做功的那部分能量lGibbs自由能是恒温恒压过程中，封闭系统由体积功以外其它广义力所做功的那部分能量第4页/共75页热力学第一定律能量守恒热力学第二定律熵增原理自由能的变化1.吸力的概念Gibbs自由能增量第5页/共75页孔隙水总势能1.吸力的概念产生土体总吸力的物理与化学作用机理就是相对于自由水状态来说，土中孔隙水势能的减少量以土中水作为研究对象，系统的非体积功主要来源有：重力水中溶质的渗透力孔隙水的持水力温度压力物理化学力土中水的总势能可以表示为第6页/共75页孔隙水总势能1.吸力的概念第7页/共75页孔隙水总势能单位化学势能反映每单位质量所包含的能量，J/mol或J/kg；压力势能反映每单位体积所包含的能量，J/m3=N/m2=Pa；水头势能反映每单位重量所包含的能量，J/N=m。势能、压力和水头的单位间具有如下关系1.吸力的概念第8页/共75页忽略温度、重力与惯性的影响，促使土体孔隙水势能降低的主要因素有：l渗透作用孔隙水中溶质溶解的结果l毛细作用水-气交界面曲率以及负孔隙水压力l吸附作用固-液（即土中孔隙水）交界面附近的电场与范德华力场作用而产生孔隙水总势能1.吸力的概念第9页/共75页孔隙水总势能1.吸力的概念osmoticelectricalcapillaryvanderWaals第10页/共75页孔隙水势能渗透作用1.吸力的概念第11页/共75页孔隙水势能渗透作用1.吸力的概念对更普遍的、非稀释溶液的情况，渗透压力可以表示为（Shaw，1992）第12页/共75页孔隙水势能渗透作用1.吸力的概念随着溶液溶质浓度的增加，渗透压力也逐渐增加，相应的孔隙溶液的化学势能却降低当土中有化学溶液输运时，土中孔隙水的化学溶度发生改变，此时渗透吸力对土的性质具有较大的影响。然而就其它情况而言可忽略吸力中的渗透部分。第13页/共75页孔隙水势能毛细作用1.吸力的概念随着交界面曲率的增大（意味着含水量降低，负孔隙水压力的数值变得更大），化学势能会显著降低第14页/共75页孔隙水势能其它势能1.吸力的概念以上两式表明：当水分子相对远离土颗粒表面时，电场和范德华力引起的化学势能的降低值要远小于水分子靠近土颗粒表面是化学势能的降低值。上述每种作用机理均造成了化学势能的下降，这些下降势能的总和就定义为土水系统的总吸力第15页/共75页土体的总吸力1.吸力的概念Suction孔隙水与土颗粒间相互作用而产生的吸力中的各个部分因溶质溶解作用而产生的吸力毛细作用吸附作用渗透作用第16页/共75页基质吸力的探讨 Matric suction源自毛细与短程吸附综合作用的吸力通常定义为基质吸力，它具有与压力一样的单位。基质，是指细小的土颗粒。基质吸力可以看作是土基质对水的吸持潜能。土基质对水的吸持潜能与土体的含水率有关。由于饱和土体所有孔隙都已经被水占据，因此没有吸水的能力了，所以，我们认为饱和土的基质吸力等于0。而非饱和土的孔隙并没有被全部占据，所以具有吸水的潜能，基质吸力大于0。将基质吸力引入到非饱和土及土体变形强度稳定的研究与分析中去是当前非饱和力学研究发展的一条基本线索1.吸力的概念第17页/共75页基质吸力的探讨土基质吸持水分的机理十分复杂，但可概括为吸附作用和毛细作用。Derjaguin等（1987）和Tuller等（1999）强调基质吸力由两部分组成，可表示为：单纯使用毛细模型来解释和表达基质吸力具有一定的局限性，但是由于毛细作用的机理比较清楚，数学处理也比较便利，所以在非饱和土力学的定量研究中常常还是将基质吸力作用等同于毛细效应。通常，基质吸力中的毛细部分可表示为：1.吸力的概念第18页/共75页基质吸力的探讨1.吸力的概念理想化的基质吸力分布 处于地下水位以上的非饱和土体的基质吸力大小与土体的深度呈现一定的关系。假定土壤地表基质吸力为So,地下水位处基质吸力为0,基质吸力从地表至地下水位处线性减小。第19页/共75页基质吸力的探讨基质吸力中的毛细部分和粘吸部分在概念上的区分是明显的，但难以通过试验手段加以区分基质吸力中的各个部分对非饱和土行为影响的机理并不相同。对于低塑性的或较高含水量下的土体，基质吸力中毛细部分占支配地位；然而对于高塑性的粘土或较低含水量下的土体，基质吸力中粘吸部分占支配地位。当基质吸力很大时，只表明液相吸附到固相的程度很大，但绝不能认为是传统意义上的负孔隙压力很大。然而，在以往的研究中，通常忽略基质吸力中粘吸部分的作用，认为基质吸力仅由毛细作用产生，致使基质吸力很大时，认为负孔隙水压力亦很大。目前针对非饱和土所建立的本构模型多是基于毛细机理，因此这些模型对高塑性的粘土或低含水量下的土体的适用性是值得探讨的。1.吸力的概念第20页/共75页三.吸力与土水特征曲线1.吸力的概念吸力的概念2.土水特征曲线土水特征曲线3.滞后现象滞后现象4.SWCC模型模型5.SWCC影响因素影响因素第21页/共75页土水特征曲线2.土水特征曲线 基质吸力与土的含基质吸力与土的含水率有关，它与含水率水率有关，它与含水率之间的关系曲线称为土之间的关系曲线称为土水特征曲线水特征曲线。Soil Water Characteristic Curve(SWCC)Water Retention Curve(WRC)基质吸力随着含水量的增大而减小 第22页/共75页土水特征曲线2.土水特征曲线可以是含水量或饱和度第23页/共75页体积含水量含水量含水量饱和度饱和度WaterAirSolidVaVwVsVvVma=0mwmsm质量质量体积体积表达式表达式:定义定义:土中水的体积与孔隙体土中水的体积与孔隙体积的比值积的比值饱和度表示饱和度表示孔隙中充满水的程度孔隙中充满水的程度Sr=0 :干土干土Sr=1 :饱和饱和土土0Sr1:非饱和土非饱和土基本变量含水量与饱和度2.土水特征曲线定义定义:土中水的质量与土粒质量之比土中水的质量与土粒质量之比,用用百分数表示百分数表示表达式表达式:第24页/共75页基本变量进气值2.土水特征曲线第25页/共75页有效饱和度有效饱和度基本变量有效饱和度2.土水特征曲线Effective degree of saturation也可以用饱和度表示第26页/共75页土水特征曲线2.土水特征曲线第27页/共75页土水特征曲线用于描述吸力与含水量之间本构关系的函数曲线2.土水特征曲线低含水量孔隙水结合水形式存在孔隙水势能较自由水低高含水量孔隙水以毛细形式存在孔隙水势能与自由水间差值相对较小第28页/共75页土水特征曲线2.土水特征曲线第29页/共75页土水特征曲线2.土水特征曲线理想化的土水特征曲线第30页/共75页三.吸力与土水特征曲线1.吸力的概念吸力的概念2.土水特征曲线土水特征曲线3.滞后现象滞后现象4.SWCC模型模型5.SWCC影响因素影响因素第31页/共75页滞后现象3.滞后现象第32页/共75页滞后现象3.滞后现象第33页/共75页不同类型土的土水特征曲线3.滞后现象第34页/共75页滞后原因Fredlund(2000)总结了SWCC出现的滞后效应的原因，主要包括:(1)孔隙尺寸分布不均匀。l在湿化过程中，水将首先进入湿锋附近的小孔隙，并将其充满，然后再充满大孔隙。这是因为在小孔隙中的孔隙水具有最低的化学势（最稳定），而在大孔隙中孔隙水化学势较高。l在干燥过程中，位于大孔隙中的孔隙水首先排出来，然后再轮到小孔隙排水，孔隙内的气体就有可能会沿着连通大孔隙形成连通的气流路径，从而阻隔了小孔隙的进一步排水，使得孔隙水在孔隙介质中呈块状分布。然而，在湿化过程中，由于小孔隙首先被充满，所以不会形成上述水流通路阻隔现象，使得孔隙水分布相对比较均匀。3.滞后现象第35页/共75页滞后原因(2)气体体积变化不同l当吸力增加或减少时孔隙中的气体的体积及其变化是不同的，并导致饱和度的变化也不同(3)触变和时间效应。3.滞后现象第36页/共75页滞后现象(4)瓶颈效应墨水瓶滞后作用l不同大小的孔隙，以及相互连通的孔隙喉道之间的尺寸差别造成了这种作用。在浸润过程中，由于孔隙以及与其连通的喉道之间存在着尺寸差异，孔隙水在涌入的过程中自然面临着瓶颈的“约束”而难以突破，导致在相同吸力下浸润时的含水量小于干燥时的含水量。5.土水特征曲线第37页/共75页滞后原因(5)接触角的影响。l在干燥与浸润过程中，水-气交界面上的接触角会有所不同。l一般干燥时接触角小，浸润时大；小的接触角对应的表面张力较大，因此对水的滞留能量较大。l接触角的大小差异决定了水的滞留特性的差别，这种现象称之为雨点效应。3.滞后现象第38页/共75页三.吸力与土水特征曲线1.吸力的概念吸力的概念2.土水特征曲线土水特征曲线3.滞后现象滞后现象4.SWCC模型模型5.SWCC影响因素影响因素第39页/共75页滞后现象3.滞后现象第40页/共75页SWCC模型目前，习惯用室内试验加理论模型的方式对SWCC的滞后性进行描述。对于考虑滞后效应的任意路径的SWCC，由于室内试验往往费时费力，进行完整的循环路径实验资料并不多。现有的滞后模型主要包括以下几种类型：l经验模型l域模型l理论外推模型l边界面模型参考文献：Pham H Q,Fredlund D G,Barbour S L.A study of hysteresis models for soil-water characteristic curves.Canadian Geotechnical Journal,2005,42(6):1548-15684.SWCC模型第41页/共75页SWCC模型4.SWCC模型(Pham,etal,2005)第42页/共75页SWCC模型4.SWCC模型(Pham,etal,2005)第43页/共75页SWCC模型经验模型l代表性模型：Scott等(1983)提出的比例缩放模型l有一些研究者在Scott模型基础之上做出了一些修正，如Kawai等(2000)，Karube等(2001)l实际上扫描曲线形状与边界面形状并不完全相符，因此该类模型的精确度不高，但是由于其简单适用，因此得到了一定的应用。4.SWCC模型曲线拟合公式第44页/共75页SWCC模型经验模型l代表性模型：Feng&Fredlund(1999)模型l此类模型主要是对浸润/干燥边界面的描述，缺乏对任意扫描线描述的功能。但是，该模型仅需少数的几个点即可得出整条曲线。另外，在简化的Feng&Fredlund模型中可以用一条边界曲线即可拟合另外一条边界线。Feng&Fredlund模型以及后续的简化模型拟合的精度非常高，而且所需标定的数据较少，因此在边界面的模拟中得到了广泛的应用。4.SWCC模型基于干燥/浸润边界的关系进行预测的经验模型第45页/共75页SWCC模型域模型域模型是一种将土视为孔隙的集合体，以每个孔隙的吸排水特性作为基本的研究单元，在统计学的基础上，通过引入孔隙水分布函数来计算土中含水量随吸力变化规律的土水特征曲线滞后模型。域模型本质上是一种利用边界滞回圈通过内插的方法计算扫描线的计算模型，早期的域模型在计算时除了需要实测两条边界曲线外，还需要一定数量的扫描线来标定参数，由于这些模型在计算时所需的实测数据较多，因此应用起来并不方便。5.土水特征曲线第46页/共75页SWCC模型域模型Mualem(1973)假定“孔隙水分布函数可表示为两个独立分布函数的乘积”，利用“相似性假定”简化后的域模型仅需实测两条边界曲线即可预测滞回圈中的扫描线。Mualem随后将他的相似性假定应用到了一系列毛细滞回循环模型中，既提高了域模型的计算精度，又在一定程度上简化了计算过程，使得域模型在工程中得到了一定的应用。域模型的优点是具备良好的理论基础，在一定程度上能够反映土水特征曲线滞后特性的物理本质；其不足之处在于，这类模型的计算过程，尤其是计算高阶扫描线时，在吸力变化历史未知的情况下确定扫描线过程十分复杂，因而限制了它在工程上的应用。5.土水特征曲线第47页/共75页SWCC模型理论外推模型Mualem模型暗示着浸润与干燥扫描线非常规则而又光滑的穿越区域边界，但是试验结果(Topp，1971表明：干燥边界曲线的斜率往往与扫描线的斜率不同。Parlange（1976）在Mualem的相似性假设的基础上，提出了理性外推模型，Hogarth等(1988)，Liu等(1995)发展了理性外推模型，使其能考虑含气量的大小。此类理性外推模型存在一些难以解决的缺陷：在含水量变化较小时对扫描曲线的描述往往比较准确，但是一旦含水量变化范围太大，或者浸润扫描线贴近于浸润边界线，这时模型的预测结果往往与实测结果有不小的差距。4.SWCC模型第48页/共75页SWCC模型边界面模型根据边界面塑性理论，可以利用加载面上的应力点与其在边界面上的映射点之间的距离来确定加载面上的塑性反应。基于这一理论，Li(2005)和Wei(2006)分别建立了模拟土水特征曲线滞回循环的计算模型。模型以边界浸润和干燥曲线作为计算的边界，以浸润-干燥的反弯点作为投影中心，建立了扫描线上的斜率与边界曲线斜率之间的关系。优缺点：能够计算高阶扫描线，每个模型中各有一个参数，标定参数时除了需要测量边界干燥和浸润曲线外，还都需实测一条一阶扫描曲线，因此利用这两个模型计算都需要提供较多的试验数据。4.SWCC模型第49页/共75页SWCC模型举例4.SWCC模型Fredlund D G,Xing A,1994.Equations for the soil-water characteristic curveJ.Canadian Geotechnical Journal,31(4):521-532.第50页/共75页SWCC模型举例4.SWCC模型VGmodelvanGenuchten(1980)FredlundandXing(1994)FengandFredlund（1999）第51页/共75页SWCC模型举例4.SWCC模型FredlundandXing(1994)第52页/共75页SWCC模型举例4.SWCC模型FredlundandXing(1994)第53页/共75页SWCC模型举例4.SWCC模型FredlundandXing(1994)第54页/共75页边界面模型4.SWCC模型Bounding Surface PlasticityRole of the bounding surface:“The bounding surface concept expresses macroscopically a bounding material state,created by the manufacturing process and/or the past loading history,in the sense that plastic response depends on how far the current stress state on loading surface is from a corresponding bounding or image stress state on bounding surface.”Dafalia(1986)pIn this case,the boundary surface is two endpoints of the bounding zone第55页/共75页边界面模型第56页/共75页边界面模型4.SWCC模型 The plastic modulus is related to the bounding plastic modulus via the distance between the stress and image stress by an expression of the general form:第57页/共75页边界面模型4.SWCC模型Formulation of the bounding curvesFeng and Fredlunds equationisusedtofitthemeasuredboundarycurves:第58页/共75页边界面模型4.SWCC模型Model formulationModel formulationModel formulationModel formulation 韦昌富（Wei，2004）李相菘（Li，2006）刘艳（2007）only only one one parameterparameter第59页/共75页第60页/共75页Topp(1971)Caribou Silt loamC=2000=4000=1.18第61页/共75页Topp(1971)Caribou Silt loamC=2000=4000=1.18第62页/共75页Poulovassilis(1970)Sand AC=250=500=2第63页/共75页Poulovassilis(1970)Sand AC=250=500=2第64页/共75页三.吸力与土水特征曲线1.吸力的概念吸力的概念2.土水特征曲线土水特征曲线3.滞后现象滞后现象4.SWCC模型模型5.SWCC影响因素影响因素第65页/共75页SWCC影响因素5.SWCC影响因素SWCC矿物成分密实程度温度孔隙结构应力状态第66页/共75页SWCC影响因素矿物成分土体矿物成分不同，土体的持水能力会不同，随着土中粘粒含量逐渐增多，土的进气值和残余体积含水量都逐渐变大，持水能力逐渐增强。5.SWCC影响因素Se10kPaSe=10kPa-100kPa砂土进气值可达几十至几百kPa粉质粘土粘性土对于确定的土样，矿物成分和温度影响可以不考虑第67页/共75页SWCC影响因素应力状态孙德安通过研究指出，SWCC与应力状态无直接关系，即使应力状态不同，只要孔隙比相近，其SWCC就相近。当然，这一结论是指某一特定土在变形过程中孔隙结构变化不大的条件下得到的。因此在孔隙结构变化不太大的土体变形过程中，可以用孔隙比变化表示孔隙结构和密实程度的变化。5.SWCC影响因素第68页/共75页SWCC影响因素孔隙结构和密实程度土体的孔隙结构对SWCC、其自身的变形，渗透系数都存在影响，土体变形将改变土的密实程度(表现为孔隙比的改变)，进而影响到土水特征曲线在“含水量-吸力”空间中的位置和形状。孔隙比的变化会改变土的进气值，从而使SWCC的位置移动；但孔隙分布指数并没有随着孔隙比的改变而单调变化，即孔隙分布指数随着孔隙比的变化没有唯一的关系，它如何变化依赖于土的类型。另外孔隙比对SWCC的影响程度还依赖于土的初始孔隙结构：孔隙比对于具有双孔结构的土的SWCC影响较大，而对于具有分散结构的土的SWCC影响则比较小。5.SWCC影响因素第69页/共75页5.SWCC影响因素SWCC影响因素不同密度时不同密度时SWCCSWCCDatafrom孙德安第70页/共75页Density dependency of water-retention curve10010110210311.5eABCDEpn=19.6kPae0=1.76w0=26%s(kPa)10010110210320406080100s(kPa)ABCDEe：大e：大e：小Sr(%)e：小SWCCNet stress：pn=cont.,Suction schange（来源：孙德安）第71页/共75页5.SWCC影响因素考虑孔隙比影响的SWCC方程第72页/共75页三维土水特征曲线5.SWCC影响因素TsiampousiA,ZdravkoviL,PottsDM,2013.Athree-dimensionalhystereticsoil-waterretentioncurveJ.Gotechnique,63(2):155-164.第73页/共75页5.SWCC影响因素考虑骨架变形影响的SWCC方程增量基本方程双线性形式的SWCC虽然可以简单地应用于非饱和土本构模型中，但是与实际情况比有较大的误差，因此可以用一些经验的SWCC方程来替代这里的边界面，比如用vanGenuchten的方程作为边界面方程，可以较好的拟合试验数据。第74页/共75页感谢您的观看！第75页/共75页
