深埋隧道围岩形变压力计算方法研究

王明年; 王志龙; 张霄; 赵思光; 刘大刚; 童建军

doi:10.11779/CJGE202001009

深埋隧道围岩形变压力计算方法研究 English Version

王明年^{1, 2,},
王志龙^{1, 2},
张霄^{1, 2},
赵思光^{1, 2},
刘大刚^{1, 2},
童建军^{1, 3, ,}

1.
西南交通大学隧道工程教育部重点实验室,四川　成都610031
2.
西南交通大学土木工程学院,四川　成都610031
3.
西南交通大学峨眉校区土木工程系,四川　峨眉山614202

基金项目:

国家自然科学基金项目 51878567

国家自然科学基金项目 51878568

国家自然科学基金项目 51578458

中国铁路总公司科技研究开发计划课题 2017G007-F

中国铁路总公司科技研究开发计划课题 2017G007-H

详细信息

作者简介:
王明年（1965— ）,男,教授,博士生导师,主要从事地下工程方面的教学与科研工作。E-mail：wangmingnian@swjtu.edu.cn

通讯作者:
童建军, E-mail：jjtong@163.com

中图分类号: TU456;U25
计量
- 文章访问数: 567
- HTML全文浏览量: 71
- PDF下载量: 320
出版历程
- 收稿日期: 2019-01-06
- 网络出版日期: 2022-12-07
- 刊出日期: 2019-12-31

Method for calculating deformation pressure of surrounding rock of deep-buried tunnels

WANG Ming-nian^{1, 2,},
WANG Zhi-long^{1, 2},
ZHANG Xiao^{1, 2},
ZHAO Si-guang^{1, 2},
LIU Da-gang^{1, 2},
TONG Jian-jun^{1, 3, ,}

1.
Key Laboratory of Transportation Tunnel Engineering, Ministry of Education, Southwest Jiaotoug University, Chengdu 610031, China
2.
School of Civil Engineering, Southwest Jiaotoug University, Chengdu 610031, China
3.
Department of Civil Engineering, Emei Branch, Southwest Jiaotong University, E’meishan 614202, China

摘要

摘要: 通过大量的深埋隧道围岩形变压力现场实测数据,采用数理统计方法,分析了围岩级别、隧道跨度、隧道埋深对围岩形变压力的影响及其变化规律。通过多元非线性回归分析方法,推导建立了竖向形变压力计算公式。通过分析可知,形变压力约为《铁路隧道设计规范》计算荷载的84%～92%,且围岩级别越高两者比值越小。形变压力侧压力系数近似服从正态分布,与《铁路隧道设计规范》中侧压力系数相比,各围岩条件下侧压力系数均有一定的增大。研究成果对我国隧道及地下工程支护结构设计具有一定的参考指导作用。
- 深埋隧道 /
- 现场实测 /
- 形变压力 /
- 隧道埋深 /
- 隧道跨度 /
- 围岩级别
Abstract: The pressure of surrounding rock is an important research content of tunnel engineering, and it includes deformation pressure. The study on the deformation pressure is of great significance to ensure the safety and economic construction of tunnels and underground engineering. The method of mathematical statistics is used to analyze the influence degrees and variation laws of grade of surrounding rock, tunnel span and tunnel depth on the deformation pressure through a large number of filed measured data. The formula for calculating the vertical deformation pressure is established by using the multivariate nonlinear regression analysis method. The deformation pressure is about 84%~92% times the value by the "Code for design of railway Tunnels" according to the analysis. As the grade of the surrounding rock increases, the ratio decreases continuously. The lateral pressure coefficient of deformation loads obeys the normal distribution. The lateral pressure coefficient of each surrounding rock condition has a certain increase compared with that in the "Code for design of railway tunnels". The research achievements have certain reference guidance for the design of tunnels and underground supporting structures in China.
- deep-buried tunnel /
- field measurement /
- deformation pressure /
- tunnel depth /
- tunnel span /
- grade of surrounding rock

HTML全文

0. 引言

锚固技术广泛应用于边坡加固、基坑支护以及公路交通等工程中^[1],各种新型锚固技术不断创新并应用到实际工程中。注浆支盘式锚杆是自主研发的一种新型锚固结构^[2],通过扩孔设备^[3]在局部位置形成支盘状空腔进行注浆,依靠锥状支盘的挤压扩张作用来增大端阻力,从而充分发挥和调用岩土体的强度达到提高抗拔承载力的目的。

对于新型锚固技术需要从锚杆的承载特性、锚固体交界面剪应力分布规律以及锚固机理等方面进行全面研究。杨坚等^[4]通过改变支盘式锚杆的埋置深度、支盘直径和双支盘间距等条件开展室内拉拔试验,从而探究支盘式锚杆的承载特性。而在锚固体荷载传递机制方面,目前已有诸多国内外学者取得了较好的成果。Nie等^[5]假设锚杆应力应变曲线服从线性强化弹塑性模型,分析锚杆拉拔荷载传递特征,但其假设锚固界面剪切滑移呈线性递增关系,未能考虑界面剪切滑移的塑性软化等特征。Ma等^[6]提出了一种非线性剪切滑移模型,推导出了锚固段轴力和界面剪应力解析解,并进行了锚杆锚固段荷载传递的全历程特征分析,Nemcik等^[7]在前者模型基础上通过 ${FLAC}^{2D}$ 软件模拟锚杆的剪切滑移关系进行了拉拔数值模拟。Ma等^[8]基于锚固界面剪切滑移呈三段线性的假设,分析了锚杆屈服前与屈服后全长锚固锚杆的荷载传递计算模型。Ni等^[9]推导了一种荷载-沉降响应的广义非线性软化荷载传递模型,该模型只需计算单参数n而且明确地考虑了表面摩擦在剪切位移作用下的硬化和退化行为。Bryson等^[10]基于现场收集的性能数据,提出岩石锚的荷载传递模型中包含峰值和残余荷载特性是至关重要的。陈建功等^[11]基于小波函数建立了能反映锚杆界面黏结-软化-滑动力学特性的非线性模型,并结合锚固体荷载传递的力学微分方程,推导出了锚杆拉拔荷载-位移曲线的解析解。黄明华等^[12]基于锚固界面的一种非线性剪切滑移模型,采用荷载传递方法分析了锚固长度对锚杆受力特性的影响,建立了锚杆临界锚固长度的计算方法。匡政等^[13]基于荷载传递理论及Kelvin位移解,推导出GFRP抗浮锚杆杆体与锚固体的剪应力、轴力沿锚固深度的荷载分布函数。但是,由于支盘式锚杆的结构为几何非线性（存在支盘状扩径体）,加上锚固段荷载传递涉及到锚固体和岩土材料的非均质性以及锚固体交界面的接触非线性等力学特性,故该新型锚杆的荷载传递机制尚不明确。

基于预制的支盘式锚杆为研究对象开展室内模型试验,分析锚杆荷载-位移曲线的变化特征和轴力传递特性,在此基础上探讨锚杆的力学模型,引入经典的VESIC扩孔理论推导支盘端阻力的理论计算方法,并结合锚杆-土交界面的双曲线剪应力-剪切位移非线性模型,采用分段变形协调迭代算法建立了荷载传递计算模型,并与室内试验实测数据进行对比验证,最后对支盘挤扩角、位移等参数进行敏感性分析。

1. 支盘式锚杆拉拔试验

1.1 试验装置

本次室内模型拉拔试验的装置包括模型箱、反力辅助设备、加载系统和应变采集系统。模型箱的尺寸为长×宽×高=3 m×2 m×1.5 m,其外框架由钢材焊接而成,两侧黏接15 mm厚的有机钢化玻璃,其表面标注锚杆埋置位置,如图1,2所示。反力辅助设备是自制的工字型钢反力架,反力架上设置的直角三棱柱（倾角15°）能够保证设备在安装及受拉过程中千斤顶油缸与锚杆在同一轴线上,提高测试精度。

图 1 室内模型示意图

Figure 1. Scheme of indoor model tests

下载: 全尺寸图片幻灯片

图 2 拉拔系统示意图

Figure 2. Scheme of drawing system

下载: 全尺寸图片幻灯片

加载系统设备是HCYL-60锚杆综合参数测定仪,包括位移传感器、负荷传感器、液压油缸和中文组态器等,如图2所示,其中负荷传感器受力产生电压信号,位移传感器随位移量的变化输出电阻信号,通过中文组态器中的20位A/D 转换器转换成数字信号,经单片机处理后由液晶显示器显示并存储压力值、位移值,其压力、位移最大测量范围分别为100 kN和50 mm,显示分辨率为0.01 kN和0.01 mm。

应变采集系统是武汉优泰生产的uT7110Y高速静态应变仪,测量量程达到±30000 με。应变片采用黄岩巨星电测生产的型号BX120-6AA（6X2）（电阻值（120±0.2）Ω,灵敏度2.08±1%）。

1.2 锚杆模型

支盘式锚固体系的施工工艺流程包括：①预成孔；②扩孔；③分序扩孔；④锚杆体安放；⑤注浆；⑥紧固锚杆。本试验为了便于进行大量的应变片粘贴、安装以及各类参数的测试,因此采用预制锚杆进行试验。

预制的支盘式锚杆由拉杆（HRB400螺纹钢筋）和浆材（PP42.5硅酸盐水泥、水）构成,其锚固浆液配比为水泥∶水=1∶0.4,将浆液倒入模具凝固养护而成,试验锚杆参数见表1。锚杆外侧周围分布有8对电阻应变片,从锚端起每隔18 cm设置一对应变片,其具体布置如图3所示。为确保应变片的存活率,在应变片表面刷上一层环氧树脂进行封装。

表 1 锚杆试验参数

Table 1. Test parameters for anchor

支盘直径R/mm	锚杆长度L/mm	锚杆直径d/mm	挤扩角θ/(°)	平行试验次数
—	1100	50	—	2
150	1100	50	60	2
200	1100	50	60	2
250	1100	50	60	2

下载: 导出CSV

| 显示表格

图 3 应变片分布图

Figure 3. Distribution of strain gauges

下载: 全尺寸图片幻灯片

1.3 边坡模型填筑

所用粉质黏土取自位于福州市晋安区的茉莉花园西北侧边坡,具体岩土物理力学性质参数见表2。边坡的制备是采用分层填筑压实的方法^[14],分层厚度为5 cm,通过打夯机压实次数控制土层的密实度。锚杆采用预埋的方式,当填土高度到达预设高度后,将锚杆模型轻置于黏土面上（与水平面呈15°）,并同钢化玻璃上的标注线保持平行,然后以放置锚杆的坡面为基准面,继续分层填筑黏土至达到预设埋置深度为止,每次试验均需重新填土夯实。

表 2 土的基本力学性质

Table 2. Basic mechanical properties of soil

密度ρ/(g·cm^-3)	含水率w/%	黏聚力c/kPa	内摩擦角 $φ$ /(°)	弹性模量E/(kN·m²)	泊松比 $ν$
1.90	25.2	21	15.6	7400	0.31

下载: 导出CSV

| 显示表格

1.4 试验过程

覆盖塑料膜养护48 h后进行拉拔试验,试验拉拔过程采用逐级加载法进行,相邻两级荷载的加载时间间隔至少为5 min,直至锚杆发生破坏。锚杆破坏的判定标准^[15]为：①单位荷载下锚固体位移的增量超过前一级荷载产生的单位荷载下位移增量的两倍；②锚固体位移不收敛,荷载不再增加。

1.5 试验结果分析

（1）荷载-位移曲线特征

图4为不同支盘直径下锚杆的荷载-位移曲线。从图4中可知：①相较于普通锚杆,支盘式锚杆的抗拔承载力得到了显著提高,分别达到了3.7,4.6,5.6倍（直径D由小到大）。挖出锚杆后发现其破坏模式是支盘达到极限强度导致破裂失效,最终发生剪切破坏,说明支盘端在拉拔的过程中承担了极大的荷载比重,也是该锚杆的显著特点之一,故支盘端的受力性状是研究重点。②结合Liu等^[16-18]基于数字照相测量技术对平板圆锚拉拔破坏机制的研究,可将荷载-位移曲线分为3个阶段（图4中的阶段划分以D = 250 mm为例）：第Ⅰ阶段为加载初期（基本为一直线）,当荷载较小,锚杆承载力主要由锚杆自重、锚固段与土的侧摩阻力承担,当位移增加到一定程度时,支盘端开始挤压端部的土体,小部分土体开始产生塑性屈服；第Ⅱ阶段为加载中期（为一曲线）,随着荷载与位移逐步增加,锚杆的锚固段在拉力作用下拉伸变形增大,支盘前端不断向周围土体挤压,使土体的塑性屈服区向外延伸,逐步进入极限屈服；第Ⅲ阶段为加载后期,随着荷载与位移持续增加,支盘前端土体屈服区向前发展,支盘达到极限强度发生破裂,锚杆与土体之间产生滑动,位移急剧增大,锚杆失效破坏。

图 4 荷载-位移曲线

Figure 4. Load-displacement curves

下载: 全尺寸图片幻灯片

（2）轴力传递特性

图5,6为支盘直径分别为200,250 mm锚杆沿轴向深度的轴力分布图。从图5,6中可知：①支盘式锚杆的轴力随着轴向深度的增加而逐渐减小,传递至底部时几乎为零,这是由于拉杆上的作用力传递到锚固段时逐步转化为侧摩阻力和盘阻力；②随着拉拔力增加,曲线分布的斜率也愈大,说明锚固段在拉力作用下拉伸变形逐步增大,剪应力不断增加；③在支盘的位置区域轴力发生了台阶状突变,说明支盘承担了较大的荷载,与上文所述破坏现象（见图4）相符,充分体现了支盘在抗拔方面的作用；④盘阻力随着荷载的增大而增大,荷载比（盘阻力与总荷载的比值）也是呈不断增长趋势（由30%至50%）。说明支盘随着位移增加而不断挤压前端土体,并逐步发挥其承载力作用,直至锚杆破坏（支盘达到极限强度）。

图 5 支盘式锚杆（D=200 mm）沿轴向深度的轴力分布图

Figure 5. Distribution of axial force of branch-type anchor (D = 200 mm) along axial depth

下载: 全尺寸图片幻灯片

图 6 支盘式锚杆（D=250 mm）沿轴向深度的轴力分布图

Figure 6. Distribution of axial force of branch-type anchor (D = 250 mm) along axial depth

下载: 全尺寸图片幻灯片

2. 支盘端阻力分析与计算

支盘式锚杆的力学模型如图7所示,锚杆的承载力主要有两部分组成：支盘端阻力与锚固段的侧阻力。在拉拔试验加载过程中,支盘与土体的相互作用与静力触探试验的锥头相似,土体挤压变形引起盘周围土体的径向位移,同时压缩挤密盘端土体（塑性区不断向外延伸发展）,提高了支盘式锚杆的抗拔承载力。

图 7 支盘式锚杆力学模型

Figure 7. Mechanical model for branch-type anchor

下载: 全尺寸图片幻灯片

因此,可借鉴静力触探试验锥阻力的计算方法,圆孔扩张理论提供了一种对锥状支盘端阻力的理论计算思路,它考虑了土的压缩性以及在拉拔过程中周围应力增加的影响,能够真实模拟锚杆支盘端在加载中的受力过程。

2.1 圆孔扩张理论

当支盘端位移达到最大值时S_max,达到极限扩孔压力p_u,根据经典VESIC^[19]扩孔理论可表示为

$p_{u} = c c o t φ (1 + \sin φ) {(\frac{R_{p}}{R_{u}})}^{\frac{2 \sin φ}{1 + \sin φ}} - c c o t φ$ ,

(1)

式中,R_p为支盘周围土体弹、塑性交界处的半径,R_u为支盘-土界面径向应力达到极限扩孔应力p_u时的圆孔半径,如图8所示,c为土体黏聚力, $φ$ 为土体内摩擦角。

图 8 圆孔扩张分析示意图

Figure 8. Analysis scheme of expansion cavity

下载: 全尺寸图片幻灯片

根据体积变形协调的关系,支盘位移产生的土体挤压变形等于弹性区加塑性区的土体体积变化,即

$π R_{u}^{2} - π R_{0}^{2} = π R_{p}^{2} - π(R_{p}^{} - u_{p}) + π(R_{p}^{2} - R_{u}^{2}) Δ$ ,

(2)

式中, $R_{0}^{}$ 为初始孔径, $u_{p}$ 为塑性区边界上的径向位移, $Δ$ 为塑性区土体平均体积应变,采用 $Δ = 0 .015$ ^[20-21]。

在弹塑性边界,当R=R_p时,则有 $σ_{r} = σ_{p}$ ,结合土体的几何方程、位移场、弹塑性应力场以及Mohr-Coulomb屈服准则,有

$1 + Δ - {(\frac{R_{0}}{R_{u}})}^{2} = {(\frac{R_{p}}{R_{u}})}^{2} Δ + 2 \frac{R_{p}}{R_{u}^{2}} \frac{1 + u}{E} c \cos φ$ ,

(3)

由式（3）整理可求得塑性区的半径R_p：

$R_{p} = R_{u} \sqrt{\frac{G (1 + Δ - {(\frac{R_{0}}{R_{u}})}^{2})}{c \cos φ + \frac{E Δ}{2 (1 + μ)}}}$ 。

(4)

把式（4）求得的R_p/R_u代入式（1）,可以求得支盘-土界面的径向挤压应力P_u：

$P_{u} = c cot φ (1 + sin φ) {[\frac{G (1 + Δ - {(\frac{R_{0}}{R_{u}})}^{2})}{c cos φ + G Δ}]}^{\frac{sin φ}{1 + sin φ}} - c cot φ$ ,

(5)

式中,G为土体的剪切模量, $G = E / (2 + 2 μ)$ 。

在上述对塑性区内体积变化影响的分析中,通常采用迭代法来修正塑性区平均体积应变 $Δ$ ^[22]：①假设初始的塑性区平均体积应变为 $Δ_{1}$ ,进而计算出塑性区的应力状态；②通过试验确定体积应变与应力的关系,确定修正的平均塑性体积应变 $Δ_{2}$ ；③设定 $| Δ_{n} - Δ_{n - 1} |$ 的预设值,重复①～②步骤直到数值收敛,从而得到塑性区平均体积应变。

2.2 挤土扩张过程分析

在锚杆拉拔过程中,随着荷载不断增加支盘端会逐步向前挤压土体,产生图9所示的虚线位移。

图 9 拉拔时支盘的位移示意图

Figure 9. Displacement scheme of branch plate during drawing

下载: 全尺寸图片幻灯片

在 $L_{pt}$ 段任意水平线上实线至虚线的距离即可表示为圆孔扩张的过程,在 $L_{1}$ 段中初始孔径是不变的,扩张后的孔径随着深度增加,故可表示为

$\frac{R_{0}}{R_{u}} = \frac{r_{0}}{r_{x}}$ ,

(6)

式中, $r_{0}$ 为锚杆直杆半径, $r_{x}$ 为支盘位移后轨迹线至中轴线的距离。

在L₂段中初始孔径与最终孔径是随深度不断变化的,但BC与B'C'是相互平行的,则说明这两条线之间的距离是始终不变的,故可表示为

$\frac{R_{0}}{R_{u}} = \frac{r_{x} - s_{m} \tan (θ)}{r_{x}}$ ,

(7)

式中, $s_{m}$ 为拉拔过程中产生的位移,θ为支盘挤扩角。

因此,将式（6）,（7）代入式（5）则可得到 $L_{pt}$ 段支盘接触面上不同深度的极限扩孔压力。

当 $r_{0} \leq r_{x} \leq (s_{m} \tan θ + r_{0})$ 时：

$P_{u} = c cot φ (1 + sin φ) {[\frac{G (1 + Δ - {(\frac{r_{0}}{r_{x}})}^{2})}{c cos φ + G Δ}]}^{\frac{sin φ}{1 + sin φ}} - c cot φ$ 。

(8)

当 $(s_{m} \tan θ + r_{0}) \leq r_{x} \leq R$ 时：

$\begin{matrix} P_{u} = c cot φ (1 + sin φ) {[\frac{G (1 + Δ - {(\frac{r_{x} - s_{m} \tan (θ)}{r_{x}})}^{2})}{c cos φ + G Δ}]}^{\frac{sin φ}{1 + sin φ}} \\ - c cot φ \end{matrix}$ 。

(9)

2.3 端阻力计算

取支盘局部段进行受力分析（见图10）,垂直和平行于盘-土界面的应力分别用 $σ_{n}$ 和 $τ_{n}$ 表示,假设土体服从莫尔库仑强度理论,

$τ_{n} = σ_{n} \tan φ + c$ 。

(10)

图 10 局部段支盘的受力分析示意图

Figure 10. Force analysis scheme of local branch plate

下载: 全尺寸图片幻灯片

利用支盘-土界面处应力的正交分解,可以得到垂直应力分量 $P_{v}$ ,该应力状态可以表示为

$P_{v} = P_{u} tan (φ + θ) + \frac{c \sec^{2} θ}{(1 - \tan φ \tan θ)}$ 。

(11)

从而对支盘上盘面进行积分,可以进一步得到盘端的估算阻力与支盘端的荷载传递函数,即

$Q = 2 \int_{0}^{L_{pt}} P_{v} r_{x} πd h = Q (s_{m})$ 。

(12)

3. 支盘式锚杆非线性变形分析

3.1 荷载传递函数

基于锚杆周围岩土体材料的非线性以及锚固交界面的剪应力 $τ_{s}$ 在拉拔过程中与剪切位移 $S$ 呈非线性^[23-25],即当剪切位移 $S$ 达到 $S_{ult}$ 时剪应力 $τ_{s}$ 会接近最大值 $τ_{f}$ 而不再增加。因此本文将采用双曲线模型（见图11）来描述支盘式锚杆界面-土之间非线性力学性状。

图 11 双曲线模型示意图

Figure 11. Scheme of hyperbolic model

下载: 全尺寸图片幻灯片

其表达式为

$τ_{s} = \frac{s}{a + b s}$ ,

(13)

式中, $a$ 表示锚-土界面初始刚度 $K$ 倒数,即曲线初始斜率； $b$ 表示锚-土界面极限剪应力值 $τ_{u}$ 的倒数。

$a = \frac{1}{K} = \frac{R_{0} \ln (\frac{R_{m}}{R_{0}})}{G}$ ,

(14)

式中, $K$ 为土层的初始刚度, $G$ 为土层的剪切模量, $R_{m}$ 为影响半径,根据Randolph等^[26]的研究成果, $R_{m} =$ $2.5 L ρ (1 - υ)$ ,其中,L为锚杆长度,ρ为非均质参数。

$b = \frac{1}{τ_{u}}$ ,

(15)

式中, $τ_{u}$ 为极限剪应力。

3.2 分段变形协调迭代算法

考虑如图12所示的支盘式锚杆荷载传递分析模型,将长为L的支盘式锚杆分为n个单元,每个单元长度为ΔL,其单元的受力形式分为直杆段与支盘段,二者的区别在于侧阻力的函数表达式和横截面面积 $A$ 不同。

图 12 支盘式锚杆荷载传递分析模型示意图

Figure 12. Scheme of load transfer analysis model for branch-type anchor

下载: 全尺寸图片幻灯片

基本假设：①锚杆的拉伸变形是线弹性的；②锚杆-土界面采用的是双曲线模型,界面不发生相对滑移,锚杆周围土体位移与锚固体位移一致；③单元的轴力是线弹性变化的。

本文的分段变形协调迭代算法的程序如图13流程图所示：

图 13 分段变形协调迭代算法流程图

Figure 13. Flow chart of sectional deformation coordinated iterative algorithm

下载: 全尺寸图片幻灯片

（1）给定锚杆一个初始拉拔力为 $P$ ,假定锚杆顶部产生了一个极小的位移 $S$ ,则锚杆上的首个单元 $N_{1}$ ： $P_{t 1} = P$ , $S_{t 1} = S$ 。

（2）假定单元N₁的平均拉拔力为 $P_{t 1}$ ,则 $N_{1}$ 的初始弹性变形为

$ε_{1} = \frac{P_{t 1} Δ L}{A E_{p}}$ 。

(16)

式中 $A$ 为锚杆直杆段的横截面面积, $A = π {(R_{0})}^{2}$ （若进入支盘端分段单元,则 $A = π {(R + R_{0})}^{2} / 2$ ）； $E_{p}$ 为锚杆的弹性模量, $E_{p} = [E_{a} A_{a} + E_{b} (A - A_{a})] / A$ ； $E_{a}$ 为杆体的弹性模量； $E_{b}$ 为注浆体的弹性模量； $A_{a}$ 为杆体的横截面面积。

（3）单元 $N_{1}$ 的中点位移：

$S_{c 1} = S_{t1} - \frac{ε_{1}}{2}$ 。

(17)

当 $N_{i} \leq N_{t}$ 或 $N_{i} > N_{t} + N_{p}$ 时,将 $S_{c i}$ 代入式（13）,得到 $N_{i}$ 段的剪应力 $τ_{s i}$ ,从而 $N_{i}$ 段的总侧摩阻力可表示为

$T_{i} = 2 π R_{0} Δ L τ_{s i}$ 。

(18)

反之则进入支盘段单元,将 $S_{c i}$ 代入式（11）,可得到 $N_{i}$ 段支盘端的侧应力 $P_{v}$ ,则 $N_{i}$ 段总侧摩阻力表示为

$T_{i} = 2 π r_{x} Δ L P_{v}$ 。

(19)

所以锚杆单元 $N_{1}$ 的底部荷载为 $P_{b 1} = P_{t 1} - T_{1}$ ,平均拉力为 $P_{1} = (P_{t 1} + P_{b 1}) / 2$ 。

（4）由此可以得到该段弹性变形的修正值 $ε^{'}$ ：

$ε^{'} = \frac{P_{1} Δ L}{E_{p} A}$ 。

(20)

（5）与初始弹性变量 $ε$ 进行收敛判定

$| ε^{'} - ε | \leq 10^{- 6}$ 。

(21)

若判定为假,则假定 $ε = ε^{'}$ ,重复（2）,（3）步骤,直到判定为真为止。反之则修正单元 $N_{1}$ 底部的位移为 $S_{b 1} = S_{t 1} - ε^{'} / 2$ ,并过度到下一个单元,将单元 $N_{1}$ 底部拉拔力和位移赋值于 $N_{2}$ 单元顶部拉拔力和位移。

（6）重复（2）～（5）步骤,依次计算完支盘式锚杆的n个单元,最后对总阻力 $T_{sum}$ 与拉拔力 $P$ 进行收敛判定(其单位为N),其中 $T_{sum} = \sum_{i = 1}^{n} T_{i}$ 。

$| T_{sum} - P | \leq 10$ 。

(22)

若为假,则修正 $S$ ,重复（2）～（5）步骤,直到拉拔力与总阻力平衡。反之则加大拉拔力 $P$ ,重复整个步骤,当达到锚杆的破坏判定（单位荷载下的锚固体位移增量超过前一级荷载产生的单位荷载下的位移增量2倍）,则结束整个程序,输出不同拉拔力 $P$ 所对应的支盘式锚杆的位移值 $S$ ,由此得到支盘锚杆的P-S曲线。

3.3 室内模型试验验证

根据拉拔力测试试验和试验材料性能试验,可确定出各材料参数取值：注浆体弹性模量为20 GPa,锚杆弹性模量为210 GPa,锚固体界面极限抗剪强度为68 kPa,土层参数与表2相同。将锚杆划分分为70个单元（支盘段为20个,直杆段为50个）,依照图13的计算流程图利用Matlab软件编写计算程序,再结合室内拉拔试验所得数据,对锚杆的荷载-位移曲线以及沿轴向深度的轴力分布特性进行对比分析,如图14,15所示。

图 14 拉拔曲线对比分析图

Figure 14. Comparative analysis of drawing curves

下载: 全尺寸图片幻灯片

图 15 支盘式锚杆（D=200 mm）轴力分布对比分析图

Figure 15. Comparative analysis of distribution of axial force of branch-type anchor (D=200 mm)

下载: 全尺寸图片幻灯片

从图14,15中可知,采用本文所述的理论模型计算得到的P-S曲线与各级荷载下的轴力值分布,跟室内试验实测值基本一致,说明了本文提出的模型有效描述了支盘式锚杆的荷载-位移曲线和轴力分布特征。

3.4 影响参数分析

（1）盘-土界面P_u分布及位移影响

图16为不同位移下盘-土界面的极限扩孔压力分布示意图。从图16中可以看出,支盘的极限扩孔压力从圆心沿着支盘的径向先快速非线性增长（图9中的 $L_{1}$ 段）,当达到峰值点后缓慢下降（图9中的 $L_{2}$ 段）。随着支盘前端土体的挤压位移量加大,其极限扩孔压力峰值也会不断增大,同时以支盘中心为圆点径向外移,而在 $L_{2}$ 段极限扩孔压力的减小幅度会逐步变缓,整体分布由“锥台状”发展为“圆柱状”。说明了支盘端对周围土体挤压越密实,极限扩孔压力越大且分布越均匀。

图 16 不同位移下盘-土界面的极限扩孔压力分布示意图

Figure 16. Scheme of ultimate pressure of cavity expansion at plate-soil interface under different displacements

下载: 全尺寸图片幻灯片

（2）支盘挤扩角的敏感性分析

图17为不同支盘挤扩角下的荷载位移曲线对比分析图。从图17中可以看出,支盘式锚杆的极限拉拔力随着支盘挤扩角的扩大而增大,当位移小于3.8 mm时,荷载位移曲线基本保持一致,当大于3.8 mm时,曲线斜率开始发生变化,挤扩角越大斜率越大,当挤扩角为25°～50°时为低敏感区,对承载性状影响不明显,斜率几乎不变,当挤扩角大于50°时,极限拉拔力显著提高。这是由于当位移大于3.8 mm时支盘端开始发挥承载力作用,挤扩角愈大盘端阻力愈大,从而越不易发生挤压位移,荷载位移曲线斜率也就越大。挤扩角对锚杆拉拔力学性状影响显著,但考虑到支盘需要足够的纵向高度保证其强度以及工程施工难度,所以建议 $θ$ 的取值范围为50°～65°之间,既能满足锚杆强度要求,也可大幅度提高抗拔承载力。

图 17 不同支盘挤扩角下的荷载位移曲线对比分析图

Figure 17. Comparative analysis of load-displacement curves under different squeezing angles

下载: 全尺寸图片幻灯片

4. 结论

（1）相较于普通锚杆,支盘式锚杆的抗拔承载力得到了显著提高,其荷载-位移曲线可划分为3个阶段,支盘端土体塑性区随荷载增加不断向外延伸发展,当支盘达到极限强度时破裂失效,位移急剧增加,轴力分布沿锚固体深度呈减小趋势且在支盘处发生台阶状突变,荷载比也不断增长。

（2）在试验模型基础上探讨支盘式锚杆力学模型,通过借鉴静力触探试验锥阻力的计算方法,引用经典的VESIC扩孔理论模拟支盘拉拔过程中的挤土效应,推导出了盘端阻力的计算理论公式。

（3）锚-土界面选用双曲线非线性模型,采用分段变形协调迭代算法建立了荷载传递计算模型,计算所得结果与实测数据基本一致。

（4）支盘的极限扩孔压力沿支盘径向先快速非线性增长后缓慢减小,其峰值随着位移增大而增大并以支盘中心为圆点径向外移,整体分布由“锥台状”发展为“圆柱状”。支盘的挤扩角大于50°时,对锚杆承载性状影响显著。

图 1 岩石应力-应变曲线及围岩破坏过程

Figure 1. Stress-strain curve of rock and failure process of surrounding rock

下载: 全尺寸图片幻灯片

图 2 支护力与围岩特性关系分析图

Figure 2. Relationship between supporting force and characteristics of surrounding rock

下载: 全尺寸图片幻灯片

图 3 多点位移计布置及现场安装图

Figure 3. Layout of multi-point displacement meter and field installation

下载: 全尺寸图片幻灯片

图 4 塑性区范围

Figure 4. Range of plastic zone

下载: 全尺寸图片幻灯片

图 5 各围岩级别样本占比

Figure 5. Percentage of grades of surrounding rock

下载: 全尺寸图片幻灯片

图 6 各隧道跨度样本占比

Figure 6. Percentage of tunnel spans

下载: 全尺寸图片幻灯片

图 7 各施工工法样本占比

Figure 7. Percentage of construction methods

下载: 全尺寸图片幻灯片

图 8 Ⅳ级围岩形变压力分布图

Figure 8. Deformation load distribution of surrounding rock Grade IV

下载: 全尺寸图片幻灯片

图 9 Ⅴ级围岩形变压力分布图

Figure 9. Deformation load distribution of Grade V surrounding rock

下载: 全尺寸图片幻灯片

图 10 形变压力测点布置及荷载模式分布图

Figure 10. Arrangement of measuring points for deformation load and loading mode

下载: 全尺寸图片幻灯片

图 11 竖向形变压力变化曲线

Figure 11. Curve of vertical deformation loads

下载: 全尺寸图片幻灯片

图 12 竖向形变压力变化曲线（幂函数）

Figure 12. Curve of vertical deformation loads (power function)

下载: 全尺寸图片幻灯片

图 13 竖向形变压力变化曲线(线性函数)

Figure 13. Curve of vertical deformation loads (linear function)

下载: 全尺寸图片幻灯片

图 14 竖向形变压力变化曲线(Ⅳ)

Figure 14. Curve of vertical deformation loads (Ⅳ)

下载: 全尺寸图片幻灯片

图 15 竖向形变压力变化曲线(Ⅴ)

Figure 15. Curve of vertical deformation loads (Ⅴ)

下载: 全尺寸图片幻灯片

图 16 竖向形变压力变化曲线(Ⅳ)

Figure 16. Curve of vertical deformation loads (Ⅳ)

下载: 全尺寸图片幻灯片

图 17 竖向形变压力变化曲线(Ⅴ)

Figure 17. Curve of vertical deformation loads (Ⅴ)

下载: 全尺寸图片幻灯片

图 18 竖向形变压力变化曲线(Ⅳ)

Figure 18. Curve of vertical deformation loads (IV)

下载: 全尺寸图片幻灯片

图 19 竖向形变压力变化曲线(Ⅴ)

Figure 19. Curve of vertical deformation loads (Ⅴ)

下载: 全尺寸图片幻灯片

图 20 IV级围岩侧压力系数分布图

Figure 20. Lateral pressure coefficient of surrounding rock Grade IV

下载: 全尺寸图片幻灯片

图 21 V级围岩侧压力系数分布图

Figure 21. Lateral pressure coefficient of surrounding rock Grade V

下载: 全尺寸图片幻灯片

图 22 形变压力随围岩级别变化曲线

Figure 22. Curves of deformation loads with grades of surrounding rock

下载: 全尺寸图片幻灯片

图 23 形变压力随隧道埋深变化曲线

Figure 23. Curves of deformation loads with tunnel depth

下载: 全尺寸图片幻灯片

图 24 形变压力随隧道跨度变化曲线

Figure 24. Curves of deformation loads with tunnel span

下载: 全尺寸图片幻灯片

表 1 形变压力随围岩级别拟合公式

Table 1 Fitting formula for deformation loads with grade of surrounding rock

形变压力	拟合公式	相关系数
竖向形变压力	$q = 15 {.86e}^{0. 51 s}$	0.961

下载: 导出CSV

表 2 形变压力随隧道跨度拟合公式

Table 2 Fitting formula for deformation loads with tunnel span

形变压力拟合公式	拟合公式	相关系数
幂函数	$q = 3.0 B^{1.40}$	0.91
线性函数	$q = 12.0 B - 45.25$	0.94

下载: 导出CSV

表 3 形变压力随隧道埋深拟合公式

Table 3 Fitting formula for deformation loads with tunnel depth

形变压力	拟合公式	相关系数
（IV）竖向形变压力	$q / B = 2. 81 \cdot H^{0.2 0}$	0.89
（V）竖向形变压力	$q / B = 4. 35 \cdot H^{0. 22}$	0.86

下载: 导出CSV

表 4 形变压力随隧道埋深拟合公式

Table 4 Fitting formula for deformation loads with tunnel depth

形变压力	拟合公式	相关系数
（IV）竖向形变压力	$q / B = 10.38 \cdot e^{\frac{- 35.2}{H}}$	0.80
（V）竖向形变压力	$q / B = 16.63 \cdot e^{\frac{- 32.03}{H}}$	0.83

下载: 导出CSV

表 5 各围岩级别正态分布函数表

Table 5 Normal distribution function for various grades of surrounding rock

围岩级别	数学期望值	标准差	概率密度函数
IV	0.47	0.10	$f (x) = \frac{1}{\sqrt{0.2 π}} \exp (- \frac{(λ - 0.47)}{0.02})$
V	0.71	0.12	$f (x) = \frac{1}{\sqrt{0.24 π}} \exp (- \frac{(λ - 0.71)}{0.0288})$

下载: 导出CSV

表 6 侧压力系数分布区间

Table 6 Distribution intervals of lateral pressure coefficient

数据区间	侧压力系数（ $λ$ ）（IV级围岩）	概率/%	侧压力系数（ $λ$ ）（V级围岩）	概率/%
$(μ - σ, μ + σ)$	（0.37,0.57）	68.26	（0.59,0.83）	68.26
$(μ - 1.96 σ, μ +1 .96 σ)$	（0.27,0.67）	95.44	（0.47,0.95）	95.44
$(μ - 1.96 σ, μ +1 .96 σ)$	（0.21,0.73）	99.73	（0.40,1.02）	99.73

下载: 导出CSV

表 7 侧压力系数表

Table 7 List of lateral pressure coefficient

围岩级别	侧压力系数	置信区间/%	《隧规》侧压力系数
Ш	<0.25	—	<0.15
Ⅳ	0.25~0.5	88.5	0.15~0.30
Ⅴ	0.5~1.0	91.0	0.30~0.50

下载: 导出CSV

表 8 竖向形变压力计算公式

Table 8 Formulas for calculating vertical deformation loads

序号	公式	考虑因素	修正系数 $α$	相关系数
1	q=α∙B^1.25∙H^0.1∙e^0.5s	s,B,H	$α$ =0.35	0.897
2	q=α∙B^1.25∙e^(0.5s-15/H)	s,B,H	$α$ =0.71	0.920
3	q=α∙B^1.25∙e^0.45s	s,B	$α$ =0.81	0.920
4	q=α∙B∙H^0.1∙e^0.5s	s,B,H	$α$ =0.69	0.907
5	q=α∙B∙e^(0.5s-15/H)	s,B,H	$α$ =1.38	0.935
6	q=α∙B∙e^0.45s	s,B	$α$ =1.55	0.915
注：s为围岩级别,B为隧道跨度,H为隧道埋深。

下载: 导出CSV

表 9 侧压力系数表

Table 9 List of lateral pressure coefficient

围岩级别	Ⅲ	Ⅳ	Ⅴ
侧压力系数	<0.25	0.25~0.5	0.5~1.0

下载: 导出CSV

参考文献(23)

[1]	关宝树. 隧道工程设计要点集[M]. 北京: 人民交通出版社, 2003. GUAN Bao-shu. Gist of tunnel engineering design[M]. Beijing: China Communications Press, 2003. (in Chinese)
[2]	王毅才. 隧道工程[M]. 北京: 人民交通出版社, 2003. WANG Yi-cai. Tunnel Engineering[M]. Beijing: China Communications Press, 2003. (in Chinese)
[3]	吴海科, 刘远明. 围岩压力中形变压力的影响因素[J]. 建筑技术开发, 2016, 43(1): 163-164. https://www.cnki.com.cn/Article/CJFDTOTAL-JZKF201601090.htm WU Hai-ke, LIU Yuan-ming. Factors rock pressure in deformation pressure[J]. Building Technology Development, 2016, 43(1): 163-164. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JZKF201601090.htm
[4]	铁路隧道设计规范：TB10003—2016[S]. 北京: 中国铁道出版社, 2016. Code for Design of Railway Tunnel: TB10003—2016[S]. Beijing: China Railway Publishing House, 2016. (in Chinese)
[5]	李鹏飞, 赵勇, 张顶立, 等. 基于现场实测数据统计的隧道围岩压力分布规律研究[J]. 岩石力学与工程学报, 2013, 32(7): 1392-1399. doi: 10.3969/j.issn.1000-6915.2013.07.014 LI Peng-fei, ZHAO Yong, ZHANG Ding-li, et al. Study of distribution laws of tunnel surrounding rock pressure based on field measured data statistics[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(7): 1392-1399. (in Chinese) doi: 10.3969/j.issn.1000-6915.2013.07.014
[6]	李鹏飞, 周烨, 伍冬. 隧道围岩压力计算方法及其适用范围[J]. 中国铁道科学, 2013, 34(6): 55-60. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGTK201306010.htm LI Peng-fei, ZHOU Ye, WU Dong. Calculation methods for surrounding rock pressure and application scopes[J]. China Railway Science, 2013, 34(6): 22-60. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGTK201306010.htm
[7]	GOODMAN R E. Introduction to rock mechanics[M]. New York: Wiley, 1978.
[8]	曾钱帮, 王恩志, 王思敬. 深埋圆形硐室围岩塑性形变压力的广义Hoek- Brown破坏准则解及其三元非线性回归模型[J]. 岩石力学与工程学报, 2009, 28(增刊2): 3543-3549. ZENG Qian-bang, WANG En-zhi, WANG Si-jing. Analytical solution for plastic deformation pressure around a deep circular opening based on generalized Hoek-Brown failure criterion and its ternary nonlinear regression model[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(S2): 3543-3549. (in Chinese)
[9]	付国彬. 巷道围岩破裂范围与位移的新研究[J]. 煤炭学报, 1995, 20(3): 304-310. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB503.015.htm FU Guo-bin. Recent investigation of extent of fractured zone and displacement of rocks around the roadways[J]. Journal of China Coal Society, 1995, 20(3): 304-310. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB503.015.htm
[10]	PARK K H, KIM Y J. Analytical solution for a circular opening in an elastic-brittle-plastic rock[J]. International Journal of Rock Mechanics and Mining Sciences, 2006, 43(4): 616-622. doi: 10.1016/j.ijrmms.2005.11.004
[11]	GUAN Z, JIANG Y, TANABASI Y. Ground reaction analyses in conventional tunnelling excavation[J]. Tunnelling and Underground Space Technology, 2007, 22(2): 230-237. doi: 10.1016/j.tust.2006.06.004
[12]	LEE Y K, PIETRUSZCZAK S. A new numerical procedure for elasto-plastic analysis of a circular opening excavated in a strain-softening rock mass[J]. Tunnelling and Underground Space Technology, 2008, 23(5): 588-599. doi: 10.1016/j.tust.2007.11.002
[13]	公路隧道设计规范：JTG D70—2004[S]. 北京: 人民交通出版社, 2004. Code for Design of Road Tunnel: JTG D70—2004[S]. Beijing: China Communications Press, 2004. (in Chinese)
[14]	伍冬. 山岭隧道围岩压力计算方法及其适用性研究[D]. 北京: 北京交通大学, 2012. WU Dong. Study on Surrounding Rock Pressure Calculation Method and Application of Mountain Tunnel[D]. Beijing: Beijing Jiaotong University, 2012(in Chinese)
[15]	刘学增, 叶康. 山岭公路隧道围岩压力统计规律分析[J]. 岩土工程学报, 2011, 33(6): 890-895. LIU Xue-zeng, YE Kang. Statistical analysis of surrounding rock pressure of mountain road tunnels[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(6): 890-895. (in Chinese)
[16]	陈进. 岩体结构渐进破坏的理论与实践[D]. 徐州: 中国矿业大学, 1988. CHEN Jin. Theory and Practice of Progressive Failure of Rock Mass Structure[D]. Xuzhou: China University of Mining and Technology, 1988. (in Chinese)
[17]	JIANG Yu-jing. Theoretical and Experimental Study on the Stability of Deep Undergrounding[D]. Fukuoka-ken: Kyushu University, 1993.
[18]	何满潮, 景海河, 孙晓明. 软岩工程力学[M]. 北京: 科学出版社, 2002. HE Man-chao, JING Hai-he, SUN Xiao-ming. Soft Rock Engineering Mechanics[M]. Beijing: Science Press, 2002. (in Chinese)
[19]	柳厚祥, 方风华. 预埋式多点位移计现场确定围岩松动圈的方法研究[J]. 矿冶工程, 2006, 26(1): 1-4. LIU Hou-xiang, FANG Feng-hua. Study on the method to Define the wall rock loose zone with embedded displacement meter[J]. Mining and Metallurgical Engineering, 2006, 26(1): 1-4. (in Chinese)
[20]	乔春生, 张清. 锚杆轴力分布与软弱岩体中隧道塑性区的关系[J]. 铁道学报, 1999, 21(2): 72-75. QIAO Chun-sheng, ZHANG Qing. Study on the method to define the wall rocks loose zone with embedded displacement mete[J]. Journal of the China Railway Society, 1999, 21(2): 72-75. (in Chinese)
[21]	杜国文. 小浪底工程多点位移计测值突变原因分析[J]. 东北水利水电, 1998(6): 26-29. DU Guo-wen. Study on the method to define the wall rocks loose zone with embedded displacement meter[J]. Water Resourse & Hydropower of Northeast, 1998(6): 26-29. (in Chinese)
[22]	张顶立, 陈立平. 隧道围岩的复合结构特性及其荷载效应[J]. 岩石力学与工程学报, 2016, 35(3): 456-469. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201603003.htm ZHANG Ding-li, CHEN Li-ping. Compound structural characteristics and load effect of tunnel surrounding rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(3): 456-469. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201603003.htm
[23]	张顶立, 孙振宇. 复杂隧道围岩结构稳定性及其控制[J]. 水力发电学报, 2018, 37(2): 1-11. https://www.cnki.com.cn/Article/CJFDTOTAL-SFXB201802001.htm ZHANF Ding-li, SUN Zhen-yu. Structural stability of complex tunnel surrounding rock and its control[J]. Journal of Hydroelectric Engineering, 2018, 37(2): 1-11. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SFXB201802001.htm

施引文献(7)

期刊类型引用(5)

1.	黄炜，简文彬，杨坚，豆红强，罗金妹. 新型多支盘土锚的承载性能研究. 工程地质学报. 2025(02): 772-782 . 百度学术
2.	谷复光，王军，常虹. 注浆支盘式锚杆承载特性数值模拟分析. 低温建筑技术. 2024(07): 134-138 . 百度学术
3.	孙世国，谢远东，宋腾飞. 锚杆扩大头不同受力方向对其承载能力影响机制探究. 北京工业职业技术学院学报. 2023(01): 1-4 . 百度学术
4.	黄炜，简文彬，杨坚，豆红强，罗金妹. 多支盘锚杆的原型试验与荷载传递特征分析. 岩土力学. 2023(02): 520-530 . 百度学术
5.	董旭光，李峥，崔自治，李宏波，王永胜. 一种边坡真空锚管研发及工作机理分析. 铁道学报. 2022(12): 114-124 . 百度学术