盾构掘进引起土体扰动的数值模拟与现场验证_顾荣华.pdf
《盾构掘进引起土体扰动的数值模拟与现场验证_顾荣华.pdf》由会员分享,可在线阅读,更多相关《盾构掘进引起土体扰动的数值模拟与现场验证_顾荣华.pdf(10页珍藏版)》请在咨信网上搜索。
1、第5卷第1期 地 基 处 理 Vol.5 No.1 2023 年 1 月 Journal of Ground Improvement Jan.2023 盾构掘进引起土体扰动的数值模拟与现场验证 顾荣华1,沈卓恒2,王益坚3,高 飞1,边学成1*,陈云敏1(1.浙江大学 软弱土与环境土工教育部重点实验室/岩土工程研究所,浙江 杭州 310058;2.杭州市城市建设投资集团有限公司,浙江 杭州 310012;3.杭州市运河综合保护开发建设集团有限责任公司,浙江 杭州 310016)摘 要:为了研究结构性软土中盾构掘进过程对土体扰动的影响,引入了修正结构性剑桥模型考虑土的结构性,在 Abaqus 中
2、编写了该模型的子程序并建立了盾构掘进的三维有限元模型,以模型中结构强度参数为评价指标建立了盾构掘进引起土体扰动评价方法。以宁波地铁 1 号线二期工程邱隘东站到五乡站左线盾构隧道项目为例,分析了盾构推进过程中的土体变形、超孔压分布及土体扰动情况。结果表明,土体结构性变化对预测结果的准确性有较大影响,当考虑土体结构性的影响时,预测结果更接近实测值;数值模拟得到隧道底部土体扰动度大致呈三角形分布,最大值为 31%,影响范围达 3.0 m,与采用静力触探锥尖阻力扰动评价结果基本一致。对盾构施工控制参数进行分析,开挖面的支护压力对土体扰动的影响较大,当增大 0.5 倍支护压力时,隧道底面扰动度从 31%
3、增加至 40%,影响范围从 3.0 m 增大到 4.0 m。关键词:结构性软土;修正结构性剑桥模型;结构强度;盾构掘进;土体扰动度 中图分类号:TU443 文献标识码:A 文章编号:20967195(2023)01003310 Numerical simulation and field verification of the soil disturbance caused by shield tunneling in structural soft soil GU Rong-hua1,SHEN Zhuo-heng2,WANG Yi-jian3,GAO Fei1,BIAN Xue-cheng1
4、*,CHEN Yun-min1(1.Key Laboratory of Soft Soils and Geoenvironmental Engineering of the Ministry of Education/Institute of Geotechnical Engineering,Zhejiang University,Hangzhou 310058,Zhejiang,China;2.China Hangzhou Urban Construction Assets Management Co.,Ltd.,Hangzhou 310012,Zhejiang,China;3.Hangzh
5、ou Canal Comprehensive Protection and Development Construction Group Co.,Ltd.,Hangzhou 310016,Zhejiang,China)Abstract:In order to study the influence of shield tunneling on soil disturbance,the Modified Structured Cam Clay Model(MSCC)subroutine is compiled to establish the analysis model of shield t
6、unneling in structural soft soil established in Abaqus.A theoretical calculation method of soil disturbance degree is proposed based on structure strength in the constitutive model.Taking the excavation process of shield tunnel from Qiuai station to the Wuxiang station in the second stage of Ningbo
7、Metro Line 1 as an example,the soil deformation,excess pore water pressure and soil disturbance degree are analyzed.The results show that the accuracy of the predicted results can be improved by considering the soil structure.Numerical results shows that the disturbance degree of the soil at the bot
8、tom of the tunnel has roughly triangular distribution,with a maximum value of 31%and an influence range of 3.0 m,which is basically consistent with the evaluation results of the static cone tip resistance disturbance.By analyzing the control parameters of shield construction,the support pressure has
9、 a greater influence on the soil disturbance degree.When the support pressure is increased by 0.5 times,the soil disturbance degree under the tunnel increases from 31%to 40%and the influence range increases from 3.0 m to 4.0 m.Key words:soil structure;Modified Structured Cam Clay Model(MSCC);structu
10、re strength;shield tunneling;soil disturbance degree 0 引 言 土体在漫长的历史沉积过程中,土颗粒间通过化学、生物和力学等多种作用下形成某种特定的结构形式或在土颗粒间形成胶结作用力,导致原状土的力学特性与重塑土相比存在较大差异,表现出结 收稿日期:2022-03-10 基金项目:国家自然科学基金项目(52125803);浙江省省重点研发项目(2019C03111)。作者简介:顾荣华(1996),男,贵州安顺人,硕士研究生,主要从事隧道工程数值分析方面的研究工作。E-mail:。*通信作者:边学成(1976),男,浙江诸暨人,博士,教授,主要
11、从事交通岩土工程和地铁保护研究工作。E-mail:。DOI:10.3785/j.issn.2096-7195.2023.01.005 34 地 基 处 理 2023年1月 构性特征1-2:当土体未达到屈服状态时,土体中胶结结构的存在提高了土骨架的稳定性,与重塑土相比,在同等应力条件下结构性土可以维持更高的孔隙比;当土体达到屈服后,结构性的破坏使土体的压缩性急剧增大3-4。文献5 指出东南沿海地区广泛分布着海相沉积的软弱黏土层,通常具有显著的结构性特征。在这样的土层中进行地下结构施工时,会对土体产生严重的扰动,使土体结构性发生破坏,导致过大的土体变形,危及地下结构的安全6-7。随着沿海地区工程建
12、设规模不断增大,建立施工扰动程度评价方法,探究受扰动土体力学性质的变化规律具有重要的工程价值及实践意义。根据原位测试与现场监测结果,许多学者采用不同的试验指标建立了土体扰动度评价方法。XU等8 在研究盾构隧道掘进过程中发现有效应力的变化对土体力学性质有明显影响,提出了将有效应力作为指标评估土体扰动度。王军9 发现土体扰动后不排水强度大为降低,提出将十字板剪切强度来定量评价现场地基扰动度,并建立了土体变形参数与扰动度之间的关系。陈云敏等10 在对某基坑坍塌后进行现场测试,结果表明,锥尖阻力与十字板强度基本呈线性关系,因此将锥尖阻力作为土体扰动度的评价指标。鲁泰山等11 以不排水抗剪强度为评价指标
13、,建立了基于孔压静力触探锥尖阻力的扰动评价方法。准确评价土体扰动度对工程变形控制十分重要,但是基于原位测试扰动评价方法普遍存在时效性差、监测困难、监测精度低等问题。因此,发展一种基于有限元模型的施工扰动评价方法十分必要。有限元模型计算的精度取决于所选的材料本构模型及参数,目前已经发展出很多能够考虑土体结构性的本构模型。CHAI 等12 将结构性土的压缩曲线方程引入到修正剑桥模型中,建立了结构性参数模型;LIU等13 和CARTER等14 在结构性土压缩曲线的基础上,提出了应力状态、应力历史、孔隙比和结构性指数的四维空间理论,以此建立了结构性本构模型。结构性软土中普遍存在胶结作用15-17,影响
14、着土的力学特性。为了考虑胶结作用力的影响,SUEBSUK 等18 将屈服面平移至与平均主应力轴负半轴相交,通过假定截距随塑性偏应变发展来考虑胶结作用的演化,提出了修正结构性剑桥模型(MSCC),并通过对比Osaka 黏土19 与 Marl 黏土20 模拟结果与试验结果,验证了该模型在不同地区的结构性软黏土的适用性。尽管结构性土的本构模型研究取得了很大的进展,但是目前还未运用到常见的有限元软件中。为此,本文首先基于 Abaqus 提供的二次开发平台编写了MSCC 模型的子程序,并以宁波地铁 1 号线二期工程邱隘东站到五乡站左线盾构隧道掘进过程为例,利用MSCC 模型在 Abaqus 中建立三维有
15、限元模型,以评估土体结构性对隧道周围土体响应的影响。最后,基于 MSCC 模型中结构强度参数提出土体扰动度计算方法,并通过与现场实测结果对比以验证该方法的合理性,分析了盾构施工控制参数对土体扰动的影响。研究结果能够为类似的工程设计施工提供参考。1 修正结构性剑桥模型 1.1 模型介绍 ROSCOE 等21 基于能量方程建立了修正剑桥模型(MCC),模型相对简单、参数易于确定,可以准确地描述重塑土的力学特性。但工程中的实际土具有明显的结构性,其强度变形特性与结构性密切相关,不同结构性土的临界状态特性亦不同。因此,用修正剑桥模型来描述结构性软土的力学反映显然存在着较大的误差。为此,SUEBSUK
16、等18 在建立定量描述土结构性的基础上,考虑土的胶结作用,提出了修正结构性剑桥模型(MSCC)。模型屈服面见图 1,屈服函数F为:22b0()()FqMpppp=-+-(1)式中:p为平均主应力;q为广义剪应力;bp为考虑胶结强度而引入的参数,称为结构强度;M为临界状态线斜率。bp0pqdp0结构性土屈服面重塑土屈服面M未扰动临界状态线完全扰动临界状态线p 图 1 结构性土屈服面和重塑土屈服面 Fig.1 Yield surface of structural soil and remolded soil HORPIBULSUK 等15 和 LIU 等16 指出结构强度bp的变化与塑性应变密切
17、相关,模型中假定结构强度bp取决于塑性偏应变pd的大小,并将其变化分为两个阶段来表示胶结结构的退化和破碎,计算公式为:结构退化阶段:()pbb0dexppp=-(2)结构破碎阶段:()ppbb,fdd,fexppp=-(3)式中:b0p为初始结构强度;pd,f为土体抗剪强度达第 1 期 顾荣华,等:盾构掘进引起土体扰动的数值模拟与现场验证 35 到峰值时对应的塑性偏应变;b,fp为pd,f对应的结构强度;为结构强度衰减指数。模型中,将重塑土作为参考状态来描述原状土的压缩特性,见图 2。当固结压力p小于初始屈服应力y,ip时,土体处于弹性状态,原状土孔隙比可根据固结压力p得到;当固结压力p大于初
18、始屈服应力y,ip时,土体发生弹塑性变形,假设原状土与重塑土之间的孔隙比差值e随应力p按指数b衰减,原状土孔隙比e的计算公式见式(4)。()*iy,i*IC/lnbeeeeeppeep=+=-(4)式中:e*为固结压力为p的重塑土的孔隙比;ie为初始屈服应力y,ip对应的孔隙比差值;b为孔隙比衰减指数;*ICe为重塑土初始孔隙比,表示重塑土在固结压力1 kPap=情况下的孔隙比;为lnep-坐标中重塑土压缩曲线斜率。y,ip*=+e ee*ee0peie*ICe1重塑土等向固结线ln p原状土等向固结线孔隙比 e 图 2 结构性土压缩特性 Fig.2 Compression behavior
19、of structural soil 与修正剑桥模型(MCC)相同,将土体理想化为各向同性材料,屈服面的大小只取决于塑性体应变,即将塑性体应变作为硬化参数,根据式(4)推导出该模型的硬化规律为:()()0p0v00d(1)dd(1)*pMb eMMe ppb eMe p-+Er(误差控制值,取 Er=105),当前增量步不收敛,则减小子步长度,取 Tnew=0.8(Er/R)1/2回到第(2)步重新计算,如果如果 REr,则移到下一步。2121():()12():()R-=+(7)(5)子步计算收敛,局部误差控制成功,更新应力分量+、硬化参数0,10,2(+)/2pp。(6)令TTT=+,新子
20、步增量同样采用Tnew=0.8(Er/R)1/2T 确定。如果new+1TT,取36 地 基 处 理 2023年1月 new1TT=-,回到第(2)步重新计算,当1T=时循环终止。(7)根据最终计算应力+计算屈服函数值F,参照文献27,当610F-时进行应力漂移修正,并更新弹塑性矩阵。(8)保存应力、应变、硬化参数及其他计算需要的量,将刚度矩阵赋值给雅克比矩阵,本次增量步结束。2 盾构掘进数值模拟 2.1 工程背景 宁波地铁 1 号线二期工程西起一期终点东环南路站,终点至霞浦站。线路全长 25.3 km,是宁波市轨道交通线网规划中的东西走向主干线。其中,邱隘东站到五乡站区间全长为 3.7 km
21、,左线地下段长 825.8 m,采用盾构法施工。隧道覆土厚度在 3.06.9 m 之间,掘进地层以软黏土地层为主。隧道外径为 6.2 m,内径为 5.5 m,管片环宽 1.2 m,开挖直径为 6.39 m,隧道中心埋深为 6.810.1 m。宁波地区软黏土具有显著的结构性,受施工扰动影响,其强度会降低、压缩性急剧增大。为了研究结构性软土地层中盾构掘进过程中土体的变形、孔压响应及扰动情况,在工程现场开展了试验测试,试验场地位于邱隘东站至五乡站区间,试样现场位置见图3。图 3 试验现场位置示意图 Fig.3 Schematic diagram of test site location 根据地质勘
22、察资料,试验场地隧道中心埋深为8.7 m,土层分布情况及土层物理力学指标见表 1。从表中可以看到,土层 1、土层 2 灵敏度都大于 4,属于高灵敏度土,结构性较强。2.2 有限元模型 利用数值模拟软件 Abaqus 对盾构隧道开挖引起的地层响应进行三维弹塑性有限元分析,模型尺寸需满足圣维南原理的要求,即隧道与整体模型的边界距离大于 35 倍的隧道直径,文献28 表明在数值模拟中,隧道长度至少要 120 m 才能充分消除边界效应。考虑模型的对称性,取地层几何模型的尺寸在 X、Y、Z 方向分别为 120 m、50 m、50 m,开挖面在 YZ 面,开挖方向沿 X 负轴方向开挖,且模型底部边界固定,
23、侧边竖向位移自由,顶面没有任何约束,开挖模型的建立及网格的划分见图 4。图 4 模型建立及网格划分 Fig.4 Modeling and mesh 有限元模型中,土体采用可考虑孔压的 C3D8P实体单元进行模拟,并假设土体完全饱和且渗透系数在相同土层中各向同性,根据地质勘察资料确定各土层渗透系数见表 1。土层 1、2 采用 MSCC 模型模拟,土层 3 灵敏度低、结构性较弱,采用修正剑桥模型(MCC)模拟,各土层的本构模型参数见表 2。为了反映不考虑土的结构性对模拟结果的影响,土层 1、2 同时采用 MCC 本构模型模拟作为参照模型,材料参数与 MSCC 模型相同,其中参数、b、ei、Pb0、
24、Pv,i取 0,取=2,土体为正常固结土。衬砌、注浆层采用 C3D8I 实体单元进行模拟。盾构机盾壳长 8.65 m,模型中将其简化为长度 8.4 m(7 个管片环宽度)、外径 6.39 m 的柱体,采用 S4R壳单元模拟。衬砌、注浆层、盾构机均采用弹性材料,材料参数取值见表 3。施工过程的模拟参照文献30,采用单元生死法模拟盾构隧道的掘进过程。表 1 土层分布及其参数 Table 1 Soil strata and soil parameters 土层编号 土层名称 层厚 相对密度/Gs 饱和重度 sat/(kN m3)含水率 孔隙比 e0 土压力 系数 k0 渗透系数 kv/(m s1)灵
25、敏度 St H/m w/%土层 1 淤泥质粉质黏土 10 2.73 17.40 48.23 1.34 0.58 6.0 10-9 4.9 土层 2 淤泥质黏土 10 2.74 17.67 44.47 1.27 0.52 2.5 10-9 7.3 土层 3 黏土 30 2.73 18.80 34.60 1.00 0.47 8.1 10-9 1.8 掌子面 盾壳 掘进方向 衬砌与注浆层 埋深 5.5 m 淤泥质粉质黏土 淤泥质黏土 黏土 50 m 120 m 50 m 50 m 10 m 10 m Z Y X 邱隘东站 地下段,全长 825.8 m 试验现场 宁波地铁 1 号线 五乡站 第 1 期
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 盾构 掘进 引起 扰动 数值 模拟 现场 验证 荣华
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。