振冲桩法碎石原料与实际成桩体抗剪强度差异试验研究

江洎洧; 王小华; 饶锡保; 丁高俊; 张雪杨; 黄帅

doi:10.11779/CJGE2020S1031

振冲桩法碎石原料与实际成桩体抗剪强度差异试验研究 English Version

1.
长江水利委员会长江科学院水利部岩土力学与工程重点实验室,湖北　武汉　430010
2.
长江勘测规划设计研究有限责任公司,湖北　武汉　430010

基金项目:

国家自然科学基金项目 51979009

中央级公益性科研院所基本科研业务费资金项目 CKSF2019191/YT

详细信息

作者简介:
江洎洧（1984— ）,男,博士,高级工程师,主要从事岩土工程数值计算及粗粒土的试验研究。E-mail：jiangjw1023@163.com

中图分类号: TU411
计量
- 文章访问数: 185
- HTML全文浏览量: 15
- PDF下载量: 73
出版历程
- 收稿日期: 2020-06-03
- 网络出版日期: 2022-12-07
- 刊出日期: 2020-10-31

Experimental study on difference of shear strength between gravel materials and actual vibroflotation piles

1.
Changjiang River Scientific Research Institute, Key Laboratory of Geotechnical Mechanics and Engineering of the Ministry of Water Resources, Wuhan 430010, China
2.
Changjiang Institute of Survey Planning Design and Research, Wuhan 430010, China

摘要

摘要: 振冲桩法通过碎石等透水性材料对软弱地基进行置换挤密而形成复合地基,以往对该复合地基抗剪强度关注相对较少。由于碎石桩和桩周软土存在挤密和置换作用,碎石原料与振冲形成的实际桩体不论颗粒组成还是变形受力特性都存在差异,基于工程现场对成桩体的检测和取样,在室内开展土工试验,分析碎石原料与实际成桩体在物理力学参数方面的差异,试验表明：①实际成桩体中的细颗粒含量显著高于碎石原料,其主要来源于振冲施工过程中桩周软土的混（挤）入；②成桩体最大干密度高于碎石原料,与碎石原料级配不良及细颗粒混入有直接关系；③以现场实测压实度备样分别对碎石原料和两种桩间距成桩体开展大型三轴试验,成桩体整体强度较碎石原料降低约10%,考虑其在应力变形特性上的差异,若以5%应变对应偏应力求取强度参数,成桩体较碎石料强度减小超过20%,桩周土混入对桩体抗剪强度及应力变形特性均有较显著的劣化影响。结合初步研究结论,提出了复合地基强度设计参数的修正建议。
- 振冲碎石桩 /
- 碎石料 /
- 级配特性 /
- 强度及应变 /
- 试验研究
Abstract: The vibroflotation pile method is used to replace and compact the soft foundation with permeable materials such as gravel to form composite foundation, and in the past, the attention to shear strength of this kind composite foundation is little. Due to the compaction and replacement of gravel and soil around piles, difference is shown between gravel materials and actual vibroflotation gravel piles, regardless of particle composition or stress-strain characteristics. Based on field detection, sampling for vibroflotation gravel piles and gravel materials and indoor geotechnical tests to analyze the difference of physical and mechanical properties between vibroflotation gravel piles and gravel materials, the researches show: (1) The content of fine particles in actual piles is significantly higher than that of gravel materials, which mainly comes from the mixing (squeezing) of soft soil around piles during vibroflotation process; (2) The maximum dry density of piles is higher than that of gravel materials, which is directly related to the poor grading of gravel materials and mixing of fine particles; (3) Refer to actual compaction degree in the field, large-scale triaxial tests are carried out on gravel materials and two kinds of piles with different pile spacings, the strength of piles is about 10% lower than that of gravel materials, considering the difference in stress-strain characteristics, if the strength parameter is taken according to 5% axial strain, the strength of piles is more than 20% lower than that of gravel materials, and mixing of soil around piles has a significant deterioration effect on the shear strength and also the stress-strain characteristics of the piles. Based on the above initial conclusions, the modified suggestions for strength design parameters of composite foundation are proposed.
- vibroflotation gravel pile /
- gravel soil /
- particle gradation /
- strength and strain /
- experimental study

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. Gravel materials of vibroflotation gravel piles

下载: 全尺寸图片幻灯片

图 2 振冲碎石桩散体状态

Figure 2. Scattered state of vibroflotation gravel pile

下载: 全尺寸图片幻灯片

图 3 不同试样的颗粒级配曲线

Figure 3. Gradation curves for different kinds of samples

下载: 全尺寸图片幻灯片

图 4 大型三轴压缩试验仪

Figure 4. Large-scale triaxial compression apparatus

下载: 全尺寸图片幻灯片

图 5 桩间距1.5 m成桩体的应力应变关系曲线

Figure 5. Stress-strain curves of vibroflotation gravel pile with 1.5 m spacing

下载: 全尺寸图片幻灯片

图 6 碎石料的应力应变关系曲线

Figure 6. Stress-strain curves of gravel materials

下载: 全尺寸图片幻灯片

图 7 5%轴向应变对应的莫尔圆及强度包络线

Figure 7. Mohr's circle and linear fitting curves (ε₁=5%)

下载: 全尺寸图片幻灯片

表 1 研究试样颗粒级配统计

Table 1 Statistics of particle size distribution for research samples (%)

材料	粒径区间/mm
材料	100~150	60~100	40~60	20~40	10~20	5~10	2~5	1~2	0.5~1	0.25~0.5	0.1~0.25	<0.1
料场碎石	31.44	40.75	17.72	6.31	0.98	0.43	0.24	0.22	0.16	0.54	0.90	0.29
桩间距2.0 m桩体	27.90	28.92	9.36	10.67	7.83	6.03	2.27	1.26	0.78	2.23	1.85	0.91
桩间距1.5 m桩体	30.54	30.15	10.09	10.60	5.54	3.82	1.59	1.32	0.61	1.92	2.99	0.84

下载: 导出CSV

表 2 成桩材料及碎石料密度测试统计表

Table 2 Densities of vibroflotation gravel piles and gravel materials

材料	现场干密度/(g·cm^-3)	最大干密度/(g·cm^-3)	密实度
2.0 m间距桩体	2.026	2.249	0.901
1.5 m间距桩体	2.031	2.264	0.897
料场碎石	—	2.130	—

下载: 导出CSV

表 3 大型三轴压缩试验成果表

Table 3 Test results of large-scale triaxial compression tests

试样名称	干密度/(g·cm^-3)	抗剪强度指标
试样名称	干密度/(g·cm^-3)	c/ kPa	φ/(°)
桩间距2.0 m桩体	2.026	36	32.8
桩间距1.5 m桩体	2.031	52	31.7
料场碎石	1.92	42	36.2

下载: 导出CSV

表 4 复合地基抗剪强度计算值对比

Table 4 Comparison of shear strength for composite foundation

材料	1.5 m间距复合地基		2.0 m间距复合地基
材料	φ/(°)	c/ kPa	φ/(°)	c/ kPa
成桩体	26.49	7.76	19.42	10.05
碎石原料	31.17	7.76	22.66	10.05

下载: 导出CSV

参考文献(12)

[1]	李进元. 振冲碎石桩法地基处理在阴坪水电站中的应用[J]. 岩石力学与工程学报, 2013, 32(增刊1): 2968-2976. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2013S1049.htm LI Jin-yuan. Application of foundation treatment with vibro-replacement stone piles to yinping hydropower station[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(S1): 2968-2976. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2013S1049.htm
[2]	陈健, 苏岩松, 张杰, 等. 碎石桩软土复合地基整体抗剪强度研究[J]. 水运工程, 2016(5): 146-148, 161. https://www.cnki.com.cn/Article/CJFDTOTAL-SYGC201605027.htm CHEN Jian, SU Yan-song, ZHANG Jie, et al. Overall shear strength of stone column-soil composite ground[J]. Port & Waterway Engineering, 2016(5): 146-148, 161. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SYGC201605027.htm
[3]	RAYAMAJHI D, ASHFORD S A, BOULANGER R W, et al. Dense granular columns in liquefiable ground: I shear reinforcement and cyclic stress ratio reduction[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2016: 04016023.
[4]	彭芝平, 杨军, 陈耀光, 等. 饱和盐渍土碎石排水桩加强夯试验研究[J]. 岩土工程学报, 2010, 32(增刊2): 136-141. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2010S2035.htm PENG Zhi-ping, YANG Jun, CHEN Yao-guang, et al. Experimental study on ground improvement of saturated saline soil with dynamic consolidation and drainage stone columns[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(S2): 136-141. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2010S2035.htm
[5]	王德咏, 梁小丛, 牛犇. 振冲密实法处理吹填地基的两个关键技术问题[J]. 水运工程, 2019(6): 163-166. https://www.cnki.com.cn/Article/CJFDTOTAL-SYGC201906029.htm WANG De-yong, LIANG Xiao-cong, NIU Ben, et al. Two key technical issues of vibroflotation compaction for hydraulic-filled foundation treatment[J]. Port & Waterway Engineering, 2019(6): 163-166. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SYGC201906029.htm
[6]	欧阳芳. 包裹碎石桩复合地基静动力特性及其设计方法研究[D]. 成都: 西南交通大学, 2017. OU Yang-fang. Studies on the Behavior and Design Method for Geosynthetic Encased Stone Columns Composite Ground Under Static and Dynamic Conditions[D]. Chengdu: Southwest Jiaotong Unversity, 2017. (in Chinese)
[7]	孙立强, 邵丹丹, 冯守中, 等. 土工合成材料约束碎石桩承载特性研究[J]. 岩土工程学报, 2019, 41(增刊2): 29-32. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2019S2009.htm SUN Li-qiang, SHAO Dan-dan, FENG Shou-zhong, et al. Bearing capacity of geosynthetics-enclosed stone columns[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(S2): 29-32. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2019S2009.htm
[8]	周建安. 碎石桩模型三轴试验及其成果应用分析[J]. 地质灾害与环境保护, 1997, 8(2): 40-45. https://www.cnki.com.cn/Article/CJFDTOTAL-DZHB702.006.htm ZHOU Jian-an. Analysis and application of triaxial test results of crushed-stone pile model[J]. Journal of Geological Hazards and Environment Preservation, 1997, 8(2): 40-45. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DZHB702.006.htm
[9]	SIAHAAN F, INDRARATNA B, NGO N T, et al. Influence of particle gradation and shape on the performance of stone columns in soft clay[J]. Geotechnical Testing Journal, 2018, 41(6): 20160234.
[10]	建筑地基处理规范：JGJ 79—2012[S]. 2012. Technical code for ground treatment of buildings: JGJ 79—2012[S]. 2011. (in Chinese)
[11]	水电水利工程振冲法地基处理技术规范：DL/T 5214—2016[S]. 2016. Technical specification for vibroflotation ground trwatment of hydroelectric and water resources projects: DL/T 5214—2016[S]. 2011. (in Chinese)
[12]	龚晓南. 地基处理手册[M]. 3版.北京: 中国建筑工业出版社, 2008. GONG Xiao-nan. Foundation Treatment Manual[M]. 3rd ed. Beijing: China Architecture & Building Press, 2008. (in Chinese)

施引文献(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 . 百度学术