一种隧道洞口段抗震设防长度计算方法技术

本发明专利技术公开了一种隧道洞口段抗震设防长度计算方法，包括：A、对地基梁单元进行分析；B、引入隧道横向变形系数η1和轴向变形系数η2，得到修正后的地震波动形状位移方程；C、计算得到所述横向变形系数η1和轴向变形系数η2；D、根据步骤B与步骤C，得到对应的最大弯曲应力方程和最大轴向应力方程；E、得到修正后的最大弯曲应力方程和修正后的最大轴向应力方程；F、确定所述修正后的最大弯曲应力方程和所述修正后的最大轴向应力方程的各个参数数值；G、根据所述各个参数数值计算隧道最大抗震设防长度。本发明专利技术的一种隧道洞口段抗震设防长度计算方法，对隧道洞口段的设防长度提出了明确的解决方法，提高了现有技术的隧道洞口段的设防长度计算的完整度和精确度。

【技术实现步骤摘要】
一种隧道洞口段抗震设防长度计算方法
本专利技术涉及隧道抗震领域，具体涉及一种隧道洞口段抗震设防长度计算方法。
技术介绍
有关隧道震害特征研究起始于本世纪初，国内外学者通过现场调查，对隧道衬砌结构受损程度和震害原因进行了初步研究。如W.L.Wang等(W.L.Wanga，T.T.Wangb，f，J.J.Sua.AssessmentofdamageinmountaintunnelsduetotheTaiwanChi-ChiEarthquake.TunnellingandUndergroundSpaceTechnology，2011，16：133-150.)根据台湾集集地震中的57个山岭隧道受损程度，率先把受损隧道分为无破坏、轻度破坏和重度破坏三个等级，并给出了隧道衬砌结构震害可能涉及的原因。国内学者吉随旺等(吉随旺，唐永建，胡德贵等.四川省汶川地震灾区干线公路典型震害特征分析.岩石力学与工程学报，2009，28(6)2009：1250-1261)依据震害调查统计资料，研究认为地震惯性力和围岩失稳作用是导致隧道结构震害的两个主要原因。而高波等(高波，王峥峥，袁松.汶川地震公路隧道震害启示.西南交通大学学报，2013，44(3)：337-343)认为，山岭隧道震害多为面波、横波、地基失效和围岩作用所致。隧道洞口段作为抗震薄弱部位，在对地震区隧道进行抗震设计时，隧道洞口地段的隧道衬砌一般应予以加强，但隧道洞口段衬砌的设防长度到底多长，结构抗震才能达到最优，是隧道工程界普遍关心的问题。比如《铁路工程抗震设计规范》(GB50111-2006)也对隧道洞口段抗震设防长度做出了规定：设防长度可根据地形、地质及设防烈度确定，并不小于2.5倍结构洞跨，该设防长度主要根据一些宏观震害调查来确定。索然序等(朱长安，高波，索然绪.强震区隧道洞口段振动台模型试验研究.现代隧道技术，2009，45(1)：48-52)从弹性介质半无限空间的波动方程出发，建立的瑞利波下的运动方程，得到了瑞利波下隧道的位移方程，进而推导出结构的地层峰值位移和地层峰值加速度表达式。在此基础上考虑隧道埋深、地层条件和波长对隧道设防的影响；分别从地层曲率、地震烈度和位移三个方面对洞口设防长度的取值进行了探讨。结果表明：洞口段埋深超过150m，在9度地震作用下，地层对隧道的影响位移已很小，不需要再设防；再根据抗防埋深确定洞口段抗震设防长度；应该延长现行隧道设计规范规定的洞口段抗震设防长度。耿萍等(李林，何川，耿萍等.浅埋偏压洞口段隧道地震响应振动台模型试验研究.岩石力学与工程学报，2011，30(12)：2540-2548)以跨度6.4m、单线140km/h的铁路隧道为例，利用数值模型分析了衬砌物理力学参和数围岩等级条件等因素对隧道洞口段在地震作用下衬砌内力变化的影响，数值分析结果表明：距离洞口段超过3倍隧道洞跨后，衬砌内力变化明显减小；振动台试验结果也证明取隧道洞口段为3倍隧道洞跨的设防长度时，减震效果显著。而高峰等(高峰，石玉成，严松宏.隧道洞口段的抗震设防长度.中国公路学报，2006，1，(3)：65-70)、苏惠等(苏惠，贾良，严松宏等.隧道洞口段结构地震响应分析.水利与建筑工程学报，2010，8(2)：156-159)采用粘-弹性人工边界，运用Newmark隐式时间积分有限元单元法，对隧道三维地震反应分析来确定隧道洞口抗震设防长度；研究结果表明：隧道洞口段的围岩岩性主要关系着隧道抗震设防长度，洞口段的围岩岩性越差，隧道抗震设防长度值就越长；隧道衬砌横断面形式和洞口段周围临空面的存在与否与隧道的设防长度取值关系不大。周德培(周德培.强震区隧道洞口段的动力特性研究.地震工程与工程振动，1998，8(1)：124-130)根据南昆线上草庵隧道和乐善村2号隧道的试验结果，研究了强震区隧道洞口段的破坏形态和抗震设防长度。李育枢等(李育枢，李天斌.高烈度地震区山岭隧道洞口减震问题的数值模拟研究.公路交通科技，2009，26(10)：100-104)以国道318线上黄草坪2号隧道进口段为参考原型，针对隧道减震措施，开展大型振动台物理模型试验分析研究，研究结果表明：一般从距模型洞口120～150cm(对应原型48～60m)后，地震反应逐渐趋于平缓等等。隧道因其埋于地下，纵向可长可短且被围岩所包围，抗震设计比较复杂。如何保证地下隧道在未来强震中安全可靠地发挥应有的功能，是隧道工程界迫切亟待解决的问题。国内外学者已对地下隧道洞口段抗震分析进行了研究，其结果分析出隧道洞口段为抗震的薄弱区域，需要进行加固减震处理。即便对隧道洞口抗震设防的研究取得了很多成果，但是洞口段的设防长度仍然没有明确的结论。
技术实现思路
为了解决这些潜在问题，本专利技术的目的在于克服现有技术中所存在的上述不足，提供一种隧道洞口段抗震设防长度计算方法。为了实现上述专利技术目的，本专利技术采用的技术方案是：一种隧道洞口段抗震设防长度计算方法，包括：A、对地基梁单元进行分析，得到在地震作用下隧道的二维运动微分方程；B、引入隧道横向变形系数η1和轴向变形系数η2，得到修正后的地震波动形状位移方程；C、根据所述修正后的地震波动形状位移方程、所述地震作用下隧道的二维运动微分方程，计算得到所述横向变形系数η1和轴向变形系数η2；D、根据步骤B与步骤C，得到对应的最大弯曲应力方程和最大轴向应力方程；E、对所述对应的最大弯曲应力方程和所述最大轴向应力方程进行修正后的，得到修正后的最大弯曲应力方程和修正后的最大轴向应力方程；F、确定所述修正后的最大弯曲应力方程和所述修正后的最大轴向应力方程的各个参数数值；G、根据所述各个参数数值计算隧道最大抗震设防长度。进一步地，所述二维运动微分方程为：其中，E为梁的弹性模量，I为梁截面的惯性矩，v(x,t)为隧道横向位移，u(x,t)为隧道的轴向位移，ρ为梁密度，s为梁横截面面积，kl为地基梁横向变形系数，ka为地基梁轴向变形系数，gl(x,t)为地基土横向位移，ga(x,t)为地基土轴向位移。进一步地，所述修正后的地震波动形状位移方程为：其中，v(x,t)为隧道横向位移，u(x,t)为地基土的横向位移，w为地震波的圆频率，λ′为地震波的波长，A′为地震波的振幅，为地震波的入射方向与隧道的轴线夹角，t为时间，x为波的传播位移。进一步地，所述横向变形系数η1和轴向变形系数η2分别为：其中，E为梁的弹性模量，I为梁截面的惯性矩，ρ为梁密度，s为梁横截面面积，kl为地基梁横向变形系数，ka为地基梁纵向变形系数，λ′为地震波的波长，w为地震波的圆频率，为地震波的入射方向与隧道的轴线夹角。进一步地，所述对应的最大弯曲应力方程和最大轴向应力方程分别为：其中，σlmax为最大弯曲应力，σamax为最大轴向应力，E为梁的弹性模量，I为梁截面的惯性矩，kl为地基梁横向变形系数，ka为地基梁轴向变形系数，ρ为梁密度，s为梁横截面面积，H为基岩至地表面的距离，Sv为速度反应谱的值，λ′为地震波的波长，为地震波的入射方向与隧道的轴线夹角，h为隧道中心埋深，D为隧道的等效外径，Vs为土层的剪切波速，Ts为土层的固有振动周期，其中，Ts＝4H/Vs。进一步地，所述修正后的最大弯曲应力方程和修正后的最大轴向应力方程分别为本文档来自技高网...

【技术保护点】
一种隧道洞口段抗震设防长度计算方法，其特征在于，包括：A、对地基梁单元进行分析，得到在地震作用下隧道的二维运动微分方程；B、引入隧道横向变形系数η1和轴向变形系数η2，得到修正后的地震波动形状位移方程；C、根据所述修正后的地震波动形状位移方程、所述地震作用下隧道的二维运动微分方程，计算得到所述横向变形系数η1和轴向变形系数η2；D、根据步骤B与步骤C，得到对应的最大弯曲应力方程和最大轴向应力方程；E、对所述对应的最大弯曲应力方程和所述最大轴向应力方程进行修正后的，得到修正后的最大弯曲应力方程和修正后的最大轴向应力方程；F、确定所述修正后的最大弯曲应力方程和所述修正后的最大轴向应力方程的各个参数数值；G、根据所述各个参数数值计算隧道最大抗震设防长度。

【技术特征摘要】
1.一种隧道洞口段抗震设防长度计算方法，其特征在于，包括：A、对地基梁单元进行分析，得到在地震作用下隧道的二维运动微分方程；B、引入隧道横向变形系数η1和轴向变形系数η2，得到修正后的地震波动形状位移方程；C、根据所述修正后的地震波动形状位移方程、所述地震作用下隧道的二维运动微分方程，计算得到所述横向变形系数η1和轴向变形系数η2；D、根据步骤B与步骤C，得到对应的最大弯曲应力方程和最大轴向应力方程；E、对所述对应的最大弯曲应力方程和所述最大轴向应力方程进行修正后的，得到修正后的最大弯曲应力方程和修正后的最大轴向应力方程；F、确定所述修正后的最大弯曲应力方程和所述修正后的最大轴向应力方程的各个参数数值；G、根据所述各个参数数值计算隧道最大抗震设防长度。2.根据权利要求1所述的一种隧道洞口段抗震设防长度计算方法，其特征在于，所述二维运动微分方程为：其中，E为梁的弹性模量，I为梁截面的惯性矩，v(x,t)为隧道横向位移，u(x,t)为隧道的轴向位移，ρ为梁密度，s为梁横截面面积，kl为地基梁横向变形系数，ka为地基梁轴向变形系数，gl(x,t)为地基土横向位移，ga(x,t)为地基土轴向位移。3.根据权利要求1所述的一种隧道洞口段抗震设防长度计算方法，其特征在于，所述修正后的地震波动形状位移方程为：其中，v(x,t)为隧道横向位移，u(x,t)为地基土的横向位移，w为地震波的圆频率，λ′为地震波的波长，A′为地震波的振幅，为地震波的入射方向与隧道的轴线夹角，t为时间，x为波的传播位移。4.根...

【专利技术属性】
技术研发人员：张军伟，杨婉，赵成林，马小杰，梅志荣，钱进，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：四川;51