Transverse isotropic rock mass discrete element model parameter checking method and system

Abstract

The invention relates to the technical field of coal mining, and provides a transverse isotropic rock mass discrete element model parameter checking method and system. According to the analysis method, a complete matrix model Matrix of a rock mass is generated based on a parallel combination model, the matrix modulus of the matrix model Matrix is calibrated by adjusting the bonding rigidity of bonding particles of the matrix model Matrix, and the compressive strength of the matrix model Matrix is calibrated by adjusting the strength parameter of the bonding particles; inserting a smooth contact interface for simulating the contact between a bedding surface and a weak surface of the layer rock mass into the matrix model Materix based on the smooth contact model, and generating a rock mass model Mass of the layer rock mass; when the dip angle of the bedding surface of the rock mass model Mass is (0 degree, 90 degrees), the compressive strength of the rock mass model Mass is adjusted by adjusting the strength parameter of the smooth contact interface, so that the calibration process of the layer rock mass model has repeatability.

Description

Technical Field

[0001] This application relates to the technical field of coal mining, and particularly relates to a method and system for calibrating parameters of a transversely isotropic rock mass discrete element model. Background Art

[0002] The application of the discrete element method (DEM) in geologic materials and rock mechanics has a rather long history. Especially in simulating particulate materials and rock crack propagation, it is widely used to simulate and study the mechanical behaviors of materials such as rocks and soils, among which the bonded particle model (BPM) is the most widely applied.

[0003] For the existing parallel bonded model (PBM) based on the discrete element method (DEM), especially the calibration of composite materials created by simultaneously adopting the constitutive model of the parallel bonded model and the smooth joint model (SJM), usually a large number of trial-and-error processes are required. The relationship between the microscopic parameters and macroscopic behaviors of the model is complex and non-linear. The calibration process requires repeatedly adjusting multiple parameters. Especially when the model involves material stiffness and strength in different directions, the number of trial-and-error times will increase significantly. Moreover, this trial-and-error method is complex and time-consuming on the one hand; on the other hand, it is difficult to be systematic. During the calibration process, the coupling relationship between parameters makes the adjustment of one parameter often affect other parameters, causing other parameters to change, increasing the complexity of calibration and unable to ensure repeatability and efficiency.

[0004] Therefore, there is an urgent need to provide a technical solution to address the deficiencies of the above-mentioned existing technologies. Summary of the Invention

[0005] The purpose of this application is to provide a method and system for calibrating parameters of a transversely isotropic rock mass discrete element model to solve or alleviate the problems existing in the above-mentioned existing technologies.

[0006] To achieve the above purpose, this application provides the following technical solutions:

[0007] This application provides a method for calibrating parameters of a transversely isotropic rock mass discrete element model, including: Step S101, generating a complete matrix model Matrix of the layered rock mass based on the parallel bonded model, calibrating the matrix modulus of the matrix model Matrix by adjusting the bond stiffness of the bonded particles of the matrix model Matrix, and calibrating the compressive strength of the matrix model Matrix by adjusting the strength parameters of the bonded particles;

[0008] Step S102: Insert a smooth contact interface for simulating the contact between the bedding plane and the weak plane of the layered rock mass into the matrix model Matrix based on the smooth contact model, and generate the rock mass model Ms of the layered rock mass;

[0009] Step S103: In response to the bedding plane inclination angle of the rock mass model Mass being (0°, 90°), adjust the compressive strength of the rock mass model Mass by adjusting the strength parameters of the smooth contact interface.

[0010] Preferably, in step S101, the calibration of the matrix modulus of the matrix model Matrix by adjusting the bond stiffness between the bonded particles in the matrix model Matrix is specifically: adjusting the normal stiffness k n,pbm and shear stiffness k s,pbm of the bonded particles in the matrix model Matrix to calibrate the effective modulus E pbm of the matrix model Matrix.

[0011] Preferably, in step S101, the calibration of the compressive strength of the matrix model Matrix by adjusting the strength parameters of the bonded particles is specifically: adjusting the tensile strength σ t,pbm , cohesion c pbm and friction angle μ pbm of the bonded particles in the matrix model Matrix to calibrate the compressive strength of the matrix model Matrix.

[0012] Preferably, in step S103, the adjustment of the compressive strength of the rock mass model Mass in response to the bedding plane inclination angle of the rock mass model Mass being (0°, 90°) by adjusting the strength parameters of the smooth contact interface is specifically: adjusting the normal stiffness k n,sj and tangential stiffness k s,sj of the smooth contact interface to adjust the matrix modulus E0 in the direction parallel to the bedding plane and the rock mass modulus E 90 in the direction perpendicular to the bedding plane of the rock mass model Mass;

[0013] Adjust the tensile strength σ t,sj , cohesion c sj and friction angle μ sj of the smooth contact interface to the uniaxial compressive strength in the rock mass model Mass.

[0014] Preferably, the adjustment of the normal stiffness k n,sj and tangential stiffness k s,sj, adjust the matrix modulus E0 in the direction parallel to the bedding plane and the rock mass modulus E in the direction perpendicular to the bedding plane in the rock mass model Mass, including: in response to the error between the matrix modulus E0 in the direction parallel to the bedding plane in the rock mass model Mass and the matrix modulus E0 in the direction parallel to the bedding plane in the layered rock mass being greater than a preset first threshold, adjust the tangential stiffness k 90 of the smooth contact interface; s,sj

[0015] In response to the error between the rock mass modulus E in the direction perpendicular to the bedding plane in the rock mass model Mass and the rock mass modulus E in the direction perpendicular to the bedding plane in the layered rock mass 90 being greater than a preset second threshold, adjust the normal stiffness k 90 of the smooth contact interface. n,sj

[0016] Preferably, when adjusting the tensile strength σ t,sj , cohesion c sj and friction angle μ sj of the smooth contact interface for the uniaxial compressive strength in the rock mass model Mass, including: in response to the error between the uniaxial compressive strength UCS0 in the direction parallel to the loading direction in the rock mass model Mass and the uniaxial compressive strength UCS0 in the direction parallel to the loading direction in the layered rock mass being greater than a preset third threshold, keep the friction angle μ sj of the smooth contact interface in the rock mass model Mass unchanged, and adjust the tensile strength σ t,si , cohesion c sj .

[0017] The embodiment of the present application also provides a transverse isotropic rock mass discrete element model parameter calibration system, including:

[0018] A matrix model and a matrix stiffness adjustment unit, configured to generate a matrix model Matrix of the layered rock mass based on the parallel bond model, calibrate the matrix modulus of the matrix model Matrix by adjusting the bond stiffness of the bonded particles in the matrix model Matrix, and calibrate the compressive strength of the matrix model Matrix by adjusting the strength parameters of the bonded particles;

[0019] A rock mass model generation unit, configured to generate a rock mass model Mass of the layered rock mass by inserting a smooth contact interface for simulating the contact between the bedding plane and the weak plane of the layered rock mass into the matrix model Matrix based on the smooth contact model;

[0020] The bedding plane stiffness adjustment unit is configured to adjust the compressive strength of the rock model Mass by adjusting the strength parameter of the smooth contact interface in response to the bedding plane inclination angle of the rock model Mass being (0°, 90°).

[0021] Beneficial effects:

[0022] In the parameter calibration method of the transversely isotropic rock mass discrete element model provided in the embodiment of the present application, first, a complete matrix model Matrix of the layered rock mass is generated based on the parallel bonding model, and the matrix modulus of the matrix model Matrix is calibrated by adjusting the bonding stiffness of the bonding particles of the matrix model Matrix, and the compressive strength of the matrix model Matrix is calibrated by adjusting the strength parameters of the bonding particles; then, a smooth contact interface simulating the contact between the bedding plane and the weak plane of the layered rock mass is inserted into the matrix model Matrix based on the smooth contact model to generate a rock mass model Mass of the layered rock mass; when the inclination angle of the bedding plane of the rock mass model Mass is (0°, 90°), the compressive strength of the rock mass model Mass is adjusted by adjusting the strength parameters of the smooth contact interface.

[0023] Therefore, in the process of layered rock mass simulation, the parallel bonding model and the smooth contact model are combined to divide the calibration process of the layered rock mass model parameters into multiple independent calibration steps, and in each independent calibration step, different parameters of the layered rock mass model are independently calibrated in sequence. The calibration process has multi-scale and multi-directional adaptability, so that the calibration of the layered rock mass model is no longer affected by the size of the model; at the same time, it effectively avoids the parameter coupling during the calibration of the layered rock mass model in the existing technology, which not only simplifies the parameter adjustment process of the model and improves the model calibration efficiency, but also makes the calibration process of the layered rock mass model repeatable. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] The drawings in the specification, which constitute a part of this application, are used to provide further understanding of this application. The illustrative embodiments of this application and their descriptions are used to explain this application and do not constitute improper limitations on this application.

[0025] in:

[0026] Figure 1 A schematic flow chart of a parameter calibration method for a transversely isotropic rock mass discrete element model provided according to some embodiments of the present application;

[0027] Figure 2 A schematic diagram of a mechanical constitutive model of a smooth contact surface in a transversely isotropic rock mass provided according to some embodiments of the present application;

[0028] Figure 3Schematic diagram of the smooth contact surface model in the transversely isotropic rock mass discrete element model provided according to some embodiments of the present application;

[0029] Figure 4 Schematic diagram of the structure of a parameter calibration system for a transversely isotropic rock mass discrete element model provided according to some embodiments of the present application. Specific embodiments

[0030] The present application will be described in detail below with reference to the accompanying drawings and in conjunction with embodiments. Each example is provided by way of explanation of the present application rather than limitation of the present application. In fact, those skilled in the art will appreciate that modifications and variations can be made to the present application without departing from the scope or spirit thereof. For example, features shown or described as part of one embodiment can be used in another embodiment to yield yet another embodiment. All other embodiments obtained by those of ordinary skill in the art based on the embodiments in the present invention shall fall within the scope of protection of the embodiments of the present invention.

[0031] The existing calibration process for combining particle (i.e., bonded particle) models in geological material simulation is only applicable to a single contact model (such as a parallel bond model or a soft contact model, etc.). However, transversely isotropic materials (such as shale) usually require combining multiple contact models to simulate complex mechanical behaviors, such as the strength and stiffness in the layer plane. And when the model simultaneously includes multiple contact models (such as a parallel bond model and a smooth contact model), the existing numerical simulation techniques can usually only be applied to a single contact model, and it is impossible to simultaneously consider the different mechanical properties of the layer plane and the collective material. Especially when dealing with transversely isotropic materials (such as shale, coal seam, etc.), it is difficult to accurately calibrate multiple parameters in the model when dealing with the isotropy of the material.

[0032] Based on this, the embodiments of the present application provide a method for calibrating parameters of a transversely isotropic rock mass discrete element model. By dividing the calibration process into multiple independent steps and adjusting different parameters respectively, it is ensured that the influence of parameter changes in each step on the parameters in other steps is minimized, simplifying the process of adjusting parameters of the model and improving the calibration efficiency. As Figures 1 to 3 shown, the method for calibrating parameters of the transversely isotropic rock mass discrete element model includes:

[0033] Step S101: Generate a complete matrix model Matrix of the layered rock mass based on the parallel bond model, calibrate the matrix model of the matrix model Matrix by adjusting the bond stiffness of the bonded particles in the matrix model Matrix, and calibrate the compressive strength of the matrix model Matrix by adjusting the strength parameters of the bonded particles.

[0034] In the parallel bond model, the simulated material of the layered rock mass is connected by bonded particles to simulate the formation and propagation of rock cracks. However, in the parallel bond model, only single contact can be simulated, and the characteristics of the material matrix and bedding plane cannot be simulated and calibrated simultaneously. Here, using the parallel bond model to simulate the layered rock mass, the elastic modulus and stiffness of the material matrix are simulated through the bonding relationship between particles (simulated materials), and a complete matrix model Matrix of the layered rock mass is generated. Then, the matrix model of the matrix model Matrix is calibrated by adjusting the bonding stiffness of the bonded particles in the matrix model Matrix, while retaining the advantages of the matrix stiffness in the parallel bond model.

[0035] Specifically, adjust the normal stiffness k of the bonded particles in the matrix model Matrix n,pbm and the shear stiffness k s,pbm , and calibrate the effective modulus E pbm of the matrix model Matrix, laying a foundation for the subsequent calibration of the bedding plane. At the same time, by adjusting the bonding particle strength parameters (tensile strength σ t,pbm , cohesion c pbm and friction angle μ pbm ), the compressive strength of the matrix model Matrix is calibrated. That is to say, by adjusting the tensile strength σ t,pbm , cohesion c pbm and friction angle μ pbm of the bonded particles, the uniaxial compressive strength of the matrix model Matrix is calibrated.

[0036] Among them, first set the tensile strength σ t,pbm , cohesion c pbm and friction angle μ pbm of the bonded particles as the laboratory theoretical data. On this basis, compare the compressive strength of the matrix model Matrix with the actual laboratory strength. If the compressive strength of the matrix model Matrix is greater than the actual laboratory strength, then according to the ratio of the compressive strength of the matrix model Matrix to the actual laboratory strength, the values of the tensile strength σ t,pbm , cohesion c pbm and friction angle μ pbm of the bonded particles are reduced in the same proportion; conversely, if the compressive strength of the matrix model Matrix is less than the actual laboratory strength, then according to the ratio of the compressive strength of the matrix model Matrix to the actual laboratory strength, the values of the tensile strength σ t,pbm , cohesion c pbm and friction angle μ pbm of the bonded particles are increased in the same proportion.

[0037] Step S102: Insert a smooth contact interface that simulates the contact between the bedding plane and the weak plane of the layered rock mass in the matrix model Matrix based on the smooth contact model, and generate the rock mass model Mass of the layered rock mass.

[0038] In this application, in the parallel bond model, the elastic model and strength of the material matrix are simulated through the bond relationship between the bonded particles. The adjustment of the matrix modulus is achieved through the parallel bond model, and the advantage of the matrix stiffness in the parallel bond model is retained. Then, a smooth contact model is inserted into the parallel bond model, and the smooth contact model is used to simulate the contact between the bedding plane and the weak plane in the layered rock mass, so as to handle the normal stiffness and shear stiffness on the weak plane of the layered rock mass, enabling the generated rock mass model Mass to more accurately reflect the isotropic behavior of the material, thereby achieving precise control of the mechanical properties in different directions (parallel to the bedding plane direction and perpendicular to the bedding plane direction).

[0039] Step S103: In response to the bedding plane inclination angle of the rock mass model Mass being in the range of (0°, 90°), adjust the compressive strength of the rock mass model Mass by adjusting the strength parameters of the smooth contact interface.

[0040] When the bedding plane inclination angle of the rock mass model Mass is in the range of (0°, 90°), by adjusting the normal stiffness k n,sj and the tangential stiffness k s,sj , adjust the matrix modulus E0 parallel to the bedding plane direction and the rock mass modulus E 90 perpendicular to the bedding plane direction in the rock mass model Mass. Among them, when the error between the matrix modulus E0 parallel to the bedding plane direction in the rock mass model Mass and the matrix modulus E0 parallel to the bedding plane direction in the layered rock mass is greater than the preset first threshold, then adjust the tangential stiffness k s,sj of the smooth contact interface; when the error between the rock mass modulus E 90 perpendicular to the bedding plane direction in the rock mass model Mass and the rock mass modulus E 90 perpendicular to the bedding plane direction in the layered rock mass is greater than the preset second threshold, then adjust the normal stiffness k n,sj of the smooth contact interface.

[0041] Therefore, the interaction between the normal stiffness and the shear stiffness is fully considered in the calibration process of the model, enabling the model to not only perform well in the bedding plane directions of 0° (parallel to the bedding plane direction) and 90° (perpendicular to the bedding plane direction), but also exhibit more accurate mechanical behavior in other bedding plane inclination angles (such as 45°, 60°), effectively improving the prediction performance and practicality of the model, and being more reliable in engineering applications.

[0042] In addition, when the bedding plane of the rock mass model Mass is inclined within the range of (0°, 90°), by adjusting the tensile strength σ t,sj , cohesion c sj and friction angle μ sj of the smooth contact interface, the uniaxial compressive strength in the rock mass model Mass is determined. Specifically, when the error between the uniaxial compressive strength UCS0 parallel to the loading direction in the rock mass model Mass and the uniaxial compressive strength UCS0 parallel to the loading direction in the layered rock mass is greater than a preset third threshold, the friction angle μ sj of the smooth contact interface in the rock mass model Mass is kept unchanged, and the tensile strength σ t,si , cohesion c sj are adjusted.

[0043] In this application, by means of the bonding relationship between the bonded particles in the parallel bonding model, the elastic modulus and strength of the material matrix are simulated, and the matrix stiffness of the matrix model of the layered rock mass is adjusted; by inserting a smooth contact interface in the matrix model Matrix through the smooth contact model to simulate the contact between the bedding plane and the weak plane of the layered rock mass, and then, through the normal stiffness and shear stiffness on the weak plane at the smooth contact interface, the bedding plane stiffness of the rock mass model of the layered rock mass is adjusted. This makes the adjustment of parameters more targeted, reduces the coupling influence of the interaction between parameters in the complex model calibration process, improves the transparency of parameter adjustment, reduces unnecessary parameter iteration, and improves the stability of the model. By simultaneously adjusting the matrix parameters in the parallel bonding model and the weak plane parameters in the smooth contact model, the performance of the model under different directions and stress conditions becomes more consistent.

[0044] Combining the parallel bonding model and the smooth contact model, the calibration process of the layered rock mass model parameters is divided into multiple independent calibration steps, and in each independent calibration step, different parameters of the layered rock mass model are independently calibrated in sequence. The calibration process has multi-scale and multi-directional adaptability, so that the calibration of the layered rock mass model is no longer affected by the model scale size. At the same time, it effectively avoids the parameter coupling during the calibration of the layered rock mass model in the prior art, not only simplifies the parameter adjustment process of the model, improves the model calibration efficiency, but also makes the calibration process of the layered rock mass model repeatable, greatly enhances the adaptability of the model. Compared with the traditional trial-and-error method, it shortens the calibration time, improves the calibration efficiency, and makes the calibration of the model repeatable and standardized, reducing the influence of human experience.

[0045] The embodiment of this application also provides a parameter verification system for a transversely isotropic rock mass discrete element model, as Figure 4 shown. This system includes:

[0046] A matrix model and a matrix stiffness adjustment unit 401 are configured to generate a matrix model Matrix of a layered rock mass based on a parallel bond model, calibrate the matrix modulus of the matrix model Matrix by adjusting the bond stiffness of the bonded particles in the matrix model Matrix, and calibrate the compressive strength of the matrix model Matrix by adjusting the strength parameters of the bonded particles;

[0047] A rock mass model generation unit 402 is configured to insert a smooth contact interface for simulating the contact of bedding planes and weak planes of the layered rock mass into the matrix model Matrix based on a smooth contact model to generate a rock mass model Mass of the layered rock mass;

[0048] A bedding plane stiffness adjustment unit 403 is configured to adjust the compressive strength of the rock mass model Mass by adjusting the strength parameters of the smooth contact interface in response to the bedding plane dip angle of the rock mass model Mass being in the range of (0°, 90°).

[0049] The transverse isotropic rock mass discrete element model parameter calibration system provided by the embodiments of the present application can implement the steps and processes of the transverse isotropic rock mass discrete element model parameter calibration method in any of the above embodiments and achieve the same technical effects, which will not be elaborated here one by one.

[0050] In the description of the present invention, it should be understood that the terms "first" and "second" are only used for descriptive purposes and cannot be construed as indicating or implying relative importance or implicitly specifying the quantity of the indicated technical features. Thus, the features defined with "first" and "second" may explicitly or implicitly include at least one of such features. In the description of the present invention, the meaning of "a plurality" is at least two, such as two, three, etc., unless otherwise specifically defined.

[0051] In the present invention, the terms "an embodiment", "some embodiments", "example", "specific example", or "some examples", etc. mean that the specific features, structures, materials, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described can be combined in a suitable manner in any one or more embodiments or examples.

[0052] The above are only the preferred embodiments of the present application and are not used to limit the present application. For those skilled in the art, the present application can have various changes and modifications. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for calibrating parameters of a transversely isotropic rock mass discrete element model, characterized in that, Including: Step S101: Generate a complete matrix model Matrix of the layered rock mass based on the parallel bond model, calibrate the matrix modulus of the matrix model Matrix by adjusting the bond stiffness of the bonded particles in the matrix model Matrix, and calibrate the compressive strength of the matrix model Matrix by adjusting the strength parameters of the bonded particles; Step S102: Insert a smooth contact interface for simulating the contact between the bedding plane and the weak plane of the layered rock mass into the matrix model Matrix based on the smooth contact model to generate the rock mass model Mass of the layered rock mass; Step S103: In response to the bedding plane dip angle of the rock mass model Mass being (0°, 90°), adjust the compressive strength of the rock mass model Mass by adjusting the strength parameters of the smooth contact interface.

2. The method for calibrating parameters of the transversely isotropic rock mass discrete element model according to claim 1, characterized in that In step S101, the calibration of the matrix modulus of the matrix model Matrix by adjusting the bond stiffness between the bonded particles in the matrix model Matrix is specifically as follows: Adjust the normal stiffness k of the bonded particles in the Matrix model n,pbm and the shear stiffness k s,pbm to calibrate the effective modulus E of the Matrix model pbm .

3. The brittle and shear failure analysis method of the layered rock mass according to claim 1, characterized in that In step S101, the calibration of the compressive strength of the matrix model Matrix by adjusting the strength parameters of the bonded particles is specifically as follows: Adjust the tensile strength σ of the bonded particles in the Matrix model t,pbm , the cohesion c pbm and the friction angle μ pbm to calibrate the compressive strength of the Matrix model.

4. The method for checking the parameters of the transversely isotropic rock mass discrete element model according to claim 1, characterized in that, In step S103, the adjustment of the compressive strength of the rock mass model Mass by adjusting the strength parameters of the smooth contact interface in response to the bedding plane dip angle of the rock mass model Mass being (0°, 90°) is specifically as follows: By adjusting the normal stiffness k of the smooth contact interface n,j and the tangential stiffness k s,sj , the matrix modulus E0 in the direction parallel to the bedding plane and the rock mass modulus E in the direction perpendicular to the bedding plane in the rock mass model Ma are adjusted; 90 ​ Adjust the tensile strength σ of the smooth contact interface t,sj , the cohesion c sj and the friction angle μ sj , for the uniaxial compressive strength in the rock mass model Mass.

5. The method for calibrating the parameters of the transversely isotropic rock mass discrete element model according to claim 4, wherein By adjusting the normal stiffness k of the smooth contact interface n,sj and the tangential stiffness k s,sj , the matrix modulus E0 in the direction parallel to the bedding plane and the rock mass modulus E in the direction perpendicular to the bedding plane in the rock mass model Mass are adjusted, including: 90 ​ If the error between the matrix modulus E0 in the rock mass model Mass parallel to the bedding plane direction and the matrix modulus E0 in the layered rock mass parallel to the bedding plane direction is greater than a preset first threshold, then the tangential stiffness k of the smooth contact interface s,sj is adjusted; In response to the modulus of rock mass E in the direction perpendicular to the bedding plane in the rock mass model Mass 90 and the modulus of rock mass E in the direction perpendicular to the bedding plane in the layered rock mass 90 whose error is greater than a preset second threshold, the normal stiffness k of the smooth contact interface n,sj is adjusted.

6. The method for calibrating the parameters of the transversely isotropic rock mass discrete element model according to claim 4, characterized in that, Adjusting the tensile strength σ of the smooth contact interface t,sj , the cohesion c sj and the friction angle μ sj , for the uniaxial compressive strength in the rock mass model Mass, including: If the error between the uniaxial compressive strength UCS0 parallel to the loading direction in the rock mass model Mass and the uniaxial compressive strength UCS0 parallel to the loading direction in the layered rock mass is greater than a preset third threshold, then the friction angle μ of the smooth contact interface in the rock mass model Mass is maintained sj unchanged, and the tensile strength σ t,si , cohesion c sj is adjusted.

7. A parameter calibration system for a transversely isotropic rock mass discrete element model, characterized in that, Including: A matrix model and a matrix stiffness adjustment unit, configured to generate a matrix model Matrix of the layered rock mass based on the parallel bond model, calibrate the matrix modulus of the matrix model Matrix by adjusting the bond stiffness of the bonded particles in the matrix model Matrix, and calibrate the compressive strength of the matrix model Matrix by adjusting the strength parameters of the bonded particles; A rock mass model generation unit, configured to insert a smooth contact interface for simulating the contact between the bedding plane and the weak plane of the layered rock mass into the matrix model Matrix based on the smooth contact model to generate the rock mass model Mass of the layered rock mass; A bedding plane stiffness adjustment unit, configured to adjust the compressive strength of the rock mass model Mass by adjusting the strength parameters of the smooth contact interface in response to the bedding plane dip angle of the rock mass model Mass being (0°, 90°).

Images