Evolution of the radiated seismic wave energy of underground explosion in visco-elastic solids
-
摘要: 地下爆炸与介质的能量耦合和介质中的波传播机制是理解地下爆炸源物理的重要基础。为研究地下爆炸辐射地震波能量的传播衰减规律，分析了黏弹性介质中地下爆炸地震波能量的组成。基于无限介质中黏弹性球面波理论，给出了速度、位移、应力、应变等物理量Laplace域的理论解。利用Laplace数值逆求解方法，建立了黏弹性介质中地下爆炸辐射地震波场的计算方法。以干黄土作为典型黏弹性材料，计算给出了地震波能量的传播特征，分析了地下爆炸辐射能量的传播衰减规律。结果表明：（1）在黏弹性介质中，某球面处流入的能量随半径增加而逐渐降低。在理想弹性介质中，某球面处流入的能量在几倍弹性半径外即可稳定到某一定值；（2）在某一固定的有限观测区域内，当观测时间足够长时，势能和耗散能均趋于某一定值，辐射动能趋于零；（3）当有限的观测区域能容纳一个完整波长的地震波时，地震波辐射动能的稳态值随波传播距离的增大而减小，总体上可以用指数函数和幂函数进行分段拟合。Abstract: The energy coupling between the underground explosion and the medium and the wave propagation mechanism in the medium are important bases for understanding the physics of the underground explosion source. In order to study the law of the propagation and attenuation for the seismic wave energy of underground explosion, the composition of the radiated energy of underground explosion in viscoelastic medium was analyzed, and the formulas for calculating inflow energy, outflow energy, radiated kinetic energy, potential energy and dissipation energy in a limited observation region were given. Based on the theory of viscoelastic spherical wave in infinite medium, the theoretical solutions of velocity, displacement, stress and strain in the Laplace domain were given by using the exponential attenuation pressure model and the generalized Maxwell viscoelastic model. The numerical solutions of velocity, displacement, stress and strain were given by using the Laplace numerical inverse method, and the inflow energy, outflow energy, radiated kinetic energy, potential energy and dissipation energy were calculated by these numerical results. The numerical results of different components of seismic wave energy are consistent with the theoretical results, and the correctness of this method is proved. Using the dry loess as typical viscoelastic material, the radial stress and particle velocity at different radii were calculated, and the relationship between the inflow energy at different radii and the wave propagation distance was obtained. The spatial distributions of the radiated kinetic energy of seismic wave were calculated by using the spatial distributions of the radial particle velocity at different times, and the propagation law of the radiated kinetic energy was obtained. The changes of the inflow energy and the radiated kinetic energy with the propagation distance in the limited observation area were analyzed, and the results show as follows: (1) In a viscoelastic medium, the energy flowing into a sphere surface decreases gradually with the increase of radius. In an ideal elastic medium, the energy flowing into a sphere surface at the elastic radius of about several times can be stabilized to a constant value. (2) The potential energy and the dissipative energy tend to constant values when the observation time is long enough in a fixed limited observation region, and the radiated kinetic energy tends to zero. (3) When a limited observation area can hold the seismic waves with complete wave length, the steady-state value of the radiated kinetic energy of seismic waves decreases with the increase of the wave propagation distance. In general, exponential function and power function can be used for piecewise fitting of the attenuation law for the steady-state value of the radiated kinetic energy of the seismic wave.
-
Key words:
- underground explosion /
- spherical wave /
- viscoelasticity /
- seismic wave energy /
- radiated energy
-
材料和结构在经受爆炸等冲击载荷时，压力-冲量（p-I）曲线是评价其损伤程度的重要指标。起初，p-I曲线被用来评价建筑物受爆炸载荷时的破坏程度和人员的伤亡程度；近年来，它逐渐被用来评价结构和材料的损伤程度[1-2]。
美国对不同截面形状和材料的梁进行了大量的爆炸毁伤实验[3]，得到了满足一定关系的p-I曲线。汪明[1]、Zhang等[4]分别对不同爆炸载荷下的钢柱和夹层玻璃窗进行了研究，得到了相应的p-I曲线，其爆炸损伤参数与梁结构的爆炸损伤参数不同。p-I曲线只有压力和时间2个量纲，在相同载荷、相同材料结构工况下，可认为爆炸损伤参数与时间无关，从而通过微积分理论将其转化成压力-时间公式。
自19世纪以来，尤其是二战时期，为了满足海洋军事的需要，科学家们对水下爆炸问题进行了大量的研究[5-6]。在军事上，水下爆炸研究大多是关于炸药爆炸对舰船、潜艇等结构的损伤效应[6]。水下爆炸的损伤效应主要体现在冲击波和气泡脉动两方面[7]，而在冲击波的理论研究中，压力衰减规律一直是研究的重点内容。
关于炸药水下爆炸冲击波的压力分布及其时间变化规律，诸多学者基于冲击波传播理论和实验研究给出大量的模型。Cole[8]基于冲击波传播理论推导出了冲击波的压力衰减指数形式；Орленко[9]基于Cole的指数衰减公式，给出了冲击波压力在衰减后半段的反比形式。Cole和Орленко理论都是冲击波衰减规律的常用理论，前者适用于压力衰减前半段，后者适用于压力衰减后半段。Kiciński等[10]总结了诸多学者的冲击波压力衰减模型，提出当炸药的尺寸以及炸药所处深度不同时，经验公式的参数值会有所不同。Stiepanow对Cole公式做了进一步研究，得到了可模拟TNT炸药水下爆炸冲击波压力衰减全过程的经验公式，包括指数衰减段、倒数衰减段和倒数衰减后段3个部分[11]。Keil[12]提出了冲击波的压力衰减表达式，并结合实验研究给出了经验参数。Ming等[13]、Reid[14]、Rajendran等[15]在Keil[12]的研究基础上，对经验参数进行了修正。Geers等[16]为了提高精度、减小模型与实验之间的误差，提出了一种双指数衰减模型，但其形式复杂，参数较多。
前面诸多研究主要基于冲击波传播理论，结合实验研究得到经验公式，对水下爆炸压力的衰减特性进行研究，然而大多忽视了爆炸的损伤效应，较少从爆炸损伤的角度分析压力的衰减特性。本文中，基于材料的冲击波损伤机理，推导一种新的压力衰减规律的数学表达式，采用新表达式对实验数据进行拟合，比较拟合、Cole理论、Орленко理论和实验得到的比冲量以及比冲击波能，分析新公式的计算精度。
1. 理论推导
在材料的损伤效应研究中，采用p-I图研究材料爆炸损伤的重要参数，材料的p-I公式[2-4]可表达为：
(p−a)(I−b)=c (1) 式中：a为准静态载荷作用下材料达到某一损伤程度的临界压力，b为动态载荷作用下材料达到损伤破坏的临界冲量，c 为与损伤程度和材料相关的常数[17-19]。
将式(1)应用在水下爆炸研究中时，由于水的不可压缩性，可认为a、b、c均为与时间无关的常数。水承受负压时，其稳定性会发生变化，当负压达到某个阈值时，其内部会出现蒸汽泡，并发生空化[20-21]。水受到冲击载荷时，其状态会发生短暂改变，之后又很快恢复。因此，对于水来说，a为准静态载荷作用下水达到与冲击载荷作用后相同状态的临界压力，通常为负值，数值的大小与药量、距离和水变化后的状态相关；爆炸载荷为强冲击载荷，b近似等于冲击波冲量。
在水下爆炸的参数计算中，I可表示为：
I=∫pdt (2) 式中：t为时间。
将式(1)等号两边同时对时间t进行微分，结合式(2)可得：
p=−c(p−a)2d(p−a)dt (3) 将式(3)两边同时乘以dt/p，并对t积分可得：
t=ca2[ln(1−ap)+1p/a−1]+C (4) 式中：C为积分常数。
将式(4)转化成无量纲量的形式：
k=ta2c=f(p/a)=ln(1−ap)+1p/a−1 (5) 在Cole和Орленко的理论中，冲击波压力衰减表达式为：
p={pme−tθ0≤t＜θpm0.368θtθ≤t≤(5∼10)θ (6) 式中：pm为冲击波峰值压力，
θ 为压力从pm下降到pm/e时所需要的时间。通常pm≫ a，当p=pm时，式(4)中ln(1−a/p)+1/(p/a−1)≈0 ，t ≈ 0，因此在后续拟合中，可认为C=0。比较式(5)和式(6)可以看出，式(5)的第1项为压力关于时间的指数表达式，第2项为反比表达式。 式(5)结合了Орленко和Cole理论的优点，在描述冲击波压力衰减规律时具有更高的精度。
2. 水下爆炸实验
2.1 实验材料
进行水下爆炸实验所需的材料和设备包括：黑索今炸药（RDX）、雷管、PCB压力传感器、恒流源、示波器、起爆器。
2.2 实验装置
基于水下爆炸罐开展水下爆炸实验，水下爆炸罐是一个上端开口的圆柱体，直径5.0 m，壁厚5 mm，高5.0 m。药包和PCB传感器置于水面以下3.0 m处，传感器与药包的距离r取1.0或1.5 m，实验装置见图1。
采用式(4)拟合实验结果时，p代表入射波的压力，不需要添加结构，测出冲击波在水中的压力时程曲线即可。水下爆炸实验采用的炸药分别为9 g RDX、12 g RDX和19 g RDX。
2.3 实验结果
水下爆炸实验采集的数据受噪声、杂波、起爆信号等干扰，为此，先采用1 MHz的低通滤波获取有效信号；基于采样点的频率，使用低通FFT（fast Fourier transform）滤波器，截止频率设置为1 MHz[22]。图2显示了滤波后9 g RDX、12 g RDX和19 g RDX的水下爆炸压力时程曲线，同一质量炸药的实验重复2次。从图2可直接读取水下爆炸冲击波的峰值压力[23]。
图 2 不同实验工况下水下爆炸压力时程曲线Figure 2. Underwater explosion pressure-time curves under different experimental conditions冲击波的冲量I为压力对时间的积分[24]，可表示为：
I=∫6.7θ0p(t)dt (7) 比冲击波能[18]可表示为：
Es=4πr2mρwCw∫6.7θ0p2(t)dt (8) 式中：
ρw 为水的密度，取1000 kg/m3；Cw 为水中的声速，取1450 m/s；m为RDX的质量；Es为比冲击波能。3. 水下爆炸时程曲线衰减段的MATLAB拟合
3.1 水下爆炸时程曲线衰减段的模拟结果
为了验证新模型的有效性，采用MATLAB软件基于式(5)对实验的压力衰减过程进行拟合，得到拟合曲线和a、c两个常数。采用新模型、Cole模型、Орленко模型对9 g RDX、12 g RDX和19 g RDX实验进行拟合，结果如图3～5所示。
图 3 9 g RDX实验的3种模拟结果Figure 3. Three model fitting results of 9 g RDX experimental data
图 4 12 g RDX实验的3种模拟结果Figure 4. Three model fitting results of 12 g RDX experimental data
图 5 19 g RDX实验的3种模拟结果Figure 5. Three model fitting results of 19 g RDX experimental data拟合曲线的精度由可决系数R2决定：
R2=1−n∑i=1(yi−ˆyi)2n∑i=1(yi−ˉy)2 (9) 式中：
yi 为实验值，ˆyi 表示拟合结果，ˉy 为实验值的平均值，i为数据点的序号，n为数据点总数。R2的取值范围为0～1，数值越大，拟合精度越高。3种拟合方法的R2列于表1。从表1可以看出，新模型的R2均大于0.988，大多数结果大于0.99，说明新模型的模拟值与实验值接近，计算方法有效。相较于Cole和Орленко理论，新模型的R2更接近1，拟合精度更高，说明新模型能更准确地描述水下爆炸冲击波的衰减规律。
表 1 3种拟合方法的精度Table 1. Accuracy of three fitting methods炸药 组号 R2 a/kPa c/(kPa2·s) 新模型 Cole理论 Орленко理论 9 g RDX 1 0.9882 0.9282 0.9854 −0.097 0.053 2 0.9923 0.9363 0.9881 −0.089 0.043 12 g RDX 1 0.9936 0.9327 0.9861 −0.210 0.135 2 0.9963 0.9306 0.9861 −0.200 0.133 19 g RDX 1 0.9897 0.9274 0.9895 −0.180 0.131 2 0.9908 0.9286 0.9900 −0.160 0.094 3.2 冲击波参数评价
水下爆炸损伤效应的主要参数是I和Es。分别计算新模型、Cole-Орленко理论和实验的I与Es，验证新模型对冲击波参数模拟的准确性。由于Cole和Орленко理论分别适用于冲击波衰减的前半段和后半段，因此为计算完整的冲击波冲量与比冲击波能变化过程，将两种理论结合在一起计算。Cole-Орленко理论在计算冲量和比冲击波能时，将常用的Cole和Орленко理论代入式(7)～(8)，直接积分可得：
I=pmθ(1−1e+ln6.7e) (10) Es=4πθr2p2mmρwCw(12−12e2−16.7e2) (11) 新模型中，基于实验数据，通过式(7)～(8)积分得到冲击波参数。表2和表3分别列出了冲量和比冲击波能的实验与理论结果，其中δ1为新模型与实验结果的误差，δ2为Cole-Орленко模型与实验结果的误差，负数表示理论值小于实验值。
表 2 冲量的理论值与实验值Table 2. Theoretical and experimental values of impulse炸药 组号 I/(Pa·s) δ1(I)/% δ2(I)/% 实验 新模型 Cole-Орленко模型 9 g RDX 1 325.862 320.943 325.575 −1.510 −0.088 2 282.572 279.108 325.575 −1.226 15.280 12 g RDX 1 376.875 375.155 326.037 −0.456 −13.490 2 349.426 348.407 326.037 −0.292 −6.694 19 g RDX 1 439.536 439.484 431.403 −0.012 −1.850 2 400.488 400.228 431.403 −0.065 7.719 表 3 比冲击波能的理论值与实验值Table 3. Theoretical and experimental values of specific shock wave energy炸药 组号 Es/(MJ·kg−1) δ1(Es)/% δ2(Es)/% 实验 新模型 Cole-Орленко模型 9 g RDX 1 1.057 1.051 0.929 −0.578 −12.110 2 0.902 0.899 0.929 −0.370 2.993 12 g RDX 1 1.171 1.168 1.106 −0.256 −5.551 2 1.087 1.082 1.106 −0.460 1.748 19 g RDX 1 0.840 0.836 0.772 −0.429 −8.095 2 0.757 0.754 0.772 −0.396 1.982 由表2～3可以看出，新模型计算的冲量和比冲击波能具有更高的精度。从冲量来看，新模型结果与实验结果更接近，误差不超过5%，远低于Cole-Орленко模型结果；从比冲击波能来看，新模型结果与Cole-Орленко模型结果相差不大，新模型与实验结果之间的误差略小，并且误差不超过1%，而Cole-Орленко模型结果与实验结果的误差波动较大。
4. 结 论
通过对损伤效应的p-I公式进行推导，得到水下爆炸的压力时程公式。结合Cole和Орленко模型的优势，提出了一个新模型来描述水下爆炸的冲击波衰减过程，并通过水下爆炸实验得到新模型的经验参数。与Cole和Орленко模型相比，新模型与实验数据的近似程度更高，可决系数超过0.988。
-
图 1 黏弹性固体中爆炸地震波能量的组成
Figure 1. The composition of the energy of explosion seismic wave in visco-elastic solids
图 3 黄土中地下爆炸辐射地震波能量随观测区域的变化
Figure 3. The variation of the energy of the radiated seismic wave from underground explosion in loess with the observation region
图 4 黄土不同半径球面处最终流入能量
W(r,∞) 的变化Figure 4. The changes in the inflow energy
W(r,∞) at different radii of the sphere in loess
图 5 黄土中不同时刻粒子速度的空间分布
Figure 5. The spatial distribution of the particle velocity at different times in loess
图 6 黄土中不同时刻辐射动能的空间分布
Figure 6. The spatial distribution of the radiated kinetic energy at different times in loess
图 7 黄土中地震波辐射动能随波前位置的变化
Figure 7. The variation of the radiated kinetic energy of seismic wave with the position of wave front in loess
-
[1] 靳平, 徐果明, 楼为涛. 受低频动态正压力加载的椭球腔的地震矩张量表示及其在无限介质中辐射的地震波 [J]. 地震学报, 1997, 19(5): 447–456.JIN P, XU G M, LOU W T. Seismic moment tensor representation of ellipsoidal cavity loaded with low frequency dynamic positive pressure and seismic wave radiated in infinite medium [J]. Aata Seismologica Sinica, 1997, 19(5): 447–456. [2] CHOY G L, BOATWRIGHT J. Global patterns of radiated seismic energy and apparent stress [J]. Journal of Geophysical Research, 1995, 100(89): 18205–18228. DOI: 10.1029/95JB01969. [3] BOATWRIGHT J, CHOY G L. Teleseismic estimates of the energy radiated by shallow earthquakes [J]. Journal of Geophysical Research, 1986, 91(B2): 2095–2112. DOI: 10.1029/JB091iB02p02095. [4] MUELLER R A, MURPHY J R. Seismic characteristics of underground nuclear detonations: Part Ⅰ: seismic spectrum scaling [J]. Bulletin of the Seismological Society of America, 1971(61): 1675–1692. [5] MUELLER R A, MURPHY J R. Seismic characteristics of underground nuclear detonations: Part Ⅱ. Elastic energy and magnitude determinations [J]. Bulletin of the Seismological Society of America, 1971(61): 1693–1704. [6] 周钟, 王肖钧, 肖卫国, 等. 花岗岩介质中地下爆炸震源函数研究 [J]. 爆炸与冲击, 2007, 27(1): 18–25. DOI: 10.11883/1001-1455(2007)01-0018-08.ZHOU Z, WANG X J, XIAO W G, et al. Study on the main characteristics of underground explosion seismic source function in granite [J]. Explosion and Shock Waves, 2007, 27(1): 18–25. DOI: 10.11883/1001-1455(2007)01-0018-08. [7] MURPHY J R. Free-field seismic observations from underground nuclear explosions [M]. Explosion Source Phenomenology: American Geophysical Union Monograph, 1991. DOI: 10.1029/GM065p0025. [8] 袁乃荣, 刘瑞丰, 李赞, 等. 能量震级及其测定 [J]. 地震地磁观测与研究, 2018, 39(5): 1–7. DOI: 10.3969/j.issn.1003-3246.2018.05.001.YUAN N R, LIU R F, LI Z, et al. Energy magnitude and its determination [J]. Seismological and Geomagnetic Observation and Research, 2018, 39(5): 1–7. DOI: 10.3969/j.issn.1003-3246.2018.05.001. [9] 李赞, 刘瑞丰, 孔韩东, 等. 中强地震能量震级测定 [J]. 地震学报, 2019, 41(3): 289–301. DOI: 10.11939/jass.20180139.LI Z, LIU R F, KONG H D, et al. Energy magnitude determination of moderate-strong earthquakes [J]. Acta Seismologica Sinica, 2019, 41(3): 289–301. DOI: 10.11939/jass.20180139. [10] SANCHIDRIÁR J A, SEGARRA P, LÓPEM L M. Energy components in rock blasting [J]. International Journal of Rock Mechanics & Mining Sciences, 2007, 44: 130–147. [11] 田振农, 张乐文, 李世海. 岩体中爆腔内压力脉动特征和爆炸能量分布的数值模拟 [J]. 岩土工程学报, 2010, 32(8): 1247–1252. DOI: http://dspace.imech.ac.cn/handle/311007/43766.TIAN Z N, ZHANG L W, LI S H. Numerical simulation of pulsation features of pressure in explosion cavity and distribution of explosive energy in rock blasting [J]. Chinese Journal of Geotechnical Engineering, 2010, 32(8): 1247–1252. DOI: http://dspace.imech.ac.cn/handle/311007/43766. [12] 郭家豪, 范锦彪. 基于振动的钻地弹爆炸能量计算方法研究 [J]. 振动与冲击, 2020, 39(10): 180–184.GUO J H, FAN J B. Numerical calculation method for earth penetrating weapon explosion energy based on vibration [J]. Journal of Vibration and Shock, 2020, 39(10): 180–184. [13] 吴亮, 卢文波, 宗琦. 岩石中柱状装药爆炸能量分布 [J]. 岩土力学, 2006, 27(5): 734–738.WU L, LU W B, ZONG Q. Distribution of explosive energy consumed by column charge in rock [J]. Rock and Soil Mechanics, 2006, 27(5): 734–738. [14] 肖卫国, 王肖钧, 朱号锋, 等. 不同介质地下爆炸的地震耦合效应 [J]. 爆炸与冲击, 2012, 32(3): 267–272. DOI: 10.11883/1001-1455(2012)03-0267-06.XIAO W G, WANG X J, ZHU H F, et al. Experimental study on seismic coupling effects of underground explosions in different materials [J]. Explosion and Shock waves, 2012, 32(3): 267–272. DOI: 10.11883/1001-1455(2012)03-0267-06. [15] LU Q, WANG Z J. Studies of the propagation of viscoelastic spherical divergent stress waves based on the generalized Maxwell model [J]. Journal of Sound and Vibration, 2016, 371: 183–195. DOI: 10.1016/j.jsv.2016.02.034. [16] 卢强, 王占江. 标准线性固体材料中球面应力波传播特征研究 [J]. 物理学报, 2015, 64(10): 108301. DOI: 10.7498/aps.64.108301.LU Q, WANG Z J. Characteristics of spherical stress wave propagation in the standard linear solid material [J]. Acta Physica Sinica, 2015, 64(10): 108301. DOI: 10.7498/aps.64.108301. [17] LU Q, WANG Z J, Ding Y. Inversion for the complex elastic modulus of material from spherical wave propagation data in free field [J]. Journal of Sound and Vibration, 2019, 459: 1–18. DOI: 10.1016/j.jsv.2019.114851. [18] 卢强, 王占江, 李进, 等. 球面波加载下黄土线黏弹性本构关系 [J]. 岩土力学, 2012, 33(11): 3292–3298.LU Q, WANG Z J, LI J, et al. Linear viscoelastic constitutive relation of loess under spherical stress wave [J]. Rock and Soil Mechanics, 2012, 33(11): 3292–3298. [19] JAMES K G, JOHN Q E, CARGILE J D. Cavity expansion experiments with spherical explosive charges in concrete [R]. US Army Corps of Engineers, 2009. [20] WANG L L, LAI H W, WANG Z J, et al. Studies on nonlinear visco-elastic spherical waves by characteristics analyses and its application [J]. International Journal of Impact Engineering, 2013, 55: 1–10. DOI: 10.1016/j.ijimpeng.2012.12.001. [21] MILLER S A, FLORENCE A L. Laboratory particle velocity experiments on rock from a USSR underground nuclear tests site: AD-A223108 [R]. 1990. [22] 王占江, 李孝兰, 张若棋, 等. 固体介质中球形发散波的实验装置 [J]. 爆炸与冲击, 2000, 20(2): 103–109.WANG Z J, LI X L, ZHANG R Q, et al. An experimental apparatus for spherical wave propagation in solid [J]. Explosion and Shock Waves, 2000, 20(2): 103–109. [23] 卢强, 王占江, 朱玉荣, 等. 基于波传播系数构建填实爆炸下花岗岩中运动及变形场 [J]. 爆炸与冲击, 2019, 39(8): 1–10. DOI: 10.11883/bzycj-2019-0140.LU Q, WANG Z J, ZHU Y R, et al. Construction of motion and deformation field in granite under tamped explosion using wave propagation coefficient [J]. Explosion and Shock Waves, 2019, 39(8): 1–10. DOI: 10.11883/bzycj-2019-0140. -
下载:
下载:
计量
- 文章访问数: 694
- HTML全文浏览量: 262
- PDF下载量: 77
- 被引次数: 0