Simulational experiment on compression and fracture of propellant charge based on the discrete element method
-
摘要: 为了揭示发射装药破碎引起的膛炸现象，急需进行相应装药结构下发射装药挤压破碎数值模拟研究。以硝胺花边十九孔发射药为研究对象，基于离散单元法建立了发射装药挤压破碎模拟系统，同时进行了发射装药动态挤压破碎实验，通过数值模拟与实验获得了不同冲击载荷下的破碎发射装药和挤压应力；分别对获得的破碎发射装药进行了密闭爆发器数值模拟和实验。结果表明：模拟与实验获得的发射装药挤压应力时间历程、密闭爆发器压力时间曲线和起始动态活度比的一致性较好，实验验证了发射装药挤压破碎模拟系统的有效性及合理性。该模拟系统具有重大工程应用价值，为高能发射装药冲击破碎过程和发射装药发射安全性研究奠定了基础。Abstract: In order to reveal the mechanisms of gun breech-blow phenomenon caused by the fracture of propellant charge, it is urgent to carry out the simulation research on the compression and fracture of propellant charge under the corresponding charge structure. Through the analysis of the mechanical environment in the gun bore and the fracture progress of propellant charge, the discrete element method was employed to simulate the compression and fracture of propellant charge. The lace 19-hole propellant for the large caliber artillery was taken as the research object, a simulation system of compression and fracture of propellant charge was constructed using the EDEM software. And the Hertz-Mindlin contact model parameters were determined by using the drop hammer impact test of the single propellant at low temperature (–40 ℃). Then the compression andfracture simulation of the propellant charge was verified through the dynamic compression and fracture test of the propellant charge at low temperature (–40 ℃). Under the same impact load, the fracture of propellant charges and the compression stress-time curves of propellant charge were achieved by test and simulation, respectively. Using the obtained fracture propellant charge, the closed bomb simulation and test were carried out respectively. Among them, a combustion function based on the discrete element method was used to represent the gas generation law of the simulation fracture propellant charge. Finally, the initial dynamic vivacity ratio of the fracture propellant charge was processed according to the pressure-time curve. The researchresults show that the time histories of compression stress of propellant charge, the closed bomb pressure-time curves, and the initial dynamic vivacity ratios obtained by simulation and test are in good agreement with each other, indicating the designed simulation system is effective and reasonable. The research method has great engineering application value, which lays a foundation for the study of the impact fracture process of high-energy propellant charge and the launch safety of propellant charge.
-
砌体填充墙被广泛地应用于建筑物的外墙或钢筋混凝土框架内。但是, 由于它抗爆能力弱, 容易受到爆炸荷载的破坏, 且墙体的破坏与否对主体结构承受的爆炸荷载有一定影响[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时, 墙体发生飞散。
-
图 6 发射装药动态挤压破碎实验装置及原理图
Figure 6. Experiment device and schematic diagram of the dynamic compression and fracture of propellant charge
图 8 低温下实验获得的破碎发射装药
Figure 8. Fracture propellant charge obtained by experiments at low temperature
图 9 采用离散单元法模型获得的挤压破碎模拟结果
Figure 9. Compression and fracture simulation results with the discrete element model
图 10 采用黏结键模型得到的挤压破碎模拟结果
Figure 10. Compression and fracture simulation results with the bond model
图 11 不同燃烧室压力下模拟与实验挤压应力对比曲线
Figure 11. Comparison of compression stress curves between simulation and experiment under different chamber pressures
图 12 密闭爆发器数值模拟与实验压力时间对比曲线
Figure 12. Comparison of simulated and experiment pressure-time curves in the closed bomb
表 1 数值模拟参数与结果
Table 1. Model parameters and simulation results
序号 法向黏结刚度/（GN·m−3） 切向黏结刚度/（GN·m−3） 法向临界应力/MPa 切向临界应力/MPa 最大应力/MPa 1 449.55 134.87 120 36.0 87.31 2 449.55 134.87 125 37.5 90.95 3 449.55 134.87 130 39.0 95.31 4 449.55 134.87 135 40.5 98.24 5 449.55 134.87 140 42.0 102.67 6 449.55 134.87 150 45.0 110.14 
下载: 导出CSV
表 2 实验条件和结果
Table 2. Experiment conditions and results
序号 燃烧室最大压力/MPa 挤压应力峰值/MPa 1 37.46 3.96 2 40.43 11.06 3 44.81 15.35
下载: 导出CSV
表 3 数值模拟结果
Table 3. Simulation results
序号 挤压应力峰值/MPa 黏结键连接个数 1 3.92 18 498 2 9.14 15 190 3 13.34 12 233
下载: 导出CSV
表 4 数值模拟与实验起始动态活度比对比
Table 4. Comparision of initial dynamic vivacity ratios in simulation and experiment
序号 数值模拟的起始动态活度比 实验的起始动态活度比 误差/% 1 1.144 1.189 3.78 2 1.570 1.628 3.56 3 1.942 2.067 6.05
下载: 导出CSV
-
[1] 芮筱亭, 贠来峰, 王国平, 等. 弹药发射安全性导论[M]. 北京: 国防工业出版社, 2009: 1–5.RUI X T, YUN L F, WANG G P, et al. Direction to lanch safety of ammunition [M]. Beijing: National Defense Industry Press, 2009: 1–5. [2] 芮筱亭, 冯宾宾, 王燕, 等. 发射装药发射安全性评定方法研究 [J]. 兵工学报, 2015, 36(1): 1–11. DOI: 10.3969/j.issn.1000-1093.2015.01.001.RUI X T, FENG B B, WANG Y, et al. Research on evaluation method for launch safety of propellant charge [J]. Acta Armamentarii, 2015, 36(1): 1–11. DOI: 10.3969/j.issn.1000-1093.2015.01.001. [3] GAZONAS G A, FORD J C. Uniaxial compression testing of M30 and JA2 gun propellant using a statistical design strategy [J]. Experimental Mechanics, 1992, 32(2): 154–162. DOI: 10.1007/BF02324727. [4] 陈言坤, 甄建伟, 武慧恩, 等. 粒状发射药动态破碎研究进展 [J]. 爆破器材, 2014, 43(1): 43–48. DOI: 10.3969/j.issn.1001-8352.2014.01.010.CHEN Y K, ZHEN J W, WU H E, et al. Research progress of dynamic fracture of granular propellant [J]. Explosive Materials, 2014, 43(1): 43–48. DOI: 10.3969/j.issn.1001-8352.2014.01.010. [5] CUNDALL P A. A computer model for simulating progressive, large-scale movements in block rock systems [C] // Symposium of International Society of Rock Mechanics. Nancy, France, 1971. [6] JIANG S P, RUI X T, HONG J, et al. Numerical simulation of impact breakage of gun propellant charge [J]. Granular Matter, 2011, 13(5): 611–622. DOI: 10.1007/s10035-011-0276-1. [7] 王燕, 芮筱亭, 宋振东, 等. 初始堆积对发射装药底部挤压应力的影响 [J]. 爆炸与冲击, 2014, 34(5): 560–566. DOI: 10.11883/1001-1455(2014)05-0560-07.WANG Y, RUI X T, SONG Z D, et al. Effect of original packing on compression stress at the bottom of propellant bed [J]. Explosion and Shock Waves, 2014, 34(5): 560–566. DOI: 10.11883/1001-1455(2014)05-0560-07. [8] EDEM tutorial [M]. Edinburgh: DEM Solutions, 2011. [9] POTYONDY D O, CUNDALL P A. A bonded-particle model for rock [J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(8): 1329–1364. DOI: 10.1016/j.ijrmms.2004.09.011. [10] QUIST J. Cone crusher modelling and simulation [D]. Göteborg: Chalmers University of Technology, 2012: 1–45. [11] 夏露. 高能材料结构和性能的分子动力学模拟[D]. 苏州: 苏州大学, 2008: 46.XIA L. Molecular dynamics simulations of the structures and properties of highly energetic materials [D]. Suzhou: Suzhou University, 2008: 46. [12] 徐浩, 芮筱亭, 郁兆华, 等. 发射装药挤压应力与破碎规律研究 [J]. 火炸药学报, 2012, 35(4): 61–68. DOI: 10.3969/j.issn.1007-7812.2012.04.016.XU H, RUI X T, YU Z H, et al. Study on extrusion stress and fracture law of propellant charge [J]. Chinese Journal of Explosives and Propellants, 2012, 35(4): 61–68. DOI: 10.3969/j.issn.1007-7812.2012.04.016. [13] ZHANG R H, RUI X T, WANG Y, et al. Study on the change of gas generation law caused by fracture of propellant charge [J]. Journal of Energetic Materials, 2018, 36(4): 454–467. DOI: 10.1080/07370652.2018.1486895. 期刊类型引用(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)
-
下载:
百度学术
计量
- 文章访问数: 672
- HTML全文浏览量: 337
- PDF下载量: 83
- 被引次数: 23