logo

基于主成分分析的物理整形对再生粗骨料形貌特征影响的研究

轨道与基础

基于主成分分析的物理整形对再生粗骨料形貌特征影响的研究

景涛
新宇
昆林
友均
广成
明文
铁道科学与工程学报第22卷, 第3期pp.1112-1120纸质出版 2025-03-28
400

物理整形是再生粗骨料强化的重要方法之一。为了综合评价物理整形对再生粗骨料形态特征的影响,为再生粗骨料的制备及高品质利用提供研究基础,研究了物理整形对再生粗骨料物理性能的影响。通过集料图像测量系统(AIMS)采集不同整形时间的再生粗骨料典型形态特征参数,基于主成分分析方法将圆度、球度、细长度和扁平度组合成了整体形态因子,将棱角度、粗糙度和分形维数组合成表面形态因子,分析整形对再生粗骨料的形态特征的影响。研究结果表明:物理整形显著降低了再生粗骨料的空隙率、吸水率和压碎指标,改善了再生粗骨料性能。随着物理整形时间增加,再生粗骨料的整体形态因子增大,而表面形态因子先降低后增大。再生粗骨料的整体形态因子增大,则粗骨料的整体形状趋于球形,粗骨料的堆积密实程度越好,整体形态因子的增大是堆积空隙率降低的主要原因;表面形态因子降低,则粗骨料表面形态越规则。物理整形使再生粗骨料的整体形态较天然粗骨料更近似于球形,但不能使再生粗骨料表面形态比天然粗骨料更规则,且随整形时间延长,再生粗骨料可能发生二次破碎使颗粒表面更不规则,从而造成表面形态因子增大。通过物理整形强化再生粗骨料时,可根据粗骨料表面形态因子的变化来合理控制整形时间。

再生粗骨料物理整形形态特征主成分分析

在建筑固废资源化利用背景下,再生粗骨料作为混凝土的原材料,在工程中逐渐得到应用[1-2]。相比较天然碎石粗骨料,再生粗骨料的强度更低且形状更不规则,具有表面棱角复杂和质地粗糙等特点[3],通过物理整形可以改善再生粗骨料的强度和形态特征[4]。形态特征是粗骨料的基本特征,对混凝土的工作性、力学性能和耐久性能有着重要影响[5-8]。在实际过程应用中,对混凝土用天然粗骨料形态特征通常用针片状和扁平状含量表示[8],无法精确量化。随着数字图像处理技术和三维重构技术的发展,现有研究主要围绕形状、棱角度和粗糙度3个尺度表征粗骨料的形态特征[9-15],并就每个尺度提出大量计算方法,张雪芹等[8]总结了棱角度的计算方法有轮廓曲率、轮廓梯度、多边形拟合、椭球拟合和腐蚀膨胀法等近十种。除了上述3个尺度外,分形理论也被用于量化粗骨料的形态特征,常用方法包括面积-周长法和计盒维数法等[16-17]。可见,粗骨料形态特征指标的计算方法不统一且计算复杂,影响了形态特征指标的量化应用。主成分分析方法是一种从高维变量中提取相互独立的公共变量从而降低输入维度的一种多元统计学方法[18]。骆行文等[19]采用主成分分析方法对岩石的5个性能指标进行分析,最终得到了一个综合指标评价岩石的质量。肖柏林等[20]对7种矿山充填粗骨料的形状和棱角度进行评价,并通过主成分分析方法得到了2组新的指标。目前,还没有简单明确的指标能直接反映再生粗骨料的形态特征,对再生粗骨料形态特征的系统研究还需要进一步开展。因此,考虑再生粗骨料和天然粗骨料的形态特征差异,从众多形态特征指标中找出合适的指标表征再生粗骨料形态特征,这作为再生粗骨料的基础性研究对再生粗骨料进一步推广应用具有积极意义。基于此,本文对再生粗骨料进行物理整形,采集了再生粗骨料在整形过程中圆度、棱角度、扁平度、细长度、球度、粗糙度和分形维数等形态特征参数,并采用主成分分析法研究物理整形对再生粗骨料形态特征综合形态的影响并提出了评价指标,以期为再生粗骨料的生产制备及高品质应用提供基础性研究。

1 原材料及试验

1.1 原材料

试验选用原材料包括天然碎石粗骨料(G)、道路混凝土再生粗骨料(RCA1)和砖混再生粗骨料(B-RCA),其中RCA1和B-RCA的原生混凝土强度分别为C30和C20,B-RCA中砖含量约为35%。再生粗骨料均由拆除后的建筑垃圾经实验室小型颚式破碎机破碎筛分得到,粗骨料照片如图1所示。采用4.75~9.5 mm粒径粗骨料测试基本物理性能,如表1所示。

图1
试验所用粗骨料
pic
表1
粗骨料基本物理性能
粗骨料类型

表观密度/

(g∙cm-3)

堆积密度/(g∙cm-3)空隙率/%吸水率/%压碎指标/%
G2.741.6539.571.658.3
RCA12.611.3050.046.2017.8
B-RCA2.571.2451.928.9932.3
展开更多
1.2 试验方案
1.2.1 粗骨料物理整形

本文采用洛杉矶磨耗仪进行物理整形,洛杉矶磨耗仪为浙江路达机械仪器有限公司生产的DM-Ⅱ型,尺寸710 mm×510 mm,转速30~33 r/min,模拟反击式破碎机的动能冲击作用环境整形再生粗骨料[13]。每次试验称取固定质量4.75~19.5 mm的粗骨料倒入料筒中,并放入6个钢球强化整形效果,以5 min间隔逐次递增。当整形时间为15 min时,再生粗骨料出现严重破碎,最终将试验整形时间设置为0、5、10和15 min。由于再生粗骨料整形过程中粒径会发生变化,影响形态特征测试结果,因此将整形后的粗骨料筛分取样为4.75~9.5 mm粒径,并进行物理性能和形态特征测试。

1.2.2 粗骨料形态特征测试

试验粗骨料种类包括3种粗骨料(G、RCA1和B-RCA),每种粗骨料有4种整形时间(0、5、10和15 min),因此共有12种类型粗骨料。使用集料图像测量系统(AIMS)采集粗骨料的形态特征,每种粗骨料均随机选取100颗进行试验,试验过程图2所示。为从形状、棱角度和粗糙度3个尺度表征再生粗骨料,采用圆度、棱角度、扁平度、细长度、球度、粗糙度和分形维数计算粗骨料的形态特征,计算方法如表2所示。

图2
粗骨料形态特征指标采集过程
pic
表2
粗骨料粒形特征指标的计算方法
特征参数计算公式特征描述
圆度pic

圆度越大,粗骨料的形状越趋近圆形。

式中:A为粗骨料的投影面积;P为投影周长。

棱角度pic

棱角度越大,粗骨料的棱角性越明显。

式中:θ为图像边缘某点的梯度向量角度;n为边缘点的总数量;i为边缘的第i个点。

球度pic

球度越大,粗骨料形状越接近球形。

式中:DSDIDL分别为粗骨料的短轴、中轴和长轴长。

扁平度pic扁平度越小,粗骨料形状越趋于片状。
细长度pic细长度越小,粗骨料形状越趋于针状。
粗糙度pic

粗糙度越小,粗骨料表面越光滑。

式中:N为图像中细节系数的数量;ii张图像;j为小波系数;D为分解函数;x,y为转换域中的位置坐标。

分形维数pic分形维数越大,粗骨料的形状轮廓越复杂。式中:D为分形维数,k为尺度常数。
展开更多

2 试验结果

2.1 粗骨料物理性能

图3为不同整形时间时粗骨料的物理性能的变化。由图3可知,随着整形时间的增加,再生粗骨料的空隙率、吸水率和压碎指标均不断降低,当整形时间为15 min时,RCA1和B-RCA的空隙率较整形前分别降低了10.67%和10.24%、吸水率降低了23.23%和37.59%、压碎指标降低了25.84%和37.15%,表明物理整形能显著提高再生粗骨料的物理性能。这主要是由于在整形过程中,再生粗骨料表面的浆体在摩擦碰撞作用下不断剥落,附着浆体含量降低使得再生粗骨料物理性能得到提升。

图3
不同整形时间下粗骨料物理性能的变化
pic
2.2 粗骨料形态特征

表3为不同整形时间下粗骨料的粒形指标,由表3可知,随着整形时间的增加,再生粗骨料的圆度、扁平度、细长度和球度出现不同程度的增加,而棱角度、粗糙度和分形维数则出现降低。当整形时间为15 min时,RCA1、B-RCA的圆度分别增加了10.11%和7.03%;扁平度分别增加了7.43%和1.62%;细长度分别增加了14.39%和8.49%;球度分别增加了8.81%和7.38%;棱角度分别降低了24.82%和23.70%;粗糙度分别降低了13.59%和24.74%;分形维数分别降低了1.53%和1.06%。这表明再生粗骨料经过物理整形后,针片状含量减少、形状更趋于球形、表面更光滑、轮廓更规则,整形显著改善了再生粗骨料的形态特征。

表3
不同整形时间下粗骨料粒形指标的变化
粗骨料类型球磨时间/min圆度棱角度扁平率细长度球度粗糙度分形维数
G00.8302 726.20.6420.7250.688170.51.001
50.8432 404.30.6300.7310.683170.11.010
100.8432 371.60.6630.7110.688169.91.006
150.8512 305.00.6460.7310.687167.91.005
RCA100.7903 380.40.6460.6670.647212.61.046
50.8532 857.50.6680.7200.691198.91.015
100.8572 681.10.6820.7240.700207.11.010
150.8702 541.30.6940.7630.704183.71.030
B-RCA00.8073 478.70.6800.6950.682267.61.035
50.8442 792.50.6760.7250.698220.21.021
100.8552 596.70.6830.7320.705184.51.024
150.8642 654.30.6910.7540.732201.41.024
展开更多

3 主成分分析

由于每个粒形指标仅在单方面描述粗骨料的形态特征,当粒形指标种类较多时,复杂的高维信息不利于直观反映粗骨料颗粒形态。因此本文使用主成分分析方法将粒形指标重新组合,用少量几个主成分来反映多数粒形指标的信息,简洁高效地描述整形过程中再生粗骨料形态特征的变化规律。

3.1 数据预处理

主成分分析前需要通过KMO(Kaiser-Meyer-Olkin,KMO)和Bartlett球形检验方法(Bartlett’s Test of Sphericity)验证粒形指标间的相关性,使用SPSS 25.0软件(Statistical Product and Service Solutions 25.0)进行计算,得到粒形指标的KMO值为0.690(>0.5),Bartlett球形检验指标为71.361(p<0.05),这说明不同粒形指标间的相关性较好,可以使用主成分分析方法。由于不同粒形指标的计算公式不同,为了消除量纲对试验结果的影响,按式(1)采用标准差标准化方法计算表3中的原始粒形指标数据,标准化后的数据均值为0,标准差为1,计算结果如表4所示。

表4
数据标准化的粒形指标
粗骨料种类球磨时间/minZ1Z2Z3Z4Z5Z6Z7
G0-0.527-0.017-1.1690.073-0.207-0.888-1.320
50.032-0.895-1.7360.314-0.461-0.901-0.657
100.032-0.984-0.177-0.487-0.207-0.908-0.952
150.376-1.165-0.9800.314-0.258-0.977-1.025
RCA10-2.2481.766-0.980-2.250-2.2860.5661.995
50.4630.3410.059-0.127-0.0550.093-0.289
100.635-0.1400.7200.0330.4010.376-0.657
151.194-0.5211.2871.5960.604-0.4320.816
B-RCA0-1.5172.0340.626-1.129-0.5112.4661.185
50.0750.1640.4370.0730.3000.8290.153
100.549-0.3700.7680.3540.655-0.4040.374
150.936-0.2131.1461.2352.0240.1800.374
展开更多
表中:Z1Z2Z3Z4Z5Z6Z7分别为经过数据标准化后粗骨料的圆度、棱角度、扁平度、细长度、球度、粗糙度和分形维数值。
pic (1)

式中:pic为标准化后的粒形指标值,pic为原始粒形指标值,pic为第j个粒形指标均值,pic为第j个粒形指标的标准差,i=1,2,…,12,j=1,2,…,7。

3.2 粒形指标相关性分析

表4的标准化数据按式(2)计算所有粒形指标间的相关系数,计算结果如图4所示。

图4
不同粒形指标间的相关系数矩阵
pic
pic (2)

式中:rxy为2个粒形指标xy的相关系数;n为样本容量;xiyi分别为2个粒形的指标样本值;picpic分别为2个粒形指标的样本均值。

图4为不同粒形指标间的相关系数矩阵(R),由图4可知,正相关性较好的粒形指标分别有圆度和细长度、圆度和球度、细长度和球度、棱角度和粗糙度以及棱角度和分形维数。相关系数分别为0.87、0.82、0.80、0.85和0.78。负相关性较好的粒形指标有圆度和棱角度,相关系数为-0.74。这表明当粗骨料越趋于球形时,其针状和片状的含量会越少;当粗骨料表面的棱角越多时,其表面会越粗糙,轮廓会越不规则。

3.3 综合形态指标分析

将相关系数矩阵R式(3)计算特征值λi,并将结果列于表5中。

表5
各主成分特征值与贡献率
组分初始特征值提取荷载
特征值解释量/%累积量/%特征值解释量/%累积量/%
14.0557.8557.854.0557.8557.85
22.1931.2389.082.1931.2389.08
30.415.8294.90
40.182.5197.40
50.091.2898.68
60.060.8599.53
70.030.47100.00
展开更多
pic (3)

主成分贡献率反映出对原始数据信息的包含程度。按照特征值大于1的原则,提取出主成分1和主成分2,累计贡献率之和达到89.08%,说明能全面反映原始信息。为解释2个主成分的物理意义,采用最大方差法正交旋转初始载荷矩阵,计算出因子载荷矩阵,如表6所示。

表6
主成分旋转因子载荷矩阵
因子圆度棱角度扁平率细长度球度粗糙度分形维数
10.867-0.4290.7930.8600.950-0.060-0.106
2-0.4500.8700.545-0.405-0.1060.9260.852
展开更多

表6为主成分旋转因子载荷矩阵。由表6可知,因子1(F1)反映了再生骨料的圆度、扁平度、细长度和球度,载荷值分别为0.867、0.793、0.860和0.950,这4个粒形指标均侧重于描述粗骨料的整体形状特性,即判断粗骨料的形状是否接近球形,将因子1命名为整体形态因子。因子2(F2)中反映了棱角度、粗糙度和分形维数,载荷值分别为0.870、0.926和0.852。由于棱角度和粗糙度均可理解为粗骨料表面的形态规则程度,棱角度较大说明表面凸起较多,粗糙度较大说明表面凹凸不平的程度较大,将因子2命名为表面形态因子。

为了比较不同种类粗骨料的形态因子差异,使用SPSS对表6的因子载荷矩阵进行回归分析,计算得到因子得分系数矩阵,将其与标准化后的粒形指标值结合可得2个形态因子F1F2的综合得分表达式,如式(4)式(5)所示,F1越大说明粗骨料的整体形态越接近球形,F2越大说明粗骨料的表面形态越不规则。将表4中的数据代入式(4)式(5)中便可得到粗骨料2个形态因子的综合得分,结果如图5所示。

pic (4)pic (5)

式中:Z1Z2Z3Z4Z5Z6Z7分别为不同种类和整形时间下粗骨料的圆度、棱角度、扁平度、细长度、球度、粗糙度、分形维数值的数据标准化值。

图5
不同整形时间下粗骨料的形态因子综合得分
pic

图5(a)为不同种类粗骨料的整体形态因子得分。由图5(a)可知,在整形前粗骨料的整体形态因子得分排序为RCA1<G<B-RCA。粗骨料的整体形状主要受母岩强度和破碎方法影响[8]。对于颚式破碎机,当母岩强度越高时,生产出的粗骨料针片状含量也越高[13],因此原生混凝土强度的差异导致了破碎后再生粗骨料整体形态的不同,强度较高的RCA1整体形态因子得分最低。而G在实际生产过程中会经过冲击式破碎机或反击式破碎机的多次破碎加工,因此改善后的形状优于RCA1。随着整形时间的增加,粗骨料的整体形态因子均呈增长趋势,但再生粗骨料的增长幅度均大于天然粗骨料,当整形时间为5 min时,再生粗骨料的整体形态因子得分均超过了天然粗骨料。这表明相比于天然粗骨料,再生粗骨料的整体形态在整形过程中更容易改善,其原因是粗骨料的强度越低,针片状颗粒在整形过程中越容易破碎,因此再生粗骨料经过短暂整形后,整体形状比天然粗骨料更接近球形。

图5(b)为不同种类粗骨料的表面形态因子得分。由图5(b)可知,在整形前粗骨料的表面形态因子得分排序为G<RCA1<B-RCA,再生粗骨料的表面形态因子均高于天然粗骨料,这是由于再生粗骨料表面附着的浆体有较多棱角和孔隙,使其表面形态较天然粗骨料相比更不规则[3]。此外,浆体强度越低时孔隙含量会越高,强度较低的原生混凝土破碎后会生产出表面更不规则的再生粗骨料,因此B-RCA的表面形态因子高于RCA1。随着整形时间的增加,不同种类粗骨料的表面形态因子得分出现不同程度的降低,但再生粗骨料的表面形态因子得分始终大于天然粗骨料。在整形过程中,再生粗骨料附着浆体含量的减少改善了其表面形态特征,但浆体中存在的孔隙仍使再生粗骨料的表面形态比天然粗骨料更不规则,这表明物理整形方法有利于改善再生粗骨料的表面形态,但无法使其达到近似于天然粗骨料的程度。此外,由图5还可知,整形时间超过10 min时,RCA1和B-RCA的因子得分均出现了增长趋势,研究表明整形时间过长会使再生粗骨料发生二次破碎,产生了微裂纹增加了表面形态的不规则程度[4]。因此本次试验条件下的最佳整形时间为10 min,在工程实际中可结合表面形态因子的变化规律来合理控制再生粗骨料的整形时间。

相关研究表明,骨料的空隙率主要决于级配和形态特征[19]。为探究骨料形态因子和空隙率的关系,将G、RCA1和B-RCA的2个因子得分和空隙率线性拟合,图6为将G、RCA1和B-RCA的2个因子得分和空隙率进行拟合,得到的粗骨料整体形态因子和表面形态因子对空隙率的影响。由图6可知,整体形态因子、堆积密度的线性拟合表面形态因子和堆积密度的线性拟合的系数。由图6(a)可知,随着整体形态因子的增加,骨料空隙率降低,这表明经过物理整形后,再生粗骨料整体形态因子增加,堆积密实程度也因此提高。由图6(b)可知,表面形态因子和堆积密度的线性相关性,表面形态因子对空隙率的影响规律不明显。

图6
空隙率和因子得分的相关性分析
pic

4 结论

1) 物理整形能有效改善再生粗骨料的物理性能和形态特征。本文实验条件下,整形时间为15 min时,RCA1和B-RCA的空隙率较整形前降低了10.67%和10.24%,吸水率降低了23.23%和37.59%,压碎指标降低了25.84%和37.15%,且再生粗骨料经过物理整形后,针片状含量减少,形状更趋于球形,表面更光滑,轮廓更规则,整形显著改善了再生粗骨料的形态特征。

2) 采用主成分析的方法能够将再生粗骨料的多个形态特征参数进行归一化简化。圆度、球度、细长度和扁平度可以组合成整体形态因子,从而更好地表征粗骨料的整体形状特征,而棱角度、粗糙度和分形维数则可以组合成了表面形态因子,则能够更好地表征粗骨料表面的规则程度。整体形态因子增大,则粗骨料的形状越趋于球形,粗骨料的堆积密实程度会越好,空隙率降低;表面形态因子增大,则粗骨料表面形态越不规则。

3) 随着物理整形时间的增加,再生粗骨料可能发生二次破碎,将导致再生粗骨料的表面形态因子增加,从而增大粗骨料表面的不规则程度,不利于再生粗骨料的使用,因此结合本文研究建议根据表面形态特征的变化合理控制整形时间。

参考文献
1LIU Zhihui, ZHAO Yangang, MA Lie, et al.

Review on high-strength recycled aggregate concrete: mix design, properties, models and structural behaviour

[J]. Structures, 2024, 64: 106598.
百度学术谷歌学术
2马昆林, 莫文波, 莫建红, .

砖混再生细粗骨料再生砂浆性能及工程应用

[J].铁道科学与工程学报, 2021, 18(8): 2073-2080.
百度学术谷歌学术
3MA Kunlin, HUANG Xinyu, SHEN Jingtao, et al.

The morphological characteristics of brick-concrete recycled coarse aggregate based on the digital image processing technique

[J]. Journal of Building Engineering, 2021, 44: 103292.
百度学术谷歌学术
4付鲁鑫, 贾致荣, 王立志, .

研磨功对再生集料性能改善的量化评价

[J]. 科学技术与工程, 2020, 20(31): 12986-12991.
百度学术谷歌学术
5PAUZI N N M, JAMIL M, HAMID R, et al.

Influence of spherical and crushed waste Cathode-Ray Tube (CRT) glass on lead (Pb) leaching and mechanical properties of concrete

[J]. Journal of Building Engineering, 2019, 21: 421-428.
百度学术谷歌学术
6NADERI M, KABOUDAN A.

Experimental study of the effect of aggregate type on concrete strength and permeability

[J]. Journal of Building Engineering, 2021, 37: 101928.
百度学术谷歌学术
7马昆林, 刘建, 申景涛, .

骨料强化方法对再生混凝土多界面过渡区微观结构的影响

[J]. 铁道科学与工程学报, 2023, 20(10): 3809-3819.
百度学术谷歌学术
8张雪芹, 马昆林, 龙广成, .

粗骨料形态特征表征参数及其与混凝土性能关系的研究进展

[J]. 材料导报, 2024, 38(2): 118-129.
百度学术谷歌学术
9ZHANG Saisai, LI Rui, PEI Jianzhong.

Evaluation methods and indexes of morphological characteristics of coarse aggregates for road materials: a comprehensive review

[J]. Journal of Traffic and Transportation Engineering (English Edition), 2019, 6(3): 256-272.
百度学术谷歌学术
10LI Ruize, LU Wenbo, CHEN Ming, et al.

Quantitative analysis of shapes and specific surface area of blasted fragments using image analysis and three-dimensional laser scanning

[J]. International Journal of Rock Mechanics and Mining Sciences, 2021, 141: 104710.
百度学术谷歌学术
11李良, 龙广成, 谢友均, .

基于三维扫描与数值技术的粗集料形状特征与级配研究

[J]. 铁道科学与工程学报, 2022, 19(3): 714-721.
百度学术谷歌学术
12聂志红, 廖靖云, 周苏华, .

火山渣颗粒图像处理及形状定量分析

[J]. 湖南大学学报(自然科学版), 2020, 47(3): 131-139.
百度学术谷歌学术
13KAMANI M, AJALLOEIAN R.

The effect of rock crusher and rock type on the aggregate shape

[J]. Construction and Building Materials, 2020, 230: 117016.
百度学术谷歌学术
14肖柏林, 杨志强, 高谦, .

混合骨料堆积模型研究与新型模型建立

[J]. 材料导报, 2018, 32(14): 2400-2406.
百度学术谷歌学术
15WU Junqing, ZHOU Xiang, ZENG Xiaohui, et al.

Effect of aggregate morphology characteristics on the voidage of aggregate loose packing based on 3D discrete element method

[J]. Construction and Building Materials, 2022, 348: 128598.
百度学术谷歌学术
16WANG Yuting, DIAMOND S.

A fractal study of the fracture surfaces of cement pastes and mortars using a stereoscopic SEM method

[J]. Cement and Concrete Research, 2001, 31(10): 1385-1392.
百度学术谷歌学术
17YE Sheng, HAN Ding, LI Yang.

Generation and evaluation of aggregate section using the fractal dimension

[J]. Case Studies in Construction Materials, 2019, 11: e00255.
百度学术谷歌学术
18朱建平. 应用多元统计分析[M]. 北京: 北京大学出版社, 2017: 265-280.
19骆行文, 姚海林.

基于主成分分析的岩石质量综合评价模型与应用

[J]. 岩土力学, 2010, 31(S2): 452-455.
百度学术谷歌学术
20肖柏林, 杨志强, 陈得信, .

充填骨料颗粒形态参数定量评价方法

[J]. 天津大学学报(自然科学与工程技术版), 2019, 52(5): 545-553.
百度学术谷歌学术
注释

申景涛,黄新宇,马昆林等.基于主成分分析的物理整形对再生粗骨料形貌特征影响的研究[J].铁道科学与工程学报,2025,22(03):1112-1120.

SHEN Jingtao,HUANG xinyu,MA Kunlin,et al.Influence of physical shaping on morphological characteristics of recycled coarse aggregates based on principal component analysis[J].Journal of Railway Science and Engineering,2025,22(03):1112-1120.