摘 要:根据现场地质条件,以地基关键断面为研究对象建立相应的数值模型,应用自行开发RFPA2D程序就广州地铁二号线盾构开挖对新南方购物中心大楼稳定性影响问题进行了数值模拟研究·分析了隧道断面洞顶的沉降位移和建筑物地基的沉降位移,数值模拟结果再现了隧洞开挖引起的周围岩体的损伤破坏过程,得到了隧洞及地基周围岩土的应力变化图,塑性区破坏图,以及位移矢量图,预报了破碎带的未来发展趋势·数值计算结果与采用解析算法的结果进行了对比验证,结果与现场观测资料基本吻合·

　　由于城市地铁开挖引起的相应的问题(比如说地基沉降、地表沉陷等)近年来在实际工程中显得尤为关注,特别是地铁开挖引起地基沉降的位移变化情况、应力分布状态以及隧洞的稳定性是设计者和施工单位所关心的·为了限制岩土体的变形破坏,必须选用适当的开挖方法和相应的加固措施·正确预测开挖后地铁周围岩土体的位移变化、应力矢量分布以及隧洞的变形破坏和破坏范围显得十分重要·但由于工程开挖问题涉及到的地质因素、工程因素复杂[1,2],本文提出采用自行开发的RFPA2D数值模拟技术和工程分析法对隧洞盾构施工过程引起浅基础下土层的应力应变进行分析·

　　工程实例以广州地铁盾构开挖对新南方购物中心大楼稳定性影响的分析为例·

　　1 工程概况

　　广州地铁二号线客村—鹭江站区间是盾构法的施工区间·该区间隧洞穿越六层的新南方购物中心时,有一排柱的基础位于开挖隧洞的侧上方,垂直距离是7.79~8.69m·按设计单位提供的地质资料、地铁断面宽度为10m,最高拱高为6m,平均埋深为11m,围岩主要为砂岩,具体地质及隧洞相关位置如图1所示,地层由上至下分布顺次为:杂填土(1)、粉质粘土(2)、硬塑,标贯击数N为16~36,强风化泥质砂岩(3)、中风化泥质砂岩(4)、微风化泥质砂岩(5)的层状分布·

　　1.1 开挖方式

　　地铁开挖方式采用盾构法,土压平衡式盾构机内径5.4m,衬砌砼管片厚0.3m,盾构外径6.14m,纵向长度约7.8m,依靠安装在已有足够强度管片的周圈上10多个千斤顶的水平推力掘进,最大水平总推力约31000kN·当盾构向前掘进时,会对地面、建筑物、桩基产生影响·在盾构机的后部,需对已安装好1.5m长的圆环混凝土管片进行注浆,此时,圆环管片与周围的岩土体存在15cm的间隙,在土体及基底荷载作用下,基础将产生沉降·

　　2 数值模拟

　　2.1 岩石破坏过程分析系统(RFPA2D)简介

　　本文所用RFPA2D系统[3~12]是一个基于有限元应力分析原理以及弹性损伤理论开发能模拟岩石介质变形与破坏的数值模拟软件,当残余强度取值很大时,同样可以模拟理想弹塑性和脆塑性岩土介质·与其他已有的逐渐破坏模型(ProgressiveFailure/Fracture/DamageMode)一样,RFPA2D包括应力分析和破坏分析·应力分析求解器采用有限元进行,依据有效应力原理完成每一计算步在外载荷或环境因素作用下(加载,开挖和位移边界条件的改变等)的力学响应·逐渐破坏模型是在假定组成材料各单元的力学性质整体上服从Weibull分布的基础上,根据一定的破坏准则检查材料中是否有单元破坏,对破坏单元则采用刚度特性退化(分离处理)和刚度重建(接触处理)的办法进行处理·

　　考虑到岩土类材料的抗拉强度,因此本文采用了修正的库仑(Coulomb)准则(包含拉伸断裂Tensioncut-off[6])作为单元破坏的强度判据·上述分析过程可以处理裂纹萌生、扩展以及裂纹的非规则路径选择问题·

　　2.2 数值模型

　　取垂直于隧洞的横剖面作为研究对象,按照平面应变问题考虑·计算取20m×40m的分析区域,上部取到地表,左右两侧的水平位移取零,垂直方向可自由滑动,底部水平和垂直位移取零·分析在盾构施工过程中,管片衬砌与周围岩土体之间150mm缝隙引起的地基应力和位移变化·计算时,为了将实际的空间问题简化成平面问题计算,将柱底的荷载N除以3,并乘1.3的增大系数,而隧道仅在1.5m长有15cm间隙变成无限长的间隙·

　　2.3 计算参数

　　针对该工程的具体问

[1] [2] 下一页