高手帮忙看看感应加热热分析的加载和约束

wangdou

很显然结果不对，应该是在热分析的时候加载以及约束的数值或条件不对.特别是在各条线上加对流系数以及空气场温度约束的时候,希望大虾们指教.
/prep7
cyl4,0.020,,0.001
agen,10,1,,,0.003
rectng,0,0.050,0.013,0.014
rectng,0.049,0.050,0.014,0.030
aadd,11,12
rectng,0,0.049,0.014,0.025
rectng,0,0.080,-0.030,0.070
aovlap,all
numcmp,area

! Electromagnetic model
*SET,freq,25000               ! frequency (Hz.)
*SET,pi,4*atan(1)               ! pi
*SET,cond,.392e7                ! maximum conductivity  
*SET,muzero,4e-7*pi             ! free-space permeability   
*SET,mur,200                    ! maximum relative permeability
*SET,skind,sqrt(1/(pi*freq*cond*muzero*mur))  ! skin depth  
*SET,ftime,300                   ! final time
*SET,tinc,.05                   ! time increment for harmonic analysis  
*SET,time,0                     ! initialize time   
*SET,delt,.01                   ! maximum delta time step  
emunit,mks
! Electromagnetic model
et,1,53,,,1                ! PLANE53, axisymmetric, AZ dof  
et,2,53,,,1
emunit,mks                 ! set magnetic units
mp,murx,1,1                ! air relative permeability  
mp,murx,3,1                ! coil relative permeability
mp,murx,4,1
mptemp
mptemp,1,25.5,160,291.5,477.6,635,698     ! temps for relative permeability
mptemp,7,709,720.3,742,761,1000
mpdata,murx,2,1,200,190,182,161,135,104   ! steel relative permeability
mpdata,murx,2,7,84,35,17,1,1
mptemp  
mptemp,1,0,125,250,375,500,625            ! temps for resistivity   
mptemp,7,750,875,1000   
mpdata,rsvx,2,1,.184e-6,.272e-6,.384e-6,.512e-6,.656e-6,.824e-6
mpdata,rsvx,2,7,1.032e-6,1.152e-6,1.2e-6  ! steel resistivity   
mptemp  
mptemp,1,0,730,930,1000                   ! temps for conductivity  
mpdata,kxx,2,1,60.64,29.5,28,28
mptemp                                    ! temps for enthalpy  
mptemp,1,0,27,127,327,527,727   
mptemp,7,765,765.001,927
mpdata,enth,2,1,0,91609056,453285756,1.2748e9,2.2519e9,3.3396e9
mpdata,enth,2,7,3.548547e9,3.548556e9,4.3520e9  
mp,emis,2,.68              ! emissivity

asel,s,,,13   
aatt,1,,2  
asel,s,,,11   
aatt,4,,2
asel,s,,,12   
aatt,2,,1  
asel,s,,,1,10,1   
aatt,3,1,2
asel,s,,,all

asel,s,,,1,10,1
smrt,1
mshape,0,2d
mshkey,0
amesh,all

asel,s,,,12
lsla,r
lesize,all,0.0002
mshape,0,2d
mshkey,0
amesh,12
asel,s,,,11
lsla,r
lesize,all,0.0002
mshape,0,2d
mshkey,0
amesh,11
asel,s,,,13
lsla,r
lesize,all,0.00002
mshape,0,2d
mshkey,0
amesh,13
nsel,s,loc,x,0
nsel,a,loc,x,0.08
nsel,a,loc,y,-0.03
nsel,a,loc,y,0.07
d,all,az,0
allsel,all

allsel,all
esel,s,mat,,3
bfe,all,js,,,,3.18e6
allsel,all

finish
/solu
antyp,harm                 ! Emag analysis options
harfrq,freq
solve
save
finish
!theraml
/prep7
etchg,mtt
mptemp  
mptemp,1,0,730,930,1000                   ! temps for conductivity  
mpdata,kxx,2,1,60.64,29.5,28,28
mptemp                                    ! temps for enthalpy  
mptemp,1,0,27,127,327,527,727   
mptemp,7,765,765.001,927
mpdata,enth,2,1,0,91609056,453285756,1.2748e9,2.2519e9,3.3396e9
mpdata,enth,2,7,3.548547e9,3.548556e9,4.3520e9  
mp,emis,2,.68              ! emissivity
mp,kxx,1,0.023              ! 空气水及铜的密度、导热系数和比热
mp,c,1,1.4
mp,dens,1,1.29
mp,kxx,3,270
mp,c,3,390
mp,dens,3,8900
mp,kxx,4,0.61
mp,c,4,4200
mp,dens,4,1000
finish
/solu
LDREAD,HGEN,,, ,0,'file','rmg',' '
allsel
asel,s,,,12
nsel,r,loc,y,0.013
sf,all,rdsf,0.68,1
allsel,all
asel,s,,,12
nsel,r,loc,x,0.050
sf,all,rdsf,0.68,1
allsel,all

stef,5.67e-8               ! Stefan-Boltzman constant
radopt,,0.01
spctemp,1,25               ! Space temperature
v2dopt,1
esel,s,mat,,2              ! select billet material
bfe,all,fvin,,1            ! define volumetric interface
allsel,all
allsel
asel,s,,,12              !铁区表面对流系数为70
lsel,r,,,all
lsel,u,loc,x,0
nsll,r,,,all
sf,all,conv,70,25
allsel
nsel,s,loc,x,0      !对称轴约束对流为0
sf,all,conv,0,25   
lsel,s,loc,y,0.025    !水上表面对流系数为0.54
nsll,r,,,all
sf,all,conv,0.54,25   
allsel
asel,s,,,1,10,1      !线圈四周对流系数为100
lsla,r,,,all
nsll,r,,,all   
sf,all,conv,100,25
allsel                 !空气场外围对流系数为0.024
lsel,s,loc,x,0.08
lsel,a,loc,y,-0.03
lsel,a,loc,y,0.07
nsel,r,,,all
sf,all,conv,0.024,25

allsel                  !空气场外围约束温度为25度
nsel,s,loc,x,0.08
nsel,a,loc,x,0
nsel,a,loc,y,-0.03
nsel,a,loc,y,0.07
d,all,temp,25
allsel
antype,trans               ! Thermal analysis options
trnopt,full
toffst,273
tunif,25                  ! initial uniform temperature
time,ftime                     
autots,on                   ! auto time-stepping
deltim,1,.5,1,on     ! time step control
kbc,1  
outres,all,all
solve
save