Deformation,fragmentation and acceleration of a controlled fragmentation charge casing

Two different finite element software, LS-DYNA and Impetus, have been evaluated to test their ability to predict the deformation, fragmentation and acceleration of a controlled fragmentation charge casing. The general-purpose program LS-DYNA was used with a multi-material ALE formulation and a mass-preserving erosion criterion coupled to a Johnson-Cook fracture criterion. In the Impetus simulations, a third order Lagrangian element formulation was used for the casing and a node-splitting element erosion treatment coupled to a Cockcroft-Latham failure criterion was used to describe casing fracture. The high-explosive gases were described by a discrete particle formalism.In order to acquire data to validate our computational tools and constitutive models, a series of experiments have been performed using a laboratory charge with an internal grooved casing. In the test series, the charge geometry was fixed except that the groove depth were varied from very shallow to very deep resulting in different deformation patterns, fracture modes and terminal velocities. Various diagnostic tools captured the different stages of the expansion and fragmentation of the casing. A high-speed framing camera depicted the deformation pattern before fragmentation and was used to determine the moment when the casing failed. Three different complementary techniques were used to follow the acceleration of the fragments; a Photon Doppler velocimetry to determine the initial acceleration of the casing, double exposed radiographs to estimate the fragment velocity after break-up and a high-speed video to determine the terminal velocity of the fragment after leaving the fireball. In addition, the fragments were soft recovered in a set of sawdust pit tests and their final shape and weight were measured. A SEM was used to characterise the fracture surfaces and to determine the modus of fracture (tensile or shear failure).Comparisons to experiments show that both software can predict the change in deformation behaviour when the groove depth increases, from tangential necking for shallow grooves to radial punching for deep groves. Both software could also reasonable well predict the acceleration of the fragments, though both overestimates the terminal velocity for the charge with the deepest grooves.     

A novel method for determination of lethal radius for high-explosive artillery projectiles

A new model has been defined that enables the estimation of the lethal radius (radius of efficiency) of HE (High Explosive) artillery projectiles against human targets. The model is made of several modules: CAD (Computer Aided Design) modeling, fragment mass distribution estimation, fragment initial velocity prediction, fragment trajectory calculation, effective fragment density estimation, and high explosive projectile lethal radius estimation. The results were compared with the experimental results obtained based on tests in the arena used in our country, and the agreement of the results was good. This model can be used in any terminal-ballistics scenario for high explosive projectiles since it is general, para-metric, fast and relatively easy to implement.     

UCAV空面多目标攻击三维轨迹规划技术

研究了单架无人作战飞机(UCAV)攻击多个地面目标的三维轨迹规划问题。首先,将问题形式化为一类特殊的旅行商问题(TSP),即带动力学约束的邻域访问TSP问题(DCTSPN)。其次,针对规划空间维度过高、搜索代价过大的问题,提出了一种基于概率路标图(PRM)的方法。该方法借鉴了基于采样的运动规划方法的思想,并结合多种组合优化技术,将原本连续状态空间中的轨迹规划问题转化为离散拓扑图上的路由问题。求解过程分为离线预处理和在线查询两个阶段。离线阶段采用Halton拟随机采样算法及Noon-Bean转换方法,将原问题转化为经典的非对称旅行商问题(ATSP);在线阶段根据战场态势的实时变化,快速更新路标图,然后采用LKH算法在线求解问题的近似最优解。为了保证生成的飞行轨迹满足平台的运动学/动力学约束,算法基于Gauss伪谱法构建了局部轨迹规划器。最后,以攻击时间最短为优化指标对算法进行了仿真实验。结果表明,本文提出的方法能够以较高的精度和在线收敛速度生成真实可行的、较优的多目标攻击轨迹。     

IMM-EKF的弹道导弹轨道实时动态预测

快速准确的导弹轨道预测是有效进行导弹防御的前提。提出一种基于IMM-EKF的弹道导弹轨道实时动态预测方法,对导弹发射点和落点分别进行动态预测研究。首先,构建不包含导弹先验信息的IMM-EKF弹道估计器;其次,定义一种模型概率累积因子,用于估计导弹的关机时刻及相应的运动状态;最后,基于导弹运动状态的实时估计,动态预测落点和发射点。仿真实验结果表明:该方法能实时动态预测导弹发射点和落点;发射点预测应利用主动段状态估计,而落点预测则在导弹自由段早期进行为宜。     

榴弹攻击下铝质结构抗爆能力的数值评估

在静爆试验的基础上,结合有限元数值方法对小口径(30 mm)燃烧榴弹舷侧接触爆炸攻击下铝质艇体结构的毁伤进行评估,着重分析了铝质艇体结构在小口径武器攻击下,弹体的攻击角度、弹体的初始侵彻速度以及爆炸产物等对艇体结构破损和毁伤范围的影响以及榴弹爆炸的破片效应.     

Reaction degree of composition B explosive with multi-layered compound structure protection subjected to detonation loading

The explosive reaction degree and protection from explosions are concerns in the military field.In this work,the reaction degree of the composition B explosive was investigated experimentally.Multi-layered compound structures were used as barriers to weaken the blast loads.A comprehensive experiment using a high-speed camera and image processing techniques,side witness plates,and bottom witness plates was presented.Using the experimental fragment velocities,fragment piercing patterns,and damage characteristics,the reaction degree of the explosive impeded by different multi-layered com-pound structures could be precisely differentiated.Reaction parameters of the explosive obstructed by compound structures were obtained by theoretical analysis and numerical simulations.Unlike the common method in which the explosive reaction degree is only distinguished based on the initial pressure amplitude transmitted into the explosive,a following shock wave reflected from the side steel casing was also considered.Different detonation growth paths in the explosive formed.Therefore,all these shock wave propagation characteristics must be considered to analyze the explosive response impeded by compound structures.     

Enhancing the explosive characteristics of a Semtex explosive by involving admixtures of BCHMX and HMX

A well-known ternary plastic explosive, Czech Semtex 1H, contains a mixture of PETN and RDX softened by SBR. In this work, BCHMX was used to replace PETN in Semtex 1H to form Sem-BC+RDX. In addition, another mixture based on BCHMX and HMX as energetic fillers bonded by the polymeric matrix of Semtex 1H (Sem-BC+HMX) was studied. The particle size distribution of each individual explosive was determined to obtain the optimum mixing conditions. Friction and impact sensitivities were determined. The velocity of detonation was reported practically and the detonation properties were calculated by EXPLO5 code. The explosive strength of each sample was measured by the ballistic mortar test. The conclusion confirms that the velocity of detonation of Sem-BC+HMX was the highest in comparison with the prepared samples. Sem-BC+RDX has the least impact and frictions sensitivities. Sem-BC+RDX has higher detonation velocity, detonation properties and explosive strength than Semtex 1H. Addition of BCHMX in Semtex 1H as a replacement for PETN is the candidate to produce a high performance advanced Czech plastic explosive.     

临近空间高超声速目标拦截弹弹道规划

利用弹道规划设计了针对临近空间高超声速飞行器的拦截弹道。分析了临近空间高超声速目标拦截问题,将其定性为临近空间的远程高超声速拦截,并提出弹道规划需求;设计了一种两级助推的拦截弹,建立了考虑地球曲率和自转的拦截弹质点平面运动模型;根据弹道规划需求设计弹道约束,以末速最大、与终点距离误差最小和全程热量最小为指标建立拦截弹弹道规划问题;采用粒子群算法求解弹道,结果表明:符合约束的规划弹道是高抛再入形式,与比例导引弹道和准最佳弹道相比,拦截弹大部分时间飞行在大气层外,有效降低了气动热效应影响和对弹体材料的性能需求,且为末制导段提供良好的初始工作环境。     

Numerical study on the case effect of a bomb air explosion

When considering the bomb explosion damage effect, the air shock wave and high-speed fragments of the bomb case are two major threats. In experiments, the air shock wave was studied by the bare ex-plosives superseding the real cased bomb; in contrast, the bomb case influence was ignored to reduce risk. The air explosion simulations of the MK84 warhead with and without the case were conducted. The numerical simulation results showed that the bomb case significantly influenced the shock wave generated by the bomb: the spatial distribution of shock wave in the near field changed, and the peak value of shock wave was reduced. Breakage of the case and kinetic energy of the fragmentation consumed 3 and 38% of the explosion energy, respectively. The increasing factors of the peak over-pressure induced by the bare explosive on the ground and in the air were 1.43-3.04 and 1.37-1.57, respectively. Four typical stages of case breakage were defined. The mass distribution of the fragments follows the Mott distribution. The initial velocity distribution of the fragments agreed well with the Gurney equation.     

基于车辆稳态动力学特性的汽车动态轨迹规划

针对高动态环境下驾驶辅助系统（ADAS）的轨迹规划问题,提出基于车辆稳态动力学特性的动态轨迹规划算法。该算法首先在交通车轨迹预估基础上建立搜索空间,接着利用车辆稳态动力学模型作为轨迹发生器,在搜索空间中对轨迹发生器产生的轨迹进行评价,最后根据驾驶意图和优化条件选择最优轨迹。文章采用基于最优加速度预瞄理论的轨迹跟踪算法,建立了七自由度车辆动力学模型,并在Simulink环境下搭建仿真平台,进行超车仿真实验。实验结果表明,本文提出的算法可行、有效。     

Residual Characteristic of High-Velocity Steel Fragments Penetrating UHMWPE FRP Laminates

The mechanical performance of ultra-high molecular weight polyethylene fiber (UHMWPE) and its composites were proposed. Penetrated properties of different thicknesses UHMWPE FRP laminates (URP) impacted by 3.3g cubic high velocity fragments were studied. According to the ballistic experimental results and theoretical analysis, the linear relation between ballistic limit vBL and area density AD was confirmed. The relative parameters of showing experientially residual velocity vr were expressed by the function of AD. In the end, versatile experiential expression between vr and AD was found. Prediction of vr and vBL using obtained expressions under the above stated condition of impacting URP was consistent with the experimentaled results. Consequently, the two experiential relations can be used to predict the residual velocity and ballistic limit of cubic high velocity fragments impacting URP. The residual characteristic of high-velocity steel fragments penetrating UHMWPE FRP laminates can be more exactly forecasted by the two derived experiential formulas.     

杀伤战斗部破片定向飞散特性研究

研究了钢预制破片壳体在四种装药爆轰驱动下的飞散特性。求得了使破片飞行路径呈平行型的预制破片壳体的临界半径及相应的破片初始速度。     

非奇异快速终端滑模及动态面控制的轨迹跟踪制导律

针对飞行器跟踪预设轨迹的问题,提出非奇异快速终端滑模和角度约束的轨迹跟踪制导律。通过引入虚拟目标点,提出参考轨迹曲率半径的期望视线角约束条件,建立带有视线角约束并考虑自动驾驶仪动态特性的轨迹跟踪数学模型。为了保证在有限时间内跟踪预设轨迹并避免出现奇异问题,采用快速非奇异终端滑模和动态面控制方法进行制导律设计。推导出视线角误差和轨迹跟踪误差之间的数学关系,并利用Lyapunov稳定性准则证明轨迹跟踪误差最终有界任意小。与弹道成型轨迹跟踪制导律进行仿真对比,仿真结果表明所提出的制导律具有良好的跟踪性能及鲁棒性。     

探月飞船初次再入段纵程的解析预测方法

探月飞船升阻比较低,为实现长纵程飞行,必须采用跳跃式再入方式。在跳跃式再入轨迹在线规划或预测制导中,如何快速准确地预测初次再入段纵程是一个非常关键的问题。针对这一问题,研究提出一种解析预测方法:利用匹配渐进展开方法得到再入纵向运动方程的闭型近似解;将初次再入段轨迹分为三段,第一段采用高度作为积分自变量,并利用复合梯形公式得到纵程,第二段和第三段分别采用二次多项式来拟合阻力加速度-能量剖面,根据近似解结果反解出多项式系数,并将得到的阻力加速度倒数-能量函数进行积分,得到第二段和第三段的纵程;对解析预测方法的精度和计算效率进行分析,结果表明该方法计算精度较高,速度快,可用于跳跃式再入轨迹的在线规划和制导。     

Performances and direct writing of CL-20 based ultraviolet curing explosive ink

A new type of explosive ink formulation that can be quickly cured was prepared with unsaturated polyester as binder,styrene as active monomer,2,4,6-trimethylbenzoyl-diphenylphosphine oxide as photoinitiator,and hexanitrohexaazaisowurtzitane (CL-20) as the main explosive.Then the explosive ink direct writing technology was used to charge the micro-sized energetic devices,the curing mechanism of the explosive ink was discussed,and the microstructure,safety performance and explosive transfer performance of the explosive ink molded samples were tested and analyzed.Results indicate that the composite material has a fast curing molding speed,its hardness can reach 2H within 8 min.The crystal form of CL-20 in the molded sample is still ε type.The CL-20 based W-curing explosive ink formulation has good compatibility,its apparent activation energy is increased by about 3.5 kJ/mol.The composite presents a significant reduction in impact sensitivity and its characteristic drop height can reach 39.8 cm,which is about 3 times higher than the raw material.When the line width of charge is 1.0 mm,the critical thickness of the explosion can reach 0.015 mm,and the explosion velocity is 7129 m/s when the charge density is 1.612 g/cm3.     

Liquid phase in-situ synthesis of LiF coated boron powder composite and performance study

Among practical metal additives, boron (B) has a high volumetric heating value, making it a promising choice as a fuel additive. Although B can theoretically yield a large amount of energy upon complete combustion, its combustion is retarded by the initial presence of B oxide, which coats the surfaces of B particle. To improve the ignition and combustion properties of B powder, LiOH and NH₄F were used as precursors to synthesize uniformly LiF-coated B composites (LiF-B) in situ. The LiF-B mixture was also prepared for comparison using a physical method. X-ray diffraction (XRD), Fourier-transform infrared (FTIR), scanning electron microscope (SEM), and energy-dispersive X-ray spectroscopy (EDS) were used to characterize the morphologies and compositions of the products. The thermal and combustion properties of the samples were characterized by thermal gravity-differential thermal gravity (TG-DTG), differential scanning calorimetry (DSC) and closed bomb experiment. The XRD, FTIR, SEM and EDS results demonstrated the successful preparation of the coated LiF-B sample. The TG-DTG and closed bomb experiment results indicated that the addition of LiF decreased the ignition temperature of B powder, and increasing its reaction efficiency. DSC results show that when LiF-B was added, the released heat of underwater explosive increased by 6727.2, 7280.4 and 3109.6 J/g at heating rates of 5, 10, and 15 °C/min, respectively. Moreover, LiF-B decreased the activation energy of secondary combustion reaction of explosive system as calculated through Kissinger's method by 28.9%, which indicated an excellent catalytic effect for the thermal decomposition of underwater explosive. The results reveal that LiF can improve the combustion efficiency of B powder, thereby increasing the total energy of explosives. The mechanical sensitivity increased slightly after adding LiF-B to the underwater explosive. Compared to the underwater explosive with added B, the mechanical sensitivity of the explosive with added LiF-B was significantly lower.     

Sensitivity Test and Data Analysis for Storage Reliability Assessment of Explosive Initiator

The explosive initiator is one kind of sensitivity products with long life and high reliability.In order to improve the storage reliability assessment,the method of storage reliability assessment for explosive initiator was proposed based on time series model using the sensitivity test data.In the method,the up and down test was used to estimate the distribution parameters of threshold.And an approach to design the up and down test was present to draw better estimations.Furthermore,the method of shrinkage estimation was introduced to get a better estimation of scale parameter by combining the sample information with prior information.The simulation result shows that the shrinkage estimation is better than traditional estimation under certain conditions.With the distribution parameters estimations,the time series models were used to describe the changing trends of distribution parameters along with storage time.Then for a fixed storage time,the distribution parameters were predicted based on the models.Finally,the confidence interval of storage reliability was obtained by fiducial inference.The illustrative example shows that the method is available for storage reliability assessment of the explosive initiator with high reliability.     

舰炮防空战法的可视化技术

舰炮防空是水面舰艇自防御的重要手段之一.文章围绕可视化和虚拟现实技术,分析了舰炮防空战法研究的内容,给出了舰炮防空战法可视化的方法,重点探讨了其中的防空弹道仿真、反舰导弹弹道仿真、海浪水文气象环境仿真、战场焰火和爆炸效果仿真及舰炮防空战法仿真中实时碰撞检测的基本方法和优化技术,提出了解决这些问题的主要方法,并在PC机上进行了实现.最后给出了舰炮防空战法的可视化结果.     

高超声速滑翔飞行器滑翔段初始状态的唯一性和最优性分析

临近空间高超声速滑翔飞行器的弹道特性主要受滑翔段初始状态和飞行器控制律影响。在飞行器控制律确定的情况下,研究了滑翔段初始状态对高超声速滑翔飞行器弹道特性的影响规律。按照滑翔弹道的不同形式,在纵向平衡滑翔条件下,通过理论推导得出飞行器状态变量的解析式,结合平衡滑翔条件分析平衡滑翔弹道滑翔段初始状态的唯一性；在纵向跳跃滑翔条件下,构建弹道性能评价指标,利用群智能算法,寻找弹道性能最优时的滑翔段初始状态。利用单因素敏感性分析方法,分别对两种滑翔弹道的滑翔段初始状态进行敏感性分析,初始状态中初始速度对弹道特性的影响最大。对高超声速滑翔飞行器初始状态唯一性与最优性的分析,可为高超声速滑翔飞行器的弹道设计、弹道跟踪、轨迹预测和轨迹优化提供借鉴。     

充气式防热罩再入轨道设计

充气式防热罩是最近几年才出现的新型航天回收技术 ,为各种轨道航天器的回收提供了简单、可靠、经济的途径 ,有望逐渐取代降落伞在回收系统中的地位。合理选择防热罩外形 ,用CFD方法计算了防热罩以亚、跨、超声速飞行时的阻力特性 ,并与工程估算方法作对比。分别采用单次和二次充气方案 ,完成了 90kg货舱安全着陆的再入轨道设计。在对两种充气减速方案的再入轨道进行比较分析之后指出 :对于回收同样重量的货舱 ,二次充气方案在保证减速效果的前提下大幅度减轻了防热系统设计的负担 ,是理想的ITS再入减速方案