|
板凳
楼主 |
发表于 2019-12-24 13:07:29
|
只看该作者
感谢雪球。这就把代码弄来
SUBROUTINE URDFIL(LSTOP,LOVRWRT,KSTEP,KINC,DTIME,TIME)
C
INCLUDE'ABA_PARAM.INC'
C
DIMENSION ARRAY(513),JRRAY(NPRECD,513),TIME(2),LRUNIT(2,1)
1,COORD(3)
EQUIVALENCE(ARRAY(1),JRRAY(1,1))
CALL POSFIL(KSTEP,KINC,ARRAY,JRCD)
OPEN(UNIT=18,FILE='E:\tempabaqus\1224WX\F11.txt')
DO 1000 K2=1,10
DO 100 K1=1,9999999
C
CALL DBFILE(0,ARRAY,JRCD)
IF(JRCD.NE.0)GO TO 110
KEY=JRRAY(1,2)
C
IF(KEY.EQ.101)THEN
WRITE(18,160)ARRAY(6)
ELSE IF(KEY.EQ.102)THEN
WRITE(18,160)ARRAY(6)
160 FORMAT(F20.10,5X)
END IF
C
100 CONTINUE
1000 CONTINUE
110 CONTINUE
CLOSE(18)
RETURN
END
SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS
1,JLTYP,SNAME)
C
INCLUDE'ABA_PARAM.INC'
C
DIMENSION TIME(2),COORDS(3)
CHARACTER*80 SNAME
C
REAL UB,UT,VB,VT,U,V
OPEN(UNIT=18,FILE='E:\tempabaqus\1224WX\F11.txt')
OPEN(UNIT=101,FILE='E:\tempabaqus\1224WX\F12.txt')
READ(18,*)UB,UT,VB,VT
U=UT-UB
V=VT-VB
F=-100000
WRITE(101,180)UB,UT,U
180 FORMAT(F20.10,5X,F20.10,5X,F20.10,5X)
CLOSE(18)
RETURN
END |
|