我代码如下:
/units,si
pi=4*atan(1)
i=300
r=1.6e-3 !钨极半径
/PREP7
k,1,0,8.3e-3,,,,,
k,2,8.3e-3,8.3e-3,,,,,
k,3,8.3e-3,0,,,,,
k,4,0,0,,,,,
k,5,0,5.5e-3,,,,,
k,6,0.84e-3,5.5e-3,,,,,
k,7,1.6e-3,8.3e-3,,,,,
k,8,0.84e-3,0,,,,,
ET,1,PLANE53
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,1
KEYOPT,1,4,0
KEYOPT,1,5,0
KEYOPT,1,7,0
ET,2,0
MPTEMP,,,,,,,, !密度kg/m3
MPTEMP,1,273
MPTEMP,2,2500
MPTEMP,3,5000
MPTEMP,4,7500
MPTEMP,5,10000
MPTEMP,6,12500
MPTEMP,7,15000
MPTEMP,8,17500
MPTEMP,9,20000
MPDATA,DENS,1,,1.67
MPDATA,DENS,1,,0.2
MPDATA,DENS,1,,0.1
MPDATA,DENS,1,,0.07
MPDATA,DENS,1,,0.04
MPDATA,DENS,1,,0.03
MPDATA,DENS,1,,0.02
MPDATA,DENS,1,,0.015
MPDATA,DENS,1,,0.014
MPTEMP,,,,,,,, !热导率J/m*s*K
MPTEMP,1,273
MPTEMP,2,2500
MPTEMP,3,5000
MPTEMP,4,7500
MPTEMP,5,10000
MPTEMP,6,12500
MPTEMP,7,15000
MPTEMP,8,17500
MPTEMP,9,20000
MPDATA,KXX,1,,0
MPDATA,KXX,1,,0.1
MPDATA,KXX,1,,0.1
MPDATA,KXX,1,,0.3
MPDATA,KXX,1,,0.6
MPDATA,KXX,1,,1.6
MPDATA,KXX,1,,2.4
MPDATA,KXX,1,,2
MPDATA,KXX,1,,2.2
MPTEMP,,,,,,,, !电阻率 欧姆.m
MPTEMP,1,273
MPTEMP,2,2500
MPTEMP,3,5000
MPTEMP,4,7500
MPTEMP,5,10000
MPTEMP,6,12500
MPTEMP,7,15000
MPTEMP,8,17500
MPTEMP,9,20000
MPDATA,RSVX,1,,1000
MPDATA,RSVX,1,,100
MPDATA,RSVX,1,,1
MPDATA,RSVX,1,,1.3e-3
MPDATA,RSVX,1,,3.3e-4
MPDATA,RSVX,1,,1.6e-4
MPDATA,RSVX,1,,1.2e-4
MPDATA,RSVX,1,,9.8e-5
MPDATA,RSVX,1,,8.5e-5
MPTEMP,,,,,,,, !相对磁导率
MPTEMP,1,0
MPDATA,MURX,1,,0.9
MPTEMP,,,,,,,, !比热J/kg.K
MPTEMP,1,273
MPTEMP,2,2500
MPTEMP,3,5000
MPTEMP,4,7500
MPTEMP,5,10000
MPTEMP,6,12500
MPTEMP,7,15000
MPTEMP,8,17500
MPTEMP,9,20000
MPDATA,c,1,,521
MPDATA,c,1,,521
MPDATA,c,1,,521
MPDATA,c,1,,600
MPDATA,c,1,,1250
MPDATA,c,1,,6000
MPDATA,c,1,,10000
MPDATA,c,1,,4000
MPDATA,c,1,,2500
MPTEMP,,,,,,,, !粘度kg/m.s
MPTEMP,1,273
MPTEMP,2,2500
MPTEMP,3,5000
MPTEMP,4,7500
MPTEMP,5,10000
MPTEMP,6,12500
MPTEMP,7,15000
MPTEMP,8,17500
MPTEMP,9,20000
MPDATA,visc,1,,2e-5
MPDATA,visc,1,,7e-5
MPDATA,visc,1,,1e-4
MPDATA,visc,1,,2.2e-4
MPDATA,visc,1,,2.7e-4
MPDATA,visc,1,,2e-4
MPDATA,visc,1,,6e-5
MPDATA,visc,1,,3e-5
MPDATA,visc,1,,2.6e-5
!钨极材料性能
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,19300
MPDATA,RSVX,2,,0.0532
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,2,,180
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,130
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,2,,1.26e-6
!建模及网格划分
LSTR, 1, 5
LSTR, 5, 4
LSTR, 3, 2
LSTR, 2, 7
LSTR, 7, 1
LSTR, 5, 6
LSTR, 6, 7
LSTR,6,8
LSTR, 4, 8
LSTR, 8, 3
al,1,6,7,5
al,6,2,9,8
al,7,8,10,3,4
aglue,all
numcmp,area
lesize,6,,,20
lesize,1,,,10
lesize,2,,,60
lccat,4,3
mshape,0,2d
mshk,1
aatt,2,,1,,,
amesh,1
aatt,1,,1,,,
asel,s,area,,2,3,1
amesh,all
allsel,all
!加载及边界条件
ri=i/(pi*r**2)
BFA,1,JS, , ,ri,0
lsel,s,,,10
lsel,a,,,3
lsel,a,,,4
lsel,a,,,5
lsel,a,,,9
sfl,all,inf
lsel,s,,,5,7,1
lsel,a,,,1
lsel,a,,,9,10
BFL,all,TEMP,10000
lsel,s,,,3,4
bfl,all,temp,6000
/solu
ANTYPE,0
/SOL
MAGSOLV,0,3,0.001, ,25, |