Processing math: 100%
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus 收录
  • 全国中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

基于超前应力释放与注浆加固的隧道围岩大变形控制分析

王俊杰 张帅

赵春发, 李煜寒, 彭也也, 杨晶, 宁晓芳, 冯洋. 超导电动悬浮列车明线气动特性及其对悬浮状态的影响[J]. 西南交通大学学报. doi: 10.3969/j.issn.0258-2724.20240470
引用本文: 王俊杰, 张帅. 基于超前应力释放与注浆加固的隧道围岩大变形控制分析[J]. 西南交通大学学报. doi: 10.3969/j.issn.0258-2724.20230706
ZHAO Chunfa, LI Yuhan, Peng Yeye, YANG Jing, NING Xiaofang, FENG Yang. Aerodynamic Characteristics of Open Wire of Superconducting Maglev Train and Its Influence on Levitation State[J]. Journal of Southwest Jiaotong University. doi: 10.3969/j.issn.0258-2724.20240470
Citation: WANG Junjie, ZHANG Shuai. Large Deformation Control of Tunnel Surrounding Rock Based on Advance Stress Release and Grouting Reinforcement[J]. Journal of Southwest Jiaotong University. doi: 10.3969/j.issn.0258-2724.20230706

基于超前应力释放与注浆加固的隧道围岩大变形控制分析

doi: 10.3969/j.issn.0258-2724.20230706
基金项目: 重庆市建设科技计划项目(城科字2022第8-12号)
详细信息
    作者简介:

    王俊杰(1977—),男,教授,博士,研究方向为岩土与地下工程,E-mail:wangjunjie@cqjtu.edu.cn

  • 中图分类号: U451.2

Large Deformation Control of Tunnel Surrounding Rock Based on Advance Stress Release and Grouting Reinforcement

  • 摘要:

    为克服超前导洞用于软岩隧道围岩大变形控制时存在的局限性,深入分析软岩大变形特征及其存在的问题,提出基于“超前应力释放 + 环向(滞后)注浆 + 加长锚杆”的变形控制方式. 首先,借助软岩峰后刚度与强度统一劣化模型与统一强度准则得到超前导洞及正洞围岩的弹塑性解;然后,利用FLAC3D有限差分软件实现软岩峰后刚度与强度统一劣化模型的本构开发,得到超前导洞及正洞围岩的变形及应力分布;最后,对软化模量、注浆参数、超前导洞半径和两掌子面间距等影响因素进行分析. 研究结果表明:超前导洞可以有效释放围岩挤压变形,岩体松动破碎是造成释放层及正洞围岩变形过大、稳定性下降的主要原因,环形(滞后)注浆可以有效控制扩挖过程中的围岩松动变形,并改善围岩应力分布,提高围岩承载能力;软化模量取值越大时超前导洞围岩变形越大,注浆参数取值越大时正洞围岩变形越小,增大超前导洞开挖半径、增加两掌子面间距(超前导洞及正洞)均可使初始地应力释放更加充分.

     

  • 超导电动悬浮列车(EDS)以其低能耗、低噪声、高速度及自稳定性等优点,近年来获得了显著发展. 磁悬浮列车的设计速度已达到600 km/h,如何实现高速稳定运行已成为磁悬浮列车空气动力学研究的关键方向之一[1-2]. 目前,日本已于山梨试验线完成超导电动悬浮列车最高603 km/h试运行,证明超导电动悬浮的有效性[3-4]. 磁浮列车行车速度增加至600 km/h恶化了高速磁浮列车的气动特性,严重影响高速磁浮列车运行平稳性及乘坐舒适性[5].

    国内外学者对磁浮列车气动特性已开展较多研究. 针对常导高速磁浮列车:丁叁叁等[6-7]指出磁浮列车底部流量直接影响磁浮列车气动升力,并提出在头尾车鼻尖增加导流板用于控制间隙流量从而控制磁浮列车气动升力;Zhou等[8]采用IDDES(improved delayed detached eddy simulation)方法计算了磁浮列车明线气动特性,给出了轨道附近列车风安全限值;孟石等[9]研究了轨道间隙对磁浮列车气动特性影响,发现轨道间隙越大,头车升力减小,尾车升力增大;此外还有学者探究流线型鼻尖对磁浮列车气动噪声的影响[10]. 针对钉扎磁浮列车:李一凡等[11]基于钉扎磁浮列车分析了底部流动对磁浮列车升力特性的影响,指出底部间隙流量减小,升力及点头力矩均减小;王潇飞等[12]则探究钉扎磁浮列车轨道结构对磁浮列车气动特性的影响. 而日本自20世纪70年代开始研究超导电动磁浮列车,山梨试验线测试结果表明,磁浮列车流线型长度增加,车体压力变化更为平缓,升力则一定程度减小[13]. Yamamoto等[14]对比了不同的列车截面发现,底部近似方形的列车截面可有效减小气动湍流,同时指出不同车体截面气动特性差异较大. Lin等[15]研究了超导电动悬浮列车在不同悬浮间隙下的气动特性,发现悬浮间隙越大,阻力变化较小,升力则减小,且低频涡流对磁浮列车升力和阻力有显著影响.

    在运行过程中,磁浮列车外流场气动特性直接反映于车体受到气动荷载. 列车起浮运行后仅受到电磁悬浮推进力及气动荷载的作用,为实现高速稳定运行,探究气动荷载作用特征的重要性逐步提升[16]. 相关研究表明,磁浮列车受到的气动力与运行速度的平方成正比,可以预见的是提速至600 km/h下气动荷载也将显著增加[17]. 刘堂红等[18]计算了磁浮列车以相对860 km/h交会时的气动特性,同时将气动荷载作用于包含车体、悬浮架详细车辆动力学模型中,结果表明,磁浮列车车体产生了较大的横向振动. 南凯威等[19]计算了TR08型磁浮列车明线及交会工况下气动荷载时程特征,结果表明,当磁浮列车速度增大到600 km/h时,悬浮间隙将超过安全极限. Wu等[20]建立了磁浮列车垂向动力学模型,通过特征值方法和直接积分法研究气动荷载对悬架稳定性的影响规律.

    上述国内外文献表明,磁浮列车与轨道之间的间隙流动直接影响磁浮列车气动特性. 但现有研究大多基于常导高速磁浮列车模型开展,对于超导电动悬浮列车的研究较少. 受限于超导电动悬浮的结构不同,流动特性存在一定差异[21]. 本文以某型超导电动悬浮列车为研究对象,采用CFD(computational fluid dynamics)数值计算方法对三编组磁浮列车明线工况气动特性进行分析,并将悬浮架视为独立部件,分析其附近的详细空气流动在此基础上通过SIMPACK和Simulink建立车辆动力学模型,并使用精细化电磁力计算方法计算电磁力,分析不同气动荷载加载方式下磁浮列车动力学响应差异.

    本研究选用的某型超导电动悬浮列车模型如图1所示,由头车、中间车和尾车三编组组成,头车与尾车外形一致,保留风挡特征. 磁浮列车车体与悬浮架非刚性连接,其简化模型如图1(c)所示. 使用磁浮列车车体高度H作为归一化长度,H为3.35 m. 车体总长为24.87H,建立简化超导电动悬浮列车悬浮架外形,悬浮架宽为0.95H,车体宽0.89H,头尾车流线型长度为4.24H,磁浮列车稳定运行悬浮高度为100 mm. 磁浮列车车体及悬浮架底部均简化为平面,不考虑车轮等部件.

    图  1  磁浮列车几何模型
    Figure  1.  Geometric model of maglev train

    本文主要研究磁浮列车明线稳定运行气动特性,采用相对运动方法模拟列车高速运行[22]. 建立计算域如图2所示,流体域长、宽、高分别为125.41H、71.64H和35.72H,磁浮列车头车鼻尖距离32.36H. 为保证尾流充分发展,尾车鼻尖距离计算域出口边界68.18H,计算域参数满足列车空气动力学计算标准TB/T 3503.4—2018.

    图  2  计算域示意
    Figure  2.  Computational domain

    为获得准确计算结果,需对计算域赋予正确的边界条件. 超导电动悬浮列车设计速度600 km/h,大于0.30 Ma,因此需要考虑空气压缩性. 入口边界选用压力远场边界,设置边界马赫数为0.49 Ma(对应600 km/h). 出口边界选用压力出口. 地面及U型轨道表面设置为滑移壁面,设置壁面滑移矢量与气体流动方向相同以避免轨道壁面效应影响. 车体表面为光滑壁面,计算域左、右侧及顶面设置为对称边界.

    网格离散方法中采用四面体网格,在头(尾)车、风挡及4个悬浮架附近设置加密区,头尾车加密区内最大体网格为100 mm,最大面网格为100 mm. 悬浮架最大面网格为50 mm,由于车体流线型特征较为复杂,曲率较大部分设置面网格尺寸为20 mm. 网格示意如图3.

    图  3  计算网格及加密区
    Figure  3.  Computational grid and area of encryption

    网格无关性验证共设置粗糙、中等及精细3类网格,具体网格尺寸如表1所示:为适应湍流模型,控制壁面第一层网格Y + 值在1附近,棱柱层统一采用0.02 mm,增长率设置为1.2,共30层. 采用无量纲压力系数Cp、升力系数CL及阻力系数CD评判磁浮列车网格无关性,如式(1)、(2)、(3)所示.

    表  1  网格参数
    Table  1.  Grid parameters
    网格质量 数量/
    (×107个)
    车体最大面网格/mm 加密区最大体网格/mm
    粗糙 4.74 100 200
    中等 6.60 100 100
    精细 7.92 80 100
    下载: 导出CSV 
    | 显示表格
    Cp=p0.5ρv2S, (1)
    CL=L0.5ρv2S, (2)
    CD=D0.5ρv2S, (3)

    式中:ρ为流体密度,25 ℃时空气取值为1.17 kg/m3p为测点压力;v为列车运行速度,取166.67 m/s;L为车体升力;D为车体阻力;S为车体截面积取9.19 m2

    提取车底中线压力系数如图4所示,其中横轴正方向为头车鼻尖. 车底压力系数在头车鼻尖处及一位悬浮架后端出现一定程度分离. 放大相关位置可以发现,尽管幅值有一定差异,但仍具有相同的规律特征.

    图  4  车底中线压力系数
    Figure  4.  Center line pressure coefficient beneath vehicle

    表2为3套网格无关性验证,其中下标h代表头车,t代表尾车,CL,h为头车升力系数,其余类似. 提取头(尾)车的升力及阻力进行计算,以精细网格计算结果为标准. 中等网格尾车升力最大相对误差达到5.85%,粗糙网格尾车阻力最大误差达到8.51%,3套网格误差均不超过10.00%,网格对计算结果影响较小. 为提高计算效率并保证准确性,计算采用中等网格策略.

    表  2  网格无关性验证
    Table  2.  Verification of grid independence
    网格 CL,h CL,t CD,h CD,t
    粗糙 0.0826 0.0665 0.0975 0.0828
    中等 0.0875 0.0755 0.0954 0.0861
    精细 0.0868 0.0713 0.0975 0.0905
    下载: 导出CSV 
    | 显示表格

    采用商业软件Fluent完成数值计算,使用定常方法计算磁浮列车稳定流动特征,采用基于剪切应力的SST k-ω湍流模型避免网格诱导分离. 考虑压缩性选用25 ℃时理想气体作为计算流体;采用基于高斯节点的方法进行梯度离散;动量方程、湍流动能方程及湍流耗散率方程采用二阶迎风离散格式.

    由于超导电动悬浮列车尚无实车试验数据,因此,采用TR08现场试验数据进行数值方法的验证[23]. 建立轨道及三编组磁浮列车气动外形模型,如图5所示. 边界条件除入口和滑移壁面速度更改为430 km/h,其余保持一致. 采用相同网格策略及边界层参数. 测点p1p2p3选择如图5所示,时均压力系数与试验相对误差如表3所示,最大误差为9.83%,满足工程误差需要,证明数值方法有效.

    图  5  TR08测点示意(单位:m)
    Figure  5.  Measuring points of TR08 (unit: m)
    表  3  时均压力系数计算方法验证
    Table  3.  Validation of time-averaged pressure coefficient calculation method
    测点 实测值 计算值 相对误差/%
    p1 0.240 0.21641 9.83
    p2 0.034 0.03294 3.13
    p3 0.016 0.01547 3.31
    下载: 导出CSV 
    | 显示表格

    以磁浮列车明线600 km/h运行为例,车体表面压力分布如图6所示. 由图6可知:头车鼻尖出现流动驻点使得鼻尖位置动压减少,静压增大,鼻尖附近出现较大的正压力;自头车鼻尖向后发展,车体上表面由正压逐渐回落至常压,并随着空气沿磁浮列车表面加速运动,车体附近流动逐步加速,车体流线型区域表面为低压;尾车流线型区域流动截面逐渐增大,引起膨胀效应产生负压,尾车鼻尖由于顶部及底部同时产生膨胀流动,正压分布区域相比头车鼻尖更大.

    图  6  磁浮列车表面压力分布
    Figure  6.  Pressure distribution in surface of maglev train

    磁浮列车车体底部主要为常压分布. 4个悬浮架底部则为不同程度负压,靠近悬浮架的车体底部呈现正压. 一位及四位悬浮架靠近磁浮列车鼻尖,流动截面变化速率更大,流动加剧,其底部负压明显大于二位及三位悬浮架.

    磁浮列车空间涡结构采用Q准则识别,如式(4).

    Q=12(Ω2FS2F) (4)

    式中:Ω为涡量张量,S为应变速率张量.

    提取车体附近Q=80 s−1等值面涡流,采用速度幅度染色车体及尾车附近涡流,如图7所示.

    图  7  磁浮列车空间涡流
    Figure  7.  Space vortex of maglev train

    车体附近主要有涡流A、B、C:其中涡流A出现在头车鼻尖附近轨道梁侧壁正上方,并沿车体向后上方延伸直到尾流区域散逸;涡流B为磁浮列车尾车上表面及四位悬浮架顶部膨胀流动产生的尾涡并向后延伸;涡流C为尾车底部产生的分离涡流,在较短距离内耗散.

    悬浮架附近主要有涡流D、E:悬浮架前端与轨道侧壁形成涡流D并向后扩散流动,逐步减慢,在风挡后方散逸消失;由于横向宽度变化,悬浮架与车体过渡段底部存在小型涡流E.

    纵截面向提取磁浮列车Y=0截面,如图8(a)所示;横截面提取任意相邻悬浮架中间截面、头尾车鼻尖前端及尾流充分发展位置,如图8(b)所示. 由图8可知:头车鼻尖前端挤压U型轨道内部气体,车体流线型表面向上导流加速,少部分气体流入车体底部;受轨道侧壁梁限制,头车压缩U型轨内部空气流动,在侧壁梁顶部逐步形成涡流A,沿流动方向,涡流A占主要并沿着车体表面逐步从侧面移动至车体顶面;尾车鼻尖产生涡流B,主要由于膨胀流动使得侧壁梁外侧空气流入后形成;涡流C为磁浮列车底部与轨道内车轮踏面形成的膨胀涡流;在−14.33H截面位置,涡流A及涡流B仍然存在,涡流C已基本耗散消失.

    图  8  不同截面涡量流线
    Figure  8.  Vortex streamlines of different sections

    通过上述对磁浮列车表面压力分布及整车流动分析发现,悬浮架外形对于车体附近流动存在一定影响. 为探究悬浮架附近详细气动特性,同时避免鼻尖对悬浮架流动的影响,以二位悬浮架为例分析悬浮架附近流动特征. 提取同一悬浮架表面压力云图及流动迹线,如图9所示.

    图  9  悬浮架附近压力及流动迹线
    Figure  9.  Pressure and flow trace near suspension frame

    图9可知:悬浮架前端将气流分为三部分,F部分气流导流至悬浮架上方,G部分气流被悬浮架外形压缩继续向后流动,此外还有H气体被导流压缩至悬浮架底部;悬浮架后端与车体过渡部分随着横向间隙增大,流速减小,出现扩张流动形成流动I.

    进一步分析悬浮架底部流动特征,以悬浮架几何中心为0点,分别提取图9中所示悬浮架前(后)端、车体过渡段中点及悬浮架三等分截面处压力及速度云图,如图10所示.

    图  10  悬浮架截面压力及速度云图
    Figure  10.  Pressure and velocity contours of suspension frame section

    图10可知:沿空气流动方向分析,车体两侧压力场及速度场分布左右对称;1.38H截面位置为悬浮架前端,此位置磁浮列车横向间隙较大,轨道内由于即将受到悬浮架外形影响产生一定的拥塞流动,压力呈现正压,轨道内流速明显降低;随着空气向后流动进入悬浮架与车体过渡位置1.02H,横向间隙开始减小,对横向间隙流动空间继续压缩,压力增大,速度减小,悬浮架底部流速逐渐回归到流场速度,压力也回归常压;0.47H、0及−0.47H截面为悬浮架导向间隙流动,由于横向间隙减小,流动速度有所增大,静压减小,底部均为不同程度负压,0位置对应风挡截面,车体上附近流速减小;随后再次进入悬浮架与车体过渡区域截面−1.02H,横向间隙开始增大,流速减小,轨道内负压逐步增大至常压;随后进入悬浮架后端−1.38H,横向间隙已增大至车体宽度,流速进一步减小,静压增大,横向间隙内主要为正压.

    有限元数值方法计算得到的磁浮列车气动荷载是由压力场的矢量和产生. 结果表明,磁浮列车车体与悬浮架附近压力分布差异较大,可以推断作用于磁浮列车车体与悬浮架的气动荷载同样有较大的差异. 同时车体结构也表明,磁浮列车车体与悬浮架之间通过二系悬挂连接,悬浮架与车体之间存在相对运动. 为更好地计算磁浮列车受到气动荷载作用下的响应,提出一种气动荷载提取方式(B),即将悬浮架外形受到的气动荷载区别于车体单独提取,并与常规仅区分车体的提取方式(A)获得的气动荷载进行对比. 将2种提取方式得到的气动荷载作用于车辆动力学模型中获得动力学响应.

    图11为不同提取方式下车体的气动荷载. 由图11可知,气动力均随速度增加而增大. 1) 整体提取方式中,头车与中间车升力差异较小,均为正升力;尾车则为负升力,气动升力幅值大小基本相同. 2) 分部件提取方式中,三车体均为升力,升力由大到小分别为头车、尾车以及中间车;分部件提取的头(尾)车车体气动升力约为整体提取方式的5倍. 3) 相比于整体提取方式,分部件提取方式时,中间车升力增大约1倍;不同提取方式下,头车及尾车阻力差异明显大于中间车;分部件提取得到的三车体阻力约为总阻力的75%,其中头车约占33%,尾车占26%,证明磁浮列车阻力主要由车体产生且头车阻力大于尾车阻力.

    图  11  不同提取方式下车体气动荷载
    Figure  11.  Aerodynamic load of vehicle body under different extraction methods

    分部件提取方式中,单独提取悬浮架所受气动荷载如图12所示. 一~四位悬浮架升力均为负值,为整车升力的−10.8倍,一位及四位悬浮架升力绝对值明显大于二位及三位悬浮架;一~四位悬浮架的阻力由大到小排列为四位、一位、二位和三位悬浮架,悬浮架外形产生的阻力均不超过4 kN.

    图  12  悬浮架气动荷载
    Figure  12.  Aerodynamic load of suspension frame

    三编组车辆动力学模型如图13所示,主要包括3个车体、4个悬浮架及超导磁体单元,其中,二位及三位悬浮架为铰接式悬浮架,同时承载2个车体. 超导磁体单元与悬浮架之间连接刚度相比二系悬挂刚度大一个数量级,因此,将超导磁体单元与构架间简化为刚性连接,并将两部件视为悬浮架整体. 以上构件均视为刚体,每个刚体考虑对应的横向、垂向、侧滚、点头及摇头5个方向自由度,共35个自由度. 由于本文磁浮列车明线运行,其横向力、侧滚力矩及摇头力矩均为0附近的极小值,因此,不再将此3个自由度荷载作为输入. 采用电磁力快速计算模型[24]对超导电动悬浮列车电磁力进行计算.

    图  13  磁浮列车动力学模型
    Figure  13.  Dynamics model of maglev train

    对应2种气动荷载提取方式(A、B),磁浮列车也提出2种气动荷载加载方式,分别为集中加载方式(a)和分部件加载方式(b),2种加载方式下磁浮列车受力分析如图14所示. 加载方式a下,整体提取的气动力Fair仅作用在车体,悬浮架则通过二系悬挂传递荷载;加载方式b中,车体及悬浮架分别受到车体气动力Fair_car及悬浮架气动力Fair_bogie作用. 2种加载方式下,磁浮列车气动荷载合力一致,但车体与悬浮架之间的受力有一定差异.

    图  14  不同加载方式示意
    Figure  14.  Different loading modes

    以悬浮架无气动荷载作用下稳定悬浮高度为基准,Y轴正向为相对于无气动荷载作用时上浮高度,一~四位悬浮架悬浮高度变化量如图15所示. 由于头车及中间车升力通过二系悬挂向下传递,带动一、二及三位悬浮架出现一定程度上浮,二位及三位悬浮架上升程度大于一位悬浮架,最大上升量达到5.50 mm,由于尾车负升力影响四位悬浮架下沉约6.97 mm.

    图  15  悬浮架悬浮高度变化量
    Figure  15.  Change in suspension height of suspension frame

    悬浮高度变化直接影响悬浮力,由于超导电动悬浮列车的被动悬浮稳定特征,2种加载方式下,悬浮架悬浮间隙变化量最大差异为1.05 mm. 四位悬浮架悬浮力变化量如图16所示.

    图  16  悬浮力变化量
    Figure  16.  Change in suspension force

    悬浮力大小与悬浮架位移呈反比例关系,四位悬浮架悬浮力增大,一、二及三位悬浮架则减小. 2种加载方式下,悬浮力最大变化了2.92 kN;一位及四位独立悬浮架悬浮力差异较小,二位及三位铰接式悬浮架差异较大;有、无气动荷载对比时悬浮力最大变化量达到15.98 kN.

    2种加载方式下,磁浮列车车体高度基本不变,悬浮架与车体之间的相对变化反映在二系悬挂空气弹簧力变化. 图17给出了2种气动荷载加载方式下磁浮列车二系空簧力变化量随速度变化曲线. 由图17可知:1) 以磁浮列车无气动荷载作用下空簧受力为基准,加载方式a作用下,尾车车体受到负升力,反映为下压,因此,四位悬浮架空簧力增大;其余位置受车体升力影响,空簧力均减小;加载方式a下,磁浮列车空簧力最大变化量为4.94 kN. 2) 加载方式b下,车体与悬浮架受力相反,空气弹簧呈减小拉伸趋势,所有空簧力均减小;一位及四位悬浮架为独立悬浮架,其前端及后端空簧变化量差异较小;二位及三位悬浮架为铰接式悬浮架,同时承载2个车体,空簧前后端存在一定差异,以二位悬浮架为例,前端头车升力较大,对应前端空簧力减少的也更多;加载方式b下,空簧力最大变化量达到18.7 kN,为加载方式a的2.86倍.

    图  17  空簧力变化量
    Figure  17.  Chang in air spring force

    基于某型磁浮列车几何模型,计算分析其明线气动特性及气动荷载,在此基础上探究2种加载方式下车辆响应,得出以下结论:

    1) U型轨包车结构直接影响磁浮列车车体附近流动特性. 头车鼻尖压缩轨道内空气并沿头车流线型表面向上导流至轨道外,在侧壁梁顶部形成涡流. 尾车后端膨胀分离流动则被一定程度限制在U型轨道内部.

    2) 超导电动悬浮列车车体和悬浮架宽度不一致,导致悬浮架附近流动规律出现明显变化. 悬浮架从前至后压力分布为正压—负压—正压,其中悬浮架底部及侧面主要以负压为主. 一位及四位悬浮架受头尾车流线型影响,底部压力远小于二位和三位底部压力.

    3) 基于气动特性差异,提出2种气动荷载提取方式. 以600 km/h速度运行为例,整体提取方式中,头车与中间车升力均为正并且相近,头车和中间车升力均大于尾车负升力;分部件提取方式中,车体受到升力均为正,头车和尾车大小相近,均大于中间车;四位悬浮架升力均为负,幅值大小依次为一位和四位升力相近,均大于二位及三位悬浮架升力. 分部件提取方式中,车体气动升力相比于整车提取方式增大约5倍,增加了气动荷载作用下的不确定性,悬浮架负升力幅值占车体升力的91.5%,悬浮架阻力约占整车的四分之一.

    4) 相比于整体气动荷载作用,分布荷载作用下磁浮列车响应主要反映为空气弹簧变化量. 悬浮架受气动荷载作用,高度变化量不超过7 mm,不同加载方式悬浮架高度变化基本一致. 分部件加载方式相比于整车加载方式,空簧力变化量增大1.86倍. 分部件加载方式可更加准确反应二悬挂空气弹簧受力特征,为车辆动力学提供更为准确的计算参考.

    致谢:

  • 图 1  “超前导洞 + 环向注浆 + 加长锚杆”支护示意

    Figure 1.  Support of “advance guide tunnel + circumferential grouting + lengthened anchor rod”

    图 2  围岩分区示意

    Figure 2.  Zoning scheme of surrounding rock

    图 3  岩石参数随塑性应变的劣化规律

    Figure 3.  Degradation of rock parameters with plastic strain

    图 4  围岩变形控制效果对比

    Figure 4.  Comparison of deformation control effects of surrounding rock

    图 5  无导洞无注浆时正洞围岩数值模拟计算结果

    Figure 5.  Numerical simulation results of main tunnel surrounding rock without guide tunnel and grouting

    图 6  超前导洞围岩数值模拟计算结果

    Figure 6.  Numerical simulation results of surrounding rock of advance guide tunnel

    图 7  正洞围岩(注浆后)数值模拟计算结果

    Figure 7.  Numerical simulation results of main tunnel surrounding rock (after grouting)

    图 8  软化模量McMφME的影响

    Figure 8.  Influence of softening modulus McMφ and ME

    图 9  注浆参数的影响

    Figure 9.  Influence of grouting parameters

    图 10  超前导洞半径r0的影响

    Figure 10.  Influence of radius r0 of advance guide tunnel

    图 11  掌子面间距的影响

    Figure 11.  Influence of distance between tunnel faces

    表  1  围岩计算参数

    Table  1.   Calculation parameters of surrounding rock

    材料 E/MPa c/MPa φ/(º) v 容重γ/kN·m−3 黏聚力软化模量Mc/MPa 内摩擦角软化模量Mφ/ MPa P0/MPa


    初始值 2 800 3.2 32
    0.2

    20.5

    400

    1 800

    12.2
    残余值 1 000 1.2 20
    注浆后 2 500 2.6 28
    下载: 导出CSV
  • [1] HE M C, SUI Q R, LI M N, et al. Compensation excavation method control for large deformation disaster of mountain soft rock tunnel[J]. International Journal of Mining Science and Technology, 2022, 32(5): 951-963. doi: 10.1016/j.ijmst.2022.08.004
    [2] 康永水,耿志,刘泉声,等. 我国软岩大变形灾害控制技术与方法研究进展[J]. 岩土力学,2022,43(8): 2035-2059.

    KANG Yongshui, GENG Zhi, LIU Quansheng, et al. Research progress on support technology and methods for soft rock with large deformation hazards in China[J]. Rock and Soil Mechanics, 2022, 43(8): 2035-2059.
    [3] 王志杰,周平,杨建民,等. 昔格达地层隧道围岩的失稳特征及变形控制工法[J]. 西南交通大学学报,2019,54(4): 757-768.

    WANG Zhijie, ZHOU Ping, YANG Jianmin, et al. Instability properties and deformation control methods of rocks surrounding xigeda strata[J]. Journal of Southwest Jiaotong University, 2019, 54(4): 757-768.
    [4] DENG H S, FU H L, SHI Y, et al. Countermeasures against large deformation of deep-buried soft rock tunnels in areas with high geostress: a case study[J]. Tunnelling and Underground Space Technology Incorporating Trenchless Technology Research, 2022, 119: 104238.1-104238.19.
    [5] SUN X M, ZHANG B, ZHANG Y, et al. Field tests and research of large deformation control measures for soft rock tunnel with high geostress[J]. IOP Conference Series: Earth and Environmental Science, 2020, 570(4): 042006.1-042006.9.
    [6] 李锦辉,张俊齐,魏强,等. 基于自感知纤维增强复合材料锚杆的隧道围岩松动圈识别[J]. 西南交通大学学报,2024,59(1): 11-19.

    LI Jinhui, ZHANG Junqi, WEI Qiang, et al. Loose Zone Identification for Surrounding Rock of Tunnels Using Self-Sensing Fiber Reinforced Plastic Anchors[J]. Journal of Southwest Jiaotong University, 2024, 59(1): 11-19.
    [7] 王明年,王志龙,桂登斌,等. 开挖方法影响下的深埋隧道形变压力计算方法[J]. 西南交通大学学报,2021,56(5): 1116-1124.

    WANG Mingnian, WANG Zhilong, GUI Dengbin, et al. Calculation Method of Deformation Load of Deep-Buried Tunnel under Influence of Excavation Method[J]. Journal of Southwest Jiaotong University, 2021, 56(5): 1116-1124.
    [8] 张梅,徐双永,张民庆,等. 高应力软岩隧道超前导洞法应力释放试验研究[J]. 现代隧道技术,2013,50(4): 68-75,83. doi: 10.3969/j.issn.1009-6582.2013.04.011

    ZHANG Mei, XU Shuangyong, ZHANG Qingmin, et al. Experimental study of stress release using a pilot heading for a highly stressed soft-rock tunnel[J]. Modern Tunnelling Technology, 2013, 50(4): 68-75,83. doi: 10.3969/j.issn.1009-6582.2013.04.011
    [9] 夏才初,金天垚,徐晨,等. 软岩隧道超前导洞应力释放力学机制及适用性[J]. 隧道建设(中英文),2020,40(增2): 1-9.

    XIA Caicu, JIN Tianqiu, Xuchen, et al. Mechanical mechanism and applicability of stress release for pilot heading in soft rock tunnel[J]. Tunnel Construction, 2020, 40(S2): 1-9.
    [10] 刘国庆. 木寨岭隧道软岩大变形段支护措施研究[J]. 现代隧道技术,2011,48(4): 135-141,147.

    LIU Guoqing. Study of support measures for large deformed section of muzhailing tunnel in soft rocks[J]. Modern Tunnelling Technology, 2011, 48(4): 135-141,147.
    [11] 周宗青,李术才,李利平,等. 导洞超前开挖施工性态空间效应分析[J]. 岩石力学与工程学报,2014,33(增1): 2611-2619.

    ZHOU Zongqing, LI Shucai, LI Liping, et al. Analysis of spatial effects of construction behavior on pilot tunnel excavation[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(S1): 2611-2619.
    [12] 马伟斌,张梅,赵有明,等. 高地应力软岩隧道超前洞室断面优化模拟研究[J]. 铁道工程学报,2011,28(9): 58-62.

    MA Weibin, ZHANG Mei, ZHAO Youming, et al. Simulation study on the section optimization of the super pilot tunnel in high geo-stress tunnel[J]. Journal of Railway Engineering Society, 2011, 28(9): 58-62.
    [13] 韩常领,夏才初,徐晨,等. 软岩隧道的超前导洞施工位置优化数值分析[J]. 地下空间与工程学报,2020,16(增1): 264-269,277.

    HAN Changling, XIA Caichu, XU Chen, et al. Numerical analysis on the optimization of the position of the pilot tunnel in soft rock[J]. Chinese Journal of Underground Space and Engineering, 2020, 16(S1): 264-269,277.
    [14] 崔光耀,魏杭杭,王明胜. 高地应力强风化炭质板岩隧道大变形控制现场试验研究[J]. 现代隧道技术,2022,59(3): 183-189,200.

    CUI Guangyao, WEI Hanghang, WANG Mingsheng. A field test study on the large deformation control for tunnels built in seriously-weathered carbonaceous slate strata with high geostress[J]. Modern Tunnelling Technology, 2022, 59(3): 183-189,200.
    [15] 王道远,刘佳,张逴,等. 高地应力深埋隧道断裂破碎带大变形控制方法现场试验研究[J]. 岩土工程学报,2020,42(4): 658-666.

    WANG Daoyuan, LIU Jia, ZHANG Chuo, et al. Field tests on large deformation control method for surrounding rock of deep tunnel in fault zone with high geostress[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(4): 658-666.
    [16] 尤显明,李沿宗. 极高地应力软岩隧道超前导洞应力释放及多层支护变形控制技术[J]. 隧道建设,2017,37(7): 832-837. doi: 10.3973/j.issn.1672-741X.2017.07.008

    YOU Xianming, LI Yanzong. Stress releasing and deformation control technology of multilayer support for pilot heading of a soft rock tunnel with extremely high ground stress[J]. Tunnel Construction, 2017, 37(7): 832-837. doi: 10.3973/j.issn.1672-741X.2017.07.008
    [17] WU Y S, TAN Z S, LI S M, Experimental study on the basic characteristics of tunnel in squeezing surrounding rock with large deformation[J]. hina Civil Engineering Journal, 2015, 48: 398-402.
    [18] 郑可跃,施成华,雷明锋,等. 考虑松动效应的高地应力构造破碎带隧道稳定性分析及大变形支护设计优化[J]. 岩石力学与工程学报,2021,40(8): 1603-1613.

    ZHENG Keyue, SHI Chenghua, LEI Mingfeng, et al. Stability analysis and support design optimization of large-deformation tunnels in structural fracture zones with high in situ stresses considering loose effect[J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(8): 1603-1613.
    [19] 李峰. 高地应力构造破碎带隧道大变形灾变机制及控制技术研究[J]. 铁道科学与工程学报,2021,18(5): 1222-1230.

    LI Feng. Study on the formation mechanism and large deformation control of tunnel in structural fracture zone with high in-situ stress[J]. Journal Of Railway Science And Engineering, 2021, 18(5): 1222-1230.
    [20] 张广泽,邓建辉,王栋,等. 隧道围岩构造软岩大变形发生机理及分级方法[J]. 工程科学与技术,2021,53(1): 1-12.

    ZHANG Guangze, DENG Jianhui, WANG Dong, et al. Mechanism and Classification of Tectonic-induced Large Deformation of Soft Rock Tunnels[J]. Advanced Engineering Sciences, 2021, 53(1): 1-12.
    [21] PAN Jiliang, GAO Zhaoning, REN Fenhua. Effect of strength criteria on surrounding rock of circular roadway considering strain softening and dilatancy[J]. Meitan Xuebao Journal of the China Coal Society, 2018, 43(12): 3293-3301.
    [22] CHEN L, MAO X B, LI M, et al. A new strain-softening constitutive model for circular opening considering plastic bearing behavior and its engineering application[J]. Mathematical Problems in Engineering, 2018, 2018: 4521492.1-4521492.13.
    [23] 陆银龙,王连国,杨峰,等. 软弱岩石峰后应变软化力学特性研究[J]. 岩石力学与工程学报,2010,29(3): 640-648.

    LU Yinlong, WANG Lianguo, YANG Feng, et al. Post-peak strain softening mechanical properties of weak rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(3): 640-648.
    [24] 郑颖人,孔亮. 岩土塑性力学[M]. 2版. 北京:中国建筑工业出版社,2019:73-79.
    [25] 徐芝纶. 弹性力学简明教程[M]. 北京:高等教育出版社,2013:70-72.
  • 期刊类型引用(3)

    1. 杨春雷,黄运华,丁军君. 曲线几何参数对不同类型货车转向架轮轨动力作用的影响分析. 振动工程学报. 2023(01): 159-169 . 百度学术
    2. 杨春雷,王开云,黄运华,丁军君. 重载货车通过曲线轨道时的悬挂特性分析. 铁道学报. 2022(06): 18-29 . 百度学术
    3. 杨春雷,黄运华,丁军君. 重载轨道曲线几何参数对轮轨耦合动力特性的影响. 交通运输工程学报. 2021(03): 215-227 . 百度学术

    其他类型引用(4)

  • 加载中
图(11) / 表(1)
计量
  • 文章访问数:  34
  • HTML全文浏览量:  21
  • PDF下载量:  5
  • 被引次数: 7
出版历程
  • 收稿日期:  2023-12-25
  • 修回日期:  2024-04-04
  • 网络出版日期:  2025-02-10

目录

/

返回文章
返回