!Zylinderspule
finish
/clear
*abbr,RobinG3Dmag2,/input,RobinG3Dmag2,txt

/prep7
!___________________________________________________________________Elementtypen
et,1,solid97          !Luft  
et,2,solid97,3          !Spule
et,3,circu124,3,0
et,4,circu124,5
!__________________________________________________________Materialeigenschaften
mp,murx,1,1          ! Luft
mp,murx,2,1          ! Spule
mp,rsvx,2,1/56e6
!_____________________________________________________________Geometrieparameter
pi = 3.14159265
!_____Zylinder
r_in = 0.09e-3             !Innenradius 
r_au = 0.1e-3              !Auáenradius
r_mi=0.095e-3
h = 0.3e-3                 !Zylinderhöhe
!_____Spulenwerte
area = h*(r_au-r_in)
strdi=1000e6
strdi=500e6                 !Stromdichte
strom = strdi*area         !Spulenstrom  
!_____Luft
Breite = 0.4e-3             !Luftblock-Breite
Tiefe = Breite
Hoehe = Breite
!____________________________________________________________Geometrieerstellung
!_____Zylinder

cylind,r_in,r_mi,h/2,-h/2,  0, 90
cylind,r_in,r_mi,h/2,-h/2, 90,180
cylind,r_in,r_mi,h/2,-h/2,180,270
cylind,r_in,r_mi,h/2,-h/2,270,360
cylind,r_mi,r_au,h/2,-h/2,  0, 90
cylind,r_mi,r_au,h/2,-h/2, 90,180
cylind,r_mi,r_au,h/2,-h/2,180,270
cylind,r_mi,r_au,h/2,-h/2,270,360

vglue,all
!_____Luftblock
block,Breite/2,-Breite/2,Hoehe/2,-Hoehe/2,Tiefe/2,-Tiefe/2        
vovlap,all            
numcmp,all
!V2:Luft   V1,3bis9 Spule
!/eof         
!_____________________________________________________________________Vernetzung
!_____Zylinderspule
esize,h/6

type,2
mat,2
!r,Nr,spulenflächen,windungszahl,volumen,x-richtung,y-richtung,z-richtung
r,1,h*(r_au-r_in),8,pi*h*(r_au*r_au-r_in*r_in),0,1,0 ! Stromrichtung definieren
real,1
vmesh,1
vmesh,3,9

!_____Luftblock
type,1    
mat,1
mshkey,0
mshape,1,3d
!esize,0.05e-3
vmesh,2

esel,s,type,,2
esel,r,real,,1
nsle
local,11,1,-x_Wick_Mit,dicke_magnetkopf/2,Y_Wick_Mit
emodif,all,esys,11      !Zylinderkoordinaten werden eingeschaltet
csys,0
allsel
!/eof
!_________________________________________________________________________Lasten
esel,s,type,,2
nsle
cp,1,curr,all     !Kopplung
cp,2,emf,all
*get,node_sp1,node,,num,max

allsel
numcmp,node
*get,node1,node,,num,max
n,node1+1,Breite+1500e-6,0,0
n,node1+2,Breite+2000e-6,500e-6,0
n,node1+3,Breite+2500e-6,0,0
!/eof
type,3
r,2,-strom/100
real,2
e,node1+1,node1+2

type,4
r,3,1
real,3
e,node1+2,node1+3,node_sp1
/icscale,,0.5

eplot
allsel
nsel,s,node,,node1+1,node1+3,2
d,all,volt,0
allsel
eplot
!------Randbedingung------
esel,s,type,,1,2
nsle
nsel,r,ext
d,all,az,0

!________________________________________________________________________Loesung
!/eof
allsel
/solu
solve
!/eof
!-------------------------------
/post1
nsel,s,loc,x,-breite/2,0
esln,s,1
esel,r,mat,,2

ples,jheat
plve,jt

!----Thermische Berechnung------

!/eof

/prep7
allsel
!___________________________________________________________________Elementtypen
et,1,0
et,3,0
et,4,0
!et,2,solid5,8
et,2,solid70
!__________________________________________________________Materialeigenschaften
!mp,kxx,1,0.2          ! Luft
mp,kxx,2,400          ! Spule
mp,dens,2,8900
mp,c,2,380
esel,s,type,,2
!LDREAD, Lab, LSTEP, SBSTEP, TIME, KIMG, Fname, Ext, --
!ldread,hgen ,      ,       ,     ,     ,RESULT,RST !Name der Result-Datei
ldread,hgen
!allsel
tunif,20
nsle
nsel,r,ext
sf,all,conv,8,20
esel,s,type,,2
nsle

!-------Lösung-------
/solu
!antype,trans
!time,2
!deltim,0.2
!outres,all,all
solve

!-------Auswertung------
/post1
plns,temp