基于CVISC模型的深埋圆硐黏弹塑性解.pdf
《基于CVISC模型的深埋圆硐黏弹塑性解.pdf》由会员分享,可在线阅读,更多相关《基于CVISC模型的深埋圆硐黏弹塑性解.pdf(9页珍藏版)》请在咨信网上搜索。
1、引用格式:鲁义强,史博然,贺飞,等.基于 CVISC 模型的深埋圆硐黏弹塑性解J.隧道建设(中英文),2023,43(增刊 1):145.LU Yiqiang,SHI Boran,HE Fei,et al.Viscous-elastic-plastic solution of deep-buried circular opening based on CVISC modelJ.Tunnel Construction,2023,43(S1):145.收稿日期:2022-10-05;修回日期:2023-05-20基金项目:国家自然科学基金(52204108);中国中铁股份重大专项(cz02-专项-
2、03)第一作者简介:鲁义强(1989),男,河南开封人,2020 年毕业于四川大学,土木工程专业,博士(后),高级工程师,主要从事 TBM 设计和隧道防灾工作。E-mail:yiqianglu007 。通信作者:袁勇,E-mail:yuany 。基于 CVISC 模型的深埋圆硐黏弹塑性解鲁义强1,2,史博然2,贺 飞1,袁 勇2,姚旭朋2,张姣龙2(1.中铁工程装备集团,河南 郑州 450000;2.同济大学土木工程学院,上海 200092)摘要:为解决高地应力围岩挤压型大变形问题,根据以往挤压型大变形隧道案例,选取适用于泥岩、砂岩、粉质砂岩等各向同性软岩,利用 Burgers 模型和 M-C
3、 塑性体串联形成的 CVISC 黏弹塑性模型,结合 CVISC 模型黏弹性元件和弹塑性元件相互独立的特性,将围岩径向位移分为瞬时弹塑性阶段位移和黏弹性阶段位移,分别推导塑性区和黏弹性区径向位移解析解,叠加组合形成围岩位移分布解析解,分析围岩径向位移场、应力场随时间变化规律,全面描述硐室变形规律。借助 FLAC3D 和相关文献理论解答验证其适用性,在此基础上进行长期变形的参数敏感性分析,结果表明各阶段最关键的 3 个参数分别是内摩擦角、黏聚力和弹性模量,蠕变参数仅在其控制阶段具有显著敏感度。关键词:岩石隧道;挤压型大变形;黏弹塑性解析;有限差分法DOI:10.3973/j.issn.2096-4
4、498.2023.S1.017中图分类号:U 45 文献标志码:A 文章编号:2096-4498(2023)S1-0145-09V Vi is sc co ou us s-E El la as st ti ic c-P Pl la as st ti ic c S So ol lu ut ti io on n o of f D De ee ep p-B Bu ur ri ie ed d C Ci ir rc cu ul la ar r O Op pe en ni in ng g B Ba as se ed d o on n C CV VI IS SC C MMo od de el lLU Yiq
5、iang1,2,SHI Boran2,HE Fei1,YUAN Yong2,*,YAO Xupeng2,ZHANG Jiaolong2(1.China Railway Engineering Equipment Group Co.,Ltd.,Zhengzhou 450000,Henan,China;2.College of Civil Engineering,Tongji University,Shanghai 200092,China)A Ab bs st tr ra ac ct t:To control the large squeezing deformation of surround
6、ing rocks with high geostress,the CVISC viscoelastic-plastic model formed by the Burgers model and M-C plastic body in series is selected for isotropic soft rocks such as mudstone,sandstone,and siltstone based on the previous tunnel cases with large squeezing deformation.The radial displacement of t
7、he surrounding rock is divided into instantaneous elastic-plastic phase displacement and viscoelastic phase displacement using the independent properties of viscoelastic and elastic-plastic elements of the CVISC model.The analytical solutions of radial displacement in the plastic and viscoelastic zo
8、nes are derived respectively,and the superimposed combination forms the analytical solution of displacement distribution in the surrounding rock.Following this,the radial displacement and stress fields of the surrounding rock are analyzed as a function of time,and the deformation law of the opening
9、is described.The applicability of model is verified with the help of FLAC3D and the theoretical solution of related literature.On this basis,parametric sensitivity analysis of long-term deformation is carried out.The results show that the three most critical parameters in each phase are internal fri
10、ction angle,cohesion,and Youngs modulus,and the creep parameter has significant sensitivity only in its control phase.K Ke ey yw wo or rd ds s:rock tunnels;large squeezing deformation;viscous-elastic-plastic analysis;finite difference method0 引言 高地应力围岩挤压型大变形贯穿隧道开挖、支护乃至运营的各个阶段,是一个随时间发展的动态、长期过程,易诱发支护结
11、构屈曲失稳和 TBM 盾壳卡机,威隧道建设(中英文)第 43 卷胁工程安全。深入研究隧道围岩挤压型大变形,应当以围岩的流变效应为基础,从隧道长期变形机制入手。随着岩石流变力学的不断发展,各国学者针对不同的流变力学行为提出了各类流变本构模型,基于黏弹性1、黏弹-塑性2-3、黏弹-黏塑性4-7流变模型的圆形硐室变形解析公式更是层出不穷。黏弹性解析解方面,赵旭峰8基于 Burgers 黏弹性平面模型,推导了考虑开挖面荷载释放的圆硐变形解析解。Wang 等9借助复变函数推导了考虑开挖边界状态随时间变化的解析解。卞跃威等10基于Burgers 黏弹性模型,推导了非静水应力场考虑实际施工过程应力释放的圆形
12、隧道黏弹性解析解。储昭飞等11基于黏弹性对应性原理,推导了隧道围岩-衬砌结构光滑接触和完全接触 2 种情况下的黏弹性通解。黏弹-塑性解析解方面,现有研究通常假定围岩的流变特性符合黏弹性方程,再引入强度准则。张良辉等12给出了一种既能够反映围岩应变软化和塑性剪胀特性,又能体现围岩流变特性的黏弹-塑性解析公式。卞跃威等13在考虑施工过程应力释放的圆硐黏弹塑性解的基础上,提出了考虑塑性软化、塑性体积膨胀和围岩应力释放的圆形隧道弹塑性解。夏才初等14推导了不考虑塑性区围岩的剪胀及峰后应变软化,弹、塑性区围岩的应变及位移表达式相同假设下,Kelvin-Voigt 黏弹-塑性模型的圆硐解析解。黏弹-黏塑性
13、解析解主要以西原模型为主,曹瑞琅等15建立了考虑支护阻力和掌子面推进对围岩应力释放影响的基于 Hoek-Brown 强度准则的西原模型圆形隧洞黏弹塑性解。夏才初等7假设黏塑性体的偏应变张量的一阶导数与瞬态偏应力张量和稳态偏应力张量之差成正比,基于 Laplace 正、逆变换,推导了基于西元模型的黏弹-黏塑性解析解。王晓飞16基于广义 Bingham 的弹黏塑性模型,推导了基于 Mohr-Coulomb 和 Hoek-Brown 强度准则的圆硐解析解。上述黏弹性解析解物理意义明确但未考虑塑性区对位移场分布的影响;现有的黏弹-塑性解析解主要建立在 Kelvin-Voigt 黏弹-塑性模型基础上,该
14、模型能够反映岩体瞬时弹塑性性态,但蠕变变形最终收敛,变形速率趋向于 0,难以反映强流变条件下围岩的长期变形;黏弹-黏塑性解析解能够综合考虑围岩的黏弹性和黏塑性特征,但其黏塑性区半径将随时间变化,系变边界问题,应力场求解复杂。此外,相关模型元件复杂,实际应用时参数获取较为困难。因此,对深埋硐室围岩变形的黏弹塑性解的探索仍然是学术界关注的热点。从室内试验结果来看17,瞬时弹塑性、黏弹性和黏性分别反映岩石瞬时性态、衰减蠕变和定常蠕变,能够准确刻画挤压型围岩(泥岩、砂岩、粉质砂岩等各向同性软岩)的普遍力学特性,其对应元件 Hook 体、Saint-Venant 体、Kelvin 体和 Newton 体
15、串联构成的CVISC 流变模型物理意义明确,参数识别简便。赵旭峰8借助反分析程式,基于简单模型辨识理论,对比了 5 种流变模型,认为 CVISC 模型对隧道现场挤压性围岩流变时效特性的模拟效果最佳。Fuente 等18、Lee 等19、Sharifzadeh 等20均选取 CVISC 模型并对隧道位移展开反分析,以获取围岩参数。但目前基于CVISC 黏弹塑性模型解析解尚不完善。本文对静水应力场基于 CVISC 模型的圆形硐室解析解进行推导。推导静水应力场深埋圆形硐室围岩长期变形公式,并利用 FLAC3D 有限差分软件和相关文献理论解答验证其适用性,在此基础上进行长期变形的参数敏感性分析,为后续
16、围岩长期变形的研究提供理论支撑。1 围岩变形的基本方程1.1 基本假定 针对深埋硐室围岩变形问题,进行如下假设:1)研究对象为处于连续、均质、各向同性岩体的深埋圆形硐室,岩体侧压力系数为 1(即静水压力场),计算时可简化为轴对称的厚壁筒模型,并忽略硐室影响范围(5 倍硐径)内的岩体自重17;2)硐室纵向尺度远大于硐径尺度,将计算简化为平面应变问题,假定硐室每一横断面力学响应相同;3)围岩流变特性,采用 CVISC 黏弹-塑性模型。1.2 CVISC 黏弹-塑性模型基本方程 CVISC 黏弹塑性流变模型由 Burgers 体和 M-C 塑性体 2 部分串联组成,其中 M-C 塑性体由一个表征Mo
17、hr-Coulomb 强度准则的开关和一个表征理想刚塑性流动法则的塑性体并联组成。对于软岩隧道工程,在反复开挖扰动和流变效应双重影响下,围岩的应变速率包含了 Burgers 模型的黏弹性应变速率和 M-C 塑性体的刚塑性应变速率,CVISC 模型可以真实地反映挤压型地层隧道开挖过程中所表现的非线性力学行为21。CVISC 元件模型示意如图 1 所示。(a)黏弹性力学行为(b)黏弹塑性力学行为图 1 CVISC 模型示意图Fig.1 Schematic of CVISC model641增刊 1鲁义强,等:基于 CVISC 模型的深埋圆硐黏弹塑性解 为元件所受应力,s为屈服极限,当 0 时,由于
18、外荷载恒定,不再发生弹塑性变形,弹塑性元件可视作为刚体,仅由黏弹性元件发生黏弹性变形。因此,根据式(36)、(37),弹性区弹性位移和黏弹性位移,塑性区弹塑性位移和黏弹性位移能够叠加,一方面变形协调使然,另一方面是将不同时段变形的叠加,本质是时间的叠加。塑性区总位移用式(38)表示。pr=dupdr=sprJ(t)=sprMt+sprGK1-exp-GKKt();(32)spr=pr-13pr+p+pz()=13A 2-1()+pi+A()rri()Kp-1()2-Kp-Kp-()-1-2()p0;(33)dup=sprJ(t)dr=sprMt+sprGK1-exp-GKKt()dr;(34)
19、up=rJ(t)3pi+A()2-Kp-Kp-()1KprKp-1()irKp-1()-rKp-1()p()+3 p0-rp();(35)up=upe+up;(36)ue=uee+ue;(37)uit,pi()=12GD1+D2ri-rpri()Kps-2Gu(rp)+D1rpri()Kp+D2rp+riJ(t)32-Kp-Kp-()pi-rp()1Kp+3 p0-rp()。(38)式中:dupv为塑性区黏弹性位移增量;srp为塑性区径向应力偏量;upe为塑性区弹性位移;uee为弹性区弹性位移。值得注意的是,硐内支护力 pi、原岩应力 p0以及围岩材料参数一定时,对于不同设计硐径 ri,塑性区
20、的相对半径 rp/ri保持恒定,硐壁处相对收敛 ui(t,pi)/ri相同,围岩的应力场和位移场分布相同,尺寸效应不会影响围岩的相对应力和相对位移。因此,任意硐径硐室的相对收敛相同,硐室绝对收敛与设计硐径成正比。3 解析解的验证 利用 FLAC3D 有限差分软件和文献14推导公式的衍生式,对解析解进行验证,证明公式的可靠性和适用性。3.1 参数选取 数值模型涉及弹性参数、强度参数、蠕变参数和隧道设计参数。文献24-25采用单轴压缩试验、三轴压缩试验和三轴蠕变试验对泥岩的力学特性进行了深入研究,本节采用其参数进行分析,如表 1 所示。3.2 数值模型的建立 考虑力学对称性,采用 1/4 模型进行
21、分析,设计硐径 1 m,根据圣维南理论,设置模型边长 10 m,厚度1 m 的平面应变模型。边界条件方面,分别约束 x、z、y方向位移;设置右边界和上边界为应力边界条件,分别施加 x 轴负方向和 z 轴负方向 8.94 MPa 压力(强度应力比 0.5),以模拟原岩应力 p0。模型示意如图 4所示。表 1 数值模型参数Table 1 Numerical model parameters参数类型参数参数值弹性参数弹性模量 E/GPa 0.650泊松比 0.130体积模量 K/GPa 0.293剪切模量 G/GPa 0.288强度参数单轴抗压强度 c/MPa 4.470黏聚力 c/MPa 2.06
22、0内摩擦角/()10.000蠕变参数Maxwell 体剪切模量 GM/MPa306.000Maxwell 体黏滞系数 M/(MPah)1.60105Kelvin 体剪切模量 GK/MPa93.000Kelvin 体黏滞系数 K/(MPah)3.73103隧 道 设计参数设计硐径 ri/m 1.00原岩应力 p0/MPa 8.943.3 数值解与解析解之间的比较 本节分为瞬时弹塑性模型和时程黏弹塑性模型 2部分进行验证,其中相对误差计算以数值解为分母。文献14公式适用于 Kelvin-Voigt 黏弹-塑性模型,本文推导其适用于 CVISC 模型的衍生式,考虑到衍生式瞬时弹塑性力学机制与本文解析
23、解一致,仅对黏弹塑性模型进行验证。3.3.1 瞬时弹塑性模型验证 瞬时弹塑性模型围岩径向和切向应力的数值解和解析解随相对半径增加的变化曲线如图 5 所示。横轴和纵轴都进行了归一化处理。941隧道建设(中英文)第 43 卷图 4 数值模型示意图(单位:MPa)Fig.4 Schematic of a numerical model(unit:MPa)图 5 不同相对半径处围岩径向和切向应力比较Fig.5 Comparison of radial and tangential stresses of surrounding rock at different relative radii塑性区半径
24、方面,解析解和数值解的塑性区半径分别为 2.46 m 和 2.47 m,相对误差为 0.4%,误差可以接受。塑性区切向应力和径向应力的最大相对误差分别为 0.75%和 4.77%,随着相对半径增加,位置靠近塑性区半径,相对误差单调递减。弹性区切向应力和径向应力的最大相对误差分别为 2.37%和 1.77%,分别出现于相对半径 10.0 和 8.3位置处,即模型应力边界处。应力相对误差均小于5%,在可接受的范围内;应力误差分布方面,弹塑性区的最大相对误差主要分布在模型两端的应力边界处,对实际分析影响微弱。围岩径向相对位移的数值解和解析解随相对半径增加的变化曲线如图 6 所示。图 6 不同相对半径
25、处围岩径向相对位移比较Fig.6 Comparison of radial relative displacement of surrounding rock at different relative radii塑性区径向相对位移的最大相对误差为 3.69%,出现于相对半径 1.0 位置处,即硐壁位置处,随着相对半径增加,相对误差单调减小,并在相对半径 1.9 处达到最小,随后继续单调递增。塑性区半径处相对误差为 2.89%。弹性区径向相对位移的相对误差随相对半径增加而单调递增,最大相对误差为 12.2%,出现于相对半径 10.0 位置处,即应力边界位置。塑性区径向相对位移的相对误差均小于
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 CVISC 模型 深埋圆硐黏弹 塑性
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。