|
本帖最后由 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
|
|