Rock mass stability analysis method, device, terminal and storage medium

Abstract

The application provides a rock mass stability analysis method, device, terminal and storage medium. The method comprises the following steps: searching all movable blocks in a target rock mass; calculating sub-contact forces of each structural surface of a target movable block; calculating a resultant force of active forces suffered by the target movable block according to the gravity of the target movable block and the sub-contact forces of each structural surface of the target movable block; calculating a safety factor of the target movable block according to the resultant force of the active forces and a preset safety factor calculation formula corresponding to a plurality of motion modes; and generating a stability analysis result of the target rock mass according to the safety factor of the target movable block. The application considers the influence of the interaction force between blocks on the translational and rotational stability of key blocks, and has a good effect on the stability analysis of the rock mass.

Description

Technical Field

[0001] This invention relates to the field of rock mass engineering technology, and in particular to a method, apparatus, terminal and storage medium for rock mass stability analysis. Background Technology

[0002] In recent years, numerous large-scale water conservancy, hydropower, and transportation projects have been constructed both domestically and internationally. Among these, the stability of rock masses in engineering structures such as slopes, tunnels, and underground chambers is one of the key issues that needs to be addressed for the safe construction and normal operation of these projects. Rock masses are natural geological formations composed of structural planes and structural elements. Due to the weakening effect of structural planes on rock mass strength and the impact of excavation disturbance on rock mass stability, rock blocks at free faces may detach, slide, or rotate and fail. Therefore, to ensure the safety of engineering structures, the primary issue is accurately analyzing rock mass stability.

[0003] Currently, there are two main categories of theoretical analysis methods for rock mass stability: one is the continuous medium analysis method, which treats the rock mass as a continuous medium and uses elastoplastic mechanics as its main basis, employing techniques such as the finite element method, finite difference method, and boundary element method to analyze the stress and displacement fields of the rock mass; the other is the discontinuous medium analysis method, which treats the rock mass as a discontinuous medium and uses techniques such as the discrete element method, block theory, discontinuous deformation analysis methods, and numerical manifold methods to analyze rock mass stability. Among the discontinuous medium analysis methods, the critical block theory does not require complex mesh generation and stress-displacement analysis; it only searches for unstable blocks on the free surface. It has the advantages of complete theory and simple and fast calculation, and therefore has been widely used.

[0004] However, current research and application of critical block theory usually do not consider the influence of inter-block forces on the translational and rotational stability of critical blocks. This results in the calculated sliding force and rotational torque being too small, and some unstable blocks being ignored, leading to poor stability analysis of rock masses and support designs that do not meet actual engineering requirements. Summary of the Invention

[0005] This invention provides a method, apparatus, equipment, and computer storage medium for rock mass stability analysis, which can improve the effectiveness of rock mass stability analysis.

[0006] In a first aspect, embodiments of the present invention provide a method for rock mass stability analysis, comprising:

[0007] Locate all movable blocks within the target rock mass;

[0008] Calculate the contact forces of each structural surface of the target movable block; where the target movable block is any one of the movable blocks, and the contact forces include normal contact forces and tangential contact forces.

[0009] Calculate the resultant force of the principal driving force on the target movable block based on the gravity of the target movable block and the contact force of each structural surface of the target movable block;

[0010] The safety factor of the target movable block is calculated based on the resultant force of the active force and the preset safety factor calculation formula corresponding to multiple motion modes; among them, the multiple motion modes include at least the double-sided sliding mode and the edge rotation mode.

[0011] Based on the safety factor of the target movable block, the stability analysis results of the target rock mass are generated.

[0012] In one possible implementation, the contact forces of each structural surface of the target movable block are calculated, including:

[0013] Displacement constraints are applied to all movable blocks on free faces within the target rock mass.

[0014] Calculate the contact forces of each structural surface of the target movable block based on all free-face movable blocks under displacement constraints.

[0015] In one possible implementation, based on all free-faced movable blocks under displacement constraints, the contact forces of each structural surface of the target movable block are calculated, including:

[0016] Calculate the contact velocity of each structural surface of the target movable block according to the preset motion equation;

[0017] Calculate the relative displacement increment of each structural surface of the target movable block based on the contact velocity of each structural surface;

[0018] The contact force of each structural surface of the target movable block is calculated based on the relative displacement increment of each structural surface of the target movable block, the area of ​​each structural surface contact point of the target movable block, the preset normal coefficient and the preset tangential coefficient.

[0019] The preset equations of motion include:

[0020]

[0021] A i B i Let be the position vectors of the centroids of blocks A and B; Let be the translational velocity vectors of block A and block B; e is the angular velocity vector of block A and block B; ijk To arrange the tensors, the indices i, j, k take values ​​of 1, 2, and 3, respectively, representing the components of the vector or tensor in the global coordinate system; C i This serves as a reference point for the contact surface.

[0022] In one possible implementation, the formula for calculating the resultant force of the principal driving force on the target movable block includes:

[0023]

[0024] Among them, F ni Indicates the normal contact force i, Let i represent the tangential contact force.

[0025] In one possible implementation, the formula for calculating the preset safety factor corresponding to the double-sided sliding mode includes:

[0026]

[0027] Among them, F s This indicates the safety factor for the double-sided sliding mode.

[0028] In one possible implementation, the formula for calculating the preset safety factor corresponding to the edge rotation mode includes:

[0029]

[0030] Among them, M fs This indicates the safety factor for the edge rotation mode.

[0031] In one possible implementation, stability analysis results of the target rock mass are generated based on the safety factor of the target movable block, including:

[0032] When M fs <[M fs And F s >[F s The target movable block rotates around the rotation axis;

[0033] When M fs <[M fs And F s <[F s The target movable block slides along the structural surface i while rotating.

[0034] Secondly, embodiments of the present invention provide a rock mass stability analysis device, comprising:

[0035] The search module is used to locate all movable blocks within the target rock mass.

[0036] The first calculation module is used to calculate the contact forces of each structural surface of the target movable block; wherein, the target movable block is any one of the movable blocks, and the contact forces include normal contact forces and tangential contact forces;

[0037] The second calculation module is used to calculate the resultant force of the main driving force on the target movable block based on the gravity of the target movable block and the contact force of each structural surface of the target movable block;

[0038] The third calculation module is used to calculate the safety factor of the target movable block based on the resultant force of the active force and the preset safety factor calculation formula corresponding to multiple motion modes; among which, the multiple motion modes include at least the double-sided sliding mode and the edge rotation mode;

[0039] The analysis module is used to generate stability analysis results for the target rock mass based on the safety factor of the target movable block.

[0040] In one possible implementation, the first computing module is specifically used for:

[0041] Displacement constraints are applied to all movable blocks on free faces within the target rock mass.

[0042] Calculate the contact forces of each structural surface of the target movable block based on all free-face movable blocks under displacement constraints.

[0043] In one possible implementation, the first computing module is further specifically used for:

[0044] Calculate the contact velocity of each structural surface of the target movable block according to the preset motion equation;

[0045] Calculate the relative displacement increment of each structural surface of the target movable block based on the contact velocity of each structural surface;

[0046] The contact force of each structural surface of the target movable block is calculated based on the relative displacement increment of each structural surface of the target movable block, the area of ​​each structural surface contact point of the target movable block, the preset normal coefficient and the preset tangential coefficient.

[0047] The preset equations of motion include:

[0048]

[0049] A k B k Let be the position vectors of the centroids of blocks A and B; Let be the translational velocity vectors of block A and block B; e is the angular velocity vector of block A and block B; ijk For permutation tensors. The indices i, j, k take values ​​from 1 to 3 and represent the components of the vector or tensor in the global coordinate system; C k This serves as a reference point for the contact surface.

[0050] In one possible implementation, the formula for calculating the resultant force of the principal driving force on the target movable block includes:

[0051]

[0052] Among them, Fni Indicates the normal contact force i, Let i represent the tangential contact force.

[0053] In one possible implementation, the formula for calculating the preset safety factor corresponding to the double-sided sliding mode includes:

[0054]

[0055] Among them, F s This indicates the safety factor for the double-sided sliding mode.

[0056] In one possible implementation, the formula for calculating the preset safety factor corresponding to the edge rotation mode includes:

[0057]

[0058] Among them, M fs This indicates the safety factor for the edge rotation mode.

[0059] In one possible implementation, the analysis module is specifically used for:

[0060] When M fs <[M fs And F s >[F s The target movable block rotates around the rotation axis;

[0061] When M fs <[M fs And F s <[F s The target movable block slides along the structural surface i while rotating.

[0062] Thirdly, embodiments of the present invention provide a terminal, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the rock mass stability analysis method as described in the first aspect or any possible implementation thereof.

[0063] Fourthly, embodiments of the present invention provide a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the rock mass stability analysis method as described in the first aspect or any possible implementation thereof.

[0064] This invention provides a rock mass stability analysis method, apparatus, terminal, and storage medium. The rock mass stability analysis method is based on the critical block theory and considers the influence of inter-block forces on the stability of fractured rock mass blocks. Combining the critical block theory and the discrete element method, it analyzes and evaluates the stability of fractured rock mass blocks under the condition of considering inter-block forces, demonstrating good computational accuracy and applicability. Furthermore, based on considering the influence of inter-block forces on the stability of fractured rock mass, it further considers the block rotational motion modes, analyzing and calculating the net sliding force and safety factor for five motion modes (translation and rotation), providing a more comprehensive analysis of the influence of inter-block forces on the stability of fractured rock mass. Attached Figure Description

[0065] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the embodiments of the present invention will be briefly introduced below. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0066] Figure 1 This is a schematic flowchart of the rock mass stability analysis method provided in the embodiments of the present invention;

[0067] Figure 2 This is a schematic diagram of the rock mass stability analysis process provided in the embodiments of the present invention.

[0068] Figure 3 This is a schematic diagram of the tunnel cutting model provided in an embodiment of the present invention;

[0069] Figure 4 This is a schematic diagram of the inter-block forces provided in an embodiment of the present invention;

[0070] Figure 5 This is a schematic diagram of the block motion mode provided in an embodiment of the present invention;

[0071] Figure 6 This is a schematic diagram of the force vector between blocks provided in an embodiment of the present invention;

[0072] Figure 7 This is a diagram showing the analysis results of key tunnel blocks provided in an embodiment of the present invention;

[0073] Figure 8 This is a schematic diagram of the rock mass stability analysis device provided in an embodiment of the present invention;

[0074] Figure 9 This is a schematic diagram of the terminal provided in an embodiment of the present invention. Detailed Implementation

[0075] The features and exemplary embodiments of various aspects of the present invention will now be described in detail. To make the objectives, technical solutions, and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are only configured to explain the present invention and are not configured to limit the present invention. For those skilled in the art, the present invention can be practiced without some of these specific details. The following description of the embodiments is merely intended to provide a better understanding of the present invention by illustrating examples of the invention.

[0076] It should be noted that, in this document, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising..." does not exclude the presence of additional identical elements in the process, method, article, or apparatus that includes said element.

[0077] To address the problems of existing technologies, embodiments of the present invention provide a method, apparatus, terminal, and storage medium for rock mass stability analysis. The rock mass stability analysis method provided by these embodiments is described below.

[0078] The main body executing the rock mass stability analysis method can be a rock mass stability analysis device, which can be a terminal with data processing capabilities, such as a server, network attached storage (NAS), or personal computer (PC). This embodiment of the invention does not make specific limitations.

[0079] Figure 1 A flowchart illustrating the implementation of a rock mass stability analysis method provided in this embodiment of the invention is detailed below:

[0080] Step 110: Locate all movable blocks in the target rock mass.

[0081] In some embodiments, the target rock mass can be a rock slope in any region. A movable block refers to a potential critical block, that is, a block that has the potential to become a critical block.

[0082] Step 120: Calculate the contact forces of each structural surface of the target movable block.

[0083] In some embodiments, the target movable block is any one of all movable blocks, and the sub-contact forces include normal sub-contact forces and tangential sub-contact forces. The stability analysis of the target rock mass can be completed by processing all movable blocks according to the following procedure for processing the target movable block.

[0084] Optionally, the specific processing in step 120 can be as follows: apply displacement constraints to all free-face movable blocks in the target rock mass; and calculate the contact force of each structural surface of the target movable block based on all free-face movable blocks under displacement constraints.

[0085] In some embodiments, the specific processing for calculating the contact force of each structural surface of the target movable block based on all free-face movable blocks under displacement constraints can be as follows: calculate the contact velocity of each structural surface of the target movable block according to the preset motion equation; calculate the relative displacement increment of each structural surface of the target movable block according to the contact velocity of each structural surface of the target movable block; calculate the contact force of each structural surface of the target movable block according to the relative displacement increment of each structural surface of the target movable block, the area of ​​each structural surface contact point of the target movable block, the preset normal coefficient and the preset tangential coefficient.

[0086] In some embodiments, the preset equation of motion includes:

[0087]

[0088] Among them, A i B i Let be the position vectors of the centroids of blocks A and B; Let be the translational velocity vectors of block A and block B; e is the angular velocity vector of block A and block B; ijk For permutation tensors. The indices i, j, k take values ​​from 1 to 3 and represent the components of the vector or tensor in the global coordinate system; C i This serves as a reference point for the contact surface.

[0089] Step 130: Calculate the resultant force of the main driving force on the target movable block based on the gravity of the target movable block and the contact force of each structural surface of the target movable block.

[0090] In some embodiments, the formula for calculating the resultant force of the principal driving force on the target movable block may include:

[0091]

[0092] Among them, F ni Indicates the normal contact force i, Let i represent the tangential contact force.

[0093] Step 140: Calculate the safety factor of the target movable block according to the preset safety factor calculation formula corresponding to the resultant force of the active force and various motion modes.

[0094] In some embodiments, the various motion modes may include a free-fall mode, a single-sided sliding mode, a double-sided sliding mode, an edge rotation mode, and a sliding rotation mode. Generally, when the sliding mode of the block is the free-fall mode, its safety factor is considered to be 0.

[0095] In some embodiments, the formula for calculating the preset safety factor corresponding to the single-sided sliding mode may include:

[0096]

[0097] In some embodiments, the formula for calculating the preset safety factor corresponding to the double-sided sliding mode may include:

[0098]

[0099] Among them, F s This indicates the safety factor for the double-sided sliding mode.

[0100] In some embodiments, the formula for calculating the preset safety factor corresponding to the edge rotation mode may include:

[0101]

[0102] Among them, M fs This indicates the safety factor for the edge rotation mode.

[0103] Step 150: Based on the safety factor of the target movable block, generate the stability analysis results of the target rock mass.

[0104] In some embodiments, after obtaining the safety factor of the target movable block, the relationship between the safety factor and a corresponding threshold, such as a critical safety factor, can be used to determine whether the target movable block is a critical block and the rotation status of the target movable block. Specifically, when M fs <[M fs And F s >[F s The target movable block rotates around the rotation axis; when M fs <[M fs And F s <[F s The target movable block slides along the structural surface i while rotating.

[0105] In this way, by following the above processing procedure for the target movable block, the safety factor of all movable blocks of the target fractured rock mass can be obtained, and then the rotation of each movable block can be obtained. Thus, by combining the rotation of all possible blocks, the stability analysis result of the target rock mass can be generated.

[0106] For a better understanding of the above rock mass stability analysis method, please refer to [link / reference]. Figure 2 This diagram illustrates a rock mass stability analysis process.

[0107] Step 210: Establish a three-dimensional fractured rock mass model.

[0108] First, investigate the engineering background, establish a rock mass model, and import it into discrete element software, such as 3DEC Rhino software, to generate a fractured rock mass model by cutting the rock mass model according to the rock mass fracture information.

[0109] Taking an underground tunnel as an example, assuming the model is 41 meters wide and 65 meters high, containing a horseshoe-shaped chamber 65 meters long, and the fractured rock mass contains 8 sets of structural planes, their occurrence information and distribution patterns are as follows: Figure 3 As shown in (a) of the model. The specific generation process of this model is as follows: (1) Use 3DEC to establish a rock mass geometric model (without joint information); (2) Import the rock mass geometric model into 3DEC, input joint information to perform block cutting, and generate a fractured rock mass block model. Figure 3 As shown in (b), the model has a total of 44,723 blocks, including blocks formed by 3DEC virtual joint cutting.

[0110] Step 220, Block Mobility Analysis: Based on the key block theory, call the block information of the 3DEC fractured rock mass model to analyze the mobility of each block on the free surface and find all movable blocks on the free surface.

[0111] The basic theory of critical block mobility analysis is as follows:

[0112] (1) There are n sets of structural surfaces, each set of structural surfaces P i The attitude, i.e., the dip angle α i and tendency β i Find P i Upward unit normal vector

[0113]

[0114] (2) Find the intersection vector I of each structural surface. ij :

[0115]

[0116] (3) Introduce directional parameters for all possible intersection combinations of structural surfaces. Calculate its edge e ij Is it a true edge of the pyramid, where ij is the structure plane ij? The edges formed by the intersection of structure planes constitute the pyramid:

[0117]

[0118] Specifically, the directional parameters form a "directional parameter matrix" [I], and the sign numbers of each JP are listed in a diagonal matrix [D]:

[0119]

[0120] Where I(a1) is the symbol label of block structure surface 1, -1 represents the lower half space of the block on the plane, +1 represents the upper half space of the block on the plane, 0 represents the plane is not the interface of the block, and ±1 represents a pair of parallel interfaces of the plane that form the block.

[0121] (4) Establish the block discrimination matrix [T] = [I]·[D]

[0122] Specifically, when [T] corresponds to I ij When all elements ① in a certain row are "0" or contain both "+1" and "-1", I ij It is not the edge of BP; ② When both are "0" and "+1" I ij For true edges; ③ When both are "0" and "-1" -I ij . represents a real edge;

[0123] Thus, by examining the [T] matrix, if any one edge is a true edge of the block cone BP, then BP is a non-empty set, meaning its corresponding block is infinitely large. If none of the edges are true edges of BP, then BP is an empty set, meaning its corresponding block is finite in size.

[0124] In summary, firstly, the unit normal vector of the free face is calculated. Then, the intersection lines of each set of structural faces (including the intersection lines with the free face), i.e., the vectors of the edges, are calculated. Next, a discrimination matrix [T] is established to determine whether the intersection lines of each set of structural faces (including the intersection lines with the free face) are true edges of the pyramid. Finally, based on the necessary and sufficient conditions for the block to be movable, i.e., JP (fissure cone) ≠ Ф; BP (block cone) = JP ∩ EP (excavation cone) = Ф, all movable blocks can be identified.

[0125] Specifically, displacement boundaries are applied to the model, and the effects of ground stress are ignored. The mobility of the tunnel free face blocks is analyzed. The structural surface information in the free face block group information file is read to determine whether the fracture cone JP is a non-hollow fracture cone. If the fracture cone is a non-empty fracture cone, then read the block structure surface and free surface information, that is, all the constituent surfaces of the free surface block group, and determine...

[0126] Step 230, Calculation of inter-block forces.

[0127] The principle for calculating the forces between blocks is as follows:

[0128] After the mobility analysis is completed, the normal and shear force vectors of all sub-contact points on the block can be calculated using the displacement constraint method proposed based on 3DEC. The specific processing flow of this method is as follows: (1) Perform displacement constraints on all free-faced movable blocks; (2) Limit their displacement to enable the 3DEC iterative calculation to converge quickly; (3) Read the sub-contact forces of each structural surface of the movable block, i.e., the forces between the blocks, such as Figure 4 As shown. The calculation principle is as follows:

[0129] The unit normal perpendicular to the contact surface is considered the contact normal (pointing from block A to block B), and it is assumed that the relative velocity through the sub-contact is determined by the velocity V associated with the sub-contact. i V and the velocity V of the corresponding point on the opposite surface i F Obtained. For a rigid block, the velocity is calculated from the equation of motion of the rigid block. The equation of translational motion of the rigid block is:

[0130]

[0131] The equation of motion for the rigid block is:

[0132]

[0133] in:

[0134] In the above formula, α is the velocity of the block's center of mass; α is the viscous damping coefficient, F i The sum of forces acting on the block, where m is the mass of the block; g i I1, I2, I3 are the principal moments of inertia of the block; ω1, ω2, ω3 are the angular velocities about the principal axis of the block; M1, M2, M3 are the torque components applied to the principal axis.

[0135] For a rigid block, the sub-contact velocity (defined as the velocity of block B relative to block A at the sub-contact position) is calculated as follows:

[0136]

[0137] In the above formula, A k B kThe position vectors of the centroids of blocks A and B; Let be the translational velocity vectors of block A and block B; e is the angular velocity vector of block A and block B; ijk For permuting tensors, the subscripts i, j, k range from 1 to 3, and represent the components of the vector or tensor in the global coordinate system; C k This serves as a reference point for the contact surface.

[0138] The velocity of the sub-contact associated with the rigid block contact is interpolated from the contact velocity, and the relative displacement increment at the rigid block sub-contact is given by the following equation:

[0139] ΔU i =V i Δt (8)

[0140] It can be decomposed into normal and tangential components along the contact surface, where:

[0141] The normal displacement increment is:

[0142] ΔU n =ΔU i n i (9)

[0143] The tangential displacement increment vector is:

[0144]

[0145] The sub-contact displacement increment is used to calculate the tangential force increment and the normal force increment. Taking the compressive force as positive, the normal force increment is:

[0146] ΔF n =-K n ΔU n A c (11)

[0147] The increment of the tangential force vector is:

[0148]

[0149] Where A c K represents the area of ​​the sub-contact. n K is the normal stiffness of the contact surface. s This represents the shear stiffness of the contact surface.

[0150] The total normal and tangential force vectors of the sub-contacts are updated as follows:

[0151] F ni =F ni +ΔF n (13)

[0152] Fsi =F si +ΔF i s (14)

[0153] The sub-contact force vector represents the action of block A on block B, given by equation (12):

[0154]

[0155] F n For the normal contact force i, Let i be the tangential contact force vector.

[0156] Step 240, Key Block Analysis:

[0157] According to formula (16), the calculated inter-block forces can be used to correct the active resultant force vector and evaluate the stability of the movable block.

[0158] The resultant force vector acting on each block can be calculated using the following equation:

[0159]

[0160] F ni Normal contact force i, Let i be the tangential sub-contact force vector, n refers to the normal direction of the block, s refers to the tangential direction, and i is the sub-contact force number.

[0161] The kinematic conditions that each motion mode of the block must satisfy are as follows: Figure 5 As shown, the net sliding force and net rotational torque are:

[0162]

[0163]

[0164]

[0165] in, The resultant force is the active force; F is the net sliding force. Let c be the internal friction angle of structural surface i; i A represents the cohesion of structural plane i; i Let i be the area of ​​the structural surface.

[0166] If the block intersects the edge line I fi If the block rotates, the forces acting on it will affect I. fi The moment at any endpoint o is

[0167]

[0168] Where, r ir is the radius vector from point o to the point of application of the contact force vector; g r is the radius vector from point o to the center of mass of the block; c Let be the radius vector from point o to the centroid of structural surface i. For example... Figure 6 As shown:

[0169] During edge rotation, the axis of rotation is the intersection of the free surface and the structural surface that satisfy the motion conditions. The block rotates around the axis of rotation I. fi The net torque is:

[0170] M fi =M o ·n fi (twenty one)

[0171] n fi =n f ×n i / |n f ×n i | (22)

[0172] Where point O is any endpoint of the intersection line between the free surface and the structural surface that satisfies the motion conditions, and Mo is the force exerted on the block on I. fi The moment n at any endpoint o fi The line of intersection between structural plane i and free plane f, n i It is the upward unit normal vector of the structural surface, n f It is the upward unit normal vector of the free plane.

[0173] The safety factor for each motion mode is (the safety factor for block detachment is 0):

[0174] Slide along one side:

[0175]

[0176] Slide along both sides:

[0177]

[0178] Edge rotation:

[0179]

[0180] When M fs <[M fs And F s >[F s The block rotates around the rotation axis.

[0181] Sliding rotation:

[0182] The block simultaneously satisfies the conditions for sliding and rotational motion, when M fs <[M fs And Fs <[F s The block slides along structural surface i while rotating.

[0183] In the formula, The driving force is the net sliding force; F is the net sliding force. Let c be the internal friction angle of structural surface i; i A represents the cohesion of structural plane i; i Let r be the area of ​​structural surface i. i r is the radius vector from point o to the point of application of the contact force vector; g r is the radius vector from point O to the center of mass of the block; c Let be the radius vector from point o to the centroid of the structural surface i.

[0184] Taking the aforementioned underground tunnel as an example, after the mobility analysis is completed and the displacement of all movable blocks on the free face is restricted, iterative calculations are performed. The calculation parameters for the tunnel surrounding rock are: density ρ = 2680 kg / m³. 3 The normal and tangential stiffness between the blocks are both 20 GPa / m, the internal friction angle is 30°, and the cohesion is 0.18 MPa.

[0185] After the iteration is completed, the normal sub-contact force vector and tangential sub-contact force vector of all sub-contacts on the block group are obtained, which means the inter-block forces acting on the movable block.

[0186] Based on formula (1), the driving resultant force is corrected by the calculated inter-block interaction force, and the motion mode of the movable block (free fall, single-sided sliding, double-sided sliding, edge rotation, sliding rotation) is determined.

[0187] Define the critical safety factor value of the block (the block below this value is the critical block). The safety factor of the free-falling block is 0. The safety factors of the block with single-sided sliding, double-sided sliding and edge rotation are respectively determined by formulas (23), (24) and (25). After the sliding and rotation conditions are met, the safety factor is the minimum value of the sliding safety factor and the rotation safety factor.

[0188] Without considering the inter-block forces, the calculated key blocks of the tunnel free face are as follows: Figure 7 As shown in (a), considering the inter-block forces, the key blocks of the tunnel free face are as follows: Figure 7 As shown in (b).

[0189] In this embodiment of the invention, a rock mass stability method considering inter-block forces is provided. This method is based on the critical block theory and considers the influence of inter-block forces on the stability of fractured rock mass blocks. Combining the critical block theory and the discrete element method, the stability of fractured rock mass blocks considering inter-block forces is analyzed and evaluated, exhibiting good computational accuracy and applicability.

[0190] Furthermore, based on the influence of inter-block forces on the stability of fractured rock masses, the rotational motion modes of the blocks are further considered, and the net sliding force and safety factor of five motion modes (translation and rotation) are analyzed and calculated, providing a more comprehensive analysis of the influence of inter-block forces on the stability of fractured rock masses.

[0191] In addition, a method for analyzing the stability of fractured rock masses based on the critical block theory, which considers the inter-block forces and the blocks themselves, can provide theoretical support for the stability study of fractured rock masses under complex loading conditions. It also has the advantages of simple operation, accurate and efficient calculation, and good applicability.

[0192] It should be understood that the sequence number of each step in the above embodiments does not imply the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.

[0193] The following are device embodiments of the present invention. For details not described in detail, please refer to the corresponding method embodiments described above.

[0194] Figure 8 A schematic diagram of the rock mass stability analysis device provided in an embodiment of the present invention is shown. For ease of explanation, only the parts related to the embodiment of the present invention are shown, and are described in detail below:

[0195] like Figure 8 As shown, the rock mass stability analysis device includes:

[0196] The search module 810 is used to find all movable blocks in the target rock mass;

[0197] The first calculation module 820 is used to calculate the contact forces of each structural surface of the target movable block; wherein, the target movable block is any one of the movable blocks, and the contact forces include normal contact forces and tangential contact forces;

[0198] The second calculation module 830 is used to calculate the resultant force of the main driving force on the target movable block based on the gravity of the target movable block and the contact force of each structural surface of the target movable block;

[0199] The third calculation module 840 is used to calculate the safety factor of the target movable block according to the resultant force of the active force and the preset safety factor calculation formula corresponding to multiple motion modes; wherein, the multiple motion modes include at least the double-sided sliding mode and the edge rotation mode.

[0200] Analysis module 850 is used to generate stability analysis results of the target rock mass based on the safety factor of the target movable block.

[0201] In one possible implementation, the first computing module is specifically used for:

[0202] Displacement constraints are applied to all movable blocks on free faces within the target rock mass.

[0203] Calculate the contact forces of each structural surface of the target movable block based on all free-face movable blocks under displacement constraints.

[0204] In one possible implementation, the first computing module is further specifically used for:

[0205] Calculate the contact velocity of each structural surface of the target movable block according to the preset motion equation;

[0206] Calculate the relative displacement increment of each structural surface of the target movable block based on the contact velocity of each structural surface;

[0207] The contact force of each structural surface of the target movable block is calculated based on the relative displacement increment of each structural surface of the target movable block, the area of ​​each structural surface contact point of the target movable block, the preset normal coefficient and the preset tangential coefficient.

[0208] The preset equations of motion include:

[0209]

[0210] A i B i Let be the position vectors of the centroids of blocks A and B; Let be the translational velocity vectors of block A and block B; e is the angular velocity vector of block A and block B; ijk For permutation tensors. The indices i, j, k take values ​​from 1 to 3 and represent the components of the vector or tensor in the global coordinate system; C i This serves as a reference point for the contact surface.

[0211] In one possible implementation, the formula for calculating the resultant force of the principal driving force on the target movable block includes:

[0212]

[0213] Among them, F ni Indicates the normal contact force i, Let i represent the tangential contact force.

[0214] In one possible implementation, the formula for calculating the preset safety factor corresponding to the double-sided sliding mode includes:

[0215]

[0216] Among them, F s This indicates the safety factor for the double-sided sliding mode.

[0217] In one possible implementation, the formula for calculating the preset safety factor corresponding to the edge rotation mode includes:

[0218]

[0219] Among them, M fs This indicates the safety factor for the edge rotation mode.

[0220] In one possible implementation, the analysis module is specifically used for:

[0221] When M fs <[M fs And F s >[F s The target movable block rotates around the rotation axis;

[0222] When M fs <[M fs And F s <[F s The target movable block slides along the structural surface i while rotating.

[0223] Figure 8 The provided rock mass stability analysis device has various modules that can achieve... Figure 1 The functions of each step in the illustrated embodiment are as follows, and they achieve the same results as... Figure 1 The same technical effect as the rock mass stability analysis method shown is not described in detail here for the sake of brevity.

[0224] Figure 9 A schematic diagram of the hardware structure of the terminal provided in an embodiment of the present invention is shown.

[0225] The terminal may include a processor 901 and a memory 902 storing computer program instructions.

[0226] Specifically, the processor 901 may include a central processing unit (CPU), an application specific integrated circuit (ASIC), or one or more integrated circuits that can be configured to implement the embodiments of the present invention.

[0227] Memory 902 may include a large-capacity memory for data or instructions. For example, and not limitingly, memory 902 may include a hard disk drive (HDD), a floppy disk drive, flash memory, optical disk, magneto-optical disk, magnetic tape, or a Universal Serial Bus (USB) drive, or a combination of two or more of these. Where appropriate, memory 902 may include removable or non-removable (or fixed) media. Where appropriate, memory 902 may be internal or external to the integrated gateway disaster recovery device. In a particular embodiment, memory 902 is a non-volatile solid-state memory. In a particular embodiment, memory 902 includes read-only memory (ROM). Where appropriate, the ROM may be a mask-programmed ROM, a programmable ROM (PROM), an erasable PROM (EPROM), an electrically erasable PROM (EEPROM), an electrically rewritable ROM (EAROM), or flash memory, or a combination of two or more of these.

[0228] The processor 901 reads and executes computer program instructions stored in the memory 902 to implement any of the rock mass stability analysis methods in the above embodiments.

[0229] In one example, the terminal may also include a communication interface 903 and a bus 910. For example, Figure 9 As shown, the processor 901, memory 902, and communication interface 903 are connected through bus 910 and complete communication with each other.

[0230] The communication interface 903 is mainly used to realize communication between various modules, devices, units and / or equipment in the embodiments of the present invention.

[0231] Bus 910 includes hardware, software, or both, that couples terminal components together. For example, and not limitingly, the bus may include an Accelerated Graphics Port (AGP) or other graphics bus, an Enhanced Industry Standard Architecture (EISA) bus, a Front Side Bus (FSB), HyperTransport (HT) interconnect, an Industry Standard Architecture (ISA) bus, an Infinite Bandwidth Interconnect, a Low Pin Count (LPC) bus, a memory bus, a Microchannel Architecture (MCA) bus, a Peripheral Component Interconnect (PCI) bus, a PCI-Express (PCI-X) bus, a Serial Advanced Technology Attachment (SATA) bus, a Video Electronics Standards Association Local (VLB) bus, or other suitable buses, or combinations of two or more of these. Where appropriate, bus 910 may include one or more buses. While specific buses are described and illustrated in embodiments of the invention, the invention contemplates any suitable bus or interconnect.

[0232] This terminal can execute the rock mass stability analysis method in this embodiment of the invention, thereby achieving the combination of Figure 1 and... Figure 8 The described rock mass stability analysis method and apparatus.

[0233] Furthermore, in conjunction with the rock mass stability analysis methods described in the above embodiments, this invention can be implemented using a computer storage medium. This computer storage medium stores computer program instructions; when these computer program instructions are executed by a processor, they implement any of the rock mass stability analysis methods described in the above embodiments.

[0234] It should be clarified that the present invention is not limited to the specific configurations and processes described above and shown in the figures. For the sake of brevity, detailed descriptions of known methods are omitted here. In the above embodiments, several specific steps are described and shown as examples. However, the method process of the present invention is not limited to the specific steps described and shown. Those skilled in the art can make various changes, modifications, and additions, or change the order of steps, after understanding the spirit of the present invention.

[0235] The functional blocks shown in the above-described structural diagram can be implemented as hardware, software, firmware, or a combination thereof. When implemented in hardware, they can be, for example, electronic circuits, application-specific integrated circuits (ASICs), appropriate firmware, plug-ins, function cards, etc. When implemented in software, the elements of this invention are programs or code segments used to perform the required tasks. The programs or code segments can be stored on a machine-readable medium or transmitted over a transmission medium or communication link via data signals carried in a carrier wave. "Machine-readable medium" can include any medium capable of storing or transmitting information. Examples of machine-readable media include electronic circuits, semiconductor memory devices, ROM, flash memory, erasable ROM (EROM), floppy disks, CD-ROMs, optical disks, hard disks, fiber optic media, radio frequency (RF) links, etc. Code segments can be downloaded via computer networks such as the Internet, intranets, etc.

[0236] It should also be noted that the exemplary embodiments mentioned in this invention describe methods or systems based on a series of steps or apparatus. However, this invention is not limited to the order of the steps described above; that is, the steps can be performed in the order mentioned in the embodiments, or in a different order, or several steps can be performed simultaneously.

[0237] The above description is merely a specific embodiment of the present invention. Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working processes of the systems, modules, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here. It should be understood that the protection scope of the present invention is not limited thereto. Any person skilled in the art can easily conceive of various equivalent modifications or substitutions within the technical scope disclosed in the present invention, and these modifications or substitutions should all be covered within the protection scope of the present invention.

Claims

1. A method for analyzing rock mass stability, characterized in that, include: Locate all movable blocks within the target rock mass; Calculate the contact forces of each structural surface of the target movable block; wherein, the target movable block is any one of the movable blocks, and the contact forces include normal contact forces and tangential contact forces; Calculate the resultant force of the main driving force on the target movable block based on the gravity of the target movable block and the contact force of each structural surface of the target movable block; The safety factor of the target movable block is calculated according to the preset safety factor calculation formula corresponding to the resultant force of the active force and the multiple motion modes; wherein, the multiple motion modes include at least the double-sided sliding mode and the edge rotation mode. Based on the safety factor of the target movable block, the stability analysis results of the target rock mass are generated; The calculation of the contact forces of each structural surface of the target movable block includes: Displacement constraints are applied to all free-face movable blocks in the target rock mass; based on all free-face movable blocks under displacement constraints, the contact forces of each structural surface of the target movable block are calculated; Among them, the calculation of the contact forces of each structural surface of the target movable block based on displacement constraints for all free-face movable blocks includes: The contact velocity of each structural surface of the target movable block is calculated according to the preset motion equation; Based on the contact velocity of each structural surface of the target movable block, calculate the relative displacement increment of each structural surface of the target movable block; The contact force of each structural surface of the target movable block is calculated based on the relative displacement increment of each structural surface of the target movable block, the area of ​​each structural surface contact point of the target movable block, the preset normal coefficient and the preset tangential coefficient. The preset motion equations include: Ai,Bi are the position vectors of the centroids of blocks A and B; Let be the translational velocity vectors of block A and block B; e is the angular velocity vector of block A and block B; ijk For arranging tensors, the indices i, j, k range from 1 to 3 and represent the components of the vector or tensor in the global coordinate system; Ck is the reference point of the contact surface.

2. The rock mass stability analysis method according to claim 1, characterized in that, The formula for calculating the resultant force of the main driving force on the target movable block includes: Among them, F ni Indicates the normal contact force i, Let i represent the tangential contact force.

3. The rock mass stability analysis method according to claim 2, characterized in that, The preset safety factor calculation formula for the double-sided sliding mode includes: Among them, F s Indicates the safety factor for the double-sided sliding mode; Let c be the internal friction angle of structural surface i. i The cohesive force of structural plane i; The internal friction angle of structural surface j; c j Let be the cohesive force of structural plane j.

4. The rock mass stability analysis method according to claim 3, characterized in that, The formula for calculating the preset safety factor for the edge rotation mode includes: n fi =n f ×n i / |n f ×n i |; Among them, M fs Indicates the safety factor for edge rotation mode; r i r is the radius vector from point o to the point of application of the contact force vector; g r is the radius vector from point O to the center of mass of the block; c v is the radius vector from point o to the centroid of structural plane i, where point o is any endpoint of the intersection line between the free plane and the structural plane; i Let n be the normal vector pointing from structural plane i to the interior of the block; i It is the upward unit normal vector of structure surface i, n f It is the upward unit normal vector of the free plane f.

5. The rock mass stability analysis method according to claim 4, characterized in that, The step of generating stability analysis results for the target rock mass based on the safety factor of the target movable block includes: When Mfs < [Mfs] and Fs > [Fs], the target movable block rotates about the rotation axis; When Mfs < [Mfs] and Fs < [Fs], the target movable block slides along structural surface i while rotating; [Mfs] represents the safety factor for the edge rotation mode allowed by the project; [Fs] represents the safety factor for the double-sided sliding mode allowed by the project.

6. A rock mass stability analysis device, characterized in that, include: The search module is used to locate all movable blocks within the target rock mass. The first calculation module is used to calculate the contact forces of each structural surface of the target movable block; wherein, the target movable block is any one of the movable blocks, and the contact forces include normal contact forces and tangential contact forces; The second calculation module is used to calculate the resultant force of the main driving force on the target movable block based on the gravity of the target movable block and the contact force of each structural surface of the target movable block; The third calculation module is used to calculate the safety factor of the target movable block according to the resultant force of the active force and the preset safety factor calculation formula corresponding to the multiple motion modes; wherein, the multiple motion modes include at least the double-sided sliding mode and the edge rotation mode; The analysis module is used to generate stability analysis results for the target rock mass based on the safety factor of the target movable block. The first calculation module is specifically used for: Displacement constraints are applied to all free-face movable blocks in the target rock mass; based on all free-face movable blocks under displacement constraints, the contact forces of each structural surface of the target movable block are calculated; Among them, the calculation of the contact forces of each structural surface of the target movable block based on displacement constraints for all free-face movable blocks includes: The contact velocity of each structural surface of the target movable block is calculated according to the preset motion equation; Based on the contact velocity of each structural surface of the target movable block, calculate the relative displacement increment of each structural surface of the target movable block; The contact force of each structural surface of the target movable block is calculated based on the relative displacement increment of each structural surface of the target movable block, the area of ​​each structural surface contact point of the target movable block, the preset normal coefficient and the preset tangential coefficient. The preset motion equations include: Ai,Bi are the position vectors of the centroids of blocks A and B; Let be the translational velocity vectors of block A and block B; e is the angular velocity vector of block A and block B; ijk For arranging tensors, the indices i, j, k range from 1 to 3 and represent the components of the vector or tensor in the global coordinate system; Ck is the reference point of the contact surface.

7. A terminal, comprising a memory and a processor, wherein the memory is used to store a computer program, and the processor is used to call and run the computer program stored in the memory, characterized in that, When the processor executes the computer program, it implements the steps of the rock mass stability analysis method as described in any one of claims 1 to 5.

8. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the steps of the rock mass stability analysis method as described in any one of claims 1 to 5.

Images