CN114065590B - 一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法 - Google Patents
一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法 Download PDFInfo
- Publication number
- CN114065590B CN114065590B CN202111409221.0A CN202111409221A CN114065590B CN 114065590 B CN114065590 B CN 114065590B CN 202111409221 A CN202111409221 A CN 202111409221A CN 114065590 B CN114065590 B CN 114065590B
- Authority
- CN
- China
- Prior art keywords
- block
- sliding surface
- force
- node
- sliding
- 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
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/20—Hydro energy
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明涉及地下洞室围岩块体稳定分析领域,具体涉及一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法,极大地提高了对地下洞室块体稳定分析的准确性。技术方案包括:建立块体三维几何模型,采用刚体极限平衡法计算块体的安全系数,获得第一安全系数;建立块体三维有限元模型,在块体三维有限元模型上,考虑块体滑面的黏聚力和摩擦系数,以及重力、预应力合支护力、非预应力合支护力、地下水以及地震力,对各滑面上的单元节点的应力,再根据刚体极限平衡法计算块体安全系数,获得第二安全系数;比较第一安全系数与第二安全系数的大小,取较小值作为块体稳定安全系数。本发明适用于对地下洞室围岩块体的稳定性分析。
Description
技术领域
本发明涉及地下洞室围岩块体稳定分析领域,具体涉及一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法。
背景技术
刚体极限平衡法是根据块体的力学平衡原理分析块体在各种破坏模式下的受力状态,以及块体滑面上的抗滑力和下滑力之间的大小关系来评价块体的稳定性。刚体极限平衡法是块体稳定性分析计算的主要方法,也是工程实践中应用最多的一种方法。
刚体极限平衡法虽然在工程中应用的很广泛,但它不能计算失稳前岩体的真实应力情况,更不能得出岩体的变形,即不能反映出岩体内真实的应力-应变关系,只是用此公式来推求岩体稳定性系数,而且对滑面上的应力作了与实际不相符的平均分布假定,所求稳定性系数相当于是滑动面上的一个平均值,而不是精确解。所以,它无法分析岩体从变形到破坏的发展过程,也无法考虑岩体渐进性破坏的具体影响。
发明内容
本发明的目的是提供一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法,极大地提高了对地下洞室块体稳定分析的准确性。
本发明采取如下技术方案实现上述目的,一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法,包括:
步骤1、建立块体三维几何模型,根据块体几何参数、材料参数、力学参数,采用刚体极限平衡法计算块体的安全系数,获得第一安全系数;
步骤2、建立块体三维有限元模型,根据块体的稳定状态和软弱结构面类型,获得块体滑面上节点和单元的应力位移结果;
步骤3、在块体三维有限元模型上,考虑块体滑面的黏聚力和摩擦系数,以及重力、预应力合支护力、非预应力合支护力、地下水以及地震力,对各滑面上的单元节点的应力,计算该节点处的法向力和切向力;
步骤4、根据块体滑面上节点和单元的应力位移结果以及节点处的法向力和切向力计算和显示块体的力系图,再根据刚体极限平衡法计算块体安全系数,获得第二安全系数;
步骤5、比较第一安全系数与第二安全系数的大小,取较小值作为块体稳定安全系数。
进一步的是,步骤2中,根据块体的稳定状态和软弱结构面类型,获得块体滑面上节点和单元的应力位移结果的具体方法包括:
对块体的稳定状态和软弱结构面类型进行判断,若块体处于稳定状态,则无需建立滑面单元,滑面两侧的实体单元边界共用滑面上的节点,进行计算分析,获得滑面上节点和单元的应力位移结果;
若块体处于不稳定状态且软弱夹层厚度小于等于阈值厚度,则不考虑结构面变形,在各滑面建立接触对单元,进行计算分析,获得滑面上节点和单元的应力位移结果;
若块体处于不稳定状态且软弱夹层厚度大于阈值厚度,则需考虑结构面变形,在各滑面建立实体单元,进行计算分析,获得滑面上节点和单元的应力位移结果。
进一步的是,块体滑面上的接触对单元在闭合状态下传递剪力时,考虑了块体软弱结构面类型,并根据软弱结构面的类型,采用不同单元模拟滑面。
进一步的是,块体滑面上的接触对单元在闭合状态下传递剪力时,考虑了块体滑面黏聚力。
进一步的是,对块体滑面上单元的接触状态判断,判断的具体方法包括:在每一个增量步开始之前检查块体结构面上所有接触对相互作用状态,以判断从属节点为脱开还是闭合,若从属节点为脱开,则释放节点约束;若从属节点为闭合,则施加节点约束;然后进行迭代计算,直至平衡收敛或达到最大迭代次数。
进一步的是,在检查平衡收敛之前,先检查从属节点上接触状态的变化,若从属节点在迭代后存在节点侵入情况,并且节点侵入距离>0,则它的状态由脱开变为闭合,然后施加节点约束,重新迭代计算;若从属节点在迭代后的节点压力<0,则它的状态则由闭合变为脱开,然后释放节点约束,重新迭代计算。
进一步的是,步骤4中,获得第二安全系数时,考虑了预应力支护和非预应力支护的相互作用机制,相互作用机制具体包括:先判断块体是否处于稳定状态,若块体处于稳定状态,则在块体未发生滑移时,将非预应力支护作为安全储备,不计入安全系数公式;当块体处于不稳定状态或发生滑移时,非预应力支护才发生作用,才将非预应力支护计入块体安全系数计算公式,作为块体极限安全分析,并评估加强支护措施。
进一步的是,所述块体包括双滑块体,所述双滑块体沿双滑面i和j滑动时,安全系数可按下列公式计算:
(1)双滑块体基本组合
安全系数K的计算公式为:
当合支护力分量Ph方向指向滑面i时:
当合支护力分量Ph方向指向滑面j时:
(2)双滑块体偶然组合
安全系数K的计算公式为:
当合支护力分量Ph方向指向滑面i时:
当合支护力分量Ph方向指向滑面j时:
式中,Ui为滑面i上的地下水扬压力,Uj为滑面j上的地下水扬压力;为滑面i的内摩擦角,为滑面j的内摩擦角;Pδ指向滑面i的合支护力P在通过滑面i和j交线的铅垂平面内的分力,且与水平面夹角为δ;γ为滑面i和j交线的倾角,δ合支护力在铅垂平面内的分力Pδ与水平面的夹角,俯角为正;ci为滑面i的黏聚力,cj滑面j的黏聚力;Ai为滑面i的面积,Aj为滑面j的面积;Uc为地下水产生的扬压力在沿滑面i和j交线方向的合力;Ph为指向滑面i的合支护力P在与滑面i和j交线垂直的平面内水平方向的分力;βi为滑面i与通过滑面i和j交线的铅垂辅助面的夹角,βj为滑面j与通过滑面i和j交线的铅垂辅助面的夹角。
进一步的是,步骤4中,计算和显示块体力系图具体方法包括:计算和显示块体的空间力系的分解图,三维空间向量方式表示块体的空间受力情况,并在块体三维有限元模型上显示空间力系的方向、大小以及滑动方向信息。
本发明对复杂地质条件下的块体,引入有限元接触力学分析方法,并考虑岩体和块体变形、采用刚体极限平衡法处理有限元计算结果,计算和评估块体稳定安全系数,并将次安全系数和在三维几何模型下采用刚体极限平衡法得出的安全系数进行比较,根据比价结果选择相应安全系数作为块体稳定的安全系数,极大地提高了对地下洞室块体稳定分析的准确性。
附图说明
图1是本发明的块体滑面接触分析流程图。
图2是本发明双滑块体沿双滑面i和j滑动时受力分析示意图。
附图中,1为滑面i,2为滑面j,3为水平参考面,4为通过滑面i和j交线的铅垂辅助面,5为与滑面i和j交线垂直的辅助面,6为其他张裂面k。
具体实施方式
本发明一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法,包括:
步骤1、建立块体三维几何模型,根据块体几何参数、材料参数、力学参数,采用刚体极限平衡法计算块体的安全系数,获得第一安全系数;
步骤2、建立块体三维有限元模型,根据块体的稳定状态和软弱结构面类型,获得块体滑面上节点和单元的应力位移结果;
步骤3、在块体三维有限元模型上,考虑块体滑面的黏聚力和摩擦系数,以及重力、预应力合支护力、非预应力合支护力、地下水以及地震力,对各滑面上的单元节点的应力,计算该节点处的法向力和切向力;
步骤4、根据块体滑面上节点和单元的应力位移结果以及节点处的法向力和切向力计算和显示块体的力系图,再根据刚体极限平衡法计算块体安全系数,获得第二安全系数;
步骤5、比较第一安全系数与第二安全系数的大小,取较小值作为块体稳定安全系数。
步骤2中,根据块体的稳定状态和软弱结构面类型,获得块体滑面上节点和单元的应力位移结果的具体方法包括:
对块体的稳定状态和软弱结构面类型进行判断,若块体处于稳定状态,则无需建立滑面单元,滑面两侧的实体单元边界共用滑面上的节点,进行计算分析,获得滑面上节点和单元的应力位移结果;
若块体处于不稳定状态且软弱夹层厚度小于等于阈值厚度,则不考虑结构面变形,在各滑面建立接触对单元,进行计算分析,获得滑面上节点和单元的应力位移结果;
若块体处于不稳定状态且软弱夹层厚度大于阈值厚度,则需考虑结构面变形,在各滑面建立实体单元,进行计算分析,获得滑面上节点和单元的应力位移结果。
块体滑面上的接触对单元在闭合状态下传递剪力时,考虑了块体软弱结构面类型,并根据软弱结构面的类型,采用不同单元模拟滑面。
块体滑面上的接触对单元在闭合状态下传递剪力时,考虑了块体滑面黏聚力。
对块体滑面上单元的接触状态判断,判断的具体流程如图1所示:在每一个增量步开始之前检查块体结构面上所有接触对相互作用状态,以判断从属节点为脱开还是闭合,若从属节点为脱开,则释放节点约束;若从属节点为闭合(即接触),则施加节点约束;然后进行迭代计算,直至平衡收敛或达到最大迭代次数;
在检查平衡收敛之前,先检查从属节点上接触状态的变化,若从属节点在迭代后存在节点侵入情况,并且节点侵入距离h>0,则它的状态由脱开变为闭合,然后施加节点约束,重新迭代计算;若从属节点在迭代后的节点压力p<0,则它的状态则由闭合变为脱开,然后释放节点约束,重新迭代计算。
步骤4中,获得第二安全系数时,考虑了预应力支护和非预应力支护的相互作用机制,相互作用机制具体包括:先判断块体是否处于稳定状态,若块体处于稳定状态,则在块体未发生滑移时,将非预应力支护作为安全储备,不计入安全系数公式;当块体处于不稳定状态或发生滑移时,非预应力支护才发生作用,才将非预应力支护计入块体安全系数计算公式,作为块体极限安全分析,并评估加强支护措施。
块体包括双滑块体,双滑块体沿双滑面i和j滑动时,其受力分析图如图2所示,安全系数可按下列公式计算:
(1)双滑块体基本组合
安全系数K的计算公式为:
当合支护力分量Ph方向指向滑面i1时:
当合支护力分量Ph方向指向滑面j2时:
(3)双滑块体偶然组合
安全系数K的计算公式为:
当合支护力分量Ph方向指向滑面i1时:
当合支护力分量Ph方向指向滑面j2时:
式中,Ui为滑面i1上的地下水扬压力,Uj为滑面j2上的地下水扬压力(单位为kN);为滑面i的内摩擦角,为滑面j的内摩擦角(单位°);Pδ指向滑面i的合支护力P在通过滑面i和j交线的铅垂平面内的分力(单位为kN),且与水平面夹角为δ;γ为滑面i和j交线的倾角(单位°),δ为合支护力在铅垂平面内的分力Pδ与水平面的夹角(单位°),俯角为正;δ'为非合支护力在铅垂平面内的分力P'δ与水平面的夹角(单位°),俯角为正;ci为滑面i的黏聚力(单位kPa),cj滑面j的黏聚力(单位kPa);Ai为滑面i的面积(m2),Aj为滑面j的面积(m2);Uc为地下水产生的扬压力在沿滑面i和j交线方向的合力(单位为kN);Ph为指向滑面i的合支护力P在与滑面i和j交线垂直的平面内水平方向的分力(单位为kN);βi为滑面i与通过滑面i和j交线的铅垂辅助面的夹角(单位°),βj为滑面j与通过滑面i和j交线的铅垂辅助面的夹角(单位°)。
其滑面i1以及滑面j2的位置参考水平面3,分析合支护力P在通过滑面i和j交线的铅垂平面内的分力时,可参考通过滑面i和j交线的铅垂辅助面4,在分析合支护力P在与滑面i和j交线垂直的平面内水平方向的分力时,可参考与滑面i和j交线垂直的辅助面5,以及在分析相关力时,可参考其他张裂面k6。
步骤4中,计算和显示块体力系图具体方法包括:计算和显示块体的空间力系的分解图,三维空间向量方式表示块体的空间受力情况,并在块体三维有限元模型上显示空间力系的方向、大小以及滑动方向信息。
本发明解决了刚体极限平衡法无法考虑不同滑面、不同区域的差异,尤其是双滑块体模式下的各结构面接触刚度差异、软弱结构面的黏聚力和变形影响等因素的适应性难题。
综上所述,本发明极大地提高了对地下洞室块体稳定分析的准确性。
Claims (7)
1.一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法,其特征在于,包括:
步骤1、建立块体三维几何模型,根据块体几何参数、材料参数、力学参数,采用刚体极限平衡法计算块体的安全系数,获得第一安全系数;
步骤2、建立块体三维有限元模型,根据块体的稳定状态和软弱结构面类型,获得块体滑面上节点和单元的应力位移结果;
步骤3、在块体三维有限元模型上,考虑块体滑面的黏聚力和摩擦系数,以及重力、预应力合支护力、非预应力合支护力、地下水以及地震力,对各滑面上的单元节点的应力,计算该节点处的法向力和切向力;
步骤4、根据块体滑面上节点和单元的应力位移结果以及节点处的法向力和切向力计算并显示块体的力系图,再根据刚体极限平衡法计算块体安全系数,获得第二安全系数;
步骤5、比较第一安全系数与第二安全系数的大小,取较小值作为块体稳定安全系数;
步骤4中,获得第二安全系数时,考虑了预应力支护和非预应力支护的相互作用机制,相互作用机制具体包括:先判断块体是否处于稳定状态,若块体处于稳定状态,则在块体未发生滑移时,将非预应力支护作为安全储备,不计入安全系数公式;当块体处于不稳定状态或发生滑移时,非预应力支护才发生作用,才将非预应力支护计入块体安全系数计算公式,作为块体极限安全分析,并评估加强支护措施;
所述块体包括双滑块体,所述双滑块体沿双滑面i和j滑动时,第二安全系数按下列公式计算:
(1)双滑块体基本组合
安全系数K的计算公式为:
当合支护力分量Ph方向指向滑面i时:
当合支护力分量Ph方向指向滑面j时:
(2)双滑块体偶然组合
安全系数K的计算公式为:
当合支护力分量Ph方向指向滑面i时:
当合支护力分量Ph方向指向滑面j时:
2.根据权利要求1所述的一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法,其特征在于,步骤2中,根据块体的稳定状态和软弱结构面类型,获得块体滑面上节点和单元的应力位移结果的具体方法包括:
对块体的稳定状态和软弱结构面类型进行判断,若块体处于稳定状态,则无需建立滑面单元,滑面两侧的实体单元边界共用滑面上的节点,进行计算分析,获得滑面上节点和单元的应力位移结果;
若块体处于不稳定状态且软弱夹层厚度小于等于阈值厚度,则不考虑结构面变形,在各滑面建立接触对单元,进行计算分析,获得滑面上节点和单元的应力位移结果;
若块体处于不稳定状态且软弱夹层厚度大于阈值厚度,则需考虑结构面变形,在各滑面建立实体单元,进行计算分析,获得滑面上节点和单元的应力位移结果。
3.根据权利要求2所述的一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法,其特征在于,块体滑面上的接触对单元在闭合状态下传递剪力时,考虑了块体软弱结构面类型,并根据软弱结构面的类型,采用不同单元模拟滑面。
4.根据权利要求2所述的一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法,其特征在于,块体滑面上的接触对单元在闭合状态下传递剪力时,考虑了块体滑面黏聚力。
5.根据权利要求2所述的一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法,其特征在于,对块体滑面上单元的接触状态判断,判断的具体方法包括:在每一个增量步开始之前检查块体结构面上所有接触对相互作用状态,以判断从属节点为脱开还是闭合,若从属节点为脱开,则释放节点约束;若从属节点为闭合,则施加节点约束;然后进行迭代计算,直至平衡收敛或达到最大迭代次数。
6.根据权利要求5所述的一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法,其特征在于,在检查平衡收敛之前,先检查从属节点上接触状态的变化,若从属节点在迭代后存在节点侵入情况,并且节点侵入距离>0,则它的状态由脱开变为闭合,然后施加节点约束,重新迭代计算;若从属节点在迭代后的节点压力<0,则它的状态则由闭合变为脱开,然后释放节点约束,重新迭代计算。
7.根据权利要求1所述的一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法,其特征在于,步骤4中,计算和显示块体力系图具体方法包括:计算和显示块体的空间力系的分解图,三维空间向量方式表示块体的空间受力情况,并在块体三维有限元模型上显示空间力系的方向、大小以及滑动方向信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111409221.0A CN114065590B (zh) | 2021-11-25 | 2021-11-25 | 一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111409221.0A CN114065590B (zh) | 2021-11-25 | 2021-11-25 | 一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114065590A CN114065590A (zh) | 2022-02-18 |
CN114065590B true CN114065590B (zh) | 2023-04-07 |
Family
ID=80276011
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111409221.0A Active CN114065590B (zh) | 2021-11-25 | 2021-11-25 | 一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114065590B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106436647A (zh) * | 2016-10-26 | 2017-02-22 | 南昌大学 | 一种装配式防波堤防护块体 |
CN107330145A (zh) * | 2017-05-26 | 2017-11-07 | 昆明理工大学 | 一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法 |
WO2019176835A1 (ja) * | 2018-03-13 | 2019-09-19 | 日本電気株式会社 | 斜面監視システム、斜面監視方法及び記録媒体 |
-
2021
- 2021-11-25 CN CN202111409221.0A patent/CN114065590B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106436647A (zh) * | 2016-10-26 | 2017-02-22 | 南昌大学 | 一种装配式防波堤防护块体 |
CN107330145A (zh) * | 2017-05-26 | 2017-11-07 | 昆明理工大学 | 一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法 |
WO2019176835A1 (ja) * | 2018-03-13 | 2019-09-19 | 日本電気株式会社 | 斜面監視システム、斜面監視方法及び記録媒体 |
Non-Patent Citations (5)
Title |
---|
刘耀儒 ; 杨强 ; 薛利军 ; 周维垣 ; .基于三维非线性有限元的边坡稳定分析方法.岩土力学.2007,(第09期),130-134. * |
张雨霆 ; 丁秀丽 ; 卢波 ; .基于数值分析的岩石块体稳定性评价一般性方法.岩石力学与工程学报.2017,(第01期),83-97. * |
曾海钊 ; 何江达 ; 谢红强 ; 左林勇 ; .基于块体理论的大型地下洞室围岩稳定性分析.贵州水力发电.2009,(第03期),25-30. * |
熊兆平 ; .某水电站地下厂房块体稳定性分析.水电与新能源.2016,(第02期),27-31. * |
连凌云 ; 侯哲生 ; .一种确定圆弧滑动面均质边坡安全系数的弹性极限平衡方法.水利科技与经济.2008,(第04期),19-22. * |
Also Published As
Publication number | Publication date |
---|---|
CN114065590A (zh) | 2022-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Explicit discontinuous deformation analysis method with lumped mass matrix for highly discrete block system | |
Lian et al. | A numerical study on toppling failure of a jointed rock slope by using the distinct lattice spring model | |
Jiang et al. | A rigorous solution for the stability of polyhedral rock blocks | |
Caballero et al. | Consistent tangent formulation for 3D interface modeling of cracking/fracture in quasi-brittle materials | |
Zhang et al. | New stability calculation method for rock slopes subject to flexural toppling failure | |
Li et al. | Soil-slope stability considering effect of soil-strength nonlinearity | |
Fu et al. | Computation of the safety factor for slope stability using discontinuous deformation analysis and the vector sum method | |
Zhang et al. | Extensions of edge-to-edge contact model in three-dimensional discontinuous deformation analysis for friction analysis | |
CN107341285A (zh) | 一种地震力作用下节理岩质边坡倾倒破坏极限承载力分析下限法 | |
Zhou et al. | Flow–stress coupled permeability tensor for fractured rock masses | |
Zhong et al. | Hydraulic fracture at the dam-foundation interface using the scaled boundary finite element method coupled with the cohesive crack model | |
Jiang et al. | A rigorous method for three-dimensional asymmetrical slope stability analysis | |
CN106126892B (zh) | 一种浆砌石挡土墙极限承载力的塑性极限分析上限法 | |
CN109800459B (zh) | 重力式挡土墙设计方法及装置 | |
CN107844650A (zh) | 基于全结构面屈服法的坝肩安全系数计算方法 | |
Zhao et al. | A novel contact algorithm based on a distance potential function for the 3D discrete-element method | |
Zhang et al. | Method for resolving contact indeterminacy in three-dimensional discontinuous deformation analysis | |
Fu et al. | Seismic stability analysis of a rock block using the block theory and Newmark method | |
Hashimoto et al. | Improvement of discontinuous deformation analysis incorporating implicit updating scheme of friction and joint strength degradation | |
Dang-Trung et al. | Numerical modeling of wing crack propagation accounting for fracture contact mechanics | |
CN114065590B (zh) | 一种融合刚体极限平衡法和有限元法的地下洞室块体稳定分析方法 | |
Fan et al. | An improved 3D DDA method considering the unloading effect of tunnel excavation and its application | |
He et al. | Investigation on strength and stability of jointed rock mass using three‐dimensional numerical manifold method | |
Liu et al. | Slope stability analysis based on a multigrid method using a nonlinear 3D finite element model | |
Asadollahi | Stability analysis of a single three dimensional rock block: effect of dilatancy and high-velocity water jet impact |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |