GPU-accelerated smoothed particle hydrodynamics modeling of jet formation and penetration capability of shaped charges

The prediction of the penetration of three-dimensional (3D) shaped charge into steel plates is a challenging task. In this paper, the smoothed particle hydrodynamics (SPH) method is applied to simulate the jet formation generated by the shaped charge detonation and its damage to steel plates. The Jones–Wilkins–Lee (JWL) equation of state (EOS), Tillotson EOS, and elastic–perfectly plastic constitutive model were incorporated into SPH for the modeling of explosive detonation and dynamic behavior of metal material. The compute unified device architecture (CUDA) parallel programming interface has been employed in SPH to improve the computational efficiency of SPH. Firstly, the constitutive models and EOSs are validated by 3D TNT slab detonation and aluminum–aluminum (Al–Al) high velocity impact. Then the jet formation of the shaped charge detonation and its penetration into the steel plates are investigated using the graphics processing unit (GPU)-accelerated SPH methodology. The numerical results of these test cases are compared against the published experimental data or analytical result, which shows that the GPU-accelerated SPH methodology is capable of tackling the 3D shaped charge detonation and penetration involving millions of particles with high computational efficiency.     

Fracture behavior of explosively loaded spherical molded steel shells

Experimental and numerical works are made to study the fracture of explosively loaded spherical molded steel shells. The first series of experiments included three sawdust recovery shots to save fragments for examination. In this series, detonation was initiated from the center of the sphere. Results of the experiments show that two types of fractures are observed in spherical shells: radial and shear as in cylindrical shells. Spall fracture is also observed in spherical shells. To assist understanding of the experimental results, a computer simulation of expanding shells is performed to provide information on the stress, strain, strain rate and position of each element of the shell wall as a function of time after detonation. For t=7.5 μs after detonation, triaxial non-uniform strain prevails in the middle of the thickness of the wall. The peak of the stress equals to 6.5 GPa and exceeds the spall strength of carbon steel. In the second series of experiments, spall fracture is suppressed by displacing the point of detonation initiation from the center to periphery of HE charge.     

爆轰加载下金属柱壳膨胀破裂过程研究

采用前照明的高速分幅照相对不同壁厚的45钢柱壳在爆轰加载下的膨胀破裂过程进行了研究,实验成功拍摄到45钢柱壳前表面裂纹生成、扩展及产物泄漏的过程图像。研究表明:膨胀应变率在104s-1附近时,随着应变率的增加,45钢柱壳的各特征应变均有所增加,断裂模式由拉剪混合向单一剪切转变;45钢柱壳的膨胀破裂有一个过程,从柱壳的外壁面可以观察到裂纹生成至产物开始泄漏经历时间历程的长短与柱壳壁厚相关,柱壳壁越厚,从外壁出现裂纹到产物开始泄漏的时间间隔越长;由回收破片厚度测量推算出的膨胀断裂应变与动态实验过程中产物开始泄漏时的断裂应变基本一致。     

壳装炸药殉爆实验和数值模拟

进行了壳装固黑铝炸药殉爆实验,通过观察残留炸药、壳体和见证板变形,判断被发炸药的爆炸情况,得到了炸药临界殉爆距离。建立了壳装炸药殉爆实验计算模型,采用非线性有限元计算方法,对壳装固黑铝炸药殉爆实验进行了数值模拟。计算中采用预设壳体单元破片方法描述主发炸药壳体破片的形成和破片对被发炸药的撞击起爆,炸药临界殉爆距离的计算结果与实验结果基本一致。主要是主发炸药中部的壳体破片撞击到被发炸药,被发炸药起爆位置也在装药中部。炸药壳体厚度主要影响破片速度和质量、被发炸药的防护性能,从而影响炸药临界殉爆距离。 更多还原     

爆轰波在阻尼管道中声吸收的实验研究

实验旨在研究气相爆轰波在阻尼管道 (管壁上衬有吸收材料 )中传播时的衰减现象。先是在光滑管壁的管道中产生稳定的具有胞格结构的爆轰波 ,然后使其通过专门设计的管壁上衬有吸收材料 (金属丝网或不锈钢纤维 )的阻尼段。利用高速摄影、压力传感器和烟迹技术等手段 ,记录和测试了阻尼段对几种混合气体爆轰波的传播速度、爆压及胞格结构产生的影响。实验分别在方管和圆管中进行。发现在某些条件下爆轰波可以被衰减成强爆燃。     

爆轰波对碰驱动下金属圆管膨胀变形特性研究

采用高速分幅摄影和激光速度干涉仪（velocity interferometer system for any reflector, VISAR）联合测试技术,对HR-2钢管在爆轰波对碰驱动下的膨胀变形特性进行了研究。实验得到了圆管膨胀变形的过程图像以及圆管特定部位径向膨胀的位移-时间曲线和速度-时间曲线。采用LSDYNA-3D程序对HR-2钢管的膨胀变形过程进行了数值模拟,计算结果与实验结果吻合较好。     

滑移爆轰驱动钢管的层裂

将VG损伤模型推广到二维情况,考虑了最大主应力方向对损伤演化的影响,并使用显式断裂算法对20钢管在GI-920炸药滑移爆轰驱动下层裂的问题进行二维数值模拟。分析了一维内爆和滑移爆轰两种加载方式下作用于钢管外表面的压力及钢管内部受力状态的区别,考察了滑移爆轰加载方式下钢管外表面的受力随炸药厚度变化的规律,进而研究了钢管内损伤的分布和演化,以及裂纹的产生和扩展现象。计算得到的层裂片初始厚度随炸药厚度的变化规律与实验结果符合较好。     

激波风洞高低压段钢膜片破裂特性研究

激波风洞是用于高超声速飞行器气动外形设计和优化的常用地面试验装置,基于爆轰驱动技术,激波风洞能够在短时间(毫秒级)内产生高温、高压的驱动气体来模拟高超声速试验气流.主膜片位于激波风洞中的爆轰驱动段和激波管段之间,试验时膜片在爆轰脉冲压力下打开,膜片的打开状态和脱落情况对激波风洞气流品质有很大的影响. 同时,膜片也是形成激波的先决条件. 传统的风洞采用铝质膜片进行试验,在激波风洞中需要承压能力更强的膜片, 此时铝质膜片不再适用, 需要采用钢质膜片.因此, 对激波风洞中的钢膜片破裂特性进行研究很有必要.将数值计算结果与试验结果进行比较, 发现数值计算结果与试验结果吻合得比较理想,计算结果具有可靠性. 基于膜片的应力-应变模型, 建立了膜片打开的动力学模型,根据CJ爆轰理论, 采用有限元软件计算模拟了膜片破裂的过程,分析总结了膜片破裂的机制和力学特性规律.采用控制变量法对不同厚度和凹槽长度的膜片进行分析研究,得到了膜片破膜压力和有效破膜时间的变化规律. 在激波风洞试验中,根据膜片总破膜时间设计了适用于JF-12复现风洞的膜片参数.      

金属圆管在炸药滑移爆轰作用下的运动姿态

采用二维Ｌａｇｒａｎｇｅ数值摸拟方法，计算在管内装填岩石硝铰炸药滑移爆轰作用下，２０＃无缝钢管的运动过程，计算结果与用斜电阻丝法测量得到的结果基本相符。     

聚奥-9C装药的传爆管殉爆

为了研究冲击波作用下引信传爆装置的响应规律,进行了以主发炸药为RDX-8701、被发装置为聚奥-9C（JO-9C）装药的传爆管（含导爆药柱）的殉爆实验。通过观测残留传爆药、壳体和见证块变形,判断传爆管的爆炸程度,分析了殉爆过程中JO-9C爆轰波的成长历程及传播规律。采用AUTODYN软件建立了殉爆实验有限元模型,计算模型中主要考虑了主发炸药产生的爆炸冲击波对传爆管的冲击响应。基于流固耦合方法,通过调整距离模拟计算得到了传爆管的临界殉爆距离和殉爆安全距离。结果表明,传爆管上端的侧角受到爆炸冲击后产生的爆轰波沿斜下方传播,使传爆药柱完全爆轰,并引起导爆药柱发生殉爆;数值模拟结果显示,JO-9C装药的传爆管临界殉爆距离为5.7 mm,殉爆安全距离为8.8 mm。     

Critical energy for direct initiation of detonation induced by laser ablation

This study describes experimental work examining the critical energy for direct initiation of detonation by laser ablation in a stoichiometric acetylene–oxygen mixture. The amount of input energy, the target material, and the surface roughness of the target were varied to study their effects on shock wave generation. Aluminum and stainless steel were used as target materials. The propagating shock wave induced by laser ablation was observed using high-speed shadow imaging. The critical energy for direct initiation of detonation was calculated using the strong blast wave theory. The critical input energy for aluminum was found to be lower than that for stainless steel. Because the thermodynamic critical temperature of aluminum is lower than that of stainless steel, aluminum caused a phase explosion more easily than stainless steel, thus resulting in direct initiation of detonation with a lower amount of input energy. The effects of surface roughness on critical input energy and shock wave generation were negligibly small. The critical initiation energy was estimated to be \(10.3 \pm 0.2\) mJ, which is in agreement with the experimental data obtained in previous work. The estimated critical initiation energy was independent of the target material. However, other predictions of the critical initiation energy by using the cell size overestimated this value because of the scatter in cell size data of an unstable cellular structure. Furthermore, interaction between plasma plumes formed by laser ablation and those formed by breakdown near the target surface might have contributed to requiring a lower amount of energy for initiating detonation.     

密闭空间内爆炸载荷抑制效应实验研究

炸药在密闭空间内爆炸时的爆炸载荷与在敞开环境中有很大差异,在密闭空间内,TNT炸药的爆炸产物能够与周围空气充分混合并发生燃烧反应进而释放额外的能量,使密闭空间内的反射冲击波及准静态压力均明显增加。为探究不同气体环境对密闭空间内爆炸载荷的抑制效应,开展了3种不同药量的TNT分别在空气、水雾和氮气环境密闭空间内的爆炸实验研究,通过理论计算和实验对比分析了密闭空间内的爆炸载荷压力、温度及受载钢板试件响应特性。结果表明,水雾和氮气均能有效降低空间内的准静态压力和温度,对准静态压力的平均降幅分别为36.0%和51.7%,对温度的平均降幅分别为42.6%和40.3%;在水雾和氮气环境中的爆炸载荷作用下,钢板试件动态响应较空气环境中显著降低,其中160 g药量下,水雾和氮气环境中钢板试件的最终变形分别减少了15.9%和23.5%,氮气的减弱效果优于水雾;水雾和氮气均能及时有效地抑制封闭空间内的爆炸载荷,降低结构的损伤程度。     

水中钢板爆破水介质对装药量的影响

为分析水介质对水中钢板爆破装药量的影响,理论推导了钢板背衬水介质条件下,钢板爆破最小装药量与空气中钢板爆破最小装药量的倍数关系:最小倍数关系为3.76。数值计算了钢板背衬空气介质和背衬水介质情形下钢板爆破的最小装药量,其倍数关系在3.5左右,与理论结果相近,表明钢板背衬水介质或空气介质是决定装药量大小的关键因素。数值计算了装药在空气介质中及水介质中钢板爆破爆轰产物对钢板的冲量大小,结果接近,表明水介质对炸药爆轰产物的约束作用是影响水中钢板装药量的次要因素。     

An experimental study on DDT distance of JP-8 mixtures

Previous research into detonation physics has mostly utilized gaseous fuels such as hydrogen, acetylene, ethylene, and propane. If these fuels were to be used for a pulse detonation engine, they have to be stored under high pressure in steel containers which increase weight safety risks. In order to increase energy density of fuel, liquid fuel was chosen. Tests were conducted on detonation initiation of JP-8/oxygen mixtures at different initial temperatures and equivalence ratios. These tests found a reduction in the rich limit with increasing initial temperature, and the minimum deflagration-to-detonation run-up distance was approximately 200 mm, which was similar to propane/oxygen mixture results. A rapid increase in deflagration-to-detonation run-up distance was observed at equivalence ratios close to the lean and rich limits. Experiments of JP-8/oxygen and propane/oxygen mixtures with nitrogen dilution were also conducted. As the nitrogen/oxygen ratio increased, the lean and rich limits decreased while the detonation wave could not be successfully initiated as the nitrogen-to-oxygen ratio was greater than 0.4.     

内部爆炸作用下钢筒变形过程的电探针测量技术

为了研究圆柱形爆炸容器在炸药爆炸作用下的变形过程,设计了电探针测量内部爆炸作用下钢筒变形的方法。采用数值模拟方法进行预估,在钢筒的中心进行了120g TNT和180g TNT当量球形装药下的爆炸加载实验,获得了爆炸不同时刻钢筒径向位移随时间的变化关系,电探针测量钢筒最大变形与实验后钢筒变形测量结果较好吻合。     

水及盛水容器对近距离爆炸载荷影响的实验研究

为研究用水包围炸药的方式对爆炸载荷的作用机理,在两端开口的钢筒内进行了水直接包覆炸药的爆炸实验,利用光纤位移干涉仪获取了钢筒外壁的径向速度和动态变形。结果表明:不同于无水爆炸,爆轰产物通过水的“裂缝”在空气中形成冲击波,造成该冲击波出现时间更晚、强度更低、持续时间更长,并要求相应的数值模拟采用二维以上的计算模型;盛水结构的材料密度越低、厚度越小对爆炸载荷的影响越小。     

爆轰波碰撞的聚能效应

利用两高爆速导爆索对称布置于药卷两侧,起爆后炸药爆轰波在对称线处汇聚碰撞,当碰撞角度达到一定值时,发生马赫反射,使爆轰压力成倍增加,形成高压、高能量密度区域的聚能效应。本文在爆轰波传播碰撞理论的基础上,通过炸药做功能力和猛度试验验证爆轰波碰撞的聚能效果。做功能力试验结果表明爆轰波碰撞能够增加炸药能量利用率;猛度试验结果表明采用对称起爆技术下的爆轰波碰撞能够改变爆轰波在特定方向上的扩散作用。试验结果与爆轰波入射角的几何关系表明,当高爆速起爆药条与主装药爆速比例在1.15倍以上时,爆轰波碰撞能够达到一定的聚能效果。     

破片撞击起爆柱面带壳装药的临界速度修正判据

为获得适用于柱面带壳装药的冲击起爆修正判据，以Picatinny工程判据为基础加入修正项进行修正。采用AUTODYN-3D软件对破片撞击柱面带钢壳的B炸药进行数值计算，获得了破片入射角、装药曲率半径对炸药临界起爆速度的影响规律；通过拟合得到修正项表达式，建立了考虑破片入射角、柱壳装药形状函数的炸药起爆临界速度修正判据。判据计算值与实验数据和数值计算值吻合较好，该判据能较好的预测柱形带壳装药的冲击起爆条件。     

一维与二维爆轰传播的时空关联特性数值研究

一维爆轰传播的理论完备、计算准确, 二维斜爆轰传播由于壁面与黏性效应, 大尺度、高精度预测还有一定难度. 利用Euler方程和H₂-Air基元反应模型, 对二维有限长楔面诱导的斜爆轰和活塞驱动一维非定常正爆轰进行计算比较研究, 从时空两个维度方面, 分析了两者在起爆过程、稀疏波传播、爆轰波面演化中的关联特性. 研究结果表明: 在过驱动度相同的条件下, 经过时空变换的活塞驱动一维爆轰传播与二维驻定斜爆轰在起爆区波系结构、波面演化特征和主要参数分布规律方面无论定性或者定量对比均符合较好, 所以, 一维非定常爆轰和二维驻定斜爆轰具有时空相关性. 两者的差异主要体现在过驱动斜爆轰受稀疏波影响过渡到近Chapman-Jouguet (C-J)爆轰状态所需的弛豫时间不同, 原因可能是起源于活塞和壁面稀疏波强度的差异. 本文提出的一维与二维爆轰传播的时空关联方法不仅有助于认知斜爆轰起爆、过驱爆轰产生、胞格爆轰演化的三阶段规律, 还可以对比揭示壁面、边界层和黏性效应的影响, 应用在斜爆轰发动机燃烧室设计中能够有效节约计算时间和成本, 并降低复杂度.