近场近地爆炸下建筑柱爆炸荷载分布规律及简化模型

喻君; 刘福余; 方秦

doi:10.11883/bzycj-2022-0366

近场近地爆炸下建筑柱爆炸荷载分布规律及简化模型

doi: 10.11883/bzycj-2022-0366

喻君^{1, 2,},
刘福余^{2, 3},
方秦⁴

1.
东南大学爆炸安全防护教育部工程研究中心，江苏南京 211189
2.
河海大学土木与交通学院，江苏南京 210098
3.
浙江精工钢结构集团有限公司，浙江绍兴 312030
4.
中国人民解放军陆军工程大学，江苏南京 210001

基金项目: 国家自然科学基金（51978246）；中央高校基本科研业务费与专项资金（2242022R10073）

详细信息

作者简介:
喻　君（1982－　），男，博士，教授，junyu@seu.edu.cn

中图分类号: O383; O384
计量
- 文章访问数: 642
- HTML全文浏览量: 104
- PDF下载量: 396
- 被引次数: 0
出版历程
- 收稿日期: 2022-07-02
- 修回日期: 2022-10-08
- 网络出版日期: 2022-10-14
- 刊出日期: 2024-01-11

Distribution pattern and simplified model of blast load for building columns under near-field near-ground explosion

YU Jun^{1, 2
,},
LIU Fuyu^{2, 3},
FANG Qin⁴

1.
Engineering Research Center of Safety and Protection of Explosion & Impact of Ministry of Education, Southeast University, Nanjing 211189, Jiangsu, China
2.
College of Civil and Transportation Engineering, Hohai University, Nanjing 210098, Jiangsu, China
3.
Jinggong Steel Building Group, Shaoxing, 312030, Zhejiang, China
4.
Army Engineering University of PLA, Nanjing 210001, Jiangsu, China

摘要

摘要: 为了快速评估近场近地爆炸荷载下建筑柱的动力响应和破坏模式，通过数值仿真方法，探究了近场近地爆炸工况下冲击波在建筑柱迎爆面的分布规律，并提供了该工况下的爆炸荷载简化模型。为此，首先利用已有实验数据验证数值模型，并建立典型近地近场爆炸工况的数值模型，然后研究比例爆距和比例爆高对建筑柱冲击波特征参数的影响规律，最后拟合出柱迎爆面反射冲量和正相超压持续时间的计算公式，将柱迎爆面各点爆炸荷载转化为等效三角形荷载模型，为工程实践中建筑柱遭受近场近地爆炸作用下的抗爆设计提供荷载输入。研究结果表明：当比例爆高小于0.3 m/kg^1/3、比例爆距在0.4～0.6 m/kg^1/3范围时，最大反射冲量沿柱高可简化为三折线分布；当比例爆距在0.6～1.4 m/kg^1/3范围时，最大反射冲量沿柱高可近似简化为双折线分布；在同一比例爆距和比例爆高工况下，随着炸药当量的增加，柱迎爆面相同比例高度处反射超压峰值保持不变而反射冲量正比于当量的立方根。
- 建筑柱 /
- 近场爆炸 /
- 近地爆炸 /
- 爆炸荷载 /
- 荷载模型
Abstract: To rapidly assess the dynamic responses and failure modes of the building columns under near-field near-ground explosions, in this paper numerical simulation method is employed to investigate the distribution pattern of the shock waves that are applied on the front face of building columns under near-field near-ground blast scenarios, and a corresponding simplified blast load model is proposed. To this end, firstly, the existing experimental data of overpressure and impulse were selected to validate the numerical model for blast load. Then, a typical numerical model under near-field near-ground blast scenarios was established to study the effects of the scaled distance and the scaled height of spherical charges on the characteristic values of the shock waves acting at the building columns. Finally, formulae for the maximum reflected impulse and the representative value of the positive overpressure duration were derived based on nonlinear regression analysis, and the blast load at each location of the column front face was represented by an equivalent triangular load model. The results indicate that when the scaled height of the charge is less than 0.3 m/kg^1/3, the distribution of the maximum reflected impulse along the column length can be represented as a trilinear model and a bilinear model for the scaled distance of 0.4−0.6 m/kg^1/3 and 0.6−1.4 m/kg^1/3, respectively. In comparison, the distribution of the shock waves in the transverse direction of a column section was approximately uniform. Moreover, under a given scaled distance and a scaled height, the peak reflected overpressure remains constant as the charge weight increases, but the maximum reflected impulse is proportional to the cubic root of the charge weight at the locations with the identical scaled height of the column.
- building column /
- near-field blast /
- near-ground blast /
- blast load /
- load model

HTML全文

图 1 实验设置

Figure 1. Experimental set-up

下载: 全尺寸图片幻灯片

图 2 参照实验对应的有限元模型

Figure 2. Finite element models for reference tests

下载: 全尺寸图片幻灯片

图 3 近场近地球形炸药爆炸荷载数值模型（单位：mm）

Figure 3. Numerical model of blast load generated by spherical charges under near-field near-ground scenarios (unit: mm)

下载: 全尺寸图片幻灯片

图 4 基于柱表面反射超压的数值模型参数敏感性分析（η为网格尺寸大小）

Figure 4. Sensitivity analysis of numerical model parameters in accordance with overpressure at the front face of the columns (η is the element size)

下载: 全尺寸图片幻灯片

图 5 沿柱高的冲击波传播过程（底部传播情况）

Figure 5. Shock wave propagation along the column length (propogation at the bottom of the column)

下载: 全尺寸图片幻灯片

图 6 沿柱高的冲击波传播过程（中部传播情况）

Figure 6. Shock wave propagation along the column length (propogation at the middle of the column)

下载: 全尺寸图片幻灯片

图 7 沿柱高的冲击波传播过程（顶部传播情况）

Figure 7. Shock wave propagation along the column length (propogation near the top of the column)

下载: 全尺寸图片幻灯片

图 8 柱迎爆面最大反射冲量分布

Figure 8. Distribution of the maximum reflected impulse at the front face of the column

下载: 全尺寸图片幻灯片

图 9 最大反射冲量分布简化模型

Figure 9. Simplified models of the maximum reflected impulse distribution

下载: 全尺寸图片幻灯片

图 10 关键冲量的拟合结果

Figure 10. Fitting results of critical impulses

下载: 全尺寸图片幻灯片

图 11 冲量沿柱截面宽度的横向分布简化

Figure 11. Simplification of impulse distribution in transverse direction along the column width

下载: 全尺寸图片幻灯片

图 12 正相持续时间t₀与最大反射冲量沿柱高的分布及t₀的拟合结果

Figure 12. Distribution of the positive phase duration (t₀) and the maximum reflected impulse (I_r) along the colume height and the fitting result of t₀ at the detonation height

下载: 全尺寸图片幻灯片

图 13 不同柱高处反射超压时程曲线

Figure 13. Time-histories of reflected overpressure at different column heights

下载: 全尺寸图片幻灯片

图 14 炸药相似比对超压、冲量和超压时程的影响（Z=0.6 m/kg^1/3, H_c/W^1/3=0.1 m/kg^1/3）

Figure 14. Influence of λ on peak overpressure, maximum reflected impulse and time history of overpressure (Z=0.6 m/kg^1/3, H_c/W^1/3=0.1 m/kg^1/3)

下载: 全尺寸图片幻灯片

图 15 炸药相似比对超压、冲量和超压时程的影响（Z=0.6 m/kg^1/3, H_c/W^1/3=0.3 m/kg^1/3）

Figure 15. Influence of λ on overpressure, impulse and overpressure time history（Z=0.6 m/kg^1/3, H_c/W^1/3=0.3 m/kg^1/3）

下载: 全尺寸图片幻灯片

图 16 双折线模型对比验证

Figure 16. Verification of the bilinear model by the maximum reflected impulse

下载: 全尺寸图片幻灯片

图 17 三折线模型对比验证

Figure 17. Verification of the trilinear model by the maximum reflected overpressure

下载: 全尺寸图片幻灯片

表 1 Woodson实验工况表^[18]

Table 1. Woodson’s experimental cases^[18]

工况	结构	m_C4/kg	m_TNT/kg	R/m	Z/(m·kg⁻³)
实验1	无填充墙	7.1	9.6	1.52	0.72
实验2	无填充墙	7.1	9.6	1.07	0.50
实验3	全填充墙	7.1	9.6	1.07	0.50
实验4	含窗填充墙	7.1	9.6	1.07	0.50
实验5	车库	7.1	9.6	1.07	0.50
注：m为质量，C4是实验所使用炸药，模拟使用TNT炸药，当量转化系数取1.35；Z为比例爆距；R为爆距。

下载: 导出CSV

表 2 Liu实验工况^[19]

Table 2. Liu’s experimental cases^[19]

工况	TNT/kg	R/m	H_c/m	Z/(m·kg⁻³)	R_so/m
1	1	3	1.05	3.00	2
2	2	3	1.05	2.38	2
3	4	2	1.05	1.89	2
注：R为炸药中心到柱等高位置处的爆距；H_c为爆心距地表高度；R_so为自由场超压传感器布置距离；Z为比例爆距。

下载: 导出CSV

表 3 Woodson实验^[18]爆炸荷载的模拟结果与实测值对比

Table 3. Comparisons of numerical values and Woodson’s experimental^[18] results

工况	测点	p_r/MPa		I_r/(kPa∙s)		误差/%		工况	测点	p_r/MPa		I_r/(kPa∙s)		误差/%
工况	测点	实测^[18]	模拟	实测^[18]	模拟	p_r	I_r	工况	测点	实测^[18]	模拟	实测^[18]	模拟	p_r	I_r
实验1	BC1	8.32	12.17	1.04	1.10	46.27	5.77	实验3	BC1	15.05	16.61	3.02	3.43	10.37	13.58
	BC2	0.40	0.56	0.24	0.19	40.00	−20.83		BC2
	BC3	9.36	6.23	1.22	1.02	−33.44	−16.39		BC3	9.57	15.67	1.71	2.28	63.74	33.33
实验2	BC1	11.52	12.53	2.77	2.38	8.77	−14.08	实验5	BC1	11.52	12.53	2.33	2.38	8.77	2.15
	BC2	0.66	0.47			−28.79			BC2	0.66	0.47			−28.79
	BC3	10.41	15.23	1.23	1.67	46.30	35.77		BC3	10.41	15.23	1.48	1.67	46.30	12.84
注：p_r为超压，I_r为冲量；所有超压和冲量模拟误差的算术平均值为16.32%和5.79%。

下载: 导出CSV

表 4 Liu实验^[19]超压模拟结果与实测值对比

Table 4. Comparison of numerical predictions and Liu’s experimental^[19] results

测点	Case 1超压/MPa		误差/%	测点	Case 2超压/MPa		误差/%	测点	Case 3超压/MPa		误差/%
测点	实测	模拟	误差/%	测点	实测	模拟	误差/%	测点	实测	模拟	误差/%
so	0.21	0.12	−42.86	ry2	0.31	0.23	−25.81	ry2	0.52	0.57	9.62
ry3	0.14	0.15	7.14	ry3	0.43	0.21	−51.16	ry3	0.77	0.84	9.09
rb3	0.06	0.03	−50.00	ry4	0.39	0.23	−41.03	ry4	0.50	0.57	14.00
				rb3	0.08	0.04	−50.00	so	0.74	0.82	10.81
平均			−28.57				−42.00				10.88
注：三个工况超压模拟误差的算术平均值为−19.90%

下载: 导出CSV

表 5 反射冲量沿柱高分布特征及简化模型

Table 5. Distribution characteristics and simplified models of reflected impulse along column length

爆心高度（H_c/W^1/3）	分布特征及简化模型
爆心高度（H_c/W^1/3）	小比例爆距（0.4≤Z/(m∙kg^−1/3)≤0.6）	大比例爆距（0.6<Z/(m∙kg^−1/3)≤1.4）
小比例爆高（≤0.3 m/kg^1/3）	“1/x函数形分布”，可简化为三折线模型	柱底“针状峰值”分布，可简化为双折线模型
大比例爆高（>0.3 m/kg^1/3）	“局部峰值”归为局部突变分布，建议通过数值方法确定荷载	线性分布，可简化为双折线模型

下载: 导出CSV

表 6 反射冲量拟合计算公式的参数取值

Table 6. Coefficients of regressed formula for reflected impulse calculation

i	a	b	c	d
I_r-b	1.9833	−1.4980	0.0062	−2.2080
I_r-m/I_r-b	0.4690	0.9709	−0.0068	−1.3546
I_r-t/I_r-b	0.3440	0.7766	−0.0491	−0.5805
I_r-a/I_r-b	0.8792	0.4335	−0.0126	−1.4213

下载: 导出CSV

表 7 不同高度处最大反射冲量横向分布（$Z=0.4\;{\mathrm{m/kg}}^{1/3},\;H_{\mathrm{c}}/W^{1/3}=0.1\;{\mathrm{m/kg}}^{1/3} $）

Table 7. Transverse distribution of the maximum impulse at different column heights ($Z=0.4\;{\mathrm{m/kg}}^{1/3},\;H_{\mathrm{c}}/W^{1/3}=0.1\;{\mathrm{m/kg}}^{1/3} $)

h/mm	I_mid/(kPa∙s)	I_edg/(kPa∙s)	I_edg/I_mid
200	6.54	6.06	0.93
700	2.25	2.11	0.94
1200	0.63	0.62	0.98
1700	0.35	0.35	0.99
2200	0.12	0.11	0.94
注：I_mid为中轴线最大反射冲量，I_edg为柱边缘最大反射冲量。

下载: 导出CSV

表 8 双折线与三折线模型冲量误差比较

Table 8. Comparison of fitting results and numerical predictions

h/m	冲量/(kPa·s)				误差/%
	双折线		三折线		误差/%
	模拟值	公式计算值	模拟值	公式计算值	双折线	三折线
0	2.52	2.41	7.09	6.37	−4.37	−10.16
0.5	1.90	1.91	3.16	2.27	0.52	−28.17
1.0	1.25	1.43	1.60	1.46	14.40	−8.75
1.5	0.94	0.94	0.88	0.66	0.00	−25.00
2.0	0.68	0.8	0.54	0.49	17.65	−9.26
平均误差					5.64	−16.27

下载: 导出CSV

参考文献(20)

[1]	YU J, RINDER T, STOLZ A, et al. Dynamic progressive collapse of an RC assemblage induced by contact detonation [J]. Journal of Structural Engineering, 2014, 140(6): 04014014. DOI: 10.1061/(ASCE)ST.1943-541X.0000959.
[2]	SHI Y C, LI Z X, HAO H. A new method for progressive collapse analysis of RC frames under blast loading [J]. Engineering Structures, 2010, 32(6): 1691–1703. DOI: 10.1016/j.engstruct.2010.02.017.
[3]	LI J, HAO H. Numerical study of structural progressive collapse using substructure technique [J]. Engineering Structures, 2013, 52: 101–113. DOI: 10.1016/j.engstruct.2013.02.016.
[4]	BRODE H L. Blast wave from a spherical charge [J]. The Physics of Fluids, 1959, 2(2): 217–229. DOI: 10.1063/1.1705911.
[5]	BAKER W E. Explosions in air [M]. Austin, Texas: University of Texas Press, 1973: 6–10.
[6]	HENRYCH J. The dynamics of explosion and its use [M]. Amsterdam: Elsevier Scientific, 1979: 161–164.
[7]	Department of Defense (DOD). UFC 3-340-02 Structures to resist the effects of accidental explosions [S]. Washington: Department of Defense, 2008.
[8]	汪维, 刘光昆, 赵强, 等. 近爆作用下方形板表面爆炸载荷分布函数研究 [J]. 中国科学: 物理学力学天文学, 2020, 50(2): 024615. DOI: 10.1360/SSPMA-2019-0188. WANG W, LIU G K, ZHAO Q, et al. Study on load distributing function of square slab surface under close-in blast loading [J]. Scientia Sinica: Physica, Mechanica & Astronomica, 2020, 50(2): 024615. DOI: 10.1360/SSPMA-2019-0188.
[9]	WU C Q, HAO H. Modeling of simultaneous ground shock and airblast pressure on nearby structures from surface explosions [J]. International Journal of Impact Engineering, 2005, 31(6): 699–717. DOI: 10.1016/j.ijimpeng.2004.03.002.
[10]	XIAO W, ANDRAE M, GEBBEKEN N. Development of a new empirical formula for prediction of triple point path [J]. Shock Waves, 2020, 30(6): 677–686. DOI: 10.1007/s00193-020-00968-7.
[11]	CUI J, SHI Y C, LI Z X, et al. Failure analysis and damage assessment of RC columns under close-in explosions [J]. Journal of Performance of Constructed Facilities, 2015, 29(5): B4015003. DOI: 10.1061/(ASCE)CF.1943-5509.0000766.
[12]	闫秋实, 杜修力. 典型地铁车站柱在爆炸荷载作用下损伤评估方法研究 [J]. 振动与冲击, 2017, 36(1): 1–7. DOI: 10.13465/j.cnki.jvs.2017.01.001. YAN Q S, DU X L. Damage evaluation for a column of a typical subway station subjectedto internal blast loading [J]. Journal of Vibration and Shock, 2017, 36(1): 1–7. DOI: 10.13465/j.cnki.jvs.2017.01.001.
[13]	CHEN L, HU Y, REN H Q, et al. Performances of the RC column under close-in explosion induced by the double-end-initiation explosive cylinder [J]. International Journal of Impact Engineering, 2019, 132: 103326. DOI: 10.1016/j.ijimpeng.2019.103326.
[14]	DUA A, BRAIMAH A, KUMAR M. Experimental and numerical investigation of rectangular reinforced concrete columns under contact explosion effects [J]. Engineering Structures, 2020, 205: 109891. DOI: 10.1016/j.engstruct.2019.109891.
[15]	YU J, YU X F, TANG J H, et al. Local damage of precast concrete columns with grout sleeve connections under contact detonation [J]. Engineering Structures, 2022, 265: 114499. DOI: 10.1016/j.engstruct.2022.114499.
[16]	HU Y, CHEN L, FANG Q, et al. Blast loading model of the RC column under close-in explosion induced by the double-end-initiation explosive cylinder [J]. Engineering Structures, 2018, 175: 304–321. DOI: 10.1016/j.engstruct.2018.08.013.
[17]	彭玉林, 吴昊, 方秦. 爆炸荷载在圆截面桥梁墩柱上的分布规律 [J]. 爆炸与冲击, 2019, 39(12): 122201. DOI: 10.11883/bzycj-2018-0317. PENG Y L, WU H, FANG Q. Blast loading distributions on the circular sectional bridge columns [J]. Explosion and Shock Waves, 2019, 39(12): 122201. DOI: 10.11883/bzycj-2018-0317.
[18]	WOODSON S C, BAYLOT J T. Structural collapse: quarter-scale model experiments [R]: Vicksburg: Engineer Research and Development Center, US Army Corps of Engineers, 1999.
[19]	LIU L, MA Z J, ZONG Z H, et al. Blast response and damage mechanism of prefabricated segmental RC bridge piers [J]. Journal of Bridge Engineering, 2021, 26(4): 04021012. DOI: 10.1061/(ASCE)BE.1943-5592.0001698.
[20]	胡志乐, 马亮亮, 吴昊, 等. 远距离近地面爆炸空气冲击波计算的网格尺寸优化与验证 [J]. 爆炸与冲击, 2022, 42(11): 114201. DOI: 10.11883/bzycj-2021-0499. HU Z L, MA L L, WU H, et al. Optimization and verification of mesh size for air shock wave from large distance and near ground explosion [J]. Explosion and Shock Waves, 2022, 42(11): 114201. DOI: 10.11883/bzycj-2021-0499.