342土质边坡德定分析一原理·方法程序 IF(ABS(XXL).LT0.0001THEN SLO=L0E-6 ELSE SLO=(YN(J2YN(J1)XXL!斜率 ENDIF IF(NEQ 1 OR SLO. LEO)GOTOll C第一条外边坡线为水平同时圆弧与其有两个交点不接受继续搜索 IF( R. LT D ORREQ. D)GOTO1无交点 NLOW=LNUM NUPPENLOW RETURN 0 IF(I NE.OGOTO14 NLOW=LNUM 11=11+I 11 CONTINUE IF(EQ 1 AND R.LTRIAND R GT R2)RETURN l=+1 IF(GT.LSLGOTO14 GOTO12 14 NUPP=LNUM 15 CONTINUE RETURN SUBROUTINE DD(X1, Y1, X2, Y2, XC,YC, D, IN) REAL XI, Y1, X2, Y2, XC,YC, D REAL XXL, AK, AKF X,YXT,YT NTEGER IN XXLEX2-XI IF(ABS(XXL). LT0.001)XXL=1.0E-6 IF(ABS((Y2-Y1)XXL). LT0.OOIRETURN AK=(Y2-YI/XXL AKF=l/AK X=(AK*Xl-Y1+AKF*XC+YC)(AK+AKF) Y=AK*XAK*X1+Y1!从圆心作垂直与线段的线获与线段的交点 X1OX-X2 YT=(YY1)*(Y-Y2) IF(XT.LT0OR YT. LT.OGOTO12 IN=11交点不在线段内 RETURN 12D=SQRT(X-XC)**2+(YYO)*2)!线段与圆心的垂直距离 RETURN END 3.例题 数据文件 CIRCLE DAT 00. 24.235-35.164-42.046 21 2 3342 土质边坡稳定分析 原理 ⋅ 方法 ⋅ 程序 IF(ABS(XXL).LT.0.0001)THEN SLO=1.0E-6 ELSE SLO=(YN(J2)-YN(J1))/XXL !斜率 ENDIF IF(IN.EQ.1.OR.SLO.LE.0)GOTO11 C 第一条外边坡线为水平,同时圆弧与其有两个交点,不接受,继续搜索 IF(R.LT.D.OR.R.EQ.D)GOTO11 !无交点 NLOW=LNUM(I) NUPP=NLOW RETURN 20 IF(I1.NE.0)GOTO14 NLOW=LNUM(I) I1=I1+1 11 CONTINUE IF(I.EQ.1.AND.R.LT.R1.AND.R.GT.R2)RETURN I=I+1 IF(I.GT.LSL)GOTO14 GOTO12 14 NUPP=LNUM(I) 15 CONTINUE RETURN END SUBROUTINE DD(X1,Y1,X2,Y2,XC,YC,D,IN) REAL X1,Y1,X2,Y2,XC,YC,D REAL XXL,AK,AKF,X,Y,XT,YT INTEGER IN D=9999. IN=0 XXL=X2-X1 IF(ABS(XXL).LT.0.001)XXL=1.0E-6 IF(ABS((Y2-Y1)/XXL).LT.0.001)RETURN AK=(Y2-Y1)/XXL AKF=1./AK X=(AK*X1-Y1+AKF*XC+YC)/(AK+AKF) Y=AK*X-AK*X1+Y1 !从圆心作垂直与线段的线,获与线段的交点 XT=(X-X1)*(X-X2) YT=(Y-Y1)*(Y-Y2) IF(XT.LT.0.OR.YT.LT.0)GOTO12 IN=1 !交点不在线段内 RETURN 12 D=SQRT((X-XC)**2+(Y-YC)**2) !线段与圆心的垂直距离 RETURN END 3. 例题 数据文件 CIRCLE.DAT 0 0. -24.235 -35.164 -42.046 21 9 1 2 3 2 3 1 3 4 1 2 5 3 5 6 2 6 7 2