土力学中的渗透力与超静孔隙水压力_李广信

1-11-2、在直剪、单剪和环剪试验中，试样的应力和应变有什么特点？（1）在直剪试验中，由于试验的破坏面是人为确定得，试样中得应力和应变不均匀且十分复杂，试样各点应力状态及应力路径不同。

在剪切面附近单元的主应力大小是变化的，方向是旋转的。

（2）在单剪试验中，仪器用一系列环形圈代替刚性盒，因而没有明显的应力应变不均匀，试样内所加的应力被认为是纯剪。

（3）环剪试验，试样是环状的，剪切沿着圆周方向旋转，所以剪切面的总面积不变，特别适用于量测大应变后土的残余强度和终极强度。

1-3、说明围压σc、土的平均粒径d50、土的级配和橡皮膜的厚度对于三轴排水试验的膜嵌入效应各有什么影响？（1）对于常规三轴压缩排水试验，由于其围压σ_c=σ_3是不变的，其对膜嵌入的影响很小。

但对于三轴不排水试验，其有效围压随孔压变化而变化，围压对膜嵌入影响较大。

一般来说，围压越大，膜嵌入越明显。

（2）土的平均有效粒径越大，则土越粗，一般而言，粗粒土膜嵌入明显，细粒土则相反。

（3）土的级配越好，膜嵌入越不明显，反之则相反。

（4）橡皮膜的越厚，膜嵌入越不明显，越薄则相反。

1-6、对土工栅格进行蠕变试验，120年后应变达到5%的荷载为70kN/m。

在n=100的土工离心机试验中，该栅格在70kN/m的荷载作用下，应变达到5%需要多少时间？土工离心机加速度am=100g，对于蠕变问题，时间的比尺因素为1，故在离心机上试验时，在相同荷载下，达到同样的应变时，其时间相等，为120年。

2-2说明土与金属材料的应力应变关系有什么主要区别（找了两种答案，自己开心就好）土体的应力应变关系主要特点是其非线性与非弹性。

如下图（没图），左边为金属材料，下图为土的材料。

金属材料开始时有一段直线。

而土体应力应变曲线显示出其很明显的非线性关系。

其应变很大一部分是塑性应变，而且土的变形为非弹性。

答：金属材料被视作线弹性材料，符合弹性力学中的五个假定：连续性、线弹性、均匀性、各向同性和微小变形假定，土体应力应变与金属材料完全不同，体现在以下几个方面：1）土体应力应变的非线性和弹塑性：金属材料的应力应变在各个阶段呈线性，在屈服强度以内呈弹性；而由于土体是由碎散的固体颗粒组成，其变形主要是由于颗粒间的错位引起，颗粒本身的变形不是主要因素，因此在不同应力水平下由相同的应力增量引起的变形增量不同，表现出应力应变关系的非线性。

第3章习题摩尔-库仑公式推导:ϕ+ϕσ+σ=σ-σcos c sin 223131 即: 231231]cos c 2sin )[()(ϕ+ϕσ+σ=σ-σ,同理有;232232]cos c 2sin )[()(ϕ+ϕσ+σ=σ-σ; 221221]cos c 2sin )[()(ϕ+ϕσ+σ=σ-σ破坏面条件:{}{}{}0]cos c 2sin )[()(]cos c 2sin )[()(]cos c 2sin )[()(221221232232231231=ϕ+ϕσ+σ=σ-σ⨯ϕ+ϕσ+σ=σ-σ⨯ϕ+ϕσ+σ=σ-σ⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧+⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧π-θ-θπ+θ=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧σσσ1112321I 31I 31I 31)6cos()sin()6cos(J 32 将该式代入上式得:0cos C J )3sin sin (cos sin I 3121=ϕ+ϕθ+θ-ϕ π平面上各轴的投影:在1σ轴上的投影:2S 2321321=σ-σ-σ在2σ轴上的投影:2S 2322312=σ-σ-σ在3σ轴上的投影:2S 2323213=σ-σ-σ如: 1σ=400kPa, 2σ=3σ=100kPa. 则在三个轴上的投影分别为: 141kPa, -71kPa, -71kPa.1、临界状态：是指土在常应力和常孔隙比下不断变形的状态。

临界孔隙比：表示土在这种密度状态下，受剪作用只产生剪应变而不产生体应变。

水力劈裂：由于孔隙水压力的升高，引起土体产生拉伸裂缝发生和发展的现象。

饱和松砂的流滑：饱和松砂在受静力剪切后，因体积收缩导致超孔压骤然升高，从而失去强度和流动的现象。

真强度理论：为了反映孔隙比对粘土抗剪强度及其指标的影响，将抗剪强度分为受孔隙比影响的粘聚分量与不受孔隙比影响的摩擦分量。

通过不同的固结历史，形成等孔隙比的试样，在不同的法向压力下剪切，试样破坏时的孔隙比相同，强度包线即为孔隙比相同的试样的强度包线，该强度称为在此孔隙比时的真强度。

2002年5月水 利 学 报SHUILI XUE BAO 第5期收稿日期:2001-12-28基金项目:国家自然科学基金委员会和长江水利委员会联合资助项目(50099620)作者简介:李广信(1941-),男,黑龙江呼兰人,博士,教授.主要从事土的本构、土工合成材料和基础工程方面的研究.文章编号:0559-9350(2002)05-0075-06渗透对基坑水土压力的影响李广信1,刘早云1,温庆博1(1.清华大学水利水电工程系　北京　100084)摘　要:基坑地基土中水的渗流不但可能引起渗透破坏,引起水压力,而且也对其土压力有重大影响,从而决定抗滑稳定性.本文作者针对有上层滞水、一般自由渗透、有承压水、基坑内排水与基坑外降水以及有超静孔压等情况对基坑支护结构物上的水土压力进行计算分析,结果表明:水土压力大小及分布与静水时的明显不同,且此时较宜于用库伦土压力理论.在有上层滞水情况下,用水土合算大体上是可以接受的.在有承压水情况下,其作为抗力的被动土压力可能丧失殆尽.基坑外人工降水与基坑内排水相比,更有利于基坑的稳定.正的超静孔压大大提高了土压力,负的超静孔压明显减少土压力.在很多有渗流的情况下,不宜用朗肯土压力计算土压力,而应当用库仑土压力理论的图解法来搜索可能滑裂面.关键词:基坑;水压力;土压力;渗流;渗透力中图分类号:TU476文献标识码:A随着我国大规模建筑基坑和地下工程的发展,支护结构设计计算中的许多问题逐步凸现出来.支护结构水土压力计算得到越来越多的重视和讨论[1～3].基坑水土压力计算常常采用朗肯和库伦土压力理论,其中朗肯理论由于其简便而被广泛使用.土中水的问题是土压力计算的难点,简单的水土合算与分算并不能解决实际工程中的复杂问题.土中水的存在状态有多种,而地下水存在的形态又有上层滞水、潜水和承压水[4].在基坑开挖过程中,基坑的水常常处于流动状态.由于朗肯理论要求墙后土应力状态为一维情况,这样,在一些土中水为静水压力或者水压力为一维渗流情况下适用,但有平面渗流的情况就不适用.库伦土压力理论由于考虑土楔体的极限平衡,因而更为适用在有渗流的情况下计算水土压力[11].当挡土墙墙后水为二维渗流时,由于渗透力方向不全是竖直方向的,故朗肯理论不适用.这时朗肯理论与库伦理论计算的结果有很大的不同.水土压力的分布还受不同土层渗透系数的影响,当土层的渗透系数由大到小,或者由小到大,考虑渗透影响,其水土压力分布有很大的不同,这点在北京某些地区含地下水的土中表现特别明显.此外,渗透力方向的不同,将影响基坑的水土压力.挡土结构后面的土中存在二维分布的超静孔压时,此时不宜用朗肯理论而应用库伦土压力理论.杨晓军、李广信等[5-8]对基坑有渗流情况的水土压力进行了一些分析,本文将通过对基坑的几种工况和不同方法的计算,对朗肯、库伦土压力理论在有渗流情况下的适用性进行分析和比较,对渗透系数不同的土层水土压力进行了计算,考虑到不同渗流方向,对同一个基坑采用基坑内外降水(轻型井点和明沟降水)对板桩墙水土压力的影响进行分析,并对挡土结构后面的土中存在二维分布的超静孔压的情况进行了计算和分析.1　一维垂直渗流的水土压力计算北京东部土层中往往有3层水,即潜水、滞水和承压水,由不同的土层组合而形成,基坑及基础常涉及到这3层水,由于不同土层组合,3层水的组合常常使水土压力的分布与静水时相比有所不同,有时甚至出现管涌流土现象.算例1,各土层的γsat均为18kN m3,摩擦角φ为30°,C值近似取为0,各土层垂直向下一维渗流,土层分布如图1所示,其中渗透系数k为相对值.考虑渗透力,土压力用朗肯土压力理论计算,计算基坑水土压力随深度变化如图2所示(图中实线为水压力,长划线为土压力,短划线为为水土总压力).在图2(a)中,粉土层中采用水土分算,但由于垂直下渗,计算结果与水土合算相同.从图2(b)可见,不同渗透系数土层的分布导致水土压力的不同分布.可见在这种上层滞水情况下,简单地水土分算,即用γ′计算土压力加上静水压力,和简单的水土合算都是不合适的.图1　土层的一维渗流情况(算例1)算例2,如图3一维渗流的情况,设饱和容重γsat均为18.5kg m3,摩擦角φ均为30°,粘土的粘聚力为5kPa.由于粘土上下均为渗透性较大的砂土,所以计算时可以考虑为沿板桩墙的垂直渗透.从图中可以看出,当上、下土层渗透系数如图所示,基坑的被动土压力将会显著的减小,抗力大大小于荷载,且被动侧很容易发生流土破坏((γ′-iγw)γw=0.533),这将对基坑的稳定极为不利.图2　不同土层的水土压力分布(算例1)图3　一维渗流水土压力分布2　二维渗流情况下基坑水土压力的计算算例3,一基坑挡土结构如图4所示,该土层为正常固结土,饱和容重γsat=18.5kN m3,粘聚力C=0,内摩擦角φ=30°,水头差■H=10m,基坑内排水,计算墙左侧每延米所受的总压力E左(E左=被动土压力P P+水压力E W1)和墙右侧所受的总压力E右(E右=主动土压力P a+水压力E w2).采用如下几种方法计算:(1)用朗肯土压力理论,忽略板桩与土间摩擦力.假设沿板桩墙水头均匀损失,即假定沿板桩墙背各点的水力坡降相等.以这种简化方法计算挡土墙所受的力.(2)用朗肯理论计算,但是渗透力用流网来计算.用有限元程序计算基坑的等势线图,再手动绘制流线进行计算.(3)不考虑板桩墙与土间的摩擦力,用水土合算的方法,即直接取饱和容重用朗肯土压力理论计算作用于板桩墙上的总压力.(4)用库伦理论对板桩墙进行滑裂面计算.墙右侧土压力E右的计算如图5.取土水混合体作为滑动的楔体,楔体自重用饱和容重计算.取直线滑裂面,在滑裂面上除作用着土体支承反力R外,还有垂直于滑裂面上的水压力P W,P W根据滑裂面切割到的流网网格叠加进行计算,其大小随滑裂面的位置变化而变化.板桩墙对土楔体的作用力E包括水压力E w和土压力P a (或P p).另一方法是滑动楔体取其骨架(浮容重)为隔离体,此时需根据一个个有限元小块计算渗流力,进行叠加计算,这将非常繁琐,结果应当是一致的.图4　板桩墙流网示意(算例3)图5　板桩墙主动土压力计算(算例3)上述4种计算方法计算结果列于表1,其中方法4(2)和方法4(1)的板桩墙与土间的摩擦角分别为20°、0°,总压力为水土压力的矢量合成,取方法4(2)为比较基准.表1　板桩墙土压力计算结果(算例3)计算方法土压力水压力总压力墙左侧P p kN墙右侧P a kN墙左侧E W1/kN墙右侧E W2/kN(Pp+EW1)kN(Pa+EW2)kN方法179.44460.270122.304422.576201.74882.850方法248.63422.494153.360579.345201.99989.810方法3————444604.333方法4(1)54.29400.470153.360579.345200.61971.820方法4(2)156.33353.540153.360579.345304.99919.550与方法4(2)比较(%)方法1-49.230.2-20.3-27.1-33.9-4.0方法2-68.919.500-33.87.6方法3————45.6-34.3方法4(1)-68.315.100-34.25.7方法4(1)的主动滑裂面角为35°,被动滑裂面角为60°,方法4(2)主动滑裂面角为37°,被动滑裂面角为60°;且方法4(2)的主动和被动总土压力与板桩墙的外法线夹角均为20°,总水土压力与板桩墙的外法线方向夹角分别为7.6°、10.1°.从表1可以看出,在考虑板桩墙与土摩擦力的情况下,方法3的主、被动总压力偏差都较大,达13以上,说明此时这种用朗肯理论水土合算已不适用.方法1和方法2的被动总土压力偏差较大,接近13,主动总压力偏差较小,较为符合.方法4(1)与方法4(2)比较,可以发现,考虑板桩墙与土间摩擦力,被动总压力偏差较大,主动总压力偏差较小.值得说明的是当板桩墙与土间摩擦角δ较大时,用库仑土压力理论计算被动土压力(即方法4(2))也有较大的误差.3　轻型井点降水与明沟降水的二维渗流计算当地下水位比较高时,常常需要采取降水措施.当降水在6m 以内,土体渗透系数较小,可以采用一级轻型井点降水[8]或者明沟降水.下面通过对两者的水土压力对比计算,分析不同方向渗流下水土压力的计算特点.图6　轻型井点降水流网等势图算例4,某一工程进行轻型井点降水,基坑深5m ,板桩墙深7m ,有锚杆,坑外坑内水位与地面平,不透水层距原地面13m ,井点设计采用干式真空泵轻型井点,井点管埋深7m ,距基坑板桩墙1m .基坑地层为正常固结粘土为主,饱和容重γsat 为18.2kN m 3,摩擦角φ为30°,C 值近似取为0,土与板桩墙的摩擦角δ=0.6φ.设基坑内外水位不变,参考文献[8～10],通过对真空区和重力区分析,取降水管进水处水位为-6m ,设自由水面不变,用有限元分析,计算等势线如图6.计算时采用库仑土压力理论,对土楔体进行受力平衡分析.板桩墙左右侧墙面所受的水压力分别为E W 1、E W 2,总压力分别为 E a 、 E P ,计算同上例方法4(2).同上工程,采用明沟排水.在基坑内侧挖沟深0.3m ,宽0.4m ,每隔30～40m 处或基坑角处挖一集水井.设基坑上和坑底水位不变,通过有限元分析,计算方法同上例方法4(2).无渗流的情况水压力按静水计算,采用库伦土压力理论进行计算.计算结果如表2所示.表2　轻型井点与明沟降水水土压力计算(算例4)计算方法土压力水压力总压力P a kN P P kN E W 1/kN E W 2/kN ( E a ) kN ( E P ) kN 轻型井点降水136.94139.557.828.51144.40147.67明沟排水92.7647.02148.0422.10237.9968.38无渗流61.7494.39240.1019.60299.43113.19跟轻型井点明沟排水-32.3-66.31793.1159.764.8-53.7相比(%)无渗流　-54.9-32.42970.3130.3107.4-23.3从表2可以看出,由于渗流的作用,轻型井点降水能够降低板桩墙的总主动水土压力,提高板桩墙的总被动土压力,从而提高基坑支护结构的稳定性.从计算中可以看出,由于渗透力的作用,库伦土压力理论计算的土楔体,其主动土楔体的开裂角度为38°,大于45°-φ 2=30°,被动土楔体的开裂角度为72°,大于45°+φ 2=60°;明沟排水的主动开裂角为38°,被动土楔体开裂角为68°;无渗流静水时,主动开裂角为33°,被动开裂角为70°.对照表1也发现主动开裂角为37°(>30°),被动开裂角为60°.这说明渗透力对主动开裂角影响较大,对被动开裂角影响较小.其原因是被动区范围较小,基坑内外土的渗透力除了有竖直方向的分量,还有水平向基坑内方向的分量.4　二维超静孔压力分布的土压力计算[11]在某些砂石材料缺乏的地方,挡土结构后面的填土常常不得不采用当地的残积土作为填土.对于这种粘性填土,当其完全饱和或饱和度较高时,常常会遇到土中存在超静孔隙水压力的情况.这种超静孔压力一般不是线性分布,当墙面是排水结构时,它又是二维分布的,且它的大小和分布也是随时间变化的.具有二维分布的超静孔压的填土,朗肯土压力理论将不再适用.这时应该根据库仑理论用图解法来计算.此时二维超静孔压力分布可以通过分解成两个一维固结问题分别计算.图7　有顶盖时掺气浓度分布(Q=0.33m3s)算例5,墙面不排水的5m高挡土墙后填土是粘性土,上下边界都是排水层.饱和度S r=90%,容重γ=20kN m3,渗透系数k=1.0×10-6cm s,有效应力强度指标C′=5kPa,φ′=30°,初始孔隙比e0=0.5,压缩系数〗a=0.2MPa-1.分10层填筑,每层施工8h.填土完成后,在墙后填土表面一次施加均布荷载q=20kPa,再过24h后墙后填土中计算的孔压分布如图7(a)所示.它是二维问题,不能用朗肯土压力理论计算,只能用库伦理论的图解法进行计算,其滑裂面与水平方向夹角不是45°+φ′2=60°,而是53°,得到的滑裂面上的超静孔压,墙上的有效主动土压力p′a及水土压力之和p a 的分布如图8.计算得到主动状态下总压力为E a=155kN.(无孔压时总主动土压力E a=104.5kN).值得指出的是,u是滑裂面上的孔压,墙面上的孔压为0,所以u实际是通过渗透力引起的土压力作用在支挡结构物上的孔压.算例6　土质及荷载与边界条件与算例5完全相同,问题变成基坑开挖的支护问题。

高等土力学主要知识点整理（李广信版）第二章土的本构关系（一）概述材料的本构关系是反映其力学性能的数学表达式，一般为应力-应变时间-强度的关系，也称本构定律、本构方程。

土的强度是土受力变形的一个阶段，即微小应力增量小，发生无限大(或不可控制）应变增量，实际是本构关系一个组成部分，是土受力变形的最后阶段。

第一应力不变量kk z y x I σσσσ=++=1第二应力不变量kk yz xz xy z y z x y x I στττσσσσσσ=---++=2222第三应力不变量22232xyz xz y yz x yz xz xy z y x I τστστστττσσσ---+= 坐标系选择使剪应力为零3211σσσ++=I ,3231212σσσσσσ++=I 3213σσσ=I 球应力张量)(31)(3131321332211σσσσσσσσ++=++==kk m 偏应力张量ii kk ij ij s δσσ31-=，其中=≠=j i j i ii 10δ，克罗内克解第一偏应力不变量01≡=kk s J 第二偏应力不变量()()()[]23123222126121σσσσσσ-+-+-==ji ij s s J 第二偏应力不变量()()()213312321322227131σσσσσσσσσ------==ki jk ij s s s J 1.土的应力应变特性：非线性（应变/加工硬化、应变/加工软化）、剪胀性、弹塑性、各向异性、结构性、流变性（蠕变、应力松弛）。

加工硬化：应力随应变增加而增加，但增加速率越来越慢，最后趋于稳定（正常固结黏土、松砂）加工软化：应力一开始随应变增加而增加，超过一个峰值后，应力随应变增加而减小，最后趋于稳定（超固结黏土、松砂）剪胀性：剪应力引起的体积变化，含剪胀和剪缩土的结构性：由土颗粒空间排列集合、土中各相和颗粒间作用力造成，可明显提高土的强度和刚度。

灵敏度：原状黏性土与重塑土的无侧限抗压强度之比土的蠕变：应力状态不变条件下，应变随时间逐渐增长的现象，随土的塑性、活动性、含水量增加而加剧土的应力松弛:维持应变不变，材料内应力随时间逐渐减小的现象压硬性：土的变形模量（指无侧限，压缩模指完全侧限）随围压而提高的现象。

再议水压率李广信(清华大学水利水电工程系)11翻开贵刊2007年11期,看到方玉树先生的反驳文章,题目赫然是/水压率理论被推翻了吗0[1],对照拙文的题目/有效应力原理能够被推翻吗0[2],可见方先生以水压力理论取有效应力原理而代之的决心,尽管他目前还谦虚地借用有效应力原理的外衣。

21所谓孔隙水压力,在5岩土工程基本术语标准6(BG/T50279-98)[3]中规定:/土中某点孔隙水承受的压力0。

既然是水受到的压力,用测管水头来表示就是天经地义的,即u=p w=C w h,不仅土力学如此,水力学、流体力学和普通物理学都是如此,没有打折扣(N)的余地。

所以/水压率原理0还是应大大方方地写成:R=R c+N C w h=R c+N u(1)这就与各种其他企图推翻有效应力原理的公式相似了。

在科学的探索过程中,在学术讨论和交流时,不能在大家公认的名词术语方面偷换概念,否则就没有讨论的共同语言了。

偷换概念的结果常常是名为创新,实为抄袭;名为发展,实为篡改。

31将饱和含水量小于流限的粘土试样放在压缩仪中,施加竖向压力p,在t=0的瞬时,会测得孔隙水压力u=p。

亦即如方先生认为不可能的/侧压力系数为10。

可见结合水是能够传递水压力的;随着时间延续,孔压消散,土体固结,密度增加,侧压力系数变为K0,可见结合水是可以流动的。

这个过程似乎没有/水压率0置啄的余地。

41本科5土力学6教材只是介绍了一些基础的知识,未必句句是真理。

何况我国目前多如过江之鲫的土力学教材,基本是互相/参考0和/借用0,所以方先生列出了5本教材以表明我的/结合水可以传递水压力0之谬,不会使分量增加了5倍。

如果方先生认为这些土力学教材都是真理,那么/水压率0之果早在20年前(5本教材中最早出版者的时间)就应结在/结合水不能传递水压力0之树上。

既然方先生已经用水压率理论对于5土力学6教材各章进行了大刀阔斧的全面改造,根本就是否定了它们的权威性,而我只是认为/结合水能够传递水压力0,就似乎违背了公理而大逆不道。

土力学笔记（清华二版）第一章土的物理性质和工程分类1.1土的形成1.土的特点碎散性岩石风化或破碎的产物，非连续体受力易变形，强度低体积变化主要是孔隙变化剪切变形主要由颗粒相对位移引起三相性固、液、气受力后由三相共同承担相间存在复杂相互作用孔隙水和孔隙气可流动天然性自然界的产物，存在自然变异性非均匀性各向异性时空变异性2.三大问题、三大特性：强度、变形、渗透1.2土的三相组成1．根据三相比例不同，将土分为：饱和土、非饱和土、干土2．粒径级配：粒径大小及不同尺寸颗粒在土中占的百分比3．巨粒土＞60mm 、粗粒土＞0.075、细粒土≤0.075；4．粗粒土：以砾石和砂石为主的土，也称无黏性土5．细粒土：以粉粒和黏粒为主的土，也称黏性土6．粒径级配分析方法：筛分法（适用于粒径大于）0.075；水分法（粒径小于0.075）7．粒径级配曲线：横坐标为土颗粒直径（mm ），纵坐标为小于某粒径的土颗粒累积含量（百分比）8．粒径级配曲线的用途：了解土的粗细程度；粒径分布的均匀程度，分布连续性程度来判断土的级配优劣。

9．重要参数：1d50：平均粒径，表粗细；d10：有效粒径，细颗粒代表值；d30：连续粒径，表连续性；d60：控制粒径，粗颗粒代表值。

2两点：某粒径范围内土颗粒所占质量百分数，陡多，缓少，平缺10．不均匀系数Cu ：Cu=d60/d10；Cu ＞5不均匀土11．曲率系数：1060230c d d d C ⨯=；Cc ＜1或＞3表示级配曲线不连续。

12．级配良好的土：土的级配不均匀Cu ≥5，且级配曲线连续（Cc=1-3），适用于填方工程。

13．土中水：自由水：毛细水，重力水；结合水：强结合水，弱结合水14．土中气：封闭气体，自由气体1.3土的物理状态15．最基本指标1土的密度：单位体积土的质量，g vργρ==;m2土粒比重（土粒相对密度）：土粒的质量与同体积纯蒸馏水在4℃时的质量之比，数值上等于土粒的密度。