!Dilatometer 2D
!Probenerwärmung
!harmonisch
!-------------------------------------------------
finish  $/clear
*abbr,dilatometer,/input,dilatometer,inp

/config,nproc,2
!------PARAMETERFESTLEGUNG------------------------
!------Geometrie----------------------------------
radpro=3e-3     !Radius Probe
lanpro=12e-3    !Länge Probe
radba=10e-3     !Radius Basis
langes=60e-3    !Länge gesamt
!-----Wicklung------
radwin=16e-3    !Radius Windung
raddra=2e-3     !Radius Draht
stei=raddra*3   !Steigung der Wicklung

radau=60e-3	!Radius außen, Luftraum
strom=200
freq=1000

!------Sonstige Kenngrößen------------------------
wikappa=56e6		!spez. el. Leitfähigkeit der Wicklungen
prokappa=10e6           !Probe
!------Elementetypen------------------------------
/prep7			!Aufruf des Preprozessors
et,1,plane53,0		!Luft
et,2,plane53,1		!Windung
et,3,plane53,0          !Probe
et,4,plane53,0          !Basis
!------Materialeigenschaften------------------------
mp,murx,1,1       	!Luft
mp,murx,2,1		!Wicklung
mp,rsvx,2,1/wikappa	!spez. el. Leitfähigkeit (Wicklung)
mp,murx,3,1000
mp,rsvx,3,1/prokappa
!------Real-Konstanten---------------------------------
!------Geometrie------------------
!------Probe---------
rectng,0,radpro,-langes/2,langes/2
rectng,0,radba,lanpro/2,langes/2
rectng,0,radba,-lanpro/2,-langes/2

!------Wicklung--------
local,11,0,radwin,stei*0-stei/2
local,12,0,radwin,stei*1-stei/2
local,13,0,radwin,stei*2-stei/2
local,14,0,radwin,stei*3-stei/2
local,15,0,radwin,stei*(-1)-stei/2
local,16,0,radwin,stei*(-2)-stei/2

wpcsys,,11  $pcirc,raddra
wpcsys,,12  $pcirc,raddra
wpcsys,,13  $pcirc,raddra
wpcsys,,14  $pcirc,raddra
wpcsys,,15  $pcirc,raddra
wpcsys,,16  $pcirc,raddra

wpcsys,,0
csys,0
!------Luft--------
pcirc,,radau,-90,90  

aovlap,all             
numcmp,area  
/pnum,area,1        
!A11:Probe, A7-10:Basis
!A1-A6:Windung,  A12:Umgebung

/eof

!-------Vernetzung Probe-------------
esize,radpro/4  	!Maschenweite
type,3	$mat,3	  	!Probe
amesh,11

!-------Vernetzung Basis-------------
esize,radpro/4  	!Maschenweite
type,4	$mat,1	  	!Basis
amesh,7,10


!-------Vernetzung Windung-------------
esize,raddra/4  	!Maschenweite
type,2	$mat,2	  	!Windung
amesh,1,6

!------Vernetzung der Umgebung-------------
esize
!mshkey,0		
!mshape,1,3D		
type,1	$mat,1	
amesh,12
	
!------Rand----------------------
nsel,s,ext
d,all,az,0


!------Lasten-----
asel,s,,,1
nsla,s,1
cp,1,volt,all
*get,n_1,node,0,num,max
f,n_1,amps,strom

asel,s,,,2
nsla,s,1
cp,2,volt,all
*get,n_2,node,0,num,max
f,n_2,amps,strom

asel,s,,,3
nsla,s,1
cp,3,volt,all
*get,n_3,node,0,num,max
f,n_3,amps,strom

asel,s,,,4
nsla,s,1
cp,4,volt,all
*get,n_4,node,0,num,max
f,n_4,amps,strom

asel,s,,,5
nsla,s,1
cp,5,volt,all
*get,n_5,node,0,num,max
f,n_5,amps,strom

asel,s,,,6
nsla,s,1
cp,6,volt,all
*get,n_6,node,0,num,max
f,n_6,amps,strom



!------LÖSUNG----------------------------------
/solu			
allsel			
antype,harm	
harfre,freq
solve

!-------------------
/post1
/plopts,minm,0
esel,s,mat,,3
etable,verlu,jheat
pletab,verlu,avg

/eof

set,1,1          !Realteil
esel,s,mat,,3
etable,strdi,jt,sum
pletab,strdi,avg      !Stromdichte

!-----3D-Darstellung-geschnitten-------
/expand,27,axis,,,10
pletab,verlu,avg