feap *Schallwelle Dehnstab*
  0 0 0 2 2 2

parameter
 n  = 101		! Anzahl der Knoten
 m  = n-1		! Anzahl der Elemente
 l  = 1000		! Laenge in mm
 b  = 100		! Breite in mm
 h  = 100		! Höhe in mm
 E  = 210000.0		! E-Modul in N/mm^2 bzw. MPa
 A  = b*h		! Querschnittsflaeche in mm^2
 nu = 0.3		! Querdehnzahl
 ro = 7.85e-9		! Dichte t/mm^3
 cs = sqrt(E/ro)	! Schallgeschwindigkeit Stahl
 t1 = l/cs		! Zeit fuer Durchlauf Schallwelle
 t2 = 4.0*t1		! Zeit fuer vier Stabdurchlaeufe = Berechnungszeit
 tc = l/m/cs            ! kritischer Zeitschritt
 pr = 0.1		! Faktor fuer Zeitschrittweite pr % des kritischen Zeitschritts
 dt = pr*tc		! Zeitschrittweite
 nt = t2/dt		! Anzahl Zeitschritte
 np = 1			! Anzahl eps Plots
 nl = nt/np		! Anzahl Gleichungsloesungen pro Plott
 v0 = 1.0e+4		! Geschwindigkeit v0=10 m/s=1e4 mm/s
 zs = 400.0		! Anzahl zu speichernder Zeitschritte
 z  = nt/zs 		! TPLOt: zu jedem z-en Zeitschritt Daten speichern
 z  = 10.0 !nt/zs 	! TPLOt: zu jedem z-en Zeitschritt Daten speichern
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Knoten und Elemente
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 

coord			! Knotenkoordinaten festlegen
  1  1  0  0		! Knotennr.          , Inkr., x-Koord, y-Koord
  n  0  l  0		! Kn.-nr. letzter Kn., Inkr., x-Koord, y-Koord

elem
  1  1  1  1  2		! Elem.-nr., Inkr.,Mat.,Kn.-nr.1, Kn.-nr.2 
  m  0  1  m  n

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Randbedingungen
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~   
  
eboun			! Einspannung: linker Rand fest
  1 0 1 1

eboun			! Fuehrung: alle dof in y fest
  2 0 0 1
  
eboun                   ! Für Verschiebungsrand: 
  1 1000 1 1            ! x-Richtung Verschiebungsrand definieren=1 und y-Richtung fest=1

edisp
1  1000  -v0  0         ! Definition des Verschiebungsrands | Funktion u(t) s. u. im Befehl PROP

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

mate 1
  truss
    elastic isotropic  E
    cross section      A
    density mass       ro

end
!*************************************************************
!**                     ENDE NETZGENERIERUNG                **
!*************************************************************

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

 batch				! Bildung der Belastungsfunktion
   prop,,1
 end
 2,1				! 2 gibt den Typ an (Tabelle); 1 Anzahl der Wertepaare pro Zeile
 0.0 ,  0.0		! Wertepaar: Zeitpunkt Temperatur
 t2  , 1.0*t2
 

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Abspeichern der Spannungen pro Element über der Zeit
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~   
 
 batch				! Herausschreiben der Daten zum z-en Zeitschritt
  tplot,,z
 end
   DISP,  coor, 1, 1000.0	! Verschiebung am freien Stabende
   STRESS,1,1			! stress, elem, component (1=sig_xx)
   STRESS,2,1
   STRESS,3,1
   STRESS,4,1
   STRESS,5,1
   STRESS,6,1
   STRESS,7,1
   STRESS,8,1
   STRESS,9,1
   STRESS,10,1
   STRESS,11,1
   STRESS,12,1
   STRESS,13,1
   STRESS,14,1
   STRESS,15,1
   STRESS,16,1
   STRESS,17,1
   STRESS,18,1
   STRESS,19,1
   STRESS,20,1
   STRESS,21,1
   STRESS,22,1
   STRESS,23,1
   STRESS,24,1
   STRESS,25,1
   STRESS,26,1
   STRESS,27,1
   STRESS,28,1
   STRESS,29,1
   STRESS,30,1
   STRESS,31,1
   STRESS,32,1
   STRESS,33,1
   STRESS,34,1
   STRESS,35,1
   STRESS,36,1
   STRESS,37,1
   STRESS,38,1
   STRESS,39,1
   STRESS,40,1
   STRESS,41,1
   STRESS,42,1
   STRESS,43,1
   STRESS,44,1
   STRESS,45,1
   STRESS,46,1
   STRESS,47,1
   STRESS,48,1
   STRESS,49,1
   STRESS,50,1
   STRESS,51,1
   STRESS,52,1
   STRESS,53,1
   STRESS,54,1
   STRESS,55,1
   STRESS,56,1
   STRESS,57,1
   STRESS,58,1
   STRESS,59,1
   STRESS,60,1
   STRESS,61,1
   STRESS,62,1
   STRESS,63,1
   STRESS,64,1
   STRESS,65,1
   STRESS,66,1
   STRESS,67,1
   STRESS,68,1
   STRESS,69,1
   STRESS,70,1
   STRESS,71,1
   STRESS,72,1
   STRESS,73,1
   STRESS,74,1
   STRESS,75,1
   STRESS,76,1
   STRESS,77,1
   STRESS,78,1
   STRESS,79,1
   STRESS,80,1
   STRESS,81,1
   STRESS,82,1
   STRESS,83,1
   STRESS,84,1
   STRESS,85,1
   STRESS,86,1
   STRESS,87,1
   STRESS,88,1
   STRESS,89,1
   STRESS,90,1
   STRESS,91,1
   STRESS,92,1
   STRESS,93,1
   STRESS,94,1
   STRESS,95,1
   STRESS,96,1
   STRESS,97,1
   STRESS,98,1
   STRESS,99,1
   STRESS,100,1
 show

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Implizite Zeitdiskretisierung zur Gleichungsloesung
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~    
 
BATCH				! starten der Berechnung
check				! Ueberpruefung der Knotennum. und Elemente
   DT,,dt			! Zeitschrittweite
   INITIAL, DISPLACEMENTS
   INITIAL, Rates
   MASS
   TRANsient,newmark		! implizite Zeitintegration mit NEWMARK-Verfahren
   LOOP,print,np                ! Schleife ueber Plott-Zustaende
      LOOP,time,nl		! Schleife ueber die Zeit
	TIME
	  TANG,,1               ! Gleichungsloeser
      NEXT time
!    DISP,all		  	! Output
!    STRE,all		  	! Output
!    REAC,all		  	! Output
!   plot post
!   plot default off
!   plot rang,off
!   plot rang,-1000,0
!   plot cont,-1400,,1400
!   plot,defa,off		! Turn off defaults for ranges
!   plot,cont,1,,1		! contour axial displacment and set ranges
!   plot defo 100 1
!   plot mesh
!   plot nodes
!   plot stre,1
!   plot wipe
!   plot post
!   plot wipe
!   plot clear
   NEXT print
END

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Postprozess
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~   

! BATCH				! Postprozess: PLOT Einstellungen
!PLOT, post
!   PLOT, pers,1
! !   PLOT, PERSpectiv
!   PLOT, axis
! !   PLOT, hide
!   PLOT, MESH
!   PLOT, NODE
! !   PLOT, ELEM
!   PLOT, BOUN
! !   PLOT, CONT
!   PLOT, LOAD
! !  PLOT, STRE 
!PLOT, post
!  END

! 1500, -1000, 500		! Gibt die Entfernung zum Objekt an für x,y,z
! 0, 0, 1			! Gibt die Achse, die nach oben zeigen soll an

inter
stop