内部爆炸载荷作用下砌体墙结构的失效规律

曹宇航; 张晓伟; 张庆明

doi:10.11858/gywlxb.20210810

内部爆炸载荷作用下砌体墙结构的失效规律

doi: 10.11858/gywlxb.20210810

北京理工大学爆炸科学与技术国家重点实验室，北京　100081

基金项目: 国家重点研发计划项目子课题（2016YFC0801204）

详细信息

作者简介:
曹宇航（1996－），男，硕士研究生，主要从事毁伤防护与技术研究. E-mail：yanxiu153@qq.com

通讯作者:
张晓伟（1982－），男，博士，副教授，主要从事冲击动力学、毁伤与防护技术研究.E-mail：mezhangxw@bit.edu.cn

中图分类号: O383
计量
- 文章访问数: 1141
- HTML全文浏览量: 611
- PDF下载量: 40
出版历程
- 收稿日期: 2021-06-10
- 修回日期: 2021-07-07

Failure Characteristics of Masonry Wall under Internal Explosion

State Key Laboratory of Explosive Science and Technology, Beijing Institute of Technology, Beijing 100081, China

摘要

摘要: 为了探究钢筋混凝土框架-砌体墙建筑物的内爆毁伤效应，开展了内爆载荷下砌体墙结构动力响应和失效规律研究。采用数值模拟方法，并结合理论分析，研究了不同当量内爆载荷下砌体墙的失效机理以及失效形式转化过程。结果表明，小药量内部爆炸时，砌体墙的主要失效形式为弯曲导致的墙面开裂。随着药量增加，墙体边界在首道冲击波反射超压作用下发生以剪切变形为主导的失效。而在大药量条件下，首道冲击波超压使砌体墙材料达到极限抗压强度而发生压溃失效。研究结果可为砌体墙结构毁伤评估与防护设计提供技术参考。
- 钢筋混凝土 /
- 砌体墙 /
- 内部爆炸载荷 /
- 失效模式 /
- 数值模拟
Abstract: In order to investigate the internal explosion damage effect of reinforced concrete frame-masonry wall, the dynamic response and failure characteristics of masonry walls under internal explosion were studied. Based on numerical simulation and theoretical analysis, the failure mechanism and failure transformation process of masonry wall under internal explosion with different TNT equivalence were analyzed. The results showed that under internal explosive with small mass of explosive, the failure of the masonry wall is mainly due to bending effect. With the increase of the mass of explosive, the failure dominated by shear deformation will happen due to the reflection of the first shock wave. For the cases with larger mass of explosive, compressive failure occurs because the overpressure of the first shock wave reaches the ultimate compressive strength of masonry wall. The research results can provide technical reference for explosive damage evaluation and protection design of masonry wall structures.
- reinforced concrete /
- masonry wall /
- internal explosion /
- failure mode /
- numerical simulation

HTML全文

烧结钕铁硼永磁材料具有优越的磁性能，广泛应用于诸多领域，如航空航天、机械加工、电子计算机等。其中基于Nd₂Fe₁₄B合金的爆炸驱动脉冲电源被认为是最有效的紧凑型脉冲电源^[1]。对于烧结钕铁硼永磁材料，人们关注的重点是其磁性能，对其力学性能的研究则相对较少，尤其是动态力学性能。然而，在现代精密仪器中经常用到烧结钕铁硼，如汽车、导弹制导、国防通讯设备，这些仪器在使用过程中常常承受外力作用，要求其具有一定的抗冲击性能，因此研究烧结钕铁硼在冲击压缩下的动态行为和断裂机理具有很高的学术意义和应用价值。李岩峰等^[2]在一级轻气炮系统中对Nd-Fe-B磁体进行了冲击压缩实验，研究了退磁现象；李巧燕等^[3]采用二级轻气炮，测量了钕铁硼永磁体在19～78 GPa压力范围内的冲击压缩特性。

本课题组在前期工作中利用MTS万能试验机和分离式霍普金森压杆（Split Hopkinson Pressure Bar，SHPB）研究了烧结钕铁硼在准静态下的抗弯强度和弹性模量^[4]以及在单轴压缩下的动态断裂^[5]；此外还开展了动态巴西圆盘试验^[6]，研究了烧结钕铁硼的断裂过程，根据弹性理论得到了其拉伸强度，并利用扫描电子显微镜（Scanning Electron Microscope，SEM）观察其微观结构。本研究将利用一级气体炮开展烧结钕铁硼在高应变率、一维平面冲击加载下的层裂实验，利用激光干涉测速技术VISAR（Velocity Interferometer System for Any Reflector）测量自由面粒子速度历史，确定其层裂行为，最后利用SEM观察断口微观形貌，分析烧结钕铁硼在冲击波作用下的断裂机理。

1. 实验方案和原理

采用宁波大学的一级气体炮进行冲击加载实验，图1为实验装置示意图。采用铝制弹托驱动飞片，飞片采用有机玻璃（PMMA）和铝合金（LY12）两种材料，直径均为42.5 mm。通过改变弹丸的速度和飞片的材质，实现较大范围的加载应力。利用测速探针测量飞片的撞击速度，利用VISAR测量样品自由面粒子速度历史。

图 1 实验装置示意图

Figure 1. Schematic of experimental setup

下载: 全尺寸图片幻灯片

层裂是材料的一种动态失效形式，是由稀疏（膨胀）波相互作用产生拉应力而形成的^[7]。层裂原理如图2所示。飞片以速度V撞击样品，在飞片和样品内产生冲击压缩波；压缩波到达飞片和样品自由面时反射稀疏波，随后稀疏波在样品内相遇，相互作用产生拉伸应力；当拉伸应力超过材料的强度极限时，样品发生破坏，即发生层裂；层裂形成后，从新产生的自由面上反射稀疏波，使拉伸应力降低；当稀疏波到达样品自由面时，下降的自由面粒子速度曲线立即回跳，这便是发生层裂行为的信号。根据带有层裂信号的自由面速度曲线，通过（1）式确定样品的层裂强度^[8]，即

图 2 层裂原理示意图

Figure 2. Schematic illustration of spall process

下载: 全尺寸图片幻灯片

${\sigma _{{\rm{sp}}}} = \frac{1}{2}{\rho _0}{c_{\rm{L}}}({u_{{\rm{f}},\max }} - {u_{{\rm{f}},\min }})$

(1)

式中：u_{f, max}和u_{f, min}是出现第1次回跳信号前自由面粒子速度的最大值和最小值， ${\sigma _{{\rm{sp}}}}$ 为层裂强度， ${\rho _0}$ 为样品的初始密度，c_L为样品的弹性纵波声速。

实验材料采用牌号为N42的烧结钕铁硼块状磁体，由宁波市三环钕铁硼磁业有限公司提供，采用未充磁的粉末冶金方法制备。通过电火花线切割方式，加工成直径为40 mm、厚度为4 mm的圆形样品，样品的厚度方向即磁化方向，也是冲击压缩方向，最后利用砂纸对样品表面进行打磨抛光。烧结钕铁硼的实测平均密度 ${\rho _0}$ 为7.544 g/cm³，泊松比 $\nu$ 取0.24，弹性模量E为155 GPa，根据

${c_{\rm L}} = {\left[ {\frac{{1 - \nu }}{{(1 + \nu )(1 - 2\nu )}}\frac{E}{{{\rho _0}}}} \right]^{1/2}}$

(2)

可得一维应变下钕铁硼样品的弹性纵波声速c_L=4.921 km/s^[7]。表1列出了两种飞片材料参数，其中 ${c_{0{\rm{i}}}}$ 、 ${\lambda _{\rm{i}}}$ 为Hugoniot关系常数， ${\rho _{0{\rm{i}}}}$ 为飞片的初始密度。

表 1 飞片材料的Hugoniot参数^[9]

Table 1. Hugoniot parameters of flyers^[9]

Material	${\rho _{0{\rm{i}}}}$ $/({\rm{g}} \cdot {\rm{c}}{{\rm{m}}^{ - 3}})$	${c_{0{\rm{i}}}}$ $/({\rm{km}} \cdot {{\rm{s}}^{ - 1}})$	${\lambda _{\rm{i}}}$
PMMA	1.186	2.65	1.54
LY12	2.784	5.37	1.29

下载: 导出CSV

| 显示表格

2. 实验结果与分析

烧结钕铁硼样品的层裂实验结果如表2所示，其中加载应力 $\sigma$ 根据飞片和样品碰撞时的动量守恒及冲击界面上连续性条件计算得到。计算时，飞片采用流体力学模型，烧结钕铁硼样品采用线弹性模型^[10]。

表 2 烧结钕铁硼层裂实验参数和结果

Table 2. Parameters and results of spall experiments for Nd-Fe-B

Exp. No.	Material of flyer	Thickness of flyer/mm	V/(m·s^–1)	Thickness of sample/mm	$\sigma$ /GPa	${\sigma _{{\rm{sp}}}}$ /GPa
01	PMMA	1.17	122.67	3.99	0.375	0.209
02	PMMA	1.16	173.99	3.99	0.550	0.249
03	PMMA	1.17	230.72	3.96	0.739	0.263
04	LY12	2.28	146.33	3.95	1.591	0.274
05	LY12	2.33	161.90	3.90	1.737	0.313
06	LY12	2.28	190.89	3.94	2.092	0.251
07	LY12	2.34	230.73	3.95	2.512	0.224

下载: 导出CSV

| 显示表格

${\sigma _{\rm{i}}} = {\rho _{0{\rm{i}}}}({c_{0{\rm{i}}}} + {\lambda _{\rm{i}}}{u_{\rm{i}}}){u_{\rm{i}}}$

(3)

${\sigma _{\rm{s}}} = {\rho _{0{\rm{s}}}}{c_{\rm{L}}}{u_{\rm{s}}}$

(4)

${\sigma _{\rm{s}}} = {\sigma _{\rm{i}}} = \sigma ,\quad {u_{\rm{s}}} + {u_{\rm{i}}} = V$

(5)

式中：u为粒子速度，V为飞片撞击速度， ${\rho _0}$ 为初始密度，下标i和s分别表示飞片和样品。

由自由面速度u_f可以得到粒子速度，按忽略熵增的方法处理，则有^[11]

${u_{\rm{f}}} = u + {u_{\rm{r}}} \approx 2u$

(6)

式中： ${u_{\rm{r}}}$ 表示由中心稀疏波引起的附加质点速度。

层裂实验实测样品的自由面速度剖面如图3所示。可以看出，烧结钕铁硼在飞片撞击下发生了层裂破坏，所有实测曲线均有明显的层裂信号，如图3中黑色箭头所示。在冲击载荷下烧结钕铁硼的层裂强度与加载应力的关系如图4所示。可见：层裂强度随着加载应力的增加先增大后减小；存在一个明显的临界应力，约为1.7 GPa，当加载应力超过该临界应力时，材料的层裂强度就会减小。该临界应力值同时反映了冲击载荷下材料是否发生压缩损伤，如图3中06号实验所测自由面速度曲线所示，在速度上升过程中，接近峰值时速度突然变缓，说明材料在压缩过程中发生压缩损伤，材料开始由弹性向非弹性转化，导致材料的层裂强度减小。

图 4 层裂强度与加载应力的关系

Figure 4. The relationship between spall strength and impact stress

下载: 全尺寸图片幻灯片

图 3 层裂实验样品的自由面速度剖面

Figure 3. Free surface velocity profile of the sample in the spall experiment

下载: 全尺寸图片幻灯片

通过SEM对冲击压缩实验后烧结钕铁硼样品的断口进行观察，如图5所示，可见晶粒表面非常光滑，断口比较平整，表现出脆性断裂的明显特征。另外，从图5（b）中可以看出，微裂纹处有明显的穿晶断裂。

图 5 断口的SEM图像

Figure 5. SEM images of the fractured surface

下载: 全尺寸图片幻灯片

对比之前工作中烧结钕铁硼的断口形貌发现，在准静态下的三点弯曲实验中烧结钕铁硼主要是沿晶断裂^[4]，在动态巴西圆盘试验中出现穿晶断裂现象^[6]，而本实验中在微裂纹处出现多处穿晶断裂。准静态压缩下断口表面极不平整，充满了尖锐的边缘，这是由于晶界的强度较弱，裂纹沿着晶粒边界传播，说明在低加载速率下沿晶断裂在材料失效过程中占主导地位。随着加载速率的升高，穿晶断裂越来越明显，说明当加载速率小于某个临界值时，微裂纹沿晶界传播，导致断口不均匀，即沿晶断裂；当加载速率大于某个临界值时，由于裂纹尖端附近释放的应变能显著增大，导致其前面的单个晶粒分裂，这时断口非常平整、光滑，即穿晶断裂。因此，冲击波加载下烧结钕铁硼的失效过程是由两者共同控制下的沿晶断裂和穿晶断裂。

3. 结　论

在一级气体炮平面冲击压缩实验中，采用VISAR测量了烧结钕铁硼样品在冲击载荷下的自由面速度历史，得到了烧结钕铁硼在0.375～2.512 GPa范围内的层裂强度；分析了层裂强度与加载应力的关系，发现存在一个临界应力阈值（约1.7 GPa），当加载应力大于该阈值时，材料发生压缩损伤，导致烧结钕铁硼的层裂强度减小。观察烧结钕铁硼断口微观形貌发现，冲击波加载下烧结钕铁硼有明显的穿晶断裂，失效过程由沿晶断裂和穿晶断裂共同控制，而在准静态下则由沿晶断裂控制。研究结果有助于进一步了解烧结钕铁硼在冲击加载下的损伤和破坏。

图 1 单个房间内部爆炸有限元模型

Figure 1. Finite element model of a single room under internal explosion

下载: 全尺寸图片幻灯片

图 2 砌体墙的材料模型

Figure 2. Material models for the masonry wall

下载: 全尺寸图片幻灯片

图 3 文献[15]中的实验装置示意图

Figure 3. Schematic diagram of experiment device in Ref.[15]

下载: 全尺寸图片幻灯片

图 4 实验砌体墙及有限元模型

Figure 4. Masonry walls for experiment and finite element model

下载: 全尺寸图片幻灯片

图 5 墙体破坏范围对比

Figure 5. Comparison of masonry wall failure range

下载: 全尺寸图片幻灯片

图 6 房间中部剖面压力云图

Figure 6. Pressure distribution at the middle cross-section of room

下载: 全尺寸图片幻灯片

图 7 墙面壁面压力峰值分布

Figure 7. Peak pressure distribution of t wall surface

下载: 全尺寸图片幻灯片

图 8 墙面中心壁面压力时程曲线

Figure 8. Pressure history curves at the center of the wall

下载: 全尺寸图片幻灯片

图 9 典型砌体墙开裂形式^[7]

Figure 9. Typical cracking modes of masonry wall^[7]

下载: 全尺寸图片幻灯片

图 10 竖向裂缝示意图

Figure 10. Schematic diagram of vertical cracking

下载: 全尺寸图片幻灯片

图 11 水平裂缝示意图

Figure 11. Schematic diagram of horizontal cracking

下载: 全尺寸图片幻灯片

图 12 斜向裂缝示意图

Figure 12. Schematic diagram of diagonal cracking

下载: 全尺寸图片幻灯片

图 13 双向固支砌体墙裂缝分布

Figure 13. Crack distribution of masonry wall with two-way fixed support

下载: 全尺寸图片幻灯片

图 14 横向载荷p与斜向裂缝角度θ的关系曲线

Figure 14. Relationship between lateral load p and diagonal cracking angle $\theta$

下载: 全尺寸图片幻灯片

图 15 6 kg TNT药量下砌体墙的失效过程

Figure 15. Failure process of masonry wall with 6 kg TNT explosion

下载: 全尺寸图片幻灯片

图 16 弯曲主导的失效模式中不同药量下墙体的失效特征（上图为开裂形式，下图为速度分布）

Figure 16. Failure characteristics of the wall dominated by bending (Upper figures correspond to cracking form, and lower figures correspond to velocity distribution.)

下载: 全尺寸图片幻灯片

图 17 50 kg TNT药量下墙体的毁伤过程（1/4模型）

Figure 17. Damage process of wall with 50 kg TNT (1/4 model)

下载: 全尺寸图片幻灯片

图 18 不同药量下墙体中心水平轴线速度分布（40 ms）

Figure 18. Velocity distribution at the center horizontal axis of the wall with different TNT equivalence (40 ms)

下载: 全尺寸图片幻灯片

图 19 100 kg TNT药量下中心水平轴线不同位置的超压曲线

Figure 19. Overpressure curves at different locations of the center horizontal axis of the wall with 100 kg TNT

下载: 全尺寸图片幻灯片

图 20 100 kg TNT药量下墙体的毁伤过程（1/4模型）

Figure 20. Damage process of the wall with 100 kg TNT (1/4 model)

下载: 全尺寸图片幻灯片

图 21 100 kg TNT药量下墙体中心水平轴线速度分布（40 ms）

Figure 21. Velocity distribution at the center horizontal axis of the wall with 100 kg TNT equivalence (40 ms)

下载: 全尺寸图片幻灯片

表 1 砖块和砂浆的基本材料参数^[10-12]

Table 1. Basic material parameters of brick and mortar^[10-12]

Material	p- $\alpha$ equation of state					RHT strength model
Material	$\,\rho$ ₀/(kg·m⁻³)	$\alpha_0$	N	p_el/MPa	p_comp/GPa	E/GPa	ν	f_c/MPa	f_t/ f_c	f_v/ f_c
Mortar	2100	1.33	3	2.4	2.5	1.44	0.20	7.2	0.08	0.15
Brick	1800	1.33	3	5.2	2.5	4.71	0.12	15.7	0.10	0.18

下载: 导出CSV

表 2 空气材料参数^[13]

Table 2. Material parameters of air^[13]

$\,\rho$ /(kg·m⁻³)	T/K	c/(kJ·kg⁻¹·K⁻¹)
1.225	288.2	0.718

下载: 导出CSV

表 3 TNT材料参数^[13]

Table 3. Material parameters of TNT^[13]

$\,\rho$ /(kg·m⁻³)	A/GPa	B/GPa	R₁	R₂	$\omega$	E_m0/(kJ·m⁻³)
1630	373.8	3.75	4.15	0.9	0.35	6.0 × 10⁶

下载: 导出CSV

表 4 数值模拟结果与实验数据对比

Table 4. Comparison of numerical simulation results and experimental data

Method	Peak overpressure of shock wave/MPa
Method	2 m	6 m	10 m	16 m	20 m	28 m	36 m
Experiment^[16]	7.684	1.624	1.057	0.721	0.605	0.467	0.387
Numerical simulation	6.755	1.430	0.969	0.751	0.613	0.457	0.368
Relative error/%	−12.10	−11.90	−8.33	4.16	1.32	−2.14	−4.91

下载: 导出CSV

参考文献(18)

[1]	SMITH P D, MAYS G C, ROSE T A, et al. Small scale models of complex geometry for blast overpressure assessment [J]. International Journal of Impact Engineering, 1992, 12(3): 345–360. doi: 10.1016/0734-743X(92)90112-7
[2]	王瑀. 典型结构内部爆炸波的传播规律与超压荷载模型[D]. 天津: 天津大学, 2008. WANG Y. Interior propagation of blast wave and overpressure model for common structures [D]. Tianjin: Tianjin University, 2008.
[3]	张晓伟, 张庆明, 施鹏, 等. 内爆条件下建筑物等效载荷研究 [J]. 北京理工大学学报, 2013, 33(2): 133–138. doi: 10.3969/j.issn.1001-0645.2013.02.005 ZHANG X W, ZHANG Q M, SHI P, et al. Investigation on the equivalent loading of constructions under internal explosion [J]. Journal of Beijing Institute of Technology, 2013, 33(2): 133–138. doi: 10.3969/j.issn.1001-0645.2013.02.005
[4]	陈鹏宇, 侯海量, 金键, 等. 舰船舱内爆炸载荷简化载荷计算模型 [J]. 舰船科学技术, 2020, 42(17): 22–29. CHEN P Y, HOU H L, JIN J, et al. Simplified calculation model for explosion loading in ship cabin [J]. Ship Science and Technology, 2020, 42(17): 22–29.
[5]	杨科之, 董军, 吴平安, 等. 砌体隔墙的抗爆炸能力分析[C]//第十三届全国结构工程学术会议论文集（第Ⅲ册）. 北京: 清华大学出版社, 2004: 146–151.
[6]	SCHNEIDER J M, RAMIN M V, STOTTMEISTER A, et al. Characterization of debris throw from masonry wall sections subjected to blast [J]. Engineering Structures, 2019, 203: 109729.
[7]	KEYS R A, CLUBLEY S K. Experimental analysis of debris distribution of masonry panels subjected to long duration blast loading [J]. Engineering Structures, 2017, 130: 229–241.
[8]	ZHOU X, HAO H, DEEKS A. Numerical modeling of response and damage of masonry walls to blast loading [J]. Transactions of Tianjin University, 2006, 12: 132–137.
[9]	李娟, 丛培宇, 刘香, 等. 内爆作用下不同类型填充墙对RC框架结构的影响 [J]. 工程爆破, 2020, 26(4): 28–34. doi: 10.3969/j.issn.1006-7051.2020.04.005 LI J, CONG P Y, LIU X, et al. Influence of different types infill walls on RC frame structure under implosion [J]. Engineering Blasting, 2020, 26(4): 28–34. doi: 10.3969/j.issn.1006-7051.2020.04.005
[10]	高连玉, 徐建. 砌体结构设计规范: GB 50003—2011 [S]. 北京: 中国建筑工业出版社, 2011.
[11]	MICHALOUDIS G, GEBBEKEN N. Modeling masonry walls under far-field and contact detonations [J]. International Journal of Impact Engineering, 2019, 123: 84–97.
[12]	范俊余, 方秦, 张亚栋, 等. 砖墙抗爆性能的数值模拟研究 [J]. 防护工程, 2011(5): 35–40. FAN J Y, FANG Q, ZHANG Y D, et al. Numerical simulation of the anti-blast properties of masonry walls [J]. Protective Engineering, 2011(5): 35–40.
[13]	Century Dynamics Inc.. Theory manual for Autodyn [M]. Texas: Century Dynamics Inc., 2005.
[14]	汪维, 刘瑞朝, 李林, 等. 不同爆炸距离下单向支撑方形钢筋混凝土板破坏模式数值模拟研究 [J]. 兵工学报, 2015, 36(Suppl 1): 233–241. WANG W, LIU R C, LI L, et al. Numerical simulation of one-way square reinforced concrete slab at different blast distances [J]. Acta Armamentarii, 2015, 36(Suppl 1): 233–241.
[15]	张玉磊, 王胜强, 袁建飞, 等. 方形坑道内爆炸冲击波传播规律 [J]. 含能材料, 2020, 28(1): 46–51. ZHANG Y L, WANG S Q, YUAN J F, et al. Experimental study on the propagation law of blast waves in a square tunnel [J]. Chinese Journal of Energetic Materials, 2020, 28(1): 46–51.
[16]	熊维. 近距离爆炸下砌体填充墙局部破坏试验研究与数值模拟 [D]. 天津: 天津大学, 2015. XIONG W. Experimental study and numerical simulation of local damage of unreinforced masonry under close-in explosion [D]. Tianjin: Tianjin University, 2015.
[17]	Department of Defense. UFC 3-340-02 structures to resist the effects of accidental explosions [M]. New York: National Institute of Building Sciences, 2008.
[18]	JONES N. 结构冲击 [M]. 2版. 许骏, 蒋平, 译. 北京: 国防工业出版社, 2018: 57–72. JONES N. Structural impact [M]. 2nd ed.Translated by XU J, JIANG P. Beijing: National Defense Industry Press, 2018: 57–72.