CN106560814A - Discrete element based three dimensional form characteristic aggregate forming method - Google Patents
Discrete element based three dimensional form characteristic aggregate forming method Download PDFInfo
- Publication number
- CN106560814A CN106560814A CN201610029558.1A CN201610029558A CN106560814A CN 106560814 A CN106560814 A CN 106560814A CN 201610029558 A CN201610029558 A CN 201610029558A CN 106560814 A CN106560814 A CN 106560814A
- Authority
- CN
- China
- Prior art keywords
- ball
- axis
- materials
- plane
- disruption
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a discrete element based three dimensional form characteristic aggregate forming method; the method comprises the following steps: randomly selecting real aggregates; collecting long axis, short axis, middle axis and fragmentation surface information of the aggregate; cutting a ball particle according to the aggregate information characteristics and certain rules; using regularly ranked small particle size balls to fill zones formed by cutting the ball particle, thus representing true aggregates. The discrete element based three dimensional form characteristic aggregate forming method can select the long axis, short axis, middle axis and fragmentation surface parameters of the aggregate according to the real aggregate statistics data and normal distribution rules, thus representing grading, corner angle and fragmentation surface characteristics of real aggregates.
Description
Technical field
The present invention relates to asphalt three-dimensional emulation technology, especially, a kind of three-dimensional configuration feature based on discrete element gather materials generation method,.
Background technology
Gathered materials using the Numerical Method Study of discrete element or asphalt mechanical property when, how to ensure numeral gather materials can represent the morphological characteristic truly gathered materials be key.The morphological characteristic gathered materials is obtained based on the conventional program by series of complex such as indoor shaping test piece, machinery section, Digital Image Processing, but the two dimension that this mode can only obtain Limited information is gathered materials, and differs greatly with practical situation.Or on the basis of molding indoors, though can obtain the image information that three-dimensional is gathered materials by CT scan technology, this equipment is sufficiently expensive, and needs the substantial amounts of test specimen of molding when required asphalt type is different, and consumes many energy.When noticeable, test specimen of the CT scan to density difference very little of gathering materials, it is difficult to differentiate the boundary characteristic for gathering materials, this undoubtedly increases the difficulty of Digital Image Processing.Therefore, with the development of computer technology, gathered materials especially with distinct element method research using technology of numerical simulation or the mechanical property of asphalt becomes hot issue.
The elementary cell provided by the software of Three Dimensional Discrete Element only has preferable ball granule, represents truly gathering materials with complex shapes such as corner angle, texture and the planes of disruption with Utopian ball, and this there will naturally be huge difference.Have scholar by substantial amounts of it is assumed that the plane of disruption of granule is limited to tetrahedron, pentahedron or hexahedron, and exist it is some be mutually perpendicular to plane, right angle etc., trial is contacted with foundation of actually gathering materials.This method can set up level with relation with truly gathering materials well, but for there is certain difference in morphological characteristic, in addition, this method is when test specimen or bitumen mixture specimen are virtually gathered materials in generation, it is difficult to reach higher degree of compaction by arranging action of gravity and applying external force effect.Even if reaching larger degree of compaction, the primary stress inside test specimen is excessive, has a strong impact on the carrying out of emulation.
The content of the invention
Goal of the invention：In order to overcome the deficiencies in the prior art; the present invention provides a kind of three-dimensional configuration feature based on discrete element and gathers materials generation method; the plane of disruption is obtained by Cutting ball granule; and fill the method for path grain ball and obtain reaction three-dimensional configuration feature and gather materials, solve the deficiencies in the prior art.
Technical scheme：For achieving the above object, the technical solution used in the present invention is：A kind of three-dimensional configuration feature based on discrete element is gathered materials generation method, it is characterised in that statistics is truly gathered materials parameter, and the initial ball granule of the simulation of generation is cut according to the parameter of gathering materials, some facets are formed；And the space between facet and initial ball granule is filled with the small particle ball of queueing discipline；The parameter of gathering materials includes major axis, short axle, axis and the plane of disruption；
The Grading feature for truly gathering materials is represented with major axis, short axle and axis, the plane of disruption is represented with the facet；Corner angle are formed between facet, or facet forms corner angle with the sphere of ball granule；Parameter of gathering materials truly is gathered materials by calculating acquisition according to what is randomly selected.
Further, the method specifically includes following steps：
1) randomly select in truly gathering materials, gathered materials at random, parameter of gathering materials described in statistical analysiss at random, statistics the major axis for gathering materials at random, short axle, axis, calculate average and corresponding standard deviation respectively；
2) matched somebody with somebody according to the level truly gathered materials, the initial ball granule of simulation is generated using normal distyribution function；
3) according to step 1) and step 2), calculate and obtain the plane of disruption；
4) according to step 3) plane of disruption that obtains cut to initial ball granule, forms multiple facets, while forming multiple corner angle after repeatedly cutting；
5) with the space between the small particle ball filling facet and initial ball granule of queueing discipline；After completing filling, initial ball is deleted, by gathering materials for the small particle ball reaction of formation three-dimensional configuration feature of filling.
Further, calculate the plane of disruption to comprise the following steps that：
2R=la+N(0,1)×Sa (1)
2M=lb+N(0,1)×Sb (2)
2N=lc+N(0,1)×Sc (3)
(x-x0)2+(y-y0)2+(z-z0)2=R2 (4)
D=N+ (M-N) × urand (6)
nx1×(x-x0)+ny1×(y-y0)+nz1×(z-z0)=0 (7)
Initially ball centre of sphere expression formula is：O(x0,y0,z0)；Radius is R；The diameter of ball is equal to major axis, and formula (5) is the normal vector of the initial ball, formula (6) for the ball centre of sphere to plane distance range between semi-minor axis and semi-mean axis, formula (7) is with (nx1,ny1,nz1) normal vector is through centre of sphere O (x0,y0,z0) plane equation；
Wherein, calculate major axis average l gathered materials at randoma, axis average lb, short axle average lc；The corresponding standard deviation of average is respectively Sa、Sb、Sc；N (0,1) is normal distribution random number；Urand is uniform random number.
Beneficial effect：A kind of three-dimensional configuration feature based on discrete element that the present invention is provided is gathered materials generation method, the parameters such as the size, the plane of disruption and the corner angle that gather materials are obtained through substantial amounts of data statistic analysis according to the three-dimensional configuration feature truly gathered materials, and are truly gathered materials with good relatedness.And the three-dimensional data for representing feature of gathering materials can be obtained by succinct formula step.
Description of the drawings
Fig. 1 is that initial ball particles generation gathers materials the schematic diagram of facet
Fig. 2 is method of the present invention flow chart
Specific embodiment
Below in conjunction with the accompanying drawings the present invention is further described.
It is illustrated in figure 1 a kind of three-dimensional configuration feature based on discrete element to gather materials generation method, it is characterised in that statistics is truly gathered materials parameter, and the initial ball granule of the simulation of generation is cut according to the parameter of gathering materials, some facets are formed；And the space between facet and initial ball granule is filled with the small particle ball of queueing discipline；The parameter of gathering materials includes major axis, short axle, axis and the plane of disruption；
The Grading feature for truly gathering materials is represented with major axis, short axle and axis, the plane of disruption is represented with the facet；Corner angle are formed between facet, or facet forms corner angle with the sphere of ball granule；Parameter of gathering materials truly is gathered materials by calculating acquisition according to what is randomly selected.
Further, the method specifically includes following steps：
1) randomly select in truly gathering materials, gathered materials at random, parameter of gathering materials described in statistical analysiss at random, statistics the major axis for gathering materials at random, short axle, axis, calculate average and corresponding standard deviation respectively；
2) matched somebody with somebody according to the level truly gathered materials, the initial ball granule of simulation is generated using normal distyribution function；
3) according to step 1) and step 2), calculate and obtain the plane of disruption；As shown in Fig. 2 calculate a plane of disruption comprising the following steps that：
2R=la+N(0,1)×Sa (1)
2M=lb+N(0,1)×Sb (2)
2N=lc+N(0,1)×Sc (3)
(x-x0)2+(y-y0)2+(z-z0)2=R2 (4)
D=N+ (M-N) × urand (6)
nx1×(x-x0)+ny1×(y-y0)+nz1×(z-z0)=0 (7)
Initially ball centre of sphere expression formula is：O(x0,y0,z0)；Radius is R；The diameter of ball is equal to major axis, formula (5) is the normal vector of the initial ball, formula (6) for the ball centre of sphere to plane distance range between semi-minor axis and semi-mean axis, the center for gathering materials to the plane of disruption distance between semi-minor axis and semi-mean axis, formula (7) is with (nx1,ny1,nz1) normal vector is through centre of sphere O (x0,y0,z0) plane equation；
Wherein, calculate major axis average l gathered materials at randoma, axis average lb, short axle average lc；The corresponding standard deviation of average is respectively Sa、Sb、Sc；N (0,1) is normal distribution random number；Urand is uniform random number.
4) according to step 3) plane of disruption that obtains cut to initial ball granule, forms multiple facets, while forming multiple corner angle after repeatedly cutting；
5) with the space between the small particle ball filling facet and initial ball granule of queueing discipline；After completing filling, initial ball is deleted, by gathering materials for the small particle ball reaction of formation three-dimensional configuration feature of filling.Can characterize the level truly gathered materials match somebody with somebody, corner angle and the plane of disruption etc. feature.
The above is only the preferred embodiment of the present invention, it should be pointed out that：For those skilled in the art, under the premise without departing from the principles of the invention, some improvements and modifications can also be made, these improvements and modifications also should be regarded as protection scope of the present invention.
Claims (3)
1. a kind of three-dimensional configuration feature based on discrete element is gathered materials generation method, it is characterised in that statistics is truly gathered materials parameter,
The initial ball granule of the simulation of generation is cut according to the parameter of gathering materials, some facets are formed；And with arrangement
Space between the small particle ball filling facet of rule and initial ball granule；The parameter of gathering materials includes major axis, short
Axle, axis and the plane of disruption；
The Grading feature for truly gathering materials is represented with major axis, short axle and axis, the plane of disruption is represented with the facet；Facet
Between form corner angle, or facet forms corner angle with the sphere of ball granule；Parameter of gathering materials is according to the true collection for randomly selecting
Material is obtained by calculating.
2. a kind of three-dimensional configuration feature based on discrete element as claimed in claim 1 is gathered materials generation method, it is characterised in that
The method specifically includes following steps：
1) randomly select in truly gathering materials, gathered materials at random, parameter of gathering materials described in statistical analysiss at random is counted simultaneously respectively
The major axis for gathering materials at random, short axle, axis, calculate average and corresponding standard deviation；
2) matched somebody with somebody according to the level truly gathered materials, the initial ball granule of simulation is generated using normal distyribution function；
3) according to step 1) and step 2), calculate and obtain the plane of disruption；
4) according to step 3) plane of disruption that obtains cut to initial ball granule, forms multiple cuttings after repeatedly cutting
Face, while forming multiple corner angle；
5) with the space between the small particle ball filling facet and initial ball granule of queueing discipline；After completing filling,
Initial ball is deleted, by gathering materials for the small particle ball reaction of formation three-dimensional configuration feature of filling.
3. a kind of three-dimensional configuration feature based on discrete element as claimed in claim 2 is gathered materials generation method, it is characterised in that
Calculate the plane of disruption to comprise the following steps that：
2R=la+N(0,1)×Sa (1)
2M=lb+N(0,1)×Sb (2)
2N=lc+N(0,1)×Sc (3)
(x-x0)2+(y-y0)2+(z-z0)2=R2 (4)
D=N+ (M-N) × urand (6)
nx1×(x-x0)+ny1×(y-y0)+nz1×(z-z0)=0 (7)
Initially ball centre of sphere expression formula is：O(x0,y0,z0)；Radius is R；The diameter of ball is equal to major axis, formula (5)
For the normal vector of the initial ball, formula (6) for the ball centre of sphere to plane distance range semi-minor axis and semi-mean axis it
Between, formula (7) is with (nx1,ny1,nz1) normal vector is through centre of sphere O (x0,y0,z0) plane equation；
Wherein, calculate major axis average l gathered materials at randoma, axis average lb, short axle average lc；Average is corresponding
Standard deviation is respectively Sa、Sb、Sc；N (0,1) is normal distribution random number；Urand is uniform random number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610029558.1A CN106560814A (en) | 2016-01-15 | 2016-01-15 | Discrete element based three dimensional form characteristic aggregate forming method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610029558.1A CN106560814A (en) | 2016-01-15 | 2016-01-15 | Discrete element based three dimensional form characteristic aggregate forming method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106560814A true CN106560814A (en) | 2017-04-12 |
Family
ID=58485515
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610029558.1A Pending CN106560814A (en) | 2016-01-15 | 2016-01-15 | Discrete element based three dimensional form characteristic aggregate forming method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106560814A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109002651A (en) * | 2018-08-22 | 2018-12-14 | 中国空气动力研究与发展中心超高速空气动力研究所 | The fly-cutting faceted search method decomposed for parallel numerical simulating area |
| CN109753667A (en) * | 2017-11-01 | 2019-05-14 | 长沙理工大学 | A kind of consideration is gathered materials the three-dimensional random number test specimen generation method of feature |
| CN109991401A (en) * | 2019-04-19 | 2019-07-09 | 长沙理工大学 | A kind of consideration is gathered materials the three-dimensional random number test specimen generation method of feature |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101354241A (en) * | 2008-07-11 | 2009-01-28 | 长安大学 | Method and system for evaluating aggregate digital image |
| CN103218480A (en) * | 2013-03-20 | 2013-07-24 | 东南大学 | Method for randomly building bituminous mixture multilayered structure simulation model |
| CN104599321A (en) * | 2015-01-24 | 2015-05-06 | 合肥工业大学 | Building method of discrete element module of true aggregate particles based on X-ray CT image |
-
2016
- 2016-01-15 CN CN201610029558.1A patent/CN106560814A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101354241A (en) * | 2008-07-11 | 2009-01-28 | 长安大学 | Method and system for evaluating aggregate digital image |
| CN103218480A (en) * | 2013-03-20 | 2013-07-24 | 东南大学 | Method for randomly building bituminous mixture multilayered structure simulation model |
| CN104599321A (en) * | 2015-01-24 | 2015-05-06 | 合肥工业大学 | Building method of discrete element module of true aggregate particles based on X-ray CT image |
Non-Patent Citations (3)
| Title |
|---|
| GAO, Y.: "Numerical Simulation of Motion Rules of Coarse Aggregates in the Compaction Process", 《JOURNAL OF TESTING AND EVALUATION》 * |
| WEIDONG LIU: "Evaluation Indices of an Asphalt-Mixture Digital Specimen Based on the Discrete-Element Method", 《JOURNAL OF TESTING AND EVALUATION》 * |
| 张东: "沥青混合料粗集料的形态特征研究和力学性能的离散元模拟", 《万方学术论文库》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109753667A (en) * | 2017-11-01 | 2019-05-14 | 长沙理工大学 | A kind of consideration is gathered materials the three-dimensional random number test specimen generation method of feature |
| CN109002651A (en) * | 2018-08-22 | 2018-12-14 | 中国空气动力研究与发展中心超高速空气动力研究所 | The fly-cutting faceted search method decomposed for parallel numerical simulating area |
| CN109002651B (en) * | 2018-08-22 | 2019-04-09 | 中国空气动力研究与发展中心超高速空气动力研究所 | The fly-cutting faceted search method decomposed for parallel numerical simulating area |
| CN109991401A (en) * | 2019-04-19 | 2019-07-09 | 长沙理工大学 | A kind of consideration is gathered materials the three-dimensional random number test specimen generation method of feature |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103218480B (en) | Method for randomly building bituminous mixture multilayered structure simulation model | |
| Nie et al. | Probabilistic characterization and simulation of realistic particle shape based on sphere harmonic representation and Nataf transformation | |
| CN106560814A (en) | Discrete element based three dimensional form characteristic aggregate forming method | |
| CN106934826B (en) | Rock slope structure refined modeling and block identification method | |
| CN109164117A (en) | Based on-bituminous concrete damage forecast the method for mortar interface Fracture of gathering materials | |
| CN108491677A (en) | Pore character statistical method based on the micro pore model for improving maximum ball | |
| CN104573705A (en) | Clustering method for building laser scan point cloud data | |
| Hongxin et al. | Study on multi-size seed-metering device for vertical plate soybean precision planter | |
| CN109241646B (en) | Multi-factor two-dimensional soil-rock mixture generation method based on elliptical stacking and random field | |
| CN104216974B (en) | The method of unmanned plane images match based on words tree Block Cluster | |
| Zheng et al. | Characterization of two-and three-dimensional morphological properties of fragmented sand grains | |
| CN109509251A (en) | Multifactor three-dimensional soil-rock mixture generation method | |
| Qing et al. | Impact factors of soil wind erosion in the center of Taklimakan Desert | |
| CN106769322A (en) | A kind of prismatical joint rock mass preparation method of sample of cylinder containing stomata | |
| CN109632429B (en) | Sample preparation method for soil-rock mixture double-shaft compression test | |
| CN104123428A (en) | CMP process simulation method and system | |
| Yan et al. | A review of the application of discrete element method in agricultural engineering: A case study of soybean | |
| CN109191571B (en) | Method for preparing mechanical test standard aggregate by applying 3D printing technology | |
| Sathyaprakash et al. | Morphology of clusters and superclusters in N-body simulations of cosmological gravitational clustering | |
| CN107845140A (en) | A kind of three-dimensional modeling method and threedimensional model | |
| Jin et al. | Bayesian analysis for exponential random graph models using the adaptive exchange sampler | |
| Hao et al. | Calibration of simulation parameters for wind erosion gas-solid two-phase flow in arid and semiarid soils | |
| Guérin et al. | Efficient modeling of entangled details for natural scenes | |
| Dun et al. | Model construction of soybean average diameter and hole parameters of seed-metering wheel based on DEM | |
| Liu et al. | A novel stochastic generation method of shale blocks in S-RM considering the particle size effect of morphological features |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170412 |
|
| RJ01 | Rejection of invention patent application after publication |