!Homogenes Magnetfeld
!isoliert

finish  
/clear

/filname,test
/title,3-D Feld
!-----------------Variablen----------
abx=0.1			
aby=0.2
abz=0.1

Bx=0
By=0
Bz=0.001

versatz=0

prfrad=0.07
prfdi=0.005
prfho=0.005
aw3=prfdi*prfho	!Wicklungsfläche Prüfling
pi=4*atan(1)
kappa=56e6

abm=0.15

!-----------Modell-------
/prep7
emunit,mks
et,1,solid97
et,2,solid97,1
!et,3,circu124,5,0

		
!-----Material------
mp,murx,1,1      	!Luft
mp,murx,2,1		!Leiter
mp,rsvx,2,1/kappa	!Leiter
!------Geometrie Luft--------------
/view,,1,2,3

!-------Püfling---------------

cylind,prfrad-prfdi/2,prfrad+prfdi/2,-prfho/2,prfho/2,5,355


!-------Luft------------------
cylind,0,abx,-(versatz+abz)/2,(abz+versatz)/2,  0, 90
cylind,0,abx,-(versatz+abz)/2,(abz+versatz)/2, 90,180
cylind,0,abx,-(versatz+abz)/2,(abz+versatz)/2,180,270
cylind,0,abx,-(versatz+abz)/2,(abz+versatz)/2,270,360
bopt,numb,off
vovlap,all
numcmp,volu




!------Vernetzung-----------------
!---Pruef
type,2
mat,2
!real,1
vmesh,1,4

!---Luft
type,1
mat,1
mshkey,0
mshape,1,3d
esize,abx/4
vmesh,5,8
eplot


!------Randbedingung--------
esel,s,mat,,1
csys,1
nsel,s,loc,x,abx
csys,0
*get,nanz,node,,count
*get,nmin,node,,num,min
narb=nmin
*do,i,1,nanz
d,narb,ax,-bz/2*ny(narb)
d,narb,ay,bz/2*nx(narb)
d,narb,az,0
narb=ndnext(narb)
*enddo
nsel,s,loc,z,-(versatz+abz)/2
nsel,a,loc,z,(versatz+abz)/2
d,all,az,0

allsel

!---Sekundärspule-----
!esel,s,type,,2
!nsle,s
!cp,1,emf,all
!cp,2,curr,all

allsel

!csys,1
!esel,s,mat,,2
!nsel,s,loc,y,5
!*get,n1,node,,num,min
!Csys,0

!allsel
!*get,knr,node,,num,max
!*get,emax,elem,,num,max
!n,knr+1,-2*abm,0
!n,knr+2,-2*abm,abm

!r,2,1
!type,3
!real,2
!e,knr+2,knr+1,n1


!------Schaltungsmasse-----
!d,knr+1,volt,0

allsel

!-----Loesung----------
/solu
antype,3
harfrq,50
solve

/post1
set,1,1
plvect,b

esel,s,mat,,2
plns,volt