高强低合金钢甲弧形弹侵彻速度与深度的计算

给出了弧形弹头侵入３０ＭｎＣｒＮｉＭｏＢ高强低合金钢过程中弹头的侵彻深度与速度的关系式。阐明了弹头初速度只有超过某临界值才能侵入钢甲。指出弹头侵入钢甲一定深度后，空穴膨胀理论才能适用，此时弹头与钢甲之间的摩擦系数取０．０２更符合实际。     

Lagrangian description of transport equations for shock waves in three-dimensional elastic solids

A set of transport equations for the growth or decay of theamplitudes of shock waves along an arbitrary propagation directionin three-dimensional nonlinear elastic solids is derived using theLagrangian coordinates.The transport equations obtained showthat the time derivative of the amplitude of a shock wave alongany propagation ray depends on (i) an unknown quantity immediatelybehind the shock wave,(ii) the two principal curvatures of theshock surface,(iii) the gradient taken on the shock surface ofthe normal shock wave speed and (iv) the inhomogeneous term.whichis related to the motion ahead of the shock surface.vanisheswhen the motion ahead of the shock surface is uniform.Severalchoices of the propagation vector are given for which the tran-sport equations can be simplified.Some universal relations,which relate the time derivatives of various jump quantities toeach other but which do not depend on the constitutive equationsof the material,are also presented.     

球形弹丸正撞击薄板防护屏碎片云特性研究

以质量、动量和能量守恒为基础,结合平面激波理论和热力学理论,对球形铝弹丸正撞击薄板铝防护屏形成的碎片云特性进行了建模,模型计算结果与实验结果吻合较好。利用该模型对多种工况下碎片云特性进行了计算,结果表明:(1)碎片云质心速度和膨胀速度随撞击速度和弹丸直径的增加而增大,随防护屏厚度增大而减小;膨胀半角随撞击速度和防护屏厚度的增大而增大,随弹丸直径的增大而减小;(2)速度和膨胀半角变化曲线具有相似性;(3)对于给定的弹丸和防护屏材料,碎片云中受到激波加载部分的材料各相态质量分数只与弹丸撞击速度有关。这些规律与实验规律吻合。     

考虑头形磨损变化的动能弹侵彻深度研究

质量损失主要发生在弹体头部,讨论残余弹头形状因子的主要影响因素,给出了更为一般化的残余弹头形状因子与初始撞击速度之间的关系. 假设在整个侵彻过程中弹体头部形状因子及质量保持瞬时速度平方的线性关系,得到了能够考虑弹体头部质量磨损的侵彻深度计算方法. 弹体质量损失将导致其侵彻深度存在上限,该修正模型可以找出达到侵彻深度上限的弹体初始撞击速度.      

弹体垂直侵彻深度工程计算模型

利用空腔膨胀理论计算各种形状弹头的法向侵彻阻力，并考虑作用在弹体头部的库仑摩擦阻力和作用在弹体侧壁粘滞摩擦阻力的影响，根据牛顿定律确定弹体在靶体介质中的运动微分方程，由运动边界条件确定弹体的侵入深度，最后与试验结果进行了比较，表明结果是可靠的。     

穿甲机理的数值研究

为了深入开展穿甲机理的研究,本文对三组穿甲问题——钨弹以1000m/s打钢靶,钢弹以1000m/s打钢靶,钢弹以1000m/s打铝靶,作了数值研究,应用流体弹塑性力学模型,二维不定常欧拉型编码CPG-LTDL code进行计算,计算结果与试验对比符合良好。 对计算结果的分析指出:(1)弹丸侵彻能力不是取决于它的整体动能,而是取决于弹靶接触面上的压力,当靶板参数和弹材固定时,可用弹丸初始面积比动能来标志,(2)铝靶在抗侵彻过程中呈流体性态,高强度铝钣的抗侵彻性能不如普通的A3钢钣,因此应致力于研制具有抗高压、耐高温的新型铝材,而不是单纯追求提高铝材的常规强度。     

Numerical simulation of flow fields induced by a supersonic projectile moving in tubes

Numerical study of the shocked flows generated by a supersonic projectile released from a launch tube into a big chamber has been performed in this paper. Based on fixed Cartesian grids, the two-dimensional axisymmetric Euler equations are solved by the fifth-order WENO scheme implemented with moving boundary conditions. Both the level set technique and ghost fluid method are used for capturing the moving interface of the projectile implicitly. The numerical results show that complex shock phenomena exist in the transient shock flow, resulting from shock-wave reflection, shock-wave focusing, shock-wave/projectile interaction and shock-wave/contact surface interactions. The relationships between the acceleration of the projectile and the transient shock flow are also discussed in detail.     

弹体高速入水特性实验研究

进行了速度在35～160 m/s的平头、卵形和截卵形弹体入水实验,利用高速相机记录了弹体入水和空泡扩展的详细过程,得到了3种弹体在入水初期的水中弹道轨迹和空泡形状,并比较分析了弹体头部形状对入水弹道稳定性的影响.结果表明,平头弹体在水中飞行有良好的弹道稳定性,截卵形弹体往往由于受力不均衡在入水后期发生偏转,卵形弹体则在...     

SIMILARITY LAWS AND SCALING OF PENETRATION TESTS

利用球形空腔膨胀理论的响应函数对靶板在弹体撞击下的侵彻进行分析，并基于量纲理论建立了侵彻相似律的一般形式.通过对侵彻混凝土靶的已有实验数据进行分析，得到了不同缩比因子下，无量纲侵彻深度曲线随侵彻速度不同而发生偏转和偏移.考虑到材料的应变率效应和损伤演化，认为几何相似律对混凝土等脆性材料在抗侵彻时不成立，有必要引入材料应变率效应和损伤演化对几何相似律的修正.     

A new,high-resolution shock-capturing hybrid scheme of flux vector splitting-Harten's TVD

The paper presents a new high-resolution hybrid scheme combining implicit flux vector splitting with Harten's TVD, which is proved suitable for shock-capturing calculation in gasdynamics. Fluxsplitting procedures are applied to discretize the implicit part of the Euler equations whereas Harten's numerical fluxes are used to calculate the residual of steady-state solutions. It ensures good shock-capturing properties and produces sharp numerical discontinuities without oscillations. It excludes expansion shocks and leads only to physically relevant solutions. The block-line-Gauss-Seidel relaxation procedure (block-LGS) is used to solve the resulting difference equations. The time step and the CFL number are much larger than those in the linearized block-alternating-direction-implicit approximate factorization method (block-ADI). Numerical experiments suggest that the hybrid scheme not only has a fairly rapid convergence rate, but also can generate a highly resolved approximation to the steady-state solution. Hence scheme seems to lead to an effective nonoscillatory shock capturing method for steady transonic flow. Project Supported by National Natural Science Foundation of China     

高速侵彻混凝土弹体的形状演化模拟

高速侵彻时, 弹靶间摩擦力所做功全部转变为热并作用于弹体, 假设弹体质量损失来源于弹体表面材料熔化脱落, 可以预测弹体总质量损失. 结合分析得到的弹靶间摩擦功与压力轴向功关系, 通过相应的弹体质量损失分析, 获得侵彻过程中弹体表面各点的回退位移(垂直于弹体轴心线方向), 实现弹体形状演化的模拟. 弹体形状、侵彻深度及质量损失的计算结果均与实验结果基本吻合.      

细长薄壁弹体的屈曲和靶体等效分析

利用金属靶开展动能深侵彻弹的穿甲屈曲实验研究。实验观察到长/短型弹体动塑性屈曲破坏分别表现为轴向皱褶型和轴向外翻撕裂型2种基本模式。不同的屈曲破坏模式与弹体几何、撞击初条件以及靶材等密切相关。区别于刚性尖头弹穿甲金属靶的韧性隧道开孔,尖头弹因屈曲破坏易变形为钝头形,导致弹体穿透靶板表现为挤凿穿甲。利用薄壁中空柱壳的弹性欧拉屈曲分析和弹体塑性屈服的极限分析给出弹体动塑性屈曲的临界条件。同时给出屈曲实验中混凝土靶和金属靶的等效条件。     

Ti-6Al-4V弹体破坏模式对冲击反应的影响研究

Ti-6Al-4V材料是武器结构轻量化时的重要替代材料，其冲击反应将可能增加战斗部毁伤威力，但目前缺乏对其冲击反应条件及反应机理的研究。本文将采用试验与理论分析方法，研究结构破坏模式对Ti-6Al-4V材料冲击反应的影响，获得其冲击反应条件及反应机理。设计并开展了钛合金弹（头部与壳体均为钛合金）与复合弹（头部碳/碳复合材料、壳体空心钛合金圆柱）正侵彻混凝土试验，撞击速度在222~1008 m/s之间。钛合金弹激发了剧烈的氧化冲击反应，但复合弹未产生冲击反应。破坏模式宏细观分析显示，钛合金弹侵彻后宏观结构基本完整，仅表面发生摩擦磨损，以细观组织剪切变形为主要失效模式，形成尺寸在微米量级至百微米量级的颗粒碎片，碎片个数可高达3×106。复合弹的钛合金空心圆柱被撕裂成块，撕裂面沿剪切带方向发展，碎块尺寸在毫米或以上量级，个数至多百余个。碎片供氧和供热的效率均与碎片尺寸成反比，而特定供氧与供热条件下，碎片尺寸足够小是Ti-6Al-4V材料发生冲击反应的必要条件，这是钛合金弹发生冲击反应而钛合金空心圆柱无法激发冲击反应的本质原因。在具备冲击反应必要条件的前提下，碎片个数越多，冲击反应烈度越高。     

基于局部相互作用理论的侵彻弹头部形状优化及仿真

以局部相互作用理论为基础,引入与弹体头部形状相关的开坑计算方法和归一化弹体头部形状方程,给出了任意头部形状弹体侵彻混凝土深度的计算模型。利用最大侵深法,得到了无量纲头部形状控制参数表达式及经典变分头部形状优化设计方法。理论计算及弹靶分离仿真模拟计算结果与实验结果吻合较好。研究结果表明:弹体头部相对半径较小时,球头锥形和球头卵形弹体优化后得到的头部形状分别为尖头锥形和尖头卵形;优化截头弹体的侵彻深度大于优化尖头弹体,而优化截锥形弹体的侵彻深度最大;弹体头部形状对弹体侵彻过载的影响显著,优化弹体头部形状可以有效地提高侵彻深度。     

Focusing shock waves generated by an accelerating projectile

Unsteady flow fields past a projectile either in parabolic motion or in accelerated motion are investigated numerically, with special attention to shock wave focusing phenomena induced by the motions. The two-dimensional Euler equations are solved by a moving-cell, MUSCL TVD finite-volume method with Steger-Warming flux-vector splitting. The pattern of shock wave focusing induced by the projectile is found to be similar to that of shock waves reflected from a concave wall observed in a laboratory experiment using a shock tube. Effects of projectile shapes and trajectories on focusing are also examined.     

Three dimensional shock wave/boundary layer interactions

An investigation into a three-dimensional, curved shock wave interacting with a three-dimensional, curved boundary layer on a slender body is presented. Three different nose profiles mounted on a cylindrical body were tested in a supersonic wind tunnel and numerically simulated by solving the Navier–Stokes equations. The conical and hemispherical nose profiles tested were found to generate shock waves of sufficient strength to separate the boundary layer on the cylinder, while the shock wave generated by the ogival profile did not separate the boundary layer. For the separated flow, separation was found to occur predominantly on the windward side of the cylinder with the lee-side remaining shielded from the direct impact of the incident shock wave. A thickening of the boundary layer on the lee-side of all the profiles was observed, and in the conical and hemispherical cases this leads to the re-formation of the incident shock wave some distance away from the surface of the cylinder. A complex reflection pattern off the shock wave/boundary layer interaction (SWBLI) was also identified for the separated flow cases. For comparative purposes, an inviscid simulation was performed using the hemispherical profile. Significant differences between the viscous and inviscid results were noted including the absence of a boundary layer leading to a simplified shock wave reflection pattern forming. The behaviour of the incident shock wave on the lee-side of the cylinder was also affected with the shock wave amalgamating on the surface of the cylinder instead of away from the surface as per the viscous case. Test data from the wind tunnel identified two separation lines present on the cylindrical surface of the hemispherical SWBLI generator. The pair of lines were not explicitly evident in the original CFD simulations run, but were later identified in a high-resolution simulation.     

Extension,torsion and expansion of an incompressible,hemitropic Cosserat circular cylinder

Constitutive equations for the stress and couple stres on an incompressible, hemitropic, constrained Cosserat material are derived, and the theory is applied to study the problem of finite extension, torsion and expansion of a circular cylinder. As in the theory of isotropic simple elastic materials, it is found that the deformation is controllable by application of only a normal force and a tosional moment at the cylinder ends. It is shown that in general the well known universal relation between the torsional stiffness and the axial force for incompressible, isotropic simple materials in the limit as the twist goes to zero does not exist for incompressible, hemitropic Cosserat materials. However, for a special and unusual class of hemitropic materials, the same universal formula is found to hold for a certain reduced torsional stiffness. The main problem is solved completely for incompressible, hemitropic, linearly elastic, Cosserat materials; and certain additional special features of the Kelvin-Poynting type, which here appear to the first order in the amount of twist of the cylinder, are derived and discussed in relation to experimentally observed composite material behavior.     

Modeling on mass abrasion of kinetic energy penetrator

An engineering model on mass abrasion of kinetic energy penetrator is presented to predict the nose shape and mass loss of the residual projectile after high-speed penetration into concrete. The experimental analysis indicates that the kinetic energy of penetrator (i.e., mass and velocity of projectile) and the hardness of aggregate of concrete significantly affect the mass abrasion of projectile. A theoretical upper limit exists for the mass loss. More general relationship between mass loss and impact function I of projectile is constructed. Graphical discussion declares that the most mass loss occurs on the nose of the projectile and the eroding nose approaches to an ogival shape with a smaller value of caliber-radius-head (CRH). A relative rate of mass abrasion on ogive-nose is further defined and analyzed. The mass loss from abrasion on kinetic energy (KE) penetrator may be evaluated through the variation of nose shape.     

Simulation of shape variation of projectile nose during high-speed penetration into concrete

A formula is developed to estimate the total mass loss of projectile, based on the assumptions that the peeling of molten surface layer in projectile nose is the primary cause of mass loss, and the frictional heat is totally absorbed by the projectile. Extrapolating this formula to predict the mass loss of local area of projectile, the receding displacement on projectile surface is obtained, which is vertical to the symmetry axis of projectile. Thereby, a finite difference method model is constructed to simulate the variation of projectile shape. The shape of residual projectile, depth of penetration of projectile and its mass loss obtained by calculation are found in good consistency with respective experimental data.     

PELE贯穿薄靶后外壳破片径向速度计算方法与影响因素分析

基于横向效应增强型弹丸(PELE)侵彻金属薄靶板过程分析,将弹体前端在撞击作用下的变形过程分解为轴向一维压缩和径向自由膨胀两个变形阶段;依据冲击波理论,给出了弹体前端的冲击波压缩势能,由功能转化原理,给出了PELE前端外壳在靶后形成破片的最大径向飞散速度计算公式。计算结果在多种工况下均与文献的实验结果较为一致。计算结果表明:PELE靶后外壳破片的最大径向飞散速度与外壳和内芯材料的体积模量和泊松比有关,且随二者的增大而增大;PELE外壳破片的最大径向飞散速度是壳体和内芯在冲击波压缩作用下共同径向膨胀的结果,且外壳膨胀能在弹体整体膨胀能中所占比例较大,计算中应当同时考虑弹体外壳和内芯材料的横向膨胀效应对弹体破片径向飞散速度的影响。