考虑线性分布应力影响的混凝土氯离子侵蚀模型

刘建文，彭立敏，雷明锋, 2



考虑线性分布应力影响的混凝土氯离子侵蚀模型

刘建文1，彭立敏1，雷明锋1, 2

(1. 中南大学 土木工程学院，湖南 长沙，410075；2. 中国建筑第五工程局有限公司，湖南 长沙，410004)

基于氯离子扩散系数与结构应力服从多项式函数关系的假定以及Fick第二定律，通过边界条件的奇延拓和扩散方程的傅里叶变换推导建立综合考虑结构截面线性分布应力的混凝土氯离子侵蚀模型，并给出其求解流程以及相关参数的确定方法，进一步结合既有混凝土氯离子侵蚀试验数据进行实例验证和对比分析。研究结果表明：所建立的氯离子侵蚀模型能够反映结构截面线性分布应力对氯离子扩散的影响，特别是在浅层区，计算的氯离子质量分数与既有成果较吻合，且更接近于试验值，模型是可靠的。

线性应力；Fick第二定律；氯离子扩散系数；傅里叶变换

任何工程结构的设计应满足耐久性的技术原则，特别是对于混凝土结构，耐久性是确保结构功能正常发挥以及达到预定使用年限的前提条件[1−2]。影响混凝土结构耐久性的因素很多，其中，最关键的是钢筋锈蚀。张誉等[3−4]认为在引起混凝土破坏的原因重要性排序中，钢筋锈蚀处于第1位。因钢筋锈蚀而导致混凝土耐久性不足所造成的建筑物破坏的比例和经济损失较大[4−5]。而造成混凝土中钢筋锈蚀的因素也很 多[4, 6]，其中，环境氯盐造成的损失占比最大，达30%以上。为此，人们针对氯盐环境中钢筋混凝土耐久性问题进行了大量研究，特别是人们针对氯离子在混凝土中的侵蚀机制问题，基于Fick定律相继提出了诸多修正、补充模型[7−9]。在这些模型中，绝大部分是通过修正氯离子扩散系数来开展研究的。如，金立兵[10]研究了氯离子扩散系数与侵蚀时间的关系，指出扩散系数随侵蚀时间的增加而减小，且服从指数函数变化规律。BITARAF等[11−12]从氯离子的扩散需要借助水通道的角度研究了相对湿度和温度对氯离子扩散系数的影响，发现氯离子扩散速度随相对湿度的增大而增大，随着相对湿度的降低而减小，并存在临界相对湿度；温度升高会加剧氯离子活动，使氯离子扩散系数增加。目前，人们针对氯盐环境下混凝土结构耐久性问题进行了大量研究，但由于该问题本身的复杂性以及研究者对耐久性理解的差异，目前仍有诸多问题有待解决[13−14]，如关于荷载作用对氯离子扩散的影响，当前的研究均考虑均布荷载作用，即在垂直氯离子侵蚀方向的应力为均匀分布，研究成果也仅仅能反映计算截面上的平均应力对氯离子扩散的影响，而不能全面表征其在截面各处的影响状态。为此，本文作者在既有研究成果的基础上，考虑结构中应力的线性分布特征，进一步探讨荷载作用下氯离子在混凝土中的扩散规律和计算方法。

1 荷载作用对氯离子扩散系数的 影响

氯离子侵蚀示意图见图1。当前，对于混凝土中氯离子扩散的研究，大多采用的是Fick第二定律，即

式中：ws为表面氯离子质量分数；w为侵蚀时间t时深度x处的氯离子质量分数；D为氯离子扩散系数。在研究过程中，D可为常数，也可为耦合温度、相对湿度、荷载及时间等因素的复杂函数。如袁承斌等[15]使用FQY25型盐雾箱、质量分数为5%的NaCl溶液，分别对无应力、压应力和拉应力这3组混凝土试件进行盐雾试验，拟合出应力与氯离子扩散系数Dσ的关系：

式中：0为无应力状态下的氯离子有效扩散系数；为混凝土应力(拉应力为正，压应力为负)；和为待定系数。

LU等[16]采用同样的试验方法得

式中：为混凝土极限压或拉应力。王元战等[17]通过试验分别研究了干湿交替、盐雾和浸泡环境下荷载与氯离子扩散系数的关系：

涂永明等[18]利用盐雾箱模拟沿海地区环境，研究了不同应力状态及水平下的氯离子传输速率，得

可见：荷载作用对氯离子在混凝土中扩散的影响不可忽略。当前的研究大多针对均布应力状态，而实际混凝土结构单纯受轴拉(压)的情况很少，绝大部分为压弯结构，即在满足平截面假定条件时，其截面应力为线性分布，见图1和下式：

式中：为侵蚀深度；为应力分布曲线斜率；为图1中1处应力。

2 考虑线性分布应力的氯离子侵蚀模型

2.1 模型推导

考虑荷载作用时，氯离子扩散系数与应力具有简单的多项式关系。为此，假设考虑荷载作用时氯离子扩散系数与应力满足式(2)，将线性分布应力函数式(6)代入式(2)并耦合时间因素得

对式(8)中的第1式求导并将式(7)代入整理，有

对式(12)进行傅里叶变换逆变换，并代入(，)，经整理得

进一步代入式(8)中第2~3式的初始条件，有

式中：erf()为高斯误差函数。

于是，将和代入式(17)并整理可得到考虑线性分布应力影响的氯离子扩散方程：

2.2 求解方法与计算流程

考察式(18)，其中待确定的参数有0，，，s和，其中时间衰减指数与混凝土配合比及其所处环境有关，一般可取为0.25~0.60[19−20]；0为无荷载作用时的氯离子扩散系数，可由无荷载作用(或荷载水平近似为0 N)的工况来确定[15]。再利用有荷载工况下的侵蚀试验结果确定，和s，并判断这些参数是否为常数，若不是，则分析该参数的影响因素(如溶液质量分数、侵蚀时间等)并确定其与相关因素的关系式；否则，表明该参数是与其他因素无关的常数。具体求解流程见图2。

图2 求解流程图

3 参数分析与确定

为进一步明确上述模型中各参数的确定方法及取值特征，本文参考文献[21]中关于盾构管片混凝土在荷载作用下的氯离子侵蚀试验结果，见表1~2。其试验过程见参考文献[21]，并以此作为基础数据确定 参数。

3.1 时间衰减指数m

时间衰减指数取决于混凝土的配合比及其所处的环境，如表3所示[22−23]。文献[21]中，盾构管片混凝土配合比中添加了粉煤灰和矿渣掺合料，管片的赋存环境应为水下区，因此，将取为0.6。

3.2 氯离子扩散系数D0

无荷载作用时的氯离子扩散系数0可由表2中参数确定，具体拟合结果见表4。从表4可见拟合相关性系数均在0.95以上，由此可将其平均值10.93作为无荷载作用时的氯离子扩散系数0的取值。

表1 氯离子质量分数分布

注：1和2位置见图1。

表2 荷载约为0 N时氯离子质量分数分布

表3 m取值

表4 D0拟合值

3.3 A，B和ws分析

确定=0.6和0=10.93 mm2/月，再利用表1中试验结果确定，及表面氯离子质量分数s，拟合结果如表5所示。

表5 A，B和ws拟合值

3.3.1的分析与确定

从表5及图3可见：在溶液质量分数一定时，不同侵蚀时间的大致相等，并且随着溶液质量分数的升高而增大，说明受溶液质量分数的影响较大，而时间对其影响很小。故可假定只是溶液质量分数的函数而与侵蚀时间和侵蚀深度无关，即=()，拟合曲线见图3。从图3可见：有随增大而趋近于某一稳定值。因此，与溶液质量分数的函数关系可用双曲线的形式来描述：

图3 溶液质量分数对A的影响

式中：为环境溶液氯离子质量分数；1和1为拟合参数，在本试验条件下分别取4.10和0.17。

3.3.2的确定与分析

由表5及图4可见：分布在0.38~0.46之间，其变化幅度较小。随溶液质量分数升高而线性递减的特性，如图6所示。由此，与溶液质量分数的函数关系式可用如下线性函数来拟合：

式中：2和2为拟合参数，在本试验条件下分别取0.46和0.47。

3.3.3 表面氯离子质量分数s的确定与分析

由表5及图5和图6可见表面氯离子质量分数s受侵蚀时间与溶液质量分数的影响较明显，且随着和的增大而递增，近似满足线性关系，故可假定s与和满足如下关系式：

图4 溶液质量分数对B的影响

式中：α3，β3和γ3为拟合参数，在本试验条件下分别取−0.52，12.40和0.30。

图6 溶液氯离子质量分数对表面氯离子质量分数的影响

4 模型验证

为验证模型的可靠性，以文献[21]中的试验结果为基础数据，采用文献[24]和本文推导的氯离子侵蚀模型进行计算，并进行对比分析。

当侵蚀时间为2月、溶液质量分数为10%，侵蚀时间为4月、溶液质量分数为15%，侵蚀时间为6月、溶液质量分数15%时，采用文献[24]和本文推导的氯离子侵蚀模型进行计算所得氯离子质量分数见图7。从图7可见：

1) 总体而言，这2种模型的氯离子质量分数计算结果较接近，曲线形态、发展趋势吻合程度较高，且与试验结果拟合性好。可见本文所推导的考虑线性分布应力的氯离子侵蚀模型是可靠的。

2) 相对而言，考虑了线性分布应力状态下混凝土氯离子侵蚀模型与试验数据符合程度更高，预测效果更准确，特别是在浅层区域(≤20 mm)。而事实上，在该区域内，结构应力也较大(见图1)，可见本文所推导的模型能够反映荷载(应力水平)对氯离子侵蚀的影响。

(a) w=10%，t=2月；(b) w=15%，t=4月；(c) w=15%，t=6月

5 结论

1) 基于荷载作用下氯离子扩散系数与结构应力满足多项式函数关系的假定，以Fick第二定律为理论依据，结合傅里叶变换，推导得到了考虑线性分布应力作用下，混凝土氯离子扩散侵蚀模型，并给出了其计算流程和相关参数的确定方法。

2) 在既有试验数据的基础上，通过与经典混凝土氯离子侵蚀模型对比验证，结果表明本文所推导的计算模型是可靠的，能反映结构线性分布应力对氯离子扩散的影响，计算结果也更接近试验值。

[1] 覃维祖. 混凝土结构耐久性的整体论[J]. 建筑技术, 2003, 34(1): 19−22. QIN Weizu. Holistic view of durability of concrete structures[J]. Architecture Technology, 2003, 34(1): 19−22.

[2] 邓敏, 唐明述. 混凝土的耐久性与建筑业的可持续发展[J]. 混凝土, 1999(2): 8−12. DENG Min, TANG Mingshu. Durability of concrete and the sustainable development of construction industry[J]. Concrete, 1999(2): 8−12.

[3] 张誉, 蒋利学. 混凝土结构耐久性概论[M]. 上海: 上海科学技术出版社, 2003: 4−5. ZHANG Yu, JIANG Lixue. Introduction to the durability of concrete structure[M]. Shanghai: Shanghai Science and Technology Press, 2003: 4−5.

[4] LEI Mingfeng, PENG Limin, SHI Chenghua. An experimental study on durability of shield segments under load and chloride environment coupling effect[J]. Tunnelling and Underground Space Technology, 2014, 42(5): 15−24.

[5] 王建强. 地铁隧道衬砌结构钢筋锈蚀及耐久性研究[D]. 西安: 长安大学公路学院, 2009: 1−2. WANG Jianqiang. Study on lining structure steel corrosion and durability of subway tunnels[D]. Xi’an: Chang’an University. School of Highway, 2009: 1−2.

[6] 张友红, 曹勇, 陈晓斌. 氯盐侵蚀隧道结构寿命预测及耐久性技术研究[J]. 湖南交通科技, 2009, 35(3): 151−155. ZHANG Youhong, CAO Yong, CHEN Xiaobin. Study on life prediction and durability technology of tunnel structure under chloride erosion[J]. Hunan Communication Science and Technology, 2009, 35(3): 151−155.

[7] 李积平, 潘德强, 田俊峰, 等. 海工高性能混凝土抗氯离子侵蚀耐久寿命预测[C]//土建结构工程的安全性与耐久性科技论坛. 北京, 2001: 1−6. LI Jiping, PAN Deqiang, TIAN Junfeng, et al. Prediction of durability of high performance concrete against chloride ion erosion[C]//Safety and Durability of Civil Engineering and Technology Forum. Beijing, 2001: 1−6.

[8] 余红发, 孙伟. 混凝土在多重因素作用下的氯离子扩散方程[J]. 建筑材料学报, 2002, 5(3): 240−247. YU Hongfa, SUN Wei. Diffusion equations of chloride ion in concrete under the combined action of durability factors[J]. Journal of Building Materials, 2002, 5(3): 240−247.

[9] 金伟良, 袁迎曙, 卫军, 等. 氯盐环境下混凝土结构耐久性理论与设计方法[M]. 北京: 科学出版社, 2011: 58. JIN Weiliang, YUAN Yingshu, WEI Jun, et al. Durability theory and design method of concrete structure in chloride salt environment[M]. Beijing: Science Press, 2011: 58.

[10] 金立兵. 多重环境时间相似理论及其在沿海混凝土结构耐久性中的应用[D]. 杭州: 浙江大学建筑工程学院, 2008: 128−135. JIN Libing. Multi-environmenta time similarity(METS) theory and its application in coastal concrete structural durability[D]. Hangzhou: Zhejiang University. School of Architectural Engineering, 2008: 128−135.

[11] BITARAF M, MOHAMMADI S. Analysis of chloride diffusion in concrete structures for prediction of initiation time of corrosion using a new meshless approach[J]. Construction & Building Materials, 2008, 22(4): 546−556.

[12] BERKE N S, HICKS M C. Predicting chloride profiles in concrete[J]. Corrosion (Houston), 1994, 50(3): 234−239.

[13] 雷明锋, 彭立敏, 施成华. 氯盐侵蚀环境和荷载耦合作用下盾构管片耐久性评价与寿命预计[J]. 中南大学学报(自然科学版), 2015, 46(8): 3092−3099. LEI Mingfeng, PENG Limin, SHI Chenghua. Durability evaluation and life prediction of shield segment under coupling effect of chloride salt environment and load[J]. Journal of Central South University (Science and Technology), 2015, 46(8): 3092−3099.

[14] 雷明锋, 彭立敏, 谢友均, 等. 侵蚀环境和荷载耦合作用下盾构管片氯离子扩散规律试验研究[J]. 现代隧道技术, 2014, 51(3): 98−104. LEI Mingfeng, PENG Limin, XIE Youjun, et al. Experimental study on the diffusion law of the chloride ion in a shield segment under coupled erosion and load[J]. Modern Tunnelling Technology, 2014, 51(3): 98−104.

[15] 袁承斌, 张德峰, 刘荣桂, 等. 不同应力状态下混凝土抗氯离子侵蚀的研究[J]. 河海大学学报(自然科学版), 2003, 31(1): 50−54. YUAN Chengbin, ZHANG Defeng, LIU Ronggui, et al. Diffusivity of chloride in concrete in different stress states[J]. Journal of Hohai University (Natural Sciences), 2003, 31(1): 50−54.

[16] LU Chunhua, WANG Hailong, JIN Weiliang. Modeling the influence of stress level on chloride transport in prestressed concrete[C]//Proceedings of the international conference on durability of concrete structures. Hangzhou, 2008: 239−245.

[17] 王元战, 田双珠, 王军, 等. 不同环境条件下考虑荷载影响的氯离子扩散模型[J]. 水道港口, 2010, 31(2): 125−131. WANG Yuanzhan, TIAN Shuangzhu, WANG Jun, et al. Chloride diffusion model of RC member in various marine environments considering loading effect[J]. Journal of Waterway and Harbor, 2010, 31(2): 125−131.

[18] 涂永明, 吕志涛. 应力状态下混凝土结构的盐雾侵蚀试验研究[J]. 工业建筑, 2004, 34(5): 1−3. TU Yongming, LÜ Zhitao. The experimental research on prestressed concrete structure under salt fog corrosion environment[J]. Industrial Construction, 2004, 34(5):1−3.

[19] 姚昌建, 金伟良, 王海龙, 等. 海工混凝土氯离子扩散系数随深度的变化规律[J]. 水利水运工程学报, 2008, 2008(4): 14−18. YAO Changjian, JIN Weiliang, WANG Hailong, et al. Change of chloride diffusion coefficient in marine concrete with depth[J]. Hydro-science and Engineering, 2008, 2008(4): 14−18.

[20] 李智, 彭立敏, 施成华. 管片混凝土变扩散系数氯离子侵蚀模型[J]. 现代隧道技术, 2014, 51(1): 105−109. LI Zhi, PENG Limin, SHI Chenghua. Model of chloride ion erosion in segment concrete considering variable diffusion coefficients[J]. Modern Tunnelling Technology, 2014, 51(1): 105−109.

[21] 雷明锋. 侵蚀环境下盾构隧道结构性能全寿命计算方法研究[D]. 长沙: 中南大学土木工程学院, 2013: 45−47. LEI Mingfeng. Research on calculation methods of life-cycle structure performance of shield tunnel in aggressive environment[D]. Changsha: Central South University. School of Civil Engineering, 2013: 45−47.

[22] COSTA A, APPLETON J. Chloride penetration into concrete in marine environment-Part II: Prediction of long term chloride penetration[J]. Materials and Structures, 1999, 32(6): 354−359.

[23] BE95-1347, General guidelines for durability design and redesign[S].

[24] MANGAT P, LIMBACHIYA M C. Effect of initial curing on chloride diffusion in concrete repair materials[J]. Materials and Structures, 1994, 27(7): 338−345.

(编辑 陈灿华)

A concrete chloride ion erosion model considering influence of linear distributed stress

LIU Jianwen1, PENG Limin1, LEI Mingfeng1, 2

(1. School of Civil Engineering, Central South University, Changsha 410075, China; 2. China Construction Fifth Engineering Division Co. Ltd., Changsha 410004, China)

Based on the assumptions of polynomial relation between chloride diffusion coefficient and stress along with the Fick’s second law, and with the boundary conditions of odd extension and diffusion equation in the Fourier transform, a concrete chloride ion erosion model, which took the structure section linear distribution of force into consideration, was established, and the solving processes and the determination method of the related parameters were given. With the existing experimental data of concrete chloride ion erosion, a practical example was given to verify the results. The results show that the chloride ion erosion model can reflect the impacts of the linear section stress distribution of structure on chloride ion diffusion, and the calculated mass fraction of chloride ion coincides well with the existing results and is closer to the test results, especially in shallow area, which verifies that the model is reliable.

linear distributed stress; Fick’s second law; chloride diffusion coefficient; Fourier transformation

10.11817/j.issn.1672−7207.2018.02.018

U456.2

A

1672−7207(2018)02−0393−08

2017−02−10；

2017−04−08

国家自然科学基金资助项目(51508575)；湖南省创新平台与科技人才项目(2015RS4006)；中南大学中央高校基本科研业务费专项资金资助项目(2017zzts523, 2017zzts154)；湖南省研究生科研创新项目(CX2016B050) (Project(51508575) supported by the National Natural Science Foundation of China; Project(2015RS4006) supported by the Innovative Talents of Science and Technology Plan of Hunan Province; Projects(2017zzts523,2017zzts154) supported by the Fundamental Research Funds for the Central Universities of Central South University; Project (CX2016B050) supported by Hunan Provincial Innovation Foundation for Postgraduate)

雷明锋，副教授，从事隧道与地下工程研究工作；E-mail：mingfenglei@aliyun.com