!eine  Wicklung (Kupfer) stromdurchflossen
!zweidimensional axialsysmmetrisch
!keine Symmetrienutzung
!Schätzing, 21.7.2010

finish  $/clear
*abbr,induk,/input,induk,txt

wiradin=30e-3  !Windungs-Radius innen
wiradau=36e-3  !Windungs-Radius außen
wiho=60e-3     !Windungs-Höhe
berad=100e-3   !Bereichs-Radius
strom=1        !Speisestrom der Wicklung
wianz=100      !Windungszahl
wifl=(wiradau-wiradin)*wiho  !Wicklungsquerschnitt
strdi=strom*wianz/wifl  !Stromdichte
rocu=57e6   
mw=(wiradau-wiradin)/2         !Maschenweite 

!------------------------------------
/prep7  
et,1,plane53,,,1 
et,2,infin110,,,1 
mp,murx,1,1        !Luft 
mp,murx,2,1        !Wicklung
mp,rsvx,2,1/rocu   !Wicklung

!----Geometrie------
pcirc,berad,,-90,90
rectng,wiradin,wiradau,-wiho/2,wiho/2
aovlap,all
!A2:Leiter A3:Bereich

!-------Vernetzung--Leiter------
esize,mw  !Maschenweite
mat,2
amesh,2
!----Vernetzung---Luft------
esize,mw*2
mat,1
amesh,3

!-------Rand-------
csys,1
nsel,s,loc,x,berad
esln
emodif,all,type,2
sf,all,inf
csys,0

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

!--------
/solu
solve
!--------------
/post1
plf2d

!------Induktivität-------
etable,energie,sene
ssum
*get,sumener,ssum,,item,energie
L_ansys=2*sumener/strom/strom

!----Näherungsbeziehung-----
pi=4*atan(1)
my0=12.56e-7
L_naeher=wianz**2*my0*pi*wiradin**2/wiho

!----Relation-----
rel=L_naeher/L_ansys