!Eine Wicklung mit Kern
!harmonisch
!axialsymmetrisch
!Vorab-Übung für Spartrafo
!Schätzing, 18.10.2008
finish  $/clear
*abbr,vorab1,/inp,vorab1,txt

/prep7
et,1,plane53,,,1         !Luft, Kern
et,2,plane53,3,,1        !Wicklung
et,3,infin110,,1,1       !Rand
et,4,circu124,0          !Widerstand
et,5,circu124,4,1        !Spannungsquelle,sinus
et,6,circu124,5          !Kopplung zur Wicklung

mp,murx,1,1         !Luft
mp,murx,2,1         !Wicklung1
mp,rsvx,2,1/56e6    !Wicklung1
mp,murx,4,200       !Kern
!--------------------------
rak=10e-3       !Radius außen Kern 
hk=50e-3        !Höhe Kern
riw=11e-3       !Radius innen Wicklung
raw=14e-3       !Radius außen Wicklung
hw=40e-3        !Höhe Wicklung
nw=1000         !Windungszahl Wicklung
!=================

aw=(raw-riw)*hw  !Fläche Wicklung
!Real für Wicklung:care,turn,leng,dirz,fill
r,1,aw,nw,,1,0.8
!-------------------------
rba=70e-3  !Radius Bereich außen
rbi=65e-3  !Radius Bereich innen

!-----Geometrie------
/color,outl,blac

!-----Kern------
rectng,0,rak,-hk/2,hk/2
!-----Wicklung-----
rectng,riw,raw,-hw/2,+hw/2
!----Luft------
pcirc,0,rba,-90,90
pcirc,0,rbi,-90,90
!-------------------
aovlap,all
!A2:Wi A6:Kern A5:Rand A7:Luft

!/eof

!-----Vernetzung-----
type,1
mat,4
esize,rak/4
amesh,6   !Kern
!--------------------
type,2
mat,2
esize,(raw-riw)/3
amesh,2   !Wicklung

!------Vernetzung-Luft------
type,1
mat,1
csys,1
lsel,s,loc,x,(rba+rbi)/2
lesize,all,,,1
esize,,36
amesh,5
nsel,s,loc,x,rba  !Randknoten
esln
emodif,all,type,3
sf,all,inf        !Infiniti-Flag
csys,0

esize
!mshkey,0
!mshape,1
amesh,7    !Luft innen

!-----Kopplung------
esel,s,mat,,2
nsle
cp,1,curr,all
cp,2,emf,all

allsel
!-----Speiseschaltung-----
!-----Knoten------
abst=20e-3  !für Knoten
*get,n_max,node,,num,max   !höchste Knotennummer
*get,e_max,elem,,num,max   !höchste Elementnummer
n,n_max+1,-abst*2,-abst
n,n_max+2,-abst*2,abst

type,6   !Koppelelement
r,6,1
real,6
e,n_max+1,n_max+2,node(riw,-hw/2,0)
!------------------

!-------Stromspeisung-------
f,n_max+2,amps,5
!-------Spannungsspeisung----
!d,n_max+2,volt,129.58
!-------Knoten auf Masse---------
d,n_max+1,volt,0

allsel
!----------------
/solu
antype,harm
harfre,50
solve
!-------------
/post1
!-------Realteil----
set,1,1
plf2d
!esel,s,,,e_wider
!presol,elem
ug_real=volt(n_max+2)

!-------Imaginärteil-----
set,1,1,,1
plf2d
ug_imag=volt(n_max+2)
!------Betragsbildung-------
ug_betr=sqrt(ug_real**2+ug_imag**2)