钢筋组合方式对混凝土保护层锈胀开裂影响

2019-09-05程旭东苏巧芝向恩泽梁晓岑王子栋张如林黄思凝

中国石油大学学报（自然科学版） 2019年4期

程旭东, 苏巧芝, 向恩泽, 梁晓岑, 王子栋, 张如林，黄思凝

(中国石油大学(华东)储运与建筑工程学院,山东青岛 266580)

对于跨海大桥、LNG储罐和海洋平台等氯盐环境中的在役钢筋混凝土结构,氯离子侵蚀引起的钢筋锈蚀是影响其耐久性的重要因素[1-3]。目前,关于钢筋锈蚀造成的混凝土开裂问题已经做了大量研究工作,但大多为针对单根钢筋的腐蚀膨胀行为进行分析[4-6]。实际上,混凝土构件中往往含有多根钢筋,并按照强度和构造要求对钢筋进行组合和分布,目前极少有实验和有限元模型分析不同钢筋组合对混凝土耐久性的影响规律。因此,笔者通过一个综合的扩散-力学模型研究多根钢筋非均匀锈蚀引起的保护层破坏行为,探讨分析不同钢筋组合对同一混凝土梁开裂路径、开裂时间、裂缝宽度和总截面损失率的影响。

1 数值模型

1.1 钢筋非均匀腐蚀模型

混凝土的氯离子侵蚀是一个复杂的物理化学过程。根据程旭东、Muthulingam等的研究[1,3],将温度传递、水分运输和氯离子扩散等过程进行耦合,各控制方程为

(1)

(2)

(3)

距离混凝土保护层最近的钢筋表面首先脱钝变为活化区,其余表面为钝化区。在活化区发生阳极铁氧化与阴极氧还原现象,相邻的阳极和阴极之间形成微电池。钝化区可以忽略阳极铁氧化,而阴极氧还原起主要作用。此外,在活化区和钝化区之间会形成腐蚀宏电池,这两种电化学腐蚀电池加速了钢筋的溶解。

在动态腐蚀演化过程中,活化区钢筋总体腐蚀电流密度为

itotal=imic+imac.

(4)

式中,imic和imac分别为微电池和宏电池电流密度,A/m2。

腐蚀电流密度的计算过程详见参考文献[7]-[9]。

由法拉第定律可得沿钢筋周向的钢筋半径减小值r1(θ,t),表达式[7]为

(5)

式中,θ为沿钢筋周向的角度,(°);AFe为铁相对原子质量;t1为钢筋腐蚀开始的时间,s;ZFe为阳极反应化合价;ρs为钢筋密度,kg/m3;F1为法拉第常数。

锈蚀产物的体积膨胀率n取值为2～6,与文献[10]中取相同的值n=3,铁锈膨胀层厚度ur(θ,t)表达式为

ur(θ,t)=(n-1)r1.

(6)

钢筋锈蚀率为

(7)

式中,R和r分别为钢筋半径和锈蚀后钢筋半径,m。

1.2 保护层开裂模型

基于损伤塑性模型描述混凝土的力学特性,该模型考虑了塑性应变引起的弹性刚度退化。混凝土采用CPE4R单元,将钢筋膨胀行为视为平面应变问题。各向同性损伤的混凝土本构方程为

σ=(1-G)D0(εt-εp)=D′(εt-εp) .

(8)

式中,G为弹性刚度损伤变量,MPa;D0为混凝土初始弹性刚度,MPa;εt为总应变;εp为塑性应变;D′为退化后的弹性刚度,MPa。

裂缝宽度计算为

(9)

式中,E0为混凝土的弹性模量,MPa;l为开裂单元尺寸,m。

氯离子扩散和钢筋非均匀锈蚀过程通过COMSOL有限元软件数值计算,采用ABAQUS自带的混凝土损伤塑性模型模拟保护层开裂,具体计算步骤为:①由式(1)、(2)和(3)可以得到氯离子浓度和钢筋周向各点的去钝化时间;②由式(5)计算钢筋半径损失量;③将由式(6)计算得到的锈胀层厚度作为位移荷载施加到钢筋周向节点上;④通过混凝土损伤塑性模型可得到保护层破坏行为。

2 模型验证及结果分析

2.1 模型验证

图1 膨胀层厚度对比试验结果Fig.1 Comparison of expansion thickness with experiment data

准确地预测腐蚀分布和膨胀层厚度是模拟混凝土开裂过程的关键,为了验证模型的合理性,该模型基于姬永生等[11]的实验参数,得到了沿钢筋周向分布的膨胀层厚度。如图1所示,模拟结果与实验结果吻合较好。按式(9)得到的裂缝宽度与Pedrosa等[12]的加速试验结果进行对比(图2),从图2中可以看出两者结果近似一致。但是当锈蚀深度比较大时,数值模拟结果要大于实验结果,原因在于实验中部分锈蚀产物会进入到孔隙和裂缝中去,使锈胀压力减小,从而延缓裂缝宽度增加。

图2 裂缝宽度对比试验结果Fig.2 Comparison of crack width with experiment data

2.2 结果分析

选一截面尺寸为200 mm×500 mm(宽×高)混凝土梁,取对称结构一半,如图3(以3根钢筋的模型为例)所示。模型上表面和两侧面均暴露在侵蚀环境中且边界自由,下表面限制竖向位移,图中,c为保护层厚度,d为钢筋直径,s为钢筋间距。混凝土弹性模量取28 GPa,混凝土抗拉强度取为1.78 MPa,泊松比为0.18。按照强度要求假定该混凝土梁的计算配筋率约为1%,在满足构造要求的情况下根据钢筋截面面积选取合适的钢筋根数和直径,保护层厚度为30 mm时,对应的钢筋组合 (钢筋总截面面积As/mm2)为2Φ25(As=982), 3Φ22(As=1 140),4Φ18(As=1 018);保护层厚度为40 mm时,对应的钢筋组合为2Φ25(As=982),3Φ20(As=942),4Φ16(As=804)。

图3 混凝土几何模型和边界条件Fig.3 Concrete geometric model and boundary condition

2.2.1 保护层开裂路径

以塑性应变的变化描述保护层的开裂路径,不同钢筋组合下的保护层开裂路径如图4所示。

图4 保护层的开裂路径Fig.4 Cracking paths of cover

从图4中可以看出,按照根数少且直径大的配筋方式时会形成4条表面裂缝,如果这些表面裂缝早早地贯穿整个保护层,很可能先出现保护层局部剥落的现象;相反,间距较小的两钢筋之间容易形成水平贯通裂缝,且只有2条侧表面裂缝,所以根数多且直径小的钢筋组合易引发保护层大面积层状剥落,破坏结构整体性。造成这种开裂路径的原因主要在于钢筋间距的变化,钢筋锈蚀膨胀会使上部混凝土变形,如果钢筋间距较小,膨胀应力会以压应力的形式作用在相邻钢筋周围的混凝土上,从而限制了保护层上表面裂缝的扩展(图5);但是混凝土为了消耗腐蚀膨胀增加的应变能,侧面裂缝和内部裂缝将开裂的更严重。通过数值模拟发现钢筋直径对开裂路径影响并不大,这与文献[13]中结论相同。另外,保护层越厚,上部混凝土需消耗更多的应变能使其开裂,从而限制表面裂缝的形成,这也是在3Φ20情况下未出现上表面裂缝的部分原因。

图5 钢筋锈胀对保护层表面开裂影响示意图Fig.5 Schematic diagram of influence of corrosion expansion on cover surface cracking

2.2.2 裂缝出现和贯通时间

裂缝出现和贯通时间是评估混凝土结构耐久性的关键。当混凝土的塑性应变大于0时,表示裂缝出现;当保护层表面的塑性应变向内扩展到钢筋或内部水平方向的塑性应变相互连通时,表示裂缝贯通。表1、2为不同保护层厚度时不同钢筋组合下裂缝出现和贯通时间,其中tin1、tsi1、tth1、tto1、tin2、tth2分别为角部钢筋周边内部裂缝出现时间、保护层侧表面裂缝出现时间、保护层侧表面裂缝贯通时间、保护层上表面裂缝出现时间、中部钢筋两侧水平裂缝出现时间和内部水平裂缝贯通时间。

从表1、2中可以看出,保护层从30 mm增大到40 mm时,裂缝出现和贯通时间延长了0.8～2 a;当保护层厚度一定时,根数少且直径大的钢筋组合也可以延缓水平裂缝和侧面裂缝的出现和贯通时间,如果间距足够大,表面裂缝很可能在内部裂缝相互贯通之前贯穿整个保护层。因此为了延长实际钢筋混凝土梁的开裂时间,建议适当提高保护层厚度并选择根数少且直径大的钢筋组合。

表1 不同钢筋组合下裂缝出现和贯通时间(c=30 mm)Table 1 Crack initiation and through time for different reinforcement arrangements(c=30 mm)

表2 不同钢筋组合下裂缝出现和贯通时间(c=40 mm) Table 2 Crack initiation and through time for different reinforcement arrangements(c=40 mm)

2.2.3 表面最大裂缝宽度

保护层表面裂缝宽度可以直观地反映混凝土破坏程度和锈裂情况。图6给出了不同钢筋组合下表面最大裂缝宽度随侵蚀时间的变化规律。

图6 不同钢筋组合下表面最大裂缝宽度Fig.6 Surface maximum crack width for different reinforcement arrangements

从图6中可以看出,侵蚀时间在18 a内,不同钢筋组合下的表面最大裂缝宽度均不会超过0.35 mm,而且裂缝宽度在2Φ25的情况下最小,在4Φ18和4Φ16的情况下较大,主要是因为较小的钢筋间距限制了上表面裂缝扩展,反而使侧面裂缝(表面最大裂缝)开裂更严重。因此根数少且直径大的钢筋组合对减小表面裂缝宽度更有利。

2.2.4 钢筋截面损失率

钢筋锈蚀不仅会引起保护层开裂,更直接的是造成钢筋有效截面面积的减少,从而降低结构的承载和变形能力。图7给出了不同钢筋组合下钢筋总截面损失率与侵蚀时间的关系,可以看出钢筋组合从2根钢筋变为4根钢筋时,截面损失率增大了约2.5倍,所以在保护层厚度和侵蚀时间一定的情况下,钢筋越少且直径越大的组合方式,其总截面损失率越小,说明对结构受力性能的影响也就越小。