二级高压驱动阵列弹珠同步弹射微型爆源的研制

郭纬; 徐小辉; 李干; 李杰; 蒋海明; 李志浩

doi:10.11883/bzycj-2021-0343

二级高压驱动阵列弹珠同步弹射微型爆源的研制

doi: 10.11883/bzycj-2021-0343

陆军工程大学爆炸冲击防灾减灾国家重点实验室，江苏南京 210007

基金项目: 国家自然科学基金（12072371）

详细信息

作者简介:
郭　纬（1998－　），男，硕士研究生, guo15965317797@163.com

通讯作者:
徐小辉（1983－　），男，博士，副教授, xuxiaohui168@126.com

中图分类号: O383; O389
计量
- 文章访问数: 289
- HTML全文浏览量: 150
- PDF下载量: 42
- 被引次数: 1
出版历程
- 收稿日期: 2021-08-16
- 修回日期: 2021-10-22
- 网络出版日期: 2022-08-01
- 刊出日期: 2022-09-09

Development of a miniature explosion device initiated by a synchronous launcher of marbles driven by two-stage high-pressure gas

State Key Laboratory of Explosion & Impact and Disaster Prevention & Mitigation, Army Engineering University of PLA, Nanjing 210007, Jiangsu, China

摘要

摘要: 针对当前大当量地下爆炸真空室模拟试验中爆源起爆方式高度依赖火药制品等问题，基于地下爆炸相似理论和二级气炮原理，自主研制了二级高压驱动阵列弹珠同步弹射微型爆源装置。装置利用二级高压气体驱动弹珠同步击碎玻璃球壳，释放球内高压气体，以模拟真实爆炸气体生成物的推出。整套爆源装置的发射参数：高压气室充气压力4 MPa，玻璃球壳内残余稳态气体压力约为3 kPa，能够用于0~20 kt当量地下爆炸成坑效应的真空室模拟。爆源适用性试验验证表明，该爆源装置的爆破机制和爆破效果满足大当量地下抛掷爆炸真空室模拟试验的功能需求，且具有较高的安全性、可控性和可操作性，为开展相关模拟试验提供了新的技术方法。
- 大当量地下爆炸 /
- 真空室模拟试验 /
- 模型爆源 /
- 二级高压驱动
Abstract: Aiming at the problem that the initiation mode of the explosion device is highly dependent on the gunpowder products in the simulation experiments of large-scale underground explosions in a vacuum chamber, and based on the similarity theory of underground explosions and the principle of the two-stage gas gun, a micro explosion device initiated by a synchronous launcher of marbles driven by two-stage high-pressure gas was developed independently. A glass enclosure with compressed gas (filled by air compressor) was used to simulate the high-pressure cavity generated at the beginning of a real underground explosion. Two-stage high-pressure gas was used to drive marbles to break the glass shell synchronously, thus releasing the high-pressure gas in the spherical shell to simulate the ejection of gas products in a real underground explosion. The pressure in the launcher chamber is 4 MPa, and the residual steady-state gas pressure in the glass enclosure is about 3 kPa. The above set of the launch parameters can be used for simulation experiments of real underground explosions with an equivalent of 0−20 kt TNT. Through high-speed imaging of the air and water blasting sphericity tests, the reliability of the explosion device and the sphericity of the blasting effect were verified. When there is a difference in the internal and external pressure of the glass spherical shell, the cracks of the shell are fully developed and the fragments are evenly distributed. The applicability test shows that the blasting mechanism and blasting effect of the explosion device can meet the requirements of the simulation experiment of large-scale underground explosions in the vacuum chamber, and the device has the characteristics of high efficiency, low pollution, convenient operation, good repeatability, good controllability and low requirements for site conditions, which can provide a novel technology for the simulation experiments of large-scale underground explosions in the vacuum chamber.
- large-scale underground explosion /
- similarity experiment in a vacuum chamber /
- explosion device model /
- two-stage gas gun

HTML全文

多胞材料在高速冲击下具有变形局部化的典型特征。Reid等^[1-2]研究了由圆环组成的一维系统的动态冲击行为，提出了“结构冲击波”的概念，并在对木材的动态压溃实验中证实了压溃带的冲击波传播特征。Zheng等^[3]、Liu等^[4]在对随机蜂窝的动态冲击中观察到3种变形模式，即均匀模式、过渡模式和冲击模式。张新春等^[5]、胡玲玲等^[6]探讨了胞元构型对蜂窝结构动态冲击性能的影响。Zou等^[7]认为蜂窝中传播的冲击波波阵面在高速冲击下具有胞元尺寸的宽度。Liao等^[8]采用了局部应变场计算方法表征了蜂窝铝中一维冲击波的传播并得到了冲击波速度的大小。一系列的一维冲击波模型被发展来表征多胞材料中塑性压溃行为，如R-PP-L模型(率无关、刚性-理想塑性-锁定假设)^{[2, 9-10]}和R-PH模型(率无关，刚性-塑性硬化假设)^[11-13]。然而，大部分的冲击波模型都是基于准静态名义应力-应变曲线得到的，直接得到动态的应力应变关系在实验上是困难的。

多胞材料在动态冲击下存在变形局部化和应力增强等现象，名义应力应变曲线失去了物理意义^[4]。Zheng等^[13]通过对泡沫铝的三维细观有限元模型的分析得到了波后应力应变状态点，并提出了D-R-PH(动态，刚性-塑性硬化)模型来表征动态应力-应变行为。Barnes等^[14]、Sun等^[15]分别对开孔泡沫和闭孔泡沫进行实验和数值研究，也得到了多胞材料的动态应力应变状态。但这些研究都是针对冲击模式下的应力应变状态，对过渡模式下的应力应变状态的认识仍然不清楚。此外，文献[13-15]中波后应力都是由撞击端的名义应力得到的，缺少对材料内部局部应力信息的直接表征。

认识多胞材料在中速过渡模式下的应力应变行为，对于完整地理解多胞材料的动态本构关系有重要意义。本文中将采用截面应力计算方法研究基于细观有限元模型的随机蜂窝结构在恒速冲击下的应力分布，分析冲击波在随机蜂窝中的传播特性，表征随机蜂窝的动态应力应变关系。

1. 数值模型

1.1 细观有限元模型

采用二维Voronoi技术构造随机蜂窝结构，详细建模过程参见文献[3]。如图 1，蜂窝试件的长和宽均为50 mm，试件的不规则度为0.3，由500个胞元组成。试件相对密度ρ₀/ρ_s= 0.1，其中ρ₀为蜂窝试样变形前的初始密度，ρ_s为基体材料的密度。蜂窝试件的平均胞元尺寸约为2.5 mm，其值定义为与平均胞元面积相等的圆的直径。

图 1 Voronoi蜂窝试件细观有限元模型

Figure 1. A cell-based finite element model of a Voronoi honeycomb specimen

下载: 全尺寸图片幻灯片

采用ABAQUS/Explicit有限元软件对蜂窝材料的恒速压溃过程进行数值模拟。胞壁材料由S4R壳单元(4节点，减缩积分、沙漏控制、有限膜应变)进行模拟。通过网格收敛性分析将模型划分了10 480个壳单元，平均单元长度约为0.2 mm。基体材料属性设置为弹性-理想塑性，材料参数设置为泊松比ν = 0.3，弹性模量E= 66 GPa，屈服应力σ_y = 175 MPa，基体密度ρ_s = 2 700 kg/m³。所有接触面均定义为通用接触，摩擦因数取为0.02。为了模拟面内变形，模型中所有节点的面外变形均被限制。

蜂窝试件的左右两端设置两个刚性面，左侧刚性面以恒定不变的速度v沿x方向向右压缩试件。在压缩过程中右侧刚性板为固支端，如图 1所示。

1.2 截面应力计算方法

对于蜂窝试件在承受匀速压缩的情况下，一维近似是适用的，因此可以通过截面工程应力来表征局部应力信息，详见文献[16]。不同的拉格朗日截面位置的工程应力由截面内力除以截面初始面积得到。截面的内力由两部分组成，一部分是通过基体材料的单元节点所传递的力，称为节点力，记为F_x⁽ⁿ⁾；另一部分是由变形中胞壁接触所引起的接触力，记为F_x^(c)。这两部分力由下式可以计算得到:

$\left\{ \begin{array}{l} F_x^{\left( {\rm{n}} \right)} = \sum\limits_{i = 1}^M {F_{xi}^{\left( {\rm{n}} \right)}} \\ F_x^{\left( {\rm{c}} \right)} = \sum\limits_{i = 1}^K {\left( {F_{xi}^{\left( {{\rm{cn}}} \right)} + F_{xi}^{\left( {{\rm{ct}}} \right)}} \right)} \end{array} \right.$

(1)

式中:M为所考察的截面所贯穿的单元在参考构型中左侧的节点数，K是从撞击端到所考察截面的总的节点数目。对于恒速压缩情形，K包含了刚性端的节点。接触力F_x^(c)的简化计算公式的算法得益于牛顿第三定律，完整的接触对的接触力大小相等方向相反，互相抵消之后剩余的部分为被截面所隔开的接触对的接触力。由于刚性面与试件撞击端位置重合，因此撞击端面上的接触力由刚性面输出接触力得到。F_xi⁽ⁿ⁾、F_xi^(cn)、F_xi^(ct)分别为任意节点编号i的节点力的x分量以及节点的法向和切向接触力的x分量，均可以直接从ABAQUS的输出文件提取得到。截面的工程应力可以由下式计算为：

$\sigma = \sigma _x^{\left( {\rm{n}} \right)} + \sigma _x^{\left( {\rm{c}} \right)}$

(2)

式中:σ_x⁽ⁿ⁾=F_x⁽ⁿ⁾/A₀，表示由节点所传递的应力；σ_x^(c)=F_x^(c)/A₀，表示由胞壁接触所引发的应力。其中A₀为截面的初始面积。

2. 结果

2.1 一维应力分布

冲击速度分别为200和50 m/s时的应力分布如图 2所示。在任意名义应变时刻下，应力会在某个位置发生陡然减小，并且这个位置随着名义应变的增大往试件右侧移动。这说明在高速和中速撞击下试件内均存在冲击波的传播，应力发生急剧变化的位置即为冲击波波阵面的位置，中速撞击下的波阵面区域相较于高速撞击更加平缓。从图 2可以明显地观察到应力增强现象，跨过冲击波阵面，应力由波前应力迅速增加到波后应力并几乎保持不变。

图 2 两种冲击速度下的应力分布图

Figure 2. One-dimensional stress distributions at two impact velocities

下载: 全尺寸图片幻灯片

2.2 基于应力分布得到冲击波速度

可以由应力分布梯度来确定冲击波波阵面Φ的位置，如图 3所示。应力分布梯度反映了一维应力在加载方向上的变化，在应力梯度绝对值达到最大的拉格朗日坐标位置，应力变化最快，这个位置即为当前时刻冲击波波阵面的位置。中速和高速压溃下波阵面位置均可以采用此方法得到。不同冲击速度下得到的波阵面位置随冲击时间变化关系如图 4所示。在恒速压缩情形下，冲击波的波阵面位置与冲击时间关系为线性关系。可通过对线性关系拟合得到的斜率来预估冲击波速度。从图 4中可以看出，冲击速度越大，斜率越大，则相应的冲击波速度也越大。

图 3 一维的应力分布及其应力梯度分布

Figure 3. One-dimensional stress distributions and the corresponding stress gradients

下载: 全尺寸图片幻灯片

图 4 不同冲击速度下冲击波位置随时间的关系

Figure 4. Variation of shock front position with impact time at different impact velocities

下载: 全尺寸图片幻灯片

2.3 一维冲击波理论

除了由有限元结果直接得到的应力信息来得到冲击波速度之外，一维冲击波理论也可以对蜂窝中的冲击波传播行为进行描述。由应力波理论^[17]，跨波阵面的质量和动量守恒关系分别为:

${v_{\rm{B}}}\left( t \right) - {v_{\rm{A}}}\left( t \right) = {v_{\rm{s}}}\left( t \right)\left( {{\varepsilon _{\rm{B}}}\left( t \right) - {\varepsilon _{\rm{A}}}\left( t \right)} \right)$

(3)

${\sigma _{\rm{B}}}\left( t \right) - {\sigma _{\rm{A}}}\left( t \right) = {\rho _0}{v_{\rm{s}}}\left( t \right)\left( {{v_{\rm{B}}}\left( t \right) - {v_{\rm{A}}}\left( t \right)} \right)$

(4)

式中:v_s(t)是冲击波速度, σ_B、ε_B、v_B为波后方的应力、应变和速度，σ_A、ε_A、v_A为波前方的应力、应变和速度。联立这两个守恒关系可以得到冲击波速度与波前波后应力、应变的关系为：

${v_{\rm{s}}}\left( t \right) = \sqrt {\frac{{{\sigma _{\rm{B}}}\left( t \right) - {\sigma _{\rm{A}}}\left( t \right)}}{{{\rho _0}\left( {{\varepsilon _{\rm{B}}}\left( t \right) - {\varepsilon _{\rm{A}}}\left( t \right)} \right)}}}$

(5)

如果能直接从有限元结果中得到波前和波后应力/应变信息，就可以由式(5)计算得到冲击波速度的大小。为了简便，这里将这种结合冲击波理论和数值模拟结果的方法叫做一维冲击波理论方法。

通过对一维应力和应变分布求平均值，可以得到波前和波后的应力和应变，如图 5所示。图中一维应变分布采用文献[8]提出的局部应变计算方法得到。在求平均的过程中，为了消除冲击波波阵面平均效应的影响，将紧邻波阵面附近1.5倍胞元半径R_c内的数据排除掉。波后应力和应变信息由冲击端到Φ-1.5R_c区域内平均得到，波前应力应变信息由Φ+1.5R_c到支撑端区域平均得到。再由式(5)就可以得到不同冲击速度下的冲击波速度。

图 5 一维应力和应变分布及其理想化

Figure 5. One-dimensional stress and strain distributions and the idealizations

下载: 全尺寸图片幻灯片

2.4 冲击波速度的比较

通过截面应力和局部应变场等有限元方法以及一维冲击波理论方法都可以得到冲击波速度。此外，也可由已有的冲击波模型得到冲击波速度的大小。由R-PP-L模型^[2]得到的冲击波速度为:

${v_{\rm{s}}} = v/{\varepsilon _{\rm{L}}}$

(6)

式中:ε_L为多胞材料的锁定应变。由R-PH模型^[13]得到的冲击波速度为:

${v_{\rm{s}}} = v + \sqrt {C/{\rho _0}}$

(7)

式中:C为多胞材料的塑性硬化参数。在本文中，ε_L=0.64，C=0.226 MPa。

由不同方法得到的冲击波速度如图 6所示。由一维冲击波理论方法以及R-PH冲击波模型得到的冲击波速度比较接近有限元方法得到的冲击波速度，并且与冲击速度的差值均为一个常值，即材料冲击参数。然而由R-PP-L模型预估得到的冲击波速度明显高于有限元结果。这是由于R-PP-L模型采用的是锁定不变的压实应变，实际上其值小于动态压溃的压实应变。由质量守恒关系可以发现，压实应变越小，冲击波速度越大。R-PH模型中的压实应变是随冲击速度变化的动态应变，因此更加适合描述多胞材料在动态压溃下的冲击波传播。

图 6 不同方法得到的冲击波速度的比较

Figure 6. Comparison of shock wave speeds obtained using different methods

下载: 全尺寸图片幻灯片

3. 分析和讨论

3.1 冲击波速度的一致近似

前文中通过应力分布得到了高速冲击下的冲击波速度与冲击速度的关系。在中等冲击速度下，仍然可以观察到冲击波传播的现象^{[4, 8]}，因此仍然采用应力梯度分布方法来确定冲击波速度的大小，结果见图 7。当冲击速度很高的时候，冲击波速度与冲击速度呈线性关系并且两者之间的差值几乎为一个常数^[13]。但随着冲击速度的减小，冲击波速度趋近于常数。因此，冲击波速度对冲击速度的一阶导数在低速和高速下分别趋于0和1，从低速到高速的完整曲线具有明显的S型特征。利用S函数(sigmoid function)，冲击波速度对冲击速度的一阶导数可以写作：

$\frac{{{\rm{d}}{v_{\rm{s}}}}}{{{\rm{d}}v}} = \frac{1}{{1 + {{\rm{e}}^{\frac{{b - v}}{a}}}}}$

(8)

图 7 冲击波速度与冲击速度的关系

Figure 7. Variations of the shock wave speed with the impact velocity

下载: 全尺寸图片幻灯片

式中:a和b为待定参数。上式的积分形式可以写作：

${v_{\rm{s}}} = a\ln \left( {1 + {{\rm{e}}^{\frac{{v - b}}{a}}}} \right) + c$

(9)

式中:c为待定参数。可用式(9)对有限元结果进行拟合，从图 7可见拟合效果较好，其中拟合参数为a=3.603 m/s, b=40.68 m/s和c=54.91 m/s。若式(8)右边的值较小，对应于低速冲击，此时并不能观察到冲击波传播的现象，因此不妨取该值为0.05，可以估算从均匀模式向过渡模式转变的临界速度v_c1=b－aln19=30.07 m/s。若式(8)右边的值较大，对应于高速冲击，不妨取该值为0.95，可以估算从过渡模式向冲击模式转变的临界速度v_c2=b+aln19=51.29 m/s。可见参数a和b的线性组合反映了模式转变速度的大小。当0＜v≪b时，冲击波速度v_s近似等于c。当v≫b时， $1 + {{\rm{e}}^{\frac{{v - b}}{a}}} \approx {{\rm{e}}^{\frac{{v - b}}{a}}}$ ，此时冲击波速度v_s=v-b+c，和R-PH模型中得到的冲击波速度与冲击速度变化关系一致，即c₁=c-b。因此，式(9)表征了冲击波速度与冲击速度的一致近似关系。

3.2 动态应力应变状态的一致近似

在本文中考察的恒速冲击情形中，若取ε_A=0, v_A=0, v_B=v和σ_A=σ₀，由跨波阵面的守恒关系式(3)和(4)有:

$\left\{ \begin{array}{l} {\varepsilon _{\rm{B}}} = v/{v_{\rm{s}}}\\ {\sigma _{\rm{B}}} = {\sigma _0} + {\rho _0}v{v_{\rm{s}}} \end{array} \right.$

(10)

将式(9)代入上式，可得：

$\left\{ \begin{array}{l} {\varepsilon _{\rm{B}}} = v/\left[ {a\ln \left( {1 + {{\rm{e}}^{\frac{{v - b}}{a}}}} \right) + c} \right]\\ {\sigma _{\rm{B}}} = {\sigma _0} + {\rho _0}v\left[ {a\ln \left( {1 + {{\rm{e}}^{\frac{{v - b}}{a}}}} \right) + c} \right] \end{array} \right.$

(11)

与有限元结果的比较如图 8所示，图中虚线对应于式(11)。可见式(11)与有限元结果存在一定的误差。事实上，在恒速冲击下，冲击波波阵面前方的应变和粒子速度不能直接忽略。

图 9给出了不同冲击速度下试件中的速度分布情况，局部粒子速度在冲击波波阵面位置附近约两个胞元区间内急剧变小。当冲击速度为中等大小时，波阵面位置处的粒子速度变化更为平缓。整个波前区域处于低速状态，近似呈线性分布，与冲击速度的大小关系不大。冲击波前方的粒子速度和变形不能忽略，可以通过图 9估算v_A≈10 m/s，初始压溃应变ε_A取准静态实验中常用的ε₀ =0.02。因此，由跨波阵面的守恒关系式(3)和(4)有:

$\left\{ \begin{array}{l} {\varepsilon _{\rm{B}}} = {\varepsilon _{\rm{0}}} + \left( {v - {v_{\rm{A}}}} \right)/{v_{\rm{s}}}\\ {\sigma _{\rm{B}}} = {\sigma _0} + {\rho _0}\left( {v - {v_{\rm{A}}}} \right)/{v_{\rm{s}}} \end{array} \right.$

(12)

图 8 波后应变和波后应力随冲击速度的变化

Figure 8. Variations of strain and stress behind shock front with impact velocity

下载: 全尺寸图片幻灯片

图 9 不同冲击速度下的速度分布

Figure 9. Velocity distributions at different impact velocities

下载: 全尺寸图片幻灯片

将式(9)代入上式，可得:

$\left\{ \begin{array}{l} {\varepsilon _{\rm{B}}} = {\varepsilon _{\rm{0}}} + \left( {v - {v_{\rm{A}}}} \right)/\left[ {a\ln \left( {1 + {{\rm{e}}^{\frac{{v - b}}{a}}}} \right) + c} \right]\\ {\sigma _{\rm{B}}} = {\sigma _0} + {\rho _0}\left( {v - {v_{\rm{A}}}} \right)/\left[ {a\ln \left( {1 + {{\rm{e}}^{\frac{{v - b}}{a}}}} \right) + c} \right] \end{array} \right.$

(13)

如所示，这个结果与有限元结果吻合较好。引入无量纲参数 $\xi = {{\rm{e}}^{\frac{{v - b}}{a}}}$ ，式(13)可以进一步写作：

$\left\{ \begin{array}{l} {\varepsilon _{\rm{B}}} = {\varepsilon _{\rm{0}}} + \frac{{\ln \xi + \beta - \alpha }}{{\ln \left( {1 + \xi } \right) + \gamma }}\\ {\sigma _{\rm{B}}} = {\sigma _0} + B\left( {\ln \xi + \beta - \alpha } \right)\left[ {\ln \left( {1 + \xi } \right) + \gamma } \right] \end{array} \right.$

(14)

式中:α=v_A/a, β =b/a, γ=c/a和B=ρ₀a²。因此，式(14)给出恒速冲击下动态应力应变状态的一致近似关系。准静态和动态应力应变关系的比较如图 10所示。在图 10中，准静态应力应变曲线及动态应力应变状态点由有限元计算得到，动态应力应变曲线由式(14)得到。在动态应力应变关系的曲线上的每个点都对应着特定的一个冲击速度。在高应变情形下，动态应力应变关系位于准静态曲线右侧，这与已有文献中的结论一致。随着应变的减小，动态的应力应变关系逼近准静态应力应变曲线。

图 10 准静态及动态应力应变关系

Figure 10. Quasi-static and dynamic stress-strain relations

下载: 全尺寸图片幻灯片

4. 结论

采用细观有限元模型研究了随机蜂窝结构在恒速冲击下的动态性能，采用截面应力计算方法得到了随机蜂窝试件内的应力分布和冲击波速度，比较了高速冲击下由不同方法得到的冲击波速度与冲击速度的关系。结果表明，R-PP-L模型高估了冲击波速度，但R-PH模型以及一维冲击波理论得到的冲击波速度与有限元模拟得到的结果接近。在冲击速度很高的时候，冲击波速度与冲击速度的关系趋于线性，但随着冲击速度的减小，冲击波速度不断减少并趋于常数。最后，发展了可以表征冲击波速度与冲击速度的关系的一致近似模型，并基于一维冲击波理论发展了动态应力应变关系的一致近似模型。这个模型可以对多胞材料在过渡模式及冲击模式下的动态压溃行为进行表征。

图 1 微型爆源装置

Figure 1. Miniature explosion device

下载: 全尺寸图片幻灯片

图 2 微型爆源装置设计原理

Figure 2. Schematic of the micro explosion device

下载: 全尺寸图片幻灯片

图 3 二级高压驱动阵列弹珠同步弹射装置设计原理

Figure 3. Principle of marble synchronous launcher driven by two-stage high-pressure gas

下载: 全尺寸图片幻灯片

图 4 炮管顶部发射弹头设计及实物

Figure 4. Schematic diagram and photo of the barrel end

下载: 全尺寸图片幻灯片

图 5 活塞及破膜

Figure 5. Piston and ruptured diaphragms

下载: 全尺寸图片幻灯片

图 6 玻璃球内稳态压力测试

Figure 6. Tests of the steady-state pressure in the glass shell

下载: 全尺寸图片幻灯片

图 7 空气中爆破球形度测试

Figure 7. Blasting sphericity tests in the air

下载: 全尺寸图片幻灯片

图 8 空气中玻璃球壳（内外无压差）爆破高速分幅摄影图像（试验S1）

Figure 8. High-speed split photos of a glass shell (without pressure difference inside and outside) blasting in air (test S1)

下载: 全尺寸图片幻灯片

图 9 空气中玻璃球壳（内外压差80 kPa）爆破高速分幅摄影图像（试验S2）

Figure 9. High-speed split photos of a glass shell (with an internal and external pressure difference of 80 kPa) blasting in air (test S2)

下载: 全尺寸图片幻灯片

图 10 玻璃球壳水中爆炸高速分幅镜头（试验S3）

Figure 10. High-speed split photos of a glass shell blasting in water (test S3)

下载: 全尺寸图片幻灯片

图 11 水中玻璃球壳爆炸与水面相互作用

Figure 11. Interaction of the glass shell in water with the water surface

下载: 全尺寸图片幻灯片

表 1 地下爆炸空腔气体势能计算表达式^[10]

Table 1. Formulas for the potential energy of the cavity gas in underground explosions^[10]

岩石特性	空腔气体势能
不含气体岩石	$A = {{0.49q}/ {\bar r_{\text{n}}^{0.84}}}$
仅含自由水的硅酸盐类岩石（适用于花岗岩、凝灰岩、冲积层等岩石）	$A = \dfrac{{0.49q}}{{\bar r_{\text{n}}^{0.84}}}\left( {1 + 5.8\eta _{\text{w}}^{0.7}} \right)$
仅含碳酸气的碳酸盐类岩石（适用于硬石膏、方解石、石灰岩等岩石）	$A = \dfrac{{0.49q}}{{\bar r_{\text{n}}^{0.84}}}\left( {1 + 1.96\eta _{{\text{c}}{{\text{o}}_{\text{2}}}}^{0.7}} \right)$
混合含气岩石	$A = \dfrac{{0.49q}}{{\bar r_{\text{n}}^{0.84}}}\left( {1 + 5.8\eta _{{\varepsilon }}^{0.7}} \right)$

下载: 导出CSV

表 2 弹珠弹射装置适用性测试试验部分结果

Table 2. Partial test results of the applicability of the marble launcher

试验序号	高压气室压力/MPa	玻璃球壳直径/cm	玻璃球壳稳态压力/kPa	试验情况
1-1	1	10	0.7	活塞未发射
1-2	1	10	1.1	活塞未发射
2-1	2	10	3.4	活塞正常发射，膜片破裂，活塞与炮管锥段齐平，入锥不充分
2-2			3.6
2-3			1.8
4-1	4	10	2.9	活塞正常发射，膜片破裂，入锥充分，玻璃球壳稳态压力波动较小
4-2			3.0
4-3			2.7
5-1	5	10	3.6	活塞正常发射，膜片破裂，入锥充分，玻璃球壳稳态压力偏高。
5-2			7.7
5-3			6.9

下载: 导出CSV

表 3 不同规模大当量地下爆炸真空室模拟试验主要参数

Table 3. Key parameters for the vacuum chamber simulation tests of large-scale underground explosions

爆炸代号	等效TNT 当量/kt^[10]	真实爆炸空腔半径/m^[10]	模拟玻璃球壳半径/cm	玻璃球壳气压/kPa	局部二级高压气体稳态压力/kPa	驱动气体体积在球内占比/%
Neptun	0.115	7.30	5.00	135	3.00	2.20
1003竖井	1.10	13.3	5.00	87.0	3.00	3.40
Palanquin	4.30	17.9	5.00	84.3	3.00	3.60
125竖井	19.0	26.7	5.00	60.4	3.00	4.96
Schooner	31.0	47.1	5.00	15.9	3.00	18.9
Sedan	100	69.6	5.00	12.8	3.00	23.4

下载: 导出CSV

表 4 爆源球形度试验

Table 4. Sphericity test for the explosion device

试验	介质	球壳中绝对气压/kPa	玻璃球壳埋深/cm	拍摄频率/kHz	高压气室压力/MPa
S1	空气	100	−	5	4
S2	空气	180	−	5	4
S3	水	180	25	3	4

下载: 导出CSV

参考文献(14)

[1]	ADUSHKIN V V, KHRISTOFOROV B D. Craters of large-scale surface explosions [J]. Combustion, Explosion, and Shock Waves, 2004, 40(6): 674–678. DOI: 10.1023/B:CESW.0000048270.62239.01.
[2]	SADOVSKII M A, ADUSHKIN V V, RODIONOV V N, et al. A method of modeling large cratering explosions [J]. Combustion, Explosion, and Shock Waves, 1969, 3(1): 73–79. DOI: 10.1007/BF00741616.
[3]	ADUSHKIN V V, PERNIK L M. Design of linear charges for caving a slope by a large-scale explosion [J]. Soviet Mining, 1989, 25(6): 505–509. DOI: 10.1007/BF02528298.
[4]	ADUSHKIN V V, KAMALYAN R Z, KOROLEV K D. Character of the interaction of concentrated excavation charges detonated at different times [J]. Combustion, Explosion, and Shock Waves, 1989, 25(6): 782–785. DOI: 10.1007/BF00758751.
[5]	BLINOV I M, VAKHRAMEEV Y S. The method for modelling large-scale outburst explosions by microexplosions of explosive charges [J]. Fizika Goreniya I Vzryva, 1995, 31(2): 102–109.
[6]	VAKHRAMEEV Y S. Physical foundations for approximately modeling explosions with ejecta [J]. Combustion, Explosion, and Shock Waves, 1995, 31(1): 120–125. DOI: 10.1007/BF00755969.
[7]	徐小辉, 邱艳宇, 王明洋, 等. 大当量浅埋地下爆炸抛掷成坑效应的缩比模拟实验装置 [J]. 爆炸与冲击, 2018, 38(6): 1333–1343. DOI: 10.11883/bzycj-2017-0144. XU X H, QIU Y Y, WANG M Y, et al. Development of the testing apparatus for modeling large equivalent underground cratering explosions [J]. Explosion and Shock Waves, 2018, 38(6): 1333–1343. DOI: 10.11883/bzycj-2017-0144.
[8]	徐小辉, 邱艳宇, 王明洋, 等. 大当量地下浅埋爆炸真空室模拟相似材料研究 [J]. 岩石力学与工程学报, 2018, 37(S1): 3550–3556. DOI: 10.13722/j.cnki.jrme.2016.1539. XU X H, QIU Y Y, WANG M Y, et al. Similar materials for vacuum chamber model test under large scale throw blasting [J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(S1): 3550–3556. DOI: 10.13722/j.cnki.jrme.2016.1539.
[9]	王明洋, 徐小辉, 邱艳宇, 等. 一种用于模拟爆炸效应的爆源装置: CN201710295250.6 [P]. 2018-11-13.
[10]	ADUSHKIN V V, SPIVAK A. Underground explosions: WGC-2015-03 [R]. Lexington: Weston Geophysical Corp, 2015.
[11]	冯建宁, 彭炎午, 林俊德. 新型非火药驱动二级轻气炮内弹道诸问题的研究 [J]. 西北工业大学学报, 1994(3): 477–481. FENG J N, PENG Y W, LIN J D. On some problems of interior ballistics in a new type two-stage light gas gun with non-powder projection [J]. Journal of Northwestern Polytechnical University, 1994(3): 477–481.
[12]	TANG W Q, WANG Q, WEI B C, et al. Performance and modeling of a two-stage light gas gun driven by gaseous detonation [J]. Applied Sciences, 2020, 10(12): 4383. DOI: 10.3390/app10124383.
[13]	RINGROSE T J, DOYLE H W, FOSTER P S, et al. A hypervelocity impact facility optimised for the dynamic study of high pressure shock compression [J]. Procedia Engineering, 2017, 204: 344–351. DOI: 10.1016/j.proeng.2017.09.756.
[14]	徐坤博, 龚自正. 超高速发射技术研究进展 [C] // 中国数学力学物理学高新技术交叉研究学会第十三届学术年会. 敦煌: 中国数学力学物理学高新技术交叉研究学会, 2010.

施引文献

期刊类型引用(0)
其他类型引用(1)

资源附件(0)

访问统计

图(11) / 表(4)

计量

文章访问数: 289
HTML全文浏览量: 150
PDF下载量: 42
被引次数: 1

1. 数值模型
1.1 细观有限元模型
1.2 截面应力计算方法
2. 结果
2.1 一维应力分布
2.2 基于应力分布得到冲击波速度
2.3 一维冲击波理论
2.4 冲击波速度的比较
3. 分析和讨论
3.1 冲击波速度的一致近似
3.2 动态应力应变状态的一致近似
4. 结论

二级高压驱动阵列弹珠同步弹射微型爆源的研制

doi: 10.11883/bzycj-2021-0343

作者简介: 郭 纬（1998－ ），男，硕士研究生, guo15965317797@163.com

通讯作者: 徐小辉（1983－ ），男，博士，副教授, xuxiaohui168@126.com

计量

出版历程