矩形坑道内爆炸冲击波压力仿真计算及试验验证

针对坑道试验爆轰波传播进行研究,通过数值模拟手段获取不同爆心距处坑道侧壁超压时程曲线,并通过试验手段,验证了数值模拟模型正确性。对比数值模拟结果和试验数据,发现坑道试验将混合装药等效TNT质量后,数值模拟超压峰值与实测误差范围为8.49%～25.49%,在距爆心15 m外数值模拟误差较小;等截面坑道内,爆轰波在近爆心处峰值超压相差较为明显,反映近爆心爆轰波较为紊乱。随着距爆心越远,等截面处峰值压力趋于一致,表明爆轰波呈稳定传播状态。研究结果表明,坑道内爆试验将混合装药等效为TNT当量,采用数值模拟研究坑道不同位置峰值超压方法可行性,为坑道试验超压测试传感器量程、规格选择提供理论指导,为坑道内侵爆类弹药爆炸毁伤效应评估等提供支撑。     

采用不同化学反应源项处理方法的胞格爆轰数值研究

采用一种改进的化学非平衡流解耦方法对Euler反应流方程进行解耦处理,对流项采用五阶WENO格式进行离散,化学反应源项分别采用梯形公式和α-QSS拟稳态逼近两种方法处理,时间步进采用二阶精度的Runge-Kutta方法,对H2/O2/Ar预混气的胞格爆轰进行了二维数值模拟。两种处理方法都得到了合理的三波点结构,计算结果接近,但α-QSS拟稳态逼近处理方法的化学反应较完全,爆轰波传播速度较高;揭示了爆轰波结构的发展变化过程,横波与横波、横波与壁面碰撞引起的二次起爆对爆轰波的稳定发展至关重要;分析了数值胞格结构,胞格长宽比与参考文献基本一致,验证了本计算方法的有效性。     

多点偏心起爆对破片速度增益的影响

多点偏心起爆后,爆轰波相互作用在定向区域内产生了马赫波超压,使得该区域破片速度产生增益.利用爆轰波反射理论,建立多点偏心起爆后爆轰波相互作用特性的理论模型,并确定定向区内破片速度的增益特性.3种不同起爆方式的仿真计算和静爆试验结果对比表明,定向区内不同起爆方式得到的破片增益程度较为吻合.     

Q235钢与304不锈钢多层爆炸焊接的数值模拟研究

利用SPH-FEM耦合算法及SPH算法对Q235钢与304不锈钢多层爆炸焊接实验分别进行了三维与二维数值模拟,旨在研究多层爆炸焊接过程中动态参数分布情况及结合质量。理论计算了三维模拟中复板的碰撞压力峰值并与模拟结果进行了对比,在此基础上利用振动波及爆轰波的传播规律解释了复板压力分布特性。通过二维模拟获得了基复板波状结合界面,与实验所得波形较为吻合,验证了多层爆炸焊接实验中基复板的结合质量。模拟结果显示了相较于单层爆炸焊接,多层爆炸焊接理论上能够节省近68%的炸药量,有效节约了生产成本。     

平面飞片动态高压装置的一维差分计算

本文讨论了炸药接触爆轰推动金属飞片,与多层材料组成的靶板高速碰撞,获得高压的一维差分数值计算方法,并对中国科学院物理所设计的爆炸法合成金刚石装置作了细致计算。此外,本文还给出了用D=C_0+S_0u及Gruneison状态方程导出的冲击波压力及冲击波恒压的时间宽度等计算公式,它们被用来直接计算平面飞片增压装置中的有关参量。     

爆炸激波管管口稀疏波对试验段的影响

数值模拟了爆炸激波管不同隔离段长度时管口稀疏波对试验段超压的影响。为了准确高效地模拟试验段入口的超压曲线,采用了一种将一维球对称程序和三维程序相结合的计算方法,并在一维计算中利用爆炸相似律,采用小当量爆炸来模拟实际超压波形。计算结果表明,隔离段长度L的变化不影响超压峰值;L小于等于20m时,稀疏波的影响使得试验段超压的作用时间、比冲量减小;L大于等于30m时,稀疏波对试验段超压无影响。     

氧气鱼雷供氧系统中扰动波的建模与仿真

鱼雷发射时,对高压氧气舱中的氧气起开关作用的启动阀快速开启,从而在启动阀的两端产生激波和膨胀波等两种扰动波.考虑了高压氧气的分子间相互作用力及分子所占体积的影响,采用范德瓦耳斯公式来代替完全气体的状态方程,同时建立了扰动波的数学模型,并推导了气体在接触间断面上的速度计算模型,仿真计算得出激波波阵面上的温度、由膨胀波所引起的低温及其它物理量的数值,为分析氧气鱼雷使用高压氧气的安全性提供了参考数据.     

聚能射流形成过程的理论建模与分析

分析了聚能射流的形成过程,并对其中的各阶段进行了详细建模。在模型中考虑了炸药爆轰、金属的驱动、药型罩压垮以及射流和杵体的形成过程。采用该模型对某一聚能装药结构进行了计算,计算结果表明:药型罩顶部和底部微元的压垮速度较小,在射流头部形成反向速度梯度,与试验数据吻合较好。该模型对于多级侵彻战斗部的工程设计与侵彻参数的计算具有一定的参考价值。     

爆轰波斜冲击作用下破片飞散特性研究

为研究非对称结构战斗部的破片飞散特性,利用斜激波理论对爆轰波作用于壳体表面的过程进行研究,并利用自由面速度倍增定律对波在自由面反射后质点速度的计算进行简化,得到了破片飞散角的计算模型。利用D型战斗部试验数据对计算模型进行验证,结果表明,斜激波理论计算得到的破片飞散与试验结果吻合很好;当入射角较小时,壳体飞散角与入射角成线性关系。     

爆炸合成金刚石中温度的计算

本文讨论材料在强冲击波作用下,状态参数温度的计算。在数十万大气压的冲击波作用下,固体材料的强度及电子对比热的贡献都可忽略。这样,可采用流体力学模型,用压力代替应力,材料的本构方程即为状态方程,本文取格留乃森状态方程形式。下面结合爆炸合成金刚石课题,讨论石墨、金刚石在强击波作用下,对其温度的数值计算。当然,对于其它材料,在动高压状态下的温度也可用这个方法估算。本文最后给出了适用于爆炸合成金刚石中对温度作工程估算用的解析拟合式。     

Effects of nano-sized aluminum on detonation characteristics and metal acceleration for RDX-based aluminized explosive

Nano-sized aluminum(Nano-Al)powders hold promise in enhancing the total energy of explosives and the metal acceleration ability at the same time.However,the near-detonation zone effects of reaction between Nano-Al with detonation products remain unclear.In this study,the overall reaction process of 170 nm Al with RDX explosive and its effect on detonation characteristics,detonation reaction zone,and the metal acceleration ability were comprehensively investigated through a variety of experiments such as the detonation velocity test,detonation pressure test,explosive/window interface velocity test and confined plate push test using high-resolution laser interferometry.Lithium fluoride(LiF),which has an inert behavior during the explosion,was used as a control to compare the contribution of the reaction of aluminum.A thermochemical approach that took into account the reactivity of aluminum and ensuing detonation products was adopted to calculate the additional energy release by afterburn.Combining the numerical simulations based on the calculated afterburn energy and experimental results,the param-eters in the detonation equation of state describing the Nano-Al reaction characteristics were calibrated.This study found that when the 170 nm Al content is from 0％to 15％,every 5％increase of aluminum resulted in about a 1.3％decrease in detonation velocity.Manganin pressure gauge measurement showed no significant enhancement in detonation pressure.The detonation reaction time and reaction zone length of RDX/Al/wax/80/15/5 explosive is 64 ns and 0.47 mm,which is respectively 14％and 8％higher than that of RDX/wax/95/5 explosive(57 ns and 0.39 mm).Explosive/window interface velocity curves show that 170 nm Al mainly reacted with the RDX detonation products after the detonation front.For the recording time of about 10 μs throughout the plate push test duration,the maximum plate velocity and plate acceleration time accelerated by RDX/Al/wax/80/15/5 explosive is 12％and 2.9 μs higher than that of RDX/LiF/wax/80/15/5,respectively,indicating that the aluminum reaction energy significantly increased the metal acceleration time and ability of the explosive.Numerical simulations with JWLM explosive equation of state show that when the detonation products expanded to 2 times the initial volume,over 80％of the aluminum had reacted,implying very high reactivity.These results are significant in attaining a clear understanding of the reaction mechanism of Nano-Al in the development of aluminized explosives.     

A quasi-isentropic model of a cylinder driven by aluminized explosives based on characteristic line analysis

A quasi-isentropic study on the process of driving a cylinder with aluminized explosives was carried out to examine the influence of the aluminum (Al) reaction rate on cylinder expansion and the physical parameters of the detonation products. Based on the proposed quasi-isentropic hypothesis and relevant isentropic theories, the characteristic lines of aluminized explosives driving a cylinder were analyzed, and a quasi-isentropic model was established. This model includes the variation of the cylinder wall velocity and the physical parameters of the detonation products with the Al reaction degree. Using previously reported experimental results, the quasi-isentropic model was verified to be applicative and accurate. This model was used to calculate the physical parameters for cylinder experiments with aluminized cyclotrimethylenetrinitramine explosives with 15.0 % and 30.0 % Al content. The results show that this quasi-isentropic model can be used not only to calculate the cylinder expansion rule or Al reaction degree, but also to calculate the physical parameters of the detonation products in the process of cylinder expansion. For explosives with 15.0 % and 30.0 % Al, 24.3 % and 18.5 % of the Al was found to have reacted at 33.9 μs and 34.0 μs, respectively. The difference in Al content results in different reaction intensity, occurrence time, and duration of two forms of reaction (diffusion and kinetic) between the Al powder and the detonation products; the post-detonation burning reaction between the Al powder and the detonation products prolongs the positive pressure action time, resulting in a continuous rise in temperature after detonation.     

Study of simple plane wave generator with an air-metal barrier

Plane wave generators (PWGs) are used to accelerate flyer plates to high velocities with their generated plane waves, which are widely used in the test of dynamic properties of materials. The traditional PWG is composed of two explosives with different detonation velocities. It is difficult to implement the related fabrication processes and control the generated waves due to its complicated structures. A simple plane wave generator is presented in this paper, which is composed of two identical cylindrical high explosive (HE) charges and an air-metal barrier. A theoretical model was established based on two different paths of the propagation of detonation waves, based on which the size of air-metal barrier was calculated for a given charge. The corresponding numerical simulations were also carried out by AUTODYN-2D^® based on the calculated results, which were used to compare with the theoretical calculations. A detonation wave with a flatness of 0.039 μs within the range of 70-percent diameter of the main charge was obtained through the simulations.     

Shock tube design for high intensity blast waves for laboratory testing of armor and combat materiel

Shock tubes create simulated blast waves which can be directed and measured to study blast wave effects under laboratory conditions. It is desirable to increase available peak pressure from ∼1 MPa to ∼5 MPa to simulate closer blast sources and facilitate development and testing of personal and vehicle armors. Three methods are experimentally investigated to increase peak simulated blast pressure produced by an oxy-acetylene driven shock tube while maintaining suitability for laboratory studies. The first method is the addition of a Shchelkin spiral priming section which supports a deflagration to detonation transition. This approach increases the average peak pressure from 1.17 MPa to 5.33 MPa while maintaining a relevant pressure-time curve (near Friedlander waveform). The second method is a bottleneck between the driving and driven sections. Coupling a 79 mm diameter driving section to a 53 mm driven section increases the peak pressure from 1.17 MPa to 2.25 MPa. A 103 mm driving section is used to increase peak pressure to 2.64 MPa. The third method, adding solid fuel to the driving section with the oxy-acetylene, results in a peak pressure increasing to 1.70 MPa.     

Application of Two-dimensional Viscous CE/SE Method in Calculation of Two-phase Detonation

The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.     

Numerical modelling of non-ideal detonation in ANFO explosives applying Wood-Kirkwood theory coupled with EXPLO5 thermochemical code

Ammonium nitrate and fuel oil (ANFO) based explosive is a classic example of non-ideal high explosives. Its detonation is characterized by a strong dependence of detonation parameters on explosive charge diameter, presence and characteristics of confinement, as well as incomplete consumption of explosive at the sonic point.In this work we propose a detonation model based on the Wood-Kirkwood (WK) theory coupled with the thermochemical code EXPLO5 and supplemented with reaction rate models. Our objective is to analyze the validity of the model for highly non-ideal ANFO explosives, with emphasis on effect of reaction rate models.It was found that both single-step and two-step pressure-based models can be calibrated to reproduce experimental detonation velocity-charge radius data of ANFO at radii significantly above the failure radius (i.e. for D/D_id > ∼0.6). Single-step pressure-based model, with the pressure exponent equal to 1.4, proved to be the most accurate, even in the vicinity of the failure radius. The impact of the rate models is most evident on temporal (and spatial) distribution of flow parameters in detonation driving zone, especially when it comes to the conversion and width of detonation driving zone.     

Experimental study on the detonation process of a pulse detonation engine with ionized seeds

An experimental platform of a pulse detonation engine (PDE) was established to study the effect of different K₂CO₃ ionized seed mass contents on the detonation process. The pressure and ion concentration were detected in the detonation process of the PDE with different contents of ionized seeds. The initiation process of the PDE at different ignition frequencies was studied. The results show that the gas conductivity in the detonation process increased by adding ionized seeds to the PDE tube, and the conductivity increased with the increase in ionized seed mass content. With the increase in ionized seed mass content, the range of the conductivity decreased. The PDE was successfully ignited and formed a stable detonation wave at ignition frequencies of 5 Hz and 10 Hz, and the peak pressure of the stable detonation with the ignition frequency of 5 Hz was 17% higher than that with an ignition frequency of 10 Hz. The detonation wave intensity was weakened and degenerated to a shock wave that propagated in the tube without the fuel filled at the ignition frequency of 20 Hz.     

金属圆管内爆轰波对碰部位速度的影响因素

利用LS-DYNA动力有限元程序对张世文等人做的金属圆管内爆轰波相互作用效应实验进行数值模拟,其模拟结果与实验结果具有较好的一致性。在此基础上着重分析了装药和金属圆管壁厚对金属圆管内爆轰波对碰部位速度的影响规律,为弹丸的设计提供了参考。     

Nitro-tetrazole based high performing explosives: Recent overview of synthesis and energetic properties

Heterocyclic skeleton (Azoles) and different energetic groups containing high performing explosives are highly emerged in recent years to meet the challenging requirements of energetic materials in both military and civilian applications with improved performance. For this purpose tetrazole (Azole) is identified as an attractive heterocyclic backbone with energetic functional groups nitro (-NO₂), nitrato (-ONO₂), nitrimino (-NNO₂), and nitramino (–NH–NO₂) to replace the traditionally used high performing explosives. The tetrazole based compounds having these energetic functional groups demonstrated advanced energetic performance (detonation velocity and pressure), densities, and heat of formation (HOF) and became a potential replacement of traditional energetic compounds such as RDX. This review presents a summary of the recently reported nitro-tetrazole energetic compounds containing poly-nitro, di/mono-nitro, nitrato/nitramino/nitrimino, bridged/bis/di tetrazole and nitro functional groups, describing their preparation methods, advance energetic properties, and further applications as high-performing explosives, especially those reported in the last decade. This review aims to provide a fresh concept for designing nitro-tetrazole based high performing explosives together with major challenges and perspectives.     

Effects of total pressures and equivalence ratios on kerosene/air rotating detonation engines using a paralleling CE/SE method

In this paper, the kerosene/air rotating detonation engines(RDE) are numerically investigated, and the emphasis is laid on the effects of total pressures and equivalence ratios on the operation characteristics of RDE including the initiation, instabilities, and propulsive performance. A hybrid MPI + OpenMP parallel computing model is applied and it is proved to be able to obtain a more effective parallel performance on high performance computing(HPC) systems. A series of cases with the total pressure of 1 MPa, 1.5 MPa, 2 MPa, and the equivalence ratio of 0.9, 1, 1.4 are simulated. On one hand, the total pressure shows a significant impact on the instabilities of rotating detonation waves. The instability phenomenon is observed in cases with low total pressure (1 MPa) and weakened with the increase of the total pressure. The total pressure has a small impact on the detonation wave velocity and the specific impulse. On the other hand, the equivalence ratio shows a negligible influence on the instabilities, while it affects the ignition process and accounts for the detonation velocity deficit. It is more difficult to initiate rotating detonation waves directly in the lean fuel operation condition. Little difference was observed in the thrust with different equivalence ratios of 0.9, 1, and 1.4. The highest specific impulse was obtained in the lean fuel cases, which is around 2700 s. The findings could provide insights into the understanding of the operation characteristics of kerosene/air RDE.