feap * * Übung ebene 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 
  disp,coor,1,0   ! Ausgabe der Verschiebung bei x=0
  stress,coor,1,0 ! Ausgabe der Hauptspannung bei x=0
  disp,node,1     ! Ausgabe der Spannungen für Knoten 1
  stress,node,1   ! Ausgabe der Spannungen für Knoten 1
end

batch
 plot postscript
 plot mesh
 plot axis
 plot load
 plot boun
 plot postscript
end

batch
 plot postscript
 plot cont,3     ! Kontourplot Verschiebung u_z
 plot postscript
end

inter			! Ausgabe auf den Bildschirm
stop