基于上下限解的膨胀土边坡首次滑动区域分析

贺伟明; 石胜伟; 蔡强; 梁炯

doi:10.16030/j.cnki.issn.1000-3665.202311059

基于上下限解的膨胀土边坡首次滑动区域分析

贺伟明^{1, 2,},
石胜伟^3, ,,
蔡强^{1, 2},
梁炯^{1, 2}

1.
中国地质科学院探矿工艺研究所，四川成都　611734
2.
自然资源部地质灾害风险防控工程技术创新中心，四川成都　611734
3.
中国地质调查局昆明自然资源综合调查中心，云南昆明　650111

基金项目: 国家重点研发计划项目（2019YFC1509904）；中国地质调查局地质调查项目（DD20230451）

详细信息

作者简介:
贺伟明（1992—），男，硕士，工程师,主要从事地质灾害防治方面的研究工作。E-mail：1197037455@qq.com

通讯作者:
石胜伟（1972—），男，本科，教授级高级工程师，主要从事地质灾害形成机理与防治技术研究。E-mail：519903933@qq.com

中图分类号: TU43
计量
- 文章访问数: 266
- HTML全文浏览量: 3
- PDF下载量: 23
出版历程
- 收稿日期: 2023-11-28
- 修回日期: 2024-05-05
- 网络出版日期: 2024-12-04
- 刊出日期: 2025-01-14

First sliding area analysis of expansive soil slope based on the upper and lower bound solutions

HE Weiming^{1, 2,},
SHI Shengwei^3, ,,
CAI Qiang^{1, 2},
LIANG Jiong^{1, 2}

1.
Institute of Exploration Technology, Chinese Academy of Geological Sciences, Chengdu, Sichuan　611734, China
2.
Technology Innovation Center for Risk Prevention and Mitigation of Geohazard, Ministry of Natural Resources, Chengdu, Sichuan　611734, China
3.
Kunming General Survey of Natural Resources Center, China Geological Survey, Kunming, Yunnan　650111, China

HTML全文

摘要

摘要:
受膨胀土胀缩性、裂隙性及超固结性的影响，膨胀土滑坡灾害易发、频发，而传统的计算分析方法难以反映出膨胀土边坡“渐进性”失稳破坏的特点，并有效估算其首次滑动的区域。为揭示膨胀土边坡的成灾特性，量化膨胀土边坡的阶段性成灾规模，基于膨胀土吸水膨胀变形而引起部分土条底切向力方向可能改变的特性，根据Morgenstern-Price静力平衡分析法，严格地导出了膨胀土边坡滑动趋势临界点的位置，分析了边坡首次滑动区域的上限解与下限解，并以某工程为例进行了计算说明，同时分析了抗剪强度对上下限解的影响。结果表明：受膨胀作用影响，按经典分析方法判定为稳定的边坡仍有可能发生局部失稳渐进破坏，其首次滑动区域的上限解受有效抗剪强度的影响较小，而残余抗剪强度对下限解的影响尤为显著；上限解及其变化速率随有效黏聚力与有效内摩擦角的增大而减小，下限解及其变化速率则随残余内摩擦角与残余黏聚力的增加而增加；同时，可基于上下限解的大小关系，判定膨胀土边坡的稳定性。研究成果可为膨胀土滑坡的成灾规模评价及治理工程设计提供参考。
- 膨胀土边坡 /
- 首次滑动区域 /
- 边坡稳定性 /
- 上下限解 /
- Morgenstern-Price法
Abstract:
Expansive soil landslides are prone to frequent disasters due to the expansion, cracking, and over-consolidation characteristics of expansive soil. However, traditional analysis methods cannot adequately reflect the progressive failure characteristic of expansive soil slopes, and effectively estimate the first sliding area. To better understand the disaster characteristics of expansive soil slopes and quantify the phased disaster scale, based on the change of the tangential force direction at the bottom of some soil strips caused by water swelling of expansive soil, the upper and lower bound solutions of the first sliding area of slope are analyzed. The sliding trend critical point of expansive soil slope is strictly derived using the Morgenstern-Price static equilibrium analysis method. A case study is used to illustrate the calculation and analyze the influence of shear strength on the upper solution and lower solutions. The results show that partial and progressive failure may still occur on expansive soil slope under the influence of expansion although they are determined to be stable by classical analysis methods. The upper bound solution of the first slip area is less affected by the effective shear strength, while the influence of residual shear strength on the lower bound solution is particularly significant. The upper bound solution and its change rate decrease with the increase of effective cohesion and effective internal friction angle, while the lower bound solution and its change rate increase with the increase of residual internal friction angle and residual cohesion. In addition, the stability of expansive soil slope can be determined based on the relationship between the upper and lower limit solutions. This study provides the basic information for evaluating disaster scales and designing prevention and control measures for expansive soil landslides.
- expansive soil slope /
- first sliding area /
- stability of slope /
- upper and lower bound solution /
- Morgenstern-Price method

HTML全文

膨胀土作为我国及世界范围内广泛分布的典型特殊性黏土^[1]，具有显著的往复吸水膨胀和失水收缩变形特性^{[2 − 3]}，严重影响着相关地区的工程建设和维护^{[4 − 5]}。特别是二十世纪末以来，随着一批重点项目的谋划及实施^[6]，膨胀土工程问题日益突出，俨然已成为岩土工程领域研究的重点课题、热点问题，同时也是技术难题，尤其是膨胀土边坡的稳定性问题^[7]。

对此，相关学者立足于膨胀土的基本特性，开展了系列模型试验与数值模拟研究，并基于所获现象，分析了不同因素对边坡稳定性的影响，建立了相应的计算方法^{[8 − 9]}。在理论分析方面，不少学者以作用于土体单元的膨胀力反映膨胀性的不利影响，包括仅考虑侧面或底面的膨胀作用，以及同时考虑侧面与底面膨胀作用等情况^[10]。其中，计算时对膨胀力的分布形式也存在不同假设，常见的有三角形分布^[11]、矩形分布^[12]、抛物线型分布^[13]及折线型分布^[14]等。此外，连继峰等^[15]通过概化膨胀土边坡破坏形态为“顺坡曲面”失稳模式，采用低应力下强度非线性幂函数模型，分析了膨胀土路堤边坡浅层稳定性。另一方面，由于理论分析的复杂性及其假设的简化性，采用数值模拟方式分析膨胀土边坡的稳定性更为活跃^{[16 − 17]}。其中，采用湿度应力场理论模拟膨胀作用较为普遍，即通过吸湿膨胀将膨胀应变以附加应变的形式加入到总应变中^{[18 − 20]}。其次，考虑渗流软化作用，以抗剪强度参数与含水率的关系为输入，动态反映因土体含水率变化而引起的强度衰减^{[21 − 22]}。

可见，无论是理论分析还是数值模拟膨胀土边坡的稳定性时，考虑膨胀土的两大主要特性是十分必要的，一是裂隙引起的强度衰减；二是膨胀力的作用。其中，裂隙引起的强度衰减可通过试验有效获取^{[23 − 25]}，膨胀力也可在给定条件下通过室内试验获得^{[26 − 29]}，但其在边坡中作用的机制复杂，且作用模式尚未达成统一。因此，本文基于膨胀土边坡“渐进性”破坏的特点，导出了仅考虑膨胀力宏观影响的边坡首次滑动区域上下限解，以快速估测膨胀土边坡首次滑动的规模，为其危害评估及工程防治提供参考。同时，该方法亦可用于判定膨胀土边坡的稳定性，并界定边坡稳定分析中所施加的膨胀力的上限。

1. 土条底的切向力方向分析

1.1 经典边坡计算方法

诸如Fellenious法、Bishop法、Janbu法等经典边坡稳定分析方法，其重要特点是假设自然作用下所有土条均趋于向坡脚滑动，因而土条底的切向力方向均延续，且与滑动方向相反。

1.2 膨胀土边坡

膨胀土在吸水膨胀时，将沿无外界约束或约束作用较弱的方向产生膨胀变形。以置于水平面上的单位宽度均质膨胀土块为例，其长、高在膨胀初始时分别为 $x (t_{0})$ 、 $y (t_{0})$ ，膨胀终止时分别为 $x (t_{f})$ 、 $y (t_{f})$ ，假设在土块底边存在点s，在膨胀过程中始终处于“静止”状态，且以此点为界，被分割的土块底边长分别为 $x_{1} (t)$ 、 $x_{2} (t)$ ，如图1所示。

图 1 土块膨胀分析模型

Figure 1. Analysis model for cold swelling

下载: 全尺寸图片幻灯片

假设土块底部与其接触面的黏聚力为 $c (t)$ 、内摩擦角为 $φ (t)$ ，土体重度为 $γ (t)$ ，则在膨胀过程中有：

x_{1} (t) [c (t) + y (t) γ \tan φ (t)] = x_{2} (t) [c (t) + y (t) γ \tan φ (t)]

(1)

由式（1）可知， $x_{1} (t) = x_{2} (t)$ ，说明s点为土块底部中点，亦即在膨胀变形过程中，以该点为界，两侧土体的滑动趋势相反。因此，可称s点为滑动趋势临界点（sliding trend critical point，STCP）。

对膨胀土边坡而言，一般地，其失稳破坏多发于降雨后，此时边坡潜在滑体主要趋于临空面（含坡脚、坡顶）方向产生膨胀变形，因而其土条沿潜在滑面的运动趋势将有别于经典边坡分析的假设，具体表现为：受膨胀作用影响，靠近坡脚的部分土条趋于向坡脚滑动，而靠近坡顶的部分土条可能具有向坡顶滑动的趋势，如图2所示。因而，土条底的切向力方向可能在STCP发生突变。

图 2 降雨条件下边坡膨胀变形

注：图中数值表示孔隙水压力。

Figure 2. Swelling deformation of expansive slope under rainfall condition

下载: 全尺寸图片幻灯片

2. 边坡滑动趋势临界点的确定

采用可同时满足力与力矩平衡条件，且适用于任意滑裂面形状的Morgenstern-Price法^[30]进行分析。

2.1 基本假设

基于上述分析，建立膨胀土边坡滑动趋势临界点的分析模型（图3），并做如下假设：

图 3 膨胀土边坡滑动趋势临界点分析模型

Figure 3. Analysis model for the slip tendency critical point at the expansive soil slope

下载: 全尺寸图片幻灯片

1）边坡潜在滑面曲线为 $y_{h} (x)$ ，坡面曲线为 $y_{p} (x)$ ，两曲线相交于o、a两点；

2）s为滑动趋势临界点，以此为界，下部（Ⅰ区）滑体具有向坡脚的滑动趋势，上部（Ⅱ区）滑体反之。

2.2 静力平衡分析

切取Ⅰ区中一垂直土条，建立其力学分析模型，如图4所示。

图 4 Ⅰ区土条受力分析示意图

注：G₁、ΔG₁——土条垂直边上的总作用力及其增量/kN；G_1x、G_1y、(G₁+ΔG₁)_x、(G₁+ΔG₁)_y——各力沿x、y轴方向的分量/kN；β、Δβ——作用力G₁与水平向的夹角及其增量/（°）；y_h、Δy_h——土条底距基准线（x轴）的距离及其增量/m；y_t1、Δy_t1——G₁作用点距基准线（x轴）的距离及其增量/m； ΔQ——水平地震力/kN，其值为ΔQ=ηΔW； η——水平地震影响系数；ΔW——土条自重/kN；ΔN₁——土条底反力/kN；ΔU——土条底受到的总孔隙水压力/kN；ΔT₁——抗滑力/kN，其值为ΔT₁=c_tΔxsecα+(ΔN₁−uΔxsecα)tanφ_t；c_t、φ_t——随土体基质吸力而变化的总黏聚力/kPa及总内摩擦角/（°）；α——土条底倾角/（°）；u——孔隙水压力/（kN∙m⁻¹）。

Figure 4. Schematic diagram of force analysis for soil stripes in area Ⅰ

下载: 全尺寸图片幻灯片

根据静力平衡条件，建立x向和y向的静力平衡方程：

{\begin{aligned} G_{1} \cos β - Δ Q - Δ N_{1} \sin α + Δ T_{1} \cos α - (G_{1} + Δ G_{1}) \cos (β + Δ β) = 0 \\ G_{1} \sin β - Δ W + Δ T_{1} \sin α + Δ N_{1} \cos α - (G_{1} + Δ G_{1}) \sin (β + Δ β) = 0 \end{aligned}

(2)

消去式（2）中的 $Δ N_{1}$ ，并令 $Δ x \to 0$ ，则有微分方程：

\frac{d G_{1}}{d x} - \tan δ_{1} \frac{d β}{d x} G_{1} = P_{1} (x) \sec δ_{1}

(3)

式中：

\begin{aligned} P_{1} (x) = & (c_{t} \cos φ_{t} - u \sin φ_{t}) \sec α + \sin (φ_{t} - α) d W / d x - \\ \cos (φ_{t} - α) d Q / d x, δ_{1} = φ_{t} - α + β 。 \end{aligned}

同时，对土条底中点取矩，建立力矩平衡方程，并简化后有：

G_{1} \sin β + y_{h} \frac{d}{d x} (G_{1} \cos β) = \frac{d}{d x} (y_{t 1} G_{1} \cos β) + \frac{d Q}{d x} h_{e}

(4)

式中： $h_{e}$ ——水平地震力作用点距土条底的距离/m。

以反向的 $Δ T_{2}$ 取代 $Δ T_{1}$ ，建立Ⅱ区土条的力学计算模型，如图5所示。

图 5 Ⅱ区土条受力分析示意图

注：图中变量符号含义参照图4，下标2表示Ⅱ区土条。

Figure 5. Schematic diagram of force analysis for soil stripes in area Ⅱ

下载: 全尺寸图片幻灯片

同理，对Ⅱ区土条建立力的平衡方程，并经整理后有：

\frac{d G_{2}}{d x} + \tan δ_{2} \frac{d β}{d x} G_{2} = P_{2} (x) \sec δ_{2}

(5)

式中：

\begin{aligned} P_{2} (x) = & (u \sin φ_{t} - c_{t} \cos φ_{t}) \sec α - \sin (φ_{t} + α) d W / d x - \\ \cos (φ_{t} + α) d Q / d x, δ_{2} = φ_{t} + α - β 。 \end{aligned}

对Ⅱ区土条建立力矩平衡方程，并简化后有：

G_{2} \sin β + y_{h} \frac{d}{d x} (G_{2} \cos β) = \frac{d}{d x} (y_{t 2} G_{2} \cos β) + \frac{d Q}{d x} h_{e}

(6)

2.3 滑动趋势临界点的确定

式（3）与式（5）分别为关于 $G_{1}$ 、 $G_{2}$ 的一阶非齐次线性微分方程，对其求解，并结合边界条件 $G_{1} (o) = 0$ 、 $G_{2} (a) = 0$ ，有：

G_{1} (x) = \sec δ_{1} S_{1}^{- 1} (x) \int_{x_{o}}^{x} P_{1} (l) S_{1} (l) d l

(7)

G_{2} (x) = \sec δ_{2} S_{2}^{- 1} (x) \int_{x_{a}}^{x} P_{2} (l) S_{2} (l) d l

(8)

其中，

S_{1} (x) = \sec δ_{1} \exp (- \int_{x_{o}}^{x} \tan δ_{1} \frac{d β}{d l} d l)

S_{2} (x) = \sec δ_{2} \exp (\int_{x_{a}}^{x} \tan δ_{2} \frac{d β}{d l} d l)

式中： $x_{o}$ 、 $x_{a}$ ——滑面曲线与坡面曲线交点o、a两点的横坐标/m。

对式（4）、式（6）分别积分，可得：

y_{t 1} (x) = \frac{\int_{x_{o}}^{x} G_{1} (\sin β - \cos β \tan α) d x - \int_{x_{o}}^{x} η \frac{d W}{d x} h_{e} d x}{G_{1} \cos β} + y_{h}

(9)

y_{t 2} (x) = \frac{\int_{x_{a}}^{x} G_{2} (\sin β - \cos β \tan α) d x - \int_{x_{a}}^{x} η \frac{d W}{d x} h_{e} d x}{G_{2} \cos β} + y_{h}

(10)

上述 $β$ 为一未知变量，假定其符合某一形状分布^[31]：

\tan β = f_{0} (x) + λ f (x)

(11)

式中： $f_{0} (x)$ 、 $f (x)$ ——假定的函数分布；

$λ$ ——待求参数。

根据连续性条件，在滑动趋势临界点 $s$ 处有：

G_{1} (s) = G_{2} (s)

(12)

y_{t 1} (s) = y_{t 2} (s)

(13)

$f_{0} (x)$ 、 $f (x)$ 一旦确定，即可联立式（7）—（13），求解 $s$ 及 $λ$ 。

3. 膨胀土边坡首次滑动区域的上下限解

3.1 膨胀土边坡的滑面形态

基于某干渠22处膨胀土滑坡的资料收集与野外调研，其滑面形态及采用二次函数曲线拟合的结果如图6所示。

图 6 膨胀土滑坡滑面形态及其拟合结果

Figure 6. Sliding surface shape of expansive soil landslides and its fitting results

下载: 全尺寸图片幻灯片

由图6分析可知，18处滑面（占比为81.82%）采用二次函数曲线拟合的效果较好，其余4处（SS-2#、SS-3#、XF-4#、AT-3#）滑面的拟合效果则一般。进一步地，采用三次函数曲线对该4处滑面拟合，可获得较为理想的结果。由此可见，采用二次函数曲线可有效拟合绝大部分膨胀土滑坡的滑面，而其余部分则可采用三次函数曲线进行拟合。

在此，为便于说明，以其滑面及坡面均可拟合为二次函数曲线为例，分析边坡首次滑动区域的上下限解。

3.2 首次滑动区域的上限解

假定滑面曲线与坡面曲线分别为：

y_{h} (x) = q_{1} x^{2} + p_{1} x + q_{1} (x_{o} ⩽ x ⩽ x_{a})

(14)

y_{p} (x) = o_{2} x^{2} + p_{2} x + q_{2} (x_{o} ⩽ x ⩽ x_{a})

(15)

式中： $o_{i}$ 、 $p_{i}$ 、 $q_{i}$ （ $i = 1, 2$ ）——拟合参数。

假定 $β (x)$ 为常数^[31]，联立式（7）（14）（15）有：

G_{1} = M_{11} x^{3} + M_{12} x^{2} + M_{13} x - M_{14} ln (1 + \frac{E_{1} x}{H_{1}})

(16)

式中：

M_{11} = A_{1} / (3 E_{1})

M_{12} = B_{1} / (2 E_{1}) - A_{1} H_{1} / (2 E_{1}^{2})

M_{13} = C_{1} / E_{1} - B_{1} H_{1} / E_{1}^{2} + A_{1} H_{1}^{2} / E_{1}^{3}

M_{14} = (- D_{1} E_{1}^{3} + C_{1} H_{1} E_{1}^{2} - B_{1} E_{1} H_{1}^{2} + A_{1} H_{1}^{3}) / E_{1}^{4}

A_{1} = - 2 γ o_{1} (o_{2} - o_{1}) (\cos φ_{t} + η \sin φ_{t})

\begin{aligned} B_{1} = & 4 o_{1}^{2} c_{t} \cos φ_{t} - 2 γ o_{1} (p_{2} - p_{1}) (\cos φ_{t} + η \sin φ_{t}) + \\ γ (o_{2} - o_{1}) (\sin φ_{t} - η \cos φ_{t} - p_{1} \cos φ_{t} - p_{1} η \sin φ_{t}) \end{aligned}

\begin{aligned} C_{1} = & 4 o_{1} p_{1} c_{t} \cos φ_{t} - 2 γ o_{1} (\cos φ_{t} + η \sin φ_{t}) (q_{2} - q_{1}) + \\ γ (\sin φ_{t} - η \cos φ_{t} - p_{1} \cos φ_{t} - p_{1} η \sin φ_{t}) (p_{2} - p_{1}) \end{aligned}

\begin{aligned} D_{1} = & γ (q_{2} - q_{1}) (\sin φ_{t} - η \cos φ_{t} - p_{1} \cos φ_{t} - p_{1} η \sin φ_{t}) + \\ c_{t} \cos φ_{t} (p_{1}^{2} + 1) \end{aligned}

E_{1} = 2 o_{1} \sin (φ_{t} + β)

H_{1} = \cos (φ_{t} + β) + p_{1} \sin (φ_{t} + β)

同理，可求得：

\begin{aligned} G_{2} = & M_{21} (x^{3} - a^{3}) + M_{22} (x^{2} - a^{2}) + M_{23} (x - a) + \\ M_{24} ln (\frac{H_{2} + E_{2} a}{H_{2} + E_{2} x}) \end{aligned}

(17)

式中：

M_{21} = A_{2} / (3 E_{2})

M_{22} = B_{2} / (2 E_{2}) - A_{2} H_{2} / (2 E_{2}^{2})

M_{23} = C_{2} / E_{2} - B_{2} H_{2} / E_{2}^{2} + A_{2} H_{2}^{2} / E_{2}^{3}

M_{24} = (- D_{2} E_{2}^{3} + C_{2} H_{2} E_{2}^{2} - B_{2} E_{2} H_{2}^{2} + A_{2} H_{2}^{3}) / E_{2}^{4}

A_{2} = - 2 γ o_{1} (o_{2} - o_{1}) (\cos φ_{t} - η \sin φ_{t})

\begin{aligned} B_{2} = & - 4 o_{1}^{2} c_{t} \cos φ_{t} - 2 γ o_{1} (p_{2} - p_{1}) (\cos φ_{t} - η \sin φ_{t}) - \\ γ (o_{2} - o_{1}) (\sin φ_{t} + η \cos φ_{t} + p_{1} \cos φ_{t} - p_{1} η \sin φ_{t}) \end{aligned}

\begin{aligned} C_{2} = & - 4 o_{1} p_{1} c_{t} \cos φ_{t} - 2 γ o_{1} (\cos φ_{t} - η \sin φ_{t}) (q_{2} - q_{1}) - \\ γ (\sin φ_{t} + η \cos φ_{t} + p_{1} \cos φ_{t} - p_{1} η \sin φ_{t}) (p_{2} - p_{1}) \end{aligned}

\begin{aligned} D_{2} = & - γ (q_{2} - q_{1}) (\sin φ_{t} + η \cos φ_{t} + p_{1} \cos φ_{t} - p_{1} η \sin φ_{t}) - \\ c_{t} \cos φ_{t} (p_{1}^{2} + 1) \end{aligned}

E_{2} = - 2 o_{1} \sin (φ_{t} - β)

H_{2} = \cos (φ_{t} - β) - p_{1} \sin (φ_{t} - β)

联立式（9）、式（14）—（16）有：

\begin{aligned} y_{t 1} (x) = & \frac{J_{G 1} (\sin β - p_{1} \cos β) - 2 o_{1} \cos β J_{x G 1} - η γ J_{W 1} / 2}{G_{1} \cos β} + \\ o_{1} x^{2} + p_{1} x + q_{1} \end{aligned}

(18)

式中：

\begin{aligned} J_{G 1} = & \frac{M_{11}}{4} x^{4} + \frac{M_{12}}{3} x^{3} + \frac{M_{13}}{2} x^{2} - \\ M_{14} [(\frac{H_{1}}{E_{1}} + x) In (1 + \frac{E_{1} x}{H_{1}}) - x] \end{aligned}

\begin{aligned} J_{x G 1} = & \frac{M_{11}}{5} x^{5} + \frac{M_{12}}{4} x^{4} + \frac{M_{13}}{3} x^{3} - \\ \frac{M_{14}}{2} [(x^{2} - \frac{H_{1}^{2}}{E_{1}^{2}}) In (1 + \frac{E_{1} x}{H_{1}}) - \frac{x^{2}}{2} + \frac{H_{1}}{E_{1}} x] \end{aligned}

\begin{aligned} J_{W 1} = & \frac{{(o_{2} - o_{1})}^{2}}{5} x^{5} + \frac{(o_{2} - o_{1}) (p_{2} - p_{1})}{2} x^{4} + \\ \frac{2 (o_{2} - o_{1}) (q_{2} - q_{1}) + {(p_{2} - p_{1})}^{2}}{3} x^{3} + \\ (p_{2} - p_{1}) (q_{2} - q_{1}) x^{2} + {(q_{2} - q_{1})}^{2} x \end{aligned}

同理，可求得：

\begin{aligned} y_{t 2} (x) = & \frac{J_{G 2} (\sin β - p_{1} \cos β) - 2 o_{1} \cos β J_{x G 2} - η γ J_{W 2} / 2}{G_{2} \cos β} + \\ o_{1} x^{2} + p_{1} x + q_{1} \end{aligned}

(19)

式中：

J_{G 2} = {[\begin{matrix} \frac{M_{21}}{4} x^{4} + \frac{M_{22}}{3} x^{3} + \frac{M_{23}}{2} x^{2} - \\ M_{24} ((\frac{H_{2}}{E_{2}} + x) ln (H_{2} + E_{2} x) - x) - \\ (M_{21} a^{3} + M_{22} a^{2} + M_{23} a - M_{24} ln (H_{2} + E_{2} a)) x \end{matrix}]}_{x_{a}}^{x}

J_{x G 2} = {[\begin{matrix} \frac{M_{21}}{5} x^{5} + \frac{M_{22}}{4} x^{4} + \frac{M_{23}}{3} x^{3} - \\ M_{24} ((\frac{x^{2}}{2} - \frac{F_{2}^{2}}{2 E_{2}^{2}}) ln (H_{2} + E_{2} x) - \frac{x^{2}}{4} + \frac{H_{2}}{2 E_{2}} x) - \\ \frac{M_{21} a^{3} + M_{22} a^{2} + M_{23} a - M_{24} ln (H_{2} + E_{2} a)}{2} x^{2} \end{matrix}]}_{x_{a}}^{x}

J_{W 2} = {[\begin{matrix} \frac{{(o_{2} - o_{1})}^{2}}{5} x^{5} + \frac{(o_{2} - o_{1}) (p_{2} - p_{1})}{2} x^{4} + \\ \frac{2 (o_{2} - o_{1}) (q_{2} - q_{1}) + {(p_{2} - p_{1})}^{2}}{3} x^{3} + \\ (p_{2} - p_{1}) (q_{2} - q_{1}) x^{2} + {(q_{2} - q_{1})}^{2} x \end{matrix}]}_{x_{a}}^{x}

联立式（12）（13）、式（16）—（19）即可求得首次滑动区域的上限解，对应于该处的条间作用力亦可用于界定该处膨胀力的上限。

需要特别说明的是，按经典边坡计算方法求得的边坡整体安全系数小于1的情况，不在该分析方法所包含的范围内。

3.3 首次滑动区域的下限解

受膨胀作用影响，被推动部分土体的抗剪强度将转变为残余强度，故首次滑动区域的下限解（d）基于残余强度的平衡状态求取。仍可采用上述求取上限解的方法按式（16）计算求取，所不同的是，部分参数（c_t变为残余黏聚力c_r，φ_t变为残余内摩擦角φ_r）及边界条件的变化：

G_{1} (d) = 0

(20)

y_{t 1} (d) = y_{h} (d)

(21)

具体计算过程在此不再赘述。

4. 算例分析

某支渠膨胀土边坡如图7所示，地层岩性为：①层上更新统重粉质壤土（ ${Qp}_{3}^{a l}$ ）：夹粉质黏土，灰黄、棕黄色，硬可塑～硬塑，局部坚硬，湿，夹铁锰质结核及灰白色高岭土，中等压缩性；②层下白垩统朱巷组全风化～强风化泥质粉砂岩（ $K_{1} z$ ）：软质岩石，泥质胶结，紫红色，破碎，岩心呈土状，少数块状，揭露最大厚度3.20 m。潜在滑体均处于①层，其湿重度为19.4 kN/m³，场地地下水类型分为两类：一类为上层滞水，分布于①层重粉质壤土表层孔隙中，受大气降水和地表水渗入补给，无统一水位，水量一般较小，多随季节性降水变化而变化，勘探期间，水位位于滑坡鼓丘以下5～7 m，地下水受降水及干渠河水的补给。第二类为基岩裂隙水，分布于下伏基岩泥质粉砂岩中，微承压性，水量小。采用固结快剪试验测得潜在滑面原状土体的总黏聚力 $c_{t}$ =10 kPa，总内摩擦角 $φ_{t} = 19 °$ ；残余黏聚力c_r=2.5 kPa，残余内摩擦角 $φ_{r} = 12 °$ 。此外，水平地震影响系数（ $η$ ）取0.05。

图 7 某支渠膨胀土边坡示意图

Figure 7. Schematic diagram of expansive soil slope along a branch canal

下载: 全尺寸图片幻灯片

拟合后的滑面曲线与坡面曲线（单位均为m）分别为：

y_{h} (x) = 0.012 2 x^{2} + 0.064 x + 0.222 4 (0 ⩽ x ⩽ 23.3)

(22)

y_{p} (x) = - 0.005 6 x^{2} + 0.469 6 x + 0.731 9 (0 ⩽ x ⩽ 23.3)

(23)

经编程计算，边坡整体稳定性系数为1.62，在本次研究范围内，其首次滑动区域的上限解s=17.71 m，下限解d=11.47 m，由此可估算该膨胀土滑坡单位宽度首次滑动的体量（S）介于23.56～39.67 m³，如图8所示。

图 8 首次滑动区域的上下限解

注：此图为滑动区域的平面示意，故S采用面积单位m²。

Figure 8. Upper and lower bound solutions in the first sliding area

下载: 全尺寸图片幻灯片

5. 抗剪强度对上下限解的影响

基于上述算例，分别分析有效抗剪强度对上限解的影响及残余抗剪强度对下限解的影响。

5.1 有效抗剪强度对上限解的影响

为保障边坡不发生整体失稳，取有效黏聚力的变化范围为7～13 kPa，相应地有效内摩擦角的变化范围取13°～25°。经计算整理，得到上限解与有效黏聚力及有效内摩擦角的关系，如图9所示。

图 9 上限解与有效抗剪强度的关系曲面

Figure 9. Relationship between upper bound solution and effective shear strength

下载: 全尺寸图片幻灯片

由图可见，上限解随有效黏聚力及有效内摩擦角的增大而减小，且其减小速率也均随有效黏聚力及有效内摩擦角的增大而减小。当有效黏聚力为7 kPa，有效内摩擦角由13°增至25°时，上限解减小4.44 m；而当有效黏聚力为13 kPa，有效内摩擦角由13°增至25°时，上限解减小2.41 m。同样地，当有效内摩擦角为13°，有效黏聚力由7 kPa增至13 kPa时，上限解减小2.97 m；而当有效内摩擦角为25°，有效黏聚力由7 kPa增至13 kPa时，上限解减小仅为0.94 m。

5.2 残余抗剪强度对下限解的影响

残余黏聚力取0.5～4.5 kPa，相应地，残余内摩擦角取8°～16°。经计算整理，得到下限解与残余黏聚力及残余内摩擦角的关系，如图10所示。

图 10 下限解与残余抗剪强度的关系曲面

Figure 10. Relationship between lower bound solution and residual shear strength

下载: 全尺寸图片幻灯片

图中反映出，下限解随残余内摩擦角及残余黏聚力的增加而增加，且其增加速率也均随残余黏聚力及残余内摩擦角的增加而增加。当残余黏聚力为0.5 kPa，残余内摩擦角由8°增至16°时，下限解增加9.57 m；而当残余内摩擦角为8°，残余黏聚力由0.5 kPa增至4.5 kPa时，下限解增加6.9 m。

特别地，在残余黏聚力及残余内摩擦角均较大的区域（如： $φ_{r}$ =16°， $c_{r}$ ≥3.5 kPa），无下限解，表明该类情况下，边坡不会发生失稳破坏。

综合分析，膨胀土边坡首次滑动区域的上限解受有效抗剪强度的影响较小，而残余抗剪强度对于下限解的影响尤为显著。同时，可进一步基于上下限解的大小关系，分析膨胀土边坡的稳定性：

（1）当下限解小于上限解时，边坡可能发生局部失稳破坏，其首次滑动区域可估算确定；

（2）当下限解等于上限解时，边坡可能发生局部失稳破坏，且其首次滑动区域可准确确定；

（3）当无下限解或下限解大于上限解时，边坡处于稳定状态。

6. 结论

（1）采用二次函数可有效拟合绝大部分膨胀土滑坡的滑面，而其余部分则可进一步地采用三次函数曲线进行拟合。

（2）受膨胀作用影响，按经典分析方法判定为稳定的膨胀土边坡仍有可能发生局部失稳渐进破坏，其首次滑动区域的范围可通过上下限解确定。

（3）膨胀土边坡首次滑动区域的上限解受有效抗剪强度的影响较小，而残余抗剪强度对于下限解的影响尤为显著。上限解随有效黏聚力及有效内摩擦角的增大而减小，且其减小速率也均随有效黏聚力及有效内摩擦角的增大而减小；下限解随残余抗剪强度的变化则相反。

（4）可基于上下限解的大小关系，分析膨胀土边坡的稳定性。当下限解不大于上限解时，边坡可能发生局部失稳破坏；当无下限解或下限解大于上限解时，边坡处于稳定状态。

图 1 土块膨胀分析模型

Figure 1. Analysis model for cold swelling

下载: 全尺寸图片幻灯片

图 2 降雨条件下边坡膨胀变形

注：图中数值表示孔隙水压力。

Figure 2. Swelling deformation of expansive slope under rainfall condition

下载: 全尺寸图片幻灯片

图 3 膨胀土边坡滑动趋势临界点分析模型

Figure 3. Analysis model for the slip tendency critical point at the expansive soil slope

下载: 全尺寸图片幻灯片

图 4 Ⅰ区土条受力分析示意图

注：G₁、ΔG₁——土条垂直边上的总作用力及其增量/kN；G_1x、G_1y、(G₁+ΔG₁)_x、(G₁+ΔG₁)_y——各力沿x、y轴方向的分量/kN；β、Δβ——作用力G₁与水平向的夹角及其增量/（°）；y_h、Δy_h——土条底距基准线（x轴）的距离及其增量/m；y_t1、Δy_t1——G₁作用点距基准线（x轴）的距离及其增量/m； ΔQ——水平地震力/kN，其值为ΔQ=ηΔW； η——水平地震影响系数；ΔW——土条自重/kN；ΔN₁——土条底反力/kN；ΔU——土条底受到的总孔隙水压力/kN；ΔT₁——抗滑力/kN，其值为ΔT₁=c_tΔxsecα+(ΔN₁−uΔxsecα)tanφ_t；c_t、φ_t——随土体基质吸力而变化的总黏聚力/kPa及总内摩擦角/（°）；α——土条底倾角/（°）；u——孔隙水压力/（kN∙m⁻¹）。

Figure 4. Schematic diagram of force analysis for soil stripes in area Ⅰ

下载: 全尺寸图片幻灯片

图 5 Ⅱ区土条受力分析示意图

注：图中变量符号含义参照图4，下标2表示Ⅱ区土条。

Figure 5. Schematic diagram of force analysis for soil stripes in area Ⅱ

下载: 全尺寸图片幻灯片

图 6 膨胀土滑坡滑面形态及其拟合结果

Figure 6. Sliding surface shape of expansive soil landslides and its fitting results

下载: 全尺寸图片幻灯片

图 7 某支渠膨胀土边坡示意图

Figure 7. Schematic diagram of expansive soil slope along a branch canal

下载: 全尺寸图片幻灯片

图 8 首次滑动区域的上下限解

注：此图为滑动区域的平面示意，故S采用面积单位m²。

Figure 8. Upper and lower bound solutions in the first sliding area

下载: 全尺寸图片幻灯片

图 9 上限解与有效抗剪强度的关系曲面

Figure 9. Relationship between upper bound solution and effective shear strength

下载: 全尺寸图片幻灯片

图 10 下限解与残余抗剪强度的关系曲面

Figure 10. Relationship between lower bound solution and residual shear strength

下载: 全尺寸图片幻灯片

参考文献(31)

[1]	SADEGHI H,ALIPANAHI P. Saturated hydraulic conductivity of problematic soils measured by a newly developed low-compliance triaxial permeameter[J]. Engineering Geology,2020,278:105827. DOI: 10.1016/j.enggeo.2020.105827
[2]	JIANG Guanlu,CHEN Wweizhi,LIU Xianfeng,et al. Field study on swelling-shrinkage response of an expansive soil foundation under high-speed railway embankment loads[J]. Soils and Foundations,2018,58(6):1538 − 1552. DOI: 10.1016/j.sandf.2018.09.008
[3]	ZHAN Liangtong,CHEN Ping,NG C W W. Effect of suction change on water content and total volume of an expansive clay[J]. Journal of Zhejiang University:Science A,2007,8(5):699 − 706. DOI: 10.1631/jzus.2007.A0699
[4]	KHAN M S,HOSSAIN S,AHMED A,et al. Investigation of a shallow slope failure on expansive clay in Texas[J]. Engineering Geology,2017,219:118 − 129. DOI: 10.1016/j.enggeo.2016.10.004
[5]	AVSAR E,ULUSAY R,SONMEZ H. Assessments of swelling anisotropy of Ankara clay[J]. Engineering Geology,2009,105(1/2):24 − 31.
[6]	NG C W W,ZHAN L T,BAO C G,et al. Performance of an unsaturated expansive soil slope subjected to artificial rainfall infiltration[J]. Géotechnique,2003,53(2):143 − 157.
[7]	程展林,李青云,郭熙灵,等. 膨胀土边坡稳定性研究[J]. 长江科学院院报,2011,28(10):102 − 111. [CHENG Zhanlin,LI Qingyun,GUO Xiling,et al. Study on the stability of expansive soil slope[J]. Journal of Yangtze River Scientific Research Institute,2011,28(10):102 − 111. (in Chinese with English abstract)] DOI: 10.3969/j.issn.1001-5485.2011.10.019 CHENG Zhanlin, LI Qingyun, GUO Xiling, et al. Study on the stability of expansive soil slope[J]. Journal of Yangtze River Scientific Research Institute, 2011, 28（10）: 102 − 111. （in Chinese with English abstract） DOI: 10.3969/j.issn.1001-5485.2011.10.019
[8]	SHI Bin,JIANG Hongtao,LIU Zhibin,et al. Engineering geological characteristics of expansive soils in China[J]. Engineering Geology,2002,67(1/2):63 − 71.
[9]	ESTABRAGH A R,PARSAEI B,JAVADI A A. Laboratory investigation of the effect of cyclic wetting and drying on the behaviour of an expansive soil[J]. Soils and Foundations,2015,55(2):304 − 314. DOI: 10.1016/j.sandf.2015.02.007
[10]	张震,林宇亮,何红忠,等. 膨胀土边坡的失稳特征与稳定性分析[J]. 中南大学学报(自然科学版),2022,53(1):104 − 113. [ZHANG Zhen,LIN Yuliang,HE Hongzhong,et al. Instability characteristics and stability analysis of expansive soil slope[J]. Journal of Central South University (Science and Technology),2022,53(1):104 − 113. (in Chinese with English abstract)] DOI: 10.11817/j.issn.1672-7207.2022.01.006 ZHANG Zhen, LIN Yuliang, HE Hongzhong, et al. Instability characteristics and stability analysis of expansive soil slope[J]. Journal of Central South University （Science and Technology）, 2022, 53（1）: 104 − 113. （in Chinese with English abstract） DOI: 10.11817/j.issn.1672-7207.2022.01.006
[11]	杨果林,陈子昂,张红日,等. 干湿循环作用下平缓型膨胀土边坡失稳破坏机制研究[J]. 中南大学学报(自然科学版),2022,53(1):95 − 103. [YANG Guolin,CHEN Ziang,ZHANG Hongri,et al. Collapse mechanism of gentle expansive soil slope in drying and wetting cycles[J]. Journal of Central South University (Science and Technology),2022,53(1):95 − 103. (in Chinese with English abstract)] DOI: 10.11817/j.issn.1672-7207.2022.01.005 YANG Guolin, CHEN Ziang, ZHANG Hongri, et al. Collapse mechanism of gentle expansive soil slope in drying and wetting cycles[J]. Journal of Central South University （Science and Technology）, 2022, 53（1）: 95 − 103. （in Chinese with English abstract） DOI: 10.11817/j.issn.1672-7207.2022.01.005
[12]	杨文琦,周成,王林,等. 降雨作用下考虑膨胀推力的膨胀土边坡稳定性分析[J]. 防灾减灾工程学报,2018,38(5):874 − 880. [YANG Wenqi,ZHOU Cheng,WANG Lin,et al. Stability analysis of an expansive soil slope under rainfall using the thrust line method[J]. Journal of Disaster Prevention and Mitigation Engineering,2018,38(5):874 − 880. (in Chinese with English abstract)] YANG Wenqi, ZHOU Cheng, WANG Lin, et al. Stability analysis of an expansive soil slope under rainfall using the thrust line method[J]. Journal of Disaster Prevention and Mitigation Engineering, 2018, 38（5）: 874 − 880. （in Chinese with English abstract）
[13]	孙即超,王光谦,董希斌,等. 膨胀土膨胀模型及其反演[J]. 岩土力学,2007,28(10):2055 − 2059. [SUN Jichao,WANG Guangqian,DONG Xibin,et al. Expansive model and black analysis of expansive soils[J]. Rock and Soil Mechanics,2007,28(10):2055 − 2059. (in Chinese with English abstract)] DOI: 10.3969/j.issn.1000-7598.2007.10.009 SUN Jichao, WANG Guangqian, DONG Xibin, et al. Expansive model and black analysis of expansive soils[J]. Rock and Soil Mechanics, 2007, 28（10）: 2055 − 2059. （in Chinese with English abstract） DOI: 10.3969/j.issn.1000-7598.2007.10.009
[14]	肖世国. 膨胀土堑坡稳定性分析[J]. 岩土力学,2001,22(2):152 − 155. [XIAO Shiguo. Stability analysis of expansive clay slope[J]. Rock and Soil Mechanics,2001,22(2):152 − 155. (in Chinese with English abstract)] DOI: 10.3969/j.issn.1000-7598.2001.02.008 XIAO Shiguo. Stability analysis of expansive clay slope[J]. Rock and Soil Mechanics, 2001, 22（2）: 152 − 155. （in Chinese with English abstract） DOI: 10.3969/j.issn.1000-7598.2001.02.008
[15]	连继峰,罗强,张文生,等. 路堤边坡膨胀土强度非线性应力阈值与浅层稳定性[J]. 哈尔滨工业大学学报,2021,53(3):142 − 151. [LIAN Jifeng,LUO Qiang,ZHANG Wensheng,et al. Nonlinear stress threshold for expansive soil strength and shallow stability of embankment slope[J]. Journal of Harbin Institute of Technology,2021,53(3):142 − 151. (in Chinese with English abstract)] DOI: 10.11918/201908149 LIAN Jifeng, LUO Qiang, ZHANG Wensheng, et al. Nonlinear stress threshold for expansive soil strength and shallow stability of embankment slope[J]. Journal of Harbin Institute of Technology, 2021, 53（3）: 142 − 151. （in Chinese with English abstract） DOI: 10.11918/201908149
[16]	QI Shunchao,VANAPALLI S K. Influence of swelling behavior on the stability of an infinite unsaturated expansive soil slope[J]. Computers and Geotechnics,2016,76:154 − 169. DOI: 10.1016/j.compgeo.2016.02.018
[17]	QI Shuncaho,VANAPALLI S K. Hydro-mechanical coupling effect on surficial layer stability of unsaturated expansive soil slopes[J]. Computers and Geotechnics,2015,70:68 − 82. DOI: 10.1016/j.compgeo.2015.07.006
[18]	张良以,陈铁林,张顶立. 降雨诱发膨胀土边坡渐进破坏研究[J]. 岩土工程学报,2019,41(1):70 − 77. [ZHANG Liangyi,CHEN Tielin,ZHANG Dingli. Progressive failure of expansive soil slopes under rainfall[J]. Chinese Journal of Geotechnical Engineering,2019,41(1):70 − 77. (in Chinese with English abstract)] DOI: 10.11779/CJGE201901007 ZHANG Liangyi, CHEN Tielin, ZHANG Dingli. Progressive failure of expansive soil slopes under rainfall[J]. Chinese Journal of Geotechnical Engineering, 2019, 41（1）: 70 − 77. （in Chinese with English abstract） DOI: 10.11779/CJGE201901007
[19]	赵思奕,石振明,鲍燕妮,等. 考虑吸湿膨胀及软化的膨胀土边坡稳定性分析[J]. 工程地质学报,2021,29(3):777 − 785. [ZHAO Siyi,SHI Zhenming,BAO Yanni,et al. Stability analysis of expansive soil slope considering seepage softening and moistening swelling deformation[J]. Journal of Engineering Geology,2021,29(3):777 − 785. (in Chinese with English abstract)] ZHAO Siyi, SHI Zhenming, BAO Yanni, et al. Stability analysis of expansive soil slope considering seepage softening and moistening swelling deformation[J]. Journal of Engineering Geology, 2021, 29（3）: 777 − 785. （in Chinese with English abstract）
[20]	张坤勇,孙斌,李威,等. 膨胀土边坡失稳机理分析及基于完全软化强度的边坡极限平衡计算方法[J]. 中南大学学报(自然科学版),2022,53(1):75 − 94. [ZHANG Kunyong,SUN Bin,LI Wei,et al. Study on the failure mechanism of expansive soil slope and limit equilibrium analysis method based on fully softening strength[J]. Journal of Central South University (Science and Technology),2022,53(1):75 − 94. (in Chinese with English abstract)] DOI: 10.11817/j.issn.1672-7207.2022.01.004 ZHANG Kunyong, SUN Bin, LI Wei, et al. Study on the failure mechanism of expansive soil slope and limit equilibrium analysis method based on fully softening strength[J]. Journal of Central South University （Science and Technology）, 2022, 53（1）: 75 − 94. （in Chinese with English abstract） DOI: 10.11817/j.issn.1672-7207.2022.01.004
[21]	饶鸿,王金淑,赵志明,等. 基于有限元软件自定义本构模型的膨胀土边坡降雨入渗分析[J]. 水文地质工程地质,2021,48(1):154 − 162. [RAO Hong,WANG Jinshu,ZHAO Zhiming,et al. An analysis of rainfall infiltration of expansive soil slope based on the finite element software custom constitutive model[J]. Hydrogeology & Engineering Geology,2021,48(1):154 − 162. (in Chinese with English abstract)] RAO Hong, WANG Jinshu, ZHAO Zhiming, et al. An analysis of rainfall infiltration of expansive soil slope based on the finite element software custom constitutive model[J]. Hydrogeology & Engineering Geology, 2021, 48（1）: 154 − 162. （in Chinese with English abstract）
[22]	李晋鹏,汪磊,王俊,等. 考虑抗剪强度衰减特性的膨胀土边坡稳定性分析[J]. 中国地质灾害与防治学报,2022,33(6):29 − 36. [LI Jinpeng,WANG Lei,WANG Jun,et al. Stability analysis of expansive soil slopes considering shear strength decay characteristics[J]. The Chinese Journal of Geological Hazard and Control,2022,33(6):29 − 36. (in Chinese with English abstract)] LI Jinpeng, WANG Lei, WANG Jun, et al. Stability analysis of expansive soil slopes considering shear strength decay characteristics[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33（6）: 29 − 36. （in Chinese with English abstract）
[23]	李晶晶,孔令伟,穆坤. 膨胀土原位孔内剪切试验与强度响应特征[J]. 岩土力学,2017,38(2):453 − 461. [LI Jingjing,KONG Lingwei,MU Kun. In-situ borehole shear test on expansive soil and its strength characteristics[J]. Rock and Soil Mechanics,2017,38(2):453 − 461. (in Chinese with English abstract)] LI Jingjing, KONG Lingwei, MU Kun. In-situ borehole shear test on expansive soil and its strength characteristics[J]. Rock and Soil Mechanics, 2017, 38（2）: 453 − 461. （in Chinese with English abstract）
[24]	刘华强,殷宗泽. 裂缝对膨胀土抗剪强度指标影响的试验研究[J]. 岩土力学,2010,31(3):727 − 731. [LIU Huaqiang,YIN Zongze. Test study of influence of crack evolution on strength parameters of expansive soil[J]. Rock and Soil Mechanics,2010,31(3):727 − 731. (in Chinese with English abstract)] DOI: 10.3969/j.issn.1000-7598.2010.03.010 LIU Huaqiang, YIN Zongze. Test study of influence of crack evolution on strength parameters of expansive soil[J]. Rock and Soil Mechanics, 2010, 31（3）: 727 − 731. （in Chinese with English abstract） DOI: 10.3969/j.issn.1000-7598.2010.03.010
[25]	赵鑫,阳云华,朱瑛洁,等. 裂隙面对强膨胀土抗剪强度影响分析[J]. 岩土力学,2014,35(1):130 − 133. [ZHAO Xin,YANG Yunhua,ZHU Yingjie,et al. Analysis of impact of crack surface on shear strength of strong expansive soil[J]. Rock and Soil Mechanics,2014,35(1):130 − 133. (in Chinese with English abstract)] ZHAO Xin, YANG Yunhua, ZHU Yingjie, et al. Analysis of impact of crack surface on shear strength of strong expansive soil[J]. Rock and Soil Mechanics, 2014, 35（1）: 130 − 133. （in Chinese with English abstract）
[26]	CHEN Yonggui,LIU Lina,YE Weimin,et al. Deterioration of swelling pressure of compacted Gaomiaozi bentonite induced by heat combined with hyperalkaline conditions[J]. Soils and Foundations,2019,59(6):2254 − 2264. DOI: 10.1016/j.sandf.2019.12.008
[27]	叶为民,刘樟荣,崔玉军,等. 膨润土膨胀力时程曲线的形态特征及其模拟[J]. 岩土工程学报,2020,42(1):29 − 36. [YE Weimin,LIU Zhangrong,CUI Yujun,et al. Features and modelling of time-evolution curves of swelling pressure of bentonite[J]. Chinese Journal of Geotechnical Engineering,2020,42(1):29 − 36. (in Chinese with English abstract)] DOI: 10.11779/CJGE202001003 YE Weimin, LIU Zhangrong, CUI Yujun, et al. Features and modelling of time-evolution curves of swelling pressure of bentonite[J]. Chinese Journal of Geotechnical Engineering, 2020, 42（1）: 29 − 36. （in Chinese with English abstract） DOI: 10.11779/CJGE202001003
[28]	TANAKA Y,WATANABE Y. Modelling the effects of test conditions on the measured swelling pressure of compacted bentonite[J]. Soils and Foundations,2019,59(1):136 − 150. DOI: 10.1016/j.sandf.2018.09.009
[29]	赖小玲,叶为民,刘毅,等. 高庙子膨润土膨胀力时效性试验研究[J]. 岩土工程学报,2014,36(3):574 − 579. [LAI Xiaoling,YE Weimin,LIU Yi,et al. Experimental investigation on ageing effects on swelling pressure of unsaturated GMZ01 bentonite[J]. Chinese Journal of Geotechnical Engineering,2014,36(3):574 − 579. (in Chinese with English abstract)] DOI: 10.11779/CJGE201403022 LAI Xiaoling, YE Weimin, LIU Yi, et al. Experimental investigation on ageing effects on swelling pressure of unsaturated GMZ01 bentonite[J]. Chinese Journal of Geotechnical Engineering, 2014, 36（3）: 574 − 579. （in Chinese with English abstract） DOI: 10.11779/CJGE201403022
[30]	MORGENSTERN N R,PRICE V E. The analysis of the stability of general slip surfaces[J]. Géotechnique,1965,15(1):79 − 93.
[31]	陈祖煜. 土质边坡稳定分析——原理·方法·程序[M]. 北京:中国水利水电出版社,2003. [CHEN Zuyu. Soil slope stability analysis:Theory methods and programs [M]. Beijing:China Water & Power Press,2003. (in Chinese)] CHEN Zuyu. Soil slope stability analysis: Theory methods and programs [M]. Beijing: China Water & Power Press, 2003. （in Chinese）

施引文献

资源附件(0)

图(10)

计量

文章访问数: 266
HTML全文浏览量: 3
PDF下载量: 23
被引次数: 0

1. 土条底的切向力方向分析
1.1 经典边坡计算方法
1.2 膨胀土边坡
2. 边坡滑动趋势临界点的确定
2.1 基本假设
2.2 静力平衡分析
2.3 滑动趋势临界点的确定
3. 膨胀土边坡首次滑动区域的上下限解
3.1 膨胀土边坡的滑面形态
3.2 首次滑动区域的上限解
3.3 首次滑动区域的下限解
4. 算例分析
5. 抗剪强度对上下限解的影响
5.1 有效抗剪强度对上限解的影响
5.2 残余抗剪强度对下限解的影响
6. 结论

1. 土条底的切向力方向分析
1.1 经典边坡计算方法
1.2 膨胀土边坡
2. 边坡滑动趋势临界点的确定
2.1 基本假设
2.2 静力平衡分析
2.3 滑动趋势临界点的确定
3. 膨胀土边坡首次滑动区域的上下限解
3.1 膨胀土边坡的滑面形态
3.2 首次滑动区域的上限解
3.3 首次滑动区域的下限解
4. 算例分析
5. 抗剪强度对上下限解的影响
5.1 有效抗剪强度对上限解的影响
5.2 残余抗剪强度对下限解的影响
6. 结论

参考文献(31)

施引文献

资源附件(0)

基于上下限解的膨胀土边坡首次滑动区域分析

作者简介:
贺伟明（1992—），男，硕士，工程师,主要从事地质灾害防治方面的研究工作。E-mail：1197037455@qq.com

通讯作者:
石胜伟（1972—），男，本科，教授级高级工程师，主要从事地质灾害形成机理与防治技术研究。E-mail：519903933@qq.com

计量

First sliding area analysis of expansive soil slope based on the upper and lower bound solutions

1. 土条底的切向力方向分析

1.1 经典边坡计算方法

1.2 膨胀土边坡

2. 边坡滑动趋势临界点的确定

2.1 基本假设

2.2 静力平衡分析

2.3 滑动趋势临界点的确定

3. 膨胀土边坡首次滑动区域的上下限解

3.1 膨胀土边坡的滑面形态

3.2 首次滑动区域的上限解

3.3 首次滑动区域的下限解

4. 算例分析

5. 抗剪强度对上下限解的影响

5.1 有效抗剪强度对上限解的影响

5.2 残余抗剪强度对下限解的影响

6. 结论

计量

目录

1. 土条底的切向力方向分析

1.1 经典边坡计算方法

1.2 膨胀土边坡

2. 边坡滑动趋势临界点的确定

2.1 基本假设

2.2 静力平衡分析

2.3 滑动趋势临界点的确定

3. 膨胀土边坡首次滑动区域的上下限解

3.1 膨胀土边坡的滑面形态

3.2 首次滑动区域的上限解

3.3 首次滑动区域的下限解

4. 算例分析

5. 抗剪强度对上下限解的影响

5.1 有效抗剪强度对上限解的影响

5.2 残余抗剪强度对下限解的影响

6. 结论

友情链接：

基于上下限解的膨胀土边坡首次滑动区域分析

作者简介: 贺伟明（1992—），男，硕士，工程师,主要从事地质灾害防治方面的研究工作。E-mail：1197037455@qq.com

通讯作者: 石胜伟（1972—），男，本科，教授级高级工程师，主要从事地质灾害形成机理与防治技术研究。E-mail：519903933@qq.com

计量

出版历程

First sliding area analysis of expansive soil slope based on the upper and lower bound solutions

1. 土条底的切向力方向分析

1.1 经典边坡计算方法

1.2 膨胀土边坡

2. 边坡滑动趋势临界点的确定

2.1 基本假设

2.2 静力平衡分析

2.3 滑动趋势临界点的确定

3. 膨胀土边坡首次滑动区域的上下限解

3.1 膨胀土边坡的滑面形态

3.2 首次滑动区域的上限解

3.3 首次滑动区域的下限解

4. 算例分析

5. 抗剪强度对上下限解的影响

5.1 有效抗剪强度对上限解的影响

5.2 残余抗剪强度对下限解的影响

6. 结论

计量

出版历程

目录

1. 土条底的切向力方向分析

1.1 经典边坡计算方法

1.2 膨胀土边坡

2. 边坡滑动趋势临界点的确定

2.1 基本假设

2.2 静力平衡分析

2.3 滑动趋势临界点的确定

3. 膨胀土边坡首次滑动区域的上下限解

3.1 膨胀土边坡的滑面形态

3.2 首次滑动区域的上限解

3.3 首次滑动区域的下限解

4. 算例分析

5. 抗剪强度对上下限解的影响

5.1 有效抗剪强度对上限解的影响

5.2 残余抗剪强度对下限解的影响

6. 结论

友情链接：

作者简介:
贺伟明（1992—），男，硕士，工程师,主要从事地质灾害防治方面的研究工作。E-mail：1197037455@qq.com

通讯作者:
石胜伟（1972—），男，本科，教授级高级工程师，主要从事地质灾害形成机理与防治技术研究。E-mail：519903933@qq.com