!Ermittlung einer Feldlinie
!ein erster Versuch
!
!Forum blid
!Schätzing, 27.11.2009

finish
/clear,start

!fini  $cle
!abbr,feldlinie,inp,feldlinie.txt

rad=20e-3
sw=12e-3
spannung=1

winkel=30*3.14159/180/2			! Winkel für Partikelstartpunkt
x_start=rad*sin(winkel)
y_start=2*rad-rad*cos(winkel)

!-----------------------------------
/prep7
et,1,plane67,,,1
mp,rsvx,1,1

!---Geometrie-----
pcirc,rad,,-90,0
rectng,0,rad*2,-rad-sw,0
asba,2,1

!----vernetzung-----
!esize,sw10
esize,sw/20
amesh,all

!-----Randwertvorgaben-----
nsel,s,loc,y,-sw-rad
d,all,volt,0

csys,1
nsel,s,loc,x,rad
csys,0
d,all,volt,Spannung

allsel
!----------------------
/solu
solve          !Lösung des elektrostat. Feldes

!----Auswertung---------

/post1
!plns,volt,1
!plns,ef,sum
PLVECT,ef

!eof

!---Vereinbarungen-----

*dim,ort,array,5000,3
e_mittel=spannung*sw


!---------------------
!Startpunkt festlegen

!n_start=192

csys,1
n_start=node(rad,-80,0)
csys,0


ort(1,1)=1
ort(1,2)=nx(n_start)
ort(1,3)=ny(n_start)

n_i=n_start
n_end=0

*do,i,2,500

!*dowhile ort(i-1,3) > -rad-sw

*get,ex,node,n_i,ef,x  !Feldstärke X-Komponente
*get,ey,node,n_i,ef,y  !Feldstärke X-Komponete

!!!ort(i,1)=i
!!!ort(i,2)=ort(i-1,2)+ex/e_mittel/sw/10000
!!!ort(i,3)=ort(i-1,3)+ey/e_mittel/sw/10000

ort(i,1)=i
ort(i,2)=ort(i-1,2)+(0.5e-4)*ex/SQRT(ex*ex+ey*ey)*10
ort(i,3)=ort(i-1,3)+(0.5e-4)*ey/SQRT(ex*ex+ey*ey)*10

n_i=node(ort(i,2),ort(i,3),0)

*IF,ort(i,3),LT,(-rad-sw),EXIT
!*IF,n_end,EQ,0
!letzerort=i
!*ENDIF

*ENDDO

ort(i,3)=-sw-rad


/prep7
*do,ii,1,i
k,1000+ii,ort(ii,2),ort(ii,3)
*enddo

lplot
*do,ii,1,i-1
l,1000+ii,1001+ii
*enddo