feap * * Beispiel gevouteter Stab
  0 0 0 2 2 2
! 0 0 0 2D-3D DOFs/Knoten Knotenanzahl/Element 

parameter
el= 4		! Anzahl der Elemente
k = el+1	! Anzahl der Knoten
E = 3*10^4	! E-Modul
A1 = 10	        ! Querschnitt links
A2 = 1	        ! Querschnitt rechts
L = 100	        ! Länge
P2= 20	        ! Kraft
i = 0

coord			! Knotenkoordinaten
  1 1    0    0		! Knotennr, Inkrem, x-Koord, y-Koord 
  k 0    L    0

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! Element 1
para
 i=i+1			! Elementnummer
 x=0.5*L*(2*i-1)/el	! Position in der Elementmitte
 A=(L-x)*A1+x*A2	! Querschnittswert in Elementmitte
 A=A/L			! Querschnittswert in Elementmitte

elem                    ! Knotenzuweisung zu Elementen
i 0 i i i+1

mate i			! Materialkarte
truss
elastic isotropic  E
cross section A

boun			! Randbedingungen
i+1 0 0 1


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! Element 2


para
 i=i+1				! Elementnummer
 x=0.5*L*(2*i-1)/el	! Position in der Elementmitte
 A=(L-x)*A1+x*A2		! Querschnittswert in Elementmitte
 A=A/L			! Querschnittswert in Elementmitte

elem
i 0 i i i+1


mate i			! Materialkarte
  truss
    elastic isotropic  E
    cross section A

boun			! Randbedingungen
i+1 0 0 1

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! Element 3
para
 i=i+1			! Elementnummer
 x=0.5*L*(2*i-1)/el	! Position in der Elementmitte
 A=(L-x)*A1+x*A2	! Querschnittswert in Elementmitte
 A=A/L			! Querschnittswert in Elementmitte

elem                    ! Knotenzuweisung zu Elementen
i 0 i i i+1

mate i			! Materialkarte
truss
elastic isotropic  E
cross section A

boun			! Randbedingungen
i+1 0 0 1


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! Element 4


para
 i=i+1				! Elementnummer
 x=0.5*L*(2*i-1)/el	! Position in der Elementmitte
 A=(L-x)*A1+x*A2		! Querschnittswert in Elementmitte
 A=A/L			! Querschnittswert in Elementmitte

elem
i 0 i i i+1


mate i			! Materialkarte
  truss
    elastic isotropic  E
    cross section A

boun			! Randbedingungen
i+1 0 0 1



cboun			! Randbedingungen
node 0 0 1 1

cforce			! Einzellast	
node L 0 P2 0		! Richtung (x,y-Koord), Kraft-x, Kraft-y


end				! Ende der Netzgenerierung

batch				! Solver
  tangent
  form
  solve
end

batch				! Post
  disp,all
  stre,all
  reac,all
!  plot cont 1
!  plot load
!  plot boun
!  plot mesh
!  plot elem
!  plot nodes
end

inter				! Ausgabe auf den Bildschirm
stop