爆破振动作用下小净距隧道围岩力学行为分析.docVIP

爆破振动作用下小净距隧道围岩力学行为分析 　　摘要:本文运用数值分析计算软件,对小净距隧道中既有隧道受邻近隧道爆破震动的动力响应特性进行了分析,得出了爆破地震波作用下,隧道周边围岩应力、位移和振动速度的分布规律,为小净距隧道的设计与安全施工提供了科学依据。 　　关键词:小净距;隧道;爆破震动;数值模拟 　　 　　爆破震动作用下邻近洞室的响应是一个复杂的动力学问题。数学解析法虽然能给出准确可靠的分析结果,但数学处理上的困难制约了它的适用范围。目前,理论分析仍只能解决介质、结构线弹性,入射波为稳态波,边界条件较简单的情况。要分析更复杂的实际问题,则需借助实验或数值计算方法。计算机数值模拟技术的发展为这一问题的研究提供了很好的工具,很多学者已经作过这方面的研究。 　　本文结合某高速公路小净距隧道开挖施工,应用数值模拟软件模拟开挖过程中隧道爆破施工对临近既有隧道的影响,得到隧道围岩应力场、位移场及爆破振动速度场的分布规律。 　　1爆破振动作用下影响邻近隧道稳定性的因素 　　(1) 根据动力弹性理论,爆破时,岩体中产生的应力与质点的振动速度成正比。在爆破地震波引起相同振速的情况下,低频波比高频波对洞室围岩的影响要大。目前,许多国家在制定或修订爆破震动安全判据时,普遍考虑了振速和频率两个指标的共同影响。 　　(2) 地质条件因素。地质条件主要指相临隧道之间的地质材料与地质构造,也是地震波的传播介质。研究表明,在相同爆破地震波的作用下,岩石越硬,则围岩中引起的峰值振动速度越小。从波动力学的角度讲,爆破震动对邻近洞室的影响是地震波在传播过程中与地质介质相互作用的结果。应力波参数是两者波阻抗的函数,不同的岩石对爆破地震波的阻尼效应不同,特别是当洞室围岩存在某些特殊地质构造(如断层软弱夹层,溶洞等)时,应力波的传播更受影响。 　　(3) 隧道自身条件因素。隧道的横断面形状、被覆材料,衬砌结构、爆源的相对位置、洞径与爆破应力波长的相对尺寸等直接影响洞壁周边的应力分布,而洞径是影响围岩应力和振速的重要因素。在其它条件相同的情况下,随着隧道直径的增大,迎爆一侧主要由于应力波的反射拉伸作用,其动力集中因子增大,而背爆一侧没有明显的变化。长江水利水电研究院建议的允许爆破震动速度和破坏标准见表1。表1中的数据为重复爆破的速度值,单次爆破时,可适当加大,但不超过50%。 　　2 工程概况 　　某高速公路隧道进口端前方铁路隧道下穿,出口端基本与铁路平行。隧道K2+880~K3+250,YK2+880~YK3+260为小净距,线间距4~12.7m,其余为分离式,隧道洞轴线走向约196°,净空10.25×5.0m,其中左线隧道全长2092m,纵坡坡度为-2.8%,隧道底高程为1473.36~1686.11m;右线隧道底高程为1473.36~1684.93m,全长2124m,纵坡坡度为-2.8%。隧道最大埋深231.45m。 　　3爆破对邻近隧道影响的数值模拟分析 　　3.1数值计算模型 　　计算模型由两个截面形状为马蹄形的隧道及其周围围岩组成。分析地下开挖时,模型边界应距离开挖洞室中心大约10倍直径的长度。隧道断面上半部分为半径3 m的半圆,下半部分为3 m×10m的矩形。两硐中心相距20 m,计算模型x方向范围为-70～70 m, Y方向范围为0～30m, Z方向范围为30～-60。用FLAC3D进行动力计算,网格剖分尺寸受输入波动的最短波长控制。设网格最大尺寸为Δl, 输入波动的最短波长为λ,则Δl必须小于(1/10～1/8) 。另外,动力计算要求网格划分的尺寸要相对均匀,网格最大尺寸和最小尺寸的比例一般不超过5。因此,划分模型为44160个单元,共50094个节点。 　　3.2岩体物理力学参数及力学模型 　　不考虑岩体物理力学参数变化对计算的影响,计算模型中的岩体为单一岩石材料。参数取值为体积模量K = 7.73 GPa,剪切模量G = 5.32 GPa,抗拉强度σt=0.2459MPa,粘聚力C=0.725MPa,体积密度ρ=2720kg.m-3,内摩擦角φ=70.77°。模型采用理想的弹塑性模型,屈服准则采用Mohr-Coulomb强度准则。 　　3.3初始条件与边界条件 　　(1) 初始地应力场。初始地应力场仅按自重应力场考虑。根据弹性力学原理,竖向应力σv =γH,水平应力σh = kσv ,式中γ为岩体容重;H为埋置深度;k为侧压力系数,k =μ/(1-μ);μ为泊松比。 　　(2) 静力计算边界条件。采用位移边界条件,即模型的左右(X方向)边界,前后(Y方向)边界和底边界均施加位移约束条件,上边界取自由边界。 　　(3) 动力计算边界条件。FLAC3D求解动力问题的边界条件设置有远置人工边界条件和粘滞边界条