392土质边坡德定分析一原理·方法·程序 12.3数据文件的格式 12.3.1主要数据 STAB读入数据文件的源程序如表12.1所示 表12.1主要数据输入格式 说明 1 SUBROUTINE READI 2 COMMONAl/N,IN,PF(20,PC(20)PDW(20),PDS(2用户需确认没有使用越界的数组 0)PH,ICPH(40,3) 3 COMMON/A2/NN, NJ(20), NHY, DHRP, WL(80), BDA 80),CA(8O)FAA(80) 4#,LRU(20+B59 5 COMMON/A3/RWRUL UWL 6 COMMON/A4/LSL, LNUM(60), IUPWJUPW 7 COMMON/A6/X(80), Y(80), XN(99), YN(99),X2(80),Y 8 COMMON/A7/W(80), RU(80), C(80), F(80) 9 COMMON/A9/DIR, FQUH(80), FQUV(80), FDIS(80), I l0 COMMON/AlO/PPX(15, 15), PY(15, 15), P(15, 15),LOXO 15),LOY 11 COMMON/ l/FA, FB, FK, FX(80), FXO(80),ALF(80), NS, GF(80), MQP, KXYX, ASP 12 COMMON/A1S/UWLI, NCOUN, MLA(80), PFI(20), P CI(20,MID(20) 13 COMMON/A16/IPH2, ICPH2(40, 3), IPHl, ICPHI(40, 3) 14 COMMON/A9A/LQH, CYE(1O), COE(10), VOE(10) 15 COMMON/CRU/RUS(IO),LINES(8O), YSF(80) 16 COMMON/RELIAI/NPROB, NLP, ITM, IPBD, MEEDO 17 COMMON/SUR/NSUR, NNP, NP(IO), PN(10, 2), PT(10 18 INTEGER 4OPTION(6),OP1(6) 19 COMMON/OPP/OPTION OP1 20 COMMON/WALL/TWALL GWALL HMWEWALL 21 COMMON/SWEAK/TWEAK, KQ3(50),LWK, IC3(40, 3 22 READ(5, (A64)HEAD 标题 23 READ(5, (OPTION), I=1, 6) 读入6个功能控制码 24 IF(OPTION(3).EQ.-2)READ(, )HITE 拉力缝内充水高度由用户指定 25 IF(OPTION(6).EQ. ITHEN 在 OPTION(6)为时读入6个扩充功能 26 READ(, OPI I=1, 28 !OPI(DWALL OPl(1)不为零时进行土压力计算 29 !OPI(3)FO(X)ANDF(X)SPECIFIEDBYUSER OPl(3)=0时程序默认使用 SPENCER
392 土质边坡稳定分析 原理 ⋅ 方法 ⋅ 程序 12. 3 数据文件的格式 12. 3. 1 主要数据 STAB 读入数据文件的源程序如表 12.1 所示 表 12. 1 主要数据输入格式 编号 源程序 说明 1 SUBROUTINE READ1 2 COMMON/A1/IN,IN1,PF(20),PC(20),PDW(20),PDS(2 0),IPH,ICPH(40,3) 用户需确认没有使用越界的数组 3 COMMON/A2/NN,NJ(20),NHY,DHRP,WL(80),BDA( 80),CA(80),FAA(80) 4 #,LRU(20)+B59 5 COMMON/A3/RW,RU1,UWL 6 COMMON/A4/LSL,LNUM(60),IUPW,JUPW 7 COMMON/A6/X(80),Y(80),XN(99),YN(99),X2(80),Y 2(80),IC(80,3),N 8 COMMON/A7/W(80),RU(80),C(80),F(80) 9 COMMON/A9/IDIR,FQUH(80),FQUV(80),FDIS(80),I QUA,HBAS,HBO 10 COMMON/A10/PPX(15,15),PY(15,15),P(15,15),LOX( 15),LOY 11 COMMON/A11/FA,FB,FK,FX(80),FXO(80),ALF(80), NS,GF(80),MQP,KXYX,ASP 12 COMMON/A15/UWL1,NCOUN,MLA(80),PF1(20),P C1(20),MID(20) 13 COMMON/A16/IPH2,ICPH2(40,3),IPH1,ICPH1(40,3) 14 COMMON/A9A/LQH,CYE(10),COE(10),VOE(10) 15 COMMON/CRU/RUS(10),LINES(80),YSF(80) 16 COMMON/RELIA1/NPROB,NLP,ITM,IPBD,MEED( 40) 17 COMMON/SUR/NSUR,NNP,NP(10),PN(10,2),PT(10, 2) 18 INTEGER*4OPTION(6),OP1(6) 19 COMMON/OPP/OPTION,OP1 20 COMMON/WALL/IWALL,GWALL,HMW,EWALL 21 COMMON/SWEAK/IWEAK,KQ3(50),LWK,IC3(40,3 ) 22 READ(5,'(A64)')HEAD 标题 23 READ(5,*)(OPTION(I),I=1,6) 读入6个功能控制码 24 IF(OPTION(3).EQ.-2)READ(5,*)HITE 拉力缝内充水高度由用户指定 25 IF(OPTION(6).EQ.1)THEN 在OPTION(6)为1时读入6个扩充功能 码 26 READ(5,*)(OP1(I),I=1,6) 27 ENDIF 28 !OP1(1)WALL OP1(1)不为零时进行土压力计算 29 !OP1(3)F0(X)ANDF(X)SPECIFIEDBYUSER OP1(3) = 0时程序默认使用SPENCER 法
第12拿程序使用说明393 续表 编号源程序 说明 OPl(3)不为0时,用户输入的关信息进 行BN通用条分法计算 30 !OPI(4)WEAKSEAM OP(4)不为0时,执行软弱夹层功能 31 !OPI(5)InputxaO 32oPl(6外边坡线由用户输入 OP1(6)不为0时,外边坡线由用户输入 33 READ(S, MTE 择有效应力法或总应力法 4READ5,*儿 OWL, UWLL, ASP 坡外水位、总应力法骤降水位、滑楔 法控制码 35 NSS=0 6 IF(ABSOUWLI-UWL)GT.0.OOINSS=1 总应力法坡外水位不同于骤降水位 使用库水位骤降总应力法,此时NSS= 37 MTE=0 38 MRD=0 39 READ(, )IQUA, HBAS, HBO, IDIR 读入地震加速度信息 40 IF(IQUA. GE.O)GOTO841 当IQUA小于零时,地震系数由用户输 入 41 IF(QUA NE-2 IQUA=-2时同时输入垂直和水平地 READ(S, ) LQH,(CYE(), COED), I=1, LQH) 42 IF(IQUA EQ READ( LQH, (CYE, COE(D,VOE I=1, LQH) 否则仅输入水平地震力 43 DO8431=1, LQH 44 IF(IQUA NE-2)VOE=0 45 843CONTINUE 46 READ(5, RW 水容重 47 READ(S, ")NN, L, XND, YN I=1,NN) 节点信息 8READ(S,*)N(LPF(D,PC①,PFl(D,PCI(DPDW(①,土层信息 PDS(), 49 %LRU I=1, IN) 50 IF(MTE. GE.OGOTO279 51 CALL READ2 输入地震总应力法动强度 52279NHY=0 53DO342I=1,N 处理孔隙水压力信息 55 IF(LRUO)GE.OJGOTO323 当LRU(①)小于零时读入该层土的孔隙 水压力系数 56 READ(S, RUSO 孔隙水压力系数 57NHY=-1 58 GOTO342 59 323IF(LRU).EQO)GOTO342 61 GOTO343 62 342CONTINUE
第 12 章 程序使用说明 393 续 表 编号 源程序 说明 OP1(3)不为0时 用户输入的关信息进 行BN通用条分法计算 30 !OP1(4)WEAKSEAM OP1(4)不为0时 执行软弱夹层功能 31 !OP1(5)Inputxa() 32 !OP1(6)外边坡线由用户输入 OP1(6)不为0时 外边坡线由用户输入 33 READ(5,*)MTE 选择有效应力法或总应力法 34 READ(5,*)UWL,UWL1,ASP 坡外水位 总应力法骤降水位 滑楔 法控制码 35 NSS=0 36 IF(ABS(UWL1-UWL).GT.0.001)NSS=1 总应力法坡外水位不同于骤降水位 使用库水位骤降总应力法,此时NSS = 1 37 MTE=0 38 MRD=0 39 READ(5,*)IQUA,HBAS,HBO,IDIR 读入地震加速度信息 40 IF(IQUA.GE.0)GOTO841 当IQUA小于零时 地震系数由用户输 入 41 IF(IQUA.NE.-2) READ(5,*)LQH,(CYE(I),COE(I),I=1,LQH) IQUA = −2时同时输入垂直和水平地 震力 42 IF(IQUA.EQ.-2) READ(5,*)LQH,(CYE(I),COE(I),VOE(I),I=1,LQH) 否则仅输入水平地震力 43 DO843I=1,LQH 44 IF(IQUA.NE.-2)VOE(I)=0. 45 843CONTINUE 46 READ(5,*)RW 水容重 47 READ(5,*)NN,(L,XN(I),YN(I),I=1,NN) 节点信息 48 READ(5,*)IN,(L,PF(I),PC(I),PF1(I),PC1(I),PDW(I), PDS(I), 土层信息 49 %LRU(I),I=1,IN) 50 IF(MTE.GE.0)GOTO279 51 CALL READ2 输入地震总应力法动强度 52 279NHY=0 53 DO342I=1,IN 处理孔隙水压力信息 54 RUS(I)=0. 55 IF(LRU(I).GE.0)GOTO323 当LRU(I)小于零时读入该层土的孔隙 水压力系数 56 READ(5,*)RUS(I) 孔隙水压力系数 57 NHY=-1 58 GOTO342 59 323IF(LRU(I).EQ.0)GOTO342 60 NHY=1 61 GOTO343 62 342CONTINUE
394土质边坡德定分析一原理·方法·程序 续表 编号源程序 说明 63 343CONTUNUE 64 IF(NHY. LE.OGOTO2010 65 READ(S, LOY 输入孔压或不排水强度内插网 6 DO2011=LLOY 相应MTE=0和>0分别输入孔隙水压 67 READ(, LOX) 68 I=LOX 69 READ(5, )(PPX( ) PY(D),P(LJ=1, I1) 0 201ICONTINUE 71 IF(OPTION(4).EQ.O)GOTO1874 72 CALL READS 可靠度分析的物理力学参数 线段和浸润线总数 75 READ(5, (CPH(D),J=1,3),I=1, PH) 76 225IF(OPI(4).GT.O)THEN 读入软弱夹层线段信息 77 READ(, LWK 78READ(S,*)(C3(,J),J=1,3)=1,LWK) 80 IF(NSSEQO)GOTO431 库水位骤降情况读入降后浸氵 81READ(S,川PH1 82 READ(S, (CPHI(D),J=1,3), I=1, PHI) 83 43ICONTINUE 84 CONTINUE 85 IF(OPI(6).EQOTHEN! 读入外边坡线信息 87 CALL OUTIDEN(NNXN, YNINLIO,LSL,LNUM)调用确定外边坡线的程序 88 ELSE 90 READ(, LNUMGID), II=1, LSL) 外边坡线编号 91 ENDIF 92 IF(OPTION(5).GT.O)CALL READ3 I 第一种功能输入外荷载 93 IF(OPTION(5).LT.OCALL READ3 按第二种功能输入外荷载 94 IF(OPI(I).NEOTHEN 算主动土压力 95 READ(, GWALL, HMW,EWALL 假设的土压力值,土压力作用高度,角 度 97 CALL READ4(IDO) 读入滑裂面信息 98 RETU 99 END
394 土质边坡稳定分析 原理 ⋅ 方法 ⋅ 程序 续 表 编号 源程序 说明 63 343CONTUNUE 64 IF(NHY.LE.0)GOTO2010 65 READ(5,*)LOY 输入孔压或不排水强度内插网 66 DO2011I=1,LOY 相应MTE = 0和>0分别输入孔隙水压 和qcu 67 READ(5,*)LOX(I) 68 I1=LOX(I) 69 READ(5,*)(PPX(I,J),PY(I,J),P(I,J),J=1,I1) 70 2011CONTINUE 71 IF(OPTION(4).EQ.0)GOTO1874 72 CALL READ5 读入可靠度分析的物理力学参数 73 1874READ(5,*)IN1,IPH 线段和浸润线总数 74 IF(IPH.EQ.0)GOTO225 75 READ(5,*)((ICPH(I,J),J=1,3),I=1,IPH) 76 225IF(OP1(4).GT.0)THEN 读入软弱夹层线段信息 77 READ(5,*)LWK 78 READ(5,*)((IC3(I,J),J=1,3),I=1,LWK) 79 ENDIF 80 IF(NSS.EQ.0)GOTO431 库水位骤降情况读入降后浸润线 81 READ(5,*)IPH1 82 READ(5,*)((ICPH1(I,J),J=1,3),I=1,IPH1) 83 431CONTINUE 84 CONTINUE 85 IF(OP1(6).EQ.0)THEN! 读入外边坡线信息 86 LSL=0 87 CALL OUTIDEN(NN,XN,YN,IN1,IC,LSL,LNUM) 调用确定外边坡线的程序 88 ELSE 89 READ(5,*)LSL 外边坡线总数 90 READ(5,*)(LNUM(II),II=1,LSL) 外边坡线编号 91 ENDIF 92 IF(OPTION(5).GT.0)CALL READ3_1 按第一种功能输入外荷载 93 IF(OPTION(5).LT.0)CALL READ3 按第二种功能输入外荷载 94 IF(OP1(1).NE.0)THEN 计算主动土压力 95 READ(5,*)GWALL,HMW,EWALL 假设的土压力值,土压力作用高度,角 度 96 ENDIF 97 CALL READ4(IDO) 读入滑裂面信息 98 RETURN 99 END
第12章程序使用说明395 12.3.2滑裂面信息 表122为STAB读入的圆弧和任意形状滑裂面。 表12.2滑裂面信息数据输入格式 编号源程序 说明 2 COMMON/A4/LSL, LNUM(60), IUPWJUPW COMMON/AS/NLOW, NUPP, NDP, NLOW1, CX, CY, D COMMON/A6X(80,Y(80)XN(99),YN99),X2(80,Y 2(80),C(80,3),N 5 COMMON/A8/MYK KXSX, MTE, KL 6 COMMON/All/FA, FB, FK, FX(80), FXO(80),ALF(80), 7 COMMON/A12/DT, BET(80), DQ, DM, AB1, ALAM, D B, DF, NO, NOO(20), LMC 8 COMMON/A13/ID. IWRLIWR2IWR3 9 COMMON/A23NS1KQ2(50),X1(50)Y1(50) 10 COMMON/A24/LNO, LOO(60) 11 COMMON/BI/ZB, NL(50),AL(S0), XE(50), YE(50), ZS( 12 COMMON/B2/MPU, MPUA(60), MPL, MPLA(60) 13 COMMON/B4/NOPT, IRE,TL 14 COMMON/B5 V/NED/ISP/ISPE, IBOUN, FAB, FBB XF UN(20), FUN(20), IFUN 5 COMMON/RN/LRND, IX, CONR, MRS/TW/TWR4, WR 5,WR6 16 DIMENSION WS(20) 17 COMMON/SER/SERIAL ISIM 18 INTEGER 4OPTION(6),OPI(6) 19 COMMON/OPP/OPTION, OPI/SWEAK/TWEAK, KQ3 (50),LWK,IC3(40,3) 20 COMMON/WALL/TWALL, GWALL, HMW,EWALL ETA 21 COMMON/AI/N, INI, PF(20), PC(20), PDW(20), PDS(2 0)PH,ICPH(40,3) 22 IF(MYK EQ0OR MYK. LT.OGOTO22 3 READ(, LLO, LLL, LL2 读入圆弧滑裂面信息 .6 IF(LLO EQ.OGOTO102 27 READ(S, NGRID, NS, CCD, BCD 读入固定圆弧上、下交点时的信息 28 IF(LLO GT.O)READ(, )XLl, YLl, XL2, YL2 29 IF(LLO.LTOREAD(S, )XLL, YLL, LDN, DS 30 RETURN 32C 33 102READ(S, ")CCX, CCY, CCD, BCX, BCY, BCD 按网格布置圆心位置和滑弧深度 4 READ(, NGRIX, NGRIY, NGRID, NS
第 12 章 程序使用说明 395 12. 3. 2 滑裂面信息 表 12.2 为 STAB 读入的圆弧和任意形状滑裂面 表 12. 2 滑裂面信息数据输入格式 编号 源程序 说明 1 SUBROUTINE READ4(IDO) 2 COMMON/A4/LSL,LNUM(60),IUPW,JUPW 3 COMMON/A5/NLOW,NUPP,NDP,NLOW1,CX,CY,D S 4 COMMON/A6/X(80),Y(80),XN(99),YN(99),X2(80),Y 2(80),IC(80,3),N 5 COMMON/A8/MYK,KXSX,MTE,KL 6 COMMON/A11/FA,FB,FK,FX(80),FXO(80),ALF(80), NS,GF(80),MQP,KXYX,ASP 7 COMMON/A12/DT,BET(80),DQ,DM,AB1,ALAM,D B,DF,NO,NOO(20),LMC 8 COMMON/A13/IID,IWR1,IWR2,IWR3 9 COMMON/A23/NS1,KQ2(50),X1(50),Y1(50) 10 COMMON/A24/LNO,LOO(60) 11 COMMON/B1/ZB,NL(50),AL(50),XE(50),YE(50),ZS( 50),ZP(50) 12 COMMON/B2/MPU,MPUA(60),MPL,MPLA(60) 13 COMMON/B4/NOPT,IRE,TL 14 COMMON/B5V/NED/ISP/ISPE,IBOUN,FAB,FBB,XF UN(20),FUN(20),IFUN 15 COMMON/RN/LRND,IX,CONR,MRS/IW/IWR4,IWR 5,IWR6 16 DIMENSION WS(20) 17 COMMON/SER/SERIAL,ISIM 18 INTEGER*4OPTION(6),OP1(6) 19 COMMON/OPP/OPTION,OP1/SWEAK/IWEAK,KQ3 (50),LWK,IC3(40,3) 20 COMMON/WALL/IWALL,GWALL,HMW,EWALL, ETA 21 COMMON/A1/IN,IN1,PF(20),PC(20),PDW(20),PDS(2 0),IPH,ICPH(40,3) 22 IF(MYK.EQ.0.OR.MYK.LT.0)GOTO22 23 READ(5,*)LLO,LL1,LL2 读入圆弧滑裂面信息 26 IF(LLO.EQ.0)GOTO102 27 READ(5,*)NGRID,NS,CCD,BCD 读入固定圆弧上 下交点时的信息 28 IF(LLO.GT.0)READ(5,*)XL1,YL1,XL2,YL2 29 IF(LLO.LT.0)READ(5,*)XL1,YL1,LDN,DS 30 RETURN 31 C 32 C 33 102READ(5,*)CCX,CCY,CCD,BCX,BCY,BCD 按网格布置圆心位置和滑弧深度 34 READ(5,*)NGRIX,NGRIY,NGRID,NS
396土质边坡德定分析一原理·方法·程序 续表 编号源程序 说明 35 IF(OPTION(3).NEOTHEN 36 READ(, ) YTENSION 输入圆弧滑裂面的拉力缝深度 37 ENDIF 38 RETURN 40C 41 22 READ(, *)NSI, NOPT,IRE 滑裂面上控制点的总数,优化方 迭代次数 46READ(5,*)KQ2(,X1(D),Y1(,NLI),AL(1 控制点的条块界面编号,x、y坐标值, 自由度,移动方向 47 259 CONTINUE 51DOI=1,NS1-1 52Ws(I)=0 53 ENDDO 55 IF(OPI(4).GT.O)READ(S, CWSI=I, NSI-1) OI=1.NS1-1 57 IF(WS).EQ.OCYCLE 58 II=wS KQ3(1+1)=Ⅱ 61 NQL=wS( 621=IC(NLQ,2) 63Y1(=(x1(XN()*S1+YN(I1) 将代表软弱夹层的控制点的y坐标调 整到紧贴软弱夹层 64AL(I=ATAN(S)*180/3.14159 将控制点的移动方向调整到与软弱夹 层平行 65 ENDDO 66 IF(NOPT. LT. 10GOTO511 67 READ(5, ) IX 读入随机搜索参数 68 READ(S, )(KK, ZP(), I=1, NS1) 滑裂面上控制点号和从该点定义的半 带宽 69 51ICONTINUE 70 READ(S, LNO 71 IF(LNO. LT NS1OR LNO GT.0) 读入滑裂面上形状为直线的段的编号 72 #READ(, (LOO),I=1, LNO) 当LNO为零时,滑裂面为一光滑曲线, 没有直线段。当LNO为一大于NSl-1 的数时,则程序默认为全部为直线段 74 CONTINUE 这两种情况都无须填写LOO(I 5 CONTINUE
396 土质边坡稳定分析 原理 ⋅ 方法 ⋅ 程序 续 表 编号 源程序 说明 35 IF(OPTION(3).NE.0)THEN 36 READ(5,*)YTENSION 输入圆弧滑裂面的拉力缝深度 37 ENDIF 38 RETURN 39 C 40 C 41 22 READ(5,*)NS1,NOPT,IRE 滑裂面上控制点的总数 优化方法 迭代次数 44 C 45 DO 259 I=1,NS1 46 READ(5,*)KQ2(I),X1(I),Y1(I),NL(I),AL(I) 控制点的条块界面编号 x y坐标值 自由度 移动方向 47 259 CONTINUE 48 C 49 C 50 C 51 DO I=1,NS1-1 52 WS(I)=0 53 ENDDO 54 C 55 IF(OP1(4).GT.0) READ(5,*)(WS(I),I=1,NS1-1) 56 DO I=1,NS1-1 57 IF(WS(I).EQ.0)CYCLE 58 II=WS(I) 59 KQ3(I)=II 60 KQ3(I+1)=II 61 NQL=WS(I) 62 I2=IC(NLQ,2) 63 Y1(I)=(X1(I)-XN(I1))*S1+YN(I1) 将代表软弱夹层的控制点的 y 坐标调 整到紧贴软弱夹层 64 AL(I)=ATAN(S1)*180./3.14159 将控制点的移动方向调整到与软弱夹 层平行 65 ENDDO 66 IF(NOPT.LT.10)GOTO511 67 READ(5,*)IX 读入随机搜索参数 68 READ(5,*)(KK,ZP(I),I=1,NS1) 滑裂面上控制点号和从该点定义的半 带宽 69 511CONTINUE 70 READ(5,*)LNO 71 IF(LNO.LT.NS1.OR.LNO.GT.0) 72 #READ(5,*)(LOO(I),I=1,LNO) 73 C 74 CONTINUE 75 CONTINUE 读入滑裂面上形状为直线的段的编号 当LNO为零时 滑裂面为一光滑曲线 没有直线段 当LNO为一大于NS1−1 的数时 则程序默认为全部为直线段 这两种情况都无须填写LOO(I)