爆破振动下小净距隧道围岩的力学行为分析

2023-05-19 来源：个人技术集锦

桥

RIDGEAND

TUNNELENGINEERING

隧工程

西部交通科技

WesternChinaCommunications

Science&Technology

文章编号:16734874(2010)02、030058005

爆破振动下小净距隧道围岩的力学行为分析

柴国辉,柴　旺

(11贵州省公路桥梁工程总公司,贵州　贵阳　550001;21贵州高速公路开发总公司,贵州　贵阳　550001)

摘　要:文章运用数值分析计算软件,对小净距隧道中既有隧道受邻近隧道爆破震动影响的动力响应特性进行分析,得出了爆破地震波作用下,隧道周边围岩应力、位移和振动速度的分布规律,为小净距隧道的设计与安全施工提供科学依据。关键词:小净距;隧道;爆破震动;数值模拟中图分类号:U452　　　　　　文献标识码:A

;

爆破震动作用下邻近洞室的响应是一个复杂的动力学问题。数学解析法虽然能给出准确可靠的分析结果,但数学处理上的困难制约了它的适用范围。目前,理论分析仍只能解决介质、结构线弹性,入射波为稳态波,边界条件较简单的情况。要分析更复杂的实际问题,则需借助实验或数值计算方法。计算机数值模拟技术的发展为这一问题的研究提供了很好的工具,很多学者已经做过这方面的研究。

作者简介

柴国辉(1975—),男,仡佬族,贵州铜仁人,工程师,主要从事公路工程施工管理工作。

本文结合贵州某高速公路小净距隧道开挖施工,应用数值模拟软件模拟开挖过程中隧道爆破施工对临近既有隧道的影响,得到隧道围岩应力场、位移场及爆破振动速度场的分布规律。

　西部交通科技58

1　爆破振动作用下影响邻近隧道稳定性的2　工程概况

某高速公路隧道进口端前方铁路隧道下穿,出口

端基本与铁路平行。隧道K2+880～K3+250,YK2+880～YK3+260为小净距,线间距4～1217m,其余为分离式,隧道洞轴线走向约196°,净空(10125×510)m,其中左线隧道全长2092m,纵坡坡度为-218%,隧道底高程为1473136～1686111m;右线隧道底高程为1473136～1684193m,全长2124m,纵坡坡度为-218%。隧道最大埋深231145m。

因素

(1)根据动力弹性理论,爆破时,岩体中产生的

应力与质点的振动速度成正比。在爆破地震波引起相同振速的情况下,低频波比高频波对洞室围岩的影响要大。目前,许多国家在制定或修订爆破震动安全判据时,普遍考虑了振速和频率两个指标的共同影响。

(2)地质条件因素。地质条件主要指相临隧道之间的地质材料与地质构造,是地震波的传播介质。研究表明,在相同爆破地震波的作用下,岩石越硬,则围岩中引起的峰值振动速度越小。从波动力学的角度讲,爆破震动对邻近洞室的影响是地震波在传播过程中与地质介质相互作用的结果。应力波参数是两者波阻抗的函数,不同的岩石对爆破地震波的阻尼效应不同,特别是当洞室围岩存在某些特殊地质构造(如断层软弱夹层,溶洞等)时,应力波的传播更受影响。

(3)隧道自身条件因素。隧道的横断面形状、被覆材料,衬砌结构、爆源的相对位置、洞径与爆破应力波长的相对尺寸等直接影响洞壁周边的应力分布,而洞径是影响围岩应力和振速的重要因素。在其它条件相同的情况下,随着隧道直径的增大,迎爆一侧主要由于应力波的反射拉伸作用,其动力集中因子增大,而背爆一侧没有明显的变化。长江水利水电研究院建议的允许爆破震动速度和破坏标准见表1。表1中的数据为重复爆破的速度值,单次爆破时,可适当加大,但≯50%。

表1　地下结构允许的爆破振动速度和破坏标准表

质点振动速度

(cm/s)510～10101010～20102010～30103010～60106010～9010

3　爆破对邻近隧道影响的数值模拟分析311数值计算模型计算模型由两个截面形状为马蹄形的隧道及其周围围岩组成。分析地下开挖时,模型边界应距离开挖洞室中心大约10倍直径的长度。隧道断面上半部分为半径3m的半圆,下半部分为3m×10m的矩形。两洞室中心相距20m,计算模型X方向范围为-70～70m,Y方向范围为0～30m,Z方向范

围为30～-60m。用FLAC进行动力计算,网格剖分尺寸受输入波动的最短波长控制。设网格最大尺寸为Δl,输入波动的最短波长为λ,则Δl必须小于(1/10～1/8)。另外,动力计算要求网格划分的尺寸要相对均匀,网格最大尺寸和最小尺寸的比例一般≯5。因此,划分模型为44160个单元,共50094个节点。

地下结构破坏现象

未衬砌的松散洞体有小的掉块

隧道原有裂缝有时扩大;破碎岩体有掉块;管道接头有细

微变位

隧道有大掉块,有时有小的塌落;岩柱有掉块衬砌出现裂缝;管道变形;顶板出现塌方地下建筑物或衬砌开裂;硬岩体裂缝严重扩张

图1　计算模型图

2010年第2、3期总第32期

312　岩体物理力学参数及力学模型

率。

315　模拟计算结果分析

不考虑岩体物理力学参数变化对计算的影响,计算模型中的岩体为单一岩石材料。参数取值为体积模量K=7173GPa,剪切模量G=5132GPa,抗拉强度σt=012459MPa,粘聚力C=01725MPa,体积密度ρ=2720kg/mCoulomb强度准则。

313　初始条件与边界条件

(1)初始地应力场。初始地应力场仅按自重应

-3

计算分为静力和动力计算。动力计算以静力计算为基础,考虑到实际施工过程,静力开挖一次就应施加一次动荷载,但计算量太大。因此,计算过程进行了简化,先完成静力开挖,再一次性施加动荷载。静力分析时,整个开挖分5步完成。沿Y轴方向开挖,每一步开挖6m,并都计算至平衡状态,最后将速度时程施加于开挖隧道的侧壁,进行动力计算,得到动力计算结果。31511　应力分析

,内摩擦角φ=70177°。模

型采用理想的弹塑性模型,屈服准则采用Mohr2

力场考虑。根据弹性力学原理,竖向应力σv=γH,水

σv,式中γ为岩体容重;H为埋置深度;平应力σh=k

μ/(1-μ);μ为泊松比。k为侧压力系数,k=

(2)静力计算边界条件。采用位移边界条件,

静力开挖和动力开挖均会在隧道顶底板和两侧壁产生拉应力集中。静力开挖后,最大拉应力区在底板,其拉应力最大值为01028MPa。施加动荷载后,迎爆侧壁和底板拉应力集中程度增加,尤其是迎爆侧受影响最明显,最大拉应力达到0104MPa,说明动荷载施加后,既有隧道围岩应力发生了重新分布。动力开挖后,最小主应力发生重新分布,在隧道两底角的应力集中程度增加,在迎爆一侧的圆弧区域内产生压应力集中,最大压应力从2119MPa增加到213MPa。

图2是静力开挖后的最大主应力分布图,可以看出:邻近隧道开挖后,既有隧道的两侧壁、底板和顶板均产生一定程度的拉应力集中,底板的拉应力集中程度最大,最大拉应力达到010346MPa,洞室两侧壁的拉应力集中值较小,几乎为0。图3是动荷载施加012s以后,隧道周围的最大主应力分布,可以看出:既有隧道的两侧壁、底板和顶板的拉应力集中程度均不同程度的增大,尤其是邻近开挖隧道一侧,拉应力集中程度增大最明显,最大拉应力达到了010425MPa,比静力开挖后底板的最大拉应力增大了2315%。

图5是动荷载施加012s以后,隧道周围的压应力分布情况,可以看出:应力同样出现了重分布,静力开挖时最大压应力出现在远离开挖隧道一侧的底角,而施加动荷载后,则最大压主应力除了出现在远离开挖隧道一侧的底角附近外,还出现在了邻近开挖隧道一侧的圆弧位置,最大值也由2147MPa增大

即模型的左右(X方向)边界,前后(Y方向)边界和底边界均施加位移约束条件,上边界取自由边界。(3)动力计算边界条件。FLAC3D求解动力问题的边界条件设置有远置人工边界条件和粘滞边界条件。远置人工边界就是在动力计算中沿用静力计算中使用的边界条件,但这样的边界对波起着完全反射的作用,计算中要求将模型的范围取得足够大,从而给计算造成很大的负荷。粘滞边界条件通过在边界的法线方向和切线方向设置独立的粘壶得以实现,以便吸收来自模型内部的入射波。采用粘滞边界条件,即在模型的左右边界,前后边界和下边界均施加粘滞边界条件,上边界作为自由边界。

314　爆破震动荷载的输入及阻尼的选取

FLAC计算中动荷载的输入方式有加速度时

程、速度时程、应力时程和位移时程等4种。输入

λsin(2π/T)t,最大幅的速度时程为正弦压缩波ν=

值λ取012m/s,振动频率取100Hz,地震波的持续时间为013s。爆破产生的地震波以动压力形式均布作用于隧道洞壁上,方向垂直于壁面。

FLAC中,采用了瑞雷阻尼和局部阻尼两种形

式的阻尼,计算中多采用瑞雷阻尼。设置瑞雷阻尼时,需要确定阻尼比和中间频率两个参数,其中阻

尼比按照经验取值(一般取0105);中间频率则取自振频率或输入波的主频,这里取输入波的主频

　西部交通科技60

到216MPa,同样是邻近开挖隧道一侧的应力变化最明显。

比较图6和图7,可以看出两者均在隧道周围产生明显剪应力集中,并呈反对称的形式分布,但是动力开挖的剪应力集中程度要明显大于静力开挖,尤其是右侧底角和右半圆弧部分两个位置,无论从剪应力集中的区域还是最大剪应力值来看均有明

2010年第2、3期总第32期

显的增大。静力开挖后最大剪应力值为01709MPa,动力开挖后,右侧底角的剪应力值为0182MPa,右侧圆弧部分的剪应力为1118MPa,分别增大1516%和6614%,增大程度明显。

计算结果表明,迎爆侧一面的振动速度峰值均比背爆侧一面的振动速度峰值大很多。同一截面中的最大振速峰值出现在迎爆侧直墙中上部分区域,该区域的振动速度峰值均达到10cm/s以上。根据表1的判别标准,围岩会发生较小的破坏,该区域是隧道截面中最危险的区域。迎爆侧起拱线以上部分和底角为次峰值振动区,也是次危险区。背爆侧峰值振动速度普遍较小,最大值只有318cm/s,≯5cm/s,不会发生破坏,因此,背爆侧比较安全。

在动力计算过程中分别对爆破震动荷载作用时间为0101s、0103s、0105s、0108s、011s、012s时的应力分布情况进行记录。记录结果显示

最大拉应力分布规律是:在邻近开挖隧道一侧直墙的中上部均会产生明显的拉应力集中现象,并且是隧道周围最大拉应力值出现的位置,其应力程度要明显大于远离开挖隧道一侧;随着作用时间的增长,最大拉应力值有增大的趋势,在0103s时最大拉应力值为010388MPa,在011s时为010392MPa,在012s时为01425MPa。

4　结论

(1)动力开挖使既有隧道周围围岩发生应力重最小主应力分布规律是:在隧道的两个底角附近和隧道的右半圆面范围内会产生明显的压应力集中,最大压应力随时间的变化不大,在0103s时为216MPa,在011s时为2168MPa,在012s时为217MPa。31512　振速分析

分布。动荷载作用后,隧道破坏的可能性增加。(2)隧道截面危险部位依次为:迎爆侧直墙中部到起拱线部分、迎爆侧起拱线以上部分和迎爆侧底角、背爆侧。既有隧道最大振动点、最大拉应力区以及最大位移位置均出现在迎爆侧直墙的中上部。迎爆侧的拉应力区是由于爆破地震波的反射作用造成的,而背爆侧则主要受爆破地震波的绕射作用,因此,受影响较小。

此次模拟对既有隧道周围围岩X向振动速度进行了计算,得出了既有隧道周围质点振动速度的变化情况。输入的是正弦地震荷载,各计算点的振动速度随时间的变化也基本成一种简谐波的形式。实际工程中是以围岩质点的峰值振动速度作为判别洞室稳定性的根据。计算得出的各点峰值振动速度见表2。

表2　截面各点的振速峰值Vmax表(cm/s)

y=0

y=3

Vmax

y=6

Vmax

y=9

Vmax

y=12

Vmax

参考文献

[1]刘春玲,祁生文,童立强,等.利用FALC3D分析某边坡

地震稳定性[J].岩石力学与工程学报,2004,23(16):

2730-2733.

[2]易长平,卢文波.开挖爆破对临近隧道的震动影响研究

[J].工程爆破,2004,10(1):1-5.

[3]杨年华,刘　慧.近距离爆破引起的隧道周边振动场

[J].工程爆破,2000,6(2):6-10.

[4]毕继红,钟建辉.邻近隧道爆破震动对既有隧道影响的

坐标(-815,0,216)(-714,0,115)(-7,0,0)(-7,0,-1)(-7,0,-2)(-7,0,-3)

坐标(-815,3,216)

坐标坐标坐标

Vmax

51661771476415

918(-815,6,216)1118(-815,9,216)1215(-815,12,216)1315

(-714,3,115)1218(-714,6,115)1512(-714,9,115)1615(-714,12,115)1713(-7,3,0)(-7,3,-1)(-7,3,-2)(-7,3,-3)

1415141118811

(-7,6,0)(-7,6,-1)(-7,6,-2)(-7,6,-3)

171616141012

(-7,9,0)(-7,9,-1)(-7,9,-2)(-7,9,-3)

1812181511

(-7,12,0)(-7,12,-1)(-7,12,-2)(-7,12,-3)

191812151611

研究[J].工程爆破,2004,10(4):69-73.

[5]刘　波,韩彦辉.FLAC原理、实例与应用指南[M].北

京:人民交通出版社,2005.

[6]潘晓马,张成满,温向东,等.新建隧道施工对邻近既有

(-1115,0,216)-215(-11153216)(-1216,0,115)115(-13,0,0)(-13,0,-1)(-13,0,-2)(-13,0,-3)

11511251116

(-12163115)(-1330)(-133-1)(-133-2)(-133-3)

316(-1115,6,216)415(-1115,9,216)2211119116214

5(-1115,12,2112)513

3314219312315

(-1216,6,115)215(-1216,9,115)218(-1216,12,115)(-13,6,0)(-13,6,-1)(-13,6,-2)(-13,6,-3)

2182142219

(-13,9,0)(-13,9,-1)(-13,9,-2)(-13,9,-3)

313217213313

(-13,12,0)(-13,12,-1)(-13,12,-2)(-13,12,-3)

隧道安全性影响数值分析[J].铁道建筑技术,2002,

(1):29-31.

收稿日期:2010202204

　西部交通科技62

因篇幅问题不能全部显示，请点此查看更多更全内容

查看全文

全部栏目

爆破振动下小净距隧道围岩的力学行为分析