基于Mein-Larson入渗模型的浅层降雨滑坡稳定性研究

降雨入渗模型与边坡稳定性方法有机结合是研究降雨诱发滑坡稳定性的有效方法。在深入分析降雨入渗基本理论的基础上,将Mein-Larson入渗模型与基于饱和土、非饱和土以及近似非饱和土的无限边坡稳定性分析方法有机结合,综合考虑了高强度短历时和低强度长历时两种降雨情形,拓展并改进了基于Green-Ampt入渗模型的饱和土和近似非饱和土无限边坡稳定性模型,并依托典型边坡实例对各评价模型的计算结果和适用性等进行了对比分析。     

考虑降雨入渗的非饱和边坡稳定分析

分析了降雨入渗的边坡稳定问题,通过正交实验确定饱和-非饱和渗流计算参数。采用SEEP/W计算饱和-非饱和渗流,利用SLOPE/W模块计算树坪滑坡在降雨入渗情况下的安全系数。     

降雨诱发型残坡积土滑坡形成机理分析

以马家村滑坡为例,用有限元方法分析了降雨入渗条件下残坡积土斜坡的暂态非饱和渗流场,用极限平衡方法分析了残坡积土斜坡的渗流稳定性。分析了斜坡结构、临空面、降雨对斜坡变形与稳定性的影响,分析了降雨诱发残坡积土滑坡的过程。降雨诱发残坡积土滑坡的机理为:降雨入渗使残坡积土斜坡中形成软弱带;坡脚开挖使斜坡坡面出现拉张裂隙;降雨使坡面裂隙中充满水,导致斜坡土体饱和,抗剪强度降低,同时产生动水压力,不利于斜坡的稳定;当斜坡的下滑力大于抗滑力时发生滑动。     

降雨作用下考虑膨胀推力的膨胀土边坡稳定性分析

为研究膨胀力作用下膨胀土边坡的稳定性,通过合理假设,在传统剩余滑坡推力法中引入膨胀力项,得出了适用于膨胀土边坡的滑坡稳定性系数计算公式。基于该公式对某膨胀土坡的离心模型降雨试验进行了模拟和计算,通过计算结果和试验结果的对比验证了该法的可行性;同时分析了膨胀力施加与否、不同裂隙分布位置和深度、不同降雨历时和降雨强度下膨胀土边坡的稳定性,得出了安全系数的变化规律。结果表明,降雨条件下考虑膨胀力时膨胀土边坡的安全系数相对于未考虑膨胀力时大幅下降,符合膨胀土边坡降雨破坏的实际情况;裂隙位于坡中时膨胀土边坡较裂隙位于坡顶时更危险;降雨强度和降雨历时也对膨胀土边坡的稳定性造成了影响,但超过一定程度时这种影响会受到限制。该方法为工程中膨胀土边坡稳定性的计算提供了参考。     

渗流作用下节理型黄土开挖边坡塌滑破坏分析

黄土节理在黄土地区普遍发育.通过野外现场地质调查。发现调查区域黄土开挖边坡垂直节理较多,不稳定的黄土边坡体被黄土垂直节理分割后,与其下部缓倾角黄土节理组合而成的裂隙滑动面在降雨作用下往往能形成塌滑破坏。建立了调查区域露天煤矿某黄土开挖边坡的概化三维数值分析模型,模拟分析了在3种降雨条件(20 mm/d,40 mm/d,100 mm/d)下降雨2 d后概化模型的孔隙水压力、位移、塑性区面积情况,提出此类黄土开挖边坡应注意黄土垂直节理与缓倾角黄土节理组合而成的裂隙滑动面的防治,为工程地质条件地区相同的黄土开挖边坡塌滑变形破坏分析以及防治提供了一种可行的分析方法。     

降雨入渗条件下内河航道岸坡稳定性分析

降雨入渗是引起内河航道岸坡失稳的重要因素。目前降雨入渗条件下岸坡稳定性数值分析中,关于岸坡竖向剖面的饱和度与孔压分布规律,主要针对不同的雨型分布开展研究,而不同的降雨强度、降雨历时和土体饱和渗透系数也会显著影响饱和度与孔压随深度的分布,从而影响岸坡的稳定性;降雨作用下的许多岸坡失稳现象发生在降雨停止后的某一时刻,目前鲜见针对降雨作用停止后岸坡稳定性变化规律的数值模拟研究;当降雨强度与岸坡土体饱和渗透系数较为接近时,对岸坡的稳定性最为不利,已有研究中没有结合降雨强度考虑土体饱和渗透系数变化对岸坡稳定性的影响。基于此,通过建立ABAQUS有限元数值模型,分析降雨作用对岸坡饱和度、孔隙水压力、位移变形和安全系数的影响。结果表明:降雨入渗会使岸坡浅层土体的饱和度与孔压急剧升高;雨水渗流及基质吸力会使岸坡土体产生变形,强降雨会使浅层土体发生严重塑性破坏;降雨强度一定时,随着降雨时长增加,岸坡稳定性降低,土体渗透系数较小时,岸坡安全系数最小值出现在降雨结束后某一时刻;土体渗透系数对于岸坡稳定性的影响需同降雨强度结合起来考虑。     

间歇型降雨作用下缓倾滑坡稳定性数值模拟

为研究间歇型降雨对缓倾滑坡稳定性的影响,以洛带樱桃沟滑坡为例,通过设置不同的间歇工况对洛带樱桃沟滑坡进行数值模拟,结果表明:非雨季坡体稳定性受间歇型降雨影响显著,且雨停后坡体稳定性恢复能力较弱;雨季同等雨量下,间歇周期越长,坡体越不稳定;长期间歇后的首次降雨对坡体强度弱化影响显著;降雨入渗易在缓倾坡体基岩面上储存,形成暂态地下水位和坡向渗流场,入渗水主要起参数弱化作用,饱和区的渗流推力对坡体安全系数影响较小;随着坡体水位的抬升,坡体安全系数对降雨的敏感性逐渐被削弱;降雨易在缓坡处聚集,形成局部高孔隙水压力区域,增强了坡体的储水能力;缓倾滑坡稳定性系数波动幅度小,在间歇型降雨作用下易发生蠕变滑动。     

不同含水率下考虑筋土界面刚度软化的加筋边坡地震响应∗

加筋土结构具有较好的整体变形协调能力以及抗震性能，被广泛应用于边坡工程中。基于不同含水率下残积土及筋土界面室内直剪试验数据，在有限差分软件中考虑不同含水率下土体和筋土界面抗剪强度参数，并通过 FISH 语言编程在数值计算中实现筋土界面剪切刚度动态变化过程，分析加筋边坡在不同含水率、筋材长度、筋材间距及地震峰值加速度工况下的位移、应力、加速度响应。模型计算结果表明：相比不考虑筋土界面刚度软化情况，考虑筋土界面剪切刚度软化情况下含水率对加筋边坡水平位移和加速度放大系数的影响较大；坡面侧向位移、 坡顶沉降和筋材拉应力与回填土含水率呈正相关，加速度放大系数与含水率呈负相关；在一定范围内增加筋材长度能够有效降低加筋边坡坡面的水平位移；筋材间距对含水率越低的加筋边坡坡面侧向位移影响越大；地震峰值加速度越大，含水率对坡面水平位移影响越大。     

三维土工网防护边坡整体稳定性分析

针对三维土工网护坡技术的前期阶段不断出现的边坡整体失稳现象,基于三维土工网护坡技术的作用机理,采用强度系数折减法,利用ABAQUS建立土工网防护边坡的三维有限元模型,以坡肩特征点的位移突变情况作为边坡失稳判断的标准,分别计算和比较了不同降雨条件下土工网防护边坡和不含土工网边坡的安全系数,探究了边坡失稳时坡体内部塑性区的发展情况和位移分布规律,并利用室内边坡抬升破坏试验分析比较了两种边坡的破坏抬升角度和最终破坏情况,与数值模拟分析结果进行比较。结果表明:土工三维网在一定程度上可以提高边坡的整体稳定性,但作用影响较小;降雨强度的大小对边坡稳定性的影响较为明显;在土工网和铆钉共同作用下,土工网和其覆盖的土体共同作用,呈现较好的整体性。     

碎石土边坡破坏机理的敏感性分析

在碎石土类边坡中常常发育稳定的地下水管网排泄系统,它们对控制地下水水位上升,保持边坡稳定十分重要。以官家滑坡为例,通过对滑坡稳定性系数有关的各因素敏感性分析,发现地下水是影响边坡变形破坏以及复发破坏的最主要因素。当坡脚开挖或坡体堆载时,会破坏管道状地下水排泄系统,降雨入渗导致地下水水位升高,从而引起坡体内孔隙压力比、水头高度和水力坡度增大,使潜在滑面上的孔隙水压升高,影响碎石土边坡的稳定性；同时,地下水位的升高降低了土体的内摩擦角,而因素敏感性分析发现,内摩擦角对边坡失稳具有极其重要的影响。在官家滑坡的后期治理中,治水作为主要工程措施的理念已经得到很好的贯彻,效果十分明显。     

非饱和土边坡三维稳定性分析方法及考例

为了同时考虑边坡的饱和状态和三维尺度,根据非饱和土Barcelona模型进行了三维有限元数值模拟,确定了非饱和土边坡的应力位移分布。结合稳定性分析理论,将吸力对强度的贡献嵌入安全系数公式中,研究了非饱和土边坡的三维稳定性分析方法,并通过标准算例边坡的稳定性分析,确定了吸力对安全系数的贡献,并证明了该方法能够很好地退化到常规稳定性分析。     

面板砂砾石坝异常状态渗流计算及稳定分析

运用非饱和土渗流与强度理论,把渗流计算的水力条件用于边坡稳定分析之中,对沟后面板砂砾石坝进行了渗流场分析与边坡稳定计算,得到一些新的结论,可补充以前相关的研究结果。计算结果表明:坝料填筑分离造成的渗透各向异性尽管使浸润线升高,但是由于水平排水能力的提高反而使坝体的孔隙水压力总体上降低,从而也就提高了边坡的抗滑稳定安全性;当顶部面板与坝体脱开并导致异常状态渗流时,坝体顶部有一饱和区并且渗流出逸点很高,但坝体中部仍未饱和,异常状态渗流增加的孔隙水压力有限,对坝坡的抗滑稳定影响也很小,因此,按照最高浸润线的位置确定孔隙水压力进行抗滑稳定计算会导致错误的结果;坝顶异常状态渗流增加的孔隙水压力并不会导致下游坝坡滑动失稳,而坝体顶部渗透破坏会导致局部坍塌,最终发生溃坝事故。     

水对某库岸边坡稳定性的影响分析

针对某水电站库岸边坡，分析了地质条件，采用ADINA软件建立了有限元二维模型，用强度折减法对比分析了该边坡在自重，蓄水，降雨入渗及水位骤降的情况下的位移、应力、塑性区及安全系数的变化，判断分析了该边坡在以上各工况下的稳定性，对比分析了库水位变化及降雨情况下水的作用对边坡稳定性的不利影响。     

基于局部强度阶梯折减法的边坡渐进破坏研究∗

通过强度折减法获取边坡的潜在滑裂带是目前数值计算广泛使用的方法,研究边坡滑裂带形成的渐进过程及坡体内部土体的应力状态具有重要的理论意义和工程价值。基于边坡破坏的渐进性思想,针对均质边坡提出了一种局部强度阶梯折减的方法。该方法以屈服接近度(YAI)作为判别土体是否存在破损区的评价标准,将局部破损区按不同破损次序分配予不同的折减系数,利用通用软件ABAQUS对局部坡体进行强度折减计算,研究边坡渐进变形过程,直至边坡发生整体破坏。计算结果表明,坡体内部初始破损区参与了潜在滑裂面的形成过程,对边坡后期渐进变形具有贡献;有效解决了传统强度折减法潜在滑裂面上折减区和非折减区边界处的塑性应变值的跳跃性问题,合理地演化了边坡的渐进失稳进程。     

基于可滑性分析的土质边坡局部稳定性研究

传统的局部安全系数计算方法通常是基于某种屈服准则,通过判断坡体内部是否发生塑性屈服来评价边坡局部稳定性。实际上,发生塑性屈服的区域通常远大于滑动面的范围,仅以是否发生塑性屈服作为局部稳定性评价指标偏于保守,为此进行坡体可滑性分析,即考虑坡体内各点具有的形成滑动面的可能性,并定义该指标为滑动因子。同时,将坡体内各点的塑性屈服程度用剪切塑性屈服系数来表示。将两者结合起来共同评价坡体的局部稳定性。研究表明,边坡局部稳定性不仅与坡体内各点沿潜在滑动面的屈服程度有关,还取决于该点形成滑动面的可能性,可表达为剪切塑性屈服系数与滑动因子的乘积。用ACADS边坡考题验证了本文方法的正确性,为边坡稳定性分析提供了思路。     

降雨强度对泥石流起动影响的模型试验研究

利用自制小比例模型槽,采用可控降雨强度降雨模拟器,进行了人工降雨诱发泥石流的室内模型试验,研究降雨强度对泥石流起动的影响,并确定泥石流起动的临界降雨强度。基于临界降雨强度下泥石流的起动过程,对孔隙水压力和坡体内部含水率的变化规律进行了研究,并与现场试验进行了对比。运用Geodog软件对泥石流形成前坡体的位移场进行了分析,研究坡体破坏规律。研究结果表明,泥石流起动的临界降雨强度范围为43~50 mm/h。泥石流发生前存在一个孔隙水压力增加的过程,在泥石流形成后,孔隙水压力急剧下降。降雨强度影响雨水在坡体内部的渗透,是泥石流起动的主要原因。泥石流起动与坡体含水率和饱和度的变化有关,当坡体饱和度超过80%时,产生整体滑动破坏,并迅速转化为泥石流。     

全风化花岗岩滑坡稳定性与降雨关系分析∗

降雨是滑坡诱发的主要因素,当前我国防灾工作理念从注重灾后救助向注重灾前预防转变、从减少灾害损失向减轻灾害风险转变,因此明晰降雨诱发机理并建立合理的预警机制,对做好地质灾害防灾减灾工作至关重要。针对降雨型全风化花岗岩滑坡稳定性问题,以青岛崂山风景区全风化花岗岩返岭滑坡为实例,进行了不同含水率原状土剪切试验与大型物理相似模型试验。揭示了边坡的降雨响应规律,探究了全风化花岗岩滑坡稳定性与降雨关系,并拟合相关量化公式。同时采用ABAQUS 建立边坡流固耦合三维数值模型,基于强度折减法验证了公式的合理性。研究结果表明：（1）降雨型全风化花岗岩滑坡的破坏模式分为浸润侵蚀→表层变形→破坏加深→整体失稳4个阶段,变形期间坡体存在“ 鼓状凸起”与“片状溜滑”现象,最终发生推移式破坏;（2）边坡对降雨入渗的响应规律在水平及竖向空间上存在差异,降雨强度越大,含水率与孔压增速越大;（3）基于试验结果推导了滑坡含水率、安全系数与降雨强度、降雨持时之间的影响关系公式,数值模拟验证误差率较小,能够较好的描述全风化花岗岩滑坡稳定性与降雨之间的定量关系。研究结果为类似强~全风化花岗岩地区滑坡的预警与防治提供参考。     

下蜀土边坡地震稳定性的大型振动台试验研究(Ⅱ)——试验结果及分析

选用天然下蜀土,设计并完成了一几何相似常数为20、坡高为0.5 m、坡角为45°的模型边坡1 g大型振动台试验,通过对模型边坡动力特性、动力反应及宏观变形等三方面的研究,揭示了下蜀土边坡的地震反应机理。试验结果表明,下蜀土边坡有着稳定的一、二阶自振频率和阻尼比,随着地震动加载幅值(PGA)逐渐增大,其一阶自振频率逐渐减小,一阶阻尼比逐渐增大;坡高范围内,PGA放大系数均大于1,且自下而上逐渐增大,并在坡肩处达到最大,其在坡脚处的值小于水平方向其他位置处的值,输入El Centro波(经时间相似常数4.47压缩)时的值分别大于输入El Centro原波和Kobe波(经时间相似常数4.47压缩)时的值;随着输入PGA的增大,土体非线性特征逐渐增强,边坡也逐渐趋于失稳破坏状态;震后模型边坡坡脚、坡顶等部位出现了不同程度的裂缝,说明边坡中出现了通过坡脚的(潜在)滑动面。     

基于重力加载比例法的岩质边坡稳定性分析

根据岩质边坡的稳定性主要由边坡内部的结构面所控制的特点,采用有限元重力加载比例法对含有一组平行节理面和含两组节理面的岩质边坡进行数值模拟,分析表明:与强度折减法相比,采用有限元重力加载比例法能更加快速地计算出岩质边坡的安全系数,节省计算耗时和人工干预的工作量,并能通过塑性应变的发展分析其破坏过程,确定潜在的滑动面。     

下蜀土边坡地震稳定性的大型振动台试验研究(Ⅰ)——模型试验设计

为研究下蜀土边坡的地震稳定性,选用天然下蜀土,按照Meymand的相似法则,对一坡高为0.5m、宽度为1.05m、坡角为45°的模型边坡进行了几何相似常数为20的1g大型振动台模型试验详细设计。首先,运用数值模拟的方法确定了模型边坡尺寸适当的边界范围,模型边坡地基厚度、坡脚前缘和坡肩后缘长度均取为0.5m;其次,按照确定的模型几何尺寸设计了符合试验要求的模型箱,其内壁净尺寸(长×宽×高)为1500mm×1050mm×1100mm,属刚性模型箱,且不会与模型边坡发生共振;最后,制定了包括加速度和位移在内的数据量测策略,运用拟静力法求得下蜀土边坡平均屈服加速度为0.561g,综合确定了地震动峰值(PGA)按0.1g、0.3g、0.6g和0.9g的顺序逐级加载的原则,并制定了详细的试验步骤,从而可保障试验得到可靠的结果。