CN116754587A - 三维层面的再生骨料表面残余砂浆层分布的量化表征方法 - Google Patents
三维层面的再生骨料表面残余砂浆层分布的量化表征方法 Download PDFInfo
- Publication number
- CN116754587A CN116754587A CN202310742620.1A CN202310742620A CN116754587A CN 116754587 A CN116754587 A CN 116754587A CN 202310742620 A CN202310742620 A CN 202310742620A CN 116754587 A CN116754587 A CN 116754587A
- Authority
- CN
- China
- Prior art keywords
- recycled aggregate
- particles
- natural
- residual
- formula
- 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.)
- Pending
Links
- 239000004570 mortar (masonry) Substances 0.000 title claims abstract description 32
- 238000012512 characterization method Methods 0.000 title claims abstract description 16
- 238000009826 distribution Methods 0.000 title claims abstract description 13
- 239000002245 particle Substances 0.000 claims abstract description 165
- 238000000034 method Methods 0.000 claims abstract description 46
- 239000002002 slurry Substances 0.000 claims abstract description 39
- 230000011218 segmentation Effects 0.000 claims abstract description 8
- 229910001385 heavy metal Inorganic materials 0.000 claims abstract description 7
- 238000012856 packing Methods 0.000 claims abstract description 7
- 238000002591 computed tomography Methods 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 7
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 239000002612 dispersion medium Substances 0.000 abstract 1
- 239000010410 layer Substances 0.000 description 17
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 11
- 239000000741 silica gel Substances 0.000 description 9
- 229910002027 silica gel Inorganic materials 0.000 description 9
- AIYUHDOJVYHVIT-UHFFFAOYSA-M caesium chloride Chemical compound [Cl-].[Cs+] AIYUHDOJVYHVIT-UHFFFAOYSA-M 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 238000007581 slurry coating method Methods 0.000 description 6
- 239000004567 concrete Substances 0.000 description 5
- 230000010354 integration Effects 0.000 description 5
- 150000002500 ions Chemical class 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 239000002699 waste material Substances 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 239000002253 acid Substances 0.000 description 4
- 238000010276 construction Methods 0.000 description 4
- 238000004090 dissolution Methods 0.000 description 4
- GAISRYMZAXLPHD-UHFFFAOYSA-N 2-hydroxybenzoic acid;methanol Chemical compound OC.OC(=O)C1=CC=CC=C1O GAISRYMZAXLPHD-UHFFFAOYSA-N 0.000 description 3
- OKKJLVBELUTLKV-UHFFFAOYSA-N Methanol Chemical compound OC OKKJLVBELUTLKV-UHFFFAOYSA-N 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- VEXZGXHMUGYJMC-UHFFFAOYSA-N Hydrochloric acid Chemical compound Cl VEXZGXHMUGYJMC-UHFFFAOYSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000010426 asphalt Substances 0.000 description 2
- 239000004568 cement Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000011978 dissolution method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- YGSDEFSMJLZEOE-UHFFFAOYSA-N salicylic acid Chemical compound OC(=O)C1=CC=CC=C1O YGSDEFSMJLZEOE-UHFFFAOYSA-N 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- VEXZGXHMUGYJMC-UHFFFAOYSA-M Chloride anion Chemical compound [Cl-] VEXZGXHMUGYJMC-UHFFFAOYSA-M 0.000 description 1
- 238000003763 carbonization Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 239000011247 coating layer Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000009440 infrastructure construction Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000011326 mechanical measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- FJKROLUGYXJWQN-UHFFFAOYSA-N papa-hydroxy-benzoic acid Natural products OC(=O)C1=CC=C(O)C=C1 FJKROLUGYXJWQN-UHFFFAOYSA-N 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 229960004889 salicylic acid Drugs 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000002791 soaking Methods 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/401—Imaging image processing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/421—Imaging digitised image, analysed in real time (recognition algorithms)
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种三维层面的再生骨料表面残余砂浆层分布的量化表征方法，包括：1采用重金属离子溶液浸泡再生骨料；2采用分散介质将浸泡后的再生骨料固定于塑料管内；3采用工业CT系统对再生骨料样品进行断面扫描；4基于数字图像技术得到区分天然颗粒和残余砂浆的灰度阈值；5通过阈值分割法得到总再生骨料和其中天然颗粒的表面体元的实际空间坐标；6基于球谐函数对总再生骨料和其中天然颗粒进行重构，并基于此计算残余浆体包裹体积比和残余浆体包裹表面积比等量化参数。本发明能准确方便地定量再生骨料残余砂浆，实现对再生骨料特征的精确描述，从而有利于再生骨料的合理和高效应用。
Description
技术领域
本发明属于再生骨料混凝土领域，具体涉及一种在三维层面上的再生骨料表面残余砂浆层分布的量化表征方法。
背景技术
随着我国城市化进行的加快发展，建筑垃圾排放量逐年增长，建筑垃圾中可再生的组成比例也不断的提高。大部分的建筑垃圾没有经过任何处理，直接运往郊外或城市周边进行简单填埋或露天堆放，这不仅浪费了土地和资源，而且污染了环境。另一方面，随着建筑业和基础设施建设对砂石材料的需求居高不下，长期的砂石原材料开采造成资源枯竭、山体滑坡、河床改道，严重破坏了自然环境。生产和再利用再生骨料对于节约资源、保护环境和实现建筑业的可持续发展具有重要意义
废弃混凝土仅通过简单的破碎、筛分工艺制备的再生骨料表面会残留硬化的水泥浆体，从而导致再生骨料存在空隙率大、吸水率大、堆积密度小、压碎值指标高等缺点。采用此类再生骨料制备的再生混凝土存在多重过渡区界面，严重影响其硬化强度、弹性模量、抗渗性、抗冻性、抗碳化能力和抗氯离子渗透等性能的主要影响因素。确定再生骨料表面浆体的含量及其分布情况，定量化描述再生骨料特点，对再生骨料的应用十分重要和必要。现有方法通常采用酸溶处理，仅能测得残余砂浆的含量，无法测定砂浆层的厚度、包裹位置等颗粒特征；鉴于此，CN110823060A采用沥青包裹再生骨料，利用残余砂浆与再生骨料、沥青胶浆的纳米力学参数的差异，利用原位纳米力学测量系统进行硬度测试，从而确定特征位置处的残余砂浆厚度，但此方法对制样打磨精度要求高，且仅能实现二维局部厚度的测试，不能准确表征真实再生骨料表面残余浆体的特征。截止到目前，尚没有在三维层面针对再生骨料表面残余浆体含量及其分布情况的量化评价方法。
发明内容
本发明是为了解决上述现有技术存在的不足之处，提出一种三维层面的再生骨料表面残余砂浆层分布的量化表征方法，以期能准确方便地定量再生骨料残余砂浆，实现对再生骨料特征的精确描述，从而有利于再生骨料的合理和高效应用。
本发明为达到上述发明目的，采用如下技术方案：
本发明一种三维层面的再生骨料表面残余砂浆层分布的量化表征方法的特点在于，包含以下步骤：
步骤1：采用重金属离子溶液浸泡再生骨料后，得到用于CT扫描的再生骨料颗粒；
步骤2：采用分散介质将浸泡后的再生骨料分散地固定于圆形塑料管内，得到用于CT扫描的再生骨料样品；
步骤3：采用工业CT系统对所述再生骨料样品进行断面扫描，得到样品的断面灰度图序列；
步骤4：从样品的断面灰度图序列中提取单个再生骨料颗粒的灰度图序列，并采用最大类间方差法确定灰度图序列的三值化最佳阈值，作为区分单个再生骨料颗粒中天然颗粒和残余砂浆的灰度阈值；
步骤5：根据灰度阈值，采用阈值分割法对所述灰度图序列进行处理，分别得到单个再生骨料颗粒的体素集合及单个再生骨料颗粒中的天然颗粒的体素集合；
以灰度图序列中的每张灰度图的左下角为二维笛卡尔坐标系的原点，以与所述原点相邻的两条边分别为x轴和y轴，并以厚度方向为z轴正方向，建立全局三维笛卡尔坐标系；
识别对应体素集合的表面体元，并根据表面体元的像素位置和像素分辨率分别得到在所述全局三维笛卡尔坐标系中总再生骨料表面体元的实际空间坐标以及所述单个再生骨料颗粒中天然颗粒表面体元的实际空间坐标；
步骤6：将所述全局三维笛卡尔坐标系的原点移至单个再生骨料颗粒中天然颗粒的几何中心位置，得到局部三维笛卡尔坐标系；
计算在所述局部三维笛卡尔坐标系中单个再生骨料颗粒表面体元的坐标数据以及所述单个再生骨料颗粒中天然颗粒表面体元的坐标数据，并分别换算为球坐标系下的坐标数据；且所述球坐标系下的原点与局部三维笛卡尔坐标系的原点重合；
基于球谐函数对单个再生骨料颗粒和其中天然颗粒的轮廓进行重构，从而利用式(1)和式(2)分别计算残余浆体包裹体积比RV和残余浆体包裹表面积比RS；
式(1)和式(2)中，VRCA为单个再生骨料颗粒的体积，VNA为单个再生骨料颗粒中天然颗粒的体积；Sw-NA为单个再生骨料颗粒中天然颗粒被残余浆体包裹的表面积，St-NA为单个再生骨料颗粒中天然颗粒的总表面积。
本发明所述的量化表征方法的特点也在于，所述步骤6中是分别采用式(3)和式(4)分别对单个再生骨料颗粒和其中天然颗粒的轮廓进行球谐重构：
式(3)和式(4)中，N1为球谐展开阶数；为球谐函数的n阶m次基，anm和bnm分别为单个再生骨料颗粒和其中天然骨料的轮廓对应于/>的球谐系数；/>和分别代表单个再生骨料颗粒和其中天然颗粒的表面体元在极角为θ和方位角为/>时的半径值；/>和/>分别为单个再生骨料颗粒和其中天然颗粒表面轮廓的重构解析式。
所述步骤6中是采用式(4)和式(5)分别计算单个再生骨料颗粒的体积VRCA和其中天然颗粒的体积VNA：
所述步骤6中是残余浆体包裹表面积Sw-NA和总表面积St-NA是按如下步骤获得：
步骤a：分别将极角θ在[0,2π]和方位角在[0,π]范围内等分为N份，在单个再生骨料颗粒和其中天然颗粒的表面上均得到N2个表面点；
步骤b：采用Delaunay三角剖分法分别对上述两个轮廓表面的N2个表面点进行三角剖分，绘制三角面片后得到总再生骨料和其中天然颗粒的近似表面轮廓；
步骤c：分别对比在N2个角度下的和/>的值，若在某个角度下两者存在差值，则代表相应角度下单个再生骨料颗粒中的天然颗粒存在残余浆体包裹，反之，则代表无残余浆体包裹；
步骤d：累计计算单个再生骨料颗粒中天然颗粒表面所有三角面片的面积，从而得到再生骨料中天然颗粒的总表面积St-NA；
步骤e：删除单个再生骨料颗粒中天然颗粒上未被残余浆体包裹的三角面片，并累计计算剩余三角面片的面积之和，从而得到单个再生骨料颗粒中天然颗粒被残余浆体包裹的表面积Sw-NA。
本发明一种电子设备，包括存储器以及处理器的特点在于，所述存储器用于存储支持处理器执行所述量化表征方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
本发明一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序的特点在于，所述计算机程序被处理器运行时执行所述所述量化表征方法的步骤。
与现有技术相比，本发明的有益效果在于：
1、本发明采用重金属离子浸泡再生骨料以放大残余浆体和天然骨料的密度差异，然后利用工业CT扫描技术和数字图像技术提取得到再生骨料和其中天然骨料的表面体元球坐标数据；进一步基于球谐基函数对进行再生骨料和其中天然骨料的表面进行重构得到各表面的数学解析式，最后通过数值积分法和三角面片求和法可准确计算再生骨料中残余浆体包裹体积比RV和残余浆体包裹表面积比RS等参数的目的。此类量化指标能更准确全面地表征再生骨料表面残余砂浆的分布情况，对于分析不同破碎方式对再生骨料残余浆体剥离率的影响提供基础。
2、本发明相较于现有采用纳米力学技术在二维层面测量残余砂浆层的局部包裹厚度的方法，采用CT扫描、数字图像和球谐重构等技术，可进一步从三维层面识别并量化表征再生骨料中残余浆体的包裹情况，有助于在再生骨料混凝土离散元模拟时生成再生骨料几何模型。
3、本发明相较于现有的酸溶法测量残余砂浆含量的方法，通过对再生骨料和其中天然骨料表面的数学重构，采用数值积分法和三角面片求和法，可进一步量化表征再生骨料中残余砂浆分布特征，对建立残余浆体不同包裹参数和再生骨料混凝土性能间的关系提供基础。
附图说明
图1为本发明量化表征方法的流程图；
图2为本发明工业CT扫描后某层再生骨料的三维视图；
图3为本发明提取的单个颗粒的三维CT视图；
图4为本发明按最大类间方差法三值化结果示意图；
图5a为本发明阈值分割法得到的总再生骨料的形状示意图；
图5b为本发明阈值分割法得到总再生骨料中天然颗粒的形状示意图；
图6a为本发明经球谐重构后总再生骨料的形状示意图；
图6b为本发明经球谐重构后总再生骨料中天然颗粒的形状示意图；
图7为本发明三角面片求和法计算表面积的示意图。
具体实施方式
本实施例中，如图1所示，一种三维层面的再生骨料表面残余砂浆层分布的量化表征方法，包含以下步骤：
步骤1：由于残余砂浆吸水率远大于天然骨料，故采用CsCl等重金属离子溶液浸泡再生骨料，可达到放大再生骨料中残余砂浆和天然骨料的密度差异，得到用于CT扫描的再生骨料颗粒；其中采用的CsCl重金属离子溶液为的质量分数一般不低于30％；
步骤2：为防止再生骨料颗粒相互接触以及扫描过程中的颗粒晃动，需采用硅胶等分散介质将浸泡后的再生骨料分散地固定于圆形塑料管内，得到用于CT扫描再生骨料样品；其中采用硅胶类等低密度基质的原因是为了后期便于采用阈值分割法提取再生骨料相。
步骤3：采用工业CT系统对再生骨料样品进行断面扫描，得到样品的断面灰度图序列，如图2所示；为准确识别再生骨料表面的残余浆体包裹层，需控制样品尺寸或扫描参数以保证像素分辨率不大于30～40μm；
步骤4：从样品的断面灰度图序列中提取单个再生骨料颗粒的灰度图序列(如图3所示)，采用最大类间方差法(Otsu)确定灰度图序列的三值化最佳阈值，作为区分单个再生骨料颗粒中天然颗粒和残余砂浆的灰度阈值，图4显示了某再生颗粒某断面基于Otsu进行三值化的结果；
步骤5：基于阈值分割法获得单个再生骨料颗粒和其中天然颗粒表面体元的空间坐标：
步骤5-1：根据灰度阈值，采用阈值分割法对灰度图序列进行处理，分别得到单个再生骨料颗粒的体素集合及单个再生骨料颗粒中的天然颗粒的体素集合，如图5a和图5b所示；
步骤5-2：以灰度图序列中的每张灰度图的左下角为二维笛卡尔坐标系的原点，以与原点相邻的两条边分别为x轴和y轴，并以厚度方向为z轴正方向，建立全局三维笛卡尔坐标系；
步骤5-3：采用Matlab中的内置函数bwperim识别对应体素集合的表面体元，并根据表面体元的像素位置和像素分辨率分别得到在全局三维笛卡尔坐标系中总再生骨料表面体元的实际空间坐标以及单个再生骨料颗粒中天然颗粒表面体元的实际空间坐标；
步骤6：对单个再生骨料颗粒和其中天然颗粒表面体元坐标进行球谐重构和残余砂浆层分布参数的量化计算：
步骤6-1：为方便坐标数据向球坐标系下的转换，首先将全局三维笛卡尔坐标系的原点移至单个再生骨料颗粒中天然颗粒的几何中心位置，得到局部三维笛卡尔坐标系；
步骤6-2：计算在局部三维笛卡尔坐标系中单个再生骨料颗粒表面体元的坐标数据以及单个再生骨料颗粒中天然颗粒表面体元的坐标数据，并按照式(1)将局部坐标系下的数据分别换算到球坐标系下；且球坐标系下的原点与局部三维笛卡尔坐标系的原点重合；
式(1)中，θ为颗粒中心到各表面体元的有向线段与x轴正方向的夹角，即极角；从正z轴看逆时针为正；为颗粒中心到各表面体元的有向线段与与z轴正向的夹角，即方位角；/> 为极角为颗粒表面体元的到质心的距离，即半径。x,y和z分别为颗粒表面体元在局部三维笛卡尔坐标系中的坐标值。
步骤6-3：按照式(2)和式(3)，基于球谐函数对单个再生骨料颗粒和其中天然颗粒的轮廓进行重构，重构后的轮廓形状如图6a和图6b所示；
式(2)和式(3)中，N1为球谐展开阶数，应不低于15以保证重构精度；为球谐函数的n阶m次基，anm和bnm分别为总再生骨料和其中天然骨料轮廓对应于/>的球谐系数，其计算可采用Gauss积分法进行计算，其中在极角和方位角上的高斯积分点数均应不少于120个，以保证积分的精确度；/>和/>分别代表总再生骨料和其中天然颗粒表面体元在极角为θ和方位角为/>时的半径值；/>和/>分别为总再生骨料和其中天然颗粒表面的重构解析式。
步骤6-4：分别采用式(4)和式(5)分别计算总再生骨料的体积VRCA和其中天然颗粒的体积VNA：
步骤6-5：计算单个再生骨料中天然颗粒的总表面积St-NA和残余浆体包裹表面积Sw-NA。首先分别将极角θ在[0,2π]和方位角在[0,π]范围内等分为N份，其中，为保证三角面片求和法计算表面积的精确度，极角和方位角的等分份数N均不应小于256份，在单个再生骨料颗粒和其中天然颗粒的表面上均得到N2个表面点；
采用Delaunay三角剖分法分别对上述两个轮廓表面的N2个表面点进行三角剖分，绘制三角面片后得到总再生骨料和其中天然颗粒的近似表面轮廓；
分别对比在N2个角度下的和/>的值，考虑到重构精度的影响，可以设置一个截断误差。若在某角度下两者存在的差值大于此阶段误差，则代表相应角度下单个再生骨料颗粒中的天然颗粒存在残余浆体包裹，反之，则代表无残余浆体包裹；
累计计算单个再生骨料颗粒中天然颗粒表面所有三角面片的面积，从而得到再生骨料中天然颗粒的总表面积St-NA；
删除单个再生骨料颗粒中天然颗粒上未被残余浆体包裹的三角面片，并累计计算剩余三角面片的面积之和，从而得到单个再生骨料颗粒中天然颗粒被残余浆体包裹的表面积Sw-NA；其中各三角面片面积的计算原理如图7和式(6)所示。
式(6)中，Si为第i个三角面片的面积；pi1,pi2和pi3分别代表第i个三角面片的三个顶点；θi代表第i个三角面片中向量和/>的夹角。
步骤6-6：利用式(7)和式(8)分别计算残余浆体包裹体积比RV和残余浆体包裹表面积比RS；
式(7)和式(8)中，VRCA为总再生骨料(Recycled Aggregate)的体积，VNA为再生骨料中天然颗粒(Natural Aggregate)的体积；Sw-NA为再生骨料中天然颗粒被残余浆体包裹的表面积，St-NA为再生骨料中天然颗粒的总表面积。
步骤7：将再生骨料颗粒从硅胶等分散基质中剥离出来，烘干后使用研钵进行捣碎，由于钙质骨料也会被盐酸溶解掉一部分，所以采用水杨酸甲醇溶液溶解再生骨料中的残余浆体，采用量筒测定酸溶前后再生骨料的体积变化可得到再生骨料残余浆体包裹体积比RV，以验证基于CT扫描技术和球谐函数计算的残余浆体包裹体积比RV的准确性。
实施例：采用以下步骤测量和计算再生骨料残余砂浆特征参数：
步骤1：重金属离子浸泡再生骨料：
将15g的CsCl粉末溶于35g水中，制备质量分数为30％的CsCl溶液，取数颗粒径在2.36～4.75mm范围内的再生细骨料置于其中，浸泡12h；
步骤2：硅胶基质分散固定再生骨料：
将加成型有机硅胶分层倒入内径为30mm的塑料圆管内，待其固化强度达到可支撑再生颗粒不沉降时，将再生骨料分散地摆放在有机硅胶基质内，需确保各颗粒相互不接触，最后继续倒入硅胶没过颗粒；
步骤3：工业CT断面扫描试样：
待硅胶固化后，将装有再生骨料的塑料圆管进行CT断层扫描，获取的图像分辨率为1980×1980×253，平面及轴向扫描切片图像层间距为0.0197mm；
步骤4：提取单个颗粒、确定区分浆体和颗粒的灰度阈值：
提取单个颗粒的灰度图序列，如图3所示，采用最大类间方差法确定三值化最佳阈值为68，作为区分天然颗粒和残余砂浆的灰度阈值；
步骤5：获得总再生骨料颗粒及其中天然颗粒的表面体元坐标：
设置不同阈值，分别得到总再生骨料和其中天然颗粒的体素集合，统计各自的体素个数乘以图像像素大小的立方即可得到颗粒的扫描体积，VRCA和VNA分别为40.69和36.56mm3；借助Matlab中的bwperim内置函数即可得到总再生骨料和其中天然颗粒的表面体元的像素坐标，乘以图像分辨率即可得到实际空间坐标；
步骤6：基于球谐基函数进行轮廓重构以及参数计算：
将三维笛卡尔坐标系的原点平移至单个再生骨料中天然颗粒区域中心，并将平移后的表面体元坐标转化到球坐标系下，接着按式(2)和式(3)将单个再生骨料和其中天然颗粒进行球谐重构，其中采用15阶球谐函数进行重构。
基于球谐重构后的表面轮廓的数学解析式，按式(4)和式(5)计算总再生骨料VRCA和其中天然颗粒的体积VNA，分别为39.85mm3和36.98mm3，和CT结果的体素叠加法吻合良好。
分别将极角和方位角等分为256份，基于Delaunay三角剖分法绘制各表面三角面片，累加再生骨料中天然颗粒的各三角面片面积即可得到其总表面积St-NA为68.23mm2。
遍历N2个角度，设置截断误差为0.005mm，对比各角度下和/>的差值以判断此角度下是否存在包裹；删除再生骨料中天然颗粒上未被残余浆体包裹的三角面片，累计计算剩余三角面片的面积和，即可得到再生骨料中天然颗粒被残余浆体包裹的表面积Sw-NA为15.32mm2；
进一步按式(7)计算得到再生骨料残余浆体包裹体积比RV为7.76％；按式(8)计算得到再生骨料残余浆体包裹表面积比RS为22.06％。
步骤7：基于酸溶法对计算结果进行验证：
将CT扫描的再生骨料颗粒从硅胶基质中剥离出来，选择提取的再生骨料，首先利用量程为0.2ml的高精度量筒(测量精度可达0.003ml)基于排水法测量其体积，结果为39mm3；然后将再生颗粒在研钵中进行捣碎，取14g水杨酸溶于80ml无水甲醇中制成水杨酸甲醇溶液，将捣碎后的再生颗粒倒入上述水杨酸甲醇溶液中溶解再生骨料中的残余浆体，待溶解2h后过滤并测量剩余体积，结果为36mm3；根据溶解前后的体积变化计算得到再生骨料残余浆体体积比RV为7.69％，和基于CT和球谐重构方法计算得到的RV相差较小，即可验证基于CT技术和球谐函数计算的残余浆体包裹体积比RV的准确性。
本实施例中，一种电子设备，包括存储器以及处理器，该存储器用于存储支持处理器执行上述方法的程序，该处理器被配置为用于执行该存储器中存储的程序。
本实施例中，一种计算机可读存储介质，是在计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时执行上述方法的步骤。
Claims (6)
1.一种三维层面的再生骨料表面残余砂浆层分布的量化表征方法，其特征在于，包含以下步骤：
步骤1：采用重金属离子溶液浸泡再生骨料后，得到用于CT扫描的再生骨料颗粒；
步骤2：采用分散介质将浸泡后的再生骨料分散地固定于圆形塑料管内，得到用于CT扫描的再生骨料样品；
步骤3：采用工业CT系统对所述再生骨料样品进行断面扫描，得到样品的断面灰度图序列；
步骤4：从样品的断面灰度图序列中提取单个再生骨料颗粒的灰度图序列，并采用最大类间方差法确定灰度图序列的三值化最佳阈值，作为区分单个再生骨料颗粒中天然颗粒和残余砂浆的灰度阈值；
步骤5：根据灰度阈值，采用阈值分割法对所述灰度图序列进行处理，分别得到单个再生骨料颗粒的体素集合及单个再生骨料颗粒中的天然颗粒的体素集合；
以灰度图序列中的每张灰度图的左下角为二维笛卡尔坐标系的原点，以与所述原点相邻的两条边分别为x轴和y轴，并以厚度方向为z轴正方向，建立全局三维笛卡尔坐标系；
识别对应体素集合的表面体元，并根据表面体元的像素位置和像素分辨率分别得到在所述全局三维笛卡尔坐标系中总再生骨料表面体元的实际空间坐标以及所述单个再生骨料颗粒中天然颗粒表面体元的实际空间坐标；
步骤6：将所述全局三维笛卡尔坐标系的原点移至单个再生骨料颗粒中天然颗粒的几何中心位置，得到局部三维笛卡尔坐标系；
计算在所述局部三维笛卡尔坐标系中单个再生骨料颗粒表面体元的坐标数据以及所述单个再生骨料颗粒中天然颗粒表面体元的坐标数据，并分别换算为球坐标系下的坐标数据；且所述球坐标系下的原点与局部三维笛卡尔坐标系的原点重合；
基于球谐函数对单个再生骨料颗粒和其中天然颗粒的轮廓进行重构，从而利用式(1)和式(2)分别计算残余浆体包裹体积比RV和残余浆体包裹表面积比RS；
式(1)和式(2)中，VRCA为单个再生骨料颗粒的体积，VNA为单个再生骨料颗粒中天然颗粒的体积；Sw-NA为单个再生骨料颗粒中天然颗粒被残余浆体包裹的表面积，St-NA为单个再生骨料颗粒中天然颗粒的总表面积。
2.根据权利要求1所述的量化表征方法，其特征在于，所述步骤6中是分别采用式(3)和式(4)分别对单个再生骨料颗粒和其中天然颗粒的轮廓进行球谐重构：
式(3)和式(4)中，N1为球谐展开阶数；为球谐函数的n阶m次基，anm和bnm分别为单个再生骨料颗粒和其中天然骨料的轮廓对应于/>的球谐系数；/>和/>分别代表单个再生骨料颗粒和其中天然颗粒的表面体元在极角为θ和方位角为/>时的半径值；/>和/>分别为单个再生骨料颗粒和其中天然颗粒表面轮廓的重构解析式。
3.根据权利要求2所述的量化表征方法，其特征在于，所述步骤6中是采用式(4)和式(5)分别计算单个再生骨料颗粒的体积VRCA和其中天然颗粒的体积VNA：
4.根据权利要求2所述的量化表征方法，其特征在于，所述步骤6中是残余浆体包裹表面积Sw-NA和总表面积St-NA是按如下步骤获得：
步骤a：分别将极角θ在[0,2π]和方位角在[0,π]范围内等分为N份，在单个再生骨料颗粒和其中天然颗粒的表面上均得到N2个表面点；
步骤b：采用Delaunay三角剖分法分别对上述两个轮廓表面的N2个表面点进行三角剖分，绘制三角面片后得到总再生骨料和其中天然颗粒的近似表面轮廓；
步骤c：分别对比在N2个角度下的和/>的值，若在某个角度下两者存在差值，则代表相应角度下单个再生骨料颗粒中的天然颗粒存在残余浆体包裹，反之，则代表无残余浆体包裹；
步骤d：累计计算单个再生骨料颗粒中天然颗粒表面所有三角面片的面积，从而得到再生骨料中天然颗粒的总表面积St-NA；
步骤e：删除单个再生骨料颗粒中天然颗粒上未被残余浆体包裹的三角面片，并累计计算剩余三角面片的面积之和，从而得到单个再生骨料颗粒中天然颗粒被残余浆体包裹的表面积Sw-NA。
5.一种电子设备，包括存储器以及处理器，其特征在于，所述存储器用于存储支持处理器执行权利要求1-4中任一所述量化表征方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
6.一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，其特征在于，所述计算机程序被处理器运行时执行权利要求1-4中任一所述所述量化表征方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310742620.1A CN116754587A (zh) | 2023-06-21 | 2023-06-21 | 三维层面的再生骨料表面残余砂浆层分布的量化表征方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310742620.1A CN116754587A (zh) | 2023-06-21 | 2023-06-21 | 三维层面的再生骨料表面残余砂浆层分布的量化表征方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116754587A true CN116754587A (zh) | 2023-09-15 |
Family
ID=87947522
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310742620.1A Pending CN116754587A (zh) | 2023-06-21 | 2023-06-21 | 三维层面的再生骨料表面残余砂浆层分布的量化表征方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116754587A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116663378A (zh) * | 2023-05-23 | 2023-08-29 | 合肥工业大学 | 一种考虑再生砂颗粒形貌的级配优化方法 |
-
2023
- 2023-06-21 CN CN202310742620.1A patent/CN116754587A/zh active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116663378A (zh) * | 2023-05-23 | 2023-08-29 | 合肥工业大学 | 一种考虑再生砂颗粒形貌的级配优化方法 |
| CN116663378B (zh) * | 2023-05-23 | 2024-06-04 | 合肥工业大学 | 一种考虑再生砂颗粒形貌的级配优化方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Coral aggregate concrete: Numerical description of physical, chemical and morphological properties of coral aggregate | |
| Erdogan et al. | Three-dimensional shape analysis of coarse aggregates: New techniques for and preliminary results on several different coarse aggregates and reference rocks | |
| Yio et al. | Representative elementary volume (REV) of cementitious materials from three-dimensional pore structure analysis | |
| Tian et al. | An investigation on the three-dimensional transport of chloride ions in concrete based on X-ray computed tomography technology | |
| Xu et al. | Microstructural characterization of fresh cement paste via random packing of ellipsoidal cement particles | |
| Koenig | Analysis of air voids in cementitious materials using micro X-ray computed tomography (µXCT) | |
| CN116754587A (zh) | 三维层面的再生骨料表面残余砂浆层分布的量化表征方法 | |
| CN110702536A (zh) | 一种岩体结构面剪切破坏表面形态演化规律研究方法及系统 | |
| CN102944505A (zh) | 一种定量表征水泥基材料孔隙率增加量空间分布的方法 | |
| Ren et al. | Quantitative analysis of relationships between crack characteristics and properties of soda-saline soils in Songnen Plain, China | |
| Hu et al. | Shape characterization of concrete aggregate | |
| WO2021087842A1 (zh) | 一种混凝土表面三维粗糙度的测量方法 | |
| Liu et al. | Magnification effects on the interpretation of SEM images of expansive soils | |
| Song et al. | A 3D petrographic analysis for concrete freeze-thaw protection | |
| CN110363855B (zh) | 堆石坝透明化建模方法 | |
| Deliang et al. | Random field model of corroded steel plate surface in neutral salt spray environment | |
| Liu et al. | Multiscale structural characterizations of anisotropic natural granite residual soil | |
| Ye et al. | Characterization of the development of microstructure and porosity of cement-based materials by numerical simulation and ESEM image analysis | |
| Kenny et al. | Characterization of the interfacial transition zone around steel rebar by means of the mean shift method | |
| Dong et al. | Review on the mesoscale characterization of cement-stabilized macadam materials | |
| CN110599477A (zh) | 一种植草混凝土有效孔隙率的检测方法及系统 | |
| Wang et al. | Microstructure and unsaturated geotechnical properties of net-like red soils in Xuancheng, China | |
| CN114283133A (zh) | 一种基于最邻近分析的沥青混合料均匀性评价方法 | |
| Zhang et al. | Study on the non-uniformity of rust layer and residual rebar caused by pitting corrosion in mortar | |
| Hu et al. | Ballast deterioration inspection and quantification with 3D form method based on particle inscribed ellipsoid |
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 |