一种计算凸形多面体粒子堆积体系孔径分布的方法
【专利摘要】本发明公开了一种计算凸形多面体粒子堆积体系孔径分布的方法，对水泥浆体的内部结构进行合理的假设，形成三维空间凸形多面体粒子堆积体系，凸形多面体粒子堆积间隙即为内部孔隙结构。将体系进行单元划分，判断每个单元格的n个点与凸形多面体粒子的包含关系，从而得出该单元格的分类，单元格的种类有固相、孔相和混合相。对其中的混合相单元格进行更小的单元划分，最终得到全部孔径的大小。通过本发明计算得到模拟孔径分布准确，并且计算速度快。
【专利说明】一种计算凸形多面体粒子堆积体系孔径分布的方法
【技术领域】
[0001] 本发明涉及到一种计算凸形多面体粒子堆积体系孔径分布的方法。
【背景技术】
[0002] 水泥浆体的孔隙结构是关系到土木工程材料的渗透性，强度和耐久性的关键因 素。当材料受到外部荷载，冻融循环和各种盐类侵蚀时，材料的内部结构会发生严重变化。 为了提高材料的各项性能，使材料在受到损伤时具有足够的抵抗能力，需要对材料的内部 结构有一定深入的了解，孔径分布是材料结构的重要参数之一。
[0003] 通过压汞法，氮气吸脱附法，微观扫描电镜和热孔计法等试验可以得到水泥浆体 内部的孔径分布情况，但是这些试验结果往往受制于试验条件和试验手段，且试验过程中 会对水泥浆体内部的孔隙结构有破坏作用。为把人力和资源从大量试验中解脱出来，从本 质上把握材料的结构性能变化规律，用于模拟材料微细观结构的计算机模型研究成为当前 的一个研究热点。水泥水化后，内部结构非常复杂，包含大量未水化的水泥颗粒，内部和外 部的水化产物。这些固体之间相互交错连接，中间形成大量的孔隙。为了研究这些孔隙的大 小，需要将固体或孔隙进行简化，即将水泥浆体简化为规则的三维立体体系，而这些固体简 化为凸多面体粒子，杂乱无序的排列在三维体系中，而凸多面体粒子间隙即为所要研究的 孔相，其大小即为孔径。基于凸多面体粒子的蒙特卡罗法堆积及相应的水化反应，得到混凝 土的模拟微观及细观结构【van Breugel K. Simulation of hydration and formation of structure in hardening cement-based materials. Ph. D thesis,Delft University of Technology，Delft，The Netherlands，199L 】【Xu W.X·，Chen H.S.，Lv Z. An overlapping detection algorithm for random sequential packing of elliptical particles. Physica A，390(2011)2452-67.】。对于孔径分布的模型分析，已提出的有丢球算法，通过 向模拟结构中依次丢入由大到小的球体，判断球体是否能够放入来计算模拟结构的孔径 分布【Ye G. Experimental study and numerical simulation of the development of the microstructure and permeability of cementitious materials, PhD thesis, Delft University of Technology, Delft, 2003.】。丢球法仅适用于球形粒子堆积模拟的微观结 构，且计算量大，计算时间较长，尤其是丢入较小粒径的球体时。本发明的计算方法不仅可 以计算凸多面体粒子堆积的孔径分布，且适用于凸多面体粒子堆积的模拟微观结构，而且 相对已提出的算法而言，计算时间大大缩短，提高了计算效率。
【发明内容】
[0004] 发明目的：本发明所要解决的技术问题在于克服了原有孔径分布计算算法计算时 间长，效率低，适用性窄等问题，而且大大提高了计算孔径的效率，能够合理的计算出水泥 浆体的孔径分布情况。
[0005] 技术方案：本发明提出一种计算凸形多面体粒子堆积体系孔径分布的方法，通过 将混合相单元格不断地进行划分最后计算求得孔径，包括以下步骤：
[0006] 步骤1、将水泥浆体假定为三维凸形多面体粒子堆积体系，其微观结构是由固相和 毛细孔相组成的两相结构体系；
[0007] 步骤2、在所述三维凸形多面体粒子堆积体系中确定每一个凸形多面体的顶点坐 标和每一个面的空间方程；
[0008] 步骤3、在整体坐标系中，将三维凸形多面体粒子堆积体系按照边长屯进行单元划 分，并且确定出每一个单元格的η个合适点的坐标吨=L/N，N取大于等于2的整数，L为 所述三维凸形多面体粒子堆积体系的边长，且屯< 2 μ m，2 μ m为水泥净浆的最大孔径；
[0009] 步骤4、利用单元格η个点的坐标、凸形多面体的顶点坐标和每个平面的空间方 程，判断该单元格与凸形多面体粒子的关系，从而得出单元格的种类；
[0010] 步骤5、将其中的混合相单元格按照边长d2进行划分，d2 =屯/N'，Ν'取大于等于 2的整数，不断将这个三维连续结构体系进行划分，最终可得出全部孔隙的直径。
[0011] 作为优选，所述微观结构的固相为凸形多面体粒子，所述凸形多面体粒子包括正 四面体粒子、正六面体粒子、正八面体粒子、正十二面体粒子和正二十面体粒子。
[0012] 作为优选，将所包含的固相和毛细孔相的三维体系全部划分。
[0013] 作为优选，所述η个合适点的坐标包括8个顶点坐标和1个中心点坐标，利用每个 单元格9个点的坐标、凸形多面体粒子的顶点坐标和每个平面的空间方程，判断单元格与 凸形多面体粒子的关系。
[0014] 作为优选，所述单元格的种类包含：固相、孔相和混合相。
[0015] 作为优选，仅将其中的混合相单元格不断进行划分，最终得到全部孔径的大小。
[0016] 作为优选，取 Ν' = 2,即 d2 = d/2。
[0017] 使用时，具体操作如下：
[0018] 步骤1、将水泥浆体假定为三维凸形多面体粒子堆积体系，其微观结构是由固相 (包括未水化水泥颗粒和水化产物）和毛细孔相组成的两相结构体系，是由正四面体粒子、 正六面体粒子、正八面体粒子、正十二面体粒子、正二十面体粒子组成的堆积体系；
[0019] 步骤2、在三维凸形多面体粒子堆积体系中选择全局坐标系为参照坐标系，并且根 据所设定的坐标系，计算出每个凸形多面体粒子的顶点坐标和每个平面的空间方程，如Si 的平面方程如下：
[0020] Α (χ-χ^+B (y-y^+C (z-z^ =0
[0021] 其中，yp Zi为Si上的Pi顶点坐标。
[0022] 步骤3、将三维体系按照边长屯进行单元划分，并且确定出每一个单元格八个顶点 和一个中心的坐标。屯=L/N，N取大于等于2的整数，L为所述三维凸形多面体粒子堆积 体系的边长，且为水泥净楽的最大孔径。
[0023] 步骤4、对于平面Si，利用单元格的九个点的坐标和不在Si上的其他凸形多面体粒 子的顶点坐标，判断该单元格与凸形多面体粒子的关系。判别过程如下：
[0024] 将点V和不在Si上的其他凸形多面体粒子的顶点P4, P5，…，Pk的坐标分别带入 Si的平面方程，观察所得到的结果是否同号。
[0025] 若同号，则说明V点与不在Si上的其他凸形多面体粒子的顶点均在平面Si的同一 侦k对于凸形多面体粒子的任意一平面，若V点与不在该平面上的其他凸形多面体粒子的 顶点均在该平面的同侧，则说明V点在该凸形多面体粒子内部。
[0026] 对于正四面体凸形粒子，需判定V点与其四个表面的关系，同理，对于正六面体凸 形粒子，需判定V点与其六个表面的关系；对于正八面体凸形粒子，需判定V点与其八个表 面的关系；对于正十二面体凸形粒子，需判定V点与其十二个表面的关系；对于正二十面体 凸形粒子，需判定V点与其二十个表面的关系；
[0027] 若单元格的九个顶点全部符合上述条件，则说明该单元格在凸形多面体粒子的内 部，称这种单兀格为固相单兀格；
[0028] 若单元格的九个顶点只有部分顶点满足上述条件，其他顶点不满足条件，则说明 该单元格只有部分空间在凸形多面体粒子内部，称这种单元格为混合相单元格，混合相单 兀格是研究孔径分布的关键；
[0029] 若单元格的全部顶点均不满足上述方程，则说明该单元格在该凸形多面体粒子外 部。若该单元格在三维体系内所有凸形多面体粒子外部，则称这种单元格为孔相单元格。
[0030] 步骤5、选取其中所有的混合相单元格，将混合相单元格按照边长d2进行划分，一 般取d 2 = 0. δΑ，不断重复步骤4将这个三维连续结构体系进行划分，最终可得出孔隙的半 径。
【专利附图】
【附图说明】
[0031] 图1为正四面体、正六面体、正八面体、正十二面体和正二十面体的堆积体系结 构；
[0032] 图2为在设定坐标系下，凸形多面体粒子（正四面体）的堆积体系结构；
[0033] 图3为正四面体粒子的几何形状与表示方法；
[0034] 图4为将三维体系进行单元格的划分；
[0035] 图5为固相单元格；
[0036] 图6为混合相单元格；
[0037] 图7为孔相单元格；
[0038] 图8为混合相单元格的划分。
【具体实施方式】
[0039] 实施例：本发明的计算凸形多面体粒子堆积体系孔径分布的方法，具体包括以下 步骤：
[0040] 步骤1、水泥水化后，内部结构非常复杂，包含大量未水化的水泥颗粒，内部和外 部的水化产物（参见文献[刘琳.静荷载和冰冻荷载耦合作用下水泥基材料的劣化研究 [D].东南大学，南京：东南大学，2012])。这些固体之间相互交错连接，中间形成大量的 孔隙。为了研究这些孔隙的大小，需要将固体和孔隙进行简化，即将水泥浆体简化为规则的 三维立体体系，而这些固体简化为凸形多面体粒子，杂乱无序的排列在三维体系中，而凸形 多面体粒子间隙即为所要研究的孔相，其大小即为孔径。将水泥浆体假定为三维凸形多面 体粒子堆积体系后，其微观结构是由固相（包括未水化水泥颗粒和水化产物）和毛细孔相 组成的两相结构体系，包括由如图1所示的正四面体粒子、正六面体粒子、正八面体粒子、 正十二面体粒子、正二十面体粒子组成的堆积体系；
[0041] 步骤2、如图2和图3所示，计算正四面体粒子随机堆积体系的孔径分布。在全局 坐标下计算出每个正四面体粒子的顶点坐标和每个平面的空间方程，假定平面Si的三个凸 多面体顶点的坐标为Pi (Xu yp Z)，P2 (x2, y2, z2)，P3(x3, y3, z3)，则Si的平面方程如下：
【权利要求】
1. 一种计算凸形多面体粒子堆积体系孔径分布的方法，通过将混合相单元格不断地进 行划分最后计算求得孔径，其特征在于包括以下步骤： 步骤1、将水泥浆体假定为三维凸形多面体粒子堆积体系，其微观结构是由固相和毛细 孔相组成的两相结构体系； 步骤2、在所述三维凸形多面体粒子堆积体系中确定每一个凸形多面体的顶点坐标和 每一个面的空间方程； 步骤3、在整体坐标系中，将三维凸形多面体粒子堆积体系按照边长屯进行单元划分， 并且确定出每一个单元格的η个合适点的坐标deL/UN取大于等于2的整数，L为所述三 维凸形多面体粒子堆积体系的边长，且屯< 2Mm，2Mm为水泥净浆的最大孔径； 步骤4、利用单元格η个点的坐标、凸形多面体的顶点坐标和每个平面的空间方程，判 断该单元格与凸形多面体粒子的关系，从而得出单元格的种类； 步骤5、将其中的混合相单元格按照边长d2进行划分，d^di/N'，Ν'取大于等于2的整 数，不断将这个三维连续结构体系进行划分，最终可得出全部孔隙的直径。
2. 根据权利要求1所述的计算凸形多面体粒子堆积体系孔径分布的方法，其特征在 于：所述微观结构的固相为凸形多面体粒子，所述凸形多面体粒子包括正四面体粒子、正六 面体粒子、正八面体粒子、正十二面体粒子和正二十面体粒子。
3. 根据权利要求1所述的计算凸形多面体粒子堆积体系孔径分布的方法，其特征在 于：将所包含的固相和毛细孔相的三维体系全部划分。
4. 根据权利要求1所述的计算凸形多面体粒子堆积体系孔径分布的方法，其特征在 于：所述η个合适点的坐标包括8个顶点坐标和1个中心点坐标，利用每个单元格9个点的 坐标、凸形多面体粒子的顶点坐标和每个平面的空间方程，判断单元格与凸形多面体粒子 的关系。
5. 根据权利要求1所述的计算椭球粒子堆积体系孔径分布的方法，其特征在于：所述 单元格的种类包含：固相、孔相和混合相。
6. 根据权利要求1所述的计算凸形多面体粒子堆积体系孔径分布的方法，其特征在 于：仅将其中的混合相单元格不断进行划分，最终得到全部孔径的大小。
7. 根据权利要求1所述的计算凸形多面体粒子堆积体系孔径分布的方法，其特征在 于：取Ν' =2,即dfcV〗。
【文档编号】G01N15/08GK104089864SQ201410307872
【公开日】2014年10月8日 申请日期:2014年6月30日 优先权日:2014年6月30日
【发明者】刘琳, 赵晓光, 沈德建 申请人:河海大学