A method and system for calculating the stability of a slope under the action of a coal-rock combination

Abstract

The application discloses a kind of coal rock combination under the action of slope stability calculation method and system, it is related to slope stability analysis technical field, method includes: analysis research area engineering geological characteristics, determine the form parameters of slope and coal rock mass in research area and the physical and mechanical parameters of each stratum rock mass;Based on the form parameters of most dangerous sliding surface searched by Geo-Studio geotechnical analysis software;Based on the form parameters of most dangerous sliding surface, stratum tendency distribution situation will main sliding surface vertical strip column be divided into micro strip column, analyze the three-dimensional mechanical effect of micro strip column bottom sliding surface;Equivalent three-dimensional mechanical effect of micro strip column bottom sliding surface to sliding body main sliding surface, the superimposed sliding force of each row micro strip column is obtained, superimposed anti-sliding force is divided into stages, determine equivalent physical and mechanical parameters;Equivalent physical and mechanical parameters are introduced into Morgenstern-Price method, set the bulk density λ and the slope stability coefficient F s Under the action of coal rock combination, f (x) expression satisfies boundary condition, then output F s The present application accurately predicts the slope stability under the action of coal rock combination.

Description

Technical Field

[0001] This invention relates to the field of slope stability analysis technology, and in particular to a method and system for calculating slope stability under the action of coal and rock composites. Background Technology

[0002] During open-pit mining, it is crucial to ensure the stability of mine slopes to guarantee safe production. Complex slopes formed under the influence of coal and rock masses exhibit significant three-dimensional effects, directly threatening the safety of miners and machinery. Therefore, stability analysis is essential. Existing three-dimensional stability analysis methods primarily rely on two-dimensional calculations. For example, the application No. 20211140557, "A Slope Stability Prediction Method," first determines geometric prediction coefficients based on the slope angle and height, then determines mechanical prediction coefficients based on unit weight, cohesion, friction angle, and void pressure ratio, and finally predicts the slope's stability using these parameters. However, this two-dimensional calculation method is generally only applicable to solving the stability of homogeneous slopes with regular shapes and cannot meet the needs of slope stability calculations under the influence of coal and rock masses in actual production. Therefore, a method for calculating slope stability coefficients under the influence of coal and rock masses is urgently needed to provide a theoretical basis for determining the stability coefficients of complex open-pit coal mine slopes. Summary of the Invention

[0003] The purpose of this invention is to provide a method and system for calculating slope stability under the action of coal and rock composites. By establishing a three-dimensional model of the slope under the action of coal and rock composites, the three-dimensional mechanical effects of the sliding surface at the bottom of the micro-strip columns and the superimposed mechanical effects of each row of micro-strip columns are analyzed in different regions, and the accurate slope stability coefficient under the action of coal and rock composites is obtained.

[0004] To achieve the above objectives, the present invention provides the following solution:

[0005] A method for calculating slope stability under the action of coal and rock composites includes the following steps:

[0006] S1: Analyze the regional engineering geological characteristics and determine the morphological parameters of slopes and coal and rock masses in the region, as well as the physical and mechanical parameters of each stratum and rock mass.

[0007] S2: Based on Geo-Studio geotechnical analysis software, search for potential sliding surfaces, analyze the morphology of potential sliding surfaces, and determine the morphological parameters of the most dangerous sliding surface;

[0008] S3: Based on the morphological parameters and stratigraphic dip distribution of the most dangerous slip surface, the main slip surface is vertically divided into several micro-slip surfaces, and the three-dimensional mechanical effects of the slip surface at the bottom of the micro-slip surfaces are analyzed.

[0009] S4: Superimpose the three-dimensional mechanical effects of the bottom sliding surface of each micro-strip column onto the main sliding surface of the sliding body, solve the superimposed sliding force X and superimposed anti-slip force T of each row of micro-strip columns in stages, and determine the equivalent physical and mechanical parameters;

[0010] S5: Introduce equivalent physical and mechanical parameters into the Morgenstern-Price method, setting the unit weight λ and the slope stability coefficient F under the action of the coal-rock composite. s The initial value is set according to the elasticity theory, and the expression of the function f(x) is set accordingly;

[0011] S6: Define the f(x) function and the assumptions λ, F s Substitute the initial value into the boundary conditions and determine whether the boundary conditions are met. If not, return to S5; if met, output F. s .

[0012] Furthermore, the analysis of the engineering geological characteristics of the area, the determination of the morphological parameters of the slopes and coal and rock masses within the area, and the determination of the physical and mechanical parameters of the strata and rock masses specifically include:

[0013] S1.1: Approximating the potential slip surface of the slope under the action of the coal-rock composite as an ellipsoid, a three-dimensional model under the action of the coal-rock composite is established, and the slip surface equation z is obtained. h as follows:

[0014]

[0015] In the formula, (x0, z0) are the coordinates of the center of the circular arc segment on the main sliding surface of the slope; R a R b These are the major axis radius and minor axis radius of the ellipsoid, respectively.

[0016] The slope equation z is obtained by dividing the region according to the shape of the sliding body boundary. p as follows:

[0017]

[0018] In the formula, H and h are Z p1 With Z p3 The distance from the slope to the x-axis; θ and η are the distances from the slope to the x-axis, respectively. p2 With Z p4 The angle between the slope and the x-axis; d is the Z-axis. p3 The length of the slope along the x-axis;

[0019] The slip surface of the slope is a combination of an ellipsoid and a weak layer. The dip angles of the weak layer in the x-axis and y-axis directions are β and δ, respectively. In this case, the slip surface equation of the slope is as follows:

[0020]

[0021] S1.2: Based on the morphological parameters of the weak layer, determine the spatial morphological equation of the weak layer, and combine the morphological parameters of the internal spoil heap and the lithological stratigraphic model to construct a three-dimensional stability calculation model of the slope under the action of the coal-rock composite.

[0022] Furthermore, the physical parameters include cohesion c and internal friction angle. The bulk density γ; the morphological parameters include the slope height H / m of the inner spoil heap, the slope angle θ / (°) of the inner spoil heap, the coal pillar height h / m above the weak surface, the bottom angle η / (°) of the coal pillar, the dip angle δ / (°) of the weak surface along the dip direction, the coordinates of the center, and the minor axis length R of the sliding body. a / m, length of the major axis of the sliding body R b / m.

[0023] Furthermore, S3: Based on the morphological parameters and stratigraphic dip distribution of the most dangerous slip surface, the main slip surface is vertically divided into several micro-slip columns, and the three-dimensional mechanical effects of the slip surface at the bottom of the micro-slip columns are analyzed, specifically including:

[0024] S3.1: Treat any microstrip column as a hexahedron with mutually perpendicular sides, and perform mechanical analysis on any microstrip column based on the three-dimensional limit equilibrium slice method;

[0025] S3.2: Ignore all tangential and normal forces between the micro-strips, and make the tangential force on the bottom sliding surface parallel to the main sliding direction. Assume that the point of application of the normal force at the bottom of the micro-strip is located at the midpoint of the bottom of the strip.

[0026] S3.3: Based on the force balance condition, the normal force and sliding force of the bottom sliding surface of the micro-strip column are obtained, and the anti-slip force of the micro-strip column is obtained based on the Mohr-Coulomb strength criterion.

[0027] S3.4: The volume force of the micro-strip column is obtained by analyzing the spatial morphology of the slope and the sliding surface.

[0028] Furthermore, in step S4: the three-dimensional mechanical effect of the bottom sliding surface of each micro-strip column is equivalently applied to the main sliding surface of the sliding body, and the superimposed sliding force X and superimposed anti-slip force T of each row of micro-strip columns are solved in stages to determine the equivalent physical and mechanical parameters, specifically including:

[0029] By combining the equivalent mechanical effects on the main sliding surface with the superimposed mechanical effects of each row of micro-strips, the equivalent unit weight γ of the sliding body at each stage can be calculated. d Equivalent cohesion c d Equivalent friction coefficient

[0030] Furthermore, in S5: equivalent physical and mechanical parameters are introduced into the Morgenstern-Price method, setting the unit weight λ and the slope stability coefficient F under the action of the coal-rock composite. sThe initial value of f(x) is set according to elasticity theory, and the expression of the function f(x) is specifically included as follows:

[0031] Assume the normal force E' of the column i Tangential force X of the bar column i The relationship between them is:

[0032] X i =λf(x)E' i

[0033] The boundary conditions are:

[0034]

[0035] A system for calculating the slope stability coefficient under the action of a coal-rock composite, applied to the method for calculating the slope stability coefficient under the action of the coal-rock composite, the system comprising:

[0036] Slope data acquisition module: Analyzes regional engineering geological characteristics, determines the morphological parameters of slopes and coal and rock masses within the region, and the physical and mechanical parameters of various strata and rock masses;

[0037] Most dangerous slip surface analysis module: Input the morphological parameters of the slope and coal and rock mass in the area and the physical and mechanical parameters of the strata and rock mass in each area into Geo-Studio geotechnical analysis software, use Geo-Studio geotechnical analysis software to search for potential slip surfaces, analyze the morphology of potential slip surfaces, and determine the morphological parameters of the most dangerous slip surface;

[0038] Strip column division module: Based on the morphological parameters of the most dangerous slip surface and the distribution of stratigraphic dip, the main slip surface is vertically divided into strip columns to obtain several micro-strip columns, and the three-dimensional mechanical effects of the slip surface at the bottom of the micro-strip columns are analyzed.

[0039] Equivalent physical and mechanical parameters solution module: The three-dimensional mechanical effects of the bottom sliding surface of each micro-strip column are superimposed onto the main sliding surface of the sliding body, and the superimposed sliding force X and superimposed anti-slip force T of each row of micro-strip columns are solved in stages to determine the equivalent physical and mechanical parameters;

[0040] Initial value setting module: Introduces equivalent physical and mechanical parameters into the Morgenstern-Price method, setting the unit weight λ and the slope stability coefficient F under the action of coal and rock composite. s The initial value is set according to the elasticity theory, and the expression of the function f(x) is set accordingly;

[0041] Stability coefficient output module: outputs the set f(x) function and the assumptions λ and F. s Substitute the initial values ​​into the boundary conditions, continuously adjust the values ​​of the function f(x) and λ, and determine whether the boundary conditions are satisfied. Continue until the boundary conditions are satisfied, then output F. s .

[0042] According to specific embodiments provided by the present invention, the present invention discloses the following technical effects: The slope stability calculation method under the action of coal and rock composite provided by the present invention, based on establishing the spatial morphology equation of the sliding body, determining the mechanical model of the strip pillars, the strip pillar division method and mechanical assumptions, and combining the spatial morphology of the sliding body itself, divides it into 5 calculation stages, and uses the Mohr-Coulomb shear strength theory to characterize the three-dimensional mechanical effect of the strip pillars in stages; starting from the essence of the three-dimensional mechanical effect of the strip pillars, the equivalent physical and mechanical parameters of each strip on the main sliding surface of the sliding body are analyzed in stages using the cumulative equivalent method, and introduced into the Morgenstern-Price method to construct a three-dimensional stability calculation method for the slope of the internal spoil heap under coal pillar support. It can obtain accurate stability coefficients and is simple to operate, and can provide a theoretical basis for determining the stability coefficient of complex slopes in open-pit coal mines. Attached Figure Description

[0043] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0044] Figure 1 This is a flowchart of the slope stability calculation method under the action of coal and rock composites according to the present invention;

[0045] Figure 2 This is a schematic diagram of the weak-level spatial morphology in an embodiment of the present invention;

[0046] Figure 3 This is a schematic diagram of the cross-sectional morphology of the weak layer along the strike in an embodiment of the present invention, wherein A represents the calculated length of the coal-rock composite along the strike of approximately 293m;

[0047] Figure 4 This is a schematic diagram of the cross-sectional morphology of the weak layer along the dip direction in an embodiment of the present invention;

[0048] Figure 5 This is a schematic diagram of a three-dimensional slope model under the action of coal and rock composite in an embodiment of the present invention;

[0049] Figure 6 This is a three-dimensional stability calculation model of the slope under the action of coal and rock combination in the embodiments of the present invention;

[0050] Figure 7 This is a schematic diagram of the main smooth noodle block division in an embodiment of the present invention;

[0051] Figure 8 This is a schematic diagram of the force analysis of the micro-strip column model in an embodiment of the present invention;

[0052] Figure 9 This is a schematic diagram of the force analysis model and elements of the micro-strip column bottom interface in an embodiment of the present invention;

[0053] Figure 10 This is a schematic diagram (top view) of the computational region division in an embodiment of the present invention;

[0054] Figure 11 This is a schematic diagram illustrating the number of columns divided into regions and the conversion of superposition angles in an embodiment of the present invention;

[0055] Figure 12 This is a schematic diagram of the superposition of mechanical effects of the main sliding surface and the force analysis of the strip in an embodiment of the present invention (taking stage V as an example);

[0056] Figure 13 The results of the three-dimensional stability calculation of the slope under the action of coal and rock combination in the embodiments of the present invention are shown. Among them, a, b, and c refer to the three-dimensional stability calculation results of the slope of the inner spoil heap under the conditions that the top width of the coal pillar is 20m, 30m, and 38m, respectively. Detailed Implementation

[0057] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0058] The purpose of this invention is to provide a method and system for calculating slope stability under the action of coal and rock composites. By establishing a three-dimensional model of the slope under the action of coal and rock composites, the three-dimensional mechanical effects of the sliding surface at the bottom of the micro-strip columns and the superimposed mechanical effects of each row of micro-strip columns are analyzed in different regions, and the accurate slope stability coefficient under the action of coal and rock composites is obtained.

[0059] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0060] This embodiment takes an open-pit coal mine located in Inner Mongolia as an example. The weak layer between the two coal seams has a dip angle of 4°, and the coal-rock composite has a calculated length of approximately 293m along the strike. The spatial morphology of the weak layer is as follows: Figure 2 As shown, the cross-section along the strike is as follows Figure 3 As shown, the weak layer is along the dip section as follows Figure 4 As shown in Table 1, the physical and mechanical parameters of each selected soil and rock layer are as follows.

[0061] Table 1 Physical and mechanical parameters

[0062]

[0063] like Figure 1 As shown in the figure, the slope stability calculation method under the action of coal and rock composite provided in this embodiment of the invention includes the following steps:

[0064] S1: Analyze the regional engineering geological characteristics and determine the morphological parameters of slopes and coal and rock masses in the region, as well as the physical and mechanical parameters of each stratum and rock mass.

[0065] S2: Based on Geo-Studio geotechnical analysis software, potential sliding surfaces are searched, their morphologies are analyzed, and the morphological parameters of the most dangerous sliding surface are determined; the landslide model is set to a combined sliding with a circular arc as the side interface and weak interlayers between coal seams as the base, such as... Figure 7 As shown, the morphological parameters are shown in Table 2;

[0066] Table 2 Slip surface morphology parameters

[0067]

[0068] S3: Based on the morphological parameters and stratigraphic dip distribution of the most dangerous slip surface, the main slip surface is vertically divided into several micro-slip surfaces, and the three-dimensional mechanical effects of the slip surface at the bottom of the micro-slip surfaces are analyzed.

[0069] S4: Superimpose the three-dimensional mechanical effects of the bottom sliding surface of each micro-strip column onto the main sliding surface of the sliding body, solve the superimposed sliding force X and superimposed anti-slip force T of each row of micro-strip columns in stages, and determine the equivalent physical and mechanical parameters;

[0070] S5: Introduce equivalent physical and mechanical parameters into the Morgenstern-Price method, setting the unit weight λ and the slope stability coefficient F under the action of the coal-rock composite. s The initial value is set according to the elasticity theory, and the expression of the function f(x) is set accordingly;

[0071] S6: Define the f(x) function and the assumptions λ, F s Substitute the initial value into the boundary conditions and determine whether the boundary conditions are met. If not, return to S5; if met, output F. s .

[0072] In this embodiment, S1: Analyzing the regional engineering geological characteristics, determining the morphological parameters of slopes and coal and rock masses within the region, and determining the physical and mechanical parameters of each stratum and rock mass, specifically includes:

[0073] S1.1: Approximate the potential sliding surface of the slope under the action of the coal-rock composite as an ellipsoid, and establish a three-dimensional model of the slope under the action of the coal-rock composite (e.g., Figure 5 As shown), the equation of the slip surface z is obtained. h as follows:

[0074]

[0075] In the formula, (x0, z0) are the coordinates of the center of the circular arc segment on the main sliding surface of the slope; R a R b These are the major axis radius and minor axis radius of the ellipsoid, respectively, in meters (m).

[0076] The region is divided according to the boundary morphology of the sliding body, and the slope equation z p The regions can be divided as follows:

[0077]

[0078] In the formula, H and h are Z p1 With Z p3 The distance from the slope to the x-axis; θ and η are the distances from the slope to the x-axis, respectively. p2 With Z p4 The angle between the slope and the x-axis; d is the Z-axis. p3 The length of the slope along the x-axis;

[0079] The slip surface of the slope is a combination of an ellipsoid and a weak layer. The dip angles of the weak layer in the x-axis and y-axis directions are β and δ, respectively. In this case, the slip surface equation of the slope is as follows:

[0080]

[0081] S1.2: In 3D modeling software, based on the morphological parameters of the weak layer, the spatial morphological equation of the weak layer is determined. Combined with the morphological parameters of the internal spoil heap and the lithological stratigraphic model, a 3D slope stability calculation model under the action of coal-rock composite is constructed, such as... Figure 6 As shown.

[0082] In this embodiment, the physical parameters include cohesion c and internal friction angle. The bulk density γ; the morphological parameters include the slope height H / m of the inner spoil heap, the slope angle θ / (°) of the inner spoil heap, the coal pillar height h / m above the weak surface, the bottom angle η / (°) of the coal pillar, the dip angle δ / (°) of the weak surface along the dip direction, the coordinates of the center, and the minor axis length R of the sliding body. a / m, length of the major axis of the sliding body R b / m.

[0083] In this embodiment, S3: Based on the morphological parameters and stratigraphic dip distribution of the most dangerous slip surface, the main slip surface is vertically divided into several micro-slip surfaces, and the three-dimensional mechanical effects of the slip surface at the bottom of the micro-slip surfaces are analyzed, specifically including:

[0084] S3.1: Treat any microstrip column as a hexahedron with mutually perpendicular sides, and perform mechanical analysis on any microstrip column based on the three-dimensional limit equilibrium slice method;

[0085] Based on the three-dimensional limit equilibrium slice method, mechanical analysis of any micro-slab reveals that each surface of the micro-slab, except the top surface, is subjected to three forces. The lateral surfaces are subjected to a normal force E and a tangential force T. The normal force E is indicated by a subscript representing the surface on which it acts, such as E... y The normal force T represents the force acting on a surface with the y-normal vector. The first subscript of the tangential force T indicates the surface on which it acts, and the second subscript indicates the direction of the force. For example, Tyx represents a tangential force acting in the x-direction on a surface with the y-normal vector. The normal force d acting on the bottom surface... N Anti-slip force d T , volume force dW; α x α y Let be the angles between the lines of intersection of the bottom sliding surface with planes xoz and yoz, and the horizontal plane; let d be the length of the x-direction. x Let the length in the y-direction be d. y The stress analysis model of the column is as follows: Figure 8 As shown;

[0086] S3.2: ① Ignore all tangential and normal forces between the micro-strip columns, and ensure that the tangential force on the bottom sliding surface is parallel to the principal sliding direction, i.e., ignore the tangential force in the zoy plane. The anti-slip force in the principal sliding direction is expressed as dT. x Indicates: ① Assuming the point of application of the normal force at the bottom of the micro-strip column is located at the midpoint of the bottom of the column; the force analysis model and elements of the micro-strip column are shown in the figure below. Figure 9 As shown.

[0087] S3.3: Based on the force balance condition, the normal force dN and the sliding force dX of the bottom sliding surface of the micro-strip column are obtained, and the anti-slip force of the micro-strip column is obtained based on the Mohr-Coulomb strength criterion.

[0088] The normal force dN on the sliding surface at the bottom of the micro-strip column is as follows: dN=dWcosα x ,

[0089] The expression for the sliding force dX on the bottom surface of the micro-strip column is as follows: dX=dWsinα x ,

[0090] Anti-slip force dT of micro-strip column x :

[0091] In the formula, Let be the internal friction angle of the microstrip column (°); c is the cohesion of the microstrip column, in kPa; the analytical formulas are: normal force dN, sliding force dX, and anti-slip force dT. x From these three equations, we can see that sinα x cosα x dA and dW are both unknown parameters. If we want to solve for the anti-slip force dT on the bottom sliding surface of the micro-strip column, xTo determine the sliding force dX, it is necessary to analyze the geometric relationships between various parameters and then indirectly solve for one or more known parameters. A schematic diagram of the geometric elements of the micro-strip column is shown below. Figure 9 As shown.

[0092] S3.4: First, using the trigonometric relationships and the definition of partial derivatives, we get:

[0093]

[0094] Where x0 is the X-axis coordinate of the center of the arc segment on the main sliding surface of the slope.

[0095] Then, based on the geometric relationships of the parameters of the sliding surface, the unknown parameter sinα can be determined. x cosα x dA is as follows:

[0096]

[0097]

[0098]

[0099] In the formula, α x α y These are the angles between the lines of intersection of the bottom sliding surface with plane xoz and plane yoz, and the horizontal plane, respectively.

[0100] S3.5: The volumetric force dW of the micro-strip column can be obtained by analyzing the slope and slip surface morphology.

[0101] dW=γ(z p -z h )dxdy

[0102] In the formula, γ is the unit weight of the microstrip column, in kN / m³. 3 It should be noted that due to the irregular shape of the spatial volume formed by the slope and slip surface, calculations need to be performed in seven regions, named A to G. The division of the calculation regions is illustrated below. Figure 10 As shown;

[0103] Therefore, the volume force dW of the micro bar column in the i-th row and j-th column is... ij It can be represented as follows:

[0104]

[0105] In this embodiment, S4: The three-dimensional mechanical effect of the bottom sliding surface of each micro-strip column is equivalent to the main sliding surface of the sliding body, and the superimposed sliding force X and superimposed anti-slip force T of each row of micro-strip columns are solved in stages to determine the equivalent physical and mechanical parameters, including:

[0106] By combining the equivalent mechanical effects on the main sliding surface with the superimposed mechanical effects of each row of micro-strips, the equivalent unit weight γ of the sliding body at each stage can be calculated. d Equivalent cohesion c d Equivalent friction coefficient Force analysis of the strip column on the main sliding surface at each stage is as follows: Figure 12 As shown in the figure. The calculated physical and mechanical parameters, position and morphological parameters of each block are shown in Table 3.

[0107] Table 3 Calculation parameters for each block on the main sliding surface

[0108]

[0109] In a further embodiment, step S4: the three-dimensional mechanical effect of the bottom sliding surface of each micro-strip column is equivalently applied to the main sliding surface of the sliding body, and the superimposed sliding force X and superimposed anti-slip force T of each row of micro-strip columns are solved in stages to determine the equivalent physical and mechanical parameters, specifically including:

[0110] S4.1: The mechanical effects are superimposed into five stages, namely stages I to V. Let stage I have n bars along the x-direction, and region A in row r have r1 bars. The diagram showing the number of bars in each region and the conversion of superposition angles is as follows. Figure 11 As shown. The superimposed downward force X of the r-th row of bars in stage I. r With anti-slip force T r as follows:

[0111]

[0112]

[0113] In the formula,

[0114]

[0115] S4.2: Let k be the number of columns along the x-direction in stage II, and u1 be the number of columns in region B of row u. Then the superimposed sliding force X of the columns in row u is... u With anti-slip force T u as follows:

[0116]

[0117]

[0118] In the formula, cosω uj =cos(α) xuj -θ)

[0119] S4.3: Let the number of columns along the x-direction in stage III be p, and the number of columns in regions B and C of row s be s1+s3 and s2 respectively. The superimposed sliding force X of the columns in row s of stage III is... s With anti-slip force T s as follows:

[0120]

[0121]

[0122] In the formula, cosω sj =cos(α) xsj -θ)

[0123] S4.4: Let the number of columns along the x-direction in stage IV be q, and the number of columns in regions D and E of the m-th row be m1+m3 and m2, respectively; the superimposed sliding force X of the columns in the m-th row of stage IV. m With anti-slip force T m :

[0124]

[0125]

[0126] In the formula, cosω mj =cos(α) xmj -θ)

[0127] S4.5: Let the number of columns along the x-direction in stage V be i, and the number of columns in regions F and G of row t be t1+t3 and t2, respectively. The superimposed sliding force X of the columns in row t of stage V is... t With anti-slip force T t :

[0128]

[0129]

[0130] In the formula, cosω tj =cos(α) xtj -θ)

[0131] S4.6: Transform the three-dimensional mechanics problem into a plane strain problem. When dividing the columnar sections, set the base of the columnar section to a unit length in the y-direction. By combining like terms of the equivalent mechanical effects on the main sliding surface and the superimposed mechanical effects of each row of micro-slabs, calculate the equivalent unit weight γ of the sliding body at each stage. d Equivalent cohesion c d Equivalent friction coefficient Force analysis of the strip column on the main sliding surface at each stage is as follows: Figure 11As shown. Let X be the sliding force and anti-slip force of the bottom interface after the two-dimensional equivalent of the three-dimensional mechanical effect in stage I. Ⅰd T Ⅰd Then, for the r-th bar on the main sliding surface, we should have:

[0132]

[0133]

[0134] In the formula, y r Let be the length of the sliding boundary corresponding to the bar in row r.

[0135] According to X Ⅰdr =X r T Ⅰdr =T r have to:

[0136]

[0137]

[0138]

[0139] S4.7: Let X be the sliding force and anti-slip force of the bottom interface after the two-dimensional equivalent of the three-dimensional mechanical effect in stage II. Ⅱd T Ⅱd Then, for the u-th bar on the main sliding surface, we should have:

[0140]

[0141]

[0142] According to X Ⅱdu =X u T Ⅱdu =T u have to:

[0143]

[0144]

[0145]

[0146] S4.8: Let X be the sliding force and anti-slip force of the bottom interface after the two-dimensional equivalent of the three-dimensional mechanical effect in stage III. Ⅲd T Ⅲd Then, for the s-th bar on the main sliding surface, we should have:

[0147] X Ⅲds =dW s sinδ=γ Ⅲds (zp2 -z h22 )sinδdx s y s

[0148]

[0149] In the formula, z h22 = -xtanδ.

[0150] According to X Ⅲds =X s T Ⅲds =T s have to:

[0151]

[0152]

[0153]

[0154] S4.9: Let X be the sliding force and anti-slip force of the bottom interface after the two-dimensional equivalent of the three-dimensional mechanical effect in stage IV. Ⅳd T Ⅳd Then, for the m-th bar on the main sliding surface, we should have:

[0155] X Ⅳdm =dW m sinδ=γ Ⅳdm (z p3 -z h22 )sinδdx m y m

[0156]

[0157] According to X Ⅳdm =X m T Ⅳdm =T m have to:

[0158]

[0159]

[0160]

[0161] S4.10: Let X be the sliding force and anti-slip force of the bottom interface after the two-dimensional equivalent of the three-dimensional mechanical effect in stage V. Ⅴd T Ⅴd Then, for the t-th bar on the main sliding surface, we should have:

[0162] X Ⅴdt =dW tsinδ=γ Ⅴdt (z p4 -z h22 )sinδdx t y t

[0163]

[0164] According to X Ⅴdt =X t T Ⅴdt =T t have to:

[0165]

[0166]

[0167]

[0168] The equivalent physical and mechanical parameters of each block on the main sliding surface were calculated, as shown in Table 4.

[0169] Table 4 Equivalent physical and mechanical parameters of each block on the main sliding surface

[0170]

[0171]

[0172] In this embodiment, step S5 involves introducing equivalent physical and mechanical parameters into the Morgenstern-Price method, setting the unit weight λ and the slope stability coefficient F under the action of the coal-rock composite. s The initial value of f(x) is set according to elasticity theory, and the expression of the function f(x) is specifically included as follows:

[0173] Assume the normal force E' of the column i Tangential force X of the bar column i The relationship between them is:

[0174] X i =λf(x)E' i

[0175] The boundary conditions are expressed as follows:

[0176]

[0177] Where, x n Let M(x) be the coordinate of the center point of the nth column on the x-axis; n ) represents the bending moment of the column.

[0178] In this embodiment, step S6: setting the f(x) function and assuming λ, F sSubstitute the initial value into the boundary conditions and determine whether the boundary conditions are met. If not, return to S5; if met, output F. s Specifically, it includes:

[0179] The calculation process terminates when no external forces act on the slope sliding body (boundary condition). Based on elasticity theory, an expression for the function f(x) is assumed, and λ and F are continuously adjusted. s The value is substituted into the boundary conditions for calculation. It is then determined whether the boundary conditions are met. If not, S5 is returned; if they are met, the corresponding F value is output. s This is the slope stability coefficient under the action of the coal and rock composite, such as... Figure 13 As shown, this embodiment calculates the slope stability of the coal-rock composite body under different designed top widths.

[0180] Example 2

[0181] Based on the slope stability calculation method under coal-rock composite action provided in Embodiment 1 of the present invention, this embodiment provides a slope stability calculation system under coal-rock composite action, the system comprising:

[0182] Slope data acquisition module: Analyzes regional engineering geological characteristics, determines the morphological parameters of slopes and coal and rock masses within the region, and the physical and mechanical parameters of various strata and rock masses;

[0183] Most dangerous slip surface analysis module: Input the morphological parameters of the slope and coal and rock mass in the area and the physical and mechanical parameters of the strata and rock mass in each area into Geo-Studio geotechnical analysis software, use Geo-Studio geotechnical analysis software to search for potential slip surfaces, analyze the morphology of potential slip surfaces, and determine the morphological parameters of the most dangerous slip surface;

[0184] Strip column division module: Based on the morphological parameters of the most dangerous slip surface and the distribution of stratigraphic dip, the main slip surface is vertically divided into strip columns to obtain several micro-strip columns, and the three-dimensional mechanical effects of the slip surface at the bottom of the micro-strip columns are analyzed.

[0185] Equivalent physical and mechanical parameters solution module: The three-dimensional mechanical effects of the bottom sliding surface of each micro-strip column are superimposed onto the main sliding surface of the sliding body, and the superimposed sliding force X and superimposed anti-slip force T of each row of micro-strip columns are solved in stages to determine the equivalent physical and mechanical parameters;

[0186] Initial value setting module: Introduces equivalent physical and mechanical parameters into the Morgenstern-Price method, setting the unit weight λ and the slope stability coefficient F under the action of coal and rock composite. s The initial value is set according to the elasticity theory, and the expression of the function f(x) is set accordingly;

[0187] Stability coefficient output module: outputs the set f(x) function and the assumptions λ and F. sSubstitute the initial values ​​into the boundary conditions, continuously adjust the values ​​of the function f(x) and λ, and determine whether the boundary conditions are satisfied. Continue until the boundary conditions are satisfied, then output F. s .

[0188] In summary, the slope stability calculation method and system provided by this invention, starting from the essence of the three-dimensional mechanical effect of the strip column, accumulates and equates the three-dimensional mechanical effect of the sliding body to the main sliding surface of the sliding body. Combining the spatial morphological characteristics of the sliding body itself, the equivalent physical and mechanical parameters of each strip on the main sliding surface are analyzed in stages. Combined with the Morgenstern-Price method, the accurate slope stability system coefficient under the action of the coal-rock combination can be calculated. Moreover, this method is simple and easy to operate.

[0189] The remaining technical features in this embodiment can be flexibly selected by those skilled in the art to meet different specific practical needs. However, it is obvious to those skilled in the art that these specific details are not necessary to implement the present invention. In other instances, to avoid obscuring the present invention, well-known components, structures, or parts are not specifically described, and all are within the scope of technical protection defined by the claims of the present invention.

[0190] Modifications and variations made by those skilled in the art that do not depart from the spirit and scope of this invention should be within the protection scope of the appended claims. In the above description, numerous specific details have been set forth to provide a thorough understanding of the invention. However, it will be apparent to those skilled in the art that these specific details are not necessary to practice the invention. In other instances, to avoid obscuring the invention, well-known techniques, such as specific construction details, operating conditions, and other technical conditions, have not been specifically described.

[0191] This document uses specific examples to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. Furthermore, those skilled in the art will recognize that, based on the ideas of the present invention, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for calculating the slope stability coefficient under the action of coal and rock composites, characterized in that, Includes the following steps: S1: Analyze the engineering geological characteristics of the study area and determine the morphological parameters of the slopes and coal and rock masses and the physical and mechanical parameters of the strata and rock masses in the study area; S2: Input the morphological parameters of the slope and coal and rock mass in the area and the physical and mechanical parameters of the strata and rock mass in each area into Geo-Studio geotechnical analysis software, use Geo-Studio geotechnical analysis software to search for potential sliding surfaces, analyze the morphology of potential sliding surfaces, and determine the morphological parameters of the most dangerous sliding surface; S3: Based on the morphological parameters and stratigraphic dip distribution of the most dangerous slip surface, the main slip surface is vertically divided into several micro-slip surfaces, and the three-dimensional mechanical effects of the slip surface at the bottom of the micro-slip surfaces are analyzed. S4: Superimpose the three-dimensional mechanical effects of the bottom sliding surface of each micro-strip column onto the main sliding surface of the sliding body, solve the superimposed sliding force X and superimposed anti-slip force T of each row of micro-strip columns in stages, and determine the equivalent physical and mechanical parameters; S5: Introduce equivalent physical and mechanical parameters into the Morgenstern-Price method, setting the unit weight λ and the slope stability coefficient F under the action of the coal-rock composite. s The initial value is set according to the elasticity theory, and the expression of the function f(x) is set accordingly; S6: Define the f(x) function and the assumptions λ, F s Substitute the initial value into the boundary conditions and determine whether the boundary conditions are met. If not, return to S5; if met, output F. s .

2. The method for calculating the slope stability coefficient under the action of coal and rock composite as described in claim 1, characterized in that, S1: Analyze the regional engineering geological characteristics, determine the morphological parameters of slopes and coal and rock masses within the region, and determine the physical and mechanical parameters of each stratum and rock mass, specifically including: S1.1: Approximating the potential slip surface of the slope under the action of the coal-rock composite as an ellipsoid, a three-dimensional model under the action of the coal-rock composite is established, and the slip surface equation z is obtained. h as follows: In the formula, (x0, z0) are the coordinates of the center of the circular arc segment on the main sliding surface of the slope, and R a R b These are the major axis radius and minor axis radius of the ellipsoid, respectively. The slope equation z is obtained by dividing the region according to the shape of the sliding body boundary. p as follows: In the formula, H and h are Z p1 With Z p3 The distance from the slope to the x-axis; θ and η are the distances from the slope to the x-axis, respectively. p2 With Z p4 The angle between the slope and the x-axis; d is the Z-axis. p3 The length of the slope along the x-axis; The slip surface of the slope is a combination of an ellipsoid and a weak layer. The dip angles of the weak layer in the x-axis and y-axis directions are β and δ, respectively. In this case, the slip surface equation of the slope is as follows: S1.2: Based on the morphological parameters of the weak layer, determine the spatial morphological equation of the weak layer, and combine the morphological parameters of the internal spoil heap and the lithological stratigraphic model to construct a three-dimensional stability calculation model of the slope under the action of the coal-rock composite.

3. The method for calculating the slope stability coefficient under the action of coal and rock composite as described in claim 1, characterized in that, The physical and mechanical parameters include cohesion c and internal friction angle. The bulk density γ; the morphological parameters include the slope height H / m of the inner spoil heap, the slope angle θ / (°) of the inner spoil heap, the coal pillar height h / m above the weak surface, the bottom angle η / (°) of the coal pillar, the dip angle δ / (°) of the weak surface along the dip direction, the coordinates of the center, and the minor axis length R of the sliding body. a / m, length of the major axis of the sliding body R b / m.

4. The method for calculating the slope stability coefficient under the action of coal and rock composite as described in claim 1, characterized in that, S3: Based on the morphological parameters and stratigraphic dip distribution of the most dangerous slip surface, the main slip surface is vertically divided into several micro-slip surfaces, and the three-dimensional mechanical effects of the slip surface at the bottom of the micro-slip surfaces are analyzed, specifically including: S3.1: Treat any microstrip column as a hexahedron with mutually perpendicular sides, and perform mechanical analysis on any microstrip column based on the three-dimensional limit equilibrium slice method; S3.2: Ignore all tangential and normal forces between the micro-strips, and make the tangential force on the bottom sliding surface parallel to the main sliding direction. Assume that the point of application of the normal force at the bottom of the micro-strip is located at the midpoint of the bottom of the strip. S3.3: Based on the force balance condition, the normal force and sliding force of the bottom sliding surface of the micro-strip column are obtained, and the anti-slip force of the micro-strip column is obtained based on the Mohr-Coulomb strength criterion. S3.4: The volume force of the micro-strip column is obtained by analyzing the spatial morphology of the slope and the sliding surface.

5. The method for calculating the slope stability coefficient under the action of coal and rock composite as described in claim 1, characterized in that, S4: The three-dimensional mechanical effect of the bottom sliding surface of each micro-strip column is equivalent to the main sliding surface of the sliding body. The superimposed sliding force X and superimposed anti-slip force T of each row of micro-strip columns are solved in stages to determine the equivalent physical and mechanical parameters, specifically including: By combining the equivalent mechanical effects on the main sliding surface with the superimposed mechanical effects of each row of micro-strips, the equivalent unit weight γ of the sliding body at each stage can be calculated. d Equivalent cohesion c d Equivalent friction coefficient 6. The method for calculating the slope stability coefficient under the action of coal and rock composite as described in claim 1, characterized in that, S5: Introducing equivalent physical and mechanical parameters into the Morgenstern-Price method, setting the unit weight λ and the slope stability coefficient F under the action of the coal-rock composite. s The initial value of f(x) is set according to elasticity theory, and the expression of the function f(x) is specifically included as follows: Assume the normal force E' of the column i Tangential force X of the bar column i The relationship between them is: X i =λf(x)E' i The boundary conditions are: Where, x n Let M(x) be the coordinate of the center point of the nth column on the x-axis; n ) represents the bending moment of the column.

7. A system for calculating the slope stability coefficient under the action of a coal-rock composite, applied to the method for calculating the slope stability coefficient under the action of a coal-rock composite as described in any one of claims 1-6, characterized in that, include: Slope data acquisition module: Analyzes regional engineering geological characteristics, determines the morphological parameters of slopes and coal and rock masses within the region, and the physical and mechanical parameters of various strata and rock masses; Most dangerous slip surface analysis module: Input the morphological parameters of the slope and coal and rock mass in the area and the physical and mechanical parameters of the strata and rock mass in each area into Geo-Studio geotechnical analysis software, use Geo-Studio geotechnical analysis software to search for potential slip surfaces, analyze the morphology of potential slip surfaces, and determine the morphological parameters of the most dangerous slip surface; Strip column division module: Based on the morphological parameters of the most dangerous slip surface and the distribution of stratigraphic dip, the main slip surface is vertically divided into strip columns to obtain several micro-strip columns, and the three-dimensional mechanical effects of the slip surface at the bottom of the micro-strip columns are analyzed. Equivalent physical and mechanical parameters solution module: The three-dimensional mechanical effects of the bottom sliding surface of each micro-strip column are superimposed onto the main sliding surface of the sliding body, and the superimposed sliding force X and superimposed anti-slip force T of each row of micro-strip columns are solved in stages to determine the equivalent physical and mechanical parameters; Initial value setting module: Introduces equivalent physical and mechanical parameters into the Morgenstern-Price method, setting the unit weight λ and the slope stability coefficient F under the action of coal and rock composite. s The initial value is set according to the elasticity theory, and the expression of the function f(x) is set accordingly; Stability coefficient output module: outputs the set f(x) function and the assumptions λ and F. s Substitute the initial values ​​into the boundary conditions, continuously adjust the values ​​of the function f(x) and λ, and determine whether the boundary conditions are satisfied. Continue until the boundary conditions are satisfied, then output F. s .

Images