CN114297860A - 延时引信弹药与跳弹式防护结构碰撞分析方法 - Google Patents

延时引信弹药与跳弹式防护结构碰撞分析方法 Download PDF

Info

Publication number
CN114297860A
CN114297860A CN202111649114.5A CN202111649114A CN114297860A CN 114297860 A CN114297860 A CN 114297860A CN 202111649114 A CN202111649114 A CN 202111649114A CN 114297860 A CN114297860 A CN 114297860A
Authority
CN
China
Prior art keywords
collision
projectile
impact
speed
bullet
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111649114.5A
Other languages
English (en)
Other versions
CN114297860B (zh
Inventor
秦有权
吴应祥
陶西贵
马媛媛
张伟锋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Academy of Defense Engineering of PLA Academy of Military Science
Original Assignee
National Academy of Defense Engineering of PLA Academy of Military Science
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by National Academy of Defense Engineering of PLA Academy of Military Science filed Critical National Academy of Defense Engineering of PLA Academy of Military Science
Priority to CN202111649114.5A priority Critical patent/CN114297860B/zh
Publication of CN114297860A publication Critical patent/CN114297860A/zh
Application granted granted Critical
Publication of CN114297860B publication Critical patent/CN114297860B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明提出一种延时引信弹药与跳弹式防护结构碰撞分析方法，预设初始飞行状态下弹体的攻角为零且弹体沿铅垂方向入射，弹体轴线、弹体偏转力、弹体速度均在入射平面内，将三维问题转化为二维问题，通过弹体与异形结构靶板碰撞后的运动轨迹与碰撞后的运动参数的分析，计算二次防弹障碍物碰撞点的位置，可为跳弹式防钻地弹结构的研发设计以及对周边目标的毁伤分析提供理论依据。

Description

延时引信弹药与跳弹式防护结构碰撞分析方法
技术领域
本发明涉及侵彻防护技术领域，特别是涉及一种延时引信弹药与跳弹式防护结构碰撞分析方法。
背景技术
延期触发引信,简称延期引信，其装有延期元件或延期装置,与目标碰撞后能延期一段时间起作用引爆。当前很多钻地战斗部都采用延期引信，该钻地战斗部依靠自身坚硬的弹头钻过土壤、岩石、混凝土等介质，进入地下工程内部或穿透多层楼板进入建筑物下层精确延时爆炸，这类高效能武器对战场目标的防护以及其他重要目标的生存安全构成了极大的威胁。因此，现有技术中进行了很多针对这类武器的研究。合理的安全防护结构是这些研究试验中需要关注的问题之一，例如，为减少钻地弹的侵彻深度，地下工程中的遮弹层设计或表面异形结构这类遮弹层的研究和检测也需要做弹靶斜碰试验。同时，合理的安全防护结构也是部队实弹训练所面临的现实问题，实弹打靶训练中经常会因为弹靶斜碰造成跳弹事故。现有防护措施是在靶体周围设置强度一致的墙板式防弹仓，这类没有重点、四面平均的设防形式十分浪费材料，并且由于缺少科学的弹道预测技术无法完全避免跳弹引起的安全事故。因此，如何有重点地设置防弹仓的角度和强度是训练和试验中均需要解决的问题。
异形靶板在防护工程中的应用十分广泛。现阶段防护工程领域中异形靶板多为固定靶板。现有的异形靶板有半球形、两种不同直径半球组合形式、半球和其他异形结构组合形式。其中两种不同直径半球组合形式为大球和小球组合结构，即在两个大的球壳之间的防护薄弱处采用小球壳进行辅助以增强整体防护效果。现有技术中半球形异形靶板的结构设计是基于半穿甲弹侵彻遮弹层模型试验得到的，并在分析实验结果的基础上确定了各个结构的材料；然后根据侵彻经验公式进行计算，结果表明采用组合形式结构比单一结构在抗侵彻方面更有优势。现有技术中还提出了对异形结构，例如，半球、圆锥以及正三棱锥，导偏模型的分析方法，但这种分析方法是基于MAT-JOHNSOM-COOK材料模型以及GRUNEISEN状态方程通过数值模拟进行的仿真实验，之后，通过现场试验对数值仿真模型进行验证，给出该类异形靶板的结构形状参数，以便结构参数对导偏性能的影响；这种仿真模型的对象为全金属材料，且其防护对象为小口径弹体，并不能直接确定这种针对小口径弹体的设计方法在大口径弹体的防护结构中的适用性。
发明内容
针对现有技术中存在的问题，本发明的目的在于提出一种延时引信弹药与跳弹式防护结构碰撞分析方法，提出了弹体的第一次碰撞模型和第二次碰撞模型，根据首次碰撞后的修正弹体角速度ω'2以及首次撞击后弹体质心速度为V'C所在的方向与首次撞击的碰撞面MN的夹角β，建立第二次碰撞模型；得出弹体与首次碰撞的异形靶板相邻的异形靶板的碰撞面上的碰撞点B的位置。本发明解决了炮击或侵彻试验过程中缺乏针对弹体跳弹或发生大角度偏转情况下的试验安全防护措施的问题。通过本发明所建立的碰撞模型和碰撞分析方法，可以对碰撞后弹体的飞行轨迹进行预测，指导试验场中安全防护设施，例如，收弹仓的设计。现有技术中收弹仓的壁厚是直接借鉴人防工程设计规范中防护层的厚度，本发明的碰撞模型可以从理论上对收弹仓的壁厚进行计算避免试验过程中的危险。
本发明中所提出一种延时引信弹药与跳弹式防护结构碰撞分析方法，预设初始飞行状态下弹体的攻角为零且弹体沿铅垂方向入射，弹体轴线、弹体偏转力、弹体速度均在入射平面内，将三维问题转化为二维问题，根据首次碰撞后弹体的运动轨迹，确定二次碰撞点的位置，可为跳弹式防钻地弹结构的研发设计以及对周边目标的毁伤分析提供理论依据。
本发明的技术方案如下：
一种延时引信弹药与跳弹式防护结构碰撞分析方法，其具体包括如下步骤：
步骤一：预设弹体为线性杆状，初始飞行状态下所述弹体的攻角为零且防侵彻结构的碰撞面光滑；
步骤二：建立第一次碰撞模型；
步骤三：求出碰撞后弹头的速度v'A以及首次碰撞后弹头速度在Y轴方向的分量投影v'Ay；
步骤四：确定材料恢复系数e；得到首次碰撞后弹头速度在Y轴方向的分量投影v'Ay与材料恢复系数e、首次撞击前的弹体的质心速度vc以及弹体的入射方向与首次撞击面之间的夹角θ之间的关系；
步骤五：计算首次碰撞后弹体质心速度在Y方向上的分量v'cy；
步骤六：最终得到首次碰撞后弹体的角速度ω2以及首次撞击后弹体质心速度v'c所在的方向与首次撞击的碰撞面MN的夹角β；
Figure BDA0003446056100000031
Figure BDA0003446056100000032
其中，e为材料恢复系数；
θ为弹体的入射方向与首次撞击面之间的夹角；
L为弹体长度；
vc为首次撞击前的弹体的质心速度；
v'cy为首次碰撞后弹体质心速度在Y方向上的分量；
v'cx为首次撞击后弹体质心速度在X轴方向的分量；
步骤七：对弹体的质心位置进行修正；修正后，其弹体质心与弹头顶部的距离L'；
步骤八：得到修正后的首次碰撞后弹体角速度ω'2；
步骤九：根据首次碰撞后的修正弹体角速度ω'2以及首次撞击后弹体质心速度为v'c所在的方向与首次撞击的碰撞面MN的夹角β，建立第二次碰撞模型；
步骤十：根据异形靶板的几何尺寸，确定弹体与首次碰撞的异形靶板相邻的异形靶板的碰撞面上的碰撞点B的位置。
优选地，步骤二中，异形靶板的第一尖顶部的两个侧面为异形靶板的撞击面，通过所述撞击面进行导偏，根据撞击前弹体的质心速度、撞击面上的弹着点，弹体入射角度建立第一次碰撞模型。
优选地，步骤三中，所述首次碰撞后弹头速度在Y轴方向的分量投影v'Ay带动所述弹体进行偏转；
所述碰撞后弹头的速度v'A如下所示：
v'A＝v'c+v'AC (5)
式中，vA为首次撞击前弹头的速度；
v'A为首次撞击后弹头的速度；
v'c为首次撞击后弹体质心的速度；
v'AC为首次碰撞后弹头相对于弹体质心的速度；
首次碰撞后弹头速度在Y轴方向的分量投影v'Ay如下式所示：
Figure BDA0003446056100000041
其中，v'cy为首次碰撞后弹体质心速度在Y方向上的分量；
θ为弹体的入射方向与首次撞击面之间的夹角；
L为弹体长度。
优选地，步骤四中，首次碰撞后弹头速度在Y轴方向的分量投影v'Ay与材料恢复系数e、首次撞击前的弹体的质心速度vc以及弹体的入射方向与首次撞击面之间的夹角θ之间的关系如下；
Figure BDA0003446056100000051
式中，y'Ay为首次碰撞后弹头速度在Y轴方向的分量投影；
vAy为初始状态下弹头速度沿着Y轴方向的分量投影；
vc为首次撞击前的弹体的质心速度；
θ为弹体的入射方向与首次撞击面之间的夹角。
优选地，首次碰撞后弹体角速度ω'2的修正后的表达式为：
Figure BDA0003446056100000052
式中，θ为弹体的入射方向与首次撞击面之间的夹角；
L'为弹体质心与弹头顶部的距离；
m为弹体的质量；
J'c为对弹体进行质心位置修正后，弹体的转动惯量。
优选地，步骤五中，首次碰撞后弹体质心速度在Y方向上的分量v'cy与首次碰撞后弹体的角速度ω2之间的关系式；
Figure BDA0003446056100000053
其中，L为弹体长度；
θ为弹体的入射方向与首次撞击面之间的夹角；
vc为首次撞击前的弹体的质心速度。
优选地，首次碰撞模型的建立方法如下，其具体包括如下步骤：
S1：建立初始坐标系，预设表面MN方向为X轴，垂直于MN的方向Y轴，在初始坐标系中，所述弹着点为A，所述弹体的质心为C，撞击面为MN所在平面，弹体的入射方向与首次撞击面之间的夹角为θ；
S2：设置弹体的参数；
预设弹体的质量为m，弹体的长度为L；首次撞击前的弹体的质心速度为vC；首次撞击后弹体的质心速度为v'c；
S3：建立弹体首次碰撞模型，其表达式如下所示：
mv'cx-mvcx＝∑Ix (1)
mv'cy-mvcy＝∑Iy (2)
JCω2-JCω1＝∑MC(I(e)) (3)
其中，vcx为撞击前弹体质心速度在X轴方向的分量；
vcy为撞击前弹体质心速度在Y轴方向的分量；
v'cx为首次撞击后弹体质心速度在X轴方向的分量；
v'cy为首次撞击后弹体质心速度在Y轴方向的分量；
Ix为首次碰撞时弹体受到X轴方向的碰撞冲量；
Iy为首次碰撞时弹体受到Y轴方向的碰撞冲量；
JC为弹体绕其质心的转动惯量；
ω1为首次碰撞前弹体的角速度；
ω2为首次碰撞后弹体的角速度；
∑MC(I(e))为首次碰撞冲量对弹体质心之矩的矢量和；
S4：弹体在碰撞之前角速度ω1＝0，弹体首次碰撞时弹体受到Y轴方向的碰撞冲量Iy，Ix＝0，首次碰撞后弹体质心速度在X轴方向的分量与首次撞击前弹体质心速度在X轴方向的分量相等，对首次碰撞模型进行修正；
首次碰撞后弹体质心速度在X轴方向的分量的表达式如下所示：
v'cx＝vcx＝vccosθ (4)
其中，vcx为撞击前弹体质心速度在X轴方向的分量；v'cx为首次撞击后弹体质心速度在X轴方向的分量；vc为首次撞击前的弹体的质心速度；θ为弹体的入射方向与首次撞击面之间的夹角。
优选地，第二次碰撞模型的建立步骤如下：
第一步，预设直线AD与首次碰撞后弹体质心速度v'c的方向平行；以直线AD为X轴，垂直于直线AD的方向为Y轴，重新建立第二次碰撞模型的定坐标系；同时，以t时刻后的弹体质心C1为原点建立动坐标系；其中t为首次碰撞至第二次碰撞之间的时间间隔；
第二步，得到弹头顶点的相对运动轨迹(x',y')；
Figure BDA0003446056100000071
其中，L'为弹体质心与弹头顶部的距离；
θ为弹体的入射方向与首次撞击面之间的夹角；
ω'2为首次碰撞后弹体角速度；
t为首次碰撞至第二次碰撞之间的时间间隔；
第三步，建立弹头顶点的牵连运动轨迹(xo',yo')，其表达式如下所示：
Figure BDA0003446056100000072
其中，φ为X轴到X'的转角；
t为首次碰撞至第二次碰撞之间的时间间隔；
L'为弹体质心与弹头顶部的距离；
vc为首次撞击前的弹体的质心速度；
第四步，经过坐标变换，得到弹头顶点的绝对运动轨迹(x，y)，其表达式如下所示：
Figure BDA0003446056100000073
其中，t为首次碰撞至第二次碰撞之间的时间间隔；
L'为弹体质心与弹头顶部的距离；
vc为首次撞击前的弹体的质心速度；
ω'2为首次碰撞后弹体角速度；
θ为弹体的入射方向与首次撞击面之间的夹角；
优选地，确定二次撞击点B的位置的方法,具体包括如下步骤：
A)确定第一次撞击点的位置A；
B)添加异形靶板的边线；
首次碰撞的异形靶板的边线MN和第二次碰撞的异形靶板的边线PN；
C)得到第二次碰撞的撞击面方程(xPN，yPN)；
D)当发生第二次碰撞时，(xPN，yPN)与(x，y)重合，得到第二次撞击点B的位置。
优选地，二次碰撞防弹障碍物或二次碰撞防弹结构物为直线结构PG，其与水平X轴夹角为
Figure BDA0003446056100000081
与X轴交点为K，弹体入射点A与K之间的距离AK＝s1。
与现有技术相比，本发明的优点在于：
本发明的延时引信弹药与跳弹式防护结构碰撞分析方法，预设初始飞行状态下弹体的攻角为零且弹体沿铅垂方向入射，弹体轴线、弹体偏转力、弹体速度均在入射平面内，将三维问题转化为二维问题，通过弹体与异形结构靶板碰撞后的运动轨迹与碰撞后的运动参数的分析，计算二次防弹障碍物碰撞点的位置，可为跳弹式防钻地弹结构的研发设计以及对周边目标的毁伤分析提供理论依据。
附图说明
本发明上述和/或附加方面的优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是根据本发明的延时引信弹药与跳弹式防护结构碰撞分析方法中异形靶板的模型示意图。
图2是根据本发明的延时引信弹药与跳弹式防护结构碰撞分析方法中第一次碰撞模型的流程图。
图3是根据本发明的延时引信弹药与跳弹式防护结构碰撞分析方法中第二次碰撞模型的流程图。
图4是根据本发明的延时引信弹药与跳弹式防护结构碰撞分析方法中优选确定二次撞击点B的位置的方法流程图。
图5是根据本发明的延时引信弹药与跳弹式防护结构碰撞分析方法的第一次碰撞模型的分析模型图，也就是说，弹体在首次与异形靶板碰撞后，如图5所示，在外力作用下会发生沿v'c方向的质心运动以及弹体本身绕着质心C的转动。
图6是根据本发明的延时引信弹药与跳弹式防护结构碰撞分析方法中v'C的方向的平行线AD为X轴，垂直于该方向为Y轴，建立定坐标系；以t时刻后的弹体质心C1为原点建立动坐标系，得到的分析模型示意图。
图7是如图6所示的根据本发明的延时引信弹药与跳弹式防护结构碰撞分析方法的分析模型图中添加异形靶板边线的分析模型示意图。
图8是根据本发明的延时引信弹药与跳弹式防护结构碰撞分析方法中二次防弹障碍物设置示意图。
图9是根据本发明的另一实施例的延时引信弹药与跳弹式防护结构碰撞分析方法的二次防弹障碍物设置位置理论结果示意图。
图10至图12为现场实验图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。
图1至图12中示出了基于矢量防护的跳弹式固定靶板碰撞模型及碰撞分析方法中弹体首次撞击固定异形靶板的情形以及弹体二次撞击固定异形靶板的情形，可根据本申请中的碰撞模型和碰撞方法进行矢量防护结构的防护模块的设计。本申请中的异形靶板包括第一顶部；第一顶部为第一尖顶部，所述第一顶部的两个侧面为碰撞面；
其具体包括如下步骤：
步骤一：预设所述弹体为线性杆状，初始飞行状态下所述弹体的攻角为零且所述防侵彻结构的碰撞面光滑；
优选地，预设所述弹体为刚体，预设弹体与异形靶板撞击时的变形为零；
优选地，预设所述弹体为线性杆状，在所述弹体初始飞行状态下，其速度与弹轴重合，即初始飞行状态下所述弹体的攻角为零，所述弹体仅受偏转力作用；
优选地，当弹体速度＞200m/s的范围内，弹体与防护结构上的弹着点之间的摩擦力忽略不计；也就是说，预设所述防护结构的碰撞面光滑；
优选地，所述弹体的重力、所述弹体在初始飞行状态下的空气阻力和空气阻力矩相对于所述弹体的冲量可忽略不计；优选地，弹体轴线、偏转力、速度均在入射平面内，可忽略弹体绕轴旋转的阻力、阻力矩，及其对攻角和弹道的影响。
步骤二：所述第一尖顶部两侧连接的两个侧面为异形靶板的撞击面，通过所述撞击面进行导偏，根据撞击前弹体的质心速度、撞击面上的弹着点，弹体入射角度建立第一次碰撞模型。
具体地，第一次碰撞模型，即首次碰撞模型的建立方法如下，其具体包括如下步骤：
S1：建立初始坐标系，预设表面MN方向为X轴，垂直于MN的方向Y轴，在初始坐标系中，所述弹着点为A，所述弹体的质心为C，撞击面为MN所在平面，弹体的入射方向与首次撞击面之间的夹角为θ；
S2：设置弹体的参数；
优选地，预设弹体的质量为m，弹体的长度为L；首次撞击前的弹体的质心速度为vc；首次撞击后弹体的质心速度为v'c；
S3：建立弹体首次碰撞模型，其表达式如下所示：
mv'cx-mvcx＝∑Ix (1)
mv'cy-mvcy＝∑Iy (2)
JCω2-JCω1＝∑MC(I(e)) (3)
其中，vcx为撞击前弹体质心速度在X轴方向的分量；
vcy为撞击前弹体质心速度在Y轴方向的分量；
v'cx为首次撞击后弹体质心速度在X轴方向的分量；
v'cy为首次撞击后弹体质心速度在Y轴方向的分量；
Ix为首次碰撞时弹体受到X轴方向的碰撞冲量；
Iy为首次碰撞时弹体受到Y轴方向的碰撞冲量；
JC为弹体绕其质心的转动惯量；
ω1为首次碰撞前弹体的角速度；
ω2为首次碰撞后弹体的角速度；
∑MC(I(e))为首次碰撞冲量对弹体质心之矩的矢量和；
S4：弹体在碰撞之前角速度ω1＝0，且由于碰撞面光滑，弹体首次碰撞时弹体受到Y轴方向的碰撞冲量Iy，Ix＝0，首次碰撞后弹体质心速度在X轴方向的分量与首次撞击前弹体质心速度在X轴方向的分量相等，对首次碰撞模型进行修正；
优选地，首次碰撞后弹体质心速度在X轴方向的分量的表达式如下所示：
v'cx＝vcx＝vccosθ (4)
其中，vcx为撞击前弹体质心速度在X轴方向的分量；v'cx为首次撞击后弹体质心速度在X轴方向的分量；vc为首次撞击前的弹体的质心速度；θ为弹体的入射方向与首次撞击面之间的夹角。
步骤三：根据碰撞后弹头的速度v'A以及首次碰撞后弹头速度在Y轴方向的分量投影v'Ay，得到首次碰撞后弹头速度在Y轴方向的分量投影v'Ay与弹体长度L、弹体的入射方向与首次撞击面之间的夹角θ，首次碰撞后弹体的角速度ω2以及首次碰撞后弹体质心速度在Y方向上的分量v'cy之间的关系；其中，所述首次碰撞后弹头速度在Y轴方向的分量投影v'Ay带动所述弹体进行偏转；
优选地，所述碰撞后弹头的速度v'A如下所示：
v'A＝v'c+v'AC (5)
式中，vA为首次撞击前弹头的速度；
v'A为首次撞击后弹头的速度；
v'c为首次撞击后弹体质心的速度；
v'AC为首次碰撞后弹头相对于弹体质心的速度；
优选地，首次碰撞后弹头速度在Y轴方向的分量投影v'Ay如下式所示：
Figure BDA0003446056100000121
其中，v'cy为首次碰撞后弹体质心速度在Y方向上的分量；
θ为弹体的入射方向与首次撞击面之间的夹角；
L为弹体长度；
步骤四：确定材料恢复系数e；得到首次碰撞后弹头速度在Y轴方向的分量投影v'Ay与材料恢复系数e、首次撞击前的弹体的质心速度vc以及弹体的入射方向与首次撞击面之间的夹角θ之间的关系；
Figure BDA0003446056100000122
式中，v'Ay为首次碰撞后弹头速度在Y轴方向的分量投影；
vAy为初始状态下弹头速度沿着Y轴方向的分量投影；
vc为首次撞击前的弹体的质心速度；
θ为弹体的入射方向与首次撞击面之间的夹角；
步骤五：得到首次碰撞后弹体质心速度在Y方向上的分量v'cy与首次碰撞后弹体的角速度ω2之间的关系式；
Figure BDA0003446056100000131
其中，L为弹体长度；
θ为弹体的入射方向与首次撞击面之间的夹角；
vc为首次撞击前的弹体的质心速度；
步骤六：最终得到首次碰撞后弹体的角速度ω2以及首次撞击后弹体质心速度为v'c所在的方向与首次撞击的碰撞面MN的夹角β；
Figure BDA0003446056100000132
Figure BDA0003446056100000133
其中，e为材料恢复系数；
θ为弹体的入射方向与首次撞击面之间的夹角；
L为弹体长度；
vc为首次撞击前的弹体的质心速度；
v'cy为首次碰撞后弹体质心速度在Y方向上的分量；
v'cx为首次撞击后弹体质心速度在X轴方向的分量；
步骤七：对弹体的质心位置进行修正；
优选地，初始状态下，弹体简化为细长杆，其绕质心的转动惯量按照
Figure BDA0003446056100000134
修正后，其弹体质心与弹头顶部的距离L'。
步骤八：得到首次碰撞后弹体角速度ω'2的修正后的表达式为：
Figure BDA0003446056100000135
式中，θ为弹体的入射方向与首次撞击面之间的夹角；
L'为弹体质心与弹头顶部的距离；
m为弹体的质量；
J'c为对弹体进行质心位置修正后，弹体的转动惯量。
步骤九：根据首次碰撞后的修正弹体角速度ω'2以及首次撞击后弹体质心速度为v'c所在的方向与首次撞击的碰撞面MN的夹角β，建立第二次碰撞模型。
优选地，弹体和异形靶板在首次碰撞过程中均未发生破坏，也就是说，在首次碰撞过程中仅发生弹塑性变形。
优选地，第二次碰撞模型的建立步骤如下所示：
第一步，预设直线AD与首次碰撞后弹体质心速度v'c的方向平行；以直线AD为X轴，垂直于直线AD的方向为Y轴，重新建立第二次碰撞模型的定坐标系；同时，以t时刻后的弹体质心C1为原点建立动坐标系；其中t为首次碰撞至第二次碰撞之间的时间间隔；
第二步，得到弹头顶点的相对运动轨迹(x',y')；
Figure BDA0003446056100000141
其中，L'为弹体质心与弹头顶部的距离；
θ为弹体的入射方向与首次撞击面之间的夹角；
ω'2为首次碰撞后弹体角速度；
t为首次碰撞至第二次碰撞之间的时间间隔；
第三步，建立弹头顶点的牵连运动轨迹(xo',yo')，其表达式如下所示：
Figure BDA0003446056100000142
其中，φ为X轴到X'的转角；
t为首次碰撞至第二次碰撞之间的时间间隔；
L'为弹体质心与弹头顶部的距离；
vc为首次撞击前的弹体的质心速度；
第四步，经过坐标变换，得到弹头顶点的绝对运动轨迹(x，y)，其表达式如下所示：
Figure BDA0003446056100000151
其中，t为首次碰撞至第二次碰撞之间的时间间隔；
L'为弹体质心与弹头顶部的距离；
vc为首次撞击前的弹体的质心速度；
ω'2为首次碰撞后弹体角速度；
θ为弹体的入射方向与首次撞击面之间的夹角；
步骤十：根据异形靶板的几何尺寸，确定弹体与首次碰撞的异形靶板相邻的异形靶板的碰撞面上的碰撞点B的位置；
确定二次撞击点B的位置的方法,具体包括如下步骤：
A)确定第一次撞击点的位置A；
优选地，所述第一次撞击点A与异形靶板的顶点M之间的距离MA为a，所述第一次撞击点A与异形靶板的侧棱部之间的距离AN为b，如图1所示；
B)添加异形靶板的边线；
优选地，首次碰撞的异形靶板的边线MN和第二次碰撞的异形靶板的边线PN，
C)得到第二次碰撞的撞击面方程(xPN，yPN)，
yPN-b sinβ-L sin(θ+β)＝tan(180°-β-2θ)[xPN-b cosβ-L cos(θ+β)] (15)
其中，a为所述第一次撞击点A与异形靶板的顶点M之间的距离；
b为所述第一次撞击点A与异形靶板的侧棱部之间的距离；
θ为弹体的入射方向与首次撞击面之间的夹角；
β为首次撞击后弹体质心速度为v'c所在的方向与首次撞击的碰撞面MN的夹角；
L为弹体的长度。
D)当发生第二次碰撞时，(xPN，yPN)与(x，y)重合，得到第二次撞击的位置，也就是二次撞击点B的位置。
优选地，各异形靶板垂直于其长度方向的截面也可以为五边形，所述异形靶板包括第一侧壁、第二侧壁、第三侧壁、第四侧壁和底面；所述第一侧壁和所述第二侧壁形成异形靶板的顶部，该顶部为第一顶部；所述第一顶部均为尖顶部；所述第三侧壁与所述第四侧壁分别与所述底面连接，所述第一侧壁与所述第三侧壁的连接处形成第一侧棱部，所述第二侧壁与所述第四侧壁的连接处形成第二侧棱部；安装时，相邻两个防护单元中第一个防护单元的第一侧棱部与第二防护单元的第二侧棱部对齐。
3二次碰撞防弹障碍物设置探索
如图8中所示，假设首次导偏结构的尺寸已知，弹体首次撞击位置已知，则有MA＝a,AN＝b。现在需要设置二次碰撞防弹障碍物或二次碰撞防弹结构物。假设二次碰撞防弹障碍物或二次碰撞防弹结构物为直线结构PG，其与水平X轴夹角为
Figure BDA0003446056100000161
其中
Figure BDA0003446056100000162
与X轴交点为K，弹体入射点A与K之间的距离AK＝s1，则可以得到直线PG的方程：
Figure BDA0003446056100000163
进而可得，
Figure BDA0003446056100000164
得到t与
Figure BDA0003446056100000165
s1的关系式，简写为
Figure BDA0003446056100000166
根据
Figure BDA0003446056100000167
s1的取值规定进行试算，得到t，可求出的解为B点位置坐标。如果二次碰撞防弹障碍物或二次碰撞防弹结构物PK的方程与弹头顶点的绝对运动轨迹方程不存在交点，则需调整s1、
Figure BDA0003446056100000168
二次碰撞防弹障碍物或二次碰撞防弹结构物高度PK，直至方程有解。如果二次碰撞防弹障碍物或二次碰撞防弹结构物PK的方程与弹头顶点的绝对运动轨迹方程存在交点；那么，B点二次撞击前的速度是由与质心C速度大小和方向一致的
Figure BDA0003446056100000171
以及垂直于C1A1的转动速度
Figure BDA0003446056100000172
如果二次碰撞防弹障碍物或二次碰撞防弹结构物PK的方程与弹头顶点的绝对运动轨迹方程存在一个交点，则该交点即为二次撞击点B(xB，yB)。如果二次碰撞防弹障碍物或二次碰撞防弹结构物PK的方程与弹头顶点的绝对运动轨迹方程存在多个交点，则弹头顶点的绝对运动轨迹与二次碰撞防弹障碍物或二次碰撞防弹结构物的首次接触点即为二次撞击点B(xB，yB)。
试算过程中，优先选择
Figure BDA0003446056100000173
从90度开始。同时试验中有可能入射角出现偏差，为了保证安全，应该是的G点位于AN的延长线上。假设KG的长度为s2，那么，可得如下关系式：
Figure BDA0003446056100000174
应设置PB＝2BG以上，由此防弹结构的整体尺寸PG可以的得4算例假设弹径d＝30mm，弹速v＝400m/s，弹体长度l＝0.12m，弹体质心位于弹体中心，即l’＝0.06m。弹体质量m＝0.466kg，e＝0.2，θ＝60度。
则ω2＝5938.5度/秒，β＝-28.56°
则弹头顶点绝对运动方程：
Figure BDA0003446056100000175
根据试验现场设置，弹体与靶板的碰撞点和测孔板之间的水平距离为51.25cm。且测孔板垂直于靶板入射面。假设测孔板为直线，则改直线方程为：
y＝1.8372*x-1.1612
通过计算可得弹头轨迹曲线与测孔板方程的交点为R(0.6709，0.0714)。图9中，∠GAR即为弹体的偏转角，其等于34.66度。
图10、图11、图12为现场试验图，弹速分别为398m/s和411m/s，弹体与靶板的碰撞点和测孔板之间的水平距离平均值为51.25cm。最终测得弹体偏转角1在水平面上的投影分别为34.22度和33.94度。误差分别为1.27％和2.08％，与理论计算结果基本一致。
可以理解的是，以上实施方式仅仅是为了说明本发明的原理而采用的示例性实施方式，然而本发明并不局限于此。本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或全部技术特征进行等同替换；而这些修改或替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
在本发明的描述中，需要说明的是，术语“上”、“下”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性。
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“连通”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接连通，也可以通过中间媒介间接连通，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。此外，在本发明的描述中，除非另有说明，“至少三个”的含义是两个或两个以上。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种延时引信弹药与跳弹式防护结构碰撞分析方法，其特征在于，其具体包括如下步骤：
步骤一：预设初始飞行状态下所述弹体的攻角为零；
步骤二：建立第一次碰撞模型；
步骤三：基于第一次碰撞模型得到碰撞后弹头的速度v'A以及首次碰撞后弹头速度，确定首次碰撞后弹头速度在Y轴方向的分量投影v'Ay；
步骤四：确定材料恢复系数e；
步骤五：计算首次碰撞后弹体质心速度在Y方向上的分量v'cy；
步骤六：基于首次碰撞后弹体质心速度在Y方向上的分量v'cy，得到首次碰撞后弹体的角速度ω2以及首次撞击后弹体质心速度v'c所在的方向与首次撞击的碰撞面MN的夹角β；
步骤七：对弹体的质心位置进行修正；修正后，其弹体质心与弹头顶部的距离L'；
步骤八：得到修正后的首次碰撞后弹体角速度ω'2；
步骤九：根据首次碰撞后的修正弹体角速度ω'2以及首次撞击后弹体质心速度为v'c所在的方向与首次撞击的碰撞面MN的夹角β，建立第二次碰撞模型；
步骤十：根据异形靶板的几何尺寸，确定弹体与首次碰撞的异形靶板相邻的异形靶板的碰撞面上的碰撞点B的位置。
2.如权利要求1所述的延时引信弹药与跳弹式防护结构碰撞分析方法，其特征在于，步骤二中，异形靶板的第一尖顶部的两个侧面为异形靶板的撞击面，通过所述撞击面进行导偏，根据撞击前弹体的质心速度、撞击面上的弹着点，弹体入射角度建立第一次碰撞模型。
3.如权利要求2所述的延时引信弹药与跳弹式防护结构碰撞分析方法，其特征在于，步骤三中，所述首次碰撞后弹头速度在Y轴方向的分量投影v'Ay带动所述弹体进行偏转。
4.如权利要求3所述的延时引信弹药与跳弹式防护结构碰撞分析方法，其特征在于，步骤四中，确定材料恢复系数e，其表达式如下；
Figure FDA0003446056090000021
式中，v'Ay为首次碰撞后弹头速度在Y轴方向的分量投影；
vAy为初始状态下弹头速度沿着Y轴方向的分量投影；
vc为首次撞击前的弹体的质心速度；
θ为弹体的入射方向与首次撞击面之间的夹角。
5.如权利要求4所述的延时引信弹药与跳弹式防护结构碰撞分析方法，其特征在于，首次碰撞建模的初始坐标系为设表面MN方向为X轴，垂直于MN的方向Y轴，在初始坐标系中，所述弹着点为A，所述弹体的质心为C，撞击面为MN所在平面，弹体的入射方向与首次撞击面之间的夹角为θ。
6.如权利要求5所述的延时引信弹药与跳弹式防护结构碰撞分析方法，其特征在于，弹体在碰撞之前角速度ω1＝0，弹体首次碰撞时弹体受到Y轴方向的碰撞冲量Iy，Ix＝0，首次碰撞后弹体质心速度在X轴方向的分量与首次撞击前弹体质心速度在X轴方向的分量相等，对首次碰撞模型进行修正。
7.如权利要求6所述的延时引信弹药与跳弹式防护结构碰撞分析方法，其特征在于，第二次碰撞建模中，预设直线AD与首次碰撞后弹体质心速度v'c的方向平行；以直线AD为X轴，垂直于直线AD的方向为Y轴，重新建立第二次碰撞模型的定坐标系；同时，以t时刻后的弹体质心C1为原点建立动坐标系；其中t为首次碰撞至第二次碰撞之间的时间间隔；根据弹头顶点的相对运动轨迹(x',y')和建立弹头顶点的牵连运动轨迹(xo',yo')弹头顶点的绝对运动轨迹(x，y)。
8.如权利要求7所述的延时引信弹药与跳弹式防护结构碰撞分析方法，其特征在于，添加首次碰撞的异形靶板的边线MN和第二次碰撞的异形靶板的边线PN；得到第二次碰撞的撞击面方程(xPN，yPN)。
9.如权利要求6所述的延时引信弹药与跳弹式防护结构碰撞分析方法，其特征在于，二次碰撞防弹障碍物或二次碰撞防弹结构物为直线结构PG，其与水平X轴夹角为
Figure FDA0003446056090000031
其与X轴交点为K，弹体入射点A与K之间的距离AK＝s1。
CN202111649114.5A 2021-12-30 2021-12-30 延时引信弹药与跳弹式防护结构碰撞分析方法 Active CN114297860B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111649114.5A CN114297860B (zh) 2021-12-30 2021-12-30 延时引信弹药与跳弹式防护结构碰撞分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111649114.5A CN114297860B (zh) 2021-12-30 2021-12-30 延时引信弹药与跳弹式防护结构碰撞分析方法

Publications (2)

Publication Number Publication Date
CN114297860A true CN114297860A (zh) 2022-04-08
CN114297860B CN114297860B (zh) 2022-11-22

Family

ID=80972603

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111649114.5A Active CN114297860B (zh) 2021-12-30 2021-12-30 延时引信弹药与跳弹式防护结构碰撞分析方法

Country Status (1)

Country Link
CN (1) CN114297860B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115183632A (zh) * 2022-08-18 2022-10-14 北京汽车集团越野车有限公司 B柱防弹结构和车辆

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004024A (zh) * 2010-10-22 2011-04-06 中国工程物理研究院总体工程研究所 一种用于实现弹体倾斜撞击的模型实验系统
US20140358505A1 (en) * 2013-05-31 2014-12-04 The Board Of Trustees Of The University Of Illinois Collision impulse derived discrete element contact force determination engine, method, software and system
CN108549733A (zh) * 2017-12-22 2018-09-18 彩虹无人机科技有限公司 一种导弹在轨运动的动力学建模方法
CN109145386A (zh) * 2018-07-24 2019-01-04 张瑞杰 一种用于工程结构的地震碰撞反应分析方法
CN111553019A (zh) * 2020-04-20 2020-08-18 哈尔滨工程大学 一种回转体撞击船体板架横梁剩余速度与转角计算方法
US20210102781A1 (en) * 2018-03-26 2021-04-08 Korea Military Academy R&Db Foundation Point-of-impact analysis apparatus for improving accuracy of ballistic trajectory and point of impact by applying shooting environment of real personal firearm to virtual reality, and virtual shooting training simulation using same

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004024A (zh) * 2010-10-22 2011-04-06 中国工程物理研究院总体工程研究所 一种用于实现弹体倾斜撞击的模型实验系统
US20140358505A1 (en) * 2013-05-31 2014-12-04 The Board Of Trustees Of The University Of Illinois Collision impulse derived discrete element contact force determination engine, method, software and system
CN108549733A (zh) * 2017-12-22 2018-09-18 彩虹无人机科技有限公司 一种导弹在轨运动的动力学建模方法
US20210102781A1 (en) * 2018-03-26 2021-04-08 Korea Military Academy R&Db Foundation Point-of-impact analysis apparatus for improving accuracy of ballistic trajectory and point of impact by applying shooting environment of real personal firearm to virtual reality, and virtual shooting training simulation using same
CN109145386A (zh) * 2018-07-24 2019-01-04 张瑞杰 一种用于工程结构的地震碰撞反应分析方法
CN111553019A (zh) * 2020-04-20 2020-08-18 哈尔滨工程大学 一种回转体撞击船体板架横梁剩余速度与转角计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
秦有权: "矢量防护工程理论探讨", 《防护工程》 *
韩永要等: "长管体斜侵彻半无限均质靶简化模型", 《工程力学》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115183632A (zh) * 2022-08-18 2022-10-14 北京汽车集团越野车有限公司 B柱防弹结构和车辆

Also Published As

Publication number Publication date
CN114297860B (zh) 2022-11-22

Similar Documents

Publication Publication Date Title
KR20110126770A (ko) 피탄충격 시험용 충격탄 및 이를 이용한 피탄충격 모의 시험장치
US8563910B2 (en) Systems and methods for targeting a projectile payload
EP2040025A1 (en) Shooting target system for automatic determination of point of impact
US8468926B2 (en) Ballistic armor system
CN114297860B (zh) 延时引信弹药与跳弹式防护结构碰撞分析方法
CN106052491A (zh) 一种模拟爆炸冲击波和高速破片群联合载荷作用试验方法
CN112035780B (zh) 一种导弹末制导阶段杀伤效果计算方法
CN113406966A (zh) 一种无人机的反制方法及无人机的反制系统
US8863666B2 (en) Method and system for electronically shaping detonated charges
CN114297871B (zh) 基于弹靶斜碰的跳弹轨迹预测模型
Ji et al. Dynamic behaviors of multi-layered steel targets with air gaps subjected to the impact of EFP simulants
EP3615884B1 (en) Projectile with selectable angle of attack
CN114297861B (zh) 弹体与可动类防弹障碍物碰撞分析方法
CN116797032A (zh) 一种坦克打击装甲类目标的毁伤效能评估方法
RU2294526C1 (ru) Способ определения дальности срабатывания неконтактного взрывателя при проведении испытаний управляемого снаряда
Żochowski et al. Numerical analysis of effectiveness for vehicle net systems protecting against shaped charge projectiles
CN113405411A (zh) 一种战斗部破片质量分布的试验方法及系统
KR101925690B1 (ko) 기갑차량 방어장치
Jiang et al. Dynamic multiprojectile attack and killing effects of detonation warheads
RU2315939C1 (ru) Способ наведения на цель ракет, управляемых по лучу
RU2230284C2 (ru) Кассетный снаряд "хорс"
RU2237230C1 (ru) Осколочный снаряд направленного действия "стрибог"
Zhang et al. Research progress of anti-penetration yaw technology for concrete protective structures
JPH07159095A (ja) 射撃模擬装置
Shi et al. The Research on Damage Effect of Target under the Action of Fragmentation Warhead

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant