Experimental study on failure characteristics of pressured under-reamed anchors under horizontal loads
-
摘要: 通过自行研制的模型试验装置,利用PIV技术研究水平荷载作用下压力型扩大头锚杆力学特性以及周围土体位移发展机理。试验结果表明,压力型扩大头锚杆的受力过程可分为侧摩阻力阶段、过渡阶段和端部阻力阶段,在每个阶段转变之间荷载位移曲线会出现“拐点”。基于数字图像分析,提出压力型扩大头锚杆扩大头锚固段近端土体临界位移场影响范围呈“灯泡型”,且当扩大头锚固段长径比超过3∶1后,随着继续增大长径比对扩大头锚固段近端的位移场范围基本无影响。Abstract: Through a self-developed model test device and PIV technology, the mechanical characteristics of pressured under-reamed anchors under horizontal loads and the development mechanism of the surrounding soil displacement are studied. The test results show that the stress process of the pressured under-reamed anchors can be divided into lateral friction stage, transition stage and end resistance stage, and there will be an "inflection point" in the load-displacement curve between the transitions of each stage. According to the image analysis, it is proposed that the influence range of critical displacement field of soil around the under-reamed anchorage segment of pressured anchor is "bulb-shaped". When the length-diameter ratio of the enlarged head anchoring section exceeds 3∶1, the displacement field at the proximal end of the under-reamed anchoring section is basically unaffected with the continuous increase of length-diameter ratio.
-
0. 引言
膨胀土具有较强的亲水性,对于处于季节冻土区的渠道工程、边坡工程等,反复冻融循环对膨胀土的劣化作用会极大影响工程稳定性[1-2]。已有的研究表明,冻融作用对土体的性质改变相对复杂,部分特征并无统一认识[3-4],典型的如内摩擦角,变大、变小、基本不变的结论均有。由于冻融循环作用下,不同土体响应机制的复杂性,针对膨胀土特定土体的冻融分析在工程实践中具有更好的指导意义。
目前针对膨胀土的冻融循环研究已有了较大进展,蔡正银等[2]以高寒渠道膨胀土为研究对象,在干湿冻融耦合的试验条件下系统讨论了冻融作用下膨胀土的裂隙演化规律。Lu等[5]、Tang等[6]则讨论了不同状态下的膨胀土与循环冻融次数相关的劣化规律,研究表明,随着冻融循环次数的增加,膨胀土的破坏强度、弹性模量、抗剪强度等力学指标劣化加剧,且第一次冻融循环下其劣化最为明显。
目前有关冻融循环下的损伤模型,更多是基于连续介质理论展开。损伤力学中的损伤因子D部分研究需要借助其他测试手段进行确定,如卢再华等[7]通过CT扫描技术定义损伤因子,其构建的损伤模型包含23个参数,均需测试确定。目前主流的损伤因子定义是在已有的应力应变测试成果上进行,其中假设土体微元满足weibull分布,基于应变等价原理开展损伤模型构建获得了广泛应用。杨明辉等[8]、李向东等[9]进行了基质吸力和外部荷载耦合状态下的非饱和土损伤模型构建。Tang等[6]基于多种技术手段综合研究了损伤程度与冻融循环次数之间的关系,并构建了冻融循环作用下的统一损伤模型。总结目前的土的损伤模型,其损伤定义时,主要考虑了弹性模量的退化,核心均集中到应力应变曲线的模拟,但是对于土体本身的性质如压硬性等考虑不足。
传统弹性力学中的应变是微应变,但是土工试验中的土体变形量较高,因此完全按照传统应变去研究其土体性质无法充分考虑土体的内在特质。基于此,通过引入自然应变的概念、weibull统计分布规律、应变等价原理等,推导构建适用于膨胀土的损伤模型,并以损伤模型为基础进行冻融损伤分析等工作。
1. 基于自然应变的冻融损伤模型
1.1 自然应变适用性
以图1中所示的岩土体为分析模型,其满足均质各向同性的基本假定,受荷时,假定无新的损伤发生,其应力应变关系用HOOK定律进行描述,
dσ=Kdε, (1) 式中,
σ 为岩土材料截面应力(压为正),K为岩土体的刚度,ε 为当前应力状态下的应变量。根据弹性力学理论,岩土体的应变可为
dε=−dH/H, (2) 式中,H为当前应力状态下的试件高度。
式(2)计算的应变,Freed将其定义为自然应变或真实应变[10]。具体到应用阶段通过下式计算岩土体应变量更为广泛:
dε=−dH/H0, (3) 式中,H0为未受荷时的试件高度。对于高刚度材料,通常变形量较原来的形状较为微小,数量级相差在2个以上,此时式(2),(3)基本等效。式(3)由于更简洁的表达形式,在岩土工程建设中具有广泛的应用场景。但对于部分软岩、松散土体、断裂结构、多孔介质等,全部采用式(3)去分析其力学特征,则可能会阻碍对岩土本质特征的正确认识。
对于土体而言,土颗粒是高刚度材料,由于空气和水的存在,其形成的多介质土体则不完全满足式(3)的适用条件。尤其是对于亲水性矿物含量高的膨胀土而言,饱和时,会极大的削弱土体的刚度,干湿循环、冻融循环的条件下,会进一步导致其刚度的降低。因此在分析其力学特征时,式(2)具有更好的适用性。
1.2 基于应变等价原理的冻融损伤变量
岩土体的损伤产生主要有受荷和环境影响2部分。针对产生损伤岩土体,目前常用的损伤变量定义方法是基于应变等价原理实现,即可将损伤状态的材料在有效应力作用下的应变与无损状态下的应变进行等效,其无损状态的材料所受的应力称为名义应力。原理图如图1所示,表达式为
ε=σE′=σ′E, (4) σA=σ′A′, (5) 式中,
σ 为损伤状态下的有效应力,A为损伤状态下有效承载面积,E'为损伤状态下的岩土体模量,σ'等效无损状态下的名义应力,A'等效无损状态下有效承载面积,E为等效无损状态下的岩土体模量。损伤变量可用岩土体模量来表示,
D=1−E'/E。 (6) 更一般的情况下,岩土体的损伤可用岩土体变形模量的损伤来定义,目前多数文献均采用此定义[6, 8-9]。
基于式(2)的自然应变概念,可进一步建立岩土体刚度和传统方法变形模量的关系,
E=σ1−exp(−σ/K)。 (7) 将式(7)代入式(6),可得自然应变下的岩土体损伤变量为
D=1−K'/K, (8) 式中,K'为损伤状态下的土体刚度。式(6),(8)的形式基本一致,但是代表的是两种不同应变概念下的损伤定义方式。对于高刚度材料,E=K,传统的应变定义是自然应变定义下的一种特殊情况。
对于环境损伤Df(冻融循环、干湿循环)和受荷损伤Ds耦合的情况,总的损伤变量D可表示为
D=1−(1−Df)(1−Ds)。 (9) 1.3 受荷损伤模型构建
假设土体内部微元服从概率分布,可用Weibull分布定量描述微元强度的分布特征[11],
φ(F)=mF0(FF0)m−1e−(FF0)m=dDsdF, (10) 式中,φ(F)为土体微元强度概率分布函数,F0,m为Weibull分布的参数,F为微元的强度。
对式(10)进行积分可得,受荷损伤的损伤变量Ds与土体微元强度F的关系,
Ds=1−exp[−(F/F0)m]。 (11) 土体微元的强度准则采用D-P准则进行计算[12]。土体微元的强度F为
F=αI′1+√J′2, (12) 式中,
α 为材料参数,I′1 为有效应力状态下的应力第一不变量,J2'为有效应力状态下的偏应力第二不变量,α=2sinφ√3(3−sinφ) ,I′1=σ′1+2σ′3 ,J′2=13(σ′1−σ′3)2 ,} (13) 其中,φ土体微元内摩擦角,
σ′1 ,σ′3 为常规三轴条件下的竖向有效应力和有效围压。通过自然应变的概念,有效应力作用下的最大主应力方向的应变
ε′1 为ε′1=1−exp(−σ′1−σ′3K')(σ′1−σ′3)(σ′1−2μσ′3), (14) 式中,K'为土体的刚度,μ为泊松比。常规三轴试验时的固结过程,测得应力应变曲线主要是偏应力
(σ1−σ3) 和除去初始应变之后的应变ε1 ,基于应变等价原理,将有效应力进行转换,可将式(14)改写为σ1−σ3=−K0(1−Ds)ln(1−ε1), (15) 式中,K0为偏应力为0时的初始刚度。设孔隙水压力为PW,则联立式(12)~(15),土体基于自然应变的损伤本构模型如下:
σ1−σ3=−K0ln(1−ε1)exp[−(FF0)m], (16) 其中,
F={2sinφ√3(3−sinφ)[(σ1−σ3)+3(σ3−PW)]+1√3(σ1−σ3)}⋅−K0ln(1−ε1)σ1−σ3。 (17) 2. 基于遗传算法的模型参数确定
从式(16),(17)中可以发现,对于损伤参量无法直接测量的参数主要有初始刚度K0、土体微元内摩擦角φ、Weibull分布的参数F0和m。由于土体结构的复杂性,为避免误差累积,降低模型本身的精度,具体参数确定时采用智能算法遗传算法实现。遗传算法本质上属于生物进化算法,通过模拟生物进化中的优胜劣汰,对于具体的问题给出最优解[13]。对于给定的应力应变曲线,限定参数的取值范围,通过数值计算软件MATLAB编程,进行参数自动搜索,全局寻优,具体原理如图2所示。
3. 损伤模型验证
3.1 冻融循环试验数据
为进一步研究冻融循环后,膨胀土的损伤特性,利用Tang等[6]。的试验数据进行进一步的分析。试验土样充分饱水后进行冻融循环,-10℃保持12 h后,10℃保持12 h为1次冻融循环。试验采用常规三轴压缩试验,试验条件为固结排水,因此可不考虑超静孔隙水压力影响,试验围压为100 kPa。文献[6]采用的循环冻融次数NFT分别为0,1,3,7,11。试验结果如图3所示。
3.2 损伤模型验证
通过MATLAB编程,遗传算法可实现损伤模型参数的自动获取。通过实测的应变应变数据模拟,得到到实测值和模拟值的对比如图4所示,模拟参数如表1所示。从图4中可明显的发现,基于自然应变推导的损伤模型可以很好的模拟膨胀土在不同冻融循环下的应力应变关系,模拟精度较高。
表 1 损伤模型拟合参数Table 1. Fitting parameters for damage modelNFT K0/MPa φ/(°) F0 m 0 15.96 35.31 154.05 0.3174 1 14.89 35.83 153.06 0.3205 3 12.96 34.18 146.31 0.3200 7 12.42 35.30 145.52 0.3202 11 11.90 34.86 143.63 0.3209 3.3 损伤模型参数分析
从表1可以明显发现,4个参数中初始刚度K0、参数F0呈现明显的劣化特征,土体微元内摩擦角φ,m相对变化不大。Weibull分布中的参数F0和m通常代表着宏观平均强度和均布系数。由于土体结构本身的劣化,对土体形成一定的损伤,导致F0的减小。由于冻融过程中,土颗粒的团聚效应,均布系数m有略微增加的趋势,验证了论文中本构模型对于分析冻融损伤的合理性。冻融损伤产生的一个重要原因就是孔隙的增加,结构劣化,因此膨胀土的冻融损伤同样表现在初始刚度K0的明显降低。
选取受冻融影响明显的初始刚度K0和参数F0作为分析对象,绘制其与冻融循环次数的关系曲线如图5所示。从图5中可以明显看出,初始刚度K0和参数F0的劣化规律均满足一定的负指数关系,随着冻融次数的增加,劣化速度减慢,逐渐趋于一稳定值。
4. 讨论
通过建立的自然应变损伤模型,高精度的模拟了膨胀土冻融循环后的应力应变曲线,以及分析了冻融损伤对模型参数的影响。为进一步研究冻融损伤的损伤演化过程,利用式(9)绘制损伤因子D的变化曲线如图6所示。通过图6可以看出,在应力应变曲线的后期,不同冻融循环下的损伤因子基本上趋于一致。在应力应变曲线前期,随着冻融循环次数的增加,冻融损伤速度也进一步加快。但是这种加快的趋势随循环次数的增加开始减弱。
事实上,基于自然应变的膨胀土损伤模型优势不仅在于模拟精度上的提升,而且可以解决传统应变模式下对于土体本身性质认识的局限性。通过自然应变损伤模型,将非损伤变形和损伤变形进行一定区分,绘制图7的应力变形曲线。
通过图7可以发现,非损伤变形部分呈现逐渐密实的状态,即对应加载过程中不存在土体损伤时的变形状态,变形曲线特征基本上与土体的侧限压缩试验对应。三轴压缩时,由于损伤变形贡献了主要的变形量,因此曲线特征基本上与损伤变形曲线一致。
5. 结论
通过引入自然应变的概念,基于Weibull统计理论和应变等价原理,构建了损伤模型。在模型验证的基础上进一步分析膨胀土的冻融损伤特征。通过研究得出以下结论。
(1)构建的自然应变损伤模型可高精度的模拟膨胀土不同冻融循环下的应力应变曲线特征。遗传算法在模型验证中自动化程度高,可实现参数的自动寻优。
(2)膨胀土初始刚度K0和参数F0的冻融劣化明显,与冻融循环次数呈负指数关系。内摩擦角呈现出一定的波动下降趋势,均布系数m有略微增加的趋势。
(3)从不同冻融循环下的损伤因子演化曲线可以看出,冻融损伤主要发生在应力应变曲线前期,且随着冻融循环次数增加,应力应变曲线前期的损伤速度有加快趋势。应力应变后期损伤过程则基本一致。
(4)绘制非损伤变形和损伤变形,非损伤变形呈现逐渐压密的状态,损伤变形贡献了更多的变形量。其变形特征与侧限压缩试验曲线和三轴压缩曲线具有较好的对应性。
-
表 1 标准砂物理力学参数表
Table 1 Physical and mechanical parameters of standard sand
重度 γ /(kN·m-3)密度 ρ /(g·cm-3)含水率/% 黏聚力c/kPa 内摩擦角 φ /(°)土粒相对密度 Gs 不均匀系数 Cu 17.25 1.76 0.037 0 33.2 2.67 4.77 表 2 扩大头锚杆物理力学参数表
Table 2 Physical and mechanical parameters of under-reamed anchors
类型 重度 γ /(kN·m-3) 弹性模量E/GPa 抗压强度/MPa 泊松比 ν 扩大头锚固段 木块 6.9 7.5 24 0.47 橡胶 12.5 4.4 — 0.49 钢螺纹杆 79.3 206 — 0.30 表 3 锚固段对锚杆承载力的影响
Table 3 Effects of water content on bearing capacity of anchors
类型 锚固段长度/cm 第一阶段承载力/N 第二阶段承载力/N 第三阶段承载力/N 木制 10 160.15 206.72 296.61 15 116.09 196.66 323.93 20 163.77 232.29 331.34 橡胶 10 61.02 198.90 321.48 15 76.30 200.93 367.66 20 81.21 215.98 350.31 表 4 临界位移场影响范围统计
Table 4 Influence ranges of critical displacement field
类别 长度/cm D1/D D2/D D3/D β1 /(°)β2 /(°)木制锚固段 10 3.6 3.0 2.1 34 28 15 4.3 3.0 2.0 35 28 20 4.3 3.1 2.1 25 21 橡胶锚固段 10 4.0 2.5 2.4 35 34 15 3.8 2.8 1.9 37 34 20 3.9 2.8 2.1 32 28 -
[1] 贾金青. 深基坑预应力锚杆柔性支护法的理论及实践[M]. 北京: 中国建筑工业出版社, 2014. JIA Jing-qing. Theory and Practice of Flexible Support Method of Pre-Stressed Anchor Rod in Deep Foundation Pit[M]. Beijing: China Architecture and Building Press, 2014. (in Chinese)
[2] 张乐文, 汪稔. 岩土锚固理论研究之现状[J]. 岩土力学, 2002, 23(5): 627-631. ZHANG Le-wen, WANG Ren. Research on status quo of anchorage theory of rock and soil[J]. Rock and Soil Mechanics, 2002, 23(5): 627-631. (in Chinese)
[3] 尤春安. 锚固系统应力传递机理理论及应用研究[D]. 青岛: 山东科技大学, 2004. YOU Chun-an. Study on the Theory and Application of Stress Transfer Mechanism of Anchorage System[D]. Qingdao: Shandong University of Science and Technology, 2004. (in Chinese)
[4] 程良奎. 岩土锚固研究与新进展[J]. 岩石力学与工程学报, 2005, 24(21): 3803-3810. CHENG Liang-kui. Research and new progress in ground anchorage[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(21): 3803-3810. (in Chinese)
[5] 程良奎, 胡建林, 张培文. 岩土锚固技术新发展[J]. 工业建筑, 2010, 40(1): 98-101. CHENG Liang-kui, HU Jian-lin, ZHANG Pei-wen. The new development of ground anchorage in China[J]. Industrial Construction, 2010, 40(1): 98-101. (in Chinese)
[6] 尤志嘉, 付厚利, 尤春安, 等. 土层锚固体应力传递机制[J]. 岩土力学, 2018, 39(1): 85-92. YOU Zhi-jia, FU Hou-li, YOU Chun-an, et al. Stress transfer mechanism of soil anchor body[J]. Rock and Soil Mechanics, 2018, 39(1): 85-92. (in Chinese)
[7] 孙志峰. 高精度水平锚索技术研究及在三峡工程中的应用[D]. 长春: 吉林大学, 2006. SUN Zhi-feng. 2006. Research on Technology of Horizontal Cables with High Accuracy and its Application in the Three Gorges Project[D]. Changchun: Jilin University. (in Chinese)
[8] 张永兴, 卢黎, 饶枭宇, 等. 压力型锚杆力学性能模型试验研究[J]. 岩土力学, 2010, 37(7): 2045-2050. ZHANG Yong-xing, LU Li, RAO Xiao-yu, et al. Model test research on mechanical behavior of compression type-rock bolt[J]. Rock and Soil Mechanics, 2010, 37(7): 2045-2050. (in Chinese)
[9] 贺建清, 陈清, 陈秋南. 压力型与拉力型锚杆工作性状的室内足尺模型对比试验研究[J]. 水文地质工程地质, 2013, 40(4): 50-55. HE Jian-qing, CHEN Qing, CHEN Qiu-nan. A comparative study of work characteristics of pressure-type and tensile-type anchors based on indoor full-sized model tests[J]. Hydrogeology & Engineering Geology, 2013, 40(4): 50-55. (in Chinese)
[10] LI Jian, LI Jie. Numerical analysis on a new pressure-type anchor cable with precast anchor head based on FLAC3D[C]//ACSR-Advances in Computer Science Research, 2016, 67: 237-241.
[11] 夏柏如, 方铭东. 一种扩体头型土锚的试验[J]. 现代地质, 1997, 11(1): 91-94. XIA Bo-ru, FANG Ming-dong. Test of one under-reamed ground anchor[J]. Modern Geology, 1997, 11(1): 91-94. (in Chinese)
[12] 曾庆义, 杨晓阳, 杨昌亚. 扩大头锚杆的力学机制和计算方法[J]. 岩土力学, 2010, 31(5): 1359-1367. ZENG Qing-yi, YANG Xiao-yang, YANG Chang-ya. Mechanical mechanism and calculation method of bid expanded anchor rods[J]. Rock and Soil Mechanics, 2010, 31(5): 1359-1367. (in Chinese)
[13] 张慧乐, 刘钟, 徐龚鑫, 等. 扩体锚杆尺寸效应模型试验研究[J]. 岩土工程学报, 2010, 32(增刊1): 468-472. ZHANG Hui-le, LIU Zhong, XU Gong-xin, et al. Model tests on size effect of under-reamed ground anchors[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(S1): 468-472. (in Chinese)
[14] 郭钢, 刘钟, 邓益兵, 等. 砂土中扩体锚杆承载特性模型试验研究[J]. 岩土力学, 2012, 33(12): 3645-3652. GUO Gang, LIU Zhong, DENG Yi-bing, et al. Model test research on bearing capacity characteristics of underreamed ground anchor in sand[J]. Rock and Soil Mechanics, 2012, 33(12): 3645-3652. (in Chinese)
[15] LIU J, LIU M, ZHU Z. Sand Deformation around an uplift plate anchor[J]. Journal of Geotechnical Geoenvironmental Engineering, 2012, 138(6): 728-737.
[16] 王哲, 王乔坎, 马少俊. 扩大头可回收预应力锚索极限抗拔力计算方法研究[J]. 岩土力学, 2018, 39(增刊2): 202-208. WANG Zhe, WANG Qiao-kan, MA Shao-jun. A method for calculating ultimate pullout force of recoverable under-reamed prestressed anchor cable[J]. Rock and Soil Mechanics, 2018, 39(S2): 201-208. (in Chinese)
-
期刊类型引用(5)
1. 谢柏涵,孔元元,张学飞,孙东彦,王智慧. 冻融循环作用下盐渍土损伤模型的建立及验证. 人民长江. 2025(02): 183-189+200 . 百度学术
2. 孙超,谷鸿雁,郭浩天. 冻融循环次数对粉质黏土强度劣化影响的试验研究. 吉林建筑大学学报. 2024(05): 66-71 . 百度学术
3. 吕思清,朱杰兵,汪斌,祝永锁. 冻融荷载耦合作用下含开口裂隙砂岩宏细观损伤模型研究. 岩石力学与工程学报. 2023(05): 1124-1135 . 百度学术
4. 李新宇,凌贤长,曲娜. 考虑温度效应的冻结膨胀土统计损伤模型. 吉林大学学报(工学版). 2023(08): 2339-2349 . 百度学术
5. 张雁,薛晓雨,杨晓蕴. 低温状态下纤维加筋膨胀土数值模拟试验. 中国科技论文. 2022(12): 1351-1357 . 百度学术
其他类型引用(5)