不同埋置深度山岭隧道地震响应分析.docVIP

不同埋置深度的山岭隧道地震响应分析 蒋树屏 方 林 林 志 （重庆交通科研设计院重庆400067） 摘要：根据大量震害资料的调查分析发现，地下结构埋埋置深度对其地震破坏程度影响很大。通过有限元方法计算八种不同埋置深度条件下的山岭隧道地震响应，并对计算模型的地震输入方法进行了验证，证明地震波输入处理方式的合理性。讨论埋深对结构动力响应的影响，提取衬砌关键节点的竖直向和水平向加速度、位移峰值，分析随着埋深的增大，加速度和位移峰值的变化情况；并以拱顶为例，计算每个埋深变化段的加速度和位移峰值变化率，得到了一定的规律性。此外，还分析了埋深增大对衬砌结构内力峰值的影响。最后提出在高烈度地震区修建隧道时其埋置深度尽可能不小于50m，这为相关工程的修建提供了参考依据。 关键词：埋置深度；山岭隧道；地震响应 0 引言 近些年，地下结构的抗震研究越来越被重视，大量的震害资料也反映出，地下结构在地震中会产生不同程度的破坏。对于山岭隧道，在几次典型的地震中均有遭受严重破坏的实例[1],[2],[3]。特别是1995年发生的日本阪神大地震，阪神地震震级7.3级，灾区内有10%的山岭隧道受到了破坏，改变了以往山岭隧道不会遭受地震破坏的传统观点。 埋置深度是隧道地震响应的重要影响因素之一，在高烈度地震区修建隧道，选择合适的埋深，既能将可能发生的地震对隧道的影响降到最小，又能在此基础上满足经济性的要求。Sunil Sharma等对132例地下工程震害资料进行统计后提出，深度大于50m时破坏程度明显减小，在300m以下没有严重的破坏[4]。李天斌对5.12汶川大地震中遭受震害的隧道进行分析后认为，在硬质岩隧道埋深大于50m时震害程度为中等～轻微，埋深大于100m后隧道几乎没有震害或震害轻微[5]。2008年5月12日后，作者对汶川地震震灾进行了调查。调查表明：此次地震对灾区的房屋、公路、桥梁和隧道等工程结构破坏严重，甚至倒塌（如照片1所示），隧道出现衬砌开裂、错位、剥落、垮塌，山体发生滑动、滚石等，尤其在隧道洞口段或断裂带地层较为严重[]。穿越断层隧道振动台模型试验埋置深度隧道地震响应试验研究[]。 照片1 地震对公路的损情况(2008.5.16) 照片振动台模型试验 在调查的基础上，通过有限元数值计算，分析了在岩质地层条件下不同埋置深度山岭隧道的地震响应规律及内力变化特征，为今后实际工程应用提供依据。 1 计算模型及工况 1.1 计算模型 本文计算采用的有限元模型均为二维模型，横向宽度为160m，拱顶以上部分为隧道埋置深度D，拱顶至下边界的距离取60m。为简化计算模型，便于进行结果分析，均假设隧道位于均质岩层中。模型单元划分长度按2m控制，其中围岩采用四节点平面单元，衬砌采用梁单元模拟。有限元模型如图1所示。 图1 计算模型 隧道采用公路三车道标准断面，内轮廓净宽为15.17m，净高9.4m，衬砌厚度50cm，衬砌材料为C30钢筋混凝土。 阻尼耗散能量，使自由振动的振幅随时间衰减，限制其频率与固有频率一致的加载引起振动的幅值。本文计算均采用瑞利阻尼，阻尼系数通过前两阶自振频率计算得到。 1.2 模型边界及输入方式 通过很多的计算分析可知[], [9]，粘弹性人工边界相对于粘性边界、固定边界等具有更高的精度，因此本文计算采用粘弹性人工边界，通过弹簧、阻尼器等实现，弹簧和阻尼器按下式进行赋值。 ， （1） ， （2） 式中：KBN 、KBN——分别为法向与切向弹簧刚度； CBT 、CBT——分别为法向与切向阻尼器的阻尼系数； R——波源至人工边界点的距离； cp、cs——分别为P波和S波波速； G——剪切模量； ρ——介质密度。 实现地震波动合理输入是决定波动模拟是否成功的关键，直接影响到计算结果的精度及可信度。本文采用与一致粘弹性人工边界相对应的地震动输入方法，即在人工边界上施加如下等效荷载： （3） 式中：FB——施加的等效荷载； KB 、CB——分别为弹簧刚度和阻尼系数；——已知位移场； τ0——在原连续介质中由位移产生的应力。 1.3 计算工况及参数 整个计算过程按隧道不同的埋置深度共分为8个工况，分别对应埋深（D）为5m、20m、50m、100m、200m、350m、500m和1000m。围岩条件按级围岩考虑，各计算参数取值如表1所示。 表1 材料力学参数 参数 弹性模量,GPa 重度,kN/m3 泊松比 围岩 10.0 23.0 0.28 衬砌 28.0 25.0 0.2 1.4 地震波 在模型底部同时输入水平和竖直向地震波，地震波采用1999年台湾集集地震记录，持时为30s。其中，水平向地震峰值加速度为-4.94m/s2，对应时刻10.38s；竖直向地震峰值加速度为5.20m/s2，对应时刻9.28s。地震波波形如图2