将网上一段电热代码纠正后产生的温度场，欢迎使用，想学APDL编程的小伙伴，可以更改建模的坐标，建立任何想要的电热物理模型。加载方式KBC,1或KBC,0 、是否考虑辐射、电阻是否随温度变化都可以自己添加。
代码1
/FILNAME,Ex Eletric Thermal PolySi /PREP7
! Define problem parameters !
E3=155e9 P3=0.22
K3=150 !导热系数 A3=2.9e-6
R3=2.3e-3 !电阻率 C3=700 !热容 D3=2330 !密度 U=5
! Choose Element and define MP ET,1,SOLID69 KEYOPT,1,1,0 MP,RSVX,1,R3 MP,KXX,1,K3
! Finite element model
BLOCK,0,5E-2,0,5E-2,0,0.5E-2
BLOCK,0,4.5E-2,0.5E-2,4.5E-2,0,0.5E-2 /PNUM,VOLU,1 VPLOT VSBV,1,2 SAVE
LESIZE,ALL,0.5e-2 TYPE,1 MAT,1
SMRTsize,1 MSHAPE,1,3D MSHKEY,0 VMESH,all
/FACET,WIRE /VIEW,1,1,1,1 numcmp,all /pnum,volu,1
vplot
EPLOT SAVE
! Define Load and boundary conditions /PNUM,AREA,1 APLOT toffst,273 tunif,25 /SOL antype,4 allsel,all ASEL,S,,,7 nsla,s,1
DA,ALL,VOLT,0
*get,n_5,node,,num,min
f,n_5,amps,2.2*1!施加f集中电流载荷 allsel,all ASEL,S,,,8 nsla,s,1
DA,ALL,VOLT,U
ALLSEL ASEL,S,,,1 NSLA,S,1
SFA,ALL,1,CONV,12.5,20
ASEL,S,,,2 NSLA,S,1
SFA,ALL,1,CONV,12.5,20
ASEL,S,,,3 NSLA,S,1
SFA,ALL,1,CONV,12.5,20
ASEL,S,,,3 NSLA,S,1
SFA,ALL,1,CONV,12.5,20
ASEL,S,,,4 NSLA,S,1
SFA,ALL,1,CONV,12.5,20
ASEL,S,,,5 NSLA,S,1
SFA,ALL,1,CONV,12.5,20
ASEL,S,,,6 NSLA,S,1
SFA,ALL,1,CONV,12.5,20
ASEL,S,,,7 NSLA,S,1
SFA,ALL,1,CONV,12.5,20
ASEL,S,,,8 NSLA,S,1
SFA,ALL,1,CONV,12.5,20
ASEL,S,,,9 NSLA,S,1
SFA,ALL,1,CONV,12.5,20
ASEL,S,,,10 NSLA,S,1
SFA,ALL,1,CONV,12.5,20
ALLSEL
TRNOPT,FULL LUMPM,0 time,3535
NSUBST,25,25,10 OUTRES,ALL,ALL SOLVE /post1 Set,last
plnsol,temp,,0
温度场1
代码2
/FILNAME,Ex Eletric Thermal PolySi /PREP7
! Define problem parameters !
E3=155e9 P3=0.22
K3=150 !导热系数 A3=2.9e-6
R3=2.3e-3 !电阻率 C3=700 !热容 D3=2330 !密度 U=5
! Choose Element and define MP ET,1,SOLID69 KEYOPT,1,1,0 MP,RSVX,1,R3 MP,KXX,1,K3
! Finite element model
BLOCK,0,5E-2,0,5E-2,0,0.5E-2
BLOCK,0,4.5E-2,0.5E-2,4.5E-2,0,0.5E-2 /PNUM,VOLU,1 VPLOT VSBV,1,2 SAVE
LESIZE,ALL,0.5e-2 TYPE,1 MAT,1
SMRTsize,1 MSHAPE,1,3D MSHKEY,0 VMESH,all
/FACET,WIRE /VIEW,1,1,1,1 numcmp,all /pnum,volu,1 vplot
EPLOT
Eletric Thermal PolySi
