Effects of Truncated Ovate Nose Diameter of the Penetration Warhead on the Ballistic Deflection
-
摘要: 针对截卵形头部弹体斜侵彻靶标时弹道发生初始偏转的问题,建立了分析截卵平台直径对初始弹道偏转影响的理论和数值仿真计算模型,计算了相同侵彻条件下不同截卵平台直径时,头部侵彻产生的偏转函数和偏转角速度。结果表明:截卵平台头部侵彻产生的偏转力矩会减小弹轴与靶标法线之间的夹角,且随着截卵平台直径的增大,偏转力矩增大,偏转角速度增大,当截卵平台直径增大到1.5倍时,偏转力矩增大到约1.2倍,当截卵平台直径增大到2.0倍时,偏转力矩增大到约2倍;相同截卵平台直径下,随着头部形状系数的减小,偏转力矩和偏转角速度增大。Abstract: In order to study the problem of preliminary ballistic deflection when a penetration warhead with truncated ovate nose penetrates target obliquely, the theoretical model and simulation model are set up to analyze the effects of truncated ovate nose diameter on the preliminary ballistic deflection of penetration. The defection function and deflection angular velocity of warhead penetrating target are calculated in the same penetration condition with different truncated ovate nose diameters. It is concluded that the angle between axis of warhead and normal of target decreases with the action of deflection moment caused by the warhead of truncated ovate nose penetrating target. Both the deflection moment and the deflection angular velocity increase with the increase of truncated ovate nose diameter. When the truncated ovate nose diameter increases to 1.5 times, the deflection moment will increase to about 1.2 times. When the truncated ovate nose diameter increases to 2.0 times, the deflection moment will increase to about 2 times. Under the condition of the same truncated ovate nose diameter, the defection moment and the deflection angular velocity are increased with the shape coefficient of warhead.
-
Key words:
- penetration /
- ballistic /
- ballistic deflection /
- deflection moment
-
侵彻弹道偏转是考核战斗部侵彻性能的主要参数之一,尤其斜侵彻多层钢筋混凝土靶时,侵彻弹道偏转控制尤为重要。侵彻过程可分为弹头侵彻和弹体侵彻两个阶段:弹头率先侵彻,因弹头上下表面受力不均匀而产生偏转力矩,导致战斗部发生初始偏转;弹体侵彻时,受力不均匀产生的偏转力矩使战斗部在初始偏转的基础上进一步偏转,偏转程度可能减小也可能增强,直至战斗部完成侵彻。偏转力矩的大小与弹头形状、弹体外型直接相关,因此研究战斗部侵彻弹道的偏转性能时,分析弹头形状和弹体外型的影响至关重要。
截卵形是一种典型的侵彻战斗部头部形状。关于截卵形头部弹体的侵彻性能,我国已开展了大量的研究工作,如异形头部弹体的侵彻机理[1]、不同头部形状弹体侵彻混凝土性能[2]、斜侵彻弹道偏转[3]、截卵形弹丸侵彻混凝土性能[4]、截卵形弹体正侵彻加强筋结构靶[5]、非对称作用下的动态响应[6]、截卵形弹入水性能[7]、考虑弹头形状[8]和头部变形[9]的混凝土靶侵彻深度等。同时,研究人员还开展了相关侵彻仿真计算[10],分析了尖卵形和截卵形头部结构对侵彻深度、侵彻弹道偏转的影响规律,建立了特定条件下侵彻特定靶标时弹靶相互作用的理论模型和仿真模型,得出了相对于尖卵形头部,截卵形头部的侵彻弹道稳定性较好而侵彻深度相对降低的定性结论,为战斗部设计提供了参考,但是对于战斗部结构优化设计,尤其在重点考虑弹道偏转控制和不同侵彻条件时,如何设计截卵形头部平台尚未开展深入研究。基于此,本研究拟通过理论分析和数值仿真计算,重点探讨截卵形头部不同平台直径对初始弹道偏转的影响,为该类型战斗部头部结构优化设计提供参考。
1. 理论模型
以截卵平台中心
O 为原点,沿弹轴方向和垂直弹轴方向建立Oxy 坐标系。如图1所示,设弹轴与靶面法线的夹角为φ ;弹体圆弧段称为头部,弹体头部弧线延长线与弹轴的夹角为θ ;截卵平台直径为d。侵彻着靶时率先接触靶面的头部侧面称为下表面,通过弹轴与之对称的头部侧面称为上表面。为使问题得到简化,作如下假设:
(1) 弹体为刚体,头部侵彻靶标时不发生变形;
(2) 不同截卵平台直径下,侵彻条件相同;
(3) 相同的侵彻条件下,侵彻过程中压垮单位面积混凝土产生的阻力相同;
(4) 侵彻速度方向与弹轴的方向相同,即攻角为0°;
(5) 理论计算中,进行对比的两个截卵平台直径变化幅度较小,即计算步长较小,以减小因不均衡力的力臂按照平均力臂计算所带来的误差。
侵彻时,存在弹轴与靶面法线的夹角
φ ,弹体头部上下表面受到的阻力不同,导致沿平行弹轴方向和垂直弹轴方向产生不均衡力,该不均衡力以弹体质心为支点,产生偏转力矩;另外截卵平台的存在也会产生偏转力矩,进而引起弹体偏转。本研究只探讨截卵形头部侵彻产生的偏转,暂不分析弹体侵彻产生的偏转问题。1.1 截卵平台直径对初始偏转的影响
根据基本假设,弹体头部侵彻时,分为截卵平台侵彻和弧面段侵彻。截卵平台侵彻时产生的平行弹轴的力Fpp为
Fpp=18F0πd2icos2φ (1) 式中:
F0 为侵彻过程中压垮单位面积混凝土产生的阻力,di 为截卵平台直径(i = 1, 2, 3, …)。垂直弹轴的力Fpv为零。截卵平台直径
d1 、d2 (d2>d1 )之间的弧段侵彻时,下端面产生的平行弹轴的力Fup为Fup=12F0π[12d2cosφ+tanφ(d2−d1)2tanθ]2sinθ (2) 下端面产生的垂直弹轴的力Fuv为
Fuv=12F0π[12d2cosφ+tanφ(d2−d1)2tanθ]2cosθ (3) 上端面产生的平行弹轴的力Fdp为
Fdp=12F0π[12d2cosφ−tanφ(d2−d1)2tanθ]2sinθ (4) 上端面产生的垂直弹轴的力Fdv为
Fdv=12F0π[12d2cosφ−tanφ(d2−d1)2tanθ]2cosθ (5) 截卵平台侵彻时产生的平行弹轴的力臂按
di/4 (i=1, 2)计算,在截卵平台直径d1 、d2 之间的弧段侵彻时产生的平行弹轴方向不均衡力的平均力臂按(d1+d2)/4 计算,垂直弹轴方向的不均衡力的平均力臂按L计算,其中L为弹体质心距截卵平台中心的距离。设减小弹轴与靶标法线夹角的力矩方向为正方向,增大弹轴与靶标法线夹角的力矩方向为负方向。截卵平台直径分别为
d2 、d1 的弹体头部侵彻靶标时,平行弹轴方向产生的偏转力矩差ΔMp为ΔMp=132F0π(d32−d31)cos2φ+18F0πsinθ(d1+d2){[12d2cosφ+tanφ(d2−d1)2tanθ]2−[12d2cosφ−tanφ(d2−d1)2tanθ]2} (6) 垂直弹轴方向产生的偏转力矩差ΔMv为
ΔMv=12F0πcosθL{[12d2cosφ+tanφ(d2−d1)2tanθ]2−[12d2cosφ−tanφ(d2−d1)2tanθ]2} (7) 通过式(6)和式(7)可以得出,在
F0 、φ 、θ 、L 确定的条件下,平行弹轴方向产生的偏转力矩差和垂直弹轴方向产生的偏转力矩差为d1 、d2 的函数,即平行弹轴方向产生的偏转力矩Mp和垂直弹轴方向产生的偏转力矩Mv分别可表达为Mp=F0πf1(d1,d2) (8) Mv=F0πLf2(d1,d2) (9) 式中:
f1(d1,d2) 和f2(d1,d2) 称为偏转函数。在
θ 、φ 不同的工况(见表1)下,设d2=kd1(k>1) ,d1=40mm ,计算得出f1(d1,d2) 和f2(d1,d2) 的曲线,如图2、图3所示。表 1 计算工况Table 1. Conditions of calculationCase No. θ/(°) φ/(°) Case No. θ/(°) φ/(°) 1 25 10 5 35 10 2 25 15 6 35 15 3 25 20 7 35 20 4 25 30 8 35 30 从图2中可以得出:随着截卵平台直径的增大,侵彻时平行弹轴方向产生的偏转力矩为正,且逐渐增大;随着
φ 的增大,平行弹轴方向产生的偏转力矩增大;随着θ 的增大,即头部形状系数(CRH)减小,在相同的侵彻条件下,平行弹轴方向产生的偏转力矩为正,且随着截卵平台直径的增大而增大,但增大幅度呈现减小趋势。从图3中可以得出:随着截卵平台直径的增大,侵彻时垂直弹轴方向产生的偏转力矩为正,且逐渐增大;随着
φ 的增大,垂直弹轴方向产生的偏转力矩增大;随着θ 的增大,即头部形状系数减小,在相同的侵彻条件下,垂直弹轴方向产生的偏转力矩为正,且随着截卵平台直径的增大而增大,但增大幅度呈现减小趋势。综上所述,随着截卵平台直径的增大,侵彻时产生的平行和垂直弹轴方向的偏转力矩均为正,且均逐渐增大,即产生的偏转力矩均使弹轴与靶标法线夹角减小,且随着截卵平台直径的增大,减小弹轴与靶标法线夹角的幅度增大,偏转角速度增大;随着头部角度
θ 的增大,即头部变钝,在相同的侵彻条件下,产生的偏转力矩减小,致使减小弹轴与靶标法线夹角的幅度减小,偏转角速度减小。2. 仿真模型
2.1 模型建立
应用ANSYS/LS-DYNA程序建立仿真模型。对于壳体,采用与应变率相关的塑性随动硬化模型描述,并考虑失效;对于混凝土材料,采用JHC累计损伤本构模型描述,钢筋采用beam梁单元描述,钢筋与混凝土之间设置耦合关系;弹体与靶板之间采用面面侵蚀接触算法。本计算采用cm-g-μs单位制,每20 μs输出一次结果文件。弹体和靶标网格均采用拉格朗日网格算法,利用TRUEGRID前处理软件完成有限元模型建立。考虑到模型结构具有对称性,建立1/2计算模型,计算到弹体头部完全侵入靶标为止。
侵彻模型中,壳体材料选用G50高强度钢,壳体材料参数见表2[11],其中:
ρ 为密度,E为弹性模量,μ 为泊松比,σ 为动态屈服应力,Et 为塑性剪切模量,β 为各向同性硬化系数,fs 为失效应变系数。混凝土的JHC本构模型参数和材料参数见表3[11],其中:G为剪切模量,A为归一化黏聚强度系数,B为归一化压力硬化系数,C为应变率影响系数,N为压力硬化指数,fc为靶标抗压强度,T为抗拉强度,˙ε0 为参考应变率,σfmin 为最小断裂应变,Sfmax为归一化最大强度,pc为压碎压力,μc 为压碎体积应变,pL为压实压力,μL 为压实体积应变,D1、D2为损伤常数,K1、K2、K3为混凝土材料常数。2.2 模型校核
结合某240 kg弹体(截卵形头部)以842 m/s的速度、0°攻角、15°着角侵彻5层C40钢筋混凝土靶(厚度分别为300、180、180、180 、180 mm,每层间隔2.7 m)试验结果,对模型进行校核。该试验中,出靶后5 m处的余速为570 m/s,弹体向下偏转22.6°,试验结果如图4所示。弹体的结构参数:弹长1150 mm;弹径280 mm;头部形状系数3.0;截卵平台直径40 mm;质心距头部顶点距离640 mm;x、y、z轴的转动惯量Ix = 2.68 kg·m2,Iy = 43.90 kg·m2,Iz = 43.90 kg·m2。
应用建立的数值仿真模型参数,对侵彻试验进行了仿真计算,着重计算侵彻前后弹体的偏转和速度变化情况,并与试验结果进行对比。仿真计算结果如图5和图6所示。
仿真计算结果显示,战斗部以850 m/s的速度、15°着角侵彻5层钢筋混凝土靶时,出靶后余速为568 m/s,弹体向下偏转24°。仿真计算结果与试验结果基本一致,表明仿真模型及参数准确。
2.3 计算结果
应用校核过的模型及参数,对3个算例进行仿真计算,着重研究不同截卵平台直径的弹体头部侵彻靶标时产生的偏转情况。计算从截卵平台撞击靶标瞬间开始,到弹体头部完全侵入靶标为止。仿真计算模型如图7所示。计算所用的弹体结构参数见表4,其中:Lp为弹体长度,dp为弹体外径,mp为弹体质量,d为截卵直径,S为质心与战斗部头部顶端面的距离。不同截卵平台直径的弹体头部侵彻靶标时产生的偏转角速度和偏转角度仿真计算结果如图8和图9所示。
表 4 仿真算例参数Table 4. Parameters of simulation examplesCase Lp/mm dp/mm mp/kg θ/(°) d/mm S/mm Ix/(kg·m2) Iy/(kg·m2) Iz/(kg·m2) A 1534.4 300 400 25 40 828.5 5.63 80.52 80.52 A1 1534.4 300 400 35 40 828.1 5.64 80.60 80.60 B 1535.3 300 400 25 60 828.9 5.63 80.39 80.39 C 1532.0 300 400 25 80 828.9 5.64 79.85 79.85 通过图8和图9可以得出:在相同的侵彻条件下,随着截卵平台直径的增加,偏转角速度随之增大;在相同的截卵平台直径和相同的侵彻条件下,随着头部形状系数的减小,即头部变钝,偏转角速度随之增大。
从图2和图3中的工况2可得,当k = 1.5时(d1 = 40 mm,d2 = 60 mm),f1 = 0.05,f2 = 0.06,综合偏转力矩峰值(平行弹轴方向产生的偏转力矩和垂直弹轴方向产生的偏转力矩矢量和)约为40 mm截卵平台直径时的1.5倍;当k = 2.0时(d1 = 40 mm,d2 = 80 mm),f1 = 0.494,f2 = 0.362,综合偏转力矩峰值约为40 mm截卵平台直径时的2倍。由图8可得:截卵平台直径为60 mm时的综合偏转力矩峰值约为平台直径40 mm时的1.2倍,截卵平台直径为80 mm时的综合偏转力矩峰值约为平台直径40 mm时的2倍。
从图2和图3中的工况2和工况6可得,对于d1 = 40 mm:当k = 1.1时,工况2与工况6的f1比值约为1.0,f2的比值约为1.1,工况2的综合偏转力矩约为工况6的综合偏转力矩的1.5倍;当k = 1.5时,工况2与工况6的f1比值约为1.0,f2的比值约为1.2,工况2的综合偏转力矩约为工况6的综合偏转力矩的1.5倍。由图9可得,
θ 分别为25°和35°时,综合偏转力矩峰值的比值约为1.4。对比理论计算和数值仿真结果得出,两种计算结果得到的趋势相同,理论计算结果偏高。
3. 结 束 语
通过理论分析和数值仿真计算,分析了在相同的侵彻条件下不同截卵平台直径对初始侵彻弹道偏转的影响。仿真结果与理论分析结果得到的趋势一致,数值模拟结果基本相同。结果表明:截卵形头部侵彻产生的偏转力矩均为减小弹轴与靶标法线夹角,且截卵平台直径越大,弹头侵彻时产生的偏转力矩越大,偏转角速度增大;相同截卵平台直径下,随着头部形状系数的减小,偏转力矩和偏转角速度随之增大。
由于理论研究时作出部分假设,与实际情况存在一定偏差,如弹头侵彻是动态变化过程,而非单一状态过程,因此理论计算结果偏高。下一步工作重点是考虑弹头动态变化的理论模型精确描述。
-
表 1 计算工况
Table 1. Conditions of calculation
Case No. θ/(°) φ/(°) Case No. θ/(°) φ/(°) 1 25 10 5 35 10 2 25 15 6 35 15 3 25 20 7 35 20 4 25 30 8 35 30 ρ/(g·cm−3) E/GPa μ σ/MPa Et/MPa β fs 7.85 210 0.2 1900 0 1 0.8 ρ/(g·cm−3) G/GPa A B C N fc/MPa 2.44 11.147 0.79 1.6 0.007 0.61 40 T/MPa ˙ε0/s−1 σfmin Sfmax pc/MPa μc pL/MPa 3.28 1 0.01 7 9.34 0.008 800 μL D1 D2 K1/GPa K2/GPa K3/GPa fs 0.116 0.04 1 85 −171 208 0.1 表 4 仿真算例参数
Table 4. Parameters of simulation examples
Case Lp/mm dp/mm mp/kg θ/(°) d/mm S/mm Ix/(kg·m2) Iy/(kg·m2) Iz/(kg·m2) A 1534.4 300 400 25 40 828.5 5.63 80.52 80.52 A1 1534.4 300 400 35 40 828.1 5.64 80.60 80.60 B 1535.3 300 400 25 60 828.9 5.63 80.39 80.39 C 1532.0 300 400 25 80 828.9 5.64 79.85 79.85 -
[1] 刘坚成, 黄风雷, 皮爱国, 等. 异型头部弹体增强侵彻性能机理研究 [J]. 爆炸与冲击, 2014, 34(4): 409–414. doi: 10.11883/1001-1455(2014)04-0409-06LIU J C, HUANG F L, PI A G, et al. On enhanced penetration performance of modified nose projectiles [J]. Explosion and Shock Waves, 2014, 34(4): 409–414. doi: 10.11883/1001-1455(2014)04-0409-06 [2] 彭永, 方秦, 吴昊, 等. 不同头部形状弹体侵彻混凝土靶体的终点弹道参数分析 [J]. 兵工学报, 2014, 35(Suppl 2): 128–134.PENG Y, FANG Q, WU H, et al. Theoretical analyses for terminal ballistic of the projectiles with different nose geometries penetrating into concrete targets [J]. Acta Armamentarii, 2014, 35(Suppl 2): 128–134. [3] 葛超, 董永香, 陆志超, 等. 弹丸头部对斜侵彻弹道偏转影响研究 [J]. 兵工学报, 2015, 36(2): 255–262. doi: 10.3969/j.issn.1000-1093.2015.02.010GE C, DONG Y X, LU Z C, et al. Ballistic deflection on oblique penetration of projectiles with different noses [J]. Acta Armamentarii, 2015, 36(2): 255–262. doi: 10.3969/j.issn.1000-1093.2015.02.010 [4] 贺登高, 贺虎成. 基于Teland模型的截卵形弹丸侵彻混凝土计算分析 [J]. 兵器材料科学与工程, 2010, 33(5): 78–80. doi: 10.3969/j.issn.1004-244X.2010.05.021HE D G, HE H C. Calculation and analysis of rigid truncated-ogive-nosed projectile penetrating into concrete based on Teland force model [J]. Ordnance Material Science and Engineering, 2010, 33(5): 78–80. doi: 10.3969/j.issn.1004-244X.2010.05.021 [5] 展婷变, 吕淑芳, 黄德雨. 截卵形弹体正侵彻加强筋结构靶的理论分析 [J]. 弹道学报, 2012, 24(1): 52–57. doi: 10.3969/j.issn.1004-499X.2012.01.011ZHAN T B, LÜ S F, HUANG D Y. Theoretical analysis on normal penetration of truncated oval-nosed projectile into stiffened plate [J]. Journal of Ballistics, 2012, 24(1): 52–57. doi: 10.3969/j.issn.1004-499X.2012.01.011 [6] 杨华伟, 王志华, 晋小超. 高速侵彻混凝土弹体在横向非对称作用下的动态响应 [J]. 高压物理学报, 2017, 31(1): 73–80. doi: 10.11858/gywlxb.2017.01.011YANG H W, WANG Z H, JIN X C. Dynamic response of high-speed projectiles penetrating into concrete target by asymmetric load [J]. Chinese Journal of High Pressure Physics, 2017, 31(1): 73–80. doi: 10.11858/gywlxb.2017.01.011 [7] 郭子涛, 张伟, 郭钊, 等. 截卵形弹水平入水的速度衰减及空泡扩展特性 [J]. 爆炸与冲击, 2017, 37(4): 727–733. doi: 10.11883/1001-1455(2017)04-0727-07GUO Z T, ZHANG W, GUO Z, et al. Characteristics of velocity attenuation and cavity expansion induced by horizontal water-entry of truncated-ogive nosed projectiles [J]. Explosion and Shock Waves, 2017, 37(4): 727–733. doi: 10.11883/1001-1455(2017)04-0727-07 [8] 吴昊, 方秦, 龚自明. 考虑刚性弹弹头形状的混凝土(岩石)靶体侵彻深度半理论分析 [J]. 爆炸与冲击, 2012, 32(6): 573–580. doi: 10.11883/1001-1455(2012)06-0573-08WU H, FANG Q, GONG Z M. Semi-theoretical analyses for penetration depth of rigid projectiles with different nose geometries into concrete (rock) targets [J]. Explosion and Shock Waves, 2012, 32(6): 573–580. doi: 10.11883/1001-1455(2012)06-0573-08 [9] 赵军, 陈小伟, 金丰年, 等. 考虑头形磨损变化的动能弹侵彻深度研究 [J]. 力学学报, 2010, 42(2): 212–218. doi: 10.6052/0459-1879-2010-2-2009-009ZHAO J, CHEN X W, JIN F N, et al. Studying on the penetration depth of penetrator with including the effect of mass abrasion [J]. Chinese Journal of Theoretical and Applied Mechanics, 2010, 42(2): 212–218. doi: 10.6052/0459-1879-2010-2-2009-009 [10] 邓佳杰, 张先锋, 葛贤坤, 等. 基于局部相互作用理论的侵彻弹头部形状优化及仿真 [J]. 爆炸与冲击, 2017, 37(4): 611–620. doi: 10.11883/1001-1455(2017)04-0611-10DENG J J, ZHANG X F, GE X K, et al. Nose-shape optimization and simulation of projectiles penetrating into concrete target based on local interaction theory [J]. Explosion and Shock Waves, 2017, 37(4): 611–620. doi: 10.11883/1001-1455(2017)04-0611-10 [11] 张丁山, 谷鸿平, 吕永柱, 等. 战斗部后端盖结构强度的数值仿真及应力波分析方法 [J]. 探测与控制学报, 2018, 40(5): 21–25.ZHANG D S, GU H P, LÜ Y Z, et al. Numerical simulation of structure strength calculation and stress waves analysis [J]. Journal of Detection and Control, 2018, 40(5): 21–25. -