恳请看下编写的一段子程序有啥问题UTEMP

806940059 · 发表于 2014-4-1 20:34:33

本帖最后由 806940059 于 2014-4-1 20:41 编辑

用fortran编写abaqus调用的UTEMP子程序，调入计算结果全部为零，不知道是啥原因啊！求指教：
   SUBROUTINE UTEMP(TEMP,NSECPT,KSTEP,KINC,TIME,NODE,COORDS)
   INCLUDE 'ABA_PARAM.INC'
   DIMENSION TEMP(NSECPT), TIME(2), COORDS(3)
   INTEGER NODE

   INTEGER I,J,TFLAG,TAGBEAM
   REAL*8 T,X,Y,Z,POINT
   REAL*8 A1(5),A_COEF1(6,5),A_COEF2(6,5),A_COEF3(6,5),A_COEF4(6,5)
   REAL*8 A_COEF5(6,5)
   REAL*8 TTEMP

      T = TIME(1)
      IF(T.LT.0.0)THEN
      T = 0.0
      ENDIF

      TTEMP=30.0+345.0*log10(8.0*T/60.0+1.0)

         SELECT CASE(NODE)
         CASE(1:5)
         TEMP(1) = 0.7*TTEMP
         TEMP(2) = 0.1*TTEMP
         TEMP(3) = 0.4*TTEMP
         TEMP(4) = 0.7*TTEMP
         TEMP(5) = 0.1*TTEMP
         CASE(6:11)
         TEMP(1) = 0.7*TTEMP
         TEMP(2) = 0.1*TTEMP
         TEMP(3) = 0.4*TTEMP
         TEMP(4) = 0.7*TTEMP
         TEMP(5) = 0.1*TTEMP
      END SELECT

         DO I=1, NSECPT
               TEMP(I)=I*TTEMP
         ENDDO

            WRITE (*,*) NSECPT

      DO I=1, NSECPT
               WRITE (*,*) TEMP(I)
      ENDDO

      RETURN
      END