• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST、EBSCO、DOAJ收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

半空间内含有部分脱胶的椭圆夹杂及圆孔对SH波的散射

齐辉 陈洪英 张希萌 赵元博 项梦

樊保龙, 白春华, 王博, 高康华, 李斌. 大尺寸密闭容器内天然气的爆炸超压场[J]. 爆炸与冲击, 2018, 38(2): 404-408. doi: 10.11883/bzycj-2016-0191
引用本文: 齐辉, 陈洪英, 张希萌, 赵元博, 项梦. 半空间内含有部分脱胶的椭圆夹杂及圆孔对SH波的散射[J]. 爆炸与冲击, 2018, 38(6): 1344-1352. doi: 10.11883/bzycj-2017-0142
FAN Baolong, BAI Chunhua, WANG Bo, GAO Kanghua, LI Bin. Explosion overpressure field of natural gas in a large-scaled confined vessel[J]. Explosion And Shock Waves, 2018, 38(2): 404-408. doi: 10.11883/bzycj-2016-0191
Citation: QI Hui, CHEN Hongying, ZHANG Ximeng, ZHAO Yuanbo, XIANG Meng. Scattering of SH-wave by elliptical inclusion with partial debond curve and circular cavity in half space[J]. Explosion And Shock Waves, 2018, 38(6): 1344-1352. doi: 10.11883/bzycj-2017-0142

半空间内含有部分脱胶的椭圆夹杂及圆孔对SH波的散射

doi: 10.11883/bzycj-2017-0142
基金项目: 

黑龙江省自然科学基金项目 A201404

详细信息
    作者简介:

    齐辉(1963-), 男, 博士, 教授, 博导

    通讯作者:

    陈洪英, hongyingchen@hrbeu.edu.cn

  • 中图分类号: O343.4

Scattering of SH-wave by elliptical inclusion with partial debond curve and circular cavity in half space

  • 摘要: 采用复变函数法,结合"保角映射"技术及Green函数法,研究SH波作用下半空间内含有部分脱胶的椭圆夹杂以及圆形孔洞的散射问题。首先,利用"保角映射"技术将椭圆夹杂映射为圆夹杂,求出散射波位移场,同时,利用Green函数法与"虚设点源"的方法,求出半空间内椭圆夹杂以及圆孔的位移及应力场;然后,根据椭圆夹杂周围位移、应力连续、圆孔周围应力自由的边界条件,建立无穷线性代数方程组,求解出波函数中的未知系数;最后,在脱胶部分施加大小相等、方向相反的应力,构造出"脱胶模型",得到半空间内含有部分脱胶的椭圆夹杂以及圆形孔洞的总位移场。数值算例表明,入射角度、入射波频率、缺陷之间的距离、夹杂埋深及脱胶角度等对动应力集中因子有较大影响。
  • 当可燃气体、可燃液体的蒸气(或可燃粉尘)与空气混合并达到一定浓度时,遇到火源就会发生爆炸。这个能够发生爆炸的浓度范围,叫做爆炸极限,通常用可燃气体、蒸气或粉尘在空气中的体积分数来表示。可燃气体在空气(氧气)中的爆炸极限范围是众多学者关心的焦点,但关于作为描述可燃气体爆炸后毁伤效果的典型参数,如爆压、爆温、爆速等的报道较少[1-4]。这是因为目前对爆炸极限的研究多局限于两个方面:一方面是因为偏于实际生产应用的缘故,多数测试中只需知道可燃气体在空气中的爆炸极限范围、临界氧含量或者获得爆炸三角形图即可,而对具体的燃烧爆炸过程及其结果并不关心;另一方面是由于测试仪器的局限性,多采用小尺寸容器,如小型激波管、20 L球、圆柱形爆炸罐等,其测试手段较单一,多数只在耐压容器壁面安置一个压力传感器或温度传感器,所得数据有限[5-7]。正因为如此,目前对于气体爆炸的研究多采用放宽条件,改变初始温度、初始压力、当量体积分数等手段,以获得某种可燃气体或气体混合物较全面的爆炸特性参数[8-10]。大尺寸密闭空间内可燃气体爆炸过程更贴合于实际。但由于爆炸容器尺寸大,操作复杂,开展相关研究难度较大。本文中,通过在大尺寸密闭容器中开展天然气爆炸超压场的研究,以期获得大尺寸密闭空间内天然气爆炸超压的发展规律,丰富目前天然气燃爆威力的测试数据,为密闭空间内天然气爆炸危害的预防及毁伤能力评估提供数据支持。

    天然气-空气混合物的燃爆过程在容积为10 m3的爆炸罐内进行,爆炸罐示意图见图 1。为更清楚地记录爆炸罐内天然气-空气混合物爆炸后超压的发展过程,在爆炸罐内沿罐体轴向典型位置布置4个压力传感器,沿罐体径向典型位置布置3个压力传感器,以获得天然气-空气混合物燃爆后,爆炸罐内部的超压场状态,传感器布置见图 2。点火装置选用高能放电器,单次点火,点火能量为40 J;超压测试系统包括压力传感器、信号调理器、数据采集系统、信号线等。压力传感器为PCB公司的113B系列;信号调理器为PCB公司的信号调理器;数据采集系统为南汇科技虚拟仪器系统;高速摄像系统为Photron公司的NX100相机,实验中采样频率为1 000 s-1

    图  1  10 m3爆炸罐结构示意图
    Figure  1.  Schematic structure of 10 m3 explosion vessel
    图  2  罐体内压力传感器布置示意图
    Figure  2.  Layout of pressure transducers in the vessel

    选用华北地区的工业天然气为实验样品,其密度为728.9 g/m3, 高位发热量为40.38 MJ/m3。该工业天然气的组分及其体积分数分别为:氧气,0.05%;氮气,1.28%;天然气,92.40%;二氧化碳,1.72%;乙烷,3.62%;丙烷,0.65%;异丁烷,0.11%;正丁烷,0.11%;异戊烷,0.04%;正戊烷,0.02%。

    将不同量的天然气在容积为10 m3的爆炸罐内与空气进行均匀混合,测量各点天然气的体积分数,达到罐体内各位置处天然气体积分数相对均匀时,进行点火操作,混合过程见文献[11],利用压力测试系统记录相关测试点的压力数据。

    超压状态场通常被用来评估受限空间内可燃气体点爆过程中的爆炸效果。对大尺寸密闭空间来说,其超压状态场与可燃气体的体积分数存在对应关系[11]。本次研究主要从近爆炸下限(5.4%)的天然气点爆过程入手,分析不同体积分数下的天然气爆炸超压状态场。

    经过系列实验测试后发现,在实际天然气体积分数接近天然气爆炸下限(5.4%)处,开展天然气点爆实验时,能够获得3种超压曲线,且曲线状态随天然气实际体积分数的不断升高发生一系列的变化。在天然气爆炸下限附近选取5.5%、5.8%和6.5%等3种天然气体积分数进行实验,以轴向第一个传感器的信号为典型信号,来研究不同体积分数条件下压力传感器获取的不同超压曲线,如图 3所示。

    图  3  天然气体积分数不同时的超压时程曲线
    Figure  3.  Overpressure-time curves at different volume fractions of natural gas

    图 3可以看出,图 3(a)是天然气-空气混合物点火后典型的冲击波压力曲线,由于点火位置处天然气的体积分数较低,在40 J点火能量的作用下,天然气中可燃组分与空气发生化学反应的速率较低,从时间坐标可以看出整个反应持续了十多秒,是典型的缓慢燃烧反应。当点火位置处天然气体积分数上升至5.8%时,由图 3(b)可以看出,测得的压力-时间曲线分成两部分:蓝色椭球框内的初始冲击波压力突跃以及后续的持续燃烧过程。蓝色椭球框内的压力曲线对应着图 3(a)中压力曲线的发展状态,不同的是,当初始冲击波过去后的一个豫驰时间后(约4 s),图 3(a)中的压力曲线并无继续增长趋势,而是持续下降,而图 3(b)中的反应被进一步加速,造成了后续大范围持续燃烧的过程。

    当天然气体积分数进一步升高,达到6.5%时,超压时程曲线如图 3(c)所示。此时,从作用初期,已无法捕捉到初始冲击波的作用曲线,也无法观察到一个明显的豫驰时间,天然气爆炸后压力直接上升至最高值,前期的冲击波作用和后期的燃烧波发展已形成一个整体。产生这种现象的主要原因是:气体的爆燃过程也是一种化学反应过程,在初始环境参数不变的情况下,可燃气体体积分数越高,单位空间内的可燃气体分子越多,可燃气体分子发生有效碰撞的几率越大,反应速率越高。在点火的瞬间,当可燃气体体积分数较低时,点火源周围局部的可燃气体分子在外界能量的输入下发生反应,但由于可燃气体分子少,反应没有完全传播下去,造成了图 3(a)所示的现象,在宏观上表现为点火后产生了前导冲击波,但前导冲击波没有得到能量支持继续发展;当可燃气体体积分数较高时,化学反应速率很高,使得点火瞬间气体分子的反应从局部很快发展到整个空间,宏观上表现为前导冲击波波后气体产物运动速度追上或超过前导冲击波发展速度,使二者形成一个整体,表现为图 3(c)所示的形式。而在这二者之间,存在前导冲击波缓慢发展最终形成燃烧波的过程,如图 3(b)所示,即存在一定的豫驰时间[11]

    以容积为10 m3的爆炸罐为研究对象,对其轴向的4个压力传感器(距爆源由近至远分别命名为OP1~OP4)的压力数据进行分析,典型结果如图 4(a)所示。当天然气体积分数为5.5%时,接近实验测得的爆炸下限(5.4%),因此,此爆炸超压发展曲线图为近爆炸极限时的临界压力发展趋势图。由图 4(a)可看出,经过滤波处理后,轴向上的压力传感器随着距离点火位置的远近,其超压峰值分别为82.5、32.9、23.4和15.1 kPa。距爆源最近的传感器测得的压力曲线较接近典型的冲击波超压曲线,其他3个传感器所测得的压力曲线都接近于燃烧波的压力曲线。这主要是由于点火点处天然气的体积分数较低,接近爆炸下限,点火初期,点火位置处的天然气-空气混合物被点燃,初始冲击波产生,但由于能量支持不够,未继续发展,使得后续的轴向传感器测得的压力信号较弱且随距离呈递减趋势。

    图  4  天然气体积分数为5.5%时压力传感器测得的爆炸超压曲线
    Figure  4.  Overpressure-time curves obtained by the pressure transducers at the volume fraction 5.5% of natural gas

    仍以容积为10 m3的爆炸罐为研究对象,对其径向的3个压力传感器(距轴线由近至远分别为OP4~OP6)的压力数据进行分析,典型结果如图 4(b)所示。由图 4(b)可以看出,传感器距离爆源中心轴线越远时,爆炸超压峰值越大,但增幅不大。这是由于天然气-空气燃爆发生并沿爆炸罐体轴向传播的同时,也沿爆炸罐体径向传播,呈体积性发展趋势,距爆源一定距离后,整个燃爆体系传播过程已成整体化趋势,同一波阵面的压力数据基本相当,但由于壁面反射的影响,偏离轴线处压力可能略有升高。

    在对天然气体积分数为5.5%的天然气-空气混合物燃爆超压曲线分析后,针对3种天然气体积分数情况下的天然气-空气混合物燃爆发展进行研究,得到其燃爆参数随轴向和径向的发展规律,如图 5所示。由图 5可以看出:当天然气体积分数接近爆炸下限时,天然气-空气混合物燃爆的最高压力即为前导冲击波的超压峰值,其值相对较低;随着天然气体积分数的升高,天然气燃爆的最高压力为前导冲击波过后燃烧波的峰值压力,且此压力值随着初始天然气体积分数的升高而增大。从空间发展角度来看,距爆源距离对天然气爆炸超压峰值影响不大。而对于豫驰时间来说,与超压发展规律类似,燃爆豫驰时间与距爆源距离关系不大,初始天然气体积分数是决定性因素。

    图  5  天然气燃爆参数沿轴向和径向的发展规律
    Figure  5.  Development of natural gas deflagration parameters along the axis and radial directions

    天然气爆炸下限附近存在3种典型的超压状态:(1)当可燃气体的体积分数接近爆炸下限时,点火后只存在点火点周边气体燃烧产生的前导冲击波;(2)当可燃气体体积分数略高于爆炸下限时,点火后前导冲击波和后续燃烧波共存;(3)当可燃气体体积分数高于爆炸下限一定程度后,前导冲击波与后续燃烧波重叠。

    经过系统实验发现,爆炸下限附近的爆炸超压峰值及燃爆豫驰时间主要取决于初始天然气体积分数,而与距爆源距离关系不大,这主要是由于大尺寸密闭容器内气体燃爆过程的体积性效果。

  • 图  1  半空间含有部分脱胶的椭圆夹杂及圆形孔洞模型

    Figure  1.  Model of elliptical inclusion with a partially debond curve and circular cavity in half space

    图  2  含椭圆形弹性夹杂及圆孔的半空间作用一出平面线源荷载模型

    Figure  2.  Half-space model of elliptical inclusion and circular cavity impacted by an out-plane source load

    图  3  SH波垂直入射时τθz*随着k1a的分布

    Figure  3.  Distribution of DSCF with k1a for vertically incident SH wave

    图  4  SH波水平入射时τθz*随着θ1*的变化情况

    Figure  4.  Distribution of DSCF with θ1* for horizontally incident SH wave

    图  5  SH波水平入射时随不同角度α在不同入射波频率k1a条件下τθz*的分布

    Figure  5.  Distribution of DSCF versus the frequency of incident wave at different incident angles for horizontally incident SH wave

    图  6  SH波水平入射时τθz*随着k1*的分布

    Figure  6.  Distribution of DSCF with k1* for horizontally incident SH wave

    图  7  SH波水平入射时τθz*随着h/a的变化情

    Figure  7.  Distribution of DSCF with h/a for horizontally incident SH wave

    图  8  SH波水平入射时τθz*随着d/r的变化情况

    Figure  8.  Distribution of DSCF with d/r for horizontally incident SH wave

    图  9  SH波水平入射时τθz*随着k1*μ1*的分布

    Figure  9.  Distribution of DSCF with k1* and μ1* for horizontally incident SH wave

  • [1] PAO Y H, MOW C C, Achenbach J D. Diffraction of elastic waves and dynamic stress concentrations[J]. Journal of Applied Mechanics, 1973, 40(4):213-219. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=JJ0220516407
    [2] COUSSY O. Scattering of SH-waves by a cylindrical inclusion presenting an interface crack[J]. CR Acad Sci Paris, 1982, 295:1043-1046. https://www.researchgate.net/publication/309224500_Scattering_of_SH_Wave_by_a_Cylindrical_Inclusion_and_a_Semi-Cylindrical_Hollow_Near_Vertical_Interface_Crack_in_the_Bi-Material_Half_Space
    [3] COUSSY O. Scattering of elastic waves by an inclusion with an interface crack[J]. Wave Motion, 1983, 6(3):223-236. doi: 10.1016-0165-2125(84)90026-X/
    [4] 汪越胜, 王铎.SH波对有部分脱胶衬砌的圆形孔洞的散射[J].力学学报, 1994, 26(4):462-469. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400323665

    WANG Yuesheng, WANG Duo. Diffraction of SH-waves by a circular cavity with a particular debonded liner[J]. Acta Mechanica Sinica, 1994, 26(4):462-469. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400323665
    [5] 汪越胜, 王铎.剪切波作用下圆弧形界面裂纹的动强度因子[J].固体力学学报, 1993, 14(4):362-367. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000000699056

    WANG Yuesheng, WANG Duo. Dynamic stress intensity factors of circular arcshaped interface cracks subjected to shear waves[J]. Acta Solida Mechanics Sinica, 1993, 14(4):362-367. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000000699056
    [6] LEE V W, KARL J. Diffraction of SV waves by underground, circular, cylindrical cavities[J]. Soil Dynamics and Earthquake Engineering, 1992, 11(8):445-456. doi: 10.1016/0267-7261(92)90008-2
    [7] QI Hui, YANG Jie, SHI Yong. Scattering of SH-wave by cylindrical inclusion near interface in bimaterial half-space[J]. Journal of Mechanics, 2011, 27(1):37-45. doi: 10.1017/jmech.2011.5
    [8] 齐辉, 丁晓浩, 赵元博.双相介质弹性半空间垂直界面附近椭圆形夹杂对SH波的散射[J].振动与冲击, 2015, 34(24):76-81. http://d.old.wanfangdata.com.cn/Periodical/zdycj201524013

    QI Hui, DING Xiaohao, ZHAO Yuanbo. Scattering of SH-wave by an elliptic inclusion near vertical interface in a bi-material half-space[J]. Journal of Vibration and Shock, 2015, 34(24):76-81. http://d.old.wanfangdata.com.cn/Periodical/zdycj201524013
    [9] 杨杰.半空间垂直界面裂纹及圆夹杂对SH波的散射[D].哈尔滨: 哈尔滨工程大学, 2012: 22-34. http://cdmd.cnki.com.cn/Article/CDMD-10217-1012518460.htm
    [10] LOEBER J F, SIH G C. Diffraction of antiplane shear waves by a finite ccrack[J]. Journal of the Acoustical Society of America, 1967, 44(1):90-98. doi: 10.1121/1.1911091
    [11] 陈冬妮, 齐辉.SH波入射覆盖层半空间圆孔及圆夹杂的动力分析[J].哈尔滨工程大学学报, 2014, 35(2):171-176. http://d.old.wanfangdata.com.cn/Periodical/hebgcdxxb201402006

    CHEN Dongni, QI Hui. Dynamic analysis for circular cavity and inclusion impacted by SH-wave in a layered half-space[J]. Journal of Harbin Engineering University, 2014, 35(2):171-176. http://d.old.wanfangdata.com.cn/Periodical/hebgcdxxb201402006
    [12] 张希萌, 齐辉, 项梦.半空间双相压电介质垂直边界附近圆孔对SH波的散射[J].爆炸与冲击, 2017, 37(4):591-599. http://www.bzycj.cn/CN/abstract/abstract9758.shtml

    ZHANG Ximeng, QI Hui, XIANG Meng. Scattering of SH-wave by a circular cavity near the vertical boundary[J]. Explosion and Shock Waves, 2017, 37(4):591-599. http://www.bzycj.cn/CN/abstract/abstract9758.shtml
    [13] 杨在林, 许华南, 黑宝平.半空间椭圆夹杂与裂纹对SH波的散射[J].振动与冲击, 2013, 32(11):56-61. doi: 10.3969/j.issn.1000-3835.2013.11.012

    YANG Zailin, XU Huanan, HEI Baoping. Interaction of elliptical and crack under incident SH-wave in a half-space[J]. Journal of Vibration and Shock, 2013, 32(11):56-61. doi: 10.3969/j.issn.1000-3835.2013.11.012
    [14] 赵嘉喜, 齐辉, 杨在林.含有部分脱胶的浅埋圆夹杂对SH波的散射[J].岩土力学, 2009, 30(5):1297-1302. doi: 10.3969/j.issn.1000-7598.2009.05.018

    ZHAO Jiaxi, QI Hui, YANG Zailin. Scattering of SH-waves by a shallow buried cylindrical inclusion with a partially debonded curve in half space[J]. Rock and Soil Mechanics, 2009, 30(5):1297-1302. doi: 10.3969/j.issn.1000-7598.2009.05.018
    [15] 林宏, 刘殿魁.半无限空间中圆形孔洞周围SH波的散射[J].地震工程与工程振动, 2002, 22(2):9-16. doi: 10.3969/j.issn.1000-1301.2002.02.002

    LIN Hong, LIU Diankui. Scattering of SH-wave around circular cavity in half-space[J]. Earthquake Engineering and Engineering Vibration, 2002, 22(2):9-16. doi: 10.3969/j.issn.1000-1301.2002.02.002
    [16] 赵元博, 齐辉, 丁晓浩, 等.SH波作用下地表软覆盖层中圆形夹杂的动应力分析[J].爆炸与冲击, 2017, 37(6):982-989. http://www.bzycj.cn/CN/abstract/abstract9823.shtml

    ZHAO Yuanbo, QI Hui, DING Xiaohao, et al. Dynamic analysis for shallow buried circular inclusion impacted by SH-wave in a softlayered half-space[J]. Explosion and Shock Waves, 2017, 37(6):982-989. http://www.bzycj.cn/CN/abstract/abstract9823.shtml
    [17] 陈志刚.平面SH波在浅埋椭圆柱形孔洞上的散射与应力集中[J].暨南大学学报, 2005, 26(3):324-330. doi: 10.3969/j.issn.1000-9965.2005.03.009

    CHEN Zhigang. Scattering of SH-wave by a shallow buried elliptic cylindrical cavity and dynamic stress concentration[J]. Journal of Jinan University, 2005, 26(3):324-330. doi: 10.3969/j.issn.1000-9965.2005.03.009
  • 期刊类型引用(1)

    1. 王晔,白春华,刘文杰. 装置参数对125 kg燃料成雾性能的影响. 爆炸与冲击. 2025(05): 33-41 . 本站查看

    其他类型引用(5)

  • 加载中
推荐阅读
城市地下排水管道中燃气爆炸及气-液两相耦合作用规律
周刚 等, 爆炸与冲击, 2024
基于高压气体驱动的爆炸波模拟激波管冲击波衰减历程控制方法
程帅 等, 爆炸与冲击, 2024
Cf3i和co2抑制甲烷-空气爆炸实验研究
程方明 等, 爆炸与冲击, 2022
氢气-甲烷-乙醇混合燃料的爆炸压力特性
郭宏展 等, 爆炸与冲击, 2023
液化天然气罐式集装箱安全阀起跳预防和应急措施
常心洁 等, 集装箱化, 2021
H型巷道内采用不同布置方式的双爆源瓦斯爆炸传播特性
叶青 等, 高压物理学报, 2024
爆炸作用下综合管廊燃气舱内同舱管道失效的影响因素
米红甫 等, 高压物理学报, 2022
Cd36 as a gatekeeper of myocardial lipid metabolism and therapeutic target for metabolic disease
Glatz, Jan F. C., PHYSIOLOGICAL REVIEWS, 2024
Study on the minimum safe thickness of water inrush prevention in karst tunnel under the coupling effect of blasting power and water pressure
TUNNELLING AND UNDERGROUND SPACE TECHNOLOGY, 2024
Study on internal rise law of fracture water pressure and progressive fracture mechanism of rock mass under blasting mpact
TUNNELLING AND UNDERGROUND SPACE TECHNOLOGY
Powered by
图(9)
计量
  • 文章访问数:  5786
  • HTML全文浏览量:  1470
  • PDF下载量:  33
  • 被引次数: 6
出版历程
  • 收稿日期:  2017-05-02
  • 修回日期:  2017-08-02
  • 刊出日期:  2018-11-25

目录

/

返回文章
返回