感应加热 温度上不去

xulang1013

/batch,list
/title, Induction heating of a solid cylinder billet
/prep7
shpp,off
/com,
/com,
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
freq=10000                 ! frequency (Hz.)
pi=4*atan(1)               ! pi
cond=1.9e6                 ! 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
emunit,mks                 ! set magnetic units
mp,murx,1,1                ! air relative permeability
mp,murx,3,1                ! coil relative permeability
mptemp,1,20,100,200,300,400     ! temps for relative permeability
mpdata,murx,2,1,100,97.5,93.3,89,83.6   ! steel relative permeability
mptemp
mptemp,1,20,100,200,300,400            ! temps for resistivity
mpdata,rsvx,2,1,.52e-6,.59e-6,.65e-6,.71e-6,.8e-6
mptemp
mptemp,1,20,100,200,300,400                  ! temps for conductivity
mpdata,kxx,2,1,32.2,30.9,29.9,29.5,29.2
mptemp                                   
mptemp,1,20,100,200,300,400
mpdata,c,2,1,460,469,480,502,534
mptemp                                   
mptemp,1,20,100,200,300,400
mpdata,dens,2,1,7800,7800,7800,7800
mp,emis,2,.68              ! emissivity
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
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
mshape,0,2d
mshk,1
amesh,1                    ! mesh billet area
lsel,s,loc,y,0
lsel,a,loc,y,t
lsel,u,loc,x,row/2
lesize,all,.001
lsel,all
amesh,all                  ! mesh remaining areas
nsel,s,loc,x
d,all,az,0                 ! apply flux-normal b.c.
nsel,all
esel,s,mat,,3
bfe,all,js,,,,1.87e6         ! apply current density to coil
alls
! 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
alls
spctemp,1,20               ! Space temperature
v2dopt,1
radopt,,0.01
stef,5.67e-8               ! Stefan-Boltzman constant
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,10000
outres,all,all
tunif,20
mfcm,1,                    ! Write Emag analysis options
mfclear,solu               ! Clear analysis options
antype,trans               ! Thermal analysis options
toffst,273
tunif,20                  ! 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
finish
/post26
file,field2,rth
esel,s,type,,4
nsle,s
nsel,r,loc,x,row           ! get node at outer radius
nsel,r,loc,y,0
*get,nor,node,,num,min
nsle,s
nsel,r,loc,x,0
nsel,r,loc,y,0             ! get node at centerline
*get,nir,node,,num,min
nsol,2,nor,temp,,outerR    ! Outer radius
nsol,3,nir,temp,,inner     ! Inner radius (centerline)
/axlab,y,Temperature
plvar,2,3                  ! plot temperature
prvar,2,3
finish
/post1
file,field2,rth
set,last                   ! Solution at 3 seconds
esel,s,type,,4             ! select thermal elements
plns,temp                  ! plot temperature
finish
就把帮助的参数载荷换成了自己的了，算了3秒温度上去不 才3度，请问怎么回事啊，谢谢。

weizumin

3秒应该也差不多就升3度左右吧，你把时间加长计算看能不能升啊？如果只是改了电流密度的加载应该是没问题的。

xulang1013

我把材料的参数也换成自己的了，这样在网格划分上是不会出问题是

清小溪

你好！！你这个温度 不升的问题解决没有啊？我的感应加热也是出现这个问题，温度升不上去