Loading [MathJax]/jax/output/HTML-CSS/jax.js
  • ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

中心点火火焰在药床中传播规律的试验研究

薛绍 陶如意 王浩 程申申

黄勇, 鲁长波, 安高军, 熊春华, 解立峰. 柴油爆炸性能外场实验研究[J]. 爆炸与冲击, 2015, 35(4): 482-488. doi: 10.11883/1001-1455(2015)04-0482-07
引用本文: 薛绍, 陶如意, 王浩, 程申申. 中心点火火焰在药床中传播规律的试验研究[J]. 爆炸与冲击, 2021, 41(11): 112101. doi: 10.11883/bzycj-2021-0030
Huang Yong, Lu Chang-bo, An Gao-jun, Xiong Chun-hua, Xie Li-feng. Experimental research on explosion performance of diesel fuel in the external field[J]. Explosion And Shock Waves, 2015, 35(4): 482-488. doi: 10.11883/1001-1455(2015)04-0482-07
Citation: XUE Shao, TAO Ruyi, WANG Hao, CHENG Shenshen. Experimental research on the law of flame spreading in the charge bed of a central ignition tube[J]. Explosion And Shock Waves, 2021, 41(11): 112101. doi: 10.11883/bzycj-2021-0030

中心点火火焰在药床中传播规律的试验研究

doi: 10.11883/bzycj-2021-0030
详细信息
    作者简介:

    薛 绍(1991- ),男,博士研究生,xue18761686700@163.com

    通讯作者:

    陶如意(1978- ),女,博士,副教授,tao801801@163.com

  • 中图分类号: O381

Experimental research on the law of flame spreading in the charge bed of a central ignition tube

  • 摘要: 为了研究中心点火管火焰在药床中的传播规律,设计了可视化模拟试验平台,开展了不同点火药量、不同装药结构的中心点传火试验。采用高速图像采集系统记录了中心点火管火焰在药床中的传播过程,采用瞬态压力记录仪记录膛内压力的时空变化。结果表明,点火药量为20 g时,出火时间为0.6 ms;点火药量为30 g时,出火时间为1.5 ms;杆状装药床的传火时间平均为2.2 ms,粒状装药床的传火时间平均为3.4 ms,而杆粒混装药床的传火时间为3.1 ms。可见,点火药量对药床出火时间影响显著,较大的点火药量导致药床出火时间延长;不同装药床结构传火性能差异较大,单一杆状装药床传火性能优于单一粒状装药和杆粒混装药床,并且粒状装药床易形成气体壅塞,膛内会出现明显的压力波动现象;根据火焰传播时序位置点,利用一阶指数衰减函数拟合建立了火焰传播过程数学模型,拟合优度大于0.98。
  • 薄壁圆管凭借其自身结构简单,变形、破坏模式稳定,轴压峰值力与平均力的可设计性以及能量吸收可控等优点,被广泛地应用于抗冲击领域,如飞机坠撞防护、缓冲器以及其他的工程领域等[1-2]。薄壁圆管在轴向冲击压缩载荷作用下的屈曲问题长期以来备受关注。Alexander等[3]针对理想刚塑性材料制成的薄壁圆管,应用具有塑性铰的简化轴对称变形模式,推导出了轴对称条件下的轴压载荷下薄壁圆管的平均力。该变形模型获得了学者的广泛认可并得到了大量的推广应用。Abramowicz等[4]基于Alexander等[3]的研究结果,对薄壁圆管轴对称变形时的折皱长度公式进行了修正,试验结果表明,修正后的公式可以较好地预测轴对称变形以及金刚石变形模式下的薄壁圆管的折皱长度。Karagiozova等[5-6]对受轴向冲击的弹塑性材料薄壁圆管的峰值载荷进行了研究,分别获得了Mises屈服条件和Tresca屈服条件下的峰值载荷。Lu等[7]在Alexander等[3]的基础上,提出了一种新的运动位移场,求解了更加合理的薄壁圆管折皱长度公式,并改进了轴对称变形下的载荷位移曲线。这些重要的研究为理解轴向冲击圆管的变形行为奠定了理论基础。

    结构的冲击承载安全性是结构设计的重要组成部分,但是大结构的全尺寸冲击试验极其耗时耗力,需要付出大量的人力物力和财力,因此一般使用比例模型试验来替代结构全尺寸原型[2, 8-9]。Jones[1]通过量纲分析,建立了结构冲击问题的经典相似理论,由于该理论基于简单的几何相似关系,因而广受重视。Oshiro等[10-11] 采用新的量纲分析体系克服经典相似律解决材料相似时的理论困难,提出一种通过修正冲击速度或密度的方式,来适应应变率敏感性以及比例模型与原型使用不同材料时的畸变问题的相似方案。王帅等[12-13]通过建立应变率区间上比例模型预测的流动屈服应力与原型流动屈服应力的最佳逼近关系,获得了速度比例因子的直接相似关系,从而扩展了相似律的应用技术。李肖成等[14]通过调整薄板结构的厚度,修正了比例模型与原型采用不同材料时引起的畸变问题。秦健等[15]基于VSG量纲体系,研究得到了加筋板结构在水下爆炸冲击波作用下模型和原型的动态响应相似关系。

    然而,在实际工程应用中还存在加工制造等难题,通常,薄壁圆管的厚度不一定能够和半径以及长度按照相同的比例因子进行缩放,这就会产生畸变模型,从而导致原有的相似律失效。Alves等[16]提出单独考虑几何畸变以及应变率效应,并通过经验方法不断迭代几何畸变对缩比相似影响所占的比重来修正速度比例因子。Mazzariol等[17]通过塑性弯矩的无量纲数改进了经验法,提出了一种解决薄板和梁结构的几何畸变问题的新方法。Wang等[18-19]基于定向量纲的分析方法,从理论上克服了几何畸变的相似律难题,针对薄板冲击模型系统表达了几何畸变的相似体系。尽管这些研究在解决几何畸变相似律的问题上取得了重要的理论进步,但是针对薄壁圆管的几何畸变修正方法的研究分析并未开展,而且在实际的受轴向冲击载荷作用的薄壁圆管缩比相似分析中厚度畸变几乎不可避免。

    本文中针对受轴向冲击载荷作用的薄壁圆管,基于方程分析与量纲分析,以载荷、位移和能量为关键设计参数,推导薄壁圆管几何畸变的动态响应相似律,此方法可推广至考虑材料应变率以及应变硬化效应和几何畸变问题耦合时的比例模型与原型物理量对应关系的描述。

    在结构冲击相似分析中,一般使用基于一组无量纲数的比例因子来预测比例模型的力学行为,如VSG量纲体系[10]或DLV量纲体系[18],其比例因子如表1所示。其中,比例模型和原型相关的物理量分别用下标m和p表示,βK=Km/Kp表示比例模型和原型相关物理量K的比值,例如βv=vm/vp表示比例模型和原型速度的比值,比例模型和原型的几何比例因子β=Lm/Lp

    表  1  纯几何相似律比例因子[10-11]
    Table  1.  Scaling factors by pure geometric similarity[10-11]
    变量比例因子变量比例因子
    长度L[10] β=Lm/Lp 位移δ[10] βδ=β
    密度ρ[10] βρ=ρm/ρp 应力σd[10] βσd=βρβ2v
    速度v[10] βv=vm/vp 应变ε[10] βε=1
    质量m[10] βm=βρβ3 应变率˙ε[10] β˙ε=βv/β
    时间t[10] βt=β/βv 载荷P[11] βP=βρβ2β2v
    加速度a[10] βa=β2v/β 动能Ek[11] βEk=βρβ3β2v
    下载: 导出CSV 
    | 显示表格

    基于表1中的相似关系,目前针对受轴向冲击载荷作用的薄壁圆管的相似律研究对象均为纯几何缩放的比例模型,也就是说薄壁圆管的厚度h、长度L和半径R服从同一个几何比例因子,亦即βh = βR = βL = β

    但对于薄壁圆管而言,由于其自身几何尺寸的特性,厚度方向的尺寸与半径和长度方向的尺寸完全不在一个量级,厚度方向若要按照长度和半径的几何比例因子进行缩放,会导致比例模型加工的成本远远超过加工一个原型试验件的费用。并且,受到加工限制,其厚度方向的尺寸h并不一定能够与长度L和半径R按照同一个几何比例因子进行缩放,这时就会产生几何畸变模型,这时,βh = βR = βL = β这一基础条件会遭到破坏,将会导致表1中的相似律不能全部满足,这就需要重新推导薄壁圆管的相似律。

    考虑到薄壁圆管的结构响应特性,以薄壁圆管压溃过程中的平均载荷、位移和能量这3个响应变量作为相似设计的关键参数[20]。影响薄壁圆管压溃特性的参量有:厚度h、半径R、长度L、冲击速度v、冲击质量M、材料的密度ρ、流动屈服应力σd

    根据薄壁圆管理论,在推导过程中做如下假设:

    (1)受轴向冲击载荷作用的薄壁圆管的变形模式服从Alexander所给出的轴对称变形模式,如图1所示;

    图  1  轴向受压薄壁圆管的轴对称压溃模式
    Figure  1.  Axisymmetric crushing mode of the thin-walled cylindrical shell under axial compression

    (2)薄壁圆管壁厚与半径及长度相比为小量,服从薄壳理论;

    (3)与弹性变形相比,塑性变形占据主导地位,弹性变形可以忽略不计;

    (4)薄壁圆管使用理想弹塑性材料,材料的流动屈服应力为一定值;

    (5)薄壁圆管受到重撞击物的低速冲击;

    (6)定义基本几何比例因子:β = βR = βL

    当薄壁圆管发生图1所示的准静态压溃模式时,Alexander[3]给出了该模式下的薄壁圆管每产生一个完整折皱时的平均压溃力P和折皱长度l的表达式:

    {P=2(πh)3/2R1/2σd/31/4l=(πRh/3)1/2 (1)

    式中:σd为材料的流动屈服应力,对于理想弹塑性材料,其为一常值。

    当薄壁圆管受重撞击物的低速冲击时,可忽略惯性力的影响,薄壁圆管的动态压溃过程可近似视为准静态,此时式(1)能够适用于动态载荷下的薄壁圆管的轴对称压溃问题[1, 4, 21]

    薄壁圆管的质量m为:

    m=2ρπRhL (2)

    当薄壁圆管受到轴向冲击载荷作用时,忽略摩擦以及温度软化导致的能量损失,那么,在整个薄壁圆管压溃过程中的能量转化就是结构的动能与塑性耗散能的相互转化。根据能量守恒定律,塑性耗散能可等效视为外力P在整个薄壁圆管的压缩行程δ中所做的功,也就是说外力P所做的功WP应等于薄壁圆管所受到的冲击能量Ek [22]

    WP=EkPδ=12Mv02 (3)

    式中:M为冲击块的质量,v0为撞击物的初始速度。

    将式(2)代入式(3)可得:

    Pδ=πρRhLv20Mm (4)

    将式(1)代入式(4),则式(4)进一步可变形为:

    δL=31/42πρv20σd(Rh)0.5Mm (5)

    通过式(5)可以得到,发生轴对称变形的薄壁圆管受轴向冲击载荷作用时的无量纲位移、质量和应力分别为:

    Πδ=δL,Πm=Mm,Πσd=ρv2σd(Rh)0.5 (6)

    为了使比例模型可以准确预测原型的动态响应,原型与比例模型的无量纲数应相等,也就是说:

    (Πδ)m(Πδ)p=βδβ=1βδ=β (7)
    (Πm)m(Πm)p=βMβm=1βM=βm=βρβ2βh (8)
    (Πσd)m(Πσd)p=βρβ2vβσd(ββh)1/2=1βv=βρβσd(ββh)1/2 (9)

    对于遵循图1所示的压溃模式的薄壁圆管来说,在折皱的轴对称部分ab段,距离a点径向距离为x的微元伸长量为:

    ds=2π(R+x)2πR=2πx (10)

    对应该微元处的环向应变为:

    εθ=2πx2πR=xR (11)

    不难发现,图1ab段各点处的环向应变与x成正比,那么对于一个完整折皱,x∈[0, l],其平均环向应变为:

    εθ=120+lR=l2R (12)

    平均环向应变率为:

    ˙εθ=εθt=l2Rt (13)

    对发生轴对称变形的薄壁圆管,与环向应变相比,其轴向应变十分小,而薄壁圆管壁厚很薄,可近似认为其径向应变为零,因此可以用其环向的平均应变εθ和平均应变率˙εθ替代等效应变ε与等效应变率˙ε,即:

    εθε,˙εθ˙ε (14)

    将式(1)代入式(12)~(13),可得应变和应变率的比例因子为:

    {βε=(εθ)m/(εθ)p=βlβ=βhββ=βh/ββ˙ε=(˙εθ)m/(˙εθ)p=βlββt=βlββvβ=(βhβ)1/2βvβ2=βvββhβ (15)

    综上所述,得到受轴向冲击载荷的薄壁圆管的缩放比例因子,见表2,当比例模型与原型满足表2中的比例因子时,比例模型可以准确地预测原型的动态响应。

    表  2  受轴向冲击的理想弹塑性薄壁圆管比例因子
    Table  2.  Scaling factors of the elastic-ideal plastic thin-walled cylindrical shell under axial impact loading
    变量比例因子变量比例因子
    长度L β=Lm/Lp位移δ βδ=β
    密度ρ βρ=ρm/ρp应力σd βσd=βρβ2vβ/βh
    速度v βv=vm/vp应变ε βε=βh/β
    质量m βm=βρβ2βh应变率˙ε β˙ε=(βv/β)βL/βh
    时间t βt=β/βv载荷P βP=βρββhβ2v
    加速度 a βa=β2v/β动能Ek βEk=βρβhβ2β2v
    下载: 导出CSV 
    | 显示表格

    实际工程中,理想弹塑性的材料是几乎不存在的,绝大多数金属材料的塑性阶段都存在着应变率效应和应变硬化效应[23],如图2所示,材料的流动屈服应力是应变和应变率的函数而非一固定值,即:

    图  2  金属材料应力-应变曲线
    Figure  2.  Stress-strain curves of metal materials
    σd=f(ε,˙ε) (16)

    将式(15)代入式(16),得到流动屈服应力的比例因子为:

    βσd=(σd)m(σd)p=f(εm,˙εm)f(εp,˙εp)=f(εpβε,˙εpβ˙ε)f(εp,˙εp)=f(εpβh/β,˙εpβvβh/β3)f(εp,˙εp) (17)

    可见,式(17)所得到的 βσd在整个薄壁圆管的压溃过程中是一个不断变化的值。但是若要使比例模型准确地预测原型的动态响应,原型与比例模型的流动屈服应力在整个冲击过程中的应变区间[εa,εb]与应变率区间[˙εa,˙εb]中应始终满足(σd)m=(σd)p。而由于材料的应变率和应变硬化效应引起的缩比相似的畸变问题,使这一条件几乎无法满足。

    对此,可在应变率区间[˙εa,˙εb]和应变区间[εa,εb]上,通过(σd)m(σd)p的误差优化来减小两者的误差,从而减小应变率和应变硬化引起的误差[12]

    如果原型的应变区间为[εa,εb],应变率˙εp[˙εa,˙εb]上变化,那么比例模型预测的流动屈服应力(σd)m与原型流动屈服应力(σd)p存在一个最小误差,即:(σd)m(σd)p2=minβvR(σd)m(σd)p2成立。那么就有:

    fm(εpβhβ,˙εpβvββhβ)1βρβv2βhβfp(εp,˙εp)2=minβvRfm(εpβhβ,˙εpβvββhβ)1βρβv2βhβfp(εp,˙εp)2 (18)

    式中:βv为比例模型所预测的流动屈服应力(σd)m与原型的流动屈服应力(σd)p在整个应变区间和应变率区间上误差最小时的速度比例因子,(σd)m为速度比例因子选取为 βv时的比例模型预测的原型的流动屈服应力。·为2范数,其中应变区间[εa,εb]与应变率区间[˙εa,˙εb]可根据试验确定,当难以进行试验时可根据经验粗略选取。

    至此,我们将理想弹塑性假设下的薄壁圆管缩比相似分析拓展到了同时伴有应变率效应和应变硬化效应的一般材料薄壁圆管缩比相似分析的应用中。

    为了验证本文所推导的相似律对于结构缩比相似分析中的应变率和应变硬化效应畸变问题以及几何畸变问题的修正能力,在ABAQUS中建立了如图3所示的薄壁圆管模型进行验证分析。薄壁圆管原型的长度L=240 mm,中性面半径R=60 mm,壁厚h=1.8 mm。薄壁圆管材料选取应变率较为敏感的1006 钢,采用Johnson-Cook本构模型进行数值模拟,材料参数见表3[24]

    图  3  受轴向冲击的薄壁圆管示意图
    Figure  3.  Schematic diagram of a thin-walled cylindrical shell under axial impact
    表  3  1006 钢的材料参数[24]
    Table  3.  Material parameters of 1006 steel[24]
    ρ/(g∙cm−3)E/GPaμA/MPaB/MPaCn˙ε0/s1
    7.892070.33502750.0220.361
    下载: 导出CSV 
    | 显示表格

    薄壁圆管采用CAX4R轴对称单元进行仿真,锤头和地面则使用刚体单元。薄壁圆管自由放置于刚性地面上,顶部受到质量M=128 kg、速度v0 =12 m/s的质量块冲击。薄壁圆管与地面之间的摩擦因数为0.2,圆管内壁、外壁之间设置为无摩擦。

    分别建立了基本几何比例因子β=1/10,厚度方向的比例因子βh 分别为 0.12、0.15、0.18和0.20的比例模型。并按照表2所给出的相似律获得速度修正比例因子,定义几何畸变度η = βh /β。由表2可得理想弹塑性材料的薄壁圆管几何畸变比例模型的各动态变量的比例因子,见表4。同时为了验证本文对于应变率和应变硬化效应与几何畸变效应耦合后的修正方法,建立了考虑应变率效应和应变硬化效应后的薄壁圆管几何畸变比例模型,并按照式(18)所给出的方法获得应变率区间和应变硬化区间上的最佳速度比例因子,通过表2即可获得比例模型与原型各动态变量的比例因子,见表5。其中,应变率区间和应变区间在大量数值模拟经验的基础上分别选取为[0, 500]和[0, 0.6]。

    表  4  理想弹塑性模型几何畸变比例因子
    Table  4.  Scaling factors of geometrically-distorted models of elastic-ideal plastic material
     βη βh βM βv βt βP βδ βε
    0.11.20.120.121.04660.09550.01310.11.0954
    0.11.50.150.151.10670.09040.01840.11.2247
    0.11.80.180.181.15830.08630.02410.11.3416
    0.12.00.200.201.18920.08410.02830.11.4142
    下载: 导出CSV 
    | 显示表格
    表  5  考虑应变率效应和应变硬化效应几何畸变模型比例因子
    Table  5.  Scaling factors of geometrically-distorted models considering strain-rate sensitivity and strain hardening
     βη βh βM βv βt βP βδ βε
    0.11.20.120.121.07560.09300.01390.11.0954
    0.11.50.150.151.14420.08740.01960.11.2247
    0.11.80.180.181.20380.08310.02610.11.3416
    0.12.00.200.201.23960.08070.03060.11.4142
    下载: 导出CSV 
    | 显示表格
    2.2.1   理想弹塑性几何畸变模型结果分析

    图4为使用本文提出的受轴向冲击的薄壁圆管几何畸变修正方法后,获得的修正后理想弹塑性薄壁圆管的几何畸变比例模型与原型结构动态响应时间曲线的对比。图4(a)为锤头的速度-时间响应曲线,图4(b)为薄壁圆管首个折皱处塑性铰位置的等效塑性应变-时间响应,图4(c)、(d)分别为锤头的位移和能量-时间响应曲线,图4(e)为薄壁圆管几何畸变比例模型与原型载荷时间-曲线的对比。表6为修正后的薄壁圆管几何畸变模型所预测原型的峰值位移和平均载荷与原型响应之间的相对误差。

    图  4  原型与比例模型的动态响应-时间曲线
    Figure  4.  Dynamic response-time curves of the scale models and the prototype
    表  6  比例模型与原型的峰值位移和平均载荷相对误差
    Table  6.  Relative errors in the peak displacement and average force between the scale models and prototype
    模型δ/βδ)/mm相对误差/%P/βP/kN相对误差/%
    原型124.25871.95
    β=1/10, η=1.2121.6692.08473.432.057
    β=1/10, η=1.5120.6512.90369.483.433
    β=1/10, η=1.8118.6424.52075.855.420
    β=1/10, η=2.0117.6785.29572.250.417
    下载: 导出CSV 
    | 显示表格

    图4可以看出,修正前的不同几何畸变度的理想弹塑性薄壁圆管比例模型的响应曲线基本无法重合,各条曲线上的相应点所对应的时间与值均不相等,而修正后的理想弹塑性薄壁圆管几何畸变比例模型与原型的速度、能量、位移、等效塑性应变和载荷-时间曲线重合性明显改善,修正后的比例模型的速度和应变-时间曲线与原型的响应曲线完全重合。对于本文推导过程中的3个关键设计参数:能量Ek、位移δ和平均载荷P,修正后的比例模型的位移和能量曲线与原型的响应曲线重合性也非常好,但是修正后的比例模型与原型的载荷-时间曲线重合度不是很高。这是因为本文推导的相似律以能量Ek、位移δ和平均载荷P作为关键相似设计参数,薄壁圆管的折皱长度l与位移δ并不服从同一个比例因子,因此几何畸变比例模型与原型的折皱数量并不一致,导致修正后的比例模型和原型的载荷-时间曲线无法完全重合,但是相较修正前,其误差已经大大减小,且其整体的时间趋势一致性比较良好。因此我们又分析了比例模型与原型的峰值位移和平均载荷的误差,从表6可以看到,即使几何畸变度达到2,比例模型所预测的结果与原型的峰值位移和平均载荷的相对误差均不超过6%。

    综上所述,使用本文所推导的相似律,修正后的理想弹塑性薄壁圆管几何畸变比例模型的动态响应可以准确地预测原型的动态响应。

    2.2.2   考虑应变率效应和应变硬化效应几何畸变模型结果分析

    使用本文提出的受轴向冲击的一般材料薄壁圆管几何畸变修正方法,同时考虑应变率效应和应变硬化效应与几何畸变问题,获得的修正后比例模型预测结果与原型动态响应时间曲线的对比如图5所示。图5(a)为锤头的速度-时间响应曲线,图5(b)为薄壁圆管首个折皱处塑性铰位置的等效塑性应变-时间响应,图5(c)、(d)分别为锤头的位移和能量-时间响应曲线,图5(e)为修正前后的薄壁圆管几何畸变比例模型与原型载荷时间-曲线的对比。表7为考虑应变率效应和应变硬化效应,修正后的薄壁圆管几何畸变模型所预测原型的峰值位移与平均载荷与原型响应之间的相对误差。

    图  5  修正后的比例模型与原型的动态响应-时间曲线
    Figure  5.  Dynamic response-time curves of the modified scaled models and the prototype
    表  7  考虑应变率与应变硬化效应的比例模型位移与平均载荷相对误差
    Table  7.  Relative errors in the peak displacement and average force of the scaled models considering strain-rate sensitivity and strain hardening
    模型(δ/βδ) /mm相对误差/%(P/βP)/kN相对误差/%
    原型90.74097.79
    β=1/10, η=1.289.6191.23598.420.644
    β=1/10, η=1.589.2221.67399.471.718
    β=1/10, η=1.886.2015.002101.063.344
    β=1/10, η=2.086.2324.968105.147.516
    下载: 导出CSV 
    | 显示表格

    图5可以看出,修正前的不同几何畸变度的一般材料薄壁圆管比例模型的响应曲线基本无法重合,各条曲线上的相应点所对应的时间与值均不相等,而修正后的一般材料薄壁圆管几何畸变比例模型与原型的速度、能量、位移、等效塑性应变和载荷-时间曲线重合度明显提高。修正后的比例模型的速度和应变-时间曲线与原型的响应曲线基本完全重合。对于本文推导过程中的3个关键设计参数:能量Ek、位移δ和平均载荷P,修正后的比例模型的位移和能量曲线与原型的响应曲线重合度非常高,但修正后的几何畸变比例模型的载荷-时间曲线与原型重合度不是很高,原因如2.2.1节所述。然而随着几何畸变度的增大,修正后的比例模型的动态响应变量曲线与原型的动态变量响应曲线的重合程度逐渐降低。对比修正后的比例模型与原型的峰值位移和平均载荷,从表7中可以看到,即便几何畸变度增大至2,比例模型所预测的结果与原型的峰值位移相对误差依然不超过6%,平均载荷的误差亦不超过8%。

    可见,通过使用本文所推导出的薄壁圆管相似律,对考虑应变率效应和应变硬化效应的一般材料薄壁圆管几何畸变比例模型进行修正后,比例模型可以准确地预测原型的动态响应。

    以受轴向冲击载荷的薄壁圆管为研究对象,基于量纲分析与方程分析,推导了薄壁圆管结构的几何畸变相似律。

    (1)通过能量守恒定律,直接建立了结构几何特征与冲击速度和冲击质量的关系,解决了理想弹塑性材料薄壁圆管的壁厚h、半径R和长度L非等比例缩放的不相似问题。

    (2)建立应变区间和应变率区间上比例模型预测的流动屈服应力与原型流动屈服应力的最佳逼近关系,解决了考虑材料应变率和应变硬化效应与薄壁圆管几何畸变问题耦合时的不相似问题。

    (3)通过数值模拟验证了该方法的有效性,研究表明,使用该方法可大幅降低几何畸变模型预测结果与原结构实际响应的误差,特别是轴向冲击问题备受关注的载荷和能量,相似误差最大不超过8%。

    但该相似律只适用于几何畸变比例模型满足薄壳假设且弹性变形能占比可忽略时的薄壁圆管几何畸变相似性分析,对于不满足该条件的薄壁圆管相似性分析还需进一步研究。

  • 图  1  模拟试验系统流程图

    Figure  1.  Flow chart of the simulation experimental system

    图  2  可视化模拟装置

    Figure  2.  The visualization device for experimental investigation of propellant charge ignition

    图  3  中心点火管结构示意图

    Figure  3.  Schematic diagram of the central ignition tube

    图  4  测试现场布置

    Figure  4.  Arrangement of the test site

    图  5  不同装药结构

    Figure  5.  Different charging configurations

    图  6  不同试验方案火焰传播序列照片

    Figure  6.  High-speed photography of flame spreading through propellant charge in different cases

    图  7  方案1、3、5火焰传播位置

    Figure  7.  Flame positions along the propellant chamber during charge ignition of cases 1, 3 and 5

    图  8  各方案P1、P2、P3测点的压力曲线

    Figure  8.  Pressure-time curves measured at measuring points P1, P2 and P3 in each case

    表  1  各方案装填参数

    Table  1.   Charging parameters in each case

    方案电底火点火药量/g仿真发射药装药结构仿真发射药量/kg仿真发射药装填密度/(g·cm3)
    13#20单一杆状发射药2.70.58
    23#30单一杆状发射药2.70.58
    33#20单一粒状发射药3.00.64
    43#30单一粒状发射药3.00.64
    53#20杆状-粒状混装发射药1.8(杆状),0.9(粒状)0.58
    下载: 导出CSV

    表  2  各方案传火时间

    Table  2.   Flame speeding time of each case

    方案tL/mstR/mstt/ms
    10.81.62.4
    20.61.42.0
    31.02.63.6
    41.02.23.2
    50.82.33.1
    下载: 导出CSV

    表  3  方案1、3、5火焰传播位置函数拟合系数

    Table  3.   Exponential decay function coefficients in cases 1, 3 and 5

    方案火焰传播方向A/mmt1/msX0/mm
    1 960.73748 5.48834−746.77192
    −2140.6712610.281582129.96097
    3 337.19979 0.96082 −65.56418
    −670.86944 3.58306 675.19170
    5 449.90302 1.92159−216.47406
    −567.63017 1.26204 458.43930
    下载: 导出CSV

    表  4  各方案的最高压力和破膜压力

    Table  4.   The highest pressure and membrane-broken pressure in each case

    方案最高压力/MPa破膜压力/MPa
    12.510.98
    22.820.96
    32.530.85
    42.911.00
    52.631.10
    下载: 导出CSV
  • [1] 金志明, 翁春生. 高等内弹道学 [M]. 北京: 高等教育出版社, 2003: 237−239.
    [2] 翁春生, 王浩. 计算内弹道学 [M]. 北京: 国防工业出版社, 2006: 95−97.
    [3] 王升晨, 周彦煌, 刘千里, 等. 膛内多相燃烧理论及应用 [M]. 北京: 兵器工业出版社, 1994: 84−86.
    [4] ACHARYA R. Effect of different flash-tube vent-hole patterns on interior ballistic processes of ignition cartridge of 120-mm mortar system [J]. International Journal of Energetic Materials and Chemical Propulsion, 2008, 7(5): 383–397. DOI: 10.1615/IntJEnergeticMaterialsChemProp.v7.i5.30.
    [5] HOUIM R W, KUO K K. Understanding interior ballistic processes in a flash tube [J]. Journal of Applied Mechanics, 2010, 77(5): 051403. DOI: 10.1115/1.4001285.
    [6] MOORE J D, ACHARYA R, FERRARA P J. Effect of flash-tube vent-hole patterns on the combustion product discharge rate [J]. International Journal of Energetic Materials and Chemical Propulsion, 2009, 8(3): 199–220. DOI: 10.1615/IntJEnergeticMaterialsChemProp.v8.i3.30.
    [7] JIANG S P, RUI X T, HONG J, et al. Numerical simulation of impact breakage of gun propellant charge [J]. Granular Matter, 2011, 13(5): 611. DOI: 10.1007/s10035-011-0276-1.
    [8] 张瑞华, 芮筱亭, 赵宏立, 等. 基于离散单元法的发射装药挤压破碎模拟实验 [J]. 爆炸与冲击, 2021, 41(6): 062301. DOI: 10.11883/bzycj-2020-0157.

    ZHANG R H, RUI X T, ZHAO H L, et al. Simulational experiment on compression and fracture of propellant charge based on the discrete element method [J]. Explosion and Shock Waves, 2021, 41(6): 062301. DOI: 10.11883/bzycj-2020-0157.
    [9] 陈言坤, 罗兴柏, 甄建伟, 等. 发射药动态力学性能的研究进展 [J]. 含能材料, 2013, 21(5): 675–680. DOI: 10.3969/j.issn.1006-9941.2013.05.023.

    CHEN Y K, LUO X B, ZHEN J W, et al. Review on dynamical mechanical properties of propellants [J]. Chinese Journal of Energetic Materials, 2013, 21(5): 675–680. DOI: 10.3969/j.issn.1006-9941.2013.05.023.
    [10] 王浩, 梁世超, 张莺, 等. 火焰在传火管装药床中的传输特性研究 [J]. 爆炸与冲击, 1999, 19(1): 66–71.

    WANG H, LIANG S C, ZHANG Y, et al. Study of flame propagation characteristics in fire transfering tube [J]. Explosion and Shock Waves, 1999, 19(1): 66–71.
    [11] 王珊珊, 张玉成, 王浩, 等. 大长径比点火管高密实火药床点传火过程两相流的数值模拟 [J]. 爆炸与冲击, 2013, 33(4): 444–448. DOI: 10.11883/1001-1455(2013)04-0444-05.

    WANG S S, ZHANG Y C, WANG H, et al. Two-phase flow in ignition process of consolidated charge bed within a large length-to-diameter ratio igniter tube [J]. Explosion and Shock Waves, 2013, 33(4): 444–448. DOI: 10.11883/1001-1455(2013)04-0444-05.
    [12] 王珊珊, 王浩, 黄明, 等. 药床透气性对大长径比点传火管点传火性能的影响研究 [J]. 弹道学报, 2013, 25(4): 95–99. DOI: 10.3969/j.issn.1004-499X.2013.04.019.

    WANG S S, WANG H, HUANG M, et al. Research on influence of charge permeability on ignition characteristics on large length-diameter ratio igniter tube [J]. Journal of Ballistics, 2013, 25(4): 95–99. DOI: 10.3969/j.issn.1004-499X.2013.04.019.
    [13] 韩博, 张晓志, 邢浴仁, 等. 大口径火炮发射装药点传火模拟试验装置的研究 [J]. 兵工学报, 2008, 29(3): 262–265. DOI: 10.3321/j.issn:1000-1093.2008.03.002.

    HAN B, ZHANG X Z, XING Y R, et al. Research on simulation testing device of the ignition system for large caliber propellant charging [J]. Acta Armamentarii, 2008, 29(3): 262–265. DOI: 10.3321/j.issn:1000-1093.2008.03.002.
    [14] 程诚, 张小兵. 某制导炮弹二维两相流内弹道性能分析与数值模拟研究 [J]. 兵工学报, 2015, 36(1): 58–63. DOI: 10.3969/j.issn.1000-1093.2015.01.009.

    CHENG C, ZHANG X B. Two-dimensional numerical simulation on two-phase flow interior ballistic performance of a guided projectile [J]. Acta Armamentarii, 2015, 36(1): 58–63. DOI: 10.3969/j.issn.1000-1093.2015.01.009.
    [15] MIURA H, MATSUO A, NAKAMURA Y. Three-dimensional simulation of pressure fluctuation in a granular solid propellant chamber within an ignition stage [J]. Propellants, Explosives, Pyrotechnics, 2011, 36(3): 259–267. DOI: 10.1002/prep.201000058.
    [16] JANG J S, Oh S H, ROH T S. Development of three-dimensional numerical model for combustion-flow in interior ballistics [J]. Journal of Mechanical Science and Technology, 2016, 30(4): 1631–1637. DOI: 10.1007/s12206-016-0319-y.
    [17] JENARO G, OTON-MARTINEZ R A, MORATILLA D, et al. Analysis of pressure waves generation in gun combustion chambers during the propellant ignition process [J]. Propellants, Explosives, Pyrotechnics, 2016, 41(2): 228–237. DOI: 10.1002/prep.201500173.
    [18] JARAMAZ S, MICKOVIĆ D, ELEK P. Two-phase flows in gun barrel: theoretical and experimental studies [J]. International Journal of Multiphase Flow, 2011, 37(5): 475–487. DOI: 10.1016/j.ijmultiphaseflow.2011.01.003.
    [19] COLLETTI A, BOYER E, KUO K K, et al. Characterization of black-powder igniter jets into a granular bed through analysis of condensed-phase product deposition and regions of penetration [J]. International Journal of Energetic Materials and Chemical Propulsion, 2010, 9(6): 505–522. DOI: 10.1615/IntJEnergeticMaterialsChemProp.2011001422.
  • 期刊类型引用(9)

    1. 刘志宇,王昌建,汪泉,马鸿盛,张佳庆,黄玉彪. 变压器油爆炸作用下的消防设备抗爆实验研究. 实验力学. 2023(06): 762-770 . 百度学术
    2. 葛超,余庆波,卢冠成,郑元枫,王海福. 活性芯体子弹对柴油油箱引燃效应及机理研究. 北京理工大学学报. 2020(10): 1072-1080+1087 . 百度学术
    3. 黄勇,张红伟,申荣艳,解立峰,陈群. 抑爆柴油池火燃烧特性试验研究. 中国安全科学学报. 2019(08): 55-60 . 百度学术
    4. 宋刚,张文铖,封晓杰,李斌,黄勇,解立峰. 油箱内部强起爆条件下安全柴油燃爆性能的试验研究. 爆破器材. 2018(02): 28-32 . 百度学术
    5. 刘健,姚箭,宋述忠,李斌,解立峰,王永旭. 柴油的烤燃燃爆性能实验. 爆炸与冲击. 2018(03): 534-540 . 本站查看
    6. 李俊,鲁长波,安高军,熊春华,解立峰. 抑爆高闪点喷气燃料的抑爆特性. 高压物理学报. 2017(03): 328-334 . 百度学术
    7. 刘健,李斌,解立峰,姚箭,王永旭. 油箱明火烤燃下燃爆危险性分析与安全对策. 安全与环境工程. 2017(03): 144-149 . 百度学术
    8. 刘健,李斌,解立峰,姚箭,王永旭. 微乳化柴油燃爆性能研究. 中国安全生产科学技术. 2017(05): 140-144 . 百度学术
    9. 李俊,鲁长波,安高军,熊春华,解立峰,苗建敏. 含水柴油国内外研究进展. 化工进展. 2016(S1): 109-115 . 百度学术

    其他类型引用(4)

  • 加载中
图(8) / 表(4)
计量
  • 文章访问数:  443
  • HTML全文浏览量:  193
  • PDF下载量:  47
  • 被引次数: 13
出版历程
  • 收稿日期:  2021-01-21
  • 修回日期:  2021-03-19
  • 网络出版日期:  2021-09-27
  • 刊出日期:  2021-11-23

目录

/

返回文章
返回