混凝土中钢筋时常发生锈蚀,不仅会减小钢筋的有效承载面积,而且还会对钢筋的屈服强度、极限强度、屈服应变、极限应变等力学和变形性能产生显著影响。因此,准确、全面地分析锈蚀钢筋的力学和变形性能,对于在役钢筋混凝土结构的耐久性分析和安全性评估具有重要意义。
围绕锈蚀钢筋的力学和变形性能已开展了大量试验研究。其中,文献[1-6]分析了质量锈蚀率与锈蚀钢筋的屈服强度、极限强度、极限应变、屈服平台长度等力学和变形性能参数之间的关系;Kashani等[7]分析了质量锈蚀率对受拉钢筋延性的影响;张伟平等[8-9]分析了平均截面锈蚀率与锈蚀钢筋的力学与变形性能之间的关系。由于上述质量锈蚀率和平均截面锈蚀率只能从平均意义上描述钢筋的锈蚀程度,所以,只适用于均匀锈蚀的情况。然而,混凝土属于典型的非均质材料,其内部随机分布的骨料、孔隙和裂缝等导致混凝土中钢筋的空间锈蚀分布形态具有显著的非均匀性[10-11]。因此,利用最大截面锈蚀率、最小剩余截面积和点蚀因子等特征参数能够更加合理地描述钢筋的非均匀锈蚀程度。对此,罗小勇等[12]、Lu等[13]分析了最大截面锈蚀率与锈蚀钢筋的力学和变形性能之间的关系;曲福来等[14]分析了各微段的最大质量锈蚀率与锈蚀钢筋的屈服强度、极限强度和伸长率之间的关系;孙晓燕等[15]建立了最小剩余截面积相对值与锈蚀钢筋的屈服强度相对值、极限强度相对值和伸长率之间的关系;Zhu等[16]分析了最大截面锈蚀率对锈蚀钢筋延性的影响;Sheng等[17]研究了锈坑形状、深度和分布对钢筋的极限强度、屈服强度和伸长率的影响;Tang等[18]分析了最小剩余截面积与锈蚀钢筋的力学和变形性能之间的关系;Lu等[19]、Imperatore等[20]分析了面积点蚀因子对锈蚀钢筋的屈服强度、极限强度和伸长率的影响;Balestra等[21]分析了最小直径对锈蚀钢筋的屈服强度、极限强度和伸长率的影响。由此可见,目前的研究成果大多局限于单一或有限的锈蚀特征参数分析,难以准确、全面地描述钢筋的均匀/非均匀锈蚀形态对力学和变形性能的影响。因此,有必要进一步开展考虑均匀/非均匀锈蚀特征参数的钢筋力学和变形性能分析研究。
笔者首先开展了45个混凝土试件中135根加速锈蚀钢筋的单轴拉伸试验,然后,基于试验数据分别建立了均匀锈蚀特征参数(质量锈蚀率)和非均匀锈蚀特征参数(包括最大截面锈蚀率、最大锈蚀深度相对值、最小剩余截面积相对值和面积点蚀因子)与锈蚀钢筋的力学性能参数(包括名义屈服强度相对值和名义极限强度相对值)和变形性能参数(包括屈服应变相对值、极限应变相对值和屈服平台长度相对值)之间的定量关系,进而分别建立了基于上述锈蚀特征参数的锈蚀钢筋应力-应变关系模型,并与试验数据进行了验证分析。
本文试验与余波等[22]试验基于同一批试件开展,但二者的试验内容和研究侧重点存在明显区别。其中,文献[22]重点关注混凝土中钢筋的空间锈蚀形态测试装置和空间锈蚀特征参数分析方法,而本文重点开展锈蚀钢筋的单轴拉伸试验。
如图 1所示,共制作45个长、宽和高分别为500、d+30和150 mm的混凝土试件[22]。其中,d为钢筋直径。混凝土中共埋设135根钢筋,包括直径为8和12 mm的HPB300级光圆钢筋各18根,直径为12 mm的HRB400E级带肋钢筋18根,以及直径为16 mm的HRB400E级带肋钢筋81根。每个试件沿高度方向等间距放置3根钢筋。混凝土强度等级为C30,混凝土配合比和钢筋的性能参数信息见文献[22]。待锈蚀钢筋的长度为450 mm,钢筋在加工完后经过酸洗、烘干、称重获取钢筋锈蚀前的质量M0。在混凝土试件的两端表面涂3层环氧树脂,以防止NaCl溶液从混凝土端部渗入混凝土内部,待环氧树脂凝固后,将试件放入温度为20±2 ℃、湿度为95%的标准养护室中养护7 d。试验采用外加直流电的方法加速混凝土中的钢筋锈蚀,通过控制通电时间可以得到不同锈蚀程度的钢筋[22]。通电结束后,对混凝土试件进行破型,取出锈蚀钢筋,经过酸洗、干燥、称重等步骤,获取钢筋锈蚀后的质量M1。
利用自主研发的连续获取钢筋全角度锈蚀形态图像的测试装置[23]来获取钢筋的非均匀锈蚀形态,该装置主要由底座、钢筋固定装置、钢筋驱动装置、调节支架和数码相机等组成。利用该装置可以获得钢筋的非均匀锈蚀形态图像(识别精度为0.1 mm),并获得钢筋全角度锈蚀轮廓的上、下边缘线到轴线的距离rik和rn+1+ik及其对应的角度αi和αn+1+i(i=0, 1, 2, ..., n,为旋转角度的编号;k=1, 2, 3, ..., m,为沿钢筋纵向的控制截面编号;n和m分别为旋转总次数和控制截面总数),以及锈蚀钢筋控制截面的剩余截面积Sk,进而可以计算钢筋的均匀锈蚀特征参数(质量锈蚀率ηm, a)和非均匀锈蚀特征参数(包括最大截面锈蚀率ηs, max、最小剩余截面积Smin、最大锈蚀深度Dmax和面积点蚀因子Ps)[22]。
式中:S0和r0分别为未锈蚀钢筋的横截面积和半径;Smin、S和rmin分别为锈蚀钢筋的最小剩余截面积、平均剩余截面积和最小剩余半径。
根据金属材料拉伸试验规范[24]的操作规定,采用WAW-600型微机控制电液伺服万能试验机开展锈蚀钢筋的轴向拉伸试验。用钢筋打点机每隔10 mm对待测钢筋进行打点处理,并用油性笔在钢筋锈蚀最严重部位的5倍公称直径范围两端做标记。通过轴向拉伸试验,可以获得不同锈蚀程度钢筋的荷载-位移曲线,如图 2所示。编号P和R分别表示HPB300级钢筋和HRB400E级钢筋,后面的数字代表钢筋直径,Mi(i=1, 2, 3, …9)表示9种混凝土试件的配合比,A、B、C分别表示在每种配合比条件下的3组混凝土试件,-Ai(i=1, 2, 3)表示每个混凝土试件中不同位置钢筋的编号,-i(i=1, 2, 3, …5)表示未锈蚀钢筋的编号。由图 2可知,随着质量锈蚀率的增加,钢筋的屈服荷载和极限荷载及其对应的位移都逐渐降低,屈服平台也随着质量锈蚀率的增加而逐渐减小,甚至消失。此外,质量锈蚀率对锈蚀钢筋荷载-位移曲线在弹性阶段的斜率影响相对较小。
根据图 2中不同锈蚀程度钢筋的荷载-位移曲线,可以确定钢筋锈蚀前后的应力-应变曲线,如图 3所示。基于Park-Paulay模型[25],可以将锈蚀钢筋的应力-应变曲线划分为弹性段、屈服平台段及强化段,对应的应力-应变曲线为
式中:$m=\frac{\left(f_{\mathrm{uc}} / f_{\mathrm{yc}}\right)(30 r+1)^{2}-60 r-1}{15 r^{2}}$;r=εuc-εhc;ζ=εs-εhc;Es为锈蚀钢筋的弹性模量;σs和εs分别为锈蚀钢筋的应力和应变;fyc和fuc分别为锈蚀钢筋的名义屈服强度与名义极限强度。
式中:fy0和fu0分别为未锈蚀钢筋的名义屈服强度与名义极限强度;αyc和αuc分别为锈蚀钢筋的名义屈服强度相对值和名义极限强度相对值;Fy0和Fu0分别为未锈蚀钢筋的屈服荷载和极限荷载;A0为未锈蚀钢筋的公称截面积。
式(3)中εyc、εuc、εhc和Δεc分别为锈蚀钢筋的屈服应变、极限应变、强化应变和屈服平台长度应变,可以分别描述为
式中:αεy、αεu和Δε分别为锈蚀钢筋的屈服应变相对值、极限应变相对值和屈服平台长度相对值;εy0、εu0和Δε0分别为未锈蚀钢筋的屈服应变、极限应变和屈服平台长度应变。
由式(3)~式(6)可知,要确定不同锈蚀程度钢筋的应力-应变曲线,首先需要确定锈蚀钢筋的力学性能参数(名义屈服强度fyc与名义极限强度fuc)和变形性能参数(屈服应变εyc、极限应变εuc、强化应变εhc和屈服平台长度应变Δεc),其核心和关键是建立钢筋锈蚀特征参数(包括质量锈蚀率ηm, a、最大截面锈蚀率ηs, max、最小剩余截面积Smin、最大锈蚀深度Dmax和面积点蚀因子Ps)与名义屈服强度相对值αyc、名义极限强度相对值αuc、屈服应变相对值αεy、极限应变相对值αεu和屈服平台长度相对值Δε之间的定量关系。
分析了钢筋的均匀锈蚀特征参数ηm, a和非均匀锈蚀特征参数(包括ηs, max、Ps、最大锈蚀深度相对值αd, max=Dmax/d、最小剩余截面积相对值αs, min=Smin/S0)与锈蚀钢筋的名义屈服强度相对值αyc和名义极限强度相对值αuc之间的相关性。分析结果表明,是否对光圆钢筋和带肋钢筋进行分类分析对结果影响不明显,但分类后需要针对两类钢筋分别建立不同的预测模型,不便于工程应用。综合考虑计算精度和简便性,本文暂不区分钢筋类别的影响。如图 4和图 5所示,随着ηm, a、ηs, max和αd, max的增加,αyc和αuc呈线性减小趋势。其中,αyc和αuc与ηm, a之间预测方程的斜率绝对值相对较大,说明随着ηm, a的增加,αyc和αuc下降趋势较为明显。同时,随着αs, min和Ps的增加,αyc和αuc呈线性增加趋势。其中,与αs, min相比,αyc和αuc与Ps之间预测方程的斜率相对较大,说明随着Ps的增加,αyc和αuc增大趋势更为明显。
此外,由于钢筋锈蚀特征参数与αyc和αuc之间的预测方程非常接近,为了便于使用,可以建立钢筋锈蚀特征参数与αyc和αuc之间的统一预测方程,如图 6所示。图 6中,k为名义屈服强度和名义极限强度在钢筋锈蚀前后的相对值(统称为名义强度相对值)。
与力学性能分析类似,暂不考虑钢筋类别的影响,均匀/非均匀锈蚀特征参数(包括ηm, a、ηs, max、αd, max、αs, min和Ps)与锈蚀钢筋的屈服应变相对值αεy和极限应变相对值αεu之间的相关性分别如图 7和图 8所示。由图 7和图 8可知,随着ηm, a、ηs, max和αd, max的增加,αεy和αεu呈线性减小趋势;随着αs, min和Ps的增加,αεy和αεu呈线性增加趋势。其中,αεy和αεu与ηm, a预测方程的斜率绝对值相对较大,说明随着ηm, a的增加,αεy和αεu的变化趋势更为明显。此外,与αεy相比,锈蚀特征参数(包括ηm, a、ηs, max、αd, max、αs, min和Ps)与αεu之间预测方程的斜率绝对值相对更大,说明锈蚀特征参数对αεu的影响更为显著。
钢筋的均匀/非均匀锈蚀特征参数(包括ηm, a、ηs, max、αd, max、αs, min和Ps)与锈蚀钢筋的屈服平台长度相对值Δε之间的相关性如图 9所示。由图 9可知,Δε随着ηm, a、ηs, max和αd, max的增加而逐渐降低,随着αs, min和Ps的增加而逐渐增大。锈蚀钢筋的屈服平台长度相对值Δε与锈蚀特征参数之间的定量关系,可以确定当屈服平台基本消失(即Δε=0)时各锈蚀特征参数的临界值,即当ηm, a、ηs, max、αd, max分别大于20%、60%和40%左右或αs, min和Ps分别小于40%和60%左右时,锈蚀钢筋的屈服平台基本消失。由此可见,随着钢筋锈蚀程度的增加,锈蚀钢筋的屈服平台会出现明显退化现象,对于锈蚀严重的钢筋,其屈服平台可能完全消失。
基于图 6~图 9中建立的均匀/非均匀锈蚀特征参数(包括ηm, a、ηs, max、αd, max、αs, min和Ps)与锈蚀钢筋的名义屈服强度相对值αyc、名义极限强度相对值αuc、屈服应变相对值αεy、极限应变相对值αεu和屈服平台长度相对值Δε之间的定量关系,结合式(4)~式(6),可以计算不同锈蚀程度钢筋的力学性能参数(名义屈服强度fyc和名义极限强度fuc)和变形性能参数(屈服应变εyc、极限应变εuc、强化应变εhc和屈服平台长度应变Δεc),见表 1。进而根据式(3)可以分别以上述锈蚀特征参数作为控制参数,建立锈蚀钢筋的应力-应变关系模型。从135根试验钢筋中任意选取5根光圆钢筋(如M1P12-A2、M1P12-A3、M1P8-C2、M1P12-B2、M2P12-A3)和5根带肋钢筋(如M9R16-C1、M3R16-C1、M2R16-B3、M1R12-A2、M1R16-C1)进行对比分析,应力-应变关系曲线的模型计算值与实验测试值的对比分析如图 10所示。由图 10可知,当钢筋锈蚀程度较低(如低于10%)时,可以利用均匀锈蚀特征参数(ηm, a)和非均匀锈蚀特征参数(ηs, max、αd, max、αs, min和Ps)作为控制参数描述锈蚀钢筋的应力-应变关系曲线;而当钢筋发生重度锈蚀(如ηm, a大于25%)时,利用非均匀锈蚀特征参数(ηs, max、αd, max、αs, min和Ps)作为控制参数可以更好地描述锈蚀钢筋的应力-应变关系曲线。同时,从文献[1, 14]中任意选取了5组以质量锈蚀率为特征参数的试验数据(编号为SD-10、SG-2、SG-6、R16-D和P12-E),结合文献[8]的应力-应变关系模型,对本文模型开展对比验证分析,如图 11所示。图 11中,P1表示本文型,P2表示文献[8]模型。由图 11可知,文献[8]模型和本文模型的预测结果与实验数据总体吻合较好,但是本文模型的应力应变曲线与实验数据更加接近,说明本文模型具有较好的预测精度和适用性。
1) 锈蚀钢筋的力学和变形性能参数与锈蚀特征参数之间的相关性较大,且随着质量锈蚀率、最大截面锈蚀率和最大锈蚀深度相对值的增加呈现线性减小趋势,而随着最小剩余截面积相对值和面积点蚀因子的增大呈现线性增大趋势。与锈蚀钢筋的力学性能参数相比,锈蚀钢筋的变形性能参数与锈蚀特征参数之间的离散性相对更大。
2) 随着钢筋锈蚀程度的增加,锈蚀钢筋的屈服平台会出现明显退化现象。当质量锈蚀率ηm, a、最大截面锈蚀率ηs, max、最大锈蚀深度相对值αd, max分别大于20%、60%和40%左右或最小剩余截面积相对值αs, min和面积点蚀因子Ps分别小于40%和60%左右时,锈蚀钢筋的屈服平台基本消失。
3) 当钢筋锈蚀程度较低(如ηm, a低于10%)时,可以利用均匀锈蚀特征参数ηm, a和非均匀锈蚀特征参数(ηs, max、αs, min、Ps和αd, max)作为控制参数描述锈蚀钢筋的应力-应变关系曲线;而当钢筋发生重度锈蚀(如ηm, a大于25%)时,利用非均匀锈蚀特征参数(ηs, max、αs, min、Ps和αd, max)可以更好地描述锈蚀钢筋的应力-应变关系曲线。