FEAP *Stand. lin. Festkörper Relaxation*
0 0 0 2 2 2

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Parameter und Variablen
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 

PARAMETER
dt = 0.01       ! Zeitschritt in s
t2 = 20	        ! Berechnungszeit in s
nt = t2/dt      ! Anzahl der Berechnungsschritte
u  = 1          ! Vorgabe der Endverschiebung in mm
E0 = 7000.0     ! Elastizitaetsmodul N/mm^2
E1 = 5000.0     ! Elastizitaetsmodul N/mm^2
E  = E0+E1      ! Elastizitaetsmodul N/mm^2
v1 = 7000.0     ! Viskositaet in Ns/mm^2
b1 = E1/E       ! Beteiligungsfaktor
t1 = v1/E1      ! Relaxationszeit in s
L  = 100        ! Laenge in mm
A  = 2          ! Querschnittsflaeche in mm^2

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Element- und Materialdefinition
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~   

MATERIAL, 1
 TRUSS
  ELASTIC ISOTROPIC     E           ! Elastizitaetmodul
  CROSS SECTION         A           ! Querschnittsflaeche
  VISCOELASTIC term1    b1 t1       ! Beteiligungsfaktor / Relaxationszeit
 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Knoten und Elemente
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  

COORDINATES
 1  0  0  0  0      ! Knotennr., Inkr., x, y, z
 2  0  L  0  0      ! Knotennr., Inkr., x, y, z

ELEMENTS
 1  0  1  1  2      ! Elementnr., Inkr., Material, Knotennr., Knotennr.

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Randbedingungen
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

BOUNDARY
 1  0  1  1    	! Knotennr., Inkr., Randbed.-x, Randbed.-y. Randbed.-z
 2  0  1  0        ! Knotennr., Inkr., Randbed.-x, Randbed.-y. Randbed.-z

DISPLACEMENT
 2  0  u  0        ! Knotennr., Inkr., Verschiebung-x, Verschiebung-y. Verschiebung-z

END

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Belastungs-Zeitfunktion
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  

BATCH
 PROP,,1        ! Belastungs-Zeit-Funktion
END 
 2 1            ! 1: Funktion / 2: Tabelle, Anzahl Wertepaare
 0	1           ! Zeit, Skalierungsfaktor
 t2	1           ! Zeit, Skalierungsfaktor

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Abspeichern der Ergebnisse pro Zeitschritt
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 

BATCH
 TPLOT,,1
END
 DISP 2 1       ! Verschiebung, Knotennr., Achse
 STRESS 1 1 	! Spannung, Elementnr., Achse
 STRESS 1 2 	! Spannung, Elementnr., Achse
 STRESS 1 3     ! Spannung, Elementnr., Achse
SHOW

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Explizite Zeitdiskretisierung zur Gleichungsloesung
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 

BATCH
 DT,,dt             ! Zeitinkrement
 LOOP TIME nt		! Schleifendurchlaeufe
  TIME              ! Aufaddieren der Zeit
   LOOP ITER 10		! Anzahl der Newtoniterationen
     TANG,,1		! Steifigkeitsmatrix bilden, residual und loesen
   NEXT ITER		! Naechster Iterationsschritt
  DISP,all          ! Output
  STRE,all          ! Output
  REAC,all          ! Output
 NEXT TIME          ! Naechster Zeitschritt
END
 
INTERactive
STOP