! -- Parameter -------------------------------
strom			= 10 		! [A]
windungen		= 100 		! []

ajhoehe			= 20e-2		! [m]
ajbreite		= 40e-2		! [m]
ajtiefe			= 4e-2		! [m]
dicke			= ajtiefe	! [m]
ijbreite		= dicke		! [m]
lspalthoehe		= 3e-2		! [m]
ijhoehe			= (ajhoehe-2*dicke)/2

charSpuleDY		= ijhoehe/10
charSpuleDZ		= ijhoehe/2


/PREP7
ET,1,SOLID96			! 3-D Magnetic Scalar Solid
ET,2,SOURC36			! Element für Spule
ET,3,PLANE13			! 2-D Coupled-Field Solid


MP,MURX,1,1 			! Relative Permeabilität für Material 1 (Luft); Magnetic relative permeabilities (also MURY, MURZ).
MP,MURX,2,2000 			! Relative Permeabilität für Material 2 (Eisen)
R,2,1,strom*windungen,charSpuleDY,charSpuleDZ ! Spuleneigenschaften (Preprocessor > real constants > Add/Edit/Delete > Add, Type 2 Sourc36
						! 1 TYPE Source type 
						! 2 CUR Total current through source 
						! 3 DY Characteristic Y dimension 
						! 4 DZ Characteristic Z dimension 


/PNUM,NODES,0
/PNUM,KP,0
/PNUM,AREAS,1
/PNUM,VOLUM,1
/PNUM,LINES,0

! -- Geometrie -------------------------------
! Aufbau 2D, dann extrudieren und zweimal spiegeln
K,1,0,lspalthoehe/2,0
K,2,-ijbreite/2,lspalthoehe/2,0
K,3,-ijbreite/2,ajhoehe/2-dicke,0
K,4,-ajbreite/2+dicke,ajhoehe/2-dicke,0
K,5,-ajbreite/2+dicke,0,0
K,6,-ajbreite/2,0,0
K,7,-ajbreite/2,ajhoehe/2,0
K,8,0,ajhoehe/2,0

! Verbinden der KPs zu Linien
*DO,numkp,1,7
	L,numkp,numkp+1
*ENDDO
L,numkp+1,1

! Fläche erstellen
AL,ALL


! Luftspalt mit Material füllen
K,10,0,0,0
K,11,-ijbreite/2,0,0

L,1,10
L,10,11
L,11,2
L,2,1

ksel,none
ksel,a,kp,,1
ksel,a,kp,,2
ksel,a,kp,,10
ksel,a,kp,,11

LSLK,S,1 ! Alle Linien an den KPs auswählen

AL,ALL

ALLSEL
APLO


! Meshen Luftspalt
ESIZE,0.3

! Linie Luftspalt unten
LSEL,S,LOC,Y,0
LSEL,R,LOC,X,-ijbreite/2,0
*GET,numl,LINE,,NUM,MAX
lesi,numl,,,10
! Linie Luftspalt oben
LSEL,S,LOC,Y,lspalthoehe/2
LSEL,R,LOC,X,-ijbreite/2,0
*GET,numl,LINE,,NUM,MAX
lesi,numl,,,10


! Linie Luftspalt rechts
LSEL,S,LOC,X,0
LSEL,R,LOC,Y,0,lspalthoehe/2
*GET,numl,LINE,,NUM,MAX
lesi,numl,,,10
! Linie Luftspalt links
LSEL,S,LOC,X,-ijbreite/2
LSEL,R,LOC,Y,0,lspalthoehe/2
*GET,numl,LINE,,NUM,MAX
lesi,numl,,,10



TYPE,3
AMESH,2  ! Nummer noch anpassen!

ALLSEL

lsel,s,line,,2
lsel,a,line,,3
lsel,a,line,,4
!LCCAT,all
lsel,s,line,,6
lsel,a,line,,7
lsel,a,line,,8
!LCCAT,all

ALLSEL
LESIZE,1, , ,10, , , , ,1 
LESIZE,5, , ,10, , , , ,1 
LESIZE,2, , ,15
LESIZE,3, , ,15
LESIZE,4, , ,15
LESIZE,6, , ,15
LESIZE,7, , ,15
LESIZE,8, , ,15
AMESH,1

! 3D Element aktivieren, Anker
TYPE,1

! Erstellte Geometrie mit Mesh extrudieren, Anker
EXTOPT,ESIZE,10,0,   
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0   
MAT,2   
!REAL,1  
!ESYS,0  
VEXT,1, , ,0,0,ajtiefe,,,,   


! 3D Element aktivieren, Luft
TYPE,1

! Erstellte Geometrie mit Mesh extrudieren, Luft
EXTOPT,ESIZE,10,0,   
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0   
MAT,1   
!REAL,1  
!ESYS,0  
VEXT,2, , ,0,0,ajtiefe,,,,  

! -- Spulen -----------------------------
! Erstellen der Spule oben
! KS definieren
LOCAL,11,,0,ajhoehe/2-dicke-ijhoehe/2+charSpuleDY,ajtiefe/2

N,100000,-ijbreite
N,100001,,,-ijbreite
N,100002,0,0,0
TYPE,2
MAT,3
REAL,2
E,100000,100001,100002

! Erstellen der Spule unten
! KS definieren
LOCAL,12,,0,-(ajhoehe/2-dicke-ijhoehe/2+charSpuleDY),ajtiefe/2

N,100003,-ijbreite
N,100004,,,-ijbreite
N,100005,0,0,0
TYPE,2
MAT,3
REAL,2
E,100003,100004,100005


/ESHA,1
EPLO

! Reset CS
CSYS,0

! Plane ET löschen
ETDELE,3

ALLSEL

! Volumen spiegeln
VSYMM,Y,1, , , ,0,0 
VSYMM,Y,2, , , ,0,0 

VSYMM,X,ALL, , , ,0,0 


/SOLU
ALLSSEL

MAGSOLV,3

/POST1
SET,LAST
/RESET

PLVECT,B, , , ,VECT,ELEM,ON,0   
/VSCALE,1,1,0