滑坡安全系数 Bishop 简化法模拟器

参数设置

斜坡高度 H

斜坡角度 β

粘聚力 c'

kPa

有效应力法中的粘聚力

内部摩擦角 φ'

单位体积重量 γ

kN/m³

孔隙水压比 r_u

降雨和地下水的影响。0.5 时水位上升至斜坡表面

地震系数 k_h

拟静力分析的水平震度。第一级地震动为0.10～0.15

破坏形式

本工具主要进行Bishop圆弧滑动分析

计算结果

—

安全系数 Fs (Bishop)

—

限界斜坡角 (Fs=1) (°)

—

临界 r_u (Fs=1)

—

抵抗力矩

—

驱动力矩

—

风险评估

—

斜坡断面·圆弧滑动面·条分分割

蓝线=圆弧滑动面，虚线=条分边界，浅蓝色带=地下水位，箭头=驱动力和抵抗力。颜色表示 Fs 状态（绿色=稳定／橙色=需注意／红色=破坏危险）。

降雨影响 — Fs vs 孔隙水压比 r_u

地震影响 — Fs vs 地震系数 k_h

理论·主要公式

$$F_s = \frac{\sum \left[c'b + (W - ub)\tan\phi'\right]/m_\alpha}{\sum W\sin\alpha + k_h W\cos\alpha},\quad m_\alpha = \cos\alpha + \frac{\sin\alpha \tan\phi'}{F_s}$$

Bishop简化法的迭代求解公式。c'=有效粘聚力 [kPa]，φ'=有效内部摩擦角 [°]，W=条分重量 [kN/m]，u=条分底面孔隙水压 [kPa]，b=条分宽度 [m]，α=条分底面倾斜角，k_h=水平震度。由于 Fs 作为未知数出现在右边，需要通过3～5次迭代收敛。

$$u = r_u \cdot \gamma \cdot h, \qquad r_{u,\text{crit}} \approx \frac{c'/(\gamma H) + \tan\phi'\cos^2\beta - \sin\beta\cos\beta}{\tan\phi'\cos^2\beta}$$

孔隙水压比 r_u 的定义和使 Fs=1 的临界 r_u（无限斜坡近似）。降雨使地下水位上升时，u 增加，有效应力减少，摩擦阻力随之降低。

滑坡安全系数 Fs — Bishop 简化法的斜坡稳定性分析

🙋

在新闻中经常听到"大雨导致斜坡的安全系数下降"之类的说法，具体是怎样计算的呢？

🎓

简单来说，安全系数 Fs = 抵抗力 / 驱动力，衡量斜坡抵抗下滑的能力。Fs > 1 表示抵抗力还占优势，相对安全；Fs < 1 表示驱动力已占上风，随时可能崩塌。实际工程中，我们假设斜坡沿着一个圆形的"滑动面"整体滑下来，这叫"圆弧滑动"，用Bishop简化法（1955年）计算是最普遍的。

🙋

Bishop先生的方法啊。与普通的力的平衡有什么区别呢？

🎓

斜坡被切分成许多"条分"（这个工具用10条），逐条分计算重量、水压和摩擦力。Bishop法的特点是，它不仅考虑条分侧面的水平力，还把安全系数 Fs 本身作为未知数迭代求解。公式右边包含 Fs 作为 m_α，所以先猜 Fs=1.5，算一下，得到新的 Fs，再用新的 Fs 继续算...这样反复3～5次就收敛了。本工具自动做10次迭代，非常稳定。

🙋

我把滑块上的"孔隙水压比 r_u"从0.3改到0.6，Fs就从1.5暴跌到1.0了！这就是"雨水引发滑坡"吧？

🎓

完全正确。降雨会抬升地下水位，滑动面上的孔隙水压 u 增加。Bishop公式里有个(W - ub)·tan φ'项，这是"有效应力造成的摩擦阻力"，u越大这一项就越小，Fs就急速下降。日本滑坡事件的70～80%都在梅雨季和台风季集中发生，原因就在这儿。r_u = 0.5 大约对应"地下水位升至斜坡表面"，这时原本 Fs=1.5 的斜坡很容易跌破 Fs=1.0。所以降雨预警和避难指令都是基于这个 r_u 来制定阈值的。

🙋

还有一个"地震系数 k_h"。我把它设为0.15，驱动力矩马上飙升，Fs也跌了。地震确实也会引发滑坡呢。

🎓

这就是"拟静力分析"。真实的地震是一个随时间变化的波形，但在工程设计中，我们简化为"水平方向恒定施加相当于重力 k_h 倍的惯性力"。日本《道路土工指南》推荐第一级地震动 k_h=0.10～0.15，第二级 0.15～0.20。美国USACE也是0.10～0.15。用本工具试试 k_h=0.15，驱动力矩增加30～40%，Fs 降至静力值的70%左右。虽然不如 FLAC、PLAXIS Dynamic 那样的详细动力分析精确，但作为设计初期的筛选判定足够了。Vajont大坝（意大利，1963年，约2000人遇难）和台湾八八水灾（2009年，小林村被毁）都是储水位上升加暴雨／台风的复合触发因素。

🙋

Fs=1.5就算安全了吧？安全系数的标准是多少呢？

🎓

日本《宅地造成等规制法》和《斜坡稳定性基准》的一般目标是：常时（无地震、常规水位）Fs ≥ 1.5，地震时 Fs ≥ 1.0，降雨时 Fs ≥ 1.2。但要注意，土质参数（c'、φ'、γ）的波动很大。室内试验测出 c'=30kPa，但实地斜坡可能从20～50kPa都有。所以即使 Fs=1.5，如果参数往坏了偏20%，Fs 就能掉到1.0以下。实务中还要做"可靠性分析"（基于破坏概率来评估），对策工也得"从Fs=1.2提升到1.5"这样有余量地设计。

常见问题

Bishop简化法（1955年）是求圆弧滑动面安全系数 Fs 的极限平衡法之一。与Fellenius法（简易分割法）只考虑各条分力矩不同，Bishop法考虑条分侧面的水平力，并将 Fs 作为未知数进行迭代求解。通常Bishop法给出的 Fs 比Fellenius法大5～20%，对于圆弧滑动与Spencer法或Morgenstern-Price法的结果基本相同。本工具使用10条分进行10次迭代，在典型条件下3～5次迭代即可收敛。

孔隙水压比 r_u 是条分底面的孔隙水压 u 与上覆土的竖向全应力 γh 的比值，是一个无量纲量（r_u = u / (γh)）。地下水位上升至斜坡表面，在无流动的静水压条件下的极限情况下，r_u ≈ 0.5；在有流动时可达0.5～0.7。在本工具的默认条件下（H=20m, β=30°, c'=30kPa, φ'=25°），当 r_u = 0.3 时 Fs ≈ 1.5，当 r_u = 0.6（暴雨后）时 Fs ≈ 1.0 左右。这就是"降雨引发滑坡"的物理原因。

在拟静力分析中，地震动引起的水平加速度 a 与重力 g 的比值 k_h = a/g 被作为水平惯性力加在各条分上。日本《道路土工指南》（NEXCO）中，第一级地震动推荐 k_h = 0.10～0.15，第二级 0.15～0.20。美国USACE推荐0.10～0.15。在本工具中输入 k_h = 0.15 时，驱动力矩增加约30～40%，Fs 降至静力值的70%左右。虽然不如详细的动态响应分析（FLAC、PLAXIS Dynamic）精确，但对设计初期的拟静力判定足够使用。

Fs = 1.0 是"刚好平衡"的状态，在实际工程中绝不认为是安全的。日本《宅地造成等规制法》和《斜坡稳定性基准》的一般目标值为：常时（无地震、常规水位）Fs ≥ 1.5，地震时 Fs ≥ 1.0，降雨时 Fs ≥ 1.2。本工具的"风险评估"也采用以下标准：Fs > 1.5 为稳定，1.2～1.5 为需要注意，1.0～1.2 为警戒，< 1.0 为破坏危险。此外，土质参数具有很大的不确定性，即使 Fs = 1.5，如果 c' 或 φ' 的实际值比设计值低20%，Fs 也会迅速跌破1.0，这一点需要特别注意。

实际应用

陡峭地的防灾和灾害地图：日本根据《土砂災害防止法》将全国陡峭地划分为"土砂災害警戒区域（黄色区域）"和"特别警戒区域（红色区域）"。这些区域的划定和住宅开发的许可，都以本工具所用的Bishop简化法 Fs 计算为基本评估指标。国土地理院的土砂灾害灾害地图和各都道府县公开的Web灾害地图背后，都在进行包含地形、土质数据的数千～数万个地点的 Fs 评估。

道路和铁路的路堑和填土设计：NEXCO（高速公路）和各JR公司的土构造物中，路堑斜坡的 Fs 必须满足：常时 ≥ 1.5，地震时 ≥ 1.0，降雨时 ≥ 1.2。当 Fs 不足时，对策工程从轻到重依次是：(1)排水工（横向钻孔、暗管）以降低 r_u，(2)表面处理（砂浆喷涂、植被），(3)锚杆工和抑止桩，(4)抑止填土。用本工具模拟"将 r_u 从0.5降至0.3"就能定量看出排水工的 Fs 恢复效果。

大坝和水库的斜坡稳定性：水位的快速上升和快速下降（Rapid Drawdown）是典型的滑坡触发因素。1963年意大利Vajont大坝滑坡事件（约2亿立方米岩石滑入水库，约2000人遇难）就是这个现象的典范。大坝上下游的斜坡必须在4种工况下进行 Fs 审查：稳定运行时、蓄水期、地震时和快速下降时。其中快速下降时，孔隙水压的消散跟不上水位下降速度，Fs 会达到最小值。本工具中r_u保持高位而增大β的感觉，对应的就是快速下降时的不稳定化。

滑坡监测和早期预警：近年来，在斜坡上安装IoT倾斜计、GNSS、伸缩计等实时监测设备日益普遍。将观测到的位移速度和降雨量与本工具这样的 Fs 计算相结合，制定"累积降雨达到 X mm 且位移速度超过 Y mm/h 时发布避难勧告"这样的阈值。2014年广岛土砂灾害、2021年热海伊豆山土石流之后，这类AI整合预警系统的投资急剧增加。

常见误区和注意点

最大的陷阱是"Fs计算隐藏了土质参数的不确定性"。输入c'=30kPa到本工具就能得出一个唯一的Fs，但实地斜坡因地点、深度、季节等因素，c'可能从10～60kPa波动，室内三轴试验与原位测试差异可达2～3倍。即使计算出Fs=1.5，如果参数向坏处偏低20%，Fs就会跌破1.0。实务中应该用"可靠性分析"（基于破坏概率的Monte Carlo仿真）并行，对每个参数标注敏感性分析和影响度，不能仅凭一个确定性的Fs值判断"OK"。

其次是"误认为Bishop法能解所有斜坡"。Bishop是圆弧滑动专用，地层倾斜的岩石边坡的平面滑动、沿断层或软弱层的楔形崩塌、流动性强的土石流就不适用。地层与斜坡表面平行堆积的风化岩和崩积土斜坡，平面滑动或Janbu广义分割法更合适。复杂地质需要FEM/FDM（PLAXIS、FLAC）进行弹塑性详细分析。本工具的模式选项虽然有"平面滑动"和"楔形"，但其实都还是用Bishop模型跑的（供教学参考足够）。

最后是"绝对化拟静力分析的k_h"。用k_h=0.15和Fs≥1.0来判定"抗震OK"虽然便利，但真实地震是数十秒的反复动荷载，会产生累积位移、液化、残余强度下降等多重效应。1995年阪神淡路大地震、2011年东日本大地震、2018年北海道胆振东部地震中，拟静力Fs>1的斜坡也观测到了由累积位移引发的破坏。重要构造物邻近或大规模边坡应并行Newmark法（位移法）或动态FEM分析。本工具的k_h仅作"初期筛选"定位。

滑坡安全系数 Fs — Bishop 简化法的斜坡稳定性分析

常见问题

实际应用

常见误区和注意点

使用指南

具体计算例

实务中的注意事项