基于改进最小势能法的水平–竖向立体筋边坡稳定性分析研究

蒋建国; 陈媛; 刘孙光; 赵发嘉

doi:10.11779/CJGE202105008

基于改进最小势能法的水平–竖向立体筋边坡稳定性分析研究 English Version

中南大学土木工程学院,湖南　长沙　410075

详细信息

作者简介:
蒋建国(1973—),男,湖南洞口人,副教授,主要从事路基路面工程、岩土工程方面的研究工作。E-mail: yuanchencsu@126.com

中图分类号: TU473
计量
- 文章访问数: 257
- HTML全文浏览量: 23
- PDF下载量: 132
出版历程
- 收稿日期: 2020-05-10
- 网络出版日期: 2022-12-04
- 刊出日期: 2021-04-30

Stability analysis of horizontal-vertical ribs-reinforced slopes based on improved minimum potential energy method

School of Civil Engineering, Central South University, Changsha 410075, China

摘要

摘要: 基于改进的最小势能算法提出了一种应用于水平–竖向（horizontal-vertical,简称H-V）立体筋边坡稳定性分析的新方法。在建立边坡支护系统的势能函数时考虑水平筋的摩擦势能、线弹性应变能以及竖筋变形势能的影响。通过单元体虚位移方向的静力平衡方程分析加筋边坡剪切势能,构建了适用于施加水平筋和H-V立体筋边坡的势能计算模型。通过计算系统势能最小时的虚位移进而求解边坡安全系数,简化了加筋边坡的分析过程,提高了计算效率,经与已有文献计算结果进行对比分析,证明计算方法的可行性。分析了不同土体强度参数黏聚力和内摩擦角对加筋边坡稳定性的影响,为立体加筋技术在工程中应用和推广提供参考。
- 改进最小势能法 /
- 边坡稳定性 /
- 水平–竖向体筋（H-V筋） /
- 安全系数
Abstract: Based on the improved minimum potential energy algorithm, a new method for the stability analysis of horizontal-vertical (H-V) ribs-reinforced slopes is proposed. In this method, the frictional potential energy of horizontal reinforcement, linear elastic strain energy and deformation potential energy of vertical reinforcement are considered when the potential energy function for the slope-supporting system is established. A potential energy model suitable for the studied slope is established through the shear potential energy of the reinforced slope according to the static equilibrium equation for virtual displacement. The analysis process of the reinforced slope is simplified, and the calculation efficiency is improved by solving the safety factor of the slope with the virtual displacement of the minimum potential energy of the studied slope. By comparing and analyzing the calculated results with those of the published papers, the feasibility of the proposed method is proved. The effect of strength parameters of soils such as cohesive force and internal friction angle on the slope stability is also analyzed, which provides reference for the application and promotion of three-dimensional reinforced slope technology in engineering.
- improved minimum potential energy method /
- slope stability /
- horizontal-vertical ribs /
- safety factor

HTML全文

0. 引言

自1960年Henri Vidal^[1]首次提出了“加筋土”的概念和设计理论以来,土工合成材料广泛应用于岩土工程领域^[1]。在边坡中铺设土工合成材料能够优化土体应力分布,提高边坡的稳定性和承载力^[2-4]。有别于二维土工合成材料如土工格栅、土工布。张孟喜^[5]首先提出H-V立体筋的概念。作为一种新型加筋方式,H-V立体筋在筋–土相互作用上展现出更好的优势,除了具有水平筋作用外,竖筋能约束土体的变形并在竖筋间形成“加强土楔体”或“加强土柱”,增强土体的抗拉强度和抗剪强度^[6],具有较好的应用前景。

常用的加筋边坡稳定性分析方法主要有极限平衡法、有限元法等。改进最小势能法作为边坡稳定性分析方法之一,研究成果表明,相对极限平衡法,该方法计算方便快捷,具有不需要条块划分、不需要反复迭代求解等优点^[7-9]。目前,改进最小势能法主要应用于有加固措施如锚杆、锚索、抗滑桩的边坡稳定性分析中。邹治来^[10]考虑土工格栅与土体接触面之间摩擦势能和土工格栅弹性拉伸势能,建立土工格栅加筋边坡的最小势能稳定性分析方法。但是分析过程中,由于土工格栅与土体之间摩擦阻力计算过程复杂,参数众多,需要进行拉拔试验才能确定参数值,实现有很大难度。采用最小势能法分析了普通土工格栅加筋边坡稳定性,李铀等^[11]、孙加平等^[12]提出了考虑滑床剪切势能影响的近似计算模型,改进最小势能对锚杆（索）加固边坡稳定分析方法。温树杰等^[13]分析桩体变形特征,考虑抗滑桩与滑体相互接触处的弹性压缩势能、挠曲势能提出基于最小势能法加固边坡稳定性分析方法。未检索到最小势能法应用于H-V立体筋加筋边坡稳定性分析。

H-V立体筋边坡稳定性分析主要采用条分法,张孟喜等^[14]采用瑞典条分法对立体加筋边坡的稳定性进行探讨,在计算过程中未考虑条间力的作用,容易导致计算结果偏小。马学宁等^[15]建立水平条间力的假设条件,推导出水平条分瑞典法和简布法分析H-V立体加筋边坡的稳定性,该方法需要在水平筋材处于土条的中线上的条件下进行,且水平条分简布法需要进行反复迭代,计算量较大。其它对于H-V立体筋边坡稳定性分析文献尚不多见。

因此,本文基于改进最小势能法,对H-V立体筋作用的基本原理和受力特性进行分析,构建出水平筋和竖筋的弹性势能方程,并对加筋边坡单元体虚位移方向的剪切势能进行推导,通过计算边坡最危险滑裂面系统势能最小时的虚位移,求解出边坡安全系数,并通过算例对该方法的适用性和可行性进行验证。

本文利用改进最小势能法从能量的角度对H-V立体筋边坡的稳定性展开分析,在系统中涉及的能量种类主要有滑体储存的弹性势能,H-V筋因水平筋线弹性变形产生弹性势能、筋土之间因摩擦储存摩擦势能、竖筋产生侧向变形储存弹性势能和滑面剪切势能,概念清晰,可直接计算加筋边坡稳定性分析结果,相对于计算复杂的有限元法和需反复迭代、需对超静定问题进行多次假定的极限平衡法更有实用价值,为H-V立体筋边坡稳定性分析提供新的思路。

1. 滑体储存的弹性势能

如图1所示的H-V立体筋均质边坡,设潜在滑裂面为AC,其曲线方程表示为 $y = f (x)$ ,土体黏聚力为 $c$ ,内摩擦角为 $φ$ ,重度为 $γ$ ,土体的重力为 $W$ ,第 $i$ 排筋条与水平方向的夹角为 $α_{i}$ ,水平方向布置则 $α_{i} = 0$ ；第 $i$ 排筋条与竖直方向的夹角为 $β_{i}$ ,每排筋条的单位方向向量记为 $\vec{T} = (\cos α_{i}, \cos β_{i})$ 。为合理建立力学模型,做如下说明：①假设滑动体整体为刚性体,在滑裂面上存在弹性压缩变形和剪切变形,储存的弹性压缩势能为 $V_{e}$ ,剪切势能为 $V_{τ}$ ；②在合外力 $\vec{R}$ 的作用下加筋边坡与滑体产生相同的虚位移 $\vec{d} =$ $(d_{1}, d_{2})$ ；③微分段的长度为 $d l$ ,弹性压缩变形可用刚度系数为 $K_{i}$ 的弹簧模拟,其大小与微段 $d l$ 呈正比,即 $K_{i} = m_{i} d l$ , $m_{i}$ 为地基系数,可通过试验测定或查阅相关资料获得。

图 1 H-V加筋均质边坡计算图示

Figure 1. Calculation of H-V ribs-reinforced homogeneous slope

下载: 全尺寸图片幻灯片

滑裂面方程为 $y = f (x)$ ,则求得曲线上任意切线上法线的斜率为 $k = - 1 / f^{'} (x)$ 。滑裂面外法线单位向量和滑面切向单位向量可表示为

$\vec{n} (x) = [n_{1} (x), n_{2} (x)] = {\frac{f^{'} (x)}{\sqrt{1 + {[f^{'} (x)]}^{2}}}, \frac{- 1}{\sqrt{1 + {[f^{'} (x)]}^{2}}}}$ ,

(1)

$\vec{t} (x) = [t_{1} (x), t_{2} (x)] = {\frac{1}{\sqrt{1 + {[f^{'} (x)]}^{2}}}, \frac{f^{'} (x)}{\sqrt{1 + {[f^{'} (x)]}^{2}}}}$ 。

(2)

假设在各种力的作用下,滑体发生虚位移,使系统势能最小,滑体此时的势能即是系统达到极限平衡状态时的最小势能,则任意微段储存的弹性势能为

$V_{e1} = \sum_{i = 1}^{N} \frac{1}{2} K_{i} {(\vec{d} \cdot \vec{n})}^{2} = \frac{1}{2} \int_{l} m {(\vec{d} \cdot \vec{n})}^{2} d l$ 。

(3)

2. H-V立体筋势能构建

H-V立体筋结构在传统水平筋的基础上设置竖筋,竖筋宽度可与水平筋相同,其截面形状可根据工程需要采用不同厚度变化的矩形、多边形等,形成一种空间形式的三维组合加筋结构,如图2所示。

图 2 H-V立体筋三维结构示意图

Figure 2. Three-dimensional structural diagram of H-V ribs

下载: 全尺寸图片幻灯片

H-V立体筋结构的基本原理：将一定量具有抗拉强度的H-V立体筋铺设在土体中,水平筋的拉力作用会逐渐延伸到路堤边坡自由面的筋条中,水平筋土之间摩擦力提供锚固力；竖筋的作用则表现为增加筋材拔出过程中竖筋的抗阻力,使土体结合更紧密,将整个H-V筋锚固在了土体中,提高边坡承载力。相对普通均质边坡,H-V筋的这种拉力状态也会使上部土体传递来的荷载较为集中,相对更均匀地分布到筋材下部的土体,使下部土体受力更合理。

2.1 筋条特征坐标的确定

构建H-V立体筋势能函数时需要确定水平筋与滑面交点坐标和竖筋有效数目,因此首先需要确定筋条的特征坐标。加筋边坡在发生失稳破坏时,根据坡体的滑动状态可分为滑动区和锚固区,如图3所示,设坡面倾角为 $θ$ ,筋材各层的垂直间距为 $h_{i}$ ,第一层水平筋距离坡脚的竖直距离为 $h_{1}$ 。则第 $i$ 层水平筋与滑动面的交点纵坐标 $y_{i}$ 可表示为

图 3 筋条的特征坐标计算简图

Figure 3. Simplified calculation of characteristic coordinates of geogrid

下载: 全尺寸图片幻灯片

$y_{i} = y_{1} + h_{i} (i - 1)$ 。

(4)

将式（4）与滑动面方程 $y = f (x)$ 联立即可确定各排筋条与滑面的交点坐标,若方程组有解,则筋条产生变形储存了势能,确认筋条储存势能且取 $x > 0$ 时,第 $i$ 层水平筋与滑动面的交点横坐标 $x_{i}$ 可表示为

$x_{i} = f^{'} (y_{i}) (x > 0)$ 。

(5)

土工格栅长度为 $l$ ,则第 $i$ 层水平筋右末端横坐标 ${x^{'}}_{i}$ 表示为

${x^{'}}_{i} = x_{1} + l + h_{i} (i - 1) / \tan θ$ 。

(6)

由式（5）,（6）可得第 $i$ 层土工格栅的锚固区长度为

$l_{i} = {x^{'}}_{i} - x_{i}$ 。

(7)

设竖筋间距为D,则锚固区内第 $i$ 层内竖筋的个数 $n_{i}$ 可表示为

$n_{i} = l_{i} / D + 1$ 。

(8)

则可计算出锚固区内竖筋的总个数 $n = \sum_{i = 1}^{m} n_{i}$ 。

2.2 水平筋弹性势能

在土体滑动的过程中,水平筋与土体之间摩擦力作功会产生摩擦势能,水平筋自身也会由于拉伸变形产生线弹性应变能。

对水平筋进行作用力分析时,由于筋土之间接触复杂,水平筋和填土分开考虑计算筋材的摩擦力参数较多,实现起来有很大的难度。加筋路堤稳定性计算时,土工合成材料应用手册将水平筋与路堤填土分开考虑,所以在计算稳定性系数时水平筋的作用力一般直接取其抗拉强度值。本文将水平筋与其附近土体考虑为一种复合体,所以筋条的抗拉力应该分散到其作用范围以内的土体中^[14]。假设筋条的抗拉力平均分配到每层筋条的土体中,若路堤中每层筋条水平间距为 $S_{x}$ ,竖向间距为 $S_{y}$ （图4）。

图 4 水平筋复合抗拉力计算示意图

Figure 4. Schematic diagram of calculation of horizontal ribs

下载: 全尺寸图片幻灯片

则第 $i$ 个土条中单位厚度复合土体中抗拉力为

$T_{st} = \frac{T_{h}}{S_{x} \cdot S_{y}} l_{i} \sin θ_{i}$ ,

(9)

式中, $T_{h}$ 为单根水平筋条抗拉力, $S_{x}$ 为筋条水平间距, $S_{y}$ 为筋条竖向间距, $l_{i}$ 为第 $i$ 条底弧长度, $θ_{i}$ 为第 $i$ 条底弧仰角。

如图5所示,对单位长度为 $l_{i}$ 的土体单元进行分析,水平筋的摩擦阻力 $T_{k i}$ 与水平筋的作用力 $T_{st}$ 平衡,所以满足 $T_{k i} = T_{st}$ ,由于土工格栅所受摩擦阻力与位移呈钝角关系,所以摩擦所储存的势能为力所做功的负值。

图 5 土体单元水平筋受力分析图

Figure 5. Analysis diagram of horizontal reinforcement of soil element

下载: 全尺寸图片幻灯片

$V_{e2} = - \sum_{i = 1}^{n} T_{k i} (\vec{d} \cdot \vec{T})$ ,

(10)

式中, $T_{k i}$ 为第 $i$ 排水平筋摩擦阻力, $\vec{d}$ 为滑体虚位移, $\vec{T}$ 为水平筋水平轴向单位向量。

根据以往文献资料结果显示,基于筋-土界面加筋摩擦机理的圆弧滑动分析方法可以计算加筋后的安全系数,在抗滑项中加入界面摩擦阻力,但结果可能偏小^[16]。本文将水平筋与其附近土体考虑为一种复合体,因为筋土间的界面被忽略,在计算水平筋所受摩擦阻力时缺少界面摩擦阻力,抗滑力中缺少界面摩擦项,且土工格栅摩擦阻力是在抗拉强度条件下计算的,因此可能导致安全系数计算结果偏大,具体影响还需要进一步研究。

分析水平筋的线应变势能时,滑体在向最小势能方向移动时会产生虚位移,因筋材均必须要伸入到较深坡体内部,因此可假设筋材与右端与滑裂面的交点g处是固定的^[11]。设布置了 $m$ 排筋材,筋条的刚度为 $K^{'}$ ,如图6所示,在合外力作用下滑体发生虚位移 $\vec{d}$ ,水平筋的坡面端头由 $e$ 移至 $e^{'}$ ,伸长量为 $Δ_{k^{'}}$ ,假设 $e$ 坐标为（ $x_{e}$ , $y_{e}$ ）, $g$ 的坐标为（ $x_{g}$ , $y_{g}$ ）。可得

图 6 立体筋材结构受力示意图

Figure 6. Schematic stress diagram of three-dimensional rib structure

下载: 全尺寸图片幻灯片

$\begin{matrix} Δ_{k^{'}} = | g e^{'} - g e | \\ = \sqrt{{(d_{1} + v_{k^{'} x})}^{2} + {(d_{2} + v_{k^{'} y})}^{2}} - \sqrt{v_{k^{'} x}^{2} + v_{k^{'} y}^{2}} \end{matrix}$ ,

(11)

式中, $v_{k^{'} x} = x_{g} - x_{e}, v_{k^{'} y} = y_{g} - y_{e}$ 。

水平筋的线应变能为

$\begin{matrix} V_{e3} = \sum_{i = 1}^{m} [\frac{1}{2} K^{'} Δ_{k^{'}}^{2}] \\ = \frac{K^{'}}{2} \sum_{i = 1}^{m} [{(d_{1} + v_{k^{'} x})}^{2} + {(d_{2} + v_{k^{'} y})}^{2} + (v_{k^{'} x}^{2} + v_{k^{'} y}^{2})] - \\ K^{'} \sum_{i = 1}^{m} [\sqrt{{(d_{1} + v_{k^{'} x})}^{2} + {(d_{2} + v_{k^{'} y})}^{2}} \cdot \sqrt{v_{k^{'} x}^{2} + v_{k^{'} y}^{2}}] \end{matrix}$ 。

(12)

2.3 竖筋弹性势能

当H-V立体筋在荷载作用下竖筋首先发挥作用,在其自身刚度下对土体产生侧阻力,阻挡土体向两侧运动,这相当于给土体增加多个侧向约束,竖筋对土体的势能主要是侧阻力对竖筋作用产生的变形势能。

对竖筋的作用力进行分析,竖筋的侧向抵抗作用在竖筋间将形成“土体加强区”,以有效改变加筋土的受力状况^[6]。如图7所示,在竖筋两侧分别作用有主动土压力 $σ_{a}$ 和被动土压力 $σ_{p}$ ,沿竖筋表面作用有剪切应力 $τ_{v}$ , $τ_{h}$ ,在竖筋顶部作用有竖向压力 $σ_{b}$ 。由于竖筋表面积与厚度及其顶部面积相对路堤尺寸可忽略,因此竖筋的主要作用力为竖向筋条两侧的主动土压力及被动土压力,在稳定性分析时,忽略竖筋表面的剪应力作用,竖筋的侧阻力可由被动土压力与主动土压力之差计算得到^[14]。

图 7 竖筋单元受力图

Figure 7. Force diagram of three-dimensional tendons

下载: 全尺寸图片幻灯片

则主动、被动土压力计算公式为

$T_{a} = (γ H_{i} K_{a} - 2 c \sqrt{K_{a}}) A$ ,

(13)

$T_{P} = (γ H_{i} K_{P} + 2 c \sqrt{K_{P}}) A$ 。

(14)

竖筋作用力为 $T_{sh}$ ,

$\begin{matrix} T_{s h} = T_{p} - T_{a} = [T_{p} - T_{a}] \\ = [γ H_{i} (K_{p} - K_{a}) + 2 c (\sqrt{K_{p}} + \sqrt{K_{a}})] A \end{matrix}$ ,

(15)

式中, $c$ 为土体黏聚力, $γ$ 为土体重度, $H_{i}$ 为竖筋上覆土高度, $φ$ 为内摩擦角, $K_{a}$ 为主动土压力系数, $K_{p}$ 为被动土压力系数,且 $K_{a} = \tan^{2} (45^{°} - φ / 2)$ , $K_{a} =$ $\tan^{2} (45^{°} + φ / 2)$ , $A$ 为竖筋侧面积, $A = B h$ ,其中B为竖筋宽度, $h$ 为竖筋高度。

由于竖筋的侧向变形对H-V筋中竖筋作用发挥起重要作用,所以尽管竖筋高度相对边坡尺寸较小,也需要进行势能分析。在进行H-V加筋边坡稳定性分析时,考虑到竖筋的排列会对微分圆弧作用面的受力产生影响,所以在计算抗滑力时只计算锚固区（即破裂面外部）的有效竖筋阻力,锚固区内竖筋总个数为 $n$ ,则竖筋在侧阻力作用下产生的变形势能为

$V_{e4} = \sum_{i = 1}^{n} T_{sh} Δ_{x} = \sum_{i = 1}^{n} T_{sh} d_{1}$ ,

(16)

式中, $T_{sh}$ 为竖筋作用力, $Δ_{x}$ 为竖筋水平方向的变形量,大小为虚位移 $\vec{d}$ 在水平方向的分量。

通过式（15）可看出竖筋作用力 $T_{sh}$ 的作用点位于竖筋高度的中部,但由于竖筋高度 $h$ 较小,在计算立体筋合力时,竖筋作用力的方向可近似认为作用在水平筋上。则立体筋的合力为竖筋作用力 $T_{sh}$ 和水平筋作用力 $T_{st}$ 的和,即

$T_{i} = T_{s t} + n T_{sh}$ 。

(17)

3. 滑面剪切势能的构建

边坡的失稳往往是沿着土体抵抗力最小的方向滑动,即为使坡体系统势能最小的虚位移方向。由文献[8]可知采用弹性势能以及合外力做功的方式表示边坡的剪切势能,可大大简化计算量,同时便于程序的实现。为便于势能公式的表达,令 $Δ = \sqrt{R_{1}^{2} + R_{2}^{2}}$ , $Δ^{'} = \sqrt{d_{1}^{2} + d_{2}^{2}}$ ,则合外力 $\vec{R}$ 和虚位移 $\vec{d}$ 的方向向量可表示为

$\begin{array}{l} {\vec{R}}_{r} = \frac{1}{Δ} (R_{1}, R_{2}), \\ {\vec{d}}_{d} = \frac{1}{Δ^{'}} (d_{1}, d_{2}) 。 \end{array}}$

(18)

从图1中任取土条作为研究对象,土条的受力分析如图8所示。

图 8 边坡土条受力图

Figure 8. Force diagram of slope strips

下载: 全尺寸图片幻灯片

作用在土条上的力有法向力 $N_{i}$ ,剪应力 $τ_{i}$ ,合外力为 $R_{i}$ ,H-V立体筋合力为 $T_{i}$ ,土条左右面上的条间力为 $E_{i}$ , $E_{i + 1}$ 。由于边坡失稳沿着岩土体抵抗力最小的方向,即虚位移 $\vec{d}$ ,因此各个土条左右面上的条间力作用方向与虚位移方向一致,在计算过程中,将第 $i$ 个土条条间力合力记为 $E_{i}$ , $E_{i + 1}$ ,将第 $i + 1$ 个土条条间力合力记为 ${E^{'}}_{i + 1}$ , $E_{i + 2}$ , $Δ E = E_{i + 1} - E_{i}$ ,则各个土条沿虚位移方向满足平衡方程：

$N_{i} \cos α + τ_{i} l_{i} \cos β + T_{i} \cos ω + Δ E_{i} = R_{i} \cos ϕ$ 。

(19)

式中 $α$ 为任意微面上法向力与虚位移方向的夹角,且 $\cos α = \vec{n} \cdot (- {\vec{d}}_{d})$ ； $l_{i}$ 为滑裂面微段长度； $β$ 为剪应力与虚位移的夹角,且 $\cos β = \vec{t} \cdot (- {\vec{d}}_{d})$ ； $ω$ 为水平筋与虚位移的夹角, $\cos ω = \vec{T} \cdot (- {\vec{d}}_{d})$ ； $ϕ$ 为合外力与虚位移的夹角, $\cos ϕ = {\vec{R}}_{r} \cdot {\vec{d}}_{d}$ 。

由式（19）可知滑面上剪应力的解析解为

$τ_{i} = \frac{\frac{R_{i}}{Δ} {\vec{R}}_{i} \cdot {\vec{d}}_{d} + T_{i} \vec{T} \cdot {\vec{d}}_{d} - N_{i} {\vec{n}}_{i} \cdot {\vec{d}}_{d} - Δ E_{i}}{- {\vec{d}}_{d} \cdot \vec{t} l_{i}}$ 。

(20)

为得到线弹性范围内剪切应变势能产生的函数,需要先构建剪应变和剪应力的力学求解模型,土体由于剪切变形发生的剪应变为δ_i,在剪切力 $τ_{i}$ 的作用下滑体沿着滑裂面切向会产生相对位移,即由滑裂面上的点 $f$ 移动至点 $f^{'}$ ,受剪切变形影响的土体深度为 $h_{i}$ ,剪切位移记为 $d^{'}$ ,则

$δ_{i} = \frac{d^{'}}{h_{i}}$ 。

(21)

滑面上的剪切势能为

$\begin{array}{l} V_{τ} = \frac{1}{2} \sum_{i = 1}^{n} τ_{i} δ_{i} h_{i} \\ = - \frac{1}{2} \sum_{i = 1}^{n} \frac{\frac{R_{i}}{△} {\vec{R}}_{i} \cdot {\vec{d}}_{d} + T_{i} \vec{T} \cdot {\vec{d}}_{d} - N_{i} {\vec{η}}_{i} \cdot {\vec{d}}_{d} - Δ E_{i}}{- {\vec{d}}_{d} \cdot \vec{t} l_{i}} \frac{\vec{d} \cdot \vec{t}}{h_{i}} h_{i} l_{i} \\ = - \frac{1}{2} \sum_{i = 1}^{n} (\frac{{\vec{R}}_{i}}{△} \vec{R} \cdot \vec{d} + T_{i} \vec{T} \cdot \vec{d} - N_{i} {\vec{η}}_{i} \cdot \vec{d} - Δ E_{i}) \\ = V_{e 1} - \frac{1}{2} \vec{R} \cdot \vec{d} - \frac{1}{2} \sum_{i = 1}^{n} T_{i} \vec{T} \cdot \vec{d} + \frac{1}{2} \sum_{i = 1}^{n} Δ E_{i} \end{array}$ ,

(22)

式中, $\sum_{i = 1}^{n} Δ E_{i}$ 的物理意义为各个条块条间力的矢量和,由于各个条块之间的作用力属于作用力与反作用力,因此 $\sum_{i = 1}^{n} Δ E_{i} = 0$ ,则上式可化简为

$V_{τ} = V_{e1} - \frac{1}{2} \vec{R} \cdot \vec{d} - \frac{1}{2} \sum_{i = 1}^{n} T_{i} \vec{T} \cdot \vec{d}$ 。

(23)

综合式（3）、（10）、（12）、（17）、（23）可得H-V立体筋边坡系统的总势能为

$\begin{matrix} V = V_{e1} + V_{e2} + V_{e3} + V_{e4} + V_{τ} - \vec{R} \cdot \vec{d} \\ = 2 V_{e1} + V_{e 2} + V_{e 3} + V_{e4} - \frac{3}{2} \vec{R} \cdot \vec{d} - \frac{1}{2} \sum_{i = 1}^{n} T_{i} \vec{T} \cdot \vec{d} \end{matrix}$ 。

(24)

4. 安全系数的计算

安全系数 $F_{s}$ 为沿着虚位移 $\vec{d}$ 方向投影得到的抗滑力 $F_{m}$ 与下滑力 $F_{n}$ 的比值^[8],即

$F_{s} = \frac{F_{m}}{F_{n}},$

(25)

式中,抗滑力主要由微分圆弧 $d l$ 的极限抗滑力和H-V立体筋复合抗拉力产生的抗滑力提供,下滑力由合外力（包括滑体自重）提供。本文的抗滑力 $F_{m}$ 及下滑力 $F_{n}$ 均与虚位移相关,因此需先求解虚位移,再进行安全系数求解。

4.1 虚位移 $\vec{d}$ 的求解

根据最小势能原理可知,H-V立体筋边坡系统势能最小时,总势能函数在虚位移的一阶偏导处需满足

$\begin{array}{l} \frac{\partial V}{\partial d_{1}} = 0, \\ \frac{\partial V}{\partial d_{2}} = 0 。 \end{array}}$

(26)

将式（24）代入式（26）,即可求出系统最小势能时坡体的虚位移 $\vec{d}$ 。其中令

$a_{1} = \sum_{i = 1}^{n} k_{i} η_{x}^{2} = m \int_{x 1}^{x 2} \frac{{[f^{'} (x)]}^{2}}{\sqrt{1 + {[f^{'} (x)]}^{2}}} d x$ ,

(27)

$a_{2} = a_{3} = \sum_{i = 1}^{n} k_{i} η_{x}^{2} = m \int_{x 1}^{x 2} \frac{- f^{'} (x)}{\sqrt{1 + {[f^{'} (x)]}^{2}}} d x$ ,

(28)

$a_{4} = \sum_{i = 1}^{n} k_{i} η_{y}^{2} = m \int_{x 1}^{x 2} \frac{1}{\sqrt{1 + {[f^{'} (x)]}^{2}}} d x$ ,

(29)

$A = K^{'} [d_{1} + v_{k^{'} x} - \frac{2 (v_{k^{'} x}^{2} + v_{k^{'} y}^{2}) (d_{1} + v_{k^{'} x})}{\sqrt{[{(d_{1} + v_{k^{'} x})}^{2} + {(d_{2} + v_{k^{'} y})}^{2}] \cdot (v_{k^{'} x}^{2} + v_{k^{'} y}^{2})}}],$

(30)

$B = K^{'} [d_{2} + v_{k^{'} y} - \frac{2 (v_{k^{'} x}^{2} + v_{k^{'} y}^{2}) (d_{2} + v_{k^{'} x})}{\sqrt{[{(d_{1} + v_{k^{'} x})}^{2} + {(d_{2} + v_{k^{'} y})}^{2}] \cdot (v_{k^{'} x}^{2} + v_{k^{'} y}^{2})}}],$

(31)

$C = (\frac{1}{2} \sum_{i = 1}^{n} T_{i} + \sum_{i = 1}^{n} T_{k i} - \sum_{i = 1}^{n} T_{sh}) \cos α_{i} + \frac{3}{2} R_{x}$ ,

(32)

$D = (\frac{1}{2} \sum_{i = 1}^{n} T_{i} + \sum_{i = 1}^{n} T_{k i}) \cos β_{i} + \frac{3}{2} R_{y}$ 。

(33)

则式（26）可表示为

$\begin{array}{l} 2 (a_{1} d_{1} + a_{2} d_{2}) + A = C, \\ 2 (a_{4} d_{2} + a_{3} d_{1}) + B = D 。 \end{array}}$

(34)

根据方程组计算得出 $d_{1}$ , $d_{2}$ 。

4.2 抗滑力 $F_{m}$ 的求解

抗滑力 $F_{m}$ 由滑裂面上各微段的法向力N_i、极限剪切力τ_il_i以及水平加筋等效应力或H-V立体加筋应力沿着虚位移方向提供的抗滑力投影之和组成。本文研究的土工格栅对抗滑力的贡献主要体现在增加了滑体的剪切力。整个滑裂面上任意微段的法向力 $N_{i}$ 提供的下滑力为

$\begin{matrix} F_{m1} = \sum_{i = 1}^{n} N_{i} \cos α = \sum_{i = 1}^{n} m (\vec{d} \cdot \vec{n}) l_{i} \cos α \\ = - \sum_{i = 1}^{n} m [\frac{d_{1} f^{'} (x) - d_{2}}{\nabla}] {\vec{d}}_{d} \cdot \vec{n} l_{i} \end{matrix}$ ,

(35)

式中, $T_{s h}$ 为任意微面上法向力与虚位移方向的夹角,且 $\cos α = \vec{n} \cdot (- {\vec{d}}_{d})$ 。

由莫尔-库仑破坏准则,微分圆弧 $d l$ 滑动面上极限抗剪力在虚位移方向上所能提供的抗滑力 $F_{m2}$ 为

$\begin{matrix} F_{m2} = \int_{x 1}^{x 2} d F_{m2} \\ = \int_{x 1}^{x 2} [c + m (\vec{d} \cdot \vec{n}) \tan φ] d l \cos β \end{matrix}$ ,

(36)

式中, $β$ 为剪应力与虚位移的夹角, $\cos β = \vec{t} \cdot (- {\vec{d}}_{d})$ 。

任意微段上的H-V立体筋复合抗拉力提供的抗滑力 $F_{m3}$ 为

$F_{m3} = \sum_{i = 1}^{n} T_{i} \cdot \cos ω = \sum_{i = 1}^{n} (T_{st} + T_{sh}) \cdot \cos ω$ ,

(37)

式中, $ω$ 为水平筋与虚位移的夹角, $\cos ω = \vec{T} \cdot (- {\vec{d}}_{d})$ 。

在虚位移方向上总的抗滑力 $F_{m}$ 为

$F_{m} = F_{m1} + F_{m2} + F_{m3}$ 。

(38)

4.3 下滑力 $F_{n}$ 的求解

作用在滑体上的所有力中,只有合外力 $\vec{R}$ 沿着虚位移 $c = 10$ 方向上提供下滑力 $F_{n}$ ,则总下滑力为

$F_{n} = \sum_{i = 1}^{n} {\overset{\leftarrow}{R}}_{i} \cos ϕ = | \vec{R} | \frac{R_{1} d_{1} + R_{2} d_{2}}{Δ Δ^{'}}$ 。

(39)

5. 算例

如图9所示的H-V立体筋均质边坡^[15],边坡坡度为1∶1,边坡高度6 m,土体黏聚力 $c = 10$ kPa,重度 $γ = 18$ kN/m³,内摩擦角 $φ = 20^{\circ}$ 。该边坡采用的H-V立体筋材料为单向拉伸高密度聚乙烯拉筋,筋材的的抗拉刚度为10000 kN/m,土体的地基系数为100 kN/m³,共铺设5层,每层筋材长为6 m,第一层距坡顶0.6 m,最后一层距坡底0.6 m,竖筋间距D=20 cm,筋材之间的间距为1.2 m,竖筋高度 $H$ =10 mm。

图 9 H-V立体筋边坡示意图

Figure 9. Schematic diagram of reinforced slope

下载: 全尺寸图片幻灯片

建立如图9所示直角坐标系,通过4.5H法搜索滑动面圆心线,通过改变边坡坡度,得到不同的滑动面的方程,利用MATLAB编程试算求解出安全系数的数值最小的滑裂面即为最危险滑动面。则试算得出的边坡最危险滑动面圆心 $o$ 的坐标为（-1.02,10.72）,半径R=10.77 m,滑动面方程AC为

${(x + 1.02)}^{2} + {(y - 10.72)}^{2} = {10.77}^{2}$

。

5.1 算例计算结果分析

本文采用改进最小势能法计算未加筋、加水平筋及加H-V立体筋边坡的安全系数,并与文献[14,15]的水平条分瑞典法的计算结果进行比较,计算结果如表1所示。

表 1 与文献计算结果对比

Table 1. Comparison between calculated results and those in published papers

项目	改进最小势能法	文献[14]的方法	文献[15]的水平条分瑞典法
未加筋	1.178	1.204	1.215
加水平筋	1.563	1.576	1.749
加立体筋	1.728	1.719	2.121

下载: 导出CSV

| 显示表格

由表1可见,对未加筋边坡,本文采用的改进最小势能法计算的安全系数与文献[14,15]的计算结果相差最大值为0.037,相差最小值0.026,相对误差小于3.05%,证明本文的计算方法是合理可行的。对于加水平筋,加立体筋边坡稳定性提高了14%,说明立体筋加筋效果更明显。对比文献[15]的水平条分瑞典法计算结果,本文的计算结果明显偏小,原因在于对加筋边坡稳定性分析时考虑了单个条块水平条间合力为0和整体力矩平衡的条件,考虑了竖向条间力作用,使计算结果偏大。

5.2 不同土体参数算例计算结果对比

目前,文献[15]采用的水平条分法还是一种相对新的计算方法,为了更好地验证本文计算方法的可靠性,针对算例在表2中设置不同的参数与极限平衡法（Bishop法和Janbu法）在无筋（不加筋）、加筋（加水平筋和H-V立体筋）状态下的计算结果进行对比。

表 2 算例参数

Table 2. Parameters of examples

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

下载: 导出CSV

| 显示表格

根据以往文献资料表明,当边坡为单一地层且不考虑加固措施时,通过最小势能法计算边坡稳定性系数过程中可将地基系数m约去^[8],即单一均质边坡稳定性不受地基系数的影响。因为极限平衡法分析边坡稳定性参数中没有地基系数m,为避免地基系数测试结果取值误差产生干扰,以及为后文验证H-V立体筋对提高边坡稳定性的作用,设置表3进行改进最小势能法和极限平衡法无筋边坡计算结果对比。

表 3 无筋边坡算例计算结果对比

Table 3. Comparison of calculated results of examples of unreinforced slopes

对比分组	表2编号	改进最小势能法	极限平衡法		备注
对比分组	表2编号	改进最小势能法	Bishop	Janbu	备注
A	1	1.178	1.117	1.098	$φ$ 渐增
	2	1.180	1.144	1.137
	3	1.182	1.179	1.160
	4	1.184	1.207	1.192
B	5	1.117	1.242	1.231	c渐增
	6	1.238	1.337	1.329
	7	1.344	1.415	1.411
	8	1.594	1.588	1.584

下载: 导出CSV

| 显示表格

从表3计算结果可知,其他参数保持不变,设置两组对比方式,A组内摩擦角 $φ$ 逐渐增大,B组黏聚力c逐渐增大,从计算结果可以看出无筋边坡利用改进最小势能法与极限平衡法计算结果相对误差保持在10%范围以内,证明本文分析计算边坡稳定性的方法是合理可行的,且改进最小势能法计算结果偏大于Janbu法,偏小于Bishop法。

由表4可看出：①改进最小势能法计算结果跟极限平衡法对比,施加水平筋相对误差小于4.7%,施加立体筋相对误差小于2.2%,证明改进最小势能法与极限平衡方法的计算结果具有很好的一致性。②加筋边坡的安全系数均随土体强度指标的增加呈现递增的趋势,其他参数保持不变,对照组C内摩擦角 $φ$ 逐渐增大,对照组D黏聚力c逐渐增大,分析数值变化规律可知时D组加筋边坡安全系数值变化更大,可知黏聚力c对边坡稳定性影响较大。③综合表3,4计算结果还可分析得出,对比无筋边坡,施加水平筋边坡的安全系数提高33.4%,对比施加水平筋,施加立体筋边坡的安全系数提高16.5%,说明施加立体筋在提高边坡稳定性方面的效果更显着。在不同的条件下,采用改进最小势能法的计算结果对比极限平衡法（Bishop法和Janbu法）还存在微小差异,主要因为本文的计算方法与极限平衡法采用的计算模型不一样,且对安全系数的定义和计算仍存在差异性。

表 4 加筋边坡算例计算结果对比

Table 4. Comparison of calculated results of examples of reinforced slope

对比分组	表2编号	改进最小势能法		极限平衡法				备注
		改进最小势能法		Bishop		Janbu
		水平筋	H-V筋	水平筋	H-V筋	水平筋	H-V筋
C	1	1.563	1.749	1.585	1.735	1.559	1.723	φ渐增
	2	1.565	1.761	1.613	1.762	1.558	1.758
	3	1.567	1.784	1.645	1.791	1.561	1.782
	4	1.570	1.792	1.678	1.830	1. 576	1.804
D	5	1.630	1.919	1.646	1.927	1.627	1.930	$c$ 渐增
	6	1.824	2.078	1.844	2.105	1.815	2.111
	7	1.879	2.125	1.902	2.173	1.880	2.123
	8	1.979	2.221	2.020	2.250	1.972	2.218

下载: 导出CSV

| 显示表格

6. 结论

本文在改进最小势能边坡稳定性分析基础上,提出H-V立体筋边坡稳定性分析新方法,并通过算例验证该方法的有效性和合理性,得出以下两点结论。

（1）本文基于改进的最小势能算法,考虑水平筋的摩擦势能、线弹性应变能以及竖筋变形势能的影响,进一步分析了加筋边坡的剪切势能,并将其表示为弹性势能的函数,简化了计算过程。通过计算边坡最危险滑裂面沿虚位移方向静力平衡方程,得到适用于施加水平筋和H-V立体筋的边坡稳定性分析方法。

（2）利用改进最小势能法求解文献算例,将其计算结果与极限平衡法（Bishop法和Janbu法）的计算的结果相比较,相对误差在10%以内,表明本文分析方法是可行且合理的。相比施加水平筋,施加立体筋时边坡安全系数提高了14%,说明H-V立体筋在提高边坡稳定性方面效果显着。

图 1 H-V加筋均质边坡计算图示

Figure 1. Calculation of H-V ribs-reinforced homogeneous slope

下载: 全尺寸图片幻灯片

图 2 H-V立体筋三维结构示意图

Figure 2. Three-dimensional structural diagram of H-V ribs

下载: 全尺寸图片幻灯片

图 3 筋条的特征坐标计算简图

Figure 3. Simplified calculation of characteristic coordinates of geogrid

下载: 全尺寸图片幻灯片

图 4 水平筋复合抗拉力计算示意图

Figure 4. Schematic diagram of calculation of horizontal ribs

下载: 全尺寸图片幻灯片

图 5 土体单元水平筋受力分析图

Figure 5. Analysis diagram of horizontal reinforcement of soil element

下载: 全尺寸图片幻灯片

图 6 立体筋材结构受力示意图

Figure 6. Schematic stress diagram of three-dimensional rib structure

下载: 全尺寸图片幻灯片

图 7 竖筋单元受力图

Figure 7. Force diagram of three-dimensional tendons

下载: 全尺寸图片幻灯片

图 8 边坡土条受力图

Figure 8. Force diagram of slope strips

下载: 全尺寸图片幻灯片

图 9 H-V立体筋边坡示意图

Figure 9. Schematic diagram of reinforced slope

下载: 全尺寸图片幻灯片

表 1 与文献计算结果对比

Table 1 Comparison between calculated results and those in published papers

项目	改进最小势能法	文献[14]的方法	文献[15]的水平条分瑞典法
未加筋	1.178	1.204	1.215
加水平筋	1.563	1.576	1.749
加立体筋	1.728	1.719	2.121

下载: 导出CSV

表 2 算例参数

Table 2 Parameters of examples

编号	黏聚力 $c$ $c$ /kPa	重度 $γ$ $γ$ /(kN·m^-3)	内摩擦角 $φ$ $φ$ /(°)	地基系数 $m$ $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

下载: 导出CSV

表 3 无筋边坡算例计算结果对比

Table 3 Comparison of calculated results of examples of unreinforced slopes

对比分组	表2编号	改进最小势能法	极限平衡法		备注
对比分组	表2编号	改进最小势能法	Bishop	Janbu	备注
A	1	1.178	1.117	1.098	$φ$ $φ$ 渐增
	2	1.180	1.144	1.137
	3	1.182	1.179	1.160
	4	1.184	1.207	1.192
B	5	1.117	1.242	1.231	c渐增
	6	1.238	1.337	1.329
	7	1.344	1.415	1.411
	8	1.594	1.588	1.584

下载: 导出CSV

表 4 加筋边坡算例计算结果对比

Table 4 Comparison of calculated results of examples of reinforced slope

对比分组	表2编号	改进最小势能法		极限平衡法				备注
		改进最小势能法		Bishop		Janbu
		水平筋	H-V筋	水平筋	H-V筋	水平筋	H-V筋
C	1	1.563	1.749	1.585	1.735	1.559	1.723	φ渐增
	2	1.565	1.761	1.613	1.762	1.558	1.758
	3	1.567	1.784	1.645	1.791	1.561	1.782
	4	1.570	1.792	1.678	1.830	1. 576	1.804
D	5	1.630	1.919	1.646	1.927	1.627	1.930	$c$ $c$ 渐增
	6	1.824	2.078	1.844	2.105	1.815	2.111
	7	1.879	2.125	1.902	2.173	1.880	2.123
	8	1.979	2.221	2.020	2.250	1.972	2.218

下载: 导出CSV

参考文献(16)

[1]	VIDAL H. The principles of reinforced earth[J]. Highway Research Record. 1969, 28(2): 1-16.
[2]	SONG Fei, CHEN Ru-yi, MA Li-qiu, et al. A new method for the stability analysis of geosynthetic-reinforced slopes[J]. Journal of Mountain Science, 2016, 13(11): 2069-2078. doi: 10.1007/s11629-016-4001-8
[3]	Bhattacherjee DIPANKANA, VISWANADHAM B V S. Effect of geocomposite layers on slope stability under rainfall condition[J]. Indian Geotechnical Journal, 2018, 48(2): 316-326. doi: 10.1007/s40098-017-0280-4
[4]	EVİRGEN Burak, Mustafa TUNCAN, TUNCAN Ahmet. Modelling study on the geotextile, geogrid and steel strip reinforced slopes[J]. Çukurova Üniversitesi Mühendislik- Mimarlık Fakültesi Dergisi, 2017: 227-240.
[5]	张孟喜. 立体加筋土: ZL200520043841.7[P]. 2009-09. ZHANG Meng-xi. Three-dimensional reinforced soil: ZL200520043841.7[P]. 2009-09. (in Chinese)
[6]	张孟喜, 张石磊. H-V加筋土性状的颗粒流细观模拟[J]. 岩土工程学报, 2008, 30(5): 626-631. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200805001.htm ZHANG Meng-xi, ZHANG Shi-lei. Meso-simulation of particle flow characteristics of H-V reinforced soil[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(5): 626-631. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200805001.htm
[7]	SONG Fei, CHEN Ru-yi, MA Li-qiu, et al. A new method for the stability analysis of geosynthetic-reinforced slopes[J]. Journal of Mountain Science, 2016, 13(11): 2069-2078. doi: 10.1007/s11629-016-4001-8
[8]	孙加平. 改进最小势能边坡稳定性分析方法及程序开发[D]. 南昌: 江西理工大学, 2016. SUN Jia-ping. Analysis Method and Procedure Developing of Slope Stability Based on Improved Minimum Potential Energy[D]. Nanchang: Jiangxi University of Technology, 2016. (in Chinese)
[9]	SONG Fei, CHEN Ru-yi. Stability analysis of reinforced slope based on limit equilibrium method[J]. Tehnicki Vjesnick- Technical Gazette, 2018, 25(1): 224-229.
[10]	邹冶来. 土工格栅加筋边坡的最小势能稳定性分析方法程序开发[D]. 南昌: 江西理工大学, 2019. ZOU Ye-lai. Analysis Method and Program Development of Minimum Potential Energy Stability of Geogrid-Reinforced Slope[D]. Nanchang: Jiangxi University of Science and Technology, 2019. (in Chinese)
[11]	李铀, 陆洋. 锚杆(索)加固边坡的最小势能稳定分析方法研究[J]. 岩土力学, 2008, 29(9): 2330-2334. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200809006.htm LI You, LU Yang. Research on stability analysis method of minimum potential energy of slope reinforced by anchor (cable)[J]. Rock and Soil Mechanics: 2008, 29(9): 2330-2334. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200809006.htm
[12]	孙加平, 顾厚煜. 锚杆作用下边坡稳定性分析中改进最小势能法研究[J]. 岩土力学, 2017, 38(1): 292-298. SUN Jia-ping, GU Hou-yu. An improved minimum potential energy method for slope stability analysis under the action of bolt[J]. Rock and Soil Mechanics, 2017, 38(1): 292-298. (in Chinese)
[13]	温树杰, 宋亮亮, 梁超, 等. 考虑桩体变形的最小势能边坡稳定性分析方法[J]. 工业建筑, 2017, 47(6): 90-94. https://www.cnki.com.cn/Article/CJFDTOTAL-GYJZ201706019.htm WEN Shu-jie, SONG Liang-liang, LIANG Chao, et al. Analysis method of slope stability with minimum potential energy considering pile deformation[J]. Industrial Construction, 2017, 47(6): 90-94. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GYJZ201706019.htm
[14]	张孟喜, 张石磊. H-V加筋路堤稳定性分析及机理研究[D]. 上海: 上海大学, 2010. ZHANG Meng-xi, ZHANG Shi-lei. The Investigation on Stability and Mechanism of Embankment Reinforced with H-V inclusions[D]. Shanghai: Shanghai University: 2010. (in Chinese)
[15]	马学宁, 吴培元, 王旭, 等. 加筋土边坡稳定性分析水平条分简化计算方法[J]. 铁道学报, 2017, 39(9): 156-160. https://www.cnki.com.cn/Article/CJFDTOTAL-TDXB201709022.htm MA Xue-ning, WU Pei-yuan, WANG Xu, et al. Simplified calculation methods of the horizontal slice for the stability analysis of reinforced soil slope[J]. Journal of the China Railway Society, 2017, 39(9): 156-160. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TDXB201709022.htm
[16]	刘动, 陈晓平, 常学宁. 基于筋-土界面摩擦特性的加筋边坡分析[J]. 岩土工程学报, 2013, 35(增刊2): 217-223. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2013S2037.htm LIU Dong, CHEN Xiao-ping, CHANG Xue-ning. Analysis of reinforced slope based on interface properties of geosynthetics and soil[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(S2): 217-223. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2013S2037.htm

施引文献(7)

期刊类型引用(4)

1.	邢亮. 基于判别分析法和机器学习对边坡稳定性的预测研究. 广东土木与建筑. 2024(12): 25-29+50 . 百度学术
2.	郝广杰. 考虑土层变形特征的边坡稳定性分析方法研究. 粉煤灰综合利用. 2023(02): 62-67 . 百度学术
3.	温树杰，赖光甜，孙自立. 改进的边坡最小势能法及潜在滑动方向研究. 铁道科学与工程学报. 2022(08): 2249-2258 . 百度学术
4.	何军. 深层岩土开采中边坡水土流失控制方法研究. 环境科学与管理. 2022(09): 148-153 . 百度学术

其他类型引用(3)

资源附件(0)

图(9) / 表(4)

计量

文章访问数: 257
HTML全文浏览量: 23
PDF下载量: 132
被引次数: 7

0. 引言
1. 滑体储存的弹性势能
2. H-V立体筋势能构建
2.1 筋条特征坐标的确定
2.2 水平筋弹性势能
2.3 竖筋弹性势能
3. 滑面剪切势能的构建
4. 安全系数的计算
4.1 虚位移d→ $\vec{d}$ 的求解
4.2 抗滑力Fm $F_{m}$ 的求解
4.3 下滑力Fn $F_{n}$ 的求解
5. 算例
5.1 算例计算结果分析
5.2 不同土体参数算例计算结果对比
6. 结论

0. 引言
1. 滑体储存的弹性势能
2. H-V立体筋势能构建
2.1 筋条特征坐标的确定
2.2 水平筋弹性势能
2.3 竖筋弹性势能
3. 滑面剪切势能的构建
4. 安全系数的计算
4.1 虚位移 $\vec{d}$ 的求解
4.2 抗滑力 $F_{m}$ 的求解
4.3 下滑力 $F_{n}$ 的求解
5. 算例
5.1 算例计算结果分析
5.2 不同土体参数算例计算结果对比
6. 结论

参考文献(16)

施引文献(7)

资源附件(0)

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

编号	黏聚力 $c$ $c$ /kPa	重度 $γ$ $γ$ /(kN·m^-3)	内摩擦角 $φ$ $φ$ /(°)	地基系数 $m$ $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

编号	黏聚力 $c$ $c$ /kPa	重度 $γ$ $γ$ /(kN·m^-3)	内摩擦角 $φ$ $φ$ /(°)	地基系数 $m$ $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

基于改进最小势能法的水平–竖向立体筋边坡稳定性分析研究 English Version

作者简介: 蒋建国(1973—),男,湖南洞口人,副教授,主要从事路基路面工程、岩土工程方面的研究工作。E-mail: yuanchencsu@126.com

计量

出版历程

Stability analysis of horizontal-vertical ribs-reinforced slopes based on improved minimum potential energy method

0. 引言

1. 滑体储存的弹性势能

2. H-V立体筋势能构建

2.1 筋条特征坐标的确定

2.2 水平筋弹性势能

2.3 竖筋弹性势能

3. 滑面剪切势能的构建

4. 安全系数的计算

4.1 虚位移→d<math xmlns="http://www.w3.org/1998/Math/MathML" id="M163"><mover accent="true"><mi>d</mi><mo>→</mo></mover></math>的求解

4.2 抗滑力Fm<math xmlns="http://www.w3.org/1998/Math/MathML" id="M176"><mrow><msub><mi>F</mi><mtext>m</mtext></msub></mrow></math>的求解

4.3 下滑力Fn<math xmlns="http://www.w3.org/1998/Math/MathML" id="M193"><mrow><msub><mi>F</mi><mi>n</mi></msub></mrow></math>的求解

5. 算例

5.1 算例计算结果分析

5.2 不同土体参数算例计算结果对比

6. 结论

期刊类型引用(4)

其他类型引用(3)

计量

出版历程

目录

0. 引言

1. 滑体储存的弹性势能

2. H-V立体筋势能构建

2.1 筋条特征坐标的确定

2.2 水平筋弹性势能

2.3 竖筋弹性势能

3. 滑面剪切势能的构建

4. 安全系数的计算

4.1 虚位移d→<math xmlns="http://www.w3.org/1998/Math/MathML" id="M163"><mover accent="true"><mi>d</mi><mo>→</mo></mover></math>的求解

4.2 抗滑力Fm<math xmlns="http://www.w3.org/1998/Math/MathML" id="M176"><mrow><msub><mi>F</mi><mtext>m</mtext></msub></mrow></math>的求解

4.3 下滑力Fn<math xmlns="http://www.w3.org/1998/Math/MathML" id="M193"><mrow><msub><mi>F</mi><mi>n</mi></msub></mrow></math>的求解

5. 算例

5.1 算例计算结果分析

5.2 不同土体参数算例计算结果对比

6. 结论

作者简介:
蒋建国(1973—),男,湖南洞口人,副教授,主要从事路基路面工程、岩土工程方面的研究工作。E-mail: yuanchencsu@126.com

4.1 虚位移 $\vec{d}$ 的求解

4.2 抗滑力 $F_{m}$ 的求解

4.3 下滑力 $F_{n}$ 的求解

4.1 虚位移 $\vec{d}$ 的求解

4.2 抗滑力 $F_{m}$ 的求解

4.3 下滑力 $F_{n}$ 的求解

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

编号	黏聚力 $c$ $c$ /kPa	重度 $γ$ $γ$ /(kN·m^-3)	内摩擦角 $φ$ $φ$ /(°)	地基系数 $m$ $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200