如何利用abaqus用户子程序和应用程序输出单元PE文件？

xyt_lee

求指导！！！:'(:'(
各位大佬们，帮忙看一下下面子程序为何没法输出单元塑性应变？
  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