ABAQUS岩土工程实例：deal.for文件操作

~未来可期，人生值得~

          DIMENSION NODE(1701),COORD(1701,3),DIS(10,1701,3),CDIS(1701,3)

DIMENSION FDIS(1701,3)

OPEN(1,FILE="U1.RPT")

OPEN(2,FILE="U2.RPT")

OPEN(3,FILE="U3.RPT")

OPEN(4,FILE="U4.RPT")

OPEN(5,FILE="U5.RPT")

OPEN(6,FILE="U6.RPT")

OPEN(7,FILE="U7.RPT")

OPEN(8,FILE="U8.RPT")

OPEN(9,FILE="U9.RPT")

OPEN(10,FILE="U10.RPT")

OPEN(11,FILE="OUTU.DAT")

DO 20 K=1,10

DO 10 I=1,1701

READ(k,*)NODE(I),COORD(I,1),COORD(I,2),

     1 COORD(I,3),DIS(K,I,1),DIS(K,I,2),DIS(K,I,3)

10 CONTINUE

20 CONTINUE

C 读入节点编号及相应的坐标、位移值

DO 30 I=1,1701

NFILL=ABS((COORD(I,3)-0.001))/10+1

C 确定当前节点属于哪一级加载

FEIH=10.0

C 当前填土高度，本例中各级填土均为10.0m

FEIZ=NFILL*10-COORD(I,3)

C 确定当前节点在分级填土中的埋深（从分级填土表面算起）

CDIS(I,1)=2*FEIZ/(FEIH+FEIZ)*DIS(NFILL,I,1)

CDIS(I,2)=2*FEIZ/(FEIH+FEIZ)*DIS(NFILL,I,2)

CDIS(I,3)=2*FEIZ/(FEIH+FEIZ)*DIS(NFILL,I,3)

C 进行位移修正

FDIS(I,1)=CDIS(I,1)+DIS(10,I,1)-DIS(NFILL,I,1)

FDIS(I,2)=CDIS(I,2)+DIS(10,I,2)-DIS(NFILL,I,2)

FDIS(I,3)=CDIS(I,3)+DIS(10,I,3)-DIS(NFILL,I,3)

C 累计位移

30 CONTINUE

DO 40 I=1,1701

WRITE(11,*)COORD(I,1),COORD(I,3),FDIS(I,1)*100,FDIS(I,3)*100

C 输出结果，将位移单位从m转化为cm

40 CONTINUE

END