!WS Sensor in Seitenansicht auf Stahldraht
!erfordert Rand infin110

fini  $/cle
*abbr,pro_3_4,/inp,pro_3_4.inp
/CWD,'E:\Holgis_bunte_ANSYS_Welt\Ansys_Ergebnisse\Ferrittest\mue1'
! Materialvorgaben:
!------------------
/prep7		  !Aufruf Preprozessor
/units,si

et,1,plane53,0         	!Luft, Leiter
et,2,plane53,3         	!Kupferwicklung durch CIRCU124 gespeist,mit Induktionsspannung
et,3,infin110        	!Randelement, Boundary
et,4,124,3,1            !Stromquelle, Sinus, fuer Speisung
et,5,124,5              !Induktor fuer Speisung
!--------------------------------------------------------------------

mp,murx,1,1       	!Permeabilitaet Luft
mp,rsvx,1,1E8   	!Resistivitaet  Luft
!--------------------------------------------------------------------

mp,murx,2,1      	!Permeabilitaet Kupfer
mp,rsvx,2,1/56E6   	!Resistivität   Kupfer
!--------------------------------------------------------------------

MPTEMP,R5.0, 2, 1,  20.0000000    ,  60.0000000    ,
MPDATA,rsvx, 3, 1,  2.500000000E-07, 1.000000000E-06, !Resistivitaet Stahl

 TB,BH  ,3   ,   1,  45	!BH-Kurve Stahl
 TBTEM,  20.00000000    ,   1
 TBPT,,  50.0000000    , 1.000000000E-03
 TBPT,,  194.000000    , 3.300000000E-03
 TBPT,,  309.000000    , 1.000000000E-02
 TBPT,,  540.000000    , 3.090000000E-02
 TBPT,,  826.000000    , 5.420000000E-02
 TBPT,,  1069.00000    , 0.105100000    
 TBPT,,  1264.00000    , 0.231300000    
 TBPT,,  1439.00000    , 0.427200000    
 TBPT,,  1555.00000    , 0.484700000    
 TBPT,,  1749.00000    , 0.920700000    
 TBPT,,  2090.00000    ,  1.22390000    
 TBPT,,  2381.00000    ,  1.34790000    
 TBPT,,  3402.00000    ,  1.55260000    
 TBPT,,  4422.00000    ,  1.62450000    
 TBPT,,  5199.00000    ,  1.65900000    
 TBPT,,  5783.00000    ,  1.68100000    
 TBPT,,  6366.00000    ,  1.69640000    
 TBPT,,  7629.00000    ,  1.72630000    
 TBPT,,  8309.00000    ,  1.73850000    
 TBPT,,  9087.00000    ,  1.75290000    
 TBPT,,  9719.00000    ,  1.76280000    
 TBPT,,  10448.0000    ,  1.77720000    
 TBPT,,  11420.0000    ,  1.78940000    
 TBPT,,  12003.0000    ,  1.79940000    
 TBPT,,  13364.0000    ,  1.81600000    
 TBPT,,  14287.0000    ,  1.82590000    
 TBPT,,  15259.0000    ,  1.83810000    
 TBPT,,  16134.0000    ,  1.84810000    
 TBPT,,  17154.0000    ,  1.85580000    
 TBPT,,  17689.0000    ,  1.86240000    
 TBPT,,  18515.0000    ,  1.86910000    
 TBPT,,  19049.0000    ,  1.87570000    
 TBPT,,  20313.0000    ,  1.88570000    
 TBPT,,  21722.0000    ,  1.89670000    
 TBPT,,  22256.0000    ,  1.89900000    
 TBPT,,  23957.0000    ,  1.91000000    
 TBPT,,  25366.0000    ,  1.91890000    
 TBPT,,  27067.0000    ,  1.93110000    
 TBPT,,  28817.0000    ,  1.93990000    
 TBPT,,  29692.0000    ,  1.94430000    
 TBPT,,  30858.0000    ,  1.95100000    
 TBPT,,  32607.0000    ,  1.95760000    
 TBPT,,  33822.0000    ,  1.96310000    
 TBPT,,  34405.0000    ,  1.96540000    
 TBPT,,  35037.0000    ,  1.96870000    

!--------------------------------------------------------------------

MPTEMP, 1,  20.0000000    ,  60.0000000    ,
MPDATA,RSVX, 4, 1,  3.000000000E-02, 4.000000000E-02,

 TB,BH  ,4   ,   1,  16!38	!BH-Kurve Ferrit
 TBTEM,  20.0000000    ,   1
 !TBPT,,  1.00000000    , 0.145000000    
 !TBPT,,  10.0000000    , 0.199000000    
 !TBPT,,  13.0000000    , 2.000000000E-03
 !TBPT,,  14.0000000    , 2.200000000E-02
 !TBPT,,  18.0000000    , 4.200000000E-02
 !TBPT,,  20.0000000    , 6.500000000E-02
 !TBPT,,  23.0000000    , 7.500000000E-02
 !TBPT,,  25.0000000    , 9.000000000E-02
 !TBPT,,  26.0000000    , 0.105000000    
 !TBPT,,  28.0000000    , 0.115000000    
 !TBPT,,  31.0000000    , 0.125000000    
 !TBPT,,  33.0000000    , 0.133000000    
 !TBPT,,  35.0000000    , 0.192000000    
 !TBPT,,  38.0000000    , 0.217000000    
 !TBPT,,  40.0000000    , 0.226000000    
 !TBPT,,  43.0000000    , 0.242000000    
 !TBPT,,  46.0000000    , 0.260000000    
 !TBPT,,  48.0000000    , 0.269000000    
 !TBPT,,  50.0000000    , 0.280000000    
 !TBPT,,  53.0000000    , 0.294000000    
 !TBPT,,  56.0000000    , 0.302000000    
 !TBPT,,  60.0000000    , 0.310000000    
 TBPT,,  66.0000000    , 0.323000000    
 TBPT,,  70.0000000    , 0.332000000    
 TBPT,,  73.0000000    , 0.338000000    
 TBPT,,  76.0000000    , 0.346000000    
 TBPT,,  81.0000000    , 0.352000000    
 TBPT,,  85.0000000    , 0.358000000    
 TBPT,,  90.0000000    , 0.365000000    
 TBPT,,  103.000000    , 0.375000000    
 TBPT,,  110.000000    , 0.380000000    
 TBPT,,  132.000000    , 0.391000000    
 TBPT,,  156.000000    , 0.397000000    
 TBPT,,  161.000000    , 0.398000000    
 TBPT,,  190.000000    , 0.401000000    
 TBPT,,  285.000000    , 0.406000000    
 TBPT,,  357.000000    , 0.409000000    
 TBPT,,  396.000000    , 0.410000000    

!---------------------------------------------------------------------------------------
!---- Geometrie erstellen ----

BFE 	 = 2.000000000E-03
BSCH	 = 1.900000000E-03
BSt	 = 5.000000000e-03
BW       = 1.500000000E-03
LFE      = 9.000000000E-03
LSCH     = 2.050000000E-03
LSP      = 5.000000000E-05
LST      = 5.400000000E-02
LW       = 1.800000000E-03
MST      = 5.000000000E-03
SGH      = 1.800000000E-02
DM	 = 5.000000000E-03
ABS	 = 7.500000000E-03
HM	 = 5.000000000E-03

radau=55e-3             !Radius außen
radmi=50e-3      	!Radius innen
winkel=180       	!Winkel des Geometriesegments
!-----------------------!Variablenvereinbarung------------------------------------------
strom=0.092   		!(65mA_eff) Strom im Leiter (Scheitelwert)
N=125.0			!Gesamtwindungszahl pro Schenkel
/prep7		 	!Aufruf Preprozessor

!zeichnen der Geometrie

cyl4,radau,0,radau,0,,winkel				              !Außenkreis							                                !Zeichnen des aeusseren 4tel-Kreises
cyl4,radau,0,radmi,0,,winkel					      !Innenkreis								                                !Zeichnen des inneren 4tel-Kreises
rectng,(Radau*2-LST)/2,(Radau*2-LST)/2+Lst,0,Bst/2						                                !Zeichnen des Stahls innen
rectng,(Radau*2-LST)/2,(Radau*2-LST)/2+Lst,Bst/2,3*Bst/2					                                !Zeichnen des Stahls aussen

rectng,(Radau*2-LST)/2+Lst/2-LFE/2,(Radau*2-LST)/2+Lst/2-LFE/2+BSCH,3*Bst/2+LSP,3*Bst/2+LSCH+LSP    				!Zeichnen des Ferrit
rectng,(Radau*2-LST)/2+Lst/2+LFE/2,(Radau*2-LST)/2+Lst/2+LFE/2-BSCH,3*Bst/2+LSP,3*Bst/2+LSCH+LSP    				!Zeichnen des Ferrit
rectng,(Radau*2-LST)/2+Lst/2-LFE/2,(Radau*2-LST)/2+Lst/2+LFE/2,3*Bst/2+LSP+LSCH,3*Bst/2+LSP+2*LSCH 				!Zeichnen des Ferrit

rectng,(Radau*2-LST)/2+Lst/2-LFE/2-BW-LSP/2,(Radau*2-LST)/2+Lst/2-LFE/2-LSP/2,3*Bst/2+LSP+LSCH-Lw,3*Bst/2+LSP+LSCH    		!Zeichnen der Spulen
rectng,(Radau*2-LST)/2+Lst/2-LFE/2+BSCH+LSP/2,(Radau*2-LST)/2+Lst/2-LFE/2+BSCH+LSP/2+BW,3*Bst/2+LSP+LSCH-Lw,3*Bst/2+LSP+LSCH    !Zeichnen der Spulen

rectng,(Radau*2-LST)/2+Lst/2+LFE/2-Bsch-BW-LSP/2,(Radau*2-LST)/2+Lst/2+LFE/2-Bsch-LSP/2,3*Bst/2+LSP+LSCH-Lw,3*Bst/2+LSP+LSCH    !Zeichnen der Spulen
rectng,(Radau*2-LST)/2+Lst/2+LFE/2+LSP/2,(Radau*2-LST)/2+Lst/2+LFE/2+BW+LSP/2,3*Bst/2+LSP+LSCH-Lw,3*Bst/2+LSP+LSCH              !Zeichnen der Spulen

/Pnum,area,1

/prep7		  !Aufruf Preprozessor

aovlap,all	  !alle Flaechen ueberlappen

aadd,5,6,15       !Verbinde Ferritflaechen

!Zuordnung von Material und ElementTyp zu den einzelnen Flaechen und Vernetzung

!Plane53,AZ      : Elem.Typ1 Luft, Ferrit Eisen
!Plane53,AZ Curr : Elem.Typ2 Kupfer
!Infin110,AZ     : Elem.Typ3 Boundary

!A8,A9,A10,A11   :Leiter
!A1		 :Ferritkern
!A13,A14	 :Stahl
!A16		 :Luft innen A12:luft aussen					

!------------------------Vernetzung Luft aussen--------------------

asel,s,area,,12    			!Flaeche 12 sei Mat.1,Type4
aatt,1,,3

csys,1
lsel,s,loc,x,(radau+radmi)/2		!Festlegen des Linienschwerpunktes 
lesize,all,,,1,1/3			!nichtlineare Teilung
mshape,0,2d				!nur Rechtecke
mshkey,1
k=Radau-radmi				!mapped mashing
esize,k!,11				!Radiale Unterteilung der Flaeche
type,3  $mat,1
amesh,12				!Vermaschen der Flaeche 4
allsel

!-------------------------Vernetzung Ferrit------------------------

asel,s,area,,1    			 !Flaeche 1 sei Mat.4,Type1
aatt,4,,1

esize,Lst/120 ! 240
mshkey,2
type,1  $mat,2
amesh,1		   			 !Vermaschen der Flaeche 1

!------Vernetzung Kupfer--------
asel,s,area,,8,11  			 !Flaeche 8-11 sei Mat.2,Type2
aatt,2,,2

esize,Lst/120 ! 240
mshkey,1
type,2  $mat,2
amesh,8,11				 !Vermaschen der Flaeche 8-11


!------Vernetzung Stahl1---------

asel,s,area,,13,14 			 !Flaeche 13-14 sei Mat.3,Type1
aatt,3,,1

esize,Lst/120 ! 240
mshkey,0
type,1  $mat,2
amesh,13				!Vermaschen der Flaeche 13

!------Vernetzung Stahl2---------

esize,Lst/80
mshkey,0
type,1  $mat,2
amesh,14				!Vermaschen der Flaeche 14

!------Vernetzung Luft innen----

asel,s,area,,16    			!Flaeche 16 sei Mat.1,Type1
aatt,1,,1
mshape,1,2d
esize
mat,1
amesh,16	  			!Vermaschen der Flaeche 16

allsel
!-----------Festlegung der Spulenkonstanten-----------------------------------

/PREP7  

!------Festlegung der Erregungsart für jedes Spulensegment (Flächen A9-A13)----

!cpdele,all

!-------1--------

aplot   
asel,s, , , 8 
esla,r  
eplot 
 
cm,WICK8,elem
nsle,s 
nplot  
CP,next,CURR,all
CP,next,EMF,all
*get,mast8,node,,num,min

allsel
!-------2--------

aplot   
asel,s, , , 9 
esla,r  
eplot
cm,WICK9,elem
nsle,s 
nplot  
CP,next,CURR,all
CP,next,EMF,all
*get,mast9,node,,num,min    ! ##### fuer Kopplung mit CIRCU124

allsel

!-------3--------

aplot   
asel,s, , , 10 
esla,r  
eplot  
cm,WICK10,elem
nsle,s 
nplot  
CP,next,CURR,all
CP,next,EMF,all
*get,mast10,node,,num,min

allsel

!-------4--------

aplot   
asel,s, , , 11 
esla,r  
eplot  
cm,WICK11,elem
nsle,s 
nplot  
CP,next,CURR,all
CP,next,EMF,all
*get,mast11,node,,num,min

allsel

! ############ Reals fuer Wicklungen
fuellfac = 1.0        ! Spulen-Fuellfaktor
asel,s,area,,8
asum
*get,coilarea,area,8,area
r,8,coilarea,N,1, 1,fuellfac   ! DIRZ = 1/-1 gibt die Stromrichtung an
r,9,coilarea,N,1,-1,fuellfac
r,10,coilarea,N,1,-1,fuellfac
r,11,coilarea,N,1, 1,fuellfac
cmsel,s,WICK8
emodif,all,real,8
cmsel,s,WICK9
emodif,all,real,9
cmsel,s,WICK10
emodif,all,real,10
cmsel,s,WICK11
emodif,all,real,11
allsel

!CIRCU124 fuer Speisung und Induktion
freq = 1000            ! Hz
r,4,0,Strom,freq       ! Sinuserregung
r,5,1

csys,0
*get,nmax,node,,num,max
mm = 0.001
n,nmax+1,20*mm,-10*mm
n,nmax+2,20*mm,-20*mm
n,nmax+3,25*mm,-10*mm
n,nmax+4,25*mm,-20*mm
n,nmax+5,40*mm,-10*mm
n,nmax+6,40*mm,-20*mm
n,nmax+7,45*mm,-10*mm
n,nmax+8,45*mm,-20*mm
type,4 				!124er Stromquelle
real,4
e,nmax+1,nmax+2
e,nmax+3,nmax+4
e,nmax+5,nmax+6
e,nmax+7,nmax+8
type,5
real,5
e,nmax+1,nmax+2,mast8		!gekoppelte Knoten
e,nmax+3,nmax+4,mast9
e,nmax+5,nmax+6,mast10
e,nmax+7,nmax+8,mast11

!--------Randknoten setzen------

CSYS,1 
lsel,s,radius,,radau
nsll,s,1
esln
eplot
emodif,all,type,3
sf,all,inf
allsel

!------------ Symmetrieachse festlegen---------------
csys,0
nsel,s,loc,y,0   !Alle Knoten der x-Achse erfassen
D,all,az,0
allsel
!-----------Eingabe der Erregung----------------------
/cont,,7,0.001,1.2

/solu
! ########## Parameter freq s.o.
antype,trans
trnopt,full

KBC,1
OUTRES,ALL,ALL,
cnvtol,csg
cnvtol,amps
cnvtol,vltg

! Dauermagnetfeld statisch aufbauen
timint,off
TIME,1e-6
nsubst,1
solve
!-----------------------------------------------------
! transient
perioden = 1
timint,on
time,perioden/freq
autots,1
DELTIM,1/freq/100,1/freq/100,1/freq/10
solve

/post1 !post26
SET,FIRST
/EFACET,,,,,1   
PLNSOL, B,SUM, 0
!set,list
!plvect,jt,z
!set,,,,,0.25e-3
/cont,1,10,0,,0.25
/auto,1
/ZOOM,1,SCRN,0.3333,0.01,0.8,-0.425 
/ZOOM,1,SCRN,0.3333,0.01,0.8,-0.425 
/ZOOM,1,SCRN,0.3333,0.01,0.8,-0.425 
  
/FOC,1,,-0.3,,1 
/REP,FAST   
/FOC,1,,-0.3,,1 
/REP,FAST   
/FOC,1,,-0.3,,1 
/REP,FAST   
/FOC,1,,-0.3,,1 
/REP,FAST   
/FOC,1,,-0.3,,1 
/REP,FAST   
/FOC,1,,-0.3,,1 
/REP,FAST   
/ 
/DIST,1,1.37174211248,1 
/REP,FAST   
/DIST,1,1.37174211248,1 
/REP,FAST 
/replot

/COLOR,SMIN,WHITE
/COLOR,CNTR,BLUE,1  
/COLOR,CNTR,CBLU,2  
/COLOR,CNTR,CYAN,3  
/COLOR,CNTR,GCYA,4  
/COLOR,CNTR,GREE,5  
/COLOR,CNTR,YGRE,6  
/COLOR,CNTR,YELL,7  
/COLOR,CNTR,ORAN,8  
/COLOR,CNTR,RED,9   

/rgb, index, 100, 30, 0, 13,
/Tlabel,1.5,-0.91,Tesla
!PLVECT,B, , , ,VECT,ELEM,ON,0   
PLNS,B,SUM  
ANTIME,70,0.1, ,0,2,0,0.1E-02
/title,Modell: Wirbelstromsensor mit magnetischer Arbeitspunktstabilisierung 
/replot