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

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

有限元强度折减法计算边坡稳定的对比分析

采用目前边坡稳定性分析比较流行的强度折减法,对比研究了MIDAS/GTS、FLAC、ANSYS配合Drucker-Prager（简称D-P）屈服准则和Mohr-Coulomb（简称M-C）屈服准则时软黏土、硬黏土、弱膨胀土3种工况下计算结果的偏差。软黏土工况下D-P准则和M-C准则计算结果的偏差相对较小,当边坡土体为硬黏土时,采用D-P准则与采用M-C准则计算结果的偏差明显增加。3种软件2种屈服准则下的计算结果都反映出,硬黏土的滑动面比弱膨胀土和软黏土的滑动面浅,而且同等情况下MIDAS计算得到的滑动面比ANSYS计算得到的滑动面浅;坡度较小时FLAC（M-C）计算的安全系数比MIDAS（M-C）计算得到的大,坡度较大时则相反;坡度较小时计算过程中先出现塑性区贯通,后出现计算不收敛;坡度较大时计算过程中先出现计算不收敛,后出现塑性区贯通。坡度较小时计算不收敛时的折减系数与出现塑性区贯通时的折减系数差别较大;坡度较大时这一差别较小,甚至计算到不收敛时塑性区仍未贯通,在用MIDAS计算时这一现象反映得更加明显。     

基于强度折减法的多级边坡稳定性分析

结合工程算例,采用强度折减有限元法,对边坡治理前后的稳定性进行了分析。结果表明,有限元强度折减法能更加直观的得到坡体的实际破坏形式,求得的边坡稳定系数更接近边坡的实际稳定状态,显示出其在边坡稳定性分析中的一定优势。     

一种考虑变形参数和强度参数协调折减的强度折减法研究

对于一种确定的岩体材料而言,岩体的变形参数与强度参数存在一种必然的内在联系。从材料性质来看,强度折减法的本质就是在当前边坡构型和岩体重度条件下,寻找另一种合适的岩体材料使边坡恰好达到临界状态。显然,若只对黏聚力和内摩擦角进行折减而其他参数保持不变,折减后的材料参数则可能违背岩体强度与变形参数之间存在的必然联系,即新找到的岩体材料客观上不存在。因此,在对黏聚力和内摩擦角进行折减时,有必要对岩体的其他力学参数进行调整。基于以上认识,提出了一种考虑变形与强度参数协调折减的强度折减法,探讨了变形参数（弹性模量和泊松比）、抗拉强度对安全系数、塑性区的影响,结果表明,同时考虑变形参数与强度参数协调折减的方案所得的安全系数最小,与极限平衡法的计算结果最为接近,且该方案所得的剪切塑性区主要集中在临空面附近,边坡模型中下部的剪切塑性区范围较小,同时,在边坡顶部区域出现张拉屈服区,此种塑性区分布范围最接近边坡实际破坏状态。     

非均质边坡强度折减法折减范围研究

强度折减法分析边坡稳定的研究多针对均质简单边坡,而当涉及到非均质边坡时,就存在选择局部区域还是选择所有区域进行强度折减的问题。本文基于FLAC3D的FISH语言自编整体强度折减法和局部强度折减法,在此基础上,对非均质边坡分别按照整体强度折减与局部强度折减进行分析,结果表明两者所得安全系数并不总是一致,整体强度折减法给出的安全系数更为合理,因此建议采用整体强度折减法分析非均质边坡的稳定性。     

基于强度折减法的土石混合体边坡长期稳定性研究

土石混合体边坡是自然界中常见的一种边坡类型,具有明显的不均匀性和不连续性。土体的蠕变是造成土石混合体边坡变形及失稳的主要原因之一,而目前关于土石混合体边坡稳定性的研究几乎均未考虑土体的蠕变性。首先运用数字图像处理技术对我国某水电站库区的一个土石混合体边坡进行建模,而后利用FLAC3D软件中的强度折减法对其进行稳定性分析,并着重研究了土体的蠕变特征及其蠕变参数对土石混合体边坡变形及稳定性的影响。计算结果表明,土体的蠕变特征会明显降低土石混合体边坡的安全系数,增大边坡变形,进而对边坡的稳定性造成不利影响;其中蠕变黏性系数对边坡的长期稳定性具有较大影响,黏性系数越大,则土石混合体边坡的安全系数越低。     

关于强度折减有限元方法中边坡失稳判据的讨论

在边坡稳定性分析中采用强度折减弹塑性有限元方法时,所得到的总体安全系数在一定程度上依赖于所采用的失稳评判标准,通常以数值计算的收敛性作为边坡失稳判据。而有限元计算的数值收敛性受多种因素的影响,因而由此所得到的安全系数的合理性及其唯一性受到了质疑。为了考查目前各种失稳判据的合理性及其适用性,分别依据计算的收敛性、特征部位位移的突变性和塑性区的贯通性等3个失稳判据,针对某一典型边坡算例,采用强度折减弹塑性有限元方法进行稳定性分析,并与Spencer极限平衡法所得到的总体安全系数进行了对比。对比分析表明,以有限元数值计算的收敛性作为失稳判据在某些情况下所得到的安全系数可能误差较大,而采用特征部位位移的突变性或塑性区的贯通性作为失稳判据所得到的边坡安全系数与Spencer极限平衡法的计算结果比较接近,考虑到实用性与简便性,建议在边坡稳定性分析的强度折减有限元方法中联合采用特征部位位移的突变性和塑性区的贯通性作为边坡的失稳判据。     

有限元强度折减法边坡失稳判据的统一性

目前有限元强度折减法边坡稳定性分析中广泛使用的边坡失稳判据主要有3类,即有限元数值迭代不收敛判据、特征部位位移突变判据、广义塑性应变或等效塑性应变贯通判据。对于上述3类判据的选用,一直以来存在较大的分歧。而从本质上看,3类边坡失稳判据具有一致性和统一性。之所以使用不同判据时计算结果表现出差异,是由于人为误判和有限元数值计算误差造成的。在具有足够计算精度的情况下,3类判据所得结果一致。     

一种基于强度折减法的次级滑动面分析方法研究

传统强度折减法通常对整个边坡区域进行折减,因而只能得到最小安全系数及对应的最危险临界滑动面。但在工程实践中,边坡的治理范围不能局限于最危险滑动面所包围的区域,不满足规范要求的次级滑动面范围内的坡体亦应得到治理。为了克服传统强度折减法在搜索次级滑动面方面的缺陷,提出了一种新的强度折减法。一般而言,边坡的屈服区域主要集中在滑动面两侧附近,形成一个沿滑动面分布的剪切破坏带,因此,在用强度折减法进行边坡稳定性分析时,只需对剪切破坏带范围内的局部坡体进行折减,则可得到该滑动面的安全系数和滑动面。基于以上认识,假设剪切破坏带沿对数螺旋曲线分布,不断调整对数螺旋曲线的形状以得到不同范围的剪切破坏带,对各种不同的剪切破坏带范围内的坡体进行局部折减计算,即可得到不同的安全系数及对应的滑动面。通过两个边坡算例（单台阶、双台阶）验证了该方法的可行性,并最后讨论了剪切破坏带的宽度对结果的影响。两个算例的计算结果表明,该方法不仅可以得到最危险临界滑动面,亦可得到任意安全系数对应的滑动面。     

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

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

基于双强度折减策略的边坡稳定性分析方法探讨

基于能耗分析理论,假定黏聚力和内摩擦角按不同折减系数进行强度折减,采用不同边坡综合安全系数定义方式,根据虚功原理推导了各自综合安全系数的目标函数表达式。采用序列二次规划法和内点法编制了非线性规划迭代程序,对综合安全系数目标函数进行优化求解。依据算例结果探讨了双强度折减技术评价方法对边坡安全系数计算的影响规律。研究结果表明,边坡稳定性影响因素众多,抗剪强度参数c与φ的折减系数之间不存在惟一确定的函数关系,现有预先假定黏聚力和内摩擦角按某一比例进行强度折减的处理方法尚难以完全合理解释;通过强度折减的最短路径来定义边坡综合安全系数的方法物理意义更为明确。研究结论可为工程技术人员加深采用强度折减技术进行边坡稳定性分析的认识提供有益参考。     

基于岩体抗压强度折减的边坡稳定性分析方法

在岩体边坡稳定性分析方法中,基于非线性Hoek-Brown准则的强度折减法,存在折减方案不统一或计算太复杂等问题。针对这一问题,提出了一种基于岩体抗压强度的广义Hoek-Brown准则强度折减法,比直接折减材料参数,更能体现强度衰减的物理意义。方案以岩体单轴抗压强度作为折减依据,对岩石单轴抗压强度s_ci和反映岩体结构特征的参数组合s^a（s、a 均为模型参数）进行同比折减,分析该方案下抗剪切强度和抗拉强度的折减比例,并基于平均抗剪强度定义安全系数。应用该方案,对两个经典边坡算例进行稳定性分析,将所得安全系数、临界滑动面与局部线性化强度折减法、极限平衡法进行对比。结果表明,对于两算例,本方法所得安全系数与其他两种方法结果接近,相对误差低于2%;在算例1中,本方法得到的临界滑动面位置与其他两种方法结果一致性很高;在算例2中,本方案滑动面位置更接近于局部线性化法,且能同时反映剪切和拉伸破坏。以上结果验证了本强度折减方案和安全系数定义法的可靠性和合理性。     

基于局部强度折减法的边坡多滑面分析方法及应用研究

边坡最危险滑动面的搜索方法一直是研究的热点,但边坡内部次级滑动面也可能不满足安全设计要求,因此,考虑边坡多条滑动面的分析方法亦应得到关注。在传统强度折减法中,对边坡整个区域的抗剪强度参数进行折减,此方法仅可得到一个临界滑动面和最小安全系数。提出了一种基于局部强度折减法的多滑面分析方法,即首先定义单元安全系数的概念,并且计算边坡每个单元的安全系数,然后自动搜索出单元安全系数处于不同范围内的单元集合,对各个单元集合的强度参数进行折减计算,即可得到不同安全系数对应的滑动面。通过单台阶和双台阶边坡算例验证了该方法的可行性,结果表明,随着安全系数的增大,潜在滑面的深度和潜在滑动区域亦增大。最后把该方法应用到某隧道进口仰坡的稳定性评价中, 通过该方法得到的多级滑动面与现场监测数据吻合较好。     

极限分析有限元法讲座--Ⅱ有限元强度折减法中边坡失稳的判据探讨

边坡失稳,滑体滑出,滑体由稳定静止状态变为运动状态,同时产生很大的且无限发展的位移,这就是边坡破坏的特征。有限元中通过强度折减使边坡达到极限破坏状态,滑动面上的位移和塑性应变将产生突变,且此位移和塑性应变的大小不再是一个定值,有限元程序无法从有限元方程组中找到一个既能满足静力平衡又能满足应力-应变关系和强度准则的解,此时,不管是从力的收敛标准,还是从位移的收敛标准来判断有限元计算都不收敛。塑性区从坡脚到坡顶贯通并不一定意味着边坡破坏,塑性区贯通是破坏的必要条件,但不是充分条件,还要看是否产生很大的且无限发展的塑性变形和位移,有限元计算中表现为塑性应变和位移产生突变。在突变前计算收敛,突变之后计算不收敛,表征滑面上土体无限流动,因此可把有限元静力平衡方程组是否有解,有限元计算是否收敛作为边坡破坏的依据。-     

基于FLAC强度折减法的江苏田湾核电站岩质高边坡稳定性分析

为了评价田湾核电站岩质高边坡的稳定性,应用有限差分强度折减法,对该核电站三、四期工程人工边坡的稳定性进行了计算分析,并与传统的极限平衡法的计算结果进行了对比。结果表明,由于有限差分强度折减法克服了极限平衡法未能考虑材料本构关系的不足,因而更能贴近实际模拟分析边坡的稳定状态,求得的边坡稳定系数更为合理。     

基于FLAC强度折减法的江苏田湾核电站岩质高边坡稳定性分析

为了评价田湾核电站岩质高边坡的稳定性,应用有限差分强度折减法,对该核电站三、四期工程人工边坡的稳定性进行了计算分析,并与传统的极限平衡法的计算结果进行了对比。结果表明,由于有限差分强度折减法克服了极限平衡法未能考虑材料本构关系的不足,因而更能贴近实际模拟分析边坡的稳定状态,求得的边坡稳定系数更为合理。