CN113239510A - Method for calculating optimal excavation slope rate of mine side slope - Google Patents
Method for calculating optimal excavation slope rate of mine side slope Download PDFInfo
- Publication number
- CN113239510A CN113239510A CN202110291622.4A CN202110291622A CN113239510A CN 113239510 A CN113239510 A CN 113239510A CN 202110291622 A CN202110291622 A CN 202110291622A CN 113239510 A CN113239510 A CN 113239510A
- Authority
- CN
- China
- Prior art keywords
- slope
- excavation
- mine
- optimal
- side slope
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000009412 basement excavation Methods 0.000 title claims abstract description 240
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000005457 optimization Methods 0.000 claims abstract description 19
- 230000005484 gravity Effects 0.000 claims abstract description 17
- 231100000817 safety factor Toxicity 0.000 claims description 29
- 230000009977 dual effect Effects 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000010008 shearing Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/23—Dune restoration or creation; Cliff stabilisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Operation Control Of Excavators (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for calculating the optimal excavation slope rate of a mine side slope, which takes the mine side slope as a research object, establishes a multi-objective optimization nonlinear mathematical programming model of the stability safety and the excavation economy of the mine side slope, and takes the maximum value of a safety coefficient and the maximum value of excavation benefit as objective functions at the same time; the excavation slope rate, the shearing force and the normal force of the structural surface are used as decision variables, the balance equation constraint condition of the side slope slide body and the yield condition of the structural surface are used as constraint conditions, and meanwhile, the gravity function constructed in the constraint conditions can change along with the change of the excavation slope rate, so that the constraint conditions of the limit balance of the side slope slide body can well relate to the safety and the economy of the mine side slope, and the actual mine side slope excavation requirements can be better met; the method can directly obtain the optimal excavation slope rate of the mine side slope, so that the stability and safety factor are the maximum while the excavation benefit is the highest, and further the unity of the excavation economy and the stability and safety is achieved.
Description
Technical Field
The invention relates to a method for calculating the optimal excavation slope rate of a mine side slope, and belongs to the technical field of analysis of stability of the mine side slope.
Background
In recent years, the problem of slope stability caused by excavation of mine slopes is receiving wide attention of a great number of geotechnical engineers. Two problems are often involved in performing mine slope excavation: firstly, the stability and safety of the mine side slope after excavation and secondly, the excavation economy caused by the excavation volume. The stability and safety of the excavated mine side slope mainly depend on the property and spatial combination of the structural surface, and the bedding surface, the weak interlayer, the fault and the joint surface in the mine side slope are all main reasons for causing instability and damage of the side slope in the excavation process; larger excavation volumes result in higher excavation economics and vice versa. In actual engineering, excavation slope rate design is usually performed by adopting an excavation slope rate method, which is mainly used for determining the excavation slope rate of an excavation slope according to the suggestions of geological specialties and similar engineering experiences, such as: the recommended excavation slope rate of the mine slope is given by the technical Specification of building slope engineering (GB 50330-2013). The proposed excavation slope rate is an empirical value, and the safety factor and the excavation benefit of the mine side slope cannot be guaranteed to be optimal at the same time, namely the purposes of safety and economy are achieved.
The current mine side slope excavation design has the following problems:
(1) in actual engineering, the excavation slope rate of a side slope is determined according to the suggestions of geological specialties and similar engineering experiences, and relevant principles and theories are still to be studied;
(2) the mine side slope suggested excavation slope rate given by the specification is adopted, and the stability, the safety and the excavation economy are not guaranteed to be optimal;
(3) the mine side slope excavation design is a multi-target problem of seeking the mine side slope excavation benefit as large as possible and seeking the mine side slope safety factor as high as possible at the same time, and no good method is considered comprehensively at present;
based on the current situation, a method for calculating the optimal excavation slope rate of the mine side slope is needed to satisfy stable safety and excavation economy.
Disclosure of Invention
The invention provides a method for calculating the optimal excavation slope rate of a mine side slope, which can obtain the optimal stable safety coefficient, optimal excavation benefit and optimal excavation slope rate which simultaneously meet the requirements of stable safety and optimal excavation economy.
The technical scheme of the invention is as follows: a method for calculating the optimal excavation slope rate of a mine side slope comprises the following steps:
step 1, drawing up geometric parameters of a drawn-up mine side slope;
step 2, drawing up physical and mechanical parameters of the mine side slope;
step 3, analyzing the stress of the mine slope slide body;
step 4, establishing a stable safety objective function and an excavation economy objective function of the mine side slope;
step 5, establishing a constraint condition of the limit balance of the mine slope landslide;
step 6, establishing a multi-objective optimization nonlinear mathematical programming model of mine slope stability safety and excavation economy;
and 7, solving a multi-objective optimization nonlinear mathematical programming model of mine slope stability safety and excavation economy, and obtaining stability safety factor, excavation benefit and optimal excavation slope rate when the stability safety and the excavation economy are optimal.
The drawing up of the geometric parameters of the drawn-up mine slope comprises the following steps: height of mine side slope H and distance H from bottom end of structural surface to x axis1(ii) a The natural slope rate beta of the mine side slope; and (4) the excavation slope rate theta of the mine side slope.
The physical and mechanical parameters of the drawn up mine slope comprise: volume weight gamma of mine slope slide, inclination angle alpha of structural surface, cohesion c of structural surface and friction angle of structural surface
The stress analysis of the mine side slope slide body is specifically as follows: gravity G acting on the centroid of the mine slope slide body, and shear force F acting on the mine slope slide body structural surfaceSAnd normal force FNIn which F isSTurning in the counterclockwise direction negative, FNThe pressure is positive.
The establishment of the stable safety objective function and the excavation economy objective function of the mine side slope specifically comprises the following steps:
firstly, establishing a stable safety objective function of a mine side slope
Taking the stability safety coefficient of the mine slope as a stability safety objective function, seeking the maximum value of the stability safety coefficient of the mine slope, and establishing the stability safety objective function as follows:
Maximize:K
in the formula: maximize means "max"; k is the mine slope stability safety coefficient;
establishing an excavation economy objective function of the mine side slope
Taking the excavation benefit, which is the product of the price of the unit volume of the excavation slope sliding body and the excavation volume, as an excavation economic objective function, seeking the maximum value of the excavation benefit, and establishing the objective function of the excavation economic as follows:
Maximize:J
in the formula: maximize means "max"; j is mine side slope excavation benefit, and when a side slope sliding body with unit width is taken for calculation:
in the formula: k is a radical of1The unit volume selling price of the mine side slope sliding body is H, the height of the mine side slope, theta, the excavation slope rate of the mine side slope and beta, the natural slope rate of the mine side slope.
The constraint conditions for establishing the limit balance of the mine slope landslide body are as follows:
the constraint conditions of the limit balance of the mine slope landslide body comprise: the balance equation constraint condition of the mine slope landslide body and the yield condition of the structural plane;
balance equation constraint conditions of mine slope landslide:
the balance equation of the mine slope sliding body in the horizontal direction is established as follows:
FS cosα-FN sinα=0
the balance equation in the vertical direction of the mine slope sliding body is established as follows:
G+FS sinα+FN cosα=0
wherein:g is the gravity acting on the mine slope slide centroid, FSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, wherein FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, and alpha rotates in the counterclockwise direction to be positive; h is the height of the mine side slope1The distance from the bottom end of the structural plane to the x axis is used, and theta is the excavation slope rate of the mine side slope; gamma is the volume weight of the mine slope sliding body;
secondly, the yield condition of the mine side slope slide body structure surface is as follows:
in the formula: k is the mine slope stability safety coefficient; c.respectively the cohesion and friction angle of the structural surface; l is the length of the structural plane.
The establishment of the multi-objective optimization nonlinear mathematical programming model for the mine slope stability safety and the excavation economy is as follows:
the method comprises the following steps of establishing a multi-objective optimization nonlinear mathematical programming model of mine slope stability safety and excavation economy, which takes mine slope stability safety factors and excavation benefits as dual objective function and simultaneously meets constraint conditions of mine slope slide mass balance equations and yield conditions of structural surfaces, as follows:
in the formula: k is the mine slope stability safety coefficient; j, mine side slope excavation benefit; fSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, wherein FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, and alpha rotates in the counterclockwise direction to be positive; g is the gravity acting on the centroid of the mine slope slide body; c.respectively the cohesion and friction angle of the structural surface; l is the length of the structural plane; k is a radical of1The unit volume selling price of the mine slope slide body is obtained; h is the height of the mine side slope; theta is the excavation slope rate of the mine side slope; beta is the natural slope rate of the mine side slope.
Solving a multi-objective optimization nonlinear mathematical programming model of mine slope stability safety and excavation economy to obtain stability safety factor, excavation benefit and optimal excavation slope rate when the stability safety and the excavation economy are optimal, specifically:
solving a nonlinear mathematical programming model with optimal stability and safety by adopting a penalty function method, wherein the nonlinear mathematical programming model with optimal stability and safety is as follows:
in the formula:the stability and safety factor of the side slope when the excavation stability and safety are optimal is obtained; fSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, wherein FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, and alpha rotates in the counterclockwise direction to be positive; g is the gravity acting on the centroid of the mine slope slide body; c.respectively the cohesion and friction angle of the structural surface; l is the length of the structural plane; k is a radical of1The unit volume selling price of the mine slope slide body is obtained; h is the height of the mine side slope; j. the design is a square1The side slope excavation benefit is achieved when the excavation stability and safety are optimal; thetaKExcavating slope rate for the side slope with optimal excavation stability and safety; beta is the natural slope rate of the mine side slope;
solving the nonlinear mathematical programming model with the optimal excavation economy by adopting a penalty function method, wherein the nonlinear mathematical programming model with the optimal excavation economy is as follows:
in the formula:the side slope excavation benefit is achieved when the excavation economy is optimal; k1The safety coefficient is stabilized for the side slope when the excavation economy is optimal; thetaJExcavating slope rate for the side slope with optimal excavating economy;
and thirdly, constructing an objective function with optimal stability, safety and excavation economy as follows:
in the formula: z is an objective function when the stability safety and the excavation economy are optimal;the method is used for realizing the slope excavation benefit when the stability safety and the excavation economy are optimal;the slope stability safety factor is the slope stability safety factor when the stability safety and the excavation economy are both optimal;
solving the nonlinear mathematical programming model of the mine slope excavation economy and the stable safety by adopting a penalty function method, wherein the nonlinear mathematical programming model of the mine slope excavation economy and the stable safety is as follows:
in the formula: thetaoptThe optimal excavation slope rate is achieved when the stability, the safety and the excavation economy are optimal.
The invention has the beneficial effects that:
the method takes the mine side slope as a research object, establishes a multi-objective optimization nonlinear mathematical programming model of the stability safety and the excavation economy of the mine side slope, and simultaneously takes the maximum value of the safety coefficient and the maximum value of the excavation benefit as objective functions; the excavation slope rate, the shearing force and the normal force of the structural surface are used as decision variables, the balance equation constraint condition of the side slope slide body and the yield condition of the structural surface are used as constraint conditions, and meanwhile, the gravity function constructed in the constraint conditions can change along with the change of the excavation slope rate, so that the constraint conditions of the limit balance of the side slope slide body can well relate to the safety and the economy of the mine side slope, and the actual mine side slope excavation requirements can be better met; the model is different from the traditional method that the model is established by combining theory and practice according to the suggestions of geological specialties and similar engineering experiences, and the optimal excavation slope rate of the mine side slope can be directly obtained by the method, so that the stable safety factor is the maximum while the excavation benefit is the highest, and the unification of the excavation economy and the stable safety is further achieved.
The method is suitable for calculating the optimal excavation slope rate of the mine with comprehensive consideration of stable safety and excavation economy, and has the characteristics of clear concept, high calculation precision and simple and convenient engineering application.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a diagram of an excavation area of a mine slope according to the present invention;
FIG. 3 is a schematic view of a mine slope of the present invention;
FIG. 4 is a schematic view of a mine slope according to an embodiment of the present invention;
fig. 5 is a schematic view of the force applied to the side slope sliding body according to the embodiment of the invention.
Detailed Description
Example 1: as shown in fig. 1 to 5, a method for calculating an optimal excavation slope rate of a mine side slope includes:
step 1, drawing up geometric parameters of a drawn-up mine side slope;
step 2, drawing up physical and mechanical parameters of the mine side slope;
step 3, analyzing the stress of the mine slope slide body;
step 4, establishing a stable safety objective function and an excavation economy objective function of the mine side slope;
step 5, establishing a constraint condition of the limit balance of the mine slope landslide;
step 6, establishing a multi-objective optimization nonlinear mathematical programming model of mine slope stability safety and excavation economy;
step 7, solving a multi-objective optimization nonlinear mathematical programming model of mine slope stability safety and excavation economy, and obtaining a stability safety coefficient when the stability safety and the excavation economy are optimalBenefits of excavationAnd optimum excavation slope rate θopt。
Further, the invention provides the following implementation processes:
(1) drawing up the geometric parameters of mine side slope
As shown in fig. 4, the mine slope of the embodiment is established by a coordinate system: and (3) taking an angular point A on the outer side of the slope top of the mine slope as a vertical line, intersecting the slope bottom with a point O, taking a straight line where BO is located as an x-axis, taking the x-axis as a right positive, taking a straight line where OA is located as a y-axis, and taking the y-axis as a positive, and establishing a coordinate system as shown in the figure. Drawing up the geometric parameters of the side slope according to the actual situation of the mine side slope, comprising the following steps: the height H of the mine side slope is 50 m; the distance from the bottom end of the structural surface to the x axis is 8.49 m; AB is a natural face of the mine side slope, and the natural slope rate beta is 60 degrees; OC is the face-to-air surface after the excavation of the mine side slope, and the excavation slope rate is theta.
(2) Drawing up physical and mechanical parameters of mine side slope
The physical and mechanical parameters of the side slope are drawn up according to the actual situation of the mine side slope, and the drawing comprises the following steps: the volume weight gamma of the side slope slide is 19.00kN/m3The inclination angle of the structural plane is alpha, the cohesion c of the structural plane is 40kPa, and the friction angle of the structural planeIs 40 deg..
(3) Stress analysis of mine slope slide
As shown in fig. 5, which is a schematic diagram of a mine slope slide body, the stress analysis is performed on the slope slide body, the centroid of the slope slide body acts on the gravity G in the y direction, and the shearing force F acts on the structural surface of the slope slide bodySAnd normal force FNIn which F isSTurning in the counterclockwise direction negative, FNThe pressure is positive.
(4) Establishing stable safety objective function and excavation economy objective function of mine side slope
Firstly, establishing a stable safety objective function of a mine side slope
The obtained strength reserve coefficient is the stability safety coefficient, and the stability safety coefficient is higher. Taking the mine slope stability safety coefficient as a stability safety objective function, seeking the maximum value of the mine slope stability safety coefficient, and establishing the stability safety objective function as follows:
Maximize:K (1)
in the formula: maximize means "max"; k is the mine slope stability safety coefficient.
Establishing an excavation economy objective function of the mine side slope
The excavation economy of mine side slope embodies through the size of excavation volume, and the excavation volume is big more, and the excavation economy is higher, will excavate the product of the price of side slope landslide unit volume and excavation volume, excavate the benefit promptly as excavation economy objective function to seek the maximum value of excavation benefit, then the objective function of excavation economy can write:
Maximize:J (2)
in the formula: maximize means "max"; j, the slope excavation benefit;
fig. 2 is a mine slope excavation area diagram, where the excavation benefit is the product of the unit volume selling price of the slope sliding body and the excavation volume, and when the slope sliding body with unit width is taken for calculation:
in the formula: k is a radical of1The unit volume selling price of the side slope sliding body is 1000, and theta is the excavation slope rate of the mine side slope;
(5) constraint condition for establishing limit balance of mine slope landslide
The constraint conditions of the limit balance of the mine slope landslide body comprise: the balance equation constraint condition of the mine slope landslide body and the yield condition of the structural plane;
the mine side slope slide body balance equation constraint conditions are as follows:
the mine slope slide body is kept balanced by gravity on the centroid, shearing force on the structural surface and normal force on the structural surface.
The balance equation of the mine slope sliding body in the horizontal direction is established as follows:
FS cosα-FN sinα=0 (4)
the balance equation in the vertical direction of the mine slope sliding body is established as follows:
G+FS sinα+FN cosα=0 (5)
wherein:g is the gravity of the slide centroid of the side slope along the y direction, FSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, wherein FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, alpha rotates in the counterclockwise direction as positive, and gamma is the volume weight of the mine slope slide.
Secondly, the yield condition of the mine side slope slide body structure surface is as follows:
assuming that the mine slope sliding body is a rigid body block which cannot generate any damage, the damage of the mine slope only occurs on a structural plane, the structural plane meets the Mohr-Coulomb yield condition, and the yield condition on the structural plane is written as follows:
K|FS|-FN tan40°-40×l≤0 (6)
in the formula: fSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, wherein FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, and alpha rotates in the counterclockwise direction to be positive; l is the length of the structural plane; c.respectively the cohesion and friction angle of the structural surface; k is the mine slope stability safety coefficient.
(6) Multi-objective optimization nonlinear mathematical programming model for establishing mine slope stability safety and excavation economy
The method comprises the following steps of establishing a multi-objective optimization nonlinear mathematical programming model of mine slope stability safety and excavation economy, which takes mine slope stability safety factors and excavation benefits as dual objective function and simultaneously meets constraint conditions of mine slope slide mass balance equations and yield conditions of structural surfaces, as follows:
in the formula: k is a slope stability safety coefficient; j, the slope excavation benefit; theta is the slope rate of slope excavation(ii) a G is the gravity acting on the centroid of the side slope slide body; fSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, wherein FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, and alpha rotates in the counterclockwise direction to be positive; l is the length of the structural plane; c.respectively, the cohesion and the friction angle of the structural surfaces.
The formula (7) is a multi-target nonlinear mathematical programming model, and the stability safety factor when the stability safety and the excavation economy are optimal cannot be obtained simultaneouslyBenefits of excavationAnd optimum excavation slope rate θoptTherefore, the invention firstly solves the nonlinear mathematical programming model with optimal stability and safety to obtainThen solving a nonlinear mathematical programming model with optimal excavation economy to obtainAnd then according to the solutionConstructing an objective function with optimal stability safety and optimal excavation economy so as to obtain the stability safety coefficient when the stability safety and the excavation economy are optimalBenefits of excavationAnd optimum excavation slope rate θopt(ii) a In this wayThe deficiency that ideal solutions cannot be obtained simultaneously is avoided (an ideal solution can be determined for a nonlinear mathematical programming model with optimal excavation economy)An ideal solution can be determined for a nonlinear mathematical programming model that optimizes stability and safetyGenerally speaking of'i≠α″iThis makes it impossible to obtain ideal solutions of the multi-objective functions simultaneously, and therefore by solving them according toConstructing a target function with optimal stability safety and excavation economy, thereby realizing finding an alphaiSo that J (alpha)i) Andand K (alpha)i) Andas close as possible). The method comprises the following specific steps:
(7) solving a nonlinear mathematical programming model for multi-objective optimization of mine slope stability safety and excavation economy to obtain a stability safety coefficient when the stability safety and the excavation economy are both optimalBenefits of excavationAnd optimum excavation slope rate θopt。
step 1: nonlinear mathematical programming model for solving optimal stability and safety
According to the target function formula (1) and the constraint condition formulas (3), (4), (5) and (6), obtaining a nonlinear mathematical programming model with optimal stability and safety of mine slope excavation:
in the formula:the stability and safety factor of the side slope when the excavation stability and safety are optimal is obtained; j. the design is a square1The side slope excavation benefit is achieved when the excavation stability and safety are optimal; thetaKExcavating slope rate for the side slope with optimal excavation stability and safety; g is the gravity of the slide centroid of the side slope along the y direction; fSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, wherein FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, and alpha rotates in the counterclockwise direction to be positive; l is the length of the structural plane.
The invention discloses a non-linear mathematical programming model for solving the optimal safety of mine slope stability in a formula (8) by adopting a penalty function method.
step 2: nonlinear mathematical programming model for solving optimal excavation economy
According to the target function formula (2) and the constraint condition formulas (3), (4), (5) and (6), obtaining a nonlinear mathematical programming model with optimal mine slope excavation economy as follows:
in the formula: k1The safety coefficient is stabilized for the side slope when the excavation economy is optimal;the side slope excavation benefit is achieved when the excavation economy is optimal; thetaJExcavating slope rate for the side slope with optimal excavating economy; g is the gravity of the slide centroid of the side slope along the y direction; fSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, whichMiddle FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, and alpha rotates in the counterclockwise direction to be positive; l is the length of the structural plane.
And (3) solving the nonlinear mathematical programming model with the optimal economic performance of the mine slope excavation by adopting a penalty function method (9).
step 3: objective function with optimal construction stability safety and optimal excavation economy
The objective function with optimal construction stability safety and excavation economy is as follows:
in the formula: z is an objective function when the stability safety and the excavation economy are both optimal,for the side slope excavation benefit when the excavation economy is optimal,the method is used for realizing the slope excavation benefit when the stability safety and the excavation economy are optimal;for the stability and safety factor of the side slope when the stability and safety of excavation are optimal,the slope stability safety factor is the slope stability safety factor when the stability safety and the excavation economy are both optimal;
step 4: nonlinear mathematical programming model for establishing mine slope excavation economy and stable safety
Converting the nonlinear mathematical programming model of the multi-objective optimization of mine slope stability safety and excavation economy into a nonlinear mathematical programming model of single-objective optimization represented by the following formula according to the objective function formula (10) and the constraint condition formulas (3), (4), (5) and (6):
in the formula: z is an objective function when the stability safety and the excavation economy are both optimal,for the side slope excavation benefit when the excavation economy is optimal,the method is used for realizing the slope excavation benefit when the stability safety and the excavation economy are optimal;for the stability and safety factor of the side slope when the stability and safety of excavation are optimal,the slope stability safety factor is the slope stability safety factor when the stability safety and the excavation economy are both optimal; thetaoptThe optimal excavation slope rate is obtained when the stability safety and the excavation economy are optimal; g is the gravity of the slide centroid of the side slope along the y direction; fSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, wherein FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, and alpha rotates in the counterclockwise direction to be positive; l is the length of the structural plane.
And (3) solving a nonlinear mathematical programming model of the mine slope excavation economy and stable safety in the formula (14) by adopting a penalty function method.
The calculation results obtained when the excavation economy is optimal under different structural surface inclination angles alpha are listed in the table 1; the results of calculations obtained when the stability and safety are optimal for different structural surface inclination angles α are listed in table 2; the stability factor obtained with the process of the invention is shown in Table 3Benefits of excavationAnd optimum excavation slope rate θopt(ii) a Table 4 compares the stability safety factor and excavation benefit obtained by the method of the present invention when the excavation economy is optimal under different structural surface inclination angles α; table 5 compares the optimum stability and safety for different structural surface inclination angles alpha with the stability and safety factor and excavation benefit obtained by the method of the present invention.
TABLE 1 example calculation results for mine slope when excavation economics are optimal
Table 2 example calculation results for mine slope when stability and safety are optimal
Table 3 example mine slope calculated results obtained by the method of the invention
TABLE 4 comparison of excavation economics optimal with stability safety factor and excavation benefit of the method of the present invention
TABLE 5 comparison of optimal stability and safety with stability, safety factor and excavation benefit of the method of the invention
As shown in table 4, the excavation efficiency obtained by the method of the present invention was found to be equal to 45 ° in the structural plane inclination angle αCompared with the excavation benefit when the excavation economy is optimalThe stable safety factor can be realized only by reducing 2.00 percentCompared with stable safety factor K when the excavation economy is optimal1The increase of 56.84 percent achieves the unification of stable safety and economy.
As shown in Table 5, the stability factor obtained by the process according to the invention is higher for a structural surface inclination α of 50 °Stability factor of safety compared to stability factor of safety when optimalThe excavation benefit can be realized only by reducing 0.36 percentExcavation benefit J compared with stable safety when optimal1The increase of 2.82 percent, and the unification of safety and economy is achieved.
While the present invention has been described in detail with reference to the embodiments shown in the drawings, the present invention is not limited to the embodiments, and various changes can be made without departing from the spirit of the present invention within the knowledge of those skilled in the art.
Claims (8)
1. A method for calculating the optimal excavation slope rate of a mine side slope is characterized by comprising the following steps: the method comprises the following steps:
step 1, drawing up geometric parameters of a drawn-up mine side slope;
step 2, drawing up physical and mechanical parameters of the mine side slope;
step 3, analyzing the stress of the mine slope slide body;
step 4, establishing a stable safety objective function and an excavation economy objective function of the mine side slope;
step 5, establishing a constraint condition of the limit balance of the mine slope landslide;
step 6, establishing a multi-objective optimization nonlinear mathematical programming model of mine slope stability safety and excavation economy;
and 7, solving a multi-objective optimization nonlinear mathematical programming model of mine slope stability safety and excavation economy, and obtaining stability safety factor, excavation benefit and optimal excavation slope rate when the stability safety and the excavation economy are optimal.
2. The method for calculating the optimal excavation slope rate of the mine side slope according to claim 1, characterized in that: the drawing up of the geometric parameters of the drawn-up mine slope comprises the following steps: height of mine side slope H and distance H from bottom end of structural surface to x axis1(ii) a The natural slope rate beta of the mine side slope; and (4) the excavation slope rate theta of the mine side slope.
3. The method for calculating the optimal excavation slope rate of the mine side slope according to claim 1, characterized in that: the physical and mechanical parameters of the drawn up mine slope comprise: volume weight gamma of mine slope slide, inclination angle alpha of structural surface, cohesion c of structural surface and friction angle of structural surface
4. The method for calculating the optimal excavation slope rate of the mine side slope according to claim 1, characterized in that: the stress analysis of the mine side slope slide body is specifically as follows: gravity G acting on the centroid of the mine slope slide body, and shear force F acting on the mine slope slide body structural surfaceSAnd normal force FNIn which F isSTurning in the counterclockwise direction negative, FNThe pressure is positive.
5. The method for calculating the optimal excavation slope rate of the mine side slope according to claim 1, characterized in that: the establishment of the stable safety objective function and the excavation economy objective function of the mine side slope specifically comprises the following steps:
firstly, establishing a stable safety objective function of a mine side slope
Taking the stability safety coefficient of the mine slope as a stability safety objective function, seeking the maximum value of the stability safety coefficient of the mine slope, and establishing the stability safety objective function as follows:
Maximize:K
in the formula: maximize means "max"; k is the mine slope stability safety coefficient;
establishing an excavation economy objective function of the mine side slope
Taking the excavation benefit, which is the product of the price of the unit volume of the excavation slope sliding body and the excavation volume, as an excavation economic objective function, seeking the maximum value of the excavation benefit, and establishing the objective function of the excavation economic as follows:
Maximize:J
in the formula: maximize means "max"; j is mine side slope excavation benefit, and when a side slope sliding body with unit width is taken for calculation:
in the formula: k is a radical of1The unit volume selling price of the mine side slope sliding body is H, the height of the mine side slope, theta, the excavation slope rate of the mine side slope and beta, the natural slope rate of the mine side slope.
6. The method for calculating the optimal excavation slope rate of the mine side slope according to claim 1, characterized in that: the constraint conditions for establishing the limit balance of the mine slope landslide body are as follows:
the constraint conditions of the limit balance of the mine slope landslide body comprise: the balance equation constraint condition of the mine slope landslide body and the yield condition of the structural plane;
balance equation constraint conditions of mine slope landslide:
the balance equation of the mine slope sliding body in the horizontal direction is established as follows:
FS cosα-FN sinα=0
the balance equation in the vertical direction of the mine slope sliding body is established as follows:
G+FS sinα+FN cosα=0
wherein:g is the gravity acting on the mine slope slide centroid, FSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, wherein FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, and alpha rotates in the counterclockwise direction to be positive; h is the height of the mine side slope1The distance from the bottom end of the structural plane to the x axis is used, and theta is the excavation slope rate of the mine side slope; gamma is the volume weight of the mine slope sliding body;
secondly, the yield condition of the mine side slope slide body structure surface is as follows:
7. The method for calculating the optimal excavation slope rate of the mine side slope according to claim 1, characterized in that: the establishment of the multi-objective optimization nonlinear mathematical programming model for the mine slope stability safety and the excavation economy is as follows:
the method comprises the following steps of establishing a multi-objective optimization nonlinear mathematical programming model of mine slope stability safety and excavation economy, which takes mine slope stability safety factors and excavation benefits as dual objective function and simultaneously meets constraint conditions of mine slope slide mass balance equations and yield conditions of structural surfaces, as follows:
in the formula: k is the mine slope stability safety coefficient; j, mine side slope excavation benefit; fSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, wherein FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, and alpha rotates in the counterclockwise direction to be positive; g is the gravity acting on the centroid of the mine slope slide body; c.respectively the cohesion and friction angle of the structural surface; l is the length of the structural plane; k is a radical of1The unit volume selling price of the mine slope slide body is obtained; h is the height of the mine side slope; theta is the excavation slope rate of the mine side slope; beta is the natural slope rate of the mine side slope.
8. The method for calculating the optimal excavation slope rate of the mine side slope according to claim 1, characterized in that: solving a multi-objective optimization nonlinear mathematical programming model of mine slope stability safety and excavation economy to obtain stability safety factor, excavation benefit and optimal excavation slope rate when the stability safety and the excavation economy are optimal, specifically:
solving a nonlinear mathematical programming model with optimal stability and safety by adopting a penalty function method, wherein the nonlinear mathematical programming model with optimal stability and safety is as follows:
in the formula:for excavation stabilizingThe slope stability safety coefficient when the comprehensiveness is optimal; fSShear forces acting on structural surfaces, in which FSNegative in counterclockwise rotation; fNIs a normal force acting on a structural surface, wherein FNThe pressure is taken as positive; alpha is a positive included angle between the structural plane and the x axis, and alpha rotates in the counterclockwise direction to be positive; g is the gravity acting on the centroid of the mine slope slide body; c.respectively the cohesion and friction angle of the structural surface; l is the length of the structural plane; k is a radical of1The unit volume selling price of the mine slope slide body is obtained; h is the height of the mine side slope; j. the design is a square1The side slope excavation benefit is achieved when the excavation stability and safety are optimal; thetaKExcavating slope rate for the side slope with optimal excavation stability and safety; beta is the natural slope rate of the mine side slope;
solving the nonlinear mathematical programming model with the optimal excavation economy by adopting a penalty function method, wherein the nonlinear mathematical programming model with the optimal excavation economy is as follows:
in the formula:the side slope excavation benefit is achieved when the excavation economy is optimal; k1The safety coefficient is stabilized for the side slope when the excavation economy is optimal; thetaJExcavating slope rate for the side slope with optimal excavating economy;
and thirdly, constructing an objective function with optimal stability, safety and excavation economy as follows:
in the formula: z is an objective function when the stability safety and the excavation economy are optimal;the method is used for realizing the slope excavation benefit when the stability safety and the excavation economy are optimal;the slope stability safety factor is the slope stability safety factor when the stability safety and the excavation economy are both optimal;
solving the nonlinear mathematical programming model of the mine slope excavation economy and the stable safety by adopting a penalty function method, wherein the nonlinear mathematical programming model of the mine slope excavation economy and the stable safety is as follows:
in the formula: thetaoptThe optimal excavation slope rate is achieved when the stability, the safety and the excavation economy are optimal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110291622.4A CN113239510B (en) | 2021-03-18 | 2021-03-18 | Method for calculating optimal excavation slope rate of mine side slope |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110291622.4A CN113239510B (en) | 2021-03-18 | 2021-03-18 | Method for calculating optimal excavation slope rate of mine side slope |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113239510A true CN113239510A (en) | 2021-08-10 |
CN113239510B CN113239510B (en) | 2022-05-20 |
Family
ID=77130322
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110291622.4A Active CN113239510B (en) | 2021-03-18 | 2021-03-18 | Method for calculating optimal excavation slope rate of mine side slope |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113239510B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102002955A (en) * | 2010-11-09 | 2011-04-06 | 河海大学 | Slope deformation monitoring and arranging method |
CN102789516A (en) * | 2012-05-21 | 2012-11-21 | 河北钢铁集团矿业有限公司 | Stability numerical analysis and optimization design method based on monitoring information in slope construction process |
CN106049511A (en) * | 2016-06-20 | 2016-10-26 | 中钢集团马鞍山矿山研究院有限公司 | Combined reinforcement method suitable for non-coal strip mine side slope landslide control |
KR101894495B1 (en) * | 2017-08-17 | 2018-09-04 | 서울대학교산학협력단 | Prediction Method and System of Disaster of Debris-flow |
CN110688747A (en) * | 2019-09-19 | 2020-01-14 | 昆明理工大学 | Rock slope anchoring calculation method based on safety and economy multi-objective optimization |
CN111368458A (en) * | 2020-03-23 | 2020-07-03 | 青岛理工大学 | Method for calculating foundation pit excavation slope safety coefficient |
CN111611728A (en) * | 2020-06-29 | 2020-09-01 | 中国有色金属工业昆明勘察设计研究院有限公司 | Stability calculation method for Xigeda stratum rock slope trailing edge slip fracture surface position |
CN111912953A (en) * | 2020-07-31 | 2020-11-10 | 青岛理工大学 | Deep-well mining slope stability determination method based on excavation amount monitoring |
-
2021
- 2021-03-18 CN CN202110291622.4A patent/CN113239510B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102002955A (en) * | 2010-11-09 | 2011-04-06 | 河海大学 | Slope deformation monitoring and arranging method |
CN102789516A (en) * | 2012-05-21 | 2012-11-21 | 河北钢铁集团矿业有限公司 | Stability numerical analysis and optimization design method based on monitoring information in slope construction process |
CN106049511A (en) * | 2016-06-20 | 2016-10-26 | 中钢集团马鞍山矿山研究院有限公司 | Combined reinforcement method suitable for non-coal strip mine side slope landslide control |
KR101894495B1 (en) * | 2017-08-17 | 2018-09-04 | 서울대학교산학협력단 | Prediction Method and System of Disaster of Debris-flow |
CN110688747A (en) * | 2019-09-19 | 2020-01-14 | 昆明理工大学 | Rock slope anchoring calculation method based on safety and economy multi-objective optimization |
CN111368458A (en) * | 2020-03-23 | 2020-07-03 | 青岛理工大学 | Method for calculating foundation pit excavation slope safety coefficient |
CN111611728A (en) * | 2020-06-29 | 2020-09-01 | 中国有色金属工业昆明勘察设计研究院有限公司 | Stability calculation method for Xigeda stratum rock slope trailing edge slip fracture surface position |
CN111912953A (en) * | 2020-07-31 | 2020-11-10 | 青岛理工大学 | Deep-well mining slope stability determination method based on excavation amount monitoring |
Non-Patent Citations (2)
Title |
---|
张小艳 等: "基于塑性极限分析上限法理论的土质边坡可靠度分析", 《岩土力学》 * |
陈秋茹,梁敏: "挖方边坡优化设计方法研究", 《重庆交通学院学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113239510B (en) | 2022-05-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107330146B (en) | Rock slope ultimate bearing capacity analysis upper limit method considering rotation effect | |
CN110688747B (en) | Rock slope anchoring calculation method based on safety and economy multi-objective optimization | |
CN111611728B (en) | Stability calculation method for Xigeda stratum rock slope trailing edge slip fracture surface position | |
CN110083883B (en) | Three-dimensional slope safety coefficient solving method based on limit analysis lower limit theorem | |
CN103775080A (en) | Construction method of attitude angle regulating model of manless working face coalcutter for thin seam | |
CN108446251B (en) | Method for calculating active rock pressure of rock slope | |
CN103266617A (en) | Method for computing optimal anchoring angle of rock slope wedge | |
CN113239510B (en) | Method for calculating optimal excavation slope rate of mine side slope | |
CN108170898B (en) | Lower limit method for reliability analysis of jointed rock slope | |
CN112069574A (en) | Circular foundation pit stability analysis method based on limit analysis upper limit theorem | |
Kongshu et al. | Force transmission characteristics for the non-equidistant arrangement thrust systems of shield tunneling machines | |
CN109214050A (en) | A kind of polar region deviation of plumb line is without the detailed computation model of singularity and its modeling method | |
CN106815445A (en) | A kind of shallow buried coal seam section coal pillar LOAD FOR model and its computational methods | |
CN108763833B (en) | Method for calculating deflection of foundation pit supporting pile in consideration of soil resistance sudden change | |
Wang et al. | Research on the calculation model and control method of initial supporting force for temporary support in the underground excavation roadway of coal mine | |
CN101761136B (en) | Structure joint simultaneously meeting constraint requirements of space in different directions | |
CN103761367B (en) | Displacement control objective based analytic calculating method of string beam structural internal force | |
CN113094905B (en) | Calculation method and system suitable for multi-middle-section continuous empty area support key points | |
Huang et al. | Study of prevention and control technology for roadway excavation under the soft and extra-thick coal roof in luling coal mine | |
CN113486500B (en) | Method for obtaining optimal excavation load reduction depth of side slope | |
CN113468766A (en) | Mechanical characteristic research and optimization design for jacking construction of rectangular curve jacking pipe | |
CN203411981U (en) | Foundation pit horizontal support component deformation control structure | |
Peng et al. | Calculation of passive earth pressure on rigid retaining wall affected by seepage and soil arching | |
CN102194251B (en) | Technology for controlling deformation of slope rock mass | |
CN114970122B (en) | Method and system for calculating unsteady optimal supporting force of shield tunnel driving surface |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |