Dynamic responses of hollow steel pipes directly buried in high-saturated clay to blast waves
-
摘要: 设计和实施了系列爆炸波作用下钢管的动态响应实验,获得了管道应变、振动速度和加速度及地表振动速度的时间历程。分析实验数据可知:在爆炸波近中场,峰值应变大小同管土相对刚度因数负相关,并随比例距离呈幂函数形式衰减,且在爆炸波不同分区衰减指数不同;地表峰值粒子速度、各管的峰值振动速度和测点峰值应变之间具有良好的线性关系。对各测试量做FFT (Fast Fourier Transformation)频谱分析得:振动信号主要集中在低频段,质心频率在10~60 Hz范围内,频谱同天然地震波谱有明显差异,动应变谱质心频率随药量增加亦呈幂函数形式衰减;考虑了爆腔因素以后,管道和地表振速频谱的质心频率和比例距离取对数后具有线性衰减关系。实验所得数据可应用于相似条件下管道的抗震计算,一些结论可作为深入研究爆炸波作用下埋地管道冲击振动机理的依据。Abstract: A series of experiments were designed and implemented to explore dynamic responses of steel pipes to blast waves The time histories of dynamic strains, vibration velocities and accelerations of the steel pipes were obtained, and the vibration velocity-time curves were gained. It is known by analyzing the experimental data that in the near and middle fields of the blast wave, the peak dynamic strains are negatively correlated with the relative stiffness coefficient of the pipe and the soil, and they follow the attenuation law of the power function with scaled distance; that the attenuation indexes are different for the blast wave propagating in the different sections of the field. The peak particle vibration velocities of the ground and the pipes have good linear correlations with the peak strains of the measuring points at the pipes. On the basis of the spectrum analysis by the fast Fourier transform on each test quantity, it is found that the spectrum energy of each test quantity is mainly concentrated in the low-frequency band, the centroid frequency is in 10−60 Hz, but there are obvious differences compared to the spectra of natural seismic waves. The centroid frequency of the dynamic strain spectrum is decayed in the exponent form of a power function with the increase of the explosive charge. By taking logarithm, there is a linear attenuation relationship between the scaled distance and the centroid frequencies of velocity spectra that the blasting cavity factor is considered. The test data can be directly applied to the seismic calculation under similar conditions, and some conclusions can be used as the theoretical basis for further study of the impact damage mechanism of buried pipelines.
-
Key words:
- blast wave /
- spectrum analysis /
- buried pipe /
- dynamic strain /
- vibration velocity /
- vibration acceleration
-
港口码头是交通物流的重要枢纽,也是国家海洋经济的主要依托,由于其重要的经济战略意义和开放的地理环境,易受到多形式打击破坏。水下爆炸是主要打击形式之一,具有毁伤作用强、毁伤范围大、隐蔽性高的特点。沉箱码头工作性能和耐久性好,在我国应用广泛,开展水下爆炸下沉箱码头毁伤效应研究,具有重要的研究价值和工程指导意义。
水下爆炸结构毁伤的研究最先开始于海军装备的研究,早期以舰船、潜艇等为研究对象。Rajendran[1]、吴林杰等[2]、Wang等[3]简化舰船结构,对钢板、钢梁水下爆炸破坏过程进行了系统的试验及数值研究,得到了钢板、钢梁的破坏形态和毁伤过程。Wardlaw等[4]分析了近场水下爆炸流固耦合作用,发现气泡脉动是结构毁伤的重要因素,特别是水下接触爆炸气泡脉动不容忽视。周章涛等[5]结合试验和数值模拟软件,考虑结构表面空化和气泡膨胀等因素研究了水下接触和近距爆炸作用下钢板的加载机理和破坏过程。针对混凝土水工结构,徐强等[6]、孙金山等[7]对大坝、桥梁开展结构水下抗爆研究,张社荣等[8]、王高辉等[9]通过数值模拟软件考虑水深等因素,研究了混凝土坝的毁伤机理和破坏形态。基于三峡围堰拆除工程,刘美山等[10]进行了一系列水下爆破试验,获得了混凝土水下爆炸的破坏形态。国内外就水下爆炸混凝土结构毁伤效应已开展一些研究,但关于沉箱码头水下接触爆炸的相关研究较少,为探究水下接触爆炸作用下沉箱码头的破坏过程和毁伤机理,数值模拟是必不可少的研究手段。
本文中在试验研究基础上,通过LS-DYNA有限元软件揭示水下接触爆炸作用下沉箱码头毁伤过程,分析冲击波和气泡对沉箱码头的荷载规律,研究炸药深度对沉箱码头毁伤效应的影响,探析沉箱码头的毁伤机理和毁伤特征。
1. 有限元模型
依据沉箱码头模型水下爆炸试验[11],图1所示为沉箱码头水下接触爆炸试验示意图,试验场地为直径8 m、深8 m的钢筋混凝土空心圆柱爆坑,内衬2 cm厚钢板,设计水深1.8 m,模型四面临水,底部无人为约束。采用1 kg TNT当量的圆柱状PETN炸药,药包轴线垂直于迎爆面,在0.9 m水深接触沉箱码头爆炸。沉箱码头模型长2.98 m、宽1.62 m、高2.19 m,其中沉箱结构高1.8 m,由6个长86 cm、宽65 cm仓格构成,仓格内填满饱和砂。沉箱底板厚25 cm,沉箱上部依次为管沟和面板,高39 cm,沉箱码头具体尺寸和配筋情况如表1所示。沉箱混凝土抗压强度为35.0 MPa,码头上部结构混凝土抗压强度为28.2 MPa,钢筋采用HPB335型号。考虑码头结构对称性,采用LS-DYNA软件建立沉箱码头对称模型如图2所示,有限元模型包括沉箱钢筋混凝土结构、炸药、水和空气,其中混凝土、钢筋和远场仓格土C2~C6采用Lagrange单元,网格尺寸为2~3 cm;考虑炸药附近流场和仓格土体发生大变形,炸药、水、空气和近场仓格C1土采用Euler单元,网格尺寸为1.5 cm~2 cm,小于Lagrange单元网格,两者均在爆心加密向四周渐变划分。通过关键字Ale_Coupling_Nodal_Constraint定义混凝土、钢筋之间的粘结,关键字Constrained_Lagrange_In_Solid定义流场与结构之间的流固耦合作用。考虑水下爆炸气泡脉动受重力和浮力影响,通过关键字Initial_Stress_Depth和Load_Body_Z初始化静水压力和重力,并在对称面施加对称约束,流场边界采用环境单元(Ambient)保证流场压力流出。在仓格C1外墙迎爆面设置测点,水平方向测点为H1和H2,垂直方向测点为V1~V7,测点之间的距离为20 cm,仓格C1内墙中心位置设置测点V8。
表 1 主要部位混凝土厚度及配筋情况Table 1. Concrete thickness and matching bar conditions of main parts位置 混凝土厚度/cm 配筋情况 保护层厚度/cm 仓格外墙 12 双层双向配筋,钢筋直径1.2 cm,间距18 cm 2.0 仓格内隔墙 8 双层双向配筋,钢筋直径0.8 cm,间距9 cm 1.5 沉箱底板 25 双层双向配筋,钢筋直径2.0 cm,间距18 cm 4.0 管沟底板 13 双层双向配筋,钢筋直径0.6 cm,间距15 cm 2.0 管沟外壁 12 双层双向配筋,钢筋直径0.6 cm,间距15 cm 1.5 面板 6 管沟上部面板单层双向配筋,其他部位不配筋 1.5 封仓板 6 不配筋 采用正确的材料参数对数值结果的准确性十分重要,混凝土采用Concrete_Damage_Rel3模型[12],该材料模型引入初始屈服面、极限强度面和残余强度面这3种强度面,考虑偏应力不变量对强度破坏面的影响。Mat_Concrete_Dmage_Rel3混凝土材料模型中,当应力达到初始屈服面,但未达到极限强度面时,通过初始屈服面和极限强度面线性插值表示;当应力达到极限强度面,但未达到残余强度面时,通过极限强度面和残余强度面线性插值表示:
Δσf=ηΔσm+(1−η)Δσyλ≤λm(后继屈服面) (1) Δσh=ηΔσm+(1−η)Δσrλ>λm(后继软化面) (2) 式中:λ是损伤变量,是等效塑性应变的函数;λm表示损伤转折点,是强化段和软化段的边界;η是损伤变量λ的函数,λ<λm时为强化段由0增至1,λm>λ为软化段由1减至0。
沉箱混凝土的抗压强度为35.0 MPa,码头上部结构混凝土的抗压强度为28.2 MPa,密度为2550 kg/m3,泊松比为0.2,通过动载增大系数κ曲线设定混凝土应变率效应,混凝土抗压动载增大系数κc和抗拉动载增大系数κt分别为:[13-14]。
κc=fdcfc={1˙ε<˙εstat(˙ε˙εstat)1.026α˙εstat<˙ε≤30s−1γ(˙ε˙εstat)1/330s−1<˙ε<300s−1 (3) κt=fdtft={(˙ε˙εstat2)β˙εstat2<˙ε≤1s−1ξ(˙ε˙εstat2)1/31s−1<˙ε<160s−1 (4) 式中:
˙ε 为动载应变率,s−1;˙εstat 为静态应变率,取3×10−5 s−1;fdc为混凝土在动载应变率为˙ε 下的动态抗压强度,MPa;fc为混凝土在静态应变率˙εstat 下的抗压强度,MPa;α=15+0.9fc,lgγ=6.156α−2 ;˙εstat2 为静态应变率,取10−6 s−1;fdt为混凝土在动载应变率为˙ε 下的动态抗拉强度,MPa;ft为混凝土在静态应变率˙εstat2 下的抗拉强度,MPa;β=11+0.8fc,lgξ=6β−2 ,30 MPa≤fc≤70 MPa。钢筋采用Mat_Plastic_Kinematic弹塑性模型,采用HRB335强度,泊松比为0.3,应变率参数C和P分别取40.4和5。空气采用线性多项式状态方程,C0~C6为状态方程参数,E为初始质量内能;水采用Grüneisen状态方程,C、S1~S3为状态方程参数,γ为Grüneisen常数;炸药采用标准JWL方程,A、B、ω、R1、R2为状态方程参数,黏土采用线弹性模型,E为弹性模量,G为剪切模量,参数设置如表2所示。
表 2 材料参数Table 2. Material parameters空气 ρ/(kg·m−3) C0−C3 C4 C5 C6 E/(J·kg−1) 1.29 0 0.4 0.4 0 250000 水 ρ/(kg·m−3) C S1 S2 S3 γ 1000 1480 2.56 −1.986 0.2268 0.5 炸药 ρ/(kg·m−3) A B ω R1 R2 1630 3.74×1011 7.33×109 0.3 4.15 0.95 黏土 ρ/(kg·m−3) E/MPa G/MPa 1800 16 8 2. 计算结果及分析
2.1 冲击波阶段
水下爆炸载荷主要由两个阶段构成:冲击波阶段和气泡脉动阶段。图3为冲击波阶段沉箱码头迎爆面区域流场与码头结构的压力云图,炸药引爆后在水中形成超高压冲击波,直接作用在沉箱外墙,根据压力云图可以发现:透射冲击波在仓格土中的衰减速度远大于初始冲击波在水中的衰减速度,仓格C1内墙迎爆面(测点V8)受到的压力峰值仅7.56 MPa,因此冲击波阶段仓格内墙损伤很小。冲击波沿外墙迎爆面向水域四周传播,t=650 μs时在水面、水底发生反射,t=1100 μs时水底反射冲击波传播到气泡位置再次发生反射透射,根据压力云图可以发现此时流场冲击波压力衰减很大。图4所示为沉箱外墙毁伤现象,在冲击波和爆轰产物的共同作用下外墙迎爆面混凝土压缩失效形成爆坑,初始冲击波在混凝土内部转变为压缩波传播,当冲击波传播到外墙背爆面时反射形成稀疏拉伸波,由于混凝土抗拉强度较弱,造成外墙背爆面拉伸破坏,因此沉箱外墙形成内外大、中间小的锥形爆坑,并且外墙迎爆面在压缩波的切向压缩和径向扩张下出现环状损伤裂缝。沉箱外墙顶部(外墙与上部管沟连接处)和底部(外墙与底板连接处)混凝土出现局部冲切损伤,码头上部面板未出现明显损伤。
设定混凝土单元的比例损伤变量δ:
δ=2λ/(λ+λm) (5) 式中:0<δ<1时,混凝土屈服,进入强化阶段;1<δ<2时,混凝土进入软化阶段。
选取沉箱外墙迎爆面炸药附近测点(V4、H1、V5)有效应力时程曲线如图5所示,在t=250 μs时有效应力达到峰值,分别为114.64、133.57、154.62 MPa,炸药下方测点峰值应力更大,这是由于炸药并不在外墙迎爆面中心,而是距离沉箱底板更近,因此炸药下方混凝土受到的冲击作用更大。炸药右侧测点(H1)由于临近仓格侧墙,应力波在外墙边界反射,应力在下降段出现明显振荡。V4、H1、V5测点速度时程曲线如图6所示,外墙迎爆面测点在冲击作用后迅速加速达到第一个峰值速度,外墙发生变形产生局部空化,空化的发生截断了压力载荷,外墙失去了加载载荷并开始减速,随着空化闭合出现第2个速度峰值,然后逐渐减速,迎爆面空化造成明显的二次加载现象。该现象对垂向测点(V4、V5)来说比较明显,这是由于冲击波作用下外墙中垂线方向变形明显,水平测点H1处变形较小,二次加载现象较弱。
查看冲击波阶段迎爆面测点混凝土的比例损伤变量δ如图7所示,由于接触爆炸超高的爆轰压力,炸药周围混凝土瞬间屈服并进入强化阶段,同时由于应变率效应,混凝土并没有立即失效,距离炸药最近的V4率先失效,随后H2、V5相继失效,其中V4、V5失效是由于近场混凝土压坏,而侧面H2失效则因为测点位于外墙和仓格侧墙连接处,混凝土发生剪切破坏失效。
2.2 气泡脉动阶段
图8为水下接触爆炸气泡与沉箱外墙耦合过程的损伤云图,其中灰色单元为已失效单元。冲击波阶段沉箱码头外墙为主要损伤区域,迎爆面形成内外大、中间小的破口(t=0.5 ms)。由于破口的产生,一部分爆轰产物涌入仓格内,直接作用于仓格内土体,外墙严重凹陷土体挤压导致仓格顶板向上膨胀,仓格与上部管沟连接处混凝土冲切破坏,管沟被直接顶起,码头面板由于反射拉伸应力波振荡,混凝土拉伸失效。气泡贴紧外墙膨胀,由于气泡上浮形成水冢,靠近外墙一侧水无法及时补充,t=20 ms时外墙附近水面形成空压区,靠近墙面水面凹陷。在气泡脉动阶段,由于气泡膨胀挤压沉箱外墙,沉箱外墙破口进一步增大,仓格顶板被向上顶起,在仓格变形和拉伸应力的共同作用下,靠近外墙的面板被直接掀起。
Cole[15]通过大量试验数据得到气泡半径经验公式:
Rmax=3.36(WZ)1/3 (6) 式中:Rmax为气泡膨胀最大半径,m;W为炸药质量,kg;Z为代表炸药位置的流体静压的等效水深,m,Z=H0+γ,H0为水面大气压的等效水深,H0=10.33;γ为炸药水深,m。
1 kg TNT在0.85 m水深的理论最大直径为3.01 m,而试验水深为1.7 m,因此气泡在t=23 ms冲出水面,气泡内压释放外界气体涌入,气泡溃灭,t=40 ms时码头毁伤基本不再发展。根据经验公式:
T=2.08W1/3Z5/6 (7) 式中:T为第一次气泡脉动周期,s。气泡第1次脉动周期为0.28 s,计算时长约为0.14倍的气泡第一次脉动周期,大于文献[2]评估舷侧仓外板的0.05倍气泡脉动周期。
3. 试验现象验证
接触爆炸下沉箱码头模型的毁伤现象如图9所示,码头迎爆面、侧面和顶面均出现一定破坏。沉箱迎爆面严重内凹炸药附近混凝土破坏严重,形成外径1.10 m×0.96 m、内径0.92 m×0.74 m的椭球形爆坑,数值模拟的破口外径1.07 m×0.95 m、内径1.03 m×0.86 m,结果较为吻合,外墙爆坑下部毁伤程度大于上部,这是由于炸药位置偏于外墙中心下方。外墙钢筋网架整体内凹,试验模型节点最大挠度为13 cm,数值模型为14.08 cm,基本一致。
沉箱顶部靠近炸药侧面板被掀起,沿着管沟方向形成通长裂缝并贯穿面板。码头侧面管沟下方形成45°斜裂缝,迎爆仓格侧墙上部与管沟连接处形成裂缝并分离,数值模拟结果的破坏形态与试验结果基本一致。
4. 炸深对毁伤效应的影响
水下接触爆炸时,炸药深度对气泡脉动形态和码头毁伤效应影响很大,调整炸深分别为0.5和1.3 m,分别距水面、水底0.5 m。比较气泡脉动形态和码头毁伤现象,如图10所示。近水面和近水底接触爆炸气泡分别在t=2 ms和t=48 ms冲出水面溃散,近水面爆炸时码头面板的毁伤增强,迎爆外墙上部形成椭圆形爆坑,毁伤由上向下逐渐衰减为环状裂缝,爆坑外径为0.59 m×0.91 m,沉箱外墙底部和码头侧面毁伤较弱。近水底爆炸时,由于沉箱底部为25 cm厚的钢筋混凝土底板,沉箱外墙毁伤效应较小,爆坑外径为0.40 m×0.74 m,迎爆仓格毁伤变形较小,码头面板沿管沟形成通长裂缝,没有完全掀起。近水面和近水底爆炸码头沉箱外墙的破口大小和侧面的毁伤现象均小于水域中部爆炸,近水面爆炸时码头面板的毁伤现象最严重。
根据水下接触爆炸作用下沉箱码头结构内能变化,可以得到不同炸深条件下码头各部分吸收能量占爆轰能量的比值,如表3所示。水下接触爆炸沉箱码头迎爆仓格外墙和钢筋相较远场仓格和码头面板吸收能量最多,近水面接触爆炸时码头面板吸收能量增加,其余部分吸收能量减少,近水底接触爆炸沉箱除迎爆仓格的其余仓格吸收能量轻微增加,其余部分吸收能量减少;各个工况下码头结构各部分吸收能量比例与码头破坏形态显示一致。因此,水域中部接触爆炸对沉箱的毁伤效应最强,近水面爆炸对码头上部面板的毁伤效应更强,近水底爆炸对沉箱码头的毁伤效应最弱。
表 3 码头结构各部分吸收能量Table 3. Energy absorption of different parts工况 爆轰能量/kJ 迎爆仓格外墙吸能/kJ 沉箱其余仓格吸能/kJ 码头面板吸能/kJ 钢筋吸能/kJ 水域中部爆炸 2560 89.68(3.5%) 27.77(1.09%) 0.47(0.2%) 124.85(4.88%) 近水面爆炸 2560 65.45(2.56%) 18.18(0.71%) 0.62(0.2%) 80.08(3.13%) 近水底爆炸 2560 68.63(2.68%) 28.29(1.11%) 0.30(0.1%) 95.38(3.73%) 注:括号内为吸收能量占总能量的百分比。 5. 结 论
基于LS-DYNA有限元软件建立水下接触爆炸沉箱码头全耦合数值模型,对水下接触爆炸下沉箱码头毁伤特性开展了研究,通过与试验结果对比,验证了有限元模型的准确性。分别从冲击波阶段和气泡脉动阶段研究了码头模型的破坏过程和毁伤机理,通过改变炸深,分析了不同炸深条件下沉箱码头的破坏形态,得到如下结论:
(1)通过有限元计算与模型试验结果的对比,码头各部位毁伤现象的数值计算结果与试验结果吻合良好,水下接触爆炸沉箱码头破坏过程可分为冲击波阶段和气泡脉动阶段,接触爆炸爆轰能量主要由沉箱外墙和钢筋吸收。
(2)水下接触爆炸冲击波阶段,沉箱外墙迎爆面出现空化及其闭合再加载的荷载特性,沉箱外墙形成内外大,中间小的锥形爆坑,外墙与上部管沟和沉箱底部的连接处混凝土出现不同程度的剪切破坏,仓格土对沉箱内部仓格有很好的缓冲作用。沉箱外墙毁伤在冲击波阶段已经基本形成,码头上部结构基本没有损伤。
(3)水下接触爆炸气泡脉动阶段,随着气泡膨胀爆轰产物直接涌入仓格,仓格外墙破口大小增加,仓格顶板向上膨胀并顶起管沟,导致上部码头面板破坏,仓格和沉箱侧墙产生拉伸裂缝,码头面板和沉箱侧墙的损伤在气泡脉动阶段逐步形成。
(4)随着炸深增加,气泡冲出水面溃散的时间延长,沉箱码头毁伤完成时间增加。水域中部接触爆炸对沉箱的毁伤作用最强,近水面接触爆炸对码头面板的毁伤作用更大,近水底接触爆炸对沉箱码头毁伤最弱。
-
表 1 管道几何参数
Table 1. Geometrical parameters of pipes
管道编号 D/mm δ/mm L/m S1 110 4.24 4.6 S2 160 4.70 4.6 S3 300 4.40 4.6 表 2 实验方案
Table 2. Experimental schemes
H/m R/m Q/g 2.0 2.2 50 75 100 125 150 175 200 2.7 50 75 100 125 150 175 200 1.5 2.2 50 75 100 125 150 175 200 2.7 50 75 100 125 150 175 200 3.2 50 75 100 125 150 175 200 1.0 2.7 50 75 100 125 150 175 200 表 3 实验参数
Table 3. Experimental parameters
H/m R/m Q/g ˉR/(m·kg−1/3) 管道编号 1.0 2.7 75 6.402 S3 1.0 3.2 125 6.4 S3 表 4 校验结果
Table 4. Results of check
测点编号 相关系数 可信概率 测点编号 相关系数 可信概率 S323H 0.836 0.872 S332Z 0.855 0.893 S311H 0.786 0.834 S341Z 0.897 0.910 表 5 S3管不同测点的最大峰值应变(R=2.2 m)
Table 5. The maximum peak strains at different measured points of pipe S3 (R=2.2 m)
Q/g 最大峰值应变/10−6 S321H S341Z S323H S323Z S331Z S332H S332Z S333H S333Z S334Z 50 −31.1 35.1 −50.7 −31.7 104.3 53.8 −66.3 48.3 44.5 −58.5 75 −32.1 39.4 −51.8 −35.8 103.9 62.7 −75.4 43.8 50.4 −65.9 100 −30.7 52.9 −63.7 −35.9 97.3 61.1 −77.8 45.2 52.9 −61.6 125 −26.4 44.3 −45.2 −33.3 122.5 74.5 −81.2 60.3 54.3 −58.2 150 −42.7 60.7 −98.1 56.5 174.9 90.6 −117.7 69.6 76.4 −78.1 175 −78.5 72.2 −71.7 −64.3 187.8 93.6 −118.6 79.9 79.2 −79.3 200 −53.9 86.3 −89.2 −84.1 228.9 103.4 −131.3 100.5 90.4 −87.4 注:压应变按绝对值最大给出,下同。 表 6 S3管不同测点的最大峰值应变(R=2.7 m)
Table 6. The maximum peak strains at different measured points of pipe S3 (R=2.7 m)
Q/g 最大峰值应变/10−6 S321H S341Z S323H S323Z S331Z S332H S332Z S333H S333Z S334Z 50 −16.5 7.2 −33.7 −14.4 −26.1 44.5 14.6 21.0 7.6 −29.8 75 −22.4 10.4 −52.0 −20.2 −34.9 61.1 19.4 28.1 11.2 −28.2 100 −28.2 14.9 −67.2 −28.3 −49.4 80.7 28.8 36.2 17.9 −31.3 125 −30.5 17.6 −69.6 −30.9 −58.3 88.7 33.6 39.1 23.9 −42.0 150 −36.5 21.0 −91.0 −37.1 −71.3 102.9 42.9 48.2 30.1 −48.4 175 −38.2 33.2 −103.2 −48.5 −94.1 121.3 52.3 60.1 34.5 −57.0 200 −42.1 28.9 −103.7 −44.6 −85.7 118.9 49.3 54.5 34.3 −48.5 表 7 S3管不同测点的最大峰值应变(R=3.2 m)
Table 7. The maximum peak strains at different measured points of pipe S3 (R=3.2 m)
Q/g 最大峰值应变/10−6 S321H S341Z S323H S323Z S331Z S332H S332Z S333H S333Z S334Z 50 −30.2 18.6 −36.9 −22.1 −49.1 39.3 −26.9 27.6 −15.5 34.4 75 −16.2 17.4 −32.0 −22.3 −45.9 31.8 −25.5 26.2 −16.8 32.8 100 −18.6 22.7 −36.6 −26.8 −52.8 35.4 −32.9 26.5 −19.6 33.3 125 −19.1 22.5 −38.9 −32.7 −34.1 45.3 −30.7 31.1 −20.9 37.9 150 −19.3 24.7 −41.5 −37.5 −61.1 52.3 −35.1 33.0 −22.1 38.2 175 −23.0 28.0 −43.0 −36.6 −71.8 53.0 −41.3 34.8 −28.6 40.9 200 −25.9 29.5 −47.7 −38.7 −90.8 58.1 −43.0 38.7 −38.1 42.5 表 8 S3管道振动的峰值粒子速度和地表振动的峰值粒子速度(R=2.7 m, H=2.0 m)
Table 8. Peak particle velocities of S3 pipe and ground vibrations (R=2.7 m, H=2.0 m)
Q/g 峰值粒子速度/(cm∙s−1) 比值/% Q/g 峰值粒子速度/(cm∙s−1) 比值% S3管道振动 地表振动 S3管道振动 地表振动 50 4.3 13.8 31 150 9.4 25.1 37 75 3.7 11.4 32 175 12.1 30.6 39 100 4.9 15.0 33 200 14.5 35.2 41 125 5.2 13.7 38 表 9 管道S3和地面振动速度质心频率(R=2.7 m, H=2.0 m)
Table 9. Centroid frequencies of pipe S3 and ground vibration velocities (R=2.7 m, H=2.0 m)
Q/g fc/Hz Q/g fc/Hz 管道 地表 管道 地表 X方向 Y方向 Z方向 X方向 Y方向 Z方向 X方向 Y方向 Z方向 X方向 Y方向 Z方向 50 24.3 18.0 27.5 18.3 22.3 18.3 150 22.6 26.9 28.8 29.7 16.0 29.8 75 25.1 19.7 25.8 14.9 19.6 16.4 175 20.4 23.0 22.7 32.9 16.5 39.9 100 25.4 22.4 28.8 19.1 18.0 15.5 200 33.2 24.5 26.9 36.3 15.7 45.8 125 26.2 26.2 36.8 28.6 17.0 15.2 表 10 S2管道振动的峰值加速度(H=1.5 m, R=2.7 m)
Table 10. Peak particle vibration acceleration of pipe S2 (H=1.5 m, R=2.7 m)
Q/g 加速度/(m·s−2) X方向 Y方向 Z方向 合成 Max Min Max Min Max Min 50 4.7 −11.1 3.7 0.3 6.7 −3.7 13.5 75 6.7 −12.9 4.0 0.0 8.1 −3.5 15.8 100 4.1 −10.4 2.2 −1.2 7.7 −2.5 13.2 125 4.7 −16.1 1.5 −2.9 8.4 −7.0 17.8 150 3.1 −15.5 0.9 −3.1 6.7 −6.5 17.1 175 1.7 −14.7 1.0 −3.3 6.8 −7.1 16.6 200 3.0 −18.2 0.3 −4.0 6.9 −7.4 20.1 表 11 加速度质心频率(H=1.5 m, R=2.7 m)
Table 11. Centroid frequencies of acceleration (H=1.5 m, R=2.7 m)
Q/g fc/Hz X方向 Y方向 Z方向 50 55.9 40.1 32.9 100 41.3 41.9 33.0 200 41.1 39.1 36.8 表 12 S3管加速度、速度和位移峰值(H=2.0 m, R=2.2 m)
Table 12. Peak acceleration, peak velocity and peak displacement of S3 (H=2.0 m, R=2.2 m)
Q/g dmax/mm vmax/(cm·s−1) amax/(m·s−2) (vmax⋅d−1max)/Hz (amax⋅v−1max)/Hz 50 0.99 7.45 12.73 75.60 170.85 75 0.82 5.89 11.23 71.69 190.78 100 1.02 7.21 15.17 70.50 210.32 125 1.33 9.33 23.59 70.40 252.82 150 1.57 10.19 20.53 64.92 201.37 175 1.88 11.85 25.31 62.91 213.55 200 2.54 16.68 41.92 65.58 251.26 -
[1] 李强, 陈德利, 屈洋. 爆破对输气管道本体影响的监测 [J]. 油气储运, 2015, 34(2): 190–194. DOI: 10.6047/j.issn.1000-8241.2015.02.017.LI Q, CHEN D L, QU Y. Monitoring on effects of blasting on gas pipeline body distortion [J]. Oil & Gas Storage and Transportation, 2015, 34(2): 190–194. DOI: 10.6047/j.issn.1000-8241.2015.02.017. [2] 程围峰, 梁旭, 王振宇. 隧道爆破施工对邻近输油管道的影响评价 [J]. 石油工程建设, 2011, 37(4): 44–46; 9. DOI: 10.3969/j.issn.1001-2206.2011.04.013.CHENG W F, LIANG X, WANG Z Y. Assessment of blasting impact in tunnel excavation on nearby oil pipeline [J]. Petroleum Engineering Construction, 2011, 37(4): 44–46; 9. DOI: 10.3969/j.issn.1001-2206.2011.04.013. [3] 张紫剑, 赵昌龙, 张黎明, 等. 埋地管道爆破振动安全允许判据试验探究 [J]. 爆破, 2016, 33(2): 12–16. DOI: 10.3963/j.issn.1001-487X.2016.02.003.ZHANG Z J, ZHAO C L, ZHANG L M, et al. Experimental investigation of blasting vibration safety criterion on buried pipeline [J]. Blasting, 2016, 33(2): 12–16. DOI: 10.3963/j.issn.1001-487X.2016.02.003. [4] 戴联双, 张海珊, 孟国忠, 等. 在役油气管道周边爆破作业风险分析 [J]. 油气储运, 2012, 31(11): 801–803; 887. DOI: 10.6047/j.issn.1000-8241.2012.11.001.DAI L S, ZHANG H S, MENG G Z, et al. Risk analysis of blasting operations around in-service oil and gas pipelines [J]. Oil & Gas Storage and Transportation, 2012, 31(11): 801–803; 887. DOI: 10.6047/j.issn.1000-8241.2012.11.001. [5] 王振洪, 侯雄飞, 边明, 等. 爆破对天然气长输管道振动影响的安全判据 [J]. 油气储运, 2016, 35(8): 813–818. DOI: 10.60 47/j.issn.1000-8241.2016.08.003.WANG Z H, HOU X F, BIAN M, et al. Safety criterion for impacts of bursting on long-distance gas pipeline [J]. Oil & Gas Storage and Transportation, 2016, 35(8): 813–818. DOI: 10.60 47/j.issn.1000-8241.2016.08.003. [6] SISKIND D E, STAGG M S, WIEGAND J E, et al. Surface mine blasting near pressurized transmission pipelines [R]. Minneapolis, MN: US Department of the Interior Bureau of Mines, 1994: 2−3; 22−23. [7] DOWDING C H. Blast vibration monitoring and control [M]. New Jersey: Prentice-Hall, 1985: 167−171. [8] ESPARZA E D, WESTINE P S, WENZEL A B. Pipeline response to buried explosive detonations: Vol. 2 [R]. Arlington, USA: American Gas Association, 1981: 146−152. [9] KOURETZIS G P, BOUCKOVALAS G D, GANTES C J. Analytical calculation of blast-induced strains to buried pipelines [J]. International Journal of Impact Engineering, 2007, 34(10): 1683–1704. DOI: 10.1016/j.ijimpeng.2006.08.008. [10] ABEDI A S, HATAF N, GHAHRAMANI A. Analytical solution of the dynamic response of buried pipelines under blast wave [J]. International Journal of Rock Mechanics and Mining Sciences, 2016, 88: 301–306. DOI: 10.1016/j.ijrmms.2016.07.014. [11] 钟冬望, 龚相超, 涂圣武, 等. 高饱和黏土中爆炸波作用下直埋聚乙烯管的动力响应 [J]. 爆炸与冲击, 2019, 39(3): 033102. DOI: 10.11883/bzycj-2017-0334.ZHONG D W, GONG X C, TU S W, et al. Dynamic responses of PE pipes directly buried in high saturated clay to blast wave [J]. Explosion and Shock Waves, 2019, 39(3): 033102. DOI: 10.11883/bzycj-2017-0334. [12] 冯双, 马郧, 徐光黎. 武汉地区基坑一般黏性土的力学特性与其物理指标间的相关性研究 [J]. 安全与环境工程, 2016, 23(5): 149–154. DOI: 10.13578/j.cnki.issn.1671-1556.2016.05.025.FENG S, MA Y, XU G L. Correlation between mechanical properties and physical indicators of general cohesive soil in foundation pits in Wuhan area [J]. Safety and Environmental Engineering, 2016, 23(5): 149–154. DOI: 10.13578/j.cnki.issn.1671-1556.2016.05.025. [13] 库特乌佐夫. 爆破工程师手册 [M]. 刘清泉, 译. 北京: 煤炭工业出版社, 1992: 137–138 [14] 卢文波, 张乐, 周俊汝, 等. 爆破振动频率衰减机制和衰减规律的理论分析 [J]. 爆破, 2013, 30(2): 1–6; 11. DOI: 10.3963/j.issn.1001-487X.2013.02.001.LU W B, ZHANG L, ZHOU J R, et al. Theoretical analysis on decay mechanism and law of blasting vibration frequency [J]. Blasting, 2013, 30(2): 1–6; 11. DOI: 10.3963/j.issn.1001-487X.2013.02.001. [15] 陈勇, 杨庆华. 埋地管道在地震波作用下的变形分析 [J]. 世界地震工程, 2015, 31(4): 275–279.CHEN Y, YANG Q H. Deformation analysis of buried pipelines under seismic waves [J]. World Earthquake Engineering, 2015, 31(4): 275–279. 期刊类型引用(6)
1. 吉新博,卢卫建,吕琛,王玺,廉政,何丽. 不同冲击载荷下典型路面结构的动态弯沉响应. 兵工学报. 2024(01): 299-311 . 百度学术
2. 王继伟,曲占庆,郭天魁,陈铭,吕明锟,王旭东. 基于CDEM的页岩甲烷原位燃爆压裂数值模拟. 中国石油大学学报(自然科学版). 2023(01): 106-115 . 百度学术
3. 牛犇,张新伟,周玉,李婧,徐兴全,张一鸣. 基于连续-非连续元降雨工况三维边坡稳定性分析. 山东大学学报(工学版). 2023(01): 92-99 . 百度学术
4. 王涛,程垒明,项远铠,承宁,王博,周航. 砾石影响下裂缝扩展规律数值模拟. 新疆石油天然气. 2023(01): 42-48 . 百度学术
5. 路扬,冯春,程鹏达,张一鸣. 车轮-地面耦合动力学行为的CDEM-DEM分析. 光学精密工程. 2023(05): 746-756 . 百度学术
6. 张震东,高原,马大为,朱忠领,王玺. 多次模拟弹射载荷作用下水泥混凝土路面动态响应. 兵工学报. 2022(07): 1706-1717 . 百度学术
其他类型引用(3)
-