CN108387495B - 一种多孔混凝土孔隙率计算和孔隙参数表征方法 - Google Patents

一种多孔混凝土孔隙率计算和孔隙参数表征方法 Download PDF

Info

Publication number
CN108387495B
CN108387495B CN201810059436.6A CN201810059436A CN108387495B CN 108387495 B CN108387495 B CN 108387495B CN 201810059436 A CN201810059436 A CN 201810059436A CN 108387495 B CN108387495 B CN 108387495B
Authority
CN
China
Prior art keywords
test piece
pore
porous concrete
images
image
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
Application number
CN201810059436.6A
Other languages
English (en)
Other versions
CN108387495A (zh
Inventor
范宏
赵铁军
张鹏
刘树明
王鹏刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong Zhengyue Intellectual Property Co.,Ltd.
Wanning Yongji Commercial Concrete Co ltd
Original Assignee
Qingdao University of Technology
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 Qingdao University of Technology filed Critical Qingdao University of Technology
Priority to CN201810059436.6A priority Critical patent/CN108387495B/zh
Publication of CN108387495A publication Critical patent/CN108387495A/zh
Application granted granted Critical
Publication of CN108387495B publication Critical patent/CN108387495B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/38—Registration of image sequences
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N2015/0846—Investigating permeability, pore-volume, or surface area of porous materials by use of radiation, e.g. transmitted or reflected light
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30108—Industrial image inspection
    • G06T2207/30132—Masonry; Concrete

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Geometry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • Dispersion Chemistry (AREA)
  • Computer Hardware Design (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Sampling And Sample Adjustment (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明属于混凝土孔隙参数分析技术领域，涉及一种多孔混凝土孔隙率计算和孔隙参数表征方法，采用优化设计的脆性材料剖面磨削箱室配合磨削刀具对多相材料试件进行固定间距的逐层干磨，拍摄多相材料试件的连续剖面图像，进而利用图像处理软件对多相材料试件连续剖面图像进行配准、标准化处理和预处理操作获取标准可用的多相材料试件连续剖面图像，在连续剖面图像获取时考虑了剖面图像在Z轴方向的对准，实现了较小层间距和较高校准度的连续剖面图像的获取，具备“模拟CT”图像的特点，获得的真彩色图像克服了CT扫描的灰度图像信息量丢失的缺点；其图像获取的方式简单、快捷，图像获取成本低，获取的图像信息量完整，便于应用推广。

Description

一种多孔混凝土孔隙率计算和孔隙参数表征方法
技术领域：
本发明属于混凝土孔隙参数分析技术领域，涉及一种多孔混凝土孔隙率计算和孔隙参数表征方法，基于三维重建的多孔混凝土模型计算多孔混凝土的孔隙率、表征多孔混凝土的孔隙参数。
背景技术：
混凝土(简称为砼)是指由胶凝材料将骨料胶结成整体的工程复合材料的统称，通常讲的混凝土一词是指用水泥作胶凝材料，砂、石作骨料；与水(可含外加剂和掺合料)按一定比例配合，经搅拌而得的水泥混凝土，也称普通混凝土，它广泛应用于土木工程。多孔混凝土作为多孔质轻的新型混凝土材料越来越受到人们的青睐，基于多孔混凝土的透气、透水、高散热等性能能够有效地缓解城市热岛效应，其被越来越多地将其应用于城市道路及河道护坡工程中。
目前，对于多孔混凝土的研究多集中于孔隙率与抗压强度、透水系数等宏观方面的研究，只有较少的研究涉及孔隙结构形态特征与空间分布特征的影响，关于孔隙的研究，研究者均集中于混凝土切片进行二维平面内孔隙面积、等效直径、轮廓周长的研究，没有定量地对孔隙的三维空间分布进行研究；多孔混凝土孔隙结构形态及空间分布特征影响着多孔混凝土的透水性能，然而对多孔混凝土孔隙结构形态及空间分布的研究，尚未有较好的评价方法和借鉴依据。
多孔混凝土中的孔隙按其连通情况可分为连通孔隙、半连通孔隙和封闭孔隙三种孔隙形态，三者之和即为总孔隙率，其中封闭孔隙因其孤立不连通的特性起不到排水效果，故称之为无效孔隙；连通孔隙有较好的透水排水性能，为有效孔隙；半连通孔隙是为一端开口一端封闭孔隙，虽然没有透水排水功能，但其具有一定的蓄水能力，对路面降水方面是有效的，也称之为有效孔隙。孔隙率按有效使用情况分为总孔隙率和有效孔隙率两种，多孔混凝土孔隙率的计算方法包括传统实验法和图像分析法：
1、传统实验法包括重量法和体积法两种，《透水混凝土及其应用技术》(宋中南，石云兴.透水混凝土及其应用技术，中国建筑工业出版社，2011年8月)中分别给出了重量法和体积法的计量多孔混凝土孔隙率的过程，其中重量法对孔隙率的测定使用的是电子天平，分别称量试件烘干后的重量和在水中的重量，两者之差即为试件因孔隙被水所填充而实际收到的浮力，假设试件无孔隙，用理论上所受到的浮力减去实际收到的浮力即可得孔隙率P的公式，其计算公式为：
Figure GDA0002359215980000021
式中，m1为试件在水中的重量(g)，m2为试件在烘箱中烘烤24小时后的重量(g)，v为试件的体积；体积法所用的计量设备为CoreLok真空密封仪，利用CoreLok真空密封仪测量压密试件的毛体积密度和最大表观密度，根据毛体积密度和最大表观密度计算出试件的连通孔隙率，具体操作步骤为：(1)、计算经真空密封的试件连同密封袋的密度ρ1，(2)、在水下剪开该试件包裹的密封袋，计算得到其水中密度ρ2，连通孔隙率是将ρ1和ρ2代入公式：
Figure GDA0002359215980000022
计算得到的；
2、图像分析法是利用数学模型结合图像处理的技术分析和提取得到孔隙率的方法，中国专利201510061021.9公开的一种基于图像处理的透水混凝土孔隙特征的检测方法包括把透水混凝土试件进行横向切割，以暴露出完整和平滑的切割平面为实验面，获取经过砖粉处理的透水混凝土实验面的正面照片，从获取的正面照片中截取出面积最大的形状规则的试件横截面图片，并计算其总像素值，对获得的图片黑白化处理，提取经过黑白处理后的图片的孔隙图像，计算各个孔隙的像素值，计算孔隙总像素值，根据得到的各个孔隙的像素值和孔隙总像素值来计算总孔隙率及内部孔径的分布情况，重复以上步骤，获得同一透水混凝土试件多个切割面上的总孔隙率以及孔径的分布情况，各个面总孔隙率再取平均值，得到透水混凝土试件总孔隙率，获取透水混凝土的孔隙特征；中国专利201310408977.2公开的一种基于CT技术的多孔混凝土孔结构表征方法，通过图像处理法分析CT扫描二维断层图像，得到多孔混凝土平面孔隙率、有效孔径和孔径分布等孔结构参数，具体包括以下步骤：首先利用CT扫描仪获取多孔混凝土二维断层图像，并通过计算机对图像进行黑白二值化处理得到清晰的孔隙分布图，然后利用图像处理软件计算出代表孔隙区域的面积，并把此区域转化为等面积的圆，最后统计分析孔径大于1mm的数据得到多孔混凝土平面孔隙率、有效孔径和孔径分布等孔结构参数；上述两种方法均是利用混凝土二维剖面图像进行孔隙率计算和孔径分布研究，其孔隙率计算的具体过程如下：(1)、横向切割多孔混凝土试件，以切割面为实验面，(2)、截取混凝土截面，并计算其总像素值A，(3)、将截面图片进行黑白化处理，使孔隙成为统一纯黑色，(4)、提取孔隙图像，并计算孔隙总像素值A1，(5)、根据像素值计算二维平面总孔隙率，试件截面总孔隙率表达式为：
Figure GDA0002359215980000031
(6)、计算各个面总孔隙率，取其平均值为试件总孔隙率，多孔混凝土试件总孔隙率表达式为：
Figure GDA0002359215980000032
其中n为试件截面二维图像的数量。
综合分析以上孔隙率计算方法可知，重量法和体积法主要基于实验原理进行的孔隙率测量，这两种方法忽略了封闭孔隙的影响，所测得的孔隙率为有效孔隙率(连通孔隙与半连通孔隙之和)；图像分析法是基于混凝土连续断层图像的图像分析原理进行测量，对每层剖面图像进行面孔率测量，取其平均值作为混凝土总孔隙率值，这种方法是基于平面图像的处理，所描述的是平面孔隙的分布，没有涉及多孔混凝土三维空间分布研究，对于孔隙结构的描述较为抽象，所计算孔隙率为平面总孔隙率，(连通孔隙，半连通孔隙与封闭孔隙之和)；以上方法都能够计算得到多孔混凝土的孔隙率，但是没有区分有效孔隙率与总孔隙率，对孔隙率的计算较为单一，并且没有涉及多孔混凝土的数值模拟研究领域，没有形象的说明多孔混凝土孔隙结构特点。
另外，孔隙迂曲度、连通度、孔喉比、喉道宽度、配位数等都是多孔混凝土孔隙的重要空间参数，其对多孔混凝土的强度、透水性能均有重要影响，但是基于多孔混凝土孔隙的复杂性和无规律性，利用现有技术对孔隙空间参数进行研究还存在较大困难。因此，针对多孔混凝土孔隙结构和透水性能研究存在的不足，孔隙空间形态结构的获取难度大导致孔隙结构研究进程缓慢的问题，研发设计一种基于数值模型的多孔混凝土孔隙率计算和孔隙参数的表征方法。
发明内容：
本发明的目的在于克服现有技术存在的多孔混凝土孔隙结构分析和孔隙率测量方法的局限性，提供一种基于数值模型的虚拟测量方式进行孔隙率的计算，并通过数值模型进行孔隙结构参数的表征。
为了实现上述目的，本发明涉及的多孔混凝土孔隙率计算和孔隙参数表征方法的具体工艺过程包括获取多孔混凝土连续剖面图像、标准化处理连续剖面图像、重建多孔混凝土模型、计算孔隙率和表征孔隙特征参数共五个步骤：
(一)获取多孔混凝土连续剖面图像：将磨削刀具与脆性材料剖面磨削箱室固定连接，把预制的或在多孔混凝土道路上取芯的多孔混凝土试件装设并固定在脆性材料剖面磨削箱室的试件夹持机构上，多孔混凝土试件跟随试件夹持箱门旋转，试件夹持箱门的最大旋转角度为90°，多孔混凝土试件截面中心距地面高度为105cm；将固定在三脚架上的数码相机的镜头高度调整为105cm使其与多孔混凝土试件的截面中心点在同一条直线上，将相机镜头与多孔混凝土试件截面的水平直线距离调整为30cm；开启磨削刀具沿多孔混凝土试件的Z轴方向对多孔混凝土试件进行打磨，通过磨削刀具的精密螺丝的旋转周数控制多孔混凝土试件剖面的打磨厚度和间距，精密螺丝的螺纹间距为1mm，精密螺丝旋转一周磨头的打磨进深间距为1mm；每层剖面打磨完成后进行截面拍照，每层剖面拍照后精密螺丝旋转一圈，打磨下一层剖面，以相同的拍照参数重复打磨—拍照过程，获取层间距为1mm的多孔混凝土试件连续剖面图像；
连续剖面图像的层间距越小，构建的三维模型越精确，通过精密螺丝的控制获得层间距为1mm连续剖面图像提高了三维模型重建的精准度，大大减少了图像层间差值运算量；
(二)标准化处理连续剖面图像：采用三锚点对准方法对连续剖面图像进行校准，由于混凝土试件分层打磨时，拍摄的每一层剖面图像不能保证在层间配准，由于图像拍摄时受到各种因素的影响，如摄像机的角度、距离和方向、物体位置的移动以及其他因素，导致图像出现误差，需要对图像进行配准，选取试件夹持机构上靠近多相材料试件的三个控制点作为基准点，三个基准点不在一条直线上，连续剖面图像成像时，每个剖面图像都包含这三个基准点，选择其中一幅标准图像作为基准图像，其余剖面图像为待配准图像，标准图像是连续剖面图像中的角度、距离、方向和物体位置均符合设定要求的图像，图像配准是寻求待配准图像和标准图像间一对一的映射的过程，是将两幅图像中对应于空间同一位置的点联系起来，待配准图像中的三个控制点与基准图像中的三个基准点应重合，根据控制点的位置来推算空间映射的关系，再用空间映射的关系对待配准图像进行几何变换，在MATLAB(矩阵实验室)中，使用cpselect函数交互工具对连续剖面图像进行Z轴方向的配准，获得配准结果；拍摄时包含了不属于多孔混凝土试件剖面图像的锚点，利用Photoshop(图像处理软件)、image proplus(图像分析软件)对配准后的连续剖面图像进行统一的图像剪裁，裁剪掉连续剖面图像周围包括锚点的无用区域，保留需要分析的混凝土试件打磨区域的连续剖面图像并将连续剖面图像调整为相同尺寸，完成连续剖面图像的标准化处理；
基于连续剖面图像为全景图像，拍摄和保存图像的成像设备参数的不同，成像设备与多孔混凝土试件剖面相对位置和角度的不同，使得各层剖面图像之间存在位置偏差，需要对获得的连续剖面图像进行标准化处理，以满足实验分析的需求；利用Matlab软件或image proplus软件对连续剖面图像进行的校准，能够增加图像的对准度，对连续剖面图像进行的图像增强操作，能够增加图像孔隙介质与固体介质的对比度，突显出孔隙轮廓的边界；
(三)重建多孔混凝土模型：为保证图像像素(pixel)与真实骨料尺度的比例一致，对连续剖面图像进行图像像素与真实骨料尺寸的转换推算，在连续剖面图像导入Mimics(医学影像控制系统)三维重建软件的过程中输入X轴、Y轴和Z轴的扫描分辨率，图像上的点坐标(x’，y’)对应真实切片上的相应点坐标(x，y)，然后对真彩色连续剖面图像进行灰度转化处理，得到的连续剖面图像中不同相的图像灰度具有不同的阈值范围；进一步对连续剖面图像进行阈值分割，分别选取多孔混凝土试件的连续剖面图像中的骨料、胶凝材料和孔隙三相并保存为各相蒙板，对骨料、胶凝材料和孔隙三相蒙板进行3D计算，分别将骨料、胶凝材料和孔隙三相蒙板重建为三维矩阵模型；
(四)计算孔隙率：利用mimics软件的虚拟测量功能对骨料、胶凝材料和孔隙三相模型进行体积计量，骨料模型体积为V1，胶凝材料模型体积为V2，总孔隙模型体积为V3，有效孔隙体积为V4，则透水混凝土总孔隙率计算公式为：
Figure GDA0002359215980000061
有效孔隙率计算公式为
Figure GDA0002359215980000062
其中多孔混凝土有效孔隙体积V4的测量方法是通过区域增长功能，分别选取多孔混凝土连续剖面图像中第一层与最后一层的孔隙区域，其余各层分别对与边界相交的孔隙进行区域增长，经区域增长后选取的孔隙区域为与边界相连的包括连通孔隙和半连通孔隙的有效孔隙区域，舍弃中间封闭的孔隙，将区域增长选取的有效孔隙区域保存为新的蒙板，对新的蒙板进行3D计算重建为有效孔隙的三维模型，再对有效孔隙三维模型进行虚拟测量计算得到有效孔隙模型体积V4；
(五)表征孔隙特征参数：通过有效孔隙三维模型能够直观的观察到孔隙空间分布特征，通过mimics软件孔隙分析模块能够对孔隙特征参数：平均孔径、孔径分布、孔喉比和孔隙迂曲度进行计量分析，其中孔隙迂曲度τ是渗流通道的实际长度与穿过渗流介质的宏观长度的比值，孔隙迂曲度τ的计算公式为：
Figure GDA0002359215980000063
其中，Lε为流体介质渗流通道的实际长度，L为渗流介质宏观长度，流体介质渗流通道的实际长度的计量过程是：选取有效孔隙三维模型从二维平面图像中生成的孔隙边界轮廓线，基于孔隙边界轮廓线拟合孔隙中心线，通过软件输出功能输出孔隙拟合中心线的长度，孔隙拟合中心线的长度即为流体介质渗流通道的实际长度Lε，流体介质点穿越介质单位距离时质点在孔道中运动轨迹的真实长度即为渗流介质宏观长度L。
本发明步骤(一)涉及的脆性材料剖面磨削箱室的主体结构包括箱体、抽屉、箱盖、试件夹持箱门、铰链机构、试件夹持机构、随动密封机构、观察箱门、合页、门锁机构、窗口和把手；内空式矩形结构的箱体的下部设置有矩形结构的抽屉，箱体的顶部设置有内凹式矩形结构的箱盖，箱体与箱盖铰接式连接，箱体的前侧面设置有矩形结构的试件夹持箱门，箱体与试件夹持箱门通过2-4个铰链机构连接，试件夹持箱门的内侧面设置有试件夹持机构，箱体的后侧面设置有随动密封机构，箱体与随动密封机构焊接式连接，箱体的左侧面设置有矩形结构的观察箱门，箱体与观察箱门通过2-4个合页连接；箱体与试件夹持箱门和箱体与观察箱门分别通过门锁机构实现开闭，箱盖的底面和观察箱门的中部分别开设有矩形板状结构的窗口，抽屉的面板上和箱盖的窗口上分别设置有圆弧形结构的把手，窗口与把手螺栓式连接；使用时，将脆性材料剖面磨削箱室与磨削机具连接，磨削机具的主动杆轴由随动密封机构进入箱体的内部，主动杆轴与试件夹持箱门垂直，打开试件夹持箱门，将待磨削的材料试件固定于试件夹持机构上并调整材料试件的截面位置，设定试件夹持机构的上夹持钳和下夹持钳的角点处为校准锚点，将安置有相机的三脚架固定在试件夹持箱门的正前方位置，调整相机的垂直高度使相机的中心高度与试件剖面的中心高度均为H1，按照设定的拍照距离H2调整相机与试件剖面之间的水平距离，相机镜头对准材料试件的正截面，通过门锁机构关闭试件夹持箱门，使箱体密封，开启磨削机具，按照设定的深度对材料试件进行磨削，磨削完成后，开启试件夹持箱门，试件夹持箱门带动试件夹持机构和材料试件跟随铰链机构的转动而旋转，试件夹持箱门旋转角度达到90°时，控制杆的低端的凸起与基座中孔道内的凹槽嵌合并保持稳定，支座与基座呈90°角的稳定状态，试件夹持箱门停止旋转并固定，打开相机对磨削好的材料试件进行拍照获取材料试件的剖面图像，拍摄完毕后再次关闭试件夹持箱门对材料试件进行下一深度的磨削，依此进行多次磨削和拍照，得到设定深度范围内材料试件的连续剖面图像，以锚点为基准点对连续剖面图像进行校准，以提高连续剖面图像的对准度，实现增加三维重建模型精准性的目的；试件夹持箱门每次打开后都处于同一固定位置，使相机镜头轴线与材料试件每层剖面的中心法线位于同一直线上，保证了材料试件连续剖面图像参数的一致性。
本发明涉及的铰链机构的主体结构包括支座、连接杆、控制杆、凸起、基座、孔道、凹槽和连接凸起；板状结构的支座的一端设置有内空式圆柱形结构的连接杆，支座的上表面中部横向设置有一端高另一端低的条状结构的控制杆，控制杆的一端的上底面和下底面分别设置有半球形伸缩式结构的凸起，板状结构的基座的中部开设有矩形结构的孔道，孔道的上表面和下表面分别设置有半球形结构的凹槽，基座的一端顶部和底部分别设置有半球形结构的连接凸起，连接杆的外壁与控制杆的外侧壁平齐，支座与试件夹持箱门螺栓式连接，基座与箱体螺栓式连接，连接杆嵌于连接凸起之间将支座与基座连接，凸起嵌合在凹槽中。
本发明涉及的试件夹持机构的主体结构包括连接片、扭簧、上夹持钳、上开口、连接件、下夹持钳和下开口；矩形片状结构的连接片与试件夹持箱门的上部螺栓式连接，连接片与试件夹持箱门之间设置有螺旋式结构的扭簧，扭簧的下端与块状结构的上夹持钳固定连接，扭簧压缩时，上夹持钳向上运动，扭簧伸展时，上夹持钳向下运动，上夹持钳的底部开设有倒V形结构的上开口，U形片状结构的连接件与试件夹持箱门螺栓式连接，连接件的U形口槽中卡接有块状结构的下夹持钳，下夹持钳在连接件的U形槽中上下移动，下夹持钳的顶部开设有V形结构的下开口。
本发明涉及的随动密封机构的主体结构包括底座板、上滑槽、下滑槽、第一层薄片、第一通孔、第二层薄片、第二通孔、第三层薄片和第三通孔；矩形板状结构的底座板的上部设置有U形槽状结构的上滑槽，底座板的下部设置有U形槽状结构的下滑槽，上滑槽与下滑槽之间设置有三层矩形板状结构的薄片，外层为第一层薄片，第一层薄片的中部开设有圆形结构的第一通孔，中间层为第二层薄片，第二层薄片的中部开设有胶囊形结构的第二通孔，内层为第三层薄片，第三层薄片的中部开设有胶囊形结构的第三通孔；底座板焊接在箱体的侧面，底座板的长度大于上滑槽和下滑槽的长度，上滑槽的长度与下滑槽的长度相等，上滑槽和下滑槽的长度大于第三层薄片的长度，第三层薄片的长度大于第二层薄片的长度，第二层薄片的长度大于第一层薄片的长度，第一层薄片、第二层薄片和第三层薄片在上滑槽和下滑槽之间移动，第一通孔的直径小于第二通孔的长度，第二通孔的长度小于第三通孔的长度。
本发明涉及的门锁机构的主体结构包括底板、卡板、控制板、基板和挡板；圆形板状结构的底板的上表面与板状结构的卡板转动式连接，卡板能够围绕卡板与底板的连接点旋转360°，卡板的上表面与板状结构的控制板固定连接，矩形板状结构的基板的上表面设置有L形板状结构的挡板，底板分别与试件夹持箱门和观察箱门螺栓式连接，基板与箱体螺栓式连接，旋转控制板，控制板带动卡板旋转旋进或旋出挡板之间实现试件夹持箱门和观察箱门的闭合和开启。
本发明涉及的箱体、箱盖和合页的材质均为不锈钢；抽屉用于收集材料试件磨削过程中产生的粉末碎屑，便于集中处理；试件夹持箱门在开合过程中带动试件夹持机构旋转，试件夹持箱门的转动角度范围为0°-90°；铰链机构是旋转角度可控的钢铰链；试件夹持机构能够稳固的夹持不同规格尺寸的材料试件；随动密封机构随着磨削刀具主动杆轴的往复运动发生相对移动，在主动杆轴运动过程中，保证箱体的密封性，并且有效防止粉尘和噪声污染；观察箱门便于观察材料试件的磨削过程；门锁机构能够将试件夹持箱门和观察箱门分别与箱体紧密闭合，并能轻松开启试件夹持箱门和观察箱门；窗口的材质为玻璃，在磨削材料试件的过程中，便于肉眼观察箱体内部的状况；把手便于抽屉和箱盖的开合操作。
本发明与现有技术相比，采用优化设计的脆性材料剖面磨削箱室配合磨削刀具对多孔混凝土试件进行固定间距的逐层干磨，拍摄多孔混凝土试件的连续剖面图像，进而利用图像处理软件对多孔混凝土试件连续剖面图像进行标准化处理获取标准可用的多孔混凝土试件连续剖面图像，在连续剖面图像获取时考虑了剖面图像在Z轴方向的对准，实现了较小层间距和较高校准度的连续剖面图像的获取，进而对连续剖面图像进行多孔混凝土试件模型的三维重建，基于多孔混凝土试件三维模型进行总孔隙率和有效孔隙率的计算以及孔隙参数的表征；其图像获取的方式简单、快捷，图像获取成本低，获取的图像信息量完整，能够直观的观察多孔混凝土孔隙结构的分布状况，便于应用推广。
附图说明：
图1为本发明的工艺流程框图。
图2为本发明步骤(一)涉及的脆性材料剖面磨削箱室的主体结构原理示意图。
图3为本发明步骤(一)涉及的脆性材料剖面磨削箱室的外观结构原理示意图。
图4为本发明步骤(一)涉及的铰链机构的结构原理示意图。
图5为本发明步骤(一)涉及的支座的结构原理示意图。
图6为本发明步骤(一)涉及的基座的结构原理示意图。
图7为本发明步骤(一)涉及的试件夹持机构的结构原理示意图。
图8为本发明步骤(一)涉及的随动密封机构的结构原理示意图。
图9为本发明步骤(一)涉及的三层薄片的主体结构分散原理示意图。
图10为本发明步骤(一)涉及的门锁机构的结构原理示意图。
图11为本发明步骤(一)涉及的拍照状态示意图。
图12为本发明发明内容步骤(三)涉及的骨料三维模型图。
图13为本发明发明内容步骤(三)涉及的胶凝材料三维模型图。
图14为本发明发明内容步骤(三)涉及的孔隙三维模型图。
图15为本发明发明内容步骤(四)涉及的有效孔隙三维模型图。
图16为本发明发明内容步骤(五)涉及的有效孔隙拟合中心线示意图。
图17为本发明实施例步骤(四)涉及的重量法、平面图像法、虚拟测量法计算的孔隙率对比线图。
图18为本发明实施例步骤(五)涉及的孔隙率15％的多孔混凝土三维模型图。
图19为本发明实施例步骤(五)涉及的孔隙率20％的多孔混凝土三维模型图。
图20为本发明实施例步骤(五)涉及的孔隙率25％的多孔混凝土三维模型图。
图21为本发明实施例步骤(五)涉及的孔隙率15％的多孔混凝土三维模型的孔隙等效直径分布饼图。
图22为本发明实施例步骤(五)涉及的孔隙率20％的多孔混凝土三维模型的孔隙等效直径分布饼图。
图23为本发明实施例步骤(五)涉及的孔隙率25％的多孔混凝土三维模型的孔隙等效直径分布饼图。
图24为本发明实施例步骤(五)涉及的孔隙率与连通孔隙平均迂曲度关系示意图。
具体实施方式：
下面通过实施例并结合附图对本发明作进一步说明。
实施例1：
本实施例涉及的多孔混凝土孔隙率计算和孔隙参数表征方法的具体工艺过程包括获取多孔混凝土连续剖面图像、标准化处理连续剖面图像、重建多孔混凝土模型、计算孔隙率和表征孔隙特征参数共五个步骤：
(一)获取多孔混凝土连续剖面图像：按照表1多孔混凝土配合比制备孔隙率分别为15％、20％和25％的多孔混凝土试件，将磨削刀具与脆性材料剖面磨削箱室固定连接，分别把养护好的多孔混凝土试件装设并固定在脆性材料剖面磨削箱室的试件夹持机构6上，多孔混凝土试件跟随试件夹持箱门4旋转，试件夹持箱门4的最大旋转角度为90°，多孔混凝土试件截面中心距地面高度为105cm；将固定在三脚架上的数码相机的镜头高度调整为105cm使其与多孔混凝土试件的截面中心点在同一条直线上，将相机镜头与多孔混凝土试件截面的水平直线距离调整为30cm；开启磨削刀具沿多孔混凝土试件的Z轴方向对多孔混凝土试件进行打磨，通过磨削刀具的精密螺丝的旋转周数控制多孔混凝土试件剖面的打磨厚度和间距，精密螺丝的螺纹间距为1mm，精密螺丝旋转一周磨头的打磨进深间距为1mm；每层剖面打磨完成后进行截面拍照，每层剖面拍照后精密螺丝旋转一圈，打磨下一层剖面，以相同的拍照参数重复打磨—拍照过程80次，获取层间距为1mm的80张多孔混凝土试件连续剖面图像；
表1多孔混凝土配合比
Figure GDA0002359215980000111
Figure GDA0002359215980000121
(二)标准化处理连续剖面图像：采用三锚点对准方法对连续剖面图像进行校准，由于混凝土试件分层打磨时，拍摄的每一层剖面图像不能保证在层间配准，由于图像拍摄时受到各种因素的影响，如摄像机的角度、距离和方向、物体位置的移动以及其他因素，导致图像出现误差，需要对图像进行配准，选取试件夹持机构6上靠近多相材料试件的三个控制点作为基准点，三个基准点不在一条直线上，连续剖面图像成像时，每个剖面图像都包含这三个基准点，选择其中一幅标准图像作为基准图像，其余剖面图像为待配准图像，标准图像是连续剖面图像中的角度、距离、方向和物体位置均符合设定要求的图像，图像配准是寻求待配准图像和标准图像间一对一的映射的过程，是将两幅图像中对应于空间同一位置的点联系起来，待配准图像中的三个控制点与基准图像中的三个基准点应重合，根据控制点的位置来推算空间映射的关系，再用空间映射的关系对待配准图像进行几何变换，在MATLAB(矩阵实验室)中，使用cpselect函数交互工具对连续剖面图像进行Z轴方向的配准，获得配准结果；拍摄时包含了不属于多孔混凝土试件剖面图像的锚点，利用Photoshop(图像处理软件)、image proplus(图像分析软件)对配准后的连续剖面图像进行统一的图像剪裁，裁剪掉连续剖面图像周围包括锚点的无用区域，保留需要分析的混凝土试件打磨区域的连续剖面图像并将连续剖面图像调整为相同尺寸，使每一像素点对应的图像距离为0.25mm，完成连续剖面图像的标准化处理；
(三)重建多孔混凝土模型：为保证图像像素(pixel)与真实骨料尺度的比例一致，对连续剖面图像进行图像像素与真实骨料尺寸的转换推算，在连续剖面图像导入Mimics(医学影像控制系统)三维重建软件的过程中输入X轴、Y轴和Z轴的扫描分辨率，图像上的点坐标(x’，y’)对应真实切片上的相应点坐标(x，y)，根据推算公式
Figure GDA0002359215980000135
和
Figure GDA0002359215980000131
计算得到x＝0.25mm，y＝0.25mm，计算结果表示图像每一像素所包含的图像尺寸，x轴和y轴的扫描分辨率均为0.25mm，z轴扫描分辨率为剖面图像的间距1mm；然后对真彩色连续剖面图像进行灰度转化处理，得到的连续剖面图像中不同相的图像灰度具有不同的阈值范围；进一步对连续剖面图像进行阈值分割，选取多孔混凝土试件的连续剖面图像中的孔隙并保存为孔隙蒙板，其余部分作为固体介质保存为固体介质蒙板，对孔隙蒙板和固体介质蒙板进行3D计算，分别将固体介质蒙板和孔隙蒙板重建为多孔混凝土和孔隙的三维矩阵模型；
(四)计算孔隙率：利用mimics软件的虚拟测量功能分别对孔隙率为20％的多孔混凝土模型进行体积计量，骨料模型体积为V1，胶凝材料模型体积为V2，总孔隙模型体积为V3，有效孔隙体积为V4，则透水混凝土总孔隙率计算公式为：
Figure GDA0002359215980000132
有效孔隙率计算公式为
Figure GDA0002359215980000133
其中多孔混凝土有效孔隙体积V4的测量方法是通过区域增长功能，分别选取多孔混凝土连续剖面图像中第一层与最后一层的孔隙区域，其余各层分别对与边界相交的孔隙进行区域增长，经区域增长后选取的孔隙区域为与边界相连的包括连通孔隙和半连通孔隙的有效孔隙区域，舍弃中间封闭的孔隙，将区域增长选取的有效孔隙区域保存为新的蒙板，对新的蒙板进行3D计算重建为有效孔隙的三维模型，再对有效孔隙三维模型进行虚拟测量计算得到有效孔隙模型体积V4；
多孔混凝土模型的各相体积(mm3)测量结果见表2：
骨料体积V1 胶凝材料体积V2 总孔隙体积V3 有效孔隙体积V4
171344.68 181344.36 69479.74 62778.65
分别采用称重法和平面图像法对孔隙率为15％、20％和25％的多孔混凝土试件进行多孔混凝土孔隙率的计算，并将计算得到的多孔混凝土孔隙率与本方法涉及的虚拟测量法计算的多孔混凝土孔隙率进行对比分析，结果见表3：
Figure GDA0002359215980000134
从上表得知基于多孔混凝土三维模型计算的总孔隙率比重量法计算的孔隙率偏大，与平面图像法计算的总孔隙率的误差在5％以内，基于多孔混凝土三维模型计算的有效孔隙率与重量法计算的孔隙率的误差在5％以内，以上误差是多孔混凝土连续剖面图像获取和阈值分割过程中灰度值范围选取导致的，表明基于多孔混凝土三维模型计算孔隙率的方法是有效的；
(五)表征孔隙特征参数：通过有效孔隙三维模型能够直观的观察到孔隙空间分布特征，通过mimics软件孔隙分析模块能够对孔隙特征参数：平均孔径、孔径分布、孔喉比、孔隙迂曲度和连通孔隙平均迂曲度进行计量分析，其中孔隙迂曲度τ是渗流通道的实际长度与穿过渗流介质的宏观长度的比值，孔隙迂曲度τ的计算公式为：
Figure GDA0002359215980000141
其中，Lε为流体介质渗流通道的实际长度，L为渗流介质宏观长度，流体介质渗流通道的实际长度的计量过程是：选取有效孔隙三维模型从二维平面图像中生成的孔隙边界轮廓线，基于孔隙边界轮廓线拟合孔隙中心线，通过软件输出功能输出孔隙拟合中心线的长度，孔隙拟合中心线的长度即为流体介质渗流通道的实际长度Lε，流体介质点穿越介质单位距离时质点在孔道中运动轨迹的真实长度即为渗流介质宏观长度L；连通孔隙平均迂曲度的计算公式为：
Figure GDA0002359215980000142
其中，n为连通孔隙拟合曲线的条数，随着多孔混凝土孔隙率的增大，连通孔隙平均迂曲度减小，水流路径变短，所以，孔隙率大的多孔混凝土的透水系数随着多孔混凝土孔隙率的增大而增加；对三维重建的孔隙率为15％、20％和25％的真实多孔混凝土模型分别进行空间孔隙等效直径的提取与分析，结果显示，随着多孔混凝土数值模型孔隙率的增大，试件孔隙的等效直径主要分布范围由0-9mm增加为0-11mm，且孔隙不同等效直径所占比重也发生了变化，孔隙率为15％的试件孔隙的等效直径主要分布在1.5-4.5mm，占比为57.8％，孔隙率为20％的试件孔隙的等效直径主要分布在3-6mm，占比为56.9％，孔隙率为25％的试件孔隙的等效直径主要分布范围为4.5-9mm，占比为56.9％，4.5-10.5mm范围内的等效直径孔隙的数量明显增加，表明三维模型建立的多孔混凝土试件能够真实地反映多孔混凝土孔隙的分布特征。
本实施例步骤(一)涉及的脆性材料剖面磨削箱室的主体结构包括箱体1、抽屉2、箱盖3、试件夹持箱门4、铰链机构5、试件夹持机构6、随动密封机构7、观察箱门8、合页9、门锁机构10、窗口11和把手12；内空式矩形结构的箱体1的下部设置有矩形结构的抽屉2，箱体1的顶部设置有内凹式矩形结构的箱盖3，箱体1与箱盖3铰接式连接，箱体1的前侧面设置有矩形结构的试件夹持箱门4，箱体1与试件夹持箱门4通过2-4个铰链机构5连接，试件夹持箱门4的内侧面设置有试件夹持机构6，箱体1的后侧面设置有随动密封机构7，箱体1与随动密封机构7焊接式连接，箱体1的左侧面设置有矩形结构的观察箱门8，箱体1与观察箱门8通过2-4个合页9连接；箱体1与试件夹持箱门4和箱体1与观察箱门8分别通过门锁机构10实现开闭，箱盖3的底面和观察箱门8的中部分别开设有矩形板状结构的窗口11，抽屉2的面板上和箱盖3的窗口11上分别设置有圆弧形结构的把手12，窗口11与把手12螺栓式连接；使用时，将脆性材料剖面磨削箱室与磨削机具连接，磨削机具的主动杆轴由随动密封机构7进入箱体1的内部，主动杆轴与试件夹持箱门4垂直，打开试件夹持箱门4，将待磨削的材料试件固定于试件夹持机构6上并调整材料试件的截面位置，设定试件夹持机构6的上夹持钳63和下夹持钳66的角点处为校准锚点，将安置有相机的三脚架固定在试件夹持箱门4的正前方位置，调整相机的垂直高度使相机的中心高度与试件剖面的中心高度均为H1，按照设定的拍照距离H2调整相机与试件剖面之间的水平距离，相机镜头对准材料试件的正截面，通过门锁机构10关闭试件夹持箱门4，使箱体1密封，开启磨削机具，按照设定的深度对材料试件进行磨削，磨削完成后，开启试件夹持箱门4，试件夹持箱门4带动试件夹持机构6和材料试件跟随铰链机构5的转动而旋转，试件夹持箱门4旋转角度达到90°时，控制杆53的低端的凸起54与基座55中孔道56内的凹槽57嵌合并保持稳定，支座51与基座55呈90°角的稳定状态，试件夹持箱门4停止旋转并固定，打开相机对磨削好的材料试件进行拍照获取材料试件的剖面图像，拍摄完毕后再次关闭试件夹持箱门4对材料试件进行下一深度的磨削，依此进行多次磨削和拍照，得到设定深度范围内材料试件的连续剖面图像，以锚点为基准点对连续剖面图像进行校准，以提高连续剖面图像的对准度，实现增加三维重建模型精准性的目的；试件夹持箱门4每次打开后都处于同一固定位置，使相机镜头轴线与材料试件每层剖面的中心法线位于同一直线上，保证了材料试件连续剖面图像参数的一致性。
本实施例涉及的铰链机构5的主体结构包括支座51、连接杆52、控制杆53、凸起54、基座55、孔道56、凹槽57和连接凸起58；板状结构的支座51的一端设置有内空式圆柱形结构的连接杆52，支座51的上表面中部横向设置有一端高另一端低的条状结构的控制杆53，控制杆53的一端的上底面和下底面分别设置有半球形伸缩式结构的凸起54，板状结构的基座55的中部开设有矩形结构的孔道56，孔道56的上表面和下表面分别设置有半球形结构的凹槽57，基座55的一端顶部和底部分别设置有半球形结构的连接凸起58，连接杆52的外壁与控制杆53的外侧壁平齐，支座51与试件夹持箱门4螺栓式连接，基座55与箱体1螺栓式连接，连接杆52嵌于连接凸起58之间将支座51与基座55连接，凸起54嵌合在凹槽57中。
本实施例涉及的试件夹持机构6的主体结构包括连接片61、扭簧62、上夹持钳63、上开口64、连接件65、下夹持钳66和下开口67；矩形片状结构的连接片61与试件夹持箱门4的上部螺栓式连接，连接片61与试件夹持箱门4之间设置有螺旋式结构的扭簧62，扭簧62的下端与块状结构的上夹持钳63固定连接，扭簧62压缩时，上夹持钳63向上运动，扭簧62伸展时，上夹持钳63向下运动，上夹持钳63的底部开设有倒V形结构的上开口64，U形片状结构的连接件65与试件夹持箱门4螺栓式连接，连接件65的U形口槽中卡接有块状结构的下夹持钳66，下夹持钳66在连接件65的U形槽中上下移动，下夹持钳66的顶部开设有V形结构的下开口67。
本实施例涉及的随动密封机构7的主体结构包括底座板71、上滑槽72、下滑槽73、第一层薄片74、第一通孔75、第二层薄片76、第二通孔77、第三层薄片78和第三通孔79；矩形板状结构的底座板71的上部设置有U形槽状结构的上滑槽72，底座板71的下部设置有U形槽状结构的下滑槽73，上滑槽72与下滑槽73之间设置有三层矩形板状结构的薄片，外层为第一层薄片74，第一层薄片74的中部开设有圆形结构的第一通孔75，中间层为第二层薄片76，第二层薄片76的中部开设有胶囊形结构的第二通孔77，内层为第三层薄片78，第三层薄片78的中部开设有胶囊形结构的第三通孔79；底座板71焊接在箱体1的侧面，底座板71的长度大于上滑槽72和下滑槽73的长度，上滑槽72的长度与下滑槽73的长度相等，上滑槽72和下滑槽73的长度大于第三层薄片78的长度，第三层薄片78的长度大于第二层薄片76的长度，第二层薄片76的长度大于第一层薄片74的长度，第一层薄片74、第二层薄片76和第三层薄片78在上滑槽72和下滑槽73之间移动，第一通孔75的直径小于第二通孔77的长度，第二通孔77的长度小于第三通孔79的长度。
本实施例涉及的门锁机构10的主体结构包括底板101、卡板102、控制板103、基板104和挡板105；圆形板状结构的底板101的上表面与板状结构的卡板102转动式连接，卡板102能够围绕卡板102与底板101的连接点旋转360°，卡板102的上表面与板状结构的控制板103固定连接，矩形板状结构的基板104的上表面设置有L形板状结构的挡板105，底板101分别与试件夹持箱门4和观察箱门8螺栓式连接，基板104与箱体1螺栓式连接，旋转控制板103，控制板103带动卡板102旋转旋进或旋出挡板105之间实现试件夹持箱门4和观察箱门8的闭合和开启。
本实施例涉及的箱体1、箱盖3和合页9的材质均为不锈钢；抽屉2用于收集材料试件磨削过程中产生的粉末碎屑，便于集中处理；试件夹持箱门4在开合过程中带动试件夹持机构6旋转，试件夹持箱门4的转动角度范围为0°-90°；铰链机构5是旋转角度可控的钢铰链；试件夹持机构6能够稳固的夹持不同规格尺寸的材料试件；随动密封机构7随着磨削刀具主动杆轴的往复运动发生相对移动，在主动杆轴运动过程中，保证箱体1的密封性，并且有效防止粉尘和噪声污染；观察箱门8便于观察材料试件的磨削过程；门锁机构10能够将试件夹持箱门4和观察箱门8分别与箱体1紧密闭合，并能轻松开启试件夹持箱门4和观察箱门8；窗口11的材质为玻璃，在磨削材料试件的过程中，便于肉眼观察箱体1内部的状况；把手12便于抽屉2和箱盖3的开合操作。

Claims (7)

1.一种多孔混凝土孔隙率计算和孔隙参数表征方法，其特征在于借助脆性材料剖面磨削箱室实施，具体工艺过程包括获取多孔混凝土连续剖面图像、标准化处理连续剖面图像、重建多孔混凝土模型、计算孔隙率和表征孔隙参数共五个步骤：
(一)获取多孔混凝土连续剖面图像：将磨削刀具与脆性材料剖面磨削箱室固定连接，把预制的或在多孔混凝土道路上取芯的多孔混凝土试件装设并固定在脆性材料剖面磨削箱室的试件夹持机构上，多孔混凝土试件跟随试件夹持箱门旋转，试件夹持箱门的最大旋转角度为90°，多孔混凝土试件截面中心距地面高度为105cm；将固定在三脚架上的相机的镜头高度调整为105cm使其与多孔混凝土试件的截面中心点在同一条直线上，将相机镜头与多孔混凝土试件截面的水平直线距离调整为30cm；开启磨削刀具沿多孔混凝土试件的Z轴方向对多孔混凝土试件进行打磨，通过磨削刀具的精密螺丝的旋转周数控制多孔混凝土试件剖面的打磨厚度和间距，精密螺丝的螺纹间距为1mm，精密螺丝旋转一周磨头的打磨进深间距为1mm；每层剖面打磨完成后进行截面拍照，每层剖面拍照后精密螺丝旋转一圈，打磨下一层剖面，以相同的拍照参数重复打磨—拍照过程，获取层间距为1mm的多孔混凝土试件连续剖面图像；
(二)标准化处理连续剖面图像：采用三锚点对准方法对连续剖面图像进行校准，选取试件夹持机构上靠近混凝土试件的三个控制点作为基准点，三个基准点不在一条直线上，连续剖面图像成像时，每个剖面图像都包含这三个基准点，选择其中一幅标准图像作为基准图像，其余剖面图像为待配准图像，标准图像是连续剖面图像中的角度、距离、方向和物体位置均符合设定要求的图像，图像配准是寻求待配准图像和标准图像间一对一的映射的过程，是将两幅图像中对应于空间同一位置的点联系起来，待配准图像中的三个控制点与基准图像中的三个基准点应重合，根据控制点的位置来推算空间映射的关系，再用空间映射的关系对待配准图像进行几何变换，在MATLAB中，使用cpselect函数交互工具对连续剖面图像进行Z轴方向的配准，获得配准结果；拍摄时包含了不属于多孔混凝土试件剖面图像的锚点，利用Photoshop、image proplus对配准后的连续剖面图像进行统一的图像剪裁，裁剪掉连续剖面图像周围包括锚点的无用区域，保留需要分析的混凝土试件打磨区域的连续剖面图像并将连续剖面图像调整为相同尺寸，完成连续剖面图像的标准化处理；
(三)重建多孔混凝土模型：为保证图像像素与真实骨料尺度的比例一致，对连续剖面图像进行图像像素与真实骨料尺寸的转换推算，在连续剖面图像导入Mimics三维重建软件的过程中输入X轴、Y轴和Z轴的扫描分辨率，图像上的点坐标(x’，y’)对应真实切片上的相应点坐标(x，y)，然后对真彩色连续剖面图像进行灰度转化处理，得到的连续剖面图像中不同相的图像灰度具有不同的阈值范围；进一步对连续剖面图像进行阈值分割，分别选取多孔混凝土试件的连续剖面图像中的骨料、胶凝材料和孔隙三相并保存为各相蒙板，对骨料、胶凝材料和孔隙三相蒙板进行3D计算，分别将骨料、胶凝材料和孔隙三相蒙板重建为三维矩阵模型；
(四)计算孔隙率：利用mimics软件的虚拟测量功能对骨料、胶凝材料和孔隙三相模型进行体积计量，骨料模型体积为V1，胶凝材料模型体积为V2，总孔隙模型体积为V3，有效孔隙体积为V4，则透水混凝土总孔隙率计算公式为：
Figure FDA0002316545550000021
有效孔隙率计算公式为
Figure FDA0002316545550000022
其中多孔混凝土有效孔隙体积V4的测量方法是通过区域增长功能，分别选取多孔混凝土连续剖面图像中第一层与最后一层的孔隙区域，其余各层分别对与边界相交的孔隙进行区域增长，经区域增长后选取的孔隙区域为与边界相连的包括连通孔隙和半连通孔隙的有效孔隙区域，舍弃中间封闭的孔隙，将区域增长选取的有效孔隙区域保存为新的蒙板，对新的蒙板进行3D计算重建为有效孔隙三维模型，再对有效孔隙三维模型进行虚拟测量计算得到有效孔隙模型体积V4；
(五)表征孔隙参数：通过有效孔隙三维模型能够直观的观察到孔隙空间分布特征，通过mimics软件孔隙分析模块能够对孔隙参数：平均孔径、孔径分布、孔喉比和孔隙迂曲度进行计量分析，其中孔隙迂曲度τ是渗流通道的实际长度与穿过渗流介质的宏观长度的比值，孔隙迂曲度τ的计算公式为：
Figure FDA0002316545550000031
其中，Lε为流体介质渗流通道的实际长度，L为渗流介质宏观长度，流体介质渗流通道的实际长度的计量过程是：选取有效孔隙三维模型从二维平面图像中生成的孔隙边界轮廓线，基于孔隙边界轮廓线拟合孔隙中心线，通过软件输出功能输出孔隙拟合中心线的长度，孔隙拟合中心线的长度即为流体介质渗流通道的实际长度Lε，流体介质点穿越介质单位距离时质点在孔道中运动轨迹的真实长度即为渗流介质宏观长度L。
2.根据权利要求1所述的多孔混凝土孔隙率计算和孔隙参数表征方法，其特征在于所述步骤(一)涉及的脆性材料剖面磨削箱室的主体结构包括箱体、抽屉、箱盖、试件夹持箱门、铰链机构、试件夹持机构、随动密封机构、观察箱门、合页、门锁机构、窗口和把手；内空式矩形结构的箱体的下部设置有矩形结构的抽屉，箱体的顶部设置有内凹式矩形结构的箱盖，箱体与箱盖铰接式连接，箱体的前侧面设置有矩形结构的试件夹持箱门，箱体与试件夹持箱门通过2-4个铰链机构连接，试件夹持箱门的内侧面设置有试件夹持机构，箱体的后侧面设置有随动密封机构，箱体与随动密封机构焊接式连接，箱体的左侧面设置有矩形结构的观察箱门，箱体与观察箱门通过2-4个合页连接；箱体与试件夹持箱门和箱体与观察箱门分别通过门锁机构实现开闭，箱盖的底面和观察箱门的中部分别开设有矩形板状结构的窗口，抽屉的面板上和箱盖的窗口上分别设置有圆弧形结构的把手，窗口与把手螺栓式连接；使用时，将脆性材料剖面磨削箱室与磨削机具连接，磨削机具的主动杆轴由随动密封机构进入箱体的内部，主动杆轴与试件夹持箱门垂直，打开试件夹持箱门，将待磨削的材料试件固定于试件夹持机构上并调整材料试件的截面位置，设定试件夹持机构的上夹持钳和下夹持钳的角点处为校准锚点，将安置有相机的三脚架固定在试件夹持箱门的正前方位置，调整相机的垂直高度使相机的中心高度与试件剖面的中心高度均为H1，按照设定的拍照距离H2调整相机与试件剖面之间的水平距离，相机镜头对准材料试件的正截面，通过门锁机构关闭试件夹持箱门，使箱体密封，开启磨削机具，按照设定的深度对材料试件进行磨削，磨削完成后，开启试件夹持箱门，试件夹持箱门带动试件夹持机构和材料试件跟随铰链机构的转动而旋转，试件夹持箱门旋转角度达到90°时，控制杆的低端的凸起与基座中孔道内的凹槽嵌合并保持稳定，支座与基座呈90°角的稳定状态，试件夹持箱门停止旋转并固定，打开相机对磨削好的材料试件进行拍照获取材料试件的剖面图像，拍摄完毕后再次关闭试件夹持箱门对材料试件进行下一深度的磨削，依此进行多次磨削和拍照，得到设定深度范围内材料试件的连续剖面图像，以锚点为基准点对连续剖面图像进行校准，以提高连续剖面图像的对准度，实现增加三维重建模型精准性的目的；试件夹持箱门每次打开后都处于同一固定位置，使相机镜头轴线与材料试件每层剖面的中心法线位于同一直线上，保证了材料试件连续剖面图像参数的一致性。
3.根据权利要求2所述的多孔混凝土孔隙率计算和孔隙参数表征方法，其特征在于所述铰链机构的主体结构包括支座、连接杆、控制杆、凸起、基座、孔道、凹槽和连接凸起；板状结构的支座的一端设置有内空式圆柱形结构的连接杆，支座的上表面中部横向设置有一端高另一端低的条状结构的控制杆，控制杆的一端的上底面和下底面分别设置有半球形伸缩式结构的凸起，板状结构的基座的中部开设有矩形结构的孔道，孔道的上表面和下表面分别设置有半球形结构的凹槽，基座的一端顶部和底部分别设置有半球形结构的连接凸起，连接杆的外壁与控制杆的外侧壁平齐，支座与试件夹持箱门螺栓式连接，基座与箱体螺栓式连接，连接杆嵌于连接凸起之间将支座与基座连接，凸起嵌合在凹槽中。
4.根据权利要求1所述的多孔混凝土孔隙率计算和孔隙参数表征方法，其特征在于所述试件夹持机构的主体结构包括连接片、扭簧、上夹持钳、上开口、连接件、下夹持钳和下开口；矩形片状结构的连接片与试件夹持箱门的上部螺栓式连接，连接片与试件夹持箱门之间设置有螺旋式结构的扭簧，扭簧的下端与块状结构的上夹持钳固定连接，扭簧压缩时，上夹持钳向上运动，扭簧伸展时，上夹持钳向下运动，上夹持钳的底部开设有倒V形结构的上开口，U形片状结构的连接件与试件夹持箱门螺栓式连接，连接件的U形口槽中卡接有块状结构的下夹持钳，下夹持钳在连接件的U形槽中上下移动，下夹持钳的顶部开设有V形结构的下开口。
5.根据权利要求2所述的多孔混凝土孔隙率计算和孔隙参数表征方法，其特征在于所述随动密封机构的主体结构包括底座板、上滑槽、下滑槽、第一层薄片、第一通孔、第二层薄片、第二通孔、第三层薄片和第三通孔；矩形板状结构的底座板的上部设置有U形槽状结构的上滑槽，底座板的下部设置有U形槽状结构的下滑槽，上滑槽与下滑槽之间设置有三层矩形板状结构的薄片，外层为第一层薄片，第一层薄片的中部开设有圆形结构的第一通孔，中间层为第二层薄片，第二层薄片的中部开设有胶囊形结构的第二通孔，内层为第三层薄片，第三层薄片的中部开设有胶囊形结构的第三通孔；底座板焊接在箱体的侧面，底座板的长度大于上滑槽和下滑槽的长度，上滑槽的长度与下滑槽的长度相等，上滑槽和下滑槽的长度大于第三层薄片的长度，第三层薄片的长度大于第二层薄片的长度，第二层薄片的长度大于第一层薄片的长度，第一层薄片、第二层薄片和第三层薄片在上滑槽和下滑槽之间移动，第一通孔的直径小于第二通孔的长度，第二通孔的长度小于第三通孔的长度。
6.根据权利要求2所述的多孔混凝土孔隙率计算和孔隙参数表征方法，其特征在于所述门锁机构的主体结构包括底板、卡板、控制板、基板和挡板；圆形板状结构的底板的上表面与板状结构的卡板转动式连接，卡板能够围绕卡板与底板的连接点旋转360°，卡板的上表面与板状结构的控制板固定连接，矩形板状结构的基板的上表面设置有L形板状结构的挡板，底板分别与试件夹持箱门和观察箱门螺栓式连接，基板与箱体螺栓式连接，旋转控制板，控制板带动卡板旋转旋进或旋出挡板之间实现试件夹持箱门和观察箱门的闭合和开启。
7.根据权利要求2所述的多孔混凝土孔隙率计算和孔隙参数表征方法，其特征在于所述箱体、箱盖和合页的材质均为不锈钢；抽屉用于收集材料试件磨削过程中产生的粉末碎屑，便于集中处理；试件夹持箱门在开合过程中带动试件夹持机构旋转，试件夹持箱门的转动角度范围为0°-90°；铰链机构是旋转角度可控的钢铰链；试件夹持机构能够稳固的夹持不同规格尺寸的材料试件；随动密封机构随着磨削刀具主动杆轴的往复运动发生相对移动，在主动杆轴运动过程中，保证箱体的密封性，并且有效防止粉尘和噪声污染；观察箱门便于观察材料试件的磨削过程；门锁机构能够将试件夹持箱门和观察箱门分别与箱体紧密闭合，并能轻松开启试件夹持箱门和观察箱门；窗口的材质为玻璃，在磨削材料试件的过程中，便于肉眼观察箱体内部的状况；把手便于抽屉和箱盖的开合操作。
CN201810059436.6A 2018-01-22 2018-01-22 一种多孔混凝土孔隙率计算和孔隙参数表征方法 Active CN108387495B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810059436.6A CN108387495B (zh) 2018-01-22 2018-01-22 一种多孔混凝土孔隙率计算和孔隙参数表征方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810059436.6A CN108387495B (zh) 2018-01-22 2018-01-22 一种多孔混凝土孔隙率计算和孔隙参数表征方法

Publications (2)

Publication Number Publication Date
CN108387495A CN108387495A (zh) 2018-08-10
CN108387495B true CN108387495B (zh) 2020-03-31

Family

ID=63077103

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810059436.6A Active CN108387495B (zh) 2018-01-22 2018-01-22 一种多孔混凝土孔隙率计算和孔隙参数表征方法

Country Status (1)

Country Link
CN (1) CN108387495B (zh)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109001236A (zh) * 2018-09-03 2018-12-14 东南大学 一种基于纳米ct技术表征水泥浆体中孔隙率分布的方法
CN109374624B (zh) * 2018-12-12 2020-11-06 山东大学 透水路面堵塞检测方法及系统
CN109543350B (zh) * 2018-12-20 2023-03-24 哈尔滨工业大学（深圳） 基于细观模型重构的透水混凝土性能预测方法及系统
CN109993786B (zh) * 2019-03-08 2021-05-18 中国石油大学(北京) 迂曲度获取方法、装置、设备以及存储介质
CN109781607A (zh) * 2019-03-25 2019-05-21 南通烟滤嘴有限责任公司 烟支疏松度的测量装置及其测量方法
CN109883928A (zh) * 2019-03-29 2019-06-14 上海市城市建设设计研究总院（集团）有限公司 基于ct扫描和三维重构的水泥混凝土渗透能力模拟方法
CN110838117B (zh) * 2019-11-14 2022-05-13 中国科学院武汉岩土力学研究所 一种基于孔壁图像的岩石面孔隙度识别方法
CN112881253B (zh) * 2019-11-29 2022-11-25 华北理工大学 一种确定岩石铸体薄片图片孔隙区域的方法、面孔率计算方法及其应用
CN111060440B (zh) * 2020-01-13 2022-11-08 湖南工程学院 用于泡沫混凝土抗压强度的检测方法
CN111721925B (zh) * 2020-06-11 2022-06-10 广西交科集团有限公司 一种透水混凝土抗压强度与透水系数均匀性评价方法
CN111982743B (zh) * 2020-08-31 2022-06-24 长春工程学院 一种基于物联网的火山岩鉴定方法、系统、终端及介质
CN114199725A (zh) * 2020-09-17 2022-03-18 长沙理工大学 一种新旧沥青间扩散融合程度的测试与分析方法
CN113191017B (zh) * 2021-05-21 2022-12-27 成都数之联科技股份有限公司 虚拟量测模型构建方法、虚拟量测方法、系统、电子设备及计算机可读存储介质
CN113405965B (zh) * 2021-06-08 2022-06-28 浙江广天构件集团股份有限公司 一种水泥基材料颗粒堆积体系孔隙连通性分析方法
CN113763351B (zh) * 2021-09-06 2022-09-13 深圳大学 岩土体孔隙率测量方法、装置及智能终端
CN116223333A (zh) * 2023-01-06 2023-06-06 河北农业大学 一种土壤孔隙结构信息自动检测系统
CN116612114B (zh) * 2023-07-17 2023-10-20 安徽省交通规划设计研究总院股份有限公司 图像处理视域下的透水混凝土胶石比确定方法
CN116858751B (zh) * 2023-09-05 2023-11-17 安徽省交通规划设计研究总院股份有限公司 基于代表车型声音处理的透水混凝土路面孔隙率测算方法
CN117867938B (zh) * 2024-03-13 2024-05-28 安徽省交通规划设计研究总院股份有限公司 一种图像处理透水混凝土施工稳定性控制方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9002088B2 (en) * 2012-09-07 2015-04-07 The Boeing Company Method and apparatus for creating nondestructive inspection porosity standards
EA201500703A1 (ru) * 2012-12-25 2015-10-30 Общество С Ограниченной Ответственностью "Интровижн Ресерч Энд Девелопмент" Способ получения характеристической трехмерной модели образца пористого материала для исследования свойств проницаемости
CN105372168B (zh) * 2015-12-08 2018-03-23 清华大学 一种混凝土孔隙三维结构的构建方法
CN105806765A (zh) * 2016-04-13 2016-07-27 南京大学（苏州）高新技术研究院 一种显微ct扫描土体空间孔隙结构的精细化表征方法

Also Published As

Publication number Publication date
CN108387495A (zh) 2018-08-10

Similar Documents

Publication Publication Date Title
CN108387495B (zh) 一种多孔混凝土孔隙率计算和孔隙参数表征方法
CN108364350B (zh) 一种混凝土细观结构模型三维重建方法
CN108387524B (zh) 一种多相材料连续剖面图像获取方法
Rattanasak et al. Pore structure of cement/pozzolan composites by X-ray microtomography
CN106644637B (zh) 基于ct扫描与3d重构的透水混凝土试件制作方法
CN106827545B (zh) 基于ct扫描透明透水混凝土试件3d打印方法
Du Plessis et al. Standard method for microCT-based additive manufacturing quality control 2: Density measurement
Lux et al. 3D characterization of wood based fibrous materials: an application
CN109883928A (zh) 基于ct扫描和三维重构的水泥混凝土渗透能力模拟方法
CN106353357A (zh) 一种渗流作用下砂土介质细观结构变化的监测装置及方法
Wang et al. Construction of dynamic three-dimensional microstructure for the hydration of cement using 3D image registration
CN108663287B (zh) 一种利用ct图像精确计算煤岩密度的方法
Mac et al. 3D imaging techniques for characterising microcracks in cement-based materials
JP2003116825A (ja) レンダリングされたボリュームの画像に対する変換関数を非管理下で生成するための方法及びシステム
Wevers et al. X-ray computed tomography for nondestructive testing
Lin et al. Cone beam x-ray microtomography—a new facility for three-dimensional analysis of multiphase materials
CN201007646Y (zh) 一种液体辅助断层扫描三维形状测量装置
Veerle et al. Multi-disciplinary characterisation of a sandstone surface crust
CN109297881A (zh) 一种粗粒土渗透试验边壁处理层厚度的确定方法
CN109781605B (zh) 一种基于透明土非饱和流动可视化实验方法
CN207717475U (zh) 一种脆性材料剖面磨削箱室
CN115639131A (zh) 不同压力下散装粮堆孔隙表征方法
CN113640324B (zh) 一种结合断层扫描成像与增强造影的局部孔隙率表征方法
Bresson et al. High‐resolution bulk density images, using calibrated X‐ray radiography of impregnated soil slices
Werner et al. Quantitative image analysis of masonry mortar microstructure

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
TR01 Transfer of patent right

Effective date of registration: 20230811

Address after: 571500 Dong'aopo Third Titanium Selection Plant, Dong'ao Town, Wanning City, Hainan Province

Patentee after: Wanning Yongji Commercial Concrete 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: 20230811

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: QINGDAO TECHNOLOGICAL University