徐爱民 ，柳群义, ，朱自强，鲁光银

(1. 湖南省高速公路管理局 湖南 长沙，410003；

2. 中南大学 地球科学与信息物理学院 湖南 长沙，410083；3. 中国地质科学院 矿产资源研究所，北京，100037)

抗滑桩在边坡加固工程中广泛使用，许多研究者对此进行了大量研究[1-4]，但他们主要采用极限平衡法研究边坡稳定性以及抗滑桩受力情况，该方法需对桩土受力以及边坡滑动面进行假设，无法清晰反映边坡-抗滑桩的耦合效应以及边坡的滑移特征，并且由于桩和岩土介质具有不同的弹性模量、重度和强度，势必引起应力分布的非均匀化和进入塑性极限的不同步，给边坡的稳定分析带来了较大的困难[5]。近年来，随着计算机技术的不断发展，数值方法在岩土工程分析中逐渐得到广泛应用[6-9]，该法不受岩土体几何形状不规则和材料不均匀性的限制，如Jeong等[10]利用极限平衡法和有限差分软件FLAC3D 分析了边坡-抗滑桩系统的稳定性。但是，这些研究对象主要是针对土质边坡或者不含结构面的岩体边坡，对于层状岩体中抗滑桩加固效果的研究较少，在此，本文作者利用FLAC3D数值计算方法，探讨层状岩体边坡抗滑桩的加固效应，分析抗滑桩位置和桩长对于边坡变形和应力的影响。

1 桩-岩相互作用特征

桩单元包括3个弹簧单元，其中2个弹簧单元平行于桩轴线以提供轴力和轴向剪力，另一个垂直于轴线以提供横向抗力，如图1所示。桩单元被分为许多小段，通过这些小段的积分得到桩整体的变形和应力状态[11]。……