第11章程序设计339 11.3边坡稳定分析的源程序 11.3.1概述 本章介绍边坡稳定分析程序STAB的核心步骤和程序。程序 CIRCLE. FOR具体实施 1122节的计算步骤,对给定的边坡剖面信息,构造一个圆弧滑裂面。程序 SPLINE FOR具 体实施112、3节的计算步骤,对给定的滑面控制点信息构造一个任意形状滑裂面。程序 PROFILE FOR具体实施11.2.1节的计算步骤,计算土条重量和物理力学参数。程序SFOR 计算安全系数,包含了第2章和第3章介绍的各种分析方法的源程序,即极限平衡法的通用 条分法、瑞典法、毕肖普法、陆军工程师团法和简化法。将上述几个子程序串联起来,就是 一个计算复杂几何形状边坡稳定的安全系数的程序 在介绍这些程序时,将配以图114所示的一个例子。这一例子是将在114节中介绍的 澳大利亚 ACADS程序考核的EX1(c)。该例已在第2章和第4章中多次使用。读者在移植这 些程序时,可通过这些考核题进行考核 图11.4滑裂面和条分 11.3.2构造圆弧滑裂面的程序 CIRCLE. FOR 1.说明 本子程序对给定的边坡剖面信息,构造一个圆弧滑裂面(图114)。通过子程序SEAR找 到圆弧滑裂面与边坡的上、下交点所在线段。子程序DIVI计算上、下交点的坐标,并划分 土条。 2.源程序 REAL YTENSION, DS, CX, CY, XN(99), YN(99), X(80), Y(80), ALF(80) INTEGER INI, LSL, N, LNUM(80), IC(80,3), TWRS, L, NN, NLOW, NUPP OPEN(S, FILE=, STATUS=UNKNOWN) OPEN(6, FILE=,STATUS=UNKNOWN ILOW=0 NUPP=0 READGS, )IWRS, YTENSION, DS, CX,CY IWR5=0,1分别相应顶部有、无拉力缝; YTENSION=拉力缝底y坐标 DS. CX.CY=滑弧深度,圆心坐标 READ(5,*)N土条总数 READ(,*)INl!线段总数 DO 3021=L.INI READ5,*)(C(D)J=1,3)线段信息 302 CONTINUE第 11 章 程序设计 339 11. 3 边坡稳定分析的源程序 11. 3. 1 概述 本章介绍边坡稳定分析程序 STAB 的核心步骤和程序 程序 CIRCLE.FOR 具体实施 11.2.2 节的计算步骤 对给定的边坡剖面信息 构造一个圆弧滑裂面 程序 SPLINE.FOR 具 体实施 11.2.3 节的计算步骤 对给定的滑面控制点信息构造一个任意形状滑裂面 程序 PROFILE.FOR 具体实施 11.2.1 节的计算步骤 计算土条重量和物理力学参数 程序 S.FOR 计算安全系数 包含了第 2 章和第 3 章介绍的各种分析方法的源程序 即极限平衡法的通用 条分法 瑞典法 毕肖普法 陆军工程师团法和简化法 将上述几个子程序串联起来 就是 一个计算复杂几何形状边坡稳定的安全系数的程序 在介绍这些程序时 将配以图 11.4 所示的一个例子 这一例子是将在 11.4 节中介绍的 澳大利亚 ACADS 程序考核的 EX1(c) 该例已在第 2 章和第 4 章中多次使用 读者在移植这 些程序时 可通过这些考核题进行考核 图 11. 4 滑裂面和条分 11. 3. 2 构造圆弧滑裂面的程序 CIRCLE.FOR 1. 说明 本子程序对给定的边坡剖面信息 构造一个圆弧滑裂面(图 11.4) 通过子程序 SEAR 找 到圆弧滑裂面与边坡的上 下交点所在线段 子程序 DIVI 计算上 下交点的坐标 并划分 土条 2. 源程序 REAL YTENSION, DS, CX, CY, XN(99), YN(99), X(80), Y(80), ALF(80) INTEGER IN1, LSL, N, LNUM(80), IC(80,3), IWR5, L, NN, NLOW, NUPP OPEN(5,FILE='',STATUS='UNKNOWN') OPEN(6,FILE='',STATUS='UNKNOWN') NLOW=0 NUPP=0 READ(5,*)IWR5,YTENSION,DS,CX,CY C IWR5=0,1 分别相应顶部有 无拉力缝 YTENSION=拉力缝底 Y 坐标 C DS,CX,CY=滑弧深度 圆心坐标 READ(5,*)N!土条总数 READ(5,*)IN1!线段总数 DO 302 I=1,IN1 READ(5,*)(IC(I,J),J=1,3)!线段信息 302 CONTINUE