Damage Evolution in Concrete Interfacial Transition Zone with Ultrasonic Dynamic Load
-
摘要: 混凝土是由粗骨料、水泥砂浆以及过渡区域组成的三相非均质复合材料。混凝土过渡区域(Interfacial transition zone, ITZ)是三相中最薄弱的环节,且难以观测,对混凝土的宏观力学性能有着重要影响。基于ABAQUS的Dynamic/Explicit模块,建立了能反映混凝土基质、骨料形状和 ITZ等真实细观结构的有限元模型,并应用该模型研究了过渡区域对混凝土损伤破坏的影响。研究结果表明:粗骨料的形状对混凝土损伤性能有一定影响,当骨料形状为凸多边形时,其抗损伤性能最弱;混凝土的抗损伤性能随着ITZ强度的降低而减弱,当ITZ的强度高于砂浆的60%时,抗损伤性能逐渐增强;随着ITZ厚度的增加,混凝土的抗损伤能力减弱。Abstract: Concrete is a three-phase material composed of coarse aggregate, cement mortar and interfacial transition zone (ITZ). The ITZ is the weakest of the three phases and difficult to observe, but it has a significant impact on the efficiency of concrete crushing. In order to study the impact of ITZ on the damage performance of concrete crushing, the finite element model that reflects real mesoscopic structure of concrete matrix, aggregate shape, and ITZ was established on the Dynamic/Explicit model in ABAQUS. The results showed that the shape of coarse aggregate has a certain influence on the damage performance of concrete, and when the shape is convex polygonal, its damage resistance is the weakest. The damage resistance ability of concrete decreases with the decrease of ITZ strength. When ITZ strength is higher than 60% of mortar, the damage resistance ability gradually increases. As the thickness of the ITZ area increases, the damage resistance ability decreases.
-
Key words:
- concrete /
- interfacial transition zone /
- damage performance /
- ultrasonic crushing
-
随着社会高速发展,特别是我国“双碳目标”的提出,对清洁能源的需求与日俱增,氢能作为一种环保、高效的新能源受到广泛关注。考虑到氢气管道建设成本高,利用现有的天然气管道,向甲烷中掺入一定量的氢气进而提升混合燃料热值,已成为一种经济可行的解决方案[1]。值得关注的是,近年来天然气输运事故屡有发生,若向其掺入氢气,有可能进一步加剧爆炸的危险性。倪靖等[2]讨论了不同掺氢比对甲烷-氧气爆轰特性的影响,发现掺氢后能够提高爆轰波的传播速度和爆轰敏感性。余明高等[3]通过实验研究了障碍物对甲烷-氢气爆炸特性的影响,发现最大爆炸超压和火焰传播速度随着障碍物阻塞率以及氢气体积分数的增大而增大。Yu等[4]研究了掺氢对甲烷-空气预混火焰传播特性的影响,发现随着氢气含量的增加,火焰前沿速度和爆炸超压显著升高。
鉴于可燃气体在运输过程中的危险性,大量学者对可燃气体的抑爆问题开展了研究。Li等[5]发现CO2对甲烷爆炸的抑制效果好于N2。Luo等[6]研究BC粉对氢气-甲烷-空气预混气的抑制后发现,BC粉对掺氢比低的预混气体有较好的抑制作用。陈鹏等[7]发现当金属丝网层数大于3时,甲烷-空气预混气爆炸火焰经过金属丝网时会发生淬熄。徐海顺等[8]通过研究发现铝镁合金泡沫对甲烷-空气预混气体爆炸有较好的抑制作用,材料对传播火焰的影响机制主要体现在湍流促进和冷却抑制两方面。袁必和等[9] 研究了新型多孔聚丙烯复合材料对瓦斯爆炸的影响,发现多孔材料的抑爆性能受填充位置、材料内径以及填充长度的影响。An等[10]通过对比球形非金属材料和铝合金材料对可燃气体的抑爆结果,发现由于球形非金属材料具有较高的结构强度,因而抑爆效果较好。邵继伟等[11]研究多孔材料对可燃气体的抑爆效果后发现,组合型多孔材料在密闭容器管道系统内的抑爆效果更为突出。综合以上研究发现,非金属多孔材料对抑制甲烷爆炸具有良好的效果,但甲烷掺氢的爆炸特性较甲烷出现明显改变,并且前人对非金属多孔材料在甲烷掺氢抑爆方面的研究较少。传统的非金属材料存在易燃、爆炸后熔融物易黏附管道内壁等缺陷,严重限制了多孔非金属材料在可燃气体抑爆中的实际应用。球形抑爆材料中空多孔结构能够有效增大比表面积,扩大热损失,并且耐高温耐火焰,具有良好的阻火隔爆性质[12-13]。
本研究通过实验探究掺氢比对甲烷掺氢爆炸特性的影响,比较单一球形多孔非金属材料和组合球形多孔非金属材料的抑爆效能,并通过改变材料填充长度,分析最佳阻火抑爆性能的搭配参数,为球形多孔非金属材料在混合燃料阻隔防爆领域的应用提供理论支撑与实验依据。
1. 实验装置与方案
1.1 实验装置与流程
为研究掺氢比对甲烷掺氢预混气体爆炸的影响并测试球形多孔材料的抑爆性能,自主搭建了气体管道爆炸实验平台,如图1所示。实验平台主要包括爆炸实验管道、点火系统、配气系统、同步控制、高速摄像系统(拍摄速度4 000帧每秒)以及数据采集系统。
实验管道由3节内径60 mm、管壁厚度10 mm、外径80 mm、长度2 m的水平圆形管道以及1节内径60 mm、长0.5 m、配备可视窗的方形管道组成,各个管道之间通过法兰-螺栓连接,实验管道为封闭状态。实验管道共安装2台压力传感器(CYG508微型高压传感器),分别距离管道左端1.2和6.3 m。点火系统(距离管道右端0.5 m处)由高压点火器以及2根钨丝组成,实验中点火能采用20 J。配气系统通过气瓶连接配气仪进行配气,预混气体中可燃气体的体积分数为10%。通过高速摄像系统采集火焰图像。实验过程如下:首先,清洁管道,连接管道和仪器;然后,检查管道气密性,采用抽真空法配气,向管道中通入预混气体;接着点火触发,数据采集系统对实验数据进行采集;最后,排出废气,重复实验,每组实验重复多次。
1.2 实验材料及实验方案
本实验的主要材料为球形多孔非金属材料和聚氨酯(polyurethane,PU)材料,如图2所示。球形多孔非金属材料的主要材质为聚偏氟乙烯,化学性质稳定,具有很强的抗冲击性和耐高温性,并且球形材料为多孔结构,孔隙率大,比表面积大。组合球形内部填充材料主要为聚氨酯,聚氨酯材料的可塑性较强,对爆燃火焰的传播具有一定的抑制作用[14]。
最大爆炸压力(pmax)是评价爆炸强度以及材料抑爆性能的重要指标。为了判断掺氢比对可燃气体爆炸的影响以及单一球形多孔非金属材料和组合球形非金属材料的抑爆性能,设计如下实验方案:首先研究不同掺氢比x对pmax的影响,分别对掺氢比为0%、5%、10%的甲烷掺氢预混气体进行测试,记录管道内爆炸压力并进行对比分析;在此基础上,研究多孔抑爆材料对pmax的影响,选取pmax最大的预混气体,分别对单一球形材料和组合球形材料进行抑爆实验,填充于距管道左端1.5 m的位置,填充长度分别为20、30、40 cm,对不同填充长度的同一材料以及相同填充长度的不同材料进行对比分析,判断其抑爆性能。
掺氢比x表示为
(1) 式中:
、 分别为甲烷-氢气混合气体中甲烷和氢气的体积分数。2. 实验结果及分析
2.1 掺氢比对爆炸压力的影响
管道中不同掺氢比对预混气体爆炸压力随时间的变化曲线如图3所示。从图3可以看出,1号压力传感器所采集的压力峰值随着氢气体积分数的增大而增大。不同掺氢比对火焰图像的影响如图4所示。对比发现,随着掺氢比的提高,火焰的亮度明显变强,甲烷-空气预混气体掺氢后火焰传播特性及爆炸特性发生显著改变,氢气对甲烷爆炸具有促进作用[15]。
预混气体爆炸压力是反映爆炸破坏效应的重要指标,爆炸压力越大,爆炸造成的后果往往越严重。掺氢比对预混气体最大爆炸压力pmax和最大爆炸压力上升速率(dp/dt)max的影响如图5所示。当掺氢比为0%、5%、10%时,1号压力传感器处的pmax分别为200、235、245 kPa,(dp/dt)max分别为2 250、2 875、3 250 kPa/s。与未掺氢(掺氢比为0%)相比,掺氢比为10%时,1号压力传感器处的pmax提高22.50%,(dp/dt)max提高44.44%,2号传感器处的(dp/dt)max提高26.67%,说明掺氢比较低时,pmax和(dp/dt)max均随着氢气含量的增加而增大,爆炸强度增大。这是由于氢气的活性高于甲烷,反应比甲烷剧烈,发生爆炸反应的时间更短,随着掺氢比的增加,单位体积内氢气的量增多,氢气比例的升高增强了能量释放的集中程度,气体燃烧速率加快,缩短了热量损失时间,导致爆炸压力上升速率提高[16]。
2.2 单一球形材料对爆炸压力的影响
为研究球形非金属材料对甲烷掺氢后预混气体爆炸的抑制作用,采用爆炸压力较大的掺氢比为10%的预混气体进行对照。空管道与填充不同长度的抑爆材料对爆炸压力的影响如图6所示,其中0 cm代表无抑爆材料。由图6可知,填充抑爆材料后,各个工况下的pmax均小于无抑爆材料。1号、2号压力传感器测得无抑爆材料的pmax分别为245、200 kPa。当抑爆材料的填充长度分别为20、30、40 cm时,1号压力传感器处的pmax分别为205、185、175 kPa,对比无抑爆材料时的pmax分别降低了16.33%、24.49%、28.57%,爆炸压力得到有效抑制。这是由于球形抑爆材料是多孔中空塑料球形结构,有较大的比表面积,并且球形材料把管道空间切割成多个小空间,增大了反应面与球形材料表面的接触面积,促进了球形材料与火焰的热交换,增加了热损失,并且冲击波通过球形材料时,被切割的反应面与后面的球形材料发生碰撞导致能量损失,进而使爆炸压力降低[12]。2号压力传感器处各个工况下的pmax相较于无抑爆材料分别降低5.00%、5.00%、12.50%,抑制效果较弱。这是由于随着火焰传播,火焰通过球形抑爆材料后再无障碍物,燃烧速率加快,火焰到达2号压力传感器时已充分反应,致使降压效果不明显。
不同填充长度对(dp/dt)max的影响如图7所示。在低填充长度下,(dp/dt)max有所上升,当填充长度为20 cm时,1号压力传感器处的(dp/dt)max提高7.69%,2号压力传感器处提高36.84%。重复实验均显示(dp/dt)max提高,这是由于当球形材料的填充长度较低时,火焰在穿越过程中的湍流度增大,燃烧传质传热进程加快,燃烧强度增大,火焰传播速度进一步加快[8]。随着填充长度的增加,火焰被分离成多个离散的湍流火焰,最大爆炸压力上升速率逐渐降低,填充长度为40 cm时,相较于无抑爆材料,1号压力传感器处的(dp/dt)max降低33.85%。
2.3 组合球形材料对爆炸压力的影响
不同长度的组合球形材料对爆炸压力的影响曲线如图8所示。从图8可以看出:与无抑爆材料相比,填充组合球形材料后pmax明显降低,说明组合球形材料对管道内可燃气体爆炸有较好的抑制效果。当填充长度分别为20、30、40 cm时,1号压力传感器测得的pmax分别为160、155、120 kPa,与无抑爆材料相比,分别衰减了34.69%、36.73%、51.02%,相较于单一球形非金属材料,抑爆性能分别提升了112.43%、49.98%、78.58%;2号压力传感器测得的pmax为135、135、110 kPa,与无抑爆材料相比,分别衰减了32.50%、32.50%、45.00%。可以看到,当填充长度为20和30 cm时,抑爆性能差别不大,这是多孔材料冷却抑制与障碍物加压共同作用的结果: 当火焰经过球形多孔抑爆材料时,材料会吸收热量从而对爆炸产生抑制效果,材料的填充长度越长,抑制效果越好;球形材料在管道中相当于障碍物,火焰经过障碍物时层流会转变成湍流,增加障碍物数量能够明显增大火焰传播过程中的湍流强度,导致压力上升[3]。当填充长度较低时,材料冷却抑制效果与障碍物加压效果相互抵消,最终导致抑制作用较接近。
从图9中可以看出,(dp/dt)max也得到了较大程度的抑制,与无抑爆材料相比,填充长度为20、30、40 cm时,1号压力传感器处的(dp/dt)max分别衰减34.62%、46.15%、53.85%,2号压力传感器处的(dp/dt)max最高衰减了47.37%,抑制效果明显。与单一球形材料相比,组合球形材料的抑爆性能得到大幅增强。这是由于爆炸穿越填充材料时,被球形材料多孔结构离散后的湍流火焰和冲击波能量与球形多孔材料中的聚氨酯材料接触,聚氨酯材料粗糙的壁面消耗了链式反应中的自由基,阻碍燃烧链式反应的进行,并使部分火焰热量转移到聚氨酯材料中,在火焰与聚氨酯材料发生碰撞的过程中,聚氨酯材料与冲击波产生摩擦,导致部分冲击波能量转化为热量并消散掉[14]。
3. 结 论
基于自主搭建的气体爆炸管道平台,研究了不同掺氢比条件下甲烷掺氢的爆炸特性以及不同类型球形抑爆材料的抑爆性能,得到以下结论。
(1) 氢气对甲烷-空气爆炸具有一定的促进作用。在低掺氢比下,随着可燃气体中氢气体积分数的增加,火焰传播速率加快,最大爆炸压力和最大爆炸压力上升速率升高。
(2) 单一球形多孔非金属材料对可燃气体爆炸具有两方面作用:一方面为抑制作用,体现在填充材料后,球形材料会吸收能量,导致最大爆炸压力降低,并且填充长度越长,抑制效果越好;另一方面为促进作用,体现在填充长度较低时,会导致燃烧传质传热进程加快,火焰传播加速,导致最大爆炸压力上升速率上升。
(3) 组合球形多孔非金属材料与单一球形多孔非金属材料相比,抑爆性能大幅提升,降压效果明显,当填充长度为40 cm时,最大爆炸压力衰减51.02%。球形材料抑爆性能受填充长度影响,在冷却抑制与障碍物增压的共同作用下,填充20和30 cm的组合球形多孔非金属材料的抑爆效果相差不大。
-
表 1 混凝土各相材料参数
Table 1. Concrete material parameters of each phase
Material Elastic model/GPa Poisson′s ratio μ Bulk density t/m3 Expansion angle/(°) Flow potential offset ε Ratio of double and uniaxial compressive strength αf Constant stress ratio Kc Aggregate 30.00 0.167 2.6 Mortar 10.66 0.167 2.1 30 0.1 1.16 0.666 7 ITZ 0.167 2.1 30 0.1 1.16 0.666 7 -
[1] 陈惠苏, 孙伟, STROEVEN P. 水泥基复合材料集料与浆体界面研究综述(二):界面微观结构的形成、劣化机理及其影响因素 [J]. 硅酸盐学报, 2004, 32(1): 70–79. doi: 10.3321/j.issn:0454-5648.2004.01.013CHEN H S, SUN W, STROEVEN P. Interfacial transition zone between aggregate and paste in cementitious composites (Ⅱ): mechanism of formation and degradation of interfacial transition zone microstructure, and its influence factors [J]. Journal of the Chinese Ceramic Society, 2004, 32(1): 70–79. doi: 10.3321/j.issn:0454-5648.2004.01.013 [2] 徐晶, 王先志. 纳米二氧化硅对混凝土界面过渡区的改性机制及其多尺度模型 [J]. 硅酸盐学报, 2018, 46(8): 1053–1058.XU J, WANG X Z. Effect of nano-silica modification on interfacial transition zone in concrete and its multiscale modelling [J]. Journal of the Chinese Ceramic Society, 2018, 46(8): 1053–1058. [3] YANG C C, CHO S W. Approximate migration coefficient of percolated interfacial transition zone by using the accelerated chloride migration test [J]. Cement and Concrete Research, 2005, 35: 344–350. doi: 10.1016/j.cemconres.2004.05.038 [4] AQUINO M J, LI Z, SHAH S P. Mechanical properties of the aggregate and cement interface [J]. Advanced Cement Based Materials, 1995, 2(6): 211–223. doi: 10.1016/1065-7355(95)90040-3 [5] LEE K M, PARK J H. A numerical model for elastic modulus of concrete considering interfacial transition zone [J]. Cement and Concrete Research, 2008, 38(3): 396–402. doi: 10.1016/j.cemconres.2007.09.019 [6] 王怀亮, 宋玉普. 多轴应力状态下混凝土的动态强度准则 [J]. 哈尔滨工业大学学报, 2014, 46(4): 93–97.WANG H L, SONG Y P. A dynamic strength criterion of concrete under multiaxial stress state [J]. Journal of Harbin Institute of Technology, 2014, 46(4): 93–97. [7] 杜修力, 金浏. 考虑过渡区界面影响的混凝土宏观力学性质研究 [J]. 工程力学, 2012, 29(12): 72–79. doi: 10.6052/j.issn.1000-4750.2011.04.0216DU X L, JIN L. Research on the influence of interfacial transition zone on the macro-mechanical properties of concrete [J]. Engineering Mechanics, 2012, 29(12): 72–79. doi: 10.6052/j.issn.1000-4750.2011.04.0216 [8] 王哲. 沿应变路径准静态加载时混凝土的极限状态现象 [J]. 北京交通大学学报, 2010, 34(1): 30–34. doi: 10.3969/j.issn.1673-0291.2010.01.007WANG Z. Phenomena of concrete limit state under quasi-static loading along strain paths [J]. Journal of Beijing Jiaotong University, 2010, 34(1): 30–34. doi: 10.3969/j.issn.1673-0291.2010.01.007 [9] GUINEA G V, EL-SAYED K, ROCCO C G, et al. The effect of the bond between the matrix and the aggregates on the cracking mechanism and fracture parameters of concrete [J]. Cement and Concrete Research, 2002, 32(12): 1961–1970. doi: 10.1016/S0008-8846(02)00902-X [10] 马巍, 任建伟, 胡俊, 等. 基于不同加载制度的轻骨料混凝土动态冲击性能 [J]. 硅酸盐通报, 2019, 38(4): 974–982.MA W, REN J W, HU J, et al. Dynamical shocking property of light-weighting aggregates concrete based on impact loading regimes [J]. Bulletin of the Chinese Ceramic Society, 2019, 38(4): 974–982. [11] 阮欣, 李越, 金泽人, 等. 混凝土二维细观骨料建模方法综述 [J]. 同济大学学报(自然科学版), 2018, 46(5): 604–612.RUAN X, LI Y, JIN Z R, et al. Review of two-dimensional meso-modeling methods of concrete aggregate [J]. Journal of Tongji University (Natural Science), 2018, 46(5): 604–612. [12] 刘建南, 张昌锁. 过渡区界面对混凝土劈裂性能影响的试验与数值模拟 [J]. 科学技术与工程, 2018, 18(18): 269–274. doi: 10.3969/j.issn.1671-1815.2018.18.044LIU J N, ZHANG C S. Experiment and numerical simulation on the influence of interfacial transition zone on concrete splitting performance [J]. Science Technology and Engineering, 2018, 18(18): 269–274. doi: 10.3969/j.issn.1671-1815.2018.18.044 [13] YANG C C. Effect of the interfacial transition zone on the transport and the elastic properties of mortar [J]. Magazine of Concrete Research, 2003, 55(4): 305–312. doi: 10.1680/macr.2003.55.4.305 [14] 过镇海, 李卫. 混凝土在不同应力-温度途径下的变形试验和本构关系 [J]. 土木工程学报, 1993, 26(5): 58–69. doi: 10.3321/j.issn:1000-131X.1993.05.001GUO Z H, LI W. Deformation testing and constitutive relationship of concrete under different stress-temperature paths [J]. China Civil Engineering Journal, 1993, 26(5): 58–69. doi: 10.3321/j.issn:1000-131X.1993.05.001 [15] 刘海峰, 韩莉. 二维骨料随机分布混凝土的动态力学性能数值模拟 [J]. 高压物理学报, 2016, 30(3): 191–199. doi: 10.11858/gywlxb.2016.03.003LIU H F, HAN L. Numerical simulation of dynamic mechanical behavior of concrete with two-dimensional random distribution of coarse aggregate [J]. Chinese Journal of High Pressure Physics, 2016, 30(3): 191–199. doi: 10.11858/gywlxb.2016.03.003 [16] 蒋橙炜, 陈启东, 顾泽堃. 超声破碎混凝土的力学模型与仿真分析 [J]. 机械制造与自动化, 2019, 48(2): 84–88.JIANG C W, CHEN Q D, GU Z K. Mechanical model and simulation analysis of ultrasonic crushed concrete [J]. Machine Building & Automation, 2019, 48(2): 84–88. [17] WEIBULL W. A statistical distributions function of wide applicability [J]. Journal of Applied Mechanics, 1951, 18: 293–297. -