finish
/cle
! titel
/TITLE,Leiterschleife
*abbr,schleife,/input,leiterschleife,inp


/PREP7    ! Start Preprocessor-Phase 

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! element definition
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

et,1,plane53,,,1   ! Luft
et,2,plane53,,,1   ! Leiter
et,3,infin110,,,1  ! Luft außen


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! material properties
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

mp,murx,1,1        ! Luft 
mp,murx,2,1.07     ! NdFeB


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Geometrie Parameter
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
radin=0           ! Leiterschleifenradius innen
radau=1.5e-3      ! Leiterschleifenradius außen 
h=1e-3            ! Leiterschleife Höhe
beradin=19e-3     ! Bereichsradius innen
beradau=20e-3     ! Bereichsradius außen
stromd=8.83694e8  ! Stromdichte in Leiterschleife
strom=8.83694e5   ! Strom in Leiterschleife

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Geometrie
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! Luft
pcirc,0,beradau,-90,90
pcirc,0,beradin,-90,90

! Magnet
rectng,0,radau,0,h

! Meshverfeinerung
rectng,,radau*1.2,0,h             
aovlap,all

! A5: Luft außen, A6: Leiter, A7: Meshverfeinerung, A8: Luft innen
! /eof


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Mesh
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! Leiter
mshape,0,2d    ! Vernetzen mit Vierecken
mshkey,1       ! mapped Mesh
esize,1e-4     ! Maschenweite
type,2  $mat,2 !$real,2
amesh,6
! /eof

! Meshverfeinerung
mshape,0,2d
mshkey,1
esize,1e-4
type,1   $mat,1
amesh,7         
! /eof


! Luft außen
mshape,0,2d
mshkey,1
esize,1e-3
type,1   $mat,1
amesh,5
! /eof

! Luft innen
mshape,0,2d
mshkey,0
esize,1e-4
type,1  $mat,1
amesh,8
! /eof

! Rand
csys,1
nsel,s,loc,x,beradau
esln
emodif,all,type,3
sf,all,inf
csys,0
! /eof

! Strom
! lsel,s,line,,8
! nsll
! cp,1,volt,all
! *get,n1,node,,num,min
! f,n1,amps,strom
! /eof

! Stromdichte
lsel,s,line,,8
nsll
cp,1,volt,all
bf,all,js,,,stromd
! /eof


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Solution
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

/solu
allsel
solve


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Postprocessing
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

/post1
! /expand,2,polar,half,,90
! plf2d
! /eof

path,hoch,2
ppath,1,,0,h
ppath,2,,0,h+3e-3
/axlab,x,z-richtung
pdef,bsum,b,sum
plpath,bsum
paget,achse,table
/eof