新型结构炮口制退器的膛口冲击波数值研究与性能分析

余海伟 袁军堂 汪振华 葛苗冉 罗跃

卢景瑞, 寇自力, 刘腾, 张雷雷, 丁未, 张强, 王强, 杨鸣, 龚红霞, 贺端威. 亚微米级聚晶金刚石的高温高压合成[J]. 高压物理学报, 2018, 32(2): 023302. doi: 10.11858/gywlxb.20170574
引用本文: 余海伟, 袁军堂, 汪振华, 葛苗冉, 罗跃. 新型结构炮口制退器的膛口冲击波数值研究与性能分析[J]. 高压物理学报, 2020, 34(6): 065102. doi: 10.11858/gywlxb.20200568
LU Jingrui, KOU Zili, LIU Teng, ZHANG Leilei, DING Wei, ZHANG Qiang, WANG Qiang, YANG Ming, GONG Hongxia, HE Duanwei. Sub-Micron Polycrystalline Diamond Synthesis under High Temperature and High Pressure[J]. Chinese Journal of High Pressure Physics, 2018, 32(2): 023302. doi: 10.11858/gywlxb.20170574
Citation: YU Haiwei, YUAN Juntang, WANG Zhenhua, GE Miaoran, LUO Yue. Muzzle Blast Wave Investigation and Performance Analysis of New-Structure Muzzle Brake Based on Numerical Simulation[J]. Chinese Journal of High Pressure Physics, 2020, 34(6): 065102. doi: 10.11858/gywlxb.20200568

新型结构炮口制退器的膛口冲击波数值研究与性能分析

doi: 10.11858/gywlxb.20200568
基金项目: 国防预研联合基金(6141B0104-2)
详细信息
    作者简介:

    余海伟(1990-),男,博士研究生,主要从事流体仿真与轻量化研究. E-mail:heoiwee@foxmail.com

    通讯作者:

    袁军堂(1962-),男,博士,教授,主要从事先进加工工艺与装备研究. E-mail:mc106@njust.edu.cn

  • 中图分类号: O347.5; TJ303.2

Muzzle Blast Wave Investigation and Performance Analysis of New-Structure Muzzle Brake Based on Numerical Simulation

  • 摘要: 面向炮口制退器综合性能设计需求,基于增材制造技术优势,提出了一种叠加冲击式内壁与反作用式外孔特征的新型小口径轻金属炮口制退器方案。基于三维无黏Euler方程,建立了后效期火药气体排空过程的有限元仿真模型,通过膛口流场仿真与流固耦合分析,研究了膛口流场发展、膛口冲击波与超压分布特征、制退器效率与其强度性能。结果表明:相比于传统结构,采用同种轻质钛合金材料的新型制退器具有较高的制退效率,对于机载平台具有较优的膛口超压分布特征和低冲击波危害效应,并且结构强度满足使用要求。

     

  • 聚晶金刚石(Polycrystalline Diamond, PCD)是人造金刚石的一种。利用高温高压条件合成的PCD通常分为两种:一种是含黏结剂的PCD; 另一种是不含黏结剂的纯相PCD。含黏结剂的PCD因需要的合成条件较低(5.0~6.0GPa、1600℃)而普遍应用于工业生产中,纯相PCD的合成则需要很高的温压条件(约15.0GPa、2300℃)。在研究探索新型超硬材料[1-9]的同时,人们也非常关注在工业上大量使用的PCD材料。自20世纪50年代以来,利用Co作为黏结剂,采用高温高压的方法合成的PCD已经广泛应用于工业生产中[10-14]。Hall[15]首次报道在压力为6.5GPa、温度为1726~2226℃条件下烧结得到了多晶金刚石。Katzman等[16]以1~5μm的金刚石微粉为原料,金属钴作为黏结剂,在压力为6.2GPa、温度为1600℃条件下烧结20min成功制得多晶金刚石材料。

    PCD在工业上的应用越来越广泛,其合成工艺也基本定型,人们已经习惯使用触媒在压力为5.0~5.5GPa、温度为1600℃的条件下左右进行PCD的合成,合成出的PCD硬度约为50GPa[17]。但是,目前市场上流通的合成PCD往往存在一些问题,例如:当含Co量不同时,烧结所得PCD的性能也不同; 添加的黏结剂容易聚集在一起,使合成出的PCD存在一定的缺陷; 直接混合黏结剂烧结时,存在混料不均匀的问题; 采用熔渗法合成PCD时,尺寸在微米以下的黏结剂不容易渗透进来等[18-19]。这些问题的存在在很大程度上影响了PCD产业的发展,再加上国外合成PCD的大公司对于技术的保密,使得PCD行业每研发一种新产品都要花费很长的时间,既浪费了大量的人力物力,又在一定程度上打击了整个行业的积极性。因此,迫切需要发展一种既能够满足对合成条件要求较低,又能够获得性能优异的PCD材料的方法。

    本研究采用分层组装[20]的方法对亚微米级PCD的合成方法进行探索,成功合成出性能优异的亚微米级PCD,并对合成PCD样品进行系统的表征和分析,以期获得既能保证合成样品优异性能、又具备可重复性的亚微米级PCD烧结条件。

    初始金刚石粉末购买自河南亚龙超硬材料有限公司,该粉末通过爆轰法制得,并经过多次粒度筛选。金刚石微粉的粒度分别为0.5、10μm,纯度为99.99%。基底选取WC-16%Co(16%为质量分数)。金刚石微粉在生产过程中虽然已经经过了提纯处理,但微粉不可避免地会与空气接触,因此颗粒表面可能还会附着金属触媒、油脂等物质。特别是对于爆轰法合成的金刚石微粉而言,由于生产过程中涉及到爆轰环节,会使产物中可能残留原料石墨,同时也会混入新的杂质。因此,对金刚石微粉原料进行净化处理是非常必要的。首先采用强酸(HF、HCl)处理,除去混入金刚石微粉原料中的金属、硅及其氧化物等杂质,并用去离子水洗涤至中性,然后通过烘箱将金刚石粉体烘干。此外,金刚石粉体表面会吸附大量有害气体,这些含有其他元素的气体对PCD的合成影响很大,可以采用真空高温的方法(真空度至4.0×10-3Pa,温度为800℃,并在800℃下保温1.5h后降温至室温)除去有害气体。

    对经过酸和真空处理的金刚石粉末进行X射线衍射(XRD)和扫描电镜(SEM)分析。采用丹东方圆DX2500型X射线衍射仪,工作电压40kV,工作电流25mA, Cu Kα辐射(λ=0.15406nm),步长为0.02°,采样时间设置为3~5s,扫描角度为10°~100°。图 1为金刚石微粉的XRD图,图 2为金刚石微粉的SEM图像。从图 1中可以看出,处理后的金刚石微粉展示出立方结构的金刚石相,其衍射峰与PDF卡片库中的立方结构金刚石PDF#06-0675的(111)、(220)和(311)衍射峰均吻合很好。初始粉末的SEM图像表明,粒径大小分布均匀,经过晶粒统计软件计算之后可以确认平均晶粒大小和标识的一样,即为0.5和10μm。

    图  1  初始粉末的XRD图
    Figure  1.  XRD pattern of starting material
    图  2  初始粉末的SEM图
    ((a)0.5μm金刚石粉, (b)10μm金刚石粉)
    Figure  2.  SEM image of starting material
    ((a)0.5μm diamond powder, (b)10μm diamond powder)

    本实验采用国产铰链式6×1400t大腔体静高压装置,六面顶压机形成一个封闭的正方体大压腔,每个顶锤的锤面为正方形,边长为23.5mm。选用叶蜡石(一种含水的铝硅酸盐:Al2(Si4O10)(OH)2)作为传压介质,白云石作为内层保温材料,碳管碳片作为加热装置,盐管盐片形成一个类静水压腔包裹着样品。图 3为实验装置示意图。

    图  3  实验分层组装
    Figure  3.  Bi-layered assembly of experiment

    图 4为实验温压变化曲线,其中压力p为六面顶压机油压。在进行实验之前,选取Bi、Tl、Ba、ZnTe材料,根据其相变点对样品腔进行压力标定实验[21],测出样品腔内实际压力与压机油压的关系并拟合成曲线。同理,利用B型热电偶标出加热功率(P)和实际温度之间的关系并拟合成曲线。样品腔内取粒度为0.5μm的金刚石粉0.3g,粒度为10μm的金刚石粉0.2g,WC-16%Co基底厚度为3mm。实验中选取油压为45MPa,加热功率为2.1kW,对应的实际压力和实际温度分别为5.5GPa和1500℃。经过升压-保压-加热-保温-降温-卸压等过程之后,从压机中取出合成好的实验样品。

    图  4  实验温压条件曲线
    Figure  4.  Temperature and pressure condition

    使用金刚石磨盘对合成的PCD样品进行表面抛光处理,然后对其进行硬度测试。实验所用硬度计为日产FV-700型维氏硬度计,加载压力49.0N,压头压入样品之后保压时间为15s。每个样品均选取边缘4个点和中心1个点进行测试。

    实验合成的PCD样品直径约为8mm,高约为5mm(其中亚微米层高约为0.8mm)。对合成的亚微米级PCD层进行系统的分析测试,根据检测结果表征合成PCD样品的性能,并由此反向探索合成PCD样品的最佳工艺条件。

    硬度是反映超硬材料性能的重要指标。对不同温压条件下合成的PCD样品进行维氏硬度测试后发现,当压力为5.5GPa、温度为1500℃,保温15min时,合成样品的维氏硬度测试结果表现优异。在此条件下重复进行多次实验,结果显示,加载力为49.0N时,合成样品的硬度稳定在(57.0±2.0)GPa,说明此条件为稳定、可重复的合成工艺条件。

    表 1为不同压力、温度和保温时间条件下测得所得合成PCD样品的维氏硬度。从表 1可以看出,在实验温压条件确定的情况下,保温时间对所合成样品硬度的改善有着重要作用。

    Hv=1854.4F/L2
    (1)
    表  1  合成PCD样品的维氏硬度
    Table  1.  Vickers hardness of PCD samples
    T/℃ t/min Hv/GPa
    p=5.0GPa p=5.5GPa
    1300 5 15-22 15-22
    15 15-22 15-22
    1400 5 53.03 32.14
    15 47.08 50.33
    1450 5 19.22 18.36
    15 30.84 52.56
    1500 5 21.22 42.97
    15 39.95 57.03
    1550 5 23.40 39.14
    15 48.40 43.97
    下载: 导出CSV 
    | 显示表格

    式中:Hv为维氏硬度,F为加载力,L为压痕对角线长度。硬度测试过程中存在读数误差,通过(1)式并结合误差分析理论可以计算硬度测试结果存在的误差。

    表 1中压力为5.5GPa、保温时间为15min、不同温度条件下测得的维氏硬度绘制于图 5中。从图 5可以看出:在一定的压力条件下,温度对于PCD的合成也起到了重要的作用; 当温度为1500℃时,合成PCD样品的硬度达到最大值。

    图  5  不同温度下合成PCD样品的维氏硬度
    Figure  5.  Vickers hardness of sintered PCD samples at different temperatures

    对实验合成的PCD样品进行镶嵌抛光处理,然后对样品进行XRD分析。图 6为合成PCD样品的XRD图谱。由图 6可知,合成PCD样品中存在立方相金刚石、碳化钨、钴3种成分。该样品合成压力为5.5GPa,温度为1500℃,保温时间为15min。图 7为实验合成PCD样品的断面图,其中可以清晰地看到3层不同的颜色区域。由图 3可知:区域1为0.5μm金刚石微粉,其厚度约为0.8mm; 区域2为10μm金刚石微粉,厚度约为1.4mm; 区域3为WC-16%Co基底,厚度约为2.9mm。

    图  6  5.5GPa、1500℃下合成PCD样品的XRD图谱
    Figure  6.  XRD spectrum of sintered PCD at 5.5GPa and 1500℃
    图  7  合成PCD样品断面图
    Figure  7.  Fracture surface of the sintered PCD sample

    D─D键键合情况对于判定烧结PCD样品质量的好坏有着重要作用[22]。为了观察合成PCD样品中D─D键的键合情况,首先对实验合成的PCD样品进行24h表面酸处理,然后对酸处理后的样品进行SEM测试,如图 8所示。对图 8(a)进行粒径统计可得,合成的PCD样品晶粒尺寸平均为326nm。此外,图 8中可以看到许多的细小颗粒,说明在高压的条件下部分金刚石颗粒发生了破碎。进一步分析图 8可知,合成PCD样品中存在大量的D─D键。对样品进行电阻测量显示,在5.5GPa、1500℃、保温时间为15min条件下,合成PCD样品的电阻可达1.553Ω,由此也可以判定合成PCD样品中存在着大量D─D键键合的情况。

    图  8  5.5GPa、1500℃下合成亚微米PCD样品经过酸处理后的SEM图
    Figure  8.  SEM image of sub-micron PCD sample sintered at 5.5GPa and 1500℃ after acid treatment

    Co作为黏结剂存在于样品中,其分布情况对于烧结样品的好坏具有非常重要的参考价值。对未经酸处理的PCD样品进行EBSD分析,以观察Co在样品中的分布情况。图 9为合成PCD的背散射图像,图像中白色部分为Co,可以看出,Co在样品中分布比较均匀,以孔形存在于金刚石颗粒之间,没有出现大量聚集在一起的情况。对图片放大后(见图 9(b))发现,Co存在于金刚石颗粒之间,并充满金刚石颗粒之间形成的圆形孔隙。

    图  9  5.5GPa、1500℃下合成亚微米PCD样品的EBSD图
    Figure  9.  EBSD image of sub-micron PCD sample sintered at 5.5GPa and 1500℃

    为了获得Co在合成PCD样品中的含量,对合成的PCD样品进行EDS能谱分析。图 10为合成PCD样品的表面能谱图,其中表面白色条状物为抛光未彻底钼杯剩余的残屑。表 2为通过EDS分析得到的PCD样品中表层各元素的含量,其中wχ分别为质量分数和物质的量分数,由表 2可知,Co的质量分数为13.84%。

    图  10  5.5GPa、1500℃下合成亚微米PCD样品EDS图
    Figure  10.  EDS image of sub-micron PCD sample sintered at 5.5GPa and 1500℃
    表  2  5.5GPa、1500℃合成的亚微米PCD样品中各元素的含量
    Table  2.  Element content of sub-micron PCD sample at 5.5GPa and 1500℃
    Element w/% χ/%
    C (K) 77.70 92.80
    O (K) 3.69 3.31
    Co (K) 13.84 3.37
    Mo (L) 2.08 0.31
    W (M) 2.69 0.21
    下载: 导出CSV 
    | 显示表格

    以0.5、10.0μm金刚石微粉为初始粉料,以WC-16%Co为基底材料,利用国产铰链式6×1400t大腔体静高压装置,以分层组装的方法对PCD的合成工艺进行了探索性实验,并对合成的样品进行了一系列的分析和检测,主要结论如下。

    (1) 性能优异的亚微米级PCD稳定、可重复的合成工艺为:压力5.5GPa,温度1500℃,保温时间15min。实验合成的亚微米级PCD样品硬度达到57GPa,远超工业产品性能。

    (2) 使用分层组装的方法可以使Co均匀地分散在PCD样品中,并充满金刚石颗粒之间形成的圆形孔隙,从而保证样品具备均匀、优异的性能。

    (3) 温度和保温时间对合成PCD样品性能具有很大的影响,通过选取合理的温度和保温时间可以合成出性能较好的PCD样品。

  • 图  新型结构炮口制退器

    Figure  1.  New-structure muzzle brake

    图  冲击式炮口制退器

    Figure  2.  Impulsive type muzzle brake

    图  轴向截面膛口压力云图

    Figure  3.  Pressure nephogram at axial plane

    图  径向截面膛口压力云图

    Figure  4.  Pressure nephogramat at radial cross section

    图  膛口流场速度矢量图

    Figure  5.  Velocity vector diagram of muzzle flow field

    图  不同时刻炮轴中心线压力分布曲线

    Figure  6.  Pressure distributions along central axis of gun at different moments

    图  不同时刻侧孔径向轴线压力分布曲线

    Figure  7.  Pressure distributions along radial axis of side holes at different moments

    图  膛口流场超压监测点阵

    Figure  8.  Overpressure monitoring points of muzzle flow filed

    图  膛口超压峰值分布对比

    Figure  9.  Comparison of overpressure peak distribution

    图  10  制退器腔内压力监测点

    Figure  10.  Monitoring points inside muzzle brake

    图  11  腔内各点压力随时间变化

    Figure  11.  Variation of overpressure inside muzzle brake with time

    图  12  后效期初期制退器受力曲线

    Figure  12.  Curve of muzzle brake force in the early gas ejection period

    图  13  膛底、制退器及总受力曲线

    Figure  13.  Curves of chamber bottom force, muzzle brake force and total thrust

    图  14  制退器结构的等效应力随时间变化曲线

    Figure  14.  Variation of equivalent stress with time of muzzle brake

    图  15  制退器结构的变形随时间变化曲线

    Figure  15.  Variation of deformation with time of muzzle brake

    图  16  0.25 ms时制退器结构的等效应力云图

    Figure  16.  Equivalent stress distribution of the muzzle brake at 0.25 ms

    图  17  1.10 ms时制退器结构的变形云图

    Figure  17.  Deformation distribution of the muzzle brake at 1.10 ms

    表  1  新型炮口制退器结构参数

    Table  1.   Structural parameters of new-structure muzzle brake

    StructureBackward angle/(°)Area/mm2
    Inner side hole 15392.96
    Inner side hole 210470.34
    Inner side hole 310392.96
    Outer side hole078.50
    Front hole146.54
    下载: 导出CSV

    表  2  3D打印钛合金材料参数

    Table  2.   Material parameters of titanium alloy

    MaterialDensity/(g·cm−3)Young’s modulus/GPaYield strength/MPaTensile strength/MPaPoisson’s ratio
    TC44.4311894410580.3
    下载: 导出CSV

    表  3  不同制退器的超压峰值监测结果

    Table  3.   Overpressure peak values of different muzzle brakes

    Monitor pointOverpressure peak/kPaΔnew/kPaΔimpact/kPa
    Smooth muzzleNew-structure muzzle brakeImpact muzzle brake
    P01117.84362.54200.83244.7082.99
    P0246.3165.4653.2019.156.89
    P1183.8079.0344.58−4.76−39.21
    P1228.0723.6119.56−4.46−8.51
    P2148.8429.2830.90−19.56−17.93
    P2219.7611.4515.40−8.31−4.36
    P3122.2918.8550.87−3.4528.57
    P3210.546.9915.30−3.554.76
    P4128.2726.4541.75−1.8213.48
    P4211.559.549.73−2.01−1.82
    P5134.7544.6840.099.935.34
    P5215.7114.7913.58−0.91−2.13
    下载: 导出CSV

    表  4  制退器效率对比

    Table  4.   Comparison of efficiency for different muzzle brakes

    StructureM/kgI/(N·s)v/(m·s−1)vmax/(m·s−1)E/Jη/%
    Smooth muzzle7.09−123.68−47.00−64.4514719.0
    New-structure muzzle brake8.37−29.40−40.23−43.748003.245.62
    Impact muzzle brake8.10−26.15−41.47−44.708091.345.03
    下载: 导出CSV
  • [1] 张相炎, 郑建国, 袁人枢. 火炮设计理论[M]. 北京: 北京理工大学出版社, 2014: 153–157.

    ZHANG X Y, ZHENG J G, YUAN R S. Theory of artillery gun design [M]. Beijing: Beijing Institute of Technology Press, 2014: 153–157.
    [2] 刘凯, 赵俊利, 郭利强, 等. 钛合金在炮口制退器上的应用 [J]. 兵工自动化, 2016, 35(6): 94–96. doi: 10.7690/bgzdh.2016.06.022

    LIU K, ZHAO J L, GUO L Q, et al. Application of titanium alloy in muzzle brake [J]. Ordnance Industry Automation, 2016, 35(6): 94–96. doi: 10.7690/bgzdh.2016.06.022
    [3] 吴喜富, 郑建国. 基于流固耦合的复合结构炮口制退器强度分析 [J]. 兵工自动化, 2016, 35(7): 19–22. doi: 10.7690/bgzdh.2016.07.006

    WU X F, ZHENG J G. Strength analysis of muzzle brake with composite structure by method of fluid solid coupling [J]. Ordnance Industry Automation, 2016, 35(7): 19–22. doi: 10.7690/bgzdh.2016.07.006
    [4] 岳明凯, 刘欣宁. 炮口制退器现状及其发展趋势 [J]. 兵工自动化, 2015, 34(3): 1–6. doi: 10.7690/bgzdh.2015.03.001

    YUE M K, LIU X N. Situation and development of muzzle brake [J]. Ordnance Industry Automation, 2015, 34(3): 1–6. doi: 10.7690/bgzdh.2015.03.001
    [5] 葛苗冉, 袁军堂, 汪振华, 等. 基于正交试验的炮口制退器结构设计与性能分析 [J]. 兵器装备工程学报, 2019, 40(12): 160–164. doi: 10.11809/bqzbgcxb2019.12.032

    GE M R, YUAN J T, WANG Z H, et al. Structural design and performance analysis of muzzle brake based on orthogonal test [J]. Journal of Ordnance Equipment Engineering, 2019, 40(12): 160–164. doi: 10.11809/bqzbgcxb2019.12.032
    [6] 谭添, 戴劲松, 王茂森, 等. 封闭反射膨胀装置流场仿真分析 [J]. 火炮发射与控制学报, 2019, 40(4): 1–5. doi: 10.19323/j.issn.1673-6524.2019.04.001

    TAN T, DAI J S, WANG M S, et al. Flow field simulation of closed reflective inflation device [J]. Journal of Gun Launch & Control, 2019, 40(4): 1–5. doi: 10.19323/j.issn.1673-6524.2019.04.001
    [7] 李萌蘖, 李闯, 李绍宏. TC4合金增材制造的研究现状 [J]. 昆明理工大学学报(自然科学版), 2018, 43(6): 20–27.

    LI M N, LI C, LI S H. Current status of TC4 alloy additive manufacturing [J]. Journal of Kunming University of Science and Technology (Natural Science), 2018, 43(6): 20–27.
    [8] 张焕好, 陈志华, 姜孝海, 等. 高速弹丸穿越不同制退器时的膛口流场波系结构研究 [J]. 兵工学报, 2012, 33(5): 623–629.

    ZHANG H H, CHEN Z H, JIANG X H, et al. Investigation on the blast wave structures of a high-speed projectile flying through different muzzle brakes [J]. Acta Armamentarii, 2012, 33(5): 623–629.
    [9] ZHANG H H, CHEN Z H, JIANG X H, et al. Investigations on the exterior flow field and the efficiency of the muzzle brake [J]. Journal of Mechanical Science and Technology, 2013, 27(1): 95–101. doi: 10.1007/s12206-012-1223-8
    [10] 代淑兰, 许厚谦, 肖忠良. 带制退器的膛口燃烧流场并行数值模拟 [J]. 弹道学报, 2009, 21(4): 84–87.

    DAI S L, XU H Q, XIAO Z L. Numerical simulation of muzzle combustion flow field with brake by parallel computation [J]. Journal of Ballistics, 2009, 21(4): 84–87.
    [11] LEI H X, WANG Z J, ZHAO J L. Stress analysis of muzzle brake by using fluid-solid coupled method [J]. Journal of Engineering Science and Technology Review, 2016, 9(4): 48–55. doi: 10.25103/jestr.094.07
    [12] CHATURVEDI E, DWIVEDI R K. Computer aided design and analysis of a tunable muzzle brake [J]. Defence Technology, 2019, 15(1): 89–94. doi: 10.1016/j.dt.2018.06.011
    [13] 吴彦霖. 基于SLM制备的钛合金三维点阵结构的力学性能研究 [D]. 重庆: 重庆大学, 2016.

    WU Y L. An investigation into the mechanical properties of Ti6Al4V lattice structures manufactured using selective laser melting [D]. Chongqing: Chongqing University, 2016.
  • 加载中
图(17) / 表(4)
计量
  • 文章访问数:  6890
  • HTML全文浏览量:  3461
  • PDF下载量:  38
出版历程
  • 收稿日期:  2020-06-12
  • 修回日期:  2020-07-03

目录

/

返回文章
返回