间隙对爆轰加载下金属飞片运动特征影响的模拟分析

刘军; 孙致远; 张凤国; 殷建伟

doi:10.11883/bzycj-2022-0239

间隙对爆轰加载下金属飞片运动特征影响的模拟分析

doi: 10.11883/bzycj-2022-0239

北京应用物理与计算数学研究所，北京 100094

基金项目: 国家自然科学基金（12101062）；于敏基金（TCYM1820-02）

详细信息

作者简介:
刘　军（1981－　），男，硕士，副研究员，caepcfd@126.com

中图分类号: O383
计量
- 文章访问数: 557
- HTML全文浏览量: 177
- PDF下载量: 112
- 被引次数: 2
出版历程
- 收稿日期: 2022-05-31
- 修回日期: 2022-08-25
- 网络出版日期: 2022-09-09
- 刊出日期: 2023-04-05

Simulation analysis of the effect of clearance on motion characteristic of metal flyer under detonation loading

Institute of Applied Physics and Computational Mathematics, Beijing 100094, China

摘要

摘要: 在炸药爆轰驱动含间隙双层钢飞片情况下，两层钢飞片间的间隙会影响外层飞片的首次、二次入射波波形及强度，进而影响外层飞片自由面速度等。为了更好地认识爆轰加载条件下金属飞片的运动特征，需要深入研究间隙对该动力学过程的影响规律。首先，开展了爆轰驱动含初始间隙双层钢飞片的简化建模及数值模拟，通过模拟与实验结果的对比，验证了简化建模的合理性；然后，对该模型的加载动力学过程进行了深入分析，给出了首次、二次加载的来源；最后，开展了不同间隙厚度对该动力学过程影响的模拟分析研究。自由面速度结果表明，随着间隙厚度由0.1 mm增加至1 mm以上，外层飞片自由面的首次起跳速度峰值先逐渐降低后基本保持恒定、二次起跳速度峰值由逐渐增加至基本不变。动力学分析结果表明，可将不同间隙大小的影响分为两个阶段，其分界判据是在爆轰加载后内层金属飞片是否能够在间隙部位发展为脱体层裂片：在间隙较小的情况下，内层飞片在间隙一侧无法发展为层裂片，在此阶段内，随着间隙厚度的增加，外层飞片的首次加载峰值压力降低、二次加载峰值压力增加；在间隙较大的情况下，内层飞片在间隙一侧可以形成厚度不变、速度稳定的层裂片，在此阶段内，随着间隙厚度的增加，外层飞片的首次加载与二次加载的峰值压力均基本不变，但首次与二次加载之间的时间间隔缩短。研究结果对爆轰驱动含间隙飞片的自由面速度曲线的解读具有指导意义，从而能够更好地认识工程实验中由间隙造成的一些非预期物理现象。
- 爆轰驱动 /
- 金属飞片 /
- 间隙 /
- 二次冲击
Abstract: Clearance of certain thickness often exists between two stacked metal flyers. When a double-layer metal flyer with clearance is loaded by detonation, the closing of the clearance may affect the form and shock intensity of the first and second loading waves inside of the outer flyer, and then affects the free surface velocity of the outer flyer. In order to better grasp the motion characteristics under detonation loading, the effect of clearance on the dynamic process needs to be studied. Firstly, a detonation driven two-layer steel flyers model is presented, in which a clearance of certain thickness is assumed to exist between two steel flyers. In this model, the free surface of the outer flyer is loaded twice. By comparing the simulation results and experimental results of free surface velocity at different positions, it is confirmed that the simulation can correctly catch the dynamic process. Then, the sources of the first and second loading in the outer flyer are given by the analysis of the simulated dynamic process. The first loading wave in the outer flyer comes from the clearance closing collision, and the second loading wave mainly comes from the sustained high pressure loading of detonation products. Finally, the simulation with various clearance thicknesses is carried out, and the effect of clearance thickness change is summarized. The simulated results of free surface velocity show that with the increase of clearance thickness from 0.1 mm to more than 1 mm, the peak value of the first take-off free surface velocity first decreases and then remains unchanged, and the peak value of the second take-off free surface velocity first increases and then remains unchanged. The dynamic analysis shows that the size of the clearance thickness directly affects whether the inner steel flyer has enough time to develop into spallation on the clearance side after detonation loading. If the size of clearance is small, the inner flyer cannot develop into a spallation on clearance side, and the first loading wave formed in the outer flyer has a triangular like pulse. In this stage, with the increase of the clearance thickness, the first loading peak pressure decreases and the second loading peak pressure increases. If the size of clearance is large, the inner flyer can form a spallation with constant thickness and stable velocity on clearance side, and the first loading wave formed in the outer flyer is an approximate square wave. In this stage, with the increase of clearance thickness, the peak pressures of the first and second loading remain basically unchanged, but the time interval between the first and second loading decreases. The understanding has guiding significance for the interpretation of the free surface velocity measurement results in experiments, and some unexpected physical phenomena caused by clearance in practical problems could be better understood, too.
- detonation driving /
- metal flyer /
- clearance /
- second loading

HTML全文

在炸药爆轰驱动金属飞片时，由于装配精度^[1]、应力导致飞片变形^[2]和设计需求^[3]等原因，各个构件间不是紧密贴合的，往往存在一定厚度的间隙。在两层金属飞片之间存在间隙情况下，当应力波传播至间隙部位时，会增加间隙部位闭合压实的动力学过程。该动力学过程会直接影响金属飞片内的加载压力波形及强度，进而间接影响飞片自由面速度和终态动能等。

在以往的研究中，李涛等^[1]设计了钢-钢间含固定厚度初始间隙的爆轰驱动实验，该实验自由面测速结果验证了初始间隙的存在对飞片自由面速度具有较大影响。贺年丰等^[2]在爆轰加载钢-锡双层飞片实验中，在钢-锡飞片间设置了不同厚度的间隙，自由面测速和X光实验结果均表明，初始间隙厚度不均匀会引起飞片界面扰动，且该扰动幅值较大、不可忽略。已有的相关实验研究^[4-8]表明，此类金属飞片界面扰动在强冲击加载下可能诱发界面不稳定性的产生，且随着扰动的不断发展甚至会引起飞片断裂。此外，在爆轰驱动问题中，金属飞片间的初始间隙亦会对其他一些方面造成影响，如爆炸焊接效能^[3,9-10]及炸药推动飞片的做功能力^[11]等。为了更好地认识由间隙因素造成的一些非预期物理现象，需要深入研究间隙对爆轰驱动金属飞片动力学过程的影响规律。

本文中将首先对平面爆轰驱动含初始间隙的双层钢飞片实验进行简化建模和模拟，然后基于该简化模型进行间隙厚度的小范围外推，最后给出间隙对此类爆轰驱动金属飞片问题影响规律的初步认识。

1. 爆轰驱动含间隙双层钢飞片实验的模拟

李涛等^[1]在开展的含初始装配间隙的平面爆轰实验中，使用平面波透镜驱动 ${\varnothing }$ 100 mm×20 mm的RHT-901炸药爆轰，炸药爆轰波驱动两层 ${\varnothing }$ 100 mm×3.5 mm的45钢飞片，两层钢飞片间存在厚0.5 mm的装配间隙（见图1）。实验测量获得了 ${\varnothing }$ 30 mm与 ${\varnothing }$ 60 mm位置处的飞片自由面速度。

图 1 平面爆轰驱动含间隙双层钢飞片实验^[1]示意图

Figure 1. Schematic diagram of the plane detonation driven double-layer flyers experiment^[1]

下载: 全尺寸图片幻灯片

使用有限元程序DYNA2D^[12]进行二维轴对称建模及模拟。考虑到平面波透镜真实建模可能引入的额外不确定性^[13-14]，采用孙承纬^[15]提出的近似方法将平面波透镜简化为等效厚度12 mm的RHT-901炸药。简化后的计算模型如图2所示，其中R方向为柱面方向、Z方向为平面方向。在简化模型中，炸药起爆方式为Z=−12 mm的底面同步时间起爆。炸药-内层飞片、内层飞片-外层飞片之间的滑移线均使用开穴算法^[16]。并根据实验实际测量值^[1]设置简化模型中的间隙厚度为0.48 mm。

图 2 实验模型的简化建模及含间隙飞片的局部网格划分

Figure 2. Simplified modeling result of experiment setup and the mesh generation of flyers with clearance

下载: 全尺寸图片幻灯片

在模拟中，金属使用如下形式的Mie-Grüneisen状态方程：

$\left. p = {\rho _0}{c^2}\mu \left[ {1 + \left(1 - \frac{{{\gamma _0}}}{2}\right)\mu - \frac{\lambda }{2}{\mu ^2}} \right]\right/{\left[ {1 - ({S_1} - 1)\mu - {S_2}{\mu ^2}/(\mu + 1) - {S_3}{\mu ^3}/{{(\mu + 1)}^2}} \right]^2} + \left( {{\gamma _0} + \lambda \mu } \right)e{\rho _0}$

(1)

式中： $\ \mu = {\rho }/{{{\rho _0}}} - 1$ ，其中 ${\ \rho }_{0}$ 为初始密度， ${c}_{0}$ 为零压声速， ${S}_{1}\mathrm{、}{S}_{2}\mathrm{、}{S}_{3}$ 、 ${\gamma }_{0}$ 、 $\lambda$ 为状态方程参数。

金属使用适用于较高应变率的Steinberg-Guinan (SG)本构模型：

$G = {G_0}\left[ {1 + bp{{\left( {\frac{\rho }{{{\rho _0}}}} \right)}^{1/3}} - h\left( {T - 300} \right)} \right]$

(2)

$Y = {Y_0}{\left( {1 + \beta {\varepsilon _{\text{p}}}} \right)^n}\left[ {1 + bp{{\left( {\frac{\rho }{{{\rho _0}}}} \right)}^{1/3}} - h\left( {T - 300} \right)} \right]\qquad {Y_0}{\left( {1 + \beta {\varepsilon _{\text{p}}}} \right)^n} {\text{≤}} {Y_{\max }}$

(3)

$T = {{\left( {E - {\varepsilon _0}} \right)}/ {3R}}\qquad {\varepsilon _0} = {\varepsilon _{00}} + {\varepsilon _{01}}x + {\varepsilon _{02}}{x^2} + {\varepsilon _{03}}{x^3} + {\varepsilon _{04}}{x^4},{\text{ }}x = 1 - {\rho _0}/\rho$

(4)

${T_{\text{m}}} = {T_{{\text{m}}0}}{\left( {\frac{{{\rho _0}}}{\rho }} \right)^{2/3}}\exp \left[ {2{\gamma _0}\left( {1 - V} \right)} \right] \qquad 若 T {\text{＞}} {T_{\text{m}}}，取Y = G = 0$

(5)

式中： ${G}_{0}$ 为剪切模量，b为剪切模量硬化系数，h为温度软化系数， $T$ 为材料温度， ${Y}_{0}$ 为屈服强度， $\beta$ 为硬化系数，ε_p为等效塑性应变， $n$ 为硬化指数， ${Y}_{\mathrm{m}\mathrm{a}\mathrm{x}}$ 为最大屈服强度， ${T}_{\mathrm{m}}$ 为熔化温度， ${T}_{\mathrm{m}0}$ 为初始熔化温度，V为相对体积。

金属使用Void-Growth (VG)损伤模型^[17]。在VG损伤模型中，材料损伤度定义为D，具体计算形式如下：

$\alpha = \frac{1}{{1 - D}}$

(6)

$\Delta p = p - \frac{{{\alpha_{\rm{s}}}}}{\alpha }\ln \left( {\frac{\alpha }{{\alpha - 1}}} \right)$

(7)

$\dot{\alpha }=\left\{\begin{array}{ll}0 & \Delta p{\text{≤}} 0\\ \dfrac{{\left({\alpha }_{0}-1\right)}^{2/3}}{\eta }\alpha {\left(\alpha -1\right)}^{1/3}\Delta p &\Delta p {\text{＞}} 0\end{array}\right.$

(8)

$\dot D = {\dot \alpha }/ {{\alpha ^2}}$

(9)

式中： $\alpha$ 为材料孔隙度， ${\alpha }_{0}$ 为初始孔隙度， ${\alpha }_{\mathrm{s}}\mathrm{、}\eta$ 为材料参数。定义 ${D}_{\mathrm{c}\mathrm{u}\mathrm{t}}$ 为临界损伤度，当 $D {\text{＞}} {D}_{\mathrm{c}\mathrm{u}\mathrm{t}}$ 时，认为材料断裂，并以材料被再次压缩到临界断裂状态点作为材料再压实判据。

炸药使用Jones-Wikins-Lee (JWL)状态方程来描述：

$p = A\left( {1 - \frac{\omega }{{{R_1}V}}} \right){{\text{e}}^{ - {R_1}V}} + B\left( {1 - \frac{\omega }{{{R_2}V}}} \right){{\text{e}}^{ - {R_2}V}} + \frac{{\omega E}}{V}$

(10)

式中： $A\mathrm{、}B\mathrm{、}{R}_{1}\mathrm{、}{R}_{2}\mathrm{、}\omega$ 为JWL状态方程参数， $E$ 为体积爆热。具体计算参数参见表1～3，其中 ${p}_{\mathrm{C}\mathrm{J}}\mathrm{、}{D}_{\mathrm{C}\mathrm{J}}$ 为炸药爆压和爆速。

表 1 RHT-901炸药JWL参数^[18]

Table 1. The JWL parameters of RHT-901^[18]

ρ₀/(g·cm⁻³)	p_CJ/GPa	D_CJ/(km·s⁻¹)	A/GPa	B/GPa	R₁	R₂	ω	E/(kJ·cm⁻³)
1.7	27.0	7.8	755.975	13.175	4.8	0.58	0.34	8.5

下载: 导出CSV

| 显示表格

表 2 45钢材料模型参数^[19-21]

Table 2. The material parameters of steel 45^[19-21]

ρ₀/(g·cm⁻³)	c₀/(km·s⁻¹)	S₁	S₂	S₃	γ₀	λ	G₀/GPa	Y₀/GPa	Y_max/GPa	T_m0/K
7.85	4.57	1.49	0	0	2.17	0.43	81.8	0.355	2.0	2380

下载: 导出CSV

| 显示表格

表 3 45钢的VG损伤模型参数

Table 3. The VG parameters of steel 45

${\alpha }_{0}$	${\alpha }_{\mathrm{s}}/\mathrm{G}\mathrm{P}\mathrm{a}$	$\eta$ /（Pa·s）	${D}_{\mathrm{c}\mathrm{u}\mathrm{t}}$
1.0001	0.3	10	0.10

下载: 导出CSV

| 显示表格

图3给出了简化模型模拟得到的两层钢飞片在不同时刻的密度图。在模拟观测时间内， ${\varnothing }$ 30 mm与 ${\varnothing }$ 60 mm位置上的自由面速度曲线重合，简化模型在约 ${\varnothing }$ 65 mm范围内为近似一维加载区。由于模拟中对起爆装置进行了简化，爆轰波初始入射钢飞片的时刻在模拟与实验中存在相对偏移。在将模拟与实验中的首次速度起跳峰值时刻对齐后，图4(a)给出了 ${\varnothing }$ 30 mm自由面速度的实验结果和模拟结果对比。图4(b)给出了R=0轴线上钢飞片的Z轴坐标随时间的变化图像，其中伪彩标识了对应时间、空间位置上钢飞片的内部压力值。由模拟与实验的对比可知，模拟结果能够很好地符合实验的首次高度、二次加载时间间隔及二次平台高度等关键物理信息，简化建模及模拟能够正确反映爆轰冲击含间隙双层钢飞片动力学过程。

图 3 模拟得到的含初始间隙的两层钢飞片在不同时刻的密度图

Figure 3. The density of double-layer steel flyer with initial clearance at different times

下载: 全尺寸图片幻灯片

图 4 爆轰驱动含间隙双层钢板的自由面速度对比及R=0轴线上的压力随时间变化情况

Figure 4. The free surface velocity of the double-layer flyers with clearance driven by detonation and the pressure with time at R = 0

下载: 全尺寸图片幻灯片

由图4可知，该爆轰驱动两层飞片模型的动力学过程及波系关系相对简单。

(1) 首先，首次加载波在约5.6 µs时刻（以炸药起爆时刻作为零时）作用于外飞片自由面，形成了图4(a)中的速度曲线首次起跳。其中，首次加载压力峰值约32 GPa、自由面速度首次起跳高度约1.43 km/s。

(2) 之后，二次加载波在约6.5 µs处作用于外飞片自由面，形成了图4(a)中的速度曲线二次起跳。其中，首次与二次加载速度峰值的时间间隔约0.94 µs、二次加载后的速度平台高度约1.61 km/s。由图4(b)可以清晰地看到，二次加载波是间隙闭合碰撞形成的反射冲击波与爆轰产物持续加载波对碰所形成的（图4(b)中红色菱形区域）。但是，由于二次加载波需要经过首次加载波在界面反射的稀疏波区后才作用于外飞片自由面，此时其峰值已降至约7.5 GPa。

(3) 最后，在外飞片后界面反射的稀疏波与前界面反射的稀疏波的共同拉伸作用下，在约7.3 µs时刻外飞片内形成层裂现象。层裂现象在图4(a)中体现为在7.3 µs后自由面速度基本不变，在图4(b)中体现为在7.3 µs后外飞片内部拉伸出的真空缝隙。

2. 间隙厚度变化对飞片加载动力学过程的影响

在不同的工程问题中间隙的影响不尽相同，本文中仅针对上述爆轰驱动双层钢飞片这一特定工程模型进行讨论。通过前面的模拟分析，对该模型中的波系关系、自由面首次和二次加载的来源等有了较清晰的认识。下面将通过在小范围内改变间隙大小，来进一步分析间隙厚度变化对该模型动力学过程的影响。

首先，间隙厚度不同影响的最直观体现就是飞片自由面速度。在保持图2中计算模型及参数不变的情况下，在初始间隙0.48 mm附近的小范围内改变间隙厚度，考察间隙厚度对自由面速度的影响。7个分析模型中的间隙厚度分别设置为W = 0.10, 0.25, 0.37, 0.48, 0.75, 1.00, 1.25 mm。通过对应模型的模拟得到的 ${\varnothing }$ 30 mm位置自由面速度曲线如图5(a)所示。图5(b)给出了图5(a)中首次及二次加载自由面速度峰值随间隙厚度的变化关系。可以看到，外层钢飞片在一定时间后均出现了层裂现象，但小间隙情况下为单次加载后层裂（二次加载不明显），而在较大间隙情况下出现了二次加载后层裂的现象。间隙厚度变化对飞片自由面速度的影响较大：随着间隙的增加，首次速度峰值逐渐降低至基本不变，在考察范围内，首次加载峰值速度的影响幅值约0.09 km/s；随着间隙的增加，二次速度峰值逐渐增大至基本不变，在考察范围内，二次加载峰值速度的影响幅值达到0.5 km/s。显然，就飞片自由面速度而言，间隙对此类问题的影响不可忽略。

图 5 不同间隙厚度下爆轰驱动双飞片模型的自由面速度及首次和二次加载自由面速度峰值

Figure 5. Free surface velocity and the peak value of free surface velocity under first and second loading with different clearance thickness

下载: 全尺寸图片幻灯片

之后，针对自由面速度所体现出的随着间隙增大二次加载更明显这一现象，开展首次和二次加载动力学过程分析。图6给出了间隙较小和较大（W = 0.10, 1.25 mm）两种情况下的钢飞片Z轴坐标随时间变化的图像，其中伪彩标识了对应时间、空间位置上的压力值。

图 6 爆轰驱动双飞片模型在R=0轴线上的压力随时间变化图像

Figure 6. Pressure varied with time at R = 0 in the detonation driven double-layer flyer model

下载: 全尺寸图片幻灯片

由图6(a)可知，小间隙对爆轰加载形成的初始入射波的影响较小，飞片外界面的首次加载波与初始入射波的形貌基本相同（均为类三角波）。显然，当间隙为零时，该过程退化为初始入射波在材料内部的正常衰减过程，飞片外界面的首次加载波就是衰减后的初始入射波。而且，在间隙较小的情况下，自由面的稀疏效应发展时间较短，两层飞片碰撞时的速度差异亦较小，所以后期在飞片内部无法追赶碰撞形成高压二次加载。由图6(b)可知，初始入射波在内层钢飞片的外界面形成了层裂（因为间隙较大，拉伸稀疏区有足够时间发展为层裂片）。因此外层钢飞片的首次加载波是由内层钢飞片的层裂片撞击形成的。在图6(b)中可以清晰地看到，外层钢飞片的首次加载为近似方波加载，这就使得其自由面速度出现与类三角波加载不同的物理现象。图6(b)中的二次加载来源与图4(b)中基本相同：在内层钢飞片与外层钢飞片撞击形成的反射激波与爆轰产物加载入射波碰撞后成了高压二次加载。对比图4(b)和图6(b)可知，随着间隙增大，二次加载高压区覆盖范围逐渐增加。

因此，随着间隙逐渐增加，初始入射压力波形由类三角波逐渐变为近似方波；二次加载波更容易形成。为了深入理解上述认识并剖析原因，进一步对图5中自由面速度曲线进行分类对比。图7(a)给出了W = 0.10, 0.25, 0.37, 0.48 mm等4种情况下的自由面速度局部细节对比；图7(b)给出了W = 0.75, 1.00, 1.25 mm等3种情况下的自由面速度局部细节对比。

图 7 不同尺度间隙下爆轰驱动双飞片模型自由面速度曲线

Figure 7. Free surface velocity of detonation driven double-layer flyer model with different clearance thickness

下载: 全尺寸图片幻灯片

将间隙厚度以0.5 mm作为分界厚度，分为图7(a)～(b)中两个阶段进行对比，得到如下认识。

(1) 图7(a)中的4种间隙厚度较小，内飞片在间隙一侧未能形成层裂片，在此情况下，内飞片与外飞片发生碰撞。因此自由面首次加载波仍旧为类三角波，首次加载后其自由面卸载历程基本相同（在图7(a)中体现为几条曲线速度下降段基本重合）。由于内飞片在间隙一侧未形成速度稳定不变的层裂片，因此，在此阶段内，随着间隙厚度的增加，内飞片与外飞片的碰撞速度增大、二次加载压力增高（在图7(a)中体现为二次加载的自由面速度峰值增大）。

(2) 图7(b)中的3种间隙厚度较大，内飞片在间隙一侧形成了层裂片。由于内飞片的层裂片形成后其厚度、速度均不变，因此该层裂片与外飞片碰撞形成的（近似方波）波形、强度基本不变。所以，此阶段内，随着间隙厚度的增大，其首次起跳速度的平台高度和下降速率基本相同（在图7(b)中体现为首次加载后速度曲线的近似平移变化）。同样，由于内飞片形成的稳定层裂片的速度不变，因此，在不同间隙厚度下，层裂片与外层钢飞片碰撞形成的首次、二次波峰值压力基本不变（在图7(b)中体现为首次、二次加载后的速度峰值基本相同）。

就上述爆轰驱动双层钢飞片模型动力学过程而言，可将不同间隙大小影响分为两个阶段，其分界判据是在爆轰加载后内层金属飞片是否能够在间隙部位发展为脱体层裂片。在间隙较小的情况下，内层飞片在间隙一侧无法发展为层裂片，在此阶段内，随着间隙的增大，外层飞片的首次加载峰值压力降低、二次加载的峰值压力升高；在间隙较大的情况下，内层飞片在间隙一侧可以形成厚度不变、速度稳定的层裂片，在此阶段内，随着间隙的增大，外层飞片的首次加载与二次加载的峰值压力均基本不变，但首次与二次加载的时间间隔缩短。以此认识为指导，可以在模型中通过精确控制间隙大小来达到调整首次波形、控制二次波强度等目的，也可以在未知飞片内部间隙大小的情况下，通过测量自由面速度波形初步判断两层飞片间的间隙情况。

3. 结　论

在爆轰冲击含间隙双层钢飞片实验研究的基础上，开展了不同间隙厚度影响的模拟分析研究，得到以下结论。

(1) 随着间隙的增大，飞片自由面首次起跳速度峰值由逐渐降低至基本不变、二次起跳速度峰值由逐渐增加至基本不变。

(2) 在对间隙大小进行分类对比后，认为间隙大小直接影响的是在爆轰加载后内层钢飞片是否有足够时间在间隙部位发展为脱体层裂片：间隙较小时，内飞片在间隙一侧无法发展为层裂片，在此阶段内，随着间隙增大，首次加载自由面峰值速度降低、二次加载自由面峰值速度升高；间隙较大时，内飞片在间隙一侧可以形成厚度不变、速度稳定的层裂片，在此阶段内，随着间隙厚度增加，首次与二次加载的自由面峰值速度均基本不变、首次与二次加载的时间间隔缩短。

需要说明的是，本文中间隙影响的规律认识是建立在具有一定置信度模拟结果的小范围外推基础上的，较大间隙下模拟结果的可靠性仍需要实验验证。将本文中真空间隙填充为常压空气后得到的模拟结果基本无差异，常压空气间隙造成的影响规律不变。内飞片材料与间隙之间的关系更为密切，内飞片材料的改变可能影响本文中部分规律认识。

图 1 平面爆轰驱动含间隙双层钢飞片实验^[1]示意图

Figure 1. Schematic diagram of the plane detonation driven double-layer flyers experiment^[1]

下载: 全尺寸图片幻灯片

图 2 实验模型的简化建模及含间隙飞片的局部网格划分

Figure 2. Simplified modeling result of experiment setup and the mesh generation of flyers with clearance

下载: 全尺寸图片幻灯片

图 3 模拟得到的含初始间隙的两层钢飞片在不同时刻的密度图

Figure 3. The density of double-layer steel flyer with initial clearance at different times

下载: 全尺寸图片幻灯片

图 4 爆轰驱动含间隙双层钢板的自由面速度对比及R=0轴线上的压力随时间变化情况

Figure 4. The free surface velocity of the double-layer flyers with clearance driven by detonation and the pressure with time at R = 0

下载: 全尺寸图片幻灯片

图 5 不同间隙厚度下爆轰驱动双飞片模型的自由面速度及首次和二次加载自由面速度峰值

Figure 5. Free surface velocity and the peak value of free surface velocity under first and second loading with different clearance thickness

下载: 全尺寸图片幻灯片

图 6 爆轰驱动双飞片模型在R=0轴线上的压力随时间变化图像

Figure 6. Pressure varied with time at R = 0 in the detonation driven double-layer flyer model

下载: 全尺寸图片幻灯片

图 7 不同尺度间隙下爆轰驱动双飞片模型自由面速度曲线

Figure 7. Free surface velocity of detonation driven double-layer flyer model with different clearance thickness

下载: 全尺寸图片幻灯片

表 1 RHT-901炸药JWL参数^[18]

Table 1. The JWL parameters of RHT-901^[18]

ρ₀/(g·cm⁻³) p_CJ/GPa D_CJ/(km·s⁻¹) A/GPa B/GPa R₁ R₂ ω E/(kJ·cm⁻³)

1.7 27.0 7.8 755.975 13.175 4.8 0.58 0.34 8.5

下载: 导出CSV

表 2 45钢材料模型参数^[19-21]

Table 2. The material parameters of steel 45^[19-21]

ρ₀/(g·cm⁻³) c₀/(km·s⁻¹) S₁ S₂ S₃ γ₀ λ G₀/GPa Y₀/GPa Y_max/GPa T_m0/K

7.85 4.57 1.49 0 0 2.17 0.43 81.8 0.355 2.0 2380

下载: 导出CSV

表 3 45钢的VG损伤模型参数

Table 3. The VG parameters of steel 45

${\alpha }_{0}$ ${\alpha }_{\mathrm{s}}/\mathrm{G}\mathrm{P}\mathrm{a}$ $\eta$ /（Pa·s） ${D}_{\mathrm{c}\mathrm{u}\mathrm{t}}$

1.0001 0.3 10 0.10

下载: 导出CSV

参考文献(21)

[1]	李涛, 刘明涛, 王晓燕, 等. 装配垫层与间隙对爆轰加载下金属飞片运动特征的影响 [J]. 高压物理学报, 2018, 32(4): 044202. DOI: 10.11858/gywlxb.20170576. LI T, LIU M T, WANG X Y, et al. Effects of explosive device with foam cushion and air clearance on kinetic characteristic of steel flyer under detonation loading [J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 044202. DOI: 10.11858/gywlxb.20170576.
[2]	贺年丰, 张绍龙, 洪仁楷, 等. 间隙对金属锡爆轰加载过程的影响 [J]. 爆炸与冲击, 2021, 41(1): 012101. DOI: 10.11883/bzycj-2020-0054. HE N F, ZHANG S L, HONG R K, et al. Effects of gap on the explosive loading process of tin [J]. Explosion and Shock Waves, 2021, 41(1): 012101. DOI: 10.11883/bzycj-2020-0054.
[3]	李雪交, 马宏昊, 沈兆武. 铝合金与槽型界面钢板的爆炸焊接 [J]. 爆炸与冲击, 2016, 36(5): 640–647. DOI: 10.11883/1001-1455(2016)05-0640-08. LI X J, MA H H, SHEN Z W. Explosive welding of interface between aluminum alloy and steel plate with dovetail grooves [J]. Explosion and Shock Waves, 2016, 36(5): 640–647. DOI: 10.11883/1001-1455(2016)05-0640-08.
[4]	何长江, 周海兵, 杭义洪. 爆轰驱动金属铝界面不稳定性的数值分析 [J]. 中国科学 G辑: 物理学力学天文学, 2010, 53(2): 195–198. DOI: 10.1007/s11433-009-0261-4. HE C J, ZHOU H B, HANG Y H. A numerical study on Rayleigh-Taylor instability of aluminum plates driven by detonation [J]. Science China Physics, Mechanics and Astronomy, 2010, 53(2): 195–198. DOI: 10.1007/s11433-009-0261-4.
[5]	殷建伟, 潘昊, 吴子辉, 等. 爆轰驱动Cu界面的Richtmyer-Meshkov扰动增长稳定性 [J]. 物理学报, 2017, 66(20): 204701. DOI: 10.7498/aps.66.204701. YIN J W, PAN H, WU Z H, et al. Stability analysis of interfacial Richtmyer-Meshkov flow of explosion-driven copper interface [J]. Acta Physica Sinica, 2017, 66(20): 204701. DOI: 10.7498/aps.66.204701.
[6]	刘军, 冯其京, 周海兵. 柱面内爆驱动金属界面不稳定性的数值模拟研究 [J]. 物理学报, 2014, 63(15): 155201. DOI: 10.7498/aps.63.155201. LIU J, FENG Q J, ZHOU H B. Simulation study of interface instability in metals driven by cylindrical implosion [J]. Acta Physica Sinica, 2014, 63(15): 155201. DOI: 10.7498/aps.63.155201.
[7]	张维岩, 叶文华, 吴俊峰, 等. 激光间接驱动聚变内爆流体不稳定性研究 [J]. 中国科学: 物理学力学天文学, 2014, 44(1): 1–23. DOI: 10.1360/SSPMA2013-00039. ZHANG W Y, YE W H, WU J F, et al. Hydrodynamic instabilities of laser indirect-drive inertial-confinement-fusion implosion [J]. Scientia Sinica Physica, Mechanica & Astronomica, 2014, 44(1): 1–23. DOI: 10.1360/SSPMA2013-00039.
[8]	王涛, 汪兵, 林健宇, 等. 内爆加载下界面不稳定性和湍流混合数值模拟研究 [J]. 中国科学: 物理学力学天文学, 2020, 50(10): 104704. DOI: 10.1360/SSPMA-2019-0420. WANG T, WANG B, LIN J Y, et al. Numerical investigations of interface instability and turbulent mixing driven by implosion [J]. Scientia Sinica Physica, Mechanica & Astronomica, 2020, 50(10): 104704. DOI: 10.1360/SSPMA-2019-0420.
[9]	曾翔宇, 李晓杰, 曹景祥, 等. 材料强度对爆炸焊接结合界面的影响 [J]. 爆炸与冲击, 2019, 39(5): 137–143. DOI: 10.11883/bzycj-2018-0400. ZENG X Y, LI X J, CAO J X, et al. Interface characteristics of explosive welding for different strength plates [J]. Explosion and Shock Waves, 2019, 39(5): 137–143. DOI: 10.11883/bzycj-2018-0400.
[10]	张登霞, 李国豪, 周之洪, 等. 材料强度在爆炸焊接界面波形成过程中的作用 [J]. 力学学报, 1984, 16(1): 73–80. ZHANG D X, LI G H, ZHOU Z H, et al. Effect of material strength on forming process of explosive welding interface wave [J]. Acta Mechanica Sinica, 1984, 16(1): 73–80.
[11]	虞德水, 赵锋, 谭多望, 等. JOB-9003和JB-9014炸药平面爆轰驱动飞片的对比研究 [J]. 爆炸与冲击, 2006, 26(2): 140–144. DOI: 10.11883/1001-1455(2006)02-0140-05. YU D S, ZHAO F, TAN D W, et al. Experimental studies on detonation driving behaviorof JOB-9003 and JB-9014 slab explosives [J]. Explosion and Shock Waves, 2006, 26(2): 140–144. DOI: 10.11883/1001-1455(2006)02-0140-05.
[12]	WHIRLEY R G, ENGELMANN B E. DYNA2D: a nonlinear, explicit, two-dimensional finite element code for solid mechanics: user manual: UCRL-MA-110630 [R]. Washington: Lawrence Livermore National Laboratory, 1992.
[13]	袁帅, 周显明, 李平. 液体炸药-惰性材料组合式平面波透镜数值模拟 [J]. 兵工学报, 2010, 31(S1): 159–161. YUAN S, ZHOU X M, LI P. Simulation of plane wave lens with the combinition of liquid explosive and inert material [J]. Acta Armamentarii, 2010, 31(S1): 159–161.
[14]	韩勇, 魏智勇, 黄毅民, 等. 平面波透镜实验与数值模拟 [J]. 含能材料, 2007, 15(5): 471–473. DOI: 10.3969/j.issn.1006-9941.2007.05.008. HAN Y, WEI Z Y, HUANG Y M, et al. Experiment and numerical simulation of plane wave lens [J]. Chinese Journal of Energetic Materials, 2007, 15(5): 471–473. DOI: 10.3969/j.issn.1006-9941.2007.05.008.
[15]	孙承纬. 炸药平面波透镜的有效药量 [C]//爆轰研究论文集, 第3卷. 绵阳: 中国工程物理研究院流体物理研究所, 1998: 307-316.
[16]	林文洲, 林忠, 刘全. 非结构多边形网格滑移线开穴算法 [J]. 计算物理, 2017, 34(3): 273–282. DOI: 10.3969/j.issn.1001-246X.2017.03.004. LIN W Z, LI Z, LIU Q. An open void method of slide line on unstructured N-polygon grids [J]. Chinese Journal of Computational Physics, 2017, 34(3): 273–282. DOI: 10.3969/j.issn.1001-246X.2017.03.004.
[17]	JOHNSON JN. Dynamic fracture and spallation in ductile solids [J]. Journal of Applied Physics, 1981, 52(4): 2812–2825. DOI: 10.1063/1.329011.
[18]	董海山, 周芬芬. 高能炸药及相关物性能[M]. 北京: 科学出版社, 1989.
[19]	STEINBERG D J. Equation of state and strength properties of selected materials: UCRL-MA-106439 [R]. Washington: Lawrence Livermore National Laboratory, 1996.
[20]	张林, 张祖根, 秦晓云, 等. D6A、921和45钢的动态破坏与低压冲击特性 [J]. 高压物理学报, 2003, 17(4): 305–310. DOI: 10.11858/gywlxb.2003.04.011. ZHANG L, ZHANG Z G, QIN X Y, et al. Dynamic fracture and mechanical property of D6A, 921 and 45 steels under low shock pressure [J]. Chinese Journal of High Pressure Physics, 2003, 17(4): 305–310. DOI: 10.11858/gywlxb.2003.04.011.
[21]	胡昌明, 贺红亮, 胡时胜. 45号钢的动态力学性能研究 [J]. 爆炸与冲击, 2003, 23(2): 188–192. HU C M, HE H L, HU S S. A study on dynamic mechancial behaviors of 45 steel [J]. Explosion and Shock Waves, 2003, 23(2): 188–192.

施引文献

期刊类型引用(2)

1.	刘军，殷建伟，张凤国. 基于简易冲击分解模型的爆轰驱动硅橡胶数值模拟及实验解读. 爆炸与冲击. 2025(01): 125-134 . 本站查看
2.	周蕊，李理，田保林. 爆轰驱动多介质问题的Lagrange多分区自适应数值模拟研究. 力学学报. 2023(11): 2675-2692 . 百度学术