Cl-在钢筋混凝土板一维双向扩散的分离变量法解答
岳著文, 李镜培, 杨博, 邵伟, 吕韬
同济大学地下建筑与工程系 岩土及地下工程教育部重点实验室 上海 200092
通讯作者:岳著文, E-mail:yuezhuwen@126.com

岳著文,男,1986年生,博士生,研究方向为桩基工程及其耐久性研究

摘要

针对钢筋混凝土板双侧Cl-侵蚀,建立其一维扩散方程,采用分离变量法进行求解,得到广义解答,通过与误差函数和数值解答进行对比,证明该解答是古典解。通过算例分析得出,双向海水侵蚀导致Cl-扩散较快,但并未产生浓度叠加现象;混凝土板一侧表面Cl-浓度升高与另一侧产生浓度梯度,增加扩散动力,加快Cl-扩散;对于扩散系数不同导致的非对称扩散,Cl-浓度梯度在扩散系数变化处发生突变,但扩散浓度连续。

关键词: 混凝土; Cl-; 双向扩散; 分离变量法
中图分类号:TU528
Solve Chloride Ions Diffusion Problem by Separation Variable Method for Reinforced Concrete Slab in Marine Environment
YUE Zhuwen, LI Jingpei, YANG Bo, SHAO Wei, LV Tao
Key Laboratory of Geotechnical and Underground Engineering of Ministry of Education, Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China
Abstract

Concrete slab as a part of marine engineering structure should surly be suffered from chloride ions corrosion. The relevant chloride ions diffusion problem in a concrete slab immersed in seawater was solved by means of separation variable method in this paper. In comparison with the solutions obtained by error function and numerical calculation methods, we proved this generalized solution was the correct answer for the bilateral chloride ion diffusion in a concrete slab. According to calculation and analysis, we found that the bilateral diffusion would cause greater harm than one side diffusion, but concentration superposition would not occur; chloride ion diffusion rate would speed up when concentration gradient increase caused by the increasing chloride ion concentration at one side surface. A mutation of chloride ion concentration gradient might occur where diffusion coefficient changed, but the chloride ion concentration distribution kept continuous.

Keyword: concrete; Cl-; bilateral diffusion; separation variable method
1 前言

钢筋混凝土板应用于海洋工程中时,将受到双侧海水中Cl-同时侵蚀。Cl-在混凝土中输运方式主要是扩散作用,当Cl-扩散到钢筋表面并达到一定的阀值时[ 1],钢筋表面去钝化,钢筋发生锈蚀,在此过程中Cl-充当催化剂,并未被消耗,因此,钢筋锈蚀将持续进行,对钢筋混凝土产生极大破坏。对于Cl-在混凝土中的扩散行为,Callepari等[ 2, 3]首次提出采用Fick定律进行描述并给出一维封闭解答。Tang等[ 4]给出扩散系数依赖时间变化条件下的扩散解答。余红发等[ 5]建立多因素扩散模型并给出一维解析解。Suryavanshi等[ 6]提出扩散系数为常系数时的二维解答,以上解答均基于单一侧Cl-扩散和半无限空间假设。杨绿峰等[ 7]采用边界元法对矩形域邻边Cl-扩散进行解答。李秀梅等[ 8]采用分离变量法给出对称扩散情况下混凝土中Cl-解答,然而双向Cl-扩散往往并非对称扩散。国内外对单侧Cl-扩散研究较成熟,而对于双向Cl-扩散问题鲜有研究。

基于Fick定律,本文针对混凝土板双边界侵蚀条件,尝试采用分离变量法对其Cl-扩散进行解答,并结合算例分析,探讨混凝土板双向Cl-扩散特征。

2 混凝土Cl-扩散问题的传统解答

Cl-在混凝土中的输运机理非常复杂,本文主要考虑扩散输运。假定边界处Cl-浓度 Cs不变,初始时刻试件内部Cl-浓度为 C,用Fick第二定律描述该一维单向扩散问题为:

Ct=D2Cx2(0x<+)C(0x<+,t=0)=C0C(x=0,t)=Cs (1)(2)(3)

其中, C为扩散浓度,%; t为扩散时间,s; x为扩散深度,mm; D为扩散系数,mm2/s。

基于半无限大物体假设上述问题解答为[ 9]

C(x,t)=C0+(CS-C0)1-erf(x2Dt) (4)

其中, erf(x2Dt)为误差函数,可查表获得。

3 混凝土板Cl-扩散分离变量解答
3.1 混凝土板Cl-双侧海水侵蚀Cl-扩散描述

混凝土板受双侧海水侵蚀作用,将其看做一维扩散问题,Cl-扩散用Fick第二定律描述该问题:

Ct=D2Cx2C(0<x<l,t=0)=C0C(x=0,t)=Cin,C(x=l,t)=Cout (5)(6)(7)

其中, l为混凝土板壁厚,mm; Cin Cout分别为内外壁表面Cl-浓度,%。

对于壁厚较小的混凝土板,半无限假设显然不再适用,本文使用分离变量法对混凝土板Cl-扩散问题进行解答[ 10]

3.2 混凝土板Cl-扩散分离变量解答

先将边界条件齐次化,设

C(x,t)=V(x,t)+v(x) (8)

其中, v(x)可取 Cin -Cin-Coutlx,满足 v(0)=Cinv(l)=Cout。则方程变化为:

Vt=D2Vx2V(x,0)=C0-v(x)V(x=0,t)=0,V(x=l,t)=0 (9)(10)(11)

分离变量设 V(x,t)=X(x)H(t)代入式 (9) 中整理,并令其等于 -λ得:

H(t)DH(t)=X(x)X(x)=-λ (12)

于是上式可写成两个方程:

H(t)+λDH(t)=0X(x)+λX(x)=0 (13)(14)

先求解式 (14),当 λ0时,只有零解,则 λ>0时,令 λ=β2 ( β>0),式 (14) 的解答为:

X(x)=B1cos(βx)+B2sin(βx) (15)

由边界条件式 (11) 可得 X(0)=0X(l)=0,代入上式中,得:

0=B10=B2sin(βl) (16)(17)

为得到非零解答,令 sin(βl)=0,求得: βl=nπβn=nπln=1,2。因此得到相应的特征值和特征函数:

λ=λn=βn2=(nπl)2 (18)

Xn(x)=Bnsin(nπlx) (19)

将特征值代入式 (13) 中,求得:

Hn(t)=Anexp-(nπl)2Dt (20)

由式 (18)~(20) 得:

=Cnexp-(nπl)2Dtsin(nπlx) Vn(x,t)=Xn(x)Hn(t) (21)

应用叠加原理:

=n=1Cnexp-(nπl)2Dtsin(nπlx) V(x,t)=n=1Vn(x,t) (22)

其中, Cn=AnBn;将上式代入初始条件式 (10) 中,令 C0-v(x)=φ(x)得:

φ(x)=V(x,0)=n=1Cnsin(nπlx) (23)

式 (22) 两边分别乘以 sin(mπlx)m=1,2。利用正交性,求得:

Cn=2l0lφ(x)sin(nπlx)dx (24)

将式 (24) 代入式 (22) 中,

=n=12l0lφ(x)sin(nπlx)dx? V(x,t)=n=1Vn(x,t)

exp-(nπl)2Dtsin(nπlx) (25)

φ(x)=C0-v(x)=C0-Cin+Cin-Coutlx代入上式,得到解答:

C(x,t)=V(x,t)+v(x)=

n=12n2π2[(C0-Cin)nπ-(C0-Cout)nπ

cos(nπ)+(Cin-Cout)sin(nπ)]

exp-(nπl)2Dtsin(nπlx)+

Cin-Cin-Coutlx (26)

对称扩散条件下,内外侧壁表面Cl-浓度相同,即 Cin= Cout,并假定初始时刻混凝土板内Cl-浓度为零,即 C0=0代入式 (26),得到对称条件下混凝土板Cl-浓度解答:

C(x,t)=n=12Cin[cos(nπ)-1]nπ

exp-(nπl)2Dtsin(nπlx)+Cin (27)

以上解答为广义解答,衰减很快, n取100即可保证精度,但解答是否是古典解有待进一步验证。

4 有限元数值计算

本文采用ABAQUS6.10有限元软件进行数值计算[ 11],采用其自带的质量扩散分析 (mass diffusion analysis) 进行Cl-扩散计算,在不考虑温度和压力的条件下,计算控制方程为:

dCdt+xJ=0 (28)

其中, J为扩散通量,g/mm2s;

J=-sDϕx (29)

其中, s为溶解度; ϕ为归一化浓度, ϕ= C/ s,用其表示边界浓度。

计算单元采用热传导单元,数值分析的优点在于可以设定复杂的边界条件,并进行Cl-扩散计算。

5 计算参数选取
5.1 表面Cl-浓度

施惠生等[ 12]对杭州湾地区调查结果表明,浪溅区和水位变动区的混凝土表面Cl-浓度大约在0.6%~0.7%的水平,而大气区的混凝土Cl-浓度基本在0.4%以下,结合Bamforth[ 13]和Maage等[ 14]的研究成果,确定了混凝土表面Cl-浓度取值范围,见表1

Table1 Range of surface Cl- concentration at concrete surface (relative to concrete / %) 表1 表面Cl-浓度取值范围

对于实际工程计算,表面Cl-浓度可取大值,计算结果保守,偏于安全。因此,本文计算统一取表面Cl-浓度为0.9%。

5.2 扩散系数

Hobbs[ 15]给出海洋环境混凝土扩散系数的计算公式:

D=0.04(1166w/c)×10-6mm2/s (30)

其中, w/c为水灰比。

我国混凝土耐久性设计规范[ 16]中规定,对海工条件水下区混凝土,服役寿命要求50 a时,最大允许水灰比为0.42,采用上式计算出扩散系数为7.8×10-7 mm2/s。

5.3 初始Cl-含量

我国混凝土耐久性指南[ 17]中对新拌混凝土中Cl-含量做了相应规定,见表2

Table 2 Limited content of Cl- in green concrete (relative to binder / %) 表2 新拌混凝土Cl-含量限值

规范规定范围内,初始Cl-含量值已经很低,计算时可不予考虑。

6 算例分析

以混凝土板 (壁厚70 mm) 为例计算其双向海水侵蚀作用下Cl-扩散情况。扩散系数取7.8×10-7 mm2/s,计算分为对称扩散和非对称扩散两种情况,对称扩散条件下 Cin= Cout=0.9%;非对称条件下 Cin=0.9%, Cout=0.5%。

6.1 3种计算方法比较

对称扩散Cl-浓度扩散情况,见图1。当扩散时间为1 a时,双侧Cl-没有交汇,3种方法均精确拟合,因此,对于Cl-浓度未交汇前3种方法均计算准确。经过5和10 a后双侧Cl-发生交汇,此时传统解答结果偏离另外两种计算方法,偏离量随深度、时间呈增大趋势,而数值计算和分离变量法始终较为精确拟合,说明基于无限长物体假设的误差函数解不适宜浓度交汇的混凝土板Cl-扩散计算。

Fig.1 Comparisons of Cl- concentrations obtained by three different methods under the condition of symmetric diffusion图1 3种方法得到的对称扩散结果

6.2 对称扩散

6.2.1 对称扩散特征 对称扩散是理想状态下的混凝土板双向侵蚀状态,假定两侧海水环境不发生改变,可以看做是对称扩散情况。用分离变量法 (式 (27)) 分别计算1,5,10,20和30 a后的Cl-扩散情况,结果见图2。可以看出,随时间增加,Cl-扩散浓度不断增大,但浓度交汇前 (1 a) 和交汇后 (5,10,20和30 a) 均未出现浓度叠加现象。分析认为,扩散动力源于浓度梯度,对称扩散条件下,双侧Cl-扩散浓度在中间处等值交汇,此处并未产生浓度梯度,因此不会发生浓度叠加现象。

Fig.2 Symmetric solutions of separation of variables method图2 分离变量法对称解答

6.2.2 对称扩散与单侧扩散比较 对称扩散和单侧扩散计算均设定边界表面浓度为0.9%,取壁厚一半,计算结果见图3。可知,当Cl-未出现交汇 (扩散1 a) 时,对称扩散与单侧扩散计算结果精确符合,但当Cl-浓度交汇后 (扩散5,10和30 a),两者计算结果产生偏离,且随扩散时间延长,偏离量呈增大趋势。双侧扩散情况下较单侧扩散Cl-浓度累积的快,且时间越长,两者差异越大。

Fig.3 Cl- concentrations obtained under conditions of symmetric diffusion and unilateral diffusion图3 对称扩散与单侧扩散

6.3 非对称扩散

实际工程环境中,混凝土板外侧受到阳光照射造成内外侧产生温度差异,或者受潮汐的干湿循环作用导致混凝土板两侧Cl-扩散边界条件不同均可造成混凝土板两侧Cl-非对称扩散。

6.3.1 边界条件不同导致的非对称扩散 假定混凝土板一侧Cl-浓度0.9%,而另一侧海水Cl-浓度为0.5%,按此假定进行计算,结果见图4。这种边界条件引起的非对称扩散直接导致两侧产生浓度梯度,增加扩散动力,从这点讲,将产生比对称扩散更加快的Cl-累积,更早造成钢筋锈蚀损害。

Fig.4 Asymmetric diffusion caused by different boundary conditions图4 边界条件不同导致的非对称扩散

6.3.2 扩散系数不同导致的非对称扩散 混凝土板背阳一侧常年处于阴凉、潮湿环境,而向阳一侧受到日照作用造成混凝土板两侧温差较大。文献[ 18]显示,温度从20 ℃上升到30 ℃,扩散系数可增加一倍;温度从20 ℃下降到10 ℃,扩散系数则可减小一半。由此可见,扩散系数对温度变化极为敏感。温度对Cl-扩散的影响主要通过对扩散系数的修正来体现[ 12],其计算表达式为:

kT=exp[UR(1T0-1T)] (31)

式中: kT为影响系数; U为扩散过程的活化能,35000 J/mol; R为气体常数,8.314 J/(molK); T为温度,K;常取 T=293.15 K (20 ℃)。

现假定混凝土板从中间处一分为二,两部分具有较大的平均温度差异,从而导致不同的扩散系数,假定20 ℃时扩散系数 D=7.8×10-7 mm2/s,内侧温度10 ℃,外侧温度30 ℃。采用式 (18) 计算出相应的修正后扩散系数 Din=4.68×10-7 mm2/s, Dout=1.25×10-6 mm2/s。这种情况下,式 (5) 变为分段函数,分离变量法较难解答,现采用数值计算进行解答,计算结果见图5。由图可知,由于扩散系数不同导致两侧Cl-浓度扩散差异较大,扩散系数较大的外侧在中间处明显高于内侧,值得一提的是,Cl-浓度变化率 (浓度梯度) 将会在扩散系数变化处产生突变,且突变量随扩散时间增加呈减小趋势,但扩散系数变化处浓度依然连续。

Fig.5 Asymmetric diffusion caused by different diffusion coefficients图5 扩散系数不同导致的非对称扩散

7 结论

(1) 通过与数值解答和传统解答对比,证明分离变量法所得到的广义解是古典解,并证明基于无限大物体假设的误差函数解答不适用受双侧海水侵蚀的混凝土Cl-扩散计算。

(2) 扩散动力源于浓度梯度,双向扩散条件下,混凝土中Cl-浓度在交汇处不发生浓度叠加现象,而是向低浓度一侧继续扩散。

(3) 双侧扩散较单侧扩散情况下Cl-浓度累积的快,且时间越长,浓度差异越大。因此,实际工程中应避免双侧海水侵蚀,减少Cl-累积,提高混凝土板耐久性。

(4) 混凝土板两侧表面Cl-浓度不同时,将产生浓度差,增加扩散动力,加快Cl-累积,会更早造成钢筋锈蚀。

(5) 因温度等因素造成扩散系数差异导致的非对称扩散,Cl-浓度梯度将会在扩散系数变化处产生突变,造成该处两侧较大的浓度梯度差异,增加扩散速率,但在扩散系数变化处浓度连续。

参考文献
[1] Ann K Y, Song H. Chloride threshold level for corrosion of steel in concrete[J]. Corros. Sci. , 2007, 49(11): 4113-4133 [本文引用:1] [JCR: 3.615]
[2] Collepardi M, Marcialis A, Turriziani R. Penetration of chloride ions into cement pastes and concretes[J]. J. Am. Ceram. Soc. , 1972, 55(10): 534-535 [本文引用:1] [JCR: 2.107]
[3] Collepardi M, Marcialis A, Turriziani R. The kinetics of penetration of chloride ions into the concrete[J]. Il Cemento, 1970, 67(4): 157-164 [本文引用:1]
[4] Tang L P, Gulikers J. On the mathematics of time-dependent apparent chloride diffusion coefficient in concrete[J]. Cem. Concr. Res. , 2007, 37(4): 589-595 [本文引用:1]
[5] Yu H F, Sun W, Yan L H, et al. Study on prediction of concrete service life I—Theoretical model[J]. J. Chin. Ceram. Soc. , 2002, 30(6): 686-690
(余红发, 孙伟, 鄢良慧. 混凝土使用寿命预测方法的研究I—理论模型[J]. 硅酸盐学报, 2002, 30(6): 686-690) [本文引用:1]
[6] Suryavanshi A K, Swamy R N, Cardew G E. Estimation of diffusion coefficient for chloride ion penetration into structural concrete[J]. ACI Mater. J. , 2002, 99(4): 441-449 [本文引用:1] [JCR: 0.826]
[7] Yang L F, Chen Z, Wang Y, et al. Boundary element method for analysis of two-dimensional chloride diffusion in concrete[J]. J. Chin. Ceram. Soc. , 2009, 37(7): 1110-1117
(杨绿峰, 陈正, 王燚. 混凝土中氯离子二维扩散分析的边界元法[J]. 硅酸盐学报, 2009, 37(7): 1110-1117) [本文引用:1]
[8] Li X M, Wu F, Huang Z H. Analytical solution to chloride diffusion equation on concrete[J]. Concrete, 2009, (10): 30-33
(李秀梅, 吴锋, 黄哲华. 混凝土中氯离子扩散方程的解析解[J]. 混凝土. 2009, (10): 30-33) [本文引用:1] [CJCR: 0.587]
[9] Jin W L, Yuan Y S, Wei J, et al. Concrete Structure Durability Theory and Design Method in a Chloride Environment[M]. Beijing: Science Press, 2011: 122-124
(金伟良, 袁迎曙, 卫军. 氯盐环境下混凝土结构耐久性理论与设计方法[M]. 北京: 科学出版社, 2011: 122-124) [本文引用:1]
[10] Chen C S. Mathematical Physical Equation[M]. Beijing: Science and Education Press, 2008: 70-108
(陈才生. 数学物理方程[M]. 北京: 科学教育出版社, 2008: 70-108) [本文引用:1]
[11] Dassault Systèmes Simulia Corp. Abaqus 6. 10Analysis User's Manual, 2010 [本文引用:1]
[12] Shi H S, Wang Q. Research on service life prediction of marine concrete[J]. J. Building Mater. , 2004, 7(2): 161-167
(施惠生, 王琼. 海工混凝土使用寿命预测研究[J]. 建筑材料学报, 2004, 7(2): 161-167) [本文引用:1]
[13] Bamforth P B. The derivation of input data for modelling chloride ingress from eight-year UK coastal exposure trials[J]. Mag. Concr. Res. , 1999, 51(2): 87-96 [本文引用:1]
[14] Maage M, Helland S T, Carlsen J E. Chloride penetration in high performance concrete exposed to marine environment [A]. International RILEM Workshop on Durability of High Performance Concrete [C]. Ann Arbor : RILEM Publications SARL, 1994: 194-207 [本文引用:1]
[15] Hobbs D W. Aggregate influence on chloride ion diffusion into concrete[J]. Cem. Concr. Res. , 1999, 29(12): 1995-1998 [本文引用:1]
[16] GB/T 50476-2008. Code for durability design of concrete structures[S]
(GB/T 50476-2008. 混凝土结构耐久性设计规范[S]) [本文引用:1]
[17] CCES 01-2004 (2005 revised edition). Guide to durability design and construction of reinforced structures[S]
(CCES 01-2004 (2005年修订版). 混凝土结构耐久性设计与施工指南[S]) [本文引用:1]
[18] Oh B H, Jang S Y. Effects of material and environmental parameters on chloride penetration profiles in concrete structures[J]. Cem. Concr. Res. , 2007, 37: 47-53 [本文引用:1]