扰流板结构参数对冲击波衰减特性影响的仿真分析

赵蓓蕾; 赵继广; 崔村燕; 柳宁远; 王岩; 辛腾达; 王亚琦

doi:10.11858/gywlxb.20170585

扰流板结构参数对冲击波衰减特性影响的仿真分析

doi: 10.11858/gywlxb.20170585

赵蓓蕾^{1, 2,},
赵继广³,
崔村燕²,
柳宁远^{1, 2},
王岩^{1, 2},
辛腾达^{1, 2},
王亚琦⁴

1.
装备学院研究生管理大队，北京 101416
2.
装备学院航天装备系，北京 101416
3.
装备学院装备发展战略研究所，北京 101416
4.
北京空间信息中继传输技术研究中心，北京 100094

详细信息

作者简介:
赵蓓蕾(1991—), 女，博士研究生，主要从事航天任务分析与设计研究.E-mail:1178685385@qq.com

中图分类号: O383.1;V551.12
计量
- 文章访问数: 6646
- HTML全文浏览量: 2874
- PDF下载量: 246
出版历程
- 收稿日期: 2017-05-26
- 修回日期: 2017-06-26

Simulation Analysis of Influence of Spoiler Structural Parameters on Shock Wave Attenuation Characteristics

ZHAO Beilei^{1, 2
,},
ZHAO Jiguang³,
CUI Cunyan²,
LIU Ningyuan^{1, 2},
WANG Yan^{1, 2},
XIN Tengda^{1, 2},
WANG Yaqi⁴

1.
Company of Postgraduate Management, Equipment Academy, Beijing 101416, China
2.
Department of Aerospace Equipment, Equipment Academy, Beijing 101416, China
3.
Equipment Development Strategy Research Institute, Equipment Academy, Beijing 101416, China
4.
Beijing Space Information Ralay and Transmission Technology Research Center, Beijing 100094, China

摘要

摘要: 在坑道内设置扰流板是加速爆炸冲击波衰减的有效方法。为探究扰流板结构参数对冲击波衰减特性的影响，采用ANSYS/LS-DYNA有限元软件进行仿真研究。基于流固耦合算法，建立Schardin实验等比模型，所得仿真结果与实验结果具有良好的一致性，验证了仿真模型的有效性。以某发射井坑道为研究对象，在矩形扰流板宽度一定的情况下，研究了扰流板厚度、倾角及间距对冲击波衰减规律的影响。结果表明：保持其他参数不变，随着扰流板厚度增加，冲击波超压衰减越来越明显；扰流板厚度为40 cm时，倾角为105°、间距为6 m最有利于冲击波衰减。研究结果可以为坑道防护设计提供有价值的参考。
- 扰流板 /
- 爆炸冲击波 /
- 衰减 /
- 数值仿真
Abstract: Setting spoilers in the tunnel is an effective way to accelerate its shock wave attenuation.In order to investigate the influence of the spoiler structural parameters on its shock wave's attenuation, numerical simulation was carried out using the ANSYS/LS-DYNA finite element software.First, based on the fluid-solid coupling algorithm, the equal scale model of Schardin's experiment was established.The simulation results were found to be in good agreement with the experimental results, which verified the validity of the simulation model.Then the numerical method was used to investigate the influence of the spoiler thickness, inclination angle and interval on the shock wave attenuation in a silo tunnel when the rectangular spoiler width was a constant.The results show that when the other parameters remain the same, with the increase of the spoiler thickness, the overpressure attenuation of the shock wave becomes increasingly more obvious.When the spoiler thickness is 40 cm, the spoiler with the inclination angle of 105° and the interval of 6 m is the most beneficial condition for the shock wave attenuation.These results can provide valuable reference for the design of the tunnel protection.
- spoiler /
- shock wave /
- attenuation /
- numerical simulation

HTML全文

当冲击波通过两种不同物质的界面时，界面扰动在冲击加载作用下会增长演化，即Richtmyer-Meshkov（RM）不稳定性^[1]。RM不稳定性是目前惯性约束核聚变研究的关键基础问题，其挑战主要来源于极端热动力学条件（压力远大于100 GPa，温度远高于10⁴ K）^[2-3]。目前，尚未发现针对极端条件的新RM不稳定性理论模型，进行相关实验需要高昂的费用，并且数据测量极为困难，所获数据十分有限。

数值模拟方面，对于温稠密物质及高密度等离子体，常规的流体动力学模拟缺少有效的模型和输运参数。不过，从原子尺度，采用分子动力学的建模和模拟方法正成为一种有效的技术手段。这类方法从描述原子或核子的基本相互作用势出发，不要求预设任何宏观模型、本构关系、宏观输运系数等，通过计算大量原子或核子的运动，直接模拟各种能量密度范围内的物质运动，可以涉及非平衡态、非连续效应、多组分/多相效应、电离以及与材料混合相关的随机反应等各种复杂物理机制。目前，在极端热力学条件下材料混合问题的理论研究中已经普遍采用分子动力学方法进行计算^[4-8]，特别是基于量子力学方程的第一性原理分子动力学方法（First Principles Molecular Dynamics，FPMD）^[9-19]，且正有计划地发展大规模分子动力学程序^[4]。然而，利用经典势函数的分子动力学模拟受到势函数的限制，一般只能模拟常规条件下或非电离状态下的RM不稳定性；基于密度泛函理论的QMD（Quantum Molecular Dynamics）方法可以计算极端条件下的物性参数和状态方程，但是由于计算量巨大，无法直接模拟RM不稳定过程。Su等^[20-21]提出将单电子显式处理成波包（Wave Packet，有位置，大小可变）的模型，从而发展出一种新的具有第一性原理特质的分子动力学方法—电子力场（Electron Force Field，eFF）分子动力学。近年来，该方法得到不断发展，已经具备处理10⁶量级离子和电子及长时间尺度的能力^[22]。

近期本研究团队利用该方法^[23-24]，从微观角度开展了不同热动力学条件下（3～30 km/s高速冲击加载）金属锂和气态氢材料混合的分子动力学数值模拟，结果发现极端加载条件下RM不稳定性的演化呈现新的特征，即振幅增长的线性段加长以及非线性段振幅增长速度高于无电离效应工况，这一结果得到了Nova实验结果^[25-26]的印证。其机理在于：极端条件下激波压缩导致物质电离，由于电子和离子扩散存在差异，电子将扩散到激波之前，形成电荷分离区，进而诱导产生附加电场，由附加电场引起的界面物质附加加速度会显著影响界面演化。为深入理解极端条件下RM不稳定性的演化特征，特别是附加电场对界面物质的加速效应，本研究继续采用eFF MD方法，模拟强激波作用下电荷分离情况，对电荷沿激波和界面运动方向的分布进行统计分析，得到电荷分离区强度和宽度随激波强度变化的规律，并进一步获得附加电场影响下界面附近金属侧物质产生的加速度，通过拟合，获得附加电场引起的界面附加加速度的计算模型。

1. 方法与模型

1.1 MD方法

MD方法的本质是通过牛顿第二定律求解每个原子（或分子）位置、速度、受力等基本物理量随时间的变化，其控制方程为

${{{F}}_i} = - {\nabla _i}U = - \left( {\frac{\partial }{{\partial {x_i}}}{{i}} + \frac{\partial }{{\partial {y_i}}}{{j}} + \frac{\partial }{{\partial {z_i}}}{{k}}} \right)U$

(1)

式中：i代表粒子编号；F_i为系统中其他粒子施加在粒子i上的力的总和，可以由粒子间的势函数 $U$ 求梯度得到。势函数的表达形式为

$U = U\left( {{{{r}}_1},{{{r}}_2}, \cdot \cdot \cdot ,{{{r}}_i}, \cdot \cdot \cdot ,{{{r}}_n}} \right)$

(2)

式中：r_i为第i个原子的坐标。获得施加在粒子i上的力之后，它相应的加速度a_i、速度v_i以及位置都可以由下列力学式计算得到

$\left\{ \begin{array}{l} {{{F}}_i} = {m_i}{{{a}}_i}\\ {{{v}}_i} = {{{v}}_i^0} + {{{a}}_i}t\\ {{{r}}_i} = {{{r}}_i^0} + {{{v}}_i^0}t + \dfrac{1}{2}{{{a}}_i}{t^2} \end{array} \right.$

(3)

式中：m_i为粒子质量，t为时间。物质其他相关的物理量，如体积、温度、压力等，均可以通过数值统计方法得到。

1.2 电子力场势函数

在目前的研究中，包含电子力场势函数的模型将原子核处理成点电荷^[20-21]，而将电子处理成满足高斯分布的波包。整个系统的能量是由原子核位置R、电子位置r和相对距离s构成的函数，表达形式如下

$U\left( {R,r,s} \right) = E{}_{{\rm{NN}}}\left( R \right) + {E_{{\rm{Ne}}}}\left( {R,r,s} \right) + {E_{{\rm{ee}}}}\left( {r,s} \right) + {E_{{\rm{KE}}}}\left( {r,s} \right) + {E_{{\rm{pauli}}}}\left( { \uparrow \downarrow ,s} \right)$

(4)

${E_{{\rm{NN}}}} = \frac{1}{{4{\text{π}} {\varepsilon _0}}}\sum\limits_{i < j} {\frac{{{Z_i}{Z_j}}}{{{R_{ij}}}}}$

(5)

${E_{{\rm{Ne}}}} = - \frac{1}{{4{\text{π}}{\varepsilon _0}}}\sum\limits_{i,j} {\frac{{{Z_i}}}{{{R_{ij}}}}} {\rm{erf}}\left( {\frac{{\sqrt 2 {R_{ij}}}}{{{s_j}}}} \right)$

(6)

${E_{{\rm{ee}}}} = \frac{1}{{4{\text{π}}{\varepsilon _0}}}\sum\limits_{i < j} {\frac{1}{{{r_{ij}}}}} {\rm{erf}}\left( {\frac{{\sqrt 2 {r_{ij}}}}{{\sqrt {{s_i}^2 + {s_j}^2} }}} \right)$

(7)

${E_{{\rm{KE}}}} = \frac{{{h^2}}}{{{m_\varepsilon }}}\sum\limits_i {\frac{3}{{2{s_i}^2}}}$

(8)

${E_{{\rm{pauli}}}} = \sum\limits_{{\sigma _i} = {\sigma _j}} {E{{\left( { \uparrow \downarrow } \right)}_{ij}}} + \sum\limits_{{\sigma _i} \ne {\sigma _j}} {E{{\left( { \uparrow \downarrow } \right)}_{ij}}}$

(9)

式中：E_NN表示原子核之间的库仑作用力，E_Ne表示原子核与电子之间的库仑力，E_ee表示电子之间的库仑力，这3项能量与经典分子动力学下的相一致；E_KE表示电子的动能，E_pauli表示Pauli排斥能，这两项能量反映了量子影响； ${\varepsilon _0}$ 和 $\sigma$ 是经典Lennard-Jones势函数下的参量； $Z$ 表示原子核的带电量； ${R_{ij}}$ 和 ${r_{ij}}$ 表示相互作用的距离；m_ε为电子质量，h为普朗克常数。单个电子被处理成满足高斯分布的波包，因而波包的位置x和大小s会随时间变化。按照运动和静止体系，其波函数形式为

$\psi \left( r \right) \propto \exp \left[ { - \left( {\frac{1}{{{s^2}}} - \frac{{2{P_s}}}{s}\frac{i}{h}} \right){{\left( { { r} - { x} } \right)}^2}} \right]\exp \left( {\frac{i}{h} {{{ P}_x}} \cdot { r}} \right)$

(10)

式中：P_x和P_s分别表示高斯波包的平动动量和径向伸缩动量。将（10）式代入时间相关的薛定谔方程，得到简化Hamiltonian的运动方程

$\dot x = \frac{{{\rm{d}}x}}{{{\rm{d}}t}} = \frac{{{P_x}}}{{{m_{\rm{e}}}}}$

(11)

${\dot P_x} = \frac{{{\rm{d}}{P_x}}}{{{\rm{d}}t}} = - \frac{{\partial V}}{{\partial x}}$

(12)

$\dot s = \frac{{{\rm{d}}s}}{{{\rm{d}}t}} = \frac{4}{3}\frac{{{P_{\rm s}}}}{{{m_{\rm{e}}}}}$

(13)

${\dot P_{\rm s}} = \frac{{{\rm{d}}{P_{\rm s}}}}{{{\rm{d}}t}} = - \frac{{\partial V}}{{\partial s}}$

(14)

式中：m_e为电子质量，V为系统势能。这些方程说明电子的位置和大小遵循经典分子动力学，并且独立地随时间变化。因此，eFF方法可以用来模拟长时间、大空间尺度各种极端条件下材料的非绝热动力学过程。

1.3 计算模型

如图1所示，计算模型的大小为123.760 nm×3.536 nm×1.768 nm (x×y×z)，由金属Li和气体H₂组成。Li和H₂之间的界面为光滑平面，无初始扰动。于是，冲击压缩过程可近似成一维流动。e⁻表示电子，±1代表电子的自旋方向。

图 1 计算模型

Figure 1. Computational model

下载: 全尺寸图片幻灯片

界面冲击激波首先由Li区域右侧刚性壁模拟活塞向−x方向以恒速v_p冲击压缩产生。当它经过界面位置时，会形成向左的透射激波和向右的反射激波，电荷分离将在此冲击压缩过程中产生。

2. 结果与讨论

2.1 电荷分离

为了清晰地看到电荷分离情况，选取初始压缩速度v_p=22.50 km/s的工况，如图2（a）所示。可见，强激波冲击运动过程中，激波压缩产生的高温导致电离，电离出自由电子和离子。由于激波压缩区的密度高，电子和离子浓度均高于激波未到达区。根据等离子体动力学，电子和离子将向低浓度区域扩散。由于电子质量极小，碰撞截面小，其扩散速度极高，且未电离原子对其扩散阻力小，因此电子很快跨越激波，扩散到激波未到达区；而离子与原子的质量基本相同，且碰撞截面大，扩散速度低，受其他原子碰撞后扩散阻力大，无法跨越激波。因此，在激波附近必然出现电子-离子分离区。一旦出现电子-离子分离区，等离子体的空间局部中性电场将不复存在，而随着空间电荷的分布出现局部附加电场。附加电场的建立将阻止电子进一步向低电荷密度区扩散，直至达到动态平衡为止，最终形成如图2（a）所示的台阶状电荷分布状态。

图2（b）进一步显示了电荷密度随时间的变化。可见，电荷分离现象伴随整个激波压缩过程，且随着激波向前运动而移动。这是因为运动激波不断压缩导致其所经过区域物质电离，电子和离子不断扩散，建立附加电场，并形成新的扩散平衡。可以确定，电荷分离区将伴随运动激波，以相同的速度运动。

2.2 电荷分离区的强度与宽度

根据图2所描述的电荷一维分布曲线，构建了如图3（a）所示的一维电荷分布模型。它由3部分构成：区域I₁电荷量为负，表示电子分布区域；区域I₃电荷量为正，表示离子区域；中间的区域I₂是从负（电子区域）到正（离子区域）的过渡区，可近似地看成线性分布，它随着激波不断向激波运动方向（x负方向）移动。

图 3 电荷密度分布一维模型以及电荷分离区宽度和强度的时间历程 (v_p=22.50 km/s)

Figure 3. 1D model of charge density distribution and time histories of the width and intensity of the electron-ion separation zone for v_p =22.50 km/s

下载: 全尺寸图片幻灯片

图 2 v_p=22.50 km/s时电荷密度分布一维统计及其随时间的变化

Figure 2. 1D statistic of charge density distribution for v_p=22.50 km/s case and its patterns at different moments

下载: 全尺寸图片幻灯片

过渡区宽度定义为电荷分离区的宽度，过渡区两端的电荷密度差（q₃−q₁）定义为电荷分离区的强度。图3（b）和图3（c）分别显示v_p=22.50 km/s的工况下电荷分离区宽度和强度随时间变化的曲线。可以看出，这两个新定义的量在激波运动过程中基本保持恒定，宽度维持在4.65 nm左右，强度在4.55×10⁸ C/m³附近波动。这是由于：分离区诱导的附加电场抑制了电子和离子的进一步扩散，达到一种局部平衡状态。而扩散动力来自于电离强度，当压缩激波强度一定时，其波后温度和压力一定，电离强度趋于恒定，则扩散动力恒定，因而分离区的强度和宽度必然趋于恒定。

表1进一步统计了其他加载速度下电荷分离区的强度和宽度数据。可见，随着冲击速度的增加，电荷分离区的宽度和强度持续增加。这是因为随着冲击速度的增加，激波压缩强度增加，电离强度增加，分离区宽度和强度自然相应增加。不过从增加的程度来看，宽度增加程度较小，而强度增加程度较高，且均呈现非线性增加特征。

表 1 不同冲击加载速度下分离区的宽度和强度统计

Table 1. Statistics of the width and intensity of the electron-ion separation zone at different impact velocities

v_p/(km·s^–1)	Width/nm	Intensity/(10⁸ C·m^–3)
22.50	4.65	4.55
45.00	6.40	6.11
67.50	8.65	10.25

下载: 导出CSV

| 显示表格

2.3 附加电场对界面附近金属物质的附加加速度模型

2.3.1 附加电场和附加加速度

如前所述，在极端加载条件下，电荷分离会诱导附加电场的产生。本研究根据模拟结果构建了如图3（a）所示的一维电荷分布模型。为验证该模型的合理性，对模型进行电势积分获得一维分布的电场，如图4中红实线所示。

由于电荷分离，激波前后的物质不再呈中性状态，而是分别带部分负电和正电。在附加电场中必然产生加速度，其计算表达式如下

$a = \frac{{{\rho _{\rm{q}}}E}}{{{\rho _{\rm{m}}}}}\quad\quad\quad\quad\quad\quad\quad\quad\quad(15)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad$

(15)

式中：E为附加电场强度，ρ_q和ρ_m分别为单位体积下的物质带电量和物质密度，可通过分子动力学模拟结果的统计得到。

同时根据积分获得的一维分布的电场，计算各处带电物质获得的一维加速度分布，如图4中蓝色实线所示。图4还给出了根据eFF MD模拟结果一维统计获得的电场分布（红色虚线）和物质加速度分布（蓝色虚线），可见二者基本吻合。说明本研究建立的一维电荷分布模型基本合理。

图 4 一维附加电场和附加加速度(t=2.75 ps，v_p=22.50 km/s)的eFF MD统计分析结果和基于1D电荷密度模型的计算结果对比

Figure 4. Comparison of 1D extra electric field and extra acceleration of materials between results of eFF MD statistics and 1D charge density model (t=2.75 ps of v_p =22.50 km/s case)

下载: 全尺寸图片幻灯片

2.3.2 金属/气体界面附近金属物质的加速度模型

由于RM不稳定性的主要演化特征表现为界面物质扰动的振幅增长及物质混合，因此在极端条件下，激波诱导电荷分离导致附加电场的产生，进而引起界面物质的附加加速效应，可能显著影响界面演化。为深入理解极端条件下RM不稳定性的演化特征，特别是附加电场对界面物质的加速效应，对金属/气体界面附近金属侧物质（Li和H₂物质的量之比0.6～0.9，厚度在0～2.0 nm区间）的加速度进行跟踪统计（以激波到达界面位置为统计起点），获得了不同冲击加载速度v_p下界面附近金属物质的附加加速度随时间变化曲线，如图5所示。图5中虚线为根据eFF MD模拟结果统计分析获得，由于电子和离子的动态运动效应，统计值有一定涨幅落差，因此虚线呈现一定的振荡特征。对结果采用Rayleigh模型（（16）式）进行拟合，如图5中红实线所示。综合来看，可以获得以下几点认识。

图 5 不同冲击加载速度v_p下界面金属侧物质附加加速度的变化历程

Figure 5. Acceleration time histories of material located at the metal side of the interface at different v_p

下载: 全尺寸图片幻灯片

（1）界面附近金属物质获得的加速度呈脉冲形式。当激波刚到达界面位置时，电荷分离区跨越气/固界面，界面附近金属侧物质处于分离区中心区域，净电荷接近零；尽管此时分离区电场强度较高，但目标区受到的电场力较弱，加速度较低。随着激波往前推进，电荷分离区亦往前推进，分离区的电荷扩散处于动态平衡，目标区的净电荷逐渐偏正，并迅速向峰值增加；此时目标区受到的电场力开始增强，加速度迅速向峰值挺近，形成脉冲尖峰。随后，由于激波速度大于界面运动速度，随着激波继续向前推进，电荷分离区逐渐远离界面；此时目标区还带有较强的正电荷，但该区电场强度迅速衰减，对应的加速度也迅速衰减，形成完整的加速脉冲效应。

（2）随着冲击加载速度的增加（从22.50 km/s增加至78.75 km/s），目标区域的峰值加速度迅速提高（从1.7×10¹⁷ m/s²增加至1.2×10¹⁸ m/s²），同时脉冲持续时间逐渐减小（从1.5 ps降至1.0 ps）。这是因为随着冲击加载速度的提高，压缩程度增强，电离强度增加，电荷分离区强度相应提高，导致峰值加速度提升；但同时激波运动速度加快，激波跨越界面的时间相应缩短，脉冲持续时间自然有所降低。

（3）需要指出的是，曲线后期呈现较大幅度的振荡，这是因为激波到达左侧边界发生了Re-shock现象，Re-shock导致电荷分离呈现较大幅度的振荡。不过这一过程不是本研究关注内容，在此不作讨论。

本研究拟合所采用的模型公式为

$a\left( t \right) = {a_{\max }}{{\rm{e}}^{{1 / 2}}}\frac{t}{{{\tau _{{\rm{rise}}}}}}\exp \left( { - \frac{{{t^2}}}{{2{\tau _{{\rm{rise}}}}}}} \right)$

(16)

式中：a_max和τ_rise分别代表加速度峰值和上升时间，拟合参数如表2所示，拟合曲线如图6所示。

表 2 不同加载速度下Rayleigh模型的拟合参数

Table 2. The fitted parameters of Rayleigh model at different v_p

v_p/(km·s^–1)	a_max/(10¹⁷ m·s^–2)	τ_rise/ps
22.50	1.7	0.50
33.75	1.9	0.38
45.00	3.5	0.30
56.25	5.3	0.24
67.50	7.5	0.20
78.75	12.0	0.17

下载: 导出CSV

| 显示表格

图 6 Rayleigh模型的a_max和τ_rise与v_p的拟合关系

Figure 6. Fitted correlations of a_max and τ_rise with v_p in Rayleigh model

下载: 全尺寸图片幻灯片

从图6中可以明显看出a_max和τ_rise分别随初始加载速度呈指数上升和指数下降趋势。这反映了拟合模型的合理性，即峰值加速度及脉冲上升时间与冲击强度的非线性关系。进一步对两个参数进行拟合，得到如下表达式

${a_{\max }} = 6.233\;20 \times {10^{16}} \times {{\rm{e}}^{0.037\; 44{v_{\rm{p}}}}}\quad\quad\quad(17)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\;\quad\quad\quad\quad\quad$

(17)

${\tau _{{\rm{rise}}}} = 0.773\;69 \times {{\rm{e}}^{ - 0.020\;35{v_{\rm{p}}}}}\quad\quad\quad\quad(18)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\;\quad\quad\quad\;\;\;$

(18)

将（17）式和（18）式代入(16)式，可以获得界面附近金属侧物质受到激波压缩电离诱导形成附加电场的加速度经验公式

$a\left( {t,{v_{\rm{p}}}} \right) = 6.233\;20 \times {10^{16}} \times {{\rm{e}}^{0.037\;44{v_{\rm{p}}} + 0.5}} \times \frac{t}{{0.773\;69 \times {{\rm{e}}^{ - 0.020\;35{v_{\rm{p}}}}}}} \times \exp \left( { - \frac{{{t^2}}}{{1.547\;38 \times {{\rm{e}}^{ - 0.020\;35{v_{\rm{p}}}}}}}} \right)$

(19)

这是在本计算模型条件下冲击加载速度在20～80 km/s范围内激波压缩电离诱导附加电场引起的界面金属侧物质附加加速度的经验模型，揭示了极端加载条件下物质电离引起的附加电场对界面物质运动影响的基本特征。从积分效应来看，虽然其脉冲时间在皮秒量级，但所导致的界面物质的附加速度增量可达百千米每秒量级，显然不可忽视。当然，这只是界面物质加速的一方面，其他作用力耦合后，最终效应取决于综合效果。

3. 结　论

采用基于电子力场势函数的分子动力学方法，模拟了极端加载条件下金属/气体(Li/H₂)平滑界面运动过程，对激波压缩电离过程及电荷扩散在激波附近诱导的电荷分离效应进行了理论分析和结果统计，发现了电荷分离诱导附加电场及其对界面物质加速作用的现象，进一步开展了理论建模和经验拟合分析，获得了电离引起界面附近金属侧物质附加加速度的经验模型，主要获得如下结论：

（1）极端冲击加载下，激波诱导附加电场的基本机理是，极端冲击使材料电离出较轻的电子和较重的离子，电子和离子的扩散差异导致在激波附近形成电荷分离区并达到动态扩散平衡，诱导出附加电场；

（2）统计分析发现，电荷分离区的宽度和强度在激波运动过程中基本保持稳定，它们仅与初始冲击加载速度相关，即与激波强度相关，呈非线性关系；

（3）由附加电场引起的界面物质的附加加速效应可能显著影响界面演化，统计分析发现，界面处金属侧的加速度满足Rayleigh模型，呈脉冲形态，揭示了极端加载条件下物质电离引起的附加电场对界面物质运动影响的基本特征，从积分效应来看，虽然其脉冲时间在皮秒量级，但其影响却不可忽视。

需要指出的是，本研究结果基于微观分子动力学模拟，其尺度在百纳米量级，相关结论尽管揭示了内在机理，但在一定程度上还受尺度效应的限制。后续工作将进一步加大模拟尺度并开展理论建模，进而推广到宏观尺度范围，使之更贴合工程实际工况。

图 1 Schardin问题的几何模型

Figure 1. Geometric model of Schardin problem

下载: 全尺寸图片幻灯片

图 2 实验结果与仿真结果对比

Figure 2. Comparison of experimental results and simulation results

下载: 全尺寸图片幻灯片

图 3 设置扰流板的坑道模型

Figure 3. Tunnel model with spoilers

下载: 全尺寸图片幻灯片

图 4 坑道Section Ⅰ段的俯视图

Figure 4. Top view of tunnel Section Ⅰ

下载: 全尺寸图片幻灯片

图 5 6.5 ms时不同厚度扰流板附近的压力云图

Figure 5. Pressure nephogram near spoilers with different thickness at 6.5 ms

下载: 全尺寸图片幻灯片

图 6 波阵面峰值超压随距离变化的曲线

Figure 6. Curves of peak overpressure versus distance

下载: 全尺寸图片幻灯片

图 7 设置倾斜扰流板的坑道模型

Figure 7. Tunnel model with inclined spoiler

下载: 全尺寸图片幻灯片

图 8 单元A₁和单元B₁处超压比(a)和冲量比(b)与倾角的关系

Figure 8. Overpressure ratio (a) and impulse ratio (b) vesus the inclination angle at position A₁ and B₁

下载: 全尺寸图片幻灯片

图 9 单元C₁处超压比(a)和冲量比(b)与倾角的关系

Figure 9. Overpressure ratio (a) and impulse ratio (b) vesus the inclination angle at position C₁

下载: 全尺寸图片幻灯片

图 10 α=60°和α=120°扰流板附近的压力云图

Figure 10. Pressure nephogram of spoilers with α=60° and α=120°

下载: 全尺寸图片幻灯片

图 11 不同时刻坑道内压力流场

Figure 11. Pressure flow field in tunnel at different times

下载: 全尺寸图片幻灯片

图 12 单元超压随扰流板间距和距离变化的三维曲线

Figure 12. Three dimensional curve of over pressure changing with interval and distance

下载: 全尺寸图片幻灯片

表 1 空气的材料参数

Table 1. Material parameters of air

ρ/(kg·m^-3)	C₁	C₂	C₃	C₄	C₅	C₆	E₀/MPa	V₀
1.29	0	0	0	0.4	0.4	0	0.25	1.0

下载: 导出CSV

表 2 三角楔的材料参数

Table 2. Material parameters of triangular wedge

ρ/(g·cm^-3)	G/MPa	E/MPa	ν	A_J-C/MPa	B_J-C/MPa	C_J-C	N	T_m/K
7.83	0.30	0	0.25	496	434	0.014	0.26	1 788

下载: 导出CSV

表 3 壁面的材料参数

Table 3. Material parameters of the wall surface

ρ/(g·cm^-3)	E/GPa	ν	R_m/MPa	σ_c/MPa	η/(Pa·s)	σ_s/MPa
2.70	25	0.30	3.10	14.48	0	28.96

下载: 导出CSV

表 4 不同扰流板厚度下单元的峰值超压

Table 4. Peak overpressure of units with different spoiler thickness

d_s/cm	p_A₁/MPa	p_B₁/MPa	p_C₁/MPa
10	0.457	0.323	0.597
20	0.451	0.320	0.588
30	0.442	0.317	0.581
40	0.434	0.315	0.577
50	0.425	0.313	0.570
60	0.421	0.310	0.564
70	0.416	0.309	0.562
80	0.413	0.307	0.558

下载: 导出CSV

参考文献(12)

[1]	于文华, 张亚栋.爆炸冲击波在坑道内传播规律研究[J].爆破器材, 2013, 42(3):1-4. http://www.cnki.com.cn/Article/CJFDTotal-BPQC201303002.htm YU W H, ZHANG Y D.Study on propagation laws of explosion shock wave in tunnels[J]. Explosive Materials, 2013, 42(3):1-4. http://www.cnki.com.cn/Article/CJFDTotal-BPQC201303002.htm
[2]	刘晶波, 闫秋实, 伍俊.坑道内爆炸冲击波传播规律的研究[J].振动与冲击, 2013, 28(6):8-11. http://d.wanfangdata.com.cn/Periodical_zdycj200906003.aspx LIU J B, YAN Q S, WU J.Analysis of blast wave propagation inside tunnels[J]. Journal of Vibration and Shock, 2013, 28(6):8-11. http://d.wanfangdata.com.cn/Periodical_zdycj200906003.aspx
[3]	覃彬, 张奇, 向聪, 等.数值模拟研究分叉巷道中冲击波传播规律[J].含能材料, 2008, 16(6):741-744. http://d.wanfangdata.com.cn/Periodical_hncl200806024.aspx QIN B, ZHANG Q, XIANG C, et al.Numerical simulation research on shock wave propagation characteristics inside branch tunnel[J]. Chinese Journal of Energetic Materials, 2008, 16(6):741-744. http://d.wanfangdata.com.cn/Periodical_hncl200806024.aspx
[4]	王启睿, 张晓忠, 孔福利, 等.多级穿廊结构坑道口部内爆炸冲击波传播规律的实验研究[J].爆炸与冲击, 2011, 31(5):449-454. http://www.cnki.com.cn/Article/CJFDTOTAL-BZCJ201105002.htm WANG Q R, ZHANG X Z, KONG F L, et al.Shock wave propogation of an explosion inside its entrance to a multilevel gallery tunnel[J]. Explosion and Shock Wave, 2011, 31(5):449-454. http://www.cnki.com.cn/Article/CJFDTOTAL-BZCJ201105002.htm
[5]	杨科之, 周布奎, 王安宝, 等.坑道中扰流板对化爆冲击波的衰减作用[J].防护工程, 2010, 32(5):58-61. http://www.cqvip.com/QK/96961X/200802/27562137.html YANG K Z, ZHOU B K, WANG A B, et al.Attenuation effect of the interference plate on explosion shock wave inside a tunnel[J]. Protective Engineering, 2010, 32(5):58-61. http://www.cqvip.com/QK/96961X/200802/27562137.html
[6]	SHA S, CHEN Z H, JIANG X H.Influence of obstacle geometries on shock wave attenuation[J]. Shock Waves, 2014, 24:573-582. doi: 10.1007/s00193-014-0520-9
[7]	TAKAYAMA K. Attenuation of shock waves propagating over arrayed spheres[C]//Proceedings of 24th International Congress High-Speed Photography and Photonic. Sendai, Japan, 2000: 582-588.
[8]	BERGER S, SADOT O, BEN-DOR G.Experimental investigation on the shock wave load attenuation by geometrical means[J]. Shock Waves, 2010, 20:29-40. doi: 10.1007/s00193-009-0237-3
[9]	朱建, 王曙光, 陆伟东.坑道内挡板对化爆冲击波消减作用的数值模拟[J].爆破, 2008, 25(2):26-29. http://www.cqvip.com/QK/96961X/200802/27562137.html ZHU J, WANG S G, LU W D.Numerical simulation of the pressure attenuation of chemical shock wave in tunnel with baffles[J]. Blasting, 2008, 25(2):26-29. http://www.cqvip.com/QK/96961X/200802/27562137.html
[10]	沙莎, 陈志华, 张焕好, 等.Schardin问题的数值研究[J].物理学报, 2012, 61(6):1-8. http://d.wanfangdata.com.cn/Periodical_wlxb201206055.aspx SHA S, CHEN Z H, ZHANG H H, et al.Numerical investigations on the Schardin's problem[J]. Acta Phyisca Sinica, 2012, 61(6):1-8. http://d.wanfangdata.com.cn/Periodical_wlxb201206055.aspx
[11]	Livemore Software Technology Corporation.LS-DYNA keyword use's manual[M]. California:Livemore Software Technology Corporation, 2003.
[12]	李祥春, 聂百胜, 杨春丽, 等.封闭空间内瓦斯浓度对瓦斯爆炸反应动力学特性的影响[J].高压物理学报, 2017, 31(2):135-144. doi: 10.11858/gywlxb.2017.02.005 LI X C, NIE B S, YANG C L, et al.Effect of gas concentration on kinetic characteristics of gas explosion in confined space[J]. Chinese Journal of High Pressure Physics, 2017, 31(2):135-144. doi: 10.11858/gywlxb.2017.02.005

施引文献

资源附件(0)

访问统计

点击查看大图

图(12) / 表(4)

计量

文章访问数: 6646
HTML全文浏览量: 2874
PDF下载量: 246

1. 方法与模型
1.1 MD方法
1.2 电子力场势函数
1.3 计算模型
2. 结果与讨论
2.1 电荷分离
2.2 电荷分离区的强度与宽度
2.3 附加电场对界面附近金属物质的附加加速度模型
3. 结　论

扰流板结构参数对冲击波衰减特性影响的仿真分析

doi: 10.11858/gywlxb.20170585

作者简介: 赵蓓蕾(1991—), 女，博士研究生，主要从事航天任务分析与设计研究.E-mail:1178685385@qq.com

计量

出版历程

Simulation Analysis of Influence of Spoiler Structural Parameters on Shock Wave Attenuation Characteristics

1. 方法与模型

1.1 MD方法

1.2 电子力场势函数

1.3 计算模型

2. 结果与讨论

2.1 电荷分离

2.2 电荷分离区的强度与宽度

2.3 附加电场对界面附近金属物质的附加加速度模型

2.3.1 附加电场和附加加速度

2.3.2 金属/气体界面附近金属物质的加速度模型

3. 结 论

计量

出版历程

目录

1. 方法与模型

1.1 MD方法

1.2 电子力场势函数

1.3 计算模型

2. 结果与讨论

2.1 电荷分离

2.2 电荷分离区的强度与宽度

2.3 附加电场对界面附近金属物质的附加加速度模型

3. 结 论

作者简介:
赵蓓蕾(1991—), 女，博士研究生，主要从事航天任务分析与设计研究.E-mail:1178685385@qq.com

3. 结　论

3. 结　论