CN113405965B - 一种水泥基材料颗粒堆积体系孔隙连通性分析方法 - Google Patents
一种水泥基材料颗粒堆积体系孔隙连通性分析方法 Download PDFInfo
- Publication number
- CN113405965B CN113405965B CN202110637103.9A CN202110637103A CN113405965B CN 113405965 B CN113405965 B CN 113405965B CN 202110637103 A CN202110637103 A CN 202110637103A CN 113405965 B CN113405965 B CN 113405965B
- Authority
- CN
- China
- Prior art keywords
- pore
- pores
- tetrahedral
- tetrahedron
- connectivity
- 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
- 239000011148 porous material Substances 0.000 title claims abstract description 346
- 239000002245 particle Substances 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims abstract description 60
- 239000000463 material Substances 0.000 title claims abstract description 21
- 239000004568 cement Substances 0.000 title claims abstract description 14
- 239000012798 spherical particle Substances 0.000 claims abstract description 59
- 230000008569 process Effects 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 24
- 238000009825 accumulation Methods 0.000 claims description 35
- 239000011159 matrix material Substances 0.000 claims description 23
- 239000012071 phase Substances 0.000 claims description 21
- 239000007790 solid phase Substances 0.000 claims description 19
- 238000004891 communication Methods 0.000 claims description 16
- 230000000737 periodic effect Effects 0.000 claims description 11
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 9
- 238000012856 packing Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000000903 blocking effect Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 abstract description 3
- 230000005540 biological transmission Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 239000004567 concrete Substances 0.000 description 4
- 230000003628 erosive effect Effects 0.000 description 4
- 230000001788 irregular Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000001179 sorption measurement Methods 0.000 description 2
- 239000011800 void material Substances 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003763 carbonization Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011496 digital image analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000010257 thawing Methods 0.000 description 1
- 230000001550 time effect Effects 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种水泥基材料颗粒堆积体系孔隙连通性分析方法。目前测试孔隙的手段仅针对连通孔和部分半连通孔隙进行测试，无法探测非连通孔隙，无法准确评估材料的孔隙连通性。本发明所述球形颗粒堆积体系孔隙连通性分析方法为首先建立球体颗粒堆积模型，其次对球心点集执行Delaunay四面体剖分以把堆积模型划分为由多个四面体组成的集合后，逐一获取集合中四面体的孔隙体积以及判断每个四面体孔隙的连通性，最后进行球形颗粒堆积模型中连通孔隙率的统计计算过程。
Description
技术领域
本发明涉及一种水泥基材料颗粒堆积体系孔隙连通性分析方法，属于多孔材料孔径分布分析技术领域。
背景技术
水泥基复合材料具有多孔的性质，孔对材料的力学性能及耐久性有着重要的影响，其中耐久性如冻融、碳化、钢筋锈蚀等多为侵蚀介质通过孔隙结构的传输进入材料内部，逐步导致内部结构劣化、性能丧失。侵蚀介质向水泥基多相多孔材料中的传输速率除受环境温度影响外，主要取决于多孔材料自身的开口连通孔隙率。准确评估多孔材料的孔隙连通性是预测多孔材料抵抗侵蚀介质传输能力的关键，也是后续预测水泥基材料耐久性的宝贵依据，具有十分重要的工程实际意义。如今，测量多孔材料孔隙结构常用的实验方法主要有压汞法和吸附法。由于受试验技术和方法原理的限制，压汞法和吸附法所能测试的孔隙仅包括连通孔和部分半连通孔隙，无法探测非连通孔隙，故而无法准确评估材料的孔隙连通性。
得益于迅速发展的数字建模技术，依赖于计算机图像分析技术，多孔介质材料内部的微观三维孔隙结构可以被数值模型表征。通过对多孔材料的三维模型（根据扫描电镜结构重构或根据材料结构重新建立）进行像素化处理，即将模型划分为若干全等立方体像素，像素的中心点为固体则视该像素为固体，像素的中心点为孔隙则视该像素为孔隙，可以分析多孔材料的孔隙结构信息，当然也包括连通孔隙率。但是众所周知，多孔材料的连通性在多数情况下主要受小孔影响，而基于像素化处理后的多孔材料三维模型所能识别的最小孔隙却受分辨率的限制影响。越高的分辨率下多孔材料三维模型孔结构展示的也越精细，计算结果也越精确，但与此同时运行处理的内容也随之增加，效率下降；低分辨率下虽然能快速运算获取结果，但粗糙的多孔材料三维模型带来的却是准确度上的大幅降低。对不规则多孔材料采取像素化方法乃无奈之举，而对于具有特殊结构的球体堆积类多孔材料而言，目前未有不受分辨率限制的孔隙连通性分析方法，导致对不规则多孔材料中非连通孔隙的孔隙连通性难以准确评估。
发明内容
本发明的目的是提供一种水泥基材料颗粒堆积体系孔隙连通性分析方法。
本发明为解决上述技术问题采取的技术方案是：
一种水泥基材料颗粒堆积体系孔隙连通性分析方法，所述球体堆积体系孔隙连通性分析方法为首先建立球体颗粒堆积模型，其次对球心点集执行Delaunay四面体剖分以把堆积模型划分为由多个四面体组成的集合后，逐一获取集合中四面体的孔隙体积以及判断每个四面体孔隙的连通性，最后进行球体颗粒堆积模型中连通孔隙率的统计计算过程。
作为优选方案：建立球体颗粒堆积模型的过程为：
对于指定球体颗粒粒径分布和孔隙率的体系，首先确定立方体填充空间的尺寸大小，然后把球体颗粒按尺寸从大到小依次投入立方体填充空间中，每个球体颗粒的坐标由蒙特卡洛法随机生成，直至所有颗粒全部投放完毕，形成球体颗粒堆积模型。
作为优选方案：基于球体颗粒堆积模型，对球心点集执行Delaunay四面体剖分的操作过程为：
根据球体颗粒堆积模型信息，把所有球体中心的坐标进行整理形成二维矩阵，二维矩阵有三列，分别为球心的横坐标、纵坐标和列坐标，矩阵中每一行代表一个球体的球体坐标，采用MATLAB中的Delaunay Triangulation函数作用该二维矩阵，即实现对球心点集的Delaunay四面体剖分；在二维情况下，矩阵只有两列，分别为圆的横坐标集和列坐标集，对矩阵作用Delaunay Triangulation函数以执行对圆心点集的Delaunay三角剖分。
作为优选方案：对球心点集执行Delaunay四面体剖分后，堆积模型被划分为由多个四面体组成的集合，逐个确定集合中四面体孔隙体积的过程为：
经Delaunay四面体剖分后的颗粒堆积体系由相互独立、互不重叠的多个四面体组成，其中每个四面体由四个顶角球体颗粒占据的固相以及孔隙相组成，固相为三个角点颗粒在表面上截交的扇形区域总和，孔隙相即为四面体减去四个顶角球体颗粒所占据固相的余相，四面体的孔隙相即为四面体孔隙，孔隙相对应的体积即为该四面体的孔隙体积；
对允许重叠存在的球体堆积体系，因其结构复杂，采取下述方案进行计算：
在四面体中随机放入N个粒子，粒子坐标由蒙特卡洛法生成，统计落在孔隙中的粒子个数n，则四面体中的孔隙体积Vp可由公式一计算：
Vp = (n / N) * Vt （1）。
作为优选方案：判断每个四面体孔隙的连通性的过程为首先遍历判断相邻四面体孔隙连通关系后，再判断所有四面体孔隙的连通性；
四面体具有四个三角形表面，每个三角形表面由固相和孔隙相组成，其中表面固相为三个角点颗粒在表面上截交的扇形区域总和，表面孔隙相即为该四面体表面减去相应固相的余相；表面孔隙相是连接该四面体孔隙与其临近四面体孔隙的通道最窄处的横截面，即为咽喉横截面；
当任意两个相邻四面体孔隙之间的咽喉横截面被球体颗粒完全填充时，表明该两个相邻四面体孔隙之间的连接关系为相互阻绝关系，即不存在连通关系；当任意两个相邻四面体孔隙之间的咽喉横截面未被球体颗粒完全填充，存在空隙，表明该两个相邻四面体孔隙之间的连接关系为连通关系，依次类推，遍历判断相邻四面体孔隙连通关系即可；
判断所有四面体孔隙的连通性的过程为：
采用连通性算法判断四面体孔隙为连通孔隙还是非连通孔隙，具体步骤为：
步骤一：为每个四面体编号并将按序储存在一个数组中；
步骤二：遍历寻找所有相邻四面体孔隙之间的连通关系；
步骤三：遍历输入相邻四面体孔隙之间的连通关系后，当两个相邻四面体孔隙连通，则把该两个相邻四面体孔隙以及其它与该两个相邻四面体同编号的四面体孔隙的编号修改为同一个值；
步骤四：根据四面体孔隙的编号，把多个四面体孔隙划分为不同的孔隙集，处于同一个孔隙集中的所有四面体孔隙的编号相同；
步骤五：对非周期性重复球体堆积体系，逐一检查每个孔隙集中是否存在处于对立面的两个孔隙，当一个孔隙集中存在处于对立面的两个孔隙时，表明该孔隙集中的孔隙为连通孔隙；当一个孔隙集中不存在处于对立面的两个孔隙时，表明该孔隙集中的孔隙为非连通孔隙；
步骤六：对服从周期性重复边界的球体堆积体积，逐一检查多个孔隙集中是否存在满足周期性重复条件的两个孔隙，当一个孔隙集中存在满足周期性重复条件的两个孔隙时，表明该孔隙集中的所有孔隙为连通孔隙；当一个孔隙集中不存在满足周期性重复条件的两个孔隙时，表明该孔隙集中的所有孔隙为非连通孔隙，依次类推从而完成所有Delaunay四面体孔隙的连通性判断过程。
作为优选方案：球体颗粒堆积模型的连通孔隙率计算过程如下：
对所有四面体孔隙的连通性判断完毕后，将判断为连通孔隙的四面体孔隙的体积进行累加求和，连通孔隙的四面体孔隙的体积总和为Vcp；球体颗粒堆积模型的连通孔隙率即为Vcp除以球体颗粒堆积模型的总体积；球体颗粒堆积模型的孔隙连通度即为Vcp除以球体颗粒堆积模型的总孔隙体积，从而得到球体颗粒堆积模型的连通孔隙率以及孔隙连通度。
本发明相对于现有技术具有以下有益效果：
一、本发明为一种基于Delaunay三角剖分的球体堆积体系孔隙连通性分析方法。本发明的计算原理科学合理，步骤简单，对操作人员的操作经验无要求，获取结果时效快且稳定准确。
二、本发明立足于球体堆积的本征几何关系，避免了像素化处理方法中依靠高分辨率来获取高精度的做法，能够快速划分四面体孔隙区域，识别连通孔隙和非连通孔隙，算法执行速度快，计算效率高。
三、本发明能够处理多级配不等径球体堆积体系，填补现有对多级配不等径球体堆积体系的内部连通孔隙和非连通孔隙获取难以准确的空白，通过连通孔隙和非连通孔隙计算得到球体颗粒堆积模型的连通孔隙率以及孔隙连通度，为后续多级配不等径球体堆积体系的进一步研究提供相关数据的有效获取方法。
四、本发明适用于计算低孔隙率堆积或高孔隙率堆积的多种颗粒堆积情况，计算原理合理且周全，计算结果准确，更加符合实际需求。
五、本发明能够进一步拓展至不规则多面体颗粒堆积体系的孔径分布计算，具有广泛的应用场景和未来发展优化的潜力空间。
六、本发明提出的孔隙连通性分析方法的计算结果能够作为判断水泥基材料抵抗侵蚀介质侵入能力的依据，为后续预测水泥基材料耐久性奠定基础。
附图说明
下面结合附图对本发明作进一步的说明
图1是本发明的流程框图；
图2是三维Delaunay三角剖分后四个球形颗粒组成的单个四面体的立体结构示意图；
图3a为相邻四面体孔隙为相互连通关系的结构示意图；
图3b为相邻四面体孔隙为非相互连通关系的结构示意图；
图4是四面体孔隙连通性判断算法示意图；
图5是单粒径重叠球体堆积连通孔隙率和总孔隙率关系图；
图6是单粒径重叠球体堆积孔隙连通度和总孔隙率关系图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一个实施特例，而非全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
在此，还需要说明的是，为避免因不必要的细节而模糊本发明，在附图中仅仅展示了与本发明的方案密切相关的处理步骤，而省略了与本发明关系不大的细枝末节。
具体实施方式一：结合图1、图2、图3、图3a、图3b、图4、图5和图6说明本实施方式，本实施方式中所述球体堆积体系孔隙连通性分析方法为首先建立球体颗粒堆积模型，其次对球心点集执行Delaunay四面体剖分以把堆积模型划分为由多个四面体组成的集合后，逐一获取集合中四面体的孔隙体积以及判断每个四面体孔隙的连通性，最后进行球体颗粒堆积模型中连通孔隙率的统计计算过程。
本实施方式中球体颗粒堆积模型为可重叠随机球体堆积三维模型。
本实施方式中多个四面体组成的集合为Delaunay四面体剖分球心点集的结果。
具体实施方式二：本实施方式为具体实施方式一的进一步限定，所述球体堆积体系孔隙连通性分析方法中建立球体颗粒堆积模型的过程为：
对于指定球体颗粒粒径分布和孔隙率的体系，首先确定立方体填充空间的尺寸大小，然后根据球体颗粒粒径尺寸从大到小依次投入立方体填充空间中，每个球体颗粒的坐标由蒙特卡洛法随机生成，直至所有颗粒全部投放完毕，形成球体颗粒堆积模型。
具体实施方式三：本实施方式为具体实施方式一的进一步限定，本实施方式中把所有球体中心的坐标进行整理形成二维矩阵的过程为现有计算方法，Delaunay四面体剖分的处理原理与现有的Delaunay四面体剖分的处理原理相同。
步骤S1根据已知的球体颗粒粒径分布信息，利用蒙特卡洛算法实现球体颗粒的随机分布建模的具体过程为：将需要生成的N个球体按半径从大到小排序，当i时，随机生成第i个球体，而后进行下一步操作，即令i=i+1，当i小于或等于N时，则再重新随机生成第i个球体，重复检测操作，当i大于N时，所有球体生成完毕。
具体实施方式四：本实施方式为具体实施方式一、二或三的进一步限定，基于球体颗粒堆积模型，对球心点集执行Delaunay四面体剖分的操作过程为：
根据球体颗粒堆积模型信息，把所有球体中心的坐标进行整理形成二维矩阵，二维矩阵有三列，分别为球心的横坐标、纵坐标和列坐标，矩阵中每一行代表一个球体的球体坐标，采用MATLAB中的Delaunay Triangulation函数作用该二维矩阵，即对球心点集执行Delaunay四面体剖分；在二维情况下，矩阵只有两列，分别为圆的横坐标集和列坐标集，对矩阵作用Delaunay Triangulation函数以执行对圆心点集的Delaunay三角剖分。
本实施方式中球心点集为所有球形颗粒的球心坐标组合成的集合。
本实施方式中二维Delaunay三角剖分的结果示意如图2所示，图中圆代表颗粒，直线段集是对圆心点集的划分结果。
具体实施方式五：本实施方式为具体实施方式一、二、三或四的进一步限定，对球心点集执行Delaunay四面体剖分后，确定每个四面体对应的孔隙孔径和孔隙体积的步骤为：
经Delaunay四面体剖分后的颗粒堆积体系由相互独立、互不重叠的多个四面体组成，其中每个四面体由两部分构成：顶角球体颗粒占据的固相部分以及除固相部分以外的孔隙相。四面体中孔隙相可称为四面体孔隙，其对应的体积即为该四面体的孔隙体积。
具体实施方式六：本实施方式为具体实施方式一、二、三、四或五的进一步限定，对球心点集执行Delaunay四面体剖分后，堆积模型被划分为由多个四面体组成的集合，逐个确定Delaunay四面体剖分球心点集中四面体孔隙体积的过程为：
经Delaunay四面体剖分后的颗粒堆积体系由相互独立、互不重叠的多个四面体组成，其中每个四面体由顶角球体颗粒占据的固相容积以及孔隙容积，孔隙容积为四面体总容积减去顶角球体颗粒占据的固相容积的差值，孔隙容积即为四面体孔隙，孔隙容积对应的体积即为该四面体孔隙的孔隙体积；
对允许重叠存在的球体堆积体系，因其结构复杂，采取下述方案进行计算：
在四面体中随机放入N个粒子，粒子坐标由蒙特卡洛法生成，统计落在孔隙中的粒子个数n，则四面体中的孔隙体积Vp可由公式一计算：
Vp = (n / N) * Vt （1）
由此完成了遍历计算Delaunay四面体孔隙体积。
具体实施方式七：本实施方式为具体实施方式一、二、三、四、五或六的进一步限定，遍判断每个四面体孔隙的连通性的过程为首先遍历判断相邻四面体孔隙连通关系后，再判断所有四面体孔隙的连通性；
遍历判断相邻四面体孔隙连通关系的过程为：四面体具有四个三角形表面，每个三角形表面由固相和孔隙相组成，其中表面固相为三个角点颗粒在表面上截交的扇形区域总和，表面孔隙相即为该四面体表面减去相应固相的余相；表面孔隙相是连接该四面体孔隙与其临近四面体孔隙的通道最窄处的横截面，即为咽喉横截面；
当任意两个相邻四面体孔隙之间的咽喉横截面被球体颗粒完全填充时，表明该两个相邻四面体孔隙之间的连接关系为相互阻绝关系，即不存在连通关系；当任意两个相邻四面体孔隙之间的咽喉横截面未被球体颗粒完全填充，存在空隙，表明该两个相邻四面体孔隙之间的连接关系为连通关系，依次类推，遍历判断相邻四面体孔隙连通关系即可；
判断所有四面体孔隙的连通性的过程为：
采用连通性算法判断四面体孔隙为连通孔隙还是非连通孔隙，具体步骤为：
步骤一：为每个四面体编号并将按序储存在一个数组中；
步骤二：遍历寻找所有相邻四面体孔隙之间的连通关系；
步骤三：遍历输入相邻四面体孔隙之间的连通关系后，当两个相邻四面体孔隙连通，则该两个相邻四面体孔隙以及其它与该两个相邻四面体同编号的四面体孔隙的编号修改为同一个值；
步骤四：根据四面体孔隙的编号，把多个四面体孔隙划分为不同的孔隙集，处于同一个孔隙集中的所有四面体孔隙的编号相同；
步骤五：对非周期性重复球体堆积体系，逐一检查多个孔隙集中是否存在处于对立面的两个孔隙，当一个孔隙集中存在处于对立面的两个孔隙时，表明该孔隙集中的孔隙为连通孔隙；当一个孔隙集中不存在处于对立面的两个孔隙时，表明该孔隙集中的孔隙为非连通孔隙；
步骤六：对服从周期性重复边界的球体堆积体积，逐一检查每个孔隙集中是否存在满足周期性重复条件的两个孔隙，当一个孔隙集中存在满足周期性重复条件的两个孔隙时，表明该孔隙集中的所有孔隙为连通孔隙；当一个孔隙集中不存在满足周期性重复条件的两个孔隙时，表明该孔隙集中的所有孔隙为非连通孔隙，依次类推从而完成所有Delaunay四面体孔隙的连通性判断过程。
具体实施方式八：本实施方式为具体实施方式一、二、三、四、五、六或七的进一步限定，球体颗粒堆积模型的连通孔隙率计算过程如下：
对所有四面体孔隙的连通性判断完毕后，将判断为连通孔隙的四面体孔隙的体积进行累加求和，连通孔隙的四面体孔隙的体积总和为Vcp；球体颗粒堆积模型的连通孔隙率即为Vcp除以球体颗粒堆积模型的总体积；球体颗粒堆积模型的孔隙连通度即为Vcp除以球体颗粒堆积模型的总孔隙体积，从而得到球体颗粒堆积模型的连通孔隙率以及孔隙连通度。
本实施方式中球体颗粒堆积模型的连通孔隙率即为堆积体系连通孔隙率，球体颗粒堆积模型的孔隙连通度即为堆积体系孔隙连通度。
本发明为一种水泥基材料颗粒堆积体系孔隙连通性分析方法，主要应用于球形颗粒，其他不规则形状的颗粒也可简化为球形颗粒后利用本发明的颗粒堆积体系孔隙连通性分析方法进行处理。
具体实施方式九：本实施方式为具体实施方式一、二、三、四、五、六、七或八的进一步限定，本实施方式是基于建立的球体堆积模型，通过对球心点集执行Delaunay四面体剖分，然后采取提出的连通性判断算法确定Delaunay四面体孔隙的连通性，从而获取孔隙结构的连通孔隙率，具体实施步骤如下：
S1、可重叠随机球体堆积三维模型的建立：首先设定颗粒的立方体填充空间的尺寸大小，然后在立方体空间中依次投入球体颗粒，其坐标由蒙特卡洛法随机生成，直至所有颗粒全部投放完毕或体系到达目标孔隙率。
S2、球心点集的Delaunay四面体剖分：把所有堆积球体的球心坐标整理成二维矩阵，矩阵有三列，分别代表球心的横坐标、纵坐标和列坐标，矩阵中每一行代表一个球体的球心坐标，采用MATLAB中的Delaunay Triangulation函数作用该二维坐标矩阵，实现对球心点集的Delaunay四面体剖分。
经Delaunay四面体剖分后的颗粒填充空间由相互不重叠的四面体组成，每个四面体由四个顶角球体颗粒和孔隙部分构成；图2是Delaunay四面体的组成示例，图中灰色渐变填充的是球体颗粒，四面体中除球体颗粒以外的部分是空隙，该部分对应的体积是孔隙体积。
S3、四面体孔隙的体积计算：对允许重叠存在的球体堆积体系中的Delaunay四面体，采用几何关系计算四面体体积Vt后在四面体中随机放入N个粒子，粒子坐标由蒙特卡洛法随机生成，统计落在孔隙中的粒子个数n，则四面体孔隙体积Vp可由下式计算：
Vp = (n / N) * Vt
S4、相邻四面体孔隙的连通关系：Delaunay四面体四个表面的空隙部分对应于孔隙四个咽喉的横截面，是连接该四面体孔隙与周围四个四面体孔隙的通道；若某两个相邻四面体孔隙之间的咽喉横截面通道被球体颗粒完全填充，则它们之间被相互阻绝，不存在连通关系，如图3a示意；反之，若两个相邻四面体孔隙之间的咽喉横截面通道未被球体颗粒完全填充，留有部分空隙，则它们之间存在连通关系，如图3b示意。
S5、所有四面体孔隙的连通性判断：采用连通性算法判断Delaunay四面体孔隙为连通孔隙还是非连通孔隙，算法示意见图4，具体步骤解析如下：
①为每个四面体编号并将按序储存在一个数组中；
②遍历寻找所有相邻四面体孔隙之间的连通关系；
③遍历输入相邻四面体孔隙之间的连通关系，若两个相邻四面体孔隙连通，则把这两个相邻四面体孔隙以及其它与它们同编号的四面体孔隙的编号修改为同一个值；
④根据四面体孔隙的编号，把它们划分为不同的孔隙集，同一个孔隙集中的所有四面体孔隙的编号相同；
⑤对非周期性重复球体堆积体系，检查孔隙集中是否存在处于对立面的两个孔隙，若存在，则该孔隙集中的孔隙为连通孔隙，反之则为非连通孔隙；
⑥对服从周期性重复边界的球体堆积体积，检查孔隙集中是否存在满足周期性重复条件的两个孔隙，若存在，则该孔隙集中的所有孔隙为连通孔隙，反之则为非连通孔隙；
⑦至此，所有Delaunay四面体孔隙的连通性判断完毕。
S6、堆积体系的连通孔隙率计算：通过对S4中判断为连通孔隙的四面体孔隙的体积进行累加求和，记为Vcp；堆积体系的连通孔隙率即为Vcp除以堆积体系的总体积；堆积体系的孔隙连通度即为Vcp除以堆积体系的总孔隙体积。
结合本发明的有益效果以及说明书附图1至6说明以下实施例：
实施例一：
选择允许重叠单粒径球体堆积体系，采用蒙特卡洛法建立不同孔隙率的球体堆积三维模型，随后对每个模型都进行Delaunay四面体划分、四面体的孔隙体积遍历计算、四面体孔隙的连通性判断以及连通孔隙率和孔隙连通度计算，最后进行整理，分析孔隙率变化对体系中孔隙连通性的影响，以对算法的运行和计算结果进行展示：
堆积体系的连通孔隙率和总孔隙率关系如图5所示，图5展示了不同孔隙率堆积体系的连通孔隙率计算结果，孔隙率-连通孔隙率线条趋于斜率为1的直线，表明对于大多数孔隙率的堆积体系而言，内部孔隙均为连通孔隙；但在极低孔隙率的情况下，即孔隙率低于0.1时，直线略有弯折，呈现外凸趋势，显示此时的堆积体系中已经存在部分非连通孔隙。
堆积体系的孔隙连通度和总孔隙率关系如图6所示，图6展示了不同孔隙率堆积体系的孔隙连通度计算计算，对孔隙率大于0.1的单粒径可重叠球体颗粒堆积体系，其孔隙连通度均为1，表明其内部所有孔隙均为连通孔隙，具被可预见的良好传输性能；对孔隙率小于0.1的单粒径可重叠球体颗粒堆积体系，随着孔隙率的降低，体系中的孔隙连通度迅速下降，内部原先的连通孔隙逐渐转换为非连通孔隙，到孔隙率为0.015时，体系的孔隙连通度已然为0，此时，所有孔隙均变为非连通孔隙，体系失去物质传输能力。
在本说明书的描述中，参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。
Claims (2)
1.一种水泥基材料颗粒堆积体系孔隙连通性分析方法，球体堆积体系孔隙连通性分析方法为首先建立球体颗粒堆积模型，其次对球心点集执行Delaunay四面体剖分以把堆积模型划分为由多个四面体组成的集合后，逐一获取集合中四面体的孔隙体积以及判断每个四面体孔隙的连通性，最后进行球体颗粒堆积模型中连通孔隙率的统计计算过程；
建立球体颗粒堆积模型的过程为：
对于指定球体颗粒粒径分布和孔隙率的体系，首先确定立方体填充空间的尺寸大小，然后把球体颗粒按尺寸从大到小依次投入立方体填充空间中，每个球体颗粒的坐标由蒙特卡洛法随机生成，直至所有颗粒全部投放完毕，形成球体颗粒堆积模型；
对球心点集执行Delaunay四面体剖分的操作过程为：
根据球体颗粒堆积模型信息，把所有球体中心的坐标进行整理形成二维矩阵，二维矩阵有三列，分别为球心的横坐标、纵坐标和列坐标，矩阵中每一行代表一个球体的球体坐标，采用MATLAB中的Delaunay Triangulation函数作用该二维矩阵，即实现对球心点集的Delaunay四面体剖分；在二维情况下，矩阵只有两列，分别为圆的横坐标集和列坐标集，对矩阵作用Delaunay Triangulation函数以执行对圆心点集的Delaunay三角剖分；
对球心点集执行Delaunay四面体剖分后，堆积模型被划分为由多个四面体组成的集合，逐个确定集合中四面体孔隙体积的过程为：
经Delaunay四面体剖分后的颗粒堆积体系由相互独立、互不重叠的多个四面体组成，其中每个四面体由四个顶角球体颗粒占据的固相以及孔隙相组成，固相为三个角点颗粒在表面上截交的扇形区域总和，孔隙相即为四面体减去四个顶角球体颗粒所占据固相的余相，四面体的孔隙相即为四面体孔隙，孔隙相对应的体积即为该四面体的孔隙体积；
对允许重叠存在的球体堆积体系，因其结构复杂，采取下述方案进行计算：
在四面体中随机放入N个粒子，粒子坐标由蒙特卡洛法生成，统计落在孔隙中的粒子个数n，Vt为四面体体积，则四面体中的孔隙体积Vp可由公式一计算：
Vp = (n / N) * Vt（1）；
其特征在于：判断每个四面体孔隙的连通性的过程为首先遍历判断相邻四面体孔隙连通关系后，再判断所有四面体孔隙的连通性；
四面体具有四个三角形表面，每个三角形表面由固相和孔隙相组成，其中表面固相为三个角点颗粒在表面上截交的扇形区域总和，表面孔隙相即为该四面体表面减去相应固相的余相；表面孔隙相是连接该四面体孔隙与其临近四面体孔隙的通道最窄处的横截面，即为咽喉横截面；
当任意两个相邻四面体孔隙之间的咽喉横截面被球体颗粒完全填充时，表明该两个相邻四面体孔隙之间的连接关系为相互阻绝关系，即不存在连通关系；当任意两个相邻四面体孔隙之间的咽喉横截面未被球体颗粒完全填充，存在空隙，表明该两个相邻四面体孔隙之间的连接关系为连通关系，依次类推，遍历判断相邻四面体孔隙连通关系即可；
判断所有四面体孔隙的连通性的过程为：
采用连通性算法判断四面体孔隙为连通孔隙还是非连通孔隙，具体步骤为：
步骤一：为每个四面体编号并将按序储存在一个数组中；
步骤二：遍历寻找所有相邻四面体孔隙之间的连通关系；
步骤三：遍历输入相邻四面体孔隙之间的连通关系后，当两个相邻四面体孔隙连通，则把该两个相邻四面体孔隙以及其它与该两个相邻四面体同编号的四面体孔隙的编号修改为同一个值；
步骤四：根据四面体孔隙的编号，把多个四面体孔隙划分为不同的孔隙集，处于同一个孔隙集中的所有四面体孔隙的编号相同；
步骤五：对非周期性重复球体堆积体系，逐一检查每个孔隙集中是否存在处于对立面的两个孔隙，当一个孔隙集中存在处于对立面的两个孔隙时，表明该孔隙集中的孔隙为连通孔隙；当一个孔隙集中不存在处于对立面的两个孔隙时，表明该孔隙集中的孔隙为非连通孔隙；
步骤六：对服从周期性重复边界的球体堆积体积，逐一检查多个孔隙集中是否存在满足周期性重复条件的两个孔隙，当一个孔隙集中存在满足周期性重复条件的两个孔隙时，表明该孔隙集中的所有孔隙为连通孔隙；当一个孔隙集中不存在满足周期性重复条件的两个孔隙时，表明该孔隙集中的所有孔隙为非连通孔隙，依次类推从而完成所有Delaunay四面体孔隙的连通性判断过程。
2.根据权利要求1所述的一种水泥基材料颗粒堆积体系孔隙连通性分析方法，其特征在于：球体颗粒堆积模型的连通孔隙率计算过程如下：
对所有四面体孔隙的连通性判断完毕后，将判断为连通孔隙的四面体孔隙的体积进行累加求和，连通孔隙的四面体孔隙的体积总和为Vcp；球体颗粒堆积模型的连通孔隙率即为Vcp除以球体颗粒堆积模型的总体积；球体颗粒堆积模型的孔隙连通度即为Vcp除以球体颗粒堆积模型的总孔隙体积，从而得到球体颗粒堆积模型的连通孔隙率以及孔隙连通度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110637103.9A CN113405965B (zh) | 2021-06-08 | 2021-06-08 | 一种水泥基材料颗粒堆积体系孔隙连通性分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110637103.9A CN113405965B (zh) | 2021-06-08 | 2021-06-08 | 一种水泥基材料颗粒堆积体系孔隙连通性分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113405965A CN113405965A (zh) | 2021-09-17 |
| CN113405965B true CN113405965B (zh) | 2022-06-28 |
Family
ID=77676938
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110637103.9A Active CN113405965B (zh) | 2021-06-08 | 2021-06-08 | 一种水泥基材料颗粒堆积体系孔隙连通性分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113405965B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116337708B (zh) * | 2022-12-06 | 2023-11-17 | 武汉理工大学 | 提取多孔材料孔隙特征和流体分布的方法与系统 |
| CN116343958B (zh) * | 2023-03-03 | 2024-05-10 | 湖北工业大学 | 多相体系连通性的判断方法及系统 |
| CN116740087B (zh) * | 2023-05-31 | 2024-05-24 | 湖北工业大学 | 用于多相体系连通性判断的区域分割高效搜索方法及系统 |
| CN116884540B (zh) * | 2023-06-15 | 2024-05-10 | 湖北工业大学 | 多相材料三维模型生成与连通性判断方法 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1733329A4 (en) * | 2004-03-31 | 2015-07-29 | Exxonmobil Upstream Res Co | METHOD FOR SIMULATING AND ESTIMATING THE PROPERTIES OF GRES |
| CN104237103B (zh) * | 2014-09-23 | 2017-10-17 | 中国石油天然气股份有限公司 | 一种孔隙连通性定量表征方法及装置 |
| CN107063961A (zh) * | 2016-12-15 | 2017-08-18 | 中国矿业大学 | 一种水泥水化三维孔隙结构特征的表征方法 |
| CN108387495B (zh) * | 2018-01-22 | 2020-03-31 | 青岛理工大学 | 一种多孔混凝土孔隙率计算和孔隙参数表征方法 |
| CN109242985B (zh) * | 2018-10-29 | 2020-06-05 | 中国科学院力学研究所 | 一种从三维图像确定孔隙结构关键参数的方法 |
| CN109543350B (zh) * | 2018-12-20 | 2023-03-24 | 哈尔滨工业大学（深圳） | 基于细观模型重构的透水混凝土性能预测方法及系统 |
| CN109977469B (zh) * | 2019-02-22 | 2023-04-14 | 中国石油大学(华东) | 一种基于Voronoi图的二维多孔介质模型构建方法 |
| CN110599477B (zh) * | 2019-09-16 | 2020-11-03 | 南昌工程学院 | 一种植草混凝土有效孔隙率的检测方法及系统 |
| CN112765868B (zh) * | 2020-12-29 | 2022-11-25 | 哈尔滨工业大学 | 一种早期水泥石毛细水中可冻水含量的预判方法 |
| CN112903555B (zh) * | 2021-01-14 | 2023-05-26 | 中国石油天然气股份有限公司 | 考虑孔隙各向异性的多孔介质渗透率计算方法及装置 |
-
2021
- 2021-06-08 CN CN202110637103.9A patent/CN113405965B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN113405965A (zh) | 2021-09-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113405965B (zh) | 一种水泥基材料颗粒堆积体系孔隙连通性分析方法 | |
| Lu et al. | Critical assessment of two approaches for evaluating contacts between super-quadric shaped particles in DEM simulations | |
| Øren et al. | Process based reconstruction of sandstones and prediction of transport properties | |
| Van Marcke et al. | An improved pore network model for the computation of the saturated permeability of porous rock | |
| CN111222235B (zh) | 一种随机生成多结构层沥青路面离散元模型的方法 | |
| CN104834783A (zh) | 点蚀随机分布的圆柱壳数值模型的参数化构建方法 | |
| CN110806371B (zh) | 单重煤岩孔隙结构模型的构建方法 | |
| CN109979009A (zh) | 一种具有功能梯度的胶结颗粒材料三维细观数字模型重构方法 | |
| CN113281149B (zh) | 一种节理岩体的表征单元体积尺度综合取值方法 | |
| CN113971718B (zh) | 一种对三维点云模型进行布尔运算的方法 | |
| CA3062555A1 (en) | Porous medium permeability prediction method based on machine image intelligent learning | |
| CN116259375A (zh) | 沥青混合料三维细观模型的构建方法、系统及存储介质 | |
| Walker et al. | Analysis of grain-scale measurements of sand using kinematical complex networks | |
| CN111474574A (zh) | 基于压缩感知的地震采集观测系统的生成方法及装置 | |
| CN111751403A (zh) | 一种热障涂层数值重构模型测试方法和装置 | |
| CN113392570B (zh) | 一种水泥基材料颗粒堆积体系孔隙结构均质度的评估方法 | |
| Stroeven et al. | Porosimetry by random node structuring in virtual concrete | |
| CN116429295A (zh) | 一种评估岩体结构面接触应力分布的方法 | |
| Bloom et al. | Measurement of porosity in granular particle distributions using adaptive thresholding | |
| CN113405966B (zh) | 一种水泥基材料颗粒堆积体系孔径分布计算方法 | |
| Agrawal et al. | Voronoi tessellation-based algorithm for determining rigorously defined classical and generalized geometric pore size distributions | |
| CN115221747A (zh) | 任意形状壳体结构中随机点蚀的建模方法 | |
| CN111104728B (zh) | 结构表面磨损的仿真预测方法和装置 | |
| Lim et al. | Statistical models to describe the structure of porous ceramic membranes | |
| Byholm et al. | The application of morphological algorithms on 3-dimensional porous structures for identifying pores and gathering statistical data. |
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 |