爆炸容器（explosive containment vessels，ECVs）是一种能够将爆炸产生的冲击波、破片、有毒气体等限制在有限空间内，密切保护人员和实验仪器，回收爆炸试验碎片，并防止环境污染的限域装置。目前，爆炸容器已被广泛应用于民用和军用领域，特别在民用领域，随着近年来反恐战争形势的日益严峻，在公共场所，尤其是车站和机场，迫切需要轻质、高效的爆炸容器用于处理可疑爆炸物。

爆炸容器内壁爆炸载荷的特征和分布规律是研究容器动态响应特性的基础，同时也是抗爆容器设计的重要参考量。爆炸容器内冲击波的传播是一个极其复杂的过程，不仅涉及到冲击波在空气中的传播、冲击波在容器内部的反射、追赶叠加，同时还涉及到容器壳体响应与冲击波之间的相互作用等多个方面^[1-7]。霍宏发等^[8]通过开展椭球封头圆柱形爆炸容器中心点爆炸试验研究发现，在封闭容器内，爆炸冲击载荷具有明显的脉冲激励和频域上的超带宽特征。饶国宁等^[9]通过数值模拟方法对平板封头圆柱形爆炸容器的内壁爆炸载荷特征进行了研究，发现由于距离爆源最近，在中心爆炸加载条件下，容器的中环面和封头所受首次载荷最大，而由于汇集作用，封头处的后续压力与中环面处明显不同，其后续压力峰值高于第一个反射峰值。张亚军等^[10]通过数值模拟方法对爆炸容器内部冲击波系的演化规律进行了研究，发现对确定尺寸和结构的爆炸容器，当药量不同时，容器内的爆炸场不具有相似性，即药量不同，爆炸所得冲击波系结构和演化过程差别较大。

容器壳体在受内部径向冲击载荷作用下的动态响应特性一直是爆炸容器领域的研究热点，然而对于产生壳体响应的内壁爆炸载荷本身的大小、特征、分布规律及产生机理、载荷和与响应之间的关系等尚缺乏深入系统的研究，相关的试验研究成果更鲜有报道。本文中，通过试验与数值模拟相结合的方法，对圆柱形爆炸容器内炸药爆炸产生的内壁爆炸载荷进行系统的研究，为圆柱形爆炸容器的设计和安全评估提出建议。

1.   试验简介

自行设计的爆炸容器如图1所示，为组合式圆柱形容器，中间圆柱壳部分厚22 mm，内径为0.8 m，长度为0.845 m（不含法兰），并对称开有多个圆孔（含1个起爆孔7个测试孔），开孔处焊接有加强管如图2所示，上下为相同的椭球形端盖，厚度为22 mm，椭球长轴0.8 m，长轴与短轴长度之比为2∶1，下端盖与圆柱壳部分焊接为一整体，上端盖与圆柱壳部分通过法兰上的螺栓连接。

图  1  爆炸容器

Figure  1.  Explosion containment vessel

下载: 全尺寸图片 幻灯片

图  2  测点分布

Figure  2.  Distribution of measurement points

下载: 全尺寸图片 幻灯片

为研究装药在圆柱形爆炸容器内部爆炸时容器内壁所受载荷的特征及分布规律，在容器中选择了P₁～P₅共5个典型位置作为内壁反射超压测点，如图2所示。P₁测点位于装药中心所在的容器中心横截面上，P₂、P₃、P₄测点与P₁具有相同的环向坐标，在轴向上与P₁分别相距150、300、450 cm，P₅测点位于椭圆形下端盖的极点处。测试系统由PCB公司的CA102B03型压力传感器、低噪声电缆和Elsys公司的TraNET PPC型一体式便携高速数据采集仪构成。

为适应容器壳体为薄壁曲面的结构特点，消除内壁面超压测试过程中容器壳体的振动等对测试信号的寄生干扰，同时满足爆炸容器的密封要求，本文中对传感器的安装方式进行了设计，如图3所示。由图2和图3(a)可以看出，加强管焊接在容器壁开孔位置，传感器通过自身的螺纹固定于基座上，传感器与基座之间放置的紫铜垫片主要起密封作用，防止容器内部的高压气体通过传感器和基座之间的螺纹间隙发生泄漏。基座通过后端的压紧螺栓固定于加强管内，橡胶垫圈可以有效地防止容器内部的高压气体通过基座与加强管之间的间隙发生泄漏。在环向上，基座与加强管之间具有微间隙，并隔有O型橡胶圈；在轴向上，基座与加强管之间隔了一层橡胶垫片，因此可以有效地抑制容器壳体的振动对反射冲击波信号的干扰，同时传感器本身具有加速度补偿功能，也能很好地抑制冲击和壳体振动等带来的噪声。图3(b)为传感器安装在基座上的实物，传感器敏感面相比基座端面凹陷约为0.6 mm，实际测试时，在传感器敏感面涂上一层约0.6 mm厚的耐高温油脂，以阻止爆炸瞬间产生的高温等对冲击波信号的干扰，而基座端面与容器内壁面平齐。

图  3  传感器安装设计

Figure  3.  Design of the sensor installation

下载: 全尺寸图片 幻灯片

为确保容器在实验过程中不会发生塑性变形，采用逐步加载的方法，试验过程中，依次使用质量约27、64、100和150 g的球形压装TNT装药（不包括传爆药质量）进行试验，每种工况至少进行2次试验以确保测试结果的可靠性。

2.   试验结果

2.1   内壁爆炸载荷特征及分布规律

图4展示了64 g规格球形装药在容器中心位置爆炸时两次重复试验的实测结果。图5展示了64 g装药容器中心爆炸时各测点所受载荷的完整时间历程曲线。由图4可知：（1）2次试验中，各测点的压力时程曲线一致性非常好，表明试验可重复性好，测试结果准确可靠；（2）容器内壁所受爆炸载荷具有明显的多脉冲特征，随着测点位置的变化，脉冲载荷的到达时间和大小表现出明显差异，即载荷在容器内壁的分布具有明显的不均匀性；（3）圆柱部分与端盖极点附近容器壳体所受载荷存在较大差异。

图  4  装药为64 g的容器中心爆炸时各测点所受载荷压力曲线

Figure  4.  Pressure curves of blast loading on each measurement point generated by a spherical charge of 64 g in the center of the container

下载: 全尺寸图片 幻灯片

由图5可知：（1）从整体上看，容器内壁爆炸载荷大致可分为3个作用阶段，即初始冲击阶段Ⅰ（第1个压力脉冲载荷）、冲击波在容器内壁与中心之间来回反射的冲击反射阶段Ⅱ以及最后容器内部压力趋于均匀的准静态阶段Ⅲ；（2）对于圆柱壳体部分，初始冲击阶段的持续时间在0.6～0.7 ms之间，随着与容器中心轴向距离的增大，大致呈减小趋势，对于椭球端盖极点，初始冲击阶段的持续时间约为0.22 ms；（3）对整个容器而言，从最早的第1个脉冲载荷到达容器壳体，即到达P₁测点算起，经历约21 ms，容器内部压力整体趋于均匀，各测点压力虽仍处于震荡状态，但震荡幅度较小。

图  5  典型冲击波压力曲线（第10次试验结果）

Figure  5.  Typical pressure curves of shock waves (results of test 10)

下载: 全尺寸图片 幻灯片

将压力时程曲线上具有相对完整且较为显著的上升和下降阶段的一次作用周期看作一个压力脉冲。由图5可知：（1）在容器中建立相对均匀的压力分布之前，P₁～P₄测点约受到3～4次脉冲作用，P₅测点则受到多次脉冲作用，在初次脉冲作用后，对圆柱壳体部分，反射脉冲作用持续时间在1.2～2.4 ms之间，对椭球端盖极点位置，反射脉冲作用持续时间约10.7 ms；（2）对整个容器而言，第3次压力脉冲后，距离容器中心越远的测点，所测载荷曲线震荡越剧烈，特别是对于端盖极点位置，后续反射脉冲的峰值压力出现与圆柱壳部分所受峰值压力相当的情况，这表明越靠近容器中心，壳体所受载荷越早趋于均匀；（3）对于P₁～P₄测点，即容器圆柱壳部分，在经过1次反射之后，壳体所受脉冲载荷峰值压力整体呈衰减趋势，且衰减迅速，从第2个脉冲开始，载荷峰值压力基本在1 MPa以下，这与Duffey等^[11]的研究结论一致，第1个脉冲具有比后续脉冲高得多的峰值压力，同时圆柱壳体不同部位所受脉冲载荷峰值压力之间的差异大幅减小，壳体受载荷作用相对均匀；（4）对于P₅测点来说，即椭球端盖极点位置，第3个脉冲的峰值超压是最高的，这表明与圆柱壳体部分不同，椭球端盖极点位置所受峰值载荷并不由装药爆炸的直达冲击波作用产生。

表1为不同装药条件下容器内壁各测点所受载荷峰值压力的统计结果（对P₁～P₄测点，即首脉冲的峰值压力），表中：装药位置“O”代表容器中心，经验计算值为根据林俊德公式（适用范围为比例距离$\bar{r}$满足条件${0.2}\;\text{m/}{\text{kg}}^{{1/3}}{\text{≤}}\bar{{r}}{\text{≤}3}\;\text{m/}{\text{kg}}^{{1/3}}$）计算而来^[12]；表2为测点各次脉冲的比冲量统计，表中“1^st”、“2^nd”等分别代表脉冲到达的顺序，空格表示没有测到有效数据，比冲量值为单次脉冲正压段压力对时间的积分；表3为不同装药条件下容器内壁P₁～P₄测点所受载荷首脉冲正压作用时间统计，表中空格表示正压作用时间无法判读。由表1可知，容器内壁各测点所受爆炸载荷峰值压力实测结果与依据林俊德公式求解的经验计算值之间存在一定误差。为提高密闭爆炸容器内壁爆炸载荷经验计算公式的求解精度，根据表1中P₁～P₄测点首脉冲峰值压力实测结果对计算公式进行拟合。对于冲击波峰值超压${\Delta }{p}$（单位为MPa），根据经典指数形式公式^[3]，设：

表  1  峰值压力试验结果

Table  1.  Test results of the peak pressure

试验 编号	装药 质量/g	装药 位置	峰值压力/MPa
			P1测点				P2测点				P3测点				P4测点				P5测点
			实测	计算	误差		实测	计算	误差		实测	计算	误差		实测	计算	误差		实测	计算	误差
1	26.7	O	2.64	2.11	0.25		1.83	1.78	0.03		1.08	1.18	0.09		0.61	0.73	0.16		2.72	0.34	7.05
2	27.1	O	2.36	2.14	0.10		1.56	1.80	0.13		0.97	1.20	0.19		0.61	0.74	0.17		2.89	0.34	7.44
3	26.7	O	2.41	2.11	0.14		1.61	1.78	0.09		0.98	1.18	0.17		0.59	0.73	0.19		3.01	0.34	7.91
6	63.7	O	5.33	4.48	0.19		4.71	3.78	0.25		2.58	2.51	0.03		1.21	1.55	0.22		6.93	0.72	8.65
7	64.8	O	5.42	4.55	0.19		4.83	3.83	0.26		2.77	2.55	0.09		1.23	1.57	0.22		14.48	0.73	18.87
8	64.4	O	5.76	4.52	0.27		4.80	3.81	0.26		2.33	2.53	0.08		1.20	1.56	0.23		9.18	0.73	11.67
9	64.7	O	5.11	4.54	0.13		4.75	3.83	0.24		2.68	2.54	0.05		1.13	1.57	0.28		14.26	0.73	18.59
10	64.1	O	5.58	4.51	0.24		4.61	3.80	0.21		2.66	2.52	0.05		1.29	1.56	0.17		13.41	0.72	17.58
15	105.7	O	7.34	6.95	0.06		5.93	5.86	0.01		4.34	3.89	0.12		2.30	2.40	0.04		19.07	1.11	16.20
16	105.2	O	7.73	6.92	0.12		5.87	5.83	0.01		4.26	3.88	0.10		1.93	2.39	0.19			1.11
21	154.5	O	11.44	9.66	0.18		8.70	8.14	0.07		6.50	5.41	0.20		3.13	3.34	0.06		13.77	1.55	7.90
22	154.3	O	11.49	9.65	0.19		9.10	8.13	0.12		6.14	5.40	0.14		3.33	3.33	0		13.69	1.55	7.86
23	155.2	O	11.17	9.70	0.15		10.28	8.17	0.26		5.99	5.43	0.10		2.85	3.35	0.15		18.30	1.55	10.78
24	154.2	O	10.88	9.64	0.13		9.92	8.13	0.22		5.86	5.40	0.09		2.91	3.33	0.13			1.55
25	154.7	O	11.52	9.67	0.19		9.23	8.15	0.13		5.76	5.41	0.06		2.91	3.34	0.13			1.55

下载: 导出CSV 

| 显示表格

表  2  测点各次脉冲的比冲量

Table  2.  Specific impulses of each pulse at each measurement points

试验 编号	装药 质量/g	比冲量/（Pa∙s）
		P1测点				P2测点					P3测点					P4测点					P5测点
		1st	2nd	3rd		1st	2nd	3rd	4th		1st	2nd	3rd	4th		1st	2nd	3rd	4th		1st	2nd	3rd	4th	5th	6th
1	26.7	156	48	127		138	51	65			107	75	123			89	117	115			53	393	95	62	536	65
2	27.1	130	50	80		143	35	81			106	61	132			91	94	125			56	426	116	78	540	93
3	26.7	150	64	95		120	64	61			109	69	119			83	83	105			55	394	112	42	495	69
6	63.7	267	92	200		244	108	156			221	107	203	235		155	151	139	152		104	717	304	163	862	342
7	64.8	299	166	206		227	122	149			215	126	223	213		146	127	140	138		88	577	250	88	777	290
8	64.4	295	175	214		256	131	154			199	116	197	253		157	129	147	132		79	688	180	97	758	359
9	64.7	247	154	246		265	129	175			219	107	240	213		147	114	166	147		77	620	150	34	633	270
10	64.1	292	156	212		244	121	154			190	112	210	241		147	120	153	143		69	604	218	82	674	270
15	105.7	418	223	476		379	194	169			297	142	178	281		206	70	103	179		339	705	199	218	883	339
16	105.2	410	263	495		350	185	155			286	156	257	305		184	88	77	162
21	154.5	437	211	571		485	187	158	253		356	149	281	180		265	63	127	297		159	459	795	462	1084	357
22	154.3	593	288	559		517	211	158	290		412	161	205			280	62	141	280		194	405	866	524	1071	445
23	155.2	494	250	577		460	195	143	270		367	152	280	175		249	85	97	279		162	429	713	492	1054	234
24	154.2	630	327	707		474	225	142	335		367	189	308			248	107	81	256
25	154.7	557	294	624		461	225	162	321		357	170	297	189		236	90	88	246

下载: 导出CSV 

| 显示表格

表  3  测点首脉冲正压作用时间

Table  3.  Positive pressures action time of the first pulse at each measurement points

试验 编号	正压作用时间/ms
	P1测点				P2测点				P3测点				P4测点
	试验	计算	误差		试验	计算	误差		试验	计算	误差		试验	计算	误差
1	0.349	0.376	0.072		0.361	0.396	0.089		0.439	0.449	0.022		0.507	0.520	0.024
2	0.345	0.377	0.084		0.402	0.397	0.013		0.465	0.449	0.035		0.552	0.520	0.061
3	0.357	0.376	0.051		0.411	0.396	0.037		0.453	0.449	0.009		0.538	0.520	0.035
6	0.404	0.400	0.010		0.436	0.421	0.035		0.439	0.477	0.080		0.584	0.552	0.057
7	0.406	0.400	0.014		0.412	0.422	0.023		0.516	0.478	0.081		0.572	0.553	0.034
8	0.403	0.400	0.007		0.439	0.422	0.041		0.475	0.477	0.005		0.580	0.553	0.049
9	0.411	0.400	0.027		0.392	0.422	0.070		0.483	0.477	0.012		0.535	0.553	0.032
10	0.397	0.400	0.008		0.388	0.421	0.079		0.448	0.477	0.061		0.553	0.553	0.001
15	0.420	0.414	0.014		0.436	0.436	0.001		0.478	0.494	0.033			0.572
16	0.410	0.414	0.010		0.424	0.436	0.028		0.496	0.494	0.004			0.572
21	0.420	0.425	0.013		0.411	0.448	0.083		0.521	0.507	0.027			0.588
22	0.455	0.425	0.070		0.418	0.448	0.067		0.528	0.507	0.041			0.588
23	0.438	0.426	0.029		0.391	0.448	0.128		0.482	0.508	0.050			0.588
24	0.432	0.425	0.016		0.436	0.448	0.027		0.506	0.507	0.003			0.588
25	0.414	0.425	0.027		0.439	0.448	0.021		0.506	0.508	0.003			0.588

下载: 导出CSV 

| 显示表格

式中：k、α为待定系数；r为壳体与装药间的距离，m；m为装药质量，kg。将式(1)整理成如下形式：

令${y=}{\ln}{\Delta }{p}$，$x=\ln\left({r}/\sqrt[{3}]{{m}}\right)$，$b={\ln}{k}$，则有${y=\alpha{x}+b}$。利用Matlab自带的一元线性回归分析工具对表1中的试验结果进行拟合，图6为初步拟合结果。由图6可知，试验结果中存在奇异点，将奇异数据点剔除，重新进行拟合，直到残差图中奇异数据点不再出现，最终拟合结果如图7所示。此时，最终回归模型的回归系数$\alpha =-{2.923\;4}$，常数${b}{=1.610\;2}$，判定系数${{{R}}_{1}^{2}}{=0.987\;3}$，统计量值${F}{=7\;532}$，与F对应的概率${P}{=9.307\,4}{\,\rm{e}}^{-{94}}$，估计标准误差${{S}}_{{y}}{=0.006\;4}$。由于概率P＜0.05，判定系数非常接近1且估计标准误差较小，因此可以认为回归方程非常显著且拟合效果较好。因此有：

图  6  测点P₁～P₄首脉冲峰值压力初步拟合结果

Figure  6.  Preliminary fitting results of shock wave peak pressure of the first pulse from measurement points P₁–P₄

下载: 全尺寸图片 幻灯片

图  7  测点P₁～P₄首脉冲峰值压力最终拟合结果

Figure  7.  Final fitting results of shock wave peak pressures of the first pulse from measurement points P₁–P₄

下载: 全尺寸图片 幻灯片

则拟合的冲击波反射超压峰值经验计算公式为：

根据表2所示试验结果，使用相同的方法拟合可得P₁～P₄测点首脉冲比冲量计算公式为：

由表3可知，P₁～P₄测点首脉冲正压作用时间的林俊德公式计算结果与试验结果吻合较好，因此正压作用时间的计算公式为：

由表1～3可知：（1）对容器圆柱壳体部分而言，P₁～P₄测点首脉冲正压作用时间与测点距爆源的距离以及装药量都成正相关，测点距爆源的距离对正压作用时间的影响稍大于装药量；（2）离爆源最近的P₁测点所受载荷峰值压力及单次脉冲比冲量峰值最大，P₁～P₄测点比冲量的峰值并不总是在第1个脉冲出现，由于作用时间更长，后续第3或第4个脉冲的比冲量常常与第1个脉冲的比冲量基本相当甚至更大，又反射脉冲作用持续时间在1.2～2.4 ms之间，因此在确定容器圆柱壳体部分响应时不能简单地只考虑首脉冲的作用；（3）对整个容器而言，P₅测点的峰值压力及单次脉冲比冲量峰值总是所有测点中最大的，且其峰值压力远大于林俊德公式计算值。将P₅测点所受载荷峰值压力（p₅）与相同条件下容器圆柱壳部分所受载荷峰值压力即爆心所在环面P₁测点峰值压力（p₁）进行对比，结果如表4所示。由表4可知：（1）P₅测点所受载荷峰值压力测量结果的一致性并不理想，这主要与装药位置的设置误差以及椭球端盖极点附近内壁的不平整有关；（2）椭球端盖极点位置所受载荷峰值压力与容器圆柱壳部分所受载荷峰值压力的比值在1.03～2.79之间，因此在容器设计及安全评估时，端盖极点位置应予以足够重视。

表  4  椭球端盖所受载荷峰值压力与圆柱壳体部分所受载荷峰值压力的比较

Table  4.  Comparison between the peak pressure of the load on the ellipsoid end cap and the peak pressure of the load on the cylindrical shell

试验编号	装药质量/g	装药位置	试验/MPa		p5/p1		试验编号	装药质量/g	装药位置	试验/MPa		p5/p1
			p5	p1						p5	p1
1	26.7	O	2.72	2.64	1.03		15	105.7	O	19.07	7.34	2.47
2	27.1	O	2.89	2.36	1.23		16	105.2	O		7.73
3	26.7	O	3.01	2.41	1.25		21	154.5	O	13.77	11.44	1.20
6	63.7	O	6.93	5.33	1.30		22	154.3	O	13.69	11.49	1.19
7	64.8	O	14.48	5.42	2.67		23	155.2	O	18.30	11.17	1.64
8	64.4	O	9.18	5.76	1.59		24	154.2	O		10.88
9	64.7	O	14.26	5.11	2.79		25	154.7	O		11.52
10	64.1	O	13.41	5.58	2.40

下载: 导出CSV 

| 显示表格

图8展示了装药在容器中心爆炸时容器内壁所受载荷峰值压力随壳体与装药间轴向距离的变化。由图8可知：（1）对于容器圆柱壳体部分，随着与装药间轴向距离的增大，内壁所受载荷峰值压力逐渐减小，且随着装药质量的增大，峰值压力衰减的速率加快；（2）从整体来看，椭球端盖极点所受载荷峰值压力总是大于圆柱壳部分，且随着装药质量的增大，两者之间的差距呈先增大后减小的趋势。

图  8  载荷峰值压力的分布

Figure  8.  Distributions of the peak pressures of the blast loading

下载: 全尺寸图片 幻灯片

表5给出了不同平均装药质量下准静态阶段容器内壁所受准静态压力的平均值统计结果。基于理想气体假设，容器内部准静态压力p为：

表  5  容器内部准静态压力平均值

Table  5.  Average quasi-static pressures inside the vessel

平均装药质量/g	准静态平均压力/MPa
26.74	0.16
64.04	0.36
105.05	0.54
154.59	0.65

下载: 导出CSV 

| 显示表格

式中：${\rho}$为气体密度；R为理想气体参数，${R=}{{R}}_{{0}}{/M}$，R₀为普适气体参数，M为气体摩尔质量；T为气体绝对温度。

假设炸药爆炸以及容器内气体的温升是一个绝热过程，并且忽略容器内原有气体的质量以及假设炸药爆炸释放的能量全部用于容器内气体的温升，因此容器内气体的温升为：

式中：E为炸药爆炸释放的能量，Q为装药的TNT当量，c_V为气体的定容比热容^[5]。

容器内的准静态气体压力为：

化简可得：

式中：E₀为单位质量TNT爆炸释放的能量，T₀为气体初始温度即室温，V₀为容器体积。由于对于固定的炸药，R₀、E₀、M、c_V、T₀均为常数，因此式(10)可化简成如下形式：

式中：$K=\dfrac{{R}_{0}{E}_{0}}{M{c}_{V}}+\dfrac{{R}_{0}{T}_{0}}{M}$。

上述分析表明，准静态阶段容器内壁所受准静态压力近似与装药量成正比，与容器体积成反比。对表5的试验结果进行拟合，本文中圆柱形爆炸容器的准静态压力有：

2.2   机理分析

为了说明圆柱形爆炸容器内壁所受载荷的特征和分布规律的形成机理，考虑到圆柱形爆炸容器相对于过P₁测点的截面上下近似对称，图9给出了容器下半部分冲击波传播路径的示意图。同时利用LS-DYNA软件建立了容器1/8简化数值模型，如图10所示。图11为150 g球型TNT炸药在容器中心爆炸时数值模拟所得不同时刻容器内部爆轰场的压力云图。

图  9  冲击波传播路径示意图

Figure  9.  Schematic diagram of shock wave propagation path

下载: 全尺寸图片 幻灯片

图  10  容器1/8简化数值模型

Figure  10.  1/8 simplified numerical model of the cylindrical explosion containment vessel

下载: 全尺寸图片 幻灯片

图  11  容器内部爆轰场压力云图

Figure  11.  Pressure cloud map of the internal detonation field inside the vessel

下载: 全尺寸图片 幻灯片

由图9可知，装药爆炸后爆轰波向四周传播形成冲击波，一定时间后冲击波到达容器内壁面并同时开始发生反射，对应压力时程曲线上的初始冲击阶段；反射冲击波向容器内部传播，如图中路径1所示，冲击波到达P₁测点后发生正反射，反射波沿原路向容器中心传播，一定时间后与从容器内壁相同环面反射的冲击波在容器中心位置相遇，随后改变方向，沿相反方向继续传播，如此重复；对于圆柱壳上的其他测点P₂、P₃、P₄，类似图9中路径4所示，爆炸产生的直达冲击波首先到达测点，之后如路径2所示，从容器内壁其他位置反射的冲击波在容器内部相遇后再次向外反射到达测点，同时从容器内壁其他某些位置反射的冲击波在内壁、容器内部和椭球端盖之间来回的反射最后到达测点，反映在压力时程曲线上即表现为多个脉冲；由于随着传播距离的逐渐增大，冲击波逐渐衰减，因此随着与爆源之间轴向距离的增加，内壁所受载荷峰值压力（即首脉冲峰值压力）逐渐减小，P₁～P₄测点的后续几个脉冲峰值压力逐渐衰减，且明显小于第1个脉冲的峰值压力。

对于椭球端盖极点位置的P₅测点，如图9中路径8所示，爆炸产生的直达冲击波首先到达测点，之后如路径6、7所示，爆炸产生的第1道波与从容器内壁反射的冲击波在容器内部相遇并发生复杂的相互作用，如图11(a)所示。相遇后冲击波压力逐渐升高，当传播到椭球端盖时（图9中A₁、A₂位置），由于其入射角度大于马赫反射的临界角，因此会在端盖内壁面发生马赫反射，此时峰值压力再次增大，如图11(b)所示；马赫波沿端盖内壁传播同时逐渐向极点汇聚，如图11(c)、(d)所示，并最后在端盖极点附近完成汇聚，如图11(e)所示。汇聚后的冲击波会再次向端盖四周反射，如此反复，直到最后容器内部达到一种准静止状态。因此在P₅测点的压力时程曲线上出现了许多个脉冲，且在冲击波完成汇聚时刻测点的压力达到最大值。

3.   结　论

通过对圆柱形爆炸容器内壁爆炸载荷的试验和数值模拟研究，可以得到以下结论。

（1）容器内壁所受爆炸载荷具有明显的多脉冲特征，从整体上大致可分为3个作用阶段，即初始冲击阶段、冲击波在容器内壁和中心之间来回反射的冲击反射阶段以及最后容器内部压力趋于均匀的准静态阶段；对圆柱壳体部分，初始冲击阶段的持续时间在0.6～0.7 ms之间，对椭球端盖极点，初始冲击阶段的持续时间约为0.22 ms；对整个容器而言，经历约21 ms后，容器内部压力整体趋于均匀。

（2）在内部建立均匀分布的压力之前，容器圆柱壳上的测点所受载荷包含3～4个显著的完整脉冲，初次脉冲作用后，反射脉冲作用持续时间在1.2～2.4 ms之间，且后续第3或第4个脉冲的比冲量常常与第1个脉冲的比冲量基本相当甚至更大，因此在确定容器圆柱壳体部分响应时不能简单地只考虑首脉冲的作用；椭球端盖极点处的P₅测点受到多次脉冲的作用，在初次脉冲作用后，反射脉冲作用持续时间约10.7 ms，测点峰值压力及单次脉冲比冲量峰值总是所有测点中最大的，峰值压力最高可达圆柱壳所受最大压力的2.79倍，应予以足够重视。

（3）对整个容器而言，P₅测点所受载荷峰值压力最大，数值模拟表明，产生这种现象的主要原因是从容器内壁反射的冲击波与爆炸产生的第1道波在经过复杂的相互作用后于端盖内壁发生了马赫反射，马赫波沿端盖内壁传播并最终在端盖极点汇聚。

（4）试验得到当比例距离为${0.74\text{≤}}{r}/\sqrt[{3}]{{m}}{\text{≤}2.70}$时，容器圆柱壳体所受载荷首脉冲的峰值压力、正压作用时间和比冲量的经验计算公式分别为${\Delta }{p}{=5}{\left({r}/\sqrt[{3}]{{m}}\right)}^{{-2.92}}$，${\Delta }{T}{=}{\left({r}/\sqrt[{3}]{{m}}\right)}^{{0.79}}{{m}}^{{1}/{3}}$，$I{=658.52}{\left({r}/\sqrt[{3}]{{m}}\right)}^{{-1.3}}{{m}}^{{1}/{3}}$；容器内的准静态气体压力经验计算公式为${{p}}_{\rm{s}}{=3.468}\dfrac{{m}}{{{V}}_{{0}}}$、${0.036\text{≤}}\dfrac{{m}}{{{V}}_{{0}}}{\text{≤}0.207}$。