fini
	/clear
!Vorgabe von Ax, Ay, Az
!für ein gewünschtes Bz
!hier soll Bz vorgebbar sein
!------Parameter--------
H=1												!Magnetische Feldstäke in Tesla
mur=1											!Perameabilität des Eisens (gemittelt)
pi=3.141582654
mu0=pi*4e-7
abx=1 												!Abmessung in x-Richtung
aby=1 												!Abmessung in y-Richtung
abz=1 												!Abmessung in z-Richtung
Bx=2 												!Flussdichte in x-Richtung
By=2 												!Flussdichte in y-Richtung
Bz=2		 										!Flussdichte in z-Richtung
!------Modellerstellung-----------------
/prep7
	et,1,solid97
	et,2,solid97
	et,3,solid97
	mp,murx,2,1									!Permeabilität - Luft
	mp,murx,3,50000								!Permeabilität - Eisen
	!----Geometrie-----
	/view,,1,2,3
	block,-abx,abx,-abx,abx,-abx,abx,				!Luft
	block,-abx/2,abx/2,-abx/2,abx/2,-abx/2,abx/2,	!Probekörper
	vovlapa,all
	
	
	!----Vernetzung----------------
	type,3
	mat,3
	esize,abx/8
	mshkey,1
	mshape,0,3d
	vmesh,2
	MSHKEY,0		! Use mapped meshing
	MSHAPE,1 		! Use hexahedra 
!	LESIZE,ALL 	
!	LPLOT
	mat,2
	type,2
	esize,abx/6
	vmesh,3
	!----Randbedingungen----
	nsel,s,loc,x,-abx
	nsel,a,loc,x,abx
	nsel,a,loc,y,-abx
	nsel,a,loc,y,abx
	nsel,a,loc,z,-abx
	nsel,a,loc,z,abx
	nplot
	*get,nanz,node,,count
	*get,nmin,node,,num,min
	narb=nmin
	*do,i,1,nanz
		d,narb,ax,1*(-bz/2)*nz(narb)
		d,narb,az,1*(+bz/2)*nx(narb)
		d,narb,ay,0
		narb=ndnext(narb)
	*enddo
	esel,s,mat,,3 				!Elemente des Blocks
	cm,klotz,elem 				!Komponentendefinition
	fmagbc,'klotz' 				!für Kraftauswertung
	allsel
!------SOLVER----------------------------------------------------------	
/solu
	solve
!----------------------
esel,s,mat,,2
/post1
	!allsel
	!plvect,h,,,,vect,,on
	!----Flussdichte--------
	!plns,b,sum
	!-----Kraftermittlung------
	!esel,all
	!etable,kraft_y1,fmag,y 		!Maxwell
	!etable,kraft_y2,nmisc,4	 	!virtual work
	!ssum
	!*get,kraft1,ssum,,item,kraft_y1
	!*get,kraft2,ssum,,item,kraft_y2
	!/wait,4
	!plvect,fmag,,,,vect,,on
	
	esel,s,mat,,3
	nsel,s,ext
	esln,s,,1
	!nsle,s,,1
	*get,nanz,node,,count
	*get,nmin,node,,num,min
	narb=nmin
	*dim,nkraft,table,nanz,4
	*do,i,1,nanz
		*get,krx,node,narb,fmag,x
		*get,kry,node,narb,fmag,y
		*get,krz,node,narb,fmag,z
		nkraft(i,1)=narb
		nkraft(i,2)=krx
		nkraft(i,3)=kry
		nkraft(i,4)=krz
		narb=ndnext(narb)
	*enddo
	
!*ABBR,Einlesen,/input,D:\UNI\Diplomarbeit\ANSYS-Modell\probe.txt