CN111415346A - Particle model damage characterization method based on gridding - Google Patents
Particle model damage characterization method based on gridding Download PDFInfo
- Publication number
- CN111415346A CN111415346A CN202010205959.4A CN202010205959A CN111415346A CN 111415346 A CN111415346 A CN 111415346A CN 202010205959 A CN202010205959 A CN 202010205959A CN 111415346 A CN111415346 A CN 111415346A
- Authority
- CN
- China
- Prior art keywords
- damage
- grid
- model
- grid unit
- unit
- 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
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- 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/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a particle model damage characterization method based on gridding, which comprises the following steps: (1) establishing a rock particle model in particle flow software and calculating; (2) establishing a polygonal mesh covering the rock particle model; (3) grouping contacts and Fratures in the model by using the established grid units, and determining the number of the contacts and Fratures in each grid unit; (4) calculating damage values in each grid unit; (5) exporting the center coordinates of each grid unit and the corresponding damage value to post-processing software to obtain the damage distribution of the model; (6) and counting the number of the grid units in each damage interval. The invention realizes the quantitative characterization of the damage size and distribution, and further expands the application range of the particle flow method.
Description
Technical Field
The invention relates to a damage characterization method, in particular to a particle model damage characterization method based on gridding, and belongs to the field of rock mechanics.
Background
The damage is taken as a parameter for quantitatively reflecting the macroscopic damage condition of the rock mass, has definite meaning, is convenient to understand and is widely used in representing the damage of the engineering rock mass.
The particle flow method uses a large number of particles to construct a complex model, and can accurately simulate the microscopic damage characteristics and the fracture development condition of the model through the contact transfer force and displacement among the particles. The calculation result of the particle flow method is a direct representation of the micro-interaction among particles, and is not a quantitative characterization of model macroscopic damage, so that the method is difficult to quantitatively characterize the distribution and the size of the damage.
How to quantitatively characterize the distribution and the size of damage on the basis of a particle flow calculation result enables a particle flow method to not only characterize the microscopic damage of a model, but also characterize the macroscopic damage of an engineering rock mass, and the method has important significance for further expanding the application range of the method.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the problem that the existing method is difficult to quantitatively characterize model damage, the invention provides a gridding-based particle model damage characterization method.
The technical scheme is as follows: the invention relates to a particle model damage characterization method based on gridding, which comprises the following steps:
(1) establishing a rock particle model in particle flow software and carrying out numerical calculation;
(2) establishing a polygonal mesh covering the particle model;
(3) using a mesh unit of a polygonal mesh to group the Contact and the frame of the particle model, and determining the number of the Contact and the frame in each mesh unit;
(4) calculating damage values in each grid unit;
(5) exporting the center coordinates of each grid unit and the corresponding damage value to post-processing software to obtain the damage distribution of the model;
(6) and counting the number of the grid units in each damage interval so as to represent the damage condition of the model.
Preferably, in the step (2), the method for establishing the polygon mesh includes:
A. setting the shape of grid units in the polygonal grid and the number of grid units in the transverse direction and the longitudinal direction according to the size of the model, and determining the side length of each grid unit;
B. determining the vertex coordinates of each grid unit according to the side length of each grid unit and the position of each grid unit in the model;
C. respectively establishing corresponding grid units in different Geometry sets by using the vertex coordinates of each grid unit;
D. and B-C is executed in a circulating mode until the whole model is covered by the polygon mesh.
In the step (3), the method for grouping the Contact and the frame of the grain model by using the grid unit comprises the following steps: dividing each grid unit into a group, traversing each grid unit, and when the number of intersection points of a ray pointing to any direction from the central point of the Contact or the frame and one grid unit is odd, enabling the corresponding Contact or the frame to be in the grid unit and dividing the grid unit into the group where the grid unit is located; when the number of the intersection points is even, the corresponding Contact or frame is not in the grid cell.
In step (4), the damage value D in a certain grid cell can be calculated according to the following formula:
in the formula, NcIs the number of contacts, N, in the grid cellfThe number of frames in the grid cell.
In the step (6), specifically, a numerical interval may be set to establish a numerical interval of the damage value, the number of the grid cells in each damage interval is counted and a damage distribution histogram is drawn according to the magnitude of the damage value of each grid cell in the model, and the damage degree of the model and the distribution of the damage of different degrees are quantitatively represented.
Has the advantages that: compared with the prior art, the invention has the advantages that: the method calculates the damage of different positions of the model by a gridding method, leads the obtained damage size and position to post-processing software to obtain the damage distribution of the whole model, and finally counts the number of damage units in each damage interval, thereby realizing the quantitative characterization of the damage size and distribution in the particle model and playing an active role in further expanding the application range of the particle flow method.
Drawings
FIG. 1 is a flow chart of a particle model damage characterization method based on gridding according to the present invention;
FIG. 2 is a rock drilling model established in the example;
FIG. 3 shows the development of the fracture in the model of the example;
FIG. 4 is a diagram of a polygon mesh built for a model in an embodiment;
FIG. 5 is a cloud of lesion distributions of the model in an example;
FIG. 6 is a histogram of distribution of damage values of the model in the example.
Detailed Description
The technical scheme of the invention is further explained by combining the attached drawings.
Referring to fig. 1, the method for characterizing particle model damage based on gridding of the present invention includes the following steps:
(1) establishing a rock particle model in particle flow software and carrying out numerical calculation;
(2) establishing a polygonal mesh covering the particle model;
the specific establishment method comprises the following steps:
A. setting the shape of grid units in the polygonal grid and the number of grid units in the transverse direction and the longitudinal direction according to the size of the model, and determining the side length of each grid unit;
B. determining the vertex coordinates of each grid unit according to the side length of each grid unit and the position of each grid unit in the model;
C. respectively establishing corresponding grid units in different Geometry sets by using the vertex coordinates of each grid unit;
D. and B-C is executed in a circulating mode until the whole model is covered by the polygon mesh.
(3) Using a mesh unit of a polygonal mesh to group the Contact and the frame of the particle model, and determining the number of the Contact and the frame in each mesh unit;
the grouping method comprises the following steps: dividing each grid unit into a group, traversing each grid unit, and when the number of intersection points of a ray pointing to any direction from the central point of the Contact or the frame and one grid unit is odd, enabling the corresponding Contact or the frame to be in the grid unit and dividing the grid unit into the group where the grid unit is located; when the number of the intersection points is even, the corresponding Contact or frame is not in the grid cell.
And after grouping is finished, determining the quantity of contacts and frames in each grid unit according to the grouping condition.
(4) Calculating damage values in each grid unit;
the damage value D in a certain grid cell is calculated according to the following formula:
in the formula, NcIs the number of contacts, N, in the grid cellfThe number of frames in the grid cell.
(5) Exporting the center coordinates of each grid unit and the corresponding damage value to post-processing software to obtain the damage distribution of the model;
(6) and counting the number of the grid units in each damage interval so as to represent the damage condition of the model.
And establishing a numerical interval of the damage values according to a certain numerical interval, counting the number of the grid units in each damage interval according to the damage value of each grid unit in the model, drawing a damage distribution histogram, and quantitatively representing the damage degree of the model and the distribution condition of the damage of different degrees.
Examples
By taking the case of obtaining the damage condition of the rock under the impact of the rock drilling rod as an example, the method for characterizing the damage of the particle model based on gridding is explained in detail.
1. And establishing a particle model of the rock drilling rod impacting the rock in particle flow software, and performing numerical calculation. The particle model is shown in fig. 2, and the fracture development of the rock under the impact of a rock drilling rod is shown in fig. 3;
2. setting the number of grid units in the transverse and longitudinal directions, and establishing a polygonal grid by combining the size of the rock particle model, wherein the established polygonal grid is as shown in figure 4; and then grouping the contacts and the frame in the model by using the polygonal mesh, namely grouping the contacts and the frame in the model by using mesh units of 40 meshes in the horizontal direction and 40 meshes in the vertical direction of 1600 meshes in FIG. 4.
3. Calculating damage values in each grid unit by counting the number of contacts and Fratures in each grid unit, and exporting the center coordinates of each grid unit and the corresponding damage values to post-processing software according to a format which can be read by the post-processing software to obtain the damage distribution of the model, as shown in FIG. 5; the present example derives the results into Tecplot, and thus in a readable format of Tecplot.
4. The number of grid cells in each damage interval is counted to obtain a damage distribution histogram, as shown in fig. 6. The figure shows that the number of grids with higher damage values is less, which indicates that the fracture of the rock is mainly developed by cracks and is locally crushed under the impact action of the rock drilling rod.
Claims (5)
1. A particle model damage characterization method based on gridding is characterized by comprising the following steps:
(1) establishing a rock particle model in particle flow software and carrying out numerical calculation;
(2) establishing a polygonal mesh covering the particle model;
(3) using a mesh unit of a polygonal mesh to group the Contact and the frame of the particle model, and determining the number of the Contact and the frame in each mesh unit;
(4) calculating damage values in each grid unit;
(5) exporting the center coordinates of each grid unit and the corresponding damage value to post-processing software to obtain the damage distribution of the model;
(6) and counting the number of the grid units in each damage interval so as to represent the damage condition of the model.
2. The mesh-based particle model damage characterization method according to claim 1, wherein in the step (2), the polygonal mesh is established by:
A. setting the shape of grid units in the polygonal grid and the number of grid units in the transverse direction and the longitudinal direction according to the size of the model, and determining the side length of each grid unit;
B. determining the vertex coordinates of each grid unit according to the side length of each grid unit and the position of each grid unit in the model;
C. respectively establishing corresponding grid units in different Geometry sets by using the vertex coordinates of each grid unit;
D. and B-C is executed in a circulating mode until the whole model is covered by the polygon mesh.
3. The method for characterizing particle model damage based on gridding according to claim 1, wherein in step (3), the method for grouping Contact and frame of the particle model by using the grid cells comprises:
dividing each grid unit into a group, traversing each grid unit, and when the number of intersection points of a ray pointing to any direction from the central point of the Contact or the frame and one grid unit is odd, enabling the corresponding Contact or the frame to be in the grid unit and dividing the grid unit into the group where the grid unit is located; when the number of the intersection points is even, the corresponding Contact or frame is not in the grid cell.
4. The method for characterizing particle model damage based on gridding according to claim 1, wherein in step (4), the damage value D in a certain grid unit is calculated according to the following formula:
in the formula, NcIs the number of contacts, N, in the grid cellfThe number of frames in the grid cell.
5. The method for characterizing particle model damage based on gridding according to claim 1, wherein in the step (6), a numerical interval is set to establish a numerical interval of damage values, the number of grid cells in each damage interval is counted according to the damage value of each grid cell in the model, a damage distribution histogram is drawn, and the damage degree of the model and the distribution of damage of different degrees are quantitatively characterized.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010205959.4A CN111415346B (en) | 2020-03-23 | 2020-03-23 | Particle model damage characterization method based on gridding |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010205959.4A CN111415346B (en) | 2020-03-23 | 2020-03-23 | Particle model damage characterization method based on gridding |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111415346A true CN111415346A (en) | 2020-07-14 |
CN111415346B CN111415346B (en) | 2022-02-11 |
Family
ID=71491342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010205959.4A Active CN111415346B (en) | 2020-03-23 | 2020-03-23 | Particle model damage characterization method based on gridding |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111415346B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113720992A (en) * | 2021-07-12 | 2021-11-30 | 河海大学 | Method for simulating influence of rainfall effect on rock-soil body by raindrop falling method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016003080A1 (en) * | 2014-07-01 | 2016-01-07 | 한국화학연구원 | Real-time inhalation toxicity testing device using lung model |
US20160202164A1 (en) * | 2004-03-06 | 2016-07-14 | Michael Trainer | Methods and apparatus for determining characteristics of particles from scattered light |
CN109241216A (en) * | 2018-07-31 | 2019-01-18 | 国网河南省电力公司南阳供电公司 | The method for realizing grid-based management using power grid hierarchical storage database schema |
CN110415241A (en) * | 2019-08-02 | 2019-11-05 | 同济大学 | A kind of surface of concrete structure quality determining method based on computer vision |
CN110442922A (en) * | 2019-07-15 | 2019-11-12 | 郑州大学 | Stainless armored concrete carefully sees the method for building up of numerical model |
-
2020
- 2020-03-23 CN CN202010205959.4A patent/CN111415346B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160202164A1 (en) * | 2004-03-06 | 2016-07-14 | Michael Trainer | Methods and apparatus for determining characteristics of particles from scattered light |
WO2016003080A1 (en) * | 2014-07-01 | 2016-01-07 | 한국화학연구원 | Real-time inhalation toxicity testing device using lung model |
CN109241216A (en) * | 2018-07-31 | 2019-01-18 | 国网河南省电力公司南阳供电公司 | The method for realizing grid-based management using power grid hierarchical storage database schema |
CN110442922A (en) * | 2019-07-15 | 2019-11-12 | 郑州大学 | Stainless armored concrete carefully sees the method for building up of numerical model |
CN110415241A (en) * | 2019-08-02 | 2019-11-05 | 同济大学 | A kind of surface of concrete structure quality determining method based on computer vision |
Non-Patent Citations (2)
Title |
---|
HE CAN 等: "A New Model for Structural Damage Assessment Using Adaptive Mutation Particle Swarm Optimization and Support Vector Machine", 《IEEE》 * |
蔡改贫 等: "基于多尺度内聚颗粒模型的挤压分级破碎研究", 《岩石力学与工程学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113720992A (en) * | 2021-07-12 | 2021-11-30 | 河海大学 | Method for simulating influence of rainfall effect on rock-soil body by raindrop falling method |
CN113720992B (en) * | 2021-07-12 | 2022-08-02 | 河海大学 | Method for simulating influence of rainfall effect on rock-soil body by raindrop falling method |
Also Published As
Publication number | Publication date |
---|---|
CN111415346B (en) | 2022-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105139444B (en) | Three dimensional particles structural remodeling method based on rock core Particles in Two Dimensions image | |
CN107767453B (en) | Building LIDAR point cloud reconstruction optimization method based on rule constraint | |
CN109658431B (en) | Rock mass point cloud plane extraction method based on region growth | |
CN110675508A (en) | BIM model geometric figure simplification method | |
CN112819962B (en) | Non-uniform grid division and local grid density method in digital image correlation | |
CN103701466A (en) | Scattered point cloud compression algorithm based on feature reservation | |
CN108399290B (en) | Fracture flow MODFLOW simulation method based on fracture network connectivity | |
CN112507581B (en) | Rock mesoscopic numerical model generation method and system capable of avoiding malformation computing unit | |
CN111415346B (en) | Particle model damage characterization method based on gridding | |
CN104166776A (en) | Transmission line conductor form-finding method based on ANSYS | |
CN103106305A (en) | Space grid structure model step-by-step correction method based on actual measurement mode | |
CN106644868A (en) | Measuring method for peripheral interface concentration of two-dimensional non-convex random aggregates | |
CN112069704A (en) | Method for constructing multi-degree-of-freedom lattice model of real concrete specimen | |
CN108460832A (en) | A kind of shell extraction method based on Building Information Model | |
CN103177165A (en) | Coach body structure design system, coach side overturning simulation test system and method | |
CN109446573A (en) | A method of building multidimensional road surface simulation model | |
CN114186434A (en) | Asphalt mixture three-dimensional microscopic structure model construction method based on step-by-step division | |
CN113673186A (en) | Cartesian grid rapid generation method based on STL file | |
CN111291497A (en) | Improved digital core modeling method | |
CN113343423A (en) | Random fracture network generation method based on intensity spatial variability | |
CN104794745A (en) | 3D and isogeometric mixed unit modeling method of rifling barrel | |
Menzel | Evolvable free-form deformation control volumes for evolutionary design optimization | |
CN106934729A (en) | Building Testing and appraisal method and device | |
CN103714192A (en) | Adaptive R-tree based large-data-volume three-dimensional railway design model rendering method | |
CN109992841B (en) | Blasting block space omnibearing partition coupling numerical simulation method |
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 |