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

fini  $/cle
*abbr,feldlinie,/inp,feldlinie.txt

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

!-----------------------------------
/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,sw/10
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
!plve,ef

!/eof

!---Vereinbarungen-----
*dim,ort,array,100,3
e_mittel=spannung/sw

!---------------------
!Startpunkt festlegen
csys,1
n_start=node(rad,-45,0)
csys,0
ort(1,1)=1
ort(1,2)=nx(n_start)
ort(1,3)=ny(n_start)

n_i=n_start

*do,i,2,40

*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/10
ort(i,3)=ort(i-1,3)+ey/e_mittel*sw/10

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

*enddo



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

lplot
*do,i,1,39
l,1000+i,1001+i
*enddo