CN107330145A - 一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法 - Google Patents
一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法 Download PDFInfo
- Publication number
- CN107330145A CN107330145A CN201710381897.0A CN201710381897A CN107330145A CN 107330145 A CN107330145 A CN 107330145A CN 201710381897 A CN201710381897 A CN 201710381897A CN 107330145 A CN107330145 A CN 107330145A
- Authority
- CN
- China
- Prior art keywords
- structural plane
- slope
- msub
- mrow
- mover
- 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.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 138
- 230000000694 effects Effects 0.000 title claims abstract description 59
- 238000013519 translation Methods 0.000 title claims abstract description 53
- 238000004458 analytical method Methods 0.000 title claims abstract description 38
- 230000033001 locomotion Effects 0.000 claims abstract description 48
- 230000006378 damage Effects 0.000 claims abstract description 38
- 230000003068 static effect Effects 0.000 claims abstract description 20
- 238000010008 shearing Methods 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims description 32
- 239000000463 material Substances 0.000 claims description 14
- 238000011160 research Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 230000001413 cellular effect Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000003780 insertion Methods 0.000 claims description 3
- 230000037431 insertion Effects 0.000 claims description 3
- 230000015556 catabolic process Effects 0.000 claims description 2
- 239000002131 composite material Substances 0.000 claims description 2
- 238000006731 degradation reaction Methods 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 13
- 238000005516 engineering process Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 238000005452 bending Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011439 discrete element method Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Civil Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Structural Engineering (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Paleontology (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Soil Sciences (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
- Devices Affording Protection Of Roads Or Walls For Sound Insulation (AREA)
Abstract
本发明涉及一种同时考虑岩体平动和转动效应的承载力分析下限法，属于岩质边坡稳定性分析技术领域。本发明采用刚性块体单元将其离散为刚性岩块+结构面的几何系统，以结构面的法向力、剪力、弯矩为未知量，构建同时满足刚性块体单元考虑平动、转动效应的平衡方程、结构面剪切屈服条件、结构面拉伸屈服条件和结构面转动屈服条件以及静力边界条件的静力许可应力场，以超载系数为目标函数，建立下限法线性数学规划模型；采用内点算法对线性数学规划模型求解，同时采用迭代算法求解强度储备系数，获得节理岩质边坡的超载系数或强度储备系数的下限解和对应的屈服区。本发明概念明确、计算精度高，可应用于节理岩质边坡平动或转动破坏的承载力分析。
Description
技术领域
本发明是一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法，特别涉及一种同时考虑岩体结构面剪切破坏、拉伸破坏和转动破坏的承载力分析法。属于岩质边坡稳定性分析技术领域。
背景技术
节理岩质边坡的失稳是一种常见的地质灾害。在国内外公路、铁路、矿山和水利水电等工程中均碰到大量节理岩质边坡失稳破坏的灾难性事故，如国内天生桥二级水电站边坡倾倒破坏、金川露天矿边坡失稳破坏等等。节理岩质边坡的承载力和稳定性评价方法已成为了岩土工程界必须解决的重要技术问题。
节理岩质边坡由于内部存在着大量的节理、裂隙，使得其力学特性非常复杂。节理岩质边坡一般认为是由岩块和结构面组成的复杂系统，岩块的运动方式一般可分为平动和转动，结构面的破坏模式也呈现多样化。岩块平动滑移破坏问题相对简单，岩块转动倾倒破坏问题则比较复杂，一般岩块的转动倾倒破坏可分为块状倾倒、弯曲倾倒和块状弯曲倾倒等类型。自20世纪50年代以来，岩体平动破坏和转动破坏的研究就已经被重视。为了解决岩体平动、转动倾倒的稳定性分析问题，基于极限平衡原理的边坡稳定计算方法被提了出来，并且被广泛应用于工程中。这种方法一般对岩块变形失稳模式做一些假定，把岩质边坡的滑体分为稳定区、倾倒区或滑动区，但对一些特殊问题其合理性有待验证。其后，很多数值方法被提出来用于岩体平动、转动破坏失稳的分析计算，比如有限单元法、离散单元法、块体单元法、DDA等。同时，岩质边坡平动、转动的破坏物理模型试验也被广泛应用于工程实践。
就目前的研究成果来说，同时考虑平动、转动效应的节理岩质边坡的承载力分析领域还存在一些有待于进一步探讨的问题。具体不足表现在：
(1)很多分析方法不能同时考虑岩体的平动和转动效应，比如刚体极限平衡法、有限元法等一般只能考虑岩体的平动效应，但是工程实践表明岩体的转动效应是非常重要、不可或缺的，比如：对于顺层节理岩质边坡而言转动破坏是其主要破坏模式。
(2)节理岩质边坡中的结构面破坏的模式也非常多样，一般认为包括三种：滑移剪切破坏、拉伸破坏和转动破坏，但现有的分析方法很难同时考虑这三种破坏模式。
(3)在极限分析领域，特别是下限法分析领域，研究岩体同时考虑平动、转动效应的研究成果较少。
本发明基于国家自然科学基金项目(资助编号：51564026)的研究工作，将塑性极限分析下限法理论、刚性块体单元离散技术以及数学规划手段结合起来，提出一种同时考虑平动、转动效应的岩质边坡极限承载力分析的下限法。
发明内容
本发明的目的是提供一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法，为节理岩质边坡稳定性分析提供一种新方法。
本发明的基本原理是：以节理岩质边坡为研究对象，采用刚性块体单元将其离散为刚性岩块+结构面的几何系统，以结构面的法向力、剪力、弯矩为未知量，基于塑性极限分析下限法理论，构建同时满足：刚性块体单元考虑平动、转动效应的平衡方程、结构面剪切屈服条件、结构面拉伸屈服条件和结构面转动屈服条件以及静力边界条件的静力许可应力场，然后以荷载的超载系数为目标函数，建立同时考虑平动和转动效应的节理岩质边坡承载力分析的下限法线性数学规划模型；然后采用内点算法对线性数学规划模型进行求解，同时可采用迭代算法求解强度储备系数，最终可获得节理岩质边坡的超载系数或强度储备系数的下限解和对应的屈服区。
本发明的技术路线如图1所示。
本发明的同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法的技术方案依次按以下步骤进行：
一、拟定节理岩质边坡的计算参数
根据节理岩质边坡的实际情况，拟定下限法分析需要的计算参数，主要包括：岩质边坡的几何参数、节理参数、岩体和节理材料的参数(容重、凝聚力、摩擦角、抗拉强度等)、边坡荷载参数信息。
二、采用刚性块体单元离散节理岩质边坡
本发明采用刚性块体单元离散节理岩质边坡，刚性块体单元受力如图2所示，相邻单元的结构面受力如图3所示。其中(x,y)为总体坐标系；(nk,sk)为相邻块体单元i、j之间结构面k的局部坐标系，nk为结构面的外法线方向，sk为结构面的切向；刚性块体单元i形心Ci上作用有等效荷载力向量fxi为沿x方向的外力，fyi为沿y方向的外力，mi为外力弯矩(规定逆时针为正向)；相邻块体单元i、j之间的结构面k形心Pk上作用有内力向量另一侧作用有其反作用力向量Nk是沿结构面外法线nk方向的法向力(拉正压负)，Vk是沿结构面切向sk方向的剪力(使块体产生逆时针转动为正)，Mk为结构面的内力弯矩(规定逆时针为正向)。
采用刚性块体单元离散节理岩质边坡以后，边坡变成刚性块体单元+结构面的几何系统，本发明以刚性块体单元结构面的形心Pk上作用的内力向量为未知量构建岩质边坡的静力许可应力场。为了简化计算，本发明作如下假设：(1)假设岩块为刚体，岩质边坡的破坏只会发生在相邻块体单元之间的结构面上；(2)假设结构面材料为理想刚塑性材料，(3)假设结构面有可能发生剪切破坏、拉伸破坏或转动破坏，(4)岩块同时发生平动或转动。
三、建立同时考虑岩体平动、转动效应的节理岩质边坡承载力分析下限法数学规划模型
本发明的目的是求解节理岩质边坡超载系数或强度储备系数的下限解，需要建立静力许可应力场。根据下限定理，要建立边坡的静力许可应力场需满足平衡方程、屈服条件和边界条件。
1、目标函数
本发明方法将节理岩质边坡的超载定义为：
K1＝Fc/Fa (1)
上式中，Fc是极限荷载，Fa是边坡结构当前实际施加的外荷载，对于岩质边坡，求极限状态时其超载的外荷载可以是自重也可以是面力或集中力。
根据下限定理，岩质边坡在达到极限状态时，需求超载系数的最大值。本发明将超载系数设为目标函数如下：
Maximize:K1 (2)
2、刚性块体单元平衡方程
设整个节理岩质边坡被离散为n个刚性块体单元，为任意刚性块体单元i上的结构面数量。对于图2所示边坡中的任意刚性块体单元i，其受各种力的作用而保持平衡状态，本发明同时考虑块体单元的平动、转动效应，因此同时考虑力的平衡和力矩的平衡。对于单一刚性块体单元i其力和力矩平衡方程为：
上式中：n是边坡中块体的数量；是刚性块体单元块体i上的结构面数量；是结构面k形心Pi上作用的内力向量；是刚性块体单元i形心Ci作用的等效荷载力向量；Tk是结构面k的局部坐标(nk,sk)到总体坐标系(x,y)的转换矩阵，αk是nk轴与x轴的夹角，逆时针为正；δ1是图2中结构面k形心Dk到点Pk的距离；δ2是图2中岩块i中Dk点到形心Ci到的距离；K1为超载系数；K1为超载系数。
对于含有n个块体的二维岩质边坡，式(3)包括2n个力的平衡方程和n个力矩平衡方程，可采用矩阵、向量的形式简写如下：
上式中：C是考虑平动、转动效应的整体平衡矩阵，是所有结构面的内力向量，是所有块体单元形心的等效荷载力向量。
3、结构面的屈服条件
本发明假设节理岩质边坡的破坏仅仅发生在结构面上。在同时考虑岩块的平动、转动力学效应时，任意相邻块体单元之间的结构面会发生三种破坏模式(如图4所示)，包括剪切破坏，拉伸破坏和转动破坏。
对于任意一条结构面k(如图3所示),在局部坐标系(nk,sk)中其形心上作用有内力向量则结构面的剪切破坏Mohr-Coulomb屈服条件为：
上式中：ns是整个边坡中结构面的数量，lk是结构面k的长度，c是结构面的凝聚力，是结构面的内摩擦角，Nk是结构面k的法向力(规定以拉力为正)，Vk是结构面k的剪力。
结构面的抗拉屈服条件为：
Nk-σTlk≤0,k＝(1,…,ns) (6)
上式中：ns为整个边坡中结构面的数量，lk为结构面k的长度，Nk为结构面k的法向力(规定以拉力为正)，σT为结构面的抗拉强度，由于贯通的节理一般不能承受拉应力，因此在本发明中结构面的抗拉强度取σT＝0。
结合结构面的剪切屈服条件式(5)和抗拉屈服条件式(6)，则可得到结构面的拉剪破坏的修正Mohr-Coulomb屈服条件如图5所示。
当相邻块体单元的结构面发生转动破坏时，块体单元发生转动极限状态下的受力如图6所示，由于内力弯矩Mk的方向有两种可能的方向(顺时针和逆时针)，因此块体单元有可能分别绕结构面的两个端点P1、P2发生两种转动破坏。结构面的转动屈服条件如图7所示，可写为：
上式中：ns是整个边坡中结构面的数量，lk是结构面k的长度,Nk是结构面k的法向力(规定以拉力为正),Mk是结构面k的内力弯矩。
综合式(5)、式(6)、式(7)，得到考虑块体单元平动、转动力学效应的广义屈服条件：
上式中：k＝(1,…,ns)，为以结构面内力向量表示的结构面k的广义屈服条件。
式(9)共计为5ns个约束条件，5ns个约束均为线性约束条件。将广义屈服条件式(8)简写如下：
上式中：f为广义屈服条件，为所有的结构面内力向量。
4、静力边界条件
根据下限定理，静力许可应力场必须满足已知的静力边界条件。考虑节理岩质边坡中已知的边界条件为的界面b，其边界条件表达式为：
上式中：是已知边界条件的结构面的内力向量，是已知边界条件结构面上作用的力和力矩。
5、同时考虑岩体平动、转动效应的节理岩质边坡承载力分析下限法数学规划模型
根据下限定理，结合目标函数式(2)、平衡方程式(4)、广义屈服条件式(10)、边界条件式(11)，则可得到同时考虑岩体平动、转动效应的节理岩质边坡承载力分析下限法线性数学规划模型：
上式为一个线性数学规划模型，通过优化算法求解可到节理岩质边坡的超载系数和对应的屈服区。
四、优化算法求解节理岩质边坡的极限承载力
通过上述步骤建立了同时考虑岩体平动、转动效应的节理岩质边坡极限承载力分析的下限法数学规划模型(式(11))，其是一个大规模线性数学规划模型。本发明采用Fortran 12.0编写下限法计算程序，计算程序包括前处理、构建下限法数学规划模型、优化求解、后处理等4个部分。本发明采用内点算法对线性数学规划模型进行求解。最终可获得节理岩质边坡稳定性的超载系数的下限解以及对应屈服区。
对于节理岩质边坡，有时候除了超载系数还需要求解强度储备系数，本发明的将强度储备系数定义为：
上式中：c,为结构面的凝聚力和摩擦角，c',为结构面进行强度折减以后的凝聚力和摩擦角，σ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-3) | 黏聚力/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节理岩质边坡极限荷载计算结果
实施例2：实施例2为使用式(12)求解一个含有6个岩块的节理岩质边坡的强度储备系数。
(一)拟定节理岩质边坡的计算参数。
图11为实施例2节理岩质边坡的几何形状示意图。此边坡底面倾角25°，边坡包含6个岩块，每个岩块的几何尺寸详见图11。边坡材料的物理力学参数详见表3。6个岩块中，岩块1为固定岩块、不能移动，岩块2～岩块6可能发生平动或转动。本实施例的目的是求解岩质边坡的强度储备系数。
表3实施例2节理岩质边坡材料的物理力学参数表
| 材料名称 | 重度/(kN·m-3) | 黏聚力/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 (3)
1.一种同时考虑平动和转动效应的节理岩质边坡极限承载力分析下限法，其特征在于：以节理岩质边坡为研究对象，采用刚性块体单元将其离散为刚性岩块+结构面的几何系统，以结构面的法向力、剪力、弯矩为未知量，基于塑性极限分析下限法理论，构建同时满足刚性块体单元考虑平动、转动效应的平衡方程、结构面剪切屈服条件、结构面拉伸屈服条件和结构面转动屈服条件以及静力边界条件的静力许可应力场，然后以荷载的超载系数为目标函数，建立同时考虑平动和转动效应的节理岩质边坡承载力分析的下限法线性数学规划模型；然后采用内点算法对线性数学规划模型进行求解，同时采用迭代算法求解强度储备系数，最终获得节理岩质边坡的超载系数或强度储备系数的下限解和对应的屈服区。
2.根据权利要求1所述的同时考虑平动和转动效应的节理岩质边坡极限承载力分析下限法，其特征在于：具体步骤如下：
(一)拟定节理岩质边坡的计算参数
根据节理岩质边坡的实际情况，拟定下限法分析需要的计算参数，包括：岩质边坡的几何参数、节理参数、边坡荷载参数信息及岩体和节理材料的参数，所述岩体和节理材料的参数容重、凝聚力、摩擦角、抗拉强度；
(二)采用刚性块体单元离散节理岩质边坡
采用刚性块体单元离散节理岩质边坡，，边坡变成刚性块体单元+结构面的几何系统，(x,y)为总体坐标系；(nk,sk)为相邻块体单元i、j之间结构面k的局部坐标系，nk为结构面的外法线方向，sk为结构面的切向；刚性块体单元i形心Ci上作用有等效荷载力向量fxi为沿x方向的外力，fyi为沿y方向的外力，mi为外力弯矩，规定逆时针为正向；相邻块体单元i、j之间的结构面k形心Pi上作用有内力向量另一侧作用有其反作用力向量Nk是沿结构面外法线nk方向的法向力，拉正压负，Vk是沿结构面切向sk方向的剪力，使块体产生逆时针转动为正，Mk为结构面的内力弯矩，规定逆时针为正向；
以刚性块体单元结构面的形心Pk上作用的内力向量为未知量构建岩质边坡的静力许可应力场；为了简化计算，作如下假设：(1)假设岩块为刚体，岩质边坡的破坏只会发生在相邻块体单元之间的结构面上(2)假设结构面材料为理想刚塑性材料(3)假设结构面有可能发生剪切破坏、拉伸破坏或转动破坏(4)岩块同时发生平动或转动；
(三)建立同时考虑岩体平动和转动效应的节理岩质边坡承载力分析下限法数学规划模型
根据下限定理，要建立边坡的静力许可应力场需满足平衡方程、屈服条件和边界条件；
(1)目标函数
将超载系数设为目标函数如下：Maximize:K1，节理岩质边坡的超载系数定义为：K1＝Fc/Fa；
上式中，Fc是极限荷载，Fa是边坡结构当前实际施加的外荷载，对于岩质边坡，求极限状态时其超载的外荷载为自重、面力或集中力；
(2)刚性块体单元平衡方程
设整个节理岩质边坡被离散为n个刚性块体单元，为任意刚性块体单元i上的结构面数量，对于边坡中的任意刚性块体单元i，其受各种力的作用而保持平衡状态，同时考虑力的平衡和力矩的平衡，对于单一刚性块体单元i其力和力矩平衡方程为：
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msubsup>
<mi>m</mi>
<mi>i</mi>
<mi>k</mi>
</msubsup>
</munderover>
<msub>
<mover>
<mi>T</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>Q</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mover>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
上式中：n是边坡中块体的数量；是刚性块体单元块体i上的结构面数量；是结构面k形心Pk上作用的内力向量；是刚性块体单元i形心Ci作用的等效荷载力向量；Tk是结构面k的局部坐标(nk,sk)到总体坐标系(x,y)的转换矩阵，αk是nk轴与x轴的夹角，逆时针为正；δ1是图2中结构面k形心Dk到点Pk的距离；δ2是图2中岩块i中Dk点到形心Ci到的距离；K1为超载系数；
对于含有n个块体的二维岩质边坡，平衡方程采用矩阵、向量的形式简写如下：
<mrow>
<mi>C</mi>
<mover>
<mi>Q</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mover>
<mi>F</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>=</mo>
<mn>0</mn>
</mrow>
上式中：C是考虑平动和转动效应的整体平衡矩阵，是所有结构面的内力向量，是所有块体单元形心的等效荷载力向量；
(3)结构面的屈服条件
假设节理岩质边坡的破坏仅仅发生在结构面上，在同时考虑岩块的平动和转动力学效应时，任意相邻块体单元之间的结构面会发生三种破坏模式，包括剪切破坏，拉伸破坏和转动破坏；
对于任意一条结构面k,在局部坐标系(nk,sk)中其形心上作用有内力向量则结构面的剪切破坏Mohr-Coulomb屈服条件为：
上式中：ns是整个边坡中结构面的数量，lk是结构面k的长度，c是结构面的凝聚力，是结构面的内摩擦角，Nk是结构面k的法向力，规定以拉力为正，Vk是结构面k的剪力；
结构面的抗拉屈服条件为：
Nk-σTlk≤0,k＝(1,…,ns)
上式中：ns为整个边坡中结构面的数量，lk为结构面k的长度，Nk为结构面k的法向力，规定以拉力为正，σT为结构面的抗拉强度，由于贯通的节理一般不能承受拉应力，因此在结构面的抗拉强度取σT＝0；
当相邻块体单元的结构面发生转动破坏时，块体单元发生转动而处于极限状态，由于内力弯矩Mk的方向有顺时针和逆时针两种方向，因此块体单元分别绕结构面的两个端点P1、P2发生两种转动破坏，结构面的转动屈服条件写为：
<mrow>
<mo>|</mo>
<msub>
<mi>M</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mo>-</mo>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
<mfrac>
<msub>
<mi>l</mi>
<mi>k</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>&le;</mo>
<mn>0</mn>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>n</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
上式中：ns是整个边坡中结构面的数量，lk是结构面k的长度,Nk是结构面k的法向力，规定以拉力为正,Mk是结构面k的内力弯矩；
综合式以上各式，得到考虑块体单元平动和转动力学效应的广义屈服条件：
上式中：k＝(1,…,ns)，是以结构面内力向量表示的结构面k的广义屈服条件，简写如下：f为广义屈服条件，为所有的结构面内力向量；
(4)静力边界条件
根据下限定理，静力许可应力场必须满足已知的静力边界条件，考虑节理岩质边坡中已知的边界条件为的界面b，其边界条件表达式为： 是已知边界条件的结构面的内力向量，是已知边界条件结构面上作用的力和力矩；
(5)同时考虑岩体平动和转动效应的节理岩质边坡承载力分析下限法数学规划模型
根据下限定理，结合目标函数、平衡方程、广义屈服条件和边界条件，得到同时考虑岩体平动和转动效应的节理岩质边坡承载力分析下限法线性数学规划模型：
<mfenced open = "" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>M</mi>
<mi>a</mi>
<mi>x</mi>
<mi>i</mi>
<mi>m</mi>
<mi>i</mi>
<mi>z</mi>
<mi>e</mi>
<mo>:</mo>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>S</mi>
<mi>u</mi>
<mi>b</mi>
<mi>j</mi>
<mi>e</mi>
<mi>c</mi>
<mi>t</mi>
<mi> </mi>
<mi>t</mi>
<mi>o</mi>
<mo>:</mo>
<mi>C</mi>
<mover>
<mi>Q</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mover>
<mi>F</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>Q</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>Q</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>b</mi>
</msub>
<mo>=</mo>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
上式为一个线性数学规划模型，通过优化算法求解得到节理岩质边坡的超载系数和对应的屈服区；
(四)优化算法求解节理岩质边坡的极限承载力
采用内点算法对线性数学规划模型进行求解，获得节理岩质边坡稳定性的超载系数的下限解以及对应屈服区。
3.根据权利要求1所述的同时考虑平动和转动效应的岩质边坡承载力分析下限法，其特征在于：节理岩质边坡的强度储备系数为：
上式中：c,为结构面的凝聚力和摩擦角，c',为结构面进行强度折减以后的凝聚力和摩擦角，σT,σT'为岩体折减前、后的抗拉强度；采用迭代法求解强度储备系数K2。
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 true CN107330145A (zh) | 2017-11-07 |
| CN107330145B 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) |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108170898A (zh) * | 2017-12-06 | 2018-06-15 | 昆明理工大学 | 一种节理岩质边坡可靠度分析下限法 |
| CN108446251A (zh) * | 2018-01-17 | 2018-08-24 | 昆明理工大学 | 一种岩质边坡沿陡倾的外倾软弱结构面破坏时的主动岩石压力的计算方法 |
| CN108535121A (zh) * | 2018-03-07 | 2018-09-14 | 华能澜沧江水电股份有限公司 | 一种新的岩石统计损伤本构模型的构建方法 |
| CN108763697A (zh) * | 2018-05-18 | 2018-11-06 | 昆明理工大学 | 一种岩质边坡中岩块失效概率的计算方法 |
| CN109858128A (zh) * | 2019-01-23 | 2019-06-07 | 中国水利水电科学研究院 | 基于块分系统的边坡稳定性二维极限分析下限解法 |
| CN109871656A (zh) * | 2019-03-20 | 2019-06-11 | 东北大学 | 一种基于动态强度折减dda法的边坡稳定性分析系统 |
| CN109902345A (zh) * | 2019-01-23 | 2019-06-18 | 中国水利水电科学研究院 | 基于极限分析下限定理的二维边坡极限荷载求解方法 |
| CN110083883A (zh) * | 2019-04-04 | 2019-08-02 | 中国水利水电科学研究院 | 基于极限分析下限定理的三维边坡安全系数求解方法 |
| CN110866299A (zh) * | 2019-11-14 | 2020-03-06 | 中国水利水电科学研究院 | 基于有限元直接内力法的加高重力坝混凝土分载计算方法 |
| CN110889177A (zh) * | 2019-11-15 | 2020-03-17 | 宁波大学 | 基于多目标模糊优选动态规划的系列尺寸岩体锚固结构面剪切试验的多油缸布置方法 |
| CN111611728A (zh) * | 2020-06-29 | 2020-09-01 | 中国有色金属工业昆明勘察设计研究院有限公司 | 昔格达地层岩质边坡后缘滑裂面位置的稳定性计算方法 |
| CN111651902A (zh) * | 2020-06-29 | 2020-09-11 | 中国有色金属工业昆明勘察设计研究院有限公司 | 一种昔格达地层岩质边坡稳定性计算的拟动力下限法 |
| CN112069574A (zh) * | 2020-08-26 | 2020-12-11 | 盐城工学院 | 基于极限分析上限定理的圆形基坑稳定性分析方法 |
| CN112487675A (zh) * | 2020-11-12 | 2021-03-12 | 河海大学 | 一种基于非协调模型的边坡加固措施评价方法 |
| CN112541257A (zh) * | 2020-12-08 | 2021-03-23 | 中国有色金属工业昆明勘察设计研究院有限公司 | 一种基于溶蚀速率的灰岩边坡服役寿命的计算方法 |
| CN113742834A (zh) * | 2021-09-18 | 2021-12-03 | 辽宁工程技术大学 | 一种中间桥对边坡支挡效应的二维分析方法 |
| CN113899621A (zh) * | 2021-08-27 | 2022-01-07 | 北京工业大学 | 一种粗糙节理可换和节理倾角可变的边坡稳定性试验系统 |
| CN114065590A (zh) * | 2021-11-25 | 2022-02-18 | 中国电建集团成都勘测设计研究院有限公司 | 一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法 |
| CN114417482A (zh) * | 2022-02-16 | 2022-04-29 | 中信建筑设计研究总院有限公司 | 一种加固后砌体填充墙的简化力学模型的建立方法 |
| CN114491737A (zh) * | 2021-12-30 | 2022-05-13 | 北京市政路桥股份有限公司 | 一种路堑砌体挡墙的稳定性力学分析方法 |
| CN114943052A (zh) * | 2022-06-30 | 2022-08-26 | 中国长江三峡集团有限公司 | 一种顺层岩质边坡溃屈失稳的判定方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100042381A1 (en) * | 2008-08-18 | 2010-02-18 | Jennmar Corporation | Stress, Geologic, and Support Analysis Methodology for Underground Openings |
| 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 | 湖南科技大学 | 一种深埋硐室围岩压力的计算方法 |
-
2017
- 2017-05-26 CN CN201710381897.0A patent/CN107330145B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100042381A1 (en) * | 2008-08-18 | 2010-02-18 | Jennmar Corporation | Stress, Geologic, and Support Analysis Methodology for Underground Openings |
| 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)
| Title |
|---|
| 李泽 等: "基于刚性块体系统的岩质边坡稳定性下限法研究", 《岩土工程学报》 * |
Cited By (36)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108170898A (zh) * | 2017-12-06 | 2018-06-15 | 昆明理工大学 | 一种节理岩质边坡可靠度分析下限法 |
| CN108170898B (zh) * | 2017-12-06 | 2021-12-17 | 昆明理工大学 | 一种节理岩质边坡可靠度分析下限法 |
| CN108446251A (zh) * | 2018-01-17 | 2018-08-24 | 昆明理工大学 | 一种岩质边坡沿陡倾的外倾软弱结构面破坏时的主动岩石压力的计算方法 |
| CN108535121A (zh) * | 2018-03-07 | 2018-09-14 | 华能澜沧江水电股份有限公司 | 一种新的岩石统计损伤本构模型的构建方法 |
| CN108535121B (zh) * | 2018-03-07 | 2020-10-23 | 华能澜沧江水电股份有限公司 | 一种新的岩石统计损伤本构模型的构建方法 |
| CN108763697A (zh) * | 2018-05-18 | 2018-11-06 | 昆明理工大学 | 一种岩质边坡中岩块失效概率的计算方法 |
| CN108763697B (zh) * | 2018-05-18 | 2022-05-20 | 昆明理工大学 | 一种岩质边坡中岩块失效概率的计算方法 |
| CN109858128A (zh) * | 2019-01-23 | 2019-06-07 | 中国水利水电科学研究院 | 基于块分系统的边坡稳定性二维极限分析下限解法 |
| CN109902345A (zh) * | 2019-01-23 | 2019-06-18 | 中国水利水电科学研究院 | 基于极限分析下限定理的二维边坡极限荷载求解方法 |
| CN109902345B (zh) * | 2019-01-23 | 2021-02-12 | 中国水利水电科学研究院 | 基于极限分析下限定理的二维边坡极限荷载求解方法 |
| CN109871656B (zh) * | 2019-03-20 | 2023-04-25 | 东北大学 | 一种基于动态强度折减dda法的边坡稳定性分析系统 |
| CN109871656A (zh) * | 2019-03-20 | 2019-06-11 | 东北大学 | 一种基于动态强度折减dda法的边坡稳定性分析系统 |
| CN110083883A (zh) * | 2019-04-04 | 2019-08-02 | 中国水利水电科学研究院 | 基于极限分析下限定理的三维边坡安全系数求解方法 |
| CN110083883B (zh) * | 2019-04-04 | 2021-06-01 | 中国水利水电科学研究院 | 基于极限分析下限定理的三维边坡安全系数求解方法 |
| CN110866299A (zh) * | 2019-11-14 | 2020-03-06 | 中国水利水电科学研究院 | 基于有限元直接内力法的加高重力坝混凝土分载计算方法 |
| CN110866299B (zh) * | 2019-11-14 | 2021-08-17 | 中国水利水电科学研究院 | 基于有限元直接内力法的加高重力坝混凝土分载计算方法 |
| CN110889177A (zh) * | 2019-11-15 | 2020-03-17 | 宁波大学 | 基于多目标模糊优选动态规划的系列尺寸岩体锚固结构面剪切试验的多油缸布置方法 |
| CN110889177B (zh) * | 2019-11-15 | 2023-05-23 | 宁波大学 | 基于多目标模糊优选动态规划的系列尺寸岩体锚固结构面剪切试验的多油缸布置方法 |
| CN111611728A (zh) * | 2020-06-29 | 2020-09-01 | 中国有色金属工业昆明勘察设计研究院有限公司 | 昔格达地层岩质边坡后缘滑裂面位置的稳定性计算方法 |
| CN111651902A (zh) * | 2020-06-29 | 2020-09-11 | 中国有色金属工业昆明勘察设计研究院有限公司 | 一种昔格达地层岩质边坡稳定性计算的拟动力下限法 |
| CN112069574B (zh) * | 2020-08-26 | 2024-05-17 | 盐城工学院 | 基于极限分析上限定理的圆形基坑稳定性分析方法 |
| CN112069574A (zh) * | 2020-08-26 | 2020-12-11 | 盐城工学院 | 基于极限分析上限定理的圆形基坑稳定性分析方法 |
| CN112487675A (zh) * | 2020-11-12 | 2021-03-12 | 河海大学 | 一种基于非协调模型的边坡加固措施评价方法 |
| CN112541257B (zh) * | 2020-12-08 | 2023-02-10 | 中国有色金属工业昆明勘察设计研究院有限公司 | 一种基于溶蚀速率的灰岩边坡服役寿命的计算方法 |
| CN112541257A (zh) * | 2020-12-08 | 2021-03-23 | 中国有色金属工业昆明勘察设计研究院有限公司 | 一种基于溶蚀速率的灰岩边坡服役寿命的计算方法 |
| CN113899621B (zh) * | 2021-08-27 | 2023-12-29 | 北京工业大学 | 一种粗糙节理可换和节理倾角可变的边坡稳定性试验系统 |
| CN113899621A (zh) * | 2021-08-27 | 2022-01-07 | 北京工业大学 | 一种粗糙节理可换和节理倾角可变的边坡稳定性试验系统 |
| CN113742834A (zh) * | 2021-09-18 | 2021-12-03 | 辽宁工程技术大学 | 一种中间桥对边坡支挡效应的二维分析方法 |
| CN113742834B (zh) * | 2021-09-18 | 2023-09-15 | 辽宁工程技术大学 | 一种中间桥对边坡支挡效应的二维分析方法 |
| CN114065590B (zh) * | 2021-11-25 | 2023-04-07 | 中国电建集团成都勘测设计研究院有限公司 | 一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法 |
| CN114065590A (zh) * | 2021-11-25 | 2022-02-18 | 中国电建集团成都勘测设计研究院有限公司 | 一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法 |
| CN114491737B (zh) * | 2021-12-30 | 2023-07-14 | 北京市政路桥股份有限公司 | 一种路堑砌体挡墙的稳定性力学分析方法 |
| CN114491737A (zh) * | 2021-12-30 | 2022-05-13 | 北京市政路桥股份有限公司 | 一种路堑砌体挡墙的稳定性力学分析方法 |
| CN114417482A (zh) * | 2022-02-16 | 2022-04-29 | 中信建筑设计研究总院有限公司 | 一种加固后砌体填充墙的简化力学模型的建立方法 |
| CN114417482B (zh) * | 2022-02-16 | 2024-05-07 | 中信建筑设计研究总院有限公司 | 一种加固后砌体填充墙的简化力学模型的建立方法 |
| CN114943052A (zh) * | 2022-06-30 | 2022-08-26 | 中国长江三峡集团有限公司 | 一种顺层岩质边坡溃屈失稳的判定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107330145B (zh) | 2020-07-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107330145A (zh) | 一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法 | |
| CN107330146A (zh) | 一种同时考虑平动和转动效应的节理岩质边坡极限承载力分析上限法 | |
| CN107341285B (zh) | 一种地震力作用下节理岩质边坡倾倒破坏极限承载力分析下限法 | |
| CN105069191B (zh) | 充分考虑铁塔结构弹性形变的输电铁塔杆件应力计算方法 | |
| CN106570204A (zh) | 一种基于cpu+gpu异构并行计算的透平机械叶片静强度特性分析方法 | |
| CN114638046B (zh) | 一种铁路桥墩数字孪生变截面仿真计算方法 | |
| CN105930616A (zh) | 一种飞机动翼面超静定结构的载荷分配方法 | |
| CN106055764A (zh) | 基于三维壳有限元‑梁模型的风力机叶片位移计算方法 | |
| CN105568864A (zh) | 确定斜拉桥合理施工索力的综合算法 | |
| US20230342515A1 (en) | Reliability-based topology optimization design method for part structure considering bounded hybrid uncertainties | |
| CN110704953B (zh) | 一种大展弦比机翼静气弹性能设计敏度的分析方法 | |
| Zheng et al. | Improved Bell's method for the stability analysis of slopes | |
| CN105205293A (zh) | 用于获得飞机部件气动载荷的方法和系统 | |
| CN110852012A (zh) | 轴压加筋板极限强度预测方法 | |
| CN105404758A (zh) | 一种基于有限单元法的固体连续介质变形的数值模拟方法 | |
| Wu et al. | Nonlinear stability analysis of steel cooling towers considering imperfection sensitivity | |
| CN106202654B (zh) | 一种地下厂房多块体围岩变形的数值模拟方法 | |
| CN104573178B (zh) | 一种整体叶轮强度有限元计算方法 | |
| CN105893699A (zh) | 一种含混合不确定性参数的连杆机构鲁棒综合方法 | |
| CN106909755A (zh) | 一种车辆荷载下沉管隧道管节竖向位移的计算方法 | |
| CN101334803B (zh) | 基于变形修正的动力减缩算法及超单元构造方法 | |
| CN110569606B (zh) | 峡谷高应力区硬岩大型地下洞室轴线布置量化设计方法 | |
| CN107563102A (zh) | 一种承载结构的传力骨架可视化方法 | |
| CN102508978B (zh) | 一种基于群论的对称杆系结构的可动性判定方法 | |
| CN105404740A (zh) | 一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法 |
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 | ||
| 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 |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |