!leitendes Rohr im vorher homogenen Feld
!Schirmwirkung durch my-Rohr
!Schirmwirkung durch Wirbelströme
!kartesisch 3D
!

!Schätzing  1.12.2009

finish  $/clear
*abbr,rohr1,/input,rohr1,txt

!------Parameter--------
abx=0.3    !Abmessung in x-Richtung 
aby=0.3    !Abmessung in y-Richtung 
abz=0.3    !Abmessung in z-Richtung
!-------Rohr------
radin=40e-3  !Radius innen
radau=44e-3  !Radius außen
lang=120e-3  !Länge

Bz=1   !Flussdichte in z-Richtung

!------Modellerstellung-----------------
/prep7  
et,1,solid97       !Luft
et,2,solid97,1     !Rohr

mp,murx,1,1        !Luft
!------Rohr-------
mp,murx,2,1
mp,murx,2,1000.0   !Eisenwerkstoff
mp,rsvx,2,1/56e6    !Kupfer
mp,rsvx,2,1/7e6     !Eisenwerkstoff

!----Geometrie-----
/view,,1,2,3
block,-abx/2,    0,-aby/2,aby/2,-abz/2,abz/2
block,     0,abx/2,-aby/2,aby/2,-abz/2,abz/2

drehwink=30   !Verdrehwinkel für Rohr
local,11,0,,,,,drehwink
wpcsys,11
cylind,0,radin,-lang/2,lang/2
cylind,0,radau,-lang/2,lang/2,  0, 90
cylind,0,radau,-lang/2,lang/2, 90,180
cylind,0,radau,-lang/2,lang/2,180,270
cylind,0,radau,-lang/2,lang/2,270,360

vovlap,all
!V12,13,14,15:Rohr     !V8,9,10,11:innen
!V16,17:außen
!/eof

!----Vernetzung----------------
type,2  $mat,2    !Rohr
esize,(radau-radin)*3
vmesh,12,15

/eof

type,1  $mat,1    !Luft innen
mshkey,0
mshape,1
!esize,radin/6
vmesh,8,11


esize
vmesh,all
eplot

!----Randbedingungen----
nsel,s,ext
*get,nanz,node,,count
*get,nmin,node,,num,min
narb=nmin

*do,i,1,nanz
  d,narb,ax,1*(-bz/2)*ny(narb),0  
  d,narb,ay,1*(+bz/2)*nx(narb),0  
 narb=ndnext(narb)
*enddo

!-----------------------------
allsel
/solu
antyp,harm
harfre,50
!harfre,0.1
solve
!----------------------
/post1
nsel,s,loc,x,-abx/2,0
esln,s,1
/plopts,minm,0
set,1,1
/view,,1
/cont,,9,0,,1
plns,b,sum

/eof

!------------
/cont,,9
ples,jt,x   !Stromdichte
set,1,1,,1  !Imaginärteil
ples,jt,x