爆炸冲击波在防护挡墙下的削波作用研究.docx

《爆炸冲击波在防护挡墙下的削波作用研究.docx》由会员分享，可在线阅读，更多相关《爆炸冲击波在防护挡墙下的削波作用研究.docx（4页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、爆炸冲击波在防护挡墙下的削波作用研究 【摘 要】采纳多流体网格法，基于JWL状态方程，用面对对象的C+语言自行开发编制了NMMMIC通用多物质二维流体弹塑性程序，对有和无防护档墙的爆炸空气冲击波进行了数值模拟，并用VISC 2D可视化软件对计算结果进行动画演示，从中不仅可以看到爆炸过程的整体概貌，还可以清楚地显示微小过程。计算结果基本符合冲击波饶过挡墙的物理现象和规律，这说明该物理模型和数值算法较为合理，可用于对可能发生的爆炸事故进行风险评估。 【关键词】冲击波；多流体网格法；数值模拟 0.引言 在现代生活中，爆炸灾难的发生与日俱增，造成建筑物和设备的破坏及人员的伤亡。大多数的工业灾难事故是由

2、爆炸引起的。由爆炸产生的空气冲击波遇到目标，如建筑物、军事设施等会发生反射和绕射现象，探讨这些现象对于摧毁敌方的军事目标和加强我方的防卫工事，以及对危急品生产车间的防护都有很大的好用价值。本文采纳多流体的Euler型算法，运用自行编制的NM-MMIC通用程序，模拟了在爆点四周有无防护挡墙的二维爆炸场，并通过考虑防护挡墙的形态对其后爆炸场的影响，给出了一些对工程设计有参考价值的结果。 1.计算方法 1.1总体方程组 爆炸冲击波的形成过程可用二维不定常弹塑性流体力学方程组描述1。在冲击波的数值模拟中，涉及到三种介质：炸药、土和空气。炸药采纳TNT， 爆轰产物的初始条件采纳瞬时爆轰模型，状态方程分别

3、为：对于空气，近似采纳志向气体状态方程，对炸药采纳JWL状态方程，土介质所采纳的状态方程是： p=8.62610ee 式中，P为压力， pv为体积应变。 图1带防护墙的初始图 1.2数值计算方法 计算域取150150m，网格划分为600600个。网格步长取0.25 。边界条件为：左边界取对称边界，其他三边取连续边界条件，计算采纳显式差分格式。应力，密度，比内能，速度，应力偏量等均定义在网格中心，人工粘性定义在网格边界中点2。对于含两种物质的混合格，运用Youngs界面处理技术，当混合网格中有三种物质时，采纳简洁线界面技术。图1为正梯形挡墙的计算模型。 2.数值模拟结果 本文探讨了无挡墙、矩形挡

4、墙、正梯形挡墙、左斜右直挡墙等不同形态和位置的防护挡墙对爆炸空气冲击波的影响，考察了圆柱形区域的高压气体对地面和防护挡墙的冲击波流场的发展规律3。图2为正梯形挡墙的爆炸过程。 图2 正梯形挡墙的爆炸过程 通过计算可以得出如下结论： 靠近爆心一侧的防护挡墙形态对远场压力的影响较大，梯形挡墙要好于矩形挡墙；而远离爆心一侧的防护挡墙形态对远场压力的影响不明显，梯形挡墙与左斜右直的挡墙对远场的效果差别不大。 无防护挡墙的超压要比有防护挡墙时大得多，这表明防护挡墙在减弱远场压力方面还是起到了肯定作用的。 在挡墙后方高度的2倍距离处的峰值压力比无挡墙时要大得多，这是因为在此处形成马赫反射的原因。这一点在设

5、计防护挡墙时应特殊留意。 3.结论 本文针对爆炸冲击波问题，采纳多流体网格法，探讨了防护挡墙对爆炸效应的影响。针对不同位置、不同形态的挡墙，给出了爆炸场的各个物理量随时间的改变关系。在此基础上重点探讨了远场压力对防护挡墙的依靠关系，为工程设计供应了非常有用的数值分析手段。 【参考文献】 1宁建国，张文耀，郑祖国，李伟.爆炸与冲击数值模拟的可视化.兵工学报，第21卷增刊，2000，8：89-91. 2宁建国，王仲琦，赵衡阳等.冲击波绕流的数值模拟探讨J.北京理工高校学报，11019，19. 3郝莉.多物质二维流体弹塑性程序的理论基础及其应用探讨，北京理工高校博士学位论文，2003，8. 第4页 共4页第 4 页 共 4 页第 4 页 共 4 页第 4 页 共 4 页第 4 页 共 4 页第 4 页 共 4 页第 4 页 共 4 页第 4 页 共 4 页第 4 页 共 4 页第 4 页 共 4 页第 4 页 共 4 页