finish  $/clear
*abbr,zylinder,/input,zylinder1,inp
!/filname,mag_feld
!/title,Magnetfeld

!--------Variablen------------
strom=1       !Strom durch Spulen in A
windz=23     !Windungsanzahl
wira=0.174  !Spulenradius
widi=0.0025 !Wicklungsdicke
wiho=0.0025 !Wicklungshoehe
beradin=0.25 !Radius der umgebenden Sphaere
strd=strom*windz/(wiho*widi)  !Stromdichte
/prep7

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

!-----Material------
mp,murx,1,1     !Luft
mp,murx,2,1     !Spule
mp,murx,3,50000 !Zylinder

!----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

!----Spule2-----------------------
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

!-------Zylinder---------------------

cyl4,0,0,0.003,0,0.0031,90,0.03
cyl4,0,0,0.003,90,0.0031,180,0.03
cyl4,0,0,0.003,180,0.0031,270,0.03
cyl4,0,0,0.003,270,0.0031,360,0.03
vglue,9,12
numcmp,volu

!------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
esize,ndiv
vmesh,1,8
!------Zylinder--------------
mat,3
esize,ndiv
vmesh,9,12
!------Luft-----------------
mat,1
mshkey,0
mshape,1,3d
esize,ndiv
lsel,s,loc,z,0    !Für Pfad-Auswertung
lesize,all,,,12
vmesh,13,16
!-----Rand--------------
nsel,s,ext
d,all,az,0
!-----Last--------------
allsel
asel,u,area,,118,124
asel,u,area,,65
aplot

vsel,s,volu,,3
vsel,a,volu,,10,12
vplot


eslv,s
eplot


local,11,1         !lokalisches zylindrisches Koordinatensystem
emodif,all,esys,11 !Hineindrehen der Elemente
bfe,all,js,,,strd  !Vorgabe der Stromdichte
csys,0

allsel

!-----Loesung----------
/solu
solve
!--------------------
!/post1
!nsel,s,loc,x,0,beradin
!esln,s,1
!plns,h,sum


!----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

allsel
asel,u,area,,118,124
asel,u,area,,65
aplot

vsel,s,volu,,3
vsel,a,volu,,10,12
vplot


eslv,s
eplot
plnsol,h,z