Muzzle Blast Wave Investigation and Performance Analysis of New-Structure Muzzle Brake Based on Numerical Simulation
-
摘要: 面向炮口制退器综合性能设计需求,基于增材制造技术优势,提出了一种叠加冲击式内壁与反作用式外孔特征的新型小口径轻金属炮口制退器方案。基于三维无黏Euler方程,建立了后效期火药气体排空过程的有限元仿真模型,通过膛口流场仿真与流固耦合分析,研究了膛口流场发展、膛口冲击波与超压分布特征、制退器效率与其强度性能。结果表明:相比于传统结构,采用同种轻质钛合金材料的新型制退器具有较高的制退效率,对于机载平台具有较优的膛口超压分布特征和低冲击波危害效应,并且结构强度满足使用要求。Abstract: Taking advantages of the additive manufacturing technology, our work proposes a new small-caliber light metal gun muzzle brake scheme with an impact-type inner wall and reaction-type outer hole feature. On the basis of the finite element simulation model of the after-effective period gunpowder gas evacuation, which was established by means of the three-dimensional non-viscous Euler equation, the muzzle flow field was simulated, the fluid-solid interaction was analyzed, and then the development of muzzle flow field, the characteristics of the muzzle blast wave and overpressure distribution, as well as the efficiency of the muzzle brake and its strength performance were investigated. The results showed that, compared with the traditional structure, the new muzzle brake with the same kind of lightweight titanium alloy material has a higher recoil efficiency, being characteristic of better muzzle overpressure distribution and lower shock wave hazard benefit for the airborne platform. Additionally, its structural strength can meet the operating requirements.
-
Key words:
- muzzle brake /
- additive manufacturing /
- muzzle blast wave /
- overpressure
-
火炮发射时,后效期高温高压的火药气体从膛口瞬时流出,对火炮主体形成很强的后坐力,对炮架产生强冲击作用,影响火炮射击精度和机动性能,制约高性能火炮的发展。炮口制退器作为一种安装在炮口部位的排气装置,通过控制后效期火药气体的流量分配和气流速度,减小射击时火药气体作用于后坐部分的冲量,并为炮身提供制退力,减小火炮后坐动能和炮架的射击载荷,因而成为一种广泛应用的反后坐技术[1]。然而,炮口制退器的使用具有一定的负面效应。一方面,炮口制退器造成炮口区域强激波,加剧后效期火药气体流场的复杂性及膛口焰现象,加大火炮发射初速扰动及膛口冲击波超压,对火炮身管周围设备及操作人员造成不良影响;另一方面,加装制退器后火炮身管本身质量及挠度的增大也对射击精度造成不利影响。炮口制退器自身的质量和结构决定了其制退性能及负面危害表现,通过选材及结构的优化设计实现高性能、低危害的综合性能成为当前炮口制退器研究与发展的方向[2–6]。
按照结构形式的不同,传统的炮口制退器分为冲击式、反作用式及冲击反作用式。由于工艺、重量的限制,高效率与低负面效应往往难以兼顾。近年来,增材制造技术的发展为复杂异型结构的加工制造以及钛合金等难加工轻金属材料的成型制备提供了有力支撑[7],也为传统装置结构的优化与创新设计提供了更大的自由空间。随着计算机技术和计算流体力学的发展,数值模拟已成为低成本研究膛口流场及膛口装置的重要手段,例如:张焕好等[8–9]、代淑兰等[10]利用数值模拟方法,对不同类型炮口制退器的膛口流场波系、膛口激波及二次焰特征进行了详细研究,得到了制退器装置效率;Lei等[11]、Chaturvedi等[12]通过将流场仿真与流固耦合相结合,实现了炮口制退器结构应力与变形响应研究及性能评价。本研究利用增材制造技术优势,针对机载火炮小口径炮口制退器综合性能的设计需求,提出一种叠加冲击式与反作用式传统制退器结构特征与优点的钛合金新型炮口制退器结构方案,对膛口流场发展过程及特征进行数值模拟分析,并以传统的高效率冲击式炮口制退器为参照,对新型结构制退器的综合性能予以评价。
1. 制退器结构
炮口制退器的基本构型是后端与炮管封闭相接,前端收缩,中间拥有膨胀腔室、内部挡板和侧孔的圆筒形结构。冲击式炮口制退器有较大的腔室直径和侧孔面积,火药气体膨胀较为充分,大部分经侧向流出并冲击反射挡板,制退效率高,但炮口冲击波危害也较大。反作用式炮口制退器的直径一般较小,拥有多排小面积侧孔,没有或有很小的反射挡板,制退效率较低,但膛口火药气体流场相对简单,对弹丸初速的扰动小,不会影响大威力脱壳穿甲弹发射。基于上述两种制退器的结构特征和性能的经验规律,本研究以低冲击波危害和高制退效率为导向,在三维建模平台上设计一种面向机载航炮的30 mm小口径新型结构炮口制退器,如图1所示。
新型制退器内层为冲击式内壁结构,见图1,其上设置编号为1~3的3组阵列式大面积后倾侧孔,旨在引导火药气体出膛后的侧向分流,以实现类似于冲击式炮口制退器的高制退效率。此外,3组内层侧孔的后倾角度及面积采用非均等设计,参照一般冲击式炮口制退器的结构参数,对三维设计模型直接设定或导出得到内层侧孔的非均等结构参数,见表1。其中:为了限制火药气体在临近膛口位置的向后偏转角度及炮口后方冲击波的强度,第1组侧孔的后倾角度及面积均小于其他两组;根据火药气体出膛以后沿轴向流量逐渐衰减的特点,第2组侧孔的面积设置为大于其他两组。新型结构制退器外层为胡椒瓶式套筒结构,其上设置成组分布的阵列式小面积圆孔,各组圆孔与内层侧孔周向重叠,旨在借鉴反作用式炮口制退器的特点,对火药气体进行分散和疏导,从而抑制侧向强激波的产生,实现较低的冲击波危害效益。新型制退器弹孔处为收缩式前端面,并设有分流孔,以在出口处进一步减小中央弹孔流量及火药气体在制退器出口部位的膨胀扰动。内层结构与外层结构通过布置在孔间的阵列式筋板相连,实现整体结构的一体化,便于增材制造成型。同时,结构中不存在过于复杂或与轴向夹角过大的成型曲面,满足增材制造高质量成型工艺要求。设计方案选材为3D打印钛合金粉末,材料参数[13]见表2。为评价新型炮口制退器的综合性能,以图2所示的传统冲击式炮口制退器为对比分析对象。
表 1 新型炮口制退器结构参数Table 1. Structural parameters of new-structure muzzle brakeStructure Backward angle/(°) Area/mm2 Inner side hole 1 5 392.96 Inner side hole 2 10 470.34 Inner side hole 3 10 392.96 Outer side hole 0 78.50 Front hole 146.54 表 2 3D打印钛合金材料参数Table 2. Material parameters of titanium alloyMaterial Density/(g·cm−3) Young’s modulus/GPa Yield strength/MPa Tensile strength/MPa Poisson’s ratio TC4 4.43 118 944 1058 0.3 2. 数值方法与模型参数
采用无黏三维Euler方程描述带制退器的膛口流场气体流动。考虑到计算目的及效率,对炮口气流进行适当简化,忽略气固两相性、多组分、化学反应以及运动弹丸的影响,将火药气体与外界大气视为理想气体,且满足气体状态方程。有限元仿真的封闭方程组为
∂U∂t+∂F∂x+∂G∂y+∂H∂z=0 (1) e=pγ−1+12ρ(u2+v2+w2) (2) pV=nRT (3) 式中:e为单位质量气体的总能量,
p 为火药气体压强,V为气体体积,T为体系温度,U=[ρ,ρu,ρv,ρw,e]T ,F=[ρu,ρu2+p,ρuv,ρuw,(e+p)u]T ,G=[ρv,ρuv,ρv2+p,ρvw,(e+p)v]T ,H=[ρw,ρuw,ρvw,ρw2+p,(e+p)w]T ,ρ 为火药气体密度,u、v、w分别为x、y、z方向的速度分量,γ 为理想气体绝热指数,R为通用气体常数,n为气体的物质的量。方程的离散采用有限体积法,时间推进采用二阶精度Runge-Kutta法,对流项采用二阶精度Roe格式,湍流模型采用单方程Spalart-Allmaras模型,数值计算在ANSYS FLUENT 19.0中进行。后效期火药气体排空过程数值模拟的计算域设置为以炮口端面中心为原点、直径为3 m的球形区域,以使火药气体充分发展。计算初始时刻为弹丸飞离炮口时刻,以膛底位置为x=0,此时膛内气体沿炮轴方向任意一点的压力px、速度vx和温度Tx的初始分布状态表示为
px=pd[1+mw2ψ1mq(1−x2L2)] (4) vx=vdLx (5) Tx=pxMgρgR (6) ρg=mwVc+LS (7) 式中:
ρ g为膛内火药气体的平均密度,mw、mq、Mg分别为装药量、弹丸质量、火药气体摩尔质量,Vc、L、S分别为药室容积、身管行程长度、身管横截面积,ψ1 为内弹道次要功系数,vd为弹底火药气体速度,pd为膛口处弹底压力。本算例中,未加装炮口制退器时火炮身管质量M0=7.09 kg,初始时刻弹底火药气体速度vd=780 m/s,pd=27.9 MPa,膛口外为常温常压条件。计算时间步长取0.01 ms,当膛内气压降低至约0.2 MPa时视为后效期结束,终止计算。3. 结果分析
3.1 膛口冲击波
图3、图4所示为新型制退器轴向和径向截面处膛口冲击波在不同时刻的压力云图,为突出显示流场发展过程,图中隐去流场超高压及部分负压区域。图5为火药气体流场充分发展时轴向和径向截面处速度矢量图。可以看出,后效期内火药气体在流入和穿出新型制退器的过程中,与内外双层结构依次产生交互作用,形成独特的流场结构与激波形貌,并构成制退效益的基础。后效期膛口冲击波的形成过程可以分为3个阶段。(1)火药气体膨胀加速进入制退器腔室后,依次流经内层侧孔和外层小孔喷出,形成侧孔冲击波,见图3(a)和图4(a)。(2)随着火药气体的向前推进和持续补充,沿炮轴方向各外层圆孔依次形成各自的侧孔冲击波,并且伴有追赶、相交和叠加过程,最终在炮口一侧合并成为一个大的冲击波。由于膨胀速率不同,在此过程中出现了后排侧孔冲击波追赶并超过前排侧孔冲击波的现象,见图3(b)和图3(c);与此同时,周向上不同角度的三维侧孔冲击波也逐渐相交合并,见图4(b)和图4(c)。(3)火药气体流出弹孔以后,在出口处形成迅速膨胀的弹孔冲击波,并与侧孔冲击波相交,形成具有各向叠加特征的典型膛口波系结构,见图3(d)、图4(d)和图5。
3.2 超压分布
图6和图7为炮轴中心线、侧孔位置径向轴线处不同时刻压力分布曲线。可以看出,火药气体到达弹孔之前,不断提升新型制退器腔室内的压力水平。t=0.65 ms时,火药气体从中央弹孔流出并迅速提升当地压力,直至t=1.00 ms时弹孔处压力峰值超过制退器腔室内压力最大值,见图6(a)和图6(b)。与此同时,径向轴线方向的压力分布呈现随时间剧烈波动的高瞬态特征,见图7(a)。t=1.00 ms之后,火药气体流场逐渐衰减,腔室内外压力水平均逐渐下降(见图6(c)和图7(b)),同时径向轴线方向压力水平也呈现出规则的逐级下降趋势。
根据机载小口径航炮的实际工况,以炮口端面为中心、炮膛轴线为正方向,在0.5和1.0 m 的周向距离上,炮口区域侧向0°~150°范围内设置若干膛口流场超压监测点,如图8所示,以精确分析加装制退器后炮口区域的超压特征,评估炮口冲击波对机载平台的危害效应。得到的后效期内各点超压峰值监测结果与差值情况如表3所示。其中,超压值代表压力监测值高于环境压力的部分,加装新型制退器与传统冲击式炮口制退器情况下的超压峰值与光膛口情况下的超压峰值差值分别记为Δnew、Δimpact,负值代表产生减小超压效果,正值代表超压增大。同时,根据超压峰值,绘制了不同情况下炮口超压分布,如图9所示。
表 3 不同制退器的超压峰值监测结果Table 3. Overpressure peak values of different muzzle brakesMonitor point Overpressure peak/kPa Δnew/kPa Δimpact/kPa Smooth muzzle New-structure muzzle brake Impact muzzle brake P01 117.84 362.54 200.83 244.70 82.99 P02 46.31 65.46 53.20 19.15 6.89 P11 83.80 79.03 44.58 −4.76 −39.21 P12 28.07 23.61 19.56 −4.46 −8.51 P21 48.84 29.28 30.90 −19.56 −17.93 P22 19.76 11.45 15.40 −8.31 −4.36 P31 22.29 18.85 50.87 −3.45 28.57 P32 10.54 6.99 15.30 −3.55 4.76 P41 28.27 26.45 41.75 −1.82 13.48 P42 11.55 9.54 9.73 −2.01 −1.82 P51 34.75 44.68 40.09 9.93 5.34 P52 15.71 14.79 13.58 −0.91 −2.13 由表3和图9可知,相比于光膛口情况,两种类型制退器的膛口冲击波超压特点可总结如下:(1)炮口正前0°方向,0.5和1.0 m两个距离上两种制退器均造成超压峰值的明显增大,并且新型制退器的增压幅度更大;(2)自炮口侧向30°方向向后,新型制退器膛口冲击波超压峰值呈U形分布,冲击型制退器呈W形分布;(3)30°~120°侧向范围内,新型制退器均产生减小超压效果,各点超压峰值降低幅度为1.82~19.56 kPa,正侧方90°处超压峰值最小,该方位0.5和1.0 m距离处分别为18.85和6.99 kPa;(4)30°~120°侧向范围内,冲击型制退器对应的超压峰值高低不等,在正侧方90°方向0.5和1.0 m距离上均形成增大超压效果,超压峰值分别为50.87和15.30 kPa;(5)150°方向0.5 m处,新型制退器对应的超压峰值略高于冲击型制退器及光膛口情况,1.0 m处同样略高于冲击型制退器但低于光膛口情况。
综上可知,新型制退器的膛口冲击波表现出高超压集中于炮口正前方、低超压集中于正侧方(炮口与机载平台直线距离最近位置处)的分布特征,同时其后侧方的整体超压水平也较好地控制在低于或近似光膛口情况,相比传统冲击型炮口制退器,对于缓解因制退器安装导致的机身及炮口装置的冲击侵害更具优势。同时应该注意到,为了追求高制退效率,新型制退器设置了较大的侧孔后倾角度及向后偏转面积,因此在炮口后侧150°方向0.5 m距离处超压出现一定的增大,说明新型制退器的结构参数仍存在改进空间,以实现高制退效率和更低冲击波危害的进一步平衡。
3.3 制退效率及结构强度
为分析后效期内新型制退器受载变化及制退力产生机制,沿炮膛轴线在新型炮口制退器腔室内层侧孔1~3位置及前端面侧孔处分别设置压力监测点,如图10所示。图11为后效期初期阶段各监测点压力随时间的变化曲线。图12和图13分别为后效期初期制退器受力、后效期全程火炮身管各部分受力随时间的变化曲线。
结合图11和图12可知,火药气体的流出和偏转作用不断产生对火炮身管的后坐力和向前的制退力,后效期初期制退器受力曲线呈现阶梯式上升趋势。图12中点Ⅰ~点Ⅳ处的抬升显示了当火药气体依次流经制退器内层侧孔时逐级产生的冲撞和制退效果。火药气体流出中央弹孔之前,对制退器前端面的冲击作用使制退器受力曲线继续抬升,直至达到峰值点Ⅴ,此后由于火药气体对制退器外壁面产生反向冲击作用,又使制退器的正向受力水平略降至峰值以下。
依据仿真监测结果,首先由总受力F随时间的变化曲线计算得到后效期火药气体对身管的总冲量I,再结合动量守恒定律及冲量定律,计算得到后效期起始时刻身管自由后坐速度v、后效期结束时刻身管最大后坐速度vmax以及火炮后坐动能E,进而依据定义得到制退器的效率
η ,具体计算公式为v=mq+0.5mwM+mw+mqvd (8) I=∫τ0Fdt=M(vmax−v) (9) E=Mv2max2 (10) η=E0−EE0×100% (11) 式中:
τ 为后效期持续时间,E0为光膛口时火炮后坐动能,M为加装不同制退器后身管总质量。本算例中,设光膛口时火炮身管质量为7.09 kg,各种情况下的计算结果如表4所示,其中负号代表火炮发射相反方向。可见,新型制退器的制退效率达到45.62%,满足小口径炮口制退器的使用要求,并略超过以高效率为特点的传统冲击式炮口制退器性能水平。不过,新型炮口制退器情况下火药气体对膛口装置的反作用冲量相比传统冲击式炮口制退器的提升不明显,说明新型炮口制退器在内部结构参数上还有优化空间,可以进一步提升新结构形式的制退性能。表 4 制退器效率对比Table 4. Comparison of efficiency for different muzzle brakesStructure M/kg I/(N·s) v/(m·s−1) vmax/(m·s−1) E/J η/% Smooth muzzle 7.09 −123.68 −47.00 −64.45 14719.0 New-structure muzzle brake 8.37 −29.40 −40.23 −43.74 8003.2 45.62 Impact muzzle brake 8.10 −26.15 −41.47 −44.70 8091.3 45.03 制退器与后效期火药气体相互作用过程中,其本身所承受的动态载荷作用通过流固耦合技术由流体计算结果传递至制退器各壁面,并通过瞬态动力学有限元计算方法求得变形与等效应力随时间的变化曲线,如图14和图15所示。后效期初期,高温高压的火药燃气出膛后迅速膨胀,在剧烈变化的流体压力动态载荷作用下,制退器所承受的等效应力及其变形也在短时间内迅速抬升达到峰值;约2 ms以后,随着流场的充分发展并衰减,制退器结构的应力及变形转为逐渐下降趋势。
图16和图17分别为新型制退器的等效应力和变形在其峰值时刻的分布云图。t=0.25 ms时,火药气体首先对炮口端面处的侧孔产生冲击,并在该时刻造成距离炮口端面第2排外层圆孔处出现全局最大范式等效应力,为224.37 MPa。t=1.10 ms时,火药气体已从中央弹孔流出,火药气体对制退器前端面的冲击作用使得最大变形出现在前端面与中央弹孔的结合处,约为14 μm,其他部位的变形呈现从炮口端面沿炮口正向递减趋势。由以上结果可知,后效期内制退器结构的最大应力小于钛合金材料许用应力,满足结构强度的使用要求。
4. 结 论
提出了一种叠加冲击式内壁与反作用式外孔特征的新型小口径钛合金炮口制退器结构方案,基于膛口流场数值仿真、流固耦合技术等对其性能进行分析,并与传统冲击型炮口制退器进行对比。结果表明:新型结构方案满足小口径炮口制退器性能要求,制退效率相比传统冲击型制退器略有提高;在机载航炮应用背景下,其膛口冲击波超压分布特征及峰值水平相比传统冲击型制退器具有更好的机身与炮口装置防护优势和低冲击波危害效益,同时结构强度满足使用要求。新型结构方案对于新材料和新制造技术背景下炮口制退器的创新设计和发展具有参考意义。本研究所展示的新型制退器结构参数在高效率和更低负面危害的平衡上具有进一步改进与优化空间。研究过程中膛口流场数值模拟存在一定简化,仿真结果与实际情况存在一定偏差。另外,由于后效期膛口火药气体具有高温属性,一般借助防烧蚀涂层等技术强化钛合金装置的高温使用性能,因此未在结构强度分析中考虑热耦合效应。
-
表 1 新型炮口制退器结构参数
Table 1. Structural parameters of new-structure muzzle brake
Structure Backward angle/(°) Area/mm2 Inner side hole 1 5 392.96 Inner side hole 2 10 470.34 Inner side hole 3 10 392.96 Outer side hole 0 78.50 Front hole 146.54 表 2 3D打印钛合金材料参数
Table 2. Material parameters of titanium alloy
Material Density/(g·cm−3) Young’s modulus/GPa Yield strength/MPa Tensile strength/MPa Poisson’s ratio TC4 4.43 118 944 1058 0.3 表 3 不同制退器的超压峰值监测结果
Table 3. Overpressure peak values of different muzzle brakes
Monitor point Overpressure peak/kPa Δnew/kPa Δimpact/kPa Smooth muzzle New-structure muzzle brake Impact muzzle brake P01 117.84 362.54 200.83 244.70 82.99 P02 46.31 65.46 53.20 19.15 6.89 P11 83.80 79.03 44.58 −4.76 −39.21 P12 28.07 23.61 19.56 −4.46 −8.51 P21 48.84 29.28 30.90 −19.56 −17.93 P22 19.76 11.45 15.40 −8.31 −4.36 P31 22.29 18.85 50.87 −3.45 28.57 P32 10.54 6.99 15.30 −3.55 4.76 P41 28.27 26.45 41.75 −1.82 13.48 P42 11.55 9.54 9.73 −2.01 −1.82 P51 34.75 44.68 40.09 9.93 5.34 P52 15.71 14.79 13.58 −0.91 −2.13 表 4 制退器效率对比
Table 4. Comparison of efficiency for different muzzle brakes
Structure M/kg I/(N·s) v/(m·s−1) vmax/(m·s−1) E/J η/% Smooth muzzle 7.09 −123.68 −47.00 −64.45 14719.0 New-structure muzzle brake 8.37 −29.40 −40.23 −43.74 8003.2 45.62 Impact muzzle brake 8.10 −26.15 −41.47 −44.70 8091.3 45.03 -
[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.022LIU 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.006WU 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.001YUE 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.032GE 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.001TAN 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. -