轮胎滚动下的温度场求解
你好，想请问一下用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
