CN107657128A - The coarse aggregate granular discrete-element method of broken state - Google Patents

The coarse aggregate granular discrete-element method of broken state Download PDF

Info

Publication number
CN107657128A
CN107657128A CN201710946130.8A CN201710946130A CN107657128A CN 107657128 A CN107657128 A CN 107657128A CN 201710946130 A CN201710946130 A CN 201710946130A CN 107657128 A CN107657128 A CN 107657128A
Authority
CN
China
Prior art keywords
mrow
particle
coarse aggregate
msubsup
particle size
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
Application number
CN201710946130.8A
Other languages
Chinese (zh)
Other versions
CN107657128B (en
Inventor
马涛
丁珣昊
黄晓明
张玉贞
郭兆民
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201710946130.8A priority Critical patent/CN107657128B/en
Publication of CN107657128A publication Critical patent/CN107657128A/en
Application granted granted Critical
Publication of CN107657128B publication Critical patent/CN107657128B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a kind of coarse aggregate granular discrete-element method of broken state, comprise the following steps:(1) coarse aggregate particle primary dcreening operation;(2) CT scan;(3) discrete meta software PFC2D modelings;(4) single particle size coarse aggregate particle contact stiffness model parameter calibration;(5) crushing counts under single particle size coarse aggregate particle penetration resistance;(6) crushing counts under the coarse aggregate particle penetration resistance of different gradation;(7) it is based on the parallel link model parameter of coarse aggregate Particle Breakage amount amendment.The present invention breaches the difficult point that virtually coarse aggregate particle can not crush during road engineering discrete element simulation, using Cluster models in discrete meta software, the crushing statistics tested with reference to indoor penetration resistance, the micromodel parameter of different-grain diameter particle is determined to simulate granular materials shattering process under outside drag effect, the accuracy that virtual coarse aggregate particle shape reconstructs is ensured by scan image.

Description

The coarse aggregate granular discrete-element method of broken state
Technical field
The present invention relates to aggregate particle Discrete-parcel method in road engineering, more particularly to a kind of coarse aggregate of broken state Grain Discrete-parcel method.
Background technology
In road engineering numerical discretization simulated test, the broken state simulation of particle is always difficult point, turns into accurate mould Intend the obstruction of emulation, it is most of to use discrete meta software PFC2D/3D in existing analogy method, in simulation coarse aggregate particle When, clump command definitions are rigid body, and fictitious particle can not crush under the order, when realistic simulation is tested, this method into For the fictitious particle of type under big stress, big strained condition, virtual emulation result and actual difference are larger, can not portray particle and break The time of day that broken, test specimen destroys.
The content of the invention
Goal of the invention:The present invention provides a kind of coarse aggregate granular discrete-element method of broken state, breaches road engineering The difficult point that virtual coarse aggregate particle can not crush during discrete element simulation.
Technical scheme:The present invention comprises the following steps:
(1) coarse aggregate particle primary dcreening operation:Coarse aggregate particle is sieved by standard screen, it is thick that according to particle size it is divided into multi gear Aggregate particle, each grade of coarse aggregate particle screened out ensure that quantity is enough;
(2) CT scan:CT scan is carried out to each grade of coarse aggregate particle screened out, obtains each coarse aggregate particle Form contour images;
(3) discrete meta software PFC2D modelings:Coarse aggregate particle outline information is identified with reference to image processing techniques, discrete Established in meta software series of disks race reconstruct fictitious particle form, and using Cluster orders by disk race be defined as one it is whole Body, while establish the virtual test specimen model of uniaxial penetration test;
(4) single particle size coarse aggregate particle contact stiffness model parameter calibration:The coarse aggregate single shaft for carrying out different-grain diameter passes through Enter laboratory test and virtual test, determine the contact stiffness model micro-parameter of each grade of coarse aggregate particle;
(5) crushing counts under single particle size coarse aggregate particle penetration resistance:Determine under uniaxial penetration laboratory test, it is each Shelves coarse aggregate particle broken granular mass fraction because of outside penetration resistance;
(6) crushing counts under the coarse aggregate particle penetration resistance of different gradation:Determine under uniaxial penetration laboratory test, select The each grade of coarse aggregate particle that deciding grade and level is matched somebody with somebody broken granular mass fraction because of outside penetration resistance;
(7) it is based on the parallel link model parameter of coarse aggregate Particle Breakage amount amendment:Assign the fictitious particle that Cluster is defined Parallel link model, using ensure fictitious particle as it is deformable, can crush rigid body, according to single particle size, the coarse aggregate of different gradation Particle Breakage amount under particle penetration resistance, complete the parameters revision of parallel link model.
In step (1), after coarse aggregate screening, using cylindrical steel containers hold coarse aggregate particle., will in step (3) Particle outline information after CT scan and processing is directed into discrete meta software PFC2D, each coarse aggregate particle of precise forming True form, each coarse aggregate fictitious particle is made up of a series of disk races, is different from traditional Clump command methods, Fictitious particle is defined with Cluster command streams, the contact attribute of disk is defined with contact stiffness model, is determined with parallel link model The initial connection attribute of disk in adopted Cluster disks race.In step (4), different shelves single particle size coarse aggregate particles are carried out In uniaxial penetration room, virtual experimental, fitting virtual experimental data make the indoor experimental data of its unlimited approaching to reality, obtain simultaneously true The contact stiffness model micro-parameter of fixed each grade of single particle size particle, identified contact stiffness model micro-parameter foundation Grain particle diameter division, i.e. the particle contact stiffness parameter of different-grain diameter are not consistent, it is necessary to be separately determined.In step (5), to difference Particle Breakage amount after shelves single particle size coarse aggregate uniaxial penetration test is counted, the coarse aggregate uniaxial penetration of each grade of particle diameter Experiment progress is multigroup, and the coarse aggregate particle after every group of experiment is counted again and sieved successively, determines each grade of coarse aggregate The particle new composition of particle size, formula specific as follows after penetration resistance is broken:
In formula:I is the coarse aggregate particle ith uniaxial penetration test of this grade of particle diameter, i=1,2,3;J is that jth shelves gather materials Grain, wherein j=2,3,4,5,6,7;Gather materials prothyl at the beginning of ith uniaxial penetration test particle for jth shelves single particle size particle Amount;Gather materials the screen over-size after ith uniaxial penetration test on jth shelves standard screen for jth shelves single particle size particle;Respectively jth shelves single particle size particle gathers materials ith uniaxial penetration After experiment, jth -1, the screen over-size of 1 grade of standard screen of j-2, j-3 ...;ratiojGather materials for jth shelves single particle size particle by list The mass ratio of unbroken particle after axle penetration test;Δratioj-j-1、Δratioj-j-2、Δratioj-j-3........... Δratioj-1Respectively jth shelves single particle size particle gathers materials crushed after uniaxial penetration test after to be changed into particle size be j- 1st, the mass ratio of 1 grade of particles of j-2 ....Secondly, in step (5), during to single particle size Particle Breakage amount statistical analysis, Determine Particle Breakage amount of each grade of particle diameter under single particle size penetration test.In step (6), carried out using different gradings single Axle penetration test, and count the crushing of different-grain diameter particle in same test specimen;According to the different gradation of selection, in forming room Uniaxial penetration test test specimen, the color of each grade of particles need to keep different, each shelves particles determined according to gradation composition Institute's expense, cylinder specimen by dry-tamping compacting aftershaping different gradation indoor penetration test test specimen.Secondly, step (6) in, the indoor test specimen matched somebody with somebody to selected class and be molded carries out uniaxial penetration test, wherein only with the coarse aggregate shape in grading Into skeleton without fine aggregate, after the completion of uniaxial penetration test, the level of coarse aggregate particle in test specimen is counted again using standard screen With composition;The size of the broken particle initial particle is distinguished by the grain color remained on standard screen after experiment, to distinguish The broken state of different-grain diameter particle.In step (7), according to single particle size Particle Breakage amount, different gradation Particle Breakage amount, really Determine the micro-parameter of the parallel link model of Cluster fictitious particles established in PFC2D.
Beneficial effect:The advantage of the invention is that:(1) the two Dimension Numerical Value particle shape height one based on CT scan shaping Cause;(2) defined by cluster command streams and reconstruct fictitious particle model, assigned parallel link model parameter, realize coarse aggregate The broken state simulation of particle;(3) united according to the uniaxial penetration test Particle Breakage amount of different shelves single particle sizes, different gradation Meter, propose it is a kind of can crushed particles parallel link model determination method for parameter, important finger is provided for virtual simulation test Lead.
Brief description of the drawings
Fig. 1 is the flow chart of the present invention;
Fig. 2 is the virtual test model schematic of the present invention.
Embodiment
The present invention is further described below in conjunction with the accompanying drawings.
As shown in figure 1, the present invention comprises the following steps:
(1) coarse aggregate particle primary dcreening operation:Coarse aggregate particle is sieved by standard screen, it is divided into 4.75- according to particle size 9.5mm, 9.5-13.2mm, 13.2-16mm, 16-19mm, 19-26.5mm, 26.5-31.5mm, each grade of coarse aggregate screened out Particle should ensure that quantity is enough;After coarse aggregate screening, using cylindrical steel containers hold coarse aggregate particle, cylindrical steel holds Device a width of 12cm, a height of 15cm, the screening quantity of each grade of coarse aggregate particle it is ensured that fill the cylinder of 3 selected sizes enough Describe device.
(2) CT scan:CT scan is carried out to each grade of coarse aggregate particle screened out, obtains each coarse aggregate particle Form contour images.
(3) discrete meta software PFC2D modelings:As shown in Fig. 2 with reference to image processing techniques identification coarse aggregate particle outline letter Breath, the particle outline information after CT scan and processing is directed into discrete meta software PFC2D;In the built-in vertical system of discrete meta software Row disk race reconstruct fictitious particle form, the true form of each coarse aggregate particle of precise forming, each coarse aggregate are virtual Particle is made up of a series of disk races;Traditional Clump command methods are different from, are defined disk race using Cluster orders For an entirety;The virtual test specimen model of uniaxial penetration test is established simultaneously, the contact attribute of disk is defined with contact stiffness model, The initial connection attribute of disk in Cluster disks race is defined with parallel link model.
(4) single particle size coarse aggregate particle contact stiffness model parameter calibration:To different shelves single particle size coarse aggregate particles Carry out in uniaxial penetration room, virtual experimental, fitting virtual experimental data make the indoor experimental data of its unlimited approaching to reality, obtain And determine the contact stiffness model micro-parameter of each grade of single particle size particle, identified contact stiffness model micro-parameter according to Divided according to grain diameter, i.e. the particle contact stiffness parameter of different-grain diameter is not consistent, it is necessary to be separately determined.
(5) crushing counts under single particle size coarse aggregate particle penetration resistance:Determine under uniaxial penetration laboratory test, it is each Shelves coarse aggregate particle broken granular mass fraction because of outside penetration resistance;To different shelves single particle size coarse aggregate uniaxial penetrations Particle Breakage amount after experiment is counted, and the coarse aggregate uniaxial penetration test of each grade of particle diameter carries out 3 groups, and every group is tested Coarse aggregate particle afterwards is counted again to be sieved successively, determines each grade of coarse aggregate particle particle diameter after penetration resistance is broken The new composition of size, by taking crushing statistics after 26.5-31.5mm single particle size particle penetration tests as an example, such as following formula:
In formula:I is the coarse aggregate particle ith uniaxial penetration test of this grade of particle diameter, i=1,2,3;For 26.5- 31.5mm single particle size particle ith uniaxial penetration test particle initial mass;For 26.5-31.5mm single grain Screen over-size after the particle ith uniaxial penetration test of footpath on 26.5mm standard screens; Respectively 26.5-31.5mm single particle size particle ith single shaft After penetration test, the screen over-size of 19mm, 16mm, 13.2mm, 9.5mm, 4.75mm standard screen;For 26.5- After 31.5mm single particle size particle ith uniaxial penetration test, the throughput of 4.75mm standard screens;ratio26.5For 26.5- The mass ratio of 31.5mm single particle size coarse aggregate particle unbroken particle after uniaxial penetration test;ratio26.5-19、 ratio26.5-16、ratio26.5-13.2ratio26.5-9.5、ratio26.5-4.75、ratio26.5-fineRespectively 26.5-31.5mm list One particle diameter coarse aggregate particle is changed into particle size 19-26.5mm, 16-19mm, 13.2- after being crushed after uniaxial penetration test 16mm, 9.5-13.2mm, 4.75-9.5mm and the mass ratio less than 4.75mm particles.
According to above single particle size Particle Breakage amount statistical analysis, 4.75-9.5,9.5-13.2,13.2-16,16-19, 19-26.5mm single particle size particle according to above method, determines of each grade of particle diameter under single particle size penetration test Grain crushing.4.75-9.5mm particle is provided with parameter ratio4.75、ratio4.75-fine, i.e. 4.75-9.5mm particles are unbroken, broken The broken mass ratio for being less than 4.75mm for particle diameter;9.5-13.2mm it is provided with parameter ratio9.5、ratio9.5-4.75、ratio9.5 -fine, i.e., 9.5-13.2mm particles are unbroken, are broken for particle diameter as 4.75-9.5mm, the mass ratio less than 4.75mm;13.2- 16mm is provided with parameter ratio13.2、ratio13.2-9.5、ratio13.2-4.75、ratio13.2-fine, i.e. 13.2-16mm particles do not break It is broken, be broken for particle diameter for 9.5-13.2mm, 4.75-9.5mm, the mass ratio less than 4.75mm;16-19mm is provided with parameter ratio16、ratio16-13.2ratio16-9.5ratio16-4.75、ratio16-fine, i.e., 16-19mm particles are unbroken, are broken for grain Footpath is 13.2-16mm, 9.5-13.2mm, 4.75-9.5mm, the mass ratio less than 4.75mm;19-26.5mm it is provided with parameter ratio19、ratio19-16、ratio19-13.2、ratio19-9.5、ratio19-4.75、ratio19-fine, i.e., 19-26.5mm particles are not Crush, be broken for particle diameter as 16-19mm, 13.2-16mm, 9.5-13.2mm, 4.75-9.5mm, the mass ratio less than 4.75mm Example.
(6) crushing counts under the coarse aggregate particle penetration resistance of different gradation:Determine under uniaxial penetration laboratory test, select The each grade of coarse aggregate particle that deciding grade and level is matched somebody with somebody broken granular mass fraction because of outside penetration resistance;Carried out using different gradings Uniaxial penetration test, and the crushing of different-grain diameter particle in same test specimen is counted, wherein grading selection is arbitrarily not intended to limit, Principle is consistent.According to the different gradation of selection, uniaxial penetration test test specimen in forming room, before shaping test piece, by what is screened out Each grade of single particle size particle, which is cleaned, keeps particle surface without obvious miscellaneous dirt, with different dyestuffs respectively to each grade of particles Dyed, and dried, the color of each grade of particles need to keep different, and each shelves particles are determined according to gradation composition Institute's expense, using indoor penetration test test specimen of the cylinder specimen by dry-tamping compacting aftershaping different gradation.
When the indoor test specimen matched somebody with somebody to selected class and be molded carries out uniaxial penetration test, wherein only with the rough set in grading Material forms skeleton without fine aggregate, after the completion of uniaxial penetration test, counts coarse aggregate particle in test specimen again using standard screen Gradation composition, due to being made up of in test specimen the coarse aggregate particle of various particle diameters, so by being remained on standard screen after experiment Grain color distinguishes the size of the broken particle initial particle, to distinguish the broken state of different-grain diameter particle.With grading Exemplified by SMA13, when carrying out grading crushing statistical analysis, carried out according to following principle:
MM13.2=mm13.2+Δmm13.2-9.5+Δmm13.2-4.75+Δmm13.2-fine
MM9.5=mm9.5+Δmm9.5-4.75+Δmm9.5-fine
MM4.75=mm4.75+Δmm4.75-fine
mmfine=Δ mm13.2-fine+Δmm9.5-fine+Δmm4.75-fine
In formula:MM13.2、MM9.5、MM4.75Respectively before SMA13 test specimens uniaxial penetration test, 13.2-16,9.5- in test specimen 13.2nd, the initial mass of 4.75-9.5mm particles;mmfineAfter uniaxial penetration test, it is broken for being less than 4.75mm particle diameters The quality of particle;mm13.2、mm9.5、mm4.75Respectively after uniaxial penetration test, 13.2-16,9.5-13.2,4.75-9.5mm grain A unbroken quality in footpath;Δmm13.2-9.5、Δmm13.2-4.75、Δmm13.2-fineRespectively 13.2-16mm particles pass through Particle size is broken for after uniaxial test as 9.5-13.2,4.75-9.5mm and the granular mass less than 4.75mm;Δ mm9.5-4.75、Δmm9.5-fineIt is 4.75- that respectively 9.5-13.2mm particles are broken for particle size after uniaxial test 9.5mm and the granular mass less than 4.75mm;Δmm4.75-fineCrushed for 4.75-9.5mm particles after uniaxial test It is less than 4.75mm granular mass for particle size.
(7) it is based on the parallel link model parameter of coarse aggregate Particle Breakage amount amendment:Assign the fictitious particle that Cluster is defined Parallel link model, using ensure fictitious particle as it is deformable, can crush rigid body, according to single particle size, the coarse aggregate of different gradation Particle Breakage amount under particle penetration resistance, complete the parameters revision of parallel link model.According to single particle size Particle Breakage amount, Different gradation Particle Breakage amount, determine the micro-parameter of the parallel link model of Cluster fictitious particles established in PFC2D.It is first The coarse aggregate particle virtual experimental of different single particle sizes is first carried out, assigns random initial parallel link model micro-parameter, it is empty After intending off-test, the fragmentation quality of each grade of cluster particle in virtual experimental is read using the code language write, to not Parallel link model with shelves single particle size cluster particles carries out micro-parameter adjustment so that is obtained in virtual test ratio26.5-19、ratio26.5-16、ratio26.5-13.2、ratio26.5-9.5、ratio26.5-4.75、ratio26.5-fine、ratio19、 ratio19-16、ratio19-13.2、ratio19-9.5、ratio19-4.75、ratio19-fine、ratio16、ratio16-13.2、 ratio16-9.5、ratio16-4.75、ratio16-fine、ratio13.2、ratio13.2-9.5、ratio13.2-4.75、ratio13.2-fine、 ratio9.5、ratio9.5-4.75、ratio9.5-fine、ratio4.75、ratio4.75-fineIt is consistent substantially with laboratory test, error Control temporarily determines the optimal parallel link model micro-parameter of each grade of particles within 5%, by such method; Secondly the virtual simulation test matched somebody with somebody of selected class is carried out, using optimal parallel link model parameter and the contact model parameter of acquisition, Virtual test keeps material, experimental enviroment consistent with indoor grading uniaxial penetration test, reads each under the grading in virtual test Shelves particles fragmentation quality fraction, and contrasted with laboratory test results, error must be controlled within 5%.If virtual test knot Fruit error then further determines that within 5% and verifies the correctness of acquired parallel link model micro-parameter, if error More than 5%, then need to readjust the parallel link model micro-parameter of single particle size coarse aggregate particle, until virtual test knot The indoor penetration test that fruit had both met single particle size crushes situation, also complies with the indoor penetration test that selected class is matched somebody with somebody and crushes situation.
In this method in all uniaxial penetration rooms, virtual test loading velocity use 1.27mm/min, loading head is straight Footpath is 50mm, load time 1h;There is material category and the difference of different batches in coarse aggregate particle, can be according to institute of the present invention The method of offer, can crush modeling and the micromodel parameter setting of fictitious particle.

Claims (9)

1. a kind of coarse aggregate granular discrete-element method of broken state, it is characterised in that:Comprise the following steps:
(1) coarse aggregate particle primary dcreening operation:Coarse aggregate particle is sieved by standard screen, it is divided into multi gear coarse aggregate according to particle size Particle, each grade of coarse aggregate particle screened out ensure that quantity is enough;
(2) CT scan:CT scan is carried out to each grade of coarse aggregate particle screened out, obtains the form of each coarse aggregate particle Contour images;
(3) discrete meta software PFC2D modelings:Coarse aggregate particle outline information is identified with reference to image processing techniques, it is soft in discrete element Series of disks race reconstruct fictitious particle form is established in part, and disk race is defined as an entirety using Cluster orders, together The virtual test specimen model of Shi Jianli uniaxial penetration tests;
(4) single particle size coarse aggregate particle contact stiffness model parameter calibration:Carry out the coarse aggregate uniaxial penetration room of different-grain diameter Interior experiment and virtual test, determine the contact stiffness model micro-parameter of each grade of coarse aggregate particle;
(5) crushing counts under single particle size coarse aggregate particle penetration resistance:Determine under uniaxial penetration laboratory test, each shelves are thick Aggregate particle broken granular mass fraction because of outside penetration resistance;
(6) crushing counts under the coarse aggregate particle penetration resistance of different gradation:Determine under uniaxial penetration laboratory test, selected class The each grade of coarse aggregate particle matched somebody with somebody broken granular mass fraction because of outside penetration resistance;
(7) it is based on the parallel link model parameter of coarse aggregate Particle Breakage amount amendment:It is parallel to assign the fictitious particle that Cluster is defined Link model, using ensure fictitious particle as it is deformable, can crush rigid body, according to single particle size, the coarse aggregate particle of different gradation Particle Breakage amount under penetration resistance, complete the parameters revision of parallel link model.
2. the coarse aggregate granular discrete-element method of broken state according to claim 1, it is characterised in that:The step (1) in, after coarse aggregate screening, using cylindrical steel containers hold coarse aggregate particle.
3. the coarse aggregate granular discrete-element method of broken state according to claim 1, it is characterised in that:The step (3) in, the particle outline information after CT scan and processing is directed into discrete meta software PFC2D, each is thick for precise forming The true form of aggregate particle, each coarse aggregate fictitious particle are made up of a series of disk races, are different from traditional Clump lives Method is made, fictitious particle is defined with Cluster command streams, the contact attribute of disk is defined with contact stiffness model, with parallel company Connect the initial connection attribute of disk in model definition Cluster disks race.
4. the coarse aggregate granular discrete-element method of broken state according to claim 1, it is characterised in that:The step (4) in, different shelves single particle size coarse aggregate particles are carried out in uniaxial penetration room, virtual experimental, fitting virtual experimental data make The indoor experimental data of its unlimited approaching to reality, obtain and determine the microcosmic ginseng of contact stiffness model of each grade of single particle size particle Number, identified contact stiffness model micro-parameter divide according to grain diameter, i.e. the particle contact stiffness parameter of different-grain diameter It is not consistent, it is necessary to be separately determined.
5. the coarse aggregate granular discrete-element method of broken state according to claim 1, it is characterised in that:The step (5) in, the Particle Breakage amount after different shelves single particle size coarse aggregate uniaxial penetration tests is counted, each grade of particle diameter it is thick The uniaxial penetration test that gathers materials progress is multigroup, and the coarse aggregate particle after every group of experiment is counted again and sieved successively, it is determined that Each grade of coarse aggregate the particle new composition of particle size, formula specific as follows after penetration resistance is broken:
<mrow> <msubsup> <mi>M</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>=</mo> <msubsup> <mi>m</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>+</mo> <msubsup> <mi>&amp;Delta;m</mi> <mrow> <mi>j</mi> <mo>~</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> <mo>+</mo> <msubsup> <mi>&amp;Delta;m</mi> <mrow> <mi>j</mi> <mo>~</mo> <mi>j</mi> <mo>-</mo> <mn>2</mn> </mrow> <mi>i</mi> </msubsup> <mo>+</mo> <msubsup> <mi>&amp;Delta;m</mi> <mrow> <mi>j</mi> <mo>~</mo> <mi>j</mi> <mo>-</mo> <mn>3</mn> </mrow> <mi>i</mi> </msubsup> <mo>+</mo> <mo>...</mo> <mo>...</mo> <mo>...</mo> <mo>...</mo> <mo>...</mo> <mo>+</mo> <msubsup> <mi>&amp;Delta;m</mi> <mrow> <mi>j</mi> <mo>~</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> </mrow>
<mrow> <msub> <mi>ratio</mi> <mi>j</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mi>m</mi> <mi>j</mi> <mi>i</mi> </msubsup> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mi>M</mi> <mi>j</mi> <mi>i</mi> </msubsup> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>ratio</mi> <mrow> <mi>j</mi> <mo>~</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mi>&amp;Delta;m</mi> <mrow> <mi>j</mi> <mo>~</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mi>&amp;Delta;M</mi> <mrow> <mi>j</mi> <mo>~</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>ratio</mi> <mrow> <mi>j</mi> <mo>~</mo> <mi>j</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mi>&amp;Delta;m</mi> <mrow> <mi>j</mi> <mo>~</mo> <mi>j</mi> <mo>-</mo> <mn>2</mn> </mrow> <mi>i</mi> </msubsup> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mi>&amp;Delta;M</mi> <mrow> <mi>j</mi> <mo>~</mo> <mi>j</mi> <mo>-</mo> <mn>2</mn> </mrow> <mi>i</mi> </msubsup> </mrow> </mfrac> </mrow>
…………………………
<mrow> <msub> <mi>ratio</mi> <mrow> <mi>j</mi> <mo>~</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mi>&amp;Delta;m</mi> <mrow> <mi>j</mi> <mo>~</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mi>&amp;Delta;M</mi> <mrow> <mi>j</mi> <mo>~</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> </mrow> </mfrac> </mrow>
In formula:I is the coarse aggregate particle ith uniaxial penetration test of this grade of particle diameter, i=1,2,3;J is jth shelves aggregate particle, Wherein j=2,3,4,5,6,7;Gather materials ith uniaxial penetration test particle initial mass for jth shelves single particle size particle; Gather materials the screen over-size after ith uniaxial penetration test on jth shelves standard screen for jth shelves single particle size particle;Respectively jth shelves single particle size particle gathers materials ith uniaxial penetration After experiment, jth -1, the screen over-size of 1 grade of standard screen of j-2, j-3 ...;ratiojGather materials for jth shelves single particle size particle by list The mass ratio of unbroken particle after axle penetration test;
Δratioj-j-1、Δratioj-j-2、Δratioj-j-3...........Δratioj-1Respectively jth shelves single particle size Particle gather materials after uniaxial penetration test crush after be changed into particle size be 1 grade of particles of j-1, j-2 ... mass ratio Example.
6. the coarse aggregate granular discrete-element method of broken state according to claim 5, it is characterised in that:The step (5) in, during to single particle size Particle Breakage amount statistical analysis, particle of each grade of particle diameter under single particle size penetration test is determined Crushing.
7. the coarse aggregate granular discrete-element method of broken state according to claim 1, it is characterised in that:The step (6) in, uniaxial penetration test is carried out using different gradings, and count the crushing of different-grain diameter particle in same test specimen; According to the different gradation of selection, uniaxial penetration test test specimen in forming room, the color of each grade of particles need to keep different, root Institute's expense of each shelves particles, room of the cylinder specimen by dry-tamping compacting aftershaping different gradation are determined according to gradation composition Interior penetration test test specimen.
8. the coarse aggregate granular discrete-element method of broken state according to claim 1, it is characterised in that:The step (6) in, the indoor test specimen matched somebody with somebody to selected class and be molded carries out uniaxial penetration test, wherein only with the coarse aggregate shape in grading Into skeleton without fine aggregate, after the completion of uniaxial penetration test, the level of coarse aggregate particle in test specimen is counted again using standard screen With composition;The size of the broken particle initial particle is distinguished by the grain color remained on standard screen after experiment, to distinguish The broken state of different-grain diameter particle.
9. the coarse aggregate granular discrete-element method of broken state according to claim 1, it is characterised in that:The step (7) in, according to single particle size Particle Breakage amount, different gradation Particle Breakage amount, determine that the Cluster established in PFC2D is virtual The micro-parameter of particles parallel link model.
CN201710946130.8A 2017-10-12 2017-10-12 Coarse aggregate particle discrete element method in crushing state Active CN107657128B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710946130.8A CN107657128B (en) 2017-10-12 2017-10-12 Coarse aggregate particle discrete element method in crushing state

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710946130.8A CN107657128B (en) 2017-10-12 2017-10-12 Coarse aggregate particle discrete element method in crushing state

Publications (2)

Publication Number Publication Date
CN107657128A true CN107657128A (en) 2018-02-02
CN107657128B CN107657128B (en) 2020-12-25

Family

ID=61117871

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710946130.8A Active CN107657128B (en) 2017-10-12 2017-10-12 Coarse aggregate particle discrete element method in crushing state

Country Status (1)

Country Link
CN (1) CN107657128B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108593428A (en) * 2018-04-24 2018-09-28 重庆大学 A kind of discrete material Particle Breakage test device and test method
CN108918356A (en) * 2018-05-18 2018-11-30 中建材(合肥)粉体科技装备有限公司 A kind of mixed material squeezes smashed granularity prediction technique and system
CN109191571A (en) * 2018-09-30 2019-01-11 华南理工大学 A method of gather materials using 3D printing technique preparation mechanical test standard
CN110487240A (en) * 2019-08-29 2019-11-22 哈尔滨工业大学 Surface deformation monitoring induction element and its application based on 3D printing technique
CN112131633A (en) * 2020-09-04 2020-12-25 山东大学 Fluid-solid coupling simulation method and system based on coarse graining calculation theory
CN114166631A (en) * 2021-11-30 2022-03-11 温州理工学院 Test method for identifying flow and crushing mechanism of discrete material under external load

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103163033A (en) * 2013-03-29 2013-06-19 东南大学 Numerical value measuring method for resilience modulus of graded broken stones
CN104615482A (en) * 2015-02-10 2015-05-13 长安大学 Simulation method of breakage of cold recycling mixture in compaction process
CN106096135A (en) * 2016-06-12 2016-11-09 东南大学 A kind of grating granule library file numerical generation method based on AIMS equipment
CN106950120A (en) * 2017-02-16 2017-07-14 浙江大学 A kind of asphalt mixture gradation optimization method under virtual state

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103163033A (en) * 2013-03-29 2013-06-19 东南大学 Numerical value measuring method for resilience modulus of graded broken stones
CN104615482A (en) * 2015-02-10 2015-05-13 长安大学 Simulation method of breakage of cold recycling mixture in compaction process
CN106096135A (en) * 2016-06-12 2016-11-09 东南大学 A kind of grating granule library file numerical generation method based on AIMS equipment
CN106950120A (en) * 2017-02-16 2017-07-14 浙江大学 A kind of asphalt mixture gradation optimization method under virtual state

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
TAO MA等: "Influences by air voids on fatigue life of asphalt mixture based on discrete element method", 《CONSTRUCTION AND BUILDING MATERIALS》 *
张程林 等: "砂土颗粒三维形状模拟离散元算法研究", 《岩土工程学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108593428A (en) * 2018-04-24 2018-09-28 重庆大学 A kind of discrete material Particle Breakage test device and test method
CN108918356A (en) * 2018-05-18 2018-11-30 中建材(合肥)粉体科技装备有限公司 A kind of mixed material squeezes smashed granularity prediction technique and system
CN108918356B (en) * 2018-05-18 2021-07-09 中建材(合肥)粉体科技装备有限公司 Method and system for predicting particle size of mixed material after extrusion and crushing
CN109191571A (en) * 2018-09-30 2019-01-11 华南理工大学 A method of gather materials using 3D printing technique preparation mechanical test standard
CN109191571B (en) * 2018-09-30 2023-09-12 华南理工大学 Method for preparing mechanical test standard aggregate by applying 3D printing technology
CN110487240A (en) * 2019-08-29 2019-11-22 哈尔滨工业大学 Surface deformation monitoring induction element and its application based on 3D printing technique
CN112131633A (en) * 2020-09-04 2020-12-25 山东大学 Fluid-solid coupling simulation method and system based on coarse graining calculation theory
CN112131633B (en) * 2020-09-04 2023-01-13 山东大学 Fluid-solid coupling simulation method and system based on coarse graining calculation theory
CN114166631A (en) * 2021-11-30 2022-03-11 温州理工学院 Test method for identifying flow and crushing mechanism of discrete material under external load

Also Published As

Publication number Publication date
CN107657128B (en) 2020-12-25

Similar Documents

Publication Publication Date Title
CN107657128A (en) The coarse aggregate granular discrete-element method of broken state
Shi et al. Calibration of micro-scaled mechanical parameters of granite based on a bonded-particle model with 2D particle flow code
CN110619393A (en) Traffic simulation software parameter calibration method, system and medium based on learning algorithm
CN104200217B (en) Hyperspectrum classification method based on composite kernel function
CN111523616B (en) Coarse-grained soil filler grading identification method based on convolutional neural network and application system
CN106651896A (en) Material classification method of CT images of asphalt mixture
JP6689629B2 (en) Permeability evaluation method of ground material using a test piece that reproduces ground material
CN101675453A (en) Be used to analyze the graphics command management tool and the method for the performance that changes for order before the application modification
CN108171175B (en) Deep learning sample enhancement system and operation method thereof
DE10128160A1 (en) Method and device for extracting measurement information and for adjusting specifications through 3D visualization translates human visual test analysis of an object into physically and quantitatively determinable parameters.
Groh et al. Damage simulation of brittle heterogeneous materials at the grain size level
Khosrowshahi Simulation of expenditure patterns of construction projects
CN106295351B (en) A kind of Risk Identification Method and device
DE102004063769A1 (en) Method for automatically and quantitatively determining the amount of seed or grain of required quality comprises recording the seed and grain using an imaging device and further processing
US7801705B2 (en) Mass on model
CN109754392A (en) A kind of device and method that defect image automatically analyzes
CN105354600B (en) A kind of automatic classification method of sandstone microsection
CN108875788A (en) A kind of SVM classifier parameter optimization method based on modified particle swarm optiziation
CN113176153A (en) Asphalt mixture shear strength analysis method based on discrete elements
CN112381644A (en) Credit scene risk user assessment method based on space variable reasoning
CN115730501A (en) Soil body elastic modulus prediction method based on three-dimensional discrete element model and convolutional neural network
CN107480634A (en) A kind of geographical national conditions ground mulching monitoring method based on multistage target classification
CN114186434A (en) Asphalt mixture three-dimensional microscopic structure model construction method based on step-by-step division
CN116150667A (en) XGBoost-based flight delay reason pre-judging method
DE112020003599T5 (en) TEST DEVICE, TEST PROCEDURE, TEST PROGRAM, TEACHING DEVICE, TRAINING METHOD AND TEACHING PROGRAM

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