非线性Mohr—Coulomb破坏准则下边坡可靠度上限

贺志军　曹吉　赵炼恒　瞿召乾　杨胜博

摘要：传统边坡可靠度分析往往在岩土参数服从线性Mohr-coulomb（简称线性M-c）破坏准则的假设条件下进行，并且常常采用极限平衡法或有限元法计算安全系数。然而，岩土介质破坏准则具有一定的非线性。为能更加实际地描述岩土破坏机理和得到严格精确的解，基于非线性Mohr-coulomb（简称非线性M-c）破坏准则，结合极限分析上限法和蒙特卡洛法，进行边坡可靠度上限分析。当非线性参数m=1时，与等效的线性M-c破坏准则进行对比计算，验证了方法的可行性。同时，将初始粘聚力、内摩擦角arctan（C_{o/σt）和非线性参数作为随机变量且服从截断正态分布，进行了参数变异性和敏感性影响分析。研究表明：非线性M-c破坏准则下，边坡可靠度随初始粘聚力、内摩擦角arctan（Co/σt）和非线性参数变异性的增大而减小；边坡可靠度随初始粘聚力和内摩擦角arctan（Co/σt）的增大而增大，随非线性参数的增大而减小。}

关键词：边坡；可靠度；破坏准则；极限分析；蒙特卡洛法

中图分类号：TU 457

文献标志码：A

文章编号：1674-4764（2016）06-0001-09

边坡稳定性问题一直是岩土工程的一个重要研究内容。目前，边坡稳定性分析主要有两大体系：确定性体系与不确定性体系（可靠度体系）。确定性体系使用极限平衡法、数值模拟法或极限分析法等方法分析边坡的稳定性，求得边坡最小安全系数，以此作为边坡稳定性评价指标。然而，边坡是一个极其复杂的系统，岩土参数具有明显的随机性，采用确定性体系分析边坡稳定性不符实际。文献[2-3]也指出：由于安全系数没有考虑参数随机性和离散型对结果的影响，导致实际工程中很多结构在满足安全系数的条件下依然出现了破坏现象。以概率论为基础的可靠度体系可考虑边坡系统内部的随机关系，可给出边坡稳定程度，因而可以弥补用单一安全系数分析边坡稳定性的局限性。

边坡可靠度分析主要有两大步骤：一是构建计算边坡稳定性安全系数的模型（功能函数）；二是使用可靠度分析方法计算边坡可靠度（失效概率或可靠度指标）。在构建模型方面，目前应用最广泛的是极限平衡法，该方法理论简单、易于实施，但所作假设较多，根据塑性理论可知，该方法所获解答不是严格的上下限解。另一应用广泛的方法是数值方法，该方法将有限元技术应用到边坡稳定性分析中，可以考虑土体与其中结构物的共同作用，但其所得到的极限荷载值仍不够精确。相比上述两种方法，极限分析法可以得到边坡极限荷载的严格上限解，在此基础上进行可靠度分析，可以得到严格边坡可靠度上限值，这对于分析边坡稳定性可能具有重要意义。在可靠度分析方法方面，常用的分析方法有一次二阶矩法、JC法、MonteCarlo法等。其中，Monte Carlo法被认为是一种相对精确的方法，根据大数定律，只要抽样次数足够大，其精度就能足够高。目前，众多学者应用可靠度分析理论对边坡稳定性进行研究均是在线性M-C破坏准则假设下进行的。而事实上，岩土介质服从非线性破坏准则，线性破坏准则只是一种特例：将更为符合实际的略微弯曲岩土材料强度线简化成直线形状，该方法虽简单易于分析，但无法准确表述岩土强度特性。因此，考虑岩土破坏准则为非线性情况下的边坡可靠度研究显得十分必要。1987年，Zhang等提出了幂函数非线性破坏准则，尔后，大量学者对基于该破坏准则下的边坡稳定性问题进行了深入研究，非线性M-C破坏准则得到快速发展且已较为成熟。但多年来，鲜见基于非线性M-C破坏准则下边坡可靠度的深入研究。

基于以上考虑，在非线性Mohr-Coulomb破坏准则下，结合外切线技术和强度折减技术计算边坡安全系数上限解，进一步运用蒙特卡洛法计算边坡可靠度。通过与基于线性M-C破坏准则的计算结果进行对比，验证了方法的可行性。并视初始粘聚力C_{o、内摩擦角arctan（Co/σt）与非线性参数m为随机参数且服从截断正态分布，研究非线性M-C破坏准则下各参数变异性、敏感性对均质各向同性边坡可靠度的影响。}

1.非线性M-C破坏准则及抗剪强度参数引入方法

2.基于强度折减技术的边坡极限分析上限法

边坡可靠度分析中，安全系数计算式表示了边坡安全系数与土工参数之问的关联关系，作为功能函数的构建主体至关重要。极限分析上限法基于虚功率原理推导，根据外力做功和内部耗能相等原理获得目标函数并根据能量耗散最小化原理获得极限荷载的最小值，进一步结合强度折减技术可获得严格精确的安全系数上限解。

选取符合简单边坡的对数螺旋面破坏机构作为破坏模式，以通过坡趾下的对数螺旋线旋转间断机构为例进行分析，如图2。刚性块体ABCCA绕旋转中心O相对对数螺旋面BC以下的静止材料作刚体旋转，BC'是速度问断面。ABCCA区重力做功率为外荷载功率；问断面BC上滑面耗损率为

在已知边坡几何尺寸及岩土材料参数条件下，边坡安全系数F_{a可由θCh、θCD、β、φt等4个未知参数确定，且安全系数计算公式是Fs的隱函数。将Fs作为目标函数，通过非线性数学规划方法，利用Matlab软件，可求得边坡最小安全系数Fs。}

3.基于非线性M-C破坏准则的边坡可靠度上限计算

3.1岩土随机参数的选定及其分布类型

线性M-C破坏准则下岩土参数分布类型已有大量研究成果。马建全等认为不同环境下岩土参数具有不用分布类型，同时，比较了岩土参数服从正态分布、对数正态分布等不同分布类型下的边坡可靠度大小；张继周等从概率分布类型的产生背景、所描述对象的物理意义人手，研究各分布类型对可靠度分析的影响，得出正态分布和对数正态分布较为合理的结论；其他众多边坡可靠度研究亦在岩土参数服从正态分布下进行。因此，笔者也选取正态分布作为岩土随机参数的分布类型。

线性M-C破坏准则下，岩土参数c与φ服从正态分布。由式（2）可知，当m=1时，非线性M-C破坏准则指数型表达式中的参数C_o与无量纲参数C_o/σ_t的反正切值arctan（C_o/σ_t）服从正态分布。当m≠1时，假设参数 C_o与内摩擦角arctan（C_o/σ_t）服从正态分布，并进一步考虑非线性参数的变异性，假设m亦服从正态分布。由此，选取初始粘聚力c。内摩擦角arctan（C_o/σ_t）与非线性参数m作为随机参数。

桂勇等在边坡稳定二元体系的建立中，采用同时考虑材料指标的统计分布和区间分布的方法，即截断分布；Johari等在基于毕绍普法的地震边坡稳定概率模型建立中，采用截断正态分布作为随机变量的概率分布。截断分布可以考虑岩土参数具体变化范围，更符合工程实际。选取截断的正态分布作为岩土随机参数的分布类型，随机参数概率密度函数为

4.对比分析

目前，鲜有非线性M-C破坏准则方面的边坡可靠度研究。为验证方法的可行性，基于线性M-C破坏准则将边坡可靠度上限解与Johari等、蔡宁等及吕杨的计算结果进行对比分析。参数信息如表1所示，对比结果如表2所示。

由表2可知，线性M-C破坏准则下本文结算结果与Johari等、蔡寧等及吕杨的计算结果失效概率P_F最大相差5.33%。最小相差0.90%，平均相差2.455%；可靠指标口最大相差1.77%，最小相差0.57%，平均相差1.105%，表明该方法可行。

5.算例分析

选取一典型边坡作为研究对象，如图3所示。研究依据文献确定岩土材料参数C_o和σ_t获得参数c。与内摩擦角arctan（C_o/σ_t）的均值，分别为90kN/m°与20°（σ_t为247.3 kN/m。）。并假设坡角α=60。坡高H=25m，岩土容重y=20kN/m³，非线性参数均值um=1.6。随机参数变异系数取值情况为：1）参数C_o与内摩擦角arctan（C_o/σ_t）沿用Johari等的c与ψ的变异系数，分别为0.20与0.10。2）非线性参数m对于不同岩土材料有一定变化范围，因此，非线性参数m取较小变异系数符合实际工程，取值为0.05。截断的正态分布随机参数如表3所示。

5.1参数变异性影响

为分析非线性M-C破坏准则下随机参数（以下简称参数）的变异性对边坡可靠度的影响，进行参数变异性影响分析。当参数变异系数的增加值δ=0、COV/3、2COV/3、3COV/3、4COV/3、5COV/3、6COV/3时，边坡失效概率Pf的变化情况如图4所示；边坡可靠指标β的变化情况如图5所示。对应（δ=0、δ=COV）的安全系数统计结果对比如图6～7所示；对应（δ=0、δ=COV）的可靠度相对变化如表4所示。

由图4～7与表4可知，初始粘聚力C_o的变异性对边坡可靠度具有较大影响；内摩擦角arctan（C_o/σ_t）与非线性参数m的变异性对边坡可靠度的影响较小。随着初始粘聚力变异系数的增大，边可靠度减小得较为明显，且这种变化幅度随着变异系数的增大逐渐减小；随着内摩擦角arctan（C_o/σ_t）与非线性参数变异系数的增大，边坡可靠度减小且较为缓慢。

对于该算例，岩土材料参照文献[12]选取，其初始粘聚力为90kN/m^{2初始内摩擦角arctan（C_o/σ_t）为20°。已有研究表明，结合上限理论与外切线技术后，所获得的瞬时粘聚力C_t随非线性参数的增大而增大，而瞬时内摩擦角ψ_t随非线性参数的增大而减小，可知该算例中粘聚力参数对边坡稳定性的影响占主导作用。同时，加之在计算过程中选取了较大的变异系数（20%），使得初始粘聚力的变异性对边坡可靠度影响较其它参数大得多。}

5.2敏感性分析

为研究在非线性M-C破坏准则下随机参数（以下简称参数）的变化对边坡可靠度的影响，进行参数敏感性分析。当参数均值的增加值δ=0、std/3、2std/3、3std/3、4std/3、5std/3、6std/3时，边坡失效概率Pf的变化情况如图8所示；边坡可靠指标β的变化情况如图9所示。对应（δ=0、δ=td）的安全系数统计结果对比如图10～11所示；对应（δ=0、δ=std）的可靠度相对变化如表5所示。

由图8～11和表5可知，初始粘聚力C_o的均值对边坡可靠度有较大影响；无量纲参数arctan（C_o/σ_t）与非线性参数m的均值对边坡可靠度的影响相对较小。随着初始粘聚力和无量纲参数arctan（C_o/σ_t）均值的增大，抗剪强度参数C_t和ψ_t均增大，边坡可靠度增大；随着非线性参数均值的增大，抗剪强度参数C_t和ψ_t发生变化，边坡可靠度减小。

6.结论

基于非线性M-C破坏准则，结合外切线技术和强度折减技术，采用极限分析上限法求得边坡安全系数上限解，进一步运用蒙特卡洛法计算边坡的可靠度（失效概率Pf与可靠指标β）。在线性M-C破坏准则下，通过与已有算例对比计算，验证了方法的可行性。同时，基于非线性破坏准则分析参数变异性和参数敏感性对边坡可靠度的影响，得到如下结论：

1）边坡可靠度随初始粘聚力C_o、内摩擦角arctan（C_o/σ_t）和非线性参数m变异系数的增大而减小。

2）边坡可靠度随初始粘聚力c。和内摩擦角arctan（C_o/σ_t）均值的增大而增大，随非线性参数m均值的增大而减小。

工程实际中，岩土材料参数分布形态众多，笔者研究的基于非线性M-C破坏准则的边坡可靠度上限分析仅采用截断正态分布这一种，其他分布形态的研究有待深入；岩土材料参数取值范围和变异系数差异较大，研究的基于非线性M-C破坏准则的边坡可靠度上限分析采用了已有文献中参数分析的取值，符合工程实际的其他参数取值范围和变异性的研究有待深入。