finish
/clear,start

/prep7

et,1,solid185

mp,ex,1,210000      ! E-Modul in N/mm^2
mp,prxy,1,0.3        ! Querkontraktion
mp,dens,1,7.8 ! Dichte in g/mm^3

block,0,1,0,1,0,1,  
block,0,1,0,1,1,2,  
block,0,1,0,1,2,3,

mshkey,1
vmesh,all

nummrg,node
nummrg,kp

!-----Randbedingung-----
nsel,s,loc,y,0
d,all,all
!------Lasten------
nsel,s,loc,y,1
d,all,uy,0.01


!--------Lösung-------
allsel
/solu
solve

!------Auswertung-------
/post1

allsel
!vsel,s,,,1
!vsel,s,,,2
vsel,s,,,3
nslv,r,1

*get,KNOTEN,node,,count ! Knotenanzahl wird ermittelt
*dim,STAT_KNOTEN,array,KNOTEN		! Erstellen der Maske
*dim,STRESS_KNOTEN,,KNOTEN,3		! Array zum Festhalten der Spannungen
*dim,RESULTS,,KNOTEN,3			! Array in dem Ergebnisse eingefügt werden

*get,n_min,node,,num,min
n_arb=n_min

*do,i,1,KNOTEN

STRESS_KNOTEN(i,1)=n_arb
*get,spannung_x,node,n_arb,s,x
STRESS_KNOTEN(i,2)=spannung_x
*get,spannung_y,node,n_arb,s,y
STRESS_KNOTEN(i,3)=spannung_y
n_arb=ndnext(n_arb)

*enddo


*do,I,1,KNOTEN
SIGMA_1 = STRESS_KNOTEN(I,2)
SIGMA_2 = STRESS_KNOTEN(I,3)

I_1 = SIGMA_1-SIGMA_2
I_2 = SIGMA_2/2

*vfill,RESULTS(I,1),data,STRESS_KNOTEN(I,1)
*vfill,RESULTS(I,2),data,I_1
*vfill,RESULTS(I,3),data,I_2

*enddo


!*stat,stress_knoten