!#############################################
!---------Ring aus Keypoints------------------
!#############################################

!Zylinder-Spule
!Schätzing	16.02.2001
!Görnandt	29.03.2001
!Schätzing  	3.7.2008

finish
/clear
*abbr,spirale,/inp,spirale_3d,txt

!------Kenngrößen-----------------------
rad=15e-3		!Radius der Spule innen
radi=11.75e-3
leila=0.03
wdg=1			!Spiralwindungen
winkel=360*(wdg-1)+360	!Spiralelänge in Grad
pi=3.14159265359
strom=sqrt(2)*800	!Strom (Scheitelwert)
freq=50	

!------Modellerstellung---------------
/prep7
et,1,solid117		!Element Luft
et,2,solid117,6		!Element Alustab
et,3,solid117,1		!Element Stahlring
mp,murx,1,1		!mu_relativ Luft
mp,murx,2,1		!mu_relativ Alustab
mp,murx,3,100		!mu_relativ Ring
mp,rsvx,2,1/28.64e6	!spezifischer el. Widerstand Alustab
mp,rsvx,3,1/5e6		!spezifischer el. Widerstand Ring


!------SPULE/SPIRALE--------------------------
!------Keypoints und Linie für Spirale---------
csys,1
vteile=winkel/10	!alle 10 Grad ein Keypoint
*do,i,0,vteile
!k,i+1,rad,10*i,0.001*wdg*i/vteile
k,i+1,rad,10*i,0
*enddo
*do,j,1,vteile
l,j,j+1
*enddo
lsel,s,,,1,35,1
lcomb,all
lesize,all,,,30

!------Querschnittsfläche des Leiters-----------
csys,1
K,500,rad,0.001,0
circle,1,2.013e-3,500
lsel,s,,,2,5,1
lesize,all,,,3
al,2,3,4,5

!------Volumenelemente der Spirale---------------
csys,1
vdrag,1,,,,,,1

!-------Alustab--------------------
cylind,,radi,-0.005,leila+0.005,0,90		!Nachbildung Alustab 
cylind,,radi,-0.005,leila+0.005,90,180		!Viertelzylinder, damit 
cylind,,radi,-0.005,leila+0.005,180,270		!Vernetzung mit Quadern 
cylind,,radi,-0.005,leila+0.005,270,360
vglue,2,3,4,5
numcmp,volu

!--------Lufthülle----------------
csys,0
cylind,,rad*2,-0.01,leila+0.01			!Luft

vovlap,all					!Überlappen der Volumen
numcmp,volu					!Komprimierung der Nummerierung


!-------Vernetzung des Alustabes----------------
type,2
mat,2
lsel,s,,,20,21,1
lsel,a,,,29
lsel,a,,,38
lesize,all,,,20
type,2
mat,2
esize,,4					
mshape,0,3d				!Vernetzung mit Quadern
mshkey,1				!Mapped mashing
vmesh,2,5,1

!------Vernetzung des Ringes----------------------
type,3
mat,3
esize,,4
mshape,1,3d
mshkey,5
vmesh,1

!-------Vernetzung der Luft----------------
type,1
mat,1
mshkey,0
mshape,1,3d
smrtsize,3
vmesh,6

nummrg,node
nummrg,kp

!------Randbedingungen----------------
csys,1					!Koordinatensystem
asel,s,loc,x,2*rad			!Ränder der Luft 
da,all,az,0				!mit Potential 0 belegen
asel,s,loc,z,-0.01
asel,a,loc,z,leila+0.01
da,all,az,0

!-----Last-----------------------------
esel,s,mat,,2				!Last auf Alustab aufbringen
nsle
nsel,r,loc,z,-0.005
cp,1,volt,all
*get,n1,node,,num,min
f,n1,amps,strom				!Strom einspeisen

esel,s,mat,,2
nsle
nsel,r,loc,z,leila+0.005
d,all,volt,0				!Potential Null

/solu					!Lösung
allsel
antype,harm				!harmonisch
harfrq,freq				!Frequenzvorgabe			
solve


/post1
set,1,1
esel,s,mat,,2
powerh					!Leistung Alustab

esel,s,mat,,3				!Leistung Ring
powerh