WO2023124664A1 - 岩体稳定性极限分析方法 - Google Patents
岩体稳定性极限分析方法 Download PDFInfo
- Publication number
- WO2023124664A1 WO2023124664A1 PCT/CN2022/133922 CN2022133922W WO2023124664A1 WO 2023124664 A1 WO2023124664 A1 WO 2023124664A1 CN 2022133922 W CN2022133922 W CN 2022133922W WO 2023124664 A1 WO2023124664 A1 WO 2023124664A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- stress
- unit
- rock mass
- node
- constraint equation
- Prior art date
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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
Definitions
- the limit equilibrium method is currently the most widely used method in the stability analysis of rock mass engineering. It is based on the classic Mohr-Coulomb yield criterion. For several rigid bodies, the safety factor or external load of the system can be obtained through the static force balance equation of the rigid body. Under normal circumstances, the calculation model is statically indeterminate, and the static force assumption is introduced to make the model statically indeterminate and solvable. The concept of the limit equilibrium method is clear, the calculation is simple, and it can give reasonable results in most cases, so it is widely used in engineering practice. However, the assumptions about the slip surface and the static assumptions introduced in the solution greatly damage the theoretical rigor of the method.
- the physical and mechanical parameters of the structural surface include internal friction angle and viscosity coefficient
- the present invention has the advantage of providing a rock mass stability analysis that does not need to introduce the assumption of the most dangerous sliding surface, does not need to introduce the stress-strain relationship, does not need to introduce the assumption of stress discontinuity or velocity discontinuity, and has a strict theoretical basis and strong applicability Numerical method, stability analysis conclusions are more in line with engineering practice, and have good engineering practical value.
- Figure 2 is a flow chart of the method of the present invention.
- Fig. 4 is a schematic diagram of unit node numbers, coordinates and stresses in the method of the present invention.
- Fig. 5 is a schematic diagram of known boundary conditions of nodal stresses in the method of the present invention.
- Example 1 As shown in Figure 1, the concrete gravity dam 4 is located on the rock foundation, and the concrete gravity dam 4 has a sliding force of 3.31 ⁇ 10 8 N acting on the upper surface of the rock mass 1 at the bottom section 6 of the dam, and the resistance of the rock foundation is analyzed.
- the problem of sliding stability illustrates the rock mass stability limit analysis method of the present invention, and its flow chart is as shown in Figure 2.
- the physical and mechanical parameters of rock mass include natural bulk density, dry bulk density, saturated bulk density, water content, tensile strength, compressive strength, internal friction angle, and viscosity coefficient;
- the physical and mechanical parameters of the structural surface include internal friction angle and viscosity coefficient
- the rock foundation stability analysis range takes a horizontal length of 660m, wherein the length of section 5 in front of the dam is 165m, the length of section 6 at the base of the dam is 165m, and the length of section 7 behind the dam is 330m; the vertical thickness of the rock foundation is taken as 198m, including the vertical thickness of rock mass 1 16.5m and the vertical thickness of rock mass 2 is 181.5m.
- the height of the concrete gravity dam 4 is 226m, the water depth in front of the dam is 220m, the water depth behind the dam is 11m, the concrete density is 24000N/m 3 , and the water body density is 9810N/m 3 .
- the quadrilateral unit is used for unit division.
- the model has a total of 3960 units and 4114 nodes. Nodes are numbered sequentially from top to bottom and from left to right.
- dx is the horizontal length of the unit
- dy is the vertical height of the unit
- the known boundary conditions of the nodal stress include the normal stress and shear stress on the upper surface of rock mass 1 in section 5 in front of the dam, the normal stress and sliding force on the upper surface of rock mass 1 in section 6 at the bottom of the dam, and the rock mass 1 in section 7 behind the dam. Normal stress and shear stress on the upper surface of body 1.
- the upper surface normal stress and shear stress of the section 5 rock mass 1 in front of the dam are produced by the water body in front of the dam, according to the known boundary condition constraint equation of the nodal stress described in the method of the present invention, the known boundary condition of the nodal stress can be calculated as
- Formula (16) is equivalent to the intersection of formula (17) and formula (18);
- k represents the kth node, and the value range is all nodes; Represents the normal stress in the X direction of the kth node; Represents the normal stress in the Y direction of the kth node; Represents the shear stress of the kth node; represents the internal friction angle; c represents the viscous force; ⁇ [0, ⁇ ) is the calculation parameter.
- rock mass 1 of present embodiment has
- the problem of foundation bearing capacity under strip uniformly distributed load is taken as an example to illustrate the effectiveness of the analysis method for rock mass stability limit in the present invention.
- the foundation 9 is a homogeneous rock mass, and there is no rock mass of different properties, no structural planes, and no weak layer;
- S3 carry out unit division; the total horizontal length of the foundation 9 is 60m, of which the strip uniformly distributed load 8 is 20m on the left and right, and the vertical thickness of the foundation 9 is 30m, and the unit division is carried out with quadrilateral units.
- the model is divided into 1152 units and 1225 nodes. Nodes are numbered sequentially from top to bottom and from left to right.
- Constraint equations include unit static equilibrium constraint equations, node stress known boundary condition constraint equations and node stress yield criterion constraint equations.
- Adopt aforementioned method of the present invention can obtain unit static force equilibrium constraint equation to be:
- the constraint equation of known node stress boundary conditions is constructed as follows:
- the known node stress boundary conditions include the normal stress and shear stress on the upper surface of the foundation 9 on the left side of the load 8, the normal stress and the shear stress on the upper surface of the foundation 9 at the bottom of the load 8 Shear stress, the normal stress and shear stress on the upper surface of the foundation 9 on the right side of the load 8.
- the node stress is known as the boundary condition:
- j represents the jth node
- j represents the jth node, and the value range is the node on the upper surface of the load 8 bottom foundation 9;
- the node stress is known as the boundary condition:
- the constraint equation of this embodiment is calculated and solved, and it can be obtained that under the condition that the strip uniformly distributed load 8 is not greater than 1023.4kPa, the constraint equation has a feasible solution, and it can be determined that the rock mass is stable; When it is greater than 1023.4kPa, there is no feasible solution to the constraint equation, and the rock mass can be judged to be unstable.
- the rock mass stability limit analysis method provided by this application provides a method that does not need to introduce the assumption of the most dangerous sliding surface, does not need to introduce the stress-strain relationship, does not need to introduce the assumption of stress discontinuity or velocity discontinuity, and the theory
- the numerical method of stability analysis of rock mass with strict foundation and strong applicability, the conclusion of stability analysis is more consistent with engineering practice, and has good engineering practical value.
Abstract
本发明公开了一种岩体稳定性极限分析方法，通过地质勘察，收集目标区域内不同性质岩体、构造面、软弱层的分布范围；并测定其物理力学参数，进行单元划分，以所有节点的应力分量为自变量，构建单元静力平衡约束方程、节点应力已知边界条件约束方程和节点应力屈服准则约束方程；之后求解同时满足约束方程的可行解；如果至少存在一个可行解，判定岩体为稳定；如果不存在可行解，判定岩体为不稳定。本发明优点在于提供了一种不需要引入最危险滑裂面假定、不需要引入应力-应变关系、不需要引入应力间断或速度间断假定、且理论基础严密、适用性强的岩体稳定性分析数值方法，稳定性分析结论与工程实践也更为符合，具有较好的工程实用价值。
Description
本发明涉及岩土工程稳定性分析领域，尤其是涉及岩体稳定性极限分析方法。
岩体工程稳定性分析在水利水电、公路与轨道交通、工业与民用建筑、地质矿产开发及生态环境治理等工程中广泛存在。岩体稳定性分析主要包括三类问题：岩质边坡稳定性分析、地下洞室围岩稳定性分析和岩石基础承载力分析；所采用的稳定性分析方法主要有极限平衡法、滑移线法、有限元法和极限分析法等。
极限平衡法是目前岩体工程稳定性分析中应用最为广泛的一种方法，它以经典的莫尔-库仑屈服准则为基础，首先假定一个发生剪切破坏的滑裂面，将破坏区域划分成若干刚性体，通过刚性体的静力平衡方程，求得系统的安全系数或外荷载。在通常情况下，计算模型是静不定的，求解时通过引入静力假定，使模型变得静定可解。极限平衡法概念清晰、计算简单，在多数情况下能给出合理的结果，在工程实践中应用广泛。但有关滑裂面的假定和求解时引入的静力假定，使该方法理论上的严密性受到较大损害。
滑移线法基于静力平衡方程、屈服条件和应力边界条件，通过求解塑性区的应力与位移速度分布得到相应的极限荷载。滑移线法理论基础严密，在不考虑变形、应变软化或应变硬化的情况下，对某些具有特定边界条件并服从特定屈服准则的稳定性问题，滑移线法可以获得理论上的闭合解，这些闭合解为开发新的理论体系严密、适用性强的稳定性分析方法提供了极为重要的验证手段，具有十分重要的理论意义。但由于实际问题中边界条件的复杂性等其他原因，滑移线法在工程实践中难以直接应用。
有限元法把连续体划分成有限个单元，把单元的节点作为离散点，基于平衡方程和本构关系进行求解和分析。有限元法物理概念清晰，理论体系严密，能够 成功地适应材料的非均匀、非线性、复杂的几何形状以及复杂的边界条件等，应用十分广泛。进行岩体稳定性分析时，有限元法不需要事先假定滑裂面的形状或位置，且能够提供应力与变形的全部信息，有其自身的理论优势；但实际岩体由于受到构造因素的影响，本构关系最难以弄清，有限元法需要在最不易弄清的本构关系的基础上进行求解和稳定性分析，使得有限元法在岩体稳定性分析中的应用受到了很大限制。
极限分析法基于塑性力学中的上、下限定理，将极限分析与有限单元法相结合，在单元的边界面上引入应力间断或速度间断假定，将稳定性分析问题转化为约束条件下的最优化问题，得到了较为满意的结果。但极限分析法在求解时引入的应力间断或速度间断的假定，使该方法在理论上的严密性受到一定程度的损害。
因此，设计一种不需要引入最危险滑裂面假定、不需要引入应力-应变关系、不需要引入应力间断或速度间断假定、且理论基础严密、适用性强的新的岩体稳定性分析方法，对岩质边坡稳定性分析、地下洞室围岩稳定性分析和岩石基础稳定性分析等有着重要意义。
发明内容
本发明目的在于提供一种岩体稳定性极限分析方法。
为实现上述目的，本发明采取下述技术方案：
本发明所述一种岩体稳定性极限分析方法，包括以下步骤：
S1，通过地质勘察，收集目标区域内不同性质岩体、构造面、软弱层的分布范围；
S2，测定不同性质所述岩体、所述构造面、所述软弱层的物理力学参数；
S3，根据不同性质岩体、构造面、软弱层的分布范围，进行单元划分；
S4，以所有节点的应力分量为自变量，构建单元静力平衡约束方程、节点应力已知边界条件约束方程和节点应力屈服准则约束方程；
S5，求解同时满足所述单元静力平衡约束方程、所述节点应力已知边界条件 约束方程和所述节点应力屈服准则约束方程的可行解；
S6，分析岩体的稳定性；如果至少存在一个可行解，判定岩体为稳定；如果不存在可行解，判定岩体为不稳定。
进一步地，S2步中，所述岩体物理力学参数包括天然容重、干容重、饱和容重、含水率、抗拉强度、抗压强度、内摩擦角、粘滞系数；
所述构造面的物理力学参数包括内摩擦角、粘滞系数；
所述软弱层的物理力学参数包括天然容重、干容重、饱和容重、含水率、内摩擦角、粘滞系数。
进一步地，S3步中，所述单元划分中的单元采用四边形单元。
进一步地，所述单元静力平衡约束方程为：
其中，i代表第i个单元，取值范围为所有单元；
代表第i个单元周边节点的正应力分量在X方向对该单元的合力；
代表第i个单元周边节点的切应力分量在X方向对该单元的合力；
代表第i个单元在X方向的体力；
代表第i个单元周边节点的正应力分量在Y方向对该单元的合力；
代表第i个单元周边节点的切应力分量在Y方向对该单元的合力；
代表第i个单元在Y方向的体力。
进一步地，所述节点应力已知边界条件约束方程为:
其中，j代表第j个节点，取值范围为应力已知边界范围内的节点；
代表第j个节点X方向的正应力；
代表第j个节点Y方向的正应力；
代表第j个节点的切应力；
代表第j个节点在边界方向的正应力；
代表第j个节点在边界方向的切应力；l＝cosα、m＝sinα，α∈[0,π)，为计算参量。
进一步地，所述节点应力屈服准则约束方程为：
其中，k代表第k个节点，取值范围为所有节点；
代表第k个节点X方向的正应力；
代表第k个节点Y方向的正应力；
代表第k个节点的切应力；
代表内摩擦角；c代表粘滞力；α∈[0,π)，为计算参量。
进一步地，在三维岩体稳定性分析中，所述单元采用六面体单元进行划分，并依据所述单元静力平衡约束方程、所述节点应力已知边界条件约束方程和所述节点应力屈服准则约束方程的建立思路，建立三维的单元静力平衡约束方程、节点应力已知边界条件约束方程和节点应力屈服准则约束方程，求解所述约束方程的可行解，分析岩体的稳定性。
本发明优点在于提供了一种不需要引入最危险滑裂面假定、不需要引入应力-应变关系、不需要引入应力间断或速度间断假定、且理论基础严密、适用性强的岩体稳定性分析数值方法，稳定性分析结论与工程实践也更为符合，具有较好的工程实用价值。
图1是本发明所述实施例1中混凝土重力坝与岩石基础示意图。
图2是本发明所述方法的流程图。
图3是本发明所述实施例1中单元划分示意图。
图4是本发明所述方法中单元节点编号、坐标及应力示意图。
图5是本发明所述方法的节点应力已知边界条件示意图。
图6是本发明所述实施例2中条形均布荷载作用下的地基承载力示意图。
下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1：如图1所示，混凝土重力坝4坐落在岩石基础上，混凝土重力坝4在坝底段6岩体1上表面作用有3.31×10
8N的滑动力，分析岩石基础的抗滑稳定性问题，说明本发明所述岩体稳定性极限分析方法，其流程如图2所示。
S1，通过地质勘察，收集目标区域内不同性质岩体、构造面、软弱层的分布范围；
S2，测定不同性质岩体、构造面、软弱层的物理力学参数；
其中，岩体的物理力学参数包括天然容重、干容重、饱和容重、含水率、抗拉强度、抗压强度、内摩擦角、粘滞系数；
构造面的物理力学参数包括内摩擦角、粘滞系数；
软弱层的物理力学参数包括天然容重、干容重、饱和容重、含水率、内摩擦角、粘滞系数。
在本实施例中，通过地质勘察查明混凝土重力坝4岩石基础有两层不同性质的岩体，岩体1和岩体2，岩体1与岩体2交界面处有软弱层3，岩体1、岩体2与软弱层3均呈水平分布；岩体1竖向厚度16.5m，饱和容重25300N/m
3，粘滞力1.20MPa，内摩擦角47.73°；岩体2饱和容重25600N/m
3，粘滞力1.50MPa，内摩擦角50.19°；软弱层3的粘滞力0.21MPa，内摩擦角34.99度。
S3，根据不同性质岩体、构造面、软弱层的分布范围，进行单元划分；
如图3所示，本实施例仅进行二维分析，采用四边形单元进行单元划分，在三维分析中可以采用六面体单元进行单元划分。本实施例岩石基础稳定性分析范围取水平长度660m，其中坝前段5长度165m，坝基段6长度165m，坝后段7长度330m；取岩石基础竖向厚度198m，其中包括岩体1竖向厚度16.5m和岩体2竖向厚度为181.5m。
混凝土重力坝4高度为226m，取坝前水深220m，坝后水深11m，混凝土容重24000N/m
3，水体容重9810N/m
3。
采用四边形单元进行单元划分。水平长度范围内共划分120个单元，121个 节点，单元水平长度dx＝5.5m；竖向厚度范围内共划分33个单元，34个节点，单元竖向高度dy＝5.5m。模型共有3960个单元，4114个节点。节点按照从上到下，从左到右依次编号。
S4，以所有单元节点的应力分量为自变量，构建约束方程；
对本实施例，自变量为：
其中，i代表第i个节点，取值范围为所有节点，即i＝1，2,……4114；
分别代表第i个节点X方向的正应力、Y方向的正应力和第i个节点的切应力；共有4114*3＝12342个自变量。
约束方程包括单元静力平衡约束方程、节点应力已知边界条件约束方程和节点应力屈服准则约束方程。
其中单元静力平衡约束方程为：
其中，i代表第i个单元，取值范围为所有单元；
代表第i个单元周边节点的正应力分量在X方向对该单元的合力；
代表第i个单元周边节点的切应力分量在X方向对该单元的合力；
代表第i个单元在X方向的体力；
代表第i个单元周边节点的正应力分量在Y方向对该单元的合力；
代表第i个单元周边节点的切应力分量在Y方向对该单元的合力；
代表第i个单元在Y方向的体力。
根据本发明所述的单元静力平衡约束方程，设本实施例第i个单元的节点按逆时针方向编号依次为1，2，3，4；如图4所示，相应节点坐标为
和
相应节点在X方向、Y方向的正应力，以及相应节点的切应力记为
和
其中i的取值范围为所有单元。
设相邻节点之间的正应力和切应力在单元的边界上均按直线变化，则可得到，第i个单元节点在X方向的正应力
和
第i个单元周边节点的 正应力分量在X方向对该单元的合力
为：
第i个单元节点的切应力
和
第i个单元周边节点的切应力分量在X方向对该单元的合力
为：
取单元的厚度为单位厚度，设单元在X方向上的单位体力为γ
X，在Y方向上的单位体力为γ
Y；第i个单元的面积为S
i＝dx*dy；第i个单元在X方向的体力
为：
在本实施例中，第i个单元在X方向上的体力
为0；
第i个单元节点在Y方向的正应力
和
第i个单元周边节点的正应力分量在Y方向对该单元的合力
为：
第i个单元节点在Y方向的切应力
和
第i个单元周边节点的切应力分量在Y方向对该单元的合力
为：
第i个单元在Y方向的体力
为：
其中，dx为单元水平长度，dy为单元竖向高度；
将公式(3)、公式(4)、公式5)、公式6)、公式(7)和公式(8)，以及单元水平长度dx＝5.5m和单元竖向高度为dy＝5.5m代入单元静力平衡约束方程，即公式(2)，可得到本实施例的第i个单元在X方向和Y方向力的单元静力平衡约束方程为：
其中，i的取值范围为所有单元，即i＝1，2,……3960。
如图5所示，所述节点应力已知边界条件约束方程为
其中，j代表第j个节点，取值范围为应力已知边界范围内的节点；
代表第j个节点X方向的正应力；
代表第j个节点Y方向的正应力；
代表第j个节点的切应力；
代表第j个节点在边界方向的正应力；
代表第j个节点在边界方向的切应力；l＝cosα、m＝sinα，α∈[0,π)，为计算参量。
在本实施例中，节点应力已知边界条件包括坝前段5岩体1上表面的正应力和切应力，坝底段6岩体1上表面的正应力与滑动力，及坝后段7岩体1上表面的正应力和切应力。
坝前段5岩体1上表面正应力和切应力由坝前水体产生，根据本发明方法所述的节点应力已知边界条件约束方程，可以计算其节点应力已知边界条件为
其中，j表第j个节点，取值范围为坝前段5岩体1上表面的节点；根据本实施例对节点的编号，j的具体取值为：j＝1，2,……30。
坝底段6岩体1上表面正应力与滑动力，其节点应力已知边界条件为
其中，j表第j个节点，取值范围为坝底段6岩体1上表面的节点；根据本实施例对节点的编号，j的具体取值为：j＝31，32,……61。
坝后段7岩体1上表面正应力和切应力由坝后水体产生，其节点应力已知边界条件为
其中，j表第j个节点，取值范围为坝后段7岩体1上表面的节点；根据本实施例对节点的编号，j的具体取值为：j＝62，63,……121。
上述的公式(10)、公式(11)和公式(12)即为本实施例的节点应力已知边界条件约束方程。
所述节点应力屈服准则约束方程，是根据勘察和试验得到的不同性质岩体、构造面及软弱层的物理力学参数，依据相应的岩体强度理论，建立以应力形式表达的节点应力屈服准则约束方程。
以二维问题、莫尔-库仑屈服准则为例，节点k所对应的节点应力分量为
和
该点在方向l＝cosα、m＝sinα上的正应力和切应力分别为
和
则有：
由勘察及试验测定出的岩体内摩擦角为
并记
其节点应力在任一方向上的正应力
和切应力
均应满足莫尔-库仑准则，则可得到：
其中，c代表岩体的粘滞力；同时，公式(15)等价于
若记：
公式(16)等价于公式(17)和公式(18)的交集；
将公式(13)、公式(14)代入公式(17)，则可得到：
将l＝cosα、m＝sinα代入公式(19)中，进一步可得到：
将公式(13)、公式(14)代入公式(18)，则可得到：
将l＝cosα、m＝sinα代入公式(18)中，进一步可得到：
综合公式(20)与公式(22)，可得到
即本发明方法所述的节点应力屈服准则约束方程。其中，k代表第k个节点，取值范围为所有节点；
代表第k个节点X方向的正应力；
代表第k个节点Y方向的正应力；
代表第k个节点的切应力；
代表内摩擦角；c代表粘滞力；α∈[0,π)，为计算参量。
对本实施例，应用本发明方法所述的节点应力屈服准则约束方程，即有
根据勘察和试验，对本实施例岩体1有
c＝1.2×10
6；
f＝1.1001；
k的取值范围为：k＝1，2,……363；
α的取值范围为：
对软弱层3有：
c＝0.21×10
6；
f＝0.6999；
k的取值范围为：k＝364，365,……484；
α的取值范围为：
对岩体2有：
c＝1.5×10
6；
f＝1.1998；
k的取值范围为：k＝485，486,……4114。
上述公式(9)、公式(10)、公式(11)、公式(12)和公式(24)即为本实施例的约束方程，至此完成了本实施例约束方程的构建。
S5，求解同时满足所述单元静力平衡约束方程、所述节点应力已知边界条件约束方程和所述节点应力屈服准则约束方程的可行解；
S6，分析岩体的稳定性；如果至少存在一个可行解，判定岩体为稳定；如果不存在可行解，判定岩体为不稳定。
在本实施例中，通过计算机对约束方程求解可得到，在混凝土重力坝4坝底段6岩体1上表面作用有3.31×10
8N滑动力的情况下，不存在可行解，可判定该岩体为不稳定。
进一步计算可发现，如果要使得约束方程至少存在一个可行解，则岩体1上 表面作用的滑动力应不大于2.77×10
8N。
进一步地，针对本实施例，假定岩体1与岩体2间无软弱层3，在岩体1上表面作用有3.31×10
8N滑动力的情况下，通过计算机对约束方程求解发现，约束方程存在可行解，可判定该岩体为稳定。
进一步计算可发现，在假定岩体1与岩体2间无软弱层3的情况下，在岩体1上表面作用的滑动力不大于3.58×10
8N情况下，约束方程存在可行解；在岩体1上表面作用的滑动力大于3.58×10
8N情况下，约束方程不存在可行解。
对比上述四种求解结果可以看出，由于受到岩体1与岩体2间软弱层3的影响，岩石基础可提供的最大抗滑力由3.58×10
8N下降至2.77×10
8N，下降了22.6％，软弱层3的存在使得岩石基础的抗滑稳定性不能满足工程要求。
实施例2
本实施例以条形均布荷载作用下的地基承载力问题为例，说明本发明所述岩体稳定性极限分析方法的有效性。
如图6所示，在基础9的顶部20m范围内作用有条形均布荷载8，基础9的粘滞力c为200kPa，内摩擦角
为零，体力为零。下面采用本发明所述方法来求解基础9的地基承载力，即条形均布荷载8的最大值。
S1，收集目标区域内不同性质岩体、构造面、软弱层的分布范围；对本实施例，基础9为均质岩体，在基础范围内不存在不同性质岩体、不存在构造面、不存在软弱层；
S2，获得不同性质岩体、构造面、软弱层的物理力学参数；基础9岩体的粘滞力c为200kPa，内摩擦角φ为零，岩体容重为零；
S3，进行单元划分；基础9的水平总长度为60m，其中条形均布荷载8左、右各20m，基础9竖向厚度为30m，采用四边形单元进行单元划分。水平长度范围内共划分48个单元，49个节点，单元水平长度dx＝1.25m；竖向厚度范围内共划分24个单元，25个节点，单元竖向高度为dy＝1.25m。模型共划分为1152个 单元，1225个节点。节点按照从上到下，从左到右依次编号。
S4，以所有节点的应力分量为自变量，构建约束方程；
对本实施例，自变量为：
其中，i代表第i个节点，取值范围为所有节点，即i＝1，2,……1225；
分别代表第i个节点X方向的正应力、Y方向的正应力和第i个节点的切应力；共有1225*3＝3675个自变量。
约束方程包括单元静力平衡约束方程、节点应力已知边界条件约束方程和节点应力屈服准则约束方程。
单元静力平衡约束方程构建如下：对本实施例，单元水平长度dx＝1.25m，单元竖向高度dy＝1.25m，单元在X方向和Y方向力的体力均为零。采用本发明前述方法，可得到单元静力平衡约束方程为：
其中，i的取值范围为所有单元，即i＝1，2,……1152。
节点应力已知边界条件约束方程构建如下：在本实施例中，节点应力已知边界条件包括荷载8左侧基础9上表面的正应力和切应力，荷载8底部基础9上表面的正应力和切应力，荷载8右侧基础9上表面的正应力和切应力。
荷载8左侧基础9上表面的节点，其节点应力已知边界条件为
其中，j表第j个节点，取值范围为荷载8左侧基础9上表面的的节点；根据本实施例对节点的编号，j的具体取值为：j＝1，2,……16。
荷载8底部基础9上表面的节点，其节点应力已知边界条件为
其中，j表第j个节点，取值范围为荷载8底部基础9上表面的的节点；p
k为设定的地基承载力；根据本实施例对节点的编号，j的具体取值为：j＝17,18,……33。
荷载8右侧基础9上表面的节点，其节点应力已知边界条件为
其中，j表第j个节点，取值范围为荷载8右侧基础9上表面的节点；根据本实施例对节点的编号，j的具体取值为：j＝34,35,……49。
上述的公式(27)、公式(28)和公式(29)即为本实施例的基点应力已知边界条件约束方程。
节点应力屈服准则约束方程构建如下：对本实施例，基础9岩体的粘滞力c为200kPa，内摩擦角φ为零，f＝0，体力为零，其节点应力屈服准则约束方程为
k的取值范围为：k＝1,2,……1225。
上述公式(26)、公式(27)、公式(28)、公式(29)和公式(30)即为本实施例的约束方程，至此完成了本实施例约束方程的构建。
S5，求解同时满足所述单元静力平衡约束方程、所述节点应力已知边界条件约束方程和所述节点应力屈服准则约束方程的可行解；
S6，分析岩体的稳定性；
通过计算机，对本实施例约束方程进行计算求解可得到，在条形均布荷载8 不大于1023.4kPa的情况下，约束方程存在可行解，可判定该岩体为稳定；在条形均布荷载8大于1023.4kPa的情况下，约束方程不存在可行解，可判定该岩体为不稳定。
根据上述计算结果，可得到基础9的地基承载力为1023.4kPa。同时，根据塑性力学理论，基础9的地基承载力理论解为q＝c(π+2)，即为1028.3kPa。用本发明所述的方法所得到的解与理论解的相对误差为0.48％，精度较高，能够很好地满足工程实践需要。
实施例3
在三维岩体稳定性分析中，所述单元采用六面体单元进行划分，并依据所述单元静力平衡约束方程、所述节点应力已知边界条件约束方程和所述节点应力屈服准则约束方程的建立思路，建立三维的单元静力平衡约束方程、节点应力已知边界条件约束方程和节点应力屈服准则约束方程，求解所述约束方程的可行解，分析岩体的稳定性；如果至少存在一个可行解，判定岩体为稳定；如果可行解不存在，判定岩体为不稳定。
综上，本申请所提供的岩体稳定性极限分析方法，提供了一种不需要引入最危险滑裂面假定、不需要引入应力-应变关系、不需要引入应力间断或速度间断假定、且理论基础严密、适用性强的岩体稳定性分析数值方法，稳定性分析结论与工程实践也更为符合，具有较好的工程实用价值。
Claims (7)
- 一种岩体稳定性极限分析方法，其特征在于：包括以下步骤：S1，通过地质勘察，收集目标区域内不同性质岩体、构造面、软弱层的分布范围；S2，测定不同性质所述岩体、所述构造面、所述软弱层的物理力学参数；S3，根据不同性质岩体、构造面、软弱层的分布范围，进行单元划分；S4，以所有节点的应力分量为自变量，构建单元静力平衡约束方程、节点应力已知边界条件约束方程和节点应力屈服准则约束方程；S5，求解同时满足所述单元静力平衡约束方程、所述节点应力已知边界条件约束方程和所述节点应力屈服准则约束方程的可行解；S6，分析岩体的稳定性；如果至少存在一个所述可行解，判定岩体为稳定；如果不存在可行解，判定岩体为不稳定。
- 根据权利要求1所述的岩体稳定性极限分析方法，其特征在于：S2步中，所述岩体物理力学参数包括天然容重、干容重、饱和容重、含水率、抗拉强度、抗压强度、内摩擦角、粘滞系数；所述构造面的物理力学参数包括内摩擦角、粘滞系数；所述软弱层的物理力学参数包括天然容重、干容重、饱和容重、含水率、内摩擦角、粘滞系数。
- 根据权利要求1所述的岩体稳定性极限分析方法，其特征在于：S3步中，所述单元划分中的单元采用四边形单元。
- 根据权利要求1所述的岩体稳定性极限分析方法，其特征在于：所述单元静力平衡约束方程为：
其中，i代表第i个单元，取值范围为所有单元；代表第i个单元周边节点的正应力分量在X方向对该单元的合力；
代表第i个单元周边节点的切应力分 量在X方向对该单元的合力；
代表第i个单元在X方向的体力；
代表第i个单元周边节点的正应力分量在Y方向对该单元的合力；
代表第i个单元周边节点的切应力分量在Y方向对该单元的合力；
代表第i个单元在Y方向的体力。
- 根据权利要求1所述的岩体稳定性极限分析方法，其特征在于：所述节点应力已知边界条件约束方程为:
其中，j代表第j个节点，取值范围为应力已知边界范围内的节点；代表第j个节点X方向的正应力；
代表第j个节点Y方向的正应力；
代表第j个节点的切应力；
代表第j个节点在边界方向的正应力；
代表第j个节点在边界方向的切应力；l＝cosα、m＝sinα，α∈[0,π)，为计算参量。
- 根据权利要求1所述的岩体稳定性极限分析方法，其特征在于：所述节点应力屈服准则约束方程为：
其中，k代表第k个节点，取值范围为所有节点；代表第k个节点X方向的正应力；
代表第k个节点Y方向的正应力；
代表第k个节点的切应力；
代表内摩擦角；c代表粘滞力；α∈[0,π)，为计算参量。
- 根据权利要求1-6所述的岩体稳定性极限分析方法，其特征在于：在三维岩体稳定性分析中，所述单元采用六面体单元进行划分，并依据所述单元静力平衡约束方程、所述节点应力已知边界条件约束方程和所述节点应力屈服准则约束方程的建立思路，建立三维的单元静力平衡约束方程、节点应力已知边界条件约束方程和节点应力屈服准则约束方程，求解所述约束方程的所述可行解，分析岩体的稳定性。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111610587.4A CN114254519A (zh) | 2021-12-27 | 2021-12-27 | 岩体稳定性极限分析方法 |
| CN202111610587.4 | 2021-12-27 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2023124664A1 true WO2023124664A1 (zh) | 2023-07-06 |
Family
ID=80798089
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/CN2022/133922 WO2023124664A1 (zh) | 2021-12-27 | 2022-11-24 | 岩体稳定性极限分析方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN114254519A (zh) |
| WO (1) | WO2023124664A1 (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114254519A (zh) * | 2021-12-27 | 2022-03-29 | 黄河勘测规划设计研究院有限公司 | 岩体稳定性极限分析方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109918841A (zh) * | 2019-03-28 | 2019-06-21 | 黄河勘测规划设计研究院有限公司 | 基于有限单元与极限平衡的土体稳定性评价方法 |
| CN112069574A (zh) * | 2020-08-26 | 2020-12-11 | 盐城工学院 | 基于极限分析上限定理的圆形基坑稳定性分析方法 |
| CN114254519A (zh) * | 2021-12-27 | 2022-03-29 | 黄河勘测规划设计研究院有限公司 | 岩体稳定性极限分析方法 |
| CN114662341A (zh) * | 2022-05-16 | 2022-06-24 | 黄河勘测规划设计研究院有限公司 | 岩体临界滑动面极限分析方法 |
-
2021
- 2021-12-27 CN CN202111610587.4A patent/CN114254519A/zh active Pending
-
2022
- 2022-11-24 WO PCT/CN2022/133922 patent/WO2023124664A1/zh unknown
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109918841A (zh) * | 2019-03-28 | 2019-06-21 | 黄河勘测规划设计研究院有限公司 | 基于有限单元与极限平衡的土体稳定性评价方法 |
| CN112069574A (zh) * | 2020-08-26 | 2020-12-11 | 盐城工学院 | 基于极限分析上限定理的圆形基坑稳定性分析方法 |
| CN114254519A (zh) * | 2021-12-27 | 2022-03-29 | 黄河勘测规划设计研究院有限公司 | 岩体稳定性极限分析方法 |
| CN114662341A (zh) * | 2022-05-16 | 2022-06-24 | 黄河勘测规划设计研究院有限公司 | 岩体临界滑动面极限分析方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114254519A (zh) | 2022-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yu et al. | An improved form of smoothed particle hydrodynamics method for crack propagation simulation applied in rock mechanics | |
| Khoei et al. | Modeling of dynamic cohesive fracture propagation in porous saturated media | |
| CN109505576A (zh) | 页岩水力压裂三维全耦合离散裂缝网络模拟方法及系统 | |
| WO2023124664A1 (zh) | 岩体稳定性极限分析方法 | |
| Caselunghe et al. | Structural element approaches for soil-structure interaction | |
| Chen et al. | Seismic response and damage of underground subway station in a slightly sloping liquefiable site | |
| Aron et al. | Structural element approaches for soil-structure interaction | |
| Wang et al. | Unsaturated seepage–stress–damage coupling and dynamic analysis of stability on discrete fractured rock slope | |
| Rokonuzzaman et al. | Model tests and 3D finite element simulations of uplift resistance of shallow rectangular anchor foundations | |
| Rawat et al. | Investigation on failure mechanism of nailed soil slopes under surcharge loading: testing and analysis | |
| CN106321082A (zh) | 用于求取盐膏岩地层最大套管外载荷的方法及系统 | |
| Tao et al. | Failure and stability analysis of Jinping-I arch dam based on geomechanical model test and nonlinear numerical analysis | |
| Audenaert et al. | An analytical model to determine the ultimate load on masonry arch bridges. | |
| Zhang et al. | Creep characteristics and model of key unit rock in slope potential slip surface | |
| CN114662341A (zh) | 岩体临界滑动面极限分析方法 | |
| Lu et al. | Numerical and experimental analyses for bearing capacity of rigid strip footing subjected to eccentric load | |
| Pan et al. | On the characteristics of ground motion and the improvement of the input mode of complex layered sites | |
| CN107066744A (zh) | 富水地层中地铁车站的抗震分析方法 | |
| Cheng et al. | Shaking table test and numerical verification for free ground seismic response of saturated soft soil | |
| Sheng et al. | Modified Bishop method for stability analysis of weakly sloped subgrade under centrifuge model test | |
| Cheng et al. | Seismic stability of a tunnel considering the dynamic geologic parameters of loess | |
| Zhang et al. | Exploring the Deformation and Failure Characteristics of Interbedded Anti-inclined Rock Slopes: Insights from Physical Modelling Tests | |
| Pan et al. | The forward and inversion analysis of high rock-fill dam during construction period using the node-based smoothed point interpolation method | |
| Cheng et al. | Study on the Optimization of Pile Length of Micropiles in Soil Landslides | |
| Muthukkumaran et al. | Effect of seabed slope on pile behaviour of fixed offshore platform under lateral forces |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 22913922 Country of ref document: EP Kind code of ref document: A1 |