CN107330145B - 一种考虑转动效应的岩质边坡极限承载力分析下限法 - Google Patents

一种考虑转动效应的岩质边坡极限承载力分析下限法 Download PDF

Info

Publication number
CN107330145B
CN107330145B CN201710381897.0A CN201710381897A CN107330145B CN 107330145 B CN107330145 B CN 107330145B CN 201710381897 A CN201710381897 A CN 201710381897A CN 107330145 B CN107330145 B CN 107330145B
Authority
CN
China
Prior art keywords
structural
rock
slope
lower limit
structural surface
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710381897.0A
Other languages
English (en)
Other versions
CN107330145A (zh
Inventor
李泽
彭普
周宇
张小艳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN201710381897.0A priority Critical patent/CN107330145B/zh
Publication of CN107330145A publication Critical patent/CN107330145A/zh
Application granted granted Critical
Publication of CN107330145B publication Critical patent/CN107330145B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • E—FIXED CONSTRUCTIONS
    • E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
    • E02D—FOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
    • E02D1/00—Investigation of foundation soil in situ
    • E02D1/02—Investigation of foundation soil in situ before construction work
    • E—FIXED CONSTRUCTIONS
    • E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
    • E02D—FOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
    • E02D17/00—Excavations; Bordering of excavations; Making embankments
    • E02D17/20—Securing of slopes or inclines
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/23—Dune restoration or creation; Cliff stabilisation

Abstract

本发明涉及一种同时考虑岩体平动和转动效应的承载力分析下限法，属于岩质边坡稳定性分析技术领域。本发明采用刚性块体单元将其离散为刚性岩块+结构面的几何系统，以结构面的法向力、剪力、弯矩为未知量，构建同时满足刚性块体单元考虑平动、转动效应的平衡方程、结构面剪切屈服条件、结构面拉伸屈服条件和结构面转动屈服条件以及静力边界条件的静力许可应力场，以超载系数为目标函数，建立下限法线性数学规划模型；采用内点算法对线性数学规划模型求解，同时采用迭代算法求解强度储备系数，获得节理岩质边坡的超载系数或强度储备系数的下限解和对应的屈服区。本发明概念明确、计算精度高，可应用于节理岩质边坡平动或转动破坏的承载力分析。

Description

一种考虑转动效应的岩质边坡极限承载力分析下限法
技术领域
本发明是一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法，特别涉及一种同时考虑岩体结构面剪切破坏、拉伸破坏和转动破坏的承载力分析法。属于岩质边坡稳定性分析技术领域。
背景技术
节理岩质边坡的失稳是一种常见的地质灾害。在国内外公路、铁路、矿山和水利水电等工程中均碰到大量节理岩质边坡失稳破坏的灾难性事故，如国内天生桥二级水电站边坡倾倒破坏、金川露天矿边坡失稳破坏等等。节理岩质边坡的承载力和稳定性评价方法已成为了岩土工程界必须解决的重要技术问题。
节理岩质边坡由于内部存在着大量的节理、裂隙，使得其力学特性非常复杂。节理岩质边坡一般认为是由岩块和结构面组成的复杂系统，岩块的运动方式一般可分为平动和转动，结构面的破坏模式也呈现多样化。岩块平动滑移破坏问题相对简单，岩块转动倾倒破坏问题则比较复杂，一般岩块的转动倾倒破坏可分为块状倾倒、弯曲倾倒和块状弯曲倾倒等类型。自20世纪50年代以来，岩体平动破坏和转动破坏的研究就已经被重视。为了解决岩体平动、转动倾倒的稳定性分析问题，基于极限平衡原理的边坡稳定计算方法被提了出来，并且被广泛应用于工程中。这种方法一般对岩块变形失稳模式做一些假定，把岩质边坡的滑体分为稳定区、倾倒区或滑动区，但对一些特殊问题其合理性有待验证。其后，很多数值方法被提出来用于岩体平动、转动破坏失稳的分析计算，比如有限单元法、离散单元法、块体单元法、DDA等。同时，岩质边坡平动、转动的破坏物理模型试验也被广泛应用于工程实践。
就目前的研究成果来说，同时考虑平动、转动效应的节理岩质边坡的承载力分析领域还存在一些有待于进一步探讨的问题。具体不足表现在：
(1)很多分析方法不能同时考虑岩体的平动和转动效应，比如刚体极限平衡法、有限元法等一般只能考虑岩体的平动效应，但是工程实践表明岩体的转动效应是非常重要、不可或缺的，比如：对于顺层节理岩质边坡而言转动破坏是其主要破坏模式。
(2)节理岩质边坡中的结构面破坏的模式也非常多样，一般认为包括三种：滑移剪切破坏、拉伸破坏和转动破坏，但现有的分析方法很难同时考虑这三种破坏模式。
(3)在极限分析领域，特别是下限法分析领域，研究岩体同时考虑平动、转动效应的研究成果较少。
本发明基于国家自然科学基金项目(资助编号：51564026)的研究工作，将塑性极限分析下限法理论、刚性块体单元离散技术以及数学规划手段结合起来，提出一种同时考虑平动、转动效应的岩质边坡极限承载力分析的下限法。
发明内容
本发明的目的是提供一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法，为节理岩质边坡稳定性分析提供一种新方法。
本发明的基本原理是：以节理岩质边坡为研究对象，采用刚性块体单元将其离散为刚性岩块+结构面的几何系统，以结构面的法向力、剪力、弯矩为未知量，基于塑性极限分析下限法理论，构建同时满足：刚性块体单元考虑平动、转动效应的平衡方程、结构面剪切屈服条件、结构面拉伸屈服条件和结构面转动屈服条件以及静力边界条件的静力许可应力场，然后以荷载的超载系数为目标函数，建立同时考虑平动和转动效应的节理岩质边坡承载力分析的下限法线性数学规划模型；然后采用内点算法对线性数学规划模型进行求解，同时可采用迭代算法求解强度储备系数，最终可获得节理岩质边坡的超载系数或强度储备系数的下限解和对应的屈服区。
本发明的技术路线如图1所示。
本发明的同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法的技术方案依次按以下步骤进行：
一、拟定节理岩质边坡的计算参数
根据节理岩质边坡的实际情况，拟定下限法分析需要的计算参数，主要包括：岩质边坡的几何参数、节理参数、岩体和节理材料的参数(容重、凝聚力、摩擦角、抗拉强度等)、边坡荷载参数信息。
二、采用刚性块体单元离散节理岩质边坡
本发明采用刚性块体单元离散节理岩质边坡，刚性块体单元受力如图2所示，相邻单元的结构面受力如图3所示。其中(x,y)为总体坐标系；(nk,sk)为相邻块体单元i、j之间结构面k的局部坐标系，nk为结构面的外法线方向，sk为结构面的切向；刚性块体单元i形心Ci上作用有等效荷载力向量
Figure GDA0002444805220000021
fxi为沿x方向的外力，fyi为沿y方向的外力，mi为外力弯矩(规定逆时针为正向)；相邻块体单元i、j之间的结构面k形心Pk上作用有内力向量
Figure GDA0002444805220000022
另一侧作用有其反作用力向量
Figure GDA0002444805220000023
Nk是沿结构面外法线nk方向的法向力(拉正压负)，Vk是沿结构面切向sk方向的剪力(使块体产生逆时针转动为正)，Mk为结构面的内力弯矩(规定逆时针为正向)。
采用刚性块体单元离散节理岩质边坡以后，边坡变成刚性块体单元+结构面的几何系统，本发明以刚性块体单元结构面的形心Pk上作用的内力向量Qk＝(Nk,Vk,Mk)T为未知量构建岩质边坡的静力许可应力场。为了简化计算，本发明作如下假设：(1)假设岩块为刚体，岩质边坡的破坏只会发生在相邻块体单元之间的结构面上；(2)假设结构面材料为理想刚塑性材料，(3)假设结构面有可能发生剪切破坏、拉伸破坏或转动破坏，(4)岩块同时发生平动或转动。
三、建立同时考虑岩体平动、转动效应的节理岩质边坡承载力分析下限法数学规划模型
本发明的目的是求解节理岩质边坡超载系数或强度储备系数的下限解，需要建立静力许可应力场。根据下限定理，要建立边坡的静力许可应力场需满足平衡方程、屈服条件和边界条件。
1、目标函数
本发明方法将节理岩质边坡的超载定义为：
K1＝Fc/Fa (1)
上式中，Fc是极限荷载，Fa是边坡结构当前实际施加的外荷载，对于岩质边坡，求极限状态时其超载的外荷载可以是自重也可以是面力或集中力。
根据下限定理，岩质边坡在达到极限状态时，需求超载系数的最大值。本发明将超载系数设为目标函数如下：
Maximize:K1 (2)
2、刚性块体单元平衡方程
设整个节理岩质边坡被离散为n个刚性块体单元，
Figure GDA0002444805220000036
为任意刚性块体单元i上的结构面数量。对于图2所示边坡中的任意刚性块体单元i，其受各种力的作用而保持平衡状态，本发明同时考虑块体单元的平动、转动效应，因此同时考虑力的平衡和力矩的平衡。对于单一刚性块体单元i其力和力矩平衡方程为：
Figure GDA0002444805220000031
上式中：n是边坡中块体的数量；
Figure GDA0002444805220000032
是刚性块体单元块体i上的结构面数量；
Figure GDA0002444805220000033
是结构面k形心Pi上作用的内力向量；
Figure GDA0002444805220000034
是刚性块体单元i形心Ci作用的等效荷载力向量；Tk是结构面k的局部坐标(nk,sk)到总体坐标系(x,y)的转换矩阵，
Figure GDA0002444805220000035
αk是nk轴与x轴的夹角，逆时针为正；δ1是图2中结构面k形心Dk到点Pk的距离；δ2是图2中岩块i中Dk点到形心Ci到的距离；K1为超载系数；K1为超载系数。
对于含有n个块体的二维岩质边坡，式(3)包括2n个力的平衡方程和n个力矩平衡方程，可采用矩阵、向量的形式简写如下：
Figure GDA0002444805220000041
上式中：C是考虑平动、转动效应的整体平衡矩阵，
Figure GDA0002444805220000042
是所有结构面的内力向量，
Figure GDA0002444805220000043
是所有块体单元形心的等效荷载力向量。
3、结构面的屈服条件
本发明假设节理岩质边坡的破坏仅仅发生在结构面上。在同时考虑岩块的平动、转动力学效应时，任意相邻块体单元之间的结构面会发生三种破坏模式(如图4所示)，包括剪切破坏，拉伸破坏和转动破坏。
对于任意一条结构面k(如图3所示),在局部坐标系(nk,sk)中其形心上作用有内力向量
Figure GDA0002444805220000044
则结构面的剪切破坏Mohr-Coulomb屈服条件为：
Figure GDA0002444805220000045
上式中：ns是整个边坡中结构面的数量，lk是结构面k的长度，c是结构面的凝聚力，
Figure GDA0002444805220000046
是结构面的内摩擦角，Nk是结构面k的法向力(规定以拉力为正)，Vk是结构面k的剪力。
结构面的抗拉屈服条件为：
Nk-σT/lk≤0,k＝(1,…，ns) (6)
上式中：ns为整个边坡中结构面的数量，lk为结构面k的长度，Nk为结构面k的法向力(规定以拉力为正)，σT为结构面的抗拉强度，由于贯通的节理一般不能承受拉应力，因此在本发明中结构面的抗拉强度取σT＝0。
结合结构面的剪切屈服条件式(5)和抗拉屈服条件式(6)，则可得到结构面的拉剪破坏的修正Mohr-Coulomb屈服条件如图5所示。
当相邻块体单元的结构面发生转动破坏时，块体单元发生转动极限状态下的受力如图6所示，由于内力弯矩Mk的方向有两种可能的方向(顺时针和逆时针)，因此块体单元有可能分别绕结构面的两个端点P1、P2发生两种转动破坏。结构面的转动屈服条件如图7所示，可写为：
Figure GDA0002444805220000047
上式中：ns是整个边坡中结构面的数量，lk是结构面k的长度,Nk是结构面k的法向力(规定以拉力为正),Mk是结构面k的内力弯矩。
综合式(5)、式(6)、式(7)，得到考虑块体单元平动、转动力学效应的广义屈服条件：
Figure GDA0002444805220000051
上式中：k＝(1,…,ns)，
Figure GDA0002444805220000052
为以结构面内力向量
Figure GDA0002444805220000053
表示的结构面k的广义屈服条件。
式(9)共计为5ns个约束条件，5ns个约束均为线性约束条件。将广义屈服条件式(8)简写如下：
Figure GDA0002444805220000054
上式中：f为广义屈服条件，
Figure GDA0002444805220000055
为所有的结构面内力向量。
4、静力边界条件
根据下限定理，静力许可应力场必须满足已知的静力边界条件。考虑节理岩质边坡中已知的边界条件为
Figure GDA0002444805220000056
的界面b，其边界条件表达式为：
Figure GDA0002444805220000057
上式中：
Figure GDA0002444805220000058
是已知边界条件的结构面的内力向量，
Figure GDA0002444805220000059
是已知边界条件结构面上作用的力和力矩。
5、同时考虑岩体平动、转动效应的节理岩质边坡承载力分析下限法数学规划模型
根据下限定理，结合目标函数式(2)、平衡方程式(4)、广义屈服条件式(10)、边界条件式(11)，则可得到同时考虑岩体平动、转动效应的节理岩质边坡承载力分析下限法线性数学规划模型：
Figure GDA00024448052200000510
上式为一个线性数学规划模型，通过优化算法求解可到节理岩质边坡的超载系数和对应的屈服区。
四、优化算法求解节理岩质边坡的极限承载力
通过上述步骤建立了同时考虑岩体平动、转动效应的节理岩质边坡极限承载力分析的下限法数学规划模型(式(11))，其是一个大规模线性数学规划模型。本发明采用Fortran12.0编写下限法计算程序，计算程序包括前处理、构建下限法数学规划模型、优化求解、后处理等4个部分。本发明采用内点算法对线性数学规划模型进行求解。最终可获得节理岩质边坡稳定性的超载系数的下限解以及对应屈服区。
对于节理岩质边坡，有时候除了超载系数还需要求解强度储备系数，本发明的将强度储备系数定义为：
Figure GDA0002444805220000061
上式中：
Figure GDA0002444805220000062
为结构面的凝聚力和摩擦角，
Figure GDA0002444805220000063
为结构面进行强度折减以后的凝聚力和摩擦角，σT,σT'为岩体折减前、后的抗拉强度。
如果直接求解强度储备系数K2，式(9)会变成非线性约束，从而导致式(11)为非线性规划模型。为了避免求解非线性规划模型，本发明基于式(11)采用迭代法求解强度储备系数K2，具体方法是：在式(12)中逐步降低材料强度参数并带入式(9)，然后使用式(11)循环求解K1，当K2使得K1充分接近于1时(误差小于0.001)即为所求的强度储备系数的最大值(具体迭代过程详见图1)。
本发明的特点是：本发明方法为一种求解节理岩质边坡强度储备系数的严格的下限法，该方法将下限定理、刚性块体单元离散技术、线性数学规划有效结合起来，以相邻刚性块体单元之间结构面的法向力、剪力和弯矩为未知量，构建了下限法的静力许可应力场，静力许可应力场严格满足三个条件：(1)力和力矩的平衡方程、(2)结构面剪切破坏、拉伸破坏和转动破坏的广义屈服条件、(3)静力边界条件，并采用内点算法对线性数学规划模型进行求解，得到了节理岩质边坡同时考虑岩块平动、转动效应、同时考虑结构面剪切破坏、拉伸破坏和转动破坏三种破坏模式的极限状态，获得了对应的超载系数或强度储备系数的下限解以及对应的屈服区。
本发明具有以下有益效果：
1、本发明方法为一种求解节理岩质边坡强度储备系数的严格的下限法，可同时考虑岩块的平动、转动效应，是一种较为接近节理岩质边坡实际变形破坏形态的分析方法。
2、本发明采用同时考虑了节理岩质边坡中结构面的剪切破坏、拉伸破坏以及转动破坏，可获得三种破坏模式的屈服区，对了解边坡的破坏形态更为直观。
3、本发明方法构建的下限法模型为线性规划数学模型，可直接求解得到边坡超载系数的下限解，同时可采用迭代法求解强度储备系数的下限解，数学模型优化求解简单、准确。
4、本发明方法概念明确、计算精度高，相对于传统的刚体极限平衡法更适用于复杂节理岩质边坡的稳定性分析。
附图说明
图1为本发明技术路线图；
图2为刚性块体单元受力示意图；
图3为相邻刚性块体单元之间结构面受力示意图；
图4为相邻刚性块体单元之间结构面破坏模式：剪切破坏、拉伸破坏和转动破坏；
图5为结构面剪切拉伸破坏的修正Mohr-Coulomb屈服准则；
图6为块体单元倾倒破坏时结构面极限状态受力示意图；
图7为结构面转动破坏屈服条件示意图；
图8为实施例1节理岩质边坡的几何形状示意图；
图9为实施例1节理岩质边坡刚性块体单元离散示意图；
图10为实施例1节理岩质边坡极限状态下屈服区示意图；
图11为实施例2节理岩质边坡的几何形状示意图；
图12为实施例2节理岩质边坡刚性块体单元离散示意图；
图13为实施例2节理岩质边坡极限状态下屈服区示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
实施例1：实施例1为使用式(11)求解一个含有2个岩块的节理岩质边坡的极限外荷载。
(一)拟定节理岩质边坡的计算参数。
图8为实施例1节理岩质边坡的几何形状示意图。此边坡包含两个岩块，岩块1的形状为矩形，其宽1.0m、高2.0m；岩块2的形状为三角形，高度2.0m、底宽5.495m。岩块1与岩块2之间为一贯通节理，节理面倾角为20°。边坡材料的物理力学参数详见表1。在岩块1的形心C1处作用有一个水平向左的外力F，本实施例的目的就是求解F的极限值Fc。
表1实施例1节理岩质边坡材料的物理力学参数表
材料名称 重度/(kN·m<sup>-3</sup>) 黏聚力/kPa 内摩擦角/(°) 抗拉强度/kPa
完整岩块 25 4000.0 50.0 1000
节理面 / 1.0 30.0 0.0
(二)采用刚性块体单元离散节理岩质边坡
采用刚性块体单元离散实施例1的节理岩质边坡，并以刚性块体单元结构面的内力(法向力、剪力、弯矩)为未知量。实施例1的边坡系统共离散为2个刚性块体单元(如图9所示)。
(三)建立同时考虑岩体平动、转动效应的实施例1边坡承载力分析的下限法数学规划模型
采用刚性块体单元离散的实施例1的节理岩质边坡以后，根据式(12)建立了求解同时考虑岩体平动、转动效应的边坡承载力下限法模型。
(四)优化算法求解节理岩质边坡的极限承载力
根据式(11)建立的下限法模型，采用编制的线性规划优化求解程序，求解边坡中块体1形心处的外荷载F的极限值。计算结果列于表2中。
表2中计算了实施例1在两种失稳模式下的极限荷载，第一种是只考虑块体平动效应的剪切破坏，第二种是同时块体考虑平动+转动效应，并同时考虑结构面的剪切破坏、拉伸破坏、转动破坏。在第一种失稳模式下本发明方法计算得到的极限荷载为9673.70N，解析解为9695.73N；在第二种失稳模式下本发明方法计算得到的极限荷载为5742.61N，解析解为5754.29N；本发明方法得到的结果均小于解析解，属于下限解，且两者误差非常小，验证了本发明方法的正确性。另外，第一种失稳模式的极限荷载大于第二种失稳模式的极限荷载，说明对于此岩质边坡如果仅仅考虑平动的剪切破坏，将会高估边坡的承载能力，即：转动破坏会先于剪切破坏发生。图10为实施例1在第二种失稳模式下本发明方法计算得到的屈服区示意图，两个块体交界面的节理面发生剪切屈服、同时绕P1点发生转动屈服。本发明方法同时考虑平动、转动效应，相比传统方法能更准确的评估岩质边坡的承载力。
表2实施例1节理岩质边坡极限荷载计算结果
Figure GDA0002444805220000081
实施例2：实施例2为使用式(12)求解一个含有6个岩块的节理岩质边坡的强度储备系数。
(一)拟定节理岩质边坡的计算参数。
图11为实施例2节理岩质边坡的几何形状示意图。此边坡底面倾角25°，边坡包含6个岩块，每个岩块的几何尺寸详见图11。边坡材料的物理力学参数详见表3。6个岩块中，岩块1为固定岩块、不能移动，岩块2～岩块6可能发生平动或转动。本实施例的目的是求解岩质边坡的强度储备系数。
表3实施例2节理岩质边坡材料的物理力学参数表
材料名称 重度/(kN·m<sup>-3</sup>) 黏聚力/kPa 内摩擦角/(°) 抗拉强度/kPa
完整岩块 26 5000.0 55.0 1500
节理面 / 0.0 35.0 0.0
(二)采用刚性块体单元离散节理岩质边坡
采用刚性块体单元离散实施例2的节理岩质边坡，并以刚性块体单元结构面的内力(法向力、剪力、弯矩)为未知量。实施例2的边坡系统共离散为6个刚性块体单元，如图12所示。
(三)建立同时考虑岩体平动和转动效应的实施例2边坡承载力分析的下限法数学规划模型
采用刚性块体单元离散的实施例2的节理岩质边坡以后，根据式(12)建立了求解同时考虑岩体平动、转动效应的边坡承载力下限法模型。
(四)优化算法求解节理岩质边坡的极限承载力
由于本实施例需要求解边坡的强度储备系数，因此根据式(12)建立的下限法模型，采用线性规划优化求解程序用迭代法求解强度储备系数。最终求得的强度储备系数为K2＝1.435，相对应的失稳模式是岩块3、岩块4、岩块5、岩块6发生逆时针方向的转动屈服破坏，如图13所示。由于块体1是固定不动的，因此本实施例如果只考虑岩块的平动效应，则此边坡不会发生失稳。但采用本发明方法计算的结果是：当强度储备系数取K2＝1.435时，边坡发生逆时针的倾倒破坏。

Claims (2)

1.一种考虑转动效应的岩质边坡极限承载力分析下限法，其特征在于：以节理岩质边坡为研究对象，采用刚性块体单元将其离散为刚性岩块+结构面的几何系统，以结构面的法向力、剪力、弯矩为未知量，基于塑性极限分析下限法理论，构建同时满足刚性块体单元考虑平动、转动效应的平衡方程、结构面剪切屈服条件、结构面拉伸屈服条件和结构面转动屈服条件以及静力边界条件的静力许可应力场，然后以荷载的超载系数为目标函数，建立同时考虑平动和转动效应的节理岩质边坡承载力分析的下限法线性数学规划模型；然后采用内点算法对线性数学规划模型进行求解，同时采用迭代算法求解强度储备系数，最终获得节理岩质边坡的超载系数或强度储备系数的下限解和对应的屈服区；具体步骤如下：
(一)拟定节理岩质边坡的计算参数
根据节理岩质边坡的实际情况，拟定下限法分析需要的计算参数，包括：岩质边坡的几何参数、节理参数、边坡荷载参数信息及岩体和节理材料的参数，所述岩体和节理材料的参数容重、凝聚力、摩擦角、抗拉强度；
(二)采用刚性块体单元离散节理岩质边坡
采用刚性块体单元离散节理岩质边坡，边坡变成刚性块体单元+结构面的几何系统，(x,y)为总体坐标系；(nk,sk)为相邻块体单元i、j之间结构面k的局部坐标系，nk为结构面的外法线方向，sk为结构面的切向；刚性块体单元i形心Ci上作用有等效荷载力向量
Figure FDA0002444805210000011
fxi为沿x方向的外力，fyi为沿y方向的外力，mi为外力弯矩，规定逆时针为正向；相邻块体单元i、j之间的结构面k形心Pi上作用有内力向量
Figure FDA0002444805210000012
另一侧作用有其反作用力向量
Figure FDA0002444805210000013
Nk是沿结构面外法线nk方向的法向力，拉正压负，Vk是沿结构面切向sk方向的剪力，使块体产生逆时针转动为正，Mk为结构面的内力弯矩，规定逆时针为正向；
以刚性块体单元结构面的形心Pk上作用的内力向量
Figure FDA0002444805210000014
为未知量构建岩质边坡的静力许可应力场；为了简化计算，作如下假设：(1)假设岩块为刚体，岩质边坡的破坏只会发生在相邻块体单元之间的结构面上(2)假设结构面材料为理想刚塑性材料(3)假设结构面有可能发生剪切破坏、拉伸破坏或转动破坏(4)岩块同时发生平动或转动；
(三)建立同时考虑岩体平动和转动效应的节理岩质边坡承载力分析下限法数学规划模型
根据下限定理，要建立边坡的静力许可应力场需满足平衡方程、屈服条件和边界条件；
(1)目标函数
将超载系数设为目标函数如下：Maximize:K1，节理岩质边坡的超载系数定义为：K1＝Fc/Fa；
上式中，Fc是极限荷载，Fa是边坡结构当前实际施加的外荷载，对于岩质边坡，求极限状态时其超载的外荷载为自重、面力或集中力；
(2)刚性块体单元平衡方程
设整个节理岩质边坡被离散为n个刚性块体单元，
Figure FDA0002444805210000021
为任意刚性块体单元i上的结构面数量，对于边坡中的任意刚性块体单元i，其受各种力的作用而保持平衡状态，同时考虑力的平衡和力矩的平衡，对于单一刚性块体单元i其力和力矩平衡方程为：
Figure FDA0002444805210000022
上式中：n是边坡中块体的数量；
Figure FDA0002444805210000023
是刚性块体单元块体i上的结构面数量；
Figure FDA0002444805210000024
是结构面k形心Pk上作用的内力向量；
Figure FDA0002444805210000025
是刚性块体单元i形心Ci作用的等效荷载力向量；Tk是结构面k的局部坐标(nk,sk)到总体坐标系(x,y)的转换矩阵，
Figure FDA0002444805210000026
αk是nk轴与x轴的夹角，逆时针为正；δ1是结构面k形心Dk到点Pk的距离；δ2是岩块i中Dk点到形心Ci到的距离；K1为超载系数；
对于含有n个块体的二维岩质边坡，平衡方程采用矩阵、向量的形式简写如下：
Figure FDA0002444805210000027
上式中：C是考虑平动和转动效应的整体平衡矩阵，
Figure FDA0002444805210000028
是所有结构面的内力向量，
Figure FDA0002444805210000029
是所有块体单元形心的等效荷载力向量；
(3)结构面的屈服条件
假设节理岩质边坡的破坏仅仅发生在结构面上，在同时考虑岩块的平动和转动力学效应时，任意相邻块体单元之间的结构面会发生三种破坏模式，包括剪切破坏，拉伸破坏和转动破坏；
对于任意一条结构面k,在局部坐标系(nk,sk)中其形心上作用有内力向量
Figure FDA00024448052100000210
则结构面的剪切破坏Mohr-Coulomb屈服条件为：
Figure FDA00024448052100000211
上式中：ns是整个边坡中结构面的数量，lk是结构面k的长度，c是结构面的凝聚力，
Figure FDA00024448052100000212
是结构面的内摩擦角，Nk是结构面k的法向力，规定以拉力为正，Vk是结构面k的剪力；
结构面的抗拉屈服条件为：
Nk-σTlk≤0,k＝1,…,ns
上式中：ns为整个边坡中结构面的数量，lk为结构面k的长度，Nk为结构面k的法向力，规定以拉力为正，σT为结构面的抗拉强度，由于贯通的节理一般不能承受拉应力，因此在结构面的抗拉强度取σT＝0；
当相邻块体单元的结构面发生转动破坏时，块体单元发生转动而处于极限状态，由于内力弯矩Mk的方向有顺时针和逆时针两种方向，因此块体单元分别绕结构面的两个端点P1、P2发生两种转动破坏，结构面的转动屈服条件写为：
Figure FDA0002444805210000031
上式中：ns是整个边坡中结构面的数量，lk是结构面k的长度,Nk是结构面k的法向力，规定以拉力为正,Mk是结构面k的内力弯矩；
综合式以上各式，得到考虑块体单元平动和转动力学效应的广义屈服条件：
Figure FDA0002444805210000032
上式中：k＝(1,…,ns)，
Figure FDA0002444805210000033
是以结构面内力向量
Figure FDA0002444805210000034
表示的结构面k的广义屈服条件，简写如下：
Figure FDA0002444805210000035
f为广义屈服条件，
Figure FDA0002444805210000036
为所有的结构面内力向量；
(4)静力边界条件
根据下限定理，静力许可应力场必须满足已知的静力边界条件，考虑节理岩质边坡中已知的边界条件为
Figure FDA0002444805210000037
的界面b，其边界条件表达式为：
Figure FDA0002444805210000038
Figure FDA0002444805210000039
是已知边界条件的结构面的内力向量，
Figure FDA00024448052100000310
是已知边界条件结构面上作用的力和力矩；
(5)同时考虑岩体平动和转动效应的节理岩质边坡承载力分析下限法数学规划模型
根据下限定理，结合目标函数、平衡方程、广义屈服条件和边界条件，得到同时考虑岩体平动和转动效应的节理岩质边坡承载力分析下限法线性数学规划模型：
Figure FDA00024448052100000311
上式为一个线性数学规划模型，通过优化算法求解得到节理岩质边坡的超载系数和对应的屈服区；
(四)优化算法求解节理岩质边坡的极限承载力
采用内点算法对线性数学规划模型进行求解，获得节理岩质边坡稳定性的超载系数的下限解以及对应屈服区。
2.根据权利要求1所述的考虑转动效应的岩质边坡极限承载力分析下限法 ，其特征在于：节理岩质边坡的强度储备系数为：
Figure FDA0002444805210000041
上式中：c,
Figure FDA0002444805210000042
为结构面的凝聚力和摩擦角，c',
Figure FDA0002444805210000043
为结构面进行强度折减以后的凝聚力和摩擦角，σT,σT'为岩体折减前、后的抗拉强度；采用迭代法求解强度储备系数K2。
CN201710381897.0A 2017-05-26 2017-05-26 一种考虑转动效应的岩质边坡极限承载力分析下限法 Active CN107330145B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710381897.0A CN107330145B (zh) 2017-05-26 2017-05-26 一种考虑转动效应的岩质边坡极限承载力分析下限法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710381897.0A CN107330145B (zh) 2017-05-26 2017-05-26 一种考虑转动效应的岩质边坡极限承载力分析下限法

Publications (2)

Publication Number Publication Date
CN107330145A CN107330145A (zh) 2017-11-07
CN107330145B true CN107330145B (zh) 2020-07-31

Family

ID=60192885

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710381897.0A Active CN107330145B (zh) 2017-05-26 2017-05-26 一种考虑转动效应的岩质边坡极限承载力分析下限法

Country Status (1)

Country Link
CN (1) CN107330145B (zh)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108170898B (zh) * 2017-12-06 2021-12-17 昆明理工大学 一种节理岩质边坡可靠度分析下限法
CN108446251B (zh) * 2018-01-17 2021-08-20 昆明理工大学 一种岩质边坡主动岩石压力的计算方法
CN108535121B (zh) * 2018-03-07 2020-10-23 华能澜沧江水电股份有限公司 一种新的岩石统计损伤本构模型的构建方法
CN108763697B (zh) * 2018-05-18 2022-05-20 昆明理工大学 一种岩质边坡中岩块失效概率的计算方法
CN109858128B (zh) * 2019-01-23 2020-12-22 中国水利水电科学研究院 基于块分系统的边坡稳定性二维极限分析下限解法
CN109902345B (zh) * 2019-01-23 2021-02-12 中国水利水电科学研究院 基于极限分析下限定理的二维边坡极限荷载求解方法
CN109871656B (zh) * 2019-03-20 2023-04-25 东北大学 一种基于动态强度折减dda法的边坡稳定性分析系统
CN110083883B (zh) * 2019-04-04 2021-06-01 中国水利水电科学研究院 基于极限分析下限定理的三维边坡安全系数求解方法
CN110866299B (zh) * 2019-11-14 2021-08-17 中国水利水电科学研究院 基于有限元直接内力法的加高重力坝混凝土分载计算方法
CN110889177B (zh) * 2019-11-15 2023-05-23 宁波大学 基于多目标模糊优选动态规划的系列尺寸岩体锚固结构面剪切试验的多油缸布置方法
CN111611728B (zh) * 2020-06-29 2023-03-28 中国有色金属工业昆明勘察设计研究院有限公司 昔格达地层岩质边坡后缘滑裂面位置的稳定性计算方法
CN111651902B (zh) * 2020-06-29 2023-03-31 中国有色金属工业昆明勘察设计研究院有限公司 一种昔格达地层岩质边坡稳定性计算的拟动力下限法
CN112069574A (zh) * 2020-08-26 2020-12-11 盐城工学院 基于极限分析上限定理的圆形基坑稳定性分析方法
CN112487675A (zh) * 2020-11-12 2021-03-12 河海大学 一种基于非协调模型的边坡加固措施评价方法
CN112541257B (zh) * 2020-12-08 2023-02-10 中国有色金属工业昆明勘察设计研究院有限公司 一种基于溶蚀速率的灰岩边坡服役寿命的计算方法
CN113899621B (zh) * 2021-08-27 2023-12-29 北京工业大学 一种粗糙节理可换和节理倾角可变的边坡稳定性试验系统
CN113742834B (zh) * 2021-09-18 2023-09-15 辽宁工程技术大学 一种中间桥对边坡支挡效应的二维分析方法
CN114065590B (zh) * 2021-11-25 2023-04-07 中国电建集团成都勘测设计研究院有限公司 一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法
CN114491737B (zh) * 2021-12-30 2023-07-14 北京市政路桥股份有限公司 一种路堑砌体挡墙的稳定性力学分析方法
CN114417482A (zh) * 2022-02-16 2022-04-29 中信建筑设计研究总院有限公司 一种加固后砌体填充墙的简化力学模型的建立方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106126892A (zh) * 2016-06-20 2016-11-16 昆明理工大学 一种浆砌石挡土墙极限承载力的塑性极限分析上限法
CN106126796A (zh) * 2016-06-20 2016-11-16 昆明理工大学 一种带岩桥的岩质边坡极限承载力的塑性极限分析下限法
CN106547942A (zh) * 2016-09-26 2017-03-29 昆明理工大学 一种外倾顺层节理岩体被动侧向岩石压力的计算方法
CN106557608A (zh) * 2016-09-26 2017-04-05 昆明理工大学 一种基于混合数值离散的非贯通节理岩体的塑性极限分析上限法
CN106682330A (zh) * 2016-12-30 2017-05-17 湖南科技大学 一种深埋硐室围岩压力的计算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8224631B2 (en) * 2008-08-18 2012-07-17 Fci Holdings Delaware, Inc. Stress, geologic, and support analysis methodology for underground openings

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106126892A (zh) * 2016-06-20 2016-11-16 昆明理工大学 一种浆砌石挡土墙极限承载力的塑性极限分析上限法
CN106126796A (zh) * 2016-06-20 2016-11-16 昆明理工大学 一种带岩桥的岩质边坡极限承载力的塑性极限分析下限法
CN106547942A (zh) * 2016-09-26 2017-03-29 昆明理工大学 一种外倾顺层节理岩体被动侧向岩石压力的计算方法
CN106557608A (zh) * 2016-09-26 2017-04-05 昆明理工大学 一种基于混合数值离散的非贯通节理岩体的塑性极限分析上限法
CN106682330A (zh) * 2016-12-30 2017-05-17 湖南科技大学 一种深埋硐室围岩压力的计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于刚性块体系统的岩质边坡稳定性下限法研究;李泽 等;《岩土工程学报》;20120831;第1534-1540页 *

Also Published As

Publication number Publication date
CN107330145A (zh) 2017-11-07

Similar Documents

Publication Publication Date Title
CN107330145B (zh) 一种考虑转动效应的岩质边坡极限承载力分析下限法
CN107341285B (zh) 一种地震力作用下节理岩质边坡倾倒破坏极限承载力分析下限法
CN107330146B (zh) 一种考虑转动效应的岩质边坡极限承载力分析上限法
He et al. An innovative approach for gob-side entry retaining in thick coal seam longwall mining
Pickhaver et al. An equivalent beam method to model masonry buildings in 3D finite element analysis
Han et al. Ultimate bearing capacity of the welded hollow spherical joints in spatial reticulated structures
CN106557608A (zh) 一种基于混合数值离散的非贯通节理岩体的塑性极限分析上限法
Ding et al. Face stability analysis of shallow circular tunnels driven by a pressurized shield in purely cohesive soils under undrained conditions
CN105404758B (zh) 一种基于有限单元法的固体连续介质变形的数值模拟方法
Chen et al. Global concurrent cross-scale nonlinear analysis approach of complex CFRD systems considering dynamic impervious panel-rockfill material-foundation interactions
CN104281730B (zh) 一种大转动变形的板壳结构动响应的有限元分析方法
CN106126892B (zh) 一种浆砌石挡土墙极限承载力的塑性极限分析上限法
CN108170898B (zh) 一种节理岩质边坡可靠度分析下限法
Vo-Minh et al. A stable node-based smoothed finite element method for stability analysis of two circular tunnels at different depths in cohesive-frictional soils
Gong et al. Modelling rock failure with a novel continuous to discontinuous method
CN107992687A (zh) 基于蒙特卡洛方法的岩质边坡可靠度塑性极限分析上限法
Xiong et al. Instability control of roadway surrounding rock in close-distance coal seam groups under repeated mining
Wang et al. Optimization and experimental research on treelike joints based on generative design and powder bed fusion
CN108446251B (zh) 一种岩质边坡主动岩石压力的计算方法
Korolev et al. Basic stages of creating a BIM model for transport infrastructure objects
Lan et al. Review of the double-row pile supporting structure and its force and deformation characteristics
Liu et al. A method of backfill mining crossing the interchange bridge and application of a ground subsidence prediction model
Nguyen-Xuan et al. An adaptive strategy based on conforming quadtree meshes for kinematic limit analysis
Kang et al. Development of three-dimensional numerical manifold method with cover-based contact theory
Chen et al. Study on the Deformation Mechanisms of the Surrounding Rock and Its Supporting Technology for Large Section Whole Coal Cavern Groups

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Li Ze

Inventor after: Peng Pu

Inventor after: Zhou Yu

Inventor after: Zhang Xiaoyan

Inventor before: Li Ze

Inventor before: Zhou Yu

Inventor before: Zhang Xiaoyan

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant