扩展功能
文章信息
- 张蕾, 惠新, 陈龙, 王旭东
- ZHANG Lei, HUI Xin, CHEN Long, WANG Xu-dong
- 基于图像技术的水泥稳定级配碎石加载损伤特征研究
- Study on Damage Characters of Graded Cement Stabilized Crushed Stone under Loading Based on Image Processing Technique
- 公路交通科技, 2017, 34(6): 38-44
- Journal of Highway and Transportation Research and Denelopment, 2017, 34(6): 38-44
- 10.3969/j.issn.1002-0268.2017.06.006
-
文章历史
- 收稿日期: 2016-07-15
2. 南京铁路建设投资有限责任公司, 江苏 南京 210022;
3. 重庆交通大学, 重庆 400074
2. Nanjing Railway Building & Investment Co., Ltd., Nanjing Jiangsu 210022, China;
3. Chongqing Jiaotong University, Chongqing 400074, China
水泥稳定级配碎石是我国公路建设中普遍采用的基层材料之一,与沥青混合料类似,其体积组成是影响性能的重要参数,然而,由于水泥稳定级配碎石拌和、成型过程中水泥与水发生水化反应,无法像沥青混合料一样测定或计算其理论密度,因此,也就无法通过计算获得水泥稳定级配碎石材料的孔隙率。同时,开裂是水泥稳定级配碎石材料的典型破坏,但以往的研究手段无法获得裂缝萌生、扩展直至破坏的过程,以及材料设计对损坏劣化模式的影响。
工业CT技术分析材料组成结构的基本原理是通过射线在具有不同密度的物质之间的能量衰减差异,从而反算获得混合料内部不同物质以及物质的形态、体积等信息[1-8]。基于这一原理,采用工业CT扫描水泥稳定级配碎石材料,可以将水泥稳定级配碎石材料的不同组成部分石料(2.7~2.9 kg/cm3)、水泥砂浆(2.0 kg/cm3)及空气和裂缝(1.0~1.3 kg/cm3)识别出来。
文中采用交通运输部公路科学研究院的工业CT设备及图像重构软件,对5种不同级配组成的水泥稳定级配碎石混合料进行扫描和分析,获取试件的孔隙率及分布特征,探讨了孔隙分布特征与级配的相关性;对其中两种级配的混合料进行了单轴压缩模式的逐级加载试验和同步CT扫描,通过断层图像中裂纹的识别与分析,探讨裂纹扩展规律。
1 试验部分 1.1 原材料和级配采用内蒙古某工程的石灰岩集料,各粒径集料毛体积密度、压碎值、针片状颗粒含量等均采用《公路工程集料试验规程》(JTG E42—2005)[9]的方法进行检测,结果均符合规范要求。水泥为325# 普通硅酸盐水泥。
选择筛孔4.75 mm通过率分别为30%,35%,40%,45%和50%的5种级配,级配构成方法采用幂函数间断级配[10],编号为C1-C5。按照《公路工程无机结合料稳定材料试验规程》(JTG E51—2009)[11]的方法,根据水泥稳定级配碎石击实试验和抗压强度试验结果,确定最佳含水量和水泥剂量[12]。结果见表 1。
编号 | 通过各种筛孔(mm)的通过率/% |
水泥剂
量/% |
最佳含
水率/% | ||||||||||
26.5 | 19 | 13.2 | 9.5 | 4.75 | 2.36 | 1.18 | 0.6 | 0.3 | 0.15 | 0.075 | |||
C1 | 100 | 81 | 63.8 | 50.9 | 30 | 20.7 | 14.4 | 10.3 | 7.3 | 5.4 | 4 | 7 | 5.15 |
C2 | 100 | 82.4 | 66.4 | 54.4 | 35 | 23.9 | 16.5 | 11.5 | 8 | 5.6 | 4 | 6 | 5.25 |
C3 | 100 | 83.7 | 69 | 57.9 | 40 | 27.1 | 18.5 | 12.7 | 8.6 | 5.9 | 4 | 8 | 5.24 |
C4 | 100 | 85.1 | 71.6 | 61.4 | 45 | 30.3 | 20.5 | 13.9 | 9.3 | 6.1 | 4 | 8 | 5.20 |
C5 | 100 | 86.5 | 74.1 | 64.9 | 50 | 33.6 | 22.5 | 15.1 | 9.9 | 6.4 | 4 | 7 | 5.23 |
1.2 试验方法
(1)试件成型
按照上文确定的设计参数,采用静压成型的方法成型试件,每一级配6个φ150 mm×h150 mm的圆柱型平行试件,标准条件养生7 d后,强度满足要求,钻取φ100×h100 mm的芯样。
(2)图像采集及重构
采用交通运输部公路科学研究院的工业CT扫描设备进行试件的断层图像扫描,该设备为美国某公司生产,最大电压为320 kV,理论最大密度分辨率为0.2%,空间分辨率约为0.2 mm,扫描方式为锥束射线。通过多次试验,确定了用于此次扫描的射线参数为电压270 kV,电流2.9 A,采用小焦点采样方式。根据720 次(0.5°/rad)扫描投影获得的数据进行投影反算和图像重构。
(3)孔隙及裂纹识别
运用Matlab语言编制图像处理程序[13],对所得到的CT扫描图片进行图像预处理和自动识别分析,获得孔隙和损伤后裂纹图像,在分辨孔隙和裂纹的过程中,主要基于裂纹具有长度与宽度比大的特征,当长宽比大于设定阈值时,即判定该对象为裂纹,否则判定为孔隙, 如图 1所示。
![]() |
图 1 孔隙和损伤后裂纹二值图 Fig. 1 Binary image of air voids and cracking after damage |
|
(4)加载与损伤检测
以往的研究,受限于CT试验机的结构和工作原理,多采用试件加载与损伤扫描分离的方式[14-16],即首先在加载设备上进行试件的加载,再将加载后损伤的试件移至CT试验仪内进行扫描,这样的操作模式弊端明显,对于准脆性材料,卸载导致试件的裂缝闭合,从而无法获得真实的试件损伤过程,为了克服这一弊端,自行研发的工业CT扫描机同步加载试验装置[17],在工业CT扫描设备内,通过自反力加载的方式,实现同步逐级加载和损伤扫描,设备旋转加载油缸及可透视碳纤维反力架, 如图 2所示。
![]() |
图 2 工业CT同步加载装置 Fig. 2 Industrial CT synchronous loading equipment 1—单作用旋转加载油缸;2—可透视碳纤维反力架;3—油缸连接筒;4—球形压头;5—试验试件;6—尼龙垫块;7—底板连接法兰;8—插销 |
|
2 试验结果与分析 2.1 孔隙分布特征
(1)孔隙率沿试件纵向分布规律
首先,计算每一层断面上的孔隙像素数量,分析5种水泥稳定级配碎石芯样孔隙率沿试件深度方向的变化规律。5个位置分别为扫描第100~150,300~350,500~550,700~750,850~900个断面,每一个位置统计50 个断面的累计值。结果表明:对于静压成型的水泥稳定级配碎石芯样,试件内部孔隙率并非均匀分布,而是从上至下逐渐减小,即下部密度将明显大于上部,如图 3所示。
![]() |
图 3 芯样不同深度处孔隙特征值统计结果 Fig. 3 Statistical result of air voids characteristic values of coring sample in different depths |
|
(2)级配对孔隙率的影响
统计了孔隙面积率随筛孔4.75 mm通过率的变化规律,如图 4所示,结果可见,随着筛孔4.75 mm通过率逐渐增加,混合料的孔隙特征值将首先减小,达到最小值,随后增加,即孔隙特征值并不是随着细集料数量的增加持续减小,而是到达一定填充密实程度后,孔隙反而会增加,图 5为筛孔4.75 mm混合料不同通过率的典型断层图片。进一步分析不同筛孔4.75 mm 通过率下孔隙的尺度和数量特征的变化趋势。
![]() |
图 4 不同级配孔隙特征值统计结果 Fig. 4 Statistical result of air voids characteristic values of different gradations |
|
![]() |
图 5 混合料不同通过率(筛孔4.75 mm)典型断层图片 Fig. 5 Typical slice images of 4.75 mm mixture with different passing rates |
|
(3)级配对孔隙数量的影响
统计每种级配混合料的所有切片上孔隙个数,可见,随着筛孔4.75 mm通过率的增加,孔隙个数先减小,达到最小值,而后随着细集料含量的增加,呈增加趋势。即当细集料含量(30%)低时,随细集料填充量的增加,粗集料形成的主骨架孔隙残留量将逐步减少,而后,随着细集料数量增加,孔隙被逐级填充,试件内部可见孔隙数量逐步减少。但当细集料含量进一步增加,将逐渐撑开粗集料所构成的骨架,细集料的增加使压实过程中的阻力增加,细集料之间的摩阻力增大,细集料自身的孔隙逐渐占据主导,且数量随细集料含量增加逐步增加,如图 6所示。
![]() |
图 6 孔隙个数与通过率(筛孔4.75 mm)相关关系 Fig. 6 Relationship between number of air voids and 4.75 mm passing rate |
|
(4)级配分形特征与孔隙尺度分布频率
由上述(2)和(3)可知,不同级配的混合料其空隙在总量和数量上均有显著差异,那么,说明孔隙在尺度分布上具有明显差异。因此,下面采用统计学方法对所有断面的孔隙进行分布特征的统计。
首先,尝试将分形维度的概念引入级配的分析中,并将分形维度与孔隙分布频度建立联系。5 种水泥稳定级配碎石的分形维度计算结果见表 2,结果表明:分形维度随筛孔4.75 mm通过率增加而增加。另一方面,统计了不同尺度孔隙(单个孔隙面积)所占的比例,见表 3。
筛孔4.75 mm
通过率/% | 回归方程 |
分形维度
Di |
相关系数
R2 |
30 | y=0.554 6x+2.678 9 | 2.445 4 | 0.993 6 |
35 | y=0.549x+2.750 2 | 2.451 | 0.998 7 |
40 | y=0.543 6x+2.814 6 | 2.456 4 | 0.999 3 |
45 | y=0.539 1x+2.872 5 | 2.460 9 | 0.996 6 |
50 | y=0.534 4x+2.926 8 | 2.465 6 | 0.991 6 |
筛孔 4.75 mm 通过率/% | 各尺度孔隙(pixel)占所有孔隙的比例/% | |||||||
0~ 30 |
30~ 60 |
60~ 90 |
90~ 120 |
120~ 150 |
150~ 180 |
180~ 210 | >210 | |
30 | 55.9 | 16.6 | 8.1 | 4.7 | 3.3 | 2.2 | 1.5 | 7.8 |
35 | 58.9 | 16.8 | 7.9 | 4.4 | 2.8 | 2.3 | 1.5 | 5.5 |
40 | 66.1 | 16.4 | 6.6 | 3.5 | 2 | 1.5 | 0.9 | 3 |
45 | 69.6 | 15.2 | 6.3 | 3 | 1.9 | 1.1 | 0.7 | 2.1 |
50 | 73.3 | 14.7 | 5.5 | 2.8 | 1.3 | 0.7 | 0.4 | 1.3 |
计算各尺度孔隙所占比例与级配分形维度的相关性,结果见表 4。
孔隙尺度/ pixel |
0~ 30 |
30~ 60 |
60~ 90 |
90~ 120 |
120~ 150 |
150~ 180 |
180~ 210 | >210 |
相关系数 | 0.992 | -0.887 | -0.981 | -0.984 | -0.990 | -0.947 | -0.953 | -0.989 |
相关性分析结果可见,对于不同尺度范围的孔隙,通过5种级配的分析可见,级配的分形维度与不同尺度的孔隙频度均有较好的相关性,随着筛孔4.75 mm通过率的增加,孔隙尺度范围在0~30的小尺度孔隙数量逐渐增加,其他尺度范围的孔隙随筛孔4.75 mm通过率的增加,孔隙数量所占的比重逐渐减小。
2.2 裂纹随荷载作用的扩展规律采用位移控制模式(1 mm/min加载速率)进行加载,直至设定的荷载级位后,保持变形对试件进行CT扫描后卸载,分别进行不同级位的加-卸载试验[18]。由于裂纹与孔隙密度小,因此在CT图像中表现色彩偏重,利用这一特点,通过图像处理的方法,将孔隙与裂纹从集料与砂浆中分离出来。
分析中,将试件承受20 kN荷载的图像作为损伤初始值,假定逐级加载对试件将产生不同程度的损伤,表现为微裂纹的出现、孔隙扩大以及裂纹的扩展,分别采用相同的图像分析方法,对各级荷载后试件进行扫描和分析,并与20 kN荷载的图像初始值比较,获得每一级荷载作用后,相较与初始试件的累计损伤。
分别选择筛孔4.75 mm通过率为30%和45%的C1和C4级配进行逐级加载和损伤扫描试验,课题组目前尚未进行裂纹形态的定量分析,仅就每一级荷载作用下,较之初始试件的累计损伤像素总量进行了分析,结果如图 7所示。
![]() |
图 7 损伤增量随荷载变化规律 Fig. 7 Damage increment varying with loading |
|
对于悬浮密实结构的C4混合料,当试件达到破坏临界荷载前,试件内部损伤很小,主要表现在部分集料内部存在裂纹的延展,数量微少,如图 8所示。但当材料达到破坏临界荷载后,试件内部迅速出现裂缝,多数试件的裂缝呈环状扩展分布,并且裂缝主要沿集料与砂浆的边缘发展,并相互贯通。集料的粒径与裂纹的扩展有密切关系,裂缝更容易出现在较大、较为分散的颗粒周边,并逐渐扩展。与悬浮结构的C4级配显著不同,C1级配内部具有较多孔隙,如图 9所示,这一特点给裂纹尤其是微裂纹的识别带来了较大难度。直接将最末一级荷载的图像与初始图像相比较,可以看到试件内部分布的裂缝,但裂缝的发展特点与N9级配显著不同,裂缝较为均匀地分布在整个断面,而不似N9级配以环状形态出现。
![]() |
图 8 C4混合料不同级位裂纹图 Fig. 8 Cracking images of C4 mixture under different loading levels |
|
![]() |
图 9 C1混合料不同级位裂纹图 Fig. 9 Cracking images of C1 mixture under different loading levels |
|
2.3 圆柱形试件沿纵向裂纹分布状态定性分析
在单轴压缩试验模式下,对于C4级配,圆柱形试件内部的裂纹并非均匀分布,由于试件两端存在的端部约束作用,在试件的纵深方向,可见在试件两端,裂缝宽度较小,而且主要集中于试件的最外环,而从两端逐渐向试件中部变化过程中,裂缝宽度变大,并向试件中心部分深入,如图 10所示。对于C1级配,规律并不明显,裂缝基本处于零散、随机分布,如图 11所示,但由于C1级配内部存在大量孔隙,图像分析中明显可见,由于孔隙的大量存在而导致的裂缝误判,因此,对于存在大量孔隙的级配类型,裂缝分析的图像处理方法应进一步改进,将引入模式识别,通过判断目标的形状而筛选孔隙与裂缝。
![]() |
图 10 C4混合料裂缝沿试件纵向分布 Fig. 10 Cracking distribution of C4 mixture along vertical direction |
|
![]() |
图 11 C1混合料裂缝沿试件纵向分布 Fig. 11 Cracking distribution of C1 mixture along vertical direction |
|
3 结论
文中采用x-ray 计算机断层扫描技术,分析了水泥稳定级配碎石的孔隙分布特征尤其是孔隙尺度分布特征和加载过程中的裂纹扩展特征。研究结果表明:
(1)在静压成型条件下,水泥稳定级配碎石的孔隙在试件纵向方向并非均匀分布。当粗集料含量较高时,混合料的孔隙数量较少,且尺度偏大。
(2)引入分形维度的概念分析混合料级配,计算各级配的分形维数。结果表明随着细集料的增加,分形维数逐渐增大,分形维度可以作为级配的特征参数。分形维数与不同尺度的孔隙所占比例具有良好的相关性,由此建立了级配分形维数与孔隙尺度分布频率之间的相关关系。
(3)采用课题组自行研制的工业CT同步加载装置,开展的加载过程裂缝发展规律研究,结果表明,级配对于裂纹的扩展趋势同样有显著影响。
进一步,将探讨通过裂纹定量分析对混合料的损伤变量及损伤演化方程进行研究。
[1] | 张朝宗. 工业CT技术和原理[M]. 北京: 科学出版社, 2009. ZHANG Chao-zong. Technology and Principle of Industry Computer Tomography[M]. Beijing: Science Press, 2009. |
[2] | 易富, 金艳, 苏剑. 两种压实方法下沥青混合料空隙分布特性细观分析[J]. 公路交通科技, 2014, 31(3): 35-42 YI Fu, JIN Yan, SU Jian. Mesoscopic Analysis of Air Void Distribution of Asphalt Mixture in Two Compaction Methods[J]. Journal of Highway and Transportation Research and Development, 2014, 31(3): 35-42 |
[3] | 李智, 陈思宇. 采用X射线CT技术评价沥青路面结构内部质量均匀性[J]. 公路交通科技, 2014, 31(12): 78-83 LI Zhi, CHEN Si-yu. Evaluating Homogeneity of Interval Quality of Asphalt Pavement Structure with X-ray Computed Tomography[J]. Journal of Highway and Transportation Research and Development, 2014, 31(12): 78-83 |
[4] | 汪海年, 卜胤, 周洋, 等. 基于X-Ray CT技术的沥青混合料扫描影响因素及其优化研究[J]. 公路交通科技, 2014, 31(11): 9-15 WANG Hai-nian, BU Yin, ZHOU Yang, et al. Research on Influencing Factors and Its Optimization of Scanning Asphalt Mixture Based on X-ray CT Technology[J]. Journal of Highway and Transportation Research and Development, 2014, 31(11): 9-15 |
[5] | 王宏, 刘锋. 不同成型方式泡沫沥青冷再生混合料细微观结构性能研究[J]. 公路交通科技, 2016, 33(2): 19-27 WANG Hong, LIU Feng. Research of Microstructure Performance of Foamed Asphalt Cold Recycled Mixture in Different Compaction Methods[J]. Journal of Highway and Transportation Research and Developmen, 2016, 33(2): 19-27 |
[6] | MASAD E, JANDHYALA V K, DASGUPTA N, et al. Characterization of Air Void Distribution in Asphalt Mixes using X-ray Computed Tomography[J]. Journal of Materials in Civil Engineering, 2002, 14(2): 122-129 |
[7] | WANG L, PAUL H S, HARMAN T, et al. Characterization of Aggregates and Asphalt Concrete Using X-ray X-ray Computerized Tomography[J]. Journal of the Association of Asphalt Paving Technologists, 2004, 73: 467-500 |
[8] | MASAD E, MUHUNTHAN B, SHASHIDHAR N, et al. Internal Structure Characterization of Asphalt Concrete Using Image Analysis[J]. Journal of Computing in Civil Engineering, 1999, 13(2): 88-95 |
[9] | JTG E42—2005,公路工程集料试验规程[S]. JTG E42—2005,Test Methods of Aggregate for Highway Engineering[S]. |
[10] | 王旭东, 张蕾. 基于骨架嵌挤型原理的沥青混合料均衡设计方法[M]. 北京: 人民交通出版社, 2014. WANG Xu-dong, ZHANG Lei. Multi-performance Design of Hot Mix Asphalt Based on Stone Interlock Structure Theory[M]. Beijing: China Communications Press, 2014. |
[11] | JTG E51—2009,公路工程无机结合料稳定材料试验规程[S]. JTG E51—2009, Test Methods of Materials Stabilized with Inorganic Binders for Highway Engineering [S]. |
[12] | 寇冠宇.半刚性基层材料体积参数与疲劳性能研究[D].重庆:重庆交通大学,2015. KOU Guan-yu. Study on Volume Parameter and Fatigue Performance of Semi-rigid Base Material[D].Chongqing: Chongqing Jiaotong University,2015. |
[13] | 冈萨雷斯RC. 数字图像处理(MATLAB版)[M]. 北京: 电子工业出版社, 2005. GONZALEZ R C. Digital Images Progress (MATLAB Edition)[M]. Beijing: Electric Industry Press, 2005. |
[14] | HAO Ying. Using X-ray Computed Tomography to Quantify Damage of Hot-mix Asphalt in the Dynamic Complex Modulus and Flow Number Test[D]. Baton Rouge: Louisiana State University,2010. |
[15] | 李晓军, 张肖宁, 武建民. 沥青混合料单轴重复加卸载破损CT识别[J]. 哈尔滨工业大学学报, 2005, 37(9): 1228-1231 LI Xiao-jun, ZHANG Xiao-ning, WU Jian-min. Internal Structure Discrimination of Asphalt Mix with X-ray Computerized Tomograph under Repetitive Uniaxial Loadings[J]. Journal of Harbin Institute of Technology, 2005, 37(9): 1228-1231 |
[16] | 谢涛.基于CT实时观测的沥青混合料裂纹扩展行为研究[D].成都:西南交通大学,2006. XIE Tao. Study on Crack Propagation Behavior of Asphalt Mixture Based on CT Real-time Observation [D]. Chengdu: Southwest Jiaotong University,2006. |
[17] | 张蕾,王旭东,周兴业,等.工业CT扫描机同步加载试验装置及工业CT扫描机: 中国,201520647067.4 [P].2015-08-25. ZHANG Lei, WANG Xu-dong, ZHOU Xing-ye,et al. Industry Computer Tomography Scanner Synchronized Loading Equipment and Industry Computer Tomography Scanner: China, 201520647067.4[P].2015-08-25. |
[18] | 田威, 党发宁, 陈厚群. 单轴压缩条件下混凝土细观损伤演化机理的CT试验研究[J]. 实验力学, 2011, 26(2): 54-60 TIAN Wei, DANG Fa-ning, CHEN Hou-qun. Experimental Study of Concrete Meso-Damage Evolution Mechanism Based on Uniaxial Compression CT Method[J]. Journal of Experimental Mechanics, 2011, 26(2): 54-60 |