【讨论】求电磁感应加热的例子.

banana_wly

谢谢

kingzzz

/batch,list
/filenam, induc
/prep7
shpp,off
/title, induction heating of a solid cylinder billet
/com,
et,1,13,,,1             ! PLANE13, axisymmetric, AZ dof
et,2,13,,,1
et,3,151,,,1,1,1        ! SURF151, thermal, radiation
r,3,0                   ! Real constant set for SURF151

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=150000             ! frequency (Hz.)
pi=4*atan(1)            ! pi
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

emunit,mks           ! set magnetic units
mp,murx,1,1             ! air relative permeability
mp,murx,3,1             ! coil relative permeability
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

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

n                       ! create space node for SURF151   
*get,nmax,node,,num,max
lsel,s,loc,x,row
type,3
real,3
mat,2
lmesh,all            ! mesh billet outer radius with SURF151      
*get,emax,elem,,num,max
emodif,emax,3,nmax   ! modify element to add space node for radiation
et,3,0               ! reset type 3 to null element

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
esel,all
finish
/solu
antyp,harm
harfrq,150000
physics,write,emag      ! write emag physics file
finish

/prep7
lsclear,all             ! clear all b.c.'s and options
et,1,55,,,1             ! PLANE55 thermal element, axisymmetric
et,2,0                  ! null element type for coil and air region
et,3,151,,,1,1,1        ! SURF151 element for radiation
keyopt,3,9,1
r,3,1,5.67e-8           ! form factor, Stefan-Boltzmann constant

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
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
finish

banana_wly

那个例子我看了，只是热分析部分——辐射我没有看懂。不知能否赐教。另外，哪里还有这样的例子呢？

djwndxzn

我也想要一个ehong
xiaocuijiao@163.com   谢 谢！！

djwndxzn

怎么我一直没收到到你的例子，你给我发了吗？我现在急需一个这样的例子

mcn

发给我一份，谢谢。 email: ansys@0451.com

banana_wly

ansys@0451.com 这个邮箱好像不好用，给你发了两次都发不了。真抱歉。

perch88

给我发一个吧，谢谢！
perch88@sina.com

lanmiao

能发一个给我么
xw430@126.com,谢谢！

qhlty

发给我一份，谢谢。 email: a001@vip.163.com

myadfm

我也需要!谢谢!!
myadfm@sohu.com

myadfm

给我一份吧!我也非常需要!
谢谢了,

wangxia@mail.iee.ac.cn

banana_wly

#151436 - A 3d induction heating example

The main limitiations with this analysis type are 1) you can't have a
bh curve in a harmonic analysis and 2)  BH curves are not temperature
dependent.  Given these limitations, the most accurate thing that can be
done is to define a series of effective constant permeabilities, and
assign elements different material set numbers based on their H value.  
Of course, this would be an iterative process.

!
!  3d Induction Heating Example created at ANSYS 6.0
!
/prep7
et,1,117
et,2,117,1                ! volt dof needed in eddy current region
et,11,200,7               ! quad with mid-side nodes for meshing
mp,murx,1,1
mp,murx,2,1000            ! approximate iron properties
mp,rsvx,2,.2e-6  
wprot,,90                 ! build model in xz plane for later use of
csys,1
                         ! build model in inches and scale to meters at
the end
rect,.5,1,,1
rect,.5,1,1,2
rect,.5,1,2,3
rect,1,1.2,,1
rect,1,1.2,1,2
rect,1,1.2,2,3
rect,1.2,2,,1
rect,1.2,2,1,2
rect,1.2,2,2,3
nummrg,kp

lsel,a,,,3,11,4
LESIZE,all,,,20,10,,,,1
lsel,all

lesize,1,,,20,.1,,,,1

lesize,25,,,10,5

lsel,s,,,27,35,4
lesize,all,,,10,.2
lsel,all

FLST,5,4,4,ORDE,4   
FITEM,5,6   
FITEM,5,8   
FITEM,5,18  
FITEM,5,30  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,20, , , , ,1
!*  
FLST,5,8,4,ORDE,8   
FITEM,5,2   
FITEM,5,4   
FITEM,5,10  
FITEM,5,12  
FITEM,5,14  
FITEM,5,22  
FITEM,5,26  
FITEM,5,34  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,10,5, , , ,1
!*  
FLST,5,4,4,ORDE,4   
FITEM,5,2   
FITEM,5,12  
FITEM,5,14  
FITEM,5,26  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,10,.2, , , ,1
!*  
FLST,5,4,4,ORDE,4   
FITEM,5,13  
FITEM,5,15  
FITEM,5,19  
FITEM,5,23  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,2, , , , ,1  

type,11
esha,2
amesh,all   

k,101,
k,102,,,1
type,1
esize,,2
vrotate,all,,,,,,101,102,15,1

aclear,1,9,1
etdele,11   
vsel,s,,,2  
eslv
emod,all,mat,2  
emod,all,type,2
alls
local,11,1
emod,all,esys,11      ! so current is applied tangentially

vlscale,all,,,1/39.36,1/39.36,1/39.36,,,1     ! scale to meters
FINISH  

!!! Create 10 magnetic solutions at 10 coil locations
/filn,emag
/solu   
antype,harm
harf,,10000             ! excitation at 10 khz
nsel,s,ext
d,all,az,0              ! flux parallel at model edges
nsel,all
CoilLoc=1.15            ! parameter controlling coil location
*do,iii,1,10,1
CoilLoc=CoilLoc+.05     ! move coil .05 meters each new solution
bfedel,all,js           ! delete previous curr density load
vsel,s,,,5  
eslv
nsel,s,loc,z,CoilLoc/39.36,(CoilLoc+.1)/39.36   
esln,r,1
bfe,all,js,2,1e8     ! Arbitrary current density input, 1e8
amps/meter**2
alls
save
solve   
*enddo
fini

!!! Convert to thermal model and solve for the 10 heat generation cases
/filn,thermal
/prep7  
et,1,0           ! make non steel elements null
et,2,70          ! 8 node thermal brick
mp,kxx,2,50  
mp,dens,2,7800
mp,c,2,4
FINISH  
/solu   
antype,trans
ic,all,temp,70
tval=0
*do,jjj,1,10,1
ldread,hgen,jjj,1,,,emag,rst  
tval=tval+1      ! use 1 second between solutions, velocity = .05/1
Meters/sec.
time,tval
solve   
*enddo
fini

/exit

kingzzz

谢谢上传！