!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Dies ist die gesamte Berechnung selbst geschrieben!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!Programm aus dem Forum von Petrov
!mit provisorischen Änderungen von Schätzing
!ohne Prüfung der Richtigkeit und Sinnfälligkeit
!Vernetzung nicht geeignet
!
!12.12.2011

finish
/clear
*abbr,petrov,/input,petrov,txt



/title, Erwaermung der Kupferplatte
/VIEW,all,1,1,1 !Schräge Ansicht
/VUP,all,Z !Für die Ansicht wird das globale Koordinatensystem definiert
/PREP7
/UNITS,mks !SI System (mit Celsius statt Kelvin)



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!GEOMETRIE BAUEN!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!definieren der Größen
nutbreite = 0.01
zahnbreite = 0.0066
halberzahn = zahnbreite/2
nutenanzahl = 24
nutteilung = halberzahn+nutbreite+halberzahn !Ein Element entspricht der Polteilung
laengeprim = nutenanzahl*nutteilung !Setzt sich aus 24*nutteilung zusammen
nuttiefe = 0.015
hoeheprim = 0.03          !Höhe des gesamten Primärteils
unterbau = hoeheprim-nuttiefe   !Also die Platte, auf die die Zähne kommen
luftspalt = 0.005
hoehesek = 0.005
breiteprim = 0.25
breitesek = 0.3
ueberhang = (breitesek-breiteprim)/2 !Überhang, den das Sekundärteil über steht (jeweils links und rechts)
ueberhang2 = 0.02 !Anzahl an cm, die die Luftbox den Primärteil umspannen soll (Also Abstand...)
laengestab = (2*ueberhang)+(breiteprim)+(2*ueberhang2)  !Länge eines Leiterstabes, wieviel er pro Seite rausschaut drin verarbeitet

!Erzeugen der mittlere Leiterstäbe
!Auf diese Weise teilen sich Nutstäbe und Primäreisen immer drei Flächen nach vsbv (s.u.)
*do,i,0,0.39,nutteilung    !Inkrementierung um den Wert Nutteilung, bis 0,39 überschritten wird
    block,i+halberzahn,i+halberzahn+nutbreite,0,breiteprim,unterbau,unterbau+nuttiefe !Erzeugung der Nut (die schaut nicht aus dem Primärteil raus)
    !block,i+halberzahn,i+halberzahn+nutbreite,-ueberhang-ueberhang2,breiteprim+ueberhang+ueberhang2,unterbau,unterbau+nuttiefe !Erzeugung der Nut (die schaut aus dem Primärteil raus)
*enddo

!Erzeugung des Metallblocks des Primärteils
block,0,laengeprim,0,breiteprim,0,hoeheprim

!Boolsche Funktion, um die Nutenstäbe vom Primärteil abzuziehen
!Auf diese Weise wird sichergestellt, dass die Nutenstäbe und der Eisenteil sich gemeinsame Flächen teilen
!Volumennummern: Stäbe #1-24, Eisenteil #25
vsel,s,volu,,1,25 !Alle Volumen auswählen
vsbv,25,all,,delete,keep  !alle (bis auf #25) von #25 abziehen und neu anlegen (#26). (blank) = teilen sich die Schnittflächen. #25 wegwerfen (delete) und #1-#24 behalten (keep)

!Volumennummern: Stäbe #1-24, Eisenteil #26
numcmp,volu !Nummerierung wieder "gerade rücken", damit keine leere Nummer (25) da ist
!Volumennummern: Stäbe #1-24, Eisenteil #25

!Erzeugung der Kupferplatte #26
block,0,laengeprim,-ueberhang,breiteprim+ueberhang,hoeheprim+luftspalt,hoeheprim+luftspalt+hoehesek

!Erzeugung der Leiterstab-Überhänge
*do,i,0,0.39,nutteilung    !Inkrementierung von 0 um den Wert Nutteilung, bis 0,39 überschritten wird
    block,i+halberzahn,i+halberzahn+nutbreite,0,-ueberhang-ueberhang2,unterbau,unterbau+nuttiefe
    block,i+halberzahn,i+halberzahn+nutbreite,breiteprim,breiteprim+ueberhang+ueberhang2,unterbau,unterbau+nuttiefe
*enddo
!Volumennummern: #27-74 (liegt an Nut Nummer= (Volumennummer-26)/2) (aufrunden!)

!Immer die drei Nutteile auswählen, die zusammen gehören und zusammenfügen
*do,i,1,24,1
    vsel,s,volu,,i
    vsel,a,volu,,25+(2*i),26+(2*i)
    vadd,all
*enddo

!Neue Nummerierung verpassen
numcmp,volu
!Neue Nummerierung: #1-23 & #26 Nutenstäbe (komplett mit Überhang)
!                   #24 Primärteil-Eisen
!                   #25 Kupferplatte


!/eof

!Erstellung der Umgebungsluft Box #27 (inklusive Luftspalt)
block,-ueberhang2,laengeprim+ueberhang2,-ueberhang-ueberhang2,breiteprim+ueberhang+ueberhang2,-ueberhang2,ueberhang2+hoeheprim+luftspalt+hoehesek

!Boolsche Oberation, um die festen Materialien, aus der Luftgeometrie "rauszuschneiden"
vsel,s,volu,,all
vsbv,27,all,,delete,keep !löschen der alten Luftbox, behalten der Metalle (nuten,eisen, kupferplatte)
!Neue Nummerierung für Luftbox ist #28

!Leere Volumennummer 27 stopfen
numcmp,volu

!Neue Nummerierung: #1-23 & #26 Nutenstäbe (komplett mit Überhang)
!                   #24 Primärteil-Eisen
!                   #25 Kupferplatte
!                   #27 Luftbox

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Materialeigenschaften definieren!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!Material einfach eine Nummer zugewiesen
eisen = 1
kupfer = 2
luft = 3

!/eof


et,1,236,1         !Nummer 1, Eigenschaft von 236, Freiheitsgrade(AZ&Volt) --> Wirbelströme zugelassen (mit type,1 aufzurufen)                        <-- für das Kupfer
et,2,236,1,,,,1    !Nummer 2, Eigenschaften von 236, Freiheitsgrade (AZ&Volt) ---> Keyopt(5) = 1 -> Wirbelströme unterdrückt  (mit type,2 aufzurufen) <-- für das geblechte Eisen
et,2,236,1
et,3,236,0         !Nummer 3, Eigenschaften von 236, Freiheitsgrade (AZ) ---> Dadurch AUCH keine Wirbelströme   (mit type,3 aufzurufen)               <-- für die Luft


!Luft
mp,murx,luft,1 !setze Permeabilität von Luft auf 1

!Eisen
!tb,bh,eisen,,29  !Data table für das Eisen mit 29 Einträgen
!tbpt,,400,0.464   !es handelt sich um eine BH Kurve (oben angegeben)
!tbpt,,800,0.865
!tbpt,,1200,1.145
!tbpt,,1600,1.275
!tbpt,,2000,1.355
!tbpt,,2400,1.412
!tbpt,,2800,1.458
!tbpt,,3200,1.495
!tbpt,,3600,1.522
!tbpt,,4000,1.55
!tbpt,,5000,1.604
!tbpt,,6000,1.645
!tbpt,,7000,1.678
!tbpt,,8000,1.702
!tbpt,,10000,1.75
!tbpt,,14000,1.81
!tbpt,,16000,1.863
!tbpt,,22000,1.905
!tbpt,,26000,1.934
!tbpt,,30000,1.96
!tbpt,,40000,2.015
!tbpt,,50000,2.052
!tbpt,,60000,2.08
!tbpt,,80000,2.128
!tbpt,,100000,2.165
!tbpt,,120000,2.2
!tbpt,,140000,2.23
!tbpt,,200000,2.308
!tbpt,,5000000,8.34

mp,rsvx,eisen,1/(10e6)    !electrical resistivity fürs eisen festgelegt

mp,murx,eisen,2000

!Kupfer
mp,rsvx,kupfer,1/(59e6)  !elektrischer Widerstand für das Kupfer festgelegt
mp,murx,kupfer,1         !my_r von kupfer = 1


!thermische Eigenschaften
mp,c,kupfer,385   !Wärmekapazität
mp,c,eisen,449
mp,c,luft,1.005
mp,kxx,luft,0.0259  !therm leitfähigkeit
mp,kxx,kupfer,400 
mp,kxx,eisen,80


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Erregerströme definieren!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!Phase +R
i1re = 70
i1im = 0

!Phase -S
i2re = 35
i2im = -35*sqrt(3)

!Phase +T
i3re = -35
i3im = -35*sqrt(3)

!Phase -R
i4re = -70
i4im = 0

!Phase S
i5re = -35
i5im = 35*sqrt(3)

!Phase -T
i6re = 35
i6im = 35*sqrt(3)

!/eof

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Vernetzung der Volumen und Eigenschaften von oben hinzufügen!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!Jedesmal das Volumen ausgewählt und Parameter drauf angewendet, dann neu ausgewählt und neu angewendet bis das eine fertig ist

!Neue Nummerierung: #1-23 & #26 Nutenstäbe (komplett mit Überhang)
!                   #24 Primärteil-Eisen
!                   #25 Kupferplatte
!                   #27 Luftbox

!!!!!!Kupferstäbe: Eigenschaften hinzufügen und Volumen vermeshen
vsel,s,volu,,1,23    !23 Stäbe
vsel,a,volu,,26    !24ter Stab (erste Stelle)
     alls,below,volu
     aesize,all,0.05
     lesize,all,0.05

vsel,s,volu,,1,23
vsel,a,volu,,26
     type,2        !sets the attribut pointer (keine Wirbelströme)
     mat,2         !sets the material pointer (Material: Kupfer)

vsel,s,volu,,1,23
vsel,a,volu,,26
!     mshape,1,3D   !Standardmäßig Rechteckiges Netz, hier auf Tetraeder umgestellt (0= Rechteck)
!     vmesh,all
*do,i,1,23
vsweep,i
*enddo
vsweep,26

!/eof

!!!!!!Primärteileisen: Eigenschaften hinzufügen und Volumen vermeshen
!vsel,s,volu,,24
!     alls,below,volu
!     aesize,all,0.005
!     lesize,all,0.005

vsel,s,volu,,24
     type,2    !(Keine Wirbelströme)
     mat,1     !Material Eisen
type,3  !Luft
mat,3

!vsel,s,volu,,24
!     mshape,1,3D !Auf tetraeder Mesh umstellen. Vermeshed alles, verbraucht aber mehr Rechenleistung.
!     vmesh,all

vsweep,all




!!!!!!Kupferplatte: Eigenschaften hinzufügen und Volumen vermeshen
!vsel,s,volu,,25
!     alls,below,volu
!     aesize,all,0.0015
!     lesize,all,0.0015
     
vsel,s,volu,,25
     type,1     !Wirbelströme zugelassen
     mat,2      !Material Kupfer
     
!vsel,s,volu,,25
     mshape,0,3D  !Auf Quaderförmiges Mesh umstellen
esize,hoehesek  !Elementgröße=Plattendicke
     vmesh,all

!/eof
     
!!!!!Luft: Eigenschaften hinzufügen und Volumen vermeshen
vsel,s,volu,,27
     alls,below,volu
     aesize,all,0.04
     lesize,all,0.04

vsel,s,volu,,27
     type,3
     mat,3
     
vsel,s,volu,,27
     mshape,1,3D
     vmesh,all

!/eof


!Nummerierung der Aufprägungsflächen der Nutenstäbe:
!        -Z      +Z
!Stab 1: #183 & #190
!Stab 2: #195 & #202
!Stab 3: #207 & #214
!Stab 4: #219 & #226
!Stab 5: #231 & #238
!Stab 6: #243 & #250
!Stab 7: #255 & #262
!Stab 8: #267 & #274
!Stab 9: #279 & #286
!Stab 10: #291 & #298
!Stab 11: #303 & #310
!Stab 12: #315 & #322
!Stab 13: #327 & #334
!Stab 14: #339 & #346
!Stab 15: #351 & #358
!Stab 16: #363 & #370
!Stab 17: #375 & #382
!Stab 18: #387 & #394
!Stab 19: #391 & #406
!Stab 20: #403 & #418
!Stab 21: #415 & #430
!Stab 22: #427 & #442
!Stab 23: #439 & #454
!Stab 24: #451 & #466
!Umständlich rausgesucht

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Stromaufprägung!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!Zuerst jeweils 4 Nuten zu einem Set bündeln (d.h. deren Knoten verkoppeln)


!Set 1
asel,s,area,,190    !Kupferstab Nummer 1 Kopffläche (Volumennummer 26)
asel,a,area,,202    !Kupferstab Nummer 2 Kopffläche (Volumennummer 1)
asel,a,area,,214    !Kupferstab Nummer 3 Kopffläche (Volumennummer 2)
asel,a,area,,226   !Kupferstab Nummer 4 Kopffläche (Volumennummer 3)
nsla,s,1 !Node SeLect Area = Alle Knoten (nodes) auswählen, die mit der vorher ausgewählten Fläche zusammenhängen (interior to areas, lines and at keypoints)
cp,1,volt,all  !die erste 1 ist die Nummer des Sets. Die dafür ausgewählten Knoten (nodes) bilden ein set. Also hier die Knoten der ersten vier Leiterstäbe, usw
f,ndnext(0),amps,i1re,i1im   !f,alle Knoten davon,Stromfluss einstellen,Realteil,Imaginärteil


!Set 2
asel,s,area,,238
asel,a,area,,250
asel,a,area,,262
asel,a,area,,274
nsla,s,1
cp,2,volt,all
f,ndnext(0),amps,i2re,i2im

!Set 3
asel,s,area,,286
asel,a,area,,298
asel,a,area,,310
asel,a,area,,322
nsla,s,1
cp,3,volt,all
f,ndnext(0),amps,i3re,i3im

!Set 4
asel,s,area,,334
asel,a,area,,346
asel,a,area,,358
asel,a,area,,370
nsla,s,1
cp,4,volt,all
f,ndnext(0),amps,i4re,i4im

!Set 5
asel,s,area,,382
asel,a,area,,394
asel,a,area,,406
asel,a,area,,418
nsla,s,1
cp,5,volt,all
f,ndnext(0),amps,i5re,i5im

!Set 6
asel,s,area,,430
asel,a,area,,442
asel,a,area,,454
asel,a,area,,466
nsla,s,1
cp,6,volt,all
f,ndnext(0),amps,i6re,i6im

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Spannung für die gegenüberliegenden Flächen gleich Null setzen!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!Set 1
asel,s,area,,183    !Kupferstab Nummer 1 andere Kopffläche (Volumennummer 26)
asel,a,area,,195    !Kupferstab Nummer 2 andere Kopffläche (Volumennummer 1)
asel,a,area,,207    !Kupferstab Nummer 3 andere Kopffläche (Volumennummer 2)
asel,a,area,,219    !Kupferstab Nummer 4 andere Kopffläche (Volumennummer 3)
nsla,s,1 !Node SeLect Area = Alle Knoten (nodes) auswählen, die mit der vorher ausgewählten Fläche zusammenhängen (interior to areas, lines and at keypoints)
d,all,volt,0     !An allen ausgewählten Knoten die Spannung =0 setzen


!Set 2
asel,s,area,,231
asel,a,area,,243
asel,a,area,,255
asel,a,area,,267
nsla,s,1
d,all,volt,0

!Set 3
asel,s,area,,279
asel,a,area,,291
asel,a,area,,303
asel,a,area,,315
nsla,s,1
d,all,volt,0

!Set 4
asel,s,area,,327
asel,a,area,,339
asel,a,area,,351
asel,a,area,,363
nsla,s,1
d,all,volt,0

!Set 5
asel,s,area,,375
asel,a,area,,387
asel,a,area,,399
asel,a,area,,411
nsla,s,1
d,all,volt,0

!Set 6
asel,s,area,,423
asel,a,area,,435
asel,a,area,,447
asel,a,area,,459
nsla,s,1
d,all,volt,0

allsel
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Vektorpot an den Außenflächen parallel legen (Randbedingung)!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!Die Außenflächen der Luftbox gewählt
asel,s,area,,349,350
asel,a,area,,352,353
asel,a,area,,141,142

!Alle ihre Knoten angewählt
nsla,s,1
nsla,

!Vektorpotential =0 setzen
d,all,az,0
                                 
finish!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!SOLUTION!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/solu!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!Auf harmonische Berechnung umstellen
antype,harm

!Frequenz einstellen
harfrq,56e-01
  
CNVTOL,A, ,0.01, ,    !Sets convergence values for nonlinear analyses. A = Vektorpotential
allsel
!Lösen                                                    
solve
                                                         
finish!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Postprozessor!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/post1!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!Kupferplatte auswählen
vsel,s,volu,,25
eslv
ples,jheat

/eof

esel,s,type,,2

set,1,1
ples,jt,y

set,1,1,,1
ples,jt,y