Experimental study on the cavity evolution and liquid spurt of hydrodynamic ram

In research of the characteristics of the cavity evolution, the pressure, and the liquid spurt in hydrodynamic ram, the experiment of the high-velocity fragment impacting the water-filled container had been conducted. The relationships between the above three characteristics have been researched. The evolution of the cavity can be divided into three processes according to its shape characteristics. The first liquid spurt occurred in Process Ⅱ and the rest of it occurred in Process Ⅲ. The duration of the second liquid spurt is longer than the first liquid spurt. When the impact velocity of the fragment is less than 996 m/s, the velocity of the second liquid spurt is the highest. When the velocity of the fragment is greater than 996 m/s, the velocity of the first liquid spurt is the highest. The maximum velocities of the first and second liquid spurt are 111 m/s and 94 m/s respectively. The pressure fluctuated sharply in Processes Ⅰ and Ⅲ. The maximum peak pressures in the shock and the cavity oscillation phases are 15.51 MPa and 7.96 MPa respectively. The time interval of the two adjacent pressure pulses increases with the increase of the fragment velocity.     

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.     

驱动管中柱状装药爆轰过程的数值模拟

数值模拟了爆炸驱动管中柱状装药内爆轰波的传播过程.计算采用欧拉型有限体积方法,炸药及爆轰产物均采用JWL状态方程,空气采用理想气体状态方程,采用"点火-生长"模型计算化学反应速率.计算得到了驱动管内波系结构的发展过程,爆速与经验公式符合得较好.计算表明,驱动管侧壁的压力峰值在800MPa以上,而在管底中心处,由于激波的汇聚,压力峰值高达12.4GPa.     

Numerical investigation on free air blast loads generated from center-initiated cylindrical charges with varied aspect ratio in arbitrary orientation

In current guidelines, the free air blast loads (overpressure and impulse) are determined by spherical charges, although most of ordnance devices are more nearly cylindrical than spherical in geometry. This may result in a great underestimation of blast loads in the near field and lead to an unsafe design. However, there is still a lack of systematic quantitative analysis of the blast loads generated from cylindrical charges. In this study, a numerical model is developed by using the hydrocode AUTODYN to investigate the influences of aspect ratio and orientation on the free air blast loads generated from center-initiated cylindrical charges. This is done by examining the pressure contours, the peak overpressures and impulses for various aspect ratios ranged from 1 to 8 and arbitrary orientation monitored along every azimuth angle with an interval of 5°. To characterize the distribution patterns of blast loads, three regions, i.e., the axial region, the vertex region and the radial region are identified, and the propagation of blast waves in each region is analyzed in detail. The complexity of blast loads of cylindrical charges is found to result from the bridge wave and its interaction with primary waves. Several empirical formulas are presented based on curve-fitting the numerical data, including the orientation where the maximum peak overpressure emerges, the critical scaled distance beyond which the charge shape effect could be neglected and blast loads with varied aspect ratio in arbitrary orientation, all of which are useful for blast-resistant design.     

A fast-running method for blast load prediction shielding by a protective barrier

This study presents a simplified blast load prediction method on structures behind a protective barrier. The proposed method is basically an empirical approach based on Kingery-Bulamsh (K-B) chart and finite element (FE) analysis results. To this end, this study divides the structure into three regions by three critical points. Blast loads at each critical point can be calculated based on K-B chart and an approximation according to FE analysis results. Finally, peak reflected overpressure and impulse distributed on the structure can be approximately estimated by linearly connecting blast loads at each critical point. In order to confirm a feasibility of the proposed method, a series of numerical simulations were carried out. The simulation results were compared with FE analysis results which are presented in the open literature. From such comparisons, it was found that the proposed method is applicable to predict blast loads on structures behind a protective barrier.     

Research on the influences of confining pressure and strain rate on NEPE propellant: Experimental assessment and constitutive model

In order to study the influences of confining pressure and strain rate on the mechanical properties of the Nitrate Ester Plasticized Polyether (NEPE) propellant, uniaxial tensile tests were conducted using the self-made confining pressure system and material testing machine. The stress-strain responses of the NEPE propellant under different confining pressure conditions and strain rates were obtained and analyzed. The results show that confining pressure and strain rate have a remarkably influence on the mechanical responses of the NEPE propellant. As confining pressure increases (from 0 to 5.4 MPa), the maximum tensile stress and ultimate strain increase gradually. With the coupled effects of confining pressure and strain rate, the value of the maximum tensile stress and ultimate strain at 5.4 MPa and 0.0667 s⁻¹ is 2.03 times and 2.19 times of their values under 0 MPa and 0.00333 s⁻¹, respectively. Afterwards, the influence mechanism of confining pressure on the NEPE propellant was analyzed. Finally, based on the viscoelastic theory and continuous damage theory, a nonlinear constitutive model considering confining pressure and strain rate was developed. The damage was considered to be rate-dependent and pressure-dependent. The constitutive model was validated by comparing experimental data with predictions of the constitutive model. The whole maximum stress errors of the model predictions are lower than 4% and the corresponding strain errors are lower than 7%. The results show that confining pressure can suppress the damage initiation and evolution of the NEPE propellant and the nonlinear constitutive model can describe the mechanical responses of the NEPE propellant under various confining pressure conditions and strain rates. This research can lay a theoretical foundation for analyzing the structural integrity of propellant grain accurately under working pressure loading.     

命中角对防护结构口部爆炸毁伤影响的数值分析

为研究装药命中角对防护结构口部爆炸毁伤的影响规律,对直墙拱结构进行了一系列爆炸模拟分析,得到了直墙拱结构拱顶、拱脚和直墙底3个位置的爆炸应力和位移时程曲线。通过分析这3个位置的应力峰值和位移峰值,发现随着装药命中角的增大,拱顶处竖直方向荷载减小,而拱脚处水平方向荷载增大;随着装药命中角的增大,拱顶处竖直方向位移变形减小,而拱脚处水平方向位移变形增大。模拟分析得到的命中角对防护结构口部爆炸毁伤的影响规律,可供防护工程设计时参考。     

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.     

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

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

燃烧室压力对潜入式喷管喉衬热应力的影响

为了研究燃烧室压力对固体火箭发动机潜入式喷管热应力影响规律的问题,采用商业流体软件,基于压力求解器,求解了喷管纯气相的流场,确定了燃气温度、压力、壁面对流换热系数;采用有限元软件,依据流场计算的非均布壁面压力与非均布对流换热,求解了燃烧室压力为6 MPa下的潜入式喷管热结构问题;通过地面点火试验验证了仿真模型与数值方法的有效性与准确性;采用相同计算模型与数值方法,求解了在燃烧室压力为9 MPa、12 MPa下的喷管热结构问题,揭示了燃烧室压力对喉衬热应力的影响规律。结果表明:整个工作过程,喉衬环向应力最大值为103.9 MPa,位于内表面,且随时间增大,先增大后减小;喉衬环向拉应力也随时间先增大后减小;随压力增大,对流换热系数增大,喉衬温度升高,喉衬环向拉应力增大,喉衬环向压应力减小。     

Novel aluminum-based fuel: Facile preparation to improve thermal reactions

To improve the thermal properties of aluminum (Al) in the energetic system, a coated structure with ammonium perchlorate (AP) was prepared by a facile approach. And N, N-Dimethylformamide (DMF) was chosen as an ideal solvent based on heterogeneous nucleation theory and molecular dynamics simulation. This coated structure could enlarge the contact area and improve the reaction environment to enhance the thermal properties. The addition of AP could accelerate oxidation temperature of Al with around 17.5 °C. And the heat release of 85@15 composition rises to 26.13 kJ/g and the reaction degree is 97.6% with higher peak pressure (254.6 kPa) and rise rate (1.397 MPa/s). An ideal ratio with 15 wt% AP was probed primarily. The high energy laser-induced shockwave experiment was utilized to simulate the reaction behavior in hot field. And the larger activated mixture of coated powder could release more energy to promote the growth of shockwave with higher speed up to 518.7 ± 55.9 m/s. In conclusion, 85@15 composition is expected to be applied in energetic system as a novel metal fuel.     

Numerical simulation of the blast wave of a multilayer composite charge

The numerical simulation of a blast wave of a multilayer composite charge is investigated. A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is established using the AUTODYN finite element program. Results of the near-field and far-field calculations of the shock wave respectively converge at cell sizes of 0.25–0.5 cm and 1–3 cm. The Euler––flux-corrected transport solver is found to be suitable for the far-field calculation after mapping. A numerical simulation is conducted to study the formation, propagation, and interaction of the shock wave of the composite charge for different initiation modes. It is found that the initiation mode obviously affects the shock-wave waveform and pressure distribution of the composite charge. Additionally, it is found that the area of the overpressure distribution is greatest for internal and external simultaneous initiation, and the peak pressure of the shock wave exponentially decays, fitting the calculation formula of the peak overpressure attenuation under different initiation modes, which is obtained and verified by experiment. The difference between numerical and experimental results is less than 10%, and the peak overpressure of both internal and external initiation is 56.12% higher than that of central single-point initiation.     

Relative performance of novel blast wave mitigation system to conventional system based on mitigation percent criteria

Recent researches focused on developing robust blast load mitigation systems due to the threats of terrorist attacks. One of the main embraced strategies is the structural systems that use mitigation techniques. They are developed from a combination of structural elements and described herein as conventional systems. Among the promising techniques is that redirect the waves propagation through hollow tubes. The blast wave propagation through tubes provides an efficient system since it combines many blast wave phenomena, such as reflection, diffraction, and interaction. In this research, a novel blast load mitigation system, employed as a protection fence, is developed using a technique similar to the technique of the bent tube in manipulating the shock-wave. The relative performance of the novel system to the conventional system is evaluated based on mitigation percent criteria. Performances of both systems are calculated through numerical simulation. The proposed novel system proved to satisfy high performance in mitigating the generated blast waves from charges weight up to 500 kg TNT at relatively small standoff distances (5 m and 8 m). It mitigates at least 94％ of the blast waves, which means that only 6％ of that blast impulse is considered as the applied load on the targeted structure.     

T型坑道中爆炸冲击波传播规律的数值模拟

为了确定地下坑道中防护门上的冲击波荷载,首先需要弄清爆炸冲击波在坑道中的传播衰减规律.在模型坑道爆炸实验数据验证的基础上,应用数值模拟方法研究了T型坑道中爆炸冲击波的传播衰减规律.通过大量不同模型尺寸、不同炸药药量的数值模拟,得到了T型坑道中平面冲击波的形成位置,提出了爆炸冲击波通过T型坑道拐弯时新的冲击波压力衰减系数...     

Responses of HFR-LWC beams under close-range blast loadings accompanying membrane action

The load-carrying capacities and failure patterns of reinforced concrete components can be significantly changed by membrane effects. However, limited work has been carried out to investigate the blast resistance of Hybrid Fiber Reinforced Lightweight Aggregate Concrete (HFR-LWC) members accompanying membrane action. This paper presents a theoretical approach to quantitatively depicting the membrane behavior and its contribution on the behavior of HFR-LWC beams under close-range blast loadings, and the suitability of the proposed model is validated by a series of field tests. An improved Single-Degree-of-Freedom (SDOF) model was employed to describe the dynamic responses of beam-like members under blast loadings accompanying membrane action, where the mass-load coefficient is determined according to the nonuniformly distributed load induced by close-range explosion, and the membrane action is characterized by an in-plane (longitudinal) force and a resisting moment. The elasto-plastic and recovery responses of HFR-LWC beams under the combined action of blast load and membrane force were analyzed by the promoted model. A specially built end-constrain clamp was developed to provide membrane action for the beam member when they are subjected to blast load simultaneously. It is demonstrated that the analytical displacement-time histories are in good agreement with experimental results before peak deflections and that the improved SDOF model is an acceptable tool for predicting the behavior of HFR-LWC beams under blast loadings accompanying membrane action.     

Evaluation of effectiveness of polymer coatings in reducing blast-induced deformation of steel plates

Incorporating elastomers such as polymers in protective structures to withstand high energetic dynamic loads, has gained significant interest. The main objective of this study is to investigate the influence of a Polyurea coating towards the blast-induced response in steel plates. As such, Polyurea coated steel plates were tested under near-field blast loads, produced by the detonation of 1 kg of spherical nitromethane charges, at a standoff distance of 150 mm. Mild steel (XLERPLATE 350) and high-strength steel (BIS80) plates with thicknesses of 10 mm were Polyurea coated with thicknesses of 6 mm and 12 mm on either the front (facing the charge) or the back face. The deformation profiles were measured using 3D scanning. Numerical simulations were performed using the non-linear finite element code LS-DYNA. The strain-dependent behaviour of the steel and Polyurea were represented by Johnson-cook and Money-Rivlin constitutive models, respectively. The numerical models were validated by comparing the plate deflection results obtained from the experiments and were then used in the subsequent parametric study to investigate the optimum thickness of the Polyurea coating. The results indicate that back face coating contributes towards an approximately 20% reduction in the residual deformation as well as the absence of melting of the Polyurea layer, while the front-face coating can be used a means of providing additional standoff distance to the steel plates.     

Numerical simulation of base flow with hot base bleed for two jet models

In order to improve the benefits of base bleed in base flow field, the base flow with hot base bleed for two jet models is studied. Two-dimensional axisymmetric Navier–Stokes equations are computed by using a finite volume scheme. The base flow of a cylinder afterbody with base bleed is simulated. The simulation results are validated with the experimental data, and the experimental results are well reproduced. On this basis, the base flow fields with base bleed for a circular jet model and an annulus jet model are investigated by selecting the injection temperature from 830 K to 2200 K. The results show that the base pressure of the annular jet model is higher than that of the circular jet model with the changes of the injection parameter and the injection temperature. For the circular jet model, the hot gases are concentrated in the vicinity of the base. For the annular jet model, the bleed gases flow into the shear layer directly so that the hot gases are concentrated in the shear layer. The latter temperature distribution is better for the increase of base pressure.     

螺旋桨激励水下艇体振动的试验及数值研究

为考虑螺旋桨真实激励特性,在循环水槽中开展艇尾伴流场中螺旋桨诱导艇尾脉动压力及螺旋桨激励水下艇体振动响应的测量试验。试验结果表明:脉动压力幅值在叶频处最大,且随螺旋桨负载增加而增大,随与螺旋桨之间距离的增大而减小,四叶桨脉动压力在尾翼后的高伴流区幅值较大,五叶桨则在尾翼之间的低伴流区幅值较大；大部分测点的振动响应幅值随螺旋桨负载增加而增大,但也存在叶频处幅值较小和未随负载增加而增大的情况；五叶桨激励引起的侧向振动较四叶桨有所增强而轴向振动有所减弱；特定谱峰频率处振动响应幅值呈一阶弯曲振型,其频率范围与有限元计算结果较为一致。综合采用计算流体动力学、有限元和模态叠加法建立螺旋桨激励水下艇体振动响应的数值计算方法,通过计算与试验结果的对比发现,该方法的计算结果与试验结果吻合较好,相比于采用单位简谐激励的谐响应分析方法更加接近真实情况。     

Mitigation of the blast load effects on a building structure using newly composite structural configurations

In this study, a nonlinear three-dimensional hydrocode numerical simulation was carried out using AUTODYN-3D to investigate the effect of blasting of a high explosive material (TNT) against several configurations of the composite structure. Several numerical models were carried out to study the effect of varying the thickness of the walls and the effect of adding an air layer or aluminum foam layer inside two layers of concrete in mitigating the effect of blast waves on the structure walls. The results showed that increasing the thickness of walls has a good effect on mitigating the effect of blast waves. When a layer of air was added, the effect of blast waves was exaggerated, while when a layer of aluminum foam was added the blast wave effects were mitigated with a reasonable percentage.