浅埋暗挖隧道施工性态的数值模拟与分析.docx

浅埋暗挖隧道施工性态的数值模拟与分析.docx

《浅埋暗挖隧道施工性态的数值模拟与分析.docx》由会员分享，可在线阅读，更多相关《浅埋暗挖隧道施工性态的数值模拟与分析.docx（5页珍藏版）》请在冰豆网上搜索。

浅埋暗挖隧道施工性态的数值模拟与分析

浅埋暗挖隧道施工性态的数值模拟与分析

　　摘要:

地铁隧道施工中,隧道开挖程序、施作的步骤对隧道稳定性及地表沉降影响显着。

北京地铁10号线某标段浅埋单洞双层隧道,是出入口通道与车站中洞、旁边单层侧洞相接的暗挖部分,单、双层转换,初支结构形式和施工都很复杂。

目前国内外对单洞双层隧道分析较少,与车站和通道相连接的这种单洞双层隧道比较短,对它的分析往往被忽视。

对采用3层6导洞的CRD工法施工的这种浅埋单洞双层隧道的施工性态进行数值模拟,分析该区域隧道施工引起的地表沉降,探索一次性开挖成洞及隧道分步开挖引起的地表沉降槽的变化形态,拱顶沉降及拱顶主应力的变化规律。

数值分析表明:

施工工序不同使隧道的偏挖引起的沉降槽向未开挖一侧偏移,此过程中地层最大沉降并不发生在隧道中线处,而是完成全部开挖,地层变形稳定后,其累计沉降最大值位于隧道中线处。

隧道开挖衬砌完全施作后,地表沉降变化不大,但是后续开挖步引起的结构内力应予以重视,尤其对于中隔壁支撑及拱底衬砌的支护,要及时施作拱底二次衬砌。

研究结果为单洞双层隧道分步施工控制地表变形和洞周位移提供参考依据。

关键词:

浅埋暗挖隧道;单洞双层;施工性态;数值模拟;地表沉降

引言

在地铁隧道施工中,由于隧道开挖将造成土体应力重分布,因此在一定的土体环境中,隧道的开挖方式对隧道稳定性及地表沉降影响尤为显着[1]。

不同的施工程序,在时空上相当于荷载以不同的方式施加在隧道上,尤其对于浅埋隧道,覆土厚度较小,围岩松散,自承能力差,隧道开挖引起很大的地表沉降,往往造成周边建筑物破损、倾斜甚至倒塌。

因此研究隧道施工工序对隧道稳定及地表沉降的影响具有重要意义。

本文即对北京地铁10号线某标段单洞双层通道的施工工序引起的地表沉降及拱顶沉降进行了数值模拟与分析。

1工程简介

北京地铁10号线某标段与车站结构相联部位设4个出入口通道,每个通道在中洞和侧洞之间的部分为单洞双层通道,其余部分为单层通道,均采用暗挖法施工。

东南和东北出入口通道,中洞和侧洞之间双层部分长度为　m,双层通道的开挖尺寸宽×高为6600mm×12000mm,初期支护结构厚度为300mm,二衬结构厚度为500mm,二衬净空尺寸宽×高为5000mm×10100mm。

见图1,土层参数见表1。

2地表变形分析

地表沉降理论分析

预测地面沉降量及影响范围的常用的方法有边界元法、有限元法、根据实测数据的统计方法等。

其中应用较多的是Peck提出的经验公式:

式中:

S（x）为距离隧道中线x处的地面沉降量,m;S（max）为隧道中线x=0处的地面沉降量,m;x为距隧道中线的距离,m;i为沉降槽的宽度系数,即沉陷曲线反弯点的横坐标,m。

Peck假定横向沉陷曲线为正态分布曲线。

当横向沉陷曲线为正态分布曲线时,Smax与沉降槽体积VS满足下列关系:

沉降槽的宽度系数i取决于接近地表的强度,隧道埋深和隧道半径。

通过在均匀介质中试验,得到近似的几何关系是:

式中:

Z为隧道开挖面中心至地面的距离,m;R为隧道半径,m;K,n为试验系数,K=￣,n=￣。

由以上经验理论可知,地层沉降沿深度在任一横断面呈正态分布,最大沉降位于隧道中心位置处,随着隧道埋深Z的不断增加,地表最大沉降值也不断增大,在隧道拱顶位置处沉降值达到最大。

数值模拟

本文分析采用了Mohr-Coulomb塑性模型,有限差分计算方法对单洞双层隧道开挖过程进行模拟,其中初衬、临时支护均用梁单元,小导管注浆作用通过提高土体的弹性模量及c、Ф值来进行模拟。

目前世界上单洞双层隧道并不多见,但是在隧道的出入口与车站中洞与侧洞相连的部位,单洞双层通道较多,只是长度较短,施工往往被忽视。

而此处恰是车站与侧洞单洞单层通道的过渡,对它施工所引起的地层扰动,直接影响到测洞和车站的施工。

本文完全模拟了双层通道施工的动态过程,先进行第①步拱部超前小导管注浆,台阶法开挖①步土体,铺设Ф6-150×150钢筋网,架格栅并作格栅锁脚锚杆,隧道支护参数见表2,喷混凝土作初衬。

然后依次进行第②～⑥步的施工,见图2。

分部开挖引起的地表沉降分析

在进行分步开挖计算时,取隧道中线为零点,开挖水平影响范围约为隧道跨度的15倍,以减小模型的边界约束效应,使模型更能真实地反映开挖引起的沉降,竖向影响范围为隧道开挖深度的5倍,建立了100m×60m的网格[5-6],见图3,以隧道中心点为对称点对地表18个点进行了地表竖向位移的监测和数据采集,有关分步开挖地表的沉降分布形态见图4。

第一步导洞①开挖,地表最大沉降量模拟计算值为36mm,地表沉降呈正态分布,现场实测最大沉降量为　mm,第二步开挖地表模拟计算值最大沉降量为　mm,实测值为　mm,沉降槽的宽度没有变化,水平影响范围也没有增加,在-40m～+40m处沉降曲线几乎重合,说明随着隧道进一步开挖,地表沉降曲线沉降槽的宽度没有变化,沉降量略有增加。

本文也模拟了当第二步开挖先开挖导洞③后地表变形的情况,开挖面上对应的地表出现了mm的隆起,与第一步开挖导洞地层的变形完全相反,沉降槽的宽度也明显增加,不利于地层的稳定沉降。

　　进行第三步导洞③的开挖后,地表最大沉降发生在距离隧道中线-　m处,沉降量为　mm,沉降槽右移;第四步导洞④开挖完成后,地表的最大变形发生在距离隧道中线　m处,隆起量为mm,同时在距离隧道中线-　m处发生了mm的地表最大沉降,沉降曲线表明右侧导洞④的开挖引起沉降槽向x轴负方向移动,但并不是简单的平移,以隧道中心线为中心0～40m之间的地表均有不同程度的隆起,大于20m的范围地表沉降的趋势与导洞③开挖地表沉降的趋势一致,沉降槽的最低点在中心线附近,沉降曲线表明地层变形很平缓。

说明该施工过程有利于地表沉降的控制。

第五步导洞⑤开挖后,地表最大变形发生在距离隧道中线　m处,最大沉降量为　mm,最大隆起发生在距离隧道中线-　m处,隆起值为　mm,沉降槽的宽度比第一步、第二步开挖时沉降槽的宽度增加20%左右,并且沉降槽偏移,与实际施工性态相符,此时正是进行导洞⑤的开挖;第六步导洞⑥开挖,隧道中心点附近-　m处地表发生　mm的沉降,最大隆起发生在距离隧道中线-　m处,隆起值为　mm,沉降槽的宽度有所收敛,比第一步、第二步开挖的沉降槽宽度减小5%左右,并且沉降槽向隧道中线偏移,此时单洞双层隧道开挖完毕。

从以上数值分析,可以看出分步施工中,隧道先期导洞偏挖时,沉降槽向未开挖一侧偏移;在右侧导洞③④开挖时,可以认为开挖沉降槽是新形成的,但是新的沉降槽反映出导洞③④开挖对地层的扰动没有初始开挖对地层的扰动大,并且新的沉降槽是受先期地层沉降影响的,有利于对施工引起的地表沉降进行控制。

图5是模拟第三步施工先开挖导洞⑤引起的地表沉降曲线,随着隧道开挖深度的增加,地表沉降继续发展,最大沉降量达60mm,在开挖跨度　m范围内,地表沉降均在60mm左右,说明此时衬砌受到的围岩压力很大,沉降槽的宽度收敛,表明沉降槽的宽度与形状受隧道开挖面的形状影响较大。

在开挖过程中,拱顶的沉降见图6,随着施工步的进行,拱顶的沉降呈线性,但是距离隧道中线-3m的拱顶沉降变化很小,随着施工步的进行,拱顶的沉降值逐渐变为正值,与实测相符。

拱顶由于小导管注浆加固地层,衬砌刚度比较大,因此拱顶的地层变形呈线性。

拱顶土体主应力的变化见图7,开挖初始第一、第二步,拱顶土体主应力最大值为180kPa,距离隧道中线-　m处,主应力值为零。

随着开挖步的进行,隧道初衬的施作,远离隧道中线处隧道右侧,拱顶土体主应力逐渐增加,隧道中线左侧第一开挖步15m范围内,拱顶土体主应力变化不大,基本一致。

3结语

本文以北京地铁10号线某标段车站与出入口衔接处的单洞双层隧道工程为研究对象,对该区域隧道施工引起的地表沉降进行了数值模拟分析,探索了单洞双层隧道施工一次性开挖成洞及隧道分步开挖引起的地表沉降槽的变化形态,分析了因施工工序不同偏挖引起的地表变形的规律,拱顶沉降及拱顶主应力的变化规律。

研究表明:

（1）在断面形状尺寸,土层参数、地质条件相同的条件下,单洞双层隧道施工一次性开挖成洞,会造成大范围的地表沉降,甚至洞孔坍塌。

合理地分步开挖、及时施作衬砌有效地控制了地表变形和洞周位移。

（2）施工工序不同使隧道的偏挖引起地表沉降槽向未开挖一侧偏移,表明隧道在开挖过程中地层最大沉降并不发生在隧道中线处,而在完成全部开挖,地层变形稳定后,其累计沉降最大值位于隧道中线处。

（3）沉降槽的形状受隧道开挖断面形状与开挖深度影响,隧道中线右侧导洞③④的开挖形成的地表沉降槽表明第三、四步开挖对地表变形的影响均小于第一、二步开挖对地表变形的影响。

（4）隧道开挖衬砌完全施作后,地表沉降变化不大,但是后序开挖步引起的结构内力应予以重视,尤其对于中隔壁支撑及拱底衬砌的支护,要及时施作拱底二次衬砌。

参考文献

[1]孙钧.地下工程设计理论与实践[M].上海:

上海科学技术出版社,1996:

36-42（SunJun.TheDesignTheoryandPracticeonUndergroundEngineering[M].Shanghai:

ShanghaiScientificandTechnicalPublishers,1996:

36-42（inChinese））

[2]吴波,高波,漆泰岳,等.城市地铁区间隧道洞群开挖顺序优化分析[J].中国铁道科学,2003,24（5）:

23-28（WuBo,GaoBo,QiTaiyue,etal.Optimizationanalysisofexcavationsequenceoftunnelgroupsbetweenmetrostations[J].ChinaRailwayScience,2003,24（5）:

23-28（inChinese））

[3]施仲衡,张弥,王新杰,等.地下铁道设计与施工[M].西安:

陕西科学技术出版社,1997:

69-75（ShiZhongheng,ZhangMi,WangXinjie,etal.Designandconstructionofmetro[M].Xi’an:

ShanxiScientificandTechnicalPublishers,1997:

69-75（inChinese））

[4]刘波,韩彦辉（美国）.FLAC原理、实例与应用指南[M].北京:

人民交通出版社,2005:

12-20（LiuBo,HanYanhui（）.FLACTheory,ExampleandApplicationManual

[M].Beijing:

ChinaCommunicationsPress,2005:

12-20（inChinese））

[5]ItascaConsultingGroup,Inc.FLAC（FastLagrangianAnalysisofContinua）User’sManuals,[M].Minneapolis:

ItascaConsultingGroup,Inc.,2005

[6]DonovanK,PariseauWG,CepakM.Finiteelementap-proachtocableboltinginsteeplydippingVCRslopes[M]//GeomechanicsApplicationinUndergroundHardrockMin-ing.NewYork:

SocietyofMiningEngineers,1984:

65-90