- 积分
- 0
- 注册时间
- 2016-7-25
- 仿真币
-
- 最后登录
- 1970-1-1
|
abaqus热分析新手加子程序新手。将abaqus帮助文档的umatht例子放入简单的5*5*5的立方体中,在property中输入了Conductivity、Density 和Specific heat计算没有产生内部热源。请问还需要设置什么?子程序中的PROPS(1)、PROPS(2)应该在哪里设置?在property里选择user material,类别选择Thermal后,在1中输入Conductivity的数值,2中输入Specific heat的数值,提交后显示ERROR啊>_<含泪求解答
SUBROUTINE UMATHT(U,DUDT,DUDG,FLUX,DFDT,DFDG,
1 STATEV,TEMP,DTEMP,DTEMDX,TIME,DTIME,PREDEF,DPRED,
2 CMNAME,NTGRD,NSTATV,PROPS,NPROPS,COORDS,PNEWDT,
3 NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME
DIMENSION DUDG(NTGRD),FLUX(NTGRD),DFDT(NTGRD),
1 DFDG(NTGRD,NTGRD),STATEV(NSTATV),DTEMDX(NTGRD),
2 TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3)
C
COND=PROPS(1)
SPEC=PROPS(2)
C
DUDT=SPEC
DU=DUDT*DTEMP
U=U+DU-SPEC*DTEM
C
DO I=1,NTGRD
FLUX(I)=-COND*DTEMDX(I)
DFDG(I,I)=-COND
END DO
RETURN
END
|
|