CN114611293A - 一种结合传递系数法计算滑坡体中隧道结构荷载的方法 - Google Patents
一种结合传递系数法计算滑坡体中隧道结构荷载的方法 Download PDFInfo
- Publication number
- CN114611293A CN114611293A CN202210239814.5A CN202210239814A CN114611293A CN 114611293 A CN114611293 A CN 114611293A CN 202210239814 A CN202210239814 A CN 202210239814A CN 114611293 A CN114611293 A CN 114611293A
- Authority
- CN
- China
- Prior art keywords
- tunnel
- landslide
- arch
- sliding
- force
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012546 transfer Methods 0.000 title claims abstract description 13
- 239000011435 rock Substances 0.000 claims abstract description 21
- 230000005540 biological transmission Effects 0.000 claims abstract description 13
- 239000002689 soil Substances 0.000 claims description 9
- 239000003673 groundwater Substances 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 20
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000013461 design Methods 0.000 abstract description 3
- 230000002787 reinforcement Effects 0.000 abstract description 2
- 230000006378 damage Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 1
- 239000010438 granite Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
本发明公开了一种结合传递系数法计算滑坡体中隧道结构荷载的方法,包括:(1)首先计算常规条件下,隧道结构所受的围岩压力值;(2)通过传递系数法确定滑坡剩余下滑力与剩余抗滑力;根据传递系数法的假设,滑坡分块间作用力总是以平行于滑动面的方向作用于隧道结构。本发明很大程度上改进了位于滑坡体内的隧道缺乏荷载理论解析的问题,针对隧道轴向与滑坡滑动方向正交、且隧道位于滑动面上方的几何关系,利用传递系数法计算滑坡推力,并引入在滑坡治理中提出的剩余抗滑力概念计算隧道远山侧荷载。本发明可以近似计算出滑坡体中的隧道承受滑坡推力荷载,为隧道结构的设计、加固提供了新的计算方法。
Description
技术领域
本发明涉及到隧道-滑坡相互作用研究领域,具体涉及到对位于滑坡体中的隧道结构荷载的计算方法。
背景技术
当前越来越多的隧道建立在滑坡或潜在滑坡体当中,针对隧道轴向与滑坡滑动方向呈正交几何位置关系的情况,《公路隧道设计规范》只对常规荷载做了主要计算方法,并未针对滑坡推力这一特殊荷载进行计算方法的考虑。此外,当前对于隧道滑坡体系的研究,缺少对隧道结构受力的理论计算方法,而多采用数值模拟的方法进行隧道结构受力分析。
发明内容
对于这种情况,本发明提出一种针对隧道滑坡正交情况下的隧道支护结构受力的理论计算方法,将位于滑坡体内部的隧道结构按照传递系数法计算的滑坡推力结合普氏平衡拱理论来计算隧道结构承受的荷载。
常规的滑坡推力分布形式通常有三角形、矩形和梯形三种,本发明采用符合实际情况的三角形分布,因此将作用在隧道结构上的滑坡推力以三角形荷载考虑。
(1)首先计算常规条件下,隧道结构所受的围岩压力值。假设隧道结构形成图1所示的压力平衡拱。假设取一点Q(x,y),由于平衡压力拱上点受力平衡,根据普氏理论,得到其抛物线方程为:
式中:q—拱轴线上部岩体自重产生的均布荷载(kN/m2);T—平衡压力拱拱顶截面的水平推力(kN);x—平衡拱抛物线任意一点Q的x坐标值;y—平衡拱抛物线任意一点Q的y坐标值。
当隧道发生变形后,平衡拱侧壁稳定时,平衡拱轴线的最大高度为
b1—平衡拱轴线的最大高度(m);a1—自然平衡拱轴线的最大跨度(m);f—普氏经验系数,与岩体的完整性和地下水有关;
根据朗肯土压力理论,破裂面角度与竖直夹角为:
其中:
因此,平衡拱侧壁不稳定时,拱顶继续坍塌,此时自然平衡拱轴线的最大跨度:
由式(1)—式(5)得:
Q点至Q点下方隧道距离为:
此时,可以得出Q点下方隧道支护结构处所受竖向荷载为:
式中:γ—Q点距离下方隧道之间松散岩体的加权重度平均值(kN/m3)。
为了求得隧道拱顶所受整体荷载,对其进行积分,可得:
隧道边墙承受的侧向围岩压力为通过朗肯土压力计算得出:
式中:e1—隧道侧墙所受梯形荷载上方最小值;e2—隧道侧墙所受梯形荷载下方最大值;Ph—隧道边墙所受全部侧向围岩压力。
(2)通过传递系数法确定滑坡剩余下滑力Pi与剩余抗滑力Fi:
传递系数法属于刚体极限平衡分析法,通常分为强度储备法和超载法两种,属于刚体极限平衡分析法。在使用传递系数法时,考虑隧道-滑坡体系的最不利截面,并简化为平面应变问题。传递系数法的假设条件如下:
(1)将滑坡稳定性问题视为平面应变问题;
(2)滑动力以平行于滑动面的剪应力和垂直于滑动面的正应力集中作用于滑动面上;
(3)视滑坡体为理想刚塑性材料,认为整个加载过程中,滑坡体不会发生任何变形,一旦沿滑动面剪应力达到其剪切强度,则滑坡体开始沿滑动面产生剪切变形;
(4)滑动面的破坏服从莫尔-库仑破坏准则;
(5)剩余下滑力方向与滑动面倾角一致,剩余下滑力为负值时则传递的剩余下滑力为0;
(6)沿整个滑动面满足静力的平衡条件,但不满足力矩平衡条件。
根据传递系数法的假设,滑坡分块间作用力总是以平行于滑动面的方向作用于隧道结构。
通过常规传递系数法计算可确定隧道近山侧所受滑坡推力,根据传递系数法采用逆向求解,在隧道滑坡体系中隧道结构荷载计算时,引入与剩余下滑力相对应的剩余抗滑力概念,以计算远山侧隧道结构受力。在隧道下侧分块采用逆向计算剩余抗滑力,上侧近山侧采用剩余下滑力计算剩余下滑力,即可同时确定隧道近山侧结构所受滑坡推力及远山侧所受岩土抗力。若远山侧分块剩余抗滑力大于0,则证明隧道远山侧拱墙除了常规围岩压力外,还要受到下侧分块的岩土抗力作用;若远山侧分块的剩余抗滑力小于0,则证明下侧分块跟随滑坡体进行滑动,可以认为其不对隧道边墙产生其他力的作用。
Pi-1-第i-1块滑体下滑力(kN);
Wi-1-第i-1块滑体的自重(kN/m3);
αi-1-第i-1块滑体滑面倾角(°);
ψi-2-第i-2块滑体下滑力传递系数;
ci-1-第i-1块滑体粘聚力(kN);
li-1-第i-1块滑体滑面长度(m);
下方分块的剩余抗滑力:
特别地,当Fi+2≤0时,令Fi+2=0,即:
作用在隧道结构范围内的滑坡推力为:
将作用隧道结构右侧的滑坡推力按水平和竖直方向分解:
将有可能作用在隧道结构左侧的滑坡剩余抗滑力(若Fi-1>0)按水平和竖直方向分解:
因此,在普氏拱理论计算得到的常规围岩压力基础上考虑传递系数法计算的剩余下滑推力后,隧道结构右侧水平合力为:
隧道结构左侧水平合力为:
当Fi+2>0时:qh左=Ph+qh2 式(25)
当Fi+2≤0时:qh左=Ph 式(26)
拱部竖直合力为:
与现有技术相比,本发明具有那些技术效果,请结合本发明的技术优势进行下描述。
本发明很大程度上改进了位于滑坡体内的隧道缺乏荷载理论解析的问题,针对隧道轴向与滑坡滑动方向正交、且隧道位于滑动面上方的几何关系,利用传递系数法计算滑坡推力,并引入在滑坡治理中提出的剩余抗滑力概念计算隧道远山侧荷载。本发明最终提出了统一计算理论,可以近似计算出滑坡体中的隧道承受滑坡推力荷载,为隧道结构的设计、加固提供了新的计算方法。
附图说明
图1为本发明提出的滑坡体中的隧道衬砌受力示意图。
图2为本发明提出的滑坡体中隧道衬砌受力分解示意图。
图3为围岩普氏压力拱示意图。
图4为传递系数法对各条块下滑力计算简图。
图5为本发明中各几何参数的说明图。
具体实施方式
以下结合附图和实施例对本发明进行详细说明。
本发明以1993年6月的广西资源县城南端大塘湾发生的一处滑坡为计算案例。经勘查查明,滑坡滑体面积1700m2,体积约为6800m3。滑体主要由花岗岩全风化层、砾岩全风化层、残坡积层粗砂、砾砂、碎石土构成。为计算的带滑体中隧道衬砌结构的承担的荷载,需要将其所受围岩压力和滑坡推力分别计算并叠加得到等效荷载,如图1、图2所示。
根据本发明的计算步骤,第1步,首先根据普氏拱理论,对一定埋深的隧道承受的围岩压力值进行计算,参数说明见图3。根据该工程的现场勘查资料,其相应的参数如下:
表1计算参数一览表
Table1 Table of calculation parameters
假设隧道结构位于第4计算条块之中,将以上参数代入式(9)、(12),计算得到隧道所承受的竖向荷载Pv和水平荷载ph分别为4 130.4KN和248.5KN。
第2步,对滑坡体进行分条,应用传递系数法式(13)、(15)计算得到各计算条块的剩余下滑力和剩余抗滑力,见图4。由于隧道结构位于第4计算条块之中,来自上方第3条块的剩余下滑力即p3为473.327KN,下方第5条块的剩余抗滑力为43.090KN。各条块的剩余下滑力和剩余抗滑力计算结果见表2。
表2滑坡推力计算结果
Table2 Calculation results
第3步,根据图5计算简图所示,将表2计算结果代入式(19)~(22),得出隧道结构受上侧第3条块剩余下滑力作用的等效竖直荷载qv1和等效水平荷载qh1为944 970KN和56365KN,隧道结构受下侧第5条块剩余抗滑力作用的等效竖直荷载qv2和等效水平荷载qh2为121 800KN和408 160KN。
第4步,将计算结果代入式(23)、(24)、(26),滑坡推力产生的荷载与围岩压力叠加后得到隧道结构左侧承担水平荷载qh左和右侧承担水平荷载为qh右为56 613KN和40 840KN。隧道结构拱部承担竖直荷载qv计算结果为827 300KN。
Claims (1)
1.一种结合传递系数法计算滑坡体中隧道结构荷载的方法,将位于滑坡体内部的隧道结构按照传递系数法计算的滑坡推力结合普氏平衡拱理论来计算隧道结构承受的荷载;采用三角形分布,将作用在隧道结构上的滑坡推力以三角形荷载考虑;其特征在于:该方法的实施过程如下:
(1)首先计算隧道结构所受的围岩压力值;取一点Q(x,y),由于平衡压力拱上点受力平衡,根据普氏理论,得到其抛物线方程为:
式中:q—拱轴线上部岩体自重产生的均布荷载;T—平衡压力拱拱顶截面的水平推力;x—平衡拱抛物线任意一点Q的x坐标值;y—平衡拱抛物线任意一点Q的y坐标值;
当隧道发生变形后,平衡拱侧壁稳定时,平衡拱轴线的最大高度为
b1—平衡拱轴线的最大高度;a1—自然平衡拱轴线的最大跨度;f—普氏经验系数,与岩体的完整性和地下水有关;
根据朗肯土压力理论,破裂面角度与竖直夹角为:
其中:
平衡拱侧壁不稳定时,拱顶继续坍塌,自然平衡拱轴线的最大跨度:
由式(1)—式(5)得:
Q点至Q点下方隧道距离为:
此时,得出Q点下方隧道支护结构处所受竖向荷载为:
式中:γ—Q点距离下方隧道之间松散岩体的加权重度平均值;
为了求得隧道拱顶所受整体荷载,进行积分得:
隧道边墙承受的侧向围岩压力为通过朗肯土压力计算得出:
式中:e1—隧道侧墙所受梯形荷载上方最小值;e2—隧道侧墙所受梯形荷载下方最大值;Ph—隧道边墙所受全部侧向围岩压力;
(2)通过传递系数法确定滑坡剩余下滑力Pi与剩余抗滑力Fi:
在使用传递系数法时,考虑隧道-滑坡体系的最不利截面,并简化为平面应变问题;根据传递系数法的假设,滑坡分块间作用力总是以平行于滑动面的方向作用于隧道结构;
在隧道下侧分块采用逆向计算剩余抗滑力,上侧近山侧采用剩余下滑力计算剩余下滑力,即可同时确定隧道近山侧结构所受滑坡推力及远山侧所受岩土抗力;若远山侧分块剩余抗滑力大于0,则证明隧道远山侧拱墙除了常规围岩压力外,还要受到下侧分块的岩土抗力作用;若远山侧分块的剩余抗滑力小于0,则证明下侧分块跟随滑坡体进行滑动,认为不对隧道边墙产生其他力的作用;
Pi-1—第i-1块滑体下滑力;
Wi-1—第i-1块滑体的自重;
αi-1—第i-1块滑体滑面倾角;
ψi-2—第i-2块滑体下滑力传递系数;
ci-1—第i-1块滑体粘聚力;
li-1—第i-1块滑体滑面长度;
下方分块的剩余抗滑力:
特别地,当Fi+2≤0时,令Fi+2=0,即:
作用在隧道结构范围内的滑坡推力为:
将作用隧道结构右侧的滑坡推力按水平和竖直方向分解:
将有可能作用在隧道结构左侧的滑坡剩余抗滑力按水平和竖直方向分解:
在普氏拱理论计算得到的常规围岩压力基础上考虑传递系数法计算的剩余下滑推力后,隧道结构右侧水平合力为:
隧道结构左侧水平合力为:
当Fi+2>0时:qh左=Ph+qh2 式(25)
当Fi+2≤0时:qh左=Ph 式(26)
拱部竖直合力为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210239814.5A CN114611293B (zh) | 2022-03-12 | 2022-03-12 | 一种结合传递系数法计算滑坡体中隧道结构荷载的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210239814.5A CN114611293B (zh) | 2022-03-12 | 2022-03-12 | 一种结合传递系数法计算滑坡体中隧道结构荷载的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114611293A true CN114611293A (zh) | 2022-06-10 |
CN114611293B CN114611293B (zh) | 2024-03-29 |
Family
ID=81862341
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210239814.5A Active CN114611293B (zh) | 2022-03-12 | 2022-03-12 | 一种结合传递系数法计算滑坡体中隧道结构荷载的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114611293B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116720348A (zh) * | 2023-06-02 | 2023-09-08 | 中铁十二局集团第七工程有限公司 | 基于关键点位移监测计算支架现浇转体桥不平衡力矩的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101476466A (zh) * | 2008-12-24 | 2009-07-08 | 中国科学院水利部成都山地灾害与环境研究所 | 巨厚层滑坡预应力锚索抗滑隧道施工方法及防治新结构 |
CN105136370A (zh) * | 2015-07-13 | 2015-12-09 | 北京工业大学 | 一种深埋非对称连拱隧道土压力荷载确定方法 |
CN107609319A (zh) * | 2017-10-30 | 2018-01-19 | 中国电建集团成都勘测设计研究院有限公司 | 根据隧道衬砌变形确定滑坡推力的方法 |
CN112883459A (zh) * | 2021-01-17 | 2021-06-01 | 北京市政路桥股份有限公司 | 一种隧道-滑坡体-抗滑桩体系稳定性耦合分析方法 |
WO2021189612A1 (zh) * | 2020-03-23 | 2021-09-30 | 青岛理工大学 | 一种基坑开挖边坡安全系数的计算方法 |
WO2021213280A1 (zh) * | 2020-04-21 | 2021-10-28 | 长安大学 | 一种单洞四车道公路隧道断面参数的设计方法 |
-
2022
- 2022-03-12 CN CN202210239814.5A patent/CN114611293B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101476466A (zh) * | 2008-12-24 | 2009-07-08 | 中国科学院水利部成都山地灾害与环境研究所 | 巨厚层滑坡预应力锚索抗滑隧道施工方法及防治新结构 |
CN105136370A (zh) * | 2015-07-13 | 2015-12-09 | 北京工业大学 | 一种深埋非对称连拱隧道土压力荷载确定方法 |
CN107609319A (zh) * | 2017-10-30 | 2018-01-19 | 中国电建集团成都勘测设计研究院有限公司 | 根据隧道衬砌变形确定滑坡推力的方法 |
WO2021189612A1 (zh) * | 2020-03-23 | 2021-09-30 | 青岛理工大学 | 一种基坑开挖边坡安全系数的计算方法 |
WO2021213280A1 (zh) * | 2020-04-21 | 2021-10-28 | 长安大学 | 一种单洞四车道公路隧道断面参数的设计方法 |
CN112883459A (zh) * | 2021-01-17 | 2021-06-01 | 北京市政路桥股份有限公司 | 一种隧道-滑坡体-抗滑桩体系稳定性耦合分析方法 |
Non-Patent Citations (1)
Title |
---|
夏艳华;白世伟;: "传递系数法在滑坡治理削坡方案设计中的应用", 岩石力学与工程学报, no. 1, 15 June 2008 (2008-06-15) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116720348A (zh) * | 2023-06-02 | 2023-09-08 | 中铁十二局集团第七工程有限公司 | 基于关键点位移监测计算支架现浇转体桥不平衡力矩的方法 |
CN116720348B (zh) * | 2023-06-02 | 2023-11-28 | 中铁十二局集团第七工程有限公司 | 基于关键点位移监测计算支架现浇转体桥不平衡力矩的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114611293B (zh) | 2024-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Three-dimensional failure mechanism of a rectangular cavity in a Hoek–Brown rock medium | |
Lin et al. | Three-dimensional stress-transfer mechanism and soil arching evolution induced by shield tunneling in sandy ground | |
Peila | A theoretical study of reinforcement influence on the stability of a tunnel face | |
CN104480962B (zh) | 一种有限填土挡墙的土压力分布计算方法 | |
Wang et al. | Roof collapse mechanisms for a shallow tunnel in two-layer rock strata incorporating the influence of groundwater | |
Lee et al. | Undrained stability of surface strip footings above voids | |
CN107590357B (zh) | 一种不同施工阶段隧道稳定性判别的方法 | |
Shin et al. | Interaction between tunnel supports and ground convergence—Consideration of seepage forces | |
Xiao et al. | Undrained stability of strip footing above voids in two-layered clays by finite element limit analysis | |
CN108536942B (zh) | 一种软岩隧道开挖面变形计算方法 | |
Please et al. | Fracturing of an Euler–Bernoulli beam in coal mine pillar extraction | |
CN108830014B (zh) | 一种含砂卵石地层盾构隧道开挖面最小支护力计算方法 | |
CN114611293A (zh) | 一种结合传递系数法计算滑坡体中隧道结构荷载的方法 | |
Ghanbari et al. | Improvement parameters in dynamic compaction adjacent to the slopes | |
CN113204828A (zh) | 一种考虑滑坡坡底建筑物稳定性的整体加固设计方法 | |
CN111460678A (zh) | 一种水泥胶结材料的充填速率优化方法 | |
Du et al. | Numerical analysis of the effects of rock bolts on stress redistribution around a roadway | |
CN110671146B (zh) | 一种胶结充填体宽度确定方法 | |
Baoxu et al. | Analytical solution for stress distribution in cementitious backfills considering stope inclinations | |
Xia et al. | Damage evolution mechanism of extraction structure during mining gently dipped orebody by block caving method | |
Zhan et al. | Stability analysis and reinforcement of a high-steep rock slope with faults: numerical analysis and field monitoring | |
Sun et al. | Undrained stability analysis of dual unlined horseshoe-shaped tunnels in non-homogeneous clays using lower bound limit analysis method | |
Zhang et al. | Experimental study on the stability of surrounding soft rocks of gob-side entry retaining in fully mechanized caving | |
CN113239444A (zh) | 一种考虑滑坡坡底建筑物稳定性的反压加固设计方法 | |
CN111985021A (zh) | 一种盾构开挖面的遍布节理流固耦合的安全度分析方法 |
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 |