盾构隧道掘进数值模拟
- 格式:pdf
- 大小:1.10 MB
- 文档页数:65
专业知识分享版使命:加速中国职业化进程摘 要:针对某盾构隧道下穿既有地铁暗挖隧道的施工力学行为进行了三维有限元数值模拟分析。
研究结果表明: 在盾构推进至距既有隧道边缘3 m 前,隧道会发生隆起,且在此位置时隆起量最大,之后开始沉降,在盾构将要穿出既有隧道时,沉降增量最大; 隆起量随盾构推力和既有隧道刚度增大而增大,而沉降量与之相反; 盾构下穿时,既有隧道结构横截面上会产生扭转,扭转角的大小随盾构推力增大而增大,随既有隧道刚度增大而减小。
为确保下穿过程上方隧道的结构安全和列车的正常运行,在距既有隧道边缘 3 m 时采取措施控制盾构推力和提高既有隧道周围土体的强度非常有效。
关键词:隧道 盾构 下穿 数值模拟 竖向位移 横向扭转随着城市地下轨道交通的发展,下穿既有线路的情况时有发生。
由于新线穿越既有线不可避免地会引起既有隧道结构产生附加应力和沉降,而地铁运营又对既有线的轨道变形有非常严格的控制标准,依据《上海市地铁沿线建筑施工保护地铁技术管理暂行规定》,运营隧道结构水平和沉降最大位移应 < 20 mm;根据《铁路线路维修规定》,轨道纵向每 10 m 的沉降差应 < 4 mm 。
因此这类下穿工程对既有地铁的安全形成了严峻的考验。
分析盾构下穿既有隧道的力学行为是非常必要的。
姜忻良、赵志明等[1]用理论推导的方法,提出隧道开挖时,上覆土在不同深度处的沉降计算公式; 汪洋、何川等[2]利用模型试验和数值分析的方法,并考虑了隧道纵向和横向刚度的折减得到围岩、净距、推力对上覆隧道的位移和附加应力的分布规律; 文献[3-5]利用三维数值模拟的方法提出下穿盾构的推进与上方近接隧道结构位移的关系; 陈越峰、张庆贺等[6]通过数值模拟及实测数据的反馈,找出了上覆隧道的沉降规律; 张海波、殷宗泽等[7]运用三维数值模拟的方法得出上、下隧道间的距离和相对位置对彼此的影响。
本文以某盾构隧道下穿既有暗挖隧道工程为背景,采用有限元软件建立三维数值分析模型对盾构下穿既有隧道全过程进行动态模拟。
盾构隧道急曲线段施工数值模拟分析发表时间:2019-06-20T09:39:48.747Z 来源:《防护工程》2019年第6期作者:杜亭萱[导读] 本文利用有限元软件对隧道后靠土体的稳定性进行了估算,并且对周边建筑物的沉降进行了分析,为其他类似项目提供参考。
上海市地下空间设计研究总院有限公司 200020摘要:本文采用数值模拟的方法,对大直径急曲线隧道施工过程进行了数值模拟分析。
本工程的最大难点是大直径盾构隧道在S曲线小半径急转弯处运行,平面曲线最小半径仅为500m。
在掘进曲线段过程中,内、外千斤顶的受力有一定的差别,盾构推力通过管片传递到盾构后靠土体,可能引起后靠土层的失稳。
本文利用有限元软件对隧道后靠土体的稳定性进行了估算,并且对周边建筑物的沉降进行了分析,为其他类似项目提供参考。
关键词:大直径盾构隧道;后靠土体稳定性;周边建筑沉降1工程概况1.1项目总体情况该工程I标段隧道长约8km,其中盾构段约6km,主线设工作井3座,在这三个工作井之间的隧道直线段占本区间隧道长度仅为20%。
东西盾构隧道采用单管圆形隧道,管片外径为15m,管片厚度为650mm。
盾构隧道最小曲线半径为500m。
1.2工程地质与水文条件根据岩土勘察报告,场地90m以内分布的土层自上而下的土层分别为①层为填土,②1层~⑤3层为全新世Q4沉积层,⑥层~⑨层为上更新世Q3沉积层。
地下水主要有赋存浅层中的潜水,⑤1、⑤2层中的微承压水和⑦层、⑨层中的承压水。
盾构主要穿越地层为:⑤1粉质黏土、⑤3-1粉质黏土夹粉砂、⑥粉质黏土、⑦1粉砂、⑦2粉细砂、⑧1-1黏土等。
1.3周边环境情况在分析区段内,区间隧道的转弯半径只有500m。
且周边环境较为敏感,沿线建筑众多,下穿别墅区、公寓及政府管理中心。
这些建筑均采用浅基础形式。
2盾构隧道施工有限元模拟2.1土体本构模型为了更加精确的模拟施工过程,土体采用修正摩尔库伦模型(Modified Mohr-Coulomb)。
迎坡条件下地铁盾构隧道施工数值模拟分析摘要:主要采用数值分析方法,以广州地铁为工程背景,计算采用有限差分程序FLAC3D进行施工模拟,对地层位移、隧道洞周位移、开挖面应力、开挖面附近塑性区的计算结果进行了分析。
关键词:盾构;数值分析;数值模拟随着数值计算被普遍应用于工程设计中解决各种岩石力学问题,各种数值模拟技术在岩土力学中有了很大的发展和广泛的应用。
然而,这些数值分析方法其理论本身以及采用的算法都有着各自的局限性。
快速拉格朗日分析(Fast Lagrangian Analysis of continua,简称FLAC),则是在较好地吸取上述方法的优点和克服其缺点的基础上形成的一种数值分析方法。
1盾构隧道施工过程的数值模拟方法1.1掘进过程的模拟(1)开挖的实现FLAC3D中有内置一种“零模型”,通常用来模拟岩土体被开挖或移除。
零区域的应力被设置为零,在这些区域中没有体力作用。
可以通过单元性质的改变来实现后期的回填或支护。
(2)管片安装采用零模型来模拟开挖时,结合盾构机的施工推进速度,按每环1.5米进行开挖。
管片采用弹性圆环体模型来模拟,杨氏模量按《混凝土结构设计规范》选取为35.SGPa,考虑节头效应弹性模量乘以折减系数0.65,取为18GPa,钢筋混凝土的泊松比为0.17。
土体采用D-P模型来模拟。
盾构隧道开挖过程见图1所示,管片的拼装落后于开挖面4-6m。
在FLAC3D中,岩土体的力学参数采用的是体积模量K和剪切模量G,因此,须在参数输入前进行转换。
图1盾构隧道开挖示意图(3)工作面支护力的施加盾构施工工作面支护力应在极限主动土压力和极限被动土压力之间,在数值计算中,参照现场工作面支护压力的统计结果,确定工作面支护力的具体取值(不同坡度条件下其取值不同)。
1.2注浆层的简化(1)盾尾空隙计算盾构推过过程产生的土体沉降主要是由地层损失引起的,即隧道施工中实际开挖的土体体积与竣工体积之差,横断面的地层损失可以通过空隙厚度G来表示,Loet(1984)作过相关论述:(1)其中,(2)上式中△为盾尾壳体的厚度,为安装衬砌所需的空隙厚度,为开挖面应力释放导致开挖面及其前方土体的三维运动,使得土体塌落到开挖面造成的超挖土量,为施工因素(包括盾构的纠编、叩头、后退)及操作技术的影响产生的土体损失,为考虑土体后期固结产生的土体位移。
Indust rial Const ruct ion Vol.40,Supplement,2010工业建筑 2010年第40卷增刊盾构法施工过程的三维有限元仿真研究赵凤山(中铁十一局集团第一工程有限公司,湖北襄樊 441104)摘 要:盾构施工会导致地层的沉降并会对周边的建筑产生影响,对其进行分析具有重要的工程价值。
在综合考虑了现有的有限元模拟方法的基础上,结合济邵高速公路中开挖的隧道为例,根据非线性有限元法的基本原理,利用大型非线性有限元软件ANSYS 三维动态模拟隧道开挖的过程,计算结果揭示了在盾构推进过程中地表沉降分布以及特点,得到盾构法开挖引起的地表沉降曲线,其形态与Peck 计算得出的横向地表沉降槽正态分布曲线的形态基本一致。
仿真结果揭示了隧道推进过程中的周围土体的移动规律、地表面水平位移及地面沉降规律、已建隧道管片应力变化规律。
关键词:盾构法;公路隧道;三维有限元法;地表沉降THREE 2DIMENSIONAL FINITE ELEMENT ANALYSIS OF PROCESSOF SHIELD 2DRIVEN TUNNELZhao F engshan(China Railway Bureau 11Gr oup the F irst P roject Co 1Ltd,of No 111bureau Group,China Railwa y,Xiangfan 441104,China)Abstract:A shield const ruct ion in tunnel can cause ground subsidence and impact on the buildings ar ound,whose ana lysis is very impor tant 1Based on the three dimensional nonlinear finite element model,taking a Ji 2Shao freeway shield tunneling for example ,using ANSYS ,a la rge nonlinear FEM software to do a 3D FEM model simulation in the process of the shield constr uction ,then t he curve of subsidence by shield met hod construction had been calculated 1T he ground subsidence distr ibut ing char acteristics had been revealed 1T he curve of subsidence by the disturbance of soil was basically in accordance with nor mal distribution curve by Peck 1T he r ules of soil movement ,gr ound settlement and hor izontal displacement as well as change in str ess of t he adjoining tunnel during the shield 2dr iven constr uct ion wer e revealed.Keywor ds:shield method;highway tunnel;3D finite element method;ground subsidence作 者:赵凤山,男,1972年6月出生,工程师。
第30卷,第6期 中国铁道科学Vol 130No 16 2009年11月 C HINA RA IL WA Y SCIENCENovember ,2009 文章编号:100124632(2009)0620054207盾构法隧道下穿既有结构三维数值模拟分析杨广武1,关 龙2,刘 军2,郑知斌3(1.北京交通大学,北京 100044; 2.北京市政集团,北京 100045; 3.北京市市政工程研究院,北京 100037) 摘 要:以北京地铁10号线9标段盾构法开挖隧道穿越城铁13号线芍药居车站工程为依托,采用三维有限差分软件FL AC 3D 对盾构施工过程进行数值模拟,分析盾构穿越既有结构时对其沉降的影响规律。
研究结果表明:既有线车站结构的沉降随地基变形模量的提高而减少,且沉降趋势逐渐变缓;既有结构最大沉降增大的速率比围岩荷载释放率增大的速度快;增大开挖面的控制压力,可有效减小既有结构的沉降,但过大的控制压力会使前方土体隆起,产生负地层损失,并且随着开挖面控制压力的提高,差异沉降明显增大。
通过施工参数的优化可以减小既有结构的沉降,达到保护既有结构的目的。
关键词:隧道;盾构法;既有结构;变形;数值模拟 中图分类号:U455143 文献标识码:A 收稿日期:2009205212;修订日期:2009209211 基金项目:北京市政总公司(集团)基金资助项目(科2J 206116) 作者简介:杨广武(1956年—),男,北京人,教授级高工,博士研究生。
随着城市轨道交通的不断发展,换乘节点不断增加,必然出现新建地铁隧道近距离穿越既有地铁隧道或车站的现象,即所谓的近接施工。
控制新建隧道穿越既有地铁隧道或车站所引起的变形,确保既有地铁隧道或车站的结构安全和新建隧道的顺利掘进,是地铁施工面临的重要课题。
土压平衡盾构施工较为突出的问题是盾构挤压推进对周围土体的扰动较大,合理设置和控制土压对于控制地表沉降至关重要[123]。
盾构近距离隧道叠交施工数值模拟分析同济大学廖少明2004年1月目录一、研究问题简述二、计算方法三、土工试验四、有限元计算结果及分析五、施工隧道对已建隧道影响的研究六、盾构隧道地面沉降计算经验公式研究一、研究内容主要对下列问题进行深入研究:(1)隧道推进过程中的周围土体的移动规律及地面沉降规律(2)相邻隧道轴线位移规律(3)相邻隧道变形规律(4)土压力、注浆压力等对上述规律的影响(5)软土流变性对施工后土层和隧道应力、变形的长期影响。
Back二、计算方法1、盾构施工过程的三维有限元模拟方法1.1 盾构推进的模拟1.2 有限元网格剖分1.3 计算中考虑的主要影响因素2、计算程序3、材料本构模型3.1邓肯-张E-ν非线性弹性模型3.2 殷-Graham土体三维弹粘塑性模型Back盾构推进引起周围土体的变形机理图图中1、2、3、4、5分别代表五种不同的变形机理:1——由于盾构推进中正面土压力的不平衡而导致地层下沉或隆起,以及开挖面的崩裂;2——由于盾构外壳与土体之间摩擦而导致地层隆起;3——由于盾构姿态的变化引起地层损失而导致地层下沉;4——由于盾构推进后的注浆引起的地层隆起及下沉;5——由于以上四种作用,盾构推进后使周围土体产生超孔隙水压力和受到扰动而进行固结和蠕变导致地层下沉。
1.3 计算中考虑的主要影响因素在本次分析计算中主要考虑了以下一些由于盾构推进引起隧道周围土体位移和地面沉降的影响因素:①开挖面土体应力增加或释放。
②盾构壳周围土体的剪切。
③盾尾空隙及盾构刀盘外径大于盾构外径引起的超挖。
④盾尾注浆。
⑤软粘土次固结。
Back Continue1.1 盾构推进的模拟假设盾构跳跃式向前推进,用改变单元材料类型的方法(刚度迁移法)来反映盾构的向前,同时施加相应的荷载。
当前盾构机位置下一步盾构机位置卸荷单元 卸荷单元管片单元盾尾空隙单元管片单元盾尾空隙单元浆液单元浆液单元开挖面开挖面1、盾构施工过程的三维有限元模拟方法2、计算程序计算采用的三维有限元程序为ANSYS,MARC,FLAC等计算软件。
根据盾构施工的特点,对盾构的推进、盾尾空隙与压浆等过程进行模拟,从而分析新隧道盾构对已建隧道的影响1.2 有限元网格剖分将该分析区域的分为以下四段:水平平行推进段,由水平平行向垂直平行过渡段前、后两段,垂直平行推进段。
水平平形向垂直平行过渡前半段水平平行垂直平行后半段1.2 有限元网格剖分(续1):各段三维有限元网格1.2 有限元网格剖分(续2):有限元网格过渡情况示意3、材料本构模型本研究项目的计算中,土体本构模型主要采用了邓肯-张E-ν非线性弹性模型。
混凝土管片和盾构机钢壳作为线弹性材料处理。
流变计算采用了殷-Graham 土体三维弹粘塑性模型。
3.1邓肯~张E ~ν非线性弹性模型na a f t p p K c R E ⎟⎟⎠⎞⎜⎜⎝⎛•⎥⎦⎤⎢⎣⎡+−−−=32331sin 2cos 2))(sin 1(1σϕσϕσσϕ23)1(lg A p F G a t −⎟⎠⎞⎜⎝⎛−=σν四、有限元计算结果及分析1、隧道周围土体的位移场2、地表土体位移3、已建隧道管片变形和应力4、弹粘塑性模型计算结果5、新老隧道先后推进综合影响下的地面沉降Back1、隧道周围土体的位移场水平平行推进情况下纵、横剖面位移矢量图0-10-20-30-40-5001020304050607080901001101201301401501600.00 3.46-10-20-30-40-50-80-70-60-50-40-30-20-100102030405060700.00 2.222、地表土体位移水平平行推进情况下纵、横剖面地表土体沉降曲线-30-25-20-15-10-505-80-60-40-20204060(m)s (m m )开挖面前方10m 开挖面盾尾处盾尾后20m沿新隧道轴线纵断面地表沉-30-25-20-15-10-5050.020.040.060.080.0100.0120.0140.0160.0180.0(m)s (m m )计算结果3、已建隧道管片变形和应力开挖面前方10m盾尾后方20m2.60~3.202.00~ 2.59-0.50~-0.50-0.80~-0.50-1.99~-1.50-2.50~ -2.005、新老隧道先后推进综合影响下的地面沉降-4.0-3.0-2.0-1.00.0-70-40-10205080地表点位置(m)地面沉降(cm)由老隧道引起由新隧道引起新、老隧道综合影响-4.0-3.0-2.0-1.00.0-70-50-30-1010305070地表点位置(m)地表沉降(cm)由老隧道引起由新隧道引起新、老隧道综合影响-6.0-5.0-4.0-3.0-2.0-1.00.0-75-45-15154575地地由由新-10.0-8.0-6.0-4.0-2.00.0-75-45-15154575地表点位置(m)地面沉降(cm)由老隧道引起由新隧道引起新、老隧道综合影响Back五、施工隧道对已建隧道影响的研究1、假想新建隧道向已建隧道左下部过渡1)隧道周围土体的位移场2)已建隧道管片变形和应力2、施工隧道对已建隧道影响规律研究Back1、假想新建隧道向已建隧道左下部过渡1)隧道周围土体的位移场-10-20-30-40-50-80-70-60-50-40-30-20-10010203040506070-10-20-30-40-50-60-50-40-30-20-1001020304050601)隧道周围土体的位移场(续)-10-20-30-40-50-60-50-40-30-20-1001020304050602)已建隧道管片变形和应力1.60~1.800.80~1.20-0.50~0.50-1.10~-0.70-1.49~-1.10Back-1.80~-1.502、施工隧道对已建隧道影响规律研究已建隧道管片沉降与压缩模量的关系(左右平行)-1.4-1.2-1.0-0.8-0.6-0.4-0.20.005000100001500020000压缩模量Es(kPa)沉降(c m )H=32mH=27m已建隧道管片沉降与隧道埋深的关系已建隧道管片沉降与隧道埋深的关系(左右平行)-0.55-0.50-0.45-0.40-0.35101520253035隧道上覆土层厚度(m)沉降(c m )Es=6700kPa已建隧道管片位移与隧道间夹角的关系已建隧道管片指向新建隧道方向位移与隧道间夹角的关系(间距4m )0.30.60.91.21.51.8-90-60-30306090夹角α(度)位移(cm)已建隧道管片最大位移已建隧道管片最大位移与隧道间距离的关系已建隧道管片最大位移与隧道间距离的关系(新在上)0.20.30.40.50.60.70.80.901234567隧道间净间距(m)已建隧道管片位移(c m )已建隧道管片最大位移与隧道间距离的关系已建隧道管片最大位移与隧道间距离的关系(新在下)-2.2-1.8-1.4-1-0.6-0.201234567隧道间净间距(m)已建隧道管片位移(c m )Back六、边界元分析方法Back盾构近距离穿越施工力学解析(BEM)∫=S dSBBAGAu)(),()(φ∫∫+ =1221)(),()(),(S SdSBBAGdSBBAGφφ(2-1)∫=S dSBBAFAt)(),()(φ∫∫+ =1221)(),()(),(S SdSBBAFdSBBAFφφ(2-2)∫∑∫=ΔΔ=11111)(),()(),(S mi S i i S d B B A G dS B B A G φφ (2-3)∫∑∫=ΔΔ=22122)(),()(),(S nj S j j S d B B A G dS B B A G φφ (2-4)∫∑∫=ΔΔ=11111)(),()(),(S mi S i i S d B B A F dS B B A F φφ (2-5)∫∑∫=ΔΔ=22122)(),()(),(S nj S j j S d B B A F dS B B A F φφ (2-6)地面1隧道1d2隧道2dc图.1 弹性半无限土体中的隧道截面图图.2 隧道边界单元rtg Rrr tg R tg r R L −++−+−=Δθθθ222212)(0102030405060-30-20-100102030水平距离x(m)竖向位移(m m )z=0.0mz=5.0m z=10.0m1217.6m26.5m6.4m图.7 土层竖向位移计算值(XK12+489)102030405060-1001020水平位移(mm)深度z (m )x=10.0m(离原点)x=12.0m(离原点)x=15.0m(离原点)26.4m117.6m26.5m102030405060-1001020水平位移(mm)深度z (m )图.10 土层水平位移计算值(01020304050607080-30-20-100102030水平距离x(m)竖向位移(m m )L=7.0mL=15.0m L=20.0m L=25.0m1215.0m L15.0m图.15 隧道水平位置变化引起的地表沉降槽七弹性边界影响问题-16.285Φ1400覆土2.1m12.3m9.2m Hc=20.1m-14.5851.5m-10.39地铁二号线区间隧道6.2m11.2m 14.3mHc=21.3m1.8m-17.479地铁二号线区间隧道6.2m-12.356-16.997-8.5-18.006-18.011-10.53.8图4 M4线盾构穿越M2线纵剖面图图2 已建隧道弹性边界S2刚度计算模式u u K u /σ= (10)w w K w /σ= (11)xyzoL地平面地平面Hh盾构S 2S 1已建隧道24)1()(65.012g geq g E EI D E K μ−=-0.6-0.5-0.4-0.3-0.2-0.10.00.120406080100G Wx /m间距为50 m间距为5 mz = 10.0m间距为40 m 间距为50 m间距为20 m 间距为10 m实际工程中弹性边界的存在,会改变盾构推进时引起的位移场的分布。
通过计算表明:当盾构向弹性边界推进时,土层中位移波逐渐向前推进的同时,位移波中的波峰逐渐升高,且在盾构与弹性边界的距离达到某关键值时位移达到最大值,随后波峰开始下降;另外,盾构与弹性边界之间的竖向位移场的叠加作用及弹性边界对其后竖向位移场的屏蔽作用也越来越有效果。
隆沉分析各因素影响隆起比较由此表可得结论:盾构头部对土体隆起的影响是盾构侧壁影响的两倍左右。
工况4的纠偏情况:正面顶力引起的最大隆起:侧壁摩阻力起的最大隆起:纠偏压力起的最大隆起0.29 : 3.12 :0.136正面顶力引起的最大沉降:侧壁摩阻力起的最大沉降:纠偏压力起的最大沉降1.34:14.74 :0.22.540.411.0442.290.481.123 2.040.541.102 1.310.700.911b1/b2 盾构侧壁影响b2盾构头部影响b1工况八、隧道管片结构数值分析•梁—弹簧模型和梁—接头模型荷载假设与上图的区别,主要是在地基抗力的计算上,修正惯用法中假设垂直方向上的地基抗力与地基位移无关,水平方向上的地基抗力,则是伴随衬砌向围岩方向的变形而产生。