1.一种边坡渐进破坏潜在滑动面的计算方法，其特征在于包括如下步骤：(1)对滑体材料实施剪应力--剪应变全过程曲线试验，试验获得峰值应力、应变和全过程曲线；

(2)通过峰值应力决定凝聚力C、滑面摩擦角 值，通过峰值应变决定常系数a1,a2,a3大小，通过变化曲线特征决定剪切模量G、临界法向应力σncrit、常系数ξ、α、kn；

(3)建立数值计算模型，在考虑剪破坏分布的同时，考虑拉破坏分布区；

(4)在考虑应变软化本构模型数值计算的基础上，计算现状边坡每点的破坏率、破坏面和破坏比，以组合的方式提出不同的可能破坏路径；

(5)对于单元，利用单元破坏剪应力面与最小主应力的夹角为 计算最大主应力相对竖直方向的转动角δ，从而决定滑面相对于水平面的转角 所述转动角δ，二维计算公式为tan2δ＝-2τxy/(σxx-σyy)，三维计算公式为tan2δxx＝-2τxy/(σxx-σyy)，tan2δyy＝-2τzy/(σyy-σzz)，tan2δzz＝-2τzx/(σzz-σxx)；

(6)对于可能施加的荷载或位移工况，分步施加相应的工况，对在不同工况下可能破坏模式进行搜寻，将潜在滑动面转动角连续化，计算对应的边坡稳定系数，从而决定潜在滑动面。

2.如权利要求1所述的计算方法，其特征在于：对于具有软化和硬化特征的滑面剪应力--剪应变满足如下本构方程，(7.1)滑面剪应力-剪应变方程

滑面剪应力--剪应变为四参数本构方程：τ＝Gγ[1+γq/p]ξ (7.1)式中：τ、γ分别为剪应力和剪应变，G为剪切模量，p、q、ξ为在不同法向应力下的常系数，τ、G的单位为MPa或kPa或Pa，p、q、ξ为无单位参数；并将软化和硬化行为描述如下：(7.2)软化特征

对于具有软化特征的材料行为，则有-1＜ξ≤0和1+qξ≠0；临界应变空间满足如下关系式：p+(1+qξ)γqpeak＝0 (7.2)式中：γpeak为临界应力对应的应变；

假设临界应力空间τpeak满足摩尔库伦准则：式中：C为凝聚力，σn为法向应力，C和σn的单位为MPa或kPa或Pa，为滑面摩擦角；

假设临界应变空间仅相关于法向应力，临界应变γpeak采用如下关系式：或

式中：a1,a2,a3,ζN, 为常系数；a1,a2单位为MPa或kPa或Pa，a3,ζN为无量纲系数，

2 2 2

或 的量纲为1/MPa,1/MPa或1/kPa,1/kPa或1/Pa,1/Pa ；

且G＝G0+b1σn+b2σn2 (7.5)式中：G0为法向应力σn为零值的G值，b1,b2为常系数，b1的单位为无量纲，b2的单位为1/MPa或1/kPa或1/Pa；

对于无量纲参数ξ，软化系数演化方程表示为：式中，ξ0为法向应力σn为零值的ξ值，ξc为σn等于σnc时的ξ值，为常系数；

(7.3)硬化特征

当地质材料的法向应力大于临界法向应力 时，则没有明显的峰值应力，采用以下两种计算方法：(7.3.1)方法一

将本构方程(7.1)取ξ＝-1，q＝1，则a’＝1/(Ga”)，b'＝1/(Gp)，其方程形式与邓肯—张模型一致，此时只能描述材料的弹塑性硬化行为特征；

式中：a’、b’、a”为常系数；

在峰值应力条件下，方程(7.7)变为：定义割线模量

则

对方程(7.7)求导，相应的导数为切线模量，在任一应力状态条件下，切向模量Gi表示为：利用方程(7.11)，在最大应力时的切线模量Gt则有：Gi＝a,Kcant2 (7.12)在峰值应力时，研究试验曲线的切向模量，记为Gt，假设其具有如下特征：α、kn为常系数；

方程(7.13)的特点为：

当 切线模量等于0，此时曲线呈现近似理想弹塑性模型特征，当σn达到一定值时，曲线呈现线性特征，此时的法向应力理论上试验可以决定，记为 相对应的切向模量应该等于 则具有如下方程：在法向应力 的范围内，取某一法向应力 试验确定对应的切向模量：Ga，获得如下方程：

则由方程(7.14、7.15)，决定常系数：和

在确定某一法向应力σn条件下的峰值应力切向模量Gt之后，按方程(7.12)决定a’，按方程(7.10)决定b’，至此新的邓肯--张模型各参数得以确定；

(7.3.2)方法二

取方程(7.1)的ξ＝-1，则方程变为：在峰值应力下，则有：

同理对方程(7.17)求导，相应的导数为切线模量：当峰值应力满足目前摩尔库伦准则时，则峰值应变也满足方程(7.4)，在峰值应力时，其切线模量记为Gt；

利用方程(7.18、7.19)，在峰值应力时，切线模量满足方程：对于峰值应力对应的切线模量，按照方程(7.13)求解；利用方程(7.21)求解参数q，利用方程(7.19)求解p。

3.如权利要求2所述的计算方法，其特征在于：针对条分法，潜在滑动面计算方法包括如下子步骤，(1)利用现行方法进行条块划分；

(2)竖向应力以比重与高度之积加以计算，水平应力和剪应力以剩余推力为水平方向和垂直于水平方向力的矢量和，水平方向力和垂直于水平方向力假设满足一定的应力分布条件；

(3)计算条块底边摩阻应力。

4.如权利要求3所述的计算方法，其特征在于：条分法决定潜在滑动面，按强度折减法进行，具体为对条块底边临界抗剪强度进行折减，直至位于临空面的破坏条块处于极限平衡状态。

5.如权利要求3所述的计算方法，其特征在于：条分法决定潜在滑动面，按施加荷载或位移工况方法进行，具体为对可能破坏施加相应的荷载或位移工况，直至位于临空面的破坏条块处于极限平衡状态。