CN108280290A - A kind of aggregate numerical model method for reconstructing - Google Patents

A kind of aggregate numerical model method for reconstructing Download PDF

Info

Publication number
CN108280290A
CN108280290A CN201810059964.1A CN201810059964A CN108280290A CN 108280290 A CN108280290 A CN 108280290A CN 201810059964 A CN201810059964 A CN 201810059964A CN 108280290 A CN108280290 A CN 108280290A
Authority
CN
China
Prior art keywords
aggregate
image
serial
test specimen
polishing
Prior art date
Application number
CN201810059964.1A
Other languages
Chinese (zh)
Other versions
CN108280290B (en
Inventor
范宏
赵铁军
张鹏
万小梅
刘树明
Original Assignee
青岛理工大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 青岛理工大学 filed Critical 青岛理工大学
Priority to CN201810059964.1A priority Critical patent/CN108280290B/en
Publication of CN108280290A publication Critical patent/CN108280290A/en
Application granted granted Critical
Publication of CN108280290B publication Critical patent/CN108280290B/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Abstract

The invention belongs to concrete numerical simulation study technical fields,It is related to a kind of aggregate numerical model method for reconstructing,It takes pictures to comprising monomer aggregate test specimen or carrying out polishing comprising group's aggregate test specimen,Obtain the serial section image of aggregate,Aggregate model is established based on serial section image,According to aggregate model foundation aggregate model database,Including group's aggregate test specimen can establish multiple aggregate models simultaneously,It can solve the problems, such as to be not easy to implement because of CT scan equipment costliness,And three-dimensional reconstruction is carried out based on true aggregate particle,True aggregate resemblance is simulated to the full extent,The aggregate model of foundation preserves in a grid formation,It creates conditions for the foundation of true aggregate model database,And it can be exchanged into row format with CAD software or FEA Finite Element Simulation Softwares,Convenient for the numerical simulation study of mechanical performance of concrete;The mode that its image obtains is simple, quick, and the image sequence interlamellar spacing of acquisition is small, Aligning degree is high, information content is complete.

Description

A kind of aggregate numerical model method for reconstructing

Technical field:

The invention belongs to concrete numerical simulation study technical fields, are related to a kind of aggregate numerical model reconstruction side Method can fast and accurately obtain the small spacing serial section image of heterogeneous material.

Background technology:

Concrete (referred to as concrete) refers to being led to the general designation of the cementing integral engineered composite material of aggregate by cementitious material The concrete often said is made of coarse aggregate, fine aggregate, cement-hydrate, unhydrated cement granules, hole and crackle etc.;With water (can contain additive and admixture) coordinates by a certain percentage, it is agitated obtained by cement concrete, also referred to as normal concrete, it is wide It is general to be applied to civil engineering.It is continuous with computer technology currently, the research for concrete is limited only to macro-level Develop, the research in terms of micro-mechanics studies of concrete is increasingly subject to the attention of people, can be ground after carrying out numerical simulation to concrete Study carefully aggregate in concrete, the damage -form of mortar and its boundary layer between the two, the premise using numerical method analysis is to build The finite element geometrical model of vertical concrete material microscopical structure, the performance of concrete and the shape of aggregate, grading etc. have prodigious Relationship, so, the generation of aggregate mould shapes just becomes an important topic of concrete numerical simulation.

Microscopical experiments, can under the conditions of theoretical model is reasonable and each phase material performance data of concrete is sufficiently accurate To replace part Experiment, and the influence of the objective limitation and human factor of experiment condition to its result can be avoided, for mixed The foundation of solidifying soil meso-mechanical model uses many it is assumed that these hypothesis have certain error:For example, concrete numerical simulation Critical issue be simulation to aggregate shape, in previous research, to two-dimentional aggregate multipurpose circular or ellipse Round simulation cobble or gravel aggregate simulate crushed stone aggregate with convex polygon, and mostly ovum is simulated with ball or elliposoidal in three-dimensional Stone or gravel aggregate simulate crushed stone aggregate with convex polyhedron;These simplified ways are feasible in method, but cannot be accurate Express the morphological feature of true aggregate particle, although two-dimensional analog calculate it is easy, its disadvantage be also it will be apparent that The spatial distribution and Grading feature that aggregate cannot really be restored, can not give expression to the destructive characteristics of concrete;It is three-dimensional Spherical in simulation, elliposoidal and the random of convex polyhedron launch the spatial distribution form for capableing of simulation concrete aggregate, cannot Completely instead of the shape of true aggregate, the Shortcomings in terms of the research of aggregate and mortar interface layer.It is more true in order to obtain Concrete structure and destructive characteristics, it is necessary to true aggregate is simulated using more efficient way.

The scanning to concrete cut-plane slices is realized using industry CT scanning technique, is then mixed with three-dimensional reconstruction software The reconstruction of solidifying soil model, obtains the concrete numerical model based on true shape, then analyzed on the basis of numerical model It can obtain preferable effect.But there is also shortcoming, CT scan images point in the application of concrete field for CT scan technology When resolution improves, density resolution declines, and causes the information content for including insufficient, simultaneously as CT scan technology is set by scanning Standby is affected, and easy ting produce noise and artifact interference causes the model rebuild inaccurate, and CT scan equipment price Costliness is unfavorable for extensive use of the CT scan technology in concrete field, and thus, CT scan technology is thin for true concrete Seeing structure three-dimensional reconstruction, there is also certain limitations.

A kind of aggregate structure feature based on the moulding method of section disclosed in Chinese patent 201510801544.2 is picked up Method is taken, greasy filth module is made using the concrete material for carrying out three-dimensional reconstruction, manufactured greasy filth module is cut, is swept It retouches, the greasy filth module scans image of acquisition is subjected to image procossing, and carry out the reconstruct of concrete particle structure three-dimensional spatial model, Obtained concrete particle structure three-dimensional spatial model is finally realized into aggregate into the random distribution of line direction and position Structure feature is picked up, and when carrying out test block cutting, cutting thickness is not easy to control, the excessive model for causing to rebuild of cutting spacing and reality Error between object is larger;A kind of method of CT images of concrete three-dimensional reconstruction disclosed in Chinese patent 201210015951.7, Tomoscan is carried out to concrete sample with CT machines first, secondly ray casting algorithm carries out Three-dimensional Gravity to concrete C T images It builds, obtains the true picture of aggregate, mortar, hole, finally CT images of concrete is carried out using visualization tool VTK three-dimensional vertical Weight structure can clearly distinguish aggregate, mortar and hole by reconstruction model, there are scan image by CT equipment interfere with And the problem of precision.

Based on the above reason, it is necessary to seek a kind of more accurate and easily aggregate numerical model method for reconstructing, The aggregate model database based on true aggregate shape is further established using numerical model, in the Numerical-Mode of concrete Aggregate is conveniently extracted from model database in quasi- research and is launched at random.

Invention content:

It is an object of the invention to overcome disadvantage of the existing technology, R & D design one kind being based on true aggregate The method that the serial section image of acquisition of successively polishing carries out Model Reconstruction further utilizes the aggregate model that the method is established to give birth to At aggregate model database.To increase the validity of concrete numerical simulation.

To achieve the goals above, the specific embodiment of aggregate numerical model method for reconstructing of the present invention Including make the polishing test specimen comprising target aggregate, obtain aggregate serial section image, calibration and pretreatment serial section image, Aggregate model is generated, aggregate model parameter is calculated and goes forward side by side and row format conversion and establish aggregate model database totally six steps:

(1) the polishing test specimen for including target aggregate is made:According to the requirement that aggregate model database is established, making includes Monomer aggregate test specimen includes group's aggregate test specimen;Including the manufacturing process of monomer aggregate test specimen is that single aggregate is wrapped in glue It is fabricated to polishing test specimen in gel material, cementitious material is for ease of the fast hard material of polishing, including cement mortar and building gypsum, It is easily extracted and the reconstruction of threedimensional model including monomer aggregate test specimen is convenient for aggregate;Include the making of group's aggregate test specimen Process is identical as the manufacturing process of ordinary concrete specimens, is adhered phenomenon including existing between the aggregate of group's aggregate test specimen, needs It is split extraction;

(2) aggregate serial section image is obtained:Aggregate polishing test specimen is fixed on equipment for grinding, bistrique is beaten along aggregate It grinds the Z-direction of test specimen successively to dry grind to aggregate polishing test specimen with identical polishing spacing, the polishing spacing of section passes through The amount of feeding of bistrique controls, and minimum polishing spacing is 0.01mm, after the completion of every layer of section polishing, using digital camera to section into Row shooting, obtains aggregate serial section image:

(3) calibrate and pre-process serial section image:Aggregate serial section image is carried out using three anchor point alignment methods Calibration, when due to aggregate polishing test specimen layering polishing, each layer of profile image of shooting cannot be guaranteed to be registrated in interlayer, due to figure As shooting when influenced by various factors, as the angle of video camera, distance and direction, object space movement and other because Element causes image error occur, needs to be registrated image, chooses in specimen holder mechanism close to the three of heterogeneous material test specimen A control point as datum mark, three datum marks not point-blank, when serial section image is imaged, each profile image Including these three datum marks, select a wherein width standard picture as benchmark image, remaining profile image is image subject to registration, mark Quasi- image is the image that angle, distance, direction and object space in serial section image meets sets requirement, image registration It is the process for seeking one-to-one mapping between image and standard picture subject to registration, is that will correspond to the same position in space in two images The point set connects, and three in image subject to registration control point should be overlapped with three datum marks in benchmark image, according to control The position of system point carries out geometric transformation to calculate the relationship of space reflection, then with the relationship of space reflection to image subject to registration, In MATLAB (matrix labotstory), the registration of Z-direction is carried out to serial section image using cpselect functions interactive tool, Obtain registration result;The anchor point for being not belonging to porous concrete specimens profile image is contained when shooting, utilizes Photoshop (figures As processing software), image proplus (image analysis software) carry out unified image to the serial section image after registration and cut It cuts out, crops the useless region for including anchor point around serial section image, retain the concrete sample polishing region for needing to analyze Serial section image, the serial section image gray processing after cutting out handles and is converted to bianry image, wherein white area For aggregate, black region is cementitious material, completes the calibration and pretreatment of serial section image;

It is panoramic picture based on serial section image, shoots and preserve the difference of the imaging device parameter of image, imaging is set The standby difference with porous concrete specimens section relative position and angle so that existence position deviation between each layer profile image, It needs to be standardized the serial section image of acquisition, to meet the needs of experimental analysis;

(4) aggregate model is generated:

(1) it is to ensure that image pixel (pixel) is consistent with the true ratio of aggregate scale, figure is carried out to serial section image As the conversion reckoning of pixel and true aggregate size, it is three-dimensional to import Mimics (medical image control system) in serial section image X-axis, the scanning resolution of Y-axis and Z axis are inputted during reconstruction software, point coordinates on image (x ', y ') correspondence is really cut The scanning resolution of the corresponding point coordinates (x, y) of on piece, x-axis and y-axis is the image full-size(d) corresponding to each pixel, z-axis Scanning resolution is the interlamellar spacing d of serial section image;

(2) Threshold segmentation, extracts the aggregate Area generation aggregate mask of white, and monomer aggregate test specimen directly extracts white area Domain, group's aggregate test specimen need first to judge whether white area derives from same aggregate in different two dimensional images, and deterministic process is logical Software programming realization is crossed, the centre of form distance of any amount aggregate is compared with two aggregate maximum radius, two aggregate centre of form distances Judge that two aggregates come from same aggregate less than any one aggregate radius, using intensity value ranges definition with it is manual retouch to take be combined Method extract the white area of same aggregate and save as same mask, portion is adhered to different aggregate in serial section image Divide mask to pass through mask edit operation and carries out manual pruning and segmentation;

(3) 3D calculating is carried out to the aggregate mask that the monomer aggregate test specimen in step (2) extracts, generates initial aggregate face Grid data model carries out 3D calculating respectively to the different aggregate masks that group's aggregate test specimen extracts, while generating multiple first Beginning aggregate surface grids data model, the influence due to aggregate test specimen serial section image interlamellar spacing to Z axis resolution ratio, the bone of generation Charge level grid data model is coarse primary mold;

(4) to step (3) generate primary aggregate surface grids data model carry out fairing, mesh generation and package processing with Accurate aggregate surface grids data model is obtained, accurate aggregate surface grids data model is filled, is obtained based on true The volume mesh data model of aggregate;

The principle that Mimics (medical image control system) three-dimensional reconstruction software inversely rebuilds particles of aggregates model is:It is based on The three-dimensional information for the unknown object that serial section image includes reconstructs the digital three-dimensional model of contour of object in a computer, Point P is in the coordinate value (x, y) of X-axis and Y-axis by the location determination on its horizontal section wherein on unknown object profile, and point P is in Z The coordinate value (z) of axis is by the spacing (d) cut and the relational expression of section sequence (n):Z=d × n is determined;

(5) calculate aggregate model parameter go forward side by side row format conversion:The aggregate volume mesh number directly rebuild from Mimics softwares Include model volume according to the geometric parameter extracted in model, surface area, the long axis distance of X-direction, Y direction and Z-direction, Mass cg coordinate short axle distance and determined based on model coordinate frame, can be carried out by the geometric parameter directly extracted The calculating of aggregate model indirect form parameter including aggregate form factor, circularity, needle degree and flatness further passes through Mimics softwares provide various the corresponding interface using design software can convert aggregate volume mesh data model to dxf, The format of stl, inp and igs preserved convenient for aggregate model database, the throwing of aggregate model format mutually turned convenient for aggregate model It puts and studies;

(6) aggregate model database is established:The aggregate body net of different sizes and type is established using step (1) to (five) Lattice data model carries out aggregate model parameter statistics to the aggregate volume mesh data model of different sizes and type and sorts out whole It closes, establishes aggregate model database, and carry out the selection and dispensing of aggregate in follow-up numerical simulation study.

The interlamellar spacing size for the aggregate serial section image that step (1) of the present invention is related to affects aggregate Model Reconstruction knot The interlamellar spacing of the resolution ratio of fruit in the Z-axis direction, aggregate serial section image is smaller, the external shape of true aggregate model State fitting effect is more accurate;Step (1) obtain aggregate serial section image high resolution, including information completely, can The multiple dimensioned numerical model of multiphase is rebuild, is convenient for the research of different levels, and assign the coagulation based on aggregate serial section image The larger flexibility of the three-dimensional reconstruction of native microscopical structure;The aggregate model database that step (6) is established has validity height, bone Material chooses convenient and efficient feature, and solid foundation is provided for aggregate numerical simulation study.

Compared with prior art, the present invention to carrying out polishing bat comprising monomer aggregate test specimen or comprising group's aggregate test specimen According to obtaining the serial section image of aggregate, aggregate model established based on serial section image, according to aggregate model foundation aggregate mould Type database can be solved because of CT scan equipment costliness not including group's aggregate test specimen can establish multiple aggregate models simultaneously Problem convenient to carry out, and three-dimensional reconstruction is carried out based on true aggregate particle, true aggregate is simulated to the full extent The aggregate model of resemblance, foundation preserves in a grid formation, creates conditions for the foundation of true aggregate model database, and And can be exchanged into row format with CAD software or FEA Finite Element Simulation Softwares, it is ground convenient for the numerical simulation of mechanical performance of concrete Study carefully;The mode that its image obtains is simple, quick, and image procurement cost is low, the image sequence interlamellar spacing of acquisition is small, Aligning degree is high, Information content is complete, and the aggregate model database of foundation has true and reliable characteristic, is convenient for application.

Description of the drawings:

Fig. 1 is the process flow diagram of the present invention.

Fig. 2 is the profile image before the binaryzation that 1 step (3) of the embodiment of the present invention is related to.

Fig. 3 is the profile image after the binaryzation that 1 step (3) of the embodiment of the present invention is related to.

Fig. 4 is the first established model figure of aggregate that 1 step (4) of the embodiment of the present invention is related to.

Fig. 5 is the aggregate surface grids illustraton of model that 1 step (4) of the embodiment of the present invention is related to.

Fig. 6 is the aggregate volume mesh illustraton of model that 1 step (4) of the embodiment of the present invention is related to.

Fig. 7 is the CAD aggregate illustratons of model after the conversion that 1 step (5) of the embodiment of the present invention is related to.

Fig. 8 is that the shape description that the embodiment of the present invention 2 is related to is the aggregate illustraton of model for being biased to flat class.

Fig. 9 is the aggregate illustraton of model that the shape description that the embodiment of the present invention 2 is related to is slightly flat class.

Figure 10 is the aggregate illustraton of model that the shape description that the embodiment of the present invention 2 is related to is slightly needle-shaped class.

Figure 11 is that the shape description that the embodiment of the present invention 2 is related to is the aggregate illustraton of model for being biased to needle-shaped class.

Specific implementation mode:

The invention will be further described by way of example and in conjunction with the accompanying drawings.

Embodiment 1:

The specific embodiment for the aggregate numerical model method for reconstructing that the present embodiment is related to includes making comprising mesh Mark the polishing test specimen of aggregate, obtain aggregate serial section image, calibration and pretreatment serial section image, generate aggregate model, Aggregate model parameter is calculated to go forward side by side row format conversion totally five steps:

(1) the polishing test specimen for including target aggregate is made:Cube body concrete of cast-in-place 100mm × 100mm × 100mm Cube concrete test specimen is placed in after standard curing room conserves 7 days the concrete sample for obtaining including group's aggregate by test block;

(2) aggregate serial section image is obtained:Aggregate polishing test specimen is fixed on equipment for grinding, bistrique is beaten along aggregate It grinds the Z-direction of test specimen successively to dry grind to aggregate polishing test specimen with the polishing spacing of 1mm, the polishing spacing of section passes through mill The amount of feeding control of head shoots section using digital camera after the completion of every layer of section polishing, and acquisition interlamellar spacing is 1mm Aggregate serial section image:

(3) calibrate and pre-process serial section image:Aggregate serial section image is carried out using three anchor point alignment methods Calibration, when due to aggregate polishing test specimen layering polishing, each layer of profile image of shooting cannot be guaranteed to be registrated in interlayer, due to figure As shooting when influenced by various factors, as the angle of video camera, distance and direction, object space movement and other because Element causes image error occur, needs to be registrated image, chooses in specimen holder mechanism close to the three of heterogeneous material test specimen A control point as datum mark, three datum marks not point-blank, when serial section image is imaged, each profile image Including these three datum marks, select a wherein width standard picture as benchmark image, remaining profile image is image subject to registration, mark Quasi- image is the image that angle, distance, direction and object space in serial section image meets sets requirement, image registration It is the process for seeking one-to-one mapping between image and standard picture subject to registration, is that will correspond to the same position in space in two images The point set connects, and three in image subject to registration control point should be overlapped with three datum marks in benchmark image, according to control The position of system point carries out geometric transformation to calculate the relationship of space reflection, then with the relationship of space reflection to image subject to registration, In MATLAB (matrix labotstory), the registration of Z-direction is carried out to serial section image using cpselect functions interactive tool, Obtain registration result;The anchor point for being not belonging to porous concrete specimens profile image is contained when shooting, utilizes Photoshop (figures As processing software), image proplus (image analysis software) carry out unified image to the serial section image after registration and cut It cuts out, crops the useless region for including anchor point around serial section image, retain the concrete sample polishing region for needing to analyze Serial section image, the serial section image gray processing after cutting out handles and is converted to bianry image, wherein white area For aggregate, black region is cementitious material, completes the calibration and pretreatment of serial section image;

(4) aggregate model is generated:

(1) it is to ensure that image pixel (pixel) is consistent with the true ratio of aggregate scale, figure is carried out to serial section image As the conversion reckoning of pixel and true aggregate size, it is three-dimensional to import Mimics (medical image control system) in serial section image X-axis, the scanning resolution of Y-axis and Z axis are inputted during reconstruction software, point coordinates on image (x ', y ') correspondence is really cut The corresponding point coordinates (x, y) of on piece, according to prediction equationWith X=0.25mm, y=0.25mm is calculated, result of calculation indicates the picture size that each pixel of image is included, x-axis and y-axis Scanning resolution be 0.25mm, z-axis scanning resolution is the spacing 1mm of profile image;

(2) Threshold segmentation, extracts the aggregate Area generation aggregate mask of white, and monomer aggregate test specimen directly extracts white area Domain, group's aggregate test specimen need first to judge whether white area derives from same aggregate in different two dimensional images, and deterministic process is logical Software programming realization is crossed, the centre of form distance of any amount aggregate is compared with two aggregate maximum radius, two aggregate centre of form distances Judge that two aggregates come from same aggregate less than any one aggregate radius, using intensity value ranges definition with it is manual retouch to take be combined Method extract the white area of same aggregate and save as same mask, portion is adhered to different aggregate in serial section image Divide mask to pass through mask edit operation and carries out manual pruning and segmentation;

(3) 3D calculating is carried out respectively to different masks, the aggregate zone boundary fitting under identical mask is become into matrix form Preliminary three-dimensional data model, be fitted generation preliminary three-dimensional data model be initial aggregate surface grids data model;

(4) to step (3) generate primary aggregate surface grids data model carry out fairing, mesh generation and package processing with Accurate aggregate surface grids data model is obtained, accurate aggregate surface grids data model is filled, is obtained based on true The volume mesh data model of aggregate;

(5) calculate aggregate model parameter go forward side by side row format conversion:The CAD area of computer aided provided by Mimics softwares Mapping software interface, dxf formats the answering in order to aggregate model database for converting under CAD aggregate volume mesh data model to With, or the FEA finite element analysis software interfaces that are provided by Mimics softwares, it converts aggregate volume mesh data model to limited Format needed for meta analysis carries out the analysis of aggregate mechanical property.

Embodiment 2:

The aggregate model that the present embodiment builds embodiment 1 carries out model X-direction, the short axle of Y direction and Z-direction The geometric parameter of distance, long axis distance, volume, surface area and the mass cg coordinate determined based on model coordinate frame is carried It takes, as a result such as table 1 (mm):

, aggregate form factor, the indirect form parameter of circularity, needle degree and flatness are calculated according to the geometric parameter of extraction, As a result such as table 2:

Shape description It is biased to flat class Slightly flat class Slightly needle-shaped class It is biased to needle-shaped class Form factor -0.70 0.70-1.00 1.00-1.40 1.40- Circularity -0.88 0.88-1.00 1.00-1.12 1.12- Needle degree 1.43-2.09 1.24-1.91 1.60-2.30 1.80-2.57 Flatness 0.13-0.45 0.36-0.58 0.45-0.98 0.45-0.88

, indirect form parameter integrate and sorts out the shape feature for judging aggregate;To aggregate model in Tables 1 and 2 Parameter carries out classification integration, and establishes true aggregate model database in this approach, is thrown at random for aggregate Concrete numerical model is put and established, the numerical simulation study of concrete is provided beneficial to help.

Claims (2)

1. a kind of aggregate numerical model method for reconstructing, it is characterised in that specific embodiment includes making comprising targeted bone The polishing test specimen of material obtains aggregate serial section image, calibration and pretreatment serial section image, generates aggregate model, calculates Aggregate model parameter, which is gone forward side by side, row format conversion and establishes aggregate model database totally six steps:
(1) the polishing test specimen for including target aggregate is made:According to the requirement that aggregate model database is established, it includes monomer to make Aggregate test specimen includes group's aggregate test specimen;Including the manufacturing process of monomer aggregate test specimen is that single aggregate is wrapped in gelling material It is fabricated to polishing test specimen in material, cementitious material is for ease of the fast hard material of polishing, including cement mortar and building gypsum, including Monomer aggregate test specimen is convenient for aggregate and easily extracts and the reconstruction of threedimensional model;Include the manufacturing process of group's aggregate test specimen It is identical as the manufacturing process of ordinary concrete specimens, it is adhered phenomenon including existing between the aggregate of group's aggregate test specimen, needs to carry out Segmentation extraction;
(2) aggregate serial section image is obtained:Aggregate polishing test specimen is fixed on equipment for grinding, bistrique is polished along aggregate and tried The Z-direction of part successively dry grinds to aggregate polishing test specimen with identical polishing spacing, and the polishing spacing of section passes through bistrique Amount of feeding control, minimum polishing spacing is that 0.01mm claps section using digital camera after the completion of the polishing of every layer of section It takes the photograph, obtains aggregate serial section image:
(3) calibrate and pre-process serial section image:School is carried out to aggregate serial section image using three anchor point alignment methods Standard, when due to aggregate polishing test specimen layering polishing, each layer of profile image of shooting cannot be guaranteed to be registrated in interlayer, due to image Influenced by various factors when shooting, as the angle of video camera, distance and direction, object space movement and other because Element causes image error occur, needs to be registrated image, chooses in specimen holder mechanism close to the three of heterogeneous material test specimen A control point as datum mark, three datum marks not point-blank, when serial section image is imaged, each profile image Including these three datum marks, select a wherein width standard picture as benchmark image, remaining profile image is image subject to registration, mark Quasi- image is the image that angle, distance, direction and object space in serial section image meets sets requirement, image registration It is the process for seeking one-to-one mapping between image and standard picture subject to registration, is that will correspond to the same position in space in two images The point set connects, and three in image subject to registration control point should be overlapped with three datum marks in benchmark image, according to control The position of system point carries out geometric transformation to calculate the relationship of space reflection, then with the relationship of space reflection to image subject to registration, In MATLAB, the registration of Z-direction is carried out to serial section image using cpselect functions interactive tool, obtains registration knot Fruit;The anchor point for being not belonging to porous concrete specimens profile image is contained when shooting, utilizes Photoshop, image proplus Unified image cutting-out is carried out to the serial section image after registration, is cropped useless including anchor point around serial section image Region retains the serial section image for the concrete sample polishing region for needing to analyze, the serial section image ash after cutting out Degreeization handles and is converted to bianry image, wherein white area is aggregate, and black region is cementitious material, completes serial section The calibration and pretreatment of image;
(4) aggregate model is generated:
(1) it is to ensure that the ratio of image pixel and true aggregate scale is consistent, to serial section image progress image pixel and very The conversion of real aggregate size calculates, X-axis, Y-axis are inputted during serial section image imports Mimics three-dimensional reconstruction softwares With the scanning resolution of Z axis, the corresponding point coordinates (x, y) on the corresponding true slice of point coordinates (x ', y ') on image, x-axis and y The scanning resolution of axis is the image full-size(d) corresponding to each pixel, and z-axis scanning resolution is the layer of serial section image Spacing d;
(2) Threshold segmentation extracts the aggregate Area generation aggregate mask of white, and monomer aggregate test specimen directly extracts white area, Group's aggregate test specimen needs first to judge whether white area derives from same aggregate in different two dimensional images, and deterministic process passes through soft Part programming is realized, is compared with two aggregate maximum radius to the centre of form distance of any amount aggregate, two aggregate centre of form distances are less than Any one aggregate radius judges that two aggregates come from same aggregate, is defined using intensity value ranges and retouches the side for taking and being combined with craft Method extracts the white area of same aggregate and saves as same mask, is covered to the part that is adhered of different aggregates in serial section image Plate carries out manual pruning and segmentation by mask edit operation;
(3) 3D calculating is carried out to the aggregate mask that the monomer aggregate test specimen in step (2) extracts, generates initial aggregate surface grids Data model carries out 3D calculating to the different aggregate masks that group's aggregate test specimen extracts, while generating multiple initial bones respectively Charge level grid data model, the influence due to aggregate test specimen serial section image interlamellar spacing to Z axis resolution ratio, the aggregate face of generation Grid data model is coarse primary mold;
(4) fairing, mesh generation and package processing are carried out to obtain to the primary aggregate surface grids data model that step (3) generates Accurate aggregate surface grids data model, is filled accurate aggregate surface grids data model, obtains being based on true aggregate Volume mesh data model;
(5) calculate aggregate model parameter go forward side by side row format conversion:The aggregate volume mesh data mould directly rebuild from Mimics softwares The geometric parameter extracted in type includes model volume, surface area, long axis distance, the short axle of X-direction, Y direction and Z-direction Distance and the mass cg coordinate determined based on model coordinate frame can be included by the geometric parameter directly extracted The calculating of aggregate form factor, circularity, the aggregate model indirect form parameter of needle degree and flatness further passes through Mimics softwares provide various the corresponding interface using design software can convert aggregate volume mesh data model to dxf, The format of stl, inp and igs preserved convenient for aggregate model database, the throwing of aggregate model format mutually turned convenient for aggregate model It puts and studies;
(6) aggregate model database is established:The aggregate volume mesh number of different sizes and type is established using step (1) to (five) According to model, the aggregate volume mesh data model of different sizes and type is carried out aggregate model parameter statistics and sorted out to integrate, is built Aggregate model database is found, and carries out the selection and dispensing of aggregate in follow-up numerical simulation study.
2. aggregate numerical model method for reconstructing according to claim 1, it is characterised in that the step (1) relates to And the interlamellar spacing size of aggregate serial section image affect the resolution ratio of aggregate Model Reconstruction result in the Z-axis direction, bone Expect that the interlamellar spacing of serial section image is smaller, the formalness fitting effect of true aggregate model is more accurate;Step (1) obtain aggregate serial section image high resolution, including information completely, the multiple dimensioned numerical value of multiphase can be rebuild Model is convenient for the research of different levels, and assigns the three-dimensional reconstruction of the Meso-level Structure of Concrete based on aggregate serial section image Larger flexibility;The aggregate model database that step (6) is established has validity high, and aggregate chooses convenient and efficient feature, Solid foundation is provided for aggregate numerical simulation study.
CN201810059964.1A 2018-01-22 2018-01-22 Concrete aggregate numerical model reconstruction method CN108280290B (en)

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 true CN108280290A (en) 2018-07-13
CN108280290B CN108280290B (en) 2020-12-01

Family

ID=62804566

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810059964.1A CN108280290B (en) 2018-01-22 2018-01-22 Concrete aggregate numerical model reconstruction method

Country Status (1)

Country Link
CN (1) CN108280290B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109580624A (en) * 2018-12-11 2019-04-05 广西利升石业有限公司 The method of artificial stone aggregate simulation
CN110108561A (en) * 2019-05-24 2019-08-09 四川大学 A kind of aggregate-mortar bonding strength situation of change method in acquisition concrete
CN109580624B (en) * 2018-12-11 2021-02-09 广西利升石业有限公司 Artificial stone aggregate simulation method

Citations (6)

* Cited by examiner, † Cited by third party
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
刘斐骁: "基于随机骨料模型的混凝土单轴破坏数值分析", 《居业》 *
曹承伟: "图像处理技术在混凝土微观结构测定中的应用研究", 《工程建设》 *
郭瑞奇: "一种快速生成三维混凝土骨料模型的混合实现方法", 《土木建筑与环境工程》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109580624A (en) * 2018-12-11 2019-04-05 广西利升石业有限公司 The method of artificial stone aggregate simulation
CN109580624B (en) * 2018-12-11 2021-02-09 广西利升石业有限公司 Artificial stone aggregate simulation method
CN110108561A (en) * 2019-05-24 2019-08-09 四川大学 A kind of aggregate-mortar bonding strength situation of change method in acquisition concrete

Also Published As

Publication number Publication date
CN108280290B (en) 2020-12-01

Similar Documents

Publication Publication Date Title
Fröhlich et al. Terrestrial laser scanning—new perspectives in 3D surveying
CN104766366B (en) A kind of method for building up of three-dimension virtual reality demonstration
Guidi et al. A multi-resolution methodology for the 3D modeling of large and complex archeological areas
Yue et al. Finite element modeling of geomaterials using digital image processing
CN106780619B (en) Human body size measuring method based on Kinect depth camera
Lerones et al. A practical approach to making accurate 3D layouts of interesting cultural heritage sites through digital models
Bernard et al. New trends in rapid product development
CN102855663B (en) Method for building CSG (Constructive Solid Geometry) model according to laser radar grid point cloud
Pieraccini et al. 3D digitizing of cultural heritage
CN105184250B (en) A kind of terrain classification method of electric power corridor airborne LiDAR point cloud data
Wang et al. Three-dimensional digital representation of granular material microstructure from X-ray tomography imaging
Ivanov et al. Computer stereo plotting for 3-D reconstruction of a maize canopy
Lee et al. Direct integration of reverse engineering and rapid prototyping
Zhou et al. Topology repair of solid models using skeletons
CN103644860B (en) Large space free curved face measurement method
CN1301480C (en) Computerized 3D visual color scan system and its scanning mode
CN101339670B (en) Computer auxiliary three-dimensional craniofacial rejuvenation method
CN104282040A (en) Finite element preprocessing method for reconstructing three-dimensional entity model
CN104266587B (en) Three-dimensional measurement system and method for obtaining actual 3D texture point cloud data
CN103489214A (en) Virtual reality occlusion handling method, based on virtual model pretreatment, in augmented reality system
RU2291488C2 (en) Method for stereological examination of objects structural organization
CN102760309B (en) Bituminous pavement test piece X-ray CT (computed tomography) image-based aggregate mesoscopic solid model reconstruction method
Zhang et al. An efficient approach to directly compute the exact Hausdorff distance for 3D point sets
CN107672180B (en) A kind of 3D printing accuracy checking method based on reverse Engineering Technology
Wu et al. Quality-driven poisson-guided autoscanning

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