!*********************************************
!** Calcuation of selfind conductor         **
!**					    **
!** Peters 21.2.2006			    **
!*********************************************

fini 
/clear 

*abbr,sefind,/input,selfind,txt

!**Define constants
freq=50
strom=10e-3
kappa=56e6


/prep7 
/view,,1,2,3

bloc,-50e-5,50e-5,-25e-5,25e-5,-0.75e-2,0.75e-2 !Leiter
bloc,-5e-3,5e-3,-5e-3,5e-3,-1e-2,1e-2 !Luft


vovl,all
!V1:Leiter  V2:luft

!**************
et,1,117     !Luft 
et,2,117,1   !Leiter

mp,murx,1,1 ! Mur Luft ist 1

mp,murx,2,1 ! Mur Metallblock
mp,rsvx,2,1/kappa !widerstand 1/leitfähigkeit


!**Leiter
type,2
mat,2 ! wähle mat 1
vmesh,1 ! vermasche mat 1


!**Luft
type,1
mat,1
esize,0.001
mshape,1,3D ! Tetraheder Netz
vmesh,3


!*************************
/solu
!**Magnatic BCs*********
alls
nsel,s,ext !Nodes on exterior of selected elements (ignore remaining fields)
d,all,az,0

!**BCs**************
esel,s,mat,,2
nsle
nsel,r,loc,z,-0.75e-2
cp,1,volt,all
*get,n1,node,,num,min
f,n1,amps,strom ! Strom in Leiter in z-Richtung
nsle
nsel,r,loc,z,0.75e-2
d,all,volt,0
allsel
antyp,harm
harfr,freq
solve


!*************************
/post1


set,1,1

allsel
etable,energie,sene
ssum
*get,sumenergie,ssum,,item,energie
ind_H=2*sumenergie/strom/strom
ind_nH=ind_H*1e9


/eof