Large eddy simulation of impacted obstacles' effects on premixed flame's characteristics
-
摘要: 障碍物在预混气体火焰传播过程中对其结构及传播特性造成较大影响，对火焰的加速和爆燃转爆轰过程(deflagration-to-detonation transition, DDT)起到直接的促进作用。通过障碍物条件下可视管道中甲烷/空气预混火焰传播实验，捕获其火焰微观结构变化。采用三维物理模型，采用壁面自适应局部涡黏模型(wall-adapting local eddy-viscosity, WALE)的大涡模拟(large eddy simulation, LES)，并用火焰增厚化学反应模型(thickened flame model, TFM)对实验过程进行重现。分析开口管道中预混火焰翻越障碍物后的复杂流场变化，并分析层流向湍流转变过程的特点。揭示了在障碍物影响下预混火焰扰动失稳现象的直接原因，是由障碍物引发的3个气流涡团同时作用而形成Kelvin-Helmholtz不稳定及Rayleigh-Taylor不稳定现象耦合作用所导致。Abstract: In the process of the premixed gas flame propagation, obstacles have vital influence on the flame's structure and propagation, and will enhance the flame acceleration and the DDT process. Through the methane/air premixed flame propagation experiment with obstacles in the visual pipe, the flame microstructure changes and the propagation characteristics were captured. By means of the three-dimensional physical model, a large eddy simulation mode (LES) with the WALE sub-grid scale models, and the thickened flame model (TFM) were used to repeat experiment process. The complex changes of the flow field were obtained when the premixed flame climbed over the obstacle, and the characteristics of the flow turbulence transition were analyzed. Finally, the direct cause of the premixed flame disturbance instability under the influence of the obstacles was revealed. It was induced by the coupling effect of the Kelvin-Helmholtz instability and the Rayleigh-Taylor instability, which in turn was affected by three vortexes as a result of the obstacle.
-
Key words:
- mechanics of explosion /
- thickened flame model /
- large eddy simulation /
- premixed flame /
- obstacles
-
砌体填充墙被广泛地应用于建筑物的外墙或钢筋混凝土框架内。但是, 由于它抗爆能力弱, 容易受到爆炸荷载的破坏, 且墙体的破坏与否对主体结构承受的爆炸荷载有一定影响[1], 因此, 填充墙的抗爆性能及其加固技术成为工程结构抗爆研究的重要方向。
近年来, 随着恐怖活动的增多, 为此开展了许多研究, 如:S.T.Dennis等[2]、S.C.Woodson等[3]进行了一系列比例为1/4的单向混凝土砌体墙的抗爆实验; C.D.Eamon等[4]、D.L.Thornburg[5]和J.S.Davidson等[6]进行了多次全尺寸的混凝土单向砌体墙抗爆实验; K.Martin[7]采用有限元分析方法研究了单向未加固砌体墙在面外静荷载作用下的失效模式; M.Wang等[8]运用有限元程序AUTODYN分别计算了实心砖墙和空心砖墙在外爆炸荷载作用下的破坏情况。
本文中, 先开展砌体填充墙的野外抗爆实验, 测得作用在墙上的爆炸荷载, 观测墙体的破坏形态与特征以及碎块的飞散和分布情况等, 然后利用LS-DYNA有限元软件, 对爆炸荷载作用下砌体填充墙实验模型的响应及损伤破坏进行详细分析与讨论。
1. 实验
砌体填充墙抗爆实验在野外进行, 模型采用浇注框架填充墙, 炸药在一定距离处爆炸, 测量不同装药量条件下墙体上的荷载及墙体变形和破坏情况。模型的平面布置如图 1所示, 该单层封闭框架模型为8×1跨, 层高1.5m, 在梁柱围成的框架处前面及两侧用砖墙填充, 以面向模型中心的墙体为实验墙, 两侧的砌体填充墙防止爆炸波绕射后对实验墙背面的影响。炸药在模型中心起爆, 实验模型如图 2所示。
1.1 模型制作
按设计图纸进行现场施工制作, 实验模型的基础、梁、板、柱整体浇注。砌体填充墙在模型建成进行填充砌筑, 墙的尺寸为2.00m×1.50m×0.19m。砖采用P型粘土空心砖, 尺寸为90mm×90mm×190mm, 强度等级为MU15, 砌砖上下错缝, 灰缝宽度为(10±2)mm, 砌墙砂浆采用M7.5水泥混合砂浆, 墙体与柱子间配拉结筋, 2根∅6mm沿柱高隔500mm布置, 拉筋伸入墙内500mm, 如图 3所示。
1.2 实验测量
测试内容包括墙面的压力、墙体的变形及砖墙破坏后的飞散情况。墙上的压力测点布置如图 4所示, 共8个测点, 所有的压力传感器均设置在铁管里, 铁管通过预埋砌筑在墙中。墙体的变形测点布置在墙背面, 共4个测点, 如图 5所示, 位移计固定在钢板上, 钢板锚固在两侧墙、顶板和地面上。
1.3 装药量
炸药安放在模型中心地面上, 装药为裸装TNT集团装药, 通过装药量变化实现不同比例距离上的爆炸, 各炮次的装药量和比例距离见表 1。
表 1 装药量Table 1. Explosive quantitiesNo. Q/kg Z/(m·kg-1/3) 1 0.2 10.00 2 3.9 3.72 3 3.9 3.72 4 8.6 2.87 5 34.2 1.81 6 4.8 3.48 7 4.8 3.38 8 4.8 3.48 9 4.8 3.48 10 10.7 2.66 11 21.2 2.12 12 30.0 1.89 2. 实验结果分析
2.1 墙上荷载
第1炮为试炮, 测得的墙上压力曲线如图 6(a)所示, 与TNT在空气中爆炸的经验曲线形状相同。第2炮后墙体没有破坏, 也没有发生明显变形, 图 6(b)是测得的经过滤波后的墙上压力曲线。从图 6可以看出, 冲击波先到达墙底部P5, 最后到达顶部P1。
表 2~3为第1、2炮各测点的冲击波的峰值压力、冲量。从表 2~3可以看出, 墙上的压力峰值和冲量都从墙底部到顶部、从中间到两边衰减。将实验结果与TM5-855-1[9]计算结果进行比较, 其中计算中假设墙体为刚体。可以看出, 实验测得的各测点的压力峰值、冲量与计算结果相差不大, 证实实验结果是可靠的。
表 2 第1炮各测点的压力峰值和冲量Table 2. The peak pressures and impulses of No.1测点 pm/Pa I/(MPa·ms) 实验 计算 实验 计算 P1 28.02 29.9 34.757 33.25 P2 28.80 30.3 42.193 33.75 P3 29.40 30.6 43.726 34.14 P4 30.00 31.0 45.560 34.40 P5 - 31.2 - 34.53 P6 29.41 30.6 45.706 34.09 P7 28.21 30.5 36.857 34.18 P8 29.73 30.3 39.543 33.74 表 3 第2炮各测点的压力峰值和冲量Table 3. The peak pressures and impulses of No.2测点 pm/kPa I/(MPa·ms) 实验 计算 实验 计算 P1 181.5 181 256.0 259.8 P2 - 185 - 271.8 P3 181.0 189 264.1 275.6 P4 - 191 - 278.2 P5 240.0 192 306.5 279.6 P6 231.0 188 303.4 275.2 P7 202.0 187 277.7 273.9 P8 218.5 185 235.4 271.9 2.2 墙的破坏
2.2.1第5炮
图 7是第5炮墙的破坏情况, 墙的变形明显, 迎爆面砖块表面发生剥落, 墙顶部与梁交界处全线透光, 这是因为实验墙是在框架浇注完后砌的, 墙顶部与梁交接处的砂浆是填塞进去的, 可能导致墙体顶部强度降低, 墙体的顶部破坏严重。墙体中上部砖块与砂浆分离发生内凹, 破坏模式类似于弯曲破坏, 右侧墙体整体内移2.5cm。墙体背面裂缝明显, 可以看到自上而下的裂缝, 裂缝以墙中心上下八字展开, 主要出现在砂浆处, 可见裂缝主要是由砂浆层的破坏引起的。从测点位移曲线(见图 7(c))可以看出, 最大位移发生在D2, 即最大变形发生在墙体的中部, 最大位移为88mm, 永久变形达到60mm, 其余3点的永久变形在30mm左右。
2.2.2第11炮
图 8是第11炮墙的破坏情况, 墙体破坏严重, 墙体中部内凹, 且在墙体上半部分出现30cm×30cm的孔洞, 右上部分与柱子连接的边界处出现孔洞, 拉结筋暴露。迎爆面的砖块有较多剥落散, 这是由爆炸冲击作用下墙体迎爆面砖块的局部被压碎引起的。由图 8(b)可以看出, 墙体背面砖块飞散和震塌现象严重, 抛射出的碎块以整块砖为主, 有约15块, 进一步表明墙的破坏主要是由砖与砖之间的砂浆层的破坏引起的。
2.2.3第12炮
图 9是第12炮墙的破坏情况, 墙体几乎全部破坏, 墙体迎爆面的地面上有少量砖碎片, 墙背面碎片很多, 且抛射距离很远, 最远距离达到34.2m。墙体右侧与柱之间的拉结筋暴露在外, 墙体左侧与柱子连接处发生松动。墙体施工工艺、边界条件以及重力的共同作用, 使墙体上部分破坏严重。
2.3 碎片分布
图 10为在不同药量作用下墙破坏后碎片的分布情况, 墙体的破坏形态与荷载的大小有关, 当药量较小时(见图 10(a)), 砖墙未形成孔洞, 仅背面有砖块震塌, 随着药量的增加, 墙体的破坏程度逐渐增大, 产生的碎片逐渐增多, 碎片的飞散距离也逐渐增远(见图 10(b)~(d)), 随着炮次的增多, 碎片数量增多, 飞散距离增远。实际上, 碎片的实际抛射距离跟落地后的反弹有关, 碎片与地面的冲撞问题非常复杂, 受各种不确定因素影响, 如地面刚度、平整度等。这里仅对实验得到的图片做分析。
3. 数值模拟
在野外实验分析的基础上, 利用LS-DYNA有限元软件, 建立精细化的有限元分析模型, 模拟填充墙在爆炸荷载作用下的响应和损伤破坏, 进一步分析砌体填充墙的抗爆特性、破坏形态和影响因素。
3.1 材料模型及参数的确定
在实验室对采用的砖和砂浆进行抗压测试, 得到砖的抗压强度为15.7MPa, 砌墙砂浆的抗压强度为7.2MPa。材料模型、参数以及墙的破坏准则、材料的失效准则见文献[10]。
3.2 有限元模型的建立和验证
3.2.1精细化有限元模型的建立
采用LS-DYNA有限元软件进行建模, 模型尺寸与实验墙尺寸相同, 将砖和砂浆分别采用不同的材料建立, 根据对称性取一半, 如图 11(a)所示。砖和砂浆均采用8节点6面体实体单元, 其中砖的单元尺寸为0.025m, 砂浆单元尺寸为0.005m, 拉筋采用Beam单元, 单元尺寸为0.025m, 砖、砂浆和钢筋之间均采用共节点, 整个有限元模型共有节点1 024 982个, 划分实体单元529 614个, 梁单元1 132个, 有限元网格图见图 11(b)~(c)所示。实验中TNT炸药离砖墙距离为5.871m, 通过改变药量实现比例距离的变化。
3.2.2特征点的动力响应
根据以上分析, 实验中炸药爆炸作用在墙上的荷载与TM5-855-1计算结果较好吻合, 因此采用LS-DYNA中*LOAD_BLAST关键字加载方法, 用CONWEP爆炸加载[9]模拟实验中炸药的爆炸情况。第3炮后, 墙正面顶部与梁连接处有细微裂缝, 图 12是各测点位移曲线的数值模拟结果和实验结果的对比, 墙的最大位移发生在D1, 实验测得的最大位移为4.39mm, 永久变形为2.23mm, 数值模拟的最大位移为5.09mm, 永久位移为1.3m, 两者结果基本一致。这表明, 采用*LOAD_BLAST关键字模拟墙上爆炸荷载的方法, 既能满足计算精度的要求, 也可以大大提高计算效率, 所以本文的数值模拟方法切实可行。
3.3 边界条件的影响
实验墙的左右两侧与柱子之间设置拉结筋, 通过数值模拟, 对比有无拉结筋时墙体的破坏模式的差异, 如图 13所示。墙体两侧有拉结筋时, 墙体的破坏主要集中在墙体的中部(见图 13(a))。无拉结筋的墙体破坏模式有些类似于整体弯曲破坏(见图 13(b))。有拉筋的墙体破坏区域比无拉结筋的墙体破坏区域小, 所以拉筋的存在可以增强墙的抗爆能力, 减少碎片的产生。
3.4 砌体填充墙破坏等级划分
实验时对同一面墙体进行了多次爆炸作用, 有损伤累积效应, 很难判别墙体的破坏形态在什么药量下产生。结合实验结果, 根据文献[10]中墙体的破坏等级划分, 确定不同等级实验墙的破坏药量。
通过数值模拟得出, 当药量Q < 15kg时, 墙的破坏程度为轻微破坏, 墙体仅出现细微裂缝; 当15kg < Q < 40kg时, 墙体发生中度破坏, 图 14为Q=37kg时墙的破坏情况, 由于墙体顶部与梁之间无拉结筋, 墙体会先在顶端与梁交界处出现较大横向裂缝, 再沿灰缝产生竖向不规则的裂缝, 墙体变形明显, 在墙体变形过程中的等效塑形应变分布图(见图 14(a))中, 红色区域为变形和裂缝较集中的地方, 在裂缝分布图(见图 14(b))中, 灰缝不连续部分为产生的裂缝。与对应的实验破坏情况(见图 8)比较, 数值模拟所得墙体的破坏形态与实验现象较好吻合。
数值模拟中, 当40kg < Q < 50kg时, 墙体发生严重破坏, 墙体发生多处不规则裂缝, 外界不大的扰动也有可能导致墙体的倒塌, 如图 15(a)中红色部分为裂缝。冲击波作用于墙上, 墙体中部弯矩最大, 且墙的两侧有拉筋, 因此墙体的破坏主要集中在墙的中部, 灰缝截面由于剪切滑移发生破坏, 使砖与砂浆分离, 在冲击波和重力作用下, 墙体背面部分产生砖块震塌。实验第8炮的墙体破坏如图 15(b)所示, 数值模拟结果与实验结果完全吻合。
当Q > 50kg时, 墙发生飞散。图 16为Q=58kg时墙的破坏过程, 墙体在爆炸荷载和重力联合作用下发生整体飞散破坏, 墙体中上部先发生破坏飞散, 背部砖块普遍震塌, 随着墙体中上部砖块的飞散, 墙体顶部的砖块也发生脱落和飞散, 仅留下下部分和侧面部分墙体没有倒塌。实验的破坏(见图 9)中, 墙的破坏范围很大, 数值模拟结果基本如实反映了实验观测现象。
本文中所确定不同破坏等级时的药量只是针对实验墙, 如果墙体尺寸发生变化, 针对不同的破坏等级所需的药量也会发生变化。
4. 结论
采用砌体填充墙抗爆的野外实验, 并与数值模拟相结合, 对填充墙的抗爆特性和破坏机制进行了研究, 主要结论有:
(1) 在外爆炸作用下, 墙上各点的冲击波压力峰值、冲量按照从墙底部到顶部、从中间到两边衰减, 实验结果与TM5-855-1计算结果较好吻合。
(2) 在外爆炸作用下, 墙体的裂缝不规则发展, 主要由砂浆层的破坏引起, 墙体的破坏先在背部发生震塌现象, 随着药量逐渐增大, 墙体发生飞散, 产生的碎片数量逐渐增多, 抛射距离也逐渐增远。
(3) 墙体的边界条件对墙体的破坏模式影响很大, 无拉结筋时墙体易发生整体弯曲破坏, 拉筋的存在可以增强墙的抗爆能力, 减少碎片的产生。
(4) 通过数值模拟, 确定了实验墙不同破坏程度的药量。当药量Q < 15kg时, 墙体发生轻度破坏; 当15kg < Q < 40kg时, 墙体发生中度破坏; 当40kg < Q < 50kg时, 墙体发生严重破坏; 当Q > 50kg时, 墙体发生飞散。
-
图 2 甲烷/空气预混火焰传播过程纹影图片
Figure 2. High-speed schlieren images of premixed methane/air flame propagation
图 3 LES模型计算三维预混火焰结构时间序列图
Figure 3. Sequence diagram three-dimensional premixed flame structure by LES model
图 4 火焰翻越障碍物时管道内流场结构示意图
Figure 4. The pipe flow field structure diagram when the flame climb over the obstacle
图 5 特征点处压力时程曲线实验与LES模型计算值对比
Figure 5. Comparison of pressure histories at characteristic point between experiment and LES
图 6 预混火焰翻越障碍物过程中管道内的速度矢量及流场变化图
Figure 6. Velocity vector and the flow field when premixed flame climbed over obstacle
-
[1] Ciccarelli G, Johansen C, Kellenberger M.High-speed flames and DDT in very rough-walled channels[J].Combustion and Flame, 2013, 160(1):204-211. doi: 10.1016/j.combustflame.2012.08.009 [2] Chen Xianfeng, Zhang Yin, Zhang Ying.Effect of CH4-Air ratios on gas explosion flame microstructure and propagation behaviors[J].Energies, 2012, 5(10):4132-4146. doi: 10.3390/en5104132 [3] Sklavounos S, Rigas F.Validation of turbulence models in heavy gas dispersion over obstacles[J].Journal of hazardous materials, 2004, 108(1):9-20. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=9c7b9ce98b1236bf2e9888f9b42111d6 [4] Arntzen B J.Modelling of turbulence and combustion for simulation of gas explosions in complex geometries[J].Journal of Loss Prevention in the Process Industries, 1998, 18(4/5/6):225-237. http://cn.bing.com/academic/profile?id=bd294e5486d8a62e3f3850d02e580e62&encoded=0&v=paper_preview&mkt=zh-cn [5] Sarli V D, Benedetto A D, Long E J, et al.Time-resolved particle image velocimetry of dynamic interactions between hydrogen-enriched methane/air premixed flames and toroidal vortex structures[J].International Journal of Hydrogen Energy, 2012, 37(21):16201-16213. doi: 10.1016/j.ijhydene.2012.08.061 [6] Masri A R, Ibrahim S S, Cadwallader B J.Measurements and large eddy simulation of propagating premixed flames[J].Experimental Thermal and Fluid Science, 2006, 30(7):687-702. doi: 10.1016/j.expthermflusci.2006.01.008 [7] Gubba S R, Ibrahim S S, Malalasekera W, et al.Measurements and LES calculations of turbulent premixed flame propagation past repeated obstacles[J].Combustion and Flame, 2011, 158(12):2465-2481. doi: 10.1016/j.combustflame.2011.05.008 [8] Kessler D A, Gamezo V N, Oran E S.Simulations of flame acceleration and deflagration-to-detonation transitions in methane-air systems[J].Combustion and Flame, 2010, 157(11):2063-2077. doi: 10.1016/j.combustflame.2010.04.011 [9] Gamezo V N, Ogawa T, Oran E S.Flame acceleration and DDT in channels with obstacles:Effect of obstacle spacing[J].Combustion and Flame, 2008, 155(1/2):302-315. http://www.sciencedirect.com/science/article/pii/S0010218008001934 [10] Ogawa T, Gamezo V N, Oran E S.Flame acceleration and transition to detonation in an array of square obstacles[J].Journal of Loss Prevention in the Process Industries, 2013, 26(2):355-362. doi: 10.1016/j.jlp.2011.12.009 [11] 任少峰.可燃性气体泄爆动力学机理研究[D].武汉: 武汉理工大学, 2012. http://cdmd.cnki.com.cn/Article/CDMD-10497-1012442364.htm [12] Lesieur M.Large-eddy simulations of turbulence[M].Cambridge:Cambridge University Press, 2005. [13] 肖华华.管道中氢/空气预混火焰传播动力学实验与数值模拟研究[D].合肥: 中国科学技术大学, 2013. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y2354259 期刊类型引用(10)
1. 梁冠军，刘晓东，尤超，张涛. 聚脲弹性体加固砌体墙抗爆试验安全风险评估. 佳木斯大学学报(自然科学版). 2023(03): 145-150 . 
百度学术2. 谭英华，王国梁，李阳，胡亚超，席丰. 脉冲荷载作用下承重砌体墙动力行为分析. 山东建筑大学学报. 2023(04): 10-18 . 百度学术3. 陈德，吴昊，方秦. 爆炸荷载作用下单向砌体填充墙动态响应计算方法. 建筑结构学报. 2023(10): 197-210 . 百度学术4. 陈德，吴昊，徐世林，韦建树. 单向砌体填充墙激波管试验和动力行为分析. 爆炸与冲击. 2023(08): 136-154 . 
本站查看5. 尚雨露，徐轩，张帝，杨军. CONWEP与流固耦合爆炸加载差异性及砌体墙动力响应特征. 兵工学报. 2023(12): 3897-3908 . 百度学术6. 许林峰，陈力，李展，岳承军. 聚脲加固砖填充墙抗爆性能的试验和分析方法研究. 爆炸与冲击. 2022(07): 126-137 . 本站查看7. 王丽琼，肖杰. 砌体填充墙防爆性能及临界破坏状态的数值模拟. 安全与环境学报. 2022(05): 2461-2468 . 百度学术8. 胡嘉辉，吴昊，方秦. 近区爆炸作用下砌体填充墙损伤破坏与动态响应的数值模拟. 振动与冲击. 2021(09): 1-11 . 百度学术9. 周清，齐麟，丁杰. 构造柱对承重砌体墙抗爆性能的影响. 工业建筑. 2021(03): 42-49+21 . 百度学术10. 彭培，李展，张亚栋，陈力，方秦. 燃气爆炸作用下蒸压加气混凝土砌体墙的加固性能. 爆炸与冲击. 2020(03): 110-123 . 本站查看其他类型引用(13)
-
下载:
下载:
百度学术
计量
- 文章访问数: 4932
- HTML全文浏览量: 1748
- PDF下载量: 428
- 被引次数: 23