• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

小隔板试验冲击起爆和熄灭流场的实验与Lagrange分析

蒋国平 浣石

孔祥韶, 杨豹, 周沪, 郑成, 刘芳, 吴卫国. 基于响应面法的纤维金属层合板抗弹性能优化设计[J]. 爆炸与冲击, 2022, 42(4): 043301. doi: 10.11883/bzycj-2021-0146
引用本文: 蒋国平, 浣石. 小隔板试验冲击起爆和熄灭流场的实验与Lagrange分析[J]. 爆炸与冲击, 2005, 25(5): 472-476. doi: 10.11883/1001-1455(2005)05-0472-05
KONG Xiangshao, YANG Bao, ZHOU Hu, ZHENG Cheng, LIU Fang, WU Weiguo. Optimal design of ballistic performance of fiber-metal laminates based on the response surface method[J]. Explosion And Shock Waves, 2022, 42(4): 043301. doi: 10.11883/bzycj-2021-0146
Citation: JIANG Guo-ping, HUAN Shi. Experimental investigation and Lagrange method used for small partition assay[J]. Explosion And Shock Waves, 2005, 25(5): 472-476. doi: 10.11883/1001-1455(2005)05-0472-05

小隔板试验冲击起爆和熄灭流场的实验与Lagrange分析

doi: 10.11883/1001-1455(2005)05-0472-05

Experimental investigation and Lagrange method used for small partition assay

  • 摘要: 根据小隔板试验的拉氏传感器测试结果,用拉氏分析方法求得了粗颗粒包复TNT炸药中的冲击波成长过程和熄灭过程的流场分布,用实验结果确定了点火增长模型中的待定参数。利用拉氏分析方法时,对未反应炸药和已反应产物的状态方程都采用JWL状态方程,使确定的参数可以直接用于有限元程序LSDYAN计算。通过拉氏分析得到了各迹线上的反应度,结果表明,冲击波在熄灭过程中仍能引起炸药的局部化学反应。
  • 高速侵彻过程中,弹丸与靶板之间会产生强冲击载荷,形成复杂的高温、高压和高应变率受力环境,这会导致弹体发生破碎断裂,进而严重影响弹体的侵彻能力[1],因此学者们对高速侵彻过程中弹体的破碎断裂问题进行了相关研究[2-10]。Rakvag等[3-4]研究了动能弹高速侵彻钢制装甲过程中碎片的形成和不同速度下Taylor杆的断裂失效模式,Jones等[5]和Hiermaier等[6]详细分析了实体与壳体等在冲击下的结构屈曲形态与破坏响应,李硕等[1]和肖新科[7]分别开展了35CrMnSi和38CrSi合金钢材料的弹体失效与断裂行为的研究。这些相关研究中主要集中于实心杆弹对金属靶板侵彻的断裂问题,而对空心或带模拟装药的空心弹体结构研究较为缺乏,尤其是在高速侵彻岩石靶板的情况下[8-11]

    为了研究带模拟装药弹体高速侵彻岩石靶板的破坏机理,本文设计了两种不同壁厚的弹体,进行了着速约3Ma的现场岩石侵彻试验研究。在试验的基础上,考虑到弹体发生的完全破碎,通过Autodyn-3D软件对弹体采用SPH算法和Mott失效模型进行数值模拟,分析了弹药结构的破坏过程和机理,并且讨论了小范围内不同的高着速对弹体破坏的影响。

    设计了2种不同壁厚的弹体,结构如图 1所示,其直径80 mm、长394 mm,而壁厚h分别为17和22 mm,对应质量分别为8.35和9.95 kg。弹体材料为35CrMnSi,抗拉强度约1.5 GPa。使用125 mm口径滑膛炮发射次口径弹丸撞击现场流纹岩靶体,共发射2种不同壁厚的试验弹各2枚,速度约1 000 m/s。

    图  1  弹体结构示意图
    Figure  1.  Sketch of projectile geometry

    图 2所示为试验现场的岩石靶破坏形态,可以发现岩石靶体表面在弹体高速冲击后形成明显的放射性裂纹,靶表面岩石呈较为规则的块体脱落。然而,弹体仅在靶表面撞击后形成粉碎性破坏,未能有效侵入靶体。

    图  2  弹体冲击后岩石靶体正表面
    Figure  2.  Rock target's impacted surface

    试验中,2种不同壁厚的各2发试验弹的弹体结构均发生破碎,散布于试验现场,收集到的弹片很少,如图 3所示。从图 3中可以看出,收集到的较大尺寸的弹片主要是弹体靠后的部分,尤其是弹尾(带螺纹),而越靠近弹头部分收集到的弹片越小,破碎的越严重。另外,弹尾部收集到的破片几乎可拼凑出完整的弹尾,这意味着弹壳体后段呈现四瓣或更多瓣裂开。此外未收集到其他大尺寸的弹尖头位置的碎片,可以认为弹头部由于严重塑性变形而破碎。图 3中第2、3列破片表现出了较为明显的外翻形态,证明该段壳体发生了径向位移。

    图  3  试验现场收集的弹片
    Figure  3.  Collection of projectile fragments in field test

    此外,从图 4中弹片断裂面可以看出以下两点:壳体轴向断裂面不够规则,呈现不稳定拉伸破坏;壳体环向断裂面大体上呈现45°破坏,可以断定为剪切破坏。

    图  4  两种典型破坏断面
    Figure  4.  Two typical broken sections

    值得注意的是,试验中弹体头部发生完全破碎(图 3中第4列),与文献[6-8]中薄壁弹体撞击有限厚金属靶后呈现弹头完整现象有明显区别,这是因为后者弹头与靶作用时间较短,在严重变形之前就已经穿透;试验中,靶体为半无限厚,岩石一直挤压甚至刨蚀弹头,导致弹头部破碎而无法保持完整。

    根据图 1给出的弹体几何尺寸,在Truegrid 3D中建立模型并导入Autodyn软件中。岩石靶取为圆形靶,直径1 200 mm,厚800 mm。对整个弹体(弹壳、后盖、填充物)均使用SPH方法进行离散,粒子尺寸取为2 mm,粒子数目约为100 000;对靶体使用有限元网格离散,靶体中心网格尺寸为4 mm,向外比例扩大;此外为了减小边界效应,使用20 mm厚钢圈围住岩石靶,有限元单元总数约为650 000。由于侵彻过程中弹体材料随机失效,呈现不对称性,弹体可能发生偏转,所以需要建立二分之一模型,在软件中建立的模型如图 5所示。

    图  5  计算模型
    Figure  5.  Simulation model
    2.2.1   金属材料

    使用Johnson-Cook本构模型[12]来描述弹体、铝后盖及钢圈等金属材料的屈服强度:

    σy=(A+Bεn)(1+Cln˙ε˙ε0)[1(TT0TmT0)m]
    (1)

    式中:ABnCm分别为屈服强度、硬化强度、应变硬化系数、应变率硬化系数、温度系数,可通过实验测定;σyε˙ε˙ε0对应屈服应力、等效应变、等效应变率和参考应变率;T0Tm分别为参考温度和材料融化温度,T为材料瞬时温度。材料参数如表 1所示,其中ρ0G分别为密度和剪切模量。

    表  1  金属材料主要参数
    Table  1.  Main parameters of metal material
    位置 材料 ρ0/(kg·m-3) G/GPa A/GPa B/GPa n C
    弹体 35CrMnSi 7.93 81.8 1.500 0.500 0.26 0.014
    后盖 7039铝 2.77 27.6 0.337 0.345 0.41 0.010
    钢圈 4340钢 7.83 81.8 0.792 0.510 0.26 0.014
    下载: 导出CSV 
    | 显示表格

    从试验结果分析可以看出,弹体侵彻岩石靶后呈现了破碎现象,弹片较多,所以在模拟时有必要考虑弹材的随机破坏特性。Mott随机失效模型[12-14]是依托大量的试验数据而建立起来的经验性模型,具有广泛的适用性,并且在Autodyn计算软件中可以直接应用,是当前模拟金属材料随机失效的较多的方法。Mott模型选择材料弱化点破坏概率在[0, 1]范围内,对不同的失效塑性应变取概率分布:

    P=1exp(Dγeγε)
    (2)

    式中:P为比例塑性应变为ε*时的单元失效概率,Dγ取决于材料性质的常数。γ值越大,材料均匀性越强。

    对弹材添加Mott失效模型,失效模式为主应力失效,失效阈值设为屈服强度(1.5 GPa),取γ=10,设定随机失效从破坏应力的50%开始[13-14]

    2.2.2   岩石

    JH-2模型[15-16]是目前模拟脆性类材料使用较为广泛的一种本构模型,该模型考虑了压力、应变率和损伤对材料强度的影响,可以在AUTODYN软件材料库中直接使用[17]。JH-2模型针对原始HJC模型[18-19]无法反应材料的软化特性和部分参数无法准确得到的缺点进行了改进[15],其强度模型的无量纲形式为:

    ˉσ=ˉσiD(ˉσiˉσf)
    (3)
    ˉσi=A(ˉp+ˉT)N[1+Cln(˙ε/˙ε0)]
    (4)
    ˉσf=BˉpM[1+Cln(˙ε/˙ε0)]
    (5)

    式中:ABCMN均为材料常数,σσiσf分别为无量纲等效应力、无损伤的材料强度(D=0)和破碎的材料强度(D=1),p为无量纲静水压,T为无量纲的静水拉伸强度,D为材料损伤(0≤D≤1)。

    考虑到:(1)流纹岩和花岗岩的成分基本一致,而流纹岩的HJC参数研究较少,花岗岩的研究比较丰富;(2)靶体流纹岩取样实测密度为2.66 g/cm3,单轴抗压强度178 MPa,这与文献[16]试验所用的巴利花岗岩(Barre granite)的密度(2.66 g/cm3)和单轴抗压强度(167.1 MPa)均较为相近;(3)着重探讨弹体的破坏而非高强度的靶体响应。故本文中流纹岩石JH-2模型参数重点参考文献[16]。

    2.2.3   填充物质(硫磺)

    试验中弹体内部使用硫磺作为模拟装药进行填充,Autodyn软件材料库[17]自带硫磺参数,如表3所示,使用SHOCK状态方程,无强度方程。

    表  2  硫磺材料参数
    Table  2.  Material parameters of sulfur
    材料 ρ0/(kg·m-3) Grüneisen系数 c/(m·s-1) S1
    硫磺 2.02 0 2.7 0.95
    下载: 导出CSV 
    | 显示表格
    2.3.1   壳体失效的机理分析

    图 6给出了弹体高速冲击岩石(着速1 000 m/s)的侵彻过程,展现了弹体变形的过程以及岩石靶破坏及裂纹扩展的过程,与现场试验中(图 2)岩石靶表面放射性裂纹现象吻合。从图 6中可以看出,侵彻过程中(t=300 μs)由于应力集中导致弹体薄弱位置塑性屈服,弹体头肩部位置首先发生了明显的径向膨胀随后断裂,随着侵彻过程的继续(t=300 μs),弹体圆柱段径向扩张加剧。

    图  6  弹体侵彻岩石靶过程
    Figure  6.  Simulation of penetration into rock

    图 7单独给出了弹壳的破坏过程:由于弹体的着速高和岩石靶体的强度高,弹头表面材料破坏严重,同时弹头与后段壳体过渡处由于其壁厚较小,承受的应力超出了材料的屈服极限,进而失效直至断裂。分析原因如下:内部装填的硫磺密度低,强度低,受到了靶体与破碎弹头的高度挤压,使得壳体径向膨胀并产生了明显的裂纹(见图 8),最后壳体沿轴向及环向破坏形成了外翻撕裂型破片(t=700 μs)。试验收集到的在尾部的破片也呈现四瓣(及以上)裂开模式,图中模拟结果明显可见的花瓣数为4个,与试验结果基本相符。

    图  7  弹壳破碎过程
    Figure  7.  Projectile fracture process
    图  8  t=300 μs时刻弹壳轴向裂纹
    Figure  8.  Axial crack in projectile (t=300 μs)

    计算表明本文对弹体使用的Mott随机失效模型及参数具有较高的可信性,本文中数值模拟可以真实地还原弹高速撞击岩石靶时的失效过程。

    2.3.2   装填物对弹体破坏的影响

    上小节分析中提到:壳体的轴向及径向断裂是由于内部装填物(模拟装药)挤压的作用结果,由于该分析不是通过对比试验得出,这里在其它条件不变的情况下仅去掉装填物来进行数值对比试验,验证提出的该分析并进一步探讨弹体的破碎机理。

    图 9给出了对比验证模拟,展示了无装填物时壳体的破坏过程,可以发现:没有装填物时,弹体在薄弱处断裂后变成近似圆柱壳体,在高速冲击下不断地外翻撕裂并伴随连续失效破坏,而没有表现出图 8中所示的明显裂纹,最终弹体破坏形态即为变短的筒体,也非花瓣形破片,这显然与试验中弹体最终状态不符。这进一步确认了试验中装填物对弹体破坏的重大影响。

    图  9  无装填物的弹壳破坏过程
    Figure  9.  Crush process of projectile without filling

    空心弹体内部是否含有装填物的对比数值模拟表明:(1)高速侵彻过程中,弹体发生破碎与是否携带装填物无关,但装填物对弹壳的破坏形态有着显著的影响;(2)进一步确认了弹体破坏是由于壳体头肩部薄弱处严重塑性变形至断裂,内部低阻抗装填物受到高度挤压进而使壳体径向膨胀所致。这与文献[8-10]中薄壁弹体撞击有限厚金属靶后呈现的动屈曲破坏模式存在着本质的区别。

    2.3.3   弹速对弹体破坏的影响

    图 10给出了不同弹速时的弹壳体最终破坏形态的比较,可以看出:800~1 000 m/s的着速下,弹壳体裂纹均发展到了壳体末端,着速越高壳体外翻撕裂破坏越严重。另外,800和1 000 m/s着速下,壳体呈现四瓣裂开,而900 m/s速度下壳体呈现六瓣裂开,该计算结果初步表明:800~1 000 m/s的着速下,壳体的裂开瓣数与速度大小及其变化没有直接相关性但花瓣式的裂开模式与试验结果比较一致。

    图  10  不同弹速下弹壳破坏比较
    Figure  10.  Comparison of projectile crush at different velocites

    进行着速约1 000 m/s的2种不同壁厚的弹药结构对高强度岩石靶的侵彻试验,试验表明:薄壁弹体高速冲击岩石靶后,弹头部分完全破碎,这与撞击有限厚金属靶后呈现的弹头完整现象有所区别。另外,弹体部分有明显的外翻撕裂和剪切破坏,弹尾则分为几乎完整的4瓣破片。

    在试验基础上,结合三维数值模拟分析弹药结构的破坏机理,进一步确认弹体破坏是由于壳体头肩部薄弱处严重塑性变形至断裂,内部低阻抗装填物受到高度挤压进而使壳体径向膨胀所致,结果表明:

    (1) 结合SPH和Mott随机失效的数值模型可以真实地还原弹高速撞击岩石靶时的失效过程,具有较高的可信度,可以对弹体极限设计的数值模拟提供参考;

    (2) 高速侵彻过程中,弹体发生破碎与是否携带装填物无关,但装填物对弹壳的破坏形态有着显著的影响;

    (3) 通过不同速度侵彻模拟初步表明:800~1 000 m/s速度下,弹体花瓣式的裂开模式与高速冲击岩石靶试验结果比较一致,但是裂开瓣数没有与速度的大小及其变化表现出直接的相关性。

    本文中的试验数据和具有可信度的数值模型可为进一步探讨高速侵彻岩石混凝土靶板的弹体的结构安全提供有效参考。

  • 期刊类型引用(6)

    1. 高建村,王乐,胡守涛,孙谞,周尚勇,王涛. 不同磁性金属丝对丙烷爆炸反应抑制机理研究. 中国安全生产科学技术. 2020(07): 125-130 . 百度学术
    2. 王连聪,梁运涛,罗海珠. 我国矿井热动力灾害理论研究进展与展望. 煤炭科学技术. 2018(07): 1-9 . 百度学术
    3. 路长,刘洋,于子凯,潘荣锟,刘磊,滕飞. 四氟乙烷对甲烷/空气爆炸特性的影响. 化工进展. 2017(10): 3596-3603 . 百度学术
    4. 徐永亮,王兰云,余明高,万少杰,梁栋林. 感应式荷电细水雾对瓦斯爆炸传播速度的影响. 中南大学学报(自然科学版). 2016(08): 2884-2890 . 百度学术
    5. 余明高,梁栋林,徐永亮,郑凯,纪文涛. 荷电细水雾抑制瓦斯爆炸实验研究. 煤炭学报. 2014(11): 2232-2238 . 百度学术
    6. 丁小勇,谭迎新,李媛,秦涧,许航. 管道倾斜角对瓦斯爆炸影响的实验研究. 消防科学与技术. 2012(12): 1269-1271 . 百度学术

    其他类型引用(14)

  • 加载中
计量
  • 文章访问数:  2671
  • HTML全文浏览量:  121
  • PDF下载量:  293
  • 被引次数: 20
出版历程
  • 刊出日期:  2005-09-25

目录

/

返回文章
返回