!Dilatometer 2D
!Probenerwärmung
!elektrisch harmonisch
!thermisch transient
!-------------------------------------------------
finish  $/clear
*abbr,dil_harm_trans2,/input,dil_harm_trans2,txt

/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 um später die koordinaten der Kreise zu setzen

radau=60e-3	!Radius außen, Luftraum
strom=500
freq=1000
!---------Definitionen für Schleife und UNtersuchungszeit-----
endzeit=5   !Endzeit
sw=0.1      !Zeitschritt für transiente Analyse
zeit=1e-9   !Startzeit

!------Sonstige Kenngrößen------------------------
wikappa=56e6		!spez. el. Leitfähigkeit der Wicklungen in S/m
prokappa=10e6  !spez. el. Leitfähigkeit der Probe in S/m
!------Elementetypen------------------------------
/prep7			!Aufruf des Preprozessors
et,1,plane53,0,,1		!Luft      ELEMENTYPE,NUMMER,ELEMENTTYP AUS AUS NAME,KEYOPT(1)=0,keyopt(2) ...
                 !           0 bedeutet AZ degree of freedom: static domain, induced eddy current domain
et,2,plane53,1,,1		!Windung    1 bedeutet AZ, VOLT degrees of freedom: current-fed massive conductor
                !(Stromgespeiste Windung)
et,3,plane53,0,,1          !Probe
et,4,plane53,0,,1          !Basis         Stempel weggelasse?
!------Materialeigenschaften------------------------
mp,murx,1,1                 	!Luft MP, Lab, MAT, C0, C1, C2, C3, C4
!                            Lab = MURX	 — 	MAGNETIC RELATIVE PERMEABILITIES (also MURY, MURZ).
!                            MAT = verweist auf das element hier "1" also DIE LUFT
                            !C0= Zahlenwert für das Material hier LUFT  = 1
mp,murx,2,1		!Wicklung
mp,rsvx,2,1/wikappa	  !spez. el. Leitfähigkeit (Wicklung)   RSVX = ELECT. RESISTIV ELEKTRISCHER WIDERSTAND in X = 1/leitf.
mp,murx,3,1000       !Die Probe hat eine relative Permeabilität von 1000
mp,rsvx,3,1/prokappa
!------Real-Konstanten------------------------------
!------Geometrie------------------
!------Probe---------
rectng,0,radpro,-langes/2,langes/2 !Koordinatensystem wird in die mitte der Probe gesetzt RECTNG X1, X2, Y1, Y2
!                                   zwei Punkte die nicht auf einer Linie liegen  hier hat er es fasch verstanden
                                    !denn Probe ist nur rectng,0,radpro,-lanpro/2,lanpro/2
rectng,0,radba,lanpro/2,langes/2        !Stempel Flächen oben
rectng,0,radba,-lanpro/2,-langes/2      !Stempel Flächen oben

!------Wicklung--------
local,11,0,radwin,stei*0-stei/2      !LOCAL, KCN=beliebige referenznummer ab 11, KCS=0 bedeutet Kartesisch, XC, YC,
local,12,0,radwin,stei*1-stei/2      !hier werden nur für die windungskreise in der ebene die bezogenen Koordinatensysteme definiert
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      !in den Bezugskoordinatensystemen wird ein kreis definiert mit dem Drahtradius
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  !komprimiert Flächennummern
/pnum,area,1 !/PNUM, Label, KEY=turns on numbering
!A11:Probe, A7-10:Basis
!A1-A6:Windung,  A12:Umgebung

!/eof   ! Exits the file being read.


!-------Vernetzung Probe-------------
esize,radpro/4  	!Maschenweite
type,3	$mat,3	  	!Probe
amesh,11          !amesh= generate nodes and area elem within areas area 11

!-------Vernetzung Basis-------------
esize,radpro/4  	!Maschenweite
type,4	$mat,1	  	!Basis
amesh,7,10        !area 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

!------Randbedingungen-elektrisch---------------------
nsel,s,ext                !Randknoten
d,all,az,0                !Vektorpotential auf Null


!------Lasten-auf-Wicklung---
asel,s,,,1 !select area 1
nsla,s,1   !select all nodes for the area above
           !1= Select all nodes (interior to area, interior to lines, and at keypoints)
cp,1,volt,all   !CP, NSET=nur eine nummer, Lab=nur ein name, NODE1 =welche knoten
*get,n_1,node,0,num,max !NUM	MAX, MIN	Highest or lowest node number in the selected set.
f,n_1,amps,strom    !F, NODE=welche Knoten, Lab=force label muss zulässig sein!!, VALUE=Wert hier eine variabel

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
!--Darstellung der Verlustleistung---
etable,verlu,jheat  !ETABLE, Lab=vorher nicht benutzter name, Item, Comp
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


!/eof


/prep7
lsclear,all !Lasten löschen
!--------TRANSIENTE BERECHNUNG-----
!----------------------------------
!--------------ELEMENTDEF----------
et,1,0   !Nullelement
et,2,0           !2=temp. deg free 1= axialsymm
et,3,plane77,,,1
et,4,0

!------MATERIALEIGENSCHAFTEN--------
!mp,kxx,2,60 !Wärmeleitfähigkeit der Windung  KUPFER
!mp,c,2,460,380  !spe. w.kapazität
!mp,dens,2,8900  !Dichte
mp,kxx,3,58 !Wärmeleitfähigkeit der Stahlprobe
mp,c,3,460
mp,dens,3,7860

!-----Randbedingung thermisch--------------
tumg=20    !Umgebungstemperatur
nsel,s,loc,x,radpro
nsel,r,loc,y,-lanpro/2,lanpro/2
sf,all,conv,10,tumg   !Konvektion
sf,all,conv,100,tumg
allsel

!-----LOADS from Harmonic Solution------And new solution----
/solu
antype,trans    !transiente lösung
tunif,20        !Umgebungstemperatur "Assigns a uniform temperature to all nodes."
outres,all,all
esel,s,mat,,3       !2 Kupfer und 3=Stahl !Lücke nicht notwendig wenn mat ausgewählt
ldread,hgen,,,,2,,rst   !Übernahme der Verlustleistung
endzeit=100
time,endzeit
deltim,endzeit/10
solve

!-----Auswertung-------
/post1
plns,temp

/eof

!---------------
/post26
n_mitte=node(0,0,0)
n_oben=node(0,lanpro/2,0)
n_rand=node(radpro,0,0)
nsol,2,n_mitte,temp,,Mitte
nsol,3,n_oben,temp,,oben
nsol,4,n_rand,temp,,Rand

plvar,2,3,4