CN113935155A - 一种核心土通长黏结型对拉锚杆轴力分布的计算方法 - Google Patents
一种核心土通长黏结型对拉锚杆轴力分布的计算方法 Download PDFInfo
- Publication number
- CN113935155A CN113935155A CN202111130166.1A CN202111130166A CN113935155A CN 113935155 A CN113935155 A CN 113935155A CN 202111130166 A CN202111130166 A CN 202111130166A CN 113935155 A CN113935155 A CN 113935155A
- Authority
- CN
- China
- Prior art keywords
- anchor rod
- axial force
- length
- counter
- pulling
- 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.)
- Pending
Links
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
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Piles And Underground Anchors (AREA)
Abstract
本发明提供一种核心土通长黏结型对拉锚杆轴力分布的计算方法,该方法包括:对锚杆受力微单元进行受力分析,得出相应的锚固体控制微分方程及轴力解析式;根据受力微单元的锚固体控制微分方程及轴力解析式,按受力微单元分段处变形连续且剪应力与剪位移是线弹性关系,建立通长黏结型对拉锚杆力学模型;根据力学模型,代入地层及锚杆参数,计算通长黏结式对拉锚杆轴力。本发明提出了核心土通长黏结型对拉锚杆力学模型,符合对拉锚杆实际工作状态,并能考虑双侧壁法侧导洞开挖顺序的影响,用于预测某一位置处对拉锚杆的轴力大小,为对拉锚杆的设计与施工提供参考,可以进一步得出通长黏结型对拉锚杆的轴力沿杆体全长的分布曲线,仅输入简单的地勘参数,就能够得到轴力分布曲线。
Description
技术领域
本发明涉及岩土工程测算技术领域,具体涉及一种核心土通长黏结型对拉锚杆轴力分布的计算方法。
背景技术
双侧壁法是一种应用于软弱岩层大跨度隧道的开挖工法,是指将大跨隧道划分为三个小跨断面,依次开挖左右两侧导坑,最后开挖预留核心土;在具体的施工中,左右侧导洞开挖将对核心土产生剧烈扰动,使得核心土处于单向或双向的不利受力状态;若核心土发生大范围的损伤破坏,将急剧增加大跨隧道的失稳风险,因此,双侧壁法的核心就是保证预留核心土的稳定牢固。为保证双侧壁法大跨隧道的安全修建,通常采用通长黏结型对拉锚杆(简称对拉锚杆)对核心土进行有效的支护。
通长黏结型对拉锚杆直接贯穿核心土,两端施加预应力,其布设范围一般是两侧导洞的内侧拱脚至内侧拱腰,布设间隔为0.5m-1.5m,锚杆长度为所处位置核心土宽度,按照围岩条件的不同,施加的预紧力范围为50-150kN。
但目前关于通长黏结型对拉锚杆的力学特性,特别是其沿杆体纵长轴力的分布规律尚不明确,需要将通长黏结型对拉锚杆与无黏结型锚杆、传统注浆锚杆进行区分开来,并考虑侧导洞开挖顺序,设计一种通长黏结型对拉锚杆轴力的分析方法。
发明内容
本发明的目的提供一种核心土通长黏结型对拉锚杆轴力分布的计算方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明采用如下技术方案:
一种核心土通长黏结型对拉锚杆轴力分布的计算方法,包括以下步骤:
S1、对锚杆受力微单元进行受力分析,得出相应的锚固体控制微分方程及轴力解析式;
S2、根据受力微单元的锚固体控制微分方程及轴力解析式,按受力微单元分段处变形连续且剪应力与剪位移是线弹性关系,建立通长黏结型对拉锚杆力学模型;
S3、根据力学模型,代入地层及锚杆参数,计算通长黏结式对拉锚杆轴力。
进一步地,对锚杆受力微单元进行受力分析,得出相应的锚固体控制微分方程及轴力解析式,包括:
S11、将对拉锚杆与浆体剥离,并分为混凝土段和围岩段;
S12、单独对混凝土段进行受力分析,得出对应的锚固体控制微分方程及轴力解析式:
根据对拉锚杆的实际受力处于平衡状态,可得:
其中:ds为通长黏结型对拉锚杆的直径,dr为取出微单元体的长度,σs(r)为对拉锚杆的轴应力,τ(r)为锚杆或浆体界面的剪应力;
其中:Ea是通长黏结型对拉锚杆的弹性模量,Se为锚杆的形变量;
假定锚杆或浆体界面剪应力与剪位移的关系处于线弹性状态,得:
其中:Ss为锚杆或浆体界面剪切位移;
假定锚杆与浆体的变形是协调的,得:Ss=Se (6)
联立式(4)-(6),得通长黏结型对拉锚杆锚固体系的微分控制方程:
其中:Gc为混凝土段的剪切模量;
对于通长黏结型对拉锚杆,混凝土段边界方程为:
其中:P为通长黏结型对拉锚杆两端的预紧力,A为杆体横断面的面积,且A=πds2/4;
带入混凝土段长度r=a,则混凝土段边界方程为:
其中:Pa为混凝土段和围岩段分段处的轴力;
将式(8)带入式(7),得:
其中:待定系数
带入混凝土段长度r=a,对于混凝土段,在混凝土段和围岩段连接处锚杆或浆体界面剪切位移表示为:
S13、单独对围岩段进行受力分析,得出对应的锚固体控制微分方程及轴力解析式:
根据步骤S12,得围岩段的控制微分方程:
其中:Gs为围岩段的剪切模量;
对于长度为l的锚杆,围岩段的边界方程为:
将式(14)带入式(13),得:
带入混凝土段长度r=a,对于围岩段,在混凝土段和围岩段连接处锚杆或浆体界面剪切位移表示为:
将混凝土段和围岩段合成整体进行受力分析,按混凝土段与围岩段变形连续以得出分段处轴力值:
进一步地,根据受力微单元的锚固体控制微分方程及轴力解析式,按受力微单元分段处变形连续且剪应力与剪位移是线弹性关系,建立通长黏结型对拉锚杆力学模型,包括:
S21、将式(18)带入步骤S12和S13中,得出通长黏结型对拉锚杆的轴应力解析式:
S22、将锚杆轴应力乘以杆体截面积,得到通长黏结型对拉锚杆轴力的力学模型解析式:Pr=Aσs (20)。
进一步地,基于通长黏结型对拉锚杆轴力的力学模型,结合对拉锚杆长度和对应的锚杆轴力可建立相应的轴力分布曲线。
由以上技术方案可知,本发明通过提出了核心土通长黏结型对拉锚杆力学模型,符合对拉锚杆实际工作状态,用于预测某一位置处对拉锚杆的轴力大小,为对拉锚杆的设计与施工提供参考,并能考虑双侧壁法侧导洞开挖顺序的影响,可以进一步得出通长黏结型对拉锚杆的轴力沿杆体全长的分布曲线,仅输入简单的地勘参数,就能够得到轴力分布曲线。
附图说明
图1为双侧壁法核心土通长黏结型对拉锚杆的示意图
图2为本发明的步骤流程示意图;
图3为通长黏结型对拉锚杆力学模型示意图;
图4为通长黏结型对拉锚杆现场实测布置图;
图5为某山岭隧道根据本发明计算方法建立相应的轴力分布曲线。
图中:1、左导洞;2、右导洞;3、预留核心土;4、对拉锚杆;5、锚垫板;6、初支混凝土;A为混凝土段;B为围岩段。
具体实施方式
下面结合附图对本发明的一种优选实施方式做详细的说明。
如图1所示的双侧壁法核心土通长黏结型对拉锚杆的示意图,大跨隧道采用双侧壁法进行开挖,左导洞1先行开挖,右导洞2后行开挖,预留核心土3最后开挖;通长黏结型对拉锚杆4用来支护核心土5,对拉锚杆贯穿核心土体,包括杆体和设置在左右导洞内侧的锚垫板,在两端施加有预应力P;对拉锚杆杆体材料为HRB335钢筋,其屈服延展率可以达到18%;先行导洞开挖后,即施做初支混凝土6,使得通长黏结型对拉锚杆存在两种锚固层介质,可以分为混凝土段和围岩段。
如图2所示的核心土通长黏结型对拉锚杆轴力分布的计算方法,包括以下步骤:
S1、对锚杆受力微单元进行受力分析,得出相应的锚固体控制微分方程及轴力解析式;
具体的,如图3所示,先将对拉锚杆与浆体剥离,并分为混凝土段和围岩段,且做出以下合理假定:
(1)一般认为锚杆或浆体是最危险的界面,因此本发明仅仅研究此界面的力学响应;
(2)另外根据既有实测,施做对拉锚杆,并全长注浆后,核心土受力状态得到显著改善,内部相对变形小于8mm。因此,本发明假定锚杆及浆体和岩土体的相对变形较小,且认为锚杆与浆体的变形是相互协调的;
(3)工程实践表明,对拉锚杆由于钢筋杆体延展性好,材料强度很高,通常不会进入塑性屈服阶段,因此本发明认为其处于弹性工作状态,由此得出的轴力解析公式较为简洁,方便应用与推广;
(4)采用双侧壁法施工的大跨隧道,左、右两侧导洞一般分前后依次开挖,一般滞后距离超过两倍洞径,本发明假定左导洞先行。通长黏结型对拉锚杆在后行右导洞开挖后立即施做,此时先行左侧导洞初支已经封闭成环,强度刚度较大;而后行右导洞尚处于毛洞状态,或初支施作不久,支护结构还远未到达设计强度值,将其看成普通围岩。
(5)依据上述假定及力学性能的优劣,将通长黏结型对拉锚杆的锚固介质区分为混凝土段、围岩段两段,其中混凝土段段为靠近先行左导洞较强的混凝土锚固介质,围岩段段为靠近后行右导洞较弱的围岩锚固介质,对拉锚杆长度为l,混凝土段段长度为混凝土段。混凝土段、围岩段两段可以分别称作混凝土段和围岩段,两者主要区别是两段锚杆/浆体界面能够提供的摩阻力不一样,进而造成剪切模量不同,混凝土段段剪切模量为Gc,围岩段段剪切模量为Gs。剪切模量与锚杆/浆体界面的摩阻性能与法向应力密切相关,一般通过现场拉拔测试来确定取值大小。
然后先单独对混凝土段进行受力分析,得出对应的锚固体控制微分方程及轴力解析式:
根据对拉锚杆的实际受力处于平衡状态,可得:
其中:ds为通长黏结型对拉锚杆的直径,dr为取出微单元体的长度,σs(r)为对拉锚杆的轴应力,τ(r)为锚杆或浆体界面的剪应力;
其中:Ea是通长黏结型对拉锚杆的弹性模量,Se为锚杆的形变量;
假定锚杆或浆体界面剪应力与剪位移的关系处于线弹性状态,得:
其中:Ss为锚杆或浆体界面剪切位移;
假定锚杆与浆体的变形是协调的,得:Ss=Se (6)
联立式(4)-(6),得通长黏结型对拉锚杆锚固体系的微分控制方程:
其中:Gc为混凝土段的剪切模量;
对于通长黏结型对拉锚杆,混凝土段边界方程为:
其中:P为通长黏结型对拉锚杆两端的预紧力,A为杆体横断面的面积,且A=πds2/4;
带入混凝土段长度r=a,则混凝土段边界方程为:
其中:Pa为混凝土段和围岩段分段处的轴力;
将式(8)带入式(7),得:
其中:待定系数
带入混凝土段长度r=a,对于混凝土段,在混凝土段和围岩段连接处锚杆或浆体界面剪切位移表示为:
再单独对围岩段进行受力分析,得出对应的锚固体控制微分方程及轴力解析式:
根据上述步骤,得围岩段的控制微分方程:
其中:Gs为围岩段的剪切模量;
对于长度为l的锚杆,围岩段的边界方程为:
将式(14)带入式(13),得:
带入混凝土段长度r=a,对于围岩段,在混凝土段和围岩段连接处锚杆或浆体界面剪切位移表示为:
将混凝土段和围岩段合成整体进行受力分析,按混凝土段与围岩段变形连续以得出分段处轴力值:
S2、根据受力微单元的锚固体控制微分方程及轴力解析式,按受力微单元分段处变形连续且剪应力与剪位移是线弹性关系,建立通长黏结型对拉锚杆力学模型;
具体的,将式(18)带入步骤S1中,得出通长黏结型对拉锚杆的轴应力解析式:
再将锚杆轴应力乘以杆体截面积,得到通长黏结型对拉锚杆轴力的力学模型解析式:Pr=Aσs (20)。
S3、根据力学模型,代入地层及锚杆参数,计算通长黏结式对拉锚杆轴力。
在具体的测算中,基于通长黏结型对拉锚杆轴力的力学模型,结合对拉锚杆长度和对应的锚杆轴力还可以建立相应的轴力分布曲线。
如某山岭隧道,采用的是通长黏结型对拉锚杆,其力学参数如表1所示;
表1对拉锚杆设计参数
采用本发明计算方法,将地层及对拉锚杆参数代入,能够解出通长黏结型对拉锚杆轴力的解析解。
如图4所示,采用锚杆测力计进行现场实测,在对拉锚杆上布置A1、A2、A3、A4、A5共5个钢筋测力计;如图5所示,将通长黏结型对拉锚杆轴力计算值和现场实测值进行对比分析,由本发明理论曲线与现场实测散点可知,本发明计算方法得出通长黏结型对拉锚杆的轴力分布与实测结果基本一致,显示出可靠的预测性能,能够为设计施工提供指导。
以上,本发明建立了双侧壁法核心土通长注浆对拉锚杆的力学模型,区别于黏结型或者部分黏结型锚杆,认为锚杆通过与浆体的摩阻作用来起到锚固效果,具体可以通过界面剪切模量来反映,力学模型简洁明了,便于工程应用。
同时,本发明还提出了通长黏结型对拉锚杆的双向受拉受力模式,区别于传统注浆锚杆,符合对拉锚杆实际工作状态,假定剪应力与剪位移是线弹性关系,计算公式较为简单,更容易被设计和施工人员理解和接受。
且本发明提出了通长黏结型对拉锚杆轴力的计算方法,考虑了因大跨隧道双侧壁法左右导洞施工顺序不同,而造成的对拉锚杆锚固介质的不同,该方法计算结果与实际情况较为一致,提高通长黏结型对拉锚杆设计的科学性与经济性。
以上所述实施方式仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明的权利要求书确定的保护范围内。
Claims (4)
1.一种核心土通长黏结型对拉锚杆轴力分布的计算方法,其特征在于,包括以下步骤:
S1、对锚杆受力微单元进行受力分析,得出相应的锚固体控制微分方程及轴力解析式;
S2、根据受力微单元的锚固体控制微分方程及轴力解析式,按受力微单元分段处变形连续且剪应力与剪位移是线弹性关系,建立通长黏结型对拉锚杆力学模型;
S3、根据力学模型,代入地层及锚杆参数,计算通长黏结式对拉锚杆轴力。
2.根据权利要求1所述的核心土通长黏结型对拉锚杆轴力分布的计算方法,其特征在于,对锚杆受力微单元进行受力分析,得出相应的锚固体控制微分方程及轴力解析式,包括:
S11、将对拉锚杆与浆体剥离,并分为混凝土段和围岩段;
S12、单独对混凝土段进行受力分析,得出对应的锚固体控制微分方程及轴力解析式:
根据对拉锚杆的实际受力处于平衡状态,可得:
其中:ds为通长黏结型对拉锚杆的直径,dr为取出微单元体的长度,σs(r)为对拉锚杆的轴应力,τ(r)为锚杆或浆体界面的剪应力;
其中:Ea是通长黏结型对拉锚杆的弹性模量,Se为锚杆的形变量;
假定锚杆或浆体界面剪应力与剪位移的关系处于线弹性状态,得:
其中:Ss为锚杆或浆体界面剪切位移;
假定锚杆与浆体的变形是协调的,得:Ss=Se (6)
联立式(4)-(6),得通长黏结型对拉锚杆锚固体系的微分控制方程:
其中:Gc为混凝土段的剪切模量;
对于通长黏结型对拉锚杆,混凝土段边界方程为:
其中:P为通长黏结型对拉锚杆两端的预紧力,A为杆体横断面的面积,且A=πds2/4;
带入混凝土段长度r=a,则混凝土段边界方程为:
其中:Pa为混凝土段和围岩段分段处的轴力;
将式(8)带入式(7),得:
带入混凝土段长度r=a,对于混凝土段,在混凝土段和围岩段连接处锚杆或浆体界面剪切位移表示为:
S13、单独对围岩段进行受力分析,得出对应的锚固体控制微分方程及轴力解析式:
根据步骤S12,得围岩段的控制微分方程:
其中:Gs为围岩段的剪切模量;
对于长度为l的锚杆,围岩段的边界方程为:
将式(14)带入式(13),得:
带入混凝土段长度r=a,对于围岩段,在混凝土段和围岩段连接处锚杆或浆体界面剪切位移表示为:
将混凝土段和围岩段合成整体进行受力分析,按混凝土段与围岩段变形连续以得出分段处轴力值:
4.根据权利要求1所述的核心土通长黏结型对拉锚杆轴力分布的计算方法,其特征在于,基于通长黏结型对拉锚杆轴力的力学模型,结合对拉锚杆长度和对应的锚杆轴力可建立相应的轴力分布曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111130166.1A CN113935155A (zh) | 2021-09-26 | 2021-09-26 | 一种核心土通长黏结型对拉锚杆轴力分布的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111130166.1A CN113935155A (zh) | 2021-09-26 | 2021-09-26 | 一种核心土通长黏结型对拉锚杆轴力分布的计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113935155A true CN113935155A (zh) | 2022-01-14 |
Family
ID=79276956
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111130166.1A Pending CN113935155A (zh) | 2021-09-26 | 2021-09-26 | 一种核心土通长黏结型对拉锚杆轴力分布的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113935155A (zh) |
-
2021
- 2021-09-26 CN CN202111130166.1A patent/CN113935155A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pan et al. | Safety factor assessment of a tunnel face reinforced by horizontal dowels | |
Zhao et al. | Mechanism analysis and control technology of surrounding rock failure in deep soft rock roadway | |
Zhang et al. | Stress evolution in roadway rock bolts during mining in a fully mechanized longwall face, and an evaluation of rock bolt support design | |
CN105138767B (zh) | 一种深埋非对称小间距隧道围岩压力计算方法 | |
Dias | Convergence-confinement approach for designing tunnel face reinforcement by horizontal bolting | |
Li et al. | Stability of plane strain tunnel headings in soils with tensile strength cut-off | |
Karakus et al. | Back analysis for tunnelling induced ground movements and stress redistribution | |
Li et al. | Analytical model of shear behaviour of a fully grouted cable bolt subjected to shearing | |
Song et al. | Beam–spring structural analysis for the design of a tunnel pre-reinforcement support system | |
Liu et al. | Internal force calculation and supporting parameters sensitivity analysis of side piles in the subway station excavated by Pile-Beam-Arch method | |
Shou et al. | A study of jacking force for a curved pipejacking | |
Chen et al. | Effects of rock pillar width on the excavation behavior of parallel tunnels | |
CN102129499B (zh) | 隧道掘进机刀盘驱动载荷计算方法 | |
CN106050240A (zh) | 山区软弱围岩小净距隧道中夹岩稳定性控制方法 | |
Kanik et al. | Importance of numerical analyses for determining support systems in tunneling: A comparative study from the trabzon-gumushane tunnel, Turkey | |
Xia et al. | Numerical-modeling-based assessment of the impact of two-end-type cable support on failure characteristics of yield pillars | |
Panet et al. | Convergence-confinement method for tunnel design | |
Cui et al. | Prevention of the ground subsidence by using the foot reinforcement side pile during the shallow overburden tunnel excavation in unconsolidated ground | |
Epel et al. | The influence of face and shield annulus pressure on tunnel liner load development | |
Kitchah et al. | Finite difference analysis of an advance core pre-reinforcement system for Toulon's south tube | |
CN113935155A (zh) | 一种核心土通长黏结型对拉锚杆轴力分布的计算方法 | |
Kong et al. | Analysis on the breaking law of soft and thick roof of fully mechanized top-coal caving face | |
Luo et al. | Mechanical properties of rock bolt and analysis for the full-process of sliding failure based on rock mass absolute displacement | |
CN109681251B (zh) | 全断面开挖软土隧洞掌子面玻璃纤维锚杆设计方法 | |
Wang et al. | Analytical approach and field monitoring for mechanical behaviors of pipe roof reinforcement |
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 |