- 积分
- 0
- 注册时间
- 2014-6-3
- 仿真币
-
- 最后登录
- 1970-1-1
|
求指导!!!:'(:'(
各位大佬们,帮忙看一下下面子程序为何没法输出单元塑性应变?
SUBROUTINE URDFIL(LSTOP,LOVRWRT,KSTEP,KINC,DTIME,TIME)
C
INCLUDE 'ABA_PARAM.INC'
C
INTEGER A,B,EARRAY,NARRAY
REAL C,D,E,F,G,H,SARRAY,DARRAY,MARRAY,TARRAY,CARRAY
CHARACTER D1,D2,D3,D4,D5,D6,D7
DIMENSION ARRAY(513),JRRAY(NPRECD,513),TIME(2),LRUNIT(2,1)
1 ,COORD(3),PE(3),NARRAY(200000,1),DARRAY(200000,3),EARRAY(1200000,2)
1 ,SARRAY(1200000,6),MARRAY(1200000,6),TARRAY(1200000,6)
1 ,CARRAY(200000,3)
EQUIVALENCE (ARRAY(1),JRRAY(1,1))
LOVRWRT=1
OPEN(UNIT=16,FILE='E:\ABAQUSTEMP\\nopes.inp')
C
DO 1000 K2=1,10
DO 100 K1=1,999999
C
CALL DBFILE(0,ARRAY,JRCD) !读入一条记录
IF(JRCD.NE.0)GO TO 110 !判断结果文件是否结束
KEY=JRRAY(1,2) !把RECORD KEY赋值于KEY
IF(KEY.EQ.22) THEN
KEL=JRRAY(1,3)
PE(1)=ARRAY(3)
125 FORMAT(I6,',',F20.14)
WRITE(16,125) KEL,ARRAY(3) !输出单元X方向塑性应变
END IF
100 CONTINUE
1000 CONTINUE
110 CONTINUE
CLOSE(16)
RETURN
END
|
|