高速公路改扩建工程边坡抗滑桩加固效果数值分析

罗根传，何忠明，曾铃



高速公路改扩建工程边坡抗滑桩加固效果数值分析

罗根传1, 2，何忠明1，曾铃3

(1. 长沙理工大学交通运输工程学院，湖南长沙，410004；2. 广西壮族自治区交通投资集团有限公司，广西南宁，530028；3. 长沙理工大学土木与建筑工程学院，湖南长沙，410004)

为了分析改扩建工程中使用抗滑桩对边坡加固治理时的支护稳定性，基于FLAC3D软件对某改扩建边坡利用抗滑桩加固效果进行数值模拟，得到边坡的应力应变云图、桩身前后岩体压力及桩自身内力的分布规律。研究结果表明：抗滑桩的阻挡作用使桩周岩体局部应力增加，抗滑桩处的剪应变增量增大；抗滑桩前、后的岩土压力呈三角形分布，边坡滑裂面以下岩土压力分布复杂，可近似认为属于梯形或矩形分布；抗滑桩桩体在滑裂面位置处的剪力和弯矩最大。

改扩建；抗滑桩；桩岩相互作用机理；加固效果

随着我国基础设施建设的大规模启动，边坡与滑坡工程治理引起了众多学者的重视。抗滑桩加固因具有支挡效果好、施工期短、施工灵活等优点，成为治理边坡及滑坡工程的常见手段之一[1−3]。而且随着原 有高速公路改扩建工程的不断增加，采用抗滑桩提高改扩建边坡稳定性的方法也逐渐增多，但在此过程中会出现许多关于抗滑桩支护稳定性的问题，如支护桩开挖扰动对边坡土体容易造成影响等[4]。近年来，一些研究者对抗滑桩的问题进行了研究，如：吴鹏等[5]通过分析抗滑桩的荷载位移曲线变化规律，提出了考虑桩土滑移动态调整的分析计算方法；高长胜[6]通过试验手段对稳定计算中土体强度指标、边坡变形破坏发生机制及抗滑桩与土体相互作用进行了研究，指出采用双排抗滑桩的加固效果明显好于相对应单排桩的加固效果；朱彦鹏等[7]提出集中弹簧模型建立支护桩与土相互作用的−曲线，推导出考虑非线性共同作用的预应力锚杆支护桩的挠曲微分方程；吕庆等[8]采用平面应变有限元分析模型，对桩周土应力应变变化规律进行研究，解释了土拱作用的力学机理。以上研究均对抗滑桩受力情况或桩土共同作用产生的土拱效应进行了研究，但对于桩−岩之间的耦合效应特别是基于数值分析方法来反映桩−岩之间相互作用的研究尚少。为此，本文作者利用FLAC-3D数值分析方法，对抗滑桩桩−岩之间的相互作用机理及加固效果进行分析，以便为工程实践提供参考。

1 桩−岩相互作用机理

在FLAC-3D中，桩与实体单元之间的相互作用通过耦合弹簧实现[9−12]。耦合弹簧为非线性、可滑动的连接体，能够在桩身节点和实体单元之间传递力和弯矩。桩力学机理模型如图1所示。桩单元包括3个弹簧单元：弹簧1和2为切向弹簧，可以模拟接触面之间的相对位移；弹簧3为法向弹簧，可以模拟法向荷载的作用及桩周围岩土对桩身的挤压作用。

图1 桩力学机理模型

桩岩之间的接触面单元可以通过接触面节点和实体单元之间建立联系[13−15]。接触面法向方向所受到的力由目标面方向决定。在每个时步中，首先得到接触面节点和目标面之间的绝对法向刺入量和相对剪切速度，再利用接触面本构模型计算法向力和切向力。处于弹性阶段时，+Δ时刻接触面的法向应力和切向应力通过下式计算：

n(t+Δt)=nn+n(1)

si(t+Δt)=si(t)+sΔsi(t+0.5Δt)+si(2)

式中：n(t+Δt)为法向矢量；si(t+Δt)为切向力矢量；n为接触面单元灌入目标面的绝对位移；Δsi为相对剪切位移矢量；n为初始化造成的附加法向应力；si为初始化造成的附加切向应力；n和s为接触面单元的法向和切向刚度；为单元节点代表的面积。

桩−岩之间的剪切行为如图2所示。当桩−砂浆、砂浆−岩石界面发生剪切滑移时，利用桩单元法向平均有效围压m、砂浆的黏结强度g、砂浆的剪切刚度sa、砂浆的内摩擦角g、砂浆的外周长g，可得到相应的力学性质。m根据下式得到：

图2 桩变形模型

式中：b1和b2为桩轴线的法平面应力。假设相应接触面符合Mohr-Coulomb破坏模型，桩的轴向剪力可表示为

式中：sa为砂浆体内产生的剪切力；bs为桩和岩土界面之间的相对位移；为桩长。

桩的横向抗剪能力ns为

ns=nsbn(5)

式中：ns为桩的横向剪切刚度；bn为桩的横向位移。

2 计算模型的建立

泉州—南宁高速公路柳州(鹿寨)—南宁段改扩建工程K1369+400右侧高边坡，基岩为灰白色厚层灰岩，中—微风化，岩石坚硬，岩层产状为320°∠21°，节理较发育。原有边坡坡高28.8 m，分为3级，每级坡比为1:0.75，根据边坡地质条件分析，改扩建拟采用垂直开挖坡脚方案：先对老路底部进行抗滑桩预加固，然后，拆除桩前边坡圬工防护、自上而下逐排开挖并逐排张拉锚索，与抗滑桩形成受力整体。抗滑桩长14.0 m，直径为2.0 m，桩体采用C30混凝土，边坡岩体和抗滑桩的物理力学参数如表1所示。利用ANSYS有限元软件，根据实际地形及地层坡面土进行三维数值建模，再将数值模型导入FLAC3D中，设置桩与岩体的黏结力和内摩擦角与所接触的岩体相 同，二者的切向刚度和法向刚度均为1.3×1011N/m。模型底部置刚性约束，两侧设置水平约束，模型材料选用Mohr−Coulomb弹塑性材料，屈服准则采用莫 尔−库仑准则，以最大不平衡力比率趋近于10−5为计算收敛值[16]，模型网格划分共生成节点13 629个，单元11 590个，计算模型如图3所示。

表1 材料物理力学参数

图3 数值计算模型

3 计算结果分析

3.1 无抗滑桩支护结果分析

图4和图5所示分别为无抗滑桩条件下的水平等值云图与剪切应变等量云图，采用强度折减法计算得到边坡在无抗滑桩支护时的安全系数为1.08。从图4可以看出：当无抗滑桩支护时，边坡上部的水平位移整体较大，最大值出现在开挖后的坡脚处，达到22.9 cm，表现出由坡脚向上逐渐减小的趋势。从图5可见：剪切应变增量较大值在边坡内部已具有连续贯通的趋势，且此时边坡的安全系数为1.08，小于规范所要求的最低值，可以推断此时边坡处于较危险的状态，需要及时进行加固处理。

图4 无支护时的水平位移等值云图

图5 无支护时的剪切应变增量云图

3.2 抗滑桩支护效果分析

按照模型中所设计的抗滑桩支护方案对边坡进行支护，计算所得的水平及竖向位移等值云图见图6与图7。从图6可以看出：边坡的最大水平位移出现在垂直开挖的坡脚附近，最大值为1.45 mm；抗滑桩施加后，有效抑制了边坡土体的水平位移，其与边坡形成有机整体，利用自身的强度和刚度阻止了边坡上部的持续滑移。从图7可以看出：边坡的最大竖向位移出现在桩前开挖坡体处，最大值为12.01 mm，方向竖直向上。这是由于抗滑桩的作用抑制了桩后岩体的变形，使得桩后岩体应力沿桩向下传递至新开挖路面处，造成新开挖路面呈现向上隆起的现象，这一现象通过速度矢量图(图8)得到更好体现。

图6 有支护时的水平位移等值云图

图7 有支护时的竖向位移等值云图

图8 有支护时的速度矢量图

开挖完成后边坡的剪切应变增量云图如图9所示。从图9可见：与未支护时相比，剪切应变增量分布形状不再与坡面平行，剪力通过桩底传至行车道，这大大降低了边坡滑动的风险。由于抗滑桩的阻挡作用，致使桩周岩体局部应力增加，从而导致其剪应变也相应增加，抗滑桩处的剪应变增量也较大。

图9 有支护时的剪切应变增量云图

图10所示为抗滑桩桩前桩后岩体压力分布图。从图10可见：在抗滑桩滑裂面以上，桩前桩后岩体压力大致呈三角形分布；而在滑裂面以下，岩体压力变得复杂，但可近似取为矩形或梯形分布。桩后滑坡推力和桩前岩体抗力分布及其合力作用点受滑动面形状、部位、滑坡的类型、岩层性质、抗滑桩变形等综合因素的影响，所以，在具体分析时，应根据现场监测情况确定岩体压力分布，而数值计算可作为一种验证 手段。

1—桩前岩体压力；2—桩后岩体压力

通过对安全系数进行反算，计算抗滑桩桩体的有效推力(有效推力为桩后推力与桩前抗力之差，可通过滑裂面以上桩长范围内桩前桩后土体水平向应力的积分获得[17])，进而得到桩体自身的剪力和弯矩，分别如图11和图12所示。从图11和图12可以看出：在桩身约1/3高度处，剪力发生了突变，弯矩达到最大值，这与剪切应变增量云图相符；而位于基岩中桩体的剪力和弯矩变化不明显，因此，可以在桩体上部施加锚索结构形成桩锚体系，以进一步增加支护性能，从而也可减小桩体本身内力。

图11 桩体剪力分布图

图12 桩体弯矩分布图

4 结论

1) 抗滑桩用于支档边坡后，由于抗滑桩的阻挡作用，桩周岩体局部应力增加，从而导致其剪应变增大，且抗滑桩处的剪应变增量也较大。

2) 抗滑桩桩前、后岩体压力大致呈三角形分布，而在滑裂面以下，岩体压力可近似取为矩形或梯形分布，显得较复杂。桩后滑坡推力和桩前岩体抗力分布受滑动面形状、部位、滑坡的类型、岩层性质、抗滑桩变形等综合因素的影响。

3) 抗滑桩桩体在滑裂面位置处的剪力和弯矩最大，因此，可在相应位置施加锚索结构形成桩锚体系，以进一步增强支档性能，保证边坡安全稳定。

[1] 杨智. 预应力锚索抗滑桩在治理公路滑坡中的应用[D]. 长沙: 中南大学土木工程学院, 2011: 22−23. YANG Zhi. Anti-slide piles with prestressed anchor cables in the control of application in highway landslide[D]. Changsha: Central South University School of Civil Engineering, 2011: 22−23.

[2] 徐爱民, 柳群义, 朱自强, 等. 层状岩体边坡抗滑桩加固效应的数值分析[J]. 中南大学学报(自然科学版), 2011, 42(8): 2453−2458. XU Aimin, LIU Qunyi, ZHU Ziqiang, et al. Numerical analysis for reinforcement response of pile in stratified rock slope[J]. Journal of Central South University (Science and Technology), 2011, 42(8): 2453−2458.

[3] 傅强, 孙飞飞, 唐承铁, 等.双排微型抗滑桩的抗滑效果[J].中南大学学报(自然科学版), 2013, 44(4): 1596−1602. FU Qiang, SUN Feifei, TANG Chengtie, et al. Effect of double-row anti-slide micropiles[J]. Journal of Central South University (Science and Technology), 2013, 44(4): 1596−1602.

[4] 谭积青. 改扩建高速公路桩锚支护边坡稳定性数值分析[J]. 中外公路, 2013, 33(3): 34−36. TAN Jiqing. Reconstruction and expansion of highway slope stability of pile-anchor supporting numerical analysis[J]. Journal of China & Foreign Highway, 2013, 33(3): 34−36.

[5] 吴鹏, 龚维明, 梁书亭.考虑桩土滑移和整体刚度动态调整的群桩分析方法[J].岩土工程学报, 2007, 29(2): 225−230.WU Peng, GONG Weiming, LIANG Shuting. Analytic method of pile groups considering slipping between piles and soil and dynamic adjustment of integral stiffness[J]. Chinese Journal of Geotechnical Engineering, 2007, 29(2): 225−230.

[6] 高长胜. 边坡变形破坏及抗滑桩与土体相互作用研究[D]. 南京: 南京水利科学研究院, 2007: 17−98. GAO Changsheng. Study on slope deformation and interaction between reinforced piles and soils[D]. Nanjing:Nanjing Hydraulic Research Institute, 2007: 17−98.

[7] 朱彦鹏, 魏升华. 深基坑支护桩与土相互作用的研究[J]. 岩土力学, 2010, 31(9): 2840−2841. ZHU Yanpeng, WEI Shenghua. Research on interaction between deep excavation supporting pile and soil[J]. Rock and Soil Mechanics, 2010, 31(9): 2840−2841.

[8] 吕庆, 孙红月, 尚岳全. 抗滑桩桩后土拱效应的作用机理及发育规律[J]. 水利学报, 2010, 41(4): 471−474. LÜ Qing, SUN Hongyue, SHANG Yuequan. Anti-slide piles been arch effect mechanism and growth rule[J]. Journal of Hydraulic Engineering, 2010, 41(4): 471−474.

[9] Grasselli G, Giovanni. 3D behaviour of bolted rock joints: Experimental and numerical study[J]. International Journal of Rock Mechanics & Mining Sciences, 2005, 42(1): 13−24.

[10] Jeong S, Kim B, Won J, et al. Uncoupled analysis of stabilizing piles in weathered slopes[J]. Computers and Geotechnics, 2003, 30(8): 671−682.

[11] 谭捍华. 抗滑桩预加固边坡的能量分析方法[J]. 岩土力学, 2011, 32(S2): 190−197. TAN Hanhua. Energy analysis method for pre-reinforcing slopes with anti-slide piles[J]. Rock and Soil Mechanics, 2011, 32(S2): 190−197.

[12] 陈育民, 徐鼎平. FLAC/FLAC3D基础与工程实例[M]. 北京: 中国水利水电出版社, 2009: 160−163. CHEN Yumin, XU Dingping. FLAC3D foundation and project examples[M]. Beijing: China Water Power Press, 2009: 160−163.

[13] 李邵军, 陈静. 边坡桩−土相互作用的土拱力学模型与桩间距问题[J]. 岩土力学, 2010, 31(5): 1352−1358.LI Shaojun, CHEN Jing. Mechanical model of soil arch for interaction of piles and slope and problem of pile spacing[J]. Rock and Soil Mechanics, 2010, 31(5): 1352−1358.

[14] 吴坤铭, 王建国, 谭晓慧. 基于可靠度分析确定刚性抗滑桩锚固深度[J]. 岩土工程学报, 2012, 34(2): 237−242. WU Kunming, WANG Jianguo, TAN Xiaohui. Determination of anchorage depth for rigid anti-slide piles based on reliability analysis[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(2): 237−242.

[15] 林杭, 曹平. 锚杆长度对边坡稳定性影响的数值分析[J]. 岩土工程学报, 2009, 31(3): 470−474.LIN Hang, CAO Ping. Numerical analysis for the effect of cable length to the stability of slope[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(3): 470−474.

[16] 孙书伟, 林杭, 任连伟. FLAC3D在岩土工程中的应用[M]. 北京: 中国水利水电出版社, 2011: 271−279. SUN Suwei, LIN Hang, REN Lianwei. FLAC3D application in geotechnical engineering[M]. Beijing: China Water Power Press, 2011: 271−279.

[17] Chow Y K. Analysis of piles used for slope stabilization[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1996, 20(9): 635−646.

[18] 郑颖人, 赵尚毅. 用有限元强度折减法求滑(边)坡支裆结构的内力[J]. 岩石力学与工程学报, 2004, 23(20): 3552−3558. ZHENG Yingren, ZHAO Shangyi. Calculation of inner force of support structure for landslide/slope by using strength reduction FEM[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(20): 3552−3558.

(编辑 陈灿华)

Numerical analysis of reinforcement effect of slope anti-slide pile on expressway reconstruction projection

LUO Genchuan1, 2, HE Zhongming1, ZENG Ling3

(1. School of Traffic and Transportation Engineering,Changsha University of Science &Technology, Changsha 410004, China;2. Communications Investment Group Co. Ltd of Guangxi Zhuang Autonomous Region, Nanning 530028, China;3. School of Civil Engineering and Architecture, Changsha University of Science &Technology, Changsha 410004, China)

In order to analyze the stability in the expansion project when the anti-slide plie was used on slope governance, by using numerical simulation to analyze the effect of strengthening of anti-slide plie on a expansion project which was based on FLAC3D, the stress strain contours of the slope, body stress of rock before and after pile and its internal force distribution were obtained. The results show that the anti-slide pile blocking effect makes the local stress of rock mass increase, and anti-slide pile in shear strain increment also increases. Stress distribution of rock and soil is triangular before and after anti-slide pile, pressure distribution of rock and soil is complex under the critical slip surface, and can be considered to be a trapezoid or rectangular distribution. Shear force and bending moment of anti-slide pile are the maximum on the critical slip surface.

reconstruction; anti-slide pile; mechanism of piles of rock; reinforcement effect

10.11817/j.issn.1672-7207.2015.06.036

TU751+.7

A

1672−7207(2015)06−2244−06

2014−08−03；

2014−10−12

国家自然科学基金资助项目(51278067)；广西科学研究与技术开发计划项目(桂科转1298011-2)(Project (51278067) supported by the National Natural Science Foundation of China; Project (1298011-2) supported by Science and Technology Development Plan of Guangxi)

何忠明，博士，副教授，从事道路工程教学与科研工作；E-mail：hezhongming45@126.com