5.爆破漏斗成型过程仿真模拟.pdfVIP

爆破漏斗成型过程仿真模拟 1, 2 3* 4 1 牟宏霖 ，王秋红 ，王永强 ，卢兆明 （1. 香港城市大学 土木建筑工程系，香港；2. 中国科学技术大学 火灾科学国家重点实验室，安徽 合肥 230027 ；3. 西安科技大学 能源学院 陕西 西安 710054 ；4. 西南科技大学 环境与资源学院，四川 绵阳 621010 ） 摘 要：爆破漏斗作为一种爆炸破坏的基本形式，其成型机理在爆破工程研究领域 占有重要地位。同时，理解爆破漏斗的成型过程对爆破设计和安全防护工作具有实 际意义。然而在传统的爆破测试技术中，部分操作程序比较复杂且伴随着威胁测试 员安全的危险，这使得爆破过程中一些参数不易记录或耗费过高。计算机数值模拟 技术可以有效解决这一问题，本文利用ANSYS LS-DYNA 软件对爆破漏斗成型过 程进行了仿真模拟。以一个新的方法表现了爆破漏斗的爆炸破坏过程，分析了冲击 波的压力传播和应力应变等参数，继而验证了相关爆破现象及理论。该方法可用于 指导爆破工程中的生产及实践，并为爆破设计和安全防护等工作提供更好的数据支 持。 关键词：爆破漏斗；成型过程；有限元分析；数值模拟；DYNA 0 引言 爆破是利用炸药爆炸时所释放的能量来破坏某种介质或使介质变形，从而达到一定 [1] 的工程或工艺目的的技术。爆破漏斗是爆破破坏的基本形式之一 。在只有一个水平临空面 的半无限岩土介质中，当埋置于一定深度的药包在介质中爆炸时，除了将岩土介质破碎以外， 还会将一部分破碎了的物质抛掷，在地表形成一个漏斗形的凹坑，这个凹坑就是爆破漏斗。 由于爆破过程的复杂性，其涉及到爆炸力学、爆轰学、气体动力学、土力学、岩石动力 学、固体力学等多种学科。因此，尽管当前爆破技术已有巨大发展，但从理论上对爆破漏斗 机理的研究还远没有完善。此外，炸药的爆炸是一个剧烈的过程，在通常的爆炸条件下，爆 轰波每单位平方厘米上产生的功率可达数百万千瓦。同时爆破的过程受太多复杂因素影响， 精确的理论计算和运用实验手段研究爆破漏斗机理都存在很大困难。 随着计算机科学和数值模拟计算方法的迅速发展，现已可实现在计算机上运用ANSYS LS-DYNA 软件通过直接在岩土体数值模型中埋设炸药材料的方法，来研究爆破漏斗的成型 [2] 过程 。除此之外，在数值模拟中还可选取不同的药包埋藏深度和调整其他参数以做出最优 化的应用理论研究。本文运用有限元分析方法对爆破漏斗的成型过程进行了数值模拟；同时 对其中的破坏过程现象，冲击波的压力传播和应力应变等参数进行了分析与验证。所得到的 基础数据和研究成果可用于指导爆破工程中的生产与实践过程，并对爆破设计和安全防护[3] 等工作提供新的科学方法。 1 爆破数值模型前处理 1.1 数值模拟技术路线 本文应用ANSYS LS-DYNA 软件程序对爆破漏斗进行数值模拟，其中采用了有限元算 [4] 法 。有限元法是先将连续的求解域分解成有限单元，组成离散化模型，然后求其近似的数 值解。有限元包括结构有限元和动力有限元，动力有限元适合于计算边界形状复杂或者包含 - 1 - 物质界面的强动载问题计算。其在冲击问题的模拟计算方面得到了迅速发展和广泛应用。 本文中爆破漏斗数值模拟实验的建模分析是：数值模型由炸药、土壤层和空气三个部分 [5] 组成，其中炸药、土壤和空气三种材料采用欧拉网格建模，单元使用共用节点算法 ，空气 和土壤材料间采用定义面面接触。由于是条形装药，在中心线起爆条件下，不考虑端部效应 时，可以将模型简化成平面对称问题。为了优化计算时间，本算例采用单层实体网格建模。 数值模型采用cm-g-μs 单位制，整个模拟技术路线如图 1 所示。 模型中的实体单元选用SOLID164，它是8 节点六面体单元，不含零能模式并且不需沙 漏