FINISH
/CLEAR
/CONFIG,NRES,20000
/units,si

!----PREPROZESSOR STARTEN
/PREP7
!----BESCHLEUNIGUNGSWERTE EINLESEN
!----PARAMETER
  
anz=2048 !Anzahl der Werte , abhängig vom Akzelerogramm
*DIM,avlf,array,anz,4
*VREAD,avlf(1,1),Beschleunigungszeitverlauf2,txt,,jik,4,anz,,21
(1f12.6,1f15.8,1f15.8,1f15.8)

! ----MATERIALEIGENSCHAFTEN

! Material Stahl
MP,DENS,1,7850    ! Dichte [kg/m^3] 
MP,EX,1,2.1e11    ! E-Modul in [N/m^2]
MP,GXY,1,8.1e10   ! Schermoduli [N/m^2]
MP,PRXY,1,0.3

/TITLE BEAM3RAHMEN

ET,1,BEAM3
!AREA ,IZZ ,HEIGHT ,SHEARZ ,ISTRN, ADDMAS  
R,1,0.02265,0.00141200,0.59

!R,2,A2,FTMIZZ,Hoehe2

!----KNOTEN
real,1
type,1
mat,1

K,1,0,0,0
K,2,0,10.0
K,3,12.5,11.09,0
K,4,25,10,0
K,5,25,0,0

!----ELEMENTE

L,1,2
L,2,3
L,3,4
L,4,5
TYPE,1
REAL,1
MAT,1
LESIZE,all,,,20
LMESH,all


!----AUFLAGERBED. INITIAL
D,Node(0,0,0),ux
D,Node(0,0,0),uy
D,Node(25,0,0),ux
D,Node(25,0,0),uy

/SOLU
ANTYPE,TRANS 
TRNOPT,FULL   
!----INITIAL condition herstellen
NSUB,5
ALPHAD,0.1378      !Koeff. für Masssendämpfung,Rayleigh
BETAD,0.000348     !Koeff. für Steifigkeitsdämpfung,Rayleigh     
 
OUTRES,all,all      
AUTOTS,off  
pstres,on   

!----INITIAL condition herstellen
!TIME,0.1 
!TIMINT,OFF        ! Schaltet Trägheitseffekte aus
!SSTIF,ON
!KBC,0  
NSUB,5 
!ACEL,0,9.81,0
!SOLVE
!LSWRITE

!----EIGENTLICHE BELASTUNGSPHASE
TIMINT,ON ! Schaltet Trägheitseffekte ein
*DO,N,1,100,1 !anz
 TIME,avlf(n,1)     
 !Beschleunigung [m/sec²]
 ACEL,((avlf(n,2))/100),(1.0*(avlf(n,4))/100)0, ! Horozontal- und Vertikalbeschl.
 SOLVE
*ENDDO

FINISH

!----ERGEBNISS

!siehe Auswert.mac
!auslesen aller Rahmenmomente und Verschiebungen
!und ablegen in eine Ergebnisdatei

/EOF