CN109387385B - 一种三维tbm护盾外表面围岩荷载识别方法 - Google Patents

一种三维tbm护盾外表面围岩荷载识别方法 Download PDF

Info

Publication number
CN109387385B
CN109387385B CN201811243507.4A CN201811243507A CN109387385B CN 109387385 B CN109387385 B CN 109387385B CN 201811243507 A CN201811243507 A CN 201811243507A CN 109387385 B CN109387385 B CN 109387385B
Authority
CN
China
Prior art keywords
shield
tbm
load
shield outer
country rock
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811243507.4A
Other languages
English (en)
Other versions
CN109387385A (zh
Inventor
刘鹤
刘泉声
唐旭海
刘滨
黄兴
罗慈友
桑昊旻
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan Institute of Rock and Soil Mechanics of CAS
Original Assignee
Wuhan Institute of Rock and Soil Mechanics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan Institute of Rock and Soil Mechanics of CAS filed Critical Wuhan Institute of Rock and Soil Mechanics of CAS
Priority to CN201811243507.4A priority Critical patent/CN109387385B/zh
Publication of CN109387385A publication Critical patent/CN109387385A/zh
Application granted granted Critical
Publication of CN109387385B publication Critical patent/CN109387385B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M99/00Subject matter not provided for in other groups of this subclass
    • G01M99/005Testing of complete machines, e.g. washing-machines or mobile phones
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21DSHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
    • E21D9/00Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
    • E21D9/003Arrangement of measuring or indicating devices for use during driving of tunnels, e.g. for guiding machines
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21DSHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
    • E21D9/00Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
    • E21D9/06Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21FSAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00Methods or devices for use in mines or tunnels, not covered elsewhere
    • E21F17/18Special adaptations of signalling or alarm devices

Landscapes

  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)
  • Excavating Of Shafts Or Tunnels (AREA)

Abstract

本发明公开了一种三维TBM护盾外表面围岩荷载识别方法,涉及深埋隧道护盾式TBM高效掘进和安全控制的技术领域。本方法是:①在TBM护盾内表面预设测点安装应变监测传感器;②获得测点应变实测数据εm;③建立TBM护盾有限元计算模型;④通过护盾外表面分区得到荷载识别参数;⑤通过牛顿迭代和Moore‑Penrose广义逆求解得到护盾外表面围岩荷载;⑥对荷载积分得到护盾外表面摩阻力。本发明的计算精度高,可视化效果好,适合于现场TBM护盾外表面荷载的识别,能够为开挖隧道安全施工、卡机预警和事故处理提供指导。

Description

一种三维TBM护盾外表面围岩荷载识别方法
技术领域
本发明涉及深埋隧道护盾式TBM高效掘进和安全控制的技术领域,尤其涉及一种三维TBM护盾外表面围岩荷载识别方法。
背景技术
TBM(Tunnel boring machine)工法相对传统钻爆法具有优质、高效、经济、安全和环保等优势,在深埋长大隧道建设过程中成为优先的考虑对象。但是,TBM在穿越深埋高应力软弱破碎地层时,在开挖扰动作用下围岩容易发生挤压大变形。一旦围岩变形量超过TBM护盾与围岩之间的预留间隙,围岩便对护盾造成挤压,当机器推挤系统的额定推力不足以克服护盾表面的摩阻力,即发生卡机事故。由于护盾的存在限制了对围岩状态的观察,一般只有当机器完全被卡住才能察觉卡机事故的发生。卡机事故的处理不仅延误工期,而且还会造成严重的经济损失。围岩的变形通常具有显著的时效性,护盾表面围岩荷载大小的演变也是一个逐渐发展的过程。因此,如果在掘进过程中能够有效地识别护盾表面的围岩荷载,便能够在围岩荷载发展前期及时地采取措施避免卡机,这对隧道施工安全管理、进度控制和成本控制等都有着重大意义。但是,围岩与护盾之间的间隙狭小,且存在强烈的相互作用,在护盾外表面安装传感器直接监测围岩荷载目前还存在困难。考虑到获取护盾内表面应变相对简单,可以通过数值反演来确定护盾外表面的荷载分布,进而确定护盾表面摩阻力的大小。
发明内容
本发明的目的是提供一种三维TBM护盾表面围岩荷载识别方法,为护盾式TBM开挖隧道安全施工、卡机预警和事故处理提供指导。该方法识别精度高,实现简单,适合于三维TBM护盾外表面围岩荷载识别问题,同时也能够扩展到其它结构的荷载识别问题。
本发明采用下述技术方案:
具体地说,本方法包括以下步骤:
①在TBM护盾内表面预设测点安装应变监测传感器
所述的预设测点数目应尽量多且均匀分布,避开护盾内部移动机械构件作用区域;
所述的应变监测传感器采用表面应变计或者光纤应变传感器,传感器在每个测点分别沿护盾内表面轴向和环向布置,也可以沿单个方向布置,同时需要对传感器施加保护措施;
②获得测点应变实测数据εm
所述的实测应变数据利用现有的应变数据采集仪获取,以便实现长期工作,现场需将采集仪进行防水密封处理,并配备UPS不间断电源;
③建立TBM护盾有限元计算模型
所述的有限元计算模型通过有限元建模程序或者软件生成,应考虑护盾实际结构及其内部与盾体相连的主要部件,如焊接在护盾上的圈梁、支撑等,模型采用实体单元网格划分,模型尺寸和材料属性与实际一致。
④通过护盾外表面分区得到荷载识别参数
所述的护盾外表面分区采用多边形进行分区,各多边形区域的节点值作为荷载识别参数PI,记录荷载识别参数个数d,任意区域内部的荷载通过该区域节点值线性插值或者非线性插值得到,具体插值形式视分区方式而定;
⑤通过牛顿迭代和Moore-Penrose广义逆求解得到护盾外表面围岩荷载
所述的牛顿迭代和Moore-Penrose广义逆求解流程如下:
A、给定荷载识别参数PI初始值PI 0,收敛精度δ。采用有限元计算软件进行正分析,得到测点应变的计算值εc
B、得到系统目标方程:
r(PI)=εmc (1)
C、根据牛顿迭代公式,对于第k次迭代:
式中,J为雅可比矩阵,即系统目标方程的一阶偏导数矩阵,具有如下形式:
式中,w为测点的个数,ej(j=1,2…d)表示第j个元素为1的单位向量, 为摄动量;
D、引入Moore-Penrose广义逆对方程(2)进行求解,采用下式求解:
当w<d时:
当w≥d时:
E、检查测点应变计算值是否满足收敛要求||rk||/||εm||<δ。如果满足则停止计算,以作为反演结果,通过插值得到护盾外表面围岩荷载,否则进入下一轮迭代。
⑥对荷载积分得到护盾外表面摩阻力
所述的摩阻力Ff的计算方法如下:
本发明具有下列优点和积极效果:
①能够有效获取TBM掘进过程中护盾外表面围岩荷载的三维分布及摩阻力的大小,并可以嵌入到现场操控平台实现可视化,为隧道安全施工提供保障;
②通过进一步设定一个卡机预警值Fw:Fw=S*Ft,其中,S为TBM掘进安全系数,Ft为TBM推进系统的额定推力,当计算得到的摩阻力Ff大于Fw时,便进行卡机预警,及时提醒技术人员采取措施避免发生卡机,这对于深埋隧道护盾式TBM高效掘进和安全控制具有重要意义。
总之本发明的计算精度高,可视化效果好,适合于现场TBM护盾外表面荷载的识别,能够为开挖隧道安全施工、卡机预警和事故处理提供指导。
附图说明
图1为本方法的工作流程图;
图2-1为本发明实施例中TBM护盾外表面实际荷载分布云图;
图2-2为本发明实施例中TBM护盾外表面荷载识别结果云图。
英译汉
1、TBM:Tunnel Boring Machine,隧道掘进机。是利用回转刀具开挖,同时破碎洞内围岩及掘进,形成整个隧道断面的一种新型、先进的隧道施工机械;相对于常用的方法,TBM集钻、掘进、支护于一体,使用电子、信息、遥测、遥控等高新技术对全部作业进行制导和监控,使掘进过程始终处于最佳状态。在国际上,已广泛应用于水利水电、矿山开采、交通、市政、国防等工程中。
具体实施方式
下面结合实施例和附图详细说明。
如图1,本方法具体实施步骤如下:
①在TBM护盾内表面预设测点安装应变监测传感器-101
所述的预设测点数目应尽量多且均匀分布,避开护盾内部移动机械构件作用区域;
所述的应变监测传感器采用表面应变计或者光纤应变传感器,传感器在每个测点分别沿护盾内表面轴向和环向布置,也可以沿单个方向布置,同时需要对传感器施加保护措施;
本实施例采用表面应变计,在护盾内表面均匀布置36个测点,每个测点应变计沿护盾内表面环向布置。
②获得测点应变实测数据εm-102
所述的实测应变数据利用现有的应变数据采集仪获取,以便实现长期工作,现场需将采集仪进行防水密封处理,并配备UPS不间断电源。
③建立TBM护盾有限元计算模型-103
所述的有限元计算模型通过有限元建模程序或者软件生成(本实施例中采用Abaqus),应考虑护盾实际结构及其内部与盾体相连的主要部件,如焊接在护盾上的圈梁、支撑等;模型采用实体单元网格划分,模型尺寸和材料属性与实际一致;本实施例中护盾长12m,护盾外径2.9m,护盾厚度0.05m,护盾弹性模量200GPa,密度7500kg/m3,泊松比0.2。护盾的厚度方向划分3层实体单元,同时也需要对模型施加合适的边界条件。
④通过护盾外表面分区得到荷载识别参数-104
所述的护盾外表面分区采用多边形进行分区,各多边形区域的节点值作为荷载识别参数PI,记录荷载识别参数个数d,任意区域内部的荷载通过该区域节点值线性插值或者非线性插值得到。
本实施例中采用矩形对护盾外表面进行分区,沿护盾轴向和环向分别划分m和n个区域(本实施例中m=n=9),那么荷载识别参数的个数d=90,任意区域Si,j(i=1,2…9,j=1,2…9)内部荷载通过节点值Pi,j,Pi,j+1,Pi+1,j+1和Pi+1,j线性插值的计算方法为:
其中,P(z,θ)为区域Si,j内部任意一点的荷载值,(z,θ)为该点在柱坐标系下的轴向和环向坐标。
⑤通过牛顿迭代和Moore-Penrose广义逆求解得到护盾外表面围岩荷载;
所述的牛顿迭代和Moore-Penrose广义逆求解具体计算流程如下:
A、给定荷载识别参数PI初始值PI 0,收敛精度δ。本实施例中PI 0=0,δ=0.01。由式(1)可得到物护盾外表面所有区域的分布荷载,采用有限元计算软件(本实施例中采用Abaqus有限元软件)进行求解,得到36测点环向应变的计算值εc
B、得到系统目标方程:
r(PI)=εmc (1)
C、根据牛顿迭代公式,对于第k次迭代:
式中,J为雅可比矩阵,即系统目标方程的一阶偏导数矩阵,具有如下形式:
式中,w=36,ej(j=1,2…90)表示第j个元素为1的单位向量, 为摄动量;
D、引入Moore-Penrose广义逆对方程(2)进行求解,由于w<d,采用下式求解:
当时:
E、检查测点应变计算值是否满足收敛要求||rk||/||εm||<δ。如果满足则停止计算,以作为反演结果,代入式(1)得到护盾外表面围岩荷载,否则进入下一轮迭代。
护盾外表面荷载识别结果如图2-2,与实际荷载云图(图2-1)对比可以看出,一种三维TBM护盾外表面围岩荷载识别方法在该具体实施例中能够准确地识别护盾表面的荷载。
⑥对荷载积分得到护盾外表面摩阻力。
摩阻力Ff的计算方法如下:
本实施例中护盾外表面摩擦系数μ=0.3,通过上式计算得到该实施例中护盾表面的摩阻力大小为188MPa。该实施例中TBM推进系统的额定推力为Fw=150MPa,那么可以判断机器处于卡机状态,从识别的荷载分布云图可以判断围岩与护盾的挤压部位位于护盾后方顶部,对现场事故处理具有很好的指导作用。同时也表明,在识别荷载分布和摩阻力大小的基础上,该护盾荷载识别方法可以为现场隧道安全掘进和卡机预警提供指导。

Claims (1)

1.一种三维TBM护盾表面围岩荷载识别方法,其特征在于包括以下步骤:
①在TBM护盾内表面预设测点安装应变监测传感器(101)
所述的预设测点数目应尽量多且均匀分布,避开护盾内部移动机械构件作用区域;
所述的应变监测传感器采用表面应变计或者光纤应变传感器,传感器在每个测点分别沿护盾内表面轴向和环向布置,也可以沿单个方向布置,同时需要对传感器施加保护措施;
②获得测点应变实测数据εm(102)
所述的实测应变数据利用现有的应变数据采集仪获取,以便实现长期工作,现场需将采集仪进行防水密封处理,并配备UPS不间断电源;
③建立TBM护盾有限元计算模型(103)
所述的有限元计算模型通过有限元建模程序或者软件生成,应考虑护盾实际结构及其内部与盾体相连的部件,包括焊接在护盾上的圈梁和支撑,模型采用实体单元网格划分,模型尺寸和材料属性与实际一致;
④通过护盾外表面分区得到荷载识别参数(104)
所述的护盾外表面分区采用多边形进行分区,各多边形区域的节点值作为荷载识别参数PI,记录荷载识别参数个数d,任意区域内部的荷载通过该区域节点值线性插值或者非线性插值得到,具体插值形式视分区方式而定;
⑤通过牛顿迭代和Moore-Penrose广义逆求解得到护盾外表面围岩荷载(105)
所述的牛顿迭代和Moore-Penrose广义逆求解流程如下:
A、给定荷载识别参数PI初始值PI 0,收敛精度δ,采用有限元计算软件进行正分析,得到测点应变的计算值εc
B、得到系统目标方程:
r(PI)=εmc (1);
C、根据牛顿迭代公式,对于第k次迭代:
-(Jk(PI k))·ΔPI k=rk(PI k) (2);
式中,J为雅可比矩阵,即系统目标方程的一阶偏导数矩阵,具有如下形式:
式中,w为测点的个数,ej(j=1,2…d)表示第j个元素为1的单位向量, 为摄动量;
D、引入Moore-Penrose广义逆对方程(2)进行求解,采用下式求解:
当w<d时:
ΔPI k=-(Jk(PI k))T·(Jk(PI k)·(Jk(PI k))T)-1·rk(PI k) (4);
当w≥d时:
ΔPI k=-((Jk(PI k))T·Jk(PI k))-1·Jk(PI k)·rk(PI k) (5);
E、检查测点应变计算值是否满足收敛要求||rk||/||εm||<δ,如果满足则停止计算,以PI k+1=PI k+ΔPI k作为反演结果,通过插值得到护盾外表面围岩荷载,否则进入下一轮迭代;
⑥对荷载积分得到护盾外表面摩阻力(106)
所述的摩阻力Ff的计算方法如下:
公式中:
n:沿TBM护盾外表面环向分区个数;
m:沿TBM护盾外表面轴向分区个数;
μ:TBM护盾外表面的摩擦系数;
P(z,θ):TBM护盾外表面任意一点的荷载值;
z:柱坐标系下TBM护盾外表面任意一点的轴向坐标;
θ:柱坐标系下TBM护盾外表面任意一点的环向坐标。
CN201811243507.4A 2018-10-24 2018-10-24 一种三维tbm护盾外表面围岩荷载识别方法 Active CN109387385B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811243507.4A CN109387385B (zh) 2018-10-24 2018-10-24 一种三维tbm护盾外表面围岩荷载识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811243507.4A CN109387385B (zh) 2018-10-24 2018-10-24 一种三维tbm护盾外表面围岩荷载识别方法

Publications (2)

Publication Number Publication Date
CN109387385A CN109387385A (zh) 2019-02-26
CN109387385B true CN109387385B (zh) 2019-09-03

Family

ID=65426935

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811243507.4A Active CN109387385B (zh) 2018-10-24 2018-10-24 一种三维tbm护盾外表面围岩荷载识别方法

Country Status (1)

Country Link
CN (1) CN109387385B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110006348B (zh) * 2019-03-19 2021-05-11 天津大学前沿技术研究院有限公司 一种预警式地铁盾构隧道管片接缝防水性能智能监测方法
CN110824552B (zh) * 2019-11-14 2021-10-29 山东大学 一种搭载于盾构机的不良地质探测系统及方法
CN111832821A (zh) * 2020-07-09 2020-10-27 山东大学 一种tbm卡机风险预测方法及系统
CN112177621B (zh) * 2020-09-11 2021-11-23 山东大学 Tbm掘进过程护盾围岩垮塌灾害实时预警系统及方法
CN112362520B (zh) * 2020-10-30 2022-01-04 武汉大学 一种有限元-离散元耦合数值模拟程序(fdem)输入参数快速标定方法
CN113008157A (zh) * 2021-04-26 2021-06-22 黄河勘测规划设计研究院有限公司 隧道掘进机护盾内表面变形监测方法
CN114756939B (zh) * 2022-04-21 2023-03-07 中铁二院工程集团有限责任公司 一种敞开式tbm开挖隧道的围岩荷载计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000079675A1 (de) * 1999-06-16 2000-12-28 Fronius Schweissmaschinen Produktion Gmbh & Co. Kg Verfahren zum erkennen einer an einen wechselrichter bedarfsweise anschaltbaren last und entsprechendes lasterkennungssystem
CN105952465A (zh) * 2016-05-26 2016-09-21 中国科学院武汉岩土力学研究所 一种围岩与全断面隧道掘进机护盾相互作用过程监测方法
CN107506557A (zh) * 2017-09-01 2017-12-22 宁波市交通规划设计研究院有限公司 一种隧道软弱围岩力学参数的非线性反演分析方法
CN107607082A (zh) * 2017-10-24 2018-01-19 成都理工大学 Tbm施工围岩变形监测系统
CN207879337U (zh) * 2018-02-24 2018-09-18 黄河勘测规划设计有限公司 具有围岩收敛变形监测系统的护盾式tbm

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000079675A1 (de) * 1999-06-16 2000-12-28 Fronius Schweissmaschinen Produktion Gmbh & Co. Kg Verfahren zum erkennen einer an einen wechselrichter bedarfsweise anschaltbaren last und entsprechendes lasterkennungssystem
CN105952465A (zh) * 2016-05-26 2016-09-21 中国科学院武汉岩土力学研究所 一种围岩与全断面隧道掘进机护盾相互作用过程监测方法
CN107506557A (zh) * 2017-09-01 2017-12-22 宁波市交通规划设计研究院有限公司 一种隧道软弱围岩力学参数的非线性反演分析方法
CN107607082A (zh) * 2017-10-24 2018-01-19 成都理工大学 Tbm施工围岩变形监测系统
CN207879337U (zh) * 2018-02-24 2018-09-18 黄河勘测规划设计有限公司 具有围岩收敛变形监测系统的护盾式tbm

Also Published As

Publication number Publication date
CN109387385A (zh) 2019-02-26

Similar Documents

Publication Publication Date Title
CN109387385B (zh) 一种三维tbm护盾外表面围岩荷载识别方法
Zhu et al. Challenges and development prospects of ultra-long and ultra-deep mountain tunnels
Liu et al. In situ observation of spalling process of intact rock mass at large cavern excavation
Zhu et al. Large-scale geomechanical model testing of an underground cavern group in a true three-dimensional (3-D) stress state
Zhang et al. Geo-mechanical model test on the collaborative bearing effect of rock-support system for deep tunnel in complicated rock strata
Jongpradist et al. Development of tunneling influence zones for adjacent pile foundations by numerical analyses
CN103410516B (zh) 小直径隧洞全断面开挖围岩变形预警工法
CN106595918B (zh) 一种盾构隧道管片外土压力的长期监测装置及方法
CN109357939B (zh) 盾构空推现场试验方法
CN205593565U (zh) 隧道工程全寿命监测装置
Zhao et al. Effects of water inrush from tunnel excavation face on the deformation and mechanical performance of shield tunnel segment joints
CN114322819A (zh) 光纤光栅传感器、深部巷道围岩的应变监测方法及系统
CN110907079A (zh) 一种采动应力的动态监测系统及方法
Lorenzo Longitudinal beam response of concrete segmental linings simultaneously backfilled with bicomponent grouts
Wagner et al. Direct and distributed strain measurements inside a shotcrete lining: concept and realisation
Duddeck Application of numerical analyses for tunnelling
Djizanne et al. Some aspects of the hydro-mechanical behaviour of Callovo-Oxfordian (COx) claystone around a gallery parallel to the principal horizontal minor stress
Gong et al. Tunnel boring machines (TBMs) in difficult grounds
Hesser et al. Measurements of acoustic emission and deformation in a repository of nuclear waste in salt rock
Lunardi et al. Tunnel monitoring system—a contribution for the preparation of guidelines
Zhou et al. Physical model tests of the surrounding rock deformation and fracture mechanism of mixed-face ground under TBM tunneling
Zhou et al. Dynamic deformation and fracture characteristics of a deep roadway surrounding rock based on the machine vision monitoring method
CN113030400A (zh) 一种富水断层隧道掌子面突水演化模拟试验装置及方法
Le The effect of forepole reinforcement on tunnelling-induced movements in clay
Makarov et al. Zonal failure structure near the deep openings

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant