CN108280290B - Concrete aggregate numerical model reconstruction method - Google Patents
Concrete aggregate numerical model reconstruction method Download PDFInfo
- Publication number
- CN108280290B CN108280290B CN201810059964.1A CN201810059964A CN108280290B CN 108280290 B CN108280290 B CN 108280290B CN 201810059964 A CN201810059964 A CN 201810059964A CN 108280290 B CN108280290 B CN 108280290B
- Authority
- CN
- China
- Prior art keywords
- aggregate
- images
- model
- image
- test piece
- 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.)
- Active
Links
- 239000004567 concrete Substances 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000012360 testing method Methods 0.000 claims abstract description 54
- 238000005498 polishing Methods 0.000 claims abstract description 25
- 238000004088 simulation Methods 0.000 claims abstract description 18
- 238000011160 research Methods 0.000 claims abstract description 14
- 239000011229 interlayer Substances 0.000 claims abstract description 8
- 239000000178 monomer Substances 0.000 claims abstract description 4
- 238000013499 data model Methods 0.000 claims description 31
- 239000000463 material Substances 0.000 claims description 19
- 239000010410 layer Substances 0.000 claims description 17
- 238000004519 manufacturing process Methods 0.000 claims description 11
- 238000000227 grinding Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000013507 mapping Methods 0.000 claims description 9
- 238000005520 cutting process Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000009837 dry grinding Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000011049 filling Methods 0.000 claims description 3
- 230000002452 interceptive effect Effects 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 239000011083 cement mortar Substances 0.000 claims description 2
- 229910052602 gypsum Inorganic materials 0.000 claims description 2
- 239000010440 gypsum Substances 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims description 2
- 238000002591 computed tomography Methods 0.000 abstract description 9
- 239000002245 particle Substances 0.000 abstract description 7
- 238000005516 engineering process Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 239000004570 mortar (masonry) Substances 0.000 description 4
- 239000004568 cement Substances 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 239000010802 sludge Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- -1 and generally Substances 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Processing Or Creating Images (AREA)
Abstract
The invention belongs to the technical field of concrete numerical simulation research, and relates to a method for reconstructing a concrete aggregate numerical model, polishing and photographing a test piece containing monomer aggregate or a test piece containing group aggregate to obtain a continuous section image of aggregate, establishing an aggregate model based on the continuous section image, an aggregate model database is established according to aggregate models, a plurality of aggregate models can be established simultaneously by containing aggregate test pieces, the problem that CT scanning equipment is expensive and inconvenient to implement can be solved, and three-dimensional reconstruction is carried out based on real concrete aggregate particles, the appearance characteristics of the real aggregate are simulated to the maximum extent, the established aggregate model is stored in a grid form, creating conditions for establishing a real aggregate model database, and exchanging formats with CAD software or FEA finite element simulation software, so as to facilitate numerical simulation research on the mechanical properties of concrete; the image acquisition mode is simple and quick, and the acquired image sequence has small interlayer spacing, high alignment degree and complete information quantity.
Description
The technical field is as follows:
the invention belongs to the technical field of concrete numerical simulation research, and relates to a concrete aggregate numerical model reconstruction method which can quickly and accurately acquire multi-phase material small-space continuous section images.
Background art:
concrete (concrete for short) is a general term for engineering composite materials formed by cementing aggregate into a whole by cementing materials, and generally, concrete is composed of coarse aggregate, fine aggregate, cement hydrate, unhydrated cement particles, pores, cracks and the like; the cement concrete, also called as common concrete, is obtained by mixing with water (which may contain additives and admixtures) according to a certain proportion and stirring, and is widely applied to civil engineering. At present, the research on concrete is only limited to a macroscopic level, along with the continuous development of computer technology, the research on the aspect of concrete mesomechanics is increasingly emphasized by people, the damage forms of aggregate, mortar and an interface layer between the aggregate and the mortar in the concrete can be researched after the concrete is subjected to numerical simulation, a finite element geometric model of a concrete material mesostructure is established on the premise of adopting numerical method analysis, and the performance of the concrete has a great relation with the shape, the gradation and the like of the aggregate, so the generation of the concrete aggregate model shape becomes an important subject for the concrete numerical simulation.
The mesoscopic numerical simulation can replace part of experiments under the conditions that a theoretical model is reasonable and the characteristic data of each phase of concrete material are accurate enough, can avoid the objective limit of experimental conditions and the influence of human factors on the result, adopts a plurality of assumptions for the establishment of the concrete mesoscopic model, and has certain errors: for example, the key problem of concrete numerical simulation is the simulation of the shape of concrete aggregate, in the past research, for two-dimensional concrete aggregate, circular or elliptical simulated pebble or gravel aggregate is mostly used, broken stone aggregate is simulated by convex polygon, pebble or gravel aggregate is simulated by spherical or elliptical shape in most three dimensions, and broken stone aggregate is simulated by convex polygon; the simplified modes are feasible on the method, but the morphological characteristics of real concrete aggregate particles cannot be accurately expressed, although the two-dimensional simulation is simple and convenient to calculate, the defects of the two-dimensional simulation are obvious, the spatial distribution and grading characteristics of the concrete aggregate cannot be really restored, and the damage characteristics of the concrete cannot be expressed; in the three-dimensional simulation, the random placement of the spherical, the ellipsoidal and the convex polyhedrons can simulate the spatial distribution form of concrete aggregate, can not completely replace the shape of real aggregate, and has defects in the research aspect of aggregate and mortar interface layers. In order to obtain a more realistic concrete structure and failure characteristics, it is necessary to simulate real concrete aggregates in a more efficient way.
The scanning of the concrete section slices is realized by utilizing an industrial CT scanning technology, then the reconstruction of the concrete model is carried out by using three-dimensional reconstruction software, a concrete numerical model based on a real shape is obtained, and then the analysis is carried out on the basis of the numerical model, so that a better effect can be obtained. However, the application of the CT scanning technology in the concrete field also has disadvantages, when the resolution of the CT scanning image is improved, the density resolution is reduced, which results in insufficient information content, and meanwhile, because the CT scanning technology is greatly influenced by the scanning equipment, noise and artifact interference are easily generated, which causes the reconstructed model to be inaccurate, and the CT scanning equipment is expensive, which is not beneficial to the wide application of the CT scanning technology in the concrete field, the CT scanning technology has certain limitations for the real concrete mesoscopic structure three-dimensional reconstruction.
Chinese patent 201510801544.2 discloses a concrete aggregate structural feature pickup method based on a section shaping method, which is characterized in that an oil sludge module is made of a concrete material for three-dimensional reconstruction, the made oil sludge module is cut and scanned, an obtained scanned image of the oil sludge module is subjected to image processing, a three-dimensional space model of a concrete particle structure is reconstructed, and finally the obtained three-dimensional space model of the concrete particle structure is randomly distributed in direction and position to realize the pickup of the concrete aggregate structural feature; the method for three-dimensional reconstruction of a concrete CT image disclosed in Chinese patent 201210015951.7 includes firstly carrying out tomography scanning on a concrete sample by using a CT machine, secondly carrying out three-dimensional reconstruction on a concrete C T image by using a ray projection algorithm to obtain real images of aggregates, mortar and holes, and finally carrying out three-dimensional reconstruction on the concrete CT image by using a visualization tool VTK.
Based on the reasons, a more accurate and convenient concrete aggregate numerical model reconstruction method is needed to be found, an aggregate model database based on the real concrete aggregate shape is further established by using the numerical model, and aggregates are conveniently and quickly extracted from the model database and randomly put in the model database in the numerical simulation research of concrete.
The invention content is as follows:
the invention aims to overcome the defects in the prior art, develops and designs a method for carrying out model reconstruction on the basis of continuous sectional images obtained by polishing real concrete aggregate layer by layer, and further generates an aggregate model database by utilizing an aggregate model established by the method. To increase the effectiveness of the numerical simulation of the concrete.
In order to achieve the aim, the concrete aggregate numerical model reconstruction method comprises six steps of manufacturing a polishing test piece containing target aggregate, obtaining an aggregate continuous section image, calibrating and preprocessing the continuous section image, generating an aggregate model, calculating aggregate model parameters, performing format conversion and establishing an aggregate model database:
preparing a grinding test piece containing a target aggregate: manufacturing a test piece containing a monomer aggregate or a test piece containing a group aggregate according to the requirements established by the aggregate model database; the manufacturing process of the single aggregate-containing test piece is to wrap a single aggregate in a cementing material to manufacture a polishing test piece, the cementing material is a quick-hardening material convenient for polishing, the quick-hardening material comprises cement mortar and building gypsum, and the single aggregate-containing test piece is convenient for the convenient extraction of the aggregate and the reconstruction of a three-dimensional model; the manufacturing process of the test piece containing the group aggregate is the same as that of the common concrete test piece, and the aggregates containing the group aggregate test piece are adhered and need to be segmented and extracted;
(II) acquiring an aggregate continuous section image: fix the aggregate test piece of polishing on grinding equipment, the bistrique is polished the test piece and is carried out the successive layer dry grinding to the aggregate with the same interval of polishing along the Z axle direction of aggregate test piece of polishing, and the interval of polishing of section is controlled through the feed of bistrique, and the minimum interval of polishing is 0.01mm, and every layer of section is polished and is accomplished the back, uses digital camera to shoot the section, acquires aggregate continuous section image:
(III) calibrating and preprocessing continuous section images: the method comprises the steps of calibrating aggregate continuous section images by adopting a three-anchor point alignment method, wherein when an aggregate grinding test piece is ground in a layered mode, the photographed section images of each layer cannot be registered among layers, because the images are influenced by various factors such as the angle, the distance and the direction of a camera, the movement of the position of an object and other factors when being photographed, the images need to be registered, three control points close to a multiphase material test piece on a test piece clamping mechanism are selected as reference points, the three reference points are not on the same straight line, when the continuous section images are imaged, each section image comprises the three reference points, one standard image is selected as a reference image, the rest section images are to-be-registered images, the standard images are images with the angles, the distances, the directions and the position of the object in the continuous section images meeting the set requirements, and the image registration is a process of seeking one-to-one mapping between the to-be-registered images and the standard images, linking points corresponding to the same spatial position in two images, superposing three control points in an image to be registered with three reference points in a reference image, calculating a spatial mapping relation according to the positions of the control points, carrying out geometric transformation on the image to be registered by using the spatial mapping relation, and carrying out Z-axis direction registration on a continuous section image by using a cpselect function interactive tool in an MATLAB (matrix laboratory) to obtain a registration result; the method comprises the steps of shooting anchor points which do not belong to section images of the porous concrete sample, uniformly cutting the registered continuous section images by utilizing Photoshop (image processing software) and image processing (image analysis software), cutting useless areas including the anchor points around the continuous section images, reserving the continuous section images of the concrete sample polishing area to be analyzed, carrying out gray-scale processing on the cut continuous section images and converting the cut continuous section images into binary images, wherein the white area is aggregate, the black area is gel material, and completing calibration and pretreatment of the continuous section images;
based on the fact that the continuous section images are panoramic images, the imaging equipment parameters for shooting and storing the images are different, and the relative positions and angles of the imaging equipment and the sections of the porous concrete sample are different, so that position deviation exists among the section images of each layer, and the obtained continuous section images need to be subjected to standardization processing to meet the requirements of experimental analysis;
(IV) generating an aggregate model:
(1) in order to ensure that the proportion of image pixels (pixels) to the real aggregate size is consistent, conversion calculation of the image pixels and the real aggregate size is carried out on a continuous sectional image, the scanning resolution of an X axis, a Y axis and a Z axis is input in the process of importing a Mimics (medical image control system) three-dimensional reconstruction software into the continuous sectional image, point coordinates (X ', Y') on the image correspond to corresponding point coordinates (X, Y) on a real slice, the scanning resolution of the X axis and the Y axis is the image real size corresponding to each pixel, and the scanning resolution of the Z axis is the interlayer spacing d of the continuous sectional image;
(2) threshold segmentation, wherein white aggregate areas are extracted to generate aggregate masks, white areas are directly extracted from single aggregate samples, a group of aggregate samples need to judge whether the white areas in different two-dimensional images are from the same aggregate or not, the judgment process is realized through software programming, the centroid distance of any amount of aggregates is compared with the maximum radius of the two aggregates, the centroid distance of the two aggregates is smaller than the radius of any one aggregate to judge that the two aggregates are from the same aggregate, the white areas of the same aggregate are extracted by adopting a method of combining gray value range definition and manual drawing and are stored as the same mask, and the masks of the adhered parts of the different aggregates in the continuous sectional images are manually trimmed and segmented through mask editing operation;
(3) 3D computing the aggregate mask extracted from the single aggregate test piece in the step (2) to generate an initial aggregate surface grid data model, respectively performing 3D computing on different aggregate masks extracted from aggregate test pieces of a group, and simultaneously generating a plurality of initial aggregate surface grid data models, wherein the generated aggregate surface grid data model is a coarse primary model due to the influence of the interlayer spacing of the continuous section images of the aggregate test pieces on the Z-axis resolution;
(4) carrying out fairing, meshing and wrapping treatment on the primary aggregate surface mesh data model generated in the step (3) to obtain an accurate aggregate surface mesh data model, and filling the accurate aggregate surface mesh data model to obtain a real aggregate-based volume mesh data model;
the principle of three-dimensional reconstruction software of the Mimics (medical image control system) for reversely reconstructing the aggregate particle model is as follows: on the basis of three-dimensional information of an unknown object contained in the continuous section images, a digital three-dimensional model of the object outline is reconstructed in a computer, wherein coordinate values (X, Y) of a point P on the unknown object outline on an X axis and a Y axis are determined by the position on a horizontal section of the unknown object outline, and the coordinate value (Z) of the point P on a Z axis is determined by the relation between the cutting distance (d) and the section sequence (n): z is determined as d × n;
and (V) calculating aggregate model parameters and carrying out format conversion: the geometric parameters directly extracted from the aggregate grid data model reconstructed by the Mimics software comprise model volume, surface area, long-axis distance and short-axis distance in the X-axis direction, the Y-axis direction and the Z-axis direction and mass barycentric coordinates determined based on a model coordinate frame, the aggregate model indirect shape parameters including aggregate shape factor, roundness, needle degree and flatness can be calculated through the directly extracted geometric parameters, further, the aggregate grid data model can be converted into a format which is convenient for storing an aggregate model database and is dxf, stl, inp and igs through corresponding interfaces of various application design software provided by the Mimics software, and the aggregate model format is convenient for putting and researching;
and (VI) establishing an aggregate model database: and (5) establishing aggregate grid data models with different sizes and types by utilizing the steps from (I) to (V), carrying out aggregate model parameter statistics and classification integration on the aggregate grid data models with different sizes and types, establishing an aggregate model database, and selecting and putting aggregates in subsequent numerical simulation research.
The interlayer spacing of the aggregate continuous section image related to the step (one) influences the resolution of the aggregate model reconstruction result in the Z-axis direction, and the smaller the interlayer spacing of the aggregate continuous section image is, the more accurate the fitting effect of the external form of the real concrete aggregate model is; the aggregate continuous section image obtained in the step (one) has high resolution and complete contained information, a multiphase and multi-scale numerical model can be reconstructed, different levels of research are facilitated, and the three-dimensional reconstruction of the concrete microstructure based on the aggregate continuous section image is endowed with great flexibility; the aggregate model database established in the step (six) has the characteristics of high degree of reality and convenient and quick aggregate selection, and provides a solid foundation for numerical simulation research of concrete aggregates.
Compared with the prior art, the method has the advantages that the test pieces containing the monomer aggregate or the group aggregate are polished and photographed to obtain the continuous section images of the aggregate, the aggregate model is built based on the continuous section images, the aggregate model database is built according to the aggregate model, a plurality of aggregate models can be built simultaneously for the test pieces containing the group aggregate, the problem that the CT scanning equipment is expensive and inconvenient to implement can be solved, the three-dimensional reconstruction is carried out based on real concrete aggregate particles, the appearance characteristics of the real aggregate are simulated to the maximum extent, the built aggregate model is stored in a grid form, conditions are created for building the real aggregate model database, the format exchange can be carried out with CAD software or FEA finite element simulation software, and the numerical simulation research on the mechanical property of the concrete is facilitated; the method has the advantages of simple and quick image acquisition mode, low image acquisition cost, small interlayer spacing of the acquired image sequence, high alignment degree and complete information quantity, and the established aggregate model database has the characteristics of reality and reliability and is convenient to apply and popularize.
Description of the drawings:
FIG. 1 is a block diagram of the process flow of the present invention.
Fig. 2 is a sectional image before binarization according to step (iii) of embodiment 1 of the present invention.
Fig. 3 is a binarized sectional image obtained in step (iii) of example 1 of the present invention.
FIG. 4 is a diagram showing an aggregate primary structure model according to step (IV) in example 1 of the present invention.
FIG. 5 is a schematic diagram of an aggregate surface mesh pattern according to step (IV) in example 1 of the present invention.
Fig. 6 is a model diagram of an aggregate mesh according to step (iv) in example 1 of the present invention.
FIG. 7 is a model diagram of a transformed CAD aggregate according to step (V) of example 1 of the present invention.
FIG. 8 is a schematic representation of an aggregate having a shape described as biased toward a flat type according to example 2 of the present invention.
FIG. 9 is a schematic representation of an aggregate having a shape described as somewhat flattened according to example 2 of the present invention.
FIG. 10 is a schematic representation of an aggregate having a shape described as a somewhat needle-like shape according to example 2 of the present invention.
FIG. 11 is a schematic representation of an aggregate having a shape described as biased toward acicular according to example 2 of the present invention.
The specific implementation mode is as follows:
the invention is further illustrated by the following examples in conjunction with the accompanying drawings.
Example 1:
the concrete aggregate numerical model reconstruction method related to the embodiment comprises five steps of manufacturing a grinding test piece containing target aggregate, obtaining an aggregate continuous section image, calibrating and preprocessing the continuous section image, generating an aggregate model, calculating parameters of the aggregate model and carrying out format conversion:
preparing a grinding test piece containing a target aggregate: casting a cubic concrete test block of 100mm multiplied by 100mm in situ, and placing the cubic concrete test block in a standard curing room for curing for 7 days to obtain a concrete test piece containing group aggregate;
(II) acquiring an aggregate continuous section image: fix the aggregate test piece of polishing on grinding equipment, the bistrique is polished the test piece and is carried out the successive layer dry grinding with the interval of polishing of 1mm to the aggregate along the Z axle direction of aggregate test piece of polishing, and the interval of polishing of section is controlled through the feeding of bistrique, and every layer section is polished and is accomplished the back, uses digital camera to shoot the section, acquires the interval and is the aggregate continuous section image of 1 mm:
(III) calibrating and preprocessing continuous section images: the method comprises the steps of calibrating aggregate continuous section images by adopting a three-anchor point alignment method, wherein when an aggregate grinding test piece is ground in a layered mode, the photographed section images of each layer cannot be registered among layers, because the images are influenced by various factors such as the angle, the distance and the direction of a camera, the movement of the position of an object and other factors when being photographed, the images need to be registered, three control points close to a multiphase material test piece on a test piece clamping mechanism are selected as reference points, the three reference points are not on the same straight line, when the continuous section images are imaged, each section image comprises the three reference points, one standard image is selected as a reference image, the rest section images are to-be-registered images, the standard images are images with the angles, the distances, the directions and the position of the object in the continuous section images meeting the set requirements, and the image registration is a process of seeking one-to-one mapping between the to-be-registered images and the standard images, linking points corresponding to the same spatial position in two images, superposing three control points in an image to be registered with three reference points in a reference image, calculating a spatial mapping relation according to the positions of the control points, carrying out geometric transformation on the image to be registered by using the spatial mapping relation, and carrying out Z-axis direction registration on a continuous section image by using a cpselect function interactive tool in an MATLAB (matrix laboratory) to obtain a registration result; the method comprises the steps of shooting anchor points which do not belong to section images of the porous concrete sample, uniformly cutting the registered continuous section images by utilizing Photoshop (image processing software) and image processing (image analysis software), cutting useless areas including the anchor points around the continuous section images, reserving the continuous section images of the concrete sample polishing area to be analyzed, carrying out gray-scale processing on the cut continuous section images and converting the cut continuous section images into binary images, wherein the white area is aggregate, the black area is gel material, and completing calibration and pretreatment of the continuous section images;
(IV) generating an aggregate model:
(1) in order to ensure that the proportion of image pixels (pixels) to the real aggregate size is consistent, conversion calculation of the sizes of the image pixels and the real aggregate is carried out on continuous sectional images, the scanning resolutions of an X axis, a Y axis and a Z axis are input in the process of importing the continuous sectional images into Mimics (medical image control system) three-dimensional reconstruction software, point coordinates (X ', Y') on the images correspond to corresponding point coordinates (X, Y) on real slices, and the calculation formula is used for calculating the ratio of the image pixels (pixels) to the real aggregate size according to the scanning resolutions of the X axis, the Y
And
calculating to obtain x being 0.25mm and y being 0.25mm, wherein the calculation result represents the image size contained in each pixel of the image, the scanning resolution of the x axis and the y axis is 0.25mm, and the scanning resolution of the z axis is the interval of 1mm of the section images;
(2) threshold segmentation, wherein white aggregate areas are extracted to generate aggregate masks, white areas are directly extracted from single aggregate samples, a group of aggregate samples need to judge whether the white areas in different two-dimensional images are from the same aggregate or not, the judgment process is realized through software programming, the centroid distance of any amount of aggregates is compared with the maximum radius of the two aggregates, the centroid distance of the two aggregates is smaller than the radius of any one aggregate to judge that the two aggregates are from the same aggregate, the white areas of the same aggregate are extracted by adopting a method of combining gray value range definition and manual drawing and are stored as the same mask, and the masks of the adhered parts of the different aggregates in the continuous sectional images are manually trimmed and segmented through mask editing operation;
(3) respectively carrying out 3D calculation on different masks, fitting the aggregate region boundaries under the same mask into a matrix type primary three-dimensional data model, and taking the primary three-dimensional data model generated by fitting as an initial aggregate surface grid data model;
(4) carrying out fairing, meshing and wrapping treatment on the primary aggregate surface mesh data model generated in the step (3) to obtain an accurate aggregate surface mesh data model, and filling the accurate aggregate surface mesh data model to obtain a real aggregate-based volume mesh data model;
and (V) calculating aggregate model parameters and carrying out format conversion: the aggregate grid data model is converted into a dxf format under CAD through a CAD computer aided drawing software interface provided by the Mimics software so as to be convenient for application of an aggregate model database, or the aggregate grid data model is converted into a format required by finite element analysis through an FEA finite element analysis software interface provided by the Mimics software so as to analyze the mechanical property of the aggregate.
Example 2:
in this embodiment, the aggregate model constructed in example 1 is subjected to extraction of geometric parameters of short axial distance, long axial distance, volume, surface area and mass barycentric coordinates determined based on a model coordinate frame, such as those in X-axis direction, Y-axis direction and Z-axis direction, and the results are shown in table 1 (mm):
indirect shape parameters of aggregate shape factor, roundness, conicity and flatness were calculated from the extracted geometric parameters, and the results are shown in table 2:
| description of shape | Deflection flat type | Slightly flat class | Slightly acicular class | Biased acicular species |
| Form factor | -0.70 | 0.70-1.00 | 1.00-1.40 | 1.40- |
| Roundness degree | -0.88 | 0.88-1.00 | 1.00-1.12 | 1.12- |
| Degree of penetration | 1.43-2.09 | 1.24-1.91 | 1.60-2.30 | 1.80-2.57 |
| Flatness of | 0.13-0.45 | 0.36-0.58 | 0.45-0.98 | 0.45-0.88 |
Integrating, classifying and judging the shape characteristics of the aggregate by the indirect shape parameters; classifying and integrating the aggregate model parameters in the table 1 and the table 2, and establishing a real concrete aggregate model database by the method for randomly putting concrete aggregates and establishing a concrete numerical model, thereby providing beneficial help for numerical simulation research of concrete.
Claims (1)
1. A concrete aggregate numerical model reconstruction method is characterized in that the concrete technological process comprises six steps of manufacturing a polishing test piece containing target aggregate, obtaining aggregate continuous section images, calibrating and preprocessing the continuous section images, generating an aggregate model, calculating aggregate model parameters, carrying out format conversion and establishing an aggregate model database:
preparing a grinding test piece containing a target aggregate: manufacturing a test piece containing a monomer aggregate or a test piece containing a group aggregate according to the requirements established by the aggregate model database; the manufacturing process of the single aggregate-containing test piece is to wrap a single aggregate in a cementing material to manufacture a polishing test piece, the cementing material is a quick-hardening material convenient for polishing, the quick-hardening material comprises cement mortar and building gypsum, and the single aggregate-containing test piece is convenient for the convenient extraction of the aggregate and the reconstruction of a three-dimensional model; the manufacturing process of the test piece containing the group aggregate is the same as that of the common concrete test piece, and the aggregates containing the group aggregate test piece are adhered and need to be segmented and extracted;
(II) acquiring an aggregate continuous section image: fix the aggregate test piece of polishing on grinding equipment, the bistrique is polished the test piece and is carried out the successive layer dry grinding to the aggregate with the same interval of polishing along the Z axle direction of aggregate test piece of polishing, and the interval of polishing of section is controlled through the feed of bistrique, and the minimum interval of polishing is 0.01mm, and every layer of section is polished and is accomplished the back, uses digital camera to shoot the section, acquires aggregate continuous section image:
(III) calibrating and preprocessing continuous section images: the method comprises the steps of calibrating aggregate continuous section images by adopting a three-anchor point alignment method, wherein when an aggregate grinding test piece is ground in a layered mode, the photographed section images of each layer cannot be registered between layers, the images have errors due to the influence of various factors when the images are photographed, the images need to be registered, three control points close to a multiphase material test piece on a test piece clamping mechanism are selected as reference points, the three reference points are not on the same straight line, when the continuous section images are imaged, each section image comprises the three reference points, one standard image is selected as a reference image, the rest section images are images to be registered, the standard images are images with the angles, distances, directions and object positions in the continuous section images meeting the set requirements, the image registration is a process of seeking one-to-one mapping between the images to be registered and the standard images, and is a point corresponding to the same position in space in the two images is connected, three control points in the image to be registered and three reference points in the reference image are coincided, a spatial mapping relation is calculated according to the positions of the control points, geometric transformation is carried out on the image to be registered by using the spatial mapping relation, and in MATLAB, a cpselect function interactive tool is used for carrying out Z-axis direction registration on the continuous section image to obtain a registration result; the method comprises the steps of shooting anchor points which do not belong to section images of the porous concrete sample, uniformly cutting the registered continuous section images by utilizing Photoshop and image plus, cutting useless areas including the anchor points around the continuous section images, reserving the continuous section images of the concrete sample polishing area to be analyzed, carrying out gray-scale processing on the cut continuous section images and converting the cut continuous section images into binary images, wherein the white area is aggregate, the black area is gel material, and completing calibration and pretreatment of the continuous section images;
(IV) generating an aggregate model:
(1) in order to ensure that the proportion of image pixels to the real aggregate dimension is consistent, conversion calculation of the image pixels to the real aggregate dimension is carried out on a continuous sectional image, scanning resolutions of an X axis, a Y axis and a Z axis are input in the process of importing the continuous sectional image into Mimics three-dimensional reconstruction software, point coordinates (X ', Y') on the image correspond to corresponding point coordinates (X, Y) on a real slice, the scanning resolutions of the X axis and the Y axis are the image real dimension corresponding to each pixel, and the scanning resolution of the Z axis is the interlayer spacing d of the continuous sectional image;
(2) threshold segmentation, wherein white aggregate areas are extracted to generate aggregate masks, white areas are directly extracted from single aggregate samples, a group of aggregate samples need to judge whether the white areas in different two-dimensional images are from the same aggregate or not, the judgment process is realized through software programming, the centroid distance of any amount of aggregates is compared with the maximum radius of the two aggregates, the centroid distance of the two aggregates is smaller than the radius of any one aggregate to judge that the two aggregates are from the same aggregate, the white areas of the same aggregate are extracted by adopting a method of combining gray value range definition and manual drawing and are stored as the same mask, and the masks of the adhered parts of the different aggregates in the continuous sectional images are manually trimmed and segmented through mask editing operation;
(3) 3D computing the aggregate mask extracted from the single aggregate test piece in the step (2) to generate an initial aggregate surface grid data model, respectively performing 3D computing on different aggregate masks extracted from aggregate test pieces of a group, and simultaneously generating a plurality of initial aggregate surface grid data models, wherein the generated aggregate surface grid data model is a coarse primary model due to the influence of the interlayer spacing of the continuous section images of the aggregate test pieces on the Z-axis resolution;
(4) carrying out fairing, meshing and wrapping treatment on the primary aggregate surface mesh data model generated in the step (3) to obtain an accurate aggregate surface mesh data model, and filling the accurate aggregate surface mesh data model to obtain a real aggregate-based volume mesh data model;
and (V) calculating aggregate model parameters and carrying out format conversion: the geometric parameters directly extracted from the aggregate grid data model reconstructed by the Mimics software comprise model volume, surface area, long-axis distance and short-axis distance in the X-axis direction, the Y-axis direction and the Z-axis direction and mass barycentric coordinates determined based on a model coordinate frame, the aggregate model indirect shape parameters including aggregate shape factor, roundness, needle degree and flatness can be calculated through the directly extracted geometric parameters, further, the aggregate grid data model can be converted into a format which is convenient for storing an aggregate model database and is dxf, stl, inp and igs through corresponding interfaces of various application design software provided by the Mimics software, and the aggregate model format is convenient for putting and researching;
and (VI) establishing an aggregate model database: and (5) establishing aggregate grid data models with different sizes and types by utilizing the steps from (I) to (V), carrying out aggregate model parameter statistics and classification integration on the aggregate grid data models with different sizes and types, establishing an aggregate model database, and selecting and putting aggregates in subsequent numerical simulation research.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810059964.1A CN108280290B (en) | 2018-01-22 | 2018-01-22 | Concrete aggregate numerical model reconstruction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810059964.1A CN108280290B (en) | 2018-01-22 | 2018-01-22 | Concrete aggregate numerical model reconstruction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108280290A CN108280290A (en) | 2018-07-13 |
| CN108280290B true CN108280290B (en) | 2020-12-01 |
Family
ID=62804566
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810059964.1A Active CN108280290B (en) | 2018-01-22 | 2018-01-22 | Concrete aggregate numerical model reconstruction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108280290B (en) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109580624B (en) * | 2018-12-11 | 2021-02-09 | 广西利升石业有限公司 | Artificial stone aggregate simulation method |
| CN110542634B (en) * | 2018-12-15 | 2022-05-03 | 嘉思特华剑医疗器材(天津)有限公司 | Method for measuring through hole rate of 3D-printed titanium alloy bone trabecula test piece |
| CN110108561B (en) * | 2019-05-24 | 2020-06-12 | 四川大学 | Method for obtaining aggregate-mortar bonding strength change condition in concrete |
| CN110706352B (en) * | 2019-10-10 | 2023-03-10 | 重庆交通大学 | Concrete three-phase mesoscopic model construction and internal chloride ion erosion numerical simulation method based on polygonal random aggregate |
| CN111259944B (en) * | 2020-01-10 | 2022-04-15 | 河北工业大学 | Block stone shape classification method based on rapid PCA algorithm and K-means clustering algorithm |
| CN111241643B (en) * | 2020-02-11 | 2023-05-09 | 广东三维家信息科技有限公司 | Processing method and device of polygonal cabinet body and electronic equipment |
| CN112069704B (en) * | 2020-08-21 | 2024-04-26 | 武汉大学 | Method for constructing multi-degree-of-freedom lattice model of real concrete test piece |
| CN112070734B (en) * | 2020-08-31 | 2023-08-29 | 华南农业大学 | Seed germination soil emergence path model construction method, device, equipment and storage medium |
| CN112632814A (en) * | 2020-12-04 | 2021-04-09 | 西安理工大学 | Method for constructing three-dimensional model of recycled concrete based on ellipsoid random aggregate |
| CN113252103A (en) * | 2021-05-11 | 2021-08-13 | 安徽理工大学 | Method for calculating volume and mass of material pile based on MATLAB image recognition technology |
| CN114441566A (en) * | 2022-01-05 | 2022-05-06 | 哈尔滨工业大学 | Method for evaluating morphological characteristics of recycled aggregate for construction solid waste road |
| CN116342608B (en) * | 2023-05-30 | 2023-08-15 | 首都医科大学宣武医院 | Medical image-based stent adherence measurement method, device, equipment and medium |
| CN118094967A (en) * | 2024-04-28 | 2024-05-28 | 青岛理工大学 | Random generation method of two-dimensional elliptic concrete aggregate |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102297876A (en) * | 2011-05-20 | 2011-12-28 | 河海大学 | Construction method for thermal expansion coefficient prediction model of complete size grading large aggregate concrete in hardening phase |
| CN104493980A (en) * | 2014-12-03 | 2015-04-08 | 山东连银山环保建材有限公司 | Aggregate and concrete production equipment and aggregate and concrete production process |
| CN105599160A (en) * | 2016-03-28 | 2016-05-25 | 青岛理工大学 | Adjustable clamping fixing device with V-shaped structure |
| CN106053326A (en) * | 2016-07-19 | 2016-10-26 | 青岛理工大学 | Monitoring system for monitoring durability deterioration evolution of reinforced concrete structure |
| KR101736347B1 (en) * | 2015-11-27 | 2017-05-16 | 웰트웨이브 주식회사 | Measuring apparatus of concrete neutralization and concrete crack using smart device |
| CN107024492A (en) * | 2017-06-12 | 2017-08-08 | 太原理工大学 | A kind of method that use CT technical research high-temperature damages concrete defect develops |
-
2018
- 2018-01-22 CN CN201810059964.1A patent/CN108280290B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102297876A (en) * | 2011-05-20 | 2011-12-28 | 河海大学 | Construction method for thermal expansion coefficient prediction model of complete size grading large aggregate concrete in hardening phase |
| CN104493980A (en) * | 2014-12-03 | 2015-04-08 | 山东连银山环保建材有限公司 | Aggregate and concrete production equipment and aggregate and concrete production process |
| KR101736347B1 (en) * | 2015-11-27 | 2017-05-16 | 웰트웨이브 주식회사 | Measuring apparatus of concrete neutralization and concrete crack using smart device |
| CN105599160A (en) * | 2016-03-28 | 2016-05-25 | 青岛理工大学 | Adjustable clamping fixing device with V-shaped structure |
| CN106053326A (en) * | 2016-07-19 | 2016-10-26 | 青岛理工大学 | Monitoring system for monitoring durability deterioration evolution of reinforced concrete structure |
| CN107024492A (en) * | 2017-06-12 | 2017-08-08 | 太原理工大学 | A kind of method that use CT technical research high-temperature damages concrete defect develops |
Non-Patent Citations (3)
| Title |
|---|
| 一种快速生成三维混凝土骨料模型的混合实现方法;郭瑞奇;《土木建筑与环境工程》;20171015(第5期);第100-107页 * |
| 图像处理技术在混凝土微观结构测定中的应用研究;曹承伟;《工程建设》;20160131(第1期);第48-52页 * |
| 基于随机骨料模型的混凝土单轴破坏数值分析;刘斐骁;《居业》;20171220(第12期);第128-129页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108280290A (en) | 2018-07-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108280290B (en) | Concrete aggregate numerical model reconstruction method | |
| CN106600690B (en) | Complex building three-dimensional modeling method based on point cloud data | |
| CN109087396B (en) | Mesostructure reconstruction method based on concrete CT image pixel characteristics | |
| CN108364350B (en) | Three-dimensional reconstruction method for concrete mesoscopic structure model | |
| CN110223376B (en) | Three-dimensional particle reconstruction method based on single accumulated particle material image | |
| CN106096135B (en) | A kind of gradation particle library file numerical generation method based on AIMS equipment | |
| Kersten et al. | Potential of automatic 3D object reconstruction from multiple images for applications in architecture, cultural heritage and archaeology | |
| CN108956420A (en) | A kind of concrete pore feature extracting method and system | |
| CN111028355A (en) | Three-dimensional model reconstruction method for asphalt mixture | |
| CN112069704A (en) | Method for constructing multi-degree-of-freedom lattice model of real concrete specimen | |
| Jurda et al. | Three-dimensional documentation of Dolní Věstonice skeletal remains: can photogrammetry substitute laser scanning? | |
| CN106500626A (en) | A kind of mobile phone stereoscopic imaging method and three-dimensional imaging mobile phone | |
| CN109191571B (en) | Method for preparing mechanical test standard aggregate by applying 3D printing technology | |
| Thilakarathna et al. | Aggregate geometry generation method using a structured light 3D scanner, spherical harmonics–based geometry reconstruction, and placing algorithms for mesoscale modeling of concrete | |
| Anusree et al. | Characterization of sand particle morphology: state-of-the-art | |
| Afarani et al. | Print fidelity metrics for additive manufacturing of cement-based materials | |
| CN110599477B (en) | Method and system for detecting effective porosity of grass planting concrete | |
| CN114379092B (en) | Artificial core preparation method and system | |
| Górski et al. | Application of close-range photogrammetry in reverse engineering | |
| CN114283133A (en) | Asphalt mixture uniformity evaluation method based on nearest neighbor analysis | |
| Li et al. | A study on image-processing based identification of aspect ratio of coarse aggregate | |
| Seglins et al. | GEOVIZUALISATION OF STONE MATERIAL WEATHERING DATA FOR GEOARCHAEOLOGICAL STUDIES | |
| Xu et al. | Investigation of the coarse aggregate texture property using an improved Laplacian smoothing algorithm | |
| CN113255027B (en) | Efficient concrete three-dimensional aggregate generation and delivery method-three-dimensional residual space method | |
| CN113591321B (en) | Method for generating three-dimensional high-simulation structural model of mine earth-rock mixture |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230817 Address after: 257300 Li Que Zhen Huang Dong Cun, Guangrao County, Dongying City, Shandong Province Patentee after: Dongying Yusheng Building Materials Co.,Ltd. Address before: Room 1507, Building 3, China Railway Financial Intelligence Center, No. 59 Gongye South Road, Jinan Area, China (Shandong) Pilot Free Trade Zone, Jinan City, Shandong Province, 250000 Patentee before: Shandong Zhengyue Intellectual Property Co.,Ltd. Effective date of registration: 20230817 Address after: Room 1507, Building 3, China Railway Financial Intelligence Center, No. 59 Gongye South Road, Jinan Area, China (Shandong) Pilot Free Trade Zone, Jinan City, Shandong Province, 250000 Patentee after: Shandong Zhengyue Intellectual Property Co.,Ltd. Address before: 266520 No. 2 Changjiang Road, Qingdao economic and Technological Development Zone, Shandong Patentee before: Qindao University of Technology |
|
| TR01 | Transfer of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A numerical model reconstruction method for concrete aggregates Granted publication date: 20201201 Pledgee: Shandong Guangrao Rural Commercial Bank Co.,Ltd. Pledgor: Dongying Yusheng Building Materials Co.,Ltd. Registration number: Y2024980012289 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| CP03 | Change of name, title or address |
Address after: 257000 Li Que Zhen Huang Dong Cun, Guangrao County, Dongying City, Shandong Province Patentee after: Dongying Yusheng New Materials Co.,Ltd. Country or region after: China Address before: 257300 Li Que Zhen Huang Dong Cun, Guangrao County, Dongying City, Shandong Province Patentee before: Dongying Yusheng Building Materials Co.,Ltd. Country or region before: China |
|
| CP03 | Change of name, title or address |