强度折减法在某核电站泵房边坡分析中的应用

对边坡稳定性分别采用按极限平衡法、有限元强度折减法进行二维分析,通过对比计算结果来验证强度折减法分析的准确性与可行性。针对边坡开挖与加固及隧洞开挖的三维效应,采用三维弹塑性有限元强度折减法分析计算边坡的稳定性,计算分析了开挖完成、支护完成条件下边坡安全系数、位移分布、主应力分布、等效塑性应变区和潜在滑动面。结果表明,支护加固可较大提高边坡的稳定安全系数,三维强度折减法分析在岩土工程领域具有明显优势。     

黄土高边坡的三维动力稳定性分析——以固原九龙山边坡为例

针对高烈度区九龙山黄土高边坡的动力稳定性问题,调查分析了该边坡的工程地质条件和场地所在区域构造活动与分布特征。在考察动荷载循环、往复作用下黄土反应敏感性及动强度参数,以及沿坡高确定地震惯性力反应地震作用大小的基础上,将强度参数折减与有限差分方法结合,从而形成了边坡的拟静力强度折减有限差分分析方法。通过九龙山土边坡拟静力强度折减三维有限差分法计算分析,得到了强度折减条件下边坡的位移场和应力场,边坡关键点位移与折减系数之间的关系,以及黄土高边坡的动力稳定性安全系数。其分析结果与传统的二维极限平衡分析方法确定的结果基本一致,验证了边坡拟静力强度折减有限差分方法的合理性。     

基于强度折减法的三维边坡失稳判据

为了寻求基于强度折减法的三维边坡稳定分析的合理判据,依托珍珠坝边坡采用FALC3D对位移突变判据及塑性区贯通判据进行分析,并与ANSYS计算结果和二维极限平衡结果进行比较。结果表明:在采用强度折减法对三维边坡稳定分析时,位移突变判据根据总位移增量和折减系数关系曲线较易识别突变点且特征点选在边坡坡顶滑体上较为合理;针对塑性区贯通判据的不足,提出了以贯通率增量和折减系数曲线上的突变点来判别边坡稳定性。因此,三维边坡失稳判据的合理选用应遵循以下原则:建议在滑带明确的情况下采用塑性区贯通率增量突变判据,在滑带未知的情况下采用位移增量突变判据。     

基于FLAC3D的三维边坡稳定性强度折减法计算效率改进算法及其应用

强度折减数值分析方法在边坡工程中应用广泛,其中FLAC^3D已成为三维边坡稳定性分析中的主流计算工具。FLAC^3D自带强度折减计算模块,因其计算过程复杂且耗时,一些学者提出了基于"二分法"的强度折减法,但依然不能令人满意。本文利用Python编程语言对FLAC^3D进行二次开发,实现了"二分法"强度折减法的全面优化,包括初始应力场载入优化和计算时步上限优化,建议计算时步上限取初始应力场计算时步的1.50倍,此时既可以保证计算结果的准确性又可以节省计算时间。利用优化后方法对某公路复杂三维填方路基边坡进行了三维强度折减法计算,结果表明,在保证计算准确性的同时,本文方法计算效率较其他改进的强度折减法至少提升50%,较程序自带模块提升了300%以上。     

有限元强度折减法中边坡失稳位移突变判据的改进

采用有限元强度折减法分析边坡稳定性时,计算所得安全系数依赖于所采用的失稳评判标准。目前边坡失稳判据的多样性及其不确定性,使得计算所得安全系数的合理性及其唯一性受到了质疑,因此如何选择或改进一种既能满足工程需要又易于操作的失稳判定准则显得尤为重要。针对某一典型边坡算例分别按计算收敛性、特征部位位移突变性和塑性区贯通性等3 个失稳判据考查判据的合理性及其适用性,对比研究三维与二维数值所得边坡稳定安全系数的差异,然后对位移突变失稳判据进行了一种改进,通过绘制位移突变倍率(当前分析步位移与上一分析步位移之比)跟折减系数的关系曲线来确定边坡安全系数,并通过简化的实例模型证明了其适用性。     

基于Hoek-Brown准则的三维边坡变形稳定性分析

采用Hoek-Brown准则,分析三维边坡在开挖扰动下的变形稳定性。以某露天矿边坡为工程背景,利用快速拉格朗日差分法(FLAC3D)建立三维数值分析模型,并在边坡中布置若干监测点,利用FISH语言编制相应位移插值程序,探讨边坡开挖引起的动静态位移响应,从宏观角度揭示出边坡开挖后,不同区域的变形,为工程实践提供指导;介绍了强度折减技术在Hoek-Brown准则中的实施方法,采用计算不收敛失稳判据,计算边坡安全系数,从而进一步推广Hoek-Brown强度折减法在三维边坡稳定性分析中的应用。     

基于FLAC3D强度折减法滑坡三维稳定性研究——以三峡库区白果树古滑坡群为例

利用FLAC3D的内置fish语言基于强度折剪法对其进行了二次开发,编制了滑坡体安全系数计算程序。结合工程实例,将这种方法应用于已知滑动面条件下的滑坡稳定性分析中,得到滑坡体天然状况下的二维、三维安全系数,同时对内摩擦角和粘聚力参数进行了敏感性分析,对比并分析了传统极限平衡法和FLAC3D强度折减法的计算结果,得出以下结论:①利用二维FLAC3D强度折减法计算的滑坡安全系数略高于传统极限平衡法;②利用三维FLAC3D强度折减法所得滑坡安全系数略高于二维FLAC3D的计算结果;③滑坡安全系数都是随着内摩擦角和粘聚力的增大而增大,其中内摩擦角和粘聚力的变化对三维FLAC3D强度折减法比二维求解安全系数影响更大一些。然后模拟了滑体在天然条件下处于极限平衡状态的滑动方向,并得出天然状况下滑坡整体基本稳定,仅滑坡后缘可能出现局部失稳的分析结果。同时探讨了三维FLAC3D模型建立过程中滑坡外围不动体及滑带以下基岩所取范围大小、网格剖分粗细程度对计算安全系数的影响。     

基于边坡稳定性计算方法比较分析

将强度折减法应用于边坡稳定性分析中,折减土体强度,代入有限元程序进行计算,直至计算不收敛,此时的折减系数即为安全系数。结合工程实例,将强度折减法应用于边坡稳定性的分析,利用瑞典圆弧法,结合岩土工程设计类软件天汉以及有限元分析软件ABAQUS分析边坡稳定性,并对安全系数进行对比,3种计算方法得出的安全系数差别不大,安全系数精度都能满足工程要求。     

强度折减有限元法中的单元阶次影响分析

强度折减有限元法是当前较为有效的边坡稳定性评价方法,且应用越来越广泛。但影响强度折减有限元法的因素有很多,单元阶次是其中比较重要的一个。通过3个经典算例,这些算例分别是二维地基承载力问题、二维边坡和三维边坡问题,分析了单元阶次的选择对强度折减法的影响。计算结果表明：随着单元的增多,线性单元和二次单元都从大于真实解的一侧来逼近真解;相对于二次单元,由于线性单元过“刚”,因此,会过高地估计安全系数,对于实际工程会偏于危险,且误差大,二次单元的误差是线性单元误差的1/8左右。在采用系统最大位移收敛与否的评判标准的基础上,利用二次单元来进行强度折减分析,则可以弥补这种线性单元的不足,得到更加合理的安全系数。二次单元比线性单元更适合于强度折减有限元法。     

基于强度折减和容重增加法的三维边坡稳定性分析

利用以变形变化趋势为边坡破坏判别标准的强度折减法和容重增加法,对某水电站高边坡进行了三维稳定分析,得到了三维最危险潜在滑动块体及相应的安全系数。同时与广义楔形体法二维计算结果进行了比较。计算结果表明,以边坡变形变化趋势为破坏判别标准,基于强度折减法和容重增加法的三维边坡稳定分析是合理可行的。      

基于软化特性的三维边坡强度折减有限元分析

基于现有边坡强度折减有限元的基本原理,建立了一套新的强度折减算法。该算法把强度参数降低的过程看成脆塑性应力跌落的过程,并设置一折减增量 ,当 充分小时认为系统达到收敛条件而退出循环,最终的折减系数就是边坡的安全系数。该算法很容易扩展到三维边坡的稳定性分析,得到边坡滑动面滑动和破坏的发展趋势。通过2个算例,分析了进行三维模型计算的重要性,建议当地质条件复杂时宜从三维角度分析计算边坡,不能简单地将问题当作平面应变处理。     

基于聚类分析的边坡稳定性研究

在边坡稳定性分析中,滑面搜索方法对于得到边坡的稳定性系数至关重要。现今两种常用的边坡稳定性分析方法,极限平衡法以及强度折减法,可以得到极限应力状态下边坡的潜在滑面。而为了评价边坡在当前应力场下的稳定性,提出一种基于聚类分析概念的滑面搜索方法。这种方法综合考虑了边坡各点的位移以及点安全系数,采用K均值聚类法将具有相近性质的点进行分类,把边坡体分为潜在危险区以及稳定区。潜在危险区的边界即为边坡在当前应力状态下的潜在滑面,并运用矢量和方法计算得到其稳定性系数。采用上述方法对3个算例进行验算,并与极限平衡法以及强度折减法进行对比。结果表明：（1）基于聚类分析的边坡稳定性评价方法具有较好的合理性和可靠性;（2）当边坡处于临界状态附近,通过聚类分析得到的潜在滑面与其他两种方法得到的滑面基本一致;当边坡处于较稳定状态时,当前应力状态下的潜在滑面要比强度折减法得到的滑动带范围要更深更广,但求出的稳定性系数也较高;当存在软弱夹层时,矢量和计算出的滑面更为陡峻,稳定性系数比其他两种方法的结果要小。     

一种双折减法与经典强度折减法的关系

经典强度折减法的出发点是对强度值进行折减,进而导致对强度（准则）中的参数使用同一个折减系数。使边坡处于整体临界平衡状态的折减系数恰是边坡的强度储备安全系数。一些学者陆续提出了对黏聚力和内摩擦角采用不同折减系数的所谓双折减法,但尚无严格的理论依据。提出双折减法的出发点是对强度（准则）中的参数进行折减,或简称“强度参数折减”,这样可自然导致采用不同的折减系数,但此时已无法自动得到基于强度储备的安全系数。为此,提出了定义安全系数的新框架--基于参照边坡的安全系数定义,为双折减法建立了理论基础。将新提出的一种双折减法与经典强度折减法进行了比较,发现可以将经典强度折减法纳入该双折减法的计算过程,并从理论上证明了该双折减法的安全系数几乎总是小于经典强度折减法的安全系数,数值模拟实例也证实了这个结论。这表明,经典强度折减法有可能高估了边坡的安全性,建议的双折减法可作为边坡稳定分析的备选方法之一。     

考虑地应力分布形式的坝址区初始应力场二次反演方法

在建的乌东德水电站区域构造运动强烈,断层、褶皱及节理较为发育,加之河谷纵横、沟谷深切,导致该区域岩体地应力分布极为复杂,非线性特征显著。为准确获取坝址区初始地应力场的分布特征,通过对地应力分布形式进行优化,提出了考虑地应力张量分布特征的岩体应力场二次反演方法。首先,将地应力张量分解成自重应力、构造应力及非线性应力3个部分,同时依据坝址区地形地貌、地质构造及非均匀性岩层等因素,构建FLAC3D计算模型,利用叠加法原理对坝址区初始地应力场进行一次反演。然后,主要考虑坝址区附近小规模构造及局部开挖条件的影响,构建三维离散元（3DEC）精细模型,同时从一次反演中提取精细模型的初始地应力侧压系数和自重修正系数,并利用均匀设计法进行优化,对精细模型进行二次反演。研究结果表明：一次反演计算获得的初始应力场整体上与实际情况吻合较好,但在局部构造附近与实测值相差较大;二次反演考虑局部地质构造及开挖条件的影响,同时结合一次反演的计算成果,各测点的应力计算值与实测值基本一致;通过将大型商业软件FLAC3D和3DEC相结合来反演复杂条件下工程区初始应力场,在一定程度上实现了岩体应力场在大范围内整体连续、小范围内受局部地质构造等影响导致不连续分布的特点。该研究方法可为类似工程提供借鉴。     

修正有限元强度折减法与失稳判据在边坡稳定分析中的应用

现有的有限元强度折减法用于边坡稳定性分析,在土体变形特性的认识和失稳判据上有缺陷。基于土体的应力-应变曲线和邓肯-张非线性弹性模型进行有限元强度折减法计算分析,在现有方法仅对强度参数--黏聚力c和内摩擦角φ进行折减的基础上,提出了切线弹性模量Et的折减算法,更全面地体现折减的内涵;提出了以潜在滑动区域的能量积分变化作为判据,根据能量积分与折减系数关系曲线的特征,将其拟合成两段直线,以其交点所对应的折减系数作为边坡的安全系数。比较分析了修正Duncan有限元强度折减法、失稳判据与文献的计算结果,同时分析了边坡滑动区能量积分计算范围对结果的影响,并用于南水北调东线实际工程计算分析。结果表明：修正Duncan有限元强度折减法能较全面考虑土体的变形和强度特性,新的边坡滑动失稳判据可行且具有较好的计算稳定性,能保证工程建设的安全。     

双强度折减法的量化分析

双强度折减法（DRM）采用黏聚力和内摩擦角的非等比例的折减策略,更加符合真实情况,但对于该方法本身的研究仍然缺乏定量的分析与合理的阐述;通过安全系数的定义、折减路径以及临界状态这几个概念,在强度参数坐标系中,对双折减法和传统强度折减法（SRM）进行量化分析与对比,两种方法具有统一性,表现为DRM相当于从虚拟初始强度点开始折减的SRM.边坡失稳的渐进过程本质上为强度衰减的过程,其中c、tan φ遵循正比例函数关系;同时,强度储备的意义反映为初始强度面积S向临界强度面积S'逐步地逼近的过程,从而得到安全系数表达式;强度折减法的研究应围绕强度来进行,对于强度参数的折减只是强度弱化的实现形式,通过安全系数定义的分析以及算例结果数据的回归分析,表明临界状态是边坡的固有属性,一旦边坡的几何尺寸给定,它将是一个与初始强度参数无关的量.SRM和DRM的结果的异同分析对比,实际上是对临界点的分析.     

地震作用前后极限平衡法和强度折减法的边坡稳定分析对比

分别运用极限平衡法和ANSYS强度折减法对国内某边坡地震作用前后稳定性进行了分析,考虑地震作用时采用了拟静力法施加水平和垂直两个方向的地震力.根据地震作用前后两种方法得到的边坡稳定系数、滑动面、滑动体的对比分析发现,两种方法得出的安全系数、滑动面相差不大,说明了强度折减法运用到边坡稳定分析中是合理的.强度折减法得出安全系数均比极限平衡法略大,这是由于强度折减法考虑了土体的应力-应变关系,使得计算结果更符合实际.通过地震作用前后边坡的水平位移图的对比,总结了一些规律,为地震区边坡的锚固提出相应的建议.