!aus dem Forum von herbstritt

finish
/clear
*abbr,herbstritt2,/input,herbstritt2.txt

! 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
mp,murx,2,1        !Kupfer

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

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

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

! Leiter
rectng,radin,radau,-h/2,h/2

! Meshverfeinerung
rectng,,radau*1.2,-3*h,+3*h             
aovlap,all

! A5: Luft außen, A3: Leiter, A7: Meshverfeinerung, A6: 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,3
! /eof

! Meshverfeinerung
mshape,0,2d
mshkey,
esize,2e-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,6
! /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
esel,s,mat,,2
bfe,all,js,,,,stromd



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

/solu
allsel
solve


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

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

 /eof

!----zur Kontrolle----------
*get,b_mitte,node,node(0,0,0),b,y

!----Pfad----
path,achse,2,200,200
ppath,1,,0,-beradin
ppath,2,,0,+beradin
!/axlab,x,z-richtung
pdef,by,b,y
plpath,by

/pbc,path,1
plns,b,sum


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