你好，想請問一下用abaqus子程序將生熱率賦值給模型上的各個單元，（計算出了各個單元的生熱率，每個單元生熱率都不一樣），我編寫了一下子程序，但是好像有點問題，能麻煩幫我看看嗎，謝謝

      SUBROUTINE UEXTERNALDB(LOP,LRESTART,TIME,DTIME,KSTEP,KINC)

C

      INCLUDE 'ABA_PARAM.INC'

C

      DIMENSION TIME(2)

C

      common/txt_data/ heat(2088),kele(2088),i

      if(lop.eq.0) then               !表示分析開始時正在調用用戶子程序

         open(106,file='D:\temp\A-usersubroutine\heat_jjx\heat.txt')      

         do i=1,2088

            read(106,*)a,b

            kele(i)=a

            heat(i)=b

         end do        

         close(106)  

         open(107,file='D:\temp\A-usersubroutine\heat_jjx\heat_UEX.txt')

         write(107,*)(heat(i),i=1,2088)             

         close(107)

      end if         

      RETURN

      END

C    

       SUBROUTINE SDVINI(STATEV,COORDS,NSTATV,NCRDS,NOEL,NPT,

     1 LAYER,KSPT)

C

      INCLUDE 'ABA_PARAM.INC'

C

      DIMENSION STATEV(NSTATV),COORDS(NCRDS)

C

      common/txt_data/ heat(2088),kele(2088),i

      do i=1,2088

          statev(i)=heat(i)

      end do

      do i=1,2088

          write(107,*)statev(i)

      end do

      return 

      end

C

      SUBROUTINE HETVAL(CMNAME,TEMP,TIME,DTIME,STATEV,FLUX,

     1 PREDEF,DPRED)

C

      INCLUDE 'ABA_PARAM.INC'

C

      CHARACTER*80 CMNAME

C

      DIMENSION TEMP(2),STATEV(*),PREDEF(*),TIME(2),FLUX(2),

     1 DPRED(*)

C  

      do i=1,2088

             FLUX(1)=statev(i)

      end do

      RETURN

      END