破甲弹战斗部侵彻高速运动靶板的数值模拟

2016-09-13山西中北大学机电工程学院孙国庆张国伟乔磊

河北农机 2016年2期

山西中北大学机电工程学院　孙国庆　张国伟　乔磊

破甲弹战斗部侵彻高速运动靶板的数值模拟

山西中北大学机电工程学院孙国庆张国伟乔磊

由于国内外有关破甲弹对运动体侵彻课题研究很少，可借鉴的资源有限。本文主要使用理论分析和数值模拟相结合的方法，利用有限元软件ANSYS/LS-DYNA对破甲弹战斗部侵彻运动靶板进行仿真，得出了入射角和运动体速度对侵彻射流影响的一般性规律，为今后实验研究打下基础。

破甲弹；运动靶板；有限元；数值模拟

引言

靶板运动速度对射流破甲威力的影响主要表现在靶板对射流的横向剪切力矩。相对于静止靶板，运动靶板的剪切力矩会使射流发生严重的弯曲、偏转、变形，导致射流的断裂和分散。随着靶板速度的增加，靶板剪切力矩的影响也逐渐增大，干扰影响也趋于明显，射流变形越剧烈，对其后续侵彻能力也将产生较大的影响［1］。运动靶板对射流的干扰影响还体现在对射流能量的无效消耗上。静止靶板的侵彻，射流只在轴向方向上侵彻消耗能量，形成笔直的侵彻孔；而运动靶板的侵彻，射流还要在靶板横向作用，浪费射流能量，扩大侵彻孔径，形成上粗下窄的倾斜的侵彻孔。随着靶板运动速度的增加，倾斜孔径越大，消耗射流能力越大［2］。

1　计算模型

将靶板速度取为 400m/s、600m/s、800m/s三种情况分别进行动态斜侵彻分析仿真计算。战斗部几何模型如图1所示。标靶板五种材料的破甲战斗部有限元模型，调整适宜的网格密度，设置恰当的边界约束（空气域无反射边界约束），使用正确的算法（ALE算法），确保计算结果准确。求解计算时的单位为cm-g-μs-Mbar。

图1　战斗部示意图建立药型罩、主装药、壳体、空气和目

考虑到靶板的运动情况特意将靶板偏离药型罩，以保证在整个侵彻过程中，射流尽可能地作用于靶板的中部，避免受到边界效应的影响，这也是建模时确定靶板尺寸的依据［3］。考虑到战斗部和目标靶板结构的轴对称性，简化模型为二分之一模型，主装药、药型罩和空气设置多物质ALE算法，采用Euler网格建模，靶板和壳体采用Lagrange网格建模，并且各材料之间使用耦合算法。

2　数值模拟

图2、图3、图4分别是靶板速度为v=400m/s，v=600m/s，v=800m/s时不同时刻射流侵彻过程图。

图2　v=400m/s时不同时刻射流侵彻过程

图3　v=600m/s时不同时刻射流侵彻过程

图4　v=800m/s时不同时刻射流侵彻过程

3　结果分析

运动靶板对射流侵彻体的影响，主要表现在其对射流侵彻体的横向扰动作用，使射流侵彻体弯曲和失稳，从而降低破甲能力。从图2中，我们可以很明显地看出射流侵彻过程中，由于受运动靶板的横向扰动干扰射流变形弯曲，并在侵彻孔内左右摆动。为了更细致地分析各射流段受横向干扰的影响，特将射流的侵彻过程分为三个阶段。在射流侵彻的初始阶段，即射流与靶板侧壁接触前，射流头部几乎没有弯曲变形，说明射流头部具有较高的速度，抗干扰能力相对较强，如图2中96μs以前；在射流开始出现弯曲变形阶段，即射流受到靶板的横向扰动开始弯曲，射流头部高速段和较高速段同时受到干扰变形弯曲，其产生的影响大小也因为射流段不同而不同，射流头部影响变形较小，射流中段影响变形较大，这说明速度越低对干扰越敏感，如图2中108μs左右；在射流侵彻震荡阶段，即射流在侵彻孔内碰撞，射流不断受到靶板的横向干扰弯曲变形更剧烈、更复杂，并在侵彻孔左右侧摆动，如图2中108μs以后。

运动靶板速度v=600m/s和v=800m/s时，对于同一种射流的侵彻过程，射流弯曲变形情况大致相同，如图3和图4所示。第一阶段，射流几乎没有变形；第二阶段，射流弯曲变形都发生在射流的中部低速段；第三阶段，射流震荡向下侵彻。从以上分析可以得出，运动靶板对射流的影响主要体现在对射流中部低速段的横向干扰，速度越低干扰越大，经过干扰后的射流作用于射流高速段头部，进而影响射流侵彻。

4　结论

与侵彻静止和相对静止靶板相比，射流侵彻运动靶板时的主要区别在于射流具有了一定的径向速度vjx。射流径向速度vjx的干扰作用主要体现在：其大小的变化和方向的无序性导致射流不同部位受到的干扰方向不同，加大了对射流的影响，最终导致射流容易弯曲甚至断裂。

［1］周权.聚能装药破甲侵彻深度与穿孔直径多目标优化设计研究［D］.南京：南京理工大学，2007.

［2］隋树元，王树山.终点效应学［M］.北京：国防工业出版社，2000：196-232.

［3］耿世才.小口径破甲弹的研究［D］.南京理工大学硕士论文.2010：5-15.自适应问题，需要对管道机器人通过弯管的情况进行仿真分析。假定弯管机器人重心位置为几何中心，则对弯管设定路径配合，其重心会沿管道重心轴线运动，对管壁施加压力的方向上的偏航控制为10mm。最后，检查运动干涉，确定运动中无过定义和欠定义的情况发生。