!Zylinderspule
finish
/clear
/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
h = 0.3e-3                 !Zylinderhöhe
!_____Spulenwerte
area = h*(r_au-r_in)
strom = 10e8*area !Spulenstrom  
!_____Luft
Breite = 0.4e-3             !Luftblock-Breite
Tiefe = Breite
Hoehe = Breite
!____________________________________________________________Geometrieerstellung
!_____Zylinder
cylind,r_in,r_au,h/2,-h/2,-90,90
cylind,r_in,r_au,h/2,-h/2,90,270
vadd,all
!_____Luftblock
block,Breite/2,-Breite/2,Hoehe/2,-Hoehe/2,Tiefe/2,-Tiefe/2        
vovlap,all            
numcmp,all
!/eof         
!_____________________________________________________________________Vernetzung
!_____Zylinderspule
esize,0.01e-3
mshape,1,3D
mshkey,0
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
vmesh,2
!_____Luftblock
type,1    
mat,1
mshkey,0
mshape,1,3d
esize,0.05e-3
vmesh,1

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
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
!________________________________________________________________________Loesung
!/eof
allsel
/solu
solve
/eof