同济大学地下建筑与工程系 岩土及地下工程教育部重点实验室 上海 200092
中图分类号: TU528
通讯作者:
收稿日期: 2013-07-18
修回日期: 2013-07-18
网络出版日期: --
版权声明: 2014 《中国腐蚀与防护学报》编辑部 版权所有 2014, 中国腐蚀与防护学报编辑部。使用时,请务必标明出处。
基金资助:
作者简介:
岳著文,男,1986年生,博士生,研究方向为桩基工程及其耐久性研究
展开
摘要
针对钢筋混凝土板双侧Cl-侵蚀,建立其一维扩散方程,采用分离变量法进行求解,得到广义解答,通过与误差函数和数值解答进行对比,证明该解答是古典解。通过算例分析得出,双向海水侵蚀导致Cl-扩散较快,但并未产生浓度叠加现象;混凝土板一侧表面Cl-浓度升高与另一侧产生浓度梯度,增加扩散动力,加快Cl-扩散;对于扩散系数不同导致的非对称扩散,Cl-浓度梯度在扩散系数变化处发生突变,但扩散浓度连续。
关键词:
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.
Keywords:
钢筋混凝土板应用于海洋工程中时,将受到双侧海水中Cl-同时侵蚀。Cl-在混凝土中输运方式主要是扩散作用,当Cl-扩散到钢筋表面并达到一定的阀值时[
基于Fick定律,本文针对混凝土板双边界侵蚀条件,尝试采用分离变量法对其Cl-扩散进行解答,并结合算例分析,探讨混凝土板双向Cl-扩散特征。
Cl-在混凝土中的输运机理非常复杂,本文主要考虑扩散输运。假定边界处Cl-浓度Cs不变,初始时刻试件内部Cl-浓度为C,用Fick第二定律描述该一维单向扩散问题为:
其中,C为扩散浓度,%;t为扩散时间,s;x为扩散深度,mm;D为扩散系数,mm2/s。
基于半无限大物体假设上述问题解答为[
其中,
混凝土板受双侧海水侵蚀作用,将其看做一维扩散问题,Cl-扩散用Fick第二定律描述该问题:
其中,l为混凝土板壁厚,mm;Cin和Cout分别为内外壁表面Cl-浓度,%。
对于壁厚较小的混凝土板,半无限假设显然不再适用,本文使用分离变量法对混凝土板Cl-扩散问题进行解答[
先将边界条件齐次化,设
其中,
分离变量设
于是上式可写成两个方程:
先求解式 (14),当
由边界条件式 (11) 可得
为得到非零解答,令
将特征值代入式 (13) 中,求得:
由式 (18)~(20) 得:
应用叠加原理:
其中,
式 (22) 两边分别乘以
将式 (24) 代入式 (22) 中,
将
对称扩散条件下,内外侧壁表面Cl-浓度相同,即Cin=Cout,并假定初始时刻混凝土板内Cl-浓度为零,即C0=0代入式 (26),得到对称条件下混凝土板Cl-浓度解答:
以上解答为广义解答,衰减很快,n取100即可保证精度,但解答是否是古典解有待进一步验证。
本文采用ABAQUS6.10有限元软件进行数值计算[
其中,
其中,s为溶解度;
计算单元采用热传导单元,数值分析的优点在于可以设定复杂的边界条件,并进行Cl-扩散计算。
施惠生等[
表1 表面Cl-浓度取值范围
Table1 Range of surface Cl- concentration at concrete surface (relative to concrete / %)
Composition of cementitious material | Tide zone | Splash zone | Atmospheric zone |
---|---|---|---|
Cement | 0.6~0.7 | 0.6~0.7 | 0.35~0.45 |
Cement+Slag powder | 0.7~0.9 | 0.7~0.9 | 0.45~0.55 |
Cement+Fly ash | 0.7~0.9 | 0.7~0.9 | 0.45~0.55 |
Cement+Silica fume | 0.6~0.7 | 0.6~0.7 | 0.45~0.55 |
Cement+Slag powder+ Fly ash | 0.7~0.9 | 0.7~0.9 | 0.45~0.55 |
Cement+Slag powder+ Fly ash+Silica fume | 0.7~0.9 | 0.7~0.9 | 0.45~0.55 |
对于实际工程计算,表面Cl-浓度可取大值,计算结果保守,偏于安全。因此,本文计算统一取表面Cl-浓度为0.9%。
Hobbs[
其中,
我国混凝土耐久性设计规范[
我国混凝土耐久性指南[
表2 新拌混凝土Cl-含量限值
Table 2 Limited content of Cl- in green concrete (relative to binder / %)
Service environment | Limited content of Cl-/% | ||
---|---|---|---|
Reinforced concrete | Wet environment | General engineering | 0.15 |
Bridge foundation | 0.08 | ||
Dry environment | General engineering | 0.3 | |
Bridge foundation | 0.15 | ||
Prestressed concrete | 0.06 |
规范规定范围内,初始Cl-含量值已经很低,计算时可不予考虑。
以混凝土板 (壁厚70 mm) 为例计算其双向海水侵蚀作用下Cl-扩散情况。扩散系数取7.8×10-7 mm2/s,计算分为对称扩散和非对称扩散两种情况,对称扩散条件下Cin=Cout=0.9%;非对称条件下Cin=0.9%,Cout=0.5%。
对称扩散Cl-浓度扩散情况,见图1。当扩散时间为1 a时,双侧Cl-没有交汇,3种方法均精确拟合,因此,对于Cl-浓度未交汇前3种方法均计算准确。经过5和10 a后双侧Cl-发生交汇,此时传统解答结果偏离另外两种计算方法,偏离量随深度、时间呈增大趋势,而数值计算和分离变量法始终较为精确拟合,说明基于无限长物体假设的误差函数解不适宜浓度交汇的混凝土板Cl-扩散计算。
图1 3种方法得到的对称扩散结果
Fig.1 Comparisons of Cl- concentrations obtained by three different methods under the condition of symmetric diffusion
6.2.1 对称扩散特征 对称扩散是理想状态下的混凝土板双向侵蚀状态,假定两侧海水环境不发生改变,可以看做是对称扩散情况。用分离变量法 (式 (27)) 分别计算1,5,10,20和30 a后的Cl-扩散情况,结果见图2。可以看出,随时间增加,Cl-扩散浓度不断增大,但浓度交汇前 (1 a) 和交汇后 (5,10,20和30 a) 均未出现浓度叠加现象。分析认为,扩散动力源于浓度梯度,对称扩散条件下,双侧Cl-扩散浓度在中间处等值交汇,此处并未产生浓度梯度,因此不会发生浓度叠加现象。
6.2.2 对称扩散与单侧扩散比较 对称扩散和单侧扩散计算均设定边界表面浓度为0.9%,取壁厚一半,计算结果见图3。可知,当Cl-未出现交汇 (扩散1 a) 时,对称扩散与单侧扩散计算结果精确符合,但当Cl-浓度交汇后 (扩散5,10和30 a),两者计算结果产生偏离,且随扩散时间延长,偏离量呈增大趋势。双侧扩散情况下较单侧扩散Cl-浓度累积的快,且时间越长,两者差异越大。
图3 对称扩散与单侧扩散
Fig.3 Cl- concentrations obtained under conditions of symmetric diffusion and unilateral diffusion
实际工程环境中,混凝土板外侧受到阳光照射造成内外侧产生温度差异,或者受潮汐的干湿循环作用导致混凝土板两侧Cl-扩散边界条件不同均可造成混凝土板两侧Cl-非对称扩散。
6.3.1 边界条件不同导致的非对称扩散 假定混凝土板一侧Cl-浓度0.9%,而另一侧海水Cl-浓度为0.5%,按此假定进行计算,结果见图4。这种边界条件引起的非对称扩散直接导致两侧产生浓度梯度,增加扩散动力,从这点讲,将产生比对称扩散更加快的Cl-累积,更早造成钢筋锈蚀损害。
图4 边界条件不同导致的非对称扩散
Fig.4 Asymmetric diffusion caused by different boundary conditions
6.3.2 扩散系数不同导致的非对称扩散 混凝土板背阳一侧常年处于阴凉、潮湿环境,而向阳一侧受到日照作用造成混凝土板两侧温差较大。文献[
式中:kT为影响系数;U为扩散过程的活化能,35000 J/mol;
现假定混凝土板从中间处一分为二,两部分具有较大的平均温度差异,从而导致不同的扩散系数,假定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-浓度变化率 (浓度梯度) 将会在扩散系数变化处产生突变,且突变量随扩散时间增加呈减小趋势,但扩散系数变化处浓度依然连续。
图5 扩散系数不同导致的非对称扩散
Fig.5 Asymmetric diffusion caused by different diffusion coefficients
(1) 通过与数值解答和传统解答对比,证明分离变量法所得到的广义解是古典解,并证明基于无限大物体假设的误差函数解答不适用受双侧海水侵蚀的混凝土Cl-扩散计算。
(2) 扩散动力源于浓度梯度,双向扩散条件下,混凝土中Cl-浓度在交汇处不发生浓度叠加现象,而是向低浓度一侧继续扩散。
(3) 双侧扩散较单侧扩散情况下Cl-浓度累积的快,且时间越长,浓度差异越大。因此,实际工程中应避免双侧海水侵蚀,减少Cl-累积,提高混凝土板耐久性。
(4) 混凝土板两侧表面Cl-浓度不同时,将产生浓度差,增加扩散动力,加快Cl-累积,会更早造成钢筋锈蚀。
(5) 因温度等因素造成扩散系数差异导致的非对称扩散,Cl-浓度梯度将会在扩散系数变化处产生突变,造成该处两侧较大的浓度梯度差异,增加扩散速率,但在扩散系数变化处浓度连续。
[1] |
Chloride threshold level for corrosion of steel in concrete [J].
|
[2] |
Penetration of chloride ions into cement pastes and concretes [J].
|
[3] |
The kinetics of penetration of chloride ions into the concrete [J].
|
[4] |
On the mathematics of time-dependent apparent chloride diffusion coefficient in concrete [J].
|
[5] |
Study on prediction of concrete service life I—Theoretical model [J].
混凝土使用寿命预测方法的研究I—理论模型 [J].
|
[6] |
Estimation of diffusion coefficient for chloride ion penetration into structural concrete [J].
|
[7] |
Boundary element method for analysis of two-dimensional chloride diffusion in concrete [J].
混凝土中氯离子二维扩散分析的边界元法 [J].
|
[8] |
Analytical solution to chloride diffusion equation on concrete [J].
混凝土中氯离子扩散方程的解析解 [J].
|
[9] |
|
[10] |
|
[11] |
|
[12] |
Research on service life prediction of marine concrete [J].
海工混凝土使用寿命预测研究 [J].
|
[13] |
The derivation of input data for modelling chloride ingress from eight-year UK coastal exposure trials [J].
|
[14] |
Chloride penetration in high performance concrete exposed to marine environment [A]. International RILEM Workshop on Durability of High Performance Concrete [C].
|
[15] |
Aggregate influence on chloride ion diffusion into concrete [J].
|
[16] |
|
[17] |
|
[18] |
Effects of material and environmental parameters on chloride penetration profiles in concrete structures [J].
|
/
〈 |
|
〉 |