CN113674406A - Three-dimensional simulation method for surface rock mass of open-air deep hole bench blasting pile - Google Patents
Three-dimensional simulation method for surface rock mass of open-air deep hole bench blasting pile Download PDFInfo
- Publication number
- CN113674406A CN113674406A CN202110266412.XA CN202110266412A CN113674406A CN 113674406 A CN113674406 A CN 113674406A CN 202110266412 A CN202110266412 A CN 202110266412A CN 113674406 A CN113674406 A CN 113674406A
- Authority
- CN
- China
- Prior art keywords
- pile
- sphere
- rock
- blasting
- dimensional simulation
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application discloses a three-dimensional simulation method for blasting surface rock masses of a blasting pile in an open-air deep hole step, which relates to the technical field of blasting effect evaluation and comprises determination of the shape of the blasting pile rock masses and the three-dimensional shape of the blasting pile, division and block size statistics of local areas on the surface of the blasting pile, three-dimensional simulation of the rock masses on the surface of the blasting pile and three-dimensional simulation results of the rock masses on the surface of the blasting pile. According to the method, the surface of the detonation pile is partitioned, and the rock mass of each area is subjected to statistical analysis to obtain a detonation pile block size grading accumulation curve and a large block rate of the area, so that a detonation pile surface block size local distribution model is formed; the model and the information of the explosive pile form are combined, an explosive pile surface rock block generation algorithm is formed by improving an inward arrangement scheme in a discrete unit construction algorithm, and the three-dimensional model simulation function of the explosive pile form and the surface rock block is realized by adopting a three-dimensional visualization technology.
Description
Technical Field
The application relates to the technical field of blasting effect evaluation, in particular to a three-dimensional simulation method for blasting surface rock masses of a blasting pile in an open-air deep hole step.
Background
The research on the distribution rule of the open-air step deep hole blasting lumpiness and the visual expression thereof is an important means for evaluating the blasting effect and is also an important basis for feeding back and improving the open-air step deep hole blasting design and construction method. The detection method of the open-air deep hole step blasting block degree is mainly surface block degree detection, the internal block degree is obtained by an indirect method, the research on the open-air deep hole blasting block degree distribution rule mostly focuses on the research on the surface block degree distribution rule, and the main research contents and methods comprise: (1) researching a detection method and a rapid identification method of the bench blasting blockiness of the open-air deep hole; (2) and (3) researching a mathematical model between the distribution rule of the blasting lumpiness of the open-air deep hole step and blasting parameters. Such as: the method comprises the steps that a Gene Expression Programming (GEP) algorithm is selected by Nie army et al of the university of south China, and a blasting blockiness prediction model based on GEP is established; the Schuqizhi et al of the university of Zhongnan establishes a rock blasting blockiness prediction model based on LS-SVR, and predicts the mine rock blasting blockiness under the condition of small samples; the Wanren super-et al at Tianjin university establishes a blasting blockiness neural network prediction model based on an LM algorithm.
The main problems with these research methods and contents are: (1) the mathematical model of the open-air deep hole step blasting blockiness distribution rule is too abstract, which is not beneficial to the mastering of technicians and the popularization and application of achievements; (2) the mathematical model of the open-air deep hole step blasting block size distribution rule is a statistical mathematical model, reflects an average grading model of the surface block size of the blasting pile, and cannot reflect the actual distribution condition of the surface block size of the blasting pile. In open-air deep hole bench blasting, because the engineering rock mass of the same blasting may be composed of a plurality of geologic bodies with large differences, the distribution nonuniformity of explosives in the blasting pile and the open-air deep hole blasting action mechanism are considered, and the surface block size distribution of the blasting pile has regional differences.
Therefore, a three-dimensional simulation method for blasting and blasting the surface rock of the pile by using the open deep hole steps is required.
Disclosure of Invention
The embodiment of the application provides a three-dimensional simulation method for surface rock masses of open-air deep hole bench blasting, which comprises the following steps:
step S1, determining the explosive rock block form and the explosive three-dimensional form;
step S2, dividing local areas of the surface of the blasting pile and counting the block size:
dividing the surface of the blasting pile into areas, and establishing a corresponding block size distribution curve for each area; carrying out region division and region block degree statistics according to the following method and steps to form a local block degree distribution curve;
step 1: firstly, performing area division on 6 surfaces of the trapezoidal explosive stack, and further dividing local areas according to the positions and the number of large blocks on the top surface and the slope surface of the explosive stack;
step 2: the local area dividing method comprises the following steps:
a1, carrying out image processing on the acquired image to obtain the size of the rock and the position information of the center of the rock in the image;
a2, dividing local area, searching for large rock (the large rock is unqualified rock, i.e. the size of the rock exceeds the design size (L)max) The size is determined according to actual engineering), such as when the block size is greater than LmaxIf the block size is large, the search block size is larger than LmaxAnd obtaining the position of the rock; calculating the distance between the large rock blocks, and presetting the local area range: the center of the large rock mass is taken as the circle center, and the radius is r (r is more than or equal to L)max/2,) if:
case 1: when the distance between the large rock masses is larger than 2r, the center of each large rock mass is taken as the circle center, the radius is r, the large rock masses are divided into separate areas, and other areas form a separate area;
case 2: when the distance between the large rock mass and the large rock mass is less than 2r, if the large rock mass A and the large rock mass B respectively draw a circle by taking the rock mass A and the rock mass B as the circle center and r as the radius, and two tangent lines L of the two circles are used1And L2And a semicircle A of two circles1And B2Forming a new area for the boundary, if the distance between the rock A or the rock B and other large blocks is less than 2r, processing according to the method, combining the two areas to form a new area, and processing the rest according to the method of the case 1;
and step 3: local region block degree statistics;
respectively counting the block degrees of each local area on the top surface and the slope of the blasting pile, and drawing a local block degree distribution curve; counting all the bulk degrees on the top surface and the slope surface of the blasting pile to form an integral blasting pile surface bulk degree distribution curve; except for the top surface and the slope surface, the lump distribution of other surfaces of the blasting pile is assumed to be a whole blasting pile surface lump distribution curve;
step S3, three-dimensional simulation of the surface rock of the blasting pile:
because the whole explosive stack is simplified into a trapezoidal cylinder, 12 edges and 8 vertexes of the cylinder are respectively related to two surfaces and three surfaces, and therefore preferential treatment must be carried out; since each surface of the trapezoidal column is different from that of the trapezoidal column which is simplified by the whole exploded pile, 6 surfaces and 12 sides of the exploded pile are specified;
the overall algorithm of the three-dimensional simulation of the surface of the blasting pile comprises the following steps:
b1, respectively randomly generating 8 sphere recording radiuses according to the local block size distribution model of the blasting pile, and randomly endowing the recording radiuses to 12 sides, namely the sphere center is positioned at the end point of the side;
b2, processing the 12 sides according to the directions of respective arrows, so that the residual part of the edges after the end point spherical radius is removed is filled with the sphere, namely the spherical center is positioned on the side; when the diameter of a sphere is randomly generated, judging which area the sphere is located in, and then randomly generating the diameter of the sphere according to a shot local block size distribution model;
b3, according to the rule in fig. 3, extracting the relevant edges to form 6 faces, dividing the 6 faces into regions according to the result of the region division, counting the sphere diameter distribution on the 6 faces, and using the face 1 as the current face; note that in order to form a closed chain, the direction processing needs to be performed on the relevant side according to the information in the figure, so that the spheres on the edge on the whole surface are arranged clockwise or counterclockwise, and all the spheres on the surface are assigned with values in sequence on the basis, that is, each sphere has a corresponding serial number;
b4, randomly searching a point in the current surface as an end point, such as the center point of the surface, counting the diameter of the sphere on the current surface to form the size data of the existing sphere on the current surface, and counting the diameter of the sphere in each area to form the size data of the existing sphere in the area;
b5, calculating a ball which is farthest from the terminal point among the balls forming the initial closed chain to be used as a ball 1, searching two adjacent balls according to the serial number, and using the ball which is farthest from the terminal point among the two adjacent balls to be used as a ball 2;
b6, when the ball 1 and the ball 2 are in the same region, generating a new diameter of the ball according to the block size distribution curve of the region and the existing ball size data of the region; when the ball 1 and the ball 2 are positioned in different areas, generating the diameter of a new ball according to an area block degree distribution curve of a small ball and size data of an existing ball in the area;
b7, calculating the position of the new sphere according to the radii of the sphere 1, the sphere 2 and the existing new sphere, and paying attention to the fact that the sphere center of the sphere is located on the surface;
b8, judging the relation between the new sphere and the existing sphere on the current surface, if there is coincidence, then using the original sphere diameter to subtract a certain smaller random number to form a new diameter, and returning to step B7; until the new sphere is not coincident with the existing sphere on the current surface, go to step B9;
b9, a small gap is formed among the new sphere, the sphere 1 and the sphere 2, and a small sphere can be filled in the middle gap to be tangent with the three spheres;
b10, updating the closed chain according to the relation between the new sphere and the adjacent spheres of the sphere 1 and the sphere 2, and adding the diameter of the new sphere into the size data of the existing sphere on the current surface and the size data of the existing sphere in the area respectively; if a certain sphere is removed when the closed chain is updated, counting the removed spheres to form total sphere data, and returning to the step 5;
b11, stopping the previous filling process when the previous is filled, that is, the radius of the sphere generated in step 8 is smaller than a certain smaller random number and the unsuccessful times exceed the allowed trial times;
b12, judging whether all the surfaces 6 are filled, if not, setting the next surface as the current surface, and returning to the step 4; if so, go to step B13;
b13, rotating the 6 surfaces according to the corresponding positions to form a three-dimensional blasting bulk and drawing;
and S4, three-dimensional simulation results of the rock mass on the surface of the blasting pile.
The embodiment of the application adopts the following technical scheme: in the step S1, the blasting pile rock blocks are different and irregular in shape, and the rock blocks in the whole blasting pile are difficult to describe by one shape; in order to reduce the calculation amount of numerical simulation and shorten the calculation time, the method is more suitable, and therefore, a sphere is used for replacing the explosive rock when the explosive rock is subjected to three-dimensional simulation.
The embodiment of the application adopts the following technical scheme: the surface of the open-air deep hole step blasting pile is uneven, wherein the slope surface of the blasting pile is similar to an inclined surface, the top surface of the blasting pile is similar to a paraboloid, and the slope surface and the top surface have obvious limits.
The embodiment of the application adopts the following technical scheme: for simplifying model calculation, the blasting funnel section and the slope surface of the top surface of the blasting pile are simplified, and the blasting pile is integrally similar to a trapezoidal cylinder with a trapezoid serving as a blasting pile section.
The embodiment of the application adopts the following technical scheme: in step S2, in order to make the three-dimensional simulation of the surface of the detonation pile more accurately reduce the rock mass in the detonation pile, it is proposed to perform region division on the surface of the detonation pile and establish a corresponding blockiness distribution curve for each region.
The embodiment of the application adopts the following technical scheme: in step B1, the balls at the end points of the joining edges are the same.
The embodiment of the application adopts the following technical scheme: in step B5, in searching for ball 1 and ball 2, if a plurality of farthest spheres are present, one ball is randomly drawn as the farthest sphere among the plurality of farthest spheres.
The embodiment of the application adopts the following technical scheme: in step B9, the heap is broken up with many unidentified fine particles, which can be used as a supplement.
The embodiment of the application adopts the following technical scheme: step S4 includes the steps of: and programming the whole algorithm for three-dimensional simulation of the surface of the detonation pile, and realizing the three-dimensional simulation of the surface of the detonation pile based on the obtained main parameters of the form of the detonation pile and the bulk distribution of the detonation pile.
The embodiment of the application adopts at least one technical scheme which can achieve the following beneficial effects:
according to the method, the surface of the detonation pile is partitioned, and the rock mass of each area is subjected to statistical analysis to obtain a detonation pile block size grading accumulation curve and a large block rate of the area, so that a detonation pile surface block size local distribution model is formed; the model and the information of the explosive pile form are combined, an explosive pile surface rock block generation algorithm is formed by improving an inward arrangement scheme in a discrete unit construction algorithm, and the three-dimensional model simulation function of the explosive pile form and the surface rock block is realized by adopting a three-dimensional visualization technology.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the application and together with the description serve to explain the application and not to limit the application. In the drawings:
FIG. 1 is a graph of the segmentation result in the embodiment of the present application;
FIG. 2 is a schematic diagram illustrating a division of a local area on the surface of a burst stack in the embodiment of the present application;
FIG. 3 is a specified diagram of the burst stack in the embodiment of the present application;
FIG. 4 is a simple flow chart of a three-dimensional simulation algorithm of the surface of the blasting pile in the embodiment of the application;
FIG. 5 is a three-dimensional simulation effect diagram of the surface of the blasting pile in the embodiment of the application;
fig. 6 is a diagram of a filling effect of a three-dimensional simulated gap on a surface of a burst stack in the embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the technical solutions of the present application will be described in detail and completely with reference to the following specific embodiments of the present application and the accompanying drawings. It should be apparent that the described embodiments are only some of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
The technical solutions provided by the embodiments of the present application are described in detail below with reference to the accompanying drawings.
Examples
A three-dimensional simulation method for surface rock blocks of open-air deep hole bench blasting comprises the following steps:
1. determination of blasting pile rock block form and blasting pile three-dimensional form
The blasting pile rock blocks are different and irregular in shape, and it is difficult to describe the rock blocks in the whole blasting pile by using one shape. In order to reduce the calculation amount of numerical simulation and shorten the calculation time, the method is more suitable, and therefore, a sphere is used for replacing the explosive rock when the explosive rock is subjected to three-dimensional simulation. The surface of the open-air deep hole step blasting pile is uneven, wherein the slope surface of the blasting pile is similar to an inclined surface, the top surface of the blasting pile is similar to a paraboloid, and the slope surface and the top surface have obvious limits. For simplifying model calculation, the blasting funnel section and the slope surface of the top surface of the blasting pile are simplified, and the blasting pile is integrally similar to a trapezoidal cylinder with a trapezoid serving as a blasting pile section.
2. Explosive pile surface local area division and block size statistics
In order to enable the three-dimensional simulation of the surface of the blasting pile to more accurately reduce the rock blocks in the blasting pile, the surface of the blasting pile is divided into regions, and a corresponding blockiness distribution curve is established for each region. The method comprises the following steps of carrying out region division and region block degree statistics to form a local block degree distribution curve.
Step 1: firstly, 6 surfaces of the trapezoidal blasting pile are used for carrying out region division, and then local regions are further divided according to the positions and the number of the large blocks on the top surface and the slope surface of the blasting pile.
Step 2: a local area dividing method.
(1) The acquired image is subjected to image processing to obtain the size of the rock and the position information of the center of the rock in the image (figure 1).
(2) And dividing local areas. Searching for large rock pieces (large rock pieces refer to unqualified rock pieces, i.e. the size of the rock pieces exceeds the design specification size (L) according to the result of image processingmax) The size is determined according to actual engineering), such as when the block size is greater than LmaxIf the block size is large, the search block size is larger than LmaxAnd obtaining the location of the rock mass. Calculating the distance between the large rock blocks, and presetting the local area range: the center of the large rock mass is taken as the circle center, and the radius is r (r is more than or equal to L)max/2) if
Case 1: when the distance between the large rock masses is larger than 2r, the center of each large rock mass is taken as the center of a circle, the radius is r, the large rock masses are divided into separate areas, and other areas on the surface form a separate area, as shown in fig. 2 (a).
Case 2: when the distance between the large rock mass and the large rock mass is less than 2r, if the large rock mass A and the large rock mass B respectively draw a circle by taking the rock mass A and the rock mass B as the circle center and r as the radius, and two tangent lines L of the two circles are used1And L2And a semicircle A of two circles1And B2If the distance between the rock A or rock B and other large blocks is less than 2r, the new area is formed by combining the two areas according to the method. The rest is processed as in case 1. As shown in fig. 2 (b).
And step 3: and carrying out local region block degree statistics.
And respectively counting the block degrees of each local area on the top surface and the slope of the blasting pile to form a local block degree distribution curve.
And counting all the bulk degrees on the top surface and the slope surface of the blasting pile to form an integral blasting pile surface bulk degree distribution curve.
Except for the top surface and the slope surface, the lump distribution of other surfaces of the blasting pile is assumed to be the integral blasting pile surface lump distribution curve.
3. Three-dimensional simulation of blasting surface rock mass
Since the whole of the blasting pile is simplified into a trapezoidal column, 12 sides and 8 vertexes of the column are respectively related to two faces and three faces, and therefore, priority treatment must be carried out. Since each surface of the trapezoidal column is different from that of the trapezoidal column simplified from the whole exploded stack, 6 surfaces and 12 sides of the exploded stack are defined as shown in fig. 3;
as shown in FIG. 3, the front face is face 1, the left face is face 2, the right face is face 3, the rear face is face 4, the bottom face is face 5, and the top face is face 6. The arrow is the direction of the edge. The vertices indicated in the figure are (0, 0) points of the plane in a two-dimensional coordinate system.
The overall algorithm of the three-dimensional simulation of the surface of the blasting pile comprises the following steps:
(1) according to the local block distribution model of the burst stack, 8 sphere recording radiuses are respectively generated randomly and are assigned to 12 sides randomly, namely the sphere center is located at the end point of the side. Note that the balls at the ends of the abutting edges are the same;
(2) the 12 sides are processed according to the directions of respective arrows, so that the residual parts after the end point spherical radii are removed are filled by the spheres, namely the spherical centers are positioned on the sides. When the diameter of a sphere is randomly generated, judging which area the sphere is located in, and then randomly generating the diameter of the sphere according to a shot local block size distribution model;
(3) according to the rule in fig. 3, 6 planes are formed by extracting the relevant edges, and the 6 planes are divided into regions according to the result of the region division, and the sphere diameter distribution on the 6 planes is counted. And face 1 as the current face. Note that in order to form a closed chain, the direction processing needs to be performed on the relevant side according to the information in the figure, so that the spheres on the edge on the whole surface are arranged clockwise or counterclockwise, and all the spheres on the surface are assigned with values in sequence on the basis, that is, each sphere has a corresponding serial number;
(4) randomly searching one point in the current surface as an end point, such as the center point of the surface, counting the diameter of a sphere on the current surface to form the size data of the existing sphere on the current surface, and counting the diameter of the sphere in each area to form the size data of the existing sphere in the area;
(5) and calculating a ball which is farthest from the terminal point in the balls forming the initial closed chain to be used as a ball 1, searching two adjacent balls according to the serial numbers, and using the ball which is farthest from the terminal point in the two adjacent balls to be used as a ball 2. Note that, in searching for the balls 1 and 2, if a plurality of farthest spheres are present, one ball is randomly drawn as the farthest sphere among the plurality of farthest spheres.
(6) When the ball 1 and the ball 2 are positioned in the same region, generating the diameter of a new ball according to the block size distribution curve of the region and the existing ball size data of the region; when the ball 1 and the ball 2 are located in different areas, the diameter of a new ball is generated according to the area block size distribution curve of the small ball and the existing ball size data of the area.
(7) The position of the new sphere is calculated from the radii of sphere 1 and sphere 2 and the existing new sphere. Note that the sphere center lies on the face.
(8) And (4) judging the relation between the new sphere and the existing sphere on the current surface, if the new sphere is overlapped, subtracting a certain smaller random number from the diameter of the original sphere to form a new diameter, and returning to the step (7). And (4) performing the step (9) until the new sphere is not overlapped with the existing sphere on the current surface.
(9) There will be a small gap between the new ball, ball 1 and ball 2, and a small ball can be filled in the middle gap to be tangent to the three balls. (the blasting pile will have many unidentifiable fine particles, which can be supplemented).
(10) And updating the closed chain according to the relation between the new sphere and the adjacent spheres of the sphere 1 and the sphere 2, and adding the diameter of the new sphere into the size data of the existing sphere on the current surface and the size data of the existing sphere in the area respectively. If a sphere is rejected when updating the closed chain, the statistically rejected spheres form the total sphere data. And (5) returning.
(11) When the previous filling process is completed, i.e. the radius of the sphere generated in step (8) is smaller than a certain small random number, and the number of unsuccessful attempts exceeds the allowed number of attempts, the previous filling process is stopped.
(12) And (4) judging whether the 6 surfaces are all filled, if not, setting the next surface as the current surface, and returning to the step (4). If so, go to step (13).
(13) And rotating the 6 surfaces according to the corresponding positions to form a three-dimensional blasting stack body and drawing.
4. Three-dimensional simulation result of rock mass on surface of blasting pile
And programming the algorithm, and realizing three-dimensional simulation of the surface of the detonation pile based on the obtained main parameters of the detonation pile form and the obtained detonation pile block size distribution. FIG. 5 is a three-dimensional simulation effect diagram of the surface of the blasting pile at different angles.
In fig. 5, the visual angles of (a), (b) and (c) are outside the burst stack, and the visual angle of (d) is inside the burst stack. When the details are viewed enlarged for the entire burst, the situation as in fig. 5(d) occurs because most of the voids exist even if part of the voids are filled with step (9). Considering that the whole blasting pile is a regular trapezoid cylinder, the gap on the surface of the whole blasting pile is filled in a surface mode. The padding effect is shown in fig. 6.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
The above description is only an example of the present application and is not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.
Claims (9)
1. A three-dimensional simulation method for surface rock blocks of open-air deep hole bench blasting is characterized by comprising the following steps:
step S1, determining the explosive rock block form and the explosive three-dimensional form;
step S2, dividing local areas of the surface of the blasting pile and counting the block size:
dividing the surface of the blasting pile into areas, and establishing a corresponding block size distribution curve for each area; carrying out region division and region block degree statistics according to the following method and steps to form a local block degree distribution curve;
step 1: firstly, performing area division on 6 surfaces of the trapezoidal explosive stack, and further dividing local areas according to the positions and the number of large blocks on the top surface and the slope surface of the explosive stack;
step 2: the local area dividing method comprises the following steps:
a1, carrying out image processing on the acquired image to obtain the size of the rock and the position information of the center of the rock in the image;
a2, dividing local area, searching for large rock (the large rock is unqualified rock, i.e. the size of the rock exceeds the design size (L)max) The size is determined according to actual engineering), such as when the block size is greater than LmaxIf the block size is large, the search block size is larger than LmaxAnd obtaining the position of the rock; calculating the distance between the large rock blocks, and presetting the local area range: the center of the large rock mass is taken as the circle center, and the radius is r (r is more than or equal to L)max/2,) if:
case 1: when the distance between the large rock masses is larger than 2r, the center of each large rock mass is taken as the circle center, the radius is r, the large rock masses are divided into separate areas, and other areas form a separate area;
case 2: when the distance between the large rock mass and the large rock mass is less than 2r, if the large rock mass A and the large rock mass B respectively draw a circle by taking the rock mass A and the rock mass B as the circle center and r as the radius, and two tangent lines L of the two circles are used1And L2And a semicircle A of two circles1And B2Forming a new area for the boundary, if the distance between the rock A or the rock B and other large blocks is less than 2r, processing according to the method, combining the two areas to form a new area, and forming the rest according to the situationProcessing by the method of condition 1;
and step 3: local region block degree statistics;
respectively counting the block degrees of each local area on the top surface and the slope of the blasting pile, and drawing a local block degree distribution curve; counting all the bulk degrees on the top surface and the slope surface of the blasting pile to form an integral blasting pile surface bulk degree distribution curve; except for the top surface and the slope surface, the lump distribution of other surfaces of the blasting pile is assumed to be a whole blasting pile surface lump distribution curve;
step S3, three-dimensional simulation of the surface rock of the blasting pile:
because the whole explosive stack is simplified into a trapezoidal cylinder, 12 edges and 8 vertexes of the cylinder are respectively related to two surfaces and three surfaces, and therefore preferential treatment must be carried out; since each surface of the trapezoidal column is different from that of the trapezoidal column which is simplified by the whole exploded pile, 6 surfaces and 12 sides of the exploded pile are specified;
the overall algorithm of the three-dimensional simulation of the surface of the blasting pile comprises the following steps:
b1, respectively randomly generating 8 sphere recording radiuses according to the local block size distribution model of the blasting pile, and randomly endowing the recording radiuses to 12 sides, namely the sphere center is positioned at the end point of the side;
b2, processing the 12 sides according to the directions of respective arrows, so that the residual part of the edges after the end point spherical radius is removed is filled with the sphere, namely the spherical center is positioned on the side; when the diameter of a sphere is randomly generated, judging which area the sphere is located in, and then randomly generating the diameter of the sphere according to a shot local block size distribution model;
b3, according to the rule in fig. 3, extracting the relevant edges to form 6 faces, dividing the 6 faces into regions according to the result of the region division, counting the sphere diameter distribution on the 6 faces, and using the face 1 as the current face; note that in order to form a closed chain, the direction processing needs to be performed on the relevant side according to the information in the figure, so that the spheres on the edge on the whole surface are arranged clockwise or counterclockwise, and all the spheres on the surface are assigned with values in sequence on the basis, that is, each sphere has a corresponding serial number;
b4, randomly searching a point in the current surface as an end point, such as the center point of the surface, counting the diameter of the sphere on the current surface to form the size data of the existing sphere on the current surface, and counting the diameter of the sphere in each area to form the size data of the existing sphere in the area;
b5, calculating a ball which is farthest from the terminal point among the balls forming the initial closed chain to be used as a ball 1, searching two adjacent balls according to the serial number, and using the ball which is farthest from the terminal point among the two adjacent balls to be used as a ball 2;
b6, when the ball 1 and the ball 2 are in the same region, generating a new diameter of the ball according to the block size distribution curve of the region and the existing ball size data of the region; when the ball 1 and the ball 2 are positioned in different areas, generating the diameter of a new ball according to an area block degree distribution curve of a small ball and size data of an existing ball in the area;
b7, calculating the position of the new sphere according to the radii of the sphere 1, the sphere 2 and the existing new sphere, and paying attention to the fact that the sphere center of the sphere is located on the surface;
b8, judging the relation between the new sphere and the existing sphere on the current surface, if there is coincidence, then using the original sphere diameter to subtract a certain smaller random number to form a new diameter, and returning to step B7; until the new sphere is not coincident with the existing sphere on the current surface, go to step B9;
b9, a small gap is formed among the new sphere, the sphere 1 and the sphere 2, and a small sphere can be filled in the middle gap to be tangent with the three spheres;
b10, updating the closed chain according to the relation between the new sphere and the adjacent spheres of the sphere 1 and the sphere 2, and adding the diameter of the new sphere into the size data of the existing sphere on the current surface and the size data of the existing sphere in the area respectively; if a certain sphere is removed when the closed chain is updated, counting the removed spheres to form total sphere data, and returning to the step 5;
b11, stopping the previous filling process when the previous is filled, that is, the radius of the sphere generated in step 8 is smaller than a certain smaller random number and the unsuccessful times exceed the allowed trial times;
b12, judging whether all the surfaces 6 are filled, if not, setting the next surface as the current surface, and returning to the step 4; if so, go to step B13;
b13, rotating the 6 surfaces according to the corresponding positions to form a three-dimensional blasting bulk and drawing;
and S4, three-dimensional simulation results of the rock mass on the surface of the blasting pile.
2. The three-dimensional simulation method for the surface rock of the open-air deep-hole bench blasting detonating pile according to claim 1, characterized in that in step S1, the shape of the detonating pile rock is different and irregular, and it is difficult to describe the rock in the whole detonating pile by one shape; in order to reduce the calculation amount of numerical simulation and shorten the calculation time, the method is more suitable, and therefore, a sphere is used for replacing the explosive rock when the explosive rock is subjected to three-dimensional simulation.
3. The three-dimensional simulation method for the surface rock mass of the open-air deep-hole bench blasting detonating pile according to claim 2, characterized in that the surface of the open-air deep-hole bench blasting pile is uneven, wherein the slope surface of the pile is similar to an inclined surface, the top surface of the pile is similar to a paraboloid, and the slope surface and the top surface have obvious limits.
4. The three-dimensional simulation method for the surface rock mass of the open-air deep-hole bench blasting detonating pile according to claim 3, characterized in that for the simplification of model calculation, the blasting funnel section and the slope surface of the top surface of the detonating pile are simplified, and the whole detonating pile is similar to a trapezoidal cylinder with a trapezoidal section as the detonating pile.
5. The three-dimensional simulation method for surface rock blocks of the open-air deep-hole bench blasting detonating pile according to claim 1, wherein in step S2, in order to make the three-dimensional simulation of the surface of the detonating pile more accurately restore the rock blocks in the detonating pile, it is proposed to perform area division on the surface of the detonating pile and establish a corresponding block size distribution curve for each area.
6. The three-dimensional simulation method for the surface rock block of the open-air deep-hole bench blasting explosive according to claim 1, wherein in step B1, the balls at the end points of the connecting edges are the same.
7. The three-dimensional simulation method for the surface rock block of the open-air deep-hole bench blasting explosive pile according to claim 1, wherein in the step B5, if a plurality of farthest spheres occur when searching for the sphere 1 and the sphere 2, one sphere is randomly selected from the farthest spheres as the farthest sphere.
8. The method for three-dimensional simulation of surface rock of open-air deep-hole bench blasting explosive pile according to claim 1, wherein in step B9, the explosive pile contains many unidentifiable fine particles which can be used as a supplement.
9. The three-dimensional simulation method for the surface rock block of the open-air deep-hole bench blasting explosive according to claim 1, wherein the step S4 comprises the steps of: and programming the whole algorithm for three-dimensional simulation of the surface of the detonation pile, and realizing the three-dimensional simulation of the surface of the detonation pile based on the obtained main parameters of the form of the detonation pile and the bulk distribution of the detonation pile.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110266412.XA CN113674406B (en) | 2021-03-11 | 2021-03-11 | Three-dimensional simulation method for surface rock mass of open-air deep hole bench blasting pile |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110266412.XA CN113674406B (en) | 2021-03-11 | 2021-03-11 | Three-dimensional simulation method for surface rock mass of open-air deep hole bench blasting pile |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113674406A true CN113674406A (en) | 2021-11-19 |
| CN113674406B CN113674406B (en) | 2022-06-28 |
Family
ID=78538035
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110266412.XA Active CN113674406B (en) | 2021-03-11 | 2021-03-11 | Three-dimensional simulation method for surface rock mass of open-air deep hole bench blasting pile |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113674406B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115479512A (en) * | 2022-09-20 | 2022-12-16 | 中煤科工集团沈阳研究院有限公司 | Intelligent drilling and blasting pile block rate optimization method for strip mine |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108596881A (en) * | 2018-04-09 | 2018-09-28 | 中国科学院武汉岩土力学研究所 | The intelligent image statistical method of rock grain size after a kind of explosion |
| CN109030305A (en) * | 2018-06-08 | 2018-12-18 | 武汉大学 | A kind of gradating material specific surface area based on 3 D laser scanning and image processing techniques determines method |
| CN110068573A (en) * | 2019-04-22 | 2019-07-30 | 武汉大学 | Based on photogrammetric and particle contact theory blast fragmentation size measurement method |
| US20190293400A1 (en) * | 2018-03-21 | 2019-09-26 | Accenture Global Solutions Limited | System and Method for Controlling a Drill and Blast Event |
| CN111199108A (en) * | 2020-01-07 | 2020-05-26 | 包钢集团矿山研究院（有限责任公司） | Blasting blockiness prediction method of gridding segmentation numerical simulation blasting pile model |
| CN111383174A (en) * | 2020-02-27 | 2020-07-07 | 广东锡源爆破科技股份有限公司 | Pile bursting data acquisition method for photogrammetry |
| CN112285107A (en) * | 2020-10-19 | 2021-01-29 | 中国科学院武汉岩土力学研究所 | Method for realizing three-dimensional blasting block size distribution measurement through plane photographing |
-
2021
- 2021-03-11 CN CN202110266412.XA patent/CN113674406B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190293400A1 (en) * | 2018-03-21 | 2019-09-26 | Accenture Global Solutions Limited | System and Method for Controlling a Drill and Blast Event |
| CN108596881A (en) * | 2018-04-09 | 2018-09-28 | 中国科学院武汉岩土力学研究所 | The intelligent image statistical method of rock grain size after a kind of explosion |
| CN109030305A (en) * | 2018-06-08 | 2018-12-18 | 武汉大学 | A kind of gradating material specific surface area based on 3 D laser scanning and image processing techniques determines method |
| CN110068573A (en) * | 2019-04-22 | 2019-07-30 | 武汉大学 | Based on photogrammetric and particle contact theory blast fragmentation size measurement method |
| CN111199108A (en) * | 2020-01-07 | 2020-05-26 | 包钢集团矿山研究院（有限责任公司） | Blasting blockiness prediction method of gridding segmentation numerical simulation blasting pile model |
| CN111383174A (en) * | 2020-02-27 | 2020-07-07 | 广东锡源爆破科技股份有限公司 | Pile bursting data acquisition method for photogrammetry |
| CN112285107A (en) * | 2020-10-19 | 2021-01-29 | 中国科学院武汉岩土力学研究所 | Method for realizing three-dimensional blasting block size distribution measurement through plane photographing |
Non-Patent Citations (2)
| Title |
|---|
| XIAOHUA DING ET AL: "Distribution Characteristics of Fragments Size and Optimization of Blasting Parameters Under Blasting Impact Load in Open-Pit Mine", 《IEEE ACCESS》 * |
| 李辉: "基于三维激光点云的爆堆块度分析", 《中国优秀硕士学位论文全文数据库电子期刊 工程科技I辑》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115479512A (en) * | 2022-09-20 | 2022-12-16 | 中煤科工集团沈阳研究院有限公司 | Intelligent drilling and blasting pile block rate optimization method for strip mine |
| CN115479512B (en) * | 2022-09-20 | 2023-08-22 | 中煤科工集团沈阳研究院有限公司 | Intelligent drilling and blasting pile large block rate optimization method for strip mine |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113674406B (en) | 2022-06-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110348562A (en) | The quantization strategy of neural network determines method, image-recognizing method and device | |
| CN113674406B (en) | Three-dimensional simulation method for surface rock mass of open-air deep hole bench blasting pile | |
| CN106934826B (en) | Rock slope structure refined modeling and block identification method | |
| Suchde et al. | Point cloud generation for meshfree methods: An overview | |
| CN109241646B (en) | Multi-factor two-dimensional soil-rock mixture generation method based on elliptical stacking and random field | |
| Noroozi et al. | 3D key‐group method for slope stability analysis | |
| CN112288874B (en) | Gamma radiation modeling calculation simulation method based on CAD model and Boolean operation | |
| CN111159927A (en) | Numerical modeling method for three-dimensional irregular-shaped particle throwing based on voxel matrix | |
| CN111104728B (en) | Simulation prediction method and device for structural surface wear | |
| Ren et al. | New generation and separation method for the computational modelling of the mesostructure of concrete | |
| Manak | Voronoi‐based detection of pockets in proteins defined by large and small probes | |
| CN112991539B (en) | Three-dimensional simulation method for explosive pile and block size distribution based on discrete elements | |
| CN105786512A (en) | Program generation method and dimension manager | |
| CN116228766B (en) | Intelligent regulation and control method and system for plasma processing equipment | |
| CN113284241B (en) | Three-dimensional simulation method for surface rock mass of blasting pile with pile shape distributed in Weibull mode | |
| US20170017741A1 (en) | Computational modelling for engineering analysis | |
| US11360744B2 (en) | Two-dimensional data matching method, device and logic circuit | |
| Tsang et al. | Enhanced affine invariant matching of broken boundaries based on particle swarm optimization and the dynamic migrant principle | |
| CN115147631A (en) | Blasting block degree rapid identification method and system based on three-dimensional point cloud data | |
| Sampl | Medial axis construction in three dimensions and its application to mesh generation | |
| CN112991538B (en) | Three-dimensional simulation method for explosive stack body based on Weibull distribution of discrete element pair stack shape | |
| CN114882496A (en) | Three-dimensional component similarity calculation method based on depth image | |
| Lutton et al. | A genetic algorithm with sharing for the detection of 2D geometric primitives in images | |
| CN112214927A (en) | Machine learning-based ultrahigh-speed collision fragment cloud rapid simulation method | |
| CN116958607B (en) | Data processing method and device for target damage prediction |
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 |