CN105160052A - Displacement vector filling method based on energetic material approximate spherical filling maximization - Google Patents
Displacement vector filling method based on energetic material approximate spherical filling maximization Download PDFInfo
- Publication number
- CN105160052A CN105160052A CN201510390871.3A CN201510390871A CN105160052A CN 105160052 A CN105160052 A CN 105160052A CN 201510390871 A CN201510390871 A CN 201510390871A CN 105160052 A CN105160052 A CN 105160052A
- Authority
- CN
- China
- Prior art keywords
- axis
- particle
- small
- negative direction
- move
- 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
Abstract
The invention requests to protect a displacement vector filling method based on energetic material approximate spherical filling maximization. The method comprises the following steps: randomly putting spherical particles into a compression cabin according to a radius from big to small, and entering a regulation stage after a spherical particle putting stage is finished; and dividing the regulation stage into three processes, i.e., the spherical particles moves to an x-axis negative direction from small to big according to an x-axis coordinate, moves to a y-axis negative direction from small to big according to a y-axis coordinate and moves to a z-axis negative direction from small to big according to a z-axis coordinate. After the particle movement, gaps among the particles are reduced to a high limit, and the cabin can be continuously filled to increase a filling rate.
Description
Technical field
The invention belongs to materialogy, computer science, be specifically related to a kind of energetic material compression field solve differ in size and the random spherical particles in position is filled time dispersion, excesssive gap, causes the method for the problem that compression cabin body filling rate is low,
Background technology
Energetic material and energy-containing compound, be called for short energy material, mean the material (HEDM) of high-energy-density; It is characterized by such material and has explosivity, deflagrability or other materials through the large energy of particular excitation condition meeting two-forty height output release more.Energetic material mechanical property mainly refers to that it is under varying environment (temperature, medium, humidity), bear various plus load (stretch, compression, bending, torsion, impact, alterante stress etc.) time the mechanical characteristics that shows.Its character of energetic material is subject to material internal composition, and as various factors such as material shape, volume ratios, the change of material prescription is very large on the impact of overall mechanical property; This material has explosive nature thereof simultaneously, is therefore prepared the cost and risk of experimental study to energetic material all very high.
In recent years, the numerical simulation of energetic material mechanical property has received increasing concern, due to the complicacy of material granule forming process, illustrates its compacting mechanism very difficult, particularly how simulation material distribution of particles, to experimentation and result data particularly important.This simulation adopts Univ Utah USA C-SAFE (CenterfortheSimulationofAccidentalFiresandExplosions) research centre to develop UCF (UintahComputationalFramework) large-scale parallel visual numerical simulation instrument Uintah and VisIt software.
From research conditions, how setting up spherical particles microscopical structure model is one of key issue in material granule compression process.Foreign scholar Baer
[1]deng the modeling method using molecular dynamics, establish same size and the three-dimensional blasting explosive granules microscopical structure model of position stochastic distribution, domestic Liu Qun
[2]blasting explosive granules is approximated to two-dimensional circular, and regularly hexagon arrangement.The desirable modeling of above scholar all has certain influence to experimental data.
Simply intergranular space can be caused comparatively large to ball-type energetic material random arrangement, filling rate is low, differs greatly with true powder charge process.
If can in put procedure, the distance that conservative control spherical particles produces thus reduce space, will certainly increase the filling quantity of particle, arrange more tight between particle; Even if the distance of placing between control particle, intergranular gap still exists, if again by certain method, will compress cabin body " rock ", make gap between particle less, vacate certain space simultaneously, just can continue to place, make filling rate reach the highest.
By the analysis to existing domestic data, find also lack by controlling random put procedure at present, intervene particle produce Distance geometry place after make compression cabin body " rock " thus again attempt placement method.The gap produced when the present invention is just by analyzing reality and filling, proposes a kind of simple algorithm effectively solving the low problem of spherical particles filling rate.List of references
[1]BaerMR.ModelingheterogeneousenergeticmaterialsattheMesoscale[J].Thermochimicaacta,2002,384:351-367.
[2] Liu Qun, Chen Lang, Lu Jianying, etc. the compressing numerical simulation of blasting explosive granules [J]. Chinese Journal Of High Pressure Physics 2009,23 (6) .DOI:10.3969/j.issn.1000-5773.2009.
Summary of the invention
For the deficiencies in the prior art, proposing a kind of increases filling rate, improves being similar to ball-type based on energetic material and filling maximized displacement vector fill method of the authenticity of simulated data.Technical scheme of the present invention is as follows: one is similar to ball-type based on energetic material and fills maximized displacement vector fill method, and it comprises the following steps:
101, spherical particles is put into compression cabin body from big to small at random by radius, after completing the spherical particles placement stage, enter the adjusting stage;
102, the adjusting stage is divided into three processes, set up xyz coordinate system, wherein x-axis, y-axis are transverse axis in xy plane or the longitudinal axis, and namely spherical particles moves to x-axis negative direction successively from small to large by x-axis coordinate, moves successively from small to large with y-axis coordinate to y-axis negative direction again, finally moves to z-axis negative direction successively from small to large with z-axis coordinate; Process one: with x-axis coordinate by spherical particles all in the body of cabin with 0, x-axis coordinate for initial point, move to x-axis negative direction successively from small to large by x-axis coordinate,
103, process two: again with 0, y-axis coordinate for initial point, move to y-axis negative direction successively from small to large by y-axis coordinate;
104, process three: last with 0, z-axis coordinate for initial point, move to z-axis negative direction successively from small to large by z-axis coordinate, meet mobile after do not hand over any ball-type Particle Phase.
Further, the stage of placing of the spherical particles in step 101 is specially:
If the radius of spherical particles is r
1, r
2, r
3... r
n, often place a spherical particles, all successively compare with all particles before, if crossing with arbitrary particle, abandon this placement, if be greater than d with the gap of arbitrary particle, abandon this placement equally, often abandon once placing, the placement frequency of failure adds 1; If place the frequency of failure to be greater than maximum attempts, reducing radius is r
m(1<m≤n), n is expressed as n kind radius quantity
,and number of attempt is established 0, again place, place and successfully should meet as lower inequality:
Wherein r
p1,2represent any two particle radius sizes, x
p1,2represent two particle x coordinate sizes, y
p1,2represent two particle y coordinate sizes, z
p1,2represent two particle z coordinate sizes, d represents maximal clearance between any particle.
Further, be specially to the step of x-axis negative direction movement successively from small to large by x-axis coordinate in step 102: establish in the body of compression cabin and amount to n spherical particles, the distance of first time movement is d
1, the distance knots modification of (n-1)th movement is d
i, d
ifor attempting displacement, can according to specific experiment requirements set, generally the smaller the better, meet not crossing with any particle in moving process, move mode is: namely A1, the 1st particle move to x-axis tangent by the position nearest from x-axis, and distance is d
1; A2, the 2nd to the n-th particle, with x-axis coordinate tentative mobile spherical particles successively from small to large, displacement is with d
1for benchmark, if not crossing with arbitrary particle after moving to left, with d
ifor decrement continues to move to left, otherwise stop movement, and do not hand over any ball-type Particle Phase after meeting movement, complete and move to x-axis negative direction.
Further, be specially successively from small to large in step 103 by y-axis coordinate to the step of y-axis negative direction movement: B1, move to y-axis tangent by the 1st particle, distance is d
1; B2, the 2nd to the n-th particle, with z-axis coordinate tentative mobile spherical particles successively from small to large, displacement is with d
1for benchmark, if not crossing with arbitrary particle after reach, with d
ifor decrement continues to move down, otherwise stop movement, and do not hand over any ball-type Particle Phase after meeting movement, complete and move to y-axis negative direction.
Further, be specially successively from small to large in step 104 by z-axis coordinate to z-axis negative direction move mode: C1, move to z-axis tangent by the 1st particle nearest from z-axis, distance is d
1; C2, the 2nd to the n-th particle, with z-axis coordinate tentative mobile spherical particles successively from small to large, displacement is with d
1for benchmark, if not crossing with arbitrary particle after moving down, with d
ifor decrement continues to move down, otherwise stop movement, and do not hand over any ball-type Particle Phase after meeting movement, complete and move to z-axis negative direction.
Advantage of the present invention and beneficial effect as follows:
The placement stage, often generate stochastic particle all can adjacent before the particle that generates, this stage avoids because make Particle Phase every too far away and produce the larger problem in space at random, increases particles filled quantity;
Adjusting stage, through placing the stage, space between particle, is still had to exist.Moved by the horizontal direction of particle in x-axis, y-axis, make between particle more compact, reduce space, again by moving down in z-axis, cabin body space is also vacateed in further minimizing space, capable of circulation enter the placement stage continue filler particles, thus increase filling rate.
Accompanying drawing explanation
Fig. 1 to be preferred embodiment of the present invention Fig. 1 be traditional spherical particles is according to radius random Placement process flow diagram from big to small.
Fig. 2 is the filling algorithm process flow diagram in the stage of filling after improving that the present invention proposes.
Fig. 3 is whole algorithm flow chart, comprises filling stage and adjusting stage.
Fig. 4 is for generally to fill schematic diagram.
Fig. 5 is that invention algorithm fills complete schematic diagram.
Fig. 6 is the complete schematic diagram of general filling.
Fig. 7 is that invention algorithm fills complete schematic diagram.
Embodiment
Below in conjunction with accompanying drawing, the invention will be further described:
Be analysis process of the present invention see Fig. 1, Fig. 3, as shown in the figure, one provided by the invention is similar to ball based on energetic material and fills maximized algorithm, and this algorithm comprises filling stage and adjusting stage, and step is as follows:
The filling stage:
1, spherical particles is put into compression cabin body, if radius is r from big to small at random by radius
1, r
2, r
3... r
nif maximum attempts is 1,000,000 times, maximal clearance is that d, p represent spherical particles;
2, often place a spherical particles (except first), all can with before so particle successively compares, if crossing with arbitrary particle, abandon this placement, if be greater than d with the distance of arbitrary particle, abandon this equally to place, often abandon once placing, place the frequency of failure and add 1; If place the frequency of failure to be greater than maximum attempts, changing radius is r
m(1<m≤n), and number of attempt is established 0, again place.Place and successfully should meet as lower inequality;
If 3 reach maximum placement number of times and radius is r
n, then place complete and quit a program;
Adjusting stage:
4, by spherical particles all in the body of cabin with x-axis coordinate for initial point, move to x-axis negative direction successively from small to large, if amount to n spherical particles in compression cabin body, the distance of movement is for the first time d
1, the distance knots modification of (n-1)th movement is d
i, move mode is:
4.1, the 1st particle (nearest from x-axis) moves to x-axis tangent, and distance is d
1;
4.2, the 2nd to the n-th particle, and with x-axis coordinate tentative mobile spherical particles successively from small to large, displacement is with d
1for benchmark, if not crossing with arbitrary particle after moving to left, with d
ifor decrement continues to move to left, otherwise stop movement, and do not hand over any ball-type Particle Phase after meeting movement;
5, treat that all spherical particles are by after as above process moves, then move to y-axis negative direction successively from small to large with y-axis coordinate, move mode is:
5.1, the 1st particle (nearest from y-axis) moves to y-axis tangent, and distance is d
1;
5.2, the 2nd to the n-th particle, and with y-axis coordinate tentative mobile spherical particles successively from small to large, displacement is with d
1for benchmark, if not crossing with arbitrary particle after reach, with d
ifor decrement continues reach, otherwise stop mobile, and meet mobile after do not hand over any ball-type Particle Phase;
6, after as above process moves, then move to z-axis negative direction successively from small to large with z-axis coordinate, move mode is:
5.1, the 1st particle (nearest from z-axis) moves to z-axis tangent, and distance is d
1;
5.2, the 2nd to the n-th particle, and with z-axis coordinate tentative mobile spherical particles successively from small to large, displacement is with d
1for benchmark, if not crossing with arbitrary particle after moving down, with d
ifor decrement continues to move down, otherwise stop movement, and do not hand over any ball-type Particle Phase after meeting movement;
7, after completing the adjusting stage, there will be continuous white space in the body of cabin, get back to the filling stage, continue to perform successively, the maximization that can realize compressing spherical particles in the body of cabin is placed;
These embodiments are interpreted as only being not used in for illustration of the present invention limiting the scope of the invention above.After the content of reading record of the present invention, technician can make various changes or modifications the present invention, and these equivalence changes and modification fall into the scope of the claims in the present invention equally.
Claims (5)
1. be similar to ball-type based on energetic material and fill a maximized displacement vector fill method, it is characterized in that, comprise the following steps:
101, spherical particles is put into compression cabin body from big to small at random by radius, after completing the spherical particles placement stage, enter the adjusting stage;
102, the adjusting stage is divided into three processes, set up xyz coordinate system, wherein x-axis, y-axis are transverse axis in xy plane or the longitudinal axis, and namely spherical particles moves to x-axis negative direction successively from small to large by x-axis coordinate, moves successively from small to large with y-axis coordinate to y-axis negative direction again, finally moves to z-axis negative direction successively from small to large with z-axis coordinate; Process one: with x-axis coordinate by spherical particles all in the body of cabin with 0, x-axis coordinate for initial point, move to x-axis negative direction successively from small to large by x-axis coordinate,
103, process two: again with 0, y-axis coordinate for initial point, move to y-axis negative direction successively from small to large by y-axis coordinate;
104, process three: last with 0, z-axis coordinate for initial point, move to z-axis negative direction successively from small to large by z-axis coordinate, meet mobile after do not hand over any ball-type Particle Phase.
2. be according to claim 1ly similar to ball-type based on energetic material and fill maximized displacement vector fill method, it is characterized in that, the stage of placing of the spherical particles in step 101 is specially:
If the radius of spherical particles is r
1, r
2, r
3... r
n, often place a spherical particles, all successively compare with all particles before, if crossing with arbitrary particle, abandon this placement, if be greater than d with the gap of arbitrary particle, abandon this placement equally, often abandon once placing, the placement frequency of failure adds 1; If place the frequency of failure to be greater than maximum attempts, reducing radius is r
m(1<m≤n), n is expressed as n kind radius quantity
,and number of attempt is established 0, again place, place and successfully should meet as lower inequality:
Wherein r
p1,2represent any two particle radius sizes, x
p1,2represent two particle x coordinate sizes, y
p1,2represent two particle y coordinate sizes, z
p1,2represent two particle z coordinate sizes, d represents maximal clearance between any particle.
3. be according to claim 1ly similar to ball-type based on energetic material and fill maximized displacement vector fill method, it is characterized in that, be specially to the step of x-axis negative direction movement successively from small to large by x-axis coordinate in step 102: establish in the body of compression cabin and amount to n spherical particles, the distance of first time movement is d
1, the distance knots modification of (n-1)th movement is d
i, d
ifor attempting displacement, can according to specific experiment requirements set, generally the smaller the better, meet not crossing with any particle in moving process, move mode is: namely A1, the 1st particle move to x-axis tangent by the position nearest from x-axis, and distance is d
1; A2, the 2nd to the n-th particle, with x-axis coordinate tentative mobile spherical particles successively from small to large, displacement is with d
1for benchmark, if not crossing with arbitrary particle after moving to left, with d
ifor decrement continues to move to left, otherwise stop movement, and do not hand over any ball-type Particle Phase after meeting movement, complete and move to x-axis negative direction.
4. be according to claim 1ly similar to ball-type based on energetic material and fill maximized displacement vector fill method, it is characterized in that, be specially to the step of y-axis negative direction movement successively from small to large by y-axis coordinate in step 103: B1, move to y-axis tangent by the 1st particle, distance is d
1; B2, the 2nd to the n-th particle, with y-axis coordinate tentative mobile spherical particles successively from small to large, displacement is with d
1for benchmark, if not crossing with arbitrary particle after reach, with d
ifor decrement continues to move down, otherwise stop movement moving forward, and do not hand over any ball-type Particle Phase after meeting movement, complete and move to y-axis negative direction.
5. be according to claim 1ly similar to ball-type based on energetic material and fill maximized displacement vector fill method, it is characterized in that, be specially to z-axis negative direction move mode successively from small to large by z-axis coordinate in step 104: C1, move to z-axis tangent by the 1st particle nearest from z-axis, distance is d
1; C2, the 2nd to the n-th particle, with z-axis coordinate tentative mobile spherical particles successively from small to large, displacement is with d
1for benchmark, if not crossing with arbitrary particle after moving down, with d
ifor decrement continues to move down, otherwise stop movement, and do not hand over any ball-type Particle Phase after meeting movement, complete and move to z-axis negative direction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510390871.3A CN105160052B (en) | 2015-07-03 | 2015-07-03 | Maximized displacement vector fill method is filled based on energetic material approximation ball-type |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510390871.3A CN105160052B (en) | 2015-07-03 | 2015-07-03 | Maximized displacement vector fill method is filled based on energetic material approximation ball-type |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105160052A true CN105160052A (en) | 2015-12-16 |
| CN105160052B CN105160052B (en) | 2018-06-29 |
Family
ID=54800908
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510390871.3A Active CN105160052B (en) | 2015-07-03 | 2015-07-03 | Maximized displacement vector fill method is filled based on energetic material approximation ball-type |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105160052B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105787273A (en) * | 2016-02-26 | 2016-07-20 | 重庆邮电大学 | Computer simulation method for energetic material shock initiation |
| CN106503404A (en) * | 2016-11-29 | 2017-03-15 | 重庆邮电大学 | A kind of PBX blasting explosive granules particle Random Discrete fill method |
| CN111489427A (en) * | 2020-04-07 | 2020-08-04 | 清华大学 | Sphere filling implementation method and computer equipment |
| CN113405966A (en) * | 2021-06-08 | 2021-09-17 | 浙江广天构件集团股份有限公司 | Method for calculating pore size distribution of cement-based material particle accumulation system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8170846B2 (en) * | 2009-04-14 | 2012-05-01 | The Yokohama Rubber Co., Ltd. | Computation method of mechanical material constant of composite material and volume fraction of material component in composite material, and recording medium |
| CN103310034A (en) * | 2012-03-08 | 2013-09-18 | 住友橡胶工业株式会社 | Simulation method for high polymer material |
-
2015
- 2015-07-03 CN CN201510390871.3A patent/CN105160052B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8170846B2 (en) * | 2009-04-14 | 2012-05-01 | The Yokohama Rubber Co., Ltd. | Computation method of mechanical material constant of composite material and volume fraction of material component in composite material, and recording medium |
| CN103310034A (en) * | 2012-03-08 | 2013-09-18 | 住友橡胶工业株式会社 | Simulation method for high polymer material |
Non-Patent Citations (3)
| Title |
|---|
| M.R. BAER: "《Modeling heterogeneous energetic materials at the mesoscale》", 《THERMOCHIMICA ACTA》 * |
| 刘群等: "《炸药颗粒压制成型数值模拟》", 《高压物理学报》 * |
| 黄辉: "《颗粒级配技术及其在含能材料中的应用》", 《含能材料》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105787273A (en) * | 2016-02-26 | 2016-07-20 | 重庆邮电大学 | Computer simulation method for energetic material shock initiation |
| CN105787273B (en) * | 2016-02-26 | 2018-10-12 | 重庆邮电大学 | The computer simulation method of energetic material Impact Initiation |
| CN106503404A (en) * | 2016-11-29 | 2017-03-15 | 重庆邮电大学 | A kind of PBX blasting explosive granules particle Random Discrete fill method |
| CN106503404B (en) * | 2016-11-29 | 2019-09-10 | 重庆邮电大学 | A kind of PBX blasting explosive granules particle Random Discrete fill method |
| CN111489427A (en) * | 2020-04-07 | 2020-08-04 | 清华大学 | Sphere filling implementation method and computer equipment |
| CN111489427B (en) * | 2020-04-07 | 2021-11-23 | 清华大学 | Sphere filling implementation method and computer equipment |
| CN113405966A (en) * | 2021-06-08 | 2021-09-17 | 浙江广天构件集团股份有限公司 | Method for calculating pore size distribution of cement-based material particle accumulation system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105160052B (en) | 2018-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105160052A (en) | Displacement vector filling method based on energetic material approximate spherical filling maximization | |
| Wang et al. | DEM analysis of energy dissipation in crushable soils | |
| CN102628861B (en) | Method for simulating temperature cracking value of mass concrete | |
| Nammi et al. | Finite element analysis of closed-cell aluminium foam under quasi-static loading | |
| Yu et al. | Study of the features of outburst caused by rock cross-cut coal uncovering and the law of gas dilatation energy release | |
| Liu et al. | Numerical simulations of oblique-angle penetration by deformable projectiles into concrete targets | |
| CN112199843A (en) | Rock fracture evolution mesoscopic model and method based on particle structure and electronic equipment | |
| Xie et al. | Design of the ballistic performance of shear thickening fluid (STF) impregnated Kevlar fabric via numerical simulation | |
| Lv et al. | A probabilistic method for determining the volume fraction of pre-embedded capsules in self-healing materials | |
| Yan et al. | Simulation of granular material behaviour using DEM | |
| Zhang et al. | Effects of meso shape irregularity of metal foam on yield features under triaxial loading | |
| Khude et al. | A parallel GPU implementation of the absolute nodal coordinate formulation with a frictional/contact model for the simulation of large flexible body systems | |
| CN106503404B (en) | A kind of PBX blasting explosive granules particle Random Discrete fill method | |
| Lin et al. | Size effect and material property effect of the impactor on the damage modes of the single-layer Kiewitt-8 reticulated dome | |
| Wang et al. | Explosion resistance of three-dimensional mesoscopic model of complex closed-cell aluminum foam sandwich structure based on random generation algorithm | |
| Ma et al. | A hybrid VOF and PIC multi-material interface treatment method and its application in the penetration | |
| CN111125917A (en) | Construction method of single-particle discrete element model containing inner pores | |
| Wang et al. | Study on the effects of tamping frequency to the compaction degree of trackbeds | |
| Miedzińska | New method of numerical homogenization of functionally graded materials | |
| Huang et al. | Modelling and Simulation on Penetration into Stone Concrete | |
| Stoimenov | Advanced computing for energy efficiency of milling processes | |
| Vidal et al. | Analysis of strength of brittle materials under different strain rates using ldem simulations | |
| Ma | Influence of the Second Phase Particle’s Size Distribution on Grain Growth Based on Computer Simulation | |
| Wiśniewski et al. | Numerical Simulations of Penetration of 9 mm Parabellum Bullet into Kevlar Layers: Erosion Selection in Autodyn Program | |
| Orozco et al. | Discrete-element simulation of particle breakage inside ball mills: A 2D model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |