1. 引言
近年来随着钻地武器研制和应用的迅速发展,钻地武器的侵彻深度研究得到了广泛的关注。目前世界各国针对弹体侵彻问题做了大量的试验研究,研究最多的就是关于钢筋混凝土结构的高速侵彻,而关于土介质的研究却很少,一个是因为应用的场合比较少,另一个是因为侵彻过程中,涉及到大量材料的非弹性形变、裂缝扩展等一系列难以用数学工具准确描述的现象,导致目前还没有比较完善的模型。
对于弹体侵彻土介质靶体的深度研究,可以从分析侵彻阻力入手。影响侵彻阻力的介质性质方面的因素有:介质密度
,屈服强度
,弹性模量
,粘滞系数
。弹体结构因素有:弹体质量
,弹体直径
,弹头部卵形直径
。弹体侵彻的运动条件因素有:弹速
,偏转角
,以及侵彻深度
。则弹体侵彻土介质的阻力可表示为:
(1)
实际上,由于侵彻机理的复杂性和试验、计算工具、方法的简陋落后,并因为在建立模型时做了大量的假设和简化处理,基于空腔膨胀理论的纯力学分析很难给出完全符合实际情况的、符合试验结论的方程。于是,一些人偏向于基于有限元、有限差分、离散元等的数值模拟 [1] 。目前,商业有限元程序在侵彻和高速碰撞领域的应用非常普遍 [2] ,然而由于数值模拟结果的准确性在很大程度上取决于材料性质的定义,计算结果往往与实际情况不相符。
另外一些研究人员从试验出发,通过大量的试验数据进行分析,提出了一些经验公式。其最有代表性的有:美国Sandia National Laboratories提出的Young公式 [3] ,Young公式需要通过试验不断修正相关参数,因此不具有普遍适用性,且无法对侵彻过程各物理量的相关性进行准确分析,具有很大的局限性 [4] ;Petry公式 [5] [6] 是在Poncelet公式基础上发展起来的,是一种半经验、半分析公式,它是量纲不相符的 [7] ;M.J. Forrestal等人在考虑靶体材料性质、弹头形状因素等对侵彻过程的影响情况下,提出了Forrestal侵彻方程 [8] 。由于每个经验公式都有各自的应用范围和应用条件,最终的计算结果也存在一定的差别。
本文将从力学方面进行分析,列出侵彻阻力的表达式,并将表达式的量纲考虑进去,利用量纲分析法对低速弹体撞击土靶体试验进行分析。通过量纲分析既可以把多变量的函数关系归纳为数量最少而又能反映基本关系的无量纲化的函数关系,从而为试验研究提供了方便的条件和内在本质性的规律,得出对撞击过程的全面描述。
2. 力学分析——阻力公式的困惑
根据Poncelet阻力公式,将弹体侵彻介质过程的阻力分为两个部分 [9] :(1) 固体介质对弹体侵彻的静阻力,也就是弹体在固体介质中运动时,以自身的能量去克服的侵彻孔附近介质分子之间联系的凝聚力和弹体与介质之间的摩擦力。静阻力的大小与弹体的横截面积成正比,与弹体的侵彻速度无关;(2) 介质对弹体的动阻力,即弹体在介质中运动时介质对弹体的阻力。动阻力的大小不仅与弹体的横截面积有关,还与弹体的侵彻速度平方成正比。这里采用Poncelet阻力公式的上述思想,建立弹体的侵彻阻力公式:
(2)
分析可知,
的量纲为MLT−2,
的量纲为L2,
的量纲为MLT−2,所以公式(2)等号两边的量纲不相等,是非其次的,而能够正确反映侵彻过程这一现象规律的方程,其等号两端都应该具有相同的量纲,所以下面将通过量纲分析与建模正确分析阻力与各因素变量之间的关系。
3. 量纲分析及建模
当对一个物理概念进行定量描述时,总离不开它的一些特性,比如时间、质量、密度、速度等等,这种表示不同物理特性的量,称之为具有不同的“量纲” [10] 。在物理学中,一个物理量的量纲是质量、时间和长度的结合,分别由符号M、T、L表示,任意一个量均可按一定的规则由这三个基本量导出。例如,物理量密度的量纲是质量/长度/长度/长度(ML−3)。
任何有意义的定律,对于其方程式的每一个计量单位,都必须是齐次方程式,这个认识的最终形式成为Õ定理 [11] 。白金汉Õ定理(Buckingham π theorem)表明,对于任何一个有物理意义的方程式,如果存在
个变量,其中有
个基本量,则存在
个独立的无量纲参数,即可以将
个变量组合成
个无量纲Õ数。
弹体对土介质的撞击过程是一个复杂的物理过程,影响的因素很多,很难通过单独的数学理论分析或实验的方法得出一个精确的解,利用量纲分析法研究弹体侵彻深度问题,是一种有效且可靠性较高的方法。这里基本量纲单位采用MLT国际单位制,表1列出第一小节11个参量的量纲。
根据Õ定理,取
,并取
为量纲独立变量,则在(1)式中可以找到8个无量纲的物理量,以Õ代表无量纲量,则有:
(3)
式中:
为待定系数,因为
是无量纲量,所以分子分母量纲应该相等,按照
和
的量纲可得下式:
(4)
对于
有:![](//html.hanspub.org/file/6-2950239x38_hanspub.png)
对于
有:![](//html.hanspub.org/file/6-2950239x40_hanspub.png)
对于
有:![](//html.hanspub.org/file/6-2950239x42_hanspub.png)
解得![](//html.hanspub.org/file/6-2950239x43_hanspub.png)
则:
(5)
同理,可以得到其余7个参量的无量纲表达形式。
则(1)式的无量纲关系式为:
(6)
函数
中的前三项与介质的性质有关,后三项与介质的性质无关,因此可以简化成:
![](Images/Table_Tmp.jpg)
Table 1. The dimensional parameters of projectile and soil
表1. 弹体、土的各参量量纲
(7)
式中:
:弹体横截面密度与靶体材料密度比值影响函数;
:弹体头部形状影响函数;
:弹体侵彻深度影响函数;
:弹体偏转角影响函数;
将式(1)写成:
(8)
其中:
,
。
采用上一小节阻力公式的上述思想,结合无量纲分析,建立弹体的侵彻阻力公式:
(9)
式中,
——静抗力系数;
——动抗力系数。
由运动方程得:
(10)
令:
,![](//html.hanspub.org/file/6-2950239x71_hanspub.png)
则:
(11)
因为在弹体侵彻过程中,侵彻速度是变化的,所以有:
(12)
把(12)式代入(11)式,消去
,得:
(13)
对上式两边积分,即:
(14)
从而:
(15)
上式中的
通过实验而确定。
4. 试验
试验采用质量0.48 kg,弹径30 mm,CRH = 3的42CrMo铜材料弹体进行侵彻,土介质靶体的相关参数如下:密度
、抗压强度
、弹性模量
、粘滞系数
。
土中侵彻试验结果采用文献 [12] 中的数据,如表2所示。
利用表2中的数据进行拟合,确定公式(15)中的两个系数的值,即
,
。则公式(15)可以写成:
(16)
利用(16)式计算结果如表3所示。
从上表可以看出,当弹速低于200 m/s时,用本文方法预估的深度和试验结果相近,而当弹速高于200 m/s时,误差将达40%,这是因为当弹体高速侵彻时,弹体头部会开始钝化,从而严重变形。由于存在弹体材料的物理侵蚀现象,在此对公式(16)进行修正,乘以一个靶体修正因子,此因子与弹形、弹速无关,只决定于靶体的成分与堆砌工艺。修正方法如下:
(17)
式中,
,是无量纲相对速度的函数。
利用表3数据对k进行拟合,得出修正后的深度预估公式为:
(18)
表4为应用Petry公式、Ferrostal公式、Young公式及本文提出的公式预估弹体侵彻土介质靶体的深度结果对比。从表中可以看出,各个公式的计算结果互不相同,与试验结果也有一定的差异,但本文提出的方法预估效果比较好,误差小。
5. 总结
本文介绍了目前国内外用于混凝土等介质侵彻深度计算的经验公式,在此基础上通过试验并结合国内外试验数据,利用量纲分析和力学分析方法,得到了一个新的关于侵彻土介质深度的预估公式。最后
![](Images/Table_Tmp.jpg)
Table 2. Penetration test results with soil medium
表2. 土介质侵彻试验结果
![](Images/Table_Tmp.jpg)
Table 3. Comparison of result with the method proposed in this paper
表3. 本文方法数据处理结果对比
![](Images/Table_Tmp.jpg)
Table 4. Comparison of computed result with all kinds of method
表4. 各种方法计算结果对比
的数据处理结果表明,我们可在靶场数据处理的基础上,建立自己包含靶体修正因子在内的一个侵彻深度公式。该公式不仅适用于沙土层,调整靶体修正因子即可适用于混凝土等其他介质。