三维有限有电磁热耦合瞬态分析

986787489

三维有限有电磁热耦合瞬态分析求APDL代码……

986787489

/title, tongweidou
/nopr
/VIEW,1,1,1,1
/VUP,1,Z
/SHOW,win32c     ! use /SHOW,x11c for Unix
/CONT,1,18

! *** Problem parameters

! Ferromagnetic plate:
lp=0.15          ! widen, m
hp=1             ! height, m

eszp=1           ! element size, m
sigp=6.66e6      ! conductivity, S/m
mu_fer=200       ! relative permeability
k=50             ! thermal conductivity, W/m_C   
h=100            ! convection coefficient, W/(m^2_C)

! Coil:
ri=0.23          ! inner radius, m
ro=0.23          ! outer radius, m
hc=1             ! height, m
d=0.028          ! distance between the coil and the plate,m
eszc=1           ! element  size, m
frq=50           ! frequency, Hz
I=1e5            ! rms coil current, A
mu_coil=1        ! relative permeability

! Air box:
la=0.5           ! length, m
ha1=0.05         ! lower height, m
ha2=0.35         ! upper height, m
esza=1           ! element size, m
mu_air=1         ! relative permeability
sigc=2e-8        ! conductivity, S/m
Tamb=70          ! ambient temperature,C

/com,  !实体元件
/com, *** Eddy current analysis
/com,


/PREP7
et,1,236,1       ! AZ-VOLT - ferromagnetic plate
et,2,236,1,,,,1  ! AZ-VOLT (no eddy currents) - coil
et,3,236         ! AZ - air
et,4,236
et,5,236

/PREP7
blc4,0,0,0.128,0.128,1
blc4,0,0,0.103,0.103,1
VSEL,S,,,1
VSEL,A,,,2
CM,volum,VOLU
ALLSEL,ALL
VSBV,volum,2                        ! coil
blc4,0,0,0.075,0.075,1              ! plate

!block,0,lp/2,0,lp/2,0,lp/2         ! plate
!cylinder,ri,ro,hp+d,hp+d+hc,0,1    ! coil
block,0,0.25,0,0.25,-0.05,0.35      ! air
vovlap,all
numcmp,all

mp,murx,1,mu_fer
mp,rsvx,1,1/sigp

mp,murx,2,mu_coil
mp,rsvx,2,sigc

mp,murx,3,mu_air

type,1
mat,1
esize,eszp
vmesh,1

type,2
mat,2
esize,eszc
vmesh,2

esha,1
type,3
mat,3
esize,esza
vmesh,3

asel,s,,,11
nsla,s,1
cp,1,volt,all
f,ndnext(0),amps,I*sqrt(2)  !  rms

asel,s,,,12
nsla,s,1
d,all,volt,0
nsel,all

asel,s,,,3
asel,a,,,5
nsla,s,1
d,all,volt,0
nsel,all

nsel,s,ext
d,all,az,0     ! Flux-parallel magnetic BC
nsel,all
allsel
fini

/SOLU
antyp,harm
harf,,frq
solve
fini

/POST1
set,1,1                  
esel,s,mat,,1             ! select the plate
POWERH                    ! calculate rms power loss
plve,jt,,,,vect,elem,on   ! plot eddy currents (real part)
set,1,1,,1
plve,jt,,,,vect,elem,on   ! plot eddy currents (imag part)
alls
fini

/com,
/com, Ansys computed power loss = %Pavg*4% Watts

/com,
/com, *** Thermal analysis
/com,

/PREP7
ddele,all,az
vsel,s,volu,,1
eslv
et,1,90
LDREAD,hgen
vsel,inve
vclear,all
alls
etdele,2,3

mp,kxx,1,k
nsel,s,ext
sf,all,conv,h,Tamb
nsel,all
fini

/SOLU
anty,static
solve
fini

/POST1
set,1,1
plnsol,temp
nsort,temp
*get,minTemp,sort,,min
*get,maxTemp,sort,,max
avTemp=(minTemp+maxTemp)/2
/com,
/com, Average temperature of the plate = %avTemp% (C)
/com,
fini
试试对不对大牛们

mzwpll

楼上高手

不爱洗澡的鱼

楼主，你好，我最近也在做三维电磁感应加热，希望可以一起交流一下。
qq:995794758