一个经典的电磁感应加热求温度场分布

sunclark

FINISH
/CLEAR  
/title, Induction heating
/units,si
/prep7

row=.015                   ! outer radius of workpiece
ric=.0175                  ! inner radius of coil
roc=.0200                  ! outer radius of coil
ro=.05                     ! outer radius of model
t=.001                     ! model thickness

emunit,mks                 ! set magnetic units
freq=150                   ! frequency (Hz.)
pi=4*atan(1)               ! pai
cond=.392e7                ! maximum conductivity
muzero=4e-7*pi             ! free-space permeability
mur=200                    ! maximum relative permeability
skind=sqrt(1/(pi*freq*cond*muzero*mur))  ! skin depth

ftime=3                    ! final time
tinc=.05                   ! time increment for harmonic analysis
time=0                     ! initialize time
delt=.01                   ! maximum delta time step

! Electromagnetic model
et,1,53,,,1                ! PLANE53, axisymmetric, AZ dof
et,2,53,,,1

mp,murx,1,1                ! air relative permeability
mp,murx,3,1                ! coil relative permeability
! temps for relative permeability
mptemp,1,25.5,160,291.5,477.6,635,698,709,720.3,742,761,1000
! steel relative permeability
mpdata,murx,2,1,200,190,182,161,135,104,84,35,17,1,1
! temps for resistivity
mptemp,1,0,125,250,375,500,625,750,875,1000
! steel resistivity
mpdata,rsvx,2,1,.184e-6,.272e-6,.384e-6,.512e-6,.656e-6,.824e-6,1.032e-6,1.152e-6,1.2e-6
! temps for conductivity
mptemp,1,0,730,930,1000                  
mpdata,kxx,2,1,60.6,29.5,28,28
! temps for enthalpy
mptemp,1,0,27,127,327,527,727,765,770,927
mpdata,enth,2,1,0,91609056,453285756,1.2748e9,2.2519e9,3.3396e9,3.548547e9,3.548556e9,4.3520e9
! emissivity
mp,emis,2,.68
! Stefan-Boltzman constant  
stef,5.67e-8           

!problem modeling
rectng,0,row,0,t           ! billet
rectng,row,ric,0,t         ! air-gap
rectng,ric,roc,0,t         ! coil
rectng,roc,ro,0,t          ! outer air
aglue,all
numcmp,area
/PNUM,AREA,1
save

ksel,s,loc,x,row           ! select keypoints at outer radius of workpiece
kesize,all,skind/2         ! set meshing size to 1/2 skin depth
ksel,s,loc,x,0             ! select keypoints at center
kesize,all,40*skind        ! set meshing size
lsel,s,loc,y,t/2           ! select vertical lines
lesize,all,,,1             ! set 1 division through thickness
lsel,all

asel,s,area,,1
aatt,2,1,1                 ! set attributes for billet region
asel,s,area,,3
aatt,3,1,2                 ! set attributes for coil region
asel,s,area,,2,4,2
aatt,1,1,2                 ! set attributes for air region
asel,all
save
mshape,0,2d
mshk,1
! mesh billet area
amesh,1
save                    

! mesh remaining areaslsel,s,loc,y,0
lsel,a,loc,y,t
lsel,u,loc,x,row/2
lesize,all,.001
lsel,all
amesh,all
               
nsel,s,loc,x
d,all,az,0                 ! apply flux-normal b.c.
nsel,all
esel,s,mat,,3
bfe,all,js,,,,15e6         ! apply current density to coil
ALLSEL,ALL

! Thermal model
et,4,55,,,1                ! PLANE55 thermal element, axisymmetric

agen,2,1,,,,,,,1
aatt,2,1,4
ksel,s,loc,x,row           ! select keypoints at outer radius of workpiece
kesize,all,skind/2         ! set meshing size to 1/2 skin depth
ksel,s,loc,x,0             ! select keypoints at center
kesize,all,40*skind        ! set meshing size
lsel,s,loc,y,t/2           ! select vertical lines
lesize,all,,,1             ! set 1 division through thickness
lsel,all
mshape,0,2d
mshk,1
amesh,5
asel,,,,5
allsel,below,area
nsel,r,loc,x,row
sf,all,rdsf,0.68,1         ! Radiation
ALLSEL,ALL
! Space temperature
spctemp,1,25               
v2dopt,1
radopt,,0.01
esel,s,mat,,2              ! select billet material
bfe,all,fvin,,1            ! define volumetric interface
finish

/solu
mfan,on                    ! Activate ANSYS Multi-field solver analysis
mfel,1,1,2                 ! Field #1 ET;s,  Emag
mfel,2,4                   ! Field #2 ET's, Thermal
mfor,1,2                   ! Field solution order
mfti,ftime                 ! Final time
mfdt,tinc                  ! Stagger time increment
mfco,all,1e-3              ! Convergence criteria

antyp,harm                 ! Emag analysis options
harfrq,150
outres,all,all
tunif,100
mfcm,1,                    ! Write Emag analysis options
mfclear,solu               ! Clear analysis options

antype,trans               ! Thermal analysis options
toffst,273
tunif,100                  ! initial uniform temperature
kbc,1                      ! step loads
trnopt,full
autos,on                   ! auto time-stepping
deltim,.01,.005,.01,on     ! time step control
mfcm,2,                    ! Write Thermal analysis options       

mfvo,1,1,hgen,2            ! Transfer hgen from Emag to Thermal
mfvo,1,2,temp,1            ! Transfer Temp from Thermal to Emag
solve
save

nicejoy

这就是help里面的例子啊，但是这里面考虑了热辐射……我在做的电磁-热耦合分析中，没有考虑辐射。2d分析使用comsol分析的，但是软件内存要求较高，3d模拟只能用ansys模拟了……