动载下岩石边坡稳定性研究①_邓雄武.pdf
动载下岩石边坡稳定性研究 ① 邓雄武, 余 洋, 李绍波 (云南工程建设总承包公司,云南 昆明 650000) 摘 要 采用拟动力法分析了岩石边坡在动载力作用下的极限承载力。 采用 Hoek⁃Brown 强度准则描述岩体的抗剪特性,并得到等 效的 Mohr⁃Coulomb 强度参数。 基于离散法构建破坏模型,在此基础上,通过极限分析上限法推导并计算得到了极限承载力的上限 解。 分析表明,岩石本身强度参数对其稳定性影响比较明显;而水平动载加速度系数的影响明显大于竖直动载加速度系数的影响。 初始时间差和岩体放大系数对极限承载力有显著影响。 因此,在存在坡顶荷载的边坡稳定性设计中,应考虑因边坡地形等因素导 致的动载作用放大的影响。 关键词 岩石边坡; 极限承载力; 拟动力方法; 极限分析; 离散技术 中图分类号 TU457文献标识码 Adoi10.3969/ j.issn.0253-6099.2018.05.007 文章编号 0253-6099(2018)05-0029-04 Analysis of Rock Slope Stability under Dynamic Loading DENG Xiong⁃wu, YU Yang, LI Shao⁃bo (Yunnan Construction Engineering General Contracting Company, Kunming 650000, Yunnan, China) Abstract A pseudo⁃dynamic approach was adopted to analyze the ultimate bearing capacity of slope under dynamic loading. The Hoek⁃Brown yield criterion was used to describe the shear resistance of rock mass and the equivalent Mohr⁃Coulomb strength parameter was obtained. A failure mode was constructed based on the discretization technique and the upper⁃bound solution of ultimate bearing capacity of rock slope was derived based on upper⁃bound theory of limit analysis. It is shown that the strength parameter of rock has brought an obvious impact on its stability, and horizontal dynamic load acceleration coefficient obviously has brought more impact than vertical dynamic load acceleration coefficient. While the initial time difference and amplification factor of rock mass have brought remarkable impact on the ultimate bearing capacity of slope. It is concluded that the effect by amplification of dynamic load action caused by slope topography shall be taken into consideration in the stability design for the slope with load on its crest. Key words rock slope; ultimate bearing capacity; pseudo⁃dynamic approach; limit analysis; discretization technique 在动载的影响下,边坡更易失稳。 因此在动载条 件下,对边坡进行准确地安全评估是岩土工程师们特 别关注的问题。 在已有的研究中,拟静力方法常被用 于评估动载下的边坡稳定性。 Yang[1]等基于这种方 法得到了动载条件下边坡的屈服加速度。 然而,拟静 力法假设动载加速度不随时间发生变化,也无法考虑 动载加速度的相位影响,这些都与实际情况不符,因此 拟动力方法得以推广。 过去主要通过极限平衡方 法[2-3]来得到挡土墙上的土压力。 Choudhury[4-5]等在 不同工况下对动载作用下墙后填土的主动、被动土压 力进行研究,并分析了临水挡土墙的稳定性。 极限平衡法[6]和极限分析法[7]常被用于边坡稳 定性研究,而对于极限分析法,不同的破坏机理对上限 解有显著的影响。 强度参数的变化,尤其是摩擦角的 变化对边坡破坏机理有很重要的影响。 孙志彬[8]将 离散化技术运用到评价边坡稳定性当中,这种方法能 够有效解决当土体存在空间异性时,简单的对数螺旋 线不能表示其破坏面的形状的问题。 本文提出了一种岩石边坡抗震稳定性的分析方 法,应用拟动力方法引入正弦函数形式的动载加速度 来考虑动载的动态效应。 为了便于考虑沿深度变化的 动载加速度,采用“逐点生成”的离散化技术来生成潜 在的破坏面。 在研究中,采用等效法来考虑非线性 Hoek⁃Brown 屈服准则,得到相应的 Mohr⁃Coulomb 强度 参数。 最后通过非线性优化得到动载影响下边坡极限 承载力的最小上限解。 ①收稿日期 2018-04-04 作者简介 邓雄武(1985-),男,云南昆明人,工程师,主要从事交通土建技术工作。 第 38 卷第 5 期 2018 年 10 月 矿矿 冶冶 工工 程程 MINING AND METALLURGICAL ENGINEERING Vol.38 №5 October 2018 ChaoXing 1 基本原理 1.1 岩体 Hoek⁃brown 屈服准则 在岩土工程中应用最广泛的本构关系是 Mohr⁃ Coulomb 屈服准则,这一准则既适用于土体材料也适 用于岩体材料。 但是,通过三轴试验发现,随着节理和 断裂程度的不同,岩石的破坏具备明显的非线性特征, 一些学者也进行了相应的研究[9-10]。 为了考虑岩石破 坏时主应力间的非线性关系,改进过的 Hoek⁃Brown 屈 服准则[11]可以表示为 σ1 - σ 3 = σ c mσ3 σc + s ■ ■ ■ ■ ■ ■ n (1) 式中 σ1和 σ3分别为失效状态下的大、小主应力;σc 为岩石的单轴抗压强度;m,s 和 n 是由地质强度指标 (GSI)确定的强度参数,其取值方法参照文献[12]。 基于拟合方法,以 Mohr⁃Coulomb 强度准则的表达 式为基础,得到了用 Hoek⁃Brown 的等效粘聚力 c′和内 摩擦角 φ′ c′ = σc(1 + 2n)s + (1 + n)mσ3n[](s + mσ3n)n-1 (1 + n)(1 + 2n)1 + 6nm(s + mσ3n)n-1/(1 + n)(1 + 2n)[] (2) φ′ =arcsin 6nm(s + mσ3n)n-1 2(1 + n)(1 + 2n) + 6nm(s + mσ3n)n-1 ■ ■ ■■ ■ ■ ■■(3) 1.2 拟动力方法 传统的拟静力方法引入水平和竖直方向的动载加 速度 khg 和 kvg 来考虑动载效应,该方法忽略了动载加 速度随 时 间 t 和 所 处 竖 直 方 向 位 置 yi的 影 响。 Choudhury 首次将拟动力方法运用到重力式挡土墙的 稳定性分析当中,这种方法考虑了动载加速度的动态 特性,弥补了拟静力方法的不足。 在动载条件下,P 波和 S 波的波速可以表示为 Vp= 2G(1-υ) / ρ(1-2υ)[] 0.5和 V s=(G/ ρ) 0.5,由表达式 可知,波速与剪切模量 G、泊松比 υ 和密度 ρ 有关。 假 设边坡中剪切模量不随深度变化,岩石的泊松比取为 0.3,则 Vp/ Vs=1.87。 并且考虑到动载加速度幅值随高 度的线性变化,引入岩体放大系数 f,水平和竖直加速 度可以表示为 ah= 1 + yi H (f - 1) ■ ■ ■ ■ ■ ■ ■ ■khgsin 2π t T - yi λs ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ av= 1 + yi H (f - 1) ■ ■ ■ ■ ■ ■ ■ ■kvgsin 2π t T - yi λp + t0 T ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ (4) 式中 kh,kv分别为水平和竖直动载加速度系数;g 为重 力加速度;T 为周期;λs( =TVs),λp( =TVp)分别为 P 波 和 S 波的波长;t0为考虑 P 波和 S 波到达坡底时的时 间差。 2 离散化技术 根据相关联流动法则的要求,破坏面上每一点的 速度方向与破坏面所成的夹角应该为材料的内摩擦 角。 对于一般的岩石材料,可以假定其等效内摩擦角 为一定值,在这种条件下,潜在的破坏面可以用对数螺 旋线来表示。 但是,当岩体材料表现出一定的非均匀 性时,难以用简单的函数形式来表示其破坏机理。 因 此以正向差分为原理的离散化技术可以用来生成未知 形式的破坏面。 本文研究的是通过坡趾的破坏形式,所以以坡趾 B 为原点,水平和竖直方向为坐标轴,建立平面直角坐 标系,坡体相关几何参数如图 1 所示。 为了简化问题, 本文做出了如下假设与规定① 岩石材料服从 H⁃B 屈 服准则,并通过拟合得到等效的 M⁃C 强度参数;② 每 一个无穷小的三角形均视为刚体,并围绕同一个旋转 中心旋转。 H AB y x P1 P2 P3 P4 r1 ri θi θn θ1 Pi Pi1 Pn mqH Ox0, y0 C q D β ω 逐点生成 δθ 图 1 离散法计算示意 对于任意的初始角度 θ1和初始半径 r1(具体几何 意义见图 1),可得到旋转中心 O 的坐标。 破坏面的初始点为坡趾 B,记作点 P1,其坐标为 (0,0),下一个任务是搜索点 P2。 在这之前,需要给定 两条相邻旋转半径 ri和 ri+1之间的旋转角度 δθ,这个 角度将会影响生成破坏面的精度。 下面以通过已知点 Pi搜索点 Pi+1为例,阐述整个破坏面的生成过程。 Pi点为破坏面上任意一个已知点,如图 2 所示, 其坐标为(xi,yi),与旋转中心 O 的距离为 ri,旋转角 度为 θi,根据已有的假设和几何推导,线段 PiPi+1的长 度可用下式表达 PiPi+1 → = risinδθ sin(π/2 + φ′ - δθ) (5) 03矿 冶 工 程第 38 卷 ChaoXing ri θi Pi1 x Ox0, y0 Pixi , yi δθ vi φ α i 图 2 离散微元几何关系 在向量PiPi+1 →的长度和方向均已知的情况下,点 Pi+1的坐标可以计算得到。 当 Pi+1点的纵坐标等于边 坡高度 H 时,整个搜索过程终止,最终的破坏面生成。 3 岩石边坡稳定性 运用上限定理求解极限承载力[13]时,需要得到基 于破坏体的外力功率 W 和内部能量耗散率 D。 在本 文中,外力包括重力、坡顶荷载及相应的惯性力;由于 刚性体假设,内部能量耗散仅发生在速度间断面上。 当坡顶荷载达到其极限承载力时,有 W = D(6) 3.1 外力功率 为了方便计算动载力以及相应的惯性力功率,将 破坏体划分为若干个梯形单元 PiPi+1Qi+1Qi,单元面积 为 Si,重心 Ci的坐标记为(xci,yci),重心 Ci到旋转中 心 O 的距离为 rci,在得到这些参数之后,重力功率可 表示为 WG=∑ i Siγωricos(π - θci)(7) 式中 ω 为破坏体绕 O 旋转的角速度。 同理,可以得到 水平和垂直方向的动载力功率。 对于坡顶荷载,认为只有在破坏体范围内的坡顶 荷载对边坡稳定性是不利的,而破坏体外的坡顶荷载 对边坡破坏没有影响。 在本文中,一个与 C 点距离为 mqH 的均布荷载 q 作用在坡顶平面上。 在计算其功率 时,将其等效为一个集中力,xPn是破坏面与坡顶平面 交点的横坐标。 因此,坡顶荷载引起的功率为 Wq= q xPn- Hcotβ - mqH ω x Pn - xPn+ Hcotβ - mqH 2 - x o ■ ■ ■ ■ ■ ■(8) 在计算坡顶荷载相应的惯性力功率时,动载加速 度中的位置参数 yi等于边坡的高度 H,相应的,可以 得到其功率表达式。 3.2 内部能量耗散率 内部能量耗散可通过破坏面上每一离散单元 PiPi+1的能量耗散求和得到,表达式为 D = ∑ i Di=∑ i c′Liωricosφ′(9) 式中 Li为离散单元 PiPi+1的长度;ri为点 Pi到旋转中 心的距离;φ′为等效摩擦角。 3.3 极限承载力 对于任意一组满足要求的初始半径和初始角度 r1,θ1和不同的时间 t,都能得到一个极限承载力的上 限解答,即极限承载力能够表示成上述三个参数的函 数qcr=f(r1,θ1,t)。 通过编制 Matlab 程序,结合非线 性优化,能够得到极限承载力的最小上限解。 4 稳定性分析与讨论 为便于对相关参数进行讨论分析,假定各参数的 取值如下岩石重度 γ=25 kN/ m3,边坡高度 H=10 m, 坡脚 β= 50,荷载距离系数 mq= 0.1,地震波周期 T= 0.2 s,岩石单轴抗压强度 σc= 2 MPa,S 波波速 Vs= 2 055 m/ s,P 波波速 Vp=3 843 m/ s。 4.1 强度参数对极限承载力的影响 当地质强度指标 GSI 从20 变化到40,极限承载力 随扰动系数 D 和 mi变化趋势分别如图3 和图4 所示。 GSI 600 500 400 300 200 100 0 2025303540 qcr/ kN m-1 ■ ● ▲ ▲ ■ ● ▲ ▲ ■ ● ▲ ▲ ■ ● ▲▲ ■ ● ▲▲ ■ ● ▲▲ D 0.00 D 0.05 D 0.10 D 0.15 mi 0 图 3 极限承载力 qcr随扰动系数 D 变化趋势 GSI 1000 800 600 400 200 0 2025303540 qcr/ kN m-1 ■ ● ▲ ▲ ■ ● ▲ ▲ ■ ● ▲ ▲ ■ ● ▲ ▲ ■ ● ▲ ▲ ■ ● ▲ ▲ mi 9 mi 11 mi 13 mi 15 D 0 图 4 极限承载力 qcr随扰动系数 mi变化趋势 从图 3 可以看出,随着地质强度指标 GSI 增大,极 限承载力增大;而随着扰动系数 D 增大,极限承载力 减小;随着扰动系数 mi增大,极限承载力增大。 因此 13第 5 期邓雄武等 动载下岩石边坡稳定性研究 ChaoXing 可以得出,岩石的固有性质和受扰动程度对动载影响 下边坡的极限承载力有较大影响。 4.2 水平动载加速度系数 kh对极限承载力的影响 为了分析水平动载加速度系数对动载条件下极限 承载力的影响,对相关参数取值如下GSI = 25,mi= 10,D=0, f=1.0,t0= 0,kv =k h。 当水平动载加速度系 数从 0 增加到 0.15 时,极限承载力如表 1 所示,相应 的临界破坏面如图 5 所示。 表 1 不同 kh条件下,qcr计算结果(单位kPa) kv/ kh kh 0.000.050.100.150.20 0.0388.12277.05186.49116.8661.63 0.5388.12271.53180.17111.7358.55 1.0388.12266.23174.28107.0655.04 距离/m 12 10 8 6 4 2 0 2046814121016 高度/m kh qcr 0.00 388.12 kPa kh qcr 0.10 174.28 kPa kh qcr 0.20 55.04 kPa 图 5 kh对破坏体的影响 从表 1 可以看出,随着水平动载加速度系数增大, 极限承载力降低;随着垂直加速度系数与水平加速度 系数比值增大,极限承载力虽然有所降低,但变化幅度 并不大。 因此可以得出,水平动载加速度对边坡的破 坏起主要作用,而垂直动载加速度的影响在工程中基 本可以忽略不计。 从图 5 可以看出,破坏体的大小随 着水平动载加速度系数增大而减小。 4.3 初始时间差 t0和岩体放大系数 f 对极限承载力 的影响 为了分析岩体放大系数 f 对动载条件下极限承载 力的影响,对相关参数取值如下GSI=30,mi=10,D=0, kh=0.1,kv=0.5kh。 当初始时间差从-T/2 变化到 T/2 时,同时考虑不同的岩体放大系数,极限承载力如图 6 所示。 根据图 6 所示,当初始时间差 t0从-T/2 变化到 0 时,极限承载力减小,而 t0从 0 变化到 T/2 时,极限 承载力增大,并且在一个周期之后,极限承载力又回到 了初始值,表现出一定的周期性;而随着岩体放大系数 增大,动载条件下边坡的极限承载力显著降低。 t0/T 350 300 250 200 150 -0.50-0.250.000.250.50 qcr/kPa ■▲ ▲ ■ ● ▲ ▲ f 1.00 f 1.75 ● f 1.25 f 2.00 ◆ ■ ● ▲ ▲ ◆ ■ ● ▲ ▲ ◆ ■ ● ▲ ▲ ◆ ■ ● ▲ ▲ ◆ ■ ● ▲ ▲ ◆ ■ ● ▲ ▲ ◆ ◆ f 1.50 图 6 t0/ T 对 qcr的影响 5 结 论 采用拟动力法分析了岩石边坡在动载力作用下的 极限承载力。 采用 Hoek⁃Brown 强度准则描述岩体的 抗剪特性,并得到等效的 Mohr⁃Coulomb 强度参数。 基 于离散法构建破坏模型,在此基础上,通过极限分析上 限法推导并计算得到了极限承载力的上限解。 研究结 果表明,本文方法具有良好的适用性。 参数分析表明, 岩石本身强度参数对其稳定性影响较明显,水平动载 加速度系数 kh的影响明显大于竖直动载加速度系数 kv的影响。 初始时间差 t0和岩体放大系数 f 对极限承 载力有显著影响。 因此,在存在坡顶荷载的边坡稳定 性设计中,应考虑因边坡地形等因素导致的动载作用 放大的影响。 参考文献 [1] Yang X L. Seismic displacement of rock slopes with nonlinear Hoek⁃ Brown failure criterion[J]. International Journal of Rock Mechanics and Mining Sciences, 2007,44(6)948-953. [2] 王志凯,夏唐代,陈炜昀. 刚性挡土墙地震主动土压力的拟动力学 分析[J]. 浙江大学学报(工学版), 2012,46(1)46-51. [3] 黄 睿,夏唐代,钟丽娜,等. 地震作用下挡土墙抗滑动稳定性的 拟动力分析[J]. 中南大学学报(自然科学版), 2014(9)3186- 3192. [4] Choudhury D, Nimbalkar S S. Pseudo⁃dynamic approach of seismic active earth pressure behind retaining wall[J]. Geotechnical and Geo⁃ logical Engineering, 2006, 24(5)1103-1113. [5] Choudhury D, Nimbalkar S. Seismic passive resistance by pseudo⁃dy⁃ namic method[J]. Gotechnique, 2005,55(9)699-702. [6] 罗根传,胡庆国,谭积青. 基于极限平衡法与强度折减法的高边坡 过程稳定性分析[J]. 矿冶工程, 2013,33(2)14-17. [7] 刘秀军,贺建清. 径向条分极限分析法在土坡稳定性分析中的应 用[J]. 矿冶工程, 2011,31(2)12-15. [8] 孙志彬,潘秋景,杨小礼,等. 非均质边坡上限分析的离散机构及 应用[J]. 岩石力学与工程学报, 2017,36(7)1680-1688. [9] 杨小礼. 岩石极限分析非线性理论及其应用[J]. 中南大学学报 (自然科学版), 2009,40(1)225-229. (下转第 35 页) 23矿 冶 工 程第 38 卷 ChaoXing 从表 3 可看出,使用相同原材料,工艺不同,吸水 率差异大。 常压灌注浇筑体耐静水压为 50 MPa,真空 浇注的浇筑体耐静水压为 60 MPa。 真空浇注浇筑体 的吸水率也明显低于常压灌注浇筑体。 真空浇注更能 提高胶液与粉体的浸润,减少微小气泡,提高体系堆积 系数。 故此采用真空浇注进行实验。 3.4 综合性能测试 通过以上的试验条件对比,最终选用低粘度胶液, 大小粒径混合空心微珠,采用真空浇注的方法进行浇 注,制成轻质浮体材料,测试其不同静水压力下的吸水 率,结果见表 4。 表 4 浇注体综合性能测试 初始质量 / g 静水压力 / MPa 保压时间 / h 静水压下质量 / g 吸水率 / % 4024138.240.04 4048138.270.06 138.184096138.290.07 5024138.440.18 6024138.680.36 4024259.160.32 4048259.410.42 258.324096259.540.47 5024259.600.49 6024259.630.50 该材料密度为 0.58 g/ cm3,分别通过了 40 MPa, 50 MPa 和 60 MPa 的静水压力测试,其吸水率均低于 0.5%。 在此基础上进行了其他性能测试,测得压缩强 度 55 MPa,压缩模量 1 800 MPa,拉伸强度 25 MPa,体 积模量 2 500 MPa,热变形温度 110 ℃,耐海水 30 d 增 重小于 0.5%,静浮力大于 50%。 其综合性能已达到试 验要求。 4 结 语 1) 低粘度改性环氧树脂胶液与空心玻璃微珠有 很好的浸润性能,固化后力学性能优良,为浇筑体的制 备提供了基础的保证。 2) 通过对不同粒径空心玻璃微珠的选择,确立了 最佳粉料;对浇注工艺的比较,确立其成型工艺。 研制 材料达到试验目标。 3) 研究结果表明,该生产工艺可行,可根据不同 应用领域的要求进行不同形状规格、异形件的浇筑。 该固体浮体材料能为深海 4 000 m 水下作业、深海探 矿及科考提供浮力。 材料的使用可以提高潜器有效载 荷,减少其外形尺寸,提高其水下运动性能。 它将广泛 应用于海洋探矿、水下机器人、水下管线、潜器、浮标、 救生装置等。 该浮力材料的产业化将为海洋事业发展 提供性能稳定可靠的材料保障,为深海开发技术的不 断发展和全球深海战略的推进做出重要贡献。 参考文献 [1] 王德中. 环氧树脂生产与应用(第 2 版)[M]. 北京化学工业出 版社, 2001. [2] 陈洪江,吴 敏,虞鑫海,等. 新型环氧树脂胶粘剂的固化动力学 研究[J]. 粘接, 2009,30(8)43-45. [3] 刘淑青. 高强度轻质浮力材料研究[J]. 海洋技术学报, 2007,26 (4)118-120. [4] 黄海塘,罗丹凤, 黄桂英. 海洋资料浮标体材料研究与应用进 展[J]. 数字化用户, 2017,27(16)18-19. [5] 高 昂,胡明皓,王勇智,等. 深海高强浮力材料的研究现状[J]. 材料导报, 2016,30(2)80-83. [6] NAVY MIL-S-24154A-1967, Syntactic buoyancy material for high hydrostatic pressures[S]. 引用本文 熊 利,许晓武,金 星. 深海固体浮力材料的研制及性能 探讨[J]. 矿冶工程, 2018,38(5)33-35. (上接第 32 页) [10] 杨小礼,王作伟. 非线性破坏准则下浅埋隧道围岩压力的极限 分析[J]. 中南大学学报(自然科学版), 2010,41(1)299- 302. [11] Hoek E, Brown E T. Practical estimates of rock mass strength[J]. International Journal of Rock Mechanics and Mining Sciences, 1997, 34(8)1165-1186. [12] Hoek E, Carranza⁃Torres C. Hoek⁃Brown failure criterion2002 Edition[C]∥Proceedings of the Fifth North American Rock Me⁃ chanics Symposium, 200218-22. [13] 尹 鑫,周海祚,郑 刚. 地震作用下临近边坡的条形基础极限 承载力研究[J]. 岩土工程学报, 2017,39(s2)95-98. 引用本文 邓雄武,余 洋,李绍波. 动载下岩石边坡稳定性研究[J]. 矿冶工程, 2018,38(5)29-32. 53第 5 期熊 利等 深海固体浮力材料的研制及性能探讨 ChaoXing