0

GEO5

Ayuda en línea

Tree
Settings
Programa:
Idioma:

Spencer法

Spencer法是根据极限平衡发展起来的一种通用条分法。它需要同时满足作用在各条块上的力和力矩的平衡条件。通过土条间接触面将滑面以上的土体划分为若干条块。作用在各条块上的力如下图所示。

静力图 – Spencer法

每个条块都假设承受以下作用力:

Wi

-

条块重量,包括具有重量的实体超载并考虑了竖向地震加速度系数 Kv

Kh*Wi

-

反应地震作用的水平惯性力,Kh地震时的水平加速度系数

Ni

-

滑面上的法向力

Ti

-

滑面上的剪切力

Ei,Ei+1

-

条块间的作用力,与水平面的夹角为 δ

Fxi,Fyi

-

其它作用在条块上的水平或竖向力

M1i

-

Fxi, Fyi 对点M作用的力矩,点 M i 块条块对应滑面段的中心点

Ui

-

i 块条块对应滑面段上的孔隙水压力合力值

为了求解各条块的极限作用力和力矩平衡方程组,在Spencer法中做出了如下假设:

  • 各条块间的土条间接触面是垂直的
  • 条块重量 Wi 的作用线穿过 i 块条块对应滑面段的中心点,即 M
  • 法向力 Ni 作用在 i 块条块对应滑面段的中心点处,即 M 点处
  • 所有条块受到的条块间作用力 Ei 的倾斜角度都是常量,都为 δ。只有在滑面两个端点处 δ=0

此方法中将使用到以下方程式:

(1)

(2)

(3)

(4)

(5)

方程(1)表示的是作用在滑面上的法向力有效应力值和总应力值之间的关系。方程(2)为表征某滑面段上法向力和剪切力关系的Mohr-Coulomb条件。方程(3)为 i 滑面段法线方向上的作用力平衡方程,而方程(4)表示的是沿 i 滑面段切线方向的受力平衡。SF 是安全系数,用来折减土的强度参数值。方程(5)对应的是条块上作用力关于点 M 取矩的力矩平衡方程,其中 ygi 是条块重力作用点的垂直坐标,yM 是点 M 的垂直坐标。整理方程(3)和(4)得到下列递归公式(6):

通过给定的作用力倾角值δi和安全系数值 SF,便可以根据该方程计算出条块间作用力Ei。且该方法假设坡脚条块左侧的作用力E等于零,即E1 = 0

由力矩平衡方程式(5)得到另一个递归方程(7):

对于给定的 δ 值,可以由上式计算得到所有条块间作用力的力臂值 z。其中坡脚条块左侧的力臂值z1=0

采用以下的迭代过程计算安全系数 SF

  1. δ 的初始值设为 δ = 0
  2. 假设坡顶条块右侧上的作用力 En+1 = 0,则对于给定的 δ 值,可以由方程式(6)得到安全系数 SF
  3. 使用在前面步骤中确定的 E 值,根据坡脚条块上作用的合力矩值为 0 的条件,由方程式(7)求得一个新的δ值。方程(7)并没有给出 zn+1 的值,而是假设等于零。在此取值条件下,必须满足力矩平衡方程(5)。
  4. 重复步骤2和3,直到 δ 的值不再改变为止。

为了使迭代过程收敛,就必须避免不稳定解。当表达式(6)或(7)中的分母为零时便会出现这种不稳定求解的情况。在式(7)中,当 δ = π/2δ = -π/2 时,分母为零。所以角 δ 的取值范围为(-π/2; π/2)。

当满足下式时,方程式(6)中的除数为零:

另一项防止迭代不稳定性的检查是对参数 mα 的验算 - 必须满足以下条件:

因此在迭代计算之前,必须找出满足上述条件的最小临界值(SFmin)。使用小于临界值 SFmin 的数值计算后会得到不稳定解,因此将 SF 设置为仅比 SFmin 大一微小增量的值进行迭代,从而保证迭代过程中得到的所有 SF 的值都大于 SFmin

一般而言,严格的计算方法比简单方法(Bishop、Fellenius)的收敛性要差。存在收敛问题的情况有:滑动面区段太陡,几何模型复杂,超载突变等。如果得不到计算结果,建议稍微改变一下输入的数据,比如:滑面不要太陡,在滑面上确定更多的点位等。或者使用一些简单的分析方法。

参考文献:

Spencer, E. 1967. A method of analysis of the stability of embankments assuming parallel interslice forces. Géotechnique, 17(1): 11–26.

Pruebe GEO5. Gratis.