!
!eine Wicklung über Platte
!
!zweidimensional a
!Induktionserwärmung
!Schätzing, 31.8.2012
!
finish  $/clear
*abbr,inderw,/input,inderw,txt


!------Parameter--------
brpla=60e-3      !Breite Platte
dipla=10e-3      !Dicke Platte
radwi=30e-3      !Radius Wicklung
raddra=1.58e-3   !Radius Draht/Drahtbündel
wiho=raddra+1e-3 !Höhe über Platte
wizahl=10        !Windungszahl
strom=150
strdi=strom*wizahl/(3.14*raddra)

beradin=80e-3    !Bereichs-Radius innen
beradau=100e-3   !Bereichs-Radius außen


/prep7  
et,1,plane13     !Luft, 
et,2,plane13     !Wicklung
et,3,plane13     !Platte
et,4,infin110    !Randelement

mp,murx,1,1          !Luft 
mp,murx,2,1          !Wicklung
!mp,rsvx,2,1/56e6
mp,kxx,2,350
mp,c,2,380
mp,dens,2,8900
mp,murx,3,1          !Platte
mp,rsvx,3,1/1.4e6
mp,kxx,3,52
mp,c,3,420
mp,dens,3,7860

!----Geometrie------
pcirc,0,beradau,  0, 90
pcirc,0,beradau, 90,180
pcirc,0,beradau,180,270
pcirc,0,beradau,270,360
pcirc,0,beradin


rectng,-brpla/2,brpla/2,0,-dipla

local,11,0,0,wiho
wpcsys,,11
pcirc,,raddra,  0, 90
pcirc,,raddra, 90,180
pcirc,,raddra,180,270
pcirc,,raddra,270,360



aovlap,all
numcmp,area
!A9,10:Platte A1-4:Wicklung   
!A11-14:Bereich innen A5-8:Bereich außen
aplot
!-------Vernetzung--Leiter------
!/eof
mshape,0,2d  !Vernetzen mit Vierecken
mshkey,1
esize,raddra/2  !Maschenweite
type,2  $mat,2 
amesh,1,4
!-------Platte-------
type,3  $mat,3  $amesh,9,10  !Platte

!----Vernetzung---Luft------
csys,1
lsel,s,loc,x,beradin,beradau
lsel,u,loc,x,beradin
lsel,u,loc,x,beradau
!lesize,all,,,4,1/4
lesize,all,,,1
esize,,9
type,1   $mat,1
amesh,5,8

mshkey,0
esize
type,1  $mat,1
amesh,all
!----Rand-----
!/eof
nsel,s,loc,x,beradau
esln
emodif,all,type,4
sf,all,inf
!------------
csys,0
!---Kopplung-------------
!esel,s,mat,,2
!nsle
!*get,nanf,node,,num,min
!cp,1,volt,all
!----Platte-geschlitzt----
!esel,s,mat,,3
!nsle,s,1
!cp,2,volt,all
!allsel

!/eof

!---Stromdichtevorgabe----------
!f,nanf,amps,strom*wizahl
esel,s,mat,,2
bfe,all,js,,,,strdi*100


!--------
/solu
allsel
antype,harm
harfre,360e3
outres,all,all

solve



!--------------
/post1
/plopts,minm,0
set,1,1   !Realteil
plf2d



!/eof




esel,s,mat,,2,3
ples,jheat

esel,s,mat,,3
etable,strdi,jt,sum
pletab,strdi,avg

esel,s,mat,,2,3
etable,verlu,jheat
pletab,verlu,avg

!------thermisch---------

/prep7
allsel
lsclear,all

et,1,0
et,2,0
et,3,plane13,2
et,4,0

!transiente Lösung
/solu
antype,trans
tunif,20        !Umgebungstemperatur
outres,all,all
esel,s,mat,,3
ldread,hgen,,,,2,,rst
time,10
solve

!---------------------------
/post1
esel,s,mat,,3
plns,temp