3D感应加热

zws77826

那位大虾能够上传或Email一个3D感应加热的实例啊，小第正心急如焚的做客提，将不胜感激。
Email:zws77826@163.com

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

bopor

这个例子不行啊，根本运行不下去，error太多。

robbot1980

运行很好，没问题，不知计算结果如何

xulang1013

有没有GUI的啊