混凝土真实细观模型的生成及氯离子传输的数值模拟

来源:用户上传      作者:

　　摘要：基于改进的遗传算法和8邻域边界跟踪法，对混凝土断面的数字图像进行处理，得到二值图像并提取粗骨料的边界坐标;根据提取的坐标编写程序，生成浆体集料界面过渡区（ITZ），得到真实的混凝土细观模型;将得到的混凝土细观模型导入COMSOL软件，模拟海洋水下区氯离子侵入混凝土内部过程，得到不同时刻混凝土内部氯离子浓度云图。研究结果表明：使用的建模和模拟方法所得结果与长期实海暴露混凝土实验结果一致，可以用来研究和评价海洋环境下混凝土的耐久性能;ITZ的存在会加速氯离子向混凝土内部扩散，其厚度越大，扩散过程越快，界面区厚度增大1倍，混凝土表观氯离子扩散系数增大12.3%;对比真实混凝土细观模型与参数化生成的圆形随机骨料模型中氯离子传输模拟结果发现，圆形随机骨料模型中氯离子的浓度总是小于真实细观模型中氯离子的浓度。
　　关键词：数字图像处理;混凝土细观模型;氯离子传输;数值模拟
　　中图分类号：TU528文献标志码：A 文章编号：2096-6717（2019）06-0071-09
　　在细观尺度上，混凝土可视为由水泥砂浆、骨料及浆体集料界面过渡区（ITZ）组成，在细观尺度上进行研究可以较好地描述材料组成及其非均质性对混凝土力学性能和耐久性的影响，实现从微观结构到宏观性能的过度。在细观尺度上，混凝土数值建模方法有两种：一种是将混凝土骨料粒径、级配等参数化，通过程序生成混凝土的细观结构，在二维空间，其骨料形状可分为圆形、椭圆形和多边形，而在三维空间，其骨料可以分为球形、椭球形以及多面体。Wang等建立二维随机骨料生成及随机分布算法，并将其应用到混凝土的本构关系研究中。王彩峰等提出一种基于ANSYS三维随机骨料构筑方法，该方法可以生成任意凸多面体骨料;赵蕊等采用圆形骨料，构建了骨料、浆体集料界面区的混凝土模型，并基于COMSOL软件考虑温湿度耦合作用下对混凝土氯离子传输的研究。杜成斌等通过一次性形成所有同级骨料，再对其进行随机延拓，通过三维和二维生成任意形状的随机骨料，并编写2D-RAS和3D-RAS随机骨料生成软件，避免逐个投放骨料而导致计算效率低下的问题。在总结现有随机骨料模型的基础上，钟根全等在骨料干涉的判断中引入交线法（Intersect Line），生成了圆形、椭圆形以及内接凸多边形的混凝土随机骨料模型。另一种是基于数字图像处理实现对真实混凝土的数值重构。最常用的方法是采用混凝土的CT图像进行处理并实现混凝土的数值模型重构。Boschmann等运用图像分析技术更准确地分析出混凝土中氯离子相对于水泥的含量。姜袁等基于CT图像重构了二维混凝土细观模型。匡广平通过对CT图像处理实现了三维混凝土细观模型的重构。She等基于XCT技术表征了泡沫混凝土结构，并对其热性能进行了模拟。Nguyen等。重构了泡沫混凝土的三维细观结构，并用ABAQUS有限元软件进行了原位试验仿真。但是，CT图像的获取受到试件尺寸、形状和环境等限制，成本较高。
　　为重构真实的骨料、避免混凝土CT图像获取成本高等缺点，本文从混凝土细观尺度上人手，采用数码相机拍摄混凝土数字图像，通过对混凝土截面图像处理，提取粗骨料的边界坐标，实现真实混凝土细观模型的重构;在生成细观模型基础上编写程序构造混凝土浆体集料界面区（ITZ）;基于上述数值混凝土模型，研究其在海洋水下区的氯离子传输过程。
　　1混凝土数字图像处理及其有限元（FEM）模型生成
　　1.1混凝土数字图像预处理
　　图像预处理的目的是得到高质量的二值图像。混凝土数字图像处理过程如图1所示。首先，在混凝土断面喷洒酚酞试剂以增强骨料和砂浆的对比度;利用数码相机对100×100mm2的C50混凝土断面进行拍照，得到数字图像。数字图像是RGB类型的图像，包含红、绿、蓝3种基本颜色。考虑高像素图片经过处理之后再导人有限元软件会极大地增加软件计算量。因此，选用400×400分辨率进行图像后续处理。然后，把得到的数字图像导入MATLAB。最后，对图像进行灰度处理、分割、空洞填充、去噪以及形态学处理等。
　　1.1.1图像灰度化及其分割首先，将分辨率为400×400的RGB图像灰度化得到灰度图像（像素值为0～255），如图2（b）所示。基于灰度值，混凝土中的骨料已经被清晰地显示出来。然后，将得到的灰度图像进行图像分割，得到二值图像（像素值为0和1），水泥浆体被赋予0值，显示黑色，而骨料将被赋予1值，显示白色。二值化原理见式（1）。
　　传统阈值是在画出灰度图像的灰度直方图之后，人为地在灰度直方图的“波谷”中选出，具有较大的主观性。本文采用智能算法中的遗传算法，求解最优阈值。遗传算法的流程如图3所示，与传统遗传算法相比，引入邻域搜索算子Reverse对其改进，提高算法收敛速度。
　　分割完成的图像如图4（a）所示。从图4（b）中可以看到，在骨料内部有许多空洞，这是因为骨料的矿物成分较多，有些矿物成分和水泥浆体的颜色类似。为消除这些空洞，采用内部填充的方法，填充完成之后的图像如图4（c）所示，从图4（c）中可以看出，所有骨料已经被填充好，骨料和砂浆有明显的分界。
　　1.1.2图像滤波与形态学处理
　　图像在采集和传输的过程中，外界干扰产生的噪声使图像变得模糊。通过比较各种滤波方法的特点，采用中值滤波对图片进行去噪处理。在进行中值滤波时，其窗口的形状和尺寸直接决定着滤波结果。因此，选择方形窗口，对采用不同模板大小的中值滤波后的图像进行比较，如图5所示。从图5中可以看出，随着模板的增大，尽管越来越多的杂质被去除掉，但是骨料也会重合起来，为了便于边界提取，选择了2×2模板。从图5（a）中可以看出，经过中值滤波之后仍然存在较多杂质，且骨料边缘粗糙，因此，采用形态学的处理方法繼续对二值图像进行开运算处理。采取结构元素形状为圆形的模板，半径为2的像素矩阵，得到的图像如图6所示。与图5（a）对比发现，图像中大部分杂质被去除，而且骨料的边界变得光滑。 　　1.1.3分离和去除细骨料 混凝土中直径小于4.75mm的骨料为细骨料，将其转换成像素面积去除掉，如图7（a）所示，细骨料已经被全部去除。但在识别过程中，有些骨料的边界相互渗透，连接在一起，如图7（b）红色标注所示，考虑人工识别出来之后将其分离，但人工识别、分离耗时较长，在识别出相互连接的骨料之后，采用腐蚀算法将其分离。分离之后执行图1中程序4，再一次判断是否有细骨料存在;如果有，执行图1中程序3，再一次消除细骨料，最终得到的二值图像如图7（c）所示。
　　1.2混凝土中骨料边界提取以及界面生成
　　图像的数据有矩阵、链码、拓扑和关系结构，其中，链码是用来描述目标图像的边界。链码按照标准方向的斜率分为4向链码和8向链码。编写了8邻域（8向链码）跟踪算法，提取骨料坐标并在MATLAB画出其形状，如图8（a）所示。根据提取的坐标直接在MATLAB中编写程序生成ITZ，如图8（b）所示。这样就更真实地再现了实际混凝土的细观模型，从而为混凝土中离子传输的数值模拟提供依据。
　　2饱和状态下混凝土中氯离子传输模拟
　　2.1饱和状态下混凝土中氯离子传输模型
　　海洋水下区的混凝土结构，如跨海大桥的桥墩、基础等长期浸水，一般处于完全饱和状态。氯离子在饱和状态下的输运机制以扩散为主，其传输方程可以用Fick第二定律来描述，如式（2）所示。
　　目前，高性能海工混凝土大多掺加粉煤灰，根据表1扩散系数的龄期系数n=0.69。青岛年平均气温为12.7℃，其氯离子扩散系数的活化能取3500J/mol，因此，温度影响因子KT=0.69299。粉煤灰混凝土在水下区表面氯离子浓度Cs，根据DuraCrete模型，Ac取值为10.8，混凝土W/B为0.35，计算其表面氯离子浓度Cs=0.948%。模拟100年内青岛海洋水下区混凝土中氯离子传输过程，如图10所示。从图10中可以看到，随着腐蚀龄期的增加，氯离子的渗透深度增大，混凝土表层氯离子浓度逐渐增大;由于粗骨料的氯离子扩散系数远小于砂浆及界面过渡区的扩散系数，所以，粗骨料会对氯离子的传输形成阻碍作用;氯离子通过大骨料以及骨料周边界面區传输后，在骨料后方会逐渐形成氯离子传输前锋。
　　为进一步表征氯离子随时间的侵入过程，选取混凝土左边缘为横坐标的原点，绘制氯离子浓度随时间演变线图，其结果如图11（a）所示。从图11（a）中可以明显看出，随着侵蚀时间的增加，在同一位置处的氯离子的浓度在不断增大;混凝土氯离子浓度仍然随深度增加而逐渐降低，且在某一深度逐渐趋于稳定。课题组针对C50混凝土开展的实海暴露试验结果如图11（b）所示。显然，模拟结果与实际结果具有较好的相关性。
　　2.3浆体集料界面区对混凝土中氯离子传输影响
　　在混凝土细观模型重构之后，假定浆体集料界面区厚度（ITZ）为0、30、60um，研究浆体集料界面区厚度对混凝土中氯离子传输影响。模拟海洋水下区腐蚀100a混凝土中氯离子传输过程，其氯离子浓度云图如图10、图12、图13所示。由于浆体集料界面区的扩散系数比砂浆和骨料的扩散系数大，随着浆体集料界面区厚度的增大，同一深度混凝土的氯离子浓度增大。
　　选取混凝土左边缘为纵坐标的原点，绘制100a不同界面区厚度混凝土的氯离子浓度分布曲线，如图14所示。在表面深度10mm以内，3条曲线基本重合，浆体集料界面过渡区对氯离子传输影响小;当距离表面大于10mm时，界面过渡区厚度增加，混凝土同一深度氯离子浓度明显增加，氯离子传输距离更远。
　　对图14中3条曲线按照Fick第二定律拟合得到无界面区，以及界面区厚度为30um和60um混凝土的表观氯离子扩散系数分别为0.03583×10-12、0.04575×10-12、0.0522×10-12 m2/s。由此可知，浆体集料界面区厚度越大，混凝土氯离子传输速度越快;界面过渡区提高一倍，其表观氯离子扩散系数增加了12.3%。
　　2.4两种混凝土细观模型对氯离子传输影响
　　对真实混凝土细观模型的粗骨料进行像素面积计算，得出粗骨料的体积率为41.23%。参照赵蕊等建立圆形随机骨料模型的方法，通过设定相同的浆体集料界面区厚度和骨料粒径分别为5～10mm、10～15mm、15～20mm，并控制骨料体积率为41%，编写算法生成随机骨料模型，如图15所示。导人COMSOL软件，给定与2.2节相同的计算参数，模拟100a内混凝土氯离子传输，其浓度分布云图，如图16所示。
　　比较图16和图10，采用两种细观模型获得混凝土中氯离子传输过程基本一致。因此，参数化生成的随机骨料模型也可用于模拟混凝土中氯离子传输过程。
　　为更好地对比两个混凝土细观模型对氯离子传输的影响，在腐蚀龄期相同的条件下，绘制氯离子浓度随时间演变线图，如图17所示。
　　由图17可知，当腐蚀龄期相同时，圆形随机骨料模型混凝土中氯离子浓度小于真实细观模型混凝土中氯离子的浓度。原因是：1）圆形随机骨料模型中粗骨料位置是随机的，在扩散边缘处的骨料体积越大，其粗骨料对混凝土内氯离子传输的阻碍作用越大;2）真实混凝土细观模型粗骨料比较“瘦长”，且在浇筑的时候骨料朝向随机，在相同位置处粗骨料体积率小，其对氯离子的传输阻碍作用要更小。因此，采用参数化骨料的细观模型虽然可以较好地模拟混凝土中氯离子的传输过程，但与真实骨料模型仍存在一定差异。
　　3结论
　　1）采用智能算法对真实混凝土断面的数码相机图形进行图像分割等预处理，可得到与原彩色图像吻合较好的二值图像。基于8邻域跟踪算法，对二值图像进行粗骨料坐标提取;编写算法生成粗骨料的界面过渡区（ITZ），实现真实混凝土断面三相细观模型的建立。
　　2）将真实混凝土断面细观模型导人COMSOL有限元软件，考虑龄期系数、表面氯离子浓度变化及温度对氯离子扩散的影响，实现海洋水下区混凝土氯离子传输过程模拟，模拟结果与长期实海暴露混凝土实验结果一致。
　　3）粗骨料会对混凝土氯离子的传输形成阻碍作用;氯离子通过大骨料以及骨料周边界面区传输后，在骨料后方会逐渐形成氯离子传输前锋。在距离表面小于10mm左右时，浆体集料界面区厚度对氯离子的传输影响较小;当距离大于10mm时，浆体集料界面区厚度增加将加速氯离子向混凝土内部传输;界面区厚度提高1倍，混凝土表观氯离子扩散系数提高12.3%。
　　4）通过真实混凝土细观模型与参数化生成的圆形随机骨料模型对比可以得出结论，圆形随机骨料模型可实现混凝土中氯离子传输过程模拟，但圆形随机骨料模型混凝土中氯离子浓度小于真实细观模型混凝土中氯离子的浓度。
转载注明来源:https://www.xzbu.com/8/view-15115232.htm