!---Macro zur Auswertung der Druck- und Tragbilder


!---Auswertung für Lager 1

csys, 21
rsys, 21
dsys, 0

!---axiale Auswertung ja Nadel
anz = es+1
*dim, rforce1, array, anz+2, 4, an	!array definieren in das später die Kraftwerte geschrieben werden

*do, k, 1, an
	lsel, s, , , k			!entsprechende Linien ansprechen auf denen die Federelemente sitzen
	nsll, s, 1
	esln, s, 0
	esel, r, type, , 11

nnum = 0				!Nummern auf 0 setzen um alle im set zu erfassen
*do, j, 2, anz+1
	*get, nnum, elem, nnum, nxth		!Elementnummer bestimmen und in die Tabelle schreiben
	*vfill, rforce1(j,1,k), data, nnum
*enddo

*do, i, 2, anz+1 			
	*get, koord, elem, rforce1(i,1,k), cent, z	!z-Koordinate des Federelements bestimmen
	*vfill, rforce1(i,2,k), data, koord

	*get, frk, elem, rforce1(i,1,k), smisc, 1 	!Elementkraft bestimmen
	*if, frk, lt, 0, then				!negative Kräft positiv machen und positive zu Null
		*vfill, rforce1(i,3,k), data, -frk
	*else
		*vfill, rforce1(i,3,k), data, 0
	*endif
	
	*get, dehn1, elem, rforce1(i,1,k), nmisc, 1 
	*vfill, rforce1(i,4,k), data, dehn1

*enddo

rforce1(1,1,k) = 0	!restliche Tabelle ausfüllen um grafische Darstellung zu verbessern
rforce1(1,2,k) = 0
rforce1(1,3,k) = 0
rforce1(1,4,k) = 0
rforce1(anz+2,1,k) = 0
rforce1(anz+2,2,k) = l
rforce1(anz+2,3,k) = 0
rforce1(anz+2,4,k) = 0

*enddo

!---radiale Auswertung nach Koordinate

*dim, trag1, array, an, 14	!array anlegen um die Werte zu speichern

*do, i, 1, an
	wink = (360/an)*(i-1)		!Winkelkoordinate des Wälzkörpers berechnen und eintragen
	*vfill, trag1(i,1), data, wink
*enddo


*do, h, 1, 13
	*do, k, 1, an
		*vfill, trag1(k,h+1), data, rforce1(h,3,k)	!Kraftwert je Wälzkörper auslesen und in Array speichern
	*enddo
*enddo
allsel

!---axiale Auswertung ja Nadel in alle Richtungen

anz = es+1
*dim, rforce11, array, anz+2, 6, an

*do, k, 1, an
	lsel, s, , , k		!Nodes selektieren die zu dem jeweiligen Federelement auf der Linie gehören
	nsll, s, 1
	esln, s, 0
	esel, r, type, , 11
	nsle, s, 1
	nsel, r, loc, x, dz/2

nnum = 0
*do, j, 2, anz+1
	*get, nnum, node, nnum, nxth
	*vfill, rforce11(j,1,k), data, nnum
*enddo

*do, i, 2, anz+1 			!Reaktionskräfte des jeweiligen Nodes bestimmen und abspeichern
	nsel, s, node, , rforce11(i,1,k)
	fsum, , 
	*get, frx, fsum, 0, item, fx
	*get, fry, fsum, 0, item, fy
	*get, frz, fsum, 0, item, fz 
	*get, koord, node, rforce11(i,1,k), loc, z
	!dru = druck(1,1)+(druck(2,1)*koord)
	
	*vfill, rforce11(i,2,k), data, koord
	*vfill, rforce11(i,3,k), data, frx
	*vfill, rforce11(i,4,k), data, -fry
	*vfill, rforce11(i,5,k), data, frz
	*vfill, rforce11(i,6,k), data, dru
*enddo

rforce11(1,1,k) = 0
rforce11(1,2,k) = 0
rforce11(1,3,k) = 0
rforce11(1,4,k) = 0
rforce11(anz+2,1,k) = 0
rforce11(anz+2,2,k) = l
rforce11(anz+2,3,k) = 0
rforce11(anz+2,4,k) = 0

*enddo
	
*use, schreib21.mac, nr	!Auswertung von Nadel(nr) in Datei schreiben
*use, schreib31_neu.mac,	
*use, schreib41.mac, nr

!------------------------------------------------------------------------
!---Auswerung für Lager 2

csys, 20
rsys, 20
dsys, 0

!---axiale Auswertung ja Nadel
anz = es+1
*dim, rforce2, array, anz+2, 4, an	!array definieren in das später die Kraftwerte geschrieben werden

*do, k, 1, an
	lsel, s, , , k+1000		!entsprechende Linien ansprechen auf denen die Federelemente sitzen
	nsll, s, 1
	esln, s, 0
	esel, r, type, , 11
	
nnum = 0				!Nummern auf 0 setzen um alle im set zu erfassen
*do, j, 2, anz+1
	*get, nnum, elem, nnum, nxth		!Elementnummer bestimmen und in die Tabelle schreiben
	*vfill, rforce2(j,1,k), data, nnum
*enddo

*do, i, 2, anz+1 			
	*get, koord, elem, rforce2(i,1,k), cent, z	!z-Koordinate des Federelements bestimmen
	*vfill, rforce2(i,2,k), data, koord

	*get, frk, elem, rforce2(i,1,k), smisc, 1 	!Elementkraft bestimmen
	*if, frk, lt, 0, then				!negative Kräft positiv machen und positive zu Null
		*vfill, rforce2(i,3,k), data, -frk
	*else
		*vfill, rforce2(i,3,k), data, 0
	*endif

	*get, dehn2, elem, rforce2(i,1,k), nmisc, 1 
	*vfill, rforce2(i,4,k), data, dehn2

*enddo

rforce2(1,1,k) = 0	!restliche Tabelle ausfüllen um grafische Darstellung zu verbessern
rforce2(1,2,k) = 0
rforce2(1,3,k) = 0
rforce2(1,4,k) = 0
rforce2(anz+2,1,k) = 0
rforce2(anz+2,2,k) = l
rforce2(anz+2,3,k) = 0
rforce2(anz+2,4,k) = 0

*enddo

!---radiale Auswertung nach Koordinate

*dim, trag2, array, an, 14	!array anlegen um die Werte zu speichern

*do, i, 1, an
	wink = (360/an)*(i-1)		!Winkelkoordinate des Wälzkörpers berechnen und eintragen
	*vfill, trag2(i,1), data, wink
*enddo

*do, h, 1, 13
	*do, k, 1, an
		*vfill, trag2(k,h+1), data, rforce2(h,3,k)	!Kraftwert je Wälzkörper auslesen und in Array speichern
	*enddo
*enddo
allsel

!---axiale Auswertung ja Nadel in alle Richtungen

anz = es+1
*dim, rforce12, array, anz+2, 6, an

*do, k, 1, an
	lsel, s, , , k+1000	!Nodes selektieren die zu dem jeweiligen Federelement auf der Linie gehören
	nsll, s, 1
	esln, s, 0
	esel, r, type, , 11
	nsle, s, 1
	nsel, r, loc, x, dz/2

nnum = 0
*do, j, 2, anz+1
	*get, nnum, node, nnum, nxth
	*vfill, rforce12(j,1,k), data, nnum
*enddo

*do, i, 2, anz+1 			!Reaktionskräfte des jeweiligen Nodes bestimmen und abspeichern
	nsel, s, node, , rforce12(i,1,k)
	fsum, , 
	*get, frx, fsum, 0, item, fx
	*get, fry, fsum, 0, item, fy
	*get, frz, fsum, 0, item, fz 
	*get, koord, node, rforce12(i,1,k), loc, z
	!dru = druck(1,1)+(druck(2,1)*koord)
	
	*vfill, rforce12(i,2,k), data, koord
	*vfill, rforce12(i,3,k), data, frx
	*vfill, rforce12(i,4,k), data, -fry
	*vfill, rforce12(i,5,k), data, frz
	*vfill, rforce12(i,6,k), data, dru
*enddo

rforce12(1,1,k) = 0
rforce12(1,2,k) = 0
rforce12(1,3,k) = 0
rforce12(1,4,k) = 0
rforce12(anz+2,1,k) = 0
rforce12(anz+2,2,k) = l
rforce12(anz+2,3,k) = 0
rforce12(anz+2,4,k) = 0

*enddo

*use, schreib22.mac, nr
*use, schreib32_neu.mac,	!Ausgabe in Datei
*use, schreib42.mac, nr