关键词:桥梁工程;隧道锚;极限承载力;Mindlin应力解;动荷载
中图分类号:U448.25 文献标志码:
悬索桥在面对复杂山谷河流地形时具有显著优势,而作为承载的关键构件之一的隧道锚因其较高的性价比与环境扰动较小等特点,受到了众多设计者的青睐[1-2].隧道锚的结构形式呈现为上窄下宽的锥形体,这种独特的结构形式给隧道锚带来巨大承载力的同时也决定了其承载特性复杂.尽管近年来诸多学者在静力条件下进行了大量的现场试验[3-6]与室内试验[7-8],对隧道锚从工作至破坏的过程进行了充分的论述,对其承载力进行了一系列的公式推导[9-13],实际工程中,主缆缆力受到桥面行车、风雨等荷载的影响,无时无刻不在变化,这种变化传递至锚碇后的影响是不可忽略的,因此对于隧道锚动力荷载下的研究具有一定的现实意义.
因为夹持效应的存在,悬索桥隧道锚动力研究的核心问题与静力类似,在于明确锚岩的相互作用过程,重点揭示动荷载作用下接触面应力的分布及变化过程,掌握接触面岩体状态变化情况.刘新荣等[14]通过运用数字图像技术,研究了隧道锚承载过程中的岩体裂纹萌发开展过程,对隧道锚的渐进破坏全过程进行了详细的阐述.王东英等[15]则通过有限差分软件,借助边坡滑面搜索方法,对锚岩破坏面的产生进行了追踪.李明等[16]模拟了地震作用下锚塞体对于主缆动张拉荷载的响应.颜冠峰等[17]对隧道锚在地震荷载下各处的大小主应力以及隧道锚与隧道的力学响应异同进行了研究分析.抗拔桩作为结构形式和受力与隧道锚相近的构件,其理论与研究方法具有重要的借鉴意义,范书立等[18]建立了循环荷载下斜向抗拔桩的计算模型,并最终确定了循环荷载作用下的承载力.乔东生等[19]则通过建立拖曳锚有限元模型,并引入循环强度模型概念,揭示了循环次数对极限抗拔力的影响规律.
现阶段隧道锚仍具有巨大的承载潜力待挖掘,能实现在满足工程要求的条件下尽可能地降低成本这一目的.以往研究已经表明,影响隧道锚承载力的因素一方面在于桥址的地质条件,这部分因素在选定地点后就无法改变,另一方面则在于隧道锚本身的几何参数[20-22],这些参数是可以择优选取的,因此根据不同因素的影响程度来选取适合参数是降低造价的合理手段.但目前众多学者对于锚塞体参数的研究集中在静力破坏阶段,且没有在统一标准的情况下衡量其影响程度.
以上研究分析表明,现阶段对隧道锚在动力荷载作用下破坏模式研究较少.因此本文将动荷载的施加作为研究的基础,以锚岩接触面应力的产生与变化作为研究的出发点,探究动荷载作用下,与静载相比隧道锚位移等响应的异同点;提取并分析承载后各阶段的状态变化,得到隧道锚动荷载作用下的工作过程与破坏形态;在此基础上,进一步建立各种工况来分析锚塞体长度、间距、扩展角等几何因素以及动荷载参数对承载能力的影响.
1 隧道锚动荷载作用下的应力解答
隧道锚所受载荷主要从主缆传递而来,而主缆力又由静载和动载组成,主缆静载可直接采用设计荷载代表,因此将静力荷载P0取为101202 kN,并在此基础上逐步变化.但主缆所传递的动荷载值近年很少有人研究,要想合理准确地获得隧道锚在动载作用下的响应情况,就必须选取适合的动荷载时程曲线.目前对于隧道锚所受主缆动荷载并无明确表示方法,主要存在以下两种人工合成的方式,分别为正弦函数振动形式与随机振动分布形式[16].本文主要采用第一条动荷载人工合成方式,具体动荷载求解推导过程如下.
悬索桥正常使用过程中,主缆传递至锚碇的荷载是变化的,在受到风雨荷载、桥上行车等荷载作用时,主缆缆力会产生随时间变化并在一定均值范围内的波动.因此隧道锚接触面应力不会仅仅包括上述两种,本文在前人研究基础上,对隧道锚应力分布情况进行进一步延伸,考虑到锚塞体所受动荷载影响,对其动荷载形式以及动荷载作用下的应力分布进行理论与模型研究.主缆缆力是随时间变化的波动函数,但由于悬索桥自重及二期铺装等恒载作用,且在一定时间内,桥上车辆所带来的静载都是定值,因此缆力波动的中值应与主缆缆力的静载部分保持一致,波动部分则是随时间变化的动荷载.因此本文采用正弦函数形式模拟缆力波动,由此可知动荷载不是一定值,如式(1)所示.
动荷载对于隧道锚承载能力的影响表现为两个方面,首先,动荷载的施加,使锚岩接触部位的岩体在卸荷过程中存在残余应力,在长时间经历循环荷载的过程中,岩体中的残余应力不断叠加至后续的循环过程响应中,其残余应力与应变不断累积,使得其容许强度降低,最终导致其承载能力的下降;其二,在动荷载的作用过程中,虽然动荷载作用是长时间不断变化的,但由于动荷载幅值仅约为静载幅值的15%,因此动荷载对于岩体性质影响不大,对承载力的影响是通过在原有应力基础上增加一部分动应力,而使得在较低静载作用下即达到岩体的容许强度值,进而导致隧道锚的承载失效.以上两种情况可以看作动荷载对于隧道锚承载能力影响的两种极端情况,实际承载过程中,隧道锚受动荷载影响应介于两种情况之间.因此通过有限差分软件建立模型,提取隧道锚在动荷载作用下锚岩接触面的应力等响应分布,以验证上述情况的合理性.
2 动荷载作用下锚岩系统承载特性分析
本模型依托云南普立特大桥隧道锚工程,根据表1普立特大桥隧道锚锚区岩体主要物理力学参数表,将模型中岩体看作理想弹塑性材料,进行参数等一化处理,主要反映主体岩体的性质,所得数值模型的计算参数见表2.图2(a)为所建模型图.
采用有限差分软件建立锚岩数值模型时,周围岩体采用莫尔-库仑模型,锚塞体采用弹性模型,模型计算中需设置边界条件和接触面模型,故应力边界条件采用四周自由场边界条件,底部为静态边界,实际工程中锚塞体与岩体之间的胶结面可以通过有限差分软件中的interface单元进行模拟,如图3(a)所示,为接触面的力学模型图.图中S为滑块,T 为抗拉刚度,Lm与Ln分别为与相应节点的连接长度,ks与kn分别为接触面的剪切刚度与法向刚度,可由式(7)计算得到.
图3(b)为接触面数值模型图,表3为接触面模型具体参数值,模型因锚塞体产生的附加应力仅受到荷载本身影响,因此本节仅探究动荷载对于隧道锚的影响不考虑其他因素来建立模型,其余按照普立特大桥边坡实际工程地质条件及锚岩参数设置.
2.1 动荷载作用下锚岩接触面响应分析
缆索在外部荷载影响下会产生各个方向的振动,但可以分解到x,y,z三个坐标方向进行研究,如图2(b)所示,z 方向对应主缆拉力方向,x 方向对应垂直于主缆指向侧壁方向,y方向对应垂直于主缆指向拱顶方向,其中z 方向为拉力方向,即为主要振动方向,其余两方向动荷载幅值会进行相应的折减,加载方式采用静载与动荷载分别施加的方式,即在静载基础上施加动载,由此可以看出动荷载施加后的应力增量情况.
图 4为动静荷载作用下锚岩接触面应力分布情况.在施加低幅值的荷载时,动荷载施加后的拱顶接触面处应力分布与仅施加静载相比较小,直到荷载幅值施加至7P0时,动静荷载应力分布开始出现差异,距前锚面25 m处,动荷载施加后的应力值开始大于静载时该处的应力值,通过提取隧道锚拱顶围岩塑性区发展情况,可知部分围岩此时已经进入了塑性状态,说明在动荷载作用下产生了较为明显的残余应力.继续增大荷载施加幅值,动荷载作用下产生的残余应力也随之增大,并且可以看出应力的差异由锚塞体5/6处逐渐向前锚面处扩大,说明此时的塑性区也由后锚面附近逐渐向前锚面处延伸.
通过接触面应力分布对比图,可以得到锚岩接触面的塑性区分布及开展情况,但隧道锚整体变形情况无法得到直观的体现,因此采用有限差分软件切片功能,提取出锚岩接触部位在各幅值作用下的位移极大值点,绘制动静荷载作用下位移对比图,如图5所示.图4和图5中所表示的动荷载其实表示的是考虑汽车和风荷载与静力荷载(P0)相复合的一种荷载,具体表示为复合荷载(图中动载)=静力荷载P0+动荷载0.15P0,这种复合荷载其实是静力荷载与推导的动力荷载叠加的一种荷载,为方便表述称该荷载为动载.但在大跨度悬索桥的实际工程中的运营阶段,隧道锚的受力不只是会受到缆力的静力作用,更多的是受汽车及风荷载共同作用下与缆力叠加的作用,故此图中动载表示的是一种复合作用.由于隧道锚自身结构特性及相关规范要求,需保证隧道锚在各种极端荷载作用下留有足够的安全冗余度,因此对于其计算范围选择较为广泛.根据图4和图5的分析,动载曲线在荷载幅值达到一定程度后开始偏离静载曲线.在低幅值荷载作用下,动荷载后的拱顶接触面应力分布较静载时小,导致两者位移差异不明显.然而,当荷载幅值增至7P0时,动静荷载的应力位移分布开始显现差异.在距离前锚面25 m处,动荷载后的应力值超过静载时的值,引起动荷载位移与静载位移的逐渐偏离.分析隧道锚拱顶围岩的塑性区发展情况显示,部分围岩已进入塑性状态,表明动荷载作用下产生了显著的残余应力.在7P0幅值下,动荷载的应力位移与静载作用下的应力位移出现差距,塑性区也从后锚面附近逐渐向前锚面处扩展.这一现象解释了静载与动载之间显著差异的产生,而这一差异在7P0幅值处尤为明显.其差异出现值在7P0幅值处.
因为荷载沿z 轴施加,因此z 方向位移要显著大于x 和y方向,低荷载幅值阶段,静载位移与荷载基本保持线性关系,直至11P0左右,曲线出现拐点,并开始出现非线性发展趋势,说明此时已经进入塑性状态;相较于静载作用,动荷载位移曲线在7P0处出现明显拐点,7P0之后位移随荷载幅值的增加而迅速增加,在11P0时,动荷载z 方向位移相比于静载时增加了36%,y 方向位移则增加了112%,说明隧道锚在y 方向上的动荷载敏感程度更大.
位移沿拱顶分布情况如图6所示,在各级荷载幅值作用下,拱顶接触面位移存在较大差异,为了更好地分析位移沿拱顶分布的规律,对接触面的竖向位移w(x) 与距前锚面的距离x 进行了归一化处理,得到w(x)/wm 与x/L 的关系曲线,其中,wm 为拱顶处竖向位移最大值,L 为锚塞体轴向长度.
由图6可知,在各荷载幅值作用下,拱顶处竖向位移分布特征基本一致,整体位移呈现三段式分布,从前锚面至后锚面位移幅度先增加后减小,第一个拐点出现在锚塞体1/5处,峰值点则出现在锚塞体9/10处.动荷载作用下位移分布和应力分布形式上存在一致性,且与静载作用下的位移分布基本相同.
2.2 锚岩接触面阶段状态分析
上文分析得到锚岩接触面应力及位移分布情况,可以明显看出,7P0幅值动荷载作用下,动静荷载响应出现较大差别,说明锚岩接触部位开始进入塑性状态,但无法判断这一状态产生的原因,如果在7P0静载作用下塑性区就已经产生,那么动荷载的作用仅仅就是加快残余应力的积累,这种改变是由于“荷载”本身引起的;另一种原因则是由于“动力”这一因素引起隧道锚提前进入塑性状态,因此要想准确了解动荷载对于隧道锚的影响方式,就要对隧道锚每一阶段的状态进行分析.
由莫尔-库仑强度准则可以知道,在岩土体内部某一面上的剪应力超过其本身最大抗剪强度时,岩土体就会发生剪切破坏,同时会产生较大的塑性变形,这一种岩土体状态的改变在有限差分软件中可以采用最大剪应变增量来衡量,增量较大且连成条带状区域的可以视为进入塑性状态.提取锚岩接触面处在各级荷载作用下的剪应变增量最大值点,得到动静荷载作用下荷载幅值与最大剪应变增量的曲线,如图7所示.
图7中,随着荷载幅值的增加,锚岩接触面的最大剪应变增量也是随之增加的,但较低荷载幅值(约1P0至7P0)时期,最大剪应变增量基本不变,且动荷载与静荷载作用效果一致,此外通过图4所示在7P0荷载作用之前并无残余应力产生这一现象,也可以充分说明在此之前,无论静载还是动荷载作用下,隧道锚锚岩接触面一直处于弹性阶段;当荷载幅值大于7P0之后,静载作用下,增量变化仍较小,但动荷载作用下的应变增量出现了较大变化,8P0荷载较7P0荷载,增量增加达到了156%,同时通过剪应变增量云图可得,增量大值区域发生了转变,由后锚室围岩处转移到了后锚面锚岩接触面处,说明此时这一区域率先出现了塑性,而接触面其他部位尚未出现.
图8为同一荷载幅值下动静荷载作用下的塑性区发展程度对比图.静载作用下,荷载幅值达到11P0时,右侧接触面才开始出现塑性区,而动荷载作用下7P0时就已经产生塑性,两种情况下荷载幅值相差4P0,而本文施加设计动荷载幅值约为静载幅值的0.1 倍,因此此时动荷载幅值仅为0.70P0,远达不到4P0这一限值.同时,从刚产生塑性区的发展程度可得,静载作用下,第一次出现塑性区的范围分别为距后锚面0.30 m和4.30 m,而动荷载作用下第一次出现塑性区则已经达到距后锚面8.20 m和8.50 m处,说明即便在静载基础上增加了动荷载,整体荷载的幅值相较于静载有所增加;塑性区贯通整个接触面时,静载的荷载幅值为20P0,而动荷载仅为14P0,此时动荷载幅值为2P0.
以上结果表明,在动荷载作用下,隧道锚围岩的容许承载力已经发生了改变.由于围岩内部本身存在着一些初始裂纹与空隙,而这些裂纹和空隙存在着在动荷载作用下发生屈服破坏的“门槛值”[25],在动荷载的循环往复作用下,围岩内部部分较低“门槛值”的初始裂纹与空隙率先发生屈服,而在相同幅值静载作用下尚未达到岩体的抗压强度,因此在动荷载作用下塑性区的产生和发展较静载会更加明显,进而呈现出承载能力降低这一现象.
3 动荷载承载力影响因素分析
上述内容阐述了动静载作用下锚岩响应的异同点,为进一步研究动荷载对隧道锚的影响程度,从锚碇几何特征与动荷载参数两方面分别进行探讨.
首先动荷载作用下隧道锚的几何参数对承载力的影响,从锚塞体扩展角、锚塞体长度以及锚塞体间距三个方面共6个工况进行建模分析,如表4所示.
3.1 锚塞体扩展角对承载力的影响
计算得到不同扩展角影响下静力与动力加载对比曲线,如图9 所示. 当锚塞体扩展角α 分别为1.80°、2.86°和4.00°时,静力加载下的极限承载力分别为1 720.50 MN、1 922.80 MN和2 024 MN,动力加载下的极限承载力分别为1 365 MN、1 516.80 MN 和1 567 MN.动荷载作用下的极限承载力与静载作用下类似,会随着锚塞体扩展角的增加而增加.对比动静加载下的极限承载力可以发现,与静载极限承载力相比,动荷载承载力分别降低了20.66%、21.12%、22.58%,平均降低了21.45%,说明由于动荷载的循环往复作用,锚岩接触面岩体自身容许承载力降低,且在锚塞体扩展角较小的时候降低幅度更大,达到了20.66%.
从图10可得,动力加载下,锚岩接触面位移随着荷载的增加是不断增大的,在荷载一定的情况下,接触面位移随着扩展角的增加而减小,说明承载力随着扩展角的增加是增大的.随着动荷载的增大,z方向位移增加速率随之增大,x 方向虽然整体位移小于z 方向,但速度增加比率更大.
3.2 锚塞体长度对承载力的影响
选择锚塞体长度分别为25 m、35 m、45 m三种工况,进行长度对动荷载极限承载力影响的研究.如图11所示.三种不同工况下,静力加载下的极限承载力分别为1518 MN、1922.80 MN、2327 MN,动力加载下的极限承载力分别为1315 MN、1516.80 MN、2226 MN.随着锚塞体长度的增加,无论是静力加载还是动力加载情况下的隧道锚极限承载力都是增加的,静力加载下,锚塞体承载力与长度的关系曲线近似线性,而动力加载下,承载力增加的幅度更大,最大幅度出现在长度由35 m增加至45 m这一阶段,承载力增加幅度达到46.76%.同时也可得,动荷载承载力相比静载有所降低,三种工况下分别降低了13.37%、21.12%和4.30%.
由图12可知,动力加载下,与扩展角类似,锚岩接触面位移随着荷载的增加而增加,同时在荷载一定的情况下,承载力随着锚塞体长度的增加而增大.随着荷载的增加,z 方向位移增加速率逐渐减缓.
3.3 锚塞体间距对承载力的影响
图13为不同间距下动力与静力加载对比曲线,共三种不同工况,锚塞体间距分别为16 m、21 m、26 m,计算得到静力加载极限承载力分别为1 904 MN、1 922.80 MN、1 973 MN,动力加载条件下分别为1 500 MN、1 516.80 MN、1 619 MN.改变锚塞体间距对于承载力的影响较小,静力和动力平均承载力增加分别为1%和4%,这是因为相比于锚塞体本身桩径,两锚塞体之间间距过大,由中部围岩传递的应力较小,因此群桩效应并不明显.同时可得,动静加载下极限承载力差距依旧明显,三种工况下分别降低了21.22%、21.12%和17.94%.
从图14可知,与前两种因素类似,动力加载下锚岩接触面位移随着荷载的增加而增加,在荷载一定的情况下,承载力随着锚塞体长度的增加而增大.三种工况下z 方向位移差距不大,但x 方向位移差异明显,这是因为锚塞体中部围岩的应力是水平传递的,因此间距对于水平方向应力影响较大.
3.4 几何参数归一化处理
如图15所示,为了更好地分析各因素对隧道锚承载能力的影响程度,对各工况下锚塞体的几何参数Ii 和围岩位移wi 进行了归一化处理,经过线性拟合后得到了Ii /I0 与wi /w0的关系曲线,其中,Ii 为各工况变化参数值,i = 1,2,3;I0 为初始工况各参数值.wi为各工况下围岩位移值,i = 1,2,3;w0 为初始工况围岩位移值.各几何参数对于隧道锚承载能力影响存在差异,拟合后得到三条曲线的斜率值分别为k1 =-0.76,k2 = -0.89,k3 = -0.44,说明锚塞体长度对承载能力的影响最大,其次为锚塞体扩展角,锚塞体间距对承载能力的影响最小.
3.5 动荷载参数对承载力的影响
在考虑动荷载参数对承载力的影响之前,首先要明确有哪些动荷载参数,本文将由主缆传递的动荷载简化为正弦荷载,如式(5)所示,因此主要考虑的动荷载参数为主缆振动频率f,主缆振动频率会影响单位时间内锚岩接触面受到的循环加载次数,从而影响到围岩的容许承载力,因此研究主缆振动频率的变化影响是有必要的.
如图16(a)所示,首先通过设置边界与模型本构参数,使其在仅受重力条件下自由振动,得到锚碇自振频率约为9 Hz,之后设置频率增加梯度,计算并分别提取主缆振动频率为1~11 Hz下,振动幅值为7P0(从图5 中可以看出使围岩产生损伤破坏的“门槛值”约为7P0)的锚碇振动响应,得到不同频率下隧道锚各响应随频率增加的变化曲线.
从图16(b)可得,随着动荷载振动频率的增大,围岩的最大剪应变增量及最大位移也随之增大,且基本呈线性变化,但总体来说变化量很小.这是因为围岩内部的初始裂纹与缝隙,在动荷载的循环往复作用下,初始裂缝会逐渐闭合,产生塑性的变形,且这一过程往往与加载幅值和加载时间有关,本文设置加载时间为4 s,因此在幅值一定的情况下,频率越大,产生的塑性变形就越多.同时可得,在频率达到9 Hz左右时,三种响应都发生了不同程度的突升,这说明在荷载频率接近锚碇自有频率时,会出现“共振现象”,在此频率下,隧道锚的各响应会有一定突变.
岩体受到动荷载循环作用下的损伤破坏一般分为三个阶段[26],从图16(b)中的响应随振动次数(频率与设置时间的乘积)基本呈现线性发展的趋势可以看出,隧道锚围岩基本仍处于损伤破坏的初始阶段,随着动荷载的持续施加,围岩将会在某一时间发生迅速且大规模的破坏.在实际工程中,虽然主缆引起的振动荷载基本不会达到模拟试验的频率值,但其振动时间会一直持续,因此当振动荷载由于外界因素超过“门槛值”时,振动产生的损伤破坏就会随时间积累下来,此时振动频率增大就会因持续的振动而产生较大影响,所以实际工程中防止发生围岩疲劳损伤破坏的关键在于防止动荷载超过“门槛值”,即尽可能防止超载、偏载等事件的发生[27],其次需要采取适当措施减小锚碇所受动荷载的频率值.
4 动荷载承载力分析
上文第2.2节和第3节中已经简单分析了动荷载作用对承载力的影响,为了定量分析隧道锚承载力的降低程度,在本研究中定义承载力的降低系数η为:
式中:PD 为动荷载施加后达到与P0 相同状态时的静载幅值,隧道锚所处的状态可以通过位移、应力以及塑性区分布来划分,以塑性区分布的面积作为所处状态的衡量指标,但塑性区面积通常难以准确测量,而塑性区沿垂直于锚碇轴线方向的发展随荷载增加基本不变,因此以塑性区高度代替面积作为隧道锚状态的衡量指标;k为动荷载幅值与设计荷载的比值.
降低系数η为考虑到荷载幅值影响后,动荷载作用下,隧道锚的承载力降低情况,其变化曲线如图17所示,锚塞体左右两侧承载力降低系数曲线起点并不在后锚面,起始塑性区高度分别为8.50m和8.25 m,这是因为在此之前的动力加载并未产生塑性区,当动荷载分别达到7P0和9P0时(由图8可以得到),左右两侧围岩状态开始发生转变,并在短时间内延伸至一定高度.随着动荷载幅值的增加,承载力降低系数持续增加,岩体塑性区进一步发展,直至塑性区高度达到11.30m(锚塞体1/3处),此时左右两侧承载力降低系数均达到最大值,分别为0.31 和0.43,在此之后,降低系数开始显著下降,这说明从后锚面至锚塞体1/3处的围岩所受动荷载影响更大,与隧道锚夹持作用下所产生的应力分布规律是一致的.当塑性区贯穿整个接触面时,降低系数稳定在0.22附近,左右两侧承载力降低系数平均值分别为0.28和0.26.
5 结论
本文通过理论分析和有限差分软件相结合,量化研究了动荷载对锚岩接触面位移、应力、应变增量和塑性区的影响,及其对隧道锚围岩的破坏过程,以及各参数对隧道锚承载力的影响程度,并从塑性区发展状态出发,得到了动荷载承载力降低系数η,主要结论如下:
1)锚岩接触面应力是静应力和动应力的叠加,动荷载的循环作用促进了残余应力的积累,并随荷载幅值的增加而增加,其与静动应力差异界限值为7P0.
2)动荷载作用下,锚岩接触面应力位移等响应沿锚碇轴向分布形式与静力类似,围岩夹持效应仍存在,动应力分布峰值点同样位于距后锚面1/3总长的范围内,但各响应幅值大小存在差异,随着荷载幅值的增加,静动载的差异随之增加.
3)动荷载承载力的几何参数影响程度排序为:锚塞体长度gt;锚塞体扩展角gt;锚塞体间距,且z 方向总体响应较大;频率是动荷载影响最大的参数,相同加载时间条件下,随着频率增加,围岩最大位移、剪应变最大增量和锚碇最大剪应力都线性增长.在频率达到锚碇自振频率时,出现“共振现象”.
4)静载是围岩产生塑性破坏的主因,动荷载“门槛值”为7倍静力荷载,超过此值会影响围岩承载能力,加速塑性区发展,降低隧道锚极限承载力,相较于静载,动荷载作用下承载力的平均降低幅度约为21%.
下一步工作将通过缩尺试验对动荷载作用下隧道锚承载性能和其影响因素进行对比分析,为优化隧道锚极限承载能力提供参考.