REAL AB1
INTEGER N,NS,OPTION
COMMON/AS/ASP
COMMON/CIRCLE/ID_C,CX,CY,DS
C OPTION=0, SPENCER法
C OPTION=1, BISHOP法
C OPTION=2, 瑞典法
C OPTION=3, 工程师团法
C OPTION=4, 简化法
C ID_C=0, 任意形状滑裂面; ID_C不为零,圆弧滑裂面
REAL X2(80),Y2(80),W(80),X(80),Y(80),C(80),F(80),WL(80)
$,FXO(80),FX(80),BN(80)
REAL BDA(80),FQUH(80),FDIS(80),YSF(80),RU(80),ALF(80),RW
write(6,*)'输入数据文件名'
OPEN(5,FILE=' ',STATUS='UNKNOWN')
write(6,*)'输入计算成果文件名'
OPEN(6,FILE=' ',STATUS='UNKNOWN')
RW=9.8
AB1=0.0
WRITE(6,*)'边坡稳定分析'
READ(5,*)OPTION,ID_C
IF(ID_C.NE.0)READ(5,*)CX,CY,DS
READ(5,*) NS !土条总数
READ(5,*)(X2(I),I=1,NS)
READ(5,*)(Y2(I),I=1,NS)
N=NS+1
DO 302 I=1,N
READ(5,*)W(I),RU(I),C(I),F(I),BDA(I),WL(I),FQUH(I),FDIS(I),YSF(I)
302 CONTINUE
!读入两端点和土条底中点的数据,其中涉及力的物理量均以单宽土条计
!X2,Y2=坐标,W=土条重量,RU=孔压系数,C=凝聚力=,F=磨擦系数,BDA=土容重,
!WL=浸润线Y坐标,FQUH,FDIS,YSF=水平地震力的大小,Y坐标,垂直水平地震力的大小
IF(OPTION.EQ.3)THEN
READ(5,*)ASP
WRITE(6,*)'陆军工程师团法,平均坝坡=',ASP
ASP=ASP*3.14159/180.
ENDIF
WRITE(6,709)RW
709 FORMAT(1X,'水容重=',F10.3//)
WRITE(6,702)
702 FORMAT(T25,'DATA FOR ASSUMED SIDE FORCE FUNCTION')

下一页