平面铅飞片对碰的SPH模拟的研究.pdfVIP

平面铅飞片对碰的 SPH 模拟研究 袁帅，胡海波，张崇玉，李平 （中国工程物理研究院流体物理研究所冲击波物理与爆轰物理实验室，四川绵阳 621900 ） 摘要：简要介绍了 SPH 方法及其应用领域，分析了爆轰波对碰驱动飞层的研究背景和 应用。用 SPH 方法模拟研究了两点起爆圆柱炸药驱动组合飞层（AL+Pb ）的实验模型，得 到了爆轰波对碰区铅飞层类似射流状超前凸起的计算图像，对碰凸起呈现与实验结果相同 的多分层结构：凸起头部出现散碎特征，根部相对密实。结果表明：SPH 方法对模拟研究 铅这类低熔点低强度金属在对碰区出现散碎等破坏现象具有良好前景。 关键词：SPH 组合飞层 对碰凸起 Pb 1 引言： 光滑粒子流体动力学方法SPH (Smoothed Particle Hydrodynamics)是近几十年来逐步发 展起来的一种无网格方法[1，2] ，该方法的基本思想是将连续体用相互作用的质点组来描述， 各个物质点上承载各种物理量，包括质量、速度等，通过求解质点组的动力学方程和跟踪每 个质点的运动轨迹，求得整个系统的力学行为。由于质点之间不存在网格关系，因此它可避 免极度大变形时网格扭曲而造成的精度破坏等问题，并且也能较为方便的处理不同介质的交 界面。SPH的优点还在于它是一种纯Lagrange方法，能避免Euler描述中欧拉网格与材料的界 面问题，因此特别适合于求解高速碰撞等动态大变形问题。SPH被称为下一代数模技术,它 是一种无网格质点方法./克服了高维拉氏方法中网格缠绕的问题,和欧拉法随着维数、网格数 [3] 增加而变得吃力的问题。它基于特殊的算法，对应用的限制大为减少，可用于高速侵彻 、 [4] [5] [6] 爆炸冲击 、微喷射流 、两相流 等等数值模拟。 在爆轰驱动研究中，常常会遇到两个爆轰波相互碰撞驱动壳体的情况。两个爆轰波碰撞 后在很窄的区域内产生较高的压力，约为炸药CJ压力的 2.3~2.4 倍[7]，将该区域称为爆轰波 对碰区。爆轰波对碰导致对碰区的压力温度急剧升高，可能会引起材料特殊的表面喷射、变 形失稳、断裂破坏等动力学响应。爆轰波对碰驱动材料的动载行为表现，由于其独特的力学 行为特征以及重要的工程应用背景和理论研究价值，近年来逐渐一起研究人员的关注。 Sing.M[8] [9] [10，11] 等、陈军 等、张崇玉 等先后对爆轰波对碰驱动金属圆管的膨胀特性进行 了实验研究。实验结果表明，在爆轰波对碰驱动金属圆管对碰区出现超前现象，且该部位壳 体较邻近部位易发生提前断裂。Zhiembetov A. K[12]等用铅、铟、锡、铁等材料开展对碰实 验研究发现，采用铅、铟、锡这类低强度、低熔点金属进行实验时，观察到对碰区突起部位 出现多孔断裂、崩溃以及材料的破碎、雾化等现象。ИвановАГ[13]等采用X光照相技术对R/T 不稳定性炸药对碰驱动平面铝飞层对碰区的行为特性进行了实验研究，获得了铝飞层对碰区 凸起的形貌和随时间变化的图像，并认为借助合适的材料模型和不断发展的计算机技术可以 对该类现象进行量化描述。张崇玉[14]等完成了铅飞层的对碰实验，采用光学分幅和脉冲X光 照相两种测试技术对平面铅飞层在爆轰波对碰驱动下的动载行为进行了观测，认为：铅飞层 对碰区出现多分层区、散碎雾化的凸起结构特征，除了材料强度和冲击熔化因素外，可能还 与对碰爆轰波剖面结构细节有直接的对应关系。 本研究基于张崇玉的平面铅飞层在爆轰波对碰驱动实验模型，用光滑粒子方法（SPH） 对张崇玉的爆轰波对碰驱动平面铅飞层实验进行数值模拟计算，得到了爆轰波对碰区铅飞层 类似射流状超前凸起的计算图像，对碰凸起呈现与实验结果相同的多分层结构，且头部出现 散碎特征，根部相对密实。结果表明：SPH 方法对模拟研究铅这类低熔点低强度金属在对 碰区出现散碎等破坏现象具有良好前景。 2 实验模型 实验装置如图 1 所示，铅飞层和 Ly12 铝飞层各厚 2.0 mm，炸药为 PETN ，尺寸为φ32 mm×2