Research and development on hypervelocity impact protection using Whipple shield: An overview

Whipple shield, a dual-wall system, as well as its improved structures, is widely applied to defend the hypervelocity impact of space debris (projectile). This paper reviews the studies about the mechanism and process of protection against hypervelocity impacts using Whipple shield. Ground-based experiment and numerical simulation for hypervelocity impact and protection are introduced briefly. Three steps of the Whipple shield protection are discussed in order, including the interaction between the projectile and bumper, the movement and diffusion of the debris cloud, and the interaction between the debris cloud and rear plate. Potential improvements of the protection performance focusing on these three steps are presented. Representative works in the last decade are mentioned specifically. Some prospects and suggestions for future studies are put forward.     

填充式波纹夹层结构超高速撞击特性仿真

基于航天器空间碎片被动防护需求,对一种新型填充式波纹夹层结构进行超高速撞击仿真研究,分析超高速撞击过程以及结构的穿孔破坏情况和所形成的碎片云的特性,并与相同面密度Whipple结构进行对比。其撞击现象与Whipple结构相似,但其碎片云的头部速度小于Whipple结构,而径向膨胀最大速度和膨胀半角均大于Whipple结构。随撞击初速从3 km/s~10 km/s不断增大,波纹夹层结构的撞击穿孔尺寸变大,形状也更不规则。此外,结构中的填充树脂对碎片撞击能量的吸收贡献最大,后面板所吸收的能量所占比重较大,而前面板和波纹板对碎片撞击能量的吸收贡献较小。研究结果对空间碎片防护结构的设计具有一定的参考意义。     

Characteristics structure analysis on debris cloud in the hypervelocity impact of disk projectile on thin plate

In this paper, the gauge points setting is introduced in the SPH simulation to analyze the debris cloud structure generated by the hypervelocity impact of disk projectile on thin plate. Compared with the experiments, more detailed information of the debris cloud structure can be classified from the numerical simulation. However, due to the solitary dispersion and overlap display of the particles in the SPH simulation, accurate comparison between numerical and experimental results is difficult to be performed. To track the velocity and spatial distribution of the particles in the debris cloud induced from disk and plate, gauge points are locally set in the single-layer profile in the SPH model. By analyzing the gauge points’ spatial coordinate and velocity, the location and velocity of characteristic points in the debris cloud are determined. The boundary of debris cloud is achieved, as well as the fragments distribution outside the main structure of debris cloud.     

高速/超高速碰撞效应研究

高速/超高速碰撞会产生大量破片,并伴随闪光、等离子体等现象.高速碰撞材料受冲击产生射流及破片高温会引起碰撞闪光现象,获得了碰撞等离子体电磁场强度的变化规律.采用数值仿真与理论分析相结合的方法,获得了高速碰撞破片云团分布规律.     

B4C/Al复合板中应力波行为分析（Ⅰ）

依据7.62mm穿甲燃烧弹与B4C,Al板的物理和力学特性,提出弹、板的力学模型;在此基础上,给出弹板碰撞后B4C/Al复合板的弹性动力学方程,建立应力波的传播模型,讨论了应力波传播时复合板的力学行为,最后确定了复合板承受的初始冲击应力。理论分析表明：一维应变条件的B4C和Al板的弹性极限高于一维应力条件下的值;弹、板间的作用力服从指数衰减规律;B4C板的初始应力高于430MPa时,应力波的作用使Al板在卸载过程中发生反向屈服。     

Influence of a liquid-filled compartment structure on the incoming shaped charge jet stability

Liquid-filled compartment structure consists of a bulk steel plate with matrix blind holes which are filled with liquid and a steel front plate to seal up the liquid with rings and bolts.The liquid-filled compart-ment structure can resist the shaped charge warhead effectively.This paper presents experimental and theoretical investigations of the penetration ability of the residual shaped charge jet emerging from the liquid-filled compartment structure after the penetration process at different impact angles.On the basis of shock wave propagation theory,the influence of the liquid-filled compartment structure on jet sta-bility is analysed.The interferences of the liquid backflow caused by a reflected shock wave and a back plate on jet stability under different impact angles are also examined.In addition,the range of the disturbed velocity segments of the jet at different impact angles and the penetration ability of the re-sidual jet are obtained.A theoretical model is validated against the experimental penetration depths.     

SPH算法在超高速碰撞数值模拟中的应用

采用大型动力学软件Ls-dyna中的SPH求解器,对球形弹丸超高速碰撞靶板的过程进行数值模拟。给出由于碰撞速度不同,而引起的不同的碎片云图像,以及不同的靶板破碎孔的尺寸。同时模拟了入射角度不同,而得到不同的超高速碰撞物理过程。初步分析结果表明数值模拟结果是合理的,较清晰地模拟了超高速碰撞下碎片云的生成、膨胀过程。另外,通过数值模拟结果与实验结果的比较分析,表明SPH方法得到的结果能描述相应的超高速碰撞现象。所以采用SPH算法针对超高速碰撞物理过程进行数值模拟是可行的。     

二维垂向梯度非均匀介质中弹性波传播特性

研究了二维弹性波在宏观非均匀介质中传播的问题。由于介质物理性质的非均匀性,弹性波在该介质中的传播需用包含介质物理性质的空间导数的非均匀介质波动方程来描述。采用有限元法求解波动方程,重点研究和讨论了物理性质垂向梯度渐变的非均匀介质中弹性波的传播特点。结果显示:波速在垂向上的对称的递增或递减引起平面弹性波的会聚或发散;若波速在垂向上满足双曲余割分布,则弹性波可以在对称轴上周期性聚焦。     

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.     

Shock-induced reaction behaviors of functionally graded reactive material

In this paper, the ballistic impact experiments, including impact test chamber and impact double-spaced plates, were conducted to study the reaction behaviors of a novel functionally graded reactive material (FGRM), which was composed of polytetrafluoroethylene/aluminum (PTFE/Al) and PTFE/Al/bismuth trioxide (Bi₂O₃). The experiments showed that the impact direction of the FGRM had a significant effect on the reaction. With the same impact velocity, when the first impact material was PTFE/Al/Bi₂O₃, compared with first impact material PTFE/Al, the FGRM induced higher overpressure in the test chamber and larger damaged area of double-spaced plates. The theoretical model, which considered the shock wave generation and propagation, the effect of the shock wave on reaction efficiency, and penetration behaviors, was developed to analyze the reaction behaviors of the FGRM. The model predicted first impact material of the FGRM with a higher shock impedance was conducive to the reaction of reactive materials. The conclusion of this study provides significant information about the design and application of reactive materials.     

海水中电磁波传播特性的研究

利用电磁场传播所满足的Maxwell方程组,计算和分析出电磁波在导电媒质中传播时的特征;并以海水为例,得出一些有意义的结论,为海水中通信、信号探测、引信研究等方面工作提供理论依据.     

Experimental and numerical study of the blast wave decrease using sandwich panel by granular materials core

Among the intrinsic properties of some materials, e.g., foams, porous materials, and granular materials, are their ability to mitigate shock waves. This paper investigated shock wave mitigation by a sandwich panel with a granular core. Numerical simulations and experimental tests were performed using Autodyn hydro-code software and a shock tube, respectively. The smoothed particle hydrodynamics (SPH) method was used to model granular materials. Sawdust and pumice, whose properties were determined by several compression tests, were used as granular materials in the sandwich panel core. These granular materials possess many mechanisms, including compacting (e.g., sawdust) and crushing (e.g., pumice) that mitigate shock/blast wave. The results indicated the ineffectiveness of using a core with low thickness, yet it was demonstrated to be effective with high thickness. Low-thickness pumice yielded better results for wave mitigation. The use of these materials with a core with appropriate core reduces up to 88% of the shock wave. The results of the experiments and numerical simulations were compared, suggesting a good agreement between the two. This indicates the accuracy of simulation and the ability of the SPH method to modeling granular material under shock loading. The effects of grain size and the coefficient of friction between grains have also been investigated using simulation, implying that increasing the grain size and coefficient of friction between grains both reduce overpressure.     

Modeling of the whole process of shock wave overpressure of free-field air explosion

The waveform of the explosion shock wave under free-field air explosion is an extremely complex problem. It is generally considered that the waveform consists of overpressure peak, positive pressure zone and negative pressure zone. Most of current practice usually considers only the positive pressure. Many empirical relations are available to predict overpressure peak, the positive pressure action time and pressure decay law. However, there are few models that can predict the whole waveform. The whole process of explosion shock wave overpressure, which was expressed as the product of the three factor functions of peak, attenuation and oscillation, was proposed in the present work. According to the principle of explosion similarity, the scaled parameters were introduced and the empirical formula was absorbed to form a mathematical model of shock wave overpressure. Parametric numerical simulations of free-field air explosions were conducted. By experimental verification of the AUTODYN numerical method and comparing the analytical and simulated curves, the model is proved to be accurate to calculate the shock wave overpressure under free-field air explosion. In addition, through the model the shock wave overpressure at different time and distance can be displayed in three dimensions. The model makes the time needed for theoretical calculation much less than that for numerical simulation.     

Investigation on the spatial distribution characteristics of behind-armor debris formed by the perforation of EFP through steel target

The behind-armor debris (BAD) formed by the perforation of an EFP is the main damage factor for the secondary destruction to the behind-armor components. Aiming at investigating the BAD caused by EFP, flash X-ray radiography combined with an experimental witness plate test method was used, and the FEM-SPH adaptive conversion algorithm in LS-DYNA software was employed to model the perforation process. The simulation results of the debris cloud shape and number of debris were in good agreement with the flash X-ray radiographs and perforated holes on the witness plate, respectively. Three-dimensional numerical simulations of EFP's penetration under various impact conditions were conducted. The results show that, an ellipsoidal debris cloud, with the major-to-minor axis radio (a/b) smaller than that caused by shaped charge jets, was formed behind the target. With the increase of target thickness (h) and decrease of impact velocity (v₀) and obliquity (θ), the value of a/b decreases. The number of debris ejected from target is significantly higher than that from EFP. Based on the statistical analysis of the spatial distribution of the BAD, An engineering calculation model was established considering the influence of h, v₀ and θ. The model can with reasonable accuracy predict the quantity and velocity distribution characteristics of BAD formed by EFP.     

Shock wave mitigation using zig-zag structures and cylindrical obstructions

The present study focuses on the mitigation of shock wave using novel geometric passages in the flow field. The strategy is to produce multiple shock reflections and diffractions in the passage with minimum flow obstruction, which in turn is expected to reduce the shock wave strength at the target location. In the present study the interaction of a plane shock front (generated from a shock tube) with various geometric designs such as, 1) zig-zag geometric passage, 2) staggered cylindrical obstructions and 3) zig-zag passage with cylindrical obstructions have been investigated using computational technique. It is seen from the numerical simulation that, among the various designs, the maximum shock attenuation is produced by the zig-zag passage with cylindrical obstructions which is then followed by zig-zag passage and staggered cylindrical obstructions. A comprehensive investigation on the shock wave reflection and diffraction phenomena happening in the proposed complex passages have also been carried out. In the new zig-zag design, the initial shock wave undergoes shock wave reflection and diffraction process which swaps alternatively as the shock front moves from one turn to the other turn. This cyclic shock reflection and diffraction process helps in diffusing the shock wave energy with practically no obstruction to the flow field. It is found that by combining the shock attenuation ability of zig-zag passage (using shock reflection and diffraction) with the shock attenuation ability of cylindrical blocks (by flow obstruction), a drastic attenuation in shock strength can be achieved with moderate level of flow blocking.     

双层地下室空气冲击波特性的模型试验研究

高层建筑多层地下室内人防结构顶板核爆荷载的合理确定是目前人防结构设计中亟待解决的问题。文中在对实际工程的调查基础上,研究了双层地下室的试验模型,进行了空气冲击波在地下室模型内传播的爆炸试验,分析了作用在双层地下室底层结构顶板的爆炸压力分布特性,并得出多层地下室的底层顶板核爆荷载一般要小于单层地下室顶板核爆荷载设计值的结论.     

基于水下脉冲放电的冲击波聚焦实验研究

冲击波的聚焦可以在聚焦区域形成局部较高的压力.基于能产生较强冲击波的水中脉冲放电声源和具有会聚作用的旋转椭球面反射罩建立了水下冲击波聚焦系统,开展了水下冲击波聚焦的实验研究,分析了水下冲击波的聚焦过程.对压力历史、轴线上的压力分布及峰值压力进行了分析,研究了旋转椭球面反射罩的聚焦特性.实验发现聚焦区域会产生负压,并导致局部空化.基于冲击动力学理论对这一现象进行了分析.实验结果表明,水下冲击波聚焦系统具有显著的聚焦效果.     

高速方尾水面舰船兴波问题计算方法研究

将Dawson型兴波阻力计算方法推广用于高速方尾水面舰船兴波问题计算,在计算中同时计入了尾封板0静压力、尾倾与下沉等因素的影响.以两高速方尾水面舰船为例计算了兴波问题,将兴波阻力计算值与剩余阻力模型试验值进行了比较,求出了相应的相关系数.结果表明,本方法可较好地求解高速方尾水面舰船的兴波问题,兴波阻力相关系数可用于修正理论计算结果,以提高计算精度.     

水下近距爆炸作用下弹性钢板处的空化特性研究

弹性钢板在水下近距爆炸作用下,冲击波会使其附近流体形成局部空化,脉动气泡会使流体形成锥形空化。利用平面冲击波理论对局部空化的形成特性进行了研究,理论分析了结构目标尺度的变化对空化区域形成的影响,并通过具体试验对局部空化理论进行了验证,两者符合较好;通过试验和数值仿真方法研究了气泡脉动引起的锥形空化的形成特点,初步分析了锥形空化的形成原因。结果表明,锥形空化对结构具有较大的冲击作用。     

超声导波在管道中传播的可视化模拟研究

利用有限元方法建立管道导波3-D模型,通过在管道中激励纵向超声波,模拟导波在管道中的传播形态。根据回波信号的传播时间计算了纵波传播速度,与频散曲线得到的纵波速度符合较好。从观测点取得的波信号能够反映超声波在管道中传播时存在的反射和衰减现象,动画结果能够形象直观地反映导波在管道中传播的情形,实现了导波传播的可视化模拟。该研究有利于更好地解释导波在管道无损检测中出现的频散特征。