三维感应加热问题求助，请高手指点

qiuweihuam

以下是我三维感应加热的命令流，不知为什么结果总是不对，温度老是不变，无法加热，望高手指点，不胜感激



/filenam, induc
/prep7

/title, induction heating of a solid cylinder billet

!et,1,solid62            ! PLANE13, axisymmetric, AZ dof
!et,2,solid62
!et,3,152,,,1,1,1        ! SURF151, thermal, radiation
!KEYOPT, 1, 1, 0
!KEYOPT, 1, 6, 0
!KEYOPT, 1, 5, 0

et,1,solid97            ! PLANE13, axisymmetric, AZ dof
et,2,solid97
et,3,152,,,1,1,1        ! SURF151, thermal, radiation
KEYOPT, 2, 1, 0
!KEYOPT, 2, 2, 0
!KEYOPT, 2, 5, 0

r,3,0                   ! Real constant set for SURF151
r,2,0

!r,1,0.015*0.02,1,1,0,1,0
!rmore,1,1






ftime=20            ! final time
tinc=5                ! time increment for harmonic analysis
time=0                  ! initialize time
delt=1                ! maximum delta time step

emunit,mks           ! set magnetic units
mp,murx,1,1             ! air relative permeability
mp,rsvx,1,.435e-6
mp,murx,3,1             ! coil relative permeability
mptemp,1,20,100,200,300,400,500,   ! temps for relative permeability
mptemp,7,600,700,720,742,760,800,900,1000
mpdata,murx,2,1,200,195,186.6,178.1,167.1,154.9, ! steel relative permeability
mpdata,murx,2,7,137.8,92.5,33,11,1,1,1,1
mptemp
mptemp,1,20,100,200,300,400,500,   ! temps for relative permeability
mptemp,7,600,700,720,742,760,800,900,1000
mpdata,rsvx,2,1,.198e-6,.254e-6,.339e-6,.435e-6,.541e-6,.656e-6
mpdata,rsvx,2,7,0.79e-6,0.949e-6,0.982e-6,1.019e-6,1.042e-6,1.080e-6 ! steel resistivity

cylind,0.046,0.048,0,0.162,-45,45   ! Build thermal model
cylind,0.048,0.05,0,0.162,-45,45   ! Build thermal model
cylind,0.05,0.106,0,0.162,-45,45    ! Build thermal model
numcmp,volu

vsel,s,VOLU,,3
aatt,2,2,1              ! set attributes for billet region

!vmesh,4
vsel,s,VOLU,,2
aatt,1,2,2              ! set attributes for air region

!vmesh,3
!vmesh,2
vsel,s,VOLU,,1
aatt,3,1,2              ! set attributes for coil region
!vmesh,1

csys,1
vsel,s,VOLU,,1
aslv
lsla
lsel,r,loc,x,0.046
lesize,all,,,20,,,,1    ! 20 bricks along bl
lsla
lsel,r,loc,x,0.048
lesize,all,,,20,,,,1    ! 20 bricks along bl
vmesh,1

vsel,s,VOLU,,2
aslv
lsla
lsel,r,loc,x,0.048
lesize,all,,,20,,,,1    ! 20 bricks along bl
lsla
lsel,r,loc,x,0.05
lesize,all,,,20,,,,1    ! 20 bricks along bl

vmesh,2


vsel,s,VOLU,,3
aslv
lsla
lsel,r,loc,x,0.05
lesize,all,,,40,,,,1    ! 20 bricks along bl
lsla
lsel,r,loc,x,0.0106
lesize,all,,,10,,,,1    ! 20 bricks along bl
vmesh,3


numcmp,node


n                       ! create space node for SURF151   
*get,nmax,node,,num,max
!csys,1
!asel,s,,,24,25
type,3
real,3
mat,2
n,100000,0,0,0.17
esurf,100000            ! mesh billet outer radius with SURF151  
et,3,0   


ASEL, S,,,14
NSLA, S, 1
D,ALL,aZ, 0
ALLSEL             ! apply flux-normal b.c.
nsel,all



vsel,s,VOLU,,1
nslv,s
cp,1,curr,all
!cp,2,emf,all
vsel,s,VOLU,,1
eSLv,S  
bfe,all,js,,,,100e6      ! apply current density to coil
!bfe,all,vltg,,220
ALLSEL
esel,all
!et,1,0
finish
/solu
antyp,harm
harfrq,2500
physics,write,emag      ! write emag physics file
finish


/prep7
lsclear,all             ! clear all b.c.'s and options
et,1,solid70            
et,2,0                  ! null element type for coil and air region
et,3,152,,,1,1,1        ! SURF152 element for radiation
KEYOPT, 1, 2, 0
KEYOPT, 1, 4, 0
KEYOPT, 1, 7, 0
KEYOPT, 1, 8, 0
keyopt,3,9,1

r,3,1,5.67e-8           ! form factor, Stefan-Boltzmann constant


mp,kxx,1,20             ! air relative permeability
mptemp
mptemp,1,20,100,200,300,400,500,   ! temps for relative permeability
mptemp,7,600,700,720,742,760,800,900,1000
mpdata,kxx,2,1,47.68,43.53,40.44,38.13,36.02,34.16, ! steel relative permeability
mpdata,kxx,2,7,31.98,28.66,27.49,26.2,25.14,26.49,25.92,24.02
mptemp                          ! temps for enthalpy
mptemp,1,0,27,127,327,527,727
mptemp,7,760,760.001,927
mpdata,enth,2,1,0,96526059,4.7589e8,1.2800e9,2.1607e9,3.3507e9
mpdata,enth,2,7,3.6063e9,3.6098e9,4.5639e9
mp,emis,2,.68                   ! emissivity

finish
/solu
antype,trans
toffst,273
tunif,100                      ! initial uniform temperature
d,nmax,temp,25                  ! ambient temperature
cnvtol,heat,1                   ! convergence tolerance
kbc,1                           ! step loads
trnopt,full
autos,on                        ! auto time-stepping
deltim,1e-5,1e-6,delt,on        ! time step control
outres,basic,all                ! save all load step information
physics,write,thermal           ! write thermal physics file
finish

*do,i,1,ftime/tinc   ! solution *do loop
time=time+tinc                            ! increment time
physics,read,emag    ! read emag physics file
/solu
*if,i,eq,1,then
tunif,100                    ! initial temperature
*else
ldread,temp,last,,,,,rth     ! read thermal analysis temperatures
*endif
solve                         ! solve harmonic analysis
finish
physics,read,thermal          ! read thermal physics file
/assign,esav,therm,esav       ! redirect files for use in thermal restart
/assign,emat,therm,emat
/solu
*if,i,gt,1,then
antype,trans,rest            ! thermal restart
*endif
time,time                     ! time at end of thermal run
esel,s,mat,,2                 ! select billet region
ldread,hgen,,,,2,,rmg         ! apply coupled joule heating load from emag
esel,all
solve
finish
/assign,esav                  ! reassign files to default
/assign,emat
*enddo
                      ! end of solution looping
/post1
set,,,1,,10
plnsol,temp,,0
!gettemp              !use macor,output node temperature

/post26

nsel,all


nsol,2,2,temp
plvar,2