CN104091059B - 一种确定椭球粒子堆积体系孔径分布的方法 - Google Patents
一种确定椭球粒子堆积体系孔径分布的方法 Download PDFInfo
- Publication number
- CN104091059B CN104091059B CN201410307006.3A CN201410307006A CN104091059B CN 104091059 B CN104091059 B CN 104091059B CN 201410307006 A CN201410307006 A CN 201410307006A CN 104091059 B CN104091059 B CN 104091059B
- Authority
- CN
- China
- Prior art keywords
- pore
- cell
- size distribution
- ellipsoidal particle
- phase
- 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.)
- Expired - Fee Related
Links
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Length-Measuring Devices Using Wave Or Particle Radiation (AREA)
Abstract
本发明公开了一种计算椭球粒子堆积体系孔径分布的方法,对水泥浆体的内部结构进行合理的假设,形成三维空间椭球粒子堆积体系,椭球粒子堆积间隙即为内部孔隙结构。将体系进行单元划分,判断每个单元格的n个点与椭球粒子的包含关系,从而得出该单元格的分类,单元格的种类有固相、孔相和混合相。对其中的混合相单元格进行更小的单元划分,最终得到全部孔径的大小。通过本发明计算得到的模拟孔径分布准确,并且计算速度快。
Description
技术领域
本发明涉及到一种确定椭球粒子堆积体系孔径分布的方法。
背景技术
水泥浆体的孔隙结构是关系到土木工程材料的渗透性,强度和耐久性的关键因素。当材料受到外部荷载,冻融循环和各种盐类侵蚀时,材料的内部结构会发生严重变化。为了提高材料的各项性能,使材料在受到损伤时具有足够的抵抗能力,需要对材料的内部结构有一定深入的了解,孔径分布是材料结构的重要参数之一。
通过压汞法,氮气吸脱附法,微观扫描电镜和热孔计法等试验可以得到水泥浆体内部的孔径分布情况,但是这些试验结果往往受制于试验条件和试验手段,且试验过程中会对水泥浆体内部的孔隙结构有破坏作用。为把人力和资源从大量试验中解脱出来,从本质上把握材料的结构性能变化规律,用于模拟材料微细观结构的计算机模型研究成为当前的一个研究热点。水泥水化后,内部结构非常复杂,包含大量未水化的水泥颗粒,内部和外部的水化产物。这些固体之间相互交错连接,中间形成大量的孔隙。为了研究这些孔隙的大小,需要将固体或孔隙进行简化,即将水泥浆体简化为规则的三维立体体系,而这些固体简化为球形或椭球粒子,杂乱无序的排列在三维体系中,而球形或椭球粒子间隙即为所要研究的孔相,其大小即为孔径。基于球形粒子或椭球形粒子的蒙特卡罗法堆积及相应的水化反应,得到混凝土的模拟微观及细观结构【van Breugel K.Simulation of hydration andformation of structure in hardening cement-based materials.Ph.D thesis,DelftUniversity of Technology,Delft,The Netherlands,1991.】【Xu W.X.,Chen H.S.,LvZ.An overlapping detection algorithm for random sequential packing ofelliptical particles.Physica A,390(2011)2452-67.】。对于孔径分布的模型分析,已提出的有丢球算法,通过向模拟结构中依次丢入由大到小的球体,判断球体是否能够放入来计算模拟结构的孔径分布【Ye G.Experimental study and numerical simulation ofthe development of the microstructure and permeability of cementitiousmaterials,PhD thesis,Delft University of Technology,Delft,2003.】。丢球法仅适用于球形粒子堆积模拟的微观结构,且计算量大,计算时间较长,尤其是丢入较小粒径的球体时。本发明的计算方法不仅可以计算球形粒子堆积的孔径分布,且适用于椭球形粒子堆积的模拟微观结构,而且相对已提出的算法而言,计算时间大大缩短,提高了计算效率。
发明内容
发明目的:本发明所要解决的技术问题在于开发一种克服原有孔径分布计算算法计算时间长,效率低,适用性窄等不足,提高计算孔径效率,并能够合理的计算出水泥浆体的孔径分布情况的方法。
技术方案:本发明提出一种确定椭球粒子堆积体系孔径分布的方法,通过将混合相单元格不断地进行划分最后计算求得孔径,包括以下步骤:
步骤1、将水泥浆体假定为三维椭球粒子堆积体系,其微观结构是由固相和毛细孔相组成的两相结构体系;
步骤2、设立整体坐标系O-XYZ和局部坐标系O’-X’Y’Z’,在局部坐标系中确定每一个椭球粒子的空间方程;
步骤3、在整体坐标系中,将三维椭球粒子堆积体系按照边长d1进行单元划分,并且确定出每一个单元格的n个合适点的坐标;d1=L/N,N取大于等于2的整数,L为所述三维椭球粒子堆积体系的边长,且d1≤2μm,2μm为水泥净浆的最大孔径;
步骤4、利用单元格n个点的坐标首先转换到局部坐标系中,然后带入椭球粒子的空间方程中,判断该单元格与椭球粒子的关系;
步骤5、将其中的混合相单元格按照边长d2进行划分,d2=d1/N′,N′取大于等于2的整数,不断将这个三维连续结构体系进行划分,最终可得出全部孔隙的直径。
作为优选,所述微观结构的固相为椭球粒子。
作为优选,将所包含的固相和毛细孔相的三维体系全部划分。
作为优选,n=9,包括8个顶点坐标和1个中心点坐标;将每个单元格的9个点的整体坐标转换为局部坐标后带入椭球粒子的空间方程,判断单元格与椭球粒子的关系。
作为优选,所述单元格的种类包含:固相、孔相和混合相
作为优选,仅将其中的混合相单元格不断进行划分,最终得到全部孔径的大小。
作为优选,取N′=2,即d2=d1/2。
使用时,具体操作如下:
步骤1、将水泥浆体假定为三维椭球粒子堆积体系,其微观结构是由固相(包括未水化水泥颗粒和水化产物)和毛细孔相组成的两相结构体系;
步骤2、设立整体坐标系O-XYZ和局部坐标系O’-X’Y’Z’,局部坐标系的原点在椭球粒子的中心,其中OX’沿着椭球粒子的长半轴a,OY’沿着椭球粒子的中半轴b,OZ’沿着椭球粒子的短半轴c,在局部坐标系中确定每一个椭球粒子的空间方程:
步骤3、在整体坐标系中按照边长d1进行单元划分,并且确定出每一个单元格八个顶点和一个中心的坐标V1(x1,y1,z1),V2(x2,y2,z2),…V9(x9,y9,z9)。
步骤4、将单元格的九个点的坐标转换为局部坐标后分别带入椭球粒子的空间方程中,判断该单元格与椭球粒子的关系。判别过程如下,对于点P(x’p,y’p,z’p):
若P点坐标满足上述不等式,则说明P点在椭球粒子内部或在椭球粒子上面。
若单元格的九个顶点全部符合上述方程,则说明该单元格在椭球粒子的内部,称这种单元格为固相单元格;
若单元格的九个顶点中至少有两个顶点满足上述方程,其他顶点不满足方程,则说明该单元格只有部分空间在椭球粒子内部,称这种单元格为混合相单元格,混合相单元格是研究孔径分布的关键;
若单元格的全部顶点均不满足上述方程,则说明该单元格在该椭球粒子外部,若该单元格在所有椭球粒子的外部,则称这种单元格为孔相单元格。
步骤5、选取其中所有的混合相单元格,将混合相单元格按照边长d2进行划分,一般取d2=0.5d1,不断将这个三维连续结构体系进行划分,最终可得出全部孔隙半径。
附图说明
图1为在设定坐标系下,椭球粒子的堆积体系结构;
图2为将三维体系进行单元格的划分;
图3为固相单元格;
图4为混合相单元格;
图5为孔相单元格;
图6为混合相单元格的划分;
图7为计算得到的水灰比0.4、水化程度0.69、宽径比0.8的椭球粒子堆积的水泥浆体孔径分布。
具体实施方式
实施例:本发明的计算椭球粒子堆积体系孔径分布的方法,具体包括以下步骤:
步骤1、将水泥浆体模拟为规则的边长为150μm的三维立方体,其中的未水化水泥及水化产物等固相用宽高比为0.8的椭球粒子表示,这些椭球粒子随机排列在三维体系中,而椭球粒子间隙即为所要研究的孔相,其大小即为孔径。模拟水泥浆体的水灰比为0.4,椭球粒子的等效直径Deq范围为1~50μm,粒子分布符合Rosin-Rammler分布:
其中n=1.698,b=0.04408。
椭球粒子等效直径Deq表示为:
其中,c为椭球短半轴,κ为宽径比,等于c/a。
水泥浆体的水化程度为0.69,水化过程表现为椭球粒子增大的过程,关于水泥水化模拟见文献【van Breugel K.Simulation of hydration and formation of structurein hardening cement-based materials.Ph.D thesis,Delft University ofTechnology,Delft,The Netherlands,1991.】。
步骤2、如图1所示,设立整体坐标系O-XYZ和局部坐标系O'-X'Y’Z’,局部坐标系的原点在椭球粒子的中心,其中OX’沿着椭球粒子的长半轴a,OY’沿着椭球粒子的中半轴b,OZ’沿着椭球粒子的短半轴c,在局部坐标系中确定每一个椭球粒子的空间方程:
步骤3、如图2所示,将三维体系划分成边长为d1的若干个小单元格,并且确定出每一个单元格八个顶点和一个中心的坐标V1(x1,y1,z1),V2(x2,y2,z2),…Vk(xk,yk,zk)。
步骤4、将每个单元格的九个点的坐标首先转换为局部坐标后分别带入椭球粒子的空间方程中,判断该单元格与椭球粒子的关系。判别过程如下,对于点P(x’p,y’p,z’p):
若单元格的九个顶点全部符合上述方程,则说明该单元格在椭球粒子的内部,称这种单元格为固相单元格,如图3所示;
若单元格的九个顶点至少有两个顶点满足上述方程,其他顶点不满足方程上,则说明该单元格只有部分空间在椭球粒子内部,称这种单元格为混合相单元格,混合相单元格即是继续被划分的单元格,也是研究孔径分布的关键,如图4所示;
对于所有的椭球粒子,若单元格的全部顶点和中心点均不满足上述方程,则说明该单元格在所有椭球粒子外部,称这种单元格为孔相单元格,如图5所示。
步骤5、将混合相单元格单独研究,在已设定坐标系下,混合相内的固相空间方程不变,将混合相单元格划分为若干个边长为d2的单元格,一般取d2=d1/2,确定每个小单元格的八个顶点坐标和一个中心坐标,代入固相的空间方程,参照上述标准(步骤2至步骤5),判断每个小单元格的类别,继续将其中的混合相单元格进行划分,直至孔径小于4nm,如图6所示。
步骤6、计算得到的水灰比0.4、水化程度0.69、宽径比0.8的椭球粒子堆积的水泥浆体孔径分布如图7所示。
Claims (7)
1.一种确定椭球粒子堆积体系孔径分布的方法,通过将混合相单元格不断地进行划分最后计算求得孔径,其特征在于包括以下步骤:
步骤1、将水泥浆体假定为三维椭球粒子堆积体系,其微观结构是由固相和毛细孔相组成的两相结构体系;
步骤2、设立整体坐标系O-XYZ和局部坐标系O’-X’Y’Z’,在局部坐标系中确定每一个椭球粒子的空间方程;
步骤3、在整体坐标系中,将三维椭球粒子堆积体系按照边长d1进行单元划分,并且确定出每一个单元格的n个合适点的坐标;n=9,包括8个顶点坐标和1个中心点坐标;d1=L/N,N取大于等于2的整数,L为所述三维椭球粒子堆积体系的边长,且d1≤2µm,2µm为水泥净浆的最大孔径;
步骤4、利用单元格n个点的坐标首先转换到局部坐标系中,然后带入椭球粒子的空间方程中,判断该单元格与椭球粒子的关系,从而得出单元格的种类;
步骤5、将其中的混合相单元格按照边长d2进行划分,d2= d1/N’,N’取大于等于2的整数,不断将这个三维连续结构体系进行划分,最终可得出全部孔隙的直径。
2.根据权利要求1所述的确定椭球粒子堆积体系孔径分布的方法,其特征在于:所述微观结构的固相为椭球粒子。
3.根据权利要求1所述的确定椭球粒子堆积体系孔径分布的方法,其特征在于:将所包含的固相和毛细孔相的三维体系全部划分。
4.根据权利要求1所述的确定椭球粒子堆积体系孔径分布的方法,其特征在于:将每个单元格的9个点的整体坐标转换为局部坐标后带入椭球粒子的空间方程,判断单元格与椭球粒子的关系。
5.根据权利要求1所述的确定椭球粒子堆积体系孔径分布的方法,其特征在于:所述单元格的种类包含:固相、孔相和混合相。
6.根据权利要求1所述的确定椭球粒子堆积体系孔径分布的方法,其特征在于:仅将其中的混合相单元格不断进行划分,最终得到全部孔径的大小。
7.根据权利要求1所述的确定椭球粒子堆积体系孔径分布的方法,其特征在于:取N’=2,即d2=d1/2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410307006.3A CN104091059B (zh) | 2014-06-30 | 2014-06-30 | 一种确定椭球粒子堆积体系孔径分布的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410307006.3A CN104091059B (zh) | 2014-06-30 | 2014-06-30 | 一种确定椭球粒子堆积体系孔径分布的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104091059A CN104091059A (zh) | 2014-10-08 |
CN104091059B true CN104091059B (zh) | 2017-01-18 |
Family
ID=51638775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410307006.3A Expired - Fee Related CN104091059B (zh) | 2014-06-30 | 2014-06-30 | 一种确定椭球粒子堆积体系孔径分布的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104091059B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104732096A (zh) * | 2015-03-31 | 2015-06-24 | 河海大学 | 非均质各向异性硬化粒子周围界面体积分数的计算方法 |
CN109841266B (zh) * | 2018-12-29 | 2021-03-16 | 河海大学 | 高吸水树脂内养护水泥浆体早期孔径模型构建方法及应用 |
CN110806371B (zh) * | 2019-11-20 | 2021-04-20 | 成都理工大学 | 单重煤岩孔隙结构模型的构建方法 |
CN115808376B (zh) * | 2022-12-12 | 2024-05-07 | 滨州学院 | 一种激光熔覆粉末流聚集性测量方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103440409A (zh) * | 2013-08-14 | 2013-12-11 | 河海大学 | 一种水泥基复合材料离子扩散模型的构建方法 |
-
2014
- 2014-06-30 CN CN201410307006.3A patent/CN104091059B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103440409A (zh) * | 2013-08-14 | 2013-12-11 | 河海大学 | 一种水泥基复合材料离子扩散模型的构建方法 |
Non-Patent Citations (5)
Title |
---|
Evaluation of mesostructure of particulate composites by quantitative stereology and random sequential packing model of mono-polydisperse convex polyhedral particles;W.XU, ET AL.,;《Industrial & Engineering Chemistry Research》;20130531;第52卷(第20期);6678-6693 * |
Micrometer-scale 3-D shape characterization of eight cements:Particle shape and cement chemistry,and the effect of particle shape on laser diffraction particle size measurement;S.T.Erdogan,et al.,;《Cement and Concrete Research》;20101231;第40卷;731-739 * |
Microstructural characterization of fresh cement paste via random packing of ellipsoidal cement particles;W.X.Xu,et al.,;《Materials characterization》;20121231(第66期);16-23 * |
人工模拟堆石料颗粒破碎应变的三轴试验研究;孔德志,等;《岩土工程学报》;20090331;第31卷(第3期);464-469 * |
大气环境混凝土中锈蚀钢筋性能和锈蚀模型试验研究;沈德建,等;《混凝土与水泥制品》;20040630(第3期);46-49 * |
Also Published As
Publication number | Publication date |
---|---|
CN104091059A (zh) | 2014-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104091059B (zh) | 一种确定椭球粒子堆积体系孔径分布的方法 | |
Sok et al. | Direct and stochastic generation of network models from tomographic images; effect of topology on residual saturations | |
Jivkov et al. | A novel architecture for pore network modelling with applications to permeability of porous media | |
CN105139444B (zh) | 基于岩心二维颗粒图像的三维颗粒结构重建方法 | |
Thovert et al. | Percolation in three-dimensional fracture networks for arbitrary size and shape distributions | |
US9589081B2 (en) | Karstification simulation | |
CN102521485A (zh) | 一种利用dda对粗粒土工程性质的数值仿真算法 | |
CN105425280A (zh) | 一种矿物模量和孔隙结构预测方法 | |
Sheintuch et al. | Deterministic approaches to problems of diffusion, reaction and adsorption in a fractal porous catalyst | |
Pardo-Igúzquiza et al. | A review of fractals in karst | |
CN104089864B (zh) | 一种计算凸形多面体粒子堆积体系孔径分布的方法 | |
CN105653746A (zh) | 基于嵌入式离散裂缝模型的压裂井建模及模拟方法 | |
Liu et al. | Bimodal and multimodal descriptions of soil-water characteristic curves for structural soils | |
CN111077053A (zh) | 一种裂缝型多孔介质渗流分析可视化方法 | |
Ko et al. | Monte Carlo simulation of concentrated diblock copolymers in a selective solvent: anisotropy of the diffusion | |
Miracle et al. | Clustering in interfering binary mixtures | |
Liu et al. | Continuous characterization of pore structure in cement paste | |
Hinsen et al. | Malware propagation in urban D2D networks | |
Heama et al. | Comparative effects of adjacent loaded pile row on existing tunnel by 2D and 3D simulation models | |
CN114544286A (zh) | 一种基于粒径分布的颗粒随机生成方法 | |
Imamura et al. | A practical method of nation wide Vs30 mapping with 250m grids for developing countries by automated topographical classification and published data in Japan | |
Das et al. | Quantification of quartz reefs and mafic dykes of Bundelkhand, Craton, Central India: a study based on spatial and fractal analysis | |
Xu | Unsaturated hydraulic conductivity of fractal-textured soils | |
CN114143743B (zh) | 一种基于规则网格的三维传感器部署方法和系统 | |
Patel et al. | Can a reliable prediction of cement paste transport properties be made using microstructure models |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170118 Termination date: 20190630 |
|
CF01 | Termination of patent right due to non-payment of annual fee |