蒙特卡洛剖分投影法的爆震弹破片平均比动能评估

为了评估爆震弹封装壳体破片致死、致伤半径,并削弱其平均比动能,基于蒙特卡洛剖分投影法建立破片平均比动能计算模型。通过LS-DYNA及自编程联合仿真方法对爆震弹自然破片、半预制破片比动能计算模型进行仿真求解,获得了全破片全时域的质量分布、初始速度、垂直靶分布、平均比动能阈值及安全半径等指标。结果表明:在相同装药参数下,半预制破片相比于自然破片,质量分布、初始速度阈值更低且更集中;在小于等于2.5 m范围内半预制破片平均比动能阈值比自然破片更低,大于2.5 m后半预制破片平均比动能阈值比自然破片更高。相比于自然破片,采用半预制破片可显著减小爆震弹破片的致死半径,但并不能明显增大其安全半径。     

Influence of shaped charge structure parameters on the formation of linear explosively formed projectiles

The research of LEFP (linear explosive forming projectile) is of great value to the development of new warhead due to its excellent performance. To further improve the damage ability of the shaped charge warhead, a special shell overhanging structure was designed to increase the charge based on the traditional spherical charge, in which case the crushing energy of LEFP could be guaranteed. LS-DYNA was used to simulate different charge structures obtained by changing the number of detonation points, the length of shell platform, the radius of curvature and the thickness of liner. The RSM (response surface model) between the molding parameters of LEFP and the structural parameters of charge was established. Based on RSM model, the structure of shaped charge was optimized by using multi-objective genetic algorithm. Meanwhile, the formation process of jet was analyzed by pulsed X-ray photography. The results show that the velocity, length-diameter ratio and specific kinetic energy of the LEFP were closely related to the structural parameters of the shaped charge. After the optimization of charge structure, the forming effect and penetration ability of LEPP had been significantly improved. The experimental data of jet velocity and length were consistent with the numerical results, which verifies the reliability of the numerical results.     

双球破片冲击柱壳装药的临界起爆条件

为了研究导弹战斗部(柱壳装药)在破片场中的累积毁伤问题,在现有单破片起爆平板装药的Jacobs-Roslund经验准则的基础上,分别建立了考虑破片尺寸、破片撞击角度、柱壳装药的装药半径和壳体厚度的单球破片、双球破片冲击柱壳装药临界起爆条件的工程分析模型.该模型计算结果与数值模拟结果和现有试验结果相吻,证明利用该模型能较...     

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.     

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

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

圆柱形装药驱动轴向预制破片飞散特性

为了增强柱形战斗部轴向威力,在无壳柱形战斗部底面布置单层离散的预制破片。开展圆柱形TNT装药驱动轴向预制破片飞散试验,获得预制破片的最大初速、飞散角等特征参数;运用LS-DYNA软件对装药驱动预制破片过程进行数值模拟,阐述预制破片群飞散过程;对装药驱动整体平板理论计算公式进行改进,获得预制破片的最大初速。结果表明:破片初速理论计算结果、数值计算结果和试验结果吻合良好;随着与装药底部中心距离的加大,破片初速、径向飞散角分别近似呈“抛物线”减小、增大;试验实测、理论计算得到的破片最大初速值超过2500 m/s,试验实测的径向飞散角最大约为22°,而周向飞散角则普遍较小,均值在5°以内。     

Effect of microstructure on short pulse duration shock initiation of TATB and initial response mechanism

Micro-TATB particles with different sizes and 3D nanoporous TATB architectures with different specific surface areas were prepared through recrystallization to study short pulse duration shock initiation properties by electric gun technology. For micro-TATB, the initiation threshold significantly decreases with TATB average size ranging from 79.7 μm to 0.5 μm. For 3D nanoporous TATB architecture, the initiation threshold decreases and then increases with specific surface areas increased from 9.6 m²/g to 36.2 m²/g. The lowest initiation thresholds are obtained for the micro-TATB with average sizes of 1.3 μm and 0.5 μm, and 3D nanoporous TATB architecture with specific surface area of 22.4 m²/g. The shock initiation thresholds of micro-TATB and 3D nanoporous TATB architectures show significantly decreases with the porosity increased. The decomposition reaction and thermal conductivity properties were further investigated to understand the initial response mechanism. High porosity provides more collapse sites to generate high temperature for formation of hot spots. The low thermal conductivity and decomposition temperature could enhance the formation and ignition of the hot spots, and initial decomposition reaction of TATB. The effect of the decomposition temperature is higher than that of the thermal conductivity on the shock initiation properties. The enhanced decomposition reaction could promote energy release and transfer process from the ignition to the combustion. This work offers a new insight to understand the effects of microstructure on the shock initiation properties and the initial response mechanism of TATB.     

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.     

脱壳弹弹托分离的双目视觉测量

脱壳弹的弹托分离过程对弹体的飞行稳定性和打击效能具有显著的影响。针对脱壳弹出膛后弹托相对弹体的六自由度运动过程,提出采用基于双目视觉原理的弹托分离角测量方法。通过在弹体表面和弹托表面设置标记点,采用图像处理和跟踪算法实现对标记点的识别和跟踪。结合双目空间标定参数解算出标记点的空间位置,进而获得弹托相对弹体的分离角度,同时采用实验验证了上述测量方法的精度,其精度达到2%。以实验室进行的初速度1 550 m/s和1 750 m/s的脱壳弹射击试验为例,测量并分析了不同初速度下的弹托分离轨迹。结果表明:脱壳弹分离初速度越快,弹托分离轨迹越靠近弹体。     

一种新的动态信道化接收机设计方法

针对复杂电磁环境中存在多个非均匀分布的、不同带宽的信号的情况,提出一种新的动态信道化接收机设计方法。设计一个半带滤波器,并通过内插运算压缩它的频率响应,形成周期频谱;对输入信号进行滤波,同时使原信号经过一定延迟后减去滤波后的信号,形成两路在频谱上互补的信号;用均匀分析滤波器组对得到的两路信号分别进行滤波分解处理,实现监视频段的均匀信道化;根据能量检测结果将属于同一信号的相邻子信道信号作为下一步重构处理的输入,并设计出相应的综合滤波器组提取信号。理论分析和仿真结果验证了新方法的有效性。     

圆柱壳受不对称径向冲击时的塑性动力屈曲分析

本文发展了Lindberg、Duffey等人的思想,考虑不对称冲击下圆柱壳的塑性动力屈曲问题,旨在研究内凹变形与屈曲的关系。用能量法导出了不对称基本变形的微分方程,近似地将之与对称情形下屈曲扰动控制方程联合,构成一个初值问题,然后进行数值求解。考虑了壳厚变化,放弃了屈曲过程中速度保持不变的假定(对称情形中采用)。本文的研究结果很好地解释了高速冲击时倒塌模式上无屈曲波形,屈曲波形遍布加载范围内且并不明显减弱,局部冲击临界速度高于对称情况形下临界速度等许多实验与已有的分析结果,从而明确了屈曲在塑性动响应中的地位。     

On Growing Operating Costs in the Armed Forces

Operating costs per unit of input or output are often claimed to grow faster in defence than elsewhere in the economy. In this paper, we outline several possible reasons as to why operating costs increase at higher rates in defence, including differences in the input factor mix and productivity growth, the technological complexity required to maintain the relative effect of weapon systems and a gradual reduction in the number of units and activity. We investigate whether operating costs grow at a faster rate than elsewhere in the economy, by estimating growth in real output unit costs in the Norwegian Armed Forces for the period 1994–2013, using activity as a measure of output. We find no increase beyond general inflation in structural (fixed) costs, whereas activity-based (variable) costs per unit of activity increase significantly.     

Parametric study of single confined fragment launch explosive device

In this article, parametric study of single confined fragment launch device was carried out. The configuration proposed was further studied to derive the empirical relationship for effect of fragment size, charge size, confinement thickness on fragment velocity. The simulations were carried out using ANSYS-AUTODYNE explicit solver. Fragment velocities were estimated as a function of different parametric combinations of explosive quantities, charge length to diameter ratio, fragment height to diameter ratio, confinement thickness, fragment material and fragment mass. The data was further converted to charge to metal ratio under fragment and confinement. It was observed that, increase in confinement thickness, charge quantity and decrease in fragment height increases the fragment velocity. It is also noted that, charge to metal mass ratio under fragment significantly affects the fragment velocity. At the end, an empirical relationship for fragment velocity interms of all these parameters was established. Using these relations, two velocities 1831.92 m/s and 2523.9 m/s required for NATO STANAG 4496 IM test were estimated. The design parameters for these velocities are presented. Also, the results estimated using the empirical relationship has been compared with published experimental data. Error in the predicted velocities is within the acceptable range. The empirical relationship proposed will be useful for finalization of design of the fragment launch device.     

纵桁对环肋圆柱壳体水下振动与声辐射的影响

采用结构有限元耦合流体边界元的附加质量附加阻尼算法,对具有不同数量纵桁的加筋环肋圆柱壳水下振动与声辐射噪声进行了数值计算,并对数值计算结果进行了初步的比较和分析,讨论了纵桁的数量对圆柱壳水下振动与声辐射的影响.结果表明:纵桁对圆柱壳水下的减振降噪有较好的效果,并且纵桁的数量越多,减振降噪效果越好,但增加纵桁的数量也相应增加了艇体的重量,在潜艇结构设计时须在降噪及艇体重量间综合权衡.     

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.     

铝质壳体弹丸破碎性试验破片回收率研究

针对当前没有铝质壳体弹丸破碎性试验破片回收率标准的现状,进行某型铝质壳体破甲弹破碎性试验,认为破片回收率低的原因是回收过程中破片损失和部分金属铝熔化,并从理论上计算铝质壳体在爆炸破裂过程熔化率,验证试验中破片回收率的合理性。     

Study on energy release characteristics of reactive material casings under explosive loading

Reactive Materials (RMs), a new material with structural and energy release characteristics under shock-induced chemical reactions, are promising in extensive applications in national defense and military fields. They can increase the lethality of warheads due to their dual functionality. This paper focuses on the energy release characteristics of RM casings prepared by alloy melting and casting process under explosive loading. Explosion experiments of RM and conventional 2A12 aluminum alloy casings were conducted in free field to capture the explosive fireballs, temperature distribution, peak overpressure of the air shock wave and the fracture morphology of fragments of reactive material (RM) warhead casings by using high-speed camera, infrared thermal imager temperature and peak overpressure testing and scanning electron microscope. Results showed that an increase of both the fireball temperature and air shock wave were observed in all RM casings compared to conventional 2A12 aluminum ally casings. The RM casings can improve the peak overpressure of the air shock wave under explosion loading, though the results are different with different charge ratios. According to the energy release characteristics of the RM, increasing the thickness of RM casings will increase the peak overpressure of the near-field air shock wave, while reducing the thickness will increase the peak overpressure of the far-field air shock wave.     

Fracture mechanism of a cylindrical shell cut by circumferential detonation collision

The failure mechanism of a cylindrical shell cut into fragments by circumferential detonation collision was experimentally and numerically investigated. A self-designed detonation wave regulator was used to control the detonation and cut the shell. It was found that the self-designed regulator controlled the fragment shape. The macrostructure and micro-characteristics of fragments revealed that shear fracture was a prior mechanism, the shell fractured not only at the position of detonation collision, but the crack also penetrated the shell at the first contact position of the Chapmen-Jouguet (C-J) wave. The effects of groove number and outer layer thickness on the fracture behavior were tested by simulations. When the thickness of the outer layer was 5–18 mm, it has little effect on fragmentation of the shell, and shells all fractured at similar positions. The increase of the groove number reduced the fracture possibility of the first contact position of the C-J wave. When the groove number reached 7 with a 10 mm outer layer (1/4 model), the fracture only occurred at the position of detonation collision and the fragment width rebounded.     

Field test and probabilistic analysis of irregular steel debris casualty risks from a person-borne improvised explosive device

Person-borne improvised explosive devices (PBIEDs) are often used in terrorist attacks in Western countries. This study aims to predict the trajectories of PBIED fragments and the subsequent safety risks for people exposed to this hazard. An explosive field test with a typical PBIED composed of a plastic explosive charge and steel nut enhancements was performed to record initial fragment behaviour, including positions, velocity, and trajectory angles. These data were used to predict the full trajectory of PBIED fragments using a probabilistic analysis. In the probabilistic analyses a probability of fatality or serious injury was computed. Based on the results presented, many practical conclusions can be drawn, for instance, regarding safe evacuation distances if a person were exposed to a suspected PBIED.     

基于对偶四元数的航天器六自由度相对运动输入有界控制

利用对偶四元数的理论来分析航天器六自由度的相对运动,设计了一种考虑输入有界的姿轨一体化控制器。在介绍对偶数和对偶四元数的基础上推导六自由度相对运动的姿轨耦合模型;利用双曲正切函数绝对值小于1的特性来显式地构造有界控制器,分别设计两个相互耦合的自适应调节律来动态地改变控制器的输出,并基于李雅普诺夫稳定性理论严格证明了闭环系统的全局渐近稳定性;利用数学仿真实验来验证该控制器和控制力满足给定的约束条件,能够实现航天器六自由度相对运动的精确稳定控制,并且对模型参数不确定性和外界扰动具有鲁棒性。与其他方法相比,由该控制器计算得到的控制力矩器的收敛速度更快,输出的控制力矩和控制力的最大幅值更小,且消耗的能量也更少。