边坡渐进破坏的动态强度折减法研究

边破的变形破坏是量变积累到质变的渐进发生过程,由坡体内部潜在滑动面的逐渐破损并扩展至整体滑面。基于强度折减法思想,提出模拟边坡渐进破坏的动态强度折减法。该方法利用屈服接近度指标YAI确定边坡的破损区域,将YAI<0.2的区域定为破损区,通过不断动态折减局部破损坡体的强度参数,使边坡潜在滑动面渐进演化至贯通,边坡达到极限平衡状态。算例计算表明,该方法可有效解决强度折减法中折减范围的问题,克服整体强度折减法破损区过大的缺陷,获得的位移变化趋势与破损区演化过程也具有良好的一致性。动态强度折减法真实再现了边坡渐进失稳过程,为强度折减法更有效地应用于边坡稳定性定量评价提供有效途径。     

物质点强度折减法及其在边坡中的应用

物质点法适用于模拟连续介质大变形,如边坡失稳全过程。在物质点法中应用强度折减法,用于边坡稳定性评价。与极限平衡法相比,二者安全系数计算值、滑动面位置结果基本一致;与有限元强度折减法相比,物质点法失稳评价标准的物理意义明确。利用物质点法大变形计算优势,评价边坡失稳后的破坏后果,通过算例说明其评价不同安全系数下的滑坡堆积形态及滑移距离的能力,尤其是评价滑坡对临近建筑物的影响程度的能力。物质点强度折减法可用于边坡稳定性评价及边坡破坏后果评价。     

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

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

LU Kun-lin,ZHU Da-yong,YANG Yang 均质边坡准三维安全系数实用计算曲线

针对三维极限平衡法计算边坡安全系数存在着分析过程繁琐、较难普及推广等诸多不足,建立了一套实用的三维边坡安全系数计算曲线。假设三维滑面为中间圆柱和两端为部分球面的组合滑面,采用基于滑面正应力修正的极限平衡分析方法,计算36 864个均质边坡,根据计算结果绘制该套曲线。依据该套曲线可以得到,（1）根据已知边坡体的相关指标快速获得三维安全系数,能够得到较二维分析方法更准确、更客观的稳定性评价;（2）根据实际坡体的稳定状况反分析获得滑面土体强度参数,有效地克服二维反分析得到滑面土体强度参数偏高的弊端。研究成果可为均质边坡三维稳定性初步评价提供简单、实用的参考依据。     

应变软化边坡渐进破坏及其稳定性初步研究

采用应变软化本构模型和矢量和法,提出了一种模拟应变软化边坡渐进性破坏过程的分析方法。该方法可得到应变软化边坡的坡体材料强度参数、滑面状态、稳定安全系数、边坡破坏状态的变化过程,据此可分析边坡的渐进破坏。采用所提方法对应变软化材料边坡算例进行了分析,并与3种极限平衡法进行了对比,结果表明:(1)滑面矢量和安全系数介于传统极限平衡法峰值安全系数与残余安全系数之间,滑面强度参数部分处于峰值状态、部分处于残余状态,弥补了传统极限平衡法不能得到真实的应变软化边坡安全系数的缺陷。(2)强度参数弱化的区域与边坡破坏区域位置一致,坡脚处滑面最先出现软化,随着塑性剪切破坏的发展,其由软化状态变为残余状态,邻近部分滑面开始出现新的软化,接着进入残余状态,边坡破坏由坡脚向坡顶发展,随着这个过程的持续进行,边坡的渐进性破坏不断发展,边坡安全系数不断降低,直至渐进性破坏过程结束。(3)随着残余软化参数的增加,滑面位置逐渐变浅,矢量和法安全系数逐渐增加。该方法能够较好地分析边坡的渐进性破坏过程,具有较好的应用前景。     

基坑土钉支护边坡有限元稳定性分析方法探讨

针对基坑土钉支护开挖边坡的稳定性,讨论了极限平衡方法的局限性,指出应用有限元方法分析评价基坑土体结构稳定性是近年来研究的新趋势。采用滑面应力分析法和强度折减法两种有限元方法对基坑边坡的稳定性进行了研究。对于无支护的基坑边坡,两种有限元方法得到的安全系数和滑动面结果一致,且滑动面形状随基坑土体强度的降低保持不变。对于采用土钉支护的基坑边坡,两种有限元方法得到的安全系数结果一致,不同阶段下滑动面形状和位置不同,这主要是由于在强度折减时,仅仅对土体强度参数进行折减,而未考虑土钉自身强度的降低。最后指出,在应用强度折减法时,是否考虑土钉自身强度的折减有待进一步研究。     

基于动力学和材料软化特性的边坡渐进破坏特征研究

边坡的失稳是一个从量变到质变的动态渐进破坏过程,此问题也是边坡领域研究的重点与难点之一。在考虑岩土材料软化特性和动力学求解的基础上,建立了边坡渐进破坏仿真的理论框架;利用ABAQUS软件的动力显式求解模块实现了边坡的渐进破坏仿真;根据塑性应变揭露了剪切带的扩展过程,由软化本构确定了滑面材料的分区演化规律,根据等效塑性应变确定了边坡的滑面,通过滑面位置将边坡分为滑体、滑带、滑床,并分别研究了边坡各分区内部特征点运动学变量的发展过程,从而揭示了边坡的渐进破坏过程;基于材料参数沿滑面的时空分布,利用矢量和法得到了边坡不同演化阶段的安全系数。对比该方法与Bishop法确定的滑面位置与安全系数,发现两种方法峰值和残余强度对应的安全系数比较接近,该方法搜索所得滑面位于Bishop法自动搜索的滑面之间,验证了此方法的合理性及可靠性。最后分析了材料软化特征对边坡稳定性的影响,在保持其他参数不变的条件下,增大残余黏聚力,边坡的滑面位置加深,安全系数的初始值减小,安全系数的快速减小阶段有所推迟,并且快速减小阶段经历的时间有所延长,稳定后的安全系数有所增大。保持其他参数不变,增大残余黏聚力对应的等效塑性应变阈值,边坡的滑面位置加深,安全系数的初始值减小,安全系数的快速减小阶段有所推迟,但快速减小阶段经历的时间基本不变,达到稳定的时间有所推迟,同时稳定后的安全系数略微有所增大。     

基于Hoek-Brown准则的非等比强度折减方法

针对节理岩质边坡稳定性分析问题,开展了基于Hoek-Brown破坏准则的强度折减法研究。边坡的失稳过程是一个渐进累积破坏过程,破坏过程中各强度参数的衰减程度不同,在强度折减法中对应的折减系数也不同,不同折减系数间关系的确定及其综合安全系数的定义值得深入探讨。首先从岩体材料软化（硬化）特性出发,根据岩体强度参数从峰值强度到残余强度的变化规律,推导了Hoek-Brown破坏准则中3个强度参数折减系数间的数学关系式,然后以折减前后滑面上抗滑力之比定义综合安全系数来评价边坡的稳定性,最后结合算例验证了该方法的合理性,可有效应用于节理岩体边坡的稳定性分析。     

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

强度折减法在评价边坡稳定性时不需要假设和搜索临界滑裂面,相对于其它方法具有一定优势,但失稳判据尚无明 确标准。滑移线场理论（SLFT） 可以计算得到的极限状态下的边坡坡面曲线（简称极限坡面曲线）,已有研究表明极限坡 面曲线与边坡坡面的相对位置关系可以判断边坡稳定性。基于该结论提出一种边坡强度折减法失稳判据：不同折减系数计 算得到不同的强度参数,因此滑移线场理论计算得到的极限坡面曲线将发生变化,当极限坡面曲线与边坡坡面在坡底相离 时,判断边坡为稳定状态;当极限坡面曲线与边坡坡面相交于坡脚时,判断边坡为极限平衡状态;当极限坡面曲线与边坡 坡面相交时,判断边坡为失稳状态。对标准边坡考题的计算结果表明：提出的失稳判据收敛性较好,安全系数计算结果与 标准答案和已有失稳判据分析结论相差较小;传统失稳判据收敛指标的确定没有明确标准,当位移突变曲线光滑时很难找 到突变点,很难准确计算对应的安全系数,更重要的是传统失稳判据判断边坡极限状态受到人为主观因素的影响,而提出 的失稳判据可以实现失稳判据的客观标准化;确定的临界滑裂面形状以及在边坡中的位置都与标准答案基本一致,提出的 失稳判据适用于强度折减法。     

非稳定渗流条件下非饱和土边坡稳定分析

大量的滑坡灾害表明, 水的作用是影响边坡失稳的重要原因,合理地评价边坡的稳定性,需考虑边坡基质吸力场的变化与水的渗流作用。通过稳定与非稳定渗流条件下的非饱和土堤算例,探讨了Bishop和Fredlund非饱和土抗剪强度公式应用于强度折减有限元法的差异,并与极限平衡法的结果进行对比。结果表明,Bishop和Fredlund非饱和土抗剪强度公式的本质是相同的,而从非稳定渗流条件下的分析结果来看,强度折减有限元法分析的安全系数与极限平衡法分析临水坡和背水坡安全系数的最小值相当接近,强度折减有限元法更适合分析边坡的整体稳定性。此外,采用Bishop强度公式可以通过用饱和度代替基质吸力参数来体现基质吸力对强度的影响,在分析饱和度较高的非饱和土边坡稳定性时较为方便。     

禁忌搜索在边坡稳定分析中的应用

在禁忌搜索和极限平衡分析的基础上,提出了全局确定边坡最小安全系数的禁忌搜索方法,并建立了响应的模型;该方法是一种启发式的全局搜索方法,可以搜索出全局最小安全系数,同时还能确定出与之对应的最危险滑动面。将该方法应用于不同几何形状和土层信息的边坡稳定分析中,结果表明,提出的方法能够避免局部最小解,并能搜索到全局最小安全系数。     

圆弧滑动边坡改进条分法的数值分析

结合传统极限平衡方法和现代极限平衡方法,对圆弧滑动边坡进行垂直条分,将条块间的作用力一般化,提出了基于圆弧滑动的改进条分法。建立了圆弧滑动边坡的力学计算式,确定了一般化的条间力影响系数和边坡安全系数的关系式。运用数值计算和迭代程序,得出滑体分别在力平衡和力矩平衡下的安全系数随条间力影响系数变化情况,从而定量地描述了边坡稳定状态。结果表明,圆弧滑动垂直条分改进方法反映了滑动面和垂直面相互影响情况,数值计算简便,而且精度较高。     

Seismic stability analysis of a layered rock slope

In general, the determination of the factor of safety and the location of the critical slip surface are two major challenges in seismic slope stability analysis. In this paper, a new approach for determining the factor of safety and the corresponding critical slip surface of a layered rock slope subjected to seismic excitations is presented, through a case study based on the combination of the shear strength reduction technique and distinct element method. According to this proposed method, the seismic factor of safety and the critical slip surface of the slope are estimated and compared with those obtained by the pseudo-static approach, combined with the limit equilibrium method. It is found that the factor of safety obtained from the proposed method is slightly greater than that computed by the pseudo-static analysis, with a difference of 4.2%, and that the critical slip surface obtained from the two methods is identical, which confirms the reasonability and feasibility of the proposed method.     

基于数值应力场的边坡临界滑动场

数值分析方法（有限元、有限差分等）可计算边坡应力、变形和破坏,但直接计算边坡安全系数比较困难;而极限平衡法需对边坡体内力作出假设,计算精度受到限制。利用数值分析方法计算得到的应力场,重新定义推力最大原理中推力的含义,将常规边坡临界滑动场数值模拟方法作进一步的改进,提出基于数值应力场的边坡临界滑动场方法。该方法可充分发挥两大方法的优势,能够准确、方便地搜索出边坡临界滑动面位置,并给出更合理的边坡安全系数。通过两个算例的稳定性分析,并与其他方法进行比较,验证了该方法的准确性和合理性。     

基于交替变量局部梯度法对土质边坡的三维稳定性分析

临界滑动面的获得是土质边坡稳定性分析中一个必要的过程,然而现有分析方法难以保证该临界滑面的准确性,而且部分分析方法在计算过程中会陷入局部最小值陷阱。基于单变量方法和最速梯度法,提出了交替变量局部梯度法,该方法能够突破最优化方法应用在边坡稳定性分析上的瓶颈。首先,通过斯宾塞极限平衡方法和网格搜索法（或者其他确定临界滑面的方法）,获得边坡初始椭球状临界滑面;其次,在该滑面上布置若干节点作为变量,目标函数为安全系数Fs关于空间节点坐标Zi的方程,为了使目标函数快速降低,沿负梯度方向循环优化每个节点,当满足一定精度要求时,即可获得三维土质边坡非圆弧状临界滑动面;最后,通过算例证明该方法可行,计算结果可靠。     

基于力平衡的三维边坡安全系数显式解及工程应用

将基于滑面正应力分布的二维极限平衡显式解法推广到三维边坡稳定性分析,提出适合一般形状空间滑面的三维边坡稳定性简便计算方法。首先,假设三维滑面正应力的初始分布,然后对其修正,使滑体满足所有力的平衡条件,导出关于安全系数的2次代数方程,最终得到三维边坡安全系数显式解。该方法克服了目前条柱法存在的一些难以解决的困难,不需要考虑条柱间作用力具体分布,满足主要平衡条件,计算过程简单。算例表明,该显式解与楔体三维极限平衡解析解完全一致,与对称三维边坡严格极限平衡解答接近。该方法已用于某重力坝坝肩三维滑体初步稳定性分析,计算结果有助于坝肩边坡的加固设计。     

岩质块体渐进破坏的稳定分析方法

岩桥的破裂贯通是一种非连续变形现象。强度折减技术虽然是渐进破坏过程模拟的主流技术,但其不具备描述非连续变形现象的能力。提出连通率折减法和刚度折减法来模拟含非贯通结构面岩块的渐进破坏过程,并建立考虑渐进破坏过程的块体稳定分析方法。首先,引入Goodman单元来描述共面非贯通结构面的岩桥部分和裂隙面部分,建立静力平衡方程来求解滑裂面内岩桥单元和裂隙面单元的应力。其次,岩桥单元采用Griffith准则来判别破坏,利用连通率折减法来描述其破裂;裂隙面单元采用MC准则来判别破坏,用刚度折减法来描述其屈服;通过循环迭代,模拟岩桥单元的破裂过程和裂隙面单元的应力调整过程,实现整个滑裂面的渐进破坏过程模拟。然后,定义考虑渐进破坏过程的滑裂面极限状态;通过自重超载方式将滑裂面推送至极限状态;基于极限状态设计的理论框架,计算度量块体稳定性的安全系数指标。工程实例分析结果表明,渐进破坏过程模拟结果与现场调查结果是一致的。渐进破坏过程的模拟实现,有助于岩质边坡变形破坏机制分析从定性分析阶段扩展至定量分析阶段。