!Homogenes Magnetfeld zwischen Helmholtzspulen
!dreidimensional

finish  $/clear

!/filname,3d_leer
!/title,3-D Helmholtzanordnung

!--------Variablen------------
I=5       !Strom durch Spulen in A
N=100     !Windungsanzahl
!Empfehlung: Variablen-Namen immer mehr als ein Zeichen lang
wira=0.4  !Spulenradius
widi=0.03 !Wicklungsdicke
wiho=0.03 !Wicklungshoehe
beradin=2 !Radius der umgebenden Sphaere
prfrad=0.07
prfdi=0.005
prfho=0.005
wkappa=56e6
strd=i*n/(wiho*widi)  !Stromdichte
mu0=3.1415926*4.0E-7
jx=0
jy=strd
jz=0
/prep7

emunit,mks
et,1,solid117,   !Luft und Wicklung
et,2,solid117,1	

!-----Material------
mp,murx,1,1     !Luft
mp,murx,2,1     !Spule
mp,murx,3,1
mp,rsvx,3,1/wkappa

!----Geometrie----------
/view,,1,2,3
!----Spule1-----------------------
cylind,wira-widi/2,wira+widi/2,wira/2-wiho/2,wira/2+wiho/2,  0, 90
cylind,wira-widi/2,wira+widi/2,wira/2-wiho/2,wira/2+wiho/2, 90,180
cylind,wira-widi/2,wira+widi/2,wira/2-wiho/2,wira/2+wiho/2,180,270
cylind,wira-widi/2,wira+widi/2,wira/2-wiho/2,wira/2+wiho/2,270,360
vglue,1,4
numcmp,volu
!----Spule1-----------------------
cylind,wira-widi/2,wira+widi/2,-wira/2-wiho/2,-wira/2+wiho/2,  0, 90
cylind,wira-widi/2,wira+widi/2,-wira/2-wiho/2,-wira/2+wiho/2, 90,180
cylind,wira-widi/2,wira+widi/2,-wira/2-wiho/2,-wira/2+wiho/2,180,270
cylind,wira-widi/2,wira+widi/2,-wira/2-wiho/2,-wira/2+wiho/2,270,360
vglue,5,8
numcmp,volu

!---Pruefling
cylind,prfrad-prfdi/2,prfrad+prfdi/2,-prfho/2,prfho/2,0,360

!------Luftumgebung--------------
sphere,0,beradin,  0, 90
sphere,0,beradin, 90,180
sphere,0,beradin,180,270
sphere,0,beradin,270,360
vovlap,all
numcmp,volu
!------Vernetzung---------------------
!------Spulen-----------------
type,1
mat,2
vmesh,1,8

!---Pruef
type,2
mat,3
smrtsize,2
vmesh,9,12

!------Luft-----------------
type,1
mat,1
mshkey,0
mshape,1,3d
esize,beradin/4
lsel,s,loc,z,0    !Fr Pfad-Auswertung
lesize,all,,,12
vmesh,13,16			!9,12

!-----Rand--------------
nsel,s,ext
d,all,az,0

!-----Last--------------
esel,s,mat,,2
local,11,1         !lokalisches zylindrisches Koordinatensystem
emodif,all,esys,11 !Hineindrehen der Elemente
bfe,all,js,,jx,jy,jz  !Vorgabe der Stromdichte
csys,0

esel,s,mat,,3
nsle
cp,1,volt,all
*get,nanf,node,,num,min
f,nanf,amps,0

allsel


!-----Loesung----------
/solu
antype,3
harfreq,50,50
nsubst,,
kbc,0
solve

!--------------------
/post1
set,1,1
nsel,s,loc,x,0,beradin
esln,s,1
plns,b,sum

/eof

!----Feldverlauf entlang der Z-Achse
/post1
allsel
path,zachse,3
ppath,1,,0,0,beradin
ppath,2,,0,0,0
ppath,3,,0,0,-beradin
pdef,hzachs,h,sum
/grid,1
plpath,hzachs