!eine  Wicklung (Kupfer) stromdurchflossen
!dreidimensional 
!keine Symmetrienutzung
!Prüfung von MrTorque-Ergebnis

fini  $/clear
*abbr,wicklung_3D,/input,wicklung_3D,txt

wiradin=95e-3   !Windungs-Radius innen
wiradau=105e-3  !Windungs-Radius außen
wiho=10e-3      !Windungs-Höhe
berad=240e-3    !Bereichs-Radius
strom=200
wifl=(wiradau-wiradin)*wiho  !Wicklungsquerschnitt
strdi=strom/wifl  !Stromdichte
rocu=57e6   
mw=wiho/2         !Maschenweite 
/prep7  
et,1,solid97, 
!et,2,infin110,,,1 
mp,murx,1,1        !Luft 
mp,murx,2,1        !Kupfer
mp,rsvx,2,1/rocu   !Kupfer

!----Geometrie------
local,11,1
wpcsys,,11
cylind,wiradin,wiradau,-wiho/2,wiho/2!,  0, 90
!cylind,wiradin,wiradau,-wiho/2,wiho/2, 90,180
!cylind,wiradin,wiradau,-wiho/2,wiho/2,180,270
!cylind,wiradin,wiradau,-wiho/2,wiho/2,270,360

sphere,,berad,  0, 90
sphere,,berad, 90,180
sphere,,berad,180,270
sphere,,berad,270,360

vovlap,all
numcmp,volu
!V1-4:Wicklung  v5-8:Luft

!-------Vernetzung--Leiter------
esize,mw  !Maschenweite
mat,2
vmesh,1,4
emodif,all,esys,11

!----Vernetzung---Luft------
esize,mw*16
esize,mw*4

mat,1
mshape,1,3d
mshkey,0
vmesh,all

!-------Rand-------
!aus Bequemlichkeit ohne Randelement
!umgebender Luftraum zu klein
nsel,s,ext
d,all,ax,0
d,all,ay,0
d,all,az,0

!-----Stromdichtevorgabe-------------
esel,s,mat,,2
bfe,all,js,,,strdi   !y

allsel
!--------
/solu
solve
!--------------
/post1
csys,0
nsel,s,loc,x,0,berad
esln,s,1
ples,b,sum

/wait,1

/show,,,1
/edge,1,1
plve,b





n1=node(0,0,0)
*get,bx0,node,n1,b,x
*get,by0,node,n1,b,y
*get,bz0,node,n1,b,z


!Ziel von Mrtorque: 1.26e-3
my0=4*3.1415e-7
b=strom*my0/(2*100e-3)