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

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

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

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

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.     

Analysis of the stress wave and rarefaction wave produced by hypervelocity impact of sphere onto thin plate

Shock wave is emitted into the plate and sphere when a sphere hypervelocity impacts onto a thin plate. The fragmentation and phase change of the material caused by the propagation and unloading of shock wave could result in the formation of debris cloud eventually. Propagation models are deduced based on one-dimensional shock wave theory and the geometry of sphere, which uses elliptic equations (corresponding to ellipsoid equations in physical space) to describe the propagation of shock wave and the rarefaction wave. The “Effective thickness” is defined as the critical plate thickness that ensures the rarefaction wave overtake the shock wave at the back of the sphere. The “Effective thickness” is directly related to the form of the debris cloud. The relation of the “Effective thickness” and the “Optimum thickness” is also discussed. The impacts of Al spheres onto Al plates are simulated within SPH to verify the propagation models and associated theories. The results show that the wave fronts predicted by the propagation models are closer to the simulation result at higher impact velocity. The curvatures of the wave fronts decrease with the increase of impact velocities. The predicted “Effective thickness” is consistent with the simulation results. The analysis about the shock wave propagation and unloading in this paper can provide a new sight and inspiration for the quantitative study of hypervelocity impact and space debris protection.     

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.     

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

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

SPH的两种并行计算方案及其在自由表面流动模拟中的应用

针对光滑粒子动力学(SPH)主要计算量是近邻粒子搜索这一特点,提出了一种基于粒子分解的SPH并行计算方案。利用该方案可以方便的将任意串行SPH代码并行计算,而且每一个时间步内的信息传递量只和粒子总数有关,而和粒子的分布无关,因而特别适合于自由表面流动等大变形问题的并行数值模拟。对一个粒子总数为40万的三维溃坝问题的模拟结果表明,此方案能达到的最大加速比约为16,这一结果可能比空间分解方案(不考虑动态负载均衡)更优。     

Simulation of natural fragmentation of rings cut from warheads

Natural fragmentation of warheads that detonates causes the casing of the warhead to split into various sized fragments through shear or radial fractures depending on the toughness,density,and grain size of the material.The best known formula for the prediction of the size distribution is the Mott formulae,which is further examined by Grady and Kipp by investigating more carefully the statistical most random way of portioning a given area into a number of entities.We examine the fragmentation behavior of radially expanding steel rings cut from a 25 mm warhead by using an in house smooth particle hydrodynamic(SPH) simulation code called REGULUS.Experimental results were compared with numerical results applying varying particle size and stochastic fracture strain.The numerically obtained number of fragments was consistent with experimental results.Increasing expansion velocity of the rings increases the number of fragments.Statistical variation of the material parameters influences the fragment characteristics,especially for low expansion velocities.A least square regression fit to the cumulative number of fragments by applying a generalized Mott distribution shows that the shape parameter is around 4 for the rings,which is in contrast to the Mott distribution with a shape parameter of 1/2.For initially polar distributed particles,we see signs of a bimodal cumulative fragment distribution.Adding statistical variation in material parameters of the fracture model causes the velocity numerical solutions to become less sensitive to changes in resolution for Cartesian distributed particles.     

爆炸成型弹丸贯穿靶板靶后破片空间分布模型

深空撞击载荷对探测小天体内部物质成分和结构特性具有重要意义,因此,在考虑爆炸成型弹丸(explosively formed projectile, EFP)可变截面的特性和区分靶后破片来源的基础上,建立了EFP贯穿靶板靶后破片空间分布模型。在靶板厚度30 mm至70 mm、EFP着靶速度1 650 m/s至1 860 m/s的条件下,该模型可以定量预测靶后破片云中各个破片的速度、质量、数量与空间位置的关系。结果表明,相对速度总是随相对空间位置的增加而呈线性增加,相对质量、相对数量总是随相对空间位置的增加而呈幂函数增加。     

Interface defeat studies of long-rod projectile impacting on ceramic targets

The interface defeat phenomenon always occurs when a long-rod projectile impacting on the ceramic target with certain velocity, i.e., the projectile is forced to flow radially on the surface of ceramic plates for a period of time without significant penetration. Interface defeat has a direct effect upon the ballistic performance of the armor piercing projectile, which is studied numerically and theoretically at present. Firstly, by modeling the projectiles and ceramic targets with the SPH (Smoothed Particle Hydrodynamics) particles and Lagrange finite elements, the systematic numerical simulations on interface defeat are performed with the commercial finite element program AUTODYN. Three different responses, i.e., complete interface defeat, dwell and direct penetration, are reproduced in different types of ceramic targets (bare, buffered, radially confined and oblique). Furthermore, by adopting the validated numerical algorithms, constitutive models and the corresponding material parameters, the influences of projectile (material, diameter, nose shape), constitutive models of ceramic (JH-1 and JH-2 models), buffer and cover plate (thickness, constraints, material), as well as the prestress acted on the target (radial and hydrostatic) on the interface defeat (transition velocity and dwell time) are systematically investigated. Finally, based on the energy conservation approach and taking the strain rate effect of ceramic material into account, a modified model for predicting the upper limit of transition velocity is proposed and validated. The present work and derived conclusions can provide helpful reference for the design and optimization of both the long-rod projectile and ceramic armor.     

SPH数值模拟中固壁边界的一种处理方法

在光滑粒子流体动力学 (SPH)数值模拟中尝试了一种处理固壁边界的边界力方法 ,给出了一种新的边界力形式。利用SPH方法及边界力方法对水坝坍塌和涌波进入静止水塘这两个自由表面流动问题作了数值模拟。数值模拟结果表明 ,在SPH计算中使用本文所给边界力处理固壁是行之有效的     

基于FEM/SPH算法的弹丸侵彻土体数值模拟

基于有限元法(FEM)和光滑粒子流体动力学(SPH)结合的算法,用数值模拟了钢质弹丸对钢筒约束土体的侵彻过程。基于ANSYS/LS-DYNA显式动力分析和LS-PrePost后处理软件,形象再现了钢质弹丸冲击作用下土体飞溅形成漏斗坑和直线通道的物理过程。侵彻深度计算结果与实验数据吻合较好,钢质弹丸的速度和加速度时程曲线图与理论分析一致,体现了FEM/SPH算法的精确性。研究结果表明FEM/SPH算法在模拟侵彻土体方面具有可行性和有效性。     

改进的接触算法及其在光滑粒子流体动力学中的应用

采用以黎曼解描述粒子间相互作用的接触算法,模拟一维情况下包含间断的流场。在弱间断的条件下,采用弱波近似的黎曼解来描述粒子间的相互作用;在强间断的条件下,则引入非迭代黎曼解来描述粒子间的相互作用。并引入Taylor展开思想,提高了接触算法在自由边界的计算精度。该算法无需引入人为粘性,程序结构简洁。数值计算表明,改进后的接触算法能很好地描述强间断,并有效改善了原始接触算法在自由面的计算缺陷。     

水力旋流器湍流场数值模拟

应用CFD软件FLUENT对水力旋流器进行三维数值模拟，并对网格生成技术、湍流模型、离散方法及边界条件等问题进行了探讨。数值模拟所获得的轴截面上的压力场和速度分布规律与前人的研究结果基本吻合，验证了数值模拟的可行性。为进一步研究水力旋流器的结构优化和分离性能预测提供参考。     

粉末燃料冲压发动机内镁粉尘云层流燃烧模型

对粉末燃料冲压发动机预燃室内镁粉尘云燃烧过程进行了研究,建立了镁粉尘云的一维层流预混燃烧模型。研究表明,镁粉尘云层流火焰传播很稳定,燃烧过程中火焰结构基本不变,燃烧区很薄,而预热区厚度约是燃烧区的2-3倍。粉尘云中镁颗粒的蒸发和气相镁与氧气的均相反应是产生火焰的直接原因,也是火焰得以传播的关键。预热区气相温度升高主要靠燃烧区气体的导热和扩散过来的气相镁与氧气反应释放热量,而预热区颗粒相温度升高主要靠气相对其对流传热。分析了各参数对粉尘云燃烧的影响,颗粒相对浓度对粉尘云燃烧的影响比较复杂,在浓度较低的情况下,增大颗粒相对浓度有利于粉尘云快速燃烧;而在浓度较高的情况下,增大颗粒相对浓度则不利于粉尘云快速燃烧。随颗粒粒径的增加,火焰传播速度减小,火焰温度升高,预热区厚度增大。火焰传播速度和火焰温度随粉尘云初温增加线性增长,预热区厚度随粉尘云初温增加抛物线增长。数值模拟与文献中试验结果的变化趋势相一致。     

快速傅立叶变换用于Fokker-Planck-Landau方程的数值求解

利用谱方法和FFT技术对Fokker-Planck-Landau方程进行了数值求解,研究了均匀空间条件下粒子在速度空间的分布函数随时间的演化。数值计算表明,所用计算方法能够很好地满足质量、动量和能量守恒要求,计算速度与有限差分方法相比大大加快。     

箔条云雷达回波的仿真研究

通过建立箔条干扰下的雷达回波模型,根据回波产生的机理,模拟出回波的幅度服从瑞利分布,相位服从均匀分布。根据统计学原理,推导出箔条云雷达回波的功率谱密度,服从高斯分布。结果与实际相符,说明了该仿真方法的正确性。通过仿真比较和分析了当箔条云中含有静态或动态目标时的雷达回波、幅度分布和相位分布。对箔条干扰研究具有重要意义。     

再生冷却通道结构参数对碳氢燃料流量分配影响

以液体火箭发动机再生冷却技术为背景,采用数值模拟方法对Z-type型并联通道流量分配特性进行研究,分析了通道数、通道截面形状及入口汇流结构对流量分配特性的影响。研究结果表明:在冷却通道总流通面积相同的情况下,随着再生冷却系统通道数的增加,流量分配不均匀度系数从2.59%降低至0.5%,相同流通面积下增加通道数可以有效提高流量分配的均匀程度。通道截面形状对流量分配不均匀度系数影响不大,其影响主要表现在通道内部流体速度分布及换热面积上。从换热效果看,梯形截面构型要优于矩形构型。入口汇流结构倾斜角有利于流量的均匀性分布。该研究对于再生冷却通道结构设计具有一定参考价值。     

Fragment spatial distribution of prismatic casing under internal explosive loading

Non-cylindrical casings filled with explosives have undergone rapid development in warhead design and explosion control. The fragment spatial distribution of prismatic casings is more complex than that of traditional cylindrical casings. In this study, numerical and experimental investigations into the fragment spatial distribution of a prismatic casing were conducted. A new numerical method, which adds the Lagrangian marker points to the Eulerian grid, was proposed to track the multi-material interfaces and material dynamic fractures. Physical quantity mappings between the Lagrangian marker points and Eulerian grid were achieved by their topological relationship. Thereafter, the fragment spatial distributions of the prismatic casing with different fragment sizes, fragment shapes, and casing geometries were obtained using the numerical method. Moreover, fragment spatial distribution experiments were conducted on the prismatic casing with different fragment sizes and shapes, and the experimental data were compared with the numerical results. The effects of the fragment and casing geometry on the fragment spatial distributions were determined by analyzing the numerical results and experimental data. Finally, a formula including the casing geometry parameters was fitted to predict the fragment spatial distribution of the prismatic casing under internal explosive loading.     

固体火箭发动机二次喷射控制矢量喷管流场仿真

采用基于Favre平均的三维N-S方程和k-ε湍流模型对固体火箭发动机二次喷射推力矢量喷管复杂干扰内流场进行数值模拟。空间上采用三阶精度差分格式进行求解,时间上采用隐式Jacobi点迭代方法进行迭代推进,直至流场收敛。数值模拟得到矢量喷管二次射流的激波系结构,以及复杂的主/次流干扰流动图像。二次喷射流场包含复杂的涡系结构和波系结构,还存在着边界层与激波的相互干扰、自由剪切层、激波、膨胀波和大尺寸分离。数值模拟还表明,高温燃气射流导致喷射孔附近喷管壁面处的温度相当高,需采取相应的热防护措施。