!*********************************************
!** Berechnung von B-Feld einen             **
!** Vierkanten-Induktor                     **
!** mit Ferrotronkern                       **
!**					    **
!** Segbeaya 16.06.2006			    **
!*********************************************

!fini 
!/clear
!/BATCH  
!/COM,ANSYS RELEASE  5.7.1  UP20010418       12:06:56    06/16/2006
!/input,menust,tmp     ,'',,,,,,,,,,,,,,,,1  
!/GRA,POWER
!/GST,ON
!/PLO,INFO,3
!/COL,PBAK,ON,1,BLUE 



!!
!!**Parameter deklarieren**
!!

!hs: hoehe  der spule
!bs: breite der spule
!ds: dicke  der spule
!hm: hoehe  der Magnetkern
!bm: breite der magnetkern
!bl: breite luft
!win: Anzahl windungen
!ff: fuellfaktor
!ed: eindringtief für stromfluss
!bk: breite leiter kupfer zum extr
!hk: hoehe  leiter kupfer
!gap:Anschlußöffnung
!tm: tiefe magnetkern
!mw: maschenweite

!!
!!Variablen definieren
!!

!Geometrie der Induktor (stirn seite)

*SET,bs,0.07 !alt 0.10
*SET,hs,0.03
*SET,hm,0.05
*SET,bm,bs  
*SET,ds,0.01
*SET,bl,0.33
*SET,ed,0.001
*SET,hk,0.005
*SET,bk,0.005    !in zrichtung
*SET,tm,0.01
*SET,mw,0.001
*SET,Lspule,1000e-9 !2231e-11 !induktivität spule
*SET,Rspule,2e-6   !0.2*(hs+bs) !wiederstand spule
*SET,gap,0.001

!!Stromflussflaecheberechnung

flastr=hk*bk          !Flaeche durchstromende strom
lans=2*(bs+hs)-(2*gap)       !Lanege spule
!fla2=(hk-ed)*(bk-ed) !Flaeche Innen
!flastr=fla1-fla2     !Flaeche durchstromende strom




!Wechselstrom

*SET,strom,40e3 
*SET,freq,4e3

strdi=strom/flastr !stromdichte
!--------------------
!Knoten definieren

/PREP7  
K,1,bm/2,hm/2,, 
K,2,-bm/2,hm/2,,
K,3,-bm/2,-hm/2,,   
K,4,bm/2,-hm/2,,
K,5,bs/2,hs/2,, 
K,6,-bs/2,hs/2,,
K,7,-bs/2,-hs/2,,   
K,8,bs/2,-hs/2,,
K,9,bs/2-ds,hs/2-ds,,   
K,10,-bs/2+ds,hs/2-ds,, 
K,11,-bs/2+ds,-hs/2+ds,,
K,12,+bs/2-ds,-hs/2+ds,,


!!!Gap keypoint  
 
K,13,-gap,-hm/2,,   
K,14,-gap,-(hs)/2,,  
K,15,-gap,-hs/2+ds,,  

K,16,gap,-hm/2,,   
K,17,gap,-(hs)/2,,  
K,18,gap,-hs/2+ds,,  

!Linien defineren (Ferrotron559 aussen):1
LSTR,       1,       2  
LSTR,       2,       6  
LSTR,       6,       5  
LSTR,       5,       1  
LSTR,       8,       7  
LSTR,       7,       3  
LSTR,       3,       4  
LSTR,       4,       8  

!Linien defineren (spule):2

LSTR,       6,       7  
LSTR,       5,       8  
LSTR,       9,      10  
LSTR,      10,      11  
LSTR,      11,      12  
LSTR,      12,       9  
LSTR,       5,       9  
LSTR,       6,      10  
LSTR,      11,       7  
LSTR,      12,       8  



!Linien definieren gap

LSTR,       13,      14  
LSTR,       14,     15  
LSTR,       15,     18  
LSTR,       18,     17  
LSTR,      17,      14  
LSTR,      17,      16 
LSTR,       16,     13 
LSTR,       16,     14
LSTR,       14,     18
  
!-----------------
!Area definieren (Ferrotron559 aussen):1
al,1,2,3,4
al,5,6,7,8

!Area definieren (spule):2

al,15,3,16,11
al,16,9,17,12
al,13,17,5,18
al,18,10,15,14

!Area definieren (ferrotron559 innen):3
al,14,11,12,13

!Area definieren (luft):4 und overlap 


!Area definieren (gap)
  
AL,20,21,22,23 

AL,19,23,24,25

 

allsel,all
aovlap,all
!aglue,all

allsel,all
/pnum,area,1 
aplot
!/eof

!!spule
voffst,4,bk,,   
voffst,3,bk,, 
voffst,6,-bk,, 
voffst,13,bk,, 
voffst,12,bk,, 


!!Volume gap

voffst,8,bk,, 
voffst,9,bk,, 

allsel,all
!!!Volume magnet
BLC5,,,bm,hm,-tm  
ALLSEL,ALL


!Volum Luft 
CYLIND,0,bl/6,0,-bl/12,0,360, 
vovlap,all
allsel,all
/pnum,volu,1 
vplot


!---------
!!
!!Elemente definieren
!!
ET,1,SOLID97	!Ferrotron aussen

ET,2,SOLID97	!spule
KEYOPT,2,1,4    !massive conductor
R,2,flastr      !bei care

ET,3,SOLID97	!Luft

ET,4,INFIN111   !Randelment
KEYOPT,4,1,1

!/eof
!---Matreialeigenschaften definieren

!Ferrotron559 aussen:1
MPTEMP,1,0  
MPDATA,MURX,1,,20     !rel. Permeabilität:20  
MPDATA,RSVX,1,,10e3   !spez.wiederstand: 10e3

!Spulen(kupfer):2
MPTEMP,1,0  
MPDATA,MURX,2,,1      !rel. Permeabilität:1  
MPDATA,RSVX,2,,1/56e6 !spez.wiederstand:1/56e6 

!Ferrotron559 innen:3
!MPTEMP,1,0  
!MPDATA,MURX,3,,20     !rel. Permeabilität:20  
!MPDATA,RSVX,3,,10e3   !spez.wiederstand: 10e3

!Luft herum:4
MPTEMP,1,0  
MPDATA,MURX,4,,1    !rel. Permeabilität:1  
MPDATA,RSVX,4,,1e9  !spez.wiederstand: 1e9


!/eof
!-----------
!!Vernetzung und	
!!
!!Zuweisen von Elementen und Materialien auf:
!!

!!Ferrotron magnetkern V17

esize,mw*3
mshkey,0  
mshape,1,3d !tetraeder
type,1   $mat,1  
vmesh,17

!!Spule
esize,mw*3

mshkey,0  
mshape,1,3d     !tetraeder
vsel,s,,,10,11,1
vsel,a,,,1
vsel,a,,,14
vsel,a,,,15
!allsel
type,2   $mat,2  $real,2
vmesh,all

!allsel,all

!!gap luft
esize,mw*2
mshkey,0  
mshape,1,3d !tetraeder
vsel,s,,,12,13,1
type,3   $mat,3  
vmesh,all

!!Luft V16

csys,1
asel,s,loc,x,bl/6
aesize,all,mw*6
allsel
type,3   $mat,3  
mshkey,0
mshape,1,3d
vmesh,16


!!
!!Randelemente
!!

csys,1                    !zylinderkoordinaten
type,4                    !infini
esize,,1                  !ndivision
asel,s,loc,x,bl/6
*get,anr1,area,,num,min
*get,anr2,area,,num,max
vext,anr1,anr2,,mw         !extrusion
nsel,s,loc,x,bl/6+mw/2,bl/6+mw
sf,all,inf                 !flag setzen

csys,0                      !zuruckschalten
allsel

!!
!!
!!Randbedinguen
!!
!!
!vordere Leiterqueschnittsflaeche

esel,s,mat,,2
nsle


nsel,r,loc,x,-gap
!allsel,below,area

!nsla,s,1

!nsel,r,loc,z,lans

cp,1,curr,all
cp,2,volt,all

!stromfluss in diese flaeche hinein (-1)
sf,all,mci,-1
*get,knr1,node,0,num,max
allsel
!-----hintere Leiterquerschnittsflaeche--
esel,s,mat,,2
nsle

nsel,r,loc,x,gap
!nsla,s,1


!allsel,below,area       

!nsel,r,loc,z,lans

cp,3,curr,all
cp,4,volt,all

!stromflussaus dieser stirnflaeche heraus(+1)
sf,all,mci,+1
d,all,volt,0
*get,knr2,node,0,num,max

!!Symetrie
nsel,s,loc,x,-gap
dsym,asym,x
nsel,s,loc,x,gap
dsym,asym,x


allsel
/eof
!/eof
!!Ausssenradbedingungen
!!
!!Kopplung Areas-spule mit dof emf und curr
!!

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

!allsel
!lsel,s,,,19,22,1
!DL,all,,az,0,

!allsel

!/eof





asel,s, , ,3
allsel,below,area

cp,1,curr,all
cp,2,emf,all
*get,n1,node,,num,min



asel,s, , ,4
allsel,below,area
cp,3,curr,all
cp,4,emf,all
*get,n2,node,,num,min


asel,s, , ,5
allsel,below,area
cp,5,curr,all
cp,6,emf,all
*get,n3,node,,num,min



asel,s, , ,6
allsel,below,area
cp,7,curr,all
cp,8,emf,all
*get,n4,node,,num,min


!!!festlegung der außenrandparameter
!!
!!Vektorpotenial AZ am Rand auf 0
allsel
lsel,s,,,19,22,1
DL,all,,az,0,

!allsel

!/eof

!!
!!Definition der Speisende schaltung
!!
!!element2 spule (massive conductor)
!!area3

R,11, ,  
N,30000,-0.35,-0.3,0 
N,30001,-0.15,-0.3,0 
RMOD,11,15,0,1   
ET,11,CIRCU124,6,0   
TYPE,11  
REAL,11  
MAT,1   
!*  
E,30000,30001,n1
!*  

!!area4

R,12, ,  
N,30002,0.5E-01,-0.3,0   
RMOD,12,15,0,2   
ET,12,CIRCU124,6,0   
TYPE,12  
REAL,12 
MAT,1   
!*  
E,30001,30002,n2
!*  


!!area5

R,13, ,  
N,30003,0.25,-0.3,0  
RMOD,13,15,0,3   
ET,13,CIRCU124,6,0   
TYPE,13  
REAL,13 
MAT,1   
!*  
E,30002,30003,n3
!*  


!!area6

R,14, ,  
N,30004,0.25,-0.45,0 
RMOD,14,15,0,4   
ET,14,CIRCU124,6,0   
TYPE,14 
REAL,14  
MAT,1   
!*  
E,30003,30004,n4
!*  


!!
!!Induktivitaet 
!!
!!

R,15,Lspule , ,
N,30005,-0.15,-0.45,0
RMOD,15,15,0,5   
ET,15,CIRCU124,1,0   
TYPE,15  
REAL,15  
MAT,1   
!*  
E,30004,30005 
!*  

!!
!!Widerstand
!!
!!

R,16,Rspule ,  
N,30006,-0.35,-0.45,0
RMOD,16,15,0,6   
ET,16,CIRCU124,0,0  
TYPE,16 
REAL,16 
MAT,1   
!*  
E,30005,30006 
!*  

!!
!!Stromquelle
!!
!!

R,17,strom, ,
RMOD,17,15,0,7   
ET,17,CIRCU124,3,0  
TYPE,11 
REAL,17  
MAT,1   
!*  
E,30000,30006,30006 
!*  
*get,n6,elem,0,num,max
allsel

!ja
/eof






*get,n1,node,0,num,min
f,n1,amps,strom

!allsel
cp,1,emf,all
cp,2,curr,all


asel,s,,,4
nsla,s,0
allsel,below,area


*get,n2,node,0,num,min
f,n2,amps,strom
!allsel

cp,3,emf,all
cp,4,curr,all




asel,s,,,5
nsla,s,0
allsel,below,area



*get,n3,node,0,num,min
f,n3,amps,strom
!allsel

cp,5,emf,all
cp,6,curr,all


asel,s,,,6

nsla,s,0
allsel,below,area


*get,n4,node,0,num,min
f,n4,amps,strom

!allsel

cp,7,emf,all
cp,8,curr,all







!Ferrotron aussen(1)
asel,s,,,1,2,1
aatt,1,1,1

!spule(2)
asel,s,,,3,6,1
aatt,2,2,2
!!

!Ferrotron innen(3)
asel,s,,,7
aatt,3,3,3


!Luft herum(4)
asel,s,,,9
aatt,4,4,4

allsel
!/pnum,type,1
!/replot
!---------------

!!
!! Meshing vorbereiten (lange linien ohne luft )
!!
lsel,s,,,1,3,2
lsel,a,,,9,11,1
lsel,a,,,5,7,2
lsel,a,,,13
lesize,all,,,20,,1

!!
!! Meshing vorbereiten (alle kurze linien)
!!
lsel,s,,,2,8,2
lsel,a,,,19,22,1 !luft
lsel,a,,,15,18,1
lsel,a,,,12,14,2
lesize,all,,,10,,1

!!
!! Meshing
!!
allsel
mshape,1,2D
mshkey,2
amesh,all

allsel














/eof

!!!!
Allsel
FINISH  
/SOLU   
!*  
ANTYPE,0
/STATUS,SOLU
SOLVE 
!!!!  
!!!!Mit befehlen fuer feldlinienberechnung
!!!!
FINISH  
/SOLU   
FINISH  
/POST1  
PLF2D,27,0,10,1