feap * * Scheibe * *
  0 0 0 2 2 4
! 0 0 0 2D-3D DOFs/Knoten Knotenanzahl/Element 

global		! Nur bei 2D
 plane strain	! stress

parameter
n1 = 3		! Anzahl der Elemente
n2 = 2		! Anzahl der Elemente
a  = 40		! Abmessung
b  = 50		! Abmessung
d  = 1		! Dicke
f  = 1		! Kraft F
E  = 1		! E-Modul
nu = 1/3.0		! Querdehnzahl

block			! Erzeugt einen Block mit kart-Koord
 cart n1 n2  
  1    0   0	! Pkt. x-Koord. y-Koord.
  2    a   0	! Pkt. x-Koord. y-Koord.
  3    a   b	! Pkt. x-Koord. y-Koord.
  4    0   b	! Pkt. x-Koord. y-Koord.

eboun			! Randbedingungen
  1 0 1 1		! Richtung (x,y-Koord), Länge, DOF-x, DOF-y

cforce		! Knotenlast
  NODE a b/2.0 0 f	! NODE x-Koord. y-Koord. x-Kraft y-Kraft

mate 1		! Materialkarte
  solid
    elastic isotropic  E nu
    thick section d

end			! Ende der Netzgenerierung

batch			! Solver
  tangent
  form
  solve
  disp,all		! Daten ins O-File
  stre all		! Daten ins O-File
  reac all		! Daten ins O-File
end

batch			! Post
 plot postscript
!  plot cont 2
  plot stress
  plot load
  plot boun
  plot mesh
  plot node
  plot elem
  plot axis
 plot postscript
end

batch               ! Kontur-Plot der y-Verschiebung in eps-Datei
 plot postscript
  plot cont 2       
 plot postscript
end

batch               ! Verschiebungen bei x=a
   disp,coor,1,a 
end

batch               ! Spannung am Knoten 9
   stress,node,9
end

batch               ! Spannung entlang der y-Koord.
   stress,coor,1,0
end

inter			! Ausgabe auf den Bildschirm
stop