*HEADING
 EARTHQUAKE ANALYSIS  S8R
*NODE
*NSET,NSET=ENDS
*NGEN,NSET=NALL
*ELSET,ELSET=PR
*ELEMENT,TYPE=S8R
*ELGEN,ELSET=EALL
*MATERIAL,NAME=A1
*ELASTIC
30.E6,
*DENSITY
 .000728,
*SHELL SECTION,MATERIAL=A1,ELSET=EALL
2.,
*BOUNDARY
21, 1, 6
121, 1, 6
221, 1, 6
NALL,3,5
** AMPLITUDE CURVE DATA FROM CANTILEVER_QUAKEDATA
** THIS INPUT ASSUMES THAT CANTILEVER_QUAKEDATA.INP
** HAS BEEN COPIED TO QUAKE.AMP
*AMPLITUDE,VALUE=RELATIVE,INPUT=QUAKE.AMP,NAME=EQ
*RESTART,WRITE,FREQUENCY=25
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS
10,
*EL PRINT,FREQUENCY=0
*NODE PRINT
U,
*MODAL FILE
*output,field
*node output
u,
*END STEP
*STEP
 EARTHQUAKE INPUT
*MODAL DYNAMIC
.01,2.0
*SELECT EIGENMODES,GENERATE
*MODAL DAMPING
*BASE MOTION,DOF=1,AMPLITUDE=EQ,SCALE=386.09
*PRINT,FREQUENCY=50
*EL PRINT,ELSET=PR,FREQUENCY=50
S,E
*NODE PRINT,FREQUENCY=100
U,
TU,
V,
A,
*output,field,frequency=100
*node output
u,tu,v,a
*output,history
*node output,nset=ends
U,V,A,TU
*modal output
GU,GV,GA
*NODE FILE,NSET=ENDS
U,V,A,TU
*MODAL FILE
GU,GV,GA
*END STEP