Scattering of SH-wave by elliptical inclusion with partial debond curve and circular cavity in half space
-
摘要: 采用复变函数法,结合"保角映射"技术及Green函数法,研究SH波作用下半空间内含有部分脱胶的椭圆夹杂以及圆形孔洞的散射问题。首先,利用"保角映射"技术将椭圆夹杂映射为圆夹杂,求出散射波位移场,同时,利用Green函数法与"虚设点源"的方法,求出半空间内椭圆夹杂以及圆孔的位移及应力场;然后,根据椭圆夹杂周围位移、应力连续、圆孔周围应力自由的边界条件,建立无穷线性代数方程组,求解出波函数中的未知系数;最后,在脱胶部分施加大小相等、方向相反的应力,构造出"脱胶模型",得到半空间内含有部分脱胶的椭圆夹杂以及圆形孔洞的总位移场。数值算例表明,入射角度、入射波频率、缺陷之间的距离、夹杂埋深及脱胶角度等对动应力集中因子有较大影响。Abstract: In order to study the scattering problems of SH-wave by elliptical inclusion with partial debond curve and circular cavity in half space. The method of "conformal mapping" is used to map the elliptical inclusions into circular inclusions. The displacement field and the stress field of the elliptical inclusion and circular cavity are obtained by Green's function method with the ridding of image method. Secondly, infinite system of linear equations are established by boundary conditions, which are continuous around the elliptical inclusion with the displacement and stress, and then the unknown coefficients of wave function are solved by the free around circular cavity with stress. Finally, the "partial debond curve model" is constructed, the equal stress with opposite direction is applied in the partial debond curve. By which we obtain the total displacement field of the elliptical inclusion with a partially debond curve and circular in the half space. Numerical examples show that the dynamic stress concentration factor (DSCF) is influenced by the incident angle, the frequency of incident wave, the distance of the defect, the depth of inclusion and the partial debond curve angle.
-
Key words:
- SH wave /
- scattering /
- elliptic inclusion /
- partial debond curve /
- Green's function /
- conformal mapping /
- dynamic stress factor
-
当可燃气体、可燃液体的蒸气(或可燃粉尘)与空气混合并达到一定浓度时,遇到火源就会发生爆炸。这个能够发生爆炸的浓度范围,叫做爆炸极限,通常用可燃气体、蒸气或粉尘在空气中的体积分数来表示。可燃气体在空气(氧气)中的爆炸极限范围是众多学者关心的焦点,但关于作为描述可燃气体爆炸后毁伤效果的典型参数,如爆压、爆温、爆速等的报道较少[1-4]。这是因为目前对爆炸极限的研究多局限于两个方面:一方面是因为偏于实际生产应用的缘故,多数测试中只需知道可燃气体在空气中的爆炸极限范围、临界氧含量或者获得爆炸三角形图即可,而对具体的燃烧爆炸过程及其结果并不关心;另一方面是由于测试仪器的局限性,多采用小尺寸容器,如小型激波管、20 L球、圆柱形爆炸罐等,其测试手段较单一,多数只在耐压容器壁面安置一个压力传感器或温度传感器,所得数据有限[5-7]。正因为如此,目前对于气体爆炸的研究多采用放宽条件,改变初始温度、初始压力、当量体积分数等手段,以获得某种可燃气体或气体混合物较全面的爆炸特性参数[8-10]。大尺寸密闭空间内可燃气体爆炸过程更贴合于实际。但由于爆炸容器尺寸大,操作复杂,开展相关研究难度较大。本文中,通过在大尺寸密闭容器中开展天然气爆炸超压场的研究,以期获得大尺寸密闭空间内天然气爆炸超压的发展规律,丰富目前天然气燃爆威力的测试数据,为密闭空间内天然气爆炸危害的预防及毁伤能力评估提供数据支持。
1. 实验
1.1 实验装置
天然气-空气混合物的燃爆过程在容积为10 m3的爆炸罐内进行,爆炸罐示意图见图 1。为更清楚地记录爆炸罐内天然气-空气混合物爆炸后超压的发展过程,在爆炸罐内沿罐体轴向典型位置布置4个压力传感器,沿罐体径向典型位置布置3个压力传感器,以获得天然气-空气混合物燃爆后,爆炸罐内部的超压场状态,传感器布置见图 2。点火装置选用高能放电器,单次点火,点火能量为40 J;超压测试系统包括压力传感器、信号调理器、数据采集系统、信号线等。压力传感器为PCB公司的113B系列;信号调理器为PCB公司的信号调理器;数据采集系统为南汇科技虚拟仪器系统;高速摄像系统为Photron公司的NX100相机,实验中采样频率为1 000 s-1。
1.2 实验样品
选用华北地区的工业天然气为实验样品,其密度为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%。
1.3 实验方法
将不同量的天然气在容积为10 m3的爆炸罐内与空气进行均匀混合,测量各点天然气的体积分数,达到罐体内各位置处天然气体积分数相对均匀时,进行点火操作,混合过程见文献[11],利用压力测试系统记录相关测试点的压力数据。
2. 实验结果分析
2.1 近爆炸下限处天然气-空气混合物的爆炸超压场状态
超压状态场通常被用来评估受限空间内可燃气体点爆过程中的爆炸效果。对大尺寸密闭空间来说,其超压状态场与可燃气体的体积分数存在对应关系[11]。本次研究主要从近爆炸下限(5.4%)的天然气点爆过程入手,分析不同体积分数下的天然气爆炸超压状态场。
经过系列实验测试后发现,在实际天然气体积分数接近天然气爆炸下限(5.4%)处,开展天然气点爆实验时,能够获得3种超压曲线,且曲线状态随天然气实际体积分数的不断升高发生一系列的变化。在天然气爆炸下限附近选取5.5%、5.8%和6.5%等3种天然气体积分数进行实验,以轴向第一个传感器的信号为典型信号,来研究不同体积分数条件下压力传感器获取的不同超压曲线,如图 3所示。
由图 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]。
2.2 近爆炸下限处天然气-空气混合物的爆炸超压发展过程
以容积为10 m3的爆炸罐为研究对象,对其轴向的4个压力传感器(距爆源由近至远分别命名为OP1~OP4)的压力数据进行分析,典型结果如图 4(a)所示。当天然气体积分数为5.5%时,接近实验测得的爆炸下限(5.4%),因此,此爆炸超压发展曲线图为近爆炸极限时的临界压力发展趋势图。由图 4(a)可看出,经过滤波处理后,轴向上的压力传感器随着距离点火位置的远近,其超压峰值分别为82.5、32.9、23.4和15.1 kPa。距爆源最近的传感器测得的压力曲线较接近典型的冲击波超压曲线,其他3个传感器所测得的压力曲线都接近于燃烧波的压力曲线。这主要是由于点火点处天然气的体积分数较低,接近爆炸下限,点火初期,点火位置处的天然气-空气混合物被点燃,初始冲击波产生,但由于能量支持不够,未继续发展,使得后续的轴向传感器测得的压力信号较弱且随距离呈递减趋势。
仍以容积为10 m3的爆炸罐为研究对象,对其径向的3个压力传感器(距轴线由近至远分别为OP4~OP6)的压力数据进行分析,典型结果如图 4(b)所示。由图 4(b)可以看出,传感器距离爆源中心轴线越远时,爆炸超压峰值越大,但增幅不大。这是由于天然气-空气燃爆发生并沿爆炸罐体轴向传播的同时,也沿爆炸罐体径向传播,呈体积性发展趋势,距爆源一定距离后,整个燃爆体系传播过程已成整体化趋势,同一波阵面的压力数据基本相当,但由于壁面反射的影响,偏离轴线处压力可能略有升高。
在对天然气体积分数为5.5%的天然气-空气混合物燃爆超压曲线分析后,针对3种天然气体积分数情况下的天然气-空气混合物燃爆发展进行研究,得到其燃爆参数随轴向和径向的发展规律,如图 5所示。由图 5可以看出:当天然气体积分数接近爆炸下限时,天然气-空气混合物燃爆的最高压力即为前导冲击波的超压峰值,其值相对较低;随着天然气体积分数的升高,天然气燃爆的最高压力为前导冲击波过后燃烧波的峰值压力,且此压力值随着初始天然气体积分数的升高而增大。从空间发展角度来看,距爆源距离对天然气爆炸超压峰值影响不大。而对于豫驰时间来说,与超压发展规律类似,燃爆豫驰时间与距爆源距离关系不大,初始天然气体积分数是决定性因素。
3. 结论
天然气爆炸下限附近存在3种典型的超压状态:(1)当可燃气体的体积分数接近爆炸下限时,点火后只存在点火点周边气体燃烧产生的前导冲击波;(2)当可燃气体体积分数略高于爆炸下限时,点火后前导冲击波和后续燃烧波共存;(3)当可燃气体体积分数高于爆炸下限一定程度后,前导冲击波与后续燃烧波重叠。
经过系统实验发现,爆炸下限附近的爆炸超压峰值及燃爆豫驰时间主要取决于初始天然气体积分数,而与距爆源距离关系不大,这主要是由于大尺寸密闭容器内气体燃爆过程的体积性效果。
-
-
[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=QK199400323665WANG 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=QK000000699056WANG 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/zdycj201524013QI 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/hebgcdxxb201402006CHEN 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.shtmlZHANG 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.012YANG 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.018ZHAO 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.002LIN 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.shtmlZHAO 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.009CHEN 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)
-