348土质边坡德定分析一原理·方法·程序 LPHS=IC(IN1, 3) IF(IPHEQ.O)GO TO 7 DO 73511=1.IPH MI=ICPH(Il, 1) M2=ICPH(1, 2) IF(X(K).LTXN(MD)GO TO 73: IF(X(K).GE XN(M2))GO TO 735 L(K=YN(MI)+X(K)XN(MI*(YNM2)YN(MI)(XN(M2)XN(MI) LPHS=ICPH(1, 3 73 CONTINUE CONTINUE PP=0 FT=0 l=0 DO 20711=LINI MI=IC(l, 1) M2=lC(,2) M3=I1 IF(X(K).GEXN(M2)GO TO 207 SH2=YNMl+(X(K)XN(MI))*(YN(M2YN(MIXN(M2}XNMI) 152 FORMAT(1X,**15,2F15.6) DSH2=ABS((Y(K)-SH2)) DSHI=SH2-Y(K) ASSURE THAT THE STRENGTH PARAMETERS OF THE UNDERLYING SOIL IF(DSH2. GT..00005AND DSHIGT.0.001GO TO 207 SH(KIFSH2 M4(K1FM3 207 CONTINUE IF(KI. NE IGO TO 707 GO TO 708 707LX=M4(1) HN=SH(I) DO45L1=2.K1 IF(SH(LI).GT HN)GO TO 45 H(L1) HN=SHLl) 45 CONTINUE 708Ml=C(LX,1) M2=IC(LX, 2) DELXXN(M2)-XN(MD IF(ABS(DELX)GT0.001)BBK=YN(M2)-YN(MI)(XN(M2)XN(MI X3(K)=x2(K) Y3(KFYNMI-+(X3(K)-XN(MD))BBK 139 FORMAT(IX, KI, LX, M2, MI, XM2, XMI, 415, 2F10.3) LINES(KFLX !96 12.29 47 DO 46 LI=LKI LSB=M4(LI M4(LIFIC(LSB, 3) CONTINU IF(WL(K). GT.Y(K)GO TO 734 DO 530 ISP=l K SDIF=ABS(WL(K)-SH(ISP)) IF(SDIF.LT0.00001GO TO 734348 土质边坡稳定分析 原理 ⋅ 方法 ⋅ 程序 LPHS=IC(IN1,3) IF(IPH.EQ.0)GO TO 739 DO 735 I1=1,IPH M1=ICPH(I1,1) M2=ICPH(I1,2) IF(X(K).LT.XN(M1))GO TO 735 IF(X(K).GE.XN(M2))GO TO 735 IF(ABS(XN(M2)-XN(M1)).LT.0.001)GOTO 735 WL(K)=YN(M1)+(X(K)-XN(M1))*(YN(M2)-YN(M1))/(XN(M2)-XN(M1)) LPHS=ICPH(I1,3) 735 CONTINUE 739 CONTINUE PP=0. CT=0. FT=0. K1=0 DO 207 I1=1,IN1 M1=IC(I1,1) M2=IC(I1,2) M3=I1 IF(X(K).LT.XN(M1))GO TO 207 IF(X(K).GE.XN(M2))GO TO 207 SH2=YN(M1)+(X(K)-XN(M1))*(YN(M2)-YN(M1))/(XN(M2)-XN(M1)) 152 FORMAT(1X,'**',I5,2F15.6) DSH2=ABS((Y(K)-SH2)) DSH1=SH2-Y(K) C ASSURE THAT THE STRENGTH PARAMETERS OF THE UNDERLYING SOIL IF(DSH2.GT.0.00005.AND.DSH1.GT.0.001)GO TO 207 K1=K1+1 SH(K1)=SH2 M4(K1)=M3 207 CONTINUE IF(K1.NE.1)GO TO 707 LX=M4(1) GO TO 708 707 LX=M4(1) HN=SH(1) DO 45 L1=2,K1 IF(SH(L1).GT.HN)GO TO 45 LX=M4(L1) HN=SH(L1) 45 CONTINUE 708 M1=IC(LX,1) M2=IC(LX,2) DELX=XN(M2)-XN(M1) IF(ABS(DELX).GT.0.001)BBK=(YN(M2)-YN(M1))/(XN(M2)-XN(M1)) X3(K)=X2(K) Y3(K)=YN(M1)+(X3(K)-XN(M1))*BBK 139 FORMAT(1X,'K1,LX,M2,M1,XM2,XM1',4I5,2F10.3) LINES(K)=LX !96.12.29 47 DO 46 L1=1,K1 LSB=M4(L1) M4(L1)=IC(LSB,3) 46 CONTINUE IF(WL(K).GT.Y(K))GO TO 734 DO 530 ISP=1,K1 SDIF=ABS(WL(K)-SH(ISP)) IF(SDIF.LT.0.00001)GO TO 734