CN112991539A - Three-dimensional simulation method for explosive pile and block size distribution based on discrete elements - Google Patents
Three-dimensional simulation method for explosive pile and block size distribution based on discrete elements Download PDFInfo
- Publication number
- CN112991539A CN112991539A CN202110265853.8A CN202110265853A CN112991539A CN 112991539 A CN112991539 A CN 112991539A CN 202110265853 A CN202110265853 A CN 202110265853A CN 112991539 A CN112991539 A CN 112991539A
- Authority
- CN
- China
- Prior art keywords
- sphere
- pile
- blasting
- spheres
- new
- 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/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application discloses three-dimensional simulation method of blasting pile and block size distribution based on discrete elements, relates to the technical field of blasting effect evaluation, and comprises the following steps: the three-dimensional explosive pile simulation algorithm is divided into two steps, wherein the first step is to generate rock blocks on the surface of the whole explosive pile; secondly, the surface of the whole blasting pile is used as a closed initial surface, and the inside of the whole blasting pile is filled on the basis of the closed initial surface; the first step comprises step 1 to step 12, and the second step comprises step 13 to step 20. The method and the device have the advantages that the algorithm of the three-dimensional simulation party of the explosive pile and the block size distribution is programmed, and the three-dimensional simulation of the explosive pile is realized based on the obtained main parameters of the explosive pile form and the explosive pile block size distribution; the blasting rock is replaced by a ball; simplifying the explosive pile into a trapezoidal cylinder with a trapezoidal section as the explosive pile, and providing an explosive pile shape measuring method; the three-dimensional simulation method for blasting pile and block size distribution of open-air deep hole bench blasting is firstly provided.
Description
Technical Field
The application relates to the technical field of blasting effect evaluation, in particular to a three-dimensional simulation method for blasting pile and block size distribution based on discrete elements.
Background
The research on the explosive pile formed by open-air step deep hole blasting and the block size distribution rule thereof is an important means for evaluating the blasting effect, and the intuitive expression of the explosive pile block size distribution rule by using a three-dimensional simulation method is an important development direction of open-air rock-soil blasting. Because the blasting pile formed by open-air step deep hole blasting and the block size distribution thereof have the characteristics of large scale, complex pile shape, irregular rock block shape, overlapping adhesion, large discreteness and the like, the difficulty in establishing a three-dimensional model of the blasting pile and the block size distribution thereof is high, and the research results are few. The patent relates to the relevant technology: (1) an equivalent simulation method of the form of the blasting rock mass; (2) constructing a geometric model of a three-dimensional form of the explosive pile; (3) and (3) three-dimensional simulation of the blasting rock mass. The technology (1) mainly researches the shape and size equivalence of the blasting pile rock mass, and enables irregular rock to be equivalent to regular rock mass; the technology (2) mainly establishes a model meeting the requirements of the technology (3); the technology (3) mainly adopts large-scale discrete element commercial software developed abroad, such as 3DEC discrete element software for simulating the dynamic crushing and throwing process of step blasting. And (5) performing equivalent simulation on the three-dimensional virtual rock mass of the blasting pile.
The main problems of the achievement used for establishing the three-dimensional model of the explosive pile and the lump distribution are as follows: whether the equivalent rock researched by the technology (1) is suitable for constructing a blasting pile three-dimensional simulation rock unit or not is not carried out, and related research work is not carried out; the main commercial software developed by the technology (3) and the related technology are protected by foreign patents, and have no independent intellectual property rights at home.
Therefore, the three-dimensional model of the explosive pile and the lump size distribution thereof is constructed by means of a discrete body theory and a discrete element method.
Disclosure of Invention
The embodiment of the application provides a three-dimensional simulation method for blasting and block size distribution based on discrete elements, which comprises the following steps:
the three-dimensional explosive pile simulation algorithm is divided into two steps, wherein the first step is to generate rock blocks on the surface of the whole explosive pile; secondly, the surface of the whole blasting pile is used as a closed initial surface, and the inside of the whole blasting pile is filled on the basis of the closed initial surface;
the first step is as follows:
step 1, respectively randomly generating 8 sphere recording radiuses, and randomly assigning the radiuses to 12 sides, namely, the sphere center is positioned at the end point of the side;
step 2, the 12 sides are respectively processed according to the directions of respective arrows, so that the residual part of the 12 sides after the end point sphere radius is removed is filled by a sphere, namely the sphere center is positioned on the side;
step 3, respectively extracting relevant edges to form 6 faces according to the provisions in the figure 1, counting the diameter distribution of spheres on the 6 faces, and taking 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;
step 4, randomly searching a point in the current surface as an end point, such as a center point of the surface, and counting the diameter of the sphere on the current surface to form the size data of the existing sphere;
step 5, calculating a ball which is farthest from the terminal point in the balls forming the initial closed chain to serve as a ball 1, searching two adjacent balls according to the serial number, and taking the ball which is farthest from the terminal point in the two adjacent balls to serve as a ball 2; when searching for the ball 1 and the ball 2, if a plurality of farthest balls appear, randomly drawing one ball as a farthest ball from the plurality of farthest balls;
step 6, generating the diameter of a new sphere according to the block size distribution curve and the size data of the existing sphere;
step 7, 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;
step 8, 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; until the new sphere is not coincident with the existing sphere on the current surface, go to step 9;
step 9, 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 spheres; 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;
step 10, when the current surface is filled, namely the radius of the sphere generated in the step 8 is smaller than a certain smaller random number, and the unsuccessful times exceed the allowed trial times, stopping the current surface filling process;
step 11, 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 yes, go to step 12;
step 12, rotating the 6 surfaces according to corresponding positions to form a three-dimensional pile blasting surface;
the second step is that:
step 13, calculating the center coordinate of the detonation as a central point according to the trapezoid cylinder similar to the detonation, and counting the diameters of all spheres on the surface of the detonation to form existing sphere size data;
step 14, calculating the sphere farthest from the central point among the spheres forming the initial closed surface as a sphere 1, taking the sphere farthest from the central point among the spheres adjacent to the sphere 1 as a sphere 2, and taking the sphere farthest from the central point among the spheres adjacent to both the sphere 1 and the sphere 2 as a sphere 3;
step 15, according to the block size distribution curve and the position of the existing sphere, paying attention to the fact that the sphere center of the sphere is located inside the blasting pile;
step 16, calculating the position of the new sphere according to the radii of the sphere 1, the sphere 2, the sphere 3 and the existing new sphere; note that the sphere center of the sphere is located inside the blasting pile;
step 17, judging the relationship between the new sphere and the existing sphere on the blasting pile, 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 16; until the new sphere does not coincide with the existing sphere on the current surface, go to step 18;
step 18, small gaps are formed among the new sphere, the sphere 1, the sphere 2 and the sphere 3, and a small sphere can be filled in the middle gap to be tangent to the four spheres;
step 19, updating the closed surface according to the relation between the new sphere and the adjacent spheres of the sphere 1, the sphere 2 and the sphere 3, and adding the diameter of the new sphere into the size data of the existing spheres; if a certain sphere is removed when the closed surface is updated, counting the removed spheres to form total sphere data, and returning to the step 15;
and 20, when the inside of the burst stack is filled, namely the radius of the sphere generated in the step 17 is smaller than a certain smaller random number, and the unsuccessful times exceed the allowed trial times, stopping the filling process inside the burst stack.
The embodiment of the application adopts the following technical scheme: the method also comprises the following steps before executing the pile-blasting three-dimensional simulation algorithm: and (4) acquiring main parameters of the blasting pile form.
The embodiment of the application adopts the following technical scheme: the method also comprises the following steps before the acquisition of the main parameters of the burst mode: and (5) constructing a three-dimensional explosive pile form model.
The embodiment of the application adopts the following technical scheme: the method also comprises the following steps before the construction of the three-dimensional shape model of the blasting pile: equivalent virtual of the shape of the mound rock.
The embodiment of the application adopts the following technical scheme: the method also comprises the following steps before the equivalent virtual of the shape of the explosive pile rock block: and the shape of the blasting pile rock is virtual.
The embodiment of the application adopts the following technical scheme: the detonative rock block morphology virtualization includes equivalent virtualization of rock block size and shape.
The embodiment of the application adopts the following technical scheme: in step 1, note that the balls at the end points of the abutting edges are the same.
The embodiment of the application adopts the following technical scheme: in step 12, the three-dimensional surface of the blasting pile is the closed surface in the second step.
The embodiment of the application adopts the following technical scheme: in step 18, the heap is broken up with many unidentifiable fine particles, which can be supplemented.
The embodiment of the application adopts at least one technical scheme which can achieve the following beneficial effects:
(1) replacing blasting rock blocks with spheres;
(2) simplifying the explosive pile into a trapezoidal cylinder with a trapezoidal section as the explosive pile, and providing an explosive pile shape measuring method;
(3) the three-dimensional simulation method for blasting pile and block size distribution of open-air deep hole bench blasting is firstly provided.
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 schematic diagram illustrating explosive pile specification in a three-dimensional simulation method based on discrete element explosive pile and block size distribution according to the present invention;
FIG. 2 is a simple flow chart of a three-dimensional simulation algorithm of the blasting pile in the three-dimensional simulation method based on the blasting pile and the block size distribution of the discrete elements
FIG. 3 is a three-dimensional simulation effect diagram of the blasting pile at different angles in the embodiment of the invention.
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
1. Blasting rock block shape virtual
The blasting pile rock block shape virtual comprises equivalent virtual of the size and the shape of the rock block:
(1) a blasting rock block shape;
the blasting pile rock blocks are different and irregular in shape and are influenced by the structural surface and rock properties of the rock body. Mainly comprises a sphere, an ellipsoid, a cuboid, a needle, a sheet and the like.
(2) Equivalent virtualization of the shape of the blasting rock mass;
in order to realize statistics and simulation of the explosive rock, irregular rock is equivalently virtualized into regular rock, the explosive image is identified by adopting an image identification technology, after the pixel area, the maximum chord length and the perimeter of each part (rock) are counted, the rock is generally virtualized into a circle, an ellipse, a square, a polygon and the like according to the actual shape of the rock, and the explosive rock is virtualized according to the principle that the area and the perimeter are not changed and further converted into a three-dimensional sphere, an ellipsoid, a polyhedron and the like.
Performing equivalent simulation on the three-dimensional virtual rock block of the blasting pile:
the three-dimensional virtual process of blasting is the stacking process of rock particles, and the calculation workload of the particle contact algorithm is changed in a series manner along with the change of the particle shape. In actual engineering, the length of a blasting pile of open-air step deep hole blasting reaches dozens of meters to hundreds of meters, in order to realize three-dimensional simulation of the blasting pile, the simplest particle form of a particle contact algorithm needs to be selected for simulation, and in consideration of engineering applicability, a sphere is used for replacing a blasting pile rock mass to carry out three-dimensional simulation on the surface of the blasting pile.
2. Blasting a three-dimensional shape model;
the blasting pile surface of open-air deep hole bench blasting comprises four parts: slope, top surface and the side at both ends. The slope surface of the blasting pile is similar to an inclined surface, and the top surface of the blasting pile is similar to a paraboloid due to the existence of the blasting funnel section. Since the sides of the pile are usually constrained, approximating a plane perpendicular to the ground. The blasting funnel section of the top surface of the blasting pile is simplified for the simplification of model calculation by considering four blasting pile surfaces, and the whole blasting pile is similar to a trapezoidal pile body with a trapezoidal shape as the section of the blasting pile. As shown in fig. 1.
The simplified exploding pile trapezoidal column body is composed of 12 sides and 8 vertexes, each side is related to two faces, and each vertex is related to three faces, so that simulation needs to be carried out according to the sequence of points, sides, faces and bodies. According to the three-dimensional shape model of the explosive pile, 6 surfaces and 12 sides of the explosive pile are defined, and the definition is as shown in figure 1:
as shown in FIG. 1, 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 vertex indicated in the figure is the (0,0) point of the plane in the two-dimensional coordinate system, i.e., the start point of the plane.
3. Acquiring main parameters of the explosive pile form;
the whole blasting pile is simplified into a trapezoidal column with a trapezoidal blasting pile section, as shown in figure 1. Only information about the section of the burst and the length of the burst, i.e. side 1, side 2, side 7 and side 10 in fig. 1, has to be obtained. The length of the side 1 is obtained by measuring the two sides of the detonation pile for multiple times by using RTK and then averaging. The slope bottom (side 1) of the actual pile-blasting is a curve, so that the lengths of the sides 2 are different, the measurement needs to be carried out at different positions of the slope bottom, the distance between a to-be-blasted area and a blasted area boundary line is calculated, and then the average value of the distances is calculated to be used as the length of the side 2. The distance from the boundary line of the top of the slope and the slope to the boundary line is calculated by the same method to be the edge 10. The top of the actual detonation slope is similar to a parabola, so in order to reduce the error of the detonation square quantity, the top of the detonation slope is measured for multiple times, and the average height difference from the bottom of the detonation slope is obtained as a side 7. Thereby forming the entire detonation surface.
4. A pile blasting three-dimensional simulation algorithm;
the three-dimensional pile-blasting simulation algorithm is divided into two steps. The first step is to generate rock blocks on the surface of the whole blasting pile; the second step is to fill the whole detonation pile interior based on the whole detonation pile surface as a closed initial surface.
The overall algorithm of the three-dimensional pile-blasting simulation comprises the following steps:
the first step is as follows:
and step 1, respectively randomly generating 8 sphere recording radiuses, and randomly assigning the recording radiuses to 12 sides, namely, the sphere center is positioned at the end point of the side. Note that the balls at the ends of the abutting edges are the same;
step 2, the 12 sides are respectively processed according to the directions of respective arrows, so that the residual part of the 12 sides after the end point sphere radius is removed is filled by a sphere, namely the sphere center is positioned on the side;
and 3, respectively extracting relevant edges to form 6 faces according to the specification in the figure 1, and counting the diameter distribution of the spheres on the 6 faces. 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;
step 4, randomly searching a point in the current surface as an end point, such as a center point of the surface, and counting the diameter of the sphere on the current surface to form the size data of the existing sphere;
and 5, calculating a ball which is farthest from the terminal point in the balls forming the initial closed chain to serve as a ball 1, searching two adjacent balls according to the serial numbers, and taking the ball which is farthest from the terminal point in the two adjacent balls to serve as a ball 2. When searching for the ball 1 and the 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.
Step 6, generating the diameter of a new sphere according to the block size distribution curve and the size data of the existing sphere;
and 7, calculating the position of the new sphere according to the radii of the sphere 1, the sphere 2 and the existing new sphere. Note that the sphere center lies on the face.
And 8, 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. Until the new sphere does not coincide with the existing sphere on the current surface, proceed to step 9.
And 9, 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 spheres. If a sphere is rejected when updating the closed chain, the statistically rejected spheres form the total sphere data. And returning to the step 5.
And step 10, stopping the current filling process when the current is filled, namely the radius of the sphere generated in the step 8 is smaller than a certain smaller random number and the unsuccessful times exceed the allowed trial times.
And 11, 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 12
And 12, rotating the 6 surfaces according to corresponding positions to form a three-dimensional pile blasting surface, namely a closed surface in the second step.
The second step is that:
step 13, calculating the center coordinate of the detonation as a central point according to the trapezoid cylinder similar to the detonation, and counting the diameters of all spheres on the surface of the detonation to form existing sphere size data;
step 14, calculating the sphere farthest from the central point among the spheres forming the initial closed surface as a sphere 1, taking the sphere farthest from the central point among the spheres adjacent to the sphere 1 as a sphere 2, and taking the sphere farthest from the central point among the spheres adjacent to both the sphere 1 and the sphere 2 as a sphere 3;
step 15, generating the diameter of a new sphere according to the block size distribution curve and the size data of the existing sphere;
and step 16, calculating the position of the new sphere according to the radii of the sphere 1, the sphere 2, the sphere 3 and the existing new sphere. Note that the sphere center is located inside the burst.
And step 17, judging the relationship between the new sphere and the existing sphere on the blasting pile, if the new sphere is overlapped with the existing sphere on the blasting pile, subtracting a certain smaller random number from the diameter of the original sphere to form a new diameter, and returning to the step 16. Until the new sphere does not coincide with the sphere already on the current surface, step 18 is performed.
In step 18, a small gap is formed among the new sphere, the sphere 1, the sphere 2 and the sphere 3, and a small sphere can be filled in the middle gap to be tangent to the four spheres. The blasting pile has many unidentifiable fine particles, which can be used as a supplement
And step 19, updating the closed surface according to the relation between the new sphere and the adjacent spheres of the sphere 1, the sphere 2 and the sphere 3, and adding the diameter of the new sphere into the size data of the existing spheres. If a sphere is rejected when the closed surface is updated, the statistically rejected spheres form total sphere data. Returning to step 14.
And 20, when the inside of the burst stack is filled, namely the radius of the sphere generated in the step 17 is smaller than a certain smaller random number, and the unsuccessful times exceed the allowed trial times, stopping the filling process inside the burst stack.
In summary, the following steps:
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 blasting and block size distribution based on discrete elements is characterized by comprising the following steps:
the three-dimensional explosive pile simulation algorithm is divided into two steps, wherein the first step is to generate rock blocks on the surface of the whole explosive pile; secondly, the surface of the whole blasting pile is used as a closed initial surface, and the inside of the whole blasting pile is filled on the basis of the closed initial surface;
the first step is as follows:
step 1, respectively randomly generating 8 sphere recording radiuses, and randomly assigning the radiuses to 12 sides, namely, the sphere center is positioned at the end point of the side;
step 2, the 12 sides are respectively processed according to the directions of respective arrows, so that the residual part of the 12 sides after the end point sphere radius is removed is filled by a sphere, namely the sphere center is positioned on the side;
step 3, respectively extracting relevant edges to form 6 faces according to the provisions in the figure 1, counting the diameter distribution of spheres on the 6 faces, and taking 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;
step 4, randomly searching a point in the current surface as an end point, such as a center point of the surface, and counting the diameter of the sphere on the current surface to form the size data of the existing sphere;
step 5, calculating a ball which is farthest from the terminal point in the balls forming the initial closed chain to serve as a ball 1, searching two adjacent balls according to the serial number, and taking the ball which is farthest from the terminal point in the two adjacent balls to serve as a ball 2; when searching for the ball 1 and the ball 2, if a plurality of farthest balls appear, randomly drawing one ball as a farthest ball from the plurality of farthest balls;
step 6, generating the diameter of a new sphere according to the block size distribution curve and the size data of the existing sphere;
step 7, 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;
step 8, 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; until the new sphere is not coincident with the existing sphere on the current surface, go to step 9;
step 9, 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 spheres; 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;
step 10, when the current surface is filled, namely the radius of the sphere generated in the step 8 is smaller than a certain smaller random number, and the unsuccessful times exceed the allowed trial times, stopping the current surface filling process;
step 11, 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 yes, go to step 12;
step 12, rotating the 6 surfaces according to corresponding positions to form a three-dimensional pile blasting surface;
the second step is that:
step 13, calculating the center coordinate of the detonation as a central point according to the trapezoid cylinder similar to the detonation, and counting the diameters of all spheres on the surface of the detonation to form existing sphere size data;
step 14, calculating the sphere farthest from the central point among the spheres forming the initial closed surface as a sphere 1, taking the sphere farthest from the central point among the spheres adjacent to the sphere 1 as a sphere 2, and taking the sphere farthest from the central point among the spheres adjacent to both the sphere 1 and the sphere 2 as a sphere 3;
step 15, according to the block size distribution curve and the position of the existing sphere, paying attention to the fact that the sphere center of the sphere is located inside the blasting pile;
step 16, calculating the position of the new sphere according to the radii of the sphere 1, the sphere 2, the sphere 3 and the existing new sphere; note that the sphere center of the sphere is located inside the blasting pile;
step 17, judging the relationship between the new sphere and the existing sphere on the blasting pile, 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 16; until the new sphere does not coincide with the existing sphere on the current surface, go to step 18;
step 18, small gaps are formed among the new sphere, the sphere 1, the sphere 2 and the sphere 3, and a small sphere can be filled in the middle gap to be tangent to the four spheres;
step 19, updating the closed surface according to the relation between the new sphere and the adjacent spheres of the sphere 1, the sphere 2 and the sphere 3, and adding the diameter of the new sphere into the size data of the existing spheres; if a certain sphere is removed when the closed surface is updated, counting the removed spheres to form total sphere data, and returning to the step 15;
and 20, when the inside of the burst stack is filled, namely the radius of the sphere generated in the step 17 is smaller than a certain smaller random number, and the unsuccessful times exceed the allowed trial times, stopping the filling process inside the burst stack.
2. The three-dimensional simulation method for blasting piles and bulk distribution based on discrete elements according to claim 1, which is characterized by further comprising the following steps before executing a blasting pile three-dimensional simulation algorithm: and (4) acquiring main parameters of the blasting pile form.
3. The three-dimensional simulation method for blasting pile and block size distribution based on discrete elements according to claim 2, characterized by further comprising the following steps before acquiring the main parameters of blasting pile form: and (5) constructing a three-dimensional explosive pile form model.
4. The three-dimensional simulation method for blasting piles and bulk distribution based on discrete elements according to claim 3, which is characterized by further comprising the following steps before the construction of a blasting pile three-dimensional shape model: equivalent virtual of the shape of the mound rock.
5. The three-dimensional simulation method of blasting pile and block size distribution based on discrete elements as claimed in claim 4, wherein before the equivalent virtual of the shape of the blasting pile rock mass, the method further comprises the following steps: and the shape of the blasting pile rock is virtual.
6. The three-dimensional simulation method of discrete element-based blasting pile and block size distribution according to claim 5, wherein the blasting pile rock block shape virtual comprises equivalent virtual of rock block size and shape.
7. The three-dimensional simulation method for blasting and bulk distribution based on discrete elements according to claim 1, wherein in step 1, the spheres at the end points of the adjacent edges are the same.
8. The three-dimensional simulation method for blasting and bulk distribution based on discrete elements according to claim 1, wherein in step 12, the surface of the three-dimensional blasting is the closed surface in the second step.
9. The method of claim 1, wherein in step 18, the blasting pile has a plurality of unidentifiable fine particles, which can be used as a supplement.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110265853.8A CN112991539B (en) | 2021-03-11 | 2021-03-11 | Three-dimensional simulation method for explosive pile and block size distribution based on discrete elements |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110265853.8A CN112991539B (en) | 2021-03-11 | 2021-03-11 | Three-dimensional simulation method for explosive pile and block size distribution based on discrete elements |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112991539A true CN112991539A (en) | 2021-06-18 |
| CN112991539B CN112991539B (en) | 2022-07-08 |
Family
ID=76335141
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110265853.8A Active CN112991539B (en) | 2021-03-11 | 2021-03-11 | Three-dimensional simulation method for explosive pile and block size distribution based on discrete elements |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112991539B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113743020A (en) * | 2021-09-17 | 2021-12-03 | 中铁十九局集团矿业投资有限公司 | Method and system for establishing displacement field migration path in explosive stack |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104239637A (en) * | 2014-09-16 | 2014-12-24 | 武汉大学 | Method for simulating discrete element muck pile form |
| CN109242144A (en) * | 2018-08-01 | 2019-01-18 | 京工博创（北京）科技有限公司 | A kind of surface blasting prediction technique, apparatus and system |
| CN111199108A (en) * | 2020-01-07 | 2020-05-26 | 包钢集团矿山研究院（有限责任公司） | Blasting blockiness prediction method of gridding segmentation numerical simulation blasting pile model |
| 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 CN202110265853.8A patent/CN112991539B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104239637A (en) * | 2014-09-16 | 2014-12-24 | 武汉大学 | Method for simulating discrete element muck pile form |
| CN109242144A (en) * | 2018-08-01 | 2019-01-18 | 京工博创（北京）科技有限公司 | A kind of surface blasting prediction technique, apparatus and system |
| CN111199108A (en) * | 2020-01-07 | 2020-05-26 | 包钢集团矿山研究院（有限责任公司） | Blasting blockiness prediction method of gridding segmentation numerical simulation blasting pile model |
| 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 |
|---|
| 冯春等: "基于连续-非连续单元方法的露天矿三维台阶爆破全过程数值模拟", 《爆炸与冲击》 * |
| 苏都都等: "露天台阶爆破爆堆形态的PFC模拟", 《爆破》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113743020A (en) * | 2021-09-17 | 2021-12-03 | 中铁十九局集团矿业投资有限公司 | Method and system for establishing displacement field migration path in explosive stack |
| CN113743020B (en) * | 2021-09-17 | 2024-05-24 | 中铁十九局集团矿业投资有限公司 | Method and system for establishing migration path of displacement field in detonation reactor |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112991539B (en) | 2022-07-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8903693B2 (en) | Boundary handling for particle-based simulation | |
| CN105590339B (en) | The collision checking method and system of threedimensional model | |
| CN101615305B (en) | Method and device for detecting collision | |
| CN106777822B (en) | Two-dimentional high-content particle enhancing composite material three-phase meso-mechanical model and method for building up | |
| CN106934826B (en) | Rock slope structure refined modeling and block identification method | |
| Wang et al. | Adaptive image-based intersection volume | |
| CN104504760B (en) | The method and system of real-time update 3-D view | |
| CN112991539B (en) | Three-dimensional simulation method for explosive pile and block size distribution based on discrete elements | |
| CN101937453B (en) | 3-D model (three-dimensional model) search method based on weighted topology constraint comparison | |
| CN110162849B (en) | Modeling method of hybrid fiber concrete | |
| CN103065353A (en) | Three-dimensional model feature extraction method and system and three-dimensional model retrieval method and system | |
| CN115958593A (en) | Intelligent breaking-in robot collision detection method and system based on dual vectors and computer storage medium | |
| Hermann et al. | Ray-traced collision detection for deformable bodies | |
| CN112288874B (en) | Gamma radiation modeling calculation simulation method based on CAD model and Boolean operation | |
| CN113674406B (en) | Three-dimensional simulation method for surface rock mass of open-air deep hole bench blasting pile | |
| CN113284241B (en) | Three-dimensional simulation method for surface rock mass of blasting pile with pile shape distributed in Weibull mode | |
| Guérin et al. | Efficient modeling of entangled details for natural scenes | |
| CN112991538B (en) | Three-dimensional simulation method for explosive stack body based on Weibull distribution of discrete element pair stack shape | |
| TW200823803A (en) | Method for constructing triangular grids of point clouds | |
| CN105930568A (en) | Method for constructing particle cluster discrete element model of convex polyhedral aggregate in any shape | |
| CN107610231A (en) | A kind of dynamic crash detection method | |
| CN109753726A (en) | A kind of globe mill medium movement simulating method based on Bounding Box searching method and GPU | |
| CN114186434A (en) | Asphalt mixture three-dimensional microscopic structure model construction method based on step-by-step division | |
| CN117195382B (en) | Construction method of concrete microscopic model | |
| CN111080790A (en) | Concave body generation method and device based on convex hull algorithm |
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 |