!aus dem Forum von pettaa
!Kreisbogen mit 300 Punkten
!Schätzing,  21.9.2014

finish  $/clear
*abbr,pettaa2,/input,pettaa2,txt


!-----liest Kräftefile(kraefte) in Array(loads) ein----- 
*VREAD, loads(1), kraefte,txt
*VREAD, xcoord(1), xcoord,txt
*VREAD, ycoord(1), ycoord,txt

!-----Parameter-----

anz=*get,length of loads           !Anzahl der Punkte

radius=300e-3        !Radius 

/prep7
et,1,link1   

!-----Punkte setzen-----
*do,i,1,anz
k,i,xcoord(i-1),ycoord(i-1)
*enddo

!-----Kreis malen, dient nur der Anschauung-----
k,1001,0,0
*do,q,1,anz
larc,q,q+1,1001,radius
*enddo


!-------Vernetzung--------
lesize,all,,,1
lmesh,all

!-----tangentiale Kräfte-----------
*do,v,1,anz
   px=nx(v)
   py=ny(v)

   f,v,fx, loads(v)/radius*py
   f,v,fy,-loads(v)/radius*px
*enddo

!----Summation----
fxs=0
fys=0
*do,u,1,anz
   *get,fxi,node,u,f,fx
   *get,fyi,node,u,f,fy
   fxs=fxs+fxi
   fys=fys+fyi
*enddo

fxs=fxs
fys=fys