各位ansys大神,请问我做这个热耦合分析的时候用的是shell57单元,可以计算出来温度分布,但是到热耦合分析的时候怎么计算不出来结果了,请大神指点?另外我看shell57对应的热耦合结构单元为shell63为什么在ansys运行中 我查询到是变成了shell41了?
finish
/clear
/COM,Thermal
L1=765E-3
L2=1335E-3
L3=850E-3
L4=1000E-3
L5=2500E-3
L6=2500E-3
L7=1600E-3
/PREP7
et,1,SHELL57
r,1,8E-3
T0=28
Tair=20
Kss=173
Hair=20
mp,ex,1,3e4
!mp,prxy,1,0.2
!mp,dens,1,2600e-12
MP,KXX,1,KSS
T1=210E-3
T2=230E-3
T3=420E-3
T4=530E-3
T5=540E-3
c=500E-3
b=470E-3
h=1000E-3
h0=780E-3
s=440E-3
k, 1 , 0 , 0 , h
k, 2 , T1+T2 , , h
k, 3 , T1+T3 , , h
k, 4 , T1+L1 , , h
k, 5 , T1+L1 , , 129200/858*0.001
k, 6 , T1+T2 , , 400E-3
k, 7 , , , 400E-3
!k, 8 , T1+T3 , , 780E-3
k, 9 , T1+T3 , , 267200/858*0.001
k, 1001 , 0 , s , h
k, 1002 , T1+T2 , s , h
k, 1003 , T1+T3 , s , h
k, 1004 , T1+L1 , s , h
k, 1005 , T1+L1 , s , 129200/858*0.001
k, 1006 , T1+T2 , s , 400E-3
k, 1007 , , s , 400E-3
!k, 1008 , T1+T3 , s , 780E-3
k, 1009 , T1+T3 , s , 267200/858*0.001
A, 1001 , 1002 , 2 , 1
A, 1002 , 1003 , 3 , 2
A, 1003 , 1004 , 4 , 3
A, 1004 , 1005 , 5 , 4
A, 1007 , 1006 , 6 , 7
A, 1006 , 1009 , 9 , 6
A, 1009 , 1005 , 5 , 9
A, 2 , 3 , 9 , 6
AATT,1,1,1,,
NUMMRG,all
!合并
NUMCMP,all
!压缩
ALLSEL,ALL
ESIZE,,5,
MSHKEY,
AMESH,ALL
ALLSEL,ALL
finish
/SOLU
antype,static
Da,3,temp,T0
SFA,7,1,CONV,Hair,Tair
solve
finish
/POST1
/PREP7
ETCHG,TTS !将热单元转换为结构单元
r,1,8E-3
MP,EX,1,2.06E5 !单位是
MP,PRXY,1,0.3
MP,DENS,1,7.85E-9 !单位是
MP,ALPX,1,2.36E-5
FINISH
/SOLU
DA,1,ALL
DA,5,ALL
acel,,,98
LDREAD,TEMP,,,,,,RTH
TREF,20
!ALLSEL,ALL
SOLVE
FINISH |