爆炸波在非饱和钙质砂中的传播规律

赵章泳; 王明洋; 邱艳宇; 紫民; 邢化岛

doi:10.11883/bzycj-2019-0389

爆炸波在非饱和钙质砂中的传播规律

doi: 10.11883/bzycj-2019-0389

赵章泳^{1, 2,},
王明洋^{2, 3},
邱艳宇^{2, 3, ,},
紫民⁴,
邢化岛³

1.
96911部队，北京 100010
2.
陆军工程大学爆炸冲击防灾减灾国家重点实验室，江苏南京 210007
3.
南京理工大学机械工程学院，江苏南京 210094
4.
海军91058部队，海南三亚 572000

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

详细信息

作者简介:
赵章泳（1992－　），男，博士，zhaozhangyong1@126.com

通讯作者:
邱艳宇（1978－　），女，博士，副教授，qiuyanyu78@sina.com

中图分类号: O382
计量
- 文章访问数: 4064
- HTML全文浏览量: 1394
- PDF下载量: 115
- 被引次数: 0
出版历程
- 收稿日期: 2019-10-16
- 修回日期: 2020-06-21
- 网络出版日期: 2020-07-25
- 刊出日期: 2020-08-01

The propagation laws of blast wave in unsaturated calcareous sand

ZHAO Zhangyong^{1, 2
,},
WANG Mingyang^{2, 3},
QIU Yanyu^{2, 3
, ,},
ZI Min⁴,
XING Huadao³

1.
96911 of PLA, Beijing 100010, China
2.
State Key Laboratory of Disaster Prevention and Mitigation of Explosion and Impact, The Army Engineering University of PLA, Nanjing 210007, Jiangsu, China
3.
School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu, China
4.
Navy 91058 Force, Sanya 572000, Hainan, China

摘要

摘要: 以球形TNT药包作为爆源，在密实的非饱和钙质砂中进行了一系列大尺寸爆炸模型试验。主要研究在不同药包质量、埋深及砂土试样含水率条件下，密实钙质砂中爆炸波的主要基本参数随传播距离增加而变化的规律。试验结果表明爆炸波主要以弹塑性波的形式在密实钙质砂中传播，在干燥和潮湿试样中塑性纵波波速随试样初始密度的增大或含水率的降低而增大，且范围分别为250～282 m/s和302～339 m/s。集团装药情况下，非饱和钙质砂中封闭爆炸的临界比例埋深约为2.25 m/kg^1/3。在试验范围内，密实钙质砂中爆炸波的法向应力峰值及法向比冲量的衰减均服从爆炸相似律。封闭爆炸时，干燥钙质砂中爆炸波的应力衰减指数在测点比例爆心距大于或小于0.75 m/kg^1/3处分别为2.94或1.37；潮湿钙质砂中爆炸波的应力衰减指数随含水率升高而增大，其范围为1.39～1.79。法向比冲量衰减指数随试样含水率升高而减小，其范围为0.97～1.18。
- 非饱和钙质砂 /
- 含水率 /
- 球形药包 /
- 波速 /
- 升压时间 /
- 法向应力峰值 /
- 法向比冲量
Abstract: A series of large-scale explosion model tests were carried out in dense unsaturated calcareous sand using spherical TNT explosives. The propagation laws of the blast wave in dense calcareous sand were studied under various conditions, such as explosive mass, buried depth and water content of sand sample, based on analyzing the changes of major parameters of blast wave. The results show that the blast wave propagates mainly in the form of elastic-plastic wave in dense calcareous sand. Moreover, the plastic longitudinal wave velocity increases with the increase of initial density in the dry sand sample. For the wet sand sample, the plastic longitudinal wave velocity increase with the decrease of the water content. More specifically, the corresponding longitudinal wave velocity ranges from 250 to 282 m/s in the dry sand and ranges from 302 to 339 m/s in the wet sand. In the case of concentrated charge, the critical scaled buried depth of closed explosion in unsaturated calcareous sand is about 2.25 m/kg^1/3. In the test range, the attenuations of peak normal stress and specific normal impulse of the blast wave in dense calcareous sand obey the explosion similarity law. The stress attenuation coefficient of blast wave in dry calcareous sand is 2.94 or 1.37 respectively at the measured points whose scaled distances are greater than or less than 0.75 m/kg^1/3. The stress attenuation coefficient of explosion wave in wet calcareous sand increases with the increase of water content and ranges from 1.39 to 1.79. The attenuation coefficient of the specific normal impulse decreases with the increase of the sample water content, and the range is 0.97 to 1.18.
- unsaturated calcareous sand /
- water content /
- spherical explosive /
- wave speed /
- rising time /
- peak normal stress /
- specific normal impulse

HTML全文

钙质砂是一种碳酸钙含量达50%以上的主要由海洋生物形成的粒状材料，由于在其沉积过程中大多未经长途搬运，从而保留了原生生物骨架中的细小孔隙，因此具有孔隙比高、压缩性高、颗粒易破碎等特点，并表现出与陆源砂不同的物理力学特性^[1–3]。近几十年来为保障海洋工程的建设，国内外学者已对钙质砂在低应变率下的力学特性进行了大量的研究^[4–6]，但关于钙质砂在高应变率条件下力学特性的研究还很少。研究钙质砂在爆炸荷载作用下的动力响应对岛礁及海洋工程中的防护工程问题具有重要意义，然而目前学者们仅对陆源砂抗爆性能进行过较为系统的研究，如美军的《TM5-855-1》手册提供了不同密实程度砂土中爆炸波传播的经验公式^[7]，许多学者在饱和、非饱和砂土中进行了大量爆炸试验，并建立了砂土在爆炸荷载作用下的计算模型^[8–17]。Karinski等^[18]、Yankelevsky等^[19]、Karinski等^[20]研究了土壤体积压缩特性对爆炸波在其中衰减规律的影响，并利用前人的试验结果进行验证。对于饱和钙质砂，徐学勇等^[21-22]使用电雷管作为点爆源对其进行加载，初步研究了爆炸波在其中传播所引起的土压力变化、孔隙水压力建立与消散以及砂土密实的规律。目前，关于非饱和钙质砂中爆炸波传播衰减规律的试验研究尚处于空白阶段。

本文采用球形TNT药包作为爆源，在密实钙质砂中进行了一系列大尺寸爆炸模型试验，主要研究了在不同砂土含水率、药包埋深、药包质量条件下密实钙质砂中爆炸波的传播衰减规律。所得试验结果既可用于建立适用于工程实践的经验公式，同时还可用于验证钙质砂本构模型及相应计算参数的准确性。

1. 爆炸相似律

由于非饱和砂土的成分较复杂，为完整描述其物理力学特性需要很多参数，包括密度、相对密度、孔隙比和含水率等^[23]，且非饱和砂土的结构性也对其力学特性有显著影响^[24–27]，因此目前试验研究仍为土中爆炸效应研究的最可靠手段。

爆炸波的基本参数主要包括爆炸波的波速、升压时间、法向应力峰值及对应的比冲量等。其中比冲量的定义为单位面积上爆炸波冲量的大小，并可根据下式计算得出，其量纲为[ML⁻¹T⁻¹]。

$I=\int_{0}^{\infty} \sigma(t) {\rm d} t$

(1)

对于爆炸波传播规律的主要研究对象是在爆炸波传播过程中其应力峰值、波速等基本参数随波传播距离的增加而发生的变化^[17]。目前尚未提出在较大的爆心距变化范围内，可用于求解这些基本参数的理论解，因此对这些参数变化规律的描述常采用半理论半经验的关系式。关系式的函数形式借助于爆炸相似律确定，而有关常数则是根据试验结果确定。爆炸相似律可阐述如下：两个尺寸不同、但化学成分相同的药包在任意形状介质的相同位置上爆炸时，两者的应力场和变形场在空间、时间和强度上均是相似的^[28]。

该理论的基本假设包括两点：一是介质对变形速率不敏感，二是介质中的非定常爆炸应力场和变形场只受爆炸能量的影响（不考虑重力和其他力的影响）。虽然该理论的基本假设具有局限性，但是以往研究结果表明，在相当程度上仍能利用该理论设计爆炸模型试验。

爆炸相似律的函数形式可以通过量纲分析的手段得出^[29]。根据以往岩土介质中爆炸模型试验的分析^[30–34]，可以汇总得到在集团装药情况下影响砂土中爆炸波传播规律的主要参数如表1所示。

表 1 影响砂土中爆炸波传播规律的主要参数

Table 1. Major parameters influencing the propagation of blast wave in sand

	参数名称及符号	量纲
炸药	药包质量W	[M]
	装药密度ρ_w	[ML^-3]
	单位质量炸药的化学能E_w	[L²T^-2]
	爆炸产物的膨胀指数γ	—
砂土	密度ρ₀	[ML^-3]
	塑性纵波波速c_p	[LT^-1]
	内摩擦角φ	—
试样尺寸	药包中心埋深h	[L]
试样尺寸	测点到爆心的距离R	[L]

下载: 导出CSV

| 显示表格

爆炸波的法向应力峰值σ_max及相应的折合比冲量I应当是这些控制参数的函数，即有：

$\begin{split} \sigma_{\max }=f\left(W, \rho_{\mathrm{w}}, E_{\mathrm{w}}, \gamma, \rho_{0}, c_{\mathrm{p}}, \varphi, h, R\right), \quad I=g\left(W, \rho_{\mathrm{w}}, E_{\mathrm{w}}, \gamma, \rho_{0}, c_{\mathrm{p}}, \varphi, h, R\right) \end{split}$

(2)

取药包质量W、装药密度ρ_w、单位质量炸药的化学能E_w为基本量，式(2)可化为如下的无量纲关系：

$\begin{split} &\frac{\sigma}{\rho_{\mathrm{w}} E_{\mathrm{w}}}=f\left(\gamma, \frac{\rho_{0}}{\rho_{\mathrm{w}}}, \frac{c_{\mathrm{p}}}{\sqrt{E_{\mathrm{w}}}}, \varphi, \frac{h}{\left(\dfrac{W}{\rho_{\mathrm{w}}}\right)^{\frac{1}{3}}}, \frac{R}{\left(\dfrac{W}{\rho_{\mathrm{w}}}\right)^{\frac{1}{3}}}\right) \\ &\frac{I}{W^{\frac{1}{3}} \rho_{\mathrm{w}}^{\frac{2}{3}} E_{\mathrm{w}}^{\frac{1}{2}}}=g\left(\gamma, \frac{\rho_{0}}{\rho_{\mathrm{w}}}, \frac{c_{\mathrm{p}}}{\sqrt{E_{\mathrm{w}}}}, \varphi, \frac{h}{\left(\dfrac{W}{\rho_{\mathrm{w}}}\right)^{\frac{1}{3}}}, \frac{R}{\left(\dfrac{W}{\rho_{\mathrm{w}}}\right)^{\frac{1}{3}}}\right) \end{split}$

(3)

对于采用相同种类的炸药及装药方式，在相同的钙质砂中进行的炸药质量不同的爆炸试验，则有6个控制参数保持不变，即：

${\left( {{\rho _{\rm{w}}},{E_{\rm{w}}},\gamma ,{\rho _0},{c_{\rm{p}}},\varphi } \right) = {\rm{ const}}}$

(4)

式(3)可化简为：

$\begin{array}{l} {\sigma _{\max }} = f\left( {{R^*},{h^*}} \right),\quad\quad\quad {I^*} = g\left( {{R^*},{h^*}} \right) \end{array}$

(5)

式中：f和g通常取为幂函数的形式，相应参数由模型试验得出。R^*为比例爆心距离，h^*为药包比例埋深，I^*为折合比冲量，相应的计算公式分别为：

$\begin{aligned} R^{*} =\frac{R}{\sqrt[3]{W}}, \quad\quad\quad h^{*} =\frac{h}{\sqrt[3]{W}}, \quad\quad\quad I^{*} =\frac{I}{\sqrt[3]{W}} \end{aligned}$

(6)

2. 试验设计

2.1 砂土试样

试验所用的钙质砂取自中国南海某岛，首先将钙质砂原样平摊于广场上风干，然后筛分出砂土试样中粒径大于2 mm的颗粒。按照《土工试验方法标准》（GB/T 50123—2019）^[35]得到筛分后的钙质砂试样的级配曲线如图1所示，试样的最小干密度与最大干密度分别为1.12和1.47 g/cm³，颗粒密度为2.81 g/cm³。计算可得试样的平均粒径d₅₀=0.348 mm，不均匀因数C_u＞6，因此属于不均匀中砂。

图 1 钙质砂试样级配曲线

Figure 1. Particle size distribution curve of calcareous sand

下载: 全尺寸图片幻灯片

2.2 试验装置及球形药包

试验使用如图2所示的拼装式筒体作为爆炸试验容器。容器总高度为2 000 mm，由10个拼装单元沿垂直方向连接而成。拼装单元由2片法兰盘和1个圆环焊接而成：其中圆环外径为2 000 mm、壁厚为15 mm，共有高度分别为100、200和500 mm的3种尺寸；法兰盘宽度为75 mm，厚度为16 mm。为提高法兰与圆环的连接强度，还在拼装单元侧面设置了12道加劲肋板。每一层拼装单元侧壁上均加工有螺纹孔，以方便传感器导线的引出。

图 2 爆炸模型试验所用容器

Figure 2. The container used in explosion experiments

下载: 全尺寸图片幻灯片

关于砂土中爆炸效应的试验研究中，所使用的爆源主要包括雷管^[22]、乳化炸药^[36]和块装TNT^[31-32]。当爆心距较大时，这些爆源可以近似视为点源。然而在爆心距较小的位置，炸药类型、装药密度以及炸药形状均会对试验结果带来较大的影响，Krauthammer^[37]特别强调了需要慎重使用由非球形爆源试验所获得的爆炸近区数据。为研究在爆心距较小范围内钙质砂中的爆炸波的传播规律，本文中将采用球形药包作为爆源。

图3(a)为试验所用球形药包的示意图，药包由上下2个部分粘贴而成。每部分均由散装三硝基甲苯（TNT）压装而成，压药密度为1.5 g/cm³。其中药包的下半部为完整半球，上半部则预留有直径为8 mm的雷管安装孔。将2个半球粘接为一体，通过测量可知虫胶层厚度小于0.1 mm。试验所用球形药包具有3种规格，其名义质量分别为64、216、512，相应的名义直径分别为43.6、65.2和87.0 mm。由于在药球的脱模过程中，雷管安装孔附近会有少量炸药脱落，因此每次试验前都需要称量药包的实际质量。图3 (b)为尚未粘接的216和512 g药包，图3 (c)为粘接好的64 g球形药包。由于药球的压装密度低于铸装TNT密度，因此其感度较高可以直接使用电雷管起爆。

图 3 球形药包

Figure 3. Spherical explosives

下载: 全尺寸图片幻灯片

2.3 测试系统

本文采用高频压阻式土压力传感器测量不同测点处爆炸波法向应力的时程曲线。图4所示为试验采用的土压力传感器：图左为昆山双桥公司生产的CYG1712F型土压力传感器，其尺寸为 $\varnothing$ 30 mm×16 mm；图右为河南动感电子科技公司生产的DNS-TY型土压力传感器，其尺寸为 $\varnothing$ 40 mm×10 mm。

图 4 土压力传感器

Figure 4. Soil pressure sensors

下载: 全尺寸图片幻灯片

传感器输出的模拟信号由东华测试公司的DH8302高性能动态信号测试系统调理并采集。测试采用连续采集模式。

图5为药包及传感器位置的设计示意图。其中药包位于容器的几何中心，在钙质砂试样制备过程中预留装药孔，待试样准备完毕后将捆绑好电雷管的药包放入装药孔并回填。传感器的安装平面指传感器底面所在平面，在布置传感器时需要使用此平面对传感器进行定位；测试平面指传感器顶面（也即敏感元件）所在平面。考虑到砂土的不均匀性，在每个传感器安装平面内均布置1个中心传感器与1～2个边缘传感器，且中心传感器位于药包正下方。各传感器中心之间的水平距离设置为6 cm，以减小传感器对砂土自由场的扰动。由于第一层测点距离爆心较近，因此只能设置一个中心传感器。试验中容器将被拆分为上、下两部分，以便进行夯实、取样和安装传感器等操作。

图 5 药包及传感器位置示意图

Figure 5. The schematic of the locations of explosive and sensors

下载: 全尺寸图片幻灯片

2.4 试验流程

共进行8次爆炸试验，构成了分别以药包埋深、药包质量和试样含水率为单变量的3组对照试验，各试验的设计参数如表2所示。在试验编号中，DCS表示干燥钙质砂，WCS表示潮湿钙质砂；64、216和512表示药包设计质量。其中干燥钙质砂由原始试样经风干后所得，通过室内试验测定其平均含水率约为1%；潮湿钙质砂是通过将已知含水率的钙质砂与相应质量的水使用搅拌机混合后制备获得。

表 2 爆炸模型试验参数

Table 2. Designed parameters for explosion model test

试验编号	药包设计质量W/g	设计埋深h/m	设计含水率/%	设计爆心距R/m
DCS64-1	64	0.9	1	0.1，0.2，0.3，0.4，0.6，0.8，1.0
DCS64-2	64	0.3	1	0.1，0.2，0.3，0.4，0.6，0.8，1.0
DCS64-3	64	0	1	0.1，0.2，0.3，0.4，0.6，0.8，1.0
DCS216	216	0	1	0.15，0.2，0.3，0.45，0.6，0.8，1.0
DCS512	512	0	1	0.2，0.3，0.4，0.6，0.8，1.0
WCS64-1	64	0.9	10	0.15，0.2，0.3，0.4，0.6，0.8，1.0
WCS64-2	64	0.9	20	0.1，0.2，0.3，0.4，0.6，0.8，1.0
WCS64-3	64	0.9	25	0.1，0.2，0.3，0.4，0.6，0.8，1.0

下载: 导出CSV

| 显示表格

为提高法向应力测试结果的信噪比，需预先对测点的应力峰值进行估计，以合理设置每个测试层内传感器的量程。然而目前尚未有文献给出钙质砂中爆炸波的试验数据，因此在试验设计中使用了梁霍夫^[17]关于非饱和石英砂中爆炸波的试验结果对各测点处应力幅值进行预估。由于当容器装满钙质砂试样后难以进行夯实操作，因此本节试验中最大设计埋深为0.9 m，即略小于容器上半层的高度。

考虑到模型试验的尺寸较大，按试样干密度1.4 g/cm³计算则约需8.3 t干燥钙质砂试样。为使用有限量的钙质砂进行不同试验工况条件下的对照试验，同时缩短试验周期，将按照试样含水率由低至高的顺序进行表2中所示的试验。

考虑到钙质砂试样易破碎，在每次试验结束后，均需对各测试层试样进行回收并进行颗粒分析，舍弃砂土颗粒明显破碎的试样，以保证每次试验中试样颗粒特性一致。颗粒分析结果表明，在比例爆心距R^*＞0.5 m/kg^1/3的范围内，不同含水率试样中的钙质砂颗粒几乎不发生破碎。对于药包质量为64 g的试验，此范围内砂样的体积仅为试样总体积的5‰。因此在每次试验结束后，仅需将上述范围内的试样抛弃，即可保证试样的颗粒分配特性不随系列爆炸试验的进行而发生明显变化。

准备试件时，首先使用悬挂于叉车臂上的吊秤将盛放钙质砂的吨袋吊至容器上方，称量每一袋试样的总质量。适当打开出料口，使试样缓慢漏出，砂土将在容器内形成圆锥状堆积体。手动调节叉车臂的高度使其随堆积体高度的增大而增大，控制吨袋底部与堆积体顶部的高度差约为0.5 m，此时试样约以3.2 m/s的速度落入容器内。待试样完全从吨袋中落入容器后，使用铁锹将倒入容器之中的砂土摊平，并使用电动冲击夯对试样进行夯实。夯机质量为80 kg，电机功率为3 kW，夯实时起跳高度约为5 cm，相应每次夯击输入能量为40 J，每层试样总夯击时间约为5 min。由于在夯实操作前已将部分土压力传感器埋入试样中，因此可以使用这些传感器测量夯击产生的应力波的信号。测试结果表明由夯击产生的应力波的应力幅值约为0.015 MPa，夯击周期约为100 ms。

按照图5所示进行试样的分层装填与夯实、传感器的布设及药包的埋置。在埋设好每层传感器后，使用环刀法在该层距圆心40 cm位置处沿圆周方向取3个样本，对其称重即可计算该层试样的初始密度。与干燥钙质砂相比，潮湿钙质砂具有一定程度的假黏聚力^[23]，因此其能够维持较稳定的形状。这种特性使得对于潮湿钙质砂试验，在整个容器内的试样均被装填夯实完毕后，可以使用环刀取土钻获取不同深度处的试样，进而计算不同深度处试样的密度。通过对比由试样总质量除以总体积法、环刀法及环刀取土钻法所获得的试样密度，可以确定上层试样的夯击对下层试样密度的影响程度。

2.5 实测试验参数

表3总结了各次爆炸模型试验的实际参数。对比实测结果发现，干燥钙质砂的平均含水率随试验的进行略有上升。这是因为试验周期较长，孔隙率很高的钙质砂试样从空气中吸收了部分水蒸气。

表 3 实测爆炸试验参数表

Table 3. Summary of measured parameters in explosion test

试验编号	药包质量/g	试样含水率/%	药包埋深/m	试样密度/(g·cm⁻³)
试验编号	药包质量/g	试样含水率/%	药包埋深/m	总质量除以总体积法	环刀取样法	环刀取土钻取样法
DCS64-1	63.6	0.7±0.5	0.90	1.37	1.35±0.02	−
DCS64-2	63.1	1.5±0.4	0.31	1.39	1.36±0.04	−
DCS64-3	63.3	1.9±0.5	0	1.38	1.36±0.02	−
DCS216	214.2	2.1±0.5	0	1.39	1.37±0.03	−
DCS512	510.2	2.0±0.3	0	1.38	1.37±0.02	−
WCS64-1	62.8	10.5±0.7	0.89	1.54	1.51±0.04	1.54±0.04
WCS64-2	63.3	18.9±1.5	0.92	1.59	1.55±0.07	1.58±0.06
WCS64-3	63.1	23.7±2.2	0.91	1.80	1.76±0.06	1.81±0.04

下载: 导出CSV

| 显示表格

使用环刀取土钻的方法测得的密度与使用总质量除以总体积法测得的密度差别不大，但均比环刀取样法测得的密度大0.02～ 0.05 g/cm³，这是因为在使用环刀取样后进行的夯击会略微增加已取样砂土层中试样的密度。

考虑到起爆所用电雷管折合TNT当量为1 g，试样中的实际炸药质量与设计药包质量相差不超过1%，另外药包的埋深与设计值相差不超过3%，因此在下文的计算中对于药包质量和埋深仍采用设计值。

3. 试验结果与分析

3.1 试验现象及应力时程曲线

图6 (a)～(f)分别为药包埋深等于0.9、0.3和0 m时，在爆炸试验前后试样的表面状态。由图6(a)～ (b)可知，当药包埋深为0.9 m的试验，试验结束后试件表面除有少许浮土外几乎没有任何明显变化；由图6(c)～ (d)可知，当药包埋深为0.3 m时，在爆心上方的试样被抛掷出去并形成爆破漏斗；由图6(e)～ (f)可知，当药包埋深为0 m时，在试样表面将出现深度较小、半径较大的爆炸坑。这3种药包埋深条件下对应的爆炸类型分别为封闭爆炸、浅埋爆炸和触地爆炸。

图 6 试验前后试件表面状态变化

Figure 6. Changes of sample surface state after the explosion tests

下载: 全尺寸图片幻灯片

如图6(e)所示在DCS64-2试验中形成的漏斗坑，其宽度为80 cm，深度为27 cm。试样表面布满了黑灰色的爆轰产物，说明爆腔上部的砂土被完全抛散出去，这种漏斗坑一般被称为完全漏斗坑^[38]。

图7所示为试验中传感器记录到的全部原始数据。考虑到砂土的不均匀性，在每个传感器安装平面内均布置1个中心传感器与1～2个边缘传感器。图7中所示爆心距下对应的多条曲线即表示布置的多个传感器。除部分爆心距较小位置处传感器导线受爆炸波作用损坏导致未采集到数据外，绝大多数传感器均完整记录了测点处的爆炸波法向应力时程曲线。试验结果信噪比较高，相同爆心距处传感器记录波形一致性较好。由于爆轰产物运动会产生电磁脉冲^[34]，这些电磁脉冲耦合入测试电路致使部分试验结果中出现了高频振荡的电磁噪声。除DCS512试验的第一层传感器外，这种电磁噪声基本不会对应力测试结果的分析造成任何影响。此外，由于试验操作不当，导致在WCS64-3试验中药包埋设位置比设计位置高出3 cm，因此该试验中所有测点的爆心距均增大了3 cm。

图 7 应力时程曲线

Figure 7. Stress-time curves

下载: 全尺寸图片幻灯片

3.2 爆炸波波速与升压时间

爆炸波的波速主要包括波阵面的传播速度和应力峰值的传播速度。将4种不同含水率的封闭爆炸试验中波阵面和应力峰值的到时汇总于图7中，由图7可知在试验的测试范围内，爆心距与波阵面到时、应力峰值到时均满足线性关系，因此密实钙质砂中的波速基本保持为定值。且除DCS64-1试验中第一层测点外，在其余测点爆炸波的到时和应力峰值到时有很大区别。由于波阵面到时及应力峰值到时的特征与弹塑性波传播规律一致，根据弹塑性球面波的传播理论可知，其波阵面的波速为试样的弹性纵波波速C₀，因此可以将爆心距-波阵面到时线性拟合关系的斜率作为弹性纵波波速；应力峰值的波速为试样的塑性纵波波速C_p，因此可以将爆心距-应力峰值到时线性拟合关系的斜率作为塑性纵波波速。

表4汇总了所有试验中的波速拟合值，根据表中数据可知：在密实干燥钙质砂中弹性纵波波速为291～359 m/s，塑性纵波波速为250～282 m/s；在密实潮湿钙质砂中弹性纵波波速为372～427 m/s，塑性纵波波速为302～339 m/s。干燥钙质砂中弹性纵波波速离散性较大，且与试样密度的变化不具有一致性，因此钙质砂试样在低应力下的力学特性受不限于试样密度的多种因素影响；而干燥钙质砂塑性纵波波速随试样密度的增大而增大，因此试样密度是干燥钙质砂高应力下变形特性的主要影响因素。潮湿钙质砂中的弹性和塑性纵波波速均随含水率的升高而增大，因此含水率是潮湿钙质砂在低应力和高应力条件下力学特性的主要影响因素。

表 4 密实钙质砂中的爆炸波波速

Table 4. The blast wave velocity in dense calcareous sand

试验编号	密度/(g·cm⁻³)	弹性纵波波速c₀/(m·s⁻¹)	R²	塑性纵波波速c_p/(m·s⁻¹)	R²
DCS64-1	1.37	359	0.998	250	0.995
DCS64-2	1.39	344	0.998	282	0.995
DCS64-3	1.38	291	0.998	265	0.997
DCS216	1.39	332	0.999	278	0.999
DCS512	1.38	339	0.999	274	0.999
WCS64-1	1.50	429	0.999	377	0.998
WCS64-2	1.59	387	0.998	315	0.996
WCS64-3	1.80	374	0.999	302	0.999

下载: 导出CSV

| 显示表格

由于在试验条件下，爆炸波的弹性波速和塑性波速均基本保持不变，则在爆心距为R处其上升时间t_r可近似表示为：

$t_{\mathrm{r}}=R\left(\frac{1}{c_{0}}-\frac{1}{c_{\mathrm{p}}}\right)$

(7)

该式表明，密实钙质砂中爆炸波的升压时间与爆心距成线性关系。

3.3 爆炸波法向应力峰值的衰减规律

爆炸波法向应力峰值σ_n,_max的衰减公式通常是以如式(8)所示的幂函数的形式给出：

$\sigma_{{n},{\max }}=K_{\sigma}\left(R^{*}\right)^{-\mu_{\sigma}}$

(8)

式中： $K_{\sigma}$ 为应力衰减常数，代表比例距离等于1 m/kg^1/3时介质中的应力峰值； $\mu_{\sigma}$ 为应力衰减指数，代表着爆炸波应力峰值的衰减规律，其值越大则应力峰值衰减速度越快。

图8为在对数坐标系下，所有干燥钙质砂试验中测点的比例爆心距和法向应力峰值的对应关系。由图8可知，干燥钙质砂的应力峰值衰减规律以R^*=0.75 m/kg^1/3为界可划分为两段，表5汇总了这两段衰减规律的参数值。

表 5 干燥钙质砂法向应力峰值衰减规律参数值

Table 5. The parameters of attenuation law of peak normal stress in dry calcareous sand

试验编号	$K_{\sigma}$	$\mu_{\sigma}$	R²	$K_{\sigma}$	$\mu_{\sigma}$	R²
试验编号	第一段(R*≤0.75 m/kg^1/3)			第二段(R*＞0.75 m/kg^1/3)
DCS64-1	0.944	2.94	0.998	1.474	1.37	0.976
DCS64-2	0.944	2.94	0.998	1.474	1.37	0.976
DCS64-3	0.548	3.11	0.961	0.79	1.83	0.938
DCS216
DCS512

下载: 导出CSV

| 显示表格

图 8 封闭爆炸试验中爆心距与到时关系

Figure 8. The relationship between blast center distance and arrival time in closed explosion

下载: 全尺寸图片幻灯片

根据表5的汇总结果及图9可知：在R^*＞0.33 m/kg^1/3的范围内试验结果一致性较好且符合爆炸相似律，而在R^*=0.25 m/kg^1/3的测点处不同药包质量的试验结果偏差较大。其原因可能是当应力水平较高时，应力值及相应的测试手段对介质的不均匀性较敏感；也可能是在该比例距离处，钙质砂动力特性的率敏感性将对爆炸相似率的适用性产生影响。具体是由哪种原因造成的这种偏差，尚需进一步研究。

图 9 干燥钙质砂法向应力峰值衰减的试验结果及拟合曲线

Figure 9. Experiment results and fitting curves of the peak normal stress in dry calcareous sand

下载: 全尺寸图片幻灯片

图9(b)表明比例埋深由2.25 m/kg^1/3减小至0.75 m/kg^1/3时，钙质砂中爆炸波应力峰值的衰减规律不会发生变化。当比例埋深减小至0 m/kg^1/3时，在R^*=0.25，2.0 m/kg^1/3处的应力峰值将分别下降27%和60%，且衰减指数会略微增大。

图10为对数坐标系下，所有潮湿钙质砂试验中测点的比例爆心距和法向应力峰值的对应关系，相应的衰减规律参数值汇总于表6中。可以发现在比例爆心距在0.5～2.0 m/kg^1/3的范围内，潮湿钙质砂中的应力峰值衰减曲线并没有出现如同干燥钙质砂中的分段现象，法向应力峰值随含水率的升高先减小后增大，但衰减指数μ_σ却随含水率的升高而单调递增。

表 6 潮湿钙质砂法向应力峰值衰减规律参数值

Table 6. The parameters for attenuation law of peak normal stress in moist calcareous sand

试验编号	$K_{\sigma}$	$\mu_{\sigma}$	R²
WCS64-1	1.805	1.39	0.965
WCS64-2	1.530	1.55	0.986
WCS64-3	2.200	1.79	0.978

下载: 导出CSV

| 显示表格

图 10 潮湿钙质砂法向应力峰值试验结果及拟合曲线

Figure 10. Experiment results and fitting curves of the peak normal stress in moist calcareous sand

下载: 全尺寸图片幻灯片

对于R=100 cm测点处的传感器，其被放置于容器底板之上，而底板下方为混凝土地基。由于传感器、容器底板及混凝土地基的波阻抗均远大于钙质砂的波阻抗，因此可以将该层传感器视为刚壁，其所测应力时程曲线的应力峰值即为实测反射波应力峰值σ_R,_max。使用式(8)和表5～6中的参数值，可以通过外推的方法获得在100 cm测点处的自由场应力峰值σ_I,_max，σ_R,_max与σ_I,_max之比即为非饱和钙质砂中的刚壁反射因数。

表7汇总了所有试验中的外推自由场应力峰值、实测反射波应力峰值及刚壁反射因数，可以发现：对于含水率低于等于10%的非饱和钙质砂，刚壁反射因数在1.473～1.603之间；对于含水率高于20%的非饱和钙质砂，刚壁反射因数在1.225～1.249之间。

表 7 非饱和钙质砂刚壁反射因数

Table 7. The rigid reflection coefficient of unsaturated calcareous sand

试验编号	外推自由场应力峰值σ_I,_max/MPa	实测反射波应力峰值σ_R,_max/MPa	刚壁反射因数σ_R,_max/σ_I,_max
DCS64-1	0.421	0.675	1.603
DCS64-2	0.421	0.635	1.508
DCS64-3	0.149	0.219	1.473
DCS216	0.311	0.458	1.473
DCS512	0.529	0.796	1.504
WCS64-1	0.505	0.804	1.592
WCS64-2	0.368	0.451	1.225
WCS64-3	0.430	0.537	1.249

下载: 导出CSV

| 显示表格

3.4 爆炸波法向折合比冲量的衰减规律

当爆炸波传至底板后会产生反射波，反射波会对爆心距较大的测点（主要是R=60，80 cm）处的波形造成显著影响，图11中的实线为DCS64-1试验中R=80 cm测点处的原始波形。可以发现反射波信号并不会影响爆炸波的应力峰值的读取，但却给比冲量的计算带来巨大误差。为消除反射波的影响，依据梁霍夫^[17]的关于同一测点处爆炸波法向应力随时间而衰减的研究结果，使用指数函数对未受到底板反射信号影响的波形进行拟合，然后将其外推至试验测试结束。图11中的虚线即为拟合及外推的结果。

图 11 消除反射波对应力时程曲线影响的方法

Figure 11. The method for eliminating the influence of reflected wave on stress-time curve

下载: 全尺寸图片幻灯片

将不同传感器测得的法向应力时程曲线按照图11所示的方法进行修正后进行积分，即可得到相应位置的比冲量时程曲线，按照该方法计算DCS64-1试验中各测点的比冲量时程曲线如图12所示。由图可知在时间达到10 ms后所有位置处的比冲量基本保持为一稳定值，即爆炸波脉冲已完全通过该测点，相应的比冲量稳定值即为此位置处爆炸波的总比冲量值，该值对爆炸荷载作用下地下防护结构的动力响应计算及相应的工程设计具有重要意义。

图 12 DCS64-1试验中各测点的比冲量时程曲线

Figure 12. Specific impulse time histories at each measuring point in DCS64-1 test

下载: 全尺寸图片幻灯片

将不同测点处爆炸波的总比冲量代入式(6)中第3式，即可得到各测点的折合比冲量值。折合比冲量的衰减公式通常是以如式(9)所示的幂函数的形式给出：

$I^{*}=K_{ I}\left(R^{*}\right)^{-\mu_{\rm I}}$

(9)

式中：K_I为比冲量衰减常数，μ_I为比冲量衰减指数。

在同一测试层内不同传感器测得的比冲量相差由5%～50%不等，其离散性明显大于应力峰值试验结果的离散性。这是因为在卸载过程传感器的位移远大于加载过程的位移，由于波阵面曲率的影响，将导致传感器感压面方向发生偏转，从而给应力的测量带来较大误差。因此即使在同一测试层内各传感器测得的应力峰值几乎相同，但其卸载曲线可能存在显著差异。由于卸载时间远比加载时间长，因此卸载曲线的差异将对比冲量计算造成显著影响。

图13～14所示分别为在干燥钙质砂和潮湿钙质砂中折合比冲量的试验结果及相应的拟合曲线。可以发现在对数坐标系下，干燥钙质砂或潮湿钙质砂折合比冲量的衰减均为一段直线的形式，这一点与干燥钙质砂中应力峰值的衰减规律不同。对于干燥钙质砂，比例埋深由2.25 m/kg^1/3减小至0.75 m/kg^1/3时，虽然应力峰值的大小不变，但是折合比冲量却显著降低。

图 13 干燥钙质砂折合比冲量试验结果及拟合曲线

Figure 13. Experiment results and fitting curves of the specific impulse in dry calcareous sand

下载: 全尺寸图片幻灯片

图 14 潮湿钙质砂折合比冲量试验结果及拟合曲线

Figure 14. Experimental results and fitting curves of the specific impulse in moist calcareous sand

下载: 全尺寸图片幻灯片

表8汇总了非饱和钙质砂中折合比冲量衰减公式的参数值。可以发现当含水率不变时，比例埋深仅会影响折合比冲量的值而不影响其衰减指数，即不影响其衰减规律。对于潮湿钙质砂，折合比冲量的衰减指数随含水率的增加而减小，且10%含水率试样与干燥试样的衰减指数近似相同。

表 8 非饱和钙质砂中折合比冲量衰减公式参数值

Table 8. The parameters of attenuation law of specific impulse in unsaturated calcareous sand

试验编号	K_I	μ_I	R²
DCS64-1	7.57	1.15	0.944
DCS64-2	5.54	1.14	0.975
DCS64-3	1.57	1.16	0.933
DCS216
DCS512
WCS64-1	5.33	1.18	0.959
WCS64-2	4.92	1.02	0.901
WCS64-3	4.65	0.97	0.946

下载: 导出CSV

| 显示表格

将不同含水率的封闭爆炸试验中折合比冲量的衰减拟合曲线汇总于图15中，可以发现随着含水率的升高，在相同比例爆心距位置处折合比冲量的值与衰减速度均随含水率的升高而降低。

图 15 不同含水率钙质砂中的折合比冲量衰减

Figure 15. Attenuation of reduced specific impulse in calcareous sand with different moisture content

下载: 全尺寸图片幻灯片

耦合因数f的定义是非封闭爆炸与封闭爆炸在同一介质汇总所产生的地冲击大小比值，其取值范围为0～1，通常按照下式计算求得：

$f\left(h^{*}\right)=\frac{\left\{\sigma_{n, \max} , I^{*}, v, d, a\right\}_{ {\text{非封闭爆炸}} }}{\left\{\sigma_{n, \max} , I^{*}, v, d, a\right\}_{{\text{封闭爆炸}} }}$

(10)

式中：σ_n_,max为法向应力峰值，I^*为折合比冲量，v为峰值粒子速度，d为峰值位移，a为峰值加速度。一般情况下，上述不同的量所对应的f是不相同的。在关于爆炸近区地冲击的研究中，选取较多的是法向应力峰值和比冲量，相应的耦合因数分别称为应力型耦合因数和冲量型耦合因数。由于钙质砂中爆炸波峰值应力的的衰减随比例爆心距的增大而增强，因此不适宜用式(10)的形式表示。根据计算可得，在比例埋深为0和0.75 m/kg^1/3时，钙质砂的冲量型耦合因数f分别为0.21和0.73。

4. 结　论

本文进行了不同埋深及药包质量下干燥钙质砂中的爆炸模型试验，以及不同含水率的潮湿钙质砂中的封闭爆炸模型试验。首次获得了比例爆心距在0.25～2.0 m/kg^1/3范围内，爆炸波在不同含水率非饱和钙质砂中的自由场应力时程数据，分析了密实钙质砂中药包埋深及试样含水率对爆炸波传播衰减规律的影响，主要结论包括：

（1）集团装药情况下，非饱和钙质砂中封闭爆炸的临界比例埋深约为2.25 m/kg^1/3。

（2）在密实干燥钙质砂中弹性纵波波速为291～359 m/s，塑性纵波波速为250～282 m/s。在密实潮湿钙质砂中弹性纵波波速为372～427 m/s，塑性纵波波速为302～339 m/s。试样密度及含水率是非饱和钙质砂在高应力条件下波速的主要影响因素。密实钙质砂中爆炸波的升压时间随爆心距的增加呈线性增长。

（3）封闭爆炸情况下，在比例爆心距小于0.75 m/kg^1/3时干燥钙质砂中爆炸波的应力衰减因数为2.94，在比例爆心距大于0.75 m/kg^1/3时为1.37。比例埋深减小至0.75 m/kg^1/3时，应力峰值的大小及衰减规律均不变；比例埋深减小为0 m/kg^1/3时在不同的比例爆心距处应力峰值将下降27%～60%。封闭爆炸情况下，潮湿钙质砂中爆炸波的应力衰减因数在比例爆心距为0.5～2.0 m/kg^1/3的范围内为1.39～1.79，且随含水率的升高而增大。

（4）对于含水率低于等于10%的非饱和钙质砂，其应力峰值的刚壁反射因数在1.473～1.603之间；对于含水率大于等于20%的非饱和钙质砂，其应力峰值的刚壁反射因数在1.225～1.249之间。

（5）相同比例爆心距处爆炸波的折合比冲量值随着含水率的升高而减小，其范围为0.97～1.18。10%含水率的钙质砂与干燥钙质砂的比冲量衰减指数近似相同；当含水率高于10%时，比冲量衰减指数随含水率的升高而减小。比例埋深仅会影响折合比冲量的大小，而不会影响其衰减规律。在比例埋深为0和0.75 m/kg^1/3时干燥钙质砂中冲量型耦合因数分别为0.21和0.73。

图 1 钙质砂试样级配曲线

Figure 1. Particle size distribution curve of calcareous sand

下载: 全尺寸图片幻灯片

图 2 爆炸模型试验所用容器

Figure 2. The container used in explosion experiments

下载: 全尺寸图片幻灯片

图 3 球形药包

Figure 3. Spherical explosives

下载: 全尺寸图片幻灯片

图 4 土压力传感器

Figure 4. Soil pressure sensors

下载: 全尺寸图片幻灯片

图 5 药包及传感器位置示意图

Figure 5. The schematic of the locations of explosive and sensors

下载: 全尺寸图片幻灯片

图 6 试验前后试件表面状态变化

Figure 6. Changes of sample surface state after the explosion tests

下载: 全尺寸图片幻灯片

图 7 应力时程曲线

Figure 7. Stress-time curves

下载: 全尺寸图片幻灯片

图 8 封闭爆炸试验中爆心距与到时关系

Figure 8. The relationship between blast center distance and arrival time in closed explosion

下载: 全尺寸图片幻灯片

图 9 干燥钙质砂法向应力峰值衰减的试验结果及拟合曲线

Figure 9. Experiment results and fitting curves of the peak normal stress in dry calcareous sand

下载: 全尺寸图片幻灯片

图 10 潮湿钙质砂法向应力峰值试验结果及拟合曲线

Figure 10. Experiment results and fitting curves of the peak normal stress in moist calcareous sand

下载: 全尺寸图片幻灯片

图 11 消除反射波对应力时程曲线影响的方法

Figure 11. The method for eliminating the influence of reflected wave on stress-time curve

下载: 全尺寸图片幻灯片

图 12 DCS64-1试验中各测点的比冲量时程曲线

Figure 12. Specific impulse time histories at each measuring point in DCS64-1 test

下载: 全尺寸图片幻灯片

图 13 干燥钙质砂折合比冲量试验结果及拟合曲线

Figure 13. Experiment results and fitting curves of the specific impulse in dry calcareous sand

下载: 全尺寸图片幻灯片

图 14 潮湿钙质砂折合比冲量试验结果及拟合曲线

Figure 14. Experimental results and fitting curves of the specific impulse in moist calcareous sand

下载: 全尺寸图片幻灯片

图 15 不同含水率钙质砂中的折合比冲量衰减

Figure 15. Attenuation of reduced specific impulse in calcareous sand with different moisture content

下载: 全尺寸图片幻灯片

表 1 影响砂土中爆炸波传播规律的主要参数

Table 1. Major parameters influencing the propagation of blast wave in sand

	参数名称及符号	量纲
炸药	药包质量W	[M]
	装药密度ρ_w	[ML^-3]
	单位质量炸药的化学能E_w	[L²T^-2]
	爆炸产物的膨胀指数γ	—
砂土	密度ρ₀	[ML^-3]
	塑性纵波波速c_p	[LT^-1]
	内摩擦角φ	—
试样尺寸	药包中心埋深h	[L]
试样尺寸	测点到爆心的距离R	[L]

下载: 导出CSV

表 2 爆炸模型试验参数

Table 2. Designed parameters for explosion model test

试验编号	药包设计质量W/g	设计埋深h/m	设计含水率/%	设计爆心距R/m
DCS64-1	64	0.9	1	0.1，0.2，0.3，0.4，0.6，0.8，1.0
DCS64-2	64	0.3	1	0.1，0.2，0.3，0.4，0.6，0.8，1.0
DCS64-3	64	0	1	0.1，0.2，0.3，0.4，0.6，0.8，1.0
DCS216	216	0	1	0.15，0.2，0.3，0.45，0.6，0.8，1.0
DCS512	512	0	1	0.2，0.3，0.4，0.6，0.8，1.0
WCS64-1	64	0.9	10	0.15，0.2，0.3，0.4，0.6，0.8，1.0
WCS64-2	64	0.9	20	0.1，0.2，0.3，0.4，0.6，0.8，1.0
WCS64-3	64	0.9	25	0.1，0.2，0.3，0.4，0.6，0.8，1.0

下载: 导出CSV

表 3 实测爆炸试验参数表

Table 3. Summary of measured parameters in explosion test

试验编号	药包质量/g	试样含水率/%	药包埋深/m	试样密度/(g·cm⁻³)
试验编号	药包质量/g	试样含水率/%	药包埋深/m	总质量除以总体积法	环刀取样法	环刀取土钻取样法
DCS64-1	63.6	0.7±0.5	0.90	1.37	1.35±0.02	−
DCS64-2	63.1	1.5±0.4	0.31	1.39	1.36±0.04	−
DCS64-3	63.3	1.9±0.5	0	1.38	1.36±0.02	−
DCS216	214.2	2.1±0.5	0	1.39	1.37±0.03	−
DCS512	510.2	2.0±0.3	0	1.38	1.37±0.02	−
WCS64-1	62.8	10.5±0.7	0.89	1.54	1.51±0.04	1.54±0.04
WCS64-2	63.3	18.9±1.5	0.92	1.59	1.55±0.07	1.58±0.06
WCS64-3	63.1	23.7±2.2	0.91	1.80	1.76±0.06	1.81±0.04

下载: 导出CSV

表 4 密实钙质砂中的爆炸波波速

Table 4. The blast wave velocity in dense calcareous sand

试验编号	密度/(g·cm⁻³)	弹性纵波波速c₀/(m·s⁻¹)	R²	塑性纵波波速c_p/(m·s⁻¹)	R²
DCS64-1	1.37	359	0.998	250	0.995
DCS64-2	1.39	344	0.998	282	0.995
DCS64-3	1.38	291	0.998	265	0.997
DCS216	1.39	332	0.999	278	0.999
DCS512	1.38	339	0.999	274	0.999
WCS64-1	1.50	429	0.999	377	0.998
WCS64-2	1.59	387	0.998	315	0.996
WCS64-3	1.80	374	0.999	302	0.999

下载: 导出CSV

表 5 干燥钙质砂法向应力峰值衰减规律参数值

Table 5. The parameters of attenuation law of peak normal stress in dry calcareous sand

试验编号	$K_{\sigma}$	$\mu_{\sigma}$	R²	$K_{\sigma}$	$\mu_{\sigma}$	R²
试验编号	第一段(R*≤0.75 m/kg^1/3)			第二段(R*＞0.75 m/kg^1/3)
DCS64-1	0.944	2.94	0.998	1.474	1.37	0.976
DCS64-2	0.944	2.94	0.998	1.474	1.37	0.976
DCS64-3	0.548	3.11	0.961	0.79	1.83	0.938
DCS216
DCS512

下载: 导出CSV

表 6 潮湿钙质砂法向应力峰值衰减规律参数值

Table 6. The parameters for attenuation law of peak normal stress in moist calcareous sand

试验编号	$K_{\sigma}$	$\mu_{\sigma}$	R²
WCS64-1	1.805	1.39	0.965
WCS64-2	1.530	1.55	0.986
WCS64-3	2.200	1.79	0.978

下载: 导出CSV

表 7 非饱和钙质砂刚壁反射因数

Table 7. The rigid reflection coefficient of unsaturated calcareous sand

试验编号	外推自由场应力峰值σ_I,_max/MPa	实测反射波应力峰值σ_R,_max/MPa	刚壁反射因数σ_R,_max/σ_I,_max
DCS64-1	0.421	0.675	1.603
DCS64-2	0.421	0.635	1.508
DCS64-3	0.149	0.219	1.473
DCS216	0.311	0.458	1.473
DCS512	0.529	0.796	1.504
WCS64-1	0.505	0.804	1.592
WCS64-2	0.368	0.451	1.225
WCS64-3	0.430	0.537	1.249

下载: 导出CSV

表 8 非饱和钙质砂中折合比冲量衰减公式参数值

Table 8. The parameters of attenuation law of specific impulse in unsaturated calcareous sand

试验编号	K_I	μ_I	R²
DCS64-1	7.57	1.15	0.944
DCS64-2	5.54	1.14	0.975
DCS64-3	1.57	1.16	0.933
DCS216
DCS512
WCS64-1	5.33	1.18	0.959
WCS64-2	4.92	1.02	0.901
WCS64-3	4.65	0.97	0.946

下载: 导出CSV

参考文献(38)

[1]	刘崇权, 杨志强, 汪稔. 钙质土力学性质研究现状与进展 [J]. 岩土力学, 1995(4): 74–84. DOI: 10.16285/j.rsm.1995.04.010. LIU C Q, YANG Z Q, WANG R. The present condition and development in studies of mechanical properties of calcareous soils [J]. Rock and Soil Mechanics, 1995(4): 74–84. DOI: 10.16285/j.rsm.1995.04.010.
[2]	陈海洋. 钙质砂的内孔隙研究[D]. 武汉: 中国科学院武汉岩土力学研究所, 2005: 1−11.
[3]	陈海洋, 汪稔, 李建国, 等. 钙质砂颗粒的形状分析 [J]. 岩土力学, 2005, 26(9): 1389–1392. DOI: 10.3969/j.issn.1000-7598.2005.09.008. CHEN H Y, WANG R, LI J G, et al. Grain shape analysis of calcareous soil [J]. Rock and Soil Mechanics, 2005, 26(9): 1389–1392. DOI: 10.3969/j.issn.1000-7598.2005.09.008.
[4]	张家铭. 钙质砂基本力学性质及颗粒破碎影响研究[D]. 武汉: 中国科学院武汉岩土力学研究所, 2004: 1−9.
[5]	虞海珍. 复杂应力条件下饱和钙质砂动力特性的试验研究[D]. 武汉: 华中科技大学, 2006: 1−16.
[6]	李建国. 波浪荷载作用下饱和钙质砂动力特性的试验研究[D]. 武汉: 中国科学院武汉岩土力学研究所, 2005: 1−10.
[7]	SOCHET I, GARDEBAS D, CALDERARA S. Fundamental of protective design for conventional weapons: TM5-855-1[R]. Department of The Army Technical Manual, USA Army Corps of Engineers. 1986.
[8]	LYAKHOV G M. Fundamentals of explosion dynamics in soils and liquid media [M]. Moscow: Nedra, 1964.
[9]	LYAKHOV G M, OSADCHENKO R A, POLYAKOVA N I. Plane waves in nonhomogeneous plastic media and their interaction with obstacles [J]. Journal of Applied Mechanics and Technical Physics, 1969, 10(4): 559–566. DOI: 10.1007/bf00916211.
[10]	ZAKHAROV S D, LYAKHOV G M, MIZYAKIN S D. Determination of the dynamic compressibility of soil based on the parameters of plane detonation waves [J]. Journal of Applied Mechanics and Technical Physics, 1972, 13(1): 126–130. DOI: 10.1007/BF00852370.
[11]	LYAKHOV G M, OKHITIN V N. Spherical blast waves in multicomponent media [J]. Journal of Applied Mechanics and Technical Physics, 1974, 15(2): 208–214. DOI: 10.1007/BF00850660.
[12]	LYAKHOV G M, VOVK A A, KRAVETS V G, et al. Compaction of loessal soils by detonation of surface charges [J]. Soil Mechanics and Foundation Engineering, 1976, 13(2): 121–125. DOI: 10.1007/BF01703745.
[13]	LYAKHOV G M. Waves in soils and porous multicomponent media [M]. Nauka, Moscow, 1982.
[14]	LYAKHOV G M, SALITSKAYA V I. Dissipation of blast waves and the dynamic compressibility of soils [J]. Combustion, Explosion and Shock Waves, 1983, 19(1): 90–93. DOI: 10.1007/BF00790244.
[15]	KRYMSKⅡ A V, LYAKHOV G M. Waves from an underground explosion [J]. Journal of Applied Mechanics and Technical Physics, 1984, 25(3): 361–367. DOI: 10.1007/BF00910394.
[16]	LYAKHOV G M, LUCHKO I A, PLAKSⅡ V A, et al. Spherical detonation waves in a solid multicomponent viscoplastic medium [J]. Soviet Applied Mechanics, 1986, 22(5): 490–495. DOI: 10.1007/BF00888551.
[17]	梁霍夫 Г M. 岩土中爆炸动力学基础[M]. 刘光寰, 王明洋, 译. 南京: 工程兵工程学院, 1993.
[18]	KARINSKI Y S, FELDGUN V R, YANKELEVSKY D Z. Effect of soil locking on the cylindrical shock wave’s peak pressure attenuation [J]. Journal of engineering mechanics, 2009, 135(10): 1166–1179. DOI: 10.1061/(ASCE)EM.1943-7889.0000042.
[19]	YANKELEVSKY D Z, KARINSKI Y S, FELDGUN V R. Re-examination of the shock wave’s peak pressure attenuation in soils [J]. International Journal of Impact Engineering, 2011, 38(11): 864–881. DOI: 10.1016/j.ijimpeng.2011.05.011.
[20]	KARINSKI Y S, FELDGUN V R, RACAH E, et al. Mach stem due to an underground explosion near a rigid structure buried in soil [J]. Shock Waves, 2015, 25(1): 63–76. DOI: 10.1007/s00193-014-0544-1.
[21]	徐学勇, 汪稔, 王新志, 等. 饱和钙质砂爆炸响应动力特性试验研究 [J]. 岩土力学, 2012, 33(10): 2953–2959. DOI: 10.16285/j.rsm.2012.10.005. XU X Y, WANG R, WANG X Z, et al. Experimental study of dynamic behavior of saturated calcareous sand due to explosion [J]. Rock and Soil Mechanics, 2012, 33(10): 2953–2959. DOI: 10.16285/j.rsm.2012.10.005.
[22]	徐学勇. 饱和钙质砂爆炸响应动力特性研究[D]. 武汉: 中国科学院武汉岩土力学研究所, 2009: 89−91.
[23]	谢定义. 非饱和土土力学[M]. 北京: 高等教育出版社, 2015.
[24]	MULILIS J P, ARULANANDAN K, MITCHELL J K, et al. Effects of sample preparation on sand liquefaction [J]. Journal of the Geotechnical Engineering Division, 1977, 103(2): 91–108. DOI: 10.1016/0148-9062(77)90060-2.
[25]	LADD R S. Specimen preparation and cyclic stability of sands [J]. Journal of Geotechnical and Geoenvironmental Engineering, 1977, 103(GT-6). DOI: 10.1016/0148-9062(75)91261-9.
[26]	JUANG C H, HOLTZ R D. Fabric, pore size distribution, and permeability of sandy soils [J]. Journal of Geotechnical Engineering, 1986, 112(9): 855–868. DOI: 10.1061/(ASCE)0733-9410(1986)112:9(855).
[27]	NIMMO J R, AKSTIN K C. Hydraulic conductivity of a sandy soil at low water content after compaction by various methods [J]. Soil Science Society of America Journal, 1988, 52(2): 303–310. DOI: 10.2136/sssaj1988.03615995005200020001x.
[28]	亨利奇 J. 爆炸动力学及其应用[M]. 熊建国, 译. 北京: 科学出版社, 1987.
[29]	谢多夫 Л Н. 力学中的相似方法与量纲理论[M]. 沈青, 译. 北京: 科学出版社, 1982.
[30]	马立秋, 张建民. 黏性土爆炸成坑和地冲击传播的离心模型试验研究 [J]. 岩石力学与工程学报, 2011(S1): 3172–3178. MA L Q, ZHANG J M. Centrifugal model testing study of explosion induced craters and propagation of ground shock in clay [J]. Chinese Journal of Rock Mechanics and Engineering, 2011(S1): 3172–3178.
[31]	穆朝民, 任辉启, 辛凯, 等. 变埋深条件下土中爆炸成坑效应 [J]. 解放军理工大学学报(自然科学版), 2010, 11(2): 112–116. DOI: 10.7666/j.issn.1009-3443.20100203. MU C M, REN H Q, XIN K, et al. Effects of crater formed by explosion in soils [J]. Journal of PLA University of Science and Technology (Natural Science Edition), 2010, 11(2): 112–116. DOI: 10.7666/j.issn.1009-3443.20100203.
[32]	施鹏, 邓国强, 杨秀敏, 等. 土中爆炸地冲击能量分布研究 [J]. 爆炸与冲击, 2006, 26(3): 240–244. DOI: 10.11883/1001-1455(2006)03-0240-05. SHI P, DENG G Q, YANG X M, et al. Study on ground shock energy distribution of explosion in soil [J]. Explosion and Shock Waves, 2006, 26(3): 240–244. DOI: 10.11883/1001-1455(2006)03-0240-05.
[33]	叶亚齐, 任辉启, 李永池, 等. 砂质黏土中不同深度爆炸自由场地冲击参数预计方法研究 [J]. 岩石力学与工程学报, 2011, 30(9): 1918–1923. YE Y Q, REN H Q, LI Y C, et al. Study of prediction of ground shock parameters in free field at different depths of burst in sandy clay [J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(9): 1918–1923.
[34]	奥尔连科. 爆炸物理学[M]. 孙承纬, 译. 北京: 科学出版社, 2011.
[35]	水利电力部. 土工试验方法标准: GB/T 50123-2019 [S].
[36]	贾永胜, 王维国, 谢先启, 等. 低含水率砂土和饱和砂土场地爆炸成坑特性实验 [J]. 爆炸与冲击, 2017, 37(5): 799–806. DOI: 10.11883/1001-1455(2017)05-0799-08. JIA Y S, WANG W G, XIE X Q, et al. Characterization of blast-induced craters in low-moisture and saturated sand from field experiments [J]. Explosion and Shock Waves, 2017, 37(5): 799–806. DOI: 10.11883/1001-1455(2017)05-0799-08.
[37]	KRAUTHAMMER T. Modern protective structures [M]. CRC Press, 2008.
[38]	波克罗夫斯基 Г И, 费多罗夫 И С. 在变形介质中冲击与爆破作用[M]. 刘清荣, 黄文彬, 译. 北京: 中国工业出版社, 1965.