Gas dynamics of gunpowder for clearance seal of ejection device in projectile base
-
摘要: 针对弹尾开孔安装用于侵彻数据回收的弹射装置后,存在的弹射间隙会使弹体内的测试装置受到火药气体侵蚀破坏的问题,设计了间隙密封结构,对流经间隙的火药气体进行了气体动力学的建模分析,并对密封结构进行了膛压和密封空腔压强测试。结果表明:火药燃烧产生的气体为可压缩性气体,在药室和收缩的间隙通道中为亚声速流动状态,而在扩张的密封空腔中为超声速流动状态;在弹底火药气体温度为2 166.5 K、密度为360 kg/m3、压强为242.9 MPa的条件下,经过密封装置的密封后,密封空腔内的残余气体压强为0.49 MPa。试验所得密封空腔内的最大压强为0.18 MPa,与模型计算结果基本吻合。Abstract: The present study addresses the potential damage that might be suffered by the test device of a projectile body as a result of the gun powder gas coming from the clearance opened in the projection base for installing an ejection device for the recovery of the penetration data. A seal structure of the clearance was designed and a model was built for analyzing the dynamical state of the gunpowder gas flowing in the clearance. The chamber pressure and the seal cavity pressure of the seal structure were tested by experiment. The results show that the gas in the gunpowder combustion was a compressible gas, flowing in a subsonic condition in the powder chamber and the contracting clearance channel, but in a supersonic condition in the expansion seal cavity. When the gunpowder gas in the projection base was sealed by the seal structure and kept at a temperature of 2 166.5 K, a density of 360 kg/m3 and a pressure of 242.9 MPa, the pressure for the residual gas in the seal cavity was 0.49 MPa. The maximum pressure of the test curve in the seal cavity pressure was 0.18 MPa, fairly consistent with the theoretical results. These results can serve as reference for fabricating the seal of the high temperature and high pressure gunpowder gas in the projection base.
-
Key words:
- projection base /
- gunpowder gas /
- flow parameter /
- gas dynamics /
- clearance seal /
- labyrinth seal /
- seal ring
-
当材料表面机械加工等微缺陷经历冲击波作用后,会产生速度显著高于自由面的微米级颗粒,这就是所谓的微喷射现象[1]。微喷射物质的质量、速度、空间分布会受到样品表面粗糙度、冲击波强度和加卸载历程、充气条件等一系列因素的影响[2-6],定量给出微喷射物质分布与影响因素之间的关系是研究微喷射形成机理的重要途径。基于动量守恒原理的Asay膜方法[7-8]在研究微喷射物质方面取得了巨大成功,被广泛应用,并被不断改进[9-10],相继出现了测量更高密度的Asay窗[8, 11-12]和尺寸更小的Asay膜[10]。
传统Asay膜方法测量微喷射物质的原理如图1所示,高速飞片撞击样品后在样品自由面产生微喷射,放置与样品自由面有确定间距的薄膜(Asay膜)用于收集样品自由面的微喷射物质,微喷射物质与薄膜碰撞后发生动量传递,并一起运动。采用光子多普勒速度计(photonic Doppler velocimetry, PDV)精确测量薄膜的运动过程,并按照如下4个基本假设,可获得微喷射物质的速度和质量分布。(1)瞬时喷射假设:当冲击波在自由面卸载时,全部微喷射物质是在足够短时间内产生的,即可近似认为微喷射是瞬时的,由此通过薄膜与样品自由面的距离及微喷射物质穿过该距离的时间,可计算碰撞到薄膜上的微喷射物质的速度。(2)微喷射物质分布均匀假设:即平板运动可以采用一维运动处理。(3)完全非弹性碰撞假设:微喷物质与薄膜碰撞后无反向溅射。(4)薄膜刚体运动假设:即不考虑冲击波在薄膜中的传播。
基于上述假设,可计算不同时刻到达薄膜的微喷射粒子群的速度ve[9]:
ve=(d0+∫t0vf dt)/(d0+∫t0vf dt)tt (1) 式中:d0为薄膜到样品的初始距离,t的时间零点为冲击波从自由面卸载的时刻,vf为薄膜的运动速度。
再利用动量守恒,可计算任意时刻到达薄膜的微喷射物质的质量me:
me=mf∫vf0t dvfd0−∫vf0t dvf (2) 式中:mf为薄膜的质量。对各个时刻的me积分,则得到作用于薄膜上的微喷射物质总量分布。前述各质量都为单位面积质量,即面密度。
随着对微喷射理论认识的深入,复杂加卸载条件下的微喷射现象被广泛关注,特别是样品经历连续2次冲击加载后的微喷射物质是目前较活跃的研究方向[4, 13-14]。在样品经历2次加载的情况下,自由面及微喷物质的速度随时间变化的典型状态如图2所示,首次加载后微喷射物质的分布和Asay膜常用场景是一致的。但是经历第2次加载后,首次加载后微喷粒子群速度ve1和二次加载后微喷粒子群速度ve2的值存在跳变,不再满足瞬时喷射假设和匀速运动假设,利用空腔距离和作用时间不能准确给出二次加载后某时刻微喷颗粒速度,即式(1)不成立,传统方法失效。因此,针对复杂加载微喷射物质测量,需要对Asay膜方法进行改进。
为了解决前述复杂加载条件下传统Asay膜方法应用失效的问题,本文中,提出基于Asay膜探头相同位置PDV测速信号[15]给出微喷物速度分布信息的改进Asay膜方法(以下简称新方法),针对PDV探头对微喷物质穿透深度有限、微喷物速度上下限提取存在偏差等实际应用中的问题,采用数值实验对新方法的测量偏差及原因进行分析,通过轻气炮实验对新方法与传统Asay膜方法的测量效果进行对比评估。
1. 基本方法
1.1 基于PDV测速的Asay膜数据处理
如图3所示,t0时刻微喷射物质开始作用在Asay膜上,此时自由面与Asay膜的距离用l表示,h为膜的厚度,
ˉve 为某时刻微喷物质的平均速度;t为微喷物质作用于膜上的任意时刻,其中se为微喷射物质稠密主体区的运动距离,l′ 为未作用于膜上的微喷物质的厚度,d为作用于膜上的微喷物质压缩后的厚度,sf为膜的运动距离,vf为膜的运动速度。在t时刻,假设有n个颗粒作用于Asay膜上,且质点间完全非弹性碰撞,根据动量守恒,微喷物质的速度和质量满足:
{(me+mf)vf=n∑i=1(meivei)me=n∑i=1mei (3) 测试过程中不可能追踪每个颗粒的速度和质量,假设t时刻微喷射粒子群的速度ve和质量me满足:
n∑i=1(meivei)=meve (4) 在单颗粒速度一致的情况下,式(4)严格成立。结合式(3)~(4)以及图3所示的运动关系:
Δl=l−l′=se+d−sf (5) 并假设作用在Asay膜上的微喷射物质被压实,可以给出微喷物质的总量及分布满足关系:
{me=mfvfve−vfΔl=ve(t−t0)+1ρ0mfvfve−vf−∫t0vfdt (6) 在实际测量过程中,
mf、vf和ve 都可以通过相应测试方法给出,Δl 无需初始空腔距离l即可求得,因此通过式(6)便可以获得微喷物质的总量及分布。1.2 Asay膜数值实验
为评估上述方法的数据处理不确定度,基于完全非弹性碰撞及瞬时喷射假定推导了Asay膜响应模型,如图4所示。为简化推导过程,假定单颗粒质量mp相同(可证明颗粒质量有分布,但不影响结果),忽略碰撞产生的动量矩。令喷射产生后t0时刻第1个颗粒碰膜,可获得如下的颗粒速度分布和膜速度的控制方程组分别为:
mpvp,i=(mp+mf)vf,i,tf,i=0,sf,i=0i=1 (7) {vf,i={[(i−1)mp+mf]vf,i−1+mpvp,i}/(imp+mf)dti=(t0+ti−1)(vp,i−1−vp,i)/(vp,i−vf,i−1)tf,i=tf,i−1+dtisf,i=vf,i−1dti+sf,i−1i=2,3,⋯,N (8) 式中:vp,i为第i个颗粒的速度,vf,i为第i个颗粒碰膜后的膜速,tf,i为第i个颗粒碰膜的时间,sf,i为第i个颗粒碰膜后膜的运动距离。
通过生成有随机速度分布的颗粒群,基于每个颗粒都符合匀速运动的假设,则可以获得t0时刻颗粒空间速度位置分布。将其代入式(4),则可以获得vf-t曲线,ve-t,即Asay膜数值实验数据。
1.3 基于气体炮飞片加载对比实验
为验证上述新Asay膜方法在实际应用中的效果,设计了对比实验,如图5所示,气体炮发射搭载飞片的弹丸,撞击样品产生微喷射。相较于传统的Asay膜方法,新方法增加了直接测量微喷物速度的PDV探头。由此可以获得的实验数据包括Asay膜测速结果和PDV测微喷射速度结果,结合样品与Asay膜探头的初始距离及Asay膜的面密度,即可分别采用传统Asay膜方法和本文中提出的新Asay膜方法对微喷物的面密度进行分析,对比评估新方法的测量效果。
2. 结果与讨论
2.1 数值实验
采用Matlab随机函数,生成了3种具有典型特征的颗粒群速度分布,颗粒群运动(速度方向为水平向右)t0时间后,其分布如图6所示,其中指数分布(exponential distribution, ED)和反指数分布(anti-exponential distribution, AED)都是利用Matlab的指数分布随机函数生成的,线性分布(linear distribution,LD)则是利用Matlab的均匀分布随机函数生成的。ED分布中颗粒主要集中在低速区,使用的生成函数为:vp,i=2350+random(‘exp’, 2400, 1, N)/50;而AED分布中颗粒集中在高速区,使用的生成函数为:vp,i=2650+random(‘exp’, 2400, 1, N)/50。LD分布中每种速度的颗粒数量基本一致,使用的生成函数为:vpi=2350+300×rand (N, 1)。以上生成函数中‘exp’为生成方法,N为颗粒总数,在生成过程中N=10 000,总面密度为100 mg/cm2,单颗粒等效面密度为0.01 mg/cm2。本文中讨论的碰膜时刻喷射粒子速度分布满足高速粒子在前、低速粒子在后这一基本特征,忽略了在复杂加载后较短时间内粒子速度分布不满足该特征的特殊情形,而常见应用场景中都是在距离复杂加载较长时间后进行测量的,因此本文中暂不对这一特殊情形进行讨论。
当设定膜的面密度为160 mg/cm2时,利用控制方程(8)获得的膜速曲线vf-t如图7(a)所示。图7(b)~(c)则分别展示了理论ve-t、Me-t曲线,Me代表各个时刻作用于膜上的微喷物累积面密度。
精确反演面密度分布时,准确的vf-t曲线和ve-t曲线缺一不可,而在实际复杂加载实验条件下,几乎不可能获得精确的ve-t曲线。以下将讨论使用近似ve-t曲线的处理方法。
(1)定值法:在实验测量中,直接测量微喷物速度的PDV探头至少可以给出微喷物速度分布中的某个速度值
ˉve ,使用ˉve -t进行数据处理。(2)线性插值法:假定已知速度上限(t1, ve1)和下限(t2, ve2),利用两点线性插值给出ve-t曲线。
分别采用ve取定值(2350、2450、2550、2650 m/s)及取速度上限(t1, 2650 m/s)和下限(t2, 2350 m/s),结合3种分布下的t1、t2值,线性插值给出ve-t曲线。利用式(6)对数值实验数据处理分析结果的Me-t曲线及数据处理获得的面密度值Me与理论值Met的相对偏差如图8~10所示,其中Met通过累积数值实验中每一时刻碰膜的粒子质量精确给出。
Asay膜方法常用于较大面密度(10 mg/cm2以上)微喷物质定量测量,以下将针对实际使用测量范围对新方法的测量误差进行评估分析。
当微喷颗粒呈指数分布时,除ve取值为2650 m/s时最大偏差达到15%外,其余3种取值情况下在5~100 mg/cm2面密度范围内的偏差都小于10%;当微喷颗粒呈反指数分布时,当ve取值为2350和2450 m/s时,偏差大于10.0%,分别为18.5%和11%,其余2种取值时在5~100 mg/cm2面密度范围内的偏差都小于10%;当微喷颗粒呈线性分布时,除ve取值为2350 m/s时最大偏差达到14%外,其余3种取值情况下在5~100 mg/cm2面密度范围内的偏差都小于10%。
使用线性插值方法对3种微喷物质分布的数值实验结果进行数据分析都能取得较好的效果,即在5~70 mg/cm2面密度范围内的偏差不超过5%,在5~100 mg/cm2面密度范围内的偏差不超过10%。
由上述分析可见,在可以给出微喷物质速度分布的上限和下限的情况下,无论微喷物质分布如何,采用PDV测速结合Asay膜速度进行数据分析,可以获得较高精度的微喷物质分布和累积面密度结果。然而,实验条件下,PDV测速探头可以给出微喷物质速度上限(t1, ve1),但由于PDV测速探头使用的激光对微喷物穿透深度存在极限,通常无法给出真实的速度下限。
针对PDV信号无法给出微喷物质速度下限的情况,可以先根据定值法给出初步的Me-t分布曲线,然后根据PDV探头穿透深度在Me-t曲线上给出极限穿深面密度Me, l时的t2,并取此时PDV测速信号对应的速度下限为ve2。根据PDV测速探针使用的1550 nm激光在典型微喷射物质中穿透深度约20 mg/cm2,取Me,l为10、20和40 mg/cm2,利用数值实验结果对考虑实验条件下的PDV探头穿透深度进行测量偏差分析,分析使用的(t1, ve1)、(t2, ve2)如表1所示。
表 1 不同数值实验条件下的数据处理参数Table 1. Data processing parameters under different numerical experiment conditions面密度
分布情况t1/μs t2/μs ve1/(m·s−1) ve2/(m·s−1) Me,l=10 mg/cm2 Me,l=20 mg/cm2 Me,l=40 mg/cm2 Me,l=10 mg/cm2 Me,l=20 mg/cm2 Me,l=40 mg/cm2 指数 0 2.950 3.355 3.785 2772.70 2460.83 2426.11 2393.23 反指数 0 0.046 0.101 0.240 2649.99 2644.86 2639.21 2625.85 线性 0 0.268 0.587 1.251 2649.99 2620.40 2588.02 2528.43 图11给出了假定PDV探头穿透深度分别为10、20、40 mg/cm2时数据分析结果与理论值的偏差,3种微喷射物质分布情况下,随着穿透深度的增大,数据分析结果与理论值的偏差都呈现出减小趋势,且数据分析结果较理论值都偏大;经分析,如图7(b)所示,ve-t曲线是下凹曲线,采用曲线靠前部分的两点进行线性插值后,用于计算碰膜微喷射质量的ve将被低估,且低估值越来越大,由式(6)可知,ve偏小将导致me偏大。
当微喷物质分布集中靠近Asay膜方向(ED)或样品方向(AED)时,采用线性插值法数据处理的面密度结果在5~100 mg/cm2范围内与理论值偏差小于10%;而当微喷射物质为均匀分布(LD)时,采用线性插值方法数据处理的面密度结果在5~100 mg/cm2范围内的与理论偏差小于20%。针对LD分布时偏差较其余2种分布大的现象进一步分析,通过图7(b)可知,处于插值点附近区域的ve与理论值差较小,而越远离插值区域偏差越大,ED和AED分布中主要微喷射物质与插值区域更靠近,而LD分布中各个区域微喷物都一样,因此,新的Asay膜方法对均匀分布的微喷射累积面密度测量误差较大。
2.2 轻气炮实验
2.1节中通过数值实验分析了新Asay膜方法进行微喷物面密度测量的偏差及其原因,为检验新方法在实验中的应用效果,利用轻气炮开展了验证实验。采用飞片冲击表面预置微米粉末样品开展实验,如图5所示,飞片为尺寸
∅ 13.2 mm×2 mm的304不锈钢,飞片速度为1923.1 m/s,样品为∅ 20 mm×2 mm不锈钢片表面均匀预置粒径分布为3~6 μm的球形锡颗粒,粉末初始实测面密度约44 mg/cm2。PDV探头中心和Asay膜探头中心位于同一圆周上。测得的Asay膜速度及相应的PDV测速数据如图12所示。结合测试结果及已知的初始条件,分别利用传统的空腔法(传统Asay膜法)和本文中提出的新方法(新Asay膜法)对实验数据进行了处理,相应的处理参数如表2所示。其中新方法中速度上限(t1, ve1)和下限(t2, ve2)的选取分别为:膜有速度响应时刻为t1,膜速陡升时刻为t2,t1时刻对应的频谱速度上限为ve1,t2时刻对应的速度频谱下限为ve2。
表 2 两种Asay膜数据处理方法使用的参数Table 2. Parameters used in two different Asay foil data processing methods方法 l/mm t0/μs mf/(mg·cm−2) t1/μs t2/μs ve1/(m·s−1) ve2/(m·s−1) 传统Asay膜法 19.5 28.33 156 新Asay膜法1 156 34.05 39.33 4180 2026 新Asay膜法2 156 34.05 39.33 4380 1826 新Asay膜法3 156 34.05 39.33 3980 2226 表2中新方法选取了3组数据,其中新Asay膜法1是按图12(b)中PDV频谱对应的t1、t2时刻分别提取速度上限和下限。考虑到在频谱上提取速度时引入的误差,新Asay膜法2和新Asay膜法3中考虑了2种极端情况,即频谱数据提取速度上下限存在±200 m/s的误差。图13(a)对比了分别采用传统Asay膜方法和本文中提出的新Asay膜方法的数据处理结果,在考虑PDV频谱速度上下限误差的情况下,2种方法给出的微喷物质分布及累积面密度有较好的一致性。需特别说明的是,Asay膜方法测得的面密度比预置粉末的面密度小的主要原因是粉末未完全从304钢样品表面脱落。图13(b)分析了2种方法的测量结果的相对偏差,采用新方法测得的累积面密度在初始阶段较传统方法测得的小,但随着累积面密度的增大,新Asay膜方法与传统Asay膜方法的差异逐渐减小,并且新Asay膜法1和新Asay膜法2得到的累积面密度较传统方法的大。综合来看,在累积面密度大于5 mg/cm2的情况下,在考虑PDV穿透深度、微喷物速度上下限提取偏差等因素后,使用新Asay膜方法得到的面密度空间分布和累积面密度相对传统方法的偏差都小于20%。
3. 结 论
针对复杂加载条件下微喷射物质测量这一具体问题,分析了传统Asay膜方法存在的不足,特别是在某些多次加载条件下无法给出传统Asay膜数据处理方法所必需的初始空腔距离及自由面起动时间等参数,提出了基于PDV测速结合Asay膜测速结果进行微喷射物质质量分布分析的新方法。采用传统方法和新方法对数值实验及气炮加载实验数据进行对比分析,验证了新方法的有效性,并分析了测量偏差,获得的主要认识如下。
(1)基于PDV测速结合Asay膜测速处理微喷射物质质量分布在无法获得自由面起跳时间及膜与自由面初始空腔距离时有效,扩充了Asay膜方法的适用范围;新Asay膜方法的测量值在累积面密度大于5 mg/cm2时,相对传统Asay膜方法测量偏差小于20%,即与传统Asay膜方法测试精度基本相当。
(2)在考虑PDV探头对微喷射物质穿透深度有限的条件下,新Asay膜方法对微喷射物质空间分布具有一定的敏感性,当微喷射物质集中于高速区或低速区时,微喷物累积面密度测试结果与理论值有较小的偏差;当微喷射物质均匀分布时,结果与参考值偏差较大;对于已知被测对象大致分布的情况,新Asay膜方法通过改变对ve的拟合方式,理论上测试精度仍可提升。
(3)微喷射速度测试PDV探头要求与Asay膜探头的被测对象具有相同物理状态或处于相同位置,因此,当微喷射状态空间分布极不均匀时,如果无法保证2种探头的被测对象基本一致,该方法的测试误差可能被放大,应当慎用。
本文中,仅对针对激光干涉测速信号和Asay膜信号的本身特点开展了分析,分析中不包括Asay膜探头自身的测量偏差、实际物理过程偏离基本假设等。针对此类因素,后续将开展进一步研究。本文中讨论的喷射粒子碰膜时都具有高速粒子先碰膜、低速粒子后碰膜的特点,针对第2次加载强度较低或间隔时间较长时,理论上存在高速粒子先碰膜、低速粒子后碰膜这一特殊情况,后续将进一步开展数值实验分析,拓宽本文方法的应用场景。
-
[1] 王茂林, 黄强, 张国平, 等.衬套径向间隙部位火药气体压力计算[J].火炮发射与控制学报, 2009, 30(3):49-52. doi: 10.3969/j.issn.1673-6524.2009.03.013Wang Maolin, Huang Qiang, Zhang Guoping, et al. Pressure calculation of powder gas in bush radial clearance[J]. Journal of Gun Launch and Control, 2009, 30(3):49-52. doi: 10.3969/j.issn.1673-6524.2009.03.013 [2] 张浩, 周彦煌.一种瞬态高压条件下的密封设计[J].润滑与密封, 2006, 31(5):132-134. doi: 10.3969/j.issn.0254-0150.2006.05.043Zhang Hao, Zhou Yanhuang. A seal design under transient high pressure condition[J]. Lubrication Engineering, 2006, 31(5):132-134. doi: 10.3969/j.issn.0254-0150.2006.05.043 [3] 张浩, 陆欣, 余永刚, 等.某口径埋头弹火炮的密封与装药设计[J].兵工学报, 2006, 27(4):630-633. doi: 10.3321/j.issn:1000-1093.2006.04.012Zhang Hao, Lu Xin, Yu Yonggang, et al. Design for the seal system and charge of a CAT gun[J]. Acta Armamentarii, 2006, 27(4):630-633. doi: 10.3321/j.issn:1000-1093.2006.04.012 [4] 张讯, 张相炎.新型组合式炮膛密封结构仿真研究[J].弹道学报, 2009, 21(1):63-66. http://d.old.wanfangdata.com.cn/Periodical/ddxb200901016Zhang Xun, Zhang Xiangyan. Simulation research on new combined barrel seal structure[J]. Journal of Ballistics, 2009, 21(1):63-66. http://d.old.wanfangdata.com.cn/Periodical/ddxb200901016 [5] Fu Qiang, Chen Longmiao, Qian Linfang, et al. Optimum design for seal structure of a new-type artillery based on MOGA-Ⅱ and RSM[C]//Proceedings of the 4th International Conference on Mechanical Engineering and Mechanics. Suzhou, 2011: 164-168. [6] 赵威, 何永, 陶齐冈, 等.一种高压试验装置增压密封结构设计及仿真[J].兵工自动化, 2013, 32(11):57-59. doi: 10.7690/bgzdh.2013.11.015Zhao Wei, He Yong, Tao Qigang, et al. Design and simulation for supercharged sealing structure of high voltage test device[J]. Ordnance Industry Automation, 2013, 32(11):57-59. doi: 10.7690/bgzdh.2013.11.015 [7] 陈伟, 张敬兵, 陈瑶, 等.一种新型试验装置密封结构仿真研究[J].机械制造与自动化, 2014, 43(5):23 -25, 29. doi: 10.3969/j.issn.1671-5276.2014.05.010Chen Wei, Zhang Jingbing, Chen Yao, et al. Simulation research on new seal structure of test device[J]. Machine Building and Automation, 2014, 43(5):23 -25, 29. doi: 10.3969/j.issn.1671-5276.2014.05.010 [8] 崔凯波, 秦俊奇, 狄长春, 等.变化油压冲击下火炮制退机唇形密封圈的动密封性能[J].润滑与密封, 2014, 39(8):24-30. doi: 10.3969/j.issn.0254-0150.2014.08.005Cui Kaibo, Qin Junqi, Di Changchun, et al. Dynamic seal performance of lip ring of artillery recoil mechanism under changeable hydraulic impact[J]. Lubrication Engineering, 2014, 39(8):24-30. doi: 10.3969/j.issn.0254-0150.2014.08.005 [9] 景思睿, 张鸣远.流体力学[M].西安:西安交通大学出版社, 2001:6-9. [10] 张兆顺, 崔桂香.流体力学[M].2版.北京:清华大学出版社, 2006:8-10. [11] 孔珑.工程流体力学[M].4版.北京:中国电力出版社, 2014:151-174. [12] 张连玉, 汪令羽, 苗瑞生.爆炸气体动力学基础[M].北京:北京工业学院出版社, 1987:104-129. [13] 王保国, 刘淑艳, 黄伟光.气体动力学[M].北京:北京理工大学出版社, 2005:170-184. [14] 韩占忠, 王敬, 兰小平.FLUENT:流体工程仿真计算实例与应用[M].北京:北京理工大学出版社, 2013:84-91. [15] 黄守龙, 陆家鹏, 徐诚.直通式迷宫瞬态压力特性的实验研究[J].弹道学报, 1994, 6(4):7-11. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400029597Huang Shoulong, Lu Jiapeng, Xu Cheng. Experimental study on transient features for straight-through labyrinth[J]. Journal of Ballistics, 1994, 6(4):7-11. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400029597 [16] 塔鲁达纳夫斯基.非接触密封: 间隙密封与迷宫密封的原理和应用[M].李均卿, 刁元康, 译.北京: 机械工业出版社, 1986: 77-95. [17] 吴晓铃, 袁丽娟.密封设计入门[M].北京:化学工业出版社, 2013:56-74. [18] 陆家鹏, 谭兴良, 雷志义.自动武器气体动力学[M].北京:兵器工业出版社, 1991:83-84. [19] Pérez-Ràfols F, Larsson R, Almqvist A. Modelling of leakage on metal-to-metal seals[J]. Tribology International, 2016, 94:421-427. doi: 10.1016/j.triboint.2015.10.003 [20] 蔡仁良.流体密封技术:原理与工程应用[M].北京:化学工业出版社, 2013:65-77. [21] Kima G H, Leeb Y S, Yanga H L. A new design concept of metal O-ring seal for long-term performance[J]. Vacuum, 2016, 123:54-61. doi: 10.1016/j.vacuum.2015.10.014 [22] 崔建昆.密封设计与实用数据速查[M].北京:机械工业出版社, 2010:135-165. -