黄土和砂土岩中填实爆炸辐射弹性波的对比研究*

2019-06-21王占江张景森郭志昀

爆炸与冲击 2019年5期

卢强，王占江，张景森，丁洋，李进，郭志昀

（1.西北核技术研究所强动载与效应重点实验室，陕西西安 710024；

2.西安交通大学航天航空学院机械结构强度与振动国家重点实验室，陕西西安 710049）

在地下爆炸研究中，爆源附近的介质受到高温、高压作用会发生极大的不可逆变形。随着爆炸应力波的传播衰减，距爆心一定距离r0处介质应力-应变的状态演化为弹性或黏弹性状态，r0称为弹性半径。由弹性半径处向外辐射的弹性波（黏弹性波）源，其特性与爆炸当量、介质特性及炸药与介质的耦合形式相关[1]。地下爆炸的弹性波源特性是评估地下爆炸远区震动效应的重要参数，由震源辐射出的爆炸能量及弹性应力波峰值、频率等特征共同表征了其特性。

为研究地下爆炸辐射弹性波源的特性，Antoun 等[2]提出了一种考虑介质膨胀、屈服、损伤、孔隙压实等效应的热力学相容本构关系，以PILE DRIVER 地下核爆炸实测自由场粒子速度数据[3]为基础给出了该场地花岗岩介质本构模型的相关参数，并将该模型应用于计算该爆炸场地大当量爆炸弹塑性区波的传播，以此确定地下爆炸辐射弹性应力波的特征。周钟等[4-5]则基于连续介质力学和不相融混合物理论，假定组分间无相对运动，在屈服面中引入孔隙影响因子，提出了一种多孔含水介质流固耦合的本构模型，并给出了孔隙的演化方程，对水饱和花岗岩介质中爆炸应力波的传播进行了数值模拟，较系统地研究了含水花岗岩中地下爆炸应力波的传播规律及地下爆炸震源函数的特征。刘文韬等[6]采用含损伤的弹塑性本构模型，对封闭式地下核爆炸波的传播规律展开了研究，针对PILE DRIVER 地下核试验[3]进行了数值模拟，计算的压力和速度峰值衰减、质点位移和速度波形以及爆炸空腔都与其他模拟结果比较接近。上述研究是从爆炸能量与源区介质耦合机理的角度出发，结合爆炸近区实测的介质变形、粒子运动或应力历程等，研究高压下源区介质的状态方程及本构关系，并通过数值方法计算源区介质的弹性边界，以及计算弹性边界处的应力或粒子速度时间历程，以此确定地下爆炸辐射弹性波源特征。

另一方面，研究地下爆炸辐射弹性波源的特性也可从弹性波或黏弹性波传播的角度出发，直接测量弹性区的介质变形、粒子运动或应力历程等，结合弹性波或黏弹性波传播理论，构建弹性波源模型。Larson[7]利用在弹性区传播的球面波研究了填实和空腔爆炸条件下爆炸能量耦合特性，Rodean[8]则基于理想弹性假设研究了由球形爆炸源辐射的弹性波频率、波长及衰减特性。赖华伟等[9-10]基于特征线方法分析了线黏弹性球面应力波的传播规律，Wang 等[11]基于特征线方法研究了非线性黏弹性球面发散应力波的传播演化，给出了强间断球面波的衰减规律，并提出一种反演介质黏弹性参数的方法。卢强等基于黏弹性假设给出了表征黏弹性球面波震源特性的折合位移势、折合速度势等的Laplace 解[12]，并以实测粒子速度为基础分析了爆炸辐射弹性波震动和传播特性[13]。

本文采用微米级太安颗粒整体压装的微型炸药球（0.125 g TNT 当量、直径5 mm）作为爆炸源，分别把重塑黄土和砂土岩介质作为源区介质，利用西北核技术研究所建立的塑性源区可置换的 ∅1 370 mm×1 200 mm 黄土化爆实验平台来接受源区爆炸辐射的弹性波[14]，通过预埋在黄土化爆实验平台内部的圆环型电磁粒子速度计来测量不同距离处的粒子速度[15]，分析黄土和砂土岩分别作为源区介质时填实爆炸激发的弹性波源变化及应力波传播特征的演化，建立用于评估地质材料（特别是对于不易加工为大尺寸模型的材料）的地下爆炸弹性波能量耦合强度实验方法。

1 实验简介

球面波加载实验的典型装置如图1 所示。实验过程为：（1）对雷管进行加电起爆；（2）雷管引爆装药量为0.2 g RDX/m 的 ∅1 mm 铅皮柔爆索，柔爆索传爆至 ∅5 mm 微型炸药球中心并引爆炸药球；（3）微型炸药球激发出在介质内部传播的爆炸应力波，预埋在介质内部的粒子速度计记录径向粒子速度时间历程。

由于柔爆索传爆过程中爆炸能量会向介质中辐射，对微型炸药爆炸激发的应力波造成干扰，同时柔爆索爆炸后飞溅的铅颗粒在磁场中运动会产生强电磁干扰。为抑制柔爆索爆炸产生的强电磁和力学干扰，王占江[16]设计了双层薄壁不锈钢管防护组件作为柔爆索的传爆通道，这样既有效地屏蔽住柔爆索飞散造成的电磁干扰，又减弱柔爆索爆炸产生的引力波。在上述干扰抑制措施下，线装药为0.5 g RDX/m的柔爆索爆炸干扰信号对应的粒子速度和位移峰值比0.125 g TNT 主炸药爆炸的相应值低至少一个量级[16]（正式实验时起爆炸药球用柔爆索线装药仅为0.2 g RDX/m，干扰更小）。对岩土介质来讲，柔爆索干扰被抑制到这样的程度，已处于实验介质差别和系统环境等限制造成的实验不确定度以内。

图 1 塑性区可置换的黄土实验平台Fig.1 Experimental apparatus for loess with replaceable model for plastic zone

本文采用的黄土化爆实验平台为重复使用的平台（如图2(a)所示），总质量约3.2 t，体积约为1.8 m3，预埋入黄土介质内的圆环型电磁粒子速度计（在样品中心面上，和炸药球同心）共计19 个，安置于90～640 mm 的半径范围内（0.125 g TNT 爆炸下对应比距离 R= r/Q1/3为180～1 280 mm/g1/3，Q 为爆炸当量）。可置换塑性区的重塑黄土(loess，记为L)和砂土岩(rock-like sandy soil，记为RS)样品尺寸为∅124 mm×124 mm（见图2(b)和图2(c)），样品的组分构成见表1，黄土和砂土岩的密度分别为1.80 和1.75 g/cm3，纵波波速分别为1.1 和1.2 km/s，两种材料的波阻抗之比约为0.94，近似相等。每种样品进行两次实验，分别记为L1, L2 和RS1, RS2。0.125 g TNT 炸药球加载下黄土和砂土岩中塑性区半径约为40 mm，因此爆炸应力波在可置换塑性区模型内已形成弹性波源，在可置换塑性区以外的测试区域承受的是弹性加载[14]。

图 2 黄土实验平台及源区介质模型Fig.2 The loess platform and the sample model in source area

表 1 黄土和砂土岩样品各成分组成(单位：%)Table 1 Components of loess and rock-like sandy soil (unit: %)

2 爆炸辐射弹性波时域特征分析

图3 给出了黄土作为源区介质时0.125 g TNT 炸药球填实爆炸加载下比距离180～1 280 mm/g1/3范围内实测的19 个典型粒子速度(vr)曲线。可以看出，粒子速度峰值随着传播距离的增加而逐渐减小，在靠近样品边界的几个粒子速度计测得的粒子速度波形会受到样品边界反射波的影响（图3(d)）。以砂土岩作为源区介质时粒子速度波形和图3 类似，此处不再给出。

图 3 黄土作为源区介质的填实爆炸实测粒子速度(vr)Fig.3 Measured particle velocity in a tamped explosion with taking loess as source medium (vr)

图4 和图5 分别给出了以黄土、砂土岩作为源区介质时填实爆炸条件下粒子速度峰值(vmax)和粒子位移峰值(umax)随比距离(R)的变化曲线。可以看出，砂土岩粒子速度峰值和粒子位移峰值整体比黄土要大，这说明填实爆炸条件下砂土岩耦合的向外传播的能量比黄土耦合的能量要大。

图 4 粒子速度峰值(vmax)与比距离(R)的关系Fig.4 Maximum particle velocity (vmax) varying with specific distance (R)

图 5 粒子位移峰值(umax)随比距离(R)的关系Fig.5 Maximum particle displacement (umax) varying with specific distance (R)

工程上通常采用幂指数拟合波峰值衰减规律的方法来表征波幅的衰减规律，粒子速度峰值vmax(R)的衰减满足

式中：A 为拟合系数；ζv为衰减指数，其值越大则衰减越快，可由相邻位置的粒子速度峰值确定

进一步利用图4 给出的粒子速度峰值数据，结合式（2）给出了反映粒子速度峰值变化快慢的衰减指数ζv随比距离的变化，如图6 所示。可以看出，由实验数据计算得到的衰减指数ζv具有一定的离散性，但从统计角度看，在180～1 280 mm/g1/3比距离模拟范围内应力波粒子速度峰值的衰减指数ζv表现为随比距离增加而增加的特点，这说明随着波传播距离的增加，粒子速度峰值的衰减在此区域内有逐渐变快的趋势。把ζv进行线性拟合，可以给出ζv和比距离R 满足的关系：

式中：R 的单位为mm/g1/3。ζv在180～1 280 mm/g1/3的比距离模拟范围衰减指数从1.25 增加到2.03，可见ζv不仅包含了球面波的几何发散，同时也包含了介质的黏性衰减。

同式（1）定义类似，图7 给出了位移峰值的衰减指数ζu随比距离的变化，把ζu进行线性拟合，可以给出其和比距离R 满足的关系：

图 6 速度衰减指数(ζv)随比距离(R)的变化Fig.6 Damped exponential of velocity (ζv) varying with specific distance (R)

图 7 位移衰减指数(ζu)随比距离(R)的变化Fig.7 Damped exponential of displacement (ζu) varying with specific distance (R)

式中：R 的单位为mm/g1/3。ζu在180～1 280 mm/g1/3比距离模拟范围衰减指数从1.18 增加到1.76。从图6 和图7 给出的ζv和ζu的均值曲线看，以砂土岩和黄土作为源区介质条件时，两种介质中爆炸激发的弹性波的衰减规律基本一致，这反映了作为波传播载体的米级黄土样品平台对波传播衰减的主导作用。

把粒子速度波形的第一个正向速度区间对应的时间区间定义为正向脉宽ΔT+，把粒子速度峰值一半位置对应的时间区间定义为半高宽ΔT1/2。图8 和图9 分别给出了以砂土岩和黄土作为源区介质时粒子速度波形的正向脉宽ΔT+和半高宽ΔT1/2随比距离的变化曲线，可以看出ΔT+和ΔT1/2的实验结果的具有一定离散性。从图8 可知，砂土岩的正向脉宽ΔT+和黄土没有显著的差异。从图9 可以看出，砂土岩的半高宽ΔT1/2比黄土低。

图 8 正向脉宽(ΔT+)随比距离(R)的变化Fig.8 Width of positive speed pulse (ΔT+) varying with specific distance (R)

图 9 半高宽(ΔT1/2)随比距离(R)的变化Fig.9 Full width at half maximum (ΔT1/2) varying with specific distance (R)

3 爆炸辐射弹性波频域特征分析

下面从频域角度分析一下粒子速度频率 f随着传播距离的变化。由于1 040～1 280 mm/g1/3区域内实测粒子速度受样品边界反射波的影响，本文仅考虑180～960 mm/g1/3范围内粒子速度波形的频率变化。图10 给出了黄土作为源区介质时粒子速度（对应图3）的振幅谱，可以看出粒子速度振幅谱低频部分的值随着传播距离的增加而减小。图11 给出了以砂土岩和黄土作为源区介质时填实爆炸辐射弹性波粒子速度波形的主频 f0（振幅谱峰值处对应的频率）随传播距离的变化。

整体上看，砂土岩和黄土中填实爆炸辐射弹性波的主频随着传播距离的增加均具有先增加而后减少的现象。由主频统计均值曲线看，在比距离180～280 mm/g1/3范围内主频随传播距离增加逐渐增加，在比距离高于280 mm/g1/3时主频 f0随传播距离增加逐渐减少。另外，砂土岩作为源区介质时粒子速度主频变化曲线在黄土的下方，这说明在砂土岩中填实爆炸辐射的弹性应力波频率较低。按照爆炸相似率，球形爆炸激发的波主频 f0与 Q1/3成反比关系，即

图 10 黄土作为源区介质时粒子速度的振幅谱Fig.10 Amplitude spectrum of the particle velocity with taking the loess as source medium

图 11 主频(f0)随比距离(R)的变化Fig.11 Scaled dominant frequency (f0) varying with specific distance (R)

若假定砂土岩和黄土材料性质相同，则可按照式（5）给出砂土岩和黄土作为源区介质时的相似关系：

由于砂土岩中辐射的弹性应力波主频比黄土低，按照式（6），砂土岩中产生弹性波源的等效爆炸当量比黄土要高。实际上两种介质中填实爆炸实验采用的炸药球均是0.125 g TNT，且两种材料的性质并不相同。这说明在砂土岩中爆炸耦合的向外传播的弹性波能量比黄土介质要大，此结论与砂土岩中爆炸激发的弹性波粒子速度和位移峰值比黄土大的结论一致。