1. A method for analyzing the pore connectivity of a cement-based material particle stacking system is characterized by comprising the following steps: the method for analyzing the pore connectivity of the sphere accumulation system comprises the steps of firstly establishing a sphere particle accumulation model, secondly executing Delaunay tetrahedron subdivision on a sphere center point set to divide the accumulation model into a set consisting of a plurality of tetrahedrons, then acquiring the pore volume of the tetrahedrons in the set one by one, judging the connectivity of each tetrahedron pore, and finally performing a statistical calculation process of the connected porosity in the sphere particle accumulation model.

2. The method for analyzing the pore connectivity of the cement-based material particle packing system according to claim 1, wherein the method comprises the following steps: the process of establishing the sphere particle accumulation model comprises the following steps:

for a system with specified spherical particle size distribution and porosity, firstly determining the size of a cubic filling space, then sequentially putting spherical particles into the cubic filling space from large to small according to the size, wherein the coordinate of each spherical particle is randomly generated by a Monte Carlo method until all the spherical particles are put in, and forming a spherical particle stacking model.

3. The method for analyzing the pore connectivity of the cement-based material particle packing system according to claim 2, wherein the method comprises the following steps: based on a sphere particle accumulation model, the operation process of performing Delaunay tetrahedron subdivision on the sphere center point set is as follows:

according to the information of the sphere particle accumulation model, coordinates of centers of all spheres are arranged to form a two-dimensional matrix, the two-dimensional matrix has three rows which are respectively an abscissa, an ordinate and a column coordinate of the sphere center, each row in the matrix represents a sphere coordinate of one sphere, and the two-dimensional matrix is acted by a Delaunay Triangulation function in MATLAB, so that Delaunay tetrahedron subdivision of a sphere center point set is realized; in the two-dimensional case, the matrix has only two columns, respectively set of abscissas and sets of column coordinates of the circle, on which the Delaunay Triangulation function is applied to perform a Delaunay Triangulation of a set of circle center points.

4. The method for analyzing the pore connectivity of the cement-based material particle packing system according to claim 3, wherein the method comprises the following steps: after the Delaunay tetrahedron subdivision is performed on the sphere center point set, the stacking model is divided into a set consisting of a plurality of tetrahedrons, and the process of determining the pore volumes of the tetrahedrons in the set one by one is as follows:

the particle stacking system after being subdivided by the Delaunay tetrahedron is composed of a plurality of mutually independent and non-overlapping tetrahedrons, wherein each tetrahedron is composed of a solid phase occupied by four vertex angle spherical particles and a pore phase, the solid phase is the sum of sector areas intersected by the three corner point particles on the surface, the pore phase is the residual phase of the solid phase occupied by the four vertex angle spherical particles subtracted from the tetrahedron, the pore phase of the tetrahedron is the pore of the tetrahedron, and the corresponding volume of the pore is the pore volume of the tetrahedron;

for a sphere stacking system allowing overlapping existence, the following scheme is adopted for calculation due to the complex structure:

randomly putting N particles into a tetrahedron, generating the particle coordinates by a Monte Carlo method, counting the number N of the particles falling in the pores, and determining the volume V of the pores in the tetrahedron_pCan be calculated by the formula one:

V_p＝(n/N)*V_t (1)。

5. the method for analyzing the pore connectivity of the particle packing system of the cement-based material as claimed in claim 1, 2, 3 or 4, wherein: the process of judging the connectivity of each tetrahedral pore is to firstly judge the connectivity of all the tetrahedral pores after traversing and judging the connectivity of the adjacent tetrahedral pores;

the tetrahedron is provided with four triangular surfaces, each triangular surface consists of a solid phase and a pore phase, wherein the surface solid phase is the sum of sector areas intersected by particles with three angular points on the surface, and the surface pore phase is the residual phase of the tetrahedron surface minus the corresponding solid phase; the surface pore phase is the cross section of the narrowest part of a channel connecting the tetrahedral pore and the adjacent tetrahedral pore, namely the throat cross section;

when the throat cross section between any two adjacent tetrahedral voids is completely filled with spherical particles, the connection relationship between the two adjacent tetrahedral voids is shown to be a mutual blocking relationship, that is, a connection relationship does not exist; when the throat cross section between any two adjacent tetrahedral gaps is not completely filled with spherical particles and gaps exist, the connection relationship between the two adjacent tetrahedral gaps is a communication relationship, and the like, and the communication relationship between the adjacent tetrahedral gaps is determined through traversal;

the process of judging the connectivity of all tetrahedral pores is as follows:

judging whether the tetrahedral pore is a connected pore or a non-connected pore by adopting a connectivity algorithm, and specifically comprises the following steps:

the method comprises the following steps: numbering each tetrahedron and storing the sequence in an array;

step two: traversing and searching the communication relation among all adjacent tetrahedral pores;

step three: after traversing and inputting the communication relation between the adjacent tetrahedral pores, modifying the numbers of the two adjacent tetrahedral pores and other tetrahedral pores with the same number as the two adjacent tetrahedrons into the same value when the two adjacent tetrahedral pores are communicated;

step four: dividing a plurality of tetrahedral pores into different pore sets according to the numbering of the tetrahedral pores, wherein all the tetrahedral pores in the same pore set have the same numbering;

step five: for the aperiodic repeated sphere accumulation system, checking whether two pores in opposite surfaces exist in each pore set one by one, and when two pores in opposite surfaces exist in one pore set, indicating that the pores in the pore set are communicated pores; when two pores in opposite surfaces do not exist in one pore set, the pores in the pore set are non-connected pores;

step six: checking whether two pores meeting the periodic repetition condition exist in a plurality of pore sets one by one according to the ball stacking volume obeying the periodic repetition boundary, and when two pores meeting the periodic repetition condition exist in one pore set, indicating that all pores in the pore set are communicated pores; and when two pores meeting the periodically repeated condition do not exist in one pore set, indicating that all the pores in the pore set are non-connected pores, and repeating the steps in sequence to finish the connectivity judgment process of all the Delaunay tetrahedral pores.

6. The method for analyzing the pore connectivity of the cement-based material particle packing system according to claim 5, wherein the method comprises the following steps: the connected porosity of the sphere particle accumulation model is calculated as follows:

after the connectivity of all the tetrahedral pores is judged, the volumes of the tetrahedral pores judged as the connected pores are accumulated and summed, and the total volume of the tetrahedral pores of the connected pores is V_cp(ii) a The communicated porosity of the sphere particle accumulation model is V_cpDividing by the total volume of the spherical particle accumulation model; the porosity of the sphere particle accumulation model is V_cpAnd dividing by the total pore volume of the spherical particle stacking model to obtain the connected porosity and the pore connectivity of the spherical particle stacking model.

Background

The cement-based composite material has the property of porosity, and the porosity has important influence on the mechanical property and the durability of the material, wherein the durability such as freeze thawing, carbonization, steel bar corrosion and the like mostly cause erosion media to enter the material through the transmission of a pore structure, so that the deterioration of the internal structure and the loss of the performance are gradually caused. The transmission rate of the erosion medium into the cement-based multiphase porous material is mainly determined by the open-connected porosity of the porous material itself, in addition to being affected by the ambient temperature. The accurate evaluation of the pore connectivity of the porous material is the key for predicting the transmission capability of the porous material for resisting erosion media, is a precious basis for the subsequent prediction of the durability of the cement-based material, and has very important engineering practical significance. Nowadays, the common experimental methods for measuring the pore structure of porous materials are mainly mercury pressure method and adsorption method. Due to the limitation of the test technology and method principle, the pores which can be tested by the mercury intrusion method and the adsorption method only comprise the communicating pores and part of semi-communicating pores, and the non-communicating pores cannot be detected, so that the pore connectivity of the material cannot be accurately evaluated.

Thanks to the rapidly developing digital modeling techniques, the microscopic three-dimensional pore structure inside porous media materials can be characterized by numerical models, relying on computer image analysis techniques. The three-dimensional model of the porous material (reconstructed according to the structure of a scanning electron microscope or reconstructed according to the structure of the material) is subjected to pixelization, namely, the model is divided into a plurality of congruent cubic pixels, the pixel is considered to be a solid if the central point of the pixel is a solid, and the pixel is considered to be a pore if the central point of the pixel is a pore, so that the pore structure information of the porous material can be analyzed, and the information of the pore structure of the porous material also includes the connected porosity. However, it is known that the connectivity of porous materials is mostly affected by small pores, and the minimum pore size that can be identified based on a three-dimensional model of the porous material after pixelation is affected by the limitation of resolution. The higher the resolution, the finer the porous material three-dimensional model pore structure is displayed, the more accurate the calculation result is, but the content of operation processing is increased at the same time, and the efficiency is reduced; although the result can be obtained by fast calculation under low resolution, the accuracy is greatly reduced due to the rough three-dimensional model of the porous material. The pixelation method is adopted for the irregular porous material, and for the sphere-packed porous material with a special structure, no pore connectivity analysis method which is not limited by resolution exists at present, so that the pore connectivity of the non-connected pores in the irregular porous material is difficult to accurately evaluate.

Disclosure of Invention

The invention aims to provide a method for analyzing the pore connectivity of a cement-based material particle stacking system.

The technical scheme adopted by the invention for solving the technical problems is as follows:

a method for analyzing the pore connectivity of a cement-based material particle stacking system comprises the steps of firstly establishing a spherical particle stacking model, secondly performing Delaunay tetrahedron subdivision on a spherical center point set to divide the stacking model into a set consisting of a plurality of tetrahedrons, then acquiring the pore volume of the tetrahedrons in the set one by one, judging the connectivity of each tetrahedron pore, and finally performing a statistical calculation process of the connected porosity in the spherical particle stacking model.

As a preferable scheme: the process of establishing the sphere particle accumulation model comprises the following steps:

for a system with specified spherical particle size distribution and porosity, firstly determining the size of a cubic filling space, then sequentially putting spherical particles into the cubic filling space from large to small according to the size, wherein the coordinate of each spherical particle is randomly generated by a Monte Carlo method until all the spherical particles are put in, and forming a spherical particle stacking model.

As a preferable scheme: based on a sphere particle accumulation model, the operation process of performing Delaunay tetrahedron subdivision on the sphere center point set is as follows:

according to the information of the sphere particle accumulation model, coordinates of centers of all spheres are arranged to form a two-dimensional matrix, the two-dimensional matrix has three rows which are respectively an abscissa, an ordinate and a column coordinate of the sphere center, each row in the matrix represents a sphere coordinate of one sphere, and the two-dimensional matrix is acted by a Delaunay Triangulation function in MATLAB, so that Delaunay tetrahedron subdivision of a sphere center point set is realized; in the two-dimensional case, the matrix has only two columns, respectively set of abscissas and sets of column coordinates of the circle, on which the Delaunay Triangulation function is applied to perform a Delaunay Triangulation of a set of circle center points.

As a preferable scheme: after the Delaunay tetrahedron subdivision is performed on the sphere center point set, the stacking model is divided into a set consisting of a plurality of tetrahedrons, and the process of determining the pore volumes of the tetrahedrons in the set one by one is as follows:

the particle stacking system after being subdivided by the Delaunay tetrahedron is composed of a plurality of mutually independent and non-overlapping tetrahedrons, wherein each tetrahedron is composed of a solid phase occupied by four vertex angle spherical particles and a pore phase, the solid phase is the sum of sector areas intersected by the three corner point particles on the surface, the pore phase is the residual phase of the solid phase occupied by the four vertex angle spherical particles subtracted from the tetrahedron, the pore phase of the tetrahedron is the pore of the tetrahedron, and the corresponding volume of the pore is the pore volume of the tetrahedron;

for a sphere stacking system allowing overlapping existence, the following scheme is adopted for calculation due to the complex structure:

randomly putting N particles into a tetrahedron, generating the particle coordinates by a Monte Carlo method, counting the number N of the particles falling in the pores, and determining the volume V of the pores in the tetrahedron_pCan be calculated by the formula one:

V_p＝(n/N)*V_t (1)。

as a preferable scheme: the process of judging the connectivity of each tetrahedral pore is to firstly judge the connectivity of all the tetrahedral pores after traversing and judging the connectivity of the adjacent tetrahedral pores;

the tetrahedron is provided with four triangular surfaces, each triangular surface consists of a solid phase and a pore phase, wherein the surface solid phase is the sum of sector areas intersected by particles with three angular points on the surface, and the surface pore phase is the residual phase of the tetrahedron surface minus the corresponding solid phase; the surface pore phase is the cross section of the narrowest part of a channel connecting the tetrahedral pore and the adjacent tetrahedral pore, namely the throat cross section;

when the throat cross section between any two adjacent tetrahedral voids is completely filled with spherical particles, the connection relationship between the two adjacent tetrahedral voids is shown to be a mutual blocking relationship, that is, a connection relationship does not exist; when the throat cross section between any two adjacent tetrahedral gaps is not completely filled with spherical particles and gaps exist, the connection relationship between the two adjacent tetrahedral gaps is a communication relationship, and the like, and the communication relationship between the adjacent tetrahedral gaps is determined through traversal;

the process of judging the connectivity of all tetrahedral pores is as follows:

judging whether the tetrahedral pore is a connected pore or a non-connected pore by adopting a connectivity algorithm, and specifically comprises the following steps:

the method comprises the following steps: numbering each tetrahedron and storing the sequence in an array;

step two: traversing and searching the communication relation among all adjacent tetrahedral pores;

step three: after traversing and inputting the communication relation between the adjacent tetrahedral pores, modifying the numbers of the two adjacent tetrahedral pores and other tetrahedral pores with the same number as the two adjacent tetrahedrons into the same value when the two adjacent tetrahedral pores are communicated;

step four: dividing a plurality of tetrahedral pores into different pore sets according to the numbering of the tetrahedral pores, wherein all the tetrahedral pores in the same pore set have the same numbering;

step five: for the aperiodic repeated sphere accumulation system, checking whether two pores in opposite surfaces exist in each pore set one by one, and when two pores in opposite surfaces exist in one pore set, indicating that the pores in the pore set are communicated pores; when two pores in opposite surfaces do not exist in one pore set, the pores in the pore set are non-connected pores;

step six: checking whether two pores meeting the periodic repetition condition exist in a plurality of pore sets one by one according to the ball stacking volume obeying the periodic repetition boundary, and when two pores meeting the periodic repetition condition exist in one pore set, indicating that all pores in the pore set are communicated pores; and when two pores meeting the periodically repeated condition do not exist in one pore set, indicating that all the pores in the pore set are non-connected pores, and repeating the steps in sequence to finish the connectivity judgment process of all the Delaunay tetrahedral pores.

As a preferable scheme: the connected porosity of the sphere particle accumulation model is calculated as follows:

after the connectivity of all the tetrahedral pores is judged, the volumes of the tetrahedral pores judged as the connected pores are accumulated and summed, and the total volume of the tetrahedral pores of the connected pores is V_cp(ii) a The communicated porosity of the sphere particle accumulation model is V_cpDividing by the total volume of the spherical particle accumulation model; the porosity of the sphere particle accumulation model is V_cpAnd dividing by the total pore volume of the spherical particle stacking model to obtain the connected porosity and the pore connectivity of the spherical particle stacking model.

Compared with the prior art, the invention has the following beneficial effects:

the invention relates to a method for analyzing pore connectivity of a sphere stacking system based on Delaunay triangulation. The method has the advantages of scientific and reasonable calculation principle, simple steps, no requirement on the operation experience of operators, and quick, stable and accurate time effect of the obtained result.

Secondly, the method is based on the intrinsic geometric relationship of sphere accumulation, avoids the method of obtaining high precision by means of high resolution in a pixelation processing method, can quickly divide the tetrahedral pore area and identify connected pores and non-connected pores, and is high in algorithm execution speed and calculation efficiency.

The method can process a multistage unequal-diameter sphere accumulation system, fills the blank that the acquisition of internal communicated pores and non-communicated pores of the conventional multistage unequal-diameter sphere accumulation system is difficult to obtain accurately, obtains the communicated porosity and the pore connectivity of a sphere particle accumulation model through the calculation of the communicated pores and the non-communicated pores, and provides an effective acquisition method of relevant data for the further research of the subsequent multistage unequal-diameter sphere accumulation system.

The method is suitable for calculating the accumulation condition of various particles with low porosity accumulation or high porosity accumulation, the calculation principle is reasonable and comprehensive, the calculation result is accurate, and the actual requirement is met better.

The method can be further expanded to the calculation of the pore size distribution of an irregular polyhedral particle stacking system, and has wide application scenes and potential space for future development and optimization.

And sixthly, the calculation result of the pore connectivity analysis method provided by the invention can be used as a basis for judging the erosion medium invasion resistance of the cement-based material, and a foundation is laid for the subsequent prediction of the durability of the cement-based material.

Drawings

The invention will be further described with reference to the accompanying drawings

FIG. 1 is a block flow diagram of the present invention;

FIG. 2 is a schematic diagram of the three-dimensional structure of a single tetrahedron composed of four spherical particles after three-dimensional Delaunay triangulation;

FIG. 3a is a schematic structural diagram of adjacent tetrahedral pores in a interconnected relationship;

FIG. 3b is a schematic structural diagram of adjacent tetrahedral pores being in a non-interconnected relationship;

FIG. 4 is a schematic diagram of a tetrahedral pore connectivity determination algorithm;

FIG. 5 is a graph of the relationship between single particle diameter overlapping spheres packing interconnected porosity and total porosity;

FIG. 6 is a graph of single particle diameter overlapping spheres packing porosity and total porosity.

Detailed Description

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiment is only one specific embodiment of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

It should be noted that, in order to avoid obscuring the present invention with unnecessary details, only the processing steps closely related to the solution of the present invention are shown in the drawings, and the minor details not related to the present invention are omitted.

The first embodiment is as follows: the method for analyzing pore connectivity of a sphere stacking system in the present embodiment is described with reference to fig. 1, fig. 2, fig. 3a, fig. 3b, fig. 4, fig. 5, and fig. 6, and the method for analyzing pore connectivity of a sphere stacking system in the present embodiment is to first establish a sphere particle stacking model, then perform Delaunay tetrahedron subdivision on a sphere center point set to divide the stacking model into a set composed of a plurality of tetrahedrons, then obtain pore volumes of the tetrahedrons in the set one by one and determine connectivity of pores of each tetrahedron, and finally perform a statistical calculation process of connected porosity in the sphere particle stacking model.

In this embodiment, the sphere particle stacking model is a three-dimensional model of stacking random spheres that can be overlapped.

In the present embodiment, the set of the plurality of tetrahedrons is a result of subdividing the sphere center point set by the Delaunay tetrahedron.

The second embodiment is as follows: the embodiment is further limited by the first embodiment, and the process of establishing the sphere particle stacking model in the sphere stacking system pore connectivity analysis method is as follows:

for a system with specified spherical particle size distribution and porosity, firstly determining the size of a cubic filling space, then sequentially putting the spherical particles into the cubic filling space from large to small according to the size of the spherical particles, and randomly generating coordinates of each spherical particle by a Monte Carlo method until all the spherical particles are put into the cubic filling space to form a spherical particle stacking model.

The third concrete implementation mode: the present embodiment is further limited to the first specific embodiment, a process of arranging coordinates of centers of all spheres to form a two-dimensional matrix in the present embodiment is an existing calculation method, and a processing principle of the Delaunay tetrahedron subdivision is the same as that of the existing Delaunay tetrahedron subdivision.

Step S1 is to realize the specific process of random distribution modeling of spherical particles by using the monte carlo algorithm according to the known spherical particle size distribution information: and sequencing N spheres to be generated from large to small according to the radius, randomly generating the ith sphere when i is greater than the radius, then carrying out the next operation, namely enabling i to be i +1, randomly generating the ith sphere again when i is less than or equal to N, repeating the detection operation, and finishing the generation of all spheres when i is greater than N.

The fourth concrete implementation mode: the present embodiment is further limited to the first, second, or third embodiment, and based on the sphere particle stacking model, the operation process of performing Delaunay tetrahedron subdivision on the sphere center point set is as follows:

according to the information of the sphere particle accumulation model, coordinates of centers of all spheres are arranged to form a two-dimensional matrix, the two-dimensional matrix has three rows which are respectively an abscissa, an ordinate and a column coordinate of the sphere center, each row in the matrix represents a sphere coordinate of one sphere, the two-dimensional matrix is acted by a Delaunay Triangulation function in MATLAB, and Delaunay tetrahedron subdivision is performed on a sphere center point set; in the two-dimensional case, the matrix has only two columns, respectively set of abscissas and sets of column coordinates of the circle, on which the Delaunay Triangulation function is applied to perform a Delaunay Triangulation of a set of circle center points.

In the present embodiment, the set of spherical center points is a set in which the spherical coordinates of all spherical particles are combined.

The result of the two-dimensional Delaunay triangulation in this embodiment is schematically shown in fig. 2, where a circle represents a particle, and a straight line segment set is a division result of a circle center point set.

The fifth concrete implementation mode: the embodiment is further defined as the first, second, third or fourth embodiment, and the step of determining the pore diameter and pore volume corresponding to each tetrahedron after performing Delaunay tetrahedron subdivision on the sphere center point set is as follows:

the particle stacking system after being subdivided by the Delaunay tetrahedron is composed of a plurality of mutually independent and non-overlapping tetrahedrons, wherein each tetrahedron is composed of two parts: the solid phase portion occupied by the apex angle spherical particles and the pore phase other than the solid phase portion. The tetrahedral mesoporosity may be referred to as tetrahedral porosity, and the corresponding volume is the pore volume of the tetrahedron.

The sixth specific implementation mode: the embodiment is a further limitation of the first, second, third, fourth, or fifth embodiment, after Delaunay tetrahedron subdivision is performed on the sphere center point set, the stacking model is divided into a set composed of a plurality of tetrahedrons, and the process of determining the pore volume of the tetrahedrons in the Delaunay tetrahedron subdivision sphere center point set one by one is as follows:

the particle stacking system after being subdivided by the Delaunay tetrahedron is composed of a plurality of mutually independent and non-overlapping tetrahedrons, wherein each tetrahedron comprises a solid phase volume occupied by the vertex angle spherical particles and a pore volume, the pore volume is a difference value of the solid phase volume occupied by the vertex angle spherical particles subtracted from the total volume of the tetrahedron, the pore volume is a tetrahedral pore, and the volume corresponding to the pore volume is the pore volume of the tetrahedral pore;

for a sphere stacking system allowing overlapping existence, the following scheme is adopted for calculation due to the complex structure:

randomly putting N particles into a tetrahedron, generating the particle coordinates by a Monte Carlo method, counting the number N of the particles falling in the pores, and determining the volume V of the pores in the tetrahedron_pCan be calculated by the formula one:

V_p＝(n/N)*V_t (1)

thus, the Delaunay tetrahedral pore volume is calculated in a traversal mode.

The seventh embodiment: the embodiment is further limited to the first, second, third, fourth, fifth or sixth specific embodiments, and the process of judging the connectivity of each tetrahedral pore includes first traversing and judging the connectivity of adjacent tetrahedral pores, and then judging the connectivity of all tetrahedral pores;

the process of traversing and judging the pore connectivity of adjacent tetrahedrons comprises the following steps: the tetrahedron is provided with four triangular surfaces, each triangular surface consists of a solid phase and a pore phase, wherein the surface solid phase is the sum of sector areas intersected by particles with three angular points on the surface, and the surface pore phase is the residual phase of the tetrahedron surface minus the corresponding solid phase; the surface pore phase is the cross section of the narrowest part of a channel connecting the tetrahedral pore and the adjacent tetrahedral pore, namely the throat cross section;

when the throat cross section between any two adjacent tetrahedral voids is completely filled with spherical particles, the connection relationship between the two adjacent tetrahedral voids is shown to be a mutual blocking relationship, that is, a connection relationship does not exist; when the throat cross section between any two adjacent tetrahedral gaps is not completely filled with spherical particles and gaps exist, the connection relationship between the two adjacent tetrahedral gaps is a communication relationship, and the like, and the communication relationship between the adjacent tetrahedral gaps is determined through traversal;

the process of judging the connectivity of all tetrahedral pores is as follows:

judging whether the tetrahedral pore is a connected pore or a non-connected pore by adopting a connectivity algorithm, and specifically comprises the following steps:

the method comprises the following steps: numbering each tetrahedron and storing the sequence in an array;

step two: traversing and searching the communication relation among all adjacent tetrahedral pores;

step three: after traversing and inputting the communication relation between the adjacent tetrahedral pores, when the two adjacent tetrahedral pores are communicated, modifying the numbers of the two adjacent tetrahedral pores and other tetrahedral pores with the same number as the two adjacent tetrahedrons into the same value;

step four: dividing a plurality of tetrahedral pores into different pore sets according to the numbering of the tetrahedral pores, wherein all the tetrahedral pores in the same pore set have the same numbering;

step five: checking whether two pores in opposite surfaces exist in a plurality of pore sets one by one for an aperiodic repeated sphere accumulation system, and when two pores in opposite surfaces exist in one pore set, indicating that the pores in the pore set are connected pores; when two pores in opposite surfaces do not exist in one pore set, the pores in the pore set are non-connected pores;

step six: checking whether two pores meeting the periodic repetition condition exist in each pore set one by one according to the sphere stacking volume obeying the periodic repetition boundary, and when two pores meeting the periodic repetition condition exist in one pore set, indicating that all pores in the pore set are communicated pores; and when two pores meeting the periodically repeated condition do not exist in one pore set, indicating that all the pores in the pore set are non-connected pores, and repeating the steps in sequence to finish the connectivity judgment process of all the Delaunay tetrahedral pores.

The specific implementation mode is eight: the embodiment is further defined by the first, second, third, fourth, fifth, sixth or seventh embodiment, and the calculation process of the interconnected porosity of the sphere particle packing model is as follows:

after the connectivity of all the tetrahedral pores is judged, the volumes of the tetrahedral pores judged as the connected pores are accumulated and summed, and the total volume of the tetrahedral pores of the connected pores is V_cp(ii) a The communicated porosity of the sphere particle accumulation model is V_cpDividing by the total volume of the spherical particle accumulation model; the porosity of the sphere particle accumulation model is V_cpAnd dividing by the total pore volume of the spherical particle stacking model to obtain the connected porosity and the pore connectivity of the spherical particle stacking model.

In this embodiment, the communicating porosity of the spherical particle stacking model is the stacking system communicating porosity, and the pore communication of the spherical particle stacking model is the stacking system pore communication.

The invention relates to a method for analyzing pore connectivity of a cement-based material particle stacking system, which is mainly applied to spherical particles, and other particles with irregular shapes can be simplified into spherical particles and then processed by using the method for analyzing pore connectivity of the particle stacking system.

The specific implementation method nine: the embodiment is further limited by the specific embodiment one, two, three, four, five, six, seven or eight, and is based on the established sphere stacking model, Delaunay tetrahedron subdivision is performed on a sphere center point set, and then the connectivity of the Delaunay tetrahedron pore is determined by adopting the proposed connectivity judgment algorithm, so that the connected porosity of the pore structure is obtained, and the specific implementation steps are as follows:

s1, building a three-dimensional model of the overlapped random sphere stacking: firstly, setting the size of a cubic filling space of particles, then sequentially putting spherical particles into the cubic space, and randomly generating coordinates of the spherical particles by a Monte Carlo method until all the particles are put in or the system reaches a target porosity.

S2, Delaunay tetrahedron subdivision of the sphere center point set: the sphere center coordinates of all piled spheres are arranged into a two-dimensional matrix, the matrix comprises three rows which respectively represent the abscissa, ordinate and column coordinate of the sphere center, each row in the matrix represents the sphere center coordinate of one sphere, and the two-dimensional coordinate matrix is acted by a Delaunay Triangulation function in MATLAB, so that Delaunay tetrahedron subdivision on the sphere center point set is realized.

The particle filling space after being subdivided by the Delaunay tetrahedron consists of mutually non-overlapping tetrahedrons, and each tetrahedron consists of four spherical particles with vertex angles and a pore part; fig. 2 is a composition example of a Delaunay tetrahedron, in which spherical particles are filled in a gray gradient manner, and the portion of the tetrahedron other than the spherical particles is a void, and the corresponding volume of the portion is a pore volume.

S3, calculating the volume of tetrahedral pores: for Delaunay tetrahedrons in a sphere stacking system allowing overlapping existence, after calculating the volume Vt of the tetrahedron by adopting a geometrical relationship, randomly putting N particles into the tetrahedron, randomly generating the particle coordinates by a Monte Carlo method, counting the number N of the particles falling in the pores, and calculating the volume Vp of the tetrahedral pores by the following formula:

Vp＝(n/N)*Vt

s4, connection relation of adjacent tetrahedral pores: the void parts of the four surfaces of the Delaunay tetrahedron correspond to the cross sections of the four throats of the aperture and are channels connecting the tetrahedron aperture with the four surrounding tetrahedron apertures; if the throat cross section channel between two adjacent tetrahedral gaps is completely filled with spherical particles, they are mutually isolated and there is no communication relationship, as illustrated in fig. 3 a; on the contrary, if the throat cross-section channel between two adjacent tetrahedral voids is not completely filled with spherical particles, leaving some voids, there is a communication relationship between them, as illustrated in fig. 3 b.

S5, judging the connectivity of all tetrahedral pores: judging whether the Delaunay tetrahedral pore is a connected pore or a non-connected pore by adopting a connectivity algorithm, wherein the algorithm is schematically shown in figure 4, and the specific steps are analyzed as follows:

numbering each tetrahedron and storing the tetrahedrons in an array in sequence;

traversing and searching the communication relation among all adjacent tetrahedral pores;

traversing and inputting the communication relation between the adjacent tetrahedral pores, and modifying the numbers of the two adjacent tetrahedral pores and other tetrahedral pores with the same number as the two adjacent tetrahedral pores into the same value if the two adjacent tetrahedral pores are communicated;

dividing the tetrahedral pores into different pore sets according to the serial numbers of the tetrahedral pores, wherein the serial numbers of all the tetrahedral pores in the same pore set are the same;

checking whether two pores in opposite surfaces exist in the pore set or not for the aperiodic repeated sphere accumulation system, if so, the pores in the pore set are connected pores, otherwise, the pores are not connected pores;

checking whether two pores meeting the periodic repetition condition exist in the pore set or not for the ball stacking volume obeying the periodic repetition boundary, wherein if the two pores exist, all the pores in the pore set are connected pores, and if not, the pores are not connected pores;

and seventhly, finishing the connectivity judgment of all Delaunay tetrahedron holes.

S6, calculating the interconnected porosity of the stacking system: the volume of the tetrahedral pore which is judged to be the connected pore in S4 is added up and is recorded as V_cp(ii) a The interconnected porosity of the packing system is V_cpDividing by the total volume of the stacked system; the pore connectivity of the stacking system is V_cpDivided by the total pore volume of the packing system.

The following embodiments are described in conjunction with the advantageous effects of the present invention and the accompanying drawings 1 to 6 of the specification:

the first embodiment is as follows:

selecting a sphere stacking system allowing overlapping of single particle sizes, establishing sphere stacking three-dimensional models with different porosities by adopting a Monte Carlo method, then performing Delaunay tetrahedron division, tetrahedral pore volume traversal calculation, tetrahedral pore connectivity judgment, connected porosity and pore connectivity calculation on each model, finally sorting, analyzing the influence of porosity change on pore connectivity in the system, and displaying the operation and calculation results of the algorithm:

the relationship between the interconnected porosity and the total porosity of the stacking system is shown in fig. 5, fig. 5 shows the calculation result of the interconnected porosity of the stacking system with different porosities, and the line of the porosity-interconnected porosity tends to a straight line with a slope of 1, which indicates that for most stacking systems with different porosities, the internal pores are all interconnected pores; however, in the case of extremely low porosity, i.e., porosity lower than 0.1, the straight line is slightly bent, and a convex trend is exhibited, indicating that some non-connected pores are present in the packed system at this time.

The relationship between the pore connectivity and the total porosity of the stacking system is shown in fig. 6, fig. 6 shows the calculation of the pore connectivity of the stacking system with different porosities, and the pore connectivity of the stacking system of spherical particles with single particle size and the porosity of more than 0.1, which is a stacking system of spherical particles which can be overlapped, is 1, which shows that all pores inside the stacking system are connected pores and have predictable good transmission performance; for a single-particle-size overlappable sphere particle stacking system with the porosity of less than 0.1, along with the reduction of the porosity, the porosity in the system is rapidly reduced, the original connected pores in the system are gradually converted into unconnected pores, the porosity of the system is already 0 by 0.015, all the pores are changed into the unconnected pores, and the system loses the material transmission capability.

In the description herein, references to the description of "one embodiment," "an example," "a specific example" or the like are intended to mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed.