finish
/clear

leila=100e-6
leirad=10e-6
berad=50e-6
abstand=10e-6
strom=1

!Parameter
mw=2e-6

! Elementtypen
/prep7
et,1,solid97	!Luft
et,2,solid69	!Leiter
et,3,infin111,1	!Randelement

r,2,leirad*leirad*3.1416

!Materialeigenschaften
mp,murx,1,1
mp,murx,2,1
mp,rsvx,2,1/40e6

/view,,1,2,3

! Geometrie
/triad,off
!block,-leibrei/2, leibrei/2, -leiho/2, leiho/2, 0, leila
cylind,0,leirad,0,leila
cylind,0,berad,-abstand,leila+abstand
vovlap,all

!V1: Leiter	V3: Umgebung

! Vernetzung-Leiter
esize,mw
mshkey,0	
mshape,1,3d	
type,2 $mat,2 $real,2
vmesh,1

! Vernetzung-Luft
csys,1
asel,s,loc,x,berad
aesize,all,mw*10
allsel
type,1 $mat,1
mshkey,0
mshape,1,3D
vmesh,3

! Randelemente
csys,1
type,3
esize,,1
asel,s,loc,x,berad
*get,anr1,area,,num,min	!A9
*get,anr2,area,,num,max	!A10
vext,anr1,anr2,,mw
nsel,s,loc,x,berad+mw/2,berad+mw
sf,all,inf
csys,0
allsel

!Stromeinprägung
esel,s,mat,,2
nsle
nsel,r,loc,z,leila
cp,1,volt,all
*get,n1,node,,num,min
f,n1,amps,strom

nsle
nsel,r,loc,z,0
d,all,volt,0
allsel

/solu
allsel
solve

/prep7
esel,s,type,,2
et,2,solid97
emodif,all,type,2
allsel

/solu
ldread,js
solve

!/contour,,9
  

/post1
esel,s,mat,,1
PLVECT,H, , , ,VECT,ELEM,ON,0