(b) (c) 图2滑动面构造步骤.(a)三角形网格：(b)网格节点调整：(c)滑动面与坡面相交 Fig.2 Construction steps of the slip surface:(a)Triangle mesh;(b)Mesh node adjustment()Slip surface intersects the slope surface 滑动面三角形网格中心点应力张量由有限元网格节点处应力采量插值计算，插值方法选用三维线性 插值。本文方法安全系数计算流程如图3所示。 Calculate stress by finite element Calculate the stress tensor on the sliding surface by interpolation Calculate the normal stress and shear stress vector by Eq.2 Cuate the main sliding direction by Eq.11 录用 Calculate the sliding direction of each points by Eq.10 Calculate local and global safety factors by Eq.12 and Eq.15 图3安全系数计算流程 Fig.3 Flowchart of safety factor calculation 5算例 算例中将采用式(1)剪应力比形式定义的安全系数记作F1,将采用式(15)定义的安全系数记作 F2。 5.1算例1 算例I来自Fredlund和Krahn2对二维边坡的研究，有学者将这一算例拓展到三维情况并进行了研究。 该算例包含两种工况，一种是均质边坡，滑动面为椭球形；一种是带软弱薄层的边坡，滑动面为椭球型 和位于软弱薄层的平面组合。椭圆的长轴为133.8m,有限元模型y轴方向的长度为110m,其他相关参图 2 滑动面构造步骤. (a)三角形网格；(b)网格节点调整；(c)滑动面与坡面相交 Fig.2 Construction steps of the slip surface: (a) Triangle mesh; (b) Mesh node adjustment; (c) Slip surface intersects the slope surface 滑动面三角形网格中心点应力张量由有限元网格节点处应力张量插值计算，插值方法选用三维线性 插值。本文方法安全系数计算流程如图 3 所示。 图 3 安全系数计算流程 Fig.3 Flowchart of safety factor calculation 5 算例 算例中将采用式（1）剪应力比形式定义的安全系数记作 F1，将采用式（15）定义的安全系数记作 F2。 5.1 算例 1 算例 1 来自 Fredlund 和 Krahn[25]对二维边坡的研究，有学者将这一算例拓展到三维情况并进行了研究。 该算例包含两种工况，一种是均质边坡，滑动面为椭球形；一种是带软弱薄层的边坡，滑动面为椭球型 和位于软弱薄层的平面组合。椭圆的长轴为 133.8 m，有限元模型 y 轴方向的长度为 110 m，其他相关参 (b) (c) 录用稿件，非最终出版稿