列车振动荷载作用下地层动力响应分析

1、列车振动荷载作用下地层动力响应分析高玄涛（中铁第一勘察设计院集团有限公司，西安 710043）摘要：文章以某隧道下穿高速铁路项目为背景，隧道下穿高速铁路施工为特级风险源，对地层沉降要求极为严格。对隧道下穿施工期间，在列车振动荷载作用下产生的地层动力响应进行了三维数值模拟分析。数值模拟计算结果表明，下穿隧道与高铁路基交叉点处为该工程薄弱环节，在设计与施工过程中需对此处采取加强措施，以保证施工安全及不影响列车正常运行。关键词：下穿隧道振动荷载 数值模拟动力响应Analysis of DynamicProperties of the Soil Induced by theVibration Load

2、 ofTrainsWhen the Undercrossing Tunnel Was Excavated to Different LocationsGao Xuan tao(China Railway First Survey and Design Institute Group Ltd.,Xian 710043,China)Abstract: The article was based on the engineering of a tunnel beneath a high-speed railway, the construction of the tunnel was a premi

3、um source of hazard, so the settlement of the strata was limited strictly. In this paper, the dynamic response of strata under the dynamic train loads during the construction of the tunnel beneath existing railway was computed by means of digital simulation. The simulation results indicate thatthe c

4、rossing point of the tunnel and the high-speed railways roadbed was the weakness, some measurements should be taken to ensure the safety of the construction and the normal operation of the existing high-speed railways. Key words:undercrossing tunnelvibration loaddigital simulationdynamic response近年来

5、，我国轨道交通得到大力发展。随着我国轨道交通的大规模建设，近接工程数量越来越多，施工难度也越来越大。特别是隧道下穿高速铁路施工时，不可避免地对周围土体产生扰动，引起周围地层损失及路基沉降，造成高速铁路轨道的纵、横向不平顺，从而影响列车的运营安全。本文以某隧道下穿高速铁路项目为研究背景，主要采用数值模拟的方法对下穿隧道施工至不同位置时列车动载作用下地层动力响应进行研究分析，为下穿隧道施工提供理论指导。1 工程概况该隧道下穿高速铁路A（即将投入运营）、客运专线B（已运营）、客运专线C（在建）段，需与在建的工程同步完成施工。隧道全长300m，隧道内为单面坡，隧道开挖拱顶到上部铁路轨顶面最大埋深约2m

6、。隧道下穿段与地表高速铁路A、客运专线B及客运专线C在平面位置上呈斜交，斜交角度约30。高速铁路A本段地基采用CFG桩加固, 桩径，桩长，桩底标高约，桩端一般进入强风化层不少于，桩间距，按正方形布置，桩顶设扩大桩头，桩顶设厚垫层（厚中粗砂+厚碎石）,垫层内铺设二层双向经编高强土工格栅；土工格栅极限抗拉强度不小于100kN/m。地基表层挖除换填厚C组填料后再打桩。客运专线B地基采用CFG桩加固, 桩径，桩长9.5，CFG桩桩端进入硬层不少于；桩间距，按正方形布置，桩顶设直径为扩大桩头；桩顶设厚垫层（厚中粗砂+厚碎石）,垫层内铺设二层双向经编高强土工格栅；土工格栅极限抗拉强度不小于100kN/m。

7、地基表层挖除换填厚C组填料后再打桩。客运专线C本段地基采用CFG桩加固, 桩径，桩间距，桩长9.0m（桩长嵌入全风化层内不小于），按正三角形布置。CFG桩桩顶采用1mx1m扩大桩头，桩顶设厚碎石垫层，垫层内铺设1层双向经编土工格栅，土工格栅抗拉强度不小于250kN/m。下穿隧道与地表线路路基相对位置及纵断面图分别见图1、图2。2地质概况隧道处地貌为秦淮河及其阶地、岗地、坳谷以及剥蚀丘陵，隧道洞身位于第四系地层中。根据地质资料揭露，上部为人工填土，25m，下部为第四系上更新统粉质黏土层，硬塑，厚度05m，具有弱膨胀性，下伏基岩为白垩系砂岩或侏罗系凝灰岩、砂岩，节理裂隙发育，全风化弱风化，地下水主

8、要赋存于土层与风化层、强风化与弱风化层间及岩层裂隙内，水量较丰富。隧道拱顶以上03m范围为弱风化岩，隧道所处地质根据围岩情况判定为级，根据隧道穿越地段地表建筑物的重要性及地下水情况，将隧道按为V级型式加强支护。如图3所示。图1工程平面位置（单位：m）图2 工程纵断面（单位：m）图3 地质柱状图（单位：m）3动力分析本构模型土的动力本构关系是极其复杂的，它在不同的荷载条件、土性条件下会表现出极不相同的动力本构特性。目前具体建立的动力本构模型已达数十种，大致可以分为线弹性模型、黏弹性模型、弹塑性模型、边界面模型、内时模型等。根据弹塑性理论，土体变形可以分为弹性变形和塑性变形两部分，其中弹性变形可以

9、应用广义胡克定计算，塑性变形可以应用塑性增量理论计算。塑性增量理论包含3个基本要素：屈服准则、流动法则和硬化定律，分别规定了塑性应变增量的产生条件、方向和大小。目前基于不同的基本要素已经提出了多种弹塑性本构模型，如Drucker-Prager模型、多重屈服面模型（Nested Surfaces Model）、边界面模型（Bounding Surfaces Model）等。事实上，目前还没有一种本构模型能够模拟各种情况下土的动力非线性特征，并具有实际应用中必要的简单性。因此针对具体的问题应该选择合理而简单的本构模型。ABAQUS是国际上最先进的大型非线性有限元计算分析软件之一，具有强健的非线性计

10、算功能。其提供的非线性模型有Mohr-Coulomb模型、扩展的Drucker-Prager模型、修正Drucker-Prager帽盖模型、（修正）剑桥模型。其中Mohr-Coulomb模型在岩土工程中应用最广泛，模拟结果与实际也较为吻合。因此计算中采用Mohr-Coulomb模型。4 有限元计算模型的建立采用三维有限元模型进行计算。根据工程实际及圣维南原理进行建模，隧道顶部以上覆土为，隧道底部以下为，左右边界分别离隧道中线50m，纵向延伸50m。通过ABAQUS程序建立的计算分析模型如Error! Reference source not found.所示。图4 下穿隧道计算分析模型4.1本

11、构模型土体的非线性本构模型：采用了弹塑性物理模型。衬砌单元采用强度等级为C35的钢筋混凝土，计算时采用Mohr-Coulomb本构模型并将混凝土与钢筋合为一体考虑。为了建模及划分网格方便，将计算区域不同地层简化为同一地层，土体参数进行相应转换。根据下穿隧道段地质勘察资料，计算过程中土体的材料参数取值如表所示。表1材料物理力学参数地层弹性模量/MPa泊松比/u密度/kgm-3黏结力/Pa内摩擦角/ ()地层8020500.4E+525高铁路基2700.32250E+538轨道板3.5E+40.16728504.2 边界条件本工程通过对固定边界、黏性边界条件及无限元边界条件进行计算比较分析，计算结

12、果相差不大，一方面是由于列车振动荷载不是很大，另一方面模型选取的比较大，列车振动荷载传至边界处时已衰减较多。而从计算所需时间上看，固定边界所需时间要少，故选用固定边界进行计算。固定边界条件：土层底部完全固定；左右两侧面限制水平方向的位移，竖向自由；前后两面限制轴线方向的位移，竖向自由；地面完全自由。5动力计算及结果分析 本文针对无砟轨道的结构特点，采用弹性地基梁板模型分析了CRTS型板式无砟轨道结构动力响应，给出了CRTS型板式无砟轨道结构在CRH2列车以时速350km经过时轨下压力时程曲线，确定了列车振过程中的激励荷载。如图5所示。图5 列车轮轨激振力时程曲线本文主要计算了隧道开挖到不同位置

13、时列车以时速350km通过，研究地层动力响应。计算过程中分别在地层表面、地层不同深度处及隧道拱顶设置监测点，监测该处位移变化情况。5.1地层表面监测点位移分析隧道分别开挖至15m、25m和35m时，施加列车振动荷载，地层表面各监测点在列车动载作用下的竖向位移时程曲线如图6所示。(1) 隧道开挖至15m处施加动载地表监测点位移时程曲线(2) 隧道开挖至25m处施加动载地表监测点位移时程曲线(3) 隧道开挖至35m处施加动载地表监测点位移时程曲线图6 隧道开挖至不同位置时施加动载地表监测点位移时程曲线隧道开挖至不同位置施加动载地表监测点最大竖向位移如表2所示。表2地表监测点最大竖向位移监测点最大竖

14、向位移/mm隧道开挖至15m施加动载隧道开挖至25m施加动载隧道开挖至35m施加动载z=15m处监测点z=25m处监测点z=35m处监测点通过对隧道开挖至不同部位施加动载地表监测点位移时程曲线及最大位移值分析可知：（1）在单次列车动载作用下地表监测点产生的竖向位移值与隧道开挖引起的沉降量相比很小，可以忽略不计（隧道开挖引起的最大沉降量为2.7mm）。但是根据国内外工程实际在长期高速列车振动荷载的作用下，地层产生的沉降量是不能忽视的。（2）随着一组组车轮的滚过，监测点位移时程曲线具有明显的车轮滚过效应；当车轮压在监测点时，该监测点出现位移峰值，车轮离开监测点后，位移逐渐减小；当列车驶离后，位移逐

15、渐回弹恢复。（3）当第一节车厢经过时，z=15m处监测点产生的位移最大；随后，z=25m处即下穿隧道与高铁路基交叉点处监测点位移最大，因此，在设计与施工时要注意对下穿隧道与高铁路基交叉部分加强支护。（4）隧道中点前监测点随着隧道开挖的不断进行，施加动载产生的位移值在逐渐减小；隧道中点及中点后监测点则随着隧道开挖的不断进行，施加动载产生的位移值在不断增大。5.2隧道拱顶监测点位移分析隧道分别开挖至15m、25m和35m时，施加列车振动荷载，隧道拱顶各监测点在列车动载作用下的竖向位移时程曲线如图7所示。(1) 隧道开挖至15m处施加动载隧道拱顶监测点位移时程曲线(2) 隧道开挖至25m处施加动载隧

16、道拱顶监测点位移时程曲线(3) 隧道开挖至35m处施加动载隧道拱顶监测点位移时程曲线图7 隧道开挖至不同位置时施加动载隧道拱顶监测点位移时程曲线隧道开挖至不同位置施加动载隧道拱顶监测点最大竖向位移值如表3所示。表3拱顶监测点最大位移值监测点最大沉降值/mm隧道开挖至15m施加动载隧道开挖至25m施加动载隧道开挖至35m施加动载z=15m处监测点z=25m处监测点z=35m处监测点通过对隧道开挖至不同部位施加动载隧道拱顶监测点位移时程曲线及最大位移值分析可知：（1）隧道拱顶位移时程曲线表现出和地表位移时程曲线同样的特性；（2）隧道开挖至某一处施加动载，该处拱顶产生的位移值最大，因此，隧道开挖后，

17、特别是上台阶开挖后，要及时施加支护；（3）随着开挖的不断推进，开挖面即掌子面处拱顶位移值在不断增大。5.3地层不同深度监测点位移分析为研究在列车动载作用下，地层不同深度处动力反应，取隧道开挖面即掌子面处地层不同深度监测点来进行研究，监测点布置如图8所示。隧道开挖至不同位置时施加动载，在动载作用下各监测点的位移时程曲线如图9所示。图8不同深度监测点分布示意(1) 隧道开挖至15m处施加动载地层不同深度监测点位移时程曲线(2) 隧道开挖至25m处施加动载地层不同深度监测点位移时程曲线(3) 隧道开挖至35m处施加动载地层不同深度监测点位移时程曲线图9 隧道开挖至不同位置时施加动载地层不同深度监测点

18、位移时程曲线隧道开挖至不同位置施加动载地层不同深度监测点最大位移值如表4所示。表4地层不同深度监测点最大位移值监测点最大位移值/mm隧道开挖至15m施加动载隧道开挖至25m施加动载隧道开挖至35m施加动载监测点1监测点2监测点3监测点4监测点5监测点6通过对隧道开挖至不同部位施加动载地层不同深度监测点位移时程曲线及最大位移值分析可知：（1）随着地层深度的增加，在动载作用下产生的位移不断减小，地表处监测点位移最大，隧道拱顶处位移最小。这说明振动荷载在地层中传播时，由于地层阻尼等因素使振动波能量不断减小。（2）地层动力响应从地表到隧道拱顶衰减速度逐渐减小。（3）隧道开挖至25m处施加动载，地层不同

19、深度各监测点的位移值较大。5.4地层应力响应分析为分析在列车通过时地层的竖向应力变化情况，在地层不同深度布置监测点，如图8所示，监测点位于隧道开挖断面处即掌子面处。通过对地层监测点进行监测分析，得到不同深度处地层在列车振动荷载作用下的竖向应力时程曲线，如图10所示。(1) 隧道开挖至15m处施加动载地层不同深度监测点应力时程曲线(2) 隧道开挖至25m处施加动载地层不同深度监测点应力时程曲线(3) 隧道开挖至35m处施加动载地层不同深度监测点应力时程曲线图10 隧道开挖至不同位置时施加动载地层不同深度监测点应力时程曲线注：监测点从地表到隧道拱顶编号依次为监测点16。隧道开挖至不同位置施加动载地

20、层不同深度监测点最大应力值如表5所示。表5 隧道开挖至不同位置施加动载地层不同深度监测点应力最大值监测点最大竖向应力值/Pa隧道开挖至15m施加动载隧道开挖至25m施加动载隧道开挖至35m施加动载监测点1-1250-973-924监测点2-1111-930-856监测点3-947-833-751监测点4-799-720-651监测点5-658-606-557监测点6-327-360-368通过对隧道开挖至不同部位施加动载地层不同深度监测点应力时程曲线及应力最大值分析可知：（1）随着地层深度的增加，在动载作用下产生的竖向应力不断减小，地表处监测点沉降最大，隧道拱顶处沉降最小。这说明振动荷载在地层中传播时，由于地层阻尼等因素使振动波能量不断减小。（2）地层动力响应从地表到隧道拱顶衰减速度逐渐减小，但是隧道拱顶附近衰减较快，这主要可能是由于隧道拱顶的压力拱效应造成的。（3）随着隧道开挖的不断推进，开挖面即掌子面处不同深度地层在列车动载作用下的竖向应力在不断减小，这主要是因为隧道开挖深度越深，高铁路基与下穿隧道交叉处的“临空面”面积不断增加，地层在列车动载作用下向隧道净空发生位移，使应力在一定程度上释放。6结论通过对隧道施工至不同位置时列车以时速350km通过，研究地表及地层不同位置处位移值的变化，结论如下: