摘要
对碳化作用细观求解方法进行梳理改进,利用传输和反应过程状态变化的路径依赖特点,规避了传统方法大规模矩阵求解,提出了基于扩散路径映射的碳化细观模拟方法。方法基于Dijkstra路径算法对界面层影响进行修正,建立水泥浆体扩散路径算法,得到环境表面到水泥浆体的最短路径网络;通过水泥浆体代理模型,模拟水化反应、碳化反应以及材料孔隙率变化的耦合碳化过程,得到关键指标(氢氧化钙浓度)沿扩散路径的时变分布;最终依据相同的扩散路径长度进行结果映射,实现碳化作用中关键物质浓度随时间空间分布的高效精准预测。通过材料试验对比校验,扩散路径映射方法模拟结果能够准确反应混凝土碳化深度范围;与传统偏微分求解方法对比,扩散路径映射方法能够在细观尺度精准模拟关键物质空间分布情况,且求解效率提升20倍以上,存储空间需求降低为偏微分求解的2%以内。
近年来混凝土结构病害集中爆发,其中最为普遍的是混凝土碳化问
区别于传统混凝土宏观模拟中基于均质化假定对材料单元赋予统一的单元属性,混凝土细观建模方法将材料模拟为包含骨料、水泥浆体和界面层(interfacial transition zone, ITZ)的三相复合模
在现有混凝土材料细观碳化模拟研究的基础上,本文建立了基于扩散路径进行映射的混凝土结构层面细观碳化作用的快速模拟方法,简称为扩散路径映射(diffusion path dependent mapping, DPDM)方法,通过将结构细观模型扩散路径计算以及水泥浆体碳化作用求解相结合,实现了细观模型精度和宏观模型效率的优势互补。与传统基于偏微分方程(partial difference equation, PDE)求解的碳化作用模拟方法相比,该方法保持了高精度的局部碳化作用模拟,同时大幅提升了模拟效率并降低了结果存储的空间需求。
混凝土细观碳化过程可以分解为水泥浆体中矿物成分水化过程以及水化产物的碳化过程,其反应原理与现有宏观模拟方法相似,关键的差异在于细观尺度材料介质对于求解域的影响。矿物成分水化的化学反应过程主要包括硅酸二钙()、硅酸三钙()、铝酸三钙()以及铁铝酸四钙()消耗孔隙水分形成水合硅酸钙()凝胶以及氢氧化钙()晶体。这些反应主要发生于水泥浆体与ITZ介质中,与混凝土中骨料组分并不发生反应,其具体反应过程如下式所示:
(1) |
式中:为混凝土矿物成分i细观尺度的水化机理作用,考虑到骨料介质不参与水化反应作用,因此细观尺度的水化反应空间求解域为模型内部的水泥浆体范围,是混凝土模型范围的子集,在本研究中也将作为细观模型路径影响抽离的关键参数;为矿物成分i的水化反应速率。
在不考虑矿物成分在水泥浆体中微观不均匀分布的情况下,由
(2) |
式中:为混凝土中的水分传输通量,由于骨料致密可以视为不可渗透介质,因此水分扩散过程也局限于水泥浆体求解域,将细观介质空间分布以求解域的形式进行下标表征,可以有助于后续对相同求解域进行合并处理,以降低水泥浆体空间部分复杂影响的重复计算;为水泥浆体中的水分扩散系数;为水泥浆体中的含水量,为用于求解场域梯度的Laplace算子。
水分在环境表面的传输渗透影响材料孔隙相对湿度,涉及水泥浆体中水分的吸附与解吸附过程,具体计算方法可以参照潘子超
(3) |
混凝土碳化过程主要为大气环境中的二氧化碳扩散溶解于孔隙溶液,水化产物凝胶、晶体以及部分矿物成分和与孔隙溶液中的二氧化碳发生反应,生成碳化产物碳酸钙。其具体反应过程涉及二氧化碳气体溶解、氢氧化钙溶解,酸解中和反应以及碳酸钙生成,如下式所示:
(4) |
式中:、以及分别表示化学成分的气态、溶解态以及固态状态。
根据化学反应表达式可以发现,碳化反应进程的主要控制因素包括作为反应原料的水化产物、二氧化碳、水分的空间分布,二氧化碳的渗透扩散,二氧化碳的气态与溶解态平衡,以及二氧化碳的反应速率。其中二氧化碳的扩散过程同样可以采用
(5) |
式中:为混凝土中的二氧化碳的传输通量,同样由于骨料致密不可渗透介质,因此二氧化碳扩散过程也局限于水泥浆体求解域;为水泥浆体中的二氧化碳扩散系数;为水泥浆体中的二氧化碳浓度。
对碳化过程进行总结,如
(6) |
在混凝土碳化水化和碳化过程中另一项材料参数的场域变化也同样需要考虑水泥浆体的孔隙率指标。随着水化和碳化产物的生成,水泥浆体内部的孔隙率逐渐下降,并影响水分和二氧化碳的扩散系数,由于孔隙率变化同样局限于水泥浆体内部,因此孔隙率变化过程可以表述为
(7) |
式中:为混凝土孔隙率变化,其变化过程如
(8) |
综上,对于混凝土细观碳化问题的模拟可以表述为水化反应、碳化反应以及材料孔隙率变化过程的方程联列。如
(9) |
从目前碳化反应模拟研究来看,细观尺度求解域对于反应机理作用并不构成影响,因此本研究对求解空间域的影响进行分开考虑。将传统宏观模拟结果反映到细观水泥浆体网络中,关键问题在于反映骨料随机分布对于侵蚀扩散路径的影响,而上述结论表明,碳化过程主要涉及环境物质的Fick扩散渗透以及各类反应原料的消耗和产物的生成。在不考虑化学成分浓度堆积的影响下,Fick侵蚀扩散路径一般符合最短路径原则,因此可以对细观求解域进行路径求解,同时在宏观求解过程输出路径长度信息,最终基于相同路径长度实现混凝土细观碳化模型的时间空间域求解,路径求解用表示。
(10) |
按照前述方法对混凝土细观碳化模拟过程进行拆解,骨料颗粒空间随机分布形成的水泥浆体网络以及碳化过程的迂曲分布与反应扩散的机理无关,可以通过计算机路径求解算法直接计算路径网络;而细观碳化过程中物理化学反应机理与宏观方法本质是统一的,不考虑路径空间分布对水化碳化反应进行精细求解也是可行的。面向工程分析需求,
如
(11) |
针对水泥浆体中的水化反应、碳化反应以及孔隙率变化过程,通过
(12) |

图1 混凝土细观碳化DPDM方法模拟流程
Fig.1 DPDM simulation procedure for mesoscopic concrete carbonation
基于两个模型中一致的路径长度进行结果映射,其函数意义可以表述为
(13) |
(14) |
式中:为DPDM方法计算得到的关键物质浓度在局部位置随时间变化的浓度分布
该模型的求解主要包含3个步骤:①在结构截面内建立骨料颗粒空间分布,并映射到截面网格单元中,采用节点网络最短路径计算方法求解扩散路径;②基于水泥浆体材料属性建立代理模型,基于现有研究基础模拟水化、碳化的物理化学反应过程,获得关键物质在深度分布的时变结果;③根据具体分析需求,基于同样侵蚀扩散深度,将特定时间节点的关键物质浓度映射到结构截面内,即可得到关注时刻的结构全截面内碳化作用关键物质的场域分布特征。得益于模型求解中对于结构特征和复杂材料环境特征的分离处理,该模型对于真实混凝土结构和复杂环境条件均能够提供高效精准的解决方案。
混凝土结构细观碳化路径的求解主要包含以下流程:结构关键截面的选取,截面骨料颗粒建模,以及细观最短扩散路径计算。
研究表明,混凝土材料内骨料颗粒的几何形态和随机空间分布对于局部碳化效应影响显著,因此在截面细观建模中需要考虑实际工程材料采用的骨料配比参数、类型以及浇筑方案的影响。本研究采用基于真实桥梁结构混凝土材料截面图像进行改进的三参数骨料模
截面内细观扩散路径的求解建立在模型网格划分的基础上,对模型按照所需模拟的精度需求进行自然网格划分,按照单元节点与细观模型几何位置的关系将节点分类为骨料节点、水泥浆体节点以及ITZ节点。在可以侵蚀退化的水泥浆体和ITZ节点中,采用计算机中常用的Dijkstra最短路径算法进行路径计
通过代理模型对服役过程中水泥浆体物理化学反应过程进行模拟,主要包括反应传输的传输扩散模拟,和反应原料消耗及产物生成的计算求解。该部分算法的具体实现形式是在单个时间步中,基于Fick扩散定律的偏微分方程对自由物质(水分、二氧化碳)分布进行求解,并对其他物质(矿物成分,水化产物、剩余水分、)的浓度和指标(饱和度、孔隙率)的状态进行更新。
区别于传统混凝土碳化数值模拟,本方法中水泥浆体模型仅作为碳化作用物理化学反应模拟的代理模型,与真实结构形式和尺寸并不直接相关,只需要满足结构退化深度方向上的模拟需求即可。模型可以通过建立狭长的代理模型大幅降低求解难度。一般模型宽度方向尺寸只要满足2倍以上于水泥浆体单元尺寸即可;模型长度尺寸需要略大于截面各点到环境表面的深度以满足后续映射的需求,考虑到骨料对于物质传输路径的迂曲影响,一般可以设置该代理模型长度为结构截面深度的1.2倍以上,其代理模型求解和结果提取如

图2 水泥浆体代理模型碳化模拟结果示意图
Fig.2 Illustration of carbonation simulation results of cement mortar agent model
扩散路径映射方法中最终的碳化作用全截面模拟,是基于同样的扩散路径长度和所关注的时间节点,将代理模型中的关键物质参数沿深度分布结果映射到结构整体截面中,即可得到细观碳化作用在结构全截面的分析结果。
(1)对于真实混凝土结构,由于混凝土细观截面和扩散路径求解均为一次性计算,且骨料颗粒建模和Dijkstra运算的时间和存储消耗均与模型尺寸正相关,因此混凝土结构碳化分析在面临众多截面分析时,该部分的计算成本并不会随模拟面积呈现指数增长。
(2)对于基础设施面临的复杂多变服役环境,不同时间尺度的物理化学反应则通过代理模型进行,因此能够考虑较为精细的环境气候变化影响,且不会造成模型求解和结果存储的负担。
为了校验DPDM方法的模拟准确性以及对比传统方法在计算精度和求解效率上的提升,本研究开展了两方面的对比分析:①将DPDM方法与混凝土试块碳化试验进行对比,通过对比碳化深度预测结果说明该方法的模拟准确性;②将DPDM方法与细观偏微分求解方法进行对比,通过对比碳化反应关键物质浓度分布以及计算成本说明该方法在求解性能方面的优势。
混凝土试块碳化试验采用边长100 mm的标准立方体,使用工程C50混凝土材料配合比设计(
水泥 | 粉煤灰 | 砂 | 细碎石 | 粗碎石 | 水 | 减水剂 |
---|---|---|---|---|---|---|
497 | 0 | 787 | 600 | 400 | 165 | 13.42 |
在DPDM数值模拟中,采用同样的材料参数进行混凝土细观建模,具体多边形骨料的形状指标控制采用Ruan提出的多重抽样细观建模方法。考虑到ITZ薄层厚度为30 ~ 80 μm之间,模型采用100 μm的网格划分能够精细模拟ITZ界面层以及扩散路径。而对于混凝土各相组分在水化、碳化过程中物理化学属性通过砂浆试块碳化测试以及相关研究结论推算得到。代理模型求解过程中,水化过程求解以小时为时间步长,碳化过程求解按照天为时间步长。
试验碳化作用结束后,劈裂试块,对新鲜截面进行酚酞显色反应以及磨粉化学滴定。得到截面中关键化学成分的浓度随深度分布情况,对比DPDM方法映射得到的浓度分布,其中碳化结果空间分布情况如

图3 DPDM模拟结果试验对比校验
Fig.3 Experimental comparison and verification of simulation result of DPDM
为了校验DPDM方法相较于现有混凝土细观碳化过程模拟的优势,本研究采用与前文试验一致的模型尺寸和材料性能参数。其中细观模型同样采用100 μm的网格划分,针对为期2 d的干燥水化过程采用小时为时间步长,针对为期120 d的碳化过程采用天为时间步长。将传统细观模拟和DPDM模拟结果进行对比,可以得到不同水化和碳化阶段关键化学物质在截面内的浓度分布结果。如

图4 传统PDE方程求解结果与DPDM模拟结果对比分析
Fig.4 Comparison and analysis of solution result of traditional PDE and simulation result of DPDM
通过对同样位置的计算结果差异绝对值进行计算,对比结果之间的具体差异以及偏差主要发生的位置特征,具体计算方法如下式所示:
(15) |
式中:为两个模型在局部位置、时间的结果偏差;为DPDM方法计算得到CH浓度随时间空间分布规律;为偏微分方程求解方法得到模拟结果。
对比得到相对于传统细观方法的计算偏差分布,如
基于单个细观模型可以获得混凝土材料内部非均匀碳化的局部差异,骨料颗粒导致的侵蚀路径迂曲导致了最终材料内部碳化影响深度的变异特征,这也是混凝土结构层面的局部劣化随机涌现的根源。依据上述细观模型的DPDM模拟方法,从单个模型中即可建立完整的碳化深度概率密度分布模式,这是传统简化的宏观均质模型所不具备的。
为比较DPDM方法对求解效率提升和存储空间缩减的实际效果,本研究基于同样的计算软件平台(Matlab R2020a)和硬件配置(Intel I7 7900, 32 GB),对方法计算效率和存储空间需求进行对比分析。对5个混凝土截面模型求解,所得参数如
模型 编号 | 求解时间/s | 时间降低/% | 存储空间/MB | 空间缩减/% | ||
---|---|---|---|---|---|---|
PDE方法 | DPDM 方法 | PDE 方法 | DPDM 方法 | |||
M1 | 4432 | 116 | 97.38 | 201.65 | 3.69 | 98.17 |
M2 | 4560 | 99 | 97.83 | 200.06 | 3.76 | 98.12 |
M3 | 4570 | 101 | 97.79 | 200.34 | 3.70 | 98.15 |
M4 | 4309 | 134 | 96.89 | 203.52 | 3.72 | 98.17 |
M5 | 4574 | 110 | 97.60 | 202.06 | 3.56 | 98.24 |
(1)基于混凝土细观碳化过程的扩散路径依赖特性,通过单次扩散路径计算取代全模型迭代求解来适应混凝土结构截面分析需求,通过水泥浆体代理模型求解来考虑服役环境的精细变化,因此DPDM方法在大体积、长周期混凝土结构全断面细观碳化模拟中具备优势。
(2)通过对比材料碳化试验以及传统偏微分求解结果,DPDM方法能够较为精准地模拟混凝土材料水化效应和碳化过程,其中关键物质CH的浓度随空间和时间的变异性分布也能够精准模拟预测。在100 mm边长的正方形混凝土截面模型中,DPDM方法能够保持细观碳化模拟的精度,并将求解效率提升20倍的同时降低结果存储所需空间至传统方法的2%以内。
作者贡献声明
李 越:试验设计,程序设计,数据分析,撰写初稿。
阮 欣:项目构思,课题指导,提出论文框架,论文修改与审定。
尹志逸:参与试验设计与试验结果分析。
参考文献
顾祥林, 徐宁, 黄庆华, 等. 混凝土结构时间多尺度环境作用研究[J]. 同济大学学报(自然科学版),2012, 40(1):1. [百度学术]
GU Xianglin, XU Ning, HUANG Qinghua, et al. Temporal multi-scale environment actions for concrete structures[J]. Journal of Tongji University (Natural Science), 2012, 40(1):1. [百度学术]
彭建新, 张建仁. 考虑腐蚀参数时空变异性的RC结构顺筋锈胀随机分析[J]. 土木工程学报, 2013,46(9): 46. [百度学术]
PENG Jianxin, ZHANG Jianren. Stochastic analysis of corrosion-induced cracking for RC structures considering spatial variability of corrosive parameters under atmospheric environments [J]. China Civil Engineering Journal, 2013, 46(9): 46. [百度学术]
屈文俊, 陈道普. 混凝土碳化的随机模型[J]. 同济大学学报(自然科学版), 2007, 35(5): 577. [百度学术]
QU Wenjun, CHEN Daopu. Stochastic model of concrete carbonation [J]. Journal of Tongji University (Natural Science),2007,35(5):577. [百度学术]
田飞龙, 李国平, 张春雷. 一般大气环境中混凝土桥梁长期受力性能分析[J]. 同济大学学报(自然科学版),2015,43(4):483. [百度学术]
TIAN Feilong, LI Guoping, ZHANG Chunlei. Long-term mechanical performance analysis of concrete bridges in atmospheric environment [J]. Journal of Tongji University (Natural Science),2015,43(4):483. [百度学术]
李国平, 胡皓, 任才, 等. 桥梁混凝土结构接缝的耐久性能[J]. 土木工程学报,2018, 51(7):1. [百度学术]
LI Guoping, HU Hao, REN Cai, et al. Study on durability of joints in concrete bridge structures[J]. China Civil Engineering Journal,2018, 51(7): 1. [百度学术]
余波, 成荻, 杨绿峰. 混凝土结构的碳化环境作用量化与耐久性分析[J]. 土木工程学报,2015,48(9): 51. [百度学术]
YU Bo, CHENG Di, YANG Lüfeng. Quantification of environmental effect for carbonation and durability analysis of concrete structures[J]. China Civil Engineering Journal, 2015, 48(9): 51. [百度学术]
PAN Z, XIN R, CHEN A. A 2-D numerical research on spatial variability of concrete carbonation depth at meso-scale[J]. Comput Concr, 2015,15(2):231. [百度学术]
阮欣, 李越, 金泽人, 等. 混凝土二维细观骨料建模方法综述[J]. 同济大学学报(自然科学版) ,2018,46(5):604. [百度学术]
RUAN Xin, LI Yue, JIN Zeren, et al. Review of two-dimension meso-modeling method of concrete aggregate [J]. Journal of Tongji University (Natural Science) ,2018,46(5):604 [百度学术]
王宗敏. 不均质材料(混凝土)裂隙扩展及宏观计算强度与变形[D]. 北京:清华大学, 1996. [百度学术]
WANG Zongmin. Crack growth, computer strength and deformation of nonhomogeneous material (concrete) [D]. Beijing: Tsinghua University, 1996. [百度学术]
RUAN X, LI Y, JIN Z, et al. Modeling method of concrete material at mesoscale with refined aggregate shapes based on image recognition[J]. Constr Build Mater, 2019,204(20):562. [百度学术]
潘子超, 阮欣, 陈艾荣. 基于任意级配的二维随机骨料生成方法[J]. 同济大学学报(自然科学版) ,2013,41(5):759. [百度学术]
PAN Zichao, RUAN Xin, CHEN Airong. Simulation method of random aggregate in two dimension based on arbitrary gradation [J]. Journal of Tongji University (Natural Science), 2013,41(5):759. [百度学术]
RUAN X, PAN Z. Mesoscopic simulation method of concrete carbonation process[J]. Struct Infrastruct Eng, 2012,8(2):99. [百度学术]
JIANG H, TIAN Y, JIN N, et al. Effect of aggregates spatial distribution on three-dimensional transport of chloride ions in reinforced concrete[J]. Constr Build Mater, 2020,259(1):119694. [百度学术]
PAPADAKIS V G, VAYENAS C G, FARDIS M N. Physical and chemical characteristics affecting the durability of concrete[J]. Aci Mater J, 1991,88(2):186. [百度学术]
DIJKSTRA E W. A note on two problems in connexion with graphs[J]. Numer Math, 1959,1(1):269. [百度学术]