第12章程序使用说明387 g3=16T/m 缆机垂直荷毂分布 =qx=4.4T/ 机水平荷载分布 17937 414490 185125 P=420T P=420T 图12.12输入外荷载例题 在表121第33行输入参数MTE。相应MTE的不同数值,程序执行以下分析方法 (1)有效应力法和坝体施工期总应力法(MTE=0)。根据不同运用阶段,分别采用S剪、 Q剪、R剪的强度指标。施工、竣工期采用Q剪指标进行总应力法分析,此时,不计孔隙水 压力,其处理方法实际上同有效应力法。例题EX9介绍此情况。 (2)静力总应力法(MTE>O)本书第5章论述了总应力法的原理并介绍了详细计算步骤。 STAB程序中使用静力总应力法有以下两种情况: (a)采用现场快剪强度指标。这种情况主要针对地基土具有现场十字板剪力仪或其它现 场土工试验的总强度的情况。这个数据随深度变化。程序要求输入相应现场位置x,y值和 =qa/2的网格,通过内插确定滑裂面上各点总强度。在数据文件中,用MTE>0(表12.1第33 行)和LRU①>0(表121第48行)来控制。例题EX11介绍这个情况 (b)利用固结不排水试验强度包线R或组合强度“R-S”进行库水位骤降时的边坡稳定总 应力法分析。这一步骤已在第7章通过[例7刀]作详细介绍。STAB程序严格按照美国陆军 工程师团建议的方法操作。在数据文件中,用MITE>0(表12.1第33行)以及输入不同的 降前、降后水位UwL和UⅥLI(表12.1第34行)和输入R-S组合包线来实现此功能。执 行此功能,还需输入库水位降落后坝内的浸润线,表12.1中81~83行。例题EX8介绍此情 况
第 12 章 程序使用说明 387 图 12. 12 输入外荷载例题 在表 12.1 第 33 行输入参数 MTE 相应 MTE 的不同数值 程序执行以下分析方法 (1) 有效应力法和坝体施工期总应力法(MTE=0) 根据不同运用阶段 分别采用 S 剪 Q 剪 R 剪的强度指标 施工 竣工期采用 Q 剪指标进行总应力法分析 此时 不计孔隙水 压力 其处理方法实际上同有效应力法 例题 EX9 介绍此情况 (2) 静力总应力法(MTE>0) 本书第 5 章论述了总应力法的原理并介绍了详细计算步骤 STAB 程序中使用静力总应力法有以下两种情况 (a) 采用现场快剪强度指标 这种情况主要针对地基土具有现场十字板剪力仪或其它现 场土工试验的总强度的情况 这个数据随深度变化 程序要求输入相应现场位置 x y 值和τ =qcu/2 的网格 通过内插确定滑裂面上各点总强度 在数据文件中 用 MTE>0(表 12.1 第 33 行)和 LRU(I)>0 表 12.1 第 48 行 来控制 例题 EX11 介绍这个情况 (b) 利用固结不排水试验强度包线 R 或组合强度 R−S”进行库水位骤降时的边坡稳定总 应力法分析 这一步骤已在第 7 章通过[例 7.7]作详细介绍 STAB 程序严格按照美国陆军 工程师团建议的方法操作 在数据文件中 用 MTE>0 表 12.1 第 33 行 以及输入不同的 降前 降后水位 UWL 和 UWL1 表 12.1 第 34 行 和输入 R−S 组合包线来实现此功能 执 行此功能 还需输入库水位降落后坝内的浸润线 表 12.1 中 81∼83 行 例题 EX8 介绍此情 况
88土质边被穗定分原理方法程序 (c)地震总应力法(MTE<0)。采用第7章介绍的步骤利用动三轴试验的总强度进行边坡 稳定计算的。 土的动抗剪强度vs由下式确定: t s=f(o3,K) 上二式中:t为破坏面上的剪应力;o1为固结大主应力;o3为固结小主应力 用程序进行计算,需将三轴不排水动力试验所得结果整理成式(125)形式,然后以网格 形式,逐点将σ3、K和τ输入,如图12.10所示。对某一σ3、K,程序通过内插确定τs 本程序规定具备动指标的土料种类不得超过五种,对需进行动力分析的土层依次由1 至LSD编号,不具备动指标的土层编号为零。称这组编号为“动编号”,用LSDM(I),I=1 IN表 例如某问题共包含5层土,IN=5。如果编号为1、3、4的三层土具备动指标,其相应 的τs,K数据网格的编号分别为1、2、3。第2、5层土不做地震总应力分析。那么LSD=3, 这五层土的动编号为LSDM()=1,0,2,3,0(1=1,,5)。也就是说,程序将识别,第2、5层土 无动指标,第1层土将使用下面输入的第一组动强度数据网。第3、4层土分别使用第2、3 组动强度数据网 程序在READ2(表123)中输入各层有关数据。例题EX12介绍此功能。 2.稳定计算方法 (1)简化法。程序提供第3章介绍的全部简化法的计算功能。 用于圆弧滑裂面的简化法,按网格布置多个圆弧滑裂面,程序同时算出毕肖普法、瑞典 法和改良瑞典法(罗厄法或工程师团法)的结果,分别通过LLl=0,1,2来实现(表12.2第 23行),并找出相应于某一种方法最小安全系数值的临界滑弧 罗厄法和陆军工程师团法是满足静力平衡的适用于任意形状滑裂面的简化法,用户在两 者中选一种。在数据表121第34行中,当ASP<0时,为罗厄法;ASP≥0时,为工程师团 法,此时ASP为边坡的平均坡度(以度计),此值仅在工程师团法中用。罗厄法和陆军工程 师团法也可在任意形状滑裂面条件下使用。程序先计算罗厄法或陆军工程师团法,再以此为 初值,获得通用条分法的成果 (2)边坡稳定分析通用条分法。本法仅用于任意形状滑裂面。方法的要点是输入对土条 侧向力倾角的一个分布形状,然后根据力和力矩平衡确定安全系数,已在第2章中详细介绍。 本程序默认的功能是 Spencer法。当OPl(3)>0时,程序执行第2章介绍的对β(x)的其它假定 此时,按READ4中不同处理,对fx)和fx)输入一个由用户指定的分布函数。在OPl(3)=0 时无需输入其它参数,即自动使用 Spencer法计算。在OPl(3)大于零时,通过表122的77~91 行指定fx)和fx)的数值。 OPl(3)大于零时首先读入控制码ISPE。 ISPE=-2,采用第2章2.52节介绍的侧向力的第2种假定,即固定的两端值,f(x)采用
388 土质边坡稳定分析 原理 ⋅ 方法 ⋅ 程序 (c) 地震总应力法(MTE<0) 采用第 7 章介绍的步骤利用动三轴试验的总强度进行边坡 稳定计算的 土的动抗剪强度τfs由下式确定 ( , ) (12.5) fs 3 Kc τ = f σ 3 1 σ σ Kc = (12.6) 上二式中 τfs 为破坏面上的剪应力 σ1 为固结大主应力 σ3 为固结小主应力 用程序进行计算 需将三轴不排水动力试验所得结果整理成式(12.5)形式 然后以网格 形式 逐点将σ3 Kc和τfs输入 如图 12.10 所示 对某一σ3 Kc 程序通过内插确定τfs 本程序规定具备动指标的土料种类不得超过五种 对需进行动力分析的土层依次由 1 至 LSD 编号 不具备动指标的土层编号为零 称这组编号为 动编号 用 LSDM(I) I=1 IN 表示 例如某问题共包含 5 层土 IN = 5 如果编号为 1 3 4 的三层土具备动指标 其相应 的τfs Kc数据网格的编号分别为 1 2 3 第 2 5 层土不做地震总应力分析 那么 LSD=3 这五层土的动编号为 LSDM(I)=1,0,2,3,0 (I=1,…,5) 也就是说 程序将识别 第 2 5 层土 无动指标 第 1 层土将使用下面输入的第一组动强度数据网 第 3 4 层土分别使用第 2 3 组动强度数据网 程序在 READ2 表 12.3 中输入各层有关数据 例题 EX12 介绍此功能 2. 稳定计算方法 (1) 简化法 程序提供第 3 章介绍的全部简化法的计算功能 用于圆弧滑裂面的简化法 按网格布置多个圆弧滑裂面 程序同时算出毕肖普法 瑞典 法和改良瑞典法 罗厄法或工程师团法 的结果 分别通过 LL1 = 0,1,2 来实现 表 12.2 第 23 行 并找出相应于某一种方法最小安全系数值的临界滑弧 罗厄法和陆军工程师团法是满足静力平衡的适用于任意形状滑裂面的简化法 用户在两 者中选一种 在数据表 12.1 第 34 行中 当 ASP<0 时 为罗厄法 ASP 0 时 为工程师团 法 此时 ASP 为边坡的平均坡度 以度计 此值仅在工程师团法中用 罗厄法和陆军工程 师团法也可在任意形状滑裂面条件下使用 程序先计算罗厄法或陆军工程师团法 再以此为 初值 获得通用条分法的成果 (2) 边坡稳定分析通用条分法 本法仅用于任意形状滑裂面 方法的要点是输入对土条 侧向力倾角的一个分布形状 然后根据力和力矩平衡确定安全系数 已在第 2 章中详细介绍 本程序默认的功能是 Spencer 法 当 OP1(3)>0 时 程序执行第 2 章介绍的对β(x)的其它假定 此时 按 READ4 中不同处理 对 f(x)和 f0(x)输入一个由用户指定的分布函数 在 OP1(3)=0 时无需输入其它参数 即自动使用 Spencer 法计算 在 OP1(3)大于零时 通过表 12.2 的 77∼91 行指定 f(x)和 f0(x)的数值 OP1(3)大于零时首先读入控制码 ISPE ISPE=−2 采用第 2 章 2.5.2 节介绍的侧向力的第 2 种假定 即固定的两端值 f(x)采用
第12拿程序使用说明 正弦曲线 ISPE=0,由用户指定x)和fx) 其次,读入BOUN。当BOUN为零时两端值按第2章252节提供的判据式(2.90)由程 序自动设置,不为零时,则由用户指定端值,用FAA和FBB定义。 ISPE=-2和0的两种情况的实例,参见EX27和EX29。 表127总结了本程序进行边坡稳定分析的各种功能 3.边坡稳定的可靠度分析 程序执行第10章介绍的蒙特卡洛法、一次二阶矩法、 Rothenbleuth法三种可靠度分析 的功能。输入数据在子程序READ5(表126)中执行 首先输入两个参数 PROB和NLP PROB=1:蒙特卡洛法;=2:一次二阶矩法;=3: Rosenbluth法 NLP=具有随机特性土层总数,如果 NPROB=1则需单独为蒙特卡洛法输入以下参数 IIM=蒙特卡洛法模拟次数 K=具有随机特性的变量总数 MEED=产生随机数的种子 随后对上述三种情况,输入以下参数 MTYPE(I)=1:均匀分布;=2:正态分布;=3:负指数分布;=4:对数正态分布 FMEAND①,FDEV(I), RMEAN(),RDEV()分别为c,中,r的均值和标准差 KK=0,不相关:KK=1,相关。 在KK=1时输入相关系数矩阵Rl1(,J 与确定性方法中所用的分析过程类似,边坡稳定概率分析通常包括两步即对于一指定滑 面求β;在可能的滑面中,寻找一个所谓临界滑面,其β为最小值βm 注意蒙特卡洛法无法考虑参数相关,也不能进行自动寻找临界滑面的优化功能。 4.土压力计算 当OPl(1)为2和3时进行土压力计算,此时在表12的第95行输入以下参数: GWALL=假定的土压力大小。假定值不影响最终成果,但不应与预计的最终成果偏离 过大 HMW=土压力作用位置与墙底的垂直距离 EWALL=土压力与水平线的夹角 当OPl(1)=2时进行同时满足力和力矩平衡条件的土压力计算 当OPl(1)=3时进行仅满足力平衡条件的土压力计算,即库仑土压力理论计算。此时计 算成果与输入的HMW值无关。 12.2.5应用最优化方法确定临界滑裂面 1.圆弧滑裂面 程序提供的搜索圆弧形临界滑裂面的功能有3种,分别用表12.2中23行LL2=0,1,2来 实现。在一般情况下建议直接使用LL2=2,即下述(3)的功能
第 12 章 程序使用说明 389 正弦曲线 ISPE=0 由用户指定 f0(x)和 f(x) 其次 读入 IBOUN 当 IBOUN 为零时两端值按第 2 章 2.5.2 节提供的判据式(2.90)由程 序自动设置 不为零时 则由用户指定端值 用 FAA 和 FBB 定义 ISPE=−2 和 0 的两种情况的实例 参见 EX27 和 EX29 表 12.7 总结了本程序进行边坡稳定分析的各种功能 3. 边坡稳定的可靠度分析 程序执行第 10 章介绍的蒙特卡洛法 一次二阶矩法 Rothenbleuth 法三种可靠度分析 的功能 输入数据在子程序 READ5 表 12.6 中执行 首先输入两个参数 NPROB 和 NLP PROB=1 蒙特卡洛法 =2 一次二阶矩法 =3 Rosenbleuth 法 NLP= 具有随机特性土层总数, 如果 NPROB=1 则需单独为蒙特卡洛法输入以下参数 ITM = 蒙特卡洛法模拟次数 K= 具有随机特性的变量总数 MEED(I)= 产生随机数的种子 随后对上述三种情况 输入以下参数 MTYPE(I)=1 均匀分布 =2 正态分布 =3 负指数分布 =4 对数正态分布 FMEAN(I) FDEV(I) RMEAN(I) RDEV(I)分别为 c, φ, ru的均值和标准差 KK=0 不相关 KK=1 相关 在 KK=1 时 输入相关系数矩阵 R11(I, J) 与确定性方法中所用的分析过程类似 边坡稳定概率分析通常包括两步即对于一指定滑 面求β 在可能的滑面中 寻找一个所谓临界滑面 其β为最小值βm 注意蒙特卡洛法无法考虑参数相关 也不能进行自动寻找临界滑面的优化功能 4. 土压力计算 当 OP1(1)为 2 和 3 时进行土压力计算, 此时在表 12.1 的第 95 行输入以下参数 GWALL = 假定的土压力大小 假定值不影响最终成果 但不应与预计的最终成果偏离 过大 HMW= 土压力作用位置与墙底的垂直距离 EWALL= 土压力与水平线的夹角 当 OP1(1)=2 时进行同时满足力和力矩平衡条件的土压力计算 当 OP1(1)=3 时进行仅满足力平衡条件的土压力计算 即库仑土压力理论计算 此时计 算成果与输入的 HMW 值无关 12. 2. 5 应用最优化方法确定临界滑裂面 1. 圆弧滑裂面 程序提供的搜索圆弧形临界滑裂面的功能有 3 种 分别用表 12.2 中 23 行 LL2=0,1,2 来 实现 在一般情况下建议直接使用 LL2=2 即下述(3)的功能
390土质边坡德定分析一原理·方法程序 (1)枚举法(网格法)。此时令LL0=0,程序按输入的网格步长和网格数,计算相应不 同的圆心坐标和滑弧深度的圆弧的安全系数 此法的基本思想是不断地改变圆心位置x、y和滑弧深度D的数值,逐一比较相应的 安全系数,最终找到最小的安全系数。在具体操作中,先固定一个D,然后在圆心可能的 位置中布置一个网格。如图124所示,设网格的中心坐标为CCX,CCY。相应某一固定的 滑弧深度CCD表122中33行),在左右方向,各布置了NGRⅨX格,每格宽度为BCX,在 上下方向各布置了 NGRIY格,每格宽度为BCY,则共计有(2×NGRX+1)×(2× NGRIY+1 个网格点。分别以该网格点为圆心,以D3为滑弧深度计算相应安全系数,找出最小的安全 系数。然后改变一个D3值,重复相同的步骤。在这一过程中,有可能出现圆弧和边坡不相 交的情况,程序将自动抛弃该圆弧。同样Ds也是一个起算的中心值,在其上、下各布置 NGRID 层,每层间距为BCD。这样总计计算(2×NGRX+1)×(2×NGRY+1)×(2 K NGRID+1)个圆弧。 采用此法时,LL2=0。此法的一个特例是NGRX=NGRY= NGRID=0,此时计算一个指定 圆弧。 (2)二个自由度的搜索。应用单形法,先固定Ds不动,用单形法搜索与某一滑弧深度 D相应的使安全系数F最小的滑弧的圆心x、y。不断地改变Ds可获得一系列最小的安全 系数F 本功能以数据表LL2=1来实现。 (3)三个自由度的搜索。使用此功能,程序将一次计算出最小安全系数,不再对不同滑 弧深度逐一计算 本功能以LL2=2来实现。在使用单形法自动搜索最小安全系数和临界滑裂面时,数据 输入格式和枚举法完全一样,不同之处在于枚举法LL2=1,单形法LL2=2。输入NGRX NGRIY NGRID和CCX,CCY,CCD这些参数的作用是为最优化方法提供一个初始估计的滑 弧和搜索的范围,并不影响计算成果,也可设这些值为零。详见例题EX7 2.任意形状滑裂面 (1)任意形状滑裂面的几何参数。在表122第41行输入NS1,NOPI,IRE,即滑裂面 上控制点的总数,优化方法和迭代次数 输入KQ2(),Xl(I),Y1(①),NL(I),AL(①,分别为控制点的条块界面编号,x,y坐标值, 自由度,移动方向 A1和A2的控制点的条块界面编号分别为KQ2(1)和KQ2(2),假如KQ2(1)和KQ2(2)分 别为为1和5,那么,程序将自动在A1和A2之间内插编号为2、3、4的三条条块界面。依 次类推 自由度和移动方向分别在NL()和AL(1)中存放。 应用单形法,每个离散点都将从x,y两个方向逼近临界滑裂面相应的控制点B1,B2,… Bn此时我们称该点的自由度为2,对m个离散点,总的自由度为2m 有时,需要规定某个离散点固定不动,此时自由度为0,或规定沿某一设定方向移动 此时,自由度为1 如图126所示例,A4点如沿软弱面移动,则需规定该点自由度为1,并将移动方向相
390 土质边坡稳定分析 原理 ⋅ 方法 ⋅ 程序 (1) 枚举法 网格法 此时令 LL0=0 程序按输入的网格步长和网格数 计算相应不 同的圆心坐标和滑弧深度的圆弧的安全系数 此法的基本思想是不断地改变圆心位置 xc yc 和滑弧深度 Ds 的数值 逐一比较相应的 安全系数 最终找到最小的安全系数 在具体操作中 先固定一个 Ds 然后在圆心可能的 位置中布置一个网格 如图 12.4 所示 设网格的中心坐标为 CCX CCY 相应某一固定的 滑弧深度 CCD(表 12.2 中 33 行) 在左右方向 各布置了 NGRIX 格 每格宽度为 BCX 在 上下方向各布置了 NGRIY 格 每格宽度为 BCY 则共计有(2× NGRIX +1) × (2× NGRIY+1) 个网格点 分别以该网格点为圆心 以 DS 为滑弧深度计算相应安全系数 找出最小的安全 系数 然后改变一个 DS 值 重复相同的步骤 在这一过程中 有可能出现圆弧和边坡不相 交的情况 程序将自动抛弃该圆弧 同样DS也是一个起算的中心值 在其上 下各布置NGRID 层 每层间距为 BCD 这样总计计算(2×NGRIX+1)×(2× NGRIY+1)×(2× NGRID+1)个圆弧 采用此法时 LL2=0 此法的一个特例是 NGRIX=NGRIY=NGRID=0 此时计算一个指定 圆弧 (2) 二个自由度的搜索 应用单形法 先固定 DS 不动 用单形法搜索与某一滑弧深度 DS相应的使安全系数 F 最小的滑弧的圆心 xc yc 不断地改变 DS可获得一系列最小的安全 系数 F 本功能以数据表 LL2=1 来实现 (3) 三个自由度的搜索 使用此功能 程序将一次计算出最小安全系数 不再对不同滑 弧深度逐一计算 本功能以 LL2=2 来实现 在使用单形法自动搜索最小安全系数和临界滑裂面时 数据 输入格式和枚举法完全一样 不同之处在于枚举法 LL2=1 单形法 LL2=2 输入 NGRIX NGRIY, NGRID 和 CCX CCY, CCD 这些参数的作用是为最优化方法提供一个初始估计的滑 弧和搜索的范围 并不影响计算成果 也可设这些值为零 详见例题 EX7 2. 任意形状滑裂面 (1) 任意形状滑裂面的几何参数 在表 12.2 第 41 行输入 NS1 NOPT IRE 即滑裂面 上控制点的总数 优化方法和迭代次数 输入 KQ2(I) X1(I) Y1(I) NL(I) AL(I) 分别为控制点的条块界面编号 x, y 坐标值 自由度 移动方向 A1和 A2的控制点的条块界面编号分别为 KQ2(1)和 KQ2(2) 假如 KQ2(1)和 KQ2(2) 分 别为为 1 和 5 那么 程序将自动在 A1和 A2之间内插编号为 2 3 4 的三条条块界面 依 次类推 自由度和移动方向分别在 NL(I)和 AL(I)中存放 应用单形法 每个离散点都将从 x y 两个方向逼近临界滑裂面相应的控制点 B1 B2 Bm 此时我们称该点的自由度为 2 对 m 个离散点 总的自由度为 2m 有时 需要规定某个离散点固定不动 此时自由度为 0 或规定沿某一设定方向移动 此时 自由度为 1 如图 12.6 所示例 A4 点如沿软弱面移动 则需规定该点自由度为 1 并将移动方向相
第12章程序使用说明391 对于正x轴的夹角(规定从正x轴转向正y轴为正)作为一个参数输入程序。 在优化过程中,左、右两个端点可能移至边坡外或内,如图126中的A1和Am移至B 和Bm,程序将自动找到B1B2和BmBm1与边坡的交点B1和Bn,并对滑裂面的信息作相应 的调整。滑裂面上的NS1个控制点将其分成NS-1段。自上交点至下交点编号为1,2,,NS-1。 程序要求输入此NS-1段中形状为直线的段的编号,依次存放在数组LOO①中,直线段总 数为LNO,表122第70~72行。如果LNO=0,则滑裂面上无直线段,为光滑曲线。当LNO 为一大于NS-1的数时,则程序默认为全部为直线段。这两种情况都无须填写LOO①) (2)关于软弱夹层的输入。如前所述,我们将由NS1个点组成的滑裂面视为NS-1段, 当OPl(4)=1时,在输完KQ2()、X1(D)、Y1()、N(D)、AL(D)后,紧随于后输入软弱夹层 信息WS(①,I=1,…,NS1-1。当WS()=0时,滑裂面上第I段不是软弱夹层;当WS()不为 0时,wS(I)为该段滑裂面所属软弱夹层的编号。详见例题EX25和EX29在表122中55-65 行读入和处理这些信息。 3.应用随机搜索确定安全系数的整体极值 在计算任意形状滑裂面时,有时可能因初始滑裂面离整体极值相应的滑裂面较远,因而 无法找到。此时,如第4章介绍,可使用随机搜索的方法。先在临界滑裂面可能的位置附近 划定一个随机搜索区域,程序在给定区域自动生成100~500个滑裂面,从中找到相应最小安 全系数的那个滑裂面,作为最优化计算的初值。 如图1213所示,以初始滑裂面为中轴线,对滑裂面的每个控制点,在其移动方向,设 定的一个宽度d,称半带宽,这样即可形成一个随机搜索区域。 子程序READ4(表122)中67~68行读入相应的数据。 ⅨX=为随机搜索种子,可输入一任意整型数 KK,ZP(D)=滑裂面上控制点号和从该点定义的半带宽,KK=1,2,,NS1。 半带宽是按中轴线为轴线沿移动方向AL()定义的。用户应保证其输入的KK,ZP()确 实能形成如图12.13那样的具有物理意义的随机搜索区域,否则随机搜索并不能帮忙 随机搜索仅用于任意形状滑裂面。例题EX22介绍这一功能 随机生成的滑裂面 初始滑裂面中轴线 图1213随机搜素示意图
第 12 章 程序使用说明 391 对于正 x 轴的夹角 规定从正 x 轴转向正 y 轴为正 作为一个参数输入程序 在优化过程中 左 右两个端点可能移至边坡外或内 如图 12.6 中的 A1和 Am移至 B'1 和 B'm 程序将自动找到 B'1B2和 B'mBm-1与边坡的交点 B1和 Bm 并对滑裂面的信息作相应 的调整 滑裂面上的NS1个控制点将其分成NS−1段 自上交点至下交点编号为1,2,…, NS−1 程序要求输入此 NS−1 段中形状为直线的段的编号 依次存放在数组 LOO(I)中 直线段总 数为 LNO 表 12.2 第 70~72 行 如果 LNO=0 则滑裂面上无直线段 为光滑曲线 当 LNO 为一大于 NS1−1 的数时 则程序默认为全部为直线段 这两种情况都无须填写 LOO(I) (2) 关于软弱夹层的输入 如前所述 我们将由 NS1 个点组成的滑裂面视为 NS−1 段 当 OP1(4) = 1 时 在输完 KQ2(I) X1(I) Y1(I) NL(I) AL(I)后 紧随于后输入软弱夹层 信息 WS(I) I=1, …, NS1−1 当 WS(I)=0 时 滑裂面上第 I 段不是软弱夹层 当 WS(I)不为 0 时 WS(I)为该段滑裂面所属软弱夹层的编号 详见例题 EX25 和 EX29 在表 12.2 中 55~65 行读入和处理这些信息 3. 应用随机搜索确定安全系数的整体极值 在计算任意形状滑裂面时 有时可能因初始滑裂面离整体极值相应的滑裂面较远 因而 无法找到 此时 如第 4 章介绍 可使用随机搜索的方法 先在临界滑裂面可能的位置附近 划定一个随机搜索区域 程序在给定区域自动生成 100∼500 个滑裂面 从中找到相应最小安 全系数的那个滑裂面 作为最优化计算的初值 如图 12.13 所示 以初始滑裂面为中轴线 对滑裂面的每个控制点 在其移动方向 设 定的一个宽度 di 称半带宽 这样即可形成一个随机搜索区域 子程序 READ4 表 12.2 中 67∼68 行读入相应的数据 IX= 为随机搜索种子 可输入一任意整型数 KK ZP(I)= 滑裂面上控制点号和从该点定义的半带宽 KK=1,2,…,NS1 半带宽是按中轴线为轴线沿移动方向 AL(I)定义的 用户应保证其输入的 KK ZP(I)确 实能形成如图 12.13 那样的具有物理意义的随机搜索区域 否则随机搜索并不能帮忙 随机搜索仅用于任意形状滑裂面 例题 EX22 介绍这一功能 图 12. 13 随机搜索示意图